ENVI处理高分一号PMS影像数据

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档

---

一、多光谱影像

1、数据打开

通过菜单栏File -> Open As-> China Satellites -> GF-1，在弹出的对话框，选择GF1_PMS2_E104.0_N36.0_20140724_L1A0000284766-MSS2.xml文件打开数据。

2、辐射定标

1)在Toolbox中，选择Radiometric Correction -> Radiometric Calibration，在弹出的File Selection对话框中，选择GF1_PMS2_E104.0_N36.0_20140724_L1A0000284766-MSS2.tiff，点击OK；

2)在Radiometric Calibration参数设置面板中，直接点击Apply FLAASH Settings按钮，其余参数会自动调整为FLAASH大气校正所要求的数据输入格式；

3)设置辐射定标结果输出路径和文件名；

4)点击OK开始执行

5)右下角显示进度，定标结束后结果会自动加载到Data Manager中。

6）结果

2、大气校正

1)在Toolbox中，打开/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction工具，弹出FLAASH Atmospheric Correction Model Input Parameters参数设置面板。


2)在基本参数设置面板中，上面部分主要用于设置数据输入输出，中间部分主要用于设置影像和传感器相关参数，最下部为大气模型及气溶胶反演相关参数设置。

上面部分：

Input Radiance Image：选择上一步辐射定标后的结果；

Radiance Scale Factors：选择Use single scale factor for all bands，数值保持默认1

Output Reflectance File：设置经大气校正后的地表反射率数据输出路径及文件名；

Output Directory for FLAASH Files：校正过程中生成其他文件的存储路径；

Rootname for FLAASH Files：输出文件名前缀，可不填。

至此，上面部分参数已经设置完毕
中间部分主要用于设置影像和传感器相关参数：

Scene Center Location：影像中心经纬度；

Sensor Type：传感器类型，保持默认GF-1；

Sensor Altitude(km)：传感器高度，645km(ENVI默认不会自带填入，手动输入即可)；

Ground Elevation(km)：影像对应区域地面平均高程。此处输入2，注意单位是km；

Pixel Size(m)：像元大小，8m；

Flight Date：数据获取时间，可以通过文件名或查看元数据信息获取；

中间部分已经设置完毕。

最下部为大气模型及气溶胶反演相关参数设置：

Atmospheric Model：大气模型，一般根据影像中心纬度和获取月份确定，需借助帮助文档完成。

这里选Tropical

Water Retrieval：此处保持默认NO；

Aerosol Model：气溶胶模型，有Rural、Urban、Maritime和Tropospheric四个选项可选。观察影像可以发现影像位于郊区，此处选择Rural；

Aerosol Retrieval：气溶胶反演方法，使用暗像元反射比模型估算影像气溶胶含量和平均能见度，有None、2-Band(K-T)和2-Band Over Water三个选项可选。此处选择None；

Initial Visibility(km)：初始能见度。根据影像获取时大气情况设置，如果气溶胶无法反演时，该值将作为初始值参与大气校正，此处保持默认即可；

在基本参数设置面板中，点击右下角Advanced settings…，打开高级参数设置面板；

左方参数框中的参数一般保持默认即可，对于右侧参数框中的参数：

Use Tiled Processing：是否采用分块处理。默认为Yes，推荐进行分块，分块大小(Tile Size)可根据计算机内存情况确定，可设置为安装内存的75%，默认为100M或者大于Classic中设置的缓存大小(Cache Size)。对于包含许多0值的影像，分块大小不易设置太小，避免出现分块像元值全为0而报错。此处保持默认100；

其他参数保持默认即可，参数设置完毕后点击OK。

5)在基本参数设置面板中，点击左下角Apply按钮，弹出FLAASH Atmospheric Correction面板，显示处理进度。处理结束后，会弹出一个简单的统计面板，说明大气校正结束

至此，大气校正结束。校正后的地表反射率数据存放在设置的Output Reflectance File路径下，其他结果文件存放在设置的Output Directory for FLAASH Files路径下。

3、正射校正

1）点击Toolbox→Geometric Correction→Orthorectification→RPC Orthorectification Workflow，打开正射校正流程化工具；

2）在File Selection面板中，Input File选择经过大气校正的多光谱数据（GF1_PMS2_E104.0_N36.0_20140724_MSS_reflect），DEM File会默认选择全球分辨率为900米的DEM数据，如果有更高分辨率的DEM数据，可以替换此数据，

点击Next；

3）在RPC Refinement面板中，有四个选项卡可以选择。
①切换到Advanced选项卡，修改输出像元大小Output Pixel Size为8米，其他参数保持默认；

② 切换到Export选项卡，选择输出文件格式，设置输出路径及文件名（GF1_PMS2_E104.0_N36_MSS_rpcortho.dat），点击Finish；

结果：

二、全色图像

1、打开

通过菜单栏File -> Open As-> China Satellites -> GF-1，在弹出的对话框，选择GF1_PMS2_E104.0_N36.0_20140724_L1A0000284766-PAN2.xml文件打开数据。

2、辐射定标

1)在Toolbox中，选择Radiometric Correction -> Radiometric Calibration，在弹出的File Selection对话框中，选择GF1_PMS2_E104.0_N36.0_20140724_L1A0000284766-PAN2.tiff，点击OK；

2)在Radiometric Calibration参数设置面板中，设置如下

3)设置辐射定标结果输出路径和文件名；

4)点击OK开始执行

5)右下角显示进度，定标结束后结果会自动加载到Data Manager中。

2、正射校正

1）点击Toolbox→Geometric Correction→Orthorectification→RPC Orthorectification Workflow，打开正射校正流程化工具；

2）在File Selection面板中，Input File选择经过大气校正的多光谱数据（GF1_PMS2_E104.0_N36_PAN_reflect.dat），DEM File会默认选择全球分辨率为900米的DEM数据，如果有更高分辨率的DEM数据，可以替换此数据，点击Next；

3）在RPC Refinement面板中，有四个选项卡可以选择。
①切换到Advanced选项卡，修改输出像元大小Output Pixel Size为2米

② 切换到Export选项卡，选择输出文件格式，设置输出路径及文件名（GF1_PMS2_E104.0_N36_PAN_rpcortho.dat），点击Finish；

三、融合