第四章數(shù)字圖像的校正

1、數(shù)字圖像2、輻射校正3、幾何校正4.2數(shù)字圖像校正1模擬圖像與數(shù)字圖像遙感圖像的表示：光學圖像（模擬圖像）——連續(xù)圖像——黑白照片數(shù)字圖像——離散圖像——計算機存儲模擬圖像與數(shù)字圖像之間的轉換，稱為模/數(shù)（A/D）轉換。2一、數(shù)字圖像指能夠被計算機存儲、處理和使用的圖像“離散化”、二維矩陣：每個元素的取值是圖像連續(xù)變化的灰度的離散整數(shù)值31pixel將一幅二維的圖像通過有限個離散點來表示就成為了數(shù)字圖像，其中的每個點稱為像元。4照片照片上的手指QuickBird（校園）局部放大5模擬方法表示圖像：O圖像xy(x,y)表示像元的位置；f(x,y)表示（x,y）位置上的對應地物電磁輻射強度。對于模擬圖像,x,y,f(x,y)的取值是連續(xù)的。6O圖像xy1234213456Pixel

函數(shù)f(x,y)的取值：離散整數(shù)取樣是根據(jù)需要，將灰度空間分成2n級（目前n的取值有1、4、7、8，甚至更多），然后根據(jù)方格內電磁輻射強弱取其平均值整數(shù)作為函數(shù)f(x,y)的值。對圖像進行“離散化”：7一幅模擬圖像表示為數(shù)字圖像其實質是一個數(shù)字矩陣。數(shù)字矩陣可以在計算機里進行存儲和運算。8圖像的數(shù)字化內容：（1）圖像空間位置的數(shù)字化，即圖像的空間取樣。（2）圖像灰度的數(shù)字化，即指從圖像灰度的連續(xù)變化中進行離散的采樣，目前經常使用的灰度量度有2級，64級，128級，256級。除光學圖像可以數(shù)字化為數(shù)字圖像外，更多的遙感圖像源于傳感器獲得后直接的數(shù)字產品，如MSS,TM,ETM等航天遙感器。9遙感圖像處理系統(tǒng)計算機（主機）圖像輸入輸出設備磁帶機、數(shù)字化器等；打印機、繪圖儀、激光圖像記錄儀專用處理設備圖像計算機、陣列處理機外存設備磁盤、磁帶、光盤顯示器軟件部分系統(tǒng)軟件、應用軟件（圖像處理軟件）10

圖像處理系統(tǒng)11二、輻射校正由于傳感器響應特性和大氣的吸收、散射及其它隨機因素影響，導致圖像模糊失真，造成圖像分辨率和對比度相對下降，稱為輻射畸變。這些畸變都需要通過輻射校正復原。引起輻射畸變原因：傳感器本身產生的誤差（生產單位進行校正）；大氣對輻射的影響（用戶自行校正）。12

遙感傳感器感測的信息，是地物對太陽光的反射或地物發(fā)射的電磁波經過大氣層傳輸并與大氣發(fā)生作用后的結果。由于大氣的散射作用，傳感器在接收地物輻射信息的同時，也接收了散射所造成的非地物輻射能，從而使得遙感圖像對比度下降，導致圖像猶如蒙上了一層薄紗一樣不清晰。

大氣校正大氣校正的原因：131415（1）大氣影響的定量分析地面上單位面積的輻照度為（無大氣）

假定地表面是朗伯體，其表面為漫反射，則某方向物體的亮度為

儀器系統(tǒng)增益系數(shù)因子，進入傳感器的亮度值

無大氣影響狀況j16(1)大氣影響的定量分析受入射方向透過率和反射方向透過率的影響，進入傳感器的亮度值為：漫入射輻照度為ED，反射后進入傳感器的亮度值為：散射光向上直接進入傳感器的輻射亮度為LP（程輻射度）①②③總輻射亮度:17(1)大氣影響的定量分析大氣的影響主要是減少了圖像的對比度，使原始信號和背景信號都增加了因子。18(2)大氣影響的粗略校正精確校正：找出每個波段像元亮度值與地物反射率的關系。粗略校正：通過比較簡便的方法去掉式中的LP，從而改善圖像質量。

輻射校正：對大氣影響的糾正是通過糾正輻射亮度的辦法實現(xiàn)的。19（2）大氣影響的粗略校正直方圖最小值去除法回歸分析法①②20灰度直方圖

用平面直角坐標系表示一幅灰度范圍為0-n的數(shù)字圖像像元灰度分布狀態(tài)。橫軸表示灰度級，縱軸（Pi=mi/M)表示灰度級為gi的像元個數(shù)mi占像元總數(shù)Ｍ的百分比。將2n個Pi繪于圖上，所形成的統(tǒng)計直方圖叫灰度直方圖。21通過灰度直方圖可以直觀地了解圖像的質量。22不同反差特征的圖像23直方圖校正方法前提（假設）：深海水體處（或山的陰影）等物體的亮度值為0，大氣散射導致圖像上這些物體的灰度值不為0（輻射偏置量）。從圖像像元亮度值中減去一個輻射偏置量（LP），輻射偏置量等于圖像直方圖中最小的輻射亮度值。暗物體法（Dark-objectmethod）①24反差增強25調整前直方圖調整后直方圖26回歸分析法原理：在遙感圖像上大山的陰影區(qū)或深大水體區(qū)域，各個波段的反射為零。同時，大氣散射主要影響短波部分，波長較長的波段幾乎不受影響，因此可用其校正其它波段數(shù)據(jù)。方法：在不受大氣影響的波段（如TM5或7）和待校正的某一波段圖像中，選擇一系列目標，將每個目標的兩個待比較的波段灰度值提取出來進行回歸分析，建立線性回歸方程。Lb：待校正波段的圖像亮度值La：不受大氣影響波段的圖像亮度值②27回歸分析法怎么求a,b？②28a,b計算如下：可以認為，a就是波段b程輻射值，將波段b中每個像元的亮度值減去a，來改善圖像，去掉程輻射值。29由于遙感傳感器、遙感平臺以及地球自轉等因素的影響，遙感成像時往往會引起難以避免的幾何畸變。幾何校正就是要校正成像過程中所造成的各種幾何畸變。

幾何校正可分為兩種：幾何粗校正和幾何精校正。

三、幾何校正幾何校正的原因：30

幾何變形：遙感圖像的幾何位置發(fā)生變化，產生諸如行列不均勻，像元大小與地面對應不準確，地物形狀不規(guī)則變化等。31遙感影像變形的原因（1）遙感平臺位置和運動狀態(tài)變化的影響32遙感影像變形的原因（2）地形起伏影響33遙感影像變形的原因（3）地球表面曲率的影響

像點位置的移動像元對應于地面寬度的不等34遙感影像變形的原因（3）地球表面曲率的影響

35遙感影像變形的原因（4）大氣折射的影響

36遙感影像變形的原因（5）地球自轉的影響

37幾何粗校正是針對引起畸變的原因而進行的校正。進行校正時，只需將傳感器的校準數(shù)據(jù)、遙感平臺的位置以及衛(wèi)星運行姿態(tài)等一系列測量數(shù)據(jù)代入理論校正公式即可。

幾何精校正是利用控制點進行的幾何校正。它是用一種數(shù)學模型來近似描述遙感圖像的幾何畸變過程，求得這個幾何畸變模型，然后利用此模型對遙感圖像進行幾何校正，無需考慮引起畸變的具體原因。幾何校正方法：38幾何畸變校正基本思路39幾何畸變校正

通過計算機對圖象每個象元逐個地解析糾正處理完成，可以較精確改正線性和非線性變形誤差。

包括兩個方面（基本環(huán)節(jié)）：

（1）像元坐標變換；（2）像元灰度值重新計算（重采樣）。40數(shù)字圖像幾何糾正的主要處理過程準備工作輸入原始數(shù)字影像建立糾正變換函數(shù)確定輸出影像范圍像元坐標變換像元亮度值重采樣輸出糾正后的圖像準備工作：圖像、地圖、大地測量資料、平臺軌道參數(shù)、傳感器參數(shù)、控制點的選擇。（具體內容可選）糾正變換函數(shù)建立：輸入和輸出圖像間的坐標變換關系：如多項式法、共線方程法等。41☆幾何精校正技術關鍵：地面控制點（GCP）的選取、地面控制點的數(shù)目及校正方程的次數(shù)。☆校正步驟：Ａ、調入圖像，比較觀察；Ｂ、確定二元多項式的次數(shù)Ｃ、確定GCP的個數(shù)；Ｄ、選取GCP，盡可能均勻分布在圖像上；Ｅ、再采樣計算，得到校正后的新圖像。具體步驟42數(shù)字圖像幾何校正方法(函數(shù)的建立)多項式糾正法：常用的方法基本原理：不考慮成像的空間幾何過程，而直接對圖象變形的本身進行數(shù)學模擬。把遙感圖象的總體變形看作是平移、縮放、旋轉、仿射、偏扭、彎曲及更高次的基本變形的綜合作用結果，因此，將糾正前后圖像相應點間的坐標關系用一適當?shù)亩囗検奖磉_一般取2次，可以滿足精度要求需要地面控制點：分布、數(shù)量43確定糾正后圖像的邊界范圍

糾正后圖像的邊界范圍，指在計算機存儲器中為輸出影像所開出的存儲空間大小，以及該空間邊界范圍（首行、首列、末行、末列）的地圖（或地面）坐標定義值；糾正后圖像和原始圖像的形狀、大小、方向都不一樣。所以在糾正過程的實施之前，必須首先確定新圖像的大小范圍。44

先求出原始圖像四個角點（a,b,c,d）在糾正后圖像中的對應點（a’,b’,c’,d’）的坐標（Xa’,Ya’）(Xb’,Yb’)(Xc’,Yc’)(Xd’,Yd’)；然后求出最大值和最小值。

為了將該邊界范圍轉換為計算機糾正后的存儲數(shù)組空間，須在其中劃分出網格，每個網格代表一個輸出像元；須根據(jù)精度要求確定輸出像元的地面尺寸；確定糾正后圖像的邊界范圍45確定糾正后圖像的邊界范圍X1=min(Xa’,Xb’,Xc’,Xd’)X2=max(Xa’,Xb’,Xc’,Xd’)Y1=min(Ya’,Yb’,Yc’,Yd’)Y2=max(Ya’,Yb’,Yc’,Yd’)46坐標變換確定原始圖像和糾正后圖像間的坐標變換關系（x,y）---(X,Y)(X行數(shù),Y列數(shù),均為整數(shù))有兩種方案：直接糾正方案和間接糾正方案直接糾正方案：從原始圖像陣列出發(fā)，依次對其中每一個像元P(x，y)分別計算其在輸出（糾正后）圖像的坐標P(X，Y)，并計算P(X，Y)的灰度值;間接糾正方案：從空白圖像陣列出發(fā)，依次計算每個像元P(X，Y)在原始圖象中的位置P(x，y)，然后把該點的灰度值依次計算后返送給P(X，Y).

二者間并無本質差別，互為逆變換47坐標關系

（xp,yp）（XP,YP）分別是任意一個像元在原始圖像和糾正后圖像中的坐標。直接間接48直接糾正方法與間接糾正法:二者并無本質區(qū)別，互為逆變換坐標變換49圖像灰度值的重采樣幾何校正過程中，由于校正前后圖像的像元大小可能變化、像元點位置的相對變化等，不能簡單用原圖像像元灰度值代替輸出像元灰度值。因此，需要插值（重采樣）

三種插值方法：最鄰近法、雙線性法、三次卷積法50灰度值重采樣示意圖51最鄰近法（nearestneighbor）用距離投影點(采樣點)最近像元灰度值代替輸出像元灰度值。簡單、省時，保留原始圖象的值處理后圖像的亮度具有不連續(xù)性，邊緣出現(xiàn)鋸齒狀52該最近鄰點的坐標為：K=integer（x+0.5）L=integer（y+0.5）53原始圖像糾正后圖像（最鄰近插值）54雙線性內插法取（x，y）點周圍的4個鄰點，在y方向內插兩次，再在x方向內插一次，得到（x，y）點的亮度值f(x，y)，該方法稱雙線性內插法。55雙線性內插法比起最近鄰法雖然計算量增加，但精度明顯提高，特別是對亮度不連續(xù)現(xiàn)象或線性特征的塊狀化現(xiàn)象有明顯的改善。但是這種內插法會對圖像起到平滑作用，從而使對比度明顯的分界線變得模糊。雙線性內插法56雙線性插值原始圖像糾正（雙線性插值）57三次卷積法(cubicconvolution)取與投影點（x，y）鄰近的16個像元灰度值，先在某一方向上內插，每四個值依次內插四次，得出四個結果，再根據(jù)這四個結果在另一方向計算輸出像元的灰度值。計算量大，但圖像質量更好，細節(jié)表現(xiàn)更為清楚。需注意的是，如果想用此方法獲得好的圖像效果，就要求位置校正過程更精確。5859原始圖像幾何糾正（三次卷積）60控制點的選擇地面控制點（GCP，GroundControlPoint）:一些特定的像元，其地圖坐標或其它輸出坐標為已知人工地物線性地物交叉點不易隨時間變化的目標大比例尺的圖像：道路交叉點、機場跑道