CN110568390A - design method of biplane shimming coil - Google Patents
design method of biplane shimming coil Download PDFInfo
- Publication number
- CN110568390A CN110568390A CN201910975547.6A CN201910975547A CN110568390A CN 110568390 A CN110568390 A CN 110568390A CN 201910975547 A CN201910975547 A CN 201910975547A CN 110568390 A CN110568390 A CN 110568390A
- Authority
- CN
- China
- Prior art keywords
- biplane
- magnetic field
- coil
- point
- target
- 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.)
- Pending
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/387—Compensation of inhomogeneities
Abstract
The invention discloses a design method of a biplane shim coil, which comprises the steps of firstly introducing a source point grid of the biplane shim coil, ordering grid vertexes and triangular surfaces, then respectively calculating magnetic field values of 2-order shim coils on target points by adopting a spherical harmonic function, reading each vertex and each surface of the triangular grid by adopting MATLAB software, optimizing the arrangement sequence of the vertexes and the surfaces, setting the size of a lead according to a boundary element method, calculating the contribution value of a power-on lead of a source point region to a field point, restricting the energy storage and the power consumption of the biplane shim coil to be minimum, solving the current density distribution on the upper plane shim coil, and obtaining the winding shape of the biplane shim coil by a current function method. The invention adopts the design method of the biplane shimming coil with the structure, and can eliminate 2-order and above harmonic magnetic field, thereby improving the uniformity of the main magnetic field and improving the image quality of the open MRI system.
Description
Technical Field
The invention relates to the field of magnetic resonance imaging, in particular to a design method of a biplane shim coil.
Background
magnetic Resonance imaging (magnetic Resonance imaging MRI) technology is a multi-technology crossed technology, including the subjects of electromagnetism, digital signal processing, biomedicine, atomic physics and the like, compared with other medical imaging technologies, MRI has the advantages of no radiation, high resolution, clear image quality and the like, so the MRI technology has important clinical application value in medical examination, ideally, the main magnetic field of MRI equipment only has a B0 field and no other harmonic magnetic field, however, due to the existence of processing and assembling errors, a harmonic magnetic field other than a B0 field is inevitably introduced, the introduction of the harmonic magnetic field can cause poor imaging quality or even cannot image, a gradient coil can correct a 1-order harmonic magnetic field, and a magnetic field higher than the 1-order harmonic magnetic field needs a 2-order shimming coil and more than the 2-order shimming coil to eliminate, so that a B0 field reaches the range of imaging requirements.
the invention provides a design method of a biplane shim coil in order to solve the problem of main magnetic field nonuniformity of an open MRI magnet. The magnetic field of harmonic wave of 2 orders and above can be eliminated, the uniformity of the main magnetic field is improved, and the image quality of the open type MRI system is improved.
Disclosure of Invention
the invention aims to provide a design method of a biplane shim coil, which can eliminate 2-order and above harmonic magnetic field, improve the uniformity of a main magnetic field and improve the image quality of an open MRI system.
In order to achieve the purpose, the invention provides a method for designing a biplane shim coil, which comprises the following steps.
The method comprises the following steps: and (4) introducing a biplane shimming coil grid, and dividing the main gradient coil area into three-dimensional triangular grid nodes.
Step two: dividing spheres with set diameters into 16 layers, setting a test point every 11.6 degrees, and obtaining 1488 target points by using MATLAB to obtain x, y and z coordinates of the target points to obtain field coordinate points.
Step three: and determining a target point magnetic field value according to the target point coordinate value and the 2 nd order spherical harmonic coefficient, wherein the target point magnetic field value on the spherical surface is the product of the coordinate point z coordinate value and the spherical harmonic coefficient.
Step four: according to the boundary element method and the set dimensions of the biplane shim coil wires, calculating the contribution value of the electrified wire in the source point region to the target field point: discretizing the source point region into a vertex and a triangular surface, sequencing the source point vertex and the triangular surface by adopting an MATLAB sequencing program to obtain coordinate values of the discrete vertex, and calculating the contribution value of the electrified conducting wire of the source point region to the target field point according to a Biao savart formula.
step five: and calculating the power consumption of the biplane shim coil, constraining the biplane shim coil to have minimum power consumption and energy storage through a quadrprog quadratic programming method, and calculating the current value on the minimum shim coil node.
Step six: and obtaining the actual winding shape of the biplane shim coil by a flow function method.
Step seven: and solving a magnetic field value of the coil on a target point by utilizing the Bio savart theorem according to the winding shape of the shimming coil, judging whether the magnetic field value meets the error requirement of the target magnetic field value, and if not, modifying the weight coefficients of the power consumption matrix and the energy storage matrix until the magnetic field value meets the error requirement of the target magnetic field value.
The invention has the beneficial effects that:
The biplane shimming coil designed by the invention can effectively eliminate 2-order and above harmonic magnetic field introduced by machining and assembling errors, and simultaneously improves the uniformity of a main magnetic field and the image quality.
The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.
drawings
Fig. 1 is a flow chart of a design method of the bi-planar shim coil of the present invention.
Fig. 2 is a schematic diagram of source region meshing of the shim coil of the present invention.
FIG. 3 is a schematic view of the imaging region target point of the present invention.
FIG. 4 is a schematic diagram of a three-dimensional display source region and imaging region target point of the present invention.
fig. 5 is a schematic diagram of the structure of a biplane Z20 shim coil of the present invention.
Fig. 6 is a schematic diagram of the structure of a biplane Z21 shim coil of the present invention.
Fig. 7 is a schematic diagram of the structure of a biplane Z22 shim coil of the present invention.
Detailed Description
Embodiments of the present invention will be further described with reference to the accompanying drawings.
the invention uses 3 groups of 2-order axial shimming coils: the design method of the biplane shim coil is introduced by Z20, Z21 and Z22 as design examples, the design flow chart is shown in figure 1, the design parameters of the biplane shim coil are set, the biplane shim coil is a triangular grid with the diameter of 0.42m and the distance is 500mm, each sphere on 100mm, 240mm and 360mm spheres is divided into 16 layers, a test point is set at intervals of 11.6 degrees, 1488 target points are set, the error is not more than 5%, and the specific design steps are as follows.
The method comprises the following steps: according to the design parameters of the set biplane shim coil, performing triangularization grid division on the shim coil by using blend 2.8 software, namely dispersing the shim coil into nodes and triangular surfaces, as shown in fig. 2. And exporting the grid nodes into an obj format, reading and sequencing the obj format file through MATLAB to obtain source point coordinate points S (x, y, z).
Step two: spheres with diameters of 100mm, 240mm and 360mm are divided, each sphere is divided into 16 layers, and test points are set every 11.6 degrees for 1488 target points, as shown in fig. 3 and 4. The x, y, z coordinates of these target points were found using MATLAB to obtain field coordinate points F (x1, y1, z 1).
Step three: determining a target point magnetic field value according to the target field point coordinate value and the 2 nd order spherical harmonic coefficient, wherein the target spherical harmonic coefficient is 5mT/m, and the target point magnetic field value on the spherical surface is the product of the coordinate point z coordinate value and the spherical harmonic coefficient, namely:
Qz=Q*z
In the formula, Qz is the magnetic field value of a target point in a given target area, and the unit is mT; q is a spherical harmonic coefficient of a given order, and the unit is mT/m; z is the coordinate value of the target point in the z direction and the unit is m.
Step four: and calculating the contribution value of the electrified lead of the source point region to the target field point according to a boundary element method and the set size of the biplane shim coil lead. Discretizing the source point region into a vertex and a triangular surface, and sequencing the vertex and the triangular surface of the source point by adopting an MATLAB sequencing program to obtain coordinate values S (x, y, z) of the discrete vertex.
according to the biot savart formula:
In the formulaThe contribution value of the source point lead to the field point magnetic induction intensity is obtained; mu.s0Is a vacuum magnetic conductivity; dl is the length of the electrified lead in the source region; r is the distance from the source point to the field point; i is the current value on the source point lead; theta is an included angle between the electrified lead and a connecting line of the source point and the field point.
Step five: calculating the power consumption of the biplane shim coil, calculating the current value on the node of the minimum shim coil with the power consumption and the energy storage by adopting a quadprog function in matlab,
the coil power consumption expression is as follows:
Wherein the surface S is a discrete unit surface comprising n nodes, ImAnd InCurrent values at the m-th and n-th nodes, respectively, p being the resistance of the conductor, drIs the thickness of the conductor. RmnIs a resistance matrix of the gradient coil;
The coil energy storage expression is as follows:
Formula middle surface Smand SnAre discrete triangular surfaces and belong to nodes n and m respectively, and the nodes n and m respectively contain WnAnd WmA triangle. I ismand InCurrent values, μ, at the mth and nth nodes, respectively0Is the magnetic permeability of vacuum, rmAnd rnRespectively, the coordinates of the points in the triangular plane. v. ofmaand vnbWhich are the basis functions of nodes m and n, respectively. MmnIs an energy storage matrix of the shim coils;
The quadprog function is as follows:
F=α*I*Rmn*I’+β*I*Mmn*I’
Wherein F is an objective function, I is a discrete point current value of the source region, I' is a transposed matrix of I, RmnIs a source region discrete point resistance matrix, Mmnthe method comprises the steps that a source region discrete point inductance matrix is adopted, and alpha and beta are respectively a resistance matrix and an inductance matrix weight coefficient;
The constraint conditions are as follows:
A*I<ε*Bt
In the formula, A is a magnetic field coefficient matrix which does not contain a node current value and is calculated by Bio savart at a target point through a source point; epsilon is the magnetic field error, epsilon is 0.05; b istIs the target magnetic field value; i is the node current value.
Step six: the winding shape of the biplane shim coil is calculated by a flow function method,
In the formula, S is an equipotential difference, max (I) is a maximum node current value, min (I) is a minimum node current value, and the number of turns of the N biplane shim coil windings is equal to the maximum node current value.
Step seven: and solving the magnetic field value of the coil on the target point by utilizing the Bio Saval theorem according to the winding shape of the shimming coil, and judging whether the magnetic field value meets the error requirement of the target magnetic field value. If the target magnetic field value requirement of the shimming coil is met, stopping modifying the alpha and beta weighting coefficients, otherwise, continuing the alpha and beta weighting coefficients until the linearity requirement of the gradient coil is met, and designing results are shown in fig. 5, 6 and 7.
Therefore, the invention provides a design method of the biplane shim coil, which controls the power consumption and energy storage minimization of the shim coil, restrains the magnetic field value of the shim coil on a target point, can effectively eliminate a high-order harmonic magnetic field induced by machining and assembling errors, and further improves the image quality of an open MRI system.
Finally, it should be noted that: the above embodiment is exemplified by a 2-order axial shim coil and is not limited to a 2-order shim coil, a higher-order shim coil can be designed by modifying the order in the spherical harmonic coefficient, and the above embodiment is only used for illustrating the technical solution of the present invention and is not limited thereto, although the present invention is described in detail with reference to the preferred embodiment, those skilled in the art should understand that: modifications and equivalents may be made to the invention without departing from the spirit and scope of the invention.
Claims (1)
1. A design method of a biplane shim coil is characterized in that: the method comprises the following steps:
The method comprises the following steps: leading in a biplane shimming coil grid, and dividing a main gradient coil area into three-dimensional triangular grid nodes;
step two: dividing spheres with set diameters into 16 layers, setting a test point every 11.6, and obtaining 1488 target points, and calculating x, y and z coordinates of the target points by using MATLAB to obtain a field coordinate point F (x1, y1 and z 1);
Step three: determining a target point magnetic field value according to the target point coordinate value and the 2 nd order spherical harmonic coefficient, wherein the target point magnetic field value on the spherical surface is the product of the coordinate point z coordinate value and the spherical harmonic coefficient;
step four: according to the boundary element method and the set dimensions of the biplane shim coil wires, calculating the contribution value of the electrified wire in the source point region to the target field point: discretizing the source point region into a vertex and a triangular surface, sequencing the source point vertex and the triangular surface by adopting an MATLAB sequencing program to obtain a coordinate value of the discrete vertex, and calculating a contribution value of a power-on lead of the source point region to a target field point according to a Biao savart formula;
Step five: calculating the power consumption of the biplane shimming coil, constraining the biplane shimming coil to have minimum power consumption and energy storage through a quadrprog quadratic programming method, and calculating the current value on the minimum shimming coil node;
Step six: obtaining the actual winding shape of the biplane shim coil by a flow function method;
step seven: and solving a magnetic field value of the coil on a target point by utilizing the Bio savart theorem according to the winding shape of the shimming coil, judging whether the magnetic field value meets the error requirement of the target magnetic field value, and if not, modifying the weight coefficients of the power consumption matrix and the energy storage matrix until the magnetic field value meets the error requirement of the target magnetic field value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910975547.6A CN110568390A (en) | 2019-10-14 | 2019-10-14 | design method of biplane shimming coil |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910975547.6A CN110568390A (en) | 2019-10-14 | 2019-10-14 | design method of biplane shimming coil |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110568390A true CN110568390A (en) | 2019-12-13 |
Family
ID=68785146
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910975547.6A Pending CN110568390A (en) | 2019-10-14 | 2019-10-14 | design method of biplane shimming coil |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110568390A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112444766A (en) * | 2020-11-05 | 2021-03-05 | 上海联影医疗科技股份有限公司 | Magnetic resonance system and shimming method thereof |
CN113030811A (en) * | 2021-05-19 | 2021-06-25 | 中国科学院精密测量科学与技术创新研究院 | Design method of cylindrical shimming coil |
CN113325350A (en) * | 2020-02-28 | 2021-08-31 | 河海大学 | High-performance gradient coil design method based on discrete grids |
CN113671431A (en) * | 2021-08-20 | 2021-11-19 | 宁波健信核磁技术有限公司 | Decoupling method of magnetic resonance high-order shimming coil and related device |
CN114002634A (en) * | 2021-11-15 | 2022-02-01 | 安徽工程大学 | Calibration device and method for magnetic field-voltage coefficient of multi-channel magnetocardiogram detection system |
CN114217254A (en) * | 2021-12-13 | 2022-03-22 | 北京航空航天大学 | High-linearity gradient coil design method |
CN114355263A (en) * | 2020-09-29 | 2022-04-15 | 河海大学 | Design method of high-order shimming coil |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080116894A1 (en) * | 2006-04-10 | 2008-05-22 | Markus Weiger | Method for automatic shimming for nuclear magnetic resonance spectroscopy |
CN204515110U (en) * | 2015-04-10 | 2015-07-29 | 惠仁望都医疗设备科技有限公司 | The longitudinal gradient coil of a kind of novel self-shileding for Magnetic resonance imaging |
CN104865544A (en) * | 2015-06-10 | 2015-08-26 | 武汉中科波谱技术有限公司 | Room temperature shim coil of superconducting NMR spectrometer |
CN105548925A (en) * | 2015-12-18 | 2016-05-04 | 中国科学院苏州生物医学工程技术研究所 | Holbach magnet shim coil and design method thereof |
CN107530026A (en) * | 2015-05-12 | 2018-01-02 | 海珀菲纳研究股份有限公司 | Radio-frequency coil method and apparatus |
CN108872896A (en) * | 2018-05-29 | 2018-11-23 | 河北惠仁医疗设备科技有限公司 | A kind of design method of biplane magnetic resonance image-forming system gradient coil |
-
2019
- 2019-10-14 CN CN201910975547.6A patent/CN110568390A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080116894A1 (en) * | 2006-04-10 | 2008-05-22 | Markus Weiger | Method for automatic shimming for nuclear magnetic resonance spectroscopy |
CN204515110U (en) * | 2015-04-10 | 2015-07-29 | 惠仁望都医疗设备科技有限公司 | The longitudinal gradient coil of a kind of novel self-shileding for Magnetic resonance imaging |
CN107530026A (en) * | 2015-05-12 | 2018-01-02 | 海珀菲纳研究股份有限公司 | Radio-frequency coil method and apparatus |
CN104865544A (en) * | 2015-06-10 | 2015-08-26 | 武汉中科波谱技术有限公司 | Room temperature shim coil of superconducting NMR spectrometer |
CN105548925A (en) * | 2015-12-18 | 2016-05-04 | 中国科学院苏州生物医学工程技术研究所 | Holbach magnet shim coil and design method thereof |
CN108872896A (en) * | 2018-05-29 | 2018-11-23 | 河北惠仁医疗设备科技有限公司 | A kind of design method of biplane magnetic resonance image-forming system gradient coil |
Non-Patent Citations (2)
Title |
---|
戚金凤: "Halbach核磁共振成像系统的梯度线圈设计方法", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
熊丙权: "超低场磁共振双平面有源匀场线圈设计", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113325350A (en) * | 2020-02-28 | 2021-08-31 | 河海大学 | High-performance gradient coil design method based on discrete grids |
CN114355263A (en) * | 2020-09-29 | 2022-04-15 | 河海大学 | Design method of high-order shimming coil |
CN114355263B (en) * | 2020-09-29 | 2023-10-24 | 河海大学 | Design method of high-order shimming coil |
CN112444766A (en) * | 2020-11-05 | 2021-03-05 | 上海联影医疗科技股份有限公司 | Magnetic resonance system and shimming method thereof |
CN112444766B (en) * | 2020-11-05 | 2023-09-26 | 上海联影医疗科技股份有限公司 | Magnetic resonance system and shimming method thereof |
CN113030811A (en) * | 2021-05-19 | 2021-06-25 | 中国科学院精密测量科学与技术创新研究院 | Design method of cylindrical shimming coil |
CN113030811B (en) * | 2021-05-19 | 2021-09-07 | 中国科学院精密测量科学与技术创新研究院 | Design method of cylindrical shimming coil |
CN113671431A (en) * | 2021-08-20 | 2021-11-19 | 宁波健信核磁技术有限公司 | Decoupling method of magnetic resonance high-order shimming coil and related device |
CN114002634A (en) * | 2021-11-15 | 2022-02-01 | 安徽工程大学 | Calibration device and method for magnetic field-voltage coefficient of multi-channel magnetocardiogram detection system |
CN114217254A (en) * | 2021-12-13 | 2022-03-22 | 北京航空航天大学 | High-linearity gradient coil design method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110568390A (en) | design method of biplane shimming coil | |
JP5497785B2 (en) | Gradient magnetic field coil and nuclear magnetic resonance imaging apparatus | |
US4646024A (en) | Transverse gradient field coils for nuclear magnetic resonance imaging | |
CA2728108C (en) | Magnetic assembly and method for defining a magnetic field for an imaging volume | |
CN106556813B (en) | Linear mixed optimization method of active shimming coil in magnetic resonance system | |
US4840700A (en) | Current streamline method for coil construction | |
CN102456460B (en) | There is the superconducting magnet of cold iron shimming capability | |
CN109765510B (en) | Method for designing radial superconducting shimming coil with fillet | |
WO2013166810A1 (en) | Magnetic resonance imaging superconducting magnet system and method and device for acquiring structural parameter thereof | |
EP0140259B1 (en) | Transverse gradient field coil | |
CN105445683B (en) | A kind of cylinder transverse direction self-shielded gradient coils design method | |
CN112684392A (en) | Design method of non-planar gradient coil | |
JP2009502031A (en) | Multilayer magnet | |
CN109856575B (en) | A kind of method, apparatus, equipment and the storage medium of the passive shimming of determining magnetic resonance | |
Shan et al. | Geometric distortion characterization and correction for the 1.0 T Australian MRI‐linac system using an inverse electromagnetic method | |
CN112444766A (en) | Magnetic resonance system and shimming method thereof | |
Handler et al. | Design and construction of a gradient coil for high resolution marmoset imaging | |
CN104833930B (en) | The computational methods of open type magnetic resonance system gradient coil magnetic field intensity | |
CN110456293B (en) | Design method of self-shielding gradient coil | |
JP5481216B2 (en) | Magnetic resonance imaging system | |
CN107205690A (en) | Magnetic field adjustment method | |
Tashiro et al. | Simple-Box-9 coil system: A novel approach to design of a square coil system for producing uniform magnetic fields | |
JP4996523B2 (en) | Coil pattern calculation method and gradient coil | |
CN110703170A (en) | Design method of special breast MRI gradient coil | |
CN115032579A (en) | Design method of non-equal-height self-shielding gradient coil for superconducting magnet |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20191213 |
|
WD01 | Invention patent application deemed withdrawn after publication |