csrplus (version 1.03-0)

mwin: Moving Window Approach to Quadrat Counts

Description

Calculates quadrat counts to estimate the intensity of a spatial point process through the moving window approach proposed by Bailey and Gatrell (1995). Event counts are computed within a window of a set size over a fine lattice of points within the region of observation.

Usage

mwin(xcoord, ycoord, boundaryx, boundaryy, gridx, gridy, windowsizel, windowsizew)

Arguments

xcoord

a vector containing the x-coordinates of the observed point process

ycoord

a vector containing the y-coordinates of the observed point process

boundaryx

a vector containing the x-coordinates of the boundary

boundaryy

a vector containing the y-coordinates of the boundary

gridx

integer for the number of column lattice points

gridy

integer for the number of row lattice points

windowsizel

integer for the length of the window

windowsizew

integer for the width of the window

Value

xgrid a vector of the x-coordinates of the lattice points within the boundaryygrid a vector of the y-coordinates of the lattice points within the boundaryquadrat a vector containing the number of events within each window sampled along the lattice

Source

Cressie, Noel. "Chapter 8: Spatial Point Patterns." Statistics for Spatial Data. Revised ed. New York: John Wiley and Sons, 1993. N. pag. Print. Gatrell, Anthony C. "Chapter 3: Introductory Methods for Point Patterns." Interactive Spatial Data Analysis. By Trevor C. Bailey. N.p.: Routledge, 1995. N. pag. Print.

Details

The function first constructs a rectangular space based on the maximum and minimum values of the boundary. It then places a lattice over the rectangle, with a number of points in the lattice defined by the user. The "point.in.polygon" function determines whether or not a specific lattice point falls within the region of observation. If so, the function counts the number of events in a window of a specified size centered on that point. If not, the function moves on to the next lattice point. The function returns the location of lattice points within the region of observation in addition to the quadrat count at each point.

Examples

# To load data corresponding to the location of earthquakes in California:
data(quake)

# To load data corresponding to the boundary:
data(boundary)

# To compute quadrat counts with a 40 x 40 lattice and 1 x 1 unit window:
m <- mwin(quake[,3], quake[,2], boundary[,1], boundary[,2], 40, 40, 1, 1)

# To plot the results (with the shading corresponding to the quadrat count):
layout(matrix(c(1,2), nc=2), widths = c(4, 1))
palette(rev(heat.colors((max(as.numeric(m$quadrat))-min(as.numeric(m$quadrat))))))
plot(m$xgrid, m$ygrid, col=m$quadrat, pch=15, cex=.8,
     xlab="X-Coordinates", ylab="Y Coordinates", main="Quadrat Count")
lines(boundary[,1], boundary[,2])
breaks <- seq(min(as.numeric(m$quadrat)), (max(as.numeric(m$quadrat))), by=1)
plot.new()
plot.window(xlim = c(0, 1),ylim = range(breaks),xaxs = "i", yaxs = "i")
rect(0, breaks[-length(breaks)],1, breaks[-1],
col = rev(heat.colors(length(breaks) - 1)))
axis(2)