Correcting residual aliasing in accelerated magnetic resonance imaging

Patent metadata

Field	Value
Publication number	US-10852382-B2
Application number	US-201816022348-A
Country	US
Kind code	B2
Filing date	Jun 28, 2018
Priority date	Jun 28, 2018
Publication date	Dec 1, 2020
Grant date	Dec 1, 2020

Abstract

Magnetic resonance imaging (MRI) systems and methods to determine and/or correct slice leakage and/or residual aliasing in the image domain in accelerated MRI imaging. Some implementations process one slice of MRI image domain data by input to a sensitivity encoding (SENSE) un-aliasing matrix built from predetermined RF signal reception sensitivity maps, thereby producing SENSE-decoded MRI image domain data for one pass-through image slice and at least one extra slice, and determine inter-slice leakage and/or in-plane residual aliasing based on content of the at least one extra output slice from the SENSE-decoded MRI image domain data. Some implementations correct slice leakage in reconstructed images by generating a fractional leakage matrix of inter-slice leakage measurements, and by multiplying the inverted fractional leakage matrix with uncorrected reconstructed images.

First claim

What is claimed is:

1. A magnetic resonance imaging (MRI) image domain method for processing MRI images reconstructed from k-space MRI data acquired from multiple magnetic resonance (MR) radio frequency (RF) signal receiving coils having respectively associated known RF signal reception sensitivity fields, MR signal contributions from each slice having been de-multiplexed from a plurality of multiplexed slices or aliased locations in reconstructed image domain MRI data to be processed, said method comprising:
forming, in a digital memory connected to at least one processor, a linear transformation representing said MR signal contributions from each multiplexed slice or aliased location into other slices or locations;
inverting the linear transformation formed in the digital memory; and
generating, using the at least one processor, one or more improved SENSE reconstructed images of respective ones of said plurality of slices or locations having reduced MR signal contributions between slices or locations by a composition of the one or more SENSE reconstructed MRI images and the inverted linear transformation; and
outputting the generated one or more improved reconstructed images to a display or digital storage.

2. The MRI image domain processing method as in claim 1 wherein the linear transformation is a matrix of complex-valued inter-slice signal leakage measurements.

3. The MRI image domain processing method as in claim 2 , wherein the inter-slice leakage measurements in the linear transformation are column-wise normalized prior to the inverting such that respective values on a diagonal of the matrix are each equal to or near to 1.

4. The MRI image domain processing method as in claim 2 , wherein the inter-slice leakage measurements in the linear transformation are column-wise normalized prior to the inverting by dividing values of a respective column with a square root of sum of squares of all values in the column.

5. The MRI image domain processing method as in claim 2 , wherein the inter-slice leakage measurements in the linear transformation are obtained from MRI image domain data corresponding to respective single-slice excited acquisitions.

6. The MRI image domain processing method as in claim 5 , wherein the single-slice excited acquisitions are performed in an MRI pre-scan data acquisition sequence prior to a main MRI diagnostic scan data acquisition sequence or as part of the main MRI diagnostic scan data acquisition sequence.

7. The MRI image domain processing method as in claim 1 , wherein the linear transformation comprises a matrix of elements such that non-zero values of off-diagonal elements of the matrix correspond to inter-slice leakage measurements and values of diagonal elements of the matrix correspond to shading.

8. The MRI image domain processing method as in claim 1 , wherein the generating comprises multiplying one of the reconstructed MRI images by the inverted linear transformation.

9. The MRI image domain processing method as in claim 1 , wherein the reconstructed MRI images are generated by unfolding the k-space MRI data using a sensitivity matrix corresponding to said RF signal receiving coils.

10. The MRI image domain processing method as in claim 1 , further comprising:
accessing, prior to said demultiplexing, MRI image domain data nominally representing at least one slice of a region of interest (ROI);
processing said one slice of accessed MRI image domain data by input to a sensitivity encoding (SENSE) un-aliasing matrix built from predetermined RF signal reception sensitivity maps for multiple RF MR signal receiving coils used to acquire MRI data to be processed, thereby producing as matrix output SENSE-decoded MRI image domain data for one pass-through image slice and at least one extra slice; and
determining inter-slice leakage and/or in-plane residual aliasing based on content of said at least one extra output slice from the SENSE-decoded MRI image domain data,
wherein the linear transformation is formed based on the determined inter-slice leakage and/or in-plane residual aliasing.

11. The MRI image domain processing method as in claim 10 wherein:
said un-aliasing matrix is set up to un-alias MR images acquired by a simultaneous multi-slice (SMS) MRI data acquisition sequence so that when the MR slice image data is input to the matrix at least one extra slice image is produced, and
slice leakage from an input slice to another slice of the ROI is determined based on content of said at least one extra slice image.

12. A magnetic resonance imaging (MRI) apparatus comprising:
static and gradient magnetic field generators and a plurality of radio frequency (RF) coils disposed within said magnetic field generators, the RF coils having respectively associated known RF signal reception sensitivity fields;
RF transmitter and RF receiver circuits coupled to said RF coils; and
an MRI sequence controller including a k-space memory and at least one processor configured to control said magnetic field generators and said RF transmitter and receiver circuits to:
reconstruct MRI image domain data from k-space MRI data acquired from said RF coils, MR signal contributions from one slice having been de-multiplexed from a plurality of multiplexed slices or aliased locations in reconstructed image domain MRI data;
form, in a digital memory connected to the at least one processor, a linear transformation representing said MR signal contributions from each multiplexed slice or aliased location into other slices or locations;
invert the linear transformation formed in the digital memory;
generate, using the at least one processor, one or more improved SENSE reconstructed images of respective ones of said plurality of slices or locations having reduced MR signal contributions between slices or locations by a composition of the one or more SENSE reconstructed MRI images and the inverted linear transformation; and
output the generated one or more improved reconstructed images to a display or digital storage.

13. The MRI apparatus according to claim 12 , wherein said at least one processor is further configured to form the linear transformation to include a matrix of complex-valued inter-slice signal leakage measurements.

14. The MRI apparatus according to claim 13 , wherein said at least one processor is further configured to column-wise normalize the inter-slice leakage measurements in the linear transformation prior to the inverting such that respective values on a diagonal of the matrix are each equal to or near to 1.

15. The MRI apparatus according to claim 13 , wherein said at least one processor is further configured to column-wise normalize the inter-slice leakage measurements in the linear transformation prior to the inverting by dividing values of a respective column with a square root of sum of squares of all values in the column.

16. The MRI apparatus according to claim 12 , wherein said at least one processor is further configured to obtain the inter-slice leakage measurements in the linear transformation from MRI image domain data corresponding to respective single-slice excited acquisitions.

17. The MRI apparatus according to claim 16 , wherein said at least one processor is further configured to perform the single-slice excited acquisitions in an MRI pre-scan data acquisition sequence prior to a main MRI diagnostic scan data acquisition sequence or as part of the main MRI diagnostic scan data acquisition sequence.

18. The MRI apparatus according to claim 12 , wherein said at least one processor is further configured to g

Assignees

• Canon Medical Systems Corp

Inventors

• Dannels Wayne R
• Sharma Anuj

Classifications

• G06T12/20

  Inverse problem, i.e. transformations from projection space into object space · CPC title

• G01R33/4835 Primary

  of multiple slices · CPC title

• G06T2207/10088

  Magnetic resonance imaging [MRI] · CPC title

• G06T7/0012

  Biomedical image inspection · CPC title

• G01R33/56545 Primary

  caused by finite or discrete sampling, e.g. Gibbs ringing, truncation artefacts, phase aliasing artefacts · CPC title