Reduction of metal artefacts from bilateral hip prostheses during lower extremity computed tomography angiography: an experimental phantom study.

INTRODUCTION: Image quality reduction due to metallic artefacts is a significant challenge during vascular computed tomography (CT) imaging of the lower extremities in patients with hip prostheses. This study aims to analyse various reconstruction algorithms' ability to reduce metal artefacts due to two types of hip prostheses during lower extremity CT angiography examinations. METHODS: A pelvis phantom was fabricated with the insertion of a tube filled with contrast media to simulate the femoral artery, and the phantom was then CT scanned with and without hip prostheses. Multimodal images were acquired using different kilovoltage peak (kVp) settings and reconstructed with different algorithms, such as filtered back projection (FBP), iterative reconstruction (iDose4), iterative model-based reconstruction (IMR) and orthopaedic metal artefact reduction (O-MAR). Image quality was assessed based on image noise, signal-to-noise ratio (SNR) and Hounsfield unit (HU) deviation. RESULTS: The IMR approach significantly improved image quality compared to iDose4 and FBP. For the vascular region, O-MAR improves SNR by 5 ± 1, 23 ± 5 and 42 ± 9 for FBP, iDose4 and IMR respectively, and improves HU precision towards the baseline values by 49% and 83% for FBP and IMR, respectively. The noise reduction was 71% and 89% for FBP and IMR, and 57% for iDose4. O-MAR greatly enhances SNR corrections among the most severe artefacts, with 29 ± 1 and 43 ± 4 for FBP and IMR, compared to iDose4 by 37 ± 7. CONCLUSION: IMR combined with O-MAR could improve the CT angiography of the lower extremities of patients with a hip prosthesis.

Highly accelerated parameter mapping using model-based alternating reconstruction coupling fitting.

Objective.A model-based alternating reconstruction coupling fitting, termed Model-based Alternating Reconstruction COupling fitting (MARCO), is proposed for accurate and fast magnetic resonance parameter mapping.Approach.MARCO utilizes the signal model as a regularization by minimizing the bias between the image series and the signal produced by the suitable signal model based on iteratively updated parameter maps when reconstructing. The technique can incorporate prior knowledge of both image series and parameters by adding sparsity constraints. The optimization problem is decomposed into three subproblems and solved through three alternating steps involving reconstruction and nonlinear least-square fitting, which can produce both contrast-weighted images and parameter maps simultaneously.Main results.The algorithm is applied toT2mapping with extended phase graph algorithm integrated and validated on undersampled multi-echo spin-echo data from both phantom and in vivo sources. Compared with traditional compressed sensing and model-based methods, the proposed approach yields more accurateT2maps with more details at high acceleration factors.Significance.The proposed method provides a basic framework for quantitative MR relaxometry, theoretically applicable to all quantitative MR relaxometry. It has the potential to improve the diagnostic utility of quantitative imaging techniques.

Comparison of model-based versus deep learning-based image reconstruction for thin-slice T2-weighted spin-echo prostate MRI.

PURPOSE: To compare a previous model-based image reconstruction (MBIR) with a newly developed deep learning (DL)-based image reconstruction for providing improved signal-to-noise ratio (SNR) in high through-plane resolution (1 mm) T2-weighted spin-echo (T2SE) prostate MRI. METHODS: Large-area contrast and high-contrast spatial resolution of the reconstruction methods were assessed quantitatively in experimental phantom studies. The methods were next evaluated radiologically in 17 subjects at 3.0 Tesla for whom prostate MRI was clinically indicated. For each subject, the axial T2SE raw data were directed to MBIR and to the DL reconstruction at three vendor-provided levels: (L)ow, (M)edium, and (H)igh. Thin-slice images from the four reconstructions were compared using evaluation criteria related to SNR, sharpness, contrast fidelity, and reviewer preference. Results were compared using the Wilcoxon signed-rank test using Bonferroni correction, and inter-reader comparisons were done using the Cohen and Krippendorf tests. RESULTS: Baseline contrast and resolution in phantom studies were equivalent for all four reconstruction pathways as desired. In vivo, all three DL levels (L, M, H) provided improved SNR versus MBIR. For virtually, all other evaluation criteria DL L and M were superior to MBIR. DL L and M were evaluated as superior to DL H in fidelity of contrast. For 44 of the 51 evaluations, the DL M reconstruction was preferred. CONCLUSION: The deep learning reconstruction method provides significant SNR improvement in thin-slice (1 mm) T2SE images of the prostate while retaining image contrast. However, if taken to too high a level (DL High), both radiological sharpness and fidelity of contrast diminish.

Computational wave-based photoacoustic imaging through an unknown thick aberrating layer.

We introduce a physics-based computational reconstruction framework for non-invasive photoacoustic tomography through a thick aberrating layer. Our wave-based approach leverages an analytic formulation of diffraction to beamform a photoacoustic image, when the aberrating layer profile is known. When the profile of the aberrating layer is unknown, the same analytical formulation serves as the basis for an automatic-differentiation regularized optimization algorithm that simultaneously reconstructs both the profile of the aberrating layer and the optically absorbing targets. Results from numerical studies and proof-of-concept experiments show promise for fast beamforming that takes into account diffraction effects occurring in the propagation through thick, highly-aberrating layers.

Model-based reconstruction for looping-star MRI.

PURPOSE: The aim of this study was to develop a reconstruction method that more fully models the signals and reconstructs gradient echo (GRE) images without sacrificing the signal to noise ratio and spatial resolution, compared to conventional gridding and model-based image reconstruction method. METHODS: By modeling the trajectories for every spoke and simplifying the scenario to only echo-in and echo-out mixture, the approach explicitly models the overlapping echoes. After modeling the overlapping echoes with two system matrices, we use the conjugate gradient algorithm (CG-SENSE) with the nonuniform FFT (NUFFT) to optimize the image reconstruction cost function. RESULTS: The proposed method is demonstrated in phantoms and in-vivo volunteer experiments for three-dimensional, high-resolution T2*-weighted imaging and functional MRI tasks. Compared to the gridding method, the high resolution protocol exhibits improved spatial resolution and reduced signal loss as a result of less intra-voxel dephasing. The fMRI task shows that the proposed model-based method produced images with reduced artifacts and blurring as well as more stable and prominent time courses. CONCLUSION: The proposed model-based reconstruction results shows improved spatial resolution and reduced artifacts. The fMRI task shows improved time series and activation map due to the reduced overlapping echoes and under-sampling artifacts.

Model-based super-resolution reconstruction for pseudo-continuous Arterial Spin Labeling.

Arterial spin labeling (ASL) is a promising, non-invasive perfusion magnetic resonance imaging technique for quantifying cerebral blood flow (CBF). Unfortunately, ASL suffers from an inherently low signal-to-noise ratio (SNR) and spatial resolution, undermining its potential. Increasing spatial resolution without significantly sacrificing SNR or scan time represents a critical challenge towards routine clinical use. In this work, we propose a model-based super-resolution reconstruction (SRR) method with joint motion estimation that breaks the traditional SNR/resolution/scan-time trade-off. From a set of differently oriented 2D multi-slice pseudo-continuous ASL images with a low through-plane resolution, 3D-isotropic, high resolution, quantitative CBF maps are estimated using a Bayesian approach. Experiments on both synthetic whole brain phantom data, and on in vivo brain data, show that the proposed SRR Bayesian estimation framework outperforms state-of-the-art ASL quantification.

Motion-corrected model-based reconstruction for 2D myocardial T1 mapping.

PURPOSE: To allow for T1 mapping of the myocardium within 2.3 s for a 2D slice utilizing cardiac motion-corrected, model-based image reconstruction. METHODS: Golden radial data acquisition is continuously carried out for 2.3 s after an inversion pulse. In a first step, dynamic images are reconstructed which show both contrast changes due to T1 recovery and anatomical changes due to the heartbeat. An image registration algorithm with a signal model for T1 recovery is applied to estimate non-rigid cardiac motion. In a second step, estimated motion fields are applied during an iterative model-based T1 reconstruction. The approach was evaluated in numerical simulations, phantom experiments and in in-vivo scans in healthy volunteers. RESULTS: The accuracy of cardiac motion estimation was shown in numerical simulations with an average motion field error of 0.7 ± 0.6 mm for a motion amplitude of 5.1 mm. The accuracy of T1 estimation was demonstrated in phantom experiments, with no significant difference (p = 0.13) in T1 estimated by the proposed approach compared to an inversion-recovery reference method. In vivo, the proposed approach yielded 1.3 × 1.3 mm T1 maps with no significant difference (p = 0.77) in T1 and SDs in comparison to a cardiac-gated approach requiring 16 s scan time (i.e., seven times longer than the proposed approach). Cardiac motion correction improved the precision of T1 maps, shown by a 40% reduced SD. CONCLUSION: We have presented an approach that provides T1 maps of the myocardium in 2.3 s by utilizing both cardiac motion correction and model-based T1 reconstruction.

Regularized SUPER-CAIPIRINHA: Accelerating 3D variable flip-angle T1 mapping with accurate and efficient reconstruction.

PURPOSE: To propose an acceleration method for 3D variable flip-angle (VFA) T1 mapping based on a technique called shift undersampling improves parametric mapping efficiency and resolution (SUPER). METHODS: The proposed method incorporates strategies of SUPER, controlled aliasing in volumetric parallel imaging (CAIPIRINHA), and total variation-based regularization to accelerate 3D VFA T1 mapping. The k-space sampling grid of CAIPIRINHA is internally undersampled with SUPER along the contrast dimension. A proximal algorithm was developed to preserve the computational efficiency of SUPER in the presence of regularization. The regularized SUPER-CAIPIRINHA (rSUPER-CAIPIRINHA) was compared with low rank plus sparsity (L + S), reconstruction of principal component coefficient maps (REPCOM), and other SUPER-based methods via simulations and in vivo brain T1 mapping. The results were assessed quantitatively with NRMSE and structural similarity index measure (SSIM), and qualitatively by two experienced reviewers. RESULTS: rSUPER-CAIPIRINHA achieved a lower NRMSE and higher SSIM than L + S (0.11 ± 0.01 vs. 0.19 ± 0.03, p < 0.001; 0.66 ± 0.05 vs. 0.37 ± 0.03, p < 0.001) and REPCOM (0.16 ± 0.02, p < 0.001; 0.46 ± 0.04, p < 0.001). The reconstruction time of rSUPER-CAIPIRINHA was 6% of L + S and 2% of REPCOM. For the qualitative comparison, rSUPER-CAIPIRINHA showed improvement of overall image quality and reductions of artifacts and blurring, although with a lower apparent SNR. Compared with 2D SUPER-SENSE, rSUPER-CAIPIRINHA significantly reduced NRMSE (0.11 ± 0.01 vs. 0.23 ± 0.04, p < 0.001) and generated less noisy reconstructions. CONCLUSION: By incorporating SUPER, CAIPIRINHA, and regularization, rSUPER-CAIPIRINHA mitigated noise amplification, reduced artifacts and blurring, and achieved faster reconstructions compared with L + S and REPCOM. These advantages render 3D rSUPER-CAIPIRINHA VFA T1 mapping potentially useful for clinical applications.

Adaptive model-based Magnetic Resonance.

PURPOSE: Conventional sequences are static in nature, fixing measurement parameters in advance in anticipation of a wide range of expected tissue parameter values. We set out to design and benchmark a new, personalized approach-termed adaptive MR-in which incoming subject data is used to update and fine-tune the pulse sequence parameters in real time. METHODS: We implemented an adaptive, real-time multi-echo (MTE) experiment for estimating T2 s. Our approach combined a Bayesian framework with model-based reconstruction. It maintained and continuously updated a prior distribution of the desired tissue parameters, including T2 , which was used to guide the selection of sequence parameters in real time. RESULTS: Computer simulations predicted accelerations between 1.7- and 3.3-fold for adaptive multi-echo sequences relative to static ones. These predictions were corroborated in phantom experiments. In healthy volunteers, our adaptive framework accelerated the measurement of T2 for n-acetyl-aspartate by a factor of 2.5. CONCLUSION: Adaptive pulse sequences that alter their excitations in real time could provide substantial reductions in acquisition times. Given the generality of our proposed framework, our results motivate further research into other adaptive model-based approaches to MRI and MRS.

Towards quantitative in vivo dosimetry using x-ray acoustic computed tomography.

BACKGROUND: Radiation dosimetry is essential for radiation therapy (RT) to ensure that radiation dose is accurately delivered to the tumor. Despite its wide use in clinical intervention, the delivered radiation dose can only be planned and verified via simulation. This makes precision radiotherapy challenging while in-line verification of the delivered dose is still absent in the clinic. X-ray-induced acoustic computed tomography (XACT) has recently been proposed as an imaging tool for in vivo dosimetry. PURPOSE: Most of the XACT studies focus on localizing the radiation beam. However, it has not been studied for its potential for quantitative dosimetry. The aim of this study was to investigate the feasibility of using XACT for quantitative in vivo dose reconstruction during radiotherapy. METHODS: Varian Eclipse system was used to generate simulated uniform and wedged 3D radiation field with a size of 4 cm × $ \times \ $ 4 cm. In order to use XACT for quantitative dosimetry measurements, we have deconvoluted the effects of both the x-ray pulse shape and the finite frequency response of the ultrasound detector. We developed a model-based image reconstruction algorithm to quantify radiation dose in vivo using XACT imaging, and universal back-projection (UBP) reconstruction is used as comparison. The reconstructed dose was calibrated before comparing it to the percent depth dose (PDD) profile. Structural similarity index matrix (SSIM) and root mean squared error (RMSE) are used for numeric evaluation. Experimental signals were acquired from 4 cm × $ \times \ $ 4 cm radiation field created by Linear Accelerator (LINAC) at depths of 6, 8, and 10 cm beneath the water surface. The acquired signals were processed before reconstruction to achieve accurate results. RESULTS: Applying model-based reconstruction algorithm with non-negative constraints successfully reconstructed accurate radiation dose in 3D simulation study. The reconstructed dose matches well with the PDD profile after calibration in experiments. The SSIMs between the model-based reconstructions and initial doses are over 85%, and the RMSEs of model-based reconstructions are eight times lower than the UBP reconstructions. We have also shown that XACT images can be displayed as pseudo-color maps of acoustic intensity, which correspond to different radiation doses in the clinic. CONCLUSION: Our results show that the XACT imaging by model-based reconstruction algorithm is considerably more accurate than the dose reconstructed by UBP algorithm. With proper calibration, XACT is potentially applicable to the clinic for quantitative in vivo dosimetry across a wide range of radiation modalities. In addition, XACT's capability of real-time, volumetric dose imaging seems well-suited for the emerging field of ultrahigh dose rate "FLASH" radiotherapy.

Quantitative MRI by nonlinear inversion of the Bloch equations.

PURPOSE: Development of a generic model-based reconstruction framework for multiparametric quantitative MRI that can be used with data from different pulse sequences. METHODS: Generic nonlinear model-based reconstruction for quantitative MRI estimates parametric maps directly from the acquired k-space by numerical optimization. This requires numerically accurate and efficient methods to solve the Bloch equations and their partial derivatives. In this work, we combine direct sensitivity analysis and pre-computed state-transition matrices into a generic framework for calibrationless model-based reconstruction that can be applied to different pulse sequences. As a proof-of-concept, the method is implemented and validated for quantitative T 1 $$ {\mathrm{T}}_1 $$ and T 2 $$ {\mathrm{T}}_2 $$ mapping with single-shot inversion-recovery (IR) FLASH and IR bSSFP sequences in simulations, phantoms, and the human brain. RESULTS: The direct sensitivity analysis enables a highly accurate and numerically stable calculation of the derivatives. The state-transition matrices efficiently exploit repeating patterns in pulse sequences, speeding up the calculation by a factor of 10 for the examples considered in this work, while preserving the accuracy of native ordinary differential equations solvers. The generic model-based method reproduces quantitative results of previous model-based reconstructions based on the known analytical solutions for radial IR FLASH. For IR bSFFP it produces accurate T 1 $$ {\mathrm{T}}_1 $$ and T 2 $$ {\mathrm{T}}_2 $$ maps for the National Insitute of Standards and Technology (NIST) phantom in numerical simulations and experiments. Feasibility is also shown for human brain, although results are affected by magnetization transfer effects. CONCLUSION: By developing efficient tools for numerical optimizations using the Bloch equations as forward model, this work enables generic model-based reconstruction for quantitative MRI.

Model-based reconstructions for intravoxel incoherent motion and diffusion tensor imaging parameter map estimations.

Intravoxel incoherent motion (IVIM) imaging and diffusion tensor imaging (DTI) facilitate noninvasive quantification of tissue perfusion and diffusion. Both are promising biomarkers in various diseases and a combined acquisition is therefore desirable. This comes with challenges, including noisy parameter maps and long scan times, especially for the perfusion fraction f and pseudo-diffusion coefficient D*. A model-based reconstruction has the potential to overcome these challenges. As a first step, our goal was to develop a model-based reconstruction framework for IVIM and combined IVIM-DTI parameter estimation. The IVIM and IVIM-DTI models were implemented in the PyQMRI model-based reconstruction framework and validated with simulations and in vivo data. Commonly used voxel-wise nonlinear least-squares fitting was used as the reference. Simulations with the IVIM and IVIM-DTI models were performed with 100 noise realizations to assess accuracy and precision. Diffusion-weighted data were acquired for IVIM reconstruction in the liver (n = 5), as well as for IVIM-DTI in the kidneys (n = 5) and lower-leg muscles (n = 6) of healthy volunteers. The median and interquartile range (IQR) values of the IVIM and IVIM-DTI parameters were compared to assess bias and precision. With model-based reconstruction, the parameter maps exhibited less noise, which was most pronounced in the f and D* maps, both in the simulations and in vivo. The bias values in the simulations were comparable between model-based reconstruction and the reference method. The IQR was lower with model-based reconstruction compared with the reference for all parameters. In conclusion, model-based reconstruction is feasible for IVIM and IVIM-DTI and improves the precision of the parameter estimates, particularly for f and D* maps.

Fast and accurate T2 mapping using Bloch simulations and low-rank plus sparse matrix decomposition.

PURPOSE: MRI's T2 relaxation time is one of the key contrast mechanisms for clinical diagnosis and prognosis of pathologies. Mapping this relaxation time, however, involves extensive scan times, which are needed to collect quantitative data, thereby impeding its integration into clinical routine. This study employs a low-rank plus sparse (L + S) signal decomposition approach in order to reconstruct accurate T2-maps from highly undersampled multi-echo spin-echo (MESE) MRI data. METHODS: Two new algorithms are presented: the first uses standard L + S approach, where both L and S are iteratively updated. The second technique, dubbed SPArse and fixed RanK (SPARK), uses a fixed-rank L, under the assumption that most MESE information is found in the L component and that this rank can be pre-calculated. The utility of these new techniques is demonstrated on in vivo brain and calf data at x2 to x6 acceleration factors. RESULTS: Accelerated T2 maps showed improved accuracy compared to fully sampled ground truth maps, when using L + S and SPARK techniques vis-à-vis standard GRAPPA acceleration. CONCLUSION: SPARK provides accurate T2 maps with increased robustness to the selection of reconstruction parameters making it suitable to a wide range of applications and facilitating the use of quantitative T2 information in clinical settings.

Free-breathing myocardial T1 mapping using inversion-recovery radial FLASH and motion-resolved model-based reconstruction.

PURPOSE: To develop a free-breathing myocardial T 1 $$ {\mathrm{T}}_1 $$ mapping technique using inversion-recovery (IR) radial fast low-angle shot (FLASH) and calibrationless motion-resolved model-based reconstruction. METHODS: Free-running (free-breathing, retrospective cardiac gating) IR radial FLASH is used for data acquisition at 3T. First, to reduce the waiting time between inversions, an analytical formula is derived that takes the incomplete T 1 $$ {\mathrm{T}}_1 $$ recovery into account for an accurate T 1 $$ {\mathrm{T}}_1 $$ calculation. Second, the respiratory motion signal is estimated from the k-space center of the contrast varying acquisition using an adapted singular spectrum analysis (SSA-FARY) technique. Third, a motion-resolved model-based reconstruction is used to estimate both parameter and coil sensitivity maps directly from the sorted k-space data. Thus, spatiotemporal total variation, in addition to the spatial sparsity constraints, can be directly applied to the parameter maps. Validations are performed on an experimental phantom, 11 human subjects, and a young landrace pig with myocardial infarction. RESULTS: In comparison to an IR spin-echo reference, phantom results confirm good T 1 $$ {\mathrm{T}}_1 $$ accuracy, when reducing the waiting time from 5 s to 1 s using the new correction. The motion-resolved model-based reconstruction further improves T 1 $$ {\mathrm{T}}_1 $$ precision compared to the spatial regularization-only reconstruction. Aside from showing that a reliable respiratory motion signal can be estimated using modified SSA-FARY, in vivo studies demonstrate that dynamic myocardial T 1 $$ {\mathrm{T}}_1 $$ maps can be obtained within 2 min with good precision and repeatability. CONCLUSION: Motion-resolved myocardial T 1 $$ {\mathrm{T}}_1 $$ mapping during free-breathing with good accuracy, precision and repeatability can be achieved by combining inversion-recovery radial FLASH, self-gating and a calibrationless motion-resolved model-based reconstruction.

ADEPT: Accurate Diffusion Echo-Planar imaging with multi-contrast shoTs.

PURPOSE: To introduce a novel imaging and parameter estimation framework for accurate multi-shot diffusion MRI. THEORY AND METHODS: We propose a new framework called ADEPT (Accurate Diffusion Echo-Planar imaging with multi-contrast shoTs) that enables fast diffusion MRI by allowing diffusion contrast settings to change between shots in a multi-shot EPI acquisition (i.e., intra-scan modulation). The framework estimates diffusion parameter maps directly from the acquired intra-scan modulated k-space data, while simultaneously accounting for shot-to-shot phase inconsistencies. The performance of the estimation framework is evaluated using Monte Carlo simulation studies and in-vivo experiments and compared to that of reference methods that rely on parallel imaging for shot-to-shot phase correction. RESULTS: Simulation and real-data experiments show that ADEPT yields more accurate and more precise estimates of the diffusion metrics in multi-shot EPI data in comparison with the reference methods. CONCLUSION: ADEPT allows fast multi-shot EPI diffusion MRI without significantly degrading the accuracy and precision of the estimated diffusion maps.

Quantitative Dual-Energy Imaging of Bone Marrow Edema Using Multisource Cone-Beam CT with Model-Based Decomposition.

Purpose: We investigated the feasibility of dual-energy (DE) detection of bone marrow edema (BME) using a dedicated extremity cone-beam CT (CBCT) with a unique three-source x-ray unit. The sources can be operated at different energies to enable single-scan DE acquisitions. However, they are arranged parallel to the axis of rotation, resulting in incomplete sampling and precluding the application of DE projection-domain decompositions (PDD) for beam-hardening reduction. Therefore, we propose a novel combination of a model-based "one-step" DE two-material decomposition followed by a constrained image-domain change-of-basis to obtain virtual non-calcium (VNCa) images for BME detection. Methods: DE projections were obtained using an "alternating-kV" protocol by operating the peripheral two sources of the CBCT system at low-energy (60 kV, 0.105 mAs/frame) and the central source at high-energy (100 kV, 0.028 mAs/frame), for a total of 600 frames over 216° of gantry rotation. Projections were processed with detector lag, glare and fast Monte Carlo (MC)-based iterative scatter corrections. Model-based material decomposition (MBMD) was then implemented to obtain aluminum (Al) and polyethylene (PE) volume fraction images with minimal beam-hardening. Statistical ray weights in MBMD were modified to account for regions with highly oblique sampling by the peripheral sources. To generate the VNCa maps, image-domain decomposition (IDD) constrained by the volume conservation principle (VCP) was performed to convert the Al and PE MBMD images into volume fractions of water, fat and cortical bone. Accuracy of BME detection was evaluated using physical phantom data acquired on the multi-source extremity CBCT scanner. Results: The proposed framework estimated the volume of BME with ~10% error. The MC-based scatter corrections and the modified MBMD ray weights were essential to achieve such performance - the error without MC scatter corrections was >30%, whereas the uniformity of estimated VNCa images was 3x improved using the modified weights compared to the conventional weights. Conclusions: The proposed DE decomposition framework was able to overcome challenges of high scatter and incomplete sampling to achieve BME detection on a CBCT system with axially-distributed x-ray sources.

Aesthetic Reconstruction of Onco-surgical Mandibular Defects Using Free Fibular Flap with and without CAD/CAM Customized Osteotomy Guide: A Randomized Controlled Clinical Trial.

BACKGROUND: Reconstruction of mandibular defects following ablative surgery remains a challenge even for experienced surgeons. Virtual planning and guided surgery, including computer-aided design/computer-aided manufacturing (CAD/CAM), afford optimized ways by which to plan complex surgery. This study aimed to evaluate and compare aesthetic outcome and surgical efficiency of free fibular flap (FFF) with and without CAD/CAM customized osteotomy guide (COG) for reconstruction of onco-surgical mandibular defects. METHODS: Twenty-two patients indicated for segmental mandibulectomy were randomly assigned to either CAD/CAM with COG group or that without COG- Model based reconstruction (MB group) at a 1:1 ratio. Aesthetic outcomes were evaluated by means of morphometric assessment and comparison for each differential area (DAr) and angle (DAn) in the affected side to the contralateral side of the mandible using computerized digital imaging analysis (CDIA) based on the post-operative 3D CT-scan. Subjective evaluation was performed using the Visual Analogue Scale (VAS) and Patient's Satisfaction Score (PSS). Surgical efficiency was a secondary outcome and evaluated as total operative time and ischemia time. RESULTS: The mean sagittal DAr was significantly lower in the COG group (277.28 ± 127.05 vs. 398.67 ± 139.10 mm2, P = 0.045). Although there was an improvement in the axial DAr (147.61 ± 55.42 vs. 183.68 ± 72.85 mm2), the difference was not statistically significant (P = 0.206). The mean differences (Δ) in both sagittal and coronal DAn were significantly lower in the COG group than in the MB group (6.11 ± 3.46 and 1.77 ± 1.12° vs. 9.53 ± 4.17 and 3.44 ± 2.34°), respectively. There were no statistically significant differences in the axial DAn between the two groups (P = 0.386). The PSS was significantly higher in the COG group, reflecting better aesthetic satisfaction than in the MB group (P = 0.041). The total operation and ischemia time were significantly shorter in favor of the COG group with a mean of (562.91 ± 51.22, 97.55 ± 16.80 min vs. 663.55 ± 53.43, 172.45 ± 21.87 min), respectively. CONCLUSION: The CAD/CAM with COG is more reliable and highly valuable in enhancing aesthetic outcomes and surgical efficiency of mandibular reconstruction by FFF compared to that without COG (MB reconstruction). TRIAL REGISTRATION: This trial was registered at ClinicalTrials.gov . REGISTRATION NUMBER: NCT03757273. Registration date: 28/11/2018.

3D model-based super-resolution motion-corrected cardiac T1 mapping.

Objective. To provide 3D high-resolution cardiac T1 maps using model-based super-resolution reconstruction (SRR).Approach. Due to signal-to-noise ratio limitations and the motion of the heart during imaging, often 2D T1 maps with only low through-plane resolution (i.e. slice thickness of 6-8 mm) can be obtained. Here, a model-based SRR approach is presented, which combines multiple stacks of 2D acquisitions with 6-8 mm slice thickness and generates 3D high-resolution T1 maps with a slice thickness of 1.5-2 mm. Every stack was acquired in a different breath hold (BH) and any misalignment between BH was corrected retrospectively. The novelty of the proposed approach is the BH correction and the application of model-based SRR on cardiac T1 Mapping. The proposed approach was evaluated in numerical simulations and phantom experiments and demonstrated in four healthy subjects.Main results. Alignment of BH states was essential for SRR even in healthy volunteers. In simulations, respiratory motion could be estimated with an RMS error of 0.18 ± 0.28 mm. SRR improved the visualization of small structures. High accuracy and precision (average standard deviation of 69.62 ms) of the T1 values was ensured by SRR while the detectability of small structures increased by 40%.Significance. The proposed SRR approach provided T1 maps with high in-plane and high through-plane resolution (1.3 × 1.3 × 1.5-2 mm3). The approach led to improvements in the visualization of small structures and precise T1 values.

Model-based super-resolution reconstruction with joint motion estimation for improved quantitative MRI parameter mapping.

Quantitative Magnetic Resonance (MR) imaging provides reproducible measurements of biophysical parameters, and has become an essential tool in clinical MR studies. Unfortunately, 3D isotropic high resolution (HR) parameter mapping is hardly feasible in clinical practice due to prohibitively long acquisition times. Moreover, accurate and precise estimation of quantitative parameters is complicated by inevitable subject motion, the risk of which increases with scanning time. In this paper, we present a model-based super-resolution reconstruction (SRR) method that jointly estimates HR quantitative parameter maps and inter-image motion parameters from a set of 2D multi-slice contrast-weighted images with a low through-plane resolution. The method uses a Bayesian approach, which allows to optimally exploit prior knowledge of the tissue and noise statistics. To demonstrate its potential, the proposed SRR method is evaluated for a T1 and T2 quantitative mapping protocol. Furthermore, the method's performance in terms of precision, accuracy, and spatial resolution is evaluated using simulated as well as real brain imaging experiments. Results show that our proposed fully flexible, quantitative SRR framework with integrated motion estimation outperforms state-of-the-art SRR methods for quantitative MRI.

Embedded Quantitative MRI T1ρ Mapping Using Non-Linear Primal-Dual Proximal Splitting.

Quantitative MRI (qMRI) methods allow reducing the subjectivity of clinical MRI by providing numerical values on which diagnostic assessment or predictions of tissue properties can be based. However, qMRI measurements typically take more time than anatomical imaging due to requiring multiple measurements with varying contrasts for, e.g., relaxation time mapping. To reduce the scanning time, undersampled data may be combined with compressed sensing (CS) reconstruction techniques. Typical CS reconstructions first reconstruct a complex-valued set of images corresponding to the varying contrasts, followed by a non-linear signal model fit to obtain the parameter maps. We propose a direct, embedded reconstruction method for T1ρ mapping. The proposed method capitalizes on a known signal model to directly reconstruct the desired parameter map using a non-linear optimization model. The proposed reconstruction method also allows directly regularizing the parameter map of interest and greatly reduces the number of unknowns in the reconstruction, which are key factors in the performance of the reconstruction method. We test the proposed model using simulated radially sampled data from a 2D phantom and 2D cartesian ex vivo measurements of a mouse kidney specimen. We compare the embedded reconstruction model to two CS reconstruction models and in the cartesian test case also the direct inverse fast Fourier transform. The T1ρ RMSE of the embedded reconstructions was reduced by 37-76% compared to the CS reconstructions when using undersampled simulated data with the reduction growing with larger acceleration factors. The proposed, embedded model outperformed the reference methods on the experimental test case as well, especially providing robustness with higher acceleration factors.