资源描述：
《数字图像处理》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、%题目：实现一个直方图均衡化的程序,首先将一幅256级灰度图像,通过灰度变换，%将灰度范围压缩到64到192之间，再采用直方图均衡化程序处理。closeallclcclearI=imread('e:matlab_file3.jpg');[m,n]=size(I);M=zeros(1,256);fork=0:255M(k+1)=length(find(I==k))/(m*n);endS1=zeros(1,256);sum=0;fori=1:256if(i>=64&&i<=192)S1(i)=255*M(i)/128-255*64/128/(m*n);endendS2
2、=zeros(1,256);fori=1:256forj=1:iS2(i)=S1(j)+S2(i);endendS3=zeros(1,256);fori=1:256S3(i)=floor((S2(i)*255)+0.5);endd=zeros(size(I));fori=1:mforj=1:nd(i,j)=S3(I(i,j)+1);endendd=uint8(d);Q=zeros(1,256);x=0;fork=0:255Q(k+1)=length(find(d==k))/(m*n);endimwrite(d,'E:matlab_file1.jpg');figur
3、e(1);subplot(321);imshow(I);title('原图像');subplot(322);imshow(d);title('均衡化后的图像');subplot(323);bar(0:255,M);title('原图像直方图');xlabel('灰度值');ylabel('每级灰度出现概率');subplot(324);bar(0:255,S1);title('压缩后图像直方图');xlabel('灰度值');ylabel('每级灰度出现概率');figure(2);bar(0:255,Q);title('均衡化后图像直方图');xlabel('灰度值
4、');ylabel('每级灰度出现概率');仿真结果：%题目：实现频率域的UnSharpMasking滤波。其中低通滤波使用频率域的高斯低通滤波器滤波器参数可调。clear;closeall;clc;I=imread('E:matlab_file3.jpg');figuresubplot(321)imshow(I)title('原图像')B=fft2(I);s=fftshift(B);[M,N]=size(s);n=2;d0=20;n1=floor(M/2);%对M/2进行取整n2=floor(N/2);%对N/2进行取整fori=1:Mforj=1:Nd=sqr
5、t((i-n1)^2+(j-n2)^2);h=1*exp(-1/2*(d^2/d0^2));a(i,j)=h*s(i,j);endenda=ifftshift(a);a=uint8(real(ifft2(a)));subplot(322)imshow(a)title('GLPF滤波(d0=20)')subplot(323)b=I-a;imshow(b)title('反锐化')c0=I+0.5*b;subplot(324)imshow(c0)title('反锐化图像A=0.5')c1=I+1*b;subplot(325)imshow(c1)title('反锐化图像A=1
6、')c2=I+1.5*b;subplot(326)imshow(c2)title('反锐化图像A=1.5')仿真结果：%题目：对由于运动模糊并加入噪声而退化的图像Lena,%(需要从我这里拷贝该退化图像,图像尺寸256x256,8位灰度,无格式文件),%给出采用最小均方误差滤波(Wiener滤波)的复原图像.%其中,运动模糊中的参数为:a=0.1,b=0.1,T=1.%加入的噪声为均值为0,方差为50的加性高斯噪声.%滤波中的参数K采用交互选择的方式确定.%同学也可自己选择一幅灰度图像,采用上述参数生成退化图像,然后再给出采用最小均方误差滤波的复原图像.clear;c
7、lccloseall;I=imread('E:matlab_file3.jpg');figure;subplot(2,2,1);imshow(I);title('原图像');[m,n]=size(I);F=fftshift(fft2(I));k=0.025;foru=1:mforv=1:nH(u,v)=exp((-k)*(((u-m/2)^2+(v-n/2)^2)^(5/6)));endendG=F.*H;I0=real(ifft2(fftshift(G)));subplot(2,2,2);imshow(uint8(I0));title('模糊退
