一、案例说明
三幅经典的图像(lena、peppers、barbara)分别加上了高斯噪声、乘性噪声、椒盐噪声。高斯白噪声:均值为0,方差为0.05。椒盐噪声:噪声密度0.10。 乘性噪声:这里为均匀分布的均值为0,方差为0.10的噪声。 (1)、使用均值滤波、中值滤波进行测试,得出去噪效果并分析; (2)、使用高斯低通滤波器、巴特沃斯滤波器进行测试,得出去噪效果并分析; (3)、使用PCA进行测试,得出去噪效果并分析; (4)、使用小波变换进行测试,可以采用‘db’小波、‘sym’小波等,也可以自己去选择一些小波函数,得出去噪效果并分析; (5)、自己选择一些方法、得出去噪效果并分析。
二、图像去噪及滤波简介
1 图像去噪
1.1 图像噪声定义
噪声是干扰图像视觉效果的重要因素,图像去噪是指减少图像中噪声的过程。噪声分类有三种:加性噪声,乘性噪声和量化噪声。我们用f(x,y)表示图像,g(x,y)表示图像信号,n(x,y)表示噪声。
图像去噪是指减少数字图像中噪声的过程。现实中的数字图像在数字化和传输过程中常受到成像设备与外部环境噪声干扰等影响,称为含噪图像或噪声图像。去噪是图像处理研究中的一个重点内容。在图像的获取、传输、发送、接收、复制、输出等过程中,往往都会产生噪声,其中的椒盐噪声是比较常见的一种噪声,它属于加性噪声。
1.2 图像噪声来源
(1)图像获取过程中
图像传感器CCD和CMOS采集图像过程中受传感器材料属性、工作环境、电子元器件和电路结构等影响,会引入各种噪声。
(2)图像信号传输过程中
传输介质和记录设备等的不完善,数字图像在其传输记录过程中往往会受到多种噪声的污染。
1.3 噪声分类
噪声按照不同的分类标准可以有不同的分类形式:
基于产生原因:内部噪声,外部噪声。
基于噪声与信号的关系:
加性噪声:加性噪声和图像信号强度是不相关的,这类带有噪声的图像g可看成为理想无噪声图像f与噪声n之和:
g = f + n;
乘性嗓声:乘性噪声和图像信号是相关的,往往随图像信号的变化而变化,载送每一个象素信息的载体的变化而产生的噪声受信息本身调制。在某些情况下,如信号变化很小,噪声也不大。为了分析处理方便,常常将乘性噪声近似认为是加性噪声,而且总是假定信号和噪声是互相统计独立。
g = f + f*n
按照基于统计后的概率密度函数:
是比较重要的,主要因为引入数学模型这就有助于运用数学手段去除噪声。在不同场景下噪声的施加方式都不同,由于在外界的某种条件下,噪声下图像-原图像(没有噪声时)的概率密度函数(统计结果)服从某种分布函数,那么就把它归类为相应的噪声。下面将具体说明基于统计后的概率密度函数的噪声分类及其消除方式。
1.4 图像去噪算法的分类
(1)空间域滤波
空域滤波是在原图像上直接进行数据运算,对像素的灰度值进行处理。常见的空间域图像去噪算法有邻域平均法、中值滤波、低通滤波等。
(2)变换域滤波
图像变换域去噪方法是对图像进行某种变换,将图像从空间域转换到变换域,再对变换域中的变换系数进行处理,再进行反变换将图像从变换域转换到空间域来达到去除图像嗓声的目的。将图像从空间域转换到变换域的变换方法很多,如傅立叶变换、沃尔什-哈达玛变换、余弦变换、K-L变换以及小波变换等。而傅立叶变换和小波变换则是常见的用于图像去噪的变换方法。
(3)偏微分方程
偏微分方程是近年来兴起的一种图像处理方法,主要针对低层图像处理并取得了很好的效果。偏微分方程具有各向异性的特点,应用在图像去噪中,可以在去除噪声的同时,很好的保持边缘。偏微分方程的应用主要可以分为两类:一种是基本的迭代格式,通过随时间变化的更新,使得图像向所要得到的效果逐渐逼近,这种算法的代表为Perona和Malik的方程,以及对其改进后的后续工作。该方法在确定扩散系数时有很大的选择空间,在前向扩散的同时具有后向扩散的功能,所以,具有平滑图像和将边缘尖锐化的能力。偏微分方程在低噪声密度的图像处理中取得了较好的效果,但是在处理高噪声密度图像时去噪效果不好,而且处理时间明显高出许多。
(4)变分法
另一种利用数学进行图像去噪方法是基于变分法的思想,确定图像的能量函数,通过对能量函数的最小化工作,使得图像达到平滑状态,现在得到广泛应用的全变分TV模型就是这一类。这类方法的关键是找到合适的能量方程,保证演化的稳定性,获得理想的结果。
形态学噪声滤除器将开与闭结合可用来滤除噪声,首先对有噪声图像进行开运算,可选择结构要素矩阵比噪声尺寸大,因而开运算的结果是将背景噪声去除;再对前一步得到的图像进行闭运算,将图像上的噪声去掉。据此可知,此方法适用的图像类型是图像中的对象尺寸都比较大,且没有微小细节,对这类图像除噪效果会较好。
2 中值滤波
(1)概念:
(2)原理解释:
3 均值滤波
均值滤波是指任意一点的像素值,都是周围 N \times M 个像素值的均值。例如下图中,红色点的像素值是其周围蓝色背景区域像素值之和除25,25=5\times5 是蓝色区域的大小。
均值滤波详细的计算方法如下图所示:
其中5\times5的矩阵称为核,针对原始图像内的像素点,采用核进行处理,得到结果图像,如下图所示:
提取 1/25 可以将核转换为如下形式:
4 Lee滤波
三、部分源代码
function varargout = main(varargin)% MAIN MATLAB code for main.fig%MAIN, by itself, creates a new MAIN or raises the existing%singleton*.%%H = MAIN returns the handle to a new MAIN or the handle to%the existing singleton*.%%MAIN('CALLBACK',hObject,eventData,handles,...) calls the local%function named CALLBACK in MAIN.M with the given input arguments.%%MAIN('Property','Value',...) creates a new MAIN or raises the%existing singleton*. Starting from the left, property value pairs are%applied to the GUI before main_OpeningFcn gets called. An%unrecognized property name or invalid value makes property application%stop. All inputs are passed to main_OpeningFcn via varargin.%%*See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one%instance to run (singleton)".%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help main% Last Modified by GUIDE v2.5 17-Jan- 10:47:23% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name', mfilename, ...'gui_Singleton', gui_Singleton, ...'gui_OpeningFcn', @main_OpeningFcn, ...'gui_OutputFcn', @main_OutputFcn, ...'gui_LayoutFcn', [] , ...'gui_Callback', []);if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});endif nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});elsegui_mainfcn(gui_State, varargin{:});end% End initialization code - DO NOT EDIT% --- Executes just before main is made visible.function main_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to main (see VARARGIN)% Choose default command line output for mainhandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes main wait for user response (see UIRESUME)% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line.function varargout = main_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT);% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structurevarargout{1} = handles.output;% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global I I1[filename,pathname] = uigetfile({'*.png';'*.jpg';'*.bmp'},'选择图片');str = [pathname,filename];I1 = imread(str);if (strcmpi(filename,'lena_gaussian.png')==1||strcmpi(filename,'lena_salt-pepper.png')==1||strcmpi(filename,'lena_speckle.png')==1)I=imread('lena.png');elseif (strcmpi(filename,'barbara_gaussian.png')==1||strcmpi(filename,'barbara_salt-pepper.png')==1||strcmpi(filename,'barbara_speckle.png')==1)I=imread('barbara.png');elseif (strcmpi(filename,'peppers256_gaussian.png')==1||strcmpi(filename,'peppers256_salt-pepper.png')==1||strcmpi(filename,'peppers256_speckle.png')==1)I=imread('peppers256.png');elsemsgbox('请输入规定的噪声图像','确认','error');endaxes(handles.axes1),imshow(I),title('原图');PSNR1=psnr(I1,I);t1 = sprintf('%s%f%s','加噪后图像,PSNR=',PSNR1, 'dB');axes(handles.axes2),imshow(I1),title(t1);% --- Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global I I1I2=average_filtering(I1); I3=median_filtering(I1); I4=gaussian_lowpass_filtering(I1); I5=butterworth_lowpass_filtering(I1); I6=PCA_filtering(I1); I7=wavelet_transform_filtering(I1); I8=wiener_filtering(I1); PSNR2=psnr(I2,I);t2 = sprintf('%s%f%s','均值滤波后图像,PSNR=',PSNR2, 'dB');axes(handles.axes3),imshow(I2),title(t2);PSNR3=psnr(I3,I);t3 = sprintf('%s%f%s','中值滤波后图像,PSNR=',PSNR3, 'dB');axes(handles.axes4),imshow(I3),title(t3);PSNR4=psnr(I4,I);t4 = sprintf('%s%f%s','高斯低通滤波后图像,PSNR=',PSNR4, 'dB');axes(handles.axes5),imshow(I4),title(t4);PSNR5=psnr(I5,I);t5 = sprintf('%s%f%s','巴特沃斯滤波后图像,PSNR=',PSNR5, 'dB');axes(handles.axes6),imshow(I5),title(t5);PSNR6=psnr(I6,I);t6 = sprintf('%s%f%s','PCA滤波后图像,PSNR=',PSNR6, 'dB');axes(handles.axes7),imshow(I6),title(t6);PSNR7=psnr(I7,I);t7 = sprintf('%s%f%s','小波变换后图像,PSNR=',PSNR7, 'dB');axes(handles.axes8),imshow(I7),title(t7);PSNR8=psnr(I8,I);t8 = sprintf('%s%f%s','维纳滤波后图像,PSNR=',PSNR8, 'dB');axes(handles.axes9),imshow(I8),title(t8);
四、运行结果
五、matlab版本及参考文献
1 matlab版本
a
2 参考文献
[1] 蔡利梅.MATLAB图像处理——理论、算法与实例分析[M].清华大学出版社,.
[2]杨丹,赵海滨,龙哲.MATLAB图像处理实例详解[M].清华大学出版社,.
[3]周品.MATLAB图像处理与图形用户界面设计[M].清华大学出版社,.
[4]刘成龙.精通MATLAB图像处理[M].清华大学出版社,.
[5]基于matlab的传统算法图像去噪的实现原理
【图像去噪】基于matlab GUI均值+中值+高通低通+巴特沃斯+PCA+小波+维纳滤波图像去噪【含Matlab源码 1705期】