Polarimetric neutron tomography of magnetic fields: uniqueness of solution and reconstruction

DOI10.1088/1361-6420/AB44E0MaRDI QIDQ5000574zbMATH OpenOpenAlexFDO

Authors Naeem M. Desai, William R. B. Lionheart, Morten Sales, Markus Strobl, Søren Schmidt

Publication date 14 July 2021

Published in Inverse Problems (Search for Journal in Brave)

Copyright license Creative Commons Attribution 4.0 International

Full work available at URL https://arxiv.org/abs/1903.10884

zbMATH Keywords

reconstruction algorithmuniqueness of solutionRadon inversionnon-abelian ray transformpolarimetric neutron tomography of magnetic fields

Mathematics Subject Classification ID

Radon transform (44A12) Integral geometry (53C65) Numerical methods for integral transforms (65R10) Electromagnetic theory (general) (78A25)

Abstract: We consider the problem of determination of a magnetic field from three dimensional polarimetric neutron tomography data. We see that this is an example of a non-Abelian ray transform and that the problem has a globally unique solution for smooth magnetic fields with compact support, and a locally unique solution for less smooth fields. We derive the linearization of the problem and note that the derivative is injective. We go on to show that the linearised problem about a zero magnetic field reduces to plane Radon transforms and suggest a modified Newton Kantarovich method (MNKM) type algorithm for the numerical solution of the non-linear problem, in which the forward problem is re-solved but the same derivative used each time. Numerical experiments demonstrate that MNKM works for small enough fields (or large enough velocities) and show an example where it fails to reconstruct a slice of the simulated data set. Lastly we show viewed as an optimization problem the inverse problem is non-convex so we expect gradient based methods may fail.

Recommendations

• On iterative reconstruction in the nonlinearized polarization tomography
• Stability of the non-abelian \(X\)-ray transform in dimension \(\ge 3\)
• Generic local uniqueness and stability in polarization tomography
• scientific article; zbMATH DE number 723992
• An approach to direct reconstruction of a solenoidal part in vector and tensor tomography problems

Cites work

• Convergence and Applications of Newton-type Iterations
• Integral geometry for tensor fields. Transl. from the Russian
• ON DETERMINATION OF A GAUGE FIELD ON ℝ^d FROM ITS NON-ABELIAN RADON TRANSFORM ALONG ORIENTED STRAIGHT LINES
• On non-abelian Radon transform
• The attenuated ray transform for connections and Higgs fields

Cited in

(4)

• Histogram tomography
• On iterative reconstruction in the nonlinearized polarization tomography
• The transport Oka-Grauert principle for simple surfaces
• Loop group factorization method for the magnetic and thermostatic nonabelian ray transforms

This page was built for publication: Polarimetric neutron tomography of magnetic fields: uniqueness of solution and reconstruction

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5000574)