SIPAS: A comprehensive susceptibility imaging process and analysis studio

Quantitative susceptibility mapping (QSM) is a rising MRI-based technology and quite a few QSM-related algorithms have been proposed to reconstruct maps of tissue susceptibility distribution from phase images. In this paper, we develop a comprehensive susceptibility imaging process and analysis studio (SIPAS) that can accomplish reliable QSM processing and offer a standardized evaluation system. Specifically, SIPAS integrates multiple methods for each step, enabling users to select algorithm combinations according to data conditions, and QSM maps could be evaluated by two aspects, including image quality indicators within all voxels and region-of-interest (ROI) analysis. Through a sophisticated design of user-friendly interfaces, the results of each procedure are able to be exhibited in axial, coronal, and sagittal views in real-time, meanwhile ROIs can be displayed in 3D rendering visualization. The accuracy and compatibility of SIPAS are demonstrated by experiments on multiple in vivo human brain datasets acquired from 3T, 5T, and 7T MRI scanners of different manufacturers. We also validate the QSM maps obtained by various algorithm combinations in SIPAS, among which the combination of iRSHARP and SFCR achieves the best results on its evaluation system. SIPAS is a comprehensive, sophisticated, and reliable toolkit that may prompt the QSM application in scientific research and clinical practice.

Spinet-QSM: model-based deep learning with schatten p-norm regularization for improved quantitative susceptibility mapping

OBJECTIVE

Quantitative susceptibility mapping (QSM) provides an estimate of the magnetic susceptibility of tissue using magnetic resonance (MR) phase measurements. The tissue magnetic susceptibility (source) from the measured magnetic field distribution/local tissue field (effect) inherent in the MR phase images is estimated by numerically solving the inverse source-effect problem. This study aims to develop an effective model-based deep-learning framework to solve the inverse problem of QSM.

MATERIALS AND METHODS

This work proposes a Schatten p -norm-driven model-based deep learning framework for QSM with a learnable norm parameter p to adapt to the data. In contrast to other model-based architectures that enforce the l2 -norm or l1 -norm for the denoiser, the proposed approach can enforce any p -norm ( 0 < p ≤ 2 ) on a trainable regulariser.

RESULTS

The proposed method was compared with deep learning-based approaches, such as QSMnet, and model-based deep learning approaches, such as learned proximal convolutional neural network (LPCNN). Reconstructions performed using 77 imaging volumes with different acquisition protocols and clinical conditions, such as hemorrhage and multiple sclerosis, showed that the proposed approach outperformed existing state-of-the-art methods by a significant margin in terms of quantitative merits.

CONCLUSION

The proposed SpiNet-QSM showed a consistent improvement of at least 5% in terms of the high-frequency error norm (HFEN) and normalized root mean squared error (NRMSE) over other QSM reconstruction methods with limited training data.

Maximum spherical mean value filtering for whole-brain QSM

PURPOSE

To develop a tissue field-filtering algorithm, called maximum spherical mean value (mSMV), for reducing shadow artifacts in QSM of the brain without requiring brain-tissue erosion.

THEORY AND METHODS

Residual background field is a major source of shadow artifacts in QSM. The mSMV algorithm filters large field-magnitude values near the border, where the maximum value of the harmonic background field is located. The effectiveness of mSMV for artifact removal was evaluated by comparing existing QSM algorithms in numerical brain simulation as well as using in vivo human data acquired from 11 healthy volunteers and 93 patients.

RESULTS

Numerical simulation showed that mSMV reduces shadow artifacts and improves QSM accuracy. Better shadow reduction, as demonstrated by lower QSM variation in the gray matter and higher QSM image quality score, was also observed in healthy subjects and in patients with hemorrhages, stroke, and multiple sclerosis.

CONCLUSION

The mSMV algorithm allows QSM maps that are substantially equivalent to those obtained using SMV-filtered dipole inversion without eroding the volume of interest.

Assessment for Carotid Atherosclerotic Plaque Using Vessel Wall Magnetic Resonance Imaging: A Multireader ROC Study to Determine Optimal Sequence for Detecting Vessel Wall Calcification

INTRODUCTION

We aimed to compare conventional vessel wall MR imaging techniques and quantitative susceptibility mapping (QSM) to determine the optimal sequence for detecting carotid artery calcification.

METHODS

Twenty-two patients who underwent carotid vessel wall MR imaging and neck CT were enrolled. Four slices of 6-mm sections from the bilateral internal carotid bifurcation were subdivided into 4 segments according to clock position (0-3, 3-6, 6-9, and 9-12) and assessed for calcification. Two blinded radiologists independently reviewed a total of 704 segments and scored the likelihood of calcification using a 5-point scale on spin-echo imaging, FLASH, and QSM. The observer performance for detecting calcification was evaluated by a multireader, multiple-case receiver operating characteristic study. Weighted κ statistics were calculated to assess interobserver agreement.

RESULTS

QSM had a mean area under the receiver operating characteristic curve of 0.85, which was significantly higher than that of any other sequence (p < 0.01) and showed substantial interreader agreement (κ = 0.68). A segment with a score of 3-5 was defined as positive, and a segment with a score of 1-2 was defined as negative; the sensitivity and specificity of QSM were 0.75 and 0.87, respectively.

CONCLUSION

QSM was the most reliable MR sequence for the detection of plaque calcification.

Robust breast quantitative susceptibility mapping in the presence of silicone

PURPOSE

To (a) develop a preconditioned water-fat-silicone total field inversion (wfsTFI) algorithm that directly estimates the susceptibility map from complex multi-echo data in the breast in the presence of silicone and to (b) evaluate the performance of wfsTFI for breast quantitative susceptibility mapping (QSM) in silico and in vivo in comparison with formerly proposed methods.

METHODS

Numerical simulations and in vivo multi-echo gradient echo breast measurements were performed to compare wfsTFI to a previously proposed field map-based linear total field inversion algorithm (lTFI) with and without the consideration of the chemical shift of silicone in the field map estimation step. Specifically, a simulation based on an in vivo scan and data from five patients were included in the analysis.

RESULTS

In the simulation, wfsTFI is able to significantly decrease the normalized root mean square error from lTFI without (4.46) and with (1.77) the consideration of the chemical shift of silicone to 0.68. Both the in silico and in vivo wfsTFI susceptibility maps show reduced shadowing artifacts in local tissue adjacent to silicone, reduced streaking artifacts and no erroneous single voxels of diamagnetic susceptibility in proximity to silicone.

CONCLUSION

The proposed wfsTFI method can automatically distinguish between subjects with and without silicone. Furthermore wfsTFI accounts for the presence of silicone in the QSM dipole inversion and allows for the robust estimation of susceptibility in proximity to silicone breast implants and hence allows the visualization of structures that would otherwise be dominated by artifacts on susceptibility maps.

Data adaptive regularization with reference tissue constraints for liver quantitative susceptibility mapping

PURPOSE

To improve repeatability and reproducibility across acquisition parameters and reduce bias in quantitative susceptibility mapping (QSM) of the liver, through development of an optimized regularized reconstruction algorithm for abdominal QSM.

METHODS

An optimized approach to estimation of magnetic susceptibility distribution is formulated as a constrained reconstruction problem that incorporates estimates of the input data reliability and anatomical priors available from chemical shift-encoded imaging. The proposed data-adaptive method was evaluated with respect to bias, repeatability, and reproducibility in a patient population with a wide range of liver iron concentration (LIC). The proposed method was compared to the previously proposed and validated approach in liver QSM for two multi-echo spoiled gradient-recalled echo protocols with different acquisition parameters at 3T. Linear regression was used for evaluation of QSM methods against a reference FDA-approvedR 2 $$ {R}_2 $$ -based LIC measure andR 2 ∗ $$ {R}_2^{\ast } $$ measurements; repeatability/reproducibility were assessed by Bland-Altman analysis.

RESULTS

The data-adaptive method produced susceptibility maps with higher subjective quality due to reduced shading artifacts. For both acquisition protocols, higher linear correlation with bothR 2 $$ {R}_2 $$ - andR 2 ∗ $$ {R}_2^{\ast } $$ -based measurements were observed for the data-adaptive method (r 2 = 0 . 74 / 0 . 69 $$ {r}^2=0.74/0.69 $$ forR 2 $$ {R}_2 $$ ,0 . 97 / 0 . 95 $$ 0.97/0.95 $$ forR 2 ∗ $$ {R}_2^{\ast } $$ ) than the standard method (r 2 = 0 . 60 / 0 . 66 $$ {r}^2=0.60/0.66 $$ and0 . 79 / 0 . 88 $$ 0.79/0.88 $$ ). For both protocols, the data-adaptive method enabled better test-retest repeatability (repeatability coefficients 0.19/0.30 ppm for the data-adaptive method, 0.38/0.47 ppm for the standard method) and reproducibility across protocols (reproducibility coefficient 0.28 vs. 0.53ppm) than the standard method.

CONCLUSIONS

The proposed data-adaptive QSM algorithm may enable quantification of LIC with improved repeatability/reproducibility across different acquisition parameters as 3T.

QSM Throughout the Body

Magnetic materials in tissue, such as iron, calcium, or collagen, can be studied using quantitative susceptibility mapping (QSM). To date, QSM has been overwhelmingly applied in the brain, but is increasingly utilized outside the brain. QSM relies on the effect of tissue magnetic susceptibility sources on the MR signal phase obtained with gradient echo sequence. However, in the body, the chemical shift of fat present within the region of interest contributes to the MR signal phase as well. Therefore, correcting for the chemical shift effect by means of water-fat separation is essential for body QSM. By employing techniques to compensate for cardiac and respiratory motion artifacts, body QSM has been applied to study liver iron and fibrosis, heart chamber blood and placenta oxygenation, myocardial hemorrhage, atherosclerotic plaque, cartilage, bone, prostate, breast calcification, and kidney stone.

Quantitative susceptibility mapping (QSM) of the cardiovascular system: challenges and perspectives

Quantitative susceptibility mapping (QSM) is a powerful, non-invasive, magnetic resonance imaging (MRI) technique that relies on measurement of magnetic susceptibility. So far, QSM has been employed mostly to study neurological disorders characterized by iron accumulation, such as Parkinson's and Alzheimer's diseases. Nonetheless, QSM allows mapping key indicators of cardiac disease such as blood oxygenation and myocardial iron content. For this reason, the application of QSM offers an unprecedented opportunity to gain a better understanding of the pathophysiological changes associated with cardiovascular disease and to monitor their evolution and response to treatment. Recent studies on cardiovascular QSM have shown the feasibility of a non-invasive assessment of blood oxygenation, myocardial iron content and myocardial fibre orientation, as well as carotid plaque composition. Significant technical challenges remain, the most evident of which are related to cardiac and respiratory motion, blood flow, chemical shift effects and susceptibility artefacts. Significant work is ongoing to overcome these challenges and integrate the QSM technique into clinical practice in the cardiovascular field.

Quantitative susceptibility mapping of carotid arterial tissue ex vivo: Assessing sensitivity to vessel microstructural composition

PURPOSE

To characterize microstructural contributions to the magnetic susceptibility of carotid arteries.

METHOD

Arterial vessels were scanned using high-resolution quantitative susceptibility mapping (QSM) at 7 Tesla. Models of vessel degradation were generated using ex vivo porcine carotid arteries that were subjected to several different enzymatic digestion treatments that selectively removed microstructural components (smooth muscle cells, collagen, and elastin). Magnetic susceptibilities measured in these tissue models were compared to those in untreated (native) porcine arteries. Magnetic susceptibility measured in native porcine carotid arteries was further compared to the susceptibility of cadaveric human carotid arteries to investigate their similarity.

RESULTS

The magnetic susceptibility of native porcine vessels was diamagnetic (χ_native = -0.1820 ppm), with higher susceptibilities in all models of vessel degradation (χ_{elastin-degraded} = -0.0163 ppm; χ_{collagen-degraded} = -0.1158 ppm; χ_{decellularized} = -0.1379 ppm; χ_{fixed native} = -0.2199 ppm). Magnetic susceptibility was significantly higher in collagen-degraded compared to native porcine vessels (Tukey-Kramer, P < .01) and between elastin-degraded and all other models (including native, Tukey-Kramer, P < .001). The susceptibility of fixed healthy human arterial tissue was diamagnetic, and no significant difference was found between fixed human and fixed porcine arterial tissue susceptibilities (analysis of variance, P > .05).

CONCLUSIONS

Magnetic susceptibility measured using QSM is sensitive to the microstructural composition of arterial vessels-most notably to collagen. The similarity of human and porcine arterial tissue susceptibility values provides a solid basis for translational studies. Because vessel microstructure becomes disrupted during the onset and progression of carotid atherosclerosis, QSM has the potential to provide a sensitive and specific marker of vessel disease.

Multiecho complex total field inversion method (mcTFI) for improved signal modeling in quantitative susceptibility mapping

PURPOSE

Typical quantitative susceptibility mapping (QSM) reconstruction steps consist of first estimating the magnetization field from the gradient-echo images, and then reconstructing the susceptibility map from the estimated field. The errors from the field-estimation steps may propagate into the final QSM map, and the noise in the estimated field map may no longer be zero-mean Gaussian noise, thus, causing streaking artifacts in the resulting QSM. A multiecho complex total field inversion (mcTFI) method was developed to compute the susceptibility map directly from the multiecho gradient echo images using an improved signal model that retains the Gaussian noise property in the complex domain. It showed improvements in QSM reconstruction over the conventional field-to-source inversion.

METHODS

The proposed mcTFI method was compared with the nonlinear total field inversion (nTFI) method in a numerical brain with hemorrhage and calcification, the numerical brains provided by the QSM Challenge 2.0, 18 brains with intracerebral hemorrhage scanned at 3T, and 6 healthy brains scanned at 7T.

RESULTS

Compared with nTFI, the proposed mcTFI showed more accurate QSM reconstruction around the lesions in the numerical simulations. The mcTFI reconstructed QSM also showed the best image quality with the least artifacts in the brains with intracerebral hemorrhage scanned at 3T and healthy brains scanned at 7T.

CONCLUSION

The proposed multiecho complex total field inversion improved QSM reconstruction over traditional field-to-source inversion through better signal modeling.

Diffusion tensor imaging and arterial tissue: establishing the influence of arterial tissue microstructure on fractional anisotropy, mean diffusivity and tractography

This study investigates diffusion tensor imaging (DTI) for providing microstructural insight into changes in arterial tissue by exploring how cell, collagen and elastin content effect fractional anisotropy (FA), mean diffusivity (MD) and tractography. Five ex vivo porcine carotid artery models (n = 6 each) were compared-native, fixed native, collagen degraded, elastin degraded and decellularised. Vessels were imaged at 7 T using a DTI protocol with b = 0 and 800 s/mm2 and 10 isotopically distributed directions. FA and MD were evaluated in the vessel media and compared across models. FA values measured in native (p < 0.0001), fixed native (p < 0.0001) and collagen degraded (p = 0.0018, p = 0.0016, respectively) were significantly higher than those in elastin degraded and decellularised arteries. Native and fixed native had significantly lower MD values than elastin degraded (p < 0.0001) and decellularised tissue (p = 0.0032, p = 0.0003, respectively). Significantly lower MD was measured in collagen degraded compared with the elastin degraded model (p = 0.0001). Tractography yielded helically arranged tracts for native and collagen degraded vessels only. FA, MD and tractography were found to be highly sensitive to changes in the microstructural composition of arterial tissue, specifically pointing to cell, not collagen, content as the dominant source of the measured anisotropy in the vessel wall.