规则格网上双方向性插值DEM精度研究

简要介绍双方向性插值方法的原理,并以线性插值和三次插值过程中的误差传播为基础,利用误差传播定律,推导出双方向性插值方法生成的DEM曲面平均高程误差值,并对线性插值和三次插值公式本身的精度进行分析。     

双线性插值近似网格的栅格数据投影变换

在栅格数据投影变换常用算法的基础上,介绍一种新的基于双线性内插近似网格算法的栅格数据投影变换,并详细描述其算法过程,认真分析其特点,并证明该算法在大范围栅格数据投影变换中有显著优点。     

凸包插值算法的改进与实现

研究了凸包插值算法在n维欧拉空间构造Delaunay三角网的应用,提出了其三角剖分不惟一的重大缺陷,采用了限制区域生长法解决了这一缺陷。     

Triangulated spherical splines for geopotential reconstruction

We present an alternate mathematical technique than contemporary spherical harmonics to approximate the geopotential based on triangulated spherical spline functions, which are smooth piecewise spherical harmonic polynomials over spherical triangulations. The new method is capable of multi-spatial resolution modeling and could thus enhance spatial resolutions for regional gravity field inversion using data from space gravimetry missions such as CHAMP, GRACE or GOCE. First, we propose to use the minimal energy spherical spline interpolation to find a good approximation of the geopotential at the orbital altitude of the satellite. Then we explain how to solve Laplace’s equation on the Earth’s exterior to compute a spherical spline to approximate the geopotential at the Earth’s surface. We propose a domain decomposition technique, which can compute an approximation of the minimal energy spherical spline interpolation on the orbital altitude and a multiple star technique to compute the spherical spline approximation by the collocation method. We prove that the spherical spline constructed by means of the domain decomposition technique converges to the minimal energy spline interpolation. We also prove that the modeled spline geopotential is continuous from the satellite altitude down to the Earth’s surface. We have implemented the two computational algorithms and applied them in a numerical experiment using simulated CHAMP geopotential observations computed at satellite altitude (450 km) assuming EGM96 (n _max = 90) is the truth model. We then validate our approach by comparing the computed geopotential values using the resulting spherical spline model down to the Earth’s surface, with the truth EGM96 values over several study regions. Our numerical evidence demonstrates that the algorithms produce a viable alternative of regional gravity field solution potentially exploiting the full accuracy of data from space gravimetry missions. The major advantage of our method is that it allows us to compute the geopotential over the regions of interest as well as enhancing the spatial resolution commensurable with the characteristics of satellite coverage, which could not be done using a global spherical harmonic representation. The results in this paper are based on the research supported by the National Science Foundation under the grant no. 0327577.     

Non-stationary covariance function modelling in 2D least-squares collocation

Standard least-squares collocation (LSC) assumes 2D stationarity and 3D isotropy, and relies on a covariance function to account for spatial dependence in the observed data. However, the assumption that the spatial dependence is constant throughout the region of interest may sometimes be violated. Assuming a stationary covariance structure can result in over-smoothing of, e.g., the gravity field in mountains and under-smoothing in great plains. We introduce the kernel convolution method from spatial statistics for non-stationary covariance structures, and demonstrate its advantage for dealing with non-stationarity in geodetic data. We then compared stationary and non- stationary covariance functions in 2D LSC to the empirical example of gravity anomaly interpolation near the Darling Fault, Western Australia, where the field is anisotropic and non-stationary. The results with non-stationary covariance functions are better than standard LSC in terms of formal errors and cross-validation against data not used in the interpolation, demonstrating that the use of non-stationary covariance functions can improve upon standard (stationary) LSC.     

海洋磁力测量中测线方向的选择方法

采用测线间磁异常插值精度分析了测线方向对磁测的影响;探讨了测线方向的选择方法,结合我国不同海区磁异常分布特征给出了不同海区测量时测线方向的选择详细方案。结论表明:海洋磁力测量中,测线方向对磁异常插值精度的影响不大,但是为了获得更高的磁异常插值精度,测线方向应尽量沿磁异常梯度方向布设。给出的测线方向的选择方法和我国不同海区测线方向选择详细方案可供作业人员参考。     

海底地形数据处理与DEM生成

分析了海底DEM生成经常会遇见的问题,介绍了坐标变换、接边处理、插值加密等数据预处理方法,最终生成海底规则格网DEM。     

普通Kriging插值算法研究

地统计学是进行空间变异分析的有力工具.文章说明了地统计学的基本原理及其实现算法,并应用该算法对研究区域的草句白浆土有效磷的空间变异特征进行了分析,通过计算得出的地统计学参数说明了土壤有效磷具有中等空间相关性,并生成了该养分的空间分布图.     

日本周边海域温跃层特性研究及温度剖面模型构建

使用"中国ARGO实时资料中心"提供的2005—2009年连续5 a日本南部海域的ARGO浮标资料,通过建立标准层,统计并提取各标准层的剖面温度,采用AKIMA插值方法计算各层的温度梯度。较详细分析了该海域温度垂直剖面特征及温跃层的分布规律。在此基础上,使用分段三次多项式拟合方法建立了适用于该海区的温度垂直剖面模型,并利用2010年ARGO数据对模型进行了验证,研究表明,经过模型模拟的剖面温度与实测温度之间的相关系数可达到0.99以上。     

三次样条插值在浮性频率计算中的应用

在普遍利用海洋观测站标准层的温度、盐度数据估算稳定层化海水浮性频率的传统算法中引入自然三次样条插植方法，不仅可使计算结果的精度得以提高，而且还可方便地对观测深度范围内的任一水层进行插值，从而获得总体曲率最小的温度、盐度和浮性频率随深度的分布。本文以南海（13°30’N，118°30’E；9°30’N，lll°30’E）及北太平洋（40°00’N130°00’W）3个测站为例，说明了在浮性频率的传统算法中引入自然三次样条插值法的可行性。