CN105212897A - 一种基于心磁特征确定心脏磁场极大值点的方法 - Google Patents
一种基于心磁特征确定心脏磁场极大值点的方法 Download PDFInfo
- Publication number
- CN105212897A CN105212897A CN201510526481.4A CN201510526481A CN105212897A CN 105212897 A CN105212897 A CN 105212897A CN 201510526481 A CN201510526481 A CN 201510526481A CN 105212897 A CN105212897 A CN 105212897A
- Authority
- CN
- China
- Prior art keywords
- point
- mcg
- magnetic field
- maximum point
- signals
- 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.)
- Granted
Links
Landscapes
- Measuring Magnetic Variables (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于心磁特征快速确定心脏磁场极大值点方法,该方法首先在人体胸前或背后心磁测量区域建立坐标系,在上述坐标系中,画一正方形框定人体测量范围,选择相邻两条边的中点及对应顶点构成三角形,并以此三角形的三个顶点为测量点,分别进行心磁信号测量,然后对三边所在方向迭代两次求得对应的驻点和相应极值,对比上述所求极值即可得到极大值位点,上述水平分量的心脏磁场是单一极大值模式。利用本发明可在选定的心磁测量区域中准确测量出信号极值点的位置并给出极值点的心磁信号幅度,便于获得有价值的心磁信息。
Description
技术领域
本发明是关于心脏磁场信号测量技术,具体涉及一种基于心磁特征快速确定心脏磁场极大值点的方法。
背景技术
心磁信号比常规的心电信号可以包含更丰富的信息,通过人体心脏磁场的测量进行心脏疾病的诊断,是近年来发展起来的全新的心脏病诊断技术,比常规心电图对某些心脏病的诊测更灵敏,具有无接触、无创的优势,有望在不久的将来用作临床诊断心脏疾病的重要技术手段。
在心脏磁场测量过程中,为了能够得到足够的心脏磁场分布信息,进行基于心磁的心脏物理相关研究及进行心脏疾病的临床诊断,需要在人体体表(例如胸前或背后)进行多点(如36点)的心磁测量,如图1所示。利用心磁图仪测得每个点的心磁信号后,在去除了心磁信号的直流项后,将R波峰值大小定义为心磁信号大小,如图2所示的心磁信号大小约为62pT。
目前,多点测量点位方法通常采用基于人体结构确定心磁测量覆盖范围,其覆盖范围的选择,是根据人体骨骼、胸围等解剖结构来进行,具体步骤包括:首先在人体体表选定一块与心脏对应的区域作为测量覆盖范围,并将其划分成网格,进而测量格点位置上的心磁信号随时间的变化。该方法选定的覆盖区域的大小、位置,以及网格的划分,对最后所获得的心磁图的有效信息含量及相关分析有很大的影响,由于人体心脏在胸腔中的位置及心动强弱的个体差异较大。如采用保守的大面积覆盖方式测量,极易导致边界点附近的采集信号太小,结果没有实际意义,如测量范围过小,难以覆盖心磁信号的有效区域,致使许多对诊断有潜在价值的心磁信息被遗漏。因此,以人体解剖结构来确定多点测量点位的方法对于实际工作,特别是对于临床应用非常不利,常常存在不合理之处。
发明内容
本发明提出一种基于心磁特征快速确定心脏磁场极值点方法,该方法可在心磁测量区域中准确测量出信号极值点的位置并给出极值点的心磁信号幅度,便于获得有价值的心磁信息。本发明提供的心磁测量区域极大值点的测定方法,其步骤包括:
1)在人体胸前或背后的心磁测量区域上建立一坐标系;
2)在上述坐标系中,画一正方形框定人体测量范围,选择其相邻两条边的中点以及对应顶点构成三角形;
3)选择上述三角形的三个顶点为测量点,分别得到心磁信号值T1、心磁信号值T2和心磁信号值T3;
4)分别对上述三角形的三边所在方向的目标函数采用正交偏导数求极值法得到驻点M4、M5和M6;
5)选择M4M5M6构成的三角形,重复采取步骤4),得到驻点M7、M8和M9,将驻点M7、M8和M9的坐标代入公式从而得到对应的心磁信号值T7、T8和T9;
6)将上述心磁信号值T7、T8和T9中的最大值作为该心磁测量区域极大值点的心磁信号值,该极大值点的坐标点为心磁测量区域的极大值点位置。
所述坐标系为平面直角坐标系,坐标系以左侧锁骨中线的人体躯干方向为y轴,x轴方向为与y轴垂直方向的经过人体剑突的直线,相交点O为原点。
所述正方形的中心位于坐标系y轴+6厘米~+13厘米处,边长为10~15厘米。
所述正交偏导数求值法采用数值分量的偏微分法。
本发明的技术效果是:
心磁信号具有x,y,z三个方向,根据测量的z方向的心磁信号值采用正交偏导数方法得到xy方向的心磁信号值,其心脏磁场在xy平面呈现单一极大值模式,故通过本发明所提供的方法,可以快速、方便地在选定的测量区域中获得心磁信号的极大值点。以此极大值点的位置和信号幅度为依据,在进行心磁的多点测量之前,从信号源的本征特性出发,进而合理确定进行多点心磁测量所需覆盖的有效区域,既不会遗漏最强点信号,同时又不必包含没有诊断价值的信号,完全避免根据人体结构确定测量覆盖范围所带来的各种弊端,确保测量的结果有足够的心磁研究和医学诊断价值,为后续相关研究及诊断信息的提取、促进心磁图技术的临床应用提供有利条件。
附图说明
图1为现有的心磁图测量区域划分方法示意图;
图2为现有的心磁信号波形示意图;
图3为本发明测定胸前心磁极大值点过程中各采样点示意图;
图4为本发明实施例中得到的心磁信号值T和驻点M;
图5为本发明实施例得到的极大值磁信号值T的磁场图,(a)是三维磁场图,(b)是二维等磁图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细描述:
以确定胸前一心磁测量区域的极大值点的位置与信号幅度为例:
1)在胸前心磁测量区域内建立平面直角坐标系;
2)如图3所示:本发明实施例采用的平面直角坐标系中,y轴方向大致为人体躯干方向,从下指向上,x轴方向为与y轴垂直的方向,从右指向左;
3)在上述坐标系中,选择以过人体剑突的直线为x轴,以左侧锁骨中线为y轴,相交点O为原点;
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);
5)以上述正方形ABCM1的边BC的中点M2(0,15),边AB的中点M3(-6,9),正方形顶点M1(6,3)为顶点构成三角形,如图3所示;
6)以上述三角形M1M2M3的三个顶点为测量点进行心磁信号测量,分别得到心磁信号T1(12pT)、心磁信号值T2(40pT)和心磁信号值T3(25pT);
7)根据心磁信号值T1、T2和T3,采用正交偏导数求极值法得到边M1M2,M2M3以及边M3M1所在方向心磁信号值T4,心磁信号值T5和心磁信号值T6,如图4所示。
在数学中,一个多变量的函数的偏导数,就是它关于其中一个变量的导数而保持其他变量恒定。偏导数在工程中是很有用的。在工程中遇到的函数大部分是连续的,或者能保证在所讨论的区域内取得它的极大值或极小值。如果知道可微函数的极大值或极小值一定在区域内,且函数在区域内有仅有一个驻点,则该驻点一定是极大值或极小值。
在本实施例中,心磁信号是连续的,且是矢量,分量BxByΒz相互垂直,因而是正交的。正交偏导数的定义是:对于闭合区域D={(x,y)|-6cm≤x≤6cm,3≤y≤15cm},目标函数式中,μ0是真空磁导率,表示电流密度,r是空间场中任意一点到单位体积元dv的距离。Βz(x,y,z)是测量到竖直方向的心磁信号;Bx与By是x,y方向的心磁信号;
(1)在闭合区域D,边M1M2所在方向,有方程组:和得到唯一的驻点M4(1.2,9.3),故该驻点值为极大值点;将该极大值点坐标代入式从而求的所对应的心磁信号值T4(74pT);
(2)采用上述同样的原理和方法,得到边M2M3所在方向的驻点M5(-1.4,11)和所对应的心磁信号值T5(77.5pT);
(3)采用(1)中的原理和方法,得到边M3M1所在方向的驻点M6(-0.9,8.6)和所对应的心磁信号值T6(81pT);
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);
9)上述心磁信号值T7、T8和T9,如图4所示,三者的横坐标误差在±0.003cm以内,纵坐标误差在±0.03cm以内,幅值误差在±0.05pT以内,在测量过程中,均属于允许误差,故取心磁信号值T8(81.02pT)和驻点坐标M8(-0.087,9.82)为心磁测量区域极大值点,如图5所示,a)是心磁信号T值的三维磁场图,b)是相应的二维磁场图。
上述实施例提供了一个胸前心磁测量区域的极大值点的测量方法,本发明中,心磁测量区域可以在人体心脏周围任意选定,可以是胸前,也可以是背后,只是建立坐标系时,测量点位置选择不同。如本实施例中的胸前测量范围的框定;框定的这个正方形是与人体心磁信号极值区域相对应,通过实施本方法,可以在极值区域中寻找到极大值点。
在本发明中,除采用平面直角坐标系,还可以采用其他坐标系。在本实施例中直角坐标系的y轴方向是人体躯干方向,除此之外,还可以根据测量需要任意选定。
本发明不局限于此实施例,还可以根据实际需要和设计要求做出相应的修改,例如:
测量正方形的边长可约10~15厘米;OP的距离可约为1-5cm。
以上通过详细实施例描述了本发明所提供的心磁测量区域极值点的测定方法,本领域的技术人员应当理解,在不脱离本发明实质的范围内,可以对本发明做一定的变形或修改;其实施方法也不限于实施例中所公开的内容。
Claims (4)
1.一种确定心脏磁场极大值点的方法,其步骤包括:
1)在人体胸前或背后的心磁测量区域上建立一坐标系;
2)在上述坐标系中,画一正方形框定人体测量范围,选择其相邻两条边的中点以及对应顶点构成三角形;
3)选择上述三角形的三个顶点为测量点,分别得到心磁信号值T1、心磁信号值T2和心磁信号值T3;
4)分别对上述三角形的三边所在方向的目标函数采用正交偏导数求极值法得到驻点M4、M5和M6;
5)选择M4M5M6构成的三角形,重复采取步骤4),得到驻点M7、M8和M9,将驻点M7、M8和M9的坐标代入公式从而得到对应的心磁信号值T7、T8和T9;
6)将上述心磁信号值T7、T8和T9中的最大值作为该心磁测量区域极大值点的心磁信号值,该极大值点的坐标点为心磁测量区域的极大值点位置。
2.如权利要求1所述的确定心脏磁场极大值点的方法,其特征在于,所述坐标系为平面直角坐标系,坐标系以左侧锁骨中线的人体躯干方向为y轴,x轴方向为与y轴垂直方向的经过人体剑突的直线,相交点O为原点。
3.如权利要求2所述的确定心脏磁场极大值点的方法,其特征在于,所述正方形的边长为10~15厘米。
4.如权利要求4所述的确定心脏磁场极大值点的方法,其特征在于,所述正方形的中心位于坐标系y轴+6厘米~+13厘米处。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510526481.4A CN105212897B (zh) | 2015-08-25 | 2015-08-25 | 一种基于心磁特征确定心脏磁场极大值点的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510526481.4A CN105212897B (zh) | 2015-08-25 | 2015-08-25 | 一种基于心磁特征确定心脏磁场极大值点的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105212897A true CN105212897A (zh) | 2016-01-06 |
CN105212897B CN105212897B (zh) | 2017-10-03 |
Family
ID=54982450
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510526481.4A Expired - Fee Related CN105212897B (zh) | 2015-08-25 | 2015-08-25 | 一种基于心磁特征确定心脏磁场极大值点的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105212897B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5206589A (en) * | 1990-09-30 | 1993-04-27 | Agency Of Industrial Science & Technology | Method of measuring vector magnetic flux by performing curved surface interpolation with data obtained by plural squid type vector magnetic flux meters, and apparatus for carrying out the same |
JP2004358278A (ja) * | 2004-09-28 | 2004-12-24 | Hitachi Ltd | 磁場計測装置の表示方法 |
EP1642526A1 (en) * | 2004-09-29 | 2006-04-05 | Hitachi High-Technologies Corporation | Biomagnetic measurement apparatus |
US20080161190A1 (en) * | 2006-11-14 | 2008-07-03 | Korea Research Institute Of Standards And Science | Biomagnetic field measurement apparatus |
CN101361655A (zh) * | 2008-09-12 | 2009-02-11 | 北京大学 | 一种心磁测量区域极值点的测定方法 |
CN101862183A (zh) * | 2010-06-13 | 2010-10-20 | 北京大学 | 一种基于心磁特征快速确定心磁测量范围的方法 |
CN103142211A (zh) * | 2011-12-07 | 2013-06-12 | 同济大学 | 一种基于极值圆的心脏磁信号处理方法 |
-
2015
- 2015-08-25 CN CN201510526481.4A patent/CN105212897B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5206589A (en) * | 1990-09-30 | 1993-04-27 | Agency Of Industrial Science & Technology | Method of measuring vector magnetic flux by performing curved surface interpolation with data obtained by plural squid type vector magnetic flux meters, and apparatus for carrying out the same |
JP2004358278A (ja) * | 2004-09-28 | 2004-12-24 | Hitachi Ltd | 磁場計測装置の表示方法 |
EP1642526A1 (en) * | 2004-09-29 | 2006-04-05 | Hitachi High-Technologies Corporation | Biomagnetic measurement apparatus |
US20080161190A1 (en) * | 2006-11-14 | 2008-07-03 | Korea Research Institute Of Standards And Science | Biomagnetic field measurement apparatus |
CN101361655A (zh) * | 2008-09-12 | 2009-02-11 | 北京大学 | 一种心磁测量区域极值点的测定方法 |
CN101862183A (zh) * | 2010-06-13 | 2010-10-20 | 北京大学 | 一种基于心磁特征快速确定心磁测量范围的方法 |
CN103142211A (zh) * | 2011-12-07 | 2013-06-12 | 同济大学 | 一种基于极值圆的心脏磁信号处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105212897B (zh) | 2017-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110007350B (zh) | 一种磁探测方法盲区的分析方法 | |
US9442088B2 (en) | Single coil magnetic induction tomographic imaging | |
CN102402796B (zh) | 肝脏血管系统的三维结构化描述方法 | |
KR101912233B1 (ko) | 물체의 위치를 결정하기 위한 시스템 및 방법 | |
RU2471237C2 (ru) | Создание модели приведенного порядка электромагнитного ответного сигнала от подземной структуры | |
AU2016313274A1 (en) | Hand held devices for magnetic induction tomography | |
CN103430213A (zh) | 追踪肿瘤的方法 | |
CA2940245C (en) | Coil for magnetic induction tomography imaging | |
CN110207688B (zh) | 一种基于特征矢量的磁信标快速定向方法与系统 | |
CN103969627A (zh) | 基于fdtd的探地雷达大规模三维正演模拟方法 | |
CN113051779B (zh) | 一种三维直流电阻率法数值模拟方法 | |
CN105118091A (zh) | 一种构建多精度非均匀地质网格曲面模型的方法和系统 | |
CN100581453C (zh) | 一种心磁测量区域极值点的测定方法 | |
Somphone et al. | Fast myocardial motion and strain estimation in 3D cardiac ultrasound with sparse demons | |
CN110837079B (zh) | 一种基于雷达的目标检测方法及装置 | |
CN103654775A (zh) | 可视化crbbb/clbbb心脏电活动的等效源重构方法 | |
CN101862183B (zh) | 一种基于心磁特征快速确定心磁测量范围的方法 | |
CN105212897A (zh) | 一种基于心磁特征确定心脏磁场极大值点的方法 | |
CN105193414A (zh) | 一种快速确定心脏磁场极值的方法 | |
CN102551711B (zh) | 用超导量子干涉器测量磁场并定位磁场电流源的方法 | |
CN101836860B (zh) | 一种心脏磁场最大值和最小值的点位确定方法 | |
CN114252929A (zh) | 一种矿井底板水源核磁共振探测装置及其工作方法 | |
CN102435866A (zh) | 探地雷达考古探测地面物体干扰快速识别方法 | |
Xue et al. | A new approach for the bundle adjustment problem with fixed constraints in stereo vision | |
CN110720914A (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: 20171003 Termination date: 20180825 |