如何使用 R 中的 Spatstat 获取平面点模式（ppp）数据中每个点的边界框区域？答案

【问题标题】：How to get the area for a bounding box for each point in a plannar point pattern (ppp) data using Spatstat in R?如何使用 R 中的 Spatstat 获取平面点模式（ppp）数据中每个点的边界框区域？
【发布时间】：2018-12-16 14:32:38
【问题描述】：

我有一个点模式对象 (ppp)，即包含某些位置的地理坐标的数据集。例如：

> data(bei)
> df1<-as.data.frame(bei)
> View(df1)
> head(df1)
      x     y
1  11.7 151.1
2 998.9 430.5
3 980.1 433.5
4 986.5 425.8
5 944.1 415.1
6 940.5 410.4

现在，我想在 ppp 对象中的每个点上放置一个具有一定尺寸的小边界框并计算其面积。如何使用 spatstat 做到这一点？

或者我应该像素化 ppp 对象吗？那么，如何获取像素面积呢？

【问题讨论】：

我不熟悉spatstat，但也许您只想找到一个包含所有点的最小矩形，例如head(df1) 及其面积？我之所以问，是因为您的问题似乎可以在不知道 spatstat 的情况下得到回答，在这种情况下，您更有可能得到快速回答。
好吧，我不希望一个矩形包含所有点，这可以通过创建一个窗口轻松完成。我想为每个点创建矩形或框。如何在 R 中执行此操作（即使不使用 spatstat）？
neonscience.org/field-data-polygons-centroids

标签： r geospatial spatial spatstat

【解决方案1】：

在spatstat 包中，函数boundingbox 计算包含指定空间对象的最小矩形。在您的示例中，您可以键入boundingbox(bei) 以获取包含点模式数据集bei 中所有数据点的最小矩形。使用area(boundingbox(bei)) 计算面积。如果您只有一个坐标数据框而不是点模式数据集，例如df <- as.data.frame(bei)，然后boundingbox(df) 和area(boundingbox(df)) 给你同样的信息。

要搜索有关spatstat 的信息，您可以键入help(spatstat) 以获取可用命令的概述，或使用help.search 搜索关键字。或访问the spatstat.org website。

【讨论】：