发明名称 基于脉冲发送皮层模型的医学图像融合方法
摘要 本发明公开了一种基于脉冲发放皮层模型的医学图像融合方法,具体为:首先对图像作归一化处理,然后将图像送入脉冲发放皮层模型得到每幅待融合图像的点火映射图,对两幅点火映射图中每一待考虑像素,计算以两对应像素为中心的两图像块的欧氏距离及每一图像块的像素值的总和;最后利用基于上述参数的判决算子,对两待融合图像进行选择或加权运算,从而实现图像的有效融合。本发明将脉冲发放皮层模型用于医学图像融合,同时在融合中使用了图像块的思路,能更好地反映源图像的真实信息,弥补了现有脉冲耦合神经网络用于图像融合计算复杂的不足。
申请公布号 CN102867296B 申请公布日期 2015.08.19
申请号 CN201210296082.X 申请日期 2012.08.20
申请人 华中科技大学 发明人 张旭明;吴意;王瑞;李柳;丁明跃;熊有伦;尹周平;王瑜辉
分类号 G06T5/50(2006.01)I;G06T7/00(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 华中科技大学专利中心 42201 代理人 李智
主权项 基于脉冲发送皮层模型的医学图像融合方法,具体为:步骤1 对两幅待融合图像A<sup>0</sup>和B<sup>0</sup>进行归一化处理得到图像A和B;步骤2 分别按照如下相同方式获取两图像A和B的点火映射图A′和B′:步骤21 初始化迭代次数n=1,阈值E<sub>ij</sub>(1)为零矩阵,内部活动U<sub>ij</sub>(1)为零矩阵;步骤22 计算待融合图像的点火映射图<img file="FDA0000675396510000011.GIF" wi="634" he="136" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>U</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>fU</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>S</mi><mi>ij</mi></msub><munder><mi>&Sigma;</mi><mi>kl</mi></munder><msub><mi>W</mi><mi>ijkl</mi></msub><msub><mi>Y</mi><mi>kl</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>S</mi><mi>ij</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000675396510000012.GIF" wi="917" he="112" /></maths>Y<sub>ij</sub>(n)为第n次迭代中像素点(i,j)的输出脉冲,Y<sub>kl</sub>(n‑1)为第n‑1次迭代中像素点(k,l)的输出脉冲,S<sub>ij</sub>为外部对像素点(i,j)的神经元的刺激即像素点(i,j)的灰度值,U<sub>ij</sub>(n)为第n次迭代中像素点(i,j)的神经元的内部活动,W<sub>ijkl</sub>为像素点(i,j)的神经元与像素点(k,l)间的神经元间连接的突触权重,(k,l)≠(i,j),f为神经元内部活动的衰减系数;步骤23 若n小于预定迭代次数阈值,则n=n+1,E<sub>ij</sub>(n)=gE<sub>ij</sub>(n‑1)+hY<sub>ij</sub>(n‑1),h为阈值放大系数,g为阈值衰减系数,返回步骤22,否则,进入步骤3;步骤3 利用点火映射图进行像素级的融合:步骤31 在两点火映射图A′和B′中,分别逐像素点提取以该像素点为中心的子块<img file="FDA0000675396510000013.GIF" wi="77" he="77" />和<img file="FDA0000675396510000014.GIF" wi="103" he="77" />步骤32 令z=1;步骤33 计算由子块<img file="FDA0000675396510000021.GIF" wi="73" he="77" />和<img file="FDA00006753965100000210.GIF" wi="73" he="75" />融合得到的像素点(i,j)的像素值R(i,j):<img file="FDA0000675396510000022.GIF" wi="1179" he="383" />其中,I<sub>A′</sub>(i,j)为子块<img file="FDA0000675396510000023.GIF" wi="76" he="76" />的中心像素点(i,j)对应的图像A中的像素点的像素值,I<sub>B′</sub>(i,j)为子块<img file="FDA0000675396510000024.GIF" wi="75" he="76" />的中心像素点(i,j)对应的图像B中的像素点的像素值,<img file="FDA0000675396510000025.GIF" wi="67" he="70" />和<img file="FDA0000675396510000026.GIF" wi="64" he="70" />分别为子块<img file="FDA0000675396510000027.GIF" wi="74" he="77" />和<img file="FDA0000675396510000028.GIF" wi="80" he="77" />的像素值总和,T0为子块<img file="FDA00006753965100000211.GIF" wi="73" he="75" />和<img file="FDA00006753965100000212.GIF" wi="78" he="80" />的欧式距离,th为判决算子阈值;步骤34 若z小于子块总数,则z=z+1,返回步骤33,否则,结束。
地址 430074 湖北省武汉市洪山区珞瑜路1037号