CN104899441B - 一种适用于计算磁场vrms均匀度的数值方法 - Google Patents
一种适用于计算磁场vrms均匀度的数值方法 Download PDFInfo
- Publication number
- CN104899441B CN104899441B CN201510301197.7A CN201510301197A CN104899441B CN 104899441 B CN104899441 B CN 104899441B CN 201510301197 A CN201510301197 A CN 201510301197A CN 104899441 B CN104899441 B CN 104899441B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- vrms
- tetrahedron
- uniformitys
- point
- 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
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种适用于计算磁场VRMS均匀度的数值方法,其特征是,包括如下步骤:步骤一)将待计算的区域离散成四面体网格;步骤二)将四面体内的磁场用基函数插值表示;步骤三)确定每个四面体内的插值点,并根据插值点处的磁场计算出插值基函数的系数;步骤四)将磁场基函数展开式带入积分公式计算成像区域内磁场的VRMS值。本发明所达到的有益效果:通过本发明给出的计算方法,能够精确计算任意区域的磁场均匀度,具有很高的适用范围与计算精度,可以计算包括但不限于球、椭球、圆柱、立方体等成像区域内的磁场。算法的精度可以通过调整网格的尺寸来得到。
Description
技术领域
本发明涉及一种适用于计算磁场VRMS均匀度的数值方法,属于核磁共振成像系统部件仿真与设计技术领域。
背景技术
核磁共振成像(MRI)是一种利用核磁共振原理来进行成像的医学影像新技术。近年来,随着科技的发展,MRI系统逐渐由永磁型向超导型过渡。无论是永磁型MRI还是超导型MRI,其最核心的关键部件均为磁体。而衡量磁体性能最主要的指标是磁场均匀度。过去衡量MRI主磁场均匀度通常采用峰峰值来进行计算。所谓峰峰值,即成像区域内最大场与最小场之差与中心场之比。由于整个成像区域内的最大场与最小场很难求得,通常的做法是在成像区域内取一系列的采样点,然后根据所有采样点上的最大值与最小值计算峰峰值。这样做有两个缺点,一是采样点上的最大值与最小值与整个区域内的最大值与最小值可能相差很大,二是不能描述整个成像区域内的均匀度。因此,采用体均方根(VRMS)均匀度代替峰峰值均匀度是一种更合理的做法。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种基于VRMS表示的磁场均匀度计算方法,以解决现有技术中缺少精确、有效、适用性广的磁场均匀度度量方法的问题。
为了实现上述目标,本发明采用如下的技术方案:
一种适用于计算磁场VRMS均匀度的数值方法,其特征是,包括如下步骤:
1)将待计算的区域离散成N个四面体网格,计算每个四面体的体积,第i个四面体内的区域用Ωi表示,体积用Vi表示;
2)将每个四面体内的磁场用基函数L插值表示;
3)根据每个四面体内的插值点处的磁场计算出插值基函数的系数;
4)计算成像区域内磁场的VRMS值:将基函数代入如下积分公式
这里为整个计算区域的体积,B0为主磁场的中心场,Bz(r)为磁体在r处产生的z方向磁场。
前述的一种适用于计算磁场VRMS均匀度的数值方法,其特征是,所述步骤1)中第i个四面体内的第j个点的基函数Lij定义如下:Lij(r)=Vij(r)/Vi,(1≤j≤4),这里r为四面体内任一点P的矢量坐标,Vij(r)为该点与四面体除第j个点之外的三个顶点组成的四面体的体积。
前述的一种适用于计算磁场VRMS均匀度的数值方法,其特征是,所述步骤1)中,离散网格尺寸不超过10mm。
前述的一种适用于计算磁场VRMS均匀度的数值方法,其特征是,所述步骤2)中四面体内的插值点在四面体的顶点。
前述的一种适用于计算磁场VRMS均匀度的数值方法,其特征是,所述步骤3)中,当涉及超导磁体问题时插值点处的磁场用毕奥-萨伐尔定律求得:这里μ0为自由空间的磁导率,I为超导磁体导线上的电流,为观察点的坐标矢量。
本发明所达到的有益效果:通过本发明给出的计算方法,能够精确计算任意区域的磁场均匀度,具有很高的适用范围与计算精度,可以计算包括但不限于球、椭球、圆柱、立方体等成像区域内的磁场。算法的精度可以通过调整网格的尺寸来得到。
附图说明
图1是所设计超导磁体线圈磁场峰均值为3ppm的等高线分布;
图2是所设计超导磁体线圈磁场峰均值为1ppm的等高线分布;
图3是椭球成像区域的四面体离散网格;
图4计算磁场VRMS均匀度的算法流程图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
本申请受江苏省自然科学基金青年基金项目(项目批准号:BK20130854)支持。
现有技术中基本都是采用传统峰峰值均匀度来度量MRI主磁体产生的磁场均匀度,下面举例来进行说明,首先设计一个超导磁体线圈,然后计算该线圈在空间产生的磁场均匀度分布。线圈设计要求为在45cm DSV的球形成像空间内峰峰值均匀度为6ppm,采样点为在球面上取24个高斯分布点。图1所示为所设计的超导磁体线圈磁场峰均值为3ppm的等高线分布。图中的圆弧半径为r=22.5cm。从图中可以看出,在球面一些地方不满足设计要求。图2所示为所设计的超导磁体线圈磁场峰均值为1ppm的等高线分布。图中的曲线有大幅度的凹陷。可以看出,采用峰峰值来度量均匀度的方法具有很大的局限性,不能反映整个成像空间的磁场分布。
而在本发明中采用VRMS表示成像空间内磁场均匀度的做法,适用于超导磁体。
以下结合附图和具体实施例对本发明作具体的介绍。
本实施例结合超导磁体在椭球体成像空间内产生的磁场来说明VRMS均匀度的计算。在MRI系统中,关心的是球体内或椭球体的主磁体磁场均匀度。为了计算VRMS均匀度,采用如下的步骤:
步骤一:将待计算的区域离散成四面体网格。计算每个四面体的体积,第i个四面体内的区域用Ωi表示,体积用Vi表示。
步骤二:将四面体内的磁场用基函数L插值表示:
本实施例中,第i个四面体内的第j个点基函数Lij定义如下:Lij(r)=Vij(r)/Vi,(1≤j≤4)这里r为四面体内任一点P的矢量坐标,Vij(r)为该点与除j点之外的三个顶点组成的四面体的体积。采用该基函数的好处是,在第i个四面体内,Lij(r)在顶点j处为1,在其他顶点处为零。因此其系数代表四面体顶点处的磁场值。
步骤三:确定每个四面体内的插值点,并计算每个插值点处的磁场。根据上面的分析,基函数的系数代表四面体顶点处的磁场值。这样插值点即为四个顶点。
只需求出四面体i四个顶点处的磁场即可将整个四面体内的磁场用如下的插值函数表示:
步骤四:将基函数代入如下积分公式计算成像区域内磁场的VRMS值:
这里为整个计算区域的体积,B0为主磁场的中心场。Bz(r)采用基函数展开,并对上式进行积分,即可得到整个成像区域内的VRMS均匀度。
对于超导磁体线圈产生的磁场,插值点处的磁场可以用毕奥-萨伐尔定律求得:在超导磁体线圈结构已知时,可以采用上述公式计算各个顶点处的磁场值。如果线圈结构未知,则可以采用霍尔探头测量成像区域表面采样点上的磁场值,然后通过拟合求各个顶点处的磁场。这时候,其精度受到两方面的限制,一个是网格尺寸,一个是采样点个数。
下面给出前面设计的超导磁体线圈在不同成像区域内的峰峰值均匀度以及采用本发明中的算法计算的VRMS均匀度。Rxy表示xy平面内的成像区域半径,Rz表示z方向的成像区域半径。为了计算方便,下面的计算结果中我们首先建立一个半径为100cm的球体,采用约5cm的网格尺寸进行离散,总共得到271684个四面体,网格中共包括48044个节点,如图3所示。实际的成像区域为对该球体网格进行压缩得到。最大的网格尺寸约为1.1cm。
表1.不同成像区域的磁场Vrms均匀度与峰峰值均匀度(pk to pk)
Rxy(cm) | Rz(cm) | Vrms(ppm) | pk to pk |
5 | 5 | 0.001686 | 0.007583 |
10 | 10 | 0.015343 | 0.116245 |
15 | 15 | 0.086042 | 0.686476 |
20 | 20 | 0.426208 | 2.769701 |
22.5 | 22.5 | 0.964443 | 6.015614 |
25 | 25 | 2.277163 | 21.03525 |
20 | 17.5 | 0.320211 | 1.998711 |
22.5 | 20 | 0.773082 | 5.179824 |
采用该算法的好处是算法适用性比较广。可以计算包括但不限于球、椭球、圆柱、立方体等成像区域内的磁场。算法的精度可以通过调整网格的尺寸来得到。与传统的峰峰值均匀度度量方式相比,VRMS均匀度更能反映整个成像空间内的磁场均匀度,是一种更合理的度量方法。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (5)
1.一种适用于计算磁场VRMS均匀度的数值方法,其中VRMS表示体均方根,其特征是,包括如下步骤:
1)将待计算的区域离散成N个四面体网格,计算每个四面体的体积,第i个四面体内的区域用Ωi表示,体积用Vi表示;
2)将每个四面体内的磁场用基函数L插值表示;
3)根据每个四面体内的插值点处的磁场计算出插值基函数的系数;
4)计算成像区域内磁场的VRMS值:将基函数代入如下积分公式
这里为整个计算区域的体积,B0为主磁场的中心场,Bz(r)为磁体在r处产生的z方向磁场。
2.根据权利要求1所述的一种适用于计算磁场VRMS均匀度的数值方法,其特征是,所述步骤1)中第i个四面体内的第j个点的基函数Lij定义如下:Lij(r)=Vij(r)/Vi,其中1≤j≤4,r为四面体内任一点P的矢量坐标,Vij(r)为该点与四面体除第j个点之外的三个顶点组成的四面体的体积。
3.根据权利要求1所述的一种适用于计算磁场VRMS均匀度的数值方法,其特征是,所述步骤1)中,平均离散网格尺寸不超过10mm。
4.根据权利要求1所述的一种适用于计算磁场VRMS均匀度的数值方法,其特征是,所述步骤2)中四面体内的插值点在四面体的顶点。
5.根据权利要求1所述的一种适用于计算磁场VRMS均匀度的数值方法,其特征是,所述步骤3)中,插值点处的磁场用毕奥-萨伐尔定律求得:其中,μ0为自由空间的磁导率,I为超导磁体导线上的电流,为观察点的坐标矢量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510301197.7A CN104899441B (zh) | 2015-06-05 | 2015-06-05 | 一种适用于计算磁场vrms均匀度的数值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510301197.7A CN104899441B (zh) | 2015-06-05 | 2015-06-05 | 一种适用于计算磁场vrms均匀度的数值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104899441A CN104899441A (zh) | 2015-09-09 |
CN104899441B true CN104899441B (zh) | 2018-02-27 |
Family
ID=54032103
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510301197.7A Active CN104899441B (zh) | 2015-06-05 | 2015-06-05 | 一种适用于计算磁场vrms均匀度的数值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104899441B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105718729B (zh) * | 2016-01-20 | 2018-07-24 | 河海大学 | 一种柱面轴向梯度线圈设计中磁场与电感值的计算方法 |
CN107024723B (zh) * | 2017-06-16 | 2018-09-07 | 桂林理工大学 | 一种二度体磁场数值计算方法 |
CN112036017B (zh) * | 2020-08-21 | 2022-09-23 | 河海大学 | 一种螺线管线圈磁场快速计算方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1702647A (zh) * | 2005-07-25 | 2005-11-30 | 王华国 | 插值计算的有限单元法 |
US9588200B2 (en) * | 2011-03-25 | 2017-03-07 | Hitachi, Ltd. | Method for adjusting static magnetic field homogeneity, static magnetic field generation device for magnetic resonance imaging, magnetic field adjustment system, and program |
CN102393965B (zh) * | 2011-06-01 | 2014-01-29 | 博雅智汇国际有限公司 | 一种利用四面体插值法实现测色色差计颜色校正的方法 |
CN102707250B (zh) * | 2012-05-11 | 2014-04-02 | 中国科学院电工研究所 | 一种磁共振成像超导磁体系统的设计方法 |
-
2015
- 2015-06-05 CN CN201510301197.7A patent/CN104899441B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104899441A (zh) | 2015-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108139453A (zh) | 包含海尔贝克型圆柱环的核磁共振用磁性组件 | |
Martens et al. | Insertable biplanar gradient coils for magnetic resonance imaging | |
CN104899441B (zh) | 一种适用于计算磁场vrms均匀度的数值方法 | |
Turek et al. | Magnetic field homogeneity perturbations in finite Halbach dipole magnets | |
Lu et al. | Hybrid optimal design of square highly uniform magnetic field coils | |
CN109856575B (zh) | 一种确定磁共振被动匀场的方法、装置、设备及存储介质 | |
CN105068024B (zh) | Mri系统中梯度线圈产生的梯度磁场线性度衡量方法 | |
Li et al. | Design of highly uniform three-dimensional square magnetic field coils for external magnetic shielding of magnetometers | |
CN107144802B (zh) | 一种快速测量超导薄膜平均临界电流的方法 | |
McGinley et al. | A permanent MRI magnet for magic angle imaging having its field parallel to the poles | |
CN113514789B (zh) | 一种磁传感器阵列校准方法 | |
Ilmoniemi | The triangle phantom in magnetoencephalography | |
CN107205690A (zh) | 磁场调整方法 | |
CN105184071B (zh) | 计算环形线圈在旋转对称区域内磁场vrms均匀度算法 | |
CN104714201B (zh) | 一种有效矫正磁共振成像系统的主磁场的方法 | |
Wu et al. | Algorithm Research on the Conductor Eccentricity of a Circular Dot Matrix Hall High Current Sensor for ITER | |
Değirmenci et al. | Anisotropic conductivity imaging with MREIT using equipotential projection algorithm | |
CN109188320A (zh) | 一种基于磁阻效应的流场成像系统及成像方法 | |
Wang et al. | Design of compact self-shielded uniform magnetic field coils based on multipole moment optimization | |
Oraon et al. | Simulation-based design of optimized symmetric MRI magnet modified with ferromagnetic shell | |
Kramer et al. | Mapping flux trapping in SRF cavities to analyze the impact of geometry | |
RU2307370C1 (ru) | Способ определения магнитного момента квадратной катушки с током | |
Huang et al. | A study on the effect of the condition number in the magnetic field mapping of the Air-Core solenoid | |
Aboussouan et al. | High-precision absolute positioning of medical instruments in MRI systems | |
Chen | Demagnetizing factors of long cylinders with infinite susceptibility |
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 |