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);

%==========================================================

%对图像进行增强处理

for i=1:Csize1(2)

    c1(i)=1.2*c1(i);

end

Csize2=size(c2);

for j=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('融合后增强图像');