Description Usage Arguments Examples
Produce a raster map of LAI. The resolution of the raster depends of grain.size choosed on lad.voxel() funtion.
1 | lai.raster(VOXELS_LAD, min = 1, relative.value = NULL)
|
VOXELS_LAD |
3D grid of LAD values (output of lad.voxels() function) |
min |
mix canopy height |
relative.value |
LAI map can be made in percentage of a relative lai value (indicate for effective LAI) A Leaf Area Index (LAI) |
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 27 28 29 30 31 | library(raster)
# Get the example laz file
normlas.file = system.file("extdata", "lidar_example.laz", package="leafR")
# Calculate LAD from voxelization
# use thicker grain size to avoid voxels
# without returns
VOXELS_LAD.5 = lad.voxels(normlas.file,
grain.size = 5, k=1)
#Map using absolute values
lai_raster = lai.raster(VOXELS_LAD.5)
plot(lai_raster)
#############################
## RELATIVE LAI Raster
######################
# Calculate voxels LAD with finer grain size for
# better estimation of LAI
VOXELS_LAD = lad.voxels(normlas.file,
grain.size = 2)
# Calculate the LAD profile
lad_profile = lad.profile(VOXELS_LAD)
#Calculate LAI derived from LAD profile
lidar.lai = lai(lad_profile)
#Map using relative values (%)
relative.lai_raster = lai.raster(VOXELS_LAD.5, relative.value = lidar.lai)
plot(relative.lai_raster)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.