Decision tree segmentation method to grow individual tree crowns based on 'itcSegment'

Source: R/fa_crownseg.R

Segmentation of individual tree crowns as polygons based on a LiDAR derived canopy height model. Michele Dalponte: R package itcSegment. M. Dalponte, F. Reyes, K. Kandare, and D. Gianelle, "Delineation of Individual Tree Crowns from ALS and Hyperspectral data: a comparison among four methods," European Journal of Remote Sensing, Vol. 48, pp. 365-382, 2015.

chmseg_ITC(
  chm = NULL,
  EPSG = 3064,
  movingWin = 7,
  TRESHSeed = 0.45,
  TRESHCrown = 0.55,
  minTreeAlt = 2,
  maxCrownArea = 100
)

Arguments

chm

raster*, Canopy height model in raster or SpatialGridDataFrame file format. Should be the same that was used to create the input for treepos.

EPSG

character. The EPSG code of the reference system of the CHM raster image.

movingWin

numeric. Size (in pixels) of the moving window to detect local maxima. itcSegment

TRESHSeed

numeric. seeding threshold. itcSegment

TRESHCrown

numeric. crowns threshold. itcSegment

minTreeAlt

numeric. Height threshold (m) below a pixel cannot be a local maximum. Local maxima values are used to define tree tops.itcSegment

maxCrownArea

numeric. A single value of the maximum projected tree crown area allowed. Default 100 sqm. height of treepos.

Value

sp* object

Examples

 require(uavRst)
 data(chm_seg)
##- segmentation
crownsITC<- chmseg_ITC(chm = chm_seg[[1]],
                       EPSG =3064,
                       movingWin = 7,
                       TRESHSeed = 0.45,
                       TRESHCrown = 0.55,
                       minTreeAlt = 5,
                       maxCrownArea = 50)
#> Warning: CRS object has no comment
#> Warning: CRS object has no comment
#> Warning: Discarded datum Istituto_Geografico_Militaire_1995 in Proj4 definition
#> Warning: CRS object has no comment
#> Warning: A new CRS was assigned to an object with an existing CRS:
#> +proj=utm +zone=32 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
#> without reprojecting.
#> For reprojection, use function spTransform

##- visualisation
raster::plot(crownsITC)