大气校正方法说明.docx
《大气校正方法说明.docx》由会员分享,可在线阅读,更多相关《大气校正方法说明.docx(28页珍藏版)》请在冰豆网上搜索。
![大气校正方法说明.docx](https://file1.bdocx.com/fileroot1/2023-6/7/764eb548-6b7f-49ea-9221-dc2e8e3e3989/764eb548-6b7f-49ea-9221-dc2e8e3e39891.gif)
大气校正方法说明
利用MODTRAN进行大气校正的方法说明
一.大气校正公式、原理以及所需参数
大气是介于传感器和地球表层之间由多种气体和气溶胶组成的介质层,电磁波在地物和传感器之间传输时,必然受到大气的影响。
遥感对地观测时,要想得到目标的真实信息,大气校正是不可回避的。
由卫星传感器获取的表观反射率
可由下式表出:
(1)
式中:
:
太阳天顶角,
:
太阳方位角,
:
传感器天顶角,
:
传感器方位角,
:
目标反射率,
:
大气的路径辐射项等效反射率,
:
大气的光学厚度,S:
大气的半球反照率,
:
散射透过率,
。
通过MODTRAN4对大气辐射传输进行模拟,求得大气校正所需参数,将所求的大气校正参数和传感器获得的表观反射率一并代入大气辐射传输公式
(1),便可计算出目标的真实反射率
,从而完成大气校正的任务。
在实际的工作中,我们可以用下面的公式:
(2)
式中:
是传感器接收到的辐射亮度,
是路径辐射项,
=
是太阳下行总辐射(
是大气层顶的太阳辐照度),
=
+
是传感器和目标之间的透过率(
是直射透过率,
是散射透过率)。
在已知的观测条件(太阳和传感器的几何参数,大气廓线,地表反射率等)下,设定一组
值以及相应的传感器高度,通过MODTRAN4模拟得到一组辐射亮度
,代入方程
(2),再经过简单的代数运算就可以求出大气校正所需的参数(路径辐射项、透过率、大气半球反照率和太阳下行总辐射)。
地表反射率和相应传感器高度设置见表1:
(地面高程时候传感器不受大气影响,L0项去掉;
=1表示完全透过)
模拟次数
1
2
3
4
5
地表反射率
0
0.1
0.2
0.1
0.2
传感器高度
100
100
100
地面高程
地面高程
表1地表反射率和相应的传感器高度参数设置
由
(2)式,可以解出
,
(3)
将传感器接收的辐射亮度和MODTRAN4模拟的大气校正参数代入方程(3)就可以进行大气校正。
二.MODTRAN介绍
MODTRAN——中等光谱分辨率大气透过率算法和计算模型,是由AFRL/VSBT和SpectralSciences,Inc.共同研发的一个产品。
MODTRAN以中等光谱分辨率,主要是2cm-1(紫外波段20cm-1),计算0到50000cm-1波段的大气透过率和辐射。
由于LOWTRAN系列波段模型算法的光谱分辨率和精度满足不了人们对更高分辨率和精度的需求,例如处理散射问题的二流近似限制了其计算精度,与方向有关的辐射亮度计算精度不高。
MODTRAN便应运而生。
除了LOWTRAN7分子波段模型参数化之外,MODTRAN包含了LOWTRAN7的所有功能,包括,球面折射的几何学、太阳和月亮的源函数、散射(瑞利散射,米散射,单次散射和多次散射)和默认的廓线(气体廓线,气溶胶廓线,云,雾,雨)。
目前发布的是MODTRAN4,3.1版本。
因为MODTRAN4是MODTRAN系列第一个申请专利、发布的,这个版本号意味着一些勘误表和物理模型的加入。
MODTRAN4主要的发展在于correlated-k算法的实施,大大方便了多次散射的准确计算。
这在本质上允许MODTRAN4进行“真实的比尔-朗伯”辐射传输计算,并且使衰减/层现在有了明确的物理意义。
更加准确的透过率和辐射计算极大促进高光谱影像数据的分析。
MODTRAN4提供了二向反射分布函数集,使物体表面散射逾越了朗伯体假设的局限。
correlated-K和二向反射分布函数集的结合极大改进了散射计算的精度。
三.如何利用MODTRAN得到大气校正所需参数
MODTRAN4的运行由当前目录下的tape5文件控制,tape5由一系列的CARDS(输入行)组成,输入参数由以下五个部分构成:
控制运行参数、传感器参数、大气参数、观测几何条件和地表参量。
tape5文件的主要作用是:
提供执行过程的参数(执行透过率的计算还是辐射的计算等)、提供计算所需的大气廓线(包括大气廓线的高度、压强、温度、水汽。
CO2、气溶胶等)、控制计算结果的输入(输入文件的名称,输入文件的内容以及格式等)。
由于辐射传输算法所设计的参数很多,所以tape5文件非常复杂,而且其中每个参数都有固定的格式和具体的含义,不能有丝毫的差别。
tape5文件必须严格的写入或由输入子程序生成。
设定这些参数后就可以用MODTRAN来模拟大气辐射传输过程,求解大气校正参数。
具体的参数设置和模拟结果由下面的例子case1和case2给出。
需要校正的是北京顺义地区ETM+图像,图像获取时间2001年4月17日,顺义的地理坐标为(
)。
case1:
使用的是默认的大气廓线(Mid-LatitudeSummer),采取的波段响应函数是TM7的波段相应函数(DATA目录下的tm7_filterfunction.txt)。
计算的是大气辐射。
输入文件为tape5,输出结果为channels.out文件。
tape5文件和channels.out文件的内容及说明见附录一。
表2列出了五次模拟的亮度值。
次数
波段
1
2
3
4
5
1
6.31E-06
7.65E-06
9.10E-06
2.01E-06
4.16E-06
2
4.05E-06
5.41E-06
6.84E-06
1.95E-06
4.00E-06
3
2.49E-06
3.80E-06
5.17E-06
1.76E-06
3.61E-06
4
1.00E-06
1.94E-06
2.92E-06
1.21E-06
2.47E-06
5
5.65E-08
3.29E-07
6.05E-07
3.11E-07
6.26E-07
6
9.22E-09
1.03E-07
1.97E-07
1.08E-07
2.16E-07
表2五次模拟的亮度值
把这五组辐射亮度值和设置的地表反射率、传感器高度对应起来,代入方程
(2),联立得到由五个方程组成的方程组:
(4)
将模拟的五次辐射亮度值
代入上面的方程组,不难解出大气校正参数。
(5)
计算结果如表3所示:
参数波段
路径辐射项
大气透过率
大气半球反照率
下行总辐射
1
6.31E-06
6.71E-01
3.30E-01
1.94E-05
2
4.05E-06
6.98E-01
2.78E-01
1.89E-05
3
2.49E-06
7.45E-01
2.32E-01
1.72E-05
4
1.00E-06
7.76E-01
1.68E-01
1.19E-05
5
5.65E-08
8.76E-01
6.02E-02
3.09E-06
6
9.22E-09
8.66E-01
3.31E-02
1.07E-06
表3大气校正参数
case2:
tape5文件写入的是自定义的大气廓线,由探空数据获取。
采取的波段响应函数是TM7的波段响应函数(DATA目录下的tm7_filterfunction.txt)。
计算的是大气辐射。
输入文件为tape5,输出文件为channels.out。
tape5文件和channels.out文件的内容及说明见附录二。
最后,按照表1的设置和已知数据,用MODTRAN4模拟得到五组传感器辐射亮度,见表4:
次数
波段
1
2
3
4
5
1
8.67E-06
1.17E-05
1.49E-05
4.28E-06
8.71E-06
2
5.85E-06
9.02E-06
1.24E-05
4.21E-06
8.54E-06
3
3.54E-06
6.45E-06
9.48E-06
3.66E-06
7.42E-06
4
1.47E-06
3.53E-06
5.66E-06
2.51E-06
5.07E-06
5
8.44E-08
6.26E-07
1.17E-06
5.94E-07
1.19E-06
6
1.70E-08
2.13E-07
4.10E-07
2.14E-07
4.29E-07
表4MODTRAN4模拟的五组传感器辐射亮度值
把这五组辐射亮度值和设置的地表反射率、传感器高度对应起来,代入方程
(2),联立得到由五个方程组成的方程组(4),利用方程组(5)解出大气校正所需参数,如表5所示:
参数波段
路径辐射项
大气透过率
大气半球反照率
下行总辐射
1
8.67E-06
7.22E-01
2.76E-01
4.07E-05
2
5.85E-06
7.68E-01
2.39E-01
4.03E-05
3
3.54E-06
8.06E-01
2.04E-01
3.53E-05
4
1.47E-06
8.32E-01
1.50E-01
2.44E-05
5
8.44E-08
9.14E-01
5.77E-02
5.89E-06
6
1.70E-08
9.18E-01
3.45E-02
2.13E-06
表5大气校正参数
附录一:
tape5文件格式:
TS22211111111000.0000.00
tf2f0365.0000011fttf0.000
DATA/BMP99_01.BIN
E:
\daqi\Mod4_Input\sample\DATA\tm7_filterfunction.txt
2101006.000000.000000.000000.000000.03300
100.0000.033180.0000.0000.0000.00000
121370
40.250243.5000.0000.0000.0000.0000.0000.000
3002500115rnn2aa
1
TS22211111111000.0000.10
tf2f0365.0000011fttf0.000
DATA/BMP99_01.BIN
E:
\daqi\Mod4_Input\sample\DATA\tm7_filterfunction.txt
2101006.000000.000000.000000.000000.03300
100.0000.033180.0000.0000.0000.00000
121370
40.250243.5000.0000.0000.0000.0000.0000.000
3002500115rnn2aa
1
TS22211111111000.0000.20
tf2f0365.0000011fttf0.000
DATA/BMP99_01.BIN
E:
\daqi\Mod4_Input\sample\DATA\tm7_filterfunction.txt
2101006.000000.000000.000000.000000.03300
100.0000.033180.0000.0000.0000.00000
121370
40.250243.5000.0000.0000.0000.0000.0000.000
3002500115rnn2aa
1
TS22211111111000.0000.10
tf2f0365.0000011fttf0.000
DATA/BMP99_01.BIN
E:
\daqi\Mod4_Input\sample\DATA\tm7_filterfunction.txt
2101006.000000.000000.000000.000000.03300
100.0000.033180.0000.0000.0000.00000
121370
40.250243.5000.0000.0000.0000.0000.0000.000
3002500115rnn2aa
1
TS22211111111000.0000.20
tf2f0365.0000011fttf0.000
DATA/BMP99_01.BIN
E:
\daqi\Mod4_Input\sample\DATA\tm7_filterfunction.txt
2101006.000000.000000.000000.000000.03300
0.0340.033180.0000.0000.0000.00000
121370
40.250243.5000.0000.0000.0000.0000.0000.000
3002500115rnn2aa
0
tape5文件说明:
第一行CARD1:
MODTRN,SPEED,MODEL,ITYPE,IEMSCT,IMULT,M1,M2,M3,M4,M5,M6,MDEF,IM,NOPRNT,TPTEMP,SURREF
格式(2A1,I3,12I5,F8.3A7)
Tape5文件的第一行和CARD1对应,表示采用的模型是MODTRANbandmodel, SPEED为slow(33kvalues),采用的大气模式是Mid-LatitudeSummer,计算的路径为斜程,程序以光谱热辐射和太阳/月亮辐射模式执行,计算过程中考虑多次散射,程序以常规模式运行,设定的地表反射率为0.00。
第二行CARD1A1:
DIS,DISAZM,NSTR,LSUN,ISUN,CO2MX,H2OSTR,O3STR,LSUNFL,LBMNAM,LFLTNM,H2OARE,SOLCON
FORMAT(2L1,I3,L1,I4,F10.5,2A10,4(1X,A1),2X,F10.3)
Tape5文件的第二行和CARD1A:
对应,采用DISORT二流近似离散纵标多次散射算法,CO2混合比为365ppmv,第三行从CARD1A2读取bandmodel参数,从后续的CARD1A3中读取用户定义的仪器波度响应函数。
第四行CARD1A3
FILTNM
格式(A80)
Tape5文件第四行和CARD1A3对应,读取用户自定义仪器波度响应函数。
第五行CARD2:
APLUS,IHAZE,CNOVAM,ISEASN,ARUSS,IVULCN,ICSTL,ICLD,IVSA,VIS,WSS,WHH,RAINRT,GNDALT
格式(A2,I3,A1,I4,A3,I2,3I5,5F10.5)
Tape5文件的第五行和CARD2相对应,消光类型为RURALextinction,视距6公里,对流层和平流层的季节气溶胶廓线为SPRING-SUMMER,地面高程为0.033公里。
第六行CARD3:
H1,H2,ANGLE,RANGE,BETA,RO,LENN,PHI,
格式(6F10.3,I5,5X,F10.3)
Tape5文件的第六行和CARD3对应,表示传感器的高度为100公里,目标高程为0.033公里,传感器的天顶角为180度(垂直向下观测)。
第七行CARD3A1:
IPARM,IPH,IDAY,ISOURC
格式(415)
Tape5文件的第七行和CARD3A1对应,观测时间在一年中是第137天,太空辐射源是太阳。
第八行CARD3A2:
PARM1,PARM2,PARM3,PARM4,TIME,PSIPO,ANGLEN,G
格式(8F10.3)
Tape5文件的第八行和CARD3A2对应,传感器的纬度40.25度,经度243.5度。
第九行CARD4:
V1,V2,DV,FWHM,YFLAG,XFLAG,DLIMIT,FLAGS
格式(4F10.0,2A1,A8,A7)
Tape5文件的第九行和CARD4对应,计算的是300nm到2500nm的大气辐射,步长为1nm,半高全宽度(FWHM)为15。
第十行CARD5:
IRPT
格式(I5)
Tape5文件的第十行和CARD5对应,0表示停止程序。
1表示读取后续的Tape5文件参数。
因此将多次模拟的参数写入同一个Tape5文件,就可以进行批处理。
输出文件channels.out格式:
CHANSPECTRALRADIANCECHANNELFULLCHANNELSPECTRALSPECTRALCHANNEL
NEL(WSR-1CM-2/XXXX)RADIANCEEQUIVALENTWIDTHMINIMUMMAXIMUMDESCRIPTION
NO.(PERCM-1)(PERNM)(WSR-1CM-2)(CM-1)(NM)(NM)(NM)
------------------------------------------------------------------------------------------
11.436700E-076.306379E-064.237350E-042949.363067.1915435.0000520.0000CENTER:
464NMFWHM:
44.00NM
21.267113E-074.048306E-063.141385E-042479.168577.5975500.0000624.0000CENTER:
560NMFWHM:
80.00NM
31.085361E-072.486430E-061.492019E-041374.675460.0065614.0000704.0000CENTER:
662NMFWHM:
62.00NM
46.950373E-081.003888E-061.212020E-041743.8197120.7326740.0000914.0000CENTER:
835NMFWHM:
126.00NM
51.531583E-085.646331E-081.075310E-05702.0906190.44401504.00001791.0000CENTER:
1648NMFWHM:
202.00NM
64.473961E-099.215902E-092.317407E-06517.9766251.45752001.00002386.0000CENTER:
2205NMFWHM:
282.00NM
CHANSPECTRALRADIANCECHANNELFULLCHANNELSPECTRALSPECTRALCHANNEL
NEL(WSR-1CM-2/XXXX)RADIANCEEQUIVALENTWIDTHMINIMUMMAXIMUMDESCRIPTION
NO.(PERCM-1)(PERNM)(WSR-1CM-2)(CM-1)(NM)(NM)(NM)
------------------------------------------------------------------------------------------
11.743712E-077.654005E-065.142840E-042949.363067.1915435.0000520.0000CENTER:
464NMFWHM:
44.00NM
21.691927E-075.405551E-064.194573E-042479.168577.5975500.0000624.0000CENTER:
560NMFWHM:
80.00NM
31.657733E-073.797664E-062.278845E-041374.675460.0065614.0000704.0000CENTER:
662NMFWHM:
62.00NM
41.346207E-071.944415E-062.347543E-041743.8197120.7326740.0000914.0000CENTER:
835NMFWHM:
126.00NM
58.931410E-083.292652E-076.270659E-05702.0906190.44401504.00001791.0000CENTER:
1648NMFWHM:
202.00NM
64.978463E-081.025512E-072.578728E-05517.9766251.45752001.00002386.0000CENTER:
2205NMFWHM:
282.00NM
CHANSPECTRALRADIANCECHANNELFULLCHANNELSPECTRALSPECTRALCHANNEL
NEL(WSR-1CM-2/XXXX)RADIANCEEQUIVALENTWIDTHMINIMUMMAXIMUMDESCRIPTION
NO.(PERCM-1)(PERNM)(WSR-1CM-2)(CM-1)(NM)(NM)(NM)
------------------------------------------------------------------------------------------
12.072413E-079.096831E-066.112297E-042949.363067.1915435.0000520.0000CENTER:
464NMFWHM:
44.00NM
22.141794E-076.842832E-065.309868E-042479.168577.5975500.0000624.0000CENTER:
560NMFWHM:
80.00NM
32.257958E-075.172707E-063.103960E-041374.675460.0065614.0000704.0000CENTER:
662NMFWHM:
62.00NM
42.020057E-072.917701E-063.522616E-041743.8197120.7326740.0