Flaash大气校正.docx

上机实习内容：Flaash大气校正学生姓名 王 玲 学 号 201420771院 系 城市与环境学院 专 业 地图学与地理信息系统 年 级 2014级 教务处制Flaash大气校正实验报告实验目的通过本次实验能够更深一步理解大气校正的原理、方法。并且熟练掌握Landsat8 OLI数据的大气校正的流程。实验内容辐射定标目的：将传感器记录的电压或数字量化值（DN值）转换为绝对辐射亮度值（辐射率）。原理： L=Gain*DN + Bias步骤：首先，在Envi5.1中打开辐射定标工具，Toolbox/Radiometric Correction/ Radiometric Calibration，并在File Selection对话框中选择数据，如下所示：辐射定标参数设置当选择好辐射定标的数据时，接下来需选择定标参数。其中，①Calibration Type：辐射定标类型，因Flaash校正要求输入的数据为辐亮度值，因此辐射定标类型选择辐亮度。当数据的每个波段包含Gain和Offest参数时，Envi会自动从元数据文件中获取这些参数，并按照辐射定标公式进行定标，本实验所使用的Landsat8 OLI数据的元数据中包含这两个参数。另外，Envi默认Gain和Offest参数定标单位为W/（m2*sr*μm），因此，计算得到的辐亮度值为W/（m2*sr*μm）。②Output Interleave：输出数据存储顺序，因Flaash校正要求输入的数据存储类型为BIL或BIP，但因BIL的处理速度快，故在此选择BIL。③Output Data Type：输出数据类型，辐射定标中可以选择的输出数据类型为三种，分别是：浮点型（Float）、双精度浮点型(Double)和无符号位16整型(Uint)。本实验中使用的OLI6原始数据为无符号16位整型，在进行Flaash校正时计算缩放因子是无单位型与浮点型数据之间的缩放关系，因此，该处选择浮点型（Float）。④Scale Factor：因辐射定标计算的辐亮度值单位是W/（m2*sr*μm），而FLAASH校正所要求输入数据的辐亮度单位为μW/（cm2*sr*nm），该缩放系数是这两单位间的转换系数，而点击“Apply FLAASH Settings”按钮，会自动匹配符合FLAASH大气校正要求的缩放系数。（3）辐射定标后的结果说明：上图中左边的视图框中显示的是辐射定标前的多光谱影像，也即原始的OLI多光谱影像；右边的视图框中显示的是辐射定标后的多光谱影像。结果检验查看定标前后影像的像元值，结果如下：其中，底下的Data值为定标前影像的像元值，上面的Data值为定标后影像的像元值。且定标前影像的数据类型为无符号16位整型，定标后影像为浮点型辐亮度值。Flaash大气校正目的：消除大气和光照等因素对地物反射的影响，获取地物反射率、辐射率、地表温度等真实物理模型参数，用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧等对地物反射的影响，消除大气分子和气溶胶散射的影响。原理：是基于太阳波谱范围内（不包括热辐射）和平面朗伯体（或近似平面朗伯体），在传感器处接收的像元光谱辐射亮度公式为：其中，公式右边第一部分代表了太阳辐射经大气入射到地表后又反射直接进入传感器的辐射亮度；右边第二部分为经大气散射后进入传感器的辐射亮度；第三部分为大气后向散射辐射率（大气程辐射）。进行大气校正时，在大气校正参数获取后, 根据该方程逐个波段像元计算地表平均反射率。步骤：首先，在Envi中打开Flaash大气校正工具，Toolbox/Radiometric Correction/FLAASH Atmospheric Correction，打开的Flaash校正界面如下：Input Radiance Image：输入辐亮度文件，单击该按钮选择辐射定标后的文件，这时会弹出Radiance Scale Factor对话框，它的作用是在单位未知的辐亮度值与浮点型辐亮度值之间做转换，并且Scale Factor=浮点型或整型亮度图像（单位未知）/ 浮点型辐射亮度。当各波段的辐亮度单位不一致时，选择第一项，否则选择第二项，因做辐射定标时已经使缩放系数符合Flaash校正要求，因此此处的缩放系数为默认值1。参数设置①传感器与图像信息Scene Center Location：图像中心经纬度，该值在选择了输入辐亮度文件后自动加载，无需手动填写；Sensor Type：选择多光谱里的Landsat8 OLI传感器；Sensor Altitude：传感器飞行高度，当选择了传感器类型，该值将会自动添加；Pixel Size：像素大小，当选择了传感器类型，该值将会自动添加；Flight Data与Flight Time GMT：成像日期与时间，在相