小编Dan*_*son的帖子

我已经有一段时间了,并取得了一些成功,但是当涉及到重写文件时,我没有.

我正在使用的shapefile是一个地形和都市区域的多边形shapefile,但我对大都市区不感兴趣所以我正在使用R中的shp来尝试从我的地图中删除它们.

数据源

要下载正确的文件,请确保在"Metropolitan/Micropolitan Statistical Area(2010)"下选择"一个国家/地区中的所有州"

这是我到目前为止所拥有的:

library(maptools)
met=readShapeSpatial("tl_2010_us_cbsa10.shp")
met=met@data

子集之前:

NAMELSAD10      LSAD10
Anchorage, AK Metro Area   -  M1 
Clarksdale, MS Micro Area    - M2
Richmond, VA Metro Area   -  M1
Big Spring, TX Micro Area  -   M2
Dallas-Fort Worth-Arlington, TX Metro Area  -   M1
Rio Grande City-Roma, TX Micro Area  -   M2

然后:

submet=subset(met, LSAD10 == "M1")

子集后:

NAMELSAD10 LSAD10
Anchorage, AK Metro Area   -  M1
Richmond, VA Metro Area  -   M1
Dallas-Fort Worth-Arlington, TX Metro Area   -  M1 …