图像的分形维数研究与分析

　　分形维数简介

　　分形维数被誉为大自然的几何学的分形（Fractal）理论，是现代数学的一个新分支，但其本质却是一种新的世界观和方法论。分维反映了复杂形体占有空间的有效性，它是复杂形体不规则性的量度。它与动力系统的混沌理论交叉结合，相辅相成。它承认世界的局部可能在一定条件下或过程中，在某一方面（形态，结构，信息，功能，时间，能量等）表现出与整体的相似性，它承认空间维数的变化既可以是离散的也可以是连续的，进而拓展了视野。

　　分形维数原理

　　fractal dimension主要描述分形最主要的参量，简称分维。通常欧几里德几何中，直线或曲线是1维的，平面或球面是2维的，具有长、宽、高的形体是 3 维的；然而对于分形如海岸线、科赫曲线、谢尔宾斯基海绵等的复杂性无法用维数等于 1、2、3 这样的数值来描述。科赫曲线第一次变换将1英尺的每边换成3个各长4英寸的线段，总长度变为 3×4×4/3=16 英寸；每一次变换使总长度变为乘以4/3，如此无限延续下去，曲线本身将是无限长的。这是一条连续的回线，永远不会自我相交，回线所围的面积是有限的，它小于一个外接圆的面积。因此科赫曲线以它无限长度挤在有限的面积之内，确实是占有空间的 ，它比1维要多，但不及2维图形，也就是说它的维数在1和2之间，维数是分数。同样，谢尔宾斯基海绵内部全是大大小小的空洞，表面积是无限大，而占有的 3 维空间是有限的，其维数在2和3之间。

　　图像的分形维数研究与分析

　　1 分形基础

　　1. 1 分形理论

　　Mandelbrot 曾经为分形下过两个定义：（1） D（ A）》d（ A的集合 A，称为分形集，其中 ，D（ A）为集合A的 Hausdoff 维数（或分维数），d（ A）为其拓扑维数。一般说来 ，D（ A）不是整数，而是分数;（2）部分与整体以某种形式相似的形 ，称为分形。然而，经过理论和应用的检验，人们发现这两个定义很难包括分形如此丰富的内容。实际上，对于什么是分形，到目前为止还不能给出一个确切的定义，人们通常是列出分形的一系列特性来加以说明。 分形特性有：

　　（1）分形集具有精细的结构 ，即有任意小比例的细节 ;

　　（2）分形集是如此的不规则以至它的整体和局部都不能用传统的几何语言描述 ;

　　（3）分形集通常有某种自相似的性质 ，可能是近似的或是统计的 ;

　　（4）一般的 ，分形集的“分形维数”（以某种方式定义） 大于它的拓扑维数 ;

　　（5）在大多数令人感兴趣的情形下 ，分形集以非常简单的方法定义 ，可能由迭代产生。

　　1. 2 分形维数

　　在欧几里德空间 ，所有的物体的维数都是整数维的 ，而分形理论则把维数视为分数。 维数是分形几何理论中重要基本概念之一 ，目前已经有许多维数的定义 ，主要包括 Hausdorff 维数、计盒维数（Box2counting dimension） 、修正计盒维数、填充维数等。 在实际应用中 ，通常所说的分形维数都是指计盒维数 ，而在分形理论应用研究中提出的许多维数的概念都是计盒维数的变形。

　　2 分形维数在图像分析中的应用

　　分形理论在图像分析中的应用是基于以下两点为基础的 ：一是自然界中不同种类的形态物质一般具有不同的分维数 ;二是由 A. P. Pentland 的假设 ，自然界中的分形与图像的灰度表示之间存在着一定的对应关系

　　2. 1 分形维数在图像分类中的应用

　　李厚强等在对航空遥感图像进行分类时采用基于分形维数的方法计算纹理特征 ，以此来提取和分析遥感图像的结构特征 ，并与航空遥感图像的另一重要特征 ———光谱信息相结合 ，通过综合应用纹理信息和光谱信息实现对遥感图像的分析与识别。 尺度不变性是基于分形的纹理特征的一个重要特点 ，这非常适合遥感图像分类。 现有的基于分形模型的方法是以分形维数作为纹理特征 ，在计算分形维数时 ，该文章采用 N. Sarkar 和 B. B. Chaudhuri 提出的 DBC（Differential Box Counting） 方法。 其主要思想如下。

　　将 M ×M 大小的图像分割成 s ×s 的子块（ M/ 2 Ε s 》 1 ，s 为整数） ，令 r = s/ M 。 将图像想象成三维空间中的曲面 ， x ， y 表示平面位置 ， z 轴表示灰度值。 xy 平面被分割成许多s ×s的网格 ，在每个网格上 ，是一列 s ×s ×s的盒子。 设图像灰度在第（ i ， j）网格中的最小值和最大值分别落在第 k 和第 l 个盒子中 ，则 ：

　　nr （ i ， j） = l - k + 1 （4）

　　是覆盖第（ i ， j） 网格中的图像所需的盒子数 ，而覆盖整个图像所需的盒子数 N r 为

　　针对不同的 r，计算 N r ，应用最小二乘法 ，即可求得分形维数 D. 在实际应用中 ，仅以分形维数为特征 ，很难达到好的分类效果 ，于是 ，文章引入“空隙”特征的概念 ，并将高、低灰度图像变换作为纹理特征之一 ，结合图像的色度和饱和度 ，有效地避免了不同材料计算得出相同分形维数的情况 ，可以较为准确的对航空遥感图像进行分类。

　　基于分形维数的图像分类也可以应用于图像压缩的预处理 ，基于分形维数特征 ，将图像分为不同的区域 ，而这一分类 ，对图像进行进一步的压缩编码具有重要的指导意义。曹文伦等依据图像灰度变化的平缓或剧烈程度 ，将图像分为以下几类 ：纹理、纹理边缘、纹理边缘平滑、边缘平滑和平滑 ，并且分别计算以及确定了不同图像种类的分形维数区间。 按照此区间 ，可以对待测图像基于其分形维数进行分类。 利用分类结果 ，可以分析以及初步判断待测图像的空间不规则程度 ，对在进一步进行压缩编码具有指导意义。 如果一个物体的表面是分形的 ，则由它产生的图像灰度表面也具有分形的性质。 基于此 ，通过 Keller 的盒维数算法完成原图分形维数的提取。 其主要思想是：令离散灰度图像 F是一个在三维空间中的分形曲面 ，长、宽为 R; Nε（ F） 表示边长为ε×ε×ε的包含所要估计的图像区域 ，其分形维数 DB 将由（7） 式决定 ：

　　将 R ，ε带入（8） 式 ，即可求出灰度图像 F的盒维数 DB 。 分形维数可以很好的作为图像分析与分类的参数 ，该方法所使用的采用分形维数的图像分类与人凭借主观感受对图像进行的分类相吻合。

　　李杰等在对待压缩图像进行预处理时 ，把标准图像分为图像块 ，通过计算其分形维数 ，依据图像块内灰度及其变化情况 ，把图像块分为 3 类 ：均匀块、中值块和边缘块 ，不同种类的图像块反映了不同的纹理特征。 在对图像块分形维数特征的计算中 ，本文采用的是基于形态算子的毯子覆盖法计算图像块灰度图像的分形维数。 其主要思想是灰度图像想象成一个在三维空间中的分形曲面 ，所要估计的图像区域的分形维数 D 将由下面的公式确定

　　式中 A （ r） 为分形曲面的表面积 ， r为度量时所使用的面积元尺度 ， D为曲面的分形维数 ， F是一常数。 曲面的表面积求法如下 ：把三维空间中的所有距离曲面表面为 r的点 ，用厚度为2 r的毯子覆盖 ，曲面的表面积就是被覆盖的体积除以2 r，覆盖的毯子由曲面的上表面 ur 和下表面 br 来定义 ：设曲面上点（ i ， j） 的灰度值为 g（ i ， j） ，取 r = 1 ，2 ，3 ， …， ur ， br 的值为

　　其中点（ m ， n） 是距离点（ i ， j） 不大于 1 的点 ，取于点（ i ， j） 的四邻域。 该毯子的体积为 ：

　　并在 log［ A （ r） ］2log （ r） 坐标系中做直线的最小二乘拟合 ，直线的斜率即为图像的分形维数 D.文中引入了α，β两个衍生参数作为分类辨别参数 ，以此来集中体现图像块的纹理特点。 将待测图像块的分形维数值与α，β相比较 ，并确定图像块的所属分类 ，依据分类 ，对不同的图像采用不同方式的压缩编码。 该方法基本正确的反映了图像中的纹理分布情况。 针对图像块的纹理分布情况 ，采用不同的压缩编码 ，可以得到较好的效果.