RESUMEN
Objective. To enable practical interferometry-based phase contrast CT using standard incoherent x-ray sources, we propose an imaging system where the analyzer grating is replaced by a high-resolution detector. Since there is no need to perform multiple exposures (with the analyzer grating at different positions) at each scan angle, this scheme is compatible with continuous-rotation CT apparatus, and has the potential to reduce patient radiation dose and patient motion artifacts.Approach. Grating-based x-ray interferometry is a well-studied technique for imaging soft tissues and highly scattering objects embedded in such tissues. In addition to the traditional x-ray absorption-based image, this technique allows reconstruction of the object phase and small-angle scattering information. When using conventional incoherent, polychromatic, hard x-ray tubes as sources, three gratings are usually employed. To sufficiently resolve the pattern generated in these interferometers with contemporary x-ray detectors, an analyzer grating is used, and consequently multiple images need to be acquired for each view angle. This adds complexity to the imaging system, slows image acquisition and thus increases sensitivity to patient motion, and is not dose efficient. By simulating image formation based on wave propagation, and proposing a novel phase retrieval algorithm based on a virtual grating, we assess the potential of a analyzer-grating-free system to overcome these limitations.Main results. We demonstrate that the removal of the analyzer-grating can produce equal image contrast-to-noise ratio at reduced dose (by a factor of 5), without prolonging scan duration.Significance.By demonstrating that an analyzer-free CT system, in conjuction with an efficient phase retrieval algorithm, can overcome the prohibitive dose and workflow penalties associated grating-stepping, an alternative path towards realizing clinical inteferometric CT appears possible.
Asunto(s)
Algoritmos , Interferometría , Humanos , Rayos X , Radiografía , Cintigrafía , Interferometría/métodosRESUMEN
We present an algorithm of reduced computational cost which is able to estimate kinetic model parameters directly from dynamic ECT sinograms made up of temporally inconsistent projections. The algorithm exploits the extreme degree of parameter redundancy inherent in linear combinations of the exponential functions which represent the modes of first-order compartmental systems. The singular value decomposition is employed to find a small set of orthogonal functions, the linear combinations of which are able to accurately represent all modes within the physiologically anticipated range in a given study. The reduced dimension basis is formed as the convolution of this orthogonal set with a measured input function. The Moore-Penrose pseudoinverse is used to find coefficients of this basis. Algorithm performance is evaluated at realistic count rates using MCAT phantom and clinical 99mTc-teboroxime myocardial study data. Phantom data are modelled as originating from a Poisson process. For estimates recovered from a single slice projection set containing 2.5 x 10(5) total counts, recovered tissue responses compare favourably with those obtained using more computationally intensive methods. The corresponding kinetic parameter estimates (coefficients of the new basis) exhibit negligible bias, while parameter variances are low, falling within 30% of the Cramér-Rao lower bound.