附件 1:6S 模式使用说明
一、输入参数和说明
1.1 几何参数(geometrical parameters)
参数限制:无
参数名称:igeom
取值范围:0-7
igeom=0:用户自己选择观测几何参数
所需参数有:
太阳天顶角(度) 太阳方位角(度) 卫星天顶角(度)
卫星方位角(度) 月(1-12)日(1-31)
igeom=1-7 分别代表以下卫星的观测:
igeom=1:Meteosat 卫星
所需参数有:
月 日 世界时(十进制)
列数 行数(图象最大尺度为 5000X2500 像素)
igeom=2:GOES(东)卫星
所需参数有:
月 日 世界时(十进制)
列数 行数(图象最大尺度为 17000X12000 像素)
igeom=3:GOES(西)卫星
所需参数有:
月 日 世界时(十进制)
列数 行数(图象最大尺度为 17000X12000 像素)
igeom=4:AVHRR,下午 NOAA 卫星
所需参数有:
月 日 世界时(十进制)
列数(1-2048)经度 穿越赤道时间
igeom=5:AVHRR,上午 NOAA 卫星
所需参数有:
月 日 世界时(十进制)
列数(1-2048)经度 穿越赤道时间
26
igeom=6:HRV(SPOT)
所需参数有:
月 日 世界时(十进制) 经度 纬度
igeom=7:TM(LANDSAT)
所需参数有:
月 日 世界时(十进制) 经度 纬度
注:对 HRV 和 TM,经纬度代表图象中心位置
1.2 大气模式(atmospheric model)
参数限制:无
参数名称:idatm
取值范围:0-9
idatm =0:无气体吸收
idatm =1:热带大气
idatm =2:中纬度夏大气
idatm =3:中纬度冬季
idatm =4:亚北极区夏季
idatm =5:亚北极区冬季
idatm =6:美国标准大气(62 年)
idatm =7:用户定义大气廓线(34 层无线电探空数据)
包括:高度(km ) 气压( mb ) 温度( k ) 水汽密度( g/m3) 臭氧密度(g/m3)
idatm =8:输入水汽和臭氧总含量
水汽( g/cm2 ) 臭氧 (cm-atm)
idatm =9:读入无线电探空数据文件
1.3 气溶胶类型参数
参数限制:无
参数名称:iaer
取值范围:0-12
iaer=0: 无气溶胶
iaer=1: 大陆型气溶胶
iaer=2: 海洋型气溶胶
iaer=3: 城市气溶胶
27
iaer=5: 沙漠型气溶胶
iaer=6: 生物质燃烧型
iaer=7: 平流层模式
iaer=4: 用户自己输入以下四种粒子所占体积百分比(0-1)
c(1) : 灰尘
c(2) :水溶型
c(3) :海洋型
c(4) :烟灰
iaer=8-10:用户自己按照尺度分布类型定义气溶胶模型
iaer=8:多峰对数正态分布
iaer=9:改进的 gamma 分布
iaer=10:Junge 幂指数律分布
iaer=11:按太阳光度计测量结果定义气溶胶模型
需 要 输 入 参 数 有 : 粒 子 半 径 ( µm ) 粒 径 分 布 ( d V / d (logr),
cm3/cm2/micron)和复折射指数的实部和虚部谱
iaer=12:利用事先计算的结果
给出文件名
1.4 气溶胶含量参数(concentration)
参数限制:能见度必须大于 5 公里
参数名称:v
取值范围:
v=能见度(公里)
v=0:输入 550 纳米气溶胶光学厚度
taer55=550 纳米气溶胶光学厚度
v=-1:没有气溶胶
1.5 目标高度参数(altitude of target)
参数限制:无
参数名称:xps
取值范围:
xps >=0:目标在海平面高度
xps < 0:绝对值代表目标高度(公里)
28
1.6 传感器高度参数(sensor altitude)
参数限制:无
参数名称:xpp
取值范围:
xpp= -1000:卫星观测
xpp=
0:地面观测
-100< xpp <0:飞机观测,绝对值代表飞机相对于目标的高度(公里)
对于飞机观测,必须输入飞机和地面之间的水汽,臭氧含量和 550 纳米
气溶胶光学厚度,如无数据则输入负值,水汽和臭氧根据 62 年美国标准
大气内差,气溶胶则根据 2 公里指数廓线计算
1.7 光谱参数(spectral conditions)
参数限制:虽然在整个波段计算气体透射率和散射函数,但处理强吸收
波段吸收与散射的相互作用不精确,因此不适合强吸收带
参数名称:iwave
取值范围:-2 – 70
iwave=-2 – +1, 用户自己定义光谱条件
iwave=-2:用户输入光谱范围的下限和上限(微米),滤光片函数为 1,
输出文件中给出单色结果。
iwave=-1:单色计算,用户给出单色波长(微米)
iwave=0:用户输入光谱范围的下限和上限(微米),滤光片函数为 1
iwave=1:用户输入光谱范围的下限和上限(微米)并输入滤光片函数(间
隔为 0.0025 微米。
iwave=2-70:选择下列卫星通道
2
3
4
5
6
7
8
( 0.350-1.110 )
vis band of meteosat
( 0.490-0.900 )
vis band of goes east
vis band of goes west
( 0.490-0.900 )
1st band of avhrr(noaa6) ( 0.550-0.750 )
2nd
( 0.690-1.120 )
1st band of avhrr(noaa7) ( 0.500-0.800 )
2nd
( 0.640-1.170 )
"
"
29
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
1st band of avhrr(noaa8) ( 0.540-1.010 )
"
"
"
"
"
"
"
"
"
"
2nd
( 0.680-1.120 )
1st band of avhrr(noaa9) ( 0.530-0.810 )
2nd
( 0.680-1.170 )
1st band of avhrr(noaa10 ( 0.530-0.780 )
2nd
( 0.600-1.190 )
1st band of avhrr(noaa11 ( 0.540-0.820 )
2nd
( 0.600-1.120 )
1st band of hrv1(spot1) ( 0.470-0.650 )
2nd
( 0.600-0.720 )
( 0.730-0.930 )
3rd
pan
( 0.470-0.790 )
1st band of hrv2(spot1) ( 0.470-0.650 )
( 0.590-0.730 )
2nd
( 0.740-0.940 )
3rd
pan
( 0.470-0.790 )
1st band of tm(landsat5) ( 0.430-0.560 )
( 0.500-0.650 )
2nd
3rd
( 0.580-0.740 )
( 0.730-0.950 )
4th
( 1.5025-1.890 )
5th
7th
( 1.950-2.410 )
1st band of mss(landsat5)( 0.475-0.640 )
( 0.580-0.750 )
2nd
( 0.655-0.855 )
3rd
4th
( 0.785-1.100 )
( 0.5025-0.5875)
1st band of MAS (ER2)
( 0.6075-0.7000)
2nd
3rd
( 0.8300-0.9125)
( 0.9000-0.9975)
4th
( 1.8200-1.9575)
5th
6th
( 2.0950-2.1925)
( 3.5800-3.8700)
7th
( 0.6100-0.6850)
MODIS
MODIS
( 0.8200-0.9025)
( 0.4500-0.4825)
MODIS
MODIS
( 0.5400-0.5700)
"
"
"
"
"
"
band 1
band 2
band 3
band 4
"
"
"
"
"
"
"
"
30
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
66
67
68
69
70
"
P1
P1
band 5
band 6
band 7
( 1.2150-1.2700)
( 1.6000-1.6650)
( 2.0575-2.1825)
MODIS
MODIS
MODIS
1st band of avhrr(noaa12 ( 0.500-1.000 )
2nd
( 0.650-1.120 )
1st band of avhrr(noaa14 ( 0.500-1.110 )
( 0.680-1.100 )
2nd
"
band 1
POLDER
( 0.4125-0.4775)
band 2 (non polar( 0.4100-0.5225)
POLDER
POLDER
band 3 (non polar( 0.5325-0.5950)
band 4
POLDER
( 0.6300-0.7025)
band 5 (non polar( 0.7450-0.7800)
POLDER
POLDER
band 6 (non polar( 0.7000-0.8300)
band 7
POLDER
( 0.8100-0.9200)
band 8 (non polar( 0.8650-0.9400)
POLDER
FY-1C
band 1
( 0.5310-0.7490)
( 0.7610-0.9990)
band 2
FY-1C
( 1.4950-1.7330)
band 6
FY-1C
FY-1C
band 7
( 0.4000-0.5900)
( 0.4010-0.6190)
band 8
FY-1C
( 0.4330-0.6710)
band 9
FY-1C
FY-1C
band 10
( 0.8320-1.0700)
1.8 地表反射率类型(ground reflectance type)
参数限制:用户可以选择“补丁”结构的地表情况,即输入一个半径为
rad 的圆形目标的反射率 roc 和周围环境的反射率 roe
参数名称:inhomo
取值范围:0,1
inhomo=0:均匀表面
所需参数:
idirec=0: 无方向效应
输入均匀朗伯表面的反射率 igroun(roc=roe)
idirec=1: 有方向效应
ibrdf=0: 输入太阳天顶角为 thetas 时 10 个观测天顶角(0-80
度间隔 10 度和 85 度)和 30 个观测方位角(0-360 度
31
间隔 30 度)下的反射率;
同样,输入观测天顶角为 thetav 时各太阳入射角度下
的反射率;
地表半球反射率;
在所选的观测条件下(太阳天顶角,观测天顶角和相
对方位角)的反射率;
ibrdf=1-9: 选择模式中储存的模式
ibrdf=1: hapke model
ibrdf=2: verstraete et al. model
ibrdf=3: Roujean et al. model
ibrdf=4: walthall et al. model
ibrdf=5: minnaert model
ibrdf=6: Ocean
ibrdf=7: Iaquinta and Pinty model
ibrdf=8: Rahman et al. model
ibrdf=9: Kuusk's multispectral CR model
对于上述每种地面反射率模式,还分别需要输入各自所需的参数,请参
阅主程序说明。
1.9 激活大气订正方式
参数限制:无
参数名称:rapp
取值范围:
rapp<-1: 不激活大气订正方式
rapp> 0:反演地面反射率,反演出的地面反射率满足大气层顶
的辐射亮度=rapp(w/m2/str/mic)
-1.
9
07071010.txt
5
0
43.8231 98.8277 23.53 99.03 7 7 输入几何条件参数,分别为太阳
天顶角、方位角,卫星天顶角、
方位角,月、日
选择输入探空数据,臭氧
取自中纬度夏季模式
探空数据文件
气溶胶模型为沙漠型
选择气溶胶含量输入方式为 550
纳米光学厚度
550 纳米气溶胶光学厚度
目标高度为 1.139 KM
传感器在卫星上
光谱条件为 FY-1C 第 1 通道
选择均匀表面
选择无方向性反射特征
反射率光谱参数从外部文件读入
反射率数据文件
不激活大气订正方式
0.087
-1.139
-1000
61
0
0
5
MR07079.TXT
-2
探空数据格式(举例):
Radiosound(070710.10)
RelHum
Wind
-------------------------------------------------
P(pa)
1000
925
882
850
700
609
600
500
400
300
250
200
150
100
H(m)
0.000
0.000
1.140
1.462
3.121
4.254
4.377
5.812
7.505
9.568
10.810
12.299
14.144
16.641
0
0
0
43002
52004
0
229002
272010
288012
280018
282038
270039
270037
273018
T(c)
0.0
0.0
22.4
22.0
10.9
0.0
-0.9
-8.6
-20.1
-36.8
-42.7
-48.9
-58.5
-67.4
100
100
42
26
29
38
36
29
25
23
22
20
19
999
Td(c)
0.0
0.0
8.9
1.6
-6.4
-12.6
-14.1
-23.4
-35.0
-50.2
-55.8
-61.9
-70.6
99.9
33