Real-Time Compressive Sensing MRI Reconstruction Using GPU Computing and Split Bregman Methods

Algorithm 1

GPU-based Split Bregman Compressive Sensing Reconstruction.

(1) Transfer Fourier data and k-space mask to GPU.

(2) Initialize necessary temporary arrays on the GPU.

(3) Precompute outer and inner loop-independent factors.

outer loop 𝑛  times:

 inner loop 𝑚 times:

    (4) Minimize L2 norm by direct matrix inversion.

    (5) Update solution image.

    (6) “Shrink” TV of solution.

    (7) Update Bregman parameters.

 end

  (8) Update Fourier data matrix.

  (9) Update inner loop-independent matrices.

end

(10) Transfer solution to CPU.