基于点与多边形关系的地形图图幅查询系统实现.pdfVIP

68 测 绘 通 报 2008年 第1O期 文章编号：0494~911(2008)10．0068．02 中图分类号：P284 文献标识码：B 基于点与多边形关系的地形图图幅查询系统实现 张鹏程 (广州市城市规划勘测设计研究院，广东广州510060) RealizationofTopographicMapQueryingSystem Basedon RelationshipbetweenPointandPolygon ZHANG Peng—cheng 摘要：利用点与多边形关系判定中改进的射线算法，实现区域多边形范围内地形图图幅号的自动查询和输出功能，同其他方法相 比，该方法具有易于理解、稳定可靠 、通用性强的特点和优势。 关键词：点；多边形；地形图查询系统；地形图分幅编号 2．不能处理区域多边形被某一图幅包含的情 一 、 前 言 形。如图1中的多边形HIJKL，若以Ⅲ斜率与 日点 地形图图幅查询是城市基础地理信息系统最 和图幅24-38—10各点 (点 1，2，3，4)连线的斜率比 常用和基本的功能之一。文献 [1]在分析现有地形 较，则图幅24．38．10不满足条件，但若以 斜率与， 图图幅查询系统存在 “依赖第三方开发平台”和 “仅 点和图幅24．38．10各点连线的斜率 比较，则24—38． 提供查询区域最大外接矩形的区域查询”两大缺陷 10为满足条件的图幅(／J的斜率介于 ，1和 ，2之 的基础上，提出了一种基于斜率内插的算法，该文 间)，如此存在逻辑上的矛盾。 作者利用通用开发语言VB6．0，通过读取DXF格式 事实上，笔者认为利用空间数据拓扑关系中常 的区域边界，验证了该算法，该算法不仅实现了地 用的点与多边形关系模型，完全可以实现批量查询 形图上指定区域各种 比例尺的地形图图幅编号，而 且仅输出该区域所属地形图图幅编号的功能。 且能批量查询且仅输出该区域所属地形图图幅编 D 号，算法设计巧妙，计算简便、效率高，具有较高的 N4 ／ 推广应用价值。 。 笔者认为文献 [1]仍然存在如下两个方面的 问题 ： 十e 工ffJ 1．忽略了斜率不存在的情形。区域多边形边 n—、 f Ⅻ 、 一 的斜率与多边形边界点和图廓点连线斜率的大小 关系是文献 [1]算法的核心和关键，即计算斜率的 图1 斜率不存在及图幅被多边形包含示意 大小是前提，但作者忽略了当斜率不存在时这一 重要情形。斜率不存在表现在两个方面：一是区 域多边形存在平行于 y轴的边；二是多边形边界 二、点与多边形的关系及算法 点与图廓点的连线平行于 】，轴。由于大多数区域 点与多边形的关系包括点在多边形内、外两种