本发明涉及图像处理技术领域,尤其涉及一种超声图像的分析方法及装置。
背景技术:
在临床上,借助于cad(computer-aideddiagnosis,计算机辅助诊断)系统来辅助病症的诊断,可以避免主观因素的影响以及提高诊断结果的准确性和客观性,因此计算机辅助诊断系统被越来越多地应用于临床诊断中。在计算机辅助诊断系统中,一般通过分析医学超声图像来提供诊断结果;例如,对于乳腺肿瘤超声图像,利用cad系统来对该超声图像进行分析以提供是否患乳腺肿瘤、肿瘤位置或大小等诊断结果。而对于计算机辅助诊断系统而言,如何从超声图像中分割出病灶区域(如肿瘤区)是其中的关键一环。
然而,目前对于超声图像的分割算法一般是半自动的,即需要临床医生手工选取感兴趣区域或感兴趣区域的代表点,然后通过计算机进行区域分割。半自动的分割方法降低了cad系统的自动化性能,无法满足临床上日渐产生的海量医学图像处理的需要。而目前的全自动分割算法,大多是引入了一些形状、纹理和空间相对位置等先验约束信息来实现全自动分割。但是,由于病灶的大小、形状千差万别以及边界不清晰等问题,而且即使同一患者由不同的超声设备进行扫描所得的结果也可能不相同,因此难以准确提取乳腺肿瘤的先验约束信息,从而影响分析结果的准确性。
因此,在现有的cad系统中,其针对医学超声图像所采取的分割方式存在效率与准确性难以同时兼顾的问题。
技术实现要素:
本发明实施例提供了一种超声图像的分析方法及装置,能够在保证分割结果准确性的同时,提高分割的自动化水平和速度。
本发明实施例提供了一种超声图像的分析方法,用于从所述超声图像中分割出目标区域,所述分析方法包括:对所述超声图像进行预处理;从预处理后的所述超声图像中提取出所述目标区域的初始轮廓;对所述目标区域的初始轮廓进行演化,从而获取所述目标区域的准确边界。
其中,对所述超声图像进行预处理的步骤包括:选取所述超声图像的有效区域;以及对选取的所述有效区域进行去噪处理。
其中,所述从预处理后的所述超声图像中提取出所述目标区域的初始轮廓的步骤包括:对所述预处理后的所述超声图像进行预分割,得到预分割后的图像;对所述预分割后的图像依次进行一系列处理,其中该一系列处理包括如下至少一项:形态学处理、空洞填充、以及去除与边界相连的区域;提取经过所述一系列处理后的图像中的封闭区域;根据所述提取到的封闭区域的大小,确定所述目标区域的初始轮廓。
其中,对所述预处理后的所述超声图像进行预分割的步骤包括:采用最大类间方差法,来处理所述预处理后的所述超声图像,从而得到第一阈值;根据所述第一阈值,从所述预处理后的所述超声图像中划分出前景图像;采用所述最大类间方差法,来处理所述前景图像,以得到第二阈值;利用所述第二阈值来对所述预处理后的所述超声图像进行预分割,以得到所述预分割后的图像。
其中,对所述预处理后的所述超声图像进行预分割的步骤包括:采用最大类间方差法,来处理所述预处理后的所述超声图像,从而得到第一阈值;根据所述第一阈值,来对所述预处理后的所述超声图像进行预分割,以得到所述预分割后的图像;所述对所述预分割后的图像依次进行一系列处理的步骤包括:对所述预分割后的图像多次执行所述形状学处理。
其中,根据所述提取到的封闭区域的大小,确定所述目标区域的初始轮廓的步骤包括:对于所述提取到的封闭区域,保留面积最大的n个封闭区域,其中n为大于0的整数;根据该n个封闭区域的大小关系,确定所述目标区域的初始轮廓。
其中,所述根据该n个封闭区域的大小关系,确定所述目标区域的初始轮廓的步骤包括:当n=2时,若该2个封闭区域中,较大封闭区域的面积与较小封闭区域的面积之差小于预设阈值,则将该2个封闭区域同时确定为所述目标区域的初始轮廓;若所述较大封闭区域的面积与所述较小封闭区域的面积之差大于预设阈值,则将所述较大封闭区域确定为所述目标区域的初始轮廓。
其中,对所述目标区域的初始轮廓进行演化的步骤包括:采用基于区域的活动轮廓模型,来对所述目标区域的初始轮廓进行演化。
其中,在所述基于区域的活动轮廓模型中,通过边缘指示函数检测远离演化曲线的目标边缘,以引导所述演化曲线停止在目标轮廓的边界处,从而提高收敛速度;
其中边缘指示函数为:
其中,α∈(0,1),为图像边界强度场对演化速度的控制系数,β为比例常数,r为指数加权平均比率算子得到的边缘强度。
本发明实施例提供了一种超声图像的分析装置,用于从所述超声图像中分割出目标区域,所述分析装置包括:预处理模块,用于对所述超声图像进行预处理;初始轮廓提取模块,用于从预处理后的所述超声图像中提取出所述目标区域的初始轮廓;以及演化模块,用于对所述目标区域的初始轮廓进行演化,从而获取所述目标区域的准确边界。
本发明实施例的有益效果是:
本发明实施例,能够在保证分割结果准确性的同时,提高分割的自动化水平和速度。
附图说明
图1是是本发明的超声图像的分析方法的实施例的流程示意图;
图2a是图1中的步骤101的实施例的流程示意图;
图2b是原始超声图像和有效区域的实施例的示意图;
图3是图1中的步骤102的实施例的流程示意图;
图4a和图4b分别为传统cv模型和改进的cv模型的演化结果示意图;
图5是本发明的超声图像的分析装置的实施例的结构示意图。
具体实施方式
为了使本发明所解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
如图1所示,是本发明的超声图像的分析方法的实施例的流程示意图。其中该超声图像的分析方法可以集成于cad系统中,以用来从超声图像中分割出目标区域,从而便于cad系统提供辅助诊断数据。其中,该超声图像例如可以是乳腺肿瘤超声图像,但是本发明不限制于此。如图1所示,该超声图像的分析方法包括如下步骤:
步骤101:对超声图像进行预处理。其中,超声图像例如可以由超声设备产生,然后输入至cad系统进行处理。
步骤102:从预处理后的超声图像中提取出目标区域的初始轮廓。其中,目标区域例如可以是乳腺肿瘤超声图像中的肿瘤区域。
步骤103:对目标区域的初始轮廓进行演化,从而获取目标区域的准确边界。
本实施例,通过对超声图像进行预处理、目标区域的初始轮郭的提取,以及初始轮郭的演化,从而能够在保证分割结果准确性的同时,提高分割的自动化水平和速度。
下面将依次对图1的各个步骤进行说明。
如图2a所示,是图1中的步骤101的实施例的流程示意图。其包括如下步骤:
步骤201:选取所述超声图像的有效区域。
其中,在步骤201中,可以采用手动剪切法来选取有效区域,所谓手动剪切法即手动利用鼠标划取一个方框对原始图像进行剪裁,从而去掉周边的无关区域而保留图像的中心大区域。另外,也可以采取程序设置法来选取有效区域;因为同一型号的超声设备采集到的图像大小及边框大小是固定的,因此可以通过程序设置只保留图像中特定大小区域的像素,从而忽略掉周边的无关信息。
举例而言,如图2b所示,是原始的超声图像和经过步骤201处理后,得到的有效区域的图像;其中图2b中的左图为原始的超声图像,而右图为有效区域的图像。
步骤202:对选取的有效区域进行去噪处理。
其中,在步骤202中,可以采用srad(specklereducinganisotropicdiffusion,各向异性扩散)模型对图像进行滤波,以去除斑点噪声。其中,采用srad模型既能有效滤除超声图像斑点噪声、简化图像结构、提高图像质量,又能保持图像边缘处细节信息。
如图3所示,是图1中的步骤102的实施例的流程示意图,其包括如下步骤:
步骤301:对预处理后的超声图像进行预分割,得到预分割后的图像。
其中,在步骤301中,采用两次最大类间方差法来对预处理后的超声图像进行预分割。需要说明的是,最大类间方差法的详细内容对于本领域技术人员而言是熟悉的,因此不赘述。
具体而言,在步骤301中,首先采用最大类间方差法,来处理预处理后的超声图像,以得到第一阈值(即全局阈值)。接着,根据该第一阈值,将预处理后的超声图像划分为背景图像和前景图像,其中前景图像会包含目标区域。然后,继续采用最大类间方差法来处理前景图像,以得到第二阈值(最优阈值)。最后,利用该第二阈值来对预处理后的超声图像进行预分割,以得到预分割后的图像。需要说明的是,该预分割后的图像可以为一个二值化图像。
可选地,在步骤301中,可以仅执行一次最大类间方差法,即首先采用最大类间方差法,来处理预处理后的超声图像,以得到全局阈值。接着,利用该全局阈值来对预处理后的超声图像进行预分割,以得到预分割后的图像。
需要说明的是,在步骤301中,当采用两次最大类间方差法时,可以减小误分割。而当采用一次最大类间方差法时,为了提高精度,可以在后续的一系列处理中,增加形态学处理的次数。
步骤302:对预分割后的图像依次进行一系列处理,其中该一系列处理包括如下至少一项:形态学处理、空洞填充、以及去除与边界相连的区域。
其中,形态学处理可以是指:在经过预分割后得到的二值化图像中,往往含有一些因斑点噪声引起的噪声碎片。因此为了滤除图像中的小毛刺和孤立点、切断细长连接、平滑病灶(如肿瘤)边缘,对二值化图像(即预分割后的图像)进行形态学处理。同时,为了取得更好的效果,可以在腐蚀时采用较大的结构元素,膨胀时采用较小的结构元素。
其中,由于残余噪声、病灶内部灰度不均等影响,二值化后的病灶内部会产生“空洞”,故需对封闭区域内的小孔洞进行填充操作。
其中,二值化图像中与边界相连的区域往往是由伪影引起的误分割,不是真实的病灶区域,因此为了消除其影响,去除二值化图像中与边界相连的区域。
步骤303:提取经过一系列处理后的图像中的封闭区域。
其中,经过前述步骤的处理,图像中会出现一些封闭的连通区域(一般为黑色),这些即是肿瘤候选区域。在步骤中,将所有的封闭区域提取出来,并且可以按照各封闭区域的面积大小对其进行排序。
步骤304:根据提取到的封闭区域的大小,确定目标区域的初始轮廓。
其中,在步骤304中,对于提取到的封闭区域,保留面积最大的n个封闭区域,其中n为大于0的整数;根据该n个封闭区域的大小关系,确定目标区域的初始轮廓。需要说明的是,若仅有一个封闭区域,则直接将该封闭区域作为目标区域的初始轮廓。下面以n=2为例对该过程进行说明:
首先,对于提取到的封闭区域,提取面积最大的2个封闭区域,若较大封闭区域的面积与较小封闭区域的面积之差小于预设阈值,则将该2个封闭区域同时确定为目标区域的初始轮廓;若较大封闭区域的面积与较小封闭区域的面积之差大于预设阈值,则仅将较大封闭区域确定为目标区域的初始轮廓。其中,预设阈值可以为较大封闭区域的面积的一半或者其他数值。
本实施例,通过以上处理,可以自动提取到目标区域的初始轮廓,而无需手动设置初始轮廓,因此能够提高cad系统的自动化水平。
在图1中的步骤103中,主要采用基于区域的活动轮廓模型,来对目标区域的初始轮廓进行演化,从而获取目标区域的准确边界。其中,基于区域的活动轮廓模型也称为cv模型,是由chan和vese于2001年提出的一种基于图像全局区域信息的能量模型,其核心思想是利用分段常数去逼近图像的各个部分,通过变分法引入水平集建立方程,采用差分方法进行数值计算,最终将边界轮廓的演化问题转化为能量最小化问题。
而在本发明实施例中,并非直接采用cv模型,而是对cv模型进行了改进,以提高收敛速度。具体而言,在本发明实施例中,针对传统的cv模型,本发明实施例通过增加边缘指示函数(具体而言,是采用边缘指示函数来替换传统cv模中的dirac函数)来检测远离演化曲线的目标边缘,以引导演化曲线停止在目标轮廓的边界处,从而提高收敛速度。
其中边缘指示函数可以表示为:
其中,α∈(0,1),为图像边界强度场对演化速度的控制系数,β为比例常数,r为指数加权平均比率(theratioofexponentiallyweightedaverages,roewa)算子得到的边缘强度。此处,采用roewa算子而非一般图像分割算法采用的梯度值是考虑到超声图像中的斑点噪声为乘性噪声,而roewa算子能更好地适应噪声为乘性模型的图像。其中,所构建的边缘指示函数是一个单调递减函数,当r→∞时,g(r)→0,g的取值越趋向于0,表明该处越接近于病灶的真实边界轮廓。因此,边缘指示函数本身也可以加快cv模型的收敛速度。
另外,补充一点,经过以上改进,cv模型的能量函数为:
其中,u0≥0,ν≥0,λ1,λ2>0,且为给定的参数。
在图1的步骤103中,通过对传统的cv模型进行改进,可以提高cv模型的收敛速度,从而加快处理速度。例如,如下表所示,为在不同的迭代次数下,改进后的cv模型和传统cv模型的分割时间对比表。由该表可知,改进后的cv模型的分割时间显著地快于传统cv模型。
再例如,对于同一幅超声肿瘤图像。在迭代次数均为500次的情况下,利用传统cv模型进行演化的结果如图4(a)所示,而利用改进后的cv模型进行演化后的结果如图4(b)所示。由图可见,利用改进后的cv模型进行演化后的结果与肿瘤的真实边界更接近。
以上对本发明的实施例的图像分析方法进行了详细说明,下面说明相应地该方法的装置。需要说明的是,由于相关细节已在前面详述,因此下面仅说明装置的主要架构,而忽略其细节描述。
如图5所示,是超声图像的分析装置的实施例的结构示意图。该分析装置500可以用于从超声图像中分割出目标区域。其中,该分析装置500包括:预处理模块501,用于对所述超声图像进行预处理;初始轮廓提取模块502,用于从预处理后的所述超声图像中提取出所述目标区域的初始轮廓;以及演化模块503,用于对所述目标区域的初始轮廓进行演化,从而获取所述目标区域的准确边界。
本发明实施例的超声图像的分析装置,能够在保证分割结果准确性的同时,提高分割的自动化水平和速度。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(read-onlymemory,rom)或随机存储记忆体(randomaccessmemory,ram)等。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。