Image processing based on MATLAB median + mean + Gaussian filtering + Laplacian + Sobel + Prewitt image denoising edge enhancement

clc
clear all
close all
img = imread('testPic.tif');
img=img(1:128*16.1:128*16);
imwrite(img,'Source0.tif'%% mean filter figure(1);
subplot(3.4.5)
subimage(img);
title('Original image');
Byave = filter2(fspecial('average'.3),img)/255; The % of3*3Mean filtering % spatial sharpening P1 =fspecial('sobel');
subplot(3.4.2)
add1 =imfilter(Byave,p1);
out1=add1+Byave;
subimage(out1);
title('Sobel Sharpening result (Mean filtering)');
[SNR,MSE,PSNR]=getSNR(img,out1);
fprintf('Mean Filtering + Sobel Operator Sharpening :\nSNR:%f\nMSE:%f\nPSNR:%f\n',SNR,MSE,PSNR)
subplot(3.4.3)
p2 =fspecial('prewitt');
add2 =imfilter(Byave,p2);
out2=add2+Byave;
subimage(out2);
title('Prewitt Sharpening results (Mean Filtering)');
[SNR,MSE,PSNR]=getSNR(img,out2);
fprintf('Mean Filtering + Prewitt Operator Sharpening :\nSNR:%f\nMSE:%f\nPSNR:%f\n',SNR,MSE,PSNR)
subplot(3.4.4)
p3 =fspecial('laplacian');
add3 =imfilter(Byave,p3);
out3=add3+Byave;
imwrite(out3,'Mean filtering + Laplacian operator sharpening. Tif');
subimage(out3);
title('Laplacian Sharpening result (mean filtering)');
[SNR,MSE,PSNR]=getSNR(img,out2);
fprintf('Mean filtering + Laplacian Sharpening :\nSNR:%f\nMSE:%f\nPSNR:%f\n'Bymed= medFilt2 (img,[3.3]); % spacial sharpening p1 = fSpecial ('sobel');
subplot(3.4.6)
add1 =imfilter(Bymed,p1);
out1=add1+Bymed;
subimage(out1);
title('Sobel sharpening result (median filtering)');
[SNR,MSE,PSNR]=getSNR(img,out1);
fprintf('Median filter + Sobel Sharpening :\nSNR:%f\nMSE:%f\nPSNR:%f\n',SNR,MSE,PSNR)
subplot(3.4.7)
p2 =fspecial('prewitt');
add2 =imfilter(Bymed,p2);
out2=add2+Bymed;
subimage(out2);
title('Prewitt Sharpening results (median filtering)');
[SNR,MSE,PSNR]=getSNR(img,out2);
fprintf('Median Filtering + Prewitt sharpening :\nSNR:%f\nMSE:%f\nPSNR:%f\n',SNR,MSE,PSNR)
subplot(3.4.8)
p3 =fspecial('laplacian');
add3 =imfilter(Bymed,p3);
out3=add3+Bymed;
imwrite(out3,'Median filtering + Laplacian operator sharpening. Tif');
subimage(out3);
title('Laplacian sharpening result (median filtering)');
[SNR,MSE,PSNR]=getSNR(img,out3);
fprintf('Median filtering + Laplacian operator sharpening :\nSNR:%f\nMSE:%f\nPSNR:%f\n',SNR,MSE,PSNR)'gaussian'[3.3].1); 
Bygaussian = imfilter(img, W, 'replicate');
 p1 =fspecial('sobel');
subplot(3.4.10)
add1 =imfilter(Bygaussian,p1);
out1=add1+Bygaussian;
subimage(out1);
Copy the code

Version: 2014 a