Dataset Information

Accelerated barrier optimization compressed sensing (ABOCS) reconstruction for cone-beam CT: phantom studies.

ABSTRACT: Recent advances in compressed sensing (CS) enable accurate CT image reconstruction from highly undersampled and noisy projection measurements, due to the sparsifiable feature of most CT images using total variation (TV). These novel reconstruction methods have demonstrated advantages in clinical applications where radiation dose reduction is critical, such as onboard cone-beam CT (CBCT) imaging in radiation therapy. The image reconstruction using CS is formulated as either a constrained problem to minimize the TV objective within a small and fixed data fidelity error, or an unconstrained problem to minimize the data fidelity error with TV regularization. However, the conventional solutions to the above two formulations are either computationally inefficient or involved with inconsistent regularization parameter tuning, which significantly limit the clinical use of CS-based iterative reconstruction. In this paper, we propose an optimization algorithm for CS reconstruction which overcomes the above two drawbacks.The data fidelity tolerance of CS reconstruction can be well estimated based on the measured data, as most of the projection errors are from Poisson noise after effective data correction for scatter and beam-hardening effects. We therefore adopt the TV optimization framework with a data fidelity constraint. To accelerate the convergence, we first convert such a constrained optimization using a logarithmic barrier method into a form similar to that of the conventional TV regularization based reconstruction but with an automatically adjusted penalty weight. The problem is then solved efficiently by gradient projection with an adaptive Barzilai-Borwein step-size selection scheme. The proposed algorithm is referred to as accelerated barrier optimization for CS (ABOCS), and evaluated using both digital and physical phantom studies.ABOCS directly estimates the data fidelity tolerance from the raw projection data. Therefore, as demonstrated in both digital Shepp-Logan and physical head phantom studies, consistent reconstruction performances are achieved using the same algorithm parameters on scans with different noise levels and∕or on different objects. On the contrary, the penalty weight in a TV regularization based method needs to be fine-tuned in a large range (up to seven times) to maintain the reconstructed image quality. The improvement of ABOCS on computational efficiency is demonstrated in the comparisons with adaptive-steepest-descent-projection-onto-convex-sets (ASD-POCS), an existing CS reconstruction algorithm also using constrained optimization. ASD-POCS alternatively minimizes the TV objective using adaptive steepest descent (ASD) and the data fidelity error using projection onto convex sets (POCS). For similar image qualities of the Shepp-Logan phantom, ABOCS requires less computation time than ASD-POCS in MATLAB by more than one order of magnitude.We propose ABOCS for CBCT reconstruction. As compared to other published CS-based algorithms, our method has attractive features of fast convergence and consistent parameter settings for different datasets. These advantages have been demonstrated on phantom studies.

SUBMITTER: Niu T

PROVIDER: S-EPMC3412436 | biostudies-other | 2012 Jul

REPOSITORIES: biostudies-other

ACCESS DATA

Similar Datasets

Project description:Electron paramagnetic resonance (EPR) imaging has evolved as a promising tool to provide non-invasive assessment of tissue oxygenation levels. Due to the extremely short T2 relaxation time of electrons, single point imaging (SPI) is used in EPRI, limiting achievable spatial and temporal resolution. This presents a problem when attempting to measure changes in hypoxic state. In order to capture oxygen variation in hypoxic tissues and localize cycling hypoxia regions, an accelerated EPRI imaging method with minimal loss of information is needed.We present an image acceleration technique, partial Fourier compressed sensing (PFCS), that combines compressed sensing (CS) and partial Fourier reconstruction. PFCS augments the original CS equation using conjugate symmetry information for missing measurements. To further improve image quality in order to reconstruct low-resolution EPRI images, a projection onto convex sets (POCS)-based phase map and a spherical-sampling mask are used in the reconstruction process. The PFCS technique was used in phantoms and in vivo SCC7 tumor mice to evaluate image quality and accuracy in estimating O2 concentration.In both phantom and in vivo experiments, PFCS demonstrated the ability to reconstruct images more accurately with at least a 4-fold acceleration compared to traditional CS. Meanwhile, PFCS is able to better preserve the distinct spatial pattern in a phantom with a spatial resolution of 0.6mm. On phantoms containing Oxo63 solution with different oxygen concentrations, PFCS reconstructed linewidth maps that were discriminative of different O2 concentrations. Moreover, PFCS reconstruction of partially sampled data provided a better discrimination of hypoxic and oxygenated regions in a leg tumor compared to traditional CS reconstructed images.EPR images with an acceleration factor of four are feasible using PFCS with reasonable assessment of tissue oxygenation. The technique can greatly enhance EPR applications and improve our understanding cycling hypoxia. Moreover this technique can be easily extended to various MRI applications.

Project description:In computer vision, reference datasets from simulation and real outdoor scenes have been highly successful in promoting algorithmic development in stereo reconstruction. Endoscopic stereo reconstruction for surgical scenes gives rise to specific problems, including the lack of clear corner features, highly specular surface properties and the presence of blood and smoke. These issues present difficulties for both stereo reconstruction itself and also for standardised dataset production. Previous datasets have been produced using computed tomography (CT) or structured light reconstruction on phantom or ex vivo models. We present a stereo-endoscopic reconstruction validation dataset based on cone-beam CT (SERV-CT). Two ex vivo small porcine full torso cadavers were placed within the view of the endoscope with both the endoscope and target anatomy visible in the CT scan. Subsequent orientation of the endoscope was manually aligned to match the stereoscopic view and benchmark disparities, depths and occlusions are calculated. The requirement of a CT scan limited the number of stereo pairs to 8 from each ex vivo sample. For the second sample an RGB surface was acquired to aid alignment of smooth, featureless surfaces. Repeated manual alignments showed an RMS disparity accuracy of around 2 pixels and a depth accuracy of about 2 mm. A simplified reference dataset is provided consisting of endoscope image pairs with corresponding calibration, disparities, depths and occlusions covering the majority of the endoscopic image and a range of tissue types, including smooth specular surfaces, as well as significant variation of depth. We assessed the performance of various stereo algorithms from online available repositories. There is a significant variation between algorithms, highlighting some of the challenges of surgical endoscopic images. The SERV-CT dataset provides an easy to use stereoscopic validation for surgical applications with smooth reference disparities and depths covering the majority of the endoscopic image. This complements existing resources well and we hope will aid the development of surgical endoscopic anatomical reconstruction algorithms.

Project description:ObjectivesThe aims of this study were to develop an accelerated multiparametric magnetic resonance imaging method based on 3D-quantification using an interleaved Look-Locker acquisition sequence with a T2 preparation pulse (3D-QALAS) combined with compressed sensing (CS) and to evaluate the effect of CS on the quantitative mapping, tissue segmentation, and quality of synthetic images.Materials and methodsA magnetic resonance imaging system phantom, containing multiple compartments with standardized T1, T2, and proton density (PD) values; 10 healthy volunteers; and 12 patients with multiple sclerosis were scanned using the 3D-QALAS sequence with and without CS and conventional contrast-weighted imaging. The scan times of 3D-QALAS with and without CS were 5:56 and 11:11, respectively. For healthy volunteers, brain volumetry and myelin estimation were performed based on the measured T1, T2, and PD. For patients with multiple sclerosis, the mean T1, T2, PD, and the amount of myelin in plaques and contralateral normal-appearing white matter (NAWM) were measured. Simple linear regression analysis and Bland-Altman analysis were performed for each metric obtained from the datasets with and without CS. To compare overall image quality and structural delineations on synthetic and conventional contrast-weighted images, case-control randomized reading sessions were performed by 2 neuroradiologists in a blinded manner.ResultsThe linearity of both phantom and volunteer measurements in T1, T2, and PD values obtained with and without CS was very strong (R2 = 0.9901-1.000). The tissue segmentation obtained with and without CS also had high linearity (R2 = 0.987-0.999). The quantitative tissue values of the plaques and NAWM obtained with CS showed high linearity with those without CS (R2 = 0.967-1.000). There were no significant differences in overall image quality between synthetic contrast-weighted images obtained with and without CS (P = 0.17-0.99).ConclusionsMultiparametric imaging of the whole brain based on 3D-QALAS can be accelerated using CS while preserving tissue quantitative values, tissue segmentation, and quality of synthetic images.

Project description:PurposeCurrently, four-dimensional (4D) cone-beam computed tomography (CBCT) requires a 3-4 min full-fan scan to ensure usable image quality. Recent advancements in sparse-view 4D-CBCT reconstruction have opened the possibility to reduce scan time and dose. The aim of this study is to provide a common framework for systematically evaluating algorithms for 4D-CBCT reconstruction from a 1-min scan. Using this framework, the AAPM-sponsored SPARE Challenge was conducted in 2018 to identify and compare state-of-the-art algorithms.MethodsA clinically realistic CBCT dataset was simulated using patient CT volumes from the 4D-Lung database. The selected patients had multiple 4D-CT sessions, where the first 4D-CT was used as the prior CT, and the rest were used as the ground truth volumes for simulating CBCT projections. A GPU-based Monte Carlo tool was used to simulate the primary, scatter, and quantum noise signals. A total of 32 CBCT scans of nine patients were generated. Additional qualitative analysis was performed on a clinical Varian and clinical Elekta dataset to validate the simulation study. Participants were blinded from the ground truth, and were given 3 months to apply their reconstruction algorithms to the projection data. The submitted reconstructions were analyzed in terms of root-mean-squared-error (RMSE) and structural similarity index (SSIM) with the ground truth within four different region-of-interests (ROI) - patient body, lungs, planning target volume (PTV), and bony anatomy. Geometric accuracy was quantified as the alignment error of the PTV.ResultsTwenty teams participated in the challenge, with five teams completing the challenge. Techniques involved in the five methods included iterative optimization, motion-compensation, and deformation of the prior 4D-CT. All five methods rendered significant reduction in noise and streaking artifacts when compared to the conventional Feldkamp-Davis-Kress (FDK) algorithm. The RMS of the three-dimensional (3D) target registration error of the five methods ranged from 1.79 to 3.00 mm. Qualitative observations from the Varian and Elekta datasets mostly concur with those from the simulation dataset. Each of the methods was found to have its own strengths and weaknesses. Overall, the MA-ROOSTER method, which utilizes a 4D-CT motion model for temporal regularization, had the best and most consistent image quality and accuracy.ConclusionThe SPARE Challenge represents the first framework for systematically evaluating state-of-the-art algorithms for 4D-CBCT reconstruction from a 1-min scan. Results suggest the potential for reducing scan time and dose for 4D-CBCT. The challenge dataset and analysis framework are publicly available for benchmarking future reconstruction algorithms.

Project description:Objective.A dedicated cone-beam breast computed tomography (BCT) using a high-resolution, low-noise detector operating in offset-detector geometry has been developed. This study investigates the effects of varying detector offsets and image reconstruction algorithms to determine the appropriate combination of detector offset and reconstruction algorithm.Approach.Projection datasets (300 projections in 360°) of 30 breasts containing calcified lesions that were acquired using a prototype cone-beam BCT system comprising a 40 × 30 cm flat-panel detector with 1024 × 768 detector pixels were reconstructed using Feldkamp-Davis-Kress (FDK) algorithm and served as the reference. The projection datasets were retrospectively truncated to emulate cone-beam datasets with sinograms of 768×768 and 640×768 detector pixels, corresponding to 5 cm and 7.5 cm lateral offsets, respectively. These datasets were reconstructed using the FDK algorithm with appropriate weights and an ASD-POCS-based Fast, total variation-Regularized, Iterative, Statistical reconstruction Technique (FRIST), resulting in a total of 4 offset-detector reconstructions (2 detector offsets × 2 reconstruction methods). Signal difference-to-noise ratio (SDNR), variance, and full-width at half-maximum (FWHM) of calcifications in two orthogonal directions were determined from all reconstructions. All quantitative measurements were performed on images in units of linear attenuation coefficient (1/cm).Results.The FWHM of calcifications did not differ (P > 0.262) among reconstruction algorithms and detector formats, implying comparable spatial resolution. For a chosen detector offset, the FRIST algorithm outperformed FDK in terms of variance and SDNR (P < 0.0001). For a given reconstruction method, the 5 cm offset provided better results.Significance.This study indicates the feasibility of using the compressed sensing-based, FRIST algorithm to reconstruct sinograms from offset-detectors. Among the reconstruction methods and detector offsets studied, FRIST reconstructions corresponding to a 30 cm × 30 cm with 5 cm lateral offset, achieved the best performance. A clinical prototype using such an offset geometry has been developed and installed for clinical trials.

Project description:To demonstrate the feasibility of reconstructing a cone-beam CT (CBCT) image by deformably altering a prior fan-beam CT (FBCT) image such that it matches the anatomy portrayed in the CBCT projection data set.A prior FBCT image of the patient is assumed to be available as a source image. A CBCT projection data set is obtained and used as a target image set. A parametrized deformation model is applied to the source FBCT image, digitally reconstructed radiographs (DRRs) that emulate the CBCT projection image geometry are calculated and compared to the target CBCT projection data, and the deformation model parameters are adjusted iteratively until the DRRs optimally match the CBCT projection data set. The resulting deformed FBCT image is hypothesized to be an accurate representation of the patient's anatomy imaged by the CBCT system. The process is demonstrated via numerical simulation. A known deformation is applied to a prior FBCT image and used to create a synthetic set of CBCT target projections. The iterative projection matching process is then applied to reconstruct the deformation represented in the synthetic target projections; the reconstructed deformation is then compared to the known deformation. The sensitivity of the process to the number of projections and the DRR/CBCT projection mismatch is explored by systematically adding noise to and perturbing the contrast of the target projections relative to the iterated source DRRs and by reducing the number of projections.When there is no noise or contrast mismatch in the CBCT projection images, a set of 64 projections allows the known deformed CT image to be reconstructed to within a nRMS error of 1% and the known deformation to within a nRMS error of 7%. A CT image nRMS error of less than 4% is maintained at noise levels up to 3% of the mean projection intensity, at which the deformation error is 13%. At 1% noise level, the number of projections can be reduced to 8 while maintaining CT image and deformation errors of less than 4% and 13%, respectively. The method is sensitive to contrast mismatch between the simulated projections and the target projections when the soft-tissue contrast in the projections is low.By using prior knowledge available in a FBCT image, the authors show that a CBCT image can be iteratively reconstructed from a comparatively small number of projection images, thus saving acquisition time and reducing imaging dose. This will enable more frequent daily imaging during radiation therapy. Because the process preserves the CT numbers of the FBCT image, the resulting 3D image intensities will be more accurate than a CBCT image reconstructed via conventional backprojection methods. Reconstruction errors are insensitive to noise at levels beyond what would typically be found in CBCT projection data, but are sensitive to contrast mismatch errors between the CBCT projection data and the DRRs.

Dataset Information

Accelerated barrier optimization compressed sensing (ABOCS) reconstruction for cone-beam CT: phantom studies.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets