基于栅格数据提取主骨架线的新算法.docVIP

文库下载 免费文档下载  HYPERLINK / / 　　 本文档下载自 HYPERLINK / 文库下载网，内容可能不完整，您可以点击以下网址继续阅读或下载：  HYPERLINK /doc/abb41ec5da714f.html /doc/abb41ec5da714f.html 基于栅格数据提取主骨架线的新算法 第25卷第5期2000年10月武汉测绘科技大学学报 　　ofSurveyingandMappingJournalofWuhanTechnicalUniversityVol.25No.5 Oct.2000 文章编号:1000-050X(2000)05-0432-36 基于栅格数据提取主骨架线的新算法 杜世宏1　杜道生1　樊　红1　万幼川2 (1　武汉测绘科技大学测绘遥感信息工程国家重点实验室,武汉市珞喻路129号,430079) (2　武汉测绘科技大学信息工程学院,武汉市珞喻路129号,430079) 3 摘　要　详细介绍了基于栅格数据的面状地物主骨架线的提取算法。该算法具有较快的处理速度和较强的适应性及内存需求少的优点,同时用该算法所提取的主骨架线为最优(最长)骨架线,可用于面状地物注记的自动配置。 关键词　主骨架线;数学形态学;栅格数据处理分类号　P283.7;P283.1　文献标识码　A 　　形状特征,,双线河流的简化。栅格骨架线还可以用于栅格数据的压缩。此外,骨架线还在模式识别和计算机视觉中得到应用。骨架线的提取算法通常有基于矢量数据和栅格数据两种方法。本文研究了基于节点的栅格数据骨架线的矢量化和基于端点的栅格数据的主骨架线提取算法。该算法先把面状地物栅格化,经过数学形态学的序贯减薄运算,形成地物的骨架线,然后剪去骨架线中的短支线,从而提取主骨架线。该算法所提取的主骨架线是最优(最长)的。应该指出,形态学减薄运算的时间效率受栅格数据数量的直接制约,栅格太小时,栅格数据量大,影响运算速度;栅格太大时,不易于提取小的面状地物的骨架线。本文提出基于固定栅格数目的动态栅格数据结构,能很好地解决上述矛盾。 行数和列数。这样的结构有利于对不同的多边形采取不同的精度(栅格大小),处理完后,可以方便地转换到统一的大地坐标系中去。因此,该结构既适用于整幅图形,也适用于??个多边形。栅格图像结构如下:　　typedef 　　 structtagIMAGE 33value;栅格图像灰度值:///doc/abb41ec5da714f.htmlar starow;栅格图像值的起始列数stacol;栅格图像值的起始行数 　　intcol;栅格图像值的总行数 　　introw;栅格图像值的总列数　　floatOrgPointx,OrgPointy;栅格图像的原点大地坐标 　　floatm-Sizex,m-Sizey;栅格图像的栅格大小 IMAGE,3PIMAGE; 1　栅格数据动态处理 提取面状地物的骨架线时,首先要将面状地 物栅格化。基于固定栅格数目的动态栅格数据,对于面状地物栅格化时,栅格的数目确定,然后按照数目自动确定栅格大小。这样既避免了处理大面状地物时因栅格数目太多而费时,又防止了因小面状地物栅格太大而不能精确提取骨架线。同时,该结构中保留了栅格图像的原点和栅格大小、 　收稿日期:2000-04-07. 3国家测绘局测绘科技发展基金资助项目(95014)。 2　骨架线提取算法 2.1　骨架线提取 骨架线的提取主要采用数学形态学序贯减薄 的方法,即 X○T=X/(X??T) 式中,X为栅格图像;结构元T采用结构元表中的L模板[1]。经过序贯减薄后的栅格图像为面状地物的骨架线。序贯减薄算法具有不改变原图 第5期　　　　　　　　　　　　杜世宏等:基于栅格数据提取主骨架线的新算法　　　　　　　　　　　　　　433 像各分量的拓扑邻接关系,并能保证提取的骨架线连续和单像元宽度等优点。2.2　圆环的消除及节点和端点的提取 表1　直线段的拓扑关系 Tab.1　TheTopologicalRelationofLines直线段序号起始栅格号终止栅格号 123 (0,2)(2,3)(5,8) (0,8)(4,6)(9,6) 坐标对 在栅格数据中,点的类型有断点、端点/doc/abb41ec5da714f.html、3节点、4节点。如图1所示,P为要考察的栅格,Pi为其八邻域的栅格。若Pi≥1,则Ni=1;若Pi=0,则Ni=0。当Ni=1时,Pi为目标栅格;Ni=0时,Pi为背景栅格。为方便叙述,各种点的类型定义为:当目标栅格P0时,如果∑Ni=0,则P为断点;如果∑Ni=1,则P为端点;如果∑Ni=3,则P为3节点;如果∑Ni=4,则P为4节点。本文