Calc Average Plot Height with DSM

[HaleySchuhl]

Thinking ahead to types of plantcv-geospatial workflows that we would like to support, Digital Surface Models (DSM) are a powerful tool for calculating height estimates. These estimations tend to be the most accurate to ground truth measurements when the "ground height" region is calculated per plot/Region of Interest to account for minor changes in elevation across the entire area.

Function to read in DSM
Can rasterio already support this and we just wrap the function?

Function to produce geoJSON of labeled ground regions per plot
Would it be useful to have a function that takes a plot gridcell geoJSON/shapefile and creates a complimentary geoJSON of regions, from which an analyze_dsm function would calculate the average ground elevation per plot. Option to calculate one-sided approach (e.g. define multi-polygon geoJSON that represents $5cms of soil directly to the $north of each plot) or on all sides (e.g. in situations with circular ROIs representing replicates in the experimental layout, dilate beyond the plot polygon by $5cm and consider this "picture frame" of area to calculate the average soil height per plot).

Function to calculate height per plot

• for each plot in the shapefile --> plot height = (plot elevation - ground elevation)
• scale to real world measurements when resolution information available?
• store Outputs to plantcv class object