【摘要】 本发明涉及医学图像处理，旨在提供一种CTA图像中肝脏血管增强及肝脏与血管同时分割的方法。该方法包括：对图像进行高斯卷积预处理；计算图像每个点的各向异性特征值，进而确定各点的各向异性椭圆邻域；计算各领域内的灰度直方图，初始化一块肝脏区域，并计算其内的灰度直方图；计算各点邻域内灰度直方图与肝脏区域的灰度直方图的wasserstein距离；同时根据wasserstein距离和邻域各向异性特征增强血管；肝脏及血管分割。采用本发明提供的方法，能够克服低对比度、噪声、模糊边界等带来的负面影响，大大提高血管识别分割的准确率，从而准确的获取刚脏血管的解剖结构信息。。 【专利类型】发明授权 【申请人】浙江大学 【申请人类型】学校 【申请人地址】310027 浙江省杭州市西湖区浙大路38号 【申请人地区】中国 【申请人城市】杭州市 【申请人区县】西湖区 【申请号】CN201210014249.9 【申请日】2012-01-17 【申请年份】2012 【公开公告号】CN102609913B 【公开公告日】2014-07-09 【公开公告年份】2014 【授权公告号】CN102609913B 【授权公告日】2014-07-09 【授权公告年份】2014.0 【IPC分类号】G06T5/00 【发明人】孔德兴; 彭佳林; 吴法; 王金伟 【主权项内容】一种CTA图像中肝脏血管增强的方法，其特征在于，包括下述步骤：(1)对图像进行预处理，并交互地初始化；(2)确定图像I(x)每个像素点的各向异性特征及特征邻域；(3)确定图像I(x)每个像素点特征邻域内的灰度分布特征和纹理特征；(4)根据特征邻域内的灰度分布及纹理特征特征和各向异性特征增强血管；所述步骤(1)的具体过程是：假设输入图像I(x)是三维灰度图像数据，其中x=(x1，x2，x3)，定义的区域为对输入图像I(x)进行高斯平滑和中值滤波、去除噪声，将图像显示窗位调整至[0，N]；人工交互式地确定一块仅包含肝脏实质的区域以上所述的符号表示集合包含；所述步骤(2)的具体过程是：假设图像I(x)是三维图像数据，其中x=(x1，x2，x3)，定义的区域为a、计算每一个像素点x为中心，r为半径的球形邻域B(x，r)内灰度值分布的一阶矩M1(r, x)=∫ B(x, r)I(x~)(x~-M0(r, x))⊗ (x~-M0(r, x))dx~∫ B(x, r)I(x~)dx~]]>b、计算M1(r，x)的三个特征值分别是λ1，λ2，λ3，且假定λ1≥λ2≥λ3，相应的特征向量分别是v1，v2，v3；c、利用一阶矩的特征值确定图像每个像素点的各向异性特征AS：由于在血管的区域λ1＞＞λ2，λ3≈0，而在肝脏等区域λ1≈λ2≈λ3，从而用下面的量表征各向异性特征：AS=λ 2+λ 3λ 1+λ 2+λ 3]]>在血管内部AS≈0；d、在每一个像素点确定一个椭球形特征邻域D(x，a(x)，b(x)，b(x))，椭球半径长度分别为a2(x)=τ β , b2(x)=τ β , ]]>其中β =(1-32AS)× r+32AS, ]]>椭球的三条径的方向相应地由v1，v2，v3确定；以上所述的符号：为向量的张量积，∫B(x, r)表示在区域B(x, r)内积分，符号表示体积元，M0(r，x)表示零阶矩，即重心，符号≈表示近似等于，符号＞＞表示远远大于，τ为常数，r为半径常数，β为包含有特征值信息的半径长度；所述步骤(3)的具体过程是：假设图像I(x)是三维图像数据，x=(x1，x2，x3)，定义的区域为图像灰度范围是[0，N]；a、提取特征邻域内的灰度分布特征和纹理特征，计算每一个像素点x椭球形邻域D(x，a，b，b)中的累积灰度直方图Fx和初始化区域Ω0中累积灰度直方图b、计算任意点x累积灰度直方图Fx与累积直方图之间的Wasserstein距离W~(x)=W(Fx, FΩ 00)=∫ 0N 【当前权利人】浙江大学 【当前专利权人地址】浙江省杭州市西湖区浙大路38号 【专利权人类型】公立 【统一社会信用代码】12100000470095016Q 【引证次数】3.0 【被引证次数】4 【他引次数】3.0 【被自引次数】2.0 【被他引次数】2.0 【家族引证次数】3.0 【家族被引证次数】20

未经允许不得转载：http://www.zhongzhencnc.com/1791580391.html