发明名称 基于恒边长梯度加权图切的优化雷达微弱目标检测方法
摘要 本发明公开了一种基于恒边长梯度加权图切算法的优化雷达弱目标检测方法。该方法具体思路如下:首先对相干积累后的距离-多普勒序列进行短时非相关积累;然后利用目标模型和杂波模型在混合高斯框架下,对距离-多普勒视频帧进行基于恒边长梯度加权图切算法的全局优化目标检测,得到目标视频序列;判断该视频帧是否为最后一帧,若是结束检测,反之,利用kalman预测方法对目标视频序列进行目标特征检测与跟踪完成目标模型更新,利用混合高斯更新模型算法实现背景模型的更新;最后读取新捕获到的下一帧视频信息,重复上述步骤直至所有捕获到的视频信息均被检测完为止。本发明优化和改善了对于雷达弱目标的检测结果。
申请公布号 CN103383451A 申请公布日期 2013.11.06
申请号 CN201310225376.8 申请日期 2013.06.07
申请人 杭州电子科技大学 发明人 郭春生;刘丹;董秀青;郭云飞
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 1.基于恒边长梯度加权图切的优化雷达微弱目标检测方法,其特征在于该方法的具体步骤是:步骤1.模型参数初始化:设置背景模型的均值<img file="2013102253768100001DEST_PATH_IMAGE002.GIF" wi="21" he="25" />及方差<img file="2013102253768100001DEST_PATH_IMAGE004.GIF" wi="22" he="26" />;预测模型的观测噪声<img file="DEST_PATH_IMAGE006.GIF" wi="18" he="25" />及系统噪声<img file="DEST_PATH_IMAGE008.GIF" wi="20" he="25" />;步骤2.通过相干积累及非相干积累捕捉包含雷达弱目标的视频流数据信息,并对获得的数据信息进行二值化处理,得到大小为<img file="2013102253768100001DEST_PATH_IMAGE010.GIF" wi="48" he="20" />的像素点矩阵,<img file="2013102253768100001DEST_PATH_IMAGE012.GIF" wi="22" he="17" />为矩阵宽度,<img file="2013102253768100001DEST_PATH_IMAGE014.GIF" wi="20" he="18" />为矩阵高度;步骤3. 根据可控核梯度求解方法,对步骤2中得到的像素点矩阵求解各个邻接像素的梯度值,记为<img file="2013102253768100001DEST_PATH_IMAGE016.GIF" wi="41" he="25" />,其中<img file="DEST_PATH_IMAGE018.GIF" wi="237" he="22" />;步骤4.根据能量函数,建立二值化的能量图模型;步骤5.计算流量参数<img file="DEST_PATH_IMAGE020.GIF" wi="70" he="25" />;<img file="DEST_PATH_IMAGE022.GIF" wi="175" he="82" />,<img file="DEST_PATH_IMAGE024.GIF" wi="243" he="44" />,<img file="DEST_PATH_IMAGE026.GIF" wi="204" he="45" /><img file="DEST_PATH_IMAGE028.GIF" wi="242" he="48" /><img file="DEST_PATH_IMAGE030.GIF" wi="111" he="50" />其中,<img file="DEST_PATH_IMAGE032.GIF" wi="68" he="21" />为差补矫正因子,<img file="DEST_PATH_IMAGE034.GIF" wi="17" he="18" />为恒边长值,<img file="DEST_PATH_IMAGE036.GIF" wi="17" he="20" />是目标像素点的数目;<img file="DEST_PATH_IMAGE038.GIF" wi="17" he="25" />为输入待检测的像素点,<img file="DEST_PATH_IMAGE040.GIF" wi="62" he="25" />表示待检测的像素点<img file="433844DEST_PATH_IMAGE038.GIF" wi="17" he="25" />是目标像素点的概率,<img file="DEST_PATH_IMAGE042.GIF" wi="66" he="25" />表示待检测的像素点<img file="677743DEST_PATH_IMAGE038.GIF" wi="17" he="25" />是背景像素点的概率;<img file="DEST_PATH_IMAGE044.GIF" wi="18" he="25" />为输出的观测像素点,<img file="DEST_PATH_IMAGE046.GIF" wi="56" he="25" />为将检测到的像素点<img file="818481DEST_PATH_IMAGE044.GIF" wi="18" he="25" />判定为目标像素点的概率,<img file="DEST_PATH_IMAGE048.GIF" wi="58" he="25" />表示将检测到的像素点<img file="780621DEST_PATH_IMAGE044.GIF" wi="18" he="25" />判定为背景像素点的概率;步骤6.算得流量参数后,即得边流量加权参数值<img file="DEST_PATH_IMAGE050.GIF" wi="53" he="25" />和<img file="DEST_PATH_IMAGE052.GIF" wi="21" he="25" />,将<img file="35147DEST_PATH_IMAGE050.GIF" wi="53" he="25" />和<img file="133553DEST_PATH_IMAGE052.GIF" wi="21" he="25" />代入步骤4所建立的能量图模型,对其所建图模型的边流量进行加权,最后依据最大流最小切生长树搜索原理对能量图的流量进行优化分配,实现目标与背景的分割,继而得到所要检测的雷达弱目标;步骤7.判断所得到的目标是否为最后一帧,若是,则跳到步骤9结束检测,反之进入步骤8;步骤8.参数预测及背景更新;首先,结合kalman预测,对第<img file="DEST_PATH_IMAGE054.GIF" wi="14" he="20" />帧视频,定义目标像素点数量为,目标-背景像素对数量为<img file="DEST_PATH_IMAGE056.GIF" wi="14" he="15" />,<img file="DEST_PATH_IMAGE058.GIF" wi="15" he="18" />和<img file="DEST_PATH_IMAGE060.GIF" wi="14" he="18" />分别为<img file="DEST_PATH_IMAGE062.GIF" wi="15" he="15" />和<img file="DEST_PATH_IMAGE056A.GIF" wi="14" he="15" />的速度分量,则得系统的状态向量为<img file="DEST_PATH_IMAGE065.GIF" wi="137" he="18" />,同时,定义<img file="DEST_PATH_IMAGE067.GIF" wi="21" he="16" />为第<img file="DEST_PATH_IMAGE069.GIF" wi="8" he="11" />帧视频的第一个观测量,<img file="DEST_PATH_IMAGE071.GIF" wi="21" he="16" />为第<img file="DEST_PATH_IMAGE069A.GIF" wi="8" he="11" />帧视频的第二个观测量,则得系统的观测向量为<img file="DEST_PATH_IMAGE073.GIF" wi="102" he="20" />,建立如下系统模型为:状态方程:<img file="DEST_PATH_IMAGE075.GIF" wi="106" he="15" />;观测方程:<img file="DEST_PATH_IMAGE077.GIF" wi="86" he="15" />,其中,<img file="DEST_PATH_IMAGE079.GIF" wi="14" he="11" />为系统噪声,满足<img file="DEST_PATH_IMAGE081.GIF" wi="58" he="15" />,<img file="DEST_PATH_IMAGE083.GIF" wi="82" he="18" />;<img file="DEST_PATH_IMAGE085.GIF" wi="9" he="11" />为观测噪声,满足<img file="DEST_PATH_IMAGE087.GIF" wi="54" he="15" />,<img file="DEST_PATH_IMAGE089.GIF" wi="78" he="18" />;设定系统噪声相关矩阵为<img file="DEST_PATH_IMAGE091.GIF" wi="15" he="15" />,观测噪声相关矩阵为<img file="DEST_PATH_IMAGE093.GIF" wi="14" he="15" />,则对应的方程参数分别为:<img file="DEST_PATH_IMAGE095.GIF" wi="234" he="157" />,<img file="DEST_PATH_IMAGE097.GIF" wi="90" he="42" />,<img file="DEST_PATH_IMAGE099.GIF" wi="137" he="90" />,<img file="DEST_PATH_IMAGE101.GIF" wi="119" he="42" />;其中,<img file="DEST_PATH_IMAGE103.GIF" wi="14" he="11" />为系统读取每一帧图像信息所需要的时间,<img file="DEST_PATH_IMAGE105.GIF" wi="14" he="15" />为系统噪声,<img file="DEST_PATH_IMAGE107.GIF" wi="12" he="11" />为观测噪声;根据上述系统模型,结合下述公式:Kalman一步预测增益:<img file="DEST_PATH_IMAGE109.GIF" wi="183" he="18" />,预测状态方程:<img file="DEST_PATH_IMAGE111.GIF" wi="190" he="19" />,预测均方误差:<img file="DEST_PATH_IMAGE113.GIF" wi="190" he="18" />,通过一步迭代,实现对<img file="DEST_PATH_IMAGE115.GIF" wi="28" he="11" />帧检测视频的目标像素点个数<img file="DEST_PATH_IMAGE117.GIF" wi="24" he="15" />及目标-背景像素对数<img file="DEST_PATH_IMAGE119.GIF" wi="23" he="15" />的预测,并更新步骤5所提到的流量参数<img file="DEST_PATH_IMAGE121.GIF" wi="61" he="17" />;接着,按照下列公式,更新背景参数:权值<img file="DEST_PATH_IMAGE123.GIF" wi="31" he="20" />,均值<img file="DEST_PATH_IMAGE125.GIF" wi="30" he="20" />及方差<img file="DEST_PATH_IMAGE127.GIF" wi="48" he="21" />;<img file="DEST_PATH_IMAGE129.GIF" wi="236" he="65" /><img file="DEST_PATH_IMAGE131.GIF" wi="141" he="20" /><img file="DEST_PATH_IMAGE133.GIF" wi="166" he="20" /><img file="DEST_PATH_IMAGE135.GIF" wi="270" he="21" />其中,<img file="DEST_PATH_IMAGE137.GIF" wi="11" he="7" />为选用的混合高斯模型的个数,<img file="DEST_PATH_IMAGE139.GIF" wi="109" he="21" />表示<img file="DEST_PATH_IMAGE141.GIF" wi="4" he="9" />时刻第<img file="DEST_PATH_IMAGE069AA.GIF" wi="8" he="11" />帧视频的第<img file="DEST_PATH_IMAGE144.GIF" wi="7" he="13" />个高斯模型的高斯分布概率,<img file="DEST_PATH_IMAGE146.GIF" wi="80" he="14" />,<img file="DEST_PATH_IMAGE148.GIF" wi="53" he="14" />为更新因子,<img file="DEST_PATH_IMAGE150.GIF" wi="11" he="13" />为第t时刻检测的像素点;最后,在完成参数预测及背景更新后,读入第<img file="DEST_PATH_IMAGE115A.GIF" wi="28" he="11" />帧视频图像,并返回步骤2;步骤9.结束检测。
地址 310018 浙江省杭州市下沙高教园区2号大街