遥感图像处理结课作业.docx

遥感图像处理结课作业.docx

《遥感图像处理结课作业.docx》由会员分享，可在线阅读，更多相关《遥感图像处理结课作业.docx（22页珍藏版）》请在冰豆网上搜索。

遥感图像处理结课作业

遥感图像处理实验报告

班级姓名学号

成绩评定教师签字

专题一：

DEM图像进行彩色制图

（叙述制图过程并把自己处理结果加载到本文档里）

一、处理过程：

1、ENVI彩色表的应用

ENVIColorTables选项允许对灰阶图像进行线性对比度拉伸和应用标准颜色表（密度分割）。

①在主图像窗口中，选择Tools>ColorMapping>ENVIColorTables。

出现ENVIColorTables对话框，可以使用系统默认的IDL颜色表来调整屏幕的颜色表。

该对话框包括一个灰阶wedge（或彩色wedge，若使用颜色表）和两个滑动块来控制对比度拉伸。

它也有两个下拉菜单：

File和Options。

②选择下列选项之一：

要将颜色表的任何变化自动地应用到打开的图像中，选择Options>AutoApply:

On。

要手动地将变化应用于图像：

A.选择Options>AutoApply:

Off。

B.变化后，选择Options>Apply。

注意：

AutoApply选项自动地设置为8位颜色模式。

③分别移动标签为“StretchBottom”和“StretchTop”的滑动块，来控制被显示的最小值和最大值。

向右移动StretchBottom滑动块，导致图像亮区域变暗；而向左移动StretchTop滑动块，导致图像暗区域变亮。

把拉伸底部设置为最大值，拉伸顶部设置为最小值，拉伸可以被“inverted”。

若AutoApply是打开的，新的对比度拉伸将立即应用于图像。

④通过在所需要的颜色表名上点击，把一个选定的颜色表应用到当前图像。

IDL提供许多预先保存好的颜色表。

“B-Wlinear”表提供一个灰阶图像。

“RAINBOW”颜色表提供一个从“冷”到“热”的密度分割。

其它颜色表选项允许你应用它们首选的颜色方案。

要把颜色表保存为一个ASCII文件：

A.选择File>SaveColorTabletoASCII。

B.输入一个文件名，然后点击“OK”。

要重新设置为初始的颜色表和拉伸，选择Options>ResetColorTable。

要返回到主屏幕并保留被选择的颜色表，选择File>Cancel。

2、交互式密度分割

①在主图像窗口，选择Tools>ColorMapping>DensitySlice，或选择Overlay>DensitySlice。

将出现#nDensitySlice对话框（其中“#n”是用于启动功能的显示号），在“DefinedDensitySliceRanges”下列有八个系统默认范围。

这些范围由滚动窗口计算的最小值和最大值来限定，并显示在“Min”和“Max”文本框中。

②在适当的文本框中输入所需要的最小和最大值，来改变密度分割的范围。

要重新设置数据范围到初始值，点击“Reset”。

③通过选择对话框时底部“Windows”旁所需要的复选框，来选择是否将密度分割颜色应用到图像窗口、滚动窗口或这两个窗口。

④点击“Apply”，将系统默认的范围和颜色应用于图像上。

要编辑数据范围：

A.选择一个数据范围，并点击“EditRange”来改变范围值或颜色。

B.当出现EditDensitySliceRange对话框时，输入所需要的最小和最大值，并从“Color”菜单中选择一种颜色。

C.点击“OK”，执行改变“DefinedDensitySliceRanges”列表中的范围。

D.点击“Apply”，把新的范围和颜色应用到图像上。

要从列表中删除一个范围，选择数据范围，然后点击“DeleteRange”。

要清除密度分割范围列表，点击“ClearRanges”。

File下拉菜单：

File下拉菜单中的选项可以保存所定义的密度分割范围，以及恢复原先保存的范围。

要把当前密度分割范围和颜色保存到一个文件：

在DensitySlice对话框，选择File>SaveRanges。

为了保持一致，键入一个系统默认的扩展名为.dsr的输出文件名。

要恢复以前保存的密度分割范围和颜色：

选择File>RestoreRanges。

选择所需要的输入文件。

要关闭密度分割对话框，并重新设置为初始颜色表，选择File>Cancel。

Options下拉菜单：

Options菜单选项可以添加新范围，设置默认范围，以及设置系统要用到的默认范围数。

要添加范围：

A．在DensitySlice对话框内，选择Options>AddNewRanges。

B．当出现AddDensitySliceRanges对话框时，在适当文本框内输入范围的起始值、末尾值，及所需要的范围数。

C．从“Colors”菜单中选择起始颜色。

多个范围的颜色将遵循图形颜色列表。

D．点击“OK”，返回到DensitySlice对话框。

已经判定的范围将列在对话框中，若需要可以进行编辑。

E．点击“Apply”，把密度分割颜色应用于图像。

要输入所需要使用的系统默认范围数，选择Options>SetNumberofDefaultRanges，然后输入所需要的值。

要重新设置“DefinedDensitySliceRanges”列表为系统默认的范围和颜色，这些由“Min”和“Max”值以及设置的系统默认范围数来决定：

选择Options>ApplyDefaultRanges。

点击“Apply”，把这些范围和颜色应用于显示。

3、控制RGB图像通道（ControlRGBImagePlanes）

①选择Tools>ColorMapping>ControlRGBImagePlanes。

将出现#nRGBChannels对话框，其中n为用于启动功能的显示号。

首次启动时，三个图像颜色通道全是开启的。

②选择一个或多个左边的箭头切换按钮，关闭/开启每个颜色位面。

按钮名将改变，以反映颜色位面设置（“On”或“Off”）。

要颠倒当前的拉伸（一个low-to-high拉伸变为一个high-to-low拉伸，反之亦然），选择一个或多个右边的切换按钮。

按钮名将改变，以反映图像位面设置（“Normal”或“Inverse”）。

要把RGBChannel功能同时应用到所有颜色位面，选择Options>TurnAllChannelsOff/On或InvertAllChannels。

③选择File>Cancel，退出该功能。

二、处理结果：

专题二：

TM与SPOT数据融合

（叙述该过程并处理结果加载到本文档里）

一、处理过程：

1、HSV锐化方法

（1）对lon_tm——基本工具——数据重采样（因为lon_tm分辨率较低，需要乘比例系数2.800397，2.801786）——输出“resizedtm”。

（2）对long_spot（黑白）——基本工具——数据拉伸——outputdatarange0-1之间——输出“stretchedspot”。

（3）变换——图像融合——HSV融合——selectedinputRGB——inputbands选resizedtm红绿蓝波段；HighResolutionInputFile选stretchedspot波段——输出“sharpenedhsv”。

2、彩色标准化锐化方法

（1）对lon_tm——基本工具——数据重采样（因为lon_tm分辨率较低，需要乘比例系数2.800397，2.801786）——输出“resizedtm”。

（2）对long_spot（黑白）——基本工具——数据拉伸——outputdatarange0-1之间——输出“stretchedspot”。

（3）变换——图像融合——Brovey融合——selectedinputRGB——inputforcolorbands选Display#n即可（当前显示重采样tm的窗口）；HighResolutionInputFile选stretchedspot波段——输出“sharpenedhsv”。

二、处理结果：

专题三：

航片的配准与镶嵌

（叙述该过程并处理结果加载到本文档）

一、处理过程：

1、配准：

图像-图像地面控制点（SelectGCPs：

Image-to-Image）

图像-图像配准需要两幅图像均打开。

用每一个显示的缩放窗口选择地面控制点。

可以选择Subpixel（小数的）坐标。

一旦已经选择了足够用于定义一个纠正多项式的控制点，就能预测纠正图像中的GCP位置。

GCPs可以被存储或从文件中恢复，彩色标签与GCP标记的顺序可以更改。

（1）打开和显示图像

用AvailableBandsList打开基图像和纠正图像文件，并在两个窗口显示它们。

①一旦两幅图像都已经显示，选择Map>Registration>SelectGCPs：

ImagetoImage。

②出现ImagetoImageRegistration对话框时，在“BaseImage:

”下面点击需要显示的名字，选择基图像（参照图像）。

③在WarpImage:

下方点击需要显示的名字，选择被纠正的图像。

④点击OK，出现GroundControlPointsSelection对话框。

改变对话框的大小，用鼠标左键点击任何一个角处，拖曳到需要的大小和形状。

（2）选择地面控制点（图像对图像）

选择地面控制点（GCP），通过在缩放窗口定位十字准线，在基图像和纠正图像中选择像元。

像元的左上方是整个数字坐标的原点，自这点向右和向下对应的X值和Y值不断增加。

在缩放窗口中，可得到的像元fraction与缩放系数是成比例的。

例如，缩放系数为4x，像元将被分成4个亚区。

缩放系数为10x时，定位可能只有1/10像元。

GCP标记被放在缩放窗口中，显示出亚像元（subpixel）的位置。

在选择GCPs时，亚像元被用于提高准确度。

（3）在基础和纠正图像中，选择GCP的位置。

①为每幅图像移动缩放窗口到需要的GCP区域。

②在缩放窗口的一个特定像元上点击鼠标左键，把光标定位在该像元或像元的一部分上。

注意：

从图像窗口，点击Tools>PixelLocator（或者使用右键快捷键）选择像元位置，并提供整数的GCP坐标。

在GroundControlPointsSelection对话框，被选择处的坐标按sample、line顺序分别显示在标签为“BaseX,Y”和“WarpX,Y”文本区中。

在缩放窗口处，定位提供了亚像元坐标，因此对话框中显示浮点坐标值。

③一旦两幅图像都选择了需要的像元，在GroundControlPointsSelection对话框中点击“AddPoint”，将选择的GCPs添加到已经选择的X、Y（样本，行）坐标对列表里。

他们将按基图像（第一个圆括号）、纠正图像（第二个圆括号）顺序被列出。

当已经选择了四个或更多个GCPs，对所选纠正图像的预测X、Y坐标；X、Y的误差；RMS误差将别列在GCP列表中。

当GCPs已经被添加到列表中时，在两幅图像的缩放窗口里选择的像元（或亚像元）处画一个标记。

选择的GCP数将出现在标记附近。

标记中心（在十字准线下面）显示了实际的GCP位置。

④用同样的方法添加其它的GCPs。

2、使RMSError最小化

在GroundControlPointsSelection对