《应用6S模型进行LANDSAT TM影像大气校正1》.doc

应用6S模型进行LANDSAT TM影像大气校正 ? 一、辐射校正 1、? 用定标系数将原始DN值转换为大气层顶太阳辐亮度L； L为大气层顶太阳辐亮度，Q为记录的电信号数值，，为通道增益，为偏移量，定标系数可以在头文件中获得。 表1 LANDSAT5 TM数据定标系数 Processing Data L-5 TM Post-Calibration Dyanamic Ranges for U.S.Processed NLAPS Data Processing Date From March 1st 1984 To May 4st2003 After May 5st 2003 Band Grescale Brescale Grescale Brescale 1 0.602431 -1.52 0.762824 -1.52 2 1.175100 -2.84 1.442510 -2.84 3 0.805765 -1.17 1.03988 -1.17 4 0.814549 -1.51 0.872588 -1.51 5 0.108078 -0.37 0.119882 -0.37 6 0.055158 1.2378 0.055158 1.2378 7 0.056980 -0.15 0.065294 -0.15 ? 2、由大气层顶太阳辐亮度L转换为反射率。 其中: : 行星反射率 : 传感器口径的光谱辐射值 d: 日地距离（以天文为单位） :Mean solar exoatmospheric irradiances平均太阳外大气层辐射值 : 太阳天顶角 ? 表2 TM 太阳外大气层光谱辐射值 单位 ESUN=W/（m2.um） 模型 Neckel and Labs Chance Spectrum CHKUR 波段 Landsat 4 Landsat 5 Landsat 4 Landsat 5 1 1958 1957 1957 1957 2 1828 1829 1825 1826 3 1559 1557 1557 1554 4 1045 1047 1033 1036 5 219.1 219.3 214.9 215 7 74.57 74.52 80.72 80.67 ? 表3,日地距离（以天文为单位） 天 距离 天 距离 天 距离 1 0.9832 121 1.0076 242 1.0092 15 0.9836 135 1.0109 258 1.0057 32 0.9853 152 1.014 274 1.0011 46 0.9878 166 1.0158 288 0.9972 60 0.9909 182 1.0167 305 0.9925 74 0.9945 196 1.0165 319 0.9892 91 0.9993 213 1.0149 335 0.986 106 1.0033 227 1.0128 349 0.9843 ? ? ? ? 365 0.9833 ? 二、大气校正 经过辐射校正后，象元灰度值转换为了反射率，我们使用6S模型对可见光和近红外波段进行大气校正。 1、辐射校正完的反射率是0-1之间的值，然后把它转换为0-100之间的数值； 2、再在ENVI中把第一步中的反射率（0-100）存为RAW格式； 3、然后在inputfiles中填写大气条件的输入文件、大气条件的输出文件名、待大气校正的输入文件（RAW格式）、待大气校正的图像大小； 4、? 第3步中的大气条件的输入文件需要填写以下的几项： ????????? Landsat5 geometrical conditions ????????? month,day,hh.ddd,long.,lat. ????????? tropical atmospheric mode ????????? continental aerosols model ????????? visibility in km (aerosol model concentration) ????????? target at 600 m above sea level ????????? sensor on board of satellite ????????? third band of Landsat5 ????????? the image has values of reflectance, DN is percent (actual values only 0-100, not 0-255)(-1) ????????? number of pixels of the image=number of bytes 以上的这些参数可以根据实际