Tag: maptools

简单的方法通过R中的属性来对SpatialPolygonsDataFrame进行子集(即删除多边形)

我想简单地从SpatialPolygonsDataFrame对象中删除一些基于@data数据框中的相应属性值的多边形,这样我就可以绘制一个简化/子集的shapefile。 到目前为止,我还没有find办法做到这一点。 例如,假设我想从这个世界shapefile中删除面积小于30000的所有多边形。我该如何去做这件事? 或者,同样,我怎样才能删除Antartica? require(maptools) getinfo.shape("TM_WORLD_BORDERS_SIMPL-0.3.shp") # Shapefile type: Polygon, (5), # of Shapes: 246 world.map <- readShapeSpatial("TM_WORLD_BORDERS_SIMPL-0.3.shp") class(world.map) # [1] "SpatialPolygonsDataFrame" # attr(,"package") # [1] "sp" head(world.map@data) # FIPS ISO2 ISO3 UN NAME AREA POP2005 REGION SUBREGION LON LAT # 0 AC AG ATG 28 Antigua and Barbuda 44 83039 19 29 -61.783 17.078 […]