Forest canopy height is an essential parameter in estimating forest aboveground biomass (AGB), growing stock volume (GSV), and carbon storage, and it can provide necessary information in forest management activities. Light direction and ranging (LiDAR) is widely used for estimating canopy height. Considering the high cost of acquiring LiDAR data over large areas, we took a two-stage up-scaling approach in estimating forest canopy height and aimed to develop a method for quantifying the uncertainty of the estimation result. Based on the generalized hierarchical model-based (GHMB) estimation framework, a new estimation framework named RK-GHMB that makes use of a geostatistical method (regression kriging, RK) was developed. In this framework, the wall-to-wall forest canopy height and corresponding uncertainty in map unit scale are generated. This study was carried out by integrating plot data, sampled airborne LiDAR data, and wall-to-wall Ziyuan-3 satellite (ZY3) stereo images. The result shows that RK-GHMB can obtain a similar estimation accuracy (r = 0.92, MAE = 1.50 m) to GHMB (r = 0.92, MAE = 1.52 m) with plot-based reference data. For LiDAR-based reference data, the accuracy of RK-GHMB (r = 0.78, MAE = 1.75 m) is higher than that of GHMB (r = 0.75, MAE = 1.85 m). The uncertainties for all map units range from 1.54 to 3.60 m for the RK-GHMB results. The values change between 1.84 and 3.60 m for GHMB. This study demonstrates that this two-stage up-scaling approach can be used to monitor forest canopy height. The proposed RK-GHMB approach considers the spatial autocorrelation of neighboring data in the second modeling stage and can achieve a higher accuracy.

Accurate identification of sandy land plays an important role in sandy land prevention and control. It is difficult to identify the nature of sandy land due to vegetation covering the soil in the sandy area. Therefore, HJ-2A hyperspectral data and GF-3 Synthetic Aperture Radar (SAR) data were used as the main data sources in this article. The advantages of the spectral characteristics of a hyperspectral image and the penetration characteristics of SAR data were used synthetically to carry out mixed-pixel decomposition in the “horizontal” direction and polarization decomposition in the “vertical” direction. The results showed that in the study area of the Otingdag Sandy Land, in China, the accuracy of sandy land detection based on feature-level fusion and single GF-3 data was verified to be 92% in both cases by field data; the accuracy of sandy land detection based on feature-level fusion was verified to be 88.74% by the data collected from Google high-resolution imagery, which was higher than that based on single HJ-2A (74.17%) and single GF-3 data (88.08%). To further verify the universality of the feature-level fusion method for sandy land detection, Alxa sandy land was also used as a verification area and the accuracy of sandy land detection was verified to be as high as 88.74%. The method proposed in this paper made full use of the horizontal and vertical structural information of remote sensing data. The problem of mixed pixels in sparse-vegetation scenes in the horizontal direction and the problem of vegetation covering sandy soil in the vertical direction were both well solved. Accurate identification of sandy land can be realized effectively, which can provide technical support for sandy land prevention and control.

Airborne LiDAR (ALS) data have been extensively utilized for aboveground biomass (AGB) estimation; however, the high acquisition costs make it challenging to attain wall-to-wall estimation across large regions. Some studies have leveraged ALS data as intermediate variables to amplify sample sizes, thereby reducing costs and enhancing sample representativeness and model accuracy, but the cost issue remains in larger-scale estimations. Satellite LiDAR data, offering a broader dataset that can be acquired quickly with lower costs, can serve as an alternative intermediate variable for sample expansion. In this study, we employed a three-stage up-scaling approach to estimate forest AGB and introduced a method for quantifying estimation uncertainty. Based on the established three-stage general-hierarchical-model-based estimation inference (3sGHMB), an RK-3sGHMB inference method is proposed to make use of the regression-kriging (RK) method, and then it is compared with conventional model-based inference (CMB), general hierarchical model-based inference (GHMB), and improved general hierarchical model-based inference (RK-GHMB) to estimate forest AGB and uncertainty at both the pixel and forest farm levels. This study was carried out by integrating plot data, sampled ALS data, wall-to-wall Sentinel-2A data, and airborne P-SAR data. The results show that the accuracy of CMB ($R_{adj}^2$ = 0.37, RMSE = 33.95 t/ha, EA = 63.28%) is lower than that of GHMB ($R_{adj}^2$ = 0.38, RMSE = 33.72 t/ha, EA = 63.53%), while it is higher than that of 3sGHMB ($R_{adj}^2$ = 0.27, RMSE = 36.58 t/ha, EA = 60.43%). Notably, RK-GHMB ($R_{adj}^2$ = 0.60, RMSE= 27.07 t/ha, EA = 70.72%) and RK-3sGHMB ($R_{adj}^2$ = 0.55, RMSE = 28.55 t/ha, EA = 69.13%) demonstrate significant accuracy enhancements compared to GHMB and 3sGHMB. For population AGB estimation, the precision of the proposed RK-3sGHMB (p = 94.44%) is the highest, providing that there are sufficient sample sizes in the third stage, followed by RK-GHMB (p = 93.32%) with sufficient sample sizes in the second stage, GHMB (p = 90.88%), 3sGHMB (p = 88.91%), and CMB (p = 87.96%). Further analysis reveals that the three-stage model, considering spatial correlation at the third stage, can improve estimation accuracy, but the prerequisite is that the sample size in the third stage must be sufficient. For large-scale estimation, the RK-3sGHMB model proposed herein offers certain advantages.