PDE的一个除噪程序

clc;clear all;close all; I = imread('finger.bmp');I = double(I); figure,imshow(uint8(I)),title('Orignal image'); y = I; sigma = 1;% need some modification according image rol = 3;% need some modification according image dt = 0.1; % timestep, more smaller, more better result gamma = 0.000001; % lamd