Optical microscope tomography. II. Nonnegative constraint by a gradient-projection method

In part I of this series we proposed a new principle of tomographic imaging using an optical microscope devised with a support-constraint reconstruction algorithm [ J. Opt. Soc. Am. A4, 292 ( 1987)]. In this paper we describe a new reconstruction algorithm with a nonnegative constraint that utilizes a priori knowledge of nonnegativity of absorbance or optical density of three-dimensional sample distribution. This method drastically improves the longitudinal resolution in tomography reconstruction, which was the most serious problem with the developed microscope tomography principle. Compared with conventional nonnegative-constrained iterative reconstruction algorithms, the algorithm developed here essentially reduces computation-time and memory-capacity requirements because it includes conjugate directional searching and gradient projection based on the Kuhn–Tucker condition, guaranteeing the convergence of the algorithm within a finite (the smallest) number of iterations. Results of experiments with a biological specimen, using the proposed algorithm, demonstrate an evident improvement of longitudinal resolution with a computation time only three times longer than that of the fastest nonconstrained algorithm.

[1]  B. R. Frieden Image enhancement and restoration , 1979 .

[2]  Satoshi Kawata,et al.  THE PRINCIPLE AND APPLICATIONS OF OPTICAL MICROSCOPE TOMOGRAPHY , 1986 .

[3]  K. C. Tam,et al.  The Use of Multispectral Imaging in Limited-Angle Reconstruction , 1982, IEEE Transactions on Nuclear Science.

[4]  R L Byer,et al.  Experimental optical fan beam tomography. , 1984, Applied optics.

[5]  R. Gerchberg Super-resolution through Error Energy Reduction , 1974 .

[6]  M. Sezan,et al.  Image Restoration by the Method of Convex Projections: Part 2-Applications and Numerical Results , 1982, IEEE Transactions on Medical Imaging.

[7]  R L Byer,et al.  Quantitative optical tomographic imaging of a supersonic jet. , 1986, Optics letters.

[8]  Hiroaki Koga,et al.  Photoelastic computer tomography: a novel measurement method for axial residual stress profile in optical fibers , 1986 .

[9]  D. Youla,et al.  Image Restoration by the Method of Convex Projections: Part 1ߞTheory , 1982, IEEE Transactions on Medical Imaging.

[10]  S. Minami,et al.  Constrained resolution enhancement in optical microscopic tomography , 1986, ICASSP '86. IEEE International Conference on Acoustics, Speech, and Signal Processing.

[11]  Satoshi Kawata,et al.  Optical microscope tomography. I. Support constraint , 1987 .

[12]  B. Stuck A new proposal for estimating the spatial concentration of certain types of air pollutants , 1977 .

[13]  P. Jansson,et al.  Resolution Enhancement of Spectra , 1970 .

[14]  J R Fienup,et al.  Phase retrieval algorithms: a comparison. , 1982, Applied optics.

[15]  Yoshihisa Maruyama,et al.  Measurement of the Refractive Index Distribution in the Interior of a Solid Object from Multi-Directional Interferograms , 1977 .

[16]  B R Myers,et al.  Two-dimensional spectral line emission reconstruction as a plasma diagnostic. , 1978, The Review of scientific instruments.

[17]  O. Nalcioglu,et al.  Constrained Iterative Reconstruction by the Conjugate Gradient Method , 1985, IEEE Transactions on Medical Imaging.

[18]  Michel Defrise,et al.  Regularized iterative and noniterative procedures for object restoration in the presence of noise: an error analysis , 1983 .