胡九超 周忠發(fā)

摘要：提出了一種快速獲取大量地面控制點的方法，以TM影像為基準(zhǔn)圖像選擇初始同名點，經(jīng)過自動計算誤差、調(diào)整、平移誤差消除、刪除等步驟獲取穩(wěn)定的其他同名點，并采用線性橡皮拉伸進(jìn)行幾何精校正。試驗證明，利用該方法進(jìn)行幾何精校正后的HJ-1CCD影像不但加快了選點速度而且能保證精度，可以用于遙感信息的提取以及為地理信息系統(tǒng)等提供可用的數(shù)據(jù)。

關(guān)鍵詞：貴州高原山區(qū)；HJ-1CCD影像；地面控制點；幾何精校正

中圖分類號：P231 文獻(xiàn)標(biāo)識碼：A 文章編號：0439-8114（2015）03-0703-02

DOI：10.14088/j.cnki.issn0439-8114.2015.03.052

HJ-1CCD Image Geometric Accurate Rectification based on Fast Matching and

Linear Rubber Sheeting

HU Jiu-chao1，2， ZHOU Zhong-fa1，2

（1.Institute of South China Karst，Guizhou Normal University， Guiyang 550001， China；

2. The State Key Laboratory Incubation Base for Karst Mountain Ecology Environment of Guizhou Province， Guiyang 550001， China）

Abstract：The study put forward a rapid method to access a large number of ground control points， based on TM image as the reference image selection of initial corresponding points， through the steps of automatic calculation error， adjustment， translation error elimination， and deletion， obtained other stable matching points， and used linear rubber sheeting for geometric correction. Experiments showed that， HJ-1CCD image after geometric correction not only accelerated the speed and position but also ensured accuracy， which can be used to provide data available for extracting the information of remote sensing and geographic information system.

Key words： plateau mountain area of Guizhou； HJ-1CCD image； ground control point； geometric accurate correction

進(jìn)入新世紀(jì)以來，遙感事業(yè)蓬勃發(fā)展，當(dāng)然中國也不例外。2008年9月6日，中國新一代民用光學(xué)對地遙感小衛(wèi)星——環(huán)境減災(zāi)衛(wèi)星成功發(fā)射升空。衛(wèi)星在軌組網(wǎng)，具有嚴(yán)格的成像視場匹配要求，可以在48 h內(nèi)對中國國境及周邊地區(qū)實現(xiàn)無縫覆蓋觀測[1]。隨著科學(xué)技術(shù)和社會需求的迅猛發(fā)展，遙感技術(shù)作為高新技術(shù)領(lǐng)域中的一個重要分支逐步邁向產(chǎn)業(yè)化、實用化，對遙感圖像產(chǎn)品的處理精度要求也越來越高[2]。然而，由于地球自轉(zhuǎn)、地面曲率的變化以及衛(wèi)星本身等方面的原因，遙感圖像存在無法避免的幾何畸變，而僅僅經(jīng)過系統(tǒng)粗校正的遙感圖像還無法滿足精度需求，所以，需要進(jìn)一步進(jìn)行遙感影像的幾何精校正。

遙感影像的傳統(tǒng)幾何精校正一般都通過地形圖手動采集地面控制點（Ground Control Points，GCPs）來進(jìn)行，但地形圖的測繪年代相對來說比較久遠(yuǎn)，信息量少，因而對地面控制點的采集具有很大的限制；其次，HJ-1CCD影像幅寬大，畸變復(fù)雜，誤差不一，采用全局校正模型就會有一定的局限性。綜合考慮上述問題，提出利用ERDAS IMAGINE 9.2中的AutoSync模塊快速獲取地面控制點，并采用線性橡皮拉伸校正模型對HJ-1CCD影像進(jìn)行幾何精校正，不但加快了選點速度而且有效提高了環(huán)境衛(wèi)星影像的定位精度。

1 試驗區(qū)概況及數(shù)據(jù)情況

貴州地處中國西南部，介于北緯24°37′—29°13′，東經(jīng)103°36′—109°35′，東毗湖南，南鄰廣西，西連云南，北接四川和重慶。氣候溫暖濕潤，屬亞熱帶濕潤季風(fēng)氣候區(qū)，氣溫變化小，冬暖夏涼，氣候宜人，貴州省大部分地區(qū)年平均氣溫15 ℃左右。貴州地處云貴高原東部，境內(nèi)地勢復(fù)雜多樣，地勢西高東低，自中部向北、東、南三面傾斜，平均海拔1 100 m左右。貴州省地貌可概括分為高原山地、丘陵和盆地三種基本類型，境內(nèi)山脈眾多，重巒疊峰，綿延縱橫，山高谷深，其中92.5%的面積為山地和丘陵，復(fù)雜的地形條件給影像的校正造成了干擾。

研究采用的遙感數(shù)據(jù)為HJ-1CCD影像，環(huán)境與災(zāi)害監(jiān)測預(yù)報小衛(wèi)星星座由兩顆光學(xué)小衛(wèi)星（HJ-1A，HJ-1B）和一顆合成孔徑雷達(dá)小衛(wèi)星（HJ-1C）構(gòu)成。其中HJ-1A和HJ-1B均采用準(zhǔn)太陽同步圓軌道，軌道高度為649.093 km，相位為180 b。這兩顆光學(xué)小衛(wèi)星上均裝有兩臺寬覆蓋多光譜可見光相機(jī)（簡稱CCD相機(jī)），單臺CCD相機(jī)幅寬為360 km（兩臺為710 km），共4個波段，波譜范圍為0.43～0.90 μm[3]。研究采用的數(shù)據(jù)是2011年5月28日采集的HJ-1衛(wèi)星影像以及同時相的TM影像。

2 影像的幾何校正

2.1 控制點的采集

研究采用線性橡皮拉伸模型進(jìn)行擬合，按原理而言，該方法需要大量的控制點，但傳統(tǒng)的選點方法顯然無法滿足要求，因此試驗借助分辨率為30 m的TM影像，在ERDAS IMAGINE 9.2的AutoSync模塊下進(jìn)行控制點的采集，因為該方法可以在短暫的時間內(nèi)生成大量同名控制點、且高于純手動配準(zhǔn)的精度，同時能夠處理各種各樣的數(shù)據(jù)等優(yōu)點。利用此方法采集控制點程序簡單，首先手動選取幾個明顯的地面特征點，譬如道路交叉口、河流交匯處等，且盡量分布均勻；然后利用軟件自動生成大量的同名點；緊接著對誤差大的控制點進(jìn)行調(diào)整、剔除、重選，直至符合精度要求，至于控制點稀疏區(qū)域要手動加一些控制點，盡量做到控制點在整個試驗區(qū)均勻地分布，否則，控制點分布密集的區(qū)域影像幾何校正精度高，而在稀疏區(qū)域擬合不好使得圖像產(chǎn)生變形。流程圖如圖1所示。

2.2 遙感影像重采樣

重新定位后的像元在原圖像中的分布是不均勻的，即輸出圖像像元點在輸入圖像中的行列號不是或不全是整數(shù)關(guān)系。因此需要根據(jù)輸出圖像上的各像元在輸入圖像中的位置，對原始圖像按一定規(guī)則重新采樣，進(jìn)行亮度值的插值計算，建立新的圖像矩陣[4]。一般常用的重采樣方法有三種：最鄰近像元采樣法、雙線性內(nèi)插法、三次卷積重采樣法。這三種方法各有優(yōu)劣，綜合考慮遙感影像和研究區(qū)的情況，試驗采用的是三次卷積重采樣法，雖然計算量比較大，但圖像灰度具有連續(xù)性且采樣精度高。

2.3 精度分析

RMS（均方根）誤差是GCP（地面控制點）的輸入（原）位置和逆轉(zhuǎn)換的位置之間的距離，它是在用轉(zhuǎn)換矩陣對一個GCP做轉(zhuǎn)換時所期望輸出的坐標(biāo)與實際輸出的坐標(biāo)之間的偏差[5]。殘差是某個方向上原作標(biāo)和逆變換坐標(biāo)之間的距離。每個GCP的殘差都可以顯示出來。X殘差是原X坐標(biāo)和逆變換的X坐標(biāo)之間的距離，Y殘差是原Y坐標(biāo)和逆變換Y坐標(biāo)之間的距離[6]。為了計算，每個控制點的RMS誤差都要給出，由距離公式計算：

Ri=■ 公式（1）

式中，Ri為GCPi的RMS誤差，XRi為GCPi的X殘差，YRi為GCPi的Y殘差。

Rx=■ 公式（2）

Ry=■ 公式（3）

T=■ 公式（4）

T為總的RMS誤差，n為控制點的個數(shù)。對于每次影像校正，不僅要驗證每個控制點的誤差Ri，而且要驗證總的誤差T。若精度沒有達(dá)到要求，則要反復(fù)試驗，直至符合精度為止。

3 小結(jié)與討論

用ERDAS IMAGINE 9.2中的AutoSync模塊選取地面控制點進(jìn)行幾何精校正方便快捷，不需人為干預(yù)，自動產(chǎn)生多個同名點，大大減少了幾何精校正人力的付出。選取貴州地區(qū)為研究區(qū)，并采用此法對環(huán)境衛(wèi)星影像進(jìn)行了幾何精校正，研究表明，誤差可以控制在1個像元左右，滿足要求，方法可行，同時得出以下幾點結(jié)論：

1）用TM影像輔助進(jìn)行HJ-1衛(wèi)星的幾何精校正，可以達(dá)到預(yù)期的精度，能夠滿足一般的工作要求；

2）采用AutoSync模塊選取地面控制點的方法，克服了傳統(tǒng)方法選點難、速度慢、精度低等缺點，使環(huán)境衛(wèi)星影像能夠更快、更好地與數(shù)字化地理信息相結(jié)合。

3）采用線性橡皮拉伸（Liner Rubber Sheeting）糾正法、三次卷積重采樣（Cubic Convolution）的方法對HJ-1衛(wèi)星影像的幾何精校正在該研究區(qū)范圍內(nèi)取得了較理想的效果，同時加快了幾何精校正的處理速度。

參考文獻(xiàn)：

[1] 白照廣，沈 中，王肇宇.環(huán)境減災(zāi)-1A、1B衛(wèi)星技術(shù)[J].航天器工程，2009，18（6）：1-11.

[2] 邵鴻飛，孔慶欣.遙感圖像幾何校正的實現(xiàn)[J].氣象，2000，26（2）：41-44.

[3] 任 平，楊存建，周介銘.HJ-1A/B星CCD多光譜遙感數(shù)據(jù)特征評價及應(yīng)用研究[J].遙感技術(shù)與應(yīng)用，2010，25（1）：138-142.

[4] 趙英時.遙感應(yīng)用分析原理與方法[M].北京：科學(xué)出版社，2003.

[5] 趙書河，馮學(xué)智，趙 銳.中巴資源一號衛(wèi)星南京幅數(shù)據(jù)質(zhì)量與幾何糾正評價[J].遙感技術(shù)與應(yīng)用，2000，15（3）：170-174.

[6] 劉志麗，陳 曦.基于ERDAS IMAGING軟件的TM影像幾何精校正方法初探——以塔里木河流域為例[J].干旱區(qū)地理，2001， 24（4）：353-358.