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

1、矩阵的奇异值分解在数字图像处理的应用浅析 学 院：··· 专 业：··姓 名：··学 号：··2011年11月6日 目 录一、绪 论- 1 -二、数字图像处理简介- 2 -三、矩阵的奇异值分解原理- 4 -3.1 矩阵的奇异值- 4 -3.2 矩阵的奇异值分解（SVD）- 4 -四、奇异值分解的图像性质- 5 -五、图像的奇异值分解压缩方法- 7 -5.1 奇异值分解压缩原理分析- 7 -5.2 奇异值分解压缩应用过程- 8 -六、小 结- 9 -一、绪 论目前，随着科学技术的高速发展，现实生活中有大量的

2、信息用数字进行存储、处理和传送。而传输带宽、速度和存储器容量等往往有限制，因此数据压缩就显得十分必要。数据压缩技术已经是多媒体发展的关键和核心技术。图像文件的容量一般都比较大，所以它的存储、处理和传送会受到较大限制，图像压缩就显得极其重要。当前对图像压缩的算法有很多，特点各异，类似JPEG 等许多标准都已经得到了广泛的应用。奇异值分解(Singular Value Decomposition ，SVD) 是一种基于特征向量的矩阵变换方法，在信号处理、模式识别、数字水印技术等方面都得到了应用。由于图像具有矩阵结构，有文献提出将奇异值分解应用于图像压缩2，并取得了成功，被视为一种有效的图像压缩方法

3、。本文在奇异值分解的基础上进行图像压缩。二、数字图像处理简介首先，简单介绍一下数字图像处理。人们对数字图像都应该很熟悉。我们在计算机上看到的图像，数码相机拍到的图像，雷达图像，人体MRI图像等等。数字图像处理是指用计算机对图像进行分析处理，以达到所需结果的技术。图像处理的内容十分广泛，具体而言，可以分为：图像获取、图像增强、图像复原、图像压缩、图像分割等。这些内容都是基于矩阵的处理得到的。下面举例介绍几个重要的应用。图像获取是图像处理的第一步。图像获取有很多方法，最常用的方法就是用传感器如数字摄像机、扫描仪等设备得到。数字图像处理的定义：我们可以将一幅图像定义为一个二维函数，这里和是空间坐标，

4、在空间坐标上的幅值称为该点图像的强度或灰度。对于数字图像而言，和幅值都是有限的、离散的。这样一幅图像就可以用一个二维函数来表示。模拟图像不利于计算机处理，所以我们常常将模拟图像转换为数字图像。模拟图像转化为数字图像的方式是：取样和量化。我们将坐标值离散化称为取样，将幅度值离散化称之为量化。经过取样和量化的图像是一幅数字图像。数字图像的质量很大程度上取决于取样和量化的取样数和灰度级。取样和量化的结果是一个实际的矩阵。这个矩阵可以表示为更一般的矩阵表达方式为：图像压缩是数据压缩技术在数字图像上的应用，它的目的是减少图像数据中的冗余信息从而用更加高效的格式存储和传输数据。图像数据之所以能被压缩，就是

5、因为数据中存在着冗余。图像数据的冗余主要表现为：图像中相邻像素间的相关性引起的空间冗余；图像序列中不同帧之间存在相关性引起的时间冗余；不同彩色平面或频谱带的相关性引起的频谱冗余。图像压缩可以是有损数据压缩也可以是无损数据压缩。无损图像压缩方法主要有行程长度编码、熵编码法如LZW；有损压缩方法主要有变换编码，如离散余弦变换（DCT）或者小波变换这样的傅立叶相关变换，然后进行量化和用熵编码法压缩和分形压缩（fractal compression）。图像矩阵的奇异值（Singular Value）及其特征空间反映了图像中的不同成分和特征。奇异值分解(Singular Value Decomposit

6、ion ，SVD) 是一种基于特征向量的矩阵变换方法，在信号处理、模式识别、数字水印技术等方面都得到了应用。我们主要讨论奇异值分解在图像压缩上的应用。三、矩阵的奇异值分解原理3.1 矩阵的奇异值设，是的特征值，是的特征值，它们都是实数。且设则特征值与之间的关系为：，。设， 的正特征值，的正特征值，称，是的正奇异值，简称奇异值。若是正规矩阵，则的奇异值是的非零特征向量的模长。3.2 矩阵的奇异值分解（SVD）若，是的个正奇异值，则存在阶酉矩阵和阶酉矩阵，满足：其中，为奇异对角阵。满足是对角阵，满足是对角阵。的第列为的对应于奇异值对应的左奇异向量，的第列为的对应于奇异值对应的右奇异向量。它们的每一

7、列均为单位向量，且各列之间相互正交。若，是的个正奇异值，则总有次酉矩阵，满足：，其中。奇异值分解是一种基于特征向量的矩阵变换方法。奇异值分解是现代数值的最基本和最重要的工具之一。四、奇异值分解的图像性质任意一个矩阵的奇异值是唯一的，它刻画了矩阵数据的分布特征。直观上，可以这样理解矩阵的奇异值分解：将矩阵看成是一个线性变换，它将维空间的点映射到维空间。经过奇异值分解后，这种变换被分割成3个部分，分别为、和，其中和都是标准正交矩阵，它们对应的线性变换就相当于对维和维坐标系中坐标轴的旋转变换。若为数字图像，则可视为二维时频信息，可将的奇异值分解公式写为：其中，和分别是和的列矢量，是的非零奇异值。故上

8、式表示的数字图像可以看成是个秩为1的子图叠加的结果，而奇异值为权系数。所以也表示时频信息，对应的和可分别视为频率矢量和时间矢量，因此数字图像中的时频信息就被分解到一系列由和构成的视频平面中。由矩阵范数理论， 奇异值能与向量2-范数和矩阵Frobenious-范数(F-范数)相联系。若以F-范数的平方表示图像的能量，则由矩阵奇异值分解的定义知：。也就是说，数字图像经奇异值分解后，其纹理和几何信息都集中在、之中，而中的奇异值则代表图像的能量信息。性质1：矩阵的奇异值代表图像的能量信息，因而具有稳定性。设，是矩阵的一个扰动矩阵。和的非零奇异值分别记为：和。且，是的最大奇异值。则有：。由此可知，当图像

9、被施加小的扰动时，图像矩阵的奇异值变化不会超过扰动矩阵的最大奇异值，所以图像奇异值的稳定性很好。性质2：矩阵的奇异值具有比例不变性。设，矩阵的奇异值为，矩阵()的奇异值为。则有：。性质3：矩阵的奇异值具有旋转不变性。设，矩阵的奇异值为，。若是酉矩阵，则矩阵的奇异值与矩阵的奇异值相同：。性质4：设，。若，所以可得：上式表明，在F-范数意义下，是在空间（秩为的维矩阵构成的线性空间）中的一个将秩最佳逼近。因此可根据需要保留个大于某个阈值的而舍弃其余个小于阈值的且保证两幅图像在某种意义下的近似。这就为奇异值特征矢量的降维和数据压缩等应用找到了依据。五、图像的奇异值分解压缩方法5.1 奇异值分解压缩原理

10、分析用奇异值分解来压缩图像的基本思想是对图像矩阵进行奇异值分解，选取部分的奇异值和对应的左、右奇异向量来重构图像矩阵。根据奇异值分解的图像性质1和4可以知道，奇异值分解可以代表图像的能量信息，并且可以降低图像的维数。如果表示个维向量，可以通过奇异值分解将表示为个维向量。若的秩远远小于和，则通过奇异值分解可以大大降低的维数。对于一个像素的图像矩阵，设，其中，。按奇异值从大到小取个奇异值和这些奇异值对应的左奇异向量及右奇异向量重构原图像矩阵。如果选择的，这是无损的压缩；基于奇异值分解的图像压缩讨论的是，即有损压缩的情况。这时，可以只用个数值代替原来的个图像数据。这个数据分别是矩阵的前个奇异值， 左

11、奇异向量矩阵的前列和 右奇异向量矩阵的前列元素。比率：称为图像的压缩比。显然，被选择的奇异值的个数应该满足条件，即。故在传送图像的过程中，不需要传个数据，而只需要传个有关奇异值和奇异向量的数据即可。接收端，在接收到奇异值以及左奇异向量和右奇异向量后，可以通过：重构出原图像矩阵。与的误差为：某个奇异值对图像的贡献可以定义为，对一幅图像来说，较大的奇异值对图像信息的贡献量较大，较小的奇异值对图像的贡献较小。假如接近1，该图像的主要信息就包含在之中。通常图像的奇异值都具“大L 曲线”，只有不多的一些比较大的奇异值，其它的奇异值相对较小，因此一般只需要比较小的k 就使接近1。在满足视觉要求的基础上，按

12、奇异值的大小选择合适的奇异值个数，就可以通过将图像恢复。越小，用于表示的数据量就小，压缩比就越大，而越接近，则与就越相似。在一些应用场合中，如果是规定了压缩比，则可以由式求出，这时也同样可以求出。5.2 奇异值分解压缩应用过程在对图像进行操作时，因为矩阵的维数一般较大，直接进行奇异值分解运算量大， 可以将图像分解为子块，对各子块进行奇异值分解并确定奇异值个数，将每个子块进行重构。这样操作除了因为对较小型的矩阵进行奇异值分解的计算量比较小外，另一方面是为了利用原始图像的非均匀的复杂性。如果图像的某一部分比较简单，那么只需要少量的奇异值，就可以达到满意的近似效果。为了保证图像的质量就需要较多的奇异

13、值。但是各个子块的奇异值数目， 大小各不相同， 因此可以考虑为每个子块自适应的选择适当的奇异值数目。一种简单的方法是定义奇异值贡献量的和来选择，其中是一个接近1的数。对常见的256 ×256 .bmp格式的图像（位图），划分为4×4个子块，每个子块大小为64×64。对每个子块根据来选择所需要的奇异值数目。增大的值来选择奇异值数目，可以推理得随着不断增大，视觉效果越来越好。随着不断增大， 需要的奇异值也增多， 压缩比会减小。六、小 结经过以上讨论可知，用奇异值分解进行图像压缩，肯定能取得成功，也具有较好的应用价值，但仍然需要有以下值得去思考并改善：1、对子块的划分可以采取更加有效的方法来完成。例如对规模很大的矩阵，随机抽取矩阵的某些行列得到规模较小的矩阵，计算小矩阵的奇异值，重复若干次，用这些小矩阵的奇异值逼近原始矩阵的奇异。2、影响运算速度的因素是SVD 变换运算比较大，能否找到一个快速的SVD 变换算法。另外，若已知图像矩阵的奇异值及其特征空间，一般认为较大的奇异值及其对应的奇异向量表示图像信号，而噪声反映在较小的奇异值及其对应的奇异向量