Modern Mathematics in Biomedical Imaging

Over the past decades, modern mathematics is playing an increasingly important role in the development of biomed-ical imaging, yielding groundbreaking ideas and powerful tools for clinical and preclinical applications. Nevertheless, there are not only growing obstacles between applied mathematicians and imaging engineers to communicate more effectively but also great synergy to be unlocked across boundaries between mathematics and imaging technology. Suggested by the Editor-in-Chief, we organized this special issue to address theseunique challenges by presenting some sophisticated mathematical methods in an engineer-friendly language and demonstrating some representative applications of contemporary analytic tools in image reconstruction research. Nonrigid image registration and a diffeomorphic version of the demons algorithm have been commonly used to estimate tissue deformations in highly deformable anatomies. In the paper by G. Janssens et al. in this issue, the authors further developed the concept of continuous diffeomorphic flow and proposed a diffeomorphic version of the morphons registration method. This method can be used to obtain accurate estimation of deformations between images with variable contrast and hence can be applied to radiotherapy for lung cancer patients and to 4D respiratory-correlated CT of the thorax. Bioluminescence tomography (BLT) has been developed at a fast pace in recent years in study of physiological and pathological processes at cellular and molecular levels. In practice, fine discretization is required but may cause large datasets and an increase the ill-posedness of the problem. In the paper by J. Liu et al. in this issue, the authors present a multilevel sparse reconstruction method based on a framework of the finite-element method. Empirical results showed its effectiveness and potential applications in BLT. When one considers the time-integrated X-ray flux from multiple X-ray sources to shorten the data acquisition process, a promising way is to use overlapped projections from multiple X-ray sources. In order to perform image reconstruction effectively and efficiently from overlapping projections in this configuration, H. Yu et al. developed a multisource simultaneous-algebraic reconstruction technique regularized by a sparsity-oriented constraint in the soft-threshold filtering framework in their paper in this issue. Their numerical simulation further verified the proposed algorithm and demonstrated its advantages in image reconstruction from overlapping data. Mathematical modeling in dynamic positron emission tomography (PET) can be simplified, using compartment models as a linear system. To avoid invasive arterial sampling of blood to acquire values of tracer concentration, blind methods to estimate both blood input and kinetic parameters have recently drawn attention. In …