Contourlet transform was introduced in the SAR image statistical property research work. Based on Contourlet domain Hidden Markov Trees Model, starting from image restore methods, a novel SAR image despeckling method was given. This method combined the SAR image despeckling technology with minimum mean square error (MMSE) and Bayes estimate. A new speckle variation estimate method was given based on Laplacian pyramid (LP). The results comparison with wavelet domain method were given. Experimental results show that the method represents better performance in speckle reduction and edges information detection.