当前位置:首页 > 遥感 > 正文内容

ENVI大气校正

admin4年前 (2020-04-19)遥感2719

每次找都不方便,直接搬过来 

原地址:http://blog.sina.com.cn/s/blog_764b1e9d0102v59e.html


目录

辐射定标和大气校正    

1. 概述    

2.    辐射定标    

3.    多光谱数据FLAASH大气校正    

3.1    FLAASH输入数据要求 

3.2    详细操作步骤    

3.2.1    输入文件准备    

3.2.2    基本参数设置    

3.2.3    多光谱数据参数设置    

3.2.4    高级设置    

3.2.5    处理结果浏览    

4.    高光谱数据FLAASH大气校正    

4.1    输入文件准备    

4.2    基本参数设置    

4.3    高光谱数据参数设置    

4.4    高级设置    

4.5    处理结果浏览    

5.    快速大气校正(QUAC)    

  

 

1. 概述

        太阳辐射通过大气以某种方式入射到物体表面然后再反射回传感器,由于大气气溶胶、地形和邻近地物等影像,使得原始影像包含物体表面,大气,以及太阳的信息等信息的综合。如果我们想要了解某一物体表面的光谱属性,我们必须将它的反射信息从大气和太阳的信息中分离出来,这就需要进行大气校正过程。

        本课程学习利用ENVI 大气校正扩展模块(FLAASH和QUAC)对多光谱和高光谱数据进行大气校正的过程,还包括大气校正的准备工作——辐射定标。

        FLAASH是基于MODTRAN5辐射传输模型,MODTRAN模型是由进行大气校正算法研究的领先者SpectralSpectral Sciences, Inc和美国空军实验室(Air Force Research Laboratory)共同研发。Exelis VIS公司负责集成和GUI设计。

        注:本课程需要大气校正模块(Atmospheric Correction)扩展模块的使用许可。

 

2. 辐射定标

        一般来讲,辐射定标就是将图像的数字量化值(DN)转化为辐射亮度值或者反射率或者表面温度等物理量的处理过程。辐射定标参数一般存放在元数据文件中,ENVI中的通用辐射定标工具(Radiometric Calibration)能自动从元数据文件中读取参数(如表2.1所示),从而完成辐射定标。

        注:这里以Landsat5 L1G级的数据为例,由于数据较大,并且该工具操作简单,这里不提供例子数据。

(1) 选择File>Open As>Landsat>GeoTIFF with Metadata,选择打开*_MTL.txt文件。

表2.1 Radiometric Calibration工具自动识别的数据类型


传感器

可选定标类型

打开的元数据文件

辐射亮度值

大气表观反射率

亮温

ALOS AVNIR-2 、 PRISM Level-1B2


HDR*.txt

DMC DIMAP


.dim

EO-1 ALI


使用File > Open As > EO-1 > HDF 菜单打开 *_HDF.L1G文件,*_MTL.L1G 文件必须在同一目录中。

EO-1 Hyperion


使用 File > Open As > EO1 > HDF 菜单选择 .L1R 文件.

GeoEye-1


.til

IKONOS


metadata.txt

KOMPSAT-3



*_aux.xml

Landsat TM, ETM+, Landsat-8 OLI/TIRS

*_MTL.txt, *WO.txt, *.met

OrbView-3


.pvl

Pleiades Primary or Ortho (single 或者 mosaic)


DIM*.xml

QuickBird


.til

RapidEye Level-1B 


*_metadata.xml

ResourceSat-2



.h5

SPOT DIMAP


METADATA.DIM

SSOT DIMAP



METADATA.DIM

WorldView-1 and 2


.til



 注:以下类型的数据,ENVI打开后自动定标。

      Göktürk-2: Radiance       

  • MODIS/ASTER Simulator (MASTER): Radiance
  • RASAT: Radiance
  • ASTER Level-1B
  • MODIS Level-1B through 4 (HDF-EOS)

    (2) 在Toolbox中,选择Radiometric Correction > Radiometric Calibration,对文件对话框中选择多光谱数据。打开Radiometric Calibration面板。

    (3) 在Radiometric Calibration面板中,设置以下参数:

    (4) 设置输出路径和单位名,单击OK执行辐射定标。

    • 定标类型(Calibration Type):辐射率数据Radiance
    • 单击Apply FLAASH Settings按钮,自动设置FLAASH大气校正工具需要的数据类型,包括储存顺序(Interleave):BIL或者BIP;数据类型(Data Type):Float;辐射率数据单位调整系数(Scale Factor):0.1。

图1.jpg

图:Radiometric Calibration面板

(5) 显示辐射定标结果图像,选择Display>Profiles>Spectral查看波谱曲线,看到定标后的数值主要集中在0-10范围内,单位是µW/(cm2 * sr * nm)

图2.jpg

图:辐射定标结果的波谱曲线

 

3. 多光谱数据FLAASH大气校正

 

3.1 FLAASH输入数据要求

        FLAASH对大气校正的输入图像做了一些要求,具体要求如下:

一、图像基本参数

        波段范围:卫星图像:400-2500nm,航空图像:860nm-1135nm。

        如果要执行水汽反演,光谱分辨率<=15nm,且至少包含以下波段范围中的一个:

  • 1050-1210 nm
  • 770-870 nm
  • 870-1020 nm

        像元值类型:经过定标后的辐射亮度(辐射率)数据,单位是:(μW)/(cm2*nm*sr)。

二、数据储存类型

        数据类型:浮点型(Floating Point)、32位无符号整型(Long Integer)、16位无符号和有符号整型(Integer、Unsigned Int)。

        文件类型: ENVI标准栅格格式文件,BIP或者BIL储存结构。

三、辅助信息

        中心波长:数据头文件中(或者单独的一个文本文件)包含中心波长(wavelenth)值,如果是高光谱还必须有波段宽度(FWHM),这两个参数都可以通过编辑头文件信息输入(Edit Header)。

        波谱滤波函数(波谱响应函数)文件:对于未知多光谱传感器(UNKNOWN-MSI)需要提供波谱滤波函数文件。

 

3.2 详细操作步骤

        ENVI大气校正模块的使用主要又以下7个方面组成:1、输入文件准备,2、基本参数设置,3、多光谱数据参数设置,4、高光谱数据参数设置,5、高级设置,6、输出文件,7、处理结果。下面介绍这7个方面内容。

        下面以Landsat5的多光谱数据为例介绍操作步骤。

 

表3.1 练习数据信息

文件

说明

Landsat5_beijing_radiance.dat

经过Radiometric Calibration工具定标的Landsat5辐射率图像,同时空间子区裁剪了部分区域。

Landsat5_beijing_radiance.hdr

头文件

Landsat5_template.txt

本例子FLAASH参数设置工程文件

 

3.2.1 输入文件准备

        根据前一章节中介绍的FLAASH对数据的要求准备待校正文件。由于使用了Radiometric Calibration工具辐射定标,数据类型、储存顺序、辐射率数据单位都符合FLAASH要求,Landsat5的L1G级数据包括了中心波长信息。

 

3.2.2 基本参数设置

        在Toolbox中打开FLAASH工具:/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction。启动FLAASH Atmospheric Correction Module Input Parameters面板

(1) Input Radiance Image:选择辐射定标结果数据,在打开的Radiance Scale Factors面板中,设置Single scale factor:1。

注:辐射率数据的单位已经是:(μW)/(cm2*nm*sr)

图3.jpg

图:Radiance Scale Factors面板

(2) Output Reflectance File:设置输出路径和文件名;

(3) Output Directory for FLAASH Files:设置其他文件输出目录;

(4) 传感器基本参数设置:

  • 中心点经纬度Scene Center Location:如果图像有地理坐标则自动获取;
  • 选择传感器类型Sensor Type:Landsat TM5,其对应的传感器高度以及影像数据的分辨率自动读取;
  • 设置影像区域的平均地面高程Ground Elevation:0.05km
  • 影像成像时间(格林威治时间):在layer manager中的数据图层中右键选择View Metadata,浏览time字段获取成像时间,2009年9月22号02:43:22。

        注:也可以从元文件"*_MTL.txt"中找到,具体名称:DATE_ACQUIRED 和SCENE_CENTER_TIME;

(5) 大气模型参数选择Atmospheric Model:Mid-Latitude Summer(根据成像时间和纬度信息依据下表规则选择);

表3.2数据经纬度与获取时间对应的大气模型

Latitude (°N)

Jan.

March

May

July

Sept.

Nov.

80

SAW

SAW

SAW

MLW

MLW

SAW

70

SAW

SAW

MLW

MLW

MLW

SAW

60

MLW

MLW

MLW

SAS

SAS

MLW

50

MLW

MLW

SAS

SAS

SAS

SAS

40

SAS

SAS

SAS

MLS

MLS

SAS

30

MLS

MLS

MLS

T

T

MLS

20

T

T

T

T

T

T

10

T

T

T

T

T

T

0

T

T

T

T

T

T

-10

T

T

T

T

T

T

-20

T

T

T

MLS

MLS

T

-30

MLS

MLS

MLS

MLS

MLS

MLS

-40

SAS

SAS

SAS

SAS

SAS

SAS

-50

SAS

SAS

SAS

MLW

MLW

SAS

-60

MLW

MLW

MLW

MLW

MLW

MLW

-70

MLW

MLW

MLW

MLW

MLW

MLW

-80

MLW

MLW

MLW

SAW

MLW

MLW

(6) 气溶胶模型Aerosol Model:Urban;

(7) 气溶胶反演方法Aerosol Retrieval:2-band(K-T);

注:初始能见度Initial Visibility只有在气溶胶反演方法为None时候,以及K-T方法在没有找到黑暗像元的情况下。

(8) 其他参数按照默认设置即可。

图4.jpg

图:FLAASH基本参数设置

 

3.2.3 多光谱数据参数设置

(1) 单击Multispectral Settings,打开多光谱设置面板;

(2) K-T反演选择默认模式:Defaults->Over-Land Retrieval standard(600:2100),自动选择对应的波段;

(3) 其他参数选择默认。

图5.jpg

图:多光谱设置面板

 

3.2.4 高级设置

        单击Advanced Settings打开高级设置面板。这里一般选择默认设置能符合绝大部分数据情况,在右边面板中设置:

(1) 分块处理(Use Tiled Processing):是否分块处理,选择Yes能获得较快的处理速度,Tile Size一般设为4-200m,根据内存大小设置,这里设置为100m(计算机物理内存8G)。

(2) 空间子集(Spatial Subset):可以设置输出的空间子集,这里选择默认输出全景。

(3) 重定义缩放比例系数(Re-define Scale Factors For Radiance Image):重新选择辐射亮度值单位转换系数,这里不设置。

(4) 输出反射率缩放系数(Output Reflectance Scale Factor):为了降低结果储存空间,默认反射率乘于10000,输出反射率范围变成0~10000。

(5) 自动储存工程文件(Automatically Save Template File):选择是否自动保存工程文件。

(6) 输出诊断文件(Output Diagnostic Files):选择是否输出FLAASH中间文件,便于诊断运行过程中的错误。

图6.jpg

图:高级设置面板

        如果对Modtran模型非常熟悉,可根据数据情况进行调整,如下为其余部分的参数说明。

(1) 气溶胶厚度系数(Aerosol Scale Height):用于计算邻域效应范围。一般值为1~2km,默认为1.5km。

(2) CO2混合比率(CO2 Mixing Ratio):默认为390ppm,它是依据2001测量值为370ppm,增加20ppm以得到更好的结果。

(3) Use Square Slit Function:No

(4) 使用领域纠正(Use Adjacency Correction):Yes 或者No。

(5) 使用以前的MODTRAN模型计算结果(Reuse MODTRAN Calculations):

  • No:重新计算MODRTRAN辐射传输模型。
  • Yes:执行上一次FLAASH运行获得的MODRTRAN辐射传输模型,每次运行FLAASH后,都会在根目录和临时文件夹下生成一个acc_modroot.fla。

       (6) MODTRAN模型的光谱分辨率(Modtran Resolution):越低分辨率具有较快速度而相对较低的精度,主要影响区域在2000 nm附近。高光谱数据默认为5 cm-1,多光谱数据默认为15 cm-1。

(7) MODTRAN多散射模型(Modtran Multiscatter Model):校正大气散射对成像的影响,提供三种模型供选择ISAACS,DISORT和 Scaled DISORT。默认是Scaled DISORT和streams为8。

  • Isaacs模型计算速度快,精度一般;
  • DISORT模型对于短波(小于1000nm)具有较高的精度,但是速度非常比较慢,由于散射对短波(如可见光)影响较大,长波(近红外以上)影响较小,因此当薄雾较大和短波图像时可以选择此方法;
  • Scaled DISORT提供在大气窗口内与DISORT类似的精度,速度与Isaacs类似,这模型是推荐使用的模型。当选择DISORT 或者 Scaled DISORT,需要选择streams:2、4、8、16,这个值是用来估算散射的方向,可见streams值越大速度越慢。

        (8) 观测参数

  • 天顶角(Zenith Angle):是传感器直线视线方向和天顶的夹角,范围是90~180度,其中180为传感器垂直观测。
  • 方位角(Azimuth Angle):范围是-180~180度。

 

3.2.5 处理结果浏览

(1) 设置好参数后,单击Apply执行大气校正;

(2) 完成后会得到反演的能见度和水汽柱含量。

(3) 显示大气校正结果图像,查看像元值,可以看到像元值扩大10000倍后,值在几百到几千不等。如果要得到0-1范围内的反射率数据,可以使用BandMath除以10000.0。

(4) 选择Display>Profiles>Spectral查看典型地物波谱曲线,如植被、水体等。

图7.jpg图8.jpg

图:FLAASH大气校正结果中获取的波谱曲线(左-植被,右-水体)

 

4. 高光谱数据FLAASH大气校正

        下面以AVIRIS高光谱传感器数据为例介绍操作步骤。

表4.1 练习数据信息

文件

说明

CupriteAVIRISSubset.dat

AVIRIS高光谱传感器224个波段的辐射率图像,ENVI栅格文件格式。

CupriteAVIRISSubset.hdr

头文件

AVIRIS11_gain.txt

用于重新调整辐射率数据单位的文本文件

AVIRIS_template.txt

本例子FLAASH参数设置工程文件

 

4.1 输入文件准备

        本例中的AVIRIS高光谱数据已经是辐射率数据,包括了中心波长、波段宽度(FWHM),Integer数据类型,BIP储存顺序。在layer manager中的数据图层中右键选择View Metadata,可以浏览所有的信息。

注:高光谱数据一般有一套预处理工具,可以满足FLAASH工具的数据要求。

 

4.2 基本参数设置

        打开CupriteAVIRISSubset.dat数据,在Toolbox中打开FLAASH工具:/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction。启动FLAASH Atmospheric Correction Module Input Parameters面板。

(1) Input Radiance Image:选择CupriteAVIRISSubset.dat文件,在打开的Radiance Scale Factors面板中,选择默认Read array of scale factors (1 per band) from ASCII file。

(2) 在对话框中选择AVIRIS11_gain.txt文件,在Input ASCII File对话框中,Scale Column设置1。

注: AVIRIS辐射率原始的定标结果是浮动型,为了减少储存空间,NASA/JPL 使用一组比例系数将浮点型数据变成整型数据。设定比例系数将数据由变成。在使用的时候需要将整型变回浮点型,单位变成µW/(cm2 * sr * nm)。这一步就是完成这个,文本文件中第一列就是缩放系数。

图9.jpg图10.jpg

图:Radiance Scale Factors和输入文本文件面板

(3) Output Reflectance File:设置输出路径和文件名;

(4) Output Directory for FLAASH Files:设置其他文件输出目录;

(5) 传感器基本参数设置:

  • 中心点经纬度Scene Center Location:如果图像有地理坐标则自动获取;
    • 选择传感器类型Sensor Type:Hyperspectral > AVIRIS;
    • 传感器飞行高度Sensor Altitude (km):20
  • 设置影像区域的平均地面高程Ground Elevation:0.6km
  • 图像地面分辨率Pixel Size (m):15.5
    • 影像成像时间(格林威治时间)Flight Date:Aug 6, 2011,Flight Time (GMT): 19:20:00

        注:大部分图像成像信息存放在*.info文件中。

(6) 大气模型参数选择Atmospheric Model:U.S. Standard(根据成像时间和纬度信息依据下表规则选择);

(7) 水汽反演Water Retrieval:yes

(8) 水汽吸收波长Water Absorption Feature:1135

(9) 气溶胶模型Aerosol Model:Urban;

(10) 气溶胶反演方法Aerosol Retrieval:2-band(K-T);

注:初始能见度Initial Visibility只有在气溶胶反演方法为None时候,以及K-T方法在没有找到黑暗像元的情况下。

(11) 其他参数按照默认设置即可。

图11.jpg

图:FLAASH基本参数设置

 

4.3 高光谱数据参数设置

(1) 单击Hyperspectral Settings,打开高光谱设置面板;

(2) 其他参数选择默认。

图12.jpg

图:高光谱设置面板

 

4.4 高级设置

        单击Advanced Settings打开高级设置面板。这里一般选择默认设置能符合绝大部分数据情况,在右边面板中设置:

(1) 分块处理(Use Tiled Processing):是否分块处理,选择Yes能获得较快的处理速度,Tile Size一般设为4-200m,根据内存大小设置,这里设置为100m(计算机物理内存8G)。

(2) 其他默认设置。

图13.jpg

图:高级设置面板

 

4.5 处理结果浏览

(1) 设置好参数后,单击Apply执行大气校正;

(2) 完成后会得到反演的能见度和水汽柱含量。

图14.jpg

图:FLAASH结果

(3) 显示大气校正结果图像,像元值同样扩大了10000倍。

(4) 选择Display>Profiles>Spectral查看地物波谱曲线。

注:得到的波谱曲线有几个地方断开了,这是由于FLAASH会根据反射信号的强度判断部分波段属于"Bad band",在头文件中标识这些波段。这里标识的波段包括108-112和155-166范围内的波段。

(5) 在工具栏中的Go to文本框中输入326,290,回车。得到这个像素点的波谱曲线,如下图所示。

(6) 在1300 nm 和 1500 nm除有两个低谷,这个主要由于水汽吸收造成反射率很低。这种情况可以利用"edit header"工具将这几个波段标识为"Bad band"。

图15.jpg

图:FLAASH大气校正结果中获取(326,290)点的波谱曲线

(7) 高光谱FLAASH大气校正结果还会产生两个文件:水汽含量反演图像(water.dat)和云层掩膜图像(cloudmask.dat)。

 

5. 快速大气校正(QUAC)

        快速大气校正工具(Quick Atmospheric Correction简称QUAC)自动从图像上收集不同物质的波谱信息,获取经验值完成高光谱和多光谱的快速大气校正(如下图)。它得到结果的精度近似FLAASH或者其他基于辐射传输模型的+/-15%。它支持的多光谱和高光谱波谱范围是(0.4~2.5 μm)。

图16.jpg

图:快速大气校正流程图

        QUAC的输入数据可以是辐射亮度值、表观反射率、无单位的raw数据。可以是任何数据储存顺序       (BIL/BIP/BSQ)和储存类型,多光谱和高光谱传感器数据的每个波段必须有中心波长信息。

        QUAC的操作非常简单,如下:

(1) 在Tool box中,启动/Radiometric Correction/Atmospheric Correction Module/QUick Atmospheric Correction (QUAC),在文件输入对话框中选择校正的图像文件;

注:如果待校正图像是不规则的,包括较多的背景(如0值),则需要在选择图像文件对话框中,选择/生成一个掩膜文件。

图17.jpg

图:文件选择对话框

(2) 打开QUick Atmospheric Correction Parameters面板,在Sensor Type中选择相应的传感器类型,一般会根据图像信息自动选择。

(3) 选择文件名和路径输出。

注:QUAC大气校正的结果同样是扩大了10000倍的反射率数据。

图18.jpg

图:QUick Atmospheric Correction Parameters面板

(尚未失效的地址)

数据下载:http://pan.baidu.com/s/1dDeUOl7   http://pan.baidu.com/s/1jG64UwQ

视频下载:http://pan.baidu.com/s/1pJC0xwJ

                 http://pan.baidu.com/s/1gd7WXk3

                 http://pan.baidu.com/s/14JvB4


扫描二维码推送至手机访问。

版权声明:本文由lovedm.club发布,如需转载请注明出处。

本文链接:http://lovedm.club/?id=39

分享给朋友:

“ENVI大气校正” 的相关文章

AVIRIS数据下载与预处理相关

AVIRIS数据下载与预处理相关

AVIRIS数据下载网站 https://aviris.jpl.nasa.gov/dataportal/ ,需要科学上网,数据量较大,折腾了好长时间总算是搞下来一景。处理可参考这个网站 https://www.harrisgeospatial.com/docs/Prep...

太阳天顶角、太阳高度角、太阳方位角

太阳天顶角、太阳高度角、太阳方位角

1.天顶角(Zenith Angle)天顶角指光线入射方向和天顶方向的夹角,有个英文解释地挺好的:The zenith angle is the angle between the zenith line(pointing straight up) and the direction to the...

ENVI计算光谱一阶微分

ENVI计算光谱一阶微分

首先在Display菜单中选择Spectral Library Viewer,如图:再选择Import,导入ASCII文件,选择文件后,按照自己的数据选择以下选项:点击OK后,波谱曲线会绘制出来,如下:在ToolBox中找到Spectral中的Spectral Math工具,双击打开工具,工具中输入...

GDAL-3.0.1+PROJ-6.2.1+SQLite+HDF4+HDF5+VS2019编译

这次上传百度网盘了。 http://lovedm.club/?id=40 不要用这个了 之前的没有加上hdf文件的读取,又重新编译的。链接如下: https://pan.baidu.com/s/154blFtGE6IIcUz_b1bk4xQ 提取码:tzx8此站也可下载https://www...

利用GDAL进行波段合成

顺手写的一个小工具,给出头文件和实现 主要用了GDALDataset类下面的RasterIO函数,没有用GDALBand下面的RasterIO 另外实现了一个定时器类,实际上就是之前写的一个http://lovedm.club/?id=98 输出文件格式为固定的GTiff,函数的第三个参数是G...

Python embed版本安装使用

Python embed版本安装使用

最近一个项目需要使用Python脚本后台处理数据,但是用到了几个第三方包,客户使用时部署环境比较麻烦,查找资料后可以使用Python的embed版本,也就是嵌入式版本。好处是所有用到的第三方库直接安装到目录下,部署时直接拷贝走文件夹即可,无需其它配置。下面记录一下安装使用流程。 下载Pytho...