CN108042154B - 二维超声心动图序列中心肌形状、运动和变形分析方法 - Google Patents
二维超声心动图序列中心肌形状、运动和变形分析方法 Download PDFInfo
- Publication number
- CN108042154B CN108042154B CN201711296199.7A CN201711296199A CN108042154B CN 108042154 B CN108042154 B CN 108042154B CN 201711296199 A CN201711296199 A CN 201711296199A CN 108042154 B CN108042154 B CN 108042154B
- Authority
- CN
- China
- Prior art keywords
- myocardial
- point
- cardiac cycle
- motion
- shape
- 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
- 230000002107 myocardial effect Effects 0.000 title claims abstract description 266
- 230000033001 locomotion Effects 0.000 title claims abstract description 102
- 238000004458 analytical method Methods 0.000 title claims abstract description 24
- 238000011265 2D-echocardiography Methods 0.000 title claims description 14
- 230000000747 cardiac effect Effects 0.000 claims abstract description 103
- 210000004165 myocardium Anatomy 0.000 claims abstract description 57
- 238000007781 pre-processing Methods 0.000 claims abstract description 51
- 238000012800 visualization Methods 0.000 claims abstract description 49
- 238000006073 displacement reaction Methods 0.000 claims abstract description 46
- 238000004364 calculation method Methods 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000012512 characterization method Methods 0.000 claims description 19
- 239000000126 substance Substances 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000008676 import Effects 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 4
- 230000036541 health Effects 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 8
- 238000012545 processing Methods 0.000 abstract description 6
- 230000006872 improvement Effects 0.000 description 8
- 230000002861 ventricular Effects 0.000 description 5
- 238000003759 clinical diagnosis Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 208000019622 heart disease Diseases 0.000 description 3
- 238000003709 image segmentation Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 208000029078 coronary artery disease Diseases 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 210000001174 endocardium Anatomy 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000010247 heart contraction Effects 0.000 description 1
- 230000004217 heart function Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 210000005240 left ventricle Anatomy 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 238000000718 qrs complex Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
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/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0883—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
-
- 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
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提出了一种二维超声心动图序列中心肌形状、运动和变形分析方法,包括以下步骤:首先通过预处理单元对待分析的二维超声心动图序列进行预处理,从而获取心动周期心肌壁各点坐标,并将心肌壁各点在心动周期各时刻位置的方向于可视化单元上显示;之后心肌运动形态计算单元根据心动周期心肌壁各点坐标分别计算心肌形状、位移、速度以及纵向应变,从而直观地理解心脏力学,特别是心肌形态、运动和变形。本发明对心脏周期内每个时刻的图像进行分析,从而显示心脏运动过程中心肌形态、运动和变形情况,满足临床实践对实时性处理的要求,还能够有效规避临床上某些静止图像存在偶然性的风险,便于临床医生做出正确的判断。
Description
技术领域
本发明涉及医学图像处理和智能健康管理领域,尤其是一种二维超声心动图序列中心肌形状、运动和变形分析方法。
背景技术
冠心病等心脏疾病是目前导致人类疾病死亡的主要原因之一,在心脏疾病的诊断和治疗过程中,医学影像技术发挥着重要作用。超声心动图具有无损伤性、安全、可靠等特点,广泛应用于临床上心脏病的诊断。通过二维超声心动图来研究心脏的运动,尤其是对左心室的研究,则是关注的中心。尽管可以直接从二维实时超声图像来观察心脏的运动,但对其诊断质量仍有可能进行较大的改善。由于实时观察基本属于定性的评价,依赖于操作者的经验和熟练程度,受观察者因素的制约甚多。近年来随着计算机相关技术的迅速发展以及图形图像技术的日渐成熟,医学图像处理技术对医学科研及临床实践的作用和影响日益增大,使临床医学对人体内部病变部位的观察更直接、更清晰,确认率也更高。
从研究情况来看,心脏超声图像分析技术主要分为基于像素点亮度和基于形状两种方式,前一种通过直接计算像素点的亮度变化来估计运动情况,如计算机视觉中的光流场计算等;后一种利用形变模型等方法来跟踪心脏图像中特征点的运动变化。光流理论要求图像序列中对应点的亮度应是恒定不变的,直接将传统方法运用到心脏超声图像分析上很难准确分析其运动情况;另一方面,光流较大的计算量也是一个需要解决的问题,以满足实时处理的要求。此外,物体的运动可以为图像分割提供比静止条件下更多的有用信息。心内(外)膜等运动轮廓的准确提取,对于心脏的三维重建、心腔体积的测量以及收缩过程中室壁的增厚判断都是非常重要的。国内外针对心脏图像的动态分割进行了深入研究,先后提出了活动轮廓模型、活动形体模型、活动外观模型和水平集四类方法。活动外观模型能够将图像的形状和纹理信息结合起来考虑,可以有效地提高跟踪的准确程度。统计模型的一个难点在于要求建立模型时各点之间有很好的匹配关系。
专利号为201310639848.4的发明专利《一种PET心脏图像中心肌轮廓分割方法和装置》,包括:对接收到的PET心脏图像进行粗分割,得到心肌区域的初始轮廓线;将水平集图像分割模型的函数描述中的曲线长度正则项中的导数阶数设置为分数阶;使用所述水平集图像分割模型和选取所述初始轮廓线上每一个像素点相邻的至少三个像素点对所述初始轮廓线进行曲线演化计算;将得到的曲线演化计算结果作为所述PET心脏图像中心肌区域的轮廓线。该方法存在的不足之处是:(1)PET心脏图像是一种三维图像,心脏周期内每一时刻需要处理的数据量都非常巨大,难以满足临床实践对实时性处理的要求。(2)利用PET心脏图像仅仅提取了心脏轮廓线,未解析心肌运动和变形信息,对完善临床诊断作用有限。
发明内容
本发明的目的在于提出一种二维超声心动图序列中心肌形状、运动和变形分析系统,以及通过该系统进行分析的方法。
为了解决上述技术问题,本发明提出一种二维超声心动图序列中心肌形状、运动和变形分析系统:
所述分析系统包括依次相连的预处理单元、心肌运动形态计算单元和可视化单元,以及与预处理单元相连的采集设备;
所述预处理单元与可视化单元相连;
所述采集设备用于采集二维超声心动图;
所述预处理单元从采集设备导入二维超声心动图的序列数据,并对该序列数据进行预处理;所述预处理单元输出心动周期内心肌壁各点(心肌点)坐标至心肌运动形态计算单元,输出心动周期内心肌壁各点方向至可视化单元;
所述心肌运动形态计算单元利用预处理单元输出的心肌壁各点坐标,计算完整心动周期内心肌的形状、运动(位移和速度)和变形(纵向应变)信息,并输出至可视化单元;
所述可视化单元对心肌壁各点方向以及心肌形状、运动和变形信息进行可视化显示。
作为本发明二维超声心动图序列中心肌形状、运动和变形分析系统的改进:
所述心肌运动形态计算单元包括心肌形状表征部件、心肌运动表征部件和心肌形变表征部件;心肌形状表征部件、心肌运动表征部件和心肌形变表征部件均与预处理单元和可视化单元相连;
所述心肌形状表征部件利用预处理单元输出的心动周期心肌壁各点坐标计算完整心动周期内心肌的形状信息,并通过可视化单元进行显示;
所述心肌运动表征部件利用预处理单元输出的心动周期心肌壁各点坐标计算心肌壁各点在心动周期内各时刻位置的位移和速度信息,并通过可视化单元进行显示;
所述位移包括径向位移和纵向位移,速度包括径向速度和纵向速度;
所述心肌形变表征部件利用预处理单元输出的心动周期心肌壁各点坐标计算心肌的纵向应变信息,并通过可视化单元进行显示。
为了解决上述技术问题,本发明还提出一种二维超声心动图序列中心肌形状、运动和变形分析方法,依次进行以下步骤:
S1、预处理单元导入并预处理待分析的二维超声心动图序列,并输出心动周期内心肌壁各点坐标至心肌运动形态计算单元,输出心动周期内心肌壁各点方向至可视化单元进行显示;
S2、心肌运动形态计算单元根据步骤S1中心肌壁各点坐标分别计算完整心动周期内心肌的形状、运动(位移和速度)和变形(纵向应变)信息,并通过可视化单元进行显示。
作为本发明二维超声心动图序列中心肌形状、运动和变形分析方法的改进:
所述心肌运动形态计算单元包括心肌形状表征部件、心肌运动表征部件和心肌形变表征部件;
所述步骤S2中完整心动周期内心肌的形状信息包括心动周期开始阶段、收缩末期和舒张末期的心肌形状,通过心肌形状表征部件在可视化单元上显示其的方法如下:
作为本发明二维超声心动图序列中心肌形状、运动和变形分析方法的进一步改进:
所述步骤S2中完整心动周期内心肌的形状信息还包括收缩末期和舒张末期整个心肌形状的总长度和相对差比值,通过心肌形状表征部件计算其的方法包括以下步骤:
2.1、收缩末期和舒张末期整个心肌形状的总长度的计算步骤如下:
其中,为j时刻第i-1个心肌点的x轴坐标,为j时刻第i+1个心肌点的x轴坐标;为j时刻第i-1个心肌点的y轴坐标,为j时刻第i+1个心肌点的y轴坐标;Δlongx(j,i)为j时刻第i个心肌点的x轴坐标差,Δlongy(j,i)为j时刻第i个心肌点的y轴坐标差;|Δlong(j,i)|对应两个位置差的范数;L定义为某时刻t(t=j)整个心肌形状的总长度,通过对所有点的范数进行求和获得;
2.2、收缩末期和舒张末期整个心肌形状总长度相对差比值的计算:
根据步骤2.1计算得到整个心肌形状的最大总长度LED(舒张末期)和最小总长度LES(收缩末期);
作为本发明二维超声心动图序列中心肌形状、运动和变形分析方法的进一步改进:
所述步骤S2中完整心动周期内心肌的运动信息包括心肌壁各点在心动周期内各时刻位置的径向位移、纵向位移、径向速度和纵向速度,通过心肌运动表征部件计算其的方法包括如下步骤:
3.1、计算心肌壁各点在心动周期内各时刻位置的径向位移和纵向位移;
由于
因此,心肌点径向位移ur和纵向位移ul可通过以下公式进行计算:
其中,eradial,x表示eradial在x轴方向的分量,eradial,y表示eradial在y轴方向的分量;elong,x表示elong在x轴方向的分量,elong,y表示elong在y周方向的分量;
3.2、计算步骤3.1输出的各点的径向速度和纵向速度;
所述径向速度vr和纵向速度vl从位移对时间的导数得到,计算步骤如下:
详细计算如下:
其中,t_interval为0.01秒。
作为本发明二维超声心动图序列中心肌形状、运动和变形分析方法的进一步改进:
所述步骤S1中预处理单元中所输出心动周期内心肌壁各点方向的计算方法如下:
1.1、心肌壁各点在心动周期各时刻位置的纵向方向elong计算步骤如下:
而
1.2、心肌壁各点在心动周期各时刻位置的径向方向eradial则可将对应纵向方向elong旋转90度获得,计算公式如下:
作为本发明二维超声心动图序列中心肌形状、运动和变形分析方法的进一步改进:
所述步骤S2中完整心动周期内心肌的变形信息包括心肌壁各点在心动周期内各时刻位置的纵向应变,通过心肌形变表征部件计算其所采用的计算公式如下:
作为本发明二维超声心动图序列中心肌形状、运动和变形分析方法的进一步改进:
所述步骤S1中预处理单元预处理待分析的二维超声心动图序列的方法如下:
所述预处理单元首先对导入的待分析二维超声心动图序列进行漂移校正;然后利用散斑跟踪技术对心肌壁进行自动跟踪,沿着二维超声心动图序列传播方向对心肌进行分割,并按照QRS波群出现的先后进行排序,使数据在时间上对齐。
本发明与现有技术相比,具有如下技术优势:
1、本发明建立了利用二维超声心动图序列评估心脏运动功能的数学表征模型(即,本发明二维超声心动图序列中心肌形状、运动和变形分析系统),该模型利用可视化的方式显示心肌特征信息,能十分直观地理解心脏力学,特别是心肌形态、运动和变形。
2、本发明结合了心脏力学,能够揭示心肌形态、运动(位移和速度)和变形(纵向应变)之间的关系,进一步完善临床诊断。
3、本发明采用了动态技术显示心脏运动过程中心肌形态、运动和变形情况,能够有效规避临床上某些静止图像存在偶然性的风险,便于临床医生做出正确的判断。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细说明。
图1是本发明二维超声心动图序列中心肌形状、运动和变形分析方法的原理流程图;本图还体现了本发明二维超声心动图序列中心肌形状、运动和变形分析系统的模块示意图。
图2是本发明实施例中在心动周期开始阶段心肌壁各点的纵向方向和径向方向;
图3是本发明实施例中在心动周期开始阶段、收缩末期和舒张末期的心肌形状;
图4是本发明实施例中心肌壁的三个不同点(i=6,19,31)沿整个心动周期的径向位移;
图5是本发明实施例中心肌壁的三个不同点(i=6,19,31)沿整个心动周期的纵向位移;
图6是本发明实施例中位于中位室间隔位置的心肌点沿整个心动周期的径向和纵向速度表征;
图7是本发明实施例中位于中位室间隔位置的心肌点沿整个心动周期的纵向应变。
具体实施方式
下面结合具体实施例对本发明进行进一步描述,但本发明的保护范围并不仅限于此。
实施例1、二维超声心动图序列中心肌形状、运动和变形分析系统,如图1所示,包括依次相连的采集设备1、预处理单元2、心肌运动形态计算单元3和可视化单元4,其中预处理单元2与可视化单元4相连。
采集设备1用于采集二维超声心动图。
预处理单元2用于导入采集设备1中的二维超声心动图的序列数据并对该序列数据进行预处理;本实施例中,预处理单元2输出心动周期心肌壁各点(即,心肌点)坐标至心肌运动形态计算单元3,还计算心肌壁各点在心动周期各时刻位置的方向(纵向和径向),并将计算所得的各点方向通过可视化单元4进行显示。
心肌运动形态计算单元3利用预处理单元2输出的心肌壁各点坐标,计算完整心动周期内心肌的形状、运动(位移和速度)和变形(纵向应变)信息,并输出至可视化单元4;
心肌运动形态计算单元3包括心肌形状表征部件31、心肌运动表征部件32和心肌形变表征部件33。心肌形状表征部件31、心肌运动表征部件32和心肌形变表征部件33均与预处理单元2和可视化单元4相连。
心肌形状表征部件31利用预处理单元2输出的心动周期心肌壁各点坐标计算完整心动周期内心肌的形状信息,并通过可视化单元4进行显示;
心肌运动表征部件32利用预处理单元2输出的心动周期心肌壁各点坐标计算心肌壁各点在心动周期内各时刻位置的位移和速度信息,并通过可视化单元4进行显示;位移包括径向位移和纵向位移,速度包括径向速度和纵向速度。
心肌形变表征部件33利用预处理单元2输出的心动周期心肌壁各点坐标计算心肌的纵向应变信息,并通过可视化单元4进行显示。
可视化单元4对心肌壁各点方向以及心肌的形状、运动(位移和速度)和变形(纵向应变)信息进行可视化显示,便于直观判断心脏运动功能的健康程度。
本发明还提出一种利用上述系统进行的二维超声心动图序列中心肌形状、运动和变形分析方法,如图1-7所示,首先通过预处理单元2对待分析的二维超声心动图序列进行预处理,从而获取心动周期心肌壁各点坐标,并将心肌壁各点在心动周期各时刻位置的方向于可视化单元4上显示;之后心肌运动形态计算单元3根据心动周期心肌壁各点坐标分别计算心肌形状、位移、速度以及纵向应变。
本实施例采用的测试数据来自浙江省中医院其中一位自愿者的二维超声心动图序列,含有四个腔室,总共83帧二维超声心动图像,心动周期开始阶段的帧号为1、收缩末期的帧号为27和舒张末期的帧号为83。
具体工作内容包括如下步骤:
步骤1、预处理单元2从采集设备1导入二维超声心动图,并对该二维超声心动图的序列数据(即二维超声心动图序列)并进行预处理后,数据输出至心肌运动形态计算单元3。具体工作流程如下:
(1.1)、预处理:为了使心动周期开始和结束的位置一致,对二维超声心动图序列进行漂移校正;然后利用散斑跟踪技术对心肌壁进行自动跟踪,沿着二维超声心动图序列传播方向对心肌进行分割,输出心动周期心肌壁各点坐标(即各心肌点对应的坐标矩阵)至心肌运动形态计算单元3(即,心肌形状表征部件31、心肌运动表征部件32和心肌形变表征部件33),并按照QRS波群(即,正常心电图中幅度最大的波群)出现的先后进行排序,使数据在时间上对齐。
二维超声心动图序列漂移校正方法和散斑跟踪技术均为现有技术,故在本说明书中不对其进行详细阐述。
(1.2)、计算步骤(1.1)中心肌壁各点在心动周期各时刻位置的方向(纵向和径向),并输出至可视化单元4进行显示,具体计算方法如下:
(1.2.1)、计算每个心肌点在心动周期各时刻位置的纵向方向elong;
而
(1.2.2)、计算每个心肌点在心动周期各时刻位置的径向方向eradial;
心肌壁各点在心动周期各时刻位置的径向方向eradial可通过将对应心肌点的纵向方向elong旋转90度获得,即:
如图2所示,本实施例中志愿者1的二维超声心动图序列通过预处理单元2输出至可视化单元4进行显示的心动周期开始阶段心肌壁各点的方向;如图2中箭头所示,指向中心区域的箭头分别表示对应心肌点的径向方向,沿着心肌轮廓进行指向的箭头别表示对应心肌点的纵向方向。
步骤2、根据步骤(1.1)预处理单元2输出的心动周期心肌壁各点的坐标通过心肌形状表征部件31计算心动周期开始阶段、收缩末期和舒张末期的心肌形状,以及收缩末期和舒张末期整个心肌形状的总长度和相对差比值。
(2.1)、计算心动周期开始阶段、收缩末期和舒张末期的心肌形状;
心肌形状表征部件31接收步骤(1.1)输出的心肌壁的各点坐标,心肌形状表征部件31按照心脏周期内对应时刻(即,动周期开始阶段、收缩末期和舒张末期)分别将每个心肌点在同一个坐标系中进行表示,并通过可视化单元4进行显示,从而获取心脏周期内所对应时刻的心肌形状。
其中,为每个心肌点分别在心动周期开始阶段(出现第1帧心动图时刻)、收缩末期(出现第27帧心动图时刻)和舒张末期(出现第83帧心动图时刻)相对应的x轴坐标;为每个心肌点分别在心动周期开始阶段(出现第1帧心动图时刻)、收缩末期(出现第27帧心动图时刻)和舒张末期(出现第83帧心动图时刻)相对应时刻的y轴坐标。
通过可视化单元4可以利用坐标矩阵显示在心动周期开始阶段(出现第1帧心动图时刻)、收缩末期(出现第27帧心动图时刻)和舒张末期(出现第83帧心动图时刻)的心肌形状。本实施例中志愿者1的二维超声心动图序列通过步骤(2.1)于可视化单元4显示的心肌形状如图3所示,开始阶段与舒张末期的心肌形状相重合,收缩末期与舒张末期心肌形状可以直观地进行区分,从而便于工作人员对心脏运动过程中心肌形态进行分析。
(2.2)、计算收缩末期(帧号为27)和舒张末期(帧号为83)整个心肌形状的总长度和相对差比值;
(2.2.1)、收缩末期(帧号为27)和舒张末期(帧号为83)整个心肌形状总长度的计算公式如下:
其中,为j时刻第i-1个心肌点的x轴坐标,为j时刻第i+1个心肌点的x轴坐标;为j时刻第i-1个心肌点的y轴坐标,为j时刻第i+1个心肌点的y轴坐标;Δlongx(j,i)为j时刻第i个心肌点的x轴坐标差,Δlongy(j,i)为j时刻第i个心肌点的y轴坐标差;|Δlong(j,i)|对应两个位置差的范数;L定义为某时刻t(t=j)整个心肌形状的总长度,通过对所有点的范数进行求和获得;
注:收缩末期(出现第27帧心动图时刻)心肌形状的总长度为最小总长度LES;舒张末期(出现第83帧心动图时刻)心肌形状的总长度为最大总长度LED。
(2.2.2)、收缩末期(帧数为27)和舒张末期(帧数为83)整个心肌形状总长度相对差比值的计算:
通过计算得到本实施例中收缩末期(帧数为27)的整个心肌形状的总长度为221.441mm,收缩末期(帧数为27)的整个心肌形状的总长度为198.368mm,相对差比值为0.104。
步骤3、根据步骤(1.1)输出的心动周期内心肌壁各点位置通过心肌运动表征部件32计算径向和纵向位移以及径向和纵向速度。具体工作流程如下:
(3.1)、计算心肌壁各点在心动周期内各时刻位置的径向和纵向位移;
由于:
因此,心肌点径向位移ur和纵向位移ul可通过以下公式进行计算:
其中,eradial,x表示eradial在x轴方向的分量,eradial,y表示eradial在y轴方向的分量;elong,x表示elong在x轴方向的分量,elong,y表示elong在y周方向的分量。
如图4所示,本实施例中由心肌运动表征部件32输出至可视化单元4显示的心肌壁三个不同点(即,第6、19和31个心肌点)沿整个心动周期的径向位移;如图5所示,本实施例中由心肌运动表征部件32输出至可视化单元4显示的心肌壁三个不同点(即,第6、19和31个心肌点)沿整个心动周期的纵向位移。
(3.2)、计算步骤(3.1)输出的各点的径向速度和纵向速度;
径向速度vr和纵向速度vl由位移对时间的导数得到,计算步骤如下:
详细计算如下:
其中,t_interval为0.01秒。
如图6所示,本实施例中由心肌运动表征部件32输出至可视化单元4显示的中位室间隔位置的心肌点沿整个心动周期的径向和纵向速度表征;
步骤4、根据步骤(1.1)输出的心动周期内心肌壁各点的坐标通过心肌形变表征部件33计算其纵向应变,步骤如下:
由于心肌纵向应变与心脏收缩功能关系密切,在医学领域有重要的临床意义,计算公式如下:
根据心肌壁各点的坐标通过上述两个公式计算获得对应心肌点在时刻t0至ti内的纵向应变。如图7所示,本实施例中t0表示心动周期开始的时间,ti表示心动周期结束的时间,从而获得该心肌点xi沿整个心动周期的纵向应变,实现了对心脏运动过程中变形情况的分析。
注:t0和ti可以为心脏周期内任意前后两个时刻。
如图7所示,中位室间隔位置的心肌点由心肌形变表征部件33输出至可视化单元4显示的其沿整个心动周期的纵向应变(即,在中位室间隔位置的纵向应变)。
本发明建立了利用二维超声心动图序列评估心脏功能的表征模型(本发明提出的分析系统),并利用可视化的方式显示特征信息,能够十分直观地理解心脏力学,特别是心肌形态、运动和变形。本发明通过结合心脏力学分别建立了心肌形态、运动和变形的数学表征模型(即,心肌运动形态计算单元),从而利用心脏周期内每个时刻的各心肌点坐标即可计算出相应的心肌形状、位移和速度以及纵向应变特征,有助于了解这些特征之间的相互关系,进一步完善临床诊断。此外,本发明不只是对一幅二维超声心动图进行分析,而是对心脏周期内每个时刻的图像(实施例中总共83幅图像)进行分析,每一时刻对应一副图像,采用动态技术显示心脏运动过程中心肌形态、运动和变形情况,例如图2-7所示,工作人员可根据实际工作需要对不同时刻心肌形状进行对比分析,并针对具体心肌点的运动和变形情况进行分析,从而满足临床实践对实时性处理的要求,还能够有效规避临床上某些静止图像存在偶然性的风险,便于临床医生做出正确的判断。
最后,还需要注意的是,以上列举的仅是本发明的若干个具体实施例。显然,本发明不限于以上实施例,还可以有许多变形。本领域的普通技术人员能从本发明公开的内容直接导出或联想到的所有变形,均应认为是本发明的保护范围。
Claims (2)
1.利用二维超声心动图序列中心肌形状、运动和变形分析系统进行分析的方法,其特征在于:
分析系统包括依次相连的预处理单元(2)、心肌运动形态计算单元(3)和可视化单元(4),以及与预处理单元(2)相连的采集设备(1);
所述预处理单元(2)与可视化单元(4)相连;
所述采集设备(1)用于采集二维超声心动图;
所述预处理单元(2)从采集设备(1)导入二维超声心动图的序列数据,并对该序列数据进行预处理;所述预处理单元(2)输出心动周期内心肌壁各点坐标至心肌运动形态计算单元(3),输出心动周期内心肌壁各点方向至可视化单元(4);
所述心肌运动形态计算单元(3)利用预处理单元(2)输出的心肌壁各点坐标,计算完整心动周期内心肌的形状、运动和变形信息,并输出至可视化单元(4);
所述可视化单元(4)对心肌壁各点方向以及心肌形状、运动和变形信息进行可视化显示;
所述心肌运动形态计算单元(3)包括心肌形状表征部件(31)、心肌运动表征部件(32)和心肌形变表征部件(33);心肌形状表征部件(31)、心肌运动表征部件(32)和心肌形变表征部件(33)均与预处理单元(2)和可视化单元(4)相连;
所述心肌形状表征部件(31)利用预处理单元(2)输出的心动周期心肌壁各点坐标计算完整心动周期内心肌的形状信息,并通过可视化单元(4)进行显示;
所述心肌运动表征部件(32)利用预处理单元(2)输出的心动周期心肌壁各点坐标计算心肌壁各点在心动周期内各时刻位置的位移和速度信息,并通过可视化单元(4)进行显示;
所述位移包括径向位移和纵向位移,速度包括径向速度和纵向速度;
所述心肌形变表征部件(33)利用预处理单元(2)输出的心动周期心肌壁各点坐标计算心肌的纵向应变信息,并通过可视化单元(4)进行显示;
可视化单元(4)对心肌壁各点方向以及心肌的形状、运动和变形信息进行可视化显示,便于直观判断心脏运动功能的健康程度;
首先通过预处理单元(2)对待分析的二维超声心动图序列进行预处理,从而获取心动周期心肌壁各点坐标,并将心肌壁各点在心动周期各时刻位置的方向于可视化单元(4)上显示;之后心肌运动形态计算单元(3)根据心动周期心肌壁各点坐标分别计算心肌形状、位移、速度以及纵向应变;
方法为依次进行以下步骤:
S1、预处理单元(2)导入并预处理待分析的二维超声心动图序列,并输出心动周期内心肌壁各点坐标至心肌运动形态计算单元(3),输出心动周期内心肌壁各点方向至可视化单元(4)进行显示;
预处理单元(2)中所输出心动周期内心肌壁各点方向的计算方法如下:
1.1、心肌壁各点在心动周期各时刻位置的纵向方向elong计算步骤如下:
而
1.2、心肌壁各点在心动周期各时刻位置的径向方向eradial则可将对应纵向方向elong旋转90度获得,计算公式如下:
S2、心肌运动形态计算单元(3)根据步骤S1中心肌壁各点坐标分别计算完整心动周期内心肌的形状、运动和变形信息,并通过可视化单元(4)进行显示;
步骤S2中完整心动周期内心肌的形状信息包括心动周期开始阶段、收缩末期和舒张末期的心肌形状,通过心肌形状表征部件(31)在可视化单元(4)上显示其的方法如下:
步骤S2中完整心动周期内心肌的形状信息还包括收缩末期和舒张末期整个心肌形状的总长度和相对差比值,通过心肌形状表征部件(31)计算其的方法包括以下步骤:
2.1、收缩末期和舒张末期整个心肌形状的总长度的计算步骤如下:
其中,为j时刻第i-1个心肌点的x轴坐标,为j时刻第i+1个心肌点的x轴坐标;为j时刻第i-1个心肌点的y轴坐标,为j时刻第i+1个心肌点的y轴坐标;Δlongx(j,i)为j时刻第i个心肌点的x轴坐标差,Δlongy(j,i)为j时刻第i个心肌点的y轴坐标差;|Δlong(j,i)|对应两个位置差的范数;L定义为某时刻t整个心肌形状的总长度,t=j,通过对所有点的范数进行求和获得;
2.2、收缩末期和舒张末期整个心肌形状总长度相对差比值的计算:
根据步骤2.1计算得到整个心肌形状的最大总长度LED和最小总长度LES;
S3、完整心动周期内心肌的运动信息包括心肌壁各点在心动周期内各时刻位置的径向位移、纵向位移、径向速度和纵向速度,通过心肌运动表征部件(32)计算其的方法包括如下步骤:
3.1、计算心肌壁各点在心动周期内各时刻位置的径向位移和纵向位移;
由于
因此,心肌点径向位移ur和纵向位移ul可通过以下公式进行计算:
其中,eradial,x表示eradial在x轴方向的分量,eradial,y表示eradial在y轴方向的分量;elong,x表示elong在x轴方向的分量,elong,y表示elong在y轴方向的分量;
3.2、计算步骤3.1输出的各点的径向速度和纵向速度;
所述径向速度vr和纵向速度vl从位移对时间的导数得到,计算步骤如下:
详细计算如下:
其中,t_interval为0.01秒;
S4、完整心动周期内心肌的变形信息包括心肌壁各点在心动周期内各时刻位置的纵向应变,通过心肌形变表征部件(33)计算其所采用的计算公式如下:
2.根据权利要求1所述的利用二维超声心动图序列中心肌形状、运动和变形分析系统进行分析的方法,其特征在于:
步骤S1中预处理单元(2)预处理待分析的二维超声心动图序列的方法如下:
所述预处理单元(2)首先对导入的待分析二维超声心动图序列进行漂移校正;然后利用散斑跟踪技术对心肌壁进行自动跟踪,沿着二维超声心动图序列传播方向对心肌进行分割,并按照QRS波群出现的先后进行排序,使数据在时间上对齐。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711296199.7A CN108042154B (zh) | 2017-12-08 | 2017-12-08 | 二维超声心动图序列中心肌形状、运动和变形分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711296199.7A CN108042154B (zh) | 2017-12-08 | 2017-12-08 | 二维超声心动图序列中心肌形状、运动和变形分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108042154A CN108042154A (zh) | 2018-05-18 |
CN108042154B true CN108042154B (zh) | 2021-06-15 |
Family
ID=62123607
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711296199.7A Active CN108042154B (zh) | 2017-12-08 | 2017-12-08 | 二维超声心动图序列中心肌形状、运动和变形分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108042154B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108898622A (zh) * | 2018-07-05 | 2018-11-27 | 深圳大学 | 一种心脏的运动表征方法、装置及计算机可读存储介质 |
CN115618033B (zh) * | 2022-12-15 | 2023-03-14 | 山东科技大学 | 数据采集与处理方法、系统、计算机设备以及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6638221B2 (en) * | 2001-09-21 | 2003-10-28 | Kabushiki Kaisha Toshiba | Ultrasound diagnostic apparatus, and image processing method |
CN1774646A (zh) * | 2003-04-15 | 2006-05-17 | 皇家飞利浦电子股份有限公司 | 心壁应变成像 |
CN101317773A (zh) * | 2007-06-06 | 2008-12-10 | 株式会社东芝 | 超声波图像处理装置 |
CN101380238A (zh) * | 2007-09-07 | 2009-03-11 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置及方法 |
CN101791230A (zh) * | 2009-01-30 | 2010-08-04 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置、医用图像诊断装置 |
-
2017
- 2017-12-08 CN CN201711296199.7A patent/CN108042154B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6638221B2 (en) * | 2001-09-21 | 2003-10-28 | Kabushiki Kaisha Toshiba | Ultrasound diagnostic apparatus, and image processing method |
CN1774646A (zh) * | 2003-04-15 | 2006-05-17 | 皇家飞利浦电子股份有限公司 | 心壁应变成像 |
CN101317773A (zh) * | 2007-06-06 | 2008-12-10 | 株式会社东芝 | 超声波图像处理装置 |
CN101380238A (zh) * | 2007-09-07 | 2009-03-11 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置及方法 |
CN101791230A (zh) * | 2009-01-30 | 2010-08-04 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置、医用图像诊断装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108042154A (zh) | 2018-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3725442B2 (ja) | 医用画像診断装置及びその方法 | |
US8107703B2 (en) | Quantitative real-time 4D stress test analysis | |
EP2433567A1 (en) | Medical image diagnosis device and region-of-interest setting method therefor | |
US8144956B2 (en) | Ultrasonic diagnosis by quantification of myocardial performance | |
US8538098B2 (en) | Image processing method for displaying information relating to parietal motions of a deformable 3-D object | |
US8077944B2 (en) | Method, device, and computer programme for evaluating images of a cavity | |
EP1430837B1 (en) | Ultrasonic diagnostic device | |
JP6382036B2 (ja) | 超音波診断装置及び画像処理装置 | |
US20070135705A1 (en) | System and method for image based physiological monitoring of cardiovascular function | |
US20130253319A1 (en) | Method and system for acquiring and analyzing multiple image data loops | |
CN101283911A (zh) | 一种冠状动脉血管轴线的四维重建方法 | |
CN108042154B (zh) | 二维超声心动图序列中心肌形状、运动和变形分析方法 | |
CN110811653A (zh) | 一种测量正畸过程中牙齿三维移动的方法 | |
CN115661225B (zh) | 三维血管图像重建方法、装置、计算机设备、存储介质 | |
CN116883322A (zh) | 一种三维超声模型应用于心脏参数测量与管理方法及终端 | |
US20120008833A1 (en) | System and method for center curve displacement mapping | |
Mæhle et al. | Three‐dimensional echocardiography for quantitative left ventricular wall motion analysis: A method for reconstruction of endocardial surface and evaluation of regional dysfunction | |
Bruining et al. | Three-dimensional echocardiography paves the way toward virtual reality | |
JP3378446B2 (ja) | 超音波画像の処理方法及びその装置 | |
CN100438831C (zh) | 多视点冠脉造影序列影像的时间配准方法 | |
CN111904450B (zh) | 左心室中心和感兴趣区域的提取方法、装置及系统 | |
Walimbe et al. | Fully automatic segmentation of left ventricular myocardium in real-time three-dimensional echocardiography | |
Zheng et al. | Reconstruction of coronary vessels from intravascular ultrasound image sequences based on compensation of the in-plane motion | |
Lee et al. | Comparative study of left ventricular low wall motion with scar tissue using 4d left ventricular cardiac images | |
Blancher et al. | Time-varying, 3-D echocardiography using a fast-rotating probe |
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 |