应用6S模型进行LANDSAT

1、应用6S模型进行LANDSAT TM影像大气校正一、辐射校正1、1、用定标系数将原始DN值转换为大气层顶太阳辐亮度L；人 rescale cal rescaleL为大气层顶太阳辐亮度，Q为记录的电信号数值，G爬坤，为通道增益，Bree为偏移 量，定标系数可以在头文件中获得。表1 LANDSAT5 TM数据定标系数Processing Data L-5 TM Post-Calibration Dyanamic Ranges for U.S.Processed NLAPS DataProcessing DateFrom March 1st 1984 To May 4st2003After May

2、5st 2003BandG , rescaleB , rescaleG , rescaleB , rescale10.602431-1.520.762824-1.5221.175100-2.841.442510-2.8430.805765-1.171.03988-1.1740.814549-1.510.872588-1.5150.108078-0.370.119882-0.3760.0551581.23780.0551581.237870.056980-0.150.065294-0.152、由大气层顶太阳辐亮度L转换为反射率。L -d2 兀ESUN 了 cos 0其中：P p:行星反射率L人:

3、传感器口径的光谱辐射值d:日地距离（以天文为单位）ESUN人:Mean solar exoatmospheric irradiances平均太阳外大气层辐射值0:太阳天顶角表2 TM太阳外大气层光谱辐射值单位ESUN=W/（m2.um）模型Neckel and LabsChance Spectrum CHKUR波段Landsat 4Landsat 5Landsat 4Landsat 5119581957195719572182818291825182631559155715571554410451047103310365219.1219.3214.9215774.5774.5280.7280.

4、67表3,日地距离（以天文为单位 ）天距离天距离天距离10.98321211.00762421.0092150.98361351.01092581.0057320.98531521.0142741.0011460.98781661.01582880.9972600.99091821.01673050.9925740.99451961.01653190.9892910.99932131.01493350.9861061.00332271.01283490.98433650.9833二、大气校正经过辐射校正后，象元灰度值转换为了反射率，我们使用6S模型对可见光和近红外波段 进行大气校正。1、辐射校

5、正完的反射率是0-1之间的值，然后把它转换为0-100之间的数值；2、再在ENVI中把第一步中的反射率（0-100）存为RAW格式；3、然后在inputfiles中填写大气条件的输入文件、大气条件的输出文件名、待大气校正的输入 文件（RAW格式）、待大气校正的图像大小；4、4、第3步中的大气条件的输入文件需要填写以下的几项： Landsat5 geometrical conditionsmonth,day,hh.ddd,long.,lat.tropical atmospheric modecontinental aerosols modelvisibility in km （aerosol m

6、odel concentration）target at 600 m above sea levelsensor on board of satellitethird 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以上的这些参数可以根据实际情况进行填写。三、举例说明以1987年10月26日2： 35： 52时、太阳高度角为40.413

7、4、经纬度分别为：111.173, 31.7463、的TM （第三波段）遥感影像为例进行说明：1、1、用定标系数将原始DN值转换为大气层顶太阳辐亮度L；L=DN*1.03988-1.172、2、由大气层顶太阳辐亮度L转换为反射率；L 0.98922 冗P =p 1554 cos49.58663、 3、把第2步计算出的反射率先乘以100，然后取整，并且在ENVI中把它存为RAW格 式；4、4、运用6S模型进行大气校正；inputfiles.txt:TM3_inputfiles.txtTM3.resb3.rawb3raw.raw(atmospheric conditions input text

8、file)(atmospheric conditions output text file)(input image file, in this case 400 x 400 pixels)(outpu timage file, 400 x 400 pixels)TM3_inputfiles.txt710 26 2.583574 111.173 31.7463 universal time)215-600-100027-1100,not 0 and 255)55382784 (7796*7104)(Landsat5 geometrical conditions)(month,day,hh.ddd,long.,lat.) (hh.ddd=the decimal hour in (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)(TM image wit