我想简单地基于@data数据框中的相应属性值从SpatialPolygonsDataFrame对象中删除一些多边形,以便我可以绘制一个简化的/细分的shapefile。到目前为止,我还没有找到一种方法来做到这一点。
例如,假设我要从此世界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
# 1 AG DZ DZA 12 Algeria 238174 32854159 2 15 2.632 28.163
# 2 AJ AZ AZE 31 Azerbaijan 8260 8352021 142 145 47.395 40.430
# 3 AL AL ALB 8 Albania 2740 3153731 150 39 20.068 41.143
# 4 AM AM ARM 51 Armenia 2820 3017661 142 145 44.563 40.534
# 5 AO AO AGO 24 Angola 124670 16095214 2 17 17.544 -12.296
如果我做这样的事情,该图将不会反映任何变化。
world.map@data = world.map@data[world.map@data$AREA > 30000,]
plot(world.map)
如果执行此操作,将得到相同的结果:
world.map@data = world.map@data[world.map@data$NAME != "Antarctica",]
plot(world.map)
任何帮助表示赞赏!