CN117202842A - 确定心脏壁运动的方法 - Google Patents
确定心脏壁运动的方法 Download PDFInfo
- Publication number
- CN117202842A CN117202842A CN202280026693.4A CN202280026693A CN117202842A CN 117202842 A CN117202842 A CN 117202842A CN 202280026693 A CN202280026693 A CN 202280026693A CN 117202842 A CN117202842 A CN 117202842A
- Authority
- CN
- China
- Prior art keywords
- motion
- vector
- value
- heart
- point
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 239000013598 vector Substances 0.000 claims abstract description 191
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 230000000747 cardiac effect Effects 0.000 claims description 7
- 210000001174 endocardium Anatomy 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 4
- 210000004165 myocardium Anatomy 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims 3
- 238000003384 imaging method Methods 0.000 description 18
- 230000006870 function Effects 0.000 description 10
- 238000002059 diagnostic imaging Methods 0.000 description 8
- 238000005086 pumping Methods 0.000 description 7
- 241001522301 Apogonichthyoides nigripinnis Species 0.000 description 6
- 210000005240 left ventricle Anatomy 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 230000002107 myocardial effect Effects 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 241000124008 Mammalia Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 210000005242 cardiac chamber Anatomy 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 210000005003 heart tissue Anatomy 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1126—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb using a particular sensing technique
- A61B5/1128—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb using a particular sensing technique using image analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1107—Measuring contraction of parts of the body, e.g. organ, muscle
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/023—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/02—Measuring pulse or heart rate
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Optics & Photonics (AREA)
- Quality & Reliability (AREA)
- Multimedia (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
接收以固定时间间隔在时间上连续获取的图像序列,这些图像包括心脏的至少一部分的、具有心脏壁的图像数据,在每个图像中识别心脏壁的至少一部分,并在每个图像中定义所识别的心脏壁上的至少一个矢量点。本方法还包括基于序列确定矢量点随时间的运动,接收与心脏相关的第一参考。基于对矢量点和第一参考的跟踪,确定表示矢量点相对于参考的运动的运动矢量。
Description
技术领域
各个方面及其变型及一种用于评估用于提供心脏泵送功能的心肌组织的功能的方法和装置。
背景技术
可用各种技术来评估心脏的功能以及随后评估该功能是否与参考数据相符,以使心脏病专家能够诊断任何病理问题。这些技术是对心电图的分析和对多种成像技术记录的信息的评估。
发明内容
优选提供一种用于评估心肌组织的功能和心脏泵送功能的改进技术,并为这种评估提供数据,并且可选地,优选通过电子显示屏,将数据提供给心脏病专家或其他医疗从业者。
第一方面在电子计算设备中提供了一种确定心脏壁运动的方法。该方法包括接收以固定时间间隔在时间上连续获取的图像序列,图像包括心脏的至少一部分的、具有心脏壁的图像数据,在每个图像中识别心脏壁的至少一部分,并在每个图像中定义所识别的心脏壁上的至少一个矢量点。该方法还包括基于序列确定矢量点随时间的运动,并接收与心脏相关的第一参考。基于对矢量点和第一参考的跟踪,来确定表示矢量点相对于参考的运动的运动矢量。
通过确定一个或多个运动矢量,可以确定心脏壁的一个部分或多个部分的运动,特别是可以确定特定腔(如左心室)的壁的部分相对于参考的运动。通过在腔内提供参考,可以确定壁相对于例如腔中的中心点的运动。运动的幅度,特别是在指向或来自参考的方向上,可以提供心脏(特别是特定腔)的泵送功能的质量的指示。壁的一部分可以平行于所识别的壁的切平面移动。在这种情况下,运动对泵送功能贡献很小或没有贡献。通过在模型中的腔中提供的参考,并确定指示矢量点朝向或远离参考(参考点、参考曲线或其组合)的运动的运动矢量——可以确定有助于泵送的运动部分。所识别的心脏壁可以是心脏的外壁、特定腔的内壁或另一参考(例如在外壁和内壁之间,并且优选地提供为曲线)中的一个。
第一方面可以包括识别心脏壁上的矢量点集合,将该矢量点集合的第一子集分配给心脏壁的第一区域,基于第一子集中矢量点的确定的移动,确定第一子集中矢量点随时间的移动,确定第一区域的第一区域运动矢量,该第一区域运动矢量指示第一区域相对于第一参考的运动。
根据心脏病医学标准,特定腔的心内膜可以分成特定段。按标准化区域提供运动数据为医疗从业者设置诊断提供了改进的输入。
第一方面可以进一步包括将矢量点集合的第二子集分配给心脏壁的第二区域,确定第二子集中的矢量点随时间的运动,以及基于对第二子集中的矢量点的跟踪,确定第二区域的第二区域运动矢量,该第二区域运动矢量指示第二区域相对于第一参考的运动。
通过确定多个区域的运动矢量,可以向医疗从业者提供更多信息,允许他们比较各段的运动。
第一参考可以不同于第二参考。这可以提高准确性。
该方法可以包括:基于所确定的矢量点的至少一部分的运动来确定心脏随时间的运动,以及基于所确定的心脏运动来调整第一区域运动矢量以补偿心脏的运动。
在心跳周期中,心脏在受试者体内运动。在图像获取过程中,受试者处于稳定状态,这意味着图像捕捉到了心脏的实际运动。因此,在每个采集的图像或采集的帧中,心脏的位置可能不相同。通过确定矢量点中的共同运动或参考的运动(相对于所识别的心脏壁轮廓的至少一个位置确定),可以确定心脏的运动(特别是心跳)和心脏从一帧相对于另一帧的运动。该确定的心脏运动可用于调整与所识别的心脏壁上(因此例如在所识别的心脏壁轮廓上)的矢量点相关的所确定的运动矢量。
该方法可以包括确定运动矢量的向内矢量分量,向内矢量分量是平行于第一参考点和矢量点之间的直线的运动矢量分量。如上所述,特别是由识别的心脏壁表示的心内膜的向内运动与评估心脏的泵送功能相关。
运动矢量可以基于矢量点相对于第一参考点的最大位移。心跳中血液的位移和置换的血液量可以由心脏壁上的段或点在收缩期结束和舒张期结束之间的最大位移来确定。
参考可以设置在基部和顶点之间,特别是在该轨迹的三分之一至三分之二之间,更具体地,在该轨迹的一半和通常在该轨迹上距离基部三分之一处的点之间。参考可以包括单个参考点、多个参考点、开放曲线、闭合曲线或其任意组合。闭合曲线通常定义一个区域。可以提供从靠近基部的第一位置到远离基部的第二点的开放曲线。如果参考包括多个点,则可以在闭合曲线或开放曲线上提供这些点。如果参考包括多个点,则每个参考点可以对应于所识别的心脏壁上的矢量点。
该方法还可以包括接收指示矢量点相对于参考点的参考运动值的参考运动值,基于运动矢量导出运动矢量值,计算提供运动矢量值和参考运动值之间的差的指示的运动偏差值;以及提供用于在电子显示器上显示的运动偏差值。与参考组相比,这为医疗从业者提供了与特定受试者相关的数据。
第二方面提供了一种计算机程序产品,包括计算机可执行代码,当该代码被加载到处理单元或连接到处理单元的存储器中时,使计算机的处理单元执行第一方面的方法。
第三方面提供包括第二方面的计算机程序产品的非暂时性计算机可读介质。
第四方面提供了一种电子计算设备,包括:数据采集模块,设置为接收以固定时间间隔在时间上连续采集的图像序列,这些图像包括心脏的至少一部分的、具有心脏壁的图像数据;图像处理模块,设置为在每个图像中识别心脏壁的至少一部分,在每个图像中定义心脏壁上的至少一个矢量点,基于序列随时间跟踪矢量点;运动计算模块,设置为接收与心脏相关的第一参考,基于对矢量点和参考的跟踪来确定表示矢量点相对于参考的运动的运动矢量。
附图说明
下面将结合附图进一步详细讨论其各个方面及其变化。在附图中:
图1示出了医学成像系统;
图2示出了流程图;
图3示出了左心室的壁轮廓;
图4A示出了健康心脏的壁轮廓;
图4B示出了泵送功能降低的心脏的壁轮廓;以及
图5示出了一个靶心图。
具体实施方式
图1示出了医学成像系统100。医学成像系统100包括用于控制系统100的控制装置110。控制装置110包括中央处理单元112、存储模块114、用户界面控制器118和成像控制器116。控制设备110经由用户界面控制器118连接到作为用户输入设备的键盘142和作为用户输出设备的电子显示屏144。
成像控制器116连接到成像发射器132和成像接收器134。成像发射器132向成像接收器134发出像X射线辐射的能量—电磁波,声波、特别是超声波,以及其他能量或其组合。成像发射器132和成像接收器134是医学成像系统130的一部分。医学成像系统130可以是MRI系统、CT系统、标准X射线系统、超声成像系统、另一成像系统或其组合。医学成像系统被设置为获得哺乳动物、特别是人的心脏190的成像数据。成像数据优选地包括由近似恒定的时间间隔划分的图像。
成像控制器116可以被设置为处理成像数据接收器134,以提供心脏190或身体的任何其他部分的二维表示、三维重建。成像控制器116还可以被设置为提供三维重建的二维横截面。附加地或替代地,获得的二维图像被传递到中央处理单元112。也可以由中央处理单元112代替成像控制器116或结合成像控制器116来执行一些或所有步骤。成像控制器和中央处理单元112可以包括在同一个物理实体中。
中央处理单元112包括作为壁检测单元的边缘检测子单元122、用于分割数据的分割子单元124、用于在按时间分割的帧间跟踪数据点的运动检测子单元126和用于执行矢量操作的矢量计算子单元128。各种子单元可以以任何方式(包括临时的、永久的、易失性的、非易失性的、其他的或其组合)硬连线、编程。中央处理单元112可以进一步被设置为执行根据第一方面的方法和如下所述的方法,具有或不具有指示的任何选项及其任何变型。在任何选择的实现中,该单元最终在图像上提供多个位置,其中,这些位置随时间变化。这些位置由它们的坐标值标识,并且表示其中感兴趣的地方。
存储模块114被设置为具有存储在其上的图像数据,该图像数据可以通过医学成像系统130获得,如直接获得以及经处理后获得。此外,存储模块114可以具有存储在其上的数据,该数据包括可由中央处理单元112执行以使得中央处理单元能够如上所述进行操作的代码。
图2示出了可以由系统100、特别是控制装置110执行的过程的流程图200。流程图200的各部分可以总结如下:
202启动过程
204接收图像
206识别壁
208分割壁
210在壁上定义点
212定义参考
214所有图像完成?
216识别所有帧上的点,跟踪帧到帧的点或其组合
218定义心动周期边界
220定义运动部分
222如果存在多个循环,则对所有循环重复
224确定运动矢量
226确定平行于参考和点之间的线的矢量分量
228确定每个点或一段或一组点的运动矢量
230确定朝向参考的运动矢量分量
232所有点完成?
234确定壁和参考点的整体移动
236计算调整数据
238根据调整数据来调整数据值
240获取参考值
242计算运动偏差
244显示偏差数据
246显示指向参考的矢量分量数据
248结束过程
该过程从终止器202开始,并进行到步骤204,其中从心脏190获得图像序列。优选地,以在图像获取时间之间具有近似恒定的间隔的方式来获得图像。图像可以作为具有帧的流来提供。优选地提供图像,使得仔细观察下的心腔的二维区域尽可能大。例如,如果要评估左心室的工作,则以从基部到顶点的线尽可能长的角度提供序列中的二维图像,优选地在舒张期结束时,或心跳周期内的任何其他时刻。
替代地或附加地,图像上的左心室区域在舒张期结束时尽可能大。特别是在基于三维模型提供二维图像的情况下,该选项是可行的。否则,医学成像系统130和被观察对象相对于彼此要恰当定位。在其他示例中,下面讨论的各种步骤在三维模型上执行。在后者的情况下,接收的图像是三维图像。
在步骤206中,在图像数据流或其他序列的图像中识别心脏190的壁轮廓。壁轮廓可以通过边缘检测子单元122或另一检测子单元来检测,包括接收指示轮廓的手动输入,其被设置为检测适用的壁。图3示出了心脏190的左心室的示意图300中的壁轮廓310。可以沿着心脏190的内侧和外侧之间的心肌中部、沿着心内膜、沿着心脏190的外壁、心脏组织中的其他感兴趣生理点或其组合来定义或识别壁。
在步骤208中,如此识别的壁被分割,这可以通过分割子单元124来完成。可以基于心内膜或心肌的已建立的定义区域,根据预定的方案分割壁。替代地或附加地,基于每个段的大小来定义段。在这种情况下,定义的段的总数可能取决于心脏的大小。
在步骤210中,壁点312在壁轮廓310上被定义为矢量点。壁点312可以被放置在相对于已知区段的预定位置上,例如,已经在心脏病学中明确定义的心肌或心内膜的区段。替代地或附加地,壁点312被定义为相对于彼此的预定距离。随后,如果已经定义了段,则壁点312被分配给在其中定义了它们的段。
在步骤212中,在示意图300内定义运动参考320。运动参考320可以是直线或曲线——曲线——、单个点、有限的点集、其他点或其组合。在这种情况下,曲线或直线可以解释为一个点集。点集中的点,无论是有限的还是无限的,都可以提供在曲线或直线上。
优选地,如示意图300所提供的,相对于心脏190的顶点和/或基部设置运动参考320。更具体地,运动参考320被设置在基部和顶点之间距离的三分之一至三分之二之间,更具体地,被设置在基部和顶点之间距离的一半至三分之二之间,例如基于舒张期结束时顶点和基部之间的距离。运动参考可以与生理方向一致,如顶点和基部之间的线或垂直于顶点和基部的线。
在图3中,运动参考320被描绘为参考线322,其上定义有参考点324。注意,运动参考320可以随时间变化,因为它优选地相对于在心跳周期期间心脏上的跳动的点来定义。这样的点可以是基部和顶点。
可选地,如果运动参考320不是单个点,而是作为参考点324的有限集合或参考线322或参考直线322,则特定壁点312可以被分配给特定参考点324或参考直线322或参考曲线322上的特定点。在步骤214中,检查是否已经评估了接收到的所有适用的图像或帧。可以注意到,并不是每个接收到的图像都必须被评估——例如,也可以评估每一个第二个图像。如果没有评估所有图像,则过程跳回到步骤206。
如果已经如上所述处理了所有图像,则过程继续到步骤218,在步骤218中,通过确定心动周期边界,可以根据每个心跳对图像进行分组。心动周期边界可以基于左心室的大小、基于电信号(心电图)、其他或其组合来确定。
在步骤220中,例如通过跟踪运动检测子单元126来获得壁点随时间的运动。如果图像数据按周期分组,则这可以在每个周期的基础上进行。如果不是,则在特定的时间间隔内,在全部或至少部分的采集的流上,跟踪每个壁点的运动。如果在每个周期的基础上跟踪运动,步骤222检查是否已经跟踪了所有周期。如果运动不是每个周期的,则步骤222可能是多余的。
基于点的运动,矢量计算子单元128在步骤224中计算每个壁点312或每个段的运动矢量。运动矢量可以指示基于特定参考在特定方向上的运动的度量,如心动周期的两个阶段之间的最大位移或差分位移。在特定情况下,可以提供运动测量的平均值作为一个段的一般测量,如最大值的平均值或最小值的平均值。位移的测量可以特别针对一个心跳周期来确定。可选地,运动矢量可以指示壁点运动的特定矢量分量的特定分量,该特定矢量分量平行于从壁点312到运动参考320、从收缩期结束到舒张期结束的直线。在后一种情况下,运动矢量或参考运动矢量分量总是指向运动参考或远离运动参考。
图4A示出了舒张期结束时的第一壁轮廓310和收缩期结束时的第二壁轮廓410。图4A还示出了第一运动矢量422、第二运动矢量424、第三运动矢量426、第四运动矢量428、第五运动矢量430和第六运动矢量432。在这种情况下,运动矢量指示每个壁点312在舒张期结束和收缩期结束之间的运动。图4A所示的情况与健康的心脏有关。在另一个实施方式中,运动矢量可以被确定为从壁点312到参考320(特别是从舒张期结束时的参考点324)的第一矢量和从壁点312到参考320(特别是从收缩期结束时的参考点324)的第二矢量之间的差矢量。
因此,基于壁点312在舒张期结束时的第一位置、壁点312在收缩期结束时的第一位置和参考320的位置(可选地基于在其上提供的参考点324的位置)来确定运动矢量。在一个示例中,每个壁点312与特定的、可选地唯一的参考点324相关联。这些参考点可以位于参考曲线322上。此外,这些参考点可以间隔开。
图4B示出了另一种情况。图4B示出了舒张期结束时的第一壁轮廓310’和收缩期结束时的第二壁轮廓410’。图4B还示出了第一运动矢量422’、第二运动矢量424’、第三运动矢量426’、第四运动矢量428’、第五运动矢量430’和第六运动矢量432’。同样在这种情况下,运动矢量指示每个壁点312在舒张期结束和收缩期结束之间的运动。图4B所示的情况涉及在图4B右侧心尖和下壁附近的心肌组织有问题的心脏。
基于所确定的运动矢量,在步骤226中确定运动矢量的平行于适用的壁点312和运动参考320之间的直线的分量。基于在步骤224中确定的运动矢量,可以在步骤228中确定每个段的运动矢量。这种段运动矢量可以是与适用段相关联的壁点312的壁点运动矢量的平均值或中值。基于段运动矢量或基于在步骤226中确定的矢量分量或基于两者,在步骤230中确定段运动矢量的平行于适用段和运动参考320之间的直线的段运动矢量分量。在步骤232中,检查所有段和/或所有壁点312是否已经由矢量计算子单元128评估。如果没有,则过程分支返回到步骤224。
如果矢量计算子单元128已经评估了所有段和/或所有壁点312,则可以在步骤234中确定一个周期期间心脏190或壁轮廓310的整体运动。在心跳周期中,心脏可能在体内从左向右移动,然后返回。可以基于壁点312的移动所提供的信息,例如借助于对壁点312的跟踪,并且例如通过确定所确定的运动矢量的全部或至少一部分的共同运动分量随时间的变化,来检测该移动。可以应用该全局移动来获得运动参考320的运动,或者反过来,可以从图像中直接检测运动参考320的运动,并将其视为表示整体运动。整体运动以及运动参考320的运动—一条或多条曲线、其部分、其上的点、其他点或其组合—直线可以在它们的整个矢量值中被考虑,或者它们可以被限制为平行于生理方向(例如相对于心脏方向的纵向或横向)的运动矢量的分量。基于在心跳周期期间确定的心脏190的整体运动,可以在步骤236中确定补偿因子作为调整数据,以调整在心脏190的运动周期内壁点312或段相对于运动参考320的运动。
在步骤238中,一个周期期间,可以使用调整数据,针对心脏190的运动来调整上面讨论的运动数据值中的所有、一些、一个或没有运动数据值(矢量、矢量分量和运动路径),从而产生调整后的运动数据值。
在步骤240中,获得运动值的参考值。参考值可以涉及下列至少一个:运动轨迹、壁点运动矢量、相对于运动参考的壁点运动矢量分量、段运动矢量、相对于运动参考的段运动矢量分量、其他运动值或其组合。
在步骤242中,基于上述确定和计算的参考运动值和运动值,计算运动偏差值。偏差值可以是下列至少一个:运动矢量值和参考运动值之比或其倒数、运动矢量值和参考运动值之间的差值、运动矢量值和参考运动值之间的差值的绝对值、一方面的运动矢量值和参考运动值之间的差值与另一方面的参考运动值之比或其倒数、运动矢量值和参考运动值之间的差值的绝对值与参考运动值之比或其倒数、其他或其组合。
利用经处理的运动数据、计算出的矢量和矢量分量数据以及已经计算出的偏差数据,可以在电子显示屏144上显示可用数据的至少一部分。在步骤244中,显示偏差数据。图5示出了一个靶心图500。靶心图500基于靶心结构510。靶心结构510被分成多个扇区,该扇区的数量通常在16至18之间,它可以在不同的模式之间变化,并且根据指南,每个扇区可以代表心脏190的心内膜的一段。可选地,靶心结构的每个扇区表示如上所述的壁轮廓310的一个扇区。
在靶心结构510中,偏差数据可以通过不同的颜色、不同的颜色强度、不同的阴影、不同的限制(hedgin))、其他或其组合来绘制,如图5所示。替代地或附加地,偏差数据通过字母数据和数字数据来显示。在靶心图500中,示出了第一偏离部分522、第二偏离部分524和第三偏离部分526。偏差可以是正的,也可以是负的,这取决于偏差数据是如何计算的。在图5中,与参考数据相比,没有特定数据的区域的运动数据没有偏差或偏差很小。
在步骤246中,朝向运动参考320的矢量分量数据显示在靶心图500中。如图5所示,这个矢量数据通过字母数据或数字数据的方式显示。替代地或附加地,该矢量数据可以通过颜色数据、强度数据或其他方式来显示,如上面结合步骤244所讨论的。
在靶心图500中,偏差数据被显示,使得扇区的不同部分可以提供不同的数据。并且矢量分量数据显示在靶心的每个扇区。如果同一实体的多个数据值可用于一个扇区,则可能需要进一步处理数据,例如取扇区内所有数据点的值的中值或平均值。在一个替代方案中,偏差数据可以基于每个扇区来显示。在替代方案中,矢量数据可以以每个扇区的多个数据值来显示。
在已经处理和显示至少一部分这样讨论的数据之后,该过程在步骤248结束。
各个方面和实施例涉及以下编号的示例。
Claims (34)
1.一种在电子计算设备中确定心脏壁的运动的方法,所述方法包括:
接收以固定时间间隔在时间上连续采集的图像序列,所述图像包括心脏的至少一部分的、具有心脏壁的图像数据,
在每个图像中识别所述心脏壁的至少一部分;
在每个图像中定义所识别的心脏壁上的至少一个矢量点;
基于所述序列随时间跟踪所述矢量点;
接收与所述心脏相关的第一参考;
基于对所述矢量点和所述第一参考的跟踪,来确定表示所述矢量点相对于所述参考的运动的运动矢量。
2.根据权利要求1所述的方法,所述方法还包括:
识别所述心脏壁上的矢量点集合;
将所述矢量点集合的第一子集分配给所述心脏壁的第一区域;
随时间跟踪所述第一子集中的所述矢量点;
基于对所述第一子集中的所述矢量点的所述跟踪,确定所述第一区域的第一区域运动矢量,所述第一区域运动矢量指示所述第一区域相对于所述第一参考的运动。
3.根据权利要求2所述的方法,进一步包括:
将所述矢量点集合的第二子集分配给所述心脏壁的第二区域;
随时间跟踪所述第二子集中的所述矢量点;
基于对所述第二子集中的所述矢量点的所述跟踪,确定所述第二区域的第二区域运动矢量,所述第二区域运动矢量指示所述第二区域相对于第二参考的运动。
4.根据权利要求3所述的方法,其中,所述第一参考不同于所述第二参考。
5.根据权利要求4所述的方法,其中,所述第一参考与所述第二参考间隔开。
6.根据前述权利要求中任一项所述的方法,其中,所述运动矢量被确定为在舒张期结束时从所述第一矢量点到所述参考的第一矢量和在收缩期结束时从所述第一矢量点到所述参考的第二矢量之间的差矢量。
7.根据前述权利要求中任一项所述的方法,还包括:
在每个图像中定义所述识别的心脏壁上的多个矢量点;
对于每个矢量点,接收与所述心脏相关的参考;
基于所述序列,对于每个矢量点,相对于与该矢量点相应的参考,随时间跟踪该矢量点;
基于对所述多个矢量点的跟踪,对于每个矢量点,确定表示该矢量点相对于相应的参考的运动的运动矢量。
8.根据权利要求7所述的方法,其中,所述参考在曲线上提供。
9.根据权利要求2至8中任一项在引用权利要求2时所述的方法,还包括:
基于对所述多个矢量点的至少一部分的跟踪,确定所述心脏随时间的运动;
基于所确定的所述心脏的运动来调整所述第一区域运动矢量,以补偿所述心脏的运动。
10.根据权利要求9所述的方法,其中,所述确定所述心脏的运动基于确定所述矢量点的至少一部分的运动矢量中的公共分量。
11.根据前述权利要求中任一项所述的方法,进一步包括:
基于随时间对至少一个矢量点或所述参考的跟踪,确定指示所述参考的运动的参考运动矢量;以及
基于所确定的参考运动矢量调整所述运动矢量。
12.根据权利要求11所述的方法,所述方法还包括:
确定所述参考运动矢量的方向矢量分量,所述方向矢量分量平行于与所述心脏和包括所述心脏的身体中的至少一个的几何形状相关的生理方向,可选地平行于相对于所述心脏的方向的纵向和横向中的至少之一;以及
基于所确定的方向矢量分量调整所述运动矢量。
13.根据前述权利要求中任一项所述的方法,还包括确定所述运动矢量的向内矢量分量,所述向内矢量分量是平行于至少一个第一参考点和所述矢量点之间的线的运动矢量分量。
14.根据前述权利要求中任一项所述的方法,其中,所述运动矢量基于所述矢量点相对于所述第一参考的最大位移。
15.根据权利要求14所述的方法,其中,所述运动矢量基于在单个心动周期期间所述至少一个矢量点相对于所述第一参考的位移极值,所述极值是最大值、最小值和基于最大值和最小值的导出值中的至少一个,优选地是所述最大值和所述最小值的和或差。
16.根据前述权利要求中任一项所述的方法,其中,所述第一参考位于基部附近。
17.根据前述权利要求中任一项所述的方法,其中,所述参考包括以下至少一项:
单个参考点;
多个参考点;
开放曲线;
闭合曲线。
18.根据权利要求17所述的方法,其中,在所识别的心脏壁上确定多个矢量点,所述参考包括多个参考点,其中所述壁上的每个矢量点对应于一个参考点。
19.根据权利要求17或18所述的方法,其中,所述参考点在曲线上提供。
20.根据前述权利要求中任一项所述的方法,还包括基于所述跟踪确定所述矢量点的运动路径,其中,所述运动矢量的确定进一步基于所述运动路径。
21.根据前述权利要求中任一项所述的方法,其中,所述矢量点至少在第一方向上以第一预定距离间隔开。
22.根据权利要求21所述的方法,其中,所述矢量点在第二方向上以第二预定距离间隔,其中所述第二方向通常垂直于所述第一方向。
23.根据前述权利要求中任一项所述的方法,其中,所述矢量点定义在以下至少一个上:
心肌;以及
心内膜。
24.根据前述权利要求中任一项所述的方法,还包括提供用于在电子显示器上显示的显示运动值,所述显示运动值指示所述矢量点相对于所述参考的运动量。
25.根据前述权利要求中任一项所述的方法,还包括:
接收参考运动值,所述参考运动值指示所述矢量点相对于所述参考的参考运动值;
基于所述运动矢量导出运动矢量值;
计算运动偏差值,所述运动偏差值指示所述运动矢量值和所述参考运动值之间的差;以及
提供用于在电子显示器上显示的所述运动偏差值。
26.根据权利要求25所述的方法,其中,所述运动偏差值为下述中的至少一个:
所述运动矢量值与所述参考运动值之间的比值或比值的倒数;
所述运动矢量值和所述参考运动值之间的差值;
所述运动矢量值与所述参考运动值之间的差值的绝对值;
所述运动矢量值与所述参考运动值之间的所述差值与所述参考运动值之间的比值或比值的倒数;
所述运动矢量值和所述参考运动值之间的所述差值的绝对值与所述参考运动值的比值或的比值倒数。
27.根据权利要求23至26中任一项所述的方法,其中,所述显示运动值被设置为作为数值来显示。
28.根据权利要求23至27中任一项所述的方法,其中,所述运动偏差值被设置为作为颜色值来显示。
29.根据权利要求23至28中任一项在引用权利要求3时所述的方法,还包括:
基于所述第一区域运动矢量确定第一区域运动显示值;
基于所述第二区域运动矢量确定第二区域运动显示值;以及
提供用于在电子显示器上显示的所述第一区域运动显示值和第二区域运动显示值。
30.根据权利要求29所述的方法,其中,所述第一区域运动显示值被设置为在所述电子显示器上的第一位置上显示,并且所述第二区域运动显示值被设置为在所述电子显示器上的第二位置上显示。
31.根据权利要求24至30中任一项在引用权利要求3时所述的方法,还包括:
基于所述第一区域运动矢量和第一参考区域运动值确定第一区域偏差显示值;
基于所述第二区域运动矢量和第二参考区域运动值确定第二区域偏差显示值;以及
提供用于在电子显示器上显示的所述第一区域偏差显示值和所述第二区域偏差显示值。
32.一种计算机程序产品,包括计算机可执行代码,当所述代码被加载到处理单元或连接到所述处理单元的存储器中时,使计算机的所述处理单元执行前述权利要求中任一项所述的方法。
33.一种非暂时性计算机可读介质,包括根据权利要求32所述的计算机程序产品。
34.一种电子计算设备,包括:
数据采集模块,设置为接收以固定时间间隔在时间上连续采集的图像序列,所述图像包括心脏的至少一部分的、具有心脏壁的图像数据,
图像处理模块,设置为:
在每个图像中识别所述心脏壁的至少一部分;
在每个图像中定义所述心脏壁上的至少一个矢量点;
基于所述序列随时间跟踪所述矢量点;
运动计算模块,设置为:
接收与所述心脏相关的第一参考;
基于对所述矢量点和所述参考的所述跟踪,来确定表示所述矢量点相对于所述参考的运动的运动矢量。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
NL2028172A NL2028172B1 (en) | 2021-05-07 | 2021-05-07 | Method of determining a motion of a heart wall |
NL2028172 | 2021-05-07 | ||
PCT/NL2022/050248 WO2022235162A1 (en) | 2021-05-07 | 2022-05-06 | Method of determining a motion of a heart wall |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117202842A true CN117202842A (zh) | 2023-12-08 |
Family
ID=76159981
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202280026693.4A Pending CN117202842A (zh) | 2021-05-07 | 2022-05-06 | 确定心脏壁运动的方法 |
Country Status (4)
Country | Link |
---|---|
EP (1) | EP4333697A1 (zh) |
CN (1) | CN117202842A (zh) |
NL (1) | NL2028172B1 (zh) |
WO (1) | WO2022235162A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116616817B (zh) * | 2023-07-21 | 2023-10-03 | 深圳华声医疗技术股份有限公司 | 超声心率检测方法、装置、超声设备及存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7041061B2 (en) * | 2002-07-19 | 2006-05-09 | Cardiac Pacemakers, Inc. | Method and apparatus for quantification of cardiac wall motion asynchrony |
US20070014452A1 (en) * | 2003-12-01 | 2007-01-18 | Mitta Suresh | Method and system for image processing and assessment of a state of a heart |
US20060074315A1 (en) * | 2004-10-04 | 2006-04-06 | Jianming Liang | Medical diagnostic ultrasound characterization of cardiac motion |
US9814439B2 (en) * | 2005-01-19 | 2017-11-14 | Siemens Medical Solutions Usa, Inc. | Tissue motion comparison display |
-
2021
- 2021-05-07 NL NL2028172A patent/NL2028172B1/en active
-
2022
- 2022-05-06 EP EP22722363.3A patent/EP4333697A1/en active Pending
- 2022-05-06 CN CN202280026693.4A patent/CN117202842A/zh active Pending
- 2022-05-06 WO PCT/NL2022/050248 patent/WO2022235162A1/en active Application Filing
Also Published As
Publication number | Publication date |
---|---|
NL2028172B1 (en) | 2022-11-24 |
WO2022235162A1 (en) | 2022-11-10 |
EP4333697A1 (en) | 2024-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9072489B2 (en) | Medical image diagnostic apparatus and medical image contour extraction processing method | |
US9275266B2 (en) | Apparatus and method for tracking contour of moving object, and apparatus and method for analyzing myocardial motion | |
EP1620827B1 (en) | Non-invasive left ventricular volume determination | |
US7454043B2 (en) | Image processing unit and method of associating stored images with current images | |
US8659603B2 (en) | System and method for center point trajectory mapping | |
WO2017206023A1 (zh) | 一种心脏容积识别分析系统和方法 | |
US8577441B2 (en) | System and method for image based physiological monitoring of cardiovascular function | |
CN111095429A (zh) | 心脏病确定方法和系统 | |
JP5919296B2 (ja) | イメージレジストレーション装置 | |
US8487933B2 (en) | System and method for multi-segment center point trajectory mapping | |
CN112074866B (zh) | 4d mr图像数据中的流分析 | |
JP2007167656A (ja) | 対象の運動の解析方法および断層撮影装置 | |
KR20090098839A (ko) | 메디컬 이미징 시스템 | |
US20220277420A1 (en) | Mapping binned medical data | |
US20120008833A1 (en) | System and method for center curve displacement mapping | |
CN110636798A (zh) | 用于生理功能参数确定的方法和装置 | |
CN105877743B (zh) | 使用对准的导管图像对基于阻抗的位置系统的调校 | |
CN117202842A (zh) | 确定心脏壁运动的方法 | |
EP3554341B1 (en) | Retrospective gating of mri | |
WO2018036893A1 (en) | Image processing apparatus and method for segmenting a region of interest | |
US20220117507A1 (en) | System and method for processing measurement data from electrocardiogram electrodes | |
CN116194045A (zh) | 用于提供二次医学成像源的方法 | |
CN114173649A (zh) | 使用功能信息进行的血管配准 |
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 |