空间碎片星点提取与定位方法技术

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

【技术实现步骤摘要】
空间碎片星点提取与定位方法
本专利技术涉及一种空间碎片星点提取与定位方法。
技术介绍
文献“基于能量累加的空间目标星像质心定位，《光学精密工程》，2011，第12期，p2043-3048”提出了基于能量累积的插值求解质心方法。该方法使用插值方式降低离散采样点和实际感光区域光线照度的不一致性；利用弥散光斑能量累加自适应地确定空间碎片目标的有效窗口。由于这种自适应确定空间碎片目标有效窗口的方法仅适用于单个独立的空间碎片进行定位，当多个空间碎片由于距离较近，而使得图像上的星点光斑相互黏连或重叠时就会难以分割提取和定位。
技术实现思路
为了克服现有方法对相邻空间碎片目标分割困难的不足，本专利技术提供一种空间碎片星点提取与定位方法。该方法通过MSER算法在候选空间碎片区域内，通过构建数据结构树，提取出该区域内空间碎片像素集合并确定空间碎片星点个数，为利用高斯混合模型对星点光斑建模和利用最大期望EM算法求解是的初始化提供了基础。在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)对空间碎片区域内的多个弥散光斑进行建模。GMM在质心定位的时候需要多个像素点样本。C为空间碎片区域Ω的MSER集合中包含最多像素点的MSER，对其范围内任一位置(x,y)处的像素点进行重采样，得到个样本点，这个样本点的集合记为S(x,y)，为向上取整符号。那么C中所有的样本点集合记为SC，其中n＝1,...,L表示像素点样本的序号，L表示样本点个数，即集合SC的中的成员数目，mn表示第n个样本像素点的位置。获得了样本点集合SC和星点光斑数目K，利用GMM模型和相应的EM算法对空间碎片光斑本文档来自技高网...

【技术保护点】
一种空间碎片星点提取与定位方法，其特征在于包括以下步骤：步骤一、对于星空图像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较小时分割得到的...

【技术特征摘要】
1.一种空间碎片星点提取与定位方法，其特征在于包括以下步骤：步骤一、对于星空图像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；步骤四、考虑单个空间碎片弥散光斑能够被合理假设为一个高斯函数形状的的点...

【专利技术属性】
技术研发人员：张艳宁，巩东，丁王斌，李海森，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61