Hi David,
Different deblurring algorithms estimate and remove blur based on how much knowledge you have of the PSF and noise in the image. You can find the available functions in MATLAB with examples in the MathWorks documentation:
Estimating the Noise-to-Signal Ratio (NSR) in an image is a challenging task, as it requires distinguishing between the noise and the underlying signal. While there is no one-size-fits-all algorithm for NSR estimation, a commonly used approach is Local Variance Estimation which can be implemented as follows:
blurred_noisy = im2double(imread('blurred_noisy.jpg')); % blurred_noisy is your blurred and noisy image (attached a sample image with the answer)
patchSize = 8; % Size of the patches
numPatches = 1000; % Number of patches to sample
% Extract random patches from the blurred image
patches = im2col(blurred_noisy, [patchSize patchSize], 'sliding');
patches = patches(:, randperm(size(patches, 2), numPatches));
% Convert patches to double precision
patches = double(patches);
% Compute the local variance of the patches
variance = var(patches, [], 1);
% Approximate the NSR as the ratio of average noise variance to average signal variance
averageNoiseVariance = mean(variance);
averageSignalVariance = mean(patches(:).^2);
estimated_nsr = averageNoiseVariance / averageSignalVariance
You can try this out for replacing the NSR estimation step with the MathWorks example documentation:
The Peak Signal-to-Noise Ratio (PSNR)can be used as a metric to evaluate the improvement in image quality after deblurring. PSNR measures the ratio of the maximum possible pixel intensity to the noise level, indicating how well the deblurred image approximates the original image. You can find more about "psnr()" function in the MathWorks documentation:
I hope this helps!
