问题标签 [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.
r - 如何通过在 spatstat 中锚定中间点来旋转 psp 对象?
我有一个线段图案,它与另一条线成 90 度角:
这是一个视觉表示:
现在我想旋转lines
,使它们main
在同一点交叉,但角度为 45 度。
这可行,但angles.psp
函数似乎owin
根据我的需要单独旋转窗口()而不是行。
有没有办法lines
相对于main
线旋转所有角度,使角度为 45 度,但交叉点保持不变?
r - 如何在同一设备中绘制多个 'image' 和 'plot.default' 图形?
我有一个函数,它使用spatstat
'scolourmap
函数为绘图生成色标。该colourmap
函数使用该image
函数绘制色标。这是一个例子:
我想并排绘制其中几个图,但我的尝试失败了:
我找到了一个胎面,Paul Murrell 博士说那image.plot
是不相容的layout
。不过,我相信这个问题是可以解决的。如何改进代码以使用layout
或进行绘图par
?
我想做类似以下的事情,但要分别包含每个图的色阶:
编辑: @IShouldBuyABoat 建议调查gridBase
. 混合基本图形和 ggplot 图形有一个很好的 SO 答案,但我找不到一个有效的grid
函数组合来制作绘图。这是其中一项试验:
r - spatstat 超帧中的意外行为?
我正在尝试hyperframe
从包含owin
对象(或任何 spatstat 对象)的列表超帧中制作 spatstat ,但是在使用rbind.hyperframe
命令 indo.call
将列表绑定到更大的列表时会遇到意外行为hyperframe
:
此时,列表可以正常工作。owin
对象分别位于列出的超帧的每个单元格中:
接下来我对超帧使用do.call
命令:rbind
现在圆圈是列表(它们本质上是在 spatstat 中的),但似乎所有这些都绑定在一起:
注意一个单元格的三重名称。
我的代码依赖于列表,我非常希望能够将这些列表绑定到一个更大的hyperframe
. 有没有办法解决这种行为,即如何在不绑定“圆形”列中的单元格的情况下绑定超帧?
r - 在 R{spatstat} 中查找点之间的欧几里德距离,由不规则多边形窗口限制
我试图找到由不规则多边形限制的两点之间的欧几里德距离。(即,必须将距离计算为通过给定窗口的路线)
这是一个可重现的示例:
但是,我们仅通过绘制点就知道
路线被限制在多边形内时的距离应该更大(在本例中为 5.00)。
在 {spatstat} 中是否有另一个我不知道的功能可以做到这一点?或者有人对另一个可以做到这一点的包有任何其他建议吗?
我试图找出水体中两点之间的距离,所以我的实际数据中的不规则多边形更复杂。
任何帮助是极大的赞赏!
干杯
r - 将 R 包 spatstat 函数应用于一维数据?
spatstat 中的函数主要用于 2-3-dim 数据分析。是否有可能将它们应用于一维数据?
- 2-dim 中的类 ppp 有巨大的能力。
- 对于任意维度,有一个非常通用的类 ppx - 但这就是问题所在 - 只有很少的函数可用。
- 我可以用大锤来破解将一维数据膨胀到二维数据并最终投射回一维数据的问题吗?
- 或者我应该更好地重写一维函数(rpoispp,rmpoispp,...)?
r - Voronoi diagram polygons enclosed in geographic borders
I am trying to create Voronoi polygons (aka Dirichlet tessellations or Thiessen polygons) within a fixed geographic region for a set of points. However, I am having trouble finding a method in R that will bound the polygons within the map borders. My main goal is to get accurate area calculations (not simply to produce a visual plot). For example, the following visually communicates what I'm trying to achieve:
which produces the following:
Conceptually I want to intersect counties
with vt
which should provide a set of polygons bounded by the county borders and accurate area calculations for each. Right now, vt$summary
provides area calculations for each polygon, but they are obviously overstated for all but the one interior polygon, and deldir()
appears to only accept rectangular enclosings for its rw
argument. I am new to R's geospacial capabilities, so am open to other approaches beyond what I outlined above.
r - 时空聚类 :: 不同空间中的点模式
按照Thomas 的这些说明,我使用 R 中的spatstat包创建了一个 ppp 对象。
因为我的数据集包含时间维度,所以我想扩展此分析以考虑 3D 时空进行聚类分析。
SPATSTAT VERSION 2 Preliminary Announcement 6th revision,日期为 2010 年 2 月 9 日,表明这应该可以通过新的 R 包spatstat2 实现。公告称:
点模式的新类 pp 将支持:
- 在每个点上附加多个标记
(例如,树木可以通过它们的种类、胸径和叶片毒性测定来标记)。- 不同空间中的点模式,
包括 3D 点、2D 时空、3D 时空、2D 空间中的 1D 网络、不同的距离度量- 任何类型的
标记(每个标记可以是窗口、点图案、功能等)
但是,除了公告之外,我找不到有关预期的spatstat2包的许多其他详细信息。而且我也无法从 CRAN 安装软件包:
我的问题是:
R 中可用于事件数据聚类分析的最佳方法是什么(每个事件都有地理空间坐标 X、Y 和时间戳)?
spatstat2包发生了什么?是可用的么 ?
谢谢您的帮助!
基思·赫尔弗里希
r - Voronoi 镶嵌的共享边表
我正在尝试根据库dirichlet()
函数生成的 Voronoi 镶嵌(又名 Dirichlet 镶嵌或泰森多边形)创建一个多边形邻居表。spatstat
例如,在下图中,右上角和右下角的瓷砖各有 2 个邻居,中间右侧的瓷砖有 4 个邻居,其余两个瓷砖各有 3 个邻居。我想捕获表中的邻居对,并理想地捕获它们共享的边界线的长度:例如,'Tile1'、'Tile2'、'shared_edge_length'。
最初,我尝试使用intersect.tess()
、intersect.own()
和polyclip
函数循环并比较镶嵌中的每一对多边形,但我猜这些不起作用,因为尽管共享边缘,瓷砖根据定义在区域上不重叠。是否有一个简单的功能来实现这一点(替代方法可能是遍历这些$bdry
点)?好像regeos
包有gTouches
,但我找不到类似的东西spatstat
。
这是我目前的非工作方法:
r - 在 spatstat 上绘制的不规则多边形区域
这是我第一次使用 spatstat 包,所以我想要一些建议。我正在尝试将坐标数据绘制到不规则多边形区域(格式 .shp)中,以计算像 Ripley's K 这样的空间分析。如何将不规则多边形区域添加为绘图?如何将坐标中的 .ppp 数据合并到多边形区域中?我使用了以下代码:
将坐标数据转换为 .ppp 格式
使用此功能,我将绘图区域视为坐标的最大值和最小值。我想将多边形边界作为绘图。
绘制不规则多边形区域
或者
该包接受最后一个功能,但是当我尝试将两个文件一起绘制时,似乎 .shp 文件填充了整个区域。我无法可视化坐标数据。
谢谢你,我真的很感谢你的帮助!
ps.:如果您知道任何有关这些问题的材料,请我很乐意看看
r - spatstat 上的 Ripley 的 K 函数出现错误
我想要一些有关 spatstat 上的 Kest 功能的帮助。我正在尝试为我的模拟数据计算 de Ripley's K。当我尝试运行该函数时,R 向我显示了一些错误,即:
错误:7 段不完全位于窗口内。
但是当我绘制地图时,所有点都完美地显示在观察窗口内。我正在使用以下代码:
我试图查看小插图,但我没有成功解决这个问题。
以下链接包含数据文件。如果您需要其他内容,请告诉我:https ://drive.google.com/folderview?id=0B_cWshGkkAzjbHBnemNuMzh1Q1k&usp=sharing
谢谢你,我真的很感谢你的帮助!