Tensor Tomography

Desai, Naeem M (2018) Tensor Tomography. Doctoral thesis, University of Manchester.

[thumbnail of thesis.pdf] Text

Download (2MB)


Rich tomography is becoming increasingly popular since we have seen a substantial increase in computational power and storage. Instead of measuring one scalar for each ray, multiple measurements are needed per ray for various imaging modalities. This advancement has allowed the design of experiments and equipment which facilitate a broad spectrum of applications. We present new reconstruction results and methods for several imaging modalities including x-ray diffraction strain tomography, Photoelastic tomography and Polarimetric Neutron Magnetic Field Tomography (PNMFT). We begin with a survey of the Radon and x-ray transforms discussing several procedures for inversion. Furthermore we highlight the Singular Value Decomposition (SVD) of the Radon transform and consider some stability results for reconstruction in Sobolev spaces. We then move onto define the Non-Abelian Ray Transform (NART), Longitudinal Ray Transform (LRT), Transverse Ray Transform (TRT) and the Truncated Transverse Ray Transform (TTRT) where we highlight some results on the complete inversion procedure, SVD and mention stability results in Sobolev spaces. Thereafter we derive some relations between these transforms. Next we discuss the imaging modalities in mind and relate the transforms to their specific inverse problems, primarily being linear. Specifically, NART arises in the formulation of PNMFT where we want to image magnetic structures within magnetic materials with the use of polarized neutrons. After some initial numerical studies we extend the known Radon inversion presented by experimentalists, reconstructing fairly weak magnetic fields, to reconstruct PNMFT data up to phase wrapping. We can recover the strain field tomographically for a polycrystalline material using diffraction data and deduce that a certain moment of that data corresponds to the TRT. Quite naturally the whole strain tensor can be reconstructed from diffraction data measured using rotations about six axes. We develop an innovative explicit plane-by-plane filtered back-projection reconstruction algorithm for the TRT, using data from rotations about three orthogonal axes and state the reasoning why two-axis data is insufficient. For the first time we give the first published results of TRT reconstruction. To complete our discussion we present Photoelastic tomography which relates to the TTRT and implement the algorithm discussing the difficulties that arise in reconstructing data. Ultimately we return to PNMFT highlighting the nonlinear inverse problem due to phase wrapping. We propose an iterative reconstruction algorithm, namely the Modified Newton Kantarovich method (MNK) where we keep the Jacobian (Fréchet derivative) fixed at the first step. However, this is shown to fail for large angles suggesting to develop the Newton Kantarovich (NK) method where we update the Jacobian at each step of the iteration process.

Item Type: Thesis (Doctoral)
Subjects: MSC 2010, the AMS's Mathematics Subject Classification > 45 Integral equations
Depositing User: Dr Naeem Mustaq Desai
Date Deposited: 28 Jan 2018 12:03
Last Modified: 28 Jan 2018 12:03
URI: https://eprints.maths.manchester.ac.uk/id/eprint/2617

Actions (login required)

View Item View Item