ENVI软件学习实验报告.docx

1、ENVI软件学习实验报告遥感实验报告姓名 郭美辰第一章1.1、时间：1.2、内容 ：在这一阶段，对遥感软件的基本功能加以熟练。1.3、流程图：第二章 几何校正2.1、报告人：郭美辰2.2、实验目的：把图像校正成便于解译的图像2.3、实验原理： 通过envi软件，把要校正的图像配准，校正。2.4、实验时间：2.5、实验流程2.6、实验数据： 一幅参考图像，一幅校正图像2.7、使用软件：ENVI2.8、步骤：选择两幅图像，把参考图像作为基准地图，来校正图像。1 选择MapRegistrationSelect GCPs: image to image.在对话框里选择图像，点击OK。将出现Ground

2、 Control Points Selection 对话框，分别选择控制点，最少4个，多选精度回高些。然后点击add point .又分别尝试了其他几项，比如提高精度，多选几个点，在show list 之后，选择OptionsOrder Points by Error,在列表中 删除误差较高的点，我选了19个点，删除了三个误差较大的点2 把在这些点保存，如下图所示： 图2-1 配准点3 选择Map RegistrationWarp from GCPs, 选择要被校正的图像，然后点击OK4 选择基准图像，如下图：5 选择保存路径，校正即可完成。2.9、实验结果：校正完成。校正误差在0.5以内，完

3、成情况比较好。2.10、讨论与分析： 配准的时候，选一个点就看下精度；还有自己一定要清楚哪个是被校正的图像，哪个是参考图像。随时注意误差，如果过大，应该马上重选；一开始的四个点应该分布范围较广，最好是四个边的点，这样误差会减小。用了配准中的MapRegistrationSelect GCPs: image to image.就看到在registration 下拉框中的 image to map ，应该研究下 这个校正的方法。试了一下，暂时没成功，不太会做。第三章 镶嵌3.1、报告人： 郭美辰3.2、实验目的： 完成影像的镶嵌3.3、实验原理： 当研究区超出单幅遥感图像所覆盖的范围时，通常需要将

4、两幅或多福图像拼接起来，形成一幅或一系列覆盖全区的较大图像，这个过程即图像镶嵌。3.4、实验时间： 1018-10233.5、实验流程3.6、实验数据：实验数据为多光谱数据和全色数据，多光谱数据空间分辨率为2米；全色波段空间分辨率为0.5米。3.7、使用软件：ENVI3.8、步骤：选择MapMosaickingGeoreferenced，打开Map Based Mosaic镶嵌界面。ImportImport Files选择要镶嵌的影像，定义镶嵌范围（行列号、影像拉框选择、地理坐标，引入其他文件范围等），键入镶嵌输出的维数（或尺寸），选择镶嵌所用的波段。分别选择影像右键点击Edit Entry：

5、对背景值、羽化值、匀光等进行设置。所有参数设置完成以后FileSave Template：生成临时效果影像，方便检查、调整。生成最终结果：FileApply确定输出象素大小、重采样方式、输出路径及文件名、背景值等OK。3.9、实验结果：。3.10、讨论与分析：对每幅图都可以设置其背景，羽化值等，另外在做的过程中，要考虑用不用把fethering distance 设置，有时候有可能会在镶嵌的时候，带一条边。第四章 融合4.1、报告人： 郭美辰4.2、实验目的： 将单一传感器的多波段信息或不同类别传感器所提供的信息加以综合，消除多传感器信息之间可能存在的冗余和矛盾，加以互补，改善遥感信息提取的及

6、时性和可靠性，提高数据的使用效率。4.3、实验原理：遥感分波段分图像记录了地物波谱的微弱差异，显然充分利用地物在不同波段的差异，可以更有效地识别地物。遥感影像融合是对不同空间分辨率遥感影像的融合处理，使处理后的遥感影像既具有较高的空间分辨率（高空间分辨率数据），同时又具有多光谱特征（较低分辨率数据），从而达到图像增强的目的。4.4、实验时间：1018-10234.5、实验流程4.6、实验数据：实验数据为多光谱数据和全色数据，多光谱数据空间分辨率为2米；全色波段空间分辨率为0.5米。4.7、使用软件：ENVI4.8、步骤：一般来说，遥感影像的数据融合分为预处理和数据融合两步。1、预处理：主要包括

7、遥感影像的几何纠正、大气订正、辐射校正及空间配准。（1）几何纠正、大气订正及辐射校正的目的主要在于去处透视收缩、叠掩、阴影等地形因素以及卫星扰动、天气变化、大气散射等随机因素对成像结果一致性的影响； （2）影像空间配准的目的在于消除由不同传感器得到的影像在拍摄角度、时相及分辨率等方面的差异。2、数据融合根据融合目的和融合层次智能地选择合适的融合算法，将空间配准的遥感影像数据（或提取的图像特征或模式识别的属性说明）进行有机合成，得到目标的更准确表示或估计 。对于各种算法所获得的融合遥感信息，有时还需要做进一步的处理，如“匹配处理”和“类型变换”等，以便得到目标的更准确表示或估计。ENVI中融合的

8、方法有很多种，典型的有HSV、PCT（主成分变换）、K-T（缨帽变换）、小波变换、HIS变换等。4.8.1 自动融合主菜单中选择Basic Tools-Resize Data(Spatial/Spectral)，选择低分辨率图像（多光谱影像），单击OK。在Resize Data Parameters对话框中Samples和Lines中直接输入高分辨率影像（全色波段影像）的大小值，使得两幅影像的大小完全一致，重采样的方法最近邻法，键入输出文件名。如图1。图1 重采样打开重采样之后的图像，在主菜单中选择Transform - Image Sharpening -HSV，选择重采样之后的图像，再选择

9、全色波段作为高分辨输入图像，设置重采样方式，输出融合后的影像，如图2。图2 融合时影像在两幅有相同地理坐标系统的情况下，该融合方法不需要在融合之前调整两幅图的分辨率一致，尺寸一致，ENVI会自动完成这一过程，输出图像的分辨率与高分辨率图像保持一致。4.8.2 手动融合Transform - Color Transforms - RGB to HSV，然后选择重采样过的多光谱遥感影像作为输入的RGB影像，输入要输出的文件名，如图4-3，并显示变换后的影像。Transform - Color Transforms - RGB to HSV，然后选择重采样过的多光谱遥感影像作为输入的RGB影像，输入

10、要输出的文件名。拉伸全色波段影像并代替HSV影像的亮度值波段。从ENVI主菜单中选择Basic Tools- Stretch Data，选择全色波段文件，然后点击ok。在弹出的Data Stretching对话框的output Data部分中，在Min文本框中输入0，在Max中输入1。并输入一个输出文件名，点击ok将全色波段数据拉伸为浮点型，范围在0到1之间。在ENVI主菜单选择Transform - Color Transforms - HSV to RGB，在弹出的HSV to RGB Parameters对话框中输入要输出的文件名，点击ok进行反变换，如图。图 3 经过HSV变换后的图像

11、4.9、实验结果：图4 自动融合和手动融合对比图4.10、讨论与分析：一般来说，图像融合过程分为三个层次，即预处理、信息融合与应用层。预处理就是进行几何校正，噪声消除及配准。图像融合也分为三级，分别是像素级融合，特征级融合，决策级融合。4、11 其他方法4.11.1 PC波谱锐化 Transform - Image Sharpening - PC Spectral Sharpening依次选择多光谱文件、高分辨率文件，设置重采样模式，生成图像，如图5。4.11.2 Gram-Schimidt 波谱锐化Transform - Image Sharpening - Gram-Schimidt Sp

12、ectral Sharpening 依次选择多光谱文件、高分辨率文件，设置重采样模式，生成图像，如图5。4.11.3 BroveyTransform - Image Sharpening - Color Normalized（Brovey）依次选择多光谱文件、高分辨率文件，设置重采样模式，生成图像，如图5。图 5 三种方法的融合图像讨论与分析HIS、brovey 等方法容易使融合后的图像失真，而基于Gram-Schimidt算法的图像融合方法既能使融合后影像保真度比较好，而且计算又较为简单。第五章 分类5.1、报告人：郭美辰5.2、实验目的： 分别用几种方法，完成图像分类。5.3、实验原理：将

13、图像中每个像元根据其在不同波段的光谱亮度、空间结构特征或者其他信息，按照某种规则或者其他信息，按照某种规则或算法划分为不同的类型。5.4、实验时间：5.5、实验流程5.6、实验数据：融合数据5.7、使用软件：ENVI5.8、步骤：ENVI实验步骤非监督分类ClassificationunsupervisedIsoData,设置其参数，一步就可以对其分类，对其加个图例，分类结果如下：监督分类通过人工目视解译，判定大概分类的个数。打开图像后，选择overlay ROI 开始选择感兴趣区。通过设置，可以设置每个类的颜色，名称。如图所示：选择完感兴趣区，optioncalculate ROI sepa

14、rability,计算感性趣的分离度，选择要查看的图像，出现对话框，选择要查看的感性趣区。如图所示：显示分离度图像如图所示，分离度小于1.4的时候，视为分离度很小，可以将其合并。在optionmerge ROI ,出现如下对话框，选中，合并。Classificationsupervisedmaximum likelihood.开始分类。分类结果如下：面向对象 I zoom 中的，processing下拉菜单中的feature extraction 中，进行面向对象的特征提取。这样，相同的地物就被分割开来。决策树分类主菜单classificationdecision treebuild new

15、decision tree单击node1,在node1中，输入特定的函数名称，在expression中输入具体的运算过程。如图所示： 点击OK，出现如下对话框。在对话框中单击ndvi，可以选择要分类的影像。 选择影像后，分别单击class1 和 class2，可以输入其名称，可以选择该类别的颜色。如图：选好后，可以单击某个class右键，选择add childen ，继续增加条件。一般规则的制定如下：Class1(朝北缓坡植被):NDVI0.3, slope20, aspect270Class2(非朝北缓坡植被):NDVI0.3, slope20, 90=aspect0.3, slope=20

16、Class4(水体):NDVI=0.3, 0b420Class5(裸地):NDVI=20Class6(无数据区，背景):NDVIPrinciple ComponentsForward PC Rotation6.8、实验结果：6.9、讨论与分析：图像变换主成分变换最小噪声分离颜色空间变换（HSV，HLS）缨帽变换图像拉伸去相关，饱和度，彩色直方图生成特征值的分布图像，信息主要分布在前面的波段中，并且，特征值越大，它所涵盖的信息就越大。看第一个波段的图像：再看一下最后一个波段的图像：（信息很少，大部分都是噪声了）最后合成图像如下：在主成分变换后，其实基本就可以目视解译了。第七章 图像三维可视化7.

17、1、报告人：郭美辰7.2、实验目的：进行图像三维可视化，便于识图。7.3、实验原理：通过dem，使图像成立体状。7.4、实验时间：7.5、实验数据： 一幅多波段的遥感影像7.6、使用软件： Envi7.7、步骤：加载两幅图像，其中一幅为必须dem影像，这样才能进行三维可视化。在图像的菜单栏下tools3D Surface View，弹出的对话框中要导入一幅DEM影像。选中那幅dem 的影响。在OK之后的对话框中可以对其属性，重采样的方法等进行编辑。编辑后点击OK。7.8、实验结果：7.9、讨论与分析：单击optionssurface controls 下，可以对其位置，表面的类型（纹理，特征等

18、），dem 的分布大小等进行设置。单击optionschange back grond color 修改其背景颜色。如图所示： 设置飞行，optionsmotion controls 下，可以选点，单击add, 依次重复，多选几个点。如下图所示：然后单击play sequence ,在file菜单下，可以选择save sequence as MPEG ,然后选择路径，生成一个录像。第八章 物质制图8.1、报告人：郭美辰8.2、实验目的：了解物质制图，知道如何应用。8.3、实验原理：通过一幅多波段的影像，构建多元光谱。8.4、实验时间：8.5、实验数据： 一幅多波段的遥感影像8.6、使用软件：

19、Envi8.7、步骤：打开一幅多光谱图像，构建多元多谱。打开波谱库，spectralspectral libraries spectral library Viewer ,在如下对话框中选择中下方的open，再选择spectral library 波谱库。8.8、实验结果：在波谱库下的菜单下，选一个植被的波谱库，来据此划分植物类型。打开后，点击OK。出现的对话框中，可以依次选择想要分类的植被。我选了五种，包括草地，杨树等。选择好后，单击鼠标右键plot key,感兴趣可以多选几种，可以改颜色。现在即可对其进行分类，有几种方法，最常用的就是波谱角分类，在classificationsupervi

20、sedspectral angle mapper ,对其进行分类。在此对话框中，选择要分类的图像，即原始图像。点击OK，出现一个新的对话框。因为上述过程中，plot窗口是打开的，所以在对话框中，直接导入改窗口，在importfrom plot windows。Select all，把想分类的植被都选入其中，单击OK，出现如下对话框。直接应用，apply，会出现对话框，可以分别用不同的方式对其显示。在如下图示中，可以看出这幅图像基本上只有两种种类。8.9、讨论与分析：在上面endmember collection 对话框中，apply后，可以修改其最大波谱角，从而可以是分类有所改变。第九章 地物

21、识别8.1、报告人：郭美辰8.2、实验目的：了解地物识别的过程，增强解译图像的能力。8.3、实验原理：8.4、实验时间：8.5、实验数据： 一幅多波段的遥感影像8.6、使用软件： Envi8.7、步骤：在菜单栏波谱沙漏工具spectralspectral hourglass wizard 下，出现一个有详细说明的对话框，点击next,出现如下对话框。在select input file菜单下，选择原图。在select output root name菜单下，选择要输出的路径，点击next，再点击next，会生产一幅最小特征分离的图。点击 load MNF Result to ENVI Disp

22、lay ，生成一幅图像，如下图所示。点击load animation of MNF bands，会生产一种以动画的方式表现的波段信息。点击next，选择其中的yes，next。设置参数，其中反复迭代的指数可以小点，阈值根据需求而定，内存也可以自己定，大点就快些。即可生成一幅像素纯净指数的汇总图。再next，出现如下图幅。点击右下方的start，可以观看类别的选取情况，如果不满意，可以自己调节。如果满意，就点击retrieve endmembers。当点击 plot endmembers ,会出现一幅特征曲线的光谱图。点击next，出现对话框。选择load SAM class result,生成

23、一幅自动识别地物类型的图像。其中选择MTMF，可以对混合像元进行分类，知道其一个像元所占的百分比。地物识别方法二：步骤：打开一幅多光谱图像，选择一种想要识别的地物，然后右键Z profile(spectrum)显示其波谱曲线。在spectral菜单下，SAM target finder with BandMax中，说明了一个操作的流程，点击next。在select input file 下，选择要识别的地物所在的图像。并且选择输出路径。在导入的对话框中，可以选择import from plot windows,选中，next。Next，设置背景值，在next，设置其最大光谱角，load SAM rule image，完成。9.8、实验结果：9.9、讨论与分析：该方法适用于基于波谱库的目标地物的分类，当对植被类型不是很了解的时候，这样分类会很精确，比如树的具体类型等等。