ENVI--专题五辐射定标与大气纠正.docVIP

专题五 辐射定标与大气纠正 图像处理流程分为以下几个步骤：图像的配准、重采样、定标以及大气纠正。其中，图像配准是做所有工作的前提，是图像的几何纠正。在进行定性或者定量遥感时都要进行图像配准来确定我们所要研究的目标。在进行完图像配准（几何纠正）之后，为了使得输出图像的配置与输入图像向对应，因此要进行重采样。定标以及大气纠正则是进行辐射量纠正的重要的过程，是进行定量遥感不可少的步骤。分类是图像处理的最后一步，按某种使用意图分类之后的图像，对实际应用有很大帮助。 辐射定标和大气纠正都属于图像的预处理，辐射定标的目的是把图像上的DN值转为辐亮度或者是反射率，大气纠正的目的是消除或减少大气对图像的干扰。 以下是图像预处理的流程： 辐射定标 要对图象进行辐射定标，将图象的DN值转化为表观反射率，该过程的实现是通过应用以下两个步骤来实现的： 首先将图象的DN值转化为辐亮度：radiance=gain*DN+offset (1)式 然后将图象的辐亮度转化为表观反射率：ρ=π*L*d2/（ESUN*cos（θ）） (2)式 其中ρ为表观反射率，L为表观辐亮度，d为日地距离，ESUN为太阳平均辐射强度，θ为太阳天顶角。 将以上两个步骤结合得：ρ=π*（gain*DN+offset）* d2/（ESUN*cos（θ））(3)式 ENVI中的具体实现：采用简单的波段运算 例如，我们把2002-5-22的一幅ETM图像第3波段的DN值转化为表观反射率。 第一步，查找.FST文件（是该图像经过粗纠正后的数据信息说明文件），该图像第3波段的gain和offset，可知： gain＝0.619215662339154，offset＝-5.0000000000000 第二步，查找.FST文件可知=90-62.7=27.3度，cos（θ）＝0.8886；查找Table 11.4 （Earth-Sun Distance in Astronomical Units）可知 d=1.10109天文单位；查找Table 11.3 ? （ETM+ Solar Spectral Irradiances）可知ESUN＝1551。 第三步，把这些参数的值带入（1）式：L3＝0.619*b3-5.0 代入（3）式：ρ3＝3.1415* (0.619*b3-5.0)*1.10109^2/(1551*0.4602) 运算前： 实施运算 Basic tools band math 输入运算式: 运算之后: 对于其它波段只要知道相关的参数,可以用同样的方法作简单的定标. 其中关于ETM图像的一些参数如下: Table 11.2 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 Table 11.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 Table 11.4 Earth-Sun Distance in Astronomical Units Julian Day Distance Julian Day Distance Julian Day Distance Julian Day