本文共 2442 字，大约阅读时间需要 8 分钟。

6S模型：

PRO temp  COMPILE_OPT idl2  ;几何条件  Parameterxml='***********.xml'  SOZ=90-float(getGF1PMS_Parameters(filename=Parameterxml,nodeName2='SolarZenith'))  SOA=float(getGF1PMS_Parameters(filename=Parameterxml,nodeName2='SolarAzimuth'))  SAZ=90-float(getGF1PMS_Parameters(filename=Parameterxml,nodeName2='SatelliteZenith'))  SAA=float(getGF1PMS_Parameters(filename=Parameterxml,nodeName2='SatelliteAzimuth'))  DATE=getGF1PMS_Parameters(filename=Parameterxml,nodeName2='ReceiveTime')  DATE=strsplit(DATE,'-',/extract)  MONTH=uint(DATE[1]);月  DAY=uint(DATE[2]);日  ;计算平均高程  fid=ENVIRastertoFID(CalcuRaster)  ENVI_FILE_QUERY,fid,dims=dims,ns=ns,nl=nl,nb=nb  ENVI_CONVERT_FILE_COORDINATES,fid,0,0,minX,maxY,/to_map  ENVI_CONVERT_FILE_COORDINATES,fid,ns-1,nl-1,maxX,minY,/to_map  demraster=e.openraster('GMTED2010.jp2')  rasterSub=ENVISubsetraster(demraster,spatialref=demraster.spatialref,sub_rect=[minX,minY,maxX,maxY])  statistics=ENVIRasterStatistics(rasterSub)  ground_elevation=float(statistics['MEAN']/1000)  xps=-ground_elevation  ;中心经纬度  sLatitude = (minY+ maxY) / 2.0  print,'中心纬度',sLatitude  ; 大气模式类型  if sLatitude gt -10 and sLatitude le 20 then idatm=1  if sLatitude gt 20 and sLatitude le 30 then begin    if Month gt 05 and Month le 09 then begin      idatm=1    endif else idatm=2  endif  if sLatitude gt 30 and sLatitude le 40 then begin    if Month gt 05 and Month le 09 then begin      idatm=2    endif else idatm=4  endif  if sLatitude gt 40 and sLatitude le 50 then begin    if Month gt 03 and Month le 11 then begin      idatm=4    endif else idatm=3  endif  if sLatitude gt 50 and sLatitude le 60 then begin    if Month gt 05 and Month le 09 then begin      idatm=4    endif else idatm=3  endif  if sLatitude gt 60 and sLatitude le 70 then begin    if Month gt 03 and Month le 09 then begin      idatm=3    endif else idatm=5  endif  if sLatitude gt 70 and sLatitude le 80 then begin    if Month gt 05 and Month le 09 then begin      idatm=3    endif else idatm=5  endif  print,'大气模式',idatm  ;光谱响应函数按波段定义  ;输入数据类型  ;-1
   
    <0输入的是表观反射率,rapp>0输入的是辐射亮度  rapp=-0.5    Parameter=[SOZ,SOA,SAZ,SAA,MONTH,DAY,paraID,xps,idatm,rapp]  END
   

转载地址：https://blog.csdn.net/qq_33339770/article/details/106354384 如侵犯您的版权，请留言回复原文章的地址，我们会给您删除此文章，给您带来不便请您谅解！