Прочтите полутоновое изображение, затем примените к нему сильный Гауссов шум. Отобразите шумное изображение.

I = imread('pout.tif');
Inoisy = imnoise(I,'gaussian',0,0.005);
imshow(Inoisy)
title('Noisy Image')

Оцените порог градиента и количество итераций, необходимых для выполнения анизотропной диффузионной фильтрации изображения.

[gradThresh,numIter] = imdiffuseest(Inoisy)

gradThresh = 1x5 uint8 row vector

   64   50   39   34   29

numIter = 5

Фильтрация шумного изображения с помощью анизотропной диффузии с оцененными параметрами.

Idiffuseest = imdiffusefilt(Inoisy,'GradientThreshold', ...
     gradThresh,'NumberOfIterations',numIter);

Для сравнения также фильтруйте шумное изображение с помощью анизотропной диффузии с параметрами по умолчанию. Порог градиента по умолчанию является 25,5, потому что тип данных изображения uint8, и количество итераций по умолчанию равно 5.

Idiffusedef = imdiffusefilt(Inoisy);

Визуально сравните два отфильтрованных изображения.

 montage({Idiffusedef,Idiffuseest},'ThumbnailSize',[])
 title(['Anisotropic Diffusion Filtering Using ' ...
     'Default Parameters (Left) vs. Estimated Parameters (Right)'])

На изображении, которое было отфильтровано с использованием параметров по умолчанию, остается некоторое количество шума. Шум практически полностью отсутствует на изображении, которое было отфильтровано с помощью оценочных параметров. Сохранена резкость ребер на обоих изображениях, особенно высококонтрастных, таких как шпалера и белый воротник.