Search

In this study, for feature extraction of seismic- and nonlinear trend terms of the nonlinear and nonstationary high-rate Global Navigation Satellite System (GNSS) seismic displacements, we present an adaptive complete ensemble empirical mode decomposition (CEEMD) method based on correlation coefficient. We also introduce an improved CEEMD (ICEEMD) based on multiway principal component analysis (MPCA), hereafter, ICEEMD-MPCA, denoising method for boosting the high-rate GNSS seismic displacements with/without nonlinear trend. The results show that the proposed method is suitable for extracting features of nonlinear and nonstationary high-rate GNSS seismic displacements. The ICEEMD-MPCA denoising method is able to obviously eliminate white noise at high frequencies, remove systematic errors, and effectively preserve the earthquake wave signals. The denoising method reflects that high-rate GNSS has millimeter-level accuracy in the vertical component and submillimeter-level accuracy in the horizontal component. To verify the performance characteristics of the feature extraction and denoising method, the focal mechanism inversion for the 2010 El Mayor–Cucapah earthquake is estimated using the high-rate GNSS data (5 Hz) and collocated strong motion data (50 and 200 Hz). The focal mechanism of the earthquake estimated using the denoised high-rate GNSS with/without nonlinear trend are more consistent with the teleseismic estimates than those estimated using the raw high-rate GNSS and strong motion data. For the GNSS observations, the inversion results are more affected by errors in low-frequency, which accumulate over time and display their signatures as trend in coseismic displacements. It can be concluded that with the adaptive CEEMD method and the ICEEMD-MPCA denoising method, high-rate GNSS at near-field stations only or combined with GNSS/strong-motion records can be applied to yield better information regarding the earthquake focal mechanism inversion, the lower bounds of earthquake magnitudes, and earthquake rupture parameters for small to moderate earthquakes.

We propose a fundamental theorem for eco-environmental surface modelling (FTEEM) in order to apply it into the fields of ecology and environmental science more easily after the fundamental theorem for Earth’s surface system modeling (FTESM). The Beijing-Tianjin-Hebei (BTH) region is taken as a case area to conduct empirical studies of algorithms for spatial upscaling, spatial downscaling, spatial interpolation, data fusion and model-data assimilation, which are based on high accuracy surface modelling (HASM), corresponding with corollaries of FTEEM. The case studies demonstrate how eco-environmental surface modelling is substantially improved when both extrinsic and intrinsic information are used along with an appropriate method of HASM. Compared with classic algorithms, the HASM-based algorithm for spatial upscaling reduced the root-mean-square error of the BTH elevation surface by 9 m. The HASM-based algorithm for spatial downscaling reduced the relative error of future scenarios of annual mean temperature by 16%. The HASM-based algorithm for spatial interpolation reduced the relative error of change trend of annual mean precipitation by 0.2%. The HASM-based algorithm for data fusion reduced the relative error of change trend of annual mean temperature by 70%. The HASM-based algorithm for model-data assimilation reduced the relative error of carbon stocks by 40%. We propose five theoretical challenges and three application problems of HASM that need to be addressed to improve FTEEM.

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 high-density 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.

To attenuate positional errors of LiDAR-derived datasets for constructing digital elevation models (DEMs), a feature-preserving point denoising algorithm (F-PDA) is developed in this paper. F-PDA i...

High resolution and accurate precipitation data is significantly important for numerous hydrological applications. To enhance the spatial resolution and accuracy of satellite-based precipitation products, an easy-to-use downscaling-calibration method based on spatial Random Forest (SRF) is proposed in this paper, where the spatial autocorrelation between precipitation measurements is taken into account. The proposed method consists of two main stages. Firstly, the satellite-based precipitation was downscaled by SRF with the incorporation of some high-resolution covariates including latitude, longitude, DEM, NDVI, terrain slope, aspect, relief, and land surface temperatures. Then, the downscaled precipitation was calibrated by SRF with rain gauge observations and the aforementioned high-resolution variables. The monthly Integrated MultisatellitE Retrievals for Global Precipitation Measurement (IMERG) located in Sichuan province, China from 2015 to 2019 was processed using our method and its results were compared with those of some 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). Results show that: (1) the proposed method outperforms the other methods as well as the original IMERG; (2) the monthly-based SRF estimation is slightly more accurate than the annual-based SRF fraction disaggregation method; (3) SRF-based downscaling and calibration preforms better than bilinear downscaling (Bi-SRF) and GDA-based calibration (SRF-GDA); (4) kriging seems more accurate than GWR and ANN in terms of quantitative accuracy measures, whereas its precipitation map cannot capture the detailed spatial precipitation patterns; and (5) among the predictors for calibration, the precipitation interpolated by kriging on the gauge measurements is the most important variable, indicating the significance for the inclusion of spatial autocorrelation information in gauge measurements.

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.

Airborne light detection and ranging (lidar) is becoming a widely adopted technique for capturing elevation data, which are mainly used for creating digital terrain models (DTMs). However, the large size of lidar datasets poses a serious computational challenge to the promising radial basis function (RBF) interpolation method. In this work, to reduce the huge computational cost and improve the interpolation accuracy, random Fourier features are first introduced to approximate the Gaussian kernel of RBFs in feature space, then a random features-based weighted RBF interpolation method is developed. Based on randomized Fourier features, the nonlinear kernel-based training and evaluation of the RBF method is transformed into simple linear operations in feature space, and with the help of weighted ridge regression, the negative effect of the non-Gaussian distribution of lidar datasets on DTM production is reduced. In other words, the combination of randomized Fourier features and weighted ridge regression improves the efficiency and accuracy of the RBF interpolation method. Experiments on simulated datasets indicate that the proposed method performs better than the classical or random features-based RBF method for dealing with non-Gaussian distributed samples, with the former being slightly less accurate than the iterative RBF method due to the low-dimensional random features. However, the computational cost of the new method is much lower compared with the classical or iterative RBFs. Interpolation of airborne lidar-derived points demonstrates that the new method has a computational cost similar to the inverse distance weighting and triangulated irregular network (TIN) approaches, and is significantly faster than the ordinary kriging (OK) or thin plate spline (TPS) methods. Quantitatively, for interpolation of 644,433 points, the proposed method is approximately 833 and 21 times faster than OK and TPS, respectively. Moreover, the new method avoids the surface discontinuity artifacts presented by the OK, TPS, and TIN methods.

In order to improve the possibility of correct interpretation of some geophysical phenomena such seismology, it is important to reduce systematic errors such as common mode error (CME) and multipath errors, and to improve the accuracy of high-rate (1–5 Hz) displacements. This paper focuses on accuracy enhancement and noise assessment of high-rate Global Navigation Satellite System (GNSS) position time series of regional network of GNSS stations. Complete ensemble empirical mode decomposition-based multiscale multiway principal component analysis (C-MSMPCA) denoising method is introduced to denoise the high-rate (1–5 Hz) GNSS position time series, whose performance was compared with those of modified sidereal filtering (MSF), stacking filtering, PCA and the recently developed multiscale multiway PCA (MSMPCA). It turns out that C-MSMPCA can significantly eliminate the high-frequency random noise and remove the low-frequency CME, and other systematic errors such as multipath errors. Furthermore, for denoising high-rate GNSS time series, C-MSMPCA is more precise than the other denoising methods. Irrespective of sampling rates and data processing strategies, C-MSMPCA can respectively reach the accuracy of submillimeter- and millimeter-levels on the horizontal and vertical components of high-rate GNSS displacements. Spectral analysis on the 1 Hz coordinate time series from the Sichuan GNSS Continuous Observing Network of China suggests that Gauss-Markov process is an adequate noise model for the noise characteristics of the filtered GNSS data using MSMPCA and C-MSMPCA. The CME magnitudes of MSMPCA are 2.61, 2.55, and 6.72 mm (east, north, and vertical), and those of C-MSMPCA are 2.69, 2.80, and 7.80 mm, respectively.

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.

For many geoscience applications, the quality of digital elevation models (DEMs) is the first and foremost requirement. DEM quality generally depends on the interpolation procedure, especially in forested areas with complex terrain. However, classical interpolation methods do not always preserve terrain features, owing to their isotropic and local estimation nature. Therefore, this paper presents a feature-preserving interpolation method, where a structure tensor is integrated into the radial basis function to ensure its anisotropy. This method alleviates the influence of points located on different surface patches. The performance of the proposed method was compared with those of classical interpolation methods, including the inverse distance weighting (IDW), ordinary kriging (OK), topo to raster (ANUDEM), and natural neighbour (NN) methods. The performance of these methods was tested on benchmark data provided by the International Society for Photogrammetry and Remote Sensing (ISPRS) Commission, and on one private dataset. Both datasets were collected via the airborne light detection and ranging (LiDAR) technique. The results from both datasets demonstrate that from a quantitative perspective, the proposed method produces more accurate DEMs than the classical interpolation methods. Specifically, in terms of root mean square error (RMSE), the proposed method is at least 15.8% and 7.6% more accurate when applied to the ISPRS and private data, respectively. Moreover, the proposed method produces more visually appealing surfaces with a good trade-off between terrain feature retention and noise removal.

Surface modeling with very large data sets is challenging. An efficient method for modeling massive data sets using the high accuracy surface modeling method (HASM) is proposed, and HASM_Big is developed to handle very large data sets. A large data set is defined here as a large spatial domain with high resolution leading to a linear equation with matrix dimensions of hundreds of thousands. An augmented system approach is employed to solve the equality-constrained least squares problem (LSE) produced in HASM_Big, and a block row action method is applied to solve the corresponding very large matrix equations. A matrix partitioning method is used to avoid information redundancy among each block and thereby accelerate the model. Experiments including numerical tests and real-world applications are used to compare the performances of HASM_Big with its previous version, HASM. Results show that the memory storage and computing speed of HASM_Big are better than those of HASM. It is found that the computational cost of HASM_Big is linearly scalable, even with massive data sets. In conclusion, HASM_Big provides a powerful tool for surface modeling, especially when there are millions or more computing grid cells.

Due to the penetration ability of airborne light detection and ranging (lidar) into tree crowns, data pits commonly appear in lidar-derived canopy height models (CHMs). They have a seriously negative effect on the quality of tree detection and subsequent biophysical measurements. To construct a pit-free CHM, an algorithm based on robust locally weighted regression and robust z-score was presented to remove data pits. The significant advantage of the new algorithm is parameter-free, which makes it efficient and robust for practical applications. A numerical test and a real-world example were respectively employed to assess the performance of our method for CHM construction, and its results were compared with those of three classical methods including natural neighbor interpolation of the highest point method, mean and median filters. The numerical test demonstrates that our algorithm is more accurate than the other methods for generating pit-free CHMs under the presence of data pits. The real-world example shows that compared with the classical methods, our method has a better ability of data pit removal. Moreover, our method performs better than the other methods for deriving plot-level maximum tree height from CHMs. In a word, the new method shows high potential for pit-free CHM construction.

To suppress the influence of outliers on function estimation, we propose a least absolute deviation (LAD)-based robust support vector regression (SVR). Furthermore, an efficient algorithm based on the split-Bregman iteration is introduced to solve the optimization problem of the proposed algorithm. Both artificial and benchmark datasets are employed to compare the performance of the proposed algorithm with those of least squares SVR (LS-SVR), and two weighted versions of LS-SVR with the weight functions of Hampel and Logistic, respectively. Experiments demonstrate the superiority of the proposed algorithm.

Inaccurate extrinsic parameters are one of the major errors of a mobile LiDAR system (MLS). As different types of sensors with incommensurable precisions are integrated in an MLS, the extrinsic parameters cannot be easily isolated and estimated, especially in an automatic mode. To address this issue, this paper proposes an automatic extrinsic self-calibration method for an MLS based on planar and spherical features. First, the planar and spherical features are automatically extracted from scanned point cloud data of different strips using the Random Sample Consensus algorithm, and the corresponding features are matched after the extraction. Secondly, a rigorous relationship is established between the direct geo-referencing equation, the extrinsic parameters, and the geometric constraint model. Thirdly, the extrinsic parameters are calibrated by minimizing the sum of the squares of distances from points on the feature surface to the matched reference feature. Four datasets collected by using four types of MLS were used to verify the proposed method. As a result of combining two geometric features into a single self-calibration adjustment, the experimental results show that the proposed method is superior to the conventional plane-based method in terms of the positional accuracy. The standard deviation of the distance between the check features of the four datasets collected by several mobile platforms before (after) extrinsic calibration were 0.051 m (0.024 m), 0.060 m (0.018 m), 0.029 m (0.009 m), and 0.354 m (0.070 m), which demonstrated the high compatibility and practicality of the proposed method.

Separating ground from non-ground points is a challenging and essential task before the applications of airborne Light Detection And Ranging (LiDAR) data. The classical filters generally show good results on simple landscapes, but suffer from large filtering errors on complex landscapes with steep slopes and terrain discontinuities. To produce more satisfactory results over these landscapes, a scale-irrelevant and terrain-adaptive interpolation-based filter is presented in this paper. The contributions of the proposed method include a 1D spline-based algorithm for the collection of evenly distributed ground seeds as many as possible, a scale-irrelevant interpolation for estimating the heights of unclassified points and a terrain-adaptive elevation threshold to adapt to various terrain characteristics. The performance of the proposed method was first evaluated on the International Society Photogrammetry and Remote Sensing (ISPRS) benchmark dataset. Results show that the proposed method with the average total error of 2.70% and kappa coefficient of 90.84% outperforms the existing filtering algorithms developed in recent decade (2010–2020). The performance of the proposed method was further assessed on four high-density airborne LiDAR point clouds located in urban and forested sites and compared with four state-of-the-art filters including progressive morphological filter (PMF), cloth simulation filter (CSF), progressive TIN densification (PTD) and multiresolution hierarchical filter (MHF). Results demonstrate that the proposed method is averagely more accurate than the well-known filters in terms of total error and kappa coefficient, and much faster than PTD and MHF. Moreover, the proposed method produces more satisfactory DEMs than the classical methods.