Among various kinds of image denoising methods, the Perona-Malik model is a representative Partial Differential Equation based (PDE-based) algorithm which effectively removes the noise as well as having edge enhancement simultaneously through anisotropic diffusion controlled by the diffusion coefficient. However, the unstable behavior of the Perona-Malik model introduces staircasing artifacts in the processed images. To realize less diffusion in the texture region and to get more smooth in flat region while implementing image denoising, we propose an improved Perona-Malik model based on non-local means theory, which assumes that the image contains an extensive amount of self-similarity and uses the similarity between the region around the c...