操作-大气校正,辐射定标,气溶胶反演.doc

基于RS\GIS监测洪灾变化上机操作实例 基本原理： 大气校正 遥感图像在获取过程中，受到大气吸收与散射、传感器定标、地形等因素的影响，且会随时间的不同而有所差异。利用多时相遥感图像的光谱信息检测地物变化的重要前提是要消除不变地物的辐射值差异。 大气校正的目的是消除大气和光照等因素对地物反射的影响，大多数情况下，大气校正是反演地物真实反射率的过程。 目前可以进行大气校正的模块有很多种，如最早的 MODTRAN 4+，6S (Second Simulation of the Satellite Signal in the Solar Spectrum)，ACORN，ATREM，在 ERDAS IMAGINE 8.7上的模块 ATCOR，以及 ENVI上的模块 FLAASH（基于 MODTRAN）。 FLAASH可对LANDSAT，SPOT，AVHRR，ASTER，MODIS，MERIS，AATSR， IRS等多光谱、高光谱数据、航空影像及自定义格式的高光谱影像进行快速大气校正分析。下面的大气纠正步骤，都是基于FLAASH进行的。 辐射定标 当我们拿到一幅原始影像，先要进行辐射定标，目的是把图像上的DN（Digital Number）值转为辐亮度或者是反射率。辐射定标的结果可以是表观辐亮度（L），也可以是表观反射率（ρ）。 计算表观辐亮度（L）的公式为： Radiance=（（Lmax-Lmin）/（Qcalmax-Qcalmin）*（Qcal-Qcalmin）+Lmin ① 其中：Radiance 是表观辐亮度，注意单位是W/m2·sr·μm；Qcal为像元DN值（也就是影像数据本身）； Qcalmax为传感器处最大辐亮度值所对应的DN值，一般为255；Qcalmin 为传感器处最大辐亮度值所对应的DN值，一般为0；Lmax 和Lmin是从参数表中查询，Lmin为光谱辐亮度的最小值，单位同L；Lmax为光谱辐亮度的最大值，单位同L。 计算表观反射率（ρ）的公式为： ρ =π*L*d2/（ESUN*cos（θ）） ② 其中：ρ为表观反射率；L为①式中计算出来的表观辐亮度；d为日地距离；ESUN为大气层外的太阳辐射，也可以说是传感器接收处的太阳辐射；θ为太阳天顶角（这个可以通过影像的元数据获取）。以上参数可以查询下表获得。 对于TM影像来说，当LANDSAT传感器获取影像以后（Level 0），会将其转化为 32 位浮点型的绝对辐亮度。之后进一步处理，将绝对辐亮度变为8位的DN值，这也就是我们购买后拿到的数据（Level 1）。 如果要将 L1的DN值转化为传感器处的辐亮度值（at-sensor spectral radiance），对于TM是8bit来说，Qcalmax取255，Qcalmin取0，可以将公式①简化为下面这个公式： 上面的这个公式也可以改为： 其中， Grecale即为gain，Brescale即为offset。各个波段的，以及和见下表。 需要注意的是，上述参数在2003年5月5日前后是不一致的，所以在操作时，一定要搞清楚影像获取的时间。ETM+影像与TM影像的Lmin、Lmax、Grescale、Brescale、d及ESUN的参数不一样，可以查询下列表格： ETM+ Spectral Radiance Rangewatts/(meter squared * ster * μm) Band Number Before July 1, 2000 After July 1, 2000 Low Gain High Gain Low Gain High Gain LMIN LMAX LMIN LMAX LMIN LMAX LMIN LMAX 1 -6.2 297.5 -6.2 194.3 -6.2 293.7 -6.2 191.6 2 -6.0 303.4 -6.0 202.4 -6.4 300.9 -6.4 196.5 3 -4.5 235.5 -4.5 158.6 -5.0 234.4 -5.0 152.9 4 -4.5 235.0 -4.5 157.5 -5.1 241.1 -5.1 157.4 5 -1.0 47.70 -1.0 31.76 -1.0 47.57 -1.0 31.06 6 0.0 17.04 3.2 12.65 0.0 17.04 3.2 12.65 7 -0.35 16.60 -0.35 10.932 -0.35 16.54 -0.35 10.80 8 -5.0 244.00 -5.0 158.40 -4.7 243.1 -4.7 158.3 ETM+ Solar Spectral Irradiances Band watts/(meter squared