一種地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法
【專利摘要】本發(fā)明公開了一種地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法，利用非數(shù)值計算方法，對虛擬地質環(huán)境中的對象及其之間的相互關系加以描述并構建3D模型，抽象的地質對象主要包括點、邊、面片、表面和體及其派生類的對象，具體包括地層、斷層、侵入巖、透鏡體、礦體，在構建3D模型時，需要在各個構建過程中進行數(shù)據(jù)檢測，并指導數(shù)據(jù)及模型的修正，主要包括多源數(shù)據(jù)一體化、布點檢測、表面模型檢測、體模型檢測。通過對地質體網(wǎng)格化過程中各個階段的分析，提出相應的數(shù)據(jù)檢測及修正方法，以確保網(wǎng)格化能夠正確剖分，并使得網(wǎng)格單元的精度能夠滿足地質工作者的需求。
【專利說明】
一種地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法
技術領域
[0001] 本發(fā)明涉及一種網(wǎng)格化技術，尤其涉及一種地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方 法。
【背景技術】
[0002] 網(wǎng)格化技術是一種重要的信息處理方法，可以用來表示復雜物體的空間數(shù)據(jù)。地 質體的網(wǎng)格化是將地質體分解為較小的地質單元，單元之間存在關聯(lián)關系，以離散的形式 來描述和分析巖性、斷層之類的地質特征，以便進行地質現(xiàn)象的數(shù)值計算、可視化渲染及其 它工程應用，為開展三維地質模擬奠定基礎，使地質工作者可以獲得更多的信息。
[0003] 網(wǎng)格化作為地質工程實踐應用的前處理階段，網(wǎng)格單元的質量直接影響著工程實 際問題的解決精度。在網(wǎng)格化技術發(fā)展的早期階段，網(wǎng)格模型通常是由熟練的工程師手工 劃分，之后交由算法進行自動處理。由于地質現(xiàn)象的復雜性和多樣性，導致網(wǎng)格化問題的規(guī) 模與復雜性也在日益增加。在網(wǎng)格化過程中，按照網(wǎng)格之間的連接關系可分為結構化網(wǎng)格 與非結構化網(wǎng)格。結構化網(wǎng)格本身隱含著固定的鄰接關系，以四邊形、立方體網(wǎng)格單元為 主，結構簡單，網(wǎng)格化方法包括代數(shù)法、偏微分方程法以及超單元映射法，計算結果易于收 斂，但是網(wǎng)格單元的精度難以控制與檢測。非結構化網(wǎng)格的鄰接關系并不固定，需要記錄其 鄰接關系，以四面體、六面體、棱錐體為主要網(wǎng)格單元，能較好表達地質體的邊界特征，主要 方法有Delaunay法和AFT法，具有邊界適應能力，是目前網(wǎng)格化技術的主要研究方向，其中， 如何檢測復雜地質體剖分時的邊界一致性和生成高質量網(wǎng)格單元是難點問題之一。目前， 在網(wǎng)格化過程中，無論是樣本數(shù)據(jù)的收集、解釋及其處理，還是網(wǎng)格模型的構建，仍然缺乏 有效而系統(tǒng)性的檢測手段和方法。

【發(fā)明內容】

[0004] 本發(fā)明的目的是提供一種網(wǎng)格化能夠正確剖分、網(wǎng)格單元的精度高的地質體網(wǎng)格 化過程中的數(shù)據(jù)檢測方法。
[0005] 本發(fā)明的目的是通過以下技術方案實現(xiàn)的：
[0006] 本發(fā)明的地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法，利用非數(shù)值計算方法，對虛擬地 質環(huán)境中的對象及其之間的相互關系加以描述并構建3D模型，抽象的地質對象主要包括 點、邊、面片、表面和體及其派生類的對象，具體包括地層、斷層、侵入巖、透鏡體、礦體，在構 建3D模型時，需要在各個構建過程中進行數(shù)據(jù)檢測，并指導數(shù)據(jù)及模型的修正，主要包括步 驟：
[0007] A、多源數(shù)據(jù)一體化：
[0008] 3D地質建模主要基于區(qū)域內的多源數(shù)據(jù)，該區(qū)域記為Ω，將Ω內的邊界數(shù)據(jù)、鉆孔 數(shù)據(jù)、剖面數(shù)據(jù)、3D掃描數(shù)據(jù)、以及各種解釋數(shù)據(jù)進行一體化處理，生成樣本數(shù)據(jù)集合S，作 為地質體網(wǎng)格化的基礎數(shù)據(jù)；
[0009] B、布點檢測：
[0010] 由于Ω內的樣本數(shù)據(jù)比較缺乏，需要在其中布置數(shù)據(jù)點，采用2種方法：
[0011] 方法一、特征點插入:識別Ω內的地質特征，如斷層、尖滅，并根據(jù)給定的密度閾值 ε，插入相應點，構建點集合F;
[0012] 方法二、規(guī)格化點插入：以Ω最小/大值為界，構建包圍盒，對包圍盒內進行規(guī)格化 插入點，生成點集合G;
[0013] C、表面模型檢測：
[0014] 以S為網(wǎng)格化的點數(shù)據(jù)，采用三角形或四邊形作為面片，進行網(wǎng)格化剖分，構建地 質體的表面模型，并進行下列檢測：
[0015] 約束檢測、局部單元優(yōu)化檢測、閉合性檢測、自交性檢測、流型檢測、面片方向性檢 測；
[0016] D、體模型檢測：
[0017]對已構建地質體的表面模型進行上述檢測并修正完成之后，可以進行體剖分，生 成相應的四面體或六面體網(wǎng)格模型，并對其進行下列測試：
[0018] 網(wǎng)格優(yōu)化檢測、吻合度檢測、截面檢測。
[0019] 由上述本發(fā)明提供的技術方案可以看出，本發(fā)明實施例提供的地質體網(wǎng)格化過程 中的數(shù)據(jù)檢測方法，利用非數(shù)值計算方法，對虛擬地質環(huán)境中的對象及其之間的相互關系 加以描述并構建3D模型，抽象的地質對象主要包括點、邊、面片、表面和體及其派生類的對 象，具體包括地層、斷層、侵入巖、透鏡體、礦體，在構建3D模型時，需要在各個構建過程中進 行數(shù)據(jù)檢測，并指導數(shù)據(jù)及模型的修正。通過對地質體網(wǎng)格化過程中各個階段的分析，提出 相應的數(shù)據(jù)檢測及修正方法，以確保網(wǎng)格化能夠正確剖分，并使得網(wǎng)格單元的精度能夠滿 足地質工作者的需求。
【附圖說明】
[0020] 圖1為本發(fā)明實施例提供的地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法的總體流程示意 圖。
[0021] 圖2為本發(fā)明實施例中表面模型的拓撲結構示意圖。
[0022]圖3a、圖3b、圖3c分別為本發(fā)明實施例中懸點、懸邊、懸面片的約束檢測示例。
[0023 ]圖4a、圖4b、圖4c、圖4d分別為本發(fā)明實施例中穿透性自交、相鄰性自交、穿透性自 交重構、相鄰性自交重構的自交性檢測示例。
[0024]圖5為本發(fā)明實施例中流型檢測示例。
[0025]圖6a、圖6b分別為本發(fā)明實施例共面檢測中初始前沿、前沿推進法過程示意圖。 [0026]圖7a、圖7b分別為本發(fā)明實施例的多源數(shù)據(jù)一體化處理中的一體化處理示意圖和 地層樣本數(shù)據(jù)示例。
[0027]圖8a、圖8b分別為本發(fā)明實施例的布置數(shù)據(jù)處理及檢測中的特征點插入和規(guī)格化 點插入示意圖。
[0028]圖9為本發(fā)明實施例中表面模型拓撲結構信息示例。
[0029]圖10a、圖IOb、圖IOc、圖IOd分別為本發(fā)明實施例對模型S4進行約束檢測示例中S4 中的冗余面片示例、S沖的懸點示例、S4中的懸邊示例、約束檢測并修正后的S4網(wǎng)格示例。 [0030]圖11a、圖Ilb分別為本發(fā)明實施例對S 4進行閉合性檢測中閉合性檢測失敗和閉合 性檢測成功示意圖。
[0031] 圖12為本發(fā)明實施例中多個地層模型的自交性檢測示意圖。
[0032] 圖13為本發(fā)明實施例中對S4地層模型實行流型檢測示意圖。
[0033]圖14為本發(fā)明實施例中S4的體模型。
[0034]圖15a、圖15b、圖15c、圖15d分別為本發(fā)明實施例截面檢測與一體化顯示中S4截面 示例、單截面一體化示例、多個截面及其一體化示例、多個截面及其一體化示例的局部放大 示意圖。
【具體實施方式】
[0035] 下面將對本發(fā)明實施例作進一步地詳細描述。
[0036] 本發(fā)明的地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法，其較佳的【具體實施方式】是：
[0037] 利用非數(shù)值計算方法，對虛擬地質環(huán)境中的對象及其之間的相互關系加以描述并 構建3D模型，抽象的地質對象主要包括點、邊、面片、表面和體及其派生類的對象，具體包括 地層、斷層、侵入巖、透鏡體、礦體，在構建3D模型時，需要在各個構建過程中進行數(shù)據(jù)檢測， 并指導數(shù)據(jù)及模型的修正，主要包括步驟：
[0038] A、多源數(shù)據(jù)一體化：
[0039] 3D地質建模主要基于區(qū)域內的多源數(shù)據(jù)，該區(qū)域記為Ω，將Ω內的邊界數(shù)據(jù)、鉆孔 數(shù)據(jù)、剖面數(shù)據(jù)、3D掃描數(shù)據(jù)、以及各種解釋數(shù)據(jù)進行一體化處理，生成樣本數(shù)據(jù)集合S，作 為地質體網(wǎng)格化的基礎數(shù)據(jù)；
[0040] B、布點檢測：
[0041] 由于Ω內的樣本數(shù)據(jù)比較缺乏，需要在其中布置數(shù)據(jù)點，采用2種方法：
[0042] 方法一、特征點插入:識別Ω內的地質特征，如斷層、尖滅，并根據(jù)給定的密度閾值 ε，插入相應點，構建點集合F;
[0043] 方法二、規(guī)格化點插入：以Ω最小/大值為界，構建包圍盒，對包圍盒內進行規(guī)格化 插入點，生成點集合G;
[0044] C、表面模型檢測：
[0045] 以S為網(wǎng)格化的點數(shù)據(jù)，采用三角形或四邊形作為面片，進行網(wǎng)格化剖分，構建地 質體的表面模型，并進行下列檢測：
[0046] 約束檢測、局部單元優(yōu)化檢測、閉合性檢測、自交性檢測、流型檢測、面片方向性檢 測；
[0047] D、體模型檢測：
[0048]對已構建地質體的表面模型進行上述檢測并修正完成之后，可以進行體剖分，生 成相應的四面體或六面體網(wǎng)格模型，并對其進行下列測試：
[0049] 網(wǎng)格優(yōu)化檢測、吻合度檢測、截面檢測。
[0050] 所述步驟B中：
[0051] 對于 VfieF & V_Sj.eS,其中，i = l，2, · · .m, j = l，2, · · .n〇，m為F集合的元素個數(shù)，η〇 為S的樣本初始個數(shù)，檢測F中的每個點6與Vsj的距離，如果該距離大于 ε，則將fl合并入s 中，即S = S U {fi}，最終獲得集合S，且其元素個數(shù)為m，m < m+ηο;
[0052] 對于 VggG & VsjeS.，:其中，k=l，2, · · .r, j = l，2, · · .m，r為G集合的元素個數(shù)，檢 測G中的每個點gk與Vs3的距離，如果該距離大于ε，則進行點 gk在地質體邊界區(qū)域內的判斷， 對于在Ω內的點gk，則將gk也合并入S中，即S = SU {gk}。
[0053] 所述步驟C中：
[0054] 所述約束檢測包括:對構成表面模型中的對象點、邊、及面片，分別進行冗余性和 多余性檢測，首先構建表面模型的拓撲結構，即表面模型可以由多個面片構成;面片由邊構 成，或直接由點構成;邊由點構成，通過表面模型的拓撲結構，能準確判斷它們之間的關聯(lián) 關系，在此基礎上，進行面片、邊、點的冗/多余性檢測；
[0055] 所述局部單元優(yōu)化檢測包括:在表面模型中遍歷搜索共享邊的面片單元，對兩個 共享邊的面片作最大空圓準則檢測，如果某個頂點在面片的外接圓之內，則修正對角線，完 成單元優(yōu)化過程的檢測，對于無法滿足最大空圓準則的，進行局部布點修正；
[0056] 所述閉合性檢測包括:從拓撲結構中任意一條邊ei出發(fā)，向上搜索與邊ei共享的面 片及其總數(shù)N，如果N! =2,則閉合性檢測失敗，輸出不閉合的邊或面片，否則繼續(xù)下條邊的 搜索，直至所有邊均搜索完畢，且每條邊僅被兩個面片所共享，說明閉合性檢測成功，即表 面模型具有閉合性，對于導致非閉合的邊或面片，需要修正表面模型中的數(shù)據(jù)，并重新進行 閉合性檢測；
[0057] 所述自交性檢測主要針對網(wǎng)格面片的自交性進行檢測，包括:對于Vci面片，i e U，2, ...，M}，M為面片集合的個數(shù)，在所構建表面模型的拓撲結構中，通過包圍盒檢測方 法，搜索W，j e {1，2，. . .，M}&&j矣i ;并分別計算cdPw的平面方程;再將2個平面方程求交， 并進行交點有效性判斷，獲得交點系列*P，如果P = NULU則cdPw不自交，否則返回自交信 息，對于自交的網(wǎng)格，按以下兩種方法進行網(wǎng)格重構，方法一是返回布點階段，重新插值及 修正數(shù)據(jù)，方法二是局部添加網(wǎng)格；
[0058]所述流型檢測包括:對于下列情況，即點、邊、面片都為模型的正確網(wǎng)格剖分，但是 由于多個網(wǎng)格模型的相鄰關系，使得存在某個邊被多于2個面片所共享，此時需要進行如下 檢測:設邊ek，ke{l，2, . . .，N}，N為邊集合的個數(shù)，在所構建表面模型的拓撲結構中，向上 搜索由其構成的面片，記為Ci面片，ie {1，2,...，M}，M為面片集合的個數(shù)，L為與邊％共享的 面片個數(shù)，如果L是奇數(shù)，則需要分解表面模型，直至L為偶數(shù)；
[0059] 所述面片方向性檢測包括:對于凸面體和凹面體的表面模型，首先過面片內一點 做垂直于面片的垂線，最少與模型的其它表面面片有1個交點，當交點個數(shù)大于1時，從面片 內一點開始的垂線，一端向著模型內部方向延伸，與模型有奇數(shù)個交點；另一端向著模型外 部方向延伸，與模型有偶數(shù)個交點，只要選取向著模型內部方向延伸的這端距離面片最近 的交點，就能得到正確的面片特殊內部點，之后以面片內一點為起點，以垂線與其它面片的 交點為終點做一條向量，把與面片法向量方向相同的向量對應的交點歸為一類，方向相反 的歸為一類，判斷哪一類中具有奇數(shù)個交點，從中選擇距離面片最近的交點來求特殊內部 點，再從這個特殊內部點向面片做垂線，根據(jù)垂線與面片法向量的方向關系做出調整，如果 垂線的方向與面片法向量的方向相同，則此面片法向量指向模型外部，否則交換任意兩個 頂點的位置。
[0060] 所述步驟D中：
[0061] 所述網(wǎng)格優(yōu)化檢測包括:計算形成網(wǎng)格單元兩個面片的夾角，以檢測該網(wǎng)格單元 的扭曲度;計算網(wǎng)格單元中最長邊與最短邊之比，檢測其比值是否為1-3;計算單位雅可比 值，檢測網(wǎng)格單元雅可比值是否在[0.4，1]范圍內，通過設置參數(shù)來控制生成網(wǎng)格單元的幾 何屬性，從而提高網(wǎng)格的質量，這些參數(shù)主要包括:局部網(wǎng)格尺寸、可接受的全局最大和最 小網(wǎng)格尺寸、網(wǎng)格密度、網(wǎng)格分級、元素曲率半徑、在封閉邊界上的網(wǎng)格精細化和三維網(wǎng)格 優(yōu)化的步數(shù)；
[0062] 所述吻合度檢測包括:面積計算，即分別計算體模型網(wǎng)格外邊界單元的面積、以及 表面模型單元的面積，檢測其差值是否趨近〇.〇;體積計算，即分別計算體模型和表面模型 的體積，并檢測其差值;共面檢測，這種檢測一般適用于體模型和表面模型采用相同網(wǎng)格單 元的情況，從TIN的網(wǎng)格初始前沿開始，逐層向區(qū)域內部推進，直至填滿整個區(qū)域，對于任意 一個面片V eieTIN，ie {1，2, . . .，M}，M為TIN的面片集合的個數(shù)，通過遍歷相應的TEN模型外 邊界面片，則一定存在且僅存在一個面片aupTEM，je {1，2, . . .，N}，N為TEN的面片集合的 個數(shù)，使得CdPuj-致;當檢測所有TIN中的面片與TEN模型外邊界面片中的面片存在上述一 致性關系時，返回完全吻合檢測信息；
[0063] 所述截面檢測包括:分別對體模型和表面模型進行截面計算與重構，并將它們進 行一體化顯示，通過檢測網(wǎng)格單元的面片、邊、點，進行可視化的對比分析。
[0064] 本發(fā)明通過對地質體網(wǎng)格化過程中各個階段的分析，提出相應的數(shù)據(jù)檢測及修正 方法，以確保網(wǎng)格化能夠正確剖分，并使得網(wǎng)格單元的精度能夠滿足地質工作者的需求。 [0065] 一、本發(fā)明的總體流程：
[0066] 地質體包括諸如地層、斷層、侵入巖、透鏡體、礦體復雜對象，具有十分復雜的空間 展布及其拓撲關系，缺乏使用數(shù)值計算解決問題的方法。本發(fā)明利用非數(shù)值計算方法，對虛 擬地質環(huán)境中的對象及其之間的相互關系加以描述。抽象的地質對象主要包括點、邊、面 片、表面和體，及其派生類的對象，如地層、斷層、侵入巖、透鏡體、礦體。
[0067] 在構建3D模型時，由于樣本的稀疏以及對象形狀的復雜性，需要在各個構建過程 中進行數(shù)據(jù)檢測，并指導數(shù)據(jù)及模型的修正。主要步驟如圖1所示。
[0068] 1、多源數(shù)據(jù)一體化：
[0069] 3D地質建模主要基于區(qū)域內（記為Ω )的多源數(shù)據(jù)，將Ω內的邊界數(shù)據(jù)、鉆孔數(shù)據(jù)、 剖面數(shù)據(jù)、3D掃描數(shù)據(jù)、以及各種解釋數(shù)據(jù)進行一體化處理，生成樣本數(shù)據(jù)集合S，作為地質 體網(wǎng)格化的基礎數(shù)據(jù)。
[0070] 2、布點檢測：
[0071] 通常情況下，Ω內的樣本數(shù)據(jù)比較缺乏，需要在其中布置數(shù)據(jù)點。采用2種方法，1) 特征點插入。識別Ω內的地質特征，如斷層、尖滅，并根據(jù)給定的密度閾值 ε，插入相應點，構 建點集合F; 2)規(guī)格化點插入。以Ω最小/大值為界，構建包圍盒，對包圍盒內進行規(guī)格化插 入點，生成點集合G。
[0072] 對于 Vfi&F & Vs:jeS:，其中，i = l，2, · · .m, j = l，2, · · .n〇，m為F集合的元素個數(shù)，η〇 為S的樣本初始個數(shù)。檢測F中的每個點f·1與VSj的距離，如果該距離大于 ε，則將h合并入S 中，即S = S U {fi}。最終獲得集合S，且其元素個數(shù)為m，no < m < m+ηο。
[0073] 對于VgkeG & VSjeS，其中，k=l，2, ...r, j = l，2, ...ηι，!·為 G 集合的元素個數(shù)。檢 測G中的每個點gk與Vsj的距離，如果該距離大于ε，則進行點 gk在地質體邊界區(qū)域內的判斷。 對于在Ω內的點gk，則將gk也合并入S中，即S = SU {gk}。
[0074] 3、表面模型檢測：
[0075] 以S為網(wǎng)格化的點數(shù)據(jù)，可采用三角形或四邊形作為面片，進行網(wǎng)格化剖分，構建 地質體的表面模型。為了保證表面模型的可靠性，需要進行下列檢測。
[0076] 1)約束檢測。對構成表面模型中的對象點、邊、及面片，分別進行冗余性和多余性 檢測。首先構建表面模型的拓撲結構（如圖2所示），即表面模型可以由多個面片構成;面片 由邊構成，或直接由點構成;邊由點構成。通過表面模型的拓撲結構，可以準確判斷它們之 間的關聯(lián)關系，在此基礎上，進行面片、邊、點的冗/多余性檢測。
[0077] 以邊冗/多余性檢測為例，設第k條邊ek，在拓撲結構中，首先搜索邊集合，如果發(fā) 現(xiàn)與邊改相同的邊e q，其中，9￡{1，2，...，《-{幻4為邊集合的個數(shù)，則發(fā)生冗余，刪除冗余 的邊eq，并修正與之關聯(lián)的拓撲信息，以確保信息的正確性與完整性。之后，再向上搜索％是 否構成某個或某幾個面片，如果搜索結果為NULL，則邊e k多余;否則假設說構成面片Cl，ie {1，2，...，M}，M為面片集合的個數(shù)，需要繼續(xù)檢測面片C 1的多余性；如果向上搜索結果為 NULL，則說明邊％也是多余的；最終將多余的邊％進行刪除處理。其它檢測操作類似。如果能 夠依次進行面片、邊、點的冗/多余性檢測，則各個對象的檢測僅需要在本層及向上搜索一 層即可基本完成。通過約束檢測，可以檢測并刪除圖3中的懸點p(圖3a)、懸邊e(圖3b)、以及 部分懸面片c (圖3c)情況。
[0078] 2)局部單元優(yōu)化檢測。在表面模型中遍歷搜索共享邊的面片單元，對兩個共享邊 的面片作最大空圓準則檢測，如果某個頂點在面片的外接圓之內，則修正對角線，完成單元 優(yōu)化過程的檢測。對于無法滿足最大空圓準則的，進行局部布點修正。
[0079] 3)閉合性檢測。為了真實地描述地質體，所構建的表面模型通常應該是一個閉合 的空間形狀，以便表達諸如侵入巖、透鏡體、礦體。為此需要對表面模型進行閉合性測試。從 拓撲結構中任意一條邊ei出發(fā)，向上搜索與邊ei共享的面片及其總數(shù)N，如果N! =2,則閉合 性檢測失敗，輸出不閉合的邊或面片；否則繼續(xù)下條邊的搜索，直至所有邊均搜索完畢，且 每條邊僅被兩個面片所共享，說明閉合性檢測成功，即表面模型具有閉合性。對于導致非閉 合的邊或面片，需要修正表面模型中的數(shù)據(jù)，并重新進行閉合性檢測。閉合性檢測也可以搜 索出并刪除懸面片（圖3c)情況。
[0080] 4)自交性檢測。在地質體網(wǎng)格化剖分過程中，由于地質現(xiàn)象的復雜性，以及剖分算 法的健壯性，可能導致網(wǎng)格發(fā)生自交。自交性檢測主要針對網(wǎng)格面片的自交性進行檢測。對 于Vc i面片，16{1，2，...^}，1為面片集合的個數(shù)，在所構建表面模型的拓撲結構中，通過 包圍盒檢測方法，搜索^，#{1，2，...，}&&_]_封；并分別計算(3 1和(^的平面方程;再將2個 平面方程求交，并進行交點有效性判斷，獲得交點系列*P。如果P = NULU則CdPCj不自交;否 則返回自交信息，例如圖4a穿透性自交、以及圖4b相鄰性自交。對于自交的網(wǎng)格，需要進行 網(wǎng)格重構，方法1是返回布點階段，重新插值及修正數(shù)據(jù);方法2是局部添加網(wǎng)格，如圖4c和 圖4d，分別對應圖4a和圖4b的Ci網(wǎng)格重構結果之一。
[0081] 5)流型檢測。為了進一步進行地質體的體模型構建，需要進行流型檢測，以確保體 剖分算法的穩(wěn)定性。對于懸點、懸邊、懸面片這類非流型，可以通過上述約束檢測、閉合性檢 測實現(xiàn)非流型的檢測并刪除相關數(shù)據(jù)。對于下列情況，即點、邊、面片都為模型的正確網(wǎng)格 剖分，但是由于多個網(wǎng)格模型的相鄰關系，使得存在某個邊被多于2個面片所共享，此時需 要進行如下檢測:設邊％上6{1，2，...，《4為邊集合的個數(shù)，在所構建表面模型的拓撲結 構中，向上搜索由其構成的面片，記為C 1面片，1[{1，2，...，}1為面片集合的個數(shù)兒為與 邊％共享的面片個數(shù)。如果L是奇數(shù)，則需要分解表面模型，直至L為偶數(shù)(如圖5)。
[0082] 6)面片方向性檢測。為了進行地質體的體剖分以及可視化渲染，需要對所構建表 面模型進行面片方向性檢測，以保持面片方向一致性并使其法向量指向表面模型外部。對 于凸面體和凹面體的表面模型，首先過面片內一點做垂直于面片的垂線，最少與模型的其 它表面面片有1個交點。當交點個數(shù)大于1時，從面片內一點開始的垂線，一端向著模型內部 方向延伸，與模型有奇數(shù)個交點；另一端向著模型外部方向延伸，與模型有偶數(shù)個交點。只 要選取向著模型內部方向延伸的這端距離面片最近的交點，就可以得到正確的面片特殊內 部點。之后以面片內一點為起點，以垂線與其它面片的交點為終點做一條向量，可以把與面 片法向量方向相同的向量對應的交點歸為一類，方向相反的歸為一類。判斷哪一類中具有 奇數(shù)個交點，從中選擇距離面片最近的交點來求特殊內部點。再從這個特殊內部點向面片 做垂線，根據(jù)垂線與面片法向量的方向關系做出調整，如果垂線的方向與面片法向量的方 向相同，則此面片法向量指向模型外部，否則交換任意兩個頂點的位置。
[0083] 4、體模型檢測：
[0084]對已構建地質體的表面模型進行上述檢測并修正完成之后，可以進行體剖分，生 成相應的四面體或六面體網(wǎng)格模型，并對其進行下列測試。
[0085] 1)網(wǎng)格優(yōu)化檢測。對于構建體模型的網(wǎng)格單元進行優(yōu)化檢測，優(yōu)化質量直接關系 模型的可用性與可靠性，如關系到礦體儲量計算分析的精度、工程應用問題。主要方法:計 算形成網(wǎng)格單元兩個面片的夾角，以檢測該網(wǎng)格單元的扭曲度;計算網(wǎng)格單元中最長邊與 最短邊之比，檢測其比值是否為1-3;計算單位雅可比值，檢測網(wǎng)格單元雅可比值是否在 [0.4,1]范圍內。此外，可以通過設置參數(shù)來控制生成網(wǎng)格單元的幾何屬性，從而提高網(wǎng)格 的質量。這些參數(shù)主要包括:局部網(wǎng)格尺寸、可接受的全局最大和最小網(wǎng)格尺寸、網(wǎng)格密度、 網(wǎng)格分級、元素曲率半徑、在封閉邊界上的網(wǎng)格精細化和三維網(wǎng)格優(yōu)化的步數(shù)。
[0086] 2)吻合度檢測。重點檢測體模型網(wǎng)格外邊界與表面模型邊界的吻合度，以確保體 模型能夠最好地逼近表面模型，從而真實反映地質體原貌。主要方法包括:面積計算，即分 別計算體模型網(wǎng)格外邊界單元的面積、以及表面模型單元的面積，檢測其差值是否趨近 〇.〇;體積計算，即分別計算體模型和表面模型的體積，并檢測其差值;共面檢測，這種檢測 一般適用于體模型和表面模型采用相同網(wǎng)格單元的情況，以四面體模型TEN和三角形表面 模型TIN為例，其網(wǎng)格單元都是3個點組成。從TIN的網(wǎng)格初始前沿開始（圖6a)，逐層向區(qū)域 內部推進，直至填滿整個區(qū)域（圖6b)。這里圖6a、6b為2D前沿推進法的基本過程，3D過程類 似。對于任意一個面片Vc^TIN rie {1，2,...，M}，M為TIN的面片集合的個數(shù)，通過遍歷相應 的TEN模型外邊界面片，則一定存在且僅存在一個面片3%eTEN，je {1，2,.. .，N}，N為TEN的 面片集合的個數(shù)，使得CdPuj-致;當檢測所有TIN中的面片與TEN模型外邊界面片中的面片 存在上述一致性關系時，返回完全吻合檢測信息。
[0087] 3)截面檢測。分別對體模型和表面模型進行截面計算與重構，并將它們進行一體 化顯示，通過檢測網(wǎng)格單元的面片、邊、點，進行可視化的對比分析。
[0088] 本發(fā)明實現(xiàn)的環(huán)境包括硬件環(huán)境和軟件環(huán)境：
[0089]硬件環(huán)境包括三維地震勘探設備、鉆孔、傳統(tǒng)測量設備、電子測量設備、GPS、激光 掃描儀、高性能PC、服務器。軟件環(huán)境包括地理信息系統(tǒng)、.net平臺和OpenGL圖形庫。
[0090]本發(fā)明的特點是：
[0091 ] I)本發(fā)明針對目前在網(wǎng)格化過程中缺乏有效而系統(tǒng)性的檢測手段和方法，通過對 地質體網(wǎng)格化過程中各個階段的分析，提出相應的數(shù)據(jù)檢測及修正方法。對3D地質建模區(qū) 域內的數(shù)據(jù)進行一體化處理，經(jīng)過布點檢測、表面模型檢測、體模型檢測，為地質體的網(wǎng)格 化提供準確的數(shù)據(jù)。
[0092] 2)通過構建表面模型的拓撲結構，對表面模型進行約束檢測、局部單元優(yōu)化檢測、 閉合性檢測、自交性檢測、流型檢測、面片方向性檢測一系列的系統(tǒng)性檢測，并提出指導性 的修正信息和方法，不僅保證表面模型的網(wǎng)格質量及其真實地刻畫地質體，而且為地質體 的體模型構建提供可剖分的數(shù)據(jù)。
[0093] 3)針對如何檢測復雜地質體剖分時的邊界一致性和生成高質量網(wǎng)格單元難題，本 發(fā)明提供網(wǎng)格優(yōu)化檢測、吻合度檢測、截面檢測的方法，可有效地對地質體中存在包含關 系、相鄰關系時的邊界一致性檢測及其修正，并使得網(wǎng)格單元的精度能夠滿足地質工作者 的需求。
[0094] 4)本發(fā)明所提供的地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法，不僅能夠檢測地質、礦 山、石油領域中復雜地質體網(wǎng)格化過程中的數(shù)據(jù)問題，而且也適用于醫(yī)學、生物領域中復雜 物體網(wǎng)格化過程中的數(shù)據(jù)檢測及其修正問題。
[0095] 二、具體實施例：
[0096]以位于河北一個鐵礦為例。
[0097] 具體步驟如下：
[0098]步驟1、將礦區(qū)內的邊界數(shù)據(jù)、鉆孔數(shù)據(jù)、剖面數(shù)據(jù)、3D掃描數(shù)據(jù)、以及各種解釋數(shù) 據(jù)進行一體化處理（圖7a)，生成各個地層的樣本數(shù)據(jù)集合(S1, S2，...，S9}。以如圖7b地層S4 的樣本數(shù)據(jù)為例，對其它樣本操作類似。它們將作為地質體網(wǎng)格化的基礎數(shù)據(jù)。
[0099]步驟2、從圖7a、圖7b可見，礦區(qū)內的樣本數(shù)據(jù)比較缺乏，需要在其中布置數(shù)據(jù)點。 首先可識別區(qū)內的斷層，并根據(jù)給定的密度閾值ε = 50,插入相應點，構成斷層F2的點集(圖 8a)。之后進行規(guī)格化點插入，設置諸如Χ/Υ/Ζ方向的密度、XY旋轉角度系列參數(shù)，獲得規(guī)格 化點集(圖8b)。將所生成的點合并入S 4中，以備網(wǎng)格化之用。
[0100] 步驟3、進行表面模型檢測，具體步驟如下：
[0101] 步驟301、以S4為網(wǎng)格化的點數(shù)據(jù)，進行網(wǎng)格化剖分，構建地層S4的表面模型，并生 成該表面模型的拓撲結構。如圖9表示S 4的拓撲結構在緩存中的存儲形式，分別顯示了構成 S4的其中第10個和第22個面片及其拓撲信息。
[0102] 步驟302、對構成表面模型S4中的對象點、邊、及面片，分別進行冗余性和多余性檢 測。在拓撲結構中，搜索構成S 4的面片集合，發(fā)現(xiàn)面片冗余（圖10a)，則可從S4的面片集合中 刪除此冗余面片，并修正與之關聯(lián)的拓撲信息;通過向上搜索可分別發(fā)現(xiàn)多余點（圖IOb)和 多余邊(圖l〇c)，并刪除相關的懸點和懸邊，最終可以獲得經(jīng)過約束檢測并修正后的S 4網(wǎng)格 (圖 10d)〇
[0103]步驟303、由于S4地層中包含了若干條斷層，并進行了特征點插入，在保持斷層特 征邊的前提下，對S4地層模型中的局部網(wǎng)格單元進行優(yōu)化檢測，對于無法滿足最大空圓準 則的，轉入步驟2，進行局部布點修正。
[0104]步驟304、對S4地層模型進行閉合性檢測，從拓撲結構中任意一條邊ei出發(fā)，向上搜 索與邊ei共享的面片及其總數(shù)N，如果N! =2,則閉合性檢測失敗，并輸出不閉合的相關邊 (圖11a)，轉入步驟301，對模型進行局部重構；否則繼續(xù)下條邊的搜索，直至所有邊均搜索 完畢，且每條邊僅被兩個面片所共享，則S4地層模型具有閉合性(圖I Ib)。
[0105] 步驟305、分別對S2、S4地層模型進行自交性檢測，如果檢測失敗，則轉入步驟301， 對模型進行局部修正與重構；如果檢測成功，而由于3 2是侵入巖且侵入到S4地層中，屬于地 質體中存在包含關系，則需要將2個地層模型合并，再進行自交性檢測。結果如圖12所示，網(wǎng) 格發(fā)生穿透性自交，則轉入步驟301，進行局部添加網(wǎng)格重構。
[0106] 步驟306、對S4地層模型實行流型檢測，在所構建S4模型的拓撲結構中，向上搜索并 發(fā)現(xiàn)存在邊被偶數(shù)個面片所共享(如圖13)，則不需要分解S 4模型。
[0107] 步驟307、對S4地層模型進行面片方向性檢測。遍歷S4模型所有表面面片，逐個進行 調整。先尋找面片對應的特殊內部點，再從這個特殊內部點向面片做垂線，根據(jù)垂線與面片 法向量的方向關系做出調整。具體步驟如下：
[0108] 步驟3071、選定面片內一點。假設模型的某個表面面片ABC的初始排列順序為A (xi，yi，zi)、B(X2，y2，Z2)、C(X3，y3，Z3)，取AB中點為D，取CD中點為E(x e，ye，ze)，選定面片ABC 內一點E。
[0109] 步驟3072、求過面片內一點且垂直于面片的直線。利用向量積求面片ABC的法向量 ηΞ (m，n，p) 〇
[0110] η 三 ABXBC (I)
[0111] 向量積可以通過行列式求解：
[0112]
[0113] 過點E且垂直于面片ABC的直線1^的參數(shù)方程如下，其中t為參數(shù)，因為直線L1平行 于面片ABC的法向量η，所以直線1^的方向向量可以取法向量η。
[0114]
[0115] 步驟3073、求直線與模型其它表面面片的所有交點。直線1^與除面片ABC外的其它 表面面片最少有一個交點，需要求出所有交點的坐標。依次求直線與面片所在平面的交點， 判斷交點與面片的位置關系，位置關系有兩大類:交點在面片內或外。取模型一表面面片 A 1B1C1為例說明直線與面片交點的求法。先求面片A1B 1C1所在平面的點法式方程。已知面片 AlBlCl三個頂點的坐標分別為Al = ( Xa，ya，Za )、Bi = ( Xb，yb，Zb )和Cl = ( Xc，yC，Zc )，根據(jù)式（1) 求得面片A1B1C1的法向量!!^(!^，!^，!^，又已知面片所在平面上一點此可得面片所在平 面的點法式方程如下。
[0116]
[0117] 把式⑵代入式⑶得
到參數(shù)t的值，從而求得交點丹二化^口以亦"的坐標。
[0118]
[0119]
[0120] 求得交點坐標后，判斷交點？:落在面片A1B1C1內還是外，分別連接交點PjP面片的 三個頂點，形成向量PiA 1、P1B1和P1C1，然后依次求兩兩向量的向量積P 1A1B1 = P1A1 X P1B1、 PiB1C1 = P1B1 X P1C1和P1C1A1 = P1C1 X P1A1。如果交點在面片內，則這三個向量積的方向相同， 如果交點在面片外，則有一個向量積的方向與其它兩個的方向剛好相反。利用如下的向量 積判斷PiA 1B1 A1B1C1和P1C1A 1三個向量的方向是否相同。
[0121]
[0122] 通過對交點與面片位置關系可能存在的所有情況的分析與歸納，共得到如表1所 示的五類情況。表1中條件指公式(4)中dotl、dot2和dot3計算結果的正負（正負的表示與 dotl、dot2和dot3的順序無關）。
[0123] 表1交點與面片的位置關系
[0125] 以上五類包括了所有可能存在的情況，可以判斷交點與面片的位置關系。當交點 落在面片的頂點或邊上時，一個交點可能被2個或2個以上的面片所共有，此時要注意只能 算做一個交點。
[0126] 步驟3074、選擇奇數(shù)端離面片距離最近的交點。做一條從面片上一點指向交點的 向量，根據(jù)此向量與面片法向量的關系，可以把交點分為兩類:與法向量方向相同的歸為一 類;方向相反的歸為一類。選擇包含奇數(shù)個交點的一類中離面片距離最近的交點。
[0127] 步驟3075、求面片的特殊內部點。以所選交點和面片上一點的中點作為面片的特 殊內部點。再調整面片前兩個頂點的順序，即已經(jīng)找到面片的特殊內部點，從特殊內部點向 面片做垂線，如果和面片法向量方向相同，無需調整;如果和面片法向量方向相反，則交換 面片任意兩個頂點的順序。
[0128] 步驟4、體模型檢測
[0129] 步驟401、對S4模型進行上述測試并修正完成之后，可以進行體剖分，生成相應的 四面體網(wǎng)格模型(如圖14)，并對其進行下列測試。
[0130] 步驟402、通過設置可接受的全局最大尺寸（le6)、網(wǎng)格密度（1)、網(wǎng)格分級（0)、啟 動在封閉邊界上的網(wǎng)格精細化和三維網(wǎng)格優(yōu)化的步數(shù)，來控制生成網(wǎng)格單元的幾何屬性， 對S 4的體模型進行網(wǎng)格優(yōu)化。
[0131] 步驟403、對地層SjPS2分別檢測體模型網(wǎng)格外邊界與表面模型邊界的吻合度;而 由于Si和S2為地質體中存在的相鄰關系，需要進一步進行模型合并之后的吻合度檢測。由表 2可知，檢測所有TIN中的面片與TEN模型外邊界面片中的面片存在一致性關系，返回完全吻 合檢測信息。
[0132]表2地層51和52吻合度檢測信息
[0134] 步驟404、分別對S4的體模型和表面模型進行截面計算與重構，如圖15a為任意一 個截面的計算結果，并將截面和S 4的表面模型進行一體化顯示(如圖15b)。為了進行區(qū)域整 體性的檢測，可以進行任意多截面的計算，如圖15c為分別對S 5、S6、S7的體模型和表面模型 進行任意多個截面及其一體化顯示，其中，截面的邊界用邊界點突出表達。之后，可以通過 檢測網(wǎng)格單元的面片、邊、點，進行可視化的對比分析，如圖15d是對圖15c中的局部放大結 果，以檢測截面與表面模型的數(shù)據(jù)一致性。
[0135] 步驟5、將構建并檢測通過的體模型和表面模型，進行工程應用，為解決地質復雜 現(xiàn)象提供可靠的平臺，使地質工作者可以獲得更多的地下信息。
[0136] 以上所述，僅為本發(fā)明較佳的【具體實施方式】，但本發(fā)明的保護范圍并不局限于此， 任何熟悉本技術領域的技術人員在本發(fā)明披露的技術范圍內，可輕易想到的變化或替換， 都應涵蓋在本發(fā)明的保護范圍之內。因此，本發(fā)明的保護范圍應該以權利要求書的保護范 圍為準。
【主權項】
1. 一種地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法，其特征在于，利用非數(shù)值計算方法，對虛 擬地質環(huán)境中的對象及其之間的相互關系加 W描述并構建3D模型，抽象的地質對象主要包 括點、邊、面片、表面和體及其派生類的對象，具體包括地層、斷層、侵入巖、透鏡體、礦體，在 構建3D模型時，需要在各個構建過程中進行數(shù)據(jù)檢測，并指導數(shù)據(jù)及模型的修正，主要包括 步驟： A、 多源數(shù)據(jù)一體化： 3D地質建模主要基于區(qū)域內的多源數(shù)據(jù)，該區(qū)域記為Ω，將Ω內的邊界數(shù)據(jù)、鉆孔數(shù) 據(jù)、剖面數(shù)據(jù)、3D掃描數(shù)據(jù)、W及各種解釋數(shù)據(jù)進行一體化處理，生成樣本數(shù)據(jù)集合S，作為 地質體網(wǎng)格化的基礎數(shù)據(jù)； B、 布點檢測： 由于Ω內的樣本數(shù)據(jù)比較缺乏，需要在其中布置數(shù)據(jù)點，采用巧巾方法： 方法一、特征點插入:識別Ω內的地質特征，如斷層、尖滅，并根據(jù)給定的密度闊值ε，插 入相應點，構建點集合F; 方法二、規(guī)格化點插入：ΚΩ最小/大值為界，構建包圍盒，對包圍盒內進行規(guī)格化插入 點，生成點集合G; C、 表面模型檢測： WS為網(wǎng)格化的點數(shù)據(jù)，采用Ξ角形或四邊形作為面片，進行網(wǎng)格化剖分，構建地質體 的表面模型，并進行下列檢測： 約束檢測、局部單元優(yōu)化檢測、閉合性檢測、自交性檢測、流型檢測、面片方向性檢測； D、 體模型檢測： 對已構建地質體的表面模型進行上述檢測并修正完成之后，可W進行體剖分，生成相 應的四面體或六面體網(wǎng)格模型，并對其進行下列測試： 網(wǎng)格優(yōu)化檢測、吻合度檢測、截面檢測。2. 根據(jù)權利要求1所述的地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法，其特征在于，所述步驟 Β中： 對于Vf.i.eF色Vsye-S，.其中，i = l ,2,. . .m, j = l，2,. . .n〇，m為F集合的元素個數(shù)，η〇為S的 樣本初始個數(shù)，檢測F中的每個點fi與Vs,的距離，如果該距離大于ε，則將fi合并入S中，即S =S U也}，最終獲得集合S，且其元素個數(shù)為m，no如1如+n〇; 對于V掛€狂底Vs,eS，其中，k = l，2, . . .r, j = l，2, . . .ηι,Γ為G集合的元素個數(shù)，檢測G 中的每個點gk與V成的距離，如果該距離大于ε，則進行點gk在地質體邊界區(qū)域內的判斷，對 于在Ω內的點阱，則將阱也合并入S中，即S = SU {gk}。3. 根據(jù)權利要求2所述的地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法，其特征在于，所述步驟 帥： 所述約束檢測包括:對構成表面模型中的對象點、邊、及面片，分別進行冗余性和多余 性檢測，首先構建表面模型的拓撲結構，即表面模型可W由多個面片構成;面片由邊構成， 或直接由點構成；邊由點構成，通過表面模型的拓撲結構，能準確判斷它們之間的關聯(lián)關 系，在此基礎上，進行面片、邊、點的冗/多余性檢測； 所述局部單元優(yōu)化檢測包括:在表面模型中遍歷捜索共享邊的面片單元，對兩個共享 邊的面片作最大空圓準則檢測，如果某個頂點在面片的外接圓之內，則修正對角線，完成單 元優(yōu)化過程的檢測，對于無法滿足最大空圓準則的，進行局部布點修正； 所述閉合性檢測包括:從拓撲結構中任意一條邊ei出發(fā)，向上捜索與邊ei共享的面片及 其總數(shù)N，如果N! =2,則閉合性檢測失敗，輸出不閉合的邊或面片，否則繼續(xù)下條邊的捜索， 直至所有邊均捜索完畢，且每條邊僅被兩個面片所共享，說明閉合性檢測成功，即表面模型 具有閉合性，對于導致非閉合的邊或面片，需要修正表面模型中的數(shù)據(jù)，并重新進行閉合性 檢測； 所述自交性檢測主要針對網(wǎng)格面片的自交性進行檢測，包括：對于V&面片，i E {1， 2, ...，M}，M為面片集合的個數(shù)，在所構建表面模型的拓撲結構中，通過包圍盒檢測方法，捜 索cj，j e {1，2，.. .，M}&&扣i;并分別計算Cl和C撕平面方程;再將2個平面方程求交，并進 行交點有效性判斷，獲得交點系列*P，如果P = NULL，則Cl和cj不自交，否則返回自交信息，對 于自交的網(wǎng)格，按W下兩種方法進行網(wǎng)格重構，方法一是返回布點階段，重新插值及修正數(shù) 據(jù)，方法二是局部添加網(wǎng)格； 所述流型檢測包括:對于下列情況，即點、邊、面片都為模型的正確網(wǎng)格剖分，但是由于 多個網(wǎng)格模型的相鄰關系，使得存在某個邊被多于2個面片所共享，此時需要進行如下檢 巧。:設邊ek，ke{l，2，...，N}，N為邊集合的個數(shù)，在所構建表面模型的拓撲結構中，向上捜 索由其構成的面片，記為Cl面片，1￡{1，2，...，1}，1為面片集合的個數(shù)，1^為與邊61^共享的面 片個數(shù)，如果L是奇數(shù)，則需要分解表面模型，直至L為偶數(shù)； 所述面片方向性檢測包括:對于凸面體和凹面體的表面模型，首先過面片內一點做垂 直于面片的垂線，最少與模型的其它表面面片有1個交點，當交點個數(shù)大于1時，從面片內一 點開始的垂線，一端向著模型內部方向延伸，與模型有奇數(shù)個交點；另一端向著模型外部方 向延伸，與模型有偶數(shù)個交點，只要選取向著模型內部方向延伸的運端距離面片最近的交 點，就能得到正確的面片特殊內部點，之后W面片內一點為起點，W垂線與其它面片的交點 為終點做一條向量，把與面片法向量方向相同的向量對應的交點歸為一類，方向相反的歸 為一類，判斷哪一類中具有奇數(shù)個交點，從中選擇距離面片最近的交點來求特殊內部點，再 從運個特殊內部點向面片做垂線，根據(jù)垂線與面片法向量的方向關系做出調整，如果垂線 的方向與面片法向量的方向相同，則此面片法向量指向模型外部，否則交換任意兩個頂點 的位置。4.根據(jù)權利要求3所述的地質體網(wǎng)格化過程中的數(shù)據(jù)檢測方法，其特征在于，所述步驟 帥： 所述網(wǎng)格優(yōu)化檢測包括:計算形成網(wǎng)格單元兩個面片的夾角，W檢測該網(wǎng)格單元的扭 曲度;計算網(wǎng)格單元中最長邊與最短邊之比，檢測其比值是否為1-3;計算單位雅可比值，檢 測網(wǎng)格單元雅可比值是否在[0.4,1]范圍內，通過設置參數(shù)來控制生成網(wǎng)格單元的幾何屬 性，從而提高網(wǎng)格的質量，運些參數(shù)主要包括：局部網(wǎng)格尺寸、可接受的全局最大和最小網(wǎng) 格尺寸、網(wǎng)格密度、網(wǎng)格分級、元素曲率半徑、在封閉邊界上的網(wǎng)格精細化和Ξ維網(wǎng)格優(yōu)化 的步數(shù)； 所述吻合度檢測包括:面積計算，即分別計算體模型網(wǎng)格外邊界單元的面積、W及表面 模型單元的面積，檢測其差值是否趨近0.0;體積計算，即分別計算體模型和表面模型的體 積，并檢測其差值;共面檢測，運種檢測一般適用于體模型和表面模型采用相同網(wǎng)格單元的 情況，從TIN的網(wǎng)格初始前沿開始，逐層向區(qū)域內部推進，直至填滿整個區(qū)域，對于任意一個 面片VCienN，i e {1，2，. ..，M}，M為TIN的面片集合的個數(shù)，通過遍歷相應的TE飾莫型外邊界 面片，則一定存在且僅存在一個面片3U拒征N，je{l，2,. . .，N}，N為TEN的面片集合的個 數(shù)，使得Cl和W-致;當檢測所有TIN中的面片與TEN模型外邊界面片中的面片存在上述一致 性關系時，返回完全吻合檢測信息； 所述截面檢測包括:分別對體模型和表面模型進行截面計算與重構，并將它們進行一 體化顯示，通過檢測網(wǎng)格單元的面片、邊、點，進行可視化的對比分析。
【文檔編號】G06T17/20GK105844710SQ201610221833
【公開日】2016年8月10日
【申請日】2016年4月11日
【發(fā)明人】武強, 徐華, 魏斌
【申請人】北京石油化工學院, 中國礦業(yè)大學(北京)