Source author record

William R. B. Lionheart

William R. B. Lionheart appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Recovering the second moment of the strain distribution from neutron Bragg edge data

Point by point strain scanning is often used to map the residual stress (strain) in engineering materials and components. However, the gauge volume and hence spatial resolution is limited by the beam defining apertures and can be anisotropic for very low and high diffraction (scattering) angles. Alternatively, wavelength resolved neutron transmission imaging has a potential to retrieve information tomographically about residual strain induced within materials through measurement in transmission of Bragg edges - crystallographic fingerprints whose locations and shapes depend on microstructure and strain distribution. In such a case the spatial resolution is determined by the geometrical blurring of the measurement setup and the detector point spread function. Mathematically, reconstruction of strain tensor field is described by the longitudinal ray transform; this transform has a non-trivial null-space, making direct inversion impossible. A combination of the longitudinal ray transform with physical constraints was used to reconstruct strain tensor fields in convex objects. To relax physical constraints and generalise reconstruction, a recently introduced concept of histogram tomography can be employed. Histogram tomography relies on our ability to resolve the distribution of strain in the beam direction, as we discuss in the paper. More specifically, Bragg edge strain tomography requires extraction of the second moment (variance about zero) of the strain distribution which has not yet been demonstrated in practice. In this paper we verify experimentally that the second moment can be reliably measured for a previously well characterised aluminium ring and plug sample. We compare experimental measurements against numerical calculation and further support our conclusions by rigorous uncertainty quantification of the estimated mean and variance of the strain distribution.

preprint2021arXiv

Crystalline phase discriminating neutron tomography using advanced reconstruction methods

Time-of-flight neutron imaging offers complementary attenuation contrast to X-ray computed tomography (CT), coupled with the ability to extract additional information from the variation in attenuation as a function of neutron energy (time of flight) at every point (voxel) in the image. In particular Bragg edge positions provide crystallographic information and therefore enable the identification of crystalline phases directly. Here we demonstrate Bragg edge tomography with high spatial and spectral resolution. We propose a new iterative tomographic reconstruction method with a tailored regularisation term to achieve high quality reconstruction from low-count data, where conventional filtered back-projection (FBP) fails. The regularisation acts in a separated mode for spatial and spectral dimensions and favours characteristic piece-wise constant and piece-wise smooth behaviour in the respective dimensions. The proposed method is compared against FBP and a state-of-the-art regulariser for multi-channel tomography on a multi-material phantom. The proposed new regulariser which accommodates specific image properties outperforms both conventional and state-of-the-art methods and therefore facilitates Bragg edge fitting at the voxel level. The proposed method requires significantly shorter exposure to retrieve features of interest. This in turn facilitates more efficient usage of expensive neutron beamline time and enables the full utilisation of state-of-the-art high resolution detectors.

preprint2020arXiv

Lipschitz stability at the boundary for time-harmonic diffuse optical tomography

We study the inverse problem in Optical Tomography of determining the optical properties of a medium $Ω\subset\mathbb{R}^n$, with $n\geq 3$, under the so-called diffusion approximation. We consider the time-harmonic case where $Ω$ is probed with an input field that is modulated with a fixed harmonic frequency $ω=\frac{k}{c}$, where $c$ is the speed of light and $k$ is the wave number. We prove a result of Lipschitz stability of the absorption coefficient $μ_a$ at the boundary $\partialΩ$ in terms of the measurements in the case when the scattering coefficient $μ_s$ is assumed to be known and $k$ belongs to certain intervals depending on some a-priori bounds on $μ_a$, $μ_s$.

preprint2020arXiv

Scanning electron diffraction tomography of strain

Strain engineering is used to obtain desirable materials properties in a range of modern technologies. Direct nanoscale measurement of the three-dimensional strain tensor field within these materials has however been limited by a lack of suitable experimental techniques and data analysis tools. Scanning electron diffraction has emerged as a powerful tool for obtaining two-dimensional maps of strain components perpendicular to the incident electron beam direction. Extension of this method to recover the full three-dimensional strain tensor field has been restricted though by the absence of a formal framework for tensor tomography using such data. Here, we show that it is possible to reconstruct the full non-symmetric strain tensor field as the solution to an ill-posed tensor tomography inverse problem. We then demonstrate the properties of this tomography problem both analytically and computationally, highlighting why incorporating precession to perform scanning precession electron diffraction may be important. We establish a general framework for non-symmetric tensor tomography and demonstrate computationally its applicability for achieving strain tomography with scanning precession electron diffraction data.

preprint2019arXiv

Histogram Tomography

In many tomographic imaging problems the data consist of integrals along lines or curves. Increasingly we encounter "rich tomography" problems where the quantity imaged is higher dimensional than a scalar per voxel, including vectors tensors and functions. The data can also be higher dimensional and in many cases consists of a one or two dimensional spectrum for each ray. In many such cases the data contain not just integrals along rays but the distribution of values along the ray. If this is discretized into bins we can think of this as a histogram. In this paper we introduce the concept of "histogram tomography". For scalar problems with histogram data this holds the possibility of reconstruction with fewer rays. In vector and tensor problems it holds the promise of reconstruction of images that are in the null space of related integral transforms. For scalar histogram tomography problems we show how bins in the histogram correspond to reconstructing level sets of function, while moments of the distribution are the x-ray transform of powers of the unknown function. In the vector case we give a reconstruction procedure for potential components of the field. We demonstrate how the histogram longitudinal ray transform data can be extracted from Bragg edge neutron spectral data and hence, using moments, a non-linear system of partial differential equations derived for the strain tensor. In x-ray diffraction tomography of strain the transverse ray transform can be deduced from the diffraction pattern the full histogram transverse ray transform cannot. We give an explicit example of distributions of strain along a line that produce the same diffraction pattern, and characterize the null space of the relevant transform.

preprint2015arXiv

An explicit reconstruction algorithm for the transverse ray transform of a second rank tensor field from three axis data

We give an explicit plane-by-plane filtered back-projection reconstruction algorithm for the transverse ray transform of symmetric second rank tensor fields on Euclidean 3-space, using data from rotation about three orthogonal axes. We show that in the general case two axis data is insufficient but give an explicit reconstruction procedure for the potential case with two axis data

preprint2015arXiv

Direct high-order edge-preserving regularization for tomographic image reconstruction

In this paper we present a new two-level iterative algorithm for tomographic image reconstruction. The algorithm uses a regularization technique, which we call edge-preserving Laplacian, that preserves sharp edges between objects while damping spurious oscillations in the areas where the reconstructed image is smooth. Our numerical simulations demonstrate that the proposed method outperforms total variation (TV) regularization and it is competitive with the combined TV-L2 penalty. Obtained reconstructed images show increased signal-to-noise ratio and visually appealing structural features. Computer implementation and parameter control of the proposed technique is straightforward, which increases the feasibility of it across many tomographic applications. In this paper, we applied our method to the under-sampled computed tomography (CT) projection data and also considered a case of reconstruction in emission tomography The MATLAB code is provided to support obtained results.