ENVI影像增强处理66967.docx

《ENVI影像增强处理66967.docx》由会员分享，可在线阅读，更多相关《ENVI影像增强处理66967.docx（18页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、ENVI影像增强处理66967I ENVI 初步学习和影像增加一试验目的结合 ENVI 软件中文指导书，学习和驾驭 ENVI 软件的基本操作，学会运用ENVI 软件对遥感影像进行分析增加处理，加深对遥感图像处理试验这门课的学习，为以后的实践工作打好坚实的基础二试验数据介绍该次试验所运用数据为甘肃舟曲灾前遥感影像，图像名称为 0112024.img。三试验过程. . 软件的基本操作ENVI 主工具条： 1图像的输入与输出.首先启动ENVI,选择File-Open Image File,出现Enter Data Filename对话框，选择文件的正确路径，0112024.img，再点击OK打开文件

2、。2在打开的Available Band List菜单中，可以显示图像的各个波段的基本信息，其中Gray Scale为灰色显示，RGB Color为彩色合成，并且可以选择彩色合成的波段，单击Load就可显示图像，打开的图像由三部分组成：Scroll（滚动）窗口、主图像窗口、以及缩放（Zoom）窗口,可以运用多个显示窗口组，组中每个窗口的大小都可以调整。其中菜单中的New Diaplay可以打开一个新的图像。3在Available BandsList菜单中选择Available Files List还可以显示出遥感图像的基本信息，详细如图所示：— 欢迎下载 2 4.若要保存图像，须要

3、在图像所在窗口中选择File-Save Image As-Image File，弹出Output Display to Image File 对话框。对于单波段图像，选择 8-bit Color，而多波段彩色合成图像则选择24-bit Color。图像的保存方式有两种：一是干脆保存为文件；二是选Memory，记忆在Available BandsList菜单中。. . 图像增加与变换一 空间域增加 对单个像元的灰度值进行变换进行增加处理。1. 线性变换 线性拉伸：线性拉伸的最小和最大值分别设置为0和255，两者之间的全部 其 它 值 设 置 为 中 间 的 线 性 输 出 值 ， 在 主 图 像

4、 窗 口 ， 选 择Enhance->Interactive stretch，弹出对话框：（1）线性拉伸：Stretch_Type -> Linear 选择Options->Auto Apply，打开自动应用功能。— 欢迎下载 3 （2）分段线性拉伸：分段线性对比度拉伸可以通过运用鼠标在输入直方图中放置几个点进行交互地限定。当在点之间供应线性拉伸时，线段在点处连接起来，选择Stretch_Type > Piecewise Linear。— 欢迎下载 4 2. 非线性变换 高斯拉伸：选择Stretch_Type->Gaussian。输入拉伸的

5、最小和最大值，要手动地输入所须要的标准差值，选择Options-> SetGaussian Stdv。 设置高斯标准差3. 平方根拉伸 选择Stretch_Type->Square Root。输出直方图用一条红色曲线显示平方根函数，被拉伸数据的分布呈白色叠加显示。— 欢迎下载 54. 直方图均衡化 选择Stretch_Type->Equalization。输入直方图显示未被修改的数据分布。输出直方图用一条红色曲线显示均衡化函数，被拉伸数据的分布呈白色叠加显示。 5. 直方图规定化选择Stretch_Type > Arbitrary.通过点击或按住鼠标左键绘制

6、输出直方图的线段，在Output Histogram 窗口内绘制输出直方图。随意的直方图将用绿色来显示。输出直方图用红色显示你的直方图，匹配的数据函数用白色曲线。— 欢迎下载 6二 频率域增加 1. 傅里叶变换 傅里叶变换：傅立叶分析是一种将图像分成空间上各种频率成分的数学方法。事实上，快速的傅立叶变换被用来将数据变换成一个困难的强调频率分布的图像。ENVI 中FFT滤波包括图像正向的FFT、频率滤波器的交互式建立、滤波器的应用，以及FFT向原始数据空间的逆变换。当前，FFT 处理没有用到ENVI tiling 程序，因此能被处理的图像大小受到系统可利用内存的限制。FFT图像是复数

7、数据类型，它占用了类似大小的字节图像的8倍内存。快速傅里叶变换：正向的FFT 生成的图像能显示水平和垂直空间上的频率成分。图像的平均亮度值显示在变换后图像的中心。远离中心的像元代表图像中增加的空间频率成分。这一滤波能被设计为消退特别的频率成分，并能进行逆向变换。（1）打开一个TM图像 （2）在主菜单中，选择Filters->FFT Filtering->Forward FFT。在filter中选择FFT filtering，然后选择forward FFT进行快速傅里叶变换。— 欢迎下载 7（3）定义FFT滤波器 在主菜单中，选择Filters->FFT Filte

8、ring->Filter Definition。Circular Pass 为低通滤波，Circular Cut为高通滤波器，在Radius框中输入滤波半径，50；— 欢迎下载 8 （4）方向FFT变换 在主菜单中，选择Filters->FFT Filtering->Invert FFT。选择FFT变换后的图像 选择 Filter Definition 图像— 欢迎下载 9低通滤波后的结果 2. 频率域锐化 先对影像进行FFT变换（运用高通滤波），之后对得到的进行Inverse FFT变换进而得到锐化处理的影响。高通滤波器可使信息源的高频通过而对低频加

9、以抑制。由于抑制了反映灰度骤变的边缘特征的低频信息及包含在低频中的孤立点— 欢迎下载 10 噪声，所以高通滤波起了锐化图像的处理作用。三 彩色增加 1. 伪彩色增加 将亮度值等间隔分割分别给予不同的色调，合成处理的过程。在主图像窗口，选择overlay > Density Slice。将出现#1 Density Slice对话框，在Defined Density Slice Ranges下列有八个系统默认范围。在适当的文本框中输入所须要的最小和最大值，来变更密度分割的范围。如要重新设置数据范围到初始值，点击Reset。— 欢迎下载 112. 假彩色增加 将不同波段

10、的影像分别给予不同的色调，合成处理的过程。如分别给予TM图像1、2、3波段色调R、G、B；从Available Bands List内，选择RGB Color切换按钮。在序列中点击所须要显示的红、绿和蓝波段名。将波段名称导入R.G.B后载入图像得到合成的假彩色图像。— 欢迎下载 12 321 ：真彩色合成，即 3、2、1 波段分别给予红、绿、蓝色，则获得自然彩色合成图像，图像的色调与原地区或景物的实际色调一样，适合于非遥感应— 欢迎下载 13 用专业人员运用。432 ：标准假彩色合成，即 4、3、2 波段分别给予红、绿、蓝色，获得图像植被成红色，由于突出表现了植被的特征

11、，应用非常的广泛，而被称为标准假彩色。举例：卫星遥感图像示蓝藻暴发状况我们先看一看蓝藻爆发时遥感监测机理。蓝藻暴发时绿色的藻类生物体拌随着白色的泡沫状污染物聚集于水体表面，蓝藻覆盖区的光谱特征与四周湖面有明显差异。由于所含高叶绿素 A 的作用，蓝藻区在 LandsatTM2 波段具有较高的反射率，在 TM3 波段反射率略降但仍比湖水高，在 TM4 波段反射率达到最大。因此，在 TM4（红）、3（绿）、2（蓝）假彩色合成 图像上，蓝藻区呈绯红色，与四周深蓝色、蓝黑色湖水有明显区分。此外，蓝藻暴发聚集受湖流、风向的影响，呈条带延长，在 TM 图像上呈条带状结构和絮状纹理，与四周的湖水面也有明显不同

12、。451 ：信息量最丰富的组合，TM 图像的光波信息具有 34 维结构，其物理含义相当于亮度、绿度、热度和湿度。在 TM7 个波段光谱图像中，一般第 5个波段包含的地物信息最丰富。3 个可见光波段（即第 1、2、3 波段）之间，两个中红外波段（即第 4、7 波段）之间相关性很高，表明这些波段的信息中有相当大的重复性或者冗余性。第 4、6 波段较特别，尤其是第 4 波段与其他波段的相关性得很低，表明这个波段信息有很大的独立性。计算各种组合的熵值的结果表明，由一个可见光波段、一个中红外波段及第 4 波段组合而成的彩色合成图像一般具有最丰富的地物信息，其中又常以 4，5，3 或 4，5，1 波段的组

13、合为最佳。第 7 波段只是在探测森林火灾、岩矿蚀变带及土壤粘土矿物类型等方面有特别的作用。最佳波段组合选出后，要想得到最佳彩色合成图像，还必需考虑赋色问题。人眼最敏感的颜色是绿色，其次是红色、蓝色。因此，应将绿色给予方差最大的波段。按此原则，实行 4、5、3波段分别赋红、绿、蓝色合成的图像，色调反差明显，层次丰富，而且各类地物的色调显示规律与常规合成片相像，符合过去常规片的目视判读习惯。例如把 4、5 两波段的赋色对调一下，即 5、4、3 分别给予红、绿、蓝色，则获得近似自然彩色合成图像，适合于非遥感应用专业人员运用。741 ：波段组合图像具有兼容中红外、近红外及可见光波段信息的优势，图面色调

14、丰 富，层次 感好，具有极为丰富的地质信息和地表环境信息；而且清楚度高，干扰信息少，地质可解译程度高，各种构造形迹（褶皱及断裂）显示清晰，不同类型的岩石区边界 清楚，岩石地层单元的边界、特别岩性的展布以及火山机构也显示清晰。742 ：1992 年，完成了桂东南金银矿成矿区遥感地质综合解译，利用 1：10万 TM7、4、2 假彩色合成片进行解译，共解译出线性构造 1615 条，环形影像 481 处, 并在总结了构造蚀变岩型、石英脉型、火山岩型典型矿床的遥感— 欢迎下载 14 影像特征及成 矿模式的基础上，对全区进厅成矿预料，圈定金银 A 类成矿远景区 2 处，B 类 4 处，C 类

15、5 处。为该区优选找矿靶区供应遥感依据。743：我国利用美国的陆地卫星专题制图仪图像胜利地监测了大兴安岭林火及灾后改变。这是因为 TM7 波段（2.08-2.35 微米）对温度改变敏感；TM4、TM3 波段则分别属于红外光、红光区，能反映植被的最佳波段，并有削减烟雾影响的功能；同时 TM7、TM4、TM3（分别给予红、绿、蓝色）的彩色合成图的色调接近自然彩色，故可通过 TM743 彩色合成图的分析来指挥林火扩散与限制和灾后林木的复原状况。754：对不同时期湖泊水位的改变，也可采纳不同波段，如用陆地卫星 MSS7，MSS5，MSS4 合成的标准假彩色图像中的蓝色、深蓝色等不同层次的颜色得以区分。

16、从而可用作分析湖泊水位改变的地理规律。3. 彩色变换 在图像处理中通常应用的有两种彩色坐标系统。一种是由红（R）、绿（G）、蓝（B）三原色构成的彩色空间或坐标系统，这是我们之前所讲的彩色改变和增加处理的基础；另一种是彩色空间是由色调（H）、饱和度（S）和明度（I）三个变量构成的。彩色系统变换主要是指这两种坐标系统之间的变换。详细操作如下：RGB向HSV的转换：在主菜单中点击Transforms->Color Transforms（在做正变换的选项中有RGB to HSV、RGB to HLS 和RGB to HSV（USGS Munsell）三种方法，点击RGB to HSV。接着选择已

17、经彩色合成好的遥感图像，接着在弹出的对话框中选择display1，出现RGB to HSV（HLS、USGS Munsell）Parameters 对话框，选择保存路径后，软件就自动完成了HSV变换。H(色度)— 欢迎下载 15S(饱和度) V(亮度)— 欢迎下载 16四 多图像代数运算 Band Math 功能供应了一个敏捷的图像处理工具，其中很多功能是无法在任何其它的图像处理系统中获得的。由于每个用户都有独特的需求，因此该功能允许用户自己定义处理算法，并将之应用到打开的波段或整个图像中，可以依据须要自定义简洁或困难的处理程序。选择Basic Tools > B

18、and Math.在Enter an expression 的文本框内，输入变量名和所须要的数学运算符。变量名必需以字符b或B开头，后面跟着5个以内的数字字符。1. 差值法：选择 Basic Tools > Band Math.在Enter an expression 的文本框内输入 B1-B2，点击 add to list,确定后定义 B1，B2 为波段 b1,b2,选择 memory 后确定。— 欢迎下载 172.比值法：选择 Basic Tools > Band Math.在Enter an expression的文本框内输入 float(B2)/float(B3

19、)，点击 add to list,确定后定义 B2，B3 为波段 b2,b3,选择 memory 后确定。— 欢迎下载 18 3.混合运算法：选择 Basic Tools > Band Math.在Enter an expression 的文本框内输入 float(B2)/float(B3)+B4，点击 add to list,确定后定义 B2，B3，B4 为波段 b2,b3,b4 选择 memory 后确定。— 欢迎下载 19五 多光谱图像变换 1. 主成分变换 主成分分析（PCA）用多波段数据的一个线性变换，变换数据到一个新的坐标系统，以使数据的差异达到最大。

20、这一技术对于增加信息含量、隔离噪声、削减数据维数特别有用。由于多波段数据常常是高度相关的，主成分变换找寻一个原点在数据均值的新的坐标系统，通过坐标轴的旋转来使数据的方差达到最大，从而生成互不相关的输出波段。主成分（PC）波段是原始波谱波段的线性合成，它们之间是互不相关的。可以计算输出主成分波段。第一主成分包含最大的数据方差百分比，其次主成分包含其次大的方差，以此类推，最终的主成分波段由于包含很小的方差（大多数由原始波谱的噪声引起），因此显示为噪声。由于数据的不相关，主成分波段可以生成更多种颜色的彩色合成图像。在ENVI主菜单中选择Transforms->Principal Compone

21、nts->Forward PC Rotation->Compute New Statistics and Rotate，选择0112024.img作为输入数据。出现如下Forward PC Rotation Parameters对话框：— 欢迎下载 20 该对话框参数设置如下：在 Stats X/Y Resize Factor 文本框键入小于 1 的调整系数，对计算统计值的数据进行二次抽样。键入一个小于1的调整系数，以提高统计计算的速度。例如，在统计计算时，用一个0.1的调整系数将只用到非常之一的像元。若须要，键入一个输出统计文件名。点击按钮，选择基本Covarianc

22、e Matrix或Correlation Matrix。计算主成分时，有代表性地要用到协方差矩阵。当波段之间数据范围差异较大时，要用到相关系数矩阵，并且须要标准化。选用File或Memory输出。若选择输出到File，在标有Enter Output Filename的文本框里键入要输出的文件名；或用Choose按钮选择一个输出文件名。从Output Data Type菜单里，选择须要的输出类型（字节型，整型，无符号整型，长整型，无符号长整型，浮点型，双精度型）。用下列选项，选择输出 PC 波段数。限定输出 PC 波段数，键入须要的数字，或用Number of Output PC Bands标签

23、旁边的按钮确定输出的PC波段数。默认的输出波段数等于输入的波段数。通过检查特征值，选择输出的 PC 波段数。A 点击Select Subset from Eigenvalues标签旁边的按钮，选择YES。特征值将被计算，出现 Select Output PC Bands对话框，列表显示着每一个波段和其相应的特征值。同时也为全部波段显示出每个波段中包含的数据改变的累积百分比。— 欢迎下载 21 B 在 Number of Output PC Bands 文本框里，键入一个数字或点击按钮，确定输出的波段数。特征值大的 PC 波段包含最大量的数据差异。较小的特征值包含较少的数据信息和较多

24、的噪声。有时，为存储磁盘空间，最好仅仅输出特征值大的那些波段。C 在 Select Output PC Bands 对话框里，点击OK。选择Memory 输出，其他的默认，输出7个主成分波段，自动加载到波段中同时出现一个主成分特征值的plot 图：在该对话框中选择 Edit>Data Values，查看各成分的特征值的大小：— 欢迎下载 22特征值反映了各主成分所占信息量的大小，从上到下渐渐削减，表明各波段信息量渐渐削减，可以据此计算各成分信息量占总信息量的百分比，例如，从上面 数 据 可 以 得 知 ， PC1 图 像 所 占 的 信 息 量 为 740.9951/（740

25、.9951+178.3992+29.4076+21.9403+7.6024+2.1608+0.7218）×100%=75.52%。由于主成分图像有 7 个波段，选取了有代表性的前三个波段显示如下：PC1: PC2: PC3:— 欢迎下载 23检验主成分与原始数据之间的定量关系：下面来检验各主成分究竟与原始的7个波段数据之间究竟是什么定量关系，这对说明各主成分的含义很有用。主成分变换的系数矩阵是原始数据协方差矩阵的特征向量矩阵的一个转置，我们先通过统计计算得到特征向量矩阵如下（方法参照前面所述），留意数据源选择原始数据0112024.img：— 欢迎下载 24

26、— 欢迎下载 252. 缨帽变换 缨帽变换是一种通用的植被指数，可以被用于Landsat MMS 或Landsat TM 等数据。对于Landsat MMS 数据，缨帽变换将原始数据进行正交变换，变成四维空间（包括土壤亮度指数SBI、绿色植被指数GVI、黄色成分（stuff）指数YVI，以及与大气影响亲密相关的non-such 指数NSI）。对于Landsat TM 数据，缨帽植被指数由三个因子组成——亮度、绿度与第三重量（Third）。其中的亮度和绿度相当于MSS 缨帽的SBI和GVI，第三种重量与土壤特征有关，包括水分状况。对于Landsat 7 ETM

27、 数据穗帽植被指数由六个因子组成-亮度、绿度、湿度、以及第四、第五、第六几种重量。选择Transforms > Tassled Cap。出现Tasseled Cap Transformation Input File对话框时，选择0112024.img文件。用下拉菜单，选择Input File Type（ Landsat 5TM 、LandsatMSS、Landsat 7 ETM ）。— 欢迎下载 26ENVI 将缨帽波段名输入到Available Bands List中，在那里可以用标准ENVI 灰阶或RGB 彩色合成方法显示。本文来源：网络收集与整理，如有侵权，请联系作者删除，谢谢！第18页 共18页第 18 页 共 18 页第 18 页 共 18 页第 18 页 共 18 页第 18 页 共 18 页第 18 页 共 18 页第 18 页 共 18 页第 18 页 共 18 页第 18 页 共 18 页第 18 页 共 18 页第 18 页 共 18 页