[发明专利]时空变尺度运动目标检测方法有效

专利信息
申请号: 201210594820.9 申请日: 2012-12-31
公开(公告)号: CN103077536B 公开(公告)日: 2016-01-13
发明(设计)人: 张天序;李高飞;张超;刘宽;鲁斌 申请(专利权)人: 华中科技大学
主分类号: G06T7/20 分类号: G06T7/20
代理公司: 华中科技大学专利中心 42201 代理人: 李智
地址: 430074 湖北*** 国省代码: 湖北;42
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明公开了一种时空变尺度运动目标检测方法,具体为:将原始图像转换为反差图像;根据差分强度与帧间间隔的单调递增且收敛关系,得到次优帧间间隔通过计算t0时刻两帧图像反差图像的灰度相关性得到半运动显著性图像;在半运动显著性图像空间二维四叉树分解,得到半运动显著性字块;将相邻子块合并,合并后提取运动感兴趣区域;对各感兴趣区域做时间二叉树分解,得到最优帧间间隔;在各个感兴趣区域内通过多帧累积差分提取运动目标。本发明消除背景噪音的干扰,同时也降低了光流场计算;通过指定最佳帧间间隔,移动对象可以达到最佳的运动显著状态;通过累计多帧差法能够准确地提取出运动目标。
搜索关键词: 时空 尺度 运动 目标 检测 方法
【主权项】:
时空变尺度运动目标检测方法,具体为:(1)建立当前时刻t0的帧图像的半运动显著性特征图像mf:分别计算时刻t0~t0+rt的图像的反差图像,进而分别计算图像的反差图像的半运动显著性度量值;从中找出最大半运动显著性度量值,其对应的图像与图像的帧间间隔记为次优帧间间隔通过计算图像的反差图的灰度相关性得到半运动显著性特征图像mf;所述反差图像的计算方法为:<mrow><msub><mi>c</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>&mu;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mi>&mu;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>其中,f(x,y)为原始图像,μ(x,y)为像素点(x,y)邻域内的灰度均值,经过量化后的cf(x,y)即为变换后的反差图像;(2)在半运动显著性特征图像mf内提取运动感兴趣区域:将图像mf分为四个大小相同的子块,确定各分块是否满足继续分块条件;将符合继续分块条件的子块按相同方式进一步划分子块,依此类推,直到所有的子块都不能被继续划分;所述继续分块条件为子块的长宽大于预定分块阈值,并且子块内的运动显著性度量值也大于预定显著性阈值;对于每一相邻的半显著性子块,判断其与其他相邻子块间的半运动显著性度量值的相似性,当相似性大于预设相似性阈值时将两者合并;假设当前子块的半运动显著性度量值为α,其周围领域内一个子块的半运动显著性度量值为β,这两个子块之间半运动显著性度量值的相似性计算方法如下:<mrow><mi>s</mi><mi>i</mi><mi>m</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>&beta;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>m</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>&beta;</mi><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>&beta;</mi><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>在合并后得到的M个子块内提取连通的子块得到M个运动感兴趣区域Ri,i=1,2,...,M;(3)计算各运动感兴趣区域的最优帧间间隔:计算图像的反差图像内感兴趣区域Ri的半运动显著性度量值;从中找出最大半运动显著性度量值,其对应的图像与图像的帧间间隔记为最优帧间间隔(4)在感兴趣区域内通过多帧累积差分增强感兴趣区域:增强感兴趣区域<mrow><msubsup><mi>d</mi><mi>t</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mrow><msubsup><mi>d</mi><mrow><mi>t</mi><mo>-</mo><mi>&Delta;</mi><msub><mover><mi>t</mi><mo>~</mo></mover><mi>i</mi></msub><mo>,</mo><mi>t</mi></mrow><mi>i</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&times;</mo><msubsup><mi>d</mi><mrow><mi>t</mi><mo>,</mo><mi>t</mi><mo>+</mo><mi>&Delta;</mi><msub><mover><mi>t</mi><mo>~</mo></mover><mi>i</mi></msub></mrow><mi>i</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></msqrt><mo>,</mo></mrow>其中,<mrow><msubsup><mi>d</mi><mrow><mi>t</mi><mo>-</mo><mi>&Delta;</mi><msub><mover><mi>t</mi><mo>~</mo></mover><mi>i</mi></msub><mo>,</mo><mi>t</mi></mrow><mi>i</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>&Delta;</mi><msub><mover><mi>t</mi><mo>~</mo></mover><mi>i</mi></msub></mrow></munderover><mo>|</mo><msub><mi>f</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>f</mi><mrow><mi>t</mi><mo>-</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>|</mo><mo>/</mo><mi>&Delta;</mi><msub><mover><mi>t</mi><mo>~</mo></mover><mi>i</mi></msub><mo>,</mo></mrow><mrow><msubsup><mi>d</mi><mrow><mi>t</mi><mo>,</mo><mi>t</mi><mo>+</mo><mi>&Delta;</mi><msub><mover><mi>t</mi><mo>~</mo></mover><mi>i</mi></msub></mrow><mi>i</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>&Delta;</mi><msub><mover><mi>t</mi><mo>~</mo></mover><mi>i</mi></msub></mrow></munderover><mo>|</mo><msub><mi>f</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>f</mi><mrow><mi>t</mi><mo>+</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>|</mo><mo>/</mo><mi>&Delta;</mi><msub><mover><mi>t</mi><mo>~</mo></mover><mi>i</mi></msub><mo>,</mo></mrow>(x,y)∈Ri;(5)在各增强感兴趣区域内进行OTSU阈值分割和形态学处理,得到多个运动目标:设检查到的目标数量为N,每个目标用Oj(j=1,2,...N)表示,同时得到包围目标Oj的最小外接矩形的长和宽Oj相对于其局部背景的反差cj,Oj所占区域Ωi的形心位置坐标计算每个运动目标的7维特征向量;对于某个目标Oj其中为包围目标Oj的最小外接矩形的长和宽;为Oj运动的表观速度;cj为Oj相对于其局部背景的反差;为Oj所占区域Ωi的形心位置坐标;cj在步骤(5)中目标提取的过程中计算出来;Oj的运动表观速度可以如下计算:<mrow><msub><mi>u</mi><msub><mi>o</mi><mi>j</mi></msub></msub><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mrow><msub><mi>t</mi><mn>0</mn></msub><mo>+</mo><mi>&Delta;</mi><mover><mi>t</mi><mo>~</mo></mover></mrow></msub><mo>-</mo><msub><mi>x</mi><msub><mi>t</mi><mn>0</mn></msub></msub></mrow><mrow><mi>&Delta;</mi><msub><mover><mi>t</mi><mo>~</mo></mover><mi>j</mi></msub></mrow></mfrac><mo>,</mo></mrow><mrow><msub><mi>v</mi><msub><mi>o</mi><mi>j</mi></msub></msub><mo>=</mo><mfrac><mrow><msub><mi>y</mi><mrow><msub><mi>t</mi><mn>0</mn></msub><mo>+</mo><mi>&Delta;</mi><mover><mi>t</mi><mo>~</mo></mover></mrow></msub><mo>-</mo><msub><mi>y</mi><msub><mi>t</mi><mn>0</mn></msub></msub></mrow><mrow><mi>&Delta;</mi><msub><mover><mi>t</mi><mo>~</mo></mover><mi>j</mi></msub></mrow></mfrac><mo>,</mo></mrow>其中为目标Oj在t0时刻的形心位置,为Oj时刻的形心位置,是t0时刻目标的最优帧间间隔,目标的速度等于时刻t0到时刻的平均速度;所述半运动显著性度量值的计算方法为:令两幅图像的反差图像表示为的半运动显著性度量值的计算公式为:<mrow><msub><mi>value</mi><mi>R</mi></msub><mrow><mo>(</mo><msub><mi>c</mi><msub><mi>f</mi><mi>s</mi></msub></msub><mo>,</mo><msub><mi>c</mi><msub><mi>f</mi><mi>e</mi></msub></msub><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mrow><msub><mi>&Sigma;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo><mo>&Element;</mo><mi>R</mi></mrow></msub><msub><mi>c</mi><msub><mi>f</mi><mi>s</mi></msub></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&times;</mo><msub><mi>c</mi><msub><mi>f</mi><mi>e</mi></msub></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><msqrt><mrow><msub><mi>&Sigma;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo><mo>&Element;</mo><mi>R</mi></mrow></msub><msup><mrow><mo>&lsqb;</mo><msub><mi>c</mi><msub><mi>f</mi><mi>s</mi></msub></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow></msqrt><msqrt><mrow><msub><mi>&Sigma;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo><mo>&Element;</mo><mi>R</mi></mrow></msub><msup><mrow><mo>&lsqb;</mo><msub><mi>c</mi><msub><mi>f</mi><mi>e</mi></msub></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mfrac><mo>,</mo></mrow>R表示图像的所有像素点集。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于华中科技大学,未经华中科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201210594820.9/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top