crownmetrics is a R package for computing tree crown
metrics from basic field measurements. It provides two categories of
functions:
Throughout this document the following symbols are used:
| Symbol | Variable | Unit |
|---|---|---|
| \(CW\) | Crown width (mean diameter) | meters |
| \(CL\) | Crown length | meters |
| \(H\) | Total tree height | meters |
| \(DBH\) | Diameter at breast height | centimeters |
Crown volume is approximated by fitting a geometric solid to the crown envelope. The appropriate shape depends on the species and management history (Zhu, Kleinn & Nölke, 2021).
The ellipsoid is the most widely used approximation for broadleaf and urban trees. When \(CW = CL\) the formula reduces to a sphere.
\[V_{\text{ellipsoid}} = \frac{4}{3}\,\pi \left(\frac{CL}{2}\right)\left(\frac{CW}{2}\right)^{2}\]
Suitable for conifers with narrow, pointed crowns (e.g. Pinus spp., Picea spp.).
\[V_{\text{cone}} = \frac{1}{3}\,\pi \left(\frac{CW}{2}\right)^{2} CL\]
Represents trees with a uniform crown width from base to top.
\[V_{\text{cylinder}} = \pi \left(\frac{CW}{2}\right)^{2} CL\]
An intermediate shape between the cone and the cylinder, often considered a reasonable approximation for broadleaf trees.
\[V_{\text{paraboloid}} = \frac{1}{2}\,\pi \left(\frac{CW}{2}\right)^{2} CL\]
Used for trees with wide, flat crowns (e.g. Pinus pinea). The vertical extent of the crown solid is approximated by the stem diameter at breast height converted to meters (\(DBH / 100\)).
\[V_{\text{fan}} = \frac{\pi\,CL^{2}}{4} \cdot \frac{DBH}{100}\]
The relationship between shapes is fixed for given \(CW\) and \(CL\): \(V_{\text{cylinder}} = 2\,V_{\text{paraboloid}} = 3\,V_{\text{cone}}\), with the ellipsoid falling between the paraboloid and the cylinder.
cw <- 3.5 # crown width (m)
cl <- 5.0 # crown length (m)
volumes <- c(
ellipsoid = crown_volume_ellipsoid(cw, cl),
cone = crown_volume_cone(cw, cl),
cylinder = crown_volume_cylinder(cw, cl),
paraboloid = crown_volume_paraboloid(cw, cl)
)
round(volumes, 2)
#> ellipsoid cone cylinder paraboloid
#> 32.07 16.04 48.11 24.05crown_volume() selects the model via its
shape argument, accepting "ellipsoid",
"cone", "cylinder", "paraboloid",
or "fan":
crown_volume(crown_width = 3.5, crown_length = 5.0, shape = "ellipsoid")
#> [1] 32.07043
crown_volume(crown_width = 3.5, crown_length = 5.0, shape = "cone")
#> [1] 16.03521
crown_volume(crown_length = 5.0, dbh = 22.0, shape = "fan")
#> [1] 4.31969Two types of area are available.
These represent the two-dimensional silhouette of the crown as seen from the side (McPherson & Rowntree, 1988; Zhu, Kleinn & Nölke, 2021).
Ellipse — used with the ellipsoid and paraboloid volume models:
\[A_{\text{ellipse}} = \pi \left(\frac{CL}{2}\right)\left(\frac{CW}{2}\right)\]
Triangle — used with the cone volume model:
\[A_{\text{triangle}} = \frac{CW \cdot CL}{2}\]
Rectangle — used with the cylinder volume model:
\[A_{\text{rectangle}} = CW \cdot CL\]
Fan / umbrella — used with the fan volume model:
\[A_{\text{fan}} = \frac{\pi\,CL}{4}\]
crown_area(crown_width = 3.5, crown_length = 5.0, shape = "ellipse")
#> [1] 13.74447
crown_area(crown_width = 3.5, crown_length = 5.0, shape = "triangle")
#> [1] 8.75
crown_area(crown_width = 3.5, crown_length = 5.0, shape = "rectangle")
#> [1] 17.5
crown_area(crown_length = 5.0, shape = "fan")
#> [1] 3.926991Morphometric indices are dimensionless ratios that characterise tree architecture, stability, and competitive status. All indices follow Durlo & Denardi (1998) and Burger (1939).
The fraction of total height occupied by the living crown. Higher values indicate greater photosynthetic potential and crown vitality.
\[CR = \frac{CL}{H}\]
The ratio of crown width to crown length. Values \(> 1\) indicate wide, flat crowns; values \(< 1\) indicate tall, narrow crowns.
\[CF = \frac{CW}{CL}\]
The ratio of total height to DBH (both expressed in meters). Higher values indicate more slender trees with greater susceptibility to wind and snow damage.
\[SLI = \frac{H}{DBH / 100}\]
Expresses how many times wider the crown is than the trunk diameter (both in meters). Reflects the tree’s capacity to occupy horizontal space relative to its stem size.
\[SAI = \frac{CW}{DBH / 100}\]
The ratio of crown width to total tree height. Describes the lateral competitive reach of the crown relative to tree height.
\[SCI = \frac{CW}{H}\]
The ratio of crown projection area to stem basal area. Expresses how many times more ground area the crown occupies compared to the stem cross-section. Note that the \(\pi\) terms cancel, so \(VSI = SAI^{2}\).
\[VSI = \frac{CPA}{BA} = \frac{\pi\,(CW/2)^{2}}{\pi\,(DBH/200)^{2}} = \left(\frac{CW}{DBH/100}\right)^{2}\]
| Symbol | Name | Equation | Interpretation |
|---|---|---|---|
| CR | Crown Ratio | \(CL\,/\,H\) | Crown occupancy; higher = more leaf area |
| CF | Crown Form | \(CW\,/\,CL\) | \(>1\): wide flat; \(<1\): narrow tall |
| SLI | Slenderness | \(H\,/\,(DBH/100)\) | Higher = more wind-sensitive |
| SAI | Salience Index | \(CW\,/\,(DBH/100)\) | Crown breadth relative to stem |
| SCI | Scope Index | \(CW\,/\,H\) | Lateral reach relative to height |
| VSI | Vital Space Index | \((CW/DBH_{m})^{2}\) | Growing space efficiency |
crown_metrics()For inventory datasets, crown_metrics() computes every
volume, area, and morphometric index in one pass, returning the original
data frame with all metric columns appended. Crown shape is specified
per-row via integer codes or character strings:
| Integer | Character | Solid |
|---|---|---|
| 0 | "ellipsoid" |
Ellipsoid |
| 1 | "cone" |
Cone |
| 2 | "cylinder" |
Cylinder |
| 3 | "paraboloid" |
Paraboloid |
| 4 | "fan" |
Fan |
inventory <- data.frame(
tree_id = 1:6,
species = c("Araucaria", "Pinus", "Eucalyptus",
"Araucaria", "Pinus", "Eucalyptus"),
crown_width = c(3.0, 2.0, 4.0, 3.5, 2.5, 4.5),
crown_length = c(5.0, 6.0, 3.5, 4.5, 5.5, 3.0),
total_height = c(18.0, 20.0, 12.0, 16.0, 19.0, 11.0),
dbh = c(22.0, 18.0, 25.0, 20.0, 17.0, 28.0),
shape_code = c(0L, 1L, 0L, 0L, 1L, 3L)
)
result <- crown_metrics(inventory, col_shape = "shape_code")
result[, c("tree_id", "species", "crown_volume_m3",
"crown_projection_area_m2", "slenderness", "crown_ratio")]
#> tree_id species crown_volume_m3 crown_projection_area_m2 slenderness
#> 1 1 Araucaria 23.561945 7.068583 81.81818
#> 2 2 Pinus 6.283185 3.141593 111.11111
#> 3 3 Eucalyptus 29.321531 12.566371 48.00000
#> 4 4 Araucaria 28.863383 9.621128 80.00000
#> 5 5 Pinus 8.999354 4.908739 111.76471
#> 6 6 Eucalyptus 23.856469 15.904313 39.28571
#> crown_ratio
#> 1 0.2777778
#> 2 0.3000000
#> 3 0.2916667
#> 4 0.2812500
#> 5 0.2894737
#> 6 0.2727273If no shape column is provided, all rows default to the ellipsoid model.
Burger, H. (1939). Baumkrone und Zuwachs in zwei hiebsreifen Fichtenbeständen. Mitteilungen der Schweizerischen Anstalt für das Forstliche Versuchswesen, 21, 147–176.
Durlo, M. A., & Denardi, L. (1998). Morfometria de Cabralea canjerana, em mata secundária nativa do Rio Grande do Sul. Ciência Florestal, 8(1), 55–66.
McPherson, E. G., & Rowntree, R. A. (1988). Geometric solids for simulation of tree crowns. Landscape and Urban Planning, 15(3–4), 79–83.
Sayn-Wittgenstein, L., & Aldred, A. H. (1972). Tree size from large-scale photos. Photogrammetric Engineering, 38, 971–973.
Sterba, H. (1991). Forstliche Wuchslehre. Universität für Bodenkultur, Wien.
Wink, C., Monteiro, J. S., Reinert, D. J., & Liberalesso, E. (2012). Parâmetros da copa e a sua relação com o diâmetro e altura das árvores de eucalipto em diferentes idades. Scientia Forestalis, 40(93), 57–67.
Zhu, Z., Kleinn, C., & Nölke, N. (2021). Assessing tree crown volume — a review. Forestry, 94(1), 18–35.