Matlab图像处理简易入门教程PPT学习课件.ppt

Matlab数字图像处理/1、图像的读取和显示/1、图像的读取和显示/1、图像的读取和显示/1、图像的读取和显示/2、图像的点运算7/2、图像的点运算910/2、图像的点运算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('对数变换后');/2、图像的点运算J=imadjust(I,[low_inhigh_in],[low_outhigh_out],gamma)四、灰度阈值变换及二值化OTSU算法：最大类间方差法自动单阈值分割。自动阈值分割/2、图像的点运算图像变亮后灰度均衡化/3、图像的几何变换/3、图像的几何变换/3、图像的几何变换A=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');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、图像的几何变换/4、空间域图像增强/4、空间域图像增强B=imfilter(f,w,option1,option2,…);三、滤波器设计/4、空间域图像增强/4、空间域图像增强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),[]);/5、频率域图像增强幅度谱一、傅里叶变换幅度谱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/6、彩色图像处理RGB_image=cat(3,PR,PG,PB);%将PR、PG、PB三个矩阵在第3个维度上进行级联，进行图像合成PR=RGB_image(:,:,1);%提取红色分量PG=RGB_image(:,:,2);%提取绿色分量PB=RGB_image(:,:,3);%提取蓝色分量CMY模型主要用于彩色打印机和复印机等。HSI模型色调/7、形态学图像处理/7、形态学图像处理/7、形态学图像处理/7、形态学图像处理/7、形态学图像处理原图像三维可视化效果/8、图像分割边缘检测算法：基于一阶导数：Roberts算子、Sobel算子、Prewitt算子基于二阶导数：高斯-拉普拉斯边缘检测算子Canny边缘检测算法direction指定了我们感兴趣的边缘方向，edge函数将只检测direction中指定方向的边缘，其合法值如下：3）、基于Canny算子的边缘检测Sobel算子分析三、Hough变换直线检测2、寻找峰值—houghpeaksParam合法值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,'YDa