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

2011-03-30 16:01阅读：

http://blog.sina.cn/dpool/blog/u/1412991875

一、实验数据结构
土地利用数据库常见的图层有地类图斑、线状地物、零星地物和行政区划四层数据，各层数据结构如下表：

数据层	字段名称	字段代码	字段类型	字段长度	小数位数
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
XZQH 行政区划	权属代码	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，在OutputFieldProperties的MergeRule中选择“Sum”。
3、MatchOption选择“CONTAINS”。
这一步生成了DLTB_QSDM_SJXZDW。
同上将LXDW_QSDM用SpatialJoin到DLTB_QSDM中：

其中也要注意“CONTAINS”和“SUM”这两个参数的选择。
这一步生成了DLTB_QSDM_SJXZDW_SJLXDW
在属性表中利用FieldCalculator由SHAPE_Area向JSMJ赋值，由JSMJ与XZDWMJ和LXDMJ的差向DLMJ赋值，这样就求取了DLMJ。
删除其他多余字段，保留原有属性结构。
通过上面的操作，我们填写完成了DLTB、XZDW和LXDW所有的属性数据，生成最终文件为DLTB_QSDM_SJXZDW_SJLXDW、XZDW_QSDM_KCXS和LXDW_QSDM三个，已经基本大功告成！

五、各XZQH的DLMJ汇总
最终，我们要生成各行政单位各地类的一个汇总表，这一步我们无法ArcGIS中完成，需要借助Excel中数据透视功能。
1、将最总生成的三个数据层导出dbf属性表。
2、DLTB和LXDW的属性表只保留DLBM、QSDM和DLMJ三个属性，而XZDW要保留DLBM、两个QSDM和KCMJ。
3、对XZDW的属性表进行单独处理，将两个QSDM排成一列，并补充对应的DLBM和KCMJ。
4、将DLTB、XZDW和LXDW处理好的属性表合并到一个Excel表格文件中。
5、运用数据透视功能，最终完成了DLMJ的汇总工作！

后记：
通过上面的操作基本可以达到不用专业的土地利用现状数据库软件而完成面积统计的工作，但是同样有一些缺陷，比如说没有考虑到区域总面积的控制，这一点可以在JSMJ中给各个图斑按比例进行核减，然后再进行DLTB面积的计算就可以达到目的。此外，耕地的田坎扣除也没有考虑，这需要制作坡度分级图等一系列的操作，在今后将进一步完善！

举报/Report

我的更多文章

下载客户端阅读体验更佳

APP专享

新浪博客

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

分享

我的更多文章

下载客户端阅读体验更佳

疯狂捕鱼