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

A1编程实现图象同时对比度d=zeros(32,32);%背景灰度值为0，图像尺寸为32x32d(13:20,13:20)=2;%中心方块为8x8，灰度值为2figure;%建立图形窗口subplot(221);imshow(d,[0,5]);%显示图像dd=zeros(32,32)+1;%背景灰度值为1，图像尺寸为32x32d(13:20,13:20)=2;subplot(222);imshow(d,[0,5]);?d=zeros(32,32)+3;%背景灰度值为3，图像尺寸为32x32d(13:20,13:20)=2;subplot(223);imshow(d,[0,5]);d=zeros(32,32)+5;%背景灰度值为5，图像尺寸为32x32，d(13:20,13:20)=2;subplot(224);imshow(d,[0,5]);第2页,共22页，星期六，2024年，5月A1结果第3页,共22页，星期六，2024年，5月A2编程实现马赫带效应cleard=zeros(6,40);%背景灰度值为0，图像尺寸为6x40d(1:6,1:2)=0;d(1:6,3:4)=1;d(1:6,5:6)=2d(1:6,7:8)=3;d(1:6,9:10)=4;d(1:6,11:12)=5;d(1:6,13:14)=6;d(1:6,15:16)=7;d(1:6,17:18)=8;d(1:6,19:20)=9;d(1:6,21:22)=10;d(1:6,23:24)=11;d(1:6,25:26)=12;d(1:6,27:28)=13;d(1:6,29:30)=14;d(1:6,31:32)=15;d(1:6,33:34)=16;d(1:6,35:36)=17;d(1:6,37:38)=18;d(1:6,39:40)=19;figure;%显示图形窗口imshow(d,[0,20],notruesize);%显示图形d第4页,共22页，星期六，2024年，5月B1编程实现图象读入,计算并显示图象的复利叶功率谱;进一步实现巴特沃斯高通滤波巴特沃斯低通滤波A=imread(‘shilong.jpg’);%图像读入B=rgb2gray(A);%RGB图像转换灰度图像subplot(231),imshow(B);title(灰度图像)C=fftshift(fft2(double(B)));%显示图像傅立叶频谱subplot(232),imshow(log(abs(C)),[]);title(傅立叶频谱)%显示图像傅立叶功率谱subplot(233),imshow(log(((abs(C))^2)),[]);title(傅立叶功率谱)第5页,共22页，星期六，2024年，5月[N1,N2]=size(C);%Butterworth高通滤波器n=2;d0=5;n1=fix(N1/2);n2=fix(N2/2);fori=1:N1forj=1:N2d=sqrt((i-n1)^2+(j-n2)^2);ifd==0h=0;elseh=1/(1+(d0/d)^(2*n));endD(i,j)=h*C(i,j);endendE=ifft2(ifftshift(D));F=uint8(real(E));subplot(234),imshow(F);title(Butterworth高通滤波结果第6页,共22页，星期六，2024年，5月G=imnoise(B,saltpepper,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);fori=1:N1forj=1:N2d=sqrt((i-n1)^2+(j-n2)^2);h=1/(1+0.414*(d/d0)^(2*n));I(i,j)=h*H(i,j);