Patent Application
20140070804
The present application relates to magnetic resonance arts. It finds particular application in channel-by-channel reconstruction algorithms employed in parallel imaging methods.
The introduction of multi-coil arrays in magnetic resonance (MR) systems to improve the signal-to-noise ratio (SNR) over volume coils and large surface coils has led to the introduction of multiple channel receivers included in the systems. More recently, the success of partially parallel imaging (PPI) techniques is driving the industry to develop MRI systems with a greater number of receiver channels to enable parallel imaging with higher acceleration factors and broader coverage. These techniques use spatial information contained in the component coils of an array to partially replace spatial encoding which would normally be performed using gradients, thereby greatly reducing imaging time.
Commercial MRI systems with 32 or more receiver channels are becoming popular for regular clinical use. With the increase of the number of channels, the computation time required for reconstruction has greatly increased. For reconstruction algorithms using single channel information, such as partial Fourier and k-t FOCUSS, the reconstruction times increase linearly with channel count. For channel-by-channel PPI reconstruction techniques, such as Generalized Auto-calibrating Partially Parallel Acquisitions (GRAPPA) and Data Convolution and Combination Operation (COCOA), the reconstruction time increases quadratically. Thus reconstruction with a 32-channel coil may take 16 times longer than an 8-channel coil which renders real-time and high-throughput imaging difficult.
To reduce reconstruction time without significantly degrading the final image quality, two different sets of channel reduction techniques have been proposed. One set of techniques, channel compression, reduces reconstruction time by first linearly combining the original channels using principle component analysis (PCA) into a lesser number of channels. Of those combined channels, the ones with the largest eigenvalues are used for reconstruction. The degree of channel compression is limited by the need to avoid significant signal loss. For example, to if the number of compressed channels is fewer than 12 for a 32 channel cardiac coil, image quality is sacrificed.
The other set of techniques reconstruct only one virtual composite channel using all original channels. This set of techniques includes direct virtual coil (DVC) reconstruction and synthetic target (ST) reconstruction. Because only one virtual channel is reconstructed, the reconstruction time can be greatly reduced. However, these techniques, DVC and ST, suffer from lack of a technique to efficiently optimize the phase definition to produce an optimal calibration signal. Without an accurate calibration signal, the final reconstruction can be significantly compromised or even damaged.
The present application provides a new and improved method and system which overcomes the above-referenced problems and others.
In accordance with one aspect, a method of magnetic resonance imaging is provided. A first image representation of an examination region is generated based on a plurality of partial k-space data sets. Each data is associated with at least one channel. A relative sensitivity map is generated for each of a plurality of coil elements, each coil element being associated with at least one channel. The first image representation and each relative sensitivity map are projected to generate a plurality of recreated k-space data sets. Each recreated k-space data set corresponds to one of the partial k-space data sets. Each of the partial k-space data sets is combined with the corresponding recreated k-space data set to generate substituted k-space data sets.
In accordance with another aspect, a computer-readable medium carries software for controlling one or more processors to perform the method of the preceding paragraph.
In accordance with another aspect, a magnetic resonance imaging system includes a plurality of coil elements and corresponding receivers which define a plurality of channels. Data from the channels generated during a magnetic resonance imaging sequence constitute a plurality of partial k-space data sets. One or more processors are programmed to generate a first image representation of an examination region based on the plurality of partial k-space data sets. A relative sensitivity map is generated for each of the channels. The first image representation with each of the relative sensitivity maps is projected to generate a corresponding plurality of recreated k-space data sets. Each partial k-space data set and the corresponding recreated k-space data set are combined to generate substituted k-space data sets.
In accordance with another aspect, an image reconstruction system which generates a plurality of original channels of partial k-space data includes one or more processors. The one or more processors are programmed to compress the plurality of original channels of k-space data into a fewer number of compressed channels. A channel-by-channel reconstruction algorithm is applied to the partial k-space data of one or more of the compressed channels to generate a virtual coil image. A relative sensitivity map is calculated for each of the original or compressed channels. The virtual composite image coil is projected using the relative sensitivity maps to generate a recreated sensitivity map corresponding to each original or compressed channel. Acquired data from the partial k-space data of each of the original or compressed channels is inserted into each of the corresponding recreated k-space data sets to generate substituted k-space data sets. The recreated k-space data sets are reconstructed and combined to produce a final image.
One advantage resides in that reconstruction time is reduced.
Another advantage resides in that image degradation in partially parallel imaging techniques is reduced.
Still further advantages of the present invention will be appreciated to those of ordinary skill in the art upon reading and understand the following detailed description.
The invention may take form in various components and arrangements of components, and in various steps and arrangements of steps. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention.
With reference to
A radio-frequency (RF) coil array, such as a whole-body radio frequency coil, is disposed adjacent the examination region. The RF coil array includes a plurality of individual RF coil elements 18, or may be a birdcage-type coil with the multiple elements 18 interconnected by end ring RF coil structures. The RF coil array generates radio frequency pulses for exciting magnetic resonance in aligned dipoles of the subject. In some embodiments, the radio frequency coil assembly 18 also serves to detect magnetic resonance signals emanating from the imaging region. In other embodiments, local or surface RF coils 18′ are provided in addition to or instead of the whole-body RF coil for more sensitive, localized spatial encoding, excitation, and reception of magnetic resonance signals. The individual RF coils 18 together can act a single coil, as a plurality of independent coil elements, as an array such as in a parallel transmit system, or a combination.
A scan controller 20 controls a gradient controller 22 which causes the gradient coils to apply selected phase encode gradients across the imaging region, as may be appropriate to a selected magnetic resonance imaging or spectroscopy sequence. The scan controller 20 also controls an RF transmitter 24 which causes the whole-body or local
RF coils to generate magnetic resonance excitation and manipulation B1 pulses. The scan controller also controls an RF receiver 26 which is connected to the RF coils 18, and/or a dedicated receive coil 18′ placed inside the examination region 14, to receive magnetic resonance signals therefrom. In a parallel system, the RF receiver 26 includes a plurality of receivers or a single receiver with a plurality of receive channels, each receive channel is operatively connected to at least one corresponding coil element 18 of the assembly. For example, 32 coil elements with 35 transmitters can provide 32 transmit channels and the 32 coil elements with 32 corresponding receivers can define 32 receive channels.
The received data from the receiver channels are temporarily stored in a data buffer 28 and processed by a magnetic resonance data processor 30. The magnetic resonance data processor can perform various functions as are known in the art, including image reconstruction, magnetic resonance spectroscopy processing, catheter or interventional instrument localization, and the like. Reconstructed magnetic resonance images, spectroscopy readouts, interventional instrument location information, and other processed MR data are displayed on a graphical user interface 32. The graphic user interface 32 also includes a user input device which a clinician can use for controlling the scan controller 20 to select scanning sequences and protocols, and the like.
Single channel, such as k-t FOCal Underdetermined System Solver (k-t FOCUSS), and channel-by-channel reconstruction techniques, such as Generalized Auto-calibrating Partially Parallel Acquisitions (GRAPPA), Data Convolution and Combination Operation (COCOA), Direct Virtual Coil (DVC), Synthetic Target (ST), parallel reconstruction Based on Successive Convolution Operations (BOSCO), or the like, have been developed to improve the speed of parallel imaging (PI) and partially parallel imaging (PPI) techniques. The techniques involve replacing the spatial encoding information, typically determined by the gradient controller 22 with the gradient coils 16, with spatial information contained in the RF coil array, e.g. characteristics of each coil element 18 such as position, orientation, size, and shape are incorporated into the received MR signals.
In partially parallel imaging techniques only a fraction of the phase encodes lines acquired compared to conventional parallel acquisitions. One of the mentioned reconstruction techniques is the applied to the acquired dataset to reconstruct the missing data, resulting in a full k-space dataset and thus a full FOV image. A problem with these PI and PPI reconstruction techniques and others is that as the number of RF receiver 26 channels increases so does the reconstruction time. In a single channel technique such as k-t FOCUSS the reconstruction time increases linearly with each channel count. For channel-by-channel techniques such as GRAPPA and COCOA, the reconstruction time increases quadratically. For example, a 32 channel coil array may take up to 16 times longer than an 8 channel coil array.
The magnetic resonance data processor 30 includes one or more processors which are programmed to perform the method described in conjunction with
With reference to
A k-t FOCUSS or GRAPPA or other partially parallel reconstruction technique is implemented by a multi-channel reconstruction unit, routine, or means 70. The parallel reconstruction unit, routine, or means 70 includes a unit, routine, or means 72 which combines, interpolates, and extrapolates the data from the M source channels to create a full k-space data set. The full k-space data set is reconstructed by a reconstruction unit, routine, or means 74, into M′ full images 76. M′ can be the same as M or, for a two-stage compression, M′ can be less than M. For example, 32 original receive channels can be compressed into 12 source channels and reconstructed into 8 channels or 8 images to reduce data processing time. For example, the number of k-space data lines and/or the number of data points in each k-space data line can be reduced to reduce the number of mathematical operations.
A calculation unit, routine, or means 80 generates relative sensitivity maps for the N′ combined data channels. The relative sensitivity maps can be calculated from sensitivity maps previously generated during setup for each coil element using a DC component of each sensitivity map as a reference point to determine the N relative sensitivity maps 82. In an alternate embodiment, N relative sensitivity maps are generated from the original N channels.
A Fourier-transform/projector unit, routine, or means 90 uses the N′ relative sensitivity maps to Fourier-transform of the M′ full images pointwise multiplied by each of the N′ relative sensitivity maps 82 to generate N′ recreated k-space data sets 92, corresponding to the N′ channels. A combining unit, routine, or means 94 combines each of the N′ combined data sets or channels with a corresponding one of the recreated N′ data sets. More specifically, in each of the recreated k-space data sets 92, data values which are actually available from the corresponding channel are substituted for the corresponding recreated data value. In some instances, the recreated k-space data sets may be missing. In this manner, N′ substituted full k-space data sets 96 are generated.
A reconstruction unit, routine, or means 100 reconstructs each of the N′ substituted k-space data sets to generate N′ partial images. An image combining unit, routine, or means 104 combines the N′ k-space images to generate a full image 106.
With reference to
A unit, routine, or means 110 chooses P channels of the combined channels with the largest eigenvalues, i.e. the channels with the majority of the motion artifacts, to generate P source channels 112 with large eigenvalues. A COCOA unit, routine, or means 114 applies the data convolution or combination operation motion-compensation technique using channel-by-channel k-space convolution to generate P motion-corrected k-space data sets 116. In the meantime, the N′-P channels 118 with the smaller eigenvectors are less motion affected and have unreliable relative sensitivity maps, and therefore are not corrected for motion. A combining unit, routine, or means 120 combines the P motion-corrected channels 116 with the N′-P channels with the smaller eigenvectors to generate N′ channels 122. These combined channels with the lowest eigenvalues can be used directly in the final reconstruction (100) without extra processing (120)
A unit, routine, or means 124 performs the method of
The invention has been described with reference to the preferred embodiments. Modifications and alterations may occur to others upon reading and understanding the preceding detailed description. It is intended that the invention be constructed as including all such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
| Filing Document | Filing Date | Country | Kind | 371c Date |
|---|---|---|---|---|
| PCT/IB2012/051265 | 3/16/2012 | WO | 00 | 11/26/2013 |
| Number | Date | Country | |
|---|---|---|---|
| 61453564 | Mar 2011 | US |
