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

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