GapStats {ForestGapR}

Forest Canopy Gaps Stats

Description

This function computes a series of forest canopy gap statistics

Usage

GapStats(gap_layer, chm_layer)

Arguments

gap_layer

ALS-derived gap as raster::RasterLayer object (output of getForestGaps() function). An object of the class RasterLayer.

chm_layer

ALS-derived Canopy Height Model (CHM) raster::RasterLayer used in getForestGaps() function. An object of the class raster::RasterLayer.

Value

A data.frame containing forest canopy gap statistics

List of forest gaps statistics:

• gap_id: gap id

• gap_area - area of gap (m2)

• chm_max - Maximum canopy height (m) within gap boundary

• chm_min - Minimum canopy height (m) within gap boundary

• chm_mean - Mean canopy height (m) within gap boundary

• chm_sd - Standard Deviation of canopy heights (m) within gap boundary

• chm_gini - Gini Coefficient of canopy heights (m) within gap boundary

• chm_range - Range of canopy heights (m) within gap boundary

Author(s)

Carlos Alberto Silva.

Examples

# Loading raster library
library(raster)

# ALS-derived CHM over Adolpho Ducke Forest Reserve - Brazilian tropical forest
data(ALS_CHM_CAU_2012)

# set height thresholds (e.g. 10 meters)
threshold <- 10
size <- c(5, 10^4) # m2

# Detecting forest gaps
gaps_duc <- getForestGaps(chm_layer = ALS_CHM_DUC, threshold = threshold, size = size)

# Computing basic statistis of forest gap
gaps_stats <- GapStats(gap_layer = gaps_duc, chm_layer = ALS_CHM_DUC)