合并级别，使 window 中的 n 等于级别数

Question

关于这个

答案确实解决了问题，但是如下图情节很乱。我想要的是一个大 window，里面有所有这些点 window。但正如您所见，由于有 131 个级别，该图显示每个点都有自己的 window。我认为这可以通过合并 n = 131 in 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对象可以在下图中看到。我如何将级别合并为一个大的 window 使得 n 在 window = 131（所有站点点）和 x = 列表中的所有坐标以及 y?我认为这将解决在一个大 window（整个感兴趣区域）中绘制 131 个点的问题。我希望这是有道理的。

Answer 1

泛型split操作将数据分成几组数据。参见 help(split)。

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

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

合并级别，使 window 中的 n 等于级别数

Merging levels such that n in window is equal to the number of levels

r

spatstat