矢量地形图要素边线直角化方法研究.docx

112 测 绘 通 报 2014 年 第 12 期 引文格式: 应国伟，侯华斌，刘江，等． 矢量地形图要素边线直角化方法研究［J］． 测绘通报，2014 ( 12) : 112-113． DOI: 10． 13474 / j． cnki． 11-2246． 2014．0412 矢量地形图要素边线直角化方法研究 应国伟，侯华斌，刘 江，吴 思 ( 四川省第三测绘工程院，四川 成都 610500) Ｒight Angle Method of Edge Elements in Vector Topographic Map YING Guowei，HOU Huabin，LIU Jiang，WU Si 摘要: 针对矢量地形图数据中要素边线直角化问题，分析了最小二乘法的优缺点，提出了一种新的解决方案———边线中点固定算 法，该方法的优势在于算法模型简洁、运算速度快，适用于大数据量的处理和运算。本文详细描述了该解决方案的实现过程，并分 析了其存在的不足和适用范围。 关键词: ArcGIS; 地形图; 直角化; 最小二乘法 中图分类号: P217 文献标识码: B 文章编号: 0494-0911( 2014) 12-0112-02 一、引 言 数字地形图按其数据采集方法的不同可分为全 野外数字化测图、摄影测量数字化测图和扫描矢量化 测图［1］。但无论哪种方法都会无可避免地遇到一个 问题: 实际为矩形的要素，特别是建筑物要素，采集到 地形图上时并不是规则的矩形，而是接近矩形，本文 三、直角化的方法 伪矩形的存在使得地图表达的元素与实际存在 差异，因此，对伪矩形直角化是一项非常重要的工 作。解决边线的直角化问题必须考虑的是，一方面需 要将近似矩形直角化为标准的矩形; 同时，还要使每 n i个拐点的位移量尽量最小( 如图 2 所示) ，即∑s 2 为 i 称之为伪矩形。而如何在精度许可的范围内对这些 伪矩形进行直角化处理，是本文探讨的问题。 二、伪矩形产生的原因和表现 伪矩形产生的原因一般是由扫描纠正误差、测 量误差、立体测图误差等方面引起的。其表现形式 是 DLG 图形线划发生变形和扭曲，图形的内角与多 边形的实际直角存在差异［2］，如图 1 所示。 i = 1 最小。因此，笔者可以考虑按以下方法来达到需要。 图 2 伪矩形直角化 1． 最小二乘法算法 解决边线的直角化问题，不少学者已经作过相 关研究。但不难发现，大多数学者都不约而同地想 到了最小二乘法算法。根据最小二乘法原理，设点 在直角化之前的横坐标为 x，直角化之后的位置坐 标为 a，使拟合函数 y = f( x，a) 在 xi 处的值与试验 n 图 1 伪矩形产生原因 收稿日期: 2013-11-06 数值的偏差平方和最小，即∑［f( x ，a) － y ］2 取得 i ii = i i 最小值，从而求得最小二乘解 a = ( a1 ，a2 ，…，an ) 。 最小二乘法算法的优点是技术优化、结果可靠， 基金项目: 数字制图与国土信息应用工程国家测绘地理信息局重点实验室开放研究基金( GCwd201209) ; 测绘地理信息公益性行业科研 专项经费( 201412010) ; 四川省测绘地理信息局科技计划( J2013SZ06) 作者简介: 应国伟( 1972—) ，男，四川筠连人，高级工程师，研究方向为地理信息系统。 2014 年 第 12 期 应国伟，等: 矢量地形图要素边线直角化方法研究 113 缺点是当房屋边线比较多时，则增加递归运算负荷， 程序编制难度也较大; 而且最小二乘法只考虑图形 结果的最优符合性，并不考虑图形的实际真实性。 如图 3 所示，实线为矢量采集的建筑物边线，虚线为 经过最小二乘法优化的建筑物边线，把图 3 的 A 处 放大( 如图 4 所示) 可以看出，经优化后的边线已经 完全偏离了实际边线。 3) 获得初始边 L12 的初始方位角后，从 L12 中点 M12 起，以任意长度确定一条直线 M12 K。 4) 求第 2 条边 L23 中点 M23 到直线 M12 K 的垂足 点，此点即为节点 2 的改正点 2。 5) 以 2与第 2 条边 L23 中点 M23 确定直线，求第 3 条边 L34 中点 M34 到直线 M122的垂足点，此点即为 节点 3 的改正点 3。 6) 依次求出所有节点的改正坐标。 图 3 最小二乘法直角 化结果 图 4 最小二乘法直角 化结果( 放大图) 因此，笔者提出了另外一种新的算法，并称之为 边线中点固定算法。 2． 边线中点固定算法 相对最小二乘法算法是一种绝对算法来说，边 线中点固定算法是一种近似算法。如图 5 所示，设 虚