The deconvolution of blurred and noisy satellite images is an ill-posed inverse problem, which can be regularized within a Bayesian context by using an a priori model of the reconstructed solution. Since real satellite data show spatially variant characteristics, we propose to use an inhomogeneous model. We use the Maximum Likelihood Estimator (MLE) to estimate its parameters. We demonstrate that the MLE computed on the corrupted image is not suitable for image deconvolution, because it is not robust to noise. Then we show that the estimation is correct only if it is made from the original image. As this image is unknown, we need to compute an approximation of sufficiently good quality to provide useful estimation results. Such an approximation is provided by a wavelet-based deconvolution algorithm. Thus, an hybrid method is first used to estimate the space-variant parameters from this image and second to compute the regularized solution. The obtained results on high resolution satell...