python envi影像数据 envi compute statistics

ENVI5.3影像数据辐射定标及大气校正处理

数据准备

1）、选择合适影像数据，优先云量选择<1%的遥感数据 ；

2）、在ENVI中点击File—Open External File—Landsat—Geotiff with Metadata中只需打开***********_MTL.txt即可打开所有波段数据（除band6）；

辐射定标

• ENVI自带辐射校正模块
  1、选择Basic Tools—Preprocessing—Calibration Utilities—Landsat Calibration；
  
  2、参数选择：从MTL.txt头文件中读取DATE_ACQUIRED、Sun Elevation（按照以上步骤操作，ENVI已自动读取Sun Elevation数据），Calibration Type选择Radiance；

大气校正

简单一点的大气校正可以采用ENVI的FLAASH模块，以下是FLAASH操作的步骤：

1. FLAASH 模块的进入方法是Spectral–FLAASH，或者是Basic Tools–Preprocessing–Calibration Utilities–FLAASH。
2. FLAASH 模块的操作界面分为三块：最上部设定输入输出文件；中间设定传感器的参数；下部设定大气参数。
4. 首先设定输入输出文件。FLAASH 模块要求输入辐亮度图像，输出反射率图像。之前我们进行了辐射定标，得到辐亮度图像，在这里要把BSQ 格式的图像转换为BIL 或者BIP 格式的图像，然后再Input Radiance Image 中选择转换格式后的图像。操作步骤：(Basic Tools–Convert Data(BSQ,BIL,BIP))。这里注意，当输入图像后，程序会让你选择Scale Factor，即原始辐亮度单位与ENVI 默认辐亮度单位之间的比例，二者之间转换的比例是10，因此在下图中选择Single scale factor，填写10.000。



5. 此外，如果TM 影像的头文件中没有波段的信息，在这里也要求你提供一个.txt 文件以包含此信息。那么，准备好一个.txt 文件，其中含有一列TM 每个波段中心波长的信息。如果，我们打开的是Geotiff with meta，就不用填写波段信息了。
6. 在Output Reflectance File 和Output Directory for FLAASH files 里面设定输出文件的文件名和位置。
7. 设定传感器参数。首先是Scene Center Location，即遥感图像中心的坐标，以及Flight Date, Flight Time GMT，这三者都可以在TM 的头文件中找到，填入即可。
   （1）、遥感图像中心的坐标：如果图像有地理坐标则自动获取
   （2）、影像成像时间（格林威治时间）：在layer manager中的数据图层中右键选择View Metadata，浏览time字段获取成像时间，1985年10月30号02:05:45，也可以从元文件"*_MTL.txt"中找到，具体名称：DATE_ACQUIRED 和SCENE_CENTER_TIME。
8. 在Sensor Type 菜单中根据影像数据源选择，此处填Landsat TM5，此时Sensor altitude 自动填上为705km，而Pixel Size 填为30m。
10. 根据遥感影像研究区实际情况，填写Ground Elevation。
    *Ground Elevation获取方法：
    （1）、打开需要统计区域对应的图像。
    （2）、选择File->Open World Data ->Elevation(GMTED2010),打开ENVI自带全球900米分辨率的DEM数据。
    注：如果其他版本，需要打开其他DEM数据，如90米的srtm或者30米的G-DEM数据。
    （3）、在Toolbox中，选择/Statistics/Compute Statistics，打开Compute Statistics输入文件对话框，选择GMTED2010.jp2数据。单击Stats Subset按钮，打开Select Statistics Subset对话框。



11. （4）、在Select Statistics Subset对话框中，单击File按钮，选择统计区域对应的图像，单击Ok。



12. （5）、在统计面板中，默认选择Basic Stats基本统计即可，这样就可以得到平均海拔高程36米。



13. 
    



15. 大气参数数据设置：
    a) 、Atmospheric Model( 大气模式): 共有Sub-Arctic Winter (SAW) ，Mid-Latitude Winter (MLW)，U.S. Standard (US) ，Sub-Arctic Summer(SAS), Mid-Latitude Summer (MLS) 和Tropical (T) 。根据经纬度和时间可以选定研究区的大气模式，见下图。



16. b)、 Aerosol Model（气溶胶模式）：有Rural, Urban, Maritime 和Tropospheric四种选择。根据实际情况选择即可。关于此四种模式的解释见ENVI Help。
    c) 、当我们选择TM 时，可选的参数还有Aerosol Retrieval 和Initial Visibility。
    ·Aerosol Retrieval：选择2-Band(K-T)方法；
    ·Initial Visibility：初始能见度Initial Visibility只有在气溶胶反演方法为None时候，以及K-T方法在没有找到黑暗像元的情况下。



17. FLAASH基本参数数据设置成图：



19. 多光谱数据参数设置：
    当基本设置里设置了水气反演模型和气溶胶模型时，相应的在多光谱设置框中设置参数：
    （1）、 单击Multispectral Settings，打开多光谱设置面板；
    （2）、 K-T反演选择默认模式：Defaults->Over-Land Retrieval standard（600:2100），自动选择对应的波段；
    ·设置值表格：



20. （3）、其他参数选择默认。
    多光谱参数数据设置面板图：



22. 高级设置
    单击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中间文件，便于诊断运行过程中的错误。



23. 如果对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度。
24. 输出文件及结果浏览
    (1)、设置好参数后，单击Apply执行大气校正；
    (2)、完成后会得到反演的能见度和水汽柱含量。
    (3)、显示大气校正结果图像，查看像元值，可以看到像元值扩大10000倍后，值在几百到几千不等。如果要得到0-1范围内的反射率数据，可以使用BandMath除以10000.0。
    注：以上内容通过网上各博主内容整合而成。