瑞普利K函数（二阶强度函数）Python

Question

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

我已经创建了纤维复合材料的 Voronoi 细胞分割，并需要使用 Ripley's K 和对比泊松分布的双分布函数进行分析。

无法上传图像-声望不够。

- Mehtab Pathan

---

您能否考虑通过RPy简单调用R函数？如果是这种情况，您可以使用R的spatstat包。 - Ege Rubak

谢谢！之前不知道RPy，但听起来可行。 - Mehtab Pathan

Answer 1

我的问题已经解决，以下内容为希望其他人进行相同分析的建议。

强烈推荐使用R进行空间分析。如果您想要从Python转移过来，非常简单，因为您只需要知道点模式的坐标。

1. 使用Python编写包含x、y和z点坐标的csv文件。

2. R有很好的读取csv功能，使用命令read.csv("filename")即可。请确保使用setwd命令正确设置目录。

3. 使用as.ppp命令将刚刚读取的csv转换为ppp（R所理解的点模式）。
4. 继续使用Kest、Gest等进行所需的空间分析。

祝好！

- Mehtab Pathan

Answer 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

- Ricardo Guerreiro