基于复杂网络理论的图像分割方法与流程

本发明涉及一种图像分割处理方法，特别是涉及一种包含复杂纹理的图像分割方法，应用于计算机视觉或图像处理技术领域。

背景技术：

在计算机视觉领域，图像分割是重要的图像处理任务之一，并且很具有挑战性。图像分割技术在许多领域内都具有广泛的应用，例如在无人车驾驶、物体识别、医学图像处理等领域。目前比较常见的图像分割方法有许多种，包括基于边界信息的图像分割方法、基于机器学习的图像分割方法以及基于深度学习的图像分割方法。传统的图像分割方法对包含复杂纹理的图像很难进行分割；基于机器学习的分割算法往往需要为每个像素进行提取特征，并对每个像素点进行分类，因此，计算量以及存储消耗是非常大的；深度学习可以利用卷积神经网络提取出图像高层语义特征，实现图像的分割，但是深度学习都是基于大量图像样本进行训练且要标签数据，这在一定程度上也对其产生了限制。在一些科研领域内，如材料学、医学，其图像往往都是高分辨率但数据样本比较小，对图像分割提出了具大的挑战。

近年来，伴随着复杂网络研究的热潮，复杂网络理论在处理图像中所体现出来的性能逐渐显现，越来越多基于复杂网络理论的图像处理算法开始出现。研究人员发现利用复杂网络可以结合图像局部性质和整体性质，为图像处理提供更加全面的信息，有助于一些图像处理问题的解决，比如轮廓识别、纹理分析等。复杂网络中的社区结构是重要的拓扑结构属性之一，它揭示了复杂网络中隐藏的规律和行为特征。社区结构中要求，在社区内部节点之间的联系非常紧密，而社区之间的联系则相对比较稀疏。基于复杂网络理论社区结构中这一特征，可以把它和图像分割相关联，在图像中找出由不同的像素值所代表的不同的社区来实现图像分割任务。但现有的利用复杂网络理论处理包含复杂纹理的图像时效果还不够理想。

技术实现要素：

为了解决现有技术问题，本发明的目的在于克服已有技术存在的不足，提供一种基于复杂网络理论的图像分割方法，发展了复杂网络中相关理论，无需大量数据和标签，可以实现对复杂纹理图像的分割。

为达到上述发明创造目的，本发明采用如下技术方案：

一种基于复杂网络理论的图像分割方法，步骤如下：

a.对图像进行去噪声处理，然后基于图像中像素值的分布，利用聚类和网格化方法，在图像上生成网络节点集；

b.定义网络节点之间的相似性，生成网络边集，生成图像的网络拓扑结构；

c.利用模块度函数，优化网络拓扑结构，得到最优社区结构划分结果；

d.基于生成的图像的网络拓扑结构以及最优社区结构划分结果，进行区域轮廓的提取，获得各个图像分割区域，完成区域分割任务。

本发明采用降噪方法，优选均值滤波、中值滤波或高斯滤波方法，对输入的数字图像进行降噪。本发明采用聚类算法，优选采用k-means算法、dbscan算法或meanshift算法，对图像进行像素级别上的聚类，以实现对不同像素级目标的初步分割。

作为本发明优选的技术方案，在所述步骤a中，通过降噪方法，对输入的数字图像进行降噪；

然后根据图像中像素值的分布，采用聚类算法，对图像进行像素级别上的聚类，对不同像素级目标进行初步分割；

然后对图像进行网格化，从而将图像用多个大小为d*d的方格进行划分，d为网格的方格行数或列数，计算每个方格内不同类别的像素点数，每个方格用像素点数多的类型标记成此方格的类型，并记录此方格为一个节点，方格中心点坐标表示节点的位置，把所有节点记录在一个节点集v={v1,v2,···,vn}中，每个节点都有一个坐标值，即vi(xi,yi)；通过此步骤，获得图像上的网络节点集。

作为本发明优选的技术方案，在所述步骤b中，构建的生成图像的网络拓扑结构的步骤如下：

复杂网络用节点及节点之间的边进行表示的，即g=(v,e)，其中v={v1,v2,···,vn}为节点集，e={e1,e2,···,em}为边集；

定义集合v以及e，得到网络拓扑结构，通过所述步骤a方法，获得节点集v；

两个节点vi(xi,yi)和vj(xj,yj)相似性是通过节点之间归一化后的欧式距离来度量的，其中欧式距离di,j如下所示：

,(1)

归一化后的欧式距离同样用di,j表示，公式如下：

,(2)

其中md表示任意两个节点之间距离中的最大值；

引入半径阈值r∈[0,1]，用来判断两节点之间是否应该存在连边，如下列公式所示：

(3)

为了显示网络中节点连接的紧密程度，引入连接密度来度量节点之间的连接的紧密程度:

(4)

其中，是节点度，表示与节点相连边的个数；是节点强度，表示与节点相连边上的权值之和；把连接密度归一化到区间[0,1]中，公式如下：

(5)

利用连接密度，从网络拓扑中筛选出连接紧密的节点，即连接密度大的节点；

因此，一个连接密度阈值t被定义，若两个节点的连接密度同时大于阈值时，则保留原本连接，否则断开；从而利用本步骤以上方法，构建图像的网络拓扑结构。

作为本发明优选的技术方案，在所述步骤c中，进行优化网络拓扑结构，得到最优社区结构划分结果，采用如下步骤；

在得到所述步骤b中图像的网络拓扑结构时，对于一个个密集的区域，网络是不连通的，每一个连通区域表示一个社区结构；在不同的阈值参数下，得到不同的社区划分结果；为了找到最优的r和t，定义社区模块度如下：

(6)

其中，是总强度之和，满足以下公式：

(7)

另外，和分别代表了节点和节点所属的社区，并得到以下关系：

(8)

公式（6）主要统计了和差的累加和，和越大，则说明社区划分的越好，即区域分割的越合理，因此需要找到合适的r和t，使得公式（6）中的q最大，从而得到最优社区结构划分结果。

作为本发明优选的技术方案，在所述步骤d中，进行区域轮廓的提取和区域分割时，采用如下步骤：

通过在所述步骤b中复杂网络的建立和在所述步骤c中公式（6）的计算，对图片中的拓扑结构进行分割，从一个完整的拓扑网络分割成多个子图，每个子图就是一个分割区域，然后利用填充算法，对各个子图内的空洞进行填充，最后利用边缘检测和轮廓提取，获得各个图像分割区域，完成区域分割任务。

本发明与现有技术相比较，具有如下显而易见的突出实质性特点和显著优点：

1.本发明方法能够分割带有复杂纹理的图像，能够获得复杂纹理图像中关键部位的形状；

2.本发明方法能优化网络拓扑结构中的社区结构，实现从社区结构划分到图像区域分割的映射；

3.本发明方法无需大量数据和标签，就能实现对具有复杂纹理图像的分割。

附图说明

图1为本发明各实施例基于复杂网络理论的图像分割方法的流程示意图。

具体实施方式

以下结合具体的实施例子对上述方案做进一步说明，本发明的优选实施例详述如下：

实施例一：

在本实施例中，参见图1，一种基于复杂网络理论的图像分割方法，步骤如下：

a.对图像进行去噪声处理，然后基于图像中像素值的分布，利用聚类和网格化方法，在图像上生成网络节点集；本实施例对图像进行预处理，包括图像去噪，生成目标像素区域，对图像进行网格化，最终生成网络节点集，为更有效利用复杂网络来结合图像局部性质和整体性质打好基础；

b.定义网络节点之间的相似性，生成网络边集，生成图像的网络拓扑结构；

c.利用模块度函数，优化网络拓扑结构，得到最优社区结构划分结果；本实施例使图像按社区结构的性质，生成的拓扑结构划分社区，并得到最优社区结构划分结果，本实施例采用复杂网络中的社区结构，能揭示了复杂网络中隐藏的规律和行为特征；

d.基于生成的图像的网络拓扑结构以及最优社区结构划分结果，进行区域轮廓的提取，获得各个图像分割区域，完成区域分割任务。社区结构中要求，在社区内部节点之间的联系非常紧密，而社区之间的联系则相对比较稀疏。本实施例基于复杂网络理论社区结构中这一特征，将其和图像分割相关联，在图像中找出由不同的像素值所代表的不同的社区来实现图像分割任务。

本实施例基于复杂网络理论的图像分割方法，发展了复杂网络中相关理论，实现对复杂纹理图像的分割。首先对图像进行去噪，利用聚类获得目标像素区域，生成网络节点集；然后，对图像进行网格化，获得网络的拓扑结构；之后，优化网格结构得到最优社区结构划分结果；最后，根据社区结构的最优划分，结合区域轮廓提取和区域分割实现图像分割任务。本发明中提到的图像分割方法可应用于具有复杂纹理的数字图像分割，因此，本发明可以应用在多个学科领域，如材料学、医学、微观结构学、天文学等，提高研究速度，降低研究成本，促进这些学科的研究和发展。

实施例二：

本实施例与实施例一基本相同，特别之处在于：

在本实施例中，参见图1，在所述步骤a中，通过降噪方法，并采用均值滤波、中值滤波或高斯滤波方法，对输入的数字图像进行降噪；

然后根据图像中像素值的分布，利用聚类算法，并采用k-means算法、dbscan算法或meanshift算法，对图像进行像素级别上的聚类，对不同像素级目标进行初步分割；

然后对图像进行网格化，从而将图像用多个大小为d*d的方格进行划分，d为网格的方格行数或列数，计算每个方格内不同类别的像素点数，每个方格用像素点数多的类型标记成此方格的类型，并记录此方格为一个节点，方格中心点坐标表示节点的位置，把所有节点记录在一个节点集v={v1,v2,···,vn}中，每个节点都有一个坐标值，即vi(xi,yi)；通过此步骤，获得图像上的网络节点集。本实施例通过一次扫描下来即可获得网络中的所有节点，生成网络节点集。

实施例三：

本实施例与前述实施例基本相同，特别之处在于：

在本实施例中，参见图1，在所述步骤b中，构建的生成图像的网络拓扑结构步骤如下：

复杂网络用节点及节点之间的边进行表示的，即g=(v,e)，其中v={v1,v2,···,vn}为节点集，e={e1,e2,···,em}为边集；

定义集合v以及e，得到网络拓扑结构，通过所述步骤a方法，获得节点集v；

两个节点vi(xi,yi)和vj(xj,yj)相似性是通过节点之间归一化后的欧式距离来度量的，其中欧式距离di,j如下所示：

,(1)

归一化后的欧式距离同样用di,j表示，公式如下：

,(2)

其中md表示任意两个节点之间距离中的最大值；

引入半径阈值r∈[0,1]，用来判断两节点之间是否应该存在连边，如下列公式所示：

(3)

为了显示网络中节点连接的紧密程度，引入连接密度来度量节点之间的连接的紧密程度:

(4)

其中，是节点度，表示与节点相连边的个数；是节点强度，表示与节点相连边上的权值之和；把连接密度归一化到区间[0,1]中，公式如下：

(5)

利用连接密度，从网络拓扑中筛选出连接紧密的节点，即连接密度大的节点；

因此，一个连接密度阈值t被定义，若两个节点的连接密度同时大于阈值时，则保留原本连接，否则断开；从而利用本步骤以上方法，构建图像的网络拓扑结构。

实施例四：

本实施例与前述实施例基本相同，特别之处在于：

在本实施例中，参见图1，在所述步骤c中，进行优化网络拓扑结构，得到最优社区结构划分结果，采用如下步骤；

在得到所述步骤b中图像的网络拓扑结构时，对于一个个密集的区域，网络是不连通的，每一个连通区域表示一个社区结构；在不同的阈值参数下，得到不同的社区划分结果；为了找到最优的r和t，定义社区模块度如下：

(6)

其中，是总强度之和，满足以下公式：

(7)

另外，和分别代表了节点和节点所属的社区，并得到以下关系：

(8)

公式（6）主要统计了和差的累加和，和越大，则说明社区划分的越好，即区域分割的越合理，因此需要找到合适的r和t，使得公式（6）中的q最大，从而得到最优社区结构划分结果。

本实施例得到网络拓扑时一个个密集的区域，网络是不连通的，每一个连通区域表示一个社区结构。在不同的阈值参数下可以得到不同的社区划分结果。为了找到最优的r和t，特别定义社区模块度，优化网络结构得到最优分割结果。

实施例五：

本实施例与前述实施例基本相同，特别之处在于：

在本实施例中，参见图1，在所述步骤d中，进行区域轮廓的提取和区域分割时，采用如下步骤：

通过在所述步骤b中复杂网络的建立和在所述步骤c中公式（6）的计算，对图片中的拓扑结构进行分割，从一个完整的拓扑网络分割成多个子图，每个子图就是一个分割区域，然后利用填充算法，对各个子图内的空洞进行填充，最后利用边缘检测和轮廓提取，获得各个图像分割区域，完成区域分割任务。

本实施例方法能够分割带有复杂纹理的图像，能够获得复杂纹理图像中关键部位的形状；可以优化网络拓扑结构中的社区结构，实现从社区结构划分到图像区域分割的映射；本实施例方法无需大量数据和标签，就能实现对具有复杂纹理图像的分割。社区结构中要求，在社区内部节点之间的联系非常紧密，而社区之间的联系则相对比较稀疏。本实施例方法基于复杂网络理论社区结构中这一特征，可以把它和图像分割相关联，在图像中找出由不同的像素值所代表的不同的社区来实现图像分割任务。

综上所述，本发明上述实施例基于复杂网络理论的图像分割方法的流程示意图，共分为以下步骤：首先，对图像进行预处理，包括图像去噪、生成目标像素区域，对图像进行网格化，生成网络节点集；然后，定义和计算节点之间的相似性，生成网络边集，得到网络的拓扑结构；之后，图像按社区结构的性质，生成的拓扑结构划分社区，并得到最优分割结果；最后，区域轮廓的提取和区域分割的实现。本发明上述实例方法采用图像分割方法，是针对具有复杂纹理的数字图像提出了一种基于复杂网络理论的图像分割方法，可应用于具有复杂纹理的数字图像分割，因此，本发明可以应用在多个学科领域，如材料学、医学、微观结构学、天文学等，提高研究速度，降低研究成本，促进这些学科的研究和发展。

上面对本发明实施例结合附图进行了说明，但本发明不限于上述实施例，还可以根据本发明的发明创造的目的做出多种变化，凡依据本发明技术方案的精神实质和原理下做的改变、修饰、替代、组合或简化，均应为等效的置换方式，只要符合本发明的发明目的，只要不背离本发明基于复杂网络理论的图像分割方法的技术原理和发明构思，都属于本发明的保护范围。