采用imnoise函数对图像添加不同分布、不同强度的随机噪声,比较通过空域平滑、中值滤波处理前后图像的频谱

数字图像处理实验报告

一、

实验要求

采用

imnoise

函数对图像添加不同分布、不同强度的随机噪声，比较通过空域平滑、中值滤

波处理前后图像的频谱。

二、

实验代码

clc 

clear 

warning off 

I=imread('pout.tif'); 

if isrgb(I)%

判断是否为彩色图像

I=rgb2gray(I);%

彩色图像转黑白，若原图像为黑白的，则不要这一步

end 

%

加噪声

I3=imnoise(I,'salt & pepper',0.03); 

I12=imnoise(I,'gaussian'); 

I21=imnoise(I,'speckle'); 

%

原图像求频谱

A=fft2(I);        %

傅里叶变换

A=fftshift(A);      %

使图像对称

I1=log(1+abs(A));      %

图像幅度谱，加

log

便于显示

I2=real(120*I1./abs(I1)); 

%

椒盐频谱

C=fft2(I3);        %

傅里叶变换

I4=log(1+abs(C));      %

图像幅度谱，加

log

便于显示

C=fftshift(C);      %

使图像对称

I5=real(120*I4./abs(I4)); 

%

椒盐中值

n=3;  

a=ones(n,n); 

p=size(I3);   

x1=double(I3);x2=x1; 

for i=1:p(1)-n+1 

    for j=1:p(2)-n+1 

        c=x1(i:i+(n-1),j:j+(n-1));  

        e=c(1,:); 

        for u=2:n 

            e=[e,c(u,:)];    

        end 

        mm=median(e); 

        x2(i+(n-1)/2,j+(n-1)/2)=mm;