CN104567879B - 一种组合视场导航敏感器地心方向提取方法 - Google Patents
一种组合视场导航敏感器地心方向提取方法 Download PDFInfo
- Publication number
- CN104567879B CN104567879B CN201510041030.1A CN201510041030A CN104567879B CN 104567879 B CN104567879 B CN 104567879B CN 201510041030 A CN201510041030 A CN 201510041030A CN 104567879 B CN104567879 B CN 104567879B
- Authority
- CN
- China
- Prior art keywords
- earth
- fixed star
- image
- core
- asterism
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
- G01C21/025—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
一种组合视场导航敏感器地心方向提取方法,包括以下步骤:对敏感器拍摄的地球图像进行自动阈值计算,通过阈值将图像分割为二值图像。通过区域搜索的方法确定二值图像中地球的边缘。对边缘点进行圆中心和半径的拟合,求得地心像素坐标;对敏感器拍摄的星点图像进行星点提取,根据采样点集的灰度确定星点灰度阈值,并按照此阈值进行星点提取,对满足要求的星点进行保留。然后采用质心算法计算星点质心坐标,并按照典型的星图识别方法进行恒星识别,以此得到星点像素坐标和恒星惯性空间指向信息;利用上两步计算出的地心和恒星像素坐标,计算地心方向与恒星方向在敏感器坐标系下的矢量,再根据通过星图识别得到的恒星惯性空间指向,采用最小二乘算法计算出惯性系中的地心方向矢量。
Description
技术领域
本发明属于图像信息处理技术领域。
背景技术
随着航天技术的发展,对卫星自主导航能力提出了迫切需求。卫星自主天文导航***是基于惯性系下的地心方向矢量的观测量实现的。为了实现地心方向矢量的高精度测量,设计了组合视场的导航敏感器,其视场包括中心视场和环形视场,中心视场对地球进行成像测量,环形视场对恒星进行观测,实现对不同天体目标的同光学***、共探测器成像。
组合视场导航敏感器在同一像面上获得地球圆盘图像和恒星星点图像,传统处理方法是分别对地球圆盘图像和恒星星点图像进行处理。通过地球圆盘图像处理,给出敏感器坐标系下的地心方向矢量;通过恒星星点图像处理,给出敏感器坐标系到惯性系的坐标转换矩阵,再利用坐标转换矩阵将地心方向矢量由敏感器坐标系转换到惯性系。传统处理方法中敏感器坐标系到惯性系的坐标转换矩阵的计算精度比较低。卫星自主天文导航***通常采用惯性系中的地心方向矢量作为测量数据,为便于导航***使用,并避免地心方向矢量坐标变换带来的精度损失,提出了一种组合视场敏感器地心方向直接提取算法。
发明内容
本发明要解决的技术问题是:克服现有技术的不足,利用组合视场导航敏感器对多个天体目标的成像测量,直接计算自主导航***需要的惯性系中的地心方向矢量,避免地心方向矢量的坐标转换,提高地心方向确定精度。
本发明的技术解决方案是:一种组合视场导航敏感器地心方向提取方法,步骤如下:
第一步,对组合视场导航敏感器拍摄的地球图像进行自动阈值计算,通过阈值将图像分割为二值图像;通过区域搜索的方法确定二值图像中地球的边缘;对边缘点进行圆中心的拟合,求得地心像素坐标;
第二步,在组合视场导航敏感器拍摄的星点图像中确定采样点集,根据采样点集的灰度确定星点灰度阈值,并按照此阈值进行星点提取,计算星点像素坐标,并进行恒星识别,最终得到识别出的恒星像素坐标和恒星惯性空间指向信息;
第三步,利用以上两步计算出的地心像素坐标和恒星像素坐标,计算在敏感器坐标系下的地心方向矢量与恒星方向矢量,再根据得到的恒星惯性空间指向信息,采用最小二乘算法计算出惯性系中的地心方向矢量。
所述第一步中计算的过程为:
(1)采用典型自动阈值计算方法得到阈值,通过阈值对图像进行二值分割,大于阈值的像素标称值赋为1,小于阈值的像素标称值赋为0;
(2)采用区域搜索法得到二值图像中的边缘;
假设输入的二值图像为f(x,y),x=1,2,...,M,y=1,2,...,N,图像的大小为M×N,取一个窗口mp×np(3≤mp<M/2,3≤np<N/2),在图像f(x,y)上遍历,若当前窗口内为1的像素个数大于mp×np/2,则认为该窗口中按照横向方向发生0-1或者1-0跳变的点是边缘点;
(3)对步骤(2)确定的边缘点进行圆中心的拟合,得到地球的地心像素坐标。
所述的第三步具体实现为:
(1)根据第一步和第二步的计算结果计算地球中心和恒星在敏感器坐标系下的方向矢量;假设像素坐标为(u,v),采用以下公式计算方向矢量:
其中dx,dy表示像元的横向和纵向尺寸,f表示组合视场导航敏感器的焦距,均为已知量;
即将通过恒星识别得到的n颗恒星的恒星像素坐标分别代入上述公式得到n个恒星方向矢量S1,S2…Sn;将地心像素坐标代入上述公式得到地心方向矢量E;
(2)采用最小二乘法计算地心在惯性系下的指向:
其中,SI1,SI2…SIn为恒星惯性空间指向信息。
本发明的优点是:本发明利用组合视场导航敏感器对多个天体目标的成像测量,直接计算惯性系中的地心方向矢量。一方面,省略了敏感器测量得到的地心方向矢量从敏感器坐标系到惯性系的坐标转换过程,避免了敏感器绕光轴方向较大的测量误差对地心方向矢量计算精度的影响;同时,利用单一敏感器实现了惯性系中地心方向的高精度测量,不需要其它星敏感器参与导航解算,避免了星敏感器与导航敏感器间相对安装误差的影响。
附图说明
图1为含地球的拍摄原图;
图2为二值分割后的效果;
图3为提取的地球的边缘;
图4为仅含恒星的星点拍摄图;
图5为本发明方法流程图。
具体实施方式
下面结合具体实例对本发明做详细说明,本发明如图5所示,具体实现如下:
第二步,从二值图像中,提取轮廓即边缘点,如图3所示。输入的二值图像为f(x,y),x=1,2,…,1024,y=1,2,…,1024,图像的大小为1024×1024。则取一个窗口25×25,在图像f(x,y)上遍历。若当前窗口内为1的像素个数大于312,则该窗口中按照横向方向发生0-1或者1-0跳变的点认为是边缘点(亦称轮廓点)。
第三步,图像中像素标称值为1的点就是轮廓点,序列轮廓点为(xi,yi),i=1,2,...,ne,ne表示边缘点的个数,ne≥5,根据以下公式计算地球圆形的中心(xc,yc):
利用最小二乘法求得地球圆形的中心坐标(xc,yc)。
第四步,星图阈值确定。在组合视场导航敏感器拍摄的星点图像中取出中心的512×512大小的区域(该区域为采样点集,可以根据实际情况灵活选取),计算区域内像素灰度的平均值Vmean和方差Var。
其中I(i,j)表示区域内(i,j)点的像素灰度,M×N为星点有效图像区域,将阈值设定为:Vmean+Var/2。
第五步,星点中心提取。星点图像如图4所示,对图像中灰度大于阈值的星点区域进行灰度质心提取:
其中Rk为星点k所占图像区域,Iij为有效点的灰度值,i,j分别为该点的像素坐标,为此可得到m个星点像素坐标(xk,yk),k=1,2,...m。
第六步,采用主星识别法进行恒星识别,从最亮星开始计算敏感器视场内对应角距组,然后按照搜索每一个角距的可能星组,根据交集决定主星在导航星库中位置。先以最亮星为主星1,已知星点1和星点2的像素坐标(u1,v1)和(u2,v2),则两星点间角距d可按下式计算:
其中dx,dy表示像元的横向和纵向尺寸,f表示组合视场导航敏感器的焦距,均为已知量。
取导航星库中的三颗星i,j,k组成导航三角形,其两两之间的角距表示为θ(i,j),θ(j,k)和θ(i,k)。敏感器视场内的主星1和2,3两颗星之间的角距表示为d(1,2),d(2,3)和d(1,3),角距之间进行如下判断:
|θ(i,j)-d(1,2)|≤ε
|θ(j,k)-d(2,3)|≤ε
|θ(i,k)-d(1,3)|≤ε
ε表示测量的不确定度(根据实际需要灵活确定)。满足上式且θ(i,j),θ(j,k)和θ(i,k)之间首尾相接,则认为敏感器视场内的主星1和2,3两颗星对应导航星库中的i,j,k星。如果结果不唯一,则保留主星1,重新选取敏感器视场内的4,5两颗星组成新的三角形,在多个结果中再进行筛选,直至结果唯一。之后再对敏感器视场内的其它星执行以上识别操作,准确识别最少三颗星。导航星库中存储了各个恒星在惯性系下的方向矢量信息,因此,通过星图识别可以得到敏感器视场中的恒星惯性空间指向信息。
第七步,计算地心矢量。
假设像素坐标为(u,v),计算敏感器坐标系下的地心方向矢量与恒星方向矢量的公式为:
其中dx,dy表示像元的横向和纵向尺寸,f表示组合视场导航敏感器的焦距,均为已知量。
即将通过恒星识别得到的n颗恒星的恒星像素坐标分别代入上述公式得到n个恒星方向矢量S1,S2…Sn;将地心像素坐标代入上述公式得到地心方向矢量E;
利用敏感器坐标系下地心方向矢量E、恒星方向矢量S1,S2…Sn(n表示通过星图识别得到的恒星数目)以及通过星图识别得到的敏感器视场中的恒星惯性空间指向信息SI1,SI2…SIn,采用最小二乘法直接拟合计算惯性系下地心方向矢量EI。
经过实验结果的统计表明,该方法计算得到的惯性系下地心矢量精度优于0.02度。
本发明未详细说明部分属于本领域技术人员公知常识。
Claims (2)
1.一种组合视场导航敏感器地心方向提取方法,其特征在于步骤如下:
第一步,对组合视场导航敏感器拍摄的地球图像进行自动阈值计算,通过阈值将图像分割为二值图像;通过区域搜索的方法确定二值图像中地球的边缘;对边缘点进行圆中心的拟合,求得地心像素坐标;
第二步,在组合视场导航敏感器拍摄的星点图像中确定采样点集,根据采样点集的灰度确定星点灰度阈值,并按照此阈值进行星点提取,计算星点像素坐标,并进行恒星识别,最终得到识别出的恒星像素坐标和恒星惯性空间指向信息;
第三步,利用以上两步计算出的地心像素坐标和恒星像素坐标,计算在敏感器坐标系下的地心方向矢量与恒星方向矢量,再根据得到的恒星惯性空间指向信息,采用最小二乘算法计算出惯性系中的地心方向矢量;
所述的第三步具体实现为:
(1)根据第一步和第二步的计算结果计算地球中心和恒星在敏感器坐标系下的方向矢量;假设像素坐标为(u,v),采用以下公式计算方向矢量:
其中dx,dy表示像元的横向和纵向尺寸,f表示组合视场导航敏感器的焦距,均为已知量;
即将通过恒星识别得到的n颗恒星的恒星像素坐标分别代入上述公式得到n个恒星方向矢量S1,S2…Sn;将地心像素坐标代入上述公式得到地心方向矢量E;
(2)采用最小二乘法计算地心在惯性系下的指向:
其中,SI1,SI2…SIn为恒星惯性空间指向信息;EI为惯性系下地心方向矢量。
2.根据权利要求1所述的一种组合视场导航敏感器地心方向提取方法,其特征在于:所述第一步中计算的过程为:
(1)采用典型自动阈值计算方法得到阈值,通过阈值对图像进行二值分割,大于阈值的像素标称值赋为1,小于阈值的像素标称值赋为0;所述的典型自动阈值计算方法为按照公式计算最大值,将J(t)最大时对应的t值记为阈值Dth;其中, P(i)为每个灰度的像素个数占整个图像总像素的比值,t=0~255;
(2)采用区域搜索法得到二值图像中的边缘;
假设输入的二值图像为f(x,y),x=1,2,...,M,y=1,2,...,N,图像的大小为M×N,取一个窗口mp×np,3≤mp<M/2,3≤np<N/2,在图像f(x,y)上遍历,若当前窗口内为1的像素个数大于mp×np/2,则认为该窗口中按照横向方向发生0-1或者1-0跳变的点是边缘点;
(3)对步骤(2)确定的边缘点进行圆中心的拟合,得到地球的地心像素坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510041030.1A CN104567879B (zh) | 2015-01-27 | 2015-01-27 | 一种组合视场导航敏感器地心方向提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510041030.1A CN104567879B (zh) | 2015-01-27 | 2015-01-27 | 一种组合视场导航敏感器地心方向提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104567879A CN104567879A (zh) | 2015-04-29 |
CN104567879B true CN104567879B (zh) | 2018-08-21 |
Family
ID=53084493
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510041030.1A Active CN104567879B (zh) | 2015-01-27 | 2015-01-27 | 一种组合视场导航敏感器地心方向提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104567879B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107367275A (zh) * | 2017-06-19 | 2017-11-21 | 北京控制工程研究所 | 一种高轨卫星自主导航敏感器 |
CN109949204B (zh) * | 2019-03-29 | 2023-08-15 | 江苏亿通高科技股份有限公司 | 流水线结构的星点质心提取电路 |
CN111402283B (zh) * | 2020-02-25 | 2023-11-10 | 上海航天控制技术研究所 | 基于灰度方差导数的火星图像边缘特征自适应提取方法 |
CN111323020B (zh) * | 2020-02-25 | 2021-11-16 | 上海航天控制技术研究所 | 一种基于火星边缘及中心多矢量观测的自主定轨方法 |
CN113532446B (zh) * | 2021-07-20 | 2022-07-05 | 北京控制工程研究所 | 一种基于迭代遍历的星敏感器抗杂光星点提取方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5319969A (en) * | 1992-09-21 | 1994-06-14 | Honeywell Inc. | Method for determining 3-axis spacecraft attitude |
US5837894A (en) * | 1995-02-09 | 1998-11-17 | Honeywell Inc. | Wide field of view sensor with diffractive optic corrector |
CN101236092A (zh) * | 2008-01-31 | 2008-08-06 | 北京控制工程研究所 | 紫外导航敏感器 |
CN101713655A (zh) * | 2009-11-16 | 2010-05-26 | 北京控制工程研究所 | 一种紫外导航敏感器的双模式控制方法 |
CN102519455A (zh) * | 2011-12-08 | 2012-06-27 | 北京控制工程研究所 | 基于紫外敏感器的自主导航半物理仿真试验*** |
CN102927973A (zh) * | 2012-10-24 | 2013-02-13 | 北京控制工程研究所 | 一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法 |
CN103148849A (zh) * | 2013-03-12 | 2013-06-12 | 北京控制工程研究所 | 基于地月卫星联合测距和紫外敏感器的组合导航方法 |
CN103389099A (zh) * | 2013-07-26 | 2013-11-13 | 西安电子科技大学 | 基于x射线脉冲星的航天器姿态和位置测量***及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10123A (en) * | 1853-10-11 | Improvement in cultivate ng-plows |
-
2015
- 2015-01-27 CN CN201510041030.1A patent/CN104567879B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5319969A (en) * | 1992-09-21 | 1994-06-14 | Honeywell Inc. | Method for determining 3-axis spacecraft attitude |
US5837894A (en) * | 1995-02-09 | 1998-11-17 | Honeywell Inc. | Wide field of view sensor with diffractive optic corrector |
CN101236092A (zh) * | 2008-01-31 | 2008-08-06 | 北京控制工程研究所 | 紫外导航敏感器 |
CN101713655A (zh) * | 2009-11-16 | 2010-05-26 | 北京控制工程研究所 | 一种紫外导航敏感器的双模式控制方法 |
CN102519455A (zh) * | 2011-12-08 | 2012-06-27 | 北京控制工程研究所 | 基于紫外敏感器的自主导航半物理仿真试验*** |
CN102927973A (zh) * | 2012-10-24 | 2013-02-13 | 北京控制工程研究所 | 一种用于深空探测自主导航的目标天体亚像素图像快速边缘定位方法 |
CN103148849A (zh) * | 2013-03-12 | 2013-06-12 | 北京控制工程研究所 | 基于地月卫星联合测距和紫外敏感器的组合导航方法 |
CN103389099A (zh) * | 2013-07-26 | 2013-11-13 | 西安电子科技大学 | 基于x射线脉冲星的航天器姿态和位置测量***及方法 |
Non-Patent Citations (7)
Title |
---|
一种卫星天文自主定轨定姿方法研究;李琳琳等;《空间科学学报》;20030331;第23卷(第2期);全文 * |
一种基于紫外姿态敏感器和星间相对测量的星座组合导航方法;王裙等;《航天控制》;20131031;第31卷(第5期);全文 * |
一种姿态机动辅助下的天文导航***偏差自校准方法;魏春岭等;《宇航学报》;20100131;第31卷(第1期);全文 * |
利用地球紫外和恒星可见光的卫星自主导航;孙俊等;《光学精密工程》;20130531;第21卷(第5期);全文 * |
基于紫外敏感器的地月转移轨道慢旋探测器自主导航算法;乔国栋等;《宇航学报》;20090331;第30卷(第2期);全文 * |
基于紫外敏感器的航天器自主导航;魏春岭等;《航天控制》;20040630;第22卷(第3期);正文第2节 * |
组合大视场星敏感器星光折射卫星自主导航方法及仿真;王国权等;《***仿真学报》;20050331;第17卷(第3期);正文第2-6节及图1-5 * |
Also Published As
Publication number | Publication date |
---|---|
CN104567879A (zh) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105761242B (zh) | 一种基于计算机双目视觉与惯性测量的盲人行走定位方法 | |
CN104567879B (zh) | 一种组合视场导航敏感器地心方向提取方法 | |
CN107451593B (zh) | 一种基于图像特征点的高精度gps定位方法 | |
US9453731B2 (en) | System and method for determining orientation relative to earth | |
US11468654B2 (en) | Determining compass orientation of imagery | |
CN109443359B (zh) | 一种地面全景影像的地理定位方法 | |
CN103256920A (zh) | 利用图像处理测定倾斜角和倾斜方向 | |
CN102034238A (zh) | 基于光学成像测头和视觉图结构的多摄像机***标定方法 | |
CN109540113B (zh) | 一种全站仪及其星图识别方法 | |
CN104077760A (zh) | 一种航空摄影测量的快速拼接***及其实现方法 | |
EP3132231A1 (en) | A method and system for estimating information related to a vehicle pitch and/or roll angle | |
JP2016224953A (ja) | ナビゲーション調整のためのクロススペクトル特徴相関 | |
CN108549376A (zh) | 一种基于信标的导航定位方法及*** | |
CN105806346B (zh) | 中高轨星座星间照相观测敏感器及其星间角距测量方法 | |
US20180074154A1 (en) | Position Determination by Observing a Celestial Object Transit the Sun or Moon | |
CN111238540A (zh) | 一种基于恒星拍摄的珞珈一号相机-星敏安装标定方法 | |
US8965053B2 (en) | Method for remotely determining an absolute azimuth of a target point | |
CN104133993A (zh) | 基于误差矢量匹配的小视场恒星匹配方法 | |
CN115790539B (zh) | 一种合作靶标水下摄影测量方法 | |
CN107992677B (zh) | 基于惯导信息和亮度修正的红外弱小移动目标跟踪方法 | |
JP2010009236A (ja) | 平面領域推定装置及びプログラム | |
CN109785388A (zh) | 一种基于双目摄像头的短距离精确相对定位方法 | |
CN104599283A (zh) | 一种基于深度差恢复相机高度的图像深度改进方法 | |
CN109827578B (zh) | 基于轮廓相似性的卫星相对姿态估计方法 | |
CN107607737B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |