MODTRAN使用要点精解.docx
《MODTRAN使用要点精解.docx》由会员分享,可在线阅读,更多相关《MODTRAN使用要点精解.docx(20页珍藏版)》请在冰豆网上搜索。
MODTRAN使用要点精解
MOD37使用要点
1MODTRAN文件
1.1输出文件
输入Edit|EditFile可查看运行Modtran产生的ASCII数据文件,在usr目录
下有:
Modout1、Modout2、Modout3
Modout1包含完整结果,波段平均透过率、辐亮度在最后一段。
INTEGRATEDABSORPTIONFROM625TO25000CM-1=7170.3120CM-1
AVERAGETRANSMITTANCE=0.7058
INTEGRATEDTOTALRADIANCE=2.602227E-03WATTSCM-2STER-1(FROM
625TO25000CM-1)
MINIMUMSPECTRALRADIANCE=0.000000E+00WATTSCM-2STER-1/CM-1
AT25000CM-1
MAXIMUMSPECTRALRADIANCE=1.168842E-05WATTSCM-2STER-1/CM-1
AT640CM-1
BOUNDARYTEMPERATURE=0.000
BOUNDARYEMISSIVITY=1.000
1.2输入文件
1.2.1数据库文件
文件名.LTN,包含Modtran模型参数的输入文件和输出曲线打印方式的输入
文件。
1
1.2.2NOVAM参数文件:
文件名.NPR,6000米以下海面气溶胶分布,较以前海洋气溶胶模型有明显
改进。
1.2.3滤波器函数文件:
文件名.FLT
滤波器函数可用于计算有效大气透过率或有效路径辐亮度。
滤波器函数包含
辐射源的光谱特性和辐射测量仪器的光谱带通特性。
辐射源的光谱特性用黑体温
度表示,仪器光谱特性取决于探测器光谱响应和光学元件光谱透过或反射特性,
Modtran可设置最多80个波数点的相对响应。
1.2.4扫描函数文件:
文件名.SCN
Modtran内部计算样点的光谱间隔为1cm-1,为降低输出数据的光谱分辨率,
Modtran提供了扫描功能,可将输出光谱数据与一个狭缝函数进行卷积运算。
狭
缝函数可以是三角形、矩形、Sinc.、Sinc.平方等函数。
FWHM是狭缝函数的半宽度,最小设置值2cm-1,最大设置值50cm-1。
因此,
Modtran输出数据光谱分辨率的极限为2cm-1。
输入的频率增量(步长)DV是指输出光谱数据采样间隔,用波数表示,并
非Modtran的分辨率。
如步长DV设置为100cm-1,FWHM设置为2cm-1,Modtran
输出数据的波数间隔为100cm-1,数据点的光谱分辨率为2cm-1,波数位于100
cm-1采样间隔区的大气强吸收峰或辐射峰均不予计算。
2MODTRAN概述
2.1概述
MODTRAN(中分辨率透过率的英文缩写)是美国空军研究实验室开发的
2
LOWTRAN升级产品。
实际上,MODTRAN包含了可供用户选用的LOWTRAN
完整模型。
MODTRAN可计算特定大气路径的透过率和辐亮度。
它既可作为独立的程序
运行,可可作为子程序或分立的模块运行。
MODTRAN源程序用FORTRAN77
编写,可向终端用户提供源代码。
MODTRAN计算
1)波数范围
0-50000cm-1,即最小波长0.2um。
按1cm-1波数间隔分割,即计算步长
为1cm-1。
2)大气分层
大气高度0-100km,根据气压、温度、分子吸收与消光系数对大气层进行
分层,最大可分33层。
MODTRAN与LOWTRAN的比较:
1)LOWTRAN计算的光谱范围与MODTRAN相同,但计算步长为5cm-1,
分辨率为20cm-1。
2)LOWTRAN的近似分子吸收带模型仅有吸收系数、分子密度刻度函数等
参数。
要用作工作于下层大气和地表面战术系统大气传输计算的辅助工
具,40公里以上LOWTRAN7计算精度大大降低。
3)MODTRAN的分子吸收带模型有吸收系数、谱线密度、平均谱线密度等
3个温度有关参数。
MODTRAN在对1cm-1波数间隔内吸收谱线积分时考
虑了谱线的形状。
MODTRAN适合大气路径高度30公里以上的透过率计
算,但不适合大气路径高度60公里以上的透过率计算,因为许多分子不
处于局部的热力学平衡状态,不能根据周围的温度确定分子带的辐射。
MODTRAN大气透过计算包括:
1)气体分子谱线吸收;
3
2)气体分子连续吸收;
3)气体分子散射;
4)气溶胶吸收和散射;
MODTRAN路径辐亮度计算包括:
1)大气自身辐射;
2)以单次散射方式进入路径的太阳和/或月亮辐亮度;
3)斜路径对天观察时的直接太阳辐照度;
4)以多次散射方式进入路径的太阳和/或月亮辐亮度以及地球辐亮度。
2.2大气路径类型与参数
MODTRAN计算时,要求先输入路径类型,再输入路径参数。
2.2.1水平路径
确定水平路径只需输入2个参数,即观察点高度、路径长度。
这里的水
平路径系短路径,不是长路径。
4
2.2.2斜路径
可选用的斜路径参数包括:
H1:
初始高度(观察点高度);
H2:
最终高度(辐射源高度);
HMIN:
正切高度,仅用于长程。
正切高度是路径至地球的最低高度。
r:
起点与终点之间的直线距离,由于大气折射,长路径时,实际路
径的长度大于此值。
:
天顶角,观察点的重垂线与路径的夹角;
:
最终角,源点的重垂线与路径的夹角;
:
地心对观察点、源点的张角。
确定一个斜路径只需输入3个参数,共有6种设置方法:
1)H1,,H2;
2)H1,,r;
3)H1,H2,r;
4)H2,H1,;
5)H2,H1,;
6)H2,,r。
5
2.2.3斜路径至太空
确定一个至太空的斜路径需要输入2个参数,共有3种设置方法:
1)H1,;
2)H1,HMIN;
3)HMIN,;
其中第2种方法主要用于临边探测(LimbViewing)。
2.3气溶胶模型
气溶胶是大气中悬浮颗粒的总称,它包括地面灰尘、火山灰、工业燃烧产生
烟灰、海面喷沫、雾等。
气溶胶颗粒大小变化很大,仅存在于局部区域。
气溶胶分布与海拔高度有关,MODTRAN按高度将大气分为4层,分别建立
气溶胶模型。
1)边界层(0-2km);
2)对流层顶部(2-10km);
3)同流层底部(10-30km);
4)同流层(中间层)(30-100km)。
边界层气溶胶模型与地理环境、天气一个,MODTRAN给出的边界层气溶胶
模型主要有:
农村、都市、海洋、对流层、雾等。
6
农村的气溶胶模型由占比70%的可溶物质(氨、硫酸钙及有机化合物)和占
比30%类灰尘气溶胶混合而成。
海洋的气溶胶模型由盐颗粒和背景两部分组成,
盐颗粒是海水飞沫蒸发后再凝聚核水汽形成的较大颗粒,存在于离海面10-20
米高处。
海洋的背景气溶胶与陆地农村气溶胶相仿,唯一不同是没有非常大的颗
粒。
都市气溶胶模型由占比20%的燃烧生成物或工业源形成的类烟尘气溶胶,其
余80%为农村气溶胶。
对流层气溶胶模型代表了一种非常清彻天气条件,能见度
达50公里。
对流层顶部(2-10km)的气溶胶比边界层均匀的多,大颗粒迅速减少。
同流层
底部(10-30km)的气溶胶主要受季节影响,因为对流层顶部的高度是随季节变化
的变化。
该层呈全球均匀分布,不受地理环境影响,主要颗粒是光化学反应产生
的硫酸盐颗粒和火山喷发时的火山灰。
高层大气的主要气溶胶是流星灰。
2.4大气模型
MODTRAN定义了6种标准的大气模型,也允许用户输入气象数据,自定义
大气模型。
标准的大气模型有:
1)1976年美国“标准”大气:
1976年,由美国制定,用中纬度平均值表示。
2)赤道:
北纬15°。
3)中纬度夏天:
北纬45°,7月。
4)中纬度冬天:
北纬45°,1月;
5)亚北极夏天:
北纬60°,7月;
6)亚北极冬天:
北纬60°,1月
MODTRAN的模型大气将大气非等高度地划分为34层,能给出这6种模型大
气的气压、温度以及H2O、CO2、O3等11种气体含量随高度的分布。
例如:
中纬
度夏天模型大气的气压、温度及水分子含量的高度分布如下表所示。
7
3MODTRAN屏幕输入
MODTRAN的屏幕输入参数有模型大气(ModelAtmosphere)、气溶胶
(Aerosol)、几何参数和光谱带(GeometryandSpectralBand)等3大类,参数输入
屏幕上分别用
(1)、
(2)、(3)标注。
3.1“模型大气”输入参数
屏幕标题:
ModelAtmosphere
(1)
8
3.1.1计算选择
可选MODTRAN带模型或MODTRAN所含的LOWTRAN模型。
3.1.2大气模型
大气模型包括标准大气模型和用户自定义。
标准模型大气与LOWTRAN相
同,共6种,其中1976年美国标准大气是一种标准的大气模型,其余5种模型反映
了季节、纬度对大气性质的影响。
用户自定义大气模型有两种,要求用户输入气
象数据。
3.1.3路径类型
大气路径有水平路径、斜路径、斜路径至太空3种。
需输入的路径参数如前
所述,在“几何参数和光谱带”屏幕上输入。
3.1.4运行方式
MODTRAN有透过率(Transmittance)、热辐射亮度(ThermalRadiance)、
有散射辐亮度(RadiancewithScattering)、太阳直射照度(DirectSolarIrradiance)
等四种运行模式
3.1.4.1透过率
透过率运行模式(MODE=0)可计算路径的总透过率以及气体分子带吸收、连
续吸收、气溶胶吸收等分量的路径透过率。
路径的总透过率为各个分量的路径透
过率之积。
透过率运行模式下,可输出以下结果:
1)总路径透过率;
2)H2O,CO2,O3等吸收气体的透过率;
3)分子散射的透过率;
9
4)气溶胶透过率。
3.1.4.2热辐射亮度
热辐射亮度模式(MODE=1)可计算路径大气辐射的辐亮度(即大气辐亮度)
和路径的总透过率,大气辐射主要在热红外波段,故称大气热辐射。
热辐射亮度运行模式下可输出以下结果:
1)总路径透过率(Trans.Total)
-与路径类型、路径大气吸收、散射的衰减系数有关。
2)路径热辐射(PathThermal)
-路径大气热辐射的光谱辐亮度。
3)热散射(ThermalScat)
-路径外大气热辐射经散射进入路径的光谱辐亮度,通常可忽略。
4)表面辐射(SurfaceEmission)
-特定温度、比辐射率边界层表面热辐射产生的光谱辐亮度。
5)地面反射的总辐亮度(TotalGroundReflected)
-大气热辐射经地面反射产生的光谱辐亮度,通常可忽略。
10
6)总辐亮度(TotalRadiance)
-观察点在视线方向接收到辐射的总辐亮度。
总辐亮度是2)、3)、4)、5)
四项辐亮度之和,由于热散射3)和地面反射辐亮度5)可忽略,总辐亮度
是路径热辐射2)和边界层热辐射4)辐亮度之和。
7)光学深度(OpticalDepth)
-单位:
Km-1,应是衰减系数,光学深度(大气质量)是衰减系数对路径
的积分。
边界层可以是地球、云、飞机等,MODTRAN假设边界层均为灰体,用户
须输入边界层温度和表面反射率(SurfaceAlbedo),以确定表面发射的辐亮度。
边界层温度的缺省值为0,此时边界层温度根据大气模型确定。
如不用缺省值,
用户可输入边界层的绝对温度和表面反射率。
比辐射率等于(1-表面反射率),
表面反射率的缺省值为0,比辐射率等于1,即认为是黑体。
用户可直接输入反射率的数值,MODTRAN也提供了几种典型的边界层,包
括:
雪、森林、植被、草地、海洋、沙漠等,这几种典型边界层在红外波段的光
谱反射率和光谱比辐射率都是确定的。
例如:
雪
Snowcover(fresh)
Spectralcoverageis1.4to50.micrometers(200to7143wavenumbers)
WavelengthRHOEPS
反射率比辐射率
1.40.48.52
1.46.05.95
1.55.05.95
1.85.24.76
1.93.02.98
2.08.02.98
2.25.17.83
2.35.06.94
2.45.02.98
3.00.02.98
11
3.20.07.93
5..02.98
20..02.98
50..02.98
森林
沙漠
云
12
海洋
Ocean(notgrazingangles)
Spectralcoverageis1.0to14.micrometers(714to10000wavenumbers)
WavelengthRHOEPS
1.0.03.97
1.5.03.97
2.0.02.98
2.6.01.99
2.9.05.95
3.4.02.98
4.0.02.98
6.0.02.98
8.0.01.99
10..02.98
12..04.96
14..06.94
可用热辐射亮度模式计算HY1海洋卫星以45°倾角观察海面的路径辐亮度
和透过率。
输入文件:
PcModWin37\Usrer\HY1热辐射辐亮度(全波段).ltn
输入参数:
中纬度、夏天,23公里能见度,海水温度300K,反射率0.03。
13
输出曲线如图,大气热散射、地面反射的热辐射均可忽略,观察点的总辐亮
度是路径热辐射和边界层热辐射的辐亮度之和。
卫星红外地平仪接收到的15微米
辐射几乎全部来自CO2吸收带的大气辐射,来自地球表面的辐射可以忽略。
3.1.4.3有散射辐亮度
热辐射亮度运行模式的路径辐亮度只考虑路径大气和地面的热辐射,有散射
辐亮度运行模式(MODE=2)的路径辐亮度不仅考虑路径大气和地面的热辐射,
还增加了路径外辐射源(太阳、月亮、地球等)产生的大气散射辐亮度和地面反
射辐亮度。
路径大气的散射有分子散射、气溶胶散射,分子散射主要影响可见光波段的
14
路径辐亮度,气溶胶散射对近红外、中波红外影响较大。
路径外辐射源的辐射经
大气散射进入路径将增加路径的辐亮度。
MODTRAN的单次散射仅计算太阳或
月亮光散射后进入路径,不考虑进入路径散射光再次散射,离开路径。
散射辐亮度和观察方向与太阳/月亮照射的相对位置遥感,计算时须输入太阳
/月亮几何位置(Solar/LunarGeometry)相关参数。
有散射辐亮度运行模式下,可输出以下结果:
1)总路径透过率(Trans.Total)
-与路径类型、路径大气吸收、散射的衰减系数有关。
2)路径热辐射(PathThermal)
-路径大气热辐射的辐亮度。
3)热散射(ThermalScat)
-大气热辐射经散射进入路径的辐亮度,通常可忽略。
4)表面辐射(SurfaceEmission)
-特定温度、比辐射率界面热辐射产生的辐亮度。
5)太阳散射辐亮度(SolarScatterRadiance)。
-太阳辐射经散射进入路径并到达观察点的辐亮度。
6)单次散射辐亮度(SingleScatterRadiance)
-太阳辐射经单次散射进入路径的辐亮度。
单次散射是指不考虑进入路径的
散射辐射再次离开路径,如设置为单次散射模式,太阳散射辐亮度即单次散
射辐亮度。
7)总地面反射(TotalGroundReflected)
-太阳光地面反射产生的辐亮度,包括太阳直射光地面反射和散射光地面
反射,能量主要集中在VIS,NIR波段,少量在,SW波段,MW,LW波段
可忽略。
8)直接地面反射(DirectGroundReflected)
-太阳直射的地面反射的辐亮度(用户手册p380)。
9)总辐亮度(TotalRadiance)
-观察点在视线方向接收到的辐亮度,是路径大气热辐射、地面热辐射、
15
地面反射产生的辐亮度之和。
在MW,LW波段,地面反射辐亮度可忽略,
总辐亮是路径大气热辐射、地面热辐射辐亮度之和。
在VIS、NIR、SW波
段,路径大气热辐射可忽略,总辐亮度是路径大气太阳散射和地面反射辐
亮度之和。
10)太阳反射辐亮度(ReflectedSolar)
-大气层外太阳光100%漫反射产生的辐亮度,可作为外层空间太阳辐射
基准。
11)观察点太阳反射辐亮度(SolaratObserver)
-观察点处太阳光100%漫反射产生的辐亮度,可作为观察点处太阳辐射
的基准。
12)光学深度(OpticalDepth)
-单位:
Km-1,应是衰减系数,光学深度(大气质量)是衰减系数对路径
的积分。
可用有散射辐射亮度模式计算HY1海洋卫星以45°倾角观察海面的路径辐亮度
和透过率。
输入文件:
HY1太阳散射辐亮度(全波段).ltn。
输入参数:
太阳天
顶角为0°。
中纬度、夏天,23公里能见度,海水温度300K,反射率0.03。
输出的全波段路径辐亮度曲线如图。
中波、长波波段的路径辐亮度是海面和路
径大气的热辐射辐亮度之和,大气散射、海面反射均可忽略。
16
在可见光至短波红外波段,路径辐亮度是太阳散射和海面辐射辐亮度之和,太
阳散射远强于海面反射。
3.1.4.4太阳直射照度
太阳直射照度(“DirectSolarIrradiance”)运行模式(MODE=3)可计算:
1)大气层外的太阳正入射的光谱辐照度;
2)太阳透过大气到达观察点的光谱辐照度;
3)路径透过率。
这里的大气层外太阳照度或观察点太阳照度约定为正入射照度。
太阳的入射照度和地面照度(天顶角0°,VIS-LWIR)
17
太阳的入射照度和地面照度(天顶角0°,VIS-NIR)
至海面大气透过率(天顶角0°,VIS-NIR)
3.2“气溶胶”输入参数
屏幕标题:
Aerosole
(2)
3.2.1气溶胶模型
NoAerosolAttenuation=0=noaerosoleffectsincludedinthecalculation
Rural-VIS=23km=1=RURALextinction,defaultVIS=23km
Rural-VIS=5km=2=RURALextinction,defaultVIS=5km
NavyMaritime=3=NAVYMARITIMEextinction,setsitsownVIS
Maritime-VIS=23km=4=MARITIMEextinction,defaultVIS=23km
Urban-VIS=5km=5=URBANextinction,defaultVIS=5km
18
Tropospheric-VIS=50km=6=TROPOSPHERICextinction,defaultVIS=50km
UserDefined-VIS=23km=7=USERDEFINEDextinction,defaultVIS=23km
Fogadvection-VIS=.5km=8=FOG1(advectionfog)extinction,defaultVIS=0.2km
Fogradiation-VIS=.2km=9=FOG2(radiationfog)extinction,defaultVIS=0.5km
Desertextinction=10=DESERTextinction,setsitsownVISfromWSS
3.2.2气溶胶模型季节修正
主要对季节变化引起对流层高度的变化进行修正。
3.3“几何位置和光谱带”输入参数
屏幕提示:
GeometryandSpectralBand(3)
3.3.1路径参数
须先在屏幕
(1)的“TypeofAtmospherePath”的框中选择大气路径类型(如
水平路径、斜路径等),再在屏幕(3)“PathType”的框中选择路径参数设置方法
(如高度、天顶角、路径长度等),方可输入具体路径参数。
每个路径参数均可
设置初值、终值、步长,即可得到多组路径参数的计算结果。
3.3.2光谱带参数
须输入的参数包括:
初始、终值、增量或狭缝函数的半宽度(FWHM)的波数
(cm-1)或波长(um或nm)。
3.3.3太阳/月亮几何位置
选择“有散射辐亮度”运行方式时,先要明确散射源是太阳或月亮。
然后输
入太阳/月亮的几何位置参数,输入方法有3种:
1)输入观察点和太阳的地平经度、地平纬度;
输入内容包括:
太阳日:
1-365,仅用于修正日地距离;
观察点经纬度:
地经(0°-360°)、地纬(-90°-90°)
19
太阳经纬度:
地心与太阳连线与地球交点的地经、地纬;
路径方位角:
0-360°(正北为0°,正东为90°)。
2)输入观察点经纬度、太阳日(1-365)和观察时刻的格林威治时间;
输入内容包括:
太阳日:
1-365,用于修正日地距离和太阳的视位置;
观察点经纬度:
地经(0°-360°)、地纬(-90°-90°)
时间:
观察时刻的Greenwich时间;
路径方位角:
0-360°(正北为0°,正东为90°)。
3)输入观察视线方位角(视线与日地连线地面投影的夹角)和观察点的太
阳天顶角
输入内容包括:
太阳日:
1-365,仅用于修正日地距离;
视线方位角:
-180°-180°,所谓视线方位角实际上是视线转至日地连
线在方位方向转过的角度,由北往东为正,反之为负。
太阳天顶角:
0°-90°;
sun
path
N
W
los
path
E
S
太阳光散射与散射角(太阳光线与视线之间的夹角)有关。
方法1用太
阳的地经、地纬确定太阳光的入射方向,视线方向由路径的天顶角、方位角
20
而定。
方法2的太阳光入射方向是根据太阳日、时间而定,视线方向仍由路
径的天顶角、方位角而定。
如何确定观察点太阳的天顶角、方位角可参见附
录2。
由于路径天顶角已设定,方法3要求输入太阳天顶角和观察视线与日地
连线之间的方位角,已能完全确定太阳光与视线的相对位置,不必再要求输
入太阳和视线的方位角。
4附录1:
MO