定量遙感六章大氣糾正

定量遙感

QuantitativeRemoteSensingofLandSurface光學圖像的大氣糾正主要內(nèi)容引言；單視角圖像的糾正方法；多視角圖像的糾正方法；估算整層大氣水汽含量的方法；最新進展；小結(jié)6.1引言大氣糾正：從遙感圖像獲取地表反射率的過程；兩部分：大氣參數(shù)估計和地表反射率反演；大氣效應：分子和氣溶膠的散射和氣體吸收。分子散射和臭氧、氧氣引起的吸收容易糾正，難點是從圖像上估算氣溶膠和水汽的空間分布。包括相對糾正：以某一清晰圖為參考，將其他圖糾正到該參考圖的大氣狀況上；絕對糾正：去掉大氣影響。AdiagramshowingatmospherehowtoaffectonRSdataAtmosphericpathradianceIrradiationfromsunE0TotalradiancefromdifferentsourcesAtmospherelayerSurfacereflectionSunSensorTransmittanceTransmittance

Radiativetransfer(RT)modelTheat-sensorradianceLcanberelatedtotheLambertiansurfacereflectancer

byWhere,

:atmosphericpathradiance,:twowaytransmittanceforthesun-surface-sensorpath,S:sphericalalbedooftheatmosphere,

:exoatmosphericsolarirradiance,and:solarzenithangle.

Applythesolutiontocorrectradianceand retrievereflectanceUsingLimg–,simplycalculatecorrected radianceforsatelliteimagesWithknown ,thesurface reflectancercanberetrievedusingtheRT model.

Outputreflectanceimageswithinputoforiginal radiancefromsatelliteimage.利用6S模型進行大氣訂正：模型輸入：

landsatTM第2波段（0.500~0.650um），北京，1月12日上午10：30點，能見度23km，地表朗伯，地表反射率0.5；第一步：利用6S模型進行大氣訂正：模型輸出：可以從輸出結(jié)果得到以下參數(shù)值：E427.306T0.872557S0.11953La15.254在地表真實反射率為0.5的情況下，6S計算的表觀輻亮度78.396w/m2/sr/mic

表觀反射率0.3722（經(jīng)過大氣后反射率削弱）將此時的表觀輻亮度78.396代入（1）中，求得地表真實反射率為0.5002（幾乎等于真實的地表反射率0.5）可以得到大氣糾正公式：-------------（1）第二步：6.2單視角圖像的糾正方法6.2.1不變地物法invariant-objectmethods假設：有的像元反射率不隨時間變化；一幅影像上的像元大氣影響相同。方法，有清晰參考圖J，其他圖像利用基于“不變像元”的線性回歸來標準化為圖J。6.2單視角圖像的糾正方法6.2.1不變地物法invariant-objectmethods限制：每一波段都要找到從黑到白不同亮度的不變像元，否則回歸關系不準確；如果是對TOA輻亮度的標準化=/地表反射率的標準化，因為它們的關系不是線性的：6.2單視角圖像的糾正方法6.2.1不變地物法invariant-objectmethods優(yōu)點：簡單，正向；缺點：是統(tǒng)計方法，相對校正，很難糾正不均勻的氣溶膠散射。6.2單視角圖像的糾正方法6.2.1不變地物法invariant-objectmethods另一方法：查找表（1）根據(jù)氣溶膠厚度和水汽含量用輻射傳輸包（MOTRAN/6S）建查找表r，LTOA，（2）假設氣溶膠厚度和水汽含量來確定J圖中不變像元的反射率r（3）根據(jù)查找表和I圖中不變像元的頂層輻亮度（+第2步算出的r）確定I圖大氣參數(shù)；（4）根據(jù)大氣參數(shù)計算I圖其他像元的反射率。6.2單視角圖像的糾正方法6.2.2直方圖匹配histogrammatchingmethods假設：清晰和模糊區(qū)域的地表反射率直方圖相同（意味著兩區(qū)域的各類地物比例相同?!）6.2單視角圖像的糾正方法6.2.2直方圖匹配histogrammatchingmethods優(yōu)點：簡單，erdas軟件帶有；缺點：假設太武斷！6.2單視角圖像的糾正方法6.2.3暗目標法dark-objectmethod老暗目標法：假設有反射率為0的像元，不為0是由于大氣的影響，所以各波段像元值減去反射率本該為0的像元值；6.2單視角圖像的糾正方法6.2.3暗目標法dark-objectmethod新暗目標法：濃密植被為暗目標。步驟如下：

1、閾值法找到濃密植被像元；

2、根據(jù)藍光、紅光與2.1um中紅外的統(tǒng)計關系計算紅藍波段的反射率；

3、利用查找表確定藍、紅反射率的光學厚度；

4、由藍、紅反射率的光學厚度確定光學厚度函數(shù)，進而確定其他波段的光學厚度；

5、內(nèi)插得整圖的光學厚度，進而反演地表反射率。6.2單視角圖像的糾正方法6.2.4對比度降低法contrastreductionmethod假設：氣溶膠越多，TOA反射率（輻亮度）的差異越?。：┓瓷渎瘦椓炼冗b感得到TOA反射率差異，地表反射率差異可由清晰影像估算（需要假設地表反射率不變？），由此可計算大氣厚度，進而大氣糾正。適合灰塵和大霧6.2單視角圖像的糾正方法6.2.5類型匹配法clustermatchingmethod6.2單視角圖像的糾正方法Liang法的步驟：clustermatchingmethod（1）建立查找表：用輻射傳輸方程正向計算；如果地表郎伯，大氣水平均一：6.2單視角圖像的糾正方法Liang法的步驟：（2）估算氣溶膠厚度：假設同類地物的地表反射率一樣，無論是清晰還是模糊影像。6.2單視角圖像的糾正方法Liang法的步驟：（3）地表臨近效應：大氣和陸地之間大量多次散射造成，使得受鄰近像元影響，產(chǎn)生暗像元變亮，亮像元變暗效果?？紤]臨近效應的郎伯表面TOA反射率公式：總有效地表反射率：包含臨近效應地表到傳感器散射透過率6.2單視角圖像的糾正方法Liang法的步驟：（3）地表臨近效應：計算re按像元離散化，積分變加和：6.4估算整層大氣水汽含量的方法目的：校正近紅和中遠紅外影像，必須考慮水汽吸收。對于近紅和中紅外，只需整層大氣的水汽含量；對于熱紅外，需要水汽廓線。主要2類方法：遙感系統(tǒng)外獲??；利用遙感影像自身獲取。6.4估算整層大氣水汽含量的方法6.4.1概述外部方法：（1）地表觀測法：高時間分辨率，低空間分辨率；（2）基于地表的LiDAR法：向天空發(fā)射激光，測量回波量和時間，測量各高度的大氣密度。垂直分辨率高，白天時間分辨率高；（3）GPS法：由于水汽是GPS誤差的主要來源。如果信號接收位置已知，就可以根據(jù)檢測到的信號誤差計算水汽含量；6.4估算整層大氣水汽含量的方法6.4.1概述外部方法：（4）無線電探空儀：探空氣球，垂直分辨率高，水平分辨率低，不同時間分辨率；（5）飛機：機上搭載儀器。專業(yè)機費用昂貴，商用飛行項目尚未成熟；（6）紅外探測（infraredsoundingsensors）：基于特定分子的光譜判斷大氣組成。從6700-8200nm和13300-18600nm兩個通道得到垂直廓線，這兩個通道對水汽敏感；6.4估算整層大氣水汽含量的方法6.4.1概述外部方法：（7）微波探測（microwavesensors）：提供未結(jié)冰海洋整層大氣的水汽含量，垂直分辨率很差；（8）掩日法（solaroccultationmethods）：側(cè)向觀測大氣以及被動接收折射電磁波來確定大氣廓線。通過流體力學和平衡態(tài)方程確定大氣密度，壓力，溫度廓線；6.4估算整層大氣水汽含量的方法6.4.1概述外部方法：（9）航天GPS法：信號從GPS衛(wèi)星傳輸?shù)降蛙壭l(wèi)星（LEO），穿透大氣層時利用無線電波掩蓋進行主動探測。原理是大氣可透射、彎曲以及傳播信號，已知兩個衛(wèi)星的位置，就可以確定延時進而確定大氣廓線；6.4估算整層大氣水汽含量的方法6.4.2差分吸收法直接從多光譜和高光譜傳感器估算整層大氣的水汽含量；主要思路：在水汽吸收區(qū)域（如0.94um）取一波段，并在水汽吸收外取幾個波段，2類波段間的差異顯示水汽含量。很多算法：narrow/widealgorithm；continuuminterpolationbandratioalgorithms；curvefittingalgorithm；atmosphericprecorrecteddifferentialabsorptionalgorithm；6.4估算整層大氣水汽含量的方法6.4.2差分吸收法867nm1009nm6.4估算整層大氣水汽含量的方法6.4.2差分吸收法定義量R權(quán)重6.4估算整層大氣水汽含量的方法6.4.2差分吸收法建立水汽含量和R間關系：低估了暗目標表面的水汽含量6.4估算整層大氣水汽含量的方法6.4.2差分吸收法TahlandSchonermark(1998)修正：AfurthercorrectionisrequiredtoremovetheestimationbiasesbyV/f估算相對誤差小于20%6.4估算整層大氣水汽含量的方法6.4.2差分吸收法估算誤差3.1kg/m26.4估算整層大氣水汽含量的方法6.4.2差分吸收法估算誤差1.7kg/m26.4估算整層大氣水汽含量的方法6.4.3分裂窗算法如果沒有適用差分法的波段，用此替代；原理：兩個具有不同水汽吸收率的熱紅外波段（如11um，12um），這2個波段的亮溫影像協(xié)方差-方差可以估算整層大氣水汽含量。6.4估算整層大氣水汽含量的方法6.4.3分裂窗算法假設無云的大氣水平均一，地表溫度和反射率變化不大。那么兩波段影像上相鄰像元的亮度溫度差之比滿足：對于海面，可以分為發(fā)射率為1，那么R只與t有關，而t由于吸收氣體總量有關：6.4估算整層大氣水汽含量的方法6.4.3分裂窗算法假設其他氣體不變，整體透過率t的指數(shù)函數(shù)的一次近似是線性函數(shù)，則可以認為R與水汽總量線性相關。Jedlovec（1990）用亮度溫度空間方差擴展這一方法：6.4估算整層大氣水汽含量的方法6.4.3分裂窗算法再與水汽含量W建立統(tǒng)計關系：以上2個關系都用到了AVHRR圖像中。6.4估算