一种基于最小vrms误差的梯度线圈设计方法
【技术领域】
[0001] 本发明涉及一种基于最小VRMS误差的梯度线圈设计方法，属于核磁共振成像系 统部件设计领域。
【背景技术】
[0002]无论是永磁型MRI还是超导型MRI，其最核心的关键部件均为梯度线圈。衡量梯度 线圈性能的标准，主要有场强、线性度、切换率等。而梯度线圈的设计主要是综合兼顾这三 个指标。目前主流的梯度线圈设计算法，大多是在成像区域内或表面取一系列的采样点，然 后令采样点处的梯度场与理想的梯度场之间的误差的平方和与线圈储能的和加权最小。这 样做的缺点是，在算法上只保证了采样点处的线性误差最小，而不能保证在整个成像区域 内线性误差最小。
【发明内容】
[0003] 为了解决上述技术问题，本发明提供了一种基于最小VRMS误差的梯度线圈设计 方法。
[0004] 为了达到上述目的，本发明所采用的技术方案是：
[0005] -种基于最小VRMS误差的梯度线圈设计方法，包括以下步骤，
[0006] 步骤一，预先给定梯度线圈所在的骨架尺寸以及预设参数；
[0007] 预设参数包括成像区域范围、梯度场强度和VRMS误差Evms;
[0008] VRMS误差Ev"s为，
[0009]
[0010] 其中，V为成像区域的体积，Bz为梯度线圈产生的z方向磁场，G为梯度场的斜率， v为梯度场方向的坐标；
[0011] 步骤二，构造流函数的基函数f1;
[0012] 流函数的N项展开式为，
[0013]
/=1
[0014] 其中，3;为基函数系数，i为整数，i=1，2... N ;
[0015] 步骤三，根据/ = ￥ X (gf )求得电流密度J的表达式；
[0016] 其中，g为骨架的外法向矢量；
[0017] 步骤四，建立与电流密度J有关的泛函数〇;
[0018]
[0019] 其中，Wl，w2为权重系数，W为线圈储能；
[0020]步骤五，求解使泛函数取得极值的{aj，得到流函数的表达式；
[0021] {aj为由&1构成的列矢量；
[0022] 步骤六，根据流函数表达式求得梯度线圈的形状。
[0023] 对于球形成像区域，计算VRMS误差Evms的方法为，
[0024] A1)计算梯度线圈在成像区域内的z方向磁场民的谐波系数A ^和B +
[0025] An,，Bnini 由下式确定，
[0026]
[0027] 其中，沒,炉)为成像区域内任意一点P处的球坐标，r为球坐标半径，9为P点 与z轴的夹角，供为P点在xy平面内的投影与x轴的夹角，Pnim(cos0 )为n阶m次连带勒 让德多项式；
[0028] A2)对谐波系数进行修正；
[0029] 对于X线圈，令：
[0030] An= A n+G
[0031] 对于Y线圈，令：
[0032] Bn= B n+G
[0033] 对于Z线圈，令：
[0034] A10= A10-G
[0035] 其中，An、Bn、A1iQ分别为X、Y、Z梯度线圈对应的主谐波系数，G为理想梯度磁场 的斜率；
[0036]A3)计算出成像区域内的VRMS误差：
[0037]
[0038] 其中，
R为成像区域半径。
[0039] 梯度线圈产生的磁场B通过毕奥-萨伐尔定律求得：
[0040]
[0041] 其中，y。为真空磁导率。
[0042] 本发明所达到的有益效果：1、本发明能够确保在整个成像区域内梯度场的线性误 差（线性度）最小，而非仅仅在个别采样点处的误差最小；2、本发明中将计算体积分误差 的三维积分公式简化为用谐波系数表示的公式，能够快速精确计算成像区域内磁场的VRMS 误差。
【附图说明】
[0043] 图1为本发明的流程图。
[0044] 图2为X线圈主线圈结构分布图。
[0045] 图3为X线圈屏蔽线圈结构分布图。
[0046] 图4为不同成像区域内梯度场的归一化Vrms误差。
【具体实施方式】
[0047] 下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明 的技术方案，而不能以此来限制本发明的保护范围。
[0048] 本发明受江苏省自然科学基金青年基金项目（项目批准号：BK20130854)支持。
[0049] 如图1所示，一种基于最小VRMS误差的梯度线圈设计方法，包括以下步骤：
[0050] 步骤一，预先给定梯度线圈所在的骨架尺寸以及预设参数；。
[0051] 预设参数包括成像区域范围、梯度场强度和VRMS误差Evms;
[0052] VRMS误差Ev"s为，
[0053]
[0054] 其中，V为成像区域的体积，Bz为梯度线圈产生的z方向磁场，G为梯度场的斜率， v为梯度场方向的坐标。
[0055] 对于球形成像区域，计算VRMS误差艮_的方法为
[0056] A1)计算梯度线圈在成像区域内的z方向磁场Bz的谐波系数AniB1和BniB1;
[0057] 采用毕奥-萨伐尔定律求梯度线圈产生的磁场B ;
[0058]
[0059] 其中，y。为真空磁导率；
[0060] 成像区域内的z方向磁场民的A^和B由下式确定，
[00611
[0062] 其中，〇%沒,为成像区域内任意一点P处的球坐标，r为球坐标半径，9为P点 与z轴的夹角，@为P点在xy平面内的投影与x轴的夹角，Pnim(cos0 )为n阶m次连带勒 让德多项式；
[0063] A2)对谐波系数进行修正；
[0064] 对于X线圈，令：
[0065] An= A n+G
[0066] 对于Y线圈，令：
[0067] Bn= B n+G
[0068] 对于Z线圈，令：
[0069] A10= A10-G
[0070] 其中，An、Bn、A1iQ分别为X、Y、Z梯度线圈对应的主谐波系数，G为理想梯度磁场 的斜率；
[0071] A3)计算出成像区域内的VRMS误差：
[0072]
广，、八
[0073] 其中，
R为成像区域半径。
[0074] 步骤二，构造流函数的基函数仁。
[0075] 流函数的N项展开式为，
[0076]
[0077] 其中，3;为基函数系数，i为整数，i=1，2.? ?N。
[0078] 步骤三，根据J= ▽X (nf)求得电流密度J的表达式；
[0079] 其中，嫁为骨架的外法向矢量；
[0080] 步骤四，建立与电流密度J有关的泛函数〇 ;
[0081]
[0082] 其中，Wl，w2为权重系数，W为线圈储能。
[0083]步骤五，求解使泛函数取得极值的{aj，得到流函数的表达式；
[0084] {aj为由&1构成的列矢量。
[0085] 步骤六，根据流函数表达式求得梯度线圈的形状。
[0086]根据上述方法，设计的一款带有源屏蔽的X方向梯度线圈的例子。
[0087] 线圈的内径为330mm，外径为410mm，梯度场强为58uT/m/A，采用上述算法，设计得 到的梯度线圈主线圈长度为114m，屏蔽线圈长度为77m，梯度线圈的总电感为350uH。所设 计的梯度线圈的主线圈与屏蔽线圈结构如图2与图3所示。图4中给出了该梯度线圈在不 同成像区域内的归一化VRMS误差。这里的归一化VRMS误差为用成像区域内的最大磁场对 Ewnis 进行归一化，即：Ev"y(GR)。
[0088] 通过上述方法，梯度线圈的形状可通过求解泛函的极值得到，所设计的梯度线圈 能够使泛函取得最小值，从而保证在线圈储能一定的情况下，线圈的VRMS误差最小，确保 在整个成像区域内梯度场的线性误差（线性度）最小，而非仅仅在个别采样点处的误差最 小；同时上述方法中将计算体积分误差的三维积分公式简化为用谐波系数表示的公式，能 够快速精确计算成像区域内磁场的VRMS误差；上述方法适用于设计圆柱形梯度线圈、平板 梯度线圈、不规则形状的梯度线圈等现存的任何类型的梯度线圈。
[0089] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人 员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形 也应视为本发明的保护范围。
【主权项】
1. 一种基于最小VMS误差的梯度线圈设计方法，其特征在于：包括W下步骤， 步骤一，预先给定梯度线圈所在的骨架尺寸W及预设参数； 预设参数包括成像区域范围、梯度场强度和VMS误差E"ms; VMS误差E"ms为，其中，V为成像区域的体积，氏为梯度线圈产生的Z方向磁场，G为梯度场的斜率，V为 梯度场方向的坐标； 步骤二，构造流函数的基函数fi; 流函数的N项展开式为，其中，曰1为基函数系数，i为整数，i= 1，2. . .N; 步骤S，根据二VX求得电流密度J的表达式； 其中，巧为骨架的外法向矢量； 步骤四，建立与电流密度J有关的泛函数? ;其中，Wi,W2为权重系数，W为线圈储能； 步骤五，求解使泛函数取得极值的{aiK得到流函数的表达式； {aj为由曰1构成的列矢量； 步骤六，根据流函数表达式求得梯度线圈的形状。2. 根据权利要求1所述的一种基于最小VMS误差的梯度线圈设计方法，其特征在于： 对于球形成像区域，计算VMS误差E"ms的方法为， Al)计算梯度线圈在成像区域内的Z方向磁场Bz的谐波系数A。,"和B八。,"1和6。,。由下式确定，其中，(/'，成矜)为成像区域内任意一点P处的球坐标，r为球坐标半径，0为P点与Z 轴的夹角，矜为P点在xy平面内的投影与X轴的夹角，P。,。(cos0 )为n阶m次连带勒让德 多项式； A2)对谐波系数进行修正； 对于X线圈，令： Aii=Aii+G 对于Y线圈，令： Bii=B11+G对于Z线圈，令： Aio=AIO-G 其中，All、Bii、Al,。分别为X、Y、Z梯度线圈对应的主谐波系数，G为理想梯度磁场的斜 率. A3)计算出成像区域内的VRMS误差：0, m =0 其中，，R为成像区域半径。 1, m# 03.根据权利要求2所述的一种基于最小VMS误差的梯度线圈设计方法，其特征在于： 梯度线圈产生的磁场B通过毕奥-萨伐尔定律求得：其中，y。为真空磁导率。
【专利摘要】本发明公开了一种基于最小VRMS误差的梯度线圈设计方法，其特征在于，梯度线圈的设计满足VRMS线性误差的平方与线圈储能的加权和构成的泛函最小。本发明同时给出了一种基于最小VRMS误差设计梯度线圈的算法流程。对于球形成像区域，本发明采用谐波系数展开来进行VRMS误差计算，具有很高的效率与精度。本发明适用于设计圆柱形梯度线圈、平板梯度线圈、不规则形状的梯度线圈等现存的任何类型的梯度线圈。
【IPC分类】G06F17/50
【公开号】CN105205227
【申请号】CN201510548202
【发明人】平学伟, 殷兴辉, 李黎, 李昌利, 刘海韵
【申请人】河海大学
【公开日】2015年12月30日
【申请日】2015年8月31日