CN102927973A - 一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法 - Google Patents
一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法 Download PDFInfo
- Publication number
- CN102927973A CN102927973A CN201210409261XA CN201210409261A CN102927973A CN 102927973 A CN102927973 A CN 102927973A CN 201210409261X A CN201210409261X A CN 201210409261XA CN 201210409261 A CN201210409261 A CN 201210409261A CN 102927973 A CN102927973 A CN 102927973A
- Authority
- CN
- China
- Prior art keywords
- celestial body
- rho
- grad
- target celestial
- target
- 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
- Image Analysis (AREA)
Abstract
一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法,步骤为:目标天体图像的边缘提取;对提取得到的目标成像区间计算中心位置,得到目标天体边缘点极坐标;根据极坐标进行梯度计算;根据梯度计算结果得到最优亚像素位置计算输出。本发明解决了深空探测目标天体中心位置快速、高精度测量问题。
Description
技术领域
本发明涉及一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法,属于空间光学成像敏感器。
背景技术
在深空探测中需要对火星等其他天体目标进行拍照并确定目标天体中心位置。继承紫外月球敏感器的图像处理算法,利用目标天体图像边缘点的拟合来实现精确的中心定位,因此边缘点的定位精确性将直接影响到目标天体中心位置的确定精度。因此,深空探测敏感器提升精度的方式之一就是提高边缘点的定位精度。文献1,刘力双,张铫,卢慧卿,图像的快速亚像素边缘检测方法,光电子激光,2005(8),使用Sobel算子进行边缘粗定位然后使用最小二乘拟合算法确定亚像素位置,但未能考虑到梯度方向导致精度不足;文献2艾泽潭,石庚辰,代俊,微零件图像亚像素边缘定位算法,北京理工大学学报,2011(3),二值化后使用插值方法确定亚像素边缘位置,该方法适用于二值图而不是空间灰度图;文献3,屈玉福,浦昭邦,王亚爱,视觉检测系统中亚像素边缘检测技术的对比研究,仪器仪表学报2003(z1),综述了插值法、空间矩、最小二乘的亚像素定位方法,上述方法计算量偏大。
已有文献均无法解决空间天体目标高速高精度的亚像素边缘定位问题。
发明内容
本发明的技术解决问题是:克服现有技术的不足之处,提供一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法,解决了深空探测目标天体中心矢量快速、高精度测量问题。
本发明技术解决方案:一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法,实现步骤为:
一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法,其特征在于实现步骤为:
(1)目标天体图像的边缘提取,使用阈值法进行目标天体图像分割得到目标成像区间(x1~x2,y1~y2),x代表横坐标、y代表纵坐标,从x1开始进行横向切线,第一条切线为x1,第二条切线为x1+d一直到x2;纵向切线从y1开始,按照d间隔进行纵向切线;横向、纵向切线与目标的交点为边缘点,每一个切线对应2个边缘点,边缘点记为f(x,y),f代表输入的图像阵列灰度;
(2)对步骤(1)得到的目标成像区间计算中心位置:
x0=(x1+x2)/2 y0=(y1+y2)/2
对步骤(1)输出的边缘点坐标f(x,y),按照如下公式计算极坐标:
得到目标天体边缘点极坐标表示f(ρ,θ),ρ代表极径、θ代表角度;
(3)计算f(ρ,θ)两侧的点f(ρ-1,θ)、f(ρ+1,θ)共3点的梯度值,梯度计算如下:
Grad(ρ,θ)=f(ρ+λ,θ)+f(ρ+λ-1,θ)+…+f(ρ,θ)-[f(ρ-1,θ)+f(ρ-2,θ)+…+f(ρ-λ,θ)]
其中Grad代表梯度,λ代表梯度计算长度,其中f(ρ,θ)与f(x,y)对应的位置坐标如下计算:
x=ρ·cos(θ) y=ρ·sin(θ)
得到x,y坐标值,若出现小数则按照四舍五入原则进行取整;
(4)最优亚像素位置计算输出,按照以下公式输出
xnew=ρ_new·cos(θ)ynew=ρ_new·sin(θ)。
所述步骤(1)中间隔d为1~5像元。
所述步骤(3)中的λ的取值为2~5。
本发明与现有技术相比的优点为:本发明解决了行星际探测中目标天体边缘的快速高精度定位问题,高刷新率(优于0.5Hz)、高精度(优于0.05°)导航天体敏感器的研制具有重要工程使用价值。
附图说明
图1为切线法获得边缘示意图;
图2为本发明实现流程图。
具体实施方式
下面结合附图及实施例对本发明进一步详细说明。
如图1所示,图像阈值分割后的边缘点提取示意(横向、纵向每一条切线与目标交点为边缘点,每一条切线对应2个边缘点,d为切线的间隔)
如图2所示,本发明具体实施步骤如下:
(1)目标天体图像的边缘提取,使用阈值法进行目标天体图像分割得到目标成像区间(x1~x2,y1~y2),从x1开始进行横向切线,第一条切线为x1,第二条切线为x1+d一直到x2;纵向切线从y1开始,按照d间隔进行纵向切线,间隔d一般选择1~5像元;横向、纵向切线与目标的交点为边缘点,每一个切线对应2个边缘点,边缘点坐标记为f(x,y);
(2)对步骤(1)得到的成像区间计算中心位置:
x0=(x1+x2)/2 y0=(y1+y2)/2
对步骤(1)输出的边缘点坐标f(x,y),按照如下公式计算极坐标:
得到边缘点极坐标表示f(ρ,θ);
(3)计算f(ρ,θ)两侧的点f(ρ-1,θ)、f(ρ+1,θ)共3点的梯度值,梯度计算办法如下:
Grad(ρ,θ)=f(ρ+λ,θ)+f(ρ+λ-1,θ)+…+f(ρ,θ)-[f(ρ-1,θ)+f(ρ-2,θ)+…+f(ρ-λ,θ)]
其中λ的取值一般为2~5。其中f(ρ,θ)与f(x,y)对应的位置坐标如下计算:
x=ρ·cos(θ)y=ρ·sin(θ)
得到x,y坐标值,若出现小数则按照四舍五入原则进行取整。
(4)最优亚像素位置计算输出,按照以下公式输出
亚像素边缘点位置坐标如下:
xnew=ρ_new·cos(θ) y new=ρ_new·sin(θ)
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (3)
1.一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法,其特征在于实现步骤为:
(1)目标天体图像的边缘提取,使用阈值法进行目标天体图像分割得到目标成像区间(x1~x2,y1~y2),x代表横坐标、y代表纵坐标,从x1开始进行横向切线,第一条切线为x1,第二条切线为x1+d一直到x2;纵向切线从y1开始,按照d间隔进行纵向切线;横向、纵向切线与目标的交点为边缘点,每一个切线对应2个边缘点,边缘点记为f(x,y),f代表输入的图像阵列灰度;
(2)对步骤(1)得到的目标成像区间计算中心位置:
x0=(x1+x2)/2y0=(y1+y2)/2
对步骤(1)输出的边缘点坐标f(x,y),按照如下公式计算极坐标:
得到目标天体边缘点极坐标表示f(ρ,θ),ρ代表极径、θ代表角度;
(3)计算f(ρ,θ)两侧的点f(ρ-1,θ)、f(ρ+1,θ)共3点的梯度值,梯度计算如下:
Grad(ρ,θ)=f(ρ+λ,θ)+f(ρ+λ-1,θ)+…+f(ρ,θ)-[f(ρ-1,θ)+f(ρ-2,θ)+…+f(ρ-λ,θ)]
其中Grad代表梯度,λ代表梯度计算长度,其中f(ρ,θ)与f(x,y)对应的位置坐标如下计算:
x=ρ·cos(θ) y=ρ·sin(θ)
得到x,y坐标值,若出现小数则按照四舍五入原则进行取整;
(4)最优亚像素位置计算输出,按照以下公式输出
得到亚像素边缘点位置坐标如下:
xnew=ρ_new·cos(θ)ynew=ρ_new·sin(θ)。
2.根据权利要求1所述的一种用于深空探测自主导航的目标天体图像切线法边缘点确定方法,其特征在于:所述步骤(1)中间隔d为1~5像元。
3.根据权利要求1所述的一种用于深空探测自主导航的目标天体图像切线法边缘点确定方法,其特征在于:所述步骤(3)中的λ的取值为2~5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210409261.XA CN102927973B (zh) | 2012-10-24 | 2012-10-24 | 一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210409261.XA CN102927973B (zh) | 2012-10-24 | 2012-10-24 | 一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102927973A true CN102927973A (zh) | 2013-02-13 |
CN102927973B CN102927973B (zh) | 2015-07-08 |
Family
ID=47642831
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210409261.XA Active CN102927973B (zh) | 2012-10-24 | 2012-10-24 | 一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102927973B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104567879A (zh) * | 2015-01-27 | 2015-04-29 | 北京控制工程研究所 | 一种组合视场导航敏感器地心方向提取方法 |
CN109344785A (zh) * | 2018-10-12 | 2019-02-15 | 北京航空航天大学 | 一种在深空自主光学导航中的高精度的行星中心定位方法 |
CN109631912A (zh) * | 2019-01-10 | 2019-04-16 | 中国科学院光电技术研究所 | 一种深空球形目标被动测距方法 |
CN111127501A (zh) * | 2019-12-03 | 2020-05-08 | 重庆邮电大学 | 基于多粒度遗传算法的图像分割方法 |
CN111739039A (zh) * | 2020-08-05 | 2020-10-02 | 北京控制与电子技术研究所 | 基于边缘提取的快速形心定位方法、系统及装置 |
CN115128791A (zh) * | 2021-03-26 | 2022-09-30 | 清华大学 | 光谱成像天文望远镜及天文望远镜光谱成像方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4688092A (en) * | 1986-05-06 | 1987-08-18 | Ford Aerospace & Communications Corporation | Satellite camera image navigation |
CN101676687A (zh) * | 2008-09-17 | 2010-03-24 | 北京航空航天大学 | 一种超高精度的星敏感器 |
-
2012
- 2012-10-24 CN CN201210409261.XA patent/CN102927973B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4688092A (en) * | 1986-05-06 | 1987-08-18 | Ford Aerospace & Communications Corporation | Satellite camera image navigation |
CN101676687A (zh) * | 2008-09-17 | 2010-03-24 | 北京航空航天大学 | 一种超高精度的星敏感器 |
Non-Patent Citations (1)
Title |
---|
张美静等: "改进形态学梯度的样条插值亚像素边缘检测方法", 《沈阳理工大学学报》, vol. 31, no. 3, 30 June 2012 (2012-06-30), pages 19 - 23 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104567879A (zh) * | 2015-01-27 | 2015-04-29 | 北京控制工程研究所 | 一种组合视场导航敏感器地心方向提取方法 |
CN104567879B (zh) * | 2015-01-27 | 2018-08-21 | 北京控制工程研究所 | 一种组合视场导航敏感器地心方向提取方法 |
CN109344785A (zh) * | 2018-10-12 | 2019-02-15 | 北京航空航天大学 | 一种在深空自主光学导航中的高精度的行星中心定位方法 |
CN109344785B (zh) * | 2018-10-12 | 2021-10-01 | 北京航空航天大学 | 一种在深空自主光学导航中的高精度的行星中心定位方法 |
CN109631912A (zh) * | 2019-01-10 | 2019-04-16 | 中国科学院光电技术研究所 | 一种深空球形目标被动测距方法 |
CN109631912B (zh) * | 2019-01-10 | 2022-08-23 | 中国科学院光电技术研究所 | 一种深空球形目标被动测距方法 |
CN111127501A (zh) * | 2019-12-03 | 2020-05-08 | 重庆邮电大学 | 基于多粒度遗传算法的图像分割方法 |
CN111127501B (zh) * | 2019-12-03 | 2023-05-30 | 重庆邮电大学 | 基于多粒度遗传算法的图像分割方法 |
CN111739039A (zh) * | 2020-08-05 | 2020-10-02 | 北京控制与电子技术研究所 | 基于边缘提取的快速形心定位方法、系统及装置 |
CN111739039B (zh) * | 2020-08-05 | 2020-11-13 | 北京控制与电子技术研究所 | 基于边缘提取的快速形心定位方法、系统及装置 |
CN115128791A (zh) * | 2021-03-26 | 2022-09-30 | 清华大学 | 光谱成像天文望远镜及天文望远镜光谱成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102927973B (zh) | 2015-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102927973B (zh) | 一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法 | |
CN106056605B (zh) | 一种基于图像匹配的在轨高精度图像定位方法 | |
CN106643792B (zh) | 惯性测量单元和地磁传感器整体标定装置及标定方法 | |
CN102173313B (zh) | 一种软着陆接力避障方法 | |
CN103926625B (zh) | 一种利用地磁总场对水下磁目标高精度远距离的定位方法 | |
CN102607526A (zh) | 双介质下基于双目视觉的目标姿态测量方法 | |
CN105091880B (zh) | 一种基于标量传感器阵列的追踪定位水下远距离磁性目标的方法 | |
CN103791902B (zh) | 适用于高机动载体的星敏感器自主导航方法 | |
CN103697855A (zh) | 一种基于海天线检测的船体水平姿态测量方法 | |
CN106525001A (zh) | 地球静止轨道遥感卫星相机视轴空间指向计算方法 | |
CN104634248A (zh) | 一种双目视觉下的转轴标定方法 | |
CN103136758A (zh) | 基于正交多项式拟合的快速边缘检测方法 | |
US20190287257A1 (en) | Method and system for measuring the distance to remote objects | |
CN102944227B (zh) | 一种基于fpga实现实时提取恒星星像坐标的方法 | |
Xiaoqi et al. | Centroid automatic extraction of spaceborne laser spot image | |
CN104864851A (zh) | 一种基于矩形周长和面积加权约束的单目视觉位姿测量方法 | |
CN104567879B (zh) | 一种组合视场导航敏感器地心方向提取方法 | |
CN101793517B (zh) | 一种提高机载平台定姿精度的在线快速方法 | |
CN102706348B (zh) | 一种基于三角形的重力图快速匹配方法 | |
Zou et al. | Optimal sensor configuration for positioning seafloor geodetic node | |
CN103630109B (zh) | 一种基于星光折射确定地心矢量的方法 | |
CN105222764B (zh) | 一种对惯性角速度传感器进行地球自转补偿的方法 | |
CN104391272A (zh) | 利用测向数据进行目标定位的方法与系统 | |
Yılmaztürk et al. | Measurement of deflections in buried flexible pipes by close range digital photogrammetry | |
CN102519471A (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 |