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

请注意该程序处理的图像必须是灰度图像，不能是彩图。 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(融合后增强图像); 图 3.31 CT和MRI图像融合结果图 小波图像融合的Matlab实现示例（添加图片演示080428） 默认分类 2008-06-14 20:07:21 阅读68 评论0 ??字号：大中小?订阅 前几天我们讨论了一种基于低频融合策略的小波图像融合算法，今天就来说说算法编程实现的几个细节问题。 （1）高频带的融合 ??????? 首先，高频系数融合均采用基于像素点绝对值取大的规则，这里要注意的是，比较的对象是两幅原始图像的像素点的绝对值，而融合后图像的像素点的取值是两幅原始图像像素点绝对值大的那一点的值，也就是说，要取绝对值大的那个点的实际值。故相应的Matlab程序为： % I1，I2 为两幅原始图像 y1=mywavedec2(I1,dim); y2=mywavedec2(I2,dim); [r,c]=size(y1);??? ??? ??? % 根据低频融合算法进行图像融合 for i=1:r??? ??? ??? % 首先取两幅源图像相应的小波分解系数绝对值最大者的值作为融合图像的分解系数 ??? for j=1:c ??????? if ( abs(y1(i,j)) = abs(y2(i,j)) ) ??????????? y3(i,j)=y1(i,j); ??????? elseif ( abs(y1(i,j)) abs(y2(i,j)) ) ??????????? y3(i,j)=y2(i,j); ??????? end %??? y3(i,j)=max( abs(y1(i,j)) , abs(y2(i,j)) );? % 这行代码是错误的，取了绝对值而非实际值 ??? end end （2）低频融合的对象 ??????? 低频融合的对象是N级分解后图像的LL-N子图像，相应的程序代码为： LLa=y1(1:r/(2^dim),1:c/(2^dim));??? % 调用lowfrefus函数对低频部分的小波分解系数进行融合 LLb=y