Statistical Inversion Theory in X-ray Tomography

In this thesis the X-ray tomography is discussed from the Bayesian statistical viewpoint. The unknown parameters are assumed random variables and as opposite to traditional methods the solution is obtained as a large sample of the distribution of all possible solutions. As an introduction to tomography an inversion formula for Radon transform is presented on a plane. The vastly used filtered backprojection algorithm is derived. The traditional regularization methods are presented sufficiently to ground the Bayesian approach. The measurements are foton counts at the detector pixels. Thus the assumption of a Poisson distributed measurement error is justified. Often the error is assumed Gaussian, altough the electronic noise caused by the measurement device can change the error structure. The assumption of Gaussian measurement error is discussed. In the thesis the use of different prior distributions in X-ray tomography is discussed. Especially in severely ill-posed problems the use of a suitable prior is the main part of the whole solution process. In the empirical part the presented prior distributions are tested using simulated measurements. The effect of different prior distributions produce are shown in the empirical part of the thesis. The use of prior is shown obligatory in case of severely ill-posed problem.

Gamma spectrometry and gamma and X-ray tomography of nuclear fuel

The purpose of gamma spectrometry and gamma and X-ray tomography of nuclear fuel is to determine both radionuclide concentration and integrity and deformation of nuclear fuel. The aims of this thesis have been to find out the basics of gamma spectrometry and tomography of nuclear fuel, to find out the operational mechanisms of gamma spectrometry and tomography equipment of nuclear fuel, and to identify problems that relate to these measurement techniques. In gamma spectrometry of nuclear fuel the gamma-ray flux emitted from unstable isotopes is measured using high-resolution gamma-ray spectroscopy. The production of unstable isotopes correlates with various physical fuel parameters. In gamma emission tomography the gamma-ray spectrum of irradiated nuclear fuel is recorded for several projections. In X-ray transmission tomography of nuclear fuel a radiation source emits a beam and the intensity, attenuated by the nuclear fuel, is registered by the detectors placed opposite. When gamma emission or X-ray transmission measurements are combined with tomographic image reconstruction methods, it is possible to create sectional images of the interior of nuclear fuel. MODHERATO is a computer code that simulates the operation of radioscopic or tomographic devices and it is used to predict and optimise the performance of imaging systems. Related to the X-ray tomography, MODHERATO simulations have been performed by the author. Gamma spectrometry and gamma and X-ray tomography are promising non-destructive examination methods for understanding fuel behaviour under normal, transient and accident conditions.

Construction and test of low cost X-ray tomography scanner for physical-chemical analysis and nondestructive inspections

X-ray computed tomography (CT) refers to the cross-sectional imaging of an object measuring the transmitted radiation at different directions. In this work, we describe the development of a low cost micro-CT X-ray scanner that is being developed for nondestructive testing. This tomograph operates using a microfocus X-ray source and contains a silicon photodiode as detectors. The performance of the system, by its spatial resolution, has been estimated through its Modulation Transfer Function - MTF and the obtained value at 10% of MTF is 661 μm. It was built as a general purpose nondestructive testing device. © 2009 American Institute of Physics.

Finite element 3D reconstruction of the pulmonary acinus imaged by synchrotron X-ray tomography

The alveolated structure of the pulmonary acinus plays a vital role in gas exchange function. Three-dimensional (3D) analysis of the parenchymal region is fundamental to understanding this structure-function relationship, but only a limited number of attempts have been conducted in the past because of technical limitations. In this study, we developed a new image processing methodology based on finite element (FE) analysis for accurate 3D structural reconstruction of the gas exchange regions of the lung. Stereologically well characterized rat lung samples (Pediatr Res 53: 72-80, 2003) were imaged using high-resolution synchrotron radiation-based X-ray tomographic microscopy. A stack of 1,024 images (each slice: 1024 x 1024 pixels) with resolution of 1.4 mum(3) per voxel were generated. For the development of FE algorithm, regions of interest (ROI), containing approximately 7.5 million voxels, were further extracted as a working subunit. 3D FEs were created overlaying the voxel map using a grid-based hexahedral algorithm. A proper threshold value for appropriate segmentation was iteratively determined to match the calculated volume density of tissue to the stereologically determined value (Pediatr Res 53: 72-80, 2003). The resulting 3D FEs are ready to be used for 3D structural analysis as well as for subsequent FE computational analyses like fluid dynamics and skeletonization.

Experimental Characterization and Quantification of Cement-Bentonite Interaction using Core Infiltration Techniques Coupled with X-ray Tomography

Deep geological storage of radioactive waste foresees cementitious materials as reinforcement of tunnels and as backfill. Bentonite is proposed to enclose spent fuel canisters and as drift seals. Sand/bentonite (s/b) is foreseen as backfill material of access galleries or as drift seals. The emplacement of cementitious material next to clay material generates an enormous chemical gradient in pore-water composition that drives diffusive solute transport. Laboratory studies and reactive transport modeling predicted significant mineral alteration at and near interfaces, mainly resulting in a decrease of porosity in bentonite. The goal of this thesis was to characterize and quantify the cement/bentonite interactions both spatially and temporally in laboratory experiments. A newly developed mobile X-ray transparent core infiltration device was used to perform X-ray computed tomography (CT) scans without interruption of running experiments. CT scans allowed tracking the evolution of the reaction plume and changes in core volume/diameter/density during the experiments. In total 4 core infiltration experiments were carried out for this study with the compacted and saturated cores consisting of MX-80 bentonite and sand/MX-80 bentonite mixture (s/b; 65/35%). Two different high-pH cementitious pore-fluids were infiltrated: a young (early) ordinary Portland cement pore-fluid (APWOPC; K+–Na+–OH-; pH 13.4; ionic strength 0.28 mol/kg) and a young ‘low-pH’ ESDRED shotcrete pore-fluid (APWESDRED; Ca2+–Na+–K+–formate; pH 11.4; ionic strength 0.11 mol/kg). The experiments lasted between 1 and 2 years. In both bentonite experiments, the hydraulic conductivity was strongly reduced after switching to high-pH fluids, changing eventually from an advective to a diffusion-dominated transport regime. The reduction was mainly induced by mineral precipitation and possibly partly also by high ionic strength pore-fluids. Both bentonite cores showed a volume reduction and a resulting transient flow in which pore-water was squeezed out during high-pH infiltration. The outflow chemistry was characterized by a high ionic strength, while chloride in the initial pore water got replaced as main anionic charge carrier by sulfate, originating from gypsum dissolution. The chemistry of the high-pH fluids got strongly buffered by the bentonite, consuming hydroxide and in case of APWESDRED also formate. Hydroxide got consumed by mineral reactions (saponite and possibly talc and brucite precipitation), while formate being affected by bacterial degradation. Post-mortem analysis showed reaction zones near the inlet of the bentonite core, characterized by calcium and magnesium enrichment, consisting predominately of calcite and saponite, respectively. Silica got enriched in the outflow, indicating dissolution of silicate-minerals, identified as preferentially cristobalite. In s/b, infiltration of APWOPC reduced the hydraulic conductivity strongly, while APWESDRED infiltration had no effect. The reduction was mainly induced by mineral precipitation and probably partly also by high ionic strength pore-fluids. Not clear is why the observed mineral precipitates in the APWESDRED experiment had no effect on the fluid flow. Both s/b cores showed a volume expansion along with decreasing ionic strengths of the outflow, due to mineral reactions or in case of APWESDRED infiltration also mediated by microbiological activity, consuming hydroxide and formate, respectively. The chemistry of the high-pH fluids got strongly buffered by the s/b. In the case of APWESDRED infiltration, formate reached the outflow only for a short time, followed by enrichment in acetate, indicating most likely biological activity. This was in agreement to post-mortem analysis of the core, observing black spots on the inflow surface, while the sample had a rotten-egg smell indicative of some sulfate reduction. Post-mortem analysis showed further in both cores a Ca-enrichment in the first 10 mm of the core due to calcite precipitation. Mg-enrichment was only observed in the APWOPC experiment, originating from newly formed saponite. Silica got enriched in the outflow of both experiments, indicating dissolution of silicate-minerals, identified in the OPC experiment as cristobalite. The experiments attested an effective buffering capacity for bentonite and s/b, a progressing coupled hydraulic-chemical sealing process and also the preservation of the physical integrity of the interface region in this setup with a total pressure boundary condition on the core sample. No complete pore-clogging was observed but the hydraulic conductivity got rather strongly reduced in 3 experiments, explained by clogging of the intergranular porosity (macroporosity). Such a drop in hydraulic conductivity may impact the saturation time of the buffer in a nuclear waste repository, although the processes and geometry will be more complex in repository situation.

Detailed damage mechanisms assessment in composite materials by means of X-ray tomography

This research focused on the evaluation of damage formation on ±45º carbon fiber laminates subjected to tensile tests. The damage was evaluated by means of X-ray tomography. A high density of cracks developed during the plateau of the stress-strain curve and were qualitatively analyzed, showing that the inner plies eventually developed a higher crack concentration than the outer plies. Delamination started to occur in the outermost ply interface when the slope after the plateau of the stress-strain curve began to increase.

Coding Strategies for X-ray Tomography

This work focuses on the construction and application of coded apertures to compressive X-ray tomography. Coded apertures can be made in a number of ways, each method having an impact on system background and signal contrast. Methods of constructing coded apertures for structuring X-ray illumination and scatter are compared and analyzed. Apertures can create structured X-ray bundles that investigate specific sets of object voxels. The tailored bundles of rays form a code (or pattern) and are later estimated through computational inversion. Structured illumination can be used to subsample object voxels and make inversion feasible for low dose computed tomography (CT) systems, or it can be used to reduce background in limited angle CT systems.

On the detection side, coded apertures modulate X-ray scatter signals to determine the position and radiance of scatter points. By forming object dependent projections in measurement space, coded apertures multiplex modulated scatter signals onto a detector. The multiplexed signals can be inverted with knowledge of the code pattern and system geometry. This work shows two systems capable of determining object position and type in a 2D plane, by illuminating objects with an X-ray `fan beam,' using coded apertures and compressive measurements. Scatter tomography can help identify materials in security and medicine that may be ambiguous with transmission tomography alone.

Are Human Peripheral Nerves Sensitive To X-ray Imaging?

Diagnostic imaging techniques play an important role in assessing the exact location, cause, and extent of a nerve lesion, thus allowing clinicians to diagnose and manage more effectively a variety of pathological conditions, such as entrapment syndromes, traumatic injuries, and space-occupying lesions. Ultrasound and nuclear magnetic resonance imaging are becoming useful methods for this purpose, but they still lack spatial resolution. In this regard, recent phase contrast x-ray imaging experiments of peripheral nerve allowed the visualization of each nerve fiber surrounded by its myelin sheath as clearly as optical microscopy. In the present study, we attempted to produce high-resolution x-ray phase contrast images of a human sciatic nerve by using synchrotron radiation propagation-based imaging. The images showed high contrast and high spatial resolution, allowing clear identification of each fascicle structure and surrounding connective tissue. The outstanding result is the detection of such structures by phase contrast x-ray tomography of a thick human sciatic nerve section. This may further enable the identification of diverse pathological patterns, such as Wallerian degeneration, hypertrophic neuropathy, inflammatory infiltration, leprosy neuropathy and amyloid deposits. To the best of our knowledge, this is the first successful phase contrast x-ray imaging experiment of a human peripheral nerve sample. Our long-term goal is to develop peripheral nerve imaging methods that could supersede biopsy procedures.

Teos-based consolidants for stone conservation: evaluation of the effects and testing by x-ray imaging

Natural stones have been widely used in the construction field since antiquity. Building materials undergo decay processes due to mechanical,chemical, physical and biological causes that can act together. Therefore an interdisciplinary approach is required in order to understand the interaction between the stone and the surrounding environment. Utilization of buildings, inadequate restoration activities and in general anthropogenic weathering factors may contribute to this degradation process. For this reasons, in the last few decades new technologies and techniques have been developed and introduced in the restoration field. Consolidants are largely used in restoration and conservation of cultural heritage in order to improve the internal cohesion and to reduce the weathering rate of building materials. It is important to define the penetration depth of a consolidant for determining its efficacy. Impregnation mainly depends on the microstructure of the stone (i.e. porosity) and on the properties of the product itself. Throughout this study, tetraethoxysilane (TEOS) applied on globigerina limestone samples has been chosen as object of investigation. After hydrolysis and condensation, TEOS deposits silica gel inside the pores, improving the cohesion of the grains. X-ray computed tomography has been used to characterize the internal structure of the limestone samples,treated and untreated with a TEOS-based consolidant. The aim of this work is to investigate the penetration depth and the distribution of the TEOS inside the porosity, using both traditional approaches and advanced X-ray tomographic techniques, the latter allowing the internal visualization in three dimensions of the materials. Fluid transport properties and porosity have been studied both at macroscopic scale, by means of capillary uptake tests and radiography, and at microscopic scale,investigated with X-ray Tomographic Microscopy (XTM). This allows identifying changes in the porosity, by comparison of the images before and after the treatment, and locating the consolidant inside the stone. Tests were initially run at University of Bologna, where characterization of the stone was carried out. Then the research continued in Switzerland: X-ray tomography and radiography were performed at Empa, Swiss Federal Laboratories for Materials Science and Technology, while XTM measurements with synchrotron radiation were run at Paul Scherrer Institute in Villigen.

20 movies of x-ray tomograms of different Alpinia fossils

The fossil floras described by Dieter MAI and Harald WALTHER are invaluable for understanding the past plant diversity in Europe, and provide important information on the occurrence of taxa in the fossil record that is critical for evolutionary studies. Among the taxa they recognized were seeds assigned to the extant genus Alpinia ROXB. (Zingiberaceae, Zingiberales). We reinvestigated 28 specimens that were assigned to Alpinia arnensis (CHANDLER) MAI, Alpinia cf. arnensis, and Alpinia bivascularis MAI from the Ypresian (lower Eocene) of the UK, upper Eocene of Germany, and lower Miocene of Germany using non-destructive synchrotron-based X-ray tomography to reveal internal anatomy. None of the samples studied show an anatomy consistent with extant Alpinia or even Zingiberales. The fossils lack the globose shape, often striate external surface, seed coat structure, operculum, and micropylar collar seen in all Alpinia, and lack the chalazal chamber seen in many Alpinia species. Two specimens from the lower Miocene of Germany showed the structure of fruits of Caricoidea CHANDLER (Cyperaceae) with a single-layered exocarp, thick mesocarp, and sclerified endocarp. The other specimens are recognized as Carpolithes albolutum nom. nov. (incertae sedis) from the Ypresian of the UK, C. phoenixnordensis sp. nov. (incertae sedis) from the upper Eocene of Germany, C. bivascularis comb. nov. (incertae sedis) from the lower Miocene of Germany as well as indeterminate tegmens from the lower Miocene of Germany. This reinvestigation demonstrates that there is, as yet, no confirmed fossil record for the extant genus Alpinia. Furthermore, at least four different taxa are recognized from what had been two extinct species, enhancing our understanding of these important European Cenozoic carpofloras.

X-ray computed tomography to quantify tree rooting spatial distributions

Poor root development due to constraining soil conditions could be an important factor influencing health of urban trees. Therefore, there is a need for efficient techniques to analyze the spatial distribution of tree roots. An analytical procedure for describing tree rooting patterns from X-ray computed tomography (CT) data is described and illustrated. Large irregularly shaped specimens of undisturbed sandy soil were sampled from Various positions around the base of trees using field impregnation with epoxy resin, to stabilize the cohesionless soil. Cores approximately 200 mm in diameter by 500 mm in height were extracted from these specimens. These large core samples were scanned with a medical X-ray CT device, and contiguous images of soil slices (2 mm thick) were thus produced. X-ray CT images are regarded as regularly-spaced sections through the soil although they are not actual 2D sections but matrices of voxels similar to 0.5 mm x 0.5 mm x 2 mm. The images were used to generate the equivalent of horizontal root contact maps from which three-dimensional objects, assumed to be roots, were reconstructed. The resulting connected objects were used to derive indices of the spatial organization of roots, namely: root length distribution, root length density, root growth angle distribution, root spatial distribution, and branching intensity. The successive steps of the method, from sampling to generation of indices of tree root organization, are illustrated through a case study examining rooting patterns of valuable urban trees. (C) 1999 Elsevier Science B.V. All rights reserved.

3D reconstruction and quantification of macropores using X-ray computed tomography and image analysis

Axial X-ray Computed tomography (CT) scanning provides a convenient means of recording the three-dimensional form of soil structure. The technique has been used for nearly two decades, but initial development has concentrated on qualitative description of images. More recently, increasing effort has been put into quantifying the geometry and topology of macropores likely to contribute to preferential now in soils. Here we describe a novel technique for tracing connected macropores in the CT scans. After object extraction, three-dimensional mathematical morphological filters are applied to quantify the reconstructed structure. These filters consist of sequences of so-called erosions and/or dilations of a 32-face structuring element to describe object distances and volumes of influence. The tracing and quantification methodologies were tested on a set of undisturbed soil cores collected in a Swiss pre-alpine meadow, where a new earthworm species (Aporrectodea nocturna) was accidentally introduced. Given the reduced number of samples analysed in this study, the results presented only illustrate the potential of the method to reconstruct and quantify macropores. Our results suggest that the introduction of the new species induced very limited chance to the soil structured for example, no difference in total macropore length or mean diameter was observed. However. in the zone colonised by, the new species. individual macropores tended to have a longer average length. be more vertical and be further apart at some depth. Overall, the approach proved well suited to the analysis of the three-dimensional architecture of macropores. It provides a framework for the analysis of complex structures, which are less satisfactorily observed and described using 2D imaging. (C) 2002 Elsevier Science B.V. All rights reserved.

CLASSIFYING RADIUS FRACTURES WITH X-RAY AND TOMOGRAPHY IMAGING

Introduction: This study evaluated the interobserver reliability of plain radiograpy versus computed tomography (CT) for the Universal and AO classification systems for distal radius fractures. Patients and methods: Five observers classified 21 sets of distal radius fractures using plain radiographs and CT independently. Kappa statistics were used to establish a relative level of agreement between observers for both readings. Results: Interobserver agreement was rated as moderate for the Universal classification and poor for the AO classification. Reducing the AO system to 9 categories and to its three main types reliability was raised to a ""moderate"" level. No difference was found for interobserver reliability between the Universal classification using plain radiographs and the Universal classification using computed tomography. Interobserver reliability of the AO classification system using plain radiographs was significantly higher than the interobserver reliability of the AO classification system using only computed tomography. Conclusion: From these data, we conclude that classification of distal radius fractures using CT scanning without plain radiographs is not beneficial.