（全新）大连理工数字图像处理作业.docVIP

数字图像作业 陈柯锦 电子1102 201181442 4-4 代码： a=[7,6,5,4,-4,-5,-6,-7]; b=[2,2,-5,-5,6,6,-7,-7]; subplot(211); stem(a); subplot(212); stem(b); figure; N=8; a1=[a zeros(1,N-length(a))]; b1=[b zeros(1,N-length(a))]; y1=conv(a1,b1); z1=[zeros(1,N) y1(1:(N-1))]; z2=[y1((N+1):(2*N-1)) zeros(1,N)]; z=z1(1:(2*N-1))+z2(1:(2*N-1))+y1(1:(2*N-1)); y2=z(1:N); stem((0:N-1),y2,.); C=fft(y2); C1=fft(a).*fft(b); all(C(:)==C1(:)); 圆周卷积函数： function y=yuanzhou(x1,x2,N) x1=[x1,zeros(1,N-length(x1))]; x2=[x2,zeros(1,N-length(x2))]; m=0:1:N-1; H=zeros(N,N); for n=1:1:N H(n,:)=cirshift(x2,n-1,N); end y=x1*H; end 结果： 4-9 代码： t1=imread(Lena.tif); t=rgb2gray(t1); imshow(t),title(原图像); [a,b]=size(t); [x,y]=meshgrid(-a/2:a/2-1,-b/2:b/2-1); z=sqrt(x.^2+y.^2); c=(z15); tf=fftshift(fft2(t)); figure; imshow(mat2gray(log(1+abs(tf)))); title(原图频域表示); tf1=tf.*c; figure; ty2=ifft2(tf1); imshow(abs(ty2),[]),title(低通滤波之后图像); c1=(z15); tf2=tf.*c1; ty3=ifft2(tf2); figure; imshow(abs(ty3),[]),title(高通滤波之后图像); g1=mat2gray(fspecial(gaussian,512,10)); tf3=tf.*g1; ty4=ifft2(tf3); figure; imshow(abs(ty4),[]),title(高斯低通滤波器滤波之后图像); bl=lbutter(c,15,1); tfbl=tf.*bl; ty5=ifft2(tfbl); figure; imshow(abs(ty5),[]),title(巴特沃兹低通滤波之后图像); 巴特沃兹滤波器函数设计： function out=lbutter(im,d,n) height=size(im,1); width=size(im,2); [x,y]=meshgrid(-floor(width/2):floor((width-1)/2),-floor(height/2): ... floor((height-1)/2)); out=1./(1+(sqrt(2)-1)*((x.^2+y.^2)/d^2).^n); 结果 5-6 代码： f=imread(Lena.tif); fg=rgb2gray(f); tn=imnoise(fg,salt pepper,0.05); imshow(tn,[]),title(原图像); a1=fspecial(average); tna1=filter2(a1,tn); figure; imshow(tna1,[]),title(均值滤波之后图像); tnz=medfilt2(tn,[3,3]); figure; imshow(tnz),title(中值滤波之后图像); tnout=outlier(tn,0.2); figure; imshow(tnout,[]),title(异常值判断滤波之后D=0.2); tnord=ordfilt2(tn,3,[0 1 0; 1 1 1;0 1 0]); figure; imshow(tnord,[]),title(伪中值滤波之后图像); 伪中值滤波函数： function res=outlier(im,d) f=[0.125 0.125 0.125; 0.125 0 0.125; 0.125 0.125 0.125]; imd=im2double(im); imf=filter2(f,imd); r=abs(imd-imf)-d0; re