D:\SDP_S3_Uploads\draft_products\UER_canopy_ht_33cm_v2.tif en dataset This is a map of vegetation canopy height suitable for tree detection and crown segmentation, derived from high-density airborne LiDAR scans collected in August - September 2015 and 2019. Classified LiDAR point clouds from both datasets were re-projected to a common coordinate system (EPSG:32613, WGS84 UTM Zone 13N). A 0.33m canopy height model was constructed for each dataset by triangulating ground returns and subtracting the elevation of the ground surface. A pitfree algorithm (Khosravipour et al. 2014) was then used to reconstruct canopy geometry. All point cloud processing was done using the R package lidR version 2.1. The two canopy height models from the 2015 and 2019 LiDAR datasets were mosaiced to produce a continous dataset, with the 2019 data taking priority where there was coverage (approximately 80% of the domain). Buildings and topographic artifacts were removed by setting canopy heights to zero within 1m of building footprints in the Microsoft Bing dataset, and in areas with a negative NDVI in NAIP 4-band aerial imagery. References: Khosravipour, A., Skidmore, A. K., Isenburg, M., Wang, T., & Hussin, Y. A. (2014). Generating pit-free canopy height models from airborne lidar. Photogrammetric Engineering & Remote Sensing, 80(9), 863-872. Elevation Environment Geoscientific Information Ian Breckheimer Rocky Mountain Biological Laboratory Research Scientist ikb@rmbl.org owner Creative Commons Attribution 4.0 PROJCRS["WGS 84 / UTM zone 13N",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4326]],CONVERSION["UTM zone 13N",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",-105,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",0.9996,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["unknown"],AREA["World - N hemisphere - 108°W to 102°W - by country"],BBOX[0,-108,84,-102]],ID["EPSG",32613]] +proj=utm +zone=13 +datum=WGS84 +units=m +no_defs 3097 32613 EPSG:32613 WGS 84 / UTM zone 13N utm EPSG:7030 false