创新实验评价单元及插值

评价单元制作与插值需要什么数据？目的是什么？采用什么技术？方法与操作？评价单元与插值的关系？得到什么结果？存在什么问题？创新实验系列一、基础图件检查

1.相交图层准备

2.使用相交工具

3.整理属性表

4.重建拓扑（其中3和4哪个先做？）一、基础图件检查主要检查内容1.图件比例尺、评价对象的边界（选择来源可靠、可以相互印证的图件提取边界，并作为以后制图的边界，以后合并问题）1.1土壤图边界将土壤图调入ArcMap，“数据”－“导出数据”，形成一个新文件,如“boundary.shp”，启动编辑，编辑任务选择“修改要素”，再用“要素选择”工具选择全部多边形，在编辑器中使用“合并”命令把要素全部合并形成土壤图边界。一、基础图件检查1.2土地利用现状图检查并转换为BJ54坐标系，形成比较文件，如“boundary_used.shp”，确定与土壤图具有相同的边界。如果不相同要依据原始栅格图进行必要修改和校正。1.3行政区图（WGS84）按照相同方法处理。二、评价单元生成2.使用相交工具从工具箱中选择，“AnalysisTools”－“Overlay”－“相交”，下图：二、评价单元生成二、评价单元生成打开需要求交的

2个图层，二、评价单元生成单击确定后出现以下执行界面，完成大约需要205秒。二、评价单元生成3.整理属性表（现在做还是最后做！）求交后生成文件，如“Intersect”，生成许多小块多边形，多边形数目增加？太小的多边形如何处理？整理属性表，把不需要的字段去掉，保留必要的属性（利用现状图，利用类型及其代码，土壤图，不同土属或土类及其代码）等的集合。其中：土壤图：土属或土种名称、代码，评价需要的基本养分属性数据；利用现状图：土壤利用分类代码；保留耕地（水田和旱地）、和非耕地；其它还有政区图的要素，如地名、道路和水系等，保存为单独的图层，留待出图时用。二、评价单元生成4.重建拓扑评价单元初步做完后，还要重建拓扑。可以打开已经过空间叠置的评价单元文件，会发现有的多边形相互重叠，没有按照求交结果分离成相互独立的多边形。如果不进行处理，则评价单元的属性不能获得独立值，而且在面积统计时会出现错误。二、评价单元生成

放大可见，虽然求交后多边形进行了叠加，但是原来的一个多边形仍然联合在一起，因此要重新建立拓扑。二、评价单元生成

重建拓扑的方法是：“编辑器”－“开始编辑”，“编辑器”－“更多的编辑工具”－“拓扑”，单击以下图标，用鼠标划框选择全部多边形，

此时在该图标旁边的图标可选，单击即可完成重建拓扑。

二、评价单元生成此时可见，重建拓扑后多边形各自独立.二、评价单元生成选择“编辑器”－“停止编辑”。打开属性表，可见多边形数目比重建前多了许多,再使用字段计算器计算一下，也会发现增加了许多面积很小的多边形。这说明，进行相交叠加处理后，如果不进行重建拓扑，会得到错误的结果。对于小于一定面积的多边形如何处理？合并，采用什么策略？本实验不作要求，感兴趣可以自己找方法，自己做！三、插值方法1.插值的意义和要求插值的意义：把采样点上的数据扩充到采样区域整个面上，以栅格图形式表达。因此，栅格图像中栅格的大小取决于你需要的面上数据的密度。在耕地地力调查中，每一个评价单元都需要一个插值结果，因此实际上由评价单元的最小面积决定栅格的大小。在ArcMap中打开“pingjia_unit.shp”图层，利用属性表，选择字段统计，可见总面积和最小面积的统计值，如何处理最小面积？

三、插值方法2.计算需要的栅格数量那么至少应该有多少栅格，如果评价单元都是正方形，总面积为34492534.78平方米，最小多边形的面积是10000平米，则34492534.78平方米/10000平方米，应该有3449.25个栅格。由于评价单元不会都是正方形，按照长宽比为5计算，至少要有3449.25×5＝17246.25个栅格。要划分出这么多栅格，需要多少行列，开方即可，得行列数目必须分别大于132个，据此可以进行插值设置。三、插值方法3.插值方法3.1打开插值要求的基本文件:包括：

xx.shp，采样点及其点属性;boundary.shp，目标插值区域。三、插值方法3.2启动ArcMap，添加空间分析模块；3.3设置插值参数步骤一：“空间分析”－“选项”，出现右边的设置对话框；可以选可以不选。选择时，后面操作不对边界外取值。三、插值方法插值范围设置必须保证靠近边界的样点与边界之间的评价单元都能获得数据。三、插值方法栅格大小与数目设置调整栅格大小显示栅格数目需要多少行列？三、插值方法需要的行列数目5050时大于1314因此需要使用50×50象素的栅格才能满足基本要求。至此，插值设置完成。三、插值方法3.5插值插值范围，乡镇边界内。不对采样点区域外（乡镇边界外）取值。三、插值方法选择插值方法和插值对象5种属性分别插值。选择实际的字段名插值栅格图存盘，以该属性字段名存盘。三、插值方法插值结果三、插值方法

使用同样的方法插值采样点的所有属性，以属性名称作为栅格文件名称。每插值一个属性便会得到2个文件2个目录，其中包含系列文件。此目录插值前手工建立插值后自动建立的目录三、插值方法info文件夹内容SO120201属性文件夹内容三、插值方法

按照相同方法对所有需要插值的属性进行插值，生成相应的栅格图。存入自建目录下：自建的插值属性字段名目录。三、插值方法3.6区域统计把评价单元文件调入，“pingjie_unit_xx.shp”，作为插值对象。打开属性表增加字段，以便所有点属性插值后进行链接。增加的字段作为索引使用，字段属性设置为长整型，给字段内容赋值使等于FID。三、插值方法

选择“空间分析”－“区域统计”，出现以下对话框：评价单元文件增加的索引字段选择插值生成的栅格文件从栅格文件提取的属性值，数据库格式，实例采用属性名带“—1”三、插值方法

此操作需要较长时间，属性提取完成后，弹出结果属性，总共获得24439个数据，与24441个评价单元比较，还差2个评价单元没有获得值。据此结果，可以再改变插值的栅格大小，如从50×50变为45×45。这里，仅作说明，不再修改了。注意：在给大家的插值文件中已经作了修改，提取属性都与评价单元数目相等，索引字段名称为“FID_1”；另外实际评价单元采用的是缩减了的土壤利用现状图与土壤图求交后的结果，因此评价单元要比本小结的少一些。三、插值方法使用同样的方法，把所有插值生成的栅格图通过“区域统计”提取属性，获得各自的属性文件（DBF格式），属性名带“_1”。它们之间利用新增字段作为系列属性文件的链接依据。索引，连接依据三、插值方法参见文件目录三、插值方法参见文件目录区域统计获得三、插值方法3.7创建、编辑所有提取属性的数据库文件在每个属性插值栅格生成后都有1个目录，提取属性值生成的数据库文件也选择存放在此目录中，现在要利用以前讲过的方法，合并提取的所有属性，形成一个属性文件，再连接到评价单元图上，连接的依据都是增加的ChaZhiId字段。给大家的系列插值文件都是实际应用的插值文件，连接字段是FID_1。连接方法：先打开任意一个经“区域统计”获得的属性文件数据库，另存为Excel文件，去掉不要的字段，然后逐个打开其它文件，把需要的数据拷贝到Excel文件中，更换表头。三、插值方法全部完成后，保存。如“js_all.xls”

索引，连接依据三、插值方法3.8数据连接，形成带插值属性的评价单元将整理形成的包括所有插值属性的Excel文件转变成数据库文件，js_all.dbf。上次小结已