数字图像处理实习作业1

1、数字图像处理作业 报告编号： 01 课程编号：姓 名：赵十璧 起始日期：2012-11-01 截止日期：2012-11-07 1、 基本原理本次作业采用的计算机工具为matlab。均方误差（mse），是衡量“平均误差”的一种较方便的方法，mse可以评价数据的变化程度，mse的值越小，说明预测模型描述实验数据具有更好的精确度。计算公式为：sum=sum+(a(i,j)-b(i,j)2; mse=sum/(m*n) 。信噪比（snr），等于信号与噪声的功率谱之比 ，信噪比越高表明它产生的噪声越少。计算公式为：sum=sum+a(i,j)2; snr=10*log10(sum2/mse) 。峰值信噪

2、比（psnr），反映的是图像信噪比变化情况的统计平均。通常在经过影像压缩之后，输出的影像通常都会有某种程度与原始影像一样。为了衡量经过处理后的影像品质，我们通常会参考psnr 值来认定某个处理程序够不够令人满意。 计算公式为：psnr=10*log10(2552/mse) 。平均绝对误差（mae）是指将各次测量的绝对误差取绝对值后所求的平均值，它更加具有代表性。 计算公式为：sum=sum+a(i,j)+b(i,j); mae=sum/(m*n) 。对图像进行缩小用到两种方法：其一直接消除像素进行缩放处理；其二先进行初步的图像平滑滤波，然后消除像素。平滑滤波用于模糊处理，可以出去细节并降低噪声

3、。对图像进行放大也用到两种方法：像素复制法和双线性插值法。像素复制方法的图像缩放的原理主要是对原来输入图像的整行或是整列像素进行简单的复制与删除,达到改变图像的大小的目的，基于像素复制的方法进行图像缩放是最近邻插值方法的一种,最近邻插值是最简单的插值方法,在这种算法中,每一个插值输出像素的值就是在输入图像中与其最临近的采样点的值。双线性插值放大算法中，目标图像中新创造的象素值，是由源图像位置在它附近的2*2区域4个邻近象素的值通过加权平均计算得出的。双线性内插值算法放大后的图像质量较高，不会出现像素值不连续的的情况。然而次算法具有低通滤波器的性质，使高频分量受损，所以可能会使图像轮廓在一定程度

4、上变得模糊。灰度图像的量化和直方图均衡化直接调用函数。量化是指将像素灰度转换成离散的整数值的过程。量化等级越多，所得图像层次越丰富，灰度分辨率高，图像质量好，但数据量大；量化等级越少，图像层次欠丰富，灰度分辨率低，图像质量差。直方图均衡化是指将原图像通过某种变换，得到一幅灰度直方图为均匀分布的新图像的方法。2、 操作结果task1中，所求的mse为640.6602，snr 为68.5934，psnr为 20.0645，mae为 248.1302。task2中，用直接消除像素点得到的图像如下：图2-2 滤波后的原图图2-1 原始图像图2-6 平滑滤波后缩小4倍图2-5 直接消除像素法缩小4倍图2

5、-4 平滑滤波后缩小2倍图2-2 直接消除像素法缩小2倍 图2-8 平滑滤波后缩小8倍图2-7 直接消除像素法缩小8倍 图2-10 平滑滤波后缩小16倍 图2-9 直接消除像素法缩小16倍 通过图像之间的对比可以发现，不管用哪种方法，当放大的倍数越大时，其清晰度就越低。另外，虽然原始大小图像、2倍比例及16倍比例缩小下两种缩小方法得到的图像对比并不明显，但是通过4倍和8倍缩小比例下图像的对比可以发现，经过平滑滤波处理后的图像在缩小后比原始图像缩小后得到的图像更为清晰。 task3中，分别采用了像素重复和双线性插值两种方法放大在上一个任务中缩小的图像，并分别计算原图像与放大后的图像的均方误差（m

6、se）和峰值信噪比（psnr），然后根据图像的缩小比例分别绘出mse和psnr的变化图像。变化图像如下：图3-1 mse图片缩小比例折线图 图3-2 psnr图片缩小比例折线图 由图可以看出，在1:2的缩放比例下，滤波处理缩小后用像素重复法放大得到的mse最小且psnr最大，其他比例下滤波处理缩小后用双线性插值法放大得到的mse最小且psnr最大。图4-1 原图像及其灰度分布直方图图4-2 量化步长为2的图像及其灰度分布直方图图4-3 量化步长为4的图像及其灰度分布直方图task4中，得到的图像及其灰度分布直方图如下：图4-4 量化步长为8的图像及其灰度分布直方图图4-5 量化步长为16的图像

7、及其灰度分布直方图图4-6 量化步长为32的图像及其灰度分布直方图图4-7 量化步长为64的图像及其灰度分布直方图图4-8 量化步长为128的图像及其灰度分布直方图 由图可以看出，当量化步长越大时，图像显得越细腻，其灰度分布直方图越密集。 图4-9 mse在不同量化步长下的值 由图4-9可以看出，当量化步长越大时，mse越小，及说明图像描述的精确度越高。图5-1 fig0316(1)(top_left)原图像及其灰度分布直方图图5-2 fig0316(1)(top_left)均衡化后的图像及其灰度分布直方图图5-3 fig0316(4)(bottom_left)原图像及其灰度分布直方图 tas

8、k5中，得到的各图像如下:图5-2 fig0316(4)(bottom_left)均衡化后的图像及其灰度分布直方图 我们可以发现，均衡化之后图像的对比度更加明显，兵器其灰度分布直方图分布区域更加均匀。3、 附源代码task1 f=imread(lena_gray_512.tif); g=imnoise(f,gaussian); a=double(f); b=double(g); m,n=size(f); sum=0; for i=1:m;for j=1:n;sum=sum+(a(i,j)-b(i,j)*(a(i,j)-b(i,j);end;end; mse=sum/(m*n)mse = 640

9、.6602 sum1=0; for i=1:m;for j=1:n;sum1=sum1+a(i,j)*a(i,j);end;end; snr=10*log10(sum1/mse)snr = 68.5934 psnr=10*log10(2552/mse)psnr = 20.0645 sum=0; for i=1:m;for j=1:n;sum=sum+a(i,j)+b(i,j);end;end; mae=sum/(m*n)mae = 248.1302 task2 f=imread(livingroom.tif); imshow(f); f1=f(1:2:end,1:2:end); imshow(

10、f1); f2=f(1:4:end,1:4:end); imshow(f2); f3=f(1:8:end,1;8:end); imshow(f3); f4=f(1:16:end,1:16:end); imshow(f4); w=1 1 1;1 1 1;1 1 1/9; g=imfilter(f,w,replicate); imshow(g); w=ones(31); g=imfilter(f,w,replicate); imshow(g); w=1 1 1;1 1 1;1 1 1/9; g=imfilter(f,w,replicate); imshow(g); g1=g(1:2:end,1:2

11、:end); imshow(g1); g2=g(1:4:end,1:4:end); imshow(g2); g3=g(1:8:end,1:8:end); imshow(g3); g4=g(1:16:end,1:16:end); imshow(g4); task3 f=imread(lena_frag.tif); figure; g=histeq(1,2); imhist(f); clear; f=imread(lena_frag.tif); g=histeq(f,2); figure; subplot(2,2,1); imshow(f); subplot(2,2,2); imhist(f);

12、subplot(2,2,3); imshow(g); subplot(2,2,4); imhist(g); a=double(f); g=double(g); m,n=size(f); sum=0; for i=1:m;for j=1:n;sum=sum+(a(i,j)-b(i,j)*(a(i,j)-b(i,j);end;end;undefined function b for input arguments of type double. b=double(g); m,n=size(f); sum=0; for i=1:m;for j=1:n;sum=sum+(a(i,j)-b(i,j)*(

13、a(i,j)-b(i,j);end;end; mse=sum/(m*n)mse = 9.6389e+03 g=histeq(f,4); imshow(f); figure; f=imread(lena_frag.tif); g=histeq(f,2); figure; subplot(2,2,1); imshow(f); subplot(2,2,2); imhist(f); subplot(2,2,3); imshow(g); subplot(2,2,4); imhist(g); g=histeq(f,4); figure; subplot(2,2,1); imshow(f); subplot

14、(2,2,2); imhist(f); f=imread(lena_frag.tif); subplot(1,2,1); imshow(f); subplot(1,2,2); imhist(f); g1=histeq(f,2); figure; subplot(1,2,1); imshow(g1); subplot(1,2,2); imhist(g1); g2=histeq(f,4); figure; imshow(g2); subplot(1,2,2); subplot(1,2,1); imshow(g1); imshow(g2); subplot(1,2,2); imhist(g2); g

15、3=histeq(f,8); subplot(1,2,1); imshow(g3); subplot(1,2,2); imhist(g3); g4=histeq(f,16); figure; subplot(1,2,1); imshow(g4); subplot(1,2,2); imhist(g4); g5=histeq(f,32); figure; subplot(1,2,1); imshow(g5); subplot(1,2,2); imhist(g5); g6=histeq(f,64); figure; subplot(1,2,1); imshow(g6); subplot(1,2,2)

16、; imhist(g6); g7=histeq(f,128); figure; subplot(1,2,1); imshow(g7); subplot(1,2,2); imhist(g7); m,n=size(f); a=double(f); b1=double(g1); b1=double(g1); b1=double(g1); b2=double(g2); b3=double(g3); b4=double(g4); b5=double(g5); b6=double(g6); b7=double(g7); for i=1:m;for j=1:n;end;end; sum1=0; sum2=0

17、; sum3=0sum3 = 0 sum3=0; sum4=0; sum5=0; sum6=0; sum7=0; for i=1:m;for j=1:n;sum1=sum1+(a(i,j)-b1(i,j)*(a(i,j)-b1(i,j);sum2=sum2+(a(i,j)-b2(i,j)*(a(i,j)-b2(i,j);sum3=sum3+(a(i,j)-b3(i,j)*(a(i,j)-b3(i,j);sum4=sum4+(a(i,j)-b4(i,j)*(a(i,j)-b4(i,j);sum5=sum5+(a(i,j)-b5(i,j)*(a(i,j)-b5(i,j);sum6=sum6+(a(

18、i,j)-b6(i,j)*(a(i,j)-b6(i,j);sum7=sum7+(a(i,j)-b7(i,j)*(a(i,j)-b7(i,j);end;end; mse1=sum1/(m*n)mse1 = 9.6389e+03 mse2=sum2/(m*n)mse2 = 3.5931e+03 mse3=sum3/(m*n)mse3 = 2.4078e+03 mse4=sum4/(m*n)mse4 = 2.0164e+03 mse5=sum5/(m*n)mse5 = 1.8854e+03 mse6=sum6/(m*n)mse6 = 1.8229e+03 mse7=sum7/(m*n)mse7 =

19、1.8002e+03 clear f1=imread(fig0316(1)(top_left).tif); g=im2double(f1); subplot(2,2,1); imshow(g); subplot(2,2,2); imhist(g); g1=histeq(g); subplot(2,2,3); imshow(g1); subplot(2,2,4); imhist(g1); f2=imread() f2=imread(fig0316(4)(bottom_left);error using imread (line 369)file fig0316(4)(bottom_left) d

20、oes not exist. f2=imread(fig0316(4)(bottom_left);error using imread (line 369)file fig0316(4)(bottom_left) does not exist. f2=imread(fig0316(4)(bottom_left).tif); figure; h=im2double(f2); subplot(2,2,1); imshow(h); subplot(2,2,2); subplot(2,2,1); imshow(h); subplot(2,2,1); subplot(2,2,2); imhist(h);

21、 subplot(2,2,3); h1=histeq(h); subplot(2,2,3); imshow(h1); subplot(2,2,4); imhist(h1); clear f=imread(livingroom.tif); f1=f(1:2:end,1:2:end); f2=f(1:4:end,1:4:end); f3=f(1:8:end,1:8:end); f4=f(1:16:end,1:16:end); w=1 1 1;1 1 1;1 1 1/9; g=imfilter(f,w,replicate); g1=g(1:2:end,1:2:end); g2=g(1:4:end,1

22、:4:end); g3=g(1:8:end,1:8:end); g4=g(1:16:end,1:16:end); h1=imresize(g1,2); h2=imresize(g2,2); h2=imresize(g2,4); h3=imresize(g3,8); h4=imresize(g4,16); k1=imresize(f1,2,bilinear); k2=imresize(f2,4,bilinear); k3=imresize(f3,8,bilinear); k4=imresize(f4,16,bilinear); k1=imresize(f1,2); k2=imresize(f2,

23、4); k3=imresize(f3,8); k4=imresize(f4,16); m1=imresize(f1,2,bilinear); m2=imresize(f2,4,bilinear); m3=imresize(f3,8,bilinear); m4=imresize(f4,16,bilinear); n1=imresize(g1,2,bilinear); n2=imresize(g2,4,bilinear); n3=imresize(g3,8,bilinear); n4=imresize(g4,16,bilinear); a=double(f); m,n=size(f); b1=do

24、uble(h1); b2=double(k1); b3=double(m1); b4=double(n1); sum1=0; sum2=0; sum3=0; sum4=0; for i=1:m;for j=1:n;sum1=sum1+(a(i,j)-b1(i,j)*(a(i,j)-b1(i,j);sum2=sum2+(a(i,j)-b2(i,j)*(a(i,j)-b2(i,j);sum3=sum3+(a(i,j)-b3(i,j)*(a(i,j)-b3(i,j);sum4=sum4+(a(i,j)-b4(i,j)*(a(i,j)-b4(i,j);end;end; mse1=sum1/(m*n)m

25、se1 = 123.4420 mse2=sum2/(m*n)mse2 = 138.0037 mse3=sum3/(m*n)mse3 = 132.1592 mse4=sum4/(m*n)mse4 = 137.4242 psnr1=10*log10(2552/mse1)psnr1 = 27.2162 psnr2=10*log10(2552/mse2)psnr2 = 26.7319 psnr3=10*log10(2552/mse3)psnr3 = 26.9198 psnr4=10*log10(2552/mse4)psnr4 = 26.7502 b5=double(h2); b6=double(k2)

26、; b7=double(m1); b7=double(m2); b8=double(n2); b9=double(h3); b10=double(k3); b11=double(m3); b12=double(n3); b13=double(h4); b14=double(k4); b15=double(m4); b16=double(n4); sum5=0; sum6=0; sum7=0; sum8=0; sum9=0; sum10=0; sum11=0; sum12=0; sum13=0; sum14=0; sum15=0; sum16=0; for i=1:m;for j=1:n;sum

27、5=sum5+(a(i,j)-b5(i,j)*(a(i,j)-b5(i,j);sum6=sum6+(a(i,j)-b6(i,j)*(a(i,j)-b6(i,j);sum7=sum7+(a(i,j)-b7(i,j)*(a(i,j)-b7(i,j);sum8=sum8+(a(i,j)-b8(i,j)*(a(i,j)-b8(i,j);sum9=sum9+(a(i,j)-b9(i,j)*(a(i,j)-b9(i,j);sum10=sum10+(a(i,j)-b10(i,j)*(a(i,j)-b10(i,j);sum11=sum11+(a(i,j)-b11(i,j)*(a(i,j)-b11(i,j);s

28、um12=sum12+(a(i,j)-b12(i,j)*(a(i,j)-b12(i,j);sum13=sum13+(a(i,j)-b13(i,j)*(a(i,j)-b13(i,j);sum14=sum14+(a(i,j)-b14(i,j)*(a(i,j)-b14(i,j);sum15=sum15+(a(i,j)-b15(i,j)*(a(i,j)-b15(i,j);sum16=sum16+(a(i,j)-b16(i,j)*(a(i,j)-b16(i,j);end;end; mse5=sum5/(m*n)mse5 = 332.4135 mse6=sum6/(m*n)mse6 = 424.9083

29、mse7=sum7/(m*n)mse7 = 387.9766 mse8=sum8/(m*n)mse8 = 324.0842 mse9=sum9/(m*n)mse9 = 661.7173 mse10=sum10/(m*n)mse10 = 770.8898 mse11=sum11/(m*n)mse11 = 694.3029 mse12=sum12/(m*n)mse12 = 618.3956 mse13=sum13/(m*n)mse13 = 1.0948e+03 mse14=sum14/(m*n)mse14 = 1.2116e+03 mse15=sum15/(m*n)mse15 = 1.0934e+

30、03 mse16=sum16/(m*n)mse16 = 1.0098e+03 psnr5=10*log10(2552/mse5)psnr5 = 22.9140 psnr6=10*log10(2552/mse6)psnr6 = 21.8479 psnr7=10*log10(2552/mse7)psnr7 = 22.2427 psnr8=10*log10(2552/mse8)psnr8 = 23.0242 psnr9=10*log10(2552/mse9)psnr9 = 19.9241 psnr10=10*log10(2552/mse10)psnr10 = 19.2609 psnr11=10*lo

31、g10(2552/mse11)psnr11 = 19.7153 psnr12=10*log10(2552/mse12)psnr12 = 20.2181 psnr13=10*log10(2552/mse13)psnr13 = 17.7376 psnr14=10*log10(2552/mse14)psnr14 = 17.2972 psnr15=10*log10(2552/mse15)psnr15 = 17.7428 psnr16=10*log10(2552/mse16)psnr16 = 18.0883 task4 f=imread(lena_frag.tif); subplot(1,2,1); i

32、mshow(f); subplot(1,2,2); imhist(f); g1=histeq(f,2); figure; subplot(1,2,1); imshow(g1); subplot(1,2,2); imhist(g1); g2=histeq(f,4); figure; imshow(g2); subplot(1,2,2); subplot(1,2,1); imshow(g1); imshow(g2); subplot(1,2,2); imhist(g2); g3=histeq(f,8);figure; subplot(1,2,1); imshow(g3); subplot(1,2,

33、2); imhist(g3); g4=histeq(f,16); figure; subplot(1,2,1); imshow(g4); subplot(1,2,2); imhist(g4); g5=histeq(f,32); figure; subplot(1,2,1); imshow(g5); subplot(1,2,2); imhist(g5); g6=histeq(f,64); figure; subplot(1,2,1); imshow(g6); subplot(1,2,2); imhist(g6); g7=histeq(f,128); figure; subplot(1,2,1); imshow(g7); subplot(1,2,2); imhist(g7); m,n=size(f); a=double(f); b1=double