使用ggmap绘制形状文件边界内随机分布的点

Question

是否可以使用 ggplot2 / ggmap 在 shapefile 定义的给定空间区域内随机绘制 geom_points？

我考虑过 geom_jitter，但是我需要随机分布的地块，而不是跨越空间边界。

无耻地从@matthiash 那里借来的示例数据here。

library(rgdal)
library(ggmap)

# Get shapefile with Drammen municipality borders
tmpzip<-tempfile()
tmpdir<-tempfile()
dir.create(tmpdir)
download.file("http://www.kartverket.no/Documents/Kart/N50-N5000%20Kartdata/33_N5000_shape.zip",tmpzip)
unzip(tmpzip, exdir=tmpdir)
kommune <- readOGR(dsn=tmpdir, layer="NO_AdminOmrader_pol")
kommune<-kommune[kommune$NAVN=="Drammen",]
kommune<-spTransform(kommune, CRS("+init=epsg:4326"))
dat<-fortify(kommune)

#get the base map
map <- get_map(location = "Drammen",
           maptype = "watercolor", source = "stamen", zoom = 11)

下面的代码从绘制在顶部的 shapefile 中绘制了区域 ID 为 154 的底图。

ggmap(map, extent = "normal", maprange = TRUE)+
 geom_polygon(data = dat,
           aes(long, lat, group = group),
           fill = "orange", colour = "red", alpha = 0.2)

我想做的是在 dat$id==154

定义的 shapefile 区域内随机绘制 10 个点

Answer 1

好的，我整理好了。解决方案在包 "raster".

中的 spsample()

d<-data.frame(id=NA,x=NA,y=NA)
l<-data.frame(id=154,n=10)

for (i in unique(l$id)){
 temp<-spsample(kommune[which(kommune$OBJECTID==i),],n=l[l$id==i,"n"],type="random")
 temp<-as.data.frame(temp)
 temp$id<-i
 d<-rbind(d,temp[,c("id","x","y")])
}
d<-d[-1,] #drop the first empty row

ggmap(map, extent = "normal", maprange = T)+
  geom_polygon(data = dat,
           aes(long, lat, group = group),
           fill = "blue", colour = "yellow", alpha = 0.1)+
  geom_point(aes(x = x, y = y), data = d[which(d$id==154),], alpha = .9,show.legend = T)

使用ggmap绘制形状文件边界内随机分布的点

using ggmap to plot points randomly distributed within boundaries of shape file

gis

r

shapefile

ggplot2

ggmap