Adjoint image warping using multivariate splines with application to four-dimensional computed tomography.

PURPOSE: Adjoint image warping is an important tool to solve image reconstruction problems that warp the unknown image in the forward model. This includes four-dimensional computed tomography (4D-CT) models in which images are compared against recorded projection images of various time frames using image warping as a model of the motion. The inversion of these models requires the adjoint of image warping, which up to now has been substituted by approximations. We introduce an efficient implementation of the exact adjoints of multivariate spline based image warping, and compare it against previously used alternatives. METHODS: Using symbolic computer algebra, we computed a list of 64 polynomials that allow us to compute a matrix representation of trivariate cubic image warping. By combining an on-the-fly computation of this matrix with a parallelized implementation of columnwise matrix multiplication, we obtained an efficient, low memory implementation of the adjoint action of 3D cubic image warping. We used this operator in the solution of a previously proposed 4D-CT reconstruction model in which the image of a single subscan was compared against projection data of multiple subscans by warping and then projecting the image. We compared the properties of our exact adjoint with those of approximate adjoints by warping along inverted motion. RESULTS: Our method requires halve the memory to store motion between subscans, compared to methods that need to compute and store an approximate inverse of the motion. It also avoids the computation time to invert the motion and the tunable parameter of the number of iterations used to perform this inversion. Yet, a similar and often better reconstruction quality was obtained in comparison with these more expensive methods, especially when the motion is large. When compared against a simpler method that is similar to ours in computational demands, our method achieves a higher reconstruction quality in general. CONCLUSIONS: Our implementation of the exact adjoint of cubic image warping improves efficiency and provides accurate reconstructions.

FleXCT: a flexible X-ray CT scanner with 10 degrees of freedom.

Laboratory based X-ray micro-CT is a non-destructive testing method that enables three dimensional visualization and analysis of the internal and external morphology of samples. Although a wide variety of commercial scanners exist, most of them are limited in the number of degrees of freedom to position the source and detector with respect to the object to be scanned. Hence, they are less suited for industrial X-ray imaging settings that require advanced scanning modes, such as laminography, conveyor belt scanning, or time-resolved imaging (4DCT). We introduce a new X-ray scanner FleXCT that consists of a total of ten motorized axes, which allow a wide range of non-standard XCT scans such as tiled and off-centre scans, laminography, helical tomography, conveyor belt, dynamic zooming, and X-ray phase contrast imaging. Additionally, a new software tool 'FlexRayTools' was created that enables reconstruction of non-standard XCT projection data of the FleXCT instrument using the ASTRA Toolbox, a highly efficient and open source set of tools for tomographic projection and reconstruction.