ERDAS操作技巧.docVIP

加载波段成影像：interpreter-utilities-layer stackERDAS操作小技巧：1、配准影像图：从Viewer中 打开两幅图（一幅参照，一幅配准） 从菜单栏Raster中 选Geometric correction（几何校正） 在Set Geometric Model中选Polynomial，后点击ok，打开Polynomial Model Properties对话框，在Parameters中Polynomial Order（多项式次方）中选1或2最少GCP公式：，在Projection（投影参数）中Map Units选Meters 点击Add/Change Projection 在Custom中选择所需的Projection Type，Spheroid Name，Datum Name，Scale factor at central meridian，Longitude of central meridian（可以参考参照图中Imagine info中的信息），Latitude of origin of projection，False easting（一般选500000meters），False northing（一般忽略为零），点击Ok Set Projection from GCP Tool中选择Collect Reference Point From（选择视窗采点模式）中的 Existing Viewer选项，Ok。RMS误差(均方根)=这里：xi和yi是输入的原坐标；xr和yr是逆变换后的坐标定义：RMS误差是指GCP的输入（原）位置和逆变换的位置之间的距离（或者说是在用转换矩阵对一个GCP作转换时，所期望输出的坐标与实际输出的坐标之间的偏差）。若RMS误差以坐标系统的距离来表示的话，比如，RMS误差是2，意味着参考像元与逆转换像元之间的距离是2个像元。2、在影像中要显示符合人的肉眼所见到图像模式：可以从Raster Band combinations下面的波段组合中去选择（如5，4，3或7，4，1）。3、消除金字塔（为了使影像图打开的速度增快）：在打开的图中点击工具栏中第三个image info图标（或者Utility下的Layer Info） 选菜单栏上的Edit Compute Pyramid layers 中选22（或44），然后点击ok。4、选择一小幅或多个AOI区，保存为.img影像：首先打开一幅影像图，从AOI中打开AOI工具条，选择想要的区域，如为多个，可按住shift键，点选多个AOI区，后从DataPrep中选Subset Image，输入Input File,输出Output File，点选AOI，选Viewer（或者，直接在Subset对话框中点击AOI，选择Viewer）,最后点ok。切割image影像图选择Utility下的Inquire Box，选择区域范围，然后选DataPrep下的Subset image，填写输入、输出文件名，点击From Inquire Box，Ok。5、非监督分类：打开一幅图，选择工具栏上的Classfier图标，点击Unsupervise Classfication，输入Input Raster File，Output Cluster File，Output Signature Set Filename（分类模板），Number Classes（分几类），Color Scheme Options，Maximun Iterations（最大循环次数），Convergence Threshold（循环收敛阈值0.950），ok。调出两幅图，第一幅为要被分类的图像，第二幅为即将产生的分类图像。在调第二幅时（不必要再从Viewer中打开，直接从第一幅的工具栏上选打开按钮即可），可以在Select Layer To Add中的Raster Options下面的Clear Display点击去掉，ok。（为了便于Blend、swipe、Flicker显示）。显示分类后影像图的属性：在第二幅工具栏上点击“锤子”图标以打开Raster工具面板，然后点击“表格”图标打开Raster Attribute Editor对话框（或者点击Raster下拉菜单里的Attribute，打开Raster Attribute Editor对话框），在Edit中选Column Properties，（可以在Columns一栏中调整排列顺序，通过Up、Down、Top、Bottom、New、Delete，如果选择Editable复选框，则可以在Title中修改各个字段的名字及其它内容），在 Raster Attribute Editor对话框的属性表中选择各个类别颜色。6、数据叠加显示（Blend，Swipe，Flicker）：打开一幅影像图，在该幅图的视窗中同时打开另一个文件（File Open Raster Layer Select Layer To Add对话框，其中在 Raster Options中不选Clear Display），然后在视窗菜单栏Utility的下拉框中选择Blend，Swipe，Flicker。7、在一幅图上选择指定的要素，并将之保存：打开一幅影像图，用栅格工具面板（Rastertools）中的区域增长生成多边形AOI(Region Growing AOI)选择指定要素，然后选择视窗菜单File下的Save AOI Layer As，保存一个后缀名为.aoi的文件，ok。再在该幅影像图上打开刚才保存的.aoi文件，可以进行数据叠加显示。8、在影像图上选择多个AOI区：打开影像图，在Viewer菜单栏上点击十字光标，将其移到要选择的区域，后点击AOI Seed properties，在Area栏输入适当的数值，点击Grow at inquire（或者选择AOI工具面板中的第五行第二个图标，或者Raster下tools的工具面板中的第四行第二个图标，然后在Area栏输入适当的数值，点击Grow at inquire）；如果要多选几个AOI区，可以点击AOI Invert Polygon，然后移动十字光标到另外要选的区域，点击AOI Seed properties，在Area栏输入适当的数值，点击Grow at inquire。9、将要选择的小范围区域删掉（填充）：打开Raster下面的Fill及AOI或Raster Tools工具面板，选择第二行第一个图标，圈选（小范围），在Function栏里选Majority等。可从Preview中看颜色是否相符，单击Apply即可。查看未知像元值，可选Viewer菜单栏中的Select Profile Tool图标或Raster下拉菜单下的Profile tools（图像剖面工具），ok，在图上拉一条线，即可从Spatial Profile图上看到像元值。10、添加注记（Annotation）：打开一幅影像图，选择File下的New里的Annotation Layer对话框，保存为.ovr的文件。当要进行注记时，在一幅影像图上要同时打开.img和.ovr的两个文件，选择菜单栏上Annotation下的Tools，点选A，注记（英文）。11、三维图像操作：点击Utility下的Image Drape，打开一幅三维影像图。在Utility下的Options(三维显示信息设置，包括DEM、Fog、Background)，Dump Contents to Viewer(三维视窗信息转存到二维视窗，在Image Drape菜单条：View下Link/Unlink With Viewer完成三维到二维的直接连接；在二维视窗菜单条：Utility下Selector Properties可打开Eye/Target Edit对话框)。另外Image Drape菜单条：View下Sun Position(太阳光源参数设置)、LOD Control(详细程度设置)，Position下Current Position(观测位置参数设置)。12、图像对比度调整：菜单条Raster下的Contrast（对比度调整）有8项二级菜单命令（Histogram Equation直方图均衡化处理，Standard Deviation Stretch标准差拉伸处理，General Contrast通用对比度调整，Brightness/Contrast亮度/对比度调整，Piecewise Contrast分段对比度调整，Breakpoints直方图断点处理，Load Breakpoints加载直方图断点，Save Breakpoints保存直方图断点）。13、ArcInfo多边形裁剪以及与影像图的裁剪（将一幅影像图进行掩膜）：首先将影像图进行配准（详见步骤1），然后将Vector图转为Raster图（Main下的Image Interpreter下的Utility下的Vector to Raster,填写设置参数，尤其是将x、y参数改为30，另外在Ignore Zero in Output Stats前打勾），或者，打开Vector层，选取感兴趣的区域生成一个AOI文件（copy selection to aoi），然后再掩膜Mask；直接Copy矢量Vector图层中的感兴趣区域生成新的Vector图层（打开Vector层，选取感兴趣的区域，点击菜单栏Vector下拉菜单中的Copy,再点击File下的New（新矢量图层名），Vector layer，Paste，然后新建一个窗口，打开刚刚创建的矢量图层名，Paste即可，将它配准（打开ImageInfo，点击Edit下的Change Map Info里的Projection,选择，而后，再在Edit下的Add/Change Projection修改配准设置），成为Mask图，Ok。或者将Raster与Vector进行叠加，用AOI进行圈选，并存为.aoi的格式打开Interpreter下的Utility下的Mask，设置完参数（Input原始图，Input Mask图，Output文件名），Ok。14、图像拼接（Image Mosaic）：（a）Main下的Data Preparation下的Mosaic Image，打开Mosaic Tool视窗，点击Edit下的Add Images或工具条上的Add Image图标,打开对话框，输入要拼接的图像文件，在图像拼接区域（Image Area Option）选择Compute Active Area(edge)，点击Add;同样的操作加载要拼接的图像文件，完成后点击Close。（b）在Mosaic Tool视窗，点击工具条上的第一行的第二个图标（Set Input Mode图标），进入设置输入图像模式的状态，根据需要对图像进行上下层调整，调整完成后，在Mosaic Tool视窗图形窗口点击一下，退出图像叠置组合状态。(c)选择Edit下的Image Matching或者选择工具条上第二行的第三个图标，打开图像匹配设置对话框（Matching Options）,在设置匹配方法（Matching Method）中选Overlap Areas(重叠区域匹配)，ok。点选Edit下的Set Overlap Function，打开Set Overlap Function对话框或者选择工具条上的第一行的第三个图标和第二行的第四个图标打开Set Overlap Function对话框，设置相交关系（Intersection Method）：No Cutline Exists(没有裁切线)，设置重叠区像元灰度计算（Select Function）：Average最好用Overlay，可以保持原先每一幅图的原始状态，Apply，Close。最后运行Mosaic工具（Run the Mosaic），Process下的Run Mosaic，设置参数：输出文件名，确定输出图像区域（All）,忽略输入图像值（0），输出图像背景值（0），忽略输出统计值（0），ok。 首先打开Data Preparation下的Mosaic Image的Mosaic Tool视窗；打开AOI Tool图标面板，点击Polygon图标，沿第一幅图像轮廓绘制多边形AOI，将其保存为template.aoi；然后加载Mosaic图像，通过Edit下的Add Images实现，第一步，输入第一幅图像文件，在Image Area Option:选中Template AOI，点击Set，输入template.aoi，ok，同样方法加载第二幅图像，在Image Area Option:选中Compute Active Area(edge)，ok；在Mosaic Tool视窗，点击工具条上的第一行的第二个图标（Set Input Mode图标），进入设置输入图像模式的状态，然后选择Edit下的Image Matching或者选择工具条上第二行的第三个图标，打开图像匹配设置对话框（Matching Options），设置匹配方法（Matching Method）:Overlay Areas，ok；点击Mosaic Tool视窗工具条上的Set Intersection Mode图标（第一行的第三个图标），进入设置图像关系模式，点击选择两幅图像的相交线，使其发亮；接着点击Cutline Selection Viewer图标（第二行的第九个图标）以打开截切线选择视窗（Viewer#3），在其中绘制线状AOI功能及相交区域的外轮廓线，然后点击工具条上的AOI Cutlines图标以打开Choose AOI对话框，定义AOI来源：Viewer后点击Viewer#3，接着点选Edit下的Set Overlap Function，打开Set Overlap Function对话框或者选择工具条上的第一行的第三个图标和第二行的第四个图标打开Set Overlap Function对话框，设置相交类型（Intersection Type）:Cutline Exists，设置选择功能（Select Function）:Cut/Feature，Apply,Close；点击Edit下的Output Options，打开Output Image Options对话框，定义参数：最后运行Mosaic工具（Run the Mosaic），Process下的Run Mosaic，设置参数：输出文件名，确定输出图像区域（All）,忽略输入图像值（0），输出图像背景值（0），忽略输出统计值（0），ok。矢量图层拼接（Match）：点击工具栏Vector图标，在Vector Utilities的下拉菜单下选择Mosaic Polygon layers，选择Input、Output文件名，Ok。15、在一矢量层以不同颜色组合显示：打开矢量层，选View下的polygon symbology，Automatic下的Equal Divisions选择ID、Number Of class(分类数)、 勾上Generate New Styles，Ok。要是想要将有不同颜色组合的矢量层保存，可以选择File下的Save里的Top layer as，输入文件名，就可以保存为.arcinfo的矢量文件。（如果先点击Vector下的Viewer Properties选择ID，在Attribute-based symbology选项上打勾，然后再点击菜单栏Vector下的Symbology，在打开的对话框的菜单栏中，选View下的polygon symbology，Automatic下的Unique Value，选择下一级菜单的ID选项，ok。则要一个一个修改）可以把矢量层保存为img 格式：打开矢量层，点击File下的View to image file，输入一个文件名即可。栅格转矢量时：注意将weed tolerance的容限设为图的分辨率。16、数据输入输出：点击Main下的Import/Export或点击菜单工具条上的第三个图标，打开数据输入输出对话框。输入数据操作：Import,输入数据类型：Generic Binary，输入数据媒体(Media)：File，确定输入、输出文件路径和文件名，Ok,打开Import Generic Data对话框，定义参数：数据格式：(BSQ)BIL，数据类型：Unsigned 8 Bit，图像记录长度：0，头文件字节数：0，数据文件行数：5728，数据文件列数：6920，文件波段数：1，保存参数设置（Save Options）以打开Save Options File对话框，定义参数文件名：*.gen（目的是为了输入其余波段时直接调用该参数），Ok,预览图像效果，打开一个视窗显示输入图像，如果正确，点击Ok，打开Importing Generic Data进程状态条，Ok。重复上述部分过程，依次将多个波段数据全部输入，转换为IMG文件。组合多波段数据：Main下Image Interpreter下Utilities下Layer Stack下Layer Selection and Stacking对话框（或点击Interpreter图标下的Utilities下Layer Stack下Layer Selection and Stacking对话框），依次选择并加载（Add）单波段图像：band3.img，Add,band4.img，Add，band5.img，Add，输出多波段文件：bandstack.img，输出数据类型：Unsigned 8 Bit，波段组合选择：Union，输出统计忽略零值：Ignore Zero In Stats，Ok。TIFF图像数据输入输出：只要在打开TIFF的视窗中将TIFF文件另存为（Save as）IMG文件就可以了。输出JPEG图像数据（Export JPEG Data）：点击Main下的Import/Export或点击菜单工具条上的第三个图标，打开数据输入输出对话框，输出数据操作：Export，输出数据类型（Type）为JPEG：JFIF（JPEG），输出数据媒体：File，确定输入、输出文件路径和文件名，Ok，打开Export JFIF Data对话框，设置参数：图像对比度调整：Apply Standard Deviation Stretch，标准差拉伸倍数：2，图像转换质量：100，在Export JFIF Data对话框中点击选择Export Options，打开Export Options对话框，定义参数：选择波段：4，3，2，坐标类型：Map，定义子区：ULX、ULY、LRX、LRY，Ok,返回Export JFIJ Data对话框，Ok。17、Coverage转为AOI：点击Vector图标，选择Build Vector Layer Topology，后选择AOI下Copy Selection To AOI，选中AOI,然后即可参照步骤4进行切割。18、算法编辑功能：Modeler下Model Maker，输完点击工具栏上红色闪电标志。Model Maker里Condition条件下EITHER IF ( ) OR OTHERWISE里的arg1为原图像，test为条件（如为第二幅图，假设not第二幅图，or 0 otherwise），则输出结果为扣掉第二幅图后输出的影像图。例子：EITHER 1 IF ( 文件名=1 ) OR 0 OTHERWISE，表示将文件名中值为1的信息保留下来，其余赋值为零。而Condition条件下CONDITIONAL () , () , . 里的test1为图名1，arg1为图名1的值；即是将test1赋值为arg1。以此类推。例子：CONDITIONAL (专题图名1=1) 1, (专题图名2=1) 1输出结果为：将各专题图中值为1的数据输出，合并生成一幅专题图，并将其赋值为1。（另外，分析域值，可通过Modeler下的Modeler Maker中的算法进行实现）NDVI算法实现：Model Maker，Open，选择Veg_NDVI.gmd，输入Input Raster和输出Output Raster，即可实现。19、信息依特征提取后合并：用AOI选区域（可多选），然后用Classifier下的Signature Editor进行合并，如果不合并，则在图层拼接时会因各个图层的值不同造成拼接Mosaic失败。亦可将所选区域显示的结果加以保存（View菜单下Arrange layers，选择要删除的层，在提示是否保存时，将其存为img 图层）。另外，在打开的Signature Editor对话框和另一个小框的同时，可以选择圆形进行公式编辑，亦可选择其旁边的图形进行Cell编辑；通过框选AOI，点击Signature Editor对话框工具条上的统计图标，查看其属性Cell值，然后在编辑时，添加较明显波段的层，最小值，最大值，Ok。在提取出的专题图层进行叠加时，应注意在打开该专题图层时，在Raster Option中的Display as选Pseudo color,将该专题图层中无关的类设为透明色，就可进行比较。20、矢量层中添加Projection与TM影像图叠加：一是矢量层的Projection与TM的Projection一致；二是将矢量层的Projection里的Datum设为无定义。21、Arc Tool下可进行图像格式转换：如shape到Arc coverage。22、分辨率融合：选择Main下的Image Interpreter下的Spatial Enhancement中的Resolution Merge。设置参数：融合方法（Principle Component），重采样方法（Bilinear Interpolation），输出数据选择（Stretch Unsigned 8 bit），输出波段选择（Select Layers 1:70），Ok。融合特点：对不同空间分辨率遥感图像的融合处理，使融合后的遥感图像既具有较好的空间分辨率，又具有多光谱特征，从而达到图像增强的目的。注意：融合需要相当大的空间，故最好选择小范围区域进行融合操作。*23、将专题图Thematic中的图斑进行整合处理：Interpreter下的GIS Analysis，先进行聚类统计（Clump），connected neighbors 8， 后去除分析（Eliminate），Eliminate与过滤分析（Sieve）的区别：前者按照周围大部分像元值对其进行填充，而后者则不会，是按照设定的最小值(minimum size)进行过滤 。24、矢量层选项Option的设置：点击Vector的工具面板中第八行第四个图标，打开Option选项设置，Node Snap Dist 1.0；Arc Snap Dist 1.0；Weed Dist 0.1；Grain tolerance 0.1；Intersect；Apply。 由于在新建Vector层中，勾画养殖场闭合曲线时发生严重变形，经指导后知是Vector层的选项未设置好的原因25、专题图层颜色设置：打开图层之前，选择Raster Option中Display as的下级选项中颜色选择为Pseudo Color，后即可在专题图层中进行Attribute颜色属性的修改。26、定义分类模板：（一） 打开模板编辑器并调整显示字段点击工具条上的Classifier图标，Signature Editor，可以打开View下的Columns，对有些意义不大的字段不加以显示。（全选Column字段，按住Shift键的同时分别点击不需要的字段，Apply，Close）（二） 获取分类模板信息a 应用AOI绘图工具在原始图像获取分类模板信息 打开Raster工具面板，绘制多边形AOI，在Signature Editor对话框中点击 + 图标，将多边形AOI区域加载到Signature分类模板中，可以改变Signature Name和Color，重复上述操作以多选几个区，并将其作为新的模板加入到Signature Editor当中。如果对同一个专题类型采集了多个AOI并分别生成了模板，可以将这些模板合并（选择Signature，然后点击合并图标，这时一个综合的新模板将生成，原来的多个Signature同时存在，如果必要也可以删除。） 27、保存任意一幅影像图的任意几个波段打开一幅影像图，从Interpreter下的Utility中选Layer Stack，从里面的Layer中选择需要的波段。另外，还可打开另一幅影像，添加其他的波段，再组合而成。28、保存为专题图打开一幅影像图，在DataPrep下Subset Image里面Output中选者Thematic的格式。29、在ERDAS下滤波：Radar下Radar Interpreter里的Speckle Suppression。另外，在滤波的同时，可以选择Calculate Coefficient of Variation复选框，然后从主菜单Session下的Session log里看最适合的参数。Interpreter下的Spatial Enhancement下Convolution里的Kernet中选择高通、低通等进行滤波。30、格式输出：在ERDAS的输入输出模块中，Import/Export中可以选择需要的类型（如GeoTIFF等）。31、GIS分析： （1）邻域分析（Neighborhood）：是针对分类专题图像，采用类似卷积滤波的方法对图像分类值进行多种分析，每个像元的值都参与定义的邻域范围和分析函数所进行的分析，而邻域中心像元的值将被分析结果所取代。（2）周长计算（Perimeter）：仅对做过Clump聚类的专题图像计算每个类组的周长。（3）查找分析（Search）：是对输入的分类专题图像或矢量图形进行邻近(Proximity)分析，产生一个新的输出删格文件、输出像元的属性值取决于其位置与您选择专题类型像元的接近程度和定义的接近距离，输出文件中您所选择专题类型的属性值重新编码为0，其它相邻区域属性值取决于它们与所选择专题类型像元的欧氏距离。（4）指标分析（Index）：可以将两个分类专题图像或矢量地图数据，按照定义的权重因子（Weighting Factor）进行想加，产生一个新的综合图像文件。（5）叠加分析（Overlay）：根据两个输入