图像质量评价：合成失真图像方法
1、高斯噪声
代码
%%高斯噪声
a=imread('monarch.bmp');
subplot(2,3,1);
imshow(a);
title ('原始图像');
Inoise1=imnoise(a,'gaussian',0.1,0.01);%对图像加入高斯噪声1
subplot(2,3,2);
imshow(Inoise1);
title('加入高斯噪声1');
Inoise2=imnoise(a,'gaussian',0.2,0.01);%对图像加入高斯噪声2
subplot(2,3,3);
imshow(Inoise2);
title('加入高斯噪声2');
Inoise3=imnoise(a,'gaussian',0.1,0.05);%对图像加入高斯噪声4
subplot(2,3,4);
imshow(Inoise3);
title('加入高斯噪声3');
Inoise4=imnoise(a,'gaussian',0.2,0.09);%对图像加入高斯噪声5
subplot(2,3,5);
imshow(Inoise4);
title('加入高斯噪声4');
Inoise5=imnoise(a,'gaussian',0.5,0.01);%对图像加入高斯噪声3
subplot(2,3,6);
imshow(Inoise5);
title('加入高斯噪声6');
2、运动模糊
代码
%%运动模糊
a=imread('monarch.bmp');
subplot(2,3,1);
imshow(a);
title ('原始图像');
PSF1=fspecial('motion',10,10);
blurred=imfilter(a,PSF1,'circular','conv');
subplot(2,3,2),imshow(blurred),title('运动模糊1')
PSF2=fspecial('motion',20,5);
blurred=imfilter(a,PSF2,'circular','conv');
subplot(2,3,3),imshow(blurred),title('运动模糊2')
PSF3=fspecial('motion',30,20);
blurred=imfilter(a,PSF3,'circular','conv');
subplot(2,3,4),imshow(blurred),title('运动模糊3')
PSF4=fspecial('motion',30,50);
blurred=imfilter(a,PSF4,'circular','conv');
subplot(2,3,5),imshow(blurred),title('运动模糊4')
PSF5=fspecial('motion',40,50);
blurred=imfilter(a,PSF5,'circular','conv');
subplot(2,3,6),imshow(blurred),title('运动模糊5')
3、高斯模糊
代码
a=imread('monarch.bmp');
subplot(2,3,1);
imshow(a);
title ('原始图像');
G1=fspecial('gaussian',10,10);
gaus1=imfilter(a,G1);
subplot(2,3,2),imshow(gaus1),title('高斯模糊1')
G2=fspecial('gaussian',20,10);
gaus2=imfilter(a,G2);
subplot(2,3,4),imshow(gaus2),title('高斯模糊3')
G3=fspecial('gaussian',10,20);
gaus3=imfilter(a,G3);
subplot(2,3,3),imshow(gaus3),title('高斯模糊2')
G4=fspecial('gaussian',30,10);
gaus4=imfilter(a,G4);
subplot(2,3,5),imshow(gaus4),title('高斯模糊4')
G5=fspecial('gaussian',30,20);
gaus5=imfilter(a,G5);
subplot(2,3,6),imshow(gaus5),title('高斯模糊5')
4、椒盐噪声
代码
%%椒盐噪声
a=imread('monarch.bmp');
subplot(2,3,1);
imshow(a);
title ('原始图像');
J1=imnoise(a,'salt & pepper',0.01);
subplot(2,3,2),imshow(J1),title('椒盐噪声1');
J2=imnoise(a,'salt & pepper',0.02);
subplot(2,3,3),imshow(J2),title('椒盐噪声2');
J3=imnoise(a,'salt & pepper',0.05);
subplot(2,3,4),imshow(J3),title('椒盐噪声3');
J4=imnoise(a,'salt & pepper',0.1);
subplot(2,3,5),imshow(J4),title('椒盐噪声4');
J5=imnoise(a,'salt & pepper',0.3);
subplot(2,3,6),imshow(J5),title('椒盐噪声5');
5、均值滤波
代码
%%平均值滤波
a=imread('monarch.bmp');
subplot(2,3,1);
imshow(a);
title ('原始图像');
h1 = fspecial('average',5);
hblur1=imfilter(a,h1);
subplot(2,3,2),imshow(hblur1),title('平均值滤波1')
h2 = fspecial('average',8);
hblur2=imfilter(a,h2);
subplot(2,3,3),imshow(hblur2),title('平均值滤波2')
h3 = fspecial('average',10);
hblur3=imfilter(a,h3);
subplot(2,3,4),imshow(hblur3),title('平均值滤波3')
h4 = fspecial('average',15);
hblur4=imfilter(a,h4);
subplot(2,3,5),imshow(hblur4),title('平均值滤波4')
h5 = fspecial('average',20);
hblur5=imfilter(a,h5);
subplot(2,3,6),imshow(hblur5),title('平均值滤波5')
6、圆形滤波
代码
%%圆形滤波
a=imread('monarch.bmp');
subplot(2,3,1);
imshow(a);
title ('原始图像');
blur1=fspecial('disk',5);
averageblur1=imfilter(a,blur1,'replicate');
subplot(2,3,2),imshow(averageblur1),title('圆形滤波器模糊1')
blur2=fspecial('disk',8);
averageblur2=imfilter(a,blur2,'replicate');
subplot(2,3,3),imshow(averageblur2),title('圆形滤波器模糊2')
blur3=fspecial('disk',10);
averageblur3=imfilter(a,blur3,'replicate');
subplot(2,3,4),imshow(averageblur3),title('圆形滤波器模糊3')
blur4=fspecial('disk',15);
averageblur4=imfilter(a,blur4,'replicate');
subplot(2,3,5),imshow(averageblur4),title('圆形滤波器模糊4')
blur5=fspecial('disk',18);
averageblur5=imfilter(a,blur5,'replicate');
subplot(2,3,6),imshow(averageblur5),title('圆形滤波器模糊5')
7、JPEG压缩
代码
%%圆形滤波
a=imread('monarch.bmp');
subplot(2,3,1);
imshow(a);
title ('原始图像');
imwrite(I,'monarch1.jpeg','Quality',20);
hblur=imread('monarch1.jpeg');
subplot(2,3,2),imshow(hblur),title('jpeg压缩1');
imwrite(I,'monarch2.jpeg','Quality',15);
hblur2=imread('monarch2.jpeg');
subplot(2,3,3),imshow(hblur2),title('jpeg压缩2')
imwrite(I,'monarch3.jpeg','Quality',10);
hblur3=imread('monarch2.jpeg');
subplot(2,3,4),imshow(hblur3),title('jpeg压缩3')
imwrite(I,'monarch4.jpeg','Quality',5);
hblur4=imread('monarch4.jpeg');
subplot(2,3,5),imshow(hblur4),title('jpeg压缩4')
imwrite(I,'monarch5.jpeg','Quality',1);
hblur5=imread('monarch5.jpeg');
subplot(2,3,6),imshow(hblur5),title('jpeg压缩5')
8、JPEG2000压缩
代码
%%圆形滤波
a=imread('monarch.bmp');
subplot(2,3,1);
imshow(a);
title ('原始图像');
imwrite(I,'monarch1.jp2','CompressionRatio',50);
hblur=imread('monarch1.jp2');
subplot(2,3,2),imshow(hblur),title('jpeg2000压缩1');
imwrite(I,'monarch2.jp2','CompressionRatio',100);
hblur2=imread('monarch2.jp2');
subplot(2,3,3),imshow(hblur2),title('jpeg2000压缩2')
imwrite(I,'monarch3.jp2','CompressionRatio',200);
hblur3=imread('monarch2.jp2');
subplot(2,3,4),imshow(hblur3),title('jpeg2000压缩3')
imwrite(I,'monarch4.jp2','CompressionRatio',400);
hblur4=imread('monarch4.jp2');
subplot(2,3,5),imshow(hblur4),title('jpeg2000压缩4')
imwrite(I,'monarch5.jp2','CompressionRatio',800);
hblur5=imread('monarch5.jp2');
subplot(2,3,6),imshow(hblur5),title('jpeg2000压缩5')
9、乘性噪声
代码
a=imread('monarch.bmp');
subplot(2,3,1);
imshow(a);
title ('原始图像');
Inoise1=imnoise(a,'speckle',0.1);
subplot(2,3,2);
imshow(Inoise1);
title('加入乘性噪声1');
Inoise2=imnoise(a,'speckle',0.2);
subplot(2,3,3);
imshow(Inoise2);
title('加入乘性噪声2');
Inoise3=imnoise(a,'speckle',0.1);
subplot(2,3,4);
imshow(Inoise3);
title('加入乘性噪声3');
Inoise4=imnoise(a,'speckle',0.2);
subplot(2,3,5);
imshow(Inoise4);
title('加入乘性噪声4');
Inoise5=imnoise(a,'speckle',0.5);
subplot(2,3,6);
imshow(Inoise5);
title('加入乘性噪声6');