An enhanced approach for surface flow routing over drainage‐constrained triangulated irregular networks

The accuracy and efficiency of the simulations in distributed hydrological models must depend on the proper estimation of flow directions and paths. Numerous studies have been carried out to delineate the drainage patterns based on gridded digital elevation models (DEMs). The triangulated irregular network (TIN) has been increasingly applied in hydrological applications due to the advantages of high storage efficiency and multi‐scale adaptive performance. Much of the previous literature focuses mainly on filling the depressions on gridded DEMs rather than treating the special cases in TIN structures, which has hampered its applications to hydrological models. This study proposes a triangulation‐based solution for the removal of flat areas and pits to enhance the simulation of flow routing on triangulated facet networks. Based on the drainage‐constrained TIN generated from only a gridded DEM by the compound point extraction (CPE) method, the inconsistent situations including flat triangles, V‐shape flat edges and sink nodes are respectively identified and rectified. The optimization algorithm is an iterative process of TIN reconstruction, in which the flat areas are generalized into their center points and the pits are rectified by embedding break lines. To verify the proposed algorithm and investigate the potential for flow routing, flow paths of steepest descent are derived by the vector‐based tracking algorithm based on the optimized TIN. A case study of TIN optimization and flow path tracking was performed on a real‐world DEM. The outcomes indicate that the proposed approach can effectively solve the problem of inconsistencies without a significant loss in accuracy of the terrain model.     

Mapping hourly dynamics of urban population using trajectories reconstructed from mobile phone records

Understanding the spatiotemporal dynamics of urban population is crucial for addressing a wide range of urban planning and management issues. Aggregated geospatial big data have been widely used to quantitatively estimate population distribution at fine spatial scales over a given time period. However, it is still a challenge to estimate population density at a fine temporal resolution over a large geographical space, mainly due to the temporal asynchrony of population movement and the challenges to acquiring a complete individual movement record. In this article, we propose a method to estimate hourly population density by examining the time‐series individual trajectories, which were reconstructed from call detail records using BP neural networks. We first used BP neural networks to predict the positions of mobile phone users at an hourly interval and then estimated the hourly population density using log‐linear regression at the cell tower level. The estimated population density is linearly correlated with population census data at the sub‐district level. Trajectory clustering results show five distinct diurnal dynamic patterns of population movement in the study area, revealing spatially explicit characteristics of the diurnal commuting flows, though the driving forces of the flows need further investigation.     

Extraction of Alteration Mineral Information from Moderate Remote Sensing Images Using MPS Method

Although alteration minerals related to metallogenesis is very important in mineral exploration, information of alteration mineral is weakly expressed in remote sensing imagery, which is often subject to interfering noise and sometimes limited in spectral and spatial resolutions. Because of easy access, moderate images are the main sources of alteration mineral information. Therefore, it is very important to develop alteration mineral information extraction methods from remote sensing images. In this paper, a combined method based on Mask, principal component analysis (PCA) and support vector machine method (SVM) was used to extract alteration mineral information from Enhanced thematic mapper plus remote sensing data with limited spectral and spatial resolutions. First, a mask image of the remote sensing imagery was created to remove interference information such as vegetation, shadow and water. Then, PCA was employed to collect sample data relating to iron, argillic, and carbonatization alteration. Finally, SVM was used to deal with alteration anomaly and build a feature extraction model of high accuracy. The Mask-PCA-SVM model is used to extract alteration mineral information from remote sensing images of Hatu area, Xinjiang Uygur Autonomous Regions, China. The results show that the new methods proposed in this paper can coincide well with known deposits occurrences, rate reached 86.51%. While, the consistent rate with known deposits of the ratio model, PCA model and Spectral angle mapper model were only 3.37, 65.08 and 69.05% respectively. This suggests that the proposed model can find the actual distribution of mineral deposits more effectively by reducing interference to a greater degree.     

Analysis of the Seasonal Variation of CO<Subscript>2</Subscript> Concentration in China Based on GOSAT Satellite Data

Carbon dioxide (CO₂) is one of the major gases that contribute to the global warming. Therefore, studying the distribution of CO₂ can help people understand the carbon cycle. Based on the GOSAT retrieved CO₂ products, the temporal and spatial distribution and seasonal variation of CO₂ concentration were analyzed from 2011 to 2015. CO₂ concentration has obvious seasonal variation. It was low in summer, and was high in spring, and the annual increase was about 2 ppm. Nevertheless, the annual growth rate of CO₂ concentration in summer was higher than that in spring, it was 0.5425% in summer and was 0.46% in spring. CO₂ concentration was low in the northwest and was high in the southeast. The growth rate of CO₂ was 2.8 ppm in the northwest and was 3.42 ppm in the southeast. More human’s activities made CO₂ concentration higher in the southeast than that in other regions.     

Study on Sampling Rule and Simplification of LiDAR Point Cloud Based on Terrain Complexity

It is difficult to obtain digital elevation model (DEM) in the mountainous regions. As an emerging technology, Light Detection and Ranging (LiDAR) is an enabling technology. However, the amount of points obtained by LiDAR is huge. When processing LiDAR point cloud, huge data will lead to a rapid decline in data processing speed, so it is necessary to thin LiDAR point cloud. In this paper, a new terrain sampling rule had been built based on the integrated terrain complexity, and then based on the rule a LiDAR point cloud simplification method, which was referred as to TCthin, had been proposed. The TCthin method was evaluated by experiments in which XUthin and Lasthin were selected as the TCthin’s comparative methods. The TCthin’s simplification degree was estimated by the simplification rate value, and the TCthin’s simplification quality was evaluated by Root Mean Square Deviation. The experimental results show that the TCthin method can thin LiDAR point cloud effectively and improve the simplification quality, and at 5 m, 10 m, 30 m scale levels, the TCthin method has a good applicability in the areas with different terrain complexity. This study has theoretical and practical value in sampling theory, thinning LiDAR point cloud, building high-precision DEM and so on.     

A grid-based tropospheric product for China using a GNSS network

Tropospheric delay accounts for one source of error in global navigation satellite systems (GNSS). To better characterize the tropospheric delays in the temporal and spatial domain and facilitate the safety-critical use of GNSS across China, a method is proposed to generate a grid-based tropospheric product (GTP) using the GNSS network with an empirical tropospheric model, known as IGGtrop. The prototype system generates the GTPs in post-processing and real-time modes and is based on the undifferenced and uncombined precise point positioning (UU-PPP) technique. GTPs are constructed for a grid form (\(2.0{^{\circ }}\times 2.5{^{\circ }}\) latitude–longitude) over China with a time resolution of 5 min. The real-time GTP messages are encoded in a self-defined RTCM3 format and broadcast to users using NTRIP (networked transport of RTCM via internet protocol), which enables efficient and safe transmission to real-time users. Our approach for GTP generation consists of three sequential steps. In the first step, GNSS-derived zenith tropospheric delays (ZTDs) for a network of GNSS stations are estimated using UU-PPP. In the second step, vertical adjustments for the GNSS-derived ZTDs are applied to address the height differences between the GNSS stations and grid points. The ZTD height corrections are provided by the IGGtrop model. Finally, an inverse distance weighting method is used to interpolate the GNSS-derived ZTDs from the surrounding GNSS stations to the location of the grid point. A total of 210 global positioning system (GPS) stations from the crustal movement observation network of China are used to generate the GTPs in both post-processing and real-time modes. The accuracies of the GTPs are assessed against with ERA-Interim-derived ZTDs and the GPS-derived ZTDs at 12 test GPS stations, respectively. The results show that the post-processing and real-time GTPs can provide the ZTDs with accuracies of 1.4 and 1.8 cm, respectively. We also apply the GTPs in real-time kinematic GPS PPP, and the results show that the convergence time of the PPP solutions is shortened. These results confirm that the GTPs can act as an efficient information source to augment GNSS positioning over China.     

Weak GPS signal tracking using FFT discriminator in open loop receiver

An open loop tracking architecture, which tracks GPS signals under weak and challenging conditions, is analyzed. The in-phase and quadrature-phase integration pair is regarded as a single tone complex signal. An FFT-based method is used as a frequency discriminator to estimate the Doppler frequency residual of the single tone signal. Another FFT-based method applies complex squaring to eliminate the effect of the navigation data bits polarities. The performance of the FFT-based discriminators is assessed in three criteria. Those criteria are the signal strength and dynamic range that can be tracked and the accuracy of the estimated Doppler frequency. In addition, the performance of the discriminators is analyzed to provide the theoretical and simulated peak detection probability. The results indicate that the FFT discriminator can track signals about 5 dB weaker than the signals that can be tracked by the complex squared FFT discriminator. In a quasi-static environment, the Doppler frequency residual can be assumed to be around zero, which can enable the FFT-based discriminators to track signals with approximately 2 dB less power. Moreover, the performance of the FFT-based discriminators is compared with the performance of two other frequency discriminators, namely the fast–slow and power-based. The comparison results indicate that these two frequency discriminators give higher frequency estimation accuracy, but they have a narrower dynamic range.     

Generating GPS satellite fractional cycle bias for ambiguity-fixed precise point positioning

With the development of precise point positioning (PPP), the School of Geodesy and Geomatics (SGG) at Wuhan University is now routinely producing GPS satellite fractional cycle bias (FCB) products with open access for worldwide PPP users to conduct ambiguity-fixed PPP solution. We provide a brief theoretical background of PPP and present the strategies and models to compute the FCB products. The practical realization of the two-step (wide-lane and narrow-lane) FCB estimation scheme is described in detail. With GPS measurements taken in various situations, i.e., static, dynamic, and on low earth orbit (LEO) satellites, the quality of FCB estimation and the effectiveness of PPP ambiguity resolution (AR) are evaluated. The comparison with CNES FCBs indicated that our FCBs had a good consistency with the CNES ones. For wide-lane FCB, almost all the differences of the two products were within ±0.05 cycles. For narrow-lane FCB, 87.8 % of the differences were located between ±0.05 cycles, and 97.4 % of them were located between ±0.075 cycles. The experimental results showed that, compared with conventional ambiguity-float PPP, the averaged position RMS of static PPP can be improved from (3.6, 1.4, 3.6) to (2.0, 1.0, 2.7) centimeters for ambiguity-fixed PPP. The average accuracy improvement in the east, north, and up components reached 44.4, 28.6, and 25.0 %, respectively. A kinematic, ambiguity-fixed PPP test with observation of 80 min achieved a position accuracy of better than 5 cm at the one-sigma level in all three coordinate components. Compared with the results of ambiguity-float, kinematic PPP, the positioning biases of ambiguity-fixed PPP were improved by about 78.2, 20.8, and 65.1 % in east, north, and up. The RMS of LEO PPP test was improved by about 23.0, 37.0, and 43.0 % for GRACE-A and GRACE-B in radial, tangential, and normal directions when AR was applied to the same data set. These results demonstrated that the SGG FCB products can be produced with high quality for users anywhere around the world to carry out ambiguity-fixed PPP solutions.     

Assessing Performance of Tomo-SAR and Backscattering Coefficient for Hemi-Boreal Forest Aboveground Biomass Estimation

Tomo-SAR technique has been used for hemi-boreal forest height and further forest biomass estimation through allometric equation. Backscattering coefficient especially in longer wavelength (L- or P-band) is thought as a useful parameter for hemi-boreal forest biomass retrieval. The aim of this paper is to assess the performance of vertical backscattering power and backscattering coefficient for hemi-boreal forest aboveground biomass (AGB) estimation with airborne P-band data. The test site locates in southern Sweden called Remningstorp test site, and the in-situ forest AGB ranges from 14 t/ha to 245 t/ha at stand level. Multi-baseline P-band Pol-InSAR data in repeat-path mode collected during March and May in 2007 at Remningstorp test site was used. We found that the correlation coefficient (R) between backscattering coefficient of P-band HH polarization and the in-situ forest biomass reached 0.87. The R for P-band VV backscattering power at 5 m is 0.71 and 10 m is 0.72. Backscattering coefficient in HH polarization and vertical backscattering power at 5 m and 10 m were applied to construct a model for hemi-boreal forest AGB estimation by backward step-wise regression and cross-validation approach. The results showed that the estimated forest AGB ranges from 19 to 240 t/ha, and the constructed model obtained a higher R and smaller RMSE, the value of R is 0.91, RMSE is 30.43 t/ha at Remningstorp test site.     

Curve matching approaches to waveform classification: a case study using ICESat data

Light Detection and Ranging (LiDAR) waveforms are being increasingly used in many forest and urban applications, especially for ground feature classification. However, most studies relied on either discretizing waveforms to multiple returns or extracting shape metrics from waveforms. The direct use of the full waveform, which contains the most comprehensive and accurate information has been scarcely explored. We proposed to utilize the complete waveform to test its ability to differentiate between objects having distinct vertical structures using curve matching approaches. Two groups of curve matching approaches were developed by extending methods originally designed for pixel-based hyperspectral image classification and object-based high spatial image classification. The first group is based on measuring the curve similarity between an unknown waveform and a reference waveform, including curve root sum squared differential area (CRSSDA), curve angle mapper (CAM), and Kullback–Leibler (KL) divergence. The second group assesses the curve similarity between an unknown and reference cumulative distribution functions (CDFs) of their waveforms, including cumulative curve root sum squared differential area (CCRSSDA), cumulative curve angle mapper (CCAM), and Kolmogorov–Smirnov (KS) distance. When employed to classify open space, trees, and buildings using ICESat waveform data, KL provided the highest average classification accuracy (87%), closely followed by CCRSSDA and CCAM, and they all significantly outperformed KS, CRSSDA, and CAM based on 15 randomized sample sets.