【问题标题】:Merging levels such that n in window is equal to the number of levels合并级别,使得窗口中的 n 等于级别数
【发布时间】:2021-03-03 04:01:34
【问题描述】:

关于这个question

答案确实解决了问题,但情节非常混乱,如下所示。 我想要的是一个单一的大窗口,所有这些点都在那个窗口内。但正如您所见,由于 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 对象的类型如下图所示。如何将级别合并到一个大窗口中,使window = 131 中的n(所有站点点)和x = 它们在列表中的所有坐标并且对于y 相同?我认为这将解决在一个大窗口(整个感兴趣区域)中绘制 131 个点的问题。我希望这是有道理的。

【问题讨论】:

  • 您发布了几个相关的问题,我认为是时候停下来提出一个好的问题,您可以:1. 解释您拥有的数据,2. 您想要实现的目标3. 显示你到目前为止使用的代码行。不仅仅是这里和那里的输出截图。如果您无法共享您的数据,请使用 spatstat 中的内置数据集,例如 lansing(6 种类型中的每一种都有 100 个点的多类型)、finpines(每个点有两个不同的数字标记)或您在vignette('datasets')。请考虑使用reprex 包来解决您的问题。

标签: r spatstat


【解决方案1】:

通用的split 操作将数据分成几组数据。见help(split)

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

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

【讨论】:

  • 它成功了,非常感谢你如此详尽的解释。
猜你喜欢
  • 2021-12-15
  • 2018-04-20
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 2013-11-08
  • 1970-01-01
  • 2018-03-02
  • 2011-01-03
相关资源
最近更新 更多