在spatstat上绘制的不规则多边形区域

Question

这是我第一次使用spatstat软件包，所以我需要一些建议。 我试图将坐标数据绘制到不规则的多边形区域（格式为.shp），以计算空间分析（如Ripley的K）。如何添加不规则的多边形区域作为图？ 如何将.ppp数据从坐标合并到多边形区域？ 我使用了以下代码：

将坐标数据转换为.ppp格式

library(spatstat)
library(sp)
library(maptools)

tree.simu <- read.table("simulation.txt", h=T)
tree.simu.ppp <-ppp(x=tree.simu$X,y=tree.simu$Y,window=owin(c(min(tree.simu$X),max(tree.simu$X)),c(min(tree.simu$Y),max(tree.simu$Y))))
plot(tree.simu.ppp)

通过此功能，我正在考虑将绘图区域作为坐标的最大和最小值。 我想将多边形边界作为图。

绘制不规则多边形区域

area <- readShapePoly("Area/Fragment.shp")
plot(area)
plot(tree.simu.ppp, add=T)

要么

points(tree.simu.ppp)

该软件包接受最后一个功能，但是当我尝试将两个文件一起绘制时，.shp文件似乎充满了整个区域。 我无法可视化坐标数据。

谢谢，非常感谢您的帮助！

附注：如果您对这些问题了解任何材料，请很高兴看看

Answer 1

实际上，这是由于@jlhoward的评论中推测的边界框不一致。 您的点在[273663.9，275091.45] x [7718635，7719267]中，而多边形在[-41.17483，-41.15588] x [-20.619647，-20.610134]中包含。

假设坐标确实与窗口一致，则将其放入ppp对象的正确方法是：

library(spatstat)
library(sp)
library(maptools)

area <- readShapePoly("Area/Fragment.shp")
area <- as(area, "owin")
tree.simu <- read.table("simulation.txt", h=T)
tree.simu.ppp <-ppp(x=tree.simu$X,y=tree.simu$Y,window=area)

但是，您将收到一条警告，警告您的点由于在窗口之外而被拒绝，并且该对象将不包含任何点。