CN109443381B - 一种星敏感器质心精度自适应补偿方法 - Google Patents
一种星敏感器质心精度自适应补偿方法 Download PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- 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
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)满足:
优选地,所述步骤2中计算窗口区域质心,
窗口区域内位置(x,y)的像素灰度为g(x,y),计算x方向窗口区域灰度累加和:
计算y方向窗口区域灰度累加和:
优选地,所述步骤2中计算窗口区域质心,
计算窗口区域x方向灰度质心:
计算窗口区域y方向灰度质心:
优选地,所述步骤3中计算自适应补偿参数,
计算x方向窗口质心补偿参数:
其中k为修正参数,跟星敏感器的星点像斑成像质量相关,需要通过标定测出。
优选地,所述步骤3中计算自适应补偿参数,
计算y方向窗口质心补偿参数:
其中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)满足:
第二步,计算窗口区域质心。
窗口区域内位置(x,y)的像素灰度为g(x,y),计算x方向窗口区域灰度累加和:
计算y方向窗口区域灰度累加和:
计算窗口区域x方向灰度质心:
计算窗口区域y方向灰度质心:
第三步,计算自适应补偿参数。
计算x方向窗口质心补偿参数:
计算y方向窗口质心补偿参数:
其中k为修正参数,跟星敏感器的星点像斑成像质量相关,需要通过标定测出。
第四步,修正窗口区域质心。
修正窗口区域x方向质心:
x′zx=xzx+Δxzx (8)
修正窗口区域y方向质心:
y′zx=yzx+Δyzx (9)
一种星敏感器质心精度自适应补偿方法,包括:星点像斑灰度最高像元窗口区域确定模块、窗口区域质心计算模块和自适应补偿参数计算模块;星点像斑灰度最高像元窗口区域确定模块、窗口区域质心计算模块和自适应补偿参数计算模块顺次连通;其中星点像斑灰度最高像元窗口区域确定模块确定灰度值最高的像元及其窗口区域,窗口区域质心计算模块分别计算窗口区域x方向、y方向质心,自适应补偿参数计算模块依据窗口中心邻域的灰度比例关系计算自适应补偿参数并修正各自方向质心,从而完成星点像斑质心定位及补偿,提高质心定位精度。
利用本方法对某星敏感器进行星点像斑质心定位测试,测试结果为:
仅移动x方向,y方向质心保持不变,移动范围为从一个像元中心移动到下一个像元中心,星点像斑测量误差如图2所示,测量结果如表1所示:
表1 x方向移动质心测量结果
依据表1可以看出,当仅x方向移动质心时,x方向质心误差补偿后相对于补偿前大幅度减小,y方向质心误差略有增大,是由于补偿参数计算结果因灰度噪声影响,在修正质心时引入了一个弱小的随机噪声,对于不需要精度补偿的方向精度会稍微变差。
仅移动y方向,x方向质心保持不变,移动范围为从一个像元中心移动到下一个像元中心,星点像斑测量误差如图3所示,测量结果如表2所示:
表2 y方向移动质心测量结果
依据表2可以看出,当仅y方向移动质心时,y方向质心误差补偿后相对于补偿前大幅度减小,x方向质心误差略有增大。
以相同的速率同时移动x、y方向,移动范围为从一个像元中心移动到下一个像元中心,星点像斑测量误差如图4所示,测量结果如表3所示:
表3 x、y方向同时移动质心测量结果
依据表3可以看出,当仅x、y方向移动质心时,方向质心误差补偿后相对于补偿前大幅度减小。
综上所述,本发明方法可以有效补偿传统质心定位精度,减小质心定位误差。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种星敏感器质心精度自适应补偿方法,其特征在于,包括:
步骤一,确定星点窗口区域;设星点像斑中灰度最大的像元位置为(x0,y0),以(x0,y0)为中心,以n为半径,选取星点的窗口区域记为W,则窗口区域W内的任意像素位置(x,y)满足:
步骤二,计算窗口区域质心;窗口区域内位置(x,y)的像素灰度为g(x,y),计算x方向窗口区域灰度累加和:
计算y方向窗口区域灰度累加和:
计算窗口区域x方向灰度质心:
计算窗口区域y方向灰度质心:
步骤三,计算自适应补偿参数;
计算x方向窗口区域质心补偿参数:
计算y方向窗口区域质心补偿参数:
其中k为修正参数,跟星敏感器的星点像斑成像质量相关,需要通过标定测出;
步骤四,修正窗口区域质心;
修正窗口区域x方向质心:
x’zx=xzx+Δxzx (8);
修正窗口区域y方向质心:
y’zx=yzx+Δyzx (9);
所述星敏感器质心精度自适应补偿方法还包括:星点像斑灰度最高像元窗口区域确定模块、窗口区域质心计算模块和自适应补偿参数计算模块;
所述星点像斑灰度最高像元窗口区域确定模块、所述窗口区域质心计算模块和所述自适应补偿参数计算模块顺次连通;所述星点像斑灰度最高像元窗口区域确定模块确定灰度值最高的像元及其窗口区域,窗口区域质心计算模块分别计算窗口区域x方向、y方向质心,所述自适应补偿参数计算模块依据窗口中心邻域的灰度比例关系计算自适应补偿参数并修正各自方向质心,从而完成星点像斑质心定位及补偿。
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)
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)
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 | 中国科学院国家空间科学中心 | 一种基于点扩散函数拟合的星敏感器质心定位方法 |
-
2018
- 2018-10-17 CN CN201811210853.2A patent/CN109443381B/zh active Active
Patent Citations (9)
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)
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 |