Generative adversarial networks improve interior computed tomography angiography reconstruction

In interior computed tomography (CT), the x-ray beam is collimated to a limited field-of-view (FOV) (e.g. the volume of the heart) to decrease exposure to adjacent organs, but the resulting image has a severe truncation artifact when reconstructed with traditional filtered back-projection (FBP) type algorithms. In some examinations, such as cardiac or dentomaxillofacial imaging, interior CT could be used to achieve further dose reductions. In this work, we describe a deep learning (DL) method to obtain artifact-free images from interior CT angiography. Our method employs the Pix2Pix generative adversarial network (GAN) in a two-stage process: (1) An extended sinogram is computed from a truncated sinogram with one GAN model, and (2) the FBP reconstruction obtained from that extended sinogram is used as an input to another GAN model that improves the quality of the interior reconstruction. Our double GAN (DGAN) model was trained with 10 000 truncated sinograms simulated from real computed tomography angiography slice images. Truncated sinograms (input) were used with original slice images (target) in training to yield an improved reconstruction (output). DGAN performance was compared with the adaptive de-truncation method, total variation regularization, and two reference DL methods: FBPConvNet, and U-Net-based sinogram extension (ES-UNet). Our DGAN method and ES-UNet yielded the best root-mean-squared error (RMSE) (0.03 +/- 0.01), and structural similarity index (SSIM) (0.92 +/- 0.02) values, and reference DL methods also yielded good results. Furthermore, we performed an extended FOV analysis by increasing the reconstruction area by 10% and 20%. In both cases, the DGAN approach yielded best results at RMSE (0.03 +/- 0.01 and 0.04 +/- 0.01 for the 10% and 20% cases, respectively), peak signal-to-noise ratio (PSNR) (30.5 +/- 2.6 dB and 28.6 +/- 2.6 dB), and SSIM (0.90 +/- 0.02 and 0.87 +/- 0.02). In conclusion, our method was able to not only reconstruct the interior region with improved image quality, but also extend the reconstructed FOV by 20%.Peer reviewe

T₂-weighted magnetic resonance imaging texture as predictor of low back pain:a texture analysis-based classification pipeline to symptomatic and asymptomatic cases

Abstract Low back pain is a very common symptom and the leading cause of disability throughout the world. Several degenerative imaging findings seen on magnetic resonance imaging are associated with low back pain but none of them is specific for the presence of low back pain as abnormal findings are prevalent among asymptomatic subjects as well. The purpose of this population-based study was to investigate if more specific magnetic resonance imaging predictors of low back pain could be found via texture analysis and machine learning. We used this methodology to classify T₂-weighted magnetic resonance images from the Northern Finland Birth Cohort 1966 data to symptomatic and asymptomatic groups. Lumbar spine magnetic resonance imaging was performed using a fast spin-echo sequence at 1.5 T. Texture analysis pipeline consisting of textural feature extraction, principal component analysis, and logistic regression classifier was applied to the data to classify them into symptomatic (clinically relevant pain with frequency ≥30 days and intensity ≥6/10) and asymptomatic (frequency ≤7 days, intensity ≤3/10, and no previous pain episodes in the follow-up period) groups. Best classification results were observed applying texture analysis to the two lowest intervertebral discs (L4-L5 and L5-S1), with accuracy of 83%, specificity of 83%, sensitivity of 82%, negative predictive value of 94%, precision of 56%, and receiver operating characteristic area-under-curve of 0.91. To conclude, textural features from T₂-weighted magnetic resonance images can be applied in low back pain classification

Generative adversarial networks improve interior computed tomography angiography reconstruction

Abstract In interior computed tomography (CT), the x-ray beam is collimated to a limited field-of-view (FOV) (e.g. the volume of the heart) to decrease exposure to adjacent organs, but the resulting image has a severe truncation artifact when reconstructed with traditional filtered back-projection (FBP) type algorithms. In some examinations, such as cardiac or dentomaxillofacial imaging, interior CT could be used to achieve further dose reductions. In this work, we describe a deep learning (DL) method to obtain artifact-free images from interior CT angiography. Our method employs the Pix2Pix generative adversarial network (GAN) in a two-stage process: (1) An extended sinogram is computed from a truncated sinogram with one GAN model, and (2) the FBP reconstruction obtained from that extended sinogram is used as an input to another GAN model that improves the quality of the interior reconstruction. Our double GAN (DGAN) model was trained with 10 000 truncated sinograms simulated from real computed tomography angiography slice images. Truncated sinograms (input) were used with original slice images (target) in training to yield an improved reconstruction (output). DGAN performance was compared with the adaptive de-truncation method, total variation regularization, and two reference DL methods: FBPConvNet, and U-Net-based sinogram extension (ES-UNet). Our DGAN method and ES-UNet yielded the best root-mean-squared error (RMSE) (0.03 ± 0.01), and structural similarity index (SSIM) (0.92 ± 0.02) values, and reference DL methods also yielded good results. Furthermore, we performed an extended FOV analysis by increasing the reconstruction area by 10% and 20%. In both cases, the DGAN approach yielded best results at RMSE (0.03 ± 0.01 and 0.04 ± 0.01 for the 10% and 20% cases, respectively), peak signal-to-noise ratio (PSNR) (30.5 ± 2.6 dB and 28.6 ± 2.6 dB), and SSIM (0.90 ± 0.02 and 0.87 ± 0.02). In conclusion, our method was able to not only reconstruct the interior region with improved image quality, but also extend the reconstructed FOV by 20%

An automatic regularization method:an application for 3-D X-ray micro-CT reconstruction using sparse data

Abstract X-ray tomography is a reliable tool for determining the inner structure of 3-D object with penetrating X-rays. However, traditional reconstruction methods, such as Feldkamp-Davis-Kress (FDK), require dense angular sampling in the data acquisition phase leading to long measurement times, especially in X-ray micro-tomography to obtain high-resolution scans. Acquiring less data using greater angular steps is an obvious way for speeding up the process and avoiding the need to save huge data sets. However, computing 3-D reconstruction from such a sparsely sampled data set is difficult because the measurement data are usually contaminated by errors, and linear measurement models do not contain sufficient information to solve the problem in practice. An automatic regularization method is proposed for robust reconstruction, based on enforcing sparsity in the 3-D shearlet-transform domain. The inputs of the algorithm are the projection data and a priori known expected degree of sparsity, denoted as 0 <; C pr ≤ 1. The number Cpr can be calibrated from a few dense-angle reconstructions and fixed. Human subchondral bone samples were tested, and morphometric parameters of the bone reconstructions were then analyzed using standard metrics. The proposed method is shown to outperform the baseline algorithm (FDK) in the case of sparsely collected data. The number of X-ray projections can be reduced up to 10% of the total amount 300 projections over 180° with uniform angular step while retaining the quality of the reconstruction images and of the morphometric parameters

Machine learning prediction of collagen fiber orientation and proteoglycan content from multiparametric quantitative MRI in articular cartilage

Abstract Background: Machine learning models trained with multiparametric quantitative MRIs (qMRIs) have the potential to provide valuable information about the structural composition of articular cartilage. Purpose: To study the performance and feasibility of machine learning models combined with qMRIs for noninvasive assessment of collagen fiber orientation and proteoglycan content. Study type: Retrospective, animal model. Animal model: An open-source single slice MRI dataset obtained from 20 samples of 10 Shetland ponies (seven with surgically induced cartilage lesions followed by treatment and three healthy controls) yielded to 1600 data points, including 10% for test and 90% for train validation. Field strength/sequence: A 9.4 T MRI scanner/qMRI sequences: T₁, T₂, adiabatic T1ρ and T2ρ, continuous-wave T1ρ and relaxation along a fictitious field (TRAFF) maps. Assessment: Five machine learning regression models were developed: random forest (RF), support vector regression (SVR), gradient boosting (GB), multilayer perceptron (MLP), and Gaussian process regression (GPR). A nested cross-validation was used for performance evaluation. For reference, proteoglycan content and collagen fiber orientation were determined by quantitative histology from digital densitometry (DD) and polarized light microscopy (PLM), respectively. Statistical tests: Normality was tested using Shapiro–Wilk test, and association between predicted and measured values was evaluated using Spearman’s Rho test. A P-value of 0.05 was considered as the limit of statistical significance. Results: Four out of the five models (RF, GB, MLP, and GPR) yielded high accuracy (R² = 0.68–0.75 for PLM and 0.62–0.66 for DD), and strong significant correlations between the reference measurements and predicted cartilage matrix properties (Spearman’s Rho = 0.72–0.88 for PLM and 0.61–0.83 for DD). GPR algorithm had the highest accuracy (R² = 0.75 and 0.66) and lowest prediction-error (root mean squared [RMSE] = 1.34 and 2.55) for PLM and DD, respectively. Data conclusion: Multiparametric qMRIs in combination with regression models can determine cartilage compositional and structural features, with higher accuracy for collagen fiber orientation than proteoglycan content. Evidence level: 2 Technical efficacy: Stage

Framework for photon counting quantitative material decomposition

Abstract In this work, the accuracy of material decomposition (MD) using an energy discriminating photon counting detector was studied. An MD framework was established and validated using calcium hydroxyapatite (CaHA) inserts of known densities (50 mg/cm³, 100 mg/cm³, 250 mg/cm³, 400 mg/cm³), and diameters (1.2 mm, 3.0 mm, 5.0 mm). These inserts were placed in a cardiac rod phantom that mimics a tissue equivalent heart and measured using an experimental photon counting detector cone beam computed tomography (PCD-CBCT) setup. The quantitative coronary calcium scores (density, mass, and volume) obtained from the MD framework were compared with the nominal values. Additionally, three different calibration techniques, signal-toequivalent thickness calibration (STC), polynomial correction (PC), and projected equivalent thickness calibration (PETC) were compared to investigate the effect of the calibration method on the quantitative values. The obtained MD estimates agreed well with the nominal values for density (mass) with mean absolute percent errors (MAPEs) 8 ± 11% (9 ± 15%) and 4 ± 6% (9 ± 14%) for STC and PETC calibration methods, respectively. PC displayed large MAPEs for density (27 ± 9%), and mass (25 ± 12%). Volume estimation resulted in large deviations between true and measured values with notable MAPEs for STC (40 ± 90%), PC (40 ± 80%), and PETC (40 ± 90%). The framework demonstrated the feasibility of quantitative CaHA mass and density scoring using PCD-CBCT

Assessing post-traumatic changes in cartilage using T1ρ dispersion parameters

Abstract Degeneration of cartilage can be studied non-invasively with quantitative MRI. A promising parameter for detecting early osteoarthritis in articular cartilage is T1ρ, which can be tuned via the amplitude of the spin-lock pulse. By measuring T1ρ at several spin-lock amplitudes, the dispersion of T1ρis obtained. The aim of this study is to find out if the dispersion contains diagnostically relevant information complementary to a T1ρ measurement at a single spin-lock amplitude. To this end, five differently acquired dispersion parameters are utilized; A, B, τc, T1ρ/T2, and R2 − R1ρ. An open dataset of an equine model of post-traumatic cartilage was utilized to assess the T1ρ dispersion parameters for the evaluation of cartilage degeneration. Firstly, the parameters were compared for their sensitivity in detecting degenerative changes. Secondly, the relationship of the dispersion parameters to histological and biomechanical reference parameters was studied. Parameters A, T1ρ/T2, and R2 − R1ρ were found to be sensitive to lesion-induced changes in the cartilage within sample. Strong correlations of several dispersion parameters with optical density, as well as with collagen fibril angle were found. Most of the dispersion parameters correlated strongly with individual T1ρ values. The results suggest that dispersion parameters can in some cases provide a more accurate description of the biochemical composition of cartilage as compared to conventional MRI parameters. However, in most cases the information given by the dispersion parameters is more of a refinement than complementary to conventional quantitative MRI

Iterative and discrete reconstruction in the evaluation of the rabbit model of osteoarthritis

Abstract Micro-computed tomography (µCT) is a standard method for bone morphometric evaluation. However, the scan time can be long and the radiation dose during the scan may have adverse effects on test subjects, therefore both of them should be minimized. This could be achieved by applying iterative reconstruction (IR) on sparse projection data, as IR is capable of producing reconstructions of sufficient image quality with less projection data than the traditional algorithm requires. In this work, the performance of three IR algorithms was assessed for quantitative bone imaging from low-resolution data in the evaluation of the rabbit model of osteoarthritis. Subchondral bone images were reconstructed with a conjugate gradient least squares algorithm, a total variation regularization scheme, and a discrete algebraic reconstruction technique to obtain quantitative bone morphometry, and the results obtained in this manner were compared with those obtained from the reference reconstruction. Our approaches were sufficient to identify changes in bone structure in early osteoarthritis, and these changes were preserved even when minimal data were provided for the reconstruction. Thus, our results suggest that IR algorithms give reliable performance with sparse projection data, thereby recommending them for use in µCT studies where time and radiation exposure are preferably minimized

T2* and quantitative susceptibility mapping in an equine model of post-traumatic osteoarthritis : assessment of mechanical and structural properties of articular cartilage

Objective: To investigate the potential of quantitative susceptibility mapping (QSM) and T2* relaxation time mapping to determine mechanical and structural properties of articular cartilage via univariate and multivariate analysis. Methods: Samples were obtained from a cartilage repair study, in which surgically induced full-thickness chondral defects in the stifle joints of seven Shetland ponies caused post-traumatic osteoarthritis (14 samples). Control samples were collected from non-operated joints of three animals (6 samples). Magnetic resonance imaging (MRI) was performed at 9.4 T, using a 3-D multi-echo gradient echo sequence. Biomechanical testing, digital densitometry (DD) and polarized light microscopy (PLM) were utilized as reference methods. To compare MRI parameters with reference parameters (equilibrium and dynamic moduli, proteoglycan content, collagen fiber angle and -anisotropy), depth-wise profiles of MRI parameters were acquired at the biomechanical testing locations. Partial least squares regression (PLSR) and Spearman's rank correlation were utilized in data analysis. Results: PLSR indicated a moderate-to-strong correlation (ρ = 0.49–0.66) and a moderate correlation (ρ = 0.41–0.55) between the reference values and T2* relaxation time and QSM profiles, respectively (excluding superficial-only results). PLSR correlations were noticeably higher than direct correlations between bulk MRI and reference parameters. 3-D parametric surface maps revealed spatial variations in the MRI parameters between experimental and control groups. Conclusion: Quantitative parameters from 3-D multi-echo gradient echo MRI can be utilized to predict the properties of articular cartilage. With PLSR, especially the T2* relaxation time profile appeared to correlate with the properties of cartilage. Furthermore, the results suggest that degeneration affects the QSM-contrast in the cartilage. However, this change in contrast is not easy to quantify

Evaluation of articular cartilage with quantitative MRI in an equine model of post-traumatic osteoarthritis

Chondral lesions lead to degenerative changes in the surrounding cartilage tissue, increasing the risk of developing post-traumatic osteoarthritis (PTOA). This study aimed to investigate the feasibility of quantitative magnetic resonance imaging (qMRI) for evaluation of articular cartilage in PTOA. Articular explants containing surgically induced and repaired chondral lesions were obtained from the stifle joints of seven Shetland ponies (14 samples). Three age-matched nonoperated ponies served as controls (six samples). The samples were imaged at 9.4 T. The measured qMRI parameters included T, T, continuous-wave T (CWT), adiabatic T (AdT), and T (AdT) and relaxation along a fictitious field (T). For reference, cartilage equilibrium and dynamic moduli, proteoglycan content and collagen fiber orientation were determined. Mean values and profiles from full-thickness cartilage regions of interest, at increasing distances from the lesions, were used to compare experimental against control and to correlate qMRI with the references. Significant alterations were detected by qMRI parameters, including prolonged T, CWT, and AdT in the regions adjacent to the lesions. The changes were confirmed by the reference methods. CWT was more strongly associated with the reference measurements and prolonged in the affected regions at lower spin-locking amplitudes. Moderate to strong correlations were found between all qMRI parameters and the reference parameters (ρ = −0.531 to −0.757). T, low spin-lock amplitude CWT, and AdT were most responsive to changes in visually intact cartilage adjacent to the lesions. In the context of PTOA, these findings highlight the potential of T, CWT, and AdT in evaluation of compositional and structural changes in cartilage