This paper examines a super-exponential method for blind deconvolution. Possibly non-minimal phase point spread functions (PSFs) are identified. The PSF is assumed to he low pass in nature. No other prior knowledgeof the PSF or the original image is necessary to assure convergence of the algorithm. Results are shown usingsynthetically degraded satellite images in order to demonstrate the accuracy of the PSF estimates. In addition, radiographic images are restored with no knowledge of the PSF of the x-ray imaging system. These experimentssuggest a promising application of this algorithm to a variety of blur identification probleins.Keywords: blur identification, cumulants. higher order statistics. higher order spectra, blind deconvolution 1. INTRODUCTION Blur identification is an important and necessary step in the restoration of images having linear, space invariantdegradations. Since the 1960s. research in blind deconvolution, where there is no assumed knowledge of the degra-dation system's point spread function (PSE), has centered around communication signals. Initial research dealt withsystems whose PSF could he modeled parametrically using zero crossing information preserved in the spectrum of theblurred signal. Only certain classes of blurs. such as cylindrical or i-D motion blurs, could be successfully extractedusing these methods. Stochastic methods were soon introduced. Autoregressive Moving-Average (ARMA) modelsalong with maximum likelihood based approaches became successful due to their ability to deconvolve signals withlow to moderate noise [1, 2].Recently, the use of higher order statistics/spectra. (HOS) has become popular in PSF identification in commu-nication channels [3, 4, 5, 6. 7]. HOS preserve important phase information which is lost when using second orderspectra. Another attractive feature is that given enough realizations of a signal containing additive Gaussian noise,the use of HOS allows us to asi1y filter out the noise. These properties permit us to identify non-minimum phasesystems, such as non-causal blurs.In this work we extend and apply a blur identification algorithm based on HOS, which was previously used forthe 1-D communication signal problem [8], to 2-D radiographic images. Our primary motivation in this work is toprovide fast and accurate restorations of blurred radiographs of fatigued aircraft structures. To do this, we use ablind deconvolution algorithm that converges at a super-exponential rate independent of initialization, asproven in[8]. The algorithm can be applied to the blind deconvolutioii of any type of image, such as those obtained by theHubble Space Telescope, satellites, etc.The general blind deconvolution 1)1ohle1i is illustrated in Figure 1. In this model