如何为每个数据集条目在相应的多边形中创建随机点

Question

我创建了一个多边形图层（Poly.buffer），还有一个观察的虚拟数据集（df），它有一个带有这个多边形图层的公共字段（即SiteName）。

我现在想为相应的多边形内的每个数据集条目生成一个随机点。

目标是有138个随机点，它们位于相应的多边形中，并具有数据列表的所有属性。

以下是我尝试过的内容（基于此https://gis.stackexchange.com/questions/74250/how-to-create-random-point-in-according-polygon-for-each-dataset-entry-in-a-list）：

##load packages
library(ggplot2)
library(sp)
library(rgdal)
library(rgeos)

##First create dummy dataset
SiteName=c("Site 1", "Site 2", "Site 3", "Site 4", "Site 5")
Latitude=c("-16.40407","-16.41253","-16.39207","-16.40035","-16.40457")
Longitude=c("145.8157","145.7117","145.8308","145.8368","145.8405")

Poly=data.frame(SiteName,Latitude, Longitude)
Poly$Latitude=as.numeric(as.character(Poly$Latitude))
Poly$Longitude=as.numeric(as.character(Poly$Longitude))

Site1=data.frame(SiteName="Site 1",value=1:10)
Site2=data.frame(SiteName="Site 2",value=1:3)
Site3=data.frame(SiteName="Site 3",value=1:64)
Site4=data.frame(SiteName="Site 4",value=1:13)
Site5=data.frame(SiteName="Site 5",value=1:18)
Site1x=data.frame(SiteName="Site 1",value=1:30)
df=dplyr::bind_rows(Site1,Site2,Site3,Site4,Site5,Site1x)

##Create Poly.buffer##
coordinates(Poly)=~Longitude+Latitude
proj4string(Poly)=CRS("+proj=longlat +datum=WGS84")
Poly.utm <- spTransform(Poly, CRS("+proj=utm +zone=55 +south +units=m 
+ellps=WGS84"))
Poly.buffer=gBuffer(Poly.utm,width=400,byid=TRUE)
Poly.buffer=spTransform(Poly.buffer, CRS("+proj=longlat +datum=WGS84"))

### Preparing the SpatialPointsDataFrame
spdf <- matrix(as.numeric(NA), nlevels(Poly.buffer$SiteName), 1)
spdf <- as.list(spdf)

### Sample the coordinate, match it with data in spdf.
### sample(spsample()) fix the size of the sample
for (i in seq(Poly.buffer$SiteName))
spdf[i] <- SpatialPointsDataFrame(
           sample(spsample(Poly.buffer[order(Poly.buffer$SiteName)==i,], n = 200, "random"),
table(df$SiteName)[[i]]),
           df[df$SiteName==dimnames(table(df$SiteName))[[1]][i],],
           proj4string=poly.crs, 
           match.ID=FALSE)

## Merging together the list to make on SpatialDataFrame
do.call("rbind", spdf) -> spdf

我可以运行此代码的唯一方法是选择match.ID为FALSE。但这样做会为每个数据集条目生成与相应多边形不匹配的坐标。

见下图：绘制它时的样子：

df.spdf=as.data.frame(spdf)
fPoly.buffer=fortify(Poly.buffer)
ggplot(data=df.spdf, aes(x=x,y=y)) +  geom_polygon(data=fPoly.buffer, 
aes(long, lat, group = group)) + geom_point(size=0.1,aes(color=SiteName))

当我尝试使用match.ID = TRUE运行代码时，我收到以下错误消息：

  

错误   SpatialPointsDataFrame（sample（spsample（Poly [order（Poly $ SiteName）==：   数据和坐标的row.names不匹配

我不确定如何取得进展。任何建议都将不胜感激。

见下文，看起来应该是什么样子。当我从虚拟df中排除Site1x时，代码工作正常。但是我需要它来处理它。

Answer 1

所以问题在于，当我执行do.call函数时，它根据行号合并了不同的列表组件。我通过简单地重新排序我的数据集并重新分配新的行号来解决这个问题。

df=df[order(df$SiteName),]
rownames(df) <- 1:nrow(df)