Lidar derived Trees Per Acre, Minnesota
Go to Section:- Overview
- Data Quality
- Data Organization
- Coordinate System
- Attributes
- Distribution - Get Data
- Metadata Reference
Section 1: Overview
Originator:Minnesota Department of Natural Resources (DNR)
Title:
Lidar derived Trees Per Acre, Minnesota
Abstract:
The Lidar derived Trees Per Acre (TPA) raster layer (at 20 m spatial resolution) is a modelled forest inventory attribute based on integration of Minnesota's 3DEP USGS Lidar and temporally coincident plot based inventory (PBI) data collected by the MN DNR Division of Forestry Resource Assessment office.
For more information about all metrics and models, see the Resource Assessment Glossary and Frequently Asked Questions documents, found on the Resource Assessment Program website.
Frequently Asked Questions
Resource Assessment Glossary
Resource Assessment | MN DNR
Purpose:
TPA is a spatially explicit raster output for the entire Lidar Acquisition Block (LAB). Each pixel in the raster represents TPA (in sq ft per acre) for all live trees of five inches or larger Diameter at Breast Height (DBH). The spatial inventory modeling of TPA required field inventory data and co-located lidar data from many 0.1 acre PBI plots throughout the LAB. The sample tree DBH data from the PBI plots are used to estimate basal area of individual tree and plot level basal area (response variable) is calculated by summing up the basal areas of individual trees per plot. The lidar point cloud extracted for the same plots provide more than 40 lidar metrics as the spatial predictors. A regression model fitted to the training data is extended spatially over the entire area of interest using the most significant selected lidar predictors at 20 m spatial resolution.
Time Period of Content Date:
Currentness Reference:
The 3DEP USGS Lidar data used to produce this modeled forest inventory attribute, was acquired by LAB at various time periods in the state, starting in 2019, and up to 2025.
Progress:
In work
Maintenance and Update Frequency:
None Planned
Spatial Extent of Data:
Statewide
Bounding Coordinates:
-97.23
-89.53
49.37
43.5
Place Keywords:
Minnesota
Theme Keywords:
lidar, remote sensing, forestry, 3DEP USGS
Theme Keyword Thesaurus:
Access Constraints:
None. Please see 'Distribution Info' for details.
Use Constraints:
None. Users are advised to read the dataset's metadata thoroughly and to have background knowledge in forestry to understand appropriate use and data limitations.
Browse Graphic:
None available
Associated Data Sets:
Minnesota's 3DEP USGS Lidar, LAB extent (polygon feature class)
Section 2: Data Quality
Attribute Accuracy:
Logical Consistency:
Completeness:
The data included is for the Lidar Acquision Blocks where the data is available.
Horizontal Positional Accuracy:
N/A
Vertical Positional Accuracy:
N/A
Lineage:
The Lidar inputs are from certified 3DEP USGS data as referenced by
RLB (2021) Project Report
https://rockyweb.usgs.gov/vdelivery/Datasets/Staged/Elevation/metadata/MN_RainyLake_2020_B20/USGS_MN_RainyLake_2020_B20_Project_Report.pdf
LSB (2021) Project Report
https://rockyweb.usgs.gov/vdelivery/Datasets/Staged/Elevation/metadata/MN_LakeSuperior_2021_B21/USGS_MN_LakeSuperior_2021_B21_Project_Report.pdf
PiCo (2019) Project Report
https://rockyweb.usgs.gov/vdelivery/Datasets/Staged/Elevation/metadata/WI_Oshkosh_3Rivers_2018_D19/WI_Oshkosh_3Rivers_2018_D19_WP_Report.pdf
The mask created for Lidar Derived forest inventory models is based on two parameters. First, the areas in the forest inventory models that have negative values. The negative values in a forest inventory model correspond to areas where lidar returns are absent, or are in a condition that was not sampled during field collection. Second, the mask omits areas where less than 10% of the first returns are above 3m from the ground.
TPA is modeled using multivariate relationships between field measured plot (0.1 acre) data and lidar-derived grid metrics (at 20 m resolution). In the initial step, a large number of lidar metrics (describing canopy density, elevation distributions and various strata density) are developed for each PBI location which together with the plot-level response variables are subjected to random forest model selection algorithms to determine the important metrics that contribute the most to prediction of the forest inventory variables. The predictor variables selected in this way are further subjected to stepwise backward and forward regression to have a parsimonious set of significant predictors that are finally used to fit a weighted multiple linear regression model (also called generalized least square model, ‘gls’) in the ‘nlme’ package in R statistical programming. In general, the model weight is based on a predictor that yields highest correlation with the response inventory variable. The weighted multiple linear regression model is then extended spatially (wall-to-wall) using the same parsimonious set of predictors as grid metrics at 20 m resolution.
The version of the TPA Model refers to the sequential order of running the Model on updated PBI data.
The multiple linear regression model for AGB in the following LABs is in the form:
where a0 is intercept and a1-a6 are linear model coefficients
RLB (2021) - Model v2
TPA = -a0 + (a1 * covMeanIst.tif) + (a2 * ElevCV.tif) + (a3 * ElevIQ.tif) + (a4 * ElevMADmed.tif) + (a5 * strata3.tif) + (a6 * ElevP05.tif) + (a7 * strata4.tif)
LSB (2021) - Model v2
TPA = -a0 + (a1 * cov3mIst.tif) + (a2 * strata1.tif) + (a3 * strata2.tif) + (a4 * strata3.tif) + (a5 * strata4.tif) - (a6 * ElevL3.tif) + (a7 * ElevCV.tif)
PiCo (2019) - Model v3
TPA = a0 - (a1 * ElevSkew.tif) + (a2 * cov3mIst.tif) + (a3 * strata4.tif)
Section 3: Spatial Data Organization (not used in this metadata)
Section 4: Coordinate System
Horizontal Coordinate Scheme:
Universal Transverse Mercator
UTM Zone Number:
15
Horizontal Datum:
NAD83(2011)
Horizontal Units:
meters
Vertical Datum:
Vertical Units:
Depth Datum:
Depth Units:
Cell Width:
20m
Cell Height:
20m
Section 5: Attributes
Overview:
Value range 0+ represents the TPA in number of trees per acre.
Detailed Citation:
Support for this project is provided by the State of Minnesota and the Minnesota Environment and Natural Resources Trust Fund (ENTRF) as recommended by the Legislative-Citizen Commission on Minnesota Resources (LCCMR).
Resource Assessment | MN DNR
Table Detail:
Section 6: Distribution
Publisher:
Minnesota Department of Natural Resources (DNR)
Publication Date:
01/10/2025
Contact Person Information:
Jason Langenecker,
Assistant Resource Assessment Supervisor for Field Services
Minnesota DNR - Forestry, Resource Assessment
483 Peterson rd
Grand Rapids,
MN
55744
Phone: 218-322-2519
Email: jason.langenecker@state.mn.us
Distributor's Data Set Identifier:
Ordering Instructions:
Lidar data may be accessed directly from the 3DEP USGS program via the National Map https://apps.nationalmap.gov/viewer/.
For access to the Lidar derived forest inventory raster datasets, send a message to the Distribution Contact. Online Linkage:
None available
Section 7: Metadata Reference
Metadata Date:
01/10/2025
Metadata Standard Name:
Minnesota Geographic Metadata Guidelines
Metadata Standard Version:
1.2
Metadata Standard Online Linkage: