图像预处理例子.docVIP

图像傅里叶变换 例1. 利用傅里叶变换图像，观察分析变换结果 A=imread(E:\张儒良图像处理课件\car1.bmp); B=imread(E:\张儒良图像处理课件\girl.bmp); A1=fft2(double(A)); B1=fft2(double(B)); subplot(1,4,1);imshow(A); subplot(1,4,2);imshow(A1); subplot(1,4,3);imshow(B); subplot(1,4,4);imshow(B1); 例2. 利用傅里叶变换图像，然后用逆傅里叶变换观察分析变换结果 A=imread(E:\张儒良图像处理课件\car1.bmp); B=imread(E:\张儒良图像处理课件\girl.bmp); A1=fft2(double(A)); B1=fft2(double(B)); A2=abs(ifft2(A1)); B2=abs(ifft2(B1)); subplot(2,3,1);image(A); subplot(2,3,2);imshow(A1); subplot(2,3,3);image(A2); subplot(2,3,4);image(B); subplot(2,3,5);imshow(B1); subplot(2,3,6);image(B2); 例3. 利用傅里叶变换图像，然后用逆傅里叶变换观察压缩结果 A=imread(E:\张儒良图像处理课件\apple1.bmp); m=rgb2gray(A); A1=fft2(double(m)); s=size(A1); P1=zeros(s); P1(1:s(1),1:s(2)/1.1)=A1(1:s(1),1:s(2)/1.1); E0=(ifft2(A1)); E1=abs(ifft2(P1)); subplot(2,2,1);image(m); subplot(2,2,2);imshow(A1); subplot(2,2,3);image(E0); subplot(2,2,4);image(E1); 二、离散余弦变换 例4. 利用离散余弦变换图像，然后用逆离散余弦变换 A=imread(E:\张儒良图像处理课件\car1.bmp); B=imread(E:\张儒良图像处理课件\girl.bmp); C=imread(E:\张儒良图像处理课件\car.bmp); A1=dct2(A); B1=dct2(B); C1=dct2(C); subplot(2,3,1);imshow(A); subplot(2,3,2);imshow(A1); subplot(2,3,3);imshow(B); subplot(2,3,4);imshow(B1); subplot(2,3,5);imshow(C); subplot(2,3,6);imshow(C1) 例5. 利用离散余弦变换图像，然后用逆离散余弦变换观察压缩结果 D=imread(E:\张儒良图像处理课件\car1.bmp); s=size(D) D1=dct2(D); D2=idct2(D1); P=zeros(s); P(1:100,1:100)=D1(1:100,1:100); E=idct2(P); subplot(2,2,1); imshow(D); subplot(2,2,2); imshow(D1); subplot(2,2,3); image(D2); subplot(2,2,4); image(E); 例6. 编程实现离散余弦变换图像（实在太慢） A=imread(E:\张儒良图像处理课件\car1.bmp); A=double(A); s=size(A) M=s(1); N=s(2); ap0=sqrt(1/M); ap=sqrt(2/M); aq0=sqrt(1/N); aq=sqrt(2/N); for p=1:M-1 for q=1:N-1 k=0; for m=0:M-1 for n=0:N-1 k=A(m+1,n+1)*cos(pi*(2*m+1)*p/(2*M))*cos(pi*(2*n+1)*q/(2*N))+k; end end B(p+1,q+1)=ap*aq*k; end end B 三、灰度变换 例7：灰度变换例程(p64) clc I=imread(E:\张儒良图像处理课件\car.bmp); J1=imadjust(I,[0.3,0.7],[0,1],1); J2=imadjust(I,[0.3,0.7],[0,1],0.5); J3=imadjust(I,[0.3,0.7],[0,1],1.5); J4=imadjust(I,[0.3,0.7],[1,0],1.5); subp