CN101363732A - 一种高帧频太阳敏感器及其实现方法 - Google Patents

一种高帧频太阳敏感器及其实现方法 Download PDF

Info

Publication number
CN101363732A
CN101363732A CNA2008102224887A CN200810222488A CN101363732A CN 101363732 A CN101363732 A CN 101363732A CN A2008102224887 A CNA2008102224887 A CN A2008102224887A CN 200810222488 A CN200810222488 A CN 200810222488A CN 101363732 A CN101363732 A CN 101363732A
Authority
CN
China
Prior art keywords
hot spot
image
coordinate
center
solar facula
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
Application number
CNA2008102224887A
Other languages
English (en)
Other versions
CN101363732B (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.)
Beihang University
Beijing University of Aeronautics and Astronautics
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN200810222488A priority Critical patent/CN101363732B/zh
Publication of CN101363732A publication Critical patent/CN101363732A/zh
Application granted granted Critical
Publication of CN101363732B publication Critical patent/CN101363732B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

本发明公开了一种高帧频太阳敏感器,该敏感器通过图像采集模块、计算模块、预测模块和第一控制模块,实现了基于ROI窗口动态跟踪的高数据更新率。本发明还公开了一种高帧频太阳敏感器实现方法,该方法通过利用当前太阳光斑图像,以及前一太阳光斑图像得到的同一小孔对应的光斑的质心坐标,预测下一太阳光斑图像中该小孔对应的光斑所在的位置,来实现ROI窗口动态跟踪,与基于全帧图像采集的太阳敏感器相比,大大减小了背景像素区域的数据。采用本发明所述的高帧频太阳敏感器及其实现方法,可减少获得姿态信息所需处理的数据量,提高太阳敏感器数据更新率。

Description

一种高帧频太阳敏感器及其实现方法
技术领域
本发明涉及姿态敏感器技术,尤其涉及一种高帧频太阳敏感器及其实现方法。
背景技术
太阳敏感器是姿态敏感器的一种,它通过测量太阳光线与卫星某一体轴之间的夹角,确定太阳在太阳敏感器本体坐标系中的位置。目前,高精度太阳敏感器所用的光敏元件主要有电荷耦合器件(CCD)图像传感器、互补金属氧化物半导体(CMOS)图像传感器两种。随着卫星微小型化的趋势,太阳敏感器作为卫星的组成部分也需微小型化,基于CCD的太阳敏感器难以满足卫星组成部件微小型化的这一需求,而CMOS图像传感器凭借其低成本、低功耗、小体积、高集成度、控制灵活等优点,得到了越来越广泛的应用。
但是,目前基于CMOS图像传感器的太阳敏感器是通过全帧图像采集的方式来获取太阳光斑图像,并通过对太阳光斑图像进行全帧数据处理确定太阳在太阳敏感器本体坐标系中的位置,而对太阳光斑图像进行全帧数据处理,所需处理的数据量大、处理时间长,使得太阳敏感器的数据更新率较低,最高的也不过10Hz。
发明内容
有鉴于此,本发明的主要目的在于提供一种高帧频太阳敏感器及其实现方法,可减少获得姿态信息所需处理的数据量,进而提高太阳敏感器的数据更新率。
为达到上述目的,本发明的技术方案是这样实现的:
一种高帧频太阳敏感器,该敏感器包括:
光学系统,用于形成太阳光斑;
图像采集模块,用于通过全帧图像采集得到太阳光斑图像,还用于采集ROI窗口的太阳光斑图像;
计算模块,用于计算所述太阳光斑图像内的各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角;还用于计算所述ROI窗口的太阳光斑图像内各光斑的质心坐标,根据光斑的质心坐标计算太阳姿态角并输出;
预测模块,用于根据计算模块所得到的当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑的质心坐标,预测该小孔对应的光斑在下一太阳光斑图像中的位置;
第一控制模块,用于设置ROI窗口大小,以及根据所预测出的位置和预先设置的感兴趣区域ROI窗口大小对所述图像采集模块进行ROI窗口的配置,还用于发送控制信息,命令图像采集模块执行ROI窗口的太阳光斑图像采集。
进一步地,该敏感器进一步包括:
计数模块,用于对全帧采集的图像内的光斑进行计数,得出全帧采集的图像内的光斑数目并存储,并用于对所述ROI窗口的太阳光斑图像内的光斑进行计数,得出数目所述ROI窗口的太阳光斑图像内的光斑数目;
判断模块,用于判断所述ROI窗口的太阳光斑图像内的光斑数目与全帧采集的图像所获得的光斑数目是否相等,并输出判断结果;
第二控制模块,用于根据判断模块输出的结果控制图像采集模块返回执行全帧图像采集或控制计算模块计算ROI窗口内的各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角。
进一步地,所述计算模块包括:
质心计算模块,用于利用流水线结构,计算所述太阳光斑图像内的各光斑的质心坐标;
姿态角计算模块,用于根据光斑的质心坐标计算太阳姿态角。
进一步地,所述预测模块包括:
位移计算模块,用于计算同一小孔对应的光斑在x轴方向的位移,还用于计算同一小孔对应的光斑在y轴方向的位移;
位置确定模块,用于确定所述小孔对应的光斑在下一太阳光斑图像中的位置。
本发明还提供一种高帧频太阳敏感器实现方法,该方法包括以下步骤:
A、全帧图像采集,获得太阳光斑图像,计算所获得的太阳光斑图像的各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角;
B、重复执行步骤A至少一次后,执行步骤C;
C、根据得到的当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑的质心坐标,预测该小孔对应的光斑在下一太阳光斑图像中的位置,并根据所预测出的位置和预设的ROI窗口大小进行ROI窗口的配置;
D、采集ROI窗口的太阳光斑图像,计算所采集到的ROI窗口的太阳光斑图像内各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角,返回步骤C。
进一步地,所述采集ROI窗口的太阳光斑图像后计算质心坐标前,步骤D还包括:
判断所述ROI窗口的太阳光斑图像内的光斑数目与全帧采集的图像所获得的光斑数目是否相等,如果相等,计算ROI窗口内各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角;如果不相等,返回步骤A。
进一步地,所述计算太阳光斑图像内的各光斑质心坐标包括:
a、在当前时钟周期内,计算再前一像素灰度值和坐标值的乘积与已处理的同一光斑所有像素灰度值和坐标值乘积的累加值;识别前一像素所属光斑,将前一像素作为再前一像素;并读取当前像素作为前一像素缓存;
b、重复步骤a,直到处理完整个所述太阳光斑图像后,分别计算所的到的各个光斑灰度值和坐标值乘积的累加值与灰度累加值之商,并将所得到商作为对应光斑的质心的x轴坐标,以及计算所得到的各个光斑灰度值和y轴坐标值乘积的累加值与灰度累加值之商,并将所得到的商作为对应光斑的质心的y轴坐标。
进一步地,所述根据得到的同一小孔对应的光斑质心坐标,预测该小孔对应的光斑在下一太阳光斑图像中的位置包括:
将所述当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑质心的x轴坐标之差作为光斑在x轴方向的位移,将所述当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑质心的y轴坐标之差作为光斑在y轴方向的位移;
将当前太阳光斑图像中所述小孔对应的光斑质心的x轴坐标与所述光斑在x轴方向的位移求和,将当前太阳光斑图像中所述小孔对应的光斑质心的y轴坐标与所述光斑在y轴方向的位移求和,即得到该小孔对应的光斑在下一太阳光斑图像中的位置。
进一步地,所述计算所采集到的ROI窗口的太阳光斑图像内各光斑的质心坐标包括:
a1、在当前时钟周期内,计算再前一像素灰度值和坐标值的乘积与已处理的同一光斑所有像素灰度值和坐标值乘积的累加值;识别前一像素所属光斑,将前一像素作为再前一像素;并读取当前像素并作为前一像素缓存;
b1、重复步骤a1,直到处理完整个ROI窗口区域后,分别计算所的到的各个光斑灰度值和坐标值乘积的累加值与灰度累加值之商,并将所得到商作为对应光斑的质心的x轴坐标,以及计算所得到的各个光斑灰度值和y轴坐标值乘积的累加值与灰度累加值之商,并将所得到的商作为对应光斑的质心的y轴坐标。
本发明所提供的高帧频太阳敏感器及其实现方法,对感兴趣区域(ROI)窗口进行动态跟踪,并采集ROI窗口的太阳光斑图像,与基于全帧图像采集所获得的太阳光斑图像相比,可大大减小背景像素区域,这样,读取像素时,只读取包含太阳光斑的区域,进而可减小读取数据和后续质心坐标计算、太阳姿态角计算的数据量和处理时间,进而提高太阳敏感器的帧频,即太阳敏感器的数据更新率。
本发明中的质心坐标计算引入了流水线结构,使得针对每个像素,其累加计算比识别像素所属光斑滞后一个时钟周期,识别像素所属光斑比读取像素滞后一个时钟周期,进而可在一个时钟周期内同时进行再前一像素的累加计算、前一像素所属光斑的识别、以及当前像素的读取,如此大大减小了时钟周期,将系统的最高可运行时钟频率由原来的43MHz左右到提高到65MHz左右。
附图说明
图1为本发明高帧频太阳敏感器的组成结构示意图;
图2为本发明高帧频太阳敏感器实现方法的流程示意图;
图3为本发明ROI窗口跟踪示意图;
图4为本发明质心坐标的计算流程示意图。
具体实施方式
本发明的基本思路为:利用当前太阳光斑图像,以及前一太阳光斑图像得到的同一小孔对应的光斑质心坐标,预测下一太阳光斑图像中该小孔对应的光斑所在的位置,实现对包含光斑的ROI窗口的动态跟踪,然后进行ROI窗口的太阳光斑图像的采集,计算出ROI窗口的太阳光斑图像内各光斑的质心坐标,并根据光斑的质心坐标计算出太阳姿态角;减小了读取数据和后续质心坐标计算、太阳姿态角计算的数据量和处理时间,进而提高太阳敏感器的帧频,即太阳敏感器的数据更新率。
本发明高帧频太阳敏感器的组成结构如图1所示,包括:
光学系统,包括透光的小孔阵列,用于形成太阳光斑;
图像采集模块,用于通过全帧图像采集得到太阳光斑图像,还用于采集ROI窗口的太阳光斑图像;所述图像采集模块通过Cypress公司的CMOS图像传感器IBIS5-A1300实现;
计算模块,用于计算所述太阳光斑图像内的各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角;还用于计算所述ROI窗口的太阳光斑图像内各光斑的质心坐标,根据光斑的质心坐标计算太阳姿态角并输出;所述计算模块通过现场可编程门阵列(FPGA)实现;
预测模块,用于根据计算模块所得到的当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑的质心坐标,预测该小孔对应的光斑在下一太阳光斑图像中的位置;所述预测模块通过FPGA实现;
第一控制模块,用于设置ROI窗口大小,以及根据所预测出的位置和预先设置的感兴趣区域ROI窗口大小进行ROI窗口的配置,还用于在完成ROI窗口的配置后,发送控制信息,命令图像采集模块由执行全帧图像采集转为执行ROI窗口的太阳光斑图像采集;所述第一控制模块通过FPGA实现。
该敏感器进一步包括:
计数模块,用于对全帧采集的图像内的光斑进行计数,得出全帧采集的图像内的光斑数目并存储,并用于对所述ROI窗口的太阳光斑图像内的光斑进行计数,得出数目所述ROI窗口的太阳光斑图像内的光斑数目;
判断模块,用于判断所述ROI窗口的太阳光斑图像内的光斑数目与所存储的全帧采集的图像所获得的光斑数目是否相等,并输出判断结果;
第二控制模块,用于根据判断模块输出的结果控制图像采集模块返回执行全帧图像采集或控制计算模块计算ROI窗口内的各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角。
所述计算模块包括:
质心计算模块,用于利用流水线结构,计算所述太阳光斑图像内的各光斑的质心坐标;所述流水线结构指在每个时钟周期内,计算再前一像素灰度值和坐标值的乘积与已处理的同一光斑所有像素灰度值和坐标值乘积的累加值;识别前一像素所属光斑,将前一像素作为再前一像素;并读取当前像素作为前一像素缓存;
姿态角计算模块,用于根据光斑的质心坐标计算太阳姿态角。
所述预测模块包括:
位移计算模块,用于计算同一小孔对应的光斑在x轴方向的位移,还用于计算同一小孔对应的光斑在y轴方向的位移;
所述光斑在x轴方向的位移等于所述当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑质心的x轴坐标之差;所述光斑在y轴方向的位移等于所述当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑质心的y轴坐标之差;
位置确定模块,用于确定所述小孔对应的光斑在下一太阳光斑图像中的位置;
所述位置的x轴坐标等于当前太阳光斑图像中所述小孔对应的光斑质心的x轴坐标与所述光斑在x轴方向的位移之和;所述位置的y轴坐标等于当前太阳光斑图像中所述小孔对应的光斑质心的y轴坐标与所述光斑在y轴方向的位移之和。
本发明还提出一种高帧频太阳敏感器的实现方法,该方法的实现流程如图2所示,包括以下步骤:
步骤201:设置ROI窗口的大小;
这里,ROI窗口大小的具体设置方法包括:
a、根据光斑的直径和光斑的间距,计算出可以完全包含全帧图像采集所获得的太阳光斑图像中的所有光斑的最小正方形的边长;所述正方形的上下两边平行于x轴;所述光斑的间距指相邻两光斑的质心的距离;
b、在所述边长的基础上增加一定余量,得出ROI窗口的边长;且第n行的第m个小孔对应的光斑的质心距ROI窗口左边和上边的距离分别为:
Figure A200810222488D00111
Figure A200810222488D00112
其中,S为第n行的第m个小孔对应的光斑的质心距ROI窗口左边的距离,S为第n行的第m个小孔对应的光斑的质心距ROI窗口上边的距离,d为光斑的直径,s为光斑的间距,l为余量。
可见,ROI窗口为包含光斑的区域,与全帧图像采集所获得的太阳光斑图像相比,大大减小了背景像素区域,这样,在读取像素时,只需读取包含光斑的区域,因此读取数据的时间可以大大缩短;并且,能够减小后续处理的数据量和处理时间,进而提高太阳敏感器的帧频。
通常,ROI窗口越小,所述余量越小,帧频越高;但由于本发明中所述位置的预测是基于太阳敏感器做匀速运动的假设进行的,而实际应用中,只是基于两次图像采集间隔时间短,可近似认为太阳敏感器做匀速运动,预测的位置与实际位置有一定偏差,如果余量太小,会增大跟踪失败率,因此,较佳的余量设置为77像素。
步骤202:通过全帧图像采集得到太阳光斑图像,计算所获得的太阳光斑图像内各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角;
其中,具体如何根据光斑的质心坐标计算太阳姿态角的具体实现方法,可参见专利号为200610103797.3,发明名称为“高精度APS太阳敏感器”的中国专利,这里不再赘述;所述全帧图像采集可通过Cypress公司的CMOS图像传感器IBIS5-A1300实现;所述光斑的质心坐标计算和太阳姿态角计算均可通过FPGA实现。
具体的,本发明中计算质心坐标的方法包括:
a1、在本次质心坐标计算的第一个时钟周期内,读取当前像素并作为前一像素缓存;
b1、在本次质心坐标计算的第二个时钟周期内,将缓存器内的前一像素作为再前一像素,识别该再前一像素所属光斑,并读取当前像素作为前一像素缓存;
从本次质心坐标计算的第三个时钟周期开始的每个时钟周期,重复执行一次步骤c1;直到处理完整个所述太阳光斑图像后,再分别计算所得到的各个光斑灰度值和x轴坐标值乘积的累加值与灰度累加值之商,并将所得到的商作为对应光斑的质心的x轴坐标,以及计算所得到的各个光斑灰度值和y轴坐标值乘积的累加值与灰度累加值之商,并将所得到的商作为对应光斑的质心的y轴坐标;
c1、计算再前一像素灰度值和坐标值的乘积与已处理的同一光斑所有像素灰度值和坐标值乘积的累加值;识别前一像素所属光斑,将前一像素作为再前一像素;并读取当前像素作为前一像素缓存。其中,具体如何读取像素、识别像素所属光斑、以及累加计算的可参见专利号为200610161802.6,发明名称为“快速高精度太阳光斑图像质心定位方法及装置”的中国专利,这里不再赘述。
步骤203:再次执行步骤202后,执行步骤204。
步骤204:根据得到的当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑的质心坐标,预测该小孔对应的光斑在下一帧太阳光斑图像中的位置;并根据所预测出的位置和所设置的ROI窗口的大小,进行ROI窗口的配置;
这里,太阳光斑图像指全帧图像采集所得到太阳光斑图像或ROI窗口的太阳光斑图像。
如图3中ROI窗口跟踪示意图所示,本步骤中具体的预测方法包括:
a2、将所述当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑质心的x轴坐标之差作为光斑在x轴方向的位移,将所述当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑质心的y轴坐标之差作为光斑在y轴方向的位移;其计算式表示为:
Δx=roi_x_now-roi_x_old
Δy=roi_x_now-roi_y_old
其中,Δx和Δy分别为光斑在x轴和y轴方向的位移,roi_x_now为当前太阳光斑图像中所述小孔对应的光斑质心的x轴坐标,roi_y_now为当前太阳光斑图像中所述小孔对应的光斑质心的y轴坐标,roi_x_old为前一太阳光斑图像中所述小孔对应的光斑质心的x轴坐标,roi_y_old为前一太阳光斑图像中所述小孔对应的光斑质心的y轴坐标;
b2、将当前太阳光斑图像中所述小孔对应的光斑质心的x轴坐标与所述光斑在x轴方向的位移Δx求和,将当前太阳光斑图像中所述小孔对应的光斑质心的y轴坐标与所述光斑在y轴方向的位移Δy求和,得到该小孔对应的光斑在下一太阳光斑图像中的位置;其计算式表示为:
roi_x_next=roi_x_now+Δx
roi_y_next=roi_x_now+Δy
其中,(roi_x_next,roi_y_next)为该小孔对应的光斑在下一太阳光斑图像中的位置坐标。
步骤205:采集ROI窗口的太阳光斑图像。
步骤206:判断所采集到的ROI窗口的太阳光斑图像内的光斑数目与全帧采集的图像所获得的光斑数目是否相等,如果相等,执行步骤207;如果不相等,认为跟踪失败,返回步骤202。
步骤207:计算ROI窗口内各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角,返回步骤204;其中,具体如何根据光斑的质心坐标计算太阳姿态角的具体实现方法可参见专利号为200610103797.3,发明名称为“高精度APS太阳敏感器”的中国专利,这里不再赘述;所述光斑的质心坐标计算和太阳姿态角计算均可通过FPGA实现。
这里,具体计算质心坐标的方法同步骤202中计算质心坐标的方法。
上述流程中,步骤203、204根据前几帧得到的同一小孔对应的光斑的质心坐标,预测下一帧中该小孔对应的光斑所在的位置,实现了ROI窗口动态跟踪。在实际应用中,该ROI窗口动态跟踪可通过Cypress公司的CMOS图像传感器IBIS5-A1300实现,IBIS5-A1300具有高分辨率和双快门模式,最大分辨率可达1280*1024,帧频可达到27帧/s,ROI模式下,帧频可达100帧/s以上。
上述流程中,步骤201、205中质心坐标计算和太阳姿态角计算均可通过FPGA硬件实现,保证在高帧频下能实时输出光斑的质心坐标和太阳姿态角。
在步骤201、205中,质心坐标计算引入了流水线结构,即:使得针对每个像素,其累加计算比识别像素所属光斑滞后一个时钟周期,识别像素所属光斑比读取像素滞后一个时钟周期,这样就实现了如图4所示的质心坐标的计算流程:在一个时钟周期内可同时进行再前一像素的累加计算、前一像素所属光斑的识别、以及当前像素的读取,使得一个时钟周期只需不小于数据读入所需时间、标志判断所需时间和数据累加所需时间中的最大时间即可,大大减小了时钟周期,将系统的最高可运行时钟频率由原来的43MHz左右到提高到65MHz左右,解决了现有质心坐标的计算方法在一个时钟周期内需要针对同一像素顺序完成读取像素、识别像素所属光斑、累加计算三个步骤,会导致时序紧张,系统可运行的最高时钟频率低的问题。
以上所述,仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。

Claims (9)

1、一种高帧频太阳敏感器,其特征在于,该敏感器包括:
光学系统,用于形成太阳光斑;
图像采集模块,用于通过全帧图像采集得到太阳光斑图像,还用于采集ROI窗口的太阳光斑图像;
计算模块,用于计算所述太阳光斑图像内的各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角;还用于计算所述ROI窗口的太阳光斑图像内各光斑的质心坐标,根据光斑的质心坐标计算太阳姿态角并输出;
预测模块,用于根据计算模块所得到的当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑的质心坐标,预测该小孔对应的光斑在下一太阳光斑图像中的位置;
第一控制模块,用于设置ROI窗口大小,以及根据所预测出的位置和预先设置的感兴趣区域ROI窗口大小对所述图像采集模块进行ROI窗口的配置,还用于发送控制信息,命令图像采集模块执行ROI窗口的太阳光斑图像采集。
2、根据权利要求1所述高帧频太阳敏感器,其特征在于,该敏感器进一步包括:
计数模块,用于对全帧采集的图像内的光斑进行计数,得出全帧采集的图像内的光斑数目并存储,并用于对所述ROI窗口的太阳光斑图像内的光斑进行计数,得出数目所述ROI窗口的太阳光斑图像内的光斑数目;
判断模块,用于判断所述ROI窗口的太阳光斑图像内的光斑数目与全帧采集的图像所获得的光斑数目是否相等,并输出判断结果;
第二控制模块,用于根据判断模块输出的结果控制图像采集模块返回执行全帧图像采集或控制计算模块计算ROI窗口内的各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角。
3、根据权利要求1或2所述高帧频太阳敏感器,其特征在于,所述计算模块包括:
质心计算模块,用于利用流水线结构,计算所述太阳光斑图像内的各光斑的质心坐标;
姿态角计算模块,用于根据光斑的质心坐标计算太阳姿态角。
4、根据权利要求1或2所述高帧频太阳敏感器,其特征在于,所述预测模块包括:
位移计算模块,用于计算同一小孔对应的光斑在x轴方向的位移,还用于计算同一小孔对应的光斑在y轴方向的位移;
位置确定模块,用于确定所述小孔对应的光斑在下一太阳光斑图像中的位置。
5、一种高帧频太阳敏感器实现方法,其特征在于,该方法包括以下步骤:
A、全帧图像采集,获得太阳光斑图像,计算所获得的太阳光斑图像的各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角;
B、重复执行步骤A至少一次后,执行步骤C;
C、根据得到的当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑的质心坐标,预测该小孔对应的光斑在下一太阳光斑图像中的位置,并根据所预测出的位置和预设的ROI窗口大小进行ROI窗口的配置;
D、采集ROI窗口的太阳光斑图像,计算所采集到的ROI窗口的太阳光斑图像内各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角,返回步骤C。
6、根据权利要求5所述高帧频太阳敏感器实现方法,其特征在于,所述采集ROI窗口的太阳光斑图像后计算质心坐标前,步骤D还包括:
判断所述ROI窗口的太阳光斑图像内的光斑数目与全帧采集的图像所获得的光斑数目是否相等,如果相等,计算ROI窗口内各光斑的质心坐标,并根据光斑的质心坐标计算太阳姿态角;如果不相等,返回步骤A。
7、根据权利要求5或6所述高帧频太阳敏感器实现方法,其特征在于,所述计算太阳光斑图像内的各光斑质心坐标包括:
a、在当前时钟周期内,计算再前一像素灰度值和坐标值的乘积与已处理的同一光斑所有像素灰度值和坐标值乘积的累加值;识别前一像素所属光斑,将前一像素作为再前一像素;并读取当前像素作为前一像素缓存;
b、重复步骤a,直到处理完整个所述太阳光斑图像后,分别计算所的到的各个光斑灰度值和坐标值乘积的累加值与灰度累加值之商,并将所得到商作为对应光斑的质心的x轴坐标,以及计算所得到的各个光斑灰度值和y轴坐标值乘积的累加值与灰度累加值之商,并将所得到的商作为对应光斑的质心的y轴坐标。
8、根据权利要求5或6所述高帧频太阳敏感器实现方法,其特征在于,所述根据得到的同一小孔对应的光斑质心坐标,预测该小孔对应的光斑在下一太阳光斑图像中的位置包括:
将所述当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑质心的x轴坐标之差作为光斑在x轴方向的位移,将所述当前太阳光斑图像,以及前一太阳光斑图像中同一小孔对应的光斑质心的y轴坐标之差作为光斑在y轴方向的位移;
将当前太阳光斑图像中所述小孔对应的光斑质心的x轴坐标与所述光斑在x轴方向的位移求和,将当前太阳光斑图像中所述小孔对应的光斑质心的y轴坐标与所述光斑在y轴方向的位移求和,即得到该小孔对应的光斑在下一太阳光斑图像中的位置。
9、根据权利要求5或6所述高帧频太阳敏感器实现方法,其特征在于,所述计算所采集到的ROI窗口的太阳光斑图像内各光斑的质心坐标包括:
a1、在当前时钟周期内,计算再前一像素灰度值和坐标值的乘积与已处理的同一光斑所有像素灰度值和坐标值乘积的累加值;识别前一像素所属光斑,将前一像素作为再前一像素;并读取当前像素并作为前一像素缓存;
b1、重复步骤a1,直到处理完整个ROI窗口区域后,分别计算所的到的各个光斑灰度值和坐标值乘积的累加值与灰度累加值之商,并将所得到商作为对应光斑的质心的x轴坐标,以及计算所得到的各个光斑灰度值和y轴坐标值乘积的累加值与灰度累加值之商,并将所得到的商作为对应光斑的质心的y轴坐标。
CN200810222488A 2008-09-17 2008-09-17 一种高帧频太阳敏感器及其实现方法 Expired - Fee Related CN101363732B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200810222488A CN101363732B (zh) 2008-09-17 2008-09-17 一种高帧频太阳敏感器及其实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200810222488A CN101363732B (zh) 2008-09-17 2008-09-17 一种高帧频太阳敏感器及其实现方法

Publications (2)

Publication Number Publication Date
CN101363732A true CN101363732A (zh) 2009-02-11
CN101363732B CN101363732B (zh) 2010-05-12

Family

ID=40390232

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200810222488A Expired - Fee Related CN101363732B (zh) 2008-09-17 2008-09-17 一种高帧频太阳敏感器及其实现方法

Country Status (1)

Country Link
CN (1) CN101363732B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101487699B (zh) * 2009-03-06 2010-05-12 清华大学 一种高精度aps太阳敏感器及其实现方法
CN102034091A (zh) * 2010-12-01 2011-04-27 北京航空航天大学 一种用于数字式太阳敏感器的光斑识别方法及装置
CN102128608A (zh) * 2010-12-07 2011-07-20 北京航空航天大学 一种高动态二维姿态角测量方法和系统
CN102445287A (zh) * 2011-09-16 2012-05-09 浙江中控太阳能技术有限公司 一种表面能流密度测量系统和方法
CN103323025A (zh) * 2013-05-24 2013-09-25 哈尔滨工业大学 一种太阳光斑质心及其太阳输出角的误差修正方法
CN104034353A (zh) * 2014-06-06 2014-09-10 中国科学院长春光学精密机械与物理研究所 基于探测窗口的数字太阳敏感器质心计算方法
CN108284925A (zh) * 2018-01-09 2018-07-17 燕山大学 水线的识别检测方法
CN112581374A (zh) * 2019-09-29 2021-03-30 深圳市光鉴科技有限公司 散斑亚像素中心提取方法、系统、设备及介质
CN113155755A (zh) * 2021-03-31 2021-07-23 中国科学院长春光学精密机械与物理研究所 微透镜阵列型成像光谱仪在线标定方法
CN114543790A (zh) * 2021-12-30 2022-05-27 浙江众星志连科技有限责任公司 一种微小卫星太阳敏感器的鬼影去除方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1173155C (zh) * 2001-12-07 2004-10-27 清华大学 光学非线性补偿的模拟式两轴太阳敏感器
CN100405013C (zh) * 2007-06-05 2008-07-23 北京航空航天大学 太阳敏感器的复合光学系统及其实现方法
CN100527431C (zh) * 2007-08-30 2009-08-12 北京航空航天大学 太阳敏感器光电探测组件装配结构及安装方法

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101487699B (zh) * 2009-03-06 2010-05-12 清华大学 一种高精度aps太阳敏感器及其实现方法
CN102034091B (zh) * 2010-12-01 2013-05-22 北京航空航天大学 一种用于数字式太阳敏感器的光斑识别方法及装置
CN102034091A (zh) * 2010-12-01 2011-04-27 北京航空航天大学 一种用于数字式太阳敏感器的光斑识别方法及装置
CN102128608A (zh) * 2010-12-07 2011-07-20 北京航空航天大学 一种高动态二维姿态角测量方法和系统
CN102445287B (zh) * 2011-09-16 2014-07-09 浙江中控太阳能技术有限公司 一种表面能流密度测量系统和方法
CN102445287A (zh) * 2011-09-16 2012-05-09 浙江中控太阳能技术有限公司 一种表面能流密度测量系统和方法
CN103323025A (zh) * 2013-05-24 2013-09-25 哈尔滨工业大学 一种太阳光斑质心及其太阳输出角的误差修正方法
CN103323025B (zh) * 2013-05-24 2015-06-17 哈尔滨工业大学 一种太阳光斑质心及其太阳输出角的误差修正方法
CN104034353A (zh) * 2014-06-06 2014-09-10 中国科学院长春光学精密机械与物理研究所 基于探测窗口的数字太阳敏感器质心计算方法
CN108284925A (zh) * 2018-01-09 2018-07-17 燕山大学 水线的识别检测方法
CN112581374A (zh) * 2019-09-29 2021-03-30 深圳市光鉴科技有限公司 散斑亚像素中心提取方法、系统、设备及介质
CN113155755A (zh) * 2021-03-31 2021-07-23 中国科学院长春光学精密机械与物理研究所 微透镜阵列型成像光谱仪在线标定方法
CN113155755B (zh) * 2021-03-31 2022-05-24 中国科学院长春光学精密机械与物理研究所 微透镜阵列型成像光谱仪在线标定方法
CN114543790A (zh) * 2021-12-30 2022-05-27 浙江众星志连科技有限责任公司 一种微小卫星太阳敏感器的鬼影去除方法及系统
CN114543790B (zh) * 2021-12-30 2024-05-10 浙江众星志连科技有限责任公司 一种微小卫星太阳敏感器的鬼影去除方法及系统

Also Published As

Publication number Publication date
CN101363732B (zh) 2010-05-12

Similar Documents

Publication Publication Date Title
CN101363732B (zh) 一种高帧频太阳敏感器及其实现方法
CN103487058B (zh) 一种提高aps星敏感器动态性能的方法
CN101540892B (zh) 一种基于dsp视频采集设备的统计进出口人数的方法
CN106886216B (zh) 基于rgbd人脸检测的机器人自动跟踪方法和系统
CN101739692B (zh) 实时视频目标的快速相关跟踪方法
CN117250200B (zh) 一种基于机器视觉的方管生产质量检测系统
CN107238727A (zh) 基于动态视觉传感器芯片的光电式转速传感器及探测方法
CN109637128A (zh) 一种基于Markov的灰色Verhulst短时交通流预测方法及系统
CN105721772A (zh) 一种异步时域视觉信息成像方法
CN105354863A (zh) 基于特征滤波和快速运动检测模板预测的自适应尺度图像序列目标跟踪方法
CN105300390B (zh) 障碍物体运动轨迹的确定方法及装置
CN111105093B (zh) 基于CNN和BiLSTM的临近降水预测方法
Qimin et al. A methodology of vehicle speed estimation based on optical flow
CN104515473A (zh) 一种漆包线直径的在线检测方法
CN107301628B (zh) 基于震颤像移轨迹的卫星平台震颤影像去模糊方法
CN116310999B (zh) 一种水力发电站库区大型漂浮物检测方法
CN202404111U (zh) 基于dsp的视频测速系统
CN101509775B (zh) 一种多孔阵列式太阳敏感器的图像快速识别方法及装置
CN102034091A (zh) 一种用于数字式太阳敏感器的光斑识别方法及装置
CN115436371A (zh) 基于多源信息融合的Autoformer烧结矿FeO含量软测量方法
Ruan et al. On the use of sky images for intra-hour solar forecasting benchmarking: Comparison of indirect and direct approaches
CN105780631A (zh) 一种手推式水泥混凝土路面接缝损坏检测装置
CN108490849A (zh) 采用红外线传感器实现手势识别追踪的机器人及控制方法
Wei et al. Design of a sun tracking system based on the brightest point in sky image
CN116189229B (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100512

Termination date: 20190917