实验九遥感图像预处理

实验九遥感图像预处理（三）

一、实验内容

图像融合、镶嵌（2学时）

图像裁剪（2学时）

二、实验学时4学时

二、实验原理、方法和手段

图像融合、镶嵌、裁剪原理内容在操作中进行介绍.

三、实验数据

实验数据：第四章:遥感图像预处理

四、实验步骤

1.图像融合

数据：第四章：遥感图像预处理\5-图像融合。

图像融合是将低空间分辨率的多光谱图像或高光谱数据与高空间分辨率的单波段重新采样,

生成一幅高分辨率多光谱遥感图像的图像处理技术。使得处理后的图像既有较高的空间分辨率，

又具有多光谱特征。

ENVI中提供了两种融合方法:HSV变换和Brovey变换。这两种方法均要求数据具有地理参考

或者具有相同的尺寸，RGB输入波段必须为无符号8-bit数据或从打开的彩色display中选择。

两种方法基本类似，下面介绍Brovey变换操作过程。

操作过程：

1o打开融合的两个文件：TM-30m。img和bldr_spoimg（分别在两个display窗口中显示），将

TM—30m.img以RGB格式显示在dispIay窗口中。

2o选挂主菜单TtransformTimagesharpeningTcoIornormaIized(brovey),在selectinput

RGB对话框中，有两种选择方式(如第一图)：可用波段列表中选择或display窗口中选择，选

择display#1窗口中的RGB,电击OK按钮。

3.选中相应波段，双击，进入coIo.normaIized(brovey)对话框(如第二图)，在

coIo.normaIized(brovey)对话框中，选择重采样方式(resampIing)和输入文件路径及文件名，

单击OK按钮输出结果。融合后结果如下，可以对两幅图像链接进行比较.

对于多光谱图像，ENVI可以利用以下融合技术：

Gram-Schmidt:能保持融合前后图像波谱信息的一致性。

Colornormalized:要求数据具有中心波长和FEHM

主成分(PC)变换

下面介绍参数相对较多的Gram-schmidt操作过程。

1o打开融合的两个文件：qb_bouId_msioimg和qb_bould_pan。imgo

2O选择主菜单TtransformTimagesharpening->gram-schmitspectralsharpening,

在selectlowspatiaIresolutionmultibandinputfiIe对话框中选择低分辨率多光谱

图像qb_bouId_msi.img,双击选中图像波段，弹出gram—schmitspectralsharpening对话框。

SelectHighSpatialResolutionPanInputBand区）

：

S«l«ctInputBond：BandInformation

-sqb_boMlder_panimg

oBand1

也M«ipTit£u

Q包11qb_boulder_msi.img

oB5d1(4850000)

aBand2©60.0000)

oBand3(660.0000)

oBwd4(8300000)

(±■MapInfo

0回1brovwy^换

OCNSharpR(WarpCBand5：bldr

OCNSharpGOffarp(Btnd4:bldr

[0K][Cancel][Previous1Open▼

3o在gram—schmitspectraIsharpening对话框中，选择降低高分辨率全色波段的方法。

4.选择重采样方式(resampling)和输入文件路径及文件名，单击0K按钮输出结果。融合后结果如下，可以对

两幅图像链接进行比较。

2。图像镶嵌

将没有地理坐标系或者有地理坐标系的多幅图像合并，生成一幅单一的合成图像.主要是镶

嵌颜色不一致、接边以及重叠等问题。

⑴有地理参考的图像镶嵌

第一步：启动图像镶嵌工具

在ENVI主菜单中，选择mapTmosaickingTGeoreferrnced,打开mapbasedmosaic对话框.

上述对话框由菜单条、窗口图像和文件列表和快捷菜单组成.

第二步：加载镶嵌图像

1o在mosaic对话框中，选择importTimportfiles,选择mosaicl和mosaic2镶嵌文件导

入。

2.导入的镶嵌文件显示在图像窗口以及文件列表，文件列表中的排在下面的文件，在图像

显示窗口中显示在上层。

3.在文件列表中选择需要调整顺序的文件，单击右键选择快捷菜单raise(lower)imageto

top(bottom)(提升/降低到顶层/底层)或raise(lower)imageoneposition(提高/降低一、层),

或者在图像窗口中单击右键选择快捷菜单。通过这个功能调整图像叠加顺序。

第三步：图像重叠设置

1o选择文件列表中一个文件，单击右键选择editentry.

2O在editentry对话框中，设置datavaIuetoignore:0,忽略0值；设置feathering

distance为10,羽化半径为10个像素，单击0K按钮。

第四步：切割线设置

1o在mosaic对话框中，选择fileTsavetempIate,选择输出路径和文件名，将模板文件

显示在display中。

2.在display中，选择overlay—annotation,在重叠区域绘制一条折线当做切割线;绘制一

个symbol放在切割线一旁，标示这部分将被裁剪，注意单击两次右击以完成symbol注记的绘

制；保存注记文件。

3o回到mosaic对话框中，在文件列表最下面文件单击右键，选择editentry,在entry

参数对话框中，单击selectcutIineannotationfile按钮，选择前面生成的注记文件，单击

clear按钮可以清除注记文件.

第五步：颜色平衡设置

1.在mosaic对话框中，首先确定一个图像当做基准，在文件列表中选择这个图像，单击右

键选择editentry,打开entry对话框。

2o将mosaicdispIay设置为RGB,并选择波段合成RGB图像显示；选择colorbaIancing

参数：fixed,作为基准图像；同样的方法对其他图像文件进行设置，选择colorbalancing参

数：adjust.

第六步：输出结果

1o在mosaic对话框中,选择fi在—apply,在mosaicparameters对话框中，设置输出像

元大小、重采样方式、文件路径及文件名、背景值.其中colorbalanceusing选项中，默认的

是统计重叠区的直方图，可以单击□按钮切换到统计整个基准图像的直方图用于颜色平衡。

整个镶嵌过程已经完成，显示镶嵌结果查看效果.在上述步骤中，其中第三-五步都是可选项,

根据实际情况选择.

(2)基于像素的图像镶嵌

第一步：启动图像镶嵌工具

在ENVI主菜单中，选择mapTmosaickingTpixelbased,开始进行ENVI基于像素的镶嵌操

作，pixelbasedmosaic对话框出现在屏幕上。

第二步：加载镶嵌图像

1o在mosaic对话框中，选择importTimportfiles,选择相应的镶嵌文件导入。

2.在selectmosaicsize对话框中，指定镶嵌图像的大小，这个可以通过将所有的镶嵌图

像的行列数相加，得到一个大概的范围，设置“Xsize”为1028“Ysize”为1024。

第三步：调整图像住置

1o在mosaic对话框的下方X0文本框和Y0文本框输入像素值，调整图像位置，也可以在

图像窗口中，点击并按住鼠标左键，拖拽所选图像到所需的位置，然后松开鼠标左键就可以放

置该图像了。

2O如果镶嵌区域大小不合适，选择optionsTchangemosaicsize,重新设置镶嵌区域大小.

其他步骤和有地理参考的图像镶嵌类似，这里不一一列举.

3.图像裁剪

图像裁剪的目的是将研究之外的区域去除。

常用的方法是按照行政区划边界或自然区划边界进行图像裁剪。

在基础数据生产中，还经常要进行标准分幅裁剪，按照ENVI的图像裁剪过程,可分为规则裁

剪和不规则裁剪。

(1)规则分幅裁剪

规则分幅裁剪是指裁剪图像的边界范围是一个矩形，这个矩形范围获取途径包括行列号、左

上角和右下角两点坐标、图像文件、R0I/矢量文件。

操作步骤：

1.在主菜单中，选择fileTopenimagefiIe,打开裁剪图像bhtmsat.img。

2.在主菜单中，选择file—savefiIeasTENVIstandard,弹出newfiIebuiIder对话框.

在该对话框中，单击importfile按钮，弹出creatnewfileinputfile对话框.

3.在creat.ne.fiI.inpu.fiIe对话框中，选中selec.inpu.file列表中的裁剪图像图像,

单击spatia.subset按钮（空间波段子集），在spatia.subset对话框中，单击image按钮，弹

出subse.b.image对话框，在毋选波段中进行子波段裁剪范围设置。

好SelectSpatialSubset区）

4.在subse.b.image对话框中，可以通过输入行列数确定裁剪尺寸，按住鼠标左键拖动图像

中的红色矩形框确定裁剪区域，或者直接用鼠标左键按钮红色边框拖动来裁剪尺寸以及位置，

单击0K按钮。

5.在selec.spatia.subset对话框中可以看到裁剪区域信息，单击0K按钮.....□

6O在creatnewfiIeinputfiIe对话框中，可以通过spectralsubset按钮选择输出波

段子集，单击0K按钮。

»CreateHewFileInputFile区J

7o选择输出路径及文件名，单击OK按钮，完成规则分幅裁剪过程.

⑵不规则分幅裁剪

A.不规则分幅裁剪是指裁剪对象的外边界范围是一个任意多边形。任意多边形可以是事先

生成的一个完整的闭合多边形区域，可以是一个手工绘制的R0I多边形，也可以是ENVI

支持的矢量文件.针对不同的情况采用不同的裁剪过程，下面介绍两种方法：

B.手动绘制感兴趣区

1.打开图像can_tmr。img并显示在dispIay窗口中。

2.在image窗口中选择overIay->regionofinterest.在ROITooI对话框中，单击ROI_Type

TPolyon。

绘制窗口（window）选择image,绘制一个多边形，右键结束.根据需求可以绘制若干个多边

形。

选择主菜单一basictools->subsetdataviaROIS,或者选择ROITool->fiIe->subsetdata

viaROIS,选择裁剪图像，双击左键，进入spatialsubsetdataviaROIs对话框。

SelectByFile回

[0KCancel]Previous|Open▼

在spatialsubsetdataviaROIs对话框中,设置以下参数:

在ROIs列表中(seIectinputROIs）,选择绘制的ROIS.

在amaskpixelsoutsideofROI”项中选择：Yeso

裁剪背景值(maskbackgroundvalue):0

选择输出路径及文件名，单击OK按钮，裁剪图像.（以下为其中一个窗口的裁剪结果）

C.矢量数据生成感兴趣区

在主菜单中，选择fiIeTopenvectorfiIe,打开裁剪图像所在区域的shapefiIe矢量文件,

投影参数不变，选择导入的memory。

在avaitablevectorIist对话框中，选择fiIeTexportIayertoROI,在弹出的对话框中

选择裁剪图像，单击0K按钮。

❷AvailableVectorsList[»}[D||3^|

在exportEVFlayertoROI选择对话框中，选择将所有矢量要素转成一个ROI(convertaII

recordofanEVFIayertooneROI)，单击OK按钮°

1.选择主菜单一basictoolsubsetdataviaROIs,选择裁剪图像。

在spatialsubsetviaROIsparameters中，设置以下参数:

在ROIs列表中(seIectinputROIs),