CN101963656B - 一种附加磁场检测方法和装置 - Google Patents
一种附加磁场检测方法和装置 Download PDFInfo
- Publication number
- CN101963656B CN101963656B CN 200910157757 CN200910157757A CN101963656B CN 101963656 B CN101963656 B CN 101963656B CN 200910157757 CN200910157757 CN 200910157757 CN 200910157757 A CN200910157757 A CN 200910157757A CN 101963656 B CN101963656 B CN 101963656B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- free induction
- induction decay
- phase
- change condition
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 11
- 230000000295 complement effect Effects 0.000 title abstract 6
- 230000006698 induction Effects 0.000 claims abstract description 60
- 238000000034 method Methods 0.000 claims abstract description 20
- 230000008859 change Effects 0.000 claims description 55
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 4
- 230000004069 differentiation Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical group [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种附加磁场检测方法,包括:生成被测梯度脉冲,并在所述被测梯度脉冲结束后,获取两个自由感应衰减信号;其中,第一个自由感应衰减信号同时受到主磁场和附加磁场的影响,第二个自由感应衰减信号只受到主磁场的影响;通过分析两个自由感应衰减信号的相位差,确定附加磁场的变化情况。本发明同时公开了一种附加磁场检测装置。应用本发明所述的方法和装置,能够有效地检测出附加磁场的变化情况。
Description
技术领域
本发明涉及磁共振技术,特别涉及一种磁共振系统中的附加磁场检测方法和装置。
背景技术
在磁共振系统中,由梯度脉冲所引起的磁场会附加在主磁场上,并衰减足够长的时间,因此,可将这种由梯度脉冲所引起的磁场称为附加磁场。附加磁场的存在会对磁共振系统中采集到的图像的图像质量造成严重影响,因此,希望能够检测出附加磁场的变化情况,以便根据其变化情况对所采集的图像进行调整,以得到质量更好的图像,但现有技术中还没有一种有效的附加磁场检测技术,来有效地检测出附加磁场的变化情况。
发明内容
有鉴于此,本发明的主要目的在于提供一种附加磁场检测方法,能够有效地检测出附加磁场的变化情况。
本发明的另一目的在于提供一种附加磁场检测装置,能够有效地检测出附加磁场的变化情况。
为达到上述目的,本发明的技术方案是这样实现的:
一种附加磁场检测方法,包括:
生成被测梯度脉冲,并在所述被测梯度脉冲结束后,获取两个自由感应衰减信号;其中,第一个自由感应衰减信号同时受到主磁场和附加磁场的影响,第二个自由感应衰减信号只受到主磁场的影响;
通过分析两个自由感应衰减信号的相位差,确定附加磁场的变化情况。
其中,所述获取两个自由感应衰减信号包括:
施加两个射频脉冲;
在每个射频脉冲结束后,分别采集得到一个自由感应衰减信号。
较佳地,所述采集得到自由感应衰减信号包括:利用模数转换器采集得到自由感应衰减信号。
较佳地,所述通过分析两个自由感应衰减信号的相位差,确定附加磁场的变化情况包括:
利用所述第一个自由感应衰减信号的相位计算得到由主磁场和附加磁场所组成的混合磁场的变化情况;
利用所述第二个自由感应衰减信号的相位计算得到主磁场的变化情况;
通过比较所述混合磁场的变换情况以及所述主磁场的变化情况,确定所述附加磁场的变化情况。
其中,所述利用第一个自由感应衰减信号的相位计算得到由主磁场和附加磁场所组成的混合磁场的变化情况包括:利用以下公式计算得到混合磁场的变化情况:其中,所述表示所述第一个自由感应衰减信号的相位的一阶导数,所述γ表示氢质子的磁旋比,所述Bm(t)表示混合磁场的变化情况;
所述利用第二个自由感应衰减信号的相位计算得到主磁场的变化情况包括:利用以下公式计算得到主磁场的变化情况:其中,所述表示所述第二个自由感应衰减信号的相位的一阶导数,所述γ表示氢质子的磁旋比;所述B0(t)表示主磁场的变化情况。
所述h表示模数转换器的采样间隔。
一种附加磁场检测装置,包括:
一个生成获取单元,用于生成被测梯度脉冲,并在所述被测梯度脉冲结束后,获取两个自由感应衰减信号;其中,第一个自由感应衰减信号同时受到主磁场和附加磁场的影响,第二个自由感应衰减信号只受到主磁场的影响;
一个分析单元,用于通过分析两个自由感应衰减信号的相位差,确定附加磁场的变化情况。
其中,所述生成获取单元包括:
一个生成子单元,用于生成被测梯度脉冲;
一个获取子单元,用于在所述被测梯度脉冲结束后,施加两个射频脉冲,并在每个射频脉冲结束后,利用模数转换器分别采集得到一个自由感应衰减信号。
所述分析单元包括:
一个计算子单元,用于利用所述第一个自由感应衰减信号的相位计算得到由主磁场和附加磁场所组成的混合磁场的变化情况,利用所述第二个自由感应衰减信号的相位计算得到主磁场的变化情况;
一个确定子单元,用于通过比较所述混合磁场的变换情况以及所述主磁场的变化情况,确定所述附加磁场的变化情况。
较佳地,所述计算子单元包括:
可见,采用本发明的技术方案,在利用被测梯度脉冲产生附加磁场后,获取两个独立的FID信号,其中,第一个FID信号同时受到主磁场和附加磁场的影响,而第二个FID信号只受到主磁场的影响,通过对这两个FID信号的相位进行分析,来确定附加磁场的变化情况,所述本发明所述方案能够有效地检测出附加磁场的变化情况,进而根据检测出的附加磁场变化情况改善采集所得图像的质量。
附图说明
下面将通过参照附图详细描述本发明的优选实施例,使本领域的普通技术人员更清楚本发明的上述及其它特征和优点,附图中:
图1为本发明方法实施例的流程图;
图2为本发明方法实施例中的FID信号获取过程示意图;
图3为本发明装置实施例的组成结构示意图;
图4为本发明实施例中一个ADC周期中的第一和第二个FID信号的相位示意图;
图5为本发明实施例中从行角度所显示的附加磁场变化情况示意图;
图6为本发明实施例中从列角度所显示的附加磁场变化情况示意图;
图7为图5和图6所示附加磁场变化情况的频谱示意图。
具体实施方式
在磁共振系统中,横向磁化矢量形成后,将以共振频率自由振动并指数衰减,在此过程中产生的信号被称为自由感应衰减(FID,Free Induction Decay)信号。共振频率通常仅与主磁场的强度有关,即ω=γB0(其中γ表示氢质子的磁旋比),因此,FID信号的相位可表示为:φ=ωt+φi,其中,φi为初始相位。但是,在实际应用中,附加磁场也会对FID信号的相位造成影响。本发明所述方案中,利用双FID序列来生成两个独立的FID信号,其中,第一个FID信号会同时受到主磁场和附加磁场的影响,而第二个FID信号只会受到主磁场的影响,通过对这两个FID信号的相位进行分析,来确定附加磁场的变化情况。
为了使本发明的技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
图1为本发明方法实施例的流程图。如图1所示,包括以下步骤:
步骤101:生成被测梯度脉冲。
由于附加磁场是由梯度脉冲所引起的,因此,本实施例中,利用被测梯度脉冲来模拟附加磁场的产生情况。
图2为本发明方法实施例中的FID信号获取过程示意图。如图2所示,其中,被测梯度脉冲(Gtest)的数量、极性、振幅、上升时间(Tr)以及顶部时间(Tt)等均可根据系统的实际情况而定,即使得利用被测梯度脉冲产生的附加磁场与系统实际的附加磁场相同。
步骤102:在被测梯度脉冲结束后,获取两个FID信号;其中,第一个FID信号同时受到主磁场和附加磁场的影响,第二个FID信号只受到主磁场的影响。
本步骤中,在被测梯度脉冲结束后,施加两个翻转角可调的射频(RF,RadioFrequency)脉冲,即图2所示α;在每个RF脉冲结束后,利用模数转换器(ADC,Analog-to-Digital Converter)分别采集得到一个FID。如何施加RF脉冲以及如何采集得到FID信号均为本领域公知,不再赘述。
可以看出,第一个FID信号是在梯度脉冲结束后的一个可调整的较短时间后产生的,因此该信号会同时受到主磁场和附加磁场的影响,而第二个FID信号是在梯度脉冲结束后的一个固定的较长时间后产生的,该固定较长时间的具体取值可根据实际需要而定,但需要保证此时附加磁场已经衰减为零,即第二个FID信号不会受到附加磁场的影响。所以,可将第二个FID信号作为第一个FID信号的参考基准,两个信号的相位差可用于说明附加磁场的变化情况。上述所提到的可调整的较短时间,即指图2所示延迟时间(TD,Delay Time),TD的具体取值是可调的,也就是说,第一个FID信号的起始时间是可调的,因此,可识别出一个较长时间周期内的附加磁场的变化情况。
步骤103:通过分析两个自由感应衰减信号的相位差,确定附加磁场的变化情况。
本步骤中,首先利用第一个FID信号的相位计算得到由主磁场和附加磁场所组成的混合磁场的变化情况;然后,利用第二个FID信号的相位计算得到主磁场的变化情况;最后,通过比较混合磁场的变换情况以及主磁场的变化情况,确定附加磁场的变化情况。具体实现过程如下:
其中,γ表示氢质子的磁旋比,t表示时间。需要说明的是,为了区分两个不同的时间变量t和T,公式(1)中将Bm(t)表示为了Bm(T)。
基于公式(1),可推导出混合磁场的变化情况为:
然后,将混合磁场的变化情况Bm(t)和主磁场的变化情况B0(t)进行比较,以确定附加磁场的变化情况。
基于上述方法,图3为本发明装置实施例的组成结构示意图。如图3所示,附加磁场检测装置包括:一个生成获取单元31和一个分析单元32。
生成获取单元31生成被测梯度脉冲,并在被测梯度脉冲结束后,获取两个自由感应衰减信号;其中,第一个自由感应衰减信号同时受到主磁场和附加磁场的影响,第二个自由感应衰减信号只受到主磁场的影响;分析单元32通过分析两个自由感应衰减信号的相位差,确定附加磁场的变化情况。
其中,生成获取单元31中包括:一个生成子单元311和一个获取子单元312。生成子单元311生成被测梯度脉冲;获取子单元312在被测梯度脉冲结束后,施加两个射频脉冲,并在每个射频脉冲结束后,利用模数转换器分别采集得到一个自由感应衰减信号。
分析单元32中包括:一个计算子单元321和一个确定子单元322。计算子单元321利用第一个自由感应衰减信号的相位计算得到由主磁场和附加磁场所组成的混合磁场的变化情况,并利用第二个自由感应衰减信号的相位计算得到主磁场的变化情况;确定子单元322通过比较混合磁场的变化情况以及主磁场的变化情况,确定附加磁场的变化情况。
其中,计算子单元321中又可进一步包括:一个第一处理子单元3211和一个第二处理子单元3212。其中,第一处理子单元3211利用以下公式计算得到混合磁场的变化情况:其中,表示第一个自由感应衰减信号的相位的一阶导数,γ表示氢质子的磁旋比,Bm(t)表示混合磁场的变化情况;第二处理子单元3212利用以下公式计算得到主磁场的变化情况:其中,表示所述第二个自由感应衰减信号的相位的一阶导数,γ表示氢质子的磁旋比,B0(t)表示主磁场的变化情况。
图3所示装置的具体工作流程请参照图1所示方法实施例中的相应说明,此处不再赘述。
下面通过具体的示例,对本发明所述方案作进一步地详细说明。
图4为本发明实施例中一个ADC周期中的第一个FID信号和第二个FID信号的相位示意图。由于第一个FID信号会受到附加磁场的影响,因此其相位会出现振荡,如图4中的曲线1所示;而第二个FID信号只会受到主磁场的影响,因此其相位不会出现振荡,如图4中直线2所示;两者的相位差为一平滑的曲线,如图4中的曲线3所示,其上的振荡仅由附加磁场变化所引起。
根据图4所示相位差,可方便地利用微分方法计算出附加磁场的变化情况,如公式(2)所示。由于附加磁场为一平面场,因此图5和图6分别从行和列两个不同的角度显示出了附加磁场的变化情况,其中的Be(t)即表示附加磁场。图7为图5和图6所示附加磁场变化情况的频谱示意图。从上述图5、6和7中,可以很清楚地看出附加磁场的变化情况。
总之,采用本发明的技术方案,能够有效检测出附加磁场的变化情况。
需要说明的是,上述实施例仅用于举例说明,并不用于限制本发明的技术方案。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种附加磁场检测方法,该方法包括:
生成被测梯度脉冲,并在所述被测梯度脉冲结束后,获取两个自由感应衰减信号;其中,第一个自由感应衰减信号同时受到主磁场和附加磁场的影响,第二个自由感应衰减信号只受到主磁场的影响;
通过分析两个自由感应衰减信号的相位差,确定附加磁场的变化情况,其中,确定附加磁场的变化情况包括:
利用所述第一个自由感应衰减信号的相位计算得到由主磁场和附加磁场所组成的混合磁场的变化情况;
利用所述第二个自由感应衰减信号的相位计算得到主磁场的变化情况;
通过比较所述混合磁场的变换情况以及所述主磁场的变化情况,确定所述附加磁场的变化情况。
2.根据权利要求1所述的方法,其特征在于,所述获取两个自由感应衰减信号包括:
施加两个射频脉冲;
在每个射频脉冲结束后,分别采集得到一个自由感应衰减信号。
3.根据权利要求2所述的方法,其特征在于,利用模数转换器采集自由感应衰减信号。
6.一种附加磁场检测装置,包括:
一个生成获取单元(31),用于生成被测梯度脉冲,并在所述被测梯度脉冲结束后,获取两个自由感应衰减信号;其中,第一个自由感应衰减信号同时受到主磁场和附加磁场的影响,第二个自由感应衰减信号只受到主磁场的影响;
一个分析单元(32),用于通过分析两个自由感应衰减信号的相位差,确定附加磁场的变化情况。
7.根据权利要求6所述的装置,其特征在于,所述生成获取单元(31)包括:
一个生成子单元(311),用于生成被测梯度脉冲;
一个获取子单元(312),用于在所述被测梯度脉冲结束后,施加两个射频脉冲,并在每个射频脉冲结束后,分别采集得到一个自由感应衰减信号。
8.根据权利要求6所述的装置,其特征在于,所述分析单元(32)包括:
一个计算子单元(321),用于利用所述第一个自由感应衰减信号的相位计算得到由主磁场和附加磁场所组成的混合磁场的变化情况,利用所述第二个自由感应衰减信号的相位计算得到主磁场的变化情况;
一个确定子单元(322),用于通过比较所述混合磁场的变换情况以及所述主磁场的变化情况,确定所述附加磁场的变化情况。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910157757 CN101963656B (zh) | 2009-07-24 | 2009-07-24 | 一种附加磁场检测方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910157757 CN101963656B (zh) | 2009-07-24 | 2009-07-24 | 一种附加磁场检测方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101963656A CN101963656A (zh) | 2011-02-02 |
CN101963656B true CN101963656B (zh) | 2012-12-19 |
Family
ID=43516611
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910157757 Active CN101963656B (zh) | 2009-07-24 | 2009-07-24 | 一种附加磁场检测方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101963656B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155621B (zh) * | 2014-09-11 | 2017-01-18 | 中国科学院武汉物理与数学研究所 | 一种准确测量静磁场b0分布的方法 |
CN106772648A (zh) * | 2016-12-05 | 2017-05-31 | 中国地质大学(武汉) | 提高质子旋进类传感器调谐精度和信噪比的方法及系统 |
CN109946632B (zh) * | 2019-03-06 | 2023-10-20 | 佛山瑞加图医疗科技有限公司 | 剩磁相位差测量方法及剩磁测量方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101144852A (zh) * | 2006-09-15 | 2008-03-19 | 西门子(中国)有限公司 | 磁场分布的测试方法及其装置 |
CN101328805A (zh) * | 2007-06-19 | 2008-12-24 | 普拉德研究及开发股份有限公司 | 测量自由感应衰减信号的方法和设备及其在分析中的应用 |
-
2009
- 2009-07-24 CN CN 200910157757 patent/CN101963656B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101144852A (zh) * | 2006-09-15 | 2008-03-19 | 西门子(中国)有限公司 | 磁场分布的测试方法及其装置 |
CN101328805A (zh) * | 2007-06-19 | 2008-12-24 | 普拉德研究及开发股份有限公司 | 测量自由感应衰减信号的方法和设备及其在分析中的应用 |
Also Published As
Publication number | Publication date |
---|---|
CN101963656A (zh) | 2011-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107478402B (zh) | 一种低频振动位移检测系统和方法 | |
CA2713772C (en) | A method and apparatus for phase sensitive detection of eddy current measurements | |
US10557908B2 (en) | Magnetic field monitoring of spiral echo train imaging | |
JP2018518322A5 (zh) | ||
CN112881959B (zh) | 一种用于磁共振成像的梯度涡流补偿方法及系统 | |
CN101963656B (zh) | 一种附加磁场检测方法和装置 | |
CN104995530B (zh) | 被检体信息获取装置和被检体信息获取方法 | |
RU2012157092A (ru) | Магнитно-резонансное обследование с обнаружением инструмента | |
CN108759658A (zh) | 一种感应同步器测角误差分析与补偿方法 | |
CN104155621B (zh) | 一种准确测量静磁场b0分布的方法 | |
JP2004024669A5 (zh) | ||
CN108627788B (zh) | 一种霍尔电流传感器磁滞检测的方法及系统 | |
CN108426909B (zh) | 核磁共振弛豫时间和流速同步快速测量的方法及装置 | |
US20170135590A1 (en) | Magnetic resonance imaging apparatus and blood flow drawing method | |
Marszalek | Maxwell-Wien bridge with vector voltmeter system for measurement small and rapid changes in inductive-loop sensor impedance components | |
CN107861033B (zh) | 振荡波局部放电检测系统定位误差的检定方法和系统 | |
JP2016067895A (ja) | 磁気共鳴イメージング装置及び磁気共鳴イメージング方法 | |
JP6516063B2 (ja) | 測定装置および材料試験機 | |
CN104931731B (zh) | 一种用于检测机车电涡流速度传感脉冲发生器质量的方法 | |
JP6152806B2 (ja) | 生体情報測定方法 | |
CN110133734B (zh) | 一种用于核磁共振测井的信号检测方法 | |
WO2012139680A1 (en) | A method and system for determining the flip angle in the transient phase of a gradient echo pulse sequence | |
JP2011038969A (ja) | インピーダンス測定装置 | |
CN107219393B (zh) | 一种信号功率检测方法、装置和设备 | |
JP5855382B2 (ja) | コンデンサの絶縁抵抗測定装置およびコンデンサの絶縁抵抗測定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |