【问题标题】:Ripley's K Function (Second order intensity function) PythonRipley's K 函数(二阶强度函数)Python
【发布时间】:2015-07-01 23:05:10
【问题描述】:

我正在寻找 Python 中的 Ripley 的 k 函数实现。但到目前为止,还没有找到任何在 scipy 或其他地方实现此功能的空间模块。

我创建了纤维复合材料的 Voronoi 镶嵌,需要使用 Ripley 的 K 和对分布函数与泊松分布进行分析。

无法上传图片 - 没有足够的代表。

【问题讨论】:

  • 您可以选择通过 RPy 简单地调用 R 函数吗?在这种情况下,您可以为 R 使用 spatstat 包。
  • 谢谢!以前不知道 RPy,但这听起来可行。

标签: python spatial


【解决方案1】:

解决了我的问题,这是为其他希望进行相同分析的人准备的。

绝对推荐使用 R 进行空间分析。从 python 传输很简单,因为您只需要点模式的坐标。

  1. 使用 python 编写点的 x、y 和 z 坐标的 csv

  2. R 具有良好的使用命令读取 csv 的功能 读取.csv(“文件名”)。确保使用 setwd 命令正确设置目录。

  3. 将刚才读取的csv转换成ppp(point R 理解的模式)使用 as.ppp 命令。
  4. 继续使用 用于所需空间分析的 Kest、Gest 等。

干杯。

【讨论】:

    【解决方案2】:

    现在(2021 年)有一个用于 Python 的 Ripley's K 实现,易于安装

    https://pypi.org/project/ripleyk/

    还有一个更全面的包,包括 Ripley 的 K、G、F 和 p 值的 MC 模拟,非常好。

    https://github.com/pysal/pointpats/blob/master/notebooks/distance_statistics-numpy-oriented.ipynb

    【讨论】:

      猜你喜欢
      • 1970-01-01
      • 2012-03-19
      • 1970-01-01
      • 2014-08-19
      • 2011-07-05
      • 2015-04-16
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      相关资源
      最近更新 更多