如何在 R 中的 spatstat 包中设置坐标 X 和 Y 的范围以观察 window 几何
how to set the ranges of coordinates X & Y for observation window geometry in spatstat package in R
嗨？我有一个幼苗分布数据，其中包含 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 个点被拒绝，因为它们位于指定的 window 之外”。
我猜这与这段代码中c(…), c(…)中应该指定的X和Y坐标范围有关。我检查了 X & Y 的范围，R 给了我以下范围：
for X: 368615 and 368746,
for Y: 4587355 and 4587485
当我绘制数据时，绘图的形状看起来像 "tilted rombo"。不知道有没有帮助
这里我随便选了几个号码试试：100 & 131 & 130，网上怎么设置都找不到。
所以我的问题是如何使用这些坐标范围来设置 R 中 spatstat 包中点模式的观察 window 几何？。
非常感谢您！
你没有剧情信息吗？例如。划定地块的多边形区域角的坐标？如果您有这些坐标，请将它们用作 owin 的参数 poly 中的输入。有关详细信息，请参阅 owin 的帮助文件。在缺少任何信息的情况下，您可以尝试 ripras 来估计地块的边界。
你现在做的是说你在矩形 [0,131]×[0,130] 中定义了一个点模式，然后你提供了一堆坐标在这个区域之外的点（更大的坐标值）并且它们都被丢弃了。
owin调用中的数字不是window的宽高；它们是 window.
角的 X 和 Y 坐标
由于数据点的X坐标值范围是368615到368746，所以window至少要包含这个范围。同样，Y 值的范围必须包含在 window 中。不会发出警告的最小 window 是
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))
但这只是可以接受的最小值window；为了进行正确的分析，您需要有关调查区域的信息。如果它不是矩形，那么正如 Ege 所说，您需要使用多边形顶点的坐标位置指定 owin(poly=...)。
嗨？我有一个幼苗分布数据，其中包含 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 个点被拒绝，因为它们位于指定的 window 之外”。
我猜这与这段代码中c(…), c(…)中应该指定的X和Y坐标范围有关。我检查了 X & Y 的范围，R 给了我以下范围：
for X: 368615 and 368746,
for Y: 4587355 and 4587485
当我绘制数据时，绘图的形状看起来像 "tilted rombo"。不知道有没有帮助
这里我随便选了几个号码试试：100 & 131 & 130，网上怎么设置都找不到。
所以我的问题是如何使用这些坐标范围来设置 R 中 spatstat 包中点模式的观察 window 几何？。
非常感谢您！
你没有剧情信息吗？例如。划定地块的多边形区域角的坐标？如果您有这些坐标，请将它们用作 owin 的参数 poly 中的输入。有关详细信息，请参阅 owin 的帮助文件。在缺少任何信息的情况下，您可以尝试 ripras 来估计地块的边界。
你现在做的是说你在矩形 [0,131]×[0,130] 中定义了一个点模式，然后你提供了一堆坐标在这个区域之外的点（更大的坐标值）并且它们都被丢弃了。
owin调用中的数字不是window的宽高；它们是 window.
X 和 Y 坐标
由于数据点的X坐标值范围是368615到368746，所以window至少要包含这个范围。同样，Y 值的范围必须包含在 window 中。不会发出警告的最小 window 是
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))
但这只是可以接受的最小值window；为了进行正确的分析，您需要有关调查区域的信息。如果它不是矩形，那么正如 Ege 所说，您需要使用多边形顶点的坐标位置指定 owin(poly=...)。