主权项 |
基于脉冲发送皮层模型的医学图像融合方法,具体为:步骤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>Σ</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,否则,结束。 |