GIS面裁减
李泉 陈玉进(南京跬步科技有限公司http://www.creable.cn )
裁剪是从数据集合中抽取所需信息的过程,它是地理信息系统中许多重要问题的基础。裁剪最典型的用途是开窗裁剪,即将当前窗口显示范围之外的部分切除,以至于不绘制或者不标注不需要的部分。除了开窗裁剪之外,使用硬网格索引在处理数据阶段也需要将地图数据按照网格裁切为一个个的块状区域。叠加分析需要用的任意多边形和任意多边形的裁剪。本文的目的旨在详细介绍Weiler-Atherton算法和这种算法在实际应用中容易出现的问题。
一、Weiler-Atherton算法简介
该算法能够实现凹多边形与凹多边形的裁剪,可以用一个含内孔的凹多边形去裁剪另一个也有内孔的凹多边形,且效率比较高,但该算法略显复杂。
算法将主多边形和裁剪多边形定义为顶点的环形列表,多边形的外边界取顺时针方向,内边界或孔取逆时针方向。这种约定保证了在遍历顶点列表的时候,多边形的内部总是位于前进方向的右侧。主多边形和裁剪多边形的边界可能相交,也可能不相交。若他们相交,则交点必成对地出现。其中一个交点为主多边形边进入裁剪多边形内部时的交点(通常称作入点),而另一交点为其离开时的交点(通常称作出点)。算法从入点开始,沿主多边形的外部边界按照顺时针方向向前跟踪,直至找到一个出点为止。在出点处跳跃到裁剪多边形节点序列上,沿裁剪多边形的外部边界按照顺时针方向跟踪,直至发现一个入点,再次沿主多边形的边界跟踪。继续上述过程,直至到达此次的起始点为止,并将此次经过的所有入点都标识为已读。主多边形的内边界则按照逆时针方向跟踪。再次从头查找没有读过的入点,如没有入点