土地适宜性评价上机实验

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,福建省某县耕地适宜性评价,2011,年,6,月,实验内容：福建省某县耕地适宜性评价，制作耕地适宜性评价图和耕地适宜性评价结果表,实验目的：熟悉并掌握土地适宜性评价的一般流程，重点掌握利用ArcGis叠置分析完成加权指数和法的运算。,数据准备,Dem：高程,Dimao：地貌类型,Aspect：坡向,Slope：坡度,Xiangzhen：乡镇边界,说明：本实验数据只适用于演示实验过程，并不代表真实值,一、建立技术路线,划分评价单元,空间分析,评价指标体系,评价因子权重,专题制图,二、划分评价单元,采用网格法作为评价单元的划分方法,把各参评因子图层划分为一定尺寸的网格,以网格作为评价单元。,考虑到现有数据的精度，选择50m50m的网格作为评价精度。,各参评因子层为GRID格式,格网大小统一为50m50m,数据预处理：建立slope图层的坐标系统。,三、评价指标体系,评价因子,因子分级,因子量化,地貌类型,低河漫滩、湖积平原等,100,高河漫滩、洪积平原,80,河流高阶地,40,坡度（,）,小于,3,100,35,80,58,60,815,40,1525,20,大于,25,0,高程（,m,）,小于,100,100,100200,80,200300,65,300400,50,400500,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,为评价因子数；,W,k,为第,k,个评价因子的权重；,V,k,为第,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 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,*,1000m,2,打开“某县耕地适宜性”的属性表，增加,area,字段，并计算其面积，在,field calculator,中输入：,COUNT*50*50,其中,50,*,50,表示像元大小,适宜,不适宜,参与评价区域,总面积,高度,中度,勉强,面积（,km,2,）,百分比（,%,）,表 福建省某县耕地适宜性评价结果,