不同分辨率DEM對影像正射糾正精度的影響

董玉剛

（遼寧省自然資源事務服務中心—遼寧省基礎測繪院，遼寧 錦州 121003）

0.引言

正射影像圖是按照地理坐標生成的具有幾何精度和影像特征的數字地圖，具備直觀準確地反應地形地貌、制作周期短等特點，越來越廣泛地應用在空間地理信息建設中。遙感衛星影像在成像的過程中，由于地球曲率及地形起伏、攝影軸傾斜、大氣折光等諸多因素的影響，致使影像中每個像點產生不同程度的幾何變形而失真[1]。正射糾正不僅可以減少或消除成像過程中由于各種因素導致的影像畸變，也是各種影像應用的前提與基礎，特別是對于山區等地形起伏較大，常規的幾何校正難以消除幾何變形的區域。

研究發現：數字高程模型的分辨率決定其對高分辨率影像正射糾正的精度。因此，研究不同分辨率數字高程模型對于高分辨率衛星影像正射糾正精度的影響具有實際價值[2]。重慶市地理信息中心的陳靜，袁超等人針對實際生產中1∶2 000數字高程模型對于山區高分辨率衛星影像正射校正扭曲嚴重的現象做過數字高程模型對山區高分辨率影像正射校正精度的影響[3]。劉云峰研究了數字高程模型對遙感衛星影像正射糾正精度的影響，選取了包含山地、平地、丘陵等常見地貌的SPOT圖像，并使用所在區域1∶50 000地形圖作為試驗樣區，通多對比分析以后發現25米以及50米的格網間距對正射糾正精度的影響可以忽略[4]。本文使用資源三號衛星影像作為數據源，對衛星遙感影像進行正射糾正，并通過試驗對比分析不同格網間距的高程模型對正射糾正的精度影響的情況，以此來確定最佳分辨率的數字高程模型。

1.試驗數據準備

試驗數據主正文開始要包括資源三號影像數據和四種不同分辨率的數字高程模型，四種分辨率的數字高程模型分別是15米分辨率的DSM、30米分辨率的ASTER GDEM、45米分辨率的SRTM、90米分辨率的SRTM。然后將下載的90米分辨率數字高程模型進行內插處理，得到45米分辨率的數字高程模型。

1.1 衛星影像

資源三號衛星是我國第一顆民用高分辨率光學傳輸型立體測繪衛星，它主要用于1∶2.5 000等更大比例尺地形圖部分要素的更新，還廣泛應用于1∶50 000比例尺立體測繪和數字影像制作中。資源三號衛星配備有多光譜相機和三線陣測繪相機兩種相機，采用三線陣成像方式，通過正視相機、前視相機、后視相機三種高幾何分辨率相機和多光譜相機獲取指定的某地區的立體影像、多光譜影像、高分辨率影像以及其他輔助數據[5]。本實驗采用的標準衛星影像為資源三號全色波段影像數據，經過了傳感器校正的幾何處理方式和相對輻射校正處理，影像平面精度優于10米，高程精度優于5米。

1.2 DEM數據

數字高程模型（Digital Elevation Model）是利用地面坐標及高程數據來模擬地球表面的數據集，描述一定空間范圍的地貌特征，是基礎地理信息數字產品，能夠用來糾正遙感影像在成像過程中的投影差。根據產品的不同形式通常分為以下幾種：

（1）SRTM數據

SRTM是航天飛機雷達地形測繪的產品，使用雷達干涉儀，生成小于60°緯度的全局一致的數字高程圖，獲取的雷達影像數據經過處理，制成了分辨率為30米和90米兩種數字地形高程模型產品。

（2）ASTER GDEM數據

ASTER GDEM（先進星載熱發射和反射輻射儀全球數字高程模型）是對ASTER衛星的測繪數據進行處理生成的，其垂直精度能夠達到20米，水平精度達30米，產品覆蓋全球陸地，在大范圍地形圖制作等方面，它與SRTM兩種數據的應用較為廣泛。

（3）DSM數據

DSM（數字表面模型）是比DEM更加精細的表示地形地貌的數據集，它不僅包含了自然地球表面的形狀還包含了人工建筑物產生的高度變化，是反映地表信息的重要手段之一[6]。利用資源三號立體像對建立的15米分辨率DSM已經基本覆蓋全國。

空間地理數據云平臺是由中國科學院計算機網絡信息中心的科學數據中心研發，數字高程模型數據可以從空間地理數據云免費獲取。

2.ZY-3遙感影像正射糾正試驗

2.1 技術方法

以ZY-3遙感影像的全色波段為影像數據源、準備試驗的數字高程模型以及數字表面模型等數據，分別利用15米分辨率數字表面模型、30米分辨率數字高程模型、45米分辨率數字高程模型以及90米分辨率數字高程模型對資源三號影像進行正射糾正，將正射糾正影像與高分辨率影像進行精度對比，評價正射糾正的效果。

2.2 數字高程模型預處理

2.2.1 DEM拼接

圖像拼接是將幾幅具有一定重疊區域并且帶有地理參考的影像組合成為一幅圖像的過程。圖像輸入之前必須經過幾何校正處理或者進行校正標定[7]。由于本試驗選取的資源三號遙感影像的覆蓋范圍大于一幅數字高程模型的范圍，所以要將多幅數字高程模型的柵格圖像進行拼接，確保覆蓋試驗影像。

首先，選取覆蓋試驗區域的全部DEM數據，盡可能選擇成像時間和成像條件比較接近的數據，以減輕后面的色調調整的工作量。然后，選定拼接方案，這是極為重要的步驟，拼接方案的好壞直接關系到工作時間和工作量。本試驗確定的拼接方案為：確定標準像幅，標準像幅一般位于試驗區中央，再確定拼接順序，即以標準像幅為中心，由中央向四周逐步進行拼接，接下來，無論是幾何拼接的過程還是色調調整的過程都是以標準像幅為中心，向四周逐步擴散，完成相鄰圖像的拼接。

2.2.2 DEM裁剪

拼接后的DEM能夠完全覆蓋試驗區，需要按試驗區范圍對其進行裁剪。本次試驗采用資源三號影像數據自帶的邊界文件向外做一千米緩沖區，利用新生成的矢量文件，根據不規則分幅裁剪的方式對DEM數據進行裁剪。

2.2.3 DEM內插

DEM內插是根據若干個相鄰參考點的高程值，利用一定算法求出待定點上的高程值，在數學上屬于插值問題，任意一種內插算法都是基于原始地形起伏變化的連續光滑性，或者說比較臨近的數據點之間有很大的相關性[8]。

由于本次試驗下載的DEM數據是柵格數據，所以使用影像內插方法來進行數字高程模型內插。影像內插一般有最鄰近像元法、雙線性內插、雙三次卷積法。本次試驗采用雙三次卷積方法進行處理，這種方法的內插效果最好，但是處理時間相對較長。

2.3 正射糾正處理

許多遙感圖像處理軟件都能實現遙感影像的正射糾正處理，如ERDAS、ArcGIS等。利用成熟的遙感圖像處理軟件，在新建的處理項目選取與傳感器對應的幾何模型進行自動正射糾正即可。糾正時，分別選取15米分辨率的DSM、30米分辨率的ASTER GDEM、45米分辨率的SRTM和90米分辨率的SRTM作為影像糾正的DEM數據。

2.4 精度評定

精度評定是評價一幅圖像正射糾正結果的唯一標準。本次試驗采用的精度評價方法是尋找特征點，并判斷影像上的特征點對比標準圖像上同一個特征點在圖上像素差值的中誤差。試驗中特征點提取采用SIFT算法，SIFT算法是Scale Invariant Feature Transform算法（即尺度不變特征變換匹配算法）的簡稱，它是一種用于計算機視覺識別的算法，用來偵測和描述一景影像的局部特征，在空間尺度中尋找極值點，并提取出其位置、尺度、旋轉不變量以確定兩景影像之間的相識度。本文在利用不同DEM數據源糾正后的正射影像上分別選取了城鎮、山地、丘陵、農田以及河流等多種不同地貌類型，再進行特征點提取，統計其點位中誤差。各類特征點位中誤差如表1-表4所示。

表1城鎮特征點中誤差 單位：米

從表1可以看出：15米分辨率的DSM對影像中城鎮地區的正射糾正效果最好，而其他三種分辨率的DEM對影像正射糾正的效果雖有降低，但是降低幅度比較小。

表2山地、丘陵特征點中誤差 單位：米

從表2可以看出：同樣是15米分辨率的DSM對山地、丘陵類地形的正射糾正效果最好，而且正射糾正效果隨著DEM分辨率的降低而下降，但是與表1進行比較可以看出：山地地形整體糾正效果沒有城鎮地區的糾正效果好。

表3農田特征點中誤差 單位：米

從表3可以看出：在農田種植區15米分辨率的DSM糾正效果明顯比其他幾種分辨率的DEM的糾正效果要好，而且與城鎮區域的正射糾正中誤差相差不多。

表4河流特征點中誤差 單位：米

從以上四個表可以看出：15米分辨率的DSM對資源三號影像的正射糾正效果是最好的，其他三種分辨率DEM的糾正效果隨著DEM分辨率的提高而線性降低。相同分辨率的DEM在地形起伏較大地區的糾正效果低于平坦地區的正射糾正效果。

（1）通過對比試驗可以得出：隨著使用的DEM精度的提升，資源三號影像正射糾正的精度也隨之提高，15米分辨率的數字表面模型對資源三號影像正射糾正的效果最好。

（2）試驗過程中，針對不同地貌類型的正射糾正結果可以看出：農田等地形起伏較小的地貌類型的正射糾正精度較好，山區等地形起伏較大的地貌類型糾正結果的精度有所下降，說明對地形復雜的遙感影像而言，其地形起伏直接影響正射糾正的精度，對于地形起伏較大的區域，高程模型的精度越高，其正射糾正效果越好。

（3）對于資源三號這種空間分辨率較高的遙感影像而言，使用45米分辨率的數字高程模型和90米分辨率的數字高程模型糾正的效果不好，因此不適合利用其對資源三號遙感影像做正射糾正。

3.結束語

本文以資源三號影像為基礎影像，利用四種不同分辨率的數字高程模型對其進行正射糾正處理，并對比其正射糾正的精度結果。本試驗中使用的標準影像是使用無控制點糾正的影像，如果能夠獲得GCP文件，將控制點參與到正射糾正中可以獲得更好的效果。每種地貌類別的糾正效果都是通過在糾正后的影像上手動截取一定區域然后與參照影像上手動截取相同區域提取特征點加以判斷，在手動截圖的過程中難免會產生一些偏差。有幾幅影像用SIFT算法提取特征點效果并不好，需要進行人工判斷并提取特征點，這個過程也會產生一些誤差。這些都是本次試驗中的不足。本次試驗僅對資源三號影像進行糾正分析，下一步將對多種分辨率的影像采用不同尺度的高程模型進行正射糾正，再對糾正結果進行對比分析，進一步驗證不同分辨率高程模型對影像正射糾正精度的影響。