一种确定性星座对地覆盖计算方法

陈晓宇,戴光明,2,王茂才,2,宋志明,2

(1.中国地质大学 计算机学院，武汉 430074; 2.智能地学信息处理湖北省重点实验室(中国地质大学)，武汉 430074)

一种确定性星座对地覆盖计算方法

陈晓宇1,戴光明1,2,王茂才1,2,宋志明1,2

(1.中国地质大学 计算机学院，武汉 430074; 2.智能地学信息处理湖北省重点实验室(中国地质大学)，武汉 430074)

为提高星地覆盖计算精度和效率，研究了星座对地覆盖问题.首先，分析了星座对地覆盖几何特性，着重研究了基于Delaunay三角网和Voronoi图的星座空间几何构形划分方法;其次，通过定义卫星的所属覆盖区域，将星座对地覆盖问题简化为单星覆盖问题.并在此基础上提出了星座对任意类型地面目标覆盖能力的精确、快速计算方法;最后，比较不同覆盖目标下，所提出的方法和经典的网格点划分方法的计算精度和计算效率.实验结果表明，所提出方法具有正确性和高效性，且该方法可以有效应用于卫星星座优化设计和性能分析中.

卫星星座；覆盖计算；性能分析；Delaunay三角网；Voronoi图

随着卫星技术的发展，由多颗星组成的星座系统在多个领域都得到了重要的应用[1-2].不论是遥感卫星对地成像，还是通信卫星与地面站进行数据交换都涉及到卫星对地覆盖分析.并且大多数的星座优化设计和分析都是基于卫星对地球表面的可视性计算实现的[3].

目前，卫星星座对地覆盖问题最常用的方法是网格点法.Morrison[4]在1973年就提出将区域按照某种规则划分为若干个网格，用卫星对网格点的覆盖能力来表示卫星对整个网格区域的覆盖能力.文献[5]在分析星座对地覆盖性能时，也是将目标区域进行均匀网格划分，计算星座在任意时刻对目标区域的覆盖情况.文献[6]通过将空间球面圆进行网格划分，分析了低轨凝视传感器的空间球面覆盖性能，文献[7]分别在近空间平台对地区域覆盖优化设计和考虑多项卫星配置参数下通信卫星星座设计中，将地球表面进行网格划分，分析航天器覆盖性能.文献[8]采用网格点抽样的方法进行覆盖率计算，并从理论上推导了网格点法划分大小和计算误差之间的关系.但是，通过将目标区域划分为足够多的点进行星座对地覆盖计算与覆盖性能分析，主要缺点是计算效率低，且不能保证计算结果的可靠性.

文献[9]通过对星座空间几何构形进行划分，研究了全球连续性完全覆盖下的Rosette星座优化设计方法.文献[10]采用改进网格点法进行星座对地面目标覆盖性能计算，研究了通用航空飞行器组网星座的优化设计.文献[11]提出一种半长轴分析方法进行全球连续性覆盖周期性轨道设计，该算法通过绘制卫星每次经过不同纬线圈时的升交点和降交点，然后进行统计排序，从而判断是否存在缝隙.这几种方法在计算效率上都有所提高，但只局限于对全球目标的完全覆盖分析.文献[12]通过定义覆盖状态函数和覆盖区域函数，分析了星座对区域的完全覆盖问题和连续性覆盖问题，提出了卫星在固定时刻对纬线的覆盖范围计算方法，消除了抽样对覆盖结果的可信度影响，但该方法同样要基于经度条带的划分精度.

根据卫星对地覆盖几何特性，基于球面Delaunay三角网络和Voronoi图的性质以及快速构建方法，本文定义了卫星的所属覆盖区域，并对该区域做了定性与定量描述，从而将多星覆盖问题转换为单星覆盖问题，旨在提出满足覆盖约束下，星座对任意类型地面目标区域的确定性计算方法，从而精确、快速地计算星座对目标区域的覆盖能力.

1 星座对目标区域覆盖能力分析

1.1 卫星对地覆盖几何特性

已知地球半径为Re，卫星轨道高度为h，受地面目标最小观测仰角γmin的约束，卫星对地的最大覆盖半幅宽为[13]

同时，对于一个由轨道高度相同的圆轨道卫星构成的星座，已知A,B,C是该星座中任意3颗星的星下点，要保证球面ΔABC内任意一点都能被覆盖，只需要确保ΔABC内距离3个顶点最远的点P被覆盖.可见，将1个由N颗星构成的星座，在任意状态下的空间几何构形剖分成若干个相邻且互不相交的球面三角形.计算三角网中所有球面三角形内最远点P到顶点的半径ri，如果满足[10]:

则卫星星座在当前状态下可以实现一重全球完全覆盖，若在整个星座重构周期内的任意时刻该条件都成立，则表明该卫星星座能够实现一重连续性全球完全覆盖.

1.2 星座对地覆盖分析

定义1 卫星的邻接结点.球面三角网剖分结果中，将包含卫星星下点的球面三角形的外接圆圆心记作该卫星的邻接结点.

定义2 卫星的邻接多边形.以卫星星下点为中心，依次顺序连接该点的所有邻接结点围成的多边形区域.

定义3 卫星所属覆盖区域.卫星的邻接多边形区域和目标覆盖区域的相交区域.

已知卫星星座C={S1,S2,…,SN}，通过对某一时刻星座在地球表面上投影点的空间几何构形进行球面Delaunay三角网剖分和Voronoi图划分，其性质决定了当前时刻星座中如果存在1颗卫星Si不能实现对其所属覆盖区域的完全覆盖，则对于任意的Sj都不能对Si的所属覆盖区域的漏缝进行覆盖，即该状态下，星座必然不能实现对目标区域的完全覆盖.因此，要实现对目标区域的完全覆盖，则整个星座中所有卫星均能实现对其所属覆盖区域的完全覆盖.

2 卫星对其所属覆盖区域的覆盖计算

由上述分析可见，星座对地覆盖问题可以转换为多个独立的单星对其所属覆盖区域的覆盖问题，并且生成的卫星所属覆盖区域是球面大圆弧段和小圆弧段围城的多边形区域，因此，卫星对其所属覆盖区域的覆盖计算，就可以归结为卫星实际覆盖圆与其所属覆盖区域的相交性的分析和两球面圆的相交区域计算.

2.1 卫星对其所属覆盖区域的覆盖能力计算

如果d<α且α-d

如果d≥α，即O2点在圆Θ1外部.此时，当两圆相交且满足cosd=cosα·cosr，则可得到两个球面直角三角形O1MO2和O1NO2，其中，M,N为两圆的交点.

当d-α

当acos(cosd/cosα)≤r<α+d，如图3所示，区域R如图3中阴影部分所示，则

图1 球面扇形面积计算

图2 球面两弧度围成面积计算

图3 球面两弧度围成面积计算

2.2 星座对地覆盖率计算

2.2.1 覆盖目标是全球区域

将全球区域按球面Voronoi图剖分方法进行划分，得到每颗卫星所属覆盖区域，如图4所示，此时，卫星Si的所属覆盖区域Ωi就是其所属多边形区域Ai.图中，由顶点p1,p2,p3,p4,p5围成的球面多边形区域就是卫星Si的所属覆盖区域Ωi，阴影部分是卫星对其所属覆盖区域的实际覆盖区域.

图4 全球区域覆盖计算

当卫星覆盖圆半径r=r1时，卫星实际覆盖区域的边界顶点集合Ti=Φ，覆盖区域为以O为圆心，r1为半径的球面圆.卫星Oi对其所属覆盖区域Ωi的面积为

2.2.2 覆盖目标是球面凸多边形区域

将球面凸多边形区域按球面Voronoi图剖分方法进行划分，得到每颗卫星所属覆盖区域，如图5所示.当卫星星下点在覆盖目标区域Ω内部时，卫星对其所属覆盖区域的覆盖计算和星座对全球覆盖分析过程一样.当卫星星下点在覆盖区域外部时，卫星Si的所属覆盖区域Ωi是由顶点p1,p2,p3,p4,p5围成的球面多边形区域.

图5 球面多边形区域覆盖计算

当卫星覆盖圆半径r=r1时，卫星实际覆盖区域的边界顶点集合Ti=Φ，卫星不能实现对其所属覆盖区域的覆盖，因此卫星Oi对其所属覆盖区域Ωi的面积为

2.2.3 覆盖目标是球面圆区域

将以O为圆心，α为半径的球面圆形区域按球面Voronoi图剖分方法进行划分，得到每颗卫星所属覆盖区域，如图6所示.当卫星星下点在覆盖目标区域Ω内部，即d≤α时，卫星对其所属覆盖区域的覆盖计算和星座对全球覆盖分析过程一样，与其不同的是卫星的实际覆盖区域边界中存在小圆弧段.当卫星星下点在覆盖区域外部，即d>α.

图6 球面圆区域覆盖计算

当卫星覆盖圆半径r≤d-α时，卫星实际覆盖区域的边界顶点集合Ti=Φ，卫星Oi对其所属覆盖区域Ωi的面积为

2.2.4 覆盖目标是纬度带区间

纬度带区间为[l,u].当下纬度线l≥0时，星座对纬度带区域的覆盖计算可以分解为对以(0,π)为圆心，π/2-l为半径的圆的覆盖面积Sl，和以(0,π)为圆心，π/2-u为半径的圆的覆盖面积Su，星座对纬度带区域的覆盖面积S=Sl-Su.

当下纬度线u≤0时，星座对纬度带区域的覆盖计算可以分解为对以(0,π)为圆心，π/2-l为半径的圆的覆盖面积Sl，和以(0,π)为圆心，π/2-u为半径的圆的覆盖面积Su，星座对纬度带区域的覆盖面积S=Su-Sl.

当l<0，且u>0时，星座对纬度带区域的覆盖计算可以分解为对全球的覆盖面积SA和以(0,π)为圆心，π/2-l为半径的圆的覆盖面积Sl和以(0,π)为圆心，π/2-u为半径的圆的覆盖面积Su，星座对纬度带区域的覆盖面积S=SA-Sl-Su.

3 实例分析

以参数为(T/P/F)=(36/4/1)的Walker星座为例，卫星轨道高度为H=1 300 km.定义覆盖目标1为球面多边形区域，球面多边形的顶点集合为(63°,54°),(73°,3°), (126°,3°), (156°,44°), (120°,68°).覆盖目标2为球面圆区域，球面圆中心为(θ,φ)=(63°,54°)，半径为α=20°.卫星对目标最小观测仰角约束为γmin=10°，分别计算在一个卫星轨道周期内，该星座对该球面多边形区域和球面圆形区域的覆盖能力.

3.1 星座对地覆盖率计算

首先，对星座的空间几何构形进行球面Delaunay三角网剖分，并在此基础上，计算卫星的所属覆盖区域.以覆盖目标是球面圆区域为例.

图7所示是网格划分精度分别为5.0°,1.0°,0.1°时，整个仿真周期内星座对目标区域的覆盖率计算结果与本文所提方法计算结果的对比.图中，黑色实线为采用本文提出的球面几何剖分法计算得到的星座对目标区域的覆盖率，红色实线为采用网格点划分法计算得到的星座对目标区域的覆盖率.

比较图7中两种方法的计算结果，并计算对应的覆盖率偏差，结果如图8所示.

分别对比图8的计算结果，可见网格划分精度越小，网格点法计算结果越精确，且和本文所提出方法的计算结果偏差也越小，反映了本文所提出方法的正确性.

图7 不同网格划分精度下覆盖率计算结果对比

图8 不同网格划分精度下覆盖率计算结果差别

当网格点划分精度为1.0°时，每个网格将代表地球表面上面积为(π·Re/180)2km2的区域.将目标区域按照网格点划分进行星座对地覆盖率计算，分析进行一次星地覆盖计算的计算时长随目标区域面积和划分精度的关系.当网格划分精度为0.1°，覆盖目标是全球区域时，进行一次覆盖计算的时间约为88.3 s；而采用本文提出的星座空间几何构形剖分的方法进行一次星地覆盖计算只需要约21 ms，并且目标区域的形状和大小对计算时长均没有影响，反映了本文所提出方法在时间复杂度上的高效性.

3.2 结果分析

通过对上述实验的对比分析，计算结果验证了本文基于星座空间构形几何剖分进行星座对地覆盖计算方法的正确性和高效性.尤其当目标区域较大且对星座覆盖精度要求较高时，传统的基于网格点划分的方法往往是不可行的.并且在仿真周期内，当星座中存在两颗星在某一时刻相撞时，通过本文的星座空间构形几何剖分法可以检测出来，且不影响剖分方法和计算结果.

4 结 论

1)采用经典的网格点或者条带划分法进行覆盖率计算，得到的星座对地覆盖率精度是基于网格点或者条带的划分精度的；而网格点或条带划分的越细，对应的计算时间复杂度会很高.对比之下，本文提出的基于星座空间构形几何剖分进行星座对地覆盖计算是一种确定性方法，计算结果是精确值，更加可靠，且计算效率有显著提高.

2)基于Delaunay三角网和Voronoi图特性，对任意时刻卫星星座在地球表面投影点的集合进行划分，且通过定义每颗卫星所属覆盖区域，将星座覆盖问题简化为单星覆盖问题，可以精确快速计算星座中卫星对目标区域的覆盖能力的时变特性，以及在实现对目标区域的覆盖过程中每颗星的重要性，并且该方法适用于对任意类型的地面目标区域覆盖计算.

[1] SANDAU R, BRIEß K, D’ERRICO M. Small satellites for global coverage: potential and limits[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2010, 65(6): 492-504. DOI：10.1016/j.isprsjprs.2010.09.003.

[2] ARGOUN M B. Recent design and utilization trends of small satellites in developing countries[J]. Acta Astronautica, 2012, 71: 119-128. DOI: 10.1016/j.actaastro.2011.07.024.

[3] ULYBYSHEV Y. Satellite constellation design for complex coverage[J]. Journal of Spacecraft and Rockets, 2008, 45(4): 843-849.DOI:10.2514/1.35369.

[4] MORRISON J J. A system of sixteen synchronous satellites for worldwise navigation and surveillance[M]. Fort Belvoir, VA：Defense Technical Information Center, 1973.

[5] JIANG Yong, YANG Sen, ZHANG Gengxin, et al. Coverage performances analysis on Combined-GEO-IGSO satellite constellation[J]. Journal of Electronics (China), 2011, 28(2): 228-234. DOI: 10.1007/s11767-011-0581-1.

[6] 邓勇，王春明，张中兆. 红外低轨星座凝视传感器的空间覆盖性能分析[J]. 宇航学报, 2011, 32(1): 123-128.DOI: 10.3873/j.issn.1000-1328.2011.01.019. DENG Yong, WANG Chunming, ZHANG Zhongzhao. Analysis on coverage performance of staring sensors infrared LEO constellation[J]. Journal of Astronautics, 2011, 32(1): 123-128.DOI: 10.3873/j.issn.1000-1328.2011.01.019.

[7] MORTARI D, DE SANCTIS M, LUCENTE M. Design of flower constellations for telecommunication services[J]. Proceedings of the IEEE, 2011, 99(11): 2008-2019. DOI: 10.1109/JPROC.2011.2158766.

[8] 秦睿杰，戴光明，王茂才，等. 一种计算星座区域覆盖率的高效抽样网格点法[J]. 计算机应用研究, 2015, 32(4): 1065-1068.DOI: 10.3969/j.issn.1001-3695.2015.04.025. QIN Ruijie, DAI Guangming, WANG Maocai, et al. Efficient sampling grid-point approach for calculating regional coverage of satellite constellation[J]. Application Research of Computers, 2015, 32(4): 1065-1068.DOI: 10.3969/j.issn.1001-3695.2015.04.025.

[9] BALLARD A H. Rosette constellations of earth satellites[J]. IEEE Transactions on Aerospace and Electronic Systems, 1980, 16 (5): 656-673. DOI: 10.1109/TAES.1980.308932.

[10]王绍凯，崔红正，韩潮. 超地平覆盖飞行器组网星座优化设计[J]. 哈尔滨工业大学学报, 2013, 45(7): 109-114. WANG Shaokai, CUI Hongzheng, HAN Chao. Constellation optimization design for hyper-horizon coverage vehicles[J]. Journal of Harbin Institute of Technology,2013, 45(7): 109-114.

[11]XU Ming, HUANG Li. An analytic algorithm for global coverage of the revisiting orbit and its application to the CFOSAT satellite[J]. Astrophysics and Space Science, 2014, 352(2): 497-502. DOI: 10.1007/s10509-014-1939-2.

[12]宋志明，戴光明，王茂才，等. 卫星星座对地面目标的连续性覆盖分析[J]. 华中科技大学学报(自然科学版), 2014, 42(8): 33-37.DOI: 10.13245/j.hust.140807. SONG Zhiming, DAI Guangming ,WANG Maocai, et al. Continuous coverage analysis of satellite constellation to ground target[J]. Journal of Huazhong University of Science and Technology (Nature Science Edition) , 2014, 42(8): 33-37.DOI: 10.13245/j.hust.140807.

[13]ULYBYSHEV Y. Geometric analysis of low-earth-orbit satellite communication systems: covering functions[J]. Journal of Spacecraft and Rockets, 2000, 37(3): 385-391. DOI: 10.2514/2.3572.

(编辑 张 红)

Deterministic method for coverage of constellation to ground region

CHEN Xiaoyu1, DAI Guangming1,2, WANG Maocai1,2, SONG Zhiming1,2

(1.School of Computer, China University of Geosciences,Wuhan 430074, China; 2.Hubei Key Laboratory of Intelligent Geo-Information Processing (China University of Geosciences), Wuhan 430074, China)

To improve the coverage accuracy and efficiency of the satellite, the geometric characteristic related to satellite coverage of a ground region is analyzed, and the division method for space geometric configuration of constellation is mainly studied based on the properties of the Delaunay triangulation and the Voronoi Diagram. Subsequently, the constellation coverage problem is simplified into multiple one-satellite coverage problems by defining the coverage area of each satellite. Furthermore, an accurate and effective approach to solve the satellite constellation coverage problem for any type of coverage area is put forward. Finally, the comparison of accuracy and efficiency between the proposed approach and the classic Grid-point method is given for different target areas. The test result shows that the proposed method is valid and more efficient. The proposed method can also be used in satellite constellation optimization design and performance evaluation.

satellite constellation; coverage calculation; performance analysis; delaunay triangulation; Voronoi diagram

10.11918/j.issn.0367-6234.201509008

2015-09-02

国家自然科学基金(41571403,61472375,61103144)

陈晓宇(1990—),男，博士研究生; 戴光明(1964—)，男，教授，博士生导师

戴光明，gmdai@cug.edu.cn

V474

A

0367-6234(2017)04-0055-06