你没有关于情节的信息吗?例如。划定图的多边形区域角落的坐标?如果您有这些坐标,请将它们用作参数中的输入 poly 的 owin 。请参阅帮助文件 owin 详情。缺少任何信息,你可以尝试 ripras 估计情节的边界。
poly
owin
ripras
你现在做的是说你在矩形[0,131] [0,130]中定义一个点模式,然后你提供一堆坐标在这个区域之外的点(更大的坐标值),它们都被丢弃了。
中的数字 owin 呼叫不是窗口的宽度和高度;他们是 X 和 Y 窗口角落的坐标。
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=...) 使用多边形顶点的坐标位置。
owin(poly=...)