IDL调用6S文件生成数据，及表观反射率有关二维图像

;;===============================================================================================

;;;   该程序是通过6s辐射校正模型生成查找表，以便进行气溶胶反演

;      该程序参数设置较为简单，如果要详细的参数，请自行设置

;      改进方法：1、读取每一个影像的角度信息建立查找表，这样精度更高

;                   2、利用IDL直接调用Fortran语言的6S源程序，并行生成查找表，速度提高20倍左右

;                   3、将所有参数都设置好，缩小步长，构建一个庞大的查找表，构建神经网络，一劳永逸

;;;==============================================================================================

PRO MODIS_LUT;month,day,iwave,idatm,iaer,lutpath,lutname

   igeom=0;自定义几何条件

   phi0=100;卫星方位角++++

   month=5;月份

   day=17;日期

   idatm=8;define own atmospheric model

   water=3.0;Water Vapor

   ozone=0.35;Ozone

   iaer=1;气溶胶模式大陆型

   v=0;能见度

   xps=-0.2;目标物高度

   xpp=-1000;星测

   iwave=42;自定义1输入波段范围和反射相函数42为modis的red

   inhomo=0;地表反射率均一地表

   idirect=0;无方向效应

   igroun=0;绿色植被

   surfre=[0.0,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5];地表反射率

   rapp=-1;无大气校正

   tao=[0.0,0.3,0.6,0.9,1.2,1.5,1.8,2.1,2.4,2.7,3.0];550nm气溶胶光学厚度

   asol=30;太阳天顶角

   avis=[0,7,14,21,28,35,42,49,56,63,70];卫星天顶角

   phiv=0;太阳方位角

   rf=FLTARR(11);表观反射率数组初始化

   ; CD,lutpath

   CD,'E:\，作业\大三下学期\遥感原理\6s_test\IDLdata';自己更改6s.exe所在的文件夹路径

   ;********************表观反射率随地表反射率（取值范围 0 ～ 0.5 ）的变化情况****************

   lutname='surf reflectanceoutput.txt'

   OPENW,lutlun,lutname,/get_lun

   ;设置循环过程

   FOR e=0,10 DO BEGIN;地表反射率

&nbs