南京邮电大学图像实验

通信与信息工程学院2016/2017学年第一学期实验报告实验课程名称数字图像处理与图像通信实验专业电子信息工程学生学号学生姓名指导教师谢世鹏指导单位广播工程电视系2实验一图像信号的数字化一、实验目的通过本实验了解图像的数字化参数取样频率(像素个数)、量化级数与图像质量的关系。二、实验内容1、编写并调试图像数字化程序。要求参数k，N可调。其中，k为亚抽样比例，N为量化比特数。2、可选任意图像进行处理，在显示器上观察各种数字化参数组合下的图像效果。三、实验方法及编程在数字系统中进行处理、传输和存储图像，必须把代表图像的连续信号转变为离散信号，这种变换过程称为图像信号的数字化。它包括采样和量化两种处理。本实验对数字图像进行再采样和再量化，以考察人眼对数字图像的分辨率和灰度级的敏感程度。程序的主体部分如下：functionnewbuf=Sample_Quant(oldbuf,k,n)%%函数名称：Sample_Quant()图像数字化算法函数%oldbuf原图像数组%MN原图像尺寸%k取样间隔%n量化比特值%newbuf存放处理后的图像二维数组[M,N]=size(oldbuf);oldbuf=double(oldbuf);x=1;y=1;whilexMwhileyNnew=(round(oldbuf(x,y)/(2^(8-n))))*(2^(8-n));fori=0:k-13forj=0:k-1if(x+iM)&(y+jN)newbuf(x+i,y+j)=uint8(new);endendendy=y+k;endify=Ny=1;endx=x+k;end四、实验结果及分析实验截图如下：由实验结果可以看出，亚抽样比例k和量化比特数N对都会使图像变得模糊，但两者的影响是不相同的。亚抽样比例k的大小决定了数字化图像的方块效应是否明显。当k较大时，数字化图像会有较为明显的块状出现，对于图像的视觉效果影响很大。量化比特数N则决定了图像的灰度级，量化比特数为N时，图像有个2N灰度级。所以当N较小时，图像会出现不规则的区域有着相同的灰度值的情况，但是这些区域在原图像中却有着差别较小的不同的灰度值。特别是在原图的灰度渐变的区域，这种效应会变得尤为明显。4实验二图像灰度级修正一、实验目的掌握常用的图像灰度级修正方法，即图像的灰度变换法和直方图均衡化法，加深对灰度直方图的理解。二、实验内容1、编程实现图像的灰度变换。改变图像输入、输出映射的灰度参数范围(拉伸和反比)，观看图像处理结果。2、修改可选参数gamma值，使其大于1，等于1和小于1，观看图像处理结果。3、对图像直方图均衡化处理，显示均衡前后的直方图和图像。三、实验方法及编程灰度级修正是图像空间域上的增强方法图像增强，它以图像的灰度映射变换为基础，直接对图像中的像素进行处理。常用到三种基本方式：线性，对数，幂次。程序的主体部分如下：Functionnewbuf=GrayGamma(oldbuf,r)newbuf=imadjust(oldbuf,stretchlim(oldbuf),[],r);newbuf=uint8(newbuf);functionnewbuf=GrayEqualize(oldbuf)%函数名称：%GrayEqualize()直方图均衡算法函数%参数：%oldbuf原图像数组%newbuf存放处理后的图像二维数组[M,N]=size(oldbuf);NN=M*N;sk=0;[COUNTS,X]=imhist(oldbuf,256);fori=1:Mforj=1:Nkk=double(oldbuf(i,j));fork=1:kksk=sk+COUNTS(k);5endsk=sk/NN*256;newbuf(i,j)=sk;endendnewbuf=uint8(newbuf);四、实验结果及分析实验截图如下：由实验结果可以看出，所选的测试图像的灰度值主要分布于低值的部分。经过无gamma值的灰度值变换后，直观地看出：分布于低值部分的直方图分散开来，分布于几乎所有的灰度值，但是对映于各个灰度值的像素个数分布仍是不均匀的。变换后的图像比原图显得明亮、清晰。经过第二第三幅图像可以进一步看出：gamma值的灰度值变换则明显受gamma值的影响：当gamma值大于1时，直方图有向灰度为0的一端压缩的趋势，gamma越大，这种趋势越明显。此时的图像比原图清晰，但是原本偏暗的部分更加偏黑，原本较亮的部分则变得发白，总体而言，图像偏暗的部分较多。当gamma值小于1时，其趋势与gamma值大于1相反，故整个图像显得发白。直方图均衡后的图像也显得较为清晰、均匀。它能尽量将直方图变得均衡，分布也更为均匀，各个灰度值所对应的像素个数尽可能相同。6实验三图像的平滑滤波一、实验目的图像平滑主要目的是减少噪声对图像的影响。噪声有很多种类，不同的噪声有不同的抑制措施。本实验要求用平滑线性滤波和中值滤波两种最典型、最常用的处理算法进行程序设计，学习如何对已被污染的图像进行“净化”。通过平滑处理，对结果图像加以比较，得出自己的实验结论。二、实验内容1、编写并调试窗口尺寸为m×n的平滑滤波函数。2、编写并调试窗口尺寸为m×n的中值滤波函数。三、实验方法及编程在M*N的图像f上，用m*n大小的滤波器模板进行线性滤波由这个公式给出：),(),(,nmhnymxfyxgbbnaam程序如下：functionnewbuf=AverageFilter(oldbuf,M,N,m)%函数名称：%AverageFilter()均值滤波算法函数%参数：%oldbuf噪声图像数组%MN噪声图像尺寸%m矩形平滑窗口尺寸%newbuf存放处理后的图像二维数组oldbuf=double(oldbuf);newbuf=zeros(M,N);fori=(m+1)/2:M-(m+1)/2forj=(m+1)/2:N-(m+1)/2forx=-(m-1)/2:(m-1)/2;fory=-(m-1)/2:(m-1)/2;7newbuf(i,j)=newbuf(i,j)+oldbuf(i+x,j+y)/(m*m);endendendendfunctionnewbuf=MedianFilter(oldbuf,M,N,m)%函数名称：%MedianFilter()中值滤波算法函数%参数：%oldbuf原图像数组%MN原图像尺度%m滑动窗口尺寸%newbuf存放处理后的图像数组fori=(m+1)/2:M-(m+1)/2forj=(m+1)/2:N-(m+1)/2k=1;forx=-(m-1)/2:(m-1)/2;fory=-(m-1)/2:(m-1)/2;winbuf(k)=oldbuf(i+x,j+y);k=k+1;endendnewbuf(i,j)=SeekMid(winbuf,m);endendwinsize=length(winbuf);%取窗口尺寸fori=1:winsize%编写排序函数forj=1:winsize-iif(winbuf(i)winbuf(j+i))8t=winbuf(i);winbuf(i)=winbuf(j+i);winbuf(j+i)=t;endendendmid=winbuf(ceil(m*m/2));四、实验结果及分析实验截图如下：由实验结果可以看出，对于处理椒盐噪声，中值滤波比均值滤波要好很多。当窗口大小为3×3时，均值滤波能一定程度上降低噪声的污染，但是，仍能看出在噪声点，并没有完全消除噪声，只是将噪声与周围的图像进行了平均，噪声点只是显得模糊了，并没有完全消除。相比之下，中值滤波则效果明显，同为3×3的窗口，几乎能把所有的噪声点消除，效果非常显著。当窗口大小变大时，也能将大部分噪声点消除，但是会带来严重的模糊。同时我们可以看出，两种方法都使得处理后的图像比原图要模糊，并且随着窗口的变大，图像变得越来越模糊。通过两种方法之间的比较则可以看出，同等窗口下中值滤波带来的模糊比均值滤波要轻一些。9实验四名称：图像的锐化处理一、实验目的学习用锐化处理技术来加强图像的目标边界和图像细节。对图像进行梯度算子、Roberts算子、Sobel算子边缘检测处理和Laplace算子边缘增强处理，是图像的某些特征(如边缘、轮廓等)得以进一步的增强及突出。二、实验内容1、编写梯度算子和Roberts算子滤波函数。2、编写Sobel算子滤波函数。3、编写拉普拉斯边缘增强滤波函数。三、实验方法及编程在实验中，我们对于这三种算子的算法，都有其各自对应的模板，根据这个模板我们可以方便地编写出程序的主体结构。指导书中也有对应的流程框图可供参考。程序的主体部分如下：functionnewbuf=RobFilter(oldbuf,M,N);%RobFilter()fori=1:M-1forj=1:N-1newbuf(i,j)=abs(oldbuf(i,j)-oldbuf(i+1,j+1))+abs(oldbuf(i+1,j)-oldbuf(i,j+1));endendfunctionnewbuf=SobFilter(oldbuf,M,N);%SobFilter()fori=2:M-1forj=2:N-1sx=oldbuf(i+1,j-1)+2*oldbuf(i+1,j)+oldbuf(i+1,j+1)-oldbuf(i-1,j-1)-2*oldbuf(i-1,j)-oldbuf(i-1,j+1;10sy=oldbuf(i-1,j+1)+2*oldbuf(i,j+1)+oldbuf(i+1,j+1)-oldbuf(i-1,j-1)-2*oldbuf(i,j-1)-oldbuf(i+1,j-1;newbuf(i,j)=abs(sx)+abs(sy);endendfunctionnewbuf=LapFilter(oldbuf,M,N);%LapFilter()fori=2:M-1forj=2:N-1newbuf(i,j)=5*oldbuf(i,j)-oldbuf(i-1,j)-oldbuf(i+1,j)-oldbuf(i,j-1)-oldbuf(i,j+1);endend四、实验结果及分析实验截图如下：从上面的图像可以看出：Robert梯度算子得出的图像能够得出原图的大部分边缘细节，灰度差别越11大的地方结果越大，所以显示时较为明亮。一些边缘由于灰度差值较小，在得出的结果图像中不容易分辨出来。Sobel算子得出的图像则显得明亮而粗壮。所有的边缘细节均被显示出来，特别是人物面部。由于其结果粗壮，面部细节显得非常密集。Laplace算子则用以将图像的边缘、细节增强，通过结果结果可以看出，图像的细节明显比原来突出。但是这个方法存在的弊端是，在背景区域，结果图像中有一些噪声的图样也被加强了。12实验五图像方块编码(BTC)一、实验目的通过编程实验，掌握方块编码的基本方法及压缩性能。二、实验内容编程实现子块为n×n的方块编码基本方法，分别取n=2,4,8方块尺寸进行方块编码实验，计算编码后的均方误差和压缩比。三、实验方法及编程本实验采用的方案为：（1）miiTxmxx11（2）qmiilxqmxa101Tixx（3）mqmiihxqxa111Tixx程序的主体部分如下：functionoutbuf=BtcBlock(inbuf,n)%btcblock()方块编码算法函数%inbuf方块数组%n方块尺寸%outbuf存放处理后的方块图像temp=0;temp0=0;temp1=0;q=0;m=n*n;inbuf=double(inbuf);fori=1:nforj=1:ntemp=temp+inbuf(i,j);endendxt=temp/m;13fori=1:nforj=1:nif(inbuf(i,j)=xt)q=q+1;temp1=temp1+inbuf(i,j);elsetemp0=temp0+inbuf(i,j);endendendifq~=ma0=round(temp0/(