CT与MRI图像融合matlab程序及结果

请注意该程序处理的图像必须是灰度图像，不能是彩图。Matlab程序：clear%装载原始图像001.mat;002.mat并显示。[X,map]=imread('003.Tif');X1=X;map1=map;figure(2);X1=imread('003.Tif');subplot(2,2,1);imshow(X1);title('图像001');X2=imread('004.Tif');subplot(2,2,2);imshow(X2);title('图像002');X1=double(X1);X2=double(X2);%X1=uint8(X1);%X2=uint8(X2);%=========================================================%对上述两图像进行分解[c1,I1]=wavedec2(X1,3,'sym4');[c2,I2]=wavedec2(X2,3,'sym4');%对分解系数进行融合c=c1+c2;%=========================================================%应用融合系数进行图像重构并显示XX=waverec2(c,I1,'sym4')subplot(2,2,3);%image(uint8(XX));XX=double(XX);image(XX);title('融合图像');Csize1=size(c1);%==========================================================%对图像进行增强处理fori=1:Csize1(2)c1(i)=1.2*c1(i);endCsize2=size(c2);forj=1:Csize2(2)c2(j)=0.8*c2(j);end%===========================================================%通过减小融合技术以减小图像的亮度c=0.5*(c1+c2);%===========================================================%对融合系数进行图像重构XXX=waverec2(c,I2,'sym4');%===========================================================%显示重构结果subplot(2,2,4);XXX=double(XXX);image(XXX);%image(uint8(XXX));title('融合后增强图像');图3.31CT和MRI图像融合结果图