双通道pcnn图像融合程序

function MODPCNNFUSION %%%% THIS FUNCTION IS USED TO FUSE IMAGES BASED ON %%%% "MEDICAL IMAGE FUSION USING mPCNN" BY WANGZHAOBIN %%%% REFERENCE: Medical Image Fusion Using mPCNN %%%% Z Wang, Y Ma [lzu] Information Fusion.2008 %%%% link:http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6W76-4N %%%% HV4MP-3&_user=1002898&_rdoc=1&_fmt=&_orig=search&_sort=d&view=c&_acct=C000050164&_version=1&_urlVersion=0&_userid=1002898&md5=05e1ab129356e043049a4636b1810bd7 %%%% Editor:Denghaibo Email:denghblzu07@gmail.com %%%% date:2009-4-24 IMA=double(imread('disk1.gif')); IMB=double(imread('disk2.gif')); [ROW,COL]=size(IMA); [ROW1,COL1]=size(IMB); if(ROW~=ROW1|COL~=COL1) error('input images should be the same size!'); end %%%% SET THE mPCNN PARAMEERS [7 IN TOTAL] %%%% W=M; BETA1=BETA2; DELTA; ALPHAT; VT; BETA=0.5; %%% note BETA1=BETA2; SO WE USE BETA FOR CONVENIENCE M=[0.707,1,0.707;1,0,1;0.707,1,0.707]; DELTA=1; ALPHAT=0.012; VT=4000; %%% MAKE SURE WHETHER THIS IS RIGHT. %%%%% IMAWT=conv2(IMA,M,'same'); %%% 卷积后保持源图像大小； IMBWT=conv2(IMB,M,'same'); SUMPIXEL=ROW*COL; %%% PIXELS NUMBER OF FUSED IMAGE JUDFIR=zeros(ROW,COL); %%% JUDFIR IS USED TO JUDGE WHETHER THE NEURON(I,J) IS FIRED %%%JUDFIR用于判断神经元（i，j）是否激发。每个神经元在所有的迭代中只激发一次。 THR=500*ones(ROW,COL); %%% 不知此阈值是否合适 U=zeros(ROW,COL); %%% ITERATION ITETIME=10; %%% SET ITERATION TIME for(itime=1:10) for(i=1:ROW) for(j=1:COL) if(JUDFIR(i,j)~=1) %%% THE NEURON SHOULD NOT BE FIRED YET. %%% 对单点进行迭代计算 HA=IMAWT(i,j)+IMA(i,j); HB=IMBWT(i,j)+IMB(i,j); U(i,j)=(1+BETA*HA)*(1+BETA*HB)+DELTA; %%%判断是否激发 if (U(i,j)>THR(i,j)) %% 满足条件则激发 JUDFIR(i,j)=1; else %% 不激发，对阈值进行衰减 THR(i,j)=THR(i,j)*exp(-ALPHAT); end end %%% end of JUDFIR end %%% end j end %%% end i if(sum(sum(JUDFIR))==SUMPIXEL) fprintf(1,'fused successfully'); break end end