任江峰 穆志杰

（自然資源部第二地形測量隊，陜西 西安 710054）

0.引言

近年來，自然資源部持續(xù)推進實景三維中國建設(shè)，并在北京、上海、武漢、西安等開展了一系列的項目試點和技術(shù)攻關(guān)工作，為多源數(shù)據(jù)采集和三維模型數(shù)據(jù)生產(chǎn)等方面積累了寶貴的技術(shù)經(jīng)驗。同時，各級地方政府也在全面響應(yīng)實景三維中國建設(shè)的各項政策和要求，積極擁抱新技術(shù)，探索新的應(yīng)用場景。其中，深圳基于傾斜攝影和激光雷達等技術(shù)[1]，初步建成了城市級實景三維模型，有力支撐了地方社會經(jīng)濟發(fā)展和自然資源管理等工作。2022年，山東青島作為全國典范，全面實現(xiàn)了面向大范圍復(fù)雜城市場景的優(yōu)視航測多機協(xié)同實景三維數(shù)據(jù)采集實施驗證[2]。

目前，各行各業(yè)對傳統(tǒng)基礎(chǔ)測繪地圖產(chǎn)品的需求仍然巨大，特別是大比例尺地圖。在城市級實景三維建設(shè)過程中，可以基于基礎(chǔ)地理實體數(shù)據(jù)派生出大比例尺地圖產(chǎn)品，相應(yīng)的技術(shù)路線和生產(chǎn)工藝處于研究和探索之中，基于地理實體圖元無法完全派生出大比例尺地圖，還需要使用高程點、注記和等高線等其他傳統(tǒng)測繪成果來參與大比例尺地圖測圖和生產(chǎn)過程。由于整個過程費用高、生產(chǎn)周期長、效率低，亟需基于實景三維中國建設(shè)過程中積累的全息測繪和實景三維建模等手段，開展1∶500—1∶2 000大比例尺測圖，以滿足城市規(guī)劃、建設(shè)和精細化治理等過程中不斷出現(xiàn)的需求，服務(wù)于城市管理、水利水務(wù)、生態(tài)保護、應(yīng)急指揮、交通出行、文化旅游、疫情防控等應(yīng)用場景，實現(xiàn)社會治理智能化、交通出行便捷化、生態(tài)文明可持續(xù)化，助力國家高質(zhì)量發(fā)展。

1.總體技術(shù)路線

本文采用全息測繪等手段，采集高精度傾斜實景三維模型，并基于實景模型和裸眼3D觀測等技術(shù)進行大比例尺測圖，從而生產(chǎn)1∶500—1∶2 000范圍的大比例尺DLG矢量?；趯嵕叭S模型的大比例尺測圖技術(shù)路線如圖1所示。

圖1 基于實景三維模型的大比例尺測圖技術(shù)路線圖

1.1 實景三維建模

實景三維模型構(gòu)建包含數(shù)據(jù)源采集和實景三維模型自動化生產(chǎn)兩部分內(nèi)容：

首先，相對于傳統(tǒng)RTK+全站儀數(shù)據(jù)采集方式，大比例尺測圖對數(shù)據(jù)采集效率和數(shù)據(jù)源地面分辨率提出了較高的要求。為了方便后續(xù)輸出1∶500—1∶2 000大比例尺地圖成果，此階段可以通過載人飛機或者無人機搭載傾斜航攝儀或者五鏡頭相機等設(shè)備并低空飛行，對地面進行全息測繪，從而高效地獲取地面的高精度傾斜影像。

其次，基于上述低空傾斜航飛獲取的五鏡頭影像，利用實景三維建模工具，基于集群算力和AI自動化算法，快速提取地面景物的Mesh格網(wǎng)模型數(shù)據(jù)，數(shù)據(jù)格式通常為.osgb。

1.2 大比例尺測圖

傳統(tǒng)測圖是基于帶有外方位元素的立體像對來開展的，需要采購昂貴的數(shù)字攝影測量工作站或航空立體采編平臺，如武漢大學(xué)的VirtuoZo和美國鷹圖公司的ImageStation等，對人員技能要求也較高，采集效率普遍不高；而基于實景三維模型的大比例尺測圖主要是基于傾斜或者點云等三維地理場景進行“所見即所得”方式的地圖采編，上手門檻較低，大幅度降低了項目成本，提升了大比例尺測圖的生產(chǎn)效率。

2.實景三維模型構(gòu)建

本文基于外業(yè)航空攝影來獲取影像和GPS/IMU數(shù)據(jù)，并在業(yè)內(nèi)通過攝影測量工具生產(chǎn)實景三維模型，具體內(nèi)容包括：原始航飛、數(shù)據(jù)分析和預(yù)處理、區(qū)域網(wǎng)空中三角測量、三維重建等。實景三維模型數(shù)據(jù)生產(chǎn)的主要技術(shù)路線如圖2所示。

圖2 實景三維模型構(gòu)建技術(shù)流程圖

為了進行大比例尺測圖，需要首先獲取實景三維模型。

2.1 外業(yè)數(shù)據(jù)采集

利用機載航攝儀或者無人機搭載五鏡頭相機傾斜攝影系統(tǒng)進行中低空航攝，獲取地面分辨率優(yōu)于6 cm的航空數(shù)碼像片，同時采集少量像片控制點。

2.1.1 傾斜航飛

采用五目相機，在航線規(guī)劃軟件中進行合理的航飛任務(wù)和航線規(guī)劃，并提前人工布設(shè)靜態(tài)GPS基站或者采用連續(xù)運行參考站（CORS），然后按照規(guī)劃的航線進行航攝作業(yè)。

2.1.2 像控點采集

為了保證后續(xù)傾斜空中三角測量順利開展，需要提前進行像控點采集，包括采集一定數(shù)量的平高點和高程點，采集的點均應(yīng)滿足1∶500大比例尺測圖國家規(guī)范要求，同時需充分考慮后續(xù)分區(qū)空三和測區(qū)合并工作的順利開展。在分區(qū)的外邊界和內(nèi)部邊界應(yīng)設(shè)置合理數(shù)量的像控點，每平方千米布設(shè)3～5個控制點。應(yīng)選擇地面上固定且清晰可見的目標來標記像控點，如房屋的四角、斑馬線的端點、井蓋的中心點等。

2.2 內(nèi)業(yè)模型生產(chǎn)

利用全自動三維建模軟件進行空中三角測量并自動重建三維模型，在傾斜模型基礎(chǔ)上提取實景三維模型。

2.2.1 數(shù)據(jù)分析和預(yù)處理

對提交的航攝成果進行二次檢查分析，確保數(shù)據(jù)完整、元數(shù)據(jù)信息完備。保證POS記錄信息與影像成果對應(yīng)一致，無缺漏和不匹配現(xiàn)象；保證圖像質(zhì)量和文件格式正確無誤，并且相機參數(shù)信息齊全。

2.2.2 空中三角測量

為了確定項目中每個圖像在照片或圖像曝光時的位置和方向正確，需要開展空中三角測量[3]，得到的參數(shù)稱為外方位元素。為了精確估計外方位元素，無論項目中包含多少圖像，整個區(qū)塊至少有三個地面控制點。

為了進行空中三角測量，首先需要提前收集傾斜航飛的影像數(shù)據(jù)、POS數(shù)據(jù)，以及像控點數(shù)據(jù)[4]。其次，采用改進SIFT的傾斜無人機影像匹配方法生產(chǎn)可靠的同名點[5]。一旦生成大量的地面同名點，則可以使用插值算法加密出大量的地面控制點，從而大大減少空中三角測量對外業(yè)像控點的依賴，間接降低項目成本。目前用于空中三角測量的商業(yè)軟件有Bently旗下的Context Capture、?？怂箍礗MAGINE Photogrammetry、街景工廠Street Factory等。

2.2.3 自動化建模

空中三角測量為傾斜影像數(shù)據(jù)生成了精確的外方位元素，這樣每張像片就具有了精確的位置?？梢曰诖碎_展自動化建模工作，通過人工智能和遙感圖像智能化識別技術(shù)來提取地物實景三維模型，主要過程包含點云匹配、TIN三角網(wǎng)構(gòu)建和優(yōu)化、自動紋理映射、質(zhì)檢和編輯、坐標轉(zhuǎn)換等步驟。常用的自動化建模工具：大勢智慧重建大師、Skyline PhotoMesh、大疆智圖、天際航DP-Smart全自動建模平臺、瞰景Smart3D等。

（1）點云匹配：根據(jù)空中三角測量所計算出的精確的外方位元素，對傾斜影像進行密集點云匹配，生成海量的DSM數(shù)據(jù)。

（2）TIN三角網(wǎng)構(gòu)建和優(yōu)化：基于DSM中的密集點，采用內(nèi)插算法，全自動構(gòu)建TIN三角網(wǎng)，并對三角網(wǎng)進行平滑等處理。

（3）自動紋理映射：經(jīng)過點云匹配和三角網(wǎng)構(gòu)建，已經(jīng)生成了實景三維白模，需要利用紋理匹配算法從海量圖像中選取質(zhì)量最優(yōu)的影像進行紋理映射處理，從而生成具備紋理色彩的實景三維模型。

（4）質(zhì)檢和編輯：由于高樓和樹木之間存在遮擋、實景三維模型在水域存在破洞等問題，需要進行人工質(zhì)檢，標記出問題區(qū)域并手工編輯、修復(fù)。

（5）坐標轉(zhuǎn)換：為了滿足后續(xù)大比例尺測圖對坐標系的個性化需求，需要對傾斜實景三維模型進行坐標轉(zhuǎn)換，如，從WGS84坐標系轉(zhuǎn)換到2 000地方坐標系。

3.大比例尺測圖

在構(gòu)建完實景三維模型后，在實景三維智能測圖工具中，導(dǎo)入高精度的實景三維模型數(shù)據(jù)，便可正式開展大比例尺測圖。與傳統(tǒng)的基于立體像對的測圖模式有所不同，由于此處加載的是實景三維模型，所以無需佩戴3D眼鏡，而直接采用裸眼3D方式進行觀測，大大提高了測圖的便捷性，節(jié)約了成本。大比例尺測圖的主要流程如圖3所示。

圖3 大比例尺測圖技術(shù)流程圖

基于裸眼3D方式，可以采集1∶500—1∶2 000大比例尺下的各種地物矢量要素，包括房屋、道路、水系、植被、行政區(qū)界線、高程點、等高線等。在人工作業(yè)的過程中，可以結(jié)合包含直角繪制、直線繪制、智能繪制等在內(nèi)的智能化手段來輔助3D采集，從而大大提高生產(chǎn)效率。

4.實驗流程與精度評價

為了論證本文提出的實景三維模型在大比例尺測圖工作中的可行性，精心設(shè)計以下實驗并通過實驗數(shù)據(jù)對精度進行檢驗。

4.1 測區(qū)概況

測區(qū)位于陜西省西安市某區(qū)域，測區(qū)總面積約為0.2 km2。測區(qū)主要以農(nóng)房為主，房屋和樹木比較密集，對GNSS衛(wèi)星信號具有一定影響且通視條件一般。

4.2 實驗設(shè)備

本項目采用國產(chǎn)軟硬件，主要包括：

（1）智航SF700A四旋翼無人機：最大續(xù)航約80 min，起飛重量為5.7 kg，支持仿地飛行／斷點續(xù)飛。

（2）五鏡頭T53P傾斜相機：傳感器尺寸為23.5×15.6 mm，下視焦距為25 mm，側(cè)視焦距為35 mm，圖像分辨率為6 000×4 000。

（3）南方銀河6 RTK測量系統(tǒng)：獲取檢查點的三維坐標。

（4）天云三維一體機：用于運行傾斜空三和實景三維重建。

（5）Ufeature3D：三維立體測圖平臺。

（6）實景三維GIS平臺：用于檢驗大比例尺DLG效果和精度。

4.3 實驗過程

根據(jù)1∶500大比例尺測圖要求，本次測區(qū)高差約20 m，滿足無人機低空攝影測量技術(shù)規(guī)范。因此決定采取無人機低空傾斜攝影測量的技術(shù)手段獲取測區(qū)的實景三維模型，并基于實景模型采集1∶500的DLG數(shù)字線劃圖，詳細步驟如下：

（1）采用智航SF700A四旋翼無人機+T53P傾斜相機獲取測區(qū)的傾斜影像。相機最優(yōu)拍照間隔為0.9 s，保證航測航向重疊度80%，旁向重疊度75%，設(shè)置飛行速度為9m/s，飛行高度300 m，飛行時間46 min，共拍照9280張，場地內(nèi)共布設(shè)像控點8個，布設(shè)房角檢查點12個。

（2）基于天云三維一體機對測區(qū)進行區(qū)域網(wǎng)空中三角測量，導(dǎo)出空三加密成果，精度較好，具體為：平面中誤差為0.015 m，高程中誤差為0.021 m。

（3）在天云三維一體機中導(dǎo)入空三加密成果，進行自動化三維重建，構(gòu)建實景三維模型。

（4）將實景三維模型導(dǎo)入Ufeature3D軟件中，進行智能化三維立體測圖，并采集房屋要素，以滿足房地一體的應(yīng)用需求。

4.4 結(jié)果分析

將1∶500房屋DLG數(shù)據(jù)導(dǎo)入實景三維GIS平臺：

（1）效果評估：房屋DLG數(shù)據(jù)與正射影像底圖套合完美，效果較好如圖4所示。

圖4 大比例尺測圖效果圖

（2）精度評估

利用南方銀河6 RTK測量系統(tǒng)隨機采集DLG矢量對應(yīng)的屋頂檢查點20個，并將檢查點坐標導(dǎo)入實景三維GIS平臺，然后軟件自動計算檢查點與對應(yīng)DLG相同位置坐標的限差：平面最大限差為0.048 m、0.032 m；平均平面中誤差為0.037 m，平均高程中誤差為0.025 m，優(yōu)于1∶500地形圖（1∶500地形圖平面中誤差應(yīng)不超過0.25 m，高程中誤差應(yīng)不大于0.17 m）。

5.結(jié)束語

基于以上研究得出以下結(jié)論：

以全息測繪等智能化方式快速構(gòu)建的實景三維模型為基礎(chǔ)，進行1∶500—1∶2 000大比例尺裸眼3D測圖的技術(shù)路線完全可行。該技術(shù)手段可以充分利用各地實景三維項目中積累的存量實景三維模型數(shù)據(jù)，有效地降低了項目成本，成圖效果較好，同時也滿足大比例尺地形圖成圖精度要求，具有較高的應(yīng)用價值。