为什么网络上的点模式尺寸很大 - spatstat答案

【问题标题】：why large size of point pattern on network - spatstat为什么网络上的点模式尺寸很大 - spatstat
【发布时间】：2017-03-10 04:41:47
【问题描述】：

我正在尝试从 Baddeley 等人的书中学习在 R 中使用 spatstat 包。我已将我的 shapefile 转换为 psp 对象（这里说丢弃了 6 列数据），然后我使用 as.linnet 将其转换为适合 spatstat 的线性网络。（这里说网络未连接）

然后我想使用rpoislpp 在这个网络上创建一个点模式，我使用了 2 的泊松强度。

> abc<-rpoislpp(2,final_roads)
> head(abc)
Point pattern on linear network
3 points
Linear network with 4296 vertices and 4475 lines
Enclosing window: rectangle = [30093.5, 278045.11] x [308520.5, 606556.7] units
> abc
Point pattern on linear network
10190733 points
Linear network with 4296 vertices and 4475 lines
Enclosing window: rectangle = [30093.5, 278045.11] x [308520.5, 606556.7] units

这创建了一个巨大的 1GB 文件。我不明白怎么做。我从书中的理解是泊松强度是线性网络总长度的点数。我的网络大约5000公里（一个省的路网），所以我只是在整个网络中请求2个点。

我的问题是

我不明白怎么可能是1000万点。我只是想在整个网络中创建 2 个点并绘制它们。
如何绘制这些点并获得它们的 x,y

我做错了吗？任何建议都会很棒。谢谢

【问题讨论】：

我意识到问题1的错误。但我仍在寻找问题2的答案。

标签： spatstat

【解决方案1】：

Q1：在泊松过程中，强度 lambda 是每单位长度的预期点数。如果网络的总长度为L，那么强度为 lambda 的泊松过程将包含大约L * lambda 个点。在您的示例中，打印输出表明坐标以米为单位，因此您有L = 5 000 000 metres，如果lambda = 2，您将获得大约 1000 万个点。

要获得固定数量的随机点，请使用runiflpp 而不是rpoislpp。所以runiflpp(2, final_roads) 会在你的道路网络上放置 2 个随机点。

最佳实践：使用unitname 指定长度单位，并考虑使用rescale 将坐标从米重新调整为公里。

Q2：如果X 是线性网络上的点模式，那么要在线上绘制点，只需使用plot(X)。要仅提取点位置，请使用Y <- as.ppp(X)，然后仅绘制点位置，plot(Y)。要提取x,y 坐标，请使用xy <- coords(Y) 或xy <- coords(X, local=FALSE)。

更多详情请参阅the spatstat book 的第 17 章。

【讨论】：

谢谢阿德里安。我意识到我的错误，也感谢你的 runiflpp 提示。我使用了 write csv，它也可以给出点的坐标。我会尝试你的建议。 PS。我这周买了你的书，我很喜欢。
请看，我不明白为什么它说网络未连接（问题中的粗体文字）。我寻找悬挂物，但没有。我应该纠正什么吗？

【解决方案2】：

我无权访问你的数据，所以我以内置数据集chicago为例：

abc <- unmark(chicago)
plot(abc) # Plot the network and the points on the network
X <- as.ppp(abc) # Convert to planar point pattern without network
plot(X) # Plot the points in the plane
co <- coords(X) # Extract the Cartesian coordinates

这就是你要找的吗？也可以直接使用coords(abc)，得到笛卡尔坐标和网络中的局部坐标。

【讨论】：