一种星载图像辅助导航方法技术

一种星载图像辅助导航方法，属于航天器姿态确定技术领域。本发明专利技术是为了解决现有的图像处理辅助姿态确定方法中，由于数据量较大造成处理速度慢，硬件资源消耗较大等问题。本发明专利技术所述的方法，首先，对在轨图像进行压缩并计算图像梯度值；加载模板图像，根据压缩过程中获取的黑塞矩阵，进行多尺度的黑塞矩阵特征值求取以及多尺度金字塔运算，并根据最大值准则得到特征点；然后，根据在轨图像的特征点和模板图像的特征点，采用随机抽样一致性算法进行特征点匹配，从而得到两个图像之间的仿射变换矩阵；最后建立相机运动模型，将图像偏移量转换为姿态角变化量。本发明专利技术可应用于微小卫星对地观测系统。

A satellite borne image aided navigation method

A satellite borne image aided navigation method belongs to the technology field of spacecraft attitude determination. The purpose of the invention is to solve the problem of the existing image processing aided attitude determination method, which is caused by the large amount of data, slow processing speed and large hardware resource consumption. Methods the first on orbit image compression and image gradient calculation; loading the template image, according to the Hesse matrix in the process of compression, strike and multi-scale Pyramid operations Hesse matrix eigen multi-scale, and the maximum value is obtained according to the quasi feature points; then, according to the feature points the feature points on orbit image and template image, using the random consistency algorithm for feature point matching, so as to obtain the affine transformation matrix between two images; finally a camera motion model, the image offset for attitude change. The invention can be applied to the earth observation system of small satellite.

【技术实现步骤摘要】
一种星载图像辅助导航方法
本专利技术涉及航天器姿态确定
，具体涉及星载图像辅助导航方法。
技术介绍
随着遥感卫星分辨率越来越高，高分辨率光学载荷搭载在微小卫星上能有效降低卫星发射成本，并进一步实现卫星组网、组件更换，实现商业化运作等。微小遥感卫星相较于现有的传统卫星同时存在一些新的问题，比如卫星在轨对地指向精度较低，姿态稳定性较差等。这些问题会降低遥感图像的分辨率，甚至影响卫星的正常工作。针对这些问题目前的解决方案是：1)利用高精度的姿态传感器，比如基于磁流体效应的角速度传感器以及基于流体旋转差动效应的角位移传感器。但是此类传感器成本较高且技术难度较大。2)利用低成本的惯性传感器和星敏感器进行姿态融合和姿态确定。但是现有传感器短时精度无法达到要求，无法修正卫星拍摄过程中因为卫星颤振引起的图像抖动和模糊等问题。3)利用高精度遥感图像包含丰富的地理信息，通过在轨图像处理算法可实现低成本高效率的姿态辅助确定。然而由于遥感图像数据量较大，星上计算能力有限，造成处理速度慢，硬件资源消耗较大，因此需要一种改进的匹配策略实现在轨姿态确定。
技术实现思路
本专利技术是为了解决现有的图像处理辅助姿态确定方法中，由于数据量较大造成处理速度慢，硬件资源消耗较大等问题，现提供一种星载图像辅助导航方法。本专利技术所述一种星载图像辅助导航方法，通过以下技术方案实现：步骤一、对在轨图像进行图像压缩，在压缩过程中计算图像梯度值；步骤二、加载模板图像，根据压缩过程中获取的不同尺度下的图像梯度值的黑塞矩阵，进行多尺度的黑塞矩阵特征值求取以及多尺度金字塔运算，并通过最大值准则得到特征点；步骤三、根据在轨图像的特征点和模板图像的特征点，采用随机抽样一致性算法进行特征点匹配，得到两个图像之间的仿射变换矩阵；然后建立相机运动模型，将图像偏移量转换为姿态角变化量。本专利技术与现有技术相比较，最为突出的特点和显著的有益效果是：本专利技术所述的星载图像辅助导航方法，是一种通过获取多尺度下图像特征点进行图像配准获得姿态偏移信息的方法，充分利用了压缩图像提供的图像信息，相比于现有的SURF算法，该方法无需进行多尺度缩放操作，无需单独求解梯度矩阵，因此，本专利技术的方法简单，效率高，硬件资源消耗较低，节省计算资源，且易于实现，经过计算，在轨道变化可以忽略的情况下像素偏移0.3个单位的像素会导致姿态角偏移量为1.125μrad。附图说明图1为改进的JPEG2000压缩示意图；图2为多尺度图像特征点获取示意图；图3为本专利技术所述的利用在轨图像辅助导航方法的流程图。具体实施方式具体实施方式一：参照图3具体说明本实施方式，本实施方式给出的一种星载图像辅助导航方法，具体是按照以下步骤制备的：步骤一、首先，对在轨图像进行图像压缩，在压缩过程中计算图像梯度值；步骤二、加载事先保存好的模板图像，根据压缩过程中获取的不同尺度下的图像梯度值的黑塞矩阵，进行多尺度的黑塞矩阵特征值求取以及多尺度金字塔运算，并通过最大值准则(即特征点的获取是根据不同层之间一些点的最大值得到的)得到特征点；步骤三、根据在轨图像的特征点和模板图像的特征点，采用随机抽样一致性算法进行特征点匹配，从而得到两个图像之间的仿射变换矩阵(一种二维坐标到二维坐标之间的线性变换)；然后建立相机运动模型，将图像偏移量转换为姿态角变化量，同时考虑轨道变化对随机抽样一致性算法的影响。具体实施方式二：参照图1说明本实施方式，本实施方式与具体实施方式一不同的是：步骤一的具体过程如下：在FPGA在轨图像处理平台中，在图像经过电子学处理送入平台进行压缩时，进行5/3提升式小波变换，同时在小波变换时实施图像梯度值的求取；为进一步降低图像噪声对梯度值的影响，提高梯度值的稳定性，利用中心点以及与其上下左右邻域的5个值进行图像梯度值的求取，一维变换公式如下：其中y(n)代表压缩后的下一级图像像素值，x(n)指当前级图像像素值，Dx代表图像横向的梯度值，Dy的表达式与Dx相同，只是在计算方向上不同，Dy代表图像纵向的梯度值；Dxy代表了图像在45度方向上的梯度值，y1是中间变量，n指图像像素的索引值；而公式(1)与公式(3)因为具有相似的结构，因此公式(3)在FPGA平台的计算中可借用公式(1)的计算资源；从而有效提高了算法的运行效率；因梯度图像同原图像一致，因此在第一级图像压缩之后，获得的梯度值与图像高频分量一同输入到公共存储区中进行保存。为提高后续特征点匹配的效率，梯度信息将不参与图像压缩的量化以及编码运算。具体实施方式三：参照图2说明本实施方式，本实施方式与具体实施方式二不同的是：步骤二中对压缩图像进行多尺度特征点提取的具体步骤为：针对获取的图像梯度值，根据特征值原理，通过对梯度特征值的求取进行特征点判决；计算在轨图像上的像素点I在压缩图像第m层的黑塞矩阵为：其中，Dx(I,m)是在轨图像上的像素点I在压缩图像第m层横向的梯度值，Dy(I,m)是在轨图像上的像素点I在压缩图像第m层纵向的梯度值Dxy(I,m)是在轨图像上的像素点I在压缩图像第m层在45度方向上的梯度值；利用黑塞矩阵H的行列式进行特征值的判别：det(H)＝DxDy-wDxy2(7)其中w为近似系数，用于补偿判决式误差；det是一个计算机函数，该函数用于求矩阵的行列式；根据行列式性质，公式(7)取最大值时可以认为对应的像素点是像素点阵附近的特征点，为进一步提高特征点获取稳定性，采用如下判决式进行特征点筛选：S＝max(det(H,m),Threshold)+max(det(H,m-1),det(H,m),det(H,m+1))(8)即，首先利用预设极值进行筛选，小于预设极值的点都被抛弃，之后再利用周围的8个点以及相邻上下两层的周围9个点进行判决，满足最大值的点才会被认为强特征点；Threshold指的是预设极值，S指的是阈值；图2填充区域表现了筛选过程，其中深色填充像素代表目标像素，而浅色填充像素代表周围要比较的像素点。另外，由于在轨遥感图像会存在进光量不足，导致特征点较小的问题，因此这里的阈值应为自适应阈值，而将特征点数目作为调整目标值，当提取的特征点少于预设的特征点数目S*时，通过调整阈值以增加特征点，直到总的特征点数目符合要求；公式表述如下：Threshold＝argmin(S>＝S*,Threshold)(9)。其它步骤及参数与具体实施方式二相同。具体实施方式四：本实施方式与具体实施方式三不同的是：步骤三中采用随机抽样一致性算法进行特征点匹配的具体步骤为：根据在轨图像和模板图像之间存在的仿射变换，即针对模板图像中的任一像素点I′＝(a′,b′)和对于在轨图像中的像素点I＝(a,b)，存在如下对应关系：其中(a′,b′)为像素点I′的平面坐标，(a,b)为像素点I的平面坐标，A代表旋转矩阵，B代表平移矩阵，可用如下参数表示：式中θ表示旋转角度，b1、b2表示平移的像素大小；在(10)式中，最少需要三组(a,b)像素点才能解出仿射变换矩阵；然后，利用随机抽样一致性算法解出仿射矩阵。具体实施方式五：本实施方式与具体实施方式四不同的是：步骤三中随机抽样一致性算法解出仿射变换矩阵的具体步骤为：1)计算特征向量描述符的欧式距离，之后利用最小值比次小值原理进行特征点匹配本文档来自技高网...

【技术保护点】
一种星载图像辅助导航方法，其特征在于：所述星载图像辅助导航方法包括以下步骤：步骤一、对在轨图像进行图像压缩，在压缩过程中计算图像梯度值；步骤二、加载模板图像，根据压缩过程中获取的不同尺度下的图像梯度值的黑塞矩阵，进行多尺度的黑塞矩阵特征值求取以及多尺度金字塔运算，并通过最大值准则得到特征点；步骤三、根据在轨图像的特征点和模板图像的特征点，采用随机抽样一致性算法进行特征点匹配，得到两个图像之间的仿射变换矩阵；然后建立相机运动模型，将图像偏移量转换为姿态角变化量。

【技术特征摘要】
1.一种星载图像辅助导航方法，其特征在于：所述星载图像辅助导航方法包括以下步骤：步骤一、对在轨图像进行图像压缩，在压缩过程中计算图像梯度值；步骤二、加载模板图像，根据压缩过程中获取的不同尺度下的图像梯度值的黑塞矩阵，进行多尺度的黑塞矩阵特征值求取以及多尺度金字塔运算，并通过最大值准则得到特征点；步骤三、根据在轨图像的特征点和模板图像的特征点，采用随机抽样一致性算法进行特征点匹配，得到两个图像之间的仿射变换矩阵；然后建立相机运动模型，将图像偏移量转换为姿态角变化量。2.根据权利要求1所述的一种星载图像辅助导航方法，其特征在于：步骤一的具体过程如下：在FPGA在轨图像处理平台中，在图像经过电子学处理送入平台进行压缩时，进行5/3提升式小波变换，同时在小波变换时实施图像梯度值的求取；利用中心点邻域的值进行图像梯度值的求取：其中y(n)为压缩后的下一级图像像素值，x(n)指当前级图像像素值，Dx是图像横向的梯度值，Dy是图像纵向的梯度值；Dxy是图像在45度方向上的梯度值，y1是中间变量，n指图像像素的索引值；在第一级图像压缩之后，获得的梯度值与图像高频分量一同输入到公共存储区中进行保存。3.根据权利要求2所述的一种星载图像辅助导航方法，其特征在于：步骤二中对压缩图像进行多尺度特征点提取的具体步骤为：针对获取的图像梯度值，根据特征值原理，通过对梯度特征值的求取进行特征点判决；计算在轨图像上的像素点I在压缩图像第m层的黑塞矩阵为：其中，Dx(I,m)是在轨图像上的像素点I在压缩图像第m层横向的梯度值，Dy(I,m)是在轨图像上的像素点I在压缩图像第m层纵向的梯度值Dxy(I,m)是在轨图像上的像素点I在压缩图像第m层在45度方向上的梯度值；利用黑塞矩阵H的行列式进行特征值的判别：det(H)＝DxDy-wDxy2(7)其中w为近似系数，公式(7)取最大值时对应的像素点是像素点阵附近的特征点，采用如下判决式进行特征点筛选：S＝max(det(H,m),Threshold)+max(det(H,m-1),det(H,m),det(H,m+1))(8)Threshold指的是预设极值，S指的是阈值；阈值为自适应的，将特征点数目作为调整目标值，当提取的特征点少于预设的特征点数目S*时，通过调整阈值S以增加特征点，直到总的特征点数目符合要求；公式表述如下：Threshold＝argmin(S>＝S*,Threshold)(9)。4.根据权利要求3所述的一种星载图像辅助导航方法，其特征在于：步骤三中采用随机抽样一致性算法进行特征点匹配的具体步骤为：根据在轨图像和模板图像之间存在的仿射变换，即针对模板图像中的任一像素点I′＝(a′,b′)和对于在轨图像中的像素点I＝(a,b)，存在如下对应关系：其中(a′,b′)为I′的平面坐标，(a,b)为I的平面坐标，A代表旋转矩阵，B代表平移矩阵，用如下参数表示：式中θ表示旋转角度，b1、b2表示平移的像素大小；然后用随机抽样一致性算法解出仿射变换矩阵。5.根据权利要求4所述的一种星载图像辅助导航方法，其特征在于：步骤三中用随机抽样一致性算法...

【专利技术属性】
技术研发人员：徐国栋，张兆祥，刘明，邢雷，王梓霖，张光宇，朱晏辰，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江,23