MATLAB在数字图像处理中的应用.ppt

B2=blkproc(B,[8 8],P1.*x,mask); I2=blkproc(B2,[8 8],P1*x*P2,T,T); imshow(I),figure,imshow(I2) M=I2-I; figure,imshow(mat2gray(M)) 原始图像 经压缩、解压后的图像 相减所得图像 2.图像相加运算 设有图像f(x,y)和h(x,y)，它们的和为 g(x,y)=f(x,y)+h(x,y) 图像相加的一个重要应用是对同一场景的多幅图像求平均值。这种方法经常用来有效地降低加性随机躁声的影响。图像相加也可以用来将一幅图像的内容叠加到另一幅图像上去，以达到二次暴光的效果。 例1 实现图像相加。 I=imread(saturn.tif); figure,imshow(I) [m,n]=size(I); J(m,n)=0; for i=1:3 X=imnoise(I,gaussian); figure,imshow(X) Y=double(X); J=J+Y/3; end figure,imshow(mat2gray(J)) 图6.4.1 索引图像 图6.4.2 灰度图像 图6.4.3 真彩图像 图6.4.4 二值图像 三. 边缘检测 常用的边缘检测算子有Robert算子、Sobel算子、Prewitt算子、LOG算子、Canny算子，以下分别进行介绍。 1. Robert算子 Robert算子是一种利用局部差分算子寻找边缘的算子，其模板如图所示。 ? 1 0 0 -1 0 1 -1 0 Robert算子 在MATLAB中可以由函数实现，语法格式如下： BW=edge(I,’roberts’) BW=edge(I,’roberts’,thresh) BW=edge(I,’roberts’)自动选择阈值用Robert算子进行边缘检测。 BW=edge(I,’roberts’,thresh)根据所指定的敏感度阈值thresh用Robert算子进行边缘检测，它忽略了所有小于阈值的边缘。当thresh为空([])时，自动选择阈值。 2.Sobel算子 Sobel算子的两个卷积计算核如图所示，图像中的每个点都用这两个核做卷积，第一个核对水平边缘响应最大，第二个核对垂直边缘响应最大。 ? 1 2 1 0 0 0 -1 -2 -1 1 0 -1 2 0 -2 1 0 -1 edge函数实现的语法格式如下： BW=edge(I,’sobel’) BW=edge(I,’sobel’,thresh) BW=edge(I,’sobel’,thresh,direction) BW=edge(I,’sobel’)自动选择阈值用Sobel算子进行边缘检测。 BW=edge(I,’sobel’,thresh)根据所指定的敏感度阈值thresh用Soble算子进行边缘检测，它忽略了所有小于阈值的边缘。当thresh为空([])时，自动选择阈值。 BW=edge(I,’sobel’,thresh,direction)根据所指定的敏感度阈值thresh，在所指定的方向direction上，用Sobel算子进行边缘检测。direction可取的字符串值为’horizontal’（水平方向）、’vertial’（垂直方向）或’both’（两个方向）。 3. Prewitt算子 Prewitt算子的两个卷积计算核如图所示。 1 0 -1 1 0 -1 1 0 -1 1 1 1 0 0 0 -1 -1 -1 edge函数实现的语法格式如下： BW=edge(I,’prewitt’) BW=edge(I,’prewitt’,thresh) BW=edge(I,’prewitt’,thresh,direction) BW=edge(I,’prewitt’)自动选择阈值用Prewitt算子进行边缘检测。 BW=edge(I,’prewitt’,thresh)根据所指定的敏感度阈值thresh用Prewitt算子进行边缘检测，它忽