Ripley's K Function (Second Order Intensity Function) Python

Ripley's K Function (Second Order Intensity Function) Python

I am looking for Ripley's k function implementation in Python. But so far haven't been able to find any spatial modules implementing this in scipy or elsewhere.

I have created Voronoi tessellation of a fibre composite and need to perform analysis using Ripley's K and pair distribution functions compared to a Poisson distribution.

Cannot upload images-not enough rep.

2

2 Answers

Solved my problem, this is for others looking to do same analysis.

Definitely recommend using R for spatial analysis. A transfer from python is simple because all you need is coordinates of your point pattern.

  1. Write a csv of x,y and z coordinates of your points using python

  2. R has good functionality of reading csv using command read.csv("filename"). Make sure the directory is set properly using setwd command.

  3. Convert the csv you just read to a ppp(point pattern which R understands) using as.ppp command.
  4. Continue to use Kest, Gest etc for required spatial analysis.

Cheers.

Now (2021) there is a Ripley's K implementation for Python, easy to install

And also a more comprehensive package that includes Ripley's K, G, F and includes MC simulations for p-values, very nice.

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Sarah Jenkins
Author

Sarah Jenkins

Sarah Jenkins is a veteran tech journalist with over 12 years of experience covering artificial intelligence, mobile innovations, and digital ethics. Her insights have appeared in leading technology publications worldwide.