Tables for
Volume B
Reciprocal space
Edited by U. Shmueli

International Tables for Crystallography (2006). Vol. B. ch. 2.5, p. 318   | 1 | 2 |

Section The method of back-projection

B. K. Vainshteinc The method of back-projection

| top | pdf |

This method is also called the synthesis of projection functions. Let us consider a two-dimensional case and stretch along [\tau_{\psi_{i}}] each one-dimensional projection [L^{i}] (Fig.[link]) by a certain length b; thus, we obtain the projection function [L^{i} ({\bf x}) = {1 \over b} L^{i} (x_{i}) \cdot 1(\tau_{i}). \eqno(] Let us now superimpose h functions [L^{i}] [{\textstyle\sum\limits_{i = 1}^{h}} L^{i} ({\bf x}) = \Sigma_{2} ({\bf x}). \eqno(] The continuous sum over the angles of projection synthesis is [\eqalignno{ \Sigma_{2} ({\bf x}) &= {\textstyle\int\limits_{0}^{\pi}} L(\psi, {\bf x}) \ \hbox{d}\psi = \rho_{2} ({\bf x}) * |{\bf x}|^{-1} &\cr &\simeq {\textstyle\sum\limits_{i = 1}^{h}} L^{i} = \rho_{2} ({\bf x}) + B(1)\hbox{;} &(\cr}] this is the convolution of the initial function with a rapidly falling function [|{\bf x}|^{-1}] (Vainshtein, 1971b[link]). In ([link], the approximation for a discrete set of h projections is also written. Since the function [|{\bf x}|^{-1}] approaches infinity at [x = 0], the convolution with it will reproduce the initial function [\rho ({\bf x})], but with some background B decreasing around each point according to the law [|{\bf x}|^{-1}]. At orthoaxial projection the superposition of cross sections [\varphi_{2} ({\bf x}, z_{k})] arranged in a pile gives the three-dimensional structure [\varphi_{3}].


Figure | top | pdf |

(a) Formation of a projection function; (b) superposition of these functions.

Radon operator. Radon (1917[link]; see also Deans, 1983[link]) gave the exact solution of the problem of reconstruction. However, his mathematical work was for a long time unknown to investigators engaged in reconstruction of a structure from images; only in the early 1970s did some authors obtain results analogous to Radon's (Ramachandran & Lakshminarayanan, 1971[link]; Vainshtein & Orlov, 1972[link], 1974[link]; Gilbert, 1972a[link]).

The convolution in ([link] may be eliminated using the Radon integral operator, which modifies projections by introducing around each point the negative values which annihilate on superposition the positive background values. The one-dimensional projection modified with the aid of the Radon operator has the form [\tilde{L} (x_{\psi}) = {1 \over 2\pi^{2}} \int\limits_{0}^{\infty} {2L(x_{\psi}) - L(x_{\psi} + x'_{\psi}) - L(x_{\psi} - x'_{\psi}) \over x_{\psi}^{'2}} \hbox{ d}x'_{\psi}. \eqno(] Now [\varphi_{2} ({\bf x})] is calculated analogously to ([link], not from the initial projections L but from the modified projection [\tilde{L}]: [\varphi_{2} ({\bf x}) = {\textstyle\int\limits_{0}^{\pi}}\; \tilde{L} (\psi, {\bf x}) \ \hbox{d}\psi \simeq {\textstyle\sum\limits_{i = 1}^{k}}\; \tilde{L}_{i} (\psi_{i}, {\bf x}). \eqno(]

The reconstruction of high-symmetry structures, in particular helical ones, by the direct method is carried out from one projection making use of its equivalence to many projections. The Radon formula in discrete form can be obtained using the double Fourier transformation and convolution (Ramachandran & Lakshminarayanan, 1971[link]).


First citation Deans, S. R. (1983). The Radon transform and some of its applications. New York: John Wiley.Google Scholar
First citation Gilbert, P. F. C. (1972a). The reconstruction of a three-dimensional structure from projections and its application to electron microscopy. II. Direct methods. Proc. R. Soc. London Ser. B, 182, 89–102.Google Scholar
First citation Radon, J. (1917). Über die Bestimmung von Funktionen durch ihre Integralwerte längs gewisser Mannigfaltigkeiten. (On the determination of functions from their integrals along certain manifolds). Ber. Verh. Saechs. Akad. Wiss. Leipzig Math. Phys. Kl. 69, 262–277.Google Scholar
First citation Ramachandran, G. N. & Lakshminarayanan, A. V. (1971). Three-dimensional reconstruction from radiographs and electron micrographs: application of convolutions instead of Fourier transforms. Proc. Natl Acad. Sci. USA, 68(9), 2236–2240.Google Scholar
First citation Vainshtein, B. K. (1971b). Finding the structure of objects from projections. Sov. Phys. Crystallogr. 15, 781–787.Google Scholar
First citation Vainshtein, B. K. & Orlov, S. S. (1972). Theory of the recovery of functions from their projections. Sov. Phys. Crystallogr. 17, 213–216.Google Scholar
First citation Vainshtein, B. K. & Orlov, S. S. (1974). General theory of direct 3D reconstruction. Proceedings of International Workshop, Brookhaven National Laboratory, pp. 158–164.Google Scholar

to end of page
to top of page