Premium
Fast image reconstruction for optical absorption tomography in media with radially symmetric boundaries
Author(s) -
Hampel Uwe,
Freyer Richard
Publication year - 1998
Publication title -
medical physics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.473
H-Index - 180
eISSN - 2473-4209
pISSN - 0094-2405
DOI - 10.1118/1.598160
Subject(s) - iterative reconstruction , inverse problem , tomography , singular value decomposition , optical tomography , imaging phantom , diffuse optical imaging , integral equation , mathematical analysis , physics , mathematics , optics , computer science , algorithm , computer vision
In this paper we present a reconstruction algorithm to invert the linearized problem in optical absorption tomography for objects with radially symmetric boundaries. This is a relevant geometry for functional volume imaging of body regions that are sensitive to ionizing radiation, e.g., breast and testis. From the principles of diffuse light propagation in scattering media we derive the governing integral equations describing the effects of absorption variations on changes in the measurement data. Expansion of these equations into a Neumann series and truncation of higher‐order terms yields the linearized forward imaging operator. For the proposed geometry we utilize an invariance property of this operator, which greatly reduces the problem dimensionality. This allows us to compute the inverse by singular value decomposition and consequently to apply regularization techniques based on the knowledge of the singular value spectrum. The inversion algorithm is highly efficient computing slice images as fast as convolution‐backprojection algorithms in computed tomography (CT). To demonstrate the capacity of the inversion scheme we present reconstruction results for synthetic and phantom measurement data.