ENVI试验参考指导书专业资料

遥感图像解决实验指引实验一图像解决软件ENVI4功能简介实验二影像配准实验三图像融合实验四波段组共计算及图像增强实验五图像分类实验六地图制图实验一图像解决软件ENVI4功能简介1ENVI4支持文献格式ENVI使用是通用栅格数据格式，包括一种简朴二进制文献和一种有关ASCII头文献。读该文献格式容许ENVI使用几乎有影像文献，涉及那些自身嵌入头文献信息影像文献。通用栅格数据都会存储为二进制字节流，普通它将以BSQ（按波段顺序）、BIP（波段按像元交叉）或者BIL（波段按行交叉）方式进行存储。ENVI软件支持数据类型，涉及字节型，整型，无符号整型，长整型，无符号长整型，浮点型，双精度浮点型，复数型，双精度复数型，64位整型，以及无符号64位整型。2ENVI窗口和显示使用ENVI软件时，打开图像普通有三个窗口，即主影像窗口（MainImageWindow）、滚动窗口（ScrollWindow）和缩放窗口（ZoomWindow），见图1.1。图1.1ENVI显示窗口组：主影像窗口、滚动窗口和缩放窗口2.1ENVI主菜单ENVI主菜单位于显示窗口顶部，在默认状况下，它显示于主影像窗口顶部，并提供交互式影像显示和分析功能，见图1.2。如果所选显示组中不涉及主影像窗口，那么菜单将会出当前滚动窗口或者缩放窗口顶部，同ENVI其她菜单同样，从该菜单中可以选取任意菜单选项。此外，在任何一种显示窗口点击右键会弹出一种快捷菜单。通过这个菜单可以访问显示函数，也可以变化显示设立。图1.2主影像窗口中overlay菜单和快捷菜单2.2可用波段列表ENVI可以打开影像文献或者这些文献中单个波段。可用波段列表（AvailableBandsList）是一种特殊ENVI对话框，它涉及了因此被打开文献中可用影像波段，以及与此有关地图信息列表，见图1.3。使用波段列表可以把彩色和灰阶影像加载到一种显示窗口中，可以直接打开一种新显示窗口，也可以点击对话框底部Display按钮，并从中选取显示窗标语，打开相应显示窗口，然后点击GrayScale或者RGB单选按钮，再从列表中点击波段名，选取所需要波段，加载显示影像。提示：如果加载是单波段影像，可以直接在波段上双击来完毕。图1.3可用波段列表对话框3ENVI基本功能打开影像文献一方面启动ENVI软件，双击软件图标即可。要打开一种影像文献可以分如下几种环节：1选取File→OpenImagefile。在屏幕上弹出文献选取对话框“EnterInputDataFile”。2选取进入实验课程所提供数据目录envidata中，双击打开文献夹，选取can_tm文献夹，从文献夹中选取can_tmr.img文献，然后点击ok。3在可用波段列表中就会显示刚才选取文献can_tmr.img，使用鼠标左键点击选取影像显示方式（GrayScale或RGB）。选取某个影像波段或三个波段进行显示，所选波段在标有“SelectedBand”显示出来。4点击LoadBand或LoadRGB，将影像加载到一种新显示窗口。图像加载后，可以熟悉一下窗口各个界面，以及对图像主窗口、滚动窗口和缩放窗口进行有关操作。显示光标位置要显示鼠标光标位置和值，可以从主影像窗口菜单栏或者从ENVI主菜单栏中选取Window→CursorLocation/Value，或者在主影像窗口中点击鼠标右键，从弹出快捷菜单中选取CursorLocation/Value。接着在屏幕上浮现CursorLocation/Value对话框将显示出光标在主影像窗口。滚动窗口或者缩放窗口中位置，见图1.4。该对话框还显示了十字丝光标所相应哪个象素屏幕颜色值和实际数据值。要关闭这个对话框，可以在CursorLocation/Value对话框顶部菜单中，选取file→cancel。一旦CursorLocation/Value对话框打开后来，要隐藏或显示该对话框，可以在主影像窗口中双击鼠标左键即可。图1.4CursorLocation/Value对话框显示影像剖面廓线可以交互式选取和显示X轴（水平）、Y轴（竖直）和Z轴（波谱）剖面廓线图。这些剖面廓线图显示了穿过影像横线X，纵线Y或者波谱波段Z数据值。1从主影像窗口菜单栏中，选取Tools→Profiles→Xprofle，将会打开一种绘制窗口，该窗口将依照影像中所选取行，绘制出一幅数值与列号（samplenumber）之间关系曲线图，见图1.5。2重复上述过程，选取Yprofile，绘制出一幅数据值与行号（linenumber）之间关系曲线图。选取Zprofile，绘制出相应波谱面廓线图。提示：也可以通过任意影像窗口中快捷菜单，来打开Z轴波谱剖面廓线图。3选取Window→MouseButtonDescriptions，可以查看鼠标键在剖面廓线显示窗口中按键功能描述。4放置好这三个剖面廓线窗口，以便能同步看到它们。5在影像上移动移动十字丝，查看这三幅影像剖面廓线图在新位置上是如何显示新数据值。6从绘制图每个窗口中，选取File→Cancel，来关闭绘制窗口。图1.5水平轴（X）、Y轴（竖直）和Z轴（波谱）剖面廓线绘制曲线图进行迅速对比拉伸使用主影像窗口、缩放窗口或者滚动窗口中默认参数和数据来进行迅速对比度拉伸。选取主影像窗口菜单栏中Enhance菜单，可以进行各种各样对比度拉伸（线性拉伸，0－255间线性拉伸，2％线性拉伸、高斯拉伸、均衡化拉伸以及平方根拉伸）。1使用主影像窗口、缩放窗口或者滚动窗口中影像作为拉伸数据源，尝试进行各种类型拉伸变换。2比较在窗口显示组中线性、高斯、均衡及平方根拉伸后成果。显示交互式散点图ENVI可以绘制出两个所选影像数值关系图，即分别选定这两个波段为X、Y轴，在平面坐标上绘制两者散点图。1在主影像窗口菜单栏中，选取Tools→2DScatterPlots。接着ScatterPlotBandChoice对话框就会出当前屏幕上，在该对话框中选取要进行比较两个影像波段。2选取其中一种波段作为X轴，另一种波段作为Y轴，然后点击OK，即可生成散点图。3一旦打开了散点图绘制窗口（图1.6）就可以将鼠标光标放在主影像窗口中任意位置，并可以按住鼠标左键来拖动光标，此时，十字丝光标周边10×10范畴内像素在散点图中所相应点将会用红色突出显示出来。提示：选取Window→MouseButtonDescriptions来显示不同鼠标键在ScatterPlot显示窗口中功能。4在主影像窗口上移动鼠标光标，观测所产生跳跃像素（dancingpixels）效果。5可以使用散点图来突出显示主影像窗口中具有相应波谱值像素，将鼠标光标放在散点图窗口上，点击并拖动鼠标中间键。6在散点图菜单栏中，选取file→cancel来关闭ScatterPlot窗口。图1.6一种交互式散点图，对波段2和波段3进行比较加载一幅彩色图像1如果屏幕上没有显示可用波段列表对话框，那么在ENVI主菜单栏中，选取Window→AvailableBandslist，打开该对话框，如图1.3。注意：如果没有打开影像，那么在ENVI主菜单中，选取File→OpenImageFile，然后选取一种要打开影像。接着该影像波段就加载到可用波段列表中。2点击RGB单选按钮，然后在另一种显示窗口中，加载一幅彩色影像。3通过点击列表中波段名，从列表中为每种颜色（红、绿、蓝）选取一种波段。当点击列表中波段名时，指定RGB单选按钮会自动向前选取。4当三种颜色均有波段与其相应时，点击Display＃1下拉式菜单按钮，并从中选取NewDisplay。5点击LoadRGB按钮，就可以将一幅彩色影像加载到一种新显示窗口中。链接两个显示窗口将两个显示窗口链接在一起进行比较。当把两个显示窗口链接在一起后，在一种显示窗口中（滚动、缩放等）所进行任何操作，都会在与其相链接显示窗口中产生相似响应。要在两个显示窗口中链接在一起，操作环节如下：1从主影像窗口菜单中，选取Tools→link→LinkDisplays，或者在影像中点击鼠标右键，在弹出快捷菜单中，选取LinkDisplays。LinkDisplays对话框就会打开。2在LinkDisplays对话框中，点击ok，建立链接。3当前尝试在其中一种窗口进行滚动或者缩放操作，观测另一种窗口反映。选取感兴趣区域（RegionofInterest）ENVI容许在影像中定义感兴趣区（regionofinterest，ROIs）。感兴趣区重要被用于提取分类信息记录、生成掩膜以及其她某些操作。1从主影像窗口中选取Overlay→Regionofinterest，或者在影像中点击鼠标右键，在弹出菜单对话框中，选取ROITool。接着与该影像显示窗口相相应ROITool对话框就会出当前屏幕上（图1.7）。2绘制感兴趣区多边形，在ROI_Type中选取ROI_Type类型，涉及多边形、矩形、圆形等。详细绘制环节如下：在主影像窗口中点击鼠标左键，建立感兴趣多边形第一种顶点。再依次点击鼠标左键，按顺序选取更多边线点，点击鼠标右键来闭合该多边形。鼠标中键可以删除最新定义点，或者删除整个多边形。再一次点击鼠标右键，固定多边形位置。通过选取ROITool对话框顶部相应单选按钮，感兴趣区也可以在缩放窗口中被定义。当完毕一种感兴趣区定义后，该区域会在AvailableRegionsofInterest列表中显示出来，同步显示尚有感兴趣区名称、颜色和所包括像素数。3要定义另一新感兴趣区，在ROITool对话框中点击NewRegion即可定义。效果见图1.7。图1.7定义感兴趣区对话框图1.8定义了三个感兴趣区图像4、保存ROIROITOOL窗口菜单：File/SaveROIs…5、合并ROIROITOOL窗口菜单：Options/MergeRegions6、裁剪ROIROITOOL窗口菜单：File/SubsetDataviaROIs…或：主菜单，BasicTools/SubsetDataviaROIs…对影像进行注记在ENVI软件中，可以向地图和影像中添加某些文本、符号、色标条及某些其她符号注记。1要对一幅图像进行注记，可以从主影像菜单选项中选取Overlay→Annotation。接着与主影像窗口相相应Annotation：Text对话框就会出当前屏幕上（见图1.9）图1.9文本（Text）模式Annotation对话框2要对绘制图、3D表面以及相似对象进行注记，可以从绘制窗口菜单栏中选取Options→Annotation。注记类型：Annotation：Text对话框中容许添加不同类型注记。这些不同类型都可以从object下拉菜单中进行选取，注记类型重要有:文本、符号、矩形、椭圆、多边形、折线、箭头、地图比例尺、三北方向图标（DeclinationDiagrams）、地图图例、颜色标注及影像注记。默认状况下，注记类型为文本（Text）。对话框中其她选项用来控制文本注记字体（Font）、大小（Size）、颜色（color）、位置和角度。放置注记：例如向主影像窗口中添加文本注记：1在Annotation：Text对话框中，添加要输入文本（forexample，CenterSouthUniversityofForestry&Technology）在对话框相相应菜单和参数设立中选取文本字体、颜色和大小，然后在主影像窗口适当位置上，点击鼠标左键。接着，输入文本会显示在所选点位置上，如图1.10。图1.10带注记影像注意：选取Window→MouseButtonDescription，描述在注记对话框中鼠标按键各项功能。2使用鼠标左键，拖动文本注记小圆柄，在窗口中放置文本注记。可以在对话框中变化相应区域设立值，或者按下鼠标左键拖动文本或符号注记，以此来变化注记属性和位置。3完毕文本注记设立，可以点击鼠标右键来锁定注记位置。保存和恢复注记：1从Annotation：Text对话框菜单栏中，选取File→SaveAnnotation，来保存影像注记。2打开一种OutputAnnotationFilename对话框。要保存影像注记，在该对话框中指定要保存途径以及保存文献名，该文献扩展名为.ann。注意：如果没有将影像注记保存到文献中，那么关闭Annotation：Text对话框时，注记也将丢失。3在Annotation：Text对话框中，选取File→Restoreannotation，就可以恢复保存过注记文献。修改先前放置注记要对先前已经放置好注记进行修改，可按如下环节进行修改：1在Annotation：Text对话框中，选取Object→Selection/Edit。2用鼠标点击并托拽出一种矩形框来包括要修改注记对象。3当小圆柄浮现后，点击拖动注记及小圆柄，修改其相应属性，这就像设立新注记对象时同样。添加网格1要在影像中添加公里网格信息，可以从主影像窗口中选取Overlay→GridlineParameters对话框。注意：当给影像添加公里网格时，影像边框会自动添加进来。2在GridlineParameters对话框，设立公里网线宽、颜色和公里网间隔，来修改公里网属性。图1.11添加注记和网格影像3在EditPixelAttributes对话框中，可以变化公里网标注、网格线、矩形边框和交叉角颜色、宽度以及公里网间隔。完毕了这些属性设定之后，点击EditPixelAttributes对话框中ok按钮，将所做更改应用到这些影像当中。4当对所加公里网满意之后，点击GridlineParameters对话框中Apply按钮。保存和输出影像对添加标注和进行图像解决图像进行保存，可以将影像保存为ENVI影像文献格式，或者保存为几种通用图像格式，然后打印或者导入到其她软件当中，也可以直接通过打印机进行打印输出。将影像保存为ENVI影像格式1从主影像窗口菜单栏种，选取File→SaveImageFile，OutputDisplayToimagefile对话框就会出当前屏幕上。2选取24比特彩色或者8比特灰阶进行输出，然后再选取其她图形选项（涉及注记和公里网）以及边框设立。如果将注记和公里网都添加到彩色影像显示中了，那么注记和公里网就会自动列到graphics选项中，也可以选取其她注记文献应用于输出图像上。3使用所需单选按钮，将影像成果输出到Memory或者file中。如果选取了file，就要输入一种输出文献名。4点击ok保存影像。实验二影像配准本专项旨在简介如何在ENVI中对影像进行地理校正，添加地理坐标，以及如何使用ENVI进行影像到影像几何校正。1图像文献头文献修改打开并显示SPOT数据一方面打开软件，打开软件自带数据库中SPOT数据，环节如下：1从ENVI主菜单中，选取file→openimagefile。2当enterdatafilename文献选取对话框浮现后，选取进入envidata目录下bldr_sp.img文献。3点击ok。4当可用波段列表对话框浮现后，点击Greyscale单选按钮，使用鼠标左键，点击相应波段名，从对话框顶部所列波段中选中SPOT波段。所选取波段名显示在selectBand：字段区域中。5点击LoadBand按钮，加载这幅影像到一种新显示窗口中。修改ENVI头文献中地图信息1在可用波段列表中，右键点击bldr_sp.img文献名下Mapinfo图标，从弹出快捷菜单中选取EditMapInformation。EditMapInformation对话框出当前屏幕上。图2.1这个对话框列出了在ENVI中添加地理坐标所用地理信息。可以调节ENVI使用MagicPixel（作为地图坐标系统起始点）相相应影像坐标。由于ENVI可以从相应头文献信息和地图投影文献中，辨认出地图投影，像元大小以及地图投影参数，因此用它能计算出影像中任意像元地理坐标，既可以输入地图坐标，也可以输入地理坐标（经纬度）。2点击Projection/Datum文本旁边箭头切换按钮，显示UTMZone13North地图投影纬度/经度坐标。ENVI在解决过程中才进行转换。3点击当前DMS或者DDEG按钮，分别在度分秒（Degree－Minutes－Second）和十进制度（DecimalDegrees）之间进行切换。4点击Cancel，推出EditMapInformation对话框。图2.1EditMapInformation对话框5修改图像pixelsize信息，添加公里网格和地图标注。6保存图像。在主影像窗口中，选取file→saveimageas→imagefile。选取输出途径和填写输出文献名称。2影像对影像几何配准运用SPOT图像校正LandsatTM。打开并显示LandsatTM图像1从ENVI主菜单中，选取file→openimagefile。2当enterdatafilenames对话框浮现后，选取进入envidata目录下bldr_reg子目录，从列表中选取bldr_tm.img文献。3在文献选取对话框中，点击open，把TM影像波段加载到可用波段列表中。4在列表中选取band3，点击display#1按钮，并从下拉式菜单中选取newdisplay。5点击LoadBand按钮，把TMband3波段影像加载到一种新显示窗口中。显示光标位置/值要打开一种显示主影像窗口，滚动窗口，或者缩放窗口中光标位置信息对话框，可以按如下环节进行操作。1从主影像窗口菜单栏中，选取tools→CursorLocation/value。2在主影像窗口、滚动窗口和缩放窗口TM影像上，移动鼠标光标。注意坐标是以像素为单位给出，这是由于这个影像是基于像素坐标，它不同于上面带有地理坐标SPOT影像。3选取file→Cancel，关闭CursorLocation/value对话框。开始进行影像配准1从ENVI主菜单栏中，选取Map→Registration→SelectGCPs：ImagetoImage。2在imagetoImageRegistration对话框中，点击并选取display＃1（spot影像），作为BaseImage。点击display＃2（TM影像）作为Warpimage。图2.2imagetoImageRegistration对话框3点击ok，启动自动配准程序。通过讲光标放置在两幅影像相似地物点上，来添加单独地面控制点。注意：同名地物点，需要在两张图上均有，普通选取河流和道路交叉点或者拐角点。交叉点或拐角点选取要根据图像空间辨别率大小来拟定，如Quickbird空间辨别率为0.61米，则对Quickbird进行图像校正时，道路交叉点宽度应在1米以内为好。4启动GroundcontrolPointsSelection对话框后，一方面在两张图上找到同名地物点，如当前SPOT图像上找到一点，将鼠标移动至该点，同样办法在TM影像上找到该点。5在两个缩放窗口（Zoom）中，查看光标所在位置，把十字丝移到相应位置，记录每点X,Y值。图2.3缩放窗口（Zoom）对话框6GroundControlPointsSelection对话框baseX和Y文本框中，分别输入286和255，将SPOT影像中光标位置移动到相应位置点上。使用同样办法在WarpX和Y文本框中，分别输入133和264，将TM影像中光标移动到相应位置。7GroundControlPointsSelection对话框中，点击Addpoint，把该地面控制点添加到列表中。点击Showlist查看地面控制点列表，见图2.4。然后依同样方式寻找别的同名地物点。图2.4用来进行影像配准GroundControlPointsSelection对话框和Showlist列表注意：在缩放窗口中支持亚像元（sub-pixel）级定位，缩放比例越大，地面控制点精度就越好。地面控制点选取除了通过移动光标之外还可以直接在同名地物点上通过鼠标进行点击完毕。在showlist对话框中，一旦已经选取了4个以上地面控制点后，RMS误差就会显示出来。8如果选取控制点中，某点误差很大，删除该点，重新寻找新点来代替。在showlist对话框中用鼠标点击该点，选取delete即可删除该点。如果对所有选取控制都不满意话，可以通过GroundControlPointsSelection对话框，选取Options→ClearAllPoints，可以清除掉所有已选取地面控制点。9若选取地面控制点比较满意，则在对图像进行校正之前要对控制点信息进行保存。详细办法为在GroundControlPointsSelection对话框中，选取file→SaveGCPstoASCII…，选取保存途径和文献名，对控制点信息进行保存，以备下次使用。图2.5影像配准中所用到ImagetoImageGCPList对话框校正影像在这里咱们校正影像为TM影像一种波段，咱们也可以同步校正多波段影像中所有波段。1从GroundControlPointsSelection对话框中，选取Options→WarpDisplayedBand。2在RegistrationParameters对话框中WarpMethod按钮菜单中，选取RST。在Resampling按钮菜单中选取NearestNeighbor重采样法。图2.6RegistrationParameters对话框3输入文献名bldr_tm1.wrp,点击ok。4重复环节1和环节2，还是使用RST校正法，但是要相应选取Bilinear和CubicConvolution重采样法。5将成果分别输出为bldr_tm2.wrp和bldr_tm3.wrp文献中。6再一次重复环节1和环节2，这一次选取一次多项式Polynomial校正法，并使用Cubicconvolution重采样法。然后再选取Delaunay三角网Triangulation校正法，相应使用Cubicconvolution重采样法。7将上述成果分别输出到bldr_tm4.wrp和bldr_tm5.wrp文献中。按照同样办法可以对TM影像别的波段进行校正，与上述环节有区别地方是，对别的波段进行校正，不需要重新找点，只需要将保存控制点信息直接调用即可。在GroundControlPointsSelection对话框中，选取file→RestoreGCPsfromASCII，调用保存控制点信息即可对图像进行校正，校正后图像进行合成。8合成环节为，在ENVI主菜单栏中选取basictools→laystacking，通过laystacking选取校正后各个波段进行图层叠加，对叠加后图像进行保存，保存文献名为bldrtm_m.img。比较校正成果使用动态链接来比较校正成果：1在可用波段列表中，点击原始SPOT波段影像名bldr_sp.img，然后从菜单中，选取Greyscale显示方式来显示图像。2在可用波段列表中，选取bldr_tm1.wrp文献，在Display＃下拉式按钮中选取NewDisplay，点击LoadBand将该文献加载到一种新显示框中。3在主影像窗口中，点击鼠标右键，选取tools→Link→linkdisplays。4在linkdisplays对话框中，点击ok，把SPOT影像和已添加了地理坐标bldr_tm1.wrp影像链接起来。5使用动态链接功能，与带地理坐标SPOT影像进行比较。6取消动态链接功能，选取tools→Link→unlinkdisplays。查看地图坐标1从主影像窗口菜单栏中，选取Tools→cursorLocation/Value。2浏览带地理坐标数据集，注意不同采样办法和校正办法对数据值所产生效果。3选取File→Cancel，关闭该对话框。3影像到地图配准这一某些将逐渐地演示影像到地图配准解决过程。许多环节同影像到影像配准环节相似，因而这些环节将不会被详细地讨论。从带地理坐标SPOT影像中获取地图坐标以及一种矢量数字线划图（DLG）都将被作为基准数据，然后对基于像素坐标LandsatTM影像进行校正，以匹配相应地图数据3、1打开并显示LandsatTM影像文献：从ENVI主菜单中，选取File→OpenImageFile。当EnterDataFilenames对话框浮现后，选取进入相应文献夹目录下列表中lxmbldr_tm.img文献。点击OK。TM影像波段被加载到可用波段列表中，同步一幅彩色影像被加载到一种新显示窗口中。在可用波段列表中，点击GrayScale按钮，选取波段3。点击LoadBand按钮，把TM影像第3波段加载到已打开显示窗口中。3、2选取影像到地图配准并恢复控制点坐标:从ENVI主菜单中，选取Map>Registration>SelectGCPs:ImagetoMap。如果打开了各种影像显示窗口，那么就在ImagetoMapRegistration对话框中，点击选取包括该灰阶影像那个显示窗口显示号。从投影列表中选取UTM，并在Zone文本框中输入13。设立像素大小为30m，点击OK，启动配准程序。在要校正影像中，把光标移动到一种已知地图坐标地面点上（可以从一幅地图或者ENVI矢量文献中[见下一某些]读取所需地图坐标），来添加单个地面控制点。GroundControlPointsSelection对话框中E（东向）和N（北向）文本框中，手动地输入已知地图坐标，然后点击AddPoint来添加新地面控制点。在GroundControlPointsSelection对话框中，选取File>RestoreGCPsfromASCII，打开lxmPZhun.pts文献。在GroundControlPointsSelection对话框中，点击ShowList按钮。可以在ImagetoMapGCPList对话框中查看影像地图坐标、实际影像点和预测点坐标以及RMS误差。3、3使用矢量显示数字线划图（DLGs）来添加地图控制点：从ENVI主菜单中，选取File>OpenVectorFile>USGSDLG。在文献选取对话框中，选取lxmbldr_rd.dlg文献。在ImportOptionalDLGFileParameters对话框中，选取Memory单选按钮，点击OK，读入所需数字线划图（DLG）数据。在可用矢量列表中高亮选取ROADSANDTRAILS：BOULDER，CO文献，点击LoadSelected按钮。在LoadVector对话框中，点击NewVectorWindow。把该矢量加载到一种新矢量显示窗口中。在VectorWindow#1窗口中，点击并拖曳鼠标左键，激活一种十字形光标。光标处地图坐标会在VectorWindow#1窗口底部列出。在主影像显示窗口中，选取Tools>PixelLocator，并输入402和418，然后点击Apply，将影像光标移动到道路交叉口相应点上去。注意，在缩放窗口中，同样可以获取到亚像元（sub-pixel）级定位精度。带十字形光标矢量窗口，并且显示出了地图坐标(上图)在矢量窗口中，用鼠标左键点击并拖曳矢量光标，当十字形光标位于所需道路交叉口时，松开鼠标左键，把矢量光标放置在道路交叉口上，其坐标为477593.74，4433240.0（北纬40度3分3秒，东经-105度15分45秒）。在矢量窗口中，点击鼠标右键，并从弹出快捷菜单中，选取ExportMapLocation。新地图坐标就会出当前GroundControlPointsSelection对话框中。在GroundControlPointsSelection对话框中，点击AddPoint，添加该地图坐标/影像像素对，并观测RMS误差变化。3、4RST和三次卷积校正：GroundControlPointsSelection对话框中，选取Options>WarpFile。在InputWarpImage对话框中，高亮选取文献名zhumuqian_tm.img，点击OK，对TM6个波段进行校正。在浮现RegistrationParametersdialog对话框中，将WarpMethod选为RST，将Resampling设立为CubicConvolution。把background值改为255。在outputfile文本框中，输入输出文献名zhumuqiantm_m.img。点击OK，开始进行影像到地图校正。3、5显示成果并进行评价：使用光标位置/值（CursorLocation/Value）来对校正后彩色影像进行评价。在可用波段列表中，点击GrayScale单选按钮，接着点击校正影像波段Wrap(Band3：zhumuqiantm_m.img)。从Display#下拉式菜单按钮中，选取NewDisplay。点击LoadBand，来加载这幅校正后TM影像。注意到校正影像是倾斜，这是由于消除了LandsatTM轨道方向影响因素。此时这个影像已经带有地理坐标，但注意它空间辨别率是30米，而SPOT影像为10米辨别率。实验三图像融合数据融合是将多幅影像组合到单一合成影像解决过程。它普通使用高空间辨别率全色波段与低空间辨别率多光谱影像进行融合，达到增强光谱影像空间辨别率问题。进行数据融合影像文献必要具有地理坐标，在这种状况下，空间采样会自动运营。如果进行数据融合影像文献没有地理坐标，则影像就必要覆盖同一地理区域，并且有相似像素大小，影像大小以及相似方位。1无地理坐标系统影像融合读取并显示ERMapper影像LondonTM影像和SPOT影像都是二进制文献，具有ERMapper头文献。它们可以自动被ENVIERMapper读取程序读取。1选取file→openexternalfile→IPsoftware→ERMapper，然后进入lontmsp子目录，选取lon_tm.ers。2在可用波段列表中，单击RGBcolor单选按钮，并依次选用红、绿、蓝字段相应波段，点击LoadRGB来显示一幅彩色TM影像。3选取file→openexternalfile→IPsoftware→ERMapper，然后进入lontmsp子目录，选取lon_spot.ers文献。4在可用波段列表中，单击greyscale单选按钮，选用pseudolayer波段，从display下拉菜单中选取newdisplay，点击loadband来显示一幅灰阶SPOT影像。调节影像大小1在可用波段列表中单击SPOT影像可以发现其空间尺寸为2820×1569，用同样办法可以懂得TM影像空间尺寸为1007×560。TM影像像素大小为28米，而SPOT影像像素大小为10米。TM影像必要以2.8倍率来调节大小，以产生与SPOT影像匹配10米大小像元。2选取Basictools→ResizeData（spatial/Spectral）选取lon_tm.ers并单击ok。在resizedataparameters对话框xfac文本框中输入2.8，在yfac文本框中输入2.8009，选取保存途径和文献名lontm，点击ok来调节TM影像大小。或者直接在samples文本框中输入2820，在lines文本框中输入1569，选取保存途径和文献名lon_tm，点击ok来调节TM影像大小。图3.1ResizeData对话框3显示调节过影像，选取Tools→link→LinkDisplays将调节过大小TM影像和SPOT全色影像链接起来，使用动态链接来分析比较这两幅影像。图3.2TM彩色合成影像和SPOT全色影像进行HIS数据融合为了理解整个解决过程，选取进行手工数据融合。第一步，TM彩色影像转换到色度、饱和度和数值(也称亮度值)（hue-saturation-value(intensity)）彩色空间。将高辨别率SPOT影像替代数值（value）波段，并将其拉伸到0至1之间以满足对的数据范畴，再将从TM影像获取色度、饱和度以及SPOT影像中获取数值进行反转换，转回到红-绿-蓝彩色空间。这个过程将产生出一幅输出影像，其中包括了从TM影像中获取颜色信息以及从SPOT影像中获取空间辨别率信息。HSV正变换1从ENVI主菜单选取Transform→colorTransforms→RGBtoHSV，然后选取调节过大小TM数据作为输入RGB影像。输入要输出文献名，点击ok执行变换。2作为灰阶影像或RGB彩色影像，来显示色度、饱和度和数值影像。拉伸SPOT影像并替代TM数值波段1从ENVI主菜单选取BasicTools→StretchData，单击lon_spot文献，然后点击ok。2在Datastretching对话框OutputData中，在min文本框中输入0，Max文本框中输入1.0，并输入一种输出文献名。单击ok将SPOT影像数据拉伸为浮点型，范畴为0到1.0。图3.3Datastretching对话框HSV反变换1从ENVI主菜单中选取Transform→ColorTransforms→HSVtoRGB，选取通过变换TM影像Hue和Saturation波段作为变换H和S波段。2选取拉伸过SPOT影像作为V波段，点击ok。在HSVtoRGBParameters对话框中输入要输出文献名，点击ok进行反变换。图3.4HSVtoRGBParameters对话框显示成果1在可用波段列表中点击RGBColor单选按钮，并按顺序点击反变换后RGB波段。再点击LoadRGB按钮来显示一幅通过融合TM/SPOT彩色影像。2显示融合后影像，选取Tools→Link→LinkDisplays将融合后影像，同调节过大小TM影像以及SPOT全色影像链接起来。使用动态分析比较这些影像。进行自动HSV变换融合1在ENVI主菜单选取Transform→ImageSharpening→HSV。2如果调节过大小TM彩色影像已经显示在窗口中，则可以在SelectInputRGB对话框中直接选取适当影像显示窗口。否则，就要在SelectInputRGB对话框中，选取“RedLayer”，“GreenLayer”，“Bluelayer”所相应调节过大小TM影像波段，然后点击ok。3HighResolutionInputFile对话框中选取SPOT影像，点击ok。4输入输出文献名：lontmsp.img，在HSVSharpeningParameters对话框中点击ok。成果显示、链接和比较1显示融合后彩色影像，在可用波段列表中，选取RGBColor单选按钮，再在RGB波段中，选取融合影像中相应波段，点击LoadRGB。2通过在主影像窗口菜单中选取Tools→Link→LinkDisplays，将原始TM彩色影像，SPOT影像以及自动融合彩色影像进行链接比较。图3.5TM彩色影像、SPOT影像以及融合后彩色影像链接比较2带地理坐标系统影像融合打开要进行融合图像1选取file→openimagefile然后进入envidata目录，双击p123r043ETM文献夹，选取p123r043_7p19991224_z49_multi.hdr，点击打开，影像就出当前可用波段例表中。2在可用波段列表中，点击RGBColor，选取R（4）G（5）B（3）三个波段来进行波段显示。3选取file→openimagefile然后进入envidata目录，双击p123r043ETM文献夹，选取p123r043_7p19991224_z49_pan.hdr，点击打开，影像就出当前可用波段例表中。4在可用波段列表中，点击Greyscale，选取p123r043_7p19991224_z49_pan.hdr，点击Display＃1下拉菜单，点击Newdisplay，点击Loadband。ColorNormalized（Brovey）sharpening（原则化彩色变换融合）1在ENVI主菜单选取Transform→ImageSharpening→ColorNormalized（Brovey）。2在selectInputforcolorBands对话框中选取AvailableBandsList，点击ok。图3.6selectInputforcolorBands对话框3在随后浮现可用波段列表对话框中，从p123r043_7p19991224_z49_multi.hdr文献中选取三个波段相应RGB选项，点击ok。4在highresolutionInputfile对话框中，选取p123r043_7p19991224_z49_pan.hdr文献波段，点击ok。图3.7highresolutionInputfile对话框5在ColorNormalizedsharpening对话框中，Resampling重采样办法选取NearestNeighbor，输入输出文献名和保存途径，点击ok，就会进行自动融合，融合成果就会出当前可用波段例表当中。图3.8ColorNormalizedsharpening对话框6查看波段融合成果，在可用波段列表对话框中，选取RGBcolor显示方式，选取R（4）G（5）B（3），点击LoadRGB，显示融合图像。Gram-schmidtspectralsharpening（GS变换融合）1在ENVI主菜单选取Transform→ImageSharpening→Gram-schmidtspectralsharpening。2在selectlowspatialresolutionmultibandinputfile对话框中选取p123r043_7p19991224_z49_multi文献，点击ok。4在highresolutionInputfile对话框中，选取p123r043_7p19991224_z49_pan文献波段，点击ok。5在Gram-schmidtspectralsharpening对话框中，Resampling重采样办法选取NearestNeighbor，输入输出文献名和保存途径，点击ok，就会进行自动融合，融合成果就会出当前可用波段例表当中。6查看波段融合成果，在可用波段列表对话框中，选取RGBcolor显示方式，选取R（4）G（5）B（3），点击LoadRGB，显示融合图像。图3.9selectlowspatialresolutionmultibandinputfile对话框PCspectralsharpening(主成分变换融合)1在ENVI主菜单选取Transform→ImageSharpening→PCspectralsharpening。2在selectlowspatialresolutionmultibandinputfile对话框中选取p123r043_7p19991224_z49_multi文献，点击ok。4在highresolutionInputfile对话框中，选取p123r043_7p19991224_z49_pan文献波段，点击ok。5在PCspectralsharpening对话框中，Resampling重采样办法选取NearestNeighbor，输入输出文献名和保存途径，点击ok，就会进行自动融合，融合成果就会出当前可用波段例表当中。6查看波段融合成果，在可用波段列表对话框中，选取RGBcolor显示方式，选取R（4）G（5）B（3），点击LoadRGB，显示融合图像。注意：变化图像融合重采样办法，看看与否会对图像显示产生影响。从图像清晰度和色彩两个方面来比较不同图像融合办法，对图像影响。实验四波段组共计算及图像增强在本专项中，以ETM数据为对象，简介在图像解决过程中波段组合方式，波段之间运算方式，数据拉伸及增强解决过程。打开影像1选取file→openimagefile然后进入envidata目录，双击p123r043ETM文献夹，选取p123r043_7p19991224_z49_multi.hdr，点击打开，影像就出当前可用波段例表中。2在可用波段列表中，点击RGBColor，选取R（4）G（5）B（3）三个波段来进行波段显示。3点击LoadRGB，一幅假彩色图像就显示在影像窗口中。波段组合1在显示影像中，只用了p123r043_7p19991224_z49_multi.hdr文献三个波段，而文献有6个波段，在这里，可以尝试从6个波段当中，选用3个波段来组合，并比较不同波段组合之间图像显示效果。2在可用波段列表中，点击RGBColor，选取R（7）G（4）B（1）三个波段来进行波段显示。图4.1不同波段组合间影像显示效果3在可用波段列表中，点击Display＃1，在下拉菜单中，选取Newdisplay，点击LoadRGB，一幅真彩色图像就显示在影像窗口中。4动态链接比较先后两幅影像效果，在主影像窗口中，选取tools→link→linkdisplays，在linkdisplays对话框中，点击ok，两幅影像就链接起来了。5在可用波段列表中，选取别的波段组合方式，并加以比较。波段运算Bandmath功能为顾客提供了一种灵活图像解决工具，在Bandmath对话框中，可以实现不同波段之间加减乘除等运算。在这里还是以p123r043_7p19991224_z49_multi.hdr文献为对象进行波段运算。1在主菜单栏中，选取file→openimagefile然后进入envidata目录，双击p123r043ETM文献夹，选取p123r043_7p19991224_z49_multi.hdr，点击打开，影像就出当前可用波段例表中。2在主菜单栏中，选取BasicTools→Bandmath。将浮现bandmath对话框。图4.2Bandmath对话框3在enteranexpression文本框中，输入需要进行波段计算IDL（InteractiveDataLanguage）表达式，使用变量代替波段名或文献名，变量名必要以字符“b”或“B”开头，背面跟着5个以内数字字符。注意：例如，若想计算三个波段平均值，可以将数学表达式：（float(b1)+float(b2)+float(b3)）/3输入到文本框中，该表达式中使用三个变量b1表达第一种变量，b2表达第二个变量，b3表达第三个变量，在变量前用浮点型字节来防止计算时浮现字节溢出错误，输入有效表达式之后，点击Addtolist，表达式将添加在对话框中，点击ok。4变量赋值，在variablestoBandsPairings对话框中，容许从一种输入波段列表中，把波段赋给“enteranexpression”文本框中包括变量。图4.3variablestoBandsPairings对话框要把一种值给原先实例中变量“b1”；在variablestoBandsPairings对话框中，在标签为“Variablesusedinexpression”文本框内，点击变量“B1”。在标签为“AvailableBandsList”列表中，点击所需要波段，数据集随后显示在变量名之后，注意一旦第一种波段被选中，只有那些具备相似空间尺寸波段被显示在波段列表中，按照同样办法，为“B2”，“B3”等赋予一种值。5输出波段运算成果，三个变量都选定了相应波段之后，点击ok对运算成果进行输出，在VariablestoBandsPairings对话框中，选取“file”或“Memory”切换按钮，运算成果将显示在可用波段列表中。提示：除了波段均值计算之外，还可以尝试波段之间比值、加减或者乘法等一系列运算。波段运算普通满足四个条件：1）必要使用IDL语言书写波段运算表达式；2）输入波段必要具备相似维数；3）表达式中所有变量都必要用Bn（bn）命名；4）成果波段必要与输入波段维数相似。交互式数据拉伸（InteractiveStretching）打开要进行拉伸数据p123r043_7p19991224_z49_multi.hdr，在主影像窗口，选取Enhance→interactivestretching，在交互式对比度拉伸对话框中将显示一种输入直方图何一种输出直方图，它们表白了当前输入数据以及分别应用拉伸。两条垂线（虚线）表白了当前拉伸所用到最小值和最大值。对于彩色影像来说，直方图颜色与所选波段颜色一致（系统默认显示红色波段直方图）。交互式直方图窗口底部会列出拉伸类型和直方图来源。要浏览一幅彩色图绿波段和蓝波段直方图，点击窗口中“G”、“B”切换按钮。要把任何拉伸或直方图变换自动应用到图像，选取Options→AutoApply。普通状况下只有点击“Apply”按钮才干把任何变换应用到图像，通过取消选中Options→AutoApply来关闭自动功能。图4.4交互式对比度拉伸对话框设立拉伸参数在交互式对比度拉伸对话框中，使用“Options”菜单和鼠标交互功能可觉得拉伸设立参数，并同步浏览交互式直方图窗口中信息。更改拉伸最小值和最大值更改拉伸最小值和最大值，在任意一条垂直虚线上点击鼠标左键，然后拖放到一种新位置，或者直接在对话框顶部文本框中输入所需值。提示：当在文本框中指定最大和最小值时，输入所需要实际值（例如，最小值37，最大值68）或者数据比例（例如，最小值2％，最大值98％）。可以设立数据范畴之外值，但是只依照输入直方图数据范畴进行拉伸。输入数据后，输出直方图会自动更新，以反映输入直方图中变化，并应用到新拉伸显示数据分布。锁定拉伸条（lockingstretchbars）要锁定最小值和最大值拉伸条（垂直虚线）间距离，选取Options→lockingstretchbars。拉伸条间距离将被锁定，可以同步移动两个拉伸条。要解除锁定，取消选中lockingstretchbars即可。拉伸类型菜单在交互式对比度拉伸对话框中Stretch_type菜单中包括一种所有可供选取交互式拉伸类型列表。线性对比度拉伸线性对比度拉伸是系统默认交互式拉伸。线性拉伸最小值和最大值分别设立为0和255，两者之间所有其他值设立为中间线性输出值。在交互式对比度拉伸对话框中，选取Stretch_type→linearcontraststretch。要限定最小和最大输入值，使用鼠标左键，移动输入直方图中垂直线（白色虚线）到所需要位置，或者在stretch文本框中输入所需要DN值或一种数据比例（例如，最小值5％，最大值95％）。要重新显示本来拉伸，选取Options→ResetStretch。分段线性对比度拉伸（piecewiselinear）分段线性对比度拉伸可以通过使用鼠标在输入直方图中放置几种点进行交互地限定。对于各个点之间某些采用线性拉伸。选取Stretch_type→piecewiselinear。一种转换函数（初始为一条白色直线）将被绘制在输入直方图。在输入直方图地任何位置点击鼠标中键，从而为转换函数增长一种节点，绘制线段将把端点和绘制节点标记连接起来。图4.5分段线性对比度拉伸对话框要移动一种点位置，在标记上点击鼠标左键，然后把它拖放到一种新位置。要删除点，在标记上点击鼠标右键。要手动地键入输入和输出值，选取Options→Editpiecewiselinear。成果将被绘制在输出直方图中，该直方图显示应用新拉伸后数据分布状况。点击Apply，将拉伸应用于显示数据。3）高斯对比度拉伸（Gaussian）系统默认Gaussian拉伸使用均值DN127和相应于0-255以正负3为原则差值进行拉伸。选取Stretch_type→Gaussian。输入拉伸最小值和最大值。输出直方图用一条红色曲线显示被选取Gaussian变换函数。被拉伸数据分布呈白色，并叠加显示在红色Gaussian函数上。注意：要手动地输入所需要原则差，选取Options→SetGaussianStdv，点击Apply，把拉伸应用于显示数据。4）直方图均衡化对比度拉伸要自动缩放数据，从而使每个直方图bin中DN数相均衡；选取Stretch_type→Equalization。输入直方图显示未被修改数据分布。输出直方图用一条红色曲线显示均衡化函数，被拉伸数据分布呈白色叠加显示。点击Apply，将拉伸应用于显示数据。5）平方根对比度拉伸（SquareRoot）计算输入直方图平方根，然后应用线性拉伸，选取Stretch_type→SquareRoot，输入直方图显示未被修改数据分布。输出直方图用一条红色曲线显示平方根函数，被拉伸数据分布呈白色叠加显示，点击Apply，将拉伸应用于显示数据。6）自定义对比度拉伸和直方图匹配（Arbitrary）Arbitrary选型容许在输出直方图顶部绘制任何形状直方图，或与另一种图像直方图相匹配。选取Stretch_type→Arbitrary。输入直方图未被修改数据分布，通过点击或按住并拖放鼠标左键，可以在OutputHistogram窗口绘制输出直方图，自定义直方图将用绿色来显示。点击Apply，将拉伸应用于显示数据。注意：要删除直方图，点击鼠标中键。点击鼠标右键来接受绘制输出直方图，并把数据记录资料匹配到图中。输出直方图用红色显示自己绘制直方图，匹配数据函数用白色曲线绘制。7）直方图匹配使用自定义对比度拉伸功能也可以把一幅图像直方图与另一幅图像直方图进行匹配。要从一幅图像中获得输入或输出直方图，在图表顶部“InputHistogram”或“OutputHistogram”文本标签上点击鼠标左键拖放即可。把名字拖放到其他自定义输出直方图中，然后释放按钮。被导入直方图将用红色绘制，输出直方图将被拉伸，与导入直方图匹配。点击Apply，将拉伸应用于显示数据。8）定义查找表拉伸（DefiningLook-Up-Table）一种顾客自定义查找表可以把每个输入DN值拉伸到一种输出值。它可以恢复或交互式定义。选取Stretch_type→UserDefinedLUT。选取Options→EditUserDefinedLUT。当浮现编辑对话框时，一种包括输入DN值和相应拉伸输出值列表将显示在“EditUserDefinedLUT”标签下，这些值反映了当前拉伸状况。在值上点击进行编辑。当它出当前“EditSelectitem:”文本框中时，输入所需值，然后按回车键。注意：当该功能被启动后，要重新设立查找表为初始值，点击“Reset”。点击ok，点击Apply，将拉伸应用于显示数据。图4.6EditUserDefinedLUT对话框图像增强与滤波图像增强功能在主影像窗口中enhance下拉菜单中，有线性增强、高斯增强、直方图均衡增强和平方根对比增强几种类型，每种增强类型有三个选取，可以对主影像窗口、卷积窗口和缩放窗口分别进行增强解决。滤波普通通过特定空间频率来使图像增强。空间频率普通描述亮度或DN值与距离方差，图像涉及各种不同空间频率，消除一幅图像高频信息可以使图像显得平滑。在ENVI软件中进行图像滤波，是在主菜单栏中filter下进行操作，重要涉及卷积滤波（ConvolutionsandMorphology）、形态学滤波、纹理滤波、自适应滤波和频率滤波。卷积滤波是指在空间域对图像进行滤波。形态学滤波以形态为基本对图像进行解决。纹理滤波从图像中提取纹理信息。自适应滤波器在抑制噪声同步保存了图像尖锐信息和细节。傅立叶滤波在频率域对图像进行滤波。卷积滤波在主菜单中，选取filter→ConvolutionsandMorphology。将浮现ConvolutionsandMorphologyTool对话框。图4.7ConvolutionsandMorphologyTool对话框选取滤波器类型，在ConvolutionsandMorphologyTool对话框中，点击Convolutions，在下拉菜单中可以选取滤波器类型（高通、低通、拉普拉斯、直通、高斯高通、高斯低通、中值、sobel、Roberts、顾客自定义滤波），系统默认滤波器类型为高通类型。编辑变换核大小，在ConvolutionsandMorphologyTool对话框中，用鼠标点击“KernelSize”增减箭头时，变换核尺寸将以2为单位增减；当鼠标中键点击“KernelSize”增减箭头时，变换核尺寸将以10为单位增减；当用鼠标右键点击“KernelSize”增减箭头时，变换核尺寸将恢复到默认尺寸。变换核被默认设立为正方形。要将变换核变为非正方形核，选取Option→SquareKernel：No。默认ImageAddBack为0。纹理滤波器（基于概率记录滤波）使用Texture选项可以应用于基于概率记录或二阶概率记录纹理滤波。许多图像包括区域以亮度变化为特性，而不但局限于亮度值。纹理是指图像色调作为级别函数在空间上变化。被定义为纹理清晰区域，灰度级别相对于不同纹理地区一定是比较接近。在主菜单中，选取filter→Texture→occurrenceMeasures。在Textureinputfile对话框中，选取要进行滤波文献，点击ok。图4.8Textureinputfile对话框在OccurrenceTextureParameter对话框中，选取输出文献途径核文献名，别的选项采用默认值，点击ok，文献即进行纹理滤波解决。所选纹理图像将计算出来，并被放置在可用波段例表中。图4.9OccurrenceTextureParameter对话框自适应增强型Lee滤波器自适应滤波运用环绕每个像元原则差来计算一种新像元值。普通来说，原始像元值将被基于周边有效像元（即那些符合原则差像元）所计算新值代替。增强型Lee滤波器是Lee滤波器改进，依照单独滤波窗口中计算出记录（方差系数）对数据进行滤波。每个像元将被分到3个类型中：相似像元（homogenous）、差别像元（heterogeneous）和指向目的像元（pointtarget）。在主菜单中，选取filter→Adaptive→enhanceLee。在EnhanceLeefilterInputfile对话框中，选取要进行滤波文献，点击ok，在EnhanceLeefilterParameters对话框中，输入filterSize大小、DampingFactor及HomogenousAreas等参数，也可以采用默认方式。参数填写完毕，选取输出途径和输出文献名，点击ok即可。图4.10EnhanceLeefilterParameters对话框4）迅速傅立叶变换滤波器FFTFiltering（FastFourierTransformFiltering迅速傅立叶变换滤波）可以将图像变换成为显示不同空间频率成分合成输出图像。涉及如下操作：频率域滤波器交互式建立；滤波器应用；以及FFT向原始数据空间逆变换。正向FFT滤波打开影像，在ENVI主菜单栏中选取Filters→FFTFiltering→ForwardFFT。将浮现ForwardFFTInputFile对话框，选取要进行滤波文献，点击ok。在ForwardFFTParameters对话框中选取输出到“file”或“memory”。点击ok开始FFT计算。定义FFT滤波器使用FilterDefinition选项可以交互式直接定义滤波器，也可以通过在显示正向变换图像中绘制来定义滤波器。滤波类型涉及：circlepassandcut以及user-drawnpassandcut。在ENVI主菜单栏中选取Filters→FFTFiltering→FilterDefinition。将浮现FilterDefinition选取对话框。从下列选项中选取：如果当前显示中包括一幅FFT影像，选取包括FFT图像显示号，选取“NoDisplay”指出该滤波器不与某个特定图像有关。如果当前显示中不含FFT图像，则在“Samples”和“Lines”文本框中键入数值，以限定滤波器大小。图4.11FilterDefinition对话框在FilterDefinition选取对话框中，选取Filter_type即滤波器类型，对于“Circularpass”或“Circularcut”滤波器，需要在“Radius”文本框中，以像元为单位输入滤波半径。对于“BandPass”或“BandCut”滤波器，在“InnerRadius”和“OuterRadius”文本框中，以像元为单位键入所需值。对于“UserDefined”和“UserDefinedCut”滤波器，容许将ENVI注记导入滤波器（普通在FFT影像上，使用annotation功能）。使用箭头按钮来选取“NumberofBorderPixels”参数，用于细化滤波器。值为0表达没有平滑。选取输出到“file”或“memory”。点击apply构建FFT滤波器。反向FFT变换反向FFT变换程序实际涉及两个环节，先应用FFT滤波，然后将FFT图像反变换回原始数据空间。选取Filter→FFTFiltering→InverseFFT，将浮现InverseFFTInputFile对话框。选取一种输入文献或波段，若需要，选用任意子集。点击ok即可对图像进行反变换。实验五图像分类本专项简介遥感影像分类操作解决，比较监督分类和非监督分类后影像，并对分类后影像进行相应后解决。后解决涉及聚合（clump）解决，筛选（sieve）解决、并类（combine）解决，以及精度评估。打开显示分类影像1从ENVI主菜单中，选取File→OpenImagefile，选取进入envidata文献夹目录，点击can_tm子目录，从列表中选取can_tm.img文献，点击ok。2在可用波段列表中，选取RGBColor单选按钮，然后用鼠标左键，顺次点击波段4、波段5和波段3。所选波段就会在对话框中恰当文本框中显示出来。3点击LoadRGB按钮，把该影像加载到一种新显示窗口中。在主影像窗口对影像进行简朴区别和目视判读工作。非监督分类在ENVI主菜单中，选取Classification→Unsupervised→选取K-means或者IsoData。K均值（K-means）分类法K－均值非监督分类器使用了聚类分析办法，它需要一方面选定所需分类个数，随机地查找聚类中心位置，然后迭代地重新配备它们，直到达到最优化波谱分类。1在ENVI主菜单中，选取Classification→Unsupervised→K-means，在ClassificationInputfile对话框中，选取分类影像和分类波段数，点击ok。2在K-meansparameters对话框中，拟定分类数为6类，别的参数采用默认方式，选取文献输入名称和输出途径，点击ok，影像进行分类。3在可用波段列表中，添加一种新图像显示窗口，从Display＃1下拉菜单中，选取NewDisplay。4双击K-means（can_tm.img）分类图像，分类后图像就显示在新影像窗口中。5从主影像窗口菜单中，选取Tools→Link→LinkDisplays，然后在对话框中，点击ok。图5.1K-meansparameters对话框6使用鼠标左键，在影像上点击并拖动动态叠加显示区域，将K-means分类成果同原始彩色合成影像进行比较。7当解决完毕后，选取Tools→Link→UnlinkDisplay，关闭动态链接。IsoData（迭代自组织数据分类技术）IsoData非监督分类法将计算数据空间中均匀分布类均值，然后用最小距离规则将剩余像元进行迭代聚合。每次迭代都要重新计算均值，且依照所得新均值，对像元进行再分类。这一解决过程将持续到每一类像元数变化少于所选像元变化阈值或者达到了迭代最大次数。1在ENVI主菜单中，选取Classification→Unsupervised→IsoData，在ClassificationInputfile对话框中，选取分类影像和分类波段数，点击ok。2在IsoDataparameters对话框中，拟定分类数为6类，别的参数采用默认方式，选取文献输入名称和输出途径，点击ok，影像进行分类。3在可用波段列表中，添加一种新图像显示窗口，从Display＃1下拉菜单中，选取NewDisplay。4双击IsoData（can_tm.img）分类图像，分类后图像就显示在新影像窗口中。5从主影像窗口菜单中，选取Tools→Link→LinkDisplays，然后在对话框中，点击ok。6使用鼠标左键，在影像上点击并拖动动态叠加显示区域，将IsoData分类成果同原始彩色合成影像进行比较。按住鼠标左键同步，再次点击鼠标右键，触发第三幅影像动态叠加显示。将IsoData分类成果和K-means分类成果进行比较。7当解决完毕后，选取Tools→Link→UnlinkDisplay，关闭动态链接。监督分类法监督分类是一种边训练边学习分类办法。在进行监督分类时，一方面要训练分类样本，拟定分类数，对样本进行保存。然后选取监督分类办法，分类结束后进行分类后解决工作。1选取用于非监督分类影像can_tm.img文献。从主影像窗口中，选取Overlay→RegionofInterest。接着相应于显示窗口ROITool对话框。图5.2ROITool对话框2在ROITool对话框中，选取ROI_type为polygon，然后在主影像窗口中，绘制一种多边形，该多边形就代表了新创立感兴趣区。要完毕这一步，请按下面环节进行：在主影像窗口中，点击鼠标左键，建立感兴趣多边形第一种点。再次点击鼠标左键，按顺序选取更多边线点，点击鼠标右键来闭合该多边形。鼠标中键可以被用来删除最新定义点，或者删除已经闭合了整个多边形。再一次点击鼠标右键，固定多边形位置。通过选ROIControls对话框顶部相应单选按钮，感兴趣区可以在缩放窗口中被定义。感兴趣区域定义完后，它就会在对话框可用区域列表中显示出来，同步显示尚有感兴趣区域名字，颜色以及所包括像素总数。3要定义一种新感兴趣区，点击NewRegion按钮。点击Edit按钮，在打开EditROIParameters对话框中，可以输入感兴趣名字，选取感兴趣区颜色和填充方式。按照上述环节，定义新感兴趣区。对分类影像定义6个地类。图5.3感兴趣对话框及训练样本4对训练样本进行保存，在ROITool对话框中，选取file→saveROIs。5在SaveROIstofile对话框中，选取SelectAllItems，点击choose选取保存途径和文献名。6点击ok保存文献。图5.4SaveROIstofile对话框最小距离法分类（MinimumDistance）最小距离分类法使用了每个感兴趣区均值矢量，来计算每一种未知像元到每一类均值矢量欧氏距离。除非顾客指定了原则差和距离阈值，否则所有像元都将分类到感兴趣区中最接近那一类。1在ENVI主菜单中，选取