矩阵的奇异值分解在数字图像处理的应用要点.doc

矩阵的奇异值分解在数字图像处理的应用浅析 学 院： 专 业： 姓 名： 学 号： 目 录 一、绪 论 - 1 - 二、数字图像处理简介 - 2 - 三、矩阵的奇异值分解原理 - 4 - 3.1 矩阵的奇异值 - 4 - 3.2 矩阵的奇异值分解（SVD） - 4 - 四、奇异值分解的图像性质 - 5 - 五、图像的奇异值分解压缩方法 - 7 - 5.1 奇异值分解压缩原理分析 - 7 - 5.2 奇异值分解压缩应用过程 - 8 - 六、小 结 - 9 - 一、绪 论 目前，随着科学技术的高速发展，现实生活中有大量的信息用数字进行存储、处理和传送。而传输带宽、速度和存储器容量等往往有限制，JPEG 等许多标准都已经得到了广泛的应用。奇异值分解 Singular Value Decomposition ，SVD [2]，并取得了成功，被视为一种有效的图像压缩方法。本文在奇异值分解的基础上进行图像压缩。 二、数字图像处理简介 首先，简单介绍一下数字图像处理。人们对数字图像都应该很熟悉。我们在计算机上看到的图像，数码相机拍到的图像，雷达图像，人体MRI图像等等。数字图像处理用计算机对图像进行分析，以达到所需结果技术，这里和是空间坐标，在空间坐标上的幅值称为该点图像的强度或灰度。对于数字图像而言，和幅值都是有限的、离散的。这样一幅图像就可以用一个二维函数来表示。模拟图像不利于计算机处理，所以我们常常将模拟图像转换为数字图像。模拟图像转化为数字图像的方式是：取样和量化。我们将坐标值离散化称为取样，将幅度值离散化称之为量化。经过取样和量化的图像是一幅数字图像。数字图像的质量很大程度上取决于取样和量化的取样数和灰度级。取样和量化的结果是一个实际的矩阵。这个矩阵可以表示为 更一般的矩阵表达方式为： 图像压缩是数据压缩技术在数字图像上的应用，它的目的是减少图像数据中的冗余信息从而用更加高效的格式存储和传输数据。图像数据之所以能被压缩，就是因为数据中存在着冗余。图像数据的冗余主要表现为：图像中相邻像素间的相关性引起的空间冗余；图像序列中不同帧之间存在相关性引起的时间冗余；不同彩色平面或频谱带的相关性引起的频谱冗余。 图像压缩可以是有损数据压缩也可以是无损数据压缩。无损图像压缩方法主要有行程长度编码、熵编码法如LZW；有损压缩方法主要有变换编码，如离散余弦变换（DCT）或者小波变换这样的傅立叶相关变换，然后进行量化和用熵编码法压缩和分形压缩（fractal compression）。 图像矩阵的奇异值（Singular Value）及其特征空间反映了图像中的不同成分和特征。奇异值分解 Singular Value Decomposition ，SVD 是一种基于特征向量的矩阵变换方法，在信号处理、模式识别、数字水印技术等方面都得到了应用讨论奇异值分解在图像压缩上的应用。 三、矩阵的奇异值分解原理 3.1 矩阵的奇异值 设，，是的特征值，是的特征值，它们都是实数。且设 则特征值与之间的关系为：，。 设， 的正特征值，的正特征值，称， 是的正奇异值，简称奇异值。若是正规矩阵，则的奇异值是的非零特征向量的模长。 若，是的个正奇异值，则存在阶酉矩阵和阶酉矩阵，满足： 其中，，为奇异对角阵。满足是对角阵，满足是对角阵。的第列为的对应于奇异值对应的左奇异向量，的第列为的对应于奇异值对应的右奇异向量。它们的每一列均为单位向量，且各列之间相互正交。 若，是的个正奇异值，则总有次酉矩阵，满足：，其中。 奇异值分解是一种基于特征向量的矩阵变换方法。奇异值分解是现代数值的最基本和最重要的工具之一。 任意一个矩阵的奇异值是唯一的，它刻画了矩阵数据的分布特征。直观上，可以这样理解矩阵的奇异值分解：将矩阵看成是一个线性变换，它将维空间的点映射到维空间。经过奇异值分解后，这种变换被分割成3个部分，分别为、和，其中和都是标准正交矩阵，它们对应的线性变换就相当于对维和维坐标系中坐标轴的旋转变换。 若为数字图像，则可视为二维时频信息，可将的奇异值分解公式写为： 其中，和分别是和的列矢量，是的非零奇异值。故上式表示的数字图像可以看成是个秩为1的子图叠加的结果，而奇异值为权系数。所以也表示时频信息，对应的和可分别视为频率矢量和时间矢量，因此数字图像中的时频信息就被分解到一系列由和构成的视频平面中。 由矩阵范数理论， 奇异值能与向量2-范数和矩阵Frobenious-范数 F-范数 相联系。 若以F-范数的平方表示图像的能量，则由矩阵奇异值分解的定义知： 。 也就是说，数字图像经奇异值分解后，其纹理和几何信息都集中在、之中，而中的奇异值则代表图像的能量信息。 性质1：矩阵的奇异值代表图像的能量信息，因而具有稳定性。 设，，是矩阵的一个扰动矩阵。和的非零奇异值分别记为：和。