专利名称：基于心脏单电流偶极子模型确定心磁测量范围的方法
技术领域：
本发明是关于心脏磁场信号测量技术，具体涉及一种基于心脏单电流偶极子模型 确定心脏磁场测量范围的方法。
背景技术：
人体心磁信号是由心脏电生理活动激发的生物电流所产生的，是心磁图仪在体 表测量的结果。相较于常规的心电信号，它包含更丰富的信息，且对某些心脏病的诊断更 灵敏，是一种非接触性、无侵袭性的检测方法，可以用于心脏功能疾病的早期诊断和辅助治 疗，因此，有望在不久的将来用作临床诊断心脏疾病的重要技术手段。目前，心磁信号的多点测量点位及测量覆盖范围的确定通常是采用基于人体结构 确定心磁测量覆盖范围的方法，其覆盖范围的选择，是根据人体骨骼、胸围等解剖结构来进 行，但是这种依据人体生理结构确定心磁测量范围的方法通常存在不合理之处，容易导致 所选测量范围过大或过小，没有将心磁偶极子信号区域合理覆盖。测量范围过大时，极易导 致边界点附近信号太小，信号采集结果没有实际意义等一系列问题，使得获取的部分心磁 数据没有参考和研究价值；测量范围过小时，难以覆盖心磁信号的有效区域，致使许多对诊 断有潜在价值的心磁信息被遗漏。为克服这些问题，往往只有在多次多点测量结果之上，修 正一个合理的测量覆盖范围后，才能获得准确、合理、有使用价值的测量结果。而这对于实 际工作、特别是对于临床应用非常不利。在实际操作过程中，仪器的实用性和高效性至关重 要，但是多次测量心磁和修改测量范围，过于麻烦，不能满足医生和患者的实际需要。
发明内容
本发明提出一种基于心脏单电流偶极子模型快速确定心脏磁场测量范围的方法。本发明的技术方案是一种基于心脏单电流偶极子模型快速确定心脏磁场测量范围的方法，其步骤包 括第一步，在人体胸前或背后建立一个三维直角坐标系，该坐标系的坐标原点为人 体剑突对应的位置，y轴方向为人体躯干纵向平行方向，X轴方向与1轴方向垂直，XOy面近 似与人体胸前或背后的心磁测量平面相重合，Z轴方向与X、y轴满足右手螺旋定则，从人体 背后指向胸前；测得的心磁信号为垂直于测量平面的Z分量，设定沿Z轴正向的心磁信号为 心磁正值，反之，为心磁负值。第二步，通过计算得到心脏磁场的最大值Tmax、最小值Tmin，以及对应的坐标点 max、Xmax，^max ‘ ^ Λ min、Xmin， ymin，0)。其中，定义心磁信号正值区域(正向区域)的绝对极大
值为心磁最大值，心磁信号负值区域(负向区域)的绝对极大值为心磁最小值。第三步，确定心脏单电流偶极子的位置。其中，单电流偶极子的χ、y方向坐标分 别取心磁最大值点和心磁最小值点X、y方向坐标的算术平均值，这种单电流偶极子位置 的选取能够协调电流偶极子产生的磁场与实际测量心脏磁场的范围，有利于覆盖绝大部分的心磁极大值和极小值区域；单电流偶极子的深度d(z方向坐标的负值)是心磁最大 值和最小值点间距D(D= |Mmax-Mmin|)的1/7^，即 / = /)/^。所以，单电流偶极子的坐标
((JC +Λ . )/2,(v +J . )/2，-乃/Λ/^)。
ν maxmm ^maxmm ^下面解释如何得到心脏单电流偶极子深度d与心磁最大值和最小值点间距D的关 系:d = D /芯。若只测量心磁的ζ分量，单电流偶极子在测量平面上激发的磁场ζ分量表达式
^ u q (y-yJ-Q ^JB =^--y^ -----(1)
ζ+(J-J0)2+(ζ-ζ0)2Γ其中，qx、qy分别是单电流偶极子在x、y方向的幅值，(x，y，ζ)是测量点的坐标， 按照所建坐标系，ζ = 0，而(X(1，y0, Z0)是单电流偶极子的坐标，Z0 = -d, μ ο是真空中的磁 导率。这样，Bz可以表示成二元函数的形式/(λ:, ν) = B (x,y) = —---- (2)
4π [(x-xj+(y-yJ+dT2此二 元函数满足的边界和极值条件是
Iim f(x, j) = 0， Iim f(x, j) = 0，由此可知在距离单电流偶极子为d的平面上f (X，y)存在最值，通过多元函数的 极值解法fx(x，y) = 0，fy(x, y) = 0，可以解出最大值和最小值的坐标分别为
,qydqdqd
(X ——,，Vn + I=，0)， (χ + ,=，
λ/Κ) Λ/Κ)Λ/Κ)
q d
凡-ι 巧，0)，所以，第四步，确定心脏单电流偶极子的幅值。利用第二步获得的心脏磁场最大值Tmax和 最小值Tmin，以及对应的坐标点Mmax (Xmax，ymax, 0)、Mmin(xmin，ymin，0)，和第三步确定的单电流偶 极子坐标，求解方程组(3)、(4)
u q (v -Vn)-<7^
___Λ^Ο1 χ maxO / Jvv- maxO /rT--= 1
a Q (y -yJ-Q^
___Λ^Ο1 χ mmO / -i ν V mmO /rT--= 1
^[(xmm-xJ+(ymm-yJ+dT^ ■由上两式可以解出单电流偶极子的幅值邓和qy。第五步，以心脏单电流偶极子为心脏电流源在心磁测量平面上模拟激发磁场，如 图5所示。计算的依据是电流偶极子在测量平面上激发的磁场ζ分量表达式(1)。其中，单
电流偶极子的坐标选为((JC +Λ. )/2, (ν +V . )/2, -乃/λ/^)，单电流偶
ν maxmm ^maxmm ^
(3)
(4)极子的幅值由第四步得到。第六步，由单电流偶极子激发的心磁地图，如图5所示，计算心磁幅值降为某一较 小幅度时的，或降为心磁极值某一百分比时的等磁线与χ、y方向直线的切点。其中，极大 值区域的等磁线与X方向直线I相切于点A(xa，ya，0)，与y方向直线II相切于点B(xb，yb， 0)；极小值区域的等磁线与χ方向直线III相切于点C(x。，y。，0)，与y方向直线IV相切于 点 D (xd, yd,0)。考虑到实际心磁测量范围的需要和心磁测量仪器目前的探测能力，若人体心磁信 号较弱，直线I，II，III，IV应选择与小幅度的心磁等磁线(例如，绝对值小于等于IOpT的 等磁线)相切；若人体心磁信号较强，直线I，II，III，IV应选择与心磁极值某一百分比的 等磁线(例如，小于等于1/3心磁极值的等磁线)相切。第七步，取上述四条切线(直线I，II，III，IV)构成的矩形为心磁测量范围，如图 5所示。该矩形边长沿χ轴为k-xd|，沿y轴为k-yj。第二步中，利用插值计算得到心脏磁场的最大值点Mmax和最小值点Mmin，具体包 括1)在上述坐标系中，分别以剑突和锁骨肩峰端为圆心，半径为小于5cm画圆，在上 述两圆内任意选择两个测量点，分别为测量点M1和测量点M5 ；2)将测量点M1和测量点M5连成一直线，在该直线上位于M1点与M5点之间确定三 个心磁信号测量点，分别是测量点M2、测量点M3和测量点M4,上述五个测量点间距相等；3)对测量点M1、测量点M2、测量点M3、测量点M4和测量点M5进行心磁信号测量，分 别得到心磁信号值T1、心磁信号值T2、心磁信号值T3、心磁信号值T4和心磁信号值T5,对上 述5个心磁信号值做插值，同时得到两个极值点极大值点Tltl，对应的坐标点Mltl，极小值点 T11,对应的坐标点M11 ；4)在经过坐标点Mltl，垂直于上述直线的一条直线上，另选择两个测量点，分别是 测量点M6和测量点M7,对测量点M6和测量点M7进行心磁信号测量，分别得到心磁信号值T6 和心磁信号值T7,对上述心磁信号值Tltl、心磁信号值T6和心磁信号值T7做插值得到一极值 Tmax，对应坐标点Mmax，Mfflax就是心磁信号的最大值的点位；5)在经过坐标点M11,垂直于上述直线的一条直线上，另选择两个测量点，分别是 测量点M8和测量点M9,对测量点M8和测量点M9进行心磁信号测量，分别得到心磁信号值T8 和心磁信号值T9,对上述心磁信号值T11、心磁信号值T8和心磁信号值T9做插值得到一极值 Tmin，对应坐标点Mmin，Mfflin就是心磁信号的最小值的点位。所述测量点M6和测量点M7可分别位于坐标点Mltl的两侧，测量点M6和坐标点Mltl 之间距离与测量点M7和坐标点Mltl之间的距离相等，分别相距2 6厘米。所述测量点M8和测量点M9可分别位于坐标点M11的两侧，测量点M8和坐标点M11 之间距离与测量点M9和坐标点M11之间的距离相等，分别相距2 6厘米。所述插值选用三次样条插值函数。本专利凭借由心磁最大值和最小值及其点位确定的心脏单电流偶极子来模拟实 测心脏磁场，然后依据单电流偶极子场的衰减来确定心磁测量范围，这种方法能够方便地 确定合理的心磁测量所需覆盖的范围。单电流偶极子模型把心脏内的生物电流综合看成位 于心脏内某点的一小段电流，是一种简易的心脏功能模型，但能够较好地表征心动周期中，特别是QRS和ST-T波段的心脏电生理活动，并且它所激发的模拟心磁与实测心磁有很大的 相似性，所以，可以通过确定单电流偶极子的磁场范围来选择实际的心磁测量范围。与现有技术相比，本发明的有益效果是参考图2，人体的心磁信号分为正向区域和反向区域，这两个区域通常位于剑突水 平线和胸骨上缘水平线之间。本发明利用插值计算可得到心脏磁场的最大值Tmax和最小值 Tmin，以及对应的坐标点Mmax和Mmin，利用最大值点和最小值点来确定心脏单电流偶极子的位 置，通过模拟电流偶极子在心磁测量平面上激发的磁场来计算心磁幅值的衰减程度，进而 得到覆盖了心磁正向和反向区域的矩形区域，此即为进行多点心磁测量的合理覆盖范围， 在此矩形范围之内，可以很方便地进行测量格点的划分并作多点心磁测量。本发明在充分考虑到心脏(即信号源)的个体差异情况下，利用信号源的本征特 性，确定心磁测量范围，不但能通过简单的几步得到，避免按照人体结构确定测量范围时不 易有效覆盖心磁正向和反向区域的缺点，而且可以根据具体测量需要，灵活地确定多点心 磁测量覆盖范围的大小，有利于尽可能多地获取有用的心磁信息，并且在此范围上进行测 量后，绘制心磁地图的等磁线与测量区域边框近似相切，为后续相关研究及诊断信息的提 取提供有利条件。
图1为现有的基于人体结构确定心磁测量覆盖范围方法示意图，其中(a)为胸前 心磁测量区域，(b)为背后心磁测量区域；图2为胸前心磁各测量格点的心磁信号构成的心磁地图，“ + ”表示心磁信号的正 向区域，“_”表示心磁信号的反向区域；图3为实际测量心磁地图与人体胸前位置对应图；图4为本发明插值计算得到心磁最大值点Mmax和心磁最小值点Mmin的示意图；图5本发明利用心脏单电流偶极子确定心脏磁场测量范围的示意图。图6本发明实施例中的利用心脏单电流偶极子确定心脏磁场测量范围的示意图。
具体实施例方式下面结合附图和具体实施方式
对本发明作进一步详细描述参考图2，以胸前心磁测量区域的极值点的位置与信号幅度为例1)在人体胸前的心磁测量区域上建立一个三维坐标系，本发明采用三维平面直角 坐标系，该坐标系的坐标原点为人体剑突对应的位置，y轴方向为人体躯干纵向平行方向， 从人体下部指向上部，χ轴方向与y轴方向垂直，从人体右部指向左部，xOy面近似与人体胸 前的心磁测量平面相重合，ζ轴方向与x、y轴满足右手螺旋定则，从人体背后指向胸前。2)上述坐标系中，分别以人体剑突和锁骨左肩峰端为圆心，半径为3cm画圆，在上 述两圆内任意选择两个测量点，分别为测量点M1 (0，0，0) cm和测量点M5 (20，20,0) cm。本实 施例选择M1点为人体剑突对应的位置，M5点为锁骨左肩峰端对应的位置，参考图4。3)将测量点M1和测量点M5连成直线，在此直线上位于M1点与M5点之间再确定三 个心磁信号测量点，分别是测量点M2 (5，5，0) cm、测量点M3 (10，10,0) cm和测量点M4 (15，15， 0) cm,上述五个测量点间距相等。
4)对测量点M1、测量点M2、测量点M3、测量点M4和测量点M5进行心磁信号测量，分 别得到心磁信号值T1 = 10. 9pT、心磁信号值T2 = 30. 7ρΤ、心磁信号值T3 = -14. 8ρΤ、心磁 信号值T4 = -2. 4ρΤ和心磁信号值T5 = -0. 3ρΤ，对上述5个心磁信号值做插值，得到一个 极大值Tltl = 38. ΟρΤ，对应的坐标点Mltl (3. 9，3.9，0) cm，和一个极小值T11 = -17. 2pT，对应 的坐标点Μη(11. 3，11. 3，0)cm，如图3所示。5)在过极大值点Mltl，且垂直于上述直线的方向上选取两个测量点M6(6. 9,0. 9,0) cm、M7 (0. 9,6. 9,0) cm,使其分别位于Mltl点两侧，且与Mltl点等间距，间距为2 6厘米，分 别得到心磁信号值T6 = 18. 3pT、心磁信号值T7 = 10. 5ρΤ。对点Μ6、Μ7、Μ1(1对应的心磁信号 测量值T6、T7、T10做插值，得到一个极大值Tmax = 39. 9ρΤ，对应坐标点Mmax (xmax, ymax，zmax)= (4. 1,3. 7，0)cm。Tmax就是整个心磁测量区域心磁信号的最大值，Tmax对应的坐标点Mmax为心 磁测量区域的最大值的点位，如图4所示。6)同理，在过极小值点M11,垂直于上述直线方向上选取两个测量点M8(14. 3,8. 3， 0)cm、M9 (8. 3，14. 3，0) cm，使其分别位于M11点两侧，且与M11点等间距，间距为2 6厘米， 分别得到心磁信号值T8 = -12. 7pT、心磁信号值T9 = -2. 5ρΤ。对点M8、M9, M11对应的心磁 信号测量值T8、T9、T11做插值，得到一个极小值Tmin = -18. 6ρΤ，对应坐标点Mmin(xmin，ymin， Zfflin) = (12. 3，10. 3，0) cm。Tmin就是整个心磁测量区域心磁信号的最小值，Tmin对应的坐标 点Mmin为心磁测量区域的最小值的点位，如图4所示。7)禾Ij用心磁最大值、最小值、及其点位确定
单电流偶极子的位置和幅值。单电流偶极子的位置
(( Χ +Λ . )/2, (v. )/2，-D/λ/^) = (8.2，7.0，-6.33) cm,单电流
ν maxmm ^maxmm ^
偶极子的幅值(Px，Py) = (-1.225,2. 791) μ Am。8)已知心脏单电流偶极子后，在测量平面上模拟激发磁场，此磁场最大值和最小 值的绝对值均为29. 3ρΤ。磁场的等磁线幅值从极值向四周逐渐衰减，根据实际需要，选择心 磁幅值衰减为IOpT作为确定心磁测量范围的标准。χ方向直线I、y方向直线II与心磁极 大值区域的IOpT心磁等磁线相切(IOpT大致为心磁最大值的1/3)，x方向直线III、y方向 直线IV与心磁极小值区域的-IOpT心磁等磁线相切(-IOpT大致为心磁最小值的1/3)，如 图6所示，这四条直线围成30X 25cm2的矩形心磁测量范围。通过以上步骤就确定了心脏磁场的测量区域。本发明中，心磁测量区域可以是胸前，也可以是背后，在背后坐标系中，M1和礼点 的选取分别以与人体剑突和锁骨左肩峰端对应的点为圆心画圆。若人体是右位心，应根据左右心位的对称关系，M1和M5点的选取分别以人体剑突 和锁骨右肩峰端对应的点为圆心画圆，采用相同的方法，进行心磁测量范围的确定。以上通过详细实施例描述了本发明所提供的心磁测量区域的确定方法，本领域的 技术人员应当理解，在不脱离本发明实质的范围内，可以对本发明做一定的变形或修改；其 实施方法也不限于实施例中所公开的内容。
权利要求
一种基于心脏单电流偶极子模型确定心脏磁场测量范围的方法，其包括如下步骤第一步，在人体胸前或背后建立一个三维直角坐标系，该坐标系的坐标原点为人体剑突对应的位置，y轴方向为人体躯干纵向平行方向，x轴方向与y轴方向垂直，z轴方向与x、y轴满足右手螺旋定则；第二步，在上述坐标系中计算得到心脏磁场的最大值Tmax、最小值Tmin，以及对应的坐标点；第三步，根据心磁最大值点和心磁最小值点坐标，确定心脏电流偶极子的位置和幅值；第四步，以心脏电流偶极子为心脏电流源模拟激发磁场，形成心磁地图；第五步，利用心磁幅值的衰减程度，在上述心磁地图上确定一合理的心磁测量范围。
2.如权利要求1所述的方法，其特征在于，确定xOy面近似与人体胸前或背后的心磁测 量平面相重合，利用插值计算在该xOy面上得到心脏磁场的最大值点Mmax和最小值点Mmin， 具体包括1)在上述坐标系中，分别以剑突和锁骨肩峰端为圆心，半径为小于5cm画圆，在上述两 圆内任意选择两个测量点，分别为测量点M1和测量点M5 ；2)将测量点M1和测量点M5连成一直线，在该直线上位于M1点与M5点之间确定三个心 磁信号测量点，分别是测量点M2、测量点M3和测量点M4,上述五个测量点间距相等；3)对测量点M1、测量点M2、测量点M3、测量点M4和测量点M5进行心磁信号测量，分别得 到心磁信号值T1、心磁信号值T2、心磁信号值T3、心磁信号值T4和心磁信号值T5,对上述5 个心磁信号值做插值，同时得到两个极值点极大值点Tltl，对应的坐标点Mltl，极小值点T11, 对应的坐标点M11 ；4)在经过坐标点Mltl，垂直于上述直线的一条直线上，另选择两个测量点，分别是测量 点M6和测量点M7,对测量点M6和测量点M7进行心磁信号测量，分别得到心磁信号值T6和心 磁信号值T7,对上述心磁信号值Tltl、心磁信号值T6和心磁信号值T7做插值得到一极值Tmax， 对应坐标点Mmax，Mfflax就是心磁信号的最大值的点位；5)在经过坐标点M11,垂直于上述直线的一条直线上，另选择两个测量点，分别是测量 点M8和测量点M9,对测量点M8和测量点M9进行心磁信号测量，分别得到心磁信号值T8和心 磁信号值T9,对上述心磁信号值T11、心磁信号值T8和心磁信号值T9做插值得到一极值Tmin， 对应坐标点Mmin，Mfflin就是心磁信号的最小值的点位。
3.如权利要求2所述的方法，其特征在于，所述测量点M6和测量点M7可分别位于坐标 点Mltl的两侧，测量点M6和坐标点Mltl之间距离与测量点M7和坐标点Mltl之间的距离相等， 分别相距2 6厘米。
4.如权利要求2所述的方法，其特征在于，所述测量点M8和测量点M9可分别位于坐标 点M11的两侧，测量点M8和坐标点M11之间距离与测量点M9和坐标点M11之间的距离相等， 分别相距2 6厘米。
5.如权利要求2所述的方法，其特征在于，所述插值选用三次样条插值函数。
6.如权利要求1所述的方法，其特征在于，步骤三中，所述电流偶极子的x、y方向坐标 分别取心磁最大值点和心磁最小值点χ、y方向坐标的算术平均值，电流偶极子在ζ方向坐 标的负值是心磁最大值和最小值点间距D的1/7^。
7.如权利要求1所述的方法，其特征在于，步骤五中，若人体心磁信号较弱，上述心磁 地图上以绝对值小于等于IOpT的等磁线划定测量范围，若人体心磁信号较强，则选择心磁 极值小于等于1/3心磁极值的等磁线划定测量范围。
全文摘要
本发明公开了一种基于心脏单电流偶极子模型确定心脏磁场测量范围的方法，该方法首先在人体胸前或背后建立一个三维直角坐标系，在坐标系中计算得到心脏磁场的最大值Tmax、最小值Tmin，以及对应的坐标点；根据心磁最大值和心磁最小值及其坐标，确定心脏单电流偶极子的位置和幅值；以心脏单电流偶极子为心脏电流源在心磁测量平面上模拟激发磁场，形成心磁地图，利用心磁幅值的衰减程度，在上述心磁地图上确定一合理的心磁测量范围。本发明可快速确定心磁测量范围，同时根据具体测量需要，灵活地确定多点心磁测量覆盖范围的大小，有利于尽可能多地获取有用的心磁信息。
文档编号A61B5/00GK101897576SQ201010229258
公开日2010年12月1日 申请日期2010年7月18日 优先权日2010年7月18日
发明者刘政豪, 张辰, 戴远东, 杨涛, 王倩, 王守证, 王越, 蔡文艳, 谷红芳, 马平, 高智, 魏玉科 申请人:北京大学