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

最新推荐文章于 2024-08-23 09:44:47 发布

qq_41627642

最新推荐文章于 2024-08-23 09:44:47 发布

阅读量2.7k

点赞数 3

分类专栏：遥感地信数据采集

本文链接：https://blog.csdn.net/qq_41627642/article/details/106187868

版权

遥感地信数据采集专栏收录该内容

14 篇文章 18 订阅

订阅专栏

论文1、

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

1、基本原理

在这里插入图片描述

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

modis_lut.pro创建6S查找表

;构建气溶胶反演6s辐射传输模型。

PRO MODIS_LUT;,month,day,iwave,idatm,iaer,lutpath,lutname
  igeom=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];；卫星天顶角   ；这里设置为0
  phiv=[0,24,48,72,96,120,144,168,180];；太阳方位角(卫星方位角为0，即相对方位角为O．180)；（作者卫星相对方位角为130）
  ; CD,lutpath
  cd,'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_lun
          PRINTF,lun,igeom
          PRINTF,lun,asol[c],phiv[e],avis[d],phi0,month,day
          PRINTF,lun,idatm
          PRINTF,lun,iaer
          PRINTF,lun,v
          PRINTF,lun,tao[b]
          PRINTF,lun,xps
          PRINTF,lun,xpp
          PRINTF,lun,iwave
          PRINTF,lun,inhomo
          PRINTF,lun,idirect
          PRINTF,lun,igroun
          PRINTF,lun,rapp
          FREE_LUN,lun
          SPAWN,'6s.exe<in.txt>out.txt',/hide ;调用6s <>符号为dos系统下的重定向符号 <从文件读取命令输入>将输出结果写入文件
          txtname='out.txt'
          OPENR,lun,txtname,/get_lun
          temp=STRARR(1,120)
          READF,lun,temp
          tt=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]
        ENDFOR
      ENDFOR
    ENDFOR
  ENDFOR
  ;endfor
  FREE_LUN,lutlun