所以我有shapefile (在页面底部名为波瓦蒂 )。我想用我的和弦填充特定的多边形。
比方说,我想填写整个地区: lat:52.599427 lon:20.7572137
我知道如何在地图上设置点和弦,但如何在地图中填充整个多边形?
发布于 2015-02-06 16:08:05
该shapefile中有重复的区域名称,因此您必须使用多边形数字id来填充:
library(rgdal)
library(rgeos)
library(ggplot2)
pow <- readOGR("POWIATY.shp", "POWIATY")
plot(pow)
where <- over(SpatialPoints(cbind(20.7572137, 52.599427)), pow, TRUE)
reg <- data.frame(id=rownames(where[[1]]))
map <- fortify(pow)
gg <- ggplot()
gg <- gg + geom_map(map=map, data=map,
aes(x=long, y=lat, map_id=id),
fill="white", color="black", size=0.25)
gg <- gg + geom_map(data=reg, map=map,
aes(fill=id, map_id=id), color="steelblue")
gg <- gg + coord_map()
gg <- gg + theme_bw()
gg

https://stackoverflow.com/questions/28368279
复制相似问题