CN1948085A - 一种基于星场的星敏感器校准方法 - Google Patents
一种基于星场的星敏感器校准方法 Download PDFInfo
- Publication number
- CN1948085A CN1948085A CN 200510112552 CN200510112552A CN1948085A CN 1948085 A CN1948085 A CN 1948085A CN 200510112552 CN200510112552 CN 200510112552 CN 200510112552 A CN200510112552 A CN 200510112552A CN 1948085 A CN1948085 A CN 1948085A
- Authority
- CN
- China
- Prior art keywords
- axle
- overbar
- centerdot
- cos
- star
- 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
Links
Images
Landscapes
- Studio Devices (AREA)
Abstract
本发明属于航天测量技术,涉及对星敏感器校准方法的改进。本发明的步骤是:建立星敏感器姿态转换矩阵;建立星敏感器畸变模型;进行最小二乘参数估计。本发明方法同时将姿态估计过程和星敏感器参数校准过程结合起来,从而消除了姿态误差引入对星敏感器参数的影响;适用广泛,不仅可以应用于在轨校准,也可以应用于地面的夜空观测校准和实验室内部的星场校准;有着很好数值稳定性。
Description
技术领域
本发明属于航天测量技术,涉及对星敏感器校准方法的改进。
背景技术
星敏感器是一种利用恒星观测,为空间飞行器提供高精度姿态信息的航天测量仪器。其工作原理为:星敏感器前端摄像头单元利用CCD(或CMOS)图像传感器拍摄得到星图图像,经过图像处理程序得到恒星像点的质心坐标和亮度的信息,然后星图识别程序利用这些信息在导航星库中找到对应的恒星,最后计算出星敏感器的三轴姿态。
星敏感器校准一般分为地面校准和在轨校准两种方式。在飞行器发射前,星敏感器要首先进行地面校准和测试,在实验室的条件下标定出星敏感器的内部参数,如主点、焦距和畸变系数等。在轨校准是在飞行器发射升空后,星敏感器根据发射时的冲击和工作环境的变化情况,如重力、大气和温度等都会不同于地面情况,需要在地面校准参数的基础上适时地对这些参数进行修正。
如图1所示,П表示星敏感器靶面,星敏感器主点为O,图像坐标为(x0,y0),焦距为f,表示镜头中心到靶面的距离。以主点为原点建立星敏感器靶面坐标系,Z轴为主光轴,X轴和Y轴分别对应星敏感器图像采集的行方向和列方向。设星点i和j的成像点为Pi(xi,yi),Pj(xj,yj)。星库星向量夹角成为星内角,设为θij,同时测量成像点向量夹角,称为测量星内角,设为θ′ij。根据星向量夹角的正交变换不变性,理想小孔成像模型时,θij=θ′ij。
主点和焦距的校准。最常见的一种主点和焦距的在轨校准方法是利用星库星向量间夹角和测量星向量间夹角相等原理(如果测量参数出现偏差,则该夹角也将会出现偏差),采用非线性参数优化的方法,求解出主点和焦距偏差,使得最终的星向量夹角误差在最小二乘意义下最小。
畸变系数的校准。在轨校准畸变系数的方法常见的有两种,一种是利用姿态信息的校准方法,另一种则是不依赖姿态信息的校准方法。利用姿态信息的校准方法是:跟据主点和焦距建立理想的小孔成像模型,然后利用姿态信息计算星场星向量在靶面上的理想投影位置,进而计算它和实际测量的星点质心位置偏差。利用2维多项式模型来拟合这一偏差,完成对畸变的校准。这种方法的一个明显问题是由于对姿态信息的依赖而引入误差。由于星敏感器本身是空间飞行器上精度最高的姿态测量仪器,飞行器在星敏感器校准前只能利用陀螺或者其它姿态测量仪器获得当前姿态信息,这些外部姿态信息带有明显的误差。于是该校准方法将会把这些姿态误差引入星敏感器的参数估计过程中,引起星敏感器的最终工作误差。
不依赖姿态信息的校准方法是:和校准主点和焦距的方法相似,利用星内角不变原理,计算测量星内角和星库星内角的偏差。利用2维多项式建立测量星向量模型,由于不依赖姿态信息,需要将两个测量星向量相减来表达星向量之间的相对值,然后利用2维多项式模型来拟合该相对值,求出最小二乘解。这种方法的第一个缺点是由于星向量相减而损失了0阶项,使得误差拟合不完全;第二个缺点是算法对2维多项式的建模方法有很强的依赖性,模型不好会导致数值计算不稳定,误差会很大。
发明内容
本发明的目的是:针对上述星敏感器在轨校准方法存在的问题,提出一种基于星场的星敏感器校准方法。该方法对星敏感器姿态参数和内部参数同时进行估计运算,并利用非线性最小二乘法和共线性公式,迭代计算出星敏感器姿态参数和内部参数。该方法不仅可以用在星敏感器的在轨校准过程中,也可以用在夜空拍摄和实验室校准,因此该方法有着广泛是适用性。
本发明的技术方案是:一种基于星场的星敏感器校准方法,其特征在于,
1、建立星敏感器姿态转换矩阵;
1.1、建立天球坐标系;以地球中心O′为坐标原点,过春分点的轴线为Xn轴,过北极点的轴线为Zn轴,Yn轴则定义为垂直于XnZn平面的直线;
1.2、建立星敏感器坐标系O-XYZ,以主点为原点建立星敏感器靶面坐标系,Z轴为主光轴,X轴和Y轴分别对应星敏感器图像采集的行方向和列方向;星敏感器的姿态角由赤经α0、赤纬β0、和滚转角φ0组成,α0为Z轴在XnYn面上的投影与Xn轴的夹角,从Xn轴起逆时针计算;β0为Z轴与它在XnYn面上的投影之间的夹角,从投影逆时针计算;φ0为Zn轴在XY面上的投影与Yn轴的夹角,从投影起顺时针计算;φ0为天球坐标系子午面和像平面的交线和像面Y轴之间的夹角;
1.3、建立姿态转换矩阵;天球坐标系经过三次旋转即可到达星敏感器坐标系,旋转过程为:第一次绕Zn轴旋转
使得Xn轴和子午面相垂直;第二次绕旋转后的Xn轴旋转
使得Zn轴与Z轴重合;第三次绕两次旋转后的Zn轴旋转φ0,则天球坐标系O′-XnYnZn和星敏感器坐标系O-XYZ重合;设姿态转换矩阵为M,则有:
这里:
a1=sinα0cos0-cosα0sinβ0sin0;
a2=-sinα0sin0-cosα0sinβ0cos0;
a3=-cosα0cos0;
b1=-cosα0cos0-sinα0sinβ0sin0;
b2=cosα0sin0-sinα0sinβ0cos0;
b3=-sinα0cos0;
c1=cosβ0sin0;
c2=cosβ0cos0;
c3=-sinβ0;
设星场中有n颗星,第i颗星的惯性坐标为(αi,βi),i=1,...,n,则第i颗星的方向矢量为:
这里n1,n2,n3为星光向量在天球坐标系3个轴的投影;
假设焦距为fc,星敏感器成像的共线性公式为:
这里
(x0,y0)为主点坐标,(x,y)为像平面坐标;
2、建立星敏感器畸变模型;
假设dx,dy为x方向和y方向的畸变偏差,有:
这里,
r2=
x2+
y2;
q1,q2,q3为径向畸变系数;
p1,p2,p3为偏心畸变系数;
这里内部参数有9个,即为x0,y0,fc,q1,q2,q3,p1,p2,p3,整体成像的共线性公式为:
3、最小二乘参数估计;
主点偏差采用地面自准直校准结果;在星场校准过程中,假设主点位置为已知,这样总共需要考虑的校准参数有10个,即为α0,β0,φ0,fc,q1,q2,q3,p1,P2,p3,用参数向量来表示为:
根据共线性公式有:
假设星点采集个数为m,联合x和y方向的偏差和敏感矩阵,假设:
这里p是由x和y方向残留偏差组成的向量,M为A和B两个敏感矩阵组成的整体敏感矩阵;
于是有迭代方程为:
这里k为迭代序号,k取5~20,迭代结束后得到稳定的数据值,这时的参数即为最后的校准结果。
本发明的优点是:第一,该基于星场的方法同时将姿态估计过程和星敏感器参数校准过程结合起来,从而消除了姿态误差引入对星敏感器参数的影响;第二,该方法适用广泛,不仅可以应用于在轨校准,也可以应用于地面的夜空观测校准和实验室内部的星场校准。第三,该方法有着很好数值稳定性。
附图说明
图1是星敏感器星场成像示意图。
图2是本发明方法中星敏感器姿态角示意图。
图3是模拟星场示意图。
具体实施方式
下面对本发明做进一步详细说明。星场校准的方法在星敏感器的在轨校准、夜空拍摄校准和实验室利用星场模拟器进行校准的场合都是需要的,因此研究一种通用的基于星场的校准方法有着现实意义。本发明方法星敏感器姿态参数和内部参数同时进行估计运算,并利用非线性最小二乘法和共线性公式,迭代计算出星敏感器姿态参数和内部参数。它不仅可以用在星敏感器的在轨校准过程中,也可以用在夜空拍摄和实验室校准,因此本发明方法有着广泛是适用性。同时,由于该方法同时将星敏感器的姿态参数和内部参数同时进行校准,有效消除了姿态误差引入带来的星敏感器内部参数估计误差。本发明方法的步骤如下:
1、建立星敏感器姿态转换矩阵。
1.1、建立天球坐标系。在球面天文学中,为了与人们的直观感觉相适应,把天空假想成一个巨大的球面,即为天球。天球坐标系以地球中心O′为坐标原点,过春分点的轴线为Xn轴,过北极点的轴线为Zn轴,Yn轴则定义为垂直于XnZn平面的直线,如图2所示。
1.2、建立星敏感器坐标系。图2中O-XYZ为星敏感器坐标系。星敏感器的姿态角由赤经α0、赤纬β0、和滚转角φ0组成。这里α0为Z轴在XnYn面上的投影与Xn轴的夹角,从Xn轴起逆时针计算;β0为Z轴与它在XnYn面上的投影之间的夹角,从投影逆时针计算;φ0为Zn轴在XY面上的投影与Yn轴的夹角,从投影起顺时针计算。
1.3、建立姿态转换矩阵。φ0为天球坐标系子午面和像平面的交线和像面Y轴之间的夹角。天球坐标系经过三次旋转即可到达星敏感器坐标系。旋转过程为:第一次绕Zn轴旋转
使得Xn轴和子午面相垂直;第二次绕旋转后的Xn轴旋转
使得Zn轴与Z轴重合;第三次绕两次旋转后的Zn轴旋转φ0,则天球坐标系O′-XnYnZn和星敏感器坐标系O-XYZ重合。设姿态转换矩阵为M,则有:
[1]
这里:
a1=sinα0cos0-cosα0sinβ0sin0;
a2=-sinα0sin0-cosα0sinβ0cos0;
a3=-cosα0cos0;
b1=-cosα0cos0-sinα0sinβ0sin0;
b2=cosα0sin0-sinα0sinβ0cos0;
b3=-sinα0cos0;
c1=cosβ0sin0;
c2=cosβ0cos0;
c3=-sinβ0。
设星场中有n颗星,第i颗星的惯性坐标为(αi,βi),i=1,...,n,则第i颗星的方向矢量为:
这里n1,n2,n3为星光向量在天球坐标系3个轴的投影。
假设焦距为fc,星敏感器成像的共线性公式为:
这里
(x0,y0)为主点坐标,(x,y)为像平面坐标。
2、建立星敏感器畸变模型。
假设dx,dy为x方向和y方向的畸变偏差,有:
这里,
r2=
x2+
y2;
q1,q2,q3为径向畸变系数;
p1,p2,p3为偏心畸变系数。
这里内部参数有9个,即为(x0,y0,fc,q1,q2,q3,p1,p2,p3),整体成像的共线性公式为:
3、最小二乘参数估计。
主点偏差对于星敏感器的工作精度影响不大,可以采用地面自准直校准结果。在该星场校准过程中,假设主点位置为已知,这样总共需要考虑的校准参数有10个,即为(α0,β0,φ0,fc,q1,q2,q3,p1,p2,p3),用参数向量来表示为:
根据共线性公式有:
假设星点采集个数为m,联合x和y方向的偏差和敏感矩阵,假设:
这里p是由x和y方向残留偏差组成的向量,M为A和B两个敏感矩阵组成的整体敏感矩阵。
于是有迭代方程为:
这里k为迭代序号。k取5~20,例如取10,迭代结束后得到稳定的数据值,这时的参数即为最后的校准结果。
仿真及结果分析。
仿真的星敏感器基本参数为:
视场:12度×12度;
像素阵列:1024×1024;
像素尺寸:0.015mm×0.015mm;
焦距:73.6059mm。
假设星点质心噪声为0均值,标准差为0.05像素的高斯噪声。在实验室内部利用星场模拟器,可以进行多次测量求平均值的方法降低质心算法噪声水平。但是如果数据来源是在轨校准或者夜空拍摄,那么数据只可能是一次性的,噪声的影响无法通过这种方法来降低。
假设图像中心即为主点位置,即:
x0=512×0.015mm,
y0=512×0.015mm,
星敏感器的姿态角为:
赤经ra=0度,估计误差为0.5度;
赤纬dec=0度,估计误差为0.4度;
滚转roll=0度,估计误差为0.3度;
这时模拟星场如图3所示:
设焦距偏差为0.2mm,径向畸变参数为:
q1=2e-4,q2=-4e-7,q3=1e-8,
p1=2e-4,p2=2e-4,p3=4e-6,
畸变参数的估计值均为0。
为了验证算法,首先不添加噪声,根据上面提出的算法,经过10次迭代计算可以得到:
α0 | β0 | φ0 | Δf | q1 | q2 | q3 | p1 | p2 | p3 | |
1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2 | -0.50406 | -0.39608 | -0.33367 | 0.16732 | 0.000291 | -2.56E-06 | 2.37E-08 | 4.15E-05 | 0.000324 | 0 |
3 | -0.50001 | -0.40001 | -0.29952 | 0.20123 | 0.0002 | -4.03E-07 | 1.01E-08 | 0.000202 | 0.000202 | -3.81E-05 |
4 | -0.5 | -0.4 | -0.30001 | 0.20001 | 0.0002 | -4.00E-07 | 1.00E-08 | 0.0002 | 0.0002 | 3.64E-06 |
5 | -0.5 | -0.4 | -0.3 | 0.2 | 0.0002 | -4.00E-07 | 1.00E-08 | 0.0002 | 0.0002 | 4.00E-06 |
6 | -0.5 | -0.4 | -0.3 | 0.2 | 0.0002 | -4.00E-07 | 1.00E-08 | 0.0002 | 0.0002 | 4.00E-06 |
7 | -0.5 | -0.4 | -0.3 | 0.2 | 0.0002 | -4.00E-07 | 1.00E-08 | 0.0002 | 0.0002 | 4.00E-06 |
8 | -0.5 | -0.4 | -0.3 | 0.2 | 0.0002 | -4.00E-07 | 1.00E-08 | 0.0002 | 0.0002 | 4.00E-06 |
9 | -0.5 | -0.4 | -0.3 | 0.2 | 0.0002 | -4.00E-07 | 1.00E-08 | 0.0002 | 0.0002 | 4.00E-06 |
10 | -0.5 | -0.4 | -0.3 | 0.2 | 0.0002 | -4.00E-07 | 1.00E-08 | 0.0002 | 0.0002 | 4.00E-06 |
从上表可以看出经过4次迭代,参数已经收敛于当初的假设值。
下面添加均值为0,均方差为0.05像素的质心高斯噪声。
α0 | β0 | φ0 | Δf | q1 | q2 | q3 | p1 | p2 | p3 | |
1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2 | -0.50419 | -0.39676 | -0.33283 | 0.18359 | 0.000279 | -2.45E-06 | 2.37E-08 | 4.54E-05 | 0.000313 | 0 |
3 | -0.50053 | -0.39993 | -0.30256 | 0.23552 | 0.000174 | -9.97E-08 | 9.87E-09 | 0.000213 | 0.000224 | -1.66E-03 |
4 | -0.50035 | -0.40043 | -0.30016 | 0.22547 | 0.000175 | -3.16E-08 | 8.69E-09 | 0.000216 | 0.000197 | -7.38E-04 |
5 | -0.50023 | -0.40052 | -0.29988 | 0.22108 | 0.000181 | -1.17E-07 | 8.92E-09 | 0.000209 | 0.000193 | -2.74E-04 |
6 | -0.50027 | -0.40051 | -0.29999 | 0.22201 | 0.00018 | -1.04E-07 | 8.91E-09 | 0.00021 | 0.000193 | -3.46E-04 |
7 | -0.50026 | -0.40052 | -0.29996 | 0.22179 | 0.00018 | -1.07E-07 | 8.92E-09 | 0.00021 | 0.000193 | -3.25E-04 |
8 | -0.50026 | -0.40052 | -0.29996 | 0.22183 | 0.00018 | -1.06E-07 | 8.92E-09 | 0.00021 | 0.000193 | -3.29E-04 |
9 | -0.50026 | -0.40052 | -0.29996 | 0.22182 | 0.00018 | -1.07E-07 | 8.92E-09 | 0.00021 | 0.000193 | -3.28E-04 |
10 | -0.50026 | -0.40052 | -0.29996 | 0.22182 | 0.00018 | -1.06E-07 | 8.92E-09 | 0.00021 | 0.000193 | -3.28E-04 |
同样噪声水平下,根据估计参数产生一组数据,对应于测量获得的真实数据,计算统计偏差得到x方向0.0737像素,y方向为0.0744像素。姿态估计误差角的统计值为2.4角秒,由于仿真星图中星点分布不均匀,同时星点个数不多,导致该误差角偏大,实际的姿态误差估计精度在1个多角秒。对比现有的在轨校准方法,本校准方法的精度提高了约1倍。
Claims (1)
1、一种基于星场的星敏感器校准方法,其特征在于,
1.1、建立星敏感器姿态转换矩阵;
1.1.1、建立天球坐标系;以地球中心O′为坐标原点,过春分点的轴线为Xn轴,过北极点的轴线为Zn轴,Yn轴则定义为垂直于XnZn平面的直线;
1.1.2、建立星敏感器坐标系O-XYZ,以主点为原点建立星敏感器靶面坐标系,Z轴为主光轴,X轴和Y轴分别对应星敏感器图像采集的行方向和列方向;星敏感器的姿态角由赤经α0、赤纬β0、和滚转角φ0组成,α0为Z轴在XnYn面上的投影与Xn轴的夹角,从Xn轴起逆时针计算;β0为Z轴与它在XnYn面上的投影之间的夹角,从投影逆时针计算;φ0为Zn轴在XY面上的投影与Yn轴的夹角,从投影起顺时针计算;φ0为天球坐标系子午面和像平面的交线和像面Y轴之间的夹角;
1.1.3、建立姿态转换矩阵;天球坐标系经过三次旋转即可到达星敏感器坐标系,旋转过程为:第一次绕Zn轴旋转
使得Xn轴和子午面相垂直;第二次绕旋转后的Xn轴旋转
使得Zn轴与Z轴重合;第三次绕两次旋转后的Zn轴旋转φ0,则天球坐标系O′-XnYnZn和星敏感器坐标系O-XYZ重合;设姿态转换矩阵为M,则有:
这里:
a1=sinα0cos0-cosα0sinβ0sin0;
a2=-sinα0sin0-cosα0sinβ0cos0;
a3=-cosα0cos0;
b1=-cosα0cos0-sinα0sinβ0sin0;
b2=cosα0sin0-sinα0sinβ0cos0;
b3=-sinα0cos0;
c1=cosβ0sin0;
c2=cosβ0cos0;
c3=-sinβ0;
设星场中有n颗星,第i颗星的惯性坐标为(αi,βi),i=1,...,n,则第i颗星的方向矢量为:
这里n1,n2,n3为星光向量在天球坐标系3个轴的投影;
假设焦距为fc,星敏感器成像的共线性公式为:
这里
(x0,y0)为主点坐标,(x,y)为像平面坐标;
1.2、建立星敏感器畸变模型;
假设dx,dy为x方向和y方向的畸变偏差,有:
这里,
r2=
x2+
y2;
q1,q2,q3为径向畸变系数;
p1,p2,p3为偏心畸变系数;
这里内部参数有9个,即为x0,y0,fc,q1,q2,q3,p1,p2,p3,整体成像的共线性公式为:
1.3、最小二乘参数估计;
主点偏差采用地面自准直校准结果;在星场校准过程中,假设主点位置为已知,这样总共需要考虑的校准参数有10个,即为α0,β0,φ0,fc,q1,q2,q3,p1,p2,p3,用参数向量来表示为:
根据共线性公式有:
这里A、B为敏感矩阵,求出
的各项偏导数即可获得;
假设星点采集个数为m,联合x和y方向的偏差和敏感矩阵,假设:
这里p是由x和y方向残留偏差组成的向量,M为A和B两个敏感矩阵组成的整体敏感矩阵;
于是有迭代方程为:
这里k为迭代序号,k取5~20,迭代结束后得到稳定的数据值,这时的参数即为最后的校准结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005101125522A CN100348460C (zh) | 2005-10-12 | 2005-10-12 | 一种基于星场的星敏感器校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005101125522A CN100348460C (zh) | 2005-10-12 | 2005-10-12 | 一种基于星场的星敏感器校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1948085A true CN1948085A (zh) | 2007-04-18 |
CN100348460C CN100348460C (zh) | 2007-11-14 |
Family
ID=38017736
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005101125522A Expired - Fee Related CN100348460C (zh) | 2005-10-12 | 2005-10-12 | 一种基于星场的星敏感器校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100348460C (zh) |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101231169B (zh) * | 2008-01-31 | 2010-06-09 | 北京控制工程研究所 | 紫外月球敏感器自主积分时间调节方法 |
CN101957203A (zh) * | 2010-06-07 | 2011-01-26 | 哈尔滨工业大学 | 一种星敏感器高精确度的星跟踪方法 |
CN101699222B (zh) * | 2009-11-02 | 2011-04-20 | 中国人民解放军国防科学技术大学 | 一种星敏感器校准装置及应用于高精度星敏感器校准的方法 |
CN102521506A (zh) * | 2011-12-09 | 2012-06-27 | 中国人民解放军第二炮兵装备研究院第五研究所 | 数字天顶仪旋转轴解算方法 |
CN102514734A (zh) * | 2011-10-27 | 2012-06-27 | 北京航空航天大学 | 基于日地系统Halo轨道探测器构型与姿态指向的姿态递推方法 |
CN102607597A (zh) * | 2012-03-08 | 2012-07-25 | 清华大学 | 星敏感器的三轴精度表述与测量方法 |
CN101758934B (zh) * | 2010-01-29 | 2012-08-22 | 航天东方红卫星有限公司 | 基于任务规划的星敏感器安装角度确定方法 |
CN103308074A (zh) * | 2013-06-28 | 2013-09-18 | 上海新跃仪表厂 | 一种基于双星敏在轨数据的精度分析方法 |
CN103344872A (zh) * | 2013-06-28 | 2013-10-09 | 上海新跃仪表厂 | 一种星敏安装极性的测试方法 |
CN103438907A (zh) * | 2013-09-11 | 2013-12-11 | 哈尔滨工业大学 | 一种星敏感器六自由度像平面误差的在轨标定方法 |
CN107367235A (zh) * | 2016-05-13 | 2017-11-21 | 南京理工大学 | 一种红外面阵扫描系统的轴系误差标定方法 |
CN107402023A (zh) * | 2017-08-10 | 2017-11-28 | 上海航天控制技术研究所 | 地球静止轨道卫星星敏慢变误差预报校正的方法 |
CN108106634A (zh) * | 2017-11-21 | 2018-06-01 | 中国人民解放军国防科技大学 | 一种直接观星的星敏感器内参数标定方法 |
CN108225276A (zh) * | 2018-01-11 | 2018-06-29 | 中国人民解放军战略支援部队航天工程大学 | 一种单星成像目标运动特性反演方法及系统 |
CN109506645A (zh) * | 2018-12-13 | 2019-03-22 | 上海航天控制技术研究所 | 一种星敏感器安装矩阵地面精测方法 |
CN111044074A (zh) * | 2019-12-03 | 2020-04-21 | 上海航天控制技术研究所 | 基于外场观星的星敏感器标定装置与观星标定方法 |
CN113218577A (zh) * | 2021-04-15 | 2021-08-06 | 中国科学院光电技术研究所 | 一种星敏感器星点质心位置精度的外场测量方法 |
CN114030646A (zh) * | 2021-11-22 | 2022-02-11 | 中国空间技术研究院 | 一种具有指向自感知能力的星载设备支架 |
CN114370866A (zh) * | 2021-11-23 | 2022-04-19 | 上海航天控制技术研究所 | 一种星敏感器主点和主距测量系统及方法 |
CN114030646B (zh) * | 2021-11-22 | 2024-06-11 | 中国空间技术研究院 | 一种具有指向自感知能力的星载设备支架 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102279002B (zh) * | 2011-06-27 | 2013-06-05 | 哈尔滨工业大学 | 星敏感器测量坐标系与载体坐标系转换矩阵的标定方法 |
CN102288201B (zh) * | 2011-07-06 | 2012-08-22 | 清华大学 | 用于星敏感器的精度测量方法 |
CN102410844B (zh) * | 2011-08-12 | 2013-12-11 | 北京航空航天大学 | 一种高动态星敏感器图像非均匀校正方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6108593A (en) * | 1997-07-09 | 2000-08-22 | Hughes Electronics Corporation | Method and apparatus for estimating attitude sensor bias in a satellite |
US6691033B1 (en) * | 2000-07-26 | 2004-02-10 | Hughes Electronics Corporation | System and method for calibrating inter-star-tracker misalignments in a stellar inertial attitude determination system |
CN100386593C (zh) * | 2004-11-23 | 2008-05-07 | 清华大学 | 双视场星敏感器及利用其进行星图识别的方法 |
-
2005
- 2005-10-12 CN CNB2005101125522A patent/CN100348460C/zh not_active Expired - Fee Related
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101231169B (zh) * | 2008-01-31 | 2010-06-09 | 北京控制工程研究所 | 紫外月球敏感器自主积分时间调节方法 |
CN101699222B (zh) * | 2009-11-02 | 2011-04-20 | 中国人民解放军国防科学技术大学 | 一种星敏感器校准装置及应用于高精度星敏感器校准的方法 |
CN101758934B (zh) * | 2010-01-29 | 2012-08-22 | 航天东方红卫星有限公司 | 基于任务规划的星敏感器安装角度确定方法 |
CN101957203A (zh) * | 2010-06-07 | 2011-01-26 | 哈尔滨工业大学 | 一种星敏感器高精确度的星跟踪方法 |
CN102514734A (zh) * | 2011-10-27 | 2012-06-27 | 北京航空航天大学 | 基于日地系统Halo轨道探测器构型与姿态指向的姿态递推方法 |
CN102514734B (zh) * | 2011-10-27 | 2013-11-27 | 北京航空航天大学 | 基于日地系统Halo轨道探测器构型与姿态指向的姿态递推方法 |
CN102521506A (zh) * | 2011-12-09 | 2012-06-27 | 中国人民解放军第二炮兵装备研究院第五研究所 | 数字天顶仪旋转轴解算方法 |
CN102521506B (zh) * | 2011-12-09 | 2015-01-07 | 中国人民解放军第二炮兵装备研究院第五研究所 | 数字天顶仪旋转轴解算方法 |
CN102607597B (zh) * | 2012-03-08 | 2014-05-21 | 清华大学 | 星敏感器的三轴精度表述与测量方法 |
CN102607597A (zh) * | 2012-03-08 | 2012-07-25 | 清华大学 | 星敏感器的三轴精度表述与测量方法 |
CN103344872B (zh) * | 2013-06-28 | 2015-11-18 | 上海新跃仪表厂 | 一种星敏安装极性的测试方法 |
CN103344872A (zh) * | 2013-06-28 | 2013-10-09 | 上海新跃仪表厂 | 一种星敏安装极性的测试方法 |
CN103308074B (zh) * | 2013-06-28 | 2015-11-18 | 上海新跃仪表厂 | 一种基于双星敏在轨数据的精度分析方法 |
CN103308074A (zh) * | 2013-06-28 | 2013-09-18 | 上海新跃仪表厂 | 一种基于双星敏在轨数据的精度分析方法 |
CN103438907A (zh) * | 2013-09-11 | 2013-12-11 | 哈尔滨工业大学 | 一种星敏感器六自由度像平面误差的在轨标定方法 |
CN103438907B (zh) * | 2013-09-11 | 2016-01-20 | 哈尔滨工业大学 | 一种星敏感器六自由度像平面误差的在轨标定方法 |
CN107367235B (zh) * | 2016-05-13 | 2019-12-27 | 南京理工大学 | 一种红外面阵扫描系统的轴系误差标定方法 |
CN107367235A (zh) * | 2016-05-13 | 2017-11-21 | 南京理工大学 | 一种红外面阵扫描系统的轴系误差标定方法 |
CN107402023A (zh) * | 2017-08-10 | 2017-11-28 | 上海航天控制技术研究所 | 地球静止轨道卫星星敏慢变误差预报校正的方法 |
CN107402023B (zh) * | 2017-08-10 | 2020-06-23 | 上海航天控制技术研究所 | 地球静止轨道卫星星敏慢变误差预报校正的方法 |
CN108106634B (zh) * | 2017-11-21 | 2020-06-16 | 中国人民解放军国防科技大学 | 一种直接观星的星敏感器内参数标定方法 |
CN108106634A (zh) * | 2017-11-21 | 2018-06-01 | 中国人民解放军国防科技大学 | 一种直接观星的星敏感器内参数标定方法 |
CN108225276A (zh) * | 2018-01-11 | 2018-06-29 | 中国人民解放军战略支援部队航天工程大学 | 一种单星成像目标运动特性反演方法及系统 |
CN109506645A (zh) * | 2018-12-13 | 2019-03-22 | 上海航天控制技术研究所 | 一种星敏感器安装矩阵地面精测方法 |
CN109506645B (zh) * | 2018-12-13 | 2020-09-18 | 上海航天控制技术研究所 | 一种星敏感器安装矩阵地面精测方法 |
CN111044074A (zh) * | 2019-12-03 | 2020-04-21 | 上海航天控制技术研究所 | 基于外场观星的星敏感器标定装置与观星标定方法 |
CN113218577A (zh) * | 2021-04-15 | 2021-08-06 | 中国科学院光电技术研究所 | 一种星敏感器星点质心位置精度的外场测量方法 |
CN114030646A (zh) * | 2021-11-22 | 2022-02-11 | 中国空间技术研究院 | 一种具有指向自感知能力的星载设备支架 |
CN114030646B (zh) * | 2021-11-22 | 2024-06-11 | 中国空间技术研究院 | 一种具有指向自感知能力的星载设备支架 |
CN114370866A (zh) * | 2021-11-23 | 2022-04-19 | 上海航天控制技术研究所 | 一种星敏感器主点和主距测量系统及方法 |
CN114370866B (zh) * | 2021-11-23 | 2024-01-30 | 上海航天控制技术研究所 | 一种星敏感器主点和主距测量系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100348460C (zh) | 2007-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1948085A (zh) | 一种基于星场的星敏感器校准方法 | |
US8934721B2 (en) | Microscopic vision measurement method based on adaptive positioning of camera coordinate frame | |
CN107564057B (zh) | 顾及大气折光校正的高轨面阵光学卫星在轨几何标定方法 | |
CN108198219B (zh) | 用于摄影测量的相机标定参数的误差补偿方法 | |
CN1851408A (zh) | 基于多天体路标的星际巡航自主导航方法 | |
CN100376883C (zh) | 一种基于像素频率的星敏感器校准方法 | |
CN108765484B (zh) | 基于两台高速摄像机的活体昆虫运动采集及数据重构方法 | |
CN1851752A (zh) | 三维重构系统中双摄像机标定方法 | |
CN109341720A (zh) | 一种基于恒星轨迹的遥感相机几何标定方法 | |
CN1939807A (zh) | 基于weng模型的星敏感器在轨校准方法 | |
CN113538595B (zh) | 利用激光测高数据辅助提升遥感立体影像几何精度的方法 | |
CN112232319A (zh) | 一种基于单目视觉定位的扫描拼接方法 | |
CN114758011B (zh) | 融合离线标定结果的变焦相机在线标定方法 | |
CN107481281B (zh) | 相对位姿计算方法和装置及航天飞行器交会对接系统 | |
CN115546289A (zh) | 一种基于机器人的复杂结构件三维形貌测量方法 | |
JP2018531448A5 (zh) | ||
Pi et al. | On-orbit geometric calibration using a cross-image pair for the linear sensor aboard the agile optical satellite | |
Han et al. | A direct determination of the orientation parameters in the collinearity equations | |
CN103438907A (zh) | 一种星敏感器六自由度像平面误差的在轨标定方法 | |
CN111998823B (zh) | 基于双目异光源测距装置的目标测距方法 | |
JP5349273B2 (ja) | Rpc計算装置 | |
CN1851618A (zh) | 单目视觉半实物仿真系统及方法 | |
CN111815712A (zh) | 一种高精度的相机-单激光仪联合标定方法 | |
CN117078769A (zh) | 一种电子后视镜cms摄像头外参数标定方法及系统 | |
CN1923621A (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20071114 Termination date: 20101012 |