Precipitation data with high resolution and high accuracy are significantly important in numerous hydrological applications. To enhance the spatial resolution and accuracy of satellite-based precipitation products, an easy-to-use downscaling-calibration method based on a spatial random forest (SRF-DC) is proposed in this study, where the spatial autocorrelation of precipitation measurements between neighboring locations is considered. SRF-DC consists of two main stages. First, the satellite-based precipitation is downscaled by the SRF with the incorporation of high-resolution variables including latitude, longitude, normalized difference vegetation index (NDVI), digital elevation model (DEM), terrain slope, aspect, relief and land surface temperatures. Then, the downscaled precipitation is calibrated by the SRF with rain gauge observations and the aforementioned high-resolution variables. The monthly Integrated MultisatellitE Retrievals for Global Precipitation Measurement (IMERG) over Sichuan Province, China, from 2015 to 2019 was processed using SRF-DC, and its results were compared with those of classical methods including geographically weighted regression (GWR), artificial neural network (ANN), random forest (RF), kriging interpolation only on gauge measurements, bilinear interpolation-based downscaling and then SRF-based calibration (Bi-SRF), and SRF-based downscaling and then geographical difference analysis (GDA)-based calibration (SRF-GDA). Comparative analyses with respect to root mean square error (RMSE), mean absolute error (MAE) and correlation coefficient (CC) demonstrate that (1) SRF-DC outperforms the classical methods as well as the original IMERG; (2) the monthly based SRF estimation is slightly more accurate than the annually based SRF fraction disaggregation method; (3) SRF-based downscaling and calibration perform better than bilinear downscaling (Bi-SRF) and GDA-based calibration (SRF-GDA); (4) kriging is more accurate than GWR and ANN, whereas its precipitation map loses detailed spatial precipitation patterns; and (5) based on the variable-importance rank of the RF, the precipitation interpolated by kriging on the rain gauge measurements is the most important variable, indicating the significance of incorporating spatial autocorrelation for precipitation estimation.

Airborne light detection and ranging (LiDAR) technology has become the mainstream data source in geosciences and environmental sciences. Point cloud filtering is a prerequisite for almost all LiDAR-based applications. However, it is challenging to select a suitable filtering algorithm for handling high-density point clouds over complex landscapes. Therefore, to determine an appropriate filter on a specific environment, this paper comparatively assessed the performance of five representative filtering algorithms on six study sites with different terrain characteristics, where three plots are located in urban areas and three in forest areas. The representative filtering methods include simple morphological filter (SMRF), multiresolution hierarchical filter (MHF), slope-based filter (SBF), progressive TIN densification (PTD) and segmentation-based filter (SegBF). Results demonstrate that SMRF performs the best in urban areas, and compared to MHF, SBF, PTD and SegBF, the total error of SMRF is reduced by 1.38%, 48.21%, 48.25% and 31.03%, respectively. MHF outperforms the others in forest areas, and compared to SMRF, SBF, PTD and SegBF, the total error of MHF is reduced by 1.98%, 35.87%, 45.11% and 9.42%, respectively. Moreover, both SMRF and MHF keep a good balance between type I and II errors, which makes the produced DEMs much similar to the references. Overall, SMRF and MHF are recommended for urban and forest areas, respectively, and MHF averagely performs slightly better than SMRF on all areas with respect to kappa coefficient.

Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) inherently suffers from various errors. Many previous works employed Geoscience Laser Altimeter System onboard the Ice, Cloud, and land Elevation Satellite (ICESat/GLAS) data to assess and enhance SRTM DEM accuracy. Nevertheless, data coregistration between the two datasets was commonly neglected in their studies. In this paper, an automated and simple three dimensional (3D) coregistration method (3CM) was introduced to align the 3-arc-second SRTM (SRTM3) DEM and ICESat/GLAS data over Jiangxi province, China. Then, accuracy evaluation of the SRTM3 DEM using ICESat/GLAS data with and without data coregistration was performed on different classes of terrain factors and different land uses, with the purpose of evaluating the importance of data coregistration. Results show that after data coregistration, the root mean square error (RMSE) and mean bias of the SRTM3 DEM are reduced by 14.4% and 97.1%, respectively. Without data coregistration, terrain aspects with a sine-like shape are strongly related to SRTM3 DEM errors; nevertheless, this relationship disappears after data coregistration. Among the six land uses, SRTM3 DEM produces the lowest accuracy in forest areas. Finally, by incorporating land uses, terrain factors and ICESat/GLAS data into the correction models, the SRTM3 DEM was enhanced using multiple linear regression (MLR), back propagation neural network (BPNN), generalized regression NN (GRNN), and random forest (RF), respectively. Results exhibit that the four enhancement models with data coregistration obviously outperform themselves without the coregistration. Among the four models, RF produces the best result, and its RMSE is about 3.1%, 2.7% and 11.3% lower than those of MLR, BPNN, and GRNN, respectively. Moreover, 146 Global Navigation Satellite System (GNSS) points over Ganzhou city of Jiangxi province were used to assess the accuracy of the RF-derived SRTM3 DEM. It is found that the DEM quality is improved and has a similar error magnitude to that relative to the ICESat/GLASS data.

Airborne light detection and ranging (LiDAR) datasets with a large volume pose a great challenge to the traditional interpolation methods for the production of digital terrain models (DTMs). Thus, a fast, global interpolation method based on thin plate spline (TPS) is proposed in this paper. In the methodology, a weighted version of finite difference TPS is first developed to deal with the problem of missing data in the grid-based surface construction. Then, the interpolation matrix of the weighted TPS is deduced and found to be largely sparse. Furthermore, the values and positions of each nonzero element in the matrix are analytically determined. Finally, to make full use of the sparseness of the interpolation matrix, the linear system is solved with an iterative manner. These make the new method not only fast, but also require less random-access memory. Tests on six simulated datasets indicate that compared to recently developed discrete cosine transformation (DCT)-based TPS, the proposed method has a higher speed and accuracy, lower memory requirement, and less sensitivity to the smoothing parameter. Real-world examples on 10 public and 1 private dataset demonstrate that compared to the DCT-based TPS and the locally weighted interpolation methods, such as linear, natural neighbor (NN), inverse distance weighting (IDW), and ordinary kriging (OK), the proposed method produces visually good surfaces, which overcome the problems of peak-cutting, coarseness, and discontinuity of the aforementioned interpolators. More importantly, the proposed method has a similar performance to the simple interpolation methods (e.g., IDW and NN) with respect to computing time and memory cost, and significantly outperforms OK. Overall, the proposed method with low memory requirement and computing cost offers great potential for the derivation of DTMs from large-scale LiDAR datasets.

To overcome the huge volume problem of light detection and ranging (LiDAR) data for the derivation of digital terrain models (DTMs), a least squares compactly supported radial basis function (CSRBF) interpolation method is proposed in this paper. The proposed method has a limited support radius and fewer RBF centers than the sample points, selected by a newly developed surface variation-based algorithm. Those make the linear system of the proposed method not only much sparser but also efficiently solvable. Tests on a synthetic dataset demonstrate that the proposed method is comparable to the smoothing RBF, and far superior to the exact RBF. Moreover, the first is much faster than the others. The proposed method with the RBF centers selected by the surface variation-based algorithm obviously outperforms that with the random selection of equal number. Real-world examples on one private and ten public datasets show that the surfaces of simple interpolation methods including inverse distance weighting, natural neighbor, linear and bicubic suffer from the problems of roughness, peak-cutting, discontinuity and subtle terrain feature loss, respectively. By contrast, the proposed method produces visually appealing results, keeping a good tradeoff between noise removal and terrain feature preservation. Additionally, the new method compares favorably with ordinary kriging (OK) for the generation of high-resolution DTMs in terms of interpolation accuracy, yet the former is much more robust to spatial resolution variation and terrain characteristics than the latter. More importantly, our method is about 4 times faster than OK. In conclusion, the proposed method has high potential for the interpolation of a large LiDAR dataset, especially when both interpolation accuracy and computational cost are taken into account.

Over the past decades, plenty of filtering algorithms have been presented to distinguish ground and non-ground points from airborne LiDAR point clouds. However, with the existing methods, it is difficult to derive satisfactory filtering results on rugged terrains with dense vegetation due to the low-level penetration ability of laser pulses. Therefore, a multi-level interpolation-based filter is developed in this paper. The novelty of the algorithm lies within its usage of multi-scale morphological operations and robust z-score to correctly select ground seeds as more as possible, and a terrain-adaptive residual threshold to adapt to various terrain characteristics. Rural samples provided by International Society for Photogrammetry and Remote Sensing (ISPRS) were employed to assess the performance of the proposed method and its results were compared with 15 filtering algorithms developed in recent 5 years (2015-2019). Results show that the proposed method with the optimized parameters produces the best accuracy with the average total error and kappa coefficient of 1.89% and 87.88%, respectively. We further filtered highdensity point clouds in six forested areas with different vegetation covers and terrain slopes. Results demonstrate that the proposed algorithm is more accurate than the well-known filtering methods including morphological-based filter, progressive TIN densification filter (PTD), improved PTD and cloth simulation filter, with the average total error decreased by 26.2%, 19.9%, 3.8% and 40.4%, respectively. Moreover, the DEMs of the proposed method have lower average root mean square errors than the four classical filters. Therefore, the proposed method can be considered as an effective ground filtering algorithm for airborne LiDAR point clouds in forested areas.