Spatiotemporal filtering of regional GNSS network’s position time series with missing data using principle component analysis

The existing spatiotemporal analysis methods suppose that the involved time series are complete and have the same data interval. However missing data inevitably occur in the position time series of Global Navigation Satellite Systems networks for many reasons. In this paper, we develop a modified principal component analysis to extract the Common Mode Error (CME) from the incomplete position time series. The principle of the proposed method is that a time series can be reproduced from its principle components. The method is equivalent to the method of Dong et al. (J Geophys Res 111:3405–3421, 2006) in case of no missing data in the time series and to the extended ‘stacking’ approach under the assumption of a uniformly spatial response. The new method is first applied to extract the CME from the position time series of the Crustal Movement Observation Network of China (CMONOC) over the period of 1999–2009 where the missing data occur in all stations with the different gaps. The results show that the CMEs are significant in CMONOC. The size of the first principle components for the North, East and Up coordinates are as large as 40, 41 and 37 % of total principle components and their spatial responses are not uniform. The minimum amplitudes of the first eigenvectors are only 41, 15 and 29 % for the North, East and Up coordinate components, respectively. The extracted CMEs of our method are close to the data filling method, and the Root Mean Squared error (RMS) values computed from the differences of maximum CMEs between two methods are only 0.31, 0.52 and 1.55 mm for North, East and Up coordinates, respectively. The RMS of the position time series is greatly reduced after filtering out the CMEs. The accuracies of the reconstructed missing data using the two methods are also comparable. To further comprehensively test the efficiency of our method, the repeated experiments are then carried out by randomly deleting different percentages of data at some stations. The results show that the CMEs can be extracted with high accuracy at the non missing-data epochs. And at the missing-data epochs, the accuracy of extracted CMEs has a strong dependence on the number of stations with missing data.     

基于无几何电离层滤波模型的北斗三号系统相位周跳与中断修复方法

高精度位置服务依赖模糊度正确固定的载波相位观测值,而复杂环境的相位周跳和中断给模糊度固定和数据处理带来了巨大挑战。北斗三号系统向全球播发四频和五频信号,为提升周跳处理性能提供了机遇。本文建立了探测多频超宽巷周跳的伪距-相位无几何组合模型,以及探测窄巷周跳的相位-相位无几何组合模型,并以周跳探测成功概率最大为准则确定了多频最优超宽巷和窄巷组合。由于电离层延迟制约了数据中断和电离层变化剧烈情况下的窄巷周跳探测,进而建立了顾及电离层延迟的窄巷周跳探测滤波模型。试验结果表明,对30 s采样间隔,基于无几何组合模型探测周跳的成功率均优于96%;对于3 min的数据中断,基于无几何组合模型的周跳探测成功率仅为70%,但补偿电离层延迟的滤波模型将成功率提升至95%以上。     

时间相关观测的几种动态数据处理方法

现代大地测量数据中往往存在时间相关性,忽略观测值时间相关性会影响参数估值的精度和可靠性。因此,本文研究了几种时间相关观测的动态数据处理方法。首先,对观测值时间相关性的处理方法进行了扩展和统一。利用极大验后估计原理,提出了第3种思路,并在无历元公用参数和有历元公用参数两类模型下,详细推导了去相关变换法、差分变换法和极大验后估计法这3种思路的动态解,并对它们的特点和等价性进行了论证。其次,为了平衡实际应用中的解算效率且同时有效顾及时间相关性,利用自相关函数表达时间相关性,导出了相应退化形式。最后,通过GPS实测数据验证了本文理论公式的正确性和实用性。     

纳伪概率可控的四舍五入法及其在RTK模糊度固定中的应用

RTK模糊度固定通常采用序贯四舍五入方法,该方法只固定能可靠固定的部分模糊度,因此控制四舍五入的纳伪概率是实现RTK可靠性的关键。从含多个备选假设的假设检验理论出发,研究纳伪概率可控的四舍五入方法,根据风险水平和实数解的质量自适应地确定四舍五入取整区域,从而有效地控制整数固定的纳伪概率。并结合GNSS应用中两种特殊情况,发展两种简单实用的应用策略。采用基线长94.6 km的CORS站GPS双频数据进行网络RTK模糊度固定试验,分析本文纳伪概率可控的四舍五入方法的效果。结果表明,纳伪概率可控的四舍五入方法的固定率与传统四舍五入方法的固定率相当,且较传统方法更能有效地控制纳伪概率,提高整数解的可靠性。     

Optimal Doppler-aided smoothing strategy for GNSS navigation

GPS Solutions - Carrier-phase-smoothed code (CPSC), i.e., smoothing of the code using carrier phases, has widely been used to reduce the code noise in GNSS applications. However, the efficiency of...     

利用PPP反演海潮负荷位移参数

选取高纬地区的3个CORS站以及中纬地区的SHAO测站,利用PPP浮点解计算测站未经海潮负荷位移改正的坐标时间序列,采用傅里叶变换方法提取分潮的周期信息,进而反演其中8个主要分潮的海潮负荷位移参数,将其与5个全球海潮模型DTU10、EOT11a、GOT4.8、TPXO.7.2、FES2004计算得到的海洋负荷位移参数比较。结果表明,PPP反演得到的8个主要分潮的海潮负荷位移参数与海潮模型计算结果具有较好的一致性,两者均方根差异为mm级,说明利用PPP反演海潮负荷位移参数是可行的。     

上海城市二氧化碳排放空间特征

构建上海1 km CO₂排放网格,分析市域（UB1）、市辖区（UB2）、建成区（UB3）和城区（UB4）4个城市范围的CO₂排放特征。上海市域CO2排放空间格局是以中心城区为核心,排放水平向外递减,形成了3个梯度。空间自相关分析表明,排放在空间上存在显著的集聚效应,部分地区高强度的经济活动和能源活动对周边区域的排放有显著影响。UB4是上海城市的最佳表征,2007年UB4内CO₂排放达到1.89亿t,人均排放12.04 t;UB4排放占UB1排放的75.40%,UB1人均排放比UB4人均排放高12%。上海城市化和工业化在空间上的高度重合,导致高排放源集聚于UB4内,形成UB4的高排放特征。个别网格的排放量已经占到区域或者城市总排放量的10%～20%;前10高和前100高排放网格,其累积排放总量分别占据了3个城市范围（UB1、UB3和UB4）总排放量的60%和80%以上。     

多类重力场模型的精度分析及联合确定GPS点正常高的方法

针对模型尺度参数ae对高程异常计算的影响进行数值分析;利用各模型位系数方差信息及区域高精度GPS/水准数据分析国际上新近公布的ITG-GRACE2010S、AIUB-GRACE03S、EIGEN-6C、GOCO02S、DIR_R3、TIM_R3、SPW_R2、gif48以及EGM2008多类重力场模型对应不同谱域的位系数精度;提出采用简单谱组合方法和加权谱组合方法来提高高程异常的计算精度.在试验区域,计算结果表明,各模型的中低阶(约135阶)位系数整体精度要优于EGM2008模型,gif48和EIGEN-6C模型的整体精度均优于同阶次的EGM2008模型,简单谱组合方法和加权谱组合方法均能提高区域模型高程异常精度.