Author Contributions
Conceptualization, C.W. and A.J.E.; Funding acquisition, C.W. and M.A.C.; Methodology, C.W. and A.J.E.; Resources, C.W., A.J.E., I.N., M.A.C., S.L. and C.R.H.; Software, C.W.; Supervision, A.J.E.; Validation, C.W. and C.R.H.; Writing—original draft, C.W.; Writing—review and editing, C.W., A.J.E., I.N., M.A.C., S.L., C.R.H., Y.Z., Y.L. and Y.T. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Flowchart of the FPSF-CH and its performance analysis. The RF-CH, RH95, and C-CH are canopy heights derived from RF, GEDI waveform, and the FPSF-CH. R and W are the nearest interpolation dataset and weight for each GEDI data point, respectively. CHM is the canopy height model.
Figure 1.
Flowchart of the FPSF-CH and its performance analysis. The RF-CH, RH95, and C-CH are canopy heights derived from RF, GEDI waveform, and the FPSF-CH. R and W are the nearest interpolation dataset and weight for each GEDI data point, respectively. CHM is the canopy height model.
Figure 2.
Locations of study areas. The black rectangles show the locations of Area 1 and Area 2 with 1° × 1° tiles in the ecoregion domain of the North American eastern forests. Three red rectangles represent the locations of three NEON sites, including the BART, HOPB, and HARV. The tree cover map was downloaded from the Global Land Analysis and Discovery website (
https://glad.umd.edu/dataset, accessed on 27 April 2020). The canopy height models (CHMs) in BART, HOPB, and HARV are at the same scale with data ranging from 0 to 40 m.
Figure 2.
Locations of study areas. The black rectangles show the locations of Area 1 and Area 2 with 1° × 1° tiles in the ecoregion domain of the North American eastern forests. Three red rectangles represent the locations of three NEON sites, including the BART, HOPB, and HARV. The tree cover map was downloaded from the Global Land Analysis and Discovery website (
https://glad.umd.edu/dataset, accessed on 27 April 2020). The canopy height models (CHMs) in BART, HOPB, and HARV are at the same scale with data ranging from 0 to 40 m.
Figure 3.
The relationship between canopy heights of C-CH (A1,B1) and GDAL RH95 (A2,B2) compared to high-quality GEDI validation data in Area 1 (A) and Area 2 (B), respectively.
Figure 3.
The relationship between canopy heights of C-CH (A1,B1) and GDAL RH95 (A2,B2) compared to high-quality GEDI validation data in Area 1 (A) and Area 2 (B), respectively.
Figure 4.
Comparison of canopy height spatial distributions in Area 1 (B1,B2) and Area 2 (B3,B4). (A1,C3) are true color images, while (A2,C4) are the topographic slopes of the enlarged rectangle in Area 1 (C1,C2) and Area 2 (A3,A4), respectively. (B1,B2), and (B3,B4) are C-CH and GDAL RH95 in Area 1 and Area 2, respectively. The true color images are composited with the blue, green, and red bands of Sentinel-2. Masked areas are displayed in white.
Figure 4.
Comparison of canopy height spatial distributions in Area 1 (B1,B2) and Area 2 (B3,B4). (A1,C3) are true color images, while (A2,C4) are the topographic slopes of the enlarged rectangle in Area 1 (C1,C2) and Area 2 (A3,A4), respectively. (B1,B2), and (B3,B4) are C-CH and GDAL RH95 in Area 1 and Area 2, respectively. The true color images are composited with the blue, green, and red bands of Sentinel-2. Masked areas are displayed in white.
Figure 5.
The relationship of canopy heights of C-CH (A1–A3) and GDAL RH95 (B1–B3) compared with NEON RH95 in BART (A1,B1), HARV (A2,B2), and HOPB (A3,B3), respectively.
Figure 5.
The relationship of canopy heights of C-CH (A1–A3) and GDAL RH95 (B1–B3) compared with NEON RH95 in BART (A1,B1), HARV (A2,B2), and HOPB (A3,B3), respectively.
Figure 6.
Comparison of canopy height spatial distributions among NEON CHM, C-CH, and GDAL RH95. The (A1–A3,C1–C3) are the sub-regions of ‘1’, and ‘2’ rectangles in (B1–B3). The (B1–B3) represent the NEON CHM, C-CH, and GDAL RH95, respectively. To be comparative, the canopy heights have the same scale, ranging from 0 to 40 m. Masked areas are displayed in white.
Figure 6.
Comparison of canopy height spatial distributions among NEON CHM, C-CH, and GDAL RH95. The (A1–A3,C1–C3) are the sub-regions of ‘1’, and ‘2’ rectangles in (B1–B3). The (B1–B3) represent the NEON CHM, C-CH, and GDAL RH95, respectively. To be comparative, the canopy heights have the same scale, ranging from 0 to 40 m. Masked areas are displayed in white.
Figure 7.
Canopy height residuals of C-CH (A,B) and GDAL RH95 (C,D) plotted as a function of slope (A1–A3,C1–C3) and tree canopy cover (B1–B3,D1–D3) in BART, HARV, and HOPB. To be readable, the median and average of residuals are shown by the orange line and a green triangle in the box, respectively. The negative median and average values indicate under-and over-estimation, respectively. The upper and lower quartiles are indicated by the end of the lines.
Figure 7.
Canopy height residuals of C-CH (A,B) and GDAL RH95 (C,D) plotted as a function of slope (A1–A3,C1–C3) and tree canopy cover (B1–B3,D1–D3) in BART, HARV, and HOPB. To be readable, the median and average of residuals are shown by the orange line and a green triangle in the box, respectively. The negative median and average values indicate under-and over-estimation, respectively. The upper and lower quartiles are indicated by the end of the lines.
Figure 8.
Comparison of distributions between C-CH (C1,C3) and RF-CH (C2,C4) over space and height histograms. (A1,A3) are true color images, while (A2,A4) are topographic slopes for the enlarged rectangle in Area 1 (B1,B2) and Area 2 (B3,B4). (D1–D4) represents histograms of canopy height of C-CH (D1,D3) and RF-CH (D2,D4) in Area 1 (D1,D2) and Area 2 (D3,D4), respectively. The red circles in (D1–D4) show the obvious difference between the canopy heights of RF-CH and C-CH. The true color images are composited with the blue, green, and red bands of Sentinel-2. Masked areas are displayed in white.
Figure 8.
Comparison of distributions between C-CH (C1,C3) and RF-CH (C2,C4) over space and height histograms. (A1,A3) are true color images, while (A2,A4) are topographic slopes for the enlarged rectangle in Area 1 (B1,B2) and Area 2 (B3,B4). (D1–D4) represents histograms of canopy height of C-CH (D1,D3) and RF-CH (D2,D4) in Area 1 (D1,D2) and Area 2 (D3,D4), respectively. The red circles in (D1–D4) show the obvious difference between the canopy heights of RF-CH and C-CH. The true color images are composited with the blue, green, and red bands of Sentinel-2. Masked areas are displayed in white.
Figure 9.
Densities of GEDI observations over space (A1,A2) and amounts (B1,B2) within a 6 km width window for each object pixel in RF-CH in Area 1 and Area 2, respectively. Masked areas in (A1,A2) are displayed in white.
Figure 9.
Densities of GEDI observations over space (A1,A2) and amounts (B1,B2) within a 6 km width window for each object pixel in RF-CH in Area 1 and Area 2, respectively. Masked areas in (A1,A2) are displayed in white.
Figure 10.
Performance of the FPSF-CH using a different number of good-quality GEDI observations (A) and weighting strategies (B) for calibration. The Wos and Wsen represent weights indicated by the DI difference between an object pixel in RF-CH and the corresponding GEDI data and waveform sensitivity, respectively.
Figure 10.
Performance of the FPSF-CH using a different number of good-quality GEDI observations (A) and weighting strategies (B) for calibration. The Wos and Wsen represent weights indicated by the DI difference between an object pixel in RF-CH and the corresponding GEDI data and waveform sensitivity, respectively.
Table 1.
Geographic characteristics of the BART, HOPB, and HARV sites. The abbreviations of MAT, MAP, and MCH refer to mean annual temperature, mean annual precipitation, and mean canopy height, respectively.
Table 1.
Geographic characteristics of the BART, HOPB, and HARV sites. The abbreviations of MAT, MAP, and MCH refer to mean annual temperature, mean annual precipitation, and mean canopy height, respectively.
Sites | Latitude (°) | Longitude (°) | Mean Elevation (m) | MAT (°C) | MAP (mm) | MCH (m) |
---|
BART | 44.0639 | −71.2874 | 274 | 6.20 | 1325 | 23 |
HARV | 42.5369 | −72.1727 | 348 | 7.40 | 1199 | 26 |
HOPB | 42.4719 | −72.3295 | 203 | 7.90 | 1368 | 26 |
Table 2.
The descriptive statistics of training, validation, and good-quality GEDI waveforms. The N is the number of observations with a mean and standard deviation (SD) of canopy height.
Table 2.
The descriptive statistics of training, validation, and good-quality GEDI waveforms. The N is the number of observations with a mean and standard deviation (SD) of canopy height.
Study Area | Year | Training Data | Validation Data | Good-Quality Data |
---|
N | Mean ± SD (m) | N | Mean ± SD (m) | N | Mean ± SD (m) |
---|
Area 1 | 2019 | 2220 | 21.09 ± 6.64 | 2220 | 20.75 ± 5.59 | 52,425 | 17.86 ± 6.89 |
2020 | 4843 | 20.06 ± 11.34 | 4843 | 19.82 ± 11.54 | 65,853 | 17.26 ± 8.77 |
2021 | 93 | 17.75 ± 4.88 | 93 | 16.59 ± 5.49 | 26,517 | 17.93 ± 7.50 |
Area 2 | 2019 | 2130 | 23.22 ± 5.80 | 2130 | 23.41 ± 5.97 | 30,936 | 20.31 ± 7.99 |
2020 | 3882 | 23.06 ± 5.47 | 3882 | 23.09 ± 5.45 | 71,539 | 20.24 ± 7.95 |
2021 | 2235 | 23.21 ± 5.35 | 2235 | 23.17 ± 4.99 | 41,769 | 20.26 ± 7.82 |
Table 3.
Accuracies of the FPSF-CH for canopy height mapping (C-CH) and GDAL RH95 compared to high-quality GEDI validation data (GEDI RH95). The negative bias represents an underestimation compared to GEDI RH95. The asterisk represents a significant correlation (p < 0.05).
Table 3.
Accuracies of the FPSF-CH for canopy height mapping (C-CH) and GDAL RH95 compared to high-quality GEDI validation data (GEDI RH95). The negative bias represents an underestimation compared to GEDI RH95. The asterisk represents a significant correlation (p < 0.05).
| Study Area | r | Bias (m) | MAE (m) | RMSE (m) | rBias (%) | rRMSE (%) |
---|
C-CH | Area 1 | 0.58 * | 0.61 | 3.48 | 4.47 | 3.06 | 22.34 |
Area 2 | 0.36 * | −0.20 | 2.64 | 3.22 | −0.88 | 13.91 |
GDAL RH95 | Area 1 | 0.55 * | −1.28 | 4.60 | 6.11 | −6.41 | 30.53 |
Area 2 | 0.24 * | −2.72 | 3.94 | 5.17 | −11.74 | 22.32 |
Table 4.
Accuracies of the FPSF-CH for canopy height mapping (C-CH) and GDAL RH95 compared with NEON RH95. The negative bias represents an under-estimation compared with NEON RH95. The asterisk represents a significant correlation (p < 0.05).
Table 4.
Accuracies of the FPSF-CH for canopy height mapping (C-CH) and GDAL RH95 compared with NEON RH95. The negative bias represents an under-estimation compared with NEON RH95. The asterisk represents a significant correlation (p < 0.05).
Metrics | C-CH | GDAL RH95 |
---|
BART | HARV | HOPB | BART | HARV | HOPB |
---|
r | 0.58 * | 0.26 * | 0.16 * | 0.43 * | 0.22 * | 0.22 * |
Bias (m) | −0.84 | −1.62 | −0.04 | −1.75 | −4.45 | −4.12 |
MAE (m) | 2.96 | 3.16 | 2.66 | 3.8 | 5.01 | 4.51 |
RMSE (m) | 3.82 | 4.05 | 3.48 | 5.1 | 6.3 | 5.71 |
rBias (%) | −4.15 | −6.49 | −0.16 | −8.59 | −17.83 | −16.38 |
rRMSE (%) | 18.77 | 16.24 | 13.81 | 25.06 | 25.25 | 22.71 |
Table 5.
Accuracy of RF-CH over the two study areas and three NEON sites compared with high-quality GEDI RH95 and NEON RH95, respectively. The asterisk represents a significant correlation (p < 0.05).
Table 5.
Accuracy of RF-CH over the two study areas and three NEON sites compared with high-quality GEDI RH95 and NEON RH95, respectively. The asterisk represents a significant correlation (p < 0.05).
Datasets | Sites | r | Bias (m) | MAE (m) | RMSE (m) | rBias (%) | rRMSE (%) |
---|
GEDI RH95 | Area 1 | 0.55 * | 0.91 | 3.50 | 4.66 | 4.57 | 23.26 |
Area 2 | 0.24 * | −0.39 | 2.93 | 3.73 | −1.69 | 16.10 |
NEON RH95 | BART | 0.54 * | −0.29 | 3.33 | 4.46 | −1.43 | 21.9 |
HARV | 0.29 * | −2.06 | 3.5 | 4.54 | −8.25 | 18.19 |
HOPB | 0.22 * | −3.01 | 3.8 | 4.79 | −11.95 | 19.05 |
Table 6.
Accuracy comparison of canopy heights between good-quality GEDI data and RF-CH, compared with NEON RH95. r is the Pearson r correlation, while the asterisk represents a significant correlation (p < 0.05).
Table 6.
Accuracy comparison of canopy heights between good-quality GEDI data and RF-CH, compared with NEON RH95. r is the Pearson r correlation, while the asterisk represents a significant correlation (p < 0.05).
Datasets | r | Bias (m) | MAE (m) | RMSE (m) |
---|
Good-quality GEDI data | 0.81 * | −2.11 | 2.63 | 3.55 |
RF-CH | 0.54 * | −1.69 | 3.71 | 4.77 |