2015-04-23 15 views
5

Sto cercando l'implementazione della funzione k di Ripley in Python. Ma finora non sono stati in grado di trovare alcun modulo spaziale che lo implementasse in scipy o altrove.Ripley's K Function (funzione di intensità del secondo ordine) Python

Ho creato una tassellatura Voronoi di un composito di fibre e ho bisogno di eseguire analisi utilizzando le funzioni di distribuzione di Ripley e K rispetto a una distribuzione di Poisson.

Impossibile caricare immagini, non abbastanza rep.

+0

Potrebbe essere un'opzione per voi semplicemente chiamare una funzione R tramite RPy? In questo caso puoi usare il pacchetto 'spatstat' per R. –

+0

Grazie! Non sapevo di RPy prima ma sembra fattibile. –

risposta

2

risolto il mio problema, questo è per gli altri che cercano di fare la stessa analisi.

Consiglio vivamente di utilizzare R per l'analisi spaziale. Un trasferimento da Python è semplice perché tutto ciò di cui hai bisogno sono le coordinate del tuo pattern punto.

  1. Mandaci csv di coordinate x, yez dei punti di utilizzo pitone

  2. R ha una buona funzionalità di lettura csv utilizzando il comando read.csv ("nome file"). Assicurati che la directory sia impostata correttamente usando il comando setwd.

  3. Convertire il csv appena letto in un ppp (punto modello che R comprende) utilizzando il comando as.ppp.
  4. Continua a utilizzare Kest, Gest ecc. Per l'analisi spaziale richiesta.

Cheers.

Problemi correlati