An ECG-prep scan is used to set an optimum time phase in both systole and diastole of the heart. At each of the different time phases, an imaging scan is started to acquire a plurality of sets of echo data. An artery/vein visually separated blood flow image is produced from the echo data. The imaging scan uses a half-Fourier technique, for example. This provides high-quality blood flow images with shorter scan time, without injecting a contrast medium. Additionally, with a readout gradient pulse applied substantially parallel with a direction of slowly flowing blood, a scan is performed in synchronism with an optimally determined cardiac time phase. The readout gradient pulse has a dephasing pulse for enhancing differences in a flow void effect depending on blood flow velocities. This enables slow-speed flows, such as blood flows in the inferior limb, to be depicted without fail.
|
16. An mr imaging method of obtaining an image relating to fluid within a region to be imaged of an object, comprising:
setting two different cardiac time phases falling into a systole and a diastole of a cardiac cycle of the object;
performing toward the region to be imaged of the object, an mr imaging scan starting in turn at each of the two cardiac time phases to acquire two sets of echo data, the mr imaging scan comprising a first scan starting at one of the two cardiac time phases falling in the systole and a second scan starting at the other of the two cardiac time phases falling in the diastole, both of the first scan and the second scan being based on a half-Fourier technique; and
producing, from the two sets of acquired echo data, the image relating to the fluid.
40. An mr imaging method of obtaining an image relating to fluid within a region to be imaged of an object, comprising:
setting a cardiac time phase of an object;
performing, toward the region to be imaged of the object, a scan at the cardiac time phase with use of a pulse sequence including a readout gradient pulse of which applied direction is substantially in accordance with a moving direction of fluid in motion within the object, so that an echo signal is acquired; and
producing, from the echo signal, the image relating to the fluid,
wherein the readout gradient pulse has a main pulse to read out the echo signal and at least one of a dephase pulse and a rephase pulse responsible for dephasing and rephasing phases of magnetic spins of the fluid, respectively, the at least one pulse being added to the main pulse along a time axis of the main pulse.
45. An MRI system for obtaining an image relating to fluid within a region to be imaged of an object, comprising:
time phase setting means for setting two different cardiac time phases falling into a systole and a diastole of a cardiac cycle of the object;
scanning means for performing, toward the region to be imaged of the object, an mr imaging scan starting in turn at each of the two cardiac time phases set by the time phase setting means to acquire two sets of echo data, the mr imaging scan comprising a first scan starting at one of the two cardiac time phases falling in the systole and a second scan starting at the other of the two cardiac time phases falling in the diastole, both the first scan and the second scan being based on a half-Fourier technique; and
image producing means for producing, from the two sets of echo data acquired by the scanning means, the image relating to the fluid.
1. An MRI system for obtaining an image relating to fluid within a region to be imaged of an object, comprising:
a time phase setting unit configured to set two different cardiac time phases falling into a systole and a diastole of a cardiac cycle of the object;
a scanning unit configured to perform, toward the region to be imaged of the object, an mr imaging scan starting in turn at each of the two cardiac time pahses set by the time phase setting unit to acquire two sets of echo data, the mr imaging scan comprising a first scan starting at one of the two cardiac time phases falling in the systole and a second scan starting at the other of the two cardiac time phases falling in the diastole, both the first scan and the second scan being based on a half-Fourier technique; and
an image producing unit configured to produce, from the two sets of echo data acquired by the scanning unit, the image relating to the fluid.
17. An MRI system for obtaining an image relating to fluid within an object, in which the object placed in a static magnetic field is subjected to a scan based on a pulse sequence including a readout gradient pulse, comprising:
a time phase setting unit configured to set a cardiac time phase of the object;
a scanning unit configured to perform the scan at the cardiac time phase to acquire an echo signal from the object under a condition that an applied direction of the readout gradient pulse is substantially in accordance with a moving direction of the fluid in motion within the object; and
an image producing unit configured to produce, from the echo signal, the image relating to the fluid,
wherein the readout gradient pulse has a main pulse for reading out the echo signal and a control pulse for controlling behaviors of magnetic spins of the fluid concerning a phase of the magnetic spins, the control pulse being added to the main pulse along a time axis of the main pulse.
26. An MRI system for obtaining an image relating to fluid within an object, in which the object placed in a static magnetic field is subjected to a scan based on a pulse sequence including a readout gradient pulse, comprising:
a time phase setting unit configured to set a cardiac time phase of the object;
a scanning unit configured to perform the scan at the cardiac time phase to acquire an echo signal from the object under a condition that an applied direction of the readout gradient pulse is substantially in accordance with a moving direction of the fluid in motion within the object; and
an image producing unit configured to produce, from the echo signal, the image relating to the fluid,
wherein the time phase setting unit is configured to set two cardiac time phases of the object,
the scanning unit is configured to acquire data comprising two sets of echo signals by applying first and second scans to the object at the two cardiac time phases, respectively; and
the image producing unit is configured to produce an image of the fluid from the data.
0. 46. An MRI system for obtaining an image relating to fluid region to be imaged of an object, comprising:
a cardiac phase setting component configured to set a first cardiac phase and a second cardiac phase of a cardiac cycle of the object as first and second timings, the first and the second cardiac phase corresponding to a systole and a diastole, respectively;
a scanning component configured to perform a first 3d scan at the first timing to acquire a first echo data set and a second 3d scan at the second timing to acquire a second echo data set;
a computer configured to generate a pulse sequence such that both of the first 3d scan and the second 3d scan are performed with a gradient magnetic field including a dephase pulse so as to reduce an intensity of a signal from an artery in the systole by promoting dephasing effect, while giving reduced dephasing effect to both the artery and a vein in the diastole; and
an image producing component configured to produce a subtraction image from a first image and a second image to extract an image of the artery, the first image being generated based on the first echo data set, and the second image being generated based on the second echo data set.
0. 55. An MRI system for obtaining an image relating to fluid region to be imaged of an object, comprising:
a time phase setting component configured to set a first time phase and a second time phase of a pulse wave of the object as first and second timings, the first and the second time phase corresponding to a systole and a diastole, respectively;
a scanning component configured to perform a first 3d scan at the first timing to acquire a first echo data set from an inferior limb and a second 3d scan at the second timing to acquire a second echo data set from the inferior limb;
a computer configured to generate pulse sequences such that the first 3d scan is performed with a gradient magnetic field including a dephase pulse so as to reduce an intensity of a signal from an artery by promoting dephasing effect, while giving reduced dephasing effect to a vein in the systole, and the second 3d scan is performed with a gradient magnetic field including a rephase pulse so as to increase an intensity of signals from both the artery and the vein in the diastole; and
an image producing component configured to produce a subtraction image from a first image and a second image to extract an image of the artery, the first image being generated based on the first echo data set, the second image being generated based on the second echo data set.
0. 54. An mr imaging method of obtaining an image relating to fluid within a region to be imaged of an object, comprising:
performing, with respect to the region, a preparation scan for acquiring echo data sets at cardiac phases at a same slice position;
producing preparation image sets at the cardiac phases based on the echo data sets;
selecting a first preparation image and a second preparation image from the preparation image sets so that an image in which a vein is depicted alone in the highest intensity is determined to be the first preparation image and an image in which an artery and a vein are both depicted in the highest intensities is determined to be the second preparation image;
performing, toward the region, a first 3d imaging scan at a first cardiac phase which corresponds to a cardiac phase for the first preparation image and a second 3d imaging scan at a second cardiac phase which corresponds to a cardiac phase for the second preparation image, the first and the second cardiac phase corresponding to a systole and a diastole, respectively, a first echo data set being acquired at the first cardiac phase, a second echo data set being acquired at the second cardiac phase, both the first 3d imaging scan and the second 3d imaging scan being performed with a gradient magnetic field including a dephase pulse so as to reduce an intensity of a signal from an artery in the systole by promoting dephasing effect, while giving reduced dephasing effect to both the artery and a vein in the diastole; and
producing a subtraction image from a first image and a second image to extract an image of the artery, the first image being generated based on the first echo data set, the second image being generated based on the second echo data set.
41. An MRI system for obtaining an image relating to fluid within a region to be imaged of an object, comprising:
a magnet for generating a static magnetic field in which an object is placed;
an rf coil device through which an rf magnetic field is transmitted to the object and an echo signal emanated from the object is received;
a transmitter for transmitting the rf magnetic field to the object through the rf coil device, the rf magnetic field being based on a pulse sequence;
a gradient power supply for applying a gradient based on the pulse sequence to the object through a gradient coil;
a receiver for receiving the echo signal through the rf coil device, the echo signal being generated in response to performance of the pulse sequence;
a calculating unit for producing the echo signal received by the receiver into the image; and
a controller for controlling operations of the transmitter, receiver and gradient power supply in conformity with the pulse sequence,
wherein the controller controls the operations of transmitter, receiver and gradient power supply so that two different cardiac time phases falling into a systole and a diastole of a cardiac cycle of the object are set, and, as the pulse sequence, an imaging scan is executed in synchronism with each of the two different cardiac time phases in turn to acquire two sets of the echo signal, the imaging scan comprising a first scan starting at one of the two cardiac time phases falling in the systole and a second scan starting at the other of the two cardiac time phases falling in the diastole, both of the first scan and the second scan being based on a half-Fourier technique, and
the calculating unit produces the image relating to the fluid within the region to be imaged of the object from the two sets of the echo signal acquired correspondingly to each of the two different cardiac time phases.
42. An MRI system for obtaining an image relating to fluid within a region to be imaged of an object, comprising:
a magnet for generating a static magnetic field in which an object is placed;
an rf coil device through which an rf magnetic field is transmitted to the object and an echo signal emanated from the object is received;
a transmitter for transmitting the rf magnetic field to the object through the rf coil device, the rf magnetic field being based on a pulse sequence;
a gradient power supply for applying a gradient based on the pulse sequence to the object through a gradient coil;
a receiver for receiving the echo signal through the rf coil device, the echo signal being generated in response to performance of the pulse sequence;
a calculating unit for producing the echo signal received by the receiver into the image; and
a controller for controlling operations of the transmitter, receiver and gradient power supply in conformity with the pulse sequence,
wherein the controller controls the operations of transmitter, receiver and gradient power supply so that a cardiac time phase of the object is set and, as the pulse sequence, a pulse sequence for an imaging scan is executed in synchronism with the cardiac time phase, the imaging-scan pulse sequence including a readout gradient pulse of which applied direction being substantially in accordance with a moving direction of fluid in motion within the object,
the calculating unit produces the image relating to the fluid within the object from the echo signal acquired through the receiver correspondingly to performance of the imaging-scan pulse sequence, and
wherein the readout gradient pulse has a main pulse to read out the echo signal and a control pulse to control behaviors of magnetic spins of the fluid concerning a phase of the magnetic spins, the control pulse being added to the main pulse along a time axis of the main pulse.
0. 56. An mr imaging method of obtaining an image relating to fluid within a region to be imaged of an object, comprising:
performing, with respect to the region, a preparation scan for acquiring echo data sets from an inferior limb at time phases of a pulse wave of the object at a same slice position;
producing preparation image sets at the time phases based on the echo data sets;
selecting a first preparation image and a second preparation image from the preparation image sets so that an image in which a vein is depicted alone in the highest intensity is determined to be the first preparation image and an image in which an artery and vein are both depicted in the highest intensities is determined to be the second preparation image;
performing, toward the region, a first 3d imaging scan at a first time phase which corresponds to a cardiac phase for the first preparation image and a second 3d imaging scan at a second time phase which corresponds to a cardiac phase for the second preparation image, the first time phase the second time phase corresponding to a systole and a diastole, respectively, a first echo data set being acquired from the inferior limb at the first time phase, a second echo data set being acquired from the inferior limb at the second time phase, the first 3d scan being performed with a gradient magnetic field including a dephase pulse so as to reduce an intensity of a signal from an artery by promoting dephasing effect, while giving reduced dephasing effect to a vein in the systole, and the second 3d scan being performed with a gradient magnetic field including a rephase pulse so as to increase an intensity of signals from both the artery and the vein in the diastole; and
producing a subtraction image from a first image and a second image to extract an image of the artery, the first image being generated based on the first echo data set, the second image being generated based on the second echo data set.
2. The MRI system of
3. The MRI system of
the second scan is carried out using a pulse sequence generating an echo signal to map echo data in both a central region and one of both end regions other than the central region of a second k-space for producing the image, the central region corresponding to a lower-frequency region in a phase-encode direction of the second k-space and both of the end regions corresponding to a higher-frequency region in the phase-encode direction of the second k-space.
4. The MRI system of
5. The MRI system of
6. The MRI system of
7. The MRI system of
8. The MRI system of
9. The MRI system of
10. The MRI system of
11. The MRI system of
12. The MM system of
13. The MM system of
14. The MM system of
15. The MRI system of
18. The MRI system of
19. The MRI system of
20. The MRI system of
21. The MRI system of
22. The MRI system of
23. The MRI system of
25. The MRI system of
27. The MRI system of
28. The MRI system of
29. The MRI system of
the image producing unit has artery/vein image producing means that produces images in which the artery and vein are visually separated from each other.
30. The MRI system of
31. The MRI system of
the second scan is carried out using a pulse sequence generating an echo signal to map echo data in both of a central region and one of both end regions other than the central region of a second k-space for producing the image, the central region corresponding to a lower-frequency region in a phase-encode direction of the second k-space and both of the end regions corresponding to a higher-frequency region in the phase-encode direction of the second k-space.
32. The MRI system of
33. The MRI system of
34. The MRI system of
35. The MRI system of
36. The MRI system of
37. The MRI system of
38. The MRI system of
39. The MRI system of
43. The MRI system of
44. The MRI system of
the imaging scan consists of a first scan and a second scan made to start at the two cardiac time phases, respectively.
0. 47. The MRI system of claim 46, wherein slice encodes for the first and second 3d scans have the same amount.
0. 48. The MRI system of claim 46, wherein the second 3d scan is performed based on a pulse sequence, which has a flow dephasing pulse in a read-out direction.
0. 49. The MRI system of claim 46, wherein the second 3d scan is performed based on a pulse sequence which has a readout gradient pulse for reading out the second echo data set and a flow dephasing as the dephase pulse added to the readout gradient pulse.
0. 50. The MRI system of claim 46, wherein the first 3d scan is performed based on a pulse sequence which has a readout gradient pulse for reading out the first echo data set and a control pulse for controlling magnetic spins of the fluid.
0. 51. The MRI system of claim 46, wherein the first and second 3d scans are performed under a condition in which the direction of the readout gradient pulse substantially corresponds to a blood vessel direction to be imaged.
0. 52. The MRI system of claim 46, wherein the first and second 3d scans are performed based on FSE (Fast Spin Echo) technique.
0. 53. The MRI system of claim 46, wherein the subtraction is a weighted subtraction.
|
1. Field of the Invention
The present invention relates to magnetic resonance imaging (MRI) for internally imaging an object to be examined on the basis of a magnetic resonance phenomenon of nuclear spins of the object, particularly, to an MRI (magnetic resonance imaging) system and an MR (magnetic resonance) method capable of acquiring artery/vein visually separated images of the object without using a contrast medium.
2. Description of the Related Art
Magnetic resonance imaging is based on an imaging technique for magnetically exciting nuclear spins of an object located in a static magnetic field by applying a radio-frequency (RF) signal of a Larmor frequency and reconstructing an image from MR signals induced by the excitation.
For clinically obtaining blood flow images of the pulmonary field or abdomen of a patient by magnetic resonance imaging, MR angiography has been put in practical use, in which a contrast medium is injected into the object to highlight blood flows. However, this contrast MR angiography needs an invasive treatment to inject the contrast medium. First of all, mental and physical burdens on patients become large. Second, examination cost of the contrast MR angiography is still expensive. Third, there are some cases where a contrast medium cannot be injected into patients due to patient's physical characteristics.
In cases the contrast medium cannot be injected or is not used, imaging techniques, such as time-of-flight (TOF) and phase contrast (PC) techniques are used alternatively.
The time-of-flight method and phase contrast techniques utilize an effect of flows such as blood flows. The effect of flows is attributed to either of two natures possessed by spins in motion. One is that spins simply move their positions due to flows, while the other results from phase shifts of transverse magnetization caused when spins move in a gradient field. The nature of the position movement is used for the TOF technique and the nature of the phase shifts is used for the phase contrast technique.
However, when the TOF technique or phase contrast technique is used for obtaining MR images of a patient's pulmonary field or abdomen which depict flows of large vessels, such as the aorta, in their superior-inferior directions, it is required to scan slices located vertically to the flowing direction. That is, axial images should be acquired with a slice direction of those axial images set to the superior-inferior direction. Thus, in the case that two-dimensional slice imaging is performed to acquire such axial images, it is impossible to obtain an image in which blood flows are directly reflected. Three-dimensional image data spatially containing blood flows are therefore needed, but the number of slices increases which will cause an entire imaging time to be longer.
A novel MR imaging technique, known as an FBI (Fresh Blood Imaging) technique, has been proposed to overcome the foregoing inconveniences. In MR imaging on the FBI technique, an optimum time delayed from an R-wave of an ECG signal is predetermined, and ECG-synchronized MR scanning is performed at the delay time, thus well tracing a fresh and stable high-velocity blood flow ejected from the heart every appearance of the R-wave. In the FBI technique, three-dimensional scanning is additionally performed under the condition that signal intensities from parenchyma are actively suppressed by employing imaging conditions that include setting of a shorter repetition time TR (this causes the longitudinal relaxation time of parenchyma at rest to be insufficient) and applying an IR (Inversion Recovery) pulse or fat-suppression pulse (i.e., suppressing signals to be emanated from fat), thereby the blood flow being depicted. This eliminates the necessity of using a contrast medium, and blood flow images can be provided within a relatively shorter scan time.
For obtaining artery/vein visually separated blood flow images by using the FBI technique, a three-dimensional scan should be performed two times at different ECG-synchronized timings, and two sets of three-dimensional echo data acquired by the two-time three-dimensional scans or two sets of three-dimensional image data individually formed from the two sets of three-dimensional echo data should undergo weighted subtraction between the two sets of data.
In other words, even when the FBI technique is used, there remain some drawbacks that should be resolved. One drawback is that a longer scan time is still needed in total, because a three-dimensional scan should be performed two times. Another is that registration may be mistaken if the position of a patient's body moves between two times of scans, which is apt to deteriorate quality of blood flow images which will be produced by the subtraction.
On the other hand, the foregoing TOF and phase contrast techniques are based on the effect of flows of fluid such as blood. Although depending on the characteristics of an MRI system, it is general that either of the TOF or phase contrast method depicts only blood flows whose flowing speed is 2 to 3 cm/s or more. Blood flowing slower than this speed is scarcely detected. For example, peripheral veins, lymphatic vessels, CSF (cerebrospinal fluid), pancreatic duct, and others of a patient (human being) are slower in flow speed, and their flow speeds are approximately 1 cm/s or lower in general. Additionally, there may occur influence of positional shifts due to heartbeats, it was almost impossible to detect such slower-speed fluid flows by the conventional techniques.
The present invention has been made to break through the foregoing current situations. A first object of the present invention is to, therefore, provide an MR imaging technique for producing high-quality blood flow images in a shorter scan time, without using a contrast medium.
A second object of the present invention is to provide an MR imaging technique, in addition to the above first object, which is capable of obtaining different types of blood flow images from echo data acquired by the same scanning, thus enriching pieces of information to be provided about blood flows.
A third object of the present invention is to depict such slower-speed flows as peripheral blood flows in a steady manner, with no contrast medium injected.
A fourth object of the present invention is to depict such slower-speed flows as peripheral blood flows in a shorter period of time in a steady and high-quality manner, with no contrast medium injected.
In order to accomplish the above first and second objects, by an MRI system and an MR imaging method according to one aspect of the present invention, a plurality of different cardiac time phases of an object are set, an MR imaging scan is performed to start at the thus-set plural different time phases so that a plurality of sets of echo data are acquired successively, and a blood flow image is produced from the plurality of sets of echo data.
Preferably, the plural different time phases are two time phases falling into the systole and diastole of one cardiac cycle of the object. Still preferably, in the scan, a first scan which starts at the time phase present in the systole and a second scan which starts at the time phase present in the diastole are performed by separated pulse sequences toward the same slice of the object or the same slice encoding for the object.
Still preferably, echo data or image data thereof resultant from the first scan and echo data or image data thereof resultant from the second scan are subject to mutual subtraction, thereby producing echo data or image data thereof repenting an arterial phase image. For example, the subtraction is weighted subtraction.
Further, one example of setting time phases is directed to detection of a signal indicative of the cardiac time phases of the object. A preparing MR sequence is performed on a region to be imaged at each of different times from cyclically-appearing heartbeat reference waves of the detected signal, a plurality of times in total, so that a plurality of frames of MR images are obtained. From the plurality of frames of MR images thus obtained, two cardiac time phases, i.e., two timings in a cardiac cycle are determined. For example, the signal indicative of the cardiac time phase is an ECG signal of the object and the heartbeat reference wave is R-waves of the ECG signal.
This provides higher-quality blood flow images in a shorter scan time, without injection of a contrast medium. Moreover, from echo data acquired in the same imaging, blood flow images such as an arterial phase image and a venous phase image, which are different in types, can be produced in a simple manner. It is therefore possible to enrich blood flow information that can be provided through one time of imaging.
Additionally, in order to accomplish the foregoing third and fourth objects, another aspect of the present invention provides an MRI system and an MR imaging method, in which a scan is performed for an object placed in a static magnetic field, using a pulse sequence including a readout gradient pulse. For performing this scan, a cardiac time phase of the object is set, and the readout gradient pulse is applied to the object in a manner that its applied direction is substantially parallel to a flowing direction of blood in the object. In these states, the scan is performed in synchronism with the cardiac time phase that has been set, with echo signals acquired. An image of either a blood flow or a parenchymal region influenced by the blood flow is produced from the echo signals.
Preferably, the readout gradient pulse has a main pulse to read out an echo signal and a control pulse, which is added to the main pulse, to control behaviors of phase of magnetic spins present in blood. For instance, the control pulse is a pulse to dephase or rephase magnetic spins. By way of example, the cardiac time phase to be set is two in total, one for a systole and one for a diastole. At each of the two cardiac time phases, the scan is performed, so that data consisting of two sets of echo signals are acquired.
Since an applied direction of the readout gradient pulse is almost made to agree with a flowing direction of blood and a dephasing or rephasing pulse is added to the readout gradient pulse, a slower-speed flow, such as blood flows in an inferior limb, can be depicted with accuracy, without using a contrast medium. Particularly, a high-quality image in which arteries and veins are visually separated can be depicted within a short period of time.
In the accompanying drawings:
Hereinafter, embodiments of the present invention will be explained.
(1) First Embodiment
Refereeing to
(1.1) Configuration of System
The MRI system comprises a patient couch on which a patient P lies down, static magnetic field generating components for generating a static magnetic field, magnetic field gradient generating components for appending positional information to a static magnetic field, transmitting/receiving components for transmitting and receiving radio-frequency signals, control and operation components responsible for controlling the whole system and reconstructing images, and electrocardiogram components for acquiring an ECG signal of a patient, the ECG signal being employed as a signal indicative of cardiac time phases of the patient.
The static magnetic field generating components include a magnet 1 that is of, for example, a superconducting type and a static power supply 2 for supplying current to the magnet 1, and generates a static magnetic field H0 in an axial direction (Z-axis direction) in a cylindrical bore (diagnostic space) into which a patient P is inserted. The magnet unit includes shim coils 14. Current used to homogenize the static magnetic field is supplied from a shim coil power supply 15 to the shim coils 14 under the control of a host computer to be described later. The couch top of the patient couch on which the patient P lies down can be inserted into the bore of the magnet 1 so that the couch top can be withdrawn.
The magnetic field gradient generating components includes a gradient coil unit 3 incorporated in the magnet 1. The gradient coil unit 3 includes three pairs (kinds) of x-, y- and z-coils 3x to 3z used to generate magnetic field gradients changing in strength in X-axis, Y-axis, and Z-axis directions that are mutually orthogonal. The magnetic field gradient generating components further includes a gradient power supply 4 for supplying currents to the x-, y-, and z-coils 3x to 3z. The gradient power supply 4 supplies pulsated currents used to generate magnetic field gradients to the x-, y-, and z-coils 3x to 3z under the control of a sequencer that will be described later.
The pulsated currents supplied from the gradient power supply 4 to the x-, y-, and z-coils 3x to 3z are controlled, whereby magnetic field gradients changing in the three X-, Y-, and Z-directions (physical axis directions) are synthesized. Thus, a magnetic field gradient Os in a slice direction, a magnetic field gradient GE in a phase-encode direction, and a magnetic field gradient GR in a read-out direction (frequency-encoding direction), which are mutually orthogonal and logic axis directions, can be specified and changed arbitrarily. The gradients generated in the slice, phase-encode, and read-out directions are superposed on the static magnetic field H0.
The transmitting/receiving components includes an RF coil 7 located in the vicinity of a patient P in the diagnostic space inside the magnet 1, and a transmitter 8T and a receiver 8R both connected to the coil 7. The transmitter 8T and receiver 8R operate under the control of a sequencer 5 described later. The transmitter 8T supplies to the RF coil 7 RF current pulses of a Larmor frequency, which are used to induce the nuclear magnetic resonance (NMR). The receiver 8R takes in MR signals (radio-frequency signals) received by the RF coil 7, carries out various kinds of signal processing, such as pre-amplification, intermediate-frequency conversion, phase detection, low-frequency amplification, and filtering, on the echo signals, and carries out an A/D conversion on the processed echo signals so that digital data (original data) of the MR signals are produced.
Furthermore, the control and operation components include a sequencer 5 (also referred to as a sequence controller), host computer 6, calculation unit 10, storage unit 11, display unit 12, input device 13, and voice generator 16. Of them, the host computer 6 provides the sequencer 5 with pulse sequence information and manages the operation of the entire system according to not-shown installed software procedures
One feature of the MRI system is that it is able to perform an MR scan based on an electrocardiogram-synchronized technique depending on previously selected one or two synchronization timings (cardiac time phases). In the case that the synchronization timings are two in number, one is set to an optimum time phase residing in a diastole and the other is set to an optimum time phase residing in a systole, respectively.
In a main program not shown, the host computer 6 performs, as shown in
One execution routine of the ECG-prep scan is exemplified in
Optimum electrocardiogram synchronization timings are decided through the ECG-prep scan, before a scan for echo data acquisition is executed at the electrocardiogram synchronization timings. This permits blood flow to be traced in a steady manner and fresh blood outputted by the heart to be scanned at any time.
The sequencer 5, which has a CPU and memories, stores pulse-sequence information sent from the host computer 6, controls the operations of the gradient power supply 4, transmitter 8T, and receiver 8R according to the stored information, and temporarily receives digital data corresponding to MR signals outputted from the receiver 8R so as to transmit them to the calculation unit 10. The pulse-sequence information includes all information required for operating the gradient power supply 4, transmitter 8T, and receiver 8R according to a series of pulse sequences. For example, such information includes information about the strength, duration, and application timing of pulsed currents to be applied to the x-, y-, and z-coil 3x to 3z.
As to the pulse sequence, a two-dimensional (2D) scan or three-dimensional (3D) scan can be used, as long as a Fourier transform method is adopted. Particularly, the three-dimensional scan has a greater advantage in shortening a scan time. As pulse trains to those scans, various types of pulse trains based on a fast SE method, EPI (Echo Planar Imaging) method, FASE (Fast Asymmetric SE) method (that is, an imaging technique in which both of the fast SE and half-Fourier methods are combined), and others are available.
The calculation unit 10 receives digital data (also known as original data or raw data) sent from the receiver 8R via the sequencer 5, maps the original data in a Fourier space (also known as a k-space or frequency space) formed in its incorporated memory, and reconstructs the mapped original data into an image in the real space through a two-dimensional or three-dimensional Fourier transform for each set of data. Moreover, the calculation unit performs synthesis and subtraction (weighted subtraction is included) with data of images, according to its necessity. The synthesis includes pixel-by-pixel addition of image data of a plurality of frames and maximum intensity projection (MIP) processing of a plurality of frames of images. Another example of the synthesis is a method by which original data of a plurality of frames are synthesized into a single frame of original data, as they are, with the axes of the frames matched in the Fourier space. Additionally, the addition includes simple addition, averaging, or weighted addition.
The storage unit 11 is able to preserve image data produced by the synthesis or subtraction as well as the reconstructed image data. The display unit 12 displays an image. By using the input device 13, an operator is able to provide with the host computer 6 parameter information for selecting desired synchronization timing, scan conditions, a pulse sequence, and information about processing image synthesis and subtraction.
The voice generator 16 is capable of uttering voice messages informing a patient of the start and end of breath hold in response to instructions sent from the host computer 6.
Furthermore, the electrocardiogram components comprise an ECG sensor 17 attached to a patient body to detect an electric ECG signal and an ECG unit 18 performing various types of processing including digitization with the detected ECG signal and sending it to both the host computer 6 and the sequencer 5. The sequencer 5 uses this measured ECG signal when performing each of the ECG-prep scan and the ECG-synchronized imaging scan. This enables optimum setting of synchronization timing based on the ECG-synchronized method, and data acquisition can be done by the ECG-synchronized imaging scan on the basis of the set synchronization timing.
(1.2) ECG-prep Scan
Referring to
The host computer 6, which is in operation for a given main program not shown, responds to a command from the input device 13 and commences to execute an ECG-prep scan shown in
First, the host computer 6 reads from the input device 13 scan conditions and information about parameters both required to perform the ECG-prep scan (step S1 in the figure). The scan conditions include the type of a scan, the type of a pulse sequence, and a phase-encode direction. The parameter information includes an initial time T0 (herein, defined as an elapsing time from an R-wave peak in the ECG signal) to determine an ECG-synchronized timing (time phase), a time increment Δt, and an upper limit of a numbering counter CNT. An operator can properly set these parameters.
The initial time T0, time increment Δt, and the upper limit of the numbering counter CNT are set to amounts so that, for example, a range from a diastole to a systole in a period of “1 R-R” is almost thoroughly covered in time. In the diastole, arterial and venous phases are both depicted, while in the systole, only a venous phase is depicted. By way of example, the initial time T0 may be set to T0=0.
The host computer 6 then initializes the numbering counter CNT counting the execution times of the sequence and a time increment parameter Tinc determining the synchronization timing (CNT=0, Tinc=0; step S2). After this, the host computer 6 sends massage data to the voice generator 16 to generate breath-hold instructions, such as “Hold your breath, please.” toward an object (patient) (step S3). It is preferred that the breath hold is performed for suppression of body motions of a patient that may be caused during the ECG-prep scan. However, in some occasions, the ECG-prep scan may be performed with no breath hold
After having completed the above preparation, the host computer 6 sequentially executes processes shown after step 4. This execution permits the scan with the ECG-synchronized timing changed.
Specifically, an expression of TDL=T0+Tinc is calculated to obtain a delay time TDL from the peak time instant of an R-wave (step S4) An ECG signal that has experienced the signal processing in the ECO unit 18 is then read, and it is determined whether or not the R-wave peak value has appeared in the signal (step S5). This determination will be repeated until the R-wave appears. When the R-wave appears (Yes at step S5), it is then determined whether or not the delay time TDL calculated at step S4 has elapsed since the appearance of the R-wave peak time (step S6). This determination will also be repeated until the delay time TDL elapses.
When the time has passed by the delay time TDL since the R-wave peak time instant (Yes at step S6), the sequencer 5 is ordered to start a pulse sequence of each time (step S7, refer to
After the above sequence has been instructed to start, the numbering counter CNT is incremented such that CNT=CNT+1 (step S8), then the time increment parameter Tinc is computed such that Tinc=ΔT×CNT (step S9). In other words, every time when the pulse sequence is ordered to be executed, the count of the counter CNT increases by one and the increment parameter Tinc for adjusting the synchronization timing increases in proportion to the count.
Then, a standby state continues until a period of predetermined time (for example, approx. 500 to 1000 msec) necessary for the execution of the pulse sequence of each time passes (step S10). Then, whether the count of the numbering counter CNT reaches the preset upper limit or not is determined (step S11). In cases where, for example, five two-dimensional images are produced with the delay time TDL changed into various amounts for the purpose of optimizing the synchronization timing, the count in the counter CNT is set to “5.” If the count has not yet reached the upper limit (No at step S11), the processing is returned to step S5 to repeat the above processing. In contrast, the count of the counter CNT equals the upper limit (Yes at step S11), a command to release the breath hold is sent to the voice generator 16 (step S12), and the processing returns to the main program. A voice message to release the patient from the breath hold is such that “you may breathe.”
Executing the above processes sequentially leads to the execution of the preparing pulse sequence of which timing is exemplified in
Therefore, when the first R-wave peak appears after the instructions of the breath hold, the first scan IMGprep1 based on, for example, a two-dimensional FASE method is executed so that it lasts for a certain period of time (for example, approx. 500 to 1000 msec) from a certain time instant when the delay time TDL(=T0) has elapsed after the R-wave peak appearance, whereby echo signals being acquired. Whenever the next R-wave may appear during the continuation of this sequence, the foregoing waiting process at step S10 in
Unless the count of the numbering counter CNT has yet reached its limit, steps S5 to S11 will be executed again. Thus, in the embodiment as shown in
The echo data are sent to the calculation unit 10 via the receiver 8R and then the sequencer 5 in turn. The calculation unit 10 reconstructs image data mapped in the k-space (frequency space) into image data in the real space by means of a two-dimensional Fourier transform. The reconstructed image data are stored in the storage unit 11 as blood flow image data. The host computer 6 responds to, for example, operation signals from the input device 13 so that images of blood flow are sequentially displayed in a dynamic (CINE) mode.
As pictorially shown in
As described above, from a plurality of blood flow images scanned as the delay time TDL was changed dynamically, the operator decides delay times TDL (for example, two delay times TDL1 and TDL2) serving as optimum synchronization timings for the systole and diastole by visual observation. And the operator carries out a command, for example, by hand, for reflecting the decided delay times TDL into an imaging scan which will follow.
As to reflecting the decided delay times into the imaging scan, a further configuration can also be realized by using software. The configuration is that when images that have been determined through visual observation are specified, delay times TDL assigned to the specified images are automatically memorized as optimum synchronization timings, and the timings are automatically read out in performing the imaging scan. This makes it possible to specify ECG-synchronized timings in an automatic fashion.
By the way, in the foregoing ECG-prep scan, the phase-encode direction is positively made to agree with the running direction (i.e., body-axis direction) of a blood flow such as the aorta. Compared to cases where the phase-encode direction is set to other directions, this setting of the phase-encode direction leads to clear images in which information about blood flow directions (directional performance) is avoided from being dropped, providing its superior depiction capability.
(1.3) Imaging Scan
Referring to
The host computer 6 executes the processes shown in
Specifically, first, the host computer 6 reads from, for example, the input device 13 the two optimum delay times TDL determined by the operator through the foregoing ECG-prep scan (step S20). The delay times TDL are an optimum delay time TDL1 given to a systole and an optimum delay time TDL2(>TDL1) given to a diastole. Information in relation to the optimum delay times TDL1 and TDL2 may be stored in advance in the storage unit 11.
Then, the host computer 6 inputs information about scan conditions, image processing techniques, and others, which are specified by the operator using the input device 13, processes the information including the delay times TDL1 and TDL2 into control data, and outputs the control data to both sequencer 5 and operation unit 10 according to their necessity (step S21). The scan conditions include a phase-encode direction, an image size, the number of scans, a waiting time between scans, and a pulse sequence depending on a region to be scanned. The image processing techniques include a subtraction method and its weighting factors, an addition method (simple, averaging, or weighted addition method), and/or a maximum intensity projection (MIP) method.
When it is determined that an instruction indicating the completion of scan preparations has been issued (step S22), a command indicating the start of a breath hold is output to the voice generator 14 (step S23). This causes the voice generator 14 to utter a voice message saying, “Hold your breath, please.” like the ECG-prep scan. In response to this message, a patient is to hold breathing (refer to
After this, the host computer 6 instructs the sequencer 5 to start the imaging scan (step S24).
When having received instructions to start the imaging scan (step S24-1 in
When the specified n-th R-wave has appeared, processing to wait for the delay time TDL1 determined for a specified time phase in the systole is executed first (step S24-4). The delay time TDL1 is, as explained before, optimized through the ECG-prep scan such that echo signal intensities become the highest in imaging objective venous flows in a systole, providing a superior depiction capability of the entity.
The sequencer 5 begins to perform the imaging scan for a systole at a time when this optimum delay time TDL1 has passed, the time being regarded as an optimum ECG-synchronized timing (step S24-5). Practically, the transmitter 8T and the gradient power supply 4 are driven based on the pulse sequence information memorized beforehand. Thus a scan (the first scan) SNsys1 is performed based on the first slice-encode amount SE1 defined by, for example, a three-dimensional FASE pulse sequence according to the ECG-synchronized technique, as shown in
In the pulse sequence used for the scan SNsysn in the systole, the number of echoes is decreased as seen in
Thus, under the first slice-encode amount SE1, echo signals are acquired during a short scan time of about several hundreds msec from a three-dimensional imaging region Rima set to, for example, the hypogastrium as shown in the
The sequencer 5 then proceeds to scan control for the diastole. Specifically, processing to wait for the delay time TDL2 that is determined for a specific time phase in the diastole is performed (step S24-6). As described before, the ECG-prep scan permits the delay time TDL2 to be optimized to an amount that produces echo signals into the highest intensity in imaging targeted arterial and venous flows in a diastole, providing a superior depiction capability of the entities.
The sequencer 5 begins to perform the imaging scan for the diastole at a time when this optimum delay time TDL2 has passed, the time being regarded as an optimum ECG-synchronized timing (step S24-7). Practically, the transmitter 8T and the gradient power supply 4 are driven based on the pulse sequence information memorized beforehand. Thus a scan SNsys2 is performed based on the first slice-encode amount SE1 defined by, for example, a three-dimensional FASE pulse sequence according to the ECG-synchronized technique, as shown in
The pulse sequence used for the scan SNdian for the diastole is determined, as shown in
Thus, under the first slice-encode amount SE1, echo signals are acquired during a scan time of about 600 msec from the three-dimensional imaging region Rima set to the hypogastrium as shown in the
On having completed the first imaging scan, the sequencer 5 determines if the final imaging scan has been completed or not (step S24-8). In the case of NO at this determination (the final scan has not been completed yet), with monitoring the ECG signal, waiting is done until a shortly set period of time (for example, 2 heartbeats (2R-R) from the R-wave used in the imaging scan) passes. This results in that the recovery of longitudinal magnetization of spins in the stationary parenchyma is actively suppressed (step S24-9).
After waiting for a period of time corresponding to, for example, 2R-R, when the third R-wave appears (YES at step S2497), the sequencer 5 returns its processing to the foregoing step S24-4.
Hence, at a time instant when the specified delay time TDL1 has passed since the third R-wave peak, a second scan SNsys2 for the systole is commenced in the same way as the above under the next slice-encode amount SE2. Echo signals are therefore acquired from the three-dimensional imaging region Rima (steps S24-4 and S24-5). Then, at a time instant when the specified delay time TDL2 has passed since the third R-wave peak, a second scan SNdia2 for the diastole is commenced in the same way as the above under the next slice-encode amount SE2. Echo signals are therefore acquired from the three-dimensional imaging region Rima (steps S24-6, 7).
Likewise, echo signals are acquired for each of the systole and diastole until the final slice-encode amount SEn (for example n=8).
On having completed the final scan SNsysn and SNdian under the slice-encode amount SEn, the determination at step S24-8 becomes YES, thus a notification stating the completion of the imaging scan is sent from the sequencer 5 to the host computer 6 (step S24-10). Accordingly the processing is returned to the host computer 6.
When receiving the notification stating the completion of the imaging scan from the sequencer 5 (
Therefore, as pictorially shown in
Echo signals emanated from the patient P are received scan by scan by the RF coil 7, then sent to the receiver 8R The receiver 8R processes the echo signals with various kinds of preprocessing to convert them into digital quantities. The digital echo data are sent via to the sequencer 5 to the calculation unit 10, where they are mapped in each of two three-dimensional k-spaces Ksys and Kdia formed by memories, correspondingly to the phase-encode amounts and slice-encode amounts.
(1.4) Data Processing and Image Display
After the echo data acquisition, the host computer 6 instructs the calculation unit 10 to execute the processing shown in
As shown in
The calculations unit 10 then reconstructs images by performing a three-dimensional Fourier transform in relation to the k-space Ksys for systole and k-space Kdia for diastole, respectively (steps S33 and S34). Accordingly, as shown in
Hence, in order to obtain an arterial-phase image IAR, the calculation unit 10 performs the subtraction of “IMdia−IMsys” pixel by pixel (step S35). Here, β is a weighting factor. By designating an appropriate amount to the weighting factor β, this subtraction allows image data of the vein VE to be nearly diminished to zero, thus providing three-dimensional image data of an arterial-phase image IMAR containing only the artery AR, as pictorially shown in
Further, to obtain a venous-phase image IMVE, subtraction of “IMdia−IMAR” is performed on the pixel basis (step S36). The image data IMAR used in this calculation have been derived by the weighted-subtraction described above. This causes image data of the artery AR to reduce to zero, thereby providing three-dimensional image data of a venous-phase image IMVE containing only the vein VE. This subtraction can also be replaced by weighted-subtraction.
After such subtraction, the calculation unit 10 performs MIP (maximum intensity projection) processing with each of the arterial-phase image IMAR and venous-phase image IMVE. This produces data of a two-dimensional image (for example, coronal image) viewed along a desired observing direction of blood vessels (step S37). These two-dimensional images of the arterial- and venous-phases are displayed on the display unit 12 as shown in
(1.5) Operations and Advantages
As descried above, in the MRI system of the present embodiment, appropriate scan start timings (i.e., delay times from the R-wave) are determined for each of the systole and diastole. Then, at the timings, the scans of two shots for the systole and diastole are individually and sequentially performed under each slice-encode amount. In addition, the systolic-use scan that always precedes in each cardiac cycle is shorter in data acquisition time (the number of echoes) so that it does not overhang the following diastolic-use scan in time. Echo data acquired by the systolic-use scan are mapped in a lower-frequency region of the systolic-use k-space that is most significant in improving contrast. The remaining data in the systolic-use k-space, which have not been acquired, are duplicated from part of data acquired by the following diastolic-use san that is allowed to acquire echoes during a longer period of time. Moreover, the scans for the systole and diastole adopt the half-Fourier technique in order to set the scan time as short as possible.
Accordingly, the two-shots of scans for the systole and diastole, which are performed under one slice-encode amount, usually remain within an interval of two heartbeats. Consecutively repeating such scans enables a three-dimensional scan such that echo data of each of systolic and diastolic blood flows are acquired at appropriate timings during one time of breath hold duration in three-dimensional scanning. In other words, three-dimensional image data of blood flows in each of the systole and diastole are obtained by one time of imaging performed at proper timings. The acquired data are then subjected to the reconstruction and subtraction described before, thereby arterial-phase and venous-phase images being provided.
In the present embodiment, it is therefore unnecessary to perform imaging scans (i.e., in total, two times of scans) for the systole and diastole in a separated fashion, and one time of scan is enough. Therefore, the scan time is reduced largely and throughput of patients is improved. Particularly, an effect of shortening the san time becomes noticeable in three-dimensional scanning. Furthermore, misregistration due to patient's body motions can be lessened markedly, thereby images of higher quality being provided. Further, using echo data acquired at two cardiac phases through one time of imaging, blood flow images (MRA images) in which arterial and venous phases are visually separated can be obtained. Thus, the imaging can be done with improved efficiency, while blood flow information supplied can be enriched.
Because most optimum ECG-synchronized timings for the systole and diastole are determined in advance through the ECG-prep scan, it is possible that target blood flows are traced in a steady manner at each phase during the systole and diastole. This provides blood flow images of which signal intensity is higher, blood flow contrast is improved, and S/N is superior. On the other hand, previously setting appropriate ECG-synchronized timings eliminates the need of re-performing scans in most cases, thus relieving operational burdens on operators and physical and mental loads on patients.
Further, in addition to shortening the repetition time TR and echo train spacing, the phase-encode direction can be made to approximately agree with a blood-flow direction and the slice direction can be set along the front/rear direction of a patient. Compared to the TOF technique that requires slices to be set perpendicular to a blood flow, the entire scan time can therefore be shortened. Further, since the number of applications of the slice encode is reduced by an amount corresponding to a shortened scan range in the slice direction, the whole scan time is reduced largely, compared to the conventional TOF or phase encode technique. This will lead to a reduced burden on patients and an improved throughput of patients.
Additionally, because it is unnecessary to inject a contrast medium into an object, non-invasive imaging can be provided. This will also result in largely lessened mental and physical burdens on a patient. Cumbersome operations inherent to the contrast technique, such as paying attention to a contrast effect of the contrast medium, are not required as well. Thanks to those advantages, differently from the contrast technique, the imaging can be performed repeatedly if necessary.
Further, because the phase-encode direction is made to agree or nearly agree with a running direction of vessels, blurs of pixels can be utilized positively. This provides a remarkable depiction capability in the running direction of vessels. Further, by changing the phase-encode direction according to the vessel-running direction in a portion to be imaged, various portions of an object can be imaged with ease.
Because the present embodiment uses the pulse sequence based on the fast SE technique, the imaging is advantageous in susceptibility and contour distortion,
(1.6) Modifications of First Embodiment
The present invention is not limited to the configuration described in the above embodiment, but it can be modified into various ways and practiced into various applications.
For example, in the foregoing embodiment, the configuration has been made such that both arterial-phase and venous-phase images are present. In this respect, only an arterial-phase image can be produced by the subtraction and displayed. By this modification, the subtraction step S36 in
In the foregoing embodiment, each of the scans for the systole and diastole has been conducted by scanning based on the half-Fourier technique, but the scanning will not always be confined to a technique based on the half-Fourier technique. In such a case, the scan for the diastole is performed to acquire data mapped entirely into the k-space, and echo data mapped in its high-frequency regions positioned at both end sides in the slice-encode direction are duplicated to corresponding regions in the k-space for the systole.
The three-dimensional scanning has been employed by the foregoing embodiment, but the two-dimensional scanning can be used as well. Pulse sequences that are possible to be employed are not limited to the FASE method, but other pulse sequences derived from the FSE or EPI method may be used as well.
Moreover, the post-processing of echo data in the foregoing embodiment has been configured such that the echo data are once converted to image data in the real space, then the subtraction is performed to obtain the arterial-phase and venous-phase images. Alternatively, the subtraction may be done such that it is performed at the stage of echo data on the k-spaces Ksys and Kdia of which matrix sizes are the same to each other. Echo data resulted from the subtraction then undergo a reconstruction process to provide blood flow images.
The technique to obtain arterial-phase and/or venous-phase images is not limited to the subtraction between data acquired at the two different cardiac phases described in the embodiment. Alternatively, there may be provided a technique of performing subtraction between images of which echo train spacing are different from each other, or a technique of performing subtraction between images of which effective TE times are different from each other. Differences in the echo train spacing give changes to sensitivity in detecting speed of blood flows. This allows acquisition of echo data into which differences in blood flow speed due to inherence of the artery and vein are reflected. Therefore, the subtraction that will be performed in the similar way to the above can provide blood flow images of the artery and vein. In addition, differences in the effective TE time enables echo data to be acquired as the artery and vein of which T2 time differs from each other are differentiated. Thus the similar subtraction to the above can provide blood flow images of the artery and vein.
(2) Second Embodiment
Referring to the foregoing drawings and
MR imaging according to the second embodiment is characterized in that a dephasing or rephasing pulse is added to the read-out gradient pulse GR in order to depict slow-speed blood flow such as blood flow in the inferior limb. The present embodiment uses an MRI system that is the same or identical in both hardware configuration and ECG-prep scan as or to those in the first embodiment.
(2.1) Imaging Scan
In the present embodiment, as shown in
Referring to
Now suppose that the first time of imaging scan (imaging) is assigned to the systole. In this case, the host computer 6 first read, for example, from the input device 13, an optimum delay time TDL (=TDL1 or TDL2>TDL1) dedicated to the systole, which is determined by an operator through the foregoing ECG-prep scan (step S70).
The host computer 6 then inputs scan conditions and information about image processing that the operator specified from the input device 13, before processes those inputs into control data to be outputted to both of the sequencer 5 and calculation unit 10 (step S71). The scan conditions include an applied direction of the read-out gradient pulse, an image size, the number of times of scanning, an interval between scans, a pulse sequence according to a portion to be scanned, and others. The image processing information includes information indicative of MIP processing and/or subtraction. In the case of subtraction, the information shows that the subtraction is simple subtraction, weighted subtraction, or addition. Further, the control data include the delay time TDL.
Like the first embodiment, when it is determined that the preparation for scanning has been notified (step S72), a command of starting a breath hold (step S73), a command of starting the scans (step S74), a determination whether the scans have been completed or not (step S75), a command of releasing the breath hold (step S76), and instructions to image processing and display (step S76) are performed in turn.
At step S74, to be specific, the host computer 6 instructs the sequencer 5 to start the first-(or second-)time imaging scan.
On receiving such instructions to start the imaging scan (step S74-1,
When the optimum delay time TDL1 (or TDL2) has passed, it is regarded that an optimum ECG-synchronized timing has come. The sequencer 5 therefore performs the first-time imaging scan (step S74-5). Specifically, the sequencer 5 drives the transmitter 8T and gradient power supply 4 according to information concerning a pulse sequence stored already. As a result, the first-time imaging scan (i.e., imaging) based on a pulse sequence using the three-dimensional FASE technique is executed with the ECG-synchronized technique, as shown in
In this imaging on this pulse sequence, the read-out gradient pulse GR is applied in a direction RO, as shown in
The read-out gradient pulse GR included in the pulse sequence is provided, as shown in
In contrast, the dephasing pulses Pdephase gives few dephasing effect to stationary or almost stationary magnetic spine. It is therefore significant that the read-out gradient pulse GR be applied substantially in a direction along which fluid (blood and lymph) to be imaged moves.
Preferably, the dephasing pulses Pdephase are changeable or controllable in its intensity in accordance with flow speeds of lymph or blood, which is fluid to be imaged.
When fluid (such as blood flow) to be imaged is relatively high in speed, a total of two rephasing pulses Prephase are continuously added to the temporal forward and backward ends of the pulse body Pbody, as shown in
In this first embodiment, the first and second (later described) imaging scans employ the read-out gradient pulse GR to which either dephasing pulse Pdephase or rephasing pulse Prephase is added.
Executing the three-dimensional FASE pulse sequence makes it possible that echo signals stimulated by both of the excitation 90-degrees RF pulse and the refocusing 180-degrees RP pulses are acquired at each phase-encode amount assigned to each slice-encode amount. In each echo signal, dephasing of phase of magnetic spins caused by the dephasing pulse Pdephase or rephasing of phase of magnetic spins caused by the rephasing pulse Prephase is reflected.
This reflection will also be detailed together with a display operation described later, but the summary thereof is as follows.
To fluid that is flowing along an applied direction of the read-out gradient pulse, a dephasing effect caused by the dephasing pulse Pdephase promotes the flow void effect. Thus the dephasing pulse reduces the intensity of an echo signal. In the case that the fluid hardly flow along such direction, only a smaller amount of promotion of the flow void effect is available due to the dephasing pulse Pdephase, thereby the intensity of an echo signal being not so much reduced.
In the case of the rephasing pulse Prephase, its rephasing effect is able to control degrees of dephasing according to flow speeds of fluid.
Echo train spacing in the foregoing pulse sequence is shortened to an amount as small as approximately 5 msec. Under the first slice-encode amount SE1, echo signals are acquired from a three-dimensional imaging region Rima directed to, for example, the inferior limb as shown in
When the scanning has been completed at the first slice encoding, the sequencer 5 determines whether or not the scanning under the final slice encoding has been completed (step S74-6, in
In this way, waiting for a period of time corresponding to, for example, 2 R-R, is done until the appearance of, for example, the third R wave, and when such R wave appears (YES at step S74-7), the sequencer 5 returns the processing to the foregoing step S74-4. Accordingly, at a time when a specified delay time TDL1 has passed from the reception of an ECG trigger signal synchronized with the peak value timing of the third R wave, scanning at the next slice encode amount SE2 starts in a similar manner to the foregoing one, so that echo signals are acquired from the three-dimensional imaging region Rima again (steps S74-4 and S74-5). Hereafter, the processing will be continued until echo signals at the last slice encode amount SEn (for instance, n=8) are acquired.
When the last scanning at the slice-encode amount SEn has been completed, the determination at step S74-6 becomes YES, the sequencer 5 informs the host computer 6 of the completion of the first-time (or second-time) imaging scan (steps S74-8). Then the processing is handed to the host computer 6.
According to the above procedures, the first-time (or second-time) imaging scan (imaging) employing the ECG synchronization technique is performed every a period of time of 2 R-R on the basis of, for example, 3D-FASE method.
Each echo signal emanated from the patient P is received by the RF coil 7 and sent to the receiver 8R, for each slice encode amount supplied by the slice gradient pulse GS. The receiver 8R performs various types of pre-processing on the echo signal and converts it into digital amount of data. The digital echo data thus produced are sent to the calculation unit 10 through the sequencer 5, and mapped at given positions in a three-dimensional k-space formed by a memory, accordingly to the encoded amounts given to the echo signal.
Then, as shown in
As a result, for this second-time imaging scan, as shown in
As a result, the second-time imaging scan is able to provide image data in the diastole, which are influenced by the spin control of either dephasing pulses Pdephase or rephasing pulses Prephase added to the read-out gradient pulse GR, similarly to the first-time imaging scan.
(2.2) Data Processing and Image Display
When the echo data acquisition has been finished, the host computer 6 obliges the calculation unit 10 to execute the processing shown in
As shown therein, responsively to instructions from the host computer 6, the calculation unit 6 calculates echo data using the half Fourier technique from both of the systole-use k-space and the diastole-use k-space (step S81). That is, echo data that should be mapped in a remaining region of each k-space, but has been left with no data acquisition are calculated from the complex conjugate relationship, and mapped therein. This calculation completely fills up both k-spaces with echo data.
The calculation unit 10 reconstructs echo data in each of the k-paces for the systole and the diastole into image data through a three-dimensional Fourier transform, space by space (steps S82 and S83). Like the foregoing
A vein VE is only reflected in the systolic image IMsys, but image data of an artery VR is hardly included in the systolic image IMsys. On the other hand, in the diastolic image IMdia, both of an artery AR and a vein VE are reflected, though degrees of the reflection are different from each other.
The principle of obtaining such systolic image IMsys and diastolic image IMdia will now be detailed in terms of the applied direction of the foregoing read-out gradient pulse GR and the function of the dephasing pulses Pdephase.
The phase of magnetic spins of an object, such as blood, that flows in the applied direction of the read-out gradient pulse makes it easier to be dephased more quickly on account of the dephasing pulses applied. In other words, to an object that is flowing, this is equivalent to the fact that the flow void effect provided from the flow itself is promoted. In contrast, the rephasing pulses give a rephasing function to the phase of magnetic spins of such blood flow.
For example, the inferior limb of an object to be examined will now be exemplified. In the case of the inferior limb, even if the artery is measured in the systole, its flow speed is slow and normally less than 1 cm/sec. Moreover, as for the vein measured in the systole and the artery and vein measured in the diastole, the blood moves at extremely slow speeds that can be regarded as if the blood is stationary. As shown in
Magnetic spins of the artery and vein are excited by those imaging scans to acquire echo signals. In this acquisition, flow speeds of the artery and vein differ from each other, although the difference might be rather small. Hence the difference in flow speed is reflected into promotion of a flow void effect based on the rephasing dephasing pulses, providing relative changes between intensities of echo signals.
To be specific, the systole will now be explained as follows. Since flowing at extremely slow speeds, the vein, when observed during the systole, is less in the flow void effect and depicted as bright blood with relative higher signal intensities, through though it suffers a slight decrease in echo signal intensity due to the dephasing pulses. By contrast, the artery, when observed during the diastole systole, flows at larger speeds than those of the vein, so that the promotion of the flow void effect caused by the dephasing pulses is larger than that of the vein. This causes a larger decrease in the signal intensity of the artery, which depicts the artery as black blood. This state can be pictorially shown in a similar way to
On the other hand, since both artery and vein only move at extremely lower speeds during the diastole, they are depicted as bright blood, though they experience slight reductions in signal intensity because of the dephasing pulses. This condition is pictorially shown in a similar manner to
Backing to the explanation of
Further, to obtain a venous phase image IMVE, a subtraction of “IMdia−IMAR” is performed pixel by pixel (step S85). The image data IMAR has already been calculated by the foregoing weighed subtraction. In a similar way to
After the above subtractions, the calculation unit 10 proceeds to perform MIP (maximum intensity projection) processing for each of the arterial phase image IMAR and the venous phase image IMVE. This produces data of a two-dimensional image (e.g., coronal image) obtained by observing, along a desired direction, blood vessels residing in each of the images (step S86).
The two-dimensional images IMAR and IMVE for the arterial and venous phases are displayed on the display unit 12 as shown in
In addition to displaying the arterial and venous images IMAR and IMVE, the systolic and diastolic images IMsys and IMdia may be displayed on the same screen to those for the arterial and venous images or on the screens of different monitors from the arterial and venous images.
(2.3) Advantages
As described, the MRI system of this embodiment employs the imaging in which the applied direction of the read-out gradient pulse GR is made to almost agree with a flow direction of fluid (such as blood) of which flow speed is lower, as can be observed in the inferior limb. Concurrently, the dephasing pulses Pdephase or rephasing pulses Prephase added to the read-out gradient pulse GR.
Thus, the dephasing pulses Pdephase or rephasing pulses Prephase are able to enhance relative differences of signal intensity between a first fluid that flows and a second fluid that flows at a slower speed than the first fluid. Therefore, even if blood vessels in the inferior limb, which are slower in flow speed than the abdomen and thorax, are imaged using, for example, the dephasing pulses, the relative differences of signal intensity are able to provide an image as shown in
According the inventors' recognition, the above technique that the read-out gradient pulse is applied in the substantially same direction as the flow direction of fluid and the flow void effect is controlled by the positive use of dephasing and rephasing of magnetic spins has been newly developed. This technique can give relative differences to signal intensity between the artery and vein.
Additionally, in this embodiment, because the ECG-prep scan is used to previously determine the optimum ECG-synchronized timing for the systole and the diastole, blood flows targeted at each time phase during each of the systole and the diastole can be traced without fail. Previously conducted appropriate setting of the ECO-synchronized timing eliminates the necessity of repeating the same imaging. Operational work on operators and physical and mental burdens on patients are therefore reduced largely.
Further, it is possible to specify the slice or slice-encode direction in directions other than the superior-inferior direction of a patient, the entire scan time can be shortened, compared to imaging methods, such as the TOF technique, that require scanning to advance in the superior-inferior direction. This also lowers patient's burdens and increases throughput of patients.
It is unnecessary to inject a contrast medium into a patient, which leads to non-invasive imaging. This also remarkably reduces physical and mental burdens on patients. Also, troublesome operations inherent to the contrast technique, such as timing of a contrast effect should be measured, is also unnecessary. Differently from the contrast technique, the imaging technique according to the present embodiment can be repeated if necessary.
(Modifications of Second Embodiment)
The above embodiment uses both of the first-time and second-time imaging cans involving the read-out gradient pulse GR to which either of the dephasing pulses Pdephase or the rephasing pulses Prephase are added (refer to
As a modification of this embodiment, the dephasing pulses Pdephase may be added in the first-time imaging scan conducted at a time phase during the diastole systole, as shown in
In other words, between the systole and the diastole, the type of the control pulses to additionally control behaviors of magnetic spins is changed. This makes it possible to reflect the more effect of rephasing (i.e., flow compensation) in signal intensity in the diastole, thus increasing the signal intensity to improve a signal-to-noise ratio.
(Third Embodiment)
Referring to the foregoing figures and
In the third embodiment, the first-time and two-time imaging scans, that is, two times of imaging scans which have been conducted in the second embodiment are conducted as one-time imaging scan. In this scan, the foregoing dephasing and rephasing pulses are used according to the systole and diastole in each cardiac cycle.
The artery and vein in the inferior limb will now be employed as fluid of a slower speed and an artery/vein visually separated image thereof will now be obtained. Similarly to the sequence shown in
Then, the imaging scan is conducted in the form of a one-time imaging scan on the basis of the ECG-synchronized technique involving delay times TDL1 and TDL2. The procedures of this imaging scan, which are similar to those in
(3.1) Imaging Scan
During performance of a not-shown main program, the host computer 6 also performs the processing shown in
Specifically, first of all, the host computer 6 reads two optimum delay times TDL via the input device 13, for example (step S120). The delay times TDL, which are previously determined through the foregoing ECG-prep scan by an operator, are composed of an optimum delay time TDL1 for the systole and an optimum delay time TDL2 (>TDL1) for the diastole, as described above. Information about those optimum delay times TDL1 and TDL2 may previously be determined and stored in, for example, the storage unit 11.
Then, the host computer 6 inputs information about scan conditions, an image processing method, and others, and process the information including the delay times TDL1 and TDL2 into control data. The control data are outputted to both of the sequencer 5 and the calculation unit 10 according to necessity (step S121).
It is then determined by the host computer 6 if the preparation before scanning has been completed or not in the similar manner to that in first embodiment. When completion of the preparation is determined, breath hold is instructed, then the imaging scan is instructed to start (steps S123 and S124).
On receiving instructions of starting the imaging scan (step S124-1 at
When the appearance of the n-th R-wave is realized, the sequencer 5 waits for the delay time TDL1 set to a specific time phase in the systole (step S124-4).
A time when the optimum delay time TDL1 has passed is considered to be an optimum ECG-synchronized timing. Hence, the sequencer 5 begins to execute scanning for the systole at that time (step S124-5).
Specifically, according to information in relation to a pulse sequence memorized in advance, the transmitter 8T and gradient power supply 4 are driven. By this drive, a first scan SNsys1 is performed based on the ECG-synchronized technique as shown in
In this first san SNsys1 the read-out gradient pulse GR applied to the patient's body axis direction substantially in parallel with the artery and vein in the patient's inferior limb. Additionally, dephasing pulses Pdephase to dephase the phases of magnetic spins added to the temporal forward and backward parts of the read-out gradient pulse GR without temporal gaps. The echo train spacing used in this pulse sequence is shortened to approximately 5 msec.
The pulse sequence used for the first scan SNsysn assigned to the systole adopts a less number of echoes that consecutively continue only during a shorter period of time after the start of the scan within one heartbeat, as shown in
Therefore, at the first slice-encode amount SE1, echo signals are acquired from a three-dimensional imaging region Rima (refer to
The sequencer 5 then proceeds to scan control for the diastole. Specifically, the sequencer 5 waits for the delay time TDL2 set to a specific time phase during the diastole (step S124-6).
A time when the optimum delay time TDL2 has passed is considered to be an optimum ECG-synchronized timing. Hence, the sequencer 5 executes a second scan for the diastole (step S124-7). Specifically, according to information in relation to a pulse sequence memorized in advance, the transmitter 8T and gradient power supply 4 are driven. By this drive, the first scan SNdia1 is performed based on the ECG-synchronized technique as shown in
In this second san SNdia1, the read-out gradient pulse GR is also applied to the patient's body axis direction substantially in parallel with the artery and vein in the patient's inferior limb. Additionally, rephasing pulses Prephase to rephase the phases of magnetic spins added as shown to the temporal forward and backward parts of the read-out gradient pulse GR without temporal gaps. The echo train spacing used in this pulse sequence is also shortened to approximately 5 msec.
The pulse sequence used for the second scan SNdian assigned to the diastole is set, as shown in
Therefore, at the first slice-encode amount SE1, echo signals are acquired from the three-dimensional imaging region Rima (refer to
When the first-time imaging scan is completed, the sequencer 5 determines whether or not the last scan has been completed (step S124-8). If determined to be NO (the last scan has not been ended yet), waiting will be continued, with the ECG signal monitored, until a predetermined shorter interval of time pass. This waiting permits the longitudinal magnetization of spins in stationary parenchyma to be positively suppressed from being restored (step S124-9). Such shorter interval of time for waiting is, for example, “2 R-R” from the R-wave used for the imaging scan.
For example, when the third R-wave measured with including that used to start the scanning appears after the waiting for an interval of, for instance, “2 R-R” (YES at step S124-9), the sequencer 5 proceeds to the processing at step S124-4.
Therefore, at a time when the delay time TDL1 passes after the appearance of peak value of the third R-wave, the second-time first scan SNsys2 for the systole is performed again at the next slice-encode amount SE2 in the similar manner to the last one. As a result, echo signals are acquired from the three-dimensional imaging region Rima (steps S124-4 and -5). Further, at another time when the delay time TDL2 passes after the appearance of peak value of the third R-wave, the second-time second scan SNdia2 for the diastole is performed again at the next slice-encode amount SE2 in the similar manner to the last one. As a result, echo signals are acquired from the three-dimensional imaging region Rima (steps S124-6 and S124-7).
Hereafter, echo signals are repeatedly acquired for each of the systole and diastole until the last slice-encode amount SEn (for example, n=8).
When the last scans SNsysn and SNdian at the slice-encode amount SEn have been completed, the determination at step S124-8 becomes YES, and the notification of completion of the imaging scans is issued from the sequencer 5 to the host computer 6 (step S124-10). Thus the processing returns to the host computer 6.
On receiving such notification from the sequencer 5 (step S125 in
Accordingly, as pictorially shown in
The echo data acquired from the patient P are converted into digital echo data in a similar manner to the second embodiment. The echo data are sent to the calculation unit 10 via the sequencer 5, in which they are selectively mapped in three-dimensional systole-use and diastole-use k-spaces Ksys and Kdia both of which are formed by memories, correspondingly to each phase-encode amount and each slice-encode amount.
(3.2) Data Processing and Image Display
After the acquisition of the echo data, the host computer 6 instructs the calculation unit 10 to execute the processing shown in
As shown in
Then, the calculation unit 10 performs a three-dimensional Fourier transform on each of the k-spaces Ksys and Kdia for reconstructing images. As a result, as shown in
Considering those facts, performed in turn by the calculation unit 10 are a subtraction of “IMdia−β·IMsys” for producing an arterial phase image IMAR, a subtraction of “IMdia−·IMAR” for producing a venous phase image IMVE, MIP (maximum intensity projection) processing for each of the arterial phase image IMAR and venous phase image IMVE, two-dimensional display of the arterial phase and venous phase images, and storage of data of those images.
(3.3) Advantages
As described above, for performing the imaging scan with the MRI system of the present embodiment, the read-out gradient pulse GR is applied in a direction substantially in parallel with a flow direction of a blood vessel in the inferior limb. Concurrently, the dephasing pulses Pdephase are added to the read-out gradient pulse GR applied during the systole, whilst the rephasing pulses Prephase are added to the read-out gradient pulse GR applied during the diastole.
Therefore, similarly to control of behaviors of magnetic spins described in the second embodiment, such addition is able to reduce signal intensity by promoting the flow void effect caused in blood flowing in the systole, in particular, in the artery. In contrast, such addition is able to give an effect of flow compensation to the vein and artery flowing in the diastole.
Hence, relative differences of signal intensity between blood that flows at a certain speed and blood that flows at a speed smaller than the certain speed can be made distinctly. So even if the blood vessels in the inferior limb are imaged, which are slower in flow speed than those in the abdomen and thorax, the artery and vein can be visually separated with clearness and displayed with higher depiction capability.
According to the MRI system of the present invention, the optimum scan start timing (delay time from the R-wave) is assigned to each of the systole and diastole in one cardiac cycle. Two shots of scans for the systole and diastole at one slice encode are performed in turn in the one time of imaging scan in an alternating fashion. Additionally, the scan for the systolic, which comes to first in one cardiac cycle, is shortened in time not to overlap with the following scan for the diastole by reducing its data acquisition time (corresponding to the number of echoes). The echo data acquired by such scan are mapped in the lower-frequency region of the k-space for the systole, such region being most significant in terms of improvement in the contrast of images. Short data in the k-space for the systole can be obtained by duplicating data acquired by the following scan for the diastole, which is capable of acquiring echoes over a relative longer period of time. The scans for the systole and diastole use the half Fourier technique to reduce the scan time as short as possible.
Thus, the two shots of scans for the systole and diastole at one slice-encode amount usually remain within an interval of about two heartbeats, Sequentially and alternately repeating such scans makes it possible to acquire echo data of blood flow for the systole and diastole during a breath hold duration of one time. Namely, three-dimensional image data of blood flow for each of the systole and diastole are acquired at its optimum timing through one time of imaging.
There is therefore no need to perform each imaging scan for each of the systole and diastole (that is, in total, two times of imaging), but only one time of imaging is enough. Hence, the imaging time is lessened largely and throughput of patients is improved. Particularly, a reduction in the imaging time becomes noticeable when three-dimensional imaging is carried out. Further, misregistration due to patient's motions can also be reduced greatly, thereby improving image quality. Artery-phase/venous-phase visually separated blood flow images (MRA images) can be obtained from echo data at acquired by one time of imaging, so the imaging is excellent in efficiency. The other operations and advantages obtained in the second embodiment are also available to this third embodiment.
(Modifications of Third Embodiment)
In the foregoing third embodiment, the first-time and second-time imaging scans use, as shown in
(Modifications Common to First to Third Embodiments)
The first to third embodiments can be practiced in other various modified forms
For instance, the foregoing embodiments are configured to present both of the arterial phase and venous phase images, but only the arterial phase image may be produced by a subtraction and displayed. In this case, the step S36 in the processing of
In the foregoing embodiments, the half Fourier technique is used for each scan for each of the systole and diastole, the half Fourier technique may be replaced by other techniques. In that case, one preferred example is that the scan for the diastole acquires echo data that can be mapped in the entire k-space and echo data present in both ended regions (high-frequency regions) are individually duplicated into corresponding regions of the k-space for the systole.
Moreover, the foregoing embodiments have been described using the three-dimensional scan. Instead, the two-dimensional scan can be applied to the embodiments similarly. The pulse sequence to be used is not limited to the FASE technique itself, but pulse sequences based on an FSE technique using an inversion recovery (IR) pulse or an FASE technique modified to use the inversion recovery pulse can be available.
The post-processing of echo data in the foregoing embodiments is configured such that echo data are once converted into image data in the actual space, and then the image data undergo the subtractions to obtain the arterial phase image and venous phase image. Instead of this, the subtractions may be conducted with echo data mapped in the k-spaces Ksys and Kdia, as long as their matrix sizes are equal to each other. The subtracted echo data are then reconstructed into a blood flow image.
Further, as to the configuration to detect a signal indicative of heartbeats of an object to be imaged, the foregoing configuration detecting the ECG signal may be replaced by a PPG (peripheral gating) detector to detect a pulse wave on a finger using an optical signal.
Furthermore, in the MRI system according to each of the foregoing embodiments and modifications, image data at the two time phases are formed into one set of image data, but the present invention is not limited to this mode. For instance, the read-out gradient pulse to which the dephasing pulses or rephasing pulses are added is set so that it is applied almost in parallel with a flow direction of fluid (blood, lymph, or others). Then, an imaging scan using the read-out gradient pulse is performed one time to obtain a single image, with no relation to the cardiac time phases. The fluid is imaged in bright or black into this image, with degrees of promotion of the flow void effect in the fluid reflected. Therefore, this image also provides flow information about the fluid.
As another modification, means for controlling the intensities of the foregoing dephasing pulses and rephasing pulses according to flow speeds of fluid to be imaged can be provided. This means is composed of, for example, the input device 13, host computer 6, and/or storage unit 11. In response to information indicative of both of a region to be imaged and fluid to be imaged, which is provided by an operator via the input device 13, the host computer 6 refers to a memory table previously stored in the storage unit 11. The table memorizes pulse intensities fluid by fluid. The host computer 6 provides the sequencer 5 with the intensity of a dephasing pulse or rephasing pulse according to the reference result. Alternatively, an operator is also able to use input device 13 for directly giving the system desired pulse intensity.
Though the embodiments have been described above, the present invention is not restricted to the configurations described in the embodiments, and various modifications and adaptations will be readily apparent to those skilled in the art without departing from the substance of scope of the present invention. Those modifications and adaptations should be construed as being included into the present invention.
Miyazaki, Mitsue, Sugiura, Satoshi
Patent | Priority | Assignee | Title |
10262385, | Dec 12 2014 | SIEMENS HEALTHINEERS AG | Method and apparatus for acquiring magnetic resonance data and entering the data into k-space |
10451696, | Sep 26 2016 | Samsung Electronics Co., Ltd. | Magnetic resonance imaging apparatus and method of obtaining magnetic resonance image |
10788554, | Feb 22 2017 | SIEMENS HEALTHINEERS AG | Method for setting an operating parameter of a medical device and medical device |
9274196, | Nov 15 2011 | SIEMENS HEALTHINEERS AG | Method and magnetic resonance system to generate angiography images |
Patent | Priority | Assignee | Title |
4609872, | Aug 10 1984 | General Electric Company | NMR multiple-echo phase-contrast blood flow imaging |
4752734, | Aug 10 1984 | The General Hospital Corporation | Flow imaging by means of nuclear magnetic resonance |
4995394, | Aug 11 1989 | General Electric Company | Fast NMR cardiac profile imaging |
5101156, | Aug 09 1990 | General Electric | Rapid flow measurement using an NMR imaging system |
5122747, | Nov 21 1990 | Mayo Foundation for Medical Education and Research | Spiral three-dimensional fourier transform NMR scan |
5189369, | Mar 19 1990 | Hitachi Medical Corporation | NMR imaging method of low flow rate fluid |
5221898, | Nov 30 1990 | Hitachi Medical Corporation | Flow imaging method using an MRI apparatus |
5320099, | Aug 07 1992 | Trustees of the University of Pennsylvania | MR angiography using steady-state transport-induced adiabatic fast passage |
5417214, | Aug 07 1992 | Trustees of the University of Pennsylvania | Quantitative blood flow measurement using steady-state transport-induced adiabatic fast passage |
5519320, | Mar 31 1992 | Kabushiki Kaisha Toshiba | Nuclear magnetic resonance imaging with high speed and interactive pulse sequence control |
5565776, | Aug 26 1994 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging apparatus |
5565777, | Sep 13 1993 | Kabushiki Kaisha Toshiba | Method/apparatus for NMR imaging using an imaging scheme sensitive to inhomogeneity and a scheme insensitive to inhomogeneity in a single imaging step |
5830143, | Jan 21 1997 | Wisconsin Alumnin Research Foundation | Gated time-resolved contrast-enhanced 3D MR angiography |
5857970, | Jun 20 1997 | Siemens Medical Solutions USA, Inc | Method and apparatus for cardiac-synchronized peripheral magnetic resonance angiography |
5897496, | Sep 03 1996 | Hitachi Medical Corporation | Method and apparatus for producing magnetic resonance angiogram |
5910728, | Nov 12 1996 | BETH ISRAEL DEACONESS MEDICAL CENTER, INC | Simultaneous acquisition of spatial harmonics (SMASH): ultra-fast imaging with radiofrequency coil arrays |
6043655, | Jan 09 1997 | Toshiba Medical Systems Corporation | MR imaging utilizing the time of flight effect |
6144201, | Dec 26 1997 | Toshiba Medical Systems Corporation | MR imaging utilizing ECG gating technique |
6320377, | Mar 28 1997 | Toshiba Medical Systems Corporation | MR imaging with pre-sequence including plural saturation pulses |
6353752, | May 14 1999 | Board of Trustees of the Leland Stanford Junior University | Reduced field-of-view method for cine magnetic resonance imaging |
6493569, | Dec 30 2000 | GE Medical Systems Global Technology Company, LLC; Uniformed Services University of Health Sciences | Method and apparatus using post contrast-enhanced steady-state free precession in MR imaging |
6782286, | Apr 20 1998 | Toshiba Medical Systems Corporation | MRI system and MR imaging method |
JP10234695, | |||
JP11225992, | |||
JP11239571, | |||
JP11276453, | |||
JP2000023937, | |||
JP2000300538, | |||
JP2000342553, | |||
JP2000342555, | |||
JP243494, | |||
JP353936, | |||
JP4035647, | |||
JP4309331, | |||
JP5277085, | |||
JP8131419, |
Executed on | Assignor | Assignee | Conveyance | Frame | Reel | Doc |
Oct 04 2006 | Kabushiki Kaisha Toshiba | (assignment on the face of the patent) | / | |||
Mar 16 2016 | Kabushiki Kaisha Toshiba | Toshiba Medical Systems Corporation | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 038847 | /0108 |
Date | Maintenance Fee Events |
Date | Maintenance Schedule |
Oct 16 2015 | 4 years fee payment window open |
Apr 16 2016 | 6 months grace period start (w surcharge) |
Oct 16 2016 | patent expiry (for year 4) |
Oct 16 2018 | 2 years to revive unintentionally abandoned end. (for year 4) |
Oct 16 2019 | 8 years fee payment window open |
Apr 16 2020 | 6 months grace period start (w surcharge) |
Oct 16 2020 | patent expiry (for year 8) |
Oct 16 2022 | 2 years to revive unintentionally abandoned end. (for year 8) |
Oct 16 2023 | 12 years fee payment window open |
Apr 16 2024 | 6 months grace period start (w surcharge) |
Oct 16 2024 | patent expiry (for year 12) |
Oct 16 2026 | 2 years to revive unintentionally abandoned end. (for year 12) |