基于概率分布和连续最大流的主动脉瓣超声图像分割方法技术

本发明专利技术涉及一种基于概率分布和连续最大流的主动脉瓣超声图像分割方法，包括以下几个步骤：第一步，采集人体主动脉瓣短轴的医学超声图像数据，等间隔地提取五帧先验图像；第二步，分割五帧先验图像；第三步，构建二维灰度-距离直方图；第四步，由二维灰度-距离直方图计算获得综合概率估计函数；第五步，分别计算其各自单独的概率估计函数；第六步，对五帧先验图像，分别计算其各自最能代表前景和背景的像素灰度值；第七步，对于当前待分割图像，求取其独立的概率估计图；第八步，对上述前景区域和五帧先验图像的手动分割结果分别进行相似度测量；第九步，获得分割结果。与现有技术相比，本发明专利技术具有稳定可靠、实现方便、适用于实际临床等优点。

【技术实现步骤摘要】

本专利技术涉及一种图像处理
，尤其是涉及一种。
技术介绍
医学超声成像由于其特有的优点，如实时性、无损伤性、可重复性、灵敏度高、成本低，而受到广泛的关注。在基于医学超声图像的主动脉瓣的临床诊断和治疗中，根据病人不同的病理情况，需要提取主动脉瓣的信息，其重要的一项手段就是图像分割。图像分割结果的好坏直接影响病变组织结构的定位、定量、定性分析，三维重建等后续操作以及图像引导手术、肿瘤放射治疗等所采取治疗规划的准确性。但是，由于主动脉瓣的医学超声图像存在以下特征，使得其分割成为一个困难棘手的问题1)超声图像分辨率、对比度低；2)超声图像固有的斑点噪声的影响；3)超声图像回波纹理特性；4)主动脉瓣开放和闭合时造成的三 个瓣叶的移动；5)瓣叶和瓣环严重钙化造成的伪影。在传统的医学图像分割领域里，基于像素灰度的不连续性和相似性，其可以分为基于边缘检测的分割方法和基于区域的分割方法。前者在处理超声图像获取边缘信息时容易产生间断或者伪边缘，不能形成理想的闭合曲线，需要进一步处理以获得目标轮廓。后者中基于区域生长的方法，其分割结果与种子点的选择有很大关系，同时对噪声很敏感，容易形成孔状或不连续的区域。近年来，随着新技术的发展，医学超声图像也相继出现了新的方法，如动态规划模型、主动轮廓模型、水平集方法等。动态规划模型是为解决多节段决策问题最优化而提出的。它是将边缘检测问题看成一个优化问题，通过求解全局最优解获得最终的分割结果。但其运算量较大，同时因噪声影响会产生较多的伪边缘，从而难以得到较好的结果。主动轮廓模型即snake模型，通过构建合适的内部能量函数和外部能量函数来约束目标轮廓线的运动。其缺点是分割结果对初始曲线的位置和形状较为敏感，难以分割凹陷区域处的目标，同时容易陷入局部最优解，而无法获得理想的目标轮廓。水平集方法的核心是把η维曲线描述为(η+1)维的水平集，主要理论基础是曲线演化理论和零水平集的思想。大部分水平集模型都是通过求解偏微分方程进行数值实现的，这就造成了一个计算量较大的问题，另外难以确定合适的速度函数。目前，在临床诊断和治疗中，应用较为广泛的分割方法是基于阈值的分割方法和基于医生手动分割的方法。阈值分割方法原理简单、操作方便，但是由于上述的超声图像中固有的斑点噪声和纹理特性的影响，其分割结果十分不佳，难以为临床治疗提供精确有效的指导信息。基于医生手动分割的方法，由于医生具有长期的临床实践经验，所分割出来的结果是可以进一步用于临床诊断治疗的，但因其工作量大、耗费时间长，也不具有较好的应用价值
技术实现思路
本专利技术的目的就是为了克服上述现有技术存在的缺陷而提供一种稳定可靠，且实现方便，适用于实际临床的。本专利技术的目的可以通过以下技术方案来实现一种，其特征在于，包括以下几个步骤第一步，采集人体主动脉瓣短轴的医学超声图像数据，选择一个连续完整的心动周期，等间隔地提取五帧先验图像，其中每帧先验图像会代表一个心动周期内的不同相位；第二步，手动分割上述五帧先验图像，并计算每帧分割结果的包围盒，取其中最大的包围盒为后续过程使用；第三步，根据先验图像分割结果各自的独立中心点，计算得到先验图像的一个综 合中心点，以像素点的灰度值为横轴、以像素点与综合中心点的距离为纵轴，由上述五帧先验图像的分割结果构建一个二维灰度-距离直方图；第四步，由上述二维灰度-距离直方图计算获得综合概率估计函数，即计算每个像素点在一定的灰度以及与综合中心点一定的距离条件下属于前景的条件概率，拥有高灰度值并且距离综合中心点较近的像素点属于前景的概率高，反之，拥有高灰度值但距离综合中心点较远的像素点属于背景区域的概率高；第五步，对上述五帧先验图像，分别计算其各自单独的概率估计函数，即计算每个像素点在一定的灰度以及与自身独立中心点一定的距离条件下属于前景的条件概率，同样，拥有高灰度值并且距离综合中心点较近的像素点属于前景的概率高，反之，拥有高灰度值但距离综合中心点较远的像素点属于背景区域的概率高；第六步，根据上一步已获得的单独概率估计函数，对五帧先验图像，分别计算其各自最能代表前景和背景的像素灰度值；第七步，对于当前待分割图像，利用综合概率估计函数求取其独立的概率估计图，并采用大津算法对概率估计图进行预处理，保留前景区域；第八步，对上述前景区域和五帧先验图像的手动分割结果分别进行相似度测量，找到先验图像中与当前待分割图像相似度测量值最高的图像，记该图像为Ism ;第九步，根据Ism在第六步中已经获得的最能代表前景和背景的像素灰度值构建能量方程，确立能量方程中的数据项和平滑项，采用连续最大流算法进行图像分割，获得主动脉瓣短轴在一个心动周期内所有帧图像的分割结果。所述的第四步中的综合概率估计函数为P(x) = P(F(x) I R(x)J(x))=- Μ 尸(Λ W JW)其中，P(Fi)是像素点属于前景的概率，P(R，I)是灰度值为I、与综合中心点距离为R的像素点属于前景的概率，η是先验图像的个数，此处η = 5，F代表图像前景区域。所述的第五步中的单独的概率估计函数为 . T, P(FnR(x),I(x))Pi(X)=P(FiiX) I R(x)J(x)) = ~F(KiX)1I(X))其中，i = 1，2，3，4，5，分别代表上述五帧先验图像。所述的第六步中的最能代表各先验图像的前景和背景的像素灰度值4 权利要求1.一种，其特征在于，包括以下几个步骤 第一步，采集人体主动脉瓣短轴的医学超声图像数据，选择一个连续完整的心动周期，等间隔地提取五帧先验图像，其中每帧先验图像会代表一个心动周期内的不同相位； 第二步，手动分割上述五帧先验图像，并计算每帧分割结果的包围盒，取其中最大的包围盒为后续过程使用； 第三步，根据先验图像分割结果各自的独立中心点，计算得到先验图像的一个综合中心点，以像素点的灰度值为横轴、以像素点与综合中心点的距离为纵轴，由上述五帧先验图像的分割结果构建一个二维灰度-距离直方图； 第四步，由上述二维灰度-距离直方图计算获得综合概率估计函数，即计算每个像素点在一定的灰度以及与综合中心点一定的距离条件下属于前景的条件概率，拥有高灰度值并且距离综合中心点较近的像素点属于前景的概率高，反之，拥有高灰度值但距离综合中心点较远的像素点属于背景区域的概率高； 第五步，对上述五帧先验图像，分别计算其各自单独的概率估计函数，即计算每个像素点在一定的灰度以及与自身独立中心点一定的距离条件下属于前景的条件概率，同样，拥有高灰度值并且距离综合中心点较近的像素点属于前景的概率高，反之，拥有高灰度值但距离综合中心点较远的像素点属于背景区域的概率高； 第六步，根据上一步已获得的单独概率估计函数，对五帧先验图像，分别计算其各自最能代表前景和背景的像素灰度值； 第七步，对于当前待分割图像，利用综合概率估计函数求取其独立的概率估计图，并采用大津算法对概率估计图进行预处理，保留前景区域； 第八步，对上述前景区域和五帧先验图像的手动分割结果分别进行相似度测量，找到先验图像中与当前待分割图像相似度测量值最高的图像，记该图像为Ism ; 第九步，根据Ism在第六步中已经获得的最能代表前景和背景的像素灰度值构建能量方程，确立能量方程中的数据项和平滑项，采用连本文档来自技高网...

【技术保护点】
一种基于概率分布和连续最大流的主动脉瓣超声图像分割方法，其特征在于，包括以下几个步骤：第一步，采集人体主动脉瓣短轴的医学超声图像数据，选择一个连续完整的心动周期，等间隔地提取五帧先验图像，其中每帧先验图像会代表一个心动周期内的不同相位；第二步，手动分割上述五帧先验图像，并计算每帧分割结果的包围盒，取其中最大的包围盒为后续过程使用；第三步，根据先验图像分割结果各自的独立中心点，计算得到先验图像的一个综合中心点，以像素点的灰度值为横轴、以像素点与综合中心点的距离为纵轴，由上述五帧先验图像的分割结果构建一个二维灰度？距离直方图；第四步，由上述二维灰度？距离直方图计算获得综合概率估计函数，即计算每个像素点在一定的灰度以及与综合中心点一定的距离条件下属于前景的条件概率，拥有高灰度值并且距离综合中心点较近的像素点属于前景的概率高，反之，拥有高灰度值但距离综合中心点较远的像素点属于背景区域的概率高；第五步，对上述五帧先验图像，分别计算其各自单独的概率估计函数，即计算每个像素点在一定的灰度以及与自身独立中心点一定的距离条件下属于前景的条件概率，同样，拥有高灰度值并且距离综合中心点较近的像素点属于前景的概率高，反之，拥有高灰度值但距离综合中心点较远的像素点属于背景区域的概率高；第六步，根据上一步已获得的单独概率估计函数，对五帧先验图像，分别计算其各自最能代表前景和背景的像素灰度值；第七步，对于当前待分割图像，利用综合概率估计函数求取其独立的概率估计图，并采用大津算法对概率估计图进行预处理，保留前景区域；第八步，对上述前景区域和五帧先验图像的手动分割结果分别进行相似度测量，找到先验图像中与当前待分割图像相似度测量值最高的图像，记该图像为ISM；第九步，根据ISM在第六步中已经获得的最能代表前景和背景的像素灰度值构建能量方程，确立能量方程中的数据项和平滑项，采用连续最大流算法进行图像分割，获得主动脉瓣短轴在一个心动周期内所有帧图像的分割结果。...

【技术特征摘要】

【专利技术属性】
技术研发人员：顾力栩，聂媛媛，罗哲，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：