遥感图像辐射增强与光谱增强Word格式.docx

遥感图像辐射增强与光谱增强Word格式.docx

《遥感图像辐射增强与光谱增强Word格式.docx》由会员分享，可在线阅读，更多相关《遥感图像辐射增强与光谱增强Word格式.docx（8页珍藏版）》请在冰点文库上搜索。

●Gaussion（高斯拉伸）

●SquareRoot（平方根拉伸）

●Logarithmic（对数拉伸）

（2）主成分分析

（3）樱帽变换。

（4）色彩空间变换

（5）色彩拉伸

2实验方法

将一个多光谱图像打开并显示在视窗中，在主菜单中，选择Display→CustomStretch（或在工具箱中单击

）就可以打开交互式直方图拉伸操作面板。

在进行直方图拉伸时，可选择直方图统计的数据范围为全图统计或者当前视图范围内的数据。

多光谱图像的各个波段之间经常是高度相关的，它们的DN值以及显示出来的视觉效果往往很相似。

主成分分析（PrincipalComponentAnalysis,PCA）就是一种去除波段之间的多余信息、将多波段的图像信息压缩到比原来波段更为有效的少数几个转换波段的方法。

一般情况下，第一主成分（PCI）包含了所有波段中80%的方差信息，前三个主成分包含了所有波段中95%以上的信息量。

由于各波段之间的不相关，主成分波段可以生成更多颜色、饱和度更好的彩色合成图像。

ENVI中提供主成分正变换和主成分逆变换。

当使用主成分正变换是，ENVI可以通过计算新的统计值，或者根据已经存在的统计值进行主成分正变换。

一种特殊的主成分变换。

其中，X=[x1,x2,x3,x4,x5,x7]T,，对应于TM影像的1、2、3、4、5、7波段。

变换后得到：

Y=[y1,y2,y3,y4,y5,y7]T,这六个分量相互正交，前四个分量意义明确：

y1代表亮度分量，是TM六个波段亮度值的加权和，反映了总体的亮度变化；

y2代表绿度分量，与亮度分量垂直，是近红外与可见光波段的对比；

又称为绿度植被指数（GVI）；

y3代表黄度或湿度分量，主要与土壤湿度/水分状况有关；

y4对于MSS数据没有意义，而对于TM数据代表霾，大气辐射衰减效应。

通过对同一幅影像的不同色彩系别之间的转换，从多角度来寻找更适合人眼对颜色的敏感和识别程度，凸显或减弱某些地物的边界的对比，从而帮助对影像的分析。

ENVI支持将三波段红、绿、蓝（RGB）图像变换到一个特定的彩色空间，并且能从所选彩色空间变换回RGB。

两次变换之间包括“色度、饱和度、颜色亮度值（HSV）”、“色度、亮度、饱和度（HLS）”、HSV（USGSMunsell）。

其中，色度代表像元的颜色，取值范围为0-360；

饱和度代表颜色的纯度，取值为0-1；

颜色亮度值为颜色的亮度，取值范围为0-1；

亮度表示整个图形的明亮程度，取值范围为0-1。

其中，HSV（USGSMunsell）彩色系统被土壤科学家和地质学家用于描述土壤和岩石的颜色特征。

这套彩色系统已经被美国地质勘测部门做了修订，以描绘数字图像的颜色。

USGSMunsell变换将RGB坐标变换成了色彩坐标（色调、饱和度和颜色亮度值）。

色度变化范围为0-360，这里0和360代表蓝，120代表绿，240代表红。

饱和度变化范围为0-208，值越高代表颜色越纯。

颜色亮度值的变化范围大致为0-512，较高的值代表较量的颜色。

ENVI提供的色彩拉伸有：

去相关拉伸、Photographic拉伸和饱和度拉伸。

三、实验设备与材料

1实验设备

装有ENVI的计算机

2实验材料

资源一号02C多光谱数据

南京市Landsate5号卫星上传感器获得的一景编号为LT51200382010231BJC00的TM影像

四、实验步骤

1交互式互式直方图拉伸

（1）打开一个多光谱图像并显示在视窗中，在主菜单中，选择Display→CustomStretch（或在工具箱中单击

），打开交互式直方图拉伸操作面板。

（2）在面板右侧可以根据需要选择R，G，B3个波段之一，这里以R波段为例。

（3）在交互式直方图拉伸操作面板下方的下拉列表中，可以选择不同的拉伸方法。

图1交互式直方图拉伸操作面板

单击面板右下方

按钮，可以恢复拉伸方法到打开此面板之前的状态。

2主成分分析

1）打开图像文件

2）在Toolbox工具箱中，双击Transform/PCARotation/ForwardPCARotationNewStatisticsandRotate工具。

在PrincipalComponentsInputFile对话框中，选择图像文件

3）在ForwardPCAParameters对话框中，在StatsX/YResizeFactor文本框中键入小于等于1的调整系数，用于计算统计值时的数据二次采样。

键入小于1的调整系数，将会提高统计计算速度。

例如使用一个0.1的调整系数，在统计计算中将只用到十分之一的像元。

选择默认值为1。

图1主成分分析对话框

4）键入一个输入统计路径及文件名。

使用箭头切换按钮，选择CovarianceMatrix（协方差矩阵）或根据CorrelationMatrix（相关系数矩阵）计算主成分波段。

一般来说，计算主成分时选择使用协方差矩阵。

当波段之间数据范围差异较大时，选择相关系数矩阵。

5）选择输出的路径及文件名，输出数据类型为FloatingPoint。

6）单击SelectSubsetfromEigenvalues标签右侧的箭头切换按钮，选择“EYS”。

统计信息将被计算，并出现SeleOutputPCBands对话框，列出每个波段机器相应的特征值；

同时，也会列出每个主成分波段中包含的数据方差的累计百分比。

如果选择“NO”,则系统会计算特征值并显示供选择输出的波段数。

7）输出波段数（NumberofOutputPCBands）选择默认值（输入文件的波段数）。

8）单击OK

ENVI处理完毕后，将会出现PCEigenvalues绘图窗口。

可以看到，第一、二、三分量具有很大差异的特征值。

主成分波段被导入可用波段列表中，选择PC1,PC2,PC3合成RGB显示，色彩非常饱和。

选择PC4,PC5,PC6合成RGB显示，可以看到很多噪声。

在Toolbox工具箱中，双击Statistics/ViewStatisticsFile工具，打开主成分分析中得到的统计文件，可以得到各个波段的基本统计值、协方差矩阵、相关系数矩阵和特征向量矩阵。

在Toolbox工具箱中，双击Transform/PCARotation/InversePCARotation工具可以执行主成分逆变换。

图2主成分分析特征值窗口

图3主成分分析得到的个波段统计值

3樱帽变换

1）打开一个Landsat5TM数据

2）在Toolbox工具箱中，双击Transform/TassaledCap工具。

在TasseledCapTransformInputFile对话框中，选择输入Landsat5TM数据，单击OK

3）在TasseledCapTransformParameters面板中。

选择InputFileType选项中的”Landsat5TM”

4）选择输出的路径及文件名，单击OK.

图4TasseledCapTransformParameters面板

4色彩空间变换

色彩空间变换操作比较简单，下面以从RGB到HSV为例介绍操作过程，其他色彩空间变换操作过程类似。

1打开一个至少含3个波段的图像文件，显示RGB彩色图像。

2在Toolbox工具箱中，双击Transform/ColorTrandsform/RGBtoHSVColorTrandsform工具。

在RGBtoHSVInputBands对话框中，可以从可用波段列表中选择3个波段进行变换。

3在RGBtoHSVParameters对话框中，选择输出路径及文件名，单击OK。

、

5色彩拉伸

1）去相关拉伸

去相关拉伸处理可以消除多光谱数据各波段间的高度相关性，从而生成一幅色彩亮丽的彩色合成图像。

它首先是对图像做主成分分析，并对主成分图像进行对比度拉伸处理，然后再进行主成分逆变换，将图像恢复到RGB彩色空间，达到图像增强的目的。

在Toolbox工具箱中，双击Transform/DecorrelationStretch工具。

在DecorrelationStretchInputBands对话框中，从可用波段列表中选择3个波段作为输入。

在DecorrelationStretchParameters对话框中选择输出路径及文件名。

2）Photographic拉伸

Photographic拉伸可以对一幅真彩色输入图像进行增强，从而生成一幅与目视效果良好吻合的RGB图像。

其结果与现实彩色照片类似。

这种拉伸方法对真彩色输入图像的波段进行非线性缩放，然后将它们叠加。

在Toolbox工具箱中双击Transform/PhotographicStretch工具。

在RGBPhotographicStretchInputBands输入对话框中，从可选用波段列表中选择3个波段作为输入。

在PhotographicStretchParameters对话框中选择输出路径及文件名。

3）饱和度拉伸

饱和度拉伸是对输入的3波段图像进行彩色增强，生成具有较高颜色饱和度的波段。

输入的数据由红、绿、蓝（RGB）空间变换为色度、饱和度和颜色亮度值（HSV）空间。

对饱和度波段进行高斯拉伸，从而是数据分布到整个饱和度范围，然后逆变换回RGB空间，完成增强处理。

在Toolbox工具箱中双击Transform/SaturationStretch工具。

在SaturationStretchInputBands对话框中，从可选用波段列表中选择3个波段作为输入。

在SaturationStretchParameters对话框中选择输出路径及文件名。

WelcomeTo

Download!

!

欢迎您的下载，资料仅供参考！