ENVI5.6.2中北京三号N数据处理方法

标签：选择校正 NNDiffuse 三号 ENVI5.6 控制点图像数据处理

北京三号N卫星即北京三号国际合作星，由四颗0.3米分辨率卫星组成，首批两颗卫星分别于2021年4月29日和8月17日成功发射，除了常见的蓝、绿、红、近红外波段，还提供了深蓝（400-450nm）和红边（700-750nm）两个波段。成像宽幅为14km，还具备单轨立体和三立体成像能力。

ENVI5.6.2中完全支持北京三号N卫星4波段或者6波段图像数据处理，包括数据打开、辐射定标、大气校正、正射校正、图像融合等处理。本文介绍利用北京三号N数据生产DOM产品，包括正射校正和图像融合处理。

ENVI5.6.2软件的试用请浏览：envi.geoscene.cn/envi_license

1、数据打开

请使用易智瑞公司遥感事业部开发的“中国国产卫星支持工具”打开北京三号N卫星。

建议使用ENVI App Store进行ENVI扩展工具的安装与管理。

ENVI App Store下载地址：envi.geoscene.cn/appstore

在App Store找到中国国产卫星支持工具，点击安装进行安装。
重启 ENVI即可。

或手动安装扩展工具： envi.geoscene.cn/china

启动ENVI。选择主菜单 File > Open As > China Satellites >BJ3N。选择XML文件打开数据。
在Layer Manager中，右键选择菜单View Metadata，选择各个选项可以看到自动读取的包括RPC、定标参数、拍摄参数、成像时间等元数据信息。

说明：北京三号N卫星主要包括4波段和6波段数据产品，打开和后续处理一样。

2、正射校正

正射校正可以选择无控制点/手动选择控制点正射校正和参考影像自动选择控制点正射校正两种工具。

2.1无控制点/手动选择控制点正射校正

在Toolbox中选择/Geometric Correction/Orthorectification/ RPC Orthorectification Workflow工具，该工具是流程化操作方式操作比较简单，不选择控制点将以无控制点正射校正输出。

在RPC Orthorectification Workflow面板中，输出像元大小（Output Pixel Size）默认自动计算整个图像的平均值，可以手动修改为全色0.3米、多光谱1.2米。

注：NNDiffuse Pan Sharpening图像融合算法要求全色图像分辨率是多光谱图像的整数倍。

如下图为无控制点的正射校正结果叠加在天地图影像底图上，天地图影像放大最高级切片肉眼区分不出几何偏差。

图：无控制点正射校正精度对比

2.2参考影像自动选择控制点正射校正

参考影像应选择与待正射校正接近分辨率的影像，如0.3~2米。

在Toolbox中选择/Geometric Correction/Orthorectification/RPC Orthorectification Using Reference Image，打开基于基准影像的RPC正射校正工具；在面板上进行数据选择和参数设置：

Input Raster：选择前面打开的北京三号N全色或者多光谱图像；
Input Reference Raster：具有准确地理位置信息的参考影像；
Input DEM Raster：覆盖待校正影像区域的DEM数据；
DEM Is Height Above Ellipsoid：DEM数据是否是椭球高；
Requested Number of GCPs：需要采集的控制点最大数量；
Search Windows Size：搜索窗口大小，窗口越大找点精度越高，时间越长；
Output Coordinate System：可选项，如果不设置，则输出的正射校正结果是UTM WGS84坐标系的；
Output Pixel Size：输出像元大小，北京三号N全色为3米，多光谱为1.2米。

注：输出像元大小为可选项，如果不设置，则输出的像元大小为根据RPC自动计算的平均像元大小，建议手动设置输出数据的分辨率。

Image Resampling：重采样方法，默认双线性内插Bilinear；
Grid Spacing：像素栅格间距。保持默认10；这个值越大，校正速度越快，精度越低。
Output Raster：正射校正结果输出路径及文件名；
Output GCPs：控制点文件输出路径及文件名；

参数设置好之后，点击OK开始执行。

3、图像融合

推荐使用NNDiffuse 图像融合方法，该方法能很好的保留纹理和光谱信息。由于NNDiffuse方法对储存顺序为BIP或者BIL的多光谱图像处理效率更高，处理效率大约为BSQ储存顺序2~3倍。因此我们可以利用/Raster Management/Convert Interleave工具对储存顺序为BSQ的多光谱图像转为BIP或者BIL。

在Toolbox中选择/Image Sharpening/NNDiffuse Pan Sharpening，在NNDiffuse Pan Sharpening参数面板中分别选择全色和多光谱正射校正结果，其他参数选择默认。

Pixel Size Ratio：多光谱与全色图像的分辨率比值，不设置情况自动计算。NNDiffuse算法要求是整数倍。
Spatial Smoothness：设置空间重采样插值核参数，不填写默认为分辨率比值62。
Intensity Smoothness：平滑强度值，值越小得到的结果图像越清晰（类似高通滤波效果），可用于图像可视化。值越大得到的图像噪声越少，可用于图像分类或者定量反演。不填写默认为自动根据图像计算，可填写14.1~20之间的值。