一种全天球星图识别方法技术

一种全天球星图识别方法。首先提取拍摄星图图像中距离最远的两个星点A、B；在星图图像中除A、B两点外的其它星点中，提取与A、B距离之和最远的星点C；提取与A、B、C距离之和最远的星点D。在导航星表1中进行查找，满足与A、B、C、D三个星点距离及能量匹配的星图区域，作为初选区域。在拍摄的星图内提取出的星点，如果待识别星与图边界满足一定范围，则在该星周围若干范围内，分别取三颗星形成距离角度组合，用于识别该颗星。依次识别完所有拍摄的星图内的星点后，在识别出的星点中，选取距离其他星夹角最近的一颗星，以该星为中心，角距T度范围外的星剔除；保存剩余的识别星集合，识别过程结束。

A recognition method for all day star chart

An all day star chart recognition method. Firstly, the two farthest stars A and B are extracted from the photographed star image; in the other stars except A and B, the farthest stars C is extracted from the sum of A and B; and the farthest stars D is extracted from the sum of A, B and C. In the navigation star table 1, the star map regions matching the distance and energy of A, B, C and D are selected as the primary regions. If the boundaries of the star and the map satisfy a certain range, three stars are selected to form distance and angle combinations around the star to identify the star. After identifying all the stars in the star map in turn, a star nearest to the angle of other stars is selected to be the center of the star, and the stars outside the T degree range are removed. The remaining recognition star sets are saved and the recognition process is completed.

【技术实现步骤摘要】
一种全天球星图识别方法
本专利技术涉及一种星图识别方法，属于星敏感器设计

技术介绍
在研制恒星陀螺敏感器过程中，精度要求高，星表内恒星数量多。利用原有的全天球导航星识别方法，识别速度较慢，难以快速识别出恒星给出姿态，因此针对大导航星表，需要研制快速全天球识别算法。本研究基于多边形识别算法和网格识别算法的优点，面向导航星较多的实际情况，提出了四边形网格法、进行全天球恒星的识别。试验结果表明，该算法在DSP平台2秒内能够完成全天球恒星的识别，识别概率达到99.98％以上；移植到FPGA平台后，将在1秒内完成全天球识别。该算法具有速度快、识别概率高的优点，能够快速给空间飞行器提供姿态数据。郑胜等2004年在光学技术的《一种基于三角形几何结构的星图识别算法》中提出了三角形匹配法。该方法的基本形式是：首先在观测星和导航星中根据角距和星等等特征匹配出3颗星组成主三角形，再与这3颗星之间的角距关系来匹配星图中其它恒星，取最大匹配星组作为识别结果。该算法的优点是实现简单，占存储容量小，利用了整幅星图特征，不易出现误匹配，可以用于全天球星图识别。该方法有多种变形，如基于主边的匹配法、基于主星的匹配法和双三角形匹配法(准四边形法)。由于匹配其它恒星时，要计算这些恒星的星对角距(恒星与恒星之间与地球的球心角)，所以该方法识别时间长。1997年Mortari在JournaloftheAstronauticalSciences的《Search-LessAlgorithmforStarPatternRecognition》中提出采用k-vector方法来进行星图识别。为了降低由于星等不确定性带来的误识别，该算法只采用星对角距来识别。在研究过程中，发现星对信息表中的星对角距递增的幅值是比较平缓的，于是考虑将星对信息表的星对序号和对应的星对角距两个信息分别作为纵坐标和横坐标，拟合成某种曲线函数y＝f(x)的形式，这样，当从星敏感器拍摄到的星空图像算出一对观测星对的角距时，令其为x代入式中，计算得到星对序号y值，再根据y值来确定识别结果。该算法起到了快速初定位的作用，减少了查找次数。2001年Mortari在AASGuidance&ControlConference的《Lost-In-SpacePyramidAlgorithmforRobustStarPatternRecognition》中提出了一种金字塔算法(PyramidStarIdentification)。该算法是适当选取四颗测量星，以一颗观测星作为顶点以其余三颗观测星构成四面体。该算法以四面体为基本单位，以k-vector方法为基本算法。实验发现：随着星对信息表增加，拟合曲线y＝f(x)精度降低。而由于拟合精度降低，不能确保最佳匹配星对落在(y–δ)和(y+δ)之间的子集(δ为角距误差)。因此该算法只在局部区域内起到了快速初定位的作用。李春艳等2003年在科学通报E辑的《基于神经网络技术的星图识别方法》中提出了神经网络法。该方法利用神经网络和专家系统的优点，可以提高识别精度，且能降低由于星等测量精度引起的误差。神经网络法跟传统的算法相比，具有数据库容量低，实时性和鲁棒性好等优点，在大视场多星定位系统中表现地尤其明显。但这些算法训练需要较大计算量，要求很大的训练集合，以完成多种模式识别，因此具有识别时间长的缺点；精度易受到训练集合大小和训练时间长短的影响；潜在地需要较大存储量来存储权值，所以对硬件要求较高。
技术实现思路
本专利技术的技术解决问题是：克服现有技术的不足，提供一种全天球星图识别方法，通过本专利技术的方法进行星点识别，用于计算星敏感器的姿态，提高了星敏感器图像识别精度及姿态计算的速度。本专利技术的技术解决方案是：一种全天球星图识别方法，包括步骤如下：(14)提取拍摄的观测星的星图图像中距离最远的两个星点：星点A、星点B；(15)在星图图像中除星点A、星点B外的其它星点中，提取与星点A、星点B距离之和最大的星点C；(16)在星图图像中除星点A、星点B、星点C外的其它星点中，提取与星点A、星点B、星点C距离之和最大的星点D；(17)在导航星表1中，查找星点两端距离与星点A到星点B距离匹配且星点两端能量分别和星点A、星点B匹配的若干个星点Aˊ和星点Bˊ；(18)对于每个星点Aˊ和星点Bˊ，在导航星表1中寻找若干个星点Cˊ，所述星点Cˊ、星点Aˊ之间的距离与星点C、星点A之间的距离相匹配，星点Cˊ、星点Bˊ之间的距离与星点C、星点B之间的距离相匹配，且星点Cˊ的能量和星点C的能量相匹配；(19)对于每组星点Aˊ、星点Bˊ和星点Cˊ，在导航星表1中寻找若干个星点Dˊ，所述星点Dˊ、星点Aˊ之间的距离与星点D、星点A之间的距离相匹配，星点Dˊ、星点Bˊ之间的距离与星点D、星点B之间的距离相匹配，星点Dˊ、星点Cˊ之间的距离与星点D、星点C之间的距离相匹配，，且星点Dˊ的能量和星点D的能量匹配；上述寻找到的星点Aˊ、星点Bˊ、星点Cˊ、星点Dˊ组合构成若干个初选星空区域；(20)将从拍摄的星图内提取出的星点按照能量从大到小排序，保存到待识别星表中；(21)在待识别星表中，任意取出一颗待识别星点E，判断待识别星点E到星图图像上下左右边界的最短距离是否大于等于设定的dist像素；(22)如果待识别星点E到星图图像上下左右边界的最短距离小于设定的dist像素，则在待识别星表中取下一颗待识别星点，返回步骤(8)；如果待识别星点E到星图图像上下左右边界的最短距离都大于等于dist像素，则待识别星点E周围N1到N2个像素内，取亮度前三的三个星点F1、F2、F3；星点E周围N3到N4个像素内，取亮度前三的三个星点G1、G2、G3；星点E周围N5到N6个像素内，取亮度前四的四个星点H1、H2、H3、H4；N1、N2、N3、N4、N5、N6为正整数；N6>N5；N4>N3；N2>N1；(23)遍历星点E和星点Fi、Gj、Hk，获得星点E的标识，所述星点E的标识由5个值构成，分别为：星点Fi与星点E的距离、星点Gj与星点E的距离、星点Hk与星点E的距离、星点FiEGj的夹角、星点GiEHk的夹角；i＝1,2,3；j＝1,2,3；k＝1,2,3,4；(24)将获得的多组星点E的标识与导航星表2对应的初选星空区域中的导航星信息进行匹配，获得星点E的唯一匹配标识；(25)取下一个待识别星点，返回步骤(8)，直至识别完M颗观测星，进入步骤(13)；M为正整数；(26)在识别出的星点中，选取距离其他星最近的一颗星，以该星为中心，将角距T度范围外的识别星剔除，保存剩余的识别星集合作为最终的识别星，根据最终的识别星计算星敏感器的姿态；T为正整数。两个星点间的距离为：两个星点能量相乘后，再乘以两个星点间像素距离值；两个星点的夹角为两个星点能量相乘后再乘以两个星点物理夹角值。T在7～10范围内选取。设定的dist像素>N6。所述导航星表1中存储的星点信息为星点之间的距离和各星点的能量。所述导航星表2中存储的星点信息包括标识星与该标识星周围三颗星之间的距离、该标识星周围三颗星与标识星构成的两个夹角。本专利技术与现有技术相比的有益效果是：(1)本专利技术的整体本文档来自技高网...

【技术保护点】
1.一种全天球星图识别方法，其特征在于，包括步骤如下：(1)提取拍摄的观测星的星图图像中距离最远的两个星点：星点A、星点B；(2)在星图图像中除星点A、星点B外的其它星点中，提取与星点A、星点B距离之和最大的星点C；(3)在星图图像中除星点A、星点B、星点C外的其它星点中，提取与星点A、星点B、星点C距离之和最大的星点D；(4)在导航星表1中，查找星点两端距离与星点A到星点B距离匹配且星点两端能量分别和星点A、星点B匹配的若干个星点Aˊ和星点Bˊ；(5)对于每个星点Aˊ和星点Bˊ，在导航星表1中寻找若干个星点Cˊ，所述星点Cˊ、星点Aˊ之间的距离与星点C、星点A之间的距离相匹配，星点Cˊ、星点Bˊ之间的距离与星点C、星点B之间的距离相匹配，且星点Cˊ的能量和星点C的能量相匹配；(6)对于每组星点Aˊ、星点Bˊ和星点Cˊ，在导航星表1中寻找若干个星点Dˊ，所述星点Dˊ、星点Aˊ之间的距离与星点D、星点A之间的距离相匹配，星点Dˊ、星点Bˊ之间的距离与星点D、星点B之间的距离相匹配，星点Dˊ、星点Cˊ之间的距离与星点D、星点C之间的距离相匹配，，且星点Dˊ的能量和星点D的能量匹配；上述寻找到的星点Aˊ、星点Bˊ、星点Cˊ、星点Dˊ组合构成若干个初选星空区域；(7)将从拍摄的星图内提取出的星点按照能量从大到小排序，保存到待识别星表中；(8)在待识别星表中，任意取出一颗待识别星点E，判断待识别星点E到星图图像上下左右边界的最短距离是否大于等于设定的dist像素；(9)如果待识别星点E到星图图像上下左右边界的最短距离小于设定的dist像素，则在待识别星表中取下一颗待识别星点，返回步骤(8)；如果待识别星点E到星图图像上下左右边界的最短距离都大于等于dist像素，则待识别星点E周围N1到N2个像素内，取亮度前三的三个星点F1、F2、F3；星点E周围N3到N4个像素内，取亮度前三的三个星点G1、G2、G3；星点E周围N5到N6个像素内，取亮度前四的四个星点H1、H2、H3、H4；N1、N2、N3、N4、N5、N6为正整数；N6>N5；N4>N3；N2>N1；(10)遍历星点E和星点Fi、Gj、Hk，获得星点E的标识，所述星点E的标识由5个值构成，分别为：星点Fi与星点E的距离、星点Gj与星点E的距离、星点Hk与星点E的距离、星点Fi E Gj的夹角、星点Gi E Hk的夹角；i＝1,2,3；j＝1,2,3；k＝1,2,3,4；(11)将获得的多组星点E的标识与导航星表2对应的初选星空区域中的导航星信息进行匹配，获得星点E的唯一匹配标识；(12)取下一个待识别星点，返回步骤(8)，直至识别完M颗观测星，进入步骤(13)；M为正整数；(13)在识别出的星点中，选取距离其他星最近的一颗星，以该星为中心，将角距T度范围外的识别星剔除，保存剩余的识别星集合作为最终的识别星，根据最终的识别星计算星敏感器的姿态；T为正整数。...

【技术特征摘要】
1.一种全天球星图识别方法，其特征在于，包括步骤如下：(1)提取拍摄的观测星的星图图像中距离最远的两个星点：星点A、星点B；(2)在星图图像中除星点A、星点B外的其它星点中，提取与星点A、星点B距离之和最大的星点C；(3)在星图图像中除星点A、星点B、星点C外的其它星点中，提取与星点A、星点B、星点C距离之和最大的星点D；(4)在导航星表1中，查找星点两端距离与星点A到星点B距离匹配且星点两端能量分别和星点A、星点B匹配的若干个星点Aˊ和星点Bˊ；(5)对于每个星点Aˊ和星点Bˊ，在导航星表1中寻找若干个星点Cˊ，所述星点Cˊ、星点Aˊ之间的距离与星点C、星点A之间的距离相匹配，星点Cˊ、星点Bˊ之间的距离与星点C、星点B之间的距离相匹配，且星点Cˊ的能量和星点C的能量相匹配；(6)对于每组星点Aˊ、星点Bˊ和星点Cˊ，在导航星表1中寻找若干个星点Dˊ，所述星点Dˊ、星点Aˊ之间的距离与星点D、星点A之间的距离相匹配，星点Dˊ、星点Bˊ之间的距离与星点D、星点B之间的距离相匹配，星点Dˊ、星点Cˊ之间的距离与星点D、星点C之间的距离相匹配，，且星点Dˊ的能量和星点D的能量匹配；上述寻找到的星点Aˊ、星点Bˊ、星点Cˊ、星点Dˊ组合构成若干个初选星空区域；(7)将从拍摄的星图内提取出的星点按照能量从大到小排序，保存到待识别星表中；(8)在待识别星表中，任意取出一颗待识别星点E，判断待识别星点E到星图图像上下左右边界的最短距离是否大于等于设定的dist像素；(9)如果待识别星点E到星图图像上下左右边界的最短距离小于设定的dist像素，则在待识别星表中取下一颗待识别星点，返回步骤(8)；如果待识别星点E到星图图像上下左右边界的最短距离都大于等于dist像素，则待识别星点E周围N1到N2个像素内，取亮度前三的三个星点F1、F2、F3；星点E周围N3到N4个像素...

【专利技术属性】
技术研发人员：刘达，郝云彩，张春明，张俊，米强，程会艳，张洪健，李春艳，
申请(专利权)人：北京控制工程研究所，
类型：发明
国别省市：北京,11