Linear neural network for the solution of limited angle problems in computer-aided tomography

• Published in: Journal of electronic imaging 1998-01, Vol.7 (1), p.70-78
• Main Authors: Wong, Shun Him, Yau, Sze Fong
• Format: Article
• Language: English
• Subjects: Acoustic imaging; Algebraic reconstruction technique (ART); Clark-Palmer-Lawrence interpolation; Computer-aided tomography (CT); Convolution back-projection (CBP); Limited angle problems; Limited angle tomography; Line-integral projections; Linear neural network; Medical imaging; Missing cone problem; Nondestructive evaluation; Papoulis-Gerchberg algorithm; Projections onto convex sets (POCS); Raysums; Sinogram; Synthetic aperture radar
• ISSN: 1017-9909; 1560-229X
• DOI: 10.1117/1.482652

The problem of limited angle tomography in which a complete sinogram is not available is considered. This situation arises in many practical applications where tomographic projection over 180 deg is either physically unrealizable or infeasible. When a complete sinogram is not available, it is well known that the reconstructed images using common reconstruction algorithms, such as convolution back-projection (CBP), will have severe streak artifacts. We present a linear artificial neural network to extrapolate the missing part of the sinogram. Once the complete sinogram is obtained via extrapolation, standard reconstruction techniques such as CBP can be used to generate artifact-free reconstructions. The parameters of the neural network are designed using the sampling theory of signals with noncompact spectral support, the knowledge that complete sinograms have bow-tie-shaped spectral support and regularization. It is found that once designed, these parameters are data independent, especially for images of similar nature. For a sinogram with angular views, each having raysum per view, if views are available, the computational requirement of the neural network is only. The proposed neural network is much more efficient than other iterative algorithms, such as the method of projection onto convex sets; the Papoulis-Gerchberg algorithm; and the Clark-Palmer-Lawrence interpolation method; which requires computations of the order of where is the number of iterations. ©