Ripley's K 函数(二阶强度函数)Python

Ripley's K Function (Second order intensity function) Python

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

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

无法上传图片-代表不足

解决了我的问题,这适用于希望进行相同分析的其他人。

绝对推荐使用 R 进行空间分析。从 python 转移很简单,因为您只需要点图案的坐标。

  1. 使用 python

  2. 编写点的 x、y 和 z 坐标的 csv
  3. R 具有使用命令读取 csv 的良好功能 read.csv("filename")。使用 setwd 命令确保目录设置正确。

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

干杯。

现在 (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