一种光谱无损的遥感数据融合方法——LiIFT插值融合（二）

接上篇 《一种光谱无损的遥感数据融合方法——LiIFT插值融合（一）》

三、LiIFT融合RSD平台操作

① PANSHARP LiIFT融合的PAN/MSS自动对齐

 融合前先确认下PAN和MSS是不是配准了的。有一些软件的方案是将PAN和MSS先比较一下。但是着这有一个问题，将PAN像元分辨率降到和MSS一样，如果目视或者自动探测有一个像元的误差，融合后可能就偏离了4个像元，质量还是不太好。

RSD提供了一种PAM/MSS自动配准的功能。LiIFT融合或者其它融合时会自动监测并配准这些多光谱和全色数据集。

图1  开启PAN/MSS自动配准融合结果比较

图1是一个GF1 数据集融合结果比较。其中层3是未启用PAN/MSS自动配准的融合（图像在底层），层4是启用PAN/MSS自动配准后的融合结果。任务主窗口是利用显示图像的“蛙眼”交替显示前后两层进行比较。较亮的框内是底层，和外框图像一致的是顶层。交替闪烁的结果初看好像没有对齐，图像上下移动。但是仔细观察（比如将光标固定放在线状地物的交叉点），可见线状地物并不移动，移动的仅仅是面状的颜色。说明自动配准对面状的MSS数据相对PAN进行了移动调整，改善了融合效果。

每次融合都启动PAN/MSS自动配准通常都能保证更好的融合效果，但是这个配准也需要比较长的时间，所以RSD为此提供了一个选项，让用户选择是否进行自动配准。见图2设置窗口的选项。

图2  PANSHARP融合时PAN/MSS自动配准选项

手动选择图2中红框的值就可以启动自动配准，绝大多数情况下都可以解决PAN/MSS不齐的问题。但是数据集质量较好大多不需要PAN/MSS重新对齐时会耽误较长的时间。开始不选择PAN/MSS自动配准，对融合后检查不合格的，再次进行PAN/MSS自动对齐融合也是一种较好的解决办法。

② 菜单命令进行 GF1、2、6 全色锐化（PANSHARP）LiIFT融合

图3 GF2 PANSHARP 融合的菜单项

点击图3的菜单命令，出现一个选择文件的对话框，在里面分别选好PAN和MSS数据，点击确定。就可以等待结果了。融合后的结果是任务中一个单独的层。是否进行了PAN/MSS自动配准取决于图2的选项。

GF1 、GF6都有同样的菜单命令可以直接融合。GF1 B、C、D星数据和A星数据一样，不需要特别指定。

③ 批量的融合

批量融合与批量处理是一样的，需要预先准备一个OpenList的打开文件列表。图4。

图4 建立打开文件列表OpemList

通过图4菜单，交互选取GF1、2、6的TIFF文件（选择PAN或者MSS数据都可以，但是只选1个）。点击保存，就会保存到RSD工作路径里面一个\RsdData\ RSD_FileNameList.txt 文件。打开这个文本文件，内容大概为图5.。

图5 打开文件列表

这个文件也可以手动编辑生成。

下一步通过图6开始批量处理，我们选择处理GF2数据（文件列表中选择的是GF2）。

图6 开始批量处理GF2数据的菜单命令

点击菜单后出现图7的对话框

图7 选择处理方式

由于我们是想进行GF2数据的PANSHARP融合，所以我们选择上述图7中多光谱、全色数据正射校正+融合。如果你希望同时进行大气校正，也可以选择 多光谱数据正射、大气校正+全色数据正射校正+融合选项。

点击确定，会出现打开文件列表进行确认，点击确定，就可以等待处理结果了。结果存放在\RsdData工作目录中。

④ 菜单命令进行通用的PANSHARP  LiIFT 融合

如果已经在一个任务中加载了高分辨率的基准（依据）数据和待融合的数据，可以使用图8的菜单命令进行融合。

如果你不需要进行PAN/MSS配准，可以直接进行PANSHARP融合，建议融合步骤如下(用一个GF2的例子来演示)。

第一步：打开 MSS 数据，如：GF2_PMS2_E122.0_N41.1_20200517_L1A0004806215。正常加载后保存为一个.rsd文件：GF2_PMS2_E122.0_N41.1_20200517_L1A0004806215-MSS2.rsd

第二步：打开上述数据对应的PAN数据。这时是一个新建的任务。

第三步：从菜单 任务管理->添加RSD文件，将第一步的文件添加到当前任务中。结果见图8。注意第二层是一个ALIEN GRID 并且在主窗口可能看起来就是一些稀疏的小点子。

图8 加载了PAN数据和MSS数据

第四步：从菜单 GRID 管理->光谱数据融合->全色锐化（PANSHARP）融合，点击后出现对话框图9.。

 

图9 PANSHARP 融合参数选择对话框

图中来源数据按照红框和蓝框选择，分辨率差异倍数是4，融合处理窗口大小比分辨率差异倍数稍大一点就可以。注意这个数据在可处理范围融合程序接收这个数据，否则融合程序调节这个数据。

这个对话框中不可以指定融合前进行自动配准。点击确定就可以得到融合结果了。图13上。

如果你按照上面步骤进行了没有自动进行PAN/MSS配准的融合，融合后发现没有对齐结果不合格（图13上）。这是你可能需要自动配准的融合。步骤如下(仍然用上述GF2的例子来演示)：

第一步：打开GF2 的 PAN数据（等同上面的第二步）。或者直接使用上面例子使用的任务。

第二步： 执行命令菜单（图10） 为任务添加光谱数据。这里是添加未正射的MSS数据。(不是上面例子的.rsd文件)

图10 添加多光谱数据

注意添加以后任务增加了一个新层，这个新层不是稀疏的点子，是进行了线性插值的0.8m的MSS数据。图11.。

图11 3.2m GF２　MSS数据插值为０.８ｍ的图像。

注意：RSD并不要求融合图像有相同的空间分辨率，但是相同分辨率同样能融合。这里插值成相同分辨率的目的主要是为了和PAN自动配准。

第三步：同上述第四步，从菜单 GRID 管理->光谱数据融合->全色锐化（PANSHARP）融合，点击后出现对话框图12。图12和上述图9是一样的，仅选项不同。

图12 PANSHARP 融合自动配准参数选择对话框

注意选择 融合前自动配准 ，点击确定，就得到融合后结果，在任务中可见一个新层。图13下。

图13  自动配准前后的融合结果比较

融合前的自动配准不是专属LiIFT融合的，任何其它融合都可以使用这个功能。该功能使融合更精确，重要的是更省心，无需用户操心就能自动得到精确理想的融合结果。缺点是稍稍费了一点时间。

⑤ 菜单命令进行两个普通数据集的  LiIFT 融合

这里是一个哨兵2号10m数据和Lansat8 OLI 30m数据融合的例子。

第一步：任务中加载两个数据集，可以分辨率一致（NATIVE GRID），也可以不一致（ALIEN GRID）。

第二步：从菜单 GRID 管理->光谱数据融合->GRID层融合，点击后出现对话框图14，与图9是一样的对话框，只是选项不同。

图14 不同数据集融合参数选择对话框

图15 Landsat8 OLI 30m数据插值融合到哨兵2号的10m分辨率

 详情加企鹅群136965427，在这里解答和讨论有关遥感数据处理和RSD平台的有关技术问题。

加企鹅758461012，原来的满了。