CN112842384B - 超声心动图心肌包络量测方法、装置和存储介质 - Google Patents
超声心动图心肌包络量测方法、装置和存储介质 Download PDFInfo
- Publication number
- CN112842384B CN112842384B CN202011631107.8A CN202011631107A CN112842384B CN 112842384 B CN112842384 B CN 112842384B CN 202011631107 A CN202011631107 A CN 202011631107A CN 112842384 B CN112842384 B CN 112842384B
- Authority
- CN
- China
- Prior art keywords
- myocardial
- control point
- tracking
- control points
- envelope
- 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
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
本申请涉及超声图像处理技术领域,具体涉及一种超声心动图包络量测方法、装置和存储介质。其中,方法包括:获取超声心动图,超声心动图包括M帧超声心脏切面图像;确定超声心动图中第N帧超声心脏切面图像,N∈M;将在第N帧超声心脏切面图像中,包围特定心腔的心肌,划分为Z个心肌节段;依次确定各个心肌节段的追踪控点;根据追踪控点,确定在第N帧超声心脏切面图像中,相邻心肌节段之间的轮廓包络;装置,包括存储器和处理器,所述存储器中存储有至少一条程序指令,所述处理器通过加载并执行所述至少一条程序指令以实现上述方法。存储介质中存储有至少一条程序指令,所述至少一条程序指令被处理器加载并执行以实现上述方法。
Description
技术领域
本申请涉及超声图像处理技术领域,具体涉及一种超声心动图包络量测方法、装置和存储介质。
背景技术
心肌在心脏的心动周期中会发生收缩和舒张的形变,在心动周期中,对心肌形变状态进行准确、客观、定量的评价,有利于对各种心脏疾病在临床阶段早期症状的判定。
超声心动图是诊断复杂先天性心脏病的可靠医学依据,对明确诊断、评估病情以及预估选择手术方案具有重要的临床价值。根据超声心动图各帧的超声图像中,心肌的变化来确定心肌应变量或心肌应变率等参数,定量心肌形变,从而定量心脏室壁运动,且结果不受心脏整体旋转等运动的影响。
相关技术中,通过截取心肌兴趣区,测算兴趣区中的心肌应变状态,以反映特定心腔的应变状态。但是此种技术手段对于心肌厚度不均匀的患者,不能真实反映其特定心腔的心肌应变。
发明内容
本申请提供了一种超声心动图心肌包络量测方法、装置和存储介质,可以解决相关技术中不能真实反映其特定心腔心肌应变的问题。
为了解决上述技术问题,本申请的第一方面,提供一种超声心动图心肌包络量测方法,所述超声心动图心肌包络量测方法包括以下步骤:
获取超声心动图,所述超声心动图包括M帧超声心脏切面图像,所述M={1,2,3…m},m为大于1的整数;
确定所述超声心动图中第N帧超声心脏切面图像,所述N∈M;
将在所述第N帧超声心脏切面图像中,包围特定心腔的心肌,划分为Z个心肌节段,所述Z={1,2,3…z},z为大于1的整数;
依次确定各个所述心肌节段的追踪控点;所述追踪控点的控点类型包括:用于追踪心肌内膜的内膜控点、用于追踪心肌外膜的外膜控点,和,用于追踪心肌中膜的中膜控点;
根据所述追踪控点,确定在所述第N帧超声心脏切面图像中,相邻心肌节段之间的轮廓包络;所述轮廓包络包括:根据所述内膜控点确定的内膜包络、根据所述外膜控点确定的外膜包络,和,根据所述中膜控点确定的中膜包络。
可选的,所述将在所述第N帧超声心脏切面图像中,包围特定心腔的心肌,划分为Z个心肌节段的步骤,包括:
在所述第N帧超声心脏切面图像中,沿着所述心肌的周向,依次将包围所述特定心腔的心肌,划分为Z个心肌节段。
可选的,所述依次确定各个所述心肌节段的追踪控点的步骤中,确定第P个心肌节段的追踪控点,所述P∈Z,包括以下步骤:
获取第P个心肌节段所在图像区域各像素点的灰度值;
根据所述图像区域各像素点的灰度值,确定所述第P个心肌节段的心肌切面边缘;所述心肌切面边缘包括相对的内边缘和外边缘;
在所述内边缘处,标注至少三个追踪像素点作为内膜控点,并提取所述内膜控点的坐标;
在所述外边缘处,标注至少三个追踪像素点作为外膜控点,并提取所述外膜控点的坐标;
在所述内边缘和所述外边缘之间,标注至少三个追踪像素点作为中膜控点,并提取所述中膜控点的坐标。
可选的,确定所有所述心肌节段的追踪控点的步骤完成后,还进行:
根据所述追踪控点的坐标,确定控点集合;所述控点集合包括各个所述追踪控点的坐标信息、对应的心肌节段信息,以及控点类型信息;
判断所述控点集合中,各个所述追踪控点是否存在重复控点;
确定一追踪控点存在重复控点,获取所述重复控点对应的心肌节段和控点类型;
在与所述重复控点同一心肌节段,重新确定同一控点类型的追踪控点,并提取所述追踪控点的坐标;直至所述控点集合中不存在重复控点。
可选的,所述超声心动图中的M帧超声心脏切面图像,为按同一切面方向获取的心脏位置超声图像。
可选的,所述切面方向包括短轴切面。
可选的,在所述根据所述追踪控点,确定在所述第N帧超声心脏切面图像中,相邻心肌节段之间的轮廓包络的步骤进行前,在所述依次确定各个所述心肌节段的追踪控点的步骤完成后,还进行:
确定控点集合,所述控点集合为所有追踪控点的集合,包括起始控点和结束控点;
使得所述控点集合中不存在重复控点;
判断所述起始控点和结束控点之间的距离是否大于预设阈值;
确定所述起始控点和结束控点之间的距离大于预设阈值,在所述控点集合起始控点的周围新插入一追踪控点;
将所述新插入的所述追踪控点更新为结束控点。
本申请的第一方面,还提供一种超声心动图心肌包络量测方法,所述超声心动图心肌包络量测方法,还包括以下步骤:
根据上述超声心动图心肌包络量测方法,依次确定所述超声心动图各帧的超声心脏切面图像中,相邻心肌节段之间的轮廓包络;各条所述轮廓包络对应所述超声心动图的各帧,分别为第M轮廓包络,所述M={2,3…m},m为大于1的整数;
根据相邻X帧所述超声心脏切面图像之间对应轮廓包络的变化,计算所述特定心腔的心肌应变,所述X∈M。
本申请的第二方面,一种超声心动图心肌包络量测装置,所述超声心动图心肌包络量测装置包括存储器和处理器,所述存储器中存储有至少一条程序指令,所述处理器通过加载并执行所述至少一条程序指令以实现本申请第一方面所述的超声心动图心肌包络量测方法。
本申请的第三方面,提供一种计算机存储介质,所述计算机存储介质中存储有至少一条程序指令,所述至少一条程序指令被处理器加载并执行以实现本申请第一方面所述的方法。
本申请技术方案,至少包括如下优点:本申请通过对超声心动图中各帧超声心脏切面图像,分别进行特定心腔心肌节段划分,从而确定个节段中的追踪控点,根据追踪控点确定相邻心肌节段之间的轮廓包络,从而为真实反映其特定心腔心肌应变提供依据。
附图说明
为了更清楚地说明本申请具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本申请的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出了本申请一实施例提供的超声心动图心肌包络量测方法流程图;
图1a示出了在特定帧上,心尖四腔切面超声图像示意图;
图1b示出了按照左心室短轴切面,在特定帧上的心脏超声图像示意图;
图2a示出了超声获取的心脏房室瓣打开时,心尖四腔切面的示意图像;
图2b示出了超声获取左心室短轴切面的示意图像;
图3a示出了对图2a中各个心肌节段确定追踪控点后的示意图;
图3b示出了对图2b中各个心肌节段确定追踪控点后的示意图;
图4示出了追踪控点去重处理流程;
图5示出了左心室切面超声图像的轮廓包络调整流程;
图6示出了本申请其他实施例提供的超声心动图心肌包络量测方法流程图;
图7示出了本申请一实施例提供的超声心动图心肌包络量测装置。
具体实施方式
下面将结合附图,对本申请中的技术方案进行清楚、完整的描述,显然,所描述的实施例是本申请的一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在不做出创造性劳动的前提下所获得的所有其它实施例,都属于本申请保护的范围。
在本申请的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本申请的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电气连接;可以是直接相连,也可以通过中间媒介间接相连,还可以是两个元件内部的连通,可以是无线连接,也可以是有线连接。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本申请中的具体含义。
此外,下面所描述的本申请不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
由于心脏的左心室收缩功能是测定超声心动图检查的重要组成部分,能够用于指导制定治疗策略。例如,能够通过左心室的整体收缩功能判断疾病进展和好转情况,通过左心室局部收缩功能分析诊断冠心病和进行心脏负荷试验。因此以下实施例主要针对超声心动图中左心室心肌,阐述其轮廓包络的量测方法。
图1示出了本申请一实施例提供的超声心动图心肌包络量测方法流程图,参照图1,该超声心动图心肌包络量测方法至少包括以下步骤:
步骤S11:获取超声心动图,所述超声心动图包括M帧超声心脏切面图像,所述M={2,3…m},m为大于1的整数。
为了保证超声心动图能够更有效地反应心脏的收缩状态,优选地,该M大于8。
另外,超声声束具有一定的方向性,本实施例中该超声心动图中的M帧超声心脏切面图像,为按同一切面方向获取的心脏位置超声图像。其中该切面方向可以为心尖四腔切面,即该M帧超声心脏切面图像为M帧心尖四腔切面超声图像;该切面方向还可以为短轴切面,即M帧超声心脏切面图像为M帧短轴切面超声图像。
步骤S12:确定所述超声心动图中第N帧超声心脏切面图像,所述N∈M。
对于切面方向为心尖四腔切面,该心尖四腔切面为从左侧胸壁心尖到心底横截左、右心房和心室的平面。图1a示出了在特定帧上,心尖四腔切面超声图像示意图,参照图1a,该心尖四腔切面超声图像显示心脏的四个心腔,图1a右上方的心腔为左心室LV,呈椭圆形;左上方为右心室RV,呈三角形;右下方为左心房LA,左下方为右心房RA,心房与心室之间间隔有能够打开的房室瓣,两个心室之间间隔有室间隔,两个心房之间间隔有房间隔。
对于切面方向为短轴切面,该短轴切面超声图像的获取方式为在标准左心室长轴切面的基础上,将超声探头顺时针旋转90°,与心脏左心室的长轴垂直,在此位置发射超声声束,根据超声回波确定心脏的短轴切面超声图像。可以通过调整超声探测器的倾斜角度获取左心室短轴切面,图1b示意出了按照左心室短轴切面,在特定帧上的心脏超声图像示意图,参照图1b可以看出,该左心室短轴切面超声图像包括两个心腔,位于下侧的左心室LV和位于上侧的右心室RV。
步骤S13:将在所述第N帧超声心脏切面图像中,包围左心室的心肌,划分为Z个心肌节段,所述Z={1,2,3…z},z为大于1的整数。
图2a与图1a相对应,示出了超声获取的心脏房室瓣打开时,心尖四腔切面的示意图像,参照图2a,沿着左心室LV心肌的周向,依次将包围该左心室LV的心肌,划分为七个心肌节段。该七个心肌节段按图2a顺时针方向依次为间隔基底段211、间隔中间段212、间隔心尖段213、心尖段214、侧壁心尖段215、侧壁中间段216和侧壁基底段217。
图2b与图1b相对应,示出了超声获取左心室短轴切面的示意图像,参照图2b,沿着左心室LV心肌的周向,依次将包围该左心室LV的心肌,划分为六个心肌节段。其中,位于下侧的左心室LV和位于上侧的右心室RV之间的室间隔221划分为两个节段,左心室LV外周的左心室侧壁222划分为四个节段,且图2b所示的左心室LV六个心肌节段形成封闭的环形。
步骤S14:依次确定各个所述心肌节段的追踪控点;所述追踪控点的控点类型包括:用于追踪心肌内膜的内膜控点、用于追踪心肌外膜的外膜控点,和,用于追踪心肌中膜的中膜控点。
本实施例,在依次确定各个所述心肌节段的追踪控点的步骤中,确定第P个心肌节段的追踪控点,所述P∈Z,包括以下步骤:
步骤S141:获取第P个心肌节段所在图像区域各像素点的灰度值。
步骤S142:根据所述图像区域各像素点的灰度值,确定所述第P个心肌节段的心肌切面边缘;所述心肌切面边缘包括相对的内边缘和外边缘。
可以根据各个像素点的灰度值,通过边缘检测技术确定特定心肌切面超声图像中心肌的内边缘和外边缘,从图2a或图2b可以看出,左心室LV的心肌包括内边缘和外边缘,左心室LV的心肌内边缘为心肌内膜,左心室LV的心肌外边缘为心肌外膜,在左心室LV心肌内边缘和外边缘之间定义出心肌中膜,优选地,该心肌内膜位于特定切面方向心肌剖面的中间。
步骤S143:在所述内边缘处,标注至少三个追踪像素点作为内膜控点,并提取所述内膜控点的坐标。
步骤S144:在所述外边缘处,标注至少三个追踪像素点作为外膜控点,并提取所述外膜控点的坐标。
步骤S145:在所述内边缘和所述外边缘之间,标注至少三个追踪像素点作为中膜控点,并提取所述中膜控点的坐标。
参照图3a,其示出了对图2a中各个心肌节段确定追踪控点后的示意图。图3a中每个心肌节段的各个膜层分别确定了三个追踪控点,以心肌节段侧壁基底段217为例,沿着该节段的心肌外膜依次确定了三个外膜控点201,沿着该节段的心肌内膜依次确定了三个内膜控点202,沿着该节段的心肌中膜确定了三个中膜控点203。其他心肌节段也同样确定各自的追踪控点,并存储各个追踪控点的控点类型信息、坐标信息和所处心肌节段信息,形成控点集合。
参照图3b,其示出了对图2b中各个心肌节段确定追踪控点后的示意图。图3b中每个心肌节段的各个膜层分别确定了三个追踪控点,对于一特定心肌节段来说,沿着该节段的心肌外膜依次确定了三个外膜控点201,沿着该节段的心肌内膜依次确定了三个内膜控点202,沿着该节段的心肌中膜确定了三个中膜控点203。其他心肌节段也同样确定各自的追踪控点,并存储各个追踪控点的控点类型信息、坐标信息和所处心肌节段信息,形成控点集合。
步骤S15:根据所述追踪控点,确定在所述第N帧超声心脏切面图像中,相邻心肌节段之间的轮廓包络;所述轮廓包络包括:根据所述内膜控点确定的内膜包络、根据所述外膜控点确定的外膜包络,和,根据所述中膜控点确定的中膜包络。
继续参照图3a,依次连接同一控点类型相邻两个追踪控点,从而确定与该种控点类型对应的轮廓包络。即按图3a所示顺时针或逆时针依次连接外膜控点201,确定外膜包络A,依次连接内膜控点202,确定内膜包络B,依次连接中膜控点203,确定中膜包络C。
继续参照图3b,依次连接同一控点类型相邻两个追踪控点,从而确定与该种控点类型对应的轮廓包络。即按图3b所示顺时针或逆时针依次连接外膜控点201,确定外膜包络A,依次连接内膜控点202,确定内膜包络B,依次连接中膜控点203,确定中膜包络C。由于图3b所示的左心室LV心肌为封闭环形,因此对于左心室短轴切面超声图像中的左心室心肌轮廓包络应当为封闭曲线。
在其他实施例中,在步骤S14完成后,在步骤S15进行前,还需对所有追踪控点进行去重处理,以去除重复的追踪控点,即步骤S14完成后,在步骤S15进行前还可以进行步骤:确定控点集合,所述控点集合为所有追踪控点的集合,包括起始控点和结束控点;使得所述控点集合中不存在重复控点:
作为一种实施例可以采用图4所示的追踪控点去重处理流程,依次进行的以下步骤,使得所述控点集合中不存在重复控点:
步骤S41:根据所述追踪控点的坐标,确定控点集合;所述控点集合包括各个所述追踪控点的坐标信息、对应的心肌节段信息,以及控点类型信息。
步骤S42:判断所述控点集合中,各个所述追踪控点是否存在重复控点。本实施例中,当至少两个追踪控点的坐标信息相同,即其中一个追踪控点存在其他重复控点。
步骤S43:确定一追踪控点存在重复控点,获取所述重复控点对应的心肌节段和控点类型。
步骤S44:在与所述重复控点同一心肌节段,重新确定同一控点类型的追踪控点,并提取所述追踪控点的坐标;直至所述控点集合中不存在重复控点。
对于左心室切面超声图像的轮廓包络测量过程中,通常会出现轮廓包络不封闭的问题,为了避免上述问题,在确定心肌轮廓包络前,还可参照图5所示左心室切面超声图像的轮廓包络调整流程,进行以下步骤:
步骤S51:确定控点集合,所述控点集合为所有追踪控点的集合,包括起始控点和结束控点。
步骤S52:使得所述控点集合中不存在重复控点。
其中,步骤S51和步骤S52可以采用图4所示的追踪控点去重处理。
步骤S53:判断所述起始控点和结束控点之间的距离是否大于预设阈值。
步骤S54:确定所述起始控点和结束控点之间的距离大于预设阈值,在所述控点集合起始控点的周围新插入一追踪控点。
可以采用B-样条插值技术,在控点集合起始控点的周围新插入一追踪控点。
步骤S55:将所述新插入的所述追踪控点更新为结束控点。
若确定起始控点和结束控点之间的距离大于预设阈值,一方面会使得在进行步骤S15后所形成的轮廓包络未封闭,另一方面还可能影响所形成的轮廓包络的准确性。本实施例通过在起始控点的周围新插入一追踪控点,能够缩短轮廓包络中起始控点和结束控点之间的距离,使得最终确定的心肌轮廓包络封闭,能够提高轮廓包络的准确性。
在其他实施例中,对于各个追踪控点可通过选中、拖拽改变其位置,并在控点集合中对应更新其信息。即可以先,选中特定追踪控点;然后,确定所选中追踪控点的坐标信息、控点类型信息和所处心肌节段信息;再,拖拽所选中的追踪控点移动至特定位置;最后,更新该追踪控点移动至特定位置的坐标信息、控制类型信息和所处心肌节段信息。
需要注意的是对于内膜包络,其拖拽放大后的位置不可超过心肌外膜边缘,对于外膜包络,对其拖拽缩小后的位置不可超过心肌内膜的边缘。且拖拽后还需再次进行图4所示的追踪控点去重处理。
在其他实施例中,对于各条轮廓包络,可通过选中、拖拽缩放其大小,对应改变确定该轮廓包络的追踪控点的位置,并在控点集合中对应更新该条轮廓包络上各个追踪的信息。即可以先选中特定轮廓包络;然后,确定该条轮廓包络上的追踪控点,以及各个追踪控点的坐标信息、控点类型信息和所处心肌节段信息;再,拖拽所选中的轮廓包络进行等比例缩小或放大,在轮廓包络进行等比例缩小或放大时,该轮廓包络上的各个追踪控点对应移动;最后,确定缩小或放大后的该条轮廓包络上各个追踪控点,并更新各个追踪控点的坐标信息、控制类型信息和所处心肌节段信息。
需要注意的是对于内膜,其拖拽后的位置不可超过心肌内膜边缘,对于外膜控点,对其拖拽后的位置不可超过心肌外膜的边缘。且拖拽后还需再次进行图4所示的追踪控点去重处理。
在图1至图5中任一幅图的基础上,图6还提供其他实施例超声心动图心肌包络量测方法流程图,包括:
步骤S61:根据图1至图5中任一幅图,依次确定所述超声心动图各帧的超声心脏切面图像中,相邻心肌节段之间的轮廓包络;各条所述轮廓包络对应所述超声心动图的各帧,分别为第M轮廓包络,所述M={2,3…m},m为大于1的整数。
步骤S61:根据相邻X帧所述超声心脏切面图像之间对应轮廓包络的变化,计算所述左心室的心肌应变,所述X∈M。
本实施例根据依次确定的各帧轮廓包络,能够确定在M帧过程中左心室心肌的应变状态,有利于辅助医护人员对心肌形变状态进行准确、客观、定量的评价,有利于对各种心脏疾病在临床阶段早期症状的判定。
图7示出了本申请一实施例提供的超声心动图心肌包络量测装置,参照图7,至少一个处理器71,例如CPU(Central Processing Unit,中央处理器),至少一个通信接口73,存储器74,至少一个通信总线72。其中,通信总线72用于实现这些组件之间的连接通信。其中,通信接口73可以包括显示屏、键盘,可选通信接口73还可以包括标准的有线接口、无线接口。存储器74可以是高速RAM存储器,也可以是非不稳定的存储器,例如至少一个磁盘存储器。存储器74可选的还可以是至少一个位于远离前述处理器71的存储装置。其中处理器71可以结合图5所描述的装置,存储器74中存储应用程序,且处理器71调用存储器74中存储的程序代码,以用于执行上述任一方法步骤。
其中,通信总线72可以是外设部件互连标准总线或扩展工业标准结构总线等。通信总线72可以分为地址总线、数据总线、控制总线等。为便于表示,图7中仅用一条粗线表示,但并不表示仅有一根总线或一种类型的总线。
其中,存储器74可以包括易失性存储器,例如随机存取存储器;存储器也可以包括非易失性存储器,例如快闪存储器,硬盘或固态硬盘;存储器74还可以包括上述种类的存储器的组合。
其中,处理器71可以是中央处理器,网络处理器,或者中央处理器和网络处理器的组合。
其中,处理器71还可以进一步包括硬件芯片。上述硬件芯片可以是专用集成电路,可编程逻辑器件或其组合。可选地,存储器74还用于存储程序指令。
本申请还提供一种计算机存储介质,所述计算机存储介质中存储有至少一条程序指令,所述至少一条程序指令被处理器加载并执行以实现图1至图6中任一幅图所示的方法。
显然,上述实施例仅仅是为清楚地说明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本申请创造的保护范围之中。
Claims (8)
1.一种超声心动图心肌包络量测方法,其特征在于,所述超声心动图心肌包络量测方法包括以下步骤:
获取超声心动图,所述超声心动图包括M帧超声心脏切面图像,所述M={1,2,3…m},m为大于1的整数;
确定所述超声心动图中第N帧超声心脏切面图像,所述N∈M;
将在所述第N帧超声心脏切面图像中,包围特定心腔的心肌,划分为Z个心肌节段,所述Z={1,2,3…z},z为大于1的整数;
依次确定各个所述心肌节段的追踪控点;所述追踪控点的控点类型包括:用于追踪心肌内膜的内膜控点、用于追踪心肌外膜的外膜控点,和,用于追踪心肌中膜的中膜控点;
根据所述追踪控点,确定在所述第N帧超声心脏切面图像中,相邻心肌节段之间的轮廓包络;所述轮廓包络包括:根据所述内膜控点确定的内膜包络、根据所述外膜控点确定的外膜包络,和,根据所述中膜控点确定的中膜包络;
确定所有所述心肌节段的追踪控点的步骤完成后,还进行:
根据所述追踪控点的坐标,确定控点集合;
判断所述控点集合中,各个所述追踪控点是否存在重复控点;
确定一追踪控点存在重复控点,获取所述重复控点对应的心肌节段和控点类型;
在与所述重复控点同一心肌节段,重新确定同一控点类型的追踪控点,并提取所述追踪控点的坐标;直至所述控点集合中不存在重复控点;
在所述根据所述追踪控点,确定在所述第N帧超声心脏切面图像中,相邻心肌节段之间的轮廓包络的步骤进行前,在所述依次确定各个所述心肌节段的追踪控点的步骤完成后,还进行:
确定控点集合,所述控点集合为所有追踪控点的集合,包括起始控点和结束控点;
使得所述控点集合中不存在重复控点;
判断所述起始控点和结束控点之间的距离是否大于预设阈值;
确定所述起始控点和结束控点之间的距离大于预设阈值,在所述控点集合起始控点的周围新插入一追踪控点;
将所述新插入的所述追踪控点更新为结束控点。
2.如权利要求1所述的超声心动图心肌包络量测方法,其特征在于,所述将在所述第N帧超声心脏切面图像中,包围特定心腔的心肌,划分为Z个心肌节段的步骤,包括:
在所述第N帧超声心脏切面图像中,沿着所述心肌的周向,依次将包围所述特定心腔的心肌,划分为Z个心肌节段。
3.如权利要求1所述的超声心动图心肌包络量测方法,其特征在于,所述依次确定各个所述心肌节段的追踪控点的步骤中,确定第P个心肌节段的追踪控点,所述P∈Z,包括以下步骤:
获取第P个心肌节段所在图像区域各像素点的灰度值;
根据所述图像区域各像素点的灰度值,确定所述第P个心肌节段的心肌切面边缘;所述心肌切面边缘包括相对的内边缘和外边缘;
在所述内边缘处,标注至少三个追踪像素点作为内膜控点,并提取所述内膜控点的坐标;
在所述外边缘处,标注至少三个追踪像素点作为外膜控点,并提取所述外膜控点的坐标;
在所述内边缘和所述外边缘之间,标注至少三个追踪像素点作为中膜控点,并提取所述中膜控点的坐标。
4.如权利要求1所述的超声心动图心肌包络量测方法,其特征在于,所述超声心动图中的M帧超声心脏切面图像,为按同一切面方向获取的心脏位置超声图像。
5.如权利要求4所述的超声心动图心肌包络量测方法,其特征在于,所述切面方向包括短轴切面。
6.一种超声心动图心肌包络量测方法,其特征在于,所述超声心动图心肌包络量测方法,还包括以下步骤:
根据权利要求1至5中任一项所述的超声心动图心肌包络量测方法,依次确定所述超声心动图各帧的超声心脏切面图像中,相邻心肌节段之间的轮廓包络;
根据相邻X帧所述超声心脏切面图像之间对应轮廓包络的变化,计算所述特定心腔的心肌应变,所述X∈M。
7.一种超声心动图心肌包络量测装置,其特征在于,所述超声心动图心肌包络量测装置包括存储器和处理器,所述存储器中存储有至少一条程序指令,所述处理器通过加载并执行所述至少一条程序指令以实现如权利要求1至5任一项所述的超声心动图心肌包络量测方法。
8.一种计算机存储介质,其特征在于,所述计算机存储介质中存储有至少一条程序指令,所述至少一条程序指令被处理器加载并执行以实现如权利要求1至5任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011631107.8A CN112842384B (zh) | 2020-12-30 | 2020-12-30 | 超声心动图心肌包络量测方法、装置和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011631107.8A CN112842384B (zh) | 2020-12-30 | 2020-12-30 | 超声心动图心肌包络量测方法、装置和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112842384A CN112842384A (zh) | 2021-05-28 |
CN112842384B true CN112842384B (zh) | 2023-05-30 |
Family
ID=75999788
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011631107.8A Active CN112842384B (zh) | 2020-12-30 | 2020-12-30 | 超声心动图心肌包络量测方法、装置和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112842384B (zh) |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6674879B1 (en) * | 1998-03-30 | 2004-01-06 | Echovision, Inc. | Echocardiography workstation |
US7248725B2 (en) * | 2004-01-07 | 2007-07-24 | Ramot At Tel Avia University Ltd. | Methods and apparatus for analyzing ultrasound images |
RU2448649C2 (ru) * | 2006-05-25 | 2012-04-27 | Конинклейке Филипс Электроникс Н.В. | Количественная оценка и отображение утолщения стенки камеры сердца |
EP2529666B1 (en) * | 2010-01-29 | 2017-01-18 | Hitachi, Ltd. | Ultrasonic diagnosis device and method used therefor to track measurement point |
WO2017206023A1 (zh) * | 2016-05-30 | 2017-12-07 | 深圳迈瑞生物医疗电子股份有限公司 | 一种心脏容积识别分析系统和方法 |
CN111012377B (zh) * | 2019-12-06 | 2020-11-03 | 北京安德医智科技有限公司 | 超声心动图心脏参数计算以及心肌应变测量方法、装置 |
-
2020
- 2020-12-30 CN CN202011631107.8A patent/CN112842384B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112842384A (zh) | 2021-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mor-Avi et al. | Real-time 3-dimensional echocardiography: an integral component of the routine echocardiographic examination in adult patients? | |
RU2653274C2 (ru) | Связанная сегментация в стандартных и контрастных ультразвуковых 3d-изображениях | |
US8265363B2 (en) | Method and apparatus for automatically identifying image views in a 3D dataset | |
EP2030042B1 (en) | Quantification and display of cardiac chamber wall thickening | |
WO2017206023A1 (zh) | 一种心脏容积识别分析系统和方法 | |
US20060173317A1 (en) | Ultrasound diagnostic system for automatically detecting a boundary of a target object and method of using the same | |
BR112015025074B1 (pt) | Sistema de imageamento por ultrassom e método para gerar e avaliar vistas bidimensionais padrão a partir de dados de volume ultrassônico tridimensional | |
BR112016009829B1 (pt) | Sistema de imageamento diagnóstico ultrassônico que captura uma pluralidade de planos de vista padrão do coração e método para capturar uma pluralidade de planos de vista padrão do coração | |
JP2005218796A (ja) | 医用画像処理装置および医用画像処理方法 | |
JP2009540911A (ja) | 三次元超音波イメージングのための方法、装置及びコンピュータプログラム | |
JP6393698B2 (ja) | 3dイメージングシステムのためのデータ表示及び処理アルゴリズム | |
US10918357B2 (en) | Methods and systems for automatically determining an anatomical measurement of ultrasound images | |
JP7149516B2 (ja) | 検査支援プログラム、検査支援方法および検査支援装置 | |
EP4061231B1 (en) | Intelligent measurement assistance for ultrasound imaging and associated devices, systems, and methods | |
KR20200115488A (ko) | 심장 벽의 운동을 기록하기 위한 이미지 분석 | |
JP6039876B2 (ja) | 左心室のねじれを測定するシステム及び方法 | |
Yue et al. | Speckle tracking in intracardiac echocardiography for the assessment of myocardial deformation | |
CN112842384B (zh) | 超声心动图心肌包络量测方法、装置和存储介质 | |
US20230346339A1 (en) | Systems and methods for imaging and measuring epicardial adipose tissue | |
WO2017069699A1 (en) | Medical image processing methods and systems for assessing right ventricular function | |
Aggeli et al. | Real-time three-dimensional echocardiography: never before clinical efficacy looked so picturesque | |
JP6498316B2 (ja) | 超音波撮像装置、および画像処理方法 | |
US20240099692A1 (en) | Guided acquisition of a 3d representation of an anatomical structure | |
US20200323454A1 (en) | Cardiac trajectory curve analysis for clinical decision making and analysis | |
JP7165541B2 (ja) | ボリュームデータ処理装置、方法及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |