GPS天线相位中心校正对GRACE卫星精密定轨的影响

使用GRACE卫星星载GPS观测数据,研究地面获取的先验PCV模型和利用残差法估计的在轨PCV模型对低轨卫星精密定轨的影响,并采用GFZ精密轨道对比和SLR检核手段对其进行评估。结果表明,使用先验PCV模型会降低GRACE卫星定轨精度,相反利用在轨PCV模型可以提高定轨精度,提升数量级可达mm级。     

����Kalman�˲���GEO�����˶�ѧ����

?????GEO????????????????????????????????????Kalman???????GEO?????????????????????????????????????η????????????????????????????????????????????????????????????????????????????????????з?????????????????????????????????????????Kalman???????????н????????????????к???????????????????????,????????????????Kalman????????????????????????????,???????Kalman???????????????????????????????????????????Kalman?????????????????????GEO?????????????     

�й���½�ٶȳ�ģ�ͽ��������о�

???й?????????????????????????????????5?????壬????????????????С???????÷?????С?????????????????????3????????????????о??й?????????????????????????????????????????????????????С???????÷????????壨??????????????????壩????????????????????|?????????????????????????????????????????????????С???????÷?????й????1°×1°?????????????????????????????????С??     

具有先验信息的地面三维激光扫描地形测量数据去噪算法

为防止三维激光扫描数据中的噪声对实景重构造成扭曲,基于消噪的相关理论,提出了3种具有先验信息的去噪算法,利用高程注记点所提供的先验信息作为去噪依据,分别对比高程插值、最小二乘拟合平面及趋势面,选取适当阈值消去噪声。利用Timble GX200地面三维激光扫描系统获取的地形测量数据进行实验分析,结果表明:高程注记点可以作为先验信息为点云去噪提供依据,并确保去噪后基本符合实体特征。     

基于GIS的“一带一路”地区气温插值方法比较研究

“一带一路”倡议是新时期中国为加强对外开放提出的全球化合作倡议,资源环境的优化配置对全球化发展意义重大。气温作为重要的基础数据和输入要素,对其进行空间化处理是实现大尺度区域资源环境优化配置的前提。本文基于地理信息技术（GIS）,运用距离平方反比法（IDS）、协同克里格法（CK）、回归距离平方反比法（RIDS）和回归协同克里格法（RCK）,对“一带一路”地区1980—2017年的2679个气象站点的月平均气温和年平均气温数据进行插值,获得了“一带一路”地区10 km分辨率的气温空间分布数据。交叉验证结果表明：① IDS、CK、RIDS和RCK插值法在整体上均较好地展示了“一带一路”地区气温的地理空间分布规律,4种插值方法的月均气温的均方根误差分别在1.93~2.43、1.78~2.14、1.31~2.23和1.23~1.92 ℃之间;年均气温的均方根误差分别为1.94、1.83、1.37和1.27 ℃;② 在“一带一路”地区,加入协变量分析的CK插值精度整体优于IDS,并且削弱了IDS的极值现象;③ RIDS和RCK对年均气温的插值精度分别较IDS和CK提高了29.4%和30.6%,表明加入地理要素并进行残差修正的插值精度得到了进一步提高。总体来看,RCK插值法对气温数据的插值精度最高,可以考虑将此方法作为“一带一路”地区温度等气象要素的插值方法。     

多波束测深中声速剖面的横向加密方法

针对声速剖面在横向空间上的低密度、分布不均匀导致水深测量误差的问题,提出了一种容易实现且有效的思路,采用分层取样、空间插值的方法,在后处理过程中人为地增加声速剖面的数量,使其密度足够大、分布更均匀,达到改善声速改正的效果和提高多波束测深精度的目的。     

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.     

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

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

局部最优估计的参数求解方法

结合了飞行器轨道拟合的实际问题,介绍了在约束条件下基于局部最优估计的参数求解,阐述了其经验和方法。