52620_mesquit.tif

Ref: J. M. Rusk Report of the Secretary of Agriculture for the year 1891 (Washington, DC: Government Printing Office, 1892)

image source: http://etc.usf.edu/clipart/52600/52620/52620_mesquit.htm

Meeting with Joel and Adam

Does McClaran et al. 2013 take into account increasing tree age as a possible source of increasing divergence from an idealized allometric state?

Divergence in the height to basal area height to canopy area vs predicted allometry with increasing age
Divergence related to disturbance (wildfires), disease, damage.

Are we focused on biomass, component biomass partitioning?

What are the limitations of the airborne lidar for measuring mesquite biomass and structure (morphology)?
What are the limitations of the airborne lidar for measuring under the mesquites (species diversity + herbaceous biomass)?
- Adam will establish the relationship using his field measurements (observed).
- Compare his model to the lidar estimated canopy area/volume (estimated).
What are the limitations of the SfM?

Do we want to expand to include species diversity (α, β, γ)?

We need to build a Google.Drive document.

I need to share the lidar data files for the trees.

Need to establish figures.

Figure 1: Site map
Figure 2: LiDAR examples of mesquite
Figure 3: Regression of biomass vs height
Figure 4: Regression of canopy vs height
Figure 5: Cumulative density function of biomass across area by tree size (Joel)
Figure 6: Regression of the field derived canopy volume vs the LiDAR derived canopy volume.
Figure 7: Regression of the species diversity to tree size
Figure 8: Regression of the observed diversity (Adam's model) vs the LiDAR estimated diversity

LiDAR Canopy Height Modelling and stem segmentation

First I generate a 1 foot canopy height model in FUSION from the Windows CMD command line:

cd\fusion
## Digital Surface Model
CanopyModel /ascii F:\SRER\DSM\19S15E06_1ft.dtm 1 f f 2 0 2 2 F:\SRER\LAS\19S15E06_LDRY11.las
## Canopy Height Model
CanopyModel /outlier:-1,50 /ground:F:\SRER\DTM\19S15E06_LDRY11.dtm /ascii F:\SRER\CHM\19S15E06_1ft.dtm 1 f f 2 0 2 2 F:\SRER\LAS\19S15E06_LDRY11.las

Next, I open Matlab and select a working directory where my VLM script is saved. I then execute a process to segment out the individual trees at multiple canopy diameter-to-stem height ratio:

% Imports canopy height model (2-ft resolution) and generates stem map
tile19s15e06=asc_import('E:\PAG2011\CHM\ASC\19S15E06_LDRY11.ASC');
% Assigns a 0.75 canopy diameter to height ratio 
tile_v=vlm(tile19s15e06,0.75,2);
vlm_ex=export_utm(tile_v,tile19s15e06);
col_header={'ID','UTME','UTMN','HT','PRED','AREA','EQDIAM','MAJAX','MINAX','MAXHT','MINHT','MEANHT'};
xlswrite('F:\SRER\STEMS\cd0p75_19s15e06.xlsx',vlm_ex,1,'A2');
xlswrite('F:\SRER\STEMS\cd0p75_19s15e06.xlsx',col_header,1,'A1');
% Assigns a 1 canopy diameter to height ratio 
tile_v=vlm(tile19s15e06,1,2);
vlm_ex=export_utm(tile_v,tile19s15e06);
col_header={'ID','UTME','UTMN','HT','PRED','AREA','EQDIAM','MAJAX','MINAX','MAXHT','MINHT','MEANHT'};
xlswrite('F:\SRER\STEMS\cd1_19s15e06.xlsx',vlm_ex,1,'A2');
xlswrite('F:\SRER\STEMS\cd1_19s15e06.xlsx',col_header,1,'A1');
% Assigns a 1.25 canopy diameter to height ratio
tile_v=vlm(tile19s15e06,1.25,2);
vlm_ex=export_utm(tile_v,tile19s15e06);
col_header={'ID','UTME','UTMN','HT','PRED','AREA','EQDIAM','MAJAX','MINAX','MAXHT','MINHT','MEANHT'};
xlswrite('F:\SRER\STEMS\cd1p25_19s15e06.xlsx',vlm_ex,1,'A2');
xlswrite('F:\SRER\STEMS\cd1p25_19s15e06.xlsx',col_header,1,'A1');
% Assigns a 1.5 canopy diameter to height ratio 
tile_v=vlm(tile19s15e06,1.5,2);
vlm_ex=export_utm(tile_v,tile19s15e06);
col_header={'ID','UTME','UTMN','HT','PRED','AREA','EQDIAM','MAJAX','MINAX','MAXHT','MINHT','MEANHT'};
xlswrite('F:\SRER\STEMS\cd1p5_19s15e06.xlsx',vlm_ex,1,'A2');
xlswrite('F:\SRER\STEMS\cd1p5_19s15e06.xlsx',col_header,1,'A1');

Allometric Models

McClaran et al. (2013) reported allometric models for mesquite that include biomass and foliar leaf area (Table 2):

f(x)	x	α(95%CI)	β	r²	CF

f(x)	x	α(95%CI)	β	r²	CF
Canopy Area (m²)	Basal Diameter	-1.46 (0.28)	1.34 (0.10	0.96	1.01
Height (m)	Basal Diameter	-0.47 (0.18)	0.52 (0.06)	0.9	1
Canopy Volume (m³)	Basal Diameter	-1.92 (0.39)	1.86 (0.13)	0.96	1.03
Height (m)	Canopy Area	0.11 (0.10)	0.38 (0.04)	0.91	1
Foliar biomass (kg)	Live aboveground biomass	2.38 (0.18)	0.77 (0.05)	0.97	1.01

Table 3 reported the structural measure to biomass models as:

f(x)	x = Basal diameter (cm)				x = Canopy area (m²)				x = Height (m)				x = Volume (m³)

f(x)	x = Basal diameter (cm)				x = Canopy area (m²)				x = Height (m)				x = Volume (m³)
	α(95%CI)	β(95%CI)	r²	CF	α(95%CI)	β(95%CI)	r²	CF	α(95%CI)	β(95%CI)	r²	CF	α(95%CI)	β(95%CI)	r²	CF
Leaf area (m²)	−2.72 (0.34)	1.64 (0.12)	0.96	1.02	−0.90 (0.31)	1.21 (0.08)	0.97	1.01	−1.05 (0.37)	2.94 (0.33)	0.92	1.1	−0.99 (0.20)	0.87 (0.05)	0.97	1.01
Carbon (kg)	−3.77 (0.45)	2.19 (0.16)	0.96	1.06	−1.32 (0.31)	1.60 (0.12)	0.96	1.06	−1.58 (0.45)	3.94 (0.40)	0.93	1.23	−1.45 (0.28)	1.16 (0.08)	0.97	1.04
Nitrogen (g)	−0.41 (0.37)	1.96 (0.13)	0.97	1.02	1.77 (0.26)	1.44 (0.10)	0.97	1.03	1.56 (0.41)	3.51 (0.36)	0.93	1.15	1.67 (0.24)	1.03 (0.06)	0.97	1.02
Total biomass (kg)	−3.02 (0.45)	2.19 (0.16)	0.96	1.06	−0.59 (0.31)	1.60 (0.12)	0.96	1.06	−0.83 (0.45)	3.93 (0.40)	0.93	1.24	−0.71 (0.28)	1.16 (0.08)	0.97	1.04
Live biomass (kg)	−3.02 (0.45)	2.12 (0.16)	0.96	1.06	−0.67 (0.31)	1.55 (0.12)	0.96	1.06	−0.92 (0.41)	3.82 (0.37)	0.94	1.16	−0.79 (0.27)	1.12 (0.07)	0.97	1.03
Foliar (kg)	−4.88 (0.36)	1.67 (0.13)	0.96	1.02	−3.03 (0.20)	1.23 (0.07)	0.97	1.01	−3.20 (0.36)	3.00 (0.32)	0.92	1.09	−3.12 (0.19)	0.89 (0.05)	0.98	1.01
Fine stem (kg)	−3.15 (0.26)	1.52 (0.09)	0.98	1.01	−1.45 (0.21)	1.11 (0.08)	0.97	1.01	−1.61 (0.33)	2.70 (0.30)	0.92	1.07	−1.54 (0.20)	0.80 (0.05)	0.97	1.01
Small stem (kg)	−2.79 (0.82)	1.45 (0.26)	0.84	1.07	−1.50 (0.54)	1.18 (0.19)	0.86	1.05	−1.38 (0.67)	2.65 (0.56)	0.79	1.12	−1.56 (0.54)	0.84 (0.13)	0.87	1.05
Mid stem (kg)	−4.93 (1.68)	2.27 (0.52)	0.79	1.44	−2.77 (1.02)	1.79 (0.34)	0.84	1.23	−2.54 (0.98)	4.04 (0.77)	0.84	1.23	−2.93 (0.88)	1.30 (0.21)	0.88	1.12
Large stem (kg)	−4.30 (2.07)	2.27 (0.59)	0.79	1.07	−2.01 (1.54)	1.74 (0.50)	0.76	1.1	−1.54 (1.79)	3.36 (1.28)	0.64	1.23	−2.32 (1.52)	1.29 (0.33)	0.8	1.06
Dead biomass (kg)	−6.80 (1.36)	2.81 (0.43)	0.88	1.41	−4.09 (1.07)	2.20 (0.37)	0.86	1.62	−3.54 (1.50)	4.68 (1.23)	0.72	7.26	−4.14 (1.12)	1.56 (0.56)	0.85	1.77

Because the LiDAR does not measure basal diameter we must estimate biomass from other canopy metrics, like height, canopy area, or canopy volume (McClaran et al. 2013 used a Cylinder).

Recently, Biederman et al. (in prep) collected information about canopy volume characteristics and basal diameters for several transects located around the Santa Rita Experimental Range (SRER)

Biederman et al. report an allometric scaling exponent α = 0.6 for their height regression, while McClaran et al. (2014) reported scaling of height with basal diameter (α = 0.52)

Notes

Figure: The variation in size between diameter measured at or below root crown near ground level and the equivalent diameter (the sum of all basal areas for each minor stem converted to diameter). The dashed black line is a 1:1 reference point. The solid line is a least square regression fit to an α equal to 1.

Figure: The canopy scaling of individual trees' canopy height. Data are from the SRER collected by Biederman and Naito. The dashed black line is a 1:1 reference point. The solid line is a least square regression fit to an α equal to 1.

Joel has already collected:

1) The belt transects were 60 m long radii from the tower in 8 directions and represent about 9% of the area of a circle with the same radius.

Based on both Joel's and Adam's measurements, the canopy diameter (CD) versus max tree height (HT) of individual mesquites is approximately: CD=HT*1.41 - for the segmentation I set the parameter to 1.5; the minimum height cut-off is 1 foot.

tile19s15e06=asc_import('E:\PAG2011\CHM\ASC\19S15E06_LDRY11.ASC');
tile_v=vlm(tile19s15e06,1.5,2);
vlm_ex=export_utm(tile_v,tile19s15e06);
col_header={'ID','UTME','UTMN','HT','PRED','AREA','EQDIAM','MAJAX','MINAX','MAXHT','MINHT','MEANHT'};
xlswrite('G:\SRER\STEMS\19s15e06.xlsx',vlm_ex,1,'A2');
xlswrite('G:\SRER\STEMS\19s15e06.xlsx',col_header,1,'A1');

Figure: Canopy height model in black to white color ramp. Naito field measured stems are yellow and VLM isolation stems are orange.

Next I generated voronoi polygons in QGIS using the r.voronoi module. The voronoi polygon was the first step in segmenting the individual point clouds using FUSION PolyClipData.

Figure: Canopy height model in black to white color ramp. Naito field measured stems are red, VLM isolation are green, and voronoi polygons are visible. Tree heights are given in units of meters.

Created a FUSION Batch script to segment the original LAS data into polygon segments. The PolyClipData tool can only take in ~32k polygons at a time, so I clipped the extent of stems for analysis to the area around where Adam sampled.

PolyClipData /multifile /shape:1,* D:\mesquite_allometry\voronoi.shp D:\mesquite_allometry\LAS\stems.las D:\mesquite_allometry\LAS\19S15E06_LDRY11.las

PolyClip was able to segment out the 7165 stems in the local area, but it took quite a long time (96 minutes) to do a partial tile.

After running PolyClip I ran the ClipData function to further isolate individual trees from within the voronoi footprint. This extra step was done to ensure that (1) trees which share canopies are split by the voronoi polygon, (2) footprints that are larger than the expected canopy area are further reduced, and (3) to normalize the elevation of the points into units of height. The ClipData function was run individually for each tree after I created a spreadsheet in Excel and copy pasted the files into a text editor (TextPad).

ClipData /shape:1 /dtm:D:\mesquite_allometry\19S15E06_LDRY11.dtm /height  D:\mesquite_allometry\LAS\stems_1455.las D:\mesquite_allometry\LAS_clipped_data\tree_1455.las 1029635.002 298975.0019 1029656.998 298996.9981

Estimating the volume of a parabolic ellipsoid cone.

The filled volume space of each mesquite includes its total canopy area, which can be simulated as a half parabolic ellipsoid. However, the angle of the branches that reach up from the ground level create an empty space beneath each mesquite (estimated angle of 30-45 degrees). This results in a leafy volume space that is only a part of the total ellipsoid, while they underlying volume is a cone.

%Volume of a 1/2 ellipsoid
V=((4/3)*pi*x*y*z1)/2
 
%Volume of a cone
V=1/3*pi*x*y*z2

Identifier	DataFile	FileTitle	Total return count	Total return count above 1.00	Return 1 count above 1.00	Return 2 count above 1.00	Return 3 count above 1.00	Return 4 count above 1.00	Return 5 count above 1.00	Return 6 count above 1.00	Return 7 count above 1.00	Return 8 count above 1.00	Return 9 count above 1.00	Other return count above 1.00	Elev minimum	Elev maximum	Elev mean	Elev mode

Identifier	DataFile	FileTitle	Total return count	Total return count above 1.00	Return 1 count above 1.00	Return 2 count above 1.00	Return 3 count above 1.00	Return 4 count above 1.00	Return 5 count above 1.00	Return 6 count above 1.00	Return 7 count above 1.00	Return 8 count above 1.00	Return 9 count above 1.00	Other return count above 1.00	Elev minimum	Elev maximum	Elev mean	Elev mode