摘要：
針對(duì)推掃式成像光譜儀獲取的狹帶影像需要經(jīng)過(guò)幾何校正才能拼接形成空間二維影像的問(wèn)題，提出了基于ENVI二次開(kāi)發(fā)的高光譜推掃圖像拼接技術(shù)。基于單應(yīng)映射建立光譜儀傾斜狀態(tài)下與正射狀態(tài)下圖像上的二維點(diǎn)之間的關(guān)系，校正由姿態(tài)變化引起的圖像畸變，結(jié)合GPS數(shù)據(jù)修正因飛行速度變化引起的狹帶重疊將校正后的狹帶影像拼接起來(lái)。在ENVI二次開(kāi)發(fā)平臺(tái)上進(jìn)行技術(shù)集成，實(shí)現(xiàn)了Resonon推掃高光譜狹帶影像的自動(dòng)校正拼接。對(duì)河北保定郊區(qū)高光譜影像的校正拼接實(shí)驗(yàn)證明，該方法與光譜儀自帶拼接軟件校正結(jié)果接近經(jīng)緯度坐標(biāo)差均在1m以內(nèi)，均方根誤差約為0.7389，能夠滿足一般高光譜遙感應(yīng)用中的地理精度要求。
研究目的：
根據(jù)單應(yīng)映射原理，建立光譜儀傾斜和正射狀態(tài)下像點(diǎn)的映射關(guān)系，利用GPS/INS組合導(dǎo)航數(shù)據(jù)校正狹帶影像中的畸變，拼接成一幅完整的影像，并在ENVI二次開(kāi)發(fā)平臺(tái)上實(shí)現(xiàn)推掃狹帶影像的自動(dòng)校正和拼接。
推掃成像畸變?cè)颍?/span>
推掃式成像是利用飛行平臺(tái)的向前運(yùn)動(dòng)，借助于與飛行方向垂直的掃描線記錄而構(gòu)成二維圖像。推掃型成像光譜儀通常采用一個(gè)垂直于運(yùn)動(dòng)方向的面陣CCD來(lái)感應(yīng)地面響應(yīng)，在飛行平臺(tái)向前運(yùn)動(dòng)中完成二維空間掃描，平行于平臺(tái)運(yùn)動(dòng)方向，通過(guò)光柵和棱鏡分光完成光譜維掃描，因此，CCD上一個(gè)點(diǎn)對(duì)應(yīng)一個(gè)譜段，一條線對(duì)應(yīng)一個(gè)譜面。CCD探測(cè)器每次成像是空間一條線上的光譜信息。為了獲得空間二維圖像，再通過(guò)機(jī)械推掃，完成整個(gè)平面的圖像和光譜數(shù)據(jù)采集。
推掃成像時(shí)，CCD探測(cè)器所記錄的高光譜圖像數(shù)據(jù)是沿著飛行方向的條幅，由于搭載光譜儀的飛行平臺(tái)在飛行過(guò)程中，不能一直保證理想的姿態(tài)正射獲取影像，速度和姿態(tài)的不穩(wěn)定導(dǎo)致飛行平臺(tái)的位置、航偏角、俯仰角和橫滾角不斷隨機(jī)變化，引起光譜儀拍攝時(shí)外方位元素也不斷隨機(jī)變化。因此，CCD曝光時(shí)每條掃描線對(duì)應(yīng)的光譜儀外方位元素不一致引起了圖像的幾何畸變：
1. 飛行平臺(tái)姿態(tài)不穩(wěn)定造成地面掃描行之間相互交錯(cuò)，圖像扭曲變形，影響后期地物目標(biāo)的解析和判別。
2. 飛行平臺(tái)速度不穩(wěn)定易造成掃描行之間的行間距忽大忽小，出現(xiàn)重疊或間隙，為了獲得地面的完整影像，通常推掃成像需保證一定的采樣率。因此，在圖像拼接時(shí)就需要借助ＧＰＳ位置信息對(duì)重疊的掃描行進(jìn)行幾何糾正和圖像融合處理。
IDL實(shí)現(xiàn)
IDL是美國(guó)ITT VIS公司推出的第四代交互式、跨平臺(tái)、面向矩陣處理的編程語(yǔ)言，具有快速的數(shù)據(jù)分析、圖像處理和強(qiáng)大的可視化功能。采用IDL語(yǔ)言調(diào)用ENVI平臺(tái)中的圖像處理函數(shù)，可以很方便地進(jìn)行二次開(kāi)發(fā)，實(shí)現(xiàn)遙感數(shù)據(jù)的快速分析和可視化。
推掃圖像的自動(dòng)拼接主要包括如下３個(gè)基本步驟：
(１)影像和GPS/INS數(shù)據(jù)讀?。哼b感影像數(shù)據(jù)包含圖像本身和頭文件，ENVI二次開(kāi)發(fā)提供了函數(shù)讀取遙感影像及其屬性。如ENVI_OPEN_FILE、ENVI_FILE_QUERY、ENVI_GET_ SLICE等。GPS/INS數(shù)據(jù)存儲(chǔ)于文本文件中，按照文本文件讀取方式即可獲得狹帶影像獲取時(shí)光譜儀的姿態(tài)和位置信息。
(２)單應(yīng)矩陣計(jì)算和單應(yīng)映射：以北東地坐標(biāo)系為地理坐標(biāo)系，依據(jù)公式計(jì)算得到單應(yīng)矩陣Ｈ。主要代碼命令如下:
H=M_inv#MATRIX_POWER（C, -1）#M；計(jì)算單應(yīng)矩陣.
所以，對(duì)于狹帶影像上的每一個(gè)二維點(diǎn) ( xｂ，yｂ ) ，都可以獲得校正后的對(duì)應(yīng)點(diǎn) (xｎ，yｎ ) ，點(diǎn) ( xｎ，yｎ ) 的灰度值即為點(diǎn) ( xｂ，yｂ ) 的灰度值。
(３)圖像拼接：校正后的每條狹帶圖像中心點(diǎn)的二維地理坐標(biāo)即光譜儀成像中心的ＧＰＳ 二維坐標(biāo)，根據(jù)光譜儀的成像地面分辨率，選定影像投影方式，可以為每條狹帶設(shè)置地理信息。主要代碼命令如下:
map_info=ENVI_MAP_INFO_CREATE( /geographic, mc=mc, ps=ps ); 為狹帶添加地理信息。
拼接后的影像被認(rèn)為是光譜儀理想姿態(tài)下獲取的正射影像，具有與GPS獲取的一致的位置信息，拼接影像點(diǎn)的高光譜曲線與原始掃描行對(duì)應(yīng)點(diǎn)的一致，能夠真實(shí)地反映地面的空間特征和光譜特征。
實(shí)驗(yàn)結(jié)果與分析：
本文選擇河北省保定市郊區(qū)的高光譜影像進(jìn)行校正拼接實(shí)驗(yàn)，影像由搭載于無(wú)人機(jī)的PikaL 高光譜成像儀拍攝獲取，PikaL高光譜成像儀由美國(guó)Resonon公司設(shè)計(jì)生產(chǎn)，光譜范圍為400-1000nm，光譜分辨率為2.1 nm，CCD掃描行寬度為900像素。飛行過(guò)程中同時(shí)搭載慣導(dǎo)系統(tǒng)實(shí)時(shí)獲取光譜儀的姿態(tài)位置信息。高光譜儀將推掃獲取的原始狹帶影像先簡(jiǎn)單拼接起來(lái)，存儲(chǔ)于固態(tài)硬盤(pán)中，此時(shí)的地理信息并未經(jīng)過(guò)糾正，圖像存在幾何畸變，圖ａ所示為原始圖像的假彩色圖像，狹帶經(jīng)過(guò)幾何校正和拼接后才能正確顯示地面目標(biāo)的特征，如圖ｂ所示。
為了能夠定量檢驗(yàn)該幾何校正方法的效果，同時(shí)采用Pika L高光譜成像儀自帶的軟件對(duì)原始影像進(jìn)行幾何校正，將兩種方法得到的校正影像進(jìn)行比較。兩種校正方法均采用ＵＴＭ 投影，以WGS-84為基準(zhǔn)面。首先在軟件校正影像中隨機(jī)選取10個(gè)均勻分布的明顯地物點(diǎn)，讀取其坐標(biāo)值，作為采樣點(diǎn)用于評(píng)定校正精度，然后從本文方法校正后影像中讀取其相應(yīng)坐標(biāo)值，經(jīng)過(guò)對(duì)10個(gè)采樣點(diǎn)殘差的計(jì)算得到如表１所示的精度檢驗(yàn)結(jié)果。
表１：北向距離均方根誤差為0.6327ｍ，東向距離均方根誤差為0.3817ｍ，
總均方根誤差為0.7389ｍ
由表１可以看出，采樣點(diǎn)在ｘ和ｙ方向上的坐標(biāo)偏移均不超過(guò)１ｍ，兩種方法得到的校正圖像地理信息較為接近；ｙ方向坐標(biāo)均方根誤差大于ｘ方向坐標(biāo)均方根誤差，即像點(diǎn)坐標(biāo)的經(jīng)度值準(zhǔn)確性高于緯度值。對(duì)于某些地理精度要求不高的航空高光譜遙感應(yīng)用來(lái)說(shuō)，本方法取得的校正效果已滿足需求。如果需要進(jìn)一步提高精度，可以通過(guò)增加地面控制點(diǎn)或與高精度地圖進(jìn)行圖像配準(zhǔn)實(shí)現(xiàn)幾何精校正。
結(jié)語(yǔ)
本文根據(jù)推掃成像和單應(yīng)映射原理，結(jié)合GPS/INS組合導(dǎo)航系統(tǒng)實(shí)時(shí)獲取光譜儀姿態(tài)角度和位置信息。在ENVI二次開(kāi)發(fā)平臺(tái)上，采用IDL語(yǔ)言實(shí)現(xiàn)了高光譜儀推掃狹帶影像的自動(dòng)校正和拼接。驗(yàn)證實(shí)驗(yàn)表明，本方法與自帶軟件校正拼接效果接近，均方根誤差基本滿足一般的高光譜遙感應(yīng)用。雖然本文方法能夠取得較為理想的校正拼接效果，但是單掃描行的校正過(guò)程耗時(shí)較長(zhǎng)，無(wú)法實(shí)時(shí)獲取校正影像，下一步將就提高校正拼接效率展開(kāi)更加深入的研究。另外，拼接過(guò)程中不同成像條件下的勻色處理同樣是后續(xù)需要研究的內(nèi)容。
原文下載連接：
基于ENVI二次開(kāi)發(fā)的高光譜圖像拼接與幾何校正評(píng)估