6S parameters

6S Model ：
https://blog.csdn.net/sam92/article/details/50390556
https://www.jianshu.com/p/6a5c63534158
https://wenku.baidu.com/view/11856e7c8f9951e79b89680203d8ce2f01666551.html
MinGW install ：
https://blog.csdn.net/u012278016/article/details/101101486
https://blog.csdn.net/ZHAOJUNWEI08/article/details/86602120

PRO temp
  COMPILE_OPT idl2
  ; Geometric conditions 
  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]); month 
  DAY=uint(DATE[2]); Japan 
  ; Calculate the average elevation 
  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
  ; Center longitude and latitude 
  sLatitude = (minY+ maxY) / 2.0
  print,' Center latitude ',sLatitude
  ;  Atmospheric model type 
  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,' Atmospheric model ',idatm
  ; The spectral response function is defined by band 
  ; Enter the data type 
  ;-1<rapp<0 The input is the apparent reflectance ,rapp>0 The input is radiance 
  rapp=-0.5
  
  Parameter=[SOZ,SOA,SAZ,SAA,MONTH,DAY,paraID,xps,idatm,rapp]
  
END