Computed tomography from X-rays: old 2-D results, new 3-D problems

We consider old results on 2-D computerized tomography methods and their relevance to new fully 3-D problems. We examine the 2-D filtered back-projection method (FBP) from several perspectives to better understand how it works. Based on that understanding, we question whether stable, reliable reconstruction algorithms can be found for wide-detector cone-beam 3-D machines with their resulting large-slant beams. We use a numerical method of "point response function fitting" to compute convolution kernels H(u, v) (for use with back-projection) for a model slant-beam problem. These kernels exhibit disturbing growing and spreading oscillations which would greatly amplify errors in the projection data

Statistical Image Reconstruction for High-Throughput Thermal Neutron Computed Tomography

Neutron Computed Tomography (CT) is an increasingly utilised non-destructive analysis tool in material science, palaeontology, and cultural heritage. With the development of new neutron imaging facilities (such as DINGO, ANSTO, Australia) new opportunities arise to maximise their performance through the implementation of statistically driven image reconstruction methods which have yet to see wide scale application in neutron transmission tomography. This work outlines the implementation of a convex algorithm statistical image reconstruction framework applicable to the geometry of most neutron tomography instruments with the aim of obtaining similar imaging quality to conventional ramp filtered back-projection via the inverse Radon transform, but using a lower number of measured projections to increase object throughput. Through comparison of the output of these two frameworks using a tomographic scan of a known 3 material cylindrical phantom obtain with the DINGO neutron radiography instrument (ANSTO, Australia), this work illustrates the advantages of statistical image reconstruction techniques over conventional filter back-projection. It was found that the statistical image reconstruction framework was capable of obtaining image estimates of similar quality with respect to filtered back-projection using only 12.5% the number of projections, potentially increasing object throughput at neutron imaging facilities such as DINGO eight-fold

Three-dimensional image reconstruction in J-PET using Filtered Back Projection method

We present a method and preliminary results of the image reconstruction in the Jagiellonian PET tomograph. Using GATE (Geant4 Application for Tomographic Emission), interactions of the 511 keV photons with a cylindrical detector were generated. Pairs of such photons, flying back-to-back, originate from e+e- annihilations inside a 1-mm spherical source. Spatial and temporal coordinates of hits were smeared using experimental resolutions of the detector. We incorporated the algorithm of the 3D Filtered Back Projection, implemented in the STIR and TomoPy software packages, which differ in approximation methods. Consistent results for the Point Spread Functions of ~5/7,mm and ~9/20, mm were obtained, using STIR, for transverse and longitudinal directions, respectively, with no time of flight information included.Comment: Presented at the 2nd Jagiellonian Symposium on Fundamental and Applied Subatomic Physics, Krak\'ow, Poland, June 4-9, 2017. To be published in Acta Phys. Pol.

Elementary test for non-classicality based on measurements of position and momentum

We generalise a non-classicality test described by Kot et al. [Phys. Rev. Lett. 108, 233601 (2010)], which can be used to rule out any classical description of a physical system. The test is based on measurements of quadrature operators and works by proving a contradiction with the classical description in terms of a probability distribution in phase space. As opposed to the previous work, we generalise the test to include states without rotational symmetry in phase space. Furthermore, we compare the performance of the non-classicality test with classical tomography methods based on the inverse Radon transform, which can also be used to establish the quantum nature of a physical system. In particular, we consider a non-classicality test based on the so-called filtered back-projection formula. We show that the general non-classicality test is conceptually simpler, requires less assumptions on the system and is statistically more reliable than the tests based on the filtered back-projection formula. As a specific example, we derive the optimal test for a quadrature squeezed single photon state and show that the efficiency of the test does not change with the degree of squeezing

Image Reconstruction Techniques and Measure of Quality: Classical vs.Modern Approaches

Mathematical methods are of central importance in the new technologies of image reconstruction. Some of the most important procedures are classified as Back-projection, Filtered Back-projection and iterative reconstruction techniques. Back- projection played an important historical role but is no longer used because of sizable artifacts. Analytical methods like Filtered Back projection excel in speed and accuracy when a large number ofprojections are available. These are extensively used in x-ray imaging. Algebraic Reconstruction Technique(ART) is more attractive when the number of views is limited and when noise is significant. For these reasons, iterative methods are widely used in imaging. Two slight variants of ART are SIRT (Simultaneous Iterative reconstruction Techniques) and SART (Simultaneous ART).A modern method of image reconstruction technique is Fast Slant Slack(FSS). This method is rapidly computable, algebraically exact, geometrically faithful and invertible. A new software known as beamlab is used for FSS image reconstruction. All these reconstruction techniques are explored in this work. Also, various tasks are performed to measure the immunity to noise and quality of the images using PSNR, MSE and Universal Image Quality Index