CN109443381B - 一种星敏感器质心精度自适应补偿方法 - Google Patents

一种星敏感器质心精度自适应补偿方法 Download PDF

Info

Publication number
CN109443381B
CN109443381B CN201811210853.2A CN201811210853A CN109443381B CN 109443381 B CN109443381 B CN 109443381B CN 201811210853 A CN201811210853 A CN 201811210853A CN 109443381 B CN109443381 B CN 109443381B
Authority
CN
China
Prior art keywords
centroid
window area
calculating
star
self
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
Application number
CN201811210853.2A
Other languages
English (en)
Other versions
CN109443381A (zh
Inventor
杨飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Remote Sensing Equipment
Original Assignee
Beijing Institute of Remote Sensing Equipment
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Remote Sensing Equipment filed Critical Beijing Institute of Remote Sensing Equipment
Priority to CN201811210853.2A priority Critical patent/CN109443381B/zh
Publication of CN109443381A publication Critical patent/CN109443381A/zh
Application granted granted Critical
Publication of CN109443381B publication Critical patent/CN109443381B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/02Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Manufacturing & Machinery (AREA)
  • Astronomy & Astrophysics (AREA)
  • Automation & Control Theory (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种星敏感器质心精度自适应补偿方法,包括:步骤一,确定星点窗口区域;步骤二,计算窗口区域质心;步骤三,计算自适应补偿参数;步骤四,修正窗口区域质心。以星点像斑灰度值最高的像元为中心获取窗口区域图像,计算窗口区域的灰度质心,通过窗口中心邻域的灰度比例关系计算自适应补偿参数,利用补偿参数修正窗口区域的灰度质心。本发明方法通过计算自适应补偿参数并修正质心,来提高星点像斑质心计算精度的问题,获取更高的质心定位精度。

Description

一种星敏感器质心精度自适应补偿方法
技术领域
本发明涉及一种质心精度自适应补偿方法,特别是一种星敏感器质心精度自适应补偿方法。
背景技术
星敏感器是天文导航技术的一种,通过对恒星星空的拍摄和测量,结合恒星在天球坐标系的位置信息,解算出星敏感器光轴在天球坐标系下的矢量方向,进而获取星敏感器所处飞行平台姿态信息。星敏感器是高精度姿态测量设备,其测量精度可达1角秒量级。星敏感器在拍摄恒星成像时,每颗恒星都会在星敏感器成像面形成星点像斑,像斑尺寸覆盖若干像元,为了精确定位星点像斑位置,获取亚像元精度,需要计算星点像斑的质心位置。星点像斑的质心定位精度直接决定了星敏感器的最终测量精度,是一项非常关键的指标。
传统的亚像元星点质心定位方法主要有:一是普通质心法,以像素灰度值作为唯一权重因子进行质心定位;二是平方加权质心法,采用灰度值的平方代替灰度值作为权重因子,突出了离中心较近的灰度值像素点对中心位置的影响;三是带阈值的质心法,将原始图像与背景阈值相减,对原始图像中大于设定阈值的像素点进行质心定位;四是曲面拟合法,用高斯曲面对星点像斑的灰度分布进行拟合,再对拟合曲线采用最小二乘法求出质心的位置坐标。这些方法中以带阈值的质心法精度最高,但是其质心计算精度和星点像斑理论质心在成像面上所处位置相关,当理论质心处于某像素正中心时,质心计算精度最高,随着理论质心移动到像素结合处,质心计算精度逐渐变低,需要某种方法进行补偿。
发明内容
本发明目的在于提供一种星敏感器质心精度自适应补偿方法,解决提高传统质心定位精度的问题。
本发明提供的是一种星敏感器质心精度自适应补偿方法,包括:
步骤一,确定星点窗口区域;
步骤二,计算窗口区域质心;
步骤三,计算自适应补偿参数;
步骤四,修正窗口区域质心。
优选地,所述步骤一中,设星点像斑中灰度最大的像元位置为(x0,y0),以(x0,y0)为中心,以n为半径,选取星点的窗口区域记为W,则窗口区域W内的任意像素位置(x,y)满足:
Figure BDA0001832411740000021
优选地,所述步骤2中计算窗口区域质心,
窗口区域内位置(x,y)的像素灰度为g(x,y),计算x方向窗口区域灰度累加和:
Figure BDA0001832411740000022
计算y方向窗口区域灰度累加和:
Figure BDA0001832411740000023
优选地,所述步骤2中计算窗口区域质心,
计算窗口区域x方向灰度质心:
Figure BDA0001832411740000031
计算窗口区域y方向灰度质心:
Figure BDA0001832411740000032
优选地,所述步骤3中计算自适应补偿参数,
计算x方向窗口质心补偿参数:
Figure BDA0001832411740000033
其中k为修正参数,跟星敏感器的星点像斑成像质量相关,需要通过标定测出。
优选地,所述步骤3中计算自适应补偿参数,
计算y方向窗口质心补偿参数:
Figure BDA0001832411740000034
其中k为修正参数,跟星敏感器的星点像斑成像质量相关,需要通过标定测出。
优选地,所述步骤4中修正窗口区域质心:
x′zx=xzx+Δxzx (8)
优选地,所述步骤4中修正窗口区域y方向质心:
y′zx=yzx+Δyzx (9)
优选地,所述第四步自适应补偿参数的计算方法。
优选地,本发明还涉及一种星点像斑质心精度自适应补偿方法,包括:星点像斑灰度最高像元窗口区域确定模块、窗口区域质心计算模块和自适应补偿参数计算模块;
所述星点像斑灰度最高像元窗口区域确定模块、所述窗口区域质心计算模块和所述自适应补偿参数计算模块顺次连通;所述星点像斑灰度最高像元窗口区域确定模块确定灰度值最高的像元及其窗口区域,窗口区域质心计算模块分别计算窗口区域x方向、y方向质心,所述自适应补偿参数计算模块依据窗口中心邻域的灰度比例关系计算自适应补偿参数并修正各自方向质心,从而完成星点像斑质心定位及补偿。
综上所述,采用了上述技术方案,本发明方法通过计算自适应补偿参数并修正质心,来提高计算精度的问题,获取更高的质心定位精度。
附图说明
图1为本发明的一个实施例的一种星敏感器质心精度自适应补偿方法的流程图;
图2为本发明的一个实施例的x方向移动质心测量误差;
图3为本发明的一个实施例的y方向移动质心测量误差;
图4为本发明的一个实施例的x、y方向同时移动质心测量误差。
具体实施方式
下面结合附图说明根据本发明的具体实施方式。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述的其他方式来实施,因此,本发明并不限于下面公开的具体实施例的限制。
图1为本发明的一个实施例的一种星敏感器质心精度自适应补偿方法的流程图;图2为本发明的一个实施例的x方向移动质心测量误差;图3为本发明的一个实施例的y方向移动质心测量误差;图4为本发明的一个实施例的x、y方向同时移动质心测量误差。
本发明为一种星敏感器质心精度自适应补偿方法的具体步骤为:
第一步,确定星点窗口区域。
设星点像斑中灰度最大的像元位置为(x0,y0),以(x0,y0)为中心,以n为半径,选取星点的窗口区域记为W,则窗口区域W内的任意像素位置(x,y)满足:
Figure BDA0001832411740000051
第二步,计算窗口区域质心。
窗口区域内位置(x,y)的像素灰度为g(x,y),计算x方向窗口区域灰度累加和:
Figure BDA0001832411740000052
计算y方向窗口区域灰度累加和:
Figure BDA0001832411740000053
计算窗口区域x方向灰度质心:
Figure BDA0001832411740000054
计算窗口区域y方向灰度质心:
Figure BDA0001832411740000055
第三步,计算自适应补偿参数。
计算x方向窗口质心补偿参数:
Figure BDA0001832411740000061
计算y方向窗口质心补偿参数:
Figure BDA0001832411740000062
其中k为修正参数,跟星敏感器的星点像斑成像质量相关,需要通过标定测出。
第四步,修正窗口区域质心。
修正窗口区域x方向质心:
x′zx=xzx+Δxzx (8)
修正窗口区域y方向质心:
y′zx=yzx+Δyzx (9)
一种星敏感器质心精度自适应补偿方法,包括:星点像斑灰度最高像元窗口区域确定模块、窗口区域质心计算模块和自适应补偿参数计算模块;星点像斑灰度最高像元窗口区域确定模块、窗口区域质心计算模块和自适应补偿参数计算模块顺次连通;其中星点像斑灰度最高像元窗口区域确定模块确定灰度值最高的像元及其窗口区域,窗口区域质心计算模块分别计算窗口区域x方向、y方向质心,自适应补偿参数计算模块依据窗口中心邻域的灰度比例关系计算自适应补偿参数并修正各自方向质心,从而完成星点像斑质心定位及补偿,提高质心定位精度。
利用本方法对某星敏感器进行星点像斑质心定位测试,测试结果为:
仅移动x方向,y方向质心保持不变,移动范围为从一个像元中心移动到下一个像元中心,星点像斑测量误差如图2所示,测量结果如表1所示:
表1 x方向移动质心测量结果
Figure BDA0001832411740000071
依据表1可以看出,当仅x方向移动质心时,x方向质心误差补偿后相对于补偿前大幅度减小,y方向质心误差略有增大,是由于补偿参数计算结果因灰度噪声影响,在修正质心时引入了一个弱小的随机噪声,对于不需要精度补偿的方向精度会稍微变差。
仅移动y方向,x方向质心保持不变,移动范围为从一个像元中心移动到下一个像元中心,星点像斑测量误差如图3所示,测量结果如表2所示:
表2 y方向移动质心测量结果
Figure BDA0001832411740000081
依据表2可以看出,当仅y方向移动质心时,y方向质心误差补偿后相对于补偿前大幅度减小,x方向质心误差略有增大。
以相同的速率同时移动x、y方向,移动范围为从一个像元中心移动到下一个像元中心,星点像斑测量误差如图4所示,测量结果如表3所示:
表3 x、y方向同时移动质心测量结果
Figure BDA0001832411740000091
依据表3可以看出,当仅x、y方向移动质心时,方向质心误差补偿后相对于补偿前大幅度减小。
综上所述,本发明方法可以有效补偿传统质心定位精度,减小质心定位误差。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种星敏感器质心精度自适应补偿方法,其特征在于,包括:
步骤一,确定星点窗口区域;设星点像斑中灰度最大的像元位置为(x0,y0),以(x0,y0)为中心,以n为半径,选取星点的窗口区域记为W,则窗口区域W内的任意像素位置(x,y)满足:
Figure FDA0003528921060000011
步骤二,计算窗口区域质心;窗口区域内位置(x,y)的像素灰度为g(x,y),计算x方向窗口区域灰度累加和:
Figure FDA0003528921060000012
计算y方向窗口区域灰度累加和:
Figure FDA0003528921060000013
计算窗口区域x方向灰度质心:
Figure FDA0003528921060000014
计算窗口区域y方向灰度质心:
Figure FDA0003528921060000015
步骤三,计算自适应补偿参数;
计算x方向窗口区域质心补偿参数:
Figure FDA0003528921060000021
计算y方向窗口区域质心补偿参数:
Figure FDA0003528921060000022
其中k为修正参数,跟星敏感器的星点像斑成像质量相关,需要通过标定测出;
步骤四,修正窗口区域质心;
修正窗口区域x方向质心:
x’zx=xzx+Δxzx (8);
修正窗口区域y方向质心:
y’zx=yzx+Δyzx (9);
所述星敏感器质心精度自适应补偿方法还包括:星点像斑灰度最高像元窗口区域确定模块、窗口区域质心计算模块和自适应补偿参数计算模块;
所述星点像斑灰度最高像元窗口区域确定模块、所述窗口区域质心计算模块和所述自适应补偿参数计算模块顺次连通;所述星点像斑灰度最高像元窗口区域确定模块确定灰度值最高的像元及其窗口区域,窗口区域质心计算模块分别计算窗口区域x方向、y方向质心,所述自适应补偿参数计算模块依据窗口中心邻域的灰度比例关系计算自适应补偿参数并修正各自方向质心,从而完成星点像斑质心定位及补偿。
CN201811210853.2A 2018-10-17 2018-10-17 一种星敏感器质心精度自适应补偿方法 Active CN109443381B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811210853.2A CN109443381B (zh) 2018-10-17 2018-10-17 一种星敏感器质心精度自适应补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811210853.2A CN109443381B (zh) 2018-10-17 2018-10-17 一种星敏感器质心精度自适应补偿方法

Publications (2)

Publication Number Publication Date
CN109443381A CN109443381A (zh) 2019-03-08
CN109443381B true CN109443381B (zh) 2022-05-20

Family

ID=65547295

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811210853.2A Active CN109443381B (zh) 2018-10-17 2018-10-17 一种星敏感器质心精度自适应补偿方法

Country Status (1)

Country Link
CN (1) CN109443381B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111105446B (zh) * 2019-11-18 2023-04-18 上海航天控制技术研究所 一种星点提取与补偿方法
CN110849354B (zh) * 2019-11-28 2021-12-10 上海航天控制技术研究所 一种星敏感器寿命末期条件下星点提取与补偿方法
CN113514054A (zh) * 2021-06-16 2021-10-19 北京遥感设备研究所 一种星敏感器星点像斑检测方法及系统
CN114543790B (zh) * 2021-12-30 2024-05-10 浙江众星志连科技有限责任公司 一种微小卫星太阳敏感器的鬼影去除方法及系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5982940A (en) * 1995-11-01 1999-11-09 Minolta Co., Ltd. Image data processing device and method of processing image data
CN102496015A (zh) * 2011-11-22 2012-06-13 南京航空航天大学 二维高斯分布光斑图像中心快速高精度定位方法
CN102564457A (zh) * 2011-12-29 2012-07-11 北京控制工程研究所 一种aps星敏感器在轨噪声自主抑制方法
CN103323027A (zh) * 2013-05-30 2013-09-25 北京控制工程研究所 基于星点重构的星敏感器动态补偿方法
CN103323025A (zh) * 2013-05-24 2013-09-25 哈尔滨工业大学 一种太阳光斑质心及其太阳输出角的误差修正方法
CN105023281A (zh) * 2015-05-25 2015-11-04 中国科学院空间科学与应用研究中心 基于点扩散函数波前修正的星点像质心计算方法
CN105761288A (zh) * 2016-02-02 2016-07-13 华中科技大学 一种基于fpga的实时星点质心定位方法及装置
CN108168541A (zh) * 2017-12-20 2018-06-15 北京遥感设备研究所 一种改进的亚像元星点质心定位方法
CN108507564A (zh) * 2018-06-29 2018-09-07 中国科学院国家空间科学中心 一种基于点扩散函数拟合的星敏感器质心定位方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5982940A (en) * 1995-11-01 1999-11-09 Minolta Co., Ltd. Image data processing device and method of processing image data
CN102496015A (zh) * 2011-11-22 2012-06-13 南京航空航天大学 二维高斯分布光斑图像中心快速高精度定位方法
CN102564457A (zh) * 2011-12-29 2012-07-11 北京控制工程研究所 一种aps星敏感器在轨噪声自主抑制方法
CN103323025A (zh) * 2013-05-24 2013-09-25 哈尔滨工业大学 一种太阳光斑质心及其太阳输出角的误差修正方法
CN103323027A (zh) * 2013-05-30 2013-09-25 北京控制工程研究所 基于星点重构的星敏感器动态补偿方法
CN105023281A (zh) * 2015-05-25 2015-11-04 中国科学院空间科学与应用研究中心 基于点扩散函数波前修正的星点像质心计算方法
CN105761288A (zh) * 2016-02-02 2016-07-13 华中科技大学 一种基于fpga的实时星点质心定位方法及装置
CN108168541A (zh) * 2017-12-20 2018-06-15 北京遥感设备研究所 一种改进的亚像元星点质心定位方法
CN108507564A (zh) * 2018-06-29 2018-09-07 中国科学院国家空间科学中心 一种基于点扩散函数拟合的星敏感器质心定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"星点质心亚像元定位的高精度误差补偿法";杨君 等;《光学精密工程》;20100430;第18卷(第4期);第1002-1010页 *

Also Published As

Publication number Publication date
CN109443381A (zh) 2019-03-08

Similar Documents

Publication Publication Date Title
CN109443381B (zh) 一种星敏感器质心精度自适应补偿方法
CN108168541B (zh) 一种改进的亚像元星点质心定位方法
CN108198219B (zh) 用于摄影测量的相机标定参数的误差补偿方法
CN108871373B (zh) 一种基于俯仰滚转台和非线性优化的星敏感器标定方法
CN110260822B (zh) 一种多目结构光系统高精度标定方法
CN110490932B (zh) 单目红外共面光标迭代优化测量行车吊臂空间位姿方法
CN109285195B (zh) 基于大尺寸靶标的单目投影系统逐像素畸变校正方法及其应用
CN106023237B (zh) 一种双目摄像机定位校准方法
US7155836B2 (en) Method and apparatus for using magnetic field
CN113776455B (zh) 一种离轴非球面反射镜零位补偿检测非线性误差校正方法
CN109099857B (zh) 一种基于surf特征匹配的子孔径拼接方法
CN109798874B (zh) 一种高铁桥梁竖向动扰度测量方法
CN115236828B (zh) 一种基于离散正交多项式的共相位误差校正方法
CN110751601A (zh) 一种基于rc光学系统的畸变校正方法
Ding et al. A robust detection method of control points for calibration and measurement with defocused images
CN110595374B (zh) 基于传递像机的大型结构件实时形变监测方法
CN110211175B (zh) 准直激光器光束空间位姿标定方法
CN109813510B (zh) 基于无人机的高铁桥梁竖向动扰度测量方法
CN110490941B (zh) 一种基于法向量的远心镜头外参数标定方法
CN111664801A (zh) 一种基于高斯曲率的光学靶标点投影规划方法
CN111044039A (zh) 基于imu的单目目标区域自适应高精度测距装置和方法
CN114608540A (zh) 一种数字摄影测量系统的测量网型确定方法
CN114140444A (zh) 检测运动平台旋转精度的方法
JP2018125706A (ja) 撮像装置
CN117173255B (zh) 多相机位置的标定方法、系统及晶圆键合检测设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant