基于高精度曲面模型的DEM构建与误差分析

引入地形表达误差(terrain representation error,Etr),选择标准曲面和甘肃省董志塬地区作为研究对象,利用窗口分析法实现Etr的提取;用统计分析法得出Etr随网格分辨率变化的回归方程;根据误差传播定律计算DEM中误差。数值结果表明,该方法能更准确的计算HASM生成的DEM精度;相同的采样数下,HASM较传统方法(IDW,Spline和Kriging)能生成更高精度和分辨率的DEM。在难以获取已知数据的地区,HASM提供了生成相对准确DEM的高效工具。     

高精度曲面模型的解算

为了提高Gass-Seidel(GS)算法的收敛速度,提出了改进的GS算法(MGS),用于解算高精度曲面模型(HASM)(HASM-MGS)。以高斯合成曲面为研究对象,将HASM-MGS与HASM-GS和Matlab提供的函数进行对比,结果表明,达到相同的模拟中误差,HASM-MGS计算时间远小于HASM-GS和Matlab提供的函数;HASM-MGS计算时间与模拟区域的网格数呈非常好的线性关系,时间复杂度比传统的方法降低两个数量级。     

A high accuracy surface modeling method based on GPU accelerated multi‐grid method

A high accuracy surface modeling method (HASM) has been developed to provide a solution to many surface modeling problems such as DEM construction, surface estimation and spatial prediction. Although HASM is able to model surfaces with a higher accuracy, its low computing speed limits its popularity in constructing large scale surfaces. Hence, the research described in this article aims to improve the computing efficiency of HASM with a graphic processor unit (GPU) accelerated multi‐grid method (HASM‐GMG). HASM‐GMG was tested with two types of surfaces: a Gauss synthetic surface and a real‐world example. Results indicate that HASM‐GMG can gain significant speedups compared with CPU‐based HASM without acceleration on GPU. Moreover, both the accuracy and speed of HASM‐GMG are superior to the classical interpolation methods including Kriging, Spline and IDW.     

An application of Coons patch to generate grid-based digital elevation models

We have presented a new method for generating grid-based digital elevation models (DEMs) based on Coons patch. A series of Gaussian synthetic surfaces with various surface complexities were employed to comparatively analyze the simulation accuracies of Coons patch and the classical interpolation methods, such as IDW, Natural neighbor, SPLINE and KRIGING. Results indicate that irrespective of terrain complexity, Coons patch is much more accurate than the classical interpolation methods, and the accuracy of Coons patch is about 327.30, 15.19, 3.66 and 40.07 times as high as the accuracies of IDW, Natural neighbor, SPLINE and KRIGING on average. A real-world test area was used to test the effect of sampling density on the simulation accuracy. The results indicate that Coons patch has a higher accuracy than the classical interpolation methods regardless of sampling density. In a word, Coons patch as an accurate interpolation method can be used for grid-based DEM construction.     

A Multi‐Grid Method of High Accuracy Surface Modeling and Its Validation

A method of high accuracy surface modeling (HASM) has been constructed to find a solution for error problems that had long troubled surface modeling in geographical information systems (GIS). It is found that when a preconditioned conjugate gradient (PCG) algorithm is used to solve the large sparse linear system, which HASM can be transferred into, HASM performs best in terms of simulation compared with all other algorithms. But its computing speed is not fast enough for all applications. A multi‐grid method is introduced into HASM to try to shorten its computing time. Both numerical and real‐world tests demonstrate that there is a range of stop error (SE). The multi‐grid method of HASM (HASM‐MG) greatly increases computing speed when SEs are within this range, compared with the PCG algorithm of HASM (HASM‐PCG). HASM‐MG is suitable for applications with a need for less accuracy and a shorter computing time. HASM‐PCG is appropriate for issues needing higher accuracy. HASM‐MG performs better than HASM‐PCG in flat areas, while HASM‐PCG does better in complex terrainm in terms of accuracy and computing time.     

An Adaptive Method of High Accuracy Surface Modeling and Its Application to Simulating Elevation Surfaces

An adaptive method is employed to speed up computation of high accuracy surface modeling (HASM), for which an error indicator and an error estimator are developed. Root mean‐square error (RMSE) is used as the error estimator that is formulated as a function of gully density and grid cell size. The error indicator is developed on the basis of error surfaces for different spatial resolutions, which are interpolated in terms of the absolute errors calculated at sampled points while paying attention to the landform characteristics. The error surfaces indicate the magnitude and distribution of errors in each step of adaptive refinement and make spatial changes to the errors in the simulation process visualized. The adaptive method of high accuracy surface modeling (HASM‐AM) is applied to simulating elevation surface of the Dong‐Zhi tableland with 27.24 million pixels at a spatial resolution of 10 m × 10 m. Test results show that HASM‐AM has greatly speeded up computation by avoiding unnecessary calculations and saving memory. In addition, HASM‐AM improves simulation accuracy.     

A total error-based multiquadric method for surface modeling of digital elevation models

Digital elevation model (DEM) source data are subject to both horizontal and vertical errors owing to improper instrument operation, physical limitations of sensors, and bad weather conditions. These factors may bring a negative effect on some DEM-based applications requiring low levels of positional errors. Although classical smoothing interpolation methods have the ability to handle vertical errors, they are prone to omit horizontal errors. Based on the statistical concept of the total least squares method, a total error-based multiquadric (MQ-T) method is proposed in this paper to reduce the effects of both horizontal and vertical errors in the context of DEM construction. In nature, the classical multiquadric (MQ) method is a vertical error regression procedure, whereas MQ-T is an orthogonal error regression model. Two examples, including a numerical test and a real-world example, are employed in a comparative performance analysis of MQ-T for surface modeling of DEMs. The numerical test indicates that MQ-T performs better than the classical MQ in terms of root mean square error. The real-world example of DEM construction with sample points derived from a total station instrument demonstrates that regardless of the sample interval and DEM resolution, MQ-T is more accurate than classical interpolation methods including inverse distance weighting, ordinary kriging, and Australian National University DEM. Therefore, MQ-T can be considered as an alternative interpolator for surface modeling with sample points subject to both horizontal and vertical errors.     

1980西安坐标系向2000国家坐标系转换省级兴趣点数据成果分析

首先介绍了1980西安高斯平面坐标系向2000国家高斯平面坐标系转换的两种常用方法:改正量插值法、二维七参数法。其次,在广西测区内,以1∶10000比例尺图幅为单元,选取16个(2组,每组8个)图幅为试验区域,选取3×3块兴趣点数据,使用两种方法进行坐标转换,以从自然资源部获取的2000坐标为真值,分析两种转换成果。结果表明,两种成果均能满足1∶10000平面精度要求,改正量插值法较二维七参数法精度更高,且各图幅精度较均匀。改正量插值法最大点位误差为0.057 m,二维七参数法最大点位误差为0.679 m,出现在北海市。     

基于广义延拓逼近法的QZSS卫星钟差内插精度分析

卫星钟差数据插值是高精度定位数据处理中的重要环节,其插值结果直接影响定位精度,但常用的插值或拟合方法具有不同缺点. 本文尝试将广义延拓逼近法应用于准天顶卫星系统(QZSS)卫星钟差数据的处理中,介绍了Lagrange插值法、切比雪夫拟合法和广义延拓逼近法的原理,以及滑动式与非滑动式的区别;然后使用QZSS钟差数据探讨三种方法的参数(组)取值与插值结果精度的关系;最后比较三种方法在各自最优参数(组)取值情况下对QZSS卫星钟差的插值精度. 结果表明:选取合理的参数组合,广义延拓逼近法完全适用于QZSS卫星钟差的插值,且插值精度明显高于其他两种方法.      

基于Coons曲面的规则格网DEM表面模型

内插是数字高程模型的核心问题。目前的内插模型主要是由离散的格网数据构建的连续曲面,直接以点推面,可能存在较大的地形误差。本文建立的Coons曲面DEM表面模型,首先利用离散的格网数据构造与格网边界相对应的地形剖面曲线的拟合曲线,再基于拟合曲线构建DEM表面模型。实验表明：Coons曲面DEM表面模型是一种高精度的DEM表面模型,其地形模拟误差比直接基于格网数据建立的双线性内插、样条函数内插和移动曲面拟合法的误差都小,实际地形模拟误差与双线性模型相比减少15％-28％,且精度随着构建边界拟合曲线所用格网点的增多而逐渐提高。     

一种综合考虑采样点水平和高程误差的DEM建模算法

为了降低采样点水平和高程误差对数字高程模型（digital elevation model,DEM）建模精度的影响,受总体最小二乘算法启发,以较高精度的多面函数（multiquadric function,MQ）为基函数,发展了整体最小二乘MQ算法（MQ-T）,并分别借助数值实验和实例分析验证模型计算精度。数值实验中,以高斯合成曲面为研究对象,设计了受不同误差分量影响的采样数据,借助MQ-T曲面建模,并将计算结果与传统MQ进行比较。结果表明,当采样点仅受高程误差分量影响时,MQ-T计算结果精度与MQ相当;当采样数据受水平误差分量影响时,MQ-T计算结果中误差小于MQ中误差。实例分析中,以全站仪获取的采样数据为研究对象,借助MQ-T构建测区DEM,并将计算结果与传统插值算法进行比较,如反距离加权（inverse distance weighted,IDW）法、克里金（Kriging）法和澳大利亚国立大学DEM专用插值软件（（Australian National University DEM,ANUDEM）法。精度分析表明,随着采样点密度降低,各种插值算法精度逐步降低;不管采样密度多少,MQ-T计算精度始终高于传统插值算法;对山体阴影图分析表明,MQ-T相比Kriging法有一定峰值削平现象。     

滑动式广义延拓插值法在GLONASS钟差插值中的应用

在对卫星钟差数据进行插值处理时,所采用插值算法的精度,直接影响到卫星钟差插值结果的精度,继而影响了卫星导航定位的精度. 因此在对卫星钟差数据进行插值时,应选择适宜的插值方法. 将Lagrange插值法和切比雪夫多项式拟合法进行滑动,利用这两种传统的插值方法和滑动式广义延拓插值法,分别对历元间隔为5 min的GLONASS卫星钟差数据插值到历元间隔为30 s的钟差数据,再与历元间隔为30 s的精密钟差数据进行对比,分析三种插值方法在GLONASS卫星钟差数据中的应用效果. 结果表明:利用这三种插值方法对GLONASS钟差数据进行插值时,插值精度均能满足要求,且滑动式广义延拓插值法的插值精度最高.      

基于映射函数网络RTK的大气误差内插估计模型

在比较分析已有内插模型的基础上,以卫星高度角、参考站高程、距离为变量,以大气映射函数为基本模型,提出了一种新的适用于高程差异大的大气误差内插模型。利用江苏连续运行参考站系统部分参考站数据,采用不同的内插模型进行处理、比较与分析。结果表明,对于高程差异大、卫星高度角小的用户站,映射函数估计方法、低阶曲面模型能够准确地估计出其与主参考站之间双差对流层误差,其最大均方差不超过2cm;映射函数估计双差电离层误差方法与低阶曲面模型、线性内插模型的估计精度相当。在以上4种方法的比较中,距离线性内插模型的估计精度最差,将近5cm。     

复杂地形景区三维导览图的设计与实现

为满足复杂地形景区对三维地图导览的需求,构建更为清晰的三维地图模型,提供良好的三维导览地图设计方案显得尤为重要。为了更好地与卫星影像进行贴合,DEM数据的采样密度要与卫星影像的分辨率一致。本文通过对比当下常用的四种空间插值方法的适用范围与运用特点,设置相关阈值及权重,直观比较了四种插值方法产生的插值结果,选择更适合复杂景区DEM插值的插值方法并对其进行精度提升,以此得到符合精度要求的DEM数据。最后利用相应地区的卫星影像进行地图投影及影像贴图,两种数据结合构建复杂景区的三维地图模型,给游客带来更为精确直观的定位信息和空间要素信息。     

网络RTK对流层延迟内插模型精度分析

流动站与参考站间双差对流层延迟的精确改正是网络RTK技术的关键。本文采用河南省地基增强系统参考站网7个参考站的观测数据进行试验,选择LIM、LCM、LSM和KRG模型为研究对象,分析了各模型用于不同高差水平流动站对流层延迟改正的效果。试验结果表明：4种常用内插模型中,LIM和KRG模型对流动站对流层延迟改正的精度较优。流动站与主参考站高差达到400 m时,对于低高度角卫星,模型内插精度降低到分米级,不满足用户流动站高精度定位需求。     

土壤重金属的空间插值方法比较分析——以阳谷县为例

空间插值通过采集少量的数据点,利用其中的空间关联,推求该区域内其他位置的属性值。本文以山东省阳谷县土壤重金属Cu采样数据为例进行空间探索性分析,分别采用了反距离权重插值和普通克里金插值两种方法进行空间分布插值模拟。结果表明,针对研究区采样数据,反距离权重插值方法生成的模型平均误差为1.97 mg/kg,总体精度为92%;普通克里金插值模型的平均误差为1.91 mg/kg,总体精度为92.35%,普通克里金插值方法更优。     

小区域GNSS高程异常拟合方法研究

选取恰当的拟合模型是提高高程异常拟合精度的关键．本文利用湖南省新化县城区全球卫星导航系统（GNSS）和水准测量数据,采用反距离加权法、多项式插值法、径向基函数法、克里金插值法等多种曲面拟合方法进行高程异常拟合研究,并进行拟合精度评定．结果表明,径向基函数插值法内插精度最高,其中误差为±0.0158 m;局部多项式插值法外推精度最高,其中误差为±0.0104 m;综合来说,局部多项式插值法在县域尺度高程拟合中精度最高,本文研究结果对小区域GNSS高程拟合方法选取工作具有一定的参考意义.      

单波束测深数据的插值精度分析

基于误差传播定律,对稀疏的离散水深点内插值进行了精度分析,建立了单波束水深内插值中误差表达的数学模型,利用我国南方某海岸带的3个试验区的进行了试验。试验结果表明,反距离加权法、Shepard法和线性插值三角网法,辅以自适应搜索半径法,内插水深值粗差比例普遍低于5%,质量与效率上为较优的插值模型;内插值精度与数据源精度有关,但与数据源密度关系不大;在给定深度测量极限误差情况下,建议在制定有关数字水深模型标准时,对水深在20 m以内的格网点水深值极限误差可设置为0.4 m。     

HASM-AD algorithm based on the sequential least squares

The HASM(high accuracy surface modeling) technique is based on the fundamental theory of surfaces,which has been proved to improve the interpolation accuracy in surface fitting.However,the integral iterative solution in previous studies resulted in high temporal complexity in computation and huge memory usage so that it became difficult to put the technique into application,especially for large-scale datasets.In the study,an innovative model(HASM-AD) is developed according to the sequential least squares on...     

Noisy Data Smoothing in DEM Construction Using Least Squares Support Vector Machines

Since spatial datasets are subject to sampling errors, a smoothing interpolation method should be employed to remove noise during DEM construction. Although least squares support vector machines (LSSVM) have been widely accepted as a classifier, their effect on smoothing noisy data is almost unknown. In this article, the smoothness of LSSVM was explored, and its effect on smoothing noisy data in DEM construction was tested. In order to improve the ability to deal with large datasets, a local method of LSSVM has been developed, where only the neighboring sampling points around the one to be estimated are used for computation. A numerical test indicated that LSSVM is more accurate than the classical smoothing methods including TPS and kriging, and its error surfaces are more evenly distributed. The real‐world example of smoothing noise inherent in lidar‐derived DEMs also showed that LSSVM has a positive smoothing effect, which is approximately as accurate as TPS. In short, LSSVM with a high efficiency can be considered as an alternative smoothing method for smoothing noisy data in DEM construction.