图像处理(一).doc

现代图像处理技术教学进度表华 东 师 范 大 学 教 学 进 度 表教学周 20 周讲 课 27 学时实 验 9 学时机 动 4 学时总 共 40 学时 2010 2011 学年第一学期数学系 数学与应用数学_专业 07_年级 教研室主任（签字） 课程名称：图像处理 任课教师：潘仁良 教材名称（版次）：自编出版社名称： 教材编者：潘仁良 出版时间： 日期周次讲课内容分章和分节的名称课堂时数课外时数习题课、课堂讨论或课堂练习内容课堂时数课外时数实验、实习或科学研究名称月日10288绪论31149图像的运算31110图像的变换31811实验一32512图像的增强312213图像的分析3914实验二31615图像的复原32316图像的形态学操作33017实验三31618复习31319复习上机32020考试3333本课程所用到的基础数学工具非常广泛。例如，线性代数中矩阵分析理论，微分方程中的常微分方程、偏微分方程及稳定性理论，复变函数中的积分变换理论，概率论中的离散随机过程理论，数值分析中的常用计算方法及更进一步人工神经网络、小波分析、形态学、系统学等理论。有些在基础课中未学过的内容在本课程进行过程中作适当补充。但与纯数学课程的演绎推理体系有较大区别，并不着意追求理论上的深度与广度，而重在介绍基本思想与方法。本课程所讲授内容尽量结合计算机仿真，并引入数学软件MATLAB平台中的图像处理工具箱Image Processing Toolbox及与图像处理有关的MATLAB其它工具箱。目的想通过这门课程的学习使学生获得数字图像处理的基本思想、基本方法及应用背景，为能解决实际问题和更深入学习信息处理、图像处理方向的后继课程（模式识别、图像识别、计算机视觉、动漫、虚拟现实等）打下良好基础。增强学生应用数学知识，提取数学问题，解决实际问题的能力。（讲）：工作平台: MATLABImage Processing Toolbox及与图像处理有关的MATLAB其它工具箱一、 引言（或介绍一些与图像处理有关的概念）1图像图（Picture）：原意指图片、图画、各种照片以及各种光照影像等，是物体反射或透视光的分布，是自然界的人、物、景的模拟，它是客观存在的。像（Image）：是客观世界通过光学系统的视觉印象。是人的视觉系统所接受到的图在人脑中所形成的印象或认识，带有主观性。图像是以各种观察系统用不同的形式和手段观察客观世界而获得的，它可以直接或间接地作用于人眼并进而产生视觉印象的实体。图像信息不仅包含光通量的分布，而且还包含人类视觉的主观感受。人类社会已经进入了信息时代，对信息的获取、处理、传输等构成了现代社会的基础性工作。科学研究和统计表明，人类所获得的外界信息约有75% 是通过人的视觉系统，也就是通过图像获得的。图像中带有大量的信息，古人云“百闻不如一见”，就充分说明了这个道理。随着计算机技术和网络技术的空前发展，即计算机计算速度和大规模存储容量以几何级数增长，网络和通信速度的飞速提高，显示系统的逐步成熟，数字图像处理技术已深入到各学科领域，学习和掌握数字图像处理技术已成为信息处理中的一项重要内容。图像分模拟（或物理）图像和数字图像我们只讨论数字图像。2数字图像数字图像处理具有精度高，处理内容丰富、方便，可进行复杂的非线性处理等，且具有非常灵活的变通能力即重复性好，便于传输和保存等。(讲) 图像最一般的数学表达式： 亮度函数 其中是空间坐标，是波长（光谱的波长与色度有关），是时间，而I是坐标点（像素）处的光的强度。上式表示活动的、彩色的、三维的视频图像。特别，对于静止图像，则与时间无关；对于单色图像，则波长是常数；对于平面图像，则与坐标无关。对于平面上静止的单色图像，其则数学表达式可简化为： 特别如果是彩色图像，则是以三基色（RGB）的明亮度作为分量的二维矢量函数来表示。即数字图像处理和分析系统：数字图像处理系统基本的三个部件是：处理图像的计算机、图像数字化仪、图像显示设备。（1） 处理图像的计算机包括计算机硬件和处理图像的计算机软件。处理图像的计算机软件有Visual C+、Photoshop、3D Studio MAX等。例如Photoshop，该软件具有良好的人机交互界面，操作简单易学，能够实现图像的显示、几何运算、着色、增亮、表面材质、增强、滤波等一系列的操作，效果非常不错。且该软件生成的图像文件格式目前已被国际标准予以应用。但图像处理Photoshop。事实上，Photoshop只是一种通用的图像处理软件，一般是用来提高视觉感受和效果的，而实际的图像处理工作主要是针对不同的应用领域以提取图像的不同信息加以研究，这是Photoshop不擅长甚至于不支持的。针对不同的应用领域的不同问题，需要不同的图像处理算法，MATLAB(讲)MATLAB的Image Processing Toolbox及与图像处理有关的MATLAB其它工具箱（2）图像数字化仪包括数码照相机、数码摄像机、扫描仪等。计算机只能接受和处理离散数据，连续信号必须离散化才能被计算机接受和处理。连续信号，以为采样间隔或周期(通常记)对进行采样，得到一个数列（时间序列）： 或令=1（一个单位）， 或从用数学表示：其中为广义的单位脉冲函数，且，且。称为移动的单位脉冲函数。对一个连续信号进行采样也就是该连续信号与一个等间隔的单位脉冲序列共同作用的结果。因为只有在时刻才能被采样其中为在时刻的值。任何一个序列可看作单位脉冲的移动加权和。推广到二维情况：有 对模拟信号采样使它离散化的程度取决于采样周期T，即T愈小，离散序列愈接近于原来的信号，但离散序列的数据量就愈大。反之T愈大，离散序列与原来的信号接近的程度就差，但离散序列的数据量就小。并且当采样周期T大到超过某个限制时，就会使离散序列和原来的信号之间产生不可逆转的误差。采样周期T的这个限制由采样定理来确定。采样定理：对一个频谱有限的（）连续信号进行采样，当采样频率满足条件时，采样函数便能无失真地恢复原来的连续信号。为信号的有效频谱的最大角频率，为采样频率，一般来说图像信号是二维空间的信号，通常记为，例如黑白图像或灰度值图像表示一幅图像在水平和垂直两个方向上的光照强度的变化。一幅图像必须要在空间和灰度上都离散化才能被计算机接受和处理。这种离散化了的图像就是数字图像，相应的过程称为图像的数字化。图像在二维空间采样时常用的方法是对进行均匀采样，取得各点的亮度值，构成一个离散函数。如果是彩色图像，则是以三基色（RGB）的明亮度作为分量的二维矢量函数来表示。即相应的离散值为例1：图像的RGB分解clear all,close allRGB=imread(C18.bmp);subplot(2,2,1);imshow(RGB);R(:,:,1)=RGB(:,:,1);R(:,:,2)=0;R(:,:,3)=0;subplot(2,2,2);imshow(R);title(R分量);G(:,:,2)=RGB(:,:,2);G(:,:,1)=0;G(:,:,3)=0;subplot(2,2,3);imshow(G);title(G分量);B(:,:,3)=RGB(:,:,3);B(:,:,1)=0;B(:,:,2)=0;subplot(2,2,4);imshow(B);title(B分量);同一维信号一样，图像信号的采样也要遵循采样定理。二维信号采样定理与一维信号采样定理类似。图像信息的采样定理：对一个频谱有限的（）图像信号进行采样，当采样频率满足，条件时，采样函数便能无失真地恢复原来的连续信号。分别为信号在两个方向的频域上的有效频谱的最大角频率；分别为二维采样频率，。模拟（或物理）图像扫描、采样、量化数字化表格或矩阵可在MATLAB平台上操作。数字图像就是经过采样和量化后的图像。由于原始图像是可以用亮度函数表示的二维函数，所以，数字图像就是被采样和量化后的二维函数。通常是用矩形网格采样并对图像像素点的亮度幅值进行均匀量化实现的，所以数字图像表现为一个二维矩阵。一幅图像由许许多多的像素点构成，每个像素点包含着反映图像在读点的明暗和颜色变化等信息，这种图像叫做位图图像（bitmap）。（3）图像显示设备包括显示器（荧屏显示）、打印机、传真机等。3图像处理技术的发展从广义上说，图像技术是各种与图像有关的技术的总称目前人们主要研究的是数字图像，通常结合计算机和其它电子设备来完成一系列工作，例如图像的采集、获取、编码、存储和传输，图像的合成和产生，图像的显示和输出，图像的变换、增强、恢复（复原）和重建，图像的分割，目标的检测、表达和描述，特征的提取和测量，序列图像的校正，3-D景物的重建复原，图像数据库的建立、索引和抽取，图像的分类、表示和识别，图像模型的建立和匹配，图像与场景的解释和理解，以及基于它们的判断决策和行为规划等等另外，图像技术还可包括为完成上述功能而进行的硬件设计及制作等方面的技术。 数字图像处理的历史可以追溯到20世纪20年代报纸业中采用的Bartlane电缆图片传输系统。当时把一幅照片的数据横跨大西洋传输需要一个星期，而Bartlane系统通过编码技术把时间减少到了3个小时。 精确地说，上述例子引用的图像处理技术并没有考虑到计算机技术，不能算是真正意义上的数字图像处理。尽管计算机技术的历史可追溯到1946年世界上第1台电子计算机的诞生，但在20世纪50年代，计算机主要还是用于数值计算，满足不了处理大数据量对象（如图像数据量）的要求。 20世纪60年代，第3代计算机的研制成功以及快速傅里叶变换算法的发现和应用使得对图像的某些计算得以实现，人们逐步开始利用计算机对图像进行加工处理真正意义上的数字图像技术的诞生是这一时期图像处理技术在空间开发项目上的应用。比较著名的例子就是1964年科学家利用计算机对“旅行者7号”发回的第一张月球照片作几何校正。与此同时，数字图像处理技术在医学图像、地球遥感检测和天文学等领域也逐渐得到广泛应用。 20世纪70年代，图像技术有了长足的进展，并且第一本重要的图像处理专著 Digital Picture Processing Rosenfeld，1976也得以出版在这个时期发明的计算机轴向断层技术（也称为计算机断层或CT）是图像处理在医学领域最重要的应用。从那个时期开始，医学图像处理成为图像处理最活跃的应用领域之一。 20世纪80年代，各种硬件的发展使得人们不仅能处理静态图像，而且可以处理动态图像。许多能获取动态图像的设备和处理分析动态图像的系统被研制成功，图像技术得到了广泛的应用。除了医学和空间项目的应用外，在地球物理学领域，图像处理技术被用来从航空和卫星传送回来的图片中研究污染模型。图像增强和复原技术用于不可修复物体的已损图像在考古领域，使用图像处理方法已经成功地复原了模糊的图片。在物理学的一些领域，图像处理技术也用于增强高能等离子和电子显微镜等产生的实验图像。此外，图像处理技术也成功地应用于天文学、生物学、核医学、国防及工业领域。 进入20世纪90年代，图像技术已逐步渗透到人类生活和社会发展的各个方面。比较有代表性的有多媒体（multi-media）技术、虚拟现实（VR）技术、动漫技术等。一个越来越值得注意的事实是，如今的图像处理技术在不断地与网络技术，特别是互联网技术相融合。 上面所叙述的应用中，图像处理的结果主要是作为人类进行图像分析的辅助手段。同样，数字图像技术也可以用来解决机器感知问题，这是一个用计算机处理的形式从图像中提取信息的过程。在机器感知中，使用图像处理技术的典型问题有自动字符识别、军事识别、指纹的自动处理、用于生产线检测和公寓检测机器视觉、用于天气预报和环境鉴定的航空与卫星图像的机器处理等。4数字图像处理研究的主要内容：（1） 图像数据的变换、编码和压缩，以便于图像的存储和传输。（2） 提高图像的视感质量，如进行图像的亮度和色彩变换，增强或抑制图像的某些成份，对图像进行几何变换等，以改善图像的视觉效果。（3） 提取图像中所包含的某些特征或某些特殊信息，如图像的频率特征、灰度或颜色特征、边界或区域特征、纹理或形态特征、拓扑结构和关系结构等。数字图像处理概括的说主要包括图像的转换和存储、图像的视觉优化和图像的理解三个层次。描述数字图像的基本参数有三个，即：图像分辨率、图像深度和图像数据容量。 1)图像分辨率 图像分辨率（resolutio）实际上是指对原始图像的采样分辨率，也就是指图像水平或垂直方向单位长度上所包含的采样点数。单位是“像素点单位长度”，例如：像素点毫米（ Pixelmm）、像素点英寸（ pixelinch）。 为了方便，在描述显示器的图像分辨率时，往往简单地以显示屏的有效范围内水平方向（扫描线）上的像素点总数和垂直方向上的扫描线总数的乘积来界定。例如， 的VGA显示器的分辨率在水平方向为640个像素点，垂直方向为480个像素点。 2)图像深度 在位图图像中，各像素点的亮度或色彩信息用二进制数位来表示，这一数据位的位数称为像素深度，也叫图像深度（ image deepness）。 例如，深度为1位的图像，当取值为0时表示黑（暗）色，取值为1时表示白（亮）色，也就是一幅黑白（二值）图像。如果一个像素点用4位二进制数位表示，若是彩色图像，每个4位的像素点就可以表示16种颜色中的任何一种颜色；如果不是彩色图像，每个4位的像素点就可以表示 16个不同黑白程度的灰度级别之一。常用的图像深度是8位，这时，每个像素点可表示256种颜色（256级灰度），16位可表示65536种颜色， 24位可表示16777216种颜色。图像深度愈深，能够表现的颜色的数量愈多，图像的色彩也就愈真实。 3）图像数据量图像数据量是一幅图像的总像素点数目与每个像素点所需字节数的乘积。注：实际上人眼只能分辨40级灰度，这意味着黑白之间的灰度范围被分为超过40个以上的等级，则相邻的灰度级对人眼来说看上去是相同的。常用图像格式：数字图像处理通常会产生大量的含有数字图像的文件，而且文件的数据量相当大。它们必须被存档，而且经常需要在不同的用户以及系统间进行交换。这就要求有一些用于数字图像文件存储和交换的标准图像文件格式。图像格式是多种多样的，有未经压缩的格式（如位图BMP格式）和压缩的格式（如GIF格式）。主要的图像数据文件的格式如下所示：bit-mapped format .BMP Microsoft Wndows格式Isoft image file fomat .PCX 通用的标准图像文件格式graphical interchange format .GIF CompuServe图形格式tagged image file fomat .TIF 扫描仪生成的图像格式Trucvision fomat .TGA 国际上的图形、图像工业广泛采用的格式Macintosh format .PICT Apple Macintosh图像格式Photoshop fomat .PSD Photoshop图像格式Portable network graphics format .PNG 新兴的Web图像格式JPEG file interchanse format JPG 静止图像压缩国际标准格式等而MATLAB图像的格式或类型：MATLAB图像处理工具箱三种不同的的数据存储类型：双精度浮点（64位）8位无符号整型数据（unit8）16位无符号整型数据（unit16）注：MATLAB主包中只有一种数据存储类型即双精度浮点（64位）若unit8、unit16的图像要在MATLAB主包中进行运算，必须用im2double函数（或语句）把unit8、unit16的图像化为双精度浮点的图像（例如下面例2中的程序）。且有四种不同的图像类型：（1） 二值图像（）每个像素取两个离散值0、1中的一个，0-1矩阵。可以是unit8，也可以是双精度的，但不可以是unit16。（2） 索引图像（X，map）包括一个数据矩阵，一个调色板矩阵，其中为调色板中颜色的数目为。对应调色板矩阵中第行，。总是双精度的，可以是unit8，可以是unit16，也可以是双精度的。（3） 灰度图像（）可以是unit8，可以是unit16，也可以是双精度的。（4） RGB图像（）可以是unit8，可以是unit16，也可以是双精度的。不同的图像类型之间可互相转换。（讲）MATLAB图像的读写和显示：读：imread 外内写：imwrite内外显示：MATLAB内部显示 imshow（讲）5MATLAB图像处理初步例2: 图像的反转显示clear all,close allI=imread(I2.bmp);I1=im2double(I); J=255-I1; %或J=imadjust(I,0 1,1 0);subplot(121), imshow(I,);subplot(122), imshow(J,) 图像的反转类似于照片的底片，将灰度值为i的像素点按255-i(二值图像按1-i)显示。J就是得到的反转图像。图中同时显示了原图像和反转图像，可以作比较。例3: 图像暴光不足(或暴光太足或图像灰度值过于集中) 的例子clear all,close allI=imread(pout.tif);J=histeq(I);subplot(221),imshow(I)subplot(222),imshow(J)subplot(223),imhist(I)subplot(224),imhist(J)为了便于观察直方图均衡化前后图像的对比度的变化，先产生一幅质量较差（偏暗或者偏亮或灰度值集中）的图像，均衡化后，原图像的对比度得到了增强，原图像的直方图相应地变得平坦，灰度值均匀地分布在0255。例4: 二维低通FIR滤波器设计（解释图像滤波为什么用FIR滤波器）一维FIR滤波器设计有：fir1(或fir2) 窗口函数法remez 等波纹最佳一致逼近firls 最小二乘逼近等频率变换方法：首先设计一个性能较好的一维FIR滤波器