LR7 / LR7_1

clc; clear; close all;

% Чтение исходного изображения
input_image = imread('image.jpg'); % путь к изображению

% Радиус расфокусировки
P = 5; % радиус сферической расфокусировки

% Создание маски расфокусировки
PSF = fspecial('disk', P);

% Применение расфокусировки
blurred_image = imfilter(input_image, PSF, 'replicate');

% Сохранение изображений
imwrite(input_image, 'original_image.jpg');
imwrite(blurred_image, 'blurred_image.jpg');

% Параметры для восстановления
noise_var = 0.01; % Дисперсия шума

% Восстановление методом Тихонова
restored_Tikhonov = deconvreg(blurred_image, PSF, [], noise_var);

% Восстановление методом Винера
restored_Wiener = deconvwnr(blurred_image, PSF, noise_var);

% Сохранение результатов
imwrite(restored_Tikhonov, 'restored_Tikhonov.jpg');
imwrite(restored_Wiener, 'restored_Wiener.jpg');