如何在R中设置spatstat软件包中观测窗口几何的X和Y坐标范围

Question

嗨?我有一个苗木分布的数据，其中包括树种类型，UTM中的X和Y坐标。我想在spatstat包中的ppp()函数的帮助下，根据它们的X&Y坐标位置创建一个点阵列。我尝试了以下两种方式：

 p.patt <- ppp(mydata$X, mydata$Y)
 p.patt <- ppp(mydata$X, mydata$Y, owin(c(100,131), c(100,130)))

但是它们都有一条“警告消息: 435点被拒绝，因为它位于指定的窗口之外”。

我猜这与应该在这段代码中用c(…)指定的X和Y坐标范围有关，c(…)。我检查了X &Y和R的范围，给出了以下范围：

for X: 368615 and 368746, 
for Y: 4587355 and 4587485 

当我绘制数据图时，图的形状看起来像“倾斜的rombo”。我不知道这是不是帮助。

在这里，我随机选择了一些尝试过的数字: 100 & 131 & 130。我找不到任何关于如何在网上设置它们的信息。

所以我的问题是如何使用这些坐标范围来设置R中spatstat包中的patterm点的观察窗口几何。

非常感谢您的提前！

Answer 1

owin调用中的数字不是窗口的宽度和高度；它们是窗口角落的X和Y坐标。

由于数据点的X坐标值的范围是从368615到368746，因此窗口至少需要包含此范围。类似地，Y值的范围必须包含在窗口中。不会给出警告的最小窗口是

p.patt <- ppp(mydata$X, mydata$Y, owin(c(368615,368746), c(4587355,4587485)))

或等同于

p.patt <- ppp(mydata$X, mydata$Y, c(368615,368746), c(4587355,4587485))

但这只是可接受的最小窗口；为了进行适当的分析，您需要有关调查区域的信息。如果它不是矩形，那么，正如Ege所说的，您需要使用多边形顶点的坐标位置来指定owin(poly=...)。

Answer 2

你没有关于剧情的信息吗？例如，定界图的多边形区域的角的坐标？如果您有这些坐标，请将它们用作owin的参数poly中的输入。有关详细信息，请参阅owin的帮助文件。在缺乏任何信息的情况下，您可以尝试使用ripras来估计图的边界。

你现在要做的是，你在矩形0,131×0,130中定义了一个点模式，然后你提供了一堆坐标在这个区域之外的点(大得多的坐标值)，它们都被丢弃了。