地物光譜數(shù)據(jù)的重建

地物光譜數(shù)據(jù)的重建1第1頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月定義由于受到太陽(yáng)位置、角度條件、大氣條件、地形影響及傳感器本身的性能的影響，傳感器所記錄的數(shù)據(jù)與目標(biāo)的光譜反射率或光譜輻射亮度值并不一致。因此，將傳感器記錄的原始輻射值（DN值）轉(zhuǎn)化為地物反射率，恢復(fù)地物光譜數(shù)據(jù)本來(lái)的面目，稱為地物光譜數(shù)據(jù)重建。2第2頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月遙感器定標(biāo)大氣校正

光照及地形校正

遙感器入瞳輻射值地表輻射值

地表真實(shí)反射值

原始數(shù)字信號(hào)地物光譜數(shù)據(jù)重建的過(guò)程：3第3頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月輻射誤差：傳感器所得到的目標(biāo)測(cè)量值與目標(biāo)的光譜反射率或光譜輻亮度等物理量之間的差值輻射誤差產(chǎn)生原因：傳感器響應(yīng)特性和外界自然環(huán)境傳感器在光電變換的過(guò)程中，對(duì)各波段的靈敏度是有差異的，也就是說(shuō)，傳感器對(duì)各波段的光譜響應(yīng)是不同的，由此造成輻射畸變。另外，傳感器的光學(xué)鏡頭的非均勻性，會(huì)引起邊緣減光，也會(huì)造成圖像輻射的畸變。

地物（目標(biāo)物）的輻射（反射）經(jīng)過(guò)大氣層時(shí)，與大氣層發(fā)生散射作用和吸收作用。吸收作用直接降低地物的輻射能量，引起輻射畸變。散射作用除降低地物的輻射能量外，大氣散射的部分輻射還會(huì)進(jìn)入傳感器，直接疊加在目標(biāo)地物的輻射能量之中，成為目標(biāo)地物的噪聲，降低了圖像的質(zhì)量。光照條件的不同也會(huì)引起輻射畸變，如太陽(yáng)高度角、地面坡度等，都會(huì)引起輻射的畸變。輻射誤差4第4頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月4.1遙感器定標(biāo)一般把消除傳感器本身影響的輻射校正過(guò)程稱之為輻射定標(biāo)，即遙感器定標(biāo)。對(duì)高光譜定量遙感具有十分重要的意義，只有經(jīng)過(guò)精確的定標(biāo)后，才能獲取真實(shí)的地物物理參量，才能與不同的高光譜遙感數(shù)據(jù)，模擬數(shù)據(jù)進(jìn)行對(duì)比分析。具體形式：建立成像光譜儀每個(gè)探測(cè)元件輸出的數(shù)字量化值（DN）與它對(duì)應(yīng)視場(chǎng)中輸出輻射亮度值之間的定量關(guān)系。5第5頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月遙感器輻射定標(biāo)是指建立遙感器每個(gè)探測(cè)元所輸出信號(hào)的數(shù)值量化值與該探測(cè)器對(duì)應(yīng)像元內(nèi)的實(shí)際地物輻射亮度值之間的定量關(guān)系。Li波段i的入瞳輻射能量DNi波段i傳感器的輸出亮度值A(chǔ)i波段i的定標(biāo)增益系數(shù)Bi

波段i的定標(biāo)偏置量6第6頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月遙感器定標(biāo)主要包括三個(gè)階段：實(shí)驗(yàn)室定標(biāo)和星上定標(biāo)，場(chǎng)地定標(biāo)（1）實(shí)驗(yàn)室定標(biāo)：是成像光譜儀在運(yùn)行前在實(shí)驗(yàn)室所接受的定標(biāo)過(guò)程。包括光譜定標(biāo)和輻射定標(biāo)光譜定標(biāo)：確定成像光譜儀每個(gè)波段的中心波長(zhǎng)和帶寬。7第7頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月輻射定標(biāo)分為兩種形式：相對(duì)定標(biāo)和絕對(duì)定標(biāo)。相對(duì)定標(biāo)：校正遙感器中各個(gè)探測(cè)器元件響應(yīng)度差異而對(duì)衛(wèi)星遙感器測(cè)量到的原始數(shù)字進(jìn)行歸一化的工作。絕對(duì)定標(biāo)：通過(guò)各種標(biāo)準(zhǔn)輻射源，在不同波譜段建立成像光譜儀入瞳處光譜輻射亮度值與成像儀輸出的數(shù)字值（DN）之間的定量關(guān)系。8第8頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月其他系統(tǒng)誤差校正常見(jiàn)的由遙感系統(tǒng)引入的輻射誤差有：隨機(jī)壞像元（散粒噪聲）行起始/終止問(wèn)題行或列部分缺失行或列條紋9第9頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月隨機(jī)壞像元(散粒噪聲)BV1BV2BV3BV8BVijkBV4BV7BV6BV5原因：一個(gè)探測(cè)器未記錄某個(gè)光譜數(shù)據(jù)，該像元成為壞像元算法步驟：計(jì)算相鄰8個(gè)像元的均值，用均值代替壞像元的灰度值。10第10頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月行或列缺失原因：掃描系統(tǒng)的某個(gè)探測(cè)器工作不正常，就有可能產(chǎn)生一整行沒(méi)有光譜信息的線。對(duì)于線陣列掃描系統(tǒng)而言，就會(huì)使得整行或者整列數(shù)據(jù)缺失。解決方法：人工處理。11第11頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月（2）機(jī)上或星上定標(biāo)：飛行之前的實(shí)驗(yàn)室定標(biāo)會(huì)隨著環(huán)境的變化而變化，由于受到發(fā)射的影響和空間工作的溫度變化，發(fā)射前的定標(biāo)關(guān)系發(fā)生變化，因此，實(shí)驗(yàn)室定標(biāo)后飛機(jī)或衛(wèi)星上的定標(biāo)也是必要的。機(jī)上定標(biāo)的內(nèi)容和形式基本上與實(shí)驗(yàn)室定標(biāo)是類似的，其中穩(wěn)定的光源是主要因素。12第12頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月（3）地面場(chǎng)定標(biāo)地面場(chǎng)定標(biāo)：在遙感器處于正常運(yùn)行條件下，選擇輻射的定標(biāo)場(chǎng)地，通過(guò)地面同步測(cè)量對(duì)遙感進(jìn)行定標(biāo)。對(duì)前面兩項(xiàng)定標(biāo)進(jìn)行檢驗(yàn)和修正。13第13頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月場(chǎng)地定標(biāo)的原理：遙感器升空后，在定標(biāo)場(chǎng)地選擇若干像元區(qū)，測(cè)量成像光譜儀對(duì)應(yīng)地物的各波段光譜反射率和大氣光譜參量，利用大氣輻射傳輸模型給出成像光譜儀入瞳處各譜帶的輻射亮度，最后確定它與成像光譜儀對(duì)應(yīng)的輸出的數(shù)字量化值的關(guān)系，求解定標(biāo)系數(shù)，并進(jìn)行誤差分析。14第14頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月4.2大氣校正大氣輻射傳輸理論：大氣散射集中在太陽(yáng)輻射能量最強(qiáng)的可見(jiàn)光區(qū)。因此，散射是太陽(yáng)輻射衰減的主要原因。大氣對(duì)電磁輻射產(chǎn)生的散射和吸收是主要的輻射誤差源，因此，遙感圖像的大氣校正主要是補(bǔ)償大氣在吸收與散射方面的瞬時(shí)影響。大氣校正方法可以大致分為：（1）基于圖像統(tǒng)計(jì)模型；（2）專用的大氣校正模型15第15頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月這是一種相對(duì)定標(biāo)方法，已廣泛地應(yīng)用在成像光譜遙感數(shù)據(jù)處理中。模型要求在處理的影像數(shù)據(jù)中存在分布均勻、一定面積的純凈地物。該區(qū)域的平均光譜沒(méi)有明顯的吸收特征，區(qū)域輻射光譜主要反映的是當(dāng)時(shí)大氣條件下的太陽(yáng)光譜。處理時(shí)，選擇其中某一樣區(qū)，并求出其像元的灰度平均值，然后對(duì)圖像中每一個(gè)像元都除以該平均灰度值，從而對(duì)整幅圖像數(shù)據(jù)進(jìn)行定標(biāo)。(1)圖像統(tǒng)計(jì)模型法1）16第16頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月內(nèi)部平均法假定一幅圖像內(nèi)部的地物充分混雜，整幅圖像的平均光譜基本代表了大氣影響下的太陽(yáng)光譜信息。因而，把圖像DN值與整幅圖像的平均輻射光譜值的比值作為相對(duì)反射率，校正后的結(jié)果為相對(duì)的反射率值，這是一種基于圖像的方法。這種方法假設(shè)地面變化使充分異構(gòu)的，光譜反射特性的空間變化會(huì)相互抵消。如果這個(gè)假設(shè)不成立，得到的反射光譜會(huì)有虛假性.2）17第17頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月3）這種方法假設(shè)地面目標(biāo)的反射率與成像光譜儀探測(cè)器所獲得的信號(hào)之間具有線性關(guān)系。DN代表的是波段b上給定像元的灰度值，ρ代表的是反射率。18第18頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月19第19頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月4）波段對(duì)比法——暗像元法該方法的理論依據(jù)在于大氣散射的選擇性，即大氣散射對(duì)短波的影響大，對(duì)長(zhǎng)波的影響小。以陸地資源衛(wèi)星TM為例，1波段收到散射影響最為嚴(yán)重，其次為2波段、3波段，而7波段受散射影響最小，為處理方便，可以把近紅外圖像看做無(wú)散射影響的標(biāo)準(zhǔn)影像，通過(guò)對(duì)不同波段的對(duì)比分析計(jì)算出大氣的干擾值。若圖像中存在亮度為零的目標(biāo)，如深海水體、陰影等，則其對(duì)應(yīng)圖像的亮度值應(yīng)為0，實(shí)際上只有在沒(méi)有受大氣影響的情況下，其亮度值才可能為零，其他目標(biāo)由于受水氣散射、輻射使得目標(biāo)的亮度值不為零。20第20頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月

1978年Ahern等利用清澈的水體影像作為暗目標(biāo)提取大氣光學(xué)信息進(jìn)行大氣輻射糾正。通過(guò)比較各波段直方圖，將相應(yīng)波段中最小值與暗目標(biāo)的差，為大氣輻射影響。顯示有關(guān)圖像的直方圖，從圖上可以得知最黑的目標(biāo)亮度為零，即第七波段圖像的最小亮度值為零，第四波段的亮度最小值為a4，則a4就是第四波段圖像的大氣校正。其它波段同理可以得到大氣校正。

21第21頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月5）波段對(duì)比法——回歸分析法

用長(zhǎng)波數(shù)據(jù)來(lái)校正短波數(shù)。在不受大氣影響的波段（如TM5）和待校正的某一波段圖像中，選擇由最亮至最暗的一系列目標(biāo)，將每一目標(biāo)的兩個(gè)待比較的波段灰度值提取出來(lái)進(jìn)行回歸分析。例如：式中，為T(mén)M5波段的亮度均值；為T(mén)M1亮度均值；22第22頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月（2）大氣輻射傳輸模型23第23頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月6s模型：secondsimulationofthesatellitesignalinthesolarspectrum根據(jù)不同的成像系統(tǒng)和大氣條件，通過(guò)大氣傳輸方程對(duì)大氣輻射進(jìn)行求取。輸入的參數(shù)包括：輸出的結(jié)果1、幾何參數(shù)1、吸收2、大氣模式2、輻射3、氣溶膠模式3、散射4、氣溶膠濃度4、校正5、地面高度6、探測(cè)器高度7、地表特性8、表觀反射率24第24頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月基于大氣物理參數(shù)的模型1、lowtran模型：大氣輻射傳輸實(shí)用軟件（美國(guó)空軍地球物理實(shí)驗(yàn)室開(kāi)發(fā)）2、modtran模型：對(duì)lowtran的改進(jìn)3、acorn模型：美國(guó)科羅拉多州ImSpec有限公司開(kāi)發(fā)的一種基于圖像處理的大氣校正軟件，實(shí)現(xiàn)圖像輻射值到表觀地表反射率的轉(zhuǎn)換。25第25頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月輻射地形校正的目的：消除地形造成的影像中同類地物亮度值的差異；提高影像解譯精度。4.3輻射地形校正26第26頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月27第27頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月C糾正法在余弦糾正法的基礎(chǔ)上發(fā)展起來(lái)，由于在入射角比較小的區(qū)域會(huì)發(fā)生糾正過(guò)度的問(wèn)題。28第28頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月(a)原始影像(c)校正結(jié)果2(b)校正結(jié)果129第29頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月4.4遙感影像的幾何校正幾何畸變：遙感圖像在幾何位置上發(fā)生了變化，產(chǎn)生行列不均勻，像元大小與地面大小對(duì)應(yīng)不準(zhǔn)確，地物形狀不規(guī)則變化。產(chǎn)生畸變的圖像給定量分析及位置配準(zhǔn)造成困難。在應(yīng)用遙感圖像之前，必須將其準(zhǔn)確投影到需要的坐標(biāo)系中，這就是幾何糾正。30第30頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月（1）遙感平臺(tái)位置和運(yùn)動(dòng)狀態(tài)變化遙感影像變形的原因：31第31頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月（2）地形起伏影響32第32頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月（3）地球表面曲率的影響

像點(diǎn)位置的移動(dòng)像元對(duì)應(yīng)于地面寬度的不等33第33頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月（4）大氣折射的影響

34第34頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月（5）地球自轉(zhuǎn)的影響

35第35頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月幾何糾正分為：粗校正和精校正兩個(gè)階段粗糾正：首先由接受部門(mén)進(jìn)行校正，這種校正往往根據(jù)遙感平臺(tái)、地球、傳感器的各種參數(shù)糾正系統(tǒng)誤差，即把與傳感器有關(guān)的測(cè)定的校正數(shù)據(jù)帶入相應(yīng)的構(gòu)像方程。精糾正：然后，由用戶根據(jù)使用的目的不同或投影及比例尺不同，進(jìn)一步做幾何精校正。如何進(jìn)行幾何糾正？36第36頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月高光譜遙感圖像的精糾正基于地面控制點(diǎn)來(lái)消除圖像中的幾何變形，產(chǎn)生一幅符合某種地圖投影或圖形表達(dá)要求的新圖像的過(guò)程。37第37頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月基本思路：校正前的圖像由于某種幾何畸變，圖像像元對(duì)應(yīng)的地面距離并不相等。校正后的圖像也是由等間距的網(wǎng)格點(diǎn)組成的，且以地面為標(biāo)準(zhǔn)，符合某種投影的均勻分布。因此，可以用一個(gè)適當(dāng)?shù)亩囗?xiàng)式來(lái)表達(dá)糾正前后圖像中相應(yīng)點(diǎn)之間的坐標(biāo)關(guān)系。建立變換后圖像坐標(biāo)（x,y）與變換前圖像坐標(biāo)（u,v）的關(guān)系，顯然，整數(shù)的像元點(diǎn)在原圖像坐標(biāo)系中一般不在整數(shù)（u,v）點(diǎn)上。38第38頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月計(jì)算方法（以多項(xiàng)式法為例）首先建立兩圖像像元點(diǎn)之間的對(duì)應(yīng)關(guān)系，用數(shù)學(xué)關(guān)系f表示為二元N次多項(xiàng)式，記做：

通常數(shù)學(xué)關(guān)系表示為二元N次多項(xiàng)式：39第39頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月實(shí)際計(jì)算中常常采用二元二次多項(xiàng)式，其展開(kāi)式為：首先要計(jì)算出上面式的12個(gè)系數(shù)，這就需要12個(gè)方程，實(shí)際工作中發(fā)現(xiàn)，6個(gè)控制點(diǎn)是解線性方程的理論最低數(shù)，這樣少的控制點(diǎn)較正后的圖像效果差，因此在實(shí)際工作中，要增加控制點(diǎn)，控制點(diǎn)增加后計(jì)算結(jié)果有所改變，需采用最小二乘法，通過(guò)對(duì)控制點(diǎn)數(shù)據(jù)進(jìn)行曲面擬合來(lái)求系數(shù)。40第40頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月這里代表。其中，L為控制點(diǎn)的個(gè)數(shù)，上式寫(xiě)成距陣形式記為：

同樣以Y為主的距陣形式為：

41第41頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月依據(jù)最小二乘法進(jìn)行求解系數(shù)矩陣U1和U2的值，并進(jìn)行精度檢驗(yàn)：212111)()(PBAPAAUPBAPAAUTTTT--==42第42頁(yè)，課件共49頁(yè)，創(chuàng)作于2023年2月系數(shù)確定后，利用公式便可以根據(jù)每一個(gè)像元點(diǎn)的行列值，計(jì)算每一點(diǎn)的亮度值。由于計(jì)算后的像點(diǎn)多數(shù)不在原圖的像元中心處，因此必須重新計(jì)算新位置的亮度值。一般來(lái)說(shuō)，新點(diǎn)的亮度值介于鄰點(diǎn)亮度值之間，所以常常采用內(nèi)插法計(jì)算。不在采樣點(diǎn)43第43頁(yè)，課件共49頁(yè)