《土地适宜性评价上机实验.ppt》由会员分享，可在线阅读，更多相关《土地适宜性评价上机实验.ppt（25页珍藏版）》请在三一办公上搜索。
1、福建省某县耕地适宜性评价,2011年6月,实验内容：福建省某县耕地适宜性评价，制作耕地适宜性评价图和耕地适宜性评价结果表实验目的：熟悉并掌握土地适宜性评价的一般流程，重点掌握利用ArcGis叠置分析完成加权指数和法的运算。,数据准备,Dem：高程Dimao：地貌类型Aspect：坡向Slope：坡度Xiangzhen：乡镇边界说明：本实验数据只适用于演示实验过程，并不代表真实值,一、建立技术路线,划分评价单元,空间分析,评价指标体系,评价因子权重,专题制图,二、划分评价单元,采用网格法作为评价单元的划分方法把各参评因子图层划分为一定尺寸的网格,以网格作为评价单元。考虑到现有数据的精度，选择50
2、m50m的网格作为评价精度。各参评因子层为GRID格式,格网大小统一为50m50m数据预处理：建立slope图层的坐标系统。,三、评价指标体系,对slope指标进行量化处理,加载slope图层到ArcMap中，此时，slope图层的value值表示实际坡度值，其中0表示坡度为0-3度，3表示坡度3-5度，以此类推。前述评价指标体系表中已经对坡度因子进行了量化，例如坡度小于3度的量化值为100。在ArcGis中需采用表面分析工具reclassify实现指标的量化。具体步骤如下：,第一步，调用reclassify工具，选择input raster 和 reclass field。第二步：将新值（既
3、因子量化值）输入New value栏中。例如旧值0对应的新值为100。新生成图层为reclass_slope,四、评价因子权重,评价因子权重的确定方法很多（经验法、层次分析法），一般需要研究者具有较丰富的经验。本实验不涉及。下表给出了本实验使用的权重系数,五、评价模型,本实验采用加权指数和法作为耕地适宜性评价模型，评价模型的表达式如下:S为评价单元耕地适宜性综合总分值；n为评价因子数；Wk为第k个评价因子的权重；Vk为第k个评价因子的量化值。,耕地适宜性等级划分,得出每一个评价单元（50*50m）的S值后，还需要根据实际情况，依据S值的大小判断每一个评价单元是否适合于耕地。既耕地适宜性等级的划
4、分。划分依据见下表。,六、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-8
5、8对应新值4，旧值88-100对应新值4（高度适宜）。,适宜性分级,Reclass of Calculation图层,七、评价结果分析,制作耕地适宜性评价结果表。既统计该县分别属于高度、中度、勉强、不适宜耕地和未评价区域的面积。第一步：图层切割。由于此时耕地适宜性图层和xiangzhen图层并非完全重合。需要使用toolbox中的工具对耕地适宜性图层进行切割，从而便于统计该县范围内的不同耕地适宜程度的土地面积。工具路径：ArcToolBox/spatial analyst/extraction/extract by mask。打开对话框设置如下图,前,后，图层名：某县耕地适宜性,第二步：统计面
6、积首先计算该县的面积，在xiangzhen图层的属性表中完成。在属性表中，Addfield name：area；type：double。然后选中area字段，Field calculator,在field calculator对话框中，advanced前打钩，第一个框中输入一下命令：Dim Output as double Dim pArea as Iarea Set pArea=shape Output=pArea.area 第二个框中输入：output,矢量多边形面积计算,右键点击xiangzhen属性表中的area，调出statistics，可查阅相关统计值，其中sum表示该县的总面积。然后，计算该县四种耕地适宜性评价类型区域的面积以及未评价区域的面积。栅格图层面积可根据像元数目来计算，例如，某个栅格图层的像元大小为100*100m，像元（count）数为1000个，则该栅格图层区域的面积：100*100*1000m2,打开“某县耕地适宜性”的属性表，增加area字段，并计算其面积，在field calculator中输入：COUNT*50*50 其中50*50表示像元大小,表 福建省某县耕地适宜性评价结果,
