Virtual anthropology: A comparison between the performance of conventional X-ray and MDCT in investigating the trabecular structure of long bones.

Recently, modern cross-sectional imaging techniques such as multi-detector computed tomography (MDCT) have pioneered post mortem investigations, especially in forensic medicine. Such approaches can also be used to investigate bones non-invasively for anthropological purposes. Long bones are often examined in forensic cases because they are frequently discovered and transferred to medico-legal departments for investigation. To estimate their age, the trabecular structure must be examined. This study aimed to compare the performance of MDCT with conventional X-rays to investigate the trabecular structure of long bones. Fifty-two dry bones (24 humeri and 28 femora) from anthropological collections were first examined by conventional X-ray, and then by MDCT. Trabecular structure was evaluated by seven observers (two experienced and five inexperienced in anthropology) who analyzed images obtained by radiological methods. Analyses contained the measurement of one quantitative parameter (caput diameter of humerus and femur) and staging the trabecular structure of each bone. Preciseness of each technique was indicated by describing areas of trabecular destruction and particularities of the bones, such as pathological changes. Concerning quantitative parameters, the measurements demonstrate comparable results for the MDCT and conventional X-ray techniques. In contrast, the overall inter-observer reliability of the staging was low with MDCT and conventional X-ray. Reliability increased significantly when only the results of the staging performed by the two experienced observers were compared, particularly regarding the MDCT analysis. Our results also indicate that MDCT appears to be better suited to a detailed examination of the trabecular structure. In our opinion, MDCT is an adequate tool with which to examine the trabecular structure of long bones. However, adequate methods should be developed or existing methods should be adapted to MDCT.

Bi-planar 2D-to-3D registration in Fourier domain for stereoscopic X-ray motion tracking

In this paper we present a new method to track bonemovements in stereoscopic X-ray image series of the kneejoint. The method is based on two different X-ray imagesets: a rotational series of acquisitions of the stillsubject knee that will allow the tomographicreconstruction of the three-dimensional volume (model),and a stereoscopic image series of orthogonal projectionsas the subject performs movements. Tracking the movementsof bones throughout the stereoscopic image series meansto determine, for each frame, the best pose of everymoving element (bone) previously identified in the 3Dreconstructed model. The quality of a pose is reflectedin the similarity between its simulated projections andthe actual radiographs. We use direct Fourierreconstruction to approximate the three-dimensionalvolume of the knee joint. Then, to avoid the expensivecomputation of digitally rendered radiographs (DRR) forpose recovery, we reformulate the tracking problem in theFourier domain. Under the hypothesis of parallel X-raybeams, we use the central-slice-projection theorem toreplace the heavy 2D-to-3D registration of projections inthe signal domain by efficient slice-to-volumeregistration in the Fourier domain. Focusing onrotational movements, the translation-relevant phaseinformation can be discarded and we only consider scalarFourier amplitudes. The core of our motion trackingalgorithm can be implemented as a classical frame-wiseslice-to-volume registration task. Preliminary results onboth synthetic and real images confirm the validity ofour approach.

THREE DIMENSIONAL CHARACTERIZATION OF SOIL MACROPOROSITY BY X-RAY MICROTOMOGRAPHY

Analysis of the soil pore system represents an important way of characterizing soil structure. Properties such as the shape and number of pores can be determined through soil pore evaluations. This study presents a three-dimensional (3D) characterization of the shape and number of pores of a sub-tropical soil. To do so, a second generation X-ray microtomograph equipped with a plain type detector was employed. A voltage of 120 kV and current of 80 mA was applied to the X-ray tube. The soil samples analyzed were collected at three different depths (0-10, 10-20, and 20-30 cm). The results obtained allowed qualitative (images) and quantitative (3D) analyses of the soil structure, revealing the potential of the microtomographic technique, as well as the study of differences in soil macroporosity at different depths. Macroporosity was 5.14 % in the 0-10 cm layer, 5.10 % in the 10-20 cm layer, and 6.64 % in the 20-30 cm layer. The macroporosity of unclassified pores (UN) was 0.30 % (0-10 and 10-20 cm) and 0.40 % (20-30 cm), while equant pores (EQ) had values of 0.01 % at the three depths under analysis.

X-Ray systems: radiation protection program in basic research applications

Although the radiation doses involved in basic research radiology are relatively small, the increasing number of radiological procedures makes risks becoming increasingly high. Quality control techniques in radiological practice have to ensure an adequate system of protection for people exposed to radiation. These techniques belong to a quality assurance program for X-ray machines and are designed to correct problems related to equipment and radiological practices, to obtain radiological images of high quality and to reduce the unnecessary exposures.

Study of OSEM with different subsets in grating-based X-ray differential phase-contrast imaging.

Impressive developments in X-ray imaging are associated with X-ray phase contrast computed tomography based on grating interferometry, a technique that provides increased contrast compared with conventional absorption-based imaging. A new "single-step" method capable of separating phase information from other contributions has been recently proposed. This approach not only simplifies data-acquisition procedures, but, compared with the existing phase step approach, significantly reduces the dose delivered to a sample. However, the image reconstruction procedure is more demanding than for traditional methods and new algorithms have to be developed to take advantage of the "single-step" method. In the work discussed in this paper, a fast iterative image reconstruction method named OSEM (ordered subsets expectation maximization) was applied to experimental data to evaluate its performance and range of applicability. The OSEM algorithm with different subsets was also characterized by comparison of reconstruction image quality and convergence speed. Computer simulations and experimental results confirm the reliability of this new algorithm for phase-contrast computed tomography applications. Compared with the traditional filtered back projection algorithm, in particular in the presence of a noisy acquisition, it furnishes better images at a higher spatial resolution and with lower noise. We emphasize that the method is highly compatible with future X-ray phase contrast imaging clinical applications.

Investigation of discrete imaging models and iterative image reconstruction in differential X-ray phase-contrast tomography.

Differential X-ray phase-contrast tomography (DPCT) refers to a class of promising methods for reconstructing the X-ray refractive index distribution of materials that present weak X-ray absorption contrast. The tomographic projection data in DPCT, from which an estimate of the refractive index distribution is reconstructed, correspond to one-dimensional (1D) derivatives of the two-dimensional (2D) Radon transform of the refractive index distribution. There is an important need for the development of iterative image reconstruction methods for DPCT that can yield useful images from few-view projection data, thereby mitigating the long data-acquisition times and large radiation doses associated with use of analytic reconstruction methods. In this work, we analyze the numerical and statistical properties of two classes of discrete imaging models that form the basis for iterative image reconstruction in DPCT. We also investigate the use of one of the models with a modern image reconstruction algorithm for performing few-view image reconstruction of a tissue specimen.

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.

Detection of beetle damage in forests by X-ray CT image processing

Some beetle species can have devastating economic impacts on forest and nursery industries. A recent example is Anophophora glabripennis, a species of beetle known in the United States as the ''Asian Longhorrned beetle'', which has damaged many American forests, and is a threat which can unintentionally reach south American countries, including Brazil. This work presents a new method based on X-ray computerized tomography (CT) and image processing for beetle injury detection in forests. Its results show a set of images with correct identification of the location of beetles in living trees as well as damage evaluation with time.

Application of x-ray computed tomography in the evaluation of soil porosity in soil management systems

The study aimed to evaluate a methodology to quantify the porosity of the soil using computed tomography in areas under no-tillage, conventional tillage and native forest. Three soil management systems were selected for the study: forest, conventional tillage and no-tillage. In each soil management system, undisturbed soil samples were collected in the surface layer (0.0 to 0.10 m). The tomographic images were obtained using a X-ray microtomography. After obtaining the images, they were processed, and a methodology was evaluated for image conversion into numerical values. The statistical method which provided the greatest accuracy was the percentile method. The methodology used to analyze the tomographic image allowed quantifying the porosity of the soil under different soil management. The method enabled the characterization of soil porosity in a non-evasive and non-destructive way.

Articulated Model Registration of MRI/X-Ray Spine Data

This paper presents a method based on articulated models for the registration of spine data extracted from multimodal medical images of patients with scoliosis. With the ultimate aim being the development of a complete geometrical model of the torso of a scoliotic patient, this work presents a method for the registration of vertebral column data using 3D magnetic resonance images (MRI) acquired in prone position and X-ray data acquired in standing position for five patients with scoliosis. The 3D shape of the vertebrae is estimated from both image modalities for each patient, and an articulated model is used in order to calculate intervertebral transformations required in order to align the vertebrae between both postures. Euclidean distances between anatomical landmarks are calculated in order to assess multimodal registration error. Results show a decrease in the Euclidean distance using the proposed method compared to rigid registration and more physically realistic vertebrae deformations compared to thin-plate-spline (TPS) registration thus improving alignment.

An X-ray micro-tomography system optimised for the low-dose study of living organisms

An X-ray micro-tomography system has been designed that is dedicated to the low-dose imaging of radiation sensitive living organisms and has been used to image the early development of the first few days of plant development immediately after germination. The system is based on third-generation X-ray micro-tomography system and consists of an X-ray tube, two-dimensional X-ray detector and a mechanical sample manipulation stage. The X-ray source is a 50 kVp X-ray tube with a silver target with a filter to centre the X-ray spectrum on 22 keV.A 100 mm diameter X-ray image intensifier (XRII) is used to collect the two-dimensional projection images. The rotation tomography table incorporates a linear translation mechanism to eliminate ring artefact that is commonly associated with third-generation tomography systems' Developing maize seeds (Triticum aestivum) have been imaged using the system with a cubic voxel linear dimension of 100 mum, over a diameter of 25 mm and the root lengths and volumes measured. The X-ray dose to the plants was also assessed and found to have no effect on the plant root development. (C) 2003 Elsevier Science Ltd. All rights reserved.

3D scanning by multiple fan beam X-ray sources and sensors

*** Purpose – Computer tomography (CT) for 3D reconstruction entails a huge number of coplanar fan-beam projections for each of a large number of 2D slice images, and excessive radiation intensities and dosages. For some applications its rate of throughput is also inadequate. A technique for overcoming these limitations is outlined. *** Design methodology/approach – A novel method to reconstruct 3D surface models of objects is presented, using, typically, ten, 2D projective images. These images are generated by relative motion between this set of objects and a set of ten fanbeam X-ray sources and sensors, with their viewing axes suitably distributed in 2D angular space. *** Findings – The method entails a radiation dosage several orders of magnitude lower than CT, and requires far less computational power. Experimental results are given to illustrate the capability of the technique *** Practical implications – The substantially lower cost of the method and, more particularly, its dramatically lower irradiation make it relevant to many applications precluded by current techniques *** Originality/value – The method can be used in many applications such as aircraft hold-luggage screening, 3D industrial modelling and measurement, and it should also have important applications to medical diagnosis and surgery.

Estimating root–soil contact from 3D X-ray microtomographs

Adequate contact with the soil is essential for water and nutrient adsorption by plant roots, but the determination of root–soil contact is a challenging task because it is difficult to visualize roots in situ and quantify their interactions with the soil at the scale of micrometres. A method to determine root–soil contact using X-ray microtomography was developed. Contact areas were determined from 3D volumetric images using segmentation and iso-surface determination tools. The accuracy of the method was tested with physical model systems of contact between two objects (phantoms). Volumes, surface areas and contact areas calculated from the measured phantoms were compared with those estimated from image analysis. The volume was accurate to within 0.3%, the surface area to within 2–4%, and the contact area to within 2.5%. Maize and lupin roots were grown in soil (<2 mm) and vermiculite at matric potentials of −0.03 and −1.6 MPa and in aggregate fractions of 4–2, 2–1, 1–0.5 and < 0.5 mm at a matric potential of −0.03 MPa. The contact of the roots with their growth medium was determined from 3D volumetric images. Macroporosity (>70 µm) of the soil sieved to different aggregate fractions was calculated from binarized data. Root-soil contact was greater in soil than in vermiculite and increased with decreasing aggregate or particle size. The differences in root–soil contact could not be explained solely by the decrease in porosity with decreasing aggregate size but may also result from changes in particle and aggregate packing around the root.

Grazing incidence X-ray diffraction and atomic force microscopy analysis of BaBi2Ta2O9 thin films

Thin films of BaBi2Ta2O9 (BBT) composition were prepared through the metal organic decomposition method. The crystallinity, phase formation, crystallite size and morphology of the thin films were measured as a function of the type of substrate, stoichiometry of solution and process variables such as thickness and temperature. The thin films were investigated by grazing incidence X-ray diffractometry and atomic force microscopy (AFM) techniques. For the sample without excess of bismuth, diffraction peaks other than that of the BBT phase were observed. A well crystallized BBT single phase was observed for films prepared from a solution with 10% excess of bismuth, deposited on Si/Pt substrate, with a thickness up to 150 nm and sintered at temperatures of 700 degreesC. The thin BBT phase films heat-treated at 600 degreesC presented a diffraction pattern characteristic of samples with lower degree of crystallinity whereas for the thin films heat-treated at 800 degreesC, we observed the presence of other phases than the BBT. For the thin film deposited on the Sin+ substrate, we observe that the peaks corresponding to the BBT phase are broader than that observed on the samples deposited on the Pt and Si/Pt substrates. No variation of average crystallite size was observed as the excess of Bi increased from 10 to 20%. AFM images for the samples showed that the increasing the amount of bismuth promotes grain growth. The average surface roughness measured was in the range of 16-22 nm showing that the bismuth amount had no or little effect on the roughness of films. (C) 2002 Elsevier B.V. B.V. All rights reserved.

An experimental and theoretical approach to correct for the scattered radiation in an X-ray computer tomography system for industrial applications

The main problem connected to cone beam computed tomography (CT) systems for industrial applications employing 450 kV X-ray tubes is the high amount of scattered radiation which is added to the primary radiation (signal). This stray radiation leads to a significant degradation of the image quality. A better understanding of the scattering and methods to reduce its effects are therefore necessary to improve the image quality. Several studies have been carried out in the medical field at lower energies, whereas studies in industrial CT, especially for energies up to 450 kV, are lacking. Moreover, the studies reported in literature do not consider the scattered radiation generated by the CT system structure and the walls of the X-ray room (environmental scatter). In order to investigate the scattering on CT projections a GEANT4-based Monte Carlo (MC) model was developed. The model, which has been validated against experimental data, has enabled the calculation of the scattering including the environmental scatter, the optimization of an anti-scatter grid suitable for the CT system, and the optimization of the hardware components of the CT system. The investigation of multiple scattering in the CT projections showed that its contribution is 2.3 times the one of primary radiation for certain objects. The results of the environmental scatter showed that it is the major component of the scattering for aluminum box objects of front size 70 x 70 mm2 and that it strongly depends on the thickness of the object and therefore on the projection. For that reason, its correction is one of the key factors for achieving high quality images. The anti-scatter grid optimized by means of the developed MC model was found to reduce the scatter-toprimary ratio in the reconstructed images by 20 %. The object and environmental scatter calculated by means of the simulation were used to improve the scatter correction algorithm which could be patented by Empa. The results showed that the cupping effect in the corrected image is strongly reduced. The developed CT simulation is a powerful tool to optimize the design of the CT system and to evaluate the contribution of the scattered radiation to the image. Besides, it has offered a basis for a new scatter correction approach by which it has been possible to achieve images with the same spatial resolution as state-of-the-art well collimated fan-beam CT with a gain in the reconstruction time of a factor 10. This result has a high economic impact in non-destructive testing and evaluation, and reverse engineering.