Fat Quantification in Dual-Layer Detector Spectral Computed Tomography: Experimental Development and First In-Patient Validation.

OBJECTIVES: Fat quantification by dual-energy computed tomography (DECT) provides contrast-independent objective results, for example, on hepatic steatosis or muscle quality as parameters of prognostic relevance. To date, fat quantification has only been developed and used for source-based DECT techniques as fast kVp-switching CT or dual-source CT, which require a prospective selection of the dual-energy imaging mode.It was the purpose of this study to develop a material decomposition algorithm for fat quantification in phantoms and validate it in vivo for patient liver and skeletal muscle using a dual-layer detector-based spectral CT (dlsCT), which automatically generates spectral information with every scan. MATERIALS AND METHODS: For this feasibility study, phantoms were created with 0%, 5%, 10%, 25%, and 40% fat and 0, 4.9, and 7.0 mg/mL iodine, respectively. Phantom scans were performed with the IQon spectral CT (Philips, the Netherlands) at 120 kV and 140 kV and 3 T magnetic resonance (MR) (Philips, the Netherlands) chemical-shift relaxometry (MRR) and MR spectroscopy (MRS). Based on maps of the photoelectric effect and Compton scattering, 3-material decomposition was done for fat, iodine, and phantom material in the image space.After written consent, 10 patients (mean age, 55 ± 18 years; 6 men) in need of a CT staging were prospectively included. All patients received contrast-enhanced abdominal dlsCT scans at 120 kV and MR imaging scans for MRR. As reference tissue for the liver and the skeletal muscle, retrospectively available non-contrast-enhanced spectral CT data sets were used. Agreement between dlsCT and MR was evaluated for the phantoms, 3 hepatic and 2 muscular regions of interest per patient by intraclass correlation coefficients (ICCs) and Bland-Altman analyses. RESULTS: The ICC was excellent in the phantoms for both 120 kV and 140 kV (dlsCT vs MRR 0.98 [95% confidence interval (CI), 0.94-0.99]; dlsCT vs MRS 0.96 [95% CI, 0.87-0.99]) and in the skeletal muscle (0.96 [95% CI, 0.89-0.98]). For log-transformed liver fat values, the ICC was moderate (0.75 [95% CI, 0.48-0.88]). Bland-Altman analysis yielded a mean difference of -0.7% (95% CI, -4.5 to 3.1) for the liver and of 0.5% (95% CI, -4.3 to 5.3) for the skeletal muscle. Interobserver and intraobserver agreement were excellent (>0.9). CONCLUSIONS: Fat quantification was developed for dlsCT and agreement with MR techniques demonstrated for patient liver and muscle. Hepatic steatosis and myosteatosis can be detected in dlsCT scans from clinical routine, which retrospectively provide spectral information independent of the imaging mode.

Efficient Joint Estimation of Tracer Distribution and Background Signals in Magnetic Particle Imaging Using a Dictionary Approach.

Background signals are a primary source of artifacts in magnetic particle imaging and limit the sensitivity of the method since background signals are often not precisely known and vary over time. The state-of-the art method for handling background signals uses one or several background calibration measurements with an empty scanner bore and subtracts a linear combination of these background measurements from the actual particle measurement. This approach yields satisfying results in case that the background measurements are taken in close proximity to the particle measurement and when the background signal drifts linearly. In this work, we propose a joint estimation of particle distribution and background signal based on a dictionary that is capable of representing typical background signals. Reconstruction is performed frame-by-frame with minimal assumptions on the temporal evolution of background signals. Thus, even non-linear temporal evolution of the latter can be captured. Using a singular-value decomposition, the dictionary is derived from a large number of background calibration scans that do not need to be recorded in close proximity to the particle measurement. The dictionary is sufficiently expressive and represented by its principle components. The proposed joint estimation of particle distribution and background signal is expressed as a linear Tikhonov-regularized least squares problem, which can be efficiently solved. In phantom experiments it is shown that the method strongly suppresses background artifacts and even allows to estimate and remove the direct feed-through of the excitation field.

MRIReco.jl: An MRI reconstruction framework written in Julia.

PURPOSE: The aim of this work is to develop a high-performance, flexible, and easy-to-use MRI reconstruction framework using the scientific programming language Julia. METHODS: Julia is a modern, general purpose programming language with strong features in the area of signal/image processing and numerical computing. It has a high-level syntax but still generates efficient machine code that is usually as fast as comparable C/C++ applications. In addition to the language features itself, Julia has a sophisticated package management system that makes proper modularization of functionality across different packages feasible. Our developed MRI reconstruction framework MRIReco.jl can therefore reuse existing functionality from other Julia packages and concentrate on the MRI-related parts. This includes common imaging operators and support for MRI raw data formats. RESULTS: MRIReco.jl is a simple to use framework with a high degree of accessibility. While providing a simple-to-use interface, many of its components can easily be extended and customized. The performance of MRIReco.jl is compared to the Berkeley Advanced Reconstruction Toolbox (BART) and we show that the Julia framework achieves comparable reconstruction speed as the popular C/C++ library. CONCLUSIONS: Modern programming languages can bridge the gap between high performance and accessible implementations. MRIReco.jl leverages this fact and contributes a promising environment for future algorithmic development in MRI reconstruction.