Search

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 2

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.

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%.

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.

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