This application is based upon and claims the benefit of priority from prior Japanese Patent Applications No. 2006-248402, filed Sep. 13, 2006; and No. 2007-133193, filed May 18, 2007, the entire contents of both of which are incorporated herein by reference.
1. Field of the Invention
The present invention relates to a magnetic resonance image diagnostic apparatus fit for providing susceptibility-weighted images of the interior of, for example, the head, and also to a method of controlling the apparatus.
2. Description of the Related Art
T2*-weighted images acquired by the gradient echo method acutely reflect the non-uniformity of a local magnetic field. Therefore, the gradient echo method is widely used as an imaging method for visualizing a difference between regions of interests in terms of susceptibility. In order to provide T2*-weighted images of the subject's head, rephasing is generally performed by means of gradient moment nulling (GMN), nullifying the influence that the blood flow imposes on the image quality. T2*-weighted images are thereby acquired. The longer the echo time, the more acutely any T2*-weighted image will reflect the difference in susceptibility. This is why the echo time is set relatively long in most cases to provide T2*-weighted images. An imaging method that is more sensitive to changes in susceptibility than T2*-weighted imaging is proposed in Magn Reson, Med. 52:612-618, 2004 (hereinafter referred to as the reference document). In this imaging method, it is possible to enhance both effects of amplitude reduction and phase difference induced by susceptibility effects because phase-shift weighted mask is multiplied by amplitude images.
where γ is a gyromagnetic ratio of about 2π×42.6 MHz/T. G(t) is the gradient waveform vector, which corresponds to Gss, Gpe and Gro for the three axes (i.e., slice, phase encode and read-out), respectively, and r0, v0 and a0 are position, velocity and acceleration vectors of the spin at time t (t=0), respectively.
In equation (1), the respective terms indicate phase changes according to the position, velocity and acceleration, and correspond to the 0th-order, 1st-order and 2nd-order gradient moments. In equation (1), the third-order moment or any higher moment are not described. Nonetheless, the third-order moment and any higher-order moment contribute to the phase change, too.
GMN is a process that determines G(t) so that any gradient moment up to a certain order may have the minimum value such as 0 for the echo time TE. GMN is therefore also called rephasing, too. In the phase encode or slice encode, however, the 0th moment changes every time an encode step is performed. Thus, in the GMN for the phase-encode axis, the 0th moment at the echo time TE is set to a specific value assigned to each encode step, and the first-order moment and any higher-order moment are set to as small value as possible, such as 0. The zeroth-order GMN alone cannot achieve rephasing on the first-order moment and any higher-order moment in equation (1) for, for example, the blood flow that has active spins. Inevitably, the phases changes as the object moves. Consequently, the phase dispersion takes place, reducing the vector sum of the blood-flow spins, and the blood flow is represented by a low-level signal. Since the phase disperses in various ways in accordance with the velocity at which the blood flows, some blood flows may undergo insufficient phase dispersion and may not be represented by signals of sufficiently low level. In some cases, the blood flows will result in artifacts. Hence, in order to provide a T2*-weighted image of, for example, the head, which is free of influence of any fluid flow such as blood flow, the first-order GMN or any higher-order GMN must be performed. How high the order of GMN should be depends on whether G(t) can be rephased within the echo time TE at that order. The higher the order of the rephasing, the more greatly the influence of the flow can be reduced.
The reference document proposes that a T2*-weighted image of the head, obtained by the gradient echo method of performing GMN of the first order in the three axes, be phase-weighted so that the image may be further susceptibility-weighted. This proposed method will be explained below.
The original amplitude image, which will be phase-weighed, is an image that reflects a difference in susceptibility because it is free of the influence of the blood flow, owing to the first-order GMN. In the phase data, the second-order moment and any higher-order moment is not rephased in equation (1). Nonetheless, the susceptibility is considered to dominantly control the phase shift of the spin since the phase shift is almost free of the influence of a flow such as the blood flow. The phase shift can therefore be regarded as representing the difference in susceptibility between different kinds of tissue. Thus, a phase mask image that lowers the signal value in proportion to the phase shift may be generated from the phase data and then repeatedly applied to the amplitude image several times. Then, an image further phase-shift-weighed in accordance with the susceptibility can be acquired. The image acquired by this method reflects the difference in susceptibility between different kinds of tissue in the form of a difference in contrast. In the head, the difference between the venous blood having a high deoxyhemoglobin concentration and the surrounding tissue having a high oxyhemoglobin concentration is large in terms of susceptibility. An image clearly visualizing the veins can therefore be provided if phase making is applied to the T2*-weighted image of the head. This method has been reported to be useful and effective if used in BOLD venography.
The technique disclosed in the reference document can indeed provide susceptibility-weighted images that visualize the veins more clearly than is possible before. However, these images do not visualize the arteries more clearly, because the arteries exhibit less non-uniformity of susceptibility than the veins. Inevitably, the effect of any flow cannot contribute to image contrast. With the technique disclosed in the reference document, the phase masking must be repeated a plurality of times in order to enhance the vein-visualizing ability as much as desired. If the phase masking is repeated, however, the phase difference between the veins and the surrounding tissue will be emphasized. Consequently, the vascular cavities of veins may be over-evaluated or the artifacts may be emphasized due to the susceptibility effect.
The technique disclosed in the reference document can certainly generate phase masks reflecting phase changes that dominate susceptibility only. However, second-order GMN or any higher-order GMN is not performed, and only first-order GMN is performed. Hence, the technique cannot prevent artifacts that develop from second-order moments and higher-order moments generated as the blood vessels undergo a pulsating motion or a complex motion.
In view of the foregoing, it has been firstly desired that images should visualize, with high precision, the effect of flow or the effect of flow and susceptibility.
It has been secondly desired that the influence any flow imposes on the phase shift should be further reduced thereby to enhance the vein-visualizing ability.
According to a first aspect of the present invention, there is provided a magnetic resonance imaging diagnostic apparatus comprising: a generating unit which generates a slice gradient magnetic field, a phase-encode gradient magnetic field and a read-out gradient magnetic field that extend in a slice axis, a phase-encode axis and a read-out axis, respectively; a setting unit which sets a dephase amount for weighting a signal-level decrease resulting from flows in the arteries and veins present in a region of interest of a subject, with respect to at least one axis selected form the slice axis, phase-encode axis and read-out axis; and a control unit which controls the generating unit by using a pulse sequence for a gradient echo system, which includes a dephase gradient-magnetic-field pulse that corresponds to the dephase amount set by the setting unit for said at least one axis.
According to a second aspect of the present invention, there is provided a magnetic resonance imaging diagnostic apparatus comprising: a generating unit which generates a slice gradient magnetic field, a phase-encode gradient magnetic field and a read-out gradient magnetic field that extend in a slice axis, a phase-encode axis and a read-out axis, respectively; and a control unit which performs the first-order rephasing to the nth-order rephasing (n is 2 or a greater integer) on a spin with respect to at least one axis selected from the slice, phase-encode and read-out axes, and which controls the generating unit to achieve a pulse sequence for a gradient echo system in which an echo period is set to a value appropriate for canceling a phase shift of a vein present in a region of interest of a subject.
According to a third aspect of the present invention, there is provided a magnetic resonance imaging diagnostic apparatus comprising: a generating unit which generates a slice gradient magnetic field, a phase-encode gradient magnetic field and a read-out gradient magnetic field that extend in a slice axis, a phase-encode axis and a read-out axis, respectively; a setting unit which sets a dephase amount for weighting a signal-level decrease resulting from flows in the arteries and veins present in a region of interest of a subject, with respect to at least one axis selected form the slice, phase-encode and read-out axes; and a control unit which controls the generating unit by using a pulse sequence for an asymmetric spin-echo system, which includes a dephase gradient-magnetic-field pulse that corresponds to the dephase amount set by the setting unit for said at least one axis.
According to a fourth aspect of the present invention, there is provided a magnetic resonance imaging diagnostic apparatus comprising: a generating unit which generates a slice gradient magnetic field, a phase-encode gradient magnetic field and a read-out gradient magnetic field that extend in a slice axis, a phase-encode axis and a read-out axis, respectively; and a control unit which performs the 1st-order to the nth-order rephasing (n is 2 or a greater integer) on a spin with respect to at least one axis selected from the slice, phase-encode and read-out axes, and which controls the generating unit to achieve a pulse sequence for an asymmetric spin-echo system in which an echo period is set to a value appropriate for canceling a phase shift of a vein present in a region of interest of a subject.
According to a fifth aspect of the present invention, there is provided a method of controlling a magnetic resonance imaging diagnostic apparatus comprising a generating unit which generates a slice gradient magnetic field, a phase-encode gradient magnetic field and a read-out gradient magnetic field that extend in a slice axis, a phase-encode axis and a read-out axis, respectively, the method comprising: setting a dephase amount for weighting a signal-level decrease resulting from flows in the arteries and veins present in a region of interest of a subject, with respect to at least one axis selected form the slice, phase-encode and read-out axes; and controlling the generating unit by using a pulse sequence for a gradient echo system, which includes a dephase gradient-magnetic-field pulse that corresponds to the dephase amount set with respect to said at least one axis in the setting a dephase amount.
According to a sixth aspect of the present invention, there is provided a method of controlling a magnetic resonance imaging diagnostic apparatus comprising a generating unit which generates a slice gradient magnetic field, a phase-encode gradient magnetic field and a read-out gradient magnetic field that extend in a slice axis, a phase-encode axis and a read-out axis, respectively, the method comprising: performing the first-order rephasing to the nth-order rephasing (n is 2 or a greater integer) on a spin with respect to at least one axis selected from the slice, phase-encode and read-out axes, and which controls the generating unit and controlling the generating unit achieving a pulse sequence for a gradient echo system in which an echo period is set to a value appropriate for canceling a phase shift of a vein present in a region of interest of an subject.
According to a seventh aspect of the present invention, there is provided a method of controlling a magnetic resonance imaging diagnostic apparatus comprising a generating unit which generates a slice gradient magnetic field, a phase-encode gradient magnetic field and a read-out gradient magnetic field that extend in a slice axis, a phase-encode axis and a read-out axis, respectively, the method comprising: setting a dephase amount for weighting a signal-level decrease resulting from flows in the arteries and veins present in a region of interest of a subject, with respect to at least one axis selected form the slice axis, phase-encode axis and read-out axes; and controlling the generating unit by using a pulse sequence for an asymmetric spin-echo system, which includes a dephase gradient-magnetic-field pulse that corresponds to the dephase amount set with respect to said at least one axis in the setting a dephase amount.
According to a eighth aspect of the present invention, there is provided a method of controlling a magnetic resonance imaging diagnostic apparatus comprising a generating unit which generates a slice gradient magnetic field, a phase-encode gradient magnetic field and a read-out gradient magnetic field that extend in a slice axis, a phase-encode axis and a read-out axis, respectively, the method comprising: performing the first-order rephasing to the nth-order rephasing (n is 2 or a greater integer) on a spin with respect to at least one axis selected form the slice phase-encode axis, phase-encode axis and read-out axes, and controlling the generating unit achieving a pulse sequence for an asymmetric spin-echo system in which an echo period is set to a value appropriate for canceling a phase shift of a vein present in a region of interest of an subject.
According to a ninth aspect of the present invention, there is provided a magnetic resonance imaging diagnostic apparatus comprising: a collecting unit which collects magnetic resonance signals generated in a subject lying in a static magnetic field, when a radio-frequency pulse is applied, while applying a slice gradient magnetic field, a phase-encode gradient magnetic field and a read-out gradient magnetic field in a slice axis, a phase-encode axis and a read-out axis, respectively; and a main-scan control unit which controls the collecting unit by dephasing a spin about at least one axis selected from the slice, phase-encode and read-out axes, and by using a pulse sequence for a gradient echo system, while changing a dephase amount for the spin at least once during the collection of the magnetic resonance signals concerning a image.
According to a tenth aspect of the present invention, there is provided a method of controlling a magnetic resonance imaging diagnostic apparatus comprising a collecting unit which collects magnetic resonance signals generated in a subject lying in a static magnetic field, when a radio-frequency pulse is applied, while applying a slice gradient magnetic field, a phase-encode gradient magnetic field and a read-out gradient magnetic field in a slice axis, a phase-encode axis and a read-out axis, respectively; and a main-scan control unit which controls the collecting unit by dephasing a spin about at least one axis selected from the slice axis, phase-encode axis and read-out axis, and by using a pulse sequence for a gradient echo system, while changing a dephase amount for the spin at least once during the collection of the magnetic resonance signals concerning a image.
Additional objects and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objects and advantages of the invention may be realized and obtained by means of the instrumentalities and combinations particularly pointed out hereinafter.
The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate embodiments of the invention, and together with the general description given above and the detailed description of the embodiments given below, serve to explain the principles of the invention.
First to third embodiments of the present invention will be described, with reference to the accompanying drawings.
The MRI apparatus 100 comprises a bed unit, a static-magnetic-field generating unit, a gradient-magnetic-field generating unit, a receiving/transmitting unit, and a control/operating unit. A subject 200 may lie on the bed unit. The static-magnetic-field generating unit generates a static magnetic field. The gradient-magnetic-field generating unit is designed to add position information to the static magnetic field. The receiving/transmitting unit can receive and transmit a radio-frequency signal. The control/operating unit controls the other components of the system and reconstruct images. The MRI apparatus 100 has, as components of these units, a magnet 1, a static magnetic power supply 2, a shim coil 3, a shim power supply 4, a top plate 5, a gradient coil unit 6, a gradient power supply 7, an RF coil unit 8, a transmitter 9T, a receiver 9R, a sequencer (sequence controller) 10, an arithmetic unit 11, a storage unit 12, a display 13, an input device 14, an audio unit 15, and a host computer 16. To the MRI apparatus 100, an electrocardiograph unit is connected and detects an ECG signal representing the cardiac phase of the subject 200.
The static-magnetic-field generating unit includes the magnet 1 and the static magnetic power supply 2. The magnet 1 is, for example, a superconducting magnet or a normal magnet. The static magnetic power supply 2 supplies a current to the magnet 1. The static-magnetic-field generating unit therefore generates a static magnetic field B0 in a cylindrical aperture (examination space) into which the subject 200 is moved. The direction of the static magnetic field B0 virtually coincides with the axis (i.e., the Z-axis) of the examination space. The static-magnetic-field generating unit further includes the shim coil 3. The shim coil 3 generates a correction magnetic field for rendering the static magnetic field uniform when a current is supplied to it from the shim power supply 4 under the control of the host computer 16.
The bed unit moves the top plate 5, on which the subject 200 is lying, into or out of the examination space.
The gradient-magnetic-field generating unit includes the gradient coil unit 6 and the gradient power supply 7. The gradient coil unit 6 is arranged in the magnet 1. The gradient coil unit 6 has three coils 6x, 6y and 6z that can generate three gradient magnetic fields extending in mutually orthogonal X-, Y- and Z-axes, respectively. The gradient power supply 7 supplies pulse currents to the coils 6x, 6y and 6z, under the control of the sequencer 10. Supplied with the pulse current, the coils 6x, 6y and 6z generate gradient magnetic fields. The gradient-magnetic-field generating unit controls the pulse currents supplied from the gradient power supply 7 to the coils 6x, 6y and 6z. Thus, the gradient-magnetic-field generating unit synthesizes gradient magnetic fields extending in the three physical axes (i.e., the X-, Y- and Z-axes), respectively. The unit sets these magnetic fields in logical axes defined by a slice gradient magnetic field Gss, a phase-encode gradient magnetic field Gpe and a read-out (frequency-encode) gradient magnetic field Gro, respectively, which intersect at right angles with one another. The slice, phase-encode and read-out gradient magnetic fields, Gss, Gpe and Gro, are superposed on the static magnetic field B0.
The receiving/transmitting unit includes the RF coil unit 8, the transmitter 9T, and the receiver 9R. The RF coil unit 8 is arranged near the subject 200 in the examination space. The transmitter 9T and the receiver 9R are connected to the RF coil unit 8. The transmitter 9T and the receiver 9R operate under the control of the sequencer 10. The transmitter 9T supplies an RF current pulse of Lamor frequency to the RF coil unit 8 in order to induce nuclear magnetic resonance (NMR). The receiver 9R acquires an MR signal (radio-frequency signal), such as an eco signal, which the RF coil unit 8 has received. The receiver 9R then performs, on the MR signal, various processes, such as pre-amplification, intermediate-frequency conversion, phase detecting, low-frequency amplification and filtering. Finally, the receiver 9R performs analog-to-digital (A/D) conversion on the MR signal, producing digital data (raw data).
The control/operating unit includes the sequencer 10, the arithmetic unit 11, the storage unit 12, the display 13, the input device 14, the audio unit 15 and the host computer 16.
The sequencer 10 has a CPU and a memory. The sequencer 10 receives pulse sequence information from the host computer 16. The pulse sequence information is stored into the memory. The CPU of the sequencer 10 controls the gradient power supply 7, transmitter 9T and receiver 9R in accordance with the sequence information stored in the memory. The CPU also receives the raw data output from the receiver 9R and transfers the raw data to the arithmetic unit 11. Note that the sequence information is all data necessary for operating the gradient power supply 7, transmitter 9T and receiver 9R in accordance with the pulse sequence. It includes, for example, information about the intensity of the pulse current supplied to the coils 6x, 6y and 6z, the period of supplying the pulse current and the timing of supplying the pulse current. The sequencer 10 has the function of controlling the gradient-magnetic-field generating unit under the control of the host computer 16, by using a pulse sequence for a gradient echo system, which contains a dephase gradient-magnetic-field pulse that corresponds to the dephase amount set by the host computer 16. Further, the sequencer 10 has the function of controlling the gradient-magnetic-field generating unit under the control of the host computer 16, causing the gradient-magnetic-field generating unit to generate magnetic resonance signals about the same slice, while changing the dephase amount.
The arithmetic unit 11 receives the raw data output from the transmitter 9T, through the sequencer 10. The arithmetic unit 11 has an internal memory. The internal memory has k-space (also known as Fourier space or frequency space), in which the raw data input to the arithmetic unit 11 is stored. The raw data is subjected to two- or three-dimensional Fourier transform, thereby reconstructing video data for the real space. The arithmetic unit 11 can perform, if necessary, synthesis and differential operation (including weighted differentiation) on any data representing an image. The synthesis includes cumulative addition of pixel values, maximum intensity projection (MIP), minimum intensity projection (minIP), and the like. As another example of the synthesis, the axes of several frames may be aligned in a Fourier space, and the raw data items representing these frames may be synthesized, thereby generating one-frame raw data. The addition of pixel values may be simple addition, arithmetic-mean process or weighted-mean process.
The storage unit 12 stores video data reconstructed or video data subjected to the above-mentioned synthesis or differential operation.
The display 13 can display various images to the operator, under the control of the host computer 16. The display 13 is, for example, a display device such as a liquid crystal display.
The input device 14 is operated to input various data items such as scanning conditions the operator desires, the pulse sequence, data about the image synthesis and differential operation, and the like. These data items are sent from the input device 14 to the host computer 16. The input device 14 comprises, as the case may be, a pointing device such as a mouse or a track ball, a section device such as a mode-switching panel, or a device such as keyboard.
The audio unit 15 generates messages for the start and end of breath holding as sounds under the control of the host computer 16.
The host computer 16 can perform various functions when a prescribed software routine is executed. One of these functions is to manage the units constituting the MRI apparatus 100 so that the apparatus 100 performs such various operations as the existing MRI apparatus performs. Another of the above-mentioned functions is to set a dephase amount so that the level decrease of a signal, due to the flows in the arteries and veins existing in the region of interest of the subject 200, may be weighted with respect to at least one of the slice, phase-encode and read-out axes. Still another of the above-mentioned functions is to control the sequencer 10 so that magnetic resonance signals may be collected by using a gradient-echo system pulse sequence that contains dephase gradient-magnetic-field pulses corresponding to a dephase amount with respect to the axis for which the dephase amount has been set. A further one of the above-mentioned functions is to determine a part for which the level decrease of a signal should be weighted. Another of the above-mentioned functions is to control the sequencer 10 so that a plurality of magnetic resonance signals about the same slice may be generated, while the dephase amount is being changed. Still another of the functions is to control the arithmetic unit 11 so that preparation images of a region of interest may be generated from a plurality of magnetic resonance signals. Another of the above-mentioned functions is to control the display 13, causing the display 15 to display a plurality of preparation images. A further one of the above-mentioned functions is to determine which one of preparation images the operator wants. Another of the above-mentioned functions is to reduce the phase shift determined from a magnetic resonance signal, by a value determined from the known susceptibility of the tissue, thereby to find the phase shift that corresponds to the flow of a tissue.
After performing a preparation work such as positioning scan, the host computer 16 performs imaging scan. The imaging scan is a scan for collecting a set of echo data items that are necessary for reconstructing an image. In this case, the imaging scan is the two-dimensional scan. The imaging scan can be carried out, along with the ECG gate method that uses ECG signals. In some cases, the ECG gate method need not be performed together with the imaging scan.
The electrocardiograph unit includes an ECG sensor 17 and an ECG unit 18. The ECG sensor 17 may be attached to the subject 200 to detect an ECG signal from the subject 200. The ECG signal is an electrical signal (hereinafter referred to as a sensor signal). The ECG unit 18 performs various processes, including binary encoding, on the sensor signal. The sensor signal thus processed is output to the sequencer 10. The electrocardiograph unit is, for example, a vector electrocardiograph. The sequencer 10 uses the sensor signal generated by the electrocardiograph unit, when it is necessary to carry out a scan in synchronism with the cardiac phase of the subject 200.
How the MRI apparatus 100 configured as described above operates will be explained in detail.
The first embodiment will be described. The first embodiment is designed to achieve the first desire.
In Step Sa1, the host computer 16 instructs the sequencer 10 to collect data. Upon receiving the instruction, the sequencer 10 collects data as will be explained below.
As seen from
Assume that the time at which the RF pulse is applied is 0 and that the echo time TE is 3. Base on this assumption, gradient magnetic fields having intensities of 1, −2 and 1 are applied at time intervals of 1 as shown in
The rephasing controls the gradient magnetic fields to cancel out the change in the signal, which occurs due to the influence of the flow. If rephased, the phase change of the blood flow is reduced significantly. Thus, the blood-flow signal is collected as a high-level signal, not decreased in level. Conversely, the dephasing causes the influence of the flow to induce changes in the signal. In the dephasing, the gradient magnetic fields promote the phase dispersion of the spin of a flow, such as the blood flow, which is a motion. The vector sum of the blood-flow spins is small in the dephasing. In other words, the amplitude component of the flow signal is greatly attenuated. The flow signal is therefore is suppressed and eventually collected as a low-levels signal.
In Step Sa2, the host computer 16 instructs the arithmetic unit 11 to reconstruct an image from the data collected as described above. Upon receiving this instruction, the arithmetic unit 11 performs, for example, a known image-reconstructing method. This image-reconstructing method obtains an amplitude image, a phase image and three-dimensional (3D) volume data.
In Step Sa3, the host computer 16 instructs the arithmetic unit 11 to perform interpolation on the image reconstructed as described above. On receiving the instruction, the arithmetic unit 11 performs the known interpolation process. The interpolation process may not be carried out.
In Step Sa4, the host computer 16 instructs the arithmetic unit 11 to create a synthesized image. Upon receiving this instruction, the arithmetic unit 11 creates a synthesized image from the amplitude image and the phase image. The arithmetic unit 11 creates, for example, a phase mask image that represents a part greatly phase-shifted and a part slightly phase-shifted. The phase mask image consists of zero-pixels defining a part of the image, which advances (or is delayed) in phase, and one-pixels defining the other part of the image, which has no phase difference. That is, the phase mask image represents a part of high susceptibility (having a phase difference) and a part of low susceptibility. The phase mask image is a phase image that may be filtered so that the low-frequency component may become free of phase disturbance. The filtering may be accomplished by a process of passing the phase image through a high-pass filter or a process of differentiating the phase image before and after passing the phase image through a low-pass filter. The arithmetic unit 11 multiplies the phase image thus obtained, by the amplitude image, a predetermined number of times (one to several times), thus providing an SWI image in which the phase change due to the susceptibility is further weighed. This multiplication may not be carried out.
In Step Sa5, the host computer 16 instructs the arithmetic unit 11 to create an image to display. Upon receiving this instruction, the arithmetic unit 11 creates an image to display, by the known method. This process creates data that represents the blood vessel as a continuous tube. The process may be a projection such as maximum intensity projection (MIP), minimum intensity projection (minIP), or integral projection. In the first embodiment, both the arteries and the veins are represented by signals that are lower in level than those representing the surrounding tissue. In view of this, minimum intensity projection (minIP) is optimal. The data representing the surface of any blood vessel may be extracted in order to apply volume rendering or surface rendering. Alternatively, multi-planer reconstruction (MPR) may be employed to enable the doctor to observe the original image signal.
In Step Sa6, the host computer 16 causes the display 13 to display the image created in step Sa5, amplitude image or phase image, in accordance with the operator's instruction.
Thus, in the first embodiment, the flow can be visualized accurately, without over-evaluating the vascular cavities of veins or weighting artifacts resulting from the effect of susceptibility.
In the first embodiment, artifacts may develop, depending on the dephase amount. The dephase amount should therefore be adjusted to such an appropriate value as would not cause artifacts and also would weight the image of blood vessels.
The dephase amount can be defined by the offset amount of the moment referred to in conjunction with
The b-factor is defined as inner product of the integration of the waveform vectors of gradient magnetic fields, as in the following equation (2) and (3).
The b-factor represents a decrease in the signal level, which has resulted from dephasing. It is therefore related to the dephase amount since. A coefficient D is applied to a particular tissue, and the signal level is said to fall in accordance with the b-factor as expressed by the equation S(b)=S(0)·exp(−bD), where S(0) is the signal not dephased yet and S(b) is the signal dephased.
As seen from
The relation between the dephase amount and the ability of visualizing the blood vessels differs, from subject to subject. The host computer 16 therefore has the function of setting the dephase amount to a value desirable for the subject 200. The host computer 16 may set a dephase amount in accordance with the specific value the operator has input by operating the input device 14. Instead, it may set the dephase amount by selecting one of the several dephase-amount data items (i.e., a plurality of candidates). Otherwise, the computer 16 may automatically set an appropriate dephase amount.
The value the operator input to set the dephase amount is, for example, the b-factor, velocity encoding (VENC), or the velocity. To help the operator to determine which value should be input, preparation images may be displayed on the display 13. For example, the b-factor is changed several times, thus performing imaging the object several times and providing several images thereof, and these images may then be displayed on the display 13. More specifically, such images as shown in
To set the an appropriate dephase amount automatically, the host computer 16 may select an image having a contrast-to-noise ratio (CNR) or a signal-to-noise ratio (SNR) that falls within a tolerance range, from the preparation images obtained as described above. Then, the host computer 16 may set a dephase amount on the basis of the b-factor or the VENC used to acquire the preparation image thus selected. There may be several preparation images that accord with the conditions described above. If this is the case, the host computer 16 may select any one of the preparation images. However, the greater the dephase amount is, the higher the ability of visualizing the blood vessels will be. It should therefore be better for the host computer 16 to select an image acquired by using a large b-factor or a small VENC.
In the first embodiment, the ability of visualizing the flow never changes in accordance with the echo time TE. Thus, TE can be changed as is desired.
The echo time TE is related to the influence the non-uniformity of a local magnetic field imposes on the image contrast. That is, the shorter the echo time TE, the less the non-uniformity of the local magnetic field will influence the image contrast. Hence, TE may be set to such a small value as would not reflect the non-uniformity of the local magnetic field on the image contrast. Then, an image can be acquired, in which the effect of flow is dominant. Conversely, the TE may be set to such a large value as would reflect the non-uniformity of the local magnetic field on the image contrast. In this case, an image can be acquired, which reflects both the effect of flow and the effect of susceptibility.
The echo time TE is a period that extends from the RF excitation to the appearance of a peak in the echo signal. The peak appears in the echo signal when the integral value of the gradient magnetic field Gro decreases to zero. The echo time TE can therefore be adjusted by changing the time TA at which the gradient magnetic field Gro is applied after the RF excitation without changing the pattern of the gradient magnetic field Gro, as is shown in
An image acquired in a relatively short echo time TE, i.e., an image in which the effect of the flow is dominant, images the blood vessels accurately and clearly. Such the image helps a doctor to examine the blood vessels. In the first embodiment, the attenuation or phase change of the signal is used to visualize the flow. The image acquired can clearly visualizes the thin arteries and the collateral channels extending around from above and can therefore clinically useful in diagnosing, for example, cerebral infarction.
Cerebral infarction can be cured if a thrombus-dissolving chemical is applied within a few hours of the development of the infarction, particularly if the infarction is acute. If the patient is left untreated, the cerebral tissue will undergo necrosis over time. Treatment must therefore be started before the affected cerebral tissue become irreversibly damaged. The scan time can be shortened by reducing the number of phase encodes, namely by lowering the resolution.
However, thrombi, for example, must be imaged at a sufficiently high signal-to-noise ratio or a high spatial resolution and within a very short scan time, in order to attain reliable diagnosis.
Therefore, the RF coil unit 8 shown in
In principal, the signal-to-noise ratio decreases in the PI method, in inverse proportion to the increase rate of imaging speed. Nonetheless, the PI method should rather be used to shorten the time, while maintain the spatial resolution, since 3D images have a higher signal-to-noise ratio than 2D images.
The PI method is a technique of omitting the phase encode step and performing another step instead, thereby generating an image folded from the incomplete raw data and reconstructing an image from the raw data. As is known in the art, the PI method is classified into sensitivity encoding (SENSE) and simultaneous acquisition spatial harmonics (SMASH), which differ in terms of image-reconstructing process. Either type of PI method can be utilized in this embodiment. SENSE develops the folded part extending out of the field of view (FOV), by using the sensitivity distribution of the coil elements. SMASH first generates complete raw data from the raw data the coil elements, by using the sensitivity distribution of the coil elements, and then reconstructs an image from the complete raw data.
If the PI method (particularly, SENSE) is employed, phase data items will be acquired for the coil elements, respectively. The phase data items thus acquired are synthesized by means of, for example, weighted addition, providing synthesized data. From the synthesized data a phase mask image is formed for use in the phase-shift weighting that is performed in Step Sa4. Thus, an effective phase-shift weighting can be accomplished in the PI method, too.
The phase change resulting from the effect of susceptibility is influenced not only by the echo time TE, but also by the intensity of the static magnetic field. The value appropriate for the echo time TE in the first embodiment therefore changes in accordance with the intensity of the static magnetic field. The echo time TE may be set to about 40 ms at 1.5T, and to about 20 ms at 3T. Then, the blood vessels become flow-void because of the phase dispersion that has resulted from dephasing. As a result, an image will be obtained, which has a contrast fully reflecting the effect of susceptibility will be obtained. To induce a positive phase dispersion by means of dephasing, thereby to lower the level of the signal, a T2*-weighted image can be acquired within a short TE. This shortens the imaging time and increases the signal-to-noise ratio in the acquisition of T2*-weighted image.
A T2*-weighted image, in which the blood-flow spin induces a low-level signal because of the phase dispersion, can certainly be acquired by performing GMN involving the zeroth-order moment only, without performing rephasing using the first-order moment or any higher-order moment as in the conventional method. However, GMN using the zeroth-order moment only cannot achieve a sufficient phase dispersion of the blood flow. Consequently, the resultant image of the blood vessel consists of low-level signal components and high-level signal components where the phase dispersion is insufficient. Inevitably, the blood vessel cannot be less well visualized than is possible with the first embodiment of the present invention.
An image acquired by using a relatively long echo time TE, i.e., image formed by virtue of the effect of flow and the effect of susceptibility, provides information about thrombi and bleeding, thanks to the effect of susceptibility, in addition to the information about the blood vessels. This image therefore enables the doctor not only to determine whether thrombi exist or bleeding is occurring, but also to locate the thrombi or the bleeding. That is, the image gives the doctor information that is useful in planning an effective treatment.
In the conventional method, too, wherein rephasing is carried out, the longer the echo time TE, the more greatly the non-uniformity of the local magnetic field will be emphasized. The phase dispersion is therefore large in the will increase in any T2*-weighted image acquired for a long echo time TE. Hence, the T2*-weighted image has high ability of visualizing the blood vessels indeed. However, not only the blood flow, but also the parenchyma is susceptibility-weighted. As a consequence, the unevenness or distortion will be conspicuous in the image, due to the non-uniformity of the local magnetic field, or the signal-to-noise ratio will decrease or the imaging time will increase, because of the long echo time TE.
The conventional gradient echo method can indeed provide a T2*-weighted images that clearly visualizes an event, such as bleeding, which accompanies a large change of susceptibility. If rephasing is effected, nullifying the influence of the blood flow, the blood flow will generate a high-level signal. As a result, the image contrast is not so high that the cerebral parenchyma and that of the blood vessels cannot be well distinguished. In other words, the bleeding site and the position of any blood vessel cannot be clearly confirmed in the same image.
The phase-shift data acquired from an echo in the first embodiment represents a phase shift that has developed due to susceptibility and a flow. Unless the phase-shift data is appropriately processed, it cannot show how high the susceptibility is or how fast the flow is. Nevertheless, phase-shift data in which only the flow is dominant can be obtained by subtracting the phase shift specific to the susceptibility of the tissue from the phase shift due to the flow. In this regard, it should be noted that the susceptibilities of various kinds of tissue are known.
The image acquired in the first embodiment described above is a magnetic resonance angiograph (MRA) that can be obtained without the necessity of using any contrast medium. That is, the first embodiment can provide such a useful image as described above, without applying a contrast medium. Good use can therefore be made of non-invasiveness that is an advantageous feature of magnetic resonance imaging.
As is known in the art, in such susceptibility-weighted imaging (SWI) as disclosed in the reference document, the phase polarities inside and outside any blood vessel may change due to the magic angle effect, depending on the angle between the direction in which the static magnetic field extends and the direction in which the blood vessel runs. Inevitably, in the conventional SWI, any part adjacent to a vein is represented by a low-level signal when minIP is performed. In this case, one vein will appear as two parallel veins as is illustrated in
The second embodiment will be described. The second embodiment is designed to achieve the second desire.
As seen from
Thus, in the second embodiment, only the susceptibility is dominant with respect to phase changes and phase data can enhance the ability of visualizing the veins.
The first-order moment and the second-order moment need not be strictly zero. Rather, it suffices for the first-order moment to have such a value as will cancel the amplitude attenuation or phase shift resulting from the flow.
The imaging method used in the magnetic resonance imaging diagnostic apparatus 100 according to the first embodiment, i.e., the imaging method using a pulse sequence for the gradient echo system, which includes a dephase gradient-magnetic-field pulse, will be hereinafter referred to as flow-sensitive black blood (FSBB) imaging.
In the FSBB imaging, the phase of each flow such a blood flow, which is a motion, is gradually distributed, and the vector sum of the flow spins gradually decreases. That is, the amplitude component of the flow signal is gradually attenuated. As a result, the signal of the flow is suppressed and an image will be obtained, in which any blood vessel is represented by a low-level signal. Generally, a b-factor is used, which indicates a dephase amount. The b-factor is a tensor amount defined by equations (2) and (3) above.
In equation (3), γ is a gyromagnetic ratio of about 2π×42.6 MHz/T, and G(t) is a gradient-magnetic-field waveform vector. G(t) corresponds to Gss, Gpe and Gro for the axes of slice, phase encode and read-out, respectively. The b-factor indicates a signal-level decrease and defines a relation between signals, i.e., S(b)=S(0)·exp(−bD), where S(0) is a signal not dephased, S(b) is a signal dephased, and D is a coefficient for a specific tissue. When signals have this relation, the signal level may decrease in accordance with the b-factor.
In the FSBB imaging, a calcification in the brain is represented by a low-level signal. Hence, if the FSBB imaging is applied to image thin blood vessels, such as penetrated branches, the images of the calcifications in the brain will overlap the images of the blood vessels. In this case, the blood vessels may not be clearly detected in the MRI image. This influence of calcification can be minimized by, as is known in the art, reducing the echo time TE.
In the FSBB imaging, a dephase amount equal to or greater than a specific value must be applied in order to image blood vessels sufficiently. To preserve such a dephase amount, however, the gradient magnetic field must be intensified in inverse proportion to the echo time TE. Therefore, echo time TE cannot be shortened too much.
The intensity the gradient magnetic field should have in order to maintain the dephase amount at a particular value decreases as the encode step proceeds, as can be understood from, for example,
Since the dephase amount is limited, the shortening of the echo time TE is limited, too. Thus, it is difficult to reduce the decrease of TE sufficiently.
The third embodiment can therefore minimize the influence of the calcifications, successfully visualizing the blood vessels clearly.
In the third embodiment, the sequencer 10 includes the function of controlling the gradient power supply 7, transmitter 9T and receiver 9R to accomplish the FSBB imaging. More specifically, the sequencer 10 has the function of controlling the gradient power supply 7, transmitter 9T and receiver 9R by using a pulse sequence for the gradient echo system, while dephasing the spin about at least one of the slice axis, phase-encode axis and read-out axis. While performing this control to achieve this FSBB imaging, the sequencer 10 changes the dephase amount for the spin about at least one axis, at least one time, during the collection of magnetic resonance signals about one image.
The sequencer 10 further has the function of controlling the gradient power supply 7, transmitter 9T and receiver 9R, before the main scan is performed to collect magnetic resonance signals for the FSBB imaging. The preparation scan is thereby carried out several times, by using a different dephase amount each time, thus collecting magnetic resonance signals concerning the same slice.
The input device 14 has the function of inputting a preparation-image selecting instruction (described later) under the control of the host computer 16.
The host computer 16 has various functions. One of the functions is to make the arithmetic unit 11 reconstruct a plurality of preparation images on the basis of the magnetic resonance signals collected over the several preparation scans and to cause the display 13 to display the preparation images thus reconstructed. One of the above-mentioned functions is to set a target dephase amount that is used to collect the magnetic resonance signals for use in reconstructing any preparation image that the operator has selected by operating the input device 14.
How the MRI apparatus 100 operates in the third embodiment will be explained. The various imaging processes the MRI apparatus 100 performs as the existing MRI apparatuses do will not be described, nonetheless. Only the imaging process for acquiring the image of the penetrated branches will be explained here.
In the process loop consisting of Steps Sb1 to Sb4, the host computer 16 repeats Steps Sb2 and Sb3, incrementing the variable i from 1 to m, each time by one (1).
In Step Sb2, the host computer 16 sets a prescribed value b(i) as b-factor. Value b(i) is one of, for example, six values b(1) to b(6), i.e., 0.1, 1, 4, 16, 32 and 64 sec/mm2. In this case, variable m is “6.”
In Step Sb3, the host computer 16 instructs the sequencer 10 to perform a preparation scan for one slice. Upon receiving the instruction, the sequencer 10 controls the gradient power supply 7, transmitter 9T and receiver 9R. The magnetic resonance signals the receiver 9R has received are thereby collected. The preparation scan is performed, first by fixing the b-factor at the value set in Step Sb2 and then by using the pulse sequence for the gradient echo system, while dephasing the spin about at least one of the slice, phase-encode and read-out axes.
The process loop consisting of Steps Sb1 to Sb4 is thereby carried out. At this time, the b-factor is changed for the same slice m times, performing m preparation scans. Therefore, m magnetic resonance signals corresponding to the different b-factors are collected.
Next, in Step Sb5, the host computer 16 instructs the arithmetic unit 11 to reconstruct preparation images. On receiving this instruction, the arithmetic unit 11 reconstructs m preparation images on the basis of the m magnetic resonance signals, one preparation image from one signal. Therefore, m preparation images corresponding to different b-factors are acquired. The preparation images are FSBB images, because the spin about at least one axis is dephased during the preparation scan and the magnetic resonance signals for the blood vessels are at low levels.
In Step Sb6, the host computer 16 causes the display 13 to display the preparation images. At this point, all preparation images, i.e., m preparation images, may be displayed as shown in
As can be seen from
In Step Sb7, the host computer 16 instructs the input device 14 to input the data indicating that the best preparation image has been selected.
In Step Sb8, the host computer 16 sets a target b-factor as b-factor b(i) for the preparation image selected.
In the process loop consisting of Steps Sb9 to Sb15, the host computer 16 repeats Steps Sb10 to Sb14 several times, increasing variable j from 1 to n, by one (1) each time Steps Sb10 to Sb14 are performed. Note that n is the number of phase encode steps performed to collect magnetic resonance signals for one slice, and is, for example, “256.”
In Step Sb10, the host computer 16 finds the intensity Gpe(j) that the phase-encode gradient magnetic field Gpe should have to achieve the target b-factor in the jth phase encode step.
In Step Sb11, the host computer 16 determines whether the intensity Gpe(j) is equal to or higher than the maximum intensity Gpe_max for the phase-encode gradient magnetic field. The maximum intensity Gpe_max is determined by the abilities of the coils 6x, 6y and 6z and the ability of the gradient power supply 7. If the intensity Gpe(j) higher than the maximum intensity Gpe_max, the host computer 16 sets the intensity of the phase-encode gradient magnetic field Gpe in accordance with the maximum intensity Gpe_max. More precisely, the host computer 16 first calculates the b-factor for the case where the maximum intensity Gpe_max is set for the phase-encode gradient magnetic field Gpe. Then, the host computer 16 sets the b-factor calculated, as a dephase amount for the encode step. The host computer 16 calculates the intensity Gpe(j) on the basis of the b-factor and sets this b-factor as the intensity of the phase-encode gradient magnetic field Gpe. If the intensity Gpe(j) is not higher than the maximum intensity Gpe_max, the host computer 16 sets, in Step Sb13, the intensity Gpe(j) calculated in Step Sb10, as the intensity of the phase-encode gradient magnetic field Gpe.
In Step Sb14, the host computer 16 instructs the sequencer 10 to collect magnetic resonance signals concerning the jth encode step. Upon receiving the instruction, the sequencer 10 controls the gradient power supply 7, transmitter 9T and receiver 9R so that the magnetic resonance signals concerning the jth encode step may be collected. The magnetic resonance signals the receiver 9R has received are thereby collected. This collection of magnetic resonance signals is accomplished by a pulse sequence for the gradient echo system, while dephasing the spin of the phase-encode axis for which the intensity of the phase-encode gradient magnetic field Gpe has been set in Step Sb12 or Sb13.
The process loop consisting of Steps Sb9 to Sb15 is thereby carried out. At this time, all magnetic resonance signals necessary in n phase encode steps for one slice are collected. This process of collecting magnetic resonance signals is the main scan. The gradient magnetic field for attain the target b-factor may change as indicated by the broken line shown in
Hitherto, the b-factor has been usually not changed during the scanning. This is because the dephasing has been generally used in order to perform diffusion imaging in the scanning in which the b-factor is used as an input parameter, and the b-factor should better be constant to facilitate the calculation of a diffusion coefficient. In the third embodiment, however, the b-factor is changed in the main scan, unlike in the conventional MIR apparatus.
Although not shown in the drawing, it may be necessary to perform the main scan so that slice encoding may be achieved for a plurality of slices or in 3D imaging. In this case, the process loop consisting of Steps Sb9 to Sb15 is carried out for each phase encode step or slice encode step. At this point, the slice gradient magnetic field Gss is set so that the b-factor may become constant, or set to have the same intensity as the phase-encode gradient magnetic field Gpe.
In most cases, the intensity of a gradient magnetic field exceeds the upper limit in any encode step, when it changes in the slice direction or the phase-encode direction, and does not exceed the upper limit when it changes in the read-out direction. Therefore, it is desired to fix the dephase amount when the intensity of the gradient magnetic field changes in the read-out direction, and to change the dephase amount to the upper limit, only in such an encode step in which the intensity of the gradient magnetic field may exceed the upper limit.
In Step Sb16, the host computer 16 instructs the arithmetic unit 11 to reconstruct the image on the basis of the data acquired as described above. On receiving the instruction, the arithmetic unit 11 reconstructs the image by, for example, a known method. This image reconstruction provides an amplitude image, a phase image and 3D volume data.
In Step Sb17, the host computer 16 instructs the arithmetic unit 11 to perform interpolation on the image reconstructed. On receiving the instruction, the unit 11 performs the known interpolation process. The interpolation process may not be carried out.
In Step Sb18, the host computer 16 instructs the arithmetic unit 11 to synthesize an image. Upon receiving this instruction, the unit 11 creates a synthesized image from, for example, the amplitude image and the phase image. The process of creating the amplitude image may not be performed.
In Step Sb19, the host computer 16 instructs the arithmetic unit 11 to create an image to display. On receiving this instruction, the arithmetic unit 11 creates an image to display, by the known method. The image thus created shows the blood vessel as a continuous tube. This process of creating an image may be a projection such as the maximum intensity projection (MIP), the minimum intensity projection (minIP), or the integral projection. The data representing the surface of any blood vessel may be extracted in order to apply volume rendering or surface rendering. Alternatively, the multi-planer reconstruction (MPR) may be employed to enable the doctor to observe the original image signal.
In Step Sb20, the host computer 16 causes the display 13 to display the image created in step Sb19, amplitude image or phase image, in accordance with the operator's instruction.
In the third embodiment, the dephase amount is smaller than the target b-factor in the first phase encode step, but is gradually increased to the target b-factor as the phase encode process goes on. Thus, the dephase amount can be larger than in the case it is set as indicated by the broken line in
In the third embodiment, the b-factor is so changed as to become maximal in each encode step, not so simply as in the case of providing the image shown in
The present embodiments can be modified in various ways, as will be described below.
In the first embodiment, the dephase amount and the polarity need not be constant. Dephasing need not be performed in three axes. It may rather be performed in two axes or one axis. If this is the case, rephasing may be performed in the remaining axis or axes. Further, dephasing may be performed, not in all encodes, but in 0 encode only. Still further, the axes for dephasing or the encode steps may be limited.
In the second embodiment, the axes for rephasing and the encode steps may be changed. Moreover, the third-order moment and any higher-order moment may be large enough to cancel the amplitude attenuation or the phase shift, either due to the flow.
In the first and second embodiments, the pulse sequence may be a two-dimensional (2D) one, not a 3D one as exemplified above. The pulse sequence for an asymmetric spin-echo system can be used in place of the pulse sequence for the gradient echo method. As shown in
In the third embodiment, of the gradient magnetic fields for all encode steps, the most intensive one may be set to a value smaller than the upper limit.
In the third embodiment, the b-factor for any gradient magnetic field may be changed only once, in order to provide such an image as shown in
In the third embodiment, the target b-factor may be a value that the operator directly selects.
In the third embodiment, the signal-to-noise ratio (SNR) at the head may be automatically measured for each dephase amount, and the target b-factor may set to such a value as would render this SNR maximal.
In the third embodiment, any preparation scan need not be performed for the same time as the main scan. Rather, the preparation scan may be performed for a shorter time than the main scan, by reducing the number of encodes or by imaging only one slice by means of a 2D scan so that the doctor may visually determine the degree of visualizing the blood vessels.
In the third embodiment, any preparation scan may be performed by a method other than changing the dephase amount. For example, the preparation scan may be performed by the phase contrast method that can provide such an image as shown in
Additional advantages and modifications will readily occur to those skilled in the art. Therefore, the invention in its broader aspects is not limited to the specific details and representative embodiments shown and described herein. Accordingly, various modifications may be made without departing from the spirit or scope of the general inventive concept as defined by the appended claims and their equivalents.
Number | Date | Country | Kind |
---|---|---|---|
2006-248402 | Sep 2006 | JP | national |
2007-133193 | May 2007 | JP | national |