ArcGIS进行土地利用现状数据库面积量算研究..doc

ArcGIS进行土地利用现状数据库面积量算研究 一、实验数据结构 土地利用数据库常见的图层有地类图斑、线状地物、零星地物和行政区划四层数据，各层数据结构如下表： 数据层 字段名称 字段代码 字段类型 字段长度 小数位数 DLTB 地类图斑 标识码 ID Int 10 地类编码 DLBM Char 4 权属代码 QSDM Char 12 地类面积 DLMJ Float 15 2 线状地物面积 XZDWMJ Float 15 2 零星地物面积 LXDWMJ Float 15 2 计算面积 JSMJ Float 15 2 XZDW 线状地物 标识码 ID Int 10 地类编码 DLBM Char 4 权属代码1 QSDM1 Char 12 权属代码2 QSDM2 Char 12 地类面积 DLMJ Float 15 2 地物长度 DWCD Float 15 1 地物宽度 DWKD Float 15 1 扣除系数 KCXS Float 5 1 LXDW 零星地物 标识码 ID Int 10 地类编码 DLBM Char 4 权属代码 QSDM Char 12 地类面积 DLMJ Float 15 2 XZQH 行政区划 标识码 ID Int 10 权属代码 QSDM Char 12 其中，DLTB、XZDW、LXDW的DLBM，XZDW的DWKD，LXDW的DLMJ，以及XZQH的QSDM都是已知属性数据，是在土地调查中直接输入的，其他数据均可通过数据分析得出的。 建立一个新的文档，添加四层数据，用Editor编辑生成如图所示的要素，输入已知属性数据，进行地图整饰后保存地图文档TDLYXZ.mxd。 二、赋QSDM 在ArcGIS主要运用ArcToolBox—AnalysisTools—Overlay—Identity这一工具实现赋QSDM。 用XZQH对DLTB进行Identity： 用XZQH对XZDW进行Identity，注意选中“Keep Relationship”： 用XZQH对LXDW进行Identity： 利用属性表中FieldCalculator功能给将Identity后的QSDM赋到原有QSDM中，并删除其他多余字段，保留原有属性结构。其中，因为某些XZDW是XZQH分界线，所以有QSDM1和QSDM2两个权属代码。 这样就完成了对DLTB、XZDW和LXDW的QSDM赋值，接下来对DLTB_QSDM、XZDW_QSDM和LXDW_QSDM三层数据继续处理。 三、XZDW的DLMJ和KCXS求取 1、KCXS求取 这里我们先求取KCXS，因为有些XZDW是DLTB分界线，而有些不是，所以当XZDW是DLTB分界线时KCXS选取0.5，当XZDW是不是DLTB分界线时KCXS选取1.0，KCXS的判别同样采用Identity这一工具来实现。 用DLTB_QSDM对XZDW_QSDM进行Identity，注意选中“Keep Relationship”： 打开属性表对KCXS进行判断，对KCXS进行FieldCalculatorAdvanced运算： 代码为： if [RIGHT_DLTB_QSDM] = [RIGHT_DLTB_QSDM_1] then output=1 else output=0.5 end if 即判断线段两侧的DLTB是否为同一图斑，若是则KCXS为1，否则KCXS为0.5。 2、DLMJ求取 在属性表中用FieldCalculator由SHAPE_Length向DWCD赋值，由DWCD和DWKD的乘积向DLMJ赋值，这样就求取了DLMJ。 新建字段KCMJ，数据类型为Float，赋值为DLMJ和KCXS的乘积，以备下一步使用。 删除其他多余字段，除KCMJ外保留原有属性结构。接下来对XZDW_QSDM_KCXS进行操作。 四、DLTB的DLMJ求取 经过以上几步，XZDW和LXDW的属性字段已经求取完毕，剩下DLTB的DLMJ等字段没有求取，接下来要运用ArcToolBox—AnalysisTools—Overlay—SpatialJoin这一工具来实现。 将XZDW_QSDM_KCXS用SpatialJoin到DLTB_QSDM中： 其中要注意三点： 1、在FieldMapofJoinFeatures中删除SHAPE_Area后不需要的字段。 2、在FieldMapofJoinFeatures中删除XZDWMJ原有字段，替换成XZDW_QSDM_KCXS的KCMJ字段，右击XZDWMJ点击Properties，在OutputField