大视场星敏感器光学系统畸变自动标定方法技术方案

本发明专利技术涉及一种大视场星敏感器光学系统畸变自动标定方法，包括：设定星敏感器的成像分辨率，像元角分辨率、主点、分区步长初值、最大畸变和各分区内定标点个数；以S(1)为分区步长，计算中心区域内定标点的二维转台编码器角度信息，计算该级分区内所有畸变标定点的最大畸变；当D_max(1)≤D

Automatic calibration method of distortion in optical system of star sensor with large field of view

【技术实现步骤摘要】
大视场星敏感器光学系统畸变自动标定方法
本专利技术涉及一种大视场星敏感器光学系统畸变自动标定方法。
技术介绍
随着天文导航技术的发展，大视场星敏感器成为星敏感器研究中重要发展方向。众所周知，大视场星敏感器光学系统畸变较大，光学畸变的定标效果直接影响星敏感器姿态测量的精度，随着大视场光学测量系统应用领域的不断扩展，对相对光学畸变的要求甚至低于0.01％。光学畸变标定主要包括基于整幅图像的畸变标定和分区畸变标定，其中分区畸变标定能够统一标定各区域内的焦距误差、主点误差、像面旋转、倾斜误差以及光学畸变等参数，是一种有效的畸变标定方法。目前星敏感器分区光学畸变标定一般是将星敏感器成像区域按预设分区数量，以星敏感器主点为中心区域的中心，以固定间隔将成像区域均匀分成多个区域，通过拟合各分区内的定标点数据，实现分区畸变标定。通常情况下均匀分区存在如下问题：1.离主点越远的区域畸变越大，均匀分区无法保证畸变标定的整体精度；2.分区数量越多，畸变定标精度越高，对于对畸变标定精度较高的应用，虽然可以采用更多的分区提高畸变定标精度，分区数量不足时需重新采集畸变标定数据，造成大量的重复工作。
技术实现思路
有鉴于此，有必要提供一种大视场星敏感器光学系统畸变自动标定方法。本专利技术提供一种大视场星敏感器光学系统畸变自动标定方法，该方法包括如下步骤：a.设定星敏感器的成像分辨率M×M，像元角分辨率θ″、主点(X0,Y0)、分区步长初值S(1)、最大畸变D0角秒和各分区内定标点个数Xn×Xn，Xn为奇数；b.以S(1)为分区步长，以星敏感器主点(X0,Y0)为中心，计算中心区域内定标点的二维转台编码器角度信息，拟合星敏感器畸变修正系数，并计算该级分区内所有畸变标定点的最大畸变D_max(1)；c.当满足D_max(1)≤D0时，设第N级分区的步长细分份数，其中N≥2；d.计算第N级分区步长，计算各分区内二维转台的角度信息，采集平行光管星点图像并提取不同角度对应的星点坐标，拟合各区域畸变修正参数，计算该级分区内最大畸变D_max(N)，以及计算已定标区域的边长Stot；e.当满足D_max(N)≤D0时，输出该级分区内各区域畸变修正系数；f.当满足Stot≥M时，则现有定标区域边缘接近星敏感器成像区域边界，星敏感器光学系统畸变标定过程结束。其中，所述的步骤c还包括：当不满足D_max(1)≤D0时，缩小分区步长，以S(N)为新分区步长重新划分中心区域，然后转入步骤b。所述的步骤e还包括：当不满足D_max(N)≤D0时，返回步骤d。所述的步骤f还包括：当不满足Stot≥M时，N＝N+1，而后返回步骤d。所述的步骤b具体包括：以第1级分区S(1)为分区步长，以星敏感器主点(X0,Y0)为中心，划分中心区域，根据分区内定标点分辨率，设已定标区域的边长为Stot，Stot＝S(N)，计算中心区域内定标点的二维转台编码器角度信息，如下式所示：(A(N,m,z),E(N,m,z))＝(1)(A0+hA·(Stot-SN)·θ/2+i·S(N)·θ/Xn,E0+hB·(Stot-SN)·θ/2+j·S(N)·θ/Xn)式中：N为分区级数，m为该级分区内区域编号，第1级分区N＝1，m＝1，z∈{1，2，...，Xn2}，hA和hB为子分区处于中心点的方向，子区域中心点hA和hB取值为0，在编码器减小方向hA和hB取值为-1，在编码器增加方向hA和hB取值为1；将二维转台角度值依次设置公式(1)角度，采集平行光管星点图像并提取星点坐标，记为(X(N,m,z),Y(N,m,z))；根据当前二维转台编码器角度信息、主点对应的编码器信息及系统焦距计算星点成像的理论位置(Xl,Yl)，如下式所示：式中，Al，El为星点位于主点位置的编码器角度值，A0，E0是畸变标定过程中编码器角度值，Xl，Yl是畸变标定过程中星点理论成像位置，f为焦距；设星敏感器畸变修正系数为k1，k2，k3，k4，k5，k6，其中k1，k2，k3为星点坐标X方向修正系数，k4，k5，k6为星点坐标Y方向修正系数，采用最小二乘拟合方法，所有系数可由公式求解：假定采集星点坐标为X和Y，修正后的星点理论成像位置为Xp和Yp，则有修正公式：将该区域内所有星点坐标代入修正公式，可计算该级分区内所有畸变标定点的最大畸变，以实际星点坐标与理想星点坐标的欧氏距离定义最大畸变：所述的步骤d具体包括：利用公式(7)计算第N级分区步长，在(N-1)级分区周围划分(4*t+4)个区域，以公式(1)计算各分区内二维转台的角度信息(A(N,m,z),E(N,m,z))，采集平行光管星点图像并提取不同角度对应的星点坐标(X(N,m,z),Y(N,m,z))，利用公式(3)和公式(4)拟合各区域畸变修正参数，利用公式(6)计算该级分区内最大畸变D_max(N)；第N(N≥2)级分区：第N级分区初始步长：步骤d中所述计算已定标区域的边长Stot具体包括：Stot＝S(N)*(t+2)。所述的步骤e具体包括：输出该级分区内各区域畸变修正系数[k1(N,m),k2(N,m),k3(N,m),k4(N,m),k5(N,m),k6(N,m)]。本专利技术利用二维转台和平行光管，以星敏感器主点为中心，不断向成像区域边缘扩展畸变定标区域，以最大畸变为依据，自动搜索不同分级区域的步长，并采用最小二乘拟合计算畸变定标参数，实现大视场星敏感器光学畸变自动分区定标。本专利技术定标时可根据系统允许的最大畸变，自动对星敏感器成像区域进行分区，避免了以固定步长分区畸变标定方法边缘区域畸变标定精度差的问题，保证大视场星敏感器光学系统畸变标定的整体精度，同时避免了分区步长选择不合理造成的重复定标工作量，为大视场星敏感器以及大口径成像系统的高精度定标提供了一种简单有效的途径。附图说明图1为本专利技术大视场星敏感器光学系统畸变自动标定方法的流程图。具体实施方式下面结合附图及具体实施例对本专利技术作进一步详细的说明。参阅图1所示，是本专利技术大视场星敏感器光学系统畸变自动标定方法较佳实施例的作业流程图。步骤S1，设定星敏感器的成像分辨率、像元角分辨率、主点、分区步长初值、最大畸变和各分区内定标点个数。具体而言：设成像分辨率为M×M，像元角分辨率为θ″，主点为(X0,Y0)，分区步长初值为S(N)(N＝1，单位pixel)，最大畸变为D0角秒，各分区内定标点分辨率为Xn×Xn(Xn为奇数)。将所述星敏感器固定于高精度二维转台上，对准水平零度光管，通过转动二维转台，使光管星点成像于探测器主点位置，设此时二维转台编码器角度为(A0″，E0″)。步骤S2，第1级分区：以S(1)为分区步长，以星敏感器主点(X0,Y0)为中心，此时已定标区域的边长Stot＝S(1)，本文档来自技高网...

【技术保护点】
1.一种大视场星敏感器光学系统畸变自动标定方法，其特征在于，该方法包括如下步骤：/na.设定星敏感器的成像分辨率M×M，像元角分辨率θ″、主点(X

【技术特征摘要】
1.一种大视场星敏感器光学系统畸变自动标定方法，其特征在于，该方法包括如下步骤：
a.设定星敏感器的成像分辨率M×M，像元角分辨率θ″、主点(X0,Y0)、分区步长初值S(1)、最大畸变D0角秒和各分区内定标点个数Xn×Xn，Xn为奇数；
b.以S(1)为分区步长，以星敏感器主点(X0,Y0)为中心，计算中心区域内定标点的二维转台编码器角度信息，拟合星敏感器畸变修正系数，并计算该级分区内所有畸变标定点的最大畸变D_max(1)；
c.当满足D_max(1)≤D0时，设第N级分区的步长细分份数，其中N≥2；
d.计算第N级分区步长，计算各分区内二维转台的角度信息，采集平行光管星点图像并提取不同角度对应的星点坐标，拟合各区域畸变修正参数，计算该级分区内最大畸变D_max(N)，以及计算已定标区域的边长Stot；
e.当满足D_max(N)≤D0时，输出该级分区内各区域畸变修正系数；
f.当满足Stot≥M时，则现有定标区域边缘接近星敏感器成像区域边界，星敏感器光学系统畸变标定过程结束。


2.如权利要求1所述的方法，其特征在于，所述的步骤c还包括：
当不满足D_max(1)≤D0时，缩小分区步长，以S(N)为新分区步长重新划分中心区域，然后转入步骤b。


3.如权利要求2所述的方法，其特征在于，所述的步骤e还包括：
当不满足D_max(N)≤D0时，返回步骤d。


4.如权利要求3所述的方法，其特征在于，所述的步骤f还包括：
当不满足Stot≥M时，N＝N+1，而后返回步骤d。


5.如权利要求4所述的方法，其特征在于，所述的步骤b具体包括：
以第1级分区S(1)为分区步长，以星敏感器主点(X0,Y0)为中心，划分中心区域，根据分区内定标点分辨率，设已定标区域的边长为Stot，Stot＝S(N)，计算中心区域内定标点的二维转台编码器角度信息，如下式所示：
(A(N,m,z),E(N,m,z))＝(1)
(A0+hA·(Stot-SN)·θ/2+i·S(N)·θ/Xn,E0+hB·(Stot-SN)·θ/2+j·S(N)·θ/Xn)
式中：N为分区级数，m为该级分区内区域编号，第1级分区N＝1，m＝1，z∈{1,2,.....

【专利技术属性】
技术研发人员：吕游，何昕，袁航飞，王军，穆治亚，何家维，张磊，魏仲慧，
申请(专利权)人：中国科学院长春光学精密机械与物理研究所，
类型：发明
国别省市：吉林;22