数字图像处理-第二章课件.ppt

1、chart 02Fundamentals2.1 Digital Image Representation 数学图像的表示数学图像的表示2.2 Reading Images 读图像读图像2.3 Displaying Images 显示图像显示图像2.5 Data Classes (omitted,self-study)数据的类型（略，自学）数据的类型（略，自学）2.6 Image Types 图像的类型图像的类型2.7 Converting between Data Classes and Image Types 图像类型与数据类型之间的转换图像类型与数据类型之间的转换2.8 Array Ind

2、exing (omitted,self-study)数组的索引（下标）（略，自学）数组的索引（下标）（略，自学）2.9 Some Important Standard Arrays (omitted)一些重要的标准数组（略）一些重要的标准数组（略）2.10 M-Function Programming (omitted,self-study)M-函数编程（略，自学）函数编程（略，自学）An image may be defined as a two-dimensional function,f(x,y),where x and y are spatial(plane)coordinates,a

3、nd the amplitude of f at any pair of coordinates(x,y)is called the intensity or gray level of the image at that point.一幅图像被定义为一个二元（维）函数一幅图像被定义为一个二元（维）函数f(x,y),x和和y是空间（平面）是空间（平面）坐坐标，在坐标对坐坐标，在坐标对(x,y)处的幅值处的幅值f称为图像在该点处的亮度值或灰度等称为图像在该点处的亮度值或灰度等级（或灰度级）级（或灰度级）Color images are formed by a combination of ind

4、ividual 2-D images.For example,in the RGB color system,a color image consists of three(red,green,and blue)individual component images.颜色图像可由单个颜色图像可由单个2维图像组合而成。例如在维图像组合而成。例如在RGB颜色系统中，颜色系统中，颜色图像是由三个单色成份（即红、绿和蓝）的图像组成。颜色图像是由三个单色成份（即红、绿和蓝）的图像组成。For this reason,many of the techniques developed for monochr

5、ome images can be extended to color images by processing the three component images individually.据此，许多针对单色图像所开发的技术，要以扩展到彩色图像上。据此，许多针对单色图像所开发的技术，要以扩展到彩色图像上。即通过单独对三个成份图像分别处理。即通过单独对三个成份图像分别处理。2.1 Digital Image Representation实例说明I=imread(2353.jpg);imshow(I);I1=I;I2=I;I3=I;I1(:,:,1)=I(:,:,1);I1(:,:,2)=0;

6、I1(:,:,3)=0;I2(:,:,2)=I(:,:,2);I2(:,:,1)=0;I2(:,:,3)=0;I3(:,:,3)=I(:,:,3);I3(:,:,1)=0;I3(:,:,2)=0;figure;imshow(I);figure;imshow(I1);figure;imshow(I2);figure;imshow(I3);实现代码：实现代码：An image may be continuous with respect to the x-and y-coordinates,and also in amplitude.一幅图像关于一幅图像关于x和和y坐标可以是连续的（即位置上连接）

7、，幅值也坐标可以是连续的（即位置上连接），幅值也可以是连续的（即亮度上是连续的）。可以是连续的（即亮度上是连续的）。Converting such an image to digital form requires that the coordinates,as well as the amplitude,be digitized.将上述连续图像转换成数字形式，需要将坐标将上述连续图像转换成数字形式，需要将坐标x、y和幅度和幅度f数字化。数字化。Digitizing the coordinate values is called sampling;digitizing the amplitud

8、e values is called quantization.将坐标值数字化称为采样；幅值的数字化称为量化。将坐标值数字化称为采样；幅值的数字化称为量化。Thus,when x,y,and the amplitude values of f are all finite,discrete quantities,we call the image a digital image.因而，当因而，当x、y和幅度和幅度f都是有限的离散的量时，我们称图像为数字都是有限的离散的量时，我们称图像为数字图像。图像。2.1 Digital Image Representation The result of

9、sampling and quantization is a matrix of real numbers.We use two principal ways in this book to represent digital images.采样和量化的结果将得到一个实数矩阵。在本书中我们使用两个采样和量化的结果将得到一个实数矩阵。在本书中我们使用两个主要方法来表示图像。主要方法来表示图像。Assume that an image f(x,y)is sampled so that the resulting image has M rows and N columns.We say that

10、the image is of size MN.The values of the coordinates(x,y)are discrete quantities.设图像设图像f(x,y)是一采样后的图像，因此它有是一采样后的图像，因此它有M行和行和N列。我们称图像列。我们称图像的大小（或尺寸）是的大小（或尺寸）是MN，坐标，坐标(x,y)的值是离散的量。的值是离散的量。For notational clarity and convenience,we use integer values for these discete coordinates.为表示上的清晰和方便，我们使用整数值来表示这

11、些离散的坐标为表示上的清晰和方便，我们使用整数值来表示这些离散的坐标 2.1.1 Coordinate Conventions图 采样网格(a)正方形网格;(b)正六角形网格(a)(b)图2-2 采样示意图 采样行采样列像素行间隔采样间隔实例说明434348图2-3 量化示意图（a）量化；(b)量化为8 bit 连续灰度值 量化值 (整数值)灰度标度 灰度量化Zi1ZiZi1qi1qi125525412812710(a)(b)图图2-4 不同采样点数对图像质量的影响不同采样点数对图像质量的影响（a）原始图像）原始图像(256256)；（；（b）采样图像）采样图像1(128128)；（c）采样图

12、像采样图像2(6464)；（d）采样图像）采样图像3(3232)；（e）采样图像）采样图像4(1616)；（；（f）采样图像采样图像5(88)图图2-5 不同量化级别对图像质量的影响不同量化级别对图像质量的影响（a）原始图像原始图像(256色色)；（b）量化图像量化图像1(64色色)；（c）量化图像量化图像2(32色色)；（d）量化图像量化图像3(16色色)；（e）量化图像量化图像4(4色色)；（f）量化图像量化图像5(2色色)一般，当限定数字图像的大小时一般，当限定数字图像的大小时,为了得到质量较好的图为了得到质量较好的图像可采用如下原则：像可采用如下原则：（1）对缓变的图像，对缓变的图像，

13、应该细量化，应该细量化，粗采样，粗采样，以避免假以避免假轮廓。轮廓。（2）对细节丰富的图像，对细节丰富的图像，应细采样，应细采样，粗量化，粗量化，以避免以避免模糊（混叠）。模糊（混叠）。对于彩色图像，是按照颜色成分对于彩色图像，是按照颜色成分红（红（R）、绿）、绿（G）、蓝（）、蓝（B）分别采样和量化的。若各种颜色成分均按分别采样和量化的。若各种颜色成分均按8 bit量化，即每种颜色量级别是量化，即每种颜色量级别是 256，则可以处理则可以处理256256256=16 777 216种颜色。种颜色。In many image processing books,the image origin

14、is defined to be at(x,y)=(0,0).The next coordinate values along the first row of the image are(x,y)=(0,1).许多图像处理教材中，图像原点被定义为许多图像处理教材中，图像原点被定义为(x,y)=(0,0).即参考起始点即参考起始点的坐标值定义为的坐标值定义为(0,0),位置从第位置从第0行第行第0列位置开始计数，沿着图像第一列位置开始计数，沿着图像第一行的下一个坐标值是行的下一个坐标值是(x,y)=(0,1)，即第，即第0行第行第1列的位置。列的位置。It is important to ke

15、ep in mind that the notation(0,1)is used to signify the second sample along the first row.It does not mean that these are the actual values of physical coordinates when the image was sampled.需要记住的一个要点是：符号需要记住的一个要点是：符号(0,1)被用于标记沿第一行的第二个被用于标记沿第一行的第二个采样点。当图像被采样时，采样点。当图像被采样时，(0,1)等坐标值并不代表物理坐标的真实值，等坐标值并不

16、代表物理坐标的真实值，它代表的仅仅是采样点的相对位置。如它代表的仅仅是采样点的相对位置。如(i,j)坐标表示第坐标表示第i行第行第j个的采样个的采样点。点。2.1.1 Coordinate Conventions实例说明 Figure 2.1(a)shows this coordinate convention.Note that x ranges from 0 to M-1,and y from 0 to N-1,in integer increments.下图显示坐标约定。注意下图显示坐标约定。注意x的范围是的范围是0到到M-1，y的范围是的范围是N-1，按整，按整数值递增数值递增Figu

17、re 2.1(a)The coordinate convention used in the toolbox to denote arrays is different from the preceding paragraph in two minor ways.工具箱中的坐标约定和前面段落中的坐标约定，在两个小的方面工具箱中的坐标约定和前面段落中的坐标约定，在两个小的方面上有所不同。上有所不同。First,instead of using(x,y),the toolbox uses the notation(r,c)to indicate rows and columns.首先工具箱使用符号

18、首先工具箱使用符号(r,c)来表示行数和列数，而不是使用来表示行数和列数，而不是使用(x,y)The other difference is that the origin of the coordinate system is at(r,c)=(1,1);另一个不同之处就是在工具箱的坐标系统中，坐标的原点是在另一个不同之处就是在工具箱的坐标系统中，坐标的原点是在(r,c)=(1,1)，即第一行的第一个列位置上的采样点是从，即第一行的第一个列位置上的采样点是从(1,1)开始计数，也开始计数，也就是说，第一行记为第就是说，第一行记为第1行，而不是第行，而不是第0行（别的教材将图像的第一行行（别的

19、教材将图像的第一行记为第记为第0行），同样列亦如此。行），同样列亦如此。thus,r ranges from 1 to M,and c from 1 to N,in integer increments.This coordinate convention is shown in Fig.2.1(b)因此，因此，r的范围从的范围从1到到M，c的范围是从的范围是从1到到N，按整数值递增。这一，按整数值递增。这一坐标约定如图坐标约定如图2.1(b)所示所示Figure 2.1(b)Representation for a digitized image function:数字图像函数的表示数字图像

20、函数的表示 2.1.2 Image as Matrices A digital image can be represented naturally as a MATLAB matrix:在在MATLAB的的IPT中，数字图像表示成一个如下的中，数字图像表示成一个如下的MATLAB矩阵矩阵 Images are read into the MATLAB environment using function imread,whose syntax is:使用函数使用函数imread 将图像读入到将图像读入到MATLAB工作环境中，其语法如下：工作环境中，其语法如下：A=imread(,fmt)X

21、,map=imread(.).=imread().=imread(.,idx)CUR or ICO A,map,alpha=imread(.)CUR or ICO .=imread(.,idx)GIF .=imread(.,frames,idx)GIF .=imread(.,BackgroundColor,BG)PNG A,map,alpha=imread(.)PNG .=imread(.,idx)TIFF .=imread(.,PixelRegion,ROWS,COLS)TIFF2.2 Reading Images A=imread(,fmt)reads a grayscale or col

22、or image from the by the string.If the not in the current directory,or in a directory on the MATLAB path,specify the full pathname.A=imread(,fmt)从字符串参数所指定的文件中读取一幅灰度图或从字符串参数所指定的文件中读取一幅灰度图或彩色图像。如果文件不在当前路径，也不在彩色图像。如果文件不在当前路径，也不在MATLAB的搜索路径中，的搜索路径中，则应指定文件的完整路径。则应指定文件的完整路径。The return value A is an array

23、containing the image data.If the a grayscale image,A is an M-by-N array.If the a truecolor image,A is an M-by-N-by-3 array.返回值返回值A是是一个包含图像数据的矩阵（数组）。如果图像文件是是一个包含图像数据的矩阵（数组）。如果图像文件是一个灰度图，则是一个灰度图，则A为一为一MN矩阵（数组），如果图像文件是一个矩阵（数组），如果图像文件是一个真彩色图像，则真彩色图像，则A是个是个MN3的矩阵（注的矩阵（注MATLAB中数组和矩阵中数组和矩阵在内存分配上没有区别，有时称数组，

24、有时称矩阵）在内存分配上没有区别，有时称数组，有时称矩阵）2.2 Reading Images Supported formats MATLAB所支持的图像文件格式所支持的图像文件格式2.2 Reading Images For example,the command line:2.2 Reading Imagesf=imread(Fig0203(a)(chest-xray).tif);Read the image chest-xray into image array f.When,as in the preceding command line,no path information is

25、 included in,imread reads the the current directory and,if that fails,it tries to find the the MATLAB search path.在前述的命令行中，当文件名中没有路径信息时，在前述的命令行中，当文件名中没有路径信息时，imread就从当前路径就从当前路径中读取，如果读取失败，就试着从中读取，如果读取失败，就试着从MATLAB的搜索路径中读取。的搜索路径中读取。The simplest way to read an image from a specified directory is to inc

26、lude a full or relative path to that directory in.For example:从指定路径中读取图像的最简单方法就是在文件名中指定完整路径或相从指定路径中读取图像的最简单方法就是在文件名中指定完整路径或相对路径。例如：对路径。例如：f=imread(E:教学教学数字图像处理数字图像处理chart02(chest-xray).tif);Function size gives the row and column dimensions of an image:函数给出图像的行数和列数，即图像的行列维数。函数给出图像的行数和列数，即图像的行列维数。2.2

27、Reading ImagesM N=size(f)M=494N=600s=size(f)s=494 600 Images are displayed on the MATLAB desktop using function imshow,which has the basic syntax:在在MATLAB桌面环境下使用函数桌面环境下使用函数imshow显示图像，其语法如下：显示图像，其语法如下：imshow(f,G)where f is an image array,and G is the number of intensity levels used to display it.If G

28、 is omitted,it defaults to 256 levels.F是图像矩阵，是图像矩阵，G是显示图像所用的灰度等级数目。如果是显示图像所用的灰度等级数目。如果G省略，缺省为省略，缺省为256 Using the syntax:使用语法：使用语法：imshow(f,low high)displays as black all values less than or equal to low,and as white all values greater than or equal to high.The values in between are displayed as inte

29、rmediate intensity values using the default number of levels.将小于或等于将小于或等于low的灰度级显示为黑，将大于或等于的灰度级显示为黑，将大于或等于high的灰度级显示为白。的灰度级显示为白。low和和high之间的值按缺省灰度级数目显示为黑、白之间的中间灰度级。之间的值按缺省灰度级数目显示为黑、白之间的中间灰度级。2.3 Displaying Images Finally,the syntax:imshow(f,)sets variable low to the minimum value of array f and high

30、 to its maximum value.最后语法最后语法 imshow(f,)将将low和和high分别设置为矩阵分别设置为矩阵f的最小值和最大值。的最小值和最大值。This form of imshow is useful for displaying images that have a low dynamic range or that have positive and negative values.当图像的灰度动态范围较低或图像有正、负值时，这一语法形式非常有用当图像的灰度动态范围较低或图像有正、负值时，这一语法形式非常有用 For example:2.3 Displaying

31、 Imagesfchest=imread(Fig0203(a)(chest-xray).tif);subplot(1,2,1),imshow(fchest),subplot(1,2,2),imshow(fchest,)For example:2.3 Displaying Imagesfrose=imread(Fig0206(a)(rose-original).tif);imshow(frose)Images are written to disk using function imwrite,which has the following basic syntax:将图像写入硬盘的函数是将图像

32、写入硬盘的函数是imwrite，其基本语法：，其基本语法：imwrite(f,)For example,the following command writes f to a TIFF file anmed patient10_run1:例如，如下命令将例如，如下命令将f写入名为写入名为patient10_run1的的TIFF格式的文格式的文件中。件中。imwrite(f,patient10_run1,tif)or alternatively,imwrite(f,patient10_run1.tif)2.4 Writing Images A more general imwrite synta

33、x applicable only to JPEG images is:如果只针对如果只针对JPEG图像，能适合的一般语法是：图像，能适合的一般语法是：imwrite(f,quality,q)wher q is an integer between 0 and 100(the lower the number the higher the degradation due to JPEG compression)q是一个介于是一个介于0到到100的整数，其值越小，的整数，其值越小，JPEG图像压缩的降质越厉害。图像压缩的降质越厉害。In order to reduce storage and t

34、ransmission time,it is important that the images be compressed as much as possible while not degrading their visual appearance beyond a reasonable level.为了减少图像文件的存贮容量和传输时间，重要的一点就是要将图像尽为了减少图像文件的存贮容量和传输时间，重要的一点就是要将图像尽可能的压缩到合适的程度，以便不影响视觉效果。可能的压缩到合适的程度，以便不影响视觉效果。In this case “reasonable”means no percept

35、ible false contouring.For example,for q=25 the applicable syntax is:本例中的合适程度，意思是指没有可感知到的伪轮廓。如本例中的合适程度，意思是指没有可感知到的伪轮廓。如q=25时的所用时的所用的语法是：的语法是：imwrite(f,bubble25.jpg,quality,25)2.4 Writing Imagesfbubble=imread(Fig0204(a)(bubbles-q-100jpg).tif);imwrite(fbubble,bubbles0.jpg,quality,0);imwrite(fbubble,bub

36、bles5.jpg,quality,5);imwrite(fbubble,bubbles15.jpg,quality,15);imwrite(fbubble,bubbles25.jpg,quality,25);imwrite(fbubble,bubbles50.jpg,quality,50);2.4 Writing Imagesimshow(fbubble)imshow(f50)imshow(f25)imshow(f15)imshow(f5)imshow(f0)f0=imread(bubbles0.jpg);f5=imread(bubbles5.jpg);f15=imread(bubbles1

37、5.jpg);f25=imread(bubbles25.jpg);f50=imread(bubbles50.jpg);2.4 Writing ImagesOriginal imageq=50q=25q=0q=5q=152.4 Writing ImagesOriginal imageq=10q=8q=0q=4q=6The toolbox supports four types of images:工具箱支持四种图像类型工具箱支持四种图像类型 Intensity images 灰度图像（或叫亮度图像）灰度图像（或叫亮度图像）Binary images 二值图像（或叫二进值图像，或叫黑白图像）二值图

38、像（或叫二进值图像，或叫黑白图像）Indexed images 索引图像索引图像 RGB images RGB图像图像Most monochrome image processing operations are carried out using binary or intensity images,so our initial focus is on these two image types.大多数单色图像处理操作都是针对二值图像或灰度图像进行的。大多数单色图像处理操作都是针对二值图像或灰度图像进行的。因此我们首先集中在这两类图像的研究上。因此我们首先集中在这两类图像的研究上。2.6 I

39、mage Types实例说明灰度图像灰度图像实例说明灰度图像灰度图像实例说明索引图像索引图像实例说明RGB图像图像An intensity image is a data matrix whose values have been scaled to represent intensities.在在MATLAB中，灰度图像是一个数值矩阵。它的值归一化中，灰度图像是一个数值矩阵。它的值归一化到其所表示的图像。到其所表示的图像。When the elements of an intensity image are of class uint8,or class uint16,they have i

40、nteger values in the range 0,255 and 0,65535,respectively.当图像类型（即存贮各像素值的数据类型）是当图像类型（即存贮各像素值的数据类型）是uint8或或uint16时，其整数值的取值范围分别是时，其整数值的取值范围分别是0,255和和0,65535。Values of scaled,class double intensity images are in the range 0,1 by convention.通过约定，归一化的双精度类型的灰度图像的数据取值范围通过约定，归一化的双精度类型的灰度图像的数据取值范围是是 0,12.6.1

41、Image TypesA binary image is a logical array of 0s and 1s.Thus,an array of 0s and 1s whose values are of data class,say,uint8,is not considered a binary image in MATLAB.二值图像是由二值图像是由0和和1组成一个逻辑矩阵。其元素值是逻辑值。因此在组成一个逻辑矩阵。其元素值是逻辑值。因此在MATLAB中，如果由数据类（而不是图像类）的中，如果由数据类（而不是图像类）的0和和1组成的矩阵，不能将组成的矩阵，不能将其视为一个二值图像。注

42、意，当其视为一个二值图像。注意，当0和和1是数据类的逻辑值时，则可以将矩阵是数据类的逻辑值时，则可以将矩阵视为二值图像。视为二值图像。A numeric array is converted to binary using function logical.Thus,if A is a numeric array consisting of 0s and 1s,we create a logical array B using the statement:数值矩阵可通过函数数值矩阵可通过函数 logical转换成二值逻辑矩阵，因此，如果一个矩阵转换成二值逻辑矩阵，因此，如果一个矩阵A由数值由数

43、值(非逻辑非逻辑)0和和1组成，我们可以使用如下声明来创建逻辑矩阵：组成，我们可以使用如下声明来创建逻辑矩阵：B=logical(A)To test if an array is logical we use the islogical function:测试一个矩阵是否是一个逻辑型的，可以使用测试一个矩阵是否是一个逻辑型的，可以使用islogical函数：函数：islogical(C)If C is a logical array,this function returns a 1.Otherwise it returns a 0.如果如果C是一个逻辑矩阵，则该函数返回是一个逻辑矩阵，则该函

44、数返回1，否则返回，否则返回0.2.6.2 Binary ImagesImage Types and Type Conversionsdither：Convert image,increasing apparent color resolution by dithering 通过抖动技术，多外观上增加灰度的等级。通过抖动技术，多外观上增加灰度的等级。double：Convert data to double precision 将数据转换成双精度型。将数据转换成双精度型。gray2ind：Convert grayscale or binary image to indexed image 将灰

45、度图或二值图像转换成索引图。将灰度图或二值图像转换成索引图。graythresh：Global image threshold using Otsus method 使用使用Otsu方法将图像进行全局阈值化处理成二值图像。方法将图像进行全局阈值化处理成二值图像。im2bw：Convert image to binary image,based on threshold 使用阈值化方法，将图像转换成二值图像。使用阈值化方法，将图像转换成二值图像。im2double：Convert image to double precision 将图像转换成双精度型矩阵。将图像转换成双精度型矩阵。im2int

46、16：Convert image to 16-bit signed integers 将图像转换成将图像转换成int16型矩阵型矩阵im2java：Convert image to Java image 将将MATLAB图像转换成图像转换成Java类型图像类型图像im2java2d：Convert image to Java buffered image 将将MATLAB转换成转换成Java buffered图像图像2.7 Converting betwwen Data Classes and Image Typesim2single：Convert image to single preci

47、sion 将图像转换成单精度矩阵将图像转换成单精度矩阵im2uint16：Convert image to 16-bit unsigned integers 将图像转换成将图像转换成uint16型矩阵型矩阵im2uint8：Convert image to 8-bit unsigned integers 将图像转换成将图像转换成uint8型矩阵型矩阵2.7 Converting betwwen Data Classes and Image Typesind2gray：Convert indexed image to grayscale image 将图像转换成灰度图像将图像转换成灰度图像ind

48、2rgb：Convert indexed image to RGB image 将图像转换成将图像转换成rgb图像图像label2rgb：Convert label matrix into RGB image 将标签图像转换成将标签图像转换成rgb图像。图像。mat2gray：Convert matrix to grayscale image 将数据值（非图像类型）的矩阵转换成灰度图像。将数据值（非图像类型）的矩阵转换成灰度图像。rgb2gray：Convert RGB image or colormap to grayscale 将将rgb图像转换成灰度图像图像转换成灰度图像rgb2ind：

49、Convert RGB image to indexed image 将将rgb图像转换成索引图像。图像转换成索引图像。uint16：Convert data to unsigned 16-bit integers 将数据转换成将数据转换成uint16类型类型uint8：Convert data to unsigned 8-bit integers 将数据转换成将数据转换成uint8类型类型Please turn to help in MATLAB:Image Processing ToolboxFunctions-By CategoryImage Display and Exporation

50、Image Types and Type Conversions.2.7 Converting betwwen Data Classes and Image Types f=-0.5 0.5;0.75 1.5f=-0.5000 0.5000 0.7500 1.5000 g=im2uint8(f)g=0 128 191 255 mat2gray(f)ans=0 0.5000 0.6250 1.0000 A=0 0;1 1A=0 0 1 1 im2double(A)ans=0 0 1 1 lA=logical(A)lA=0 0 1 1 im2double(lA)ans=0 0 1 1 u8A=ui