Next Article in Journal
Improved Image Quality for Static BLADE Magnetic Resonance Imaging Using the Total-Variation Regularized Least Absolute Deviation Solver
Previous Article in Journal
Non-Hemorrhagic Adrenal Infarction during Pregnancy: The Diagnostic Imaging Keys
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simultaneous Multislice Brain MRI T1 Mapping with Improved Low-Rank Modeling

1
Siemens Healthineers Korea Ltd., Seoul 03737, Korea
2
Department of Brain and Cognitive Engineering, Korea University, Seoul 02841, Korea
3
Department of Computer Engineering, Chonnam National University, Gwangju 61186, Korea
4
Department of ICT Convergence System Engineering, Chonnam National University, Gwangju 61186, Korea
*
Author to whom correspondence should be addressed.
Submission received: 13 August 2021 / Revised: 20 September 2021 / Accepted: 29 September 2021 / Published: 7 October 2021

Abstract

:
To accelerate data acquisition speed in magnetic resonance imaging (MRI), multiple slices are simultaneously acquired using multiband pulses. Simultaneous multislice (SMS) imaging typically unfolds slice aliasing from the acquired collapsed slices. In this study, we extended the SMS framework to accelerated MR parameter quantification such as T1 mapping. Assuming that the slice-specific null space and signal subspace are invariant along the parameter dimension, we formulated the SMS framework as a constrained optimization problem under a joint reconstruction framework such that the noise and signal subspaces are used for slice separation and recovery, respectively. The proposed method was validated on 3T MR human brain scans. We successfully demonstrated that the proposed method outperforms competing methods in suppressing aliasing artifacts and noise at high SMS accelerations, thus leading to accurate T1 maps.

1. Introduction

In magnetic resonance imaging (MRI), quantitative parameter mapping, which includes voxelwise delineation of tissue-specific relaxation times, has been widely used in characterizing inherent tissue properties and evaluating various pathological diseases, thereby providing valuable insight into disease processes. These voxels serve as imaging markers for clinical applications such as acute stroke, epilepsy, and multiple sclerosis [1,2,3,4]. As parameter mapping in MRI requires that several repeated measurements with varying imaging parameters, it may result in prohibitively long imaging times, thus compromising imaging efficiency and possibly limiting a range of clinical applications.
In accelerated MR parameter mapping, simultaneous multislice (SMS) imaging has gained attention owing to its ability to provide the signal-to-noise ratio (SNR) benefit of volumetric signal averaging by simultaneously exciting multiple imaging slices [5,6,7]. Parallel imaging techniques have been used to unfold the overlapped slices by solving an inverse problem consisting of multicoil linear system equations [8,9,10,11,12]. To obtain better noise suppression and artifact mitigation over conventional data acquisition, a controlled aliasing (CAIPI) technique was introduced to SMS acquisition by applying phase-modulated RF excitation to best utilize the multicoil encoding power that maximizes the sensitivity variation between neighboring slices [12,13]. CAIPI acquisition, coupled with parallel imaging, provides higher SNR efficiency due to the improved conditioning of the matrix inversion. To flexibly control intra- and interslice artifacts, split slice-GRAPPA (SP-SG) imposes a slice leakage constraint by effectively passing signals in a slice of interest while nulling the leakage terms of the other slices [14,15,16]. However, as simultaneously acquired slices are close to each other, the neighboring slices may share spatial correlation to a certain extent, thus limiting the use of a single reconstruction kernel.
To enable fast MR parameter mapping for a large number of slices, an extended SMS Hankel subspace learning (HSL) framework, which is the generalization of SMS-HSL [17] that exploits the Hankel structured matrix property in k-space, was developed in this study. For slice unfolding, the null subspace of the acquired SMS data is combined with the signal subspace to recover the signal within a single slice. To this end, the parameter dimension was considered for dynamic SMS imaging by combining the Hankel matrices from different contrast images to yield a larger but more redundant low-rank matrix model that consists of signal and noise subspaces. In vivo experiments were performed using SMS single-shot EPI data sets with multiple inversion recovery times for T1 mapping to validate the capability of the extended SMS-HSL algorithm at high SMS factors.

2. Materials and Methods

2.1. Introduction to SMS-HSL

An SMS signal can be represented as a linear combination of multiple slices:
y = s N s x s + n
where y N x × N y × N c is the zero-filled measured signal, y N x × N y × N c is the sth slice signal, and n N x × N y × N c is the measurements noise. To make the best of redundant information from intra- and intercorrelations in k-space, each slice signal can be rewritten as a Hankel matrix form [18]:
H ( y ) = s N s H ( x s ) + N
where H ( ) :   N x × N y × N c ( N x r + 1 ) ( N y r + 1 ) × r 2 N c the Hankel operator that maps the k-space data into a Hankel-structured matrix by vectorizing r × r rectangular k-space data across all coils inside the window followed by stacking the vectors row-wise as the sliding window shifts to generate different rows. In multiband excitation, the above signal can be decomposed on the basis of a dual slice model into 1. the slice of interest and 2. its complementary slices:
H ( y ) = H ( x s ) + H ( x s c ) + N
where x s c is the complementary slice signal and is defined as the linear superposition of all simultaneously excited slices excluding the slice of interest. The Hankel matrices in Equation (3) are then projected onto the different subspaces spanned by the complementary null space N s c .
H ( y ) N s c = H ( x s ) N s c + H ( x s c ) N s c H ( x s ) N s c
Accordingly, the aliasing separation in the slice direction is performed by solving the following least squares problem:
x s = min x s ( H ( y ) H ( x s ) ) N s c F 2
Additionally, as the multicoil images become redundant by sharing identical information across the coil dimension, the Hankel structured data matrix becomes highly rank-deficient, thus imposing a low-rank prior under the SMS reconstruction framework.
E ( x s ) = min x s 1 2 ( H ( y ) H ( x s ) ) N s c F 2 + λ H ( x s ) *
where * is the nuclear norm defined as the sum of the singular values of the matrix.

2.2. Extension to Parameter Dimension

The simplest approach to applying SMS-HSL to parameter mapping is to independently reconstruct each parameter image. As the null space is observed to be contrast-invariant along the parameter dimension, the null space is applied to the individual parameter images, thus leading to signal nulling in a slice of no interest while passing a slice of interest for all TIs (Figure 1). To make better low-rank characteristics in multicontrast imaging, the larger Hankel matrices H P ( ) :   N x × N y × N c × N p ( N x r + 1 ) ( N y r + 1 ) N p × r 2 N c , formulated as:
H P ( x s ) = [ H ( x s , 1 ) ;   H ( x s , 2 ) ; ;   H ( x s , N p ) ]
can be constructed by concatenating the Np different contrast Hankel-structured matrices in the parameter dimension. The Hankel-structed low- rank properties can be observed by applying singular value decomposition (SVD) to a single slice, linearly combined composite SMS slice, and multicontrast block matrices, respectively. Figure 2 shows the corresponding singular value distributions. The block matrix shows that singular values drop quite rapidly compared to different matrices. This implies that the conversion of the individual parameter matrices into a larger multicontrast matrix increases the rank deficiency, thus making the matrix completion concept of SMS-HSL more suitable to MR parameter mapping in SMS imaging. This implies that the conversion of the individual parameter matrices into a larger multiparameter matrix increases the rank deficiency, thus making the matrix completion concept of SMS-HSL more suitable to MR parameter mapping in SMS imaging. In addition, the larger Hankel matrices potentially have a larger number of null space vectors as a result of the increased low rank of the Hankel structured matrix. This suggests that the null space projection term in Equation (6) improves the performance while unfolding the collapsed slices into individual slices by utilizing the increased redundancy of the null space. As a result, we computed the null space vectors using larger Hankel matrices obtained by:
(1)
concatenating the individual parameter Hankel structured matrices;
(2)
applying SVD to block Hankel matrix;
(3)
taking the right singular vectors corresponding to the low singular values.
Similar to the SMS-HSL based reconstruction formula of Equation (7), the proposed reconstruction problem can be formulated as:
E ( x s ) = min x s p = 1 N p 1 2 ( H ( y p ) H ( x s , p ) ) N s c F 2 + λ H P ( x s ) *
The first term measures how well the collapsed slices are separated into individual slices, and the second regularization term provides a powerful signal recovery constraint for each slice by converting the slice separation problem into a single slice signal recovery problem. For the sake of brevity, Equation (8) is represented by concatenating the Hankel structured matrix for each contrast:
E ( x s ) = min x s ( H P ( y ) H P ( x s ) ) N s c F 2 + λ H P ( x s ) *
where H P ( y p )   and   H P ( x s ) are defined as:
H P ( y ) = [ H ( y 1 ) ;   H ( y 2 ) ; ;   H ( y N P ) ] ,   H P ( x s ) = [ H ( x s , 1 ) ;   H ( x s , 2 ) ; ;   H ( x s , N p ) ]
As the proposed reconstruction framework of Equation (8) is highly similar to the SMS-HSL cost function, the same optimization algorithm can be used to minimize both terms. In this study, the low-rank matrix fitting approach (LMaFit) was used under the framework of an alternating direction method (ADM) [9] to solve the aforementioned formulation.

2.3. Optimization Algorithm

Increasing the number of parameter dimensions makes Equation (8) computationally intractable. Instead, we employ the LMaFit approach, which does not use SVD, to substantially mitigate the computational complexity. The algorithm factorizes the larger Hankel-structured matrices into spatial coefficients and coil basis functions.
H P ( x s ) * = min U S ,   V S 1 2 ( U S F 2 + V S F 2 )
Hence, using the nuclear norm minimization under the matrix factorization constraint, Equation (9) can be rewritten as:
E ( x s , U S , V S ) = min x s ,   U S , V S ( H P ( y ) H P ( x s ) ) N s c F 2 + λ 2 ( U S F 2 + V S F 2 )
subject to H P ( x s ) = U S V S .
In Equation (11), the object function is nonconvex in terms of the unknowns x s ,   U S , and V S , and thus the global optimal solution is not guaranteed. However, the object function is convex with respect to one variable when the other variables are held constant. To achieve nonconvex optimization, we employ an ADM algorithm to update x s ,   U S , and V S in an alternating fashion, to iteratively find suboptimal solutions, until the cost function stops decreasing [19,20,21].
Minimization with respect to x s : The cost function in Equation (12) is reduced to:
E ( x s ) = min x s ( H P ( y ) H P ( x s ) ) N s c F 2 + α 2 U S V S H P ( x s ) F 2
where α is a penalty parameter that balances slice separation and signal recovery.
Minimization with respect to U S : The cost function in Equation (12) is reduced to:
E ( U S ) = min U S λ 2 U S F 2 + α 2 U S V S H P ( x s ) F 2
Minimization with respect to V S : The cost function in Equation (11) is reduced to:
E ( V S ) = min V S λ 2 V S F 2 + α 2 U S V S H P ( x s ) F 2
Each subproblem in Equations (13)–(15) is convex, and a global optimal solution can thus be found for each subproblem. The quadratic subproblems are all solved in an alternating fashion using the nonlinear conjugate gradient (CG) algorithm.

2.4. Experimental Setup

Experimental studies were performed on a 3T PRISMA scanner (Siemens Healthineers, Erlangen, Germany) equipped with a 20-channel head coil using an inversion recovery (IR) gradient echo (GE) EPI sequence. All experimental procedures were performed under the approval of the Institutional Review Board at the Sungkyunkwan University (Suwon, Korea). An informed written consent was obtained from each volunteer prior to imaging. Fully sampled brain data were acquired with the following imaging parameters: TR/TE = 3 s/23 ms, FOV = 240 × 240, matrix size = 120 × 120, slice thickness = 2 mm, number of slices = 25, and 15 IR scans (inversion times TI = 50, 250, 450, 650, 850, 1050, 1250, 1450, 1650, 1850, 2050, 2250, 2450, 2650, and 2850 ms). All SMS image reconstructions were performed offline on a personal computer with 2.3 GHz CPU and 32 GB RAM using MATLAB (Mathworks Inc., Natick, MA, USA). For faster reconstruction, the acquired data sets were reduced to 12 channels using standard SVD coil compression techniques [22,23]. Prior to data acquisition, a FOV-matched low resolution 2D GRE image was also acquired with 32 fully sampled phase encoding lines at the center of k-space for calibration. To retrospectively emulate SMS data acquisition, slice-specific CAIPI-induced phase was added to each slice in k-space, before combining the slices into an overlapped slice. For reconstruction, a kernel size of 5 × 5 (phase encoding × readout) was used in the Hankel-structured matrix, and the low-rank property was imposed on the Hankel block matrix with a total number of 120 spatial coefficients and coil basis, respectively. The cut-off singular value was set toe 0.05 max for the null space selection. The performance of the proposed method depends on the choice of the regularization parameters λ and α . Here, we manually set ( λ , α ) to small values ( 1.0 e 7 ,   2.0 e 6 ) to minimize slice leakage artifacts without perturbing measured SMS data. For visual evaluations of image quality, error maps were produced by calculating the difference between the reference and estimated images. Quantitative comparison was also provided in terms of the normalized root mean square error (nRMSE): 1 max ( I s ) s = 1 N 1 N ( I s I ^ s ) 2 , where I s and I ^ s are reference and estimated images, respectively. N is the total number of pixels. To demonstrate the benefits of imposing multiparametric null space and low rank, we also performed kernel-based and null space-based SMS reconstructions (i.e., SP-SG and SMS-HSL). The null space rank order of SMS-HSL was the same as that of the proposed method. For SP-SG, kernel estimation was regularized with Tikhonov penalty, and kernel sizes and regularization parameter were selected by visual inspection.
S I = S I 0 [ 1 2 e ( T I T 1 ) + e ( T R T 1 ) ]
The term including TR could be negligible (TR > 5 × T1). However, it needs a long scan time because the T1 values at 3T are around 1000–2000 ms. Thus, in this work, we employed the above signal model including TR as a fitting model to estimate the exact T1 values of brain tissues. After reconstruction, T1 map were estimated using pixelwise nonlinear least square fitting.

3. Results

Figure 3 shows reconstructed images, error maps, and nRMSEs for SP-SG, SMS-HSL, and the proposed method at an SMS factor of 5 for two of the TIs (50 and 1250 ms) out of a total of 15 IR scans. The SP-SG reconstructed images exhibit significant noise amplification in the third image and some structured aliasing artifacts in the fourth image with TI = 50 ms and similarly with TI = 1250 ms (red arrows), in the middle of the FOV. The visual quality of the SP-SG breaks down at longer TIs due to the low SNR in signal nulling. The SMS-HSL case does not display strong artifacts, although close inspection does reveal that some noise is still present. This is to be expected due to the improved conditioning of the matrix inversion through generalized null space kernels, which are all independent of each other. Notably, our proposed method consistently reconstructs high-quality images that are devoid of visible artifacts for two of the TIs, attributed to the increased low rank of the Hankel-structured matrices of multi-TI data. Consistent with the results shown in the reconstructed images, the corresponding error maps tend to exaggerate contrast contamination with noises particularly in SP-SG while yielding the lowest nRMSEs in the proposed method compared to SP-SG and SMS-HSL.
Figure 4 compares the corresponding T1 maps (Figure S2) with enlarged figure in the white dotted boxes, calculated from the reconstructed images of multi-TI data sets at an SMS factor of 5. In line with the reconstruction results, the T1 values that are estimated from SP-SG deviate from the correct values, even for the same tissue, due to the severely amplified noise. In contrast, SMS-HSL provides dramatic improvement over SP-SG despite the presence of some noise. Finally, with the proposed method, the T1 maps have better SNR and the T1 values have higher precision than those of SP-SG and SMS-HSL, consistent with the reconstruction results.
For quantitative evaluation, T1 relaxation curves are shown in Figure 5. The values within manually chosen 3 × 3 pixel white and gray matter regions of interest (ROI) were averaged. The T1 values for white and gray matter follow the pattern of the aforementioned results. The proposed method follows the true signal intensity curve much more accurately than the other methods, which display some signal bias due to signal nulling during inversion recovery. The disparity between the proposed method and the other methods is more pronounced for white matter.

4. Discussion

In this work, we extend SMS-HSL to the context of dynamic MR imaging. This is achieved by concatenating the SMS-HSL matrices of different contrasts to yield a larger, more redundant low-rank dynamic imaging data matrix. This new reconstruction model utilizes the null-space-based slice unfolding followed by enforcing the low-rankness of the block Hankel matrix independently for each slice. The improved low rank is proof of the existence of null space filters in the k-space domain whereby k-space data are represented in a lower dimensional signal subspace, enabling the unfolding of each slice. Following up on our previous work on SMS reconstruction [17], we developed a new algorithm to incorporate the null and signal subspaces into a single optimization problem under a joint framework. We demonstrate the superior performance of the proposed algorithm in the accelerated recovery of SMS MR parameter data, using in vivo MR brain scans.
T1 mapping commonly relies on the inversion recovery sequence where EPI acquisitions are continuously applied during signal recovery followed by a computation of the T1 values during postprocessing. Multishot EPI acquisitions potentially yield motion mismatches and shot-to-shot phase variations between scans. These issues could be addressed by adopting SMS acceleration, which brings down the time to enable a single-shot SMS acquisition. The advanced technique presented in this work provides a drastic improvement in image quality at high SMS accelerations, as observed in Figure 3, while helping mitigate some of the involuntary motion and phase variation issues.
In terms of computation, the SVD-free reconstruction algorithm (LMaFit), which is based on matrix factorization, was employed. The proposed algorithm is still computationally demanding due to the existence of a parameter dimension. To mitigate this computational issue, the matrix inversion of the proposed method in Equation (7) can be precomputed. In addition, the proposed algorithm can be efficiently parallelized using a combination of multicore CPUs with GPUs, such that the null space estimation, which requires high computational cost, is performed on multicore CPUs for each slice while the matrix multiplication for data synthesis is implemented on GPUs with a large number of cores.
Our proposed method is not restricted to T1 quantification. It can be easily applied to the mapping of different MR parameters such as T2, T2 *, and diffusion with acquisitions based not only on single-shot EPI sequences but also on fast spin echo (FSE) or gradient echo (GRE) sequences.

5. Conclusions

In conclusion, we developed and evaluated a dynamic version of SMS-HSL in the reconstruction of SMS imaging data by simultaneously imposing noise and signal subspaces to the calculation of T1 maps. The proposed SMS-HSL method was demonstrated to have substantial advantages over competing methods in providing high quality reconstructions and high precision T1 maps with up to fivefold SMS acceleration.

Supplementary Materials

The proposed method was additionally evaluated on human brain from different data set of the healthy subject available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/tomography7040047/s1, Figure S1: reconstructed images using SP-SG, SMS-HSL and proposed method, Figure S2: the corresponding T1 maps.

Author Contributions

Conceptualization, S.K. and S.P.; methodology, S.P.; software, S.K.; validation, S.K.; writing—original draft preparation, S.K. and S.P.; supervision, S.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2021R1C1C1013603).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of Sungkyunkwan University (2018-05-005).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Additional data are available from the first/corresponding author upon reasonable request.

Acknowledgments

The authors acknowledge the facilities and the technical assistance at the Center for Neuroscience Imaging Research (CNIR), Sungkyunkwan University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Warntjes, J.B.; Leinhard, O.D.; West, J.; Lundberg, P. Rapid magnetic resonance quantification on the brain: Optimization for clinical usage. Magn. Reson. Med. 2008, 60, 320–329. [Google Scholar] [CrossRef]
  2. Alexander, A.; Hurley, S.; Samsonov, A.; Adluru, N.; Hosseinbor, A.; Mossahebi, P. Characterization of cerebral white matter properties using quantitative magnetic resonance imaging stains. Brain Connect. 2011, 1, 423–446. [Google Scholar] [CrossRef]
  3. Cheng, H.; Stikov, N.; Ghugre, N.; Wright, G. Practical medical applications of quantitative MR relaxometry. J. Magn. Reson. Imaging 2012, 36, 805–824. [Google Scholar] [CrossRef] [PubMed]
  4. Komiyama, M.; Yagura, H.; Baba, M.; Yasui, T.; Hakuba, A.; Nishimura, S. MR imaging: Possibility of tissue characterization of brain tumors using T1 and T2 values. Am. J. Neuroradiol. 1987, 8, 65–70. [Google Scholar] [PubMed]
  5. Larkman, D.; Hajnal, J.; Herlihy, A.; Coutts, G.; Young, I.; Ehnholm, G. Use of multicoil arrays for separation of signal from multiple slices simultaneously excited. J. Magn. Reson. Imaging 2001, 13, 313–317. [Google Scholar] [CrossRef]
  6. Moeller, S.; Yacoub, E.; Olman, C.; Auerbach, E.; Strupp, J.; Harel, N.; Uğurbil, K. Multiband multislice GE-EPI at 7 tesla, with 16-fold acceleration using partial parallel imaging with application to high spatial and temporal whole-brain fMRI. Magn. Reson. Med. 2010, 63, 1144–1153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Feinberg, D.; Moeller, S.; Smith, S.; Auerbach, E.; Ramanna, S.; Gunther, M.; Glasser, M.; Miller, K.; Ugurbil, K.; Yacoub, E. Multiplexed echo planar imaging for sub-second whole brain fMRI and fast diffusion imaging. PLoS ONE 2010, 5, e15710. [Google Scholar] [CrossRef] [Green Version]
  8. Pruessmann, K.; Weiger, M.; Scheidegger, M.; Boesiger, P. SENSE: Sensitivity encoding for fast MRI, Magn. Reson. Med. 1999, 42, 952–962. [Google Scholar] [CrossRef]
  9. Griswold, M.; Jakob, P.; Heidemann, R.; Nittka, M.; Jellus, V.; Wang, J.; Kiefer, B.; Haase, A. Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn. Reson. Med. 2002, 47, 1202–1210. [Google Scholar] [CrossRef] [Green Version]
  10. Zhu, K.; Kerr, A.; Pauly, J. Autocalibrating CAIPIRINHA: Reformulating CAIPIRINHA as a 3D problem. In Proceedings of the 20rd Annual Meeting of ISMRM, Melbourne, Australia, 5–11 May 2012. [Google Scholar]
  11. Zahneisen, B.; Poser, B.A.; Ernst, T.; Stenger, V.A. Three-dimensional Fourier encoding of simultaneously excited slices: Generalized acquisition and reconstruction framework. Magn. Reson. Med. 2014, 71, 2071–2081. [Google Scholar] [CrossRef] [Green Version]
  12. Setsompop, K.; Gagoski, B.; Polimeni, J.; Witzel, T.; Wedeen, V.; Wald, L. Blippedcontrolled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g-factor penalty. Magn. Reson. Med. 2012, 67, 1210–1224. [Google Scholar] [CrossRef] [Green Version]
  13. Breuer, F.; Blaimer, M.; Heidemann, R.; Mueller, M.; Griswold, M.; Jakob, P. Controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA) for multi-slice imaging. Magn. Reson. Med. 2005, 53, 684–691. [Google Scholar] [CrossRef]
  14. Cauley, S.; Polimeni, J.; Bhat, H.; Wald, L.; Setsompop, K. Interslice leakage artifact reduction technique for simultaneous multislice acquisitions. Magn. Reson. Med. 2014, 72, 93–102. [Google Scholar] [CrossRef] [Green Version]
  15. HashemizadehKolowri, S.K.; Chen, R.R.; Adluru, G.; Ying, L.; DiBella, E.V. Coil-combined split slice-GRAPPA for simultaneous multi-slice diffusion MRI. Magn. Reson. Imaging 2020, 66, 9–21. [Google Scholar] [CrossRef] [PubMed]
  16. HashemizadehKolowri, S.K.; Chen, R.R.; Adluru, G.; Dean, D.C.; Wilde, E.A.; Alexander, A.L.; DiBella, E.V. Simultaneous multi-slice image reconstruction using regularized image domain split slice-GRAPPA for diffusion MRI. Med. Image Anal. 2021, 70, 102000. [Google Scholar] [CrossRef] [PubMed]
  17. Park, S.; Park, J. SMS-HSL: Simultaneous multislice aliasing separation exploiting Hankel subspace learning. Magn. Reson. Med. 2017, 78, 1392–1404. [Google Scholar] [CrossRef] [PubMed]
  18. Lee, D.; Jin, K.; Kim, Y.; Park, S.H.; Ye, J. Acceleration of MR parameter mapping using annihilating filter-based low rank Hankel matrix (ALOHA). Magn. Reson. Med. 2016, 76, 1848–1864. [Google Scholar] [CrossRef] [PubMed]
  19. Wen, Z.; Yin, W.; Zhang, Y. Solving a low-rank factorization model for matrix completion by a nonlinear successive over-relaxation algorithm. Math. Program. Comput. 2012, 4, 333–361. [Google Scholar] [CrossRef]
  20. Shen, Y.; Wen, Z.; Zhang, Y. Augmented Lagrangian alternating direction method for matrix separation based on low-rank factorization. Optim. Methods Softw. 2014, 29, 239–263. [Google Scholar] [CrossRef] [Green Version]
  21. Lingala, S.G.; Jacob, M. Blind compressive sensing dynamic MRI. IEEE Trans. Med. Imaging 2013, 32, 1132–1145. [Google Scholar] [CrossRef] [Green Version]
  22. Huang, F.; Vijayakumar, S.; Li, Y.; Hertel, S.; Duensing, G.R. A software channel compression technique for faster reconstruction with many channels. Magn. Reson. Imaging 2008, 26, 133–141. [Google Scholar] [CrossRef]
  23. Buehrer, M.; Pruessmann, K.P.; Boesiger, P.; Kozerke, S. Array compression for MRI with large coil arrays. Magn. Reson. Med. 2007, 57, 1131–1139. [Google Scholar] [CrossRef]
Figure 1. An illustration of the use of null space for slice separation. This work hypothesizes that the null space is invariant along the parameter dimension. The null space projection effectively filters out a slice of no interest while passing through a slice of interest regardless of contrast changes.
Figure 1. An illustration of the use of null space for slice separation. This work hypothesizes that the null space is invariant along the parameter dimension. The null space projection effectively filters out a slice of no interest while passing through a slice of interest regardless of contrast changes.
Tomography 07 00047 g001
Figure 2. Utility of low rank in Hankel structured matrix. The larger Hankel matrices, which combine all parameter k-space data, are more rank-deficient than single and composite slices, which implies that the matrix completion scheme is more suitable to the application of MR parameter mapping.
Figure 2. Utility of low rank in Hankel structured matrix. The larger Hankel matrices, which combine all parameter k-space data, are more rank-deficient than single and composite slices, which implies that the matrix completion scheme is more suitable to the application of MR parameter mapping.
Tomography 07 00047 g002
Figure 3. Images, error maps, and nRMSEs generated using SP-SG, SMS-HSL, and our proposed method at an SMS factor of 5 according to different contrasts (TIs = 50 ms and 1250 ms). SP-SG suffers from severe noise (third slice) and aliasing (fourth slice) due to the ill-conditioning of the inverse problem while SMS-HSL mitigates the artifacts with generalized null space kernels. Finally, our proposed method gives the best reconstruction results, owing to the lower rank of the Hankel-structured matrix. The visual quality looks similar for different Tis (Figure S1).
Figure 3. Images, error maps, and nRMSEs generated using SP-SG, SMS-HSL, and our proposed method at an SMS factor of 5 according to different contrasts (TIs = 50 ms and 1250 ms). SP-SG suffers from severe noise (third slice) and aliasing (fourth slice) due to the ill-conditioning of the inverse problem while SMS-HSL mitigates the artifacts with generalized null space kernels. Finally, our proposed method gives the best reconstruction results, owing to the lower rank of the Hankel-structured matrix. The visual quality looks similar for different Tis (Figure S1).
Tomography 07 00047 g003
Figure 4. Corresponding T1 maps calculated from SP-SG, SMS-HSL, and our proposed method at an SMS factor of 5. The results from SP-SG exhibit uneven T1 values in the same tissue due to the presence of severe noise while SMS-HSL results display improved conditioning with generalized null space kernels although some noise is still visible. The proposed method provides the best T1 values, with high precision, due to the improved low rank of the Hankel structured matrix.
Figure 4. Corresponding T1 maps calculated from SP-SG, SMS-HSL, and our proposed method at an SMS factor of 5. The results from SP-SG exhibit uneven T1 values in the same tissue due to the presence of severe noise while SMS-HSL results display improved conditioning with generalized null space kernels although some noise is still visible. The proposed method provides the best T1 values, with high precision, due to the improved low rank of the Hankel structured matrix.
Tomography 07 00047 g004
Figure 5. T1 signal intensity curves of white and gray matter regions. The curves are generated from the reference and the reconstructed images by averaging the voxels in each region of interest. The signal curves from the proposed method strictly follow the decay curves of the reference, SP-SG and SMS-HSL.
Figure 5. T1 signal intensity curves of white and gray matter regions. The curves are generated from the reference and the reconstructed images by averaging the voxels in each region of interest. The signal curves from the proposed method strictly follow the decay curves of the reference, SP-SG and SMS-HSL.
Tomography 07 00047 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, S.; Park, S. Simultaneous Multislice Brain MRI T1 Mapping with Improved Low-Rank Modeling. Tomography 2021, 7, 545-554. https://0-doi-org.brum.beds.ac.uk/10.3390/tomography7040047

AMA Style

Kim S, Park S. Simultaneous Multislice Brain MRI T1 Mapping with Improved Low-Rank Modeling. Tomography. 2021; 7(4):545-554. https://0-doi-org.brum.beds.ac.uk/10.3390/tomography7040047

Chicago/Turabian Style

Kim, Sugil, and Suhyung Park. 2021. "Simultaneous Multislice Brain MRI T1 Mapping with Improved Low-Rank Modeling" Tomography 7, no. 4: 545-554. https://0-doi-org.brum.beds.ac.uk/10.3390/tomography7040047

Article Metrics

Back to TopTop