《数字图像处理MATLAB图像处理.pptx》由会员分享,可在线阅读,更多相关《数字图像处理MATLAB图像处理.pptx(41页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、概述1、MATLAB中的图像和图像处理工具箱MATLAB中大多数是用二维数据(矩阵)double(64位)浮点型,或uint8(8位无符号整数)、uint16(16位无符号整数)来存储的,矩阵中的每一个无素对应图像中的一个像素(pixel).例如:一幅由200行300列不同灰度值的点组成的图像可以用200*300的矩阵来存储;对于彩色图像,如RGB图像,需要用三维数组来存储;第1页/共41页图像文件格式 MATLAB支持下列图像文件格式:bmp(microsoft windows bitmap,位图)hdf(hierarchical data format,层次数据)jpeg(joint ph
2、otographic expert group,静止图像压缩标准)pcx(paintbrush,画刷格式)png(portable network graphics,可移植网络图像)tiff(tagged image file format,标记图像文件)xwd(x window dump)gif(graphics interchange format,图形交换)第2页/共41页MATLAB的图像类型二值图像(Binary):图像的取值为0或1;索引图像(Indexed,pseudocolor):图像为m*n矩阵,矩阵元素的值指向颜色映像表(colormap);灰度图像(intensity,g
3、ray scale,gray level):图为m*n矩阵,矩阵元素的值指的是灰度值。对single or double矩阵取值范围为0,1,对uint8矩阵取值范围为0,255,对uint16矩阵取值范围为0,65535,对uint16矩阵取值范围为-32768,32768;真彩色图像(RGB):图像为m*n*3矩阵;第3页/共41页二值图像第4页/共41页索引图像第5页/共41页灰度图像第6页/共41页真彩色图像第7页/共41页图像的读写与显示常用函数:imread;imwrite;imshow;imagesc;colormap(可以自定义颜色映射矩阵)例如:a=imread(D:MATL
4、ABworkcywaveletlenalena512.BMP);figure,imshow(a,),title(原始图像);或:figure,imagesc(cc);colormap(gray);第8页/共41页MATLAB函数图像类型转换函数gray2ind:Convert a grayscale image to an indexed image;grayslice:Convert a grayscale image,indexed image,or truecolor image,to a binary image;ind2gray:Convert an indexed image to
5、 a grayscale image;ind2rgb:Convert an indexed image to a true color image;mat2gray:Convert a data matrix to a grayscale image,by scaling the data;rgb2gray:Convert a true color image to a grayscale image;rgb2ind:Convert a true color image to an indexed image第9页/共41页数据类型转换函数Im2uint8 rgb2=im2uint8(rgb1
6、)im2uint16im2int16im2singleIm2double第10页/共41页图像的代数运算图像的加法 imadd图像的减法 imsubtract图像的乘法 immultiply图像的除法 imdivide图像的补 imcomplement图像的绝对值的差 imabsdiff图像的线性组合 imlincomb第11页/共41页1 1、加法运算、加法运算主要应用举例:主要应用举例:去除去除“叠加性叠加性”随机噪音随机噪音生成图像叠加效果生成图像叠加效果加法运算第12页/共41页(1)去除去除“叠加性叠加性”噪音噪音 对于原图象对于原图象f(x,y),f(x,y),有一个噪音图像集有一
7、个噪音图像集 g g i i(x,y)i=1,2,.M(x,y)i=1,2,.M其中:其中:g g i i(x,y)=f(x,y)+e(x,y)=f(x,y)+ei(i(x,y)x,y)当:噪音当:噪音ei(x,y)为互不相关,且均值为为互不相关,且均值为0 0时,上述图象均值将降低噪音的影响。时,上述图象均值将降低噪音的影响。M个图像的均值为:第13页/共41页则则 是是 的无偏估计的无偏估计 利用同一景物的多幅图像取平均、消除噪声。取利用同一景物的多幅图像取平均、消除噪声。取MM个图像相加求个图像相加求平均得到平均得到1 1幅新图像,一般选幅新图像,一般选8 8幅取平均。幅取平均。第14页
8、/共41页相加相加Addition:averaging for noise reductionM=1M=2M=4M=16第15页/共41页(2 2)生成图象叠加效果:可以得到各种图像合成的效果,生成图象叠加效果:可以得到各种图像合成的效果,也可以用于两张图片的衔接。也可以用于两张图片的衔接。第16页/共41页主要应用举例:主要应用举例:差影法差影法(检测同一场景两幅图像之间的变化检测同一场景两幅图像之间的变化)混合图像的分离混合图像的分离 将同一景物在不同时间拍摄的图像或同一景物在不同波段将同一景物在不同时间拍摄的图像或同一景物在不同波段的图像相减,这就是图像的减法运算。实际中常称为差影法。的
9、图像相减,这就是图像的减法运算。实际中常称为差影法。差值图像提供了图像间的差值信息,能用于指导动态监测、差值图像提供了图像间的差值信息,能用于指导动态监测、运动目标的检测和跟踪、图像背景的消除及目标识别等。运动目标的检测和跟踪、图像背景的消除及目标识别等。减法运算2 2、减法运算、减法运算第17页/共41页(1)检测同一场景两幅图像之间的变化检测同一场景两幅图像之间的变化设:时刻设:时刻设:时刻设:时刻1 1 1 1的图像为的图像为的图像为的图像为 T T T T1 1 1 1(x,y)(x,y)(x,y)(x,y),时刻时刻时刻时刻2 2 2 2的图像为的图像为的图像为的图像为 T T T
10、T2 2 2 2(x,y)(x,y)(x,y)(x,y)g(x,y)=T g(x,y)=T g(x,y)=T g(x,y)=T2 2 2 2(x,y)-T(x,y)-T(x,y)-T(x,y)-T1 1 1 1(x,y)(x,y)(x,y)(x,y)=-T T1 1(x,y)(x,y)T T2 2(x,y)(x,y)g(x,y)g(x,y)第18页/共41页(2 2)差影法在自动现场监测中的应用差影法在自动现场监测中的应用 1、在银行金库内,摄像头每隔一固定时间拍摄一幅图像,并与上一幅图像做差影,如果图像差别超过了预先设置的阈值,则表明可能有异常情况发生,应自动或以某种方式报警;2、用于遥感图
11、像的动态监测,差值图像可以发现森林火灾、洪水泛滥,监测灾情变化等;3、也可用于监测河口、海岸的泥沙淤积及监视江河、湖泊、海岸等的污染;4、利用差值图像还能鉴别出耕地及不同的作物覆盖情况。第19页/共41页第20页/共41页图:差影法进行混合图像的分离(a)混合图像)混合图像 (b)被减图像)被减图像 (c)差影图)差影图像像(3 3 3 3)混合图像的分离)混合图像的分离)混合图像的分离)混合图像的分离 第21页/共41页(3)消除背景影响:)消除背景影响:即去除不需要的叠加性图案即去除不需要的叠加性图案设:背景图像设:背景图像b(x,y)b(x,y),前景背景混合图像,前景背景混合图像f(x
12、,y)f(x,y)g(x,y)=f(x,y)g(x,y)=f(x,y)b(x,y)b(x,y)g(x,y)g(x,y)为去除了背景图像为去除了背景图像背景背景图像图像差值差值图像图像(a a)从病人头顶向下)从病人头顶向下拍摄的拍摄的X X光照片光照片(b)碘元素注入后拍摄的X光照片与背景图像的差值第22页/共41页乘法运算乘法运算(Multiplication)3、乘法运算乘法运算主要应用举例:主要应用举例:图像的局部显示图像的局部显示改变图像的灰度级改变图像的灰度级第23页/共41页(1)图像的局部显示第24页/共41页乘法乘法运算运算(Multiplication)(a)原图(b)乘以1
13、.2(c)乘以2图:乘法运算结果(2)改变图像的灰度级注意:使用乘法函数时,对uint8图像进行操作往往会发生溢出现象,immultiply函数将溢出的数据截取为数据类型允许的最大值。为了避免这种现象,可以考虑在进行乘法运算前将图像转换为一种数据范围更大的图像类型。第25页/共41页除法除法运算运算(Division)简单的除法运算可用于改变图像的灰度级,简单的除法运算可用于改变图像的灰度级,常用于遥感图像处理中。常用于遥感图像处理中。在四种算术运算中,减法与加法在图像增强在四种算术运算中,减法与加法在图像增强处理中最为有用。处理中最为有用。4、除法运算除法运算第26页/共41页几何变换1)简
14、单变换问题描述:图像的平移、放缩和旋转。解题思路:从易到难。工具:线性代数中的齐次坐标。第27页/共41页几何变换第28页/共41页几何变换第29页/共41页图像的几何操作图像的缩放 imresizeimresize(A,mrows ncols,method)图像缩放,A为缩放因子,method为nearest(默认)(最近邻插值)bilinear(双线性插值)bicubic(三线性插值)图像的旋转 imrotateimrotate(A,angle,method)图像逆时针旋转,angle为角度图像的剪裁 imcropimcrop(A,rect)图像剪切,其中rect为x y width he
15、ight第30页/共41页图像的缩放1)向前映射法通过输入图像像素位置,计算输出图像对应像素位置;将该位置像素的灰度值按某种方式分配到输出图像相邻四个像素.第31页/共41页图像的缩放2)向后映射法通过输出图像像素位置,计算输入图像对应像素位置;根据输入图像相邻四个像素的灰度值计算该位置像素的灰度值.第32页/共41页图像的缩放两种映射方法的对比对于向前映射:每个输出图像的灰度要经过多次运算;对于向后映射:每个输出图像的灰度只要经过一次运算。实际应用中,更经常采用向后映射法。实际应用中,更经常采用向后映射法。其中其中,根据四个相邻像素灰度值计算某位置的像素灰度根据四个相邻像素灰度值计算某位置的
16、像素灰度值即为灰度级插值。值即为灰度级插值。第33页/共41页图像的缩放常用的插值方法:(1)最近邻插值(2)双线性插值(3)三次内插法第34页/共41页图像的缩放(1)最近邻插值向后映射时,输出图像的灰度等于在输入图像中离它所映射位置最近的输入图像的灰度值。(x,y)(x,y+1)(x+1,y)(x+1,y+1)P(x0,y0)通过计算与点P(x0,y0)临近的四个点,并将与点P(x0,y0)最近的整数坐标点(x,y)的灰度值取为P(x0,y0)点灰度近似值。第35页/共41页优点:当P(x0,y0)点各相邻像素间灰度变化较小时,此法是最简单快速的插值方法缺点:当P(x0,y0)点相邻像素间
17、灰度值差异很大时,该灰度估值方法会产生较大的误差,甚至可能影响图像质量。第36页/共41页图像的缩放(2)双线性插值输出像素是它在输入图像中2*2领域采样点的平均值,它根据某像素周围4个像素的灰度值在水平和垂直两个方向上对其插值;四点确定一个平面函数,属于过约束问题;问题描述:单位正方形顶点已知,求正方形内任一点的P(x0,y0)值。具体原理第37页/共41页图像缩放第38页/共41页图像的缩放用最近邻插值和双线性插值的方法分别将老虎放大1.5倍。第39页/共41页图像的缩放采用最近邻插值放大采用最近邻插值放大1.51.5倍倍采用双线性插值放大采用双线性插值放大1.51.5倍倍第40页/共41页感谢您的观看!第41页/共41页