Fragstats与Arcgis协同分析.ppt

导入到ArcGIS中 景观分析LandscapeAnalysis Join表格后 ArcGIS中显示举例 原土地利用图 回旋半径图 景观分析LandscapeAnalysis ArcGIS中显示举例 原土地利用图 回旋半径图 景观分析LandscapeAnalysis Fragstats中帮助文档的描述 ArcGIS中显示举例 类似的，可以对class文件与源图像按照类型进行匹配来join到图像的属性表中 景观分析LandscapeAnalysis ASCII类型 ArcGIS中将栅格转换为ASCII后是不能直接在fragstats分析的，需要打开txt文件，把前面的几行描述信息删掉，只保留栅格信息，才能被Fragstats分析。 须删掉前面那6行信息 ASCII类型 景观分析LandscapeAnalysis ERDAS和8bit类型 貌似Fragstats不支持ERDAS9以后的img格式，所以会报错 ERDAS类型 其帮助文档是这样描述的： 8bit类型 虽然不可以直接用ERDAS类型，但可以选择8bit类型来对img图像进行分析（实际就可以是ERDAS9的文件） 但是这时候一定要选择分类文件，否则即使能够执行，但结果也是不正确的 景观分析LandscapeAnalysis 指定Class类型文件 这里可以指定你自己写的txt或文件或fdc文件，其格式很简单 ClassID, ClassName, Status, isBackground Status 代表是否参与计算， isBackground代表是否作为背景 景观分析LandscapeAnalysis 批处理 如果要对多个镇区分别计算其景观指数的话，需要在Fragstats用到批处理，格式如下： 景观分析LandscapeAnalysis GRID文件的批处理格式 Fragstats帮助文档的关于其他文件格式的批处理格式，基本类似 批处理会遇到的问题 Fragstats不能对一整幅图像自动分各镇区来计算对应的景观指数，所以先要把各镇区的土地利用图裁剪出来 景观分析LandscapeAnalysis 批量裁剪栅格图的方法： 1、split工具分割矢量图（这里要注意有些镇区包含多个polygon时，因为两个polygon名字相同，会导致分割出9999的错误，要修改其名字为为唯一值才可以，如在镇区名前加上FID，因为FID肯定是唯一的） 2、右键点击extract by mask工具来批量裁剪栅格图（这个需时较长） 3、在excel中批量生成Fragstats中需要的文件名 注意：文件路径都不能有空格，否则Fragstats分析时会报错 批处理会遇到的问题 Fragstats中的问题 景观分析LandscapeAnalysis Fragstats批处理的方法和问题： 1、在excel中批量生成Fragstats中需要的文件名（手动输入太麻烦） 2、Fragstats对批量的栅格图生成了结果文件后，以城市用地的斑块密度为例，有些镇区包含城市用地和非城市用地，而有些镇区却只有城市用地或只有非城市用地，这时下图的情况 只有城市用地这一行 有些polygon的面积比栅格图的一个像元还小，造成裁剪出来也是无实际值，Fragstats就无法识别，变成没有这个镇区一样 通常我们需要的只是城市用地的斑块密度，所以要把这两个得到的值分离开来 批处理会遇到的问题 Fragstats中的问题 景观分析LandscapeAnalysis 有些镇区有一行，有些有两行，有些连一行都没有，这样不规则的数据分离，可以用到excel中的宏命令来处理 Excel宏命令步骤 1、office2007为例子 2、视图 3、查看宏 4、输入自定义的宏名字，点编辑 批处理会遇到的问题 景观分析LandscapeAnalysis 关键是找先到定位单位格位置的函数，见下图 Sub forclass() 需要先清除excel内所有空格 Dim t As Integer Dim s As String Dim count As Integer Dim nomalCity_FD As String 对应LID的位置显示PD值 Dim orderCity_FD As String 去掉TYPE为0的重复城市 Dim order As String 序号 t = 1 累计城市数量 s = D:\split\1995output\extract_b 指定文件名和位置 count = 350 设定得到的文件数量 设定输出列--------- nomalCity_FD = f orderCity_FD = h order = g Range(nomalCity_FD + 1).Value = PD(city