CN104061929B - 共光路多视场星敏感器及其星体姿态测量方法 - Google Patents
共光路多视场星敏感器及其星体姿态测量方法 Download PDFInfo
- Publication number
- CN104061929B CN104061929B CN201410321684.5A CN201410321684A CN104061929B CN 104061929 B CN104061929 B CN 104061929B CN 201410321684 A CN201410321684 A CN 201410321684A CN 104061929 B CN104061929 B CN 104061929B
- Authority
- CN
- China
- Prior art keywords
- star
- cornea
- punctate opacity
- high reflection
- visual fields
- 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.)
- Expired - Fee Related
Links
- 230000000007 visual effect Effects 0.000 title claims abstract description 61
- 238000000691 measurement method Methods 0.000 title claims abstract description 14
- 210000004087 cornea Anatomy 0.000 claims abstract description 68
- 230000003287 optical effect Effects 0.000 claims abstract description 27
- 238000003384 imaging method Methods 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 15
- 238000000605 extraction Methods 0.000 claims abstract description 12
- 238000012937 correction Methods 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 22
- 239000013598 vector Substances 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 10
- 238000013461 design Methods 0.000 abstract description 6
- 210000001747 pupil Anatomy 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Eye Examination Apparatus (AREA)
- Navigation (AREA)
Abstract
一种共光路多视场星敏感器及其星体姿态测量方法,该共光路多视场星敏感器共用一套光学***、成像单元和数据处理单元,光学***同时观测三个不同天区的恒星,成像单元对三个不同天区的恒星同时成像,对成像的星点进行星图识别,数据处理单元同时对采集到的恒星星点进行质心坐标的提取,并通过对三个视场的偏移量进行非对齐修正,最终完成高精度的星体姿态角输出。本发明可以有效的解决单视场星敏感器沿光轴方向测量精度较低的问题,同时也可克服传统多视场星敏感器采用分体式的光学***和成像***而引起的重量、功耗和体积较大的缺点,在实现三轴高精度测量的同时也实现了轻小型化设计,可以满足多种卫星平台对高精度姿态测量的需求。
Description
技术领域
本发明涉及一种共光路多视场星敏感器及其星体姿态测量方法。
背景技术
随着空间中对高精度姿态敏感器的需求日益加剧,为了获取准确的姿态信息,必须采用高精度的姿态测量敏感器。星敏感器是目前空间应用中精度最高的姿态测量设备,但是其三个测量方向中沿光学***光轴方向的测量精度较其它两轴相比要低一个数量级,仅靠单一视场的星敏感器难以做到三轴同时高精度测量。而常规的提高测量精度的方法是采用多台星敏感器数据相互融合的方式,或者采用多套光学***、多套成像单元组成的多视场星敏感器***,这样虽然可以实现三轴高精度测量,但同时也带来了诸如重量、功耗较大,安装位置难以精确标定等问题。
发明内容
本发明提供一种共光路多视场星敏感器及其星体姿态测量方法,可以有效的解决单视场星敏感器沿光轴方向测量精度较低的问题,同时也可克服传统多视场星敏感器采用分体式的光学***和成像***而引起的重量、功耗和体积较大的缺点,在实现三轴高精度测量的同时也实现了轻小型化设计,可以满足多种卫星平台对高精度姿态测量的需求。
为了达到上述目的,本发明提供一种共光路多视场星敏感器,该共光路多视场星敏感器包含光学***、成像单元和电路连接成像单元的数据处理单元;
所述的光学***包含:
第一高反射镜组,其设置在星敏感器的光线入口处,该第一高反射镜组包含三个第一高反射镜,所述的第一高反射镜为平面反射镜;
第二高反射镜组,其包含三个第二高反射镜,每个第二高反射镜与所述的每个第一高反射镜对应设置,每个第二高反射镜分别设置在星敏感器内能够接收并反射对应的第一高反射镜的反射光线的位置处,所述的第二高反射镜为凹面反射镜;
第三高反射镜,其设置在星敏感器内能够接收并反射第二高反射镜组内的第二高反射镜的反射光线的位置处,所述的第三高反射镜为凸面反射镜;
所述的成像单元接收第三高反射镜的反射光线,对不同视场方向进入的星点同时成像并完成星图识别;
所述的数据处理单元根据成像单元获得的星图,完成星点质心坐标的提取和星体姿态的解算。
所述的三个第一高反射镜按照两两之间呈120°夹角的方式,均匀布置,每个第一高反射镜与水平位置的倾角为45°。
本发明还提供一种利用共光路多视场星敏感器实现的基于共光路多视场星敏感器的星体姿态测量方法,该基于共光路多视场星敏感器的星体姿态测量方法包含以下步骤:
步骤1、光学***同时观测三个不同天区的恒星,并将三个不同天区的恒星同时成像到成像单元;
步骤2、成像单元对三个不同天区的恒星同时成像,对成像的星点进行星图识别;
步骤3、数据处理单元同时对采集到的恒星星点进行质心坐标的提取,并通过对三个视场的偏移量进行非对齐修正,最终完成高精度的星体姿态角输出。
所述的步骤2包含以下步骤:
步骤21、利用视场点扩散函数来完成成像星点所属视场的确定;
步骤22、采用基于四棱锥的方法,通过识别星的棱锥布局来完成星图识别。
所述的步骤21中,确定视场就是确定点扩散函数的对称轴的斜率,点扩散函数的对称轴的斜率与点扩散函数图像的惯性张量J的特征向量一致;
点扩散函数的对称轴的方位,通过惯性张量J的本征分析得到:
其中求和运算扩展到属于点扩散函数的所有像素;
特征值λ1和λ2按如下给出:
λa+λb=λ1>λ2=λa-λb (2)
其中:
结合特征值λ1和λ2,给出特征向量:
如果特征值λ1和λ2完全分开,就可以确定出星光来自哪个视场。
所述的步骤22包含以下步骤:
步骤2.1、对CCD/CMOS探测器成像面上的星点成像数量进行统计,如果星点数量为3,进行步骤2.2,如果星点数量小于3,则无法进行星图识别,匹配失败,如果星点数量大于3,进行步骤2.3;
步骤2.2、判断两两星点之间的星间角距,并与星图库进行对比,看结果是否唯一,如果唯一,星图匹配成功,完成星点识别并返回;如果不唯一,则匹配失败;
步骤2.3、从星点中取出三颗星点i、j、k,并取出另一颗参考星点r,分别判断r-i、r-j和r-k之间的星间角距,看结果是否唯一,如果唯一,进行步骤2.4,如果不唯一,进行步骤2.5;
步骤2.4、星图识别成功,可以使用已经识别出的i、j、k三颗恒星进行其它星点的识别;
步骤2.5、判断该i、j、k是否星图中轮寻过的最后三个点,如果不是,从星图中再取出三颗星,进行步骤2.3,如果是最后的三个星点,舍去负数解,看结果是否唯一,如果唯一,则星图识别成功;如果不唯一,则星图识别失败。
所述的步骤3中,数据处理单元对采集到的恒星星点进行质心坐标的提取包含以下步骤:
步骤31、快速获取质心的粗略位置;
步骤32、根据质心的粗略位置,精确确定质心位置;
所述的步骤31中,仅扫描所有的CCD/CMOS像素一次,通过整数向量Ih和Iv记录CCD/CMOS探测器上在相应行和相应列的最大值的位置,从而快速的找到所有n个最亮的像素点;
所述的向量Ih由nh个分量组成,向量Iv由nv个分量组成,相对应于CCD/CMOS探测器水平方向和垂直方向的像元数;
所述的步骤32中,精确的确定质心位置:
该方法利用了星点的灰度分布特征,认为像素点距离实际星点中心越近,其灰度值就越高,其他像素的灰度值随着距离中心增加而降低,依据像素的灰度值对其坐标加权,将加权后的几何中心作为实际成像的星点目标的中心;
设星点图像中第i行、第j列像素的坐标为(xi,yj),灰度值为f(xi,yj),则星点的质心位置计算式如下:
其中,(x0,y0)为目标像点质心坐标估计值。
所述的步骤3中,采用最优四元数估计方法计算星体姿态:
用四元数方法表示的姿态矩阵如下式所示:
Kqopt=λmaxqopt (5)
其中,qopt是姿态估计的最优四元数,它与对称阵K特征向量的最大特征值有关,K的表达式如下所示:
其中,向量z=Σiαisi×vi,可以用非对称阵B的非对角元素来表示,即
z={b23-b32,b31-b13,b12-b21}T (7)
下面给出K矩阵的特征多项式:
λ4+aλ3+bλ2+cλ+d=0 (8)
其中,a=tr[K]=0,b、c、d可用下式表示:
在此引入(8)的辅助方程
u3-bu2-4du+4bd-c2=0 (10)
辅助方程的解为
其中,
当得到u1后,(8)式的特征值也就得到了:
其中,
很容易看出-1≤λ1≤λ2≤λ3≤λ4=λmax≤1
如果K矩阵的特征值数n=2,那么
其中,
当n>2时,姿态四元数qopt可以通过(13)式和(14)式计算得到;当n=2时,可以通过(15)式和(16)式得到。
本发明提出的共光路多视场星敏感器通过采取共用一套光学***、一套成像***和一套电子学单元的方式,可以有效的解决单视场星敏感器沿光轴方向测量精度较低的问题,同时也可克服传统多视场星敏感器采用分体式的光学***和成像***而引起的重量、功耗和体积较大的缺点,在实现三轴高精度测量的同时也实现了轻小型化设计,可以满足多种卫星平台对高精度姿态测量的需求。
附图说明
图1是本发明的结构示意图。
图2是本发明的星点成像点扩散函数PSF图。
图3是本发明的星图识别的流程图。
具体实施方式
以下根据图1~图3,具体说明本发明的较佳实施例。
如图1所示,本发明提供一种共光路多视场星敏感器,其包含光学***1、成像单元2和连接成像单元的数据处理单元3。
所述的光学***包含:
第一高反射镜组,其设置在星敏感器的光线入口(入瞳位置)处,该第一高反射镜组包含三个第一高反射镜101,所述的三个第一高反射镜101按照两两之间呈120°夹角的方式,均匀布置,每个第一高反射镜101与水平位置的倾角为45°,所述的第一高反射镜101为平面反射镜;
第二高反射镜组,其包含三个第二高反射镜102,每个第二高反射镜102与所述的每个第一高反射镜101对应设置,每个第二高反射镜102分别设置在星敏感器内能够接收并反射对应的第一高反射镜101的反射光线的位置处,所述的第二高反射镜102为凹面反射镜;
第三高反射镜103,其设置在星敏感器内能够接收并反射第二高反射镜组内的第二高反射镜102的反射光线的位置处,所述的第三高反射镜103为凸面反射镜。
所述的成像单元2接收第三高反射镜103的反射光线,对不同视场方向进入的星点同时成像并完成星图识别;该成像单元2采用CCD探测器(电荷耦合器件,Charge CoupledDevice)或者CMOS探测器(互补金属氧化物半导体,Complementary Metal-Oxide-Semiconductor Transistor)。
所述的数据处理单元3根据成像单元2获得的星图,完成星点质心坐标的提取和星体姿态的解算。
本发明还提供一种基于共光路多视场星敏感器的星体姿态测量方法,该方法包含以下步骤:
步骤1、光学***同时观测三个不同天区的恒星,并将三个不同天区的恒星同时成像到成像单元;
步骤2、成像单元对三个不同天区的恒星同时成像,对成像的星点进行星图识别;
步骤3、数据处理单元同时对采集到的恒星星点进行质心坐标的提取,并通过对三个视场的偏移量进行非对齐修正,最终完成高精度的星体姿态角输出;
所述的步骤1中,如图1所示,光学***采用施密特-卡塞格林折反式光学设计方法,可以有效的实现将多个视场的星点成像于同一个成像探测器上,也易于分辨图像中的星点来自于哪个视场。在星敏感器光线入口处按照两两之间120°夹角的方式,均匀布置三个45°倾角的第一高反射镜101,与安装在星敏感器内部的第二高反射镜102和高反射镜103共同作用,通过三次折反,将天球上的星点成像在同一片CCD/CMOS探测器上。三视场星敏感器的焦距与入瞳孔径的比值F大约为3,每个视场的视场角大约为20°。该光学***必须满足的条件是在物镜(第一高反射镜组和第二高反射镜组)前端的入瞳位置和形状要与每个视场相匹配,按照像方远心光路的设计方式,光阑104(决定通光孔径的大小)需要安装在物方焦平面位置处,出瞳应该在无穷远。
所述的步骤2中,可以利用典型的PSF瑞利等高线分布图来完成视场中的星图识别,由于入瞳分割,每个视场的PSF(点扩散函数,Point Spread Function)都有相同的关于轴线对称的形状,对称轴之间两两相差120°。
如图2所示,CCD/CMOS成像单元可以对不同视场方向进入的星点同时成像,根据PSF的形状,完成成像星点所属视场的确定,并采用基于四棱锥的星点组合方法,完成星图识别。施密特-卡塞格林的这种在入瞳处利用几何分割的解决方法会产生PSF,它们的形状与视场中所出现的恒星相关。因此,基于这个重要原则,通过PSF的形状来确定星光来自于哪个视场成为可能,因为PSF是轴对称,并把视场分割成两个对称的部分。本发明中每个视场的PSF形状都关于轴线对称,并且对称轴之间两两相差120°。因此,视场确定的问题得到了简化,转化为确定对称轴斜率的问题,这个斜率与图像的惯性张量J的特征向量一致,并且,此对称轴与最大特征值有关。这种解决方法特别适合入瞳处几何分割技术。因而,通过改变一些光学参数,可以获得不同形状的PSF,可以通过星的亮度来确定对称轴所对应的视场。目前,最为有效的方法就是得到环面形状的PSF。
进一步,对称轴(惯性轴)的方位,可以通过图像的惯性张量J的本征分析很容易得到。
其中求和运算扩展到属于PSF的所有像素;
特征值λ1和λ2按如下给出:
λa+λb=λ1>λ2=λa-λb (2)
其中:
结合特征值λ1和λ2,给出特征向量:
由此可以看出,对称轴(特征向量)是和特征值相关,可以明确地确定星光来源的视场。用这种方法,如果特征值λ1和λ2完全分开,就可以确定出星光来自哪个视场。
进一步,所述步骤2中星图识别的过程可以采用基于四棱锥的方法,通过识别星的棱锥布局来完成星图识别,因为在星图中四颗星的结构的误匹配的频率非常小(10-12),所以星识别几乎是很可靠的。如图3所示(其中的整形向量I、J和K代表星点坐标),基于四棱锥的星图识别的步骤包含:
步骤2.1、对CCD/CMOS探测器成像面上的星点成像数量进行统计,如果星点数量为3,进行步骤2.2,如果星点数量小于3,则无法进行星图识别,匹配失败,如果星点数量大于3,进行步骤2.3;
步骤2.2、判断两两星点之间的星间角距,并与星图库进行对比,看结果是否唯一,如果唯一,星图匹配成功,完成星点识别并返回;如果不唯一,则匹配失败;
步骤2.3、从星点中取出三颗星点i、j、k,并取出另一颗参考星点r,分别判断r-i、r-j和r-k之间的星间角距,看结果是否唯一,如果唯一,进行步骤2.4,如果不唯一,进行步骤2.5;
步骤2.4、星图识别成功,可以使用已经识别出的i、j、k三颗恒星进行其它星点的识别;
步骤2.5、判断该i、j、k是否星图中轮寻过的最后三个点,如果不是,从星图中再取出三颗星,进行步骤2.3,如果是最后的三个星点,舍去负数解,看结果是否唯一,如果唯一,则星图识别成功;如果不唯一,则星图识别失败。
进一步,所述的步骤3中,所述的数据处理单元完成星点质心坐标的提取和星体姿态的解算。
为了提高星敏感器的数据更新率,在星点质心提取方面需要快速计算,本发明中采用仅扫描所有的CCD/CMOS像素一次,即可找到所有需要的星斑在成像探测器上的位置信息。通过引进两个整数向量,Ih和Iv(向量Ih由nh个分量组成,向量Iv由nv个分量组成)相对应于CCD/CMOS探测器水平方向和垂直方向的像元数。这两个向量Ih和Iv记录了CCD/CMOS探测器上在相应行和相应列的最大值的位置。具体来说,倘若位于第6行的最亮的像元坐落在第143列,第67列最亮的像元坐落在第1行,那么Iv(6)=143,并且Ih(67)=1,这两个整数向量可以快速的找到所有n个最亮的像素点(很可能这些像素点就是恒星在成像探测器上成像的斑点),避免后续的CCD/CMOS扫描,及避免了时间的消耗。
进一步的,在快速获取质心的粗略位置后,还需要精确的确定质心位置,该方法利用了星点的灰度分布特征,认为像素点距离实际星点中心越近,其灰度值就越高,其他像素的灰度值随着距离中心增加而降低,依据像素的灰度值对其坐标加权,将加权后的几何中心作为实际成像的星点目标的中心;
设星点图像中第i行、第j列像素的坐标为(xi,yj),灰度值为f(xi,yj),则星点的质心位置计算式如下:
其中,(x0,y0)为目标像点质心坐标估计值。
对于质心计算的准确度的问题是如何准确的确定星光入射方向,在使用掩模轮廓处理时,需要依赖于CCD/CMOS点扩散函数的形状。因为复杂的点扩散函数的形状并不能很好的符合正方形和矩形掩模,因此,需要采用新的方法来处理这个问题,通过指定的椭圆形的掩膜来提高质心计算的准确度。相比较于标准的正方形掩膜,椭圆形的掩膜的优点是减小了远离点扩散函数中心(质心)的像元由于电子学噪声产生的坏点的对处理的影响。第二个重要的提高是有关于使用相关的递归函数(如lm迭代、牛顿迭代等等),采用迭代运算的方法提高质心提取精度。这些函数,在不同的约束下不同,允许识别任何复杂的PSF形状。
进一步的,在完成星点质心坐标的精确定位后,需要根据质心信息完成最终的姿态解算。对于多视场星敏感器而言,其探测器上成像的星点数量要远远多于普通的单视场星敏感器,如何快速完成姿态的解算成为一项最重要的问题,本发明使用了ESOQ(最优四元数估计)基本算法,通过使用一种更快速的编码方式,完成线性快速运算。用四元数方法表示的姿态矩阵如下式所示:
Kqopt=λmaxqopt (5)
其中,qopt是姿态估计的最优四元数,它与对称阵K特征向量的最大特征值有关,K的表达式如下所示:
其中,向量z=Σiαisi×vi,可以用非对称阵B的非对角元素来表示,即
z={b23-b32,b31-b13,b12-b21}T (7)
下面给出K矩阵的特征多项式:
λ4+aλ3+bλ2+cλ+d=0 (8)
其中,a=tr[K]=0,b、c、d可用下式表示:
在此引入(8)的辅助方程
u3-bu2-4du+4bd-c2=0 (10)
辅助方程的解为
其中,
当得到u1后,(8)式的特征值也就得到了:
其中,
很容易看出-1≤λ1≤λ2≤λ3≤λ4=λmax≤1
如果K矩阵的特征值数n=2,那么
其中,
当n>2时,姿态四元数qopt可以通过(13)式和(14)式计算得到;当n=2时,可以通过(15)式和(16)式得到。
本发明与现有技术相比,其优点和有益效果是:整体结构紧凑,集成度高,通过采用一套光学***、一套成像***和一套电子学单元,有效的解决了单视场星敏感器沿光轴方向测量精度较低的问题,同时也可克服传统多视场星敏感器采用分体式的光学***和成像***而引起的重量、功耗和体积较大的缺点,在实现三轴高精度测量的同时也实现了轻小型化设计。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (8)
1.一种共光路多视场星敏感器,其特征在于,该共光路多视场星敏感器包含光学***(1)、成像单元(2)和电路连接成像单元的数据处理单元(3);
所述的光学***包含:
第一高反射镜组,其设置在星敏感器的光线入口处,该第一高反射镜组包含三个第一高反射镜(101),所述的第一高反射镜(101)为平面反射镜;
第二高反射镜组,其包含三个第二高反射镜(102),每个第二高反射镜(102)与所述的每个第一高反射镜(101)对应设置,每个第二高反射镜(102)分别设置在星敏感器内能够接收并反射对应的第一高反射镜(101)的反射光线的位置处,所述的第二高反射镜(102)为凹面反射镜;
第三高反射镜(103),其设置在星敏感器内能够接收并反射第二高反射镜组内的第二高反射镜(102)的反射光线的位置处,所述的第三高反射镜(103)为凸面反射镜;
所述的成像单元(2)接收第三高反射镜(103)的反射光线,对不同视场方向进入的星点同时成像并完成星图识别;
所述的数据处理单元(3)根据成像单元(2)获得的星图,完成星点质心坐标的提取和星体姿态的解算。
2.如权利要求1所述的共光路多视场星敏感器,其特征在于,所述的三个第一高反射镜(101)按照两两之间呈120°夹角的方式,均匀布置,每个第一高反射镜(101)与水平位置的倾角为45°。
3.一种利用如权利要求1或2所述的共光路多视场星敏感器实现的基于共光路多视场星敏感器的星体姿态测量方法,其特征在于,该基于共光路多视场星敏感器的星体姿态测量方法包含以下步骤:
步骤1、光学***同时观测三个不同天区的恒星,并将三个不同天区的恒星同时成像到成像单元;
步骤2、成像单元对三个不同天区的恒星同时成像,对成像的星点进行星图识别;
步骤3、数据处理单元同时对采集到的恒星星点进行质心坐标的提取,并通过对三个视场的偏移量进行非对齐修正,最终完成高精度的星体姿态角输出。
4.如权利要求3所述的基于共光路多视场星敏感器的星体姿态测量方法,其特征在于,所述的步骤2包含以下步骤:
步骤21、利用视场点扩散函数来完成成像星点所属视场的确定;
步骤22、采用基于四棱锥的方法,通过识别星的棱锥布局来完成星图识别。
5.如权利要求4所述的基于共光路多视场星敏感器的星体姿态测量方法,其特征在于,所述的步骤21中,确定视场就是确定点扩散函数的对称轴的斜率,点扩散函数的对称轴的斜率与点扩散函数图像的惯性张量J的特征向量一致;
点扩散函数的对称轴的方位,通过惯性张量J的本征分析得到:
其中求和运算扩展到属于点扩散函数的所有像素;
特征值λ1和λ2按如下给出:
λa+λb=λ1>λ2=λa-λb (2)
其中:
结合特征值λ1和λ2,给出特征向量:
如果特征值λ1和λ2完全分开,就可以确定出星光来自哪个视场。
6.如权利要求4所述的基于共光路多视场星敏感器的星体姿态测量方法,其特征在于,所述的步骤22包含以下步骤:
步骤2.1、对CCD/CMOS探测器成像面上的星点成像数量进行统计,如果星点数量为3,进行步骤2.2,如果星点数量小于3,则无法进行星图识别,匹配失败,如果星点数量大于3,进行步骤2.3;
步骤2.2、判断两两星点之间的星间角距,并与星图库进行对比,看结果是否唯一,如果唯一,星图匹配成功,完成星点识别并返回;如果不唯一,则匹配失败;
步骤2.3、从星点中取出三颗星点i、j、k,并取出另一颗参考星点r,分别判断r-i、r-j和r-k之间的星间角距,看结果是否唯一,如果唯一,进行步骤2.4,如果不唯一,进行步骤2.5;
步骤2.4、星图识别成功,可以使用已经识别出的i、j、k三颗恒星进行其它星点的识别;
步骤2.5、判断该i、j、k是否星图中轮寻过的最后三个点,如果不是,从星图中再取出三颗星,进行步骤2.3,如果是最后的三个星点,舍去负数解,看结果是否唯一,如果唯一,则星图识别成功;如果不唯一,则星图识别失败。
7.如权利要求4所述的基于共光路多视场星敏感器的星体姿态测量方法,其特征在于,所述的步骤3中,数据处理单元对采集到的恒星星点进行质心坐标的提取包含以下步骤:
步骤31、快速获取质心的粗略位置;
步骤32、根据质心的粗略位置,精确确定质心位置;
所述的步骤31中,仅扫描所有的CCD/CMOS像素一次,通过整数向量Ih和Iv记录CCD/CMOS探测器上在相应行和相应列的最大值的位置,从而快速的找到所有n个最亮的像素点;
所述的向量Ih由nh个分量组成,向量Iv由nv个分量组成,相对应于CCD/CMOS探测器水平方向和垂直方向的像元数;
所述的步骤32中,精确的确定质心位置;
设星点图像中第i行、第j列像素的坐标为(xi,yj),灰度值为f(xi,yj),则星点的质心位置计算式如下:
其中,(x0,y0)为目标像点质心坐标估计值。
8.如权利要求7所述的基于共光路多视场星敏感器的星体姿态测量方法,其特征在于,所述的步骤3中,采用最优四元数估计方法计算星体姿态:
用四元数方法表示的姿态矩阵如下式所示:
Kqopt=λmaxqopt (5)
其中,qopt是姿态估计的最优四元数,它与对称阵K特征向量的最大特征值有关,K的表达式如下所示:
其中,向量z=Σiαisi×vi,可以用非对称阵B的非对角元素来表示,即
z={b23-b32,b31-b13,b12-b21}T (7)
下面给出K矩阵的特征多项式:
λ4+aλ3+bλ2+cλ+d=0 (8)
其中,a=tr[K]=0,b、c、d可用下式表示:
在此引入(8)的辅助方程
u3-bu2-4du+4bd-c2=0 (10)
辅助方程的解为
其中,
当得到u1后,(8)式的特征值也就得到了:
其中,
很容易看出-1≤λ1≤λ2≤λ3≤λ4=λmax≤1
如果K矩阵的特征值数n=2,那么
其中,
当n>2时,姿态四元数qopt可以通过(13)式和(14)式计算得到;当n=2时,可以通过(15)式和(16)式得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410321684.5A CN104061929B (zh) | 2014-07-08 | 2014-07-08 | 共光路多视场星敏感器及其星体姿态测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410321684.5A CN104061929B (zh) | 2014-07-08 | 2014-07-08 | 共光路多视场星敏感器及其星体姿态测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104061929A CN104061929A (zh) | 2014-09-24 |
CN104061929B true CN104061929B (zh) | 2017-03-08 |
Family
ID=51549756
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410321684.5A Expired - Fee Related CN104061929B (zh) | 2014-07-08 | 2014-07-08 | 共光路多视场星敏感器及其星体姿态测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104061929B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104501805B (zh) * | 2014-12-25 | 2017-07-14 | 中国科学院长春光学精密机械与物理研究所 | 大相对孔径高精度折反射型星敏感器光学*** |
CN105806346B (zh) * | 2014-12-31 | 2018-09-07 | 上海新跃仪表厂 | 中高轨星座星间照相观测敏感器及其星间角距测量方法 |
CN105425391B (zh) * | 2015-11-26 | 2019-08-02 | 中国航空工业第六一八研究所 | 一种机载星敏感器光学*** |
CN105547286B (zh) * | 2016-01-11 | 2018-04-10 | 常州工学院 | 一种复合三视场星敏感器星图模拟方法 |
CN107677264B (zh) * | 2017-08-15 | 2020-09-18 | 北京控制工程研究所 | 一种反射式星敏感器 |
CN108681617B (zh) * | 2018-03-29 | 2022-07-29 | 北京空间飞行器总体设计部 | 一种航天器多星敏感器布局优化设计方法 |
CN109520526B (zh) * | 2019-01-24 | 2023-04-18 | 中科院南京天文仪器有限公司 | 一种基于共光路的星模拟器标定与自准直测量***及方法 |
CN110503086B (zh) * | 2019-08-14 | 2021-10-22 | 中国科学院近代物理研究所 | 一种基于奇偶位置的目标质心定位***及方法 |
CN112212856B (zh) * | 2020-09-21 | 2022-08-19 | 中国科学院长春光学精密机械与物理研究所 | 一种近地面白天星敏感器的高精度姿态测量方法 |
CN112653828B (zh) * | 2020-12-28 | 2022-11-04 | 维沃移动通信有限公司 | 摄像模组、电子设备及其控制方法和控制装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4192015A (en) * | 1978-03-06 | 1980-03-04 | Bell Telephone Laboratories, Incorporated | Optical image sensor semiconductor apparatus |
CN101592490A (zh) * | 2009-07-06 | 2009-12-02 | 北京航空航天大学 | 一种自适应控阵星敏感器 |
CN101770072A (zh) * | 2009-12-30 | 2010-07-07 | 北京控制工程研究所 | 一种复合视场敏感器成像*** |
CN102116926A (zh) * | 2009-12-31 | 2011-07-06 | 北京控制工程研究所 | 一种恒星敏感器成像结构 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7372642B2 (en) * | 2006-02-13 | 2008-05-13 | 3M Innovative Properties Company | Three-channel camera systems with non-collinear apertures |
-
2014
- 2014-07-08 CN CN201410321684.5A patent/CN104061929B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4192015A (en) * | 1978-03-06 | 1980-03-04 | Bell Telephone Laboratories, Incorporated | Optical image sensor semiconductor apparatus |
CN101592490A (zh) * | 2009-07-06 | 2009-12-02 | 北京航空航天大学 | 一种自适应控阵星敏感器 |
CN101770072A (zh) * | 2009-12-30 | 2010-07-07 | 北京控制工程研究所 | 一种复合视场敏感器成像*** |
CN102116926A (zh) * | 2009-12-31 | 2011-07-06 | 北京控制工程研究所 | 一种恒星敏感器成像结构 |
Non-Patent Citations (2)
Title |
---|
多视场星敏感器结构布局优化;王真等;《红外与激光工程》;20111231;第40卷(第12期);2469-2473 * |
新型反射式星敏感器光学***设计;赵阳;《中国优秀硕士学位论文全文数据库工程科技辑》;20090215;C031-112 * |
Also Published As
Publication number | Publication date |
---|---|
CN104061929A (zh) | 2014-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104061929B (zh) | 共光路多视场星敏感器及其星体姿态测量方法 | |
CN104730802B (zh) | 光轴夹角的校准、对焦方法和***和双摄像头设备 | |
US6985183B2 (en) | Method for exploring viewpoint and focal length of camera | |
CN110517325B (zh) | 一种坐标变换及坐标变换的车身周围物体定位方法和*** | |
CN111537003B (zh) | 一种基于折射面共线的星光大气折射测量校正方法 | |
CN111709985B (zh) | 一种基于双目视觉的水下目标测距方法 | |
CN109540113B (zh) | 一种全站仪及其星图识别方法 | |
CN103426149A (zh) | 大视角图像畸变的校正处理方法 | |
CN105806346B (zh) | 中高轨星座星间照相观测敏感器及其星间角距测量方法 | |
US9453731B2 (en) | System and method for determining orientation relative to earth | |
CN110310338A (zh) | 一种基于多中心投影模型的光场相机标定方法 | |
CN106595702B (zh) | 一种基于天文标定的多传感器空间配准方法 | |
CN102072706A (zh) | 一种多相机定位与跟踪方法及*** | |
CN110009682A (zh) | 一种基于单目视觉的目标识别定位方法 | |
CN113324528B (zh) | 已知摄站位置的近景摄影测量目标定位方法及*** | |
CN111238540A (zh) | 一种基于恒星拍摄的珞珈一号相机-星敏安装标定方法 | |
CN104977000A (zh) | 中高轨星座星间照相观测敏感器及其星间角距测量算法 | |
CN103925922B (zh) | 一种适用于高动态条件下iccd星图的恒星识别方法 | |
CN105547286B (zh) | 一种复合三视场星敏感器星图模拟方法 | |
CN103134443B (zh) | 一种大口径大径厚比反射镜面形自准直检测装置及方法 | |
CN111998823B (zh) | 基于双目异光源测距装置的目标测距方法 | |
ES2933452T3 (es) | Procedimiento y dispositivo para la medición de la potencia de refracción local y/o de la distribución de la potencia de refracción de una lente de gafas | |
JPH02295535A (ja) | 可動部を要しない瞳孔間隔測定装置 | |
CN101285735B (zh) | 一种通过分离波面整体倾斜扩大动态范围的哈特曼传感器 | |
CN111024063B (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 |
Granted publication date: 20170308 |