数字图像处理图像平滑和锐化变换处理Word格式文档下载.docx

1、）;J=imadjust（I,0.2,0.4,）;subplot（2,2,1）,imshow（I）;subplot（2,2,2）,imshow（J）;subplot（2,2,3）,imhist（I）;subplot（2,2,4）,imhist（J）;直方图均衡程序如下： J=histeq（I）;Subplot（2,2,1）;Imshow（I）;Title（原图像Subplot（2,2,2）;Imshow（J）;直方图均衡化后的图像） ;Subplot（2,2,3） ;Imhist（I,64）;原图像直方图Subplot（2,2,4）;Imhist（J,64） ; Title（均衡变换后的直方图

2、伽马校正程序如下：A=imread（x=0:255;a=80,b=1.8,c=0.009;B=b.（c.*（double（A）-a）-1;y=b.（c.*（x-a）-1;subplot（3,2,1）;imshow（A）;subplot（3,2,2）;imhist（A）;subplot（3,2,3）;imshow（B）;subplot（3,2,4）;imhist（B）;subplot（3,2,6）;plot（x,y）;log变换程序如下：Image=imread（subplot（1,2,1）;imshow（Image）;Image=log（1+double（Image）;subplot（1,2,

3、2）;imshow（Image,）;onion.pngB=rgb2gray（A）;figure,subplot（3,2,1）;title（原始图象H=imnoise（B,gaussianimshow（H）;高斯噪声Q=imnoise（B,salt & pepperimshow（Q）;椒盐噪声M=fspecial（average,3*3）;E=imfilter（Q,M）;imshow（E）;3*3平均模板N=fspecial（,5*5）;K=imfilter（Q,N）;subplot（3,2,5）;imshow（K）;5*5平均模板Z=fspecial（,7*7）;J=imfilter（Q,Z）

4、;imshow（J）;7*7平均模板频域低通滤波处理噪声的程序如下：f1,f2=freqspace（25,meshgridHd=zeros（25,25）;d=sqrt（f1.2+f2.2）0.5;%0.5为截止半径大小Hd（d）=1;h=fsamp2（Hd）;figure（1）,freqz2（h,64,64）;RGB=imread（I=rgb2gray（RGB）;I1=imnoise（I,I2=imnoise（I,0.02）;I3=imnoise（I,speckleJ=imfilter（I,h,replicateJ1=imfilter（I1,h,J2=imfilter（I2,h,J3=imfi

5、lter（I3,h,figure（2）,subplot（221）,imshow（J）;原图像滤波后subplot（222）,imshow（J1）;高斯污染图像滤波后subplot（223）,imshow（J2）;椒盐污染图像污染后subplot（224）,imshow（J3）;乘法污染图像滤波后4、 空域锐化：A=RGB2gray（I）;figure,subplot（2,3,1）;原图hs=fspecial（sobelS=imfilter（A,hs）;hp=fspecial（prewittP=imfilter（A,hs）;A=double（A）;%双精度型H=0,1,0;1,-4,1;0,1,

6、0;%拉普拉斯算子J=conv2（A,H,sameK=A-J;subplot（2,3,2）,imshow（K）;拉普拉斯锐化图像B=edge（A,roberts,0.1）;subplot（2,3,3）,imshow（B）;罗伯特锐化图像subplot（2,3,4）,imshow（S）;sobel算子锐化图像subplot（2,3,5）,imshow（P）;prewitt算子锐化图像figure,subplot（131）,imshow（B）;B=double（B）;B=fftshift（fft2（B）;E=B;M,N=size（B）;k1=round（M/2）;k2=round（N/2）;D=5

7、0;for i=1:1:M for j=1:N juli=sqrt（i-k1）2+（j-k2）2）; K（i,j）=1-exp（-1/2）*juli2*（1/D2）;E（i,j）=K（i,j）*E（i,j）; endendE=ifft2（ifftshift（E）;E=uint8（real（E）;subplot（132）,imshow（E）;高斯高通滤波器四、实验结果及分析1、灰度变换：灰度拉伸运行结果如下：直方图均衡运行结果如下：伽马校正运行结果如下：log变换运行结果如下：2、空域平滑：3、频域平滑：用频率采样法构建的二维滤波器：figure（1）图像经过二维滤波器滤波后图像，figure（2）5、 频域锐化：