基于卡尔滤波的图像除噪处理汇编.docx

基于卡尔曼滤波的图像除噪处理摘要：图像在采集和和传输的过程中可能要受到外界的干扰，为了保证图像信息的准确性，除噪就显得尤为重要。本文介绍的就是利用Matlab软件使用卡尔曼滤波进行图像的除噪处理。关键词：图像除噪；Matlab；卡尔曼滤波一、引言目前大多数数字图像系统中，输入图像都是采用先冻结再扫描方式将多维图像变成一维电信号，再对其进行处理、存储、传输等加工变换。最后往往还要再组成多维图像信号，而图像噪声也将同样受到这样的分解和合成。在这些过程中电气系统和外界影响将使得图像噪声的精确分析变得十分复杂。因此，我们需要对图像进行除噪，得到更准确的图像信息。一、Matlab2.1 简介Matlab是美国MathWorks公司推出的数学软件，是一种面向工程和科学计算的交互式计算软件。主要用于算法开发、数据可视化、数据分析、以及数值计算，但是由于Matlab拥有许多的附加工具箱，因此它也可以应用到其他领域，譬如：图像处理、信号检测、音频分析以及金融建模和分析等等。2 Matlab在图像处理中的应用图像文件格式的读写和显示。MATLAB提供了图像文件的读入函数imread()，用来读取，如：bmp、jpg、hdf等格式的图像文件；图像写出函数为imwrite()，以及图像显示函数image()、imshow()等。图像处理的基本运算。MATLAB提供了图像的和、差等线性运算，以及卷积、滤波等非线性运算。图像变换。MATLAB提供了一维和二维离散傅里叶变换（DFT）、快速傅里叶变换（FFT）、离散余弦变换（DCT）以及连续小波变换（CWT）等等。图像的分析和增强。为了能对图像进行处理，MALTAB针对图像的统计计算提供了校正、中值滤波、直方图均衡、对比度调整以及自适应滤波等。2.3.均方误差（MSE)在相同测量条件下进行的测量称为等精度测量，例如在同样的条件下，用同一个游标卡尺测量铜棒的直径若干次，这就是等精度测量。对于等精度测量来说，还有一种表示误差的方法，就是标准误差。标准误差定义为各测量值的平方值的平方根，故又称为均方误差。设n个测量值的误差为，则这组测量值的标准误差等于：数理统计中均方误差是指参数估计值与参数真值之差平方的期望值，记为MSE。2.4 峰值信噪比PSNR（Peak Signal to Noise Ratio）表示的是到达噪声比率的顶点信号，PSNR一般是用于最大值信号和噪声之间的工程项目。为了衡量经过处理后的图像品质，我们通常会参考PSNR值来认定某个处理程序是否有效，或者说是处理的效果是怎么样的。Peak就是8bits表示法的最大值255，MES指的是Mean Square Error(均方误差，各值相差的n次方的平均值)。PSNR的计算公式如下：I是指原始图像第n个pixel值，P是指处理后的图像第n个pixel值。PSNR的单位为db，所以PSNR的值越大，就代表失真越少。卡尔曼滤波3.1卡尔曼滤波的特点卡尔曼滤波是用状态空间法描述系统的，有状态方程和量测方程组成。卡尔曼滤波是用一个状态的估计值和最近一个观测数据来估计状态变量的当前值，并以状态变量的估计值的形式给出。卡尔曼滤波具有以下的特点;算法是递推的；卡尔曼滤波使用于非平稳过程；卡尔曼滤波采取的误差准则仍为估计误差的均方值最小。3.2 状态方程和量测方程假设某系统k时刻的状态变量为，状态方程和量测方程可以表示为：其中，k表示时间,指的是第k步迭代时相应信号的取值；输入信号是一白噪声，输出信号的观测噪声也是一个白噪声；A表示状态变量之间的增益矩阵，表示第k步迭代时，增益矩阵A的取值；C表示状态变量与输出信号之间的增益矩阵，第k步迭代时取值用表示。3.3卡尔曼滤波的递推算法当不考虑观测噪声和输入信号时，状态方程和量测方程为：由于不考虑观测噪声的影响，输出信号的估计值与实际值是有误差的，用表示：为了提高状态估计的质量，用输出信号的估计误差来校正状态变量：其中，为增益矩阵（实质上是一个加权矩阵）。经过校正后的状态变量的估计误差以及其均方差分别用和表示，卡尔曼滤波要求状态变量的估计误差的均方差为最小，而且卡尔曼滤波的关键就是计算出加权矩阵的最佳值。把未经校正的状态变量的估计误差的均方差用表示：MATLAB仿真与分析4.1 原理基于MATLAB平台设计卡尔曼滤波器对一张.jpg图片进行除噪处理，先导入一张图片，然后依次给图片加高斯白噪声和椒盐噪声，再利用卡尔曼滤波器除噪。除噪之后利用Matlab输出除噪前后的均方误差和峰值信噪比。源程序：%1.显示原图中分别加入椒盐噪声和高斯白噪声clear;close all;mine = imread(C:\Users\yuerui\Desktop\hzw.jpg); %读取图像show_size = size(min