基于航天遠景MapMatrix的DEM制作研究

張新雯

摘 要：MapMatrix是基于衛星遙感、航空、外業等數據進行多源空間信息綜合處理的平臺。該文總結了ADS80數碼影像在MapMatrix下進行某市1∶2 000DEM數據生產流程，探討DEM快速生成的方法，成果通過精度評定與質量檢查，確實可靠，具有良好可行的運用前景。

關鍵詞：MapMatrix DEM 航天遠景 ADS80

中圖分類號：TD672 文獻標識碼：A 文章編號：1672-3791（2016）11（b）-0001-03

某市1∶2 000DEM以東西向鋪設條帶，采用ADS80數碼影像數據，經IPAS軟件對空中IMU、GPS數據進行數據預處理，解算出飛機（相機）的位置和姿態數據。結合解得的定位定姿參數和GPRO軟件，對原始L0級影像作定向定位糾正，得到WGS84坐標系統下的L1級影像。利用MapMatrix系統完成本地坐標系統轉換及恢復立體模型采集等工序，使用量測得出的特征點、線構TIN及影像匹配多種方式，生成區域DEM。

1 基于ADS80影像下的數據采集

1.1 ADS80工作原理

ADS80是徠卡公司于2008年在ADS40基礎上開發的多線陣CCD成像技術機載線陣攝影系統，其集成的POS系統（GPS/IMU）可在少量地面控制點甚至無控制點情況下進行空三解算。相較傳統膠片相機和框幅式數碼相機而言，ADS80可直接得到影像的外方位元素，大大縮短外業相控、空三加密解算等人工干預時間；并同時獲取同一地面的前視（27°）、中視（2°）及后視（14°）3個高分辨率無縫連續重疊的地面立體影像條帶，大大減少了相片掃描和拼接環節。根據空間前方交會得到的高程精度公式，提高高程精度的因素有3點：提高基高比b/f，提高像點坐標的量測精度（k ），提高GSD地面分辨率。

ADS80通過更小的傳感像元，降低航飛高度，形成0.87的基高比（B/H=tan41°）獲得較高的高程精度，彌補了以往框幅式像對的立體采集中高程精度不足的缺陷，其整體性能使數字攝影測量立體采集工作更為簡捷精確。

1.2 基于MapMatrix的ADS80立體測圖原理

MapMatrix是基于衛星遙感、航空、外業等數據進行多源空間信息綜合處理的平臺。其涵蓋幾個重要的模塊：ATMatrix空三加密模塊、DEMMatrix高程模型處理模塊、DOMMatrix正射影像處理模塊、FeatureMatrix立體測編模塊。利用MapMatrix進行ADS80影像的建模過程相較框幅式相機而言相對簡化，無需進行內定向、相對定向、絕對定向及核線重采樣等工序。用戶在該系統中可由原始影像出發，由相機參數文件（*.cam）、影像頭文件（*.sup）、金字塔影像（*.tif）、姿態定向參數（*.odf）生成立體模型，經過特征點、線采編及基于影像的自動匹配，最終生成處理后的DEM產品。

2 DEM的制作方法及技術

2.1 DEM多種制作方法及技術

DEM的數據采集方法主要包括3種：地形測量獲取、航空攝影測量、利用LiDAR點云分類提取。地形測量獲取方法如用GPS、全站儀野外測量采集等，該方法外業強度大、效率低、成本高，且對于生產產品級別的DEM數據而言，明顯高程點密度不足，精度質量欠佳。航空攝影測量方法即全數字攝影測量系統獲取DEM，能利用核線影像自動匹配生成DEM后再人工編輯，又能在立體環境下通過采集特征點線面后構TIN生成DEM。該方法效率較高，成本投入少，且能夠滿足產品級DEM的生產精度要求。利用LiDAR點云進行濾波分類提取DEM，其數據精度高、生產效率最快，但需要額外配置專業軟件對海量點云數據進行濾波處理。由于三維點云數據在三維建模等領域應用較廣，如僅從獲取DEM產品角度而言，航攝成本相對較高。

利用全數字攝影測量方式進行DEM制作又可以分為基于TIN直接內插DEM和物方相關原理生成DEM兩種。TIN的創建，主要依據于基礎數據的采集，即由矢量數據轉換成的特征點、線，按照物方DEM間隔規定內插規則正方形格網DEM，得到DEM數據。利用物方相關原理生成DEM需進行影像相關處理，進行影像匹配后MapMatrix系統自動生成DTM，但對于大比例尺（特別平坦地區），由于匹配點大量在樹木和建筑物上，由此自動生成的像方DEM效果很差，往往需要進行大量人工干預，在生成像方DEM后逐塊編輯，最終生成基于物方的DEM。

2.2 基于MapMatrix的DEM具體制作方法

數據采集是DEM生成的關鍵問題。數據點太稀會降低DEM的精度；數據點過密，又會增大數據量、處理的工作量和存儲量。因而在DEM數據采集之前，需要按照成果的精度要求確定合理的取樣密度，并在DEM數據采集過程中根據地形復雜程度動態調整采樣點密度。

該市地勢平坦，現代化程度高，區域房屋密集，高架及橋梁交通設施發達，如采用影像匹配自動生成DTM，由于房屋、樹木，高架眾多，自動匹配結果將混亂，因此在此次DEM生產中采用特征點、線，直接創建物方DEM，可大大減少編輯工作量。對于地貌交待清晰、易于表現的區域，可采用線編輯模式，選擇缺省線屬性對能表達地形變化區域進行立體采集。特征線具有單一立體效果，不遮蓋影像，能對地形細膩表示，檢查出各種小粗差。對于地貌凌亂、無法清晰表達的區域，可采用點編輯模式。根據CJJ100-2004的規定，城市DEM數據的基本格網尺寸應為5 m×5 m，利用MapMatrix中FeatureOne立體采集模塊設置自動采集高程點選項，以水平方向為行，順序從上至下，以垂直方向為列，順序從左至右，每隔5 m自動跳轉至下個格網，量測高程點。該方法操作簡單易行，但工作量較大，一般作為雜亂區域描述較好，步距及斷面線上點的間隔都可根據地形的復雜程度而改變。具體要素采集遵循如下原則。

（1）水域：對于靜止水面，測量水位高程并按此高程采集水岸線，整個水域范圍據此高程構建平三角形，并按此高程對DEM格網賦值。雙線河流水岸線的高程應依據上下游水位進行分段內插賦值，DEM高程值應自上而下平緩過渡，并且與周圍地形高程之間的關系正確、合理。（2）森林區域：在林區，DEM測量的是樹頂表面，在生成DEM格網時應減去平均樹高獲取地面高程。（3）特殊區域：山頭、凹地或埡口等處應內插高程特征點，狹長而緩坡的溝谷或山脊應內插特征線，避免出現不合理的平三角形；陡巖、斜坡、雙線沖溝等地貌應合理反映地形特征。（4）空白區域：空白區域是指數據源出現局部中斷等原因無法獲取高程的區域，位于空白區域的格網高程值賦予-9999。

2.3 基于MapMatrix的DEM編輯

DEM點應切準地面，等高線真實地反映地貌形態。為提高效率，在MapMatrix中一般采用面編輯的方式，使用系統提供的平滑、內插、擬合、定值及平均高程賦值等算法，逐塊編輯；對于道路等線狀要素也可以采用線編輯方式進行操作。

2.4 DEM接邊

像對間DEM設置10個格網間距的重疊區域用于接邊。系統自動對重疊區域內的格網點進行高程較差的統計分析，在2～3倍高程接邊誤差的點位應控制在4%以內，不得出現3倍以上中誤差的點。一旦發現需進行修測，符合限差要求后應進行像對DEM接邊，取平均數作為重疊區域內的數據值。

3 DEM的質量控制

DEM質量的最終檢查是根據野外高程點的平面坐標。在已建立的DEM中內插出檢測點位置的高程，利用航攝像控點庫以及外業RTK采集的高程檢查散點；在ArcGIS軟件中，檢查內插等高程點與原始高程點之間的偏離量是否在規定范圍內；相鄰存儲單元的DEM數據應平滑銜接；對于水域需檢查靜止水域內的DEM格網點高程是否保持一致，流動水域上下游DEM格網點高程是否呈梯度下降；并對高程較差進行統計分析。對較差較大部分需返回到立體模型上進行上機檢查，實時觀察生成的DEM點位是否切準地面。如果DEM與地面模型的高程差在2倍中誤差以上，則需進行重測。

DEM成果精度用格網點的高程中誤差Mz表示：

式中，v為高程較差；n為檢測點個數。根據城市基礎地理信息系統技術規范要求（CJJ100-2004），城市DEM數據的基本格網尺寸為5 m×5 m，采用5 km×5 km分幅存儲，精度等級為二級。據此規定，平地、丘陵、山地、高山地格網點高程中誤差分別為：0.7 m、1.7 m、3.3 m、6.7 m。高大林木覆蓋區、高層建筑陰影遮蓋區等困難區域的平面和高程中誤差可放寬50%。DEM高程值應取位至0.1 m。

DEM的質量檢查還需進行數據文件及數據完備性檢查。數據文件檢查包括DEM數據文件命名、數據格式、數據分幅、數據格網尺寸是否符合要求；數據完備性檢查包括檢查DEM數據覆蓋范圍有無不滿幅、數據有無遺漏等問題，相鄰存儲單元之間數據完整，不得出現漏洞，DEM數據應覆蓋整個區域范圍，接邊范圍數據應有一定的重疊。

4 實驗數據

該項目利用其他項目單位同源DLG數據外業高程檢查散點，對DEM成果高程坐標值進行檢核。該次核查初次抽查10幅圖，共核查點2 269個，符合精度的點為98.5%。利用公式求得高程坐標中誤差，經測算DEM成果精度優于5 m格網城市基礎地理信息系統規范要求的平地0.7 m，高大林木覆蓋區、高層建筑陰影遮蓋區等困難地區最大不超過1 m，丘陵1.7 m的中誤差要求。DEM較差及各圖幅中誤差分布情況如圖2所示。

5 結語

該文總結了ADS80數碼影像在MapMatrix下進行某市 1∶2 000DEM數據生產流程，探討DEM快速生成的方法，成果通過精度評定與質量檢查，確實可靠，具有良好可行的運用前景。

參考文獻

[1]劉凱，湯國安，江嶺，等.格網DEM侵蝕學坡長并行計算方法[J].武漢大學學報：信息科學版，2015（2）：274-279.

[2]熊漢江，李秀娟.一種提取山脊線和山谷線的新方法[J].武漢大學學報：信息科學版，2015，40（4）：498-502.

[3]陳潔，鄭雄偉，孫秉釗.航空攝影DG 法若干問題探討和實踐研究[J].地理空間信息，2015（2）：83-85.

[4]朱傳勇，文琳，沈瑩，等.基于MapMatrix的武漢市1∶2000DEM

制作[J].地理空間信息，2015（2）：101-103.

[5]明鋒，楊元喜，曾安敏，等.顧及有色噪聲的GPS位置時間序列中斷探測法[J].武漢大學學報：信息科學版，2016，41（6）：745-751.

[6]茍嬌嬌，羅明良，王飛.影響黃土高原集水面積閾值的地形因子主成分分析[J].武漢大學學報：信息科學版，2014（7）：83-90.