ERDAS+ArcGis论文技术汇总

硕士毕业论文过程中所遇到的技术问题及解答汇总2016.02.231：遥感影像的下载（1）原始数据下载网站：地理空间数据云（2）打开地理空间数据云网站，登陆后——高级检索——选择行政区（安徽省六安市）、时间（2005,2015年）、数据集（2005对应landsat7ETMSLC-off,2015对应landsat8OLI-TIRS）——搜索。（3）确定所需区域的地理位置，具体可以在操作——更多——显示地理位置。选择云量较小的数据进行下载。如果下载数据不能包含所有的研究区域，则需要下载多幅影像在面进行拼接。（4）对于landsat7ETMSLC-off的遥感图像需要进行带条修复模型，按照下图设置好相关的信息，进行条带修复。进行土地利用类型分析，选择345波段，其他时像数据中选择云量较小的数据最后在模型服务——我的计算任务，查看任务进度。并将数据下载下来整理文件夹。2：遥感影像波段融合（1）打开Erdas，新建一个viewer.（2）开始对之前下载的某一时段的遥感影像中的543波段进行融合。具体操作如下主菜单Main——ImageInterpreter——Utilities——Layerstark（2）添加需要融合的波段，每添加一个波段需要点击Add添加按钮，指定输出文件夹。合并完成之后点击OK退出。（3）在Viewer中打开合并好的文件，使用放大缩小图标进行查看。3：遥感影像的条带处理（1）在Erdas菜单栏上选择Main——ImageInterpreter——RadiometricEnhancement——DestripeTM,按照下图设置参数。（2）完成之后点击OK，在Viewer中查看去条带的影像。4：在GIS中创建适应于裁剪遥感图像的行政区边界及面域。（1）所需资料：大范围的shp文件（六安市乡镇边界线图层矢量图）。（2）打开ARCMAP，创建文件夹链接，以便打开储存文件。具体方法：主菜单——窗口——目录——文件夹链接——右键选择相应文件夹（3）将所需文件拖入窗口中加载，将相应图层可以编辑模式。具体做法：右键图层——编辑要素——开始编辑——出现编辑窗口选择不必要的线，删除。（4）将范围外面的线打断，具体做法：先选中一条待分割的线，点击编辑器的分割工具，在端点处打断。将多余的线删除（5）同时检查断口是否连接闭合。双击需要连接的线调整，将确定连接上的两条线并。（6）线转面：确定所有线闭合之后，地理处理——Arctoolbox——数据管理工具——要素——要素转面形成面之后，单击面域，检查是否形成单独面域。形成单独面域5：坐标系及投影变换：查看坐标系，并将gis坐标系转化成为ERDAS可用的坐标系。从GCS-WGS-1984转变为UTM/1984具体做法：查看shp文件的坐标系，右键点击目录中的shp文件查看属性——XY坐标确定需要转换的坐标系及区间，具体做法：在工具栏中单击ArcToolbox——数据管理工具——投影和变换——要素——投影，，打开投影变化框输出坐标系单击右侧的浏览按钮，进入到“空间参考属性”对话框，在对话框中单击“选择”，打开“浏览坐标系”对话框，依次找到ProjectedCoordinateSystems→UTM→WGS1984→Northern/SouthernHemisphere→WGS1984UTMZone***.prj，其中的UTMZone***代表不同的投影分度带，根据城市的经纬度坐标选择，具体选择方法如下：北半球地区，选择最后字母为“N”的带（在NorthernHemisphere文件夹中），南半球地区选择最后字母为“S”的带（在SouthernHemisphere文件夹中）；带数=经度整数位/6的整数部分+31（东经为正值，西经为负值）如：北京约在东经116°24'27.09，带数=116/6+31=50，50N，即WGS1984UTMZONE50N，选中后依次单击确定返回。六安市也是WGS1984UTMZONE50N。单击“确定”执行投影变换，自动加载到列表内容中。6：图像增强——光谱增强具体做法：点击interpreter——SpectralEnhancement——principalcompent，弹出对话框，按照对话框提示输入数据。单击OK确定。7：裁剪图像在ERDAS中裁剪遥感图像。利用行政面域矢量图shp文件，（采用增强后的影像作为待裁剪影像）（额……好吧，增强之后还不如不增强的，换回原来的）打开待裁剪的遥感图像，再打开六安市规划区矢量图，发现覆盖在遥感图像上。选中六安市矢量图，图像变为亮黄色，点击菜单栏的AOI——CopyselelacationtoAOI，再次选中六安市矢量图，出现矩形边框。将AOI保存到指定文件夹。File——save——AOIlayeras，点击主菜单中的Datareparation——subsetImage，在AOI处点击选择指定文件AOIfile。确定裁剪图像。8：监督分类（1）打开需要进行分类的图像（2）打开分类模板编辑器，具体步骤：点击主菜单上的Classifier图标——Classification——SignatureEditor——窗口主面板——AOI——tools（3）点击主菜单上的Classifer——SignatureEditor，跳出SignatureEditor面板。点击主窗口菜单栏上的AOI——Tools，用多边形工具画出一类用地的训练区，点击加入特征文件编辑面板。（4）画一块，添加一块，每类用地在重点地域和全区域选5-10块。点击，可将选择的几块地归为一类。归并后，点击Edit——delete将原有的删除。并重命名地类。（5）将全区域土地分为几种地类。保存分类特征信息：在SignatureEditor面板上，点击File——SaveAs，在弹出的对话框中选择路径并为文件命名，点击OK。（7）评价分类模板：可能性矩阵评价工具（8）执行监督分类，具体做法：ERDAS图标面板工具条:点击Classifier——Classification菜单——SupervisedClassification菜单项（9）定义阈值，确定哪些像元最没可能正确分类，进而对分类结果进一步优化。Classification——Threshold——9：建模（1）查看各类用地的属性值：点击主窗口的Raster——ProfileTools，在跳出的对话框选择Spatial，跳出SPATIALPROFILE面板，用在要查询的用地上拉条线，即可显示用地属性值。建筑用地林地农地水域（2）构建模型，具体做法：点击主菜单Modeler——ModelMaker——NewModel窗口，用NewModel窗口工具条构建运算模型，头和尾分别输入文件和输出文件，中间部分双击中间过程〇，跳出FuncitonDefinition对话框，在Funciton选择Conditional，再选择第二个条件EITHERarg1IF(test)ORarg2OTHERWISE语句。其中text部分输入文件名==前面计算的值，如EITHER1IF($n1_2005fl==9)OR0OTHERWISE，点击OK。再点击红色闪电图标运行。相同步骤输入建筑。林地、、农地、水体的值，运行完成，建模后的成果10：将栅格数据转换成矢量数据点击主菜单上的Vector——RastertoVector，在跳出的对话框中选择输入和输出的文件，点击OK，跳出RastertoARC/INFOCoverage对话框，选择Map和多边形Polygon，点击OK.生成info和四种类型的五个文件的矢量文件。11：使用gis统计数据并制作专题图（1）导入刚刚转换好的矢量数据（2）打开图层属性表，选出属性为1的数据，实际上就是选择实际的用地，去掉边框底图。（3）在属性表主菜单中中打开Options，按属性选择，选择GRID-CODE为1的数据。并显示所选数据，剔除不需要的数据。选中AREA—Statistic统计，弹出对话框以此类推各年份各种数据的属性表2015建筑2015林地2015农地2015水域（4）将矢量图转换成为shp格式输出成果加载之后可调整颜色以及边框样式（5）修改Gird—code值，便于不同地类的区分选中图层打开属性表，选中GIRD-CODE列，右击选择Fieldcalculator,（不同图层设置不同的值）。（5）合并图层，具体步骤地理处理——Arctoolbox——数据管理工具——常规——合并改变合成图像的颜色，具体步骤：右键图层——属性，调整好之后，点击应用——确定12：FRAGSTATS景观格局知识分析将矢量数据转换成为栅格数据选择arctoolbox——转换——栅格——要素转栅格选取18个景观格局指数注意：文件的读取路径中不能含有中文字符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，即去除背景值。具体步骤为GISAnalysis——Databasequery——reclass第一列是赋予的新值，第二列和第三列是旧地类编码的范围。由于栅格中的用地代码是1,2,3,4，于是写成如上图所示，背景值为-9999—0，重新赋予新值0（3）获取马尔科夫链矩阵方法:Modeling——Environmental/Simulationmodels——MAKOV（4）实现CA-MAKOV模型的模拟预测方法:Modeling——Environmental/Simulationmodels——CA-MAKOV（5）GISAnalysis——DatabaseQuery——CROSSTAB计算转移面积14将IDRISI的预测的25年文件，放在GIS中继续分析将idrisi的文件转化成为ASCII文件（1）将ASCII文件转化成为gis中可用的栅格文件15Arcgis中裁剪栅格图将栅格文件转换成为矢量文件，