问题标签 [spatstat]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
2 回答
84 浏览

r - 如何使用 R 在图像中获取颜色编码的绘图区域?

嗨,世界的 R 专家,假设我有一个点模式可以生成强度图,并且该图在像素化图像的 3 个区域中进行了颜色编码......我怎样才能获得颜色编码的区域?

这是一个使用 spatstat 的示例:

提前谢谢你

0 投票
2 回答
1283 浏览

r - 将“ppp”转换为多类型

我一直在通过一系列 spatstat 函数运行两个未标记的平面点模式数据集。现在我想用Kcross.inhom函数来描述两者之间的交互,但是Kcross只对标记的数据起作用,所以我将所有xy数据合并到一个csv文件中,并添加了一个区分两者的列。我已经建立了以下点模式对象,但不明白如何为我的目的编辑 Kcross 的后续示例。或者,也许有更好的方法?谢谢你的帮助!

0 投票
3 回答
743 浏览

r - 在 spatstat 中通过 ID 测量点之间的距离

我正在尝试测量雏鸟从原来的巢址移动到他们第二年建立的领地的距离。我想使用 spatstat 包来执行此操作,但我在输出的某些功能(例如,nncross、nndist、crossdist)时遇到了问题。这是我的代码的样子,数据在这里:https ://sites.google.com/site/datastackoverflow/shapefiles

结果如下所示:

是否可以保留点 ID,在这种情况下是 nestID 和相应的birdID,而不是行号和列号?

0 投票
2 回答
403 浏览

r - 在比较 R 中的 2 个不同数据集时找到具有相同值的点的最近邻

我有 2 个数据框(df1 和 df2),由三列组成;x 坐标,y 坐标,类别(有 5 个级别 AE)。所以我基本上有两组点数据,每个点都分配给一个类别

例如

等等......(虽然我的两个数据框都有100个点)

我想从第二个数据帧(df2)中为我的第一个数据帧(df1)中的每个点找到相同类别的最近邻居。

我在包 spatstat 中使用了 nncross 来为 df1 中的每个点与 df2 找到最近的邻居,然后列出这些距离中的每一个,如下所示;

但是我正在努力弄清楚如何使用存储在 ppp 对象中的类别标签(由标记定义)来从同一类别中找到最近的邻居。我相信它应该是相当直截了当的,但如果有人有任何建议或任何替代方法来达到相同的结果,我将不胜感激。

0 投票
2 回答
384 浏览

r - ggmap 上的 SpatStat 结果图

我想知道是否可以在从 ggmap 创建的绘图顶部添加从 spatstat ( https://cran.r-project.org/web/packages/spatstat/index.html )计算的结果。

例如,如果我计算一个 quadratcount,我如何将结果添加到 ggmap 的顶部。?

谢谢你的帮助,

阿尔诺

0 投票
1 回答
937 浏览

r - 在 R{Spatstat} 中将 IDW 插值限制为多边形

我是一个相当新的 R 用户,使用嵌入在美国都市区域内的地方的空间数据。我正在为许多变量生成 IDW 地图,但在定义 IDW 的边界时遇到了麻烦。我只想在都市区边界范围内插入数据,而我目前的方法填充了整个图。

我首先ppp从 a 创建一个对象dataframe

使用投影此ppp对象plot(x.ppp)会显示每个变量的图(上面代码中的“v1”和“v2”)。但是,绘图填充了整个矩形框。在绘制对象时,我ppp通过为其分配owin对象来限制此框架,如下所示:

现在使用它来绘制这个plot(x.ppp)图,可以让我在他们的都市区的多边形内正确地绘制出我的数据点图。虽然这是一个好的开始,但我想继续将此空间框架传递到 IDW 地图上。查阅 spatstat 手册让我相信我可以将我的owin对象作为as.mask参数传递,但是每当我这样做时,都会出错。示例代码如下:

这给了我错误“开关错误(在,像素= {:EXPR必须是长度为1的向量”。使用at="points"不起作用,因为我得到的是结果dataframe对象而不是imlist。如果我省略了as.mask参数,我会得到一个正确的 idw,但它插值超出了我之前能够使用该owin对象设置的都市区边界。

我认为这是一项相当常见的任务,但查阅 spatstat 文档让我有点不知所措。这个as.mask论点似乎可行,我不确定我哪里出错了。关于如何解决这个问题的任何意见或建议,或者我可以求助于尝试更好地理解这个问题的来源,都会有很大的帮助。谢谢!

0 投票
2 回答
350 浏览

r - R中镶嵌分割链的包,spatstat包?

我正在尝试创建漂亮的聚集点图形。是否有一个包可以在点的镶嵌之间创建分界链?理想情况下,它适合在ggplot.

这是一些示例代码:

这输出如下:http: //imgur.com/CCXeOEB

没有划分的点簇

我正在寻找的是这个:http: //imgur.com/7nmtXjo

具有分隔链的点簇

我知道存在一种算法。

有什么想法/替代方案吗?

0 投票
3 回答
1216 浏览

r - spatstat 中的单位长度

我对 Spatstat 中的 KEST 函数有一个非常简单的问题。KEST 图形输出我正在使用 Spatstat 中的 KEST 函数来评估数据集中的空间随机性。我已经上传了分布在伦敦的 lat 和 long 值,并将它们转换为 PPP 对象,使用 ripras 函数来指定空间域。当我在我的 ppp 上运行我的 KEST 分析并绘制图表时,我最终在 x 上得到了一个 r 值,但是虽然我知道这是一个距离测量,但我不知道它使用的是什么单位。我得到这个摘要输出:

x 轴上的最大 r 为 0.1 个单位,y 轴上的 K(r) 为 0.04。我如何弄清楚这些等于什么距离单位?

0 投票
1 回答
518 浏览

r - R 和 spatstat 的线密度

我使用 R 和 spatstat 包使用 pixellate 函数 (spatstat) 创建血管轨迹的密度图。我能够创建地图并绘制轨迹密度。

在此处输入图像描述

现在我想改进我的地图设计并做一些类似 Niels Willems ( Google Scholar ) 创建的地图的事情。下图举例说明了我想要实现的目标。你知道我可以做些什么来获得这种 3D 效果吗?

在此处输入图像描述

谢谢你的帮助,

阿尔诺

0 投票
1 回答
236 浏览

r - 最近邻标记的平均值(空间点模式)

我试图获得每个点中标记的平均值,计算每个点周围一定距离内的所有可用点,包括我们测量的点。

我已经将 markmean() {spatstats} 与缓冲区一起使用,但我不确定它是否正在做我想做的事情。

这是我正在尝试做的一个简单示例:

空间点看起来像这样:

带有 ID 和标记的点分布(标记在点 13=60,在点 14=40,其余点标记为 0 当然我的真实数据不遵循规则模式,所以点之间的距离是可变的.

我希望对于 sigma 1 的点 13,markmean 是:(60+40+0+0+0)/5= 20 或者对于点 7,markmean 是:(60+0+0+0 +0)/5=12。

但是按照上面的代码,我得到:

你知道我为什么会得到这些结果吗?您将如何正确计算某个缓冲区内每个点周围标记的平均值?

非常感谢您的时间和帮助!