TENSOR TOMOGRAPHY

UoM administered thesis: Phd

  • Authors:
  • Naeem Desai

Abstract

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 Polarimet- ric 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 Trans- verse Ray Transform (TTRT) where we highlight some results on the complete inver- sion procedure, SVD and mention stability results in Sobolev spaces. Thereafter we derive some relations between these transforms. Next we discuss the imaging modali- ties 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 im- age 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.

Details

Original languageEnglish
Awarding Institution
Supervisors/Advisors
Award date1 Aug 2018