Search

Light detection and ranging (lidar)-derived elevation data are commonly subjected to outliers due to the boundaries of occlusions, physical imperfections of sensors, and surface reflectance. Outliers have a serious negative effect on the accuracy of digital elevation models (DEMs). To decrease the impact of outliers on DEM construction, we propose a robust interpolation algorithm of multiquadric (MQ) based on a regularized least absolute deviation (LAD) technique. The objective function of the proposed method includes a regularization-based smoothing term and an LAD-based fitting term, respectively, used to smooth noisy samples and resist the influence of outliers. To solve the objective function of the proposed method, we develop a simple scheme based on the split-Bregman iteration algorithm. Results from simulated data sets indicate that when sample points are noisy or contaminated by outliers, the proposed method is more accurate than the classical MQ and two recently developed robust algorithms of MQ for surface modeling. Real-world examples of interpolating 1 private and 11 publicly available airborne lidarderived data sets demonstrate that the proposed method averagely produces better results than two promising interpolation methods including regularized spline with tension (RST) and gridded data-based robust thin plate spline (RTPS). Specifically, the image of RTPS is too smooth to retain terrain details. Although RST can keep subtle terrain features, it is distorted by some misclassified object points (i.e., pseudooutliers). The proposed method obtains a good tradeoff between resisting the effect of outliers and preserving terrain features. Overall, the proposed method can be considered as an alternative for interpolating lidar-derived data sets potentially including outliers. [ABSTRACT FROM AUTHOR]

Data pits commonly appear in lidar-derived canopy height models (CHMs) owing to the penetration ability of airborne light detection and ranging (lidar) into tree crowns. They have a seriously negative effect on the quality of tree detection and subsequent biophysical measurements. In this study, we propose an algorithm based on robust locally weighted regression and robust z-scores for the construction of a pit-free CHM. A significant advantage of the new algorithm is that it is parameter free, which makes it efficient and robust for practical applications. Simulated and airborne lidar-derived data sets are employed to assess the performance of the new method for CHM construction, and its results are compared to those of three classical methods, namely the natural neighbor (NN) interpolation of the highest point method (HPM), mean filter, and median filter. The results from the simulated data set demonstrate that our algorithm is more accurate compared to the three classical methods for generating pit-free CHMs in the presence of data pits. CHM construction using the lidar-derived data set shows that, compared to the classical methods, the new method has a better ability to remove data pits as well as preserving the edges, shapes, and structures of canopy gaps and crowns. Moreover, the proposed method performs better compared to the classical methods in deriving plot-level maximum tree heights from CHMs. Thus, the new method shows high potential for pit-free CHM construction. [ABSTRACT FROM AUTHOR]

High accuracy surface modeling method (HASM) has the capacity of constructing a surface with a higher accuracy than classical methods. However, surface construction with HASM requires a large amount of computing time due to the intensive computing arising from solving the large linear systems. An effective linear system solving method is essential to cutting down the computing time so that HASM could be applied in constructing large scale surface with high resolution. This paper aims to study the performances of seven linear system solution methods, including HASM-GS, HASM-MGS, HASM-DSPM, HASM-CG, HASM-PCG, HASM-GCG and HASM-GPCG, thus providing a guideline for choosing the right method for HASM under different situations of computing environment. These seven methods were tested against a real-world DEM surface in the Boyang lake basin, which owns a complex topography. Three performance metrics were used, including convergence rates, computing time and consumed memory size. We found that to attain the same accuracy HASM-DSPM method consumed the least memory, HASM-GCG and HASM-GPCG needed less computing time than other methods and HASM-PCG had the fastest convergent rate. This study could provide valuable suggestions for HASM users to make a wise choice of linear system solving method for different surface construction applications. [ABSTRACT FROM AUTHOR]

Irrigation water demand accounts for more than 95% of the total water use in the Kaidu-kongqi River Basin. Determination of the spatial and temporal trends in irrigation water demand is important for making sustainable and wise water management strategies in this highly water deficit region. In this study, the spatial and temporal trends in irrigation water demand as well as net crop irrigation water requirements for nine major crops during 1985-2009 were analyzed by combining the Penman-Monteith equation recommended by Food and Agriculture Organization (FAO) and GIS technology. The regional water stress was also evaluated based on the total irrigation water demand and river discharge at the annual and monthly scales. The results indicated that the annual irrigation water demand in this arid region showed a significant increasing trend during the past 25 years. Total irrigation water demand increased from 14.68 × 108 m³ in 1985 to 34.15 × 108 m³ in 2009. The spatial pattern of total irrigation water demand was significantly affected by the changes in cotton growing area. Due to differences in crop planting structure, the monthly average irrigation water demands in Korla City and Yuli County amounted to the peak in July, while those in other regions reached the maximum in June. Although the annual river runoff was much larger than the irrigation water demand, there was serious water deficit during the critical water use period in May and June in some dry years. The presented study provides important information for managers and planners on sustainable use of water resources in this arid region. [ABSTRACT FROM AUTHOR]

Light detection and ranging (LiDAR) technique is currently one of the most important tools for collecting elevation points with a high density in the context of digital elevation model (DEME) construction. However, the high density data always leads to serious time and memory consumption problems in data processing. In this paper, we have developed a thin plate spline (TPS)-based feature-preserving (TPS-F) method for LiDAR-derived ground data reduction by selecting a certain amount of significant terrain points and by extracting geomorphological features from the raw dataset to maintain the accuracy of constructed DEMs as high as possible, while maximally keeping terrain features. We employed four study sites with different topographies (i.e., flat, undulating, hilly and mountainous terrains) to analyze the performance of TPS-F for LiDAR data reduction in the context of DEME construction. These results were compared with those of the TPS-based algorithm without features (TPS-W) and two classical data selection methods including maximum z-tolerance (Max-Z) and the random method. Results show that irrespective of terrain characteristic, the two versions of TPS-based approaches (i.e., TPS-F and TPS-W) are always more accurate than the classical methods in terms of error range and root means square error. Moreover, in terms of streamline matching rate (SMR), TPS-F has a better ability of preserving geomorphological features, especially for the mountainous terrain. For example, the average SMR of TPS-F is 89.2% in the mountainous area, while those of TPS-W, max-Z and the random method are 56.6%, 34.7% and 35.3%, respectively. [ABSTRACT FROM AUTHOR]

Remote-sensing-derived elevation data sets often suffer from noise and outliers due to various reasons, such as the physical limitations of sensors, multiple reflectance, occlusions and low contrast of texture. Outliers generally have a seriously negative effect on DEM construction. Some interpolation methods like ordinary kriging (OK) are capable of smoothing noise inherent in sample points, but are sensitive to outliers. In this paper, a robust algorithm of multiquadric method (MQ) based on an Improved Huber loss function (MQ-IH) has been developed to decrease the impact of outliers on DEM construction. Theoretically, the improved Huber loss function is null for outliers, quadratic for small errors, and linear for others. Simulated data sets drawn from a mathematical surface with different error distributions were employed to analyze the robustness of MQ-IH. Results indicate that MQ-IH obtains a good balance between efficiency and robustness. Namely, the performance of MQ-IH is comparative to those of the classical MQ and MQ based on the Classical Huber loss function (MQ-CH) when sample points follow a normal distribution, and the former outperforms the latter two when sample points are subject to outliers. For example, for the Cauchy error distribution with the location parameter of 0 and scale parameter of 1, the root mean square errors (RMSEs) of MQ-CH and the classical MQ are 0.3916 and 1.4591, respectively, whereas that of MQ-IH is 0.3698. The performance of MQ-IH is further evaluated by qualitative and quantitative analysis through a real-world example of DEM construction with the stereo-images-derived elevation points. Results demonstrate that compared with the classical interpolation methods, including natural neighbor (NN), OK and ANUDEM (a program that calculates regular grid digital elevation models (DEMs) with sensible shape and drainage structure from arbitrarily large topographic data sets), and two versions of MQ, including the classical MQ and MQ-CH, MQ-IH has a better ability of maximally reducing the impact of outliers, while faithfully preserving terrain features. Theoretically, MQ-IH is not a promising interpolation method, and some side effects can be found from its simulation results. For example, the contours of MQ-IH are coarser than ANUDEM in some locations of the real-world study site, and its hill shade may not strictly agree with the real-world surface at rough terrain. Furthermore, the computing cost of MQ-IH is much bigger than that of the classical interpolation methods. However, compared with the classical interpolation methods, MQ-IH has significant potential for interpolating remote-sensing-derived elevation data sets. [ABSTRACT FROM AUTHOR]