We formulate computed tomography (CT) sinogram preprocessing as a statistical restoration problem in which the goal is to obtain the best estimate of the line integrals needed for reconstruction from the set of noisy, degraded measurements. CT measurement data are degraded by a number of factors-including beam hardening and off-focal radiation-that produce artifacts in reconstructed images unless properly corrected. Currently, such effects are addressed by a sequence of sinogram-preprocessing steps, including deconvolution corrections for off-focal radiation, that have the potential to amplify noise. Noise itself is generally mitigated through apodization of the reconstruction kernel, which effectively ignores the measurement statistics, although in high-noise situations adaptive filtering methods that loosely model data statistics are sometimes applied. As an alternative, we present a general imaging model relating the degraded measurements to the sinogram of ideal line integrals and propose to estimate these line integrals by iteratively optimizing a statistically based objective function. We consider three different strategies for estimating the set of ideal line integrals, one based on direct estimation of ideal "monochromatic" line integrals that have been corrected for single-material beam hardening, one based on estimation of ideal "polychromatic" line integrals that can be readily mapped to monochromatic line integrals, and one based on estimation of ideal transmitted intensities, from which ideal, monochromatic line integrals can be readily estimated. The first two approaches involve maximization of a penalized Poisson-likelihood objective function while the third involves minimization of a quadratic penalized weighted least squares (PWLS) objective applied in the transmitted intensity domain. We find that at low exposure levels typical of those being considered for screening CT, the Poisson-likelihood based approaches outperform the PWLS objective as well as a standard approach based on adaptive filtering followed by deconvolution. At higher exposure levels, the approaches all perform similarly.

We have developed a sinogram smoothing approach for low-dose computed tomography (CT) that seeks to estimate the line integrals needed for reconstruction from the noisy measurements by maximizing a penalized-likelihood objective function. The maximization is performed by an algorithm derived by use of the separable paraboloidal surrogates framework. The approach overcomes some of the computational limitations of a previously proposed spline-based penalized-likelihood sinogram smoothing approach, and it is found to yield better resolution-variance tradeoffs than this spline-based approach as well an existing adaptive filtering approach. Such sinogram smoothing approaches could be valuable when applied to the low-dose data acquired in CT screening exams, such as those being considered for lung-nodule detection.

Although relationships among the major groups of living gnathostomes are well established, the relatedness of early jawed vertebrates to modern clades is intensely debated. Here, we provide a new description of , a Middle Devonian (Givetian approx. 385-million-year-old) stem chondrichthyan from Germany, and one of the very few early chondrichthyans in which substantial portions of the endoskeleton are preserved. Tomographic and histological techniques reveal new details of the gill skeleton, hyoid arch and jaws, neurocranium, cartilage, scales and teeth. Despite many features resembling placoderm or osteichthyan conditions, phylogenetic analysis confirms as a stem chondrichthyan and corroborates hypotheses that all acanthodians are stem chondrichthyans. The unfamiliar character combination displayed by , alongside conditions observed in acanthodians, implies that pre-Devonian stem chondrichthyans are severely under-sampled and strongly supports indications from isolated scales that the gnathostome crown group originated at the latest by the early Silurian (approx. 440 Ma). Moreover, phylogenetic results highlight the likely convergent evolution of conventional chondrichthyan conditions among earliest members of this primary gnathostome division, while skeletal morphology points towards the likely suspension feeding habits of, suggesting a functional origin of the gill slit condition characteristic of the vast majority of living and fossil chondrichthyans.

Photoacoustic tomography (PAT), also known as optoacoustic or thermoacoustic tomography, is a hybrid imaging technique that possesses great potential for a wide range of biomedical imaging applications. Image reconstruction in PAT is tantamount to solving an inverse source problem, where the source represents the optical energy absorption distribution in the object that is induced by an interrogating pulsed optical waveform. In this work, we re-examine the PAT image reconstruction problem from a Fourier domain perspective by use of established time-harmonic inverse source concepts. A mathematical relationship between the photoacoustic pressure wavefield data on an aperture that encloses the object and the three-dimensional Fourier transform of the optical absorption distribution evaluated on a collection of concentric spheres is investigated. In addition to providing a framework for deriving both exact and approximate analytic reconstruction formulae, we demonstrate that this mapping provides an intuitive means of understanding certain spatial resolution characteristics of PAT.

We present a statistically principled sinogram smoothing approach for X-ray computed tomography (CT) with the intent of reducing noise-induced streak artifacts. These artifacts arise in CT when some subset of the transmission measurements capture relatively few photons because of high attenuation along the measurement lines. Attempts to reduce these artifacts have focused on the use of adaptive filters that strive to tailor the degree of smoothing to the local noise levels in the measurements. While these approaches involve loose consideration of the measurement statistics to determine smoothing levels, they do not explicitly model the statistical distributions of the measurement data. In this paper, we present an explicitly statistical approach to sinogram smoothing in the presence of photon-starved measurements. It is an extension of a nonparametric sinogram smoothing approach using penalized Poisson-likelihood functions that we have previously developed for emission tomography. Because the approach explicitly models the data statistics, it is naturally adaptive--it will smooth more variable measurements more heavily than it does less variable measurements. We find that it significantly reduces streak artifacts and noise levels without comprising image resolution.