Data truncation is a common problem in computed tomography (CT). Truncation causes cupping artifacts inside the field-of-view (FOV) and anatomical structures missing outside the FOV. Deep learning has achieved impressive results in CT reconstruction from limited data. However, its robustness is still a concern for clinical applications. Although the image quality of learning-based compensation schemes may be inadequate for clinical diagnosis, they can provide prior information for more accurate extrapolation than conventional heuristic extrapolation methods. With extrapolated projection, a conventional image reconstruction algorithm can be applied to obtain a final reconstruction. In this work, a general plug-and-play (PnP) method for truncation correction is proposed based on this idea, where various deep learning methods and conventional reconstruction algorithms can be plugged in. Such a PnP method integrates data consistency for measured data and learned prior image information for truncated data. This shows to have better robustness and interpretability than deep learning only. To demonstrate the efficacy of the proposed PnP method, two state-of-the-art deep learning methods, FBPConvNet and Pix2pixGAN, are investigated for truncation correction in cone-beam CT in noise-free and noisy cases. Their robustness is evaluated by showing false negative and false positive lesion cases. With our proposed PnP method, false lesion structures are corrected for both deep learning methods. For FBPConvNet, the root-mean-square error (RMSE) inside the FOV can be improved from 92HU to around 30HU by PnP in the noisy case. Pix2pixGAN solely achieves better image quality than FBPConvNet solely for truncation correction in general. PnP further improves the RMSE inside the FOV from 42HU to around 27HU for Pix2pixGAN. The efficacy of PnP is also demonstrated on real clinical head data.
18 January 2024 in American Journal of Neuroradiology
BACKGROUND AND PURPOSE: Potential utility of flat panel CT perfusion imaging (FPCT-PI) performed immediately after mechanical thrombectomy (MT) is unknown. We aimed to assess whether FPCT-PI obtained directly post-MT could provide additional potentially relevant information on tissue reperfusion status. MATERIALS AND METHODS: This was a single-center analysis of all patients with consecutive acute stroke admitted between June 2019 and March 2021 who underwent MT and postinterventional FPCT-PI (n = 26). A core lab blinded to technical details and clinical data performed TICI grading on postinterventional DSA images and qualitatively assessed reperfusion on time-sensitive FPCT-PI maps. According to agreement between DSA and FPCT-PI, all patients were classified into 4 groups: hypoperfusion findings perfectly matched by location (group 1), hypoperfusion findings mismatched by location (group 2), complete reperfusion on DSA with hypoperfusion on FPCT-PI (group 3), and hypoperfusion on DSA with complete reperfusion on FPCT-PI (group 4). RESULTS: Detection of hypoperfusion (present/absent) concurred in 21/26 patients. Of these, reperfusion findings showed perfect agreement on location and size in 16 patients (group 1), while in 5 patients there was a mismatch by location (group 2). Of the remaining 5 patients with disagreement regarding the presence or absence of hypoperfusion, 3 were classified into group 3 and 2 into group 4. FPCT-PI findings could have avoided TICI overestimation in all false-positive operator-rated TICI 3 cases (10/26). CONCLUSIONS: FPCT-PI may provide additional clinically relevant information in a considerable proportion of patients undergoing MT. Hence, FPCT-PI may complement the evaluation of reperfusion efficacy and potentially inform decision-making in the angiography suite.
15 June 2020 in IEEE Transactions on Medical Imaging
In tomographic imaging, anatomical structures are reconstructed by applying a pseudo-inverse forward model to acquired signals. Geometric information within this process is usually depending on the system setting only, i. e., the scanner position or readout direction. Patient motion therefore corrupts the geometry alignment in the reconstruction process resulting in motion artifacts. We propose an appearance learning approach recognizing the structures of rigid motion independently from the scanned object. To this end, we train a siamese triplet network to predict the reprojection error (RPE) for the complete acquisition as well as an approximate distribution of the RPE along the single views from the reconstructed volume in a multi-task learning approach. The RPE measures the motioninduced geometric deviations independent of the object based on virtual marker positions, which are available during training. We train our network using 27 patients and deploy a 21- 4-2 split for training, validation and testing. In average, we achieve a residual mean RPE of 0.013mm with an inter-patient standard deviation of 0.022 mm. This is twice the accuracy compared to previously published results. In a motion estimation benchmark the proposed approach achieves superior results in comparison with two state-of-the-art measures in nine out of twelve experiments. The clinical applicability of the proposed method is demonstrated on a motion-affected clinical dataset.
4 October 2024 in IEEE Transactions on Medical Imaging
Cone-beam computed tomography (CBCT) systems, with their flexibility, present a promising avenue for direct point-of-care medical imaging, particularly in critical scenarios such as acute stroke assessment. However, the integration of CBCT into clinical workflows faces challenges, primarily linked to long scan duration resulting in patient motion during scanning and leading to image quality degradation in the reconstructed volumes. This paper introduces a novel approach to CBCT motion estimation using a gradient-based optimization algorithm, which leverages generalized derivatives of the backprojection operator for cone-beam CT geometries. Building on that, a fully differentiable target function is formulated which grades the quality of the current motion estimate in reconstruction space. We drastically accelerate motion estimation yielding a 19-fold speed-up compared to existing methods. Additionally, we investigate the architecture of networks used for quality metric regression and propose predicting voxel-wise quality maps, favoring autoencoder-like architectures over contracting ones. This modification improves gradient flow, leading to more accurate motion estimation. The presented method is evaluated through realistic experiments on head anatomy. It achieves a reduction in reprojection error from an initial average of 3 mm to 0.61 mm after motion compensation and consistently demonstrates superior performance compared to existing approaches. The analytic Jacobian for the backprojection operation, which is at the core of the proposed method, is made publicly available. In summary, this paper contributes to the advancement of CBCT integration into clinical workflows by proposing a robust motion estimation approach that enhances efficiency and accuracy, addressing critical challenges in time-sensitive scenarios.
Treated cerebral aneurysms (IA) require follow-up imaging to ensure occlusion. Metal artifacts complicate radiologic assessment. Our aim was to evaluate an innovative metal-artifact-reduction (iMAR) algorithm for flat-detector computed tomography angiography (FD-CTA) regarding image quality (IQ) and detection of aneurysm residua/reperfusion in comparison to 2D digital subtraction angiography (DSA). Patients with IAs treated by endovascular coiling or clipping underwent both FD-CTA and DSA. FD-CTA datasets were postprocessed with/without iMAR algorithm (MAR+/MAR−). Evaluation of all FD-CTA and DSA datasets regarding qualitative (IQ, MAR) and quantitative (coil package diameter/CPD) parameters was performed. Aneurysm occlusion was assessed for each dataset and compared to DSA findings. In total, 40 IAs were analyzed (ncoiling = 24; nclipping = 16). All iMAR+ datasets demonstrated significantly better IQ (pIQcoiling < 0.0001; pIQclipping < 0.0001). iMAR significantly reduced the metal-artifact burden but did not affect the CPD. iMAR significantly improved the detection of aneurysm residua/reperfusion with excellent agreement with DSA (naneurysmdetectionMAR+/MAR−/DSA = 22/1/26). The iMAR algorithm significantly improves IQ by effective reduction of metal artifacts in FD-CTA datasets. The proposed algorithm enables reliable detection of aneurysm residua/reperfusion with good agreement to DSA. Thus, iMAR can help to reduce the need for invasive follow-up in treated IAs.
Purpose: Flat-detector computed tomography (FD-CT) is the standard for cochlear implant (CI) imaging. FD-CT systems differ in technical characteristics. Our aim was an evaluation of two different FD-CT generations with different protocols and hardware regarding image quality, radiation dose, and scan time. Methods: Two temporal bone specimens (− / + CI = TB0/TB1) were scanned using three different scanners: two FD-CT systems with different scanning protocols (standard FD-CT: 20 s 70 kV, 20 s 109 kV; high-speed FD-CT [HS-FD-CT]: 7 s 109 kV, 9 s 109 kV, 14 s 72 kV) and MS-CT (5 s 120 kV). Acquired datasets were evaluated in consensus reading regarding qualitative and quantitative parameters: addressing CI- and cochlea-specific parameters, cochlea delineation, lamina spiralis ossea visibility, distinction of single CI electrodes, determination of intracochlear implant position, stapes delineation, and mastoidal septation were assessed. Addressing protocol-specific parameters, radiation dose (dose-length-product/DLP), and scan time were assessed. Results: Two HS-FD-CT protocols (14 s/9 s) provide higher or equivalent diagnostic information regarding CI- and cochlea-specific parameters compared to both standard FD-CT protocols. The fastest HS-FD-CT protocol (7 s)—providing inferior diagnostic information compared to all other FD-CT protocols—still exceeds MS-CT. The highest DLP was recorded for the 14 s HS-FD-CT protocol (TB1 = 956 mGycm); the lowest DLPs were recorded for the 7 s HS-FD-CT protocol (TB0 = 188 mGycm) and for MS-CT (TB0 = 138 mGycm), respectively. HS-FD-CT allows a significant reduction of scan time compared to standard FD-CT. Conclusion: High-speed FD-CT improves visualization of temporal bone anatomy and postoperative assessment of CIs by combining excellent image quality, fast scan time, and reasonable radiation exposure.