CN104271038B - 一种基于电流密度变化评估心肌损伤的方法和设备 - Google Patents

一种基于电流密度变化评估心肌损伤的方法和设备 Download PDF

Info

Publication number
CN104271038B
CN104271038B CN201280072777.8A CN201280072777A CN104271038B CN 104271038 B CN104271038 B CN 104271038B CN 201280072777 A CN201280072777 A CN 201280072777A CN 104271038 B CN104271038 B CN 104271038B
Authority
CN
China
Prior art keywords
damage
diagnosed
phase
ventricle
ratio
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.)
Expired - Fee Related
Application number
CN201280072777.8A
Other languages
English (en)
Other versions
CN104271038A (zh
Inventor
伊利娅·阿纳托利欧维奇·柴可夫斯基
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Publication of CN104271038A publication Critical patent/CN104271038A/zh
Application granted granted Critical
Publication of CN104271038B publication Critical patent/CN104271038B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/242—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents
    • A61B5/243—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents specially adapted for magnetocardiographic [MCG] signals
    • 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 
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Cardiology (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明涉及心脏病学和用于诊断缺血性心肌损伤的方法。进行心磁图检查，在ST‑T间期和QRS波群的4个子间期中重建电流密度矢量映射，计算定量诊断指标。特征在于，对于所述时间间期，得到了所有矢量的总长度(总电流)，瞬时映射的自相关系数以及其与在T波峰值上的映射的相关性，计算出了这些曲线的多个定量指标(曲线下的面积，时间间期，它们的比值等)，并将它们的值的范围划分到3个时间间隔。结果是，根据以下规则诊断心室损伤是否存在——如果某个定量指标在所述3个时间间期中的一个的范围内或如果各定量指标的分值小于或等于7(8‑16,17和更大)，诊断为损伤不存在(轻微存在，明显存在)。

Description

一种基于电流密度变化评估心肌损伤的方法和设备
技术领域
本发明涉及医学，即心脏病学领域，并能用于电生理学损伤和其他心肌损伤的诊断。当前，大多数心脏疾病的可靠诊断仍然是当前存在的临床问题。
背景技术
例如，对于通过广泛使用的心脏电活动的分析-ECG方法进行的心肌缺血的诊断，存在几个已知的诊断标准可以使用，特征为在心动周期的不同时刻心脏EMF矢量的大小和方向。例如，在矢量-ECG方法[Non-invasive diagnosis of coronary artery diseaseusing cardiogoniometry performed at rest,W.Schupbach,B.Emese,P.Loretan et al,Swiss Medical Weekly,2008；138(15-16):230-238,www.smw.ch]中，使用了表征重新极化和去极化过程的13个大小和角度指标来进行缺血性心脏疾病(IHD)的诊断。
但是该方法具有以下几个明显的局限：
1)只分析了仅反映心肌整体情况的矢量和；结果很多局部的电过程损害没有检测到。
2)使用的信息仅仅是关于身体表面点之间电势的差异，而没有心肌内部电流密度分布(CDD)的信息，因此该信息由于心脏周围身体组织的电性能的不均一性和各向异性而被曲解。
因此，ECG方法的诊断准确性不足以用于IHD的临床诊断。因此，需要进一步开发无创的，完全安全的方法，该方法对病人无干扰，没有禁忌，并因此能在不同的临床状况下重复多次。即心磁图(MCG)是一种完全无创的并且安全的，能提供额外信息的方法。
用于心脏病理学诊断的MCG方法的应用前景是基于任何其他方法不能确定心肌中局部电流的方向和大小。只有MCG，能检测心脏电流产生的磁场，允许心脏中心脏电流密度分布的直接重建。
但是对于CDD重建，需要解决，存在于通过测得的磁场分布来计算磁场源(电流)中存在的所谓的反演静磁场问题(IMP)。
MCG提供能作为心脏中异常电生理过程标志的新的有用信息的能力特别涵盖在以下文献中：
1.WO 0217769,Ischemia identification,quantification and partiallocalization in MCG,A.Bakharev,Cardiomag Imaging Inc,USA,2002.
2.WO 002313,Machine learning for classification ofmagnetocardiograms,Sternickel K.B,Cardiomag Imaging Inc,USA,2005.3.ChaikovskyI.Magnetocardiography in unshielded location in coronary artery diseasedetection using computerized classification of current density vectors maps,Dr.Med.dissertation,University Duisburg-Essen,Germany,2006.
4.UA 90701.Method for estimation of the abnormality rate ofelectrical processes in heart's ventricles,IC NASU,MSEC ITiC,Chaikovsky I.,Budnyk M.,Bull.No.10,2010.
5.UA a 2011 07571,Method and device for estimation of the myocardialischemia degree on the base of the current density distribution mapsanalysis,I.Chaikovsky,M.Budnyk,16.06.2011.
6.PCT/UA2011/000121,Method and device for evaluation of myocardialischemia based on current density maps,I.Chaykovskyy,M.Budnyk,02.12.2011.
专利WO 0023 13考虑了使用神经网络的基于自动图像识别技术之一(直接基于内核的自组织映射)的诊断方法。就该识别而言，磁性映射对正常的或病理的都能适用。
该方法的局限在于，其不能使用静磁学反演问题(IPM)的方案，并且其只能适用于两类映射(正常的和病理的)而不能进一步细化。
在I.Chaikovsky的博士论文“Magnetocardiography in unshielded locationin coronary artery disease detection using computerized classification ofcurrent density vectors maps”中，分析了诊断方法，重建了从J点(QRS波群终点)到T波终点的CDD映射，基于最大电流密度矢量的方向评估，将CDD映射分成5个级别——从正常级别到最高病理级别。该方法的局限在于其仅分析了心室再极化并且该映射仅分了少量的几个级别。
在原型方法UA 90701技术中考虑了通过以下手段对心脏心室电过程的异常率进行评估：MCG-映射，从QRS波群起始点到T波终点分析CDD映射，测定局部解剖学参数映射上的高电流密度区域和电流漩涡，计算它们相对于正常的准偶极子映射(quasi-dipole map)的差异等级以及其在设置中的平均值。
方法的特点在于其考虑了具有低电流密度的区域，并基于此通过三进制标度——低、中、高来确定区域不均匀性(regional non-homogeneity degree，RNH)的程度。然后计算大小比(magnitude ratio)，即在R和T波的波峰上的最大电流密度矢量的R/T，并基于此通过三进制标度确定整体不均匀性(GNH)的程度：低——4<R/T<6,中——6<R/T<8,高——R/T>8。
结果，通过5个水平的标度对异常等级给出结论——低RNH等级和GNH时为“正常”，低RNH等级和中GNH等级以及中RNH等级和低GNH等级时为“低异常水平”，低RNH等级和高GNH等级以及高RNH等级和低GNH等级，或中RNH等级和中GNH等级时为“中异常水平”，中RNH等级或高GNH等级以及高RNH等级或中GNH等级时为“高异常水平”，高RNH等级和高GNH等级时为“极高异常水平”。
因此，原型使得确定由于心脏心室的电过程损伤引起的心肌CDD偏差相对于正常状态下的异常等级成为可能。然而，其也具有缺陷，因为指标是基于某一时刻——在R波和T波的波峰上——的电流矢量，而在不是基于心动周期的平均电流矢量计算的。平均指数对于生理学参数的短期干扰和人际分布是不敏感的。
专利WO 0217769考虑了IHD诊断的相关方法，其中根据在ST-段过程中在有效偶极子的垂直坐标(深度)存在变化的情况下，将IHD分为4个严重等级。据此，如果偶极子深度的变化仅发生在1/4的ST-段上(2/4,3/4或整个ST-段上)，出现第1(第2,第3或第4)IHD等级。在没有IHD的健康人中，偶极子的深度几乎是恒定值。
原型的优势在于：
1)定量参数是基于反演问题的方案，而不是基于磁场映射(MFM)分析的；
2)考虑了这些参数在某些心动周期的间期(cardiocycle interval)中的变化。
该模拟方法的局限在于在偶极近似和磁近似中解决了所述反演问题(IP)，但没有在电流偶极子中解决所述反演问题，这不符合心脏的真实电生理过程。简化了基于偶极子模型的分析，这就是在额状面中使用CDD映射形式的更复杂的模型表示(representations)的原因。CDD映射是一个矢量场，其中每个矢量代表给定的点处的电流密度，具有被较低幅值矢量环绕的局部极大值的映射的区域代表心肌的独立的电活性区。也没有分析去极化过程。
原型方法PCT/UA2011/000121(UA a 2011 07571)涵盖了基于电流密度CDD映射评估心肌缺血等级的方法，并用于电生理学紊乱和其他心肌损伤的诊断。根据该方法，进行MCG检查；重建额状面的CDD映射，并计算在ST-T间期的和QRS波群的每个映射的电流密度矢量长度。该方法的特征在于，对于ST-T和QRS的4个子区间，分配三个角度矢量方向分区——正常，中间和异常，4个子区间代表，室间隔1,左心室的前壁和心尖部2,左心室侧壁3,基底心肌区域4的去极化。计算定向在正常(LNORM),异常(LPATH)和中间(LI)分区的矢量的每个映射的总长度，瞬时异常指数Al＝100×LNORM/(LNORM+LPATH+1/2×LI)以及它们对于每个Alave间期的平均值。结果是，如果Alave落入间期70<Alave<100或40≤Alave<70(0≤Alave<40)，诊断1-4区不存在缺血或存在缺血(轻微，明显)。设备包括至少一个电子模块，进行至少一个计算阶段并能包括嵌入式微控制器，或进行至少一个使用电脑软件的阶段。
原型方法的优点是使用CDD映射而不是MFM，使用平均相对异常指数Alave代替瞬时指数，并将表征去极化过程的QRS波群分配到4个子间期。所有的这些使得能分析心室激发的每个阶段，并因此鉴定四个心肌区域中的一个缺血的等级(ischemization degree)。
原型方法的缺陷在于只考虑电流密度矢量的大小和方向，而没有考虑映射的其他定量参数。另一个缺陷是，应用的是平均指数，这不能考虑其随时间的变化。
因此，提出的用于诊断缺血性心肌损伤的方法相比于所述原型的优势在于：
1)不使用电流密度矢量的大小-角度参数，而是使用CDD映射的其他整体定量指标(总电流,相关性,自相关性)，它们大体上概括了CDD映射的总体变化。
2)另外使用概括的评分指数，该指数为不同等级(轻微，明显)心室损伤不存在或存在时的分数的总和。
3)记述基于CDD映射计算的定量指数随时间的变化。
发明内容
本发明的主要任务是，通过引入新的额外指标，使用这些参数的根据经验限定的范围，引入额外的操作来计算这些参数的数值，以及用于确定心肌损伤是否存在以及存在时的严重等级(轻微或明显)的诊断规则，来改进诊断方法。
为了使用提出的方法解决上述任务，记录磁场映射，重建CDD映射，以电流矢量的长度之和计算总电流映射，确定电流映射与T波顶部映射之间的相关系数和自相关系数，对于心室去极化和再极化的每个阶段建立这些参数的变化曲线，分别计算这些阶段的每个阶段中的描述这些曲线的某些间期的形状和持续时间的数值，基于这些数值得出结论，使用具有电子模块的设备，包括使用微处理器或电脑程序。
在基于电流密度分布(CDD)映射指数中的电流密度变化来评估心肌损伤的方法中给定的任务是通过以下步骤解决的：
1)进行心磁图检查；
2)在额状面重建电流密度分布映射，并计算整个心动周期中每个映射的电流密度矢量大小；
3)计算整个心动周期中每个映射的总电流；
4)绘制ST-T间期中的总电流时间变化曲线；
5)基于第4)项的总电流变化曲线，计算ST段持续时间(J点到T波起始点)；
6)计算第5)项的ST-段持续时间与整个ST-T间期持续时间的比值；
7)根据以下规则诊断心室损伤是否存在——如果第6)项所述间期的比值小于25％(25-35％，大于35％)，诊断为损伤不存在(轻微，明显)；
8)计算在代表整个心室的去极化的QRS波群(Q点到J点)中的总电流曲线；
9)根据以下规则诊断心室损伤是否存在——如果第8)项所述QRS波群的持续时间小于0.11s(0.11-0.14s,大于0.14s)，诊断为损伤不存在(轻微，明显)；
10)计算QRS波群的4个时间子间期的总电流曲线,所述4个时间子间期代表，室间隔(1),左心室LV的前壁和心尖部(2),左心室侧壁(3),基底心肌区域(4)的去极化；
11)计算第10项的4个子间期与电生理学已知的“适当”持续时间值的比值；
12)如果4个子间期中的一个的比值范围在0.8到1.2(0.7-0.8和1.2-1.5,低于0.7且高于15)，则诊断第10项中所述的心肌区域(1)-(4)是否存在(轻微，明显)损伤；
13)计算第11)项中计算得到的4个子间期比值的算术平均值；
14)根据以下规则诊断是否存在心室损伤——如果第13)项所述的算术平均值范围在0.8到1.2(0.7-0.8和1.2-1.5,低于0.7且高于15)，诊断为损伤不存在(轻微，明显)；
15)如果QRS波群中存在负波，计数它们的数量；
16)根据以下规则诊断心室损伤是否存在——如果QRS波群中负波的数值小于1(2-3,大于3)，诊断为损伤不存在(轻微，明显)；
17)计算在ST-T间期中瞬时映射与T波峰值上获得的映射的相关系数的相关性；
18)基于第17)项的相关曲线，计算ST段持续时间(J点到T波起始点)；
19)计算第18)项的ST段持续时间与整个ST-T间期持续时间的比值；
20)根据以下规则诊断心室损伤是否存在——如果第19)项的间期比值小于25％(25-35％,大于35％)，诊断为损伤不存在(轻微，明显)；
21)计算T波峰值的最大相关系数(定义其等于1)与其在ST-T间期第4个点处(J点之后40ms)的相关系数值之间的差值；
22)根据以下规则诊断心室损伤是否存在——如果第21)项所述的差值小于30％(30-50％，大于50％)，诊断为损伤不存在(轻微，明显)；
23)计算ST-T间期的平均相关系数；
24)根据以下规则诊断心室损伤是否存在——如果第23)项所述平均相关系数高于75％,(60-75％,低于60％)，诊断为损伤不存在(轻微，明显)；
25)计算在QRS波群中瞬时映射与R波峰值上的映射的相关系数的时间相关性；
26)根据以下规则诊断心室损伤是否存在——如果第25)项所述的相关曲线，在曲线迅速上升/下降时刻，含有2个(1个,0个)特征拐点，诊断为损伤不存在(轻微，明显)；
27)计算在QRS波群中瞬时映射的自相关系数的时间相关性；
28)根据以下规则诊断心室损伤是否存在——如果第27)项所述的自相关曲线含有2个(1个,0个)明显的局部最小值，诊断为损伤不存在(轻微，明显)；
29)对诊断“损伤不存在”给0分；“损伤轻微存在(明显存在)”给1(2)分；
30)根据以下规则诊断心室损伤是否存在——如果第7,9,12,14,16,18,20,22,24,26,28项所述的各个指标根据第29项计算的分数之和小于或等于7(8-16,17或更大)，诊断为损伤不存在(轻微，明显)；
31)使用设备对第3-30项中至少一项进行计算，所述设备包括至少一个通过嵌入式微处理器控制的电子单元；
32)第31)项所述设备的使用方式为通过外部计算机对其进行控制。
提出的方法相比于现有技术的新颖性在于：
1)使用基于IPM的新的定量CDD指标(总电流,相关性，自相关性)，以获得更可靠的对心肌损伤的存在和严重等级的评估；
2)使用不同类型的定量指标，如时间，电流，相关系数；
3)组合不同种类的指标用于基于广义的评分标准形成诊断结论；
4)计述了基于CDD映射计算的定量指数随时间的变化。
技术效果在于：
1)提高了损伤检测的可靠性；
2)能为心肌损伤的存在和严重等级提供更准确的评估；
3)能确定损伤的位置。
附图说明
图1为方法和设备的示意图。
图2为ST段(J点到T波起始点)的总电流曲线。垂直标记指的是T波出现的起始点:A)损伤不存在；B)损伤明显存在。
图3为在代表整个心室去极化的QRS波群(Q点到J点)中的总电流曲线。垂直标记指的对应T波峰值的点:A)损伤不存在；B)损伤明显存在。
图4为在ST-T间期瞬时映射与T波峰值上的映射的相关系数的时间相关曲线:A)损伤不存在；B)损伤明显存在。
图5为在QRS波群中瞬时映射与R波峰值上的映射的相关系数的时间相关曲线:A)损伤不存在；B)损伤明显存在。
图6为在QRS波群中瞬时映射的自相关系数的时间相关曲线:A)损伤不存在；B)损伤明显存在。
具体实施方式
在本发明方法的主要实施方式中，所述方法包含操作顺序，执行算法的各个步骤，各个步骤在执行时使用所提出的方法。所述方法的示意图如图1所示，其中虚线的步骤没有描述提出的方法，仅用于说明其背景。
步骤1.进行MCG检查，对磁场映射集进行初步的信号处理和计算。
提出的方法的具体实施包括，初步进行MCG检查，同时在第二标准引线中进行ECG注册并处理接收的信号。最后一步包括多个操作：MFM映射集的初步的处理(过滤，去除具有脉冲干扰的波群,选择相同类型的心动周期并求平均值)与计算，以及在20×20cm的方形面积内进行CDD映射的重建。这些操作记载在专利UA 77722,UA 77723,UA 20101,UA 21299中并且不是本发明的主题。
步骤2.在额状面重建CDD映射，并计算整个心动周期中的电流密度矢量大小
CDD映射是一个矢量场，其中各矢量放置在规则网格(通常方网格)上。每个矢量反映了局部值(local values),即在给定点处的电流密度。此处使用根据WO/2002/00108[Computer-based method for automatically processing data,especiallymagnetocardiographic data,of biomagnetic fields.S.Romanovich,F.Steinberg,2001]的反演问题解决方法。该方法之前在专利UA 83050,UA 83061,UA 90701和原型专利PCT/UA2011/000121(UAa201107571)中使用过。结果是获得了在额状面的CDD映射，其中各矢量放置在具有2cm间距的10×10的方网格上。
计算电流密度矢量的大小。在主要的实施方式中，计算映射中的100个点的矢量的正投影。使用已知的通过矢量投影来表示该矢量长度的公式计算该矢量的大小。以1ms的时间采样速率连续对所有CDD映射进行计算。
这些操作与UA90701以及原型PCT/UA2011/000121(UAa201107571)中的相似但不是本发明的主题。
步骤3.分析在ST-T间期的总电流。该步骤包括执行以下操作：
3.1计算整个心动周期中每个映射的总电流；
3.2绘制所述ST-T间期的总电流时间变化曲线。在图2上(a-损伤不存在；b-损伤明显存在)示出了在所述ST-T间期(J点到T波起始点)的总电流曲线。垂直标记表示T波的起点，其根据在电流密度迅速增加时刻的拐点来确定。
3.3基于根据第3.1项的总电流变化曲线，计算ST段持续时间(J点到T波起始点)；
3.4计算根据第3.2项的ST-段持续时间与整个ST-T间期持续时间的比值；
3.5根据以下规则诊断心室损伤是否存在——如果根据第3.4项所述间期的比值小于25％(25-35％，大于35％)，诊断为损伤不存在(轻微，明显)；
步骤4.分析QRS波群上的总电流时间特性。该步骤包括执行以下操作：
4.1.计算在代表整个心室的去极化的QRS波群(Q点到J点)中的总电流曲线，参见图3，其中垂直标记表示对应T波峰值的点(a-损伤不存在；b-损伤明显存在)；
4.2在所述QRS波群持续时间小于0.11s(0.11-0.14s,大于0.14s)时，诊断心室损伤是否存在(轻微，明显)；
4.3计算QRS波群的4个时间子间期的总电流曲线,所述4个时间子间期代表，室间隔(1),左心室LV的前壁和心尖部(2),左心室侧壁(3),基底心肌区域(4)的去极化；
4.4计算根据第4.3项的4个子间期的持续时间与电生理学已知的“适当”持续时间值的比值；
4.5如果4个子间期中的一个的比值范围在0.8到1.2(0.7-0.8和1.2-1.5,低于0.7且高于15)，诊断第4.3项所述的心肌区域(1)-(4)是否存在(轻微，明显)损伤；
步骤5.分析QRS波群上的总电流-时间特性。该步骤包括执行以下操作：
5.1.计算第4.4项中计算的4个子间期比值的算术平均值；
5.2.根据以下规则诊断心室损伤是否存在——如果第5.1项所述的算术平均值范围在0.8到1.2(0.7-0.8和1.2-1.5,低于0.7且高于15)，诊断为损伤不存在(轻微，明显)；
5.3.如果QRS波群中存在负波，计数它们的数量；
5.4.根据以下规则诊断心室损伤是否存在——如果QRS波群中负波的数目小于1(2-3,大于3)，诊断为损伤不存在(轻微，明显)；
步骤6.分析在ST-T间期的相关系数。该步骤包括执行以下操作：
6.1.计算在ST-T间期中瞬时映射与T波峰值上的映射的相关系数的时间相关性，参见图4(a-损伤不存在；b-损伤明显存在)；
6.2.基于第6.1项的相关系数曲线，计算ST段持续时间(J点到T波起始点)；
6.3计算第6.2项的ST段持续时间与整个ST-T间期持续时间的比值；
6.4根据以下规则诊断心室损伤是否存在——如果根据第6.3项所述的间期比值小于25％(25-35％,大于35％)，诊断为损伤不存在(轻微，明显)；
6.5计算T波峰值的最大相关系数(定义其等于1)与其在ST-T间期第4个点处(J点之后40ms)的相关系数值之间的差值；
6.6根据以下规则诊断心室损伤是否存在——如果根据第6.5项所述的差值小于30％(30-50％，大于50％)，诊断为损伤不存在(轻微，明显)；
6.7计算ST-T间期的平均相关系数；
6.8根据以下规则诊断心室损伤是否存在——如果根据第6.7项所述的平均相关系数高于75％(60-75％,低于60％)，诊断为损伤不存在(轻微，明显)；
步骤7.分析QRS波群的相关系数和自相关系数。该步骤包括执行以下操作：
7.1计算在QRS波群中瞬时映射与R波峰值上的映射的相关系数的时间相关性，参见图5(a-损伤不存在；b-损伤明显存在)；
7.2根据以下规则诊断心室损伤是否存在——如果第7.1项的相关系数曲线，在该曲线迅速上升或下降时刻，含有2(1,0)个特征拐点，诊断为损伤不存在(轻微，明显)；
7.3根据以下规则诊断心室损伤是否存在——如果根据第7.1项所述的差值小于41％(41-60％,大于60％)，诊断为损伤不存在(轻微，明显)；
7.4计算在QRS波群中瞬时映射的自相关系数的时间相关曲线，参见图6(a-损伤不存在；b-损伤明显存在)；
7.5根据以下规则诊断心室损伤是否存在——如果根据第7.4项的自相关系数曲线含有2(1,0)个明显的局部最小值，诊断为损伤不存在(轻微，明显)；
步骤8.基于广义的评分标准的诊断。该步骤包括执行以下操作：
8.1对诊断“损伤不存在”给0分；“损伤轻微存在(明显存在)”给1(2)分；
8.2根据以下规则诊断心室损伤是否存在——如果根据第3.5,4.2,4.5,5.2,5.4,6.4,6.6,6.8,7.2,7.3,7.5项所述的各个定量指标根据第8.1项计算的分数之和小于或等于7(8-16,17和更大)，诊断为损伤不存在(轻微，明显)；
在主要实施方式中，所有所述操作和根据步骤3-8的计算使用以下设备进行，所述设备包括至少一个通过嵌入式的微处理器控制的电子单元。在另一个实施方式中，所述电子模块中的至少一个模块使用外部计算机控制。
对本发明方法的具体实施方式的详细描述仅用于示例说明的目的。对于具有MCG数据分析经验且通常是心脏病学领域的人而言清楚的是，在实施时可引入一些改变和修饰，例如，将心动周期间期分为其他数目的子间期，使用定量指标间期的其他组合方式，通过一些参数改变对不同诊断的评分标度，增加或减少评分等级的数目，改变数值指标的求平均值方法，或引入其他对于心肌损伤的不存在，存在和严重等级的解释，等。但我们认为，所述改变和修饰以及其他改变和修饰，在与提出的本发明的本质和范围没有明显不同的情况下，都落入本申请的范围之内。

Claims (1)

1.一种基于电流密度变化评估心肌损伤的设备，应用至少一个电子单元，所述电子单元通过嵌入式微处理器或外部计算机控制；其中，所述至少一个电子单元用于接收心磁图数据以及根据所述心磁图数据执行以下步骤：在额状面重建电流密度分布映射，并计算整个心动周期中每个映射的电流密度矢量长度；
其特征在于，所述至少一个电子单元还用于执行以下步骤：
将QRS波群划分为4个子间期，所述4个子间期代表心肌区域(1)-(4)，所述心肌区域(1)-(4)代表室间隔(1)，左心室LV的前壁和心尖部(2)，左心室侧壁(3)，基底心肌区域(4)的去极化，计算所述每个映射的定量指标和评分标准，
计算在整个心动周期中每个映射的总电流，将总电流曲线绘制到ST-T间期，在J点到T波起始点计算ST段持续时间以及其与ST-T间期持续时间的比值，如果所述比值小于25％诊断为心室损伤不存在，如果所述比值为25％到35％诊断为心室损伤轻微存在，如果所述比值大于35％诊断为心室损伤明显存在，
计算在QRS波群中Q点到J点的总电流曲线，如果QRS波群持续时间小于0.11s诊断为心室损伤不存在，如果QRS波群持续时间为0.11至0.14s诊断为心室损伤轻微存在，如果QRS波群持续时间大于0.14s诊断为心室损伤明显存在，
计算QRS波群的4个子间期的总电流曲线的持续时间以及所述4个子间期的总电流曲线的持续时间与适当持续时间的比值，从而得到4个子间期的比值，如果所述4个子间期的比值的某一比值在0.8到1.2的范围诊断为心肌区域(1)-(4)中与所述比值在0.8到1.2的范围对应的子间期损伤不存在，如果所述4个子间期的比值的某一比值在0.7到0.8的范围以及1.2到1.5的范围诊断为心肌区域(1)-(4)中与所述比值在0.7到0.8的范围以及1.2到1.5的范围对应的子间期损伤轻微存在，如果所述4个子间期的比值的某一比值在低于0.7且高于1.5的范围诊断为心肌区域(1)-(4)中与所述比值在低于0.7且高于1.5的范围对应的子间期损伤明显存在；
计算4个比值的算术平均值，如果平均比值为0.8到1.2的范围，诊断为心室损伤不存在，如果平均比值为0.7到0.8的范围以及1.2到1.5的范围诊断为心室损伤轻微存在，如果平均比值为低于0.7且高于1.5的范围，诊断为心室损伤明显存在，
如果QRS波群含有负波，计数所述负波的数目，如果所述负波的数目小于1诊断为心室损伤不存在，如果所述负波的数目为2到3诊断为心室损伤轻微存在，如果所述负波的数目大于3诊断为心室损伤明显存在，
计算在ST-T间期中瞬时映射与T波峰值上获得的映射的相关系数的相关性，ST段持续时间，ST段持续时间与ST-T间期持续时间的比值，如果所述ST段持续时间与ST-T间期持续时间的比值小于25％诊断为心室损伤不存在，如果所述ST段持续时间与ST-T间期持续时间的比值为25到35％诊断为心室损伤轻微存在，如果所述ST段持续时间与ST-T间期持续时间的比值大于35％，诊断为心室损伤明显存在，
计算T波峰值的最大相关系数与其在ST-T间期第4个点处的相关系数值之间的差值，所述最大相关系数等于1，所述第4个点处为J点之后40ms，如果所述差值小于30％诊断为心室损伤不存在，如果所述差值为30％到50％诊断为心室损伤轻微存在，如果所述差值大于50％诊断为心室损伤明显存在，
计算在QRS波群中瞬时映射与R波峰值上映射的相关系数的时间相关曲线和自相关系数的时间相关曲线，如果在所述相关系数的时间相关曲线迅速上升或下降时刻所述相关系数的时间相关曲线含有2个特征拐点，或在所述自相关系数的时间相关曲线中出现2个局部最小值，诊断为心室损伤不存在；如果在所述相关系数的时间相关曲线迅速上升或下降时刻所述相关系数的时间相关曲线含有1个特征拐点，或在所述自相关系数的时间相关曲线中出现1个局部最小值，诊断为心室损伤轻微存在；如果在所述相关系数的时间相关曲线迅速上升或下降时刻所述相关系数的时间相关曲线含有0个特征拐点，或在所述自相关系数的时间相关曲线中出现0个局部最小值，诊断为心室损伤明显存在；
计算ST-T间期的平均相关系数和QRS波群的平均相关系数，如果所述ST-T间期的平均相关系数高于75％或所述QRS波群的平均相关系数低于41％，诊断为心室损伤不存在，如果所述ST-T间期的平均相关系数为60％到75％或所述QRS波群的平均相关系数为41％到60％，诊断为心室损伤轻微存在，如果所述ST-T间期的平均相关系数低于60％或所述QRS波群的平均相关系数高于60％，诊断为心室损伤明显存在，
对诊断“损伤不存在”给分0分，对诊断“损伤轻微存在”给分1分，对诊断“损伤明显存在”给分2分；
如果所有的指标对应的分数之和小于或等于7诊断为心室损伤不存在，如果所有的所述指标的分数之和为8到16诊断为心室损伤轻微存在，如果所有的所述指标的分数之和大于或等于17诊断为心室损伤明显存在。
CN201280072777.8A 2012-07-13 2012-09-07 一种基于电流密度变化评估心肌损伤的方法和设备 Expired - Fee Related CN104271038B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
UAA201208663 2012-07-13
UAA201208663A UA104073C2 (uk) 2012-07-13 2012-07-13 Спосіб оцінки ступеня ушкодження міокарда на основі аналізу змін у часі показників щільності струму
PCT/UA2012/000083 WO2014011135A1 (en) 2012-07-13 2012-09-07 Method and device for evaluation of myocardial damages based on the current density variations

Publications (2)

Publication Number Publication Date
CN104271038A CN104271038A (zh) 2015-01-07
CN104271038B true CN104271038B (zh) 2017-11-28

Family

ID=49916412

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280072777.8A Expired - Fee Related CN104271038B (zh) 2012-07-13 2012-09-07 一种基于电流密度变化评估心肌损伤的方法和设备

Country Status (6)

Country Link
US (1) US9451901B2 (zh)
EP (1) EP2872039B1 (zh)
JP (1) JP5996111B2 (zh)
CN (1) CN104271038B (zh)
UA (1) UA104073C2 (zh)
WO (1) WO2014011135A1 (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104794474B (zh) * 2015-03-26 2017-12-05 同济大学 一种基于心磁图的综合最大电流密度向量特征提取方法
EP3207866A1 (en) * 2016-02-18 2017-08-23 Biomagnetik Park GmbH Vector magnetocardiography method and vector magnetocardiographic system
CN108577825B (zh) * 2018-01-11 2021-06-18 中国科学院上海微系统与信息技术研究所 一种心肌缺血病变位置定位方法及系统、存储介质及终端
CN109512422A (zh) * 2018-09-25 2019-03-26 维灵（杭州）信息技术有限公司 一种ecg波形混乱程度计算方法
US10602940B1 (en) * 2018-11-20 2020-03-31 Genetesis, Inc. Systems, devices, software, and methods for diagnosis of cardiac ischemia and coronary artery disease
EP3685740A1 (en) * 2019-01-28 2020-07-29 Primax Medica Sp. z o.o. Method for determining electrical activity of cardiac muscle
CN110074774B (zh) * 2019-04-28 2022-04-12 漫迪医疗仪器（上海）有限公司 基于心磁图的心脏室间隔异常的分析方法、系统、介质及终端
CN113499075A (zh) * 2021-06-08 2021-10-15 苏州卡迪默克医疗器械有限公司 一种基于心磁信号的心肌缺血评估方法和装置
CN117137492B (zh) * 2023-11-01 2024-02-09 山东大学齐鲁医院 冠状动脉血流异常检测系统、存储介质及终端

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3824990A (en) * 1971-06-30 1974-07-23 Instr For Cardiac Res Inc Method and apparatus for producing sample electrocardiograms
US5685316A (en) 1996-04-08 1997-11-11 Rheo-Graphic Pte Ltd. Non-invasive monitoring of hemodynamic parameters using impedance cardiography
DE10082810D2 (de) 1999-08-28 2002-08-29 Squid Ag Computerbasiertes Verfahren zur automatischen Aufbereitung von Daten biomagnetischer Felder, insbesondere von Magnetokardiographischen Daten
JP3427811B2 (ja) * 2000-03-13 2003-07-22 株式会社日立製作所 生体磁場計測装置
WO2002017769A2 (en) 2000-08-29 2002-03-07 Cardiomag Imaging, Inc. Ischemia identification, quantification and partial localization in mcg
JP3944383B2 (ja) * 2001-11-16 2007-07-11 株式会社日立製作所 心臓磁場計測装置
JP4027746B2 (ja) 2002-08-07 2007-12-26 株式会社日立ハイテクノロジーズ 生体磁場計測装置
CN1552285A (zh) 2003-05-26 2004-12-08 中国科学院物理研究所 一种心磁波组图及其分析方法
US7742806B2 (en) 2003-07-01 2010-06-22 Cardiomag Imaging, Inc. Use of machine learning for classification of magneto cardiograms
UA74466C2 (en) * 2004-02-18 2005-12-15 Budnyk Mykola Mykolaiovych Method for diagnosing ischemic disease of heart
EP1721278A2 (en) 2004-02-27 2006-11-15 Koninklijke Philips Electronics N.V. Method and system to minimize power consumption by using staged life-threatening arrhythmia detection algorithm
US20080033312A1 (en) * 2004-06-01 2008-02-07 Kenji Nakai Cardiac Magnetic Field Diagnostic Apparatus and Evaluating Method of Three-Dimensional Localization of Myocardial Injury
JP3835805B2 (ja) * 2004-09-10 2006-10-18 独立行政法人科学技術振興機構 心臓磁界診断装置およびその作動方法
UA77722C2 (uk) 2004-06-29 2007-01-15 Інститут Кібернетики Ім. В.М.Глушкова Національної Академії Наук України Спосіб магнітокардіографічного обстеження
UA77723C2 (uk) 2004-06-29 2007-01-15 Hlushkov Inst Of Cybernetics O Спосіб магнітокардіографічного картування
KR20060117545A (ko) 2005-05-11 2006-11-17 인하대학교 산학협력단 심전도 신호의 특징추출방법
UA83050C2 (uk) 2006-01-23 2008-06-10 Институт Кибернетики Им. В.М. Глушкова Спосіб оцінки анормальності розподілу струмів у серці
UA83061C2 (uk) 2006-03-16 2008-06-10 Институт кибернетики им. В.М. Глушкова НАН Украины Спосіб оцінки анормальності процесу реполяризації шлуночків серця
UA20101U (en) 2006-06-27 2007-01-15 Ltd Liability Company Kyiv Med Method for processing magnetocardiographic signals
UA21299U (en) 2006-08-15 2007-03-15 V M Hlushkov Inst Of Cyberneti Method for recording magnetocardiogram
UA90701C2 (uk) 2007-07-27 2010-05-25 Международный Научно-Учебный Центр Информационных Технологий И Систем Нан И Мон Украины Спосіб оцінки ступеня анормальності електричних процесів у шлуночках серця
UA100929C2 (ru) 2011-06-16 2013-02-11 Илья Анатольевич Чайковский Способ оценки степени выразительности ишемии миокарда на основе анализа карт распределения плотности тока

Also Published As

Publication number Publication date
EP2872039A1 (en) 2015-05-20
US20150011862A1 (en) 2015-01-08
UA104073C2 (uk) 2013-12-25
US9451901B2 (en) 2016-09-27
JP2015525633A (ja) 2015-09-07
EP2872039B1 (en) 2017-05-17
JP5996111B2 (ja) 2016-09-21
EP2872039A4 (en) 2016-03-09
CN104271038A (zh) 2015-01-07
WO2014011135A1 (en) 2014-01-16

Similar Documents

Publication Publication Date Title
CN104271038B (zh) 一种基于电流密度变化评估心肌损伤的方法和设备
CN105997054B (zh) 一种心电信号预分析的方法
CN107358196A (zh) 一种心搏类型的分类方法、装置及心电仪
CN109117730A (zh) 心电图心房颤动实时判断方法、装置、系统及存储介质
Brecheisen et al. Parameter sensitivity visualization for DTI fiber tracking
CN110074774B (zh) 基于心磁图的心脏室间隔异常的分析方法、系统、介质及终端
CN103961089B (zh) 基于分段直线拟合的窦性心率震荡趋势检测方法
CN105147269B (zh) 一种无创连续血压测量方法
CN102194049A (zh) 诊断辅助设备和冠状动脉分析方法
JP2015512292A (ja) マルチプルイメージデータループを獲得して分析する方法及びシステム
JP4027867B2 (ja) 生体磁場計測装置
US20080004537A1 (en) Method and system for multi-channel biosignal processing
US10241180B2 (en) System and method for mapping and quantifying in-vivo blood flow stasis
WO2012173583A1 (en) Method and device for evaluation of myocardial ischemia based on current density maps
CN110327032A (zh) 一种单导心电信号pqrst波联合精准识别算法
CN107835656A (zh) 心肌兴奋补充/可视装置和心肌兴奋检测装置
CN103142211A (zh) 一种基于极值圆的心脏磁信号处理方法
US8842897B2 (en) Quantification of the image quality for the motion-synchronized capture of a cross-sectional tomographic picture of an object
JP2002028143A5 (zh)
KR101748261B1 (ko) 심전도 데이터를 그래프로 시각화하여 부정맥을 판별하는 방법 및 시스템
Kräuter et al. Automated vortical blood flow-based estimation of mean pulmonary arterial pressure from 4D flow MRI
CN109068984A (zh) 心肌兴奋判定装置
JP2002028144A5 (zh)
Zheng et al. 3D fractal dimension analysis: Prognostic value of right ventricular trabecular complexity in participants with arrhythmogenic cardiomyopathy
CN103377475B (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171128

Termination date: 20200907