System and method for magnetic resonance imaging acceleration

Patent metadata

Field	Value
Publication number	US-12442879-B2
Application number	US-202217662649-A
Country	US
Kind code	B2
Filing date	May 9, 2022
Priority date	Jun 29, 2017
Publication date	Oct 14, 2025
Grant date	Oct 14, 2025

Abstract

Systems and methods for magnetic resonance imaging acceleration. The systems may perform the methods to obtain imaging data of a subject, or a portion thereof, captured by the MRI system according to an undersampling pattern; execute a first iterative procedure; determine that the first iteration number meets a first threshold; execute, in response to the determination that the first iteration number meets the first threshold, a second iterative procedure; determine that a sum of the first iteration number and the second iteration number meets a second threshold; and generate a reconstruction image of the subject, or a portion thereof, according to the processed imaging data, wherein the first threshold is lower than the second threshold.

First claim

We claim:

1. A system, comprising:
at least one non-transitory computer-readable storage medium including a set of instructions; and
at least one processor in communication with the at least one non-transitory computer-readable storage medium, wherein when executing the instructions, the system is directed to:
obtain imaging data of a subject, or a portion thereof, captured by an MRI system according to an undersampling pattern, the undersampling pattern comprises a region of a regular undersampling pattern and a region of an irregular undersampling pattern;
execute a first iterative procedure comprising:
performing a soft-thresholding de-noising operation on a sparsifying transform domain of the imaging data; and
determining a first iteration number of the soft-thresholding de-noising operation;

determine that the first iteration number meets a first threshold;
execute, in response to the determination that the first iteration number meets the first threshold, a second iterative procedure comprising:
performing a total variation algorithm on the de-noised imaging data to provide processed imaging data; and
determining a second iteration number of the total variation operation;

determine that a sum of the first iteration number and the second iteration number meets a second threshold; and
generate a reconstruction image of the subject, or a portion thereof, according to a parallel imaging procedure coupled with a half-Fourier imaging procedure based on the processed imaging data,
wherein the first threshold is lower than the second threshold.

2. The system of claim 1 , wherein the system is further directed to:
optimize the total variation algorithm by a gradient descent operation.

3. The system of claim 1 , wherein the first iterative procedure further comprises:
filling back the acquired imaging data into a k-space obtained after de-noising operation; and
performing the de-noising operation on the sparsifying transform domain of filled-back imaging data.

4. The system of claim 1 , wherein the system is further directed to:
determine that the sum of the first iteration number and the second iteration number is lower than the second threshold;
perform a weighted synthesis operation on the processed imaging data to generate synthesized imaging data; and
perform the soft-thresholding de-noising operation on the sparsifying transform domain of the synthesized imaging data.

5. The system of claim 4 , wherein the system is further directed to:
divide a k-space into a first region and a second region;
assign the first region a full sampling pattern;
assign at least one part of the second region the irregular undersampling pattern;
determine an MRI scanning plan based on the assigned first region and the assigned second region; and
cause the MRI system to capture the imaging data according to the MRI scanning plan.

6. The system of claim 5 , wherein the system is further directed to:
divide the second region into a fifth region, a sixth region, a seventh region, and an eighth region,
wherein the first region is located in a center section of a plane of the k-space, the first region is defined by at least four sides including two sides in a first direction and two sides in a second direction different from the first direction, the fifth region is located along the two sides of the first region in the first direction, the sixth region is located along the two sides of the first region in the second direction, the sixth region is defined by at least two sides in the first direction, the seventh region is located along the two sides of the sixth region in the first direction, and
the first region has a larger sampling density than the fifth region, and the sixth region has a larger sampling density than the seventh region.

7. The system of claim 6 , wherein the irregular undersampling pattern includes a restrained undersampling pattern, wherein the restrained undersampling pattern includes a restraint region, and the restrained undersampling pattern includes:
determining a first unsampled point at a center section of the restraint region;
determining a second unsampled point adjacent to the first unsampled point along a preset first direction; and
determining a third unsampled point whose position satisfies a predetermined positional relationship with the first unsampled point, and the first unsampled point is aligned with the third unsampled point along a second preset direction different from the first preset direction;
sampling a plurality of remaining points in the restraint region with a certain probability.

8. The system of claim 7 , wherein the MRI scanning plan is related to the plane of the k-space, the k-space is a Fourier space of a three dimensional MRI image, and the plane of the k-space is based on at least one phase encoding direction of a pulse sequence.

9. The system of claim 8 , wherein the first region is located in the center section of the plane of the k-space, and the second region is located in a surrounding section of the plane of the k-space.

10. The system of claim 1 , wherein the system is further directed to:
perform a parallel imaging procedure on the imaging data.

11. The system of claim 1 , wherein the system is further directed to:
determine a reference phase of the imaging data;
perform a high-pass filtering operation on the imaging data to generate filtered imaging data;
determine that the sum of the first iteration number and the second iteration number meets the second threshold; and
perform, based on the reference phase, a phase correction operation on the filtered imaging data to generate phase corrected imaging data,

wherein the processed imaging data include the phase corrected imaging data.

12. The system of claim 11 , wherein the system is further directed to:
obtain the processed imaging data by performing a parallel imaging procedure on the phase corrected imaging data, wherein the parallel imaging procedure includes at least one of GRAPPA, SENSE, or SPIRIT.

13. The system of claim 12 , wherein the parallel imaging procedure is executed in an image domain.

14. The system of claim 11 , wherein the system is further directed to:
determine the reference phase based on a lowpass filtering operation and an inverse Fourier transform operation.

15. The system of claim 14 , wherein the lowpass filtering operation is executed by a symmetrical lowpass filter that is radially symmetric with respect to a central point of a k-space.

16. A system, comprising:
at least one non-transitory computer-readable storage medium including a set of instructions; and
at least one processor in communication with the at least one non-transitory computer-readable storage medium, wherein when executing the instructions, the system is directed to:
divide a k-space into a first region and a second region;
assign the first region a full sampling pattern;
assign at least one part of the second region an irregular undersampling pattern;
determine an MRI scanning plan based on the assigned first region and the assigned second region;
cause an MRI system to capture imaging data of a subject, or a portion thereof according to the MRI scanning plan;
obtain the imaging data captured by the MRI system according to an undersampling pattern, the undersampling pattern comprises a region of regular undersampling pattern and a region of irregular undersampling pattern;
execute a first iterative procedure comprising:
performing a soft-thresholding de-noising operation on a sparsifying transform domain of the imaging data; and

Assignees

• Shanghai United Imaging Healthcare Co Ltd

Inventors

• Li Guobin
• Zong Jinguang
• Li Zhaopeng

Classifications

• G01R33/561

  by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences · CPC title

• G01R33/5608 Primary

  Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels (image data processing or generation, in general G06T) · CPC title

• G06T5/70

  Denoising; Smoothing · CPC title

• G06T2207/10088

  Magnetic resonance imaging [MRI] · CPC title

• G01R33/5611

  Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE (structural details of arrays of sub-coils G01R33/3415) · CPC title