在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)

但是，您將收到一條警告，警告您的點由於在窗口之外而被拒絕，並且該對象將不包含任何點。