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

ref = imread('pout.tif');
A = imnoise(ref,'salt & pepper', 0.02);

Вычислите PSNR.

[peaksnr, snr] = psnr(A, ref);
  
fprintf('\n The Peak-SNR value is %0.4f', peaksnr);

 The Peak-SNR value is 22.6437

fprintf('\n The SNR value is %0.4f \n', snr);

 The SNR value is 15.5524 

Прочтите изображение в рабочую область, затем создайте неформатированный dlarray объект с данными изображения.

ref = imread("strawberries.jpg");
ref = im2single(ref);
dlref = dlarray(ref);

Добавьте соль и перец шум к изображению, затем создайте неформатированный dlarray объект с данными зашумленного изображения.

noisy = imnoise(ref,'salt & pepper');
dlnoisy = dlarray(noisy);

Вычислите пиковые ОСШ и ОСШ зашумленных данных относительно исходных данных.

[peaksnr,snr] = psnr(dlnoisy,dlref)

peaksnr = 
  1x1 single dlarray

   17.5941

snr = 
  1x1 single dlarray

   11.1265

Чтение ссылки изображения в рабочую область.

ref = imread("office_1.jpg");

Предварительно выделите два массива, которые хранят последовательность из шести изображений размера опорного изображения.

numFrames = 6;
imsOriginal = zeros([size(ref) numFrames],class(ref));
imsNoisy = zeros([size(ref) numFrames],class(ref));

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

for p = 1:numFrames
    filename = strcat("office_",num2str(p),".jpg");
    im = imread(filename);
    imsOriginal(:,:,:,p) = im;
    imsNoisy(:,:,:,p) = imnoise(im,"salt & pepper");
end

Отобразите последовательности изображений в монтаже. Первая строка показывает последовательность с оригинальным изображением данными. Вторая строка показывает последовательность с данными шумного изображения.

montage(cat(4,imsOriginal,imsNoisy),"Size",[2 numFrames])

Вычислите пик PSNR каждого шумного изображения относительно соответствующего первозданного изображения путем определения формата данных входных решеток как 'SSCB' (пространственный, пространственный, канал, пакет).

peak_psnrs = psnr(imsNoisy,imsOriginal,"DataFormat","SSCB");
peak_psnrs = squeeze(peak_psnrs)

peak_psnrs = 6×1

   16.3560
   16.9698
   17.8079
   18.1843
   18.0656
   17.1682