专利名称：Cta图像中肝脏血管增强及肝脏与血管同时分割的方法
技术领域：
本发明涉及医学图像处理领域，具体的说是CT (计算机断层扫描造影增强)图像中肝脏血管的增强以及肝实体与血管的同时分割。
背景技术：
肝脏及血管结构对于肝病诊断及肝病手术的指导具有非常重要的意义。CTA(计算机断层扫描造影增强)技术是最常用的肝脏及血管影像诊断技术之一。通过CTA，医生可以得到一系列血管增强的二维CT切片。肝脏体积是肝病手术的重要参考指标，然而腹部 CTA图像往往存在低对比度，高信噪比，边界模糊，肝脏与其他灰度相近组织粘连等不利因素，使得准确的肝脏分割成为一个难题。另外肝脏血管结构本身很复杂，医生很难直接通过二维切片推知完整的三维血管结构。有效地增强及准确地分割CTA图像中的血管已经成为计算机辅助手术的关键技术。由于在CT图像扫描之前检查者被注射了造影剂，CTA图像中血管的灰度值相对于肝组织会偏高。但是由于造影剂在血液中的不均匀分布及血管狭窄等因素造成血管出现明显的灰度不均匀以及血管与肝组织之间弱对比度等现象，这往往使得基于灰度值和基于血管边界强度的图像分割算法失效。人们提出了一系列的活动轮廓模型解决此类问题如CURVES，但由于血管结构复杂，在很多地方都非常细，而且由于前面提到灰度不均匀和低对比度问题，这类模型很难准确地分割血管，尤其是细的血管。其次，血管分叉的地方往往也是很多根据各向异性特征方法不能很好地探测的地方。此外由于腹部扫描图像中含有很多无关的组织，如胃，心脏，脾脏等等，与肝脏灰度值很接近，其中也会含有血管，这使得在复杂的图像中直接提取血管及分割出肝脏有很大困难。
发明内容
本发明要解决的技术问题是，克服现有技术中的不足，提供一种能够克服灰度不均、低对比度、噪声弱边界等不利影响的血管图像增强和分割的方法，同是提供了一种能够克服弱边界，高信噪比，轻度病变和灰度值相近的软组织干扰的肝脏分割方法。为解决该技术问题，本发明采用的技术方案如下本发明公开了一种CTA图像中肝脏血管增强的方法，包括下述步骤(I)对图像进行去噪预处理，并交互地初始化；(2)确定图像I (X)每个像素点的各向异性特征及特征邻域；(3)确定图像I (X)每个像素点特征邻域内的灰度分布特征和纹理特征；(4)根据特征邻域内的灰度分布及纹理特征特征和各向异性特征增强血管。本发明中，所述肝脏血管增强方法中步骤(I)的具体过程是首先对输入图像I(X)进行高斯平滑和中值滤波、去除噪声，将图像显示窗位调整至
，然后人工交互式地确定一块仅包含肝脏实质的区域c Ω Ci 3。以上所述符号Ω为图像I(X)定义域，符号c表示集合包含于。本发明中，所述肝脏血管增强方法中步骤(2)的具体过程是
a.计算每一个像素点χ半径为r的球形邻域B(x，r)内灰度值分布的一阶矩
权利要求
1.一种CTA图像中肝脏血管增强的方法，其特征在于，包括下述步骤(1)对图像进行预处理，并交互地初始化；(2)确定图像I(X)每个像素点的各向异性特征及特征邻域；(3)确定图像I(X)每个像素点特征邻域内的灰度分布特征和纹理特征；(4)根据特征邻域内的灰度分布及纹理特征特征和各向异性特征增强血管；所述步骤(I)的具体过程是假设输入图像I(x)是三维灰度图像数据，其中X= (Xl，x2，x3)，定义的区域为Qci 3; 对输入图像I(x)进行高斯平滑和中值滤波、去除噪声，将图像显示窗位调整至
； 人工交互式地确定一块仅包含肝脏实质的区域[ Ω ；以上所述的符号C表示集合包含；所述步骤(2)的具体过程是假设图像I(X)是三维图像数据，其中X= (Xl，χ2，χ3)，定义的区域为Qci 3;a、计算每一个像素点X为中心，r为半径的球形邻域B(x,r)内灰度值分布的一阶矩
2.—种CTA图像中肝脏与血管同时分割的方法，其特征在于，包括下述步骤(1)模型初始化；(2)优化新的变分能量模型，同时分割肝脏和血管；(3)进行后处理，分别得到肝脏轮廓及血管轮廓；所述步骤(I)具体过程是设输入图像/(X)是三维图像数据，X = (X1, x2, X3)，定义的区域为^，图像灰度范围是
;在肝脏区域内部初始化一块区域Qtl，其边界曲面记为Ctl;令尸O) =其中『
全文摘要
本发明涉及医学图像处理，旨在提供一种CTA图像中肝脏血管增强及肝脏与血管同时分割的方法。该方法包括对图像进行高斯卷积预处理；计算图像每个点的各向异性特征值，进而确定各点的各向异性椭圆邻域；计算各领域内的灰度直方图，初始化一块肝脏区域，并计算其内的灰度直方图；计算各点邻域内灰度直方图与肝脏区域的灰度直方图的wasserstein距离；同时根据wasserstein距离和邻域各向异性特征增强血管；肝脏及血管分割。采用本发明提供的方法，能够克服低对比度、噪声、模糊边界等带来的负面影响，大大提高血管识别分割的准确率，从而准确的获取刚脏血管的解剖结构信息。
文档编号G06T5/00GK102609913SQ201210014249
公开日2012年7月25日 申请日期2012年1月17日 优先权日2012年1月17日
发明者吴法, 孔德兴, 彭佳林, 王金伟 申请人:浙江大学