ENVI实习指南

ENVI实习指导 《 遥感应用分析原理与 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 》课程 ENVI实习指南 TOC \o "1-3" \h \z 1. 界面系统介绍 2 2. 文件的存取与显示 2 3. 图像预处理 3 4. 影像分析 5 影像统计分析 5 图像增强 (Image:Enhance) 5 专题信息提取 5 5. 定义感兴趣区(ROI)及分类 6 监督分类（Supervised Classification） 6 非监督分类（Unsupervised Classification） 9 6. 分类后处理（Post Classification） 9 7. 制图输出 10 图形的整饰 10 图形输出 11 8. 图像存取的举例说明 11 ENVI实习指南 ENVI（The Environment for Visualizing Images）是一套功能齐全的遥感图像处理系统，是处理、分析并显示多光谱数据、高光谱数据和雷达数据的高级工具，由美国系统研究公司（Research System INC.）开发。 IDL（Interactive Data Language）交互式数据语言是进行二维及多维数据可视化分析及应用开发的理想软件工具。 1. 界面系统介绍 1.​ 主菜单： 菜单项，File、Basic Tool、Classification、Tranform、Spectral实习所涉及的（粗略介绍） 2.​ Help 工具的使用 3.​ 主菜单设置（preferences）：内存设置 2. 文件的存取与显示 1．​ 图像显示 由一组三个不同的图像 窗口组成：主图像窗口、滚动窗口、缩放窗口。 1）​ 主图像Image窗口：（400*400） 100％显示（全分辨率显示）scroll的方框，可交互式分析、查询信息。主图像窗口内的功能菜单：在主图像窗口内点击鼠标右键，切换隐藏子菜单的开启和关闭。该 "Functions" 菜单控制所有的ENVI交互显示功能，这包括：图像链接和动态覆盖；空间和波谱剖面图；对比度拉伸；彩色制图；诸如ROI的限定、光标位置和值、散点图和 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 面图等交互特征；诸如注记、网格、图像等值线和矢量层等的覆盖（叠置）；动画以及显示特征。 2）​ 滚动Scroll窗口：全局，重采样(降低分辨率)显示一幅图像。只有要显示的图像比主图像窗口能显示的图象大时，才会出现滚动窗口。滚动窗口位置和大小最初在 envi.cfg 文件中被设置并且可以被修改。 3）​ 缩放Zoom窗口：（200*200）显示image的方框。缩放系数（用户自定义）出现在窗口标题栏的括号中。Scroll、Zoom可通过Image右键>Display Window style 重新打开。 2．​ 图像的头文件资料的获取和编辑 ENVI：File>>Edit ENVI Header，选择相应的文件。 从 Header Info 对话框里，你可以点击 Edit Attributes 下拉菜单中的选项，调用编辑特定文件头参数的独立对话框。这些参数包括波段名、波长、地图信息等。 3．​ 图像的存取 File > Open Image File. 当你打开任何文件，可用波段列表（ABL）自动地出现(或ENVI:Window>> ABL)。ABL列出该图像文件的所有波段，并允许你显示灰阶和彩色图像、启动新的显示窗口、打开新文件、关闭文件，以及设置显示边框。 要选择当前活动显示，请按以下步骤： 从ABL（Available Bands List）内，点击 “Display #X” 按钮菜单（其中 “X” 是与显示窗口标题栏内数字相对应的数字），再从列表中选择所需要的显示。 要开始一个新的显示，从按钮菜单选择 “New Display”。 点击 “Load Band” 或 “Load RGB” ，以把选定的波段导入选定的显示。 4．​ 灰度图像和彩色图像的显示ENVI：File>>Open Image>>Available Bands List（ABL）中选择Gray Scale或RGB Color模式 5．剖面和波谱图（Profiles and Spectral Plots） Image:>>Tools>>Profiles。ENVI 允许抽取水平的（X）、垂直的（Y）、波谱的（对每个像元为 Z ）以及任意的剖面图。剖面图显示在单独的图表窗口，并且X、Y 和 Z 剖面图可以同时是激活的。鼠标用来移动一个十字准线并交互地选择剖面图。图表窗口内 Options 菜单下的 Auto Scale Y-Axis非常有用。 3. 图像预处理 1．​ 图像的切割（取子区） ENVI：Basic Tools>> Resize Data >> Resize Data Input File对话框（如下图）。 ①选择需要切割的原始图像；②选择Spatial Subset或Spectral Subset方式；③若设置空间切割方式（Spatial Subset>>select Spatial Subset）点击 “Image”（另有‘map-输入图左上角、右下角坐标’，‘file-输入文件名’）；④出现Subset by Image对话框，Subset的尺寸用2种形式，移动图像上的方框或直接填写samples/lines（列/行）值；③’若设置波段范围（Spectral Subset>>File Spectral Subset），选择波段； 若要根据已选择的感兴趣区域进行切割，可用ENVI：Basic Tools>>Subset Data via ROIs。若要使用与上次输入的空间大小相同的文件的空间子集，点击 “Previous” 按钮。 注：Resize Data还可以进行图像重采样（如下），若仅仅进行子区的选择，则不要调整Output File Dimensions。* 图像左上角为原点(1.1 --- 列.行)。 “Masking”---把一个空间掩膜应用到图像的某个部分，包括统计、分类、分离（unmixing）、匹配滤波、连续删除（continuum removal）和波谱特征拟合（spectral feature fitting）。 2．​ 图像的重采样 ENVI：Basic Tools>> Resize Data >> Resize Data Input File对话框选择需要采样的原始图像—OK >>Resize Data Parameters——调整 Output File Dimensions的像元数；选择采用方法 >> 文件 输出 4. 影像分析 影像统计分析 1.​ 统计特征分析 ENVI：Basic Tools>> Statistics >> Compute Statistics >>Calculate Statistics Parameters对话框---选择统计类型（Basic、Histogram、Covariance）--- 输出统计文件.sta（file）---出现Statistics Results对话框Select Plot菜单可查看统计结果（Min/Max/Mean/ Stdev /Eigenvalue最小值、最大值、平均值、 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 差、特征值，Covariance协方差矩阵，Correlation相关矩阵，Eigenvector特征向量，Histogram 波段直方图） 2.​ 主成分分析 ENVI: Transform>>principle compents>>Forward PC Rotation > Compute New Statistics and Rotate.---- 选择输入文件，出现‘Forward PC Parameters’对话框，输出统计文件.sta（memory）---‘Select Subset from Eigenvalues’（yes）---出现‘Select Output PC…’显示PC各分量的特征值、方差百分比，选择输出PC波段数（4）---‘OK’---出现‘PC EigenValues’窗口（每个节点是PC各分量的特征值）--- 可进一步计算各波段信息量大小（通过PC1的方差百分比及PC1中各波段协方差矩阵的特征向量）。 3.​ 相关分析 ENVI：Basic Tools>> Statistics >> Compute Statistics 图像增强 (Image:Enhance) 1.​ 直方图调整 (1)​ 直方图匹配 Image: Enhance>>Histogram Matching 至少显示两幅图像，从想更改直方图的图像（如“Display#1—PC1”）中，选择Enhance>>Histogram Matching---‘Match To’选择想匹配直方图的图像“Display#2--V” ---“OK”，保存直方图匹配后的PC1’。 查看两图像（PC1’与V）直方图：点右键Interactive Stretching或 选择 Functions > Interactive Stretching显示直方图；若需‘图像替代’则要求两直方图输入值相同，—可根据两直方图输入值的关系，通过‘Band Math’使两直方图数值相同（PC1’变为PC1’’）---保存PC1’’，可为下一步PC1’’图像替代V，进行HSV-RGB反变换作准备。(如直方图输入值V为0.1-0.9，PC1’为0-255，则通过‘Band Math’b1*0.91/255，可得PC1’’) (2)​ 直方图的交互式拉伸 Image: Enhance>>Interactive Stretching ENVI 用 2% 的系统默认线性拉伸值来显示所有图像(两边均舍去信息量的2%)，经过这样处理后合成的假彩色图像层次分明、地物差异大，各类地物易于判别。 注：多在Scroll窗口对全局调整，分别调整R、G、B，使彩色更丰富（一般R、G敏感，B线性即可）。练习：分段线性拉伸（ctrl+左键）、均衡化Equalization、Gaussian等。 2.​ 彩色变换 ENVI: Transform>>Color Transform 包括‘HSV-色调、饱和度、数值’变换，‘HLS--色调、亮度、饱和度’变换等。 练习：彩色增强（RGB---HLS---调整‘S—S’---RGB）；数据融合（RGB---HLS---H+S+ PC1’’ ---RGB）。 3.​ MNF变换 最小噪声分离（Minimum Noise Fraction，MNF）变换是同主分量变换相似的一种方法，它被用来分离数据中的噪声，确定数据内在的维数，减少随后处理的计算量（Green 等人，1988；Boardman和Kruse，1994）。 Forward --- 选图 --- 保存噪声图MNF’、去噪声图MNF --- 保存MNF.sta --- 输出（memory）；Inverse --- 选MNF图象 --- 选MNT.sta --- 输出已去噪声的还原图像（与原图像相比，仅min/max/mean变化，其余统计值均不变）。 专题信息提取 1.​ NDVI的提取 ENVI: Transform>> NDVI(vegetation Index) NDVI（Normalized Difference Vegetation Index）归一化植被指数，是一个普遍应用的植被指数，将多波谱数据变换成显示植被分布的图像波段。NDVI 值指示着像元中绿色植被的数量，较高的NDVI值预示着较多的绿色植被。NDVI 变换可以用于 AVHRR、Landsat MSS、Landsat TM、SPOT 或 AVIRIS 数据，也可以输入其他数据类型的波段来使用。 2.​ 波段运算获取不同专题信息 ENVI: Basic Tools>>band math 例：Newband=band5-band4 具体操作是： 打开Band Math对话框(如右图)，在Enter an expression中键入：b5-b4，点击OK后将会出现Variables to Bands Pairings对话框。从可利用波段列表中，分别选择b5和b4代表的波段，并键入待输出的文件名，点击OK即可。 练习：湿度指数提取： 1.0*[(b2-b4)/(b2+b4)] 或 1.0*[(b2-b5)/(b2+b5)] 3.​ 缨帽变换 ENVI: Transform>>Tasseled Cap 4.​ 2D显示 Image：Tool>>2D Scatter Plots>>Scatter Plot Band Choice对话框---选择X、Y轴 --- OK ---出现2D 散点图‘Scatter Plot’。 5. 定义感兴趣区(ROI)及分类 监督分类（Supervised Classification） 监督分类：按照分类以前自定义的样本进行分类。 1.​ 训练样本的选择和优化 1）训练样本的提取（ROI区的选择） ENVI: Basic Tools >> Region Of Interest>> ROI tool调出感兴趣区工具窗口进行样本选择（注意：必须事先打开一幅图像），可以进行样本编辑（名称，颜色，填充方式等），样本选择越精确，分类结果越好。感兴趣区工具窗口的打开方式还有：Image：Overlay>>Region of Interest，或者直接在图像窗口上点击鼠标右键，再选择ROI Tool。 根据前面的背景资料和预处理结果选择视觉效果最佳的彩色合成图像，建立各类地物的训练区。各类地物的解译标志，即地物明显的影像特征---色调、纹理等，通过目视解译方法用鼠标在工作区影像图上选择其训练区，并使训练区的分布尽量均匀。在实际的工作中，由于存在“同物异谱”的情况，因此对于同一种类型可能有多种不同的特征。为此，我们可以对同一地物选择多个训练区，分类后再合并。 感兴趣工具窗口见下图。 提取训练样本的具体操作如下： （1）确定ROI的提取类型（ROI Tool: ROI_Type >>Polygon, Polyline, Point, Rectangel, Ellipse）和待操作窗口（主图像窗口、滚动窗口或缩放窗口）。 （2）在图像窗口（Image 或Zoom）上用左键画出感兴趣区，单击鼠标右键确定选择形状（此时可以拖动感兴趣区域，用Ctrl＋鼠标左键可以删除），再次单击右键确定此训练区（此时若要删除训练区，需要点击ROI Tool窗口中的Delete控键，此操作将删除所有该类型的感兴趣区域）。ROI Tool 窗口中将会显示选择区域的颜色和相关信息，其中，感兴趣区域名称（ROI Name）和色彩可以修改。可就某一类训练区选择多个感兴趣区域。 （3）该类训练区的选择完成后，点击ROI Tool窗口的New Region控键，再进行另一类训练样本的选择，其颜色将自动改变。按以上操作完成所有训练区的选择。 2）训练样本的优化和提纯ROI 上述步骤中选择的某类训练样本，可能混入了其他类型的样本，为了提高图像分类精度，需要对训练样本进行提纯。N维可视化分析器（N — Dimensional Visualizer）即是对选择的训练区像元进行提纯。当某些像元始终聚集在一起运动时，这些就是所需的最纯像元；若在运动时，像元分成了两个部分，则说明选择了两类地物的训练区，需把此训练区像元分开处理。 ROI Tool: File>>Export ROIs to n-D Visualizer >> n-D Control; n-D Visualizer 让训练区像元在n维空间内自由转动（可以控制转动速度Speed），当转到最能区分各类型训练区像元的位置时，停止转动，进行样本提纯操作。即： （1）在n-D Visualizer窗口中用鼠标选择某类训练区 --- 先点击右键选new class --- 在选择纯化的训练区内用左键画出纯像元区，未确认状态下可用‘ctrol+左键’删除，点击右键确定（可进行多次选择） ---再次单击右键>>Export Class，提纯后的训练区将出现在ROI Tool窗口中，确定后的改动，需到ROI Tool窗口，用delete。 （2）进行下一个类型训练区的提纯时，首先要在n-D Visualizer窗口中单击右键>>New Class，下面的操作如前。如此，完成所有训练区的提纯。 ROI Tools窗口：① 对优化后的原类别样本‘Hide’看效果，② 可删除旧类别样本，以优化后的替代，③ 编辑色彩，以达到较好视觉效果。 （3）训练区的保存：ROI Tool: File>>Save ROIs… 2.​ 选择分类方式 分类方式包括平行六面体法、最短距离法、马氏距离法、最大似然法、波谱角分类以及二进制编码法等，选择合适的分类方式。 1) 最大似然法（Maximum Likelihood Classification） ENVI: Classification >> supervised >>Maximum likelihood >> Classification Input File 选择分类的图像 >> Maximum likelihood Parameters 选择训练样本，设置说明 最大似然分类假定每个波段每一类统计呈均匀分布，并计算给定像元属于一特定类别的可能性。除非选择一个可能性阈值，所有像元都将参与分类。每一个像元被归到可能性最大的那一类里。 在Maximum Likelihood Parameters 对话框中设置一般分类参数，在 “Set Probability Threshold” 文本框里，键入一个阈值（0~1）。选项参数被用来控制像元准确分类的可能性。如果像元的可能性低于所有类的阈值，则它被归为“无类别”，在此，我们一般选择默认值。 2)​ 波谱角分类法（Spectral Angle Mapper — SAM） ENVI: Classification >> supervised >> Spectral Angle Mapper >> Classification Input File 选择分类的图像 >> Spectral Angle Mapper Parameters 选择训练样本，设置说明 波谱角分类法是以物理学为基础的一种分类法，通过比较终端光谱向量和每个像元的矢量在N维空间中的角度，将像元分配到相应的区间中去，角度值越小，分类越精确。 输入由上步提纯得到的像元数据，选择适宜的参数[Maximum Angle（radians）] 值，小于此值的像元将不参加分类，经多次实验。默认值是0.1（弧度）。 3.​ 分类 引入影像—>确定分类范围和波段—>选择训练样本—>给定阈值—>确定存储路径和文件名—OK。下图为最大似然法分类对话框。 非监督分类（Unsupervised Classification） 非监督分类：仅仅用统计方法对数据集中的像元进行分类，不需要样本。 方法： （1）IsoData: ENVI: Classification>>unsupervised>>IsoData>> IsoData Parameters对话框：参数设置说明 在ISODATA Parameters 对话框中，输入Number of Classes（分类数），Min（最少分类数）8、Max（最大分类数）15，Maximum Iteration（最大迭代数）10，Chang Threshold（像元变化的阀值）5.00，Minimum ＃Pixel in Class（每类中的最小像元数）1，Maximum Class Stdv （最大标准差）3.00，Minimum Class Distance （最小类间距）4.00，Maximum ＃Merge Pairs（最大合并数）2等8个基本参数（根据实际图像和先验知识更改参数的设置） （2）K-Means: ENVI: Classification>>unsupervised>> K-Means>> K-Means Parameters参数设置说明 6. 分类后处理（Post Classification） 1.​ 分类统计：ENVI:Classification>>Post Classification>>Class Statistics：包括每一类的点数、最小值、最大值、平均值以及类的每个波段的标准差等。其中每一类的最小值、最大值、平均值以及标准差可以以图的方式进行显示。可以显示出每一类的直方图，并且计算其协方差矩阵、相关矩阵、特征值和特征矢量等。 2.​ 类别集群：ENVI: Classification>>Post Classification>>Clump Classes细小块的合并，将一些碎块进行合并（平滑处理）。 注：未被选上用于聚块（clumping）的类，在输出图像上无变化。 3.​ 类别筛选：ENVI: Classification>>Post Classification>> Sieve Classes通过用斑点分组消除这些隔离的被分类的像元。该功能菜单将删除分类中的孤岛像元，并用黑像元表示，可以用成块分类功能代替黑像元。 注：在 “Group Min Threshold” 文本框里，输入一个类组需要包含的最少像元数（4或8）。任何一组小于这一数值的像元将从类中被删除。 4.​ 类的合并：ENVI : Classification>>Post Classification>>Combine Classes：将分过的类进行选择性的合并，可以合并为一类或几类。 5.​ 类的叠合：ENVI : Classification>>Post Classification>>Overlay Class：用一幅彩色合成影像或灰阶影像生成一幅影像地图，并且类的颜色叠置在一起，输出一幅3波段的 RGB图像。 6.​ 修改类的颜色：Image: Tools>>Color Mapping>>class color mapping：当一个分类后的图像被导入一个显示窗口时，每类自动呈现出不同的颜色。每类的颜色与监督分类中选择的感兴趣区的颜色或非监督分类中预先选择的每类颜色相对应，未分类区域在图像中呈黑色。此处可进行各类别的颜色、名称等修改。（？） 7.​ 分类精度评价：ENVI:Classification>>Post Classification>>Confusion Matrix： 通过分类结果图与地表真实感兴趣区（Ground Truth ROIs）相比较来计算混淆矩阵。评价结果 记录 混凝土 养护记录下载土方回填监理旁站记录免费下载集备记录下载集备记录下载集备记录下载 了总体精度、准确度、Kappa系数、混淆矩阵、commission 误差（每类中额外像元占的百分比）和冗长误差（类左边的像元占的百分比）等等。当用地表真实图像（Ground Truth Image）计算混淆矩阵时，还可以输出每类图像中没有被正确分类的那些像元。 7.​ 制图输出 图形的整饰 1.​ 图面大小（边界）设置Image: overlay>>annotation 出现 Annotation 对话框 --- Option>>Set Display Border --- 出现‘Display Border’对话框，输入左（80）、右（180）、上（150）、下（120）的图像边框宽度，设置边框颜色（白）---‘OK’，滚动窗口（Scroll）出现带白色边框的图。 2.​ 经纬网格线 Image: overlay>>Grid lines（经纬网和像元网）,若图像没有被配准到地图坐标，将不会显示出经纬网。可以通过：Grid Line Parameters:Options>>…来设置网格线的属性及图像边界。 3.​ 注记Image: overlay>>annotation(标题；图例；比例尺；南北指针) (1)​  添加注记：Image: overlay>>annotation。注记可以被放置在主图像窗口、滚动窗口或缩放窗口。通过从各自的 Options 菜单中选择 Annotation，每种图表，包括 X、Y、Z 剖面图或表面图，可以被注记。当出现 #n Annotation 对话框时（其中 “#n” 指正被注记的那个显示），选择 Object > >所需要的注记对象。 ​ 文本注记（Text）：Object>>text ，选择注记的属性（如font、size、color等），在文本框中输入待添加文本，用鼠标在图像中点击注记位置，按右键确定。 ​ 图例注记（map key）：Object>>map key可以直接将各类的图例加载上去， Image：Tools—color mapping-class color mapping可进行颜色、名称等编辑。 Object>>Map Key>>Edit Map Key Items可改变图例顺序：通过利用Add Item和Delete Item来调整顺序。 ​ 比例尺注记（Scale Bar）：Object>>Scale bar ​ 偏差图注记（Declination）：Object>>Declination，Declination 选项允许你在图像上放置一个磁偏角图表。磁偏角图表包括指向真北（True North，用星号显示）、坐标北（Grid North，GN）和磁北（Magnetic North，MN）的箭头的任意组合。 ​ 覆盖分类结果：Image:Overlay>>classification ​ 其他注记略 ​ 注记的选择或编辑：Object>>Selection/Edit 圈定 ---出现‘◇’待编辑状态 注：每一种注记添加到图像中后，单击鼠标右键进行确定。所有注记对象都有一个小的彩色菱形 “handle”，它被用来决定位置。 (2)​  保存注记文件： Annotation: File > Save Annotation. 出现对话框输出文件名，或点击 “Choose”，然后当出现 ENVI 文件选择对话框时后，输入一个文件名。 ​ 图形输出 1．分类结果图像输出：Image: File>>Save Image As >>Image File （or Postscript File）‘Output Display’--- Type：Jepg --- 输出名 2．分类结果的矢量输出：ENVI: Classification >> Post Classification >> Classification to Vector：将选择的类变换为矢量多边形进行输出。 注：输出图像名不能与其他文件名相同。 8.​ 图像存取的举例说明 ENVI中支持多种格式的图像文件，有些格式无法直接打开，需要添加头文件信息。下面以一DAT格式为例，说明如何打开其他格式的图像文件，并转换为ENVI标准格式。 如下图是.dat格式的TM图像文件，每一个文件包含一个波段信息，header为头文件信息。 （1）File >> Open Image File，选择待打开的图像，如先打开band1，将会出现头文件对话框；从头文件中查出图像行列数（若行列数填写错误将不能正确打开图像），填入Header Info对话框中；波段数选择1（注意这里的波段数是指待打开的文件包含几个波段，并不是指第几波段）。 注：若无法用File >> Open Image File则需要用File >> Open External File来选择相应的文件格式类型。Landsat-geotiff-全选—ABL（单波段） （2）点击OK后band1文件将出现在可利用波段列表ABL中。重复以上操作打开所有的波段。若要将这些波段合并到一个ENVI标准格式文件中，需进行下面的操作。 （3）单波段文件合并ENVI: File >> Save File As>>ENVI Standard ，将出现New File Builder对话框，点击Import File…，选择待添加的波段文件，确定后文件将出现在对话框中；点击Reorder Files…可改变文件顺序ctrl点击或拖；最后键入输出文件名，新文件将出现在可利用波段列表ABL中。 （4）修改文件中的波段名ABL表中 图像的重采样 不同像元大小波段的文件不能够合并到同一个文件下，此时需要重新采样。ETM中61、62 波段像元数为4428×3910，8波段像元数为17710×15640，其余波段像元数为8855×7820，重采样使像元数相同。 各波段文件合为一个标准文件