Чтение изображения и создание копии с добавленным шумом. Исходное изображение является ссылочным.

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);

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

[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