CN104978730B - 一种左心室心肌的分割方法和装置 - Google Patents
一种左心室心肌的分割方法和装置 Download PDFInfo
- Publication number
- CN104978730B CN104978730B CN201410141320.9A CN201410141320A CN104978730B CN 104978730 B CN104978730 B CN 104978730B CN 201410141320 A CN201410141320 A CN 201410141320A CN 104978730 B CN104978730 B CN 104978730B
- Authority
- CN
- China
- Prior art keywords
- left ventricle
- image
- central point
- blood pool
- myocardium
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
本发明提供了一种左心室心肌的分割方法和装置，所述方法包括如下步骤：(1)定位左心室中心点，根据所述左心室中心点，获取所述左心室中心点所在层图像中的感兴趣区域；(2)根据灰度信息和所述左心室中心点，在所述左心室中心点所在层图像中的感兴趣区域中定位出血池区域；(3)确定所述血池区域的质心，并根据所述血池位置质心，确定原始图像中除了所述左心室中心点所在层图像之外的所有层图像的血池区域，即获得左心室内膜的分割结果；(4)根据所述左心室内膜分割结果，通过动态规划的方法，得到左心室外膜的分割结果，即得到所述左心室心肌的分割结果。本技术方案能完全自动、均匀、稳定地分割左心室心肌。
Description
技术领域
本发明涉及图像处理领域，尤其涉及一种左心室心肌的分割方法和装置。
背景技术
左心室心肌分割是心功能分析的主要步骤，左心室心肌壁厚、射血量等心功能分析参数均基于左心室心肌的准确分割。左心室心肌和背景灰度类似，左心室内乳头肌、小梁等和心肌属于同一物质，还存在伪影。因此，左心室心肌分割是一个困难而又重要的任务。
目前，已有不少左心室内、外膜的分割方法。左心室内膜的分割方法有基于活动轮廓形状模型和基于动态规划的方法等。基于活动轮廓形状模型的方法，利用一些边界特征点和均匀间隔点来表征左心室内膜的形状，但是，该方法存在计算复杂等问题。基于动态规划的方法，利用左心室内膜的梯度和边界点之间的距离约束关系，来搜索左心室内膜边界，但是，该方法需要大致确定内膜边界的位置。
左心室外膜的分割方法有基于极坐标区域增长的方法和基于动态规划的方法等。基于极坐标区域增长的方法，将直角坐标图像转换为极坐标图像，然后利用区域增长的方式获得二值图像，最后利用傅里叶变换来平滑曲线，该方法最大的问题在于区域增长容易出现大量过分割，而滤波并不能平滑曲线。基于动态规划的方法，利用左心室外膜边界的梯度和边界之间的距离约束关系，来搜索左心室外膜边界，但是，该方法不能有效地抑制非外膜边界梯度，容易造成分割错误。
发明内容
本发明解决的问题是提供一种左心室心肌的分割方法和装置，能完全自动、均匀、稳定地分割左心室心肌。
为了解决上述问题，本发明提供了一种左心室心肌的分割方法，包括如下步骤：
(1)定位左心室中心点，根据所述左心室中心点，获取所述左心室中心点所在层图像中的感兴趣区域；
(2)根据灰度信息和所述左心室中心点，在所述左心室中心点所在层图像中的感兴趣区域中定位出血池区域；
(3)确定所述血池区域的质心，并根据所述血池位置质心，确定原始图像中除了所述左心室中心点所在层图像之外的所有层图像的血池区域，即获得左心室内膜的分割结果；
(4)根据所述左心室内膜分割结果，通过动态规划的方法，得到左心室外膜的分割结果，即得到所述左心室心肌的分割结果。
上述一种左心室心肌的分割方法，优选地，所述定位左心室中心点的过程包括：
定位图像中的异地相位图像中的中间层图像；
利用所述霍夫变换，检测所述中间层图像中的最大圆形区域；
以所述最大圆形区域的中心点作为所述左心室中心点。
上述一种左心室心肌的分割方法，优选地，所述获得左心室内膜的分割结果的过程包括：
(1)基于所述血池区域的质心，得到所述异地相位图像中除了所述中间层图像的各层图像的血池区域，即得到所述异地相位图像中的血池区域；
(2)在所述图像中，所述异地相位图像两侧分别按离异地相位图像由近及远的次序依次选择相位图像作为当前相位；
(3)根据所述当前相位前一相位或后一相位中左心室中心点，获取所述当前相位的中间层图像的左心室中心点，并根据所述当前相位的中间层图像的左心室中心点，定位所述当前相位中间层图像中的血池区域；
(4)根据所述当前相位中间层图像中血池区域的质心，获得所述当前相位图像的血池区域；
(5)对所述图像中各相位图像中得到的血池区域进行凸化处理以及滤波处理，得到所述左心室内膜分割结果。
上述一种左心室心肌的分割方法，优选地，得到所述异地相位图像中除了所述中间层图像的各层图像的血池区域的过程包括：
(1)在所述异地相位图像中，所述中间层图像两侧分别按离中间层图像由近及远的次序依次选择图像层作为当前层；
(2)根据所述当前层前一层或后一层中血池区域的质心，获取所述当前层的中心点，根据所述当前层中心点，获取所述当前层的感兴趣区域；
(3)根据灰度信息，将所述当前层感兴趣区域划分为三个区域，并根据所述当前层中心点，定位所述三个区域中的血池区域，即所述当前层图像中的血池区域。
上述一种左心室心肌的分割方法，优选地，得到所述左心室外膜分割结果的过程包括：
(1)根据所述左心室内膜分割结果和所述图像的信息，确定目标区域，其中，所述目标域包含所述左心室的外膜；
(2)抑制所述目标区域中的非心肌区域的灰度，得到抑制图像；
(3)在所述抑制图像中，利用动态规划的方法，得到左心室外膜的初次分割结果；
(4)对所述左心室外膜的初次分割结果进行修正，得到所述左心室外膜分割结果。
上述一种左心室心肌的分割方法，优选地，得到所述抑制图像的过程包括：
(1)利用聚类方法，将所述目标区域分为三类区域，对所述三类区域的灰度值居中的区域，进行二值化，得到心肌范围掩膜；
(2)对所述心肌范围掩膜内所有像素点的灰度就行统计，得到所述心肌范围掩膜的灰度分布参数；
(3)根据所述灰度分布参数，以所述心肌范围掩膜中像素点的灰度范围为轴，将所述第二感兴趣区域中灰度大于预定阈值的像素点的灰度进行映射。
上述一种左心室心肌的分割方法，优选地，将所述第二感兴趣区域中灰度大于预定阈值的像素点的灰度进行映射的过程包括：以所述像素点灰度均值为中心点，将所述第二感兴趣区域中灰度大于所述灰度范围最大值的像素点P，映射到小于所述灰度范围最小灰度值的范围内，其余像素点的灰度不变。
上述一种左心室心肌的分割方法，优选地，将所述第二感兴趣区域中灰度大于所述灰度范围最大值的像素点P映射到小于所述灰度范围最小灰度值的范围内的公式为：P′i=Rmin-(Pi-Rmax)，其中，P′i为所述像素点P映射后的灰度值；Rmin为所述灰度范围中的最小灰度值；Pi为所述像素点P映射前的灰度值；Rmax为所述灰度范围中的最大灰度值。
上述一种左心室心肌的分割方法，优选地，对所述左心室外膜的初次分割结果进行修正的过程包括：将所述初次分割结果中灰度值在所述灰度范围之外的像素点的灰度值赋值为0，得到第一修正结果；对所述第一修正结果进行凸化和傅里叶滤波，得到第二修正结果；在时间序列上，对所述第二修正结果进行均值滤波。
为了解决上述问题，本发明还提供了一种左心室心肌的分割装置，包括：
中心点定位单元，适于定位左心室中心点，根据所述左心室中心点，获取所述左心室中心点所在层图像中的感兴趣区域；
血池区域获取单元，适于根据灰度信息和所述左心室中心点，在所述左心室中心点所在层图像中的感兴趣区域中定位出血池区域；
左心室内膜分割单元，适于确定所述血池区域的质心，并根据所述血池位置质心，确定原始图像中除了所述左心室中心点所在层图像之外的所有层图像的血池区域，即获得左心室内膜的分割结果；
左心室外膜分割单元，适于根据所述左心室内膜分割结果，通过动态规划的方法，得到所述左心室的外膜的分割结果，即得到所述左心室心肌的分割结果。
与现有技术相比，本发明利用心脏的解剖结构属性，通过血池来确定左心室内膜，并根据左心室内膜，通过动态规划方法，得到左心室外膜，能自动、均匀、稳定地分割出左心室心肌，且将乳头肌、小梁等排除；
进一步地，只对异地相位图像的中间层图像进行霍夫变换，加快计算速度。
附图说明
图1所示为本发明实施例一种左心室心肌的分割方法的流程示意图；
图2所示为本发明实施例获得左心室内膜的分割结果的优选过程的流程示意图；
图3所示为本发明实施例得到所述异地相位图像中除了所述中间层图像的各层图像的血池区域的优选过程的流程示意图；
图4所示为本发明实施例得到所述左心室外膜分割结果的优选过程的流程示意图；
图5所示为本发明实施例一种左心室心肌的分割装置结构示意图。
具体实施方式
在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广，因此本发明不受下面公开的具体实施的限制。
其次，本发明利用示意图进行详细描述，在详述本发明实施例时，为便于说明，所述示意图只是实例，其在此不应限制本发明保护的范围。
下面结合附图和实施例对本发明一种左心室心肌的分割方法做出详细地说明。本发明实施例左心室心肌的分割方法如图1所示，首先，执行步骤S1，定位左心室中心点，根据所述左心室中心点，获取所述左心室中心点所在层图像中的感兴趣区域。其中，所述定位左心室中心点的过程包括：定位图像中的异地相位图像中的中间层图像；利用所述霍夫变换，检测所述中间层图像中的最大圆形区域；以所述最大圆形区域的中心点作为所述左心室中心点具体地，在本实施例中，以磁共振(MR)心脏图像为例，由于在采集图像的时候心脏跳动，MR对心脏连续检测，得到一系列数据，将这些采集到的数据重建得到10个相位图像，每一个相位图像中包含有20帧图像。优选地，由于在心脏跳动过程中，在舒张末期的时候，左心室容量最大，血液最多，因此，选取该时期的相位图像，即异地相位图像。需要说明的是，也可以为计算机断层扫描(CT)图像等其它影像图像，相位图像也可以为10个相位图像中其它的任一相位图像。
在异地相位图像中，选取中间层图像，即第10帧图像，在该层图像中，利用霍夫变换，检测该中间层图像中的最大圆形区域，所述最大圆形区域的中心为所述左心室中心点。并以该中心点为初始点，根据经验确定其周围的一定区域为感兴趣区域。需要说明的是，也可以通过区域增长等方法来确定感兴趣区域。
接着，执行步骤S2，根据灰度信息和所述左心室中心点，在所述左心室中心点所在层图像中的感兴趣区域中定位出血池区域。具体地，在本实施例中，根据异地相位图像中间层图像的灰度信息，通过Kmeans聚类方法对步骤S1中得到的感兴趣区域进行分类，得到三个区域，而所述左心室中心点所在的区域为血池区域，即所述左心室中心点所在层图像中的左心室内膜区域。需要说明的是，也可以将步骤S1中得到的感兴趣区域分成两个区域，只要将血池区域划分出来即可。
接着，执行步骤S3，确定所述血池区域的质心，并根据所述血池位置质心，确定原始图像中除了所述左心室中心点所在层图像之外的所有层图像的血池区域，即获得左心室内膜的分割结果。
接着，执行步骤S4，根据所述左心室内膜分割结果，通过动态规划的方法，得到左心室外膜的分割结果，即得到所述左心室心肌的分割结果。
其中，所述步骤S3中获得左心室内膜的分割结果的优选过程如图2所示，首先，执行步骤S201，基于所述血池区域的质心，得到所述异地相位图像中除了所述中间层图像的各层图像的血池区域，即得到所述异地相位图像中的血池区域。
具体地，在本实施例中，得到所述异地相位图像中除了所述中间层图像的各层图像的血池区域的优选过程如图3所示，首先，执行步骤S301，在所述异地相位图像中，所述中间层图像两侧分别按离中间层图像由近及远的次序依次选择图像层作为当前层。具体地，在本实施例中，在异地相位图像中，以第11帧图像为当前层。接着，执行步骤S302，根据所述当前层前一层或后一层中血池区域的质心，获取所述当前层的中心点，根据所述当前层中心点，获取所述当前层的感兴趣区域。具体地，在本实施例中，步骤S2中得到了中间层图像中的血池区域(即第10帧图像的血池区域)，从而可以得到该血池区域的质心。由于一个完整的位相图像包含整个心跳周期，因此，在时间序列上图像存在关联性。根据图像的这一个特性，第10帧图像血池区域质心所在位置，在第11帧图像中该位置所对应的位置即为第11帧图像的中心点。再以该中心点为初始点，根据经验确定其周围的一定区域为感兴趣区域。接着，执行步骤S303，根据灰度信息，将所述当前层感兴趣区域划分为三个区域，并根据所述当前层中心点，定位所述三个区域中的血池区域，即所述当前层图像中的血池区域。具体地，在本实施例中，根据第11帧图像的灰度信息，通过Kmeans聚类方法对步骤S302中得到的感兴趣区域进行分类，得到三个区域，而所述第11帧图像中心点所在的区域为血池区域，即所述第11帧图像中的左心室内膜区域。同样地，以第12帧图像为当前层图像，采用和第11帧图像相同的处理方法，可以得到第12帧图像中的左心室内膜区域。依此类推，可以得到异地相位图像中所有20帧图像的血池区域，所述20帧图像的血池区域连接即为异地相位图像的血池区域，即异地相位图像的左心室内膜区域。
接着，执行步骤S202，在所述图像中，所述异地相位图像两侧分别按离异地相位图像由近及远的次序依次选择相位图像作为当前相位。在本实施例中，所述异地相位图像为所述MR图像中10个相位图像中的第4个相位图像，在异地相位图像两侧，按照离异地相位图像由近及远的次序，选取第5个相位图像为当前相位图像。接着，执行步骤上203，根据所述当前相位前一相位或后一相位中左心室中心点，获取所述当前相位的中间层图像的左心室中心点，并根据所述当前相位的中间层图像的左心室中心点，定位所述当前相位中间层图像中的血池区域。具体地，根据心脏图像时间上序列上的关联性，步骤S1中得到的异地相位图像的左心室中心点的位置在第5个相位图像中的对应的位置，即为所述第5个相位图像的左心室中心点，也为第5个相位图像中间层图像(第10帧图像)的左心室中心点。采用和异地相位图像中间层图像同样的方法，即通过步骤S2的方法，得到第5个相位图像中间层图像(第10帧图像)的血池区域。
接着，执行步骤S204，根据所述当前相位中间层图像中血池区域的质心，获得所述当前相位图像的血池区域。具体地，采用和异地相位图像中其它层图像同样的方法，通过步骤S301、S302和S303可以得到第5个相位图像的血池区域。依此类推，可以得到所述MR图像中所有十个相位图像的血池区域。
接着，执行步骤S205，对所述图像中各相位图像中得到的血池区域进行凸化处理以及滤波处理，得到所述左心室内膜分割结果。具体地，将步骤S204中得到的十个相位图像的血池区域进行连接，得到所述MR图像中的血池区域，即左心室内膜分割的初步结果，对所述左心室分割的初步结果进行凸化处理，排除乳头肌和小梁等的影响，接着进行傅里叶滤波，从而得到MR图像的左心室内膜分割结果。
其中，所述步骤S4中得到所述左心室外膜分割结果的优选过程如图4所示，首先，执行步骤S401，根据所述左心室内膜分割结果和所述图像的信息，确定目标区域，其中，所述目标包含所述左心室的外膜。具体地，根据上述得到的左心室内膜分割结果，以及MR图像中的灰度信息，根据经验在所述MR图像中确定第二感兴趣区域。接着，执行步骤S402，抑制所述目标中的非心肌区域的灰度，得到抑制图像。具体地，在本实施例中，利用Kmeans聚类方法，将步骤S401中得到的目标分为三类区域，根据所述三类区域的灰度值从大到小排列，第二类区域即为包含有心肌的区域，接着根据第二类区域灰度的范围直接对第二类区域进行二值化，得到心肌范围掩膜。假设所述心肌范围掩膜中的像素点灰度符合高斯分布，那么，将所述心肌范围掩膜所在的图像区域作为统计对象，计算该区域内所有像素点的均值和方差，即得到了所述心肌范围掩膜的灰度分布参数。
接着，根据所述灰度分布参数，以所述心肌范围掩膜中像素点的灰度范围为轴，将所述第二感兴趣区域中灰度大于预定阈值的像素点的灰度进行映射。具体地，以所述像素点灰度均值为中心点，将所述第二感兴趣区域中灰度大于所述灰度范围最大值的像素点，通过公式(1)映射到小于所述灰度范围最小灰度值的范围内，其余像素点的灰度不变。其中，所述公式(1)为：
P′i=Rmin-(Pi-Rmax)，其中，P′i为所述像素点P映射后的灰度值；Rmin为所述灰度范围中的最小灰度值；Pi为所述像素点P映射前的灰度值；Rmax为所述灰度范围中的最大灰度值。通过上述步骤，得到所述抑制图像。
接着，执行步骤S403，在所述抑制图像中，利用动态规划的方法，得到左心室外膜的初次分割结果。具体地，以所述左心室内膜的边界为起始边界，根据经验值向外取一定距离的区域，得到环形区域，通过极坐标转换将所述环形区域转换到极坐标图像，根据动态规划公式，其中，所述公式中含有极坐标中每个像素点的内能(基于梯度)和外能(相邻点之间的距离约束)项，在所述极坐标图像中，以每一列为基准，计算并存储相对于每一个像素点的最小能量，并向后追溯最小能量路径，最后，将最小路径像素点坐标转换到直角坐标系中，得到左心室外膜的边界，即所述左心室外膜的初次分割结果。
接着，执行步骤S404，对所述左心室外膜的初次分割结果进行修正，得到所述左心室外膜分割结果。具体地，将所述初次分割结果中灰度值在所述心肌范围掩膜的灰度范围之外的像素点的灰度值赋值为0，得到第一修正结果；对所述第一修正结果进行凸化和傅里叶滤波，得到第二修正结果；接着，在时间序列上，对所述第二修正结果进行均值滤波。从而得到左心室外膜的分割结果。分割出左心室内膜和外膜，即分割出左心室心肌。
对应于上述左心室心肌的分割方法，本发明实施例还提供了一种左心室心肌的分割装置，如图5所示，所述分割装置包括中心点定位单元1、血池区域获取单元2、左心室内膜分割单元3和左心室外膜分割单元4。
所述中心点定位单元1适于定位左心室中心点，根据所述左心室中心点，获取所述左心室中心点所在层图像中的感兴趣区域。其中，定位图像中的异地相位图像中的中间层图像，利用所述霍夫变换，检测所述中间层图像中的最大圆形区域，所述最大圆形区域的中心为所述左心室中心点。
所述血池区域获取单元2适于根据灰度信息和所述左心室中心点，在所述左心室中心点所在层图像中的感兴趣区域中定位出血池区域。具体地，根据异地相位图像中间层图像的灰度信息，通过Kmeans聚类方法对步骤S1中得到的感兴趣区域进行分类，得到三个区域，而所述左心室中心点所在的区域为血池区域，即所述左心室中心点所在层图像中的左心室内膜区域。
所述左心室内膜分割单元3适于确定所述血池区域的质心，并根据所述血池位置质心，确定原始图像中除了所述左心室中心点所在层图像之外的所有层图像的血池区域，即获得左心室内膜的分割结果。具体地，所述左心室内膜分割单元3包括异地相位图像血池区域获取单元、当前相位确定单元、当前相位中间层图像血池区域确定单元、当前相位图像获取单元和血池区域处理单元。所述异地相位图像血池区域获取单元适于基于所述血池区域的质心，得到所述异地相位图像中除了所述中间层图像的各层图像的血池区域，即得到所述异地相位图像中的血池区域；当前相位确定单元适于在所述图像中，所述异地相位图像两侧分别按离异地相位图像由近及远的次序依次选择相位图像作为当前相位；当前相位中间层图像血池区域确定单元适于根据所述当前相位前一相位或后一相位中左心室中心点，获取所述当前相位的中间层图像的左心室中心点，并根据所述当前相位的中间层图像的左心室中心点，定位所述当前相位中间层图像中的血池区域；当前相位图像获取单元适于根据所述当前相位中间层图像中血池区域的质心，获得所述当前相位图像的血池区域；血池区域处理单元适于对所述图像中各相位图像中得到的血池区域进行凸化处理以及滤波处理，得到所述左心室内膜分割结果。
所述左心室外膜分割单元4适于据所述左心室内膜分割结果，通过动态规划的方法，得到所述左心室的外膜的分割结果，即得到所述左心室心肌的分割结果。具体地，所述左心室外膜分割单元4包括目标区域确定单元、抑制图像获取单元、动态规划单元以及修正单元。所述目标区域确定单元适于根据所述左心室内膜分割结果和所述图像的信息，确定目标区域，其中，所述目标域包含所述左心室的外膜；抑制图像获取单元适于抑制所述目标区域中的非心肌区域的灰度，得到抑制图像；动态规划单元适于在所述抑制图像中，利用动态规划的方法，得到左心室外膜的初次分割结果；修正单元适于对所述左心室外膜的初次分割结果进行修正，得到所述左心室外膜分割结果。
上述左心室心肌分割装置中各单元的配合及工作过程可以参考上述左心室心肌分割方法的说明，在此不再赘述。
本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。
Claims (10)
1.一种左心室心肌的分割方法，其特征在于，包括如下步骤：
(1)定位左心室中心点，根据所述左心室中心点，获取所述左心室中心点所在层图像中的感兴趣区域；
(2)根据灰度信息和所述左心室中心点，在所述左心室中心点所在层图像中的感兴趣区域中定位出血池区域；
(3)确定所述血池区域的质心，并根据所述血池区域的质心，确定原始图像中除了所述左心室中心点所在层图像之外的所有层图像的血池区域，即获得左心室内膜的分割结果；
(4)根据所述左心室内膜分割结果，通过动态规划的方法，得到左心室外膜的分割结果，即得到所述左心室心肌的分割结果。
2.如权利要求1所述一种左心室心肌的分割方法，其特征在于，所述定位左心室中心点的过程包括：
定位图像中的异地相位图像中的中间层图像；
利用霍夫变换，检测所述中间层图像中的最大圆形区域；
以所述最大圆形区域的中心点作为所述左心室中心点。
3.如权利要求2所述一种左心室心肌的分割方法，其特征在于，所述获得左心室内膜的分割结果的过程包括：
(1)基于所述血池区域的质心，得到所述异地相位图像中除了所述中间层图像的各层图像的血池区域，即得到所述异地相位图像中的血池区域；
(2)在所述图像中，所述异地相位图像两侧分别按离异地相位图像由近及远的次序依次选择相位图像作为当前相位；
(3)根据所述当前相位前一相位或后一相位中左心室中心点，获取所述当前相位的中间层图像的左心室中心点，并根据所述当前相位的中间层图像的左心室中心点，定位所述当前相位中间层图像中的血池区域；
(4)根据所述当前相位中间层图像中血池区域的质心，获得所述当前相位图像的血池区域；
(5)对所述图像中各相位图像中得到的血池区域进行凸化处理以及滤波处理，得到所述左心室内膜分割结果。
4.如权利要求3所述一种左心室心肌的分割方法，其特征在于，得到所述异地相位图像中除了所述中间层图像的各层图像的血池区域的过程包括：
(1)在所述异地相位图像中，所述中间层图像两侧分别按离中间层图像由近及远的次序依次选择图像层作为当前层；
(2)根据所述当前层前一层或后一层中血池区域的质心，获取所述当前层的中心点，根据所述当前层中心点，获取所述当前层的感兴趣区域；
(3)根据灰度信息，将所述当前层感兴趣区域划分为三个区域,并根据所述当前层中心点，定位所述三个区域中的血池区域,即所述当前层图像中的血池区域。
5.如权利要求1所述一种左心室心肌的分割方法，其特征在于，得到所述左心室外膜分割结果的过程包括：
(1)根据所述左心室内膜分割结果和所述图像的信息，确定目标区域，其中，所述目标区域包含所述左心室的外膜；
(2)抑制所述目标区域中的非心肌区域的灰度，得到抑制图像；
(3)在所述抑制图像中，利用动态规划的方法，得到左心室外膜的初次分割结果；
(4)对所述左心室外膜的初次分割结果进行修正，得到所述左心室外膜分割结果。
6.如权利要求5所述一种左心室心肌的分割方法，其特征在于，得到所述抑制图像的过程包括：
(1)利用聚类方法，将所述目标区域分为三类区域，对所述三类区域的灰度值居中的区域，进行二值化，得到心肌范围掩膜；
(2)对所述心肌范围掩膜内所有像素点的灰度就行统计，得到所述心肌范围掩膜的灰度分布参数；
(3)根据所述灰度分布参数，以所述心肌范围掩膜中像素点的灰度范围为轴，将第二感兴趣区域中灰度大于预定阈值的像素点的灰度进行映射。
7.如权利要求6所述一种左心室心肌的分割方法，其特征在于，将所述第二感兴趣区域中灰度大于预定阈值的像素点的灰度进行映射的过程包括：以所述像素点灰度均值为中心点，将所述第二感兴趣区域中灰度大于所述灰度范围最大值的像素点P，映射到小于所述灰度范围最小灰度值的范围内，其余像素点的灰度不变。
8.如权利要求7所述一种左心室心肌的分割方法，其特征在于，将所述第二感兴趣区域中灰度大于所述灰度范围最大值的像素点P映射到小于所述灰度范围最小灰度值的范围内的公式为：Pi′＝Rmin-(Pi-Rmax)，其中，Pi′为所述像素点P映射后的灰度值；Rmin为所述灰度范围中的最小灰度值；Pi为所述像素点P映射前的灰度值；Rmax为所述灰度范围中的最大灰度值。
9.如权利要求6所述一种左心室心肌的分割方法，其特征在于，对所述左心室外膜的初次分割结果进行修正的过程包括：将所述初次分割结果中灰度值在所述灰度范围之外的像素点的灰度值赋值为0，得到第一修正结果；对所述第一修正结果进行凸化和傅里叶滤波，得到第二修正结果；在时间序列上，对所述第二修正结果进行均值滤波。
10.一种左心室心肌的分割装置，其特征在于，包括：
中心点定位单元，适于定位左心室中心点,根据所述左心室中心点,获取所述左心室中心点所在层图像中的感兴趣区域；
血池区域获取单元，适于根据灰度信息和所述左心室中心点，在所述左心室中心点所在层图像中的感兴趣区域中定位出血池区域；
左心室内膜分割单元，适于确定所述血池区域的质心,并根据所述血池区域的质心，确定原始图像中除了所述左心室中心点所在层图像之外的所有层图像的血池区域，即获得左心室内膜的分割结果；
左心室外膜分割单元，适于根据所述左心室内膜分割结果,通过动态规划的方法,得到所述左心室的外膜的分割结果，即得到所述左心室心肌的分割结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410141320.9A CN104978730B (zh) | 2014-04-10 | 2014-04-10 | 一种左心室心肌的分割方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410141320.9A CN104978730B (zh) | 2014-04-10 | 2014-04-10 | 一种左心室心肌的分割方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104978730A CN104978730A (zh) | 2015-10-14 |
| CN104978730B true CN104978730B (zh) | 2019-04-23 |
Family
ID=54275211
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410141320.9A Active CN104978730B (zh) | 2014-04-10 | 2014-04-10 | 一种左心室心肌的分割方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104978730B (zh) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106611411B (zh) * | 2015-10-19 | 2020-06-26 | 上海联影医疗科技有限公司 | 一种医学图像中肋骨分割的方法及医学图像处理装置 |
| CN106709919B (zh) * | 2015-11-17 | 2020-08-04 | 上海联影医疗科技有限公司 | 一种心脏灌注磁共振图像的处理方法 |
| US10346993B2 (en) | 2015-11-17 | 2019-07-09 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for image processing in magnetic resonance imaging |
| CN106910194B (zh) * | 2015-12-22 | 2020-11-27 | 上海联影医疗科技股份有限公司 | 心脏心功能磁共振图像中的心内外膜的分割方法 |
| CN106910182B (zh) * | 2015-12-22 | 2020-11-27 | 上海联影医疗科技股份有限公司 | 心脏心功能磁共振图像中舒张末期图像中血池的分割方法 |
| US10290109B2 (en) | 2015-12-22 | 2019-05-14 | Shanghai United Imaging Healthcare Co., Ltd. | Method and system for cardiac image segmentation |
| CN110163876B (zh) * | 2019-05-24 | 2021-08-17 | 山东师范大学 | 基于多特征融合的左心室分割方法、系统、设备及介质 |
| CN111028254B (zh) * | 2019-12-09 | 2023-03-24 | 上海交通大学医学院附属上海儿童医学中心 | 一种基于心脏电影磁共振图像的左心肌分割方法 |
| CN111210423B (zh) * | 2020-01-13 | 2023-04-07 | 浙江杜比医疗科技有限公司 | Nir图像的乳房轮廓提取方法、系统及装置 |
| CN111904450B (zh) * | 2020-09-07 | 2023-11-07 | 北京永新医疗设备有限公司 | 左心室中心和感兴趣区域的提取方法、装置及系统 |
| CN112932535B (zh) * | 2021-02-01 | 2022-10-18 | 杜国庆 | 一种医学图像分割及检测方法 |
| CN113160116B (zh) * | 2021-02-03 | 2022-12-27 | 中南民族大学 | 左心室内外膜自动分割方法、系统及设备 |
| CN113469948B (zh) * | 2021-06-08 | 2022-02-25 | 北京安德医智科技有限公司 | 左心室节段识别方法及装置、电子设备和存储介质 |
| CN113362345B (zh) * | 2021-06-30 | 2023-05-30 | 武汉中科医疗科技工业技术研究院有限公司 | 图像分割方法、装置、计算机设备和存储介质 |
| CN114419032B (zh) * | 2022-03-14 | 2022-06-21 | 深圳科亚医疗科技有限公司 | 心脏左心室的心肌内膜和/或心肌外膜的分割方法和装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101596111A (zh) * | 2008-05-23 | 2009-12-09 | 西门子公司 | 在心脏电影磁共振成像中对左心室的自动定位 |
| WO2010027798A2 (en) * | 2008-08-26 | 2010-03-11 | Northwestern University | Ablation devices and related methods thereof |
-
2014
- 2014-04-10 CN CN201410141320.9A patent/CN104978730B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101596111A (zh) * | 2008-05-23 | 2009-12-09 | 西门子公司 | 在心脏电影磁共振成像中对左心室的自动定位 |
| WO2010027798A2 (en) * | 2008-08-26 | 2010-03-11 | Northwestern University | Ablation devices and related methods thereof |
Non-Patent Citations (1)
| Title |
|---|
| Segmentation of Left Ventricle in Cardiac Cine MRI: An Automatic Image-Driven Method;YingLi Lu et al.;《 Functional Imaging and Modeling of the Heart: Lecture Notes in Computer Science Volume 5528. Berlin》;20091231;339-347 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104978730A (zh) | 2015-10-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104978730B (zh) | 一种左心室心肌的分割方法和装置 | |
| CN107203741B (zh) | 血管提取方法、装置及其系统 | |
| CN105913432B (zh) | 基于ct序列图像的主动脉提取方法及装置 | |
| US9715637B2 (en) | Method and system for automatic aorta segmentation | |
| US7715609B2 (en) | Method for automatically determining the position and orientation of the left ventricle in 3D image data records of the heart | |
| CN110245721A (zh) | 神经网络模型的训练方法、装置和电子设备 | |
| CN105096310B (zh) | 利用多通道特征在磁共振图像中分割肝脏的方法和系统 | |
| CN101331517B (zh) | 利用可变形格网的帮助图像后处理的方法 | |
| US11793572B2 (en) | Selecting a medical device for use in a medical procedure | |
| CN103810363B (zh) | 血管造影术中的血管种子点选取方法及血管提取方法 | |
| CN104851108B (zh) | 基于ct图像的肝动脉分割方法 | |
| CN102890823A (zh) | 运动对象轮廓提取及左心室图像分割方法和装置 | |
| CN107545579B (zh) | 一种心脏分割方法、设备和存储介质 | |
| CN102968783B (zh) | 从腹部图像数据中自动分割骨骼的方法及系统 | |
| CN102068281A (zh) | 一种占位性病变超声图像的处理方法 | |
| CN106910194B (zh) | 心脏心功能磁共振图像中的心内外膜的分割方法 | |
| CN109886938A (zh) | 一种超声图像血管直径自动测量方法 | |
| DE102017203313A1 (de) | Verfahren zur Auswertung eines kontrastmittelverstärkten Magnetresonanzschichtbildes eines Herzens, Bildverarbeitungseinrichtung, Computerprogramm und elektronisch lesbarer Datenträger | |
| CN106133789A (zh) | 用于分割感兴趣区域的图像处理装置和方法 | |
| US8050470B2 (en) | Branch extension method for airway segmentation | |
| CN104915989B (zh) | 基于ct影像的血管三维分割方法 | |
| CN104463862A (zh) | 肾脏ct序列图像快速分割方法 | |
| CN109166127A (zh) | 一种可穿戴式植物表型感知系统 | |
| CN109102509A (zh) | 分割模型训练方法、装置及计算机可读存储介质 | |
| CN112037146A (zh) | 医学图像伪影自动校正方法、装置和计算机设备 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201815 Shanghai City, Jiading District Jiading Industrial Zone No. 1180 Building 8 Hing Road Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201815 Shanghai City, Jiading District Jiading Industrial Zone No. 1180 Building 8 Hing Road Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder | ||
| CP02 | Change in the address of a patent holder |
Address after: 201807 2258 Chengbei Road, Jiading District, Shanghai Patentee after: Shanghai Lianying Medical Technology Co.,Ltd. Address before: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai Patentee before: Shanghai Lianying Medical Technology Co.,Ltd. |
|
| CP02 | Change in the address of a patent holder |