[发明专利]空间碎片星点提取与定位方法

摘要

本发明公开了一种空间碎片星点提取与定位方法，用于解决现有方法对相邻空间碎片目标分割困难的技术问题。技术方案是通过MSER算法在候选空间碎片区域内过构建数据结构树，提取出该区域内空间碎片像素集合并确定空间碎片星点个数。在MSER算法得到的像素点树形结构中，通过检测该候选空间碎片区域中最大的MSER，为后续拟合GMM模型生成训练样本。在基于GMM的定位过程中，通过EM算法迭代求解混合模型中高斯成分的均值与协方差，对应的均值就是各个空间碎片星点光斑的质心。该方法结合MSER和GMM模型的特点，通过GMM模型对近距离弥散光斑质心进行共同求解，解决了相邻空间碎片目标分割困难的技术问题。

主权项

一种空间碎片星点提取与定位方法，其特征在于包括以下步骤：步骤一、对于星空图像I，首先对针对星空图像中每个像素统计图像的灰度直方图，根据像素值按升序对图像中的像素点进行排序，去除像素灰度值最小的λ1％的像素点，以减弱暗噪声对背景的影响，去除像素灰度值最高的λ2％的像素点，以减弱强噪声和恒星对背景估计的影响，其中λ1,λ2为设定的背景阈值参数；利用剩余的(100‑λ1‑λ2)％像素点，求得背景的灰度均值μ和灰度标准差σ，使用阈值μ+κσ对图像进行二值化处理，κ为需要设置的分割系数；二值化处理公式为：式中，I(x,y)表示在坐标为(x,y)处的图像I中像素的灰度值，B(x,y)就是(x,y)处像素点的二值化结果，把二值化的后的图像记为B；这一步为星图的二值化粗分割，因此每个B(x,y)＝1的连通区域中可能包含多于一个的星点；步骤二、对空间碎片区域进行粗分割；根据上一步图像二值化的结果，当B(x,y)＝1时表示该像素为空间碎片星点区域，当B(x,y)＝0时表示该像素为非有效空间碎片星点区域，即背景区域；对B中B(x,y)＝1的点进行广度优先搜索，提取到所有满足B(x,y)＝1的点的连通域，所有连通域集合为其中一个连通域表示为Ω；Ω为一个空间碎片星点区域，包含了一系列二值化结果为1的像素点坐标，并且他们相互连通，而且对于任一个Ω中的坐标为(x,y)的点，它在原始图像中的灰度值为I(x,y)；步骤三、在提取空间碎片区域后，需要对每一个连通域Ω通过MSER算法进行局部极大灰度值和最大稳定极值区域检测；MSER检测是通过不同阈值对空间碎片区域进行二值化处理，根据不同的结果对稳定极大灰度值点和区域进行检测；对于一个星点区域Ω，当令二值化阈值为t时，公式(2)中的Ω>t表示在空间碎片区域Ω中，像素灰度值大于t的像素点集合，当二值化阈值t不断变大的时候集合Ω>t的元素个数会越来越少，同理t变小的时候集合Ω>t的元素个数会变多；在MSER检测过程中通过集合中元素个数的变化速率来表示一个区域Ω>t是否是一个稳定的空间碎片弥散光斑；公式(3)表示当分割阈值变化大小为d时，光斑变化的稳定度计算方法；这里|Ω>t|，|Ω>t+d|和|Ω>t|分布表示集合Ω>t，Ω>t+d和Ω>t的元素个数，d为需要设置参数，st 为光斑的稳定度；Ω>t＝{(x,y)|I(x,y)>t,(x,y)∈Ω}         (2) 利用MSER算法在一个空间碎片星点连通区域Ω中分割提取星点的步骤为：（1）.首先对于图像中某个空间碎片星点连通区域Ω内的每个像素点灰度值按照大小排序，依次初始化他们为未访问状态再根据各个像素点的灰度值大小和他们空间位置关系，按照步骤（2）、步骤（3）的描述把他们组织成一颗树形结构，这个树的各个节点都代表了图像中的一个像素点；（2）.按照像素灰度值从大到小遍历Ω中所有的像素，针对某个(x,y)处的像素点遍历它周围相邻的8个点；假设(x′,y′)为一个与(x,y)相邻的位置，那么如果(x′,y′)处像素点未被遍历，则直接遍历下一个相邻点，如果已经被遍历那么：a)当Root(x,y)＝Root(x′,y′)相同时，则不进行任何处理，遍历(x,y)处像素的下一个相邻点，因为(x,y)处像素点与(x′,y′)处像素点已经处于同一个子树中；这里的Root(x,y)表示(x,y)处像素点在当前树结构中对应的根像素点的位置；b)当I(Root(x,y))＝I(Root(x′,y′))时，则把Root(x′,y′)处像素点作为子节点添加到Root(x,y)处像素点下；c)当I(Root(x,y))<I(Root(x′,y′))时，则把Root(x′,y′)处像素点添加为Root(x,y)处像素点的子节点；（3）.标记当前(x,y)处像素点为已被遍历像素点，并重新回到步骤（2）直到完成对所有像素点的遍历；（4）.当所有的像素点按照上面步骤遍历完毕，依照像素点灰度值大小关系的树形结构建立完成，该树中叶节点个数就等于该空间碎片区域Ω中局部极大灰度值像素点个数，也对应了区域内星点光斑数目，这个数目记为K，极值点的个数、位置将用于以下基于GMM模型和EM迭代算法的星点定位的初始化；为了使用EM算法拟合GMM模型提供样本，需要在空间碎片区域内检测最大稳定极值区域；对于空间碎片光斑区域Ω内，每当取到一个二值化分割阈值t，总是对空间碎片区域Ω进行分割得到一个子光斑为Ω>t，当该光斑的稳定度st小于tmser时，即满 足公式(4)时，该子光斑就是为Ω的一个MSER，这里tmser为判断一个子光斑是否为MSER的阈值；当调整t，就得到一系列MSER，这一系列的MSER呈现树形的包含现象，即t较小时分割得到的MSER包含t较大时分割得到的MSER，而最大的MSER包含所有其他的MSER；对区域Ω提取不同t情况下的MSER之后得到MSER集合其中包含像素点最多的MSER记为C；步骤四、考虑单个空间碎片弥散光斑可以被合理假设为一个高斯函数形状的的点状光斑，对于包含多个弥散光斑的空间碎片区域，利用混合高斯模型对空间碎片区域内的多个弥散光斑进行建模；GMM在质心定位的时候需要多个像素点样本；C为空间碎片区域Ω的MSER集合中包含最多像素点的MSER，对其范围内任一位置(x,y)处的像素点进行重采样，得到个样本点，这个样本点的集合记为S(x,y)，为向上取整符号；那么C中所有的样本点集合记为SC，其中n＝1,...,L表示像素点样本的序号，L表示样本点个数，即集合SC的中的成员数目，mn表示第n个样本像素点的位置；获得了样本点集合SC和星点光斑数目K，利用GMM模型和相应的EM算法对空间碎片光斑区域Ω内K个光斑同时进行定位；具体的模型描述和定位算法如下：（1）.定义用于拟合空间碎片区域星点光斑的混合高斯模型；由于空间碎片区域是有多个弥散光斑组成，首先将对于空间碎片区域中的任一位置点变量(x,y)记为随机变量m＝(x,y)，它服从公式(6)描述的高斯混合分布，其中πk、μk和Σk分别为第k个光斑的权重系数、均值和协方差矩阵，k＝1,2,...,K为GMM中高斯成分的序号，N(m|μk,Σk)表示μk和Σk为参数的随机变量m的高斯分布，高斯分布成分数即为区域内星体光斑数总数为K；利用EM算法根据样本点求解得到的K个高斯成分的参数πk、μk和Σk，而μk为对应的第k个星点光斑的质心；由于EM算法为迭代优化算法，所以需要先根据以下步骤（2）、步骤（3）对各参数进行初始化；（2）.初始化每个弥散光斑的中心；每个弥散光斑的质心，即K个高斯成分的均值μk,k＝1,2,......,K，被初始化化为MSER检测中得到的极大灰度值像素的中心位置；（3）.根据每个样本点到各个弥散光斑的初始化中心的欧式距离，为每个样本点选取最近中心点为类别中心，把样本点分为K类，并且把属于第k类的样本集合记为SCk；样本点集合SC被分割为K个集合，SCk中像素点样本被记为上标表示其属于SCk集合，下标为SCk集合元素序号；通过公式(7)初始化所有的πk，其中运算符|·|表示求集合的元素个数；根据公式(8)对K个高斯分布成分的协方差矩阵进行初始化（4）.在完成初始化工作后，进行质心估计，高斯混合模型的质心使用EM算法来估计，EM算法为迭代算法，每次迭代步骤分为两步，第一步为期望估计，第二步为最大化；第一次迭代为初始化，对于第α次迭代，α≥2，具体如下：a)通过初始化的参数计算第n个样本点mn由第个k个高斯模型生成的概率估计更新参数，包括方差、均值和权值系数；具体通过公式(10)计算新的均值，通过公式(11)计算新的协方差方差矩阵，通过公式(12)得到新的权重系数；其中计算方法如公式(13)所示，L为训练样本的数目；迭代执行步骤b)，对样本点集合SC，计算在不同参数下的对数似然函数值，当迭代到满足不等式(14)时，迭代结束；tgmm是判断两次对数似然值是否收敛的阈值；其中与分布表示模型参数与的似然值；其中的似然值计算方式如公式(15)所示对的计算同理；迭代停止后得到的K个高斯成分的均值μk,k＝1,...,K为空间碎片区域Ω中的K个空间碎片星点光斑的亚像素坐标定位结果；依次遍历联通区域集合中所有空间碎片区域，完成对图像I中所有空间碎片星点目标的提取与定位工作。