ArcGIS进行土地利用现状数据库面积量算研究一、实验数据结构土地利用数据库常见的图层有地类图斑、线状地物、零星地物和行政区划四层数据，各层数据结构如下表：数据层字段名称字段代码字段类型字段长度小数位数DLTB地类图斑标识码IDInt10地类编码DLBMChar4权属代码QSDMChar12地类面积DLMJFloat152线状地物面积XZDWMJFloat152零星地物面积LXDWMJFloat152计算面积JSMJFloat152XZDW线状地物标识码IDInt10地类编码DLBMChar4权属代码1QSDM1Char12权属代码2QSDM2Char12地类面积DLMJFloat152地物长度DWCDFloat151地物宽度DWKDFloat151扣除系数KCXSFloat51LXDW零星地物标识码IDInt10地类编码DLBMChar4权属代码QSDMChar12地类面积DLMJFloat152XZQH行政区划标识码IDInt10权属代码QSDMChar12其中，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，注意选中“KeepRelationship”：用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，注意选中“KeepRelationship”：打开属性表对KCXS进行判断，对KCXS进行FieldCalculatorAdvanced运算：代码为：if[RIGHT_DLTB_QSDM]=[RIGHT_DLTB_QSDM_1]thenoutput=1elseoutput=0.5endif即判断线段两侧的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，在OutputFieldProperties的MergeRule中选择“Sum”。3、MatchOption选择“CONTAINS”。这一步生成了DLTB_QSDM_SJXZDW。同上将LXDW_QSDM用SpatialJoin到DLTB_QSDM中：其中也要注意“CONTAINS”和“SUM”这两个参数的选择。这一步生成了DLTB_QSDM_SJXZDW_SJL