Software architecture for multi-bed FDK-based reconstruction in X-ray CT scanners

Most small-animal X-ray computed tomography (CT) scanners are based on cone-beam geometry with a flat-panel detector orbiting in a circular trajectory. Image reconstruction in these systems is usually performed by approximate methods based on the algorithm proposed by Feldkamp et al. (FDK). Besides the implementation of the reconstruction algorithm itself, in order to design a real system it is necessary to take into account numerous issues so as to obtain the best quality images from the acquired data. This work presents a comprehensive, novel software architecture for small-animal CT scanners based on cone-beam geometry with circular scanning trajectory. The proposed architecture covers all the steps from the system calibration to the volume reconstruction and conversion into Hounsfield units. It includes an efficient implementation of an FDK-based reconstruction algorithm that takes advantage of system symmetries and allows for parallel reconstruction using a multiprocessor computer. Strategies for calibration and artifact correction are discussed to justify the strategies adopted. New procedures for multi-bed misalignment, beam-hardening, and Housfield units calibration are proposed. Experiments with phantoms and real data showed the suitability of the proposed software architecture for an X-ray small animal CT based on cone-beam geometry.This work was partially funded by AMIT project from the CDTI CENIT program, TEC2007-64731, TEC2008-06715- C02-01, RD07/0014/2009, TRA2009 0175, RECAVA-RETIC, and RD09/0077/00087 (Ministerio de Ciencia e Inovación), and ARTEMIS S2009/DPI-1802 (Comunidad de Madrid).Publicad

Maximum Likelihood Estimation of Head Motion Using Epipolar Consistency

Open gantry C-arm systems that are placed within the interventional room enable 3-D imaging and guidance for stroke therapy without patient transfer. This can profit in drastically reduced time-totherapy, however, due to the interventional setting, the data acquisition is comparatively slow. Thus, involuntary patient motion needs to be estimated and compensated to achieve high image quality. Patient motion results in a misalignment of the geometry and the acquired image data. Consistency measures can be used to restore the correct mapping to compensate the motion. They describe constraints on an idealized imaging process which makes them also sensitive to beam hardening, scatter, truncation or overexposure. We propose a probabilistic approach based on the Student’s t-distribution to model image artifacts that affect the consistency measure without sourcing from motion

Correction of beam hardening artefacts in microtomography for samples imaged in containers

We explore the use of referenceless multi-material beam hardening correction methods, with an emphasis on maintaining data quality for real-world imaging of geologic materials with a view towards automation. In particular, we consider cases where the sample of interest is surrounded by a container of uniform material and propose a novel container-only pre-correction technique to allow automation of the segmentation process required for such correction methods. The effectiveness of the new technique is demonstrated using both simulated and experimental data

2D antiscatter grid and scatter sampling based CBCT pipeline for image guided radiation therapy

Poor tissue visualization and quantitative accuracy in CBCT is a major barrier in expanding the role of CBCT imaging from target localization to quantitative treatment monitoring and plan adaptations in radiation therapy sessions. To further improve image quality in CBCT, 2D antiscatter grid based scatter rejection was combined with a raw data processing pipeline and iterative image reconstruction. The culmination of these steps was referred as quantitative CBCT, qCBCT. qCBCT data processing steps include 2D antiscatter grid implementation, measurement based residual scatter, image lag, and beam hardening correction for offset detector geometry CBCT with a bow tie filter. Images were reconstructed with iterative image reconstruction to reduce image noise. To evaluate image quality, qCBCT acquisitions were performed using a variety of phantoms to investigate the effect of object size and its composition on image quality. qCBCT image quality was benchmarked against clinical CBCT and MDCT images. Addition of image lag and beam hardening correction to scatter suppression reduced HU degradation in qCBCT by 10 HU and 40 HU, respectively. When compared to gold standard MDCT, mean HU errors in qCBCT and clinical CBCT were 10 HU and 27 HU, respectively. HU inaccuracy due to change in phantom size was 22 HU and 85 HU in qCBCT and clinical CBCT images, respectively. With iterative reconstruction, contrast to noise ratio improved by a factor of 1.25 when compared to clinical CBCT protocols. Robust artifact and noise suppression in qCBCT images can reduce the image quality gap between CBCT and MDCT, improving the promise of qCBCT in fulfilling the tasks that demand high quantitative accuracy, such as CBCT based dose calculations and treatment response assessment in image guided radiation therapy

Artefact Reduction Methods for Iterative Reconstruction in Full-fan Cone Beam CT Radiotherapy Applications

A cone beam CT (CBCT) system acquires two-dimensional projection images of an imaging object from multiple angles in one single rotation and reconstructs the object geometry in three dimensions for volumetric visualization. It is mounted on most modern linear accelerators and is routinely used in radiotherapy to verify patient positioning, monitor patient contour changes throughout the course of treatment, and enable adaptive radiotherapy planning. Iterative image reconstruction algorithms use mathematical methods to iteratively solve the reconstruction problem. Iterative algorithms have demonstrated improvement in image quality and / or reduction in imaging dose over traditional filtered back-projection (FBP) methods. However, despite the advancement in computer technology and growing availability of open-source iterative algorithms, clinical implementation of iterative CBCT has been limited. This thesis does not report development of codes for new iterative image reconstruction algorithms. It focuses on bridging the gap between the algorithm and its implementation by addressing artefacts that are the results of imperfections from the raw projections and from the imaging geometry. Such artefacts can severely degrade image quality and cannot be removed by iterative algorithms alone. Practical solutions to solving these artefacts will be presented and this in turn will better enable clinical implementation of iterative CBCT reconstruction