Probabilistic Evaluation of Three-Dimensional Reconstructions from X-Ray Images Spanning a Limited Angle

An important part of computed tomography is the calculation of a three-dimensional reconstruction of an object from series of X-ray images. Unfortunately, some applications do not provide sufficient X-ray images. Then, the reconstructed objects no longer truly represent the original. Inside of the volumes, the accuracy seems to vary unpredictably. In this paper, we introduce a novel method to evaluate any reconstruction, voxel by voxel. The evaluation is based on a sophisticated probabilistic handling of the measured X-rays, as well as the inclusion of a priori knowledge about the materials that the object receiving the X-ray examination consists of. For each voxel, the proposed method outputs a numerical value that represents the probability of existence of a predefined material at the position of the voxel while doing X-ray. Such a probabilistic quality measure was lacking so far. In our experiment, false reconstructed areas get detected by their low probability. In exact reconstructed areas, a high probability predominates. Receiver Operating Characteristics not only confirm the reliability of our quality measure but also demonstrate that existing methods are less suitable for evaluating a reconstruction.

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.

Probabilistic electrical resistivity tomography of a CO2 sequestration analog

Electrical resistivity tomography (ERT) is a well-established method for geophysical characterization and has shown potential for monitoring geologic CO2 sequestration, due to its sensitivity to electrical resistivity contrasts generated by liquid/gas saturation variability. In contrast to deterministic inversion approaches, probabilistic inversion provides the full posterior probability density function of the saturation field and accounts for the uncertainties inherent in the petrophysical parameters relating the resistivity to saturation. In this study, the data are from benchtop ERT experiments conducted during gas injection into a quasi-2D brine-saturated sand chamber with a packing that mimics a simple anticlinal geological reservoir. The saturation fields are estimated by Markov chain Monte Carlo inversion of the measured data and compared to independent saturation measurements from light transmission through the chamber. Different model parameterizations are evaluated in terms of the recovered saturation and petrophysical parameter values. The saturation field is parameterized (1) in Cartesian coordinates, (2) by means of its discrete cosine transform coefficients, and (3) by fixed saturation values in structural elements whose shape and location is assumed known or represented by an arbitrary Gaussian Bell structure. Results show that the estimated saturation fields are in overall agreement with saturations measured by light transmission, but differ strongly in terms of parameter estimates, parameter uncertainties and computational intensity. Discretization in the frequency domain (as in the discrete cosine transform parameterization) provides more accurate models at a lower computational cost compared to spatially discretized (Cartesian) models. A priori knowledge about the expected geologic structures allows for non-discretized model descriptions with markedly reduced degrees of freedom. Constraining the solutions to the known injected gas volume improved estimates of saturation and parameter values of the petrophysical relationship. (C) 2014 Elsevier B.V. All rights reserved.

Extended Cahill-Glauber formalism for finite-dimensional spaces. II. Applications in quantum tomography and quantum teleportation

By means of a mod(N)-invariant operator basis, s-parametrized phase-space functions associated with bounded operators in a finite-dimensional Hilbert space are introduced in the context of the extended Cahill-Glauber formalism, and their properties are discussed in details. The discrete Glauber-Sudarshan, Wigner, and Husimi functions emerge from this formalism as specific cases of s-parametrized phase-space functions where, in particular, a hierarchical process among them is promptly established. In addition, a phase-space description of quantum tomography and quantum teleportation is presented and new results are obtained.

Discrete squeezed states for finite-dimensional spaces

We show how discrete squeezed states in an N-2-dimensional phase space can be properly constructed out of the finite-dimensional context. Such discrete extensions are then applied to the framework of quantum tomography and quantum information theory with the aim of establishing an initial study on the interference effects between discrete variables in a finite phase space. Moreover, the interpretation of the squeezing effects is seen to be direct in the present approach, and has some potential applications in different branches of physics.

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)

Multiresolution spherical wavelet analysis in global seismic tomography

Every seismic event produces seismic waves which travel throughout the Earth. Seismology is the science of interpreting measurements to derive information about the structure of the Earth. Seismic tomography is the most powerful tool for determination of 3D structure of deep Earth's interiors. Tomographic models obtained at the global and regional scales are an underlying tool for determination of geodynamical state of the Earth, showing evident correlation with other geophysical and geological characteristics. The global tomographic images of the Earth can be written as a linear combinations of basis functions from a specifically chosen set, defining the model parameterization. A number of different parameterizations are commonly seen in literature: seismic velocities in the Earth have been expressed, for example, as combinations of spherical harmonics or by means of the simpler characteristic functions of discrete cells. With this work we are interested to focus our attention on this aspect, evaluating a new type of parameterization, performed by means of wavelet functions. It is known from the classical Fourier theory that a signal can be expressed as the sum of a, possibly infinite, series of sines and cosines. This sum is often referred as a Fourier expansion. The big disadvantage of a Fourier expansion is that it has only frequency resolution and no time resolution. The Wavelet Analysis (or Wavelet Transform) is probably the most recent solution to overcome the shortcomings of Fourier analysis. The fundamental idea behind this innovative analysis is to study signal according to scale. Wavelets, in fact, are mathematical functions that cut up data into different frequency components, and then study each component with resolution matched to its scale, so they are especially useful in the analysis of non stationary process that contains multi-scale features, discontinuities and sharp strike. Wavelets are essentially used in two ways when they are applied in geophysical process or signals studies: 1) as a basis for representation or characterization of process; 2) as an integration kernel for analysis to extract information about the process. These two types of applications of wavelets in geophysical field, are object of study of this work. At the beginning we use the wavelets as basis to represent and resolve the Tomographic Inverse Problem. After a briefly introduction to seismic tomography theory, we assess the power of wavelet analysis in the representation of two different type of synthetic models; then we apply it to real data, obtaining surface wave phase velocity maps and evaluating its abilities by means of comparison with an other type of parametrization (i.e., block parametrization). For the second type of wavelet application we analyze the ability of Continuous Wavelet Transform in the spectral analysis, starting again with some synthetic tests to evaluate its sensibility and capability and then apply the same analysis to real data to obtain Local Correlation Maps between different model at same depth or between different profiles of the same model.

Image reconstruction in full-field Fourier-domain optical coherence tomography

Full-field Fourier-domain optical coherence tomography (3F-OCT) is a full-field version of spectral domain/swept source optical coherence tomography. A set of two-dimensional Fourier holograms is recorded at discrete wavenumbers spanning the swept source tuning range. The resultant three-dimensional data cube contains comprehensive information on the three-dimensional spatial properties of the sample, including its morphological layout and optical scatter. The morphological layout can be reconstructed in software via three-dimensional discrete Fourier transformation. The spatial resolution of the 3F-OCT reconstructed image, however, is degraded due to the presence of a phase cross-term, whose origin and effects are addressed in this paper. We present a theoretical and experimental study of the imaging performance of 3F-OCT, with particular emphasis on elimination of the deleterious effects of the phase cross-term.

Full-field Fourier-domain optical coherence tomography

Full-field Fourier-domain optical coherence tomography (3F-OCT) is a full-field version of spectraldomain/swept-source optical coherence tomography. A set of two-dimensional Fourier holograms is recorded at discrete wavenumbers spanning the swept-source tuning range. The resultant three-dimensional data cube contains comprehensive information on the three-dimensional morphological layout of the sample that can be reconstructed in software via three-dimensional discrete Fourier-transform. This method of recording of the OCT signal confers signal-to-noise ratio improvement in comparison with "flying-spot" time-domain OCT. The spatial resolution of the 3F-OCT reconstructed image, however, is degraded due to the presence of a phase cross-term, whose origin and effects are addressed in this paper. We present theoretical and experimental study of imaging performance of 3F-OCT, with particular emphasis on elimination of the deleterious effects of the phase cross-term.

Applications Of Visual Evoked Potentials And Fourier-domain Optical Coherence Tomography In Parkinson's Disease: A Controlled Study.

The goal of this cross-sectional observational study was to quantify the pattern-shift visual evoked potentials (VEP) and the thickness as well as the volume of retinal layers using optical coherence tomography (OCT) across a cohort of Parkinson's disease (PD) patients and age-matched controls. Forty-three PD patients and 38 controls were enrolled. All participants underwent a detailed neurological and ophthalmologic evaluation. Idiopathic PD cases were included. Cases with glaucoma or increased intra-ocular pressure were excluded. Patients were assessed by VEP and high-resolution Fourier-domain OCT, which quantified the inner and outer thicknesses of the retinal layers. VEP latencies and the thicknesses of the retinal layers were the main outcome measures. The mean age, with standard deviation (SD), of the PD patients and controls were 63.1 (7.5) and 62.4 (7.2) years, respectively. The patients were predominantly in the initial Hoehn-Yahr (HY) disease stages (34.8% in stage 1 or 1.5, and 55.8 % in stage 2). The VEP latencies and the thicknesses as well as the volumes of the retinal inner and outer layers of the groups were similar. A negative correlation between the retinal thickness and the age was noted in both groups. The thickness of the retinal nerve fibre layer (RNFL) was 102.7 μm in PD patients vs. 104.2 μm in controls. The thicknesses of retinal layers, VEP, and RNFL of PD patients were similar to those of the controls. Despite the use of a representative cohort of PD patients and high-resolution OCT in this study, further studies are required to establish the validity of using OCT and VEP measurements as the anatomic and functional biomarkers for the evaluation of retinal and visual pathways in PD patients.

Effects Of Misalignments In The Retinal Nerve Fiber Layer Thickness Measurements With Spectral Domain Optical Coherence Tomography.

Purpose. To investigate misalignments (MAs) on retinal nerve fiber layer thickness (RNFLT) measurements obtained with Cirrus(©) SD-OCT. Methods. This was a retrospective, observational, cross-sectional study. Twenty-seven healthy and 29 glaucomatous eyes of 56 individuals with one normal exam and another showing MA were included. MAs were defined as an improper alignment of vertical vessels in the en face image. MAs were classified in complete MA (CMA) and partial MA (PMA), according to their site: 1 (superior, outside the measurement ring (MR)), 2 (superior, within MR), 3 (inferior, within MR), and 4 (inferior, outside MR). We compared RNFLT measurements of aligned versus misaligned exams in all 4 sectors, in the superior area (sectors 1 + 2), inferior area (sectors 3 + 4), and within the measurement ring (sectors 2 + 3). Results. RNFLT measurements at 12 clock-hour of eyes with MAs in the superior area (sectors 1 + 2) were significantly lower than those obtained in the same eyes without MAs (P = 0.043). No significant difference was found in other areas (sectors 1 + 2 + 3 + 4, sectors 3 + 4, and sectors 2 + 3). Conclusion. SD-OCT scans with superior MAs may present lower superior RNFLT measurements compared to aligned exams.

Vertical Bone Measurements From Cone Beam Computed Tomography Images Using Different Software Packages.

This article aimed at comparing the accuracy of linear measurement tools of different commercial software packages. Eight fully edentulous dry mandibles were selected for this study. Incisor, canine, premolar, first molar and second molar regions were selected. Cone beam computed tomography (CBCT) images were obtained with i-CAT Next Generation. Linear bone measurements were performed by one observer on the cross-sectional images using three different software packages: XoranCat®, OnDemand3D® and KDIS3D®, all able to assess DICOM images. In addition, 25% of the sample was reevaluated for the purpose of reproducibility. The mandibles were sectioned to obtain the gold standard for each region. Intraclass coefficients (ICC) were calculated to examine the agreement between the two periods of evaluation; the one-way analysis of variance performed with the post-hoc Dunnett test was used to compare each of the software-derived measurements with the gold standard. The ICC values were excellent for all software packages. The least difference between the software-derived measurements and the gold standard was obtained with the OnDemand3D and KDIS3D (-0.11 and -0.14 mm, respectively), and the greatest, with the XoranCAT (+0.25 mm). However, there was no statistical significant difference between the measurements obtained with the different software packages and the gold standard (p> 0.05). In conclusion, linear bone measurements were not influenced by the software package used to reconstruct the image from CBCT DICOM data.

Towards A Deeper Understanding Of Structural Biomass Recalcitrance Using Phase-contrast Tomography.

The development of technological routes to convert lignocellulosic biomass to liquid fuels requires an in-depth understanding of the cell wall architecture of substrates. Essential pretreatment processes are conducted to reduce biomass recalcitrance and usually increase the reactive surface area. Quantitative three-dimensional information about both bulk and surface structural features of substrates needs to be obtained to expand our knowledge of substrates. In this work, phase-contrast tomography (PCT) was used to gather information about the structure of a model lignocellulosic biomass (piassava fibers). The three-dimensional cellular organization of piassava fibers was characterized by PCT using synchrotron radiation. This technique enabled important physical features that describe the substrate piassava fibers to be visualized and quantified. The external surface area of a fiber and internal surface area of the pores in a fiber could be determined separately. More than 96% of the overall surface area available to enzymes was in the bulk substrate. The pore surface area and length exhibited a positive linear relationship, where the slope of this relationship depended on the plant tissue. We demonstrated that PCT is a powerful tool for the three-dimensional characterization of the cell wall features related to biomass recalcitrance. Original and relevant quantitative information about the structural features of the analyzed material were obtained. The data obtained by PCT can be used to improve processing routes to efficiently convert biomass feedstock into sugars.

Screening Of Miners And Millers At Decreasing Levels Of Asbestos Exposure: Comparison Of Chest Radiography And Thin-section Computed Tomography.

Chest radiography (CXR) is inferior to Thin-section computed tomography in the detection of asbestos related interstitial and pleural abnormalities. It remains unclear, however, whether these limitations are large enough to impair CXR´s ability in detecting the expected reduction in the frequency of these asbestos-related abnormalities (ARA) as exposure decreases. Clinical evaluation, CXR, Thin-section CT and spirometry were obtained in 1418 miners and millers who were exposed to progressively lower airborne concentrations of asbestos. They were separated into four groups according to the type, period and measurements of exposure and/or procedures for controlling exposure: Group I (1940-1966/tremolite and chrysotile, without measurements of exposure and procedures for controlling exposure); Group II (1967-1976/chrysotile only, without measurements of exposure and procedures for controlling exposure); Group III (1977-1980/chrysotile only, initiated measurements of exposure and procedures for controlling exposure) and Group IV (after 1981/chrysotile only, implemented measurements of exposure and a comprehensive procedures for controlling exposure). In all groups, CXR suggested more frequently interstitial abnormalities and less frequently pleural plaques than observed on Thin-section CT (p<0.050). The odds for asbestosis in groups of decreasing exposure diminished to greater extent at Thin-section CT than on CXR. Lung function was reduced in subjects who had pleural plaques evident only on Thin-section CT (p<0.050). In a longitudinal evaluation of 301 subjects without interstitial and pleural abnormalities on CXR and Thin-section CT in a previous evaluation, only Thin-section CT indicated that these ARA reduced as exposure decreased. CXR compared to Thin-section CT was associated with false-positives for interstitial abnormalities and false-negatives for pleural plaques, regardless of the intensity of asbestos exposure. Also, CXR led to a substantial misinformation of the effects of the progressively lower asbestos concentrations in the occurrence of asbestos-related diseases in miners and millers.

Influence of intracanal post on apical periodontitis identified by cone-beam computed tomography

The determination of the success of endodontic treatment has been often discussed based on outcome obtained by periapical radiography. The aim of this study was to verify the influence of intracanal post on apical periodontitis detected by cone-beam computed tomography (CBCT). A consecutive sample of 1020 images (periapical radiographs and CBCT scans) taken from 619 patients (245 men; mean age, 50.1 years) between February 2008 and September 2009 were used in this study. Presence and intracanal post length (short, medium and long) were associated with apical periodontitis (AP). Chi-square test was used for statistical analyses. Significance level was set at p<0.01. The kappa value was used to assess examiner variability. From a total of 591 intracanal posts, AP was observed in 15.06%, 18.78% and 7.95% using periapical radiographs, into the different lengths, short, medium and long, respectively (p=0.466). Considering the same posts length it was verified AP in 24.20%, 26.40% and 11.84% observed by CBCT scans, respectively (p=0.154). From a total of 1,020 teeth used in this study, AP was detected in 397 (38.92%) by periapical radiography and in 614 (60.19%) by CBCT scans (p<0.001). The distribution of intracanal posts in different dental groups showed higher prevalence in maxillary anterior teeth (54.79%). Intracanal posts lengths did not influenced AP. AP was detected more frequently when CBCT method was used.