大范圍三維地形動(dòng)態(tài)可視化技術(shù)研究

董興鵬 陳偉 張麗 鄭雄偉 黃煒

摘要：三維地形以數(shù)字化的方式對(duì)地形表面起伏狀態(tài)和細(xì)節(jié)特征進(jìn)行模擬，廣泛應(yīng)用于飛行視景仿真、戰(zhàn)場(chǎng)態(tài)勢(shì)仿真、模擬演練、應(yīng)急指揮、災(zāi)害救援等應(yīng)用場(chǎng)景中。三維地形的真實(shí)性和可視化渲染的實(shí)時(shí)性直接決定了三維應(yīng)用場(chǎng)景的逼真度、沉浸感和有效性。文主要針對(duì)地質(zhì)遙感數(shù)據(jù)成果可視化及共享展示，解決地質(zhì)遙感數(shù)據(jù)體三維展示效果不佳，地質(zhì)調(diào)查成果數(shù)據(jù)展示組織方式欠缺，可視化系統(tǒng)性能欠佳的問題，提供了優(yōu)化的數(shù)據(jù)存儲(chǔ)渲染技術(shù)、多維一體化地質(zhì)大數(shù)據(jù)成果組織方式、綜合可視化承載平臺(tái)以及數(shù)據(jù)關(guān)聯(lián)分析功能，為地質(zhì)調(diào)查應(yīng)用系統(tǒng)中的遙感地質(zhì)可視化分系統(tǒng)提供技術(shù)模塊。

關(guān)鍵詞：三維展示，關(guān)聯(lián)分析，模擬演練

Abstract：3D terrain is used to simulate the undulating state and detailed features of terrain surface in a digital way，which is widely used in flight scene simulation，battlefield situation simulation，simulation exercise，emergency command，disaster rescue and other application scenarios. The reality of 3D terrain and the real-time of visual rendering directly determine the fidelity，immersion and effectiveness of 3D application scene. This paper mainly aims at the visualization and sharing display of geological remote sensing data，solves the problems of poor 3D display effect of geological remote sensing data volume，lack of organization mode of geological survey data display，and poor performance of visualization system，and provides optimized data storage and rendering technology，multi-dimensional integrated geological big data achievement organization mode，comprehensive visualization bearing platform and data Customs Combined analysis function provides technical modules for remote sensing geological visualization subsystem in geological survey application system.

Key Words：3D display;association analysis;Simulated drill

1前言

地球表面是人類社會(huì)賴以生存并從事一切實(shí)踐活動(dòng)的根基，人類的生產(chǎn)生活與其發(fā)生著緊密的聯(lián)系（劉學(xué)，2011），地形與地貌是表達(dá)地表的最重要自然地理要素。地理學(xué)、測(cè)繪學(xué)和地貌學(xué)等學(xué)科中一直把對(duì)地形的科學(xué)表達(dá)作為研究重點(diǎn)和熱點(diǎn)，隨著計(jì)算機(jī)圖形學(xué)的發(fā)展，將物理世界的三維特征用計(jì)算機(jī)進(jìn)行真實(shí)的再現(xiàn)已成為現(xiàn)實(shí)（朱永麗，2012），三維地形構(gòu)建技術(shù)應(yīng)運(yùn)而生?；诟呔群透叻直媛实牡匦慰臻g渲染、更細(xì)膩的地表細(xì)節(jié)呈現(xiàn)、更絢麗的地上景觀模擬、更豐富的光照渲染效果、更便捷流暢的交互操作成為虛擬現(xiàn)實(shí)三維場(chǎng)景普遍需求（郭浩然，2013），然而這些需求又反向推動(dòng)地形三維重建的發(fā)展。目前，三維地形重建以及基于地形的可視化應(yīng)用面臨挑戰(zhàn)。大范圍三維地形場(chǎng)景的真實(shí)感表達(dá)與呈現(xiàn)為三維場(chǎng)景可視化應(yīng)用處理的實(shí)時(shí)性帶來挑戰(zhàn)。在大規(guī)模的虛擬現(xiàn)實(shí)三維應(yīng)用場(chǎng)景中，準(zhǔn)確的繪制地形的起伏特征、逼真的渲染出地表要素都與場(chǎng)景的實(shí)時(shí)快速渲染相互矛盾，尋求既能滿足大規(guī)模地形渲染真實(shí)感的要求又能滿足實(shí)時(shí)可視化需求的方法，值得深入研究。

2基于Delaunay準(zhǔn)則和生長(zhǎng)法的網(wǎng)格化方法

在分析了三維空間離散點(diǎn)云網(wǎng)格化方法的基礎(chǔ)上，結(jié)合無人機(jī)序列圖像生成地面點(diǎn)云數(shù)據(jù)的特點(diǎn)，本文采用生長(zhǎng)法依據(jù) Delaunay準(zhǔn)則對(duì)獲取的地面點(diǎn)云數(shù)據(jù)進(jìn)行網(wǎng)格化。

2.1構(gòu)建可伸縮包圍盒

地形點(diǎn)云數(shù)據(jù)與其他點(diǎn)云數(shù)據(jù)相比，具有數(shù)據(jù)量大、不封閉、空間范圍廣的特點(diǎn)，采用傳統(tǒng)等大小包圍盒的剖分方法對(duì)地形點(diǎn)云數(shù)據(jù)進(jìn)行剖分會(huì)產(chǎn)生許多空包圍盒，造成空間的浪費(fèi)和搜索效率的下降。同時(shí)，考慮到地形點(diǎn)云數(shù)據(jù)的高程特性和連續(xù)性，點(diǎn)云在水平面上投影鄰近的點(diǎn)在空間中也是鄰近的，所以本文提出沿地形高程方向可伸縮的包圍盒對(duì)地形點(diǎn)云數(shù)據(jù)進(jìn)行空間剖分。

可伸縮的包圍盒是一種長(zhǎng)方體包圍盒，該包圍盒垂直于高程方向的一組面的大小固定，而沿著高程方向的長(zhǎng)度可以自由伸縮，以便包圍所有投影在該組平面內(nèi)的點(diǎn)。假設(shè)xoy為水平面，z軸為地形的高程方向，圖2-1投影區(qū)域?yàn)?的可伸縮包圍盒，包圍盒固定大小的一組面與xoy 平面平行，其沿z軸方向的長(zhǎng)度不固定。此時(shí)整個(gè)空間剖分的包圍盒個(gè)數(shù)只與點(diǎn)云數(shù)據(jù)所在平面投影區(qū)域面積的大小成正比，極大地減少了包圍盒的數(shù)量。

為了確定點(diǎn)云在xoy投影區(qū)域的大小，在讀入點(diǎn)云數(shù)據(jù)的同時(shí)，記錄下點(diǎn)云的x、y坐標(biāo)的最大值和最小值 Xmax、Xmin、Ymax、Ymin。

為了確定小正方形的邊長(zhǎng)，需要計(jì)算平均每個(gè)投影點(diǎn)的投影面積s，s按照下方公式計(jì)算，其中N為點(diǎn)云中點(diǎn)的總數(shù)量。

每個(gè)小正方形的邊長(zhǎng)e按照公式計(jì)算，m為平均每個(gè)正方形內(nèi)的投影點(diǎn)個(gè)數(shù)，m為整數(shù)（本文中m=5），s為公式的結(jié)果。

2.2計(jì)算空間鄰近點(diǎn)

可伸縮包圍盒的構(gòu)建極大地縮小了一個(gè)點(diǎn)的空間鄰近點(diǎn)的范圍。如圖2-2xoy投影平面的剖分圖所示，包圍盒6中點(diǎn)的空間鄰近點(diǎn)一定在包圍盒6和與其相鄰的8個(gè)包圍盒共同組成的空間中。若采用空間點(diǎn)鄰近的全部包圍盒中的點(diǎn)來定義該點(diǎn)的空間鄰近點(diǎn)范圍時(shí)，當(dāng)空間點(diǎn)位于自身所在包圍盒的邊緣時(shí)，存在較大的誤差無法準(zhǔn)確的描述其鄰近點(diǎn)。本文在包圍盒的基礎(chǔ)上通過與可變半徑的包圍球求交運(yùn)算來確定空間鄰近點(diǎn)。若要求任意一空間點(diǎn)P的鄰近點(diǎn)，定義V為其鄰近點(diǎn)的集合，初始為空，首先計(jì)算點(diǎn)P自身包圍盒編號(hào)Pnum以及鄰近的8各包圍盒編號(hào)，然后構(gòu)造球心為P的包圍球Q，將所有包圍盒與包圍球Q求交，將獲得的點(diǎn)加入到集合V中，如為空則增大包圍球的半徑，直至V非空或包圍球半徑大于給定閾值為止。

2.3種子三角形的確定

本文基于生長(zhǎng)法生成地形網(wǎng)格，首先需要選取三個(gè)頂點(diǎn)構(gòu)成一個(gè)種子三角形，然后依據(jù)生長(zhǎng)算法進(jìn)行網(wǎng)格的生成。同時(shí)為了后續(xù)網(wǎng)格繪制和顯示的需要，網(wǎng)格生成時(shí)按照右手法則給出生成三角面的正面，即三角形中各頂點(diǎn)的順序由右手法則確定（陳思環(huán)，2018）。圖2-2是對(duì)投影平面上的6個(gè)點(diǎn)A、B、C、D、E、F的三角網(wǎng)格化示意過程，其中△ABC 是種子三角形如圖2-2 a所示，依據(jù)右手法則確定種子三角形頂點(diǎn)的順序?yàn)?A->B->C，從而確定了種子三角形的正面;然后按照△ABC 邊的生成順序進(jìn)行搜索最鄰近點(diǎn)，依次生成△BAD、△CBE 和△CAF。三角形的生成過程如圖：

種子三角形頂點(diǎn)的確定：本文中基于生長(zhǎng)法的點(diǎn)云三角網(wǎng)網(wǎng)格化過程中三角形的生成是采用廣度優(yōu)先算法來對(duì)種子三角的邊進(jìn)行擴(kuò)展生成新三角形。種子三角形的確定對(duì)點(diǎn)云網(wǎng)格化結(jié)果影響較小，所以種子三角形可以依據(jù)實(shí)際需要指定其初始位置。本文在進(jìn)行種子三角形的第一點(diǎn)選擇時(shí)，選取所有點(diǎn)云中坐標(biāo)最小點(diǎn)作為第一個(gè)點(diǎn)，依據(jù)圖中空間鄰近點(diǎn)的計(jì)算算法，求出p1的鄰近點(diǎn)集合v1，在集合v1中取出最近鄰點(diǎn)p2，計(jì)算出p2的鄰近點(diǎn)集合v2，求集合v1和集合v2的交集v3，在求出距離上面兩個(gè)點(diǎn)距離和最小的點(diǎn)p3，同時(shí)要求點(diǎn)p3與p1、p2生成的三角形為銳角三角形，該三角形即為種子三角形。

種子三角形方向的確定：在對(duì)點(diǎn)云網(wǎng)格進(jìn)行渲染時(shí)，需要確定網(wǎng)格三角形的正面，剔除網(wǎng)格的背面來加快網(wǎng)格化渲染的速度。為了保證生成的網(wǎng)格中所有三角形的正面是目標(biāo)區(qū)域地形表面，本文采用右手法則給出生成三角面的正面，即三角形中各頂點(diǎn)的順序由右手法則確定。在所生成的三角網(wǎng)格中相鄰的兩個(gè)三角形（有2個(gè)共同的頂點(diǎn)）的公共邊方向相反，在種子三角形方向已知時(shí)，后續(xù)擴(kuò)展的三角形頂點(diǎn)順序按照傳遞性來獲取。所以在地形網(wǎng)格生成時(shí)，為了保證生成的網(wǎng)格的方向同地形的外表面一致，只需要使得初始的種子三角形p1p2p3的正面與地形的外表面一致即可。

地形有三種類型的地勢(shì)，平地、低谷和山峰，如圖2-3地形截面圖，圖中粗線以上部分是地形外表面，以下是地形內(nèi)表面。若地形外表面使用逆時(shí)針方向，則外法線方向與空間坐標(biāo)系中的Z軸所形成的夾角小于 90°，因此通過夾角來判斷生成三角形頂點(diǎn)的順序。以種子三角形為例若cosβ ≥ 0時(shí)，三角形的頂點(diǎn)順序是p1p2p3，反之三角形的頂點(diǎn)順序是p1p3p2。

2.4廣度優(yōu)先搜索過程

在確定了種子三角形后，需要按照生長(zhǎng)法來擴(kuò)展生成新三角形，最終完成網(wǎng)格的構(gòu)建，在新三角形生成過程中采用廣度優(yōu)先的搜索算法來對(duì)種子三角的邊進(jìn)行新三角形的擴(kuò)展。為了便于對(duì)三角形生成過程的描述，對(duì)網(wǎng)格生成過程中（圖2-4）的點(diǎn)和邊進(jìn)行如下定義：內(nèi)邊：鄰接三角形的公共邊。如邊AB（BA）、BC（CB）、AC（CA）、AD（DA）、BD（DB）、BH（HB）、BI（IB）、BE（EB）。

外邊：未生成鄰接三角形的邊。如邊AG、GD、DH、HI、IE、EC、CF、FA。

內(nèi)點(diǎn)：以其為頂點(diǎn)的邊都是內(nèi)邊的點(diǎn)稱為內(nèi)點(diǎn)。如點(diǎn)B。

外點(diǎn)：以其為頂點(diǎn)的邊中含有外邊的點(diǎn)稱為外點(diǎn)。如點(diǎn)A、G、D、H、I、E、C、F。

待連接點(diǎn)：不是任何三角形的頂點(diǎn)的點(diǎn)稱為待連接點(diǎn)。如點(diǎn) J。

圖中種子三角形是△ABC，當(dāng)種子三角形選定后，依據(jù)右手法則確定種子三角形的正面，即種子三角形頂點(diǎn)的順序，按照種子三角形的頂點(diǎn)順序依次生成種子三角形的三條邊。在種子三角形的基礎(chǔ)上，依照廣度優(yōu)先的搜索策略從外邊中選擇一條邊進(jìn)行新三角形的擴(kuò)展。如圖 2-4中種子三角形的外邊生成順序是AB->BC->CA，在進(jìn)行新三角形擴(kuò)展時(shí)，首先選擇外邊 AB 作為新三角形的一條邊，查找構(gòu)成新三角形的另一頂點(diǎn) D，生成外邊 AB 的鄰接三角形BAD，同時(shí)AB變?yōu)閮?nèi)邊。新生成的△BAD 的邊AD和DB成為了外邊。這時(shí)生成的網(wǎng)格中的外邊有 BC、CA、AD 和 DB，按照廣度優(yōu)先的搜索策略，外邊按照 BC->CA ->AD->DB的順序擴(kuò)展新的鄰接三角形，分別是：△CBE、△ACF、△DAG、△BDH。

2.5最優(yōu)點(diǎn)的候選點(diǎn)集確定

基于生長(zhǎng)法進(jìn)行點(diǎn)云數(shù)據(jù)網(wǎng)格化構(gòu)建的過程是一個(gè)不停的為已知的外邊選擇第三個(gè)頂點(diǎn)生成新三角形的過程，第三點(diǎn)的搜索集合稱為外邊的候選點(diǎn)集。在保證網(wǎng)格構(gòu)建效果的前提條件下，為了加快網(wǎng)格的構(gòu)建，應(yīng)盡可能的縮小第三點(diǎn)的搜索范圍即縮小候選點(diǎn)集，因此對(duì)第三點(diǎn)做如下限制來縮小候選點(diǎn)集。假設(shè)當(dāng)前需要進(jìn)行三角形擴(kuò)展的外邊是 MN，V 是外邊 MN 的候選點(diǎn)集，V 內(nèi)的點(diǎn)需滿足如下要求：

⑴最大化最小角特性限制。三角形中最小角最大，即相鄰兩三角形組成的凸四邊形互換對(duì)角線，六個(gè)內(nèi)角中的最小角沒有增大。為了使新擴(kuò)展的三角形形態(tài)更優(yōu)，防止生成窄帶三角形（如圖 2-5），新生成的三角形中最小內(nèi)角不小于一個(gè)下限閾值，在本文中的值為π6。

⑵二面角限制。為使重建的地形表面連續(xù)且平滑，假設(shè)外邊MN所在的三角形為△MQN，頂點(diǎn)P與外邊 MN 生成的新三角形為△MNP，△MQN 與△MNP之間的二面角越接近π2越好。在網(wǎng)格構(gòu)建時(shí)已經(jīng)對(duì)三角形的正面進(jìn)行了一致化，為了保證生成的網(wǎng)格的方向同地形的外表面相同，即網(wǎng)格中全部三角形對(duì)應(yīng)的法線方向都由內(nèi)表面指向外表面。因此，可以通過限制△MQN和△MNP法線的夾角來實(shí)現(xiàn)△MQN與△MNP之間的二面角的限制，夾角越小，兩個(gè)三角形的二面角就越趨近π2，如圖2-6。

通過夾角來判斷兩鄰接三角形的平滑性，以邊 MN 的兩個(gè)鄰接三角形為例，本文要求的余弦值需要滿cosβ>0，即的最大上限閾值為2。

⑶非法點(diǎn)限制。第三個(gè)頂點(diǎn)P不能是已生成網(wǎng)格的內(nèi)點(diǎn)、點(diǎn)P不能是已有三角形 MQN 的頂點(diǎn)、點(diǎn)P不與MN共線。

⑷邊限制。邊MP與NP不為內(nèi)邊。

以外邊 MN 為例，其候選點(diǎn)集合V可采用如下算法進(jìn)行確定：第一步：在外邊MN兩個(gè)頂點(diǎn)的鄰接外邊的頂點(diǎn)中進(jìn)行檢索，盡快完成外邊MN的三角拓展;若檢索到的鄰接點(diǎn)滿足上述（1）、（2）、（3）、（4）限制條件，判定該點(diǎn)為候選點(diǎn)加入V中，若集合V為空轉(zhuǎn)入下一步，反之結(jié)束搜索，圖 2-7表示了該步驟的搜索過程，從圖2-7a中可知確定點(diǎn)J和點(diǎn)K為符合條件的頂點(diǎn)，由于頂點(diǎn)J是已有三角形的點(diǎn)判定為非法點(diǎn)，點(diǎn)K為滿足條件的候選點(diǎn)，與外邊MN生成新三角形MNK。

第二步：對(duì)兩個(gè)頂點(diǎn) M、N 所在包圍盒內(nèi)的其他點(diǎn)進(jìn)行判斷，若滿足上述（1）、（2）、（3）、（4）限制條件，判定該點(diǎn)為候選點(diǎn)加入 V 中，待所有點(diǎn)判斷完成后，若集合V為空轉(zhuǎn)入下一步，反之結(jié)束搜索。

第三步：按照鄰近點(diǎn)計(jì)算方法，分別求外邊頂點(diǎn)M、N的鄰近點(diǎn)以M、N為圓心，初始半徑r=e時(shí)的鄰近點(diǎn)集（不包括已經(jīng)搜索過的點(diǎn)），對(duì)檢索到的所有點(diǎn)進(jìn)行判斷，若滿足（1）、（2）、（3）、（4）的限制條件，判定該點(diǎn)為候選點(diǎn)加入V中，待所有點(diǎn)判斷完成后，若集合V為空增大鄰域半徑r，擴(kuò)大搜索范圍，繼續(xù)第三步，直至V非空或r達(dá)到閾值R。

閾值R是針對(duì)在地形邊緣網(wǎng)格生成時(shí)給定的限定條件，避免無意義的查找。半徑r增大的步距和上限閾值R可通過實(shí)驗(yàn)獲得，在本文中依據(jù)本章后續(xù)實(shí)驗(yàn)，設(shè)r的步距e，R=2e。上述搜索過程結(jié)束時(shí)，若候選點(diǎn)集V為空，則將MN記為邊界邊;若候選點(diǎn)集V中只有一個(gè)唯一點(diǎn)，則該點(diǎn)就是MN構(gòu)成新三角形的最優(yōu)點(diǎn);若候選點(diǎn)集V中存在多個(gè)候選點(diǎn)，需要對(duì)點(diǎn)集中的候選點(diǎn)進(jìn)行量化評(píng)價(jià)，選出最優(yōu)點(diǎn)做作為生成新三角形的最優(yōu)候選點(diǎn)。

2.6候選點(diǎn)的評(píng)價(jià)函數(shù)

通過上節(jié)的算法可以求出外邊 MN 的候選點(diǎn)集 V，當(dāng)候選點(diǎn)集 V 中含有多個(gè)候選點(diǎn)時(shí)，需要定義一個(gè)評(píng)價(jià)函數(shù)對(duì)集合 V 中的點(diǎn)進(jìn)行量化評(píng)價(jià)，依據(jù)評(píng)價(jià)結(jié)果確定最優(yōu)點(diǎn)做作為生成新三角形的最優(yōu)候選點(diǎn)。多數(shù)評(píng)價(jià)函數(shù)從新生成的三角形是否接近 Delaunay 三角形和地形網(wǎng)格的平滑度方面去考慮，而忽略了候選點(diǎn)到外邊的距離對(duì)整個(gè)網(wǎng)格的影響。圖 2-8 a 中 MN 是外邊，點(diǎn) P1和點(diǎn) P2是外邊 MN 候選點(diǎn)集合中的兩個(gè)候選點(diǎn)，點(diǎn) P2與外邊 MN 生成的三角形△MNP2與其鄰接的△MQN 形成的二面角更接近 180°?！鱉NP2的最小內(nèi)角比△MNP1大。

若僅考慮地形平滑度和 Delaunay 特性的限制，點(diǎn) P2是生成新三角形的最優(yōu)點(diǎn)，生成的新三角形圖

圖2-8 b所示，而點(diǎn) P1由于不滿足相鄰三角形的二面角的限制，不能與其他任何外邊形成三角形而成為噪聲點(diǎn)。若果選擇候選點(diǎn)P1進(jìn)行新三角形的生成，可生成△MNP1同時(shí)點(diǎn)P2在后續(xù)網(wǎng)格化中可生成△MP1P2，如圖 2-8 c 所示，此時(shí)點(diǎn)P1和P2都是有用點(diǎn)。出現(xiàn)圖2-8 b這一現(xiàn)象的原因是候選點(diǎn)P1與外邊MN的距離更近，且候選點(diǎn)P1在平面MNP2上的投影點(diǎn)在△MNP2內(nèi)部，顯然 MP1 NP1 MP2NP2 。在網(wǎng)格構(gòu)建時(shí)為了避免類似情況的發(fā)生，在最優(yōu)候選點(diǎn)選擇時(shí)，應(yīng)綜合考慮地形平滑度、Delaunay 三角形最小內(nèi)角和候選點(diǎn)到外邊的距離。

3.總結(jié)

三維地形以數(shù)字化的方式對(duì)地形表面起伏狀態(tài)和細(xì)節(jié)特征進(jìn)行模擬，廣泛應(yīng)用于飛行視景仿真、戰(zhàn)場(chǎng)態(tài)勢(shì)仿真、模擬演練、應(yīng)急指揮、災(zāi)害救援等應(yīng)用場(chǎng)景中。以無人機(jī)采集的序列圖像進(jìn)行三維地形重建相比傳統(tǒng)方式進(jìn)行大范圍三維重建在成本、效率、數(shù)據(jù)現(xiàn)勢(shì)性等方面都具有巨大優(yōu)勢(shì)。本文圍繞著基于無人機(jī)采集到的序列圖像進(jìn)行大規(guī)模地形三維重建的關(guān)鍵技術(shù)展開研究，分別從序列圖像的特征提取、圖像匹配、稀疏重建、稠密重建、離散點(diǎn)云網(wǎng)格化、地形數(shù)據(jù)可視化、大規(guī)模地形實(shí)時(shí)渲染等多個(gè)方面展開，對(duì)重建過程中的關(guān)鍵問題予以解決，給出了具體的解決辦法，概括起來本文的創(chuàng)新性工作體現(xiàn)在以下方面：

針對(duì)離散點(diǎn)云數(shù)據(jù)的三維重建網(wǎng)格化方法進(jìn)行研究，結(jié)合無人機(jī)序列圖像獲取到的地形點(diǎn)云數(shù)據(jù)的特點(diǎn)，提出了一種基于 Delaunay 準(zhǔn)則和生長(zhǎng)法的網(wǎng)格化方法來完成離散地形點(diǎn)云數(shù)據(jù)的快速網(wǎng)格化方法。該方法首先采用可伸縮包圍盒對(duì)空間進(jìn)行剖分，求點(diǎn)的空間鄰近點(diǎn)（k-鄰近點(diǎn)）;然后從點(diǎn)云中選取一個(gè)帶方向的種子三角形，并設(shè)計(jì)點(diǎn)和邊的搜索策略;最后以已有的外邊為基礎(chǔ)，確定外邊的候選點(diǎn)集，并從該外邊的候選點(diǎn)集中選擇一個(gè)最優(yōu)點(diǎn)，增量式地生成新的三角形。

參考文獻(xiàn)

[1]白春紅.基于領(lǐng)域本體的語義圖像檢索研究[D] .重慶大學(xué)，2012.

[2]常方媛.基于無人機(jī)航拍圖像的三維地形重建[D].天津大學(xué)，2014.

[3]陳姣，米鴻燕，李超，等.基于無人機(jī)遙感系統(tǒng)應(yīng)用與影像有效利用的探討[J].價(jià)值工程. 2013：181-183.

[4]陳蔓，鐘勇，李振東. 基于SIFT字典學(xué)習(xí)的引導(dǎo)濾波多聚焦圖像融合[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào)，2018，50（11）：65-72.

[5]陳思環(huán).基于散亂點(diǎn)云數(shù)據(jù)的表面重建關(guān)鍵技術(shù)研究[D].中國(guó)科學(xué)院大學(xué).2018.

[6] 陳庭旺，王慶.面向大規(guī)模場(chǎng)景的三維線段檢測(cè)算法[J].計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào)，2011，23（5）：790-796.

[7]郭向坤.大規(guī)模三維地形構(gòu)建的關(guān)鍵技術(shù)研究.中國(guó)科學(xué)院沈陽計(jì)算技術(shù)研究所.2019

作者簡(jiǎn)介：董興鵬（1986- ），男，漢族，本科，研究方向：計(jì)算機(jī)應(yīng)用。

陳偉（1979- ），男，碩士，研究方向：人工智能與大數(shù)據(jù)應(yīng)用。

基金項(xiàng)目：本文受高分辨率對(duì)地觀測(cè)系統(tǒng)重大專項(xiàng)（GFZX0404130302）項(xiàng)目資助。