ERDAS+ArcGis论文技术汇总

1、硕士毕业论文过程中所遇到的技术问题及解答汇总2016.02.231：遥感影像的下载（1）原始数据下载网站：地理空间数据云（2）打开地理空间数据云网站，登陆后高级检索选择行政区（安徽省六安市）、时间（2005,2015年）、数据集（2005对应landsat7 ETM SLC-off, 2015对应landsat8 OLI-TIRS）搜索。（3）确定所需区域的地理位置，具体可以在操作更多显示地理位置。选择云量较小的数据进行下载。如果下载数据不能包含所有的研究区域，则需要下载多幅影像在面进行拼接。（4）对于landsat7 ETM SLC-off的遥感图像需要进行带条修复模型，按照下图设置好相关的

2、信息，进行条带修复。 进行土地利用类型分析，选择345波段，其他时像数据中选择云量较小的数据最后在模型服务我的计算任务，查看任务进度。并将数据下载下来整理文件夹。2：遥感影像波段融合（1） 打开Erdas，新建一个viewer. （2）开始对之前下载的某一时段的遥感影像中的543波段进行融合。具体操作如下主菜单MainImage InterpreterUtilitiesLayerstark（2） 添加需要融合的波段，每添加一个波段需要点击Add添加按钮，指定输出文件夹。合并完成之后点击OK退出。（3） 在Viewer中打开合并好的文件，使用放大缩小图标进行查看。 3：遥感影像的条带处理（1）在

3、Erdas菜单栏上选择MainImage InterpreterRadiometric EnhancementDestripe TM ,按照下图设置参数。（2）完成之后点击OK，在Viewer中查看去条带的影像。4：在GIS中创建适应于裁剪遥感图像的行政区边界及面域。（1） 所需资料：大范围的shp文件（六安市乡镇边界线图层矢量图）。（2） 打开ARCMAP，创建文件夹链接，以便打开储存文件。具体方法：主菜单窗口目录文件夹链接右键选择相应文件夹 （3） 将所需文件拖入窗口中加载，将相应图层可以编辑模式。具体做法：右键图层编辑要素开始编辑出现编辑窗口选择不必要的线，删除。（4） 将范围外面的线打

4、断，具体做法：先选中一条待分割的线，点击编辑器的分割工具，在端点处打断。将多余的线删除（5） 同时检查断口是否连接闭合。双击需要连接的线调整，将确定连接上的两条线并。（6） 线转面：确定所有线闭合之后，地理处理Arctoolbox数据管理工具要素要素转面形成面之后，单击面域，检查是否形成单独面域。形成单独面域5：坐标系及投影变换：查看坐标系，并将gis坐标系转化成为ERDAS可用的坐标系。从GCS-WGS-1984转变为UTM/1984 具体做法：查看shp文件的坐标系，右键点击目录中的shp文件查看属性XY坐标 确定需要转换的坐标系及区间，具体做法：在工具栏中单击ArcToolbox数据管理

5、工具投影和变换要素投影，打开投影变化框输出坐标系单击右侧的浏览按钮，进入到“空间参考属性”对话框，在对话框中单击“选择”，打开“浏览坐标系”对话框，依次找到Projected Coordinate SystemsUTMWGS1984Northern/ Southern HemisphereWGS 1984 UTM Zone *.prj，其中的UTM Zone *代表不同的投影分度带，根据城市的经纬度坐标选择，具体选择方法如下：   北半球地区，选择最后字母为“N”的

6、带（在Northern Hemisphere文件夹中），南半球地区选择最后字母为“S”的带（在Southern Hemisphere文件夹中）；带数=经度整数位/6的整数部分+31（东经为正值，西经为负值）如：北京约在东经116°24'27.09"， 带数=116/6+31=50，50N，即WGS 1984 UTM ZONE 50N ，选中后依次单击确定返回。六安市也是WGS 1984 UTM ZONE 50N。 单击“确定”执行投影变换，自动加载

7、到列表内容中。6：图像增强光谱增强 具体做法：点击interpreterSpectral Enhancementprincipal compent，弹出对话框，按照对话框提示输入数据。单击OK确定。7：裁剪图像在ERDAS中裁剪遥感图像。利用行政面域矢量图shp文件，（采用增强后的影像作为待裁剪影像）（额好吧，增强之后还不如不增强的，换回原来的）打开待裁剪的遥感图像，再打开六安市规划区矢量图，发现覆盖在遥感图像上。选中六安市矢量图，图像变为亮黄色，点击菜单栏的AOICopy selelacation to AOI，再次选中六安市矢量图，出现矩形边框。 将AOI保存到指定文件夹。Filesave

8、AOI layeras，点击主菜单中的Data reparationsubset Image，在AOI处点击选择指定文件AOI file。确定裁剪图像。8：监督分类（1） 打开需要进行分类的图像（2） 打开分类模板编辑器，具体步骤： 点击主菜单上的Classifier图标ClassificationSignature Editor窗口主面板AOItools（3） 点击主菜单上的ClassiferSignature Editor，跳出Signature Editor面板。点击主窗口菜单栏上的AOITools，用多边形工具画出一类用地的训练区，点击加入特征文件编辑面板。（4）画

9、一块，添加一块，每类用地在重点地域和全区域选5-10块。点击，可将选择的几块地归为一类。归并后，点击Editdelete将原有的删除。并重命名地类。 （5）将全区域土地分为几种地类。保存分类特征信息：在Signature Editor面板上，点击FileSave As，在弹出的对话框中选择路径并为文件命名，点击OK。（7） 评价分类模板：可能性矩阵评价工具（8） 执行监督分类，具体做法：ERDAS图标面板工具条:点击ClassifierClassification菜单 Supervised Classification菜单项（9） 定义阈值，确定哪些像元

10、最没可能正确分类，进而对分类结果进一步优化。ClassificationThreshold9：建模（1）查看各类用地的属性值：点击主窗口的RasterProfile Tools，在跳出的对话框选择Spatial，跳出SPATIAL PROFILE 面板，用在要查询的用地上拉条线，即可显示用地属性值。建筑用地 林地 农地 水域 （2）构建模型，具体做法：点击主菜单ModelerModel MakerNew Model窗口，用New Model窗口工具条构建运算模型，头和尾分别输入文件和输出文件，中间部分双击中间过程，跳出Funciton Definitio

11、n对话框，在Funciton选择Conditional，再选择第二个条件EITHER <arg1> IF (<test> ) OR <arg2> OTHERWISE语句。其中text部分输入文件名=前面计算的值，如EITHER 1 IF ($n1_2005fl =9) OR 0 OTHERWISE，点击OK。再点击红色闪电图标运行。 相同步骤输入建筑。林地、农地、水体的值，运行完成， 建模后的成果10：将栅格数据转换成矢量数据点击主菜单上的VectorRaster to 

12、;Vector，在跳出的对话框中选择输入和输出的文件，点击OK，跳出Raster to ARC/INFO Coverage对话框，选择Map和多边形Polygon，点击OK.生成info和四种类型的五个文件的矢量文件。11：使用gis统计数据并制作专题图（1）导入刚刚转换好的矢量数据（2）打开图层属性表，选出属性为1的数据，实际上就是选择实际的用地，去掉边框底图。（3）在属性表主菜单中中打开Options，按属性选择，选择GRID-CODE为1的数据。并显示所选数据，剔除不需要的数据。选中AREAStatistic统计，弹出对话框以此类推各年份各种数据的属性表20

13、15建筑2015林地2015农地2015水域（4） 将矢量图转换成为shp格式输出成果加载之后可调整颜色以及边框样式（5）修改Girdcode值，便于不同地类的区分 选中图层打开属性表，选中GIRD-CODE 列，右击选择Field calculator,（不同图层设置不同的值）。（5） 合并图层，具体步骤地理处理Arctoolbox数据管理工具常规合并改变合成图像的颜色，具体步骤：右键图层属性，调整好之后，点击应用确定 12：FRAGSTATS景观格局知识分析将矢量数据转换成为栅格数据选择arctoolbox转换栅格要素转栅格选取18个景观格局指数注意：文件的读

14、取路径中不能含有中文字符13：Idrisi景观格局预测（1） 数据格式的转换：矢量文件转化为栅格文件：10-11步骤将shp文件设置不同的ID字段，土地利用类型分别编号为1,2,3.。在arcgis中合并shp，将带有属性的合并矢量文件转化为栅格文件。具体做法：栅格文件转为ASCII文件：Arcgis的raster文件转为IDRISI的raster文件，需要用到Arcgis的ASCII文件进行转换,arctoolbox转换由栅格转出栅格转ASCII在IDRISI中将ASCII文件转化为IDRISI的raster文件：（2） 重分类转换后的栅格具有的属性值是土地利用类型的编码，没有属性的值即背景值是默认的-9999，因此惊喜土地利用重分类，将-9999转换成为0，即去除背景值。具体步骤为GIS AnalysisDatabase query reclass第一列是赋予的新值，第二列和第三列是旧地类编码的范围。由于栅格中的用地代码是1,2,3,4，于是写成如上图所示，背景值为-99990，重新赋予新值0（3） 获取马尔科夫链矩阵方法:ModelingEnvironmental/Simulation mo