辐射模型使用说明

PAGEPAGE346S模式使用说明一、输入参数和说明1.1几何参数（geometricalparameters）参数限制：无参数名称：igeom取值范围：0-7igeom=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）经度穿越赤道时间igeom=6：HRV(SPOT)所需参数有：月日世界时（十进制）经度纬度igeom=7：TM(LANDSAT)所需参数有：月日世界时（十进制）经度纬度注：对HRV和TM，经纬度代表图象中心位置1.2大气模式（atmosphericmodel）参数限制：无参数名称：idatm取值范围：0-9idatm=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-12iaer=0:无气溶胶iaer=1:大陆型气溶胶iaer=2:海洋型气溶胶iaer=3:城市气溶胶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）粒径分布（dV/d(logr)，cm3/cm2/micron）和复折射指数的实部和虚部谱iaer=12：利用事先计算的结果给出文件名1.4气溶胶含量参数（concentration）参数限制：能见度必须大于5公里参数名称：v取值范围：v=能见度（公里）v=0：输入550纳米气溶胶光学厚度taer55=550纳米气溶胶光学厚度v=-1：没有气溶胶1.5目标高度参数（altitudeoftarget）参数限制：无参数名称：xps取值范围：xps>=0：目标在海平面高度xps<0：绝对值代表目标高度（公里）1.6传感器高度参数（sensoraltitude）参数限制：无参数名称：xpp取值范围：xpp=-1000：卫星观测xpp=0：地面观测-100<xpp<0：飞机观测，绝对值代表飞机相对于目标的高度（公里）对于飞机观测，必须输入飞机和地面之间的水汽，臭氧含量和550纳米气溶胶光学厚度，如无数据则输入负值，水汽和臭氧根据62年美国标准大气内差，气溶胶则根据2公里指数廓线计算1.7光谱参数（spectralconditions）参数限制：虽然在整个波段计算气体透射率和散射函数，但处理强吸收波段吸收与散射的相互作用不精确，因此不适合强吸收带参数名称：iwave取值范围：-2–70iwave=-2–+1，用户自己定义光谱条件iwave=-2：用户输入光谱范围的下限和上限（微米），滤光片函数为1，输出文件中给出单色结果。iwave=-1：单色计算，用户给出单色波长（微米）iwave=0：用户输入光谱范围的下限和上限（微米），滤光片函数为1iwave=1：用户输入光谱范围的下限和上限（微米）并输入滤光片函数（间隔为0.0025微米。iwave=2-70：选择下列卫星通道2visbandofmeteosat(0.350-1.110)3visbandofgoeseast(0.490-0.900)4visbandofgoeswest(0.490-0.900)51stbandofavhrr(noaa6)(0.550-0.750)62nd"(0.690-1.120)71stbandofavhrr(noaa7)(0.500-0.800)82nd"(0.640-1.170)91stbandofavhrr(noaa8)(0.540-1.010)102nd"(0.680-1.120)111stbandofavhrr(noaa9)(0.530-0.810)122nd"(0.680-1.170)131stbandofavhrr(noaa10(0.530-0.780)142nd"(0.600-1.190)151stbandofavhrr(noaa11(0.540-0.820)162nd"(0.600-1.120)171stbandofhrv1(spot1)(0.470-0.650)182nd"(0.600-0.720)193rd"(0.730-0.930)20pan"(0.470-0.790)211stbandofhrv2(spot1)(0.470-0.650)222nd"(0.590-0.730)233rd"(0.740-0.940)24pan"(0.470-0.790)251stbandoftm(landsat5)(0.430-0.560)262nd"(0.500-0.650)273rd"(0.580-0.740)284th"(0.730-0.950)295th"(1.5025-1.890)307th"(1.950-2.410)311stbandofmss(landsat5)(0.475-0.640)322nd"(0.580-0.750)333rd"(0.655-0.855)344th"(0.785-1.100)351stbandofMAS(ER2)(0.5025-0.5875)362nd"(0.6075-0.7000)373rd"(0.8300-0.9125)384th"(0.9000-0.9975)395th"(1.8200-1.9575)406th"(2.0950-2.1925)417th"(3.5800-3.8700)42MODISband1(0.6100-0.6850)43MODISband2(0.8200-0.9025)44MODISband3(0.4500-0.4825)45MODISband4(0.5400-0.5700)46MODISband5(1.2150-1.2700)47MODISband6(1.6000-1.6650)48MODISband7(2.0575-2.1825)491stbandofavhrr(noaa12(0.500-1.000)502nd"(0.650-1.120)511stbandofavhrr(noaa14(0.500-1.110)522nd"(0.680-1.100)53POLDERband1(0.4125-0.4775)54POLDERband2(nonpolar(0.4100-0.5225)55POLDERband3(nonpolar(0.5325-0.5950)56POLDERband4P1(0.6300-0.7025)57POLDERband5(nonpolar(0.7450-0.7800)58POLDERband6(nonpolar(0.7000-0.8300)59POLDERband7P1(0.8100-0.9200)60POLDERband8(nonpolar(0.8650-0.9400)61FY-1Cband1(0.5310-0.7490)62FY-1Cband2 (0.7610-0.9990)66FY-1Cband6 (1.4950-1.7330)67FY-1Cband7 (0.4000-0.5900)68FY-1Cband8 (0.4010-0.6190)69FY-1Cband9 (0.4330-0.6710)70FY-1Cband10 (0.8320-1.0700)1.8地表反射率类型（groundreflectancetype）参数限制：用户可以选择“补丁”结构的地表情况，即输入一个半径为rad的圆形目标的反射率roc和周围环境的反射率roe参数名称：inhomo取值范围：0，1inhomo=0：均匀表面所需参数：idirec=0：无方向效应输入均匀朗伯表面的反射率igroun（roc=roe）idirec=1：有方向效应ibrdf=0：输入太阳天顶角为thetas时10个观测天顶角（0-80度间隔10度和85度）和30个观测方位角（0-360度间隔30度）下的反射率；同样，输入观测天顶角为thetav时各太阳入射角度下的反射率；地表半球反射率；在所选的观测条件下（太阳天顶角，观测天顶角和相对方位角）的反射率；ibrdf=1-9：选择模式中储存的模式ibrdf=1：hapkemodelibrdf=2：verstraeteetal.modelibrdf=3：Roujeanetal.modelibrdf=4：walthalletal.modelibrdf=5：minnaertmodelibrdf=6：Oceanibrdf=7：IaquintaandPintymodelibrdf=8：Rahmanetal.modelibrdf=9：Kuusk'smultispectralCRmodel对于上述每种地面反射率模式，还分别需要输入各自所需的参数，请参阅主程序说明。1.9激活大气订正方式参数限制：无参数名称：rapp取值范围：rapp<-1:不激活大气订正方式rapp>0：反演地面反射率，反演出的地面反射率满足大气层顶的辐射亮度=rapp（w/m2/str/mic）-1.<rapp<0：反演地面反射率，反演出的地面反射率满足表观反射率=rapp输入文件举例下面是敦煌辐射校正场同步观测大气辐射传输模拟计算的输入文件：输入参数说明0选择用户自己选择几何条件43.823198.827723.5399.0377输入几何条件参数，分别为太阳天顶角、方位角，卫星天顶角、方位角，月、日9选择输入探空数据，臭氧取自中纬度夏季模式07071010.txt探空数据文件5气溶胶模型为沙漠型0选择气溶胶含量输入方式为550纳米光学厚度0.087550纳米气溶胶光学厚度-1.139目标高度为1.139KM-1000传感器在卫星上61光谱条件为FY-1C第1通道0选择均匀表面0选择无方向性反射特征5反射率光谱参数从外部文件读入MR07079.TXT反射率数据文件-2不激活大气订正方式探空数据格式（举例）：Radiosound(070710.10)P(pa)H(m)T(c)Td(c)RelHumWind10000.0000.00.010009250.0000.00.010008821.14022.48.94208501.46222.01.626430027003.12110.9-6.429520046094.2540.0-12.63806004.377-0.9-14.1362290025005.812-8.6-23.4292720104007.505-20.1-35.0252880123009.568-36.8-50.22328001825010.810-42.7-55.82228203820012.299-48.9-61.92027003915014.144-58.5-70.61927003710016.641-67.499.99992730187018.794-62.3999.09992280075020.891-57.7999.0999940064022.302-57.2999.0999880093024.142-52.2999.099987012反射率数据文件格式（举例）：350.000.1044351.000.1046352.000.1048353.000.1055354.000.1056355.000.1057356.000.1062***.***.****第一列为波长(纳米),第二列为反射率

输出文件举例下面是对应于上述输入文件的输出结果（部分）：integratedvaluesof:apparentreflectance.1369appar.rad.(w/m2/sr/mic)26.279totalgaseoustransmittance.546说明：上述结果为表观反射率，表观辐射亮度和气体总透射率integratedvaluesof:downwardupwardtotalglobalgas.trans.:.63343.66785.54570water"":.63343.66785.54570ozone"":1.000001.000001.00000co2"":1.000001.000001.00000oxyg"":1.000001.000001.00000no2"":1.000001.000001.00000ch4"":1.000001.000001.00000co"":1.000001.00000