Applications of advanced and dual energy computed tomography in proton therapy

This thesis focuses on advanced reconstruction methods and Dual Energy (DE) Computed Tomography (CT) applications for proton therapy, aiming at improving patient positioning and investigating approaches to deal with metal artifacts. To tackle the first goal, an algorithm for post-processing input DE images has been developed. The outputs are tumor- and bone-canceled images, which help in recognising structures in patient body. We proved that positioning error is substantially reduced using contrast enhanced images, thus suggesting the potential of such application. If positioning plays a key role in the delivery, even more important is the quality of planning CT. For that, modern CT scanners offer possibility to tackle challenging cases, like treatment of tumors close to metal implants. Possible approaches for dealing with artifacts introduced by such rods have been investigated experimentally at Paul Scherrer Institut (Switzerland), simulating several treatment plans on an anthropomorphic phantom. In particular, we examined the cases in which none, manual or Iterative Metal Artifact Reduction (iMAR) algorithm were used to correct the artifacts, using both Filtered Back Projection and Sinogram Affirmed Iterative Reconstruction as image reconstruction techniques. Moreover, direct stopping power calculation from DE images with iMAR has also been considered as alternative approach. Delivered dose measured with Gafchromic EBT3 films was compared with the one calculated in Treatment Planning System. Residual positioning errors, daily machine dependent uncertainties and film quenching have been taken into account in the analyses. Although plans with multiple fields seemed more robust than single field, results showed in general better agreement between prescribed and delivered dose when using iMAR, especially if combined with DE approach. Thus, we proved the potential of these advanced algorithms in improving dosimetry for plans in presence of metal implants.

Computed Ultrasound Tomography in Echo mode (CUTE) of speed of sound for diagnosis and for aberration correction in pulse-echo sonography

Sound speed as a diagnostic marker for various diseases of human tissue has been of interest for a while. Up to now, mostly transmission ultrasound computed tomography (UCT) was able to detect spatially resolved sound speed, and its promise as a diagnostic tool has been demonstrated. However, UCT is limited to acoustically transparent samples such as the breast. We present a novel technique where spatially resolved detection of sound speed can be achieved using conventional pulse-echo equipment in reflection mode. For this purpose, pulse-echo images are acquired under various transmit beam directions and a two-dimensional map of the sound speed is reconstructed from the changing phase of local echoes using a direct reconstruction method. Phantom results demonstrate that a high spatial resolution (1 mm) and contrast (0.5 % of average sound speed) can be achieved suitable for diagnostic purposes. In comparison to previous reflection-mode based methods, CUTE works also in a situation with only diffuse echoes, and its direct reconstruction algorithm enables real-time application. This makes it suitable as an addition to conventional clinical ultrasound where it has the potential to benefit diagnosis in a multimodal approach. In addition, knowledge of the spatial distribution of sound speed allows full aberration correction and thus improved spatial resolution and contrast of conventional B-mode ultrasound. © (2014) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.

Pore formation during dehydration of polycrystalline gypsum observed and quantified in a time-series synchrotron radiation based X-ray micro-tomography experiment

We conducted an in-situ X-ray micro-computed tomography heating experiment at the Advanced Photon Source (USA) to dehydrate an unconfined 2.3 mm diameter cylinder of Volterra Gypsum. We used a purpose-built X-ray transparent furnace to heat the sample to 388 K for a total of 310 min to acquire a three-dimensional time-series tomography dataset comprising nine time steps. The voxel size of 2.2 μm3 proved sufficient to pinpoint reaction initiation and the organization of drainage architecture in space and time. We observed that dehydration commences across a narrow front, which propagates from the margins to the centre of the sample in more than four hours. The advance of this front can be fitted with a square-root function, implying that the initiation of the reaction in the sample can be described as a diffusion process. Novel parallelized computer codes allow quantifying the geometry of the porosity and the drainage architecture from the very large tomographic datasets (20483 voxels) in unprecedented detail. We determined position, volume, shape and orientation of each resolvable pore and tracked these properties over the duration of the experiment. We found that the pore-size distribution follows a power law. Pores tend to be anisotropic but rarely crack-shaped and have a preferred orientation, likely controlled by a pre-existing fabric in the sample. With on-going dehydration, pores coalesce into a single interconnected pore cluster that is connected to the surface of the sample cylinder and provides an effective drainage pathway. Our observations can be summarized in a model in which gypsum is stabilized by thermal expansion stresses and locally increased pore fluid pressures until the dehydration front approaches to within about 100 μm. Then, the internal stresses are released and dehydration happens efficiently, resulting in new pore space. Pressure release, the production of pores and the advance of the front are coupled in a feedback loop.

Morphometric analysis of the thoracic intervertebral foramen osseous anatomy in adolescent idiopathic scoliosis using low dose computerized tomography

Introduction. The dimensions of the thoracic intervertebral foramen in adolescent idiopathic scoliosis (AIS) have not previously been quantified. During posterior approach scoliosis correction surgery pedicle screws may occasionally breach into the foramen. Better understanding of the dimensions of the foramen may be useful in surgical planning. This study describes a reproducible method for measurement of the thoracic foramen in AIS using computerized tomography (CT). Methods. In 23 pre-operative female patients with Lenke 1 type AIS with right side convexity major curves confined to the thoracic spine the foraminal height (FH), foraminal width (FW), pedicle to superior articular process distance (P-SAP) and cross sectional foraminal area (FA) were measured using multiplanar reconstructed CT. Measurements were made at entrance, midpoint and exit of the thoracic foramina from T1/T2 to T11/T12. Results were correlated with potential dependent variables of major curve Cobb Angle measured on X-ray and CT, Age, Weight, Lenke classification subtype, Risser Grade and number of spinal levels in the major curve. Results. The FH, FW, P-SAP and FA dimensions and ratios are all significantly larger on the convexity of the major curve and maximal at or close to the apex. Mean thoracic foraminal dimensions change in a predictable manner relative to position on the major thoracic curve. There was no significant correlation with the measured foraminal dimensions or ratios and the potential dependent variables. The average ratio of convexity to concavity dimensions at the apex foramina for entrance, midpoint and exit respectively are FH (1.50, 1.38, 1.25), FW (1.28, 1.30, 0.98), FA (2.06, 1.84, 1.32), P-SAP (1.61, 1.47, 1.30). Conclusion. Foraminal dimensions of the thoracic spine are significantly affected by AIS. Foraminal dimensions have a predictable convexity to concavity ratio relative to the proximity to the major curve apex. Surgeons should be aware of these anatomical differences during scoliosis correction surgery.

Using low dose computerized tomography to morphometrically analyse the thoracic intervertebral foramen osseous anatomy in adolescent idiopathic scoliosis

INTRODUCTION The dimensions of the thoracic intervertebral foramen in adolescent idiopathic scoliosis (AIS) have not previously been quantified. During posterior approach scoliosis correction surgery pedicle screws may occasionally breach into the foramen. Better understanding of the dimensions of the foramen may be useful in surgical planning. This study describes a reproducible method for measurement of the thoracic foramen in AIS using computerized tomography (CT). METHODS In 23 pre-operative female patients with Lenke 1 type AIS with right side convexity major curves confined to the thoracic spine the foraminal height (FH), foraminal width (FW), pedicle to superior articular process distance (P-SAP) and cross sectional foraminal area (FA) were measured using multiplanar reconstructed CT. Measurements were made at entrance, midpoint and exit of the thoracic foramina from T1/T2 to T11/T12. Results were correlated with potential dependent variables of major curve Cobb Angle measured on X-ray and CT, Age, Weight, Lenke classification subtype, Risser Grade and number of spinal levels in the major curve. RESULTS The FH, FW, P-SAP and FA dimensions and ratios are all significantly larger on the convexity of the major curve and maximal at or close to the apex. Mean thoracic foraminal dimensions change in a predictable manner relative to position on the major thoracic curve. There was no significant correlation with the measured foraminal dimensions or ratios and the potential dependent variables. The average ratio of convexity to concavity dimensions at the apex foramina for entrance, midpoint and exit respectively are FH (1.50, 1.38, 1.25), FW (1.28, 1.30, 0.98), FA (2.06, 1.84, 1.32), P-SAP (1.61, 1.47, 1.30). CONCLUSION Foraminal dimensions of the thoracic spine are significantly affected by AIS. Foraminal dimensions have a predictable convexity to concavity ratio relative to the proximity to the major curve apex. Surgeons should be aware of these anatomical differences during scoliosis correction surgery.

A MatLAB Based Virtual Phantom for 2D Electrical Impedance Tomography (MatVP2DEIT): Studying the Medical Electrical Impedance Tomography Reconstruction in Computer

Electrical Impedance Tomography (EIT) is a computerized medical imaging technique which reconstructs the electrical impedance images of a domain under test from the boundary voltage-current data measured by an EIT electronic instrumentation using an image reconstruction algorithm. Being a computed tomography technique, EIT injects a constant current to the patient's body through the surface electrodes surrounding the domain to be imaged (Omega) and tries to calculate the spatial distribution of electrical conductivity or resistivity of the closed conducting domain using the potentials developed at the domain boundary (partial derivative Omega). Practical phantoms are essentially required to study, test and calibrate a medical EIT system for certifying the system before applying it on patients for diagnostic imaging. Therefore, the EIT phantoms are essentially required to generate boundary data for studying and assessing the instrumentation and inverse solvers a in EIT. For proper assessment of an inverse solver of a 2D EIT system, a perfect 2D practical phantom is required. As the practical phantoms are the assemblies of the objects with 3D geometries, the developing of a practical 2D-phantom is a great challenge and therefore, the boundary data generated from the practical phantoms with 3D geometry are found inappropriate for assessing a 2D inverse solver. Furthermore, the boundary data errors contributed by the instrumentation are also difficult to separate from the errors developed by the 3D phantoms. Hence, the errorless boundary data are found essential to assess the inverse solver in 2D EIT. In this direction, a MatLAB-based Virtual Phantom for 2D EIT (MatVP2DEIT) is developed to generate accurate boundary data for assessing the 2D-EIT inverse solvers and the image reconstruction accuracy. MatVP2DEIT is a MatLAB-based computer program which simulates a phantom in computer and generates the boundary potential data as the outputs by using the combinations of different phantom parameters as the inputs to the program. Phantom diameter, inhomogeneity geometry (shape, size and position), number of inhomogeneities, applied current magnitude, background resistivity, inhomogeneity resistivity all are set as the phantom variables which are provided as the input parameters to the MatVP2DEIT for simulating different phantom configurations. A constant current injection is simulated at the phantom boundary with different current injection protocols and boundary potential data are calculated. Boundary data sets are generated with different phantom configurations obtained with the different combinations of the phantom variables and the resistivity images are reconstructed using EIDORS. Boundary data of the virtual phantoms, containing inhomogeneities with complex geometries, are also generated for different current injection patterns using MatVP2DEIT and the resistivity imaging is studied. The effect of regularization method on the image reconstruction is also studied with the data generated by MatVP2DEIT. Resistivity images are evaluated by studying the resistivity parameters and contrast parameters estimated from the elemental resistivity profiles of the reconstructed phantom domain. Results show that the MatVP2DEIT generates accurate boundary data for different types of single or multiple objects which are efficient and accurate enough to reconstruct the resistivity images in EIDORS. The spatial resolution studies show that, the resistivity imaging conducted with the boundary data generated by MatVP2DEIT with 2048 elements, can reconstruct two circular inhomogeneities placed with a minimum distance (boundary to boundary) of 2 mm. It is also observed that, in MatVP2DEIT with 2048 elements, the boundary data generated for a phantom with a circular inhomogeneity of a diameter less than 7% of that of the phantom domain can produce resistivity images in EIDORS with a 1968 element mesh. Results also show that the MatVP2DEIT accurately generates the boundary data for neighbouring, opposite reference and trigonometric current patterns which are very suitable for resistivity reconstruction studies. MatVP2DEIT generated data are also found suitable for studying the effect of the different regularization methods on reconstruction process. Comparing the reconstructed image with an original geometry made in MatVP2DEIT, it would be easier to study the resistivity imaging procedures as well as the inverse solver performance. Using the proposed MatVP2DEIT software with modified domains, the cross sectional anatomy of a number of body parts can be simulated in PC and the impedance image reconstruction of human anatomy can be studied.

Extraction of fibre network architecture by X-ray tomography and prediction of elastic properties using an affine analytical model

This paper proposes a method for extracting reliable architectural characteristics from complex porous structures using micro-computed tomography (μCT) images. The work focuses on a highly porous material composed of a network of fibres bonded together. The segmentation process, allowing separation of the fibres from the remainder of the image, is the most critical step in constructing an accurate representation of the network architecture. Segmentation methods, based on local and global thresholding, were investigated and evaluated by a quantitative comparison of the architectural parameters they yielded, such as the fibre orientation and segment length (sections between joints) distributions and the number of inter-fibre crossings. To improve segmentation accuracy, a deconvolution algorithm was proposed to restore the original images. The efficacy of the proposed method was verified by comparing μCT network architectural characteristics with those obtained using high resolution CT scans (nanoCT). The results indicate that this approach resolves the architecture of these complex networks and produces results approaching the quality of nanoCT scans. The extracted architectural parameters were used in conjunction with an affine analytical model to predict the axial and transverse stiffnesses of the fibre network. Transverse stiffness predictions were compared with experimentally measured values obtained by vibration testing. © 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

The geometry and the number of contacts of monodisperse sphere packs using X-ray tomography

The three-dimensional structure of very large samples of monodisperse bead packs is studied by means of X-Ray Computed Tomography. We retrieve the coordinatesofeach bead inthe pack and wecalculate the average coordination number by using the tomographic images to single out the neighbors in contact. The results are compared with the average coordination number obtained in Aste et al. (2005) by using a deconvolution technique. We show that the coordination number increases with the packing fraction, varying between 6.9 and 8.2 for packing fractions between 0.59 and 0.64. © 2005 Taylor & Francis Group.

Microcomputed tomography imaging in a rat model of delayed union/non-union fracture

We aimed to develop a clinically relevant delayed union/non-union fracture model to evaluate a cell therapy intervention repair strategy. Histology, three-dimensional (3D) micro-computed tomography (micro-CT) imaging and mechanical testing were utilized to develop an analytical protocol for qualitative and quantitative assessment of fracture repair. An open femoral diaphyseal osteotomy, combined with periosteal diathermy and endosteal excision, was held in compression by a four pin unilateral external fixator. Three delayed union/non-union fracture groups established at 6 weeks-(a) a control group, (b) a cell therapy group, and (c) a group receiving phosphate-buffered saline (PBS) injection alone-were examined subsequently at 8 and 14 weeks. The histological response was combined fibrous and cartilaginous non-unions in groups A and B with fibrous non-unions in group C. Mineralized callus volume/total volume percentage showed no statistically significant differences between groups. Endosteal calcified tissue volume/endosteal tissue volume, at the center of the fracture site, displayed statistically significant differences between 8 and 14 weeks for cell and PBS intervention groups but not for the control group. The percentage load to failure was significantly lower in the control and cell treatment groups than in the PBS alone group. High-resolution micro-CT imaging provides a powerful tool to augment characterization of repair in delayed union/non-union fractures together with outcomes such as histology and mechanical strength measurement. Accurate, nondestructive, 3D identification of mineralization progression in repairing fractures is enabled in the presence or absence of intervention strategies. (c) 2007 Orthopaedic Research Society.

Verification of computed tomographic estimates of cochlear implant array position: A micro-CT and histological analysis

We examine the efficacy two volume spatial registration of pre and postoperative clinical computed tomography (CT) imaging to verify post-operative electrode array placement in cochlear implant (CI) patients. To measure the degree of accuracy with which the composite image predicts in-vivo placement of the array, we replicate the CI surgical process in cadaver heads. Pre-operative, post-operative, micro CT imaging and histology are utilized for verification.

Three dimensional analysis of cardiovascular development in mouse embryos using X-ray microcomputed tomography

Micro-computed tomography (μCT) has been successfully used to study the cardiovascular system of mouse embryos in situ. With the use of barium as a suitable contrast agent, blood vessels have been imaged and analysed quantitatively such as blood volume and vessel sizes on embryos of ages 14.5 to 16.5 days old. The advantage of using this imaging modality is that it has provided three dimensional information whilst leaving samples intact for further study.

Computed tomographic-dacryocystography (CT-DCG) of the normal canine nasolacrimal drainage system with three-dimensional reconstruction

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Synchrotron microComputed tomography of the mature bovine dentinoenamel junction

The mature dentinoenamel junction (DEJ) is viewed by some investigators and the current authors, not as a fossilized, sharp transition between enamel and dentin, but as a relatively broad structural transition zone including the mantle dentin and the inner aprismatic enamel. In this study, the DEJ structure in bovine incisors was studied with synchrotron microComputed Tomography (microCT) using small cubes cut parallel to the tooth surface. The reconstructions revealed a zone of highly variable punctate contrast between bulk dentin and enamel; the mean linear attenuation coefficients and their standard deviations demonstrated that this zone averaged less mineral than dentin or enamel but had more highly variable structure than either. The region with the punctuate contrast is, therefore, the mantle dentin. The thickness of the mantle dentin seen in a typical data set was about 30 mu m, and the mantle dentin-enamel interface deviated +/- 15 mu m from the average plane over a distance of 520 mu m. In the highest resolution data (similar to 1.5 mu m isotropic voxels, volume elements), tubules in the dentin could be discerned in the vicinity of the DEJ. Contrast sensitivity was high enough to detect differences in mineral content between near-surface and near-DEJ volumes of the enamel. Reconstructions before and after two cubes were compressed to failure revealed cracks formed only in the enamel and did not propagate across the mantle dentin, regardless of whether loading was parallel to or perpendicular to the DEJ. (C) 2007 Elsevier B.V. All rights reserved.

Development and applications of three-dimensional gamma ray tomography system using ray casting volume rendering techniques

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Effects of contrast material on computed tomographic measurements of lung volumes in patients with acute lung injury

Background. Intravenous injection of contrast material is routinely performed in order to differentiate nonaerated lung parenchyma from pleural effusion in critically ill patients undergoing thoracic computed tomography (CT). The aim of the present study was to evaluate the effects of contrast material on CT measurement of lung volumes in 14 patients with acute lung injury. Method. A spiral thoracic CT scan, consisting of contiguous axial sections of 10 mm thickness, was performed from the apex to the diaphragm at end-expiration both before and 30 s (group 1; n=7) or 15 min (group 2; n=7) after injection of 80 ml contrast material. Volumes of gas and tissue, and volumic distribution of CT attenuations were measured before and after injection using specially designed software (Lungview®; Institut National des Télécommunications, Evry, France). The maximal artifactual increase in lung tissue resulting from a hypothetical leakage within the lung of the 80 ml contrast material was calculated. Results. Injection of contrast material significantly increased the apparent volume of lung tissue by 83 ± 57 ml in group 1 and 102 ± 80 ml in group 2, whereas the corresponding maximal artifactual increases in lung tissue were 42 ± 52 ml and 31 ± 18 ml. Conclusion. Because systematic injection of contrast material increases the amount of extravascular lung water in patients with acute lung injury, it seems prudent to avoid this procedure in critically ill patients undergoing a thoracic CT scan and to reserve its use for specific indications.