MATLAB数字图像处理实验一

小咪咪 2023-05-21 15:16 14阅读 0赞

**一、实验要求**  
*1、直方图(读入你自己的图片并显示； 用默认、条形图、杆状图及曲线四种方式显示直方图，并将它们绘制在同一幅figure中进行比较)  
2、九种默认的滤波器并显示结果；  
3、对你自己的图添加噪声水平为0.15的椒盐噪声；用大小为3X3，5X5，7X7的统计滤波（取第8顺位）对其进行去噪处理，并将它们绘制在同一幅figure中进行比较；  
4、对你的图像进行傅里叶变换并直接显示；*

**二、实验代码**  
(1)

%%
    %
    I=imread('hei.png');
    I=rgb2gray(I);
    %im2uint8(mat2gray(log(1+double(I))));
    
    h=imhist(I);
    subplot(2,2,1),imhist(I),title('默认图');
    subplot(2,2,2),plot(h),title('曲线图');
    subplot(2,2,3),stem(h),title('杆状图');
    subplot(2,2,4),bar(h),title('条形图');

(2)

%%
    %九种滤波器
    I=imread('hei.png');
    imshow(I),title('原图'),figure;
    
    L1=imfilter(I,fspecial('average',3));
    subplot(3,3,1),imshow(L1),title('average');
    
    L2=imfilter(I,fspecial('disk',5));
    subplot(3,3,2),imshow(L2),title('disk');
    
    L3=imfilter(I,fspecial('gaussian'));
    subplot(3,3,3),imshow(L3),title('gaussian');
    
    L4=imfilter(I,fspecial('laplacian'));
    subplot(3,3,4),imshow(L4),title('laplacian');
    
    L5=imfilter(I,fspecial('log'));
    subplot(3,3,5),imshow(L5),title('log');
    
    L6=imfilter(I,fspecial('prewitt'));
    subplot(3,3,6),imshow(L6),title('prewitt');
    
    L7=imfilter(I,fspecial('motion'));
    subplot(3,3,7),imshow(L7),title('motion');
    
    L8=imfilter(I,fspecial('sobel'));
    subplot(3,3,8),imshow(L8),title('sobel');
    
    L9=imfilter(I,fspecial('unsharp'));
    subplot(3,3,9),imshow(L9),title('unsharp');

(3)

%%
    %
    I=imread('hei.tif');
    figure,imshow(I),title('原图');
    %%
    %
    I2=imnoise(I,'salt & pepper',0.05);%%注意中间空格，不然会错报
    figure,imshow(I2),title('椒盐噪声')
    %%
    %
    I=imread('hei.tif');
    figure;
    subplot(2,2,1),set(gca,'position'),imshow(I),title('原图');
    K = imfilter(I2,fspecial('average',3));
    subplot(2,2,2),set(gca,'position'),imshow(K),title('3*3');
    L = imfilter(I2,fspecial('average',5));
    subplot(2,2,3),set(gca,'position'),imshow(L),title('5*5');
    M = imfilter(I2,fspecial('average',7));
    subplot(2,2,4),set(gca,'position'),imshow(M),title('7*7');

(4)

%%傅里叶变换
    %%A基本操作
    I = imread('hei.tif');
    imshow(I);
    title('源图像');
    J1=rgb2gray(I);
    J2= fft2(J1);
    figure, imshow(J2);
    title('傅里叶变换');
    %%频移
    JSh = fftshift(J2);
    figure, imshow(JSh);
    title('傅里叶变换频移');
    %%直接傅立叶反变换
    Ji=ifft2(J2);
    figure, imshow(Ji/256);
    title('直接傅立叶变换');
    %%幅度
    JA=abs(J2);
    iJA=ifft2(JA);
    figure, imshow(JA/256);
    title('幅度博立叶反变换');
    %%相位
    JP=angle(J2);
    iJP=ifft2(JP);
    figure, imshow(abs(iJP)* 100);
    title('相位傅立叶反变换');