CN103824283B - 基于随机概率采样与多级融合的大幅面地形图分割方法 - Google Patents
基于随机概率采样与多级融合的大幅面地形图分割方法 Download PDFInfo
- Publication number
- CN103824283B CN103824283B CN201410029317.8A CN201410029317A CN103824283B CN 103824283 B CN103824283 B CN 103824283B CN 201410029317 A CN201410029317 A CN 201410029317A CN 103824283 B CN103824283 B CN 103824283B
- Authority
- CN
- China
- Prior art keywords
- matrix
- segmentation
- image
- cluster
- land form
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种基于随机概率采样与多级融合的大幅面地形图分割方法:步骤1,输入原始地形图图像;步骤2,对原始地形图图像进行随机概率采样;步骤3,确定地形图分割种类数;步骤4,计算聚类中心,得到聚类中心矩阵;步骤5,计算原始地形图的隶属矩阵,得到隶属矩阵;步骤6,利用隶属矩阵对原始地形图进行模糊分类,得到分割后的图像;步骤7,对分割后的图像进行多级融合,得到的融合后图像;步骤8,输出分割后的分版图像。该方法将随机概率采样与图像多级融合相结合以进行地形图分割,最终既能够更准确地获得地形图的分割版图,又能大幅度提高分割效率。尤其适用于大幅面的地形图分割。
Description
技术领域
本发明属于地理信息处理技术领域,具体涉及一种基于随机概率采样与多级融合的大幅面地形图分割方法。本发明应用于大幅面地形图的图像分割。
背景技术
在地形图处理领域中,为了快速准确地对地形图进行分割得到分割版图,并为地形图中地理要素的进一步提取提供有利条件,通常需要对地形图分割。目前,地形图分割方法主要采用基于颜色空间的聚类和基于样本学习的分类法。
郑华利等人在文献“郑华利,周献中,王建宇,‘基于色彩空间转换及模糊约束聚类的地形图自动分色’[J],测绘学报,32(2),183-187(2003)”中提出了一种模糊约束聚类算法。该方法首先对地形图的色彩空间进行转换,使其直方图的峰谷特性更加明显,然后对模糊C均值聚类算法进行模糊约束,从而较好地解决了地图的自动分色问题。由于经过色彩空间的转换后直方图的峰值比较明显,所以该算法对于初始聚类中心选择的依赖性较小。此外,采用了模糊约束的颜色聚类,对于那些在地形图直方图中占有较小比例的颜色信息也能正确地进行分层。该方法虽然能够得到较好的地形图分割效果,但是仍然存在的缺点是,模糊C均值聚类算法对地形图进行4分类,没有考虑地形图中存在的过渡色、渐变色和覆盖色等复杂问题,导致分割结果不准确。此外,该方法对地形图中所有数据进行一次性处理,没有考虑到大幅面地形图数据量大的问题,因此分割效率不高。
Alireza Khotanzad等人在文献“Alireza Khotanzad and Edmund Zink,‘Contour Line and Geographic Feature Extraction from USGS Color TopographicalPaper Maps’[J],IEEE Transactions on Pattern Analysis and MachineIntelligence,25(1),18-31(2003)”中提出了一种基于样本学习的分类方法对地形图进行分割。该方法首先根据地形图中不同地理要素之间的重叠覆盖情况,制作相应的模板,然后根据模板的颜色特征与地形图中像素的颜色特征之间的差异对地形图进行分割。此方法虽然能够较好的处理面状地理要素信息以及地形图中线划要素信息。但是仍然存在的缺点是,地形图中不同地理要素的模板样本难以制作齐全,此外,图像分割时仅仅使用颜色距离,难以很好的处理地形图中存在的过渡色、渐变色和覆盖色等复杂问题,这些因素都是导致最终分割结果不准确的原因,从而不能得到较好的地形图分割效果。
综上所述,基于颜色空间的聚类和基于样本学习的分类方法在地形图分割方面虽然获得了较好的效果,但这些方法都难以快速准确地进行大幅面地形图图像的分割,此外,在地形图分割种类数的确定方面没有一个准确的自适应的确定方法。所以在利用传统的基于颜色空间的聚类和基于样本学习的分类方法得到的分割版图中,存在分割不准确以及效率低下的问题。
发明内容
针对上述现有技术中基于颜色空间的聚类和基于样本学习的分类方法导致的地形图分割效果不理想和效率低下的缺点,本发明的目的在于,提出一种基于随机概率采样与多级融合的大幅面地形图分割方法,该方法中对地形图的分割主要包括获得聚类中心,模糊分类和多级融合三部分。本发明充分利用随机概率采样能够较大程度减少获得聚类中心所需要的数据量的优点,以及多级融合能够更准确地对分割后图像进行融合的特点,将随机概率采样与图像多级融合相结合以进行地形图分割,最终既能够更准确地获得地形图的分割版图,又能大幅度提高分割效率。尤其适用于大幅面的地形图分割。
为了达到上述目的,本发明采用如下技术方案予以解决:
一种基于随机概率采样与多级融合的大幅面地形图分割方法,具体包括如下步骤:
步骤1,输入原始地形图图像;
步骤2,对原始地形图图像进行随机概率采样,得到采样数据矩阵;
步骤3,确定地形图分割种类数;
步骤4,计算聚类中心,得到聚类中心矩阵;
步骤5,计算原始地形图的隶属矩阵,得到隶属矩阵;
步骤6,利用隶属矩阵对原始地形图进行模糊分类,得到分割后的图像;
步骤7,对分割后的图像进行多级融合,得到的融合后图像;
步骤8,输出分割后的分版图像。
进一步的,所述步骤2的具体步骤如下:
2a)生成一个符合正态分布的,且和原始地形图A等大小的随机矩阵RM×N作为随机概率矩阵,该矩阵中所有元素的值为0和1之间的任意随机数;其中,M为原始地形图的宽度,N为原始地形图的长度;
2b)计算概率阈值P,并利用步骤2a)生成的随机概率矩阵RM×N对原始地形图进行数据采样;概率阈值P的计算如下:
其中,M为原始地形图的宽度,N为原始地形图的长度。
数据采样过程如下:遍历随机概率矩阵RM×N,如果矩阵中某一元素的值大于等于概率阈值P,则原始地形图A中相同位置像素点的颜色信息被选择为采样信息;否则原始地图A中相同位置像素点的颜色信息为0;遍历结束后得到采样数据矩阵A';
进一步的,所述步骤3中,利用下式确定地形图的分割种类数:
Cluster_n=Cl 1+Cl 2+Cl 3+…Cl l+1. (5)
其中,Cluster_n为地形图的分割种类数,Cl 1为排列组合表达式,表示在由l种颜色绘制的地形图中的某区域有1种颜色相互覆盖所得到的颜色种类数,l≥1。
进一步的,所述步骤4的具体步骤如下:
4a)初始化隶属矩阵:初始化一个与原始地形图A相同大小的随机矩阵作为初始隶属矩阵U0;并使其满足下式:
其中,u0,k(i,j)为初始隶属矩阵U0中坐标(i,j)位置第k类元素的值;
4b)利用下式计算得到矩阵大小为Cluster_n×3的聚类中心矩阵vI,其中,I为迭代次数,当迭代次数I=1时,隶属度矩阵U=U0,uk(i,j)=u0,k(i,j);
其中,vI,k为第I次迭代中聚类中心矩阵vI的第k个聚类中心,k=1,…,Cluster_n;m=2;x(i,j)为原始地形图中坐标(i,j)位置的像素值;(uk(i,j))m为隶属矩阵U中坐标(i,j)位置第k类元素值的m次方;
4c)根据下式计算第I次迭代的总体价值函数JI:
其中,JI,k为第I次迭代的第k类价值函数;
如果总体价值函数JI小于阀值Jmin=10-5,或总体价值函数JI相对上次得到的价值函数值JI-1的改变量小于阀值Jdif_min=10-5,则步骤4的计算聚类中心结束,最终得到聚类中心矩阵v=vI,然后执行步骤5;否则,进行步骤4d);
4d)使用下式计算新的隶属度矩阵U';
其中,u'k(i,j)为新的隶属矩阵U'中坐标(i,j)位置第k类元素的值;vI,k和vI,p分别表示聚类中心矩阵vI中第k个和第p个聚类中心;
令U=U',uk(i,j)=u'k(i,j),I=I+1,返回步骤4b)。
进一步的,所述步骤5的具体步骤如下:
5a)读取步骤(4)中得到的聚类中心矩阵v;
5b)利用模糊C均值聚类算法中的隶属度计算公式,计算原始地形图A中各像素的隶属度,得到所有像素的隶属度组成隶属矩阵UF;
其中,uF,k(i,j)为隶属矩阵UF中位置为(i,j)处第k类元素的值,对应于原始地形图中坐标位置的像素属于k类的隶属度;vk和vp分别为聚类中心矩阵v的第k个和第p个聚类中心;x(i,j)为原始地形图中坐标(i,j)位置的像素值;m=2。
进一步的,所述步骤6的具体步骤如下:
利用隶属矩阵UF,按以下公式所示的模糊分类方法对原始地形图A进行图像分割得到分割后的图像:
x(i,j)∈imagek如果uF,k(i,j)≥uF,p(i,j)(p=1,2,…Cluster_n)
其中,imagek为原始地形图中像素x(i,j)被划入的第k个分割后的图像。
进一步的,所述步骤7的具体步骤如下:
7a)读取聚类中心矩阵;
7b)计算每个聚类中心和其他聚类中心之间的距离;
7c)利用聚类中心之间的距离以及多级融合方法对分割后的图像进行融合,直至融合后得到的分版图像数目和地形图的分割种类数Cluster_n相等时停止,得到融合后的分版图像;
进一步的,所述步骤7b)中的距离采用马氏距离:
其中,Distk1,k2为聚类中心矩阵v中第k1个聚类中心vk1和第k2个聚类中心vk2之间的马氏距离,()T为矩阵的转置。
进一步的,所述步骤7c)中的多级融合方法的规则:当且仅当与第i类的最小距离是第j类,同时与第j类的最小距离是第i类时,将第i类和第j类融合为一新类。
与现有技术相比,本发明具有以下优点:
第一,在获得聚类中心时,本发明应用随机概率采样的方法对原始地形图进行数据采样,能够较大程度上减少计算聚类中心所需要的数据量,因而能够快速获得聚类中心,以提高分割效率。克服了现有技术中对大幅面地形图分割效率低下的缺点,使得使用本发明中的方法能够快速获得大幅面地形图的分割结果。
第二,本发明在确定分割种类数时,采用基于地形图颜色覆盖的自适应确定方法,能够根据地形图中颜色的种类数及其覆盖情况自适应确定图像分割种类数,因而能够得到更准确的分割种类数。
第三,本发明在对分割后图像进行融合处理时,应用了多级融合的方法。由于多级融合采用了基于马氏距离的相互融合规则,克服了现有技术中融合结果不准确的问题,使得使用本发明中的方法能够获得更准确的地形图分版图像。
附图说明
图1为本发明的流程图。
图2为本发明用于对小幅面地形图进行分割的仿真图。其中,(a)是原始小幅面地形图图像;(b1)是本发明方法得到的绿版图;(b2)是基于地形图空间信息的图像分割方法得到的绿版图;(b3)是基于遗传模糊C均值聚类方法得到的绿版图;(b4)是基于改进的模糊C均值聚类方法得到的绿版图;(c1)是本发明方法得到的棕版图;(c2)是基于地形图空间信息的图像分割方法得到的棕版图;(c1’)是c1局部放大图;(c2’)是c2局部放大图;(c3’)是c3局部放大图;(c4’)是c4局部放大图;(c3)是基于遗传模糊C均值聚类方法得到的棕版图;(c4)是基于改进的模糊C均值聚类方法得到的棕版图;(d1)是本发明方法得到的黑版图;(d2)是基于地形图空间信息的图像分割方法得到的黑版图;(d1’)是d1局部放大图;(d2’)是d2局部放大图,(d3’)是d3局部放大图;(d4’)是d4局部放大图;(d3)是基于遗传模糊C均值聚类方法得到的黑版图;(d4)是基于改进的模糊C均值聚类方法得到的黑版图;(e1)是本发明方法得到的蓝版图;(e2)是基于地形图空间信息的图像分割方法得到的蓝版图;(e1’)是e1局部放大图;(e2’)是e2局部放大图;(e3’)是e3局部放大图;(e4’)是e4局部放大图;(e3)是基于遗传模糊C均值聚类方法得到的蓝版图;(e4)是基于改进的模糊C均值聚类方法得到的蓝版图。
图3为本发明用于对大幅面地形图进行分割的仿真图。其中,(a)是原始地形图图像;(b)是本发明的方法得到的黑版图;(c)是本发明方法得到的蓝版图;(d)是本发明的方法得到的棕版图;(e)是本发明的方法得到的红版图;(f)是本发明的方法得到的绿版图。
以下结合附图和具体实施方式对本发明进一步解释说明。
具体实施方式
参照图1,本发明的基于随机概率采样与多级融合的大幅面地形图分割方法,具体包括如下步骤:
步骤1,输入原始地形图图像A
在计算机中应用matlab软件读取存储在计算机硬盘空间中的原始地形图A。
步骤2,随机概率采样
为了减少聚类算法在计算聚类中心的数据量,先对原始地形图进行随机概率采样,具体包括如下步骤:
2a)利用matlab软件生成一个符合正态分布的,且和原始地形图A等大小的随机矩阵RM×N作为随机概率矩阵,该矩阵中所有元素的值为0和1之间的任意随机数。其中,M为原始地形图的宽度,N为原始地形图的长度。
2b)计算概率阈值P,并利用步骤2a)生成的随机概率矩阵RM×N对原始地形图进行数据采样。概率阈值P的计算如下:
其中,M为原始地形图的宽度,N为原始地形图的长度。
数据采样过程如下:遍历随机概率矩阵RM×N,如果矩阵中某一元素的值大于等于概率阈值P,则原始地形图A中相同位置像素点的颜色信息(即该像素点的R,G,B值)被选择为采样信息;否则原始地图A中相同位置像素点的颜色信息为0;遍历结束后得到采样数据矩阵A'。即:
其中,A'(i,j)为采样后图像中坐标(i,j)位置的像素值,A(i,j)为原始地形图中坐标(i,j)位置的像素值,R(i,j)为随机概率矩阵RM×N中坐标(i,j)位置的元素值,1≤i≤M,1≤j≤N。
步骤3,地形图分割种类数的自适应确定
利用地形图分割种类数的自适应方法确定地形图的分割种类数。
Cluster_n=Cl 1+Cl 2+Cl 3+…Cl l+1. (5)
其中,Cluster_n为地形图的分割种类数,Cl 1为排列组合表达式,表示在由l种颜色绘制的地形图中的某区域有1种颜色相互覆盖所得到的颜色种类数,l≥1,Cl 2、Cl 3、…Cl l同理依次类推。地形图中黑色和其他所有颜色混合都为黑色,因此最后加1。
步骤4,计算聚类中心,得到聚类中心矩阵
利用模糊C均值聚类算法对步骤2b)得到采样数据矩阵A'进行聚类处理,得到和分割种类数Cluster_n相同个数的聚类中心矩阵v,并存入计算机。v是一个大小为Cluster_n×3的矩阵。具体包括以下步骤:
4a)初始化隶属矩阵:利用matlab软件初始化一个与原始地形图A相同大小的随机矩阵作为初始隶属矩阵U0,初始隶属矩阵U0中所有元素的值为0和1之间的任意随机数。并使其满足下式:
其中,u0,k(i,j)为初始隶属矩阵U0中坐标(i,j)位置第k类元素的值。
4b)利用下式计算得到矩阵大小为Cluster_n×3的聚类中心矩阵vI,其中,I为迭代次数,当迭代次数I=1时,隶属度矩阵U=U0,uk(i,j)=u0,k(i,j);
其中,vI,k为第I次迭代中聚类中心矩阵vI的第k个聚类中心,k=1,…,Cluster_n;m为权值参数,本发明中m=2;x(i,j)为原始地形图中坐标(i,j)位置的像素值;(uk(i,j))m为隶属矩阵U中坐标(i,j)位置第k类元素值的m次方。
4c)根据下式计算第I次迭代的总体价值函数JI:
其中,JI,k为第I次迭代的第k类价值函数。
如果总体价值函数JI小于阀值Jmin=10-5,或总体价值函数JI相对上次得到的价值函数值JI-1的改变量小于阀值Jdif_min=10-5,则步骤4的计算聚类中心结束,最终得到聚类中心矩阵v=vI,然后执行步骤5。否则,进行步骤4d)。
4d)使用下式计算新的隶属度矩阵U';然后令U=U',uk(i,j)=u'k(i,j),I=I+1,返回步骤4b)。
其中,u'k(i,j)为新的隶属矩阵U'中坐标(i,j)位置第k类元素的值;vI,k和vI,p分别表示聚类中心矩阵vI中第k个和第p个聚类中心。
步骤5,计算原始地形图的隶属矩阵
5a)读取步骤(4)中得到的聚类中心矩阵v;
5b)利用模糊C均值聚类算法中的隶属度计算公式,计算原始地形图A中各像素的隶属度,得到所有像素的隶属度组成隶属矩阵UF。
其中,uF,k(i,j)为隶属矩阵UF中位置为(i,j)处第k类元素的值,对应于原始地形图中坐标位置的像素属于k类的隶属度;vk和vp分别为聚类中心矩阵v的第k个和第p个聚类中心;x(i,j)为原始地形图中坐标(i,j)位置的像素值;m为权值参数,本发明中m=2;
步骤6,模糊分类,得到分割后的图像。具体操作如下:
利用步骤5b)中得到的隶属矩阵UF,按以下公式所示的模糊分类方法对原始地形图A进行图像分割得到分割后的图像。
x(i,j)∈imagek如果uF,k(i,j)≥uF,p(i,j)(p=1,2,…Cluster_n)
其中,imagek为原始地形图中像素x(i,j)被划入的第k个分割后的图像。
步骤7,多级融合,得到的融合后图像。具体操作如下:
7a)读取步骤(4)中得到的聚类中心矩阵v;
7b)利用下式计算每个聚类中心和其他聚类中心之间的距离,本发明中采用马氏距离公式:
其中,Distk1,k2为聚类中心矩阵v中第k1个聚类中心vk1和第k2个聚类中心vk2之间的马氏距离,()T为矩阵的转置;
7c)利用步骤7b)中得到的聚类中心之间的马氏距离以及多级融合方法对步骤6得到的分割后的图像进行融合,直至融合后得到的分版图像数目和地形图的分割种类数Cluster_n相等时停止,得到融合后的分版图像;
定义:多级融合方法的规则是:当且仅当与第i类的最小距离是第j类,同时与第j类的最小距离是第i类时,将第i类和第j类融合为一新类。即如果同时满足以下两个条件:
Disti,j=min({Disti,C|C=1,2,...,Cluster_n})
Disti,j=min({DistC,j|C=1,2,...,Cluster_n})
则
mapc(i,j)=imagek1(i,j)+imagek2(i,j)
其中,mapc(i,j)为第c个融合后的图像中坐标(i,j)位置像素值,imagek1(i,j)和imagek2(i,j)分别为第k1和第k2个分割后图像中坐标(i,j)位置像素值。
步骤8,输出分割后的分版图像。
本发明的有益效果可以通过以下仿真进一步说明。
仿真试验1,该试验对本发明中小幅面地形图进行分割的仿真。
仿真条件:在MATLAB7.0软件下进行。
参照图2,分别使用本发明方法与基于地形图空间信息的图像分割方法、基于遗传模糊C均值聚类方法、基于改进的模糊C均值聚类方法等三种传统的地形图分割方法对大小为795×680像素的小幅面地形图进行仿真实验。图2(a)是原始小幅面地形图图像,主要使用四种颜色表示不同的地理要素信息,但是同时也存在多种渐变色,过渡色和覆盖色。
从图2(d2,e2,e4)可以看出,传统图像分割算法得到的结果图像中存在大量的本不应该属于此版图的杂色。特别是在图2(e4)中,蓝色分版图像中含有大量的黑色像素,导致在黑色分版图像中的线划要素出现较为严重的断裂的情况,而与之相对应的由本发明方法得到的分割结果则具有较好的效果。在图2(c1’,c2’,c3’,c4’),图2(d1’,d2’,d3’,d4’)和图2(e1’,e2’,e3’,e4’)中,这些图像时分割结果的细节图像,从这些图像中可以看出使用本发明方法得到的分割结果与其它方法得到的分割结果相比,结果图像中残留的噪声更少,线划要素更完整,更连续。
仿真试验2,该试验对本发明中大幅面地形图进行分割的仿真。
仿真条件:在MATLAB7.0软件下进行。
参照图3,对大小为5000×5000像素大幅面地形图进行仿真实验。从图3(d)中可以看出,本发明方法分割得到的棕色版图中基本不含其他杂色,且等高线信息完整,线划连续;从图3(f)中可以看出,本发明方法分割得到的绿色版图中面状区域平滑连续,边界清晰。此外,本发明方法对大小为5000×5000的地形图进行分割需要的时间为325.7302(s),是传统方法所需时间1513.2350(s)的1/5。因此,本发明方法能够对大幅面地形图进行快速准确的分割。
分别对传统的图像分割方法和本发明方法计算分割效率评价参数,最终数据如表1所示。
表1不同大小地形图分割方法使用时间表
地形图大小 | 传统的图像分割方法 | 本发明方法 |
500×500 | 32.1287(s) | 13.8662(s) |
600×600 | 44.1108(s) | 19.6646(s) |
700×700 | 59.6015(s) | 26.5989(s) |
800×800 | 75.8195(s) | 20.1525(s) |
900×900 | 96.3398(s) | 26.0943(s) |
1000×1000 | 122.2289(s) | 31.3653(s) |
1100×1100 | 152.0572(s) | 26.3407(s) |
1200×1200 | 178.4432(s) | 31.2745(s) |
1300×1300 | 205.7123(s) | 36.5182(s) |
1400×1400 | 240.4408(s) | 42.1408(s) |
1500×1500 | 286.1583(s) | 49.0366(s) |
2000×2000 | 492.9253(s) | 73.1945(s) |
2500×2500 | 805.2184(s) | 112.6271(s) |
从表1的客观评价测量值可以看出,本发明在地形图分割效率上优于传统的地形图分割方法,能够对大幅面地形图进行快速分割。且随着图像大小的增加,本发明的高效性更明显。
Claims (8)
1.一种基于随机概率采样与多级融合的大幅面地形图分割方法,其特征在于,具体包括如下步骤:
步骤1,输入原始地形图图像;
步骤2,对原始地形图图像进行随机概率采样,得到采样数据矩阵;
步骤3,确定地形图分割种类数;
步骤4,计算聚类中心,得到聚类中心矩阵;
步骤5,计算原始地形图的隶属矩阵,得到隶属矩阵;
步骤6,利用隶属矩阵对原始地形图进行模糊分类,得到分割后的图像;
步骤7,对分割后的图像进行多级融合,得到的融合后图像;
步骤8,输出分割后的分版图像;
所述步骤2的具体步骤如下:
2a)生成一个符合正态分布的,且和原始地形图A等大小的随机矩阵RM×N作为随机概率矩阵,该矩阵中所有元素的值为0和1之间的任意随机数;其中,M为原始地形图的宽度,N为原始地形图的长度;
2b)计算概率阈值P,并利用步骤2a)生成的随机概率矩阵RM×N对原始地形图进行数据采样;概率阈值P的计算如下:
其中,M为原始地形图的宽度,N为原始地形图的长度;
数据采样过程如下:遍历随机概率矩阵RM×N,如果矩阵中某一元素的值大于等于概率阈值P,则原始地形图A中相同位置像素点的颜色信息被选择为采样信息;否则原始地图A中相同位置像素点的颜色信息为0;遍历结束后得到采样数据矩阵A'。
2.如权利要求1所述的方法,其特征在于,所述步骤3中,利用下式确定地形图的分割种类数:
Cluster_n=Cl 1+Cl 2+Cl 3+…Cl l+1
其中,Cluster_n为地形图的分割种类数,Cl 1为排列组合表达式,表示在由l种颜色绘制的地形图中的某区域有1种颜色相互覆盖所得到的颜色种类数,l≥1。
3.如权利要求1所述的方法,其特征在于,所述步骤4的具体步骤如下:
4a)初始化隶属矩阵:初始化一个与原始地形图A相同大小的随机矩阵作为初始隶属矩阵U0;并使其满足下式:
其中,u0,k(i,j)为初始隶属矩阵U0中坐标(i,j)位置第k类元素的值;
4b)利用下式计算得到矩阵大小为Cluster_n×3的聚类中心矩阵vI,其中,I为迭代次数,当迭代次数I=1时,隶属度矩阵U=U0,uk(i,j)=u0,k(i,j);
其中,vI,k为第I次迭代中聚类中心矩阵vI的第k个聚类中心,k=1,…,Cluster_n;m=2;x(i,j)为原始地形图中坐标(i,j)位置的像素值;(uk(i,j))m为隶属矩阵U中坐标(i,j)位置第k类元素值的m次方;
4c)根据下式计算第I次迭代的总体价值函数JI:
其中,JI,k为第I次迭代的第k类价值函数;
如果总体价值函数JI小于阀值Jmin=10-5,或总体价值函数JI相对上次得到的价值函数值JI-1的改变量小于阀值Jdif_min=10-5,则步骤4的计算聚类中心结束,最终得到聚类中心矩阵v=vI,然后执行步骤5;否则,进行步骤4d);
4d)使用下式计算新的隶属度矩阵U';
其中,u'k(i,j)为新的隶属矩阵U'中坐标(i,j)位置第k类元素的值;vI,k和vI,p分别表示聚类中心矩阵vI中第k个和第p个聚类中心;
令U=U',uk(i,j)=u'k(i,j),I=I+1,返回步骤4b)。
4.如权利要求1所述的方法,其特征在于,所述步骤5的具体步骤如下:
5a)读取步骤(4)中得到的聚类中心矩阵v;
5b)利用模糊C均值聚类算法中的隶属度计算公式,计算原始地形图A中各像素的隶属度,得到所有像素的隶属度组成隶属矩阵UF;
其中,uF,k(i,j)为隶属矩阵UF中位置为(i,j)处第k类元素的值,对应于原始地形图中坐标位置的像素属于k类的隶属度;vk和vp分别为聚类中心矩阵v的第k个和第p个聚类中心;x(i,j)为原始地形图中坐标(i,j)位置的像素值;m=2。
5.如权利要求1所述的方法,其特征在于,所述步骤6的具体步骤如下:
利用隶属矩阵UF,按以下公式所示的模糊分类方法对原始地形图A进行图像分割得到分割后的图像:
x(i,j)∈imagek如果uF,k(i,j)≥uF,p(i,j) (p=1,2,…Cluster_n)
其中,imagek为原始地形图中像素x(i,j)被划入的第k个分割后的图像。
6.如权利要求1所述的方法,其特征在于,所述步骤7的具体步骤如下:
7a)读取聚类中心矩阵;
7b)计算每个聚类中心和其他聚类中心之间的距离;
7c)利用聚类中心之间的距离以及多级融合方法对分割后的图像进行融合,直至融合后得到的分版图像数目和地形图的分割种类数Cluster_n相等时停止,得到融合后的分版图像。
7.如权利要求6所述的方法,其特征在于,所述步骤7b)中的距离采用马氏距离:
其中,Distk1,k2为聚类中心矩阵v中第k1个聚类中心vk1和第k2个聚类中心vk2之间的马氏距离,()T为矩阵的转置。
8.如权利要求6所述的方法,其特征在于,所述步骤7c)中的多级融合方法的规则:当且仅当与第i类的最小距离是第j类,同时与第j类的最小距离是第i类时,将第i类和第j类融合为一新类。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410029317.8A CN103824283B (zh) | 2014-01-22 | 2014-01-22 | 基于随机概率采样与多级融合的大幅面地形图分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410029317.8A CN103824283B (zh) | 2014-01-22 | 2014-01-22 | 基于随机概率采样与多级融合的大幅面地形图分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103824283A CN103824283A (zh) | 2014-05-28 |
CN103824283B true CN103824283B (zh) | 2016-08-17 |
Family
ID=50759325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410029317.8A Active CN103824283B (zh) | 2014-01-22 | 2014-01-22 | 基于随机概率采样与多级融合的大幅面地形图分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103824283B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105989594B (zh) * | 2015-02-12 | 2019-02-12 | 阿里巴巴集团控股有限公司 | 一种图像区域检测方法及装置 |
CN108428235B (zh) * | 2018-02-07 | 2021-08-27 | 西北大学 | 彩色地形图中线要素分割方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1932882A (zh) * | 2006-10-19 | 2007-03-21 | 上海交通大学 | 基于目标检测的红外与可见光序列图像特征级融合方法 |
CN102800094A (zh) * | 2012-07-13 | 2012-11-28 | 南京邮电大学 | 一种快速彩色图像分割方法 |
-
2014
- 2014-01-22 CN CN201410029317.8A patent/CN103824283B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1932882A (zh) * | 2006-10-19 | 2007-03-21 | 上海交通大学 | 基于目标检测的红外与可见光序列图像特征级融合方法 |
CN102800094A (zh) * | 2012-07-13 | 2012-11-28 | 南京邮电大学 | 一种快速彩色图像分割方法 |
Non-Patent Citations (5)
Title |
---|
Texture image segmentation using combined features from spatial and spectral distribution;K.Muneeswaran 等;《Pattern Recognition Letters》;20060531;第27卷(第7期);第755-764页 * |
基于随机抽样和聚类特征的聚类算法;周兵 等;《西安交通大学学报》;20031231;第37卷(第12期);第1234-1237页 * |
基于随机抽样的加速K-均值聚类方法;王秀华;《计算机与现代化》;20031218(第220期);第27-29、33页 * |
基于颜色特征的地图要素分割和识别;朱文忠;《模式识别与人工智能》;19960630;第9卷(第2期);第194-200页 * |
改进的快速模糊C均值聚类的图像分割方法;李明 等;《兰州理工大学学报》;20070630;第33卷(第3期);第95-99页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103824283A (zh) | 2014-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103208001B (zh) | 结合形状自适应邻域和纹理特征提取的遥感图像处理方法 | |
CN107451607B (zh) | 一种基于深度学习的典型人物的身份识别方法 | |
CN106056155B (zh) | 基于边界信息融合的超像素分割方法 | |
CN105138970B (zh) | 基于空间信息的极化sar图像分类方法 | |
CN104102919B (zh) | 一种有效防止卷积神经网络过拟合的图像分类方法 | |
CN105335966B (zh) | 基于局域同质性指标的多尺度遥感影像分割方法 | |
CN106096561A (zh) | 基于图像块深度学习特征的红外行人检测方法 | |
CN106296638A (zh) | 显著性信息取得装置以及显著性信息取得方法 | |
CN106203430A (zh) | 一种基于前景聚集度和背景先验的显著性物体检测方法 | |
CN102982338B (zh) | 基于谱聚类的极化sar图像分类方法 | |
CN105893925A (zh) | 基于肤色的人手检测方法及装置 | |
CN103325117B (zh) | 一种基于matlab的岩心图像处理方法及系统 | |
CN104282008B (zh) | 对图像进行纹理分割的方法和装置 | |
CN103679191B (zh) | 基于静态图片的自动套牌车检测方法 | |
CN107229917A (zh) | 一种基于迭代聚类的多幅遥感影像共性显著目标检测方法 | |
CN107564022A (zh) | 基于贝叶斯融合的视频显著性检测方法 | |
CN105427309A (zh) | 面向对象高空间分辨率遥感信息提取的多尺度分层处理方法 | |
CN104217440B (zh) | 一种从遥感图像中提取建成区的方法 | |
Casanova et al. | Texture analysis using fractal descriptors estimated by the mutual interference of color channels | |
CN103955926A (zh) | 基于Semi-NMF的遥感图像变化检测方法 | |
CN106683102A (zh) | 基于脊波滤波器和卷积结构模型的sar图像分割方法 | |
CN106611423A (zh) | 基于脊波滤波器和反卷积结构模型的sar图像分割方法 | |
CN102073867B (zh) | 一种遥感图像分类方法及装置 | |
CN107403434A (zh) | 基于两阶段聚类的sar图像语义分割方法 | |
CN109447111A (zh) | 一种基于子类训练样本的遥感监督分类方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210729 Address after: 400031 unit 1, building 1, phase 3, R & D building, Xiyong micro power park, Shapingba District, Chongqing Patentee after: Chongqing Institute of integrated circuit innovation Xi'an University of Electronic Science and technology Address before: 710071 No. 2 Taibai South Road, Shaanxi, Xi'an Patentee before: XIDIAN University |