Phantom Scans Research Articles

High-quality multimodal MRI with simultaneous EEG using conductive ink and polymer-thick film nets

Objective. Combining magnetic resonance imaging (MRI) and electroencephalography (EEG) provides a powerful tool for investigating brain function at varying spatial and temporal scales. Simultaneous acquisition of both modalities can provide unique information that a single modality alone cannot reveal. However, current simultaneous EEG-fMRI studies are limited to a small set of MRI sequences due to the image quality and safety limitations of commercially available MR-conditional EEG nets. We tested whether the Inknet2, a high-resistance polymer thick film based EEG net that uses conductive ink, could enable the acquisition of a variety of MR image modalities with minimal artifacts by reducing the radiofrequency-shielding caused by traditional MR-conditional nets.Approach. We first performed simulations to model the effect of the EEG nets on the magnetic field and image quality. We then performed phantom scans to test image quality with a conventional copper EEG net, with the new Inknet2, and without any EEG net. Finally, we scanned five human subjects at 3 Tesla (3 T) and three human subjects at 7 Tesla (7 T) with and without the Inknet2 to assess structural and functional MRI image quality.Main results. Across these simulations, phantom scans, and human studies, the Inknet2 induced fewer artifacts than the conventional net and produced image quality similar to scans with no net present.Significance. Our results demonstrate that high-quality structural and functional multimodal imaging across a variety of MRI pulse sequences at both 3 T and 7 T is achievable with an EEG net made with conductive ink and polymer thick film technology.

Optimal weighting strategies for maximizing contrast-to-noise ratio in photon counting CT images.

Photon counting detectors (PCDs) with energy discrimination capabilities have the potential to generate grayscale CT images with improved contrast-to-noise ratio (CNR) through optimal weighting of their spectral measurements. This study evaluates the CNR performance of grayscale CT projections and images generated from spectral measurements of PCDs using three energy-weighting strategies: pre-log weighting, post-log weighting, and material decomposition (MD) weighting. This study provides the expressions of optimal weights and maximum achievable CNR of these energy-weighting strategies, which only require the knowledge of detected bin counts and do not require information of PCD energy responses or imaging techniques. We defined and solved a generalized eigenvalue problem to obtain the maximum achievable CNR in the projection domain for low-contrast tasks using three energy-weighting strategies: pre-log weighting (weighted sum of energy bin counts), post-log weighting (weighted sum of line integrals), and MD weighting (weighted sum of basis material thicknesses, which is equivalent to virtual monoenergetic images [VMIs]). These expressions only contain energy bin counts from PCD measurements. We used a realistic PCD energy response model to simulate the detected bin counts and conducted Monte Carlo simulations of different contrast tasks and phantoms to evaluate the projection- and image-domain CNR performance of these energy-weighting strategies. Additionally, the total counts method (a special case of pre-log weighting with unity weights) was included for comparison. We also conducted Gammex head and body phantom scans on an edge-on-irradiated silicon PCCT prototype to evaluate the image-domain CNR performance of these energy-weighting strategies. The results show that pre-log, post-log, and MD weighting strategies generate approximately equal projection-domain maximum achievable CNR, with a difference of less than 2%, and outperform the total counts method. These three energy-weighting strategies also generate approximately equal image-domain maximum CNR when the contrast task is located at the center of a homogeneous phantom. Pre-log weighting generates the highest image-domain CNR for an off-center contrast task location or inhomogeneous phantoms while also outperforming the total counts method. We derived the expression of projection-domain maximum achievable CNR using three energy-weighting strategies. Our results suggest that using pre-log weighting strategies enables fast grayscale CT image generation with high CNR from spectral PCD measurements for inhomogeneous phantoms and off-center region of interests (ROIs).

Rapid and quantitative CEST-MRI sequence using water presaturation.

Despite the significant potential for in vivo metabolic imaging in preclinical and clinical applications, CEST MRI suffers from long scan time and inaccurate quantification. This study aims to suppress the contaminations among signals under different frequencies, which could shorten the TR and thereby facilitate CEST imaging acceleration and quantification. A novel sequence is proposed by applying a water-presaturation (WPS) module at the beginning of each TR. WPS CEST quickly knocks down the residual signal from previous TRs so that the magnetization of all TRs recovers from zero, which aligns well with the formula of quasi-steady-state theorem and enables accurate quantification within shorter TR. WPS CEST was assessed by simulations, creatine phantom, and healthy human brain scans at 3 T. In simulation and phantom experiment, WPS CEST allows accurate estimation of exchange rate (ksw) using omega plot and using shorter delay time (Td) and saturation time (Ts) (e.g., 1 s/1 s) compared with the conventional CEST. Simulations further showed that WPS CEST could obtain consistent spin-lock relaxation (R1ρ) values over varied Tds and Tss. Six human scans indicated that R1ρ collected from conventional sequences showed significant differences between two groups with Td and Ts of (1 s/1 s) and (2 s/2 s) (amide: 1.721 ± 0.051 s-1 vs. 1.622 ± 0.050 s-1, p = 0.001; nuclear Overhauser enhancement: 1.792 ± 0.046 s-1 vs. 1.687 ± 0.053 s-1, p = 0.004), whereas WPS CEST scans using these 2 Td/Ts values obtained the same mean R1ρ (amide: 1.616 ± 0.053 s-1 vs. 1.616 ± 0.048 s-1, p = 0.862; nuclear Overhauser enhancement: 1.688 ± 0.064 s-1 vs. 1.684 ± 0.054 s-1, p = 0.544). WPS CEST demonstrated accurate quantitation within shorter TR compared with conventional sequences, and thereby may allow rapid quantitative CEST scans in various situations.

Particle‐based MR modeling with diffusion, microstructure, and enzymatic reactions

AbstractPurposeTo develop a novel particle‐based in silico MR model and demonstrate applications of this model to signal mechanisms which are affected by the spatial organization of particles, including metabolic reaction kinetics, microstructural effects on diffusion, and radiofrequency (RF) refocusing effects in gradient‐echo sequences.MethodsThe model was developed by integrating a forward solution of the Bloch equations with a Brownian dynamics simulator. Simulation configurations were then designed to model MR signal dynamics of interest, with a primary focus on hyperpolarized 13C MRI methods. Phantom scans and spectrophotometric assays were conducted to validate model results in vitro.ResultsThe model accurately reproduced the reaction kinetics of enzyme‐mediated conversion of pyruvate to lactate. When varying proportions of restrictive structure were added to the reaction volume, nonlinear changes in the reaction rate measured in vitro were replicated in silico. Modeling of RF refocusing effects characterized the degree of diffusion‐weighted contribution from preserved residual magnetization in nonspoiled gradient‐echo sequences.ConclusionsThese results show accurate reproduction of a range of MR signal mechanisms, establishing the model's capability to investigate the multifactorial signal dynamics such as those underlying hyperpolarized 13C MRI data.

Run-time motion and first-order shim control by expanded servo navigation.

To provide a navigator-based run-time motion and first-order field correction for three-dimensional human brain imaging with high precision, minimal calibration and acquisition, and fast processing. A complex-valued linear perturbation model with feedback control is extended to estimate and correct for gradient shim fields using orbital navigators (2.3 ms). Two approaches for sensitizing the model to gradient fields are presented, one based on finite differences with three additional navigators, and another projection-based approximation requiring no additional navigators. A mechanism for noise decorrelation of the matrix and the data is proposed and evaluated to reduce unwanted parameter biases. The rigid motion and first-order field control achieves robust motion and gradient shim corrections improving image quality in a series of phantom and in vivo experiments with varying field conditions. In phantom scans, magnet drifts, forced gradient field perturbations and field distortions from shifts of a second bottle phantom are successfully corrected. Field estimates of the magnet drifts are in good agreement with concurrent field probe measurements. For in vivo scans, the proposed method mitigates field variations from torso motions while being robust to head motion. In vivo gradient field precisions were along with single-digit micrometer and millidegree rigid precisions. The navigator-based method achieves accurate, high-precision run-time motion and field corrections with low sequence impact and calibration requirements.

Exploring synthetic datasets for computer-aided detection: a case study using phantom scan data for enhanced lung nodule false positive reduction.

Synthetic datasets hold the potential to offer cost-effective alternatives to clinical data, ensuring privacy protections and potentially addressing biases in clinical data. We present a method leveraging such datasets to train a machine learning algorithm applied as part of a computer-aided detection (CADe) system. Our proposed approach utilizes clinically acquired computed tomography (CT) scans of a physical anthropomorphic phantom into which manufactured lesions were inserted to train a machine learning algorithm. We treated the training database obtained from the anthropomorphic phantom as a simplified representation of clinical data and increased the variability in this dataset using a set of randomized and parameterized augmentations. Furthermore, to mitigate the inherent differences between phantom and clinical datasets, we investigated adding unlabeled clinical data into the training pipeline. We apply our proposed method to the false positive reduction stage of a lung nodule CADe system in CT scans, in which regions of interest containing potential lesions are classified as nodule or non-nodule regions. Experimental results demonstrate the effectiveness of the proposed method; the system trained on labeled data from physical phantom scans and unlabeled clinical data achieves a sensitivity of 90% at eight false positives per scan. Furthermore, the experimental results demonstrate the benefit of the physical phantom in which the performance in terms of competitive performance metric increased by 6% when a training set consisting of 50 clinical CT scans was enlarged by the scans obtained from the physical phantom. The scalability of synthetic datasets can lead to improved CADe performance, particularly in scenarios in which the size of the labeled clinical data is limited or subject to inherent bias. Our proposed approach demonstrates an effective utilization of synthetic datasets for training machine learning algorithms.

Asynchronous calibration of a CT scanner for bone mineral density estimation: sources of error and correction.

The estimation of BMD with CT scans requires a calibration method, usually based on a phantom. In asynchronous calibration, the phantom is scanned separately from the patient. A standardized acquisition protocol must be used to avoid variations between patient and phantom. However, variations can still be induced, for example, by temporal fluctuations or patient characteristics. Based on the further use of 739 forensic and 111 clinical CT scans, this study uses the proximal femur BMD value ("total hip") to assess asynchronous calibration accuracy, using in-scan calibration as ground truth. It identifies the parameters affecting the calibration accuracy and quantifies their impact. For time interval and table height, the impact was measured by calibrating the CT scan twice (once using the phantom scan with closest acquisition parameters and once using a phantom scan with standard values) and comparing the calibration accuracy. For other parameters such as body weight, the impact was measured by computing a linear regression between parameter values and calibration accuracy. Finally, this study proposes correction methods to reduce the effect of these parameters and improve the calibration accuracy. The BMD error of the asynchronous calibration, using the phantom scan with the closest acquisition parameters, was -1.2 ± 1.7% for the forensic and - 1.6 ± 3.5% for the clinical dataset. Among the parameters studied, time interval and body weight were identified as the main sources of error for asynchronous calibration, followed by table height and reconstruction kernel. Based on these results, a correction method was proposed to improve the calibration accuracy.

Ultra-low dose hip CT-based automated measurement of volumetric bone mineral density at proximal femoral subregions.

Forty to fifty percent of women and 13%-22% of men experience an osteoporosis-related fragility fracture in their lifetimes. After the age of 50 years, the risk of hip fracture doubles in every 10 years. x-Ray based DXA is currently clinically used to diagnose osteoporosis and predict fracture risk. However, it provides only 2-D representation of bone and is associated with other technical limitations. Thus, alternative methods are needed. To develop and evaluate an ultra-low dose (ULD) hip CT-based automated method for assessment of volumetric bone mineral density (vBMD) at proximal femoral subregions. An automated method was developed to segment the proximal femur in ULD hip CT images and delineate femoral subregions. The computational pipeline consists of deep learning (DL)-based computation of femur likelihood map followed by shape model-based femur segmentation and finite element analysis-based warping of a reference subregion labeling onto individual femur shapes. Finally, vBMD is computed over each subregion in the target image using a calibration phantom scan. A total of 100 participants (50 females) were recruited from the Genetic Epidemiology of COPD (COPDGene) study, and ULD hip CT imaging, equivalent to 18 days of background radiation received by U.S. residents, was performed on each participant. Additional hip CT imaging using a clinical protocol was performed on 12 participants and repeat ULD hip CT was acquired on another five participants. ULD CT images from 80 participants were used to train the DL network; ULD CT images of the remaining 20 participants as well as clinical and repeat ULD CT images were used to evaluate the accuracy, generalizability, and reproducibility of segmentation of femoral subregions. Finally, clinical CT and repeat ULD CT images were used to evaluate accuracy and reproducibility of ULD CT-based automated measurements of femoral vBMD. Dice scores of accuracy (n=20), reproducibility (n=5), and generalizability (n=12) of ULD CT-based automated subregion segmentation were 0.990, 0.982, and 0.977, respectively, for the femoral head and 0.941, 0.970, and 0.960, respectively, for the femoral neck. ULD CT-based regional vBMD showed Pearson and concordance correlation coefficients of 0.994 and 0.977, respectively, and a root-mean-square coefficient of variation (RMSCV) (%) of 1.39% with the clinical CT-derived reference measure. After 3-digit approximation, each of Pearson and concordance correlation coefficients as well as intraclass correlation coefficient (ICC) between baseline and repeat scans were 0.996 with RMSCV of 0.72%. Results of ULD CT-based bone analysis on 100 participants (age (mean±SD) 73.6±6.6 years) show that males have significantly greater (p<0.01) vBMD at the femoral head and trochanteric regions than females, while females have moderately greater vBMD (p=0.05) at the medial half of the femoral neck than males. Deep learning, combined with shape model and finite element analysis, offers an accurate, reproducible, and generalizable algorithm for automated segmentation of the proximal femur and anatomic femoral subregions using ULD hip CT images. ULD CT-based regional measures of femoral vBMD are accurate and reproducible and demonstrate regional differences between males and females.

A robust index for metal artifact quantification in computed tomography.

Objective assessment of metal artifact strength and the effectiveness of metal artifact reduction algorithms in computed tomography requires a quantitative metric. Metrics described in the literature are typically employed to compare the artifact strength in images reconstructed from the same raw data, but their robustness against varying scan conditions and repeated scans over time as it occurs in periodic quality assurance has not been investigated. A new robust metric for quantifying metal artifacts in computed-tomography images is proposed and compared to other commonly used metrics. The proposed artifact metric is based on the location parameter of the Gumbel distribution, described previously in the literature, but normalized to the location parameter in a background region-of-interest to obtain a noise-independent artifact metric. The metric was compared to three other quantitative metal artifact metrics (artifact-index, contrast-to-noise ratio, Gumbel-evaluation method) by evaluating metals artifacts in phantom scans and in clinical images. Robustness of the artifact metrics was evaluated using repeated scans with varying noise and against small variations in the selected regions-of-interest. The proposed artifact metric was independent of the underlying image noise and could be reproduced more consistently under slight changes of the region-of-interest within the artifact than the other investigated methods. The coefficient-of-variation was 5.7% on average with varying regions-of-interest in phantom scans and 2.5% in patient scans compared to 9.2% in phantoms scans and 9.9% in patient scans for the next-best performing noise-independent metric. Setup reproducibility was better than 5% and was comparable to the other metrics. The new metric correlated linearly with the artifact strength. The contrast-to-noise ratio, although often used in artifact quantification, was found to be an inadequate metric due to its lack of robustness against minute changes in the position, size, and pixel values of the region-of-interest chosen for calculating the metric and because it showed no correlation with the artifact strength. A new metal artifact metric has been proposed that is robust under changing scan conditions and less sensitive to user-dependent choices of the region-of-interest than other metrics. The new metric is straightforward to calculate and simple to implement in software commonly used for evaluation of medical imaging systems.

Virtual noncontrast images reveal gouty tophi in contrast-enhanced dual-energy CT: a phantom study

BackgroundDual-energy computed tomography (DECT) is useful for detecting gouty tophi. While iodinated contrast media (ICM) might enhance the detection of monosodium urate crystals (MSU), higher iodine concentrations hamper their detection. Calculating virtual noncontrast (VNC) images might improve the detection of enhancing tophi. The aim of this study was to evaluate MSU detection with VNC images from DECT acquisitions in phantoms, compared against the results with standard DECT reconstructions.MethodsA grid-like and a biophantom with 25 suspensions containing different concentrations of ICM (0 to 2%) and MSU (0 to 50%) were scanned with sequential single-source DECT using an ascending order of tube current time product at 80 kVp (16.5–220 mAs) and 135 kVp (2.75–19.25 mAs). VNC images were equivalently reconstructed at 80 and 135 kVp. Two-material decomposition analysis for MSU detection was applied for the VNC and conventional CT images. MSU detection and attenuation values were compared in both modalities.ResultsFor 0, 0.25, 0.5, 1, and 2% ICM, the average detection indices (DIs) for all MSU concentrations (35–50%) with VNC postprocessing were respectively 25.2, 36.6, 30.9, 38.9, and 45.8% for the grid phantom scans and 11.7, 9.4, 5.5, 24.0, and 25.0% for the porcine phantom scans. In the conventional CT image group, the average DIs were respectively 35.4, 54.3, 45.4, 1.0, and 0.0% for the grid phantom and 19.4, 17.9, 3.0, 0.0, and 0.0% for the porcine phantom scans.ConclusionsVNC effectively reduces the suppression of information caused by high concentrations of ICM, thereby improving the detection of MSU.Relevance statementContrast-enhanced DECT alone may suffice for diagnosing gout without a native acquisition.Key points• Highly concentrated contrast media hinders monosodium urate crystal detection in CT imaging• Virtual noncontrast imaging redetects monosodium urate crystals in high-iodinated contrast media concentrations.• Contrast-enhanced DECT alone may suffice for diagnosing gout without a native acquisition.Graphical

Efficient 3D cone trajectory design for improved combined angiographic and perfusion imaging using arterial spin labeling.

To improve the spatial resolution and repeatability of a non-contrast MRI technique for simultaneous time resolved 3D angiography and perfusion imaging by developing an efficient 3D cone trajectory design. A novel parameterized 3D cone trajectory design incorporating the 3D golden angle was integrated into 4D combined angiography and perfusion using radial imaging and arterial spin labeling (CAPRIA) to achieve higher spatial resolution and sampling efficiency for both dynamic angiography and perfusion imaging with flexible spatiotemporal resolution. Numerical simulations and physical phantom scanning were used to optimize the cone design. Eight healthy volunteers were scanned to compare the original radial trajectory in 4D CAPRIA with our newly designed cone trajectory. A locally low rank reconstruction method was used to leverage the complementary k-space sampling across time. The improved sampling in the periphery of k-space obtained with the optimized 3D cone trajectory resulted in improved spatial resolution compared with the radial trajectory in phantom scans. Improved vessel sharpness and perfusion visualization were also achieved in vivo. Less dephasing was observed in the angiograms because of the short TE of our cone trajectory and the improved k-space sampling efficiency also resulted in higher repeatability compared to the original radial approach. The proposed 3D cone trajectory combined with 3D golden angle ordering resulted in improved spatial resolution and image quality for both angiography and perfusion imaging and could potentially benefit other applications that require an efficient sampling scheme with flexible spatial and temporal resolution.

A modular motion compensation pipeline for prospective respiratory motion correction of multi-nuclear MR spectroscopy

Magnetic resonance (MR) acquisitions of the torso are frequently affected by respiratory motion with detrimental effects on signal quality. The motion of organs inside the body is typically decoupled from surface motion and is best captured using rapid MR imaging (MRI). We propose a pipeline for prospective motion correction of the target organ using MR image navigators providing absolute motion estimates in millimeters. Our method is designed to feature multi-nuclear interleaving for non-proton MR acquisitions and to tolerate local transmit coils with inhomogeneous field and sensitivity distributions. OpenCV object tracking was introduced for rapid estimation of in-plane displacements in 2D MR images. A full three-dimensional translation vector was derived by combining displacements from slices of multiple and arbitrary orientations. The pipeline was implemented on 3 T and 7 T MR scanners and tested in phantoms and volunteers. Fast motion handling was achieved with low-resolution 2D MR image navigators and direct implementation of OpenCV into the MR scanner’s reconstruction pipeline. Motion-phantom measurements demonstrate high tracking precision and accuracy with minor processing latency. The feasibility of the pipeline for reliable in-vivo motion extraction was shown on heart and kidney data. Organ motion was manually assessed by independent operators to quantify tracking performance. Object tracking performed convincingly on 7774 navigator images from phantom scans and different organs in volunteers. In particular the kernelized correlation filter (KCF) achieved similar accuracy (74%) as scored from inter-operator comparison (82%) while processing at a rate of over 100 frames per second. We conclude that fast 2D MR navigator images and computer vision object tracking can be used for accurate and rapid prospective motion correction. This and the modular structure of the pipeline allows for the proposed method to be used in imaging of moving organs and in challenging applications like cardiac magnetic resonance spectroscopy (MRS) or magnetic resonance imaging (MRI) guided radiotherapy.

Freehand 3D Ultrasound Imaging Based on Probe-mounted Vision and IMU System

ObjectivesFreehand three-dimensional (3D) ultrasound (US) is of great significance for clinical diagnosis and treatment, it is often achieved with the aid of external devices (optical and/or electromagnetic, etc.) that monitor the location and orientation of the US probe. However, this external monitoring is often impacted by imaging environment such as optical occlusions and/or electromagnetic (EM) interference. MethodsTo address the above issues, we integrated a binocular camera and an inertial measurement unit (IMU) on a US probe. Subsequently, we built a tight coupling model utilizing the unscented Kalman algorithm based on Lie groups (UKF-LG), combining vision and inertial information to infer the probe's movement, through which the position and orientation of the US image frame are calculated. Finally, the volume data was reconstructed with the voxel-based hole-filling method. ResultsThe experiments including calibration experiments, tracking performance evaluation, phantom scans, and real scenarios scans have been conducted. The results show that the proposed system achieved the accumulated frame position error of 3.78 mm and the orientation error of 0.36° and reconstructed 3D US images with high quality in both phantom and real scenarios. ConclusionsThe proposed method has been demonstrated to enhance the robustness and effectiveness of freehand 3D US. Follow-up research will focus on improving the accuracy and stability of multi-sensor fusion to make the system more practical in clinical environments.

Retrospective correction of second-order concomitant fields in 3D axial stack-of-spirals imaging on a high-performance gradient system.

MRI using 3D stack-of-spirals (SoS) readout on a high-performance gradient system is subject to strong second-order, spatially varying concomitant fields, which can lead to signal dropout and blurring artifacts that become more significant at locations farther from the gradient isocenter. A method for compensating for second-order concomitant fields in 3D axial SoS image reconstruction is described. We retrospectively correct for second-order concomitant field-induced phase error in the 3D SoS data by slice-dependent k-space phase compensation based on the nominal spiral readout trajectories. The effectiveness of the method was demonstrated in phantom and healthy volunteer scans in which 3D pseudo-continuous arterial spin labeling imaging was performed with SoS fast spin-echo readout at 3 T. Substantial reduction in blurring was observed with the proposed method. In phantom scans, blurring was reduced by about 53% at 98 mm from the gradient isocenter. In the in vivo 3D pseudo-continuous arterial spin labeling scans, differences of up to 10% were observed at 78 mm from the isocenter, especially around the white-matter and gray-matter interfaces, between the corrected and uncorrected proton density images, perfusion-weighted images, and cerebral blood flow maps. The described retrospective correction method provides a means to correct erroneous phase accruals due to second-order concomitant fields in 3D axial stack-of-spirals imaging.

An automated technique for global noise level measurement in CT image with a conjunction of image gradient.

Automated assessment of noise level in clinical computed tomography (CT) images is a crucial technique for evaluating and ensuring the quality of these images. There are various factors that can impact CT image noise, such as statistical noise, electronic noise, structure noise, texture noise, artifact noise, etc. In this study, a method was developed to measure the global noise index (GNI) in clinical CT scans due to the fluctuation of x-ray quanta. Initially, a noise map is generated by sliding a 10 × 10 pixel for calculating Hounsfield unit (HU) standard deviation and the noise map is further combined with the gradient magnitude map. By employing Boolean operation, pixels with high gradients are excluded from the noise histogram generated with the noise map. By comparing the shape of the noise histogram from this method with Christianson's tissue-type global noise measurement algorithm, it was observed that the noise histogram computed in anthropomorphic phantoms had a similar shape with a close GNI value. In patient CT images, excluding the HU deviation due the structure change demonstrated to have consistent GNI values across the entire CT scan range with high heterogeneous tissue compared to the GNI values using Christianson's tissue-type method. The proposed GNI was evaluated in phantom scans and was found to be capable of comparing scan protocols between different scanners. The variation of GNI when using different reconstruction kernels in clinical CT images demonstrated a similar relationship between noise level and kernel sharpness as observed in uniform phantom: sharper kernel resulted in noisier images. This indicated that GNI was a suitable index for estimating the noise level in clinical CT images with either a smooth or grainy appearance. The study's results suggested that the algorithm can be effectively utilized to screen the noise level for a better CT image quality control.

Hepatic steatosis modeling and MRI signal simulations for comparison of single- and dual-R2* models and estimation of fat fraction at 1.5T and 3T

Background and objectiveMagnetic resonance imaging (MRI) has emerged as a noninvasive clinical tool for assessment of hepatic steatosis. Multi-spectral fat-water MRI models, incorporating single or dual transverse relaxation decay rate(s) (R2*) have been proposed for accurate fat fraction (FF) estimation. However, it is still unclear whether single- or dual-R2* model accurately mimics in vivo signal decay for precise FF estimation and the impact of signal-to-noise ratio (SNR) on each model performance. Hence, this study aims to construct virtual steatosis models and synthesize MRI signals with different SNRs to systematically evaluate the accuracy of single- and dual-R2* models for FF and R2* estimations at 1.5T and 3.0T. MethodsRealistic hepatic steatosis models encompassing clinical FF range (0–60 %) were created using morphological features of fat droplets (FDs) extracted from human liver biopsy samples. MRI signals were synthesized using Monte Carlo simulations for noise-free (SNRideal) and varying SNR conditions (5–100). Fat-water phantoms were scanned with different SNRs to validate simulation results. Fat water toolbox was used to calculate R2* and FF for both single- and dual-R2* models. The model accuracies in R2* and FF estimates were analyzed using linear regression, bias plot and heatmap analysis. ResultsThe virtual steatosis model closely mimicked in vivo fat morphology and Monte Carlo simulation produced realistic MRI signals. For SNRideal and moderate-high SNRs, water R2* (R2*W) by dual-R2* and common R2* (R2*com) by single-R2* model showed an excellent agreement with slope close to unity (0.95–1.01) and R2 > 0.98 at both 1.5T and 3.0T. In simulations, the R2*com-FF and R2*W-FF relationships exhibited slopes similar to in vivo calibrations, confirming the accuracy of our virtual models. For SNRideal, fat R2* (R2*F) was similar to R2*W and dual-R2* model showed slightly higher accuracy in FF estimation. However, in the presence of noise, dual-R2* produced higher FF bias with decreasing SNR, while leading to only marginal improvement for high SNRs and in regions dominated by fat and water. In contrast, single-R2* model was robust and produced accurate FF estimations in simulations and phantom scans with clinical SNRs. ConclusionOur study demonstrates the feasibility of creating virtual steatosis models and generating MRI signals that mimic in vivo morphology and signal behavior. The single-R2* model consistently produced lower FF bias for clinical SNRs across entire FF range compared to dual-R2* model, hence signifying that single-R2* model is optimal for assessing hepatic steatosis.

Improving the quality of computed tomography brain images in the presence of cochlear implant induced metal artefacts through the additional use of tissue mimicking materials alongside metal artefact reduction software

IntroductionMetal artefact reduction software (MAR) can be used to improve Computed Tomography (CT) image quality in the presence of implanted metalwork; however, this software is not effective for superficial metallic structures such as cochlear implants (CI).This study aimed to investigate whether the effectiveness of MAR software could be improved for brain scans with CI present through the use of tissue mimicking materials (TMM) placed exteriorly to the implant. MethodsIn this two-part study, a CI was positioned on the surface of water and anthropomorphic phantoms and imaged using a helical CT brain protocol. Three TMM, Superflab, Sure Thermal heat packs, and Bart's Bolus, were utilised and images were acquired to assess the resulting artefact reduction in terms of CT numbers, noise and artefact index (Aind). Changes in CTDIvol were assessed for the anthropomorphic phantom scans. ResultsIn the water phantom, statistically significant reductions in CT number (p = 0.038) and noise (p = 0.033) were observed for Superflab, whilst the heat packs produced similar significant reductions in CT number (p < 0.001) and noise (p = 0.001) for the anthropomorphic phantom images. Aind values were significantly reduced through the use of Superflab (p = 0.009) and the heat packs (p < 0.001). No significant effects were observed for Bart's Bolus. CTDIvol increases of generally less than 5% were observed for scans with TMM in place. ConclusionThe additional use of TMM alongside MAR software yielded statistically significant reductions in CI induced metal artefacts on both water and anthropomorphic phantom scans with minimal dose increases. Implications for practiceThe extent of metal artefacts in clinical head scans with CI in place could be significantly reduced through combined use of TMM and MAR software, consequently providing greater diagnostic confidence in the images.

Improving the detection of hypo-vascular liver metastases in multiphase contrast-enhanced CT with slice thickness less than 5 mm using DenseNet

IntroductionThinner slices are more susceptible in detecting small lesions but suffer from higher statistical fluctuation. This work aimed to reduce image noise in multiphase contrast-enhanced CT reconstructed with slice thickness thinner than the clinical setting (i.e., 5 mm) using convolutional neural network (CNN) for enabling better detection of hypo-vascular liver metastasis. MethodsA DenseNet model was used to generate noise map for multiphase CT reconstructed with slice thickness of 2.5 mm and 1.25 mm. Image denoising was conducted by subtracting the CNN-generated noise map from CT images with reduced photon flux due to thinner slice thickness. The performance of DenseNet was evaluated on CT scans of electron density phantoms and patients with hypovascular liver metastases less than 1.5 cm in terms of Hounsfield Unit (HU) variation, statistical fluctuation, and contrast-to-noise ratio (CNR). ResultsThe phantom study demonstrated that the CNN-based denoising method was able to reduce statistical fluctuation in CT images reconstructed with slice thickness of 2.5 mm and 1.25 mm without causing significant edge blurring or variation in HU values. With regards to patient study, it was found that the denoised 2.5-mm and 1.25-mm slices had higher CNR than the conventional 5-mm slices for hypo-vascular liver metastases in all 4 phases of multiphase CT. ConclusionOur results demonstrated that the detection of hypo-vascular liver metastases in multiphase contrast-enhanced CT with slice thickness less than 5 mm could be improved by using the CNN-based denoising method. Implications for practiceReconstruction slice thickness has a strong influence on the image quality of CT imaging. A CNN-based denoising method was used in this work to reduce the image noise in multiphase contrast-enhanced CT reconstructed with slice thickness less than 5 mm.

A single-institution experience with 177Lu RPT workflow improvements and qualifying the SPECT/CT imaging for dosimetry.

Implementing any radiopharmaceutical therapy (RPT) program requires a comprehensive review of system readiness, appropriate workflows, and training to ensure safe and efficient treatment delivery. A quantitative assessment of the dose delivered to targets and organs at risk (OAR) using RPT is possible by correlating the absorbed doses with the delivered radioactivity. Integrating dosimetry into an established RPT program demands a thorough analysis of the necessary components and system fine-tuning. This study aims to report an optimized workflow for molecular radiation therapy using 177Lu with a primary focus on integrating patient-specific dosimetry into an established radiopharmaceutical program in a radiation oncology setting. We comprehensively reviewed using the Plan-Do-Check-Act (PDCA) cycle, including efficacy and accuracy of delivery and all aspects of radiation safety of the RPT program. The GE Discovery SPECT/CT 670DR™ system was calibrated per MIM protocol for dose calculation on MIM SurePlan™ MRT software. Jaszcak Phantom with 15-20 mCi of 177Lu DOTATATE with 2.5 µM EDTA solution was used, with the main energy window defined as 208 keV ±10% (187.6 to 229.2 keV); the upper scatter energy window was set to 240 keV ±5% (228 to 252 keV), while the lower scatter energy window was 177.8 keV ±5% (168.9 to 186.7 keV). Volumetric quality control tests and adjustments were performed to ensure the correct alignment of the table, NM, and CT gantry on SPECT/CT. A comprehensive end-to-end (E2E) test was performed to ensure workflow, functionality, and quantitative dose accuracy. Workflow improvements and checklists are presented after systematically analyzing over 400 administrations of 177Lu-based RPT. Injected activity to each sphere in the NEMA Phantom scan was quantified, and the MIM Sureplan MRT reconstruction images calculated activities within ±12% of the injected activity. Image alignment tests on the SPECT/CT showed a discrepancy of more than the maximum tolerance of 2.2 mm on any individual axis. As a result of servicing the machine and updating the VQC and COR corrections, the hybrid imaging system was adjusted to achieve an accuracy of <1 mm in all directions. Workflows and checklists, after analysis of system readiness and adequate training for staff and patients, are presented. Hardware and software components for patient-specific dosimetry are presented with a focus on hybrid image registration and correcting any errors that affect dosimetric quantification calculation. Moreover, this manuscript briefly overviews the necessary quality assurance requirements for converting diagnostic images into dosimetry measurement tools and integrating dosimetry for RPT based on 177Lu.

Harmonizing multisite data with the ComBat method for enhanced Parkinson's disease diagnosis via DAT-SPECT.

Dopamine transporter single-photon emission computed tomography (DAT-SPECT) is a crucial tool for evaluating patients with Parkinson's disease (PD). However, its implication is limited by inter-site variability in large multisite clinical trials. To overcome the limitation, a conventional prospective correction method employs linear regression with phantom scanning, which is effective yet available only in a prospective manner. An alternative, although relatively underexplored, involves retrospective modeling using a statistical method known as "combatting batch effects when combining batches of gene expression microarray data" (ComBat). We analyzed DAT-SPECT-specific binding ratios (SBRs) derived from 72 healthy older adults and 81 patients with PD registered in four clinical sites. We applied both the prospective correction and the retrospective ComBat correction to the original SBRs. Next, we compared the performance of the original and two corrected SBRs to differentiate the PD patients from the healthy controls. Diagnostic accuracy was assessed using the area under the receiver operating characteristic curve (AUC-ROC). The original SBRs were 6.13 ± 1.54 (mean ± standard deviation) and 2.03 ± 1.41 in the control and PD groups, respectively. After the prospective correction, the mean SBRs were 6.52 ± 1.06 and 2.40 ± 0.99 in the control and PD groups, respectively. After the retrospective ComBat correction, the SBRs were 5.25 ± 0.89 and 2.01 ± 0.73 in the control and PD groups, respectively, resulting in substantial changes in mean values with fewer variances. The original SBRs demonstrated fair performance in differentiating PD from controls (Hedges's g = 2.76; AUC-ROC = 0.936). Both correction methods improved discrimination performance. The ComBat-corrected SBR demonstrated comparable performance (g = 3.99 and AUC-ROC = 0.987) to the prospectively corrected SBR (g = 4.32 and AUC-ROC = 0.992) for discrimination. Although we confirmed that SBRs fairly discriminated PD from healthy older adults without any correction, the correction methods improved their discrimination performance in a multisite setting. Our results support the utility of harmonization methods with ComBat for consolidating SBR-based diagnosis or stratification of PD in multisite studies. Nonetheless, given the substantial changes in the mean values of ComBat-corrected SBRs, caution is advised when interpreting them.