Skip to content

mpl in spatstat

1 message · rohan sadler

#
Hi Suzanne,

It sounds like you are dealing with an Inhomogeneous Poisson Point 
Process (the simplest). The intensity function lambda(mu) of this 
process is the density of points near mu. i.e. the process is poisson 
locally, and the intensity varies over the region according to some 
underlying function or surface (such as soil chemistry)

e.g.
 > plot(rpoispp(function(x,y){300*cos(2*pi*x)},300))  # taken from 
Adrian Baddeley's notes

The first step would then be interpolate your soil characteristics into 
global functions (either through delaunay triangulation, smoothing 
method (loess), or krigging). This allows you to estimate soil 
characteristics where your trees are, and derive vectors of the same 
length (although with some error).

There may also be away of involving the idea of a global funciton into 
the actual inference process, but I haven't thought about that deeply 
and my memory cells aren't stirred to recall a suitable reference.

Hope this was of some use.

Rohan Sadler
Suzanne E. Blatt wrote: