中南民族大学数字图像处理程序及图像.doc

仅供参考 定有错误 木有绝对正确 实验1 1—1 close all; clear; f=zeros(40,40);% f(10:30,10:30)=1; subplot(1,3,1);% imshow(f);% F=fft2(f) subplot(1,3,2);% imshow(F);% D=log(1+abs(F));% subplot(1,3,3); imshow(D); 1——2 clc load imdemos saturn2 subplot(121) imshow(saturn2) i=fftshift(fft2(saturn2)); subplot(122) imshow(log(abs(i)),[]),colormap(jet(64)),colorbar; 1棗3 lw=imread(text.png); a=w(33:45,88:98); figure; imshow(w); figure; imshow(a); C=real(ifft2(fft2(w).*fft2(rot90(a,2),256,256))); figure; imshow(C,[]); max(C(:)); thresh=60; figure; imshow(Cthresh); 1棗4 eRGB=imread(autumn.tif); figure(1),imshow(RGB); I=rgb2gray(RGB); figure(2),imshow(I); J=dct2(I); figure(3),imshow(log(abs(J)),[]),colormap(jet(64)),colorbar; J(abs(J)10)=0; K=idct2(J)/255; figure(4),imshow(K); 1棗5棗2 eacr=0.125; sig=imread(cameraman.tif); sig=double(sig)/255; figure(1),imshow(sig); [m_sig,n_sig]=size(sig); sizi=8; snum=64; t=hadamard(sizi); hdcoe=blkproc(sig,[sizi sizi],P1*x*P2,t,t); coe=im2col(hdcoe,[sizi sizi],distinct); coe_temp=coe; [Y Ind]=sort(coe); [m,n]=size(coe); snum=snum-snum*cr; for i=1:n coe_temp(Ind(1:snum),i)=0; end re_hdcoe=col2im(coe_temp,[sizi sizi],[m_sig n_sig],distinct); re_sig=blkproc(re_hdcoe,[sizi sizi],P1*x*P2,t,t); re_sig=double(re_sig)/64; figure(2); imshow(re_sig); error=sig.^2-re_sig.^2; MSE=sum(error(:))/prod(size(re_sig)); 2棗1 1i=imread(pout.tif); imshow(i); figure; subplot(1,2,1); imhist(i); j=histeq(i,64); subplot(1,2,2) imhist(j) 2棗2 ef=imread(autumn.tif); gl=imadjust(f,[0 1],[1,0]); subplot(121); imshow(f); subplot(122); imshow(gl); 2棗3 af=imread(autumn.tif); gl=imadjust(f,[0.5 0.75],[1,0]); subplot(121); imshow(f); subplot(122); imshow(gl); 2棗4 ef=imread(autumn.tif); gl=imadjust(f,[ ],[ ],2); subplot(121); imshow(f); subplot(122); imshow(gl); test 1 I=imread(pout.tif); [row,col] = size(I); imshow(I); K=histeq(I,64) subplot(121) imhist(K); title(直方图均衡化) J = zeros(row,col); for i =1:row for j = 1: col if (I(i,j) 10 I(i,j) 0) J(i,j) = 0.5*I(i,j);