CN102944239B - 一种用于深空探测自主导航的目标天体图像晨昏线判断方法 - Google Patents
一种用于深空探测自主导航的目标天体图像晨昏线判断方法 Download PDFInfo
- Publication number
- CN102944239B CN102944239B CN201210409052.5A CN201210409052A CN102944239B CN 102944239 B CN102944239 B CN 102944239B CN 201210409052 A CN201210409052 A CN 201210409052A CN 102944239 B CN102944239 B CN 102944239B
- Authority
- CN
- China
- Prior art keywords
- segmental arc
- terminator
- candidate
- celestial body
- segmental
- 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.)
- Active
Links
Landscapes
- Image Analysis (AREA)
Abstract
一种用于深空探测自主导航的目标天体图像晨昏线判断方法,目标天体图像的边缘提取,使用传统边缘算子法进行边缘点提取;候选弧段选择;拟合判断;晨昏线补充判断;本发明中的空间目标的晨昏线若不进行判断并剔除则将直接影响目标天体图像中心位置的拟合精度,本发明解决了行星际探测中目标天体晨昏线判断的问题,通过晨昏线判断找到真正的边缘点进行目标天体中心位置拟合,提升了天体图像目标中心的定位精度;具有高鲁棒性、高刷新率、高精度导航天体敏感器的研制具有重要工程使用价值。
Description
技术领域
本发明涉及一种用于深空探测自主导航的目标天体图像晨昏线判断方法,属于空间光学成像敏感器领域。
背景技术
在奔火、环火弧段均需要对火星等其他深空探测目标进行拍照并确定火星矢量指向。由于太阳光照、轨道位置等制约火星图像将出现类似月球的晨昏线,此时用于火星中心矢量拟合的点若使用了晨昏线点将直接影响到定位精度,必须在拟合过程中判断出晨昏线点并且进行排除,因此必须设计晨昏线判断方法。文献1,李正,晨昏线的判定,中学政史地:教学指导版,2011(5),对晨昏线的概念、特点、分类、运动规律进行了描述和说明,但未对晨昏线成像模型及识别方法进行研究;文献2,杨艳华,巧用晨昏线判断太阳直射点,新课程(中学版),2010(3),针对地球晨昏线的位置特点,研究了太阳直射地球纬度的确定方法。
上述文献均未解决从成像角度进行晨昏线的识别与判断问题,无法用于深空探测成像敏感器目标天体中心矢量的测量。
发明内容
本发明的技术解决问题是:克服现有技术的不足之处,提供一种用于深空探测自主导航的目标天体图像晨昏线判断方法,解决了深空探测目标天体中心矢量高精度测量问题。
本发明技术方案:一种用于深空探测自主导航的目标天体图像晨昏线判断方法,实现步骤为:
(1)目标天体图像的边缘提取,使用传统边缘算子法进行边缘点提取;
(2)计算天体图像360°所有边缘点的梯度大小;
(3)从0°开始以90°弧段为单位,每间隔30°进行一次弧段选择,并计算该90°弧段内所含有边缘点梯度大小的均值,对所有弧段进行梯度均值排序,选择均值最小的弧段作为晨昏线候选弧段,同时将晨昏候选弧段正的90°弧段选择为边缘候选弧段;
(4)对上述晨昏候选弧段、边缘候选弧段分别进行圆拟合,若晨昏候选弧段拟合中心与边缘候选弧段拟合中心位置的偏差小于等于设定阈值(一般选择为0.2~3个像元),则说明目标天体成像中没有晨昏线,将目标天体图像360°所有边缘点进行重新拟合得到最终的目标天体图像中心位置;若晨昏候选弧段拟合中心与边缘候选弧段拟合中心位置的偏差大于设定阈值(一般选择为0.2~3个像元),则进行所述圆拟合误差比较,若晨昏线候选弧段圆拟合误差大于边缘候选弧段圆拟合误差,则说明目标天体图像存在晨昏线且位置在候选弧段附近,此时将边缘候选弧段在起始端点、终止端点分别延伸30°将原来90°弧段延长形成150°的边缘点重新进行拟合得到最终的目标天体图像中心位置;若晨昏线候选弧段圆拟合误差小于等于边缘候选弧段圆拟合误差,则转入步骤(5)进行晨昏线补充判断;
(5)晨昏线补充判断,按照步骤(3)的弧段划分方法,以90°弧段为单位,进行所有弧段进行圆拟合并记录拟合误差,拟合输出中心位置点坐标(Xi,Yi),对XY坐标分别进行排序选择中间值(Xm,Ym)设为O点,距离O点最远的输出中心位置点(Xj,Yj)对应的拟合弧段为晨昏候选弧段,选择其正对的弧段为边缘候选弧段,将边缘候选弧段在起始端、终止端分别扩张30°形成150°的边缘点重新进行拟合得到最终的目标天体图像中心位置。
所述步骤(2)中计算天体图像360°所有边缘点的梯度大小如下:
Gradx(x,y)=f(x+λ,y)+f(x+λ-1,y)+…+f(x,y)-[f(x-1,y)+f(x-2,y)+…+f(x-λ,y)]
Grady(x,y)=f(x,y+λ)+f(x,y+λ-1)+…+f(x,y)-[f(x,y-1)+f(x,y-2)+…+f(x,y-λ)]
Grad=sqrt(Gradx2+Grady2)
其中Gradx代表X方向梯度,Grady代表Y方向梯度,Grad代表总梯度大小,f(x,y)代表图像中(x,y)像素灰度,x为天体图像横坐标、y为天体图像纵坐标,上述λ代表梯度计算长度取值为3~5,sqrt为求均方符号。
本发明与现有技术相比的优点为:
(1)空间目标的晨昏线若不进行判断并剔除则将直接影响目标天体图像中心位置的拟合精度,本发明解决了行星际探测中目标天体晨昏线判断的问题,通过晨昏线判断找到真正的边缘点进行目标天体中心位置拟合,提升了天体图像目标中心的定位精度(由1°误差提高提升到0.05°)。
(2)本发明为高鲁棒性(适用于不同星相)、高刷新率(优于0.5Hz)、高精度(优于0.05°)导航天体敏感器的研制具有重要工程使用价值。
附图说明
图1为深空天体目标图像的示意图;
图2为本发明步骤(3)的弧段选择图;
图3为本发明实现流程图。
具体实施方式
下面结合附图及实施例对本发明进一步详细说明。
如图1所示为输入的目标天体图像,因为太阳照射及敏感器拍摄位置引起的目标天体图像(a,b均存在晨昏线c不存在晨昏线),存在晨昏线的目标天体图像显然无法直接由其360°的边缘点拟合出中心位置。
如图3所示,本发明具体实施步骤如下:
(1)目标天体图像的边缘提取,使用阈值法对目标天体图像进行分割,然后使用切线法提取天体图像的边缘点;
(2)计算边缘点梯度大小,计算目标天体图像360°所有边缘点的梯度大小如下:
Gradx(x,y)=f(x+λ,y)+f(x+λ-1,y)+…+f(x,y)-[f(x-1,y)+f(x-2,y)+…+f(x-λ,y)]
Grady(x,y)=f(x,y+λ)+f(x,y+λ-1)+…+f(x,y)-[f(x,y-1)+f(x,y-2)+…+f(x,y-λ)]
Grad=sqrt(Gradx2+Grady2)
其中Gradx代表X方向梯度,Grady代表Y方向梯度,Grad代表总梯度大小,f(x,y)代表图像中(x,y)像素灰度,x为天体图像横坐标、y为天体图像纵坐标,上述λ代表梯度计算长度取值为3~5,sqrt为求均方符号。
(3)候选弧段选择,如图2所示从0°开始以90°弧段为单位,每间隔30°进行一次弧段选择,逆时针第1个弧段为[0°~90°],第2弧段为[30°~120°]以此类推共12个弧段,其中第12个弧段为[330°~420°],正对关系为1弧段对7弧段,2弧段对8弧段,3弧段对9弧段,4弧段对10弧段,5弧段对11弧段,6弧段对12弧段。
计算该90°弧段内所含有边缘点梯度大小的均值,对所有弧段进行梯度均值排序,选择均值最小的弧段作为晨昏线候选弧段,同时将晨昏候选弧段正对的90°弧段选择为边缘候选弧段;
(4)拟合判断
对上述晨昏候选弧段、边缘候选弧段分别进行圆拟合,若晨昏候选弧段拟合中心与边缘候选弧段拟合中心位置的偏差小于等于设定阈值(一般选择为0.2~3个像元),则说明目标天体成像中没有晨昏线,将目标天体图像360°所有边缘点进行重新拟合得到最终的目标天体图像中心位置;若晨昏候选弧段拟合中心与边缘候选弧段拟合中心位置的偏差大于设定阈值(一般选择为0.2~3个像元),则进行所述圆拟合误差比较,若晨昏线候选弧段圆拟合误差大于边缘候选弧段圆拟合误差,则说明目标天体图像存在晨昏线且位置在候选弧段附近,此时将边缘候选弧段在起始端点、终止端点分别延伸30°将原来90°弧段延长形成150°的边缘点重新进行拟合得到最终的目标天体图像中心位置;若晨昏线候选弧段圆拟合误差小于等于边缘候选弧段圆拟合误差,则转入步骤(5)进行晨昏线补充判断;
(5)晨昏线补充判断
按照步骤(3)的弧段划分方法,以90°弧段为单位,进行所有弧段进行圆拟合并记录拟合误差,拟合输出中心位置点坐标(Xi,Yi),对XY坐标分别进行排序选择中间值(Xm,Ym)设为O点,距离O点最远的输出中心位置点(Xj,Yj)对应的拟合弧段为晨昏候选弧段,选择其正对的弧段为边缘候选弧段,将边缘候选弧段在起始端、终止端分别扩张30°形成150°的边缘点重新进行拟合得到最终的目标天体图像中心位置。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (3)
1.一种用于深空探测自主导航的目标天体图像晨昏线判断方法,其特征在于实现步骤为:
(1)目标天体图像的边缘提取,使用传统边缘算子法进行边缘点提取;
(2)计算天体图像360°所有边缘点的梯度大小;
(3)从0°开始以90°弧段为单位,每间隔30°进行一次弧段选择,并计算该90°弧段内所含有边缘点梯度大小的均值,对所有弧段进行梯度均值排序,选择均值最小的弧段作为晨昏线候选弧段,同时将晨昏候选弧段正对的90°弧段选择为边缘候选弧段;
(4)对上述晨昏候选弧段、边缘候选弧段分别进行圆拟合,若晨昏候选弧段拟合中心与边缘候选弧段拟合中心位置的偏差小于等于设定阈值,则说明目标天体成像中没有晨昏线,将目标天体图像360°所有边缘点进行重新拟合得到最终的目标天体图像中心位置;若晨昏候选弧段拟合中心与边缘候选弧段拟合中心位置的偏差大于设定阈值,则进行所述圆拟合误差比较,若晨昏线候选弧段圆拟合误差大于边缘候选弧段圆拟合误差,则说明目标天体图像存在晨昏线且位置在候选弧段附近,此时将边缘候选弧段在起始端点、终止端点分别延伸30°将原来90°弧段延长形成150°的边缘点重新进行拟合得到最终的目标天体图像中心位置;若晨昏线候选弧段圆拟合误差小于等于边缘候选弧段圆拟合误差,则转入步骤(5)进行晨昏线补充判断;
(5)晨昏线补充判断,按照步骤(3)的弧段划分方法,以90°弧段为单位,进行所有弧段进行圆拟合并记录拟合误差,拟合输出中心位置点坐标(Xi,Yi),对XY坐标分别进行排序选择中间值(Xm,Ym)设为O点,距离O点最远的输出中心位置点(Xj,Yj)对应的拟合弧段为晨昏候选弧段,选择其正对的弧段为边缘候选弧段,将边缘候选弧段在起始端、终止端分别扩张30°形成150°的边缘点重新进行拟合得到最终的目标天体图像中心位置;
所述步骤(2)中计算天体图像360°所有边缘点的梯度大小如下:
Gradx(x,y)=f(x+λ,y)+f(x+λ-1,y)+...+f(x,y)-[f(x-1,y)+f(x-2,y)+...+f(x-λ,y)]
Grady(x,y)=f(x,y+λ)+f(x,y+λ-1)+...+f(x,y)-[f(x,y-1)+f(x,y-2)+...+f(x,y-λ)]
Grad=sqrt(Gradx2+Grady2)
其中Gradx代表X方向梯度,Grady代表Y方向梯度,Grad代表总梯度大小,f(x,y)代表图像中(x,y)像素灰度,x为天体图像横坐标、y为天体图像纵坐标,λ为梯度计算长度,sqrt为求均方符号。
2.根据权利要求1所述的用于深空探测自主导航的目标天体图像晨昏线判断方法,其特征在于:所述步骤(4)中的设定阈值为0.2~3个像元。
3.根据权利要求1所述的用于深空探测自主导航的目标天体图像晨昏线判断方法,其特征在于:λ取值为3~5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210409052.5A CN102944239B (zh) | 2012-10-24 | 2012-10-24 | 一种用于深空探测自主导航的目标天体图像晨昏线判断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210409052.5A CN102944239B (zh) | 2012-10-24 | 2012-10-24 | 一种用于深空探测自主导航的目标天体图像晨昏线判断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102944239A CN102944239A (zh) | 2013-02-27 |
CN102944239B true CN102944239B (zh) | 2015-08-19 |
Family
ID=47727204
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210409052.5A Active CN102944239B (zh) | 2012-10-24 | 2012-10-24 | 一种用于深空探测自主导航的目标天体图像晨昏线判断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102944239B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106353032B (zh) * | 2015-10-10 | 2019-03-29 | 北京控制与电子技术研究所 | 一种欠光照条件下的天体形心快速检测方法 |
CN107389072B (zh) * | 2017-08-15 | 2020-07-07 | 北京理工大学 | 一种小天体表面导航特征区域检测方法 |
CN111402283B (zh) * | 2020-02-25 | 2023-11-10 | 上海航天控制技术研究所 | 基于灰度方差导数的火星图像边缘特征自适应提取方法 |
CN111695564B (zh) * | 2020-06-12 | 2023-11-14 | 上海航天控制技术研究所 | 一种目标识别与导航方法 |
CN111710016B (zh) * | 2020-06-19 | 2023-07-11 | 北京世冠金洋科技发展有限公司 | 一种晨昏线图层获取方法及装置 |
CN114252151A (zh) * | 2021-11-29 | 2022-03-29 | 中国科学院国家空间科学中心 | 一种星载光学仪器在轨工作模式自主调整方法及其装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101236092A (zh) * | 2008-01-31 | 2008-08-06 | 北京控制工程研究所 | 紫外导航敏感器 |
-
2012
- 2012-10-24 CN CN201210409052.5A patent/CN102944239B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101236092A (zh) * | 2008-01-31 | 2008-08-06 | 北京控制工程研究所 | 紫外导航敏感器 |
Non-Patent Citations (2)
Title |
---|
《一种成像敏感器对月定姿算法》;王立等;《宇航学报》;20070131;第28卷(第1期);第39-42页 * |
《深空探测自主导航与控制技术综述》;王大轶等;《空间控制技术与应用》;20090630;第35卷(第3期);第6-12页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102944239A (zh) | 2013-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102944239B (zh) | 一种用于深空探测自主导航的目标天体图像晨昏线判断方法 | |
EP3132231B1 (en) | A method and system for estimating information related to a vehicle pitch and/or roll angle | |
Stockdale et al. | Extracting ecological information from oblique angle terrestrial landscape photographs: Performance evaluation of the WSL Monoplotting Tool | |
CN109612438B (zh) | 一种虚拟共面条件约束下的空间目标初轨确定方法 | |
CN102128623B (zh) | 一种星敏感器快速星图识别方法 | |
CN108225307B (zh) | 一种惯性测量信息辅助的星图匹配方法 | |
CN105913435B (zh) | 一种适用于大区域的多尺度遥感影像匹配方法及系统 | |
CN102706352A (zh) | 一种航空用线状目标的矢量地图匹配导航方法 | |
US10234533B2 (en) | Position determination by observing a celestial object transit the sun or moon | |
CN106537409B (zh) | 确定影像的罗盘定位 | |
CN105160125A (zh) | 一种星敏感器四元数的仿真分析方法 | |
CN104236553B (zh) | 一种全天时的星光折射卫星自主定位方法 | |
CN103837150A (zh) | 一种ccd天顶望远镜地面快速天文定位的方法 | |
Yuan et al. | Automated building segmentation and damage assessment from satellite images for disaster relief | |
Jin et al. | Small scale crater detection based on deep learning with multi-temporal samples of high-resolution images | |
CN103017773B (zh) | 一种基于天体表面特征和天然卫星路标的环绕段导航方法 | |
Hua-Ming et al. | Design and verification of star-map simulation software based on CCD star tracker | |
CN107421504A (zh) | 月基对地观测光电影像的拍摄时间计算方法 | |
CN105631849A (zh) | 多边形目标的变化检测方法及装置 | |
CN112053402B (zh) | 利用偏振地理信息数据库获取航向角的方法 | |
Shaltout et al. | Studying the orientations of luxor ancient Egyptian temples using QuickBird images | |
Furgale et al. | A comparison of global localization algorithms for planetary exploration | |
CN111024063B (zh) | 一种基于星点重提取的大机动条件下星图识别算法 | |
Cristofaro et al. | Approaching skyscape archaeology: A note on method and fieldwork for the case study of Pompeii | |
CN103175526A (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 |