In this paper, We present a new image denoising algorithm. We assume a mixture of bivariate circular symmetric Laplacian probability density functions (pdfs) where for each wavelet coefficients may have different local parameter. This pdf characterizes simultaneously 1) the heavy-tailed nature, 2) the interscale dependencies of the wavelet coefficients and also 3) the empirically observed correlation between the coefficient amplitudes. We employ this local bivariate mixture model to derive a Bayesian image denoising technique. This proposed pdf, potentially can fits better the statistical properties of the wavelet coefficients than several other existing models. Our simulation results reveal that the proposed denoising method is among the best reported in the literature. This is justified since the accuracy of the employed distribution for noise-free data determines the denoising performance.