您好,欢迎访问三七文档
参考源码中的“main.f”介绍6S可以用来预测0.25~4.0μm间无云大气条件下传感器接收到的辐射。考虑了watervapor,carbondioxyde,oxygenandozone的吸收作用,moleculesandaerosols的散射作用。不均匀地表、二向反射影响也得到了考虑。总体需要以下参数输入:几何关系大气模型用于描述大气成分气溶胶类型光谱条件地表反射(类型与光谱差异)用户可以选择标准类型,也可以自定义数据。参数有效范围geometricalparametersatmosphericmodelaerosolmodelthevisibilitymustbebetterthan5.0kmforsmallervaluesalculationsmightbenomorevalid.Spectralconditionsthegaseoustransmittanceandthescatteringfuntionsarevalidfrom0.25to4.0micron.butthetreatmentofinteractionbetweenabsorptionandscatteringiscorrectfornottoolargeabsorption.ifyouwanttoocomputesignalwithinabsorptionbands,thisinteractionoughttobereconsideredgroundreflectance(type)youcanconsiderapatchystructure:thatisacirculartargetofradiusradandofreflectanceroc,withinanenvironmentofreflectanceroe.Foruniformsurfaceconditionsonly,youmayconsiderdirectionalreflectanceasboundaryconditions.Someanalyticalmodelareproposed,theusercanspecifyhisownvalues.Thecodeassumesthatthebrdfisspectrallyindependentgroundreflectance(spectralvariation)Fourtypicalreflectancesareproposed,definedwithingivenspectralrange.thisrangediffersaccordingtotheselectedcase.thereflectanceissetto0outsidethisrange,duetothedeficiencyofdata.usermustverifytheselimits.thatisobviouslyirrelevantforbrdf1.输入几何条件igeom取值0,1-70:用户输入太阳天顶角太阳方位角卫星天顶角卫星方位角月天以上参数按顺序用空格隔开1:meteosatobservation用户需要再输入月天十进制的小时(universaltime)列行(fullscal5000*2500)2:goeseastobservation用户需要再输入月天十进制的小时(universaltime)列行(fullscal17000*12000)3:goeswestobservation用户需要再输入月天十进制的小时(universaltime)列行(fullscal17000*12000)4:avhrr(PMnoaa)用户需要再输入月天十进制的小时(universaltime)列(1-2048)xlonanhna5:avhrr(AMnoaa)用户需要再输入月天十进制的小时(universaltime)列(1-2048)xlonanhnalong.(xlonan)andoverpasshour(hna)attheascendantnodeatequatorxlonan和hna分别是卫星升交点与赤道的角度经度与时间6:hrv(spot)用户需要再输入月天十进制的小时经度纬度7:tm(landsat)用户需要再输入月天十进制的小时经度纬度2.输入大气模型idatm取值0-80:nogaseousabsorption没有大气吸收1:tropical2:midlatitudesummer3:midlatitudewinter4:subarcticsummer5:subarcticwinter6:usstandard627:用户输入大气廓线radiosonde(无线探空仪)dataon34levels,每一个level输入以下5个数据:enteraltitude(inkm)pressure(inmb)temperature(ink)h2odensity(ing/m3)o3density(ing/m3)8:输入水汽含量和臭氧含量,大气廓线使用us62uw单位g/cm2uo3单位cm-atm3.输入气溶胶类型iaer取值-1,0~3,5~7,4,8~9,11,12-1:user-definedprofile先输入层数numberoflayers,在依次输入各层参数:高度(km),550nm光学厚度,气溶胶类型(下面的1-3,5-7)Exampleforiaer=-1:42.00.200110.00.02518.00.003180.00.00010:没有气溶胶1:大陆型2:海洋型3:城市5:沙漠backgrounddesert6:biomassburning7:平流层stratospheric4:defineyourownmodelusingbasiccomponentsc(1)=volumetric%ofdust-likec(2)=volumetric%ofwater-solublec(3)=volumetric%ofoceanicc(4)=volumetric%ofsoot煤灰,烟灰between0to1defineyourownmodelusingasizedistributionfunction8:MultimodalLog-Normaldistribution(upto4modes)9:ModifiedGammadistribution10:JungePower-Lawdistributiondefineamodelusingsun-photometermeasurements11:SunPhotometerdistribution(50valuesmax)Youhavetoenter:randdV/d(logr)whereristheradius(inmicron),Visthevolume,anddV/d(logr)isin(cm3/cm2/micron).Thenyouhavetoenter:nrandniforeachwavelengthwherenrandniarerespectivelytherealandtheimaginarypartsoftherefractiveindexoryoucanusetheresultscomputedandsavedpreviously12:ReadingofdatapreviouslysavedintoFILEyouhavetoentertheidentificationnameFILEinthenextlineofinputs.Note:Foriaer=8,9,10,and11:ResultsfromtheMIEsubroutinemaybesavedintothefileFILE.mie(Extinctionandscatteringcoefficients,singlescatteringalbedo,asymmetryparameter,phasefunctionatpredefinedwavelengths)andthencanbere-usedwiththeoptioniaer=12,whereFILEisanidentificationnameyouhavetoenter.So,ifyouselectiaer=8,9,10,or11,youwillhavetoenteriaerpaftertheinputsrequestedbyoptions8,9,10,or11:iaerp=0resultswillnotbesavediaerp=1resultswillbesavedintothefileFILE.mie.NextlineenterFILE.Exampleforiaerandiaerp8MultimodalLog-Normaldistributionselected0.001203Rmin,Rmax,3components0.4712.5120.17Rmean,Sigma,%density-1stcomponent1.531.531.531.531.531.531.531.531.531.531.531.531.5281.521.4621.41.3681.2761.221.2nrfor20wavelengths0.0080.0080.0080.0080.0080.0080.0080.0080.0080.0080.0080.0080.0080.0080.0080.0080.0080.0080.00850.011ni0.02852.2390.61Rmean,Sigma,%density-2ndcomponent1.531.531.531.531.531.531.531.531.531.531.531.531.5281.521.511.421.421.421.421.452nrfor20wavelengths0.0050.0050.0050.0050.0050.0050.00530.0060.0060.00670.0070.0070.00880.01090.01890.02180.01950.06750.0460.004ni0.01182.00.22Rmean,Sigma,%density-3rdcomponent1.751.751.751.751.751.751.751.751.751.751.751.751.751.751.771.7911.7961.8081.8151.9nrfor20wavelengths0.4650.460.45880.45570.4530.45120.4470.440.4360.4350.4330.43060.430.4330.44960.46290.4720.4880.50.57ni1ResultswillbesavedintoFILE.mieUrban_IndustIdentificationoftheoutputfile(FILE).resultswillbesavedintoUrban_Indust.mie4.输入气溶胶浓度v能见度,以km为单位。输入具体的能见度值v=exp(-log(taer55/2.7628)/0.79902)0:下一行输入550nm气溶胶光学厚度-1:如果iaer=0,输入v=-15.目标海拔高度xpsxps=0.thepressureisgiveninmb输入以mb为单位的气压值xps0.meansyouknowthealtitudeofthetargetexpressedinkmandyouputthatvalueasxps6.传感器海拔高度xpp以km为单位xpp=-1000meansthatthesensorisaboardasatellitexpp=0meansthatthesensorisatthegroundlevel对于航空飞机来讲-100xpp0表示你知道传感器的高度,并以km为单位输入,这个高度是相对于目标的高度。对于航空飞机来讲,还需要输入puw,po3,taerp,分别为飞机与地表之间的水汽、臭氧含量
本文标题:6SV参数输入
链接地址:https://www.777doc.com/doc-2174170 .html