数字图象与处理实验题解析 .pptVIP

数字图象处理;A1 编程实现图象同时对比度;A1 结果;A2 编程实现马赫带效应;B1 编程实现图象读入,计算并显示图象的复利叶功率谱; 进一步实现 巴特沃斯高通滤波 巴特沃斯低通滤波 ;[N1,N2]=size(C);%Butterworth高通滤波器 n=2; d0=5; n1=fix(N1/2); n2=fix(N2/2); for i=1:N1 for j=1:N2 d=sqrt((i-n1)^2+(j-n2)^2); if d==0 h=0; else h=1/(1+(d0/d)^(2*n)); end D(i,j)=h*C(i,j); end end E=ifft2(ifftshift(D)); F=uint8(real(E)); subplot(234),imshow(F);title(Butterworth高通滤波结果;G=imnoise(B,salt pepper,0.02);%添加椒盐噪声 subplot(235),imshow(G);title(添加椒盐噪声) H=fftshift(fft2(double(G))); [N1,N2]=size(H);%Butterworth低通滤波器 n=2; d0=50; n1=fix(N1/2); n2=fix(N2/2); for i=1:N1 for j=1:N2 d=sqrt((i-n1)^2+(j-n2)^2); h=1/(1+0.414*(d/d0)^(2*n)); I(i,j)=h*H(i,j); end end J=ifft2(ifftshift(I)); K=uint8(real(J)); subplot(236),imshow(K);title(Butterworth低通滤波结果);B1 结果;B2 编程实现图象的同态滤波;for i=1:M for j=1:N d=sqrt((i-n1)^2+(j-n2)^2); h=(Rh-Rl)*(1-exp(-c*(d.^2/d0.^2)))+Rl; g(i,j)=h*g(i,j); end end g=ifftshift(g); g=uint8(real(ifft2(g))); subplot(122),imshow(g);title(同态滤波效果);;C1 编程实现：3x3 平均值滤波及中值滤波 并比较两者的不同;C1 结果;C1 均值滤波与中值滤波不同： ;C2 编程实现RGB与HSI的转换，在此基础上 实现对彩色图像的饱和度、色调的改变 ;figure,subplot(1,3,1),imshow(H,[]),title(色度图H); set(gcf,outerposition,get(0,screensize)); set(gcf,NumberTitle,off,Name,HSI通道图像); subplot(1,3,2),imshow(S),title(饱和度图S); subplot(1,3,3),imshow(I),title(强度图I); ? %hsi2rgb HV=hsi(:,:,1)*2*pi; SV=hsi(:,:,2); IV=hsi(:,:,3); R=zeros(size(HV)); G=zeros(size(HV)); B=zeros(size(HV)); ? %RG Sector id=find((0=HV) (HV2*pi/3)); B(id)=IV(id).*(1-SV(id)); R(id)=IV(id).*(1+SV(id).*cos(HV(id))./cos(pi/3-HV(id))); G(id)=3*IV(id)-(R(id)+B(id));;%BG Sector id=find((2*pi/3=HV) (HV4*pi/3)); R(id)=IV(id).*(1-SV(id)); G(id)=IV(id).*(1+SV(id).*cos(HV(id)-2*pi/3)./cos(pi-HV(id))); B(id)=3*IV(id)-(R(id)+G(id)); %BR Sector id=find((4*pi/3=HV) (HV2*pi)); G(id)=IV(id).*(1-SV(id)); B(id)=IV(id).*(1+SV(id).*cos(HV(id)-4*pi/3)./cos(5*pi/3-HV(id))); R(id)=3*IV