CN105205227B - 一种基于最小vrms误差的梯度线圈设计方法 - Google Patents
一种基于最小vrms误差的梯度线圈设计方法 Download PDFInfo
- Publication number
- CN105205227B CN105205227B CN201510548202.4A CN201510548202A CN105205227B CN 105205227 B CN105205227 B CN 105205227B CN 201510548202 A CN201510548202 A CN 201510548202A CN 105205227 B CN105205227 B CN 105205227B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- vrms
- gradient
- gradient coil
- 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
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000003384 imaging method Methods 0.000 claims abstract description 31
- 238000004146 energy storage Methods 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 14
- 230000008571 general function Effects 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract 1
- 230000001788 irregular Effects 0.000 abstract 1
- 238000005070 sampling Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004949 mass spectrometry Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于最小VRMS误差的梯度线圈设计方法，其特征在于，梯度线圈的设计满足VRMS线性误差的平方与线圈储能的加权和构成的泛函最小。本发明同时给出了一种基于最小VRMS误差设计梯度线圈的算法流程。对于球形成像区域，本发明采用谐波系数展开来进行VRMS误差计算，具有很高的效率与精度。本发明适用于设计圆柱形梯度线圈、平板梯度线圈、不规则形状的梯度线圈等现存的任何类型的梯度线圈。
Description
技术领域
本发明涉及一种基于最小VRMS误差的梯度线圈设计方法，属于核磁共振成像系统部件设计领域。
背景技术
无论是永磁型MRI还是超导型MRI，其最核心的关键部件均为梯度线圈。衡量梯度线圈性能的标准，主要有场强、线性度、切换率等。而梯度线圈的设计主要是综合兼顾这三个指标。目前主流的梯度线圈设计算法，大多是在成像区域内或表面取一系列的采样点，然后令采样点处的梯度场与理想的梯度场之间的误差的平方和与线圈储能的和加权最小。这样做的缺点是，在算法上只保证了采样点处的线性误差最小，而不能保证在整个成像区域内线性误差最小。
发明内容
为了解决上述技术问题，本发明提供了一种基于最小VRMS误差的梯度线圈设计方法。
为了达到上述目的，本发明所采用的技术方案是：
一种基于最小VRMS误差的梯度线圈设计方法，包括以下步骤，
步骤一，预先给定梯度线圈所在的骨架尺寸以及预设参数；
预设参数包括成像区域范围、梯度场强度和VRMS误差Evrms；
VRMS误差Evrms为，
其中，V为成像区域的体积，Bz为梯度线圈产生的z方向磁场，G为梯度场的斜率，v为梯度场方向的坐标；
步骤二，构造流函数的基函数fi；
流函数的N项展开式为，
其中，ai为基函数系数，i为整数，i＝1,2...N；
步骤三，根据求得电流密度J的表达式；
其中，为骨架的外法向矢量；
步骤四，建立与电流密度J有关的泛函数Φ；
其中，w1,w2为权重系数，W为线圈储能；
步骤五，求解使泛函数取得极值的{ai}，得到流函数的表达式；
{ai}为由ai构成的列矢量；
步骤六，根据流函数表达式求得梯度线圈的形状。
对于球形成像区域，计算VRMS误差Evrms的方法为，
A1)计算梯度线圈在成像区域内的z方向磁场Bz的谐波系数An,m和Bn,m；
An,m和Bn,m由下式确定，
其中，为成像区域内任意一点P处的球坐标，r为球坐标半径，θ为P点与z轴的夹角，为P点在xy平面内的投影与x轴的夹角，Pn,m(cosθ)为n阶m次连带勒让德多项式；
A2)对谐波系数进行修正；
对于X线圈，令：
A11＝A11+G
对于Y线圈，令：
B11＝B11+G
对于Z线圈，令：
A10＝A10-G
其中，A11、B11、A1,0分别为X、Y、Z梯度线圈对应的主谐波系数，G为理想梯度磁场的斜率；
A3)计算出成像区域内的VRMS误差：
其中，R为成像区域半径。
梯度线圈产生的磁场B通过毕奥-萨伐尔定律求得：
其中，μ0为真空磁导率。
本发明所达到的有益效果：1、本发明能够确保在整个成像区域内梯度场的线性误差(线性度)最小，而非仅仅在个别采样点处的误差最小；2、本发明中将计算体积分误差的三维积分公式简化为用谐波系数表示的公式，能够快速精确计算成像区域内磁场的VRMS误差。
附图说明
图1为本发明的流程图。
图2为X线圈主线圈结构分布图。
图3为X线圈屏蔽线圈结构分布图。
图4为不同成像区域内梯度场的归一化Vrms误差。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
本发明受江苏省自然科学基金青年基金项目(项目批准号:BK20130854)支持。
如图1所示，一种基于最小VRMS误差的梯度线圈设计方法，包括以下步骤：
步骤一，预先给定梯度线圈所在的骨架尺寸以及预设参数；。
预设参数包括成像区域范围、梯度场强度和VRMS误差Evrms；
VRMS误差Evrms为，
其中，V为成像区域的体积，Bz为梯度线圈产生的z方向磁场，G为梯度场的斜率，v为梯度场方向的坐标。
对于球形成像区域，计算VRMS误差Evrms的方法为
A1)计算梯度线圈在成像区域内的z方向磁场Bz的谐波系数An,m和Bn,m；
采用毕奥-萨伐尔定律求梯度线圈产生的磁场B；
其中，μ0为真空磁导率；
成像区域内的z方向磁场Bz的An,m和Bn,m由下式确定，
其中，为成像区域内任意一点P处的球坐标，r为球坐标半径，θ为P点与z轴的夹角，为P点在xy平面内的投影与x轴的夹角，Pn,m(cosθ)为n阶m次连带勒让德多项式；
A2)对谐波系数进行修正；
对于X线圈，令：
A11＝A11+G
对于Y线圈，令：
B11＝B11+G
对于Z线圈，令：
A10＝A10-G
其中，A11、B11、A1,0分别为X、Y、Z梯度线圈对应的主谐波系数，G为理想梯度磁场的斜率；
A3)计算出成像区域内的VRMS误差：
其中，R为成像区域半径。
步骤二，构造流函数的基函数fi。
流函数的N项展开式为，
其中，ai为基函数系数，i为整数，i＝1,2...N。
步骤三，根据J＝▽×(nf)求得电流密度J的表达式；
其中，为骨架的外法向矢量；
步骤四，建立与电流密度J有关的泛函数Φ；
其中，w1,w2为权重系数，W为线圈储能。
步骤五，求解使泛函数取得极值的{ai}，得到流函数的表达式；
{ai}为由ai构成的列矢量。
步骤六，根据流函数表达式求得梯度线圈的形状。
根据上述方法，设计的一款带有源屏蔽的X方向梯度线圈的例子。
线圈的内径为330mm，外径为410mm，梯度场强为58uT/m/A，采用上述算法，设计得到的梯度线圈主线圈长度为114m，屏蔽线圈长度为77m，梯度线圈的总电感为350uH。所设计的梯度线圈的主线圈与屏蔽线圈结构如图2与图3所示。图4中给出了该梯度线圈在不同成像区域内的归一化VRMS误差。这里的归一化VRMS误差为用成像区域内的最大磁场对Evrms进行归一化，即：Evrms/(GR)。
通过上述方法，梯度线圈的形状可通过求解泛函的极值得到，所设计的梯度线圈能够使泛函取得最小值，从而保证在线圈储能一定的情况下，线圈的VRMS误差最小，确保在整个成像区域内梯度场的线性误差(线性度)最小，而非仅仅在个别采样点处的误差最小；同时上述方法中将计算体积分误差的三维积分公式简化为用谐波系数表示的公式，能够快速精确计算成像区域内磁场的VRMS误差；上述方法适用于设计圆柱形梯度线圈、平板梯度线圈、不规则形状的梯度线圈等现存的任何类型的梯度线圈。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。
Claims (2)
1.一种基于最小VRMS误差的梯度线圈设计方法，其特征在于：包括以下步骤，
步骤一，预先给定梯度线圈所在的骨架尺寸以及预设参数；
预设参数包括成像区域范围、梯度场强度和VRMS误差Evrms；
VRMS误差Evrms为，
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>v</mi>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<mi>V</mi>
</mfrac>
<munder>
<mo>&Integral;</mo>
<mi>V</mi>
</munder>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>B</mi>
<mi>z</mi>
</msub>
<mo>-</mo>
<mi>G</mi>
<mo>&CenterDot;</mo>
<mi>v</mi>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>V</mi>
<mo>&rsqb;</mo>
</mrow>
<mn>0.5</mn>
</msup>
</mrow>
其中，V为成像区域的体积，Bz为梯度线圈产生的z方向磁场，G为梯度场的斜率，v为梯度场方向的坐标；
步骤二，构造流函数的基函数fi；
流函数的N项展开式为，
<mrow>
<mi>f</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
</mrow>
其中，ai为基函数系数，i为整数，i＝1,2...N；
步骤三，根据求得电流密度J的表达式；
其中，为骨架的外法向矢量；
步骤四，建立与电流密度J有关的泛函数Φ；
<mrow>
<mi>&Phi;</mi>
<mo>=</mo>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
<mfrac>
<mn>1</mn>
<mi>V</mi>
</mfrac>
<munder>
<mo>&Integral;</mo>
<mi>V</mi>
</munder>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>B</mi>
<mi>z</mi>
</msub>
<mo>-</mo>
<mi>G</mi>
<mo>&CenterDot;</mo>
<mi>v</mi>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>V</mi>
<mo>+</mo>
<msub>
<mi>w</mi>
<mn>2</mn>
</msub>
<mi>W</mi>
</mrow>
其中，w1,w2为权重系数，W为线圈储能；
步骤五，求解使泛函数取得极值的{ai}，得到流函数的表达式；
{ai}为由ai构成的列矢量；
步骤六，根据流函数表达式求得梯度线圈的形状；
对于球形成像区域，计算VRMS误差Evrms的方法为，
A1)计算梯度线圈在成像区域内的z方向磁场Bz的谐波系数An,m和Bn,m；
An,m和Bn,m由下式确定，
其中， 为成像区域内任意一点P处的球坐标，r为球坐标半径，θ为P点与z轴的夹角，为P点在xy平面内的投影与x轴的夹角，Pn,m(cosθ)为n阶m次连带勒让德多项式；
A2)对谐波系数进行修正；
对于X线圈，令：
A11＝A11+G
对于Y线圈，令：
B11＝B11+G
对于Z线圈，令：
A10＝A10-G
其中，A11、B11、A1,0分别为X、Y、Z梯度线圈对应的主谐波系数，G为理想梯度磁场的斜率；
A3)计算出成像区域内的VRMS误差：
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>v</mi>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mn>3</mn>
<mn>2</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>&infin;</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>n</mi>
</munderover>
<mfrac>
<msup>
<mi>R</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>3</mn>
<mo>)</mo>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
<mo>!</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>m</mi>
<mo>)</mo>
<mo>!</mo>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>-</mo>
<msub>
<mi>&kappa;</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&kappa;</mi>
<mi>m</mi>
</msub>
<msubsup>
<mi>B</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
其中，R为成像区域半径。
2.根据权利要求1所述的一种基于最小VRMS误差的梯度线圈设计方法，其特征在于：梯度线圈产生的磁场B通过毕奥-萨伐尔定律求得：
其中，μ0为真空磁导率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510548202.4A CN105205227B (zh) | 2015-08-31 | 2015-08-31 | 一种基于最小vrms误差的梯度线圈设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510548202.4A CN105205227B (zh) | 2015-08-31 | 2015-08-31 | 一种基于最小vrms误差的梯度线圈设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105205227A CN105205227A (zh) | 2015-12-30 |
| CN105205227B true CN105205227B (zh) | 2018-05-08 |
Family
ID=54952906
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510548202.4A Active CN105205227B (zh) | 2015-08-31 | 2015-08-31 | 一种基于最小vrms误差的梯度线圈设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105205227B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106777900B (zh) * | 2016-11-25 | 2021-09-24 | 上海辰光医疗科技股份有限公司 | 一种高均匀度磁体的Vrms均匀度的单积分计算方法 |
| CN109885957A (zh) * | 2019-03-01 | 2019-06-14 | 中国科学院长春光学精密机械与物理研究所 | 一种用于生物医学成像系统的梯度线圈设计方法及梯度线圈 |
| CN109946630B (zh) * | 2019-03-27 | 2021-05-28 | 中国计量大学 | 最小化最大温度值的盘式梯度线圈设计方法 |
| CN109946631B (zh) * | 2019-03-27 | 2021-06-25 | 中国计量大学 | 最小化最大温度值的盘式梯度线圈 |
| CN113325350B (zh) * | 2020-02-28 | 2022-04-08 | 河海大学 | 一种基于离散网格的高性能梯度线圈设计方法 |
| CN111507039B (zh) * | 2020-04-26 | 2022-10-14 | 河海大学 | 一种切割铜皮式梯度线圈设计方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104007406A (zh) * | 2014-06-13 | 2014-08-27 | 河海大学 | 基于空间变换的椭圆柱面横向梯度线圈设计方法 |
| CN104062613A (zh) * | 2014-06-13 | 2014-09-24 | 河海大学 | 一种有源屏蔽梯度线圈及其设计方法 |
-
2015
- 2015-08-31 CN CN201510548202.4A patent/CN105205227B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104007406A (zh) * | 2014-06-13 | 2014-08-27 | 河海大学 | 基于空间变换的椭圆柱面横向梯度线圈设计方法 |
| CN104062613A (zh) * | 2014-06-13 | 2014-09-24 | 河海大学 | 一种有源屏蔽梯度线圈及其设计方法 |
Non-Patent Citations (2)
| Title |
|---|
| VRMS homogeneity definition：A Proposal;T.J.Havens etal;《proc inti soc resonmed》;20020131;全文 * |
| 超椭圆柱面梯度线圈设计;王亮等;《万方数据知识服务平台》;20150122;正文第3节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105205227A (zh) | 2015-12-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105205227B (zh) | 一种基于最小vrms误差的梯度线圈设计方法 | |
| CN104007406B (zh) | 基于空间变换的椭圆柱面横向梯度线圈设计方法 | |
| WO2011065357A1 (ja) | Ｍｒｉ装置用磁場調整方法 | |
| CN102096053A (zh) | 一种核磁共振系统中的梯度线圈设计法 | |
| CN105445683B (zh) | 一种柱面横向自屏蔽梯度线圈设计方法 | |
| Turek et al. | Magnetic field homogeneity perturbations in finite Halbach dipole magnets | |
| US20120235685A1 (en) | Gradient coil, magnetic resonance imaging device, and method for designing coil pattern | |
| CN102150222A (zh) | 用于限定用于成像体积的磁场的磁组件和方法 | |
| CN110568390A (zh) | 一种双平面匀场线圈的设计方法 | |
| CN101852843B (zh) | 一种超导磁体外磁屏蔽线圈的优化设计方法 | |
| Zhang et al. | An optimized target-field method for MRI transverse biplanar gradient coil design | |
| Zhao et al. | Analytical model for calculation of the magnetic field of triaxial uniform coils in magnetically shielded cylinder | |
| CN105068024B (zh) | Mri系统中梯度线圈产生的梯度磁场线性度衡量方法 | |
| CN112684392A (zh) | 一种非平面梯度线圈的设计方法 | |
| Puchalski et al. | Deuteron and triton magnetic moments from NMR spectra of the hydrogen molecule | |
| CN107205690A (zh) | 磁场调整方法 | |
| Liu et al. | An efficacious target-field approach to design shim coils for Halbach magnet of mobile NMR sensors | |
| CN104899441B (zh) | 一种适用于计算磁场vrms均匀度的数值方法 | |
| CN104833930A (zh) | 开放式磁共振系统梯度线圈磁场强度的计算方法 | |
| Zhang et al. | A spiral, bi-planar gradient coil design for open magnetic resonance imaging | |
| CN105184071A (zh) | 计算环形线圈在旋转对称区域内磁场vrms均匀度算法 | |
| JP6368849B2 (ja) | 磁場調整支援システムおよび磁場調整方法 | |
| Alnajjar et al. | An H-shaped low-field magnet for NMR spectroscopy designed using the finite element method | |
| Huang et al. | Study on magnetic field mapping method in the center volume of the air-core solenoid | |
| Du et al. | A novel design methodology for active shim coil |
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 |