rasterize: Calculate count, basal area or agb per quadrat.

In forestgeo/ctfs: Tools for the Analysis of Forest Dynamics

Description

Function to calculate the count, basal area or agb per quadrat by selecting quadrats of variable sizes.

Usage

rasterize(x, y, z = NULL, gridsize = 20, plotdim = c(1000, 500),
  FUN = sum, graph = FALSE)

Arguments

x, y	Coordinates for the point pattern.
z	Marks for the marked point process, can be basal area, agb, etc. Needed only for type = "marked".
gridsize	Side of the square quadrate.
plotdim	The x and y dimensions of the plot.
FUN	Function to apply to the point pattern when z is provided. By defult it sums the values of z per quadrat.
graph	Logical. If TRUE plots the heat map of raster data.
type	Type of rasterization. "point" runs a simple point pattern counting the number of individuals per quadrat; "marked" does a marked point pattern, where a function (FUN) is applied to the variable z in each quadrat, for example a sum of basal areas.

Value

A list containing a matrix with the raster data.

Author(s)

Matteo Detto and Tania Brenes

Examples

## Not run: 
onesp = subset(bciex::bci12t1mini, sp == "rinosy")  

# plots the density of the sp in the plot  
rast1 = rasterize(
  x = onesp$gx,
  y = onesp$gy,
  gridsize = 5,
  plotdim = c(1000, 500),
  type = 'point',
  graph = TRUE
)

## End(Not run)

forestgeo/ctfs documentation built on May 3, 2019, 6:44 p.m.