留言板

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

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

基于俯仰轴相交和球拟合的天线旋转中心测量方法

李宗春 路志勇 郭迎钢 张冠宇 何华 冯其强 陈少卿 王俊威

李宗春, 路志勇, 郭迎钢, 张冠宇, 何华, 冯其强, 陈少卿, 王俊威. 基于俯仰轴相交和球拟合的天线旋转中心测量方法[J]. 武汉大学学报 ● 信息科学版, 2019, 44(10): 1449-1456. doi: 10.13203/j.whugis20180064
引用本文: 李宗春, 路志勇, 郭迎钢, 张冠宇, 何华, 冯其强, 陈少卿, 王俊威. 基于俯仰轴相交和球拟合的天线旋转中心测量方法[J]. 武汉大学学报 ● 信息科学版, 2019, 44(10): 1449-1456. doi: 10.13203/j.whugis20180064
LI Zongchun, LU Zhiyong, GUO Yinggang, ZHANG Guanyu, HE Hua, FENG Qiqiang, CHEN Shaoqing, WANG Junwei. Two Measuring Methods for Antenna Rotating Center Based on Elevation Axes Intersection and Spherical Fitting[J]. Geomatics and Information Science of Wuhan University, 2019, 44(10): 1449-1456. doi: 10.13203/j.whugis20180064
Citation: LI Zongchun, LU Zhiyong, GUO Yinggang, ZHANG Guanyu, HE Hua, FENG Qiqiang, CHEN Shaoqing, WANG Junwei. Two Measuring Methods for Antenna Rotating Center Based on Elevation Axes Intersection and Spherical Fitting[J]. Geomatics and Information Science of Wuhan University, 2019, 44(10): 1449-1456. doi: 10.13203/j.whugis20180064

基于俯仰轴相交和球拟合的天线旋转中心测量方法

doi: 10.13203/j.whugis20180064
基金项目: 中国电子科技集团公司航天信息应用技术重点实验室开放基金
详细信息
    作者简介:

    李宗春, 教授, 博士生导师, 主要研究方向为精密工程测量。13838092876@139.com

  • 中图分类号: P258

Two Measuring Methods for Antenna Rotating Center Based on Elevation Axes Intersection and Spherical Fitting

Funds: The Open Foundation of CETC Key Laboratory of Aerospace Information Applications
More Information
    Author Bio:

    LI Zongchun, professor, PhD specializes in precise engineering survey. E-mail:13838092876@139.com

  • 摘要: 针对传统天线旋转中心测量时,俯仰轴生成方式不严密对其精度和可靠性产生不利影响的问题,提出了俯仰轴相交和球拟合两种新的旋转中心测量方法。对不同方位下标志点绕俯仰轴的运动轨迹进行圆拟合,得到含有方位信息的俯仰轴,多条俯仰轴求交点作为旋转中心,称之为俯仰轴相交法,较传统方法可靠性更高。对标志点的运动轨迹进行球拟合,取所有球心的均值作为旋转中心,称之为球拟合法,避免了圆拟合及求交点的过程,方法更为严密。以3台φ3 m上行阵天线为实验对象,应用全站仪和工业摄影测量技术获取标志点在规划姿态下的坐标,采用两种新方法处理观测数据,均得到了内符合精度达亚毫米级的旋转中心。
  • 图  1  传统天线旋转中心生成示意图

    Figure  1.  Sketch Map of Generating Traditional Rotation Center

    图  2  标志点绕俯仰轴运动的轨迹

    Figure  2.  Trajectory of a Target Point Around the Elevation Axis

    图  3  俯仰轴相交法得到的旋转中心

    Figure  3.  Rotation Center Obtained by Elevation Axes Intersection Method

    图  4  半球面模型

    Figure  4.  Hemispherical surface model

    图  5  球面拟合

    Figure  5.  Spherical Surface Fitting

    图  6  工程控制网及天线的分布

    Figure  6.  Engineering Control Network and Antennas Distribution

    图  7  标志及工装分布图

    Figure  7.  Distribution of Markers and Fixtures

    图  8  球度分布图

    Figure  8.  Sphericity Distribution

    图  9  公共点分布图

    Figure  9.  Distribution of Common Points

    图  10  旋转中心到俯仰轴的距离

    Figure  10.  Distances from Rotation Center to Elevation Axes

    表  1  GNSS观测数据处理结果

    Table  1.   Processing Results of GNSS Observation Data

    点名 大地坐标 中误差/mm 大地直角坐标/m
    J001 B *°16′18.259 7″ ±0.4 X *7 024.7 281
    L *°35′53.355 5″ ±0.3 Y *8 845.355 3
    H 66.411 4 m ±0.8 Z *9 209.548 0
    J002 B *°16′17.058 7″ ±0.3 X *7 012.947 2
    L *°35′52.390 5″ ±0.2 Y *8 875.978 8
    H 66.416 7 m ±0.6 Z *9 180.479 0
    J003 B *°16′17.985 0″ ±0.3 X *6 994.109 6
    L *°35′51.871 3″ ±0.3 Y *8 865.146 8
    H 66.416 6 m ±0.7 Z *9 202.901 3
    J004 B *°16′18.968 5″ ±0.0 X *7 012.004 5
    L *°35′53.034 7″ ±0.0 Y *8 836.294 7
    H 66.414 8 m ±0.0 Z *9 226.708 5
    注:*为省略数据
    下载: 导出CSV

    表  2  J004至其他各点的大地方位角及斜距

    Table  2.   Geodetic Azimuth and Slope Distance from J004 to Other Points

    起点 终点 大地方位角 斜距/m
    J004 J001 160° 21′ 49.119″ 23.205
    J002 194°53′ 29.198″ 60.933
    J003 223°00′3.620″ 41.466
    下载: 导出CSV

    表  3  控制点坐标及精度

    Table  3.   Coordinates and Precision of Control Points

    点名 X/m Y/m RMS/mm
    J001 -21.855 9 -7.797 8 ±0.10
    J002 -58.887 4 15.661 6 ±0.10
    J003 -30.327 4 28.281 6 ±0.11
    J004 0 0 ±0.00
    下载: 导出CSV

    表  4  工业摄影测量精度

    Table  4.   Precision of Industrial Photogrammetry

    精度 俯仰角
    10° 30° 40° 50° 60°
    RMS/mm ±0.023 ±0.021 ±0.025 ±0.025 ±0.021
    下载: 导出CSV

    表  5  公共点转换精度

    Table  5.   Precision of Coordinate Transformation

    项目 俯仰角
    10° 30° 40° 50° 60°
    公共点个数 5 7 5 6 5
    RMS/mm ±0.106 ±0.122 ±0.134 ±0.095 ±0.148
    下载: 导出CSV

    表  6  旋转中心坐标及精度(俯仰轴相交)

    Table  6.   Coordinates and Precision of Antennas Rotation Center (Elevation Axes Intersection)

    天线 拟合点数 X/mm Y/mm Z/mm RMS/mm
    T1 136 -13 664.21 -11 219.18 970.35 ±0.57
    T2 142 -12 746.45 -1 918.48 994.76 ±0.21
    T3 143 -13 337.59 7 206.31 981.46 ±0.18
    下载: 导出CSV

    表  7  旋转中心坐标及精度(球拟合)

    Table  7.   Coordinates and Precision of Antennas Rotation Center (Spherical Fitting)

    天线 拟合点数 X/mm Y/mm Z/mm RMS/mm
    T1 124 -13 666.87 -11 216.16 973.49 ±0.80
    T2 139 -12 748.05 -1 917.88 994.78 ±0.09
    T3 123 -13 338.09 7 205.46 980.49 ±0.10
    下载: 导出CSV

    表  8  部分姿态下公共点数量

    Table  8.   Number of Common Points

    方位角 俯仰角 T1 T2 T3
    60° 30° 3 6 7
    50° 3 5 5
    70° 3 5 6
    180° 30° 3 7 7
    50° 4 6 4
    70° 3 6 5
    300° 30° 6 7 6
    50° 6 5 5
    70° 5 6 3
    下载: 导出CSV

    表  9  交点到旋转中心的距离/mm

    Table  9.   Distances from Rotation Center to Intersection/mm

    点号 X Y Z 距离
    47号 -12 747.92 -1 918.38 994.30 0.58
    90号 -12 747.55 -1 918.21 994.21 0.72
    63号 -12 743.11 -1 918.74 994.55 5.01
    89号 -12 742.62 -1 918.80 994.34 5.52
    下载: 导出CSV

    表  10  两种新方法获取旋转中心的坐标偏差/mm

    Table  10.   Coordinates Deviation Between Rotation Centers Obtained from Two New Methods/mm

    天线 ΔX ΔY ΔZ 距离
    T1 2.70 -3.27 -2.71 5.03
    T2 1.61 -0.55 0.12 1.71
    T3 0.44 0.87 0.41 1.05
    下载: 导出CSV
  • [1] 李海涛, 李宇华, 匡乃雪.深空探测中的天线组阵技术[J].飞行器测控学报, 2004, 23(4):57-60 http://d.old.wanfangdata.com.cn/Periodical/fxqckxb200404012

    Li Haitao, Li Yuhua, Kuang Naixue. Antenna Array Forming Technology in Deep Space Exploration[J]. Journal of Spacecraft TT&C Technology, 2004, 23(4):57-60 http://d.old.wanfangdata.com.cn/Periodical/fxqckxb200404012
    [2] 于志坚, 李海涛, 李小梅.可扩充深空天线组阵技术研究与试验验证[J].遥测遥控, 2015, 36(6):1-7 doi:  10.3969/j.issn.2095-1000.2015.06.001

    Yu Zhijian, Li Haitao, Li Xiaomei. Research and Experiment Verification on Extendable Deep Space Antenna Arraying Techniques[J]. Journal of Telemetry, Tracking and Command, 2015, 36(6):1-7 doi:  10.3969/j.issn.2095-1000.2015.06.001
    [3] 姚飞, 匡麟玲, 詹亚锋, 等.深空通信天线组阵关键技术及其发展趋势[J].宇航学报, 2010, 31(10):2231-2238 doi:  10.3873/j.issn.1000-1328.2010.10.001

    Yao Fei, Kuang Linling, Zhan Yafeng, et al. Key Techniques and Development Trend of Antenna Arraying for Deep Space Communication[J]. Journal of Astronautics, 2010, 31(10):2231-2238 doi:  10.3873/j.issn.1000-1328.2010.10.001
    [4] 白岩夫.高速率深空通信中天线组阵关键技术的研究[D].西安: 西安电子科技大学, 2014 http://cdmd.cnki.com.cn/Article/CDMD-10701-1014324883.htm

    Bai Yanfu. A Study of Key Techniques for Antenna Arraying for Deep Space Communication[D]. Xi'an: Xidian University, 2014 http://cdmd.cnki.com.cn/Article/CDMD-10701-1014324883.htm
    [5] 周三文, 卢满宏, 黄建国.天线组阵全频谱合成效率分析[J].遥测遥控, 2009, 30(2):46-52 doi:  10.3969/j.issn.2095-1000.2009.02.009

    Zhou Sanwen, Lu Manhong, Huang Jianguo. Analysis on Combining Efficiency of Full -Spectrum Combining[J]. Journal of Telemetry, Tracking and Command, 2009, 30(2):46-52 doi:  10.3969/j.issn.2095-1000.2009.02.009
    [6] 何国龙, 李国民.上行天线组阵载波相位并行标校方法研究[J].飞行器测控学报, 2012, 31(1):19-22 http://d.old.wanfangdata.com.cn/Periodical/fxqckxb201201005

    He Guolong, Li Guomin. A Parallel Carrier Phase Calibration Method for Uplink Antenna Arrays[J]. Journal of Spacecraft TT&C Technology, 2012, 31(1):19-22 http://d.old.wanfangdata.com.cn/Periodical/fxqckxb201201005
    [7] 陈俊勇.国际地球参考框架2000(ITRF2000)的定义及其参数[J].武汉大学学报·信息科学版, 2005, 30(9):753-756 http://ch.whu.edu.cn/CN/abstract/abstract2264.shtml

    Chen Junyong. On the Definition and Adopted Parameters of International Terrestrial Reference Frame 2000[J]. Geomatics and Information Science of Wuhan University, 2005, 30(9):753-756 http://ch.whu.edu.cn/CN/abstract/abstract2264.shtml
    [8] Schuh H, Bbehrend D. VLBI:A Fascinating Technique for Geodesy and Astrometry[J]. Journal of Geodynamics, 2012, 61:68-80 doi:  10.1016/j.jog.2012.07.007
    [9] 董文亮, 孙付平, 朱新慧, 等.统一全球地壳垂直运动参考基准的方法研究[J].武汉大学学报·信息科学版, 2014, 39(3):358-361 http://ch.whu.edu.cn/CN/abstract/abstract2915.shtml

    Dong Wenliang, Sun Fuping, Zhu Xinhui, et al. A Method of Unifying the Global Crustal Vertical Motion Reference Datum[J]. Geomatics and Information Science of Wuhan University, 2014, 39(3):358-361 http://ch.whu.edu.cn/CN/abstract/abstract2915.shtml
    [10] Bergstrand S, Schmid R. Activities of the IERS Working Group on Site Survey and Co-location[C]. International VLBI Service for Geodesy and Astrometry 2016 General Meeting, South Africa, 2016
    [11] 孙中苗, 范昊鹏. VLBI全球观测系统(VGOS)研究进展[J].测绘学报, 2017, 46(10):1346-1353 doi:  10.11947/j.AGCS.2017.20170326

    Sun Zhongmiao, Fan Haopeng. Research Progress of VLBI Global Observing System (VGOS)[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10):1346-1353 doi:  10.11947/j.AGCS.2017.20170326
    [12] Altamimi Z, Rebischung P, Métivier L, et al. ITRF2014:A New Release of the International Terrestrial Reference Frame Modeling Nonlinear Station Motions[J]. Geophys Res Solid Earth, 2016, 121(8):6109-6131 doi:  10.1002/2016JB013098
    [13] 马下平, 沈云中, 王解先, 等.卫星激光测距与甚长基线干涉测量的天线旋转中心的直接解法[J].测绘学报, 2014, 43(3):257-262 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=chxb201403006

    Ma Xiaping, Shen Yunzhong, Wang Jiexian, et al. Direct Solution of SLR and VLBI Antenna Rotation Center[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(3):257-262 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=chxb201403006
    [14] Dawson J, Sarti P, Johnston G M, et al. Indirect Approach to Invariant Point Determination for SLR and VLBI Systems:An Assessment[J]. Journal of Geodesy, 2007, 81(6):433-441
    [15] LöSler M. New Mathematical Model for Reference Point Determination of an Azimuth Elevation Type Radio Telescope[J]. Journal of Surveying Engineering, 2009, 135(4):131-135 doi:  10.1061/(ASCE)SU.1943-5428.0000010
    [16] 李金岭, 张津维, 郭丽.甚长基线干涉测量天线参考点和轴线偏差监测方法探讨[J].测绘科学, 2014, 39(11):10-14 http://d.old.wanfangdata.com.cn/Periodical/chkx201411003

    Li Jinling, Zhang Jinwei, Guo Li. Discussion on Monitoring Method of Reference Point and Axis Offset of VLBI Antenna[J]. Science of Surveying and Mapping, 2014, 39(11):10-14 http://d.old.wanfangdata.com.cn/Periodical/chkx201411003
    [17] Gil K H, Sang H J, Yun H S, et al. A Comparative Analysis of 3D Circle Fitting Algorithms for Determination of VLBI Antenna Reference Point[J]. Journal of the Korean Society of Surveying Geodesy Photogrammetry & Cartography, 2015, 33(4):231-244 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=JAKO201530848469750
    [18] 沈云中, 陈廷武.上海天文台并址站的空间归心测量[J].同济大学学报(自然科学版), 2006, 34(2):217-220 doi:  10.3321/j.issn:0253-374X.2006.02.015

    Shen Yunzhong, Chen Tingwu. Determination of Space Coordinate Differences of Co-location Sites in Shanghai Observatory[J]. Journal of Tongji University (Natural Science) 2006, 34(2):217-220 doi:  10.3321/j.issn:0253-374X.2006.02.015
    [19] 李金岭, 乔书波, 刘鹂, 等. 2008年佘山25 m射电天线归心测量[J].武汉大学学报·信息科学版, 2010, 35(12):1387-1391 http://www.cnki.com.cn/Article/CJFDTotal-WHCH201012002.htm

    Li Jinling, Qiao Shubo, Liu Li, et al. Site Survey at Sheshan 25 m Radio Telescope in 2008[J]. Geomatics and Information Science of Wuhan University, 2010, 35(12):1387-1391 http://www.cnki.com.cn/Article/CJFDTotal-WHCH201012002.htm
    [20] 邓勇, 李宗春, 张万才, 等. 65 m射电望远镜方位轴和俯仰轴正交度测量[J].无线电工程, 2016, 46(7):67-70 doi:  10.3969/j.issn.1003-3106.2016.07.18

    Deng Yong, Li Zongchun, Zhang Wancai, et al. Measurement of 65 m Antenna Azimuth Axis and Elevation Axis Orthogonality[J]. Radio Engineering, 2016, 46(7):67-70 doi:  10.3969/j.issn.1003-3106.2016.07.18
    [21] 田镇, 杨志强, 石震, 等.国家授时中心昊平站40 m口径射电天线相位中心参考点坐标的精密测定[J].大地测量与地球动力学, 2016, 36(10):897-901 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz201610012

    Tian Zhen, Yang Zhiqiang, Shi Zhen, et al. Precise Measurement of 40 m Caliber Radio Telescope Phase Center's Reference Point Coordinates at the Haoping Station of National Time Service Center[J]. Journal of Geodesy and Geodynamic, 2016, 36(10):897-901 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz201610012
    [22] 李金岭, 张津维.利用GPS测量监测VLBI天线参考点的仿真分析[J].武汉大学学报·信息科学版, 2013, 38(12):1387-1391 http://ch.whu.edu.cn/CN/abstract/abstract2816.shtml

    Li Jinling, Zhang Jinwei. Simulation Analysis of Monitoring of VLBI Antennas Reference Point Via GPS Observation[J]. Geomatics and Information Science of Wuhan University, 2013, 38(12):1387-1391 http://ch.whu.edu.cn/CN/abstract/abstract2816.shtml
  • [1] 张柔, 胡志刚, 陶钧, 王晨, 赵齐乐.  顾及不同天线相位中心改正模型的北斗空间信号精度评估方法 . 武汉大学学报 ● 信息科学版, 2019, 44(6): 806-813. doi: 10.13203/j.whugis20180388
    [2] 张志斌, 王广利, 宋世泽, 张阿丽, 王宏.  归心测量中轮轨式VLBI天线的水平钢轨形变建模 . 武汉大学学报 ● 信息科学版, 2019, 44(10): 1498-1504. doi: 10.13203/j.whugis20180050
    [3] 梅文胜, 胡帅朋, 李谋思, 祁洪宇, 徐芳.  基于普通数码相机的旋转全景摄影测量方法 . 武汉大学学报 ● 信息科学版, 2017, 42(2): 243-249. doi: 10.13203/j.whugis20140994
    [4] 张鑫, 陈华明, 黄仰博, 张国柱, 欧钢.  利用旋转天线载波相位双差进行欺骗干扰检测技术 . 武汉大学学报 ● 信息科学版, 2016, 41(4): 529-534. doi: 10.13203/j.whugis20140360
    [5] 蒋振伟, 袁运斌, 欧吉坤, 王海涛.  基于站间单差相位绕转数据估计测站天线旋转速率的方法研究 . 武汉大学学报 ● 信息科学版, 2015, 40(10): 1334-1338. doi: 10.13203/j.whugis20130729
    [6] 马洋, 欧吉坤, 袁运斌, 霍星亮, 丁文武.  导航卫星天线相位中心变化估计及对LEO精密定轨影响 . 武汉大学学报 ● 信息科学版, 2015, 40(7): 894-900. doi: 10.13203/j.whugis20130626
    [7] 李金岭, 张津维.  利用GPS测量监测VLBI天线参考点的仿真分析 . 武汉大学学报 ● 信息科学版, 2013, 38(12): 1387-1391.
    [8] 李永江, 艾力·玉苏甫, 张正禄, 艾力·伊沙木丁.  天线轨道变形精密测量与指向偏差模型研究 . 武汉大学学报 ● 信息科学版, 2013, 38(2): 176-180.
    [9] 张小红, 李盼, 李星星, 郭斐.  天线相位中心改正模型对PPP参数估计的影响 . 武汉大学学报 ● 信息科学版, 2011, 36(12): 1470-1473.
    [10] 柯涛, 张祖勋, 张剑清.  旋转多基线数字近景摄影测量 . 武汉大学学报 ● 信息科学版, 2009, 34(1): 44-47.
    [11] 郭际明, 史俊波, 汪伟.  天线相位中心偏移和变化对高精度GPS数据处理的影响 . 武汉大学学报 ● 信息科学版, 2007, 32(12): 1143-1146.
    [12] 王解先.  工业测量中一种二次曲面的拟合方法 . 武汉大学学报 ● 信息科学版, 2007, 32(1): 47-50.
    [13] 郭金运, 徐泮林, 曲国庆.  GPS接收机天线相位中心偏差的三维检定研究 . 武汉大学学报 ● 信息科学版, 2003, 28(4): 448-451.
    [14] 孙海燕, 李瑞明, 闫利.  机载三线阵CCD摄影中心轨迹及姿态拟合的半参数法 . 武汉大学学报 ● 信息科学版, 2003, 28(6): 706-709.
    [15] 冯文灏.  近景摄影测量的控制 . 武汉大学学报 ● 信息科学版, 2000, 25(5): 453-458.
    [16] 廖祥春, 冯文灏.  圆形标志及其椭圆构像间中心偏差的确定 . 武汉大学学报 ● 信息科学版, 1999, 24(3): 235-239.
    [17] 熊春宝, 叶声华.  经纬仪工业测量系统的模型研究 . 武汉大学学报 ● 信息科学版, 1998, 23(3): 233-237.
    [18] 刘凤德, 郏卫萍, 邱峰, 刘先林.  JX-4A数字摄影测量系统 . 武汉大学学报 ● 信息科学版, 1996, 21(3): 269-272,276.
    [19] 施品浩.  GPS接收机天线装置定向标志线及平均相位中心参数的检测 . 武汉大学学报 ● 信息科学版, 1990, 15(1): 28-38.
    [20] 李德仁, 程家瑜, 许妙忠, 王树根.  SPOT影象的解析摄影测量处理 . 武汉大学学报 ● 信息科学版, 1988, 13(4): 28-36.
  • 加载中
图(10) / 表(10)
计量
  • 文章访问数:  1061
  • HTML全文浏览量:  104
  • PDF下载量:  112
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-08-12
  • 刊出日期:  2019-10-05

基于俯仰轴相交和球拟合的天线旋转中心测量方法

doi: 10.13203/j.whugis20180064
    基金项目:  中国电子科技集团公司航天信息应用技术重点实验室开放基金
    作者简介:

    李宗春, 教授, 博士生导师, 主要研究方向为精密工程测量。13838092876@139.com

  • 中图分类号: P258

摘要: 针对传统天线旋转中心测量时,俯仰轴生成方式不严密对其精度和可靠性产生不利影响的问题,提出了俯仰轴相交和球拟合两种新的旋转中心测量方法。对不同方位下标志点绕俯仰轴的运动轨迹进行圆拟合,得到含有方位信息的俯仰轴,多条俯仰轴求交点作为旋转中心,称之为俯仰轴相交法,较传统方法可靠性更高。对标志点的运动轨迹进行球拟合,取所有球心的均值作为旋转中心,称之为球拟合法,避免了圆拟合及求交点的过程,方法更为严密。以3台φ3 m上行阵天线为实验对象,应用全站仪和工业摄影测量技术获取标志点在规划姿态下的坐标,采用两种新方法处理观测数据,均得到了内符合精度达亚毫米级的旋转中心。

English Abstract

李宗春, 路志勇, 郭迎钢, 张冠宇, 何华, 冯其强, 陈少卿, 王俊威. 基于俯仰轴相交和球拟合的天线旋转中心测量方法[J]. 武汉大学学报 ● 信息科学版, 2019, 44(10): 1449-1456. doi: 10.13203/j.whugis20180064
引用本文: 李宗春, 路志勇, 郭迎钢, 张冠宇, 何华, 冯其强, 陈少卿, 王俊威. 基于俯仰轴相交和球拟合的天线旋转中心测量方法[J]. 武汉大学学报 ● 信息科学版, 2019, 44(10): 1449-1456. doi: 10.13203/j.whugis20180064
LI Zongchun, LU Zhiyong, GUO Yinggang, ZHANG Guanyu, HE Hua, FENG Qiqiang, CHEN Shaoqing, WANG Junwei. Two Measuring Methods for Antenna Rotating Center Based on Elevation Axes Intersection and Spherical Fitting[J]. Geomatics and Information Science of Wuhan University, 2019, 44(10): 1449-1456. doi: 10.13203/j.whugis20180064
Citation: LI Zongchun, LU Zhiyong, GUO Yinggang, ZHANG Guanyu, HE Hua, FENG Qiqiang, CHEN Shaoqing, WANG Junwei. Two Measuring Methods for Antenna Rotating Center Based on Elevation Axes Intersection and Spherical Fitting[J]. Geomatics and Information Science of Wuhan University, 2019, 44(10): 1449-1456. doi: 10.13203/j.whugis20180064
  • 上行天线组阵技术[1-2]利用多台天线向同一个深空探测器发射上行信号[3],具有提升系统性能、增强可操作性、降低成本和提高系统运行灵活性等明显优势[4-5]。为了实现各发射信号在目标处的同相合成,需要精确测定各天线参考点的相对位置[6]。具有长基线的天线阵可以支持相关的科学应用,如国际地球参考框架(international terrestrial reference frame,ITRF)的建立和维持[7]、深空探测、射电天文学研究等。ITRF的建立和维持需要综合利用甚长基线干涉测量(very long baseline interferometry,VLBI)、卫星激光测距(satellite laser ranging,SLR)、全球导航卫星系统(global navigation satellite system,GNSS)、多普勒轨道学与无线电定位集成卫星(Doppler orbitography and radiopositioning integrated by satellite,DORIS)等空间大地测量技术来实现[8],多种大地测量技术的综合利用要求将相应的观测台站并置[9]。并置台站的本地连接[10]误差是构建毫米级精度地球参考框架的重要影响因素[11-12],这就需要提高天线参考点测量结果的精度和可靠性。为此,本文分析了传统测量方法所存在的缺陷,并就如何获取精密、可靠的天线参考点进行了深入的探讨。

    天线旋转中心[13]是指天线方位轴与俯仰轴的交点,通常位于天线内部,不易直接观测,一般通过先观测天线上的标志点,再采用相关算法间接得到。标志点的布设、观测方法以及数据处理方法直接影响着旋转中心的精度和可靠性。

    国内外许多学者对天线旋转中心测量方法进行了研究。文献[14]分析讨论了在天线方位旋转受限的条件下,轴拟合、圆拟合等的几何条件对旋转中心精度的影响。文献[15]和文献[16]对VLBI天线参考点的自动监测进行了研究。文献[17]通过三维空间圆拟合来获取方位轴和俯仰轴,两轴相交得到VLBI天线的旋转中心。文献[18]采用旋转天线的方法,运用最小二乘曲线拟合得到旋转中心,外符合精度达到1 cm。文献[19]在上海天文台佘山φ25 m射电望远镜归心测量中,对测量结果进行拟合,取俯仰轴与方位轴交点作为旋转中心,外符合精度达到1 mm。文献[20]在上海天文台φ65 m天马望远镜的方位轴和俯仰轴正交度测量中,通过不同的拟合方式求取方位轴和俯仰轴,对比了各拟合方式的特点。文献[21]对国家授时中心昊平站φ40 m天线的旋转中心进行了测定,在2000国家大地坐标系下的点位精度优于8 mm。文献[22]对利用GNSS数据监测VLBI天线进行了仿真分析,期望对旋转中心实现1 mm级的监测。

    以上文献大多应用全站仪测量天线上标志点、拟合方位轴和俯仰轴,取两轴交点为旋转中心,如图 1所示。此方法存在两个缺陷:①标志点数目较少,易受测量误差、结构变形等干扰,测量结果不可靠;②生成俯仰轴的方式不够严密,主要表现为拟合俯仰轴时只在一个方位下进行观测,且标志点轨迹至多能覆盖1/4圆周,不利于圆心的准确拟合,而且受自重变形等因素的影响,用圆拟合标志点运动轨迹会引入模型误差。以上因素导致俯仰轴的测量结果可靠性不高,而方位轴的拟合不受这些因素的影响,这就导致方位轴和俯仰轴精度不匹配,从而影响了旋转中心的精度和可靠性。

    图  1  传统天线旋转中心生成示意图

    Figure 1.  Sketch Map of Generating Traditional Rotation Center

    为此,本文拟构建基于俯仰轴相交和球拟合两种天线旋转中心测量方法,以提高旋转中心的测定精度和可靠性。

    • 已有的天线旋转中心获取方法都只是利用某一方位下的俯仰轴与方位轴求交点,该交点的精度和可靠性不高。故可考虑在多个方位(3个及以上,在全圆内等间隔分布)下获取俯仰轴,以这些俯仰轴的交点作为旋转中心。

    • 理想情况下,当方位角固定时,天线上标志点绕俯仰轴运动的轨迹可以看作是1/4圆,如图 2所示。

      图  2  标志点绕俯仰轴运动的轨迹

      Figure 2.  Trajectory of a Target Point Around the Elevation Axis

    • 单个标志点对应的多方位下的俯仰轴相交可以得到一个旋转中心,如图 3所示。

      图  3  俯仰轴相交法得到的旋转中心

      Figure 3.  Rotation Center Obtained by Elevation Axes Intersection Method

      基于此思路,本文提出了俯仰轴相交法:连续、均匀地选取天线姿态,测量天线上的标志点,对每个标志点在既定方位下绕俯仰轴的运动轨迹进行圆拟合以得到俯仰轴,取各俯仰轴的交点作为旋转中心,并取多点所对应的旋转中心的均值作为最终结果,来提高其精度和可靠性。具体步骤如下:

      1) 获取观测数据。建立工程控制网,定义工程坐标系,在天线面板上布设标志点,利用工业摄影测量系统和全站仪测量系统进行联合测量,获取不同姿态下的标志点在工程坐标系下的坐标。

      2) 拟合圆,取过圆心的圆平面法线作为俯仰轴。空间圆的表达方式为:

      $$ \left\{\begin{array}{l}{A X+B Y+C Z+D=0} \\ {\left(X-X_{0}\right)^{2}+\left(Y-Y_{0}\right)^{2}+\left(Z-Z_{0}\right)^{2}=R^{2}} \\ {A X_{0}+B Y_{0}+C Z_{0}+D=0}\end{array}\right. $$ (1)

      式中,(X, Y, Z)为目标点坐标;(X0, Y0, Z0)为圆心O的坐标;ABCD为空间圆所在平面的方程参数;R为空间圆半径。

      俯仰轴l过圆心O且垂直于空间圆所在平面,方向为$u=\left(u_{x}, u_{y}, u_{z}\right) $,则l的方程可表示为:

      $$ \frac{x-X_{0}}{u_{x}}=\frac{y-Y_{0}}{u_{y}}=\frac{z-Z_{0}}{u_{z}}=a $$ (2)

      式中,a为比例参数,其绝对值等于轴线上点$ (x, y, z)$到圆心的距离。

      3) 俯仰轴相交得旋转中心。根据式(2)列误差方程:

      $$ \boldsymbol{V}=\boldsymbol{A} \boldsymbol{X}-\boldsymbol{L} $$ (3)

      式中,$\boldsymbol{X}=\left[\begin{array}{lllll}{x_{r}} & {y_{r}} & {z_{r}} & {a_{1}} & {a_{2} \cdots a_{n}}\end{array}\right]^{\mathrm{T}} $

      $$ \boldsymbol{L}=\left[\begin{array}{ll}{O_{1}} & {O_{2} \cdots O_{n}}\end{array}\right]^{\mathrm{T}} $$
      $$ \boldsymbol{A}=\left[\begin{array}{ccccc}{\boldsymbol{I}_{3 \times 3}} & {-\boldsymbol{u}_{1}^{\mathrm{T}}} & {\boldsymbol{0}_{3 \times 1}} & {\cdots} & {\boldsymbol{0}_{3 \times 1}} \\ {\boldsymbol{I}_{3 \times 3}} & {\boldsymbol{0}_{3 \times 1}} & {-\boldsymbol{u}_{2}^{\mathrm{T}}} & {\cdots} & {\boldsymbol{0}_{3 \times 1}} \\ {\vdots} & {\vdots} & {\vdots} & {\cdots} & {\vdots} \\ {\boldsymbol{I}_{3 \times 3}} & {\boldsymbol{0}_{3 \times 1}} & {\boldsymbol{0}_{3 \times 1}} & {\cdots} & {-\boldsymbol{u}_{n}^{\mathrm{T}}}\end{array}\right] $$

      依据最小二乘原理计算俯仰轴交点:

      $$ \boldsymbol{X}=-\left(\boldsymbol{A}^{\mathrm{T}} \boldsymbol{A}\right) \cdot\left(\boldsymbol{A}^{\mathrm{T}} \boldsymbol{L}\right) $$ (4)

      参数X的前3个分量$ \left(x_{r}, y_{r}, z_{r}\right)$即为旋转中心坐标。

      4) 取所有标志点所对应旋转中心的平均值作为最终的旋转中心。

      $$ \bar{X}=\frac{1}{n} \sum\limits_{j=1}^{n} x_{r j}, \bar{Y}=\frac{1}{n} \sum\limits_{j=1}^{n} y_{r j}, \bar{Z}=\frac{1}{n} \sum\limits_{j=1}^{n} z_{r j} $$ (5)

      式中,$\left(x_{r j}, y_{r j}, z_{r j}\right) $表示第j个目标点对应的旋转中心坐标;$ (\bar{X}, \bar{Y}, \bar{Z})$表示最终的旋转中心坐标。

      5) 精度评定。取n个旋转中心均值的中误差来评价其内符合精度,即

      $$ \begin{array}{l} m = \\ \pm \sqrt {\frac{{\sum\limits_{j = 1}^n {\left( {{{\left( {{x_n} - \bar X} \right)}^2} + {{\left( {{y_n} - \bar Y} \right)}^2} + {{\left( {{z_n} - \bar Z} \right)}^2}} \right)} }}{{n(n - 1)}}} \end{array} $$ (6)
    • 理想情况下,天线上标志点随方位、俯仰角变化的运动轨迹可以看作是半球面模型。根据该运动模型,只要获取不同姿态下标志点的三维坐标,通过球拟合获得球心,此球心即为旋转中心。此模型避免了圆拟合和求轴线交点的过程,相比于俯仰轴相交法更为严密,能进一步减小模型误差。

    • M为球面上任意一点,球心为$ {M_0}\left( {{x_0}, {y_0}, {z_0}} \right)$、半径为r的半球面如图 4所示。

      图  4  半球面模型

      Figure 4.  Hemispherical surface model

    • 基于半球面模型提出球拟合法,方法为:连续、均匀地选取天线姿态,测量其上的标志点,通过球拟合得到球心,取所有球心的均值作为旋转中心。该方法可将加工装配、温度及风荷载等因素的影响随机化,提高测量结果的精度和可靠性。具体步骤如下:

      1) 获取观测数据。与§1.2中步骤1)相同。

      2) 拟合球。球面方程可表示为:

      $$ f=\left(X-x_{0}\right)^{2}+\left(Y-y_{0}\right)^{2}+\left(Z-z_{0}\right)^{2}-r^{2} $$ (7)

      式中,(X, Y, Z)为目标点的坐标;$ \left( {{x_0}, {y_0}, {z_0}} \right)$为球心坐标;r为球的半径。

      图 5所示,将半径、球心坐标看作未知参数,以点到拟合球面的距离平方和最小为准则,对球面方程线性化后依据最小二乘原理求解球心坐标。

      图  5  球面拟合

      Figure 5.  Spherical Surface Fitting

      3) 计算球度,评价球拟合的质量。将第j个标志点的运动轨迹拟合为球,计算其球度为:

      $$ {\lambda _j} = \sqrt {\frac{{\sum\limits_{i = 1}^k {{{\left( {\left\| {{P_{ij}} - {M_{0j}}} \right\| - {r_j}} \right)}^2}} }}{k}} $$ (8)

      式中,Pij表示第j个点的第i个运动位置;M0j表示第j个点的拟合球心;k为第j个点的运动位置数。

      4) 取所有球心的均值作为天线旋转中心,公式同式(5)。

      5) 精度评定。用球拟合所得球心均值的中误差来评价旋转中心的内符合精度,按照式(6)进行计算。

    • 以3台φ3 m天线组成的上行阵为实验对象,对其旋转中心进行了测量。工程控制网及天线的分布如图 6所示,其中T1、T2、T3表示3台天线,J001~J004为工程控制网的4个控制点。

      图  6  工程控制网及天线的分布

      Figure 6.  Engineering Control Network and Antennas Distribution

    • 为了描述天线与探测器的关系,需要通过GNSS观测引入地心大地坐标系;为了满足高精度的三维坐标测量需求,通过全站仪观测构建三维控制网;为了获取可靠的测量信息,通过工业摄影测量技术获取天线面板上大量标志点的坐标。

    • 用3台Leica GPS System 500按照GNSS控制测量中的C级测量技术要求进行施测,以北京房山、上海、拉萨3个国际GPS服务机构基准站的GNSS观测数据及其大地坐标和运动速度计算J002和J004点的大地坐标,以J004点高精度计算结果为控制基准,采用Leica GEO office 8.4软件进行局域网数据处理,得到控制点坐标成果如表 1所示,J004至其他各点的大地方位角及斜距如表 2所示。

      表 1  GNSS观测数据处理结果

      Table 1.  Processing Results of GNSS Observation Data

      点名 大地坐标 中误差/mm 大地直角坐标/m
      J001 B *°16′18.259 7″ ±0.4 X *7 024.7 281
      L *°35′53.355 5″ ±0.3 Y *8 845.355 3
      H 66.411 4 m ±0.8 Z *9 209.548 0
      J002 B *°16′17.058 7″ ±0.3 X *7 012.947 2
      L *°35′52.390 5″ ±0.2 Y *8 875.978 8
      H 66.416 7 m ±0.6 Z *9 180.479 0
      J003 B *°16′17.985 0″ ±0.3 X *6 994.109 6
      L *°35′51.871 3″ ±0.3 Y *8 865.146 8
      H 66.416 6 m ±0.7 Z *9 202.901 3
      J004 B *°16′18.968 5″ ±0.0 X *7 012.004 5
      L *°35′53.034 7″ ±0.0 Y *8 836.294 7
      H 66.414 8 m ±0.0 Z *9 226.708 5
      注:*为省略数据

      表 2  J004至其他各点的大地方位角及斜距

      Table 2.  Geodetic Azimuth and Slope Distance from J004 to Other Points

      起点 终点 大地方位角 斜距/m
      J004 J001 160° 21′ 49.119″ 23.205
      J002 194°53′ 29.198″ 60.933
      J003 223°00′3.620″ 41.466

      以最长边J004-J002的大地方位角为方位基准,结合表 1可知,GNSS观测数据的处理结果精度较高,满足工程需求。

    • 用Leica MS50全站仪按照三等控制网测量技术要求施测,忽略垂线偏差的影响。以J004为原点,以J004点所在的大地经线北方向为X轴,水平向西方向为Y轴,J004点的法线向上方向为Z轴,构建工程坐标系,得到控制点的平面坐标如表 3所示。

      表 3  控制点坐标及精度

      Table 3.  Coordinates and Precision of Control Points

      点名 X/m Y/m RMS/mm
      J001 -21.855 9 -7.797 8 ±0.10
      J002 -58.887 4 15.661 6 ±0.10
      J003 -30.327 4 28.281 6 ±0.11
      J004 0 0 ±0.00

      表 3可知,控制点的平面坐标精度很高,都在0.1 mm左右。

    • 工业摄影测量需要回光反射标志的配合,在天线面板上均匀地布设了摄影测量标志。此工程中,3台φ3 m天线结构一致,均由8块相同的扇形面板组成,各面板标志布设方案相同,可按照图 7(a)所示的方案实施。为了实现工业摄影测量坐标系与工程坐标系的转换,需要一定数量的公共点,其分布如图 7(b)所示。公共点处安装有特制的靶座,既可安放球棱镜配合全站仪测量,又可安放摄影测量工装,且通过高精度加工保证了摄影测量工装与球棱镜的同心。标志和工装的实际分布如图 7(c)所示。

      图  7  标志及工装分布图

      Figure 7.  Distribution of Markers and Fixtures

      分别在方位为60°、180°和300°的位置,俯仰从10°~80°(间隔10°)对各天线进行摄影测量,得到标志点在各工业摄影测量坐标系下的坐标。通过全站仪观测,获取公共点在工程坐标系下的坐标,然后将所有标志点的坐标统一转换至工程坐标系下。

      以T2为例,在180°方位下不同俯仰角所对应的工业摄影测量精度如表 4所示。

      表 4  工业摄影测量精度

      Table 4.  Precision of Industrial Photogrammetry

      精度 俯仰角
      10° 30° 40° 50° 60°
      RMS/mm ±0.023 ±0.021 ±0.025 ±0.025 ±0.021

      同样以T2为例,在180°方位下不同俯仰角所对应的公共点转换精度如表 5所示。

      表 5  公共点转换精度

      Table 5.  Precision of Coordinate Transformation

      项目 俯仰角
      10° 30° 40° 50° 60°
      公共点个数 5 7 5 6 5
      RMS/mm ±0.106 ±0.122 ±0.134 ±0.095 ±0.148

      表 4表 5可知,工业摄影测量精度高,公共点转换精度较高,符合预期精度要求。

    • 用俯仰轴相交法得到3台天线的旋转中心在工程坐标系下的坐标,并根据式(6)计算其内符合精度,结果如表 6所示。

      表 6  旋转中心坐标及精度(俯仰轴相交)

      Table 6.  Coordinates and Precision of Antennas Rotation Center (Elevation Axes Intersection)

      天线 拟合点数 X/mm Y/mm Z/mm RMS/mm
      T1 136 -13 664.21 -11 219.18 970.35 ±0.57
      T2 142 -12 746.45 -1 918.48 994.76 ±0.21
      T3 143 -13 337.59 7 206.31 981.46 ±0.18

      观察表 6可知,俯仰轴相交法所得旋转中心的内符合精度达到了亚毫米级。

    • 以T2为例,通过球拟合求解旋转中心,按照式(8)计算球度,分布如图 8所示。

      图  8  球度分布图

      Figure 8.  Sphericity Distribution

      通过图 8可以看出,标志点对应的球度均在0.5 mm以内,表明利用球拟合计算旋转中心时,各标志点计算出的球心差异很小,能够据此获得高精度的旋转中心。

      计算得到3台天线的旋转中心在工程坐标系下的坐标,并依据式(6)计算其内符合精度,结果如表 7所示。

      表 7  旋转中心坐标及精度(球拟合)

      Table 7.  Coordinates and Precision of Antennas Rotation Center (Spherical Fitting)

      天线 拟合点数 X/mm Y/mm Z/mm RMS/mm
      T1 124 -13 666.87 -11 216.16 973.49 ±0.80
      T2 139 -12 748.05 -1 917.88 994.78 ±0.09
      T3 123 -13 338.09 7 205.46 980.49 ±0.10

      观察表 7可知,球拟合法所得旋转中心的内符合精度达到了亚毫米级。

    • 1) 两种新方法的主要影响因素。

      对比表 6表 7可以看出,T1的精度最差,T2和T3的精度相近,分析其原因,是由于观测过程中受到公共点的数量及分布等因素的影响。部分姿态下公共点数量如表 8所示,公共点分布如图 9所示。

      表 8  部分姿态下公共点数量

      Table 8.  Number of Common Points

      方位角 俯仰角 T1 T2 T3
      60° 30° 3 6 7
      50° 3 5 5
      70° 3 5 6
      180° 30° 3 7 7
      50° 4 6 4
      70° 3 6 5
      300° 30° 6 7 6
      50° 6 5 5
      70° 5 6 3

      图  9  公共点分布图

      Figure 9.  Distribution of Common Points

      结合表 8图 9进行分析,在观测T1的过程中,公共点较少且分布不均匀,导致公共点转换的过程中损失了精度,影响了旋转中心的精度;在观测T2、T3的过程中,公共点数量较多且分布均匀,保证了最终旋转中心的精度。这也表明在使用本文提出的两种新方法时,公共点数量及分布对测量结果有较大的影响,在实施测量过程中应给予重点关注。

      2) 新方法的优势。

      以T2为例,根据标志点绕俯仰轴旋转的运动轨迹拟合得到多条俯仰轴,以球拟合法获取的旋转中心为准,计算各俯仰轴到旋转中心的距离,如图 10所示。

      图  10  旋转中心到俯仰轴的距离

      Figure 10.  Distances from Rotation Center to Elevation Axes

      观察图 10可知,旋转中心到不同标志点所对应俯仰轴的距离从0.2~8.8 mm不等,差异明显。根据图 10中T2旋转中心到俯仰轴距离的变化情况,选择距离较小的47号点、90号点及距离较大的63号点、89号点,作进一步的分析。此4点按俯仰轴相交法计算交点,并计算到旋转中心的距离,如表 9所示。

      表 9  交点到旋转中心的距离/mm

      Table 9.  Distances from Rotation Center to Intersection/mm

      点号 X Y Z 距离
      47号 -12 747.92 -1 918.38 994.30 0.58
      90号 -12 747.55 -1 918.21 994.21 0.72
      63号 -12 743.11 -1 918.74 994.55 5.01
      89号 -12 742.62 -1 918.80 994.34 5.52

      表 9可知,在按照传统方法获取旋转中心时,如果标志点精度差或者稳定性低,会对旋转中心产生不利的影响。相比之下,本文提出的两种新方法均建立在大量观测数据基础之上,能够剔除不稳定或精度差的标志点。

      3) 两种新方法的对比。

      俯仰轴相交法和球拟合法获得的旋转中心虽然内符合精度都达到了亚毫米级,但二者存在一定的偏差,如表 10所示。

      表 10  两种新方法获取旋转中心的坐标偏差/mm

      Table 10.  Coordinates Deviation Between Rotation Centers Obtained from Two New Methods/mm

      天线 ΔX ΔY ΔZ 距离
      T1 2.70 -3.27 -2.71 5.03
      T2 1.61 -0.55 0.12 1.71
      T3 0.44 0.87 0.41 1.05

      表 10可以看出,T1天线对应的坐标偏差最大,达到了5 mm,与公共点的数量较少和分布情况不佳有关。在公共点数量充足、网形结构好的情况下,T2、T3旋转中心坐标偏差不超过2 mm,结果良好。至于哪一种方法的精度更高,有待于在具备实验条件的情况下通过外符合精度来进行评价。

    • 本文在分析传统天线旋转中心测量方法缺陷的基础上,提出了俯仰轴相交和球拟合两种旋转中心测量方法。这两种方法综合应用空间大地测量技术和精密工程测量技术,顾及天线的整个工作场景,连续、均匀地选取天线姿态,利用工业摄影测量技术获取天线面板上大量标志点的观测值,分别采用含有方位信息的俯仰轴相交或球拟合的处理方式,避免了方位轴、俯仰轴的拟合精度不匹配对旋转中心精度的影响,甚至避开了轴线拟合法而采用更为合理的球拟合法,得到了精确、可靠的旋转中心。以3台ϕ3 m天线组成的上行阵为实验对象对本文方法进行了验证,结果表明两种新方法获取的天线旋转中心内符合精度均能达到亚毫米级。该测量过程为更高性能天线组阵的标校和并置站的本地连接测量积累了经验。

      由于实验条件的限制,未能实现旋转中心外符合精度评定。可考虑在并置台站的本地连接测量中采用本文方法,将计算结果与空间技术解算的基线值进行对比,以评定其外符合精度。

参考文献 (22)

目录

    /

    返回文章
    返回