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

基于RS\GIS监测洪灾变化上机操作实例 基本原理： 大气校正 遥感图像在获取过程中，受到大气吸收与散射、传感器定标、地形等因素的影响，且会随时间的不同而有所差异。利用多时相遥感图像的光谱信息检测地物变化重要前提是要消除不变地物的辐射值差异。 大气校正的目的是消除大气和光照等因素对地物反射的影响，大多数情况下，大气校正是反演地物真实反射率的过程。ρ）。 计算表观辐亮度的公式： adiance=（（max-Lmin）/（calmax-Qcalmin）*（cal-Qcalmin）+min ① 其中：Radiance 是表观辐亮度W/m2·sr·μm；Qcal为像元DN值（也就是影像数据本身）；calmax为传感器处最大辐亮度值所对应的DN值，一般为255；Qcalmin ，一般为0max 和min是从参数表中查询min为光谱辐亮度的最小值，单位同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来说，calmax取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 * μm) 1 1969.000 2 1840.000 3 1551.000 4 1044.000 5 225.700 7 82.07 8 1368.000 Earth-Sun Distance in Astronomical Units Julian Day Distance Julian Day Distance Julian Day Distance Julian Day Distance Julian Day Distance 1 .9832 74 .9945 152 1.0140 227 1.0128 305 .9925 15 .9836 91 .9993 166 1.0158 242 1.0092 319 .9892 32 .9853 106 1.0033 182 1.0167 258 1.0057 335 .9860 46 .9878 121 1.0076 196 1.0165 274 1.0011 349 .9843 60 .9909 135 1.0109 213 1.