欢迎来到天天文库
浏览记录
ID:21478933
大小:26.00 KB
页数:6页
时间:2018-10-22
《基于matlab的多模态医学图像融合仿真》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、基于Matlab的多模态医学图像融合仿真 摘要多模态医学图像融合是一种集医学图像处理和医学影像诊断于一体的新技术,其迅猛发展对医学影像技术的进步和临床诊疗产生了深远的影响。本文使用几种多模态医学图像融合的典型方法,在Matlab环境下进行了仿真实验,为相关领域的专业人员和非专业关注者提供了参考。 【关键词】多模态医学图像融合Matlab医学图像处理医学影像学 1引言 随着计算机技术的飞速发展和信息时代的到来,医学成像成为了现代医疗技术不可缺少的一部分。由于各种成像设备的成像原理有所差异,不同模态的图像有各自的优点和局限性。在此背
2、景下,充分利用现有的成像设备,研究一种能整合来自不同成像设备的图像信息并将其作为一个整体加以表达的图像融合技术,受到了相关领域的高度重视。本文分析了多模态医学图像融合的典型方法,并在Matlab环境下进行了仿真实验,为相关领域的专业人员和非专业关注者提供了参考。 2多模态医学图像的融合方法 医学图像融合方法可分为像素级、特征级和决策级三个层次,目前应用比较广泛的是像素级层次融合融合方法,同时它又是后两种融合方法的基础,本文主要研究像素级层融合方法。 根据其特点、融合步骤和基本原理,又可将其分为基于空间域融合方法和基于变换域融合方法
3、两大类。这两类方法不是相互独立的关系,在许多算法中将两者结合使用可以达到更好的融合效果。 3多模态医学图像融合matlab仿真程序 本文针对以下几种融合方法,进行了多模态医学图像融合的仿真实验,实验所用为两幅已经做完配准处理的MRI和CT图像,如图1所示。 下面将详细阐述程序实现和实验结果,如图2所示。 3.1像素灰度值极大/极小融合法 %图像像素灰度值极大法 fori=1:m1 forj=1:n1 if(abs(M1(i,j))>=abs(M2(i,j))) M3(i,j)=M1(i,j); elseif(abs(
4、M1(i,j))5、 %傅里叶反变换 M3=ifft2(y3); M4=ifft2(y4); %数据类型转换 M3=im2uint8(M3); M4=im2uint8(M4); 3.4小波变换法 zt=2;%小波分解层数 wtype='db1';%使用的小波类型 %小波分解 [c0,s0]=wavedec2(M1,zt,wtype);%多尺度二维小波分解 [c1,s1]=wavedec2(M2,zt,wtype);%多尺度二维小波分解 %小波系数简单加权法 c=(c0+c1)*0.5; %高频部分系数选择绝对值极大法,低频部分6、系数采用二者求平均的方法 KK=size(c1); Coef_Fusion1=zeros(1,KK(2)); %低频系数的处理 Coef_Fusion1(1:s1(1,1))=(c0(1:s1(1,1))+c1(1:s1(1,1)))/2; %高频系数的处理 MM1=c0(s1(1,1)+1:KK(2)); MM2=c1(s1(1,1)+1:KK(2)); mm=(abs(MM1))>(abs(MM2)); Y=(mm.*MM1)+((~mm).*MM2); Coef_Fusion1(s1(1,1)+1:KK(2))7、=Y; %小波重构 Y1=waverec2(c,s0,wtype); Y2=waverec2(Coef_Fusion1,s0,wtype); 4结语 多模态医学图像融合技术是医学图像处理技术的一个重要分支,在临床诊疗、计算机辅助诊断、远程医疗、放射治疗及手术计划的制定等方面有着广泛的应用前景,对医学影像的进步和发展有重要的促进作用。在此背景下,本文分析了多模态医学图像融合的典型方法,并给出了Matlab环境下的仿真实例,对相关领域的专业人员和非专业关注者具有一定的参考价值。 参考文献 [1]李伟,朱学峰.医学图像融合技术及8、其应用[J].中国医学影像技术,2005,21(07):1126-1129. [2]田娅,饶妮妮,蒲立新.国内医学图像处理技术的最新动态[J].电子科技大学学报,2002,31(05):485-489.
5、 %傅里叶反变换 M3=ifft2(y3); M4=ifft2(y4); %数据类型转换 M3=im2uint8(M3); M4=im2uint8(M4); 3.4小波变换法 zt=2;%小波分解层数 wtype='db1';%使用的小波类型 %小波分解 [c0,s0]=wavedec2(M1,zt,wtype);%多尺度二维小波分解 [c1,s1]=wavedec2(M2,zt,wtype);%多尺度二维小波分解 %小波系数简单加权法 c=(c0+c1)*0.5; %高频部分系数选择绝对值极大法,低频部分
6、系数采用二者求平均的方法 KK=size(c1); Coef_Fusion1=zeros(1,KK(2)); %低频系数的处理 Coef_Fusion1(1:s1(1,1))=(c0(1:s1(1,1))+c1(1:s1(1,1)))/2; %高频系数的处理 MM1=c0(s1(1,1)+1:KK(2)); MM2=c1(s1(1,1)+1:KK(2)); mm=(abs(MM1))>(abs(MM2)); Y=(mm.*MM1)+((~mm).*MM2); Coef_Fusion1(s1(1,1)+1:KK(2))
7、=Y; %小波重构 Y1=waverec2(c,s0,wtype); Y2=waverec2(Coef_Fusion1,s0,wtype); 4结语 多模态医学图像融合技术是医学图像处理技术的一个重要分支,在临床诊疗、计算机辅助诊断、远程医疗、放射治疗及手术计划的制定等方面有着广泛的应用前景,对医学影像的进步和发展有重要的促进作用。在此背景下,本文分析了多模态医学图像融合的典型方法,并给出了Matlab环境下的仿真实例,对相关领域的专业人员和非专业关注者具有一定的参考价值。 参考文献 [1]李伟,朱学峰.医学图像融合技术及
8、其应用[J].中国医学影像技术,2005,21(07):1126-1129. [2]田娅,饶妮妮,蒲立新.国内医学图像处理技术的最新动态[J].电子科技大学学报,2002,31(05):485-489.
此文档下载收益归作者所有