实验十一：耕地地力评价
一、实验目的：
采用层次分析法以有机质、有效磷、基础设施、灌溉能力、高程、坡向为因子对研究区内的耕地进行综合地力评价 二、实验过程：
1、对“土壤.shp”和“土地利用.shp”和“村界1”和“村界2”进行联合； 2、把村界1和村界2合并，再以两张村界图层为裁剪要素裁剪上步结果； 3、对上一步结果选择，表达式为“Landuse3=旱地”，得到“评价单元”； 4、把“评价单元”导出再添加进来，对评价单元进行爆破处理；
5、为操作方便将“评价单元.shp”的属性表中的多余字段删除，仅保留新村名、土壤类型、Landuse3、Area；
6、利用“计算几何”重新计算Area（面积）字段，并将Area<20的记录删除 7、加载“测试点.shp”，以“有机质”为Z字段进行反距离权重处理，输出像元大小设为5，在环境中设置输出范围为“评价单元”；
8、对上步结果和“评价单元”执行以表格显示分区统计操作，统计类型为mean，得到统计表1；
9、连接统计表1和“评价单元”的属性表，并将有机质的均值字段复制到“评价单元”的属性表中，再增加字段“youjizhi”，在字段计算器中执行“youjizhi=0.8843*Ln(有机质_mean)-1.4180”，后断开连接； 10、 重复6—8步操作，操作字段为“有效磷”，字段计算器中的算式为“0.5530*Ln（有效磷_mean）-0.9080”
11、打开“测试点”属性表，利用选择工具将“农田基设施”字段属性值“不配套”、“基本配套”、“无农田建设”依次赋值为1,2,3；将“灌溉能力”字段属性值“可灌”、“能灌”、“无灌”依次赋值为1,2,3； 12、以“测试点”为输入要素创建泰森多边形；
13、选择连接操作中的“基于空间位置的另一图层的数据”，连接泰森多边形和”
测试点”属性表；
14、将上步结果分别以“农田基设施”和“灌溉能力”字段为值字段，输出像元
大小为5，转为栅格图；
15、分别对“农田基设施”和“灌溉能力”的栅格图与“评价单元”进行“以表
格显示分区统计”操作，统计类型为Majority，得到“分区统计灌溉”和“分区统计设施”两张统计表；
16、分别连接“分区统计灌溉”和“分区统计设施”与“评价单元”属性表，将“灌溉_Majority”和“设施_Majority”字段复制到“评价单元”属性表； 17、再增加字段“灌溉”和“设施”，利用选择工具依次将 “灌溉_Majority”
的1、2、3分别赋值为0.59、0.81、0.39，将“设施_Majority”的2、3分别赋值为0.36、0.8；
18、加载四张等高线图层，创建TIN，坐标系保持与“评价单元”相同——
beijing-1954-GK-Zone19；再转为栅格图——DEM；
19、对“评价单元”与“DEM”进行“以表格显示分区统计”操作，统计类型为
Mean，得到“高程”统计表，连接“评价单元”与“高程”属性表，将“Elev_mean”字段复制到“评价单元”； 20、增加字段“海拔”，在字段计算器中执行“海拔=0.0005*Elev_mean-0.125； 21、利用“DEM”得到研究区的坡向图“Aspect_DEM”,对“Aspect_DEM”进行重
分类，坡度从0-360分为1-9；
22、对“评价单元”与“Reclass_asp”进行“以表格显示分区统计”操作，统
计类型为Majority，得到“分区统计坡向”统计表；
23、连接“分区统计坡向”与“评价单元”属性表，将“坡向_Majority”复制
到“评价单元”，并增加字段“坡向”，利用选择工具将“坡向_Majority”的属性值1-9依次在“坡向”中赋值为0.4、0.63、0.74、0.93、0.99、0.73、
0.5、0.5；
24、在“评价单元”属性表中增加“综合地力指数”字段，在字段计算器中执行“综合地力指数=youjizhi*0.25+youxiaolin*0.17+设施*0.08+灌溉*0.3+海 拔*0.05+坡向*0.15”； 海拔 0.05 立地条件 坡向 0.15 有机质 0.25 肥力状况 有效磷 0.17 灌溉能力 0.3 土壤管理 农田基础设施 0.08
25、在图层属性——符号系统——数量——分级色彩中对“评价单元”进行分类处理，最终结果如下：
三、实验总结
实验十一 耕地地力评价
