融合影像处理.docVIP

将TM三个波段数据进行插值放大，使其分辨率与高分辨率图像SPOT保持一致。 1）地图配准 A、打开tmAtlanta.img和panAtlanta.img,用panAtlanta.img去配准tmAtlanta.img，打开tmAtlanta.img时，切换到Raster Options，选择Display as下拉条，选择Gray Scale，在Layer中输入2，选中Fit to Frame复选框。 B、打开Viewer菜单栏上的的Raster选项，选择Geometric Correction（几何校正）选项；在弹出的对话框中，选择Polynomial（几何校正属性），将多项次方改为2，点击Apply。打开几何校正工具窗口。 C、在地图上配准选点，选取明显的地物在两幅图上作出标记，为了便于显示，可以将点的颜色设置为鲜明颜色 D、配准好后，可以选择Point Matching设置，选择重采样，进行保存设置，如下图所示： E、将newtm.img和panAtlanta.img在一个viewer中打开，在菜单中点击AOI－Tools，选择一个实线矩形框工具，在图像上画出两幅图像共有的区域，然后在File－save－AOI Layer as将所做的AOI区域保存。 F、点击Interpreter—Utilities，在打开的面板上选择Subset..，然后将newtm.img和panAtlanta.img都用刚才的aoi文件裁成一个子图像。输出文件分别为subtm.img和subpan.img。 G、打开Viewer窗口并打开subpan.img，在Viewer菜单Utility－Layer info中查看subpan.img分辨率。 H、在Viewer中打开图像subTm.img，点击Raster－Geometric Correction。输入重采样后的文件名为subtm_res.img。重采样后的subtm_res.img与subPan.img分辨率相同，但是像元数目不同，还需要裁减成行列数相同的图像。用subtm_res.img和subPan.img重复步骤E~H，输出文件名分别为：sub_subpan.img， sub_subtm_res.img。 2、分别将两幅图像三个波段的数据采用指定的变换公式从RGB空间变换到IHS空间； 1）、在打开的RGB to HIS对话框中，选择输入输出文件名和波段。输入文件为sub_subtm_res.img,输出文件名为RGBtoIHS.img,波段选择R:3，G:2，B:1， 如下图所示： 2）、打开直方图匹配的对话框，进行直方图匹配， 并在对话框设置如下选项： Input File： sub_subpan.img Input File to Match： RGBtoIHS.img Output File：pan_histogrammatch.img Data Type Output：float single 3）、将直方图匹配后的pan_histogrammatch.img图像（1band）替换RGBtoIHS.img的I分量（第一波段）。选择Unlimites到LayerStack打开对话框，进行波段替换的选择，先在Input File中选择pan_histogrammatch.img ，然后点击下面的add，然后在Input File中选择RGBtoIHS.img，在Layer中选择2，然后点击add按钮，最后在Layer中选择3，点击add按钮。输出文件名为stack.img； 4）、对图像进行IHS转换为RGB，输入文件名stack.img,输出为stack_ihstorgb.img。 5）、直方图匹配，为了使融合效果更好，将stack.img与sub_subtm_res.img进行直方图匹配： Input File： stack.img Input File to Match： sub_subtm_res.img Output File：IHSend.img 3、结果评价 下图是处理过后的图像的效果显示：