sample_of_atm_cor教學(xué)教材

1、Good is good, but better carries it.精益求精，善益求善。sample_of_atm_cor-基于TM/ETM+遙感數(shù)據(jù)的地面相對反射率反演摘要：衛(wèi)星遙感數(shù)據(jù)的大氣輻射校正是定量遙感的基礎(chǔ)與關(guān)鍵。目前大氣輻射校正方法的局限性已經(jīng)嚴(yán)重影響了定量遙感科學(xué)進展。鑒于此，COST模型在TM/ETM數(shù)據(jù)上的實用性逐步得到學(xué)者們的共識。本文詳細(xì)闡述了COST模型的理論及其算法，模型參數(shù)的確定及其獲取途徑。COST模型在完成TM/ETM衛(wèi)星遙感數(shù)據(jù)大氣輻射糾正的同時實現(xiàn)了地表反射率反演，為TM/ETM衛(wèi)星遙感數(shù)據(jù)的應(yīng)用研究奠定了基礎(chǔ)。文章進一步借助地面光譜儀的實測數(shù)據(jù)對研究

2、區(qū)地表反射率反演結(jié)果進行了空間維和光譜維的分析驗證，結(jié)果認(rèn)為COST模型完全能夠滿足衛(wèi)星遙感數(shù)據(jù)基本應(yīng)用研究的需要。關(guān)鍵詞：大氣輻射校正；COST模型；地表相對反射率；衛(wèi)星遙感數(shù)據(jù)RETRIEVEOFEARTHREFLECTIONBASEONTHEREMOTESENSINGDATAOFTM/ETMAbstract:Atmosphericandradiometriccorrectionisthebasisandkeyforquantificationalremotesensing.Nowadaythelimitationofthemethodofatmosphericandradiometric

3、correctionhadbeenstrictlyeffectingthedevelopmentofRSscience.Hereinabovestatement,commonunderstandinghadbeenbuildingontheutilityofCOSTmodel.ThepaperexpatiatedtheCOSTmodelmethodandarithmetic,andhowtoconfirmtheparameters.WithCOSTmodelfinishingatmosphericandradiometriccorrectionofTM/ETM,retrieveofearthr

4、eflectionhadbeencarriedout.ThisworkestablishedthefoundationofapplicationstudyoftheremotesensingdataofTM/ETM.Furthermorethepaperanalysicedandtestedtheretrievedresultfromspationalandspectraldimentionsbythesurveyingdataofgroundspectrumapparatus.ItisconsideringthattheCOSTmodelcompletelymeettheneedofappl

5、icationstudyofremotesensingdataofTM/ETM.KeyWords:atmosphericandradiometriccorrection;COSTmodel;earthrelativereflection;LandsatThematicMapperdata1基于TM/ETM+遙感數(shù)據(jù)的COST模型與理論遙感圖像輻射值受來自大氣層的影響是多方面的。就被動傳感器而言，一方面大氣層的吸收和散射使來自地物目標(biāo)的輻射能量在到達(dá)傳感器之前因衰減而降低；另一方面，大氣反射和散射形成的路徑輻射和地物目標(biāo)輻射一起進入遙感探測器，導(dǎo)致遙感輻射量失真，質(zhì)量和對比度下降，猶如蒙上一層薄

6、霧。因此，大氣校正是遙感數(shù)據(jù)預(yù)處理的基礎(chǔ)和應(yīng)用的關(guān)鍵。盡管大氣輻射傳輸原理是相同的，但由于傳感器本身性能和參數(shù)的不同，不同傳感器獲得的遙感數(shù)據(jù)在大氣校正的具體方法上有所差別，就Lansant5和Landsat7（六波段除外）的數(shù)據(jù)而言，目前國內(nèi)外學(xué)者在大氣校正時應(yīng)用較多的是黑體消除法，其中Crist等和Chavez就黑體消除法提出了各自行之有效的模型和算法。Crist等假設(shè)研究區(qū)具有相同的黑體背景，利用TC（TasseledCap）變換的第四成分作為黑體輻射指標(biāo)(Crist，1986)，以此去除黑體輻射不失為一種有潛力的方法。傳統(tǒng)的黑體消除方法認(rèn)為近紅外輻射不存在散射影響，因此近紅外波段上測得

7、的任何潔凈深水體或陰影面積輻射值大于零均被認(rèn)為是大氣散射和路徑輻射的結(jié)果，即近紅外波段的灰度值可以近似于大氣層輻射值，以此來推算其它波段的大氣層輻射值和有關(guān)的大氣性質(zhì)參數(shù)，這種黑體去除常常會矯枉過正，不適于黑體不存在的影像(Chavez,1996)。鑒于此，COST模型基于以下的理論假設(shè)：一是每一波段均存在反射率為1的黑體輻射，黑體輻射值取決于大氣頂層的太陽輻照度；二是大氣性質(zhì)是均一的，傳感器每一波段的最小輻射亮度值除黑體輻射影響外，主要還有大氣分子的瑞利散射和氣溶膠的米氏散射和反射作用的影響。其概念模型概括為以下三個步驟：一是將遙感器記錄的DN值轉(zhuǎn)換為遙感器的光譜輻射值，即根據(jù)遙感器的增益與

8、偏移進行遙感器定標(biāo)；二是遙感器的光譜輻射值轉(zhuǎn)換為遙感器的相對反射值；三是消除因大氣吸收和散射造成的大氣影響，即大氣糾正，同時計算地球表面像元相對反射率（見圖1COST大氣糾正概念模型）。圖1COST大氣糾正概念模型Fig1ConceptmodelflowofatmosphericcorrectionforCOST根據(jù)太陽輻射和大氣傳輸原理與過程，TM/ETM+數(shù)據(jù)地面反射率反演的數(shù)學(xué)模型(Chavez,1996)34可綜合表達(dá)為：=D2(LsatI-LhazeI)/ESUNICOS2(SZ)（1）其中：地面相對反射率；D日地天文單位距離；LsatI傳感器光譜輻射值，即大氣頂層的輻射能量；Lha

9、zeI大氣層輻射值；ESUNl大氣頂層的太陽平均光譜輻射，即大氣頂層太陽輻照度；SZ太陽天頂角。2實驗遙感數(shù)據(jù)與模型參數(shù)的確定方法2.1實驗區(qū)遙感數(shù)據(jù)特征實驗區(qū)共采用了四個時相的多波段遙感數(shù)據(jù)，和Lansat7數(shù)據(jù)獲取時間分別為1990年08月29日和1991年06月29日（Landsat5，TM），2001年05月31日和2002年09月23日（Landsat7，ETM），遙感影像的空間分辨率均為30米。產(chǎn)品校正水平為系統(tǒng)校正，采用橫軸墨卡托投影（TramsverseMercator)和克拉索夫斯基(Krassovsky)橢球體建立的北京54基準(zhǔn)面進行幾何校正。L5和L7光譜通道特征見表3.

10、1。表1L5和L7光譜通道特征Table1TheCharacteristicsofSpectralChanlesforL5andL7波段號BandNumber波段類型BandStypesL7ETM+L5TM頻譜范圍Bandregion（m）波段寬度Bandwidth（m）空間分辨率SpationResolution（IFOV，m）頻譜范圍Bandregion（m）波段寬度Bandwidth（m）空間分辨率SpationResolution（IFOV，m）B1Blue-Green0.450-0.5150.066300.45-0.520.0730B2Green0.252-0.6050.082300

11、52-0.600.0830B3Red0.630-0.690.067300.63-0.690.0630B4NearIR0.775-0.900.128300.76-0.900.1430B5SWIR1.550-1.750.217301.55-1.750.230B6LWIR10.40-12.516010.40-12.52.1120B7SWIR2.090-2.350.252302.08-2.350.2730B8Pan0.520-0.900.28152.2模型常規(guī)參數(shù)確定模型（1）中相關(guān)參數(shù)可以通過不同途徑獲得。其中：日地天文單位距離D1-0.01674cos(0.9856(JD-4)/180)；JD為遙

12、感成像的儒略日（JulianDay），可以通過儒略日算法模塊獲得；太陽天頂角90太陽高度角；太陽高度角可以從遙感數(shù)據(jù)的頭文件中獲得；大氣頂層太陽輻照度（ESUNI）從遙感權(quán)威單位定期測定并公布的信息中獲取。實驗區(qū)上述遙感數(shù)據(jù)參數(shù)見遙感數(shù)據(jù)技術(shù)參數(shù)表（表3.2）。表2遙感數(shù)據(jù)技術(shù)參數(shù)（一）Table2TheImageDataParameters(one)數(shù)據(jù)類型DataType成像時間（年月日）Imagetime(YYMMDD)太陽高度角（度，SE）SolarElevationAngle太陽天頂角（度，SZ）SolarZenithAngle儒略日J(rèn)ulianDay日地距離（天文單位，D）Eart

13、h-SunDistanceinastronomicalunitsL5TM1990082951.460338.53972448132.51.015221641991062961.005728.99432448436.51.00146796L7ETM+2001053165.0025.00002452060.50.993417552002092347.5042.50002452540.51.01673999Note:JDhttp:/wise-obs.tau.ac.il/eran/Wise/Util/Julian_Day.html。表3遙感數(shù)據(jù)技術(shù)參數(shù)（二)大氣頂層太陽輻照度（ESUNI）Table3

14、SolarExoatmosphericSpectralIrradiances(ESUNI)(W/cm2.ster.m)BandTM1TM2TM3TM4TM5TM7L5TM195.7182.9155.7104.721.937.452L7ETM+196.9184155.1104.422.578.2072.3遙感器光譜輻射定標(biāo)由遙感器的靈敏度特征引起的輻射畸變主要由其光學(xué)系統(tǒng)或光電轉(zhuǎn)換系統(tǒng)的特征形成的，光電轉(zhuǎn)換系統(tǒng)的靈敏性特征通常很重復(fù)，其校正一般是通過定期的地面測定值進行的。Lansat5和7系列的遙感器糾正是通過飛行前實地測量，預(yù)先測出了各個波段的輻射值和記錄值之間的校正增益系數(shù)和校正偏移量。遙

15、感器光譜輻射定標(biāo)時采用以下轉(zhuǎn)換算式：LsatI=Bias+(GainDN)（2）單位：mWcm-ster-1mm-1(forLandsat)其中：Bias偏移；Gain增益；DN象元值。從頭文件中獲取的偏置和增益值見L5和L7各波段光譜通道的增益和偏置表（表3.3）。表4L5和L7各波段光譜通道的增益和偏置單位：W/m2.ster.mTable4theGainsandBiasesofL5andL7unit:W/m2.ster.m1750波段號BandNumberL5TM1990/1991L7ETM+2001/2002增益GAINS偏置BIASES增益GAINS偏置BIASESB10.99992

16、-0.010001.18070871-7.38070852B22.42430-0.023201.20984250-7.60984259B31.36344-0.007800.94251966-5.94251966B42.62901-0.019300.96929136-6.06929127B50.58771-0.008000.19122047-1.19122047B63.201070.25994B70.38674-0.004000.06649607-0.416496062.4大氣校正任何一種依賴大氣物理模型的大氣校正方法都需要先進行遙感器的輻射校準(zhǔn)。大氣對光學(xué)遙感的影響是很復(fù)雜的。Chavez等提

17、出了基于圖像本身的大氣參數(shù)估計方法的大氣校正COST模型為：Lhazel=LI,min-LI,1%（Chavez,1996）（3）其中：Lhazel大氣層光譜輻射值；LI,min遙感器每一波段最小光譜輻射值；LI,1%反射率為1%的黑體輻射值。遙感器的最小光譜輻射值LI，min的轉(zhuǎn)換算式為：LI，min=LMINI+QCAL（LMAXI-LMINI）/QCALMAX（4）其中：QCAL每一波段最小DN值（亮度值）；QCALMAX255；LMAXl、LMINl為常數(shù)，指遙感器光譜輻射值的上限和下限，從遙感數(shù)據(jù)頭文件或權(quán)威部門定期公布信息中獲取，本實驗參數(shù)見遙感器光譜輻射動力范圍及其大氣輻射定標(biāo)值

18、表（表3.4）。表5遙感器光譜輻射動力范圍及其大氣輻射定標(biāo)值Table5theDynamicRangesofSpectralRadiometricandAtmosphericCorrection波段號BandNumberLmin/LmaxLI,hazeL5TM1990/1991L7ETM+2001/2002L5TM1990L5TM1991L7ETM+2001L7ETM2002B1-0.01/253.9697-73807/292.548357.1181-4.11574.457612.4341B2-0.0232/615.7491-7.6098/299.720048.7808-2.37651.216

19、32.3936B3-0.0078/346.3060-5.9425/233.480823.4753-2.7087-3.0283-3.2058B4-0.0193/667.7493-6.0693/240.154515.8378-0.1300-4.0374-4.4553B5-0.008/149.2704-1.192/47.38352.40090.0057-1.2060-1.3146B60.2599/813.3307B7-0.004/98.2280-0.4165/16.47510.74860.3618-0.2983-0.2435表6遙感器光譜輻射動力范圍及其大氣輻射定標(biāo)值（續(xù)）Table6theDyna

20、micRangesofSpectralRadiometricandAtmosphericCorrection波段號BandNumberLI，minLI，1%L5TM1990L5TM1991L7ETM+2001L7ETM2002L5TM1990L5TM1991L7ETM+2001L7ETM2002B161.74190.98607.909818.49554.62395.10173.45226.0614B253.10232.39164.44238.05804.32144.76813.22605.6643B327.15411.3053-0.30901.56883.67884.05902.71934.7

21、746B418.31162.5994-2.2070-1.24142.47382.72941.83043.2139B52.91900.5774-0.8102-0.61980.51810.57170.39570.6948B6386.07002.4485B70.76640.3812-0.2840-0.21780.01790.01950.01430.0257黑體輻射值LI，1%的轉(zhuǎn)換算式：LI，1%=0.01*ESUNI*COS（SZ）/（*D）（5）LI，1%假設(shè)黑體反射率為1各波段的黑體輻射值（Chavez，1996；Moran等，1992）45。將各項參數(shù)代入上述公式（3）、（4）、（5）分別計

22、算得出遙感器1、2、3、4、5、7波段的最小光譜輻射值、黑體輻射值和大氣輻射定標(biāo)值（見表3.4遙感器光譜輻射動力范圍及其大氣輻射定標(biāo)值）；將（2）、（3）公式的計算結(jié)果及其它相應(yīng)參數(shù)代入（1）式便獲得TM/ETM遙感數(shù)據(jù)的地面相對反射率。3結(jié)果分析3.1波段相關(guān)性對比分析與驗證分別在1990年的TM的DN值影像和反射值影像上取6060像元樣本，求各波段間的相關(guān)系數(shù)。結(jié)果表明，反演后各波段地表相對反射率的相關(guān)系數(shù)總體上均大于遙感影像各波段DN值的相關(guān)系數(shù)（見表4.1）。其中較為獨立的第四波段與其它波段的反射率相關(guān)性顯著提高，其次是第七、第五波段。可見COST模型校正效果是顯著的。另外，將ASD-

23、FieldSpec-FRTM光譜儀野外實測的34個外業(yè)裸土樣點的波譜按Lansat5、7的波譜范圍匯總其反射值，并進行波段間的相關(guān)性分析，結(jié)果發(fā)現(xiàn)，地面光譜儀實測波段反射率之間的相關(guān)性與TM、ETM反射率反演后的波段相關(guān)性具有較高的一致性，尤其是1、2、3波段之間（見表4.2）。這可能是因為大氣層的選擇性散射強度與波長有關(guān)，波長越長，大氣層中由分子引起的瑞利散射減弱，由氣溶膠引起的米氏散射增強。一般氣候條件下，瑞利散射起主導(dǎo)作用。因此，從波段相關(guān)性對比分析角度講，COST模型對大氣層散射作用校正的效果還是較為突出的。表7TM90遙感影像各波段的泊松相關(guān)系數(shù)Table7PearsonCorrel

24、ationforBandsofTM90band1band2band3band4band5band7band110.9888*0.9898*0.7739*0.9646*0.9719*band20.9036*10.9922*0.7515*0.9507*0.9643*band30.9219*0.9376*10.7303*0.9529*0.9742*band40.0689*0.2742*0.0480*10.8256*0.7407*band50.8069*0.8384*0.7574*0.4011*10.9853*band70.8779*0.8763*0.8563*0.2012*0.9540*1*Corr

25、elationissignificantatthe0.01level(2-tailed).*Correlationissignificantatthe0.05level(2-tailed).注：對角線左下方為TM90影像DN值的相關(guān)系數(shù)；對角線右上方為TM90反射率的相關(guān)系數(shù)。表8地面光譜儀實測波譜合并后各波段反射率的相關(guān)系數(shù)Table8PearsonCorrelationforGroundSurveySpectrumsb1b2b3b4b5b7b110.9888*0.9898*0.7739*0.9646*0.9719*b20.9884*10.9922*0.7515*0.9507*0.9643*

26、b30.9842*0.9867*10.7303*0.9529*0.9742*b40.8302*0.8612*0.7937*10.8256*0.7407*b50.8701*0.8650*0.8995*0.7944*10.9853*b70.8555*0.8221*0.8945*0.6003*0.9421*1*Correlationissignificantatthe0.01level(2-tailed).注：對角線左下方為地面波譜合并后反射率的相關(guān)系數(shù)；對角線右上方為TM90反射率的相關(guān)系數(shù)。3.2影像特征值的對比分析標(biāo)準(zhǔn)差反映出影像像元值的分散程度，存在黑體的影像（如山體陰影），像元值的分散度高

27、，COST模型能夠消除這種黑體，剔除像元極值，降低分散度。分別在四個時點的DN值和反射率影像上相同位置取6060樣點，匯總分析各波段影像樣點的像元值標(biāo)準(zhǔn)差。從分析結(jié)果看出：各波段反射值標(biāo)準(zhǔn)差遠(yuǎn)遠(yuǎn)小于DN值標(biāo)準(zhǔn)差（見表4.3）?？梢?，COST大氣校正能夠收斂像元值域，起到消除黑體的作用。表9TM/ETM+影像DN值與反射率的標(biāo)準(zhǔn)差Table9theSDofDNValueandRefleotimfortheTM/ETM+時間time19900829199106292001053120020923平均average1波段Band1DN值DNvalue12.4711.3516.9134.9918.93

28、反射率Reflection1.657.703.556.824.932波段Band2DN值DNvalue11.379.8121.0737.0819.83反射率Reflection47.936.283波段Band3DN值DNvalue22.8720.0524.8337.6226.34反射率Reflection5.0511.375.267.427.284波段Band4DN值DNvalue15.1212.8525.0237.6722.67反射率Reflection5.0516.158.0811.3710.165波段Band5DN值DNvalue33.7631.0237.0739.

29、1635.25反射率Reflection7.2217.4610.9210.7811.607波段Band7DN值DNvalue24.3923.7438.5139.7331.59反射率Reflection7.3215.3110.8210.4610.983.3光譜特征對比分析與驗證首先基于知識結(jié)合野外調(diào)查人機交互判讀出影像上的土地利用類型，然后在ENVI平臺上分別對DN值和反射率影像進行地物光譜樣點采集，每種地類的樣點數(shù)為30。TM/ETM四個時期的影像上不同地類的光譜匯總結(jié)果見圖2至圖5。匯總結(jié)果表明，COST模型糾正后的地物光譜基本反映出Landsat數(shù)據(jù)的特點，即水田、旱地和林地在1、3波段光

30、譜吸收谷明顯；水田、旱地在2、4波段光譜反射峰較為顯著；沙地、旱地、林地含水量依次降低，在5波段反射率依次升高；研究區(qū)土壤發(fā)育于黃土母質(zhì)，成土作用微弱，其性狀與母質(zhì)相似，化學(xué)成分以SiO2和AL2O為主，粘土礦物以水云母為主，表土多為粉砂輕壤質(zhì)，粒狀或不穩(wěn)定的團粒結(jié)構(gòu)，有機質(zhì)含量低，因此采樣區(qū)的旱地、林地的地表巖性和粘土礦物基本相同，覆蓋類型不同但覆蓋度差別不大，故在7波段的光譜特征趨向一致；沙地表現(xiàn)為風(fēng)成沙母質(zhì)出露，粘土礦物微弱，表土多為砂質(zhì)或礫質(zhì)，基本無覆蓋或覆蓋度低于5%，故其光譜反射率較高；研究區(qū)水田多分布在河套川地，是長期利用富含泥沙的河水灌溉、淤積和耕作施肥交替作用下形成的，水云母含量高，表土多為中壤質(zhì)，有機質(zhì)含量高，團粒結(jié)構(gòu)穩(wěn)定，光譜反射率低。圖3橫山TM1991地物DN值與反射率Figure3spectralcurve