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: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: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