The present invention relates to software for adjusting magnetic homogeneity, a method for adjusting magnetic homogeneity, a magnetic device, and a magnetic resonance imaging apparatus.
Using a nuclear magnetic resonance phenomenon that occurs when a specimen placed in a homogeneous static magnetic field is irradiated with high frequency pulses, a magnetic resonance imaging (MRI) apparatus can obtain an image representing the physical and chemical properties of the specimen, and is particularly used for medical purposes. The MRI apparatus mainly includes a magnetic field generation source for applying a homogeneous static magnetic field in an imaging region into which the specimen is carried, an RF coil for irradiating the imaging region with high frequency pulses, a receiving coil for receiving a response from the imaging region, and a gradient magnetic field coil for applying a gradient magnetic field to provide the imaging region with position information of a resonance phenomenon.
In the MRI apparatus, one of the factors for improving the image quality is an improvement in the static magnetic field homogeneity in the imaging region. In designing and manufacturing of a magnetic device used for the MRI apparatus, magnetic field homogeneity adjustment is performed at the respective stages of designing, assembling, and installing in order to make a static magnetic field, which is generated in the imaging region by a magnetic field generation source, homogeneous.
Among these, the magnetic field homogeneity adjustment performed at the installing stage can be realized by adding or removing magnetic field homogeneity adjusting pieces (magnetic shims) of a magnetic material to/from a magnetic device, for example, when a magnetic field inhomogeneity component has been caused by a manufacturing error or the surrounding environment. For example, in a magnetic device of a type that forms an imaging region and a homogeneous magnetic field space thereof between magnetic field generation sources (magnetic poles) vertically facing each other, a structure is formed, in general, by providing and disposing a magnetic field homogeneity adjustment mechanism (means), which is called a shim tray, in a tray shape of a non-magnetic material in each of the spaces sandwiched by the respective magnetic poles and respective gradient magnetic field coils disposed inside with respect to the magnetic poles (namely, on the imaging region side) (for example, refer to Patent Document 1).
On the other hand, in a magnetic device of a type that incorporates a plurality of superconducting coils to be a magnetic field generation source in a double cylindrical container and forms an imaging region inside thereof and a homogeneous magnetic filed space along the axial direction of the cylinder, a structure is formed, in general, by providing a shim tray (magnetic field homogeneity adjuster) in the space sandwiched by a gradient magnetic field coil disposed on the inner circumferential side of the container and the inner circumferential surface of the container, or by incorporating a shim tray in the gradient magnetic field coil (for example, refer to Patent Document 2).
Consideration on where and how many magnetic shims are to be disposed on these shim trays is, in general, an optimization problem having an objective function for the magnetic field homogeneity in the imaging region, and disposition of magnetic shims is often determined by a linear optimization method or a method modified therefrom, using a given magnetic field distribution (for example, refer to Patent Document 3).
As a magnetic field homogeneity adjuster, a structure that allows disposition of magnetic shims (with a large volume) as many as possible per unit area makes a high magnetic field homogeneity adjustment capacity. This is because a large change in the magnetic field strength can be generated in an imaging region by a large volume of magnetic shims. In order to easily realize this, a method was considered that makes a large number of screw holes through a shim tray and screws magnetic material shims called shim bolts into these screw holes. In this method, by finely disposing screw holes, it is possible to dispose shim bolts as many as the number of the screw holes, and as a result, a large number of magnetic shims can be disposed. Further, by disposing screw holes at fine relative positions to each other, formation of a spatially fine magnetic filed distribution can be expected.
However, in order to realize a homogeneous magnetic field, in a case of individually managing a large number of screw holes and thus applying an optimization method such as a linear programming method to determine disposition of shim bolts, an elaborate work of screwing optimum shim bolts into the respective holes without an error is required. Assuming that, for example, several thousands of screw holes have been formed through one shim tray, a person to be engaged in the magnetic field homogeneity adjustment work needs to accurately array screws (shim bolts) necessary for the respective screw holes, which makes the work efficiency extremely low.
In order to improve this work efficiency, for example, a method of setting the number to a required minimum and thus decreasing the number of screw holes to be managed was considered. In this situation, it was also considered to make the diameter of screw holes to be provided through a shim tray large, however, the size of a shim bolt, in other words, the size of a screw hole cannot be made quite large, taking into account handling shim bolts for magnetic field homogeneity adjustment in a strong magnetic field generated by a magnetic device. Consequently, there is a problem that a sufficient magnetic field homogeneity adjustment capacity cannot be obtained.
Therefore, a method was considered, as a method for improvement, that divides a shim tray into several regions in advance such that each region includes a plurality of screw holes through the shim tray without a change in the diameter nor the number of holes, adds the volumes of shim bolts to be disposed at the screw holes in each region, and then displays the individual total volumes in the respective regions together. Herein, the above-described regions can be designed to attain a spatial accuracy sufficient for magnetic field homogeneity adjustment. It is necessary to set the regions to have a size matching a magnetic field homogeneity adjustment that is the finest adjustment expected at the time of adjusting the installation of a magnetic device. In such a manner, the work efficiency can be made higher than that for a case of individually managing a large number of screw holes.
However, because the size of these regions is normalized merely to ensure a sufficient spatial accuracy of the magnetic field distribution as described above, the size inevitably becomes small, in other words, the number of regions still remains large. If such region dividing is performed, even in a case of performing a sort of general (in other words, not-detailed) magnetic field homogeneity adjustment that does not require a high spatial accuracy of magnetic field distribution, it is necessary to finely dispose shim bolts in many regions, which causes a problem of a low work efficiency.
An object of the present invention is to provide software for adjusting magnetic homogeneity, a method for adjusting magnetic homogeneity, a magnetic device, and an MRI apparatus that contribute to improvement in the work efficiency in performing general magnetic field homogeneity adjustment which does not require, as described above, a significantly fine spatial accuracy of the magnetic field distribution.
To solve the above-described problems, in accordance with the present invention, the positions and volumes of magnetic members (for example, shim bolts) to be disposed on a shim tray are calculated first on a computational mesh, based on the magnetic field strength distribution in a magnetic field space, to make the magnetic field strength distribution homogeneous. Subsequently, from the distribution of the calculated positions and volumes of the magnetic members, the local maximum values and the local minimum values thereof are extracted; the volume distribution regions of the magnetic members, with the centers thereof respectively at the positions of the extracted local maximum and minimum values, are extracted; and the volumes of the magnetic members distributed in these distribution regions are added in the respective regions. Finally, results of these calculations are displayed together with the corresponding local maximum value positions or the corresponding local minimum value positions.
Preferably, this method of extraction of distribution regions is desired to be a method that establishes the regions while checking the relationship with respect to the mass between adjacent nodes on the computational mesh and sequentially expanding the regions.
Further, the mass display method, described above, preferably displays the positions and the mass to be visually recognizable on a screen that displays the shape of the shim tray.
According to the present invention, it is possible to provide software for adjusting magnetic homogeneity, a method for adjusting magnetic homogeneity, a magnetic device, and an MRI apparatus that contribute to improvement in the work efficiency in performing general magnetic field homogeneity adjustment which does not require a significantly fine spatial accuracy of the magnetic field distribution.
Embodiments in accordance with the present invention will be described in detail, with reference to the accompanying drawings.
As shown in
As shown in
Returning to
Likewise, the lower coil container 2 is formed with a cylindrical recess 16 on the surface thereof facing the magnetic field space 3. A shim tray 18 of a non-magnetic material is housed in the recess 16. A gradient magnetic field coil 20 is disposed on the magnetic field space 3 side of the shim tray 18, and an RF transmitting/receiving coil 22 is disposed between the magnetic field space 3 and the gradient magnetic field coil 20.
The superconducting coils 8, 9, 10, and 11 form an imaging region 23, which is a part of the magnetic field space 3, as a homogeneous magnetic field space. The superconducting coils (primary coils) 8, 9 generate the strongest magnetic field and form a static magnetic field along the vertical direction in the magnetic filed space 3. The shield coils 10, 11 are provided to prevent the magnetic field formed by the superconducting coils (primary coils) 8, 9 from leaking outside. Further, the gradient magnetic field coils 19, 20 form a dynamic magnetic field in the imaging region 23. The RF transmitting/receiving coils 21, 22 irradiate the imaging region 23 with an electromagnetic wave (radio wave) and receive the electromagnetic wave.
The superconducting coils 8, 9, 10, and 11 are disposed such as to generate a homogeneous magnetic field in the imaging region 23, as described above. If the superconducting coils 8, 9, 10, and 11 are insufficient to obtain necessary strength or homogeneity of the magnetic field, then ferromagnetic members (not shown), such as iron pieces (including iron alloy, the same hereinafter) or permanent magnets, are disposed (or removed from), for example, inside or outside the vacuum container 12, inside the radiation shield 13, or inside the helium container 14 to increase (or attenuate) the magnetic field strength or to improve the homogeneity. The above description has been made for a case of disposing four superconducting coils 8, 9, 10, and 11, however, more or fewer superconducting coils may be disposed.
In such a manner, the magnetic device 50 is designed such as to generate a homogeneous magnetic field, using the superconducting coils 8, 9, 10, and 11 and iron pieces or the like (not shown), however, in reality, an error magnetic field is generated in the imaging region 23 by an assembling error, effects by the installation environment, or the like. The shim trays 17, 18 are provided in order to remove this error magnetic field component.
From the surface of the upper coil container 1, the shim tray 17, the gradient magnetic field coil 19, and the RF transmitting/receiving coil 21 are disposed in this order. Likewise, from the surface of the lower coil container 2, the shim tray 18, the gradient magnetic field coil 20, and the RF transmitting/receiving coil 22 are disposed in this order. The gradient magnetic field coils 19, 20, and the RF transmitting/receiving coils 21, 22 are installed to be removable. The shim trays 17, 18 may be or may not be removable.
The shim trays 17, 18 have a shape of a disk or the like and are provided with a number of screw holes (female screw) 26 therethrough. In magnetic field homogeneity adjustment, when shim bolts 27, which are magnetic shims in a screw shape (male screw), are screwed into the screw holes 26, the shim trays 17, 18 are added with a magnetic material by the shim bolts 27. Shim bolts 27 are prepared in advance, having various volumes and shapes depending on the length and the machining method, and a worker for magnetic field homogeneity adjustment selects and uses appropriate shim bolts 27 with required volumes and shapes.
As shown in
The magnetic field homogeneity adjustment work means disposing shim bolts 27, which are necessary to make the magnetic field distribution in the imaging region 23 homogeneous, in the screw holes 26 provided through the shim trays 17, 18. Based on the measurement values of the magnetic field strength distribution in the imaging region (homogeneous magnetic field space) 23, software (software for adjusting magnetic homogeneity), which is installed on a computer, calculates at which positions and in what approximate volumes shim bolts 27 are to be disposed on the shim trays 17, 18 in order to obtain a desirable homogeneous magnetic field.
The algorithm of the software determining the disposition may be based on, for example, a numerical programming method such as a linear programming method, other optimization methods, and may be based on a method that solves an inverse problem. In the present embodiment, an algorithm according to the inverse problem solution method will be described as an example.
First, the initial measurement is, as described later, performed in a state that shim bolts 27 are not disposed on the shim trays 17, 18, then measurement is repeated while shim bolts 27 are sequentially added, and thereby a predetermined magnetic field homogeneity is obtained. As shown in
When a shim bolt 27 with a volume Vi and a magnetic charge M is disposed at a certain node i on the computational mesh of the shim tray 17, 18, the shim bolt 27 causes a magnetic field strength B (i, j) at a certain adjacent node j in the imaging region (homogeneous magnetic field space) 23, the magnetic field strength B being proportional to the volume Vi and the magnetic charge M.
Expression 1
B(i, j)∝ViM=mi (1)
The symbol mi represents the magnetic dipole moment of the shim bolt 27. Herein, the magnetic charge M is assumed to be constant. Accordingly, the distribution of the magnetic moments of the shim bolts 27 disposed at the respective nodes on the computational mesh of the shim tray 17, 18 is expressed by the following expression.
The distribution of the magnetic field strengths created, by these, at the respective nodes on the computational mesh of the imaging region (the homogeneous magnetic field space) 23 is expressed by the following expression.
Then, the relationship between the magnetic field distribution and the magnetic moment distribution is expressed by the following expression, representing the coefficient matrix by matrix A.
Expression 4
{right arrow over (b)}=A{right arrow over (m)} (4)
Applying a singular value decomposition method to the matrix A, a generalized inverse matrix A′ of the matrix A can be obtained. As a result, the following expression is obtained. The singular value decomposition method is, for example, the method described in the above-described “Non-patent Document 1”.
Expression 5
{right arrow over (m)}=A′{right arrow over (b)} (5)
That is, once the magnetic field distribution (to be generated) as a target is determined, required magnetic moment distribution can be calculated by calculating the matrix product between itself and the generalized inverse matrix A′ as shown in Expression (5). As described in the present invention, for a magnetic field homogeneity adjustment work, namely, a work for making the magnetic field distribution of the imaging region (homogeneous magnetic field space) 23 homogeneous, the target homogeneous magnetic field distribution is expressed by the following.
Expression 6
{right arrow over (b)}u (6)
Further, the measured values (or the calculated values) of the magnetic field distribution on the current imaging region (homogeneous magnetic field space) 23 are expressed by the following.
Expression 7
{right arrow over (b)}m (7)
Then, the magnetic field distribution to be generated can be calculated by the following expression.
Expression 8
{right arrow over (b)}={right arrow over (b)}
u
−{right arrow over (b)}
m (8)
If the magnetic moment distribution is obtained, then the volumes Vi of shim bolts 27 corresponding to respective magnetic moments mi can be simply calculated by the following expression from Expression (1).
Expression 9
V
i
=m
i
/M (9)
More concretely, these show an example of the distribution of the magnetic moments obtained by Expression (5), wherein
In such a manner, the magnetic moments mi are the amounts distributed at the respective nodes. Herein, a positive magnetic moment refers to a magnetic moment in the same direction as that of the magnetic field caused by the magnetic device 50, and a negative magnetic moment refers to a magnetic moment in the opposite direction to that of the magnetic field caused by the magnetic device 50.
As described above, reproducing the distribution of the magnetic moments (in other words, the volumes of the shim bolts 27) distributed at the respective nodes (or the respective screw holes 26) as it is on the shim trays 17, 18 makes the work efficiency significantly low due to a large number of nodes. In this situation, in order to increase the work efficiency, it is considered to define regions in a state of the grid lines (orthogonal grid) 28 such as to be a background shown in
Further, the volume VA of the shim bolts 27 is obtained by the following expression.
The unit (dimension) of the respective numbers is, for example, the cubic centimeter, for representation of the volume of shim bolts 27. A positive number means that a magnetic moment whose value corresponds to the value of this number is to be given in the direction reinforcing the static magnetic field caused in this region by the magnetic device 50. Concretely, shim bolts 27 of ferromagnetic pieces (iron or the like) corresponding to the volume represented by this number are to be given. A negative number means that a magnetic moment corresponding to the value of this number is to be given in the direction attenuating the static magnetic field caused in this region by the magnetic device 50. Concretely, shim bolts 27 of permanent magnets with a strength corresponding to this number are to be disposed in the opposite direction to that of the static magnetic field caused by the magnetic device 50, or, if shim bolts 27 of ferromagnetic pieces are already disposed, these are to be removed. The size (area, shape) of each section partitioned by the grid lines (orthogonal grid) 28 is appropriately determined in advance such as to have a sufficient performance for magnetic field adjustment work, and is, for example, 50 millimeters square. By this method, the workability is improved compared with disposing a shim bolt/bolts 27 at each node.
As shown in
In this situation, as shown in
Expansion of this principle to a real two-dimensional computational mesh will be as follows.
First, all peak positions Pn are extracted from the distribution of magnetic moments. The magnetic moment at a certain node i will be represented by mi, and the magnetic moment at each adjacent node j will be represented by mj. If the following expression is satisfied for every adjacent node j, the node i is the peak position.
Expression 12
mi>mj>0 or mi<mj<0 (12)
Next, while checking the values of the magnetic moments at the adjacent nodes with initial points at the respective peak positions Pn, the boundaries of the regions An are determined as follows.
Expression 13
0<mk<mo or mo<0 (13)
Further, if the peak position Pn has a negative magnetic moment, and the following expression is satisfied, then the nodes o are defined to be on the (n+1)th layer.
Expression 14
mo<mk<0 or mo>0 (14)
If any one of nodes o do not satisfy these expressions, then redefinition is made such that the node k instead of the nodes o is defined to be on the (n+1)th layer.
By the above-described procedure, regions An corresponding to the peak positions Pn are determined, as schematically shown in
When the magnetic moment and the volume of the nodes belonging to the region An are calculated by Expression (10) and Expression (11), this volume VAn is the volume of shim bolts 27 to be disposed in the region An. The worker is only required to dispose VAn in the vicinity of the peak position Pn.
Concretely, in the vicinity of a peak position to be applied with a positive magnetic moment, shim bolts 27 of ferromagnetic pieces, iron for example, are to be disposed. In the vicinity of a peak position to be applied with a negative magnetic moment, shim bolts of permanent magnets are to be disposed in the direction for applying a negative magnetic moment, or, if there are already existing shim bolts 27, these are to be removed.
The unit (dimension) of the respective numbers is, for example, cubic centimeter for representation of the volumes of shim bolts 27. In order to decrease work errors by intuitively expressing positive or negative of volumes, the symbols shown on the left side of the numbers represent positive amounts with a mark “Δ” and negative amounts with a mark “∇”. In order to recognizably notify a worker of positions to dispose shim bolts 27, it is desired that, for example, the shim trays 17, 18 are partitioned by the grid lines (orthogonal grid) 28, and coordinates, as shown in
When the above-described displays are made by magnetic field adjustment software, the worker is only required to perform disposition, for example, at seven positions in the case of the example shown in
The above-described method approximately determines regions An and the boundaries thereof. Accordingly, in order to improve the accuracy of magnetic field homogeneity adjustment, this magnetic field homogeneity adjustment is repeated plural times.
Next, the flow of magnetic field homogeneity adjustment work in accordance with the present invention will be described.
First, the magnetic field strength distribution in the imaging region (the homogeneous magnetic field space) 23 is measured (step S1). Concretely, a magnetic field distribution measurement device 60 operates, and a magnetic probe 63 obtains measurement signals (a measurement result). Based on this measurement result, a data obtaining computer 61 generates magnetic field analysis data (magnetic field distribution data 72) (described later with reference to
Then, a magnetic field homogeneity adjustment computer 62 (described later with reference to
Next, the worker disposes shim bolts 27 on the shim trays 17, 18, while having a view of the display on a display device 65 (refer to
Then, similarly to step S1, the magnetic field strength distribution of the imaging region (the homogeneous magnetic field space) 23 is measured (step S4).
Next, it is determined whether or not the specification of the homogeneous magnetic field is satisfied (step S5). That is, it is determined whether or not the magnetic field homogeneity of the imaging region (the homogeneous magnetic field space) 23 is within a predetermined value. More concretely, the magnetic field homogeneity adjustment computer 62 (refer to
In this repeated process, if the display in
The magnetic field distribution measurement device 60 is provided with the magnetic probe 63 that is inserted into the imaging region (the homogeneous magnetic field space) 23 of the magnetic device 50 and detects the magnetic distribution, and the data obtaining computer 61 that is connected to the magnetic probe 63 and has the display device 65 and a data obtaining program installed thereon.
The magnetic field homogeneity adjustment computer 62 is a computer that includes the display device 65 and an output device 64, such as a printer, and has software for adjusting magnetic homogeneity installed thereon.
The magnetic field homogeneity adjustment computer 62 includes a storage device 66, a computing device 67, a display device 65, and an output device 64. On the magnetic field homogeneity adjustment computer 62, the software for adjusting magnetic homogeneity is installed, and forms functions of an input section 73, a computation section 74, a display method generation section 75, and an output section 76.
The operations from obtaining the magnetic field distribution to displaying shim bolts 27 to be disposed are carried out as follows.
With reference to
In such a manner, as the correspondence between the positions (coordinates) to dispose shim bolts 27 and the volumes thereof is clear even without drawing the positions in the magnetic field distribution, the workability is significantly improved while keeping the principle of adding the distributed volumes on the shim tray 17 (18).
With reference to
While the magnetic device 50 (refer to
As shown in
By performing calculation which is similar to that conducted in the first embodiment, using such a computational mesh, magnetic field homogeneity adjustment work is all the same possible also on the magnetic device 51 with the structure shown in
In respective embodiments in accordance with the present invention, a worker for magnetic field homogeneity adjustment is only required to dispose shim bolts (magnetic shims) 27 in a minimum quantity at positions of a minimum requirement in respective stages of a magnetic field homogeneity adjustment work, which eliminates the necessity of managing all the positions of shim trays 17, 18 and disposing magnetic shims precisely at the respective positions, and thus the efficiency of the magnetic field homogeneity adjustment work can be significantly increased. Further, because a magnetic device 50 or the like using such a method or an MRI apparatus using such a apparatus can reduce the time for installation adjustment, it is possible to provide an inexpensive apparatus as a result.
Number | Date | Country | Kind |
---|---|---|---|
2008-123345 | May 2008 | JP | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/JP2009/058712 | 5/8/2009 | WO | 00 | 11/8/2010 |