《空间微弱运动目标检测方法.pdf》由会员分享，可在线阅读，更多相关《空间微弱运动目标检测方法.pdf（7页完整版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102081800 A (43)申请公布日 2011.06.01 CN 102081800 A *CN102081800A* (21)申请号 201110005155.0 (22)申请日 2011.01.06 G06T 7/20(2006.01) G06T 7/60(2006.01) (71)申请人 西北工业大学 地址 710072 陕西省西安市友谊西路 127 号 (72)发明人 张艳宁 段锋 姚睿 孙瑾秋 施建宇 杨涛 朱宇 张永鹏 张臻 李磊 (74)专利代理机构 西北工业大学专利中心 61204 代理人 黄毅新 (54) 发明名称 空间微弱运动目标检测方法 (。
2、57) 摘要 本发明公开了一种空间微弱运动目标检测方 法， 用于解决现有的空间微弱运动目标检测方法 检测率低的技术问题。技术方案是利用目标的灰 度、 背景的结构性等特性与目标轨迹连续性相结 合的方法进行运动目标的检测， 图像预处理阶段 采用了基于灰度与基于区域的分割， 有效地去除 了噪声的影响， 在检测阶段利用两帧间目标运动 连续性特征， 有效剔除了恒星干扰， 避免了使用 多帧累加剔除恒星造成的处理效率的降低， 以及 背景光照、 噪声影响下不能去除恒星背景的问题， 并且对整个目标轨迹进行后处理， 剔除虚假目标 ； 并对丢失目标进行填补， 在信噪比小 3 时， 检测率 由现有技术的不足 90提高。
3、到 95以上。 (51)Int.Cl. (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 2 页 说明书 4 页 CN 102081805 A1/2 页 2 1. 一种空间微弱运动目标检测方法， 其特征在于包括下述步骤 ： (a) 将空间微弱运动目标及恒星成像后的图像分割出星点， 把图像分为 SS 区域， 在 每个SS区域内设定自适应阈值+3， 并用阈值+3对图像进行二值化处理， 根据二 值分割结果， 建立无向图 G(V， E)， 将目标连通性的问题转化为求无向图 G(V， E) 的强连通森 林问题， 使用广度优先搜索策略提取连通分量做为后继处理的星点集合 C(t) ； 。
4、式中， 是 区域背景灰度均值， 是区域背景灰度标准差， V 是所有像素点， E 是结点之间的边， t 表示 第 t 帧图像， Ci(t)， i 1， 2， K 表示星点集合中的第 i 个星点， K 表示集合 C(t) 中共 有 K 个星点 ； 获得星点集合之后， 提取星点面积 N(Ck(t) 并通过 (1) 式求得星点质心 zk(t) (xc yc) 式中， (xc， yc) 表示星点 Ck(t) 质心在 x， y 方向的坐标， (m， n) 星点区域的长宽， I(x， y， t) 是第 t 帧 (x， y) 坐标处的灰度值 ； (b) 定义两帧图像间两个星点的距离 ： 式中， zi(t)x表。
5、示第 t 帧星点 Ci(t) 质心在 x 方向的坐标， zi(t)y表示第 t 帧星点 Ci(t) 质心在 y 方向的坐标 ； 对相邻两帧图像的星点集合 C(t)， C(t+1) 中的每个星点进行配对， 并且两个星点的区 域大小之差满足条件 ： N(Ci(t)-N(Cj(t+1) N(Ci(t) (3) 通过计算配对星点之间的距离， 区分出恒星与非恒星点 ； 第 t， t+1， t+2 帧图像各得到 一个候选目标集合， 其中前两帧得到的候选配对星点的集合为 TPt(Ci(t)， Cj(t+1)、 后两帧 得到的候选配对星点的集合为 TPt+1(Cj(t+1)， Ck(t+2)， 两个集合中的 。
6、Cj(t+1) 是它们公共 的元素， 通过两个配对星点的公共元素建立目标轨迹 ； (c) 遍历目标轨迹列表中的每个目标轨迹 TC(i)， 式中， i 1， 2， T， T 是检测到 的目标轨迹总数， 对其处理 ； 若目标轨迹中目标出现的帧数小于设定阈值 ， 将其删除 ； 从目标轨迹列表中取出 TC(j)， 式中， j 1， 2， T 且 j i， 若 TC(j) 的第一个目 标出现的帧序号Kpre与TC(i)的最后一个目标出现的帧序号Knext不满足条件Kpre-Knext2， 两条目标链不合并， 取出另一条目标轨迹进行对比 ； 若满足条件， 继续向下执行 ； 计算两条目标链轨迹之间的夹角， 。
7、即 TC(i) 的后两个目标点构成向量 a， TC(j) 的前 两个目标点构成向量 b， 若 a 与 b 之间的夹角大于 30， 不合并目标轨迹， 转步骤 ； 否则， 继续向下执行 ； 计算轨迹TC(i)的最后一个目标与轨迹TC(j)的第一个目标之间的距离DistChaini， j， 并计算 TC(i) 的目标平均运动速度 AvgDisti， 如果 权 利 要 求 书 CN 102081800 A CN 102081805 A2/2 页 3 式中， 等于 2 时， 不合并目标轨迹， 转步骤 ； 否则， 继续向下执行 ； 如果上述条件均满足， 合并 TC(i) 和 TC(j)， 直至所有的轨迹遍。
8、历完成 ； 在目标轨迹合并之后， 对由于恒星遮挡造成的丢失目标进行填补。 2. 一种权利要求 1 所述的空间微弱运动目标检测方法， 其特征在于 ： 所述 SS 区域的 大小是 100100。 权 利 要 求 书 CN 102081800 A CN 102081805 A1/4 页 4 空间微弱运动目标检测方法 技术领域 0001 本发明涉及一种空间目标检测方法， 特别是一种空间微弱运动目标检测方法。 背景技术 0002 文献 “深空背景弱小运动目标检测算法研究， 光学技术， 2010， Vol36(2)， p209-212” 公开了一种基于序列图像多帧累加并提取特征的深空背景下弱小运动目标的检。
9、 测算法。该方法先进行背景抑制 ； 采用序列图像多帧累加增强目标与恒星的对比度， 提取 局部星图并剔除恒星干扰 ； 最后结合候选目标的特征完成目标检测。由于在剔除恒星时需 要使用多于 5 帧图像进行累加， 降低了算法处理的效率， 并且在背景光照、 噪声有剧烈变化 时， 多帧累加并不能增强目标与恒星的对比度 ； 在获取候选目标之后， 只是了采用基于逻辑 的最近邻关联方法在相邻两帧间对目标进行关联， 并未完全利用序列图像整个目标轨迹信 息， 因此不能有效的剔除虚假目标， 在目标局部信噪比小于 3 时， 检测率不足 90， 且不能 解决目标被恒星遮挡问题。 发明内容 0003 为了克服现有的空间微弱。
10、运动目标检测方法检测率低的不足， 本发明提供一种空 间微弱运动目标检测方法。该方法依据目标灰度及背景的结构性特征用 2 帧图像去除恒 星， 提取目标， 并利用整个目标运动轨迹剔除虚假目标， 实现真实微弱运动目标的检测， 可 以提高空间微弱运动目标检测方法的检测率。 0004 本发明解决其技术问题所采用的技术方案 ： 一种空间微弱运动目标检测方法， 其 特点是包括下述步骤 ： 0005 (a) 将空间微弱运动目标及恒星成像后的图像分割出星点， 把图像分为 SS 区 域， 在每个 SS 区域内设定自适应阈值 +3， 并用阈值 +3 对图像进行二值化处理， 根据二值分割结果， 建立无向图 G(V， 。
11、E)， 将目标连通性的问题转化为求无向图 G(V， E) 的 强连通森林问题， 使用广度优先搜索策略提取连通分量做为后继处理的星点集合 C(t) ； 式 中， 是区域背景灰度均值， 是区域背景灰度标准差， V 是所有像素点， E 是结点之间的 边， t 表示第 t 帧图像， Ci(t)， i 1， 2， K 表示星点集合中的第 i 个星点， K 表示集合 C(t) 中共有 K 个星点 ； 0006 获得星点集合之后， 提取星点面积N(Ck(t)并通过(1)式求得星点质心zk(t) (xc yc) 0007 0008 式中， (xc， yc) 表示星点 Ck(t) 质心在 x， y 方向的坐标，。
12、 (m， n) 星点区域的长宽， I(x， y， t) 是第 t 帧 (x， y) 坐标处的灰度值 ； 说 明 书 CN 102081800 A CN 102081805 A2/4 页 5 0009 (b) 定义两帧图像间两个星点的距离 ： 0010 0011 式中， zi(t)x表示第 t 帧星点 Ci(t) 质心在 x 方向的坐标， zi(t)y表示第 t 帧 星点 Ci(t) 质心在 y 方向的坐标 ； 0012 对相邻两帧图像的星点集合 C(t)， C(t+1) 中的每个星点进行配对， 并且两个星点 的区域大小之差满足条件 ： 0013 N(Ci(t)-N(Cj(t+1) N(Ci(t。
13、) (3) 0014 通过计算配对星点之间的距离， 区分出恒星与非恒星点 ； 第 t， t+1， t+2 帧图像各 得到一个候选目标集合， 其中前两帧得到的候选配对星点的集合为 TPt(Ci(t)， Cj(t+1)、 后 两帧得到的候选配对星点的集合为 TPt+1(Cj(t+1)， Ck(t+2)， 两个集合中的 Cj(t+1) 是它们 公共的元素， 通过两个配对星点的公共元素建立目标轨迹 ； 0015 (c) 遍历目标轨迹列表中的每个目标轨迹 TC(i)， 式中， i 1， 2， T， T 是检 测到的目标轨迹总数， 对其处理 ； 0016 若目标轨迹中目标出现的帧数小于设定阈值 ， 将其删。
14、除 ； 0017 从目标轨迹列表中取出 TC(j)， 式中， j 1， 2， T 且 j i， 若 TC(j) 的第一 个目标出现的帧序号 Kpre与 TC(i) 的最后一个目标出现的帧序号 Knext不满足条件 Kpre-Knext 2， 两条目标链不合并， 取出另一条目标轨迹进行对比 ； 若满足条件， 继续向下执行 ； 0018 计算两条目标链轨迹之间的夹角， 即 TC(i) 的后两个目标点构成向量 a， TC(j) 的前两个目标点构成向量b， 若a与b之间的夹角大于30， 不合并目标轨迹， 转步骤 ； 否 则， 继续向下执行 ； 0019 计算轨迹 TC(i) 的最后一个目标与轨迹 TC。
15、(j) 的第一个目标之间的距离 DistChaini， j， 并计算 TC(i) 的目标平均运动速度 AvgDisti， 如果 0020 0021 式中， 等于 2 时， 不合并目标轨迹， 转步骤 ； 否则， 继续向下执行 ； 0022 如果上述条件均满足， 合并 TC(i) 和 TC(j)， 直至所有的轨迹遍历完成 ； 0023 在目标轨迹合并之后， 对由于恒星遮挡造成的丢失目标进行填补。 0024 所述 SS 区域的大小是 100100。 0025 本发明的有益效果是 ： 由于利用目标的灰度、 背景的结构性等特性与目标轨迹连 续性相结合的方法进行运动目标的检测， 图像预处理阶段采用了基于灰。
16、度与基于区域的分 割， 有效地去除了噪声的影响， 在检测阶段利用两帧间目标运动连续性特征， 有效剔除了恒 星干扰， 避免了使用多帧累加剔除恒星造成的处理效率的降低， 以及背景光照、 噪声影响下 不能去除恒星背景的问题， 并且对整个目标轨迹进行后处理， 剔除虚假目标 ； 并对丢失目标 进行填补， 在信噪比小 3 时， 检测率由现有技术的不足 90提高到 95以上。 0026 下面结合具体实施方式对本发明作详细说明。 具体实施方式 0027 1、 单帧图像预处理。 说 明 书 CN 102081800 A CN 102081805 A3/4 页 6 0028 目标及恒星成像后在图像中都表现为星点，。
17、 首先分割出星点， 把图像分为 SS 区 域， S 是区域大小， 设定为 100， 对每个区域统计直方图， 找出直方图最大值时的灰度作为区 域的背景均值 ， 然后迭代求出背景噪声的 3 倍方差 3， 以 +3 做为分割阈值， 按此 阈值对小块区域进行 0 和 1 的二值化处理。根据二值分割的结果建立一个无向图 G(V， E)， 其中， V 为所有像素点， E 为结点之间的边， 在该无向图中， 如果灰度阈值分割得到的二值图 中， 如果有两个像素点 (x1， y1) 和 (x2， y2) 相邻， 且两个像素点在该二值图上均为 1， 那么在 建立的无向图中则存在一条边 E(x1， y1)， (x2，。
18、 y2)， 其权值为 1。在该无向图建立之后， 将 目标连通性的问题转化为求无向图 G(V， E) 的强连通森林问题， 使用广度优先搜索策略求 出图的强联通森林作为后继处理的星点集合 C(t)， t 表示第 t 帧图像， Ci(t)， i 1， 2， K 表示星点集合中的第 i 个星点， K 表示集合 C(t) 中共有 K 个星点。 0029 获得星点集合之后， 提取星点特征， 做为后继处理分离目标与干扰的依据。 0030 (1) 星点面积 ： 根据分割出的亮斑的空间连通性搜索出星点区域为 (Ck(t)， 用 N(Ck(t) 表示 Ck(t) 内像素的个数， 即星点面积。 0031 (2) 星。
19、点质心位置 ： 星点图像的质心通过图像灰度函数的平衡点来计算， 星点的 质心 zk(t) (xc yc)( 即星点类别 Ck(t) 的聚类中心 ) 通过下式求得 ： 0032 0033 式中， (xc， yc) 表示星点 Ck(t) 质心在 x， y 方向的坐标， (m， n) 星点区域的长宽， I(x， y， t) 是第 t 帧 (x， y) 坐标处的灰度值。 0034 2、 候选目标检测与轨迹关联。 0035 对两帧图像中的星点的集合 C(t)， C(t+1) 做帧间星点的匹配。定义两帧间两个星 点的距离为它们质心位置的欧氏距离， 即 ： 0036 0037 式中， zi(t)x表示第 t。
20、 帧星点 Ci(t) 质心在 x 方向的坐标， zi(t)y表示第 t 帧 星点 Ci(t) 质心在 y 方向的坐标。 0038 首先前一帧中的星点 Ci(t) 在后帧的对应质心位置寻找候选星点， 并计算与候选 星点之间的欧氏距离， 找到距离最小的候选星点， 并且如果两个星点的区域大小之差满足 以下条件 ： 0039 N(Ci(t)-N(Cj(t+1) N(Ci(t) (7) 0040 则认为两个星点配对成功， 否则跳过该星点， 与其他星点继续判定 ； 按照这个步骤 依次进行， 直至两帧内的所有星体全部配对完成。 星点配对完成之后， 得到的是一个包含配 对点之间距离的一维向量， 正常情况下， 。
21、相邻两帧图像对应位置的恒星之间的距离为 0 或 一个小于 1 的数， 目标具有与恒星不同的运动方向和速度， 它们的运动距离通常会较大。根 据序列图像目标运动的特性， 设定阈值为 1， 过滤掉两帧中运动距离大于阈值的星点， 经过 过滤后剩余的星点即为候选目标。式中， 是系数， 取值为 0.5。 0041 通过连续三帧图像中前两帧和后两帧目标对的公共元素建立目标轨迹。具体如 说 明 书 CN 102081800 A CN 102081805 A4/4 页 7 下 ： 0042 (1)设当前帧为第t+2帧， 得到的候选目标对集合为TPt+1(Cj(t+1)， Ck(t+2)， 其中 Cj(t+1) 。
22、为第 t+1 帧图像中星点集合的第 j 个星点， Ck(t+2) 为第 t+2 帧图像中星点集合的 第 k 个星点， 对集合中的每个星点 Cj(t+1) 与已有目标轨迹列表 TC(i) 中的每个轨迹 TC(i) 的最后一个目标星点做匹配， 如果都匹配不上， 认为是新目标， 对TPt+1(Cj(t+1)， Ck(t+2)建 立新的目标轨迹 ； 如果有匹配的轨迹， 转入 (2) 继续执行 ； 0043 (2) 判断 Cj(t+1)， Ck(t+2) 之间的距离是否与匹配的目标轨迹中目标的平均运动 位移 ( 即目标轨迹中相邻两个目标运动距离之和的均值 ) 一致， 并且判断 Cj(t+1)， Ck(t。
23、+2) 形成的向量与轨迹 TC(i) 之间的夹角是否小于 30， 如果满足条件， 把 TPt+1(Cj(t+1)， Ck(t+2) 加入轨迹 TC(i)， 否则建立新链。 0044 3、 目标轨迹后处理。 0045 在图像序列中的所有帧进行检测处理完成之后， 对得到的目标轨迹进行后处理， 对由于恒星遮挡造成的丢失目标进行填补。算法如下 ： 0046 (1) 遍历目标轨迹列表中的每个目标轨迹 TC(i)， 其中 i 1， 2， ， T， T 为检测到 的目标轨迹总数， 对其处理 ； 0047 (2) 若目标轨迹中目标出现的帧数小于设定阈值 3， 则认为是由于噪声产生的链， 将其删除 ； 0048。
24、 (3) 从目标轨迹列表中取出 TC(j)， j 1， 2， ， T 且 j i， 若 TC(j) 的第一个目 标出现的帧序号Kpre与TC(i)的最后一个目标出现的帧序号Knext不满足条件Kpre-Knext2， 认为两条目标链不可合并， 取出另一条目标轨迹进行对比 ； 若满足条件， 继续向下执行 ； 0049 (4)计算两条目标链轨迹之间的夹角， 即TC(i)的后两个目标点构成向量a， TC(j) 的前两个目标点构成向量 b， 若 a 与 b 之间的夹角大于 30， 则不合并， 转 (3) ； 否则， 继续 向下执行 ； 0050 (5) 计算轨迹 TC(i) 的最后一个目标与轨迹 TC。
25、(j) 的第一个目标之间的距离 DistChaini， j， 并计算 TC(i) 的目标平均运动速度 AvgDisti， 如果 0051 0052 式中， 取 2， 则不合并目标轨迹， 转 (3) ： 否则， 继续向下执行 ； 0053 (6) 如果上述条件均满足， 合并 TC(i) 和 TC(j)， 直至所有的轨迹遍历完成。 0054 在目标轨迹合并之后， 对丢失的目标进行填充， 填补过程如下 ： 0055 (1) 遍历目标轨迹列表中的所有目标轨迹 ； 0056 (2) 找到目标轨迹中间断的帧， 假设为第 k 帧， 根据目标轨迹中目标的平均运动距 离预测目标在中断帧中的位置， 并根据预测的位置设置搜索框到该帧的星点集合 C(k) 中 寻找目标星点， 如果找到， 则认为找到的目标即为丢失的目标， 定位目标插入位置， 加入到 目标轨迹中， 否则把预测的目标加入到目标轨迹中 ； 直到所有的目标轨迹处理完成。 说 明 书 CN 102081800 A 。