ForestGapR (version 0.0.2)

GapsSpatPattern: Spatial Pattern Analysis of Forest Canopy Gaps

Description

This function computes second order statistics of forest canopy gaps (raster) to SpatialPointsDataFrame-class objects

Usage

GapsSpatPattern(gap_SPDF_layer,chm_layer)

Arguments

gap_SPDF_layer

A SpatialPointsDataFrame-class object of the forest canopy gaps. Output of (GapSPDF) function. An object of the classs SpatialPointsDataFrame-class

chm_layer

ALS-derived Canopy Height Model (CHM) RasterLayer (raster) object. An object of the classs RasterLayer.

Value

A plot with Ripley's K- and L-functions. Value of Clark-Evans index (R) and test for randomness (R=1), aggregation (R<1) or uniform distribution (R>1).

References

spatstat package,see Kest,Lest, and clarkevans.test.

Examples

# NOT RUN {
#Loading raster and viridis libraries
library(raster)
library(viridis)

# ALS-derived CHM from Fazenda Cauxi - Brazilian tropical forest
data(ALS_CHM_CAU_2012)
data(ALS_CHM_CAU_2014)

# set height thresholds (e.g. 10 meters)
threshold <- 10
size <- c(1,1000) # m2

# Detecting forest gaps
gaps_cau2012 <- getForestGaps(chm_layer = ALS_CHM_CAU_2012, threshold=threshold, size=size)
gaps_cau2014 <- getForestGaps(chm_layer = ALS_CHM_CAU_2014, threshold=threshold, size=size)

# Converting raster layers to SpatialPolygonsDataFrame
gaps_cau2012_spdf <- GapSPDF(gap_layer = gaps_cau2012)
gaps_cau2014_spdf <- GapSPDF(gap_layer = gaps_cau2014)

# Spatial pattern analysis of each year
gaps_cau2012_SpatPattern <- GapsSpatPattern(gaps_cau2012_spdf, ALS_CHM_CAU_2012)
gaps_cau2014_SpatPattern <- GapsSpatPattern(gaps_cau2014_spdf, ALS_CHM_CAU_2014)
# }