MATLAB与信号处理培训讲义.pptVIP

图象处理 5 基于MATLAB的图象处理 5.1 图象运算 5.2 图象变换 5.3 图象增强 5.4 图象分割 5.5 图象分析 5.1 图象运算 代数运算、几何运算、邻域操作和块操作。 1）代数运算： 指对两幅输入图象点对点的代数运算。 绝对差值函数 z=imabsdiff(x,y) I=imread( cameraman.tif ); J=uint8(filter2(fspecial(gaussian),I)); K=imabsdiff(I,J); subplot(1,2,1);imshow(I,[]); subplot(1,2,2);imshow(K,[]); 显示结果： 图象的叠加函数z=imadd(x,y) I=imread( cameraman.tif ); J=imread(rice.tif); K=imadd(I,J); imshow(K,[]); 增加亮度 I=imread(rice.tif); J=imadd(I,50); subplot(2,2,1);imshow(I,[]); subplot(2,2,2);imshow(J,[]); 图象的求补函数z=imcomplement(x) I = imread(text.tif); J = imcomplement(I); subplot(1,2,1);imshow(I,[]); subplot(1,2,2);imshow(J,[]) 灰度图象补运算 I = imread(bonemarr.tif); J = imcomplement(I); subplot(1,2,1);imshow(I,[]); subplot(1,2,2);imshow(J,[]) 其它计算： imdivide immultiply imsubtract imlincomb 2）几何操作 改变图象大小：y=imresize(x,m,method) 图象旋转：y=imrotate(x,angle,method) I = imread(circuit.tif); J = imrotate(I,-15,bilinear); imshow(I); figure; imshow(J) 3）邻域和块操作 实现一个滑动邻域操作步骤： （1）选择一个单独的像素。 （2）确定该像素的滑动邻域。 （3）对邻域中的像素值应用一个函数求值。 （4）将结果作为对应像素的值。 （5）对输入图象的每一个像素都重复上面四个步骤。 滑动邻域操作函数： y=nlfilter(x,[m n],fun，p1,p2,…) I=imread(tire.tif); f=inline(max(x(:))); J=nlfilter(I,[3 3],f); imshow(I);figure;imshow(J); 块操作函数： y=blkproc(x,[m n],fun，p1,p2,…) I=imread(cameraman.tif); f=@dct; J=blkproc(I,[8 8],f); imagesc(J); I = imread(alumgrns.tif) fun = inline(std2(x)*ones(size(x))); I2 = blkproc(I,[8 8],fun); imshow(I) figure, imshow(I2,[]) 5.2 图象变换 1）FFT 例：绘制一个二值图象矩阵，将其傅立叶变换可视化 f=zeros(40,40); f(5:24,13:17)=1; subplot(1,2,1);imshow(f); F=fft2(f); F2=log(abs(F)); subplot(1,2,2); imshow(F2,[-1 5]); 图象定位： w=imread(text.tif); a=w(59:71,81:91); imshow(w);figure;imshow(a); C=ifft2(fft2(w).*fft2(rot90(a,2),256,256)); figure;imshow(C,[]); thresh=45; figure;imshow(Cthresh); 离散余弦变换： 两种方法：dct2 dctmtx 例：对图象进行DCT变换和反变换 RGB=imread(autumn.tif); I=rgb2gray(RGB); subplot(1,2,1);imshow(I); J=dct2(I); subplot(1,2,2);imshow(J); 将上例中DCT变换后的矩阵中小于100的设