MR imaging system with interactive MR geometry prescription control

Information

  • Patent Grant
  • 6522141
  • Patent Number
    6,522,141
  • Date Filed
    Tuesday, October 16, 2001
    22 years ago
  • Date Issued
    Tuesday, February 18, 2003
    21 years ago
Abstract
A magnetic resonance (MR) imaging system equipped with real-time imaging capability and methods of interactively prescribing geometry to excitation profiles of structure of interest, are disclosed herein. The MR imaging system includes a graphical user interface for displaying and receiving prescription commands, a display screen for displaying MR images and the graphical user interface, and an input device for inputting prescription commands. The MR imaging system allows an operator to prescribe the boundary geometry of a subsequent imaging volume and to rapidly view the prescribed boundary imaging sections prior to committing to the subsequent imaging volume acquisition. The MR imaging system also allows the operator to retrieve boundary geometry of a previously prescribed imaging volume and to rapidly view the imaging sections corresponding to the retrieved boundary geometry prior to initiating the image volume acquisition.
Description




FIELD OF THE INVENTION




The present invention relates generally to magnetic resonance (MR) imaging systems and methods. More particularly, the present invention relates to a MR imaging system equipped for real-time imaging and methods for assisting the operator to interactively prescribe the geometry of the excitation profile of a structure of interest for subsequent acquisition of a MR image of the structure of interest.




BACKGROUND OF THE INVENTION




When a substance such as human tissue is subjected to a uniform magnetic field (polarizing field Bo), the individual magnetic moments of the spins in the tissue attempt to align with this polarizing field, but process about it in random order at their characteristic Larmor frequency. If the substance, or tissue, is subjected to a magnetic field (excitation field B


1


) which is the x-y plane and which is near the Larmor frequency, the net aligned moment, M


z


, may be rotated, or “tipped”, into the x-y plane to produce a net transverse magnetic moment M. A signal is emitted by the excited spins after the excitation signal B


1


is terminated and this signal may be received and processed to form an image.




When utilizing these signals to produce images, magnetic field gradients (G


x


, G


y


and G


z


) are employed. Typically, the region to be imaged is scanned by a sequence of measurement cycles in which these gradients vary according to the particular localization method being used. The resulting set of received NMR signals are digitized and processed to reconstruct the image using one of many well known reconstruction techniques.




When attempting to define the volume of coverage of an MR imaging scan, the NMR system operator may desire to quickly view a preview MR image (such as a real-time MR image) of the anatomical section within this volume of coverage. This process can be particularly useful when prescribing a three dimensional imaging volume, in which the desired high spatial resolution requires the thinnest slab possible. It is desirable to position this thin slab such that the anatomical section within the volume of coverage is complete, i.e. for example, covers the entire desired vascular network. Thus, a quick view of each side of the slab prior to initiating the three dimensional acquisition is useful for insuring that the entire anatomical section desired is within the defined volume of coverage.




Typically, two dimensional axial, sagittal and coronal “scout” images are first acquired. Such scout images are stored for later use. To use, the operator calls up the scout image and either graphically or explicitly (using geometry coordinates) prescribes the imaging volume directly on the scout images. The imaging volume may be either a two dimensional stack of slices or a three dimensional slab of the structure of interest. The drawback of this technique is that the operator does not actually see the results of the prescribed geometry until the subsequent imaging volume is acquired. Prescription errors cannot be detected nor corrected until the imaging volume acquisition is complete. Thus, when prescription errors exist, the operator is required to re-prescribe and re-acquire the imaging volume of the desired anatomical section.




SUMMARY OF THE INVENTION




One embodiment of the invention relates to a method for defining an imaging plane of an image to be acquired of a structure of interest. The method includes selecting a previously acquired three-dimensional image of the structure of interest, and determining a first geometry information associated with a first boundary plane of the previously acquired three-dimensional image. The method further includes determining a second geometry information associated with a second boundary plane of the previously acquired three-dimensional image, and selecting at least one of the first geometry information and the second geometry information to define the imaging plane. The image to be acquired is a planar image.




Another embodiment of the invention relates to a system for defining an imaging plane of an image to be acquired of a structure of interest. The system includes means for selecting a previously acquired three-dimensional image of the structure of interest, and means for determining a first geometry information associated with a first boundary plane of the previously acquired three-dimensional image. The system further includes means for determining a second geometry information associated with a second boundary plane of the previously acquired three-dimensional image, and means for selecting at least one of the first geometry information and the second geometry information to define the imaging plane. The image to be acquired is a planar image.




Still another embodiment of the invention relates to an imaging system retrieving geometry prescription information from an image volume of a structure of interest to define an imaging plane of an image to be acquired. The system includes a storage device configured to store the image volume, and an interface in communication with the storage device and configured to transmit at least one selection signal in response to an operator selecting the image volume on the interface. The system further includes a system control in communication with the storage device and the interface. The system control is configured to determine at least one of a first geometry prescription information and a second geometry prescription information associated with a first boundary plane and a second boundary plane, respectively, of the image volume. At least one of the first and second geometry prescription information defines the imaging plane of the image to be acquired.




Yet still another embodiment of the invention relates to a planar magnetic resonance (MR) image of a structure of interest provided using an MR imaging system. The image comprises an imaging plane defined by a geometry parameter associated with a boundary plane of at least one of a first previously acquired volume image of the structure of interest and a second previously acquired volume image of the structure of interest.











BRIEF DESCRIPTION OF THE DRAWINGS




The invention will become more fully understood from the following detailed description, taken in conjunction with the accompanying drawings, wherein like reference numerals refer to like parts, in which:





FIG. 1

is a block diagram of a MR imaging system which employs the present invention;





FIG. 2

is an electrical block diagram of the transceiver which forms part of the MR imaging system of

FIG. 1

; and





FIG. 3

is an illustration of the graphical user interface on the display screen of the operator console of the MR imaging system of FIG.


1


.











DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS




Referring first to

FIG. 1

, there is shown the major components of a preferred MR imaging system which incorporates the present invention. The operation of the system is controlled from an operator console


100


which includes an input device


101


, a control panel


102


and a display


104


. The console


100


communicates through a link


116


with a separate computer system


107


that enables an operator to control the production and display of images on the display


104


. The computer system


107


includes a number of modules which communicate with each other through a backplane. These include an image processor module


106


, a CPU module


108


and a memory module


113


, known in the art as a frame buffer for storing image data arrays. The computer system


107


is linked to a disk storage


111


and a tape drive


112


for storage of image data and programs, and it communicates with a separate system control


122


through a high speed serial link


115


.




The system control


122


includes a set of modules connected together by a backplane. These include a CPU module


119


and a pulse generator module


121


which connects to the operator console


100


through a serial link


125


. It is through this link


125


that the system control


122


receives commands from the operator which indicate the scan sequence that is to be performed. The pulse generator module


121


operates the system components to carry out the desired scan sequence. It produces data which indicates the timing, strength and shape of the RF pulses which are to be produced, and the timing of and length of the data acquisition window. The pulse generator module


121


connects to a set of gradient amplifiers


127


, to indicate the timing and shape of the gradient pulses to be produced during the scan. The pulse generator module


121


also receives patient data from a physiological acquisition controller


129


that receives signals from a number of different sensors connected to the patient, such as ECG signals from electrodes or respiratory signals from a bellows. And finally, the pulse generator module


121


connects to a scan room interface circuit


133


which receives signals from various sensors associated with the condition of the patient and the magnet system. It is also through the scan room interface circuit


133


that a patient positioning system


134


receives commands to move the patient to the desired position for the scan.




The gradient waveforms produced by the pulse generator module


121


are applied to a gradient amplifier system


127


comprised of G


x


, G


y


and G


z


amplifiers. Each gradient amplifier excites a corresponding gradient coil in an assembly generally designated


139


to produce the magnetic field gradients used for position encoding acquired signals. The gradient coil assembly


139


forms part of a magnet assembly


141


which includes a polarizing magnet


140


and a whole-body RF coil


152


.




A transceiver module


150


in the system control


122


produces pulses which are amplified by an RF amplifier


151


and coupled to the RF coil


152


by a transmit/receiver switch


154


. The resulting signals radiated by the excited nuclei in the patient may be sensed by the same RF coil


152


and coupled through the transmit/receive switch


154


to a preamplifier


153


. The amplified NMR signals are demodulated, filtered, and digitized in the receiver section of the transceiver


150


. The transmit/receive switch


154


is controlled by a signal from the pulse generator module


121


to electrically connect the RF amplifier


151


to the coil


152


during the transmit mode and to connect the preamplifier


153


during the receive mode. The transmit/receive switch


154


also enables a separate RF coil (for example, a head coil or surface coil) to be used in either the transmit or receive mode.




The NMR signals picked up by the RF coil


152


are digitized by the transceiver module


150


and transferred to a memory module


160


in the system control


122


. When the scan is completed and an entire array of data has been acquired in the memory module


160


, an array processor


161


operates to Fourier transform the data into an array of image data. This image data is conveyed through the serial link


115


to the computer system


107


where it is stored in the disk memory


111


. In response to commands received from the operator console


100


, this image data may be archived on the tape drive


112


, or it may be further processed by the image processor


106


and conveyed to the operator console


100


and presented on the display


104


.




Referring particularly to

FIGS. 1 and 2

, the transceiver


150


produces the RF excitation field B


1


through power amplifier


151


at a coil


152


A and receives the resulting signal induced in a coil


152


B. As indicated above, the coils


152


A and B may be separate as shown in

FIG. 2

, or they may be a single wholebody coil as shown in FIG.


1


. The base, or carrier, frequency of the RF excitation field is produced under control of a frequency synthesizer


200


which receives a set of digital signals (CF) from the CPU module


119


and pulse generator module


121


. These digital signals indicate the frequency and phase of the RF carrier signal produced at an output


201


. The commanded RF carrier is applied to a modulator and up converter


202


where its amplitude is modulated in response to a signal R(t) also received from the pulse generator module


121


. The signal R(t) defines the envelope of the RF excitation pulse to be produced and is produced in the module


121


by sequentially reading out a series of stored digital values. These stored digital values may, in turn, be changed from the operator console


100


to enable any desired RF pulse envelope to be produced.




The magnitude of the RF excitation pulse produced at output


205


is attenuated by an exciter attenuator circuit


206


which receives a digital command, TA, from the backplane


118


. The attenuated RF excitation pulses are applied to the power amplifier


151


that drives the RF coil


152


A. For a more detailed description of this portion of the transceiver


122


, reference is made to U.S. Pat. No. 4,952,877 which is incorporated herein by reference.




Referring still to

FIGS. 1 and 2

the NMR signal produced by the subject is picked up by the receiver coil


152


B and applied through the preamplifier


153


to the input of a receiver attenuator


207


. The receiver attenuator


207


further amplifies the signal by an amount determined by a digital attenuation signal (RA) received from the backplane


118


.




The received signal is at or around the Larmor frequency, and this high frequency signal is down converted in a two step process by a down converter


208


which first mixes the NMR signal with the carrier signal on line


201


and then mixes the resulting difference signal with the 2.5 MHz reference signal on line


204


. The down converted NMR signal is applied to the input of an analog-to-digital (A/D) converter


209


which samples and digitizes the analog signal and applies it to a digital detector and signal processor


210


which produces 16 bit in-phase (I) values and 16-bit quadrature (Q) values corresponding to the received signal. The resulting stream of digitized I and Q values of the received signal are output through backplane


118


to the memory module


160


where they are normalized in accordance with the present invention and then employed to reconstruct an image.




The 2.5 MHz reference signal as well as the 250 kHz sampling signal and the 5, 10 and 60 MHz reference signals are produced by a reference frequency generator


203


from a common 20 MHz master clock signal. For a more detailed description of the receiver, reference is made to U.S. Pat. No. 4,992,736 which is incorporated herein by reference.




In one embodiment of the present invention, an operator interactively prescribes geometry to define a subsequent MR imaging volume or receives geometry information from a previously defined MR imaging volume of the structure of interest, such as an anatomical structure. Such interactive geometry prescription is accomplished from the operator console


100


(also referred to as an operator interface) using the input device


101


. The input device


101


is selected from a group including, but not limited to, a mouse, a joystick, a keyboard, a track ball, a touch screen, a light wand, and a voice control. The MR imaging system of the present invention is capable of imaging in any desired orientation within the structure of interest and is equipped to perform both real-time acquisitions and non real-time acquisitions. In particular, real-time refers to continuous acquisition and reconstruction of MR image data as rapidly as it is acquired. A real-time MR image can be acquired and displayed in approximately one second or less, as constrained by MR imaging system performance.





FIG. 3

shows a graphical user interface


105


used in an embodiment of the present invention. The graphical user interface


105


and the MR image of the structure of interest is displayed on the display


104


(also referred to as an electronic display) of the MR imaging system. The operator interacts with the graphical user interface


105


using the input device


101


. The graphical user interface


105


includes a set start boundary icon


10


, a three-point start boundary geometry icon


12


, a set end boundary icon


14


, and a three-point end boundary geometry icon


16


. The three-point start and end boundary geometry icons


12


,


16


, respectively, each contain geometry coordinates defining the location of a planar section of the structure of interest in the imaging volume. These coordinates are defined in the patient right-left direction (R/L), patient anterior-posterior direction (A/P), and patient superior-inferior direction (S/I), hereafter referred to as center point RAS coordinates. The graphical user interface


105


also includes an acquire start boundary icon


18


, an acquire end boundary icon


20


, an apply location icon


22


, a retrieve location icon


24


, and a save series icon


26


.




First, to prescribe the boundary geometry of a subsequent or proposed imaging volume, it is desirable for the operator to view real-time imaging sections, preferably two dimensional planar sections, corresponding to the boundaries defining the desired subsequent imaging volume prior to committing to those imaging sections as the boundaries of the subsequent imaging volume. Typically the operator maneuvers the MR imaging system to acquire and display a real-time imaging section on display


104


directed to the structure of interest that defines one boundary of the desired subsequent imaging volume. The operator then registers this real-time imaging section as one boundary plane of the subsequent imaging volume by “clicking” on the set start boundary icon


10


on the graphical user interface


105


. A geometry representation of the scan plane of this imaging section is determined and stored (i.e. in a text buffer) as text in center point RAS coordinates. The geometry representation of the start boundary is also displayed in the three-point start boundary geometry icon


12


of the graphical user interface


105


.




Next, the operator manipulates the MR imaging system to acquire and display another real-time imaging section on display


104


directed to the structure of interest that defines another boundary of the desired subsequent imaging volume. The operator registers this current real-time imaging section as another boundary plane of the subsequent imaging volume by clicking on the set end boundary icon


14


on the graphical user interface


105


. Similar to above, a geometry representation of the scan plane of this current imaging section is determined, stored, and displayed in center point RAS coordinates in the three-point end boundary geometry icon


16


of the graphical user interface


105


.




It should be understood that non real-time imaging sections can also be utilized to set the start and end boundaries. The advantage of the real-time imaging sections is that the operator can very rapidly view multiple imaging sections of interest for the purposes of prescribing the subsequent imaging volume. Additionally, the operator can repeatedly set the start and/or end boundary planes by acquiring and displaying a new imaging section and then clicking on the set start boundary icon


10


or the set end boundary icon


14


, as desired. In this way, the present embodiment provides the operator with a finer degree of geometry prescription control.




The remaining boundary geometry defining the subsequent imaging volume can be identical to the corresponding boundaries of the current real-time imaging section, i.e., the in-plane field of view. Alternatively, the remaining boundary geometry can be defined independently with additional icons on the graphical user interface


105


using the input device


101


(not shown in FIG.


3


). Still further, in the case where the two boundary planes are not parallel to each other, the MR imaging system can apply a best fit algorithm, or other suitable algorithms, to the start and end boundaries to calculate the remaining boundary geometry.




The operator can now click on the apply location icon


22


, which transfers the start and end boundary geometry information contained in icons


12


,


16


to the subsequent imaging volume. Once the start and end boundary geometry information has been applied, the operator can click on the save series icon


26


. This signals the MR imaging system to check for a complete boundary geometry prescription and prepares the system for acquisition of the prescribed imaging volume.




Second, to retrieve the boundary geometry of a previously prescribed or defined imaging volume and to utilize the retrieved geometry information to check the prescribed boundary geometry or to use it as a starting point from which to prescribe a subsequent imaging volume, the operator starts by selecting a previously prescribed imaging volume from a list or display of one or more previously prescribed imaging volumes on display


104


(not shown in FIG.


3


). The previously prescribed imaging volumes can be, but is not limited to, previously stored real-time acquisitions, previously stored non real-time acquisitions, or previously stored graphically or explicitly (using geometry coordinates) prescribed imaging volumes from scout images. Then the operator clicks on the retrieve location icon


24


to load boundary geometry information, in center point RAS coordinates, into the buffers corresponding to icons


12


,


16


. Icons


12


,


16


displays the two boundary plane geometry information.




Using the acquire start boundary icon


18


or the acquire end boundary icon


20


, the operator commands the MR imaging system to acquire and display a real-time imaging section, typically a two-dimensional planer section, defined by the retrieved geometry information in the three-point start boundary geometry icon


12


or the three-point end boundary geometry icon


16


, respectively. Alternatively, the retrieved geometry information can be used to acquire and display a non real-time imaging section. The feature embodied in the acquire start and end boundary icons


18


,


20


are particularly useful for checking or previewing the boundaries of a previously prescribed imaging volume that has not been acquired, such as an imaging volume prescribed using scout images.




In another embodiment of the present invention, the imaging section acquired and displayed as a result of clicking the acquire start or end boundary icon


18


,


20


can be modified such that an acquisition of a new imaging section occurs and the said section is displayed (replacing the current imaging section displayed). The modification, for example, can be accomplished by graphically or explicitly (using geometry coordinates) changing the scan plane of the currently imaging section. This new imaging section, in turn, can be utilized to replace the retrieved geometry information stored in icon


12


or


16


by clicking on the set start or end boundary icon


10


or


14


, respectively. Thus, in this manner, the geometry information of a previously prescribed imaging volume can be used as a starting point from which to prescribe a subsequent imaging volume or to refine the prescription of a previously prescribed imaging volume.




It should be apparent that there has been provided in accordance with one embodiment of the present invention a method for accurately and efficiently prescribing the geometry of a subsequent imaging volume of a structure of interest using at least two two-dimensional MR imaging sections. Moreover, an embodiment of the present invention also provides a method for retrieving geometry information from a previously prescribed imaging volume and manipulating this geometry information. While the embodiments illustrated in the FIGs. and described above are presently preferred, it should be understood that these embodiments are offered by way of example only. For example, setting the start or end boundary described herein may be accomplished directly by inputting geometry coordinates rather than by displaying an imaging section and extracting or determining geometry coordinates therefrom. Accordingly, the invention is not limited to a particular embodiment, but extends to alternatives, modifications, and variations that nevertheless fall within the spirit and scope of the appended claims.



Claims
  • 1. A method for defining an imaging plane of an image to be acquired of a structure of interest, the method comprising:selecting a previously acquired three-dimensional image of the structure of interest; determining a first geometry information associated with a first boundary plane of the previously acquired three-dimensional image; determining a second geometry information associated with a second boundary plane of the previously acquired three-dimensional image; and selecting at least one of the first geometry information and the second geometry information to define the imaging plane, wherein the image to be acquired is a planar image.
  • 2. The method of claim 1, wherein at least one of the selecting steps is performed by a user using an input device selected from a group including a mouse, a joystick, a keyboard, a trackball, a touch screen, a light wand, and a voice control.
  • 3. The method of claim 1, further comprising storing the first and second geometry information.
  • 4. The method of claim 3, further comprising:selecting an another previously acquired three-dimensional image of the structure of interest; and determining an another geometry information associated with a boundary plane of the another previously acquired three-dimensional image.
  • 5. The method of claim 4, further comprising replacing at least one of the stored first and second geometry information with the another geometry information.
  • 6. The method of claim 1, further comprising acquiring the image to be acquired by at least one of real-time acquisition and non real-time acquisition.
  • 7. A system for defining an imaging plane of an image to be acquired of a structure of interest, the system comprising:means for selecting a previously acquired three-dimensional image of the structure of interest; means for determining a first geometry information associated with a first boundary plane of the previously acquired three-dimensional image; means for determining a second geometry information associated with a second boundary plane of the previously acquired three-dimensional image; and means for selecting at least one of the first geometry information and the second geometry information to define the imaging plane, wherein the image to be acquired is a planar image.
  • 8. The system of claim 7, wherein at least one of the means for selecting is performed by a user using an input device selected from a group including a mouse, a joystick, a keyboard, a trackball, a touch screen, a light wand, and a voice control.
  • 9. The system of claim 7, further comprising means for storing the first and second geometry information.
  • 10. The system of claim 9, further comprising:means for selecting an another previously acquired three-dimensional image of the structure of interest; and means for determining an another geometry information associated with a boundary plane of the another previously acquired three-dimensional image.
  • 11. The system of claim 10, further comprising means for replacing at least one of the stored first and second geometry information with the another geometry information.
  • 12. The system of claim 7, further comprising means for acquiring the image to be acquired by at least one of real-time acquisition and non real-time acquisition.
  • 13. An imaging system retrieving geometry prescription information from an image volume of a structure of interest to define an imaging plane of an image to be acquired, the system comprising:a storage device configured to store the image volume; an interface in communication with the storage device and configured to transmit at least one selection signal in response to an operator selecting the image volume on the interface; and a system control in communication with the storage device and the interface, the system control configured to determine at least one of a first geometry prescription information and a second geometry prescription information associated with a first boundary plane and a second boundary plane, respectively, of the image volume, wherein the at least one of the first and second geometry prescription information defines the imaging plane of the image to be acquired.
  • 14. The system of claim 13, wherein the image to be acquired is a planar image of the structure of interest.
  • 15. The system of claim 13, wherein the image volume is at least one of a three-dimensional image and a stack of planar images.
  • 16. The system of claim 13, further comprising an another storage device in communication with the system control, the another storage device configured to store at least one of the first and second geometry information determined by the system control.
  • 17. The system of claim 16, wherein the storage device is configured to store an another image volume.
  • 18. The system of claim 17, wherein the operator selects the another image volume, the system control determines an another geometry prescription information associated with the another image volume, and the another storage device stores the another geometry prescription information.
  • 19. The system of claim 18, further comprising an image acquisition unit coupled to the system control, the image acquisition unit configured to acquire the image to be acquired using at least one of the first geometry prescription, the second geometry prescription, and the another geometry prescription information as the imaging plane, as selected by the operator.
  • 20. The system of claim 13, wherein the image to be acquired is a planar image.
  • 21. A planar magnetic resonance (MR) image of a structure of interest provided using an MR imaging system, the image comprising an imaging plane defined by a geometry parameter associated with a boundary plane of at least one of a first previously acquired volume image of the structure of interest and a second previously acquired volume image of the structure of interest.
  • 22. The image of claim 21, wherein an operator selects at least one of the first and second previously acquired volume image and the MR imaging system determines the geometry parameter of the selected previously acquired volume image, the geometry parameter being at least one of a first boundary plane and a second boundary plane of the selected previously acquired volume image, and the first and second boundary planes being displaced from each other.
CROSS REFERENCE TO RELATED APPLICATIONS

The present application is a divisional of U.S. application Ser. No. 09/200,144 filed on Nov. 25, 1998, now U.S. Pat. No. 6,396,266 and which is incorporated in its entirety herein by reference.

US Referenced Citations (19)
Number Name Date Kind
4830012 Riederer May 1989 A
5270651 Wehrli Dec 1993 A
5280428 Wu et al. Jan 1994 A
5345176 LeRoux et al. Sep 1994 A
5451876 Sandford et al. Sep 1995 A
5498963 Schneider et al. Mar 1996 A
5512826 Hardy et al. Apr 1996 A
5512827 Hardy et al. Apr 1996 A
5514962 Cline et al. May 1996 A
5541513 Maier Jul 1996 A
5560361 Glusick Oct 1996 A
5657757 Hurd et al. Aug 1997 A
5711300 Schneider et al. Jan 1998 A
5749834 Hushek May 1998 A
5810729 Hushek et al. Sep 1998 A
6023653 Liu et al. Feb 2000 A
6108573 Debbins et al. Aug 2000 A
6166544 Debbins et al. Dec 2000 A
6396266 Debbins et al. May 2002 B1
Foreign Referenced Citations (1)
Number Date Country
WO 9824058 Jun 1998 WO
Non-Patent Literature Citations (4)
Entry
Debbins, et al. Cardiac Magnetic Resonance Fluoroscopy. Magnetic Resonance In Medicine, vol. 38, (1996), pp. 588-595.
Hangiandreou, et al. Interactive Selection of Optimal Section Orientations Using Real-Time MRI. Magnetic Resonance In Medicine, vol. 34 (1995), pp. 114-119.
Hardy, C. J. et al. Interactive Coronary MRI. Magnetic Resonance in Medicine, Jul. 1998, vol. 40, No. 1, pp. 105-111.
Wilman, A. H. et al. Aterial Phase Carotid And Vertebral Artery Imaging In 3D Contrast-Enhanced MR Angiography By Combining Fluoroscopic Triggering With An Elliptical Centric Acquisition Order, Magnetic Resonance in Medicine, Jul. 1998, vol. 40, No. 1, pp. 24-35.