发明名称 基于周长特征的暗星自主识别方法
摘要 本发明涉及一种基于周长特征的暗星自主识别方法,包括以下步骤:一、根据基本星表,采用星等阈值法构建导航星表。采用重叠球矩法将全天球划分成若干局部重叠天区,找出每个重叠天区中最亮的两颗星作为主星对,其余星作为辅助星构建三角形周长特征库,对生成的特征库按周长大小构造散列函数,进行分段存储。二、读入导航星表和导航特征库,星图提取结果按灰度值大小降序排列,以最亮的两颗星为主星对,其余星为辅助星,构造以主星对为公共边的观测三角形,计算三角形周长和以主星为公共点的三角形两边长。本发明大大减少了特征三角形的数量,相对传统基于星角距的三角形识别算法,减少了星角距比较次数,提高了识别速度。
申请公布号 CN106123891A 申请公布日期 2016.11.16
申请号 CN201610457479.0 申请日期 2016.06.22
申请人 中国人民解放军63680部队 发明人 张同双;郭敬明;周海渊;钟德安;潘良;焦宏伟;王二建
分类号 G01C21/02(2006.01)I 主分类号 G01C21/02(2006.01)I
代理机构 江阴市同盛专利事务所(普通合伙) 32210 代理人 唐纫兰;沈国安
主权项 一种基于周长特征的暗星自主识别方法,其特征在于:所述方法包括以下步骤:步骤一、导航星特征库的构建A、导航星的筛选⑴采用星等阈值法值剔除基本星表中大于星敏感器极限探测星的恒星;⑵双星处理当两颗恒星间的星角距小于弥散斑时,将双星当作一颗星来处理,并对恒星赤经 、赤纬及星等进行合成,则合成后的恒星赤经<img file="dest_path_image001.GIF" wi="17" he="16" />、赤纬<img file="232357dest_path_image002.GIF" wi="16" he="18" />及星等<img file="dest_path_image003.GIF" wi="18" he="16" />分别为:<img file="523399dest_path_image004.GIF" wi="197" he="75" />其中,弥散斑大小在3×3~5×5像元之间,(<img file="dest_path_image005.GIF" wi="19" he="20" />,<img file="240819dest_path_image006.GIF" wi="18" he="24" />)和(<img file="dest_path_image007.GIF" wi="21" he="20" />,<img file="573712dest_path_image008.GIF" wi="18" he="24" />)分别为两颗恒星的赤经和赤纬,<img file="dest_path_image009.GIF" wi="21" he="20" />和<img file="827232dest_path_image010.GIF" wi="22" he="20" />分别为两颗恒星的星等;⑶ 导航星表的生成对选定的恒星重新编号,并将每个恒星的星表号、赤经、赤纬和星等进行存储;B、导航特征库的构建结合主星算法构建特征三角形:以视场中的最亮星与次亮星作为主星对,其余星作为辅助星,分别与主星对构成特征三角形,每个特征三角形包括:主星星号、次星星号、辅助星星号、主星与次星间星角距、主星与辅助星间星角距和三角形周长这些数据,⑴ 分区星表的构建采用重叠球矩法将整个天球分割成若干个局部重叠天区,球矩大小应满足如下条件:<img file="dest_path_image011.GIF" wi="240" he="51" />式中,<img file="587377dest_path_image012.GIF" wi="24" he="25" />、<img file="dest_path_image013.GIF" wi="22" he="25" />分别为赤经范围与赤纬范围,<img file="792094dest_path_image014.GIF" wi="21" he="21" />、<img file="dest_path_image015.GIF" wi="19" he="24" />为星敏感器视轴指向,<img file="364895dest_path_image016.GIF" wi="17" he="18" />为视场半径,计算球矩大小时需对<img file="dest_path_image017.GIF" wi="24" he="25" />过零处理,赤纬最大值、最小值分别设为90°和‑90°,建立分区星表时,以星敏感器视轴为中心,计算球矩大小,将球矩内的所有恒星按星等升序排列,并存储到以分区号为文件名的文本文件中,分区索引号的计算公式如下所示,<img file="846823dest_path_image018.GIF" wi="237" he="22" />⑵ 特征库的建立读取分区星表,以选定的主星对分别与其它恒星构建特征三角形,读取相关星表号、赤经、赤纬和星等信息,则根据下式计算相关星角距<img file="dest_path_image019.GIF" wi="21" he="25" />,<img file="279335dest_path_image020.GIF" wi="208" he="30" />式中,“<img file="dest_path_image021.GIF" wi="9" he="9" />”为矢量点积运算符;“<img file="909031dest_path_image022.GIF" wi="22" he="27" />”为<img file="dest_path_image023.GIF" wi="14" he="14" />的范数,<img file="583726dest_path_image024.GIF" wi="19" he="23" />和<img file="dest_path_image025.GIF" wi="22" he="25" />分别为恒星<img file="dest_path_image027.GIF" wi="9" he="18" />、<img file="608576dest_path_image028.GIF" wi="14" he="19" />的参考矢量,构建特征三角形时,星角距<img file="38420dest_path_image019.GIF" wi="21" he="25" />应满足如下条件:<img file="dest_path_image029.GIF" wi="115" he="29" />式中,<img file="529313dest_path_image030.GIF" wi="14" he="21" />为弥散斑所对应的角度,若三角形的三条边均满足如上星角距限制条件,则按照导航特征数据结构存储该三角形;⑶ 周长散列特征库的构建判断所有主星是否计算完毕,否则重复⑵,是则对所有特征三角形按照周长大小分块存储,子块编号按下式构造散列函数<img file="dest_path_image031.GIF" wi="40" he="24" />:<img file="169549dest_path_image032.GIF" wi="203" he="24" />(1)式中,<img file="dest_path_image033.GIF" wi="17" he="18" />为三角形周长,<img file="688386dest_path_image034.GIF" wi="19" he="19" />为分块总数;步骤二、星图识别⑴ 初始化:读入导航星表、导航三角形特征库,设置星敏感器参数和星角距匹配门限<img file="dest_path_image035.GIF" wi="14" he="16" />参数;⑵ 对星图提取数据进行编号,并按灰度值降序排列;⑶ 判断是否存在3颗以上观测星,小于3颗则不能识别;⑷ 若观测星数为3,则按灰度值确定主星对与辅助星,构建观测三角形,计算观测三角形的对应星角距和周长,根据式(1)定位子块编号,在该子块中搜索以主星为顶点的两个星角距是否存在唯一满足下式(2)的三角形,是则将匹配结果存入识别矩阵中,其中第<img file="459771dest_path_image027.GIF" wi="9" he="18" />个元素表示星点<img file="188692dest_path_image027.GIF" wi="9" he="18" />的星表号,识别失败时则该元素置0;<img file="205190dest_path_image036.GIF" wi="106" he="51" />(2)式中,A、B、C分别为主星、次星及辅星,<img file="dest_path_image037.GIF" wi="25" he="22" />表示<img file="640850dest_path_image023.GIF" wi="14" he="14" />的绝对值;⑸ 若观测星大于3颗,则分别以星点1、2、3和星点1、2、4构建两个观测三角形,并按步骤⑷分别进行识别,若两个观测三角形存在唯一相同的主星对匹配结果,则匹配成功,将前4颗星的匹配结果存入识别矩阵中,然后根据已匹配观测星对星图中其它恒星进行投影识别,并利用观测星角距对识别结果进行验证;若两个观测三角形没有唯一匹配结果,则剔除星点3或星点4,读入下一颗星继续识别;⑹ 判断该帧星图是否识别完毕,否则转⑸,是则根据识别矩阵和导航星表确定观测星的各个参数,转入下一帧星图识别,直至整个识别任务完毕。
地址 214400 江苏省无锡市江阴市103信箱850号