ENVI5.3下高分二号（GF2）数据预处理

以一景2015年1月23日获取的GF2-PMS1数据为例介绍在ENVI5.3下GF2数据预处理的详细操作步骤。GF2数据预处理基本流程如下：

图：GF2数据预处理流程

说明：1. 针对不同的应用，有不同的处理流程，上图中列出了两种常用的预处理流程。流程一主要针对高精度的定量遥感应用，也就是对大气校正精度要求比较高应用，比如：植被参数定量反演等；流程二主要针对定性遥感或者对大气校正精度要求比较低的遥感应用，比如：土地利用类型分类等。本文介绍的主要是流程二的详细操作步骤，流程一的实现可参考博文：ENVI5.2下高分二号数据FLAASH大气校正；另外，中国资源卫星应用中心网站已经公布了最新的GF2数据绝对辐射定标系数和两个传感器的波谱响应函数，大家可以下载使用。2. 本例中所有操作都是在ENVI5.3版本下进行的，除NNDiffuse Pan Sharpening图像融合（ENVI5.2新增，ENVI5.1中可以使用G-S融合方法）外，其他操作在ENVI5.1/5.2下同样可以完成。

1. 数据打开

启动ENVI5.3，在菜单栏中，选择File > Open，弹出Open对话框，找到GF2数据文件夹所在位置，选中扩展名为.tiff的两个文件，点击打开。

图2 打开GF2多光谱和全色数据

在左侧图层管理Layer Manager面板中，选择多光谱或全色数据图层，右键View Metadata查看其元数据信息，可以看到ENVI很好地识别了数据的RPC信息。

图3 ENVI自动识别GF2数据RPC信息

2. 正射校正

有了RPC信息之后，下面我们就可以基于这些RPC信息分别对多光谱和全色数据进行正射校正。这里我们以多光谱数据正射校正为例，全色数据正射校正操作完全相同。

在Toolbox中，选择Geometric Correction > Orthorectification > RPC Orthorectification Workflow，打开正射校正流程化工具。

在File Selection面板中，Input File选择多光谱数据，DEM File会默认选择全球分辨率为900米的DEM数据，我们这里保持默认（如果有更高分辨率的DEM数据，可以替换此数据），点击Next。

图4 正射校正流程化工具文件选择面板

在RPC Refinement面板中，有四个选项卡可以选择。如果我们有实测的或从其他途径（比如：从地形图中采集）获取的控制点数据，可以在该面板中进行添加，添加后在Statistics选项卡中可以看到相应的误差统计信息（关于这两个选项卡的详细说明和操作可以参考博文：【ENVI入门系列】05.图像正射校正）。

切换到Advanced选项卡，修改输出像元大小Output Pixel Size为4米，重采样方法Image Resampling选择三次卷积法，其他参数保持默认。

图5 正射校正流程化工具参数设置选项卡

切换到Export选项卡，选择输出文件格式，设置输出路径及文件名，点击Finish。

至此，我们已经完成了GF2多光谱数据的正射校正，全色数据的正射校正操作与多光谱数据的正射校正完全相同，需要提醒的地方是GF2全色数据正射校正时输出像元大小需设置为1米，以便我们下面进行图像融合。

3. 图像融合

在进行图像融合之前，我们需要对正射校正后的多光谱和全色数据进行浏览，查看二者是否完全配准，如果没有完全配准，就需要对其进行配准，这里我们推荐大家使用ENVI中的自动配准流程化工具，以全色数据为基准对多光谱数据进行配准，此工具的位置在：Geometric Correction > Registration > Image Registration Workflow，关于该工具的详细说明和操作可以参考博文：【ENVI入门系列】07.图像自动配准。

图6 查看GF2多光谱和全色正射校正结果配准情况

对本次操作中正射校正后的多光谱和全色数据的对比发现，二者配准的比较好（目前，大部分高分辨率数据正射校正后多光谱和全色数据配准的均比较好），所以我们这里不进行图像配准，直接进行图像融合。图像融合方法我们使用从ENVI5.2版本起新增的NNDiffuse Pan Sharpening方法（关于该工具的介绍可以参考博文：ENVI5.2中的NNDiffuse融合方法），大家也可以使用Gram-Schmidt Pan Sharpening等其它方法进行融合。

在进行融合之前，还需要说明的是：经过我们的测试，不同的数据存储格式会对融合的速度产生影响，当多光谱数据的存储格式为BIL或BIP时，其融合速度较BSQ格式来说，可以提高三倍左右。所以，我们建议大家融合之前可以先查看其数据存储格式，如果是BSQ，可以将其转化为BIL或BIP（转换时需要用到的工具：Raster Management > Convert Interleave），以便提高融合效率，详细说明可以参考博文：一种提高图像融合及其他一些处理速度的方法。这里，我们也为大家提供了将上述步骤与融合步骤合并为一个步骤的扩展工具：NNDiffuse Pan Sharpening (BIL)。

在Toolbox中，选择Extensions > NNDiffuse Pan Sharpening (BIL)，弹出NNDiffuse Pan Sharpening面板。Input Low Resolution Raster选择上一步正射校正后的多光谱数据，Input High Resolution Raster选择上一步正射校正后的全色数据，其他参数保持默认。这里，需要说明的是NNDiffuse Pan Sharpening工具要求输入的多光谱和全色数据的空间分辨率是整数倍的（比如：本例中正射校正时分别将多光谱的全色的分辨率重采样为4米和1米，就是为了方便该工具的使用）。如果二者不是整数倍，使用此工具时需要将其重采样成整数倍；该工具对输入文件还有一些要求，平时我们使用的数据也基本满足这些要求，如果使用过程中报错可以查看帮助文档，判断输入数据是否满足相应要求。

图7 NNDiffuse Pan Sharpening参数设置面板

至此，我们已经完成了GF2数据图像融合处理，融合结果局部截图如下所示：

图8 GF2数据融合前后对比（左：融合前，右：融合后）

4. 快速大气校正（可选）

对于部分定性和定量遥感应用来说，大气校正是一个必需的预处理步骤。而对于一些对大气校正精度要求不是很高的应用来说，大气校正我们可以采用ENVI提供的快速大气校正工具进行，该工具在ENVI5.3版本中有所改进，具体改进可以查看博文：ENVI5.3新特性介绍中图像处理和界面改进部分。对于一些非定量遥感应用来说（比如：制作影像底图），该步骤可以不做。

在Toolbox中，选择Radiometric Correction > Atmospheric Correction Module > QUick Atmospheric Correction (QUAC)，弹出QUAC参数设置面板。Input Raster选择上一步图像融合后的结果，Sensor Type选择Generic/Unknown Sensor，设置输出路径和文件名，点击OK。

图9 快速大气校正参数设置面板

快速大气校正前后同一位置的植被光谱曲线如下图所示。可以看到，快速大气校正基本消除或减弱了大气的影响。

图10 快速大气校正前后植被光谱曲线对比

http://blog.sina.com.cn/s/blog_764b1e9d0102wrdb.html#cmt_564FDAEE-7F000001-6D9E5B31-8E9-8A0

ENVI5.3下高分二号（GF2）数据预处理相关推荐

Envi5.3——高分二号PMS数据预处理
首先第一步:打开数据open--以TIFF格式打开开始进行正射校正--RPC正射工作流程,我们先对多光谱进行校正输入需正射校正文件和DEM文件,在这里我们不用自带的2010年的DEM文件. 点击下 ...
第059篇：高分二号遥感影像预处理流程（ENVI5.3.1平台+ENVI App Store中最新的中国国产卫星支持工具）
今天被袁老的新闻刷屏,湖南衡水县水稻基地传出好消息: 袁隆平团队第三代杂交水稻测产,测得晚稻平均亩产为911.7 ...
高分二号(GF-2)号卫星数据的查询下载地址和方法
高分二号卫星是我国自主研制的首颗空间分辨优于1米的民用光学遥感卫星可在平台中查询到,搭载有两台高分辨率1米全色.4米多光谱相机,具有亚米级空间分辨率.高定位精度和快速姿态机动能力等特点,有效地提升了卫 ...
IDL下高分二号完整预处理代码
下面写的太LOW了,有好多值得修改的地方,修改更新后博客地址: https://blog.csdn.net/desertsTsung/article/details/84679969 今天去不了Ear ...
使用ERDAS对国产卫星影像进行控制点正射校正—以高分二号卫星（GF2）为例
** 使用ERDAS对国产卫星影像进行控制点正射校正-以高分二号卫星(GF2)为例 ** 一.研究区影像: 1.1.高分二号原始影像影像为山西介休市的一景获取时间为2019.01.22的GF2影像: ...
基于IDL的高分二号影像批量预处理程序
基于IDL的高分二号影像批量预处理程序(第二版) 上一版本插件预处理效果不太理想(特别是融合后分辨率不符合预期),上一版资源不再提供下载,但处理思路可以借鉴.目前在本机40GB RAM工作站上预处理一 ...
基于python实现高分二号遥感影像水体提取与水质反演（黑臭水体与水体富营养化）
高分二号遥感影像水体提取与水质反演水体提取函数--NDWI 基于几何约束提取河流生成shp,方便后续裁剪水体水质反演最终结果水体提取函数--NDWI 水体提取函数water.py impor ...
遥感技术及高分遥感影像在地震中的应用及高分二号获取
长期以来,地震预报监测.灾害调查.灾情信息获取主要依靠实地勘测手段,其获取的数据精度和置信度虽然较高,但存在工作量大.效率低.费用高和信息不直观等缺点.遥感技术手段可在一定程度上克服传统实地勘测手段的 ...
遥感图像语义分割——从原始图像开始制作自己的数据集(以高分二号为例)
遥感图像语义分割--从原始图像开始制作自己的数据集(以高分二号为例) 文章目录遥感图像语义分割--从原始图像开始制作自己的数据集(以高分二号为例) 1.遥感影像获取 2.遥感数据预处理(影像融合) ...

ENVI5.3下高分二号（GF2）数据预处理

ENVI5.3下高分二号（GF2）数据预处理相关推荐

最新文章

热门文章