Matlab图像处理教程及方法课件

Matlab数字图像处理1ppt课件Matlab数字图像处理1ppt课件/1、图像的读取和显示/2、图像的点运算/4、空间域图像增强/5、频率域图像增强/6、彩色图像处理/7、形态学图像处理/8、图像分割/9、特征提取内容安排/3、图像的几何变换2ppt课件/1、图像的读取和显示/2、图像的点运算/4、空间域图像增强/1、图像的读取和显示一、图像的读取A=imread(FILENAME,FMT)FILENAME指定图像文件的完整路径和文件名。如果在work工作目录下只需提供文件名。FMT为图像文件的格式对应的标准扩展名。I_1=imread('D:\10.06.08nir\TTC10377.BMP');%读入图像二、图像的写入imwrite(A，FILENAME,FMT)FILENAME参数指定文件名。FMT为保存文件采用的格式。imwrite(I6,'nirdilatedisk2TTC10373.bmp');3ppt课件/1、图像的读取和显示一、图像的读取A=imread(FIL/1、图像的读取和显示三、图像的显示imshow(I,[lowhigh])I为要显示的图像矩阵。[lowhigh]为指定显示灰度图像的灰度范围。高于high的像素被显示成白色；低于low的像素被显示成黑色；介于High和low之间的像素被按比例拉伸后显示为各种等级的灰色。

figure；imshow(I6);title('TheMainPassPartofTTC10373');figure；%创建一个新的窗口figure；subplot（m,n,p）；imshow（I）；Subplot(m,n,p)含义为：打开一个有m行n列图像位置的窗口，并将焦点位于第p个位置上。4ppt课件/1、图像的读取和显示三、图像的显示imshow(I,[lo/1、图像的读取和显示四、图像的格式转换im2bw(I,LEVEL);rgb2gray;从RGB图创建灰度图，存储类型不变。im2uint8将图像转换成uint8类型阈值法从灰度图、RGB图创建二值图。LEVEL为指定的阈值；（0,1）。im2double将图像转换成double类型5ppt课件/1、图像的读取和显示四、图像的格式转换im2bw(I,LE/2、图像的点运算灰度直方图描述了一副图像的灰度级统计信息，主要应用于图像分割和图像灰度变换等处理过程中。从数学角度来说，图像直方图描述图像各个灰度级的统计特性，它是图像灰度值的函数，统计一幅图像中各个灰度级出现的次数或概率。归一化直方图可以直接反映不同灰度级出现的比率。横坐标为图像中各个像素点的灰度级别，纵坐标表示具有各个灰度级别的像素在图像中出现的次数或概率。imhist(I);%灰度直方图I=imread(‘red.bmp’);%读入图像

figure;%打开新窗口[M,N]=size(I);%计算图像大小[counts,x]=imhist(I,32);%计算有32个小区间的灰度直方图counts=counts/M/N;%计算归一化灰度直方图各区间的值stem(x,counts);%绘制归一化直方图一、图像直方图图像直方图归一化6ppt课件/2、图像的点运算灰度直方图描述了一副图像7ppt课件7ppt课件/2、图像的点运算二、灰度的线性变换Fa>1时，输出图像的对比度将增大；Fa<1时，输出图像对比度将减小。Fa=1且Fb非零时，所有像素的灰度值上移或下移，使整个图像更暗或更亮。Fa<0，暗区变亮，亮区变暗。8ppt课件/2、图像的点运算二、灰度的线性变换Fa>1时，输出图像的对9ppt课件9ppt课件10ppt课件10ppt课件/2、图像的点运算三、灰度的对数变换c为尺度比例常数，s为源灰度值，t为变换后的目标灰度值。k为常数。灰度的对数变换可以增强一幅图像中较暗部分的细节，可用来扩展被压缩的高值图像中的较暗像素。广泛应用于频谱图像的显示中。Warning：log函数会对输入图像矩阵s中的每个元素进行操作，但仅能处理double类型的矩阵。而从图像文件中得到的图像矩阵大多是uint8类型的，故需先进行im2double数据类型转换。11ppt课件/2、图像的点运算三、灰度的对数变换c为尺度比例常数，s为源I=imread('nir.bmp');%读入图像F=fft2(im2double(I));%FFTF=fftshift(F);%FFT频谱平移F=abs(F);T=log(F+1);%频谱对数变换figure;imshow(F,[]);title('未经变换的频谱');figure;imshow(T,[]);title('对数变换后');12ppt课件I=imread('nir.bmp');%读入图像12ppt/2、图像的点运算三、灰度的Gamma变换其中，x、y的取值范围为[0,1]。esp为补偿系数，r则为Gamma系数。Gamma变换是根据r的不同取值选择性的增强低灰度区域的对比度或者高灰度区域的对比度。13ppt课件/2、图像的点运算三、灰度的Gamma变换其中，x、y的取值J=imadjust(I,[low_inhigh_in],[low_outhigh_out],gamma)I=imread('nir.bmp');figure;imshow(imadjust(I,[],[],0.75));%gamma=0.5title('Gamma0.5');原NIR图像Gamma0.514ppt课件J=imadjust(I,[low_inhigh_in四、灰度阈值变换及二值化Gamma1.5T为指定阈值BW=im2bw(I,level);%level为人工设定阈值范围为[0,1]thresh=graythresh(I);%自动设定所需的最优化阈值15ppt课件四、灰度阈值变换及二值化Gamma1.5T为指定阈值BWOTSU算法：最大类间方差法自动单阈值分割。Kapur算法：一维直方图熵阈值算法niblack算法：局部阈值分割阈值的计算公式是T=m+k*v，其中m为以该像素点为中心的区域的平均灰度值，v是该区域的标准差，k是一个系数。kittlerMet:表示kittler最小分类错误(minimumerrorthresholding)全局二值化算法。

原图像直接阈值分割0.2516ppt课件OTSU算法：最大类间方差法自动单阈值分割。Kapur算法：自动阈值分割OTSU算法分割KittlerMet算法Niblack算法Kapur算法17ppt课件自动阈值分割OTSU算法分割KittlerMet算法Nibl/2、图像的点运算五、直方图均衡化[J,T]=histeq(I);%J为输出图像，T为变换矩阵图像易受光照、视角、方位、噪声等的影响。使得同一类图像的不同变形体之间的差距有时大于该类图像与另一类图像之间的差距，影响图像识别、分类。图像归一化就是将图像转换到唯一的标准形式以抵抗各种变换，从而消除同类图像不同变形体之间的外观差异。也称为图像灰度归一化。原图像及直方图18ppt课件/2、图像的点运算五、直方图均衡化[J,T]=histeq(图像变亮后灰度均衡化图像变暗后灰度均衡化19ppt课件图像变亮后灰度均衡化图像变暗后灰度均衡化19ppt课件/3、图像的几何变换一、图像平移正变换逆变换strel%用来创建形态学结构元素translate(SE,[yx])%原结构元素SE上y和x方向平移imdilate%形态学膨胀20ppt课件/3、图像的几何变换一、图像平移正变换逆变换20ppt课件/3、图像的几何变换I=imread('nir.bmp');se=translate(strel(1),[180190]);B=imdilate(I,se);figure;subplot(1,2,1),subimage(I);title('原图像');subplot(1,2,2),subimage(B);title('平移后图像');21ppt课件/3、图像的几何变换I=imread('nir.bmp');/3、图像的几何变换二、图像镜像B=imtransform(A,TFORM,method);TFORM=makeform(transformtype,Matrix);%空间变换结构参数transformtype指定了变换的类型，常见的’affine’为二维或多维仿射变换，包括平移、旋转、比例、拉伸和错切等。Matrix为相应的仿射变换矩阵。22ppt课件/3、图像的几何变换二、图像镜像参数transformtypA=imread('nir.bmp');[height,width,dim]=size(A);tform=maketform('affine',[-100;010;width01]);B=imtransform(A,tform,'nearest');tform2=maketform('affine',[100;0-10;0height1]);C=imtransform(A,tform2,'nearest');figure;imshow(A);figure;imshow(B);imwrite(B,'nir水平镜像.bmp');figure;imshow(C);imwrite(B,'nir垂直镜像.bmp');原图像水平镜像图像垂直镜像图像23ppt课件A=imread('nir.bmp');原图像水平镜像图像垂A=imread('nir.bmp');tform=maketform('affine',[010;100;001]);B=imtransform(A,tform,'nearest');figure;imshow(A);figure;imshow(B);imwrite(B,'nir转置后图像.bmp');/3、图像的几何变换三、图像转置24ppt课件A=imread('nir.bmp');/3、图像的几何变换/3、图像的几何变换四、图像中心旋转B=imrotate(A,angle,method,’crop’);angle为旋转角度，正值为逆时针旋转。可选参数method为imrotate函数指定插值方法。‘crop’选项会裁减旋转后增大的图像，保持和原图像同样大小。

A=imread('nir.bmp');B=imrotate(A,30,'nearest','crop');figure;imshow(B);imwrite(B,'逆时针中心旋转30度.bmp');逆时针30度25ppt课件/3、图像的几何变换四、图像中心旋转angle为旋转角度，正/4、空间域图像增强一、噪声添加h=imnoise(I,type,parameters);type为噪声类型，合法值如下：添加高斯白噪声添加椒盐噪声26ppt课件/4、空间域图像增强一、噪声添加h=imnoise(I,ty/4、空间域图像增强二、空间域滤波滤波过程就是在图像f(x,y)中逐点移动模板，使模板中心和点(x,y)重合，滤波器在每一点(x,y)的响应是根据模板的具体内容并通过预先定义的关系来计算的。27ppt课件/4、空间域图像增强二、空间域滤波滤波过程就是在图像f(x,B=imfilter(f,w,option1,option2,…);f为要进行滤波操作的图像。w为滤波操作使用的模板，为一个二维数组，可自己定义。option1……是可选项，包括：1、边界选项（’symmetric’、’replicate’、’circular’）2、尺寸选项（’same’、’full’）3、模式选项（’corr’、’conv’）原图像滤波后图像28ppt课件f为要进行滤波操作的图像。原图像滤波后图像28ppt课件三、滤波器设计/4、空间域图像增强h=fspecial(type,parameters)parameters为可选项，是和所选定的滤波器类型type相关的配置参数，如尺寸和标准差等。type为滤波器的类型。其合法值如下：29ppt课件三、滤波器设计/4、空间域图像增强h=fspecial(ty/4、空间域图像增强四、中值滤波h=medfilt2(I1,[m,n]);m和n为中值滤波处理的模板大小，默认3*3中值滤波本质上是一种统计排序滤波器。中值不同于均值，是指排序队列中位于中间位置的元素的值。中值滤波并非线性滤波器。对于某些类型的随机噪声具有非常理想的降噪能力。典型的应用就是消除椒盐噪声。30ppt课件/4、空间域图像增强四、中值滤波h=medfilt2(I1,/4、空间域图像增强五、图像锐化图像锐化主要用于增强图像的灰度跳变部分，主要通过运算导数（梯度）或有限差分来实现。主要方法有：Robert交叉梯度，Sobel梯度，拉普拉斯算子，高提升滤波，高斯-拉普拉斯变换。Robert交叉梯度w1对接近正45°边缘有较强响应，w2对接近负45°边缘有较强响应。31ppt课件/4、空间域图像增强五、图像锐化图像锐化主要用于增强图像的灰Sobel交叉梯度对水平边缘有较大响应对垂直边缘有较大响应拉普拉斯算子32ppt课件Sobel交叉梯度对水平边缘有较大响应对垂直边缘有较大响应拉I=imread('nir.bmp');I=double(I);%双精度化w1=[-10;01];w2=[0-1;10];G1=imfilter(I,w1,‘corr’,‘replicate’);%正45°梯度G2=imfilter(I,w2,‘corr’,‘replicate’);%负45°梯度G=abs(G1)+abs(G2);%计算Robert梯度figure;imshow(G,[]);figure;imshow(abs(G1),[]);figure;imshow(abs(G2),[]);Robert交叉梯度I=imread('nir.bmp');Id=double(I);%双精度化h_1=fspecial('log',5,0.5);%大小为5，sigma=0.5的LOG算子I_1=imfilter(Id,h_1,'corr','replicate');figure;imshow(uint8(abs(I_1)),[]);h_2=fspecial('log',5,2);%大小为5，sigma=2的LOG算子I_2=imfilter(Id,h_2,'corr','replicate');figure;imshow(uint8(abs(I_2)),[]);高斯-拉普拉斯锐化33ppt课件I=imread('nir.bmp');Robert交叉梯度/5、频率域图像增强用傅里叶变换表示的函数特征可以完全通过傅里叶反变换进行重建而不丢失任何信息吉布斯现象Gibbsphenomenon（又叫吉布斯效应）：将具有不连续点的周期函数（如矩形脉冲)进行傅立叶级数展开后，选取有限项进行合成。当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时，该峰起值趋于一个常数，大约等于总跳变值的9%。这种现象称为吉布斯现象。34ppt课件/5、频率域图像增强用傅里叶变换表示的函数特征可以完全通过傅幅度谱相位谱功率谱35ppt课件幅度谱相位谱功率谱35ppt课件一、傅里叶变换/5、频率域图像增强I=fft2(x);%快速傅里叶变换I=fft2(x,m,n);x为输入图像；m和n分别用于将x的第一和第二维规整到指定的长度。当m和n均为2的整数次幂时算法的执行速度要比m和n均为素数时快。I1=abs(I);%计算I的幅度谱I2=angle(I);%计算I的相位谱Y=fftshift(I);%频谱平移I=ifft2(x);%快速傅里叶逆变换I=ifft2(x,m,n);36ppt课件一、傅里叶变换/5、频率域图像增强I=fft2(x);%快速幅度谱相位谱平移后的相位谱37ppt课件幅度谱相位谱平移后的相位谱37ppt课件functionout=imidealflpf(I,freq)%imidealflpf函数构造高斯频域低通滤波器[M,N]=size(I);out=ones(M,N);fori=1:Mforj=1:Nif(sqrt(((i-M/2)^2+(j-N/2)^2))>freq)out(i,j)=0;endendend理想低通functionout=imgaussflpf(I,sigma)%imgaussflpf函数构造高斯频域低通滤波器[M,N]=size(I);out=ones(M,N);fori=1:Mforj=1:Nout(i,j)=exp(-((i-M/2)^2+(j-N/2)^2)/2/sigma^2);endend高斯低通38ppt课件functionout=imidealflpf(I,f/6、彩色图像处理色彩模型：RGB模型、CMY模型、CMYK模型、HIS模型、HSV模型、YUV模型、YIQ模型。RGB模型国际照明委员会（CIE）规定以蓝（435.8nm）、绿（546.1nm）和红（700nm）作为主原色。39ppt课件/6、彩色图像处理色彩模型：RGB模型、CMY模型、CMYKRGB_image=cat(3,PR,PG,PB);%将PR、PG、PB三个矩阵在第3个维度上进行级联，进行图像合成PR=RGB_image(:,:,1);%提取红色分量PG=RGB_image(:,:,2);%提取绿色分量PB=RGB_image(:,:,3);%提取蓝色分量Matlab中一幅RGB图像可表示为一个M*N*3的3维矩阵。其中每一个彩色像素都在特定空间位置的彩色图像中对应红、绿、蓝3个分量。CMY模型CMY模型是采用（Cyan、Magenta、Yellow）青、品红、黄色3中基本原色按一定比例合成颜色。由于色彩的显示是由光线被物体吸收掉一部分之后反射回来的剩余光线产生，故CMY模型又称为减色法混色模型。当光都被吸收时成为黑色，都被反射时为白色。40ppt课件RGB_image=cat(3,PR,PG,PB);%将PRCMY模型主要用于彩色打印机和复印机等。CMYK模型CMY模型在实际使用中，青、品红和黄色等比例混合后的黑色并不纯，为产生真正的黑色，专门加入第四种颜色—黑色。得到CMYK模型。用于四色打印。cmy=imcomplement(rgb);%rgb转成cmyrgb=imcomplement(cmy);cmy转成rgb41ppt课件CMY模型主要用于彩色打印机和CMYK模型CMY模型在实际使HSI模型HIS模型是从人的视觉系统出发，直接使用颜色三要素色调（Hue）、饱和度（Saturation）和亮度（Intensity）来描述颜色。亮度指人眼感觉光的明暗程度。光的能量越大，亮度越大。色调由物体反射光线中占优势的波长决定。反映颜色的本质。饱和度指颜色的深浅和浓淡程度，饱和度越高，颜色越深。HIS色彩空间比RGB彩色空间更符合人的视觉特性。亮度和色度具有可分离特性，使得图像处理和机器视觉中大量灰度处理算法都可在HIS彩色空间中方便使用。42ppt课件HSI模型HIS模型是从人的视觉系统出发，直接使用颜色三要素色调强度饱和度其中43ppt课件色调强度饱和度其中43ppt课件/7、形态学图像处理一、二值图像腐蚀I2=imerode(I,SE);SE=strel(shape,parameters);I为原始图像，可以是二值或者灰度图像。shape指定了结构元素的形状。parameters是和输入shape有关的参数。44ppt课件/7、形态学图像处理一、二值图像腐蚀I2=imerode(I/7、形态学图像处理二、二值图像膨胀I2=imdilate(I,SE);SE=strel(shape,parameters);I为原始图像，可以是二值或者灰度图像。shape指定了结构元素的形状。parameters是和输入shape有关的参数。原图像腐蚀后膨胀后45ppt课件/7、形态学图像处理二、二值图像膨胀I2=imdilate(/7、形态学图像处理三、其他二值图像运算SE=strel(shape,parameters);I2=imopen(I,SE);%开运算I3=imclose(I,SE);%闭运算Ihm=bwhitmiss(I,SE1,SE2);%击中击不中变换四、连通分量提取[Lnum]=bwlabel(Ibw,conn);Ibw为一幅输入二值图像。conn为可选参数，指明提取连通分量是4连通还是8连通。默认为8。L为连通分量标注图像。num为二值图像Ibw中连通分量个数。46ppt课件/7、形态学图像处理三、其他二值图像运算SE=strel(s/7、形态学图像处理五、bwmorph函数Iout=bwmorph(I,operation,n)47ppt课件/7、形态学图像处理五、bwmorph函数Iout=bwmo/7、形态学图像处理六、顶帽变换解决非均匀光照问题。顶帽变换定义为图像f与图像f的开运算之差。a=imread('red.bmp');I=rgb2gray(a);figure,surf(double(I(1:8:end,1:8:end))),zlim([0255])，colormapgray;%显示图像I的三维可视化效果bg=imopen(I,strel(‘disk’,15));%开运算Itophat=imsubtract(I,bg);%顶帽变换figure,imshow(Itophat);figure,surf(double(Itophat(1:8:end,1:8:end))),zlim([0255]);I2=imadjust(Itophat);figure,imshow(I2);48ppt课件/7、形态学图像处理六、顶帽变换解决非均匀光照问题。顶帽变换原图像三维可视化效果顶帽变换后的三维可视化效果顶帽变换后图像对比度拉伸后效果49ppt课件原图像三维可视化效果顶帽变换后的三维可视化效果顶帽变换后图像/8、图像分割一、图像分割概述图像分割一般采用的方法有边缘检测（edgedetection）、边界跟踪（edgetracing）、区域生长（regiongrowing）、区域分离和聚合等。

图像分割算法一般基于图像灰度值的不连续性或其相似性。

不连续性是基于图像灰度的不连续变化分割图像，如针对图像的边缘有边缘检测、边界跟踪等算法。

相似性是依据事先制定的准则将图像分割为相似的区域，如阈值分割、区域生长等。二、边缘检测图像的边缘点是指图像中周围像素灰度有阶跃变化或屋顶变化的那些像素点，即灰度值导数较大或极大的地方。

边缘检测可以大幅度的减少数据量，并且剔除不相关信息，保留图像重要的结构属性。

边缘检测基本步骤：平滑滤波、锐化滤波、边缘判定、边缘连接。50ppt课件/8、图像分割一、图像分割概述图像分割一般采边缘检测算法：基于一阶导数：Roberts算子、Sobel算子、Prewitt算子基于二阶导数：高斯-拉普拉斯边缘检测算子Canny边缘检测算法Matlab实现1）、基于梯度算子的边缘检测BW=edge(I,type,thresh,direction,’nothinning’)thresh是敏感度阈值参数，任何灰度值低于此阈值的边缘将不会被检测到。默认值为空矩阵[]，此时算法自动计算阈值。51ppt课件边缘检测算法：Matlab实现1）、基于梯度算子的边缘检测Bdirection指定了我们感兴趣的边缘方向，edge函数将只检测direction中指定方向的边缘，其合法值如下：可选参数’nothinning’，指定时可以通过跳过边缘细化算法来加快算法运行的速度。默认是’thinning’，即进行边缘细化。2）、基于高斯-拉普拉斯算子的边缘检测BW=edge(I,’log’,thresh,sigma)sigma指定生成高斯滤波器所使用的标准差。默认时，标准差为2。滤镜大小n*n，n的计算方法为：n=ceil(sigma*3)*2+1。52ppt课件direction指定了我们感兴趣的边缘方向，edge函数将3）、基于Canny算子的边缘检测BW=edge(I,’canny’,thresh,sigma)thresh是敏感度阈值参数，默认值为空矩阵[]。此处为一列向量，为算法指定阈值的上下限。第一个元素为阈值下限，第二个元素为阈值上限。如果只指定一个阈值元素，则默认此元素为阈值上限，其0.4倍的值作为阈值下限。如阈值参数没有指定，则算法自行确定敏感度阈值上下限。b1=imread('nir.bmp');h58=fspecial('gaussian',5,0.8);b=imfilter(b1,h58);bw1=edge(b,‘sobel’);%sobel算子bw2=edge(b,‘prewitt’);%prewitt算子bw3=edge(b,‘roberts’);%roberts算子bw4=edge(b,‘log’);%log算子bw5=edge(b,‘canny’);%canny算子figure;imshow(bw1);imwrite(bw1,'nirbwsobel.bmp');figure;imshow(bw2);imwrite(bw2,'nirbwprewitt.bmp');figure;imshow(bw3);imwrite(bw3,'nirbwroberts.bmp');figure;imshow(bw4);imwrite(bw4,'nirbwlog.bmp');figure;imshow(bw5);imwrite(bw5,'nirbwcanny.bmp');53ppt课件3）、基于Canny算子的边缘检测BW=edge(I,’caSobel算子roberts算子prewitt算子log算子canny算子54ppt课件Sobel算子roberts算子prewitt算子log算子分析1、边缘定位精度方面：Roberts算子和Log算子定位精度较高。Roberts算子简单直观，Log算子利用二阶导数零交叉特性检测边缘。但Log算子只能获得边缘位置信息，不能得到边缘方向信息。2、边缘方向的敏感性：Sobel算子、Prewitt算子检测斜向阶跃边缘效果较好，Roberts算子检测水平和垂直边缘效果较好。Log算子不具有边缘方向检测功能。Sobel算子能提供最精确的边缘方向估计。3、去噪能力：Roberts算子和Log算子虽然定位精度高，但受噪声影响大。Sobel算子和Prewitt算子模板相对较大因而去噪能力较强，具有平滑作用，能滤除一些噪声，去掉一部分伪边缘，但同时也平滑了真正的边缘，降低了其边缘定位精度。总体来讲，Canny算子边缘定位精确性和抗噪声能力效果较好，是一个折中方案。

55ppt课件分析1、边缘定位精度方面：55ppt课件三、Hough变换直线检测/8、图像分割步骤：利用hough（）函数执行霍夫变换，得到霍夫矩阵；利用houghpeaks（）函数在霍夫矩阵中寻找峰值点；利用houghlines（）函数在之前两步结果的基础上得到原二值图像中的直线信息。[H,theta,rho]=hough(BW,param1,val1,param2,val2)1、霍夫变换(针对二值图像)56ppt课件三、Hough变换直线检测/8、图像分割步骤：[H,thet2、寻找峰值—houghpeakspeaks=houghpeaks(H,numpeaks,param1,val1,param2,val2)peaks是一个Q*2的矩阵，每行的两个元素分别为某一峰值点在Hough矩阵中的行、列索引，Q为找到的峰值点的数目。3、提取直线段—houghlineslines=houghlines(BW,theta,rho,peaks,param1,val1,param2,val2)57ppt课件2、寻找峰值—houghpeakspeaks=houghpe返回值lines的结构58ppt课件返回值lines的结构58ppt课件I1=imread('004.jpg');I=rgb2gray(I1);%旋转图像并寻找边缘rotI=imrotate(I,33,'crop');BW=edge(rotI,'canny');%执行Hough变换并显示Hough矩阵[H,T,R]=hough(BW);figure;imshow(H,[],'XData',T,'YData',R,'InitialMagnification','fit');xlabel('\theta'),ylabel('\rho');axison,axisnormal,holdon;%在Hough矩阵中寻找前5个大于Hough矩阵中最大值0.3倍的峰值P=houghpeaks(H,5,'threshold',ceil(0.3*max(H(:))));x=T(P(:,2));y=R(P(:,1));%由行、列索引转换成实际坐标plot(x,y,'s','color','white');%在Hough矩阵图像中标出峰值位置%找出并绘制直线lines=houghlines(BW,T,R,P,'FillGap',5,'MinLength',7);%合并距离小于5的线段，丢弃所有长度小于7的直线段figure,imshow(rotI),holdonmax_len=0;59ppt课件I1=imread('004.jpg');59ppt课件fork=1:length(lines)%依次标出各条直线段

xy=[lines(k).point1;lines(k).point2];plot(xy(:,1),xy(:,2),'x','LineWidth',2,'Color','green');

%绘制线段端点

plot(xy(1,1),xy(1,2),'x','LineWidth',2,'Color','yellow');plot(xy(2,1),xy(2,2),'x','LineWidth',2,'Color','red');%确定最长的线段

len=norm(lines(k).point1-lines(k).point2);if(len>max_len)max_len=len;xy_long=xy;endend%高亮显示最长线段plot(xy_long(:,1),xy_long(:,2),'LineWidth',2,'Color','cyan');60ppt课件fork=1:length(lines)%依次标出各条原图像Hough矩阵旋转并检测最长直线段61ppt课件原图像Hough矩阵旋转并检测最长直线段61ppt课件/9、特征提取特征提取的一般原则：选择在同类图像之间差异较小（较小的类内距），在不同类别的图像之间差异较大（较大的类间距）的图像特征。一、简单区域描绘子周长：区域边界上的像素数目。面积：区域中像素数目。致密性：（周长）^2/面积。区域的质心。灰度均值：区域中所有像素的平均值。灰度中值：区域中所有像素的排序中值。包含区域的最小矩形。最小或最大灰度级。大于或小于均值的像素数。欧拉数：区域中的对象数减去这些对象的孔洞数。62ppt课件/9、特征提取特征提取的一般原则：一、简单区域描绘子周长：区D=regionprops(L,properties);L为一个标记矩阵，通过连通1区域标注函数bwlabel得到。二、主成分分析（PrincipalComponentAnalysis，PCA）[COEFF,SCORE,latent]=pincomp(X);63ppt课件D=regionprops(L,properties);L为X为原始样本组成n*d的矩阵，其每一行是一个样本特征向量，每一列表示样本特征向量的一维。COEFF：主成分分量，也是样本协方差矩阵的本征向量。SCORE：主成分，X的低维表示。latent：一个包含着样本协方差矩阵本征值的向量。X=[1,2;3,3;3,5;4,5;5,6;6,5;8,7;9,3;7,6];X[COEFF,SCORE,latent]=princomp(X);%主成分分析COEFF%主成分分量SCORE%主成分latent%X样本协方差矩阵的本征值64ppt课件X为原始样本组成n*d的矩阵，其每一行是一个样本谢谢~v~65ppt课件谢谢~v~65ppt课件Matlab数字图像处理66ppt课件Matlab数字图像处理1ppt课件/1、图像的读取和显示/2、图像的点运算/4、空间域图像增强/5、频率域图像增强/6、彩色图像处理/7、形态学图像处理/8、图像分割/9、特征提取内容安排/3、图像的几何变换67ppt课件/1、图像的读取和显示/2、图像的点运算/4、空间域图像增强/1、图像的读取和显示一、图像的读取A=imread(FILENAME,FMT)FILENAME指定图像文件的完整路径和文件名。如果在work工作目录下只需提供文件名。FMT为图像文件的格式对应的标准扩展名。I_1=imread('D:\10.06.08nir\TTC10377.BMP');%读入图像二、图像的写入imwrite(A，FILENAME,FMT)FILENAME参数指定文件名。FMT为保存文件采用的格式。imwrite(I6,'nirdilatedisk2TTC10373.bmp');68ppt课件/1、图像的读取和显示一、图像的读取A=imread(FIL/1、图像的读取和显示三、图像的显示imshow(I,[lowhigh])I为要显示的图像矩阵。[lowhigh]为指定显示灰度图像的灰度范围。高于high的像素被显示成白色；低于low的像素被显示成黑色；介于High和low之间的像素被按比例拉伸后显示为各种等级的灰色。

figure；imshow(I6);title('TheMainPassPartofTTC10373');figure；%创建一个新的窗口figure；subplot（m,n,p）；imshow（I）；Subplot(m,n,p)含义为：打开一个有m行n列图像位置的窗口，并将焦点位于第p个位置上。69ppt课件/1、图像的读取和显示三、图像的显示imshow(I,[lo/1、图像的读取和显示四、图像的格式转换im2bw(I,LEVEL);rgb2gray;从RGB图创建灰度图，存储类型不变。im2uint8将图像转换成uint8类型阈值法从灰度图、RGB图创建二值图。LEVEL为指定的阈值；（0,1）。im2double将图像转换成double类型70ppt课件/1、图像的读取和显示四、图像的格式转换im2bw(I,LE/2、图像的点运算灰度直方图描述了一副图像的灰度级统计信息，主要应用于图像分割和图像灰度变换等处理过程中。从数学角度来说，图像直方图描述图像各个灰度级的统计特性，它是图像灰度值的函数，统计一幅图像中各个灰度级出现的次数或概率。归一化直方图可以直接反映不同灰度级出现的比率。横坐标为图像中各个像素点的灰度级别，纵坐标表示具有各个灰度级别的像素在图像中出现的次数或概率。imhist(I);%灰度直方图I=imread(‘red.bmp’);%读入图像

figure;%打开新窗口[M,N]=size(I);%计算图像大小[counts,x]=imhist(I,32);%计算有32个小区间的灰度直方图counts=counts/M/N;%计算归一化灰度直方图各区间的值stem(x,counts);%绘制归一化直方图一、图像直方图图像直方图归一化71ppt课件/2、图像的点运算灰度直方图描述了一副图像72ppt课件7ppt课件/2、图像的点运算二、灰度的线性变换Fa>1时，输出图像的对比度将增大；Fa<1时，输出图像对比度将减小。Fa=1且Fb非零时，所有像素的灰度值上移或下移，使整个图像更暗或更亮。Fa<0，暗区变亮，亮区变暗。73ppt课件/2、图像的点运算二、灰度的线性变换Fa>1时，输出图像的对74ppt课件9ppt课件75ppt课件10ppt课件/2、图像的点运算三、灰度的对数变换c为尺度比例常数，s为源灰度值，t为变换后的目标灰度值。k为常数。灰度的对数变换可以增强一幅图像中较暗部分的细节，可用来扩展被压缩的高值图像中的较暗像素。广泛应用于频谱图像的显示中。Warning：log函数会对输入图像矩阵s中的每个元素进行操作，但仅能处理double类型的矩阵。而从图像文件中得到的图像矩阵大多是uint8类型的，故需先进行im2double数据类型转换。76ppt课件/2、图像的点运算三、灰度的对数变换c为尺度比例常数，s为源I=imread('nir.bmp');%读入图像F=fft2(im2double(I));%FFTF=fftshift(F);%FFT频谱平移F=abs(F);T=log(F+1);%频谱对数变换figure;imshow(F,[]);title('未经变换的频谱');figure;imshow(T,[]);title('对数变换后');77ppt课件I=imread('nir.bmp');%读入图像12ppt/2、图像的点运算三、灰度的Gamma变换其中，x、y的取值范围为[0,1]。esp为补偿系数，r则为Gamma系数。Gamma变换是根据r的不同取值选择性的增强低灰度区域的对比度或者高灰度区域的对比度。78ppt课件/2、图像的点运算三、灰度的Gamma变换其中，x、y的取值J=imadjust(I,[low_inhigh_in],[low_outhigh_out],gamma)I=imread('nir.bmp');figure;imshow(imadjust(I,[],[],0.75));%gamma=0.5title('Gamma0.5');原NIR图像Gamma0.579ppt课件J=imadjust(I,[low_inhigh_in四、灰度阈值变换及二值化Gamma1.5T为指定阈值BW=im2bw(I,level);%level为人工设定阈值范围为[0,1]thresh=graythresh(I);%自动设定所需的最优化阈值80ppt课件四、灰度阈值变换及二值化Gamma1.5T为指定阈值BWOTSU算法：最大类间方差法自动单阈值分割。Kapur算法：一维直方图熵阈值算法niblack算法：局部阈值分割阈值的计算公式是T=m+k*v，其中m为以该像素点为中心的区域的平均灰度值，v是该区域的标准差，k是一个系数。kittlerMet:表示kittler最小分类错误(minimumerrorthresholding)全局二值化算法。

原图像直接阈值分割0.2581ppt课件OTSU算法：最大类间方差法自动单阈值分割。Kapur算法：自动阈值分割OTSU算法分割KittlerMet算法Niblack算法Kapur算法82ppt课件自动阈值分割OTSU算法分割KittlerMet算法Nibl/2、图像的点运算五、直方图均衡化[J,T]=histeq(I);%J为输出图像，T为变换矩阵图像易受光照、视角、方位、噪声等的影响。使得同一类图像的不同变形体之间的差距有时大于该类图像与另一类图像之间的差距，影响图像识别、分类。图像归一化就是将图像转换到唯一的标准形式以抵抗各种变换，从而消除同类图像不同变形体之间的外观差异。也称为图像灰度归一化。原图像及直方图83ppt课件/2、图像的点运算五、直方图均衡化[J,T]=histeq(图像变亮后灰度均衡化图像变暗后灰度均衡化84ppt课件图像变亮后灰度均衡化图像变暗后灰度均衡化19ppt课件/3、图像的几何变换一、图像平移正变换逆变换strel%用来创建形态学结构元素translate(SE,[yx])%原结构元素SE上y和x方向平移imdilate%形态学膨胀85ppt课件/3、图像的几何变换一、图像平移正变换逆变换20ppt课件/3、图像的几何变换I=imread('nir.bmp');se=translate(strel(1),[180190]);B=imdilate(I,se);figure;subplot(1,2,1),subimage(I);title('原图像');subplot(1,2,2),subimage(B);title('平移后图像');86ppt课件/3、图像的几何变换I=imread('nir.bmp');/3、图像的几何变换二、图像镜像B=imtransform(A,TFORM,method);TFORM=makeform(transformtype,Matrix);%空间变换结构参数transformtype指定了变换的类型，常见的’affine’为二维或多维仿射变换，包括平移、旋转、比例、拉伸和错切等。Matrix为相应的仿射变换矩阵。87ppt课件/3、图像的几何变换二、图像镜像参数transformtypA=imread('nir.bmp');[height,width,dim]=size(A);tform=maketform('affine',[-100;010;width01]);B=imtransform(A,tform,'nearest');tform2=maketform('affine',[100;0-10;0height1]);C=imtransform(A,tform2,'nearest');figure;imshow(A);figure;imshow(B);imwrite(B,'nir水平镜像.bmp');figure;imshow(C);imwrite(B,'nir垂直镜像.bmp');原图像水平镜像图像垂直镜像图像88ppt课件A=imread('nir.bmp');原图像水平镜像图像垂A=imread('nir.bmp');tform=maketform('affine',[010;100;001]);B=imtransform(A,tform,'nearest');figure;imshow(A);figure;imshow(B);imwrite(B,'nir转置后图像.bmp');/3、图像的几何变换三、图像转置89ppt课件A=imread('nir.bmp');/3、图像的几何变换/3、图像的几何变换四、图像中心旋转B=imrotate(A,angle,method,’crop’);angle为旋转角度，正值为逆时针旋转。可选参数method为imrotate函数指定插值方法。‘crop’选项会裁减旋转后增大的图像，保持和原图像同样大小。

A=imread('nir.bmp');B=imrotate(A,30,'nearest','crop');figure;imshow(B);imwrite(B,'逆时针中心旋转30度.bmp');逆时针30度90ppt课件/3、图像的几何变换四、图像中心旋转angle为旋转角度，正/4、空间域图像增强一、噪声添加h=imnoise(I,type,parameters);type为噪声类型，合法值如下：添加高斯白噪声添加椒盐噪声91ppt课件/4、空间域图像增强一、噪声添加h=imnoise(I,ty/4、空间域图像增强二、空间域滤波滤波过程就是在图像f(x,y)中逐点移动模板，使模板中心和点(x,y)重合，滤波器在每一点(x,y)的响应是根据模板的具体内容并通过预先定义的关系来计算的。92ppt课件/4、空间域图像增强二、空间域滤波滤波过程就是在图像f(x,B=imfilter(f,w,option1,option2,…);f为要进行滤波操作的图像。w为滤波操作使用的模板，为一个二维数组，可自己定义。option1……是可选项，包括：1、边界选项（’symmetric’、’replicate’、’circular’）2、尺寸选项（’same’、’full’）3、模式选项（’corr’、’conv’）原图像滤波后图像93ppt课件f为要进行滤波操作的图像。原图像滤波后图像28ppt课件三、滤波器设计/4、空间域图像增强h=fspecial(type,parameters)parameters为可选项，是和所选定的滤波器类型type相关的配置参数，如尺寸和标准差等。type为滤波器的类型。其合法值如下：94ppt课件三、滤波器设计/4、空间域图像增强h=fspecial(ty/4、空间域图像增强四、中值滤波h=medfilt2(I1,[m,n]);m和n为中值滤波处理的模板大小，默认3*3中值滤波本质上是一种统计排序滤波器。中值不同于均值，是指排序队列中位于中间位置的元素的值。中值滤波并非线性滤波器。对于某些类型的随机噪声具有非常理想的降噪能力。典型的应用就是消除椒盐噪声。95ppt课件/4、空间域图像增强四、中值滤波h=medfilt2(I1,/4、空间域图像增强五、图像锐化图像锐化主要用于增强图像的灰度跳变部分，主要通过运算导数（梯度）或有限差分来实现。主要方法有：Robert交叉梯度，Sobel梯度，拉普拉斯算子，高提升滤波，高斯-拉普拉斯变换。Robert交叉梯度w1对接近正45°边缘有较强响应，w2对接近负45°边缘有较强响应。96ppt课件/4、空间域图像增强五、图像锐化图像锐化主要用于增强图像的灰Sobel交叉梯度对水平边缘有较大响应对垂直边缘有较大响应拉普拉斯算子97ppt课件Sobel交叉梯度对水平边缘有较大响应对垂直边缘有较大响应拉I=imread('nir.bmp');I=double(I);%双精度化w1=[-10;01];w2=[0-1;10];G1=imfilter(I,w1,‘corr’,‘replicate’);%正45°梯度G2=imfilter(I,w2,‘corr’,‘replicate’);%负45°梯度G=abs(G1)+abs(G2);%计算Robert梯度figure;imshow(G,[]);figure;imshow(abs(G1),[]);figure;imshow(abs(G2),[]);Robert交叉梯度I=imread('nir.bmp');Id=double(I);%双精度化h_1=fspecial('log',5,0.5);%大小为5，sigma=0.5的LOG算子I_1=imfilter(Id,h_1,'corr','replicate');figure;imshow(uint8(abs(I_1)),[]);h_2=fspecial('log',5,2);%大小为5，sigma=2的LOG算子I_2=imfilter(Id,h_2,'corr','replicate');figure;imshow(uint8(abs(I_2)),[]);高斯-拉普拉斯锐化98ppt课件I=imread('nir.bmp');Robert交叉梯度/5、频率域图像增强用傅里叶变换表示的函数特征可以完全通过傅里叶反变换进行重建而不丢失任何信息吉布斯现象Gibbsphenomenon（又叫吉布斯效应）：将具有不连续点的周期函数（如矩形脉冲)进行傅立叶级数展开后，选取有限项进行合成。当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时，该峰起值趋于一个常数，大约等于总跳变值的9%。这种现象称为吉布斯现象。99ppt课件/5、频率域图像增强用傅里叶变换表示的函数特征可以完全通过傅幅度谱相位谱功率谱100ppt课件幅度谱相位谱功率谱35ppt课件一、傅里叶变换/5、频率域图像增强I=fft2(x);%快速傅里叶变换I=fft2(x,m,n);x为输入图像；m和n分别用于将x的第一和第二维规整到指定的长度。当m和n均为2的整数次幂时算法的执行速度要比m和n均为素数时快。I1=abs(I);%计算I的幅度谱I2=angle(I);%计算I的相位谱Y=fftshift(I);%频谱平移I=ifft2(x);%快速傅里叶逆变换I=ifft2(x,m,n);101ppt课件一、傅里叶变换/5、频率域图像增强I=fft2(x);%快速幅度谱相位谱平移后的相位谱102ppt课件幅度谱相位谱平移后的相位谱37ppt课件functionout=imidealflpf(I,freq)%imidealflpf函数构造高斯频域低通滤波器[M,N]=size(I);out=ones(M,N);fori=1:Mforj=1:Nif(sqrt(((i-M/2)^2+(j-N/2)^2))>freq)out(i,j)=0;endendend理想低通functionout=imgaussflpf(I,sigma)%imgaussflpf函数构造高斯频域低通滤波器[M,N]=size(I);out=ones(M,N);fori=1:Mforj=1:Nout(i,j)=exp(-((i-M/2)^2+(j-N/2)^2)/2/sigma^2);endend高斯低通103ppt课件functionout=imidealflpf(I,f/6、彩色图像处理色彩模型：RGB模型、CMY模型、CMYK模型、HIS模型、HSV模型、YUV模型、YIQ模型。RGB模型国际照明委员会（CIE）规定以蓝（435.8nm）、绿（546.1nm）和红（700nm）作为主原色。104ppt课件/6、彩色图像处理色彩模型：RGB模型、CMY模型、CMYKRGB_image=cat(3,PR,PG,PB);%将PR、PG、PB三个矩阵在第3个维度上进行级联，进行图像合成PR=RGB_image(:,:,1);%提取红色分量PG=RGB_image(:,:,2);%提取绿色分量PB=RGB_image(:,:,3);%提取蓝色分量Matlab中一幅RGB图像可表示为一个M*N*3的3维矩阵。其中每一个彩色像素都在特定空间位置的彩色图像中对应红、绿、蓝3个分量。CMY模型CMY模型是采用（Cyan、Magenta、Yellow）青、品红、黄色3中基本原色按一定比例合成颜色。由于色彩的显示是由光线被物体吸收掉一部分之后反射回来的剩余光线产生，故CMY模型又称为减色法混色模型。当光都被吸收时成为黑色，都被反射时为白色。105ppt课件RGB_image=cat(3,PR,PG,PB);%将PRCMY模型主要用于彩色打印机和复印机等。CMYK模型CMY模型在实际使用中，青、品红和黄色等比例混合后的黑色并不纯，为产生真正的黑色，专门加入第四种颜色—黑色。得到CMYK模型。用于四色打印。cmy=imcomplement(rgb);%rgb转成cmyrgb=imcomplement(cmy);cmy转成rgb106ppt课件CMY模型主要用于彩色打印机和CMYK模型CMY模型在实际使HSI模型HIS模型是从人的视觉系统出发，直接使用颜色三要素色调（Hue）、饱和度（Saturation）和亮度（Intensity）来描述颜色。亮度指人眼感觉光的明暗程度。光的能量越大，亮度越大。色调由物体反射光线中占优势的波长决定。反映颜色的本质。饱和度指颜色的深浅和浓淡程度，饱和度越高，颜色越深。HIS色彩空间比RGB彩色空间更符合人的视觉特性。亮度和色度具有可分离特性，使得图像处理和机器视觉中大量灰度处理算法都可在HIS彩色空间中方便使用。107ppt课件HSI模型HIS模型是从人的视觉系统出发，直接使用颜色三要素色调强度饱和度其中108ppt课件色调强度饱和度其中43ppt课件/7、形态学图像处理一、二值图像腐蚀I2=imerode(I,SE);SE=strel(shape,parameters);I为原始图像，可以是二值或者灰度图像。shape指定了结构元素的形状。parameters是和输入shape有关的参数。109ppt课件/7、形态学图像处理一、二值图像腐蚀I2=imerode(I/7、形态学图像处理二、二值图像膨胀I2=imdilate(I,SE);SE=strel(shape,parameters);I为原始图像，可以是二值或者灰度图像。shape指定了结构元素的形状。parameters是和输入shape有关的参数。原图像腐蚀后膨胀后110ppt课件/7、形态学图像处理二、二值图像膨胀I2=imdilate(/7、形态学图像处理三、其他二值图像运算SE=strel(shape,parameters);I2=imopen(I,SE);%开运算I3=imclose(I,SE);%闭运算Ihm=bwhitmiss(I,SE1,SE2);%击中击不中变换四、连通分量提取[Lnum]=bwlabel(Ibw,conn);Ibw为一幅输入二值图像。conn为可选参数，指明提取连通分量是4连通还是8连通。默认为8。L为连通分量标注图像。num为二值图像Ibw中连通分量个数。111ppt课件/7、形态学图像处理三、其他二值图像运算SE=strel(s/7、形态学图像处理五、bwmorph函数Iout=bwmorph(I,operation,n)112ppt课件/7、形态学图像处理五、bwmorph函数Iout=bwmo/7、形态学图像处理六、顶帽变换解决非均匀光照问题。顶帽变换定义为图像f与图像f的开运算之差。a=imread('red.bmp');I=rgb2gray(a);figure,surf(double(I(1:8:end,1:8:end))),zlim([0255])，colormapgray;%显示图像I的三维可视化效果bg=imopen(I,strel(‘disk’,15));%开运算Itophat=imsubtract(I,bg);%顶帽变换figure,imshow(Itophat);figure,surf(double(Itophat(1:8:end,1:8:end))),zlim([0255]);I2=imadjust(Itophat);figure,imshow(I2);113ppt课件/7、形态学图像处理六、顶帽变换解决非均匀光照问题。顶帽变换原图像三维可视化效果顶帽变换后的三维可视化效果顶帽变换后图像对比度拉伸后效果114ppt课件原图像三维可视化效果顶帽变换后的三维可视化效果顶帽变换后图像/8、图像分割一、图像分割概述图像分割一般采用的方法有边缘检测（edgedetection）、边界跟踪（edgetracing）、区域生长（regiongrowing）、区域分离和聚合等。

图像分割算法一般基于图像灰度值的不连续性或其相似性。

不连续性是基于图像灰度的不连续变化分割图像，如针对图像的边缘有边缘检测、边界跟踪等算法。

相似性是依据事先制定的准则将图像分割为相似的区域，如阈值分割、区域生长等。二、边缘检测图像的边缘点是指图像中周围像素灰度有阶跃变化或屋顶变化的那些像素点，即灰度值导数较大或极大的地方。

边缘检测可以大幅度的减少数据量，并且剔除不相关信息，保留图像重要的结构属性。

边缘检测基本步骤：平滑滤波、锐化滤波、边缘判定、边缘连接。115ppt课件/8、图像分割一、图像分割概述图像分割一般采边缘检测算法：基于一阶导数：Roberts算子、Sobel算子、Prewitt算子基于二阶导数：高斯-拉普拉斯边缘检测算子Canny边缘检测算法Matlab实现1）、基于梯度算子的边缘检测BW=edge(I,type,thresh,direction,’nothinning’)thresh是敏感度阈值参数，任何灰度值低于此阈值的边缘将不会被检测到。默认值为空矩阵[]，此时算法自动计算阈值。116ppt课件边缘检测算法：Matlab实现1）、基于梯度算子的边缘检测Bdirection指定了我们感兴趣的边缘方向，edge函数将只检测direction中指定方向的边缘，其合法值如下：可选参数’nothinning’，指定时可以通过跳过边缘细化算法来加快算法运行的速度。默认是’thinning’，即进行边缘细化。2）、基于高斯-拉普拉斯算子的边缘检测BW=edge(I,’log’,thresh,sigma)sigma指定生成高斯滤波器所使用的标准差。默认时，标准差为2。滤镜大小n*n，n的计算方法为：n=ceil(sigma*3)*2+1。117ppt课件direction指定了我们感兴趣的边缘方向，edge函数将3）、基于Canny算子的边缘检测BW=edge(I,’canny’,thresh,sigma)thresh是敏感度阈值参数，默认值为空矩阵[]。此处为一列向量，为算法指定阈值的上下限。第一个元素为阈值下限，第二个元素为阈值上限。如果只指定一个阈值元素，则默认此元素为阈值上限，其0.4倍的值作为阈值下限。如阈值参数没有指定，则算法自行确定敏感度阈值上下限。b1=imread('nir.bmp');h58=fspecial('gaussian',5,0.8);b=imfilter(b1,h58);bw1=edge(b,‘sobel’);%sobel算子bw2=edge(b,‘prewitt’);%prewitt算子bw3=edge(b,‘roberts’);%roberts算子bw4=edge(b,‘log’);%log算子bw5=edge(b,‘canny’);%canny算子figure;imshow(bw1);imwrite(bw1,'nirbwsobel.bmp');figure;imshow(bw2);imwrite(bw2,'nirbwprewitt.bmp');figure;imshow(bw3);imwrite(bw3,'nirbwroberts.bmp');figure;imshow(bw4);imwrite(bw4,'nirbwlog.bmp');figure;imshow(bw5);imwrite(bw5,'nirbwcanny.bmp');118ppt课件3）、基于Canny算子的边缘检测BW=ed