大气辐射校正6Sflaash实习课件

定量遥感实习定量遥感实习1-6S1-6S大气校正大气校正实习目的：实习目的：实习目的：实习目的：加深对加深对6S6S模型原理的理解，掌握模型原理的理解，掌握6S6S软件的使用方法与软件的使用方法与步骤，利用该软件进行步骤，利用该软件进行TMTM影像的大气校正。影像的大气校正。实习内容：实习内容：实习内容：实习内容：（1 1）熟悉各个步骤的物理含义。）熟悉各个步骤的物理含义。（2 2）利用）利用6S6S方法对方法对TMTM遥感影像的遥感影像的6 6个波段进行校个波段进行校正。正。实习步骤：实习步骤：实习步骤：实习步骤：见下页见下页实习思考：实习思考：实习思考：实习思考：尝试用暗像元校正方法对尝试用暗像元校正方法对TMTM图像进行校正，图像进行校正，并比较两种方法计算的结果。并比较两种方法计算的结果。定量遥感实习1-6S大气校正实习目的：1输入卫星的几何条件(Geometricalconditions):0-70-7中的数字：中的数字：0.use definition(0.use definition(用户自定义）用户自定义）1.meteosat observation1.meteosat observation2.goes east observation2.goes east observation3.goes west3.goes west4.avhrr(PM noaa)4.avhrr(PM noaa)5.avhrr(AM noaa)5.avhrr(AM noaa)6.hrv(spot)6.hrv(spot)7.tm(landsat)7.tm(landsat)1 输入卫星的几何条件(Geometrical condit2 2 根据输入对应时间，角度根据输入对应时间，角度Hour输入卫星过境时格林威治子午线的时间2 根据输入对应时间，角度Hour输入卫星过境时格林威治子午3 3 输入大气的模型输入大气的模型(Atmospheric Atmospheric model model)3 输入大气的模型(Atmospheric model)模型代号含义0无气体的吸收(no gaseous absorption)1热带的(Tropical)2中纬度夏季(midlatitude summer)3中纬度冬季(midlatitude winter)4亚北极夏季(subarctic summer)5亚北极(subarctic winter)6北美标准（us standard 62）3 3 输入大气的模型输入大气的模型(Atmospheric Atmospheric model model)模型代号含义0无气体的吸收(no gaseous absor4 4 输入气溶胶模型输入气溶胶模型（Aerosols modelAerosols model）大陆性的4 输入气溶胶模型（Aerosols model）大陆性的5 5 输入气溶胶模型的浓度输入气溶胶模型的浓度(可见度可见度)n n在集合-1,0,0的正值中选一个值，我们输入v=10km6 6 输入目标高度输入目标高度xps=0 xps=0 7 7 输入遥感器高度输入遥感器高度xpp=-1000 xpp=-10005 输入气溶胶模型的浓度(可见度)在集合-1,0,0的正8 8 输入传感器的波段代号输入传感器的波段代号8 输入传感器的波段代号9 9 输入地表状况输入地表状况n n0代表均匀n n1代表不均匀9 输入地表状况0代表均匀10 10 输入输入RAPP RAPP 的空气校正参数的空气校正参数在在-1至至0间得到的是间得到的是反射率反射率RAPP0时得到的是时得到的是亮度值，得到的此亮度亮度值，得到的此亮度值不能再根据反射率的值不能再根据反射率的公式转换为反射率，因公式转换为反射率，因为在为在6S模型中的算法模型中的算法额外加了一些参数，那额外加了一些参数，那个公式已不再适用个公式已不再适用10 输入RAPP 的空气校正参数在-1至0间得到的是反射率11 11 打开打开sixs.outsixs.out文件，找到经校正之后文件，找到经校正之后的系数，计算大气校正的表达式的系数，计算大气校正的表达式*atmospheric correction result *atmospheric correction result *-*-*input apparent reflectance :.309 *input apparent reflectance :.309 *measured radiance w/m2/sr/mic :90.000 *measured radiance w/m2/sr/mic :90.000 *atmospherically corrected reflectance:.367 *atmospherically corrected reflectance:.367 *coefficients xa xb xc :.00450 .02945 *coefficients xa xb xc :.00450 .02945 .08062 .08062 *y=xa*(measured radiance)-xb;acr=y/(1.+xc*y)*y=xa*(measured radiance)-xb;acr=y/(1.+xc*y)*y=0.00450*b1-0.02945y=0.00450*b1-0.02945Acr=Acr=（0.00450*b1-0.00450*b1-0.029450.02945）/(1+0.08062*/(1+0.08062*（0.00450*b1-0.029450.00450*b1-0.02945）)11 打开sixs.out文件，找到经校正之后的系数，计算大12 12 利用利用ENVIENVI软件进行大气校正软件进行大气校正FLAASHFLAASH12 利用ENVI软件进行大气校正FLAASH大气辐射校正6Sflaash实习课件大气辐射校正大气辐射校正一、辐射定标n nDN值图像向亮度值图像转换亮度值图像亮度单位*：10一、辐射定标DN值图像向亮度值图像转换亮度值图像10注意三点：l l增益（Gains）和偏置（Offset）取值头文件头文件USGSUSGS网站上定期公布数据网站上定期公布数据l l增益（Gains）和偏置（Offset）量纲和含义毫瓦毫瓦/（平方厘米（平方厘米*球面度）球面度）毫瓦毫瓦/（平方厘米（平方厘米*球面度球面度*微米）微米）注意三点：增益（Gains）和偏置（Offset）取值l l软件系统自带的增益和偏置 （通过练习体会）ENVIENVIPCIPCIERDASERDAS软件系统自带的增益和偏置二、辐射定标方法 1 1、ENVIENVI软件软件 1)1)由于由于ENVI 4.4 ENVI 4.4 中有专门进行辐射定标的模块，因此实际的中有专门进行辐射定标的模块，因此实际的操作十分简单。将原始操作十分简单。将原始TM TM 影像打开以后，选择影像打开以后，选择Basic ToolsPreprocessingCalibration UtilitiesBasic ToolsPreprocessingCalibration UtilitiesLandsat TMLandsat TM 2.)2.)进入下一步参数选择：根据传感器类型选择进入下一步参数选择：根据传感器类型选择Landsat 4,5 Landsat 4,5 或者或者7 7。从遥感影像的头文件中获取。从遥感影像的头文件中获取Data Acquisition Data Acquisition 的时间，的时间，Sun elevationSun elevation。如果你是用。如果你是用FileOpen External FileFileOpen External FileLandsatFast LandsatFast 的方法打开的方法打开header.dat header.dat 的话，的话，sun elevation sun elevation 就已经填好了。这里就已经填好了。这里Calibration Type Calibration Type 注意选择为注意选择为RadianceRadiance。输出文件，定标就完成了。输出文件，定标就完成了。2 2、ERDASERDAS 1)Model Maker 1)Model Maker 建立定标模型建立定标模型 2)2)注意生成的文件类型为浮点型注意生成的文件类型为浮点型二、辐射定标方法1、ENVI软件三、大气校正n n采用ENVI的FLAASH模块，基本操作步骤：三、大气校正采用ENVI的FLAASH模块，基本操作步骤：大气辐射校正6Sflaash实习课件大气辐射校正6Sflaash实习课件大气辐射校正6Sflaash实习课件