当前位置: 首页 > 专利查询>吉林大学专利>正文

一种基于多级阈值分割的多模态医学图像处理方法技术

技术编号:12996312 阅读:133 留言:0更新日期:2016-03-10 11:10
本发明专利技术公开了一种基于多级阈值分割的多模态医学图像处理方法,包括如下步骤:步骤一:获取相应的医学图像灰度图像数据;步骤二:构建高斯尺度空间模型来同时完成去噪处理;步骤三:采用各向异性非线性扩散和0范式梯度最小化模型来完成医学图像的平滑处理;步骤四:缩小OTSU图像分割算法中最佳阈值的搜算范围,以提高搜算效率;步骤五:通过有效利用空间位置信息来提取完整的分割区域,从而得出更好的分割效果。本发明专利技术在图像预处理去噪阶段建立多尺度高斯模型,具有相比单一尺度模型下更好地去噪效果,在采用多阈值OTSU分割方法这一阶段,缩小了算法阈值搜素范围,提高了分割算法运行效率。因此本发明专利技术具有实时性好,分割效果准确这两大优点。

【技术实现步骤摘要】
一种基于多级阈值分割的多模态医学图像处理方法
本专利技术涉及医学图像处理领域,特别涉及一种基于多级阈值分割的多模态医学图像处理方法。
技术介绍
作为各种医学图像处理技术如配准、融合等的基础,同时鉴于在当前的临床医学应用如临床辅助诊断、图像引导的外科手术以及放射治疗中,医学图像分割技术显示出越来越重要的临床价值。医学图像分割技术是医学图像处理的研究基础,如多模态医学图像配准、多模态医学图像融合等,其分割质量的好坏将直接决定后续配准、融合等操作能否有效的进行。在当前的临床医学研究及应用如临床辅助诊断、图像引导的外科手术和放射治疗中,医学图像分割技术发挥着极其重要的临床研究及应用价值。更为重要的是,随着医学成像技术的快速发展,医学影像已成为医护工作者们的第三只眼睛,为临床医生们提供了有效的诊断线索。当前,各种成像模式如核磁共振成像(MRI)、计算机断层成像(CT)、超声成像(US)、正电子发射计算机断层成像(PET)等已在临床医学中得到了广泛应用,因此对多模态医学图像分割技术的研究已成必然趋势。而与现有的分割技术如基于偏微分方程的分割技术、基于图论的分割技术等相比,阈值分割技术具有计算简单、高效等特点,满足临床医学中的实时性要求。
技术实现思路
本专利技术的目的是提供一种基于多级阈值分割的多模态医学图像处理方法,对多级阈值进行分割,使阈值搜素范围减少从而提高医学图像分割效率,以获得更准确的医学图像分割结果。本专利技术提供的技术方案为:一种基于多级阈值分割的多模态医学图像处理方法,包括如下步骤:步骤一:获取相应的医学图像灰度图像数据;步骤二:构建高斯尺度空间模型来完成去噪处理;步骤三:采用各向异性非线性扩散和0范式梯度最小化模型对步骤二处理后的医学图像的平滑处理;步骤四:缩小OTSU图像分割算法中最佳阈值的搜算范围;步骤五:利用空间位置信息来提取完整的分割区域,从而得出更好的分割效果。优选的是,步骤二包括如下分步骤:a、采用高斯函数与图像灰度函数的卷积得到高斯尺度空间模型:L(x,y,σ)=G(x,y,σ)*I(x,y)其中,*表示卷积操作,σ为尺度因子,G(x,y,σ)为高斯函数,I(x,y)为所处理的灰度图像,x,y分别为相应像素点医学图像的横坐标和纵坐标位置;b、计算第i层高斯尺度空间信息Li(x,y,σi)=G(x,y,σi)*I(x,y),σi+1=kσi,k>1,i=1,2,3...其中,σ0为初始尺度因子,σi为第i层尺度因子;k为平滑系控制参数;c、将高斯尺度空间模型作为背景信息,采用背景差方法,获取第i层目标图像信息:Di(x,y)=|I(x,y)-Li(x,y,σi)|d、对各层目标图像信息采用加权平均方式来获取最终目标图像信息:其中ωi为第i层目标图像信息所占权重,n为高斯尺度空间中图像总层数。优选的是,步骤d中,所述第i层目标图像信息所占权重ωi满足:优选的是,步骤三中,采用各向异性非线性扩散模型进行图像的平滑处理采用的模型满足:其中,It为平滑后图像灰度函数,t为平滑次数,为梯度算子,Δ表示拉普拉斯算子,div(·)表示散度算子,c(x,y,t)为扩散系数。优选的是,所述扩散系数c(x,y,t)满足:其中,或K为控制参数。优选的是,步骤三中,使用0范式梯度最小化模型时引入一个约束项C(S)来避免弱边缘泄露现象,所述约束项C(S)满足:S为平滑后图像,为x方向梯度,为y方向梯度,Sp为概率为p时的图像,p为对应像素点出现的概率。优选的是,步骤四中,多级最佳阀值满足:其中,为区域Ck概率总和,为区域Ck内的一阶矩总和,pi为灰度图像中灰度值为i的像素点出现的概率,l为最佳阈值个数,L为灰度级。优选的是,步骤五中,采用漫水法进行连通区域查找,同时将零碎小区域合并到周围与其最相似的大区域中,采用区域信息和边界相结合的计算策略,从而最终形成相似信息的分割区域。本专利技术的有益效果是:本专利技术提供了一种基于多级阈值分割的多模态医学图像处理方法,在图像预处理去噪阶段建立多尺度高斯模型,由于在多尺度下图像各项信息相互结合,为图像去噪和非均匀光照矫正均提供了保证,具有相比单一尺度模型下更好地去噪效果。在采用多阈值OTSU分割方法这一阶段,对OTSU图像分割算法做出改进,缩小了算法阈值搜素范围,从而提高分割算法运行效率,使分割结果成像的实时性变得更好。本专利技术通过对多机阈值分割后的结果进行空间位置信息分析与处理,从而确定有效的分割联通域,进而获取更好的分割效果。因此本专利技术具有实时性好,分割效果准确这两大优点。附图说明图1为本专利技术所述的基于多级阈值分割的多模态医学图像处理方法流程图。具体实施方式下面结合附图对本专利技术做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。如图1所示,本专利技术提供了一种基于多级阈值分割的多模态医学图像处理方法,通过对医学图像建立多尺度高斯模型,采用各向异性非线性扩散和0范式梯度最小化模型进行去噪和平滑等预处理,并对多级阈值分割OTSU算法进行改进,使阈值搜素范围减少从而提高医学图像分割效率,并对分割后的图像进行空间位置信息分析与处理来获得更准确的医学图像分割结果。具体步骤如下:步骤一:获取相应的医学图像灰度图像数据。通过核磁共振成像(MRI)、计算机断层成像(CT)、超声成像(US)、正电子发射计算机断层成像(PET)等方式得到医学图像,读取原始图像各像素点灰度值,归一化为0-255,并以I(x,y)表示相应灰度图像,x,y分别为相应像素点医学图像的横坐标和纵坐标位置。即I(x,y)表示在数据矩阵中第x行第y列的像素灰度值。步骤二S120:构建高斯尺度空间并进行去噪处理。鉴于现有的非均匀光照矫正和去噪通常仅限于单一尺度下进行,并且分开进行,带来的问题不仅是去噪效果不够理想,而且非均匀光照信息也没有得到有效矫正。为了弥补该缺陷,本专利技术提出了一种高斯多尺度背景模型。高斯函数在计算机视觉以及模式识别领域中有着十分广泛的应用,其定义见下式所示:其中,σ为尺度因子。对于灰度函数为I(x,y)的医学图像,采用高斯函数与该图像灰度函数的卷积来表示该医学图像对应的高斯尺度空间L(x,y,σ):L(x,y,σ)=G(x,y,σ)*I(x,y)其中,*表示卷积操作。记初始尺度因子为σ0,第i层尺度因子为σi,通过上述尺度空间的定义,则第i层尺度信息Li(x,y,σi)计算如下:Li(x,y,σi)=G(x,y,σi)*I(x,y)σi+1=kσi,k>1,i=1,2,3...其中,k为平滑系控制参数,用来控制平滑速率。通过将高斯尺度空间视作背景信息,则每层尺度信息Li(x,y,σ)所对应的目标信息Di(x,y)可采用背景差方法获取,如下式所示:Di(x,y)=|I(x,y)-Li(x,y,σi)|为了能够得到一个更加准确的目标图像,并且保留更多的目标细节信息,对各层得到的目标图像采用加权平均方式来获取高斯滤波图像计算过程如下:其中ωi为各层目标图像所占权重,n为高斯尺度空间中图像总层数。由于随着尺度因子σi不断增大提取出的目标图像也更加准确,则它所占权重也将越大,因此采用下式计算各层目标图像的权重值:在尺度空间建立过程中,本专利技术将根据差值图像的变化来控制尺度空间是否停止创建,另外,拟通过γ增强操作本文档来自技高网...
一种基于多级阈值分割的多模态医学图像处理方法

【技术保护点】
一种基于多级阈值分割的多模态医学图像处理方法,其特征在于,包括如下步骤:步骤一:获取医学图像灰度图像数据;步骤二:构建高斯尺度空间模型来完成去噪处理;步骤三:采用各向异性非线性扩散和0范式梯度最小化模型对步骤二处理后的医学图像的平滑处理;步骤四:缩小OTSU图像分割算法中最佳阈值的搜算范围;步骤五:利用空间位置信息来提取完整的分割区域,从而得出更好的分割效果。

【技术特征摘要】
1.一种基于多级阈值分割的多模态医学图像处理方法,其特征在于,包括如下步骤:步骤一:获取医学图像灰度图像数据;步骤二:构建高斯尺度空间模型来完成去噪处理;具体包括:a、采用高斯函数与图像灰度函数的卷积得到高斯尺度空间模型:L(x,y,σ)=G(x,y,σ)*I(x,y)其中,*表示卷积操作,σ为尺度因子,G(x,y,σ)为高斯函数,I(x,y)为所处理的灰度图像,x,y分别为相应像素点医学图像的横坐标和纵坐标位置;b、计算第i层高斯尺度空间信息Li(x,y,σi)=G(x,y,σi)*I(x,y),σi+1=kσi,k>1,i=1,2,3...其中,σ0为初始尺度因子,σi为第i层尺度因子;k为平滑系控制参数;c、将高斯尺度空间模型作为背景信息,采用背景差方法,获取第i层目标图像信息:Di(x,y)=|I(x,y)-Li(x,y,σi)|d、对各层目标图像信息采用加权平均方式来获取最终目标图像信息:其中ωi为第i层目标图像信息所...

【专利技术属性】
技术研发人员:秦俊申铉京冯云丛刘翔陈海鹏
申请(专利权)人:吉林大学
类型:发明
国别省市:吉林;22

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1