留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

顾及经纬度方向异性的电离层TEC IDW插值及精度分析

朱永兴 谭述森 明锋 崔先强

朱永兴, 谭述森, 明锋, 崔先强. 顾及经纬度方向异性的电离层TEC IDW插值及精度分析[J]. 武汉大学学报 ● 信息科学版, 2019, 44(11): 1605-1612. doi: 10.13203/j.whugis20180233
引用本文: 朱永兴, 谭述森, 明锋, 崔先强. 顾及经纬度方向异性的电离层TEC IDW插值及精度分析[J]. 武汉大学学报 ● 信息科学版, 2019, 44(11): 1605-1612. doi: 10.13203/j.whugis20180233
ZHU Yongxing, TAN Shusen, MIN Feng, CUI Xianqiang. IDW Ionospheric TEC Interpolation and Accuracy Analysis Considering Latitude and Longitude Anisotropy[J]. Geomatics and Information Science of Wuhan University, 2019, 44(11): 1605-1612. doi: 10.13203/j.whugis20180233
Citation: ZHU Yongxing, TAN Shusen, MIN Feng, CUI Xianqiang. IDW Ionospheric TEC Interpolation and Accuracy Analysis Considering Latitude and Longitude Anisotropy[J]. Geomatics and Information Science of Wuhan University, 2019, 44(11): 1605-1612. doi: 10.13203/j.whugis20180233

顾及经纬度方向异性的电离层TEC IDW插值及精度分析

doi: 10.13203/j.whugis20180233
基金项目: 

国家重点研发计划 2016YFB0501900

国家自然科学基金 41574013

国家自然科学基金 41604013

国家自然科学基金 41604024

国家自然科学基金 41674012

国家自然科学基金 41704035

国家自然科学基金 41874041

详细信息
    作者简介:

    朱永兴, 博士生, 主要从事卫星导航及数据处理研究。chxyzyx_2008@126.com

  • 中图分类号: P228

IDW Ionospheric TEC Interpolation and Accuracy Analysis Considering Latitude and Longitude Anisotropy

Funds: 

The National Key Research and Development Program of China 2016YFB0501900

the National Natural Science Foundation of China 41574013

the National Natural Science Foundation of China 41604013

the National Natural Science Foundation of China 41604024

the National Natural Science Foundation of China 41674012

the National Natural Science Foundation of China 41704035

the National Natural Science Foundation of China 41874041

More Information
    Author Bio:

    ZHU Yongxing, PhD candidate, specializes in satellite navigation and data processing. E-mail:chxyzyx_2008@126.com

  • 摘要: 反距离加权(inverse distance weighting,IDW)是一种简单实用的插值方法。以全球电离层格网(global ionospheric map,GIM)产品为样本,考虑电离层总电子含量(total electron content,TEC)经纬度方向异性,引入经纬度方向异性调节因子,设计了包含等权在内的6种电离层距离计算方案,分析表明,电离层TEC与经度方向相关性高于纬度方向,不同电离层距离计算方案均能有效提高IDW插值精度。采用最优方案IDW插值分析长期插值精度,结果表明,电离层活动剧烈区域(南北纬度20°)连续12 a“两分两至”日前后全球电离层格网(global ionospheric map,GIM)产品插值,最优方案比普通IDW插值精度提升约25%;2014年太阳活动高年“两分两至”日GIM产品插值,地方时14 h后3~5 h电离层活动剧烈时,最优方案插值精度提升明显,插值均方根误差(root mean square,RMS)最大不超过4.0 TECU。
  • 图  1  调节因子r不同取值对1/d的影响

    Figure  1.  Effect of Different Values of Adjustment Factors on 1/d

    图  2  插值范围示意图

    Figure  2.  Schematic Diagram of the Interpolation Range

    图  3  纬度0°电离层插值精度图

    Figure  3.  Ionospheric Interpolation Accuracy of 0° Latitude

    图  4  纬度20°电离层插值精度图

    Figure  4.  Ionospheric Interpolation Accuracy of 20° Latitude

    图  5  纬度40°电离层插值精度图

    Figure  5.  Ionospheric Interpolation Accuracy of 40° Latitude

    图  6  纬度60°电离层插值精度图

    Figure  6.  Ionospheric Interpolation Accuracy of 60° Latitude

    图  7  DOY 080电离层插值精度图

    Figure  7.  Ionospheric Interpolation Accuracy of DOY 080

    图  8  DOY 172电离层插值精度图

    Figure  8.  Ionospheric Interpolation Accuracy of DOY 172

    图  9  DOY 266电离层插值精度图

    Figure  9.  Ionospheric Interpolation Accuracy of DOY 266

    图  10  DOY 356电离层插值精度图

    Figure  10.  Ionospheric Interpolation Accuracy of DOY 356

    图  11  DOY 080电离层插值精度随时间变化

    Figure  11.  Ionospheric Interpolation Accuracy Changes of DOY 080

    图  12  DOY 172电离层插值精度随时间变化

    Figure  12.  Ionospheric Interpolation Accuracy Changes of DOY 172

    图  13  DOY 266电离层插值精度随时间变化

    Figure  13.  Ionospheric Interpolation Accuracy Changes of DOY 266

    图  14  DOY 356电离层插值精度随时间变化

    Figure  14.  Ionospheric Interpolation Accuracy Changes of DOY 356

    表  1  方案(2)插值精度提升/%

    Table  1.   Interpolation Accuracy Improvement of Plan (2)/%

    纬度 μ=2 μ=3 μ=4 μ=5 均值
    22.59 15.44 0.14 -12.94 6.31
    20° 23.67 22.75 1.51 11.01 18.49
    40° 28.23 31.07 28.79 26.27 28.59
    60° 27.91 30.60 29.13 27.90 28.89
    均值 25.60 24.97 18.64 13.06 20.57
    下载: 导出CSV

    表  2  方案(3)插值精度提升/%

    Table  2.   Interpolation Accuracy Improvement of Plan (3) /%

    纬度 μ=2 μ=3 μ=4 μ=5 均值
    13.42 6.62 -6.48 -18.08 -1.13
    20° 4.12 3.63 0.07 -3.29 1.13
    40° 9.30 11.15 8.87 6.47 8.95
    60° 11.10 11.40 10.52 10.31 10.83
    均值 9.49 8.20 3.24 -1.15 4.95
    下载: 导出CSV

    表  3  方案(4)插值精度提升/%

    Table  3.   Interpolation Accuracy Improvement of Plan (4)/%

    纬度 μ=2 μ=3 μ=4 μ=5 均值
    45.03 25.41 4.37 -11.20 15.90
    20° 9.61 8.10 5.30 4.44 6.87
    40° 3.42 2.29 1.19 1.04 1.99
    60° 1.08 0.73 0.43 0.37 0.65
    均值 14.79 9.13 2.82 -1.34 6.35
    下载: 导出CSV

    表  4  方案(5)插值精度提升/%

    Table  4.   Interpolation Accuracy Improvement of Plan (5)/%

    纬度 μ=2 μ=3 μ=4 μ=5 均值
    14.29 11.05 3.54 -2.51 6.59
    20° 17.22 17.59 14.67 12.74 15.56
    40° 19.53 21.96 21.44 21.53 21.11
    60° 21.43 26.06 28.10 30.14 26.43
    均值 18.12 19.16 16.94 15.47 17.42
    下载: 导出CSV

    表  5  相对方案(1)插值精度提升/%

    Table  5.   Interpolation Accuracy Improvement Compare Plan (1)/%

    纬度 μ=2 μ=3 μ=4 μ=5 均值
    45.03 25.41 4.37 -11.20 15.90
    20° 25.16 24.77 18.73 13.23 20.47
    40° 28.45 31.37 29.11 26.59 28.88
    60° 28.13 30.94 29.54 28.36 29.24
    均值 31.69 28.12 20.44 14.25 23.62
    下载: 导出CSV

    表  6  相对方案(2)插值精度提升/%

    Table  6.   Interpolation Accuracy Improvement Compare Plan (2)/%

    纬度 μ=2 μ=3 μ=4 μ=5 均值
    22.59 15.44 0.14 -12.94 6.31
    20° 23.67 22.75 16.51 11.01 18.49
    40° 28.23 31.07 28.79 26.27 28.59
    60° 27.91 30.60 29.13 27.90 28.89
    均值 25.60 24.97 18.64 13.06 20.57
    下载: 导出CSV

    表  7  20°S方案(6)相对方案(1)插值精度提升/%

    Table  7.   Interpolation Accuracy Improvement of Plan (6) Compare Plan (1) in 20°S/%

    DOY 年份 均值
    2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015
    080 31.33 30.85 32.28 35.46 31.85 30.24 29.95 23.13 26.15 28.52 31.46 31.86 30.26
    172 32.41 34.24 28.33 32.36 29.48 27.63 25.63 32.09 30.55 26.85 30.33 31.40 30.11
    266 29.03 33.55 32.57 22.41 20.24 26.45 32.36 30.16 31.01 29.85 27.65 31.23 28.88
    356 12.50 21.95 2.57 20.64 3.82 22.18 2.42 3.67 -2.34 22.08 12.50 21.55 11.96
    均值 26.32 30.15 23.94 27.72 21.35 26.62 22.59 22.26 21.34 26.83 25.49 29.01 25.30
    下载: 导出CSV

    表  8  20°N方案(6)相对方案(1)插值精度提升/%

    Table  8.   Interpolation Accuracy Improvement of Plan (6) Compare Plan (1) in 20°N/%

    DOY 年份 均值
    2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015
    080 30.32 28.07 26.17 28.36 26.33 29.30 27.20 29.25 30.71 28.15 28.33 25.27 28.12
    172 23.45 18.84 24.48 24.80 19.74 15.37 14.23 5.30 20.47 21.35 22.90 13.88 18.73
    266 25.73 27.46 29.66 24.04 16.18 29.46 29.95 26.01 29.25 26.11 26.97 25.62 26.37
    356 19.33 24.34 18.57 27.29 21.79 28.90 28.46 29.37 24.52 28.33 26.89 26.57 25.36
    均值 24.71 24.68 24.72 26.12 21.01 25.76 24.96 22.48 26.24 25.98 26.27 22.84 24.65
    下载: 导出CSV
  • [1] 袁运斌, 霍星亮, 张宝成.近年来我国GNSS电离层延迟精确建模及修正研究进展[J].测绘学报, 2017, 46(10):1 364-1 378 doi:  10.11947/j.AGCS.2017.20170349

    Yuan Yunbin, Huo Xingliang, Zhang Baocheng.Research Progress of Precise Models and Correction for GNSS Ionospheric Delay in China over Recent Years[J].Acta Geodaetica et Cartographica Sinica, 2017, 46(10):1 364-1 378 doi:  10.11947/j.AGCS.2017.20170349
    [2] Schaer S. Mapping and Predicting the Earth's Ionosphere Using the Global Positioning System[D]. Bern: The University of Bern, 1999
    [3] Bhuyan P K, Borah R R. TEC Derived from GPS Network in India and Comparison with the IRI[J]. Advances in Space Research, 2007, 39(5): 830- 840 doi:  10.1016/j.asr.2006.12.042
    [4] 王建平.中国及周边地区电离层TEC短期预报方法研究[D].西安: 西安电子科技大学, 2008

    Wang Jianping. Study on Short-term Prediction Methods of Ionospheric TEC in China and Surrounding Areas [D]. Xi?an: Xi?an University of Electronic Science and Technology, 2008
    [5] 黄玲, 章红平, 徐培亮, 等.中国区域VTEC模型Kriging算法研究[J].武汉大学学报∙信息科学版, 2016, 41(6):729-737 http://ch.whu.edu.cn/CN/abstract/abstract5456.shtml

    Huang Ling, Zhang Hongping, Xu Peiliang, et al. VTEC Modeling with Kriging Algorithm over China Area[J]. Geomantics and Information Science of Wuhan University, 2016, 41(6): 729-737 http://ch.whu.edu.cn/CN/abstract/abstract5456.shtml
    [6] 毛田, 万卫星, 孙凌峰.用Kriging方法构建中纬度区域电离层TEC地图[J].空间科学学报, 2007(4):279-285 doi:  10.3969/j.issn.0254-6124.2007.04.003

    Mao Tian, Wan Weixing, Sun Lingfeng. Central and Northern China TEC Map Using the Kriging Method[J]. Chin J Space Sci, 2007(4):279-285 doi:  10.3969/j.issn.0254-6124.2007.04.003
    [7] 樊子德, 李佳霖, 邓敏.顾及多因素影响的自适应IDW插值方法[J].武汉大学学报∙信息科学版, 2016, 41(6):842-847 http://ch.whu.edu.cn/CN/abstract/abstract5473.shtml

    Fan Zide, Li Jialin, Deng Min. An Adaptive Inverse-Distance Weighting Spatial Interpolation Method with the Consideration of Multiple Factors[J]. Geomantics and Information Science of Wuhan University, 2016, 41(6):842-847 http://ch.whu.edu.cn/CN/abstract/abstract5473.shtml
    [8] 薛树强, 杨元喜.广义IDW空间推估法[J].武汉大学学报∙信息科学版, 2013, 38(12):1 435-1 439 http://ch.whu.edu.cn/CN/abstract/abstract2825.shtml

    Xue Shuqiang, Yang Yuanxi. Generalized Inverse Distance Weighting Method for Spatial Interpolation[J]. Geomantics and Information Science of Wuhan University, 2013, 38(12):1 435-1 439 http://ch.whu.edu.cn/CN/abstract/abstract2825.shtml
    [9] 黄谟涛, 管铮, 欧阳永忠.海洋平均重力异常计算与精度估计[J].武汉大学学报∙信息科学版, 1995, 20(4):327-331 http://ch.whu.edu.cn/CN/abstract/abstract4062.shtml

    Huang Motao, Guan Zheng, Ouyang Yongzhong. Calculation and Accuracy Estimation of Marine Mean Free-Air Gravity Anomaly[J]. Geomantics and Information Science of Wuhan University, 1995, 20(4):327-331 http://ch.whu.edu.cn/CN/abstract/abstract4062.shtml
    [10] 熊雄, 叶小岭, 张颖超, 等.基于空间观测差异的地面气温资料质量控制算法研究[J].地球物理学报, 2017, 60(3):912-923 http://www.cqvip.com/QK/94718X/201703/671496927.html

    Xiong Xiong, Ye Xiaoling, Zhang Yingchao, et al. A Quality Control Method for the Surface Temperature Based on the Spatial Observation Diversity[J]. Chinese Journal of Geophysics, 2017, 60(3): 912-923 http://www.cqvip.com/QK/94718X/201703/671496927.html
    [11] 张锦明, 游雄, 万刚.DEM插值参数优选的试验研究[J].测绘学报, 2014, 43(2):178-192 http://www.cnki.com.cn/Article/CJFDTotal-CHXB201402013.htm

    Zhang Jinming, You Xiong, Wan Gang. Experimental Research on Optimization of DEM Interpolation Parameters[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(2):178-192 http://www.cnki.com.cn/Article/CJFDTotal-CHXB201402013.htm
    [12] 安家春, 章迪, 杜玉军, 等.极区电离层梯度的特性分析[J].武汉大学学报∙信息科学版, 2014, 39(1): 75-79 http://ch.whu.edu.cn/CN/abstract/abstract2847.shtml

    An Jiachun, Zhang Di, Du Yujun, et al. Ionospheric Electron Density Gradients of a Polar Region Using GPS[J]. Geomatics and Information Science of Wuhan University, 2014, 39(1): 75-79 http://ch.whu.edu.cn/CN/abstract/abstract2847.shtml
    [13] 冯来平, 朱永兴.北斗区域卫星导航系统电离层建模精度分析[J].海洋测绘, 2016, 36(5):47-50 doi:  10.3969/j.issn.1671-3044.2016.05.012

    Feng Laiping, Zhu Yongxing. Precision of Regional Ionosphere Model Based on BeiDou Navigation Satellite System[J]. Hydrographic Surveying and Charting, 2016, 36(5):47-50 doi:  10.3969/j.issn.1671-3044.2016.05.012
    [14] 顾春雷, 张毅, 徐如刚, 等.基于虚拟日变台进行地磁矢量数据日变通化方法[J].地球物理学报, 2013, 56(3):834-841 http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201303012

    Gu Chunlei, Zhang Yi, Xu Rugang, et al. The Method of Geomagnetic Diurnal Correction to Vector Geomagnetic Data by Virtual Diurnal Variation Station[J]. Chinese Journal of Geophysics, 2013, 56(3): 834-841 http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201303012
    [15] 张清华, 隋立芬, 贾小林, 等.北斗卫星导航系统空间信号误差统计分析[J].武汉大学学报∙信息科学版, 2014, 39(3):271-274 http://ch.whu.edu.cn/CN/abstract/abstract2913.shtml

    Zhang Qinghua, Sui Lifen, Jia Xiaolin, et al. SIS Error Statistical Analysis of BeiDou Satellite Navigation System[J]. Geomantics and Information Science of Wuhan University, 2014, 39(3): 271-274 http://ch.whu.edu.cn/CN/abstract/abstract2913.shtml
    [16] Lam N S.Spatial Interpolation Methods: A Review[J].American Cartographer, 1983, 10(2): 129- 149 doi:  10.1559/152304083783914958
    [17] Declercq F A N.Interpolation Methods for Scattered Sample Data:Accuracy, Spatial Patterns, Processing Time[J].Cartography and Geographic Information Systems, 1996, 23(3):128-144 doi:  10.1559/152304096782438882
    [18] Chen F W, Liu C W. Estimation of Spatial Rainfall Distribution Using Inverse Distance Weighting (IDW) in the Middle of Taiwan[J]. Paddy and Water Environment, 2012, 10(3):209-222 doi:  10.1007/s10333-012-0319-1
    [19] 赵海山, 杨力, 徐世依.太阳活动高低年电离层TEC变化特性分析[J].导航定位学报, 2017, 5(1):24-30 http://d.old.wanfangdata.com.cn/Periodical/dhdwxb201701006

    Zhao Haishan, Yang Li, Xu Shiyi. Analysis of Ionospheric TEC Variation Characteristics in Solar Activity Years[J].Journal of Navigation and Positioning, 2017, 5(1):24-30 http://d.old.wanfangdata.com.cn/Periodical/dhdwxb201701006
    [20] 李子申, 王宁波, 李敏, 等.国际GNSS服务组织全球电离层TEC格网精度评估与分析[J].地球物理学报, 2017, 60(10):3 718-3 729 doi:  10.6038/cjg20171003

    Li Zishen, Wang Ningbo, Li Min, et al. Evaluation and Analysis of the Global Ionospheric TEC Map in the Frame of International GNSS Services[J]. Chinese Journal of Geophysics, 2017, 60(10): 3 718- 3 729 doi:  10.6038/cjg20171003
  • [1] 郭斐, 吴维旺, 张小红, 刘万科.  Android智能手机实时精密单点定位软件实现及精度分析 . 武汉大学学报 ● 信息科学版, 2021, 46(7): 1053-1062. doi: 10.13203/j.whugis20200527
    [2] 朱永兴, 谭述森, 任夏, 贾小林.  GNSS全球广播电离层模型精度分析 . 武汉大学学报 ● 信息科学版, 2020, 45(5): 768-775. doi: 10.13203/j.whugis20180439
    [3] 朱明晨, 胡伍生, 王来顺.  GPT2w模型在中国区域的精度检验与分析 . 武汉大学学报 ● 信息科学版, 2019, 44(9): 1304-1311. doi: 10.13203/j.whugis20170387
    [4] 孔垚, 孙保琪, 杨旭海, 曹芬, 何战科, 杨海彦.  利用SLR数据进行北斗卫星广播星历精度分析 . 武汉大学学报 ● 信息科学版, 2017, 42(6): 831-837. doi: 10.13203/j.whugis20140856
    [5] 樊子德, 李佳霖, 邓敏.  顾及多因素影响的自适应反距离加权插值方法 . 武汉大学学报 ● 信息科学版, 2016, 41(6): 842-847. doi: 10.13203/j.whugis20140658
    [6] 贾剑钢, 栾威, 申文斌.  iGrav-007超导重力仪格值的精密测定 . 武汉大学学报 ● 信息科学版, 2015, 40(10): 1366-1370,1376. doi: 10.13203/j.whugis20130780
    [7] 段兵兵, 王解先, 王成.  融合速度和方向信息的船舶位置高精度预报 . 武汉大学学报 ● 信息科学版, 2015, 40(3): 422-426.
    [8] 唐卫明, 金蕾, 徐坤.  北斗卫星信号实时单站电离层估计算法及性能分析 . 武汉大学学报 ● 信息科学版, 2014, 39(11): 1318-1322.
    [9] 曹建峰, 胡松杰, 黄勇, 刘磊.  嫦娥二号卫星日地拉格朗日L_2点探测轨道定轨分析 . 武汉大学学报 ● 信息科学版, 2013, 38(9): 1029-1033.
    [10] 肖金群, 李志伟, 汪长城, 丁晓利.  相位补偿SAR差分干涉提取山区DEM算法及精度分析 . 武汉大学学报 ● 信息科学版, 2012, 37(3): 334-338.
    [11] 楼良盛, 刘思伟, 周瑜.  机载InSAR系统精度分析 . 武汉大学学报 ● 信息科学版, 2012, 37(1): 63-67.
    [12] 耿迅, 徐水平, 龚志辉, 袁军.  ADS40线阵影像空中三角测量数据处理与精度分析方法 . 武汉大学学报 ● 信息科学版, 2011, 36(7): 776-779.
    [13] 李一鹤, 沈云中.  GPS观测值的时间相关性对基线解算影响 . 武汉大学学报 ● 信息科学版, 2011, 36(4): 427-430.
    [14] 张卡, 盛业华, 叶春, 梁诚.  车载三维数据采集系统的绝对标定及精度分析 . 武汉大学学报 ● 信息科学版, 2008, 33(1): 55-59.
    [15] 刘国林, 郝晓光, 薛怀平, 独知行.  影响InSAR测高精度因素的相关性分析 . 武汉大学学报 ● 信息科学版, 2007, 32(1): 55-58.
    [16] 张永军, 张勇.  大重叠度影像的相对定向与前方交会精度分析 . 武汉大学学报 ● 信息科学版, 2005, 30(2): 126-130.
    [17] 孙凤华, 吴晓平, 张传定.  中国陆海任意点垂线偏差的快速确定及精度分析 . 武汉大学学报 ● 信息科学版, 2005, 30(1): 42-46.
    [18] 何平安, 王基尧, 翁兴涛, 李松.  角锥棱镜直角误差斜入射法检测的理论分析 . 武汉大学学报 ● 信息科学版, 2000, 25(5): 449-452.
    [19] 黄谟涛, 管铮, 欧阳永忠.  中国地区1°×1°点质量解算与精度分析 . 武汉大学学报 ● 信息科学版, 1995, 20(3): 257-262.
    [20] 王波, 刘路, 鄢建国, 高梧桐.  小行星光学定轨软件研制和数据处理分析 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20200195
  • 加载中
图(14) / 表(8)
计量
  • 文章访问数:  2699
  • HTML全文浏览量:  158
  • PDF下载量:  211
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-06-11
  • 刊出日期:  2019-11-05

顾及经纬度方向异性的电离层TEC IDW插值及精度分析

doi: 10.13203/j.whugis20180233
    基金项目:

    国家重点研发计划 2016YFB0501900

    国家自然科学基金 41574013

    国家自然科学基金 41604013

    国家自然科学基金 41604024

    国家自然科学基金 41674012

    国家自然科学基金 41704035

    国家自然科学基金 41874041

    作者简介:

    朱永兴, 博士生, 主要从事卫星导航及数据处理研究。chxyzyx_2008@126.com

  • 中图分类号: P228

摘要: 反距离加权(inverse distance weighting,IDW)是一种简单实用的插值方法。以全球电离层格网(global ionospheric map,GIM)产品为样本,考虑电离层总电子含量(total electron content,TEC)经纬度方向异性,引入经纬度方向异性调节因子,设计了包含等权在内的6种电离层距离计算方案,分析表明,电离层TEC与经度方向相关性高于纬度方向,不同电离层距离计算方案均能有效提高IDW插值精度。采用最优方案IDW插值分析长期插值精度,结果表明,电离层活动剧烈区域(南北纬度20°)连续12 a“两分两至”日前后全球电离层格网(global ionospheric map,GIM)产品插值,最优方案比普通IDW插值精度提升约25%;2014年太阳活动高年“两分两至”日GIM产品插值,地方时14 h后3~5 h电离层活动剧烈时,最优方案插值精度提升明显,插值均方根误差(root mean square,RMS)最大不超过4.0 TECU。

English Abstract

朱永兴, 谭述森, 明锋, 崔先强. 顾及经纬度方向异性的电离层TEC IDW插值及精度分析[J]. 武汉大学学报 ● 信息科学版, 2019, 44(11): 1605-1612. doi: 10.13203/j.whugis20180233
引用本文: 朱永兴, 谭述森, 明锋, 崔先强. 顾及经纬度方向异性的电离层TEC IDW插值及精度分析[J]. 武汉大学学报 ● 信息科学版, 2019, 44(11): 1605-1612. doi: 10.13203/j.whugis20180233
ZHU Yongxing, TAN Shusen, MIN Feng, CUI Xianqiang. IDW Ionospheric TEC Interpolation and Accuracy Analysis Considering Latitude and Longitude Anisotropy[J]. Geomatics and Information Science of Wuhan University, 2019, 44(11): 1605-1612. doi: 10.13203/j.whugis20180233
Citation: ZHU Yongxing, TAN Shusen, MIN Feng, CUI Xianqiang. IDW Ionospheric TEC Interpolation and Accuracy Analysis Considering Latitude and Longitude Anisotropy[J]. Geomatics and Information Science of Wuhan University, 2019, 44(11): 1605-1612. doi: 10.13203/j.whugis20180233
  • 电离层效应是卫星导航系统最严重的误差源之一,高精度电离层延迟信息对提高单频用户时延修正精度和双频/多频用户精密定位收敛速度具有重要作用[1-4]。获取整个空间高精度电离层延迟信息一般需要通过插值方法计算。电离层总电子含量(total electron content,TEC)插值是指基于区域内有限个电离层穿刺点(ionospheric pierce point,IPP)的TEC值,采用一定的插值算法获取区域内任意未知点的电离层TEC,它是电离层研究的重要内容[1, 5]。利用全球卫星导航系统(global navigation satellite system, GNSS)双频电离层TEC观测量构建格网电离层、缺失格网点数据补插、大尺度格网数据的区域加密等都离不开电离层插值[5-7]。反距离加权(inverse distance weighting,IDW)插值方法是一种应用广泛、原理简单、运算量小的空间插值方法,已广泛应用气象、水文、地质、地磁、重力等领域[8-10]。但是,普通IDW插值方法在电离层TEC插值中,由于精度较低导致应用受限。本文通过分析IDW插值方法的插值参数,构造一种考虑经纬度方向异性的IDW插值方法,提高这种简单实用方法的电离层TEC插值精度。

    插值参数是构成插值方法的基本元素,插值精度取决于插值参数对要素空间变异性和相关性的反映[8]。IDW插值方法的插值参数包括插值范围和插值核函数。其中,插值范围是以待插值点为中心、一定半径的邻域范围;插值核函数是反映两点间空间关系及其对待插值点贡献的函数[11]。本文考虑到电离层TEC分布的经纬度方向异性,引入经纬度方向异性调节因子r,设计了包含等权在内的6种电离层距离计算方案,选用一组全球电离层格网(global ionospheric map, GIM)产品,分析不同方案的插值精度,获取最优方案的IDW插值。之后,针对南北纬度20°电离层活动剧烈区域,采用连续12 a“两分两至”日前后GIM数据分析最优方案的插值精度提升情况,验证电离层TEC的经纬度方向异性特点及本文改进IDW方案的正确性;采用2014年太阳活动高年“两分两至”日数据,分析IDW插值精度与地方时的相关性。

    • IDW插值方法也称为距离乘方倒数法,基本原理是利用区域内有限个电离层IPP Z(libi)(i = 1,2…n)的电离层观测值Ti,求解区域内任意未知点Z(l0b0)的电离层T0,公式为:

      $$ \left\{ {\begin{array}{*{20}{l}} {{T_0} = \frac{{\sum\limits_{i = 1}^n {\frac{1}{{d_i^\mu }}} {T_i}}}{{\sum\limits_{i = 1}^n {\frac{1}{{d_i^\mu }}} }}}\\ {{d_i} = \sqrt {{{\left( {{l_i} - {l_0}} \right)}^2} + {{\left( {{b_i} - {b_0}} \right)}^2}} } \end{array}} \right. $$ (1)

      式中,libi为IPP经纬度;di为采样点和待插值点间距离;μ为衰减指数。

      IDW插值方法忽略了电离层TEC在经纬度方向的变化梯度。文献[12]表明,中纬地区磁平静期间,纬度方向梯度大于经度方向,磁暴期间梯度差异更加明显。因此,考虑电离层TEC经纬度方向梯度,本文对采样点与待插值点间距离函数d进行修正,增加r,公式为:

      $$ {{d}_{i}}=\sqrt{{{({{l}_{i}}-{{l}_{0}})}^{2}}+r{{({{b}_{i}}-{{b}_{0}})}^{2}}} $$ (2)

      IDW插值的基础是采样点和待插值点间的空间相关性,如果插值区域各向均匀,则相关性只与距离有关。由于电离层TEC经纬度方向变化梯度不同,纬度方向电离层TEC变化大,空间相关性低;而经度方向电离层TEC变化小,空间相关性高。对于与待插值点相同距离的样本点,显然当调节因子r = 1时,即为普通IDW插值;当r > 1时,纬度方向距离增大,也即该方向贡献减小;当r < 1时,纬度方向距离减小,也即该方向贡献增大。以与待插值点距离为10°的所有样本点为例,分析取值1.5和0.5时距离插值因子1/d的变化,如图 1所示。

      图  1  调节因子r不同取值对1/d的影响

      Figure 1.  Effect of Different Values of Adjustment Factors on 1/d

      图 1中半圆蓝线为与待插值点(0, 0)距离为10°的样本点分布位置,红色星号为r取值1.5时距离插值因子,对比r取值1.0(黑色加号),纬度方向样本点对待插值点的贡献随距离增大而减小(下凹);粉色倒立三角形为r取值0.5时距离插值因子,对比r取值1.0,纬度方向样本点对待插值点的贡献随距离增大而增大(上凸)。当r ≠ 1时,已知点对待插值点贡献由等权直线分布变为非等权分布。

    • 由式(1)可知,插值参数d是IDW插值方法的基本元素,决定着插值精度。通过优化IDW插值参数有望改善其插值效果。

      插值范围决定着采样点数量、分布,进而影响插值精度。文献[13]认为GNSS单站电离层IPP的覆盖范围约为30°×30°,选择10°×15°和20°×20°插值范围进行插值分析,如图 2所示。

      图  2  插值范围示意图

      Figure 2.  Schematic Diagram of the Interpolation Range

      由插值原理可知,插值核函数主要包括距离函数d和距离衰减指数μ两部分。

      1)距离函数d。主要考虑经纬度方向异性调节因子r,设lb为经纬度均值,T为电离层TEC均值。设计以下5种计算方案[5, 14-15]

      (1)取r=1。

      (2)考虑电离层TEC在经纬度方向变化梯度,计算公式为:

      $$ r=\sum\limits_{{}}^{i\ne j}{\frac{\left( {{l}_{i}}-{{l}_{j}} \right)}{\left( {{T}_{i}}-{{T}_{j}} \right)}/}\sum\limits_{{}}^{i\ne j}{\frac{\left( {{b}_{i}}-{{b}_{j}} \right)}{\left( {{T}_{i}}-{{T}_{j}} \right)}} $$ (3)

      (3)考虑电离层TEC与经纬度相关系数,计算公式为:

      $$ r=\frac{\sum\limits_{i=1}^{n}{\left( {{l}_{i}}-\bar{l} \right)\left( {{T}_{i}}-\bar{T} \right)}}{\sqrt{\sum\limits_{i=1}^{n}{{{({{l}_{i}}-\bar{l})}^{2}}}}\sqrt{\sum\limits_{i=1}^{n}{{{({{T}_{i}}-\bar{T})}^{2}}}}}/\frac{\sum\limits_{i=1}^{n}{\left( {{b}_{i}}-\bar{b} \right)\left( {{T}_{i}}-\bar{T} \right)}}{\sqrt{\sum\limits_{i=1}^{n}{{{({{b}_{i}}-\bar{b})}^{2}}}}\sqrt{\sum\limits_{i=1}^{n}{{{({{T}_{i}}-\bar{T})}^{2}}}}} $$ (4)

      (4)考虑电离层TEC与待插值点地理纬度φ的关系,计算公式为:

      $$ r = {\rm{sin}}\varphi $$ (5)

      (5)参考Kriging插值中电离层距离概念,取经验值r = 0.2。

      2)距离衰减指数μ。随着μ增大,样本点对待插值点影响敏感程度增加,表现为距离近样本点影响显著增加,反之则显著降低,一般不宜采用过大或过小值。文献[16]认为取1或2比较合适,文献[17]认为取2能取得更好实验效果,文献[18]则认为取0~5之间插值效果最优,综上,本文将μ取值范围设为2~5。

    • 选择太阳活动高年2014年,年积日(day of year, DOY)001天通用协调时(universal time coordinated,UTC)2时,国际GNSS服务组织(International GNSS Service,IGS)发布的全球电离层格网(global ionospheric map,GIM)产品(产品精度2.5~3.5总电离子含量单位(total electron content unit, TECU),时间分辨率2 h,空间分辨率2.5°×5°)作为样本,选择北半球纬度0°~ 60°、间隔20°,经度0°~ 360°、间隔15°所有格网点为待插值点,分析插值精度[19-20]。样本点选取范围是以待插值点为中心10°×15°和20°×20°的矩形区域,距离衰减指数μ取值2~5,距离函数di为§1.2的5种方案。计算待插值点插值结果与已知垂直总电子含量(vertical total electron content,VTEC)之差dvteci,统计同一纬度上所有格网点dvteci的均方根误差(root mean square,RMS)值,验证不同方案插值精度

    • 1)插值精度分析。不同插值范围、不同插值函数的插值精度如图 3~6所示:(1)高纬度区域插值精度高于中低纬度区域,纬度20°区域插值误差较大;(2)10°×15°插值范围的插值精度高于20°×20°插值范围;(3)考虑经纬度方向变化梯度,插值结果总体好于其他方案;(4)随着距离衰减指数μ增大,插值精度呈提高趋势,不同距离计算函数插值精度趋于相近。

      图  3  纬度0°电离层插值精度图

      Figure 3.  Ionospheric Interpolation Accuracy of 0° Latitude

      图  4  纬度20°电离层插值精度图

      Figure 4.  Ionospheric Interpolation Accuracy of 20° Latitude

      图  5  纬度40°电离层插值精度图

      Figure 5.  Ionospheric Interpolation Accuracy of 40° Latitude

      图  6  纬度60°电离层插值精度图

      Figure 6.  Ionospheric Interpolation Accuracy of 60° Latitude

      2)精度提升分析。针对10°×15°插值范围,分析各种方案插值精度提升(见表 1~4),总体上,各种方案相比方案(1)插值精度均有提升;采用方案(2)时,插值精度提升最高,平均提高20.57%;在纬度0°,方案(4)提升最高,达到15.90%,中高纬度区域方案最优,提升约为28%。

      表 1  方案(2)插值精度提升/%

      Table 1.  Interpolation Accuracy Improvement of Plan (2)/%

      纬度 μ=2 μ=3 μ=4 μ=5 均值
      22.59 15.44 0.14 -12.94 6.31
      20° 23.67 22.75 1.51 11.01 18.49
      40° 28.23 31.07 28.79 26.27 28.59
      60° 27.91 30.60 29.13 27.90 28.89
      均值 25.60 24.97 18.64 13.06 20.57

      表 2  方案(3)插值精度提升/%

      Table 2.  Interpolation Accuracy Improvement of Plan (3) /%

      纬度 μ=2 μ=3 μ=4 μ=5 均值
      13.42 6.62 -6.48 -18.08 -1.13
      20° 4.12 3.63 0.07 -3.29 1.13
      40° 9.30 11.15 8.87 6.47 8.95
      60° 11.10 11.40 10.52 10.31 10.83
      均值 9.49 8.20 3.24 -1.15 4.95

      表 3  方案(4)插值精度提升/%

      Table 3.  Interpolation Accuracy Improvement of Plan (4)/%

      纬度 μ=2 μ=3 μ=4 μ=5 均值
      45.03 25.41 4.37 -11.20 15.90
      20° 9.61 8.10 5.30 4.44 6.87
      40° 3.42 2.29 1.19 1.04 1.99
      60° 1.08 0.73 0.43 0.37 0.65
      均值 14.79 9.13 2.82 -1.34 6.35

      表 4  方案(5)插值精度提升/%

      Table 4.  Interpolation Accuracy Improvement of Plan (5)/%

      纬度 μ=2 μ=3 μ=4 μ=5 均值
      14.29 11.05 3.54 -2.51 6.59
      20° 17.22 17.59 14.67 12.74 15.56
      40° 19.53 21.96 21.44 21.53 21.11
      60° 21.43 26.06 28.10 30.14 26.43
      均值 18.12 19.16 16.94 15.47 17.42

      进一步分析:(1)IDW插值方法,样本点距离待插值点距离越近,相关性越大,调节衰减指数μ可以有效控制距离较远样本点对待插值点的影响;(2)电离层有明显的经纬度方向异性,考虑经纬度方向变化梯度的插值精度较高;(3)纬度20°区域属于电离层活动剧烈区域,插值误差较大。

    • §2.1.2分析表明,方案(2)并不是最优的,如在纬度0°时,方案(4)提升更高,即是r取值0时,插值效果最优。分析认为:(1)低纬度区域电离层TEC变化剧烈,有限样本不足以反映电离层TEC经纬度方向异性;(2)样本点与待插值点距离相等时,经度方向相关性高于纬度方向。为此,本节设计一种综合方案(方案(6)),该方案计算距离函数di时,r取值为方案(2)~(5)的最小值。从表 5~6可以看出:相对于方案(1),平均提升23.62%,相对于方案(2),平均提升20.57%;μ取值2~3时,精度提升较大。

      表 5  相对方案(1)插值精度提升/%

      Table 5.  Interpolation Accuracy Improvement Compare Plan (1)/%

      纬度 μ=2 μ=3 μ=4 μ=5 均值
      45.03 25.41 4.37 -11.20 15.90
      20° 25.16 24.77 18.73 13.23 20.47
      40° 28.45 31.37 29.11 26.59 28.88
      60° 28.13 30.94 29.54 28.36 29.24
      均值 31.69 28.12 20.44 14.25 23.62

      表 6  相对方案(2)插值精度提升/%

      Table 6.  Interpolation Accuracy Improvement Compare Plan (2)/%

      纬度 μ=2 μ=3 μ=4 μ=5 均值
      22.59 15.44 0.14 -12.94 6.31
      20° 23.67 22.75 16.51 11.01 18.49
      40° 28.23 31.07 28.79 26.27 28.59
      60° 27.91 30.60 29.13 27.90 28.89
      均值 25.60 24.97 18.64 13.06 20.57

      综上分析:(1)方案(6)计算r的插值精度优于其他5种方案;(2)当其他插值参数相同时,插值范围10°×15°结果好于20°×20°,即样本点充足时,可以适当缩小样本点选取范围或者通过增大距离衰减指数μ 来降低距离较远样本点的影响;(3)电离层TEC经度方向相关性高于纬度方向,样本点充足时可以适当选用经度差更小的点或通过调节因子r来降低纬度差较大的样本点的影响。

    • 选用2004—2015年共12 a“两分两至”日前后UTC 6时GIM产品(对应经度120°E的地方时为14时),插值范围10°×15°,μ 取值3,采用方案(6)的IDW插值方法,分析电离层活动剧烈区域(南北纬度20°)的插值精度,并与方案(1)插值结果对比。从图 7~10表 7~8可以看出:(1)电离层活动剧烈区域,方案(6)相对于方案(1)精度提升约25%,总体插值精度约为2.0~3.0 TECU;(2)春分日和秋分日前后电离层活动相对剧烈,电离层插值精度相对较差,采用方案(6)精度提升明显;(3)太阳活动剧烈年份,IDW插值精度相对较差,但方案(6)相对于方案(1)在不同年份精度提升稳定,认为是GIM样本精度所致。

      图  7  DOY 080电离层插值精度图

      Figure 7.  Ionospheric Interpolation Accuracy of DOY 080

      图  8  DOY 172电离层插值精度图

      Figure 8.  Ionospheric Interpolation Accuracy of DOY 172

      图  9  DOY 266电离层插值精度图

      Figure 9.  Ionospheric Interpolation Accuracy of DOY 266

      图  10  DOY 356电离层插值精度图

      Figure 10.  Ionospheric Interpolation Accuracy of DOY 356

      表 7  20°S方案(6)相对方案(1)插值精度提升/%

      Table 7.  Interpolation Accuracy Improvement of Plan (6) Compare Plan (1) in 20°S/%

      DOY 年份 均值
      2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015
      080 31.33 30.85 32.28 35.46 31.85 30.24 29.95 23.13 26.15 28.52 31.46 31.86 30.26
      172 32.41 34.24 28.33 32.36 29.48 27.63 25.63 32.09 30.55 26.85 30.33 31.40 30.11
      266 29.03 33.55 32.57 22.41 20.24 26.45 32.36 30.16 31.01 29.85 27.65 31.23 28.88
      356 12.50 21.95 2.57 20.64 3.82 22.18 2.42 3.67 -2.34 22.08 12.50 21.55 11.96
      均值 26.32 30.15 23.94 27.72 21.35 26.62 22.59 22.26 21.34 26.83 25.49 29.01 25.30

      表 8  20°N方案(6)相对方案(1)插值精度提升/%

      Table 8.  Interpolation Accuracy Improvement of Plan (6) Compare Plan (1) in 20°N/%

      DOY 年份 均值
      2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015
      080 30.32 28.07 26.17 28.36 26.33 29.30 27.20 29.25 30.71 28.15 28.33 25.27 28.12
      172 23.45 18.84 24.48 24.80 19.74 15.37 14.23 5.30 20.47 21.35 22.90 13.88 18.73
      266 25.73 27.46 29.66 24.04 16.18 29.46 29.95 26.01 29.25 26.11 26.97 25.62 26.37
      356 19.33 24.34 18.57 27.29 21.79 28.90 28.46 29.37 24.52 28.33 26.89 26.57 25.36
      均值 24.71 24.68 24.72 26.12 21.01 25.76 24.96 22.48 26.24 25.98 26.27 22.84 24.65
    • 选用2014年太阳活动高年“两分两至”日GIM产品,采用方案(6)IDW插值方法,分析电离层活动剧烈区域(南北纬20°),经度范围90°E~120°E(间隔10°,中心经度105°,对应地方时UTC 7时),统计同一时刻所有待插值点插值结果与已知值之差dvteci的RMS值,分析随地方时变化情况。插值范围为10°×15°,μ取值为3,比较方案(6)和方案(1)的插值精度。从图 11~14可以看出:(1)方案(1)和方案(6)插值精度均随时间变化明显,在地方时14 h(UTC 7时)后3~5 h(冬至日时间滞后),插值精度明显变差;(2)方案(6)相对于方案(1)在午后电离层活动剧烈时,插值精度提升明显,RMS最大不超过4.0 TECU。

      图  11  DOY 080电离层插值精度随时间变化

      Figure 11.  Ionospheric Interpolation Accuracy Changes of DOY 080

      图  12  DOY 172电离层插值精度随时间变化

      Figure 12.  Ionospheric Interpolation Accuracy Changes of DOY 172

      图  13  DOY 266电离层插值精度随时间变化

      Figure 13.  Ionospheric Interpolation Accuracy Changes of DOY 266

      图  14  DOY 356电离层插值精度随时间变化

      Figure 14.  Ionospheric Interpolation Accuracy Changes of DOY 356

    • 本文从IDW插值方法插值范围、插值核函数等入手,考虑电离层TEC分布经纬度方向异性,引入经纬度方向异性调节因子r ,设计多种电离层距离计算方案,选用一组GIM数据分析不同方案插值精度,获取最优方案IDW插值,之后,采用2004—2015年部分GIM数据,分析最优方案IDW插值精度,得到一些有益结论,具体如下:

      1)插值范围影响电离层TEC插值精度。基于GIM数据插值,当其他插值参数相同时,插值范围10°×15°结果好于20°×20°,插值范围对待插值点的影响可以通过调整样本选取范围或距离衰减指数μ来调节。

      2)电离层TEC经纬度方向异性影响插值精度。经度方向相关性高于纬度方向,当样本点充足时,可以适当选用经度差更小的样本点或通过调节因子r来降低纬度差较大的样本点的影响,提高插值精度。

      3)综合方案(方案(6))能较好地反映经纬度方向异性特点,插值精度优于其他方案。采用该方案的IDW插值分析长期插值精度,电离层活动剧烈区域(南北纬20°)连续12 a“两分两至”日前后的GIM数据插值,最优方案比普通IDW插值精度提升约25%;2014年太阳活动高年“两分两至”日GIM数据插值,在地方时14时后3~5 h电离层活动剧烈时,最优方案插值精度提升明显,插值误差RMS最大不超过4.0 TECU。

参考文献 (20)

目录

    /

    返回文章
    返回