第四章二圖像校正

第四章二圖像校正第一頁，共五十頁，編輯于2023年，星期五圖像校正從具有畸變的圖像中消除畸變的處理過程叫圖像校正，包括：輻射校正幾何校正第二頁，共五十頁，編輯于2023年，星期五4.3輻射校正4.3.1輻射畸變進入傳感器的輻射強度反映在圖像上就是亮度值(灰度值)。該值主要受兩個物理量影響：A.太陽輻射照射到地面的輻射強度;B.地物的光譜反射率。

當(dāng)太陽輻射相同時，圖像上像元亮度值的差異直接反映了地物目標光譜反射率的差異。第三頁，共五十頁，編輯于2023年，星期五輻射強度值還受到其他因素的影響而發(fā)生改變。這一改變的部分就是需要校正的部分，稱為輻射畸變。輻射校正（radiometriccorrection）：消除由于輻射畸變所引起的輻射量失真的過程。第四頁，共五十頁，編輯于2023年，星期五引起輻射畸變的原因：A.傳感器儀器本身產(chǎn)生的誤差；B.大氣對輻射的影響。第五頁，共五十頁，編輯于2023年，星期五儀器引起的誤差是由于多個檢測器之間存在差異，以及儀器系統(tǒng)工作產(chǎn)生的誤差，這導(dǎo)致了接收的圖像不均勻，產(chǎn)生條紋和“噪聲”。一般來說，這種畸變應(yīng)該在數(shù)據(jù)生產(chǎn)過程中，由生產(chǎn)單位根據(jù)傳感器參數(shù)進行校正，而不需要用戶自行校正。用戶應(yīng)該考慮的是大氣影響造成的畸變。第六頁，共五十頁，編輯于2023年，星期五4.3.2大氣影響的定量分析進入大氣的太陽輻射會發(fā)生反射、折射、吸收、散射和透射。其中對傳感器接收影響較大的是吸收和散射。

第七頁，共五十頁，編輯于2023年，星期五1.無大氣：在沒有大氣存在時，傳感器接收的輻照度，只與太陽輻射到地面的輻照度和地物反射率有關(guān)。設(shè)E0λ為波長λ的入射輻照度，θ為入射方向的天頂角，當(dāng)無大氣存在時，地面上單位面積的輻照度為：

第八頁，共五十頁，編輯于2023年，星期五假定地表面是朗伯體，其表面為漫反射，則某方向物體的亮度為：是地物反射率；是球面度(半球反射)。

第九頁，共五十頁，編輯于2023年，星期五傳感器接收信號時，受儀器的影響還有一個系統(tǒng)增益系數(shù)因子，這時進入傳感器的亮度值為：第十頁，共五十頁，編輯于2023年，星期五2.有大氣:由于大氣的存在，輻射經(jīng)過大氣吸收和散射，透過率小于1，從而減弱了原信號的強度。同時大氣的散射光也有一部分直接或經(jīng)過地物反射進入到傳感器，這兩部分輻射又增強了信號，但卻不是有用的。在入射方向有與入射天頂角θ和波長λ有關(guān)的透過率Tθλ；反射后，在反射方向上有與反射大頂角Φ和波長λ又有關(guān)的透過率TΦλ。因此進入傳感器的亮度值為

第十一頁，共五十頁，編輯于2023年，星期五大氣對輻射散射后，來自各個方向的散射又重新以漫入射的形式照射地物，其輻照度為ED，經(jīng)過地物的反射及反射路徑上大氣的吸收進入傳感器，其亮度值為（此值通常很小，有人主張忽略不計）第十二頁，共五十頁，編輯于2023年，星期五相當(dāng)部分的散射光向上通過大氣直接進入傳感器，這部分輻射稱為程輻射度，亮度為。第十三頁，共五十頁，編輯于2023年，星期五由于大氣影響的存在，實際到達傳感器的輻射亮度是前面所分析的三項之和，即第十四頁，共五十頁，編輯于2023年，星期五比較以下兩個公式：

大氣的主要影響是減少了圖像的對比度，使原始信號和背景信號都增加了因子。

第十五頁，共五十頁，編輯于2023年，星期五4.2.3大氣影響的粗略糾正精糾正：嚴格地說，去除大氣影響是將公式

中的附加項和附加因子求出，最終求出地物反射率R，從而恢復(fù)遙感影像中地面目標的真實面目。

第十六頁，共五十頁，編輯于2023年，星期五粗略糾正：

粗略校正指通過比較簡便的方法去掉式

中的Lp，即程輻射度，從而改善圖像質(zhì)量。式中還有漫入射因子ED及其他如透過率等影響，這些因子都作為地物反射率的因子出現(xiàn)，直接相減不易去除，常用比值法或其他校正方法去除。嚴格地說，程輻射度的大小與像元位置有關(guān)，隨大氣條件、太陽照射方向和時間變化而變化，但因其變化量微小而忽略。可以認為，程輻射度在同一幅圖像的有限面積內(nèi)是一個常數(shù)，其值的大小只與波段有關(guān)。

第十七頁，共五十頁，編輯于2023年，星期五大氣影響的粗略糾正方法1.直方圖最小值去除法

數(shù)字圖像直方圖第十八頁，共五十頁，編輯于2023年，星期五基本思想

一幅圖像中總可以找到某種或某幾種地物，其輻射亮度或反射率接近0，例如，地形起伏地區(qū)山地陰影處，反射率極低的深海水體處等，這時在圖像中對應(yīng)位置的像元亮度值應(yīng)為0。實測表明，這些位置上的像元亮度不為零。這個值就應(yīng)該是大氣散射導(dǎo)致的程輻射度值。

第十九頁，共五十頁，編輯于2023年，星期五具體校正方法：

1.確定程輻射度增值：即該圖像上確有輻射亮度或反射亮度應(yīng)為零的地區(qū)，則亮度最小值必定是這一地區(qū)大氣影響的程輻射度增值。

2.做減法運算：校正時，將每一波段中每個像元的亮度值都減去本波段的最小值。使圖像亮度動態(tài)范圍得到改善，對比度增強，從而提高了圖像質(zhì)量。第二十頁，共五十頁，編輯于2023年，星期五2.回歸分析法

一般來說由于程輻射度主要來自米氏散射，其散射強度隨波長的增大而減小，到紅外波段也有可能接近于零第二十一頁，共五十頁，編輯于2023年，星期五回歸分析法

假定某紅外波段，存在程輻射為主的大氣影響，且亮度增值最小，接近于零，設(shè)為波段a?，F(xiàn)需要找到其他波段相應(yīng)的最小值，這個值一定比a波段的最小值大一些，設(shè)為波段b，分別以a，b波段的像元亮度值為坐標，作二維光譜空間，兩個波段中對應(yīng)像元在坐標系內(nèi)用一個點表示。由于波段之間的相關(guān)性，通過回歸分析在眾多點中一定能找到一條直線與波段b的亮度Lb軸相交，且

第二十二頁，共五十頁，編輯于2023年，星期五是斜率：和是分別為a、b波段亮度的平均值。是波段a中的亮度為0處波段b中所具有的亮度?？梢哉J為就是波段b的程輻射度。校正的方法是將波段b中每個像元的亮度值減去，來改善圖像，去掉程輻射。同理可依次完成其他波段的校正。第二十三頁，共五十頁，編輯于2023年，星期五4.2.4幾何校正當(dāng)遙感圖像在幾何位置上發(fā)生了變化，產(chǎn)生諸如行列不均勻，像元大小與地面大小對應(yīng)不準確，地物形狀不規(guī)則變化等畸變時，即說明遙感影像發(fā)生了幾何畸變。遙感影像的總體變形（相對于地面真實形態(tài)而言）是平移、縮放、旋轉(zhuǎn)、偏扭、彎曲及其他變形綜合作用的結(jié)果。產(chǎn)生畸變的圖像給定量分析及位置配準造成困難，因此遙感數(shù)據(jù)接收后，首先由接收部門進行校正，這種校正往往根據(jù)遙感平臺、地球、傳感器的各種參數(shù)進行處理。而用戶拿到這種產(chǎn)品后，由于使用目的不同或投影及比例尺的不同，仍舊需要作進一步的幾何校正。

第二十四頁，共五十頁，編輯于2023年，星期五1.遙感影像變形的原因

遙感平臺位置和運動狀態(tài)變化的影響

地形起伏的影響地球表面曲率的影響大氣折射的影響地球自轉(zhuǎn)的影響

第二十五頁，共五十頁，編輯于2023年，星期五(1)遙感平臺位置和運動狀態(tài)變化的影響(航高、航速、俯仰、翻滾、偏航）

無論是衛(wèi)星還是飛機，運動過程中都會由于種種原因產(chǎn)生飛行姿勢的變化從而引起影像變形。航高：當(dāng)平臺運動過程中受到力學(xué)因素影響，產(chǎn)生相對于原標準航高的偏離，或者說衛(wèi)星運行的軌道本身就是橢圓的。航高始終發(fā)生變化，而傳感器的掃描視場角不變，從而導(dǎo)致圖像掃描行對應(yīng)的地面長度發(fā)生變化。航高越向高處偏離，圖像對應(yīng)的地面越寬

第二十六頁，共五十頁，編輯于2023年，星期五航速：衛(wèi)星的橢圓軌道本身就導(dǎo)致了衛(wèi)星飛行速度的不均勻，其他因素也可導(dǎo)致遙感平臺航速的變化。航速快時，掃描帶超前，航速慢時，掃描帶滯后，由此可導(dǎo)致圖像在衛(wèi)星前進方向上（圖像上下方向）的位置錯動。

第二十七頁，共五十頁，編輯于2023年，星期五俯仰：遙感平臺的俯仰變化能引起圖像上下方向的變化，即星下點俯時后移，仰時前移，發(fā)生行間位置錯動。

第二十八頁，共五十頁，編輯于2023年，星期五翻滾：遙感平臺姿態(tài)翻滾是指以前進方向為軸旋轉(zhuǎn)了一個角度?？蓪?dǎo)致星下點在掃描線方向偏移，使整個圖像的行向翻滾角引起偏離的方向錯動。

第二十九頁，共五十頁，編輯于2023年，星期五偏航：指遙感平臺在前進過程中，相對于原前進航向偏轉(zhuǎn)了一個小角度，從而引起掃描行方向的變化，導(dǎo)致圖像的傾斜畸變。第三十頁，共五十頁，編輯于2023年，星期五(2)地形起伏的影響

當(dāng)?shù)匦未嬖谄鸱鼤r，會產(chǎn)生局部像點的位移，使原來本應(yīng)是地面點的信號被同一位置上某高點的信號代替。由于高差的原因，實際像點P距像幅中心的距離相對于理想像點P0距像幅中心的距離移動了△r。高差引起的像點位移第三十一頁，共五十頁，編輯于2023年，星期五(3)地表曲率的影響

地球是球體，嚴格說是橢球體，因此地球表面是曲面。這一曲面的影響主要表現(xiàn)在兩個方面，一是像點位置的移動，當(dāng)選擇的地圖投影平面是地球的切平面時，使地面點P0相對于投影平面點P有一高差△h。像點位移第三十二頁，共五十頁，編輯于2023年，星期五

二是像元對應(yīng)于地面寬度的不等。由于傳感器通過掃描取得數(shù)據(jù)，在掃描過程中每一次取樣間隔是星下視場角的等分間隔。如果地面無彎曲，在地面瞬時視場寬度不大的清況下，L1，L2，L3，…的差別不大。但由于地球表面曲率的存在，對應(yīng)于地面的P1，P2，P3，…，顯然P3-P1＞L3-L1，像元對應(yīng)于地面寬度的不等距星下點越遠畸變越大，對應(yīng)地面長度越長。第三十三頁，共五十頁，編輯于2023年，星期五(4)大氣折射的影響大氣對輻射的傳播產(chǎn)生折射。由于大氣的密度分布從下向上越來越小，折射率不斷變化，因此折射后的輻射傳播不再是直線而是一條曲線，從而導(dǎo)致傳感器接收的像點發(fā)生位移大氣折射的影響NP第三十四頁，共五十頁，編輯于2023年，星期五(5)地球自轉(zhuǎn)的影響衛(wèi)星前進過程中，傳感器對地面掃描獲得圖像時，地球自轉(zhuǎn)影響較大，會產(chǎn)生影像偏離。因為多數(shù)衛(wèi)星在軌道運行的降段接收圖像，即衛(wèi)星自北向南運動，這時地球自西向東自轉(zhuǎn)。相對運動的結(jié)果，使衛(wèi)星的星下位置逐漸產(chǎn)生偏離。偏離方向如圖所示，所以衛(wèi)星圖像經(jīng)過校正后成為圖C的形態(tài)。地球自轉(zhuǎn)引起偏離(a)獲得圖像（b）實際對應(yīng)的地面位置（c）影像變形第三十五頁，共五十頁，編輯于2023年，星期五2.幾何畸變校正從具有幾何畸變的圖像中消除畸變的過程稱幾何校正。也可以說是定量地確定圖像上的像元坐標（圖像坐標）與目標物的地理坐標（地圖坐標等）的對應(yīng)關(guān)系（坐標變換式）。

第三十六頁，共五十頁，編輯于2023年，星期五

幾何畸變有多種校正方法，但常用的是一種通用的精校正方法，適合于在地面平坦，不需考慮高程信息，或地面起伏較大而無高程信息，以及傳感器的位置和姿態(tài)參數(shù)無法獲取的情況時應(yīng)用。有時根據(jù)遙感平臺的各種參數(shù)已做過一次校正，但仍不能滿足要求，就可以用該方法作遙感影像相對于地面坐標的配準校正，遙感影像相對于地圖投影坐標系統(tǒng)的配準校正，以及不同類型或不同時相的遙感影像之間的幾何配準和復(fù)合分析，以得到比較精確的結(jié)果。幾何校正需要利用地面控制點和多項式糾正模型。第三十七頁，共五十頁，編輯于2023年，星期五控制點選取的原則控制點的選擇要以配準對象為依據(jù)。以地面坐標為匹配標準的，叫做地面控制點（記作GCP）。有時也用地圖作地面控制點標準，或用遙感圖像（如用航空像片）作為控制點標準。無論用哪一種坐標系，關(guān)鍵在于建立待匹配的兩種坐標系的對應(yīng)點關(guān)系。一般來說，控制點應(yīng)選取圖像上易分辨且較精細的特征點，這很容易通過目視方法辨別，如道路交叉點、河流彎曲或分叉處、海岸線彎曲處、湖泊邊緣、飛機場、城廓邊緣等。圖像邊緣部分一定要選取控制點，以避免外推。此外，盡可能滿幅均勻選取，特征實在不明顯的大面積區(qū)域（如沙漠），可用求延長線交點的辦法來彌補，但應(yīng)盡可能避免這樣做，以避免造成人為的誤差。

地面控制點上的地物不隨時間而變化，以保證當(dāng)兩幅不同時段的圖像或地圖幾何糾正時，可以同時識別出來。第三十八頁，共五十頁，編輯于2023年，星期五（1）基本思路校正前的圖像看起來是由行列整齊的等間距像元點組成的，但實際上，由于某種幾何畸變，圖像中像元點間所對應(yīng)的地面距離并不相等（圖a）。校正后的圖像亦是由等間距的網(wǎng)格點組成的，且以地面為標準，符合某種投影的均勻分布（圖b），圖像中格網(wǎng)的交點可以看作是像元的中心。

校正的最終目的:A.確定校正后圖像的行列數(shù)值。B.找到新圖像中每一像元的亮度值。

第三十九頁，共五十頁，編輯于2023年，星期五（2）具體步驟:A.重采樣找到一種數(shù)學(xué)關(guān)系，建立變換前圖像坐標（x，y）與變換后圖像坐標（u，v）的關(guān)系，通過每一個變換后圖像像元的中心位置（u代表行數(shù)，v代表列數(shù)，均為整數(shù)）計算出變換前對應(yīng)的圖像坐標點（x，y）。分析得知，整數(shù)（u，v）的像元點在原圖像坐標系中一般不在整數(shù)（x，y）點上，即不在原圖像像元的中心。計算校正后圖像中的每一點所對應(yīng)原圖中的位置（x，y）。計算時按行逐點計算，每行結(jié)束后進入下一行計算，直到全圖結(jié)束。B.計算像元亮度值由于計算后的(x，y)多數(shù)不在原圖的像元中心處，因此必須重新計算新位置的亮度值。一股來說，新點的亮度值介于鄰點亮度值之間，所以常用內(nèi)插法計算。第四十頁，共五十頁，編輯于2023年，星期五（3）計算方法建立兩圖像像元點之間的對應(yīng)關(guān)系

（u,v)(x,y)，需求出12個系數(shù)，至少列出12個方程，即要找到6個已知的對應(yīng)點（控制點）如果要提高精度，必須大大增加控制點的數(shù)目，用最小二乘法進行曲面擬合求系數(shù)。實際計算時常采用2元2次多項式2元n次多項式第四十一頁，共五十頁，編輯于2023年，星期五為了確定校正后圖像上每點的亮度值，只要求出其原圖所對應(yīng)點（x，y）的亮度。通常有三種方法：最近鄰法雙向線性內(nèi)插法三次卷積內(nèi)插法。

第四十二頁，共五十頁，編輯于2023年，星期五I最近鄰法

圖像中兩相鄰點的距離為1，即行間距△x＝1，列間距△y=1，取與所計算點(x,y)周圍相鄰的4個點，比較它們與被計算點的距離，哪個點距離最近，就取哪個的亮度值作為（x,y）點的亮度值f（x，y）。設(shè)該最近鄰點的坐標為（k,l）,則k=Integer(x+0.5)l=Integer(y+0.5)f(x,y)=f(k,l)幾何位置上的精度為±0.5像元第四十三頁，共五十頁，編輯于2023年，星期五最鄰近內(nèi)插法以距內(nèi)插點最近的觀測點的像元值為所求的像元值。該方法最大可產(chǎn)生0.5個像元的位置誤差，優(yōu)點是不破壞原來的像元值，處理速度快。第四十四頁，共五十頁，編輯于2023年，星期五II雙線性內(nèi)插法

取（x，y）點周圍的4鄰點，在y方向（或x方向）內(nèi)插二次，再在x方向（或y方向）內(nèi)插一次，得到（x，y）點的亮度值f(x，y），該方法稱雙線性內(nèi)插法。設(shè)4個鄰點分別為(i,j),(i,j+1),(i+1,j),

(i+1,j+1),過(x,y)作直線與x軸平行，與4鄰點組成的邊相交于點（i,y)和(i+1,y)。先在y方向內(nèi)插,由f(i,j+1)和f(i,j)計算交點的亮度f（i,y)；由f(i+1,j+1)和f(i+1,j)計算交點的亮度f(i+1,y)。然后計算x方向，以f(i,y)和f(i+1,y)來內(nèi)插f(x,y)值。第四十五頁，共五十頁，編輯于2023年，星期五雙線性內(nèi)插法使用內(nèi)插點周圍的4個觀測點的像元值，對所求的像元值進行線性內(nèi)插。缺點是破壞了原來的數(shù)據(jù)，但具有平均化的濾波效果。第四十六頁，共五十頁，編輯于2023年，星期五III三次卷積內(nèi)插法增加鄰點來獲得最佳插值函數(shù)。取與計算點（x，y）周圍相鄰的16個點，與雙向線性內(nèi)插類似，可先在某一方向上內(nèi)插，每4個值依次內(nèi)插4次，求出f（x，j-1），f(x,j)，f(x,j+1),f(x,j+2)，再根據(jù)這四個計算結(jié)果在另一方向上內(nèi)插，得到f(x，y）。因這種三次多項式內(nèi)插過程實際上是一種