Landsat气溶胶反演论文整理（以后知识点实践）

论文1、

Ｌａｎｄｓａｔ８ＯＬＩ城市地区气溶胶光学厚度遥感反演
主要技术：
摘要：在地表反射率较高、结构复杂的城市地区，传统的浓密植被气溶胶反演算法难以适用。通过分析地物波谱库中的植被和土壤波谱信息，模拟建立归一化植被指数（ＮＤＶＩ）与红、蓝波段地表反射率之间的相关关系，提出使用ＭＯＤＩＳ植被指数产品（ＭＯＤ１２）确定地表反射率的方法，实现该类型区域气溶胶光学厚度（ＡＯＤ）反演。以Ｌａｎｄｓａｔ８ＯＬＩ数据为例，选取北京市为研究区进行反演实验，使用ＡＥＲＯＮＥＴ地基观测数据与ＭＯＤＩＳ气溶胶产品（ＭＯＤ０４）对反演结果进行验证。结果表明，当反射率较高时，ＮＤＶＩ与红、蓝波段地表反射率仍存在较高的相关性，利用该指数能够准确获取高反射率地区的地表信息，算法反演结果与实测值具有较好的一致性，总体相关系数达０．９６６，６８％的反演结果满足误差精度要求，当ＡＯＤ＜０．５时，有８２．３％的结果满足精度要求，较ＭＯＤ０４精度有了较大改善。

1、基本原理

（1）计算不同植被覆盖度下波长处混合地物反射率数据；
ENVI波普运算
ENVI编写函数进行波段运算
（2）计算红光、蓝光、近红外波段的反射率数据；

modis_lut.pro创建6S查找表

;构建气溶胶反演6s辐射传输模型。PRO MODIS_LUT;,month,day,iwave,idatm,iaer,lutpath,lutnameigeom=0;自定义几何条件phi0=0;卫星方位角++++month=1;；月份day=13;；日期idatm=2;：大气模式中纬度夏季，这个参数不用修改iaer=1;：气溶胶模式大陆型，这个参数不用修改v=0;；能见度xps=0;；目标物高度xpp=-1000;；星测iwave=42;：自定义1输入波段范围和反射相函数42为modis的red;采用的地表反射率波段，红蓝分别为0~0.25，0~0.15，步长均为1.inhomo=0;；地表反射率均一地表idirect=0;；无方向效应igroun=1;：绿色植被rapp=-2;：无大气校正tao=[0.0001,0.25,0.5,1.0,1.5,1.95];；550nm气溶胶光学厚度 ；这里也要修改asol=[0,12,24,36,48,60];；太阳天顶角  ；这里市[20,50],步长间隔为5，7各参数avis=[0,12,24,36,48,60];；卫星天顶角   ；这里设置为0phiv=[0,24,48,72,96,120,144,168,180];；太阳方位角(卫星方位角为0，即相对方位角为O．180)；（作者卫星相对方位角为130）; CD,lutpathcd,'C:\Users\ljkin\Desktop\6S\6S 程序';自己更改6s.exe所在的文件夹路径；这里要修改lutname='modis_lut.txt'  ;这里要修改OPENW,lutlun,lutname,/get_lun;设置循环过程;for a=0,2 do begin;蓝红两个通道FOR b=0,5 DO BEGIN;550nm气溶胶光学厚度FOR c=0,5 DO BEGIN;太阳天顶角FOR d=0,5 DO BEGIN;；卫星天顶角FOR e=0,8 DO BEGIN;；太阳方位角(卫星方位角为0，即相对方位角为O．180)txtname='in.txt'OPENW,lun,txtname,/get_lunPRINTF,lun,igeomPRINTF,lun,asol[c],phiv[e],avis[d],phi0,month,dayPRINTF,lun,idatmPRINTF,lun,iaerPRINTF,lun,vPRINTF,lun,tao[b]PRINTF,lun,xpsPRINTF,lun,xppPRINTF,lun,iwavePRINTF,lun,inhomoPRINTF,lun,idirectPRINTF,lun,igrounPRINTF,lun,rappFREE_LUN,lunSPAWN,'6s.exe<in.txt>out.txt',/hide ;调用6s <>符号为dos系统下的重定向符号 <从文件读取命令输入>将输出结果写入文件txtname='out.txt'OPENR,lun,txtname,/get_luntemp=STRARR(1,120)READF,lun,temptt=STRMID(temp[0,105],61,8)sA=STRMID(temp[0,111],61,8)rou=STRMID(temp[0,114],61,8)FREE_LUN,lun;依次为辐射方程中的P T S参数、太阳天顶角，卫星天顶角，相对方位角、气溶胶光学厚度PRINTF,lutlun,sA,tt,rou,asol[c],avis[d],phiv[e],tao[b]ENDFORENDFORENDFORENDFOR;endforFREE_LUN,lutlun

2018Esri开发竞赛ENVI-IDL组作品欣赏-Landsat8城市生态环境监测
[2012Esri开发竞赛ENVI/IDL组作品欣赏——“北京市气溶胶反演与发布系统”]
(http://blog.sina.com.cn/s/blog_764b1e9d010143nv.html)

论文2

利用 Landsat/TM 影像监测北京地区气溶胶的空间分布
遥感影像中的DN值、辐射亮度、表观反射率、地表反射率、亮度温度粗解
Landsat8数据打开和辐射定标处理（怎么样得到表观反射率）
表观辐亮度转化为表观反射率