辐射定标和大气校正过程参考

辐射定标和大气校正过程参考实验数据来源：使用的数据为广东省汕头市的ETM+影像，成像时间为2001年11月22日，2:28:18.000（格林威治时间）。数据处理一．辐射定标1.首先对图像进行辐射定标，将图像的DN值转化为辐亮度。每个角标中含有的参数表示波段不同则取值不同，具体参数可从卫星影像的头文件中得到。L是某个波段光谱辐射亮度；gain为增量校正系数，offset为校正偏差量，DN是图像灰度值，DNmax和DNmin为遥感器最大和最小灰度值，Lmax，Lmin分别为最大和最小灰度值所相应的辐射亮度。Band3:定标公式：L=(152.9+5)/(255-1)*b1-52.在ENVI中操作如图：定标前：定标后：二．大气校正1.将图像的辐亮度转化为表现反射率))cos(*/(**2ESUNdL其中ρ为表观反射率，L为表观辐亮度，d为日地距离，ESUN为太阳平均辐射强度，θ为太阳天顶角。ESUN的值从表3中查得。d的值根据影像成像的儒略日（在一年中所在天数）从表4查得，如实习影像成像时间是2001年11月22日，儒略日为第326天，d=0.9860天文单位。θ从头文件中读取为41.36°，cosθ=0.7506，表观反射率计算公式为：ρ=3.142*L*(0.9860)2/(1554*0.7506)。参考表格：2.在ENVI中操作如图：结果图：三．基于6s模型的大气纠正1.输入文件：input32.通过cmd.exe执行下列操作得到output3.txt文件3.找到所需数据由output3.txt可知coefficientsxaxbxa:0.005430.021450.05637。4.利用公式计算校正后的反射率其中，ρ为校正后的反射率，Li是i波段的辐射量度得到计算公式为：y=0.00543*Li-0.021455.利用ENVI计算用6s得到模型进行的大气纠正四．对比大气纠正完得到的是地表真实反射率，而辐射定标完得到的是表观发射率，二者的区别就是表观反射率经过大气校正之后得到的才是真实反射率，所以两个的值有所差别