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

第一页，共二十二页，2022年，8月28日 A1 编程实现图象同时对比度 d=zeros(32,32); %背景灰度值为0，图像尺寸为32x32 d(13:20,13:20)=2; %中心方块为8x8，灰度值为2 figure; %建立图形窗口 subplot(221); imshow(d,[0,5]); %显示图像d d=zeros(32,32)+1; %背景灰度值为1，图像尺寸为32x32 d(13:20,13:20)=2; subplot(222); imshow(d,[0,5]);? d=zeros(32,32)+3; %背景灰度值为3，图像尺寸为32x32 d(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]); 第二页，共二十二页，2022年，8月28日 A1 结果 第三页，共二十二页，2022年，8月28日 A2 编程实现马赫带效应 clear d=zeros(6,40); %背景灰度值为0，图像尺寸为6x40 d(1:6,1:2)=0; d(1:6,3:4)=1;d(1:6,5:6)=2 d(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 第四页，共二十二页，2022年，8月28日 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(傅立叶功率谱) 第五页，共二十二页，2022年，8月28日 [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高通滤波结果 第六页，共二十二页，2022年，8月28日 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));