CN113876327B - 一种基于serf原子磁强计的高空间分辨率心磁成像方法 - Google Patents
一种基于serf原子磁强计的高空间分辨率心磁成像方法 Download PDFInfo
- Publication number
- CN113876327B CN113876327B CN202111382922.XA CN202111382922A CN113876327B CN 113876327 B CN113876327 B CN 113876327B CN 202111382922 A CN202111382922 A CN 202111382922A CN 113876327 B CN113876327 B CN 113876327B
- Authority
- CN
- China
- Prior art keywords
- array
- layer
- signals
- magnetocardiogram
- atomic magnetometer
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 40
- 238000003491 array Methods 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 9
- 230000000747 cardiac effect Effects 0.000 claims description 8
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 239000010410 layer Substances 0.000 abstract description 53
- 230000017525 heat dissipation Effects 0.000 abstract description 6
- 230000000694 effects Effects 0.000 abstract description 4
- 239000002356 single layer Substances 0.000 abstract description 4
- 239000013598 vector Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000003745 diagnosis Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 239000002355 dual-layer Substances 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 210000002458 fetal heart Anatomy 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 208000031225 myocardial ischemia Diseases 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Cardiology (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measuring Magnetic Variables (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种基于SERF原子磁强计的高空间分辨率心磁成像方法,包括:排布多层SERF原子磁强计阵列模块,标定多层阵列信号投影至最底层阵列的函数关系模块,计算多层阵列生成心磁图像模块。本发明提出的方法相比于基于单层排布SERF原子磁强计阵列的心磁成像方法,具有的优点是在单位面积上可以采集更多信号、传感阵列的散热效果更好、传感阵列的串扰更小、成像空间分辨率更高。
Description
技术领域
本发明属于生物医学信号分析领域,涉及一种基于SERF(Spin ExchangeRelaxation Free,无自旋交换弛豫)原子磁强计阵列的高空间分辨率心磁成像方法。
背景技术
基于SERF原子磁强计的心磁成像是一种新型心脏检测方法,具有无创、无辐射、非接触、室温下工作、高灵敏度、高时间分辨率的特点。研究表明,基于SERF原子磁强计的心磁成像在心律失常、心肌缺血、胎儿心脏等疾病诊断方面优于心电等其他检测方法,有较高的医学研究和应用价值。
基于SERF原子磁强计的高空间分辨率心磁成像方法是基于SERF原子磁强计的心磁诊断技术的应用之一,其原理是使用多层排布的SERF原子磁强计阵列采集心磁信号,把多层阵列心磁信号都垂直投影到最底层阵列进行成像,从而实现高空间分辨率的心磁成像。相比于SQUID(Superconducting Quantum Interference Device,超导量子干涉装置)磁强计,SERF原子磁强计可以在室温工作,并且传感器自由组合实现多层排布,可以实现更高空间分辨率的心磁成像。
目前心磁成像方法具有一定的局限性,主要是成像空间分辨率低、传感阵列的散热效果差、传感阵列的串扰大、在二维图成像用于疾病诊断中包含局部细节信息少。
因此,如何提供一种基于SERF原子磁强计阵列的高空间分辨率心磁成像方法,以解决现有心磁成像的局限性,已成为本领域从业者亟待解决的技术问题。
发明内容
本发明的技术解决问题:克服现有技术的缺点,提供一种基于SERF原子磁强计的高空间分辨率心磁成像方法,用于解决现有技术中心磁成像方法的成像空间分辨率低、传感阵列的散热效果差、传感阵列的串扰大、在二维图成像用于疾病诊断中局部信息少的问题。
为实现上述目的及其他的相关目的,本发明一种基于SERF原子磁强计的高空间分辨率心磁成像方法,包括以下步骤:
步骤一,排布多层SERF原子磁强计阵列。该步骤设置SERF原子磁强计阵列的层数和相邻SERF原子磁强计之间水平方向排布间隔;根据成像空间分辨率指标要求,计算出需要的SERF原子磁强计阵列的层数和相邻SERF原子磁强计之间水平方向排布间隔;得到满足空间分辨率成像要求的多层SERF原子磁强计阵列,该步骤采集的心磁信号可用于步骤三。
步骤二,标定多层阵列信号投影至最底层阵列的函数关系。基于排布多层SERF原子磁强计阵列模块得到确定层数和相邻SERF原子磁强计之间水平方向排布间隔的多层SERF原子磁强计阵列,对多层阵列中每层阵列分别进行标定,标定方法是对每层阵列垂直移动多个高度测量心磁信号,然后拟合心磁信号与高度之间的关系;实现多层阵列中每层阵列的信号都可以根据标定结果计算出投影到最底层高度时的信号幅值,该步骤获得的结果用于步骤三。
步骤三,计算多层阵列生成心磁图像。根据步骤一采集的心磁信号,以及步骤二获得的标定结果,计算多层阵列信号垂直投影至最底层阵列的信号,得到多层阵列信号都化为最底层阵列的信号,对合成后的最底层阵列信号计算心磁图像,包括二维等磁图和二维伪电流密度图。
所述排布多层SERF原子磁强计阵列模块由多层SERF原子磁强计阵列构成,并且任一层SERF原子磁强计阵列中的每个SERF原子磁强计内部敏感源在水平方向位置相对其他层传感阵列的每个SERF原子磁强计内部敏感源都是交叉排布,没有重合排布。
所述标定多层阵列信号投影至最底层阵列的函数关系模块由多次测量每个阵列在不同高度心磁信号幅值,以及函数拟合幅值与高度之间函数关系构成,多层阵列中的每一层阵列先在现有高度采集一组心磁信号,然后垂直向下调低N次阵列高度至最底层阵列的高度,再以相同间隔垂直向上N次调高阵列高度,对调节的2N个阵列高度分别都采集一组心磁信号,最后对获取的2N+1个高度的心磁信号进行拟合,获得多层阵列中每个SERF原子磁强计信号与高度的函数关系式,N≥1。
本发明与现有技术相比具有以下有益效果:
(1)传统心磁成像是使用单层SERF原子磁强计阵列进行采集信号,其心磁成像无法解决在提高成像空间分辨率的同时,消除散热和串扰问题。换言之,传统心磁成像中,如果为了提高成像空间分辨率而采取单层传感阵列中SERF原子磁强计水平方向之间的间距减小,则会造成SERF原子磁强计之间的热量无法散出,并且SERF原子磁强计之间间距过小会有串扰,这都会导致传感阵列无法正常工作;如果为了解决SERF原子磁强计之间的散热和串扰问题,把SERF原子磁强计之间水平方向间距增大,则会导致成像空间分辨率减低。本发明中通过使用多层排布的SERF原子磁强计阵列,提高SERF原子磁强计投影到最底层后在水平面排布SERF原子磁强计的空间密度,在相同水平面获取更多的心磁信号,从而实现在提高成像空间分辨率的同时,消除散热和串扰问题,可以采集信号通道数目提高、成像空间分辨率提高、传感阵列的散热效果提高、传感阵列的串扰消除。
(2)传统心磁成像是使用单层SERF原子磁强计阵列进行采集信号,其阵列中SERF原子磁强计之间水平方向间距较大,在二维图成像中SERF原子磁强计之间部分使用插值绘图,局部真实细节信息无法获取;本发明中通过使用多层排布的SERF原子磁强计阵列,从而实现高空间分辨率的心磁二维图成像,在二维图成像用于疾病诊断中局部细节信息更加丰富。
附图说明
图1为本发明的高空间分辨率心磁成像方法流程示意图;
图2为本发明的双层排布SERF原子磁强计阵列示意图;
图3为本发明的标定双层阵列垂直投影关系示意图。
具体实施方式
下面结合附图及实施例对本发明进行详细说明。
如图1所示,本发明的高空间分辨率心磁成像方法包括以下步骤:
第一步,排布多层SERF原子磁强计阵列。
图2所示是双层排布SERF原子磁强计阵列示意图,从图中可以看出,上层传感阵列中的SERF原子磁强计垂直投影到下层传感阵列时,投影的SERF原子磁强计并不与下层传感阵列中的SERF原子磁强计位置重合,而是叉开排布的。排布三层以及更多层的SERF原子磁强计阵列类似。
双层SERF原子磁强计阵列的成像空间分辨率,实际就是上层阵列中的SERF原子磁强计内部敏感源投影到下层阵列时与下层阵列的SERF原子磁强计敏感源之间在水平方向上的距离。举例说明成像空间分辨率,在图2中SERF原子磁强计A1内部敏感源投影到下层阵列后,与相邻的下层阵列SERF原子磁强计B1敏感源之间在水平方向上的距离就是成像空间分辨率。反之,可以根据需要的高空间分辨率指标,来设计每层传感阵列中SERF原子磁强计敏感源之间距离。
所述SERF原子磁强计内部的敏感源,是指SERF原子磁强计内部气室中抽运光和检测光交叉部分。
第二步,标定多层阵列信号投影至最底层阵列的函数关系。
所述标定多层阵列信号投影至最底层阵列的函数关系可以参阅图3所示,图3所示是标定双层阵列中上层阵列信号投影关系示意图,双层传感阵列的上层阵列先在现有高度采集一组心磁信号,然后垂直向下调低N次阵列高度至最底层阵列的高度,再以相同间隔垂直向上N次调高阵列高度,对调节的2N个阵列高度分别都采集一组心磁信号,最后对获取的2N+1组高度的心磁信号进行拟合,获得上层阵列中每个SERF原子磁强计信号与高度的函数关系式。标定三层以及更多层的SERF原子磁强计阵列垂直投影关系类似。
所述拟合心磁信号和SERF原子磁强计高度之间函数关系,采用的拟合方法为多项式曲线拟合。
第三步,计算多层阵列生成的心磁图像。
根据步骤一采集的心磁信号,以及步骤二获得的标定结果,计算多层阵列信号垂直投影至最底层阵列的信号,得到多层阵列信号都化为最底层阵列的信号,对合成后的最底层阵列信号计算心磁图像,包括心磁二维等磁图和心磁二维伪电流密度图。
所述心磁二维等磁图,为同一时刻各测量点测得的心磁信号中具有相同幅值的测量点相连接,形成等高线图。
所述心磁二维伪电流密度图,根据麦克斯韦方程组对于电流矢量和磁场矢量关系的描述,可利用在平面上测得的心脏磁场分布数据,计算相应的电流源矢量,将各个位置的电流源矢量以箭头图的形式绘制出来,形成的图像即伪电流密度图。
心磁二维伪电流密度图中的电流矢量计算如公式(1)所示:
综上,本发明一种基于SERF原子磁强计的高空间分辨率心磁成像方法,依次进行排布多层SERF原子磁强计阵列模块、标定多层阵列信号投影至最底层阵列的函数关系模块、计算多层阵列生成心磁图像模块,最终实现了高空间分辨率心磁成像,具有的优点是采集信号通道数目提高、成像空间分辨率提高、传感阵列的散热效果提高、传感阵列的串扰消除、在二维图成像用于疾病诊断中局部细节信息更加丰富。所以,本发明有效克服了现有技术中的各个缺点而具有较高的医学应用价值。
上述实施例仅列示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所述技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (3)
1.一种基于SERF原子磁强计的高空间分辨率心磁成像方法,其特征在于,所述方法分为三个模块,分别是排布多层SERF原子磁强计阵列模块,标定多层阵列信号投影至最底层阵列的函数关系模块,计算多层阵列生成心磁图像模块;
所述排布多层SERF原子磁强计阵列模块,设置SERF原子磁强计阵列的层数和相邻SERF原子磁强计之间水平方向排布间隔;根据成像空间分辨率指标要求,计算出需要的SERF原子磁强计阵列的层数和相邻SERF原子磁强计之间水平方向排布间隔;得到满足空间分辨率成像要求的多层SERF原子磁强计阵列,该模块采集的心磁信号用于计算多层阵列生成心磁图像模块;
所述标定多层阵列信号投影至最底层阵列的函数关系模块,基于排布多层SERF原子磁强计阵列模块得到确定层数和相邻SERF原子磁强计之间水平方向排布间隔的多层SERF原子磁强计阵列,对多层阵列中每层阵列分别进行标定,标定方法是对每层阵列垂直移动多个高度测量心磁信号,然后拟合心磁信号与高度之间的关系;实现多层阵列中每层阵列的信号均根据标定结果计算出投影到最底层高度时的信号幅值,该模块获得的结果用于计算多层阵列生成心磁图像模块;
所述计算多层阵列生成心磁图像模块,根据排布多层SERF原子磁强计阵列模块采集的心磁信号,以及标定多层阵列信号投影至最底层阵列的函数关系模块获得的标定结果,计算多层阵列信号垂直投影至最底层阵列的信号,得到多层阵列信号都化为最底层阵列的信号,对合成后的最底层阵列信号计算心磁图像,心磁图像包括心磁二维等磁图和心磁二维伪电流密度图。
2.根据权利要求1所述的一种基于SERF原子磁强计的高空间分辨率心磁成像方法,其特征在于:所述排布多层SERF原子磁强计阵列模块由多层SERF原子磁强计阵列构成,并且任一层SERF原子磁强计阵列中的每个SERF原子磁强计内部敏感源在水平方向位置相对其他层传感阵列的每个SERF原子磁强计内部敏感源都是交叉排布,没有重合排布。
3.根据权利要求1所述的一种基于SERF原子磁强计的高空间分辨率心磁成像方法,其特征在于:所述标定多层阵列信号投影至最底层阵列的函数关系模块由多次测量每个阵列在不同高度心磁信号幅值,以及函数拟合幅值与高度之间函数关系构成,多层阵列中的每一层阵列先在现有高度采集一组心磁信号,然后垂直向下调低N次阵列高度至最底层阵列的高度,再以相同间隔垂直向上N次调高阵列高度,对调节的2N个阵列高度分别都采集一组心磁信号,最后对获取的2N+1个高度的心磁信号进行拟合,获得多层阵列中每个SERF原子磁强计信号与高度的函数关系式,N≥1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111382922.XA CN113876327B (zh) | 2021-11-22 | 2021-11-22 | 一种基于serf原子磁强计的高空间分辨率心磁成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111382922.XA CN113876327B (zh) | 2021-11-22 | 2021-11-22 | 一种基于serf原子磁强计的高空间分辨率心磁成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113876327A CN113876327A (zh) | 2022-01-04 |
CN113876327B true CN113876327B (zh) | 2023-05-26 |
Family
ID=79015908
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111382922.XA Active CN113876327B (zh) | 2021-11-22 | 2021-11-22 | 一种基于serf原子磁强计的高空间分辨率心磁成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113876327B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105640552A (zh) * | 2014-12-02 | 2016-06-08 | 精工爱普生株式会社 | 磁场测量方法以及磁场测量装置 |
CN109998519A (zh) * | 2019-04-19 | 2019-07-12 | 北京航空航天大学 | 一种基于serf原子磁强计的心磁测量及心磁图生成系统 |
CN111281370A (zh) * | 2020-02-19 | 2020-06-16 | 北京航空航天大学 | 一种基于serf原子磁强计的梯度仪配置式脑磁测量系统 |
CN112842344A (zh) * | 2021-02-05 | 2021-05-28 | 浙江工业大学 | 一种磁场检测系统及方法 |
CN112842343A (zh) * | 2021-02-01 | 2021-05-28 | 北京航空航天大学 | 一种基于serf原子磁强计的心磁三维源位置估计方法 |
CN113520399A (zh) * | 2021-07-12 | 2021-10-22 | 北京航空航天大学 | 一种基于亥姆霍兹线圈阵列的心脏磁场模拟系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8688192B2 (en) * | 2011-01-31 | 2014-04-01 | Seiko Epson Corporation | High-resolution magnetocardiogram restoration for cardiac electric current localization |
-
2021
- 2021-11-22 CN CN202111382922.XA patent/CN113876327B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105640552A (zh) * | 2014-12-02 | 2016-06-08 | 精工爱普生株式会社 | 磁场测量方法以及磁场测量装置 |
CN109998519A (zh) * | 2019-04-19 | 2019-07-12 | 北京航空航天大学 | 一种基于serf原子磁强计的心磁测量及心磁图生成系统 |
CN111281370A (zh) * | 2020-02-19 | 2020-06-16 | 北京航空航天大学 | 一种基于serf原子磁强计的梯度仪配置式脑磁测量系统 |
CN112842343A (zh) * | 2021-02-01 | 2021-05-28 | 北京航空航天大学 | 一种基于serf原子磁强计的心磁三维源位置估计方法 |
CN112842344A (zh) * | 2021-02-05 | 2021-05-28 | 浙江工业大学 | 一种磁场检测系统及方法 |
CN113520399A (zh) * | 2021-07-12 | 2021-10-22 | 北京航空航天大学 | 一种基于亥姆霍兹线圈阵列的心脏磁场模拟系统 |
Non-Patent Citations (1)
Title |
---|
基于心磁信号的心脏电流偶极子阵列成像及相关性质的研究;王倩;马平;华宁;陆宏;唐雪正;唐发宽;;物理学报(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113876327A (zh) | 2022-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5594849A (en) | Biomedical magnetism imaging apparatus and method | |
Iivanainen et al. | Measuring MEG closer to the brain: Performance of on-scalp sensor arrays | |
Oyama et al. | Dry phantom for magnetoencephalography—Configuration, calibration, and contribution | |
CN103987337B (zh) | 用于em跟踪的补偿、检测和误差校正的畸变指纹分析 | |
CN112842344B (zh) | 一种磁场检测系统及方法 | |
CN111281370B (zh) | 一种基于serf原子磁强计的梯度仪配置式脑磁测量系统 | |
US20110313274A1 (en) | Methods and apparatuses for 3D imaging in magnetoencephalography and magnetocardiography | |
JP5712640B2 (ja) | 磁気計測装置および生体磁気計測方法 | |
CN112842343B (zh) | 一种基于serf原子磁强计的心磁三维源位置估计方法 | |
CN111289926B (zh) | 一种脑磁测量头盔的设计和标定方法 | |
US20070016000A1 (en) | Three dimensional magnetic resonance motion estimation on a single image plane | |
JPH0542119A (ja) | 生体磁気計測装置 | |
Adachi et al. | Calibration of room temperature magnetic sensor array for biomagnetic measurement | |
Gu et al. | Automatic coregistration of MRI and on-scalp MEG | |
CN113876327B (zh) | 一种基于serf原子磁强计的高空间分辨率心磁成像方法 | |
CN103142211A (zh) | 一种基于极值圆的心脏磁信号处理方法 | |
CN105640538A (zh) | 一种全张量心磁图仪探头及其制造方法 | |
JP4077945B2 (ja) | 生体磁気計測装置 | |
CN106132288A (zh) | 三维心脏轮廓重构方法 | |
Throne et al. | Fusion of body surface potential and body surface Laplacian signals for electrocardiographic imaging | |
CN113842147A (zh) | 一种基于原子蒸气室阵列的心/脑磁测量装置 | |
Deng et al. | 4D MRI using 3D radial sampling with respiratory self-gating to characterize temporal phase-resolved respiratory motion in the abdomen | |
JP2795212B2 (ja) | 生体磁気計測装置の較正方法 | |
Oyama et al. | Evaluation of an isosceles-triangle-coil phantom for magnetoencephalography | |
Oyama et al. | Calibration of a coil array geometry using an X-ray computed tomography |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |