遙感圖像幾何處理

1、主要內(nèi)容：遙感傳感器的構(gòu)像方程遙感圖像幾何變形遙感圖像的幾何處理遙感圖像幾何處理的應(yīng)用遙感圖像幾何處理一、遙感傳感器的構(gòu)像方程構(gòu)像方程構(gòu)像方程n指地物點的圖像坐標(biāo)（x，y）和地面坐標(biāo)（X，Y，Z）之間的數(shù)學(xué)關(guān)系。n是遙感圖像幾何數(shù)字糾正的數(shù)學(xué)基礎(chǔ)n是參量誤差分析的基礎(chǔ)構(gòu)像方程中的三個坐標(biāo)系：構(gòu)像方程中的三個坐標(biāo)系： 地面坐標(biāo)系O-XYZ 圖像坐標(biāo)系o-xy 傳感器坐標(biāo)系S-UVW通用構(gòu)像方程n是傳感器坐標(biāo)系S-UVW與地面坐標(biāo)系O-XYZ之間的坐標(biāo)轉(zhuǎn)換關(guān)系n可據(jù)通用構(gòu)像方程導(dǎo)出各種不同傳感器的構(gòu)像方程PPPSSSpppWVUAZYXZYX1、框幅攝影機的構(gòu)像方程2、全景攝影機的構(gòu)像方程3、推掃

2、式傳感器的構(gòu)像方程4、紅外和多光譜掃描儀的構(gòu)像方程5、側(cè)視雷達(dá)圖像的構(gòu)像方程二、遙感圖像的幾何變形二、遙感圖像的幾何變形n遙感圖像的幾何變形是指原始圖像上各地物的幾何位置、形狀、尺寸、方位等特征與在參照系統(tǒng)參照系統(tǒng)中的表達(dá)要求不一致時產(chǎn)生的變形。n圖像投影的參照系統(tǒng)地圖投影系統(tǒng)n近似地圖投影：地球切平面坐標(biāo)系 原點：傳感器的星下點 X軸：東方向 Y軸：北方向 Z軸：由星下點鉛垂向上 引起變形誤差的原因傳感器成像方式傳感器外方位元素的變化地形起伏地球曲率大氣折射地球自轉(zhuǎn)三、遙感圖像幾何處理三、遙感圖像幾何處理n改正遙感圖像中的幾何變形，并將其投影到需要的地理坐標(biāo)系中。 目的是為滿足應(yīng)用中量測和定

3、位的要求、實現(xiàn)多源遙感圖像的幾何配準(zhǔn)、滿足利用遙感圖像進(jìn)行地形圖測圖或更新的要求、等。n糾正方法： 光學(xué)方法l使用現(xiàn)有的航測儀器l對動態(tài)獲取的影像只能進(jìn)行近似糾正（線性變形）l具體方法：分塊糾正、仿射糾正l正射影像：可消除地形起伏產(chǎn)生的像點位移數(shù)字方法l 使用計算機 ，處理數(shù)字圖像 l 建立在嚴(yán)格的數(shù)學(xué)基礎(chǔ)上，可以逐點（逐像素）地對圖像進(jìn)行糾正，原則上可以處理任何類型的傳感器圖像。數(shù)字幾何處理主要內(nèi)容：n數(shù)學(xué)基礎(chǔ)n幾何處理的兩個層次n數(shù)字糾正的基本原理n常用的糾正變換方法（一）數(shù)學(xué)基礎(chǔ)（一）數(shù)學(xué)基礎(chǔ)：遙感制圖：幾何處理：圖像坐標(biāo)系（出發(fā)點）地圖坐標(biāo)系（歸宿）地面坐標(biāo)地圖坐標(biāo)（地圖投影）地面坐標(biāo)

4、圖像坐標(biāo)（構(gòu)像方程）因此，幾何處理的實質(zhì)是將由構(gòu)像方程建立的關(guān)系與由地圖投影建立的關(guān)系相統(tǒng)一，進(jìn)一步滿足制圖的幾何要求。（二）幾何處理的兩個層次（二）幾何處理的兩個層次n粗加工處理 僅做系統(tǒng)誤差改正需要利用圖像的構(gòu)象方程來進(jìn)行 對傳感器的內(nèi)部畸變改正很有效 處理后仍有較大殘差n精加工處理 進(jìn)一步消除圖像中的幾何變形，產(chǎn)生符合某種地圖投影或圖形表達(dá)要求的新圖像； （三）數(shù)字糾正的基本原理（三）數(shù)字糾正的基本原理準(zhǔn)備工作輸入原始數(shù)字圖像建立糾正變換函數(shù)確定輸出圖像范圍逐像素的幾何位置變換像素亮度值重采樣輸出糾正后的圖像1.坐標(biāo)糾正變換兩種方案坐標(biāo)糾正變換兩種方案Q 直接法（需進(jìn)行像元的重新排列，要

5、求存儲空間大一倍，計算時間也長）Q 間接法（常采用）2.亮度值重采樣亮度值重采樣n幾何變換后坐標(biāo)計算值不為整數(shù)時，需要把該點周圍鄰近整數(shù)點位上的亮度值對該點的亮度貢獻(xiàn)累積起來，構(gòu)成該點位的新亮度值，即亮度值重采樣。n理想的重采樣函數(shù)是辛克（SINC）函數(shù)。n由于辛克函數(shù)使用不方便，所以常采用一些近似函數(shù)來代替它。三種常用的重采樣方法三種常用的重采樣方法1、雙三次卷積重采樣法（內(nèi)插精度較高，但計算量大）2、雙線性內(nèi)插法（實踐中常采用）3、最鄰近像元采樣（簡單計算量小、輻射保真度好，但幾何精度低）雙線性內(nèi)插法212221121121yyxxTyxpWWIIIIWWWIWIxWxWxx211yWyW

6、yy211) 1|0(|,|1)(cccxxxW3.幾何校正的類型n用地圖校正影像n用影像校正影像n混合校正方法：即用經(jīng)過幾何校正的影像來校正待糾正的影像（四）各種糾正變換方法n多項式糾正法n嚴(yán)格的物理模型共線方程糾正法n基于仿射變換的嚴(yán)格幾何模型n有理函數(shù)模型n神經(jīng)網(wǎng)絡(luò)校正1.多項式糾正法n回避成像的空間幾何過程，直接對圖像變形的本身進(jìn)行數(shù)學(xué)模擬，適用于各種類型傳感器影像的糾正。n用一個適當(dāng)?shù)亩囗検絹肀磉_(dá)糾正前后圖像相應(yīng)點之間的坐標(biāo)關(guān)系n校正精度與地面控制點的精度、分布、數(shù)量及實際地形有關(guān)。對地面相對平坦的情況，具有足夠好的糾正精度；對于地形起伏較大的地區(qū)，效果不好，特別是當(dāng)傾斜角大于10度

7、時，效果更差。n多項式模型糾正時，在控制點上擬合很好，但在其它點上可能有明顯偏離。n常用的多項式有一般多項式、勒讓德多項式以及雙變量分區(qū)插值多項式。一般多項式糾正變換公式)()()()(2542321025423210YbXYbXbYbXbbyYaXYaXaYaXaax多項式中系數(shù)的個數(shù)N與階數(shù)n的關(guān)系： 2/ )2)(1(nnN多項式系數(shù)的確定方法： 1、用可預(yù)測的圖像變形參數(shù)構(gòu)成； 2、利用已知控制點的坐標(biāo)值按最小二乘法原理求解。地面控制點要求： 在影像上為明顯的地物點，易于判讀； 在影像上均勻分布。 改進(jìn)的多項式模型minjpkkjiijkminjpkkjiijkZYXbyZYXax00

8、0000小結(jié)：n多項式法最簡單，用來糾正控制點處的平面變形，不考慮地面高差，適用于小范圍平坦地區(qū)；n2D多項式模型不能反映各種誤差來源，也不考慮地形高差引起的畸變，只限于小范圍、高差較小的地區(qū)，變形較小的如星下點獲取的或經(jīng)過幾何校正的數(shù)據(jù)；n2D多項式模型不適于遙感影像特別是高分辨率遙感影像的精校正；n3D多項式模型只有在無法獲取成像系統(tǒng)的相關(guān)參數(shù)或沒有可用的參數(shù)模型的條件下使用，3D非參數(shù)模型不需要成像系統(tǒng)的先驗知識，沒有明確的物理意義，通常是基于不同的XYZ數(shù)學(xué)模型。3D多項式模型是2D多項式模型的擴(kuò)展，適用于小范圍的影像數(shù)據(jù)，需要大量的高精度、分布良好的GCP,對輸入誤差非常敏感2.共線

9、方程糾正法（嚴(yán)格物理模型）n在范圍較大的的情況下，地球切平面投影與地圖投影之間有著較大差異，因此需要通過更嚴(yán)密的變換來建立地物的圖像坐標(biāo)與地圖坐標(biāo)之間的關(guān)系n以地心坐標(biāo)系為基礎(chǔ)的共線方程以地心坐標(biāo)系為基礎(chǔ)的共線方程：既能與傳感器坐標(biāo)系直接進(jìn)行三維空間線性變換，恢復(fù)成像光束的空間幾何形態(tài)，又能借助于大地測量學(xué)和地圖投影學(xué)的知識，方便地轉(zhuǎn)換為地理經(jīng)緯度坐標(biāo)，進(jìn)而轉(zhuǎn)換為任意一種所需要的地圖投影坐標(biāo)。(1)、以地心直角坐標(biāo)為基礎(chǔ)的共線方程ZYXWVUCBAZYXDZYXZYXDZYXZYXPPPOOOSCCCPSCCCPCCC切平面坐標(biāo)系切平面坐標(biāo)系地心坐標(biāo)系地心坐標(biāo)系令傳感器S的地心坐標(biāo)為：OOOS

10、CCCSCCCZYXDZYXZYXADM令傳感器平臺坐標(biāo)系與地心坐標(biāo)系間的旋轉(zhuǎn)矩陣為：ZYXWVUCBZYXPPPPPP地物點P在遙感平臺坐標(biāo)系中的坐標(biāo)為：PPPSCCCPCCCZYXMZYXZYX)()()()()()(: )()()()()()()(: )(332313322212332313312111CSCPCSCPCSCPCSCPCSCPCSCPppCSCPCSCPCSCPCSCPCSCPCSCPppZZmYYmXXmZZmYYmXXmfZYfyZZmYYmXXmZZmYYmXXmfZXfx傳感器平臺坐標(biāo)系與地心坐標(biāo)系間的坐標(biāo)變換表達(dá)式為：基于地心坐標(biāo)的共線方程：地心坐標(biāo)系地心坐標(biāo)系

11、地理經(jīng)緯度坐標(biāo)系地理經(jīng)緯度坐標(biāo)系 根據(jù)大地測量學(xué)的知識，空間地物點的地心坐標(biāo)、地心經(jīng)緯度、航偏角和地球曲率半徑都可以轉(zhuǎn)換為地理經(jīng)緯度的表達(dá)形式。地理經(jīng)緯度坐標(biāo)系地理經(jīng)緯度坐標(biāo)系地圖投影坐標(biāo)系地圖投影坐標(biāo)系坐標(biāo)正算:地面點經(jīng)緯度（L，B）地圖坐標(biāo)(Xm,Ym)坐標(biāo)反算:地圖坐標(biāo)(Xm,Ym) 地面點經(jīng)緯度（L，B）(2)、共線方程參數(shù)的確定、共線方程參數(shù)的確定n共線方程中的待定系數(shù)主要是：傳感器的三個姿態(tài)參數(shù)和三個位置參數(shù)。n參數(shù)選擇可按照兩種方案：n參數(shù)解算的兩種方式 1、利用可預(yù)測參數(shù)來直接構(gòu)成 2、借助于控制點按最小二乘法原理來求解。 mmmcscscsZYX,sssHBL(3)、共線方程

12、參數(shù)隨時間變化的表征函數(shù)、共線方程參數(shù)隨時間變化的表征函數(shù)n對動態(tài)傳感器，每個像素（或每條掃描線）都有各自的一套共線方程參數(shù)。n通常可以把整幅圖像成像過程中共線方程參數(shù)的變化看作時間的連續(xù)函數(shù)（共線方程參數(shù)的表征函數(shù)），用來表達(dá)任一時刻傳感器的位置和姿態(tài)。(4)、地物點高程數(shù)據(jù)的準(zhǔn)備、地物點高程數(shù)據(jù)的準(zhǔn)備n共線方程中需要每個像素的地面高程n建立數(shù)字地面模型可采用等高線數(shù)字化等方法；對分辨率要求不高時，也可在相關(guān)網(wǎng)站上下載；n糾正時，輸出圖像應(yīng)與DTM進(jìn)行幾何配準(zhǔn)DEM數(shù)據(jù)下載小結(jié)：n建立在嚴(yán)密的數(shù)學(xué)基礎(chǔ)上，是對成像空間幾何形態(tài)的直接描述。在糾正過程中引入了地面高程信息，比多項式法的應(yīng)用范圍更廣

13、。適用于各種分辨率的遙感影像糾正?，F(xiàn)在商業(yè)軟件中基本都以此為基礎(chǔ)實現(xiàn)各種來源的遙感影像校正功能。n雖有嚴(yán)密的理論基礎(chǔ)，但數(shù)學(xué)模型中參數(shù)的確定有著很強的近似性，因此其精度并不比多項式糾正的精度高。目前對該方法的研究主要是參數(shù)的解算問題。n 嚴(yán)格物理模型基于衛(wèi)星軌道、攝影測量、測地學(xué)和地圖學(xué)，模型反映了影像獲取時的幾何物理狀態(tài)，用來糾正由于衛(wèi)星、傳感器、地球和地圖投影等引起的變形。嚴(yán)格物理模型需要的控制點少，一般情況下3-6個，精度更高，對于高分辨率衛(wèi)星數(shù)據(jù)可以獲得一兩個像素的精度。n當(dāng)僅僅處理原始影像的一部分時，當(dāng)影像已經(jīng)經(jīng)過幾何校正時，或者你無法獲取衛(wèi)星軌道信息時，都不能采用嚴(yán)格物理模型。3.

14、基于仿射變換的嚴(yán)格幾何模型n用于高分辨率遙感影像的糾正方法n高分辨率衛(wèi)星傳感器的突出特征是長焦距和窄視場角。這種成像幾何關(guān)系若用共線方程描述會導(dǎo)致定向參數(shù)間很強的相關(guān)性。假設(shè)小視場角內(nèi)的中心投影近似于平行光投影，利用仿射模型求解方位參數(shù)，可克服方位參數(shù)的相關(guān)性。n對10米分辨率的SPOT影像 用于較小比例尺制圖、精度要求較低的情況下非常有效。對更高分辨率的影像，這種近似方法的有效性還需進(jìn)一步研究。二維仿射變換成像模型87654321/tan1)cos/()(1AZAYAXAyfyZZZAZAYAXAx4.有理函數(shù)模型n各種傳感器幾何模型的一種更廣義的表達(dá)形式，是對不同的傳感器模型更為精確的表達(dá)

15、形式，能適用于各類傳感器。n優(yōu)點是由于引入較多的定向參數(shù)，模擬精度很高；缺點是模型解算復(fù)雜，運算量大，并且要求控制點的數(shù)目相對較多。n是高分辨率遙感影像糾正的重要方法。 RPC模型的實質(zhì)： 是有理函數(shù)糾正模型，是將地面點大地大地坐標(biāo)坐標(biāo)與其對應(yīng)的像點坐標(biāo)像點坐標(biāo)用比值多項式關(guān)聯(lián)起來。為了增強參數(shù)求解的穩(wěn)定性，將地面坐標(biāo)和影像坐標(biāo)正則化到-1和1之間。對于一幅影像，定義如下比值多項式：),(),(),(),(4321nnnnnnnnnnnnnnZYXpZYXpcZYXpZYXpr123000mimjmkkjiijkZYXaP319318317216215214213212211102928276

16、543210XaYaZaYXaZXaXYaZYaXZaYZaZYXaXaYaZaXYaZXaZYaXaYaZaa（多項式次數(shù)最高不超過三次）在有理函數(shù)模型中n光學(xué)投影系統(tǒng)產(chǎn)生的誤差用一次項表示n地球曲率、大氣折射和鏡頭畸變等產(chǎn)生的誤差用二次項來模型化n未知的具有高階分量的誤差如相機振動等，用三次項來表示n在使用有理函數(shù)模型校正時，可以選擇使用或不使用GCP。n有理函數(shù)模型最最是為沒有GCP,只有DEM來對衛(wèi)星影像進(jìn)行幾何校正的用戶設(shè)計的，對于GCP無法獲取的偏遠(yuǎn)地區(qū)是非常有用的。n大多數(shù)的商業(yè)影像處理軟件和攝影測量軟件都支持有理函數(shù)模型，可以用來處理單片和立體像對，如ERDAS、PCI、ENV

17、I等，可以利用RPC（有理多項式系數(shù)）模型進(jìn)行正射糾正、3D特征提取、DEM生成和區(qū)域平差；n有理函數(shù)模型（RFM ）的應(yīng)用目前主要是在RPC 系數(shù)已知的情況下（與地形無關(guān)）。ERDAS 9. 1 軟件中已經(jīng)集成了目前主要的高分辨率傳感器模型RPC 參數(shù)解算方法。但是,因為RFM需要的控制點數(shù)目相對較多,而且在解算時對控制點的分布要求均勻分布,否則會導(dǎo)致法方程矩陣奇異,迭代求解可能不收斂,因而目前的商業(yè)軟件中也大都沒有該模型的地形相關(guān)方案的模塊集成。nRFM也有其局限性：處理的影像大小受限，為了達(dá)到較高的精度，需要大量的規(guī)則分布的GCP;只能校正控制點處的誤差，控制點之間的誤差并不能消除。幾種幾何校正主要方法的比較在相同的控制點和已有數(shù)據(jù)條件下：n共線方程模型不一定總能得到較高的糾正精度，一個重要的原因是傳感器的定向參數(shù)的確定精度；