Rigid motion correction of dual opposed planar projections in single photon imaging

Awake and/or freely moving small animal single photon emission imaging allows the continuous study of molecules exhibiting slow kinetics without the need to restrain or anaesthetise the animals. Estimating motion free projections in freely moving small animal planar imaging can be considered as a limited angle tomography problem, except that we wish to estimate the 2D planar projections rather than the 3D volume, where the angular sampling in all three axes depends on the rotational motion of the animal. In this study, we hypothesise that the motion corrected planar projections estimated by reconstructing an estimate of the 3D volume using an iterative motion compensating reconstruction algorithm and integrating it along the projection path, will closely match the true, motion-less, planar distribution regardless of the object motion. We tested this hypothesis for the case of rigid motion using Monte-Carlo simulations and experimental phantom data based on a dual opposed detector system, where object motion was modelled with 6 degrees of freedom. In addition, we investigated the quantitative accuracy of the regional activity extracted from the geometric mean of opposing motion corrected planar projections. Results showed that it is feasible to estimate qualitatively accurate motion-corrected projections for a wide range of motions around all 3 axes. Errors in the geometric mean estimates of regional activity were relatively small and within 10% of expected true values. In addition, quantitative regional errors were dependent on the observed motion, as well as on the surrounding activity of overlapping organs. We conclude that both qualitatively and quantitatively accurate motion-free projections of the tracer distribution in a rigidly moving object can be estimated from dual opposed detectors using a correction approach within an iterative reconstruction framework and we expect this approach can be extended to the case of non-rigid motion.