第四章遙感圖像校正

四、遙感圖像的校正

（RSImageCorrection）影響遙感圖像的因素（influentialfactorofRSImage）：1、引起幾何畸變的因素

衛(wèi)星的姿態(tài)、高度、速度變化及其前進(jìn)運(yùn)動(dòng)；MSS掃描鏡掃描速度不均，檢測(cè)器采樣延遲誤差，波段間配準(zhǔn)誤差及全景畸變；地球自轉(zhuǎn)、曲率、高程的影響等2/1/202312、引起輻射失真的因素如：大氣的吸收、散射；地面及傳感器系統(tǒng)中的儀器在接收、轉(zhuǎn)換、傳送、處理圖像信息過(guò)程中性能不穩(wěn)；非線性和非一致性故障及不可避免引入的噪聲這些畸變和失真影響了圖像的質(zhì)量和應(yīng)用，必須進(jìn)行消除——圖像恢復(fù)。

圖像恢復(fù)（ImageRestoration）

又叫圖像復(fù)原，是指改正或補(bǔ)償在成像過(guò)程中造成的輻射失真、系統(tǒng)噪聲和隨機(jī)噪聲、幾何畸變以及高頻信息的損失。圖像恢復(fù)主要包括兩大部分：輻射校正（RadiometricCorrection）

幾何校正（GeometricCorrection）圖像恢復(fù)一）遙感圖像的輻射校正1、概念：

輻射校正(RadiometricCorrection)

是指消除圖像數(shù)據(jù)中依附在輻射亮度中的各種失真的過(guò)程。2、內(nèi)容：1）由傳感器的靈敏度特性引起的畸變校正2）由太陽(yáng)高度角及地形等引起的畸變校正3）由大氣的散射和吸收引起的輻射校正4.1.1光學(xué)鏡頭的非均勻性引起的邊緣減光現(xiàn)象的改正減光現(xiàn)象（DimLightPhenomenon）

在使用透鏡的光學(xué)系統(tǒng)中，由于鏡頭光學(xué)特性的非均勻性，在其成像平面存在著邊緣部分比中間部分暗的現(xiàn)象。N’ONPSDθf(wàn)f/cosθ鏡頭的輻射畸變?nèi)缟蠄D所示，如果光線以平行于主光軸的方向通過(guò)鏡頭到達(dá)像平面O點(diǎn)的光強(qiáng)度為Eo，以與主光軸成θ角的方向通過(guò)鏡頭到達(dá)像平面P點(diǎn)的光強(qiáng)度為Ep，則：Ep＝Eo·cos4θ。利用這一性質(zhì)可以進(jìn)行邊緣減光現(xiàn)象造成的輻射畸變校正。另外，對(duì)于視場(chǎng)較大的成像光譜儀圖像在掃描方向上也存在明顯的輻射亮度不均勻的現(xiàn)象。這種輻射誤差主要是光線路徑長(zhǎng)短不同造成的，掃描角越大時(shí)，光線路徑越長(zhǎng)，大氣衰減越嚴(yán)重。

星（機(jī)）下點(diǎn)位置的地物輻射信息的光線路徑最短，大氣衰減所產(chǎn)生的影響也最小。因此輻射量失真最小。4.1.2光電變換系統(tǒng)的特性引起的輻射誤差校正

傳感器的光譜響應(yīng)特性和傳感器的輸出有直接的關(guān)系。在掃描方式的傳感器中，傳感器接收系統(tǒng)收集到的電磁波信號(hào)需經(jīng)光電轉(zhuǎn)換系統(tǒng)變成電信號(hào)記錄下來(lái)，這一過(guò)程也會(huì)引起輻射量誤差。

光電變換系統(tǒng)的靈敏度特性有很高的重復(fù)性，故可以定期在地面測(cè)量其特性，根據(jù)測(cè)量值可以對(duì)其進(jìn)行輻射畸變校正。例如：對(duì)于Landsat衛(wèi)星的MSS圖像和TM圖像可以按下式對(duì)傳感器輸出的輻射亮度（R）進(jìn)行校正。

V=[Rmax/(Rmax-Rmin)]·(R-Rmin

)

式中R是傳感器輸出的輻射亮度；

V是已校正過(guò)的數(shù)據(jù)；

Rmax和Rmin分別為探測(cè)器能夠輸出的最大和最小輻射亮度。4.1.3太陽(yáng)高度引起的輻射誤差校正太陽(yáng)高度角引起的畸變校正是將太陽(yáng)光線傾斜照射時(shí)獲取的圖像校正為太陽(yáng)光線垂直照射時(shí)獲取的圖像。太陽(yáng)高度角θ可由下式計(jì)算：

sinθ=sinφ·sinδ±cosφ·cosδ·cost

式中，φ為圖像對(duì)應(yīng)的地理緯度，δ為太陽(yáng)赤緯(成像時(shí)太陽(yáng)直射點(diǎn)的地理緯度），t為時(shí)角（地區(qū)經(jīng)度與成像時(shí)太陽(yáng)直射點(diǎn)地區(qū)經(jīng)度的經(jīng)差）。

太陽(yáng)高度角的校正是通過(guò)調(diào)整一幅圖像內(nèi)的平均灰度來(lái)實(shí)現(xiàn)的，太陽(yáng)以高度角θ斜射時(shí)得到的圖像g（x，y）與直射時(shí)得到的圖像f（x，y)有如下關(guān)系：

f（x，y）＝g（x，y)/sinθ如果不考慮天空光的影響，則各波段可采用相同的θ角進(jìn)行校正。4.1.4地形坡度引起的輻射誤差校正太陽(yáng)光線和地表作用以后再反射到傳感器的太陽(yáng)光的輻射亮度和地面傾斜度有關(guān)。如果光線垂直入射時(shí)水平地表受到的光照強(qiáng)度為Io，則光線垂直入射時(shí)傾斜角為α的坡面上入射點(diǎn)處的光線強(qiáng)度I為：

I＝Io·cosα因此，若處在坡度為α的傾斜面上的地物影像為g（x，y），則校正后圖像f（x，y）為：

f（x，y）＝g（x，y)/cosα4.1.5大氣校正

大氣校正（AtmosphericCorrection）

Image-processingprocedurethatcompensatesforeffectsofselectivityscatteredlightinmultispectralimages.

是清除衛(wèi)星遙感圖象在大氣傳輸中所引起的退化因素的一種圖像處理方法。

在大氣影響中，我們主要研究大氣散射的影響。大氣散射是大氣中分子和微小粒子多次作用的結(jié)果。它隨電磁波的波長(zhǎng)和散射體的大小不同而不同。

大氣散射的校正一般有三種方法：1）公式計(jì)算法（Formula)；2）野外波譜測(cè)試回歸分析法(SpectrumTest)；

3）波段對(duì)照法(WavebandComparison)

公式計(jì)算法需知道具體天氣條件下大氣路徑輻射率等參數(shù)，野外波譜測(cè)試回歸分析法需要到野外進(jìn)行與陸地衛(wèi)星同步的一致測(cè)試，這些與波段對(duì)照法相比，都相對(duì)困難，因此，我們一般采用波段對(duì)照法。大氣散射的校正波段對(duì)比法大氣散射校正

大氣選擇散射對(duì)短波段影響大，長(zhǎng)波段影響小。因此，就Landsat來(lái)說(shuō)，4波段受散射影響最嚴(yán)重，其次為5、6波段，而7波段受影響最小，一般認(rèn)為不受影響。1.回歸分析法2.直方圖法改正數(shù)：－127波段與其他波段比較，找出最黑區(qū)域，再抽出幾個(gè)像元進(jìn)行回歸分析，得到一條擬合直線。1.回歸分析法（RegressionAnalysis）2.直方圖法（Histogram）二）遙感圖像的幾何畸變與幾何校正幾何畸變geometricdistortion

圖像中的幾何圖形與該物體在所選定投影中幾何圖形的差異。幾何畸變主要由于遙感器姿態(tài)角的變化，物鏡系統(tǒng)的光學(xué)畸變，掃描速度不穩(wěn)定，地球自轉(zhuǎn)，地球曲率，地面起伏及地圖投影等因素所引起。幾何畸變使圖像產(chǎn)生幾何形式或位置的失真，主要表現(xiàn)為仿射、偏扭、彎曲及其它更高階的歪曲。幾何畸變表現(xiàn)為像元相對(duì)地面實(shí)際位置發(fā)生擠壓、伸展、扭曲或偏移。圖像的幾何畸變可通過(guò)幾何校正消除。幾何畸變的原因：1）遙感器的內(nèi)部畸變：由遙感器結(jié)構(gòu)引起2）遙感器的外部畸變：由圖像投影方式的幾何學(xué)引起

可分為平臺(tái)引起的畸變目標(biāo)物引起的畸變（地球自轉(zhuǎn)）3）圖像坐標(biāo)系的定義方式不同引起的畸變4）地圖投影方法的不同引起的畸變2.幾何校正幾何校正（geometriccorrection）——是指從具有幾何畸變的圖像中消除畸變的過(guò)程，即定量地確定圖像上像元坐標(biāo)（圖像坐標(biāo)）與目標(biāo)物的地理坐標(biāo)（地圖坐標(biāo)等）的對(duì)應(yīng)關(guān)系（坐標(biāo)變換式）。幾何校正的過(guò)程未校正圖像（輸入圖像）a系統(tǒng)性校正b非系統(tǒng)性校正c復(fù)合校正①確定校正方法②確定校正式③檢驗(yàn)校正方法校正式的有效性幾何校正后圖像（輸出圖像）④重采樣，內(nèi)插

從遙感衛(wèi)星地面站得到的CCT圖像一般

均是經(jīng)過(guò)幾何粗校正處理的圖像。因此，下

面重點(diǎn)研究幾何精校正。幾何校正可分為：幾何粗校正——針對(duì)幾何畸變的原因進(jìn)行的幾何校正。幾何精校正——利用地面控制點(diǎn)進(jìn)行的幾何校正。1）幾何精校正的原理與方法

A.原理：幾何精校正是利用地面控制點(diǎn)（GroundControlPoint，簡(jiǎn)稱(chēng)GCP）進(jìn)行的。它通過(guò)GCP數(shù)據(jù)對(duì)原始衛(wèi)星圖像的幾何畸變過(guò)程進(jìn)行數(shù)學(xué)模擬，建立原始的畸變圖像空間與地理制圖用的標(biāo)準(zhǔn)空間（亦即校正空間）之間的某種對(duì)應(yīng)關(guān)系，然后利用這種對(duì)應(yīng)關(guān)系把畸變圖像空間中的全部元素變換到校正圖像空間中去，從而實(shí)現(xiàn)幾何精校正。校正空間在我國(guó)為高斯——克里格投影空間。

B.方法:設(shè)原始的二維圖象空間為F，

其縱坐標(biāo)軸為x（衛(wèi)星前進(jìn)方向），橫坐標(biāo)軸為y（掃描方向）；校正后二維圖象空間為G，

其縱坐標(biāo)軸為ξ（南北方向），橫坐標(biāo)軸為η（東西方向）；L是某種變換關(guān)系；

f（X，Y）為F空間中像元F（X，Y）的亮度值，

g（ξ，η）為該像元在G空間中對(duì)應(yīng)像元G（ξ，η）的亮度值通?？捎袃煞N方法實(shí)現(xiàn)這兩個(gè)空間的轉(zhuǎn)換（1）直接轉(zhuǎn)換首先求出原始圖象空間F中各像元F（X，Y）在G空間中的對(duì)應(yīng)位置G（ξ，η），然后將相應(yīng)的亮度值搬移過(guò)去。即G（ξ，η）＝L{F（X，Y）}f（X，Y）＝＞g（ξ，η）

（2）重采樣技術(shù)首先反求出校正圖象空間中各像元G（ξ，η）在原始圖象空間中的共軛點(diǎn)位置F（X，Y）見(jiàn)下圖，然后再利用某種方法（下面詳細(xì)談到）確定這一共軛點(diǎn)位置的亮度值f（X，Y），最后再將該亮度值搬移到校正空間中取。即F（X，Y）＝L－1G（ξ，η）

求f（X，Y）f（X，Y）＝＞ｇ（ξ，η）我們希望校正后的數(shù)字圖象是一個(gè)均勻（指像元分布均勻）的二維圖象空間矩陣，這在利用計(jì)算機(jī)的寬行打印機(jī)輸出時(shí)尤其是這樣。由于原始圖象空間與校正圖象空間之間的對(duì)應(yīng)是非線性的，則幾乎不可能兩個(gè)圖象空間的網(wǎng)格像元在各自的空間都均勻分布。在原始圖象空間中一組均勻分布的網(wǎng)點(diǎn)在校正圖象空間中對(duì)應(yīng)著一組不均勻分布的網(wǎng)點(diǎn)。因此在使用直接法時(shí)要特別注意這一點(diǎn)；重采樣法由于是從校正空間中的一組均勻分布網(wǎng)點(diǎn)出發(fā)進(jìn)行的。這一問(wèn)題自然就獲得了解決。2）幾何精校正的步驟A：幾何位置轉(zhuǎn)換B：亮度值確定A：幾何位置轉(zhuǎn)換數(shù)學(xué)模型的建立（1）三角形線性法（2）二元多項(xiàng)式法（1）三角形線性法圖象的幾何失真一般是非線性的，但是在一個(gè)小區(qū)域內(nèi)可以認(rèn)為是線性失真?；谶@一假設(shè),我們可以利用GCP構(gòu)成一系列的小三角形區(qū)域,并使它們的全體覆蓋整個(gè)待校正區(qū)域;利用原始空間和校正空間每一對(duì)三角形區(qū)域的三個(gè)GCP（處于三角形的三個(gè)頂點(diǎn)，——求出每對(duì)三角形范圍內(nèi)的轉(zhuǎn)換關(guān)系，進(jìn)行轉(zhuǎn)換。（2）二元多項(xiàng)式法該方法以二元n次多項(xiàng)式來(lái)描述幾何失真的情況，在多項(xiàng)式系數(shù)求解時(shí)，一般用最小二乘法。B：亮度值確定通常有三種方法：（1）最近鄰點(diǎn)法（2）雙向線性插值法（3）三次卷積法（1）最近鄰點(diǎn)法該方法取離共軛點(diǎn)最近的像元的亮度值作為該網(wǎng)格的亮度值。特點(diǎn)：計(jì)算量小，具有一定的精度，是一種簡(jiǎn)單而有效的方法。但處理后圖像的亮度具有明顯的不連續(xù)性。最近鄰點(diǎn)法又稱(chēng)一級(jí)內(nèi)插法，它使用內(nèi)插點(diǎn)周?chē)?個(gè)觀測(cè)點(diǎn)的像元值，對(duì)所求的像元值進(jìn)行線性?xún)?nèi)插。（2）雙線性?xún)?nèi)插法特點(diǎn)：具有圖像均衡化和清晰化的效果，可得到較高的圖像質(zhì)量，但破壞了原有數(shù)據(jù)，使高頻分量受損，圖像輪廓模糊。

（3）三次卷積法利用多項(xiàng)式來(lái)逼近理論上的最佳插值函數(shù)函數(shù)。特點(diǎn)：計(jì)算量大，但可以克服前種方法的缺點(diǎn)，而且精度較高。2）幾何精校正的幾個(gè)關(guān)鍵問(wèn)題A：建立原始圖像空間坐標(biāo)系B：建立校正圖像空間坐標(biāo)系C：GCP的確定D：圖像的輸出A：建立原始圖像空間坐標(biāo)系B：建立校正圖像空間坐標(biāo)系該坐標(biāo)系的建立很簡(jiǎn)單，取要校正的數(shù)字圖象最左上角的點(diǎn)（象元）作為坐標(biāo)原點(diǎn)。掃描線號(hào)作為x軸坐標(biāo)值，采樣號(hào)數(shù)作為y軸坐標(biāo)值。校正圖象空間坐標(biāo)原點(diǎn)仍取在圖象的左上角，以ξ表示縱方向（南北方向）坐標(biāo)，η表示橫方向（東西方向）坐標(biāo)，坐標(biāo)以校正空間的網(wǎng)格為單位。網(wǎng)格大小由具體的圖象輸出設(shè)備和輸出圖象的比例尺來(lái)確定，網(wǎng)格大小就是以后重采樣的間隔。C：GCP的確定（1）GCP的對(duì)象

GCP應(yīng)是在原始CCT圖像中可尋找出來(lái)的，在地形圖上可精確定位的一些特征點(diǎn)、特征線等自然要素和人文要素。（2）GCP的檢測(cè)定位

人工目視檢測(cè)法；計(jì)算機(jī)自動(dòng)檢測(cè)定位（3）GCP的數(shù)量

適量增加可提高幾何校正精度，但不宜過(guò)多。二元n次多項(xiàng)式：L≥（n＋1）（n＋2）/2（4）GCP的分布均勻（5）GCP的精度

GCP的位置精度越高，幾何校正的效果越好；對(duì)于三角形線性法，當(dāng)GCP不落在三角形頂點(diǎn)時(shí)，影響不大。對(duì)于二次多項(xiàng)式擬合法，當(dāng)L≥（n＋1）（n＋2）/2時(shí)，影響不大。D：圖像的輸出寬行打印機(jī)：數(shù)字圖像繪圖機(jī)/電視顯示屏：模擬圖像數(shù)字圖象的優(yōu)點(diǎn)便于進(jìn)行定量分析缺點(diǎn)不直觀模擬圖象直觀醒目，常將之與地形圖進(jìn)行比較，評(píng)價(jià)校正精度。三）遙感圖像的配準(zhǔn)1、多圖像——