GapsSpatPattern: Spatial Pattern Analysis of Forest Canopy Gaps

Description Usage Arguments Value Author(s) References Examples

View source: R/GapsSpatPattern.R

Description

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

Usage

1
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).

Author(s)

Ruben Valbuena and Carlos Alberto Silva.

References

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

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
## 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)

## End(Not run)

ForestGapR documentation built on May 2, 2019, 1:26 a.m.