辐射定标和大气校正过程参考
实验数据来源:
使用的数据为广东省汕头市的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-5 2.在ENVI中操作如图:
定标前: 定标后:
二. 大气校正
1.将图像的辐亮度转化为表现反射率
???*L*d2/(ESUN*cos(?))
其中ρ为表观反射率,L为表观辐亮度,d为日地距离,ESUN为太阳平均辐射强度,θ为太阳天顶角。ESUN的值从表3中查得。d的值根据影像成像的儒略日(在一年中所在天数)从表4查得,如实习影像成像时间是2001年11月22日,儒略日为第326天,d=0.9860天文单位。θ从头文件中读取为41.36°,cosθ=0.7506,
2
表观反射率计算公式为: ρ=3.142*L*(0.9860)/(1554*0.7506)。 参考表格:
2.在ENVI中操作如图:
结果图:
三. 基于6s模型的大气纠正 1.输入文件:input3
2.通过cmd.exe执行下列操作得到output3.txt文件
3.找到所需数据
由output3.txt可知coefficients xa xb xa : 0.00543 0.02145 0.05637。 4. 利用公式计算校正后的反射率
其中,ρ为校正后的反射率,Li是i波段的辐射量度 得到计算公式为:y=0.00543*Li-0.02145
5.利用ENVI计算用6s得到模型进行的大气纠正
四.
对比
大气纠正完得到的是地表真实反射率,而辐射定标完得到的是表观发射率,二者的区别就是表观反射率经过大气校正之后得到的才是真实反射率,所以两个的值有所差别