GB2182450A - Nuclear magnetic resonance imaging apparatus - Google Patents
Nuclear magnetic resonance imaging apparatus Download PDFInfo
- Publication number
- GB2182450A GB2182450A GB08625116A GB8625116A GB2182450A GB 2182450 A GB2182450 A GB 2182450A GB 08625116 A GB08625116 A GB 08625116A GB 8625116 A GB8625116 A GB 8625116A GB 2182450 A GB2182450 A GB 2182450A
- Authority
- GB
- United Kingdom
- Prior art keywords
- image
- ratio
- brightness
- standard deviation
- mean
- 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/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/58—Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
In a method of assessing the performance of an NMR imaging apparatus which facilitates comparison of the performance of the apparatus when operating in different modes, inserts 30 of a phantom are filled with materials of known different properties, e.g. T1 or T2 (for example manganese chloride solutions of different strengths), or temperature. For two selected inserts with values of T1 or T2 at the ends of a range of interest, a figure of merit is given by the suppression <IMAGE> where P is proton density (mean image brighteners) for the respective insert, DELTA LnC is the difference of the natural logs of the concentrations in the inserts and equals (T'-T''), and sigma /M is the mean of sigma /M, where sigma is the standard deviation of brighteners over the area of the insert image and M is the mean brightness of the insert image. <IMAGE>
Description
SPECIFICATION
Nuclear magnetic resonance imaging apparatus
This invention relates to nuclear magnetic resonance (NMR) imaging apparatus.
In order to assess the performance of such apparatus use is normally made of a so-called phantom, that is to say, an object so constructed that an assessment of apparatus performance can readily be made from a visual inspection of the image of the object produced by the apparatus and/or quantitative analysis of data produced by the apparatus for use in presentation of the image. Such data typically comprises the brightness value of each pixel of the image produced by the apparatus, pixel brightness being a function of the total NMR signal detected by the apparatus as arising from the part of the object being imaged corresponding to that pixel.
A known form of NMR phantom comprises an enclosure made of a material which produces an effectively zero NMR response when imaged by the apparatus. Within the enclosure there may be disposed various suitable shaped members of the same or another zero NMR response material, the whole of the remainder of the interior of the enclosure being filled with a material of known NMR response. The material of known NMR response is required to have a highly homogeneous distribution of protons of the kind which are to be excited to resonance in use of the apparatus. Where, as is usual, the protons to be excited are hydrogen protons, a suitable material of known NMR response for use in an NMR phantom is a copper sulphate solution or magnanese chloride solution or other solution containing free hydrogen ions.
Since NMR imaging apparatuses are normally required to image a planar slice of an object, thereby to provide a cross-sectional view of the object in the plane of the slice, or in the case where imaging of a volume of an object is required, a plurality of parallel spaced planar slices disposed within the volume, NMR imaging phantoms typically have a tubular enclosure, the phantom being disposed in use so that the apparatus under test displays images of one or more planar slices of the phantom extending parallel to the plane ends of the enclosure.
The shaped members of the phantom, if present, are designed with a view to assessment of various parameters of the image. For example, to assess the spatial resolution of the image the shaped members typically comprise a number of spaced rods of different diameters extending parallel to the axis of the enclosure, and/or a series of plates disposed in planes parallel to the axis of the enclosure, the plates being of different thicknesses and different spacings.
To assess the image under different conditions of contrast, different regions of the phantom enclosure may be of different internal length i.e. so as to contain different depths of the material of known NMR response.
However, there does not exist at the present time a satisfactory method of comparing the performance of an apparatus when operating in different modes.
It is an object of the present invention to provide such a method.
According to the present invention a method of assessing the performance of an NMR imaging apparatus comprises obtaining a measure of the ratio of the contrast of the image produced by the apparatus for two different values of a chosen parameter; obtaining a measure of the standard deviation of the brightness of the image produced at said different values of said chosen parameter; and obtaining a measure of the ratio of said contrast ratio and said standard deviation.
Said chosen parameter is typically the spin to spin or spin to lattice relaxation time constant of material of a body being imaged.
One method in accordance with the invention will now be described by way of example with reference to the accompanying drawings in which:
Figure 1 is a plan view of a phantom used in the method;
Figures 2, 3 and 4 are sectional views of the phantom along the lines Il-Il, Ill-Ill, and IV-IV in
Fig. 1 respectively;
Figures 5 and 6 are graphs illustrating a first method of using the first phantom to determine the width of an imaged slice;
Figure 7 is a diagram illustrating said first method; and
Figure 8 is a diagram illustrating a second method of using the first phantom to determine the width of an imaged slice.
Referring to Figs. 1 to 4, the phantom has a short tubular housing made of a transparent plastics material, the housing comprising a cup-shaped body portion 1 and a lid 3 which is secured to the body portion by screws 5 of plastics material, the lid being sealed to the body portion by an '0' ring 6.
Within the housing there are disposed two right-angled wedge-shaped blocks 7 and 9 of plastics material. The blocks 7-and 9 are of the same length, height and wedge-angle, but one block 9 if half the width of the other block 7. The wider block 7 is secured to the base of the body portion 1 at one side of the housing with its base in contact with the base of the body portion 1. The narrower block 9 is secured to the lid 3 of the housing so that its inclined surface is in parallel spaced relationship with the inner half of the inclined surface of the wider block 7 when the lid 3 is secured to the body portion 1 of the housing.
The housing also contains a number of rods 11 of the plastics material of different diameters, the rods being secured to the lid 3 so as to be disposed to one side of the blocks 7 and 9 in parallel spaced relationship, with their axes parallel to the axis of the housing when the lid 3 is secured to the body portion 1 of the housing.
To the other side of the blocks 7 and 9, secured to the lid 3, there are four sets 13, 15, 17 and 19 of three parallel spaced plates of plastics material, the plates in each set being of a different thickness and spacing equal to their thickness.
At the side of the housing opposite the blocks 7 and 9, secured to the base of the body portion 1, there is a further, generally rectangular block 21 of plastics material. In the block 19 there are formed three holes 23, 25 and 27. One hole 23 extends centraly through the block 21 in a direction parallel to the axis of the housing. The other two holes 25 and 27 extend through the block 19 on opposite sides of the central hole 23 at the same angle, typically about 20 , to the base of the body portion 1, although the two holes are oppositely directed.
Where the holes 25 and 27 meet the surface of the block 21 the surface of the block is bevelled so as to be perpendicular to the axis of the holes.
Between the block 21 and the blocks 7 and 9 there is a set of thin wedge-shaped plates 29 of plastics material, the plates being secured to the base of the body portion 1 in a fan arrangement with their narrower ends meeting in a line extending parallel to the axis of the housing adjacent the block 21.
At each end of the blocks 7 and 9 and at each end of the block 21 there is disposed one of four identically shaped tubular inserts 30 of plastics material filled with a material of known NMR response, the materials in the four inserts 30 having different NMR responses. The inserts 30 are secured within holes in the lid 3 with their axes parallel to the axis of the housing, the inserts being fitted with removable plugs 30A to allow the material in the inserts to be changed.
The whole of the space within the housing not occupied by members of plastics materials is filled with a solution having a known NMR response, e.g. a solution of copper sulphate or manganese chloride in de-ionised water. The solution may be doped with sodium chloride so that the NMR response of the solution more nearly simulates that of a human body, the sodium chloride serving to distribute the NMR responsive free hydrogen ions in the solution.
A hole is provided in the curved wall of the body portion 1 of the housing to permit filling the housing with solution, the hole being plugged by a screw 31 of plastics material.
The plastics material is chosen to have a substantially zero NMR response compared with the solution used to fill the housing.
In use the phantom is disposed in the examination chamber of an NMR apparatus under test, and the apparatus is operated to image a planar slice extending through the interior of the housing of the phantom in a direction parallel to the flat ends of the phantom.
An assessment of the performance of the apparatus is made by visual and semi-quantitative assessment of the images of the rods 11, and plates 13, 15, 17, 19 and 29 in conventional manner. It will be appreciated in this connection that if the image indicates NMR response on a grey scale with zero NMR response corresponding to black, the plastic objects all appear black on the image surrounded by a bright display.
The narrower block 9 and the opposed part of the wider block 7 provide a visual indication of the position of the imaged slice with respect to the phantom. Thus the area of the image corresponding to the narrower block 9 and the opposed part of the block 7 will appear black except for a bright bar extending parallel to the width of the block 9 corresponding to the region where the image slice intersects the gap between the blocks 7 and 9. As best seen from Fig.
2, as the image slice is moved from the base of the housing towards the lid 3, the bright bar moves from one end of the block 7 to the other.
The part of the wider block 7 not opposed by the narrower block (see Fig. 1) is used for making a quantitative assessment of the width of the image slice using a method forming the subject of UK Patent Application No. 8504809 (Specification No. GB-A-2,155,187).
Referring now to Fig. 7, in the region where the image slice intersects the inclined surface of the block 7, the proportion of the width of the image slice which is occupied by material of the block 7 varies from 1 to O in a direction towards the thinner end of the block 7. Thus the corresponding image varies from black to maximum brightness over this region.
Slice width is assessed by taking the brightnesses P of pixels of the image on a line extending in a direction X parallel to the length of the block 7 across the above-mentioned region where the image slice intersects the inclined face of the block 7. It will be appreciated that these brightness values are preferably taken from data computed to display the image rather than by brightness measurements made on the displayed image. The brightness values are then plotted on a graph against pixel number in the direction X, as illustrated in Fig. 5 by way of example.
The slope dP/dX of the resulting graph is then plotted against pixel number in the direction X, as illustrated in Fig. 6.
Slice width may then be quantitatively assessed by measurement on the graph of Fig. 6.
Typically slice width will be measured where the value of dP/dX has half its maximum value.
The graph of Fig. 6 is used to obtain the X-direction dimension A at this level. The slice width
W is then obtained from the expression A Tan 6 where 6 is the wedge-angie of the block 7, as illustrated in Fig. 7.
It will be appreciated that the graph of Fig 6 in fact represents the image slice thickness profile.
The measured pixel brightness values represent the edge response function projection of the slice in the X-direction which is a convolution of the slice thickness profile with a ramp function determined by the wedge-angle of the block 7. Hence the slice profile is obtained when the edge response function projection is deconvolved by differentiation.
The graph of Fig. 6 also permits a quantitative assessment of the percentage of the area of the slice profile lying within the chosen limits, e.g. the half maximum dP/dX level.
It will be noted that the graph of Fig. 6 also gives a useful indication of the spurious out-ofslice NMR responses obtained with the apparatus.
In conventional methods of assessing the performance of NMR imaging apparatus slice width is assessed by Fourier analysis of the detected response signal or the RF pulses used to excite resonance in the slice. It will be appreciated that a phantom as described above enables a direct much more accurate assessment of slice width to be obtained, since measurement is made on the actual image obtained.
The phantom illustrated in Figs. 1 to 4 is suitable for assessing the performance of an NMR imaging apparatus when imaging a single slice of a body. It will be appreciated that for assessing performance when imaging a volume of a body i.e. a plurality of parallel spaced slices, a phantom of greater length is required, e.g. about 90 centimetres in length. The blocks 7 and 9 will then have a large wedge-angle e.g. about 60 instead of about 30 as shown. The position and width of each slice and the inter-slice spacing of the slices can then all be measured as described above.
The image of block 21 may be used in conventional manner as a partial volume effect gauge, i.e. to allow contrast ratio of the image to be assessed.
However, the image of the block 21 may also be used to obtain a rapid indication of slice width, as follows:
Referring to Fig. 8, so long as the slice width is greater than the dimension of the holes 25 and 27 in the direction of the width of the slice, then slice width W is given by the expression
W= /s where S is the ratio of the maximum brightness of the image of the hole 25 or 27 (i.e. at the centre of the image) to the brightness of the image of the hole 23. It will be appreciated that is given by B/cos where B is the diameter of the holes 25 and 27 and is the angle that the holes 25 and 27 make with the base of the body portion 1 of the housing.
The inserts 30 are used in accordance with the present invention to obtain a better assessment than is possible using the block 21 of the contrast ratio of the image, i.e. the difference in brightness between images of materials having known different properties, e.g. a different spin to lattice relaxation time constants T1 or spin to spin relaxation time constants T2.
To this end the materials filling the inserts are suitably manganese chloride solutions having strengths in the range 1.00x 1017Mni + ions/millilitre to 1.00 X 1015 Mn + ions/millilitre. In one embodiment used in practice the four solution are of strengths 1.0x10l7 (T1=160ms, T2=380 ms at 6.4 MHz) 4.20x1016 (T1=220ms, T2=310 ms) 1.40x1016 (T1=423 ms, T2= 160 ms) and 6.93x1015 (T1=700 ms, T2=86 ms)
A quantitative assessment of the contrast ratio of the image produced by an apparatus in respect of T1 or T2 can be obtained using the expression
PT PT PT where PT is the proton density, i.e. mean image brightness, for the material of an insert 30 with the apparatus operating in a mode in which image brightness is primarily dependent on the chosen relaxation time constant T1 or T2, and the suffixes ' and " indicate respective ones of two selected inserts 30 whose T1 or T2 values lie at the ends of a range of T1 or T2 values
of interest.
A figure of merit for an apparatus under test in respect of T1 or T2 may be obtained using the expression:
where LnC is the difference of the natural logs of the concentrations of the solutions in the two selected inserts 30 and /M is the mean of T/M for the two selected inserts 30, a being the standard deviation of the brightness over the area of the image of an insert, and M being the mean value of the brightness over the area of the image of an insert.
It will be appreciated that ALnC=T'T".
Hence, the figure of merit is the ratio of the contrast ratio of the image, expressed as a percentage of the chosen range of the selected time constant T1 or T1, and the standard deviation of the image brightness, expressed as a percentage of the mean of image brightness.
For a given pair of inserts 30 expression (1) may be reduced to:
where k is a constant.
It will be appreciated that the inverse of the term n/M gives a measure of the signal to noise ratio of the image. It will thus be understood that the higher is the 'figure of merit' obtained then the better may one expect to be able to distinguish visually between areas of materials of different relaxation time constant T1 or T2 in the image obtained.
It will be understood that whilst in the above example a figure of merit has been obtained in respect of the performance of the apparatus for distinguishing between areas of an image of different T1 or T2 value, the method may also be used to.obtain a figure of merit in respect of the performance of an apparatus for distinguishing between areas of an image having different values of another parameter e.g. temperature.
Claims (8)
1. A method of assessing the performance of an NMR imaging apparatus comprising: obtaining a measure of the ratio of the contrast of the image produced by the apparatus for two different values of a chosen parameter; obtaining a measure of the standard deviation of the brightness of the image produced at said different values of said chosen parameter; and obtaining a measure of the ratio of said contrast ratio and said standard deviation.
2. A method according to Claim 1 wherein said ratio of the contrast is the ratio of the difference of the brightnesses of the images produced for said two different values of the chosen parameter and the brightness of the image for one of said two different values.
3. A method according to Claim 1 or Claim 2 wherein said measure of the ratio of the contrast is said ratio of the contrast expressed as a fraction of the difference between said two values of said chosen parameter.
4. A method according to any one of Claims 1 to 3 wherein said measure of the standard deviation is said standard deviation expressed as a fraction of the mean of the image brightness.
5. A method according to Claim 4 wherein said measure of the standard deviation is the mean of the ratio of the standard deviation and mean of the brightness of the image over a first area of a body being imaged by the apparatus in which first area said chosen parameter has one of said two values and the ratio of the standard deviation and mean of the brightness of the image over a second area of the body being imaged by the apparatus in which second area said chosen parameter has the other of said two values.
6. A method according to any one of Claims 1 to 5 wherein said chosen parameter is the spin to spin relaxation time constant (T2) of material of a body being imaged.
7. A method according to any one of Claims 1 to 5 wherein said chosen parameter is the spin to lattice relaxation time constant (T1) of material of a body being imaged.
8. A method of assessing the performance of an NMR imaging apparatus for two different values of a chosen parameter according to Claim 1 and substantially as hereinbefore described by way of example.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| GB08625116A GB2182450B (en) | 1984-02-27 | 1986-10-20 | Nuclear magnetic resonance imaging apparatus |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| GB848405065A GB8405065D0 (en) | 1984-02-27 | 1984-02-27 | Nuclear magnetic resonance imaging apparatus |
| GB08625116A GB2182450B (en) | 1984-02-27 | 1986-10-20 | Nuclear magnetic resonance imaging apparatus |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| GB8625116D0 GB8625116D0 (en) | 1986-11-26 |
| GB2182450A true GB2182450A (en) | 1987-05-13 |
| GB2182450B GB2182450B (en) | 1987-10-14 |
Family
ID=26287375
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| GB08625116A Expired GB2182450B (en) | 1984-02-27 | 1986-10-20 | Nuclear magnetic resonance imaging apparatus |
Country Status (1)
| Country | Link |
|---|---|
| GB (1) | GB2182450B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009004297A2 (en) * | 2007-06-29 | 2009-01-08 | King's College Hospital Nhs Foundation Trust | Phantom for imaging apparatuses |
| GB2512416A (en) * | 2013-03-28 | 2014-10-01 | Elekta Ab | Markers, phantoms and associated methods for calibrating imaging systems |
-
1986
- 1986-10-20 GB GB08625116A patent/GB2182450B/en not_active Expired
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009004297A2 (en) * | 2007-06-29 | 2009-01-08 | King's College Hospital Nhs Foundation Trust | Phantom for imaging apparatuses |
| WO2009004297A3 (en) * | 2007-06-29 | 2009-03-12 | King S College Hospital Nhs Fo | Phantom for imaging apparatuses |
| GB2512416A (en) * | 2013-03-28 | 2014-10-01 | Elekta Ab | Markers, phantoms and associated methods for calibrating imaging systems |
| GB2512384A (en) * | 2013-03-28 | 2014-10-01 | Elekta Ab | Markers, Phantoms and Associated Methods for Calibrating Imaging Systems |
| GB2512416B (en) * | 2013-03-28 | 2015-12-30 | Elekta Ab | Markers, phantoms and associated methods for calibrating imaging systems |
| GB2512384B (en) * | 2013-03-28 | 2016-07-20 | Elekta Ab | Markers, Phantoms and Associated Methods for Calibrating Imaging Systems |
| US10022104B2 (en) | 2013-03-28 | 2018-07-17 | Elekta Ab (Publ) | Markers, phantoms and associated methods for calibrating imaging systems |
Also Published As
| Publication number | Publication date |
|---|---|
| GB2182450B (en) | 1987-10-14 |
| GB8625116D0 (en) | 1986-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4613819A (en) | Nuclear magnetic resonance imaging apparatus | |
| US4551678A (en) | Phantom for nuclear magnetic resonance machine | |
| Goodsitt et al. | Real‐time B‐mode ultrasound quality control test procedures. Report of AAPM Ultrasound Task Group No. 1 | |
| US4056970A (en) | Ultrasonic velocity and thickness gage | |
| EP3082603B1 (en) | Universal phantom structure for quality inspections both on computerized tomography and on magnetic resonance tomography | |
| US4888555A (en) | Physiological phantom standard for NMR imaging and spectroscopy | |
| US6370480B1 (en) | Quantitative analysis system and method for certifying ultrasound medical imaging equipment | |
| US4644276A (en) | Three-dimensional nuclear magnetic resonance phantom | |
| US4777442A (en) | NMR quality assurance phantom | |
| CA2064841A1 (en) | Methods and apparatus for quantifying tissue damage | |
| US5165050A (en) | Spherical test body for an image reconstructing apparatus | |
| GB2182450A (en) | Nuclear magnetic resonance imaging apparatus | |
| GB2155187A (en) | Nuclear magnetic resonance imaging apparatus | |
| Alasaarela et al. | Evaluation of image quality of ultrasound scanners in medical diagnostics. | |
| US11692959B2 (en) | System and method for imaging tendon cross sections for detecting voids and other deficiencies in grouted external tendons | |
| EP0026673A1 (en) | A method for measurement of temperature distribution over the surface of a sample, and a system for performing such measurement | |
| Chui et al. | Test method for MR image slice profile | |
| JPS59157547A (en) | Phantom for measuring nmr-ct tomographic surface distortion | |
| GB2182449A (en) | Nuclear magnetic resonance imaging apparatus | |
| Plewes et al. | The influence of partial volume averaging on sphere detectability in computed tomography | |
| US3379971A (en) | Microwave means and method for measuring the electrolytic content of solutions | |
| Lewellen et al. | A low-contrast phantom for daily quality control | |
| C atphan® 500 and 600 M anual | ||
| C atphan® 500 and 600 Product Guide | ||
| JPH06277199A (en) | Phantom for evaluating magnetic resonance image diagnosis apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PCNP | Patent ceased through non-payment of renewal fee |