技术领域本发明属于医学图像处理领域，涉及CTA序列图像中的肝脏血管增强、中心线提取和肝脏血管分割。
背景技术：
肝脏血管分割和三维重建有助于准确获取腹部肝脏血管组织的整体信息，是计算机辅助肝脏疾病诊断和肝脏手术规划的前提。CTA(计算机断层扫描血管造影)是一种无创伤的成像技术，具有密度分辨率高、对人体损伤小等优点，广泛应用于肝脏疾病的评估和诊断。由于肝脏血管结构复杂、血管之间相互缠绕、不同个体差异较大，肝脏血管分割面临巨大挑战。在临床应用中，为构建肝脏血管模型，放射科专家需根据血管的先验知识，对每一张CTA切片进行手动分割。由于肝脏CTA序列图像的切片数量很大(通常每位病人的CTA数据包含200～500张切片)，手动分割肝脏血管劳动强度大、耗时长、易产生错误分割结果，而且，分割的精度高度依赖于放射科专家的经验，具有很大的主观性。现有的三维肝脏血管分割方法一般可以分为基于灰度和梯度两大类。受血流速度、数据采样时间和血管狭窄处的影响，造影剂在血液中分布不均匀，导致在CTA图像中血管区域的灰度不一致。譬如，血管中心的灰度高于血管边界的灰度，主血管的灰度高于小血管的灰度。单一的基于灰度或梯度的分割方法，如三维区域生长、模糊聚类等，都无法有效地提取低对比度的肝脏血管。近年来，主动轮廓模型及其混合模型被广泛应用到三维血管分割中，但是这类模型的演化曲面容易越过血管的弱边界而导致严重的过分割，且需交互提供血管的初始区域。此外，以上方法对小血管的分割难以奏效。
技术实现要素：
针对现有方法的不足，本发明提出了一种肝脏CTA序列图像的血管分割方法，并且能够推广到其他三维血管的分割中。本发明的肝脏CTA序列图像的血管分割方法，包括以下步骤：步骤1：对原始三维肝脏图像进行预处理，提高血管的对比度，平滑CTA图像噪声。所述步骤1包括：首先，分别调节3D肝脏图像的窗宽/窗位到[400～800]/[200～500]范围，并将灰度值归一化到[0,255]，得到图像I0(x)，其中x为三维空间位置坐标。然后采用各向异性滤波算法去除噪声：∂I(x,t)∂t=div(c(x,t)·▿I)I(x,0)=I0(x)]]>其中，t为扩散时间，I(x,t)为t时刻的图像，div(·)为散度操作，为图像梯度，c(x,t)为非负的单调递减传导系数：c(x,t)=e-(|▿I|/k)2]]>其中，k为扩散门限，k为[50～80]之间的自然数。步骤2：采用多尺度OOF(OptimalOrientFlux)和OFA(OrientedFluxAnti-symmetry)算法增强血管及其相应的边界，细化血管中心。所述多尺度算法指针对肝脏主血管和小血管的半径(尺度)差异，选择不同半径的血管滤波器，提取滤波器的最大响应值得到最终结果。所述OOF滤波器指通过最小化球面通量f(x,r,ρ^)=14πr2∫∂Br((▿(Gσ0*I(x+h))·ρ^)ρ^)·n^dA=ρ^TQr,xρ^]]>寻找最佳的投影方向，构建血管增强滤波器，其中，*为卷积运算，Br表示以x为球心，r为半径的球表面，为球表面的外向单位法向量，为球表面上的位置坐标，dA为球表面无穷小的面积分量，是方差为σ0的高斯函数，I(x)为经各向异性滤波后得到的图像，为方向向量，Qr,x为3×3通量矩阵，其第i行、第j列元素为qr,xi,j=14πr2∂i,j(Gσ0(x)*br(x)*I(x))]]>其中，i,j∈{1,2,3