Fast and optimal multiframe blind deconvolution algorithm for high-resolution ground-based imaging of space objects.

We report a multiframe blind deconvolution algorithm that we have developed for imaging through the atmosphere. The algorithm has been parallelized to a significant degree for execution on high-performance computers, with an emphasis on distributed-memory systems so that it can be hosted on commodity clusters. As a result, image restorations can be obtained in seconds to minutes. We have compared and quantified the quality of its image restorations relative to the associated Cramér-Rao lower bounds (when they can be calculated). We describe the algorithm and its parallelization in detail, demonstrate the scalability of its parallelization across distributed-memory computer nodes, discuss the results of comparing sample variances of its output to the associated Cramér-Rao lower bounds, and present image restorations obtained by using data collected with ground-based telescopes.

[1]  Richard G. Lane,et al.  Blind deconvolution of noisy complex-valued image , 1989 .

[2]  J. Goodman.  Introduction to Fourier optics , 1969 .

[3]  Ivica Kopriva,et al.  Single frame blind image deconvolution by non-negative sparse matrix factorization , 2006 .

[4]  T J Holmes,et al.  Blind deconvolution of quantum-limited incoherent imagery: maximum-likelihood approach. , 1992, Journal of the Optical Society of America. A, Optics and image science.

[5]  D. Fried.  Optical Resolution Through a Randomly Inhomogeneous Medium for Very Long and Very Short Exposures , 1966 .

[6]  N Wang,et al.  Parallel-distributed blind deconvolution based on a self-organizing neural network. , 1999, Applied optics.

[7]  Jean J Dolne,et al.  Cramer-Rao bound and phase-diversity blind deconvolution performance versus diversity polynomials. , 2005, Applied optics.

[8]  Richard G. Lane,et al.  Automatic multidimensional deconvolution , 1987 .

[9]  A. Lohmann,et al.  Speckle masking in astronomy: triple correlation theory and applications. , 1983, Applied optics.

[10]  William H. Press,et al.  Numerical recipes , 1990 .

[11]  J. C. Dainty,et al.  Iterative blind deconvolution method and its applications , 1988 .

[12]  N Miura,et al.  Parallel scheme of the iterative blind deconvolution method for stellar object reconstruction. , 1993, Applied optics.

[13]  Jorge Nocedal,et al.  Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization , 1997, TOMS.

[14]  Donald R. McGaughey,et al.  Parallel multiframe blind deconvolution using wavelength diversity , 2004, SPIE Optics + Photonics.

[15]  Lewis C. Roberts,et al.  Characterization of the AEOS Adaptive Optics System , 2002 .

[16]  T. M. Cannon,et al.  Blind deconvolution through digital signal processing , 1975, Proceedings of the IEEE.

[17]  Peyman Milanfar,et al.  Statistical performance analysis of super-resolution , 2006, IEEE Transactions on Image Processing.

[18]  Daniel Yellin,et al.  Multichannel system identification and deconvolution: performance bounds , 1999, IEEE Trans. Signal Process..

[19]  A. M. Darling,et al.  Blind deconvolution for referenceless speckle imaging , 1990, Optics & Photonics.

[20]  C. Matson,et al.  Biased Cramér-Rao lower bound calculations for inequality-constrained estimators. , 2006, Journal of the Optical Society of America. A, Optics, image science, and vision.

[22]  Ehud Weinstein,et al.  Maximum likelihood and lower bounds in system identification with non-Gaussian inputs , 1994, IEEE Trans. Inf. Theory.

[23]  Roger Y. Tsai,et al.  Multiframe image restoration and registration , 1984 .

[24]  Michael W. Marcellin,et al.  Iterative multiframe superresolution algorithms for atmospheric-turbulence-degraded imagery , 1998 .

[25]  Nikolas P. Galatsanos,et al.  Projection-based blind deconvolution , 1994 .

[26]  David R Gerwe,et al.  Cramér-Rao analysis of orientation estimation: viewing geometry influences on the information conveyed by target features. , 2003, Journal of the Optical Society of America. A, Optics, image science, and vision.

[27]  Thomas L. Marzetta,et al.  Parameter estimation problems with singular information matrices , 2001, IEEE Trans. Signal Process..

[28]  Timothy J. Schulz,et al.  Multiframe blind deconvolution of astronomical images , 1993 .

[29]  Michael C. Roggemann,et al.  Cramér–Rao analysis of phase-diverse wave-front sensing , 1999 .

[30]  Deepa Kundur,et al.  Blind image deconvolution revisited , 1996 .

[31]  J. C. Dainty,et al.  Stellar Speckle Interferometry , 1976 .

[32]  B. C. McCallum.  Blind deconvolution by simulated annealing , 1990 .

[33]  Jesper Lomborg Jensen,et al.  Multi-processor system for real-time deconvolution and flow estimation in medical ultrasound , 1996, 1996 IEEE Ultrasonics Symposium. Proceedings.

[34]  S. Kay.  Fundamentals of statistical signal processing: estimation theory , 1993 .

[35]  Richard G. Paxman,et al.  Fine-resolution astronomical imaging with phase-diverse speckle , 1993 .

[36]  Charles L Matson,et al.  Primary and secondary superresolution by data inversion. , 2006, Optics express.

[37]  A. Labeyrie.  Attainment of diffraction limited resolution in large telescopes by Fourier analysing speckle patterns in star images , 1970 .

[38]  Naoshi Baba,et al.  Extended-object reconstruction with sequential use of the iterative blind deconvolution method , 1992 .

[39]  Charles L. Matson,et al.  Resolution, linear filtering, and positivity , 1998 .

[40]  Brian J. Thompson,et al.  New methods of processing speckle pattern star images , 1973 .

[41]  Boaz Porat,et al.  Digital Processing of Random Signals: Theory and Methods , 2008 .

[42]  Nikolas P. Galatsanos,et al.  Methods for choosing the regularization parameter and estimating the noise variance in image restoration and their relation , 1992, IEEE Trans. Image Process..

[43]  Yoram Bresler,et al.  Analysis and comparative evaluation of techniques for multichannel blind deconvolution , 1996, Proceedings of 8th Workshop on Statistical Signal and Array Processing.

[44]  Mario Bertero,et al.  Introduction to Inverse Problems in Imaging , 1998 .

[45]  R. Gerchberg.  A practical algorithm for the determination of phase from image and diffraction plane pictures , 1972 .

[46]  O. von der Lühe,et al.  Estimating Fried’s parameter from a time series of an arbitrary resolved object imaged through atmospheric turbulence , 1984 .