GPU-based Iterative Cone Beam CT Reconstruction Using Tight Frame Regularization

X-ray imaging dose from serial cone-beam CT (CBCT) scans raises a clinical concern in most image guided radiation therapy procedures. It is the goal of this paper to develop a fast GPU-based algorithm to reconstruct high quality CBCT images from undersampled and noisy projection data so as to lower the imaging dose. For this purpose, we have developed an iterative tight frame (TF) based CBCT reconstruction algorithm. A condition that a real CBCT image has a sparse representation under a TF basis is imposed in the iteration process as regularization to the solution. To speed up the computation, a multi-grid method is employed. Our GPU implementation has achieved high computational efficiency and a CBCT image of resolution 512\times512\times70 can be reconstructed in ~5 min. We have tested our algorithm on a digital NCAT phantom and a physical Catphan phantom. It is found that our TF-based algorithm is able to reconstrct CBCT in the context of undersampling and low mAs levels. We have also quantitatively analyzed the reconstructed CBCT image quality in terms of modulation-transfer-function and contrast-to-noise ratio under various scanning conditions. The results confirm the high CBCT image quality obtained from our TF algorithm. Moreover, our algorithm has also been validated in a real clinical context using a head-and-neck patient case. Comparisons of the developed TF algorithm and the current state-of-the-art TV algorithm have also been made in various cases studied in terms of reconstructed image quality and computation efficiency.Comment: 24 pages, 8 figures, accepted by Phys. Med. Bio

Consensus image method for unknown noise removal

Noise removal has been, and it is nowadays, an important task in computer vision. Usually, it is a previous task preceding other tasks, as segmentation or reconstruction. However, for most existing denoising algorithms the noise model has to be known in advance. In this paper, we introduce a new approach based on consensus to deal with unknown noise models. To do this, different filtered images are obtained, then combined using multifuzzy sets and averaging aggregation functions. The final decision is made by using a penalty function to deliver the compromised image. Results show that this approach is consistent and provides a good compromise between filters.This work is supported by the European Commission under Contract No. 238819 (MIBISOC Marie Curie ITN). H. Bustince was supported by Project TIN 2010-15055 of the Spanish Ministry of Science

A pilot study for the digital replacement of a distorted dentition acquired by Cone Beam Computed Tomography (CBCT)

Abstract Introduction: Cone beam CT (CBCT) is becoming a routine imaging modality designed for the maxillofacial region. Imaging patients with intra-oral metallic objects cause streak artefacts. Artefacts impair any virtual model by obliterating the teeth. This is a major obstacle for occlusal registration and the fabrication of orthognathic wafers to guide the surgical correction of dentofacial deformities. Aims and Objectives: To develop a method of replacing the inaccurate CBCT images of the dentition with an accurate representation and test the feasibility of the technique in the clinical environment. Materials and Method: Impressions of the teeth are acquired and acrylic baseplates constructed on dental casts incorporating radiopaque registration markers. The appliances are fitted and a preoperative CBCT is performed. Impressions are taken of the dentition with the devices in situ and subsequent dental models produced. The models are scanned to produce a virtual model. Both images of the patient and the model are imported into a virtual reality software program and aligned on the virtual markers. This allows the alignment of the dentition without relying on the teeth for superimposition. The occlusal surfaces of the dentition can be replaced with the occlusal image of the model. Results: The absolute mean distance of the mesh between the markers in the skulls was in the region of 0.09mm ± 0.03mm; the replacement dentition had an absolute mean distance of about 0.24mm ± 0.09mm. In patients the absolute mean distance between markers increased to 0.14mm ± 0.03mm. It was not possible to establish the discrepancies in the patient’s dentition, since the original image of the dentition is inherently inaccurate. Conclusion: It is possible to replace the CBCT virtual dentition of cadaveric skulls with an accurate representation to create a composite skull. The feasibility study was successful in the clinical arena. This could be a significant advancement in the accuracy of surgical prediction planning, with the ultimate goal of fabrication of a physical orthognathic wafer using reverse engineering

Improved visualization of X-ray phase contrast volumetric data through artifact-free integrated differential images

Artifacts arising when differential phase images are integrated is a common problem to several X-ray phase-based experimental techniques. The combination of noise and insufficient sampling of the high-frequency differential phase signal leads to the formation of streak artifacts in the projections, translating into poor image quality in the tomography slices. In this work, we apply a non-iterative integration algorithm proven to reduce streak artifacts in planar (2D) images to a differential phase tomography scan. We report on how the reduction of streak artifacts in the projections improves the quality of the tomography slices, especially in the directions different from the reconstruction plane. Importantly, the method is compatible with large tomography datasets in terms of computation time

Advanced Image Reconstruction for Limited View Cone-Beam CT

In a standard CT acquisition, a high number of projections is obtained around the sample, generally covering an angular span of 360º. However, complexities may arise in some clinical scenarios such as surgery and emergency rooms or Intensive Care Units (ICUs) when the accessibility to the patient is limited due to the monitoring equipment attached. X-ray systems used in these cases are usually C-arms that only enable the acquisition of planar images within a limited angular range. Obtaining 3D images in these scenarios could be extremely interesting for diagnosis or image guided surgery. This would be based on the acquisition of a small number of projections within a limited angular span. Reconstruction of these limited-view data with conventional algorithms such as FDK result in streak artifacts and shape distortion deteriorating the image quality. In order to reduce these artifacts, advanced reconstruction methods can be used to compensate the lack of data by the incorporation of prior information. This bachelor thesis is framed on one of the lines of research carried out by the Biomedical Imaging and Instrumentation group from the Bioengineering and Aerospace Department of Universidad Carlos III de Madrid working jointly with the Hospital General Universitario Gregorio Marañón through its Instituto de Investigación Sanitaria. This line of research is carried out in collaboration with the company SEDECAL, which enables the direct transfer to the industry. Previous work showed that a new iterative reconstruction method proposed by the group, SCoLD, is able to restore the altered contour of the object, suppress greatly the streak artifacts and recover to some extend the image quality by restricting the space of search with a surface constraint. However, the evaluation was only carried out using a simulated mask that described the shape of the object obtained by thresholding a previous CT image of the sample, which is generally not available in real scenarios. The general objective of this thesis is the designing of a complete workflow to implement SCoLD in real scenarios. For that purpose, the 3D scanner Artec Eva was chosen to acquire the surface information of the sample, which was then transformed to be usable as prior information for SCoLD method. The evaluation done in a rodent study showed high similarity between the mask obtained from real data and the ideal mask obtained from a CT. Distortions in shape and streak artifacts in the limited-view FDK reconstruction were greatly reduced when using the real mask with the SCoLD reconstruction and the image quality was highly improved demonstrating the feasibility of the proposal.Grado en Ingeniería Biomédica (Plan 2010

Sinogram Restoration for Low-Dosed X-Ray Computed Tomography Using Fractional-Order Perona-Malik Diffusion

Existing integer-order Nonlinear Anisotropic Diffusion (NAD) used in noise suppressing will produce undesirable staircase effect or speckle effect. In this paper, we propose a new scheme, named Fractal-order Perona-Malik Diffusion (FPMD), which replaces the integer-order derivative of the Perona-Malik (PM) Diffusion with the fractional-order derivative using G-L fractional derivative. FPMD, which is a interpolation between integer-order Nonlinear Anisotropic Diffusion (NAD) and fourth-order partial differential equations, provides a more flexible way to balance the noise reducing and anatomical details preserving. Smoothing results for phantoms and real sinograms show that FPMD with suitable parameters can suppress the staircase effects and speckle effects efficiently. In addition, FPMD also has a good performance in visual quality and root mean square errors (RMSE)