土地适宜性评价上机实验、课程.pptVIP

福建省某县耕地适宜性评价 2011年6月 实验内容：福建省某县耕地适宜性评价，制作耕地适宜性评价图和耕地适宜性评价结果表 实验目的：熟悉并掌握土地适宜性评价的一般流程，重点掌握利用ArcGis叠置分析完成加权指数和法的运算。 数据准备 Dem ：高程 Dimao： 地貌类型 Aspect： 坡向 Slope： 坡度 Xiangzhen： 乡镇边界 说明：本实验数据只适用于演示实验过程，并不代表真实值 一、建立技术路线 划分评价单元 空间分析 评价指标体系 评价因子权重 专题制图 二、划分评价单元 采用网格法作为评价单元的划分方法 把各参评因子图层划分为一定尺寸的网格,以网格作为评价单元。 考虑到现有数据的精度，选择50m×50m的网格作为评价精度。 各参评因子层为GRID格式,格网大小统一为50m×50m 数据预处理：建立slope图层的坐标系统。 三、评价指标体系 评价因子 因子分级 因子量化 地貌类型 低河漫滩、湖积平原等 100 高河漫滩、洪积平原 80 河流高阶地 40 坡度（°） 小于3 100 3~5 80 5~8 60 8~15 40 15~25 20 大于25 0 高程（m） 小于100 100 100~200 80 200~300 65 300~400 50 400~500 35 大于500 20 坡向 平台 100 南坡 80 东南坡、西南坡 70 东坡、西坡 60 东北坡、北坡 40 西北坡 20 对slope指标进行量化处理 加载slope图层到ArcMap中，此时，slope图层的value值表示实际坡度值，其中0表示坡度为0-3度，3表示坡度3-5度，以此类推。 前述评价指标体系表中已经对坡度因子进行了量化，例如坡度小于3度的量化值为100。在ArcGis中需采用表面分析工具reclassify实现指标的量化。具体步骤如下： 第一步，调用reclassify工具，选择input raster 和 reclass field。 第二步：将新值（既因子量化值）输入New value栏中。例如旧值0对应的新值为100。 新生成图层为reclass_slope 四、评价因子权重 评价因子权重的确定方法很多（经验法、层次分析法），一般需要研究者具有较丰富的经验。本实验不涉及。 下表给出了本实验使用的权重系数 评价因子 权重系数 地貌类型 0.3 坡度 0.3 高程 0.2 坡向 0.2 五、评价模型 本实验采用加权指数和法作为耕地适宜性评价模型，评价模型的表达式如下: S为评价单元耕地适宜性综合总分值；n为评价因子数；Wk为第k个评价因子的权重；Vk为第k个评价因子的量化值。 耕地适宜性等级划分 得出每一个评价单元（50*50m）的S值后，还需要根据实际情况，依据S值的大小判断每一个评价单元是否适合于耕地。既耕地适宜性等级的划分。划分依据见下表。 适宜 不适宜 高度 中度 勉强 S值范围 88-100 70-88 60-70 60 六、ArcGis中的评价过程 第一步：将4个评价因子图层加入ArcMap,调用raster calculator 工具 （工具路径：spatial analyst/ raster calculator） 在计算器空白栏中输入： [dimao] * 0.3 + [reclass_slope] * 0.3 + [dem] * 0.2 + [aspect] * 0.2 点击evaluate完成栅格计算，获得初步的耕地适宜性评价图层（calculation） 第二步：将耕地适宜性评价图层分为高度适宜、中度适宜、勉强适宜和不适宜四个等级。 调用reclassify工具，旧值60对应新值1（不适宜）、旧值60-70对应新值2，旧值70-88对应新值4，旧值88-100对应新值4（高度适宜）。 适宜性分级 Reclass of Calculation图层 七、评价结果分析 制作耕地适宜性评价结果表。既统计该县分别属于高度、中度、勉强、不适宜耕地和未评价区域的面积。 第一步：图层切割。由于此时耕地适宜性图层和xiangzhen图层并非完全重合。需要使用toolbox中的工具对耕地适宜性图层进行切割，从而便于统计该县范围内的不同耕地适宜程度的土地面积。 工具路径：ArcToolBox / spatial analyst / extraction / extract by mask。打开对话框设置如下图 前 后，图层名：某县耕地适宜性 第二步：统计面积 首先计算该县的面积，在xiangzhen图层的属性表中完成。 在属性表中，Addfield name：area； type：double。 然后选中area字段， Field calculato