CN110456293A - 一种自屏蔽梯度线圈的设计方法 - Google Patents
一种自屏蔽梯度线圈的设计方法 Download PDFInfo
- Publication number
- CN110456293A CN110456293A CN201910658436.2A CN201910658436A CN110456293A CN 110456293 A CN110456293 A CN 110456293A CN 201910658436 A CN201910658436 A CN 201910658436A CN 110456293 A CN110456293 A CN 110456293A
- Authority
- CN
- China
- Prior art keywords
- coil
- shielded
- imaging region
- main coil
- magnetic field
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
- G01R33/385—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/42—Screening
- G01R33/421—Screening of main or gradient magnetic field
- G01R33/4215—Screening of main or gradient magnetic field of the gradient magnetic field, e.g. using passive or active shielding of the gradient magnetic field
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Epidemiology (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种自屏蔽梯度线圈的设计方法,首先将根据成像区域目标点磁场值,通过边界元法和流函数法求出梯度线圈主线圈最大尺寸及主线圈在成像区域和屏蔽区域目标点的磁场值,根据主线圈在成像区域和屏蔽区域的磁场值来设定屏蔽线圈的目标点,从而控制屏蔽线圈在成像区域的磁场值,以实现对自屏蔽梯度线圈的功耗的控制。本发明采用上述结构的自屏蔽梯度线圈的设计方法,能有效提高自屏蔽梯度线圈的效率,降低自屏蔽梯度线圈的功耗。
Description
技术领域
本发明涉及磁共振成像领域,特别是涉及一种自屏蔽梯度线圈的设计方法。
背景技术
梯度线圈是磁共振成像设备的部件之一,其主要作用是在成像区域产生梯度磁场,实现对信号空间定位。在磁共振成像(MRI)中,要实现对信号的X、Y、Z三个方向空间定位,需要三路梯度线圈。梯度线圈在实现空间定位过程中,由于梯度电流的快速切断,会在梯度线圈附近的金属中产生涡流,涡流产生的磁场叠加到原梯度磁场中,会导致梯度线圈定位不准,产生涡流伪影和图像畸变。
目前克服梯度线圈涡流效应常用的方法有两种,一种是抗涡流极盘,另一种是自屏蔽梯度线圈。抗涡流极盘常采用双层硅钢片叠加,通过提高极盘电阻来降低梯度线圈产生的涡流,但硅钢片剩磁比较大,会使MRI成像系统中剩磁变大,会影响图像质量。第二种方法采用自屏蔽梯度线圈,通过使目标区域磁场变化为零,从而使极头上不会产生感应涡流。采用的设计方法主线圈和屏蔽线圈一体化设计方法,这种设计方法可以一次设计出主线圈和屏蔽线圈,但是不能控制屏蔽线圈在成像区域产生磁场的大小,在达到要求梯度场强条件下,一体化设计的线圈功耗会更大。
本发明专利为了解决上述自屏蔽梯度线圈的问题,提出了一种新的自屏蔽梯度线圈设计方法。既保证屏蔽区域磁场最小,又控制屏蔽线圈在成像区域产生梯度场最小,进一步降低自屏蔽梯度线圈的功耗。
发明内容
本发明的目的是提供一种自屏蔽梯度线圈的设计方法,能有效提高自屏蔽梯度线圈的效率且能降低自屏蔽梯度线圈的功耗。
为实现上述目的,本发明提供了一种自屏蔽梯度线圈的设计方法,其步骤如下:
步骤一:导入主梯度线圈网格,将主梯度线圈区域分为三维三角形网格节点。
步骤二:将成像区域直径360mm球体表面均匀划分156个目标点,采用MATLAB 软件读取三角网格的各个顶点和面的坐标值,优化顶点和面的排布顺序,确定成像区域目标点坐标值。
步骤三:根据成像区域目标点坐标值计算主线圈在成像区域理想梯度磁场值。
步骤四:根据边界元法,设置导线尺寸,计算源点区域通电导线对成像区域目标点的磁场贡献值。
步骤五:通过Quadprog二次规划方法,约束主线圈的功耗做小,计算出主线圈上电流的大小和方向。
步骤六:通过流函数法得到主线圈的实际绕线形状。
步骤七:根据主线圈的实际绕线形状,通过奥萨伐尔公式计算出成像区域目标点和屏蔽区域目标点磁场值。
步骤八:以屏蔽区域目标点磁场值作为屏蔽线圈在屏蔽区域的理想磁场值,以主线圈在理想区域设计值的1/5作为屏蔽线圈在成像区域的理想磁场值。
步骤九:导入屏蔽线圈的网格,排序网格顶点和三角面。
步骤十:根据边界元法和给定的导线尺寸,计算源点区域通电导线对成像区域目标点的磁场贡献值。
步骤十一:通过Quadprog二次规划方法,约束主线圈的功耗最小,计算出主线圈上电流的大小和方向。
步骤十二:通过流函数法得到主线圈的实际绕线形状。
本发明的有益效果为:
通过分步式设计自屏蔽梯度线圈的方法,在达到控制屏蔽区域磁场的同时,可有效控制屏蔽线圈在成像区域的梯度磁场,并降低整个自屏蔽梯度线圈的功耗。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1是本发明的自屏蔽梯度线圈的设计原理示意图。
图2是本发明的源区网格划分示意图。
图3是本发明的成像区域目标点示意图。
图4是本发明的三维显示源区和成像区域目标点示意图。
图5是本发明的纵向自屏蔽梯度线圈的z线圈结构示意图。
图6是本发明的纵向自屏蔽梯度线圈的x线圈结构示意图。
具体实施方式
下面结合附图对本发明的实施方式做进一步的说明。
参考图1所示,基于自屏蔽梯度线圈的原理结构,提出一种自屏蔽梯度线圈的设计方法,以下是本发明一种自屏蔽梯度线圈的设计方法的一个实施例:
所设计的主线圈间距为0.46m,直径0.405m;屏蔽线圈间距0.52m,直径0.435m,所设计的自屏蔽梯度线圈目标梯度磁场强度为15mT/m,线圈区域为非线性度不大于5%,具体设计步骤如下:
步骤一:根据所设定主线圈参数和目标点,采用blender2.7软件对主线圈进行三角网格划分,如图2所示,即将主线圈离散化。
步骤二:用matlab将球体分成13层,每隔30°设定一个测试点,一共 156个测试点,如图3所示,并求出这些测试点的x、y、z坐标。得到原坐标点S(x,y,z)和场坐标点F(x1,y1,z1)。
步骤三:根据目标场点坐标值确定目标点梯度磁场值。目标梯度磁场值为 15mT/m,设定主线圈目标磁场值为18.75mT/m,球面上目标点磁场值为坐标点z坐标值与梯度强度乘积,即
Gz=G*z
公式中Gz为给定目标区域内目标点的梯度磁场值,单位为mT;G为给定线性区域梯度强度,单位为mT/m;z为目标点z方向坐标值,单位为m。
步骤四:根据边界元法与设置的主线圈导线尺寸,计算源点区域通电导线对场点的贡献值,其计算方法为:源点区域离散化,离散成顶点和三角面,采用排序程序,对源点顶点和三角面进行排序,得出离散顶点的坐标值Ss(x,y,z)。
根据毕奥萨伐尔公式:
式中为源点导线对场点磁感应强度的贡献值;μ0为真空磁导率;dl为源点区域通电导线的长度;r为源点到场点的距离;I为源点导线上电流值;θ为通电导线与源点和场点连线之间夹角。
步骤五:计算主线圈的功耗和采用matlab中quadprog函数来计算功耗最小主线圈节点上电流值。电流功耗表达式如下:
式中面S为离散单元面,包括n各节点,Izm和Izn分别是主线圈第m个和第n个节点上电流值,ρ是主线圈导体的电阻,dzr是主线圈导体的厚度。Rzmn是主线圈的电阻矩阵。
quadprog函数如下:
Fz=(1/2)*Iz*Rzmn*Iz’
式中,F为主线圈目标函数,I为主线圈源区离散点电流值,I’为I的转置矩阵,Rzmn为主线圈源区离散点电阻矩阵。
步骤六:通过流函数法计算主线圈的绕线形状。
式中,Sz为主线圈等势差,max(Iz)为主线圈节点电流最大值,min(Iz)为主线圈节点电流最小值,Nz主线圈绕线匝数。
步骤七:通过主线圈的绕线形状及绕线的坐标点,采用比奥萨伐尔公式分别计算成像区域和屏蔽区域的磁场值。
步骤八:以目标场值的1/5倍作为屏蔽线圈在成像区域的目标磁场值,以主线圈在屏蔽区域产生的磁场的负数作为屏蔽线圈在屏蔽区域的目标磁场值。即屏蔽线圈的目标磁场值为
Bpt=[Bpimt,Bpt]
式中,Bpt为屏蔽线圈的目标磁场值,Bpimt为屏蔽线圈在成像区域的目标磁场值,Bpt为屏蔽线圈在屏蔽区域的目标磁场值。
步骤九:屏蔽线圈源区离散化,得出相应网格节点坐标并排序。
步骤十:计算屏蔽线圈源点对成像区域的贡献值,屏蔽线圈计算方法和主线圈计算公式相同。
屏蔽线圈的电阻矩阵计算公式如下:
式中面S为离散单元面,包括n各节点,Ipm和Ipn分别是屏蔽下圈第m个和第n个节点上电流值,ρ是屏蔽线圈导体的电阻,dpr是屏蔽线圈导体的厚度。Rpmn是屏蔽线圈的电阻矩阵。
步骤十一:采用matlab中quadprog函数来计算功耗最小屏蔽线圈节点上电流值。
Fp=(1/2)*Ip*Rpmn*Ip’
式中,F为屏蔽线圈目标函数,I为屏蔽线圈源区离散点电流值,I’为I的转置矩阵,Rpmn为源区离散点电阻矩阵。
步骤十二:通过流函数法计算屏蔽线圈的绕线形状,结果如图5、图6所示。
式中,Sp为屏蔽线圈等势差,max(Ip)为屏蔽线圈节点电流最大值,min(Ip)为屏蔽线圈节点电流最小值,Np屏蔽线圈绕线匝数。
因此,本发明提供了一种自屏蔽梯度线圈的设计方法,把自屏蔽梯度线圈分解成两种单独的梯度线圈设计,先设计出主梯度线圈,屏蔽线圈以主线圈在屏蔽区域产生的磁场值为目标点来设计屏蔽线圈,且控制屏蔽线圈在成像区域磁场值的大小,在达到设定梯度强度的条件下有效提高自屏蔽梯度线圈的效率且能降低自屏蔽梯度线圈的功耗。
最后应说明的是:以上实施例仅用以说明本发明的技术方案而非对其进行限制,尽管参照较佳实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对本发明的技术方案进行修改或者等同替换,而这些修改或者等同替换亦不能使修改后的技术方案脱离本发明技术方案的精神和范围。
Claims (1)
1.一种自屏蔽梯度线圈的设计方法,其特征在于:步骤如下:
步骤一:导入主梯度线圈网格,将主梯度线圈区域分为三维三角形网格节点;
步骤二:将成像区域直径360mm球体表面均匀划分156个目标点,采用MATLAB软件读取三角网格的各个顶点和面的坐标值,优化顶点和面的排布顺序,确定成像区域目标点坐标值;
步骤三:根据成像区域目标点坐标值计算主线圈在成像区域理想梯度磁场值;
步骤四:根据边界元法,设置导线尺寸,计算源点区域通电导线对成像区域目标点的磁场贡献值;
步骤五:通过Quadprog二次规划方法,约束主线圈的功耗做小,计算出主线圈上电流的大小和方向;
步骤六:通过流函数法得到主线圈的实际绕线形状;
步骤七:根据主线圈的实际绕线形状,通过奥萨伐尔公式计算出成像区域目标点和屏蔽区域目标点磁场值;
步骤八:以屏蔽区域目标点磁场值作为屏蔽线圈在屏蔽区域的理想磁场值,以主线圈在理想区域设计值的1/5作为屏蔽线圈在成像区域的理想磁场值;
步骤九:导入屏蔽线圈的网格,排序网格顶点和三角面;
步骤十:根据边界元法和给定的导线尺寸,计算源点区域通电导线对成像区域目标点的磁场贡献值;
步骤十一:通过Quadprog二次规划方法,约束主线圈的功耗最小,计算出主线圈上电流的大小和方向;
步骤十二:通过流函数法得到主线圈的实际绕线形状。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910658436.2A CN110456293B (zh) | 2019-07-22 | 2019-07-22 | 一种自屏蔽梯度线圈的设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910658436.2A CN110456293B (zh) | 2019-07-22 | 2019-07-22 | 一种自屏蔽梯度线圈的设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110456293A true CN110456293A (zh) | 2019-11-15 |
CN110456293B CN110456293B (zh) | 2021-07-20 |
Family
ID=68481508
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910658436.2A Active CN110456293B (zh) | 2019-07-22 | 2019-07-22 | 一种自屏蔽梯度线圈的设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110456293B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113030811A (zh) * | 2021-05-19 | 2021-06-25 | 中国科学院精密测量科学与技术创新研究院 | 一种圆筒型匀场线圈的设计方法 |
CN115047388A (zh) * | 2021-03-09 | 2022-09-13 | 宁波健信核磁技术有限公司 | 一种磁共振成像梯度线圈的制造组装方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0933645A2 (en) * | 1998-01-05 | 1999-08-04 | Picker International, Inc. | Design of a shielded gradient coil assembly for magnetic resonance imaging usingfinite element analysis |
US6049207A (en) * | 1998-11-25 | 2000-04-11 | Picker International, Inc. | Double-duty gradient coil assembly having two primary gradient coil sets and a common screening coil set |
CN101191829A (zh) * | 2006-11-29 | 2008-06-04 | 北京万东医疗装备股份有限公司 | 平面有源屏蔽梯度线圈的制作方法 |
CN101852843A (zh) * | 2010-05-27 | 2010-10-06 | 南京丰盛超导技术有限公司 | 一种超导磁体外磁屏蔽线圈的优化设计算法 |
CN101996273A (zh) * | 2010-11-29 | 2011-03-30 | 浙江大学 | Mri系统梯度线圈的有限差分设计方法 |
CN204515110U (zh) * | 2015-04-10 | 2015-07-29 | 惠仁望都医疗设备科技有限公司 | 一种用于核磁共振成像的新型自屏蔽纵向梯度线圈 |
CN104849681A (zh) * | 2015-06-10 | 2015-08-19 | 武汉中科波谱技术有限公司 | 一种核磁共振波谱仪梯度线圈 |
CN105445683A (zh) * | 2015-11-16 | 2016-03-30 | 河海大学 | 一种柱面横向自屏蔽梯度线圈设计方法 |
CN105718677A (zh) * | 2016-01-22 | 2016-06-29 | 中国科学院电工研究所 | 自屏蔽超导核磁共振成像系统梯度线圈设计方法 |
CN105799951A (zh) * | 2014-12-31 | 2016-07-27 | 上海新跃仪表厂 | 机电一体化微磁力矩器及磁矩测量方法 |
CN107831461A (zh) * | 2017-10-24 | 2018-03-23 | 中国计量大学 | 基于0‑1整数规划的纵向梯度线圈设计方法 |
CN107957565A (zh) * | 2017-12-21 | 2018-04-24 | 武汉中科牛津波谱技术有限公司 | 一种核磁共振波谱仪自屏蔽梯度线圈及其设计方法 |
CN108872896A (zh) * | 2018-05-29 | 2018-11-23 | 河北惠仁医疗设备科技有限公司 | 一种双平面磁共振成像系统梯度线圈的设计方法 |
-
2019
- 2019-07-22 CN CN201910658436.2A patent/CN110456293B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0933645A2 (en) * | 1998-01-05 | 1999-08-04 | Picker International, Inc. | Design of a shielded gradient coil assembly for magnetic resonance imaging usingfinite element analysis |
US6049207A (en) * | 1998-11-25 | 2000-04-11 | Picker International, Inc. | Double-duty gradient coil assembly having two primary gradient coil sets and a common screening coil set |
CN101191829A (zh) * | 2006-11-29 | 2008-06-04 | 北京万东医疗装备股份有限公司 | 平面有源屏蔽梯度线圈的制作方法 |
CN101852843A (zh) * | 2010-05-27 | 2010-10-06 | 南京丰盛超导技术有限公司 | 一种超导磁体外磁屏蔽线圈的优化设计算法 |
CN101996273A (zh) * | 2010-11-29 | 2011-03-30 | 浙江大学 | Mri系统梯度线圈的有限差分设计方法 |
CN105799951A (zh) * | 2014-12-31 | 2016-07-27 | 上海新跃仪表厂 | 机电一体化微磁力矩器及磁矩测量方法 |
CN204515110U (zh) * | 2015-04-10 | 2015-07-29 | 惠仁望都医疗设备科技有限公司 | 一种用于核磁共振成像的新型自屏蔽纵向梯度线圈 |
CN104849681A (zh) * | 2015-06-10 | 2015-08-19 | 武汉中科波谱技术有限公司 | 一种核磁共振波谱仪梯度线圈 |
CN105445683A (zh) * | 2015-11-16 | 2016-03-30 | 河海大学 | 一种柱面横向自屏蔽梯度线圈设计方法 |
CN105718677A (zh) * | 2016-01-22 | 2016-06-29 | 中国科学院电工研究所 | 自屏蔽超导核磁共振成像系统梯度线圈设计方法 |
CN107831461A (zh) * | 2017-10-24 | 2018-03-23 | 中国计量大学 | 基于0‑1整数规划的纵向梯度线圈设计方法 |
CN107957565A (zh) * | 2017-12-21 | 2018-04-24 | 武汉中科牛津波谱技术有限公司 | 一种核磁共振波谱仪自屏蔽梯度线圈及其设计方法 |
CN108872896A (zh) * | 2018-05-29 | 2018-11-23 | 河北惠仁医疗设备科技有限公司 | 一种双平面磁共振成像系统梯度线圈的设计方法 |
Non-Patent Citations (4)
Title |
---|
FENG JIA等: "Design multiple-layer gradient coils using least-squares finite element method", 《INDUSTRIAL APPLICATION》 * |
YASUSHI KANAI等: "Micromagnetic Model Analysis of Planar-Type Write Head Field Response and Dependence on Pole Tip, Return Yoke, and Shield Structure", 《IEEE TRANSACTIONS ON MAGNETICS,》 * |
徐文龙等: "基于目标场法的磁共振成像小尺寸射频线圈设计", 《高电压技术》 * |
轩倩倩: "基于目标场法的永磁微型核磁共振成像系统梯度线圈设计的多参", 《中国医学物理学杂志》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115047388A (zh) * | 2021-03-09 | 2022-09-13 | 宁波健信核磁技术有限公司 | 一种磁共振成像梯度线圈的制造组装方法 |
CN115047388B (zh) * | 2021-03-09 | 2023-10-17 | 宁波健信超导科技股份有限公司 | 一种磁共振成像梯度线圈的制造组装方法 |
CN113030811A (zh) * | 2021-05-19 | 2021-06-25 | 中国科学院精密测量科学与技术创新研究院 | 一种圆筒型匀场线圈的设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110456293B (zh) | 2021-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110456293A (zh) | 一种自屏蔽梯度线圈的设计方法 | |
JP2019515710A (ja) | 磁場シミングのための方法および装置 | |
CN102665543B (zh) | 倾斜磁场线圈、核磁共振成像装置以及线圈图形的设计方法 | |
CN106990373B (zh) | 一种磁共振系统的解耦轴向匀场线圈设计方法 | |
CN108872896A (zh) | 一种双平面磁共振成像系统梯度线圈的设计方法 | |
CN105718677A (zh) | 自屏蔽超导核磁共振成像系统梯度线圈设计方法 | |
CN110568390A (zh) | 一种双平面匀场线圈的设计方法 | |
CN109765510A (zh) | 一种带有圆角的径向超导匀场线圈及其设计方法 | |
CN109933911A (zh) | 密绕螺线管内金属圆筒电磁场有限元分析方法 | |
WO2013166810A1 (zh) | 磁共振成像超导磁体系统及其结构参数的获取方法与装置 | |
CN101852843B (zh) | 一种超导磁体外磁屏蔽线圈的优化设计方法 | |
CN102376410A (zh) | 一种开放式核磁共振超导磁体的设计方法 | |
WO2019221781A1 (en) | Resistive electromagnet systems | |
CN112684392A (zh) | 一种非平面梯度线圈的设计方法 | |
Handler et al. | Design and construction of a gradient coil for high resolution marmoset imaging | |
CN108107390A (zh) | 一种超导磁体外磁屏蔽线圈的优化设计方法 | |
Jeong et al. | Level-set-based topology optimization using remeshing techniques for magnetic actuator design | |
CN111060858B (zh) | 一种磁屏蔽桶内高均匀磁场与梯度复合磁场产生方法 | |
CN104198968B (zh) | 一种横向梯度线圈及其制作方法 | |
JP6259303B2 (ja) | 磁石装置および磁気共鳴撮像装置 | |
CN113971349A (zh) | 一种螺线管形超导磁体线圈数目和初始位置的获取方法 | |
US7633292B2 (en) | Method of reducing eddy currents caused by a gradient magnetic field in a magnetic resonance system | |
CN110298090B (zh) | 一种计及非对称性的螺线管线圈电磁场数值计算方法 | |
JPH0614900A (ja) | 勾配磁場コイルの製造方法及び勾配磁場用コイル単位並びに勾配磁場コイル | |
CN114264993A (zh) | 一种锥面自屏蔽梯度线圈的设计方法 |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20191115 Assignee: Hebei Wanshichun Railway Locomotive and Rolling Stock Parts Manufacturing Co.,Ltd. Assignor: HUIREN WANGDU MEDICAL EQUIPMENT SCIENCE & TECHNOLOGY CO.,LTD. Contract record no.: X2023980041858 Denomination of invention: A Design Method for Self Shielding Gradient Coil Granted publication date: 20210720 License type: Common License Record date: 20230914 |