Here, we present a method capable of jointly reconstructing the strain field and the \(d_0\)  field from a single set of neutron transmission images.
To achieve this both the strain and \(d_0\)  are modelled by a Gaussian process (see for example [1]}) and equilibrium and boundary conditions are built into the strain model [2]}.
This extends the Gaussian process approach presented by [3]}, [4]} to handle the inherently non-linear nature of this problem.
A numerically tractable algorithm based on variational inference (see for example [5]}, [6]}) is provided and the method is validated on a simulated data set.
