实习二 遥感图像的校正、镶嵌

《实习二 遥感图像的校正、镶嵌》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实习二熟悉软件，遥感图像的校正、镶嵌[实验目的]1.学会利用ENVI软件对遥感图像进行校正。掌握遥感图像几何校正的过程，了解其中控制点和方法的选择。2学会对图像进行镶嵌(mosaic)和切割(subsize)[实验内容]图像几何校正Image-to-ImageRegistration（图像-图像的配准）（1）用AvailableBandsList打开基图像和纠正图像文件，并在两个窗口显示它们。（2）一旦两幅图像都已经显示，选择Map>Registration>SelectGCPS>ImagetoImage。（3）出现I

2、magetoImageRegistration对话框时，在“BaseImage:”下面点击需要显示的名字，选择基准图像（参照图像）。（采用njtmcorrected,从头文件中查看其信息）（4）在“WarpImage:”下方点击需要显示的名字，选择被纠正的图像。(采用spot5-sx)（5）点击“OK”，出现GroundControlPointsSelection对话框。选择地面控制点：在基础和纠正图像中，选择GCP的位置。①为每幅图像移动缩放窗口到需要的GCP区域。②在缩放窗口的一个特定像元上点击鼠标左键，把光标定位

3、在该像元或像元的一部分上。在GroundControlPointsSelection对话框，被选择处的坐标分别显示在标签为“BaseX,Y”和“WarpX,Y”文本区中。③一旦两幅图像都选择了需要的像元，在GroundControlPointsSelection对话框中点击“AddPoint”，将选择的GCPs添加到已经选择的X、Y坐标对列表里。他们将按基图像、纠正图像顺序被列出。当已经选择了四个或更多个GCPs，对选择的纠正预测（Predict）的X、Y坐标将显示在随后的列中，之后显示的是X、Y的误差列表，RMS误差

4、列表显示在最后列。当GCPs已经被添加到列表中时，在两幅图像的缩放窗口里选择的像元处画一个标记。选择的GCP数将出现在标记附近。标记中心（在十字准线下面）显示了实际的GCP位置。④用同样的方法添加其它的GCPs。使用ShowList按钮可以显示所有控制点的信息。点击deletelastpoint可以删除最后添加的控制点。在GroundControlPoints列表中选择要编辑的GCP，双击，可修改当前值，以编辑基图像和纠正4图像的X和Y值。⑤控制点选择完成后，选择File->SaveGCPstoASCII保存控制点文件

5、。（6）像元重采样一旦为图像-图像或图像-地图配准选择好GCPs，用下列步骤指定纠正和重采样方法：①在GroundControlPointsSelection对话框，选择Options>WarpDisplayedBandorWarpFile。采用WarpFile②出现RegistrationParameters对话框时，从与“WarpMethod”邻近的下拉菜单中，选择需要的纠正方法。可以使用RST(旋转、缩放和平移)多项式（polynomial）和三角测量(triangulation)纠正方法。RST纠正是最简单的方

6、法，需要三个或更多的GCPs运行图像的旋转、缩放和平移。可以得到1次到n次多项式纠正。可以得到的次数依赖于选择的GCPs数，德洛内三角测量适于三角到不规则空间GCPs和内插数值到输出格网中。采用RST③从“Resampling”标签附近的下拉菜单里，选择需要的重采样方法。可利用的重采样方法包括最近邻（nearestneighbor）、双线性(bilinear)和立方体卷积(cubicconvolution)。最近邻重采样运用没有解译的像元建立纠正图像。双线性重采样用四个像元进行线性内插，以对纠正图像进行重采样。立方体卷

7、积通过三次多项式对图像进行重采样，用16个像元近似正弦函数。值得注意的是立方体卷积重采样比其它方法速度慢。采用nearestneighbor（4）在“BackgroundValue”文本框里，输入DN值，设定背景值（在纠正图像里，DN值用于填充没有图像数据显示的区域）。（5）选择输出到“File”或“Memory”.若选择输出到“File”,输入一个输出文件名，用“Choose”按钮选择一个文件名。进行保存切割(resize)和镶嵌(mosaic)一图像的切割（取子区）ENVI：BasicTools>>ResizeDa

8、ta>>ResizeDataInputFile对话框（如下图）。①选择需要切割的原始图像；(用软件自带数据can-tmr.img)②选择SpatialSubset或SpectralSubset方式；③若设置空间切割方式（SpatialSubset>>selectSpatialSubset）点击“Image”；④出现Subsetby