资源描述:
《数字图象处理实验题解析》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、数字图象处理张晓与地图学与地理信息系统20120711A1编程实现图象同时对比度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;%背景
2、灰度值为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]);A1结果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:1
3、0)=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;%显示图形窗口ims
4、how(d,[0,20],'notruesize');%显示图形dB1编程实现图象读入,计算并显示图象的复利叶功率谱;进一步实现巴特沃斯高通滤波巴特沃斯低通滤波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('傅立叶频谱
5、')%显示图像傅立叶功率谱subplot(233),imshow(log(((abs(C))^2)),[]);title('傅立叶功率谱')[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(ifftshi
6、ft(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);fori=1:N1for
7、j=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);endendJ=ifft2(ifftshift(I));K=uint8(real(J));subplot(236),imshow(K);title('Butterworth低通滤波结果')B1结果B2编程实现图象的同态滤波%滤波参数通过键盘输入(Rl=0.5、Rh=2.0)I=imread('shilong.jpg');%装入真彩图J=rgb2gray(I);%
8、将真彩图转换为灰度图subplot(121),imshow(J);title('灰度图像');J=double(J);f=fft2(J);%采用傅里叶变换g=fftshift(f);%数据矩阵平衡[M,N]=size(f);%进行同态滤波处理d0=10;Rl=input('请输入Rl=');%低频分量Rh=input('请输入Rh=');%高频分量c=4;n1=floor(M/2);n2=floor