专利名称：基于fpga的星敏感器图像处理方法
技术领域：
本发明涉及一种星敏感器图像处理方法，并基于FPGA以并行计算的体系结构实现该算法。
背景技术：
星敏感器是空间飞行器GNC系统的关键部件，是基于计算机视觉测量理论研制的空间精密仪器，主要用于空间飞行器三轴姿态测量和空间飞行器导航。星敏感器通常由光学及精密结构系统、光电探测器及信号处理电路和软件等三部分组成。基于计算机视觉测量理论，以恒星稳定的光信号作为输入，星敏感器可以得到恒星在其本体坐标系内的精确位置，再由其内部星表存储的恒星基于地心惯性坐标系的位置信息，可以得到星敏感器本体相对于地心惯性坐标系的三轴姿态信息，从而得到空间飞行器本体相对于地心惯性系的三轴姿态信息。恒星数字图像处理是星敏感器软件的重要组成部分，主要完成恒星成像点的精确定位功能，其稳定性、可靠性和运算速度，直接关系到星敏感器的可靠性、测量精度和数据更新率等性能指标。田玉龙等在《哈尔滨工业大学学报》O005，37<8> :1068-1070)上发表论文《基于局部熵的星敏感器星图提取方法》，给出了星敏感器星点质心提取的一种算法，基于PC机仿真实现，基于“先保存星图至存储器，再提取星点”的串行方式工作，串行的工作方式很难提高计算效率。姚太雷等在《长春理工大学学报》(2008，31<1> :71-73)上发表论文《适用于星敏感器的星体识别研究》，给出了基于RISC处理器的逐行扫描星点提取算法，但RISC处理器的计算过程本质上是串行的，其运算速度较FPGA慢，且很难与其他图像处理算法实现并行计算。郝雪涛等在《北京航空航天大学学报》0005，31<4> :381-384)上发表论文《CMOS 星敏感器图像驱动及实时星点定位算法》，基于FPGA实现，但只能识别单像素噪声点，且没有辅以图像高通滤波器、星点排序等算法，其方法鲁棒性方面需要改进。周国辉等在《航天控制》0006，M<1> 87-91)上发表论文《星敏感器图像处理系统的并行流水线操作研究》， 将图像处理工作分别由CPLD与DSP实现，功耗较大，且系统设计复杂。

发明内容
本发明的目的是提供基于FPGA的星敏感器图像处理方法，能够完成星点质心的精确定位。为达到上述目的，本发明是通过以下技术方案实现的，基于FPGA的星敏感器图像处理方法，其特征在于包括如下步骤(1)根据星敏感器数字图像能量分布的特征、恒星成像的特征和硬件实现数学运算的特点，设计高通滤波器，由此计算星点提取阈值，解决在强光干扰条件下影响星点提取的关键问题；(2)根据星敏感器恒星成像的特征，设计星点质心提取算法，在光电探测器逐行输出像素数据的同时，同步计算星点质心，由星点的图像特征，在运算过程中对星点与噪声点进行判别；(3)根据硬件运算的特点，采用FPGA实现基数排序算法，在完成星点提取工作后， 按星点能量值对星点集合进行排序，确保能量值高的星点子集用于全天识别和姿态解算， 以显著减小强光干扰产生的噪声对姿态计算的影响；(4)基于FPGA实现星敏感器图像处理的全部工作，并以并行运算的体系结构实现。所述的步骤(1)中，根据星敏感器数字图像能量分布的特征、恒星成像的特征和硬件实现数学运算的特点，设计16X16pix.的高通滤波器，在光电探测器逐行给出像素数据的同时，逐个区域的完成图像滤波和阈值的计算，得到阈值矩阵，使得星敏感器可以有效处理强光干扰条件下的图像。所述的步骤O)中，根据星敏感器恒星成像的特征，设计的同步提取星点质心的算法可以在逐行得到图像像素数据的同时，辨别单像素点和双像素点噪声的机制并及时舍弃，有效解决在星敏感器在实际运行过程中因噪声点过多带来的存储器溢出问题。所述的步骤(3)中，根据硬件运算的特点，设计基数排序算法的状态机，并运用于星点集合排序，在排序过程中，为减少存储空间开销，排序采用移动索引的方法进行。所述的步骤(1) (4)中，采用FPGA完成星敏感器图像处理算法的滤波、星点质心提取和排序的全部工作，并采用四级流水线方式设计并行运算体系结构。本发明采用的方法，与现有技术相比，其优点和有益效果是(1)充分利用星敏感器恒星成像及光电探测器数字图像的特点，能在强光干扰条件下仍然可以稳定运行。(2)充分利用FPGA在并行运算体系结构设计中的优势，采用四级流水组织图像处理算法中的运算模块，大幅提高运算速度和星敏感器的数据更新率。


下面将结合附图和实施例对本发明作进一步详细的说明。图1为星敏感器图像处理算法体系结构框图；图2为用于星敏感器图像处理的空间域高通滤波器计算示意图；图3为基数排序算法状态机示意图；图4为并行运算体系结构示意图。
具体实施例方式如图1所示，给出了星敏感器图像处理算法的体系结构。算法由高通滤波、星点提取和星点排序三部分组成，以并行运算方式基于FPGA实现方法的步骤。光电探测器在驱动信号集的控制下，将光信号转化为数字图像信号输出至空间域高通滤波器，用以实现图像滤波，其计算结果为阈值矩阵，作为星点质心提取的计算参数，星点质心提取的计算结果为星点集合，经星点排序后成为有序的星点集合，提供给星敏感器全天识别与姿态解算软件进行处理。如图2所示，给出了用于星敏感器图像处理的空间域高通滤波器计算示意图。星像信号对于空间位置是连续的，在足够小的邻域内，图像背景信号可以认为是近似均勻的。恒星像点、背景和噪声像点各自具有显著的视觉特征。背景信号的能量分布比较均勻，属于图像信号低频部分，星光信号的能量呈现出高频窄脉冲的特点，光噪声像素点能量值较高，但其几何位置分布分散，与星光信号能量较高，几何位置分布集中具有明显差别。设置滤波窗口大小为16X16pix.，窗口大小的确定基于如下几个因素①恒星在星敏感器上所成的像，按照星等的不同，像素数一般为数十个之内，且以弥散圆的形状分布，在 16X16pix.子区域内，星光信号相对于背景是能量较高的高频信号，易于提取；②星图的大小一般为2的整数次幂，选取大小为16X 16pix的滤波窗口，可以快速计算子区域行列索引值，避免使用除法和求余运算；③平均值可以采用移位运算，避免使用除法。采用FPGA实现高通滤波器。系统复位后，当一帧星图开始输出后，算法启动。首先由像素计数器计算出当前像素点所在子区域的行列索引值并从阈值矩阵中读取该区域的部分和，之后，将当前像素灰度值累加至子区域灰度的部分和。图2中四种颜色对应的图像子区域的灰度值累加和分别存放在阈值存储器对应颜色的地址中。在每行扫描结束后，与该行对应的阈值存储器中64个灰度部分和刷新一次，一帧扫描结束后，该帧对应的阈值矩阵同步给出。第k帧计算得到的阈值矩阵用于第k+Ι帧星点提取，第1帧采用经验值设定阈值以提取星点。由星点成像的特征，星点的定义如下。定义对于大小为HXW，最高量化值为Vmax，星点提取阈值为T的星图，若像素集合 P= Ipi(Xi^yijVi), i = 1,2, ...，N}可以作为星点，需满足①0 < Xi < W，0 < Yi < H，T < Vi < Vmax, N 彡 3 ；②对于e尺且满足Pi与Pj四邻域相邻。其中Xi，Yi为星点像素坐标，Vi为灰度值，N为星点超过阈值的像素数。算法可以实时地完成星点质心定位，并将星点质心数据保存于FPGA内部的存储器中。对于每个时刻得到的当前像素点p(i，j)，设其左侧相邻点为p(i，j-1)，左上相邻点为P(i-1，j-1)，上方相邻点为P(i-1，j)，令其对应的灰度值分别为V(i，j), v(i, j-1), v(i-l, j-1), v(i-l, j)，且令阈值为T，星点类属标志为starFlag。在算法中，需要记录每个星点类属标记对应的灰度值超过阈值的像素个数，为实现并行星点提取，根据定义1，需要记录第i_l行所有像素点的类属标记，此外还需要存储星点能量值的空间EnergyVec，存储星点各像素灰度加权坐标和的空间GrayXVec和GrayYVec，存放各星点像素数总和的空间MarSizeVec和星点类属复用标志starFlagMackOverlap，上述所有存储空间初始化为 O0表1和表2给出了星点提取算法原理。表1星点提取算法描述权利要求
1.基于FPGA的星敏感器图像处理方法，其特征在于包括如下步骤(1)根据星敏感器数字图像能量分布的特征、恒星成像的特征和硬件实现数学运算的特点，设计高通滤波器，由此计算星点提取阈值；(2)根据星敏感器恒星成像的特征，设计星点质心提取算法，在光电探测器逐行输出像素数据的同时，同步计算星点质心，由星点的图像特征，在运算过程中对星点与噪声点进行判别；(3)根据硬件运算的特点，采用FPGA实现基数排序算法，在完成星点提取工作后，按星点能量值对星点集合进行排序，确保能量值高的星点子集用于全天识别和姿态解算；(4)基于FPGA实现星敏感器图像处理的全部工作，并以并行运算的体系结构实现。
2.根据权利要求1所述的星敏感器图像处理方法，其特征在于所述的步骤I-(I)中， 根据星敏感器数字图像能量分布的特征、恒星成像的特征和硬件实现数学运算的特点，设计16X16pix.的高通滤波器，在光电探测器逐行给出像素数据的同时，逐个区域的完成图像滤波和阈值的计算，得到阈值矩阵。
3.根据权利要求1或2所述的星敏感器图像处理方法，其特征在于所述的步骤1-(2) 中，根据星敏感器恒星成像的特征，设计的同步提取星点质心的算法可以在逐行得到图像像素数据的同时，辨别单像素点和双像素点噪声的机制并及时舍弃。
4.根据权利要求3所述的星敏感器图像处理方法，其特征在于所述的步骤1-(3)中，根据硬件运算的特点，设计基数排序算法的状态机，并运用于星点集合排序，在排序过程中， 为减少存储空间开销，排序采用移动索引的方法进行。
5.根据权利要求4所述的星敏感器图像处理方法，其特征在于所述的步骤I-(I) (4)中，采用FPGA完成星敏感器图像处理算法的滤波、星点质心提取和排序的全部工作，并采用四级流水线方式设计并行运算体系结构。
全文摘要
本发明公开了基于FPGA的星敏感器图像处理方法，包括的步骤是图像高通滤波、星点质心计算、星点排序，基于FPGA实现星敏感器图像处理的全部工作，并以并行运算的体系结构实现。本发明采用的方法，与现有技术相比，其优点和有益效果是充分利用星敏感器恒星成像及光电探测器数字图像的特点，能在强光干扰条件下仍然可以稳定运行。充分利用FPGA在并行运算体系结构设计中的优势，采用四级流水组织图像处理算法中的运算模块，大幅提高运算速度和星敏感器的数据更新率。
文档编号G01C21/24GK102278990SQ20101020003
公开日2011年12月14日 申请日期2010年6月11日 优先权日2010年6月11日
发明者毛晓楠, 邱海辉, 郑循江 申请人:上海航天控制工程研究所