确定似大地水准面的Hotine-Helmert边值解算模型

空间大地测量技术的发展使大地高的观测成为可能,从而为第二大地边值问题的研究带来了新的机遇,本文对基于Helmert第二压缩法的第二边值问题(简称为Hotine-Helmert边值问题)展开研究。首先介绍了地形直接、间接影响的定义与算法,然后推导了Hotine-Helmert边值问题的解算模型。Hotine-Helmert边值理论无须计算地形压缩对重力的次要间接影响,因而较Stokes-Helmert边值理论更简单。此外,文中引入了一种低阶修正的Hotine截断核函数,该核函数较传统的截断核函数能有效地改善似大地水准面的解算精度。为了验证本文构建的Hotine-Helmert边值解算模型的有效性和实用性,本文将EIGEN-6C4模型的前360阶作为参考模型,利用Hotine-Helmert边值解算模型构建了我国中部地区6°×4°范围、1.5′×1.5′分辨率的重力似大地水准面,其精度达到±4.8 cm。     

Meissel-Stokes核函数应用于区域大地水准面分析

为提高区域大地水准面计算精度,基于EGM2008地球重力场位系数模型分析Meissel-Stokes核函数、截断误差系数以及截断误差。选取实验区,采用移去-恢复法评价Meissel-Stokes核函数计算大地水准面的精度。结果表明:Meissel-Stokes核函数及其截断误差系数收敛速度快;截断误差小且稳定。在积分半径不易扩展的情况下,应用Meissel-Stokes核函数计算区域大地水准面,比标准Stokes计算大地水准面精度略高。     

重力场模型谱组合精化区域似大地水准面

针对区域似大地水准面模型构建中地面重力数据分布不均匀且较难获取的难题，该文提出一种基于谱组合重力场模型的似大地水准面精化方法，实现了卫星重力场模型与高阶重力场模型的优化组合，结合293个实测全球卫星导航系统(GNSS)水准数据，构建了山东省范围的高精度似大地水准面模型。利用GNSS水准数据检核表明：在山东省范围，利用谱组合重力场模型构建的区域似大地水准面模型具有较高的精度，能够满足大多数工程要求，与采用地面重力数据构建的似大地水准面模型对比发现，在不使用地面重力数据的情况下也可构建厘米级精度似大地水准面模型，具有成本低、效率高的明显优势，便于及时开展区域似大地水准面的维护更新。     

北京市名山高程测量的方法与实践

介绍北京市松山、云蒙山和上方山海拔高程的确定方法,包括GPS数据解算、正常高确定和计算似大地水准面与大地水准面之间的差异,并根据实测数据对最终成果进行精度评定.     

区域重力大地水准面确定的相对精度估计

以频域解析方法,研究由地面重力数据、全球住模型确定区域重力大地水准面的相对精度估计.首先由Stokes公式的数值积分推导地面重力数据与球谐系敬的精度关系;再由"移去-恢复"方法的空域截断逼近模式和协方差函数的球谐表达,分别推导内区地面重力数据之误差、外区全球位模型之误差与区域重力大地水准面之相对精度的解析关系;为便于计算,提出将内区地面重力数据和外区全球位模型的频域截断误差合并,再按频段重新划分为两部分:①全球范围--地面重力数据对应频率以上的截断;②外区范围--介于全球位模型最高频率与地面重力数据对应频率之间的截断,以经验阶方差模型分别估计之.模拟计算显示了地面重力数据之精度、分辨率、积分半径和全球位模型之精度、分辨率与区域重力大地水准面之相时精度的具体对应关系.本文研究同样适用于区域重力似大地水准面的确定.     

矿区似大地水准面拟合模型适用性分析

为提高GPS/水准数据拟合高程异常的逼近精度,利用某矿区控制网观测的GPS/水准数据,计算控制点高程异常值,分别用加权平均法、Shepard插值、多项式曲面拟合和多面函数法4种方法对矿区似大地水准面进行拟合。对各种拟合模型的原理和拟合精度进行比较分析。在此基础上,根据似大地水准面物理和几何特性,采用组合拟合的方法进一步改善似大地水准面的拟合精度,利用该模型内插高程异常值与实测GPS/水准点的高程异常值比较,其均方根误差不超过0.02 m,结合GPS观测成果可以取代四等及以下几何水准测量。     

格尔木市高精度数字高程基准模型研究

大地水准面（数字高程基准）为国家高程基准的建立与维持提供了全新的思路。然而,受限于地形、重力数据等原因,高原地区高精度数字高程基准模型的建立一直是大地测量领域的难题。本文以格尔木地区为例,探讨了高原地区高精度数字高程基准模型的建立方法。首先,基于重力和地形数据,由第二类Helmert凝集法计算了格尔木重力似大地水准面。在计算中,考虑到高原地形对大地水准面模型的影响,采用了7.5″×7.5″分辨率和高精度的地形数据来恢复大地水准面短波部分的方法,以提高似大地水准面的精度。然后,利用球冠谐调和分析方法将GNSS水准与重力似大地水准面联合,建立了格尔木高精度数字高程基准模型。与实测的67个高精度GNSS水准资料比较,重力似大地水准面的外符合精度为3.0 cm,数字高程基准模型的内符合精度为2.0 cm。     

陆海交界区域厘米级精度似大地水准面的确定

为了得到我国某陆海交界区厘米级精度的区域（似）大地水准面,利用43个高精度GPS／水准点和1045个实测重力点数据对EGM96,WDM94和GFZ计算的局部重力（似）大地水准面进行了比较与评价。结果表明,在该测区用移去．恢复法确定重力（似）大地水准面时,EGM96应该是首选参考重力场模型。该测区处在陆海交界处,海域无GPS／水准数据。经比较发现,采用距离倒数加权平均法将该区重力似大地水准面拟合于GPS／水准数据比在大范围使用的多项式法效果更好。采用该方法计算的测区（似）大地水准面精度优于3cm。     

陆海交界区域厘米级精度似大地水准面的确定

为了得到我国某陆海交界区厘米级精度的区域(似)大地水准面,利用43个高精度GPS/水准点和1 045个实测重力点数据对EGM96,WDM94和GFZ计算的局部重力(似)大地水准面进行了比较与评价。结果表明,在该测区用移去-恢复法确定重力(似)大地水准面时,EGM96应该是首选参考重力场模型。该测区处在陆海交界处,海域无GPS/水准数据。经比较发现,采用距离倒数加权平均法将该区重力似大地水准面拟合于GPS/水准数据比在大范围使用的多项式法效果更好。采用该方法计算的测区(似)大地水准面精度优于3cm。     

几何法确定区域似大地水准面及工程应用

为了充分利用GPS观测数据中的高程信息,获得GPS观测点的正常高,本文利用EGM96全球重力场模型和某区域GPS/水准数据,采用移去-恢复技术以几何方法计算了该区域2.5'×2.5'分辨率似大地水准面模型,经过内、外精度的检验,似大地水准面模型的精度优于0.07 m。结合区域内某工程布设的E级GPS控制网数据进行正常高的计算,并和水准实测高程进行比较,结果表明几何法确定的小区域似大地水准面模型结合GPS观测信息可以代替低等级的水准测量,满足一般工程对高程测量的需求。     

The northern European geoid: a case study on long-wavelength geoid errors

 The long-wavelength geoid errors on large-scale geoid solutions, and the use of modified kernels to mitigate these effects, are studied. The geoid around the Nordic area, from Greenland to the Ural mountains, is considered. The effect of including additional gravity data around the Nordic/Baltic land area, originating from both marine, satellite and ground-based measurements, is studied. It is found that additional data appear to increase the noise level in computations, indicating the presence of systematic errors. Therefore, the Wong–Gore modification to the Stokes kernel is applied. This method of removing lower-order terms in the Stokes kernel appears to improve the geoid. The best fit to the global positioning system (GPS) leveling points is obtained with a degree of modification of approximately 30. In addition to the study of modification errors, the results of different methods of combining satellite altimetry gravity and other gravimetry are presented. They all gave comparable results, at the 6-cm level, when evaluated for the Nordic GPS networks. One dimensional (1-D) and 2-D fast Fourier transform (FFT) methods are also compared. It is shown that even though methods differ by up to 6 cm, the fit to the GPS is essentially the same. A surprising conclusion is that the addition of more data does not always produce a better geoid, illustrating the danger of systematic errors in data. Received: 4 July 2001 / Accepted: 21 February 2002     

重力异常向上延拓严密改化模型及向下延拓应用

重力异常向上延拓全球积分模型在航空重力测量数据质量评估和向下延拓迭代计算等领域具有广泛的应用。为了消除积分核函数奇异性影响,需要对该模型进行基于积分恒等式的移去-恢复转换及全球积分域的分区改化处理。在此过程中,传统改化处理方法往往忽略了全球积分过渡到局域积分引起的积分恒等式偏差影响,从而导致不必要的计算模型误差,最终影响向上延拓计算结果的可靠性,甚至影响向下延拓迭代解算结果的稳定性。针对此问题,本文开展了重力异常向上延拓积分模型改化及向下延拓应用分析研究,依据实测数据保障条件和积分恒等式适用条件要求,导出了重力异常向上延拓积分模型的分步改化公式,提出了补偿传统改化模型缺陷的修正公式,并将最终的严密改化模型应用于重力异常向下延拓迭代解算。使用超高阶地球位模型EGM2008作为标准位场开展数值计算检验,分别对重力异常向上延拓分步改化模型的计算精度及在向下延拓迭代解算中的应用效果进行了检核评估,验证了采用严密改化模型的必要性和有效性。     

Comparison of two forest types using adeos 1-polder data over west coast of India

Parametric kernel-based models for Bi-directional Reflectance Distribution Function (BRDF) of vegetated land surfaces attempt to describe the BRDF as a linear superposition of a set of kernels describing basic angular shapes. The kernels and associated model coefficients are derived as approximations to the underlying theories of radiative transfer as well as geometric-optical scattering within vegetation canopies. In the present study, five kernels viz. Ross-Thick, Li-Sparse, Ross-Thin, Li-Dense and Roujean geometric-optical (GO) kernels along with a purely empirical model were examined for their suitability in describing the angular characteristics of two forest types (semi-evergreen and deciduous) during November 1996-March 1997 using ADEOS-1 POLDER data. The statistically significant kernels and their combinations with respect to red and NIR (near infra-red) wavelengths were determined using correlation analysis. Atmospherically corrected POLDER measured reflectances at red and NIR wavelengths were regressed with different kernels and their combinations. It is found that the best fitting volume and GO kernel combinations depend on forest type, wavelength and date.     

Improved convergence rates for the truncation error in gravimetric geoid determination

 When Stokes's integral is used over a spherical cap to compute a gravimetric estimate of the geoid, a truncation error results due to the neglect of gravity data over the remainder of the Earth. Associated with the truncation error is an error kernel defined over these two complementary regions. An important observation is that the rate of decay of the coefficients of the series expansion for the truncation error in terms of Legendre polynomials is determined by the smoothness properties of the error kernel. Previously published deterministic modifications of Stokes's integration kernel involve either a discontinuity in the error kernel or its first derivative at the spherical cap radius. These kernels are generalised and extended by constructing error kernels whose derivatives at the spherical cap radius are continuous up to an arbitrary order. This construction is achieved by smoothly continuing the error kernel function into the spherical cap using a suitable degree polynomial. Accordingly, an improved rate of convergence of the spectral series representation of the truncation error is obtained. Received: 21 April 1998 / Accepted: 4 October 1999     

Wavelet evaluation of the Stokes and Vening Meinesz integrals

The wavelet transform is a powerful tool in evaluating some singular geodetic integrals. Due to its localization properties in both the time (space) and frequency (scale) domains, and because the kernels of some geodetic integrals have singular points and decay smoothly and quickly away from the singularities, many wavelet transform coefficients of the kernels become zeros or negligible, and only a small number of wavelet transform coefficients are significant. It is thus possible to significantly compress the kernels of these integrals on a wavelet basis by neglecting the zero coefficients and the small coefficients below a certain threshold. Therefore, wavelets provide a convenient way of efficiently evaluating these integrals in terms of fast computation and savings of computer memory. A modified algorithm for the wavelet evaluation of Stokes' integral is presented. The same modified algorithm is applied to the evaluation of the Vening Meinesz integral, whose kernel has a stronger singularity than does Stokes' kernel. Numerical examples illustrate the efficiency and accuracy of the wavelet methods.Acknowledgments.The author express their sincere thanks to Dr. Salamonwicz for providing his PhD thesis. E-mail correspondence between the authors and Dr. Barthelmes and Dr. Benciolini contributed to the work. R. Benciolini and the other two anonymous reviewers are thanked for their constructive comments. Support for this research was provided by research grants to Dr. Sideris from the Natural Sciences and Engineering Reserch Council of Canada (NSERC) and the Geomatics for Informed Decisions (GEOIDE) Network of Centres of Excellence. The MATLAB Wavelet Toolbox package was used as the platform to develop the software in this project.     

Composite kernels for hyperspectral image classification

This letter presents a framework of composite kernel machines for enhanced classification of hyperspectral images. This novel method exploits the properties of Mercer's kernels to construct a family of composite kernels that easily combine spatial and spectral information. This framework of composite kernels demonstrates: 1) enhanced classification accuracy as compared to traditional approaches that take into account the spectral information only: 2) flexibility to balance between the spatial and spectral information in the classifier; and 3) computational efficiency. In addition, the proposed family of kernel classifiers opens a wide field for future developments in which spatial and spectral information can be easily integrated.     

组合核支持向量回归提取高光谱影像不透水面

由于城市地表组成的复杂性,基于单核函数的支持向量回归模型很难满足精度。本文结合空间-光谱组合核函数和支持向量回归,提出了一种提取高光谱影像不透水面丰度的改进算法。首先从高光谱遥感图像上提取波谱特征和多通道灰度共生矩阵空间纹理特征,选取研究区10%像元特征数据作为训练数据,以线性加权求和核为多核组合方式,建立结合光谱信息和空间信息的组合核支持向量回归模型。然后,用生成的回归模型预测未知像元不透水面丰度值。最后,对实验结果进行评价。在模拟数据试验中,本文算法比单核回归均方根误差平均降低1.4%,决定系数比单核回归平均提高0.6%。在Hyperion数据两组试验中,该算法比单核回归均方根误差平均降低1.8%,决定系数比单核回归平均提高11.7%。模拟和真实两种高光谱数据实验中,本文算法均得到了空间形态上更准确的不透水面结果,单核回归结果存在失真现象。研究结果表明:本文算法能够有效提取城市不透水面丰度,与单核方法相比有较明显的精度提升。     

Vertical and horizontal spheroidal boundary-value problems

Vertical and horizontal spheroidal boundary-value problems (BVPs), i.e., determination of the external gravitational potential from the components of the gravitational gradient on the spheroid, are discussed in this article. The gravitational gradient is decomposed into the series of the vertical and horizontal vector spheroidal harmonics, before being orthogonalized in a weighted sense by two different approaches. The vertical and horizontal spheroidal BVPs are then formulated and solved in the spectral and spatial domains. Both orthogonalization methods provide the same analytical solutions for the vertical spheroidal BVP, and give distinct, but equivalent, analytical solutions for the horizontal spheroidal BVP. A closed-loop simulation is performed to test the correctness of the analytical solutions, and we investigate analytical properties of the sub-integral kernels. The systematic treatment of the spheroidal BVPs and the resulting mathematical equations extend the theoretical apparatus of geodesy and of the potential theory.     

On the accurate numerical evaluation of geodetic convolution integrals

In the numerical evaluation of geodetic convolution integrals, whether by quadrature or discrete/fast Fourier transform (D/FFT) techniques, the integration kernel is sometimes computed at the centre of the discretised grid cells. For singular kernels—a common case in physical geodesy—this approximation produces significant errors near the computation point, where the kernel changes rapidly across the cell. Rigorously, mean kernels across each whole cell are required. We present one numerical and one analytical method capable of providing estimates of mean kernels for convolution integrals. The numerical method is based on Gauss-Legendre quadrature (GLQ) as efficient integration technique. The analytical approach is based on kernel weighting factors, computed in planar approximation close to the computation point, and used to convert non-planar kernels from point to mean representation. A numerical study exemplifies the benefits of using mean kernels in Stokes’s integral. The method is validated using closed-loop tests based on the EGM2008 global gravity model, revealing that using mean kernels instead of point kernels reduces numerical integration errors by a factor of ~5 (at a grid-resolution of 10 arc min). Analytical mean kernel solutions are then derived for 14 other commonly used geodetic convolution integrals: Hotine, Eötvös, Green-Molodensky, tidal displacement, ocean tide loading, deflection-geoid, Vening-Meinesz, inverse Vening-Meinesz, inverse Stokes, inverse Hotine, terrain correction, primary indirect effect, Molodensky’s G1 term and the Poisson integral. We recommend that mean kernels be used to accurately evaluate geodetic convolution integrals, and the two methods presented here are effective and easy to implement.