Numerical solution of the identification problem for the attenuated Radon transform

DOI10.1088/0266-5611/15/5/314MaRDI QIDQ4719156zbMATH OpenOpenAlexFDO

Authors A. V. Bronnikov

Publication date 5 October 2000

Published in Inverse Problems (Search for Journal in Brave)

Full work available at URL https://doi.org/10.1088/0266-5611/15/5/314

zbMATH Keywords

attenuated Radon transformconsistency conditionsidentification problemmatrix representationNewton methodnumerical algorithmperformancesingle-photon emission computerized tomographysingular value decompositionTikhonov regularization

Mathematics Subject Classification ID

Radon transform (44A12) Numerical methods for integral transforms (65R10) Biomedical imaging and signal processing (92C55)

Review text

The attenuated Radon transform \[ R\mu f(\varphi,s)\equiv \int_{x\omega= s} \exp\Biggl(- \int^\infty_0 \mu(x+ t\omega^\perp) dt\Biggr) f(x) dx= g(\varphi,s) \] of a source function \(f(x)\), serves a mathematical tool for single-photon emission computerized tomography SPECT. The exponential weighting factor in this equation describes the process of linear photon attenuation, where \(\mu(x)\) is the attenuation coefficient in the medium, which plays the role of a parameter of the transform. The identification problem for the attenuated Radon transform is to find function \(\mu\) from \(g(\varphi,s)\) with unknown \(f(x)\). In this paper a new method for the solution of the identification problem for the attenuated Radon transform is considered. The approach is based on a matrix representation of the transform. The author considers a matrix representation of the transform and formulates the corresponding discrete consistency conditions in the form of the orthogonal projection of the data vector onto the orthogonal complement of the column space of the matrix. The singular value decomposition is applied to compute the orthogonal projector and its Fréchet derivative. The numerical algorithm suggested is based on the Newton method with Tikhonov regularization. To study the performance of the approach, a numerical simulation of a cardiac SPECT imaging experiment is also discussed. The experiment with an anthropomorphic torso phantom demonstrates promising results in the reconstructon of the thorax cross section characterized by relatively large variations of the attenuation map.

Recommendations

• scientific article; zbMATH DE number 3951958
• The identification problem for the exponential Radon transform
• A new approach towards simultaneous activity and attenuation reconstruction in emission tomography

Cited in

(11)

• The identification problem for the constantly attenuated Radon transform
• Adjoint state method for the identification problem in SPECT: recovery of both the source and the attenuation in the attenuated X-ray transform
• The identification problem for the attenuated X-ray transform
• scientific article; zbMATH DE number 3866535 (Why is no real title available?)
• Algebraic Reconstruction of Source and Attenuation in SPECT Using First Scattering Measurements
• A New Numerical Approach to Inverse Transport Equation with Error Analysis
• Fast numerical inversion of the attenuated Radon transform with full and partial measurements
• The identification problem for the exponential Radon transform
• SPECT with a multi-bang assumption on attenuation
• Inversion of the scalar and vector attenuated X-ray transforms in a unit disc
• Simultaneous recovery of attenuation and source density in SPECT

This page was built for publication: Numerical solution of the identification problem for the attenuated Radon transform

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4719156)