Ripley's K 函数(二阶强度函数)Python
Ripley's K Function (Second order intensity function) Python
我正在 Python 中寻找 Ripley 的 k 函数实现。但到目前为止,还没有在 scipy 或其他地方找到任何实现此功能的空间模块。
我已经创建了纤维复合材料的 Voronoi 镶嵌,需要使用 Ripley's K 和对分布函数与泊松分布进行分析。
无法上传图片-代表不足
解决了我的问题,这适用于希望进行相同分析的其他人。
绝对推荐使用 R 进行空间分析。从 python 转移很简单,因为您只需要点图案的坐标。
使用 python
编写点的 x、y 和 z 坐标的 csv
R 具有使用命令读取 csv 的良好功能
read.csv("filename")。使用 setwd 命令确保目录设置正确。
- 将刚刚读取的csv转换为ppp(point
R 理解的模式)使用 as.ppp 命令。
- 继续使用
所需空间分析的 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
我正在 Python 中寻找 Ripley 的 k 函数实现。但到目前为止,还没有在 scipy 或其他地方找到任何实现此功能的空间模块。
我已经创建了纤维复合材料的 Voronoi 镶嵌,需要使用 Ripley's K 和对分布函数与泊松分布进行分析。
无法上传图片-代表不足
解决了我的问题,这适用于希望进行相同分析的其他人。
绝对推荐使用 R 进行空间分析。从 python 转移很简单,因为您只需要点图案的坐标。
使用 python
编写点的 x、y 和 z 坐标的 csv
R 具有使用命令读取 csv 的良好功能 read.csv("filename")。使用 setwd 命令确保目录设置正确。
- 将刚刚读取的csv转换为ppp(point R 理解的模式)使用 as.ppp 命令。
- 继续使用 所需空间分析的 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