Owing to the blurring effect from atmosphere and camera system in the satellite imaging,a blind image restoration algorithm is proposed which includes the modulation transfer function(MTF) estimation and the image restoration.In the MTF estimation stage,based on every degradation process of satellite imaging-chain,a combined parametric model of MTF is given and used to fit the surface of normalized logarithmic amplitude spectrum of degraded image.In the image restoration stage,a maximum a posteriori(MAP) based edge-preserving image restoration method is presented which introduces multivariate Laplacian model to characterize the prior distribution of wavelet coefficients of original image.During the image restoration,in order to avoid solving high nonlinear equations,optimization transfer algorithm is adopted to decompose the image restoration procedure into two simple steps:Landweber iteration and wavelet thresholding denoising.In the numerical experiment,the satellite image restoration results from SPOT-5 and high resolution camera(HR) of China & Brazil earth resource satellite(CBERS-02B) ane compared,and the proposed algorithm is superior in the image edge preservation and noise inhibition.