使用包win combine maptoo处理聚合美国人口普查局空间数据.pdfVIP

结合空间数据*

RogerBivand

1引言

2检查拓扑结构

在本指南中，我们通过一个实际例子来探讨如何清理从人口普查局发布的shapefile文件读

入R的空间对象。然后，我们使用人口普查局的文本表将这些对象聚合到大都市区。

在此案例中，数据表示的是1990年北卡罗来纳州、南卡罗来纳州和弗吉尼亚州

的县界多边形，如图1所示。每个多边形的属性数据包括多边形标识符、州和

县标识符以及县名。所有空间对象都具有相同数量的属性数据列，且类型和名称相

同。文件以shapefile格式，不包含坐标参考系统；元数据用于选择CRS值。

owd-getwd()setwd(system.file(shapes,package=

maptools))library(maptools)nc90-readShapeSpatial(

co37_d90)proj4string(nc90)-CRS(+proj=longlat+datum=

NAD27)sc90-readShapeSpatial(co45_d90)

proj4string(sc90)-CRS(+proj=longlat+datum