CN102840860A - 一种基于混合粒子群算法的星图识别方法 - Google Patents
一种基于混合粒子群算法的星图识别方法 Download PDFInfo
- Publication number
- CN102840860A CN102840860A CN2012103151841A CN201210315184A CN102840860A CN 102840860 A CN102840860 A CN 102840860A CN 2012103151841 A CN2012103151841 A CN 2012103151841A CN 201210315184 A CN201210315184 A CN 201210315184A CN 102840860 A CN102840860 A CN 102840860A
- Authority
- CN
- China
- Prior art keywords
- star
- nautical
- nautical star
- asterism
- circle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Navigation (AREA)
Abstract
本发明公开了一种基于混合粒子群算法的星图识别方法,包括建立导航星库和基于混合粒子群算法进行星图识别两个步骤,圆半径r通过自适应调整方法来确定,采用混合粒子群算法进行快速路径寻优。本发明提出的基于混合粒子群算法的星图识别方法解决了大视场,高灵敏度恒星敏感器条件下,星图识别率低,对噪声鲁棒性差的缺点;实现了导航星的快速匹配识别,提高了识别率,对噪声的鲁棒性更好。
Description
技术领域
本发明涉及一种基于混合粒子群算法的星图识别方法,属于航天器导航、制导与控制领域。
背景技术
星敏感器是基于星图匹配识别技术的天文导航***的主要工作部件,星图匹配识别方法是星敏感器技术的核心,可以快速准确地获得载体的姿态信息。因此,研究一种识别速度快,识别成功率高的星图识别算法具有重要的理论和现实意义。
目前已有了许多较为成熟的星图匹配识别方法,如三角形算法,凸多边形匹配识别算法、栅格算法、基于Delaunay三角剖分识别法、基于神经网络的星图识别算法、基于遗传算法的星图识别方法以及基于蚁群算法的星图识别方法等。这些星图识别算法都达到了较为良好的效果,但它们在识别速度、识别成功率、数据库大小以及对噪声的鲁棒性上各有优劣;特别在大视场、高灵敏度的恒星敏感器条件下,由于光学***畸变大、敏感星点多、星点质心提取困难等原因,识别速度和成功率更是明显降低。
发明内容
本发明的目的是解决大视场条件下恒星敏感器的星图识别速度慢和识别成功率低的问题。本发明提出了一种基于混合粒子群算法的星图识别方法,显著地提高了星图识别的成功率和对噪声的鲁棒性。
本发明提出的一种基于混合粒子群算法的星图识别方法,包括以下几个步骤:
步骤一:建立导航星库,实现星图识别算法导航星库的构造;
步骤二:星图匹配识别,利用混合粒子群算法为导航星构造特征模式,然后采用二分查找法对导航星进行匹配识别。
本发明的优点在于:
(1)本发明公开的星图识别方法采用混合粒子群算法对星点进行路径寻优,根据最优路径长度值进行索引,并利用最优路径中导航星间的角距进行匹配识别,提高了方法的识别率,减少了匹配时间;
(2)本发明公开的星图识别方法提出了一种自适应确定圆半径r的方法,使该识别方法的寻优规模在合理范围内,保证了方法的实时性。
(3)本发明公开的星图识别算法加入了导航星验证环节,提高了识别结果的准确性,同时减少了冗余匹配。
附图说明
图1是本发明的方法流程图;
图2是本发明圆半径r取4°时,圆内星点数目不同时的导航星个数分布;
图3是本发明对r=4°圆内星点少于6的导航星,取半径r=5°时圆内星点个数分布;
图4是本发明对r=4°圆内星点多于25的导航星,取半径r=2.5°时圆内星点个数分布;
图5是本发明提出的基于混合粒子群算法的圆内星点路径寻优图;
图6是本发明提出的基于混合粒子群算法的星图识别方法流程;
图7是本发明提出的方法识别率与星点位置误差的关系图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是一种基于混合粒子群算法的星图识别方法,方法流程如图1所示,包括以下几个步骤:
步骤一:建立导航星库,实现星图识别方法导航星库的构造
针对基于混合粒子群算法的星图识别方法,其导航星库的构造方法如下:
第一步:选取导航星,构建基本导航星表
从基本星表中选取导航星,对恒星间角距小于某一值(一般取3-5个像素)的双星进行剔除,得到用于星图识别的基本导航星表,基本导航星表包含每一颗导航星的编号、星等、赤经和赤纬信息。
第二步:构造每颗导航星的特征数据集合
针对基本导航星表中的每颗导航星,以导航星为圆心,以半径r进行画圆,将圆内的所有导航星组成圆心导航星的特征数据集合。
星图识别是一种特征匹配方法,需要为每颗导航星构造唯一的特征模式,然后利用构造的特征实现导航星的匹配识别。基于混合粒子群算法的星图识别方法在构造特征模式时,需要给定圆半径r,而圆半径r的大小影响圆内导航星的个数,进而决定混合粒子群算法路径寻优的难易程度。为了提高星图识别方法的准确性和实时性,本发明提出了一种自适应调整方法来确定圆半径r,以保证混合粒子群算法的寻优规模在合理范围内,具体步骤包括:
(1)以一颗导航星为中心,取圆半径r=4°,统计圆内的星点个数,当星点个数小于6时,进入步骤(2),当星点个数大于25时,进入步骤(3),当星点个数大于等于6且小于等于25时,确定圆半径r=4°;
为保证星图识别方法的准确性和实时性,本发明对以上两种情况的圆半径r进行调整;
(2)当圆内的星点个数小于6时,增大圆半径r为5°;
(3)当圆内的星点个数大于25时,减小圆半径r为2.5°。
第三步:选取基本导航星表中的一颗导航星,获取特征数据集合中两两恒星间的角距,然后以特征数据集合中圆心处的导航星为起点,利用混合粒子群算法进行快速路径寻优,得到特征数据集合的最优路径,获取最优路径的长度,并选择最优路径中离圆心最近的导航星为最优路径的前进方向,得到最优路径的前三颗导航星。
所述混合粒子群算法的快速路径寻优步骤如下:
(1)假设特征数据集合中共有m颗导航星,首先对特征数据集合中的导航星进行编码,将粒子遍历所有导航星的顺序以自然数排列,作为粒子的解;然后初始化混合粒子群算法的参数,设定粒子个数为n,迭代次数最大值为Nmax;然后随机产生n个初始解(初始路径),并将粒子遍历所有导航星的路径长度之和作为粒子的适应度函数,然后根据粒子的当前位置计算起适应度值,设定各个粒子的个体极值,以及所有粒子的全局极值。
(2)将第i个粒子的当前解与个体极值进行交叉,具体交叉方法为:对于粒子的当前解,随机选择一个交叉区域(解中的一段数据),然后删除解中与个体极值交叉区域相同的元素,然后将个体极值交叉区域中的元素***粒子的当前解,得到第i个粒子的新解。
(3)将步骤(2)得到的新解与全局极值进行交叉,具体交叉方法为:对于粒子的当前解,随机选择一个交叉区域,删除解中与全局极值交叉区域相同的元素,然后将全局极值交叉区域中的元素***粒子的当前解,得到第i个粒子的新解。
(4)将步骤(3)得到的新解进行变异,为使路径长度之和达到最小,采用轮盘赌选择法,以较大的概率(概率大的导航星更易被选择)选取路径中相邻导航星间角距大的两个导航星,然后在它们之间***其它导航星,得到变异后的新解。
具体方法为:获取步骤(3)得到的新解中相邻导航星间的角距,分别记为l(k),k=1,2,…,m,其中,l(1)表示解中第一颗与第二颗导航星间的角距,以此类推,l(m-1)表示解中第m-1颗与第m颗导航星间的角距,l(m)表示解中第m颗与第一颗导航星间的角距,则每颗导航星被选择的概率为:
根据上式求得的概率,采用轮盘赌选择法选取导航星i1。采用轮盘赌选择法时,解中的每颗导航星类似于轮盘中的一小块扇形,扇形的面积大小与该导航星被选择的概率成正比,扇形面积越大,则该导航星被选择的概率也越大。
选取了导航星i1后,以较大的概率选取离星点i1较近的星点作为粒子的下一个遍历点。设d(i1,j)表示解中导航星i1与导航星j间的距离,则可计算得到离导航星i1最远导航星的距离为dmax,其中k=1,…,m且k≠i1,为防止遍历导航星i1后的下一个导航星为i1本身,令d(i1,i1)=dmax,则每颗导航星被选为下一个遍历导航星的概率为:
根据上式求得的概率,采用轮盘赌选择法选取星点j1,即概率大的星点更易被选择,然后把导航星j1安排到导航星i1之后,其余不变,得到变异后的新解。
(5)进行个体极值的更新,具体方法为:首先计算第i个粒子变异后新解的适应度值,并与该粒子个体极值的适应度值相比,得到第i个粒子的适应度值变化量,然后采用模拟退火机制确定是否接受该新解,如果接受,则更新粒子的个体极值。
判断所有粒子的个体极值是否都已更新,若都已更新,进入步骤(6),否则返回步骤(2)。
(6)进行全局极值的更新,根据所有粒子的个体极值,确定全局极值和全局极值路径,判断迭代次数是否大于迭代次数最大值Nmax,如果大于,则路径寻优结束,全局极值路径即为最优路径,否则,返回的步骤(2)。
第四步:获取最优路径中前三颗导航星的信息,即分别获取第一颗导航星、第二颗导航星和第三颗导航星所对应的赤经、赤纬信息以及星等,第一颗导航星与第二颗导航星间、第二颗导航星与第三颗导航星间、第一颗导航星与第三颗导航星间的角距信息。
第一颗导航星与第二颗导航星间、第二颗导航星与第三颗导航星间的角距信息,第一颗导航星、第二颗导航星和第三颗导航星分别对应的赤经、赤纬信息,作为导航星的基本识别信息。
第一颗导航星与第三颗导航星间的角距以及第一颗导航星、第二颗导航星和第三颗导航星的星等,作为星图识别的验证信息。
第五步:重复第三步与第四步,获取导航星表中每一颗导航星的特征信息,特征信息包括导航星的最优路径长度、导航星基本识别信息以及星图识别的验证信息。对于导航星半径r=5°的导航星识别信息按照最优路径长度升序排列,构成导航星库的第一部分,对于半径r=4°的导航星识别信息按照最优路径长度升序排列,构成导航星库的第二部分,对于半径r=2.5°的导航星识别信息按照最优路径长度升序排列,构成导航星库的第三部分,最后这三部分构成一个完整的导航星库。
步骤二:基于混合粒子群算法进行星图识别;
对于一幅星图,经过星图预处理、星体提取,得到星图中每个星点的位置坐标及灰度信息,然后进行星图匹配识别,如图6所示,基于混合粒子群算法的星图识别方法具体包括以下几个步骤:
第一步:确定候选识别主星,选取的识别主星应具有亮度大且离视场边缘角距均大于圆半径r的特点,选取方法如下:计算星图中所有离视场边缘角距均大于圆半径r的星点的平均灰度值,然后选取灰度值大于平均灰度值的星点作为候选识别主星,得到候选识别主星集合。
用户设置识别失败次数的最大值为N,设方法的识别次数为n,n的初值设为0。
第二步:进行星点快速路径寻优,在候选识别主星集合中,首先选取离视场中心第n+1近的星点作为圆心(即第一次识别时,选择最近的星点,第二次识别时,选择次近的星点,依次类推),然后以半径r进行画圆,半径r的确定方法与步骤一中第二步所述方法相同,将圆内的所有星点构成特征数据集合,如图5所示,其中白色星点表示该星点的灰度值大于平均灰度值。
计算特征数据集合中两两星点间的角距,然后以特征数据集合中圆心处的星点为起点,利用混合粒子群算法进行快速路径寻优,得到星点特征数据集合的最优路径,获取最优路径的长度,选择最优路径中离圆心最近的星点为最优路径的前进方向,得到最优路径中的前三个星点。
其中,利用混合粒子群算法进行快速路径寻优与步骤一第三步中所述混合粒子群算法的快速路径寻优方法相同,将星点的特征数据集合替代导航星的数据特征集合。
获取最优路径中前三个星点的信息,即分别获取第一颗星点、第二颗星点和第三颗星点所对应的星等,第一颗星点与第二颗星点间的角距d12、第二颗星点与第三颗星点间的角距d23、第一颗星点与第三颗星点间的角距d13。
第三步:进行星图匹配识别,根据星点的半径r,确定该星点在导航星库中对应半径r的部分,根据第二步得到的最优路径长度,采用二分查找法,查找导航星库中与该星点最优路径长度相近的导航星,并把它们作为候选识别星,然后,将第二步中得到的d12和d23与候选识别星的基本识别信息进行匹配,得到匹配结果。匹配结果为无导航星匹配、仅一颗导航星匹配或者存在多颗导航星匹配。
第四步:导航星识别结果验证,将第三步得到的匹配结果进行验证,具体验证方法如下:
(1)若仅一颗导航星匹配上,则利用第二步得到的角距d13与导航星的星图识别验证信息进行匹配,如果匹配上,则此幅星图识别成功,如果没有匹配上,则此次星图匹配失败,星图识别失败次数n加1;
(2)若没有导航星匹配上,则此次星图匹配失败,星图识别失败次数n加1;
(3)若存在多颗导航星匹配上,则利用第二步得到的角距d13与多颗导航星的星图识别验证信息进行匹配,得到进一步的匹配结果。如果无导航星匹配,则此次星图匹配失败,星图识别失败次数n加1;如果仅一颗导航星匹配,则此幅星图识别成功;如果仍存在多颗导航星匹配,则利用第二步最优路径中第一颗星点、第二颗星点和第三颗星点分别对应的星等与星图识别的验证信息进行匹配,若无导航星匹配或仍有多颗导航星匹配,则认为此次星图匹配失败,星图识别失败次数n加1,否则此幅星图识别成功。
第五步:当失败次数n小于识别失败次数最大值N时,返回步骤第二步,否则,此幅星图匹配失败,此幅星图识别结束。
实施例:
在选择导航星,构造导航星表时,选取导航星的最大星等为6.5,剔除恒星间角距小于某一值(取3个像素)的双星,最后得到的基本导航星表共有8996颗导航星。在确定圆半径r时,取圆半径r=4°,圆内星点数目取不同值时的导航星个数近似服从正态分布,圆内星点数小于6的导航星共有404颗,大于25的导航星共有628颗,如图2所示。对于圆内导航星数小于6时,增大圆半径r为5°,圆内不同星点个数的导航星分布如图3所示;对于圆内导航星数大于25时,减小圆半径r为2.5°,圆内不同星点个数的导航星分布如图4所示。
基于混合粒子群算法的星图识别方法,首先根据步骤一构造用于星图识别的导航星库,构造的导航星库由最优路径长度表、角距识别星表以及识别验证星表组成。最优路径长度表存储了最优路径的长度值,角距识别星表包含最优路径中第一颗导航星与第二颗导航星间的角距、第二颗导航星与第三颗导航星间的角距以及该三颗导航星所对应的赤经赤纬信息,识别验证星表包含最优路径中第一颗导航星与第三颗导航星间的角距以及最优路径中前三颗导航星的星等,导航星库中三个表的存储内容分别如表1、表2和表3所示。然后根据步骤二实现导航星的匹配识别,用户自行设定识别失败次数最大值(一般取4以内的正整数),方法的识别流程如图6所示。
表1最优路径长度表
最优路径值(°) |
... |
11.030192 |
11.453998 |
11.607923 |
12.864802 |
... |
表2星图角距识别表
角距1(°) | 角距2(°) | 赤经1(°) | 赤纬1(°) | 赤经2(°) | 赤纬2(°) | 赤经3(°) | 赤纬3(°) |
... | ... | ... | ... | ... | ... | ... | ... |
2.682591 | 3.831108 | 48.137904 | -57.321636 | 45.903517 | -59.737742 | 53.214933 | -61.016722 |
2.761853 | 5.601960 | 47.819796 | -16.025144 | 47.771804 | -13.263681 | 49.671696 | -18.559589 |
1.184397 | 1.059622 | 213.017000 | 69.432614 | 210.460967 | 68.678717 | 207.745563 | 68.314947 |
2.760843 | 1.592563 | 144.260825 | 16.437953 | 141.385337 | 16.585803 | 140.314033 | 15.371067 |
... | ... | ... | ... | ... | ... | ... | ... |
表3识别验证表
角距(°) | 星等1 | 星等2 | 星等3 |
... | ... | ... | ... |
4.516492 | 5.7 | 5.2 | 6.3 |
3.090185 | 6.3 | 6.5 | 5.8 |
2.203262 | 5.4 | 6.4 | 6.4 |
3.942676 | 5.9 | 6.3 | 6.5 |
... | ... | ... | ... |
应用本发明提出的基于混合粒子群算法的星图识别方法,可以显著提高大视场、高灵敏度恒星敏感器条件下方法的识别成功率和对噪声的鲁棒性。采用本发明的基于混合粒子群算法的星图识别方法、传统三角形算法以及改进的神经网络星图识别算法,分别对随机生成的1000幅星图进行识别,其中基于混合粒子群算法的星图识别方法选取的识别失败次数最大值为1,即仅进行一次路径寻优。进行识别时,三种方法均基于以下条件:星敏感器的视场大小为20°×20°,星敏感器面阵大小为1024×1024,像元尺寸为15μm×15μm。进行仿真时,星点位置误差取均值为0,方差为0至2.5个像素,对三种星图识别方法在不同位置误差条件下均仿真1000次,方法的识别率如图7所示。
从图7可以看出,在星点位置误差的方差小于2个像素的情况下,本发明的星图识别方法识别成功率最高,改进的神经网络算法识别率次之,传统三角形算法的识别率最低。当星点位置误差方差大于2个像素时,基于本发明的识别率开始小于改进的神经网络识别方法,但优于传统三角形算法。
对于本发明所述方法的识别速度,在Intel Core i3处理器2.39GHz PC机上,采用VS2010进行测试,本发明的方法从特征模式的构造到星图的匹配识别,所需的平均时间为0.15s。由于导航星库中的特征数据是由混合粒子群算法路径寻优得到的最优路径值,采用升序方式进行存储,因而可以利用二分查找法进行搜索与匹配。其检索速度比传统三角形算法快,与改进的神经网络星图识别算法相当,但三种识别方法构造特征的速度从小到大依次是传统三角形算法、改进的神经网络星图识别方法以及基于粒子群的星图识别方法。
Claims (4)
1.一种基于混合粒子群算法的星图识别方法,其特征在于,包括以下几个步骤:
步骤一:建立导航星库
具体包括以下几个步骤:
第一步:选取导航星,构建基本导航星表
从基本星表中选取导航星,对恒星间角距小于3~5个像素的双星进行剔除,得到用于星图识别的基本导航星表,基本导航星表包含每颗导航星的编号、星等、赤经和赤纬信息;
第二步:构造每颗导航星的特征数据集合
针对基本导航星表中的每颗导航星,分别以每颗导航星为圆心,以半径r进行画圆,将圆内的所有导航星组成圆心导航星的特征数据集合;
第三步:选取基本导航星表中的一颗导航星,获取特征数据集合中两两恒星间的角距,然后以特征数据集合中圆心处的导航星为起点,利用混合粒子群算法进行快速路径寻优,得到特征数据集合的最优路径,获取最优路径的长度,并选择最优路径中离圆心最近的导航星为最优路径的前进方向,得到最优路径中的前三颗导航星;
第四步:获取最优路径中前三颗导航星的信息,即分别获取第一颗导航星、第二颗导航星和第三颗导航星所对应的赤经、赤纬信息以及星等,第一颗导航星与第二颗导航星间、第二颗导航星与第三颗导航星间、第一颗导航星与第三颗导航星间的角距信息;
第一颗导航星与第二颗导航星间、第二颗导航星与第三颗导航星间的角距信息,第一颗导航星、第二颗导航星和第三颗导航星分别对应的赤经、赤纬信息,作为导航星的基本识别信息;
第一颗导航星与第三颗导航星间的角距以及第一颗导航星、第二颗导航星和第三颗导航星的星等,作为星图识别的验证信息;
第五步:重复第三步与第四步,获取导航星表中每一颗导航星的特征信息,特征信息包括导航星的最优路径长度、导航星的基本识别信息以及星图识别的验证信息,每颗导航星的特征信息按照最优路径长度的升序方式排列,并构成导航星库;
步骤二:基于混合粒子群算法进行星图识别
对于一幅星图,经过星图预处理、星体提取,得到星图中每个星点的位置坐标及灰度信息,然后进行星图匹配识别,基于混合粒子群算法的星图识别方法具体步骤如下:
第一步:确定候选识别主星,选取的识别主星应具有亮度大且离视场边缘角距均大于圆半径r的特点,选取方法如下:计算星图中所有离视场边缘角距均大于圆半径r的星点的平均灰度值,然后选取灰度值大于平均灰度值的星点作为候选识别主星,得到候选识别主星集合;
用户设置识别失败次数的最大值为N,设识别次数为n,n的初值设为0;
第二步:进行星点快速路径寻优,在候选识别主星集合中,首先选取离视场中心第n+1近的星点作为圆心,然后以半径r进行画圆,将圆内的所有星点构成特征数据集合;
计算特征数据集合中两两星点间的角距,然后以特征数据集合中圆心处的星点为起点,利用混合粒子群算法进行快速路径寻优,得到星点特征数据集合的最优路径,获取最优路径的长度,选择最优路径中离圆心最近的星点为最优路径的前进方向,得到最优路径中的前三个星点;
获取最优路径中前三个星点的信息,即分别获取第一颗星点、第二颗星点和第三颗星点所对应的星等,第一颗星点与第二颗星点间的角距d12、第二颗星点与第三颗星点间的角距d23、第一颗星点与第三颗星点间的角距d13;
第三步:进行星图匹配识别,根据星点的半径r,确定该星点在导航星库中对应半径r的部分,根据第二步得到的最优路径长度,采用二分查找法,查找导航星库中与该星点最优路径长度相近的导航星,并把它们作为候选识别星,然后利用第二步中得到的d12和d23与候选识别星的基本识别信息进行匹配,得到匹配结果;匹配结果为无导航星匹配、仅一颗导航星匹配或者存在多颗导航星匹配;
第四步:导航星识别结果验证,将第三步得到的匹配结果进行验证,具体验证方法如下:
(1)若仅一颗导航星匹配上,则利用第二步得到的角距d13与导航星的星图识别验证信息进行匹配,如果匹配上,则此幅星图识别成功,如果没有匹配上,则此次星图匹配失败,星图识别失败次数n加1;
(2)若没有导航星匹配上,则此次星图匹配失败,星图识别失败次数n加1;
(3)若存在多颗导航星匹配上,则利用第二步得到的角距d13与多颗导航星的星图识别验证信息进行匹配,得到进一步的匹配结果;如果无导航星匹配,则此次星图匹配失败,星图识别失败次数n加1;如果仅一颗导航星匹配,则星图识别成功;如果仍存在多颗导航星匹配,则利用第二步最优路径中第一颗星点、第二颗星点和第三颗星点分别对应的星等与星图识别的验证信息进行匹配,若无导航星匹配或仍有多颗导航星匹配,则认为此次星图匹配失败,星图识别失败次数n加1,否则此幅星图识别成功;
第五步:当失败次数n小于识别失败次数最大值N时,返回步骤第二步,否则,此幅星图匹配失败,此幅星图识别结束。
2.根据权利要求1所述的一种基于混合粒子群算法的星图识别方法,其特征在于,步骤一第二步中、步骤二第一步中,所述圆半径r的选取方法具体为:
(1)以一颗导航星为中心,取圆半径r=4°,统计圆内的星点个数,当星点个数小于6时,进入步骤(2),当星点个数大于25时,进入步骤(3),当星点个数大于等于6且小于等于25时,确定圆半径r=4°;
(2)当圆内的星点个数小于6时,增大圆半径r为5°;
(3)当圆内星点个数大于25时,减小圆半径r为2.5°。
3.根据权利要求1所述的一种基于混合粒子群算法的星图识别方法,其特征在于,步骤一第三步中、步骤二第二步中,所述的混合粒子群算法的快速路径寻优步骤如下:
(1)假设特征数据集合中共有m颗导航星,首先对特征数据集合中的导航星进行编码,将粒子遍历所有导航星的顺序以自然数排列,作为粒子的解;然后初始化混合粒子群算法的参数,设定粒子个数为n,迭代次数最大值为Nmax;然后随机产生n个初始解,并将粒子遍历所有导航星的路径长度之和作为粒子的适应度函数;然后根据粒子的当前位置计算其适应度值,设定每个粒子的个体极值,以及所有粒子的全局极值;
(2)将第i个粒子的当前解与个体极值进行交叉,具体交叉方法为:对于粒子的当前解,随机选择一个交叉区域,然后删除解中与个体极值交叉区域相同的元素,然后将个体极值交叉区域中的元素***粒子的当前解,得到第i个粒子的新解;
(3)将步骤(2)得到的新解与全局极值进行交叉,具体交叉方法为:对于粒子的当前解,随机选择一个交叉区域,删除解中与全局极值交叉区域相同的元素,然后将全局极值交叉区域中的元素***粒子的当前解,得到第i个粒子的新解;
(4)将步骤(3)得到的新解进行变异,为使路径长度之和达到最小,采用轮盘赌选择法,以较大的概率选取路径中相邻导航星间角距较大的两个导航星,然后在它们之间***其它导航星,得到变异后的新解;
具体方法为:获取步骤(3)得到的新解中相邻导航星间的角距,分别记为l(k),k=1,2,…,m,其中,l(1)表示解中第一颗与第二颗导航星间的角距,以此类推,l(m-1)表示解中第m-1颗与第m颗导航星间的角距,l(m)表示解中第m颗与第一颗导航星间的角距,则每颗导航星被选择的概率为:
根据上式求得的概率,采用轮盘赌选择法选取导航星i1;采用轮盘赌选择法时,解中的每颗导航星类似于轮盘中的一小块扇形,扇形的面积大小与该导航星被选择的概率成正比,扇形面积越大,则该导航星被选为下一个遍历导航星的概率也越大;
选取了导航星i1后,以较大的概率选取离星点i1较近的星点作为粒子的下一个遍历点;设d(i1,j)表示解中导航星i1与导航星j间的距离,则可计算得到离导航星i1最远导航星的距离为dmax,其中k=1,…,m且k≠i1,为防止遍历导航星i1后的下一个导航星为i1本身,令d(i1,i1)=dmax,则每颗导航星被选为下一个遍历导航星的概率为:
根据上式求得的概率,采用轮盘赌选择法选取星点j1,即概率大的星点更易被选择,然后把导航星j1安排到导航星i1之后,其余不变,得到变异后的新解;
(5)进行个体极值的更新,具体方法为:首先计算第i个粒子变异后新解的适应度值,并与该粒子个体极值的适应度值相比,得到第i个粒子的适应度值变化量,然后采用模拟退火机制确定是否接受该新解,如果接受,则更新该粒子的个体极值;
判断所有粒子的个体极值是否都已更新,若都已更新,进入步骤(6),否则返回步骤(2);
(6)进行全局极值的更新,根据所有粒子的个体极值,确定全局极值和全局极值路径,判断迭代次数是否大于迭代次数最大值Nmax,如果大于,则路径寻优结束,全局极值路径即为最优路径,否则,返回的步骤(2)。
4.根据权利要求1所述的一种基于混合粒子群算法的星图识别方法,其特征在于,步骤一第五步中,所述的导航星库构成方法为:对于半径r=5°的导航星识别信息按照最优路径长度升序排列,构成导航星库的第一部分;对于半径r=4°的导航星识别信息按照最优路径长度升序排列,构成导航星库的第二部分;对于半径r=2.5°的导航星识别信息按照最优路径长度升序排列,构成导航星库的第三部分,最后这三部分构成一个完整的导航星库。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210315184.1A CN102840860B (zh) | 2012-08-30 | 2012-08-30 | 一种基于混合粒子群算法的星图识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210315184.1A CN102840860B (zh) | 2012-08-30 | 2012-08-30 | 一种基于混合粒子群算法的星图识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102840860A true CN102840860A (zh) | 2012-12-26 |
CN102840860B CN102840860B (zh) | 2016-05-11 |
Family
ID=47368439
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210315184.1A Expired - Fee Related CN102840860B (zh) | 2012-08-30 | 2012-08-30 | 一种基于混合粒子群算法的星图识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102840860B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217163A (zh) * | 2013-04-02 | 2013-07-24 | 北京航空航天大学 | 一种混合冗余特征的星模式识别方法及装置 |
CN105258697A (zh) * | 2015-10-21 | 2016-01-20 | 北京航空航天大学 | 一种星光空白段神经网络连续导航方法及*** |
CN106403929A (zh) * | 2016-08-22 | 2017-02-15 | 清华大学 | 星敏感器星图识别和姿态解算的快速鲁棒校验方法 |
CN106595645A (zh) * | 2016-11-02 | 2017-04-26 | 上海航天控制技术研究所 | 一种基于星敏感器输出精度的导航星库制作方法 |
CN107816987A (zh) * | 2017-09-12 | 2018-03-20 | 北京航空航天大学 | 一种基于蜘蛛网模式和卷积神经网络的星图识别方法 |
CN107883947A (zh) * | 2017-12-28 | 2018-04-06 | 常州工学院 | 基于卷积神经网络的星敏感器星图识别方法 |
CN112489011A (zh) * | 2020-11-27 | 2021-03-12 | 上海航天控制技术研究所 | 一种星敏感器光机组件智能装调方法 |
CN112595312A (zh) * | 2020-11-09 | 2021-04-02 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种大视场星敏感器伪星目标滤除方法及*** |
CN114980136A (zh) * | 2022-05-20 | 2022-08-30 | 西安电子科技大学 | 一种高能效的地面基站对低空立体信号覆盖方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102072730A (zh) * | 2010-11-04 | 2011-05-25 | 北京航空航天大学 | 一种三角形星图识别方法 |
CN102128623A (zh) * | 2010-12-02 | 2011-07-20 | 哈尔滨工业大学 | 一种星敏感器快速星图识别方法 |
CN102313549A (zh) * | 2011-07-28 | 2012-01-11 | 北京航空航天大学 | 一种基于惯性比特征的三角形星图识别方法 |
-
2012
- 2012-08-30 CN CN201210315184.1A patent/CN102840860B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102072730A (zh) * | 2010-11-04 | 2011-05-25 | 北京航空航天大学 | 一种三角形星图识别方法 |
CN102128623A (zh) * | 2010-12-02 | 2011-07-20 | 哈尔滨工业大学 | 一种星敏感器快速星图识别方法 |
CN102313549A (zh) * | 2011-07-28 | 2012-01-11 | 北京航空航天大学 | 一种基于惯性比特征的三角形星图识别方法 |
Non-Patent Citations (1)
Title |
---|
全伟等: "《一种基于蚁群聚类算法的快速星图识别方法》", 《宇航学报》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103217163B (zh) * | 2013-04-02 | 2015-07-15 | 北京航空航天大学 | 一种混合冗余特征的星模式识别方法及装置 |
CN103217163A (zh) * | 2013-04-02 | 2013-07-24 | 北京航空航天大学 | 一种混合冗余特征的星模式识别方法及装置 |
CN105258697A (zh) * | 2015-10-21 | 2016-01-20 | 北京航空航天大学 | 一种星光空白段神经网络连续导航方法及*** |
CN105258697B (zh) * | 2015-10-21 | 2018-01-12 | 北京航空航天大学 | 一种星光空白段神经网络连续导航方法及*** |
CN106403929B (zh) * | 2016-08-22 | 2019-02-05 | 清华大学 | 星敏感器星图识别和姿态解算的快速鲁棒校验方法 |
CN106403929A (zh) * | 2016-08-22 | 2017-02-15 | 清华大学 | 星敏感器星图识别和姿态解算的快速鲁棒校验方法 |
CN106595645A (zh) * | 2016-11-02 | 2017-04-26 | 上海航天控制技术研究所 | 一种基于星敏感器输出精度的导航星库制作方法 |
CN107816987A (zh) * | 2017-09-12 | 2018-03-20 | 北京航空航天大学 | 一种基于蜘蛛网模式和卷积神经网络的星图识别方法 |
CN107816987B (zh) * | 2017-09-12 | 2020-07-31 | 北京航空航天大学 | 一种基于蜘蛛网模式和卷积神经网络的星图识别方法 |
CN107883947A (zh) * | 2017-12-28 | 2018-04-06 | 常州工学院 | 基于卷积神经网络的星敏感器星图识别方法 |
CN107883947B (zh) * | 2017-12-28 | 2020-12-22 | 常州工学院 | 基于卷积神经网络的星敏感器星图识别方法 |
CN112595312A (zh) * | 2020-11-09 | 2021-04-02 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种大视场星敏感器伪星目标滤除方法及*** |
CN112595312B (zh) * | 2020-11-09 | 2022-07-05 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种大视场星敏感器伪星目标滤除方法及*** |
CN112489011A (zh) * | 2020-11-27 | 2021-03-12 | 上海航天控制技术研究所 | 一种星敏感器光机组件智能装调方法 |
CN112489011B (zh) * | 2020-11-27 | 2023-01-31 | 上海航天控制技术研究所 | 一种星敏感器光机组件智能装调方法 |
CN114980136A (zh) * | 2022-05-20 | 2022-08-30 | 西安电子科技大学 | 一种高能效的地面基站对低空立体信号覆盖方法 |
CN114980136B (zh) * | 2022-05-20 | 2023-06-30 | 西安电子科技大学 | 一种高能效的地面基站对低空立体信号覆盖方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102840860B (zh) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102840860A (zh) | 一种基于混合粒子群算法的星图识别方法 | |
CN109145939B (zh) | 一种小目标敏感的双通道卷积神经网络语义分割方法 | |
CN102693542B (zh) | 一种影像特征匹配方法 | |
CN103295239B (zh) | 一种基于平面基准影像的激光点云数据的自动配准方法 | |
CN106599830A (zh) | 人脸关键点定位方法及装置 | |
CN103994765B (zh) | 一种惯性传感器的定位方法 | |
CN111681178B (zh) | 一种基于知识蒸馏的图像去雾方法 | |
CN110717927A (zh) | 基于深度学习和视惯融合的室内机器人运动估计方法 | |
CN102682302A (zh) | 一种基于关键帧的多特征融合的人体姿态识别方法 | |
CN103453905B (zh) | 一种基于平面三角形主成分分析的星图识别方法 | |
CN104933755A (zh) | 一种静态物体重建方法和*** | |
CN103954280B (zh) | 一种快速和高鲁棒性自主恒星识别方法 | |
CN103868510A (zh) | 一种快速全天星图自主恒星识别方法 | |
CN104504709A (zh) | 一种基于特征球的室外场景三维点云数据的分类方法 | |
CN104408760A (zh) | 一种基于双目视觉的高精度虚拟装配***算法 | |
CN103884281B (zh) | 一种基于主动结构光的巡视器障碍探测方法 | |
CN108447048B (zh) | 基于关注层的卷积神经网络图像特征处理方法 | |
CN109711444A (zh) | 一种新的基于深度学习的遥感影像配准方法 | |
CN107393004A (zh) | 一种获取输电线走廊中建筑物拆迁量的方法及装置 | |
CN106595702A (zh) | 一种基于天文标定的多传感器空间配准方法 | |
CN104778465A (zh) | 一种基于特征点匹配的目标跟踪方法 | |
CN103363987A (zh) | 一种多视场星敏感器的星图识别方法 | |
CN102663740A (zh) | 基于图切的sar图像变化检测方法 | |
CN109993108B (zh) | 一种增强现实环境下的手势纠错方法、***及装置 | |
CN109598261B (zh) | 一种基于区域分割的三维人脸识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160511 Termination date: 20170830 |
|
CF01 | Termination of patent right due to non-payment of annual fee |