发明名称 基于提升小波变换和像素区域相关性的图像融合方法
摘要 本发明公开了一种基于提升小波变换和像素区域相关性的图像融合方法,包括以下步骤:提升小波变换;低频子带融合系数的计算;高频子带融合系数的计算;低频、高频子带的系数融合;小波逆变换。本发明对于低频子带的每个系数,根据其自相关特征,通过比较其协方差相关系数来确定融合系数。对于高频子带的每个系数,根据其所在子带内系数分布的方向特征,以及其所在同一方向子带间的系数分布的四叉树结构特征,通过比较其子带内方向邻域的匹配度和子带间四叉树结构的匹配度,从而确定其融合系数。测试结果表明本发明的融合方法适用于多聚焦图像和医学图像的融合。
申请公布号 CN103065291A 申请公布日期 2013.04.24
申请号 CN201210573781.4 申请日期 2012.12.26
申请人 辽宁师范大学 发明人 王相海;周志光;宋传鸣;苏欣
分类号 G06T5/50(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 大连非凡专利事务所 21220 代理人 闪红霞
主权项 1.一种基于提升小波变换和像素区域相关性的图像融合方法,其特征在于按如下步骤进行:a. 将待融合的两幅图像分别进行相同级数的提升方案小波变换;b. 获得低频子带的所有融合系数:b.1对于每个待融合的低频系数,计算每幅待融合图像对应位置上低频系数的八邻域方差<img file="180849DEST_PATH_IMAGE001.GIF" wi="21" he="24" />:<img file="447882DEST_PATH_IMAGE003.GIF" wi="239" he="51" />b.2计算待融合系数对应的<img file="584466DEST_PATH_IMAGE005.GIF" wi="68" he="27" />和沿着水平方向的一步相关系数<img file="253344DEST_PATH_IMAGE007.GIF" wi="14" he="19" /><sub>1</sub>:<img file="184391DEST_PATH_IMAGE009.GIF" wi="420" he="49" />,<img file="938721DEST_PATH_IMAGE011.GIF" wi="144" he="30" />b.3计算待融合系数对应的<img file="878995DEST_PATH_IMAGE013.GIF" wi="68" he="27" />和沿着竖直方向的一步相关系数<img file="402380DEST_PATH_IMAGE007.GIF" wi="14" he="19" /><sub>2</sub>:<img file="504328DEST_PATH_IMAGE015.GIF" wi="143" he="30" />b.4计算待融合系数对应的<img file="683637DEST_PATH_IMAGE017.GIF" wi="64" he="25" />和沿着对角线方向的一步相关系数<img file="224339DEST_PATH_IMAGE007.GIF" wi="14" he="19" /><sub>3</sub>:<img file="539914DEST_PATH_IMAGE019.GIF" wi="299" he="36" />,<img file="140660DEST_PATH_IMAGE021.GIF" wi="143" he="33" />b.5计算待融合系数的一步自相关系数<img file="541685DEST_PATH_IMAGE006.GIF" wi="15" he="21" />:<img file="886079DEST_PATH_IMAGE022.GIF" wi="177" he="40" />b.6比较待融合系数的一步自相关系数<img file="56160DEST_PATH_IMAGE006.GIF" wi="15" he="21" />的大小,最终确定融合系数<img file="827807DEST_PATH_IMAGE023.GIF" wi="73" he="30" />:<img file="716129DEST_PATH_IMAGE025.GIF" wi="312" he="63" />b.7遍历所有低频子带系数,重复执行b.1~b.6,最终获得低频子带的所有融合系数; c. 获得高频子带的所有融合系数:c.1对于高频子带的每一个系数,按照方向模板,用下述公式计算待融合系数的方向邻域匹配度:<img file="598634DEST_PATH_IMAGE026.GIF" wi="248" he="52" />,<img file="888801DEST_PATH_IMAGE028.GIF" wi="312" he="73" />,其中,<img file="769033DEST_PATH_IMAGE029.GIF" wi="87" he="36" />表示图像<i>k</i>中以(<i>x,y</i>)为中心、沿着<img file="206967DEST_PATH_IMAGE030.GIF" wi="17" he="18" />方向的方向邻域能量,<img file="830847DEST_PATH_IMAGE032.GIF" wi="69" he="24" />表示待融合图像以<img file="37837DEST_PATH_IMAGE033.GIF" wi="54" he="25" />为中心的<img file="88970DEST_PATH_IMAGE030.GIF" wi="17" he="18" />方向邻域的匹配度;c.2对于高频子带的每一个系数,利用如下公式计算待融合系数所在四叉树结构的匹配度:<img file="14200DEST_PATH_IMAGE035.GIF" wi="228" he="46" />,<img file="176191DEST_PATH_IMAGE037.GIF" wi="311" he="71" />,其中,<img file="DEST_PATH_IMAGE039.GIF" wi="60" he="22" />表示图像<i>k</i>中<img file="458529DEST_PATH_IMAGE033.GIF" wi="54" he="25" />所在的四叉树<i>T</i>的能量,<i>M</i><sub><i>T</i></sub>(<i>x,y</i>)表示待融合图像<img file="477300DEST_PATH_IMAGE033.GIF" wi="54" he="25" />所在四叉树的匹配度;c.3对于待融合图像高频子带的同一位置的系数,比较其方向邻域匹配度和四叉树结构的匹配度的大小,进而确定融合算子,即若方向邻域的匹配度大,则转入c.4;若四叉树结构的匹配度大,则执行转入c.5;c.4设匹配度阈值为<img file="DEST_PATH_IMAGE041.GIF" wi="15" he="21" />;若<img file="DEST_PATH_IMAGE043.GIF" wi="104" he="28" />,则<img file="DEST_PATH_IMAGE045.GIF" wi="256" he="59" />;否则,<img file="DEST_PATH_IMAGE047.GIF" wi="391" he="49" />,其中,<img file="DEST_PATH_IMAGE049.GIF" wi="204" he="74" />c.5设匹配度阈值为<img file="702876DEST_PATH_IMAGE040.GIF" wi="16" he="24" />,若<img file="DEST_PATH_IMAGE051.GIF" wi="93" he="25" />,则<img file="DEST_PATH_IMAGE053.GIF" wi="247" he="54" />;否则,<img file="DEST_PATH_IMAGE055.GIF" wi="408" he="59" />,其中,<img file="DEST_PATH_IMAGE057.GIF" wi="215" he="63" />c.6按照上述步骤,重复执行b.1~b.5,最终获得高频子带的所有融合系数<i>f</i>(<i>x,y</i>);d.按低频子带融合规则对b步骤所获得的低频子带的所有融合系数进行融合,进而获得低频子带融合系数;按高频子带融合规则对c步骤所获得的高频子带的所有融合系数进行融合,获得高频子带融合系数;e. 对所得到的低频子带融合系数和高频子带融合系数进行提升方案小波的逆变换,进而获得融合图像。
地址 116000 辽宁省大连市沙河口区黄河路850号