envi利用两景影像进行云替换

很多时候，我们下载了遥感影像，但研究区内会有少量的云。在能找到日期相近且无云重叠影像的情况下，我们便能对有云区域进行替换。原始影像b1，替换影像b2.

图中只有少量云覆盖，但对接下来实验会产生影响。

提取云的步骤参考：利用USGS自带的云掩膜图像去云 RGB 3S博客

在提取出云后，分别对两景影像进行辐射定标，大气校正，裁剪。

（1）选择File > Save As > Save As (ENVI, NITF, TIFF, DTED).
（2）在文件对话框中选择b1，点击Mask按钮。
（3）选择云检测结果作为掩膜文件。
（4）选择相应的输出的格式和路径，点击OK输出。

得到的图像中有云区域的像素值变成nodata。

我们需要将nodata值变为0值后，再用b2中相应位置的值替换0值。

 使用bandmath 

finite(b1, /nan)*(9999) or (~finite(b1, /nan))*b1 将修改NaN为9999值

(b1 eq 9999)*b2+(b1 ne 9999)*b1           用b2中相应位置的值替换9999值

最后拼接便完成了。