4-地物光譜數(shù)據(jù)的重建PPT課件

1、2021/7/231第四章 地物光譜數(shù)據(jù)重建2021/7/232定義定義由于受到太陽位置、角度條件、大氣條件、地形影響及傳感器本身的性能的影響，傳感器所記錄的數(shù)據(jù)與目標的光譜反射率或光譜輻射亮度值并不一致。因此，將傳感器記錄的原始輻射值（DN值）轉(zhuǎn)化為地物反射率，恢復地物光譜數(shù)據(jù)本來的面目，稱為地物光譜數(shù)據(jù)重建。2021/7/233遙感器定標大氣校正 光照及地形校正 遙感器入瞳輻射值 地表輻射值 地表真實反射值 原始數(shù)字信號地物光譜數(shù)據(jù)重建的過程：2021/7/234l輻射誤差：傳感器所得到的目標測量值與目標的光譜反射率或光譜輻亮度等物理量之間的差值l輻射誤差產(chǎn)生原因：傳感器響應特性和外界自然

2、環(huán)境 地物（目標物）的輻射（反射）經(jīng)過大氣層時，地物（目標物）的輻射（反射）經(jīng)過大氣層時，與大氣層發(fā)生散射作用和吸收作用。吸收作用直接與大氣層發(fā)生散射作用和吸收作用。吸收作用直接降低地物的輻射能量，引起輻射畸變。散射作用除降低地物的輻射能量，引起輻射畸變。散射作用除降低地物的輻射能量外，大氣散射的部分輻射還會降低地物的輻射能量外，大氣散射的部分輻射還會進入傳感器，直接疊加在目標地物的輻射能量之中，進入傳感器，直接疊加在目標地物的輻射能量之中，成為目標地物的噪聲，降低了圖像的質(zhì)量。成為目標地物的噪聲，降低了圖像的質(zhì)量。 光照條件的不同也會引起輻射畸變，如太陽高光照條件的不同也會引起輻射畸變，如太

3、陽高度角、地面坡度等，都會引起輻射的畸變。度角、地面坡度等，都會引起輻射的畸變。輻射誤差2021/7/2354.1 4.1 遙感器定標遙感器定標l一般把消除傳感器本身影響的輻射校正過程稱之為輻射定標，即遙感器定標。l對高光譜定量遙感具有十分重要的意義，只有經(jīng)過精確的定標后，才能獲取真實的地物物理參量，才能與不同的高光譜遙感數(shù)據(jù)，模擬數(shù)據(jù)進行對比分析。l具體形式：建立成像光譜儀每個探測元件輸出的數(shù)字量化值（DN）與它對應視場中輸出輻射亮度值之間的定量關(guān)系。2021/7/236l遙感器輻射定標是指建立遙感器每個探測元所輸出信號的數(shù)值量化值與該探測器對應像元內(nèi)的實際地物輻射亮度值之間的定量關(guān)系。lL

4、i波段i的入瞳輻射能量lDNi波段i傳感器的輸出亮度值lAi波段i的定標增益系數(shù)lBi 波段i的定標偏置量BiDNiAiLi2021/7/237遙感器定標主要包括三個階段：實驗室定標和星上定標，場地定標（1 1）實驗室定標：）實驗室定標：是成像光譜儀在運行前在實驗室所接受的定標過程。包括光譜定標和輻射定標光譜定標：確定成像光譜儀每個波段的中心波長和帶寬。2021/7/238輻射定標分為兩種形式：相對定標和絕對定標。l相對定標：校正遙感器中各個探測器元件響應度差異而對衛(wèi)星遙感器測量到的原始數(shù)字進行歸一化的工作。l絕對定標：通過各種標準輻射源，在不同波譜段建立成像光譜儀入瞳處光譜輻射亮度值與成像儀

5、輸出的數(shù)字值（DN）之間的定量關(guān)系。2021/7/239其他系統(tǒng)誤差校正其他系統(tǒng)誤差校正l常見的由遙感系統(tǒng)引入的輻射誤差有： 隨機壞像元（散粒噪聲） 行起始/終止問題 行或列部分缺失 行或列條紋2021/7/2310隨機壞像元(散粒噪聲)BV1BV2BV3BV8BVijkBV4BV7BV6BV5原因：一個探測器未記錄某個光譜數(shù)據(jù)，該像元成為壞像元算法步驟：計算相鄰8個像元的均值，用均值代替壞像元的灰度值。2021/7/2311行或列缺失l原因：掃描系統(tǒng)的某個探測器工作不正常，就有可能產(chǎn)生一整行沒有光譜信息的線。對于線陣列掃描系統(tǒng)而言，就會使得整行或者整列數(shù)據(jù)缺失。l解決方法：人工處理。2021

6、/7/2312（2 2）機上或星上定標：）機上或星上定標：l飛行之前的實驗室定標會隨著環(huán)境的變化而變化，由于受到發(fā)射的影響和空間工作的溫度變化，發(fā)射前的定標關(guān)系發(fā)生變化，因此，實驗室定標后飛機或衛(wèi)星上的定標也是必要的。l機上定標的內(nèi)容和形式基本上與實驗室定標是類似的，其中穩(wěn)定的光源是主要因素。2021/7/2313（3 3）地面場定標）地面場定標l地面場定標：在遙感器處于正常運行條件下，選擇輻射的定標場地，通過地面同步測量對遙感進行定標。對前面兩項定標進行檢驗和修正。2021/7/2314l場地定標的原理：遙感器升空后，在定標場地選擇若干像元區(qū)，測量成像光譜儀對應地物的各波段光譜反射率和大氣光

7、譜參量，利用大氣輻射傳輸模型給出成像光譜儀入瞳處各譜帶的輻射亮度，最后確定它與成像光譜儀對應的輸出的數(shù)字量化值的關(guān)系，求解定標系數(shù)，并進行誤差分析。2021/7/23154.2 4.2 大氣校正大氣校正大氣輻射傳輸理論：l大氣散射集中在太陽輻射能量最強的可見光區(qū)。因此，散射是太陽輻射衰減的主要原因。l大氣對電磁輻射產(chǎn)生的散射和吸收是主要的輻射誤差源，因此，遙感圖像的大氣校正主要是補償大氣在吸收與散射方面的瞬時影響。l大氣校正方法可以大致分為：（1）基于圖像統(tǒng)計模型；（2）專用的大氣校正模型2021/7/2316這是一種相對定標方法，已廣泛地應用在成像光譜遙感數(shù)據(jù)處理中。模型要求在處理的影像數(shù)據(jù)

8、中存在分布均勻、一定面積的純凈地物。該區(qū)域的平均光譜沒有明顯的吸收特征，區(qū)域輻射光譜主要反映的是當時大氣條件下的太陽光譜。處理時，選擇其中某一樣區(qū)，并求出其像元的灰度平均值，然后對圖像中每一個像元都除以該平均灰度值，從 而對整幅圖像數(shù)據(jù)進行定標。(1) (1) 圖像統(tǒng)計模型法圖像統(tǒng)計模型法1）2021/7/2317內(nèi)部平均法假定一幅圖像內(nèi)部的地物充分混雜，整幅圖像的平均光譜基本代表了大氣影響下的太陽光譜信息。因而，把圖像DN值與整幅圖像的平均輻射光譜值的比值作為相對反射率，校正后的結(jié)果為相對的反射率值，這是一種基于圖像的方法。這種方法假設地面變化使充分異構(gòu)的，光譜反射特性的空間變化會相互抵消。

9、如果這個假設不成立，得到的反射光譜會有虛假性.2）2021/7/23183）這種方法假設地面目標的反射率與成像光譜儀探測器所獲得的信號之間具有線性關(guān)系。DN代表的是波段b上給定像元的灰度值，代表的是反射率。2021/7/23192021/7/23204 4）波段對比法波段對比法暗像元法暗像元法l該方法的理論依據(jù)在于大氣散射的選擇性，即大氣散射對短波的影響大，對長波的影響小。l以陸地資源衛(wèi)星TM為例，1波段收到散射影響最為嚴重，其次為2波段、3波段，而7波段受散射影響最小，為處理方便，可以把近紅外圖像看做無散射影響的標準影像，通過對不同波段的對比分析計算出大氣的干擾值。l若圖像中存在亮度為零的目

10、標亮度為零的目標，如深海水體、陰影等，則其對應圖像的亮度值應為0，實際上只有在沒有受大氣影響的情況下，其亮度值才可能為零，其他目標由于受水氣散射、輻射使得目標的亮度值不為零。2021/7/2321 1978年Ahern等利用清澈的水體影像作為暗目標提取大氣光學信息進行大氣輻射糾正。通過比較各波段直方圖，將相應波段中最小值與暗目標的差，為大氣輻射影響。顯示有關(guān)圖像的直方圖，從圖上可以得知最黑的目標亮度為零，即第七波段圖像的最小亮度值為零，第四波段的亮度最小值為a4，則a4就是第四波段圖像的大氣校正。其它波段同理可以得到大氣校正。 2021/7/23225 5）波段對比法波段對比法回歸分析法回歸分

11、析法 用長波數(shù)據(jù)來校正短波數(shù)。在不受大氣影響的波段（如TM5）和待校正的某一波段圖像中，選擇由最亮至最暗的一系列目標，將每一目標的兩個待比較的波段灰度值提取出來進行回歸分析。例如：XbaY11式中式中，X為為TM5波段的亮度均值；波段的亮度均值；Y為為TM1亮度均值；亮度均值；2021/7/2323（2 2）大氣輻射傳輸模型）大氣輻射傳輸模型2021/7/23246s模型：second simulation of the satellite signal in the solar spectrum 根據(jù)不同的成像系統(tǒng)和大氣條件，通過大氣傳輸方程對大氣輻射進行求取。輸入的參數(shù)包括： 輸出的結(jié)果1

12、、幾何參數(shù) 1、吸收2、大氣模式 2、輻射3、氣溶膠模式 3、散射4、氣溶膠濃度 4、校正5、地面高度6、探測器高度7、地表特性8、表觀反射率2021/7/2325基于大氣物理參數(shù)的模型1、lowtran模型：大氣輻射傳輸實用軟件（美國空軍地球物理實驗室開發(fā)）2、modtran模型：對lowtran的改進3、acorn模型：美國科羅拉多州ImSpec有限公司開發(fā)的一種基于圖像處理的大氣校正軟件，實現(xiàn)圖像輻射值到表觀地表反射率的轉(zhuǎn)換。2021/7/2326輻射地形校正的目的：輻射地形校正的目的：消除地形造成的影像中同消除地形造成的影像中同類地物亮度值的差異；提類地物亮度值的差異；提高影像解譯精度

13、。高影像解譯精度。4.3 4.3 輻射地形校正輻射地形校正2021/7/23272021/7/2328C糾正法在余弦糾正法的基礎上發(fā)展起來，由于在入射角比較小的區(qū)域會發(fā)生糾正過度的問題。2021/7/2329(a)原始影像(c)校正結(jié)果2(b) 校正結(jié)果12021/7/23304.4 4.4 遙感影像的幾何校正遙感影像的幾何校正l幾何畸變：遙感圖像在幾何位置上發(fā)生了變化，產(chǎn)生行列不均勻，像元大小與地面大小對應不準確，地物形狀不規(guī)則 變化。l產(chǎn)生畸變的圖像給定量分析及位置配準造成困難。在應用遙感圖像之前，必須將其準確投影到需要的坐標系中，這就是幾何糾正。2021/7/2331（1 1）遙感平臺位

14、置和運動狀態(tài)變化）遙感平臺位置和運動狀態(tài)變化遙感影像變形的原因：遙感影像變形的原因：2021/7/2332（2 2）地形起伏影響）地形起伏影響2021/7/2333（3 3）地球表面曲率的影響）地球表面曲率的影響 像點位置的移動像元對應于地面寬度的不等2021/7/2334（4 4）大氣折射的影響）大氣折射的影響 2021/7/2335（5 5）地球自轉(zhuǎn)的影響）地球自轉(zhuǎn)的影響 2021/7/2336幾何糾正分為：粗校正和精校正兩個階段l粗糾正：粗糾正：首先由接受部門進行校正，這種校正往往根據(jù)遙感平臺、地球、傳感器的各種參數(shù)糾正系統(tǒng)誤差，即把與傳感器有關(guān)的測定的校正數(shù)據(jù)帶入相應的構(gòu)像方程。l精糾

15、正：精糾正：然后，由用戶根據(jù)使用的目的不同或投影及比例尺不同，進一步做幾何精校正。如何進行幾何糾正？如何進行幾何糾正？2021/7/2337高光譜遙感圖像的精糾正高光譜遙感圖像的精糾正l基于地面控制點來消除圖像中的幾何變形，產(chǎn)生一幅符合某種地圖投影或圖形表達要求的新圖像的過程。2021/7/2338基本思路：基本思路： 校正前的圖像由于某種幾何畸變，圖像像元對應的地面距離并不相等。校正后的圖像也是由等間距的網(wǎng)格點組成的，且以地面為標準，符合某種投影的均勻分布。 因此，可以用一個適當?shù)亩囗検絹肀磉_糾正前后圖像中相應點之間的坐標關(guān)系。建立變換后圖像坐標（x,y）與變換前圖像坐標（u,v）的關(guān)系，顯

16、然，整數(shù)的像元點在原圖像坐標系中一般不在整數(shù)（u,v）點上。2021/7/2339計算方法（以多項式法為例）計算方法（以多項式法為例） 首先建立兩圖像像元點之間的對應關(guān)系，用數(shù)學關(guān)系f表示為二元N次多項式，記做： niinjjiiiijnijnjjiiiijvubyvuax0000),(),(vufyvufxyx通常數(shù)學關(guān)系表示為二元N次多項式：2021/7/2340實際計算中常常采用二元二次多項式，其展開式為：2022201101100020222011011000vbubuvbvbubbyvauauvavauaax首先要計算出上面式的12個系數(shù)，這就需要12個方程，實際工作中發(fā)現(xiàn)，6個控制

17、點是解線性方程的理論最低數(shù)，這樣少的控制點較正后的圖像效果差，因此在實際工作中，要增加控制點，控制點增加后計算結(jié)果有所改變，需采用最小二乘法，通過對控制點數(shù)據(jù)進行曲面擬合來求系數(shù)。2021/7/2341Ll 1這里代表。其中，L為控制點的個數(shù)，上式寫成距陣形式記為： 11BAU 4223332224323233222232222322221llllllllllllllllllllllllllllllllllllllllllllllllllllllvvuvuvvuvvuuvuvuuuvuvuvuvuvuvuvvuvuvvuvvuuvuvuuuvuvuvuA同樣以Y為主的距陣形式為：22BAU 0

18、220110110001aaaaaaU0220110110002bbbbbbU221llllllllllllvxuxvuxvxuxxB222llllllllllllvyuyvuyvyuyyB2021/7/2342依據(jù)最小二乘法進行求解系數(shù)矩陣U1和U2的值，并進行精度檢驗：212111)()(PBAPAAUPBAPAAUTTTT2021/7/2343系數(shù)確定后，利用公式便可以根據(jù)每一個像元點的行列值，計算每一點的亮度值。由于計算后的像點多數(shù)不在原圖的像元中心處，因此必須重新計算新位置的亮度值。一般來說，新點的亮度值介于鄰點亮度值之間，所以常常采用內(nèi)插法計算。不在采樣點2021/7/2344采樣方法：采樣方法： 1）最近鄰法 圖像中兩相臨點的距離為1，即行間距為1，列間距為1，取與所計算點（x,y）周圍相臨的4個點，比較他們與被計算點的距離，哪個點距離最近，就取哪個的亮度值作為（x,y）點的亮度值。設該最近鄰點的坐標為（K，L），則： K=Integer(x+0.5) L=Integer(y+0.5)式中，Integer為取整，于是點（K，L）的亮度值就作為點（x,y）的亮度值。優(yōu)缺點：這種方法簡單易用，計算量