Feldkamp-type VOI reconstruction from super-short-scan cone-beam data

FBP-type CT reconstruction algorithms for triple-source circular trajectory with different scanning radii

Multi-source computed tomography (CT) imaging has unique technical advantages not only for dynamic objects, but also for large-size objects by designing its imaging scan mode. Using the triple-source fan-beam imaging scan mode under three circular trajectories with two different radii, we in this study developed and analyzed theoretically several exact reconstruction algorithms in terms of full-scan and short-scan for three sets of truncated projection data. This triple-source scan configuration in different radii cases is easier to be simulated by a single-source scan configuration in an industrial CT system. The proposed algorithms are based on the idea of filtering-back-projection (FBP) algorithm, and can reconstruct the large-size objects under the same CT devices. The developed algorithms avoid data rebinning and can provide exact and fast image reconstruction. The results of the numerical simulation based data analysis verified that new algorithms were accurate and effective.

Evaluation of an Analytic Reconstruction Method as a Platform for Spectral Cone-beam CT

With the fast development of photon counting detection techniques, spectral CT with a photon counting detector has attracted considerable attention by increasing energy-resolution to identify and discriminate materials. The conventional analytic reconstruction algorithms can be directly applied to reconstruct spectral CT images for each spectrum or energy bin. However, a comprehensive evaluation of analytic reconstruction algorithms for spectral CT has not been reported yet. This motivates us to evaluate the analytic spiral cone-beam CT algorithms and to provide a fair comparison platform for the state-of-the-art iterative spectral CT reconstruction algorithms. Considering the fact that a narrow energy bin has high noise which degrades the imaging quality of spectral CT, an adaptive maximum a posterior (MAP) projection restoration algorithm is first used to reduce the noise, and then a 2D/3D weighted spiral Feldkamp-Davis-Kress (FDK) algorithm is implemented to reconstruct the spectral CT images at different energy bins. Finally, the principle component analysis (PCA) is employed to render the spectral reconstruction results into a color space. Our numerical results show that the analytic reconstruction approach is fast and it can provide high spatial resolution, high contrast resolution and high signal-noise-ratio (SNR) under higher helical pitches. This makes it possible to serve as a platform to evaluate the state-of-the-art iterative spectral CT algorithms.

Analytic reconstruction algorithms for triple-source CT with horizontal data truncation

PURPOSE

This paper explores a triple-source imaging method with horizontal data truncation to enlarge the field of view (FOV) for big objects.

METHODS

The study is conducted by using theoretical analysis, mathematical deduction, and numerical simulations. The proposed algorithms are implemented in c + + and matlab. While the basic platform is constructed in matlab, the computationally intensive segments are coded in c + +, which are linked via a mex interface.

RESULTS

A triple-source circular scanning configuration with horizontal data truncation is developed, where three pairs of x-ray sources and detectors are unevenly distributed on the same circle to cover the whole imaging object. For this triple-source configuration, a fan-beam filtered backprojection-type algorithm is derived for truncated full-scan projections without data rebinning. The algorithm is also extended for horizontally truncated half-scan projections and cone-beam projections in a Feldkamp-type framework. Using their method, the FOV is enlarged twofold to threefold to scan bigger objects with high speed and quality. The numerical simulation results confirm the correctness and effectiveness of the developed algorithms.

CONCLUSIONS

The triple-source scanning configuration with horizontal data truncation cannot only keep most of the advantages of a traditional multisource system but also cover a larger FOV for big imaging objects. In addition, because the filtering is shift-invariant, the proposed algorithms are very fast and easily parallelized on graphic processing units.

CONRAD--a software framework for cone-beam imaging in radiology

PURPOSE

In the community of x-ray imaging, there is a multitude of tools and applications that are used in scientific practice. Many of these tools are proprietary and can only be used within a certain lab. Often the same algorithm is implemented multiple times by different groups in order to enable comparison. In an effort to tackle this problem, the authors created CONRAD, a software framework that provides many of the tools that are required to simulate basic processes in x-ray imaging and perform image reconstruction with consideration of nonlinear physical effects.

METHODS

CONRAD is a Java-based state-of-the-art software platform with extensive documentation. It is based on platform-independent technologies. Special libraries offer access to hardware acceleration such as OpenCL. There is an easy-to-use interface for parallel processing. The software package includes different simulation tools that are able to generate up to 4D projection and volume data and respective vector motion fields. Well known reconstruction algorithms such as FBP, DBP, and ART are included. All algorithms in the package are referenced to a scientific source.

RESULTS

A total of 13 different phantoms and 30 processing steps have already been integrated into the platform at the time of writing. The platform comprises 74.000 nonblank lines of code out of which 19% are used for documentation. The software package is available for download at http://conrad.stanford.edu. To demonstrate the use of the package, the authors reconstructed images from two different scanners, a table top system and a clinical C-arm system. Runtimes were evaluated using the RabbitCT platform and demonstrate state-of-the-art runtimes with 2.5 s for the 256 problem size and 12.4 s for the 512 problem size.

CONCLUSIONS

As a common software framework, CONRAD enables the medical physics community to share algorithms and develop new ideas. In particular this offers new opportunities for scientific collaboration and quantitative performance comparison between the methods of different groups.

The feasibility of an inverse geometry CT system with stationary source arrays

PURPOSE

Inverse geometry computed tomography (IGCT) has been proposed as a new system architecture that combines a small detector with a large, distributed source. This geometry can suppress cone-beam artifacts, reduce scatter, and increase dose efficiency. However, the temporal resolution of IGCT is still limited by the gantry rotation time. Large reductions in rotation time are in turn difficult due to the large source array and associated power electronics. We examine the feasibility of using stationary source arrays for IGCT in order to achieve better temporal resolution. We anticipate that multiple source arrays are necessary, with each source array physically separated from adjacent ones.

METHODS

Key feasibility issues include spatial resolution, artifacts, flux, noise, collimation, and system timing clashes. The separation between the different source arrays leads to missing views, complicating reconstruction. For the special case of three source arrays, a two-stage reconstruction algorithm is used to estimate the missing views. Collimation is achieved using a rotating collimator with a small number of holes. A set of equally spaced source spots are designated on the source arrays, and a source spot is energized when a collimator hole is aligned with it. System timing clashes occur when multiple source spots are scheduled to be energized simultaneously. We examine flux considerations to evaluate whether sufficient flux is available for clinical applications.

RESULTS

The two-stage reconstruction algorithm suppresses cone-beam artifacts while maintaining resolution and noise characteristics comparable to standard third generation systems. The residual artifacts are much smaller in magnitude than the cone-beam artifacts eliminated. A mathematical condition is given relating collimator hole locations and the number of virtual source spots for which system timing clashes are avoided. With optimization, sufficient flux may be achieved for many clinical applications.

CONCLUSIONS

IGCT with stationary source arrays could be an imaging platform potentially capable of imaging a complete 16-cm thick volume within a tenth of a second.

Volume-of-change cone-beam CT for image-guided surgery

C-arm cone-beam CT (CBCT) can provide intraoperative 3D imaging capability for surgical guidance, but workflow and radiation dose are the significant barriers to broad utilization. One main reason is that each 3D image acquisition requires a complete scan with a full radiation dose to present a completely new 3D image every time. In this paper, we propose to utilize patient-specific CT or CBCT as prior knowledge to accurately reconstruct the aspects of the region that have changed by the surgical procedure from only a sparse set of x-rays. The proposed methods consist of a 3D-2D registration between the prior volume and a sparse set of intraoperative x-rays, creating digitally reconstructed radiographs (DRRs) from the registered prior volume, computing difference images by subtracting DRRs from the intraoperative x-rays, a penalized likelihood reconstruction of the volume of change (VOC) from the difference images, and finally a fusion of VOC reconstruction with the prior volume to visualize the entire surgical field. When the surgical changes are local and relatively small, the VOC reconstruction involves only a small volume size and a small number of projections, allowing less computation and lower radiation dose than is needed to reconstruct the entire surgical field. We applied this approach to sacroplasty phantom data obtained from a CBCT test bench and vertebroplasty data with a fresh cadaver acquired from a C-arm CBCT system with a flat-panel detector. The VOCs were reconstructed from a varying number of images (10-66 images) and compared to the CBCT ground truth using four different metrics (mean squared error, correlation coefficient, structural similarity index and perceptual difference model). The results show promising reconstruction quality with structural similarity to the ground truth close to 1 even when only 15-20 images were used, allowing dose reduction by the factor of 10-20.

Four-dimensional volume-of-interest reconstruction for cone-beam computed tomography-guided radiation therapy

PURPOSE

Data sufficiency are a major problem in four-dimensional cone-beam computed tomography (4D-CBCT) on linear accelerator-integrated scanners for image-guided radiotherapy. Scan times must be in the range of 4-6 min to avoid undersampling artifacts. Various image reconstruction algorithms have been proposed to accommodate undersampled data acquisitions, but these algorithms are computationally expensive, may require long reconstruction times, and may require algorithm parameters to be optimized. The authors present a novel reconstruction method, 4D volume-of-interest (4D-VOI) reconstruction which suppresses undersampling artifacts and resolves lung tumor motion for undersampled 1-min scans. The 4D-VOI reconstruction is much less computationally expensive than other 4D-CBCT algorithms.

METHODS

The 4D-VOI method uses respiration-correlated projection data to reconstruct a four-dimensional (4D) image inside a VOI containing the moving tumor, and uncorrelated projection data to reconstruct a three-dimensional (3D) image outside the VOI. Anatomical motion is resolved inside the VOI and blurred outside the VOI. The authors acquired a 1-min. scan of an anthropomorphic chest phantom containing a moving water-filled sphere. The authors also used previously acquired 1-min scans for two lung cancer patients who had received CBCT-guided radiation therapy. The same raw data were used to test and compare the 4D-VOI reconstruction with the standard 4D reconstruction and the McKinnon-Bates (MB) reconstruction algorithms.

RESULTS

Both the 4D-VOI and the MB reconstructions suppress nearly all the streak artifacts compared with the standard 4D reconstruction, but the 4D-VOI has 3-8 times greater contrast-to-noise ratio than the MB reconstruction. In the dynamic chest phantom study, the 4D-VOI and the standard 4D reconstructions both resolved a moving sphere with an 18 mm displacement. The 4D-VOI reconstruction shows a motion blur of only 3 mm, whereas the MB reconstruction shows a motion blur of 13 mm. With graphics processing unit hardware used to accelerate computations, the 4D-VOI reconstruction required a 40-s reconstruction time.

CONCLUSIONS

4D-VOI reconstruction effectively reduces undersampling artifacts and resolves lung tumor motion in 4D-CBCT. The 4D-VOI reconstruction is computationally inexpensive compared with more sophisticated iterative algorithms. Compared with these algorithms, our 4D-VOI reconstruction is an attractive alternative in 4D-CBCT for reconstructing target motion without generating numerous streak artifacts.

Cone-beam artifact evaluation of the factorization method

PURPOSE

The authors investigate the CB artifact behavior of the factorization approach recently suggested for image reconstruction in circular cone-beam computed tomography. This investigation is carried out in a typical C-arm geometry and involves simulated data and for the first time also phantom and clinical CB data acquired with a commercially available angiographic system.

METHODS

The CB artifact level is first measured using quantitative figures-of-merit that are computed from the reconstructions of the mathematical FORBILD head phantom and of a modified disk phantom. The authors then show reconstructions from a physical thorax phantom and clinical head data sets for a visual assessment of image quality. The performance of the factorization method is primarily compared to that of short-scan FDK, but the authors also show the results obtained with the full-scan FDK and the virtual PI-line BPF method for the simulation studies, as a benchmark.

RESULTS

Quantitatively, the FORBILD head phantom reconstructions of both FDK methods show a spatially averaged bias of up to 1.2% in the axial slices about 9 cm away from the plane of the scan, which is placed 4 cm below the central slice through the phantom. The artifact level for the short-scan FDK method and the virtual PI-line BPF method noticeably depends on the scan orientation. The factorization approach can significantly reduce both, this dependency as well as the reconstruction bias. It also shows visually an improved quality of the clinical images compared to short-scan FDK, particularly close to the spine and in the subcranial regions of the clinical data sets.

CONCLUSIONS

The factorization approach comes with noticeably lower reconstruction bias than the FDK methods and is least sensitive to the scan orientation among all considered short-scan methods. The data inconsistencies contained in the real data sets, such as scatter, beam hardening, or data truncation, show only little impact on the factorization results. Hence, in both, reconstructions from real and simulated data, the factorization method yields better image quality than short-scan FDK, albeit at the cost of some slight, directed high-frequency artifacts that are mostly visible in axial slices.

Hybrid cone-beam tomographic reconstruction: incorporation of prior anatomical models to compensate for missing data

We propose a method for improving the quality of cone-beam tomographic reconstruction done with a C-arm. C-arm scans frequently suffer from incomplete information due to image truncation, limited scan length, or other limitations. Our proposed "hybrid reconstruction" method injects information from a prior anatomical model, derived from a subject-specific computed tomography (CT) or from a statistical database (atlas), where the C-arm X-ray data is missing. This significantly reduces reconstruction artifacts with little loss of true information from the X-ray projections. The methods consist of constructing anatomical models, fast rendering of digitally reconstructed radiograph (DRR) projections of the models, rigid or deformable registration of the model and the X-ray images, and fusion of the DRR and X-ray projections, all prior to a conventional filtered back-projection algorithm. Our experiments, conducted with a mobile image intensifier C-arm, demonstrate visually and quantitatively the contribution of data fusion to image quality, which we assess through comparison to a "ground truth" CT. Importantly, we show that a significantly improved reconstruction can be obtained from a C-arm scan as short as 90° by complementing the observed projections with DRRs of two prior models, namely an atlas and a preoperative same-patient CT. The hybrid reconstruction principles are applicable to other types of C-arms as well.

Temporally targeted imaging method applied to ECG-gated computed tomography: preliminary phantom and in vivo experience

RATIONALE AND OBJECTIVES

Existing cardiac imaging methods do not allow for improved temporal resolution when considering a targeted region of interest (ROI). The imaging method presented here enables improved temporal resolution for ROI imaging (namely, a reconstruction volume smaller than the complete field of view). Clinically, temporally targeted reconstruction would not change the primary means of reconstructing and evaluating images, but rather would enable the adjunct technique of ROI imaging, with improved temporal resolution compared with standard reconstruction ( approximately 20% smaller temporal scan window). In gated cardiac computed tomography (CT) scans improved temporal resolution directly translates into a reduction in motion artifacts for rapidly moving objects such as the coronary arteries.

MATERIALS AND METHODS

Retrospectively electrocardiogram gated coronary angiography data from a 64-slice CT system were used. A motion phantom simulating the motion profile of a coronary artery was constructed and scanned. Additionally, an in vivo study was performed using a porcine model. Comparisons between the new reconstruction technique and the standard reconstruction are given for an ROI centered on the right coronary artery, and a pulmonary ROI.

RESULTS

In both a well-controlled motion model and a porcine model results show a decrease in motion induced artifacts including motion blur and streak artifacts from contrast enhanced vessels within the targeted ROIs, as assessed through both qualitative and quantitative observations.

CONCLUSION

Temporally targeted reconstruction techniques demonstrate the potential to reduce motion artifacts in coronary CT. Further study is warranted to demonstrate the conditions under which this technique will offer direct clinical utility. Improvement in temporal resolution for gated cardiac scans has implications for improving: contrast enhanced CT angiography, calcium scoring, and assessment of the pulmonary anatomy.

A factorization approach for cone-beam reconstruction on a circular short-scan

In this paper, we introduce a new algorithm for 3-D image reconstruction from cone-beam (CB) projections acquired along a partial circular scan. Our algorithm is based on a novel, exact factorization of the initial 3-D reconstruction problem into a set of independent 2-D inversion problems, each of which corresponds to finding the object density on one, single plane. Any such 2-D inversion problem is solved numerically using a projected steepest descent iteration scheme. We present a numerical evaluation of our factorization algorithm using computer-simulated CB data, without and with noise, of the FORBILD head phantom and of a disk phantom. First, we study quantitatively the impact of the reconstruction parameters on the algorithm performance. Next, we present reconstruction results for visual assessment of the achievable image quality and provide, for comparison, results obtained with two other state-of-the-art reconstruction algorithms for the circular short-scan.

A new scheme for view-dependent data differentiation in fan-beam and cone-beam computed tomography

In computed tomography, analytical fan-beam (FB) and cone-beam (CB) image reconstruction often involves a view-dependent data differentiation. The implementation of this differentiation step is critical in terms of resolution and image quality. In this work, we present a new differentiation scheme that is robust to changes in the data acquisition geometry and to coarse view sampling. Our scheme was compared to two previously suggested methods, which we call the direct scheme and the chain-rule scheme. Image reconstructions were performed from computer-simulated data of the Shepp-Logan phantom, the FORBILD thorax phantom and a modified FORBILD head phantom. For FB reconstruction, we investigated three acquisition geometries: a circular, an ellipse-shaped and a square-shaped trajectory. For CB reconstruction, the circle-plus-line trajectory was considered. Image comparison showed that the new scheme performs consistently well when varying the scenario, in both FB and CB geometry, unlike the other two schemes.

Lambda tomography with discontinuous scanning trajectories

Lambda tomography (LT) is a well-known local reconstruction technology to reduce the radiation dose or accommodate a limited imaging geometry. After a theoretical analysis of the so-called Calderon operator (CO), the necessary conditions for exact LT reconstruction are presented in terms of the 2D and 3D COs. Based on our previous results on LT, a general scheme is proposed to construct exact LT formulae in terms of the 2D CO with multiple segment trajectories. Every 2D formula has a corresponding 3D cone-beam formula in the Feldkamp framework in terms of the 2D CO which was illustrated in a triple-segment case. Our simulation results verify the correctness and demonstrate the merits of the proposed scheme.

Approximate and exact cone-beam reconstruction with standard and non-standard spiral scanning

The long object problem is practically important and theoretically challenging. To solve the long object problem, spiral cone-beam CT was first proposed in 1991, and has been extensively studied since then. As a main feature of the next generation medical CT, spiral cone-beam CT has been greatly improved over the past several years, especially in terms of exact image reconstruction methods. Now, it is well established that volumetric images can be exactly and efficiently reconstructed from longitudinally truncated data collected along a rather general scanning trajectory. Here we present an overview of some key results in this area.

Data consistency based rigid motion artifact reduction in fan-beam CT

It is well known that a rigid in-plane motion can be decomposed into a translation and a rotation around an origin. Based on our previous work, we first extend the Helgason-Ludwig consistency condition (HLCC) to cover a general rigid motion in fan-beam geometry. Then, we model the general motion by several parameters, and develop an iterative scheme for estimation of the in-plane motion parameters. This scheme determines the motion parameters by numerically minimizing an objective function constructed based on the HLCC. After the motion parameters are estimated, image reconstruction can be performed to compensate for the motion effects. Finally, we implement the algorithm and evaluate its performance in numerical simulations.

Data consistency based translational motion artifact reduction in fan-beam CT

A basic assumption in the classic computed tomography (CT) theory is that an object remains stationary in an entire scan. In biomedical CT/micro-CT, this assumption is often violated. To produce high-resolution images, such as for our recently proposed clinical micro-CT (CMCT) prototype, it is desirable to develop a precise motion estimation and image reconstruction scheme. In this paper, we first extend the Helgason-Ludwig consistency condition (HLCC) from parallel-beam to fan-beam geometry when an object is subject to a translation. Then, we propose a novel method to estimate the motion parameters only from sinograms based on the HLCC. To reconstruct the moving object, we formulate two generalized fan-beam reconstruction methods, which are in filtered backprojection and backprojection filtering formats, respectively. Furthermore, we present numerical simulation results to show that our approach is accurate and robust.

A backprojection-filtration algorithm for nonstandard spiral cone-beam CT with an n-PI-window

For applications in bolus-chasing computed tomography (CT) angiography and electron-beam micro-CT, the backprojection-filtration (BPF) formula developed by Zou and Pan was recently generalized by Ye et al to reconstruct images from cone-beam data collected along a rather flexible scanning locus, including a nonstandard spiral. A major implication of the generalized BPF formula is that it can be applied for n-PI-window-based reconstruction in the nonstandard spiral scanning case. In this paper, we design an n-PI-window-based BPF algorithm, and report the numerical simulation results with the 3D Shepp-Logan phantom and Defrise disk phantom. The proposed BPF algorithm consists of three steps: cone-beam data differentiation, weighted backprojection and inverse Hilbert filtration. Our simulated results demonstrate the feasibility and merits of the proposed algorithm.

Design, analysis and simulation for development of the first clinical micro-CT scanner

In this article, we propose to develop the first clinical micro-CT (CMCT) system for human temporal bone imaging in vivo. This CMCT system consists of medical CT and micro-CT scanners either as separate components or in a combination, a cross-modality registration mechanism such as a facial surface scanner, and associated software. This system integrates the strengths of state-of-the-art medical CT and micro-CT techniques to achieve a spatial resolution that is much higher than currently available for inner ear imaging at acceptable dose levels. Our design, analysis, and simulation results demonstrate that the CMCT system is feasible for inner ear imaging and other clinical applications. For example, the CMCT system has the potential to improve the safety of guiding cochlear implant electrodes within the inner ear and assist the placement of inner ear microcatheters for delivery of gene modification therapy or administration of neurotrophic factors. Imaging of microarchitectures of the cancellous bone would be also an important application.