朱 磊,趙 菲,王宏昌

(1. 山東省國(guó)土測(cè)繪院,山東 濟(jì)南 250102; 2. 山東省國(guó)土空間生態(tài)修復(fù)中心,山東 濟(jì)南 250014)

地理場(chǎng)景是指一定時(shí)空范圍內(nèi)的各種自然要素和人文要素相互聯(lián)系、相互作用所構(gòu)成的具有特定結(jié)構(gòu)和功能的地域綜合體。地形級(jí)地理場(chǎng)景主要反映地形特征及地表地物信息,側(cè)重于對(duì)地表高程、紋理等大范圍地理環(huán)境的三維描述,適用于現(xiàn)實(shí)世界的宏觀可視化表達(dá)。地理場(chǎng)景數(shù)據(jù)是對(duì)地理場(chǎng)景的空間位置、形態(tài)、紋理進(jìn)行數(shù)字化表達(dá)的地理信息數(shù)據(jù),是新型基礎(chǔ)測(cè)繪的一種基礎(chǔ)產(chǎn)品,是實(shí)景三維中國(guó)建設(shè)的重要組成部分[1]。

地理場(chǎng)景構(gòu)建與快速更新技術(shù)已成為測(cè)繪界專家學(xué)者研究的熱點(diǎn),并取得了豐碩的成果。文獻(xiàn)[2]以高速固定翼飛機(jī)搭載多鏡頭傾斜航攝儀獲取3 cm分辨率傾斜影像,經(jīng)空三加密、稠密點(diǎn)云生成、白模構(gòu)建、紋理映射等步驟實(shí)現(xiàn)了城市級(jí)地理場(chǎng)景快速構(gòu)建技術(shù);文獻(xiàn)[3]針對(duì)地理場(chǎng)景數(shù)據(jù)瓦片式建模的特點(diǎn),綜合模型數(shù)據(jù)生產(chǎn)組織和時(shí)空索引更新處理兩個(gè)階段,提出了適用于同類數(shù)據(jù)增量更新和不同類數(shù)據(jù)聯(lián)動(dòng)更新的技術(shù)策略和方法,研究了基于傾斜攝影模型的地理場(chǎng)景更新方法;文獻(xiàn)[4]為解決地理場(chǎng)景更新的及時(shí)性和高效性,提出了利用無(wú)人機(jī)獲取變化圖斑點(diǎn)云數(shù)據(jù),經(jīng)濾波提取地面點(diǎn),并進(jìn)行平滑處理,輸出數(shù)字模型,實(shí)現(xiàn)了地形級(jí)實(shí)景三維地理場(chǎng)景更新。

地形級(jí)地理場(chǎng)景構(gòu)建一般采用DEM疊加DOM的方式,DEM與DOM生產(chǎn)一般需獲取兩次數(shù)據(jù),且影像與點(diǎn)云時(shí)相較難一致。本文提出一種基于Leica CityMapper-2L混合航攝儀同步獲取的點(diǎn)云與影像數(shù)據(jù),聯(lián)合生產(chǎn)DEM和DOM融合制作地形級(jí)地理場(chǎng)景的技術(shù)方法,以及對(duì)低旁向重疊度航空影像連接點(diǎn)匹配等關(guān)鍵技術(shù)進(jìn)行分析,并通過試驗(yàn)對(duì)該方法進(jìn)行驗(yàn)證。

1 關(guān)鍵技術(shù)

1.1 技術(shù)流程

基于混合數(shù)據(jù)的地形級(jí)地理場(chǎng)景構(gòu)建主要技術(shù)流程為:①采用CityMapper-2L混合航攝儀同步獲取點(diǎn)密度優(yōu)于1個(gè)/m2的點(diǎn)云數(shù)據(jù)和0.2 m分辨率的航空影像數(shù)據(jù);②對(duì)點(diǎn)云數(shù)據(jù)進(jìn)行航帶匹配、濾噪、自動(dòng)分類和細(xì)分類處理,生成DEM和DSM數(shù)據(jù);③對(duì)航空影像進(jìn)行空三加密、基于密集匹配點(diǎn)云和同步LiDAR點(diǎn)云對(duì)航片進(jìn)行數(shù)字微分糾正、鑲嵌處理制作DOM;④DEM和DSM疊加DOM數(shù)據(jù),融合制作地形級(jí)地理場(chǎng)景數(shù)據(jù),制作流程如圖1所示。

圖1 地形級(jí)地理場(chǎng)景制作流程

1.2 DEM制作

通過IMU/DGPS數(shù)據(jù)解算、點(diǎn)云預(yù)處理、自動(dòng)分類、人工編輯分類等工序,提取符合要求的地面點(diǎn)云,然后構(gòu)建TIN、重采樣制作DEM。

1.3 低重疊度影像連接點(diǎn)匹配

航空影像數(shù)據(jù)處理過程中,連接點(diǎn)匹配受影像質(zhì)量、鏡頭畸變、地形、POS數(shù)據(jù)等多種因素的影響,其匹配效果直接影響空三精度[5-6]。

當(dāng)相鄰航帶的影像重疊度低于20%時(shí),影像之間的共同特征點(diǎn)數(shù)量減少,對(duì)連接點(diǎn)匹配造成一定困難。此時(shí),連接點(diǎn)匹配可采取根據(jù)特征點(diǎn)方向分量和灰度進(jìn)行匹配的方式完成,并通過最小二乘法(LSM)提高匹配精度。通過對(duì)圖像進(jìn)行多層不同程度的模糊處理,在每個(gè)金字塔層上運(yùn)行自動(dòng)的整體測(cè)區(qū)平差程序,以達(dá)到對(duì)連接點(diǎn)的強(qiáng)有力的內(nèi)部質(zhì)量控制,即使在結(jié)構(gòu)較差的地區(qū)及山區(qū)也能產(chǎn)生有效的連接點(diǎn)。當(dāng)連接點(diǎn)提取沒有達(dá)到預(yù)期的精度和連接性時(shí),再通過最小二乘法自動(dòng)優(yōu)化連接點(diǎn),進(jìn)一步優(yōu)化結(jié)果。通過組合匹配方式,可使連接點(diǎn)的匹配達(dá)0.1像元的精度。

1.4 影像數(shù)字微分糾正

航空影像數(shù)字微分糾正可基于影像密集匹配生成的點(diǎn)云或同步獲取的LiDAR點(diǎn)云兩種方式進(jìn)行[7-9]。基于密集匹配點(diǎn)云進(jìn)行數(shù)字微分糾正,優(yōu)點(diǎn)是可以利用現(xiàn)有的航空影像資源,但生成的點(diǎn)云精度可能受到影像質(zhì)量、特征點(diǎn)匹配算法等因素的影響;基于LiDAR點(diǎn)云進(jìn)行數(shù)字微分糾正,優(yōu)點(diǎn)是LiDAR點(diǎn)云具有較高的精度和穩(wěn)定性,可以提供更準(zhǔn)確的糾正結(jié)果。

LiDAR點(diǎn)云濾波方式生成的 DEM 在平坦區(qū)域微分糾正效果較好,房子基本不會(huì)變形,山區(qū)植被覆蓋率高的地方正射糾正的單片則會(huì)產(chǎn)生大量的樹木拉伸現(xiàn)象,而采用影像密集匹配點(diǎn)云進(jìn)行影像正射糾正效果較好,兩者之間形成互補(bǔ)。因此,采用基于LiDAR點(diǎn)云為主、密集匹配點(diǎn)云為輔的微分糾正方式,可取得理想的糾正效果[10-11]。

1.5 融合構(gòu)建地形級(jí)地理場(chǎng)景

以DEM為基底,疊加DOM作為地面紋理,并進(jìn)行場(chǎng)景編譯[12-13]。此種方式的優(yōu)點(diǎn)是山體、地形等立體感表現(xiàn)較好,高程精度高,可進(jìn)行精確的土方量算、淹沒分析等。缺點(diǎn)是房屋、橋梁等地物因DEM置平,表現(xiàn)不出三維立體形態(tài),陡涯等地形高差大的地區(qū)紋理不足,有拉花現(xiàn)象。

2 試驗(yàn)與分析

2.1 CityMapper-2混合航攝儀

CityMapper-2混合航攝儀是徠卡測(cè)量系統(tǒng)混合型傳感器(如圖2所示),集成了五視角傾斜相機(jī)和高性能激光雷達(dá),單次飛行即可獲取高分辨率五視角影像和高密度點(diǎn)云成果,影像與激光點(diǎn)云多源數(shù)據(jù)提供的信息具有一致性、冗余性和互補(bǔ)性,可以進(jìn)行目前航空攝影測(cè)量常見的4D產(chǎn)品和三維模型生產(chǎn)。

圖2 CityMapper混合航攝儀組成

2.2 試驗(yàn)情況

地理場(chǎng)景產(chǎn)品分級(jí)建設(shè),省級(jí)層面負(fù)責(zé)亞米級(jí)地形場(chǎng)景建設(shè),主要由0.5～2 m格網(wǎng)尺寸的DEM和0.2～1 m分辨率的DOM數(shù)據(jù)融合形成。在《實(shí)景三維山東建設(shè)總體實(shí)施方案》(2023-2025年)中,山東省級(jí)建設(shè)0.2 m分辨率的DOM和2 m格網(wǎng)的DEM[14]。

以山東某區(qū)域地形級(jí)地理場(chǎng)景建設(shè)項(xiàng)目為例,試驗(yàn)區(qū)位于膠東丘陵地區(qū),采用CityMapper-2L混合航攝儀同步獲取點(diǎn)云和影像數(shù)據(jù),相對(duì)航高為3777 m,飛行速度為300 km/h。激光視場(chǎng)角設(shè)置為40°,旁向重疊度為20%,相機(jī)下視角旁向視場(chǎng)角為41.2°,影像旁向重疊度為22%,單架次獲取優(yōu)于1點(diǎn)/m2點(diǎn)云數(shù)據(jù)和GSD為0.2 m的航空影像數(shù)據(jù)約1800 km2。

首先對(duì)LiDAR點(diǎn)云進(jìn)行航帶匹配、濾噪、濾波處理,制作1 m格網(wǎng)DEM;其次航空影像匹配采用根據(jù)特征點(diǎn)方向分量和灰度進(jìn)行匹配且與最小二乘法組合的方式,這樣克服了航帶重疊度小的難題;然后經(jīng)空三加密,基于LiDAR點(diǎn)云與密集匹配點(diǎn)云相結(jié)合的方式進(jìn)行微分糾正,以及影像鑲嵌等步驟制作DOM;最后DEM疊加DOM并進(jìn)行場(chǎng)景編譯制作省級(jí)地形級(jí)地理場(chǎng)景,如圖3所示。

圖3 地形級(jí)地理場(chǎng)景

2.3 精度分析

(1)平面精度。在試驗(yàn)區(qū)域內(nèi)均勻選擇60個(gè)檢測(cè)樣本區(qū)域,并與采用網(wǎng)絡(luò)RTK方法采集的平面坐標(biāo)進(jìn)行對(duì)比,計(jì)算模型平面位置誤差。采用高精度統(tǒng)計(jì)方法,得到模型平面位置誤差為0.498 m。

(2)高程精度。在試驗(yàn)區(qū)均勻選擇6組精度驗(yàn)證區(qū),采用網(wǎng)絡(luò)RTK方法采集了312個(gè)高程檢核點(diǎn),去除粗差后,共獲得305個(gè)高程差值,誤差值呈正態(tài)分布,中誤差為0.127 m。其中,21%的誤差值優(yōu)于4 cm;36%的誤差值優(yōu)于6 cm;74%的誤差值優(yōu)于8 cm;95%的誤差值優(yōu)于10 cm。

3 結(jié) 語(yǔ)

本文對(duì)地形級(jí)地理場(chǎng)景建設(shè)技術(shù)現(xiàn)狀進(jìn)行了分析,提出了基于同源混合遙感數(shù)據(jù)制作地形級(jí)地理場(chǎng)景的方法和技術(shù)路線,并對(duì)該技術(shù)中存在的技術(shù)難題給出了解決方案。結(jié)合省級(jí)基礎(chǔ)測(cè)繪項(xiàng)目,通過山東某地區(qū)對(duì)該技術(shù)方法進(jìn)行了驗(yàn)證。

采用該技術(shù)進(jìn)行地形級(jí)實(shí)景三維更新,具有高程精度高、紋理與模型時(shí)相一致、節(jié)省更新成本等特點(diǎn)。該技術(shù)革新了地形級(jí)地理場(chǎng)景構(gòu)建模式,在地形級(jí)地理場(chǎng)景更新,尤其是帶狀區(qū)域(新建高速鐵路、公路、運(yùn)河等)更新具有明顯優(yōu)勢(shì),且具有極高的推廣應(yīng)用價(jià)值。