合并级别，使得窗口中的 n 等于级别数

问题描述

答案确实解决了问题，但情节非常混乱，如下所示。

我想要的是一个大窗口，所有这些点都在那个窗口内。但是正如您所看到的，由于 131 个级别，该图显示每个点都有自己的窗口。我认为这可以通过合并 n = 131 中的 window 级别来解决。作为参考，sf_object 如下所示：

Simple feature collection with 131 features and 3 fields
geometry type:  MULTIPOINT
dimension:      XY
bBox:           xmin: -1932934 ymin: 4958872 xmax: -1439558 ymax: 5861173
projected CRS:  NAD83(2011) / UTM zone 16N
# A tibble: 131 x 4
# Groups:   COOPID [131]
   COOPID STATION_NA                         geometry Annual_Precipitation
 *  <dbl> <chr>                      <MULTIPOINT [m]>                <dbl>
 1      0 Ontario                ((-1899685 5335073))                 9.24
 2 100010 ABERDEEN EXPERIMNT STN ((-1610453 5091311))                12.4 
 3 100227 AMERICAN FALLS 3 NW    ((-1623401 5075011))                20.4 
 4 100282 ANDERSON dam           ((-1807106 5212322))                16.3 
 5 100347 ARBON 2 NW             ((-1606302 5034484))                10.2 
 6 100375 ARCO                   ((-1622855 5179969))                19.5 
 7 100448 ARROWROCK dam          ((-1834338 5254236))                20.1 
 8 100470 ASHTON                 ((-1458491 5179214))                37.5 
 9 100528 AVERY RS #2            ((-1678382 5654084))                25.3 
10 100667 BAYVIEW MODEL BASIN    ((-1691954 5753129))                 9.69

可以在下图中看到带有类型的 ppp 对象。如何将级别合并到一个大窗口中，使得 n 中的 window = 131（所有站点）和 x = 列表中的所有坐标和 y 相同？我认为这将解决在一个大窗口（整个感兴趣区域）中绘制 131 个点的问题。我希望这是有道理的。

解决方法

通用 split 操作将数据分成几组数据。请参阅help(split)。

split.ppp 包中的方法 spatstat 将一个点模式划分为一个点模式列表。例如，如果您有一个点模式 X 表示森林中树木的位置，并带有指示每棵树的种类的标记，那么 split(X) 将生成一个点模式列表，每个点模式都包含一种树种。

您似乎不想像您在前三个问题中提出的那样拆分您的点模式。相反，您可能想要生成具有因子值标记的单点模式，其中标记值标识分组/分类。您已经在上一个问题中通过执行 marks(X) <- factor(.....) 完成了此操作，因此只需绘制该点模式 X。

r r spatstat