[matlab数字图像处理5]对一副图像加噪声（高斯，椒盐等），进行平滑，锐化，观察其图像变化

一、 问题描述

对一副图像加噪声（高斯，椒盐等），进行平滑，锐化，观察其图像变化

二、 求解思路

首先利用matlab自带的imread函数读取选择的图片，然后将图片灰度化并显示原图像，加噪声使用的是matlab自带的imnoise函数，可直接选择高斯噪声或者椒盐噪声。对于图像的平滑与锐化，需要先创建一个掩膜函数，输入的参量为需要处理的图片的灰度化后的矩阵和掩膜矩阵，在函数处理中，需要对于边缘一圈的量，需要补充一圈“0“（默认掩膜为33或者22的矩阵），然后利用循环，将每个需要处理的区域与掩膜矩阵做卷积，加和后，将所得的结果赋予该区域的中心。对于平滑处理，使用的是全为1/9的33掩膜矩阵，对于锐化处理，使用的是[-1,0;0,1]的22矩阵，计算完成后，将输出的结果显示。

三、 程序代码与实验结果

%对一副图像加噪声（高斯，椒盐等），进行平滑，锐化，观察其图像变化；

clear;
clc;
close all
x=imread('2.jpeg');
x=rgb2gray(x);
%显示原图像
figure(1)
imshow(x)

%显示加高斯噪声和椒盐噪声后的图像

figure(2)
subplot(3,2,1)
x1=imnoise(x,'gaussian');
imshow(x1)
title('加高斯噪声后');
x2=imnoise(x,'salt & pepper');
subplot(3,2,2)
imshow(x2)
title('加椒盐噪声后');

平滑掩膜

%对加高斯噪声和椒盐噪声后的函数分别进行平滑，并作出对比图
t1=(1/9)*ones(3,3) ;  
x3=uint8(fun(x1,t1));
subplot(3,2,3)
imshow(x3)
title('加高斯噪声后平滑')
x4=uint8(fun(x2,t1));
subplot(3,2,4)
imshow(x4)
title('加椒盐噪声后平滑')

锐化掩膜

%对加高斯噪声和椒盐噪声后的函数分别进行锐化，并作出加噪声，平滑，锐化后的对比图
t2=[-1,0;0,1];
x5=uint8(fun(x1,t2));
subplot(3,2,5)
imshow(x5)
title('加高斯噪声后锐化')
x6=uint8(fun(x2,t2));
subplot(3,2,6)
imshow(x6)
title('加椒盐噪声后锐化')

figure(3)
%加高斯噪声后的锐化大图
imshow(x5)

figure(4)
%加椒盐噪声后的锐化大图
imshow(x6)

figure(5)
%将原图与原图锐化后的图叠加
imshow(x+uint8(fun(x,t2)))

%x为原图像，t为掩膜矩阵
function y=fun(x,t)
     row=size(x,1);
     col=size(x,2);
     new_x=zeros(row+2,col+2);
     for i=2:row+1
         for j=2:col+1
             new_x(i,j)=x(i-1,j-1);
         end
     end
     y=zeros(row,col);
     row_t=size(t,1);
     col_t=size(t,2);
     for i=1:row
         for j=1:col
             y_k=new_x(i:i+row_t-1,j:j+col_t-1);
             y(i,j)=sum(sum(y_k.*t));
         end
     end
end

四、 实验心得

在进行图像的掩膜处理时，应注意掩膜矩阵的大小，并且应对最周围一圈的情况做针对处理