含缺值GNSS基准站坐标序列的非插值小波分析与信号提取

doi:10.11947/j.AGCS.2020.20190163

摘要/Abstract

摘要： 受多种因素影响，GNSS基准站坐标序列通常都含有缺值，传统小波分析需要对缺值数据进行内插或补零处理。本文基于小波系数与时间序列观测数据的重构关系，提出了一种非插值的二进小波变换的最小范数解法，导出了相应的计算式，并严格证明了传统的补零处理算法与本文的最小范数解法等价。最后利用中国地壳运动观测网络一期27个基准站实测数据以及模拟数据进行了验证分析。结果表明，本文的非插值算法与插值算法提取的信号差异较小，27个基准站坐标序列的平均残差中误差仅相差2.01%（North），0.54%（East）和1.26%（Up），两种算法提取的信号之差与信号平均方差比仅相差1.16%（North），0.54%（East）和1.62%（Up）。

关键词: GNSS坐标时间序列, 缺失数据, 小波变换, 信号提取

Abstract: The GNSS position time series are often analyzed by using traditional dyadic wavelet transform, which requires that the time series must be complete. However, missing data inevitably occur in the GNSS position time series due to a variety of causes. In order to extract signals from the incomplete position time series, a modified dyadic wavelet transform algorithm is developed and the corresponding formulas are derived in this paper based on the principle that missing data can be reproduced by its wavelet coefficients. The equivalence between new algorithm and zero-padding algorithm is proved, which indicates that the zero-padding algorithm is essentially a least squares minimum norm solution. Finally, the real position time series of 27 based stations from Crustal Movement Observation Network of China (CMONOC) and simulated data are adopted to verify the validation of the new algorithm, the results show that the difference between the signals extracted by new algorithm and interpolation algorithm is small, with the differences of mean medium errors of 27 base stations are only 2.01%(North), 0.54%(East), 1.26%(Up) and the mean ratios of variance for difference of two signals to the variance for two signals are only 1.16%(North), 0.54%(East), 1.62%(Up).

Key words: GNSS coordinate time series, missing data, wavelet transform, signals extraction

中图分类号:

P207

嵇昆浦, 沈云中. 含缺值GNSS基准站坐标序列的非插值小波分析与信号提取[J]. 测绘学报, 2020, 49(5): 537-546.

JI Kunpu, SHEN Yunzhong. Dyadic wavelet transform and signal extraction of GNSS coordinate time series with missing data[J]. Acta Geodaetica et Cartographica Sinica, 2020, 49(5): 537-546.

参考文献

[1] MING Feng, YANG Yuanxi, ZENG Anmin, et al. Analysis of seasonal signals and long-term trends in the height time series of IGS sites in China[J]. Science China Earth Sciences, 2016, 59(6): 1283-1291.
[2] DONG D, FANG P, BOCK Y, et al. Anatomy of apparent seasonal variations from GPS-derived site position time series[J]. Journal of Geophysical Research: Solid Earth, 2002, 107(B4): ETG-1-ETG 9-16.
[3] RAY J, ALTAMIMI Z, COLLILIEUX X, et al. Anomalous harmonics in the spectra of GPS position estimates[J]. GPS Solutions, 2008, 12(1): 55-64.
[4] ZOU Rong, FREYMUELLER J T, DING Kaihua, et al. Evaluating seasonal loading models and their impact on global and regional reference frame alignment[J]. Journal of Geophysical Research: Solid Earth, 2014, 119(2): 1337-1358.
[5] BOS M S, BASTOS L, FERNANDES R M S. The influence of seasonal signals on the estimation of the tectonic motion in short continuous GPS time-series[J]. Journal of Geodynamics, 2010, 49(3-4): 205-209.
[6] HE Xiaoxing, MONTILLET J P, FERNANDES R, et al. Review of current GPS methodologies for producing accurate time series and their error sources[J]. Journal of Geodynamics, 2017, 106: 12-29.
[7] NIKOLAIDIS R. Observation of geodetic and seismic deformation with the global positioning system[D]. San Diego: University of California, 2002.
[8] WILLIAMS S D P, BOCK Y, FANG Peng, et al. Error analysis of continuous GPS position time series[J]. Journal of Geophysical Research: Solid Earth, 2004, 109(B3): B03412.
[9] LANGBEIN J. Noise in GPS displacement measurements from Southern California and Southern Nevada[J]. Journal of Geophysical Research: Solid Earth, 2008, 113(B5): 1154-1165.
[10] 李昭, 姜卫平, 刘鸿飞, 等. 中国区域IGS基准站坐标时间序列噪声模型建立与分析[J]. 测绘学报, 2012, 41(4): 496-503. LI Zhao, JIANG Weiping, LIU Hongfei, et al. Noise model establishment and analysis of IGS reference station coordinate time series inside China[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(4): 496-503.
[11] 姜卫平, 李昭, 刘鸿飞, 等. 中国区域IGS基准站坐标时间序列非线性变化的成因分析[J]. 地球物理学报, 2013, 56(7): 2228-2237. JIANG Weiping, LI Zhao, LIU Hongfei, et al. Cause analysis of the non-linear variation of the IGS reference station coordinate time series inside China[J]. Chinese Journal of Geophysics 2013, 56(7): 2228-2237.
[12] BENNETT R A. Instantaneous deformation from continuous GPS: contributions from quasi-periodic loads[J]. Geophysical Journal International, 2008, 174(3): 1052-1064.
[13] DAVIS J L, WERNICKE B P, TAMISIEA M E. On seasonal signals in geodetic time series[J]. Journal of Geophysical Research: Solid Earth, 2012, 117(B1): 1104-1113.
[14] CHEN Q, VAN DAM T, SNEEUW N, et al. Singular spectrum analysis for modeling seasonal signals from GPS time series[J]. Journal of Geodynamics, 2013, 72: 25-35.
[15] KUSCHE J, ILK K H, RUDOLPH S, et al. Application of spherical wavelets for regional gravity field recovery: a comparative study[M]. FORSBERG R, FEISSEL M, DIETRICH R. Geodesy on the Move. Berlin: Springer Verlag, 1998.
[16] 徐华君, 柳林涛, 许厚泽, 等. 利用小波分析重力的长期变化[J]. 地球物理学报, 2008, 51(3): 735-742. XU Huajun, LIU Lintao, XU Houze, et al. Wavelet approach to study the secular gravity variation[J]. Chinese Journal of Geophysics, 2008, 51(3): 735-742.
[17] 詹金刚, 王勇, 许厚泽, 等. 我国近海1992—2006年海平面变化的小波分析[J]. 测绘学报, 2008, 37(4): 438-443. ZHAN Jingang, WANG Yong, XU Houze, et al. The wavelet analysis of sea level change in China sea during 1992—2006[J]. Acta Geodaetica et Cartographica Sinica, 2008, 37(4): 438-443.
[18] 章浙涛. 小波分析理论及其在变形监测中的应用研究[D]. 长沙: 中南大学, 2014. ZHANG Zhetao. Theory of wavelet analysis and its application in deformation monitoring[D]. Changsha: Central South University, 2014.
[19] 章浙涛, 朱建军, 卢骏, 等. 小波变换在时间序列特征提取中的应用[J]. 测绘工程, 2014, 23(6): 21-26. ZHANG Zhetao, ZHU Jianjun, LU Jun, et al. Application of wavelet transform to extracting the time series feature[J]. Engineering of Surveying and Mapping, 2014, 23(6): 21-26.
[20] YI Tinghua, LI Hongnan, GU Ming. Experimental assessment of high-rate GPS receivers for deformation monitoring of bridge[J]. Measurement, 2013, 46(1): 420-432.
[21] BOGUSZ J. Geodetic aspects of GPS permanent station non-linearity studies[J]. Acta Geodynamica et Geomaterialia, 2015, 12(4): 180.
[22] KLOS A, BOS M S, BOGUSZ J. Detecting time-varying seasonal signal in GPS position time series with different noise levels[J]. GPS Solutions, 2018, 22: 21.
[23] WU Hao, LI Kui, SHI Wenzhong, et al. A wavelet-based hybrid approach to remove the flicker noise and the white noise from GPS coordinate time series[J]. GPS Solutions, 2015, 19(4): 511-523.
[24] RHIF M, BEN ABBES A, FARAH I R, et al. Wavelet transform application for/in non-stationary time-series analysis: a review[J]. Applied Sciences, 2019, 9(7): 1345.
[25] 符养. 中国大陆现今地壳形变与GPS坐标时间序列分析[D]. 上海: 中国科学院研究生院(上海天文台), 2002. FU Yang. Present-day crustal deformation in China and GPS-derived coordinate time series analysis[D]. Shanghai: Shanghai Astronomical Observatory, Chinese Academy of Sciences, 2002.
[26] 袁林果, 丁晓利, 陈武, 等. 香港GPS基准站坐标序列特征分析[J]. 地球物理学报, 2008, 51(5): 1372-1384. YUAN Linguo, DING Xiaoli, CHEN Wu, et al. Characteristics of daily position time series from the Hong Kong GPS fiducial network[J]. Chinese Journal of Geophysics, 2008, 51(5): 1372-1384.
[27] 田慧, 程鹏飞, 秘金钟. 不同插值方法对CORS高程时间序列的影响分析[J]. 测绘科学, 2013, 38(1): 16-17, 46. TIAN Hui, CHENG Pengfei, BI Jingzhong. Effect analysis of different interpolation methods on height time series of CORS[J]. Science of Surveying and Mapping, 2013, 38(1): 16-17, 46.
[28] SHEN Yunzhong, LI Weiwei, XU Guochang, et al. Spatiotemporal filtering of regional GNSS network’s position time series with missing data using principle component analysis[J]. Journal of Geodesy, 2014, 88(1): 1-12.
[29] SHEN Y, PENG F, LI B. Improved singular spectrum analysis for time series with missing data[J]. Nonlinear Processes in Geophysics, 2015, 22(4): 371-376.
[30] 沈云中, 陶本藻. 实用测量数据处理方法[M]. 北京: 测绘出版社, 2012. SHEN Yunzhong, TAO Benzao. Practical methods for surveying data processing[M]. Beijing: Surveying and Mapping Press, 2012.
[31] WALNUT D F. An introduction to wavelet analysis[M]. New York: Springer Science+Business Media, 2013.
[32] 王超, 朱宏平, 吴巧云, 等. 基于小波多分辨率分析的结构时变参数识别[J]. 华中科技大学学报(自然科学版), 2016, 44(10): 31-35. WANG Chao, ZHU Hongping, WU Qiaoyun, et al. Time-varying parameter identification of structure using wavelet multiresolution analysis[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2016, 44(10): 31-35.
[33] MALLAT S G. Multiresolution representations and wavelets[D]. Philadelphia: University of Pennsylvania, 1988.
[34] 张双成, 李振宇, 何月帆, 等. GNSS高程时间序列周期项的经验模态分解提取[J]. 测绘科学, 2018, 43(8): 80-84, 96. ZHANG Shuangcheng, LI Zhenyu, HE Yuefan, et al. Extracting of periodic component of GNSS vertical time series using EMD[J]. Science of Surveying and Mapping, 2018, 43(8): 80-84, 96.
[35] SHEN Y, LI W. Spatiotemporal signal and noise analysis of GPS position time series of the permanent stations in China[C]//Earth on the Edge: Science for a Sustainable Planet. Melbourne:[s.n.], 2011: 231-237.
[36] YANG Yuanxi. Resilient PNT concept frame[J]. Journal of Geodesy and Geoinformation Science, 2019, 2(3): 1-7. doi:10.11947/j.JGGS.2019.0301.