MATLAB与信号处理课件.ppt

MATLAB与数字信号处理;5 基于MATLAB的图象处理;5.1 图象运算;显示结果：;图象的叠加函数z=imadd(x,y) I=imread( cameraman.tif ); J=imread(rice.tif); K=imadd(I,J); imshow(K,[]); ;增加亮度;图象的求补函数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,[]) ;16;5.2 图象变换;18;图象定位： 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);;20;离散余弦变换： 两种方法：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的设为0，重构图象。 J(abs(j)100)=0; K=idct2(J); figure; imshow(K,[0 255]);;Dctdemo 重构图象和原始图象比较;Radon 变换: [R xp]=radon(I,theta) 例：计算正方形图象在0和45方向上的变换 I=zeros(100,100); I(25:75,25:75)=1; imshow(I); [R,xp]=radon(I,[0 45]); subplot(1,2,1);plot(xp,R(:,1)); subplot(1,2,2);plot(xp,R(:,2));;例：构造一个简单图象的投影并重建图象 P=phantom(256); subplot(2,2,1); imshow(P); theta1=0:10:170; theta2=0:5:175; theta3=0:2:178; [R1,xp]=radon(P,theta1); [R2,xp]=radon(P,theta2); [R3,xp]=radon(P,theta3);;I1=iradon(R1, theta1); I2=iradon(R2, theta2); I3=iradon(R3, the