1.一种专题地图地理面要素自动提取和矢量化方法,其步骤包括:
1)获取专题地图并判断每一专题地图是否包含地理投影和空间参考信息,若没有则对该专题地图进行地理配准;
2)对步骤1)处理后的每一专题地图进行自适应图像聚类,并对所得聚类图像进行中值滤波以去除噪声点;
3)对于每一聚类图像,获取该聚类图像中各要素的边界并扩展边界像素点的范围;创建一个与该聚类图像相同大小的新图像,遍历该聚类图像中每一边缘像素点,并对其周围像素进行灰度值统计,将出现频率最高的灰度值赋值给该新图像同位置的中心像素点;遍历结束后,对于新图像未赋值的像素点,其像素值与该聚类图像同位置的像素相同;然后对新图像进行中值滤波;
4)根据专题地图a对应的新图像a’的像素灰度值生成各类地理要素的二值图,然后根据每一类二值图的像素点的灰度值和空间位置的连续性生成表示面要素的多边形边界,并根据专题地图a的地理投影和空间参考信息为每个多边形对象创建相应的投影和空间参考,最后输出专题地图a的各类面要素的矢量文件。
2.如权利要求1所述的方法,其特征在于,对专题地图进行自适应图像聚类的方法为:
21)输入待处理专题图像img;
22)对专题图像img进行灰度转换,获得其灰度图像imgg;
23)统计灰度图像imgg的灰度特征,得到对应的灰度分布图;
24)计算图像imgg的步长s;
25)计算图像imgg中的类别个数k;
26)对图像img进行k-means聚类;
27)对每一类像素重新赋值,不同类别的像素赋予不同的灰度值,形成聚类图像imgk;
28)对聚类图像imgk进行中值滤波,去除噪声点,输出最终灰度图像imgc。
3.如权利要求2所述的方法,其特征在于,步骤3)的具体实现方法为:
31)对灰度图像imgc进行边缘检测,生成二值图像imgcanny;
32)对imgcanny进行图像膨胀操作,扩展边缘像素的范围,得到新的二值图像imgdilate;
33)获取二值图像imgdilate中各边缘像素,对imgc中相同位置的像素重新赋值,从而消除图像imgc中的干扰要素,然后进行中值滤波,获得图像imgl。
4.如权利要求3所述的方法,其特征在于,获得图像imgl的方法为:
331)获取imgdilate中表示边缘要素的像素索引,并建立imgc中同索引的像素点集合:edge={p(m,n)|dn(pi(m,n))=1,pi(m,n)∈imgdilate,p(m,n)∈imgc};其中,pi(m,n)为二值图imgdilate中索引为(m,n)的像素点,p(m,n)为灰度图imgc中索引为(m,n)的像素点,dn(pi(m,n))表示图像imgdilate中像素点pi(m,n)的灰度值;
332)创建灰度图imgn,即imgn=imgc;
333)对灰度图imgc进行滑窗统计,并将统计结果赋值给灰度图imgn的同位置中心像素点;
334)对步骤333)处理后的灰度图imgn进行中值滤波去除噪声点,获得灰度图像imgl。
5.如权利要求3所述的方法,其特征在于,对图像img进行k-means聚类的方法为:首先将图像img各像素点的rgb值作为输入集合:{p0,p1,…,pn},pn的rgb值为(rn,gn,bn),其中n代表图像img像素点的总数;然后从输入集合{p0,p1,…,pn}中随机选取k个点作为初始聚类中心(c1,c2,…,ck),分别对应着k个类别(c1,c2,…,ck),ck的rgb值为(rk,gk,bk);然后计算各像素点到聚类中心点的距离,并将像素点归类至距离最短的类别内;然后更新聚类中心点
6.如权利要求2所述的方法,其特征在于,
7.如权利要求2所述的方法,其特征在于,计算图像imgg中的类别个数k的方法为:
251)初始化k=0;
252)根据灰度值及其对应的像素点个数,对每个灰度值进行峰值判断:
253)遍历所有的dn之后,最终得到的k值即为图像imgg的类别数。
8.如权利要求1所述的方法,其特征在于,利用canny边缘检测方法获取专题地图中所有要素的边界。
9.一种服务器,其特征在于,包括存储器和处理器,所述存储器存储计算机程序,所述计算机程序被配置为由所述处理器执行,所述计算机程序包括用于执行权利要求1至8任一所述方法中各步骤的指令。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至8任一所述方法的步骤。