数字图像处理杨淑莹数字图象处理课程全部教案省公开课一等奖全国示范课微课金奖PPT课件.pptx

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,数字图像处理,计算机科学与工程系,杨淑莹,1/124,课程内容介绍,介绍各种数字图象处理算法分析及编程实现技术。,主要内容包含：,位图基础、图象显示、图像几何变换、图象灰度变换、图像平滑处理、图像锐化处理及边缘检测、图像分割及测量、图像形态学处理、图像变换域处理及应用、图像合成、24位彩色图像处理、JPEG图像压缩编码。,2,2/124,第一章 位图基础,1.1 数字图像基本概念,1.2 与设备无关位图(DIB),3,3/124,1.1 数字图像基本概念,位图与矢量图,数字图象基本概念,单色图象,灰度图象,伪彩色图像,24位真彩色图像,4,4/124,位图与矢量图,图形由指令集合组成描述，显示时需要对应软件读取这些命令，并将其转变为屏幕上所显示形状和颜色，图形统计主要内容是坐标值或坐标值序列，对普通画面内容颜色或亮度隐含且统一地描述。,图像是以栅格结构存放画面内容，栅格结构将一幅图划分为均匀分布栅格，每个栅格称为像素，显式地统计每一像素光度值(亮度或彩色)，全部像素位置按规则方式排列，像素位置坐标值却是有规则地隐含。,5,5/124,数字图象基本概念,在PC显示系统中，显示图像是由一个个像素组成，每一个像素都有自己颜色属性，每一个像素颜色由红，绿，蓝三原色组合而成。数字图象是连续图象f(x，y)一个近似表示，通惯用由采样点值所组成矩阵来表示:,f(0，0)f(0，1).f(0，M-1),f(1，0)f(1，1).f(1，M-1),f(N-1，0)f(N-1，1).f(N-1，M-1),每一个采样单元叫做一个象素(pixel)，上式中，M、N分别为数字图象在横、纵方向上象素。在计算机内通惯用二维数组来表示数字图象矩阵。,6,6/124,数字图象基本概念,图象数字化精度包含两个部分，即分辨率和颜色深度。,分辨率：有显示分辨率和图象分辨率。,图象分辨率：是数字化图象时划分图象象素密度，即单位长度内象素数，其单位是每英寸点数DIP（Dots per Inche）。图象分辨率说明了数字图象实际精细度。,显示分辨率：是把数字图象在输出设备（如显示器或打印机等）上能够显示象素数目和所显示象素之间点距。显示分辨率说明了数字图象表现精细度。,含有不一样图象分辨率数字图象在同一输出设备上显示分辨率相同。,颜色深度：是指表示每一象素颜色值二进制位数。颜色深度越大则能表示象素颜色数目越多，它们之间关系取决于数字图象采取颜色表示法。惯用颜色表示法有RGB、CMYK、YUV等。,7,7/124,单色图象,单色图象每一象素值含有以下特征：,1)图象中每一点红、绿、蓝颜色分量值都相等 即:,f red(x，y)=f green(x，y)=f blue(x，y),2)其值只有两个0或255。,单色图象含有比较简单格式，它普通由黑色区域和白色区域组成，其经典图象格式就是直线图和图表。一旦给定图象大小，这种格式内容将十分紧凑.,8,8/124,灰度图象,灰度图象每一象素值含有以下特征：,1)图象中每一点都不是彩色，即每一点红绿蓝颜色分量值都相等，即：,f red(x，y)=f green(x，y)=f blue(x，y),2)其值范围从0255。,对于灰度图象，其f(x，y)表示(x，y)位置处灰度值。,在一些单色照片中，经常要用灰度，才能准确表示图象真实视觉效果，这种图象被称为灰度图象。当需要在不具备真彩色支持图象显示卡上显示原始真彩色图象时，普通都要利用对应灰度图象对其进行近似处理，所以灰度图象是图象技术中包括范围比较广泛一个图象表示方法。,9,9/124,伪彩色图像,伪彩色图像与灰度图像相同，其存放文件中也带 有图像颜色表，伪彩色图像含有以下特征：,1)图像颜色表中红、绿、蓝颜色分量值不全相 等，即：,f red(x，y)f green(x，y)f blue(x，y),2)整幅图像仅有256种颜色，,要表示,256,种不一样颜色，像素必须由,8,位组成，每个像素值不是由每个基色分量数值直接决定，而是把像素值看成,图像颜色表,表项入口地址。把含有,256,色,图像,称为,8,位,彩色,图像。,10,10/124,24位真彩色图像,含有全彩色照片表示能力图像为,24,位彩色图像，,24位真彩色图像存放文件中不带有图像颜色表，其含有以下特征：,1)图像中每一像素由RGB三个分量组成，每个分量各占8位，,每个像素需,24,位。,2)f red(x，y)、f green(x，y)、f blue(x，y)取值范围为0255。,因为,24,位真彩色图像所需存放空间很大，处理速度较慢，,当需要存放空间不大，而且要求实时快速处理图像时，普通都要利用对应8位位图对其进行近似处理，所以8位位图是图像技术中包括范围比较广泛一个图像表示方法。,11,11/124,1.2与设备无关位图(DIB),BMP文件组成,数字图像处理应用,12,12/124,BMP文件组成,BMP文件由文件头、位图信息头、颜色信息和图像数据四部分组成。位图结构以下所表示:,位图颜色表RGBQUAD,位图像素数据,位图文件头结构BITMAPFILEHEADER,位图信息头结构BITMAPINFOHEADER,13,13/124,结构,数据域,含义,FILEHEADERBITMAP,bfType,字节“BM”,bfSize,文件总字节数,BfReserved1,0,BfReserved1,0,bfOffBits,位图数据距文件头偏移量,BITMAPINFOHEADER,biSize,本结构大小字节数,biWidth,位图宽度，单位为象素点,biHeight,位图高度，单位为象素点,biPlanes,1,biBitCount,象素位数（1，4，8，24）,14,14/124,biCompression,压缩方式（0表示不压缩）,bisizeImage,位图数据字节数,biXpelsPerMeter,水平分辨率（象素点/米）,biYpelsPerMeter,垂直分辨率（象素点/米）,biClrUsed,图象中使用颜色数,biClrImportant,图象中主要颜色数,RGBQUAD,rgbBlue,蓝色百分比,rgbGreen,绿色百分比,rgbRed,红色百分比,rgbReserved,0,15,15/124,数字图像处理应用,下面列举一些经典应用实例。,(1)遥感航天中应用,天文、太空星体探测及分析；,军事侦察、定位、指挥等；,地质、地形、地图普查及绘制；,地下矿藏勘探；,环境污染监测；,气象、天气预报合成份析。,(2)生物医学中应用,显微图像处理；,DNA显示分析；,生物进化图像分析；,手术规划；,内脏大小、形状、活动及异常检出分析；,癌细胞识别。,16,16/124,数字图像处理应用,(3)工业应用,产品无损检测、焊缝及内部缺点检测；,流水线零件自动检测识别；,生产过程监控；,交通管制、机场监控；,支票、署名辨伪及识别；,机器人视觉系统应用。,(4)军事公安领域中应用,罪犯脸形合成；,指纹自动识别；,巡航导弹地形识别；,手迹、印章判定识别；,遥控飞行器引导；,雷达目标侦察。,17,17/124,数字图像处理应用,(5)其它应用,多媒体计算机系统及应用；,图像远程通信,电视会议；,可视电话；,服装试穿显示；,现场视频管理。,18,18/124,第二章 图像显示,2.1 调色板基本应用,2.2 图像特效显示,19,19/124,2.1 调色板基本应用,调色板原理,调色板创建与实现,显示函数,20,20/124,调色板原理,调色板就是在256色显示系统中,将图像中出现256种颜色组成颜色表,对这些颜色按8位,即0至255进行编号,每一编号代表其中一个颜色,在这种颜色中颜色编号叫做颜色索引号。,256色位图文件中含有BMP文件由文件头、位图信息头、颜色信息表和图像数据四部分组成。图像像素值并不是颜色值,而是颜色索引表索引号。,21,21/124,调色板创建与实现,在调色板系统中,每一幅图像都有自己调色板,显示时必须将自己调色板载入系统调色板中,实现调色板。普通创建一个逻辑调色板需要以下五步：,(1)建立一个LOGPALETTE结构和PALETTEENTRY数组；,(2)对数组元素进行初始化并对组员变量进行设置；,(3)建立CPalette对象并使用CreatePalette函数初始化调 色板对象；,(4)使用SelectPalette函数来将设备描述表和调色板联络起来；,(5)使用CDC中RealizePalette函数使调色板生效。,22,22/124,显示函数,在Visual C+环境中用于显示主要有BitBlt函,数、StretchBlt函数和StretchDIBits函数。,BitBlt函数和StretchBlt函数属于CDC类，而,StretchDIBits函数为API函数。,23,23/124,显示函数,BitBlt()函数,该函数结构以下：,BOOL BitBlt(int x,int y,int nWidth,int nHeight,CDC*pSrcDC,int xSrc,int ySrc,DWORD dwrop);,参数说明：,x:指定绘制目标矩形左上角逻辑X轴位置。,y:指定绘制目标矩形左上角逻辑Y轴位置。,nWidth:指定绘制目标矩形和源位图宽度（按逻辑单位）。,nHeight:指定绘制目标矩形和源位图高度（按逻辑单位）。,pSrcDC:设备上下文指针。,xSrc,ySrc:指定源位图左上角逻辑X轴、Y轴位置。,dwrop:指定要执行光栅运算。其操作码定义了图形设备接口（GDI）在输出操作中怎样组合颜色。,当要将位图像素从内存显示环境复制到显示器（或打印机）设备环境中，普通会用到这个函数。,24,24/124,显示函数,StretchBlt（）函数,该函数结构以下：,BOOL StretchBlt（int x,int y,int nWidth,int nHeight,CDC*pSrcDC,int xSrc,int ySrc,int nSrcWidth,int nSrcHeight,DWORD dwrop);,参数说明：,x:指定绘制目标矩形左上角逻辑X轴位置。,y:指定绘制目标矩形左上角逻辑Y轴位置。,nWidth:指定绘制目标矩形宽度（按逻辑单位）。,nHeight:指定绘制目标矩形高度（按逻辑单位）。,pSrcDC:设备上下文指针,xSrc,ySrc:指定源位图左上角坐标（按逻辑单位）。,nSrcWidth,nSrcHeight:指定复制源位图宽度和高度（按逻辑单位）。,dwrop:指定要执行光栅运算。其操作码定义了图形设备接口（GDI）在输出操作中怎样组合颜色。,25,25/124,显示函数,StretchDIBits（）函数,该函数结构以下：,BOOL StretchDIBits(HDC hdc，int x,int y,int nWidth,int nHeight,int xSrc,int ySrc,int nSrcWidth,int nSrcHeight,CONST VOID*lpvBits,CONST BITMAPINFO*lpbmi,UNIT fuColorUse，DWORD dwrop);,参数说明：,hdc:设备上下文句柄。,x:指定绘制目标矩形左上角逻辑X轴位置。,y:指定绘制目标矩形左上角逻辑Y轴位置。,nWidth:指定绘制目标矩形宽度（按逻辑单位）。,nHeight:指定绘制目标矩形高度（按逻辑单位,xSrc,ySrc:指定源位图左上角坐标（按逻辑单位）。,nSrcWidth,nSrcHeight:指定复制源位图宽度和高度（按逻辑单位）。,lpvBits：指向DIB数据图像指针。,lpbmi：指向BITMAPINFO结构指针。,fuColorUse:指定BITMAPINFO结构中bmiColors包含真实RGB值还是调色板中索引值。,dwrop:指定要执行光栅运算。,26,26/124,2.2,图像特效显示,图像特效显示就是利用人眼视觉特征，经过对图像分块，然后以不一样次序显示出来。实现图像特效显示基本思绪是将图像分成不一样小块，按一定方向或次序、分阶段地显示或擦除图像块。其中四个关键点是：,(1)怎样划分图像块；,(2)确定图像块操作次序；,(3)显示或去除图像块；,(4)在两个图像块操作之间延时。,延时目标是减慢图像显示速度，方便能够看出特效显示效果。,27,27/124,图像特效显示,图像扫描,扫描是最基本特效显示方式，它没有划分图像块，只是次序地,一行一行或一列一列地显示图像或去除图像。,效果图：,(a)自上而下扫描,(b)自下而上扫描 (c)自左向右扫描,(d)自右向左扫描,28,28/124,图像特效显示,图像移动,移动是将图像看作一个整体，显示时必须按物理次序进行，比如从上向下平移时，必须先显示下面图像，后显示上面图像。平移是以复制方法显示图像，每显示一次，复制行数就增加一行，直至显示完成。,效果图：,（a）水平向右移,(b)垂直向上移动,29,29/124,图像特效显示,交叉飞入,交叉飞入是将图像平分成上下两部分，显示时上部分水平右移，下部分水平左移。所以交叉分入基本原理和平移是相同，其不一样之处只是将图像进行了分块。,交叉飞入效果图,30,30/124,图像特效显示,中间扩张,中间扩张是在显示时候，先将图像分成两部分，将中间分界处显示在屏幕中间，并快速向上扫描上半部分图像，然后将图像完整显示在屏幕上，这么人们因为视觉生理特点就会看到中间扩张效果。,中间扩张效果图,31,31/124,图像特效显示,中间收缩,中间收缩效果是从屏幕上下两边同时向中间扫描，其原理类似于将图像平分成上下两部分，从屏幕上下边界，对图像同时进行向上扫描和向下扫描，直到在图像中间分界相遇。,中间收缩效果图,32,32/124,图像特效显示,栅条特效,栅条特效分为水平栅条和垂直栅条，其效果像是将两手交叉过程，栅条显示原理是先将图像分成若干行（或若干列），将奇数行（或奇数列）组成一组，将偶数行（或偶数列）组成一组，在显示时候，奇数行（或奇数列）从右向左平移。,（a）水平栅条显示,（b）垂直栅条显示,33,33/124,图像特效显示,图像渐显,图像渐显是先统计下列图像每个像素点灰度值，显示时候先将屏幕置黑，将循环显示图像n次，设n从0，1，2，256。每一次显示像素灰度值n/256倍，图像像素点计算一遍后，显示一次，重复执行上述过程，直至每一个屏幕上像素点灰度值全部和统计中值相同为止。,图像渐显效果图,34,34/124,图像特效显示,百叶窗特效显示,百叶窗显示开始有若干个显示起始线，每条起始线相隔相同距离，显示过程中每次各个起始线同时向下（或向右）显示一行（或一列）图像，直到图像显示完成。百叶窗显示原理是先将图像按行方向（或列方向）分成若干个相同部分（比如10块），显示时，先将每块头一行（或头一列）显示在屏幕对应位置上，然后以它们为起始，分别向下扫描（或向右扫描）直至图像全部显示完成。,（a）百叶窗自左向右,（b）百叶窗自上而下,35,35/124,图像特效显示,马赛克效果,马赛克显示是图像被分成许多小区域，小区域以杂乱无章次序显 示在屏幕上。其编程思想是，先将图像分成大小相同小区域，计算出每一块区域首地址，并统计下来。设置一个随机数，用来产生随机显示区域次序，每取得一个随机区域，就依据其首址显示这块区域图像，直至全部区域都最少显示一次。,马赛克显示效果图,36,36/124,第三章 图像几何变换,3.1 图像平移,3.2 图像镜像变换,3.3 图像缩放,3.4 图像转置,3.5 图像旋转,37,37/124,3.1 图像平移,图像平移（translation）是将图像中全部点都按照指定平移量水平、垂直移动。设初始坐标为（x0,y0）点经过平移（tx.ty）后，则点(x0,y0)坐标变为(x1,y1),如图所表示。,坐标原点,tx,X,ty,(x1,y1),Y,(x0,y0),像素平移示意图,显然（x0,y0）和（x1,y1）关系以下：,x1=x0+tx,y1=y0+ty,38,38/124,3.2 图像镜像变换,图像镜像变换分为两种：一个是水平镜像，另一个是垂直镜像。图像水平镜像操作是以原图像垂直中轴线为中心，将图像分为左右两部分进行对称变换；图像垂直镜像操作是以原图像水平中轴线为中心，将图像分为上下两部分进行对称变换。镜像变换后图高和宽都不变。,设图像高度为Height，宽度为Width，原图中(x0,y0)经过水平镜像后坐标将变成（Width-x0,y0），即：,x1=Width-x0,y1=y0,一样，（x0,y0）经过垂直镜像后坐标将变成为（x0,Height-y0），即：,x1=x0,y1=Height-y0,39,39/124,3.2 图像镜像变换,效果对比图,(a)原图,(b)水平镜像处理效果图 (c)垂直镜像处理效果图,镜像处理,40,40/124,3.3 图像缩放,假设图像X轴方向缩放比率是kx，Y轴方向缩放比率是ky，那么原图中点（x0,y0）对应与新图中点（x1,y1）为：,x1=x0*kx,y1=y0*ky,（1）kx1且ky1时，原图像被放大。因为放大图像时产生了新像素，可经过插值算法来近似处理。当kx=ky=2时，图像放大2倍，原图中某一个像素，对应新图4个像素。如图所表示：,（a）原图中某一个像素,（b）对应新图4个像素,41,41/124,3.3 图像缩放,（2）当kx1且ky1时，原图像被缩小。比如，当kx=ky=0.5时，图像被缩到二分之一大小，原图中4个像素对应新图中一个像素。此时缩小后图像中（0，0）像素对应于原图中（0，0）、（0，1）、（1，0）、（1，1）像素；以这类推。在原图基础上，每行隔一个像素取一点，每割一行进行操作。如图所表示：,（a）原图中某4个像素 （b）对应新图1个像素,图像缩小示意图,42,42/124,3.3 图像缩放,效果对比图,（a）原图 （b）长宽缩小0.5倍效果图,（c）长宽各放大2倍效果图,43,43/124,3.4 图像转置,图像转置（Transpose）操作是将图像像素x坐标和y坐标交换。该操作将改变图像高度和宽度，转置后图像高度和宽度将交换。,即 x1=y0,y1=x0,效果对比图：,图像转置,44,44/124,3.5 图像旋转,图像旋转必须指明图像绕着什么旋转。普通图像旋转是以图像中心为原点，旋转一定角度。旋转后，图像大小普通会改变。和图像平移一样，我们既能够把转出显示区域图像截去，也能够扩大图像范围以显示全部图像。,以下列图，点（x0,y0）经过旋转度后坐标变成（x1,y1）。,y,（x0,y0）,(x1,y1),x,O,图像旋转示意图,在旋转前：,x0=cos(),y0=sin(),旋转后：,x1=cos(-)=cos()cos()+sin()sin()=x0cos()+y0sin(),y1=sin(-)=sin()cos()-cos()sin()=-x0sin()+y0cos(),45,45/124,3.5 图像旋转,（a）原图,（b）旋转45效果图,旋转处理图,46,46/124,第四章 图像灰度变换,4.1 灰度变换基本方法,4.2 二值化和阈值处理,4.3 灰度线性变换,4.4 窗口灰度变换处理,4.5 分段线性变换,4.6 灰度直方图,47,47/124,4.1 灰度变换基本方法,灰度变换是指依据某种目标条件按一定变换关系逐点改变原图像中每一个像素灰度值方法。设原图像像素灰度值D=f（x，y）,处理后图像像素灰度值D=g(x,y)，则灰度增强可表示为：,g(x,y)=Tf（x，y)或 D=T(D),函数,T(D),称为灰度变换函数，它描述了输入灰度值和输出灰度值之间转换关系。一旦灰度变换函数确定，则确定了一个详细灰度增强方法。,灰度变换主要针对独立像素点进行处理，经过改变原始图像数据所占据灰度范围而使图像在视觉上得到良好改观，没有利用像素点之间相互空间关系。,依据，能够将灰度变换分为线性变换和非线性变换。,48,48/124,4.2 二值化和阈值处理,非零元素取一法,固定阀值法,双固定阀值法,49,49/124,非零元素取一法,非零元素取一法是最基础二值化算法。非0取1就是对于灰度图像f，若某像素灰度值为0，则其灰度值不变，仍为零；对于灰度值不为零像素，将其像素值全部变为255。非零元素取一法阈值T=1。,非零元素取一法变换函数表示式以下：,f(x)=,0 x1,255,255,g,0 1 255 f,非零元素取一法灰度变换函数,50,50/124,固定阀值法,1.理论基础,固定阀值法就是为灰度图像f设定一个阀值T，把灰度值小于给定阈值T像素置为0，大于阈值T像素置为255，从而对灰度图像实现二值化变换。,固定阀值法变换函数表示式以下：,f(x)=,0 xT,255,其中T为指定阈值。,g,255,f,255,0 T,固定阀值法灰度变换函数,51,51/124,固定阀值法,2.效果对比图,(a)原图 (b)以100作为阈值二值化图,52,52/124,双固定阀值法,双固定阀值法预先设置了两个阀值T1和T2，T1T2，当对图像进行处理时，假如某个像素灰度值小于T1则置0（或者255）；假如大于T1而且小于T2时，则置255（或者置0）；假如大于T2时，则置0（或者255）。,0-255-0型固定阀值法变换函数表示式以下：,0 x,T2,其中T1、T2为指定阈值。,g g,255 255,0 T1 T2 255 f,0 T1 T2 255 f,(a)0-255-0,型灰度变换函数,(b)255-0-255,型灰度变换函数,255 T1x1时，输出图像对比度将增大；,当a1时，输出图像对比度将减小；,当a=1且b不等于0时，使全部图像灰度值上移或下移，其效果是使整个图像更暗或更亮；,当a=1，b=0时，输出图像和输入图像相同；,当a=-1，b=255时，输出图像灰度恰好反转；,当a0时，暗区域将变亮，亮区域将变暗，点运算完成了图像求补运算。,56,56/124,4.3 灰度线性变换,效果图,反色变换处理结果图,57,57/124,4.4 窗口灰度变换处理,当图像中大部分像素灰度级在L,U范围内，少部分像素分布在小于L和大于U区间内时，可用两端“截取式”变换使小于灰度级L和大于等于灰度级U像素强行压缩为0和255。这种变换叫灰度窗口变换。灰度窗口变换也是一个常见点运算。它限定一个窗口范围，该窗口中灰度值保持不变；小于该窗口下限灰度值直接设置为0；大于该窗口上限灰度值直接设置为255。窗口灰度变换处理结合了双固定阀值法，与其不一样之处于于窗口内灰度值保持不变。,灰度窗口变换变换函数表示式以下：,0 x U,式中L表示窗口下限，U表示窗口上限。,58,58/124,4.4 窗口灰度变换处理,效果对比图,(a)改变窗口变换参数 （b）窗口变换效果图,59,59/124,4.5 分段线性变换,1理论基础,将图像灰度区间分成两段乃至多段，分别做线性变换称之为分段线性变换。,它灰度变换函数 如图所表示，函数表示式以下：,y1*x/x1 x x2,变换后灰度,255 g,y2,y1,0 x1 x2 255 f,60,60/124,4.5 分段线性变换,式中（x1，y1）和（x2，y2）是图中两个转折点坐标。,该变换函数运算结果是将原图在x1和x2之间灰度拉伸到y1和y2之间。经过有选择拉伸某段灰度区间，能够愈加灵活地控制图像灰度直方图分布，以改进输出图像质量。假如一幅图像灰度集中在较暗区域而造成图像偏暗，能够用灰度拉伸功效来拉伸（斜率1）物体灰度区间以改进图像质量；一样假如图像灰度集中在较亮区域而造成图像偏亮，也能够用灰度拉伸功效来压缩（斜率T，,g(i，j)=（5-1）,f(i，j)其它。,此式表明，若某点值与其邻域平均值相差超出 T，用平均值代替，进行平均处理，可去除噪声；不然还保留原值，不进行平均处理，从而降低含糊。这种算法对抑制椒盐噪声比较有效，同时也能交好地保护仅有微小灰度差图像细节。,75,75/124,超限邻域平均法,2效果图,超限邻域平均法（阈值为100）,76,76/124,NN均值滤波器,1理论基础,在本程序中当灰度图像f中以像素f(i，j)为中心NN屏蔽窗口(N=3，5，7)内平均灰度值为a时，无条件作f(i，j)=a处理，N由用户给定，且取N值越大，噪声降低越显著。但平均是以图像含糊为代价。,77,77/124,选择式掩模平滑,1理论基础,噪声消除法和邻域平均法在消除噪声同时，都不可防止带来平均化缺憾，致使尖锐改变边缘或线变得含糊。考虑到图像中目标物体和背景普通都含有不一样统计特征，即不一样均值和方差，为保留一定边缘信息，可采取自适应局部平滑滤波。这么能够得到很好图像细节。自适应平滑法是以尽可能不含糊边缘轮廓为目标。,选择式掩模平滑法取5 5窗口，如图5-2所表示。在窗口内以中心像素f(i,j)为基准点，制作4个五边形、4个六边形、1个边长为3正方形共9种形状屏蔽窗口，分别计算每个窗口内平均值及方差，因为含有尖锐边缘区域，方差必定较平缓区域为大,所以采取方差最小屏蔽窗口进行平均化，这种方法在完成滤波操作同时，又不破坏区域边界细节。这种采取9种形状屏蔽窗口，分别计算各窗口内灰度值方差，并采取方差最小屏蔽窗口进行平均化方法，也称作自适应平滑方法。,78,78/124,选择式掩模平滑,(a)周围9邻近 (b)左7邻近 (c)上7邻近,(d)右7邻近 (e)下7邻近 (f)左上7邻近,(g)右上7邻近 (h)右下7邻近 (i)左下7邻近,图5-2 55窗口选择式掩模平滑法,79,79/124,选择式掩模平滑,计算个掩模均值（ai）及方差（ki）。,ai=,ki=,式中m=1,2,3.,Q；Q为各掩模对应像素个数；,在此基础上，对ki排序，最小方差Kimin 随对应掩模灰度级均值ai 作为f(i，j)平滑输出(凡含有尖锐边缘区域，方差必定较平缓区域为大)G(i，j)。,80,80/124,选择式掩模平滑,2.处理效果图,(a)原图,（b）7*7均值滤波 (c)选择式掩模平滑法,81,81/124,5.4 中值滤波,N*N中值滤波器,十字型中值滤波器,NN最大值滤波器,82,82/124,N*N中值滤波器,1.理论基础,中值滤波普通采取一个含有奇数个点滑动窗口，将窗口中各点灰度值中值来替换指定点（普通是窗口中心点）灰度值。假设窗口内有五点，其值依次为1，4，6，0，7，重新排序后（从小到大）为0，1，4，6，7，则Med1，4，6，0，7=4。,计算灰度图像f中以像素f(i，j)为中心NN屏蔽窗口(N=3，5，7)内灰度中值为u，无条件作f(i，j)=u处理，N由用户给定。,83,83/124,十字型中值滤波器,1.理论基础,计算灰度图像f中以像素f(i，j)为中心得十字型屏蔽窗口(十字型纵向和横向长度为N，N=3，5，7)内灰度值中值u，无条件作f(i，j)=u处理，N由用户给定。,55十字型中值滤波器,84,84/124,NN最大值滤波器,1.理论基础,本程序计算灰度图像f中以像素f(i，j)为中心NN屏蔽窗口 (N=3，5，7)内灰度最大值为u，无条件作f(i，j)=u处理，N由用户给定,2.处理效果,NN中值滤波、十字型中值滤波、N*N最大值滤波之后效果大致相同。,中值滤波图,85,85/124,5.5 产生噪声,随机噪声,椒盐噪声,86,86/124,随机噪声,1.理论基础,本程序经过计算机所产生随机数给图像加噪声。,循环取得各点像素值；取得随机数（rand()/1024）；该点像素值加上随机数；,87,87/124,椒盐噪声,1.理论基础,本程序经过计算机所产生随机数大小来给图像加噪声。,循环取得各点像素值；若产生随机数大于特定值，把该点置0；,2,.处理效果,(a)随机噪声 (b)椒盐噪声,88,88/124,第六章 图像锐化处理及边缘检测,6.1 图像边缘锐化基本方法,6.2 微分运算,6.3 梯度锐化,6.4 边缘检测,6.5 高斯-拉普拉斯算子,89,89/124,6.1 图像边缘锐化基本方法,利用计算机进行图像锐化处理有两个目标，一是与柔化处理相反，增强图像边缘，使含糊图像变得愈加清楚起来，颜色变得鲜明突出，图像质量有所改进，产生更适合人观察和识别图像，二是希望经过锐化处理后，目标物体边缘鲜明，方便于计算机提取目标物体边界、对图像进行分割、目标区域识别、区域形状提取等，为图像了解和分析打下基础。,90,90/124,6.2 微分运算,纵向微分运算,横向微分运算,双方向一次微分运算,91,91/124,6.2 微分运算,边缘是由相邻域灰度级不一样像素点组成，若想增强边缘，就应该突出相邻点间灰度级改变。以下列图所表示，不难发觉原图中左边暗，右边亮，中间存在着一条显著边界。,0 0 1 255 255 255 255,1 1 1 254 253 254 254,0 0 0 255 255 253 253,1 1 0 254 254 254 254,为了与f(x,y)表示法相一致，图像f(i,j)中i代表列，j代表行。假如用右列减去左列，即每一个像素值为：,G(i,j)=f(i,j)-f(i-1,j),92,92/124,6.2 微分运算,结果以下：,0 1 254 0 0 0,0 0 253 -1 1 0,0 0 255 0 -2 0,0-1 254 0 0 0,G(i,j)=f(i,j)-f(i-1,j)就起到了垂直边缘检测作用。,如对于上面那幅图像转置，得到以下列图像数据：,0 1 0 1,0 1 0 1,1 1 0 0,255 254 255 254,255 253 255 254,255 254 253 254,255 254 253 254,G(i,j)=-f(i,j-1)+f(i,j)，图像上得到一条很显著亮边。,93,93/124,纵向微分运算,1理论基础,对灰度图像在纵方向进行微分实现了将图像向左平移一个像素，再用原图像减去平移后图像。相减结果反应了原图像亮度改变率大小。原图像中像素值保持不变区域，相减结果为零，即像素为黑；原图像中像素灰度值改变猛烈区域，相减后得到较大改变率，对应像素很亮，而且像素灰度值差异越大，则得到像素就越亮，所以图像垂直边缘得到增强。,该算法数学表示式为：,G(i,j)=f(i,j)-f(i,j-1)+,该算法用以下卷积核：,94,94/124,纵向微分运算,2效果对比图,(a)原图 (b)纵向微分运算,95,95/124,横向微分运算,1理论基础,该算法数学表示式为：,G(i,j)=f(i,j)-f(i,j-1)+,该算法用以下卷积核：,96,96/124,横向微分运算,2效果对比图,横向微分运算,97,97/124,双方向一次微分运算,1理论基础,对灰度图像f在纵方向和横方向两个方向进行微分。该算法是同时增强水平和垂直方向边缘。该算法数学表示式为：,G(i,j)=sqrtf(i,j)-f(i,j-1)*f(i,j)-f(i,j-1)+f(i,j)-f(i-1,j)*f(i,j)-f(i-1,j),对于含小数G(i,j)可四舍五入。,该算法用以下卷积核：,使用水平方向卷积核得出像素值为m，使用垂直方向卷积核得出像素值为n，该像素边界强度为：sqrt（mm）+（nn）。该算法是同时增强水平和垂直方向边缘。,水平（i方向）垂直（j方向）,98,98/124,双方向一次微分运算,2处理结果,(a)原图 (b)纵向微分运算,(c)横向微分运算 (d)双方向一次微分运算,99,99/124,6.3 梯度锐化,直接以梯度值代替,辅以门限判断,给边缘要求一个特定灰度级,给背景要求灰度级,依据梯度二值化图像,100,100/124,直接以梯度值代替,利用双方向一次微分运算，算出梯度后让梯度值等于该点灰度值。,即f(i,j)=Gf(i,j)。,Gf(i,j)=f(i,j)-f(i-1,j)2+f(i,j)-f(i,j-1)2 1/2,或Gf(i,j)=|f(i,j)-f(i-1,j)|+|f(i,j)-f(i,j-1)|,这种方法直截了当。但在均匀区域，因梯度值Gf(i,j)很小，会表现出图像很暗特征，这在一些场所是不宜。就像在前面看到微分效果一样，除了黑色背景，几乎看不出什么边界，所以在这里我们就不作介绍了。,101,101/124,辅以门限判断,1.理论基础,门限判断梯度锐化公式以下：,Gf(i,j)=f(i,j)-f(i-1,j)2+f(i,j)-f(i,j-1),2,1/2,Gf(i,j)+100;Gf(i,j)T,g(i,j)=,f(i,j);其它,102,102/124,辅以门限判断,该方法基本上不破坏图像背景，又可增强边缘。这是因为Gf(i,j)表示是两个像素点之间灰度差大小，也就是梯度大小。对于图像而言，物体和物体之间，背景和背景之间梯度改变普通很小，灰度改变较大地方普通集中在图像边缘上，也就是物体和背景交接地方，当设定一个适当阈值T，Gf(i,j)大于T就认为该像素点处于图像边缘，对结果加100，以使边缘变亮，而对于Gf(i,j)小于T就认为像素点是同类像素（同是物体或同是背景）。这么既增亮了物体边界，同时又保留了图像背景原来状态。,103,103/124,辅以门限判断,2处理效果图,(a)原图 (b)辅以门限判断效果图,104,104/124,给边缘要求一个特定灰度级,1.理论基础,Gf(i,j)=f(i,j)-f(i-1,j)2+f(i,j)-f(i,j-1),2,1/2,La;Gf(i,j)T,g(i,j)=f(i,j);其它,La为一指定灰度值。这种处理实际上是门限锐化一个特殊形式，它将边界灰度值统一化，这么能够使边界愈加清楚显著。该方法基本上不破坏图像背景，又可找到边缘，并依据需要增强边缘。,105,105/124,给边缘要求一个特定灰度级,2效果对比图,（a）原图 （b）给边缘要求一个特定灰度级,106,106/124,给背景要求灰度级,Gf(i,j);Gf(i,j)T,g(i,j)=Lb;其它,Lb为一对背景指定灰度值。这种处理原理和上一节介绍给边缘要求一个特定灰度级原理是一样，只不过这种处理将图像背景灰度值统一化了，所以在这里不重复介绍了。,107,107/124,依据梯度二值化图像,1.理论基础,对阶跃边缘，在边缘点其一阶导数取极值。由此，我们对数字图像f(i,j)每个像素取它梯度值G(i,j)，适当取门限T作以下判断：若G(i,j)T则(i,j)点为阶跃状边缘点，G(i,j)称为梯度算子边缘图像,La;Gf(i,j)T,g(i,j)=Lb;其它,La 和 Lb 意义同上。梯度是向量，各向同性。梯度方向对应于F(i,j)最大改变率方向上，即,Q=arctan/,108,108/124,依据梯度二值化图像,梯度幅度百分比于邻像素灰度级，在灰度陡变区域，梯度值大，在灰度相同区，梯度值小，在灰度级为常数区，梯度为零。这么处理能够使图像锐化结果愈加清楚，把图像中我们关心部分突出出来，去除了我们不感兴趣部分。,109,109/124,依据梯度二值化图像,2处理效果图,(a)原图 （b）依据梯度二值化效果图,110,110/124,6.4 边缘检测,Roberts边缘检测算子,Sobel边缘检测算子,Prewitt边缘检测算子,Krisch边缘检测,111,111/124,Roberts边缘检测算子,1.理论基础,Roberts边缘算子采取是对角方向相邻两个像素之差。从图像处理实际效果来看，边缘定位准，对噪声敏感。Roberts边缘检测算子是一个利用局部差分算子寻找边缘算子，它由下式给出：,Roberts算子：Gi,j=|fi,j-fi+1,j+1|+|fi+1,j-fi,j+1|;,Gi,j=(fi,j-fi+1,j+1)2+(fi+1,j-fi,j+1),2,1/2,;,其中Gi,j表示处理后(i,j)点灰度值，fi,j表示处理前该点灰度