一种基于心磁特征确定心脏磁场极大值点的方法

摘要

本发明公开了一种基于心磁特征快速确定心脏磁场极大值点方法，该方法首先在人体胸前或背后心磁测量区域建立坐标系，在上述坐标系中，画一正方形框定人体测量范围，选择相邻两条边的中点及对应顶点构成三角形，并以此三角形的三个顶点为测量点，分别进行心磁信号测量，然后对三边所在方向迭代两次求得对应的驻点和相应极值，对比上述所求极值即可得到极大值位点，上述水平分量的心脏磁场是单一极大值模式。利用本发明可在选定的心磁测量区域中准确测量出信号极值点的位置并给出极值点的心磁信号幅度，便于获得有价值的心磁信息。

权利要求

1.一种确定心脏磁场极大值点的方法，其步骤包括：

1)在人体胸前或背后的心磁测量区域上建立一坐标系；

2)在上述坐标系中，画一正方形框定人体测量范围，选择其相邻两条边的中点以及对应顶点构成三角形；

3)选择上述三角形的三个顶点为测量点，分别得到心磁信号值T1、心磁信号值T2和心磁信号值T3；

4)分别对上述三角形的三边所在方向的目标函数采用正交偏导数求极值法得到驻点M4、M5和M6；

5)驻点M4、M5和M6构成第二三角形，分别对第二三角形的三边所在方向的目标函数采用正交偏导数求极值法得到驻点M7、M8和M9，将驻点M7、M8和M9的坐标代入公式从而得到对应的心磁信号值T7、T8和T9；

6)将上述心磁信号值T7、T8和T9中的最大值作为该心磁测量区域极大值点的心磁信号值，该极大值点的坐标点为心磁测量区域的极大值点位置。

2.如权利要求1所述的确定心脏磁场极大值点的方法，其特征在于，所述坐标系为平面直角坐标系，坐标系以左侧锁骨中线的人体躯干方向为y轴，x轴方向为与y轴垂直方向的经过人体剑突的直线，相交点O为原点。

3.如权利要求2所述的确定心脏磁场极大值点的方法，其特征在于，所述正方形的边长为10～15厘米。

4.如权利要求2所述的确定心脏磁场极大值点的方法，其特征在于，所述正方形的中心位于坐标系y轴+6厘米～+13厘米处。

说明书

一种基于心磁特征确定心脏磁场极大值点的方法

技术领域

[0001] 本发明是关于心脏磁场信号测量技术，具体涉及一种基于心磁特征快速确定心脏磁场极大值点的方法。

背景技术

[0002] 心磁信号比常规的心电信号可以包含更丰富的信息，通过人体心脏磁场的测量进行心脏疾病的诊断，是近年来发展起来的全新的心脏病诊断技术，比常规心电图对某些心脏病的诊测更灵敏，具有无接触、无创的优势，有望在不久的将来用作临床诊断心脏疾病的重要技术手段。

[0003] 在心脏磁场测量过程中，为了能够得到足够的心脏磁场分布信息，进行基于心磁的心脏物理相关研究及进行心脏疾病的临床诊断，需要在人体体表(例如胸前或背后)进行多点(如36点)的心磁测量，如图1所示。利用心磁图仪测得每个点的心磁信号后，在去除了心磁信号的直流项后，将R波峰值大小定义为心磁信号大小，如图2所示的心磁信号大小约为62pT。

[0004] 目前，多点测量点位方法通常采用基于人体结构确定心磁测量覆盖范围，其覆盖范围的选择，是根据人体骨骼、胸围等解剖结构来进行，具体步骤包括：首先在人体体表选定一块与心脏对应的区域作为测量覆盖范围，并将其划分成网格，进而测量格点位置上的心磁信号随时间的变化。该方法选定的覆盖区域的大小、位置，以及网格的划分，对最后所获得的心磁图的有效信息含量及相关分析有很大的影响，由于人体心脏在胸腔中的位置及心动强弱的个体差异较大。如采用保守的大面积覆盖方式测量，极易导致边界点附近的采集信号太小，结果没有实际意义，如测量范围过小，难以覆盖心磁信号的有效区域，致使许多对诊断有潜在价值的心磁信息被遗漏。因此，以人体解剖结构来确定多点测量点位的方法对于实际工作，特别是对于临床应用非常不利，常常存在不合理之处。

发明内容

[0005] 本发明提出一种基于心磁特征快速确定心脏磁场极值点方法，该方法可在心磁测量区域中准确测量出信号极值点的位置并给出极值点的心磁信号幅度，便于获得有价值的心磁信息。本发明提供的心磁测量区域极大值点的测定方法，其步骤包括：

[0006] 1)在人体胸前或背后的心磁测量区域上建立一坐标系；

[0007] 2)在上述坐标系中，画一正方形框定人体测量范围，选择其相邻两条边的中点以及对应顶点构成三角形；

[0008] 3)选择上述三角形的三个顶点为测量点，分别得到心磁信号值T1、心磁信号值T2和心磁信号值T3；

[0009] 4)分别对上述三角形的三边所在方向的目标函数采用正交偏导数求极值法得到驻点M4、M5和M6；

[0010] 5)选择M4M5M6构成的三角形，重复采取步骤4)，得到驻点M7、M8和M9，将驻点M7、M8和M9的坐标代入公式 从而得到对应的心磁信号值T7、T8和T9；

[0011] 6)将上述心磁信号值T7、T8和T9中的最大值作为该心磁测量区域极大值点的心磁信号值，该极大值点的坐标点为心磁测量区域的极大值点位置。

[0012] 所述坐标系为平面直角坐标系，坐标系以左侧锁骨中线的人体躯干方向为y轴，x轴方向为与y轴垂直方向的经过人体剑突的直线，相交点O为原点。

[0013] 所述正方形的中心位于坐标系y轴+6厘米～+13厘米处，边长为10～15厘米。

[0014] 所述正交偏导数求值法采用数值分量的偏微分法。

[0015] 本发明的技术效果是：

[0016] 心磁信号具有x,y,z三个方向，根据测量的z方向的心磁信号值采用正交偏导数方法得到xy方向的心磁信号值，其心脏磁场在xy平面呈现单一极大值模式，故通过本发明所提供的方法，可以快速、方便地在选定的测量区域中获得心磁信号的极大值点。以此极大值点的位置和信号幅度为依据，在进行心磁的多点测量之前，从信号源的本征特性出发，进而合理确定进行多点心磁测量所需覆盖的有效区域，既不会遗漏最强点信号，同时又不必包含没有诊断价值的信号，完全避免根据人体结构确定测量覆盖范围所带来的各种弊端，确保测量的结果有足够的心磁研究和医学诊断价值，为后续相关研究及诊断信息的提取、促进心磁图技术的临床应用提供有利条件。

附图说明

[0017] 图1为现有的心磁图测量区域划分方法示意图；

[0018] 图2为现有的心磁信号波形示意图；

[0019] 图3为本发明测定胸前心磁极大值点过程中各采样点示意图；

[0020] 图4为本发明实施例中得到的心磁信号值T和驻点M；

[0021] 图5为本发明实施例得到的极大值磁信号值T的磁场图，(a)是三维磁场图，(b)是二维等磁图。

具体实施方式

[0022] 下面结合附图和具体实施方式对本发明作进一步详细描述：

[0023] 以确定胸前一心磁测量区域的极大值点的位置与信号幅度为例：

[0024] 1)在胸前心磁测量区域内建立平面直角坐标系；

[0025] 2)如图3所示：本发明实施例采用的平面直角坐标系中，y轴方向大致为人体躯干方向，从下指向上，x轴方向为与y轴垂直的方向，从右指向左；

[0026] 3)在上述坐标系中，选择以过人体剑突的直线为x轴，以左侧锁骨中线为y轴，相交点O为原点；

[0027] 4)在上述坐标系中，在y轴选择点P，OP＝3cm，然后过P(0,3)点做平行于x轴的直线AM1，且AP＝PM1＝6cm。以边长为12cm，画出正方形ABCM1，顶点坐标分别为：A(-6,3)，B(-6,15)，C(6,15)和M1(6,3)；

[0028] 5)以上述正方形ABCM1的边BC的中点M2(0,15)，边AB的中点M3(-6,9)，正方形顶点M1(6,3)为顶点构成三角形，如图3所示；

[0029] 6)以上述三角形M1M2M3的三个顶点为测量点进行心磁信号测量，分别得到心磁信号T1(12pT)、心磁信号值T2(40pT)和心磁信号值T3(25pT)；

[0030] 7)根据心磁信号值T1、T2和T3，采用正交偏导数求极值法得到边M1M2，M2M3以及边M3M1所在方向心磁信号值T4，心磁信号值T5和心磁信号值T6，如图4所示。

[0031] 在数学中，一个多变量的函数的偏导数，就是它关于其中一个变量的导数而保持其他变量恒定。偏导数在工程中是很有用的。在工程中遇到的函数大部分是连续的，或者能保证在所讨论的区域内取得它的极大值或极小值。如果知道可微函数的极大值或极小值一定在区域内，且函数在区域内有仅有一个驻点，则该驻点一定是极大值或极小值。

[0032] 在本实施例中，心磁信号是连续的，且是矢量，分量Bx By Βz相互垂直，因而是正交的。正交偏导数的定义是：对于闭合区域D＝{(x,y)|-6cm≤x≤6cm,3≤y≤15cm}，目标函数  式中，μ0是真空磁导率， 表示电流密度，r是空间场中任意一点到单位体积元dv的距离。Βz(x,y,z)是测量到竖直方向的心磁信号；Bx与By是x,y方向的心磁信号；

[0033] (1)在闭合区域D，边M1M2所在方向，有方程组： 和 得到唯一的驻点M4(1.2,9.3)，故该驻点值为极大值点；将该极大值点
坐 标代入式 从而求的所对应的心磁信号值T4(74pT)；

[0034] (2)采用上述同样的原理和方法，得到边M2M3所在方向的驻点M5(-1.4,11)和所对应的心磁信号值T5(77.5pT)；

[0035] (3)采用(1)中的原理和方法，得到边M3M1所在方向的驻点M6(-0.9,8.6)和所对应的心磁信号值T6(81pT)；

[0036] 8)对于点M4M5M6构成的三角形，根据心磁信号值T4、T5和T6和相应的驻点坐标，得到边M4M5，M5M6以及边M6M4所在方向心磁信号值T7(80.87pT)和驻点坐标(-0.085,9.85)，心磁信号值T8(81.02pT)和驻点坐标(-0.087,9.82)，心磁信号值T9(80.93pT和驻点坐标(-0.083,9.80)；

[0037] 9)上述心磁信号值T7、T8和T9，如图4所示，三者的横坐标误差在±0.003cm以内，纵坐标误差在±0.03cm以内，幅值误差在±0.05pT以内，在测量过程中，均属于允许误差，故取心磁信号值T8(81.02pT)和驻点坐标M8(-0.087,9.82)为心磁测量区域极大值点，如图5所示，a)是心磁信号T值的三维磁场图，b)是相应的二维磁场图。

[0038] 上述实施例提供了一个胸前心磁测量区域的极大值点的测量方法，本发明中，心磁测量区域可以在人体心脏周围任意选定，可以是胸前，也可以是背后，只是建立坐标系时，测量点位置选择不同。如本实施例中的胸前测量范围的框定；框定的这个正方形是与人体心磁信号极值区域相对应，通过实施本方法，可以在极值区域中寻找到极大值点。

[0039] 在本发明中，除采用平面直角坐标系，还可以采用其他坐标系。在本实施例中直角坐标系的y轴方向是人体躯干方向，除此之外，还可以根据测量需要任意选定。

[0040] 本发明不局限于此实施例，还可以根据实际需要和设计要求做出相应的修改，例如：

[0041] 测量正方形的边长可约10～15厘米；OP的距离可约为1-5cm。

[0042] 以上通过详细实施例描述了本发明所提供的心磁测量区域极值点的测定方法，本领域的技术人员应当理解，在不脱离本发明实质的范围内，可以对本发明做一定的变形或修改；其实施方法也不限于实施例中所公开的内容。