Numerical Methods for Coupled Reconstruction and Registration in Digital Breast Tomosynthesis

TL;DR

This paper introduces a coupled mathematical framework for joint reconstruction and registration in digital breast tomosynthesis, improving accuracy over traditional sequential methods especially in limited-angle scenarios.

Contribution

The novel joint estimation approach simultaneously recovers image intensities and transformations, incorporating non-rigid B-spline models for enhanced deformation recovery.

Findings

Joint methods outperform sequential approaches in accuracy.

Simultaneous estimation improves reconstruction fidelity.

Non-rigid models effectively recover complex deformations.

Abstract

Digital Breast Tomosynthesis (DBT) provides an insight into the fine details of normal fibroglandular tissues and abnormal lesions by reconstructing a pseudo-3D image of the breast. In this respect, DBT overcomes a major limitation of conventional X-ray mammography by reducing the confounding effects caused by the superposition of breast tissue. In a breast cancer screening or diagnostic context, a radiologist is interested in detecting change, which might be indicative of malignant disease. To help automate this task image registration is required to establish spatial correspondence between time points. Typically, images, such as MRI or CT, are first reconstructed and then registered. This approach can be effective if reconstructing using a complete set of data. However, for ill-posed, limited-angle problems such as DBT, estimating the deformation is complicated by the significant artefacts associated with the reconstruction, leading to severe inaccuracies in the registration. This paper presents a mathematical framework, which couples the two tasks and jointly estimates both image intensities and the parameters of a transformation. We evaluate our methods using various computational digital phantoms, uncompressed breast MR images, and in-vivo DBT simulations. Firstly, we compare both iterative and simultaneous methods to the conventional, sequential method using an affine transformation model. We show that jointly estimating image intensities and parametric transformations gives superior results with respect to reconstruction fidelity and registration accuracy. Also, we incorporate a non-rigid B-spline transformation model into our simultaneous method. The results demonstrate a visually plausible recovery of the deformation with preservation of the reconstruction fidelity.