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

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