空间数据拓扑关系的自动生成

简单面复杂面
图3
图4
多边形一般包括：多边形号、中心点坐标、多边形属性数据、多边形的组成弧段号、多边形岛的信息，考虑到组成弧段的方向和多边形顶点序列的方向存在的可能的不一致性以及效率问题，可以改为记录下组成多边形的弧段指针和方向性信息，即弧段与多边形的方向是否一致，对于岛的信息则通过将构成变形的边线分块来处理的方式体现，比如多边形包含岛屿，则可以使多边形的外边界成为多边形的第一部分，岛屿作为多边形的二、三、四等部分的方式加以解决。多边形的数据结构可以表示如下：
List<Arc*> Arcs；//弧段集合
List<Arc*> NewArcs；//处理后的弧段集合
void BreakArc(Arc *arc)
{
Arc * first=0；//打断后的第一部分
Arc * second=0；//打断后的第一部分
bool IsSelfCross=false；
int LineCountofArc=弧段所包含的直线段数目；
拓扑关系自动生成算法的一般过程为：
1．弧段处理，使整幅图形中的所有弧段，除在端点处相交外，没有其他交点，即没有相交或自相交的弧段。
2．结点匹配，建立结点、弧段关系。
3．建立多边形，以左转算法或右转算法跟踪，生成多边形，建立多边形与弧段的拓扑关系。
4．建立多边形与多边形的拓扑关系。
5．调整弧段的左右多边形标识号。
Node * _end; //终止结点指针
Polygon * _LeftPolygon; //弧段左多边形指针
Polygon * _RightPolygon; //弧段右多边形指针
public:
Arc() {…}; //构造函数
~Arc() {…}; //析构函数
OtherMethod…; //其他公共操作
由于GIS的数据量大，造成了判断的工作量大，效率低下的弊端，在判断两条弧段的关系时，应尽可能的减少计算量，减少计算量的工作可以分两步来做，首先是判断两条弧段的最小矩形壁包（Minimum bouding rectangle，MBR）是否相交或具有包含关系，如果不相交或没有包含关系，那么可以断定两条弧段是不互相交的，如果相交或具有包含关系，则进一步判断第一条弧段的每一条组成线段是否和第二条弧段的MBR相交或被包含，如果不相交或没有被包含则可以判断这一部分线段不会和第二条弧段相交，否则可以使用这一条线段与组成第二条弧段的各个线段进行相交关系的判定来确定交点。
}
private:
long _PolygonID;
Point _CenterPoint;
vector<Part *> _Parts;
public:
Polygon() {…};
~Polygon() {…};
OtherMethod() {…};
OtherProperty() {…};
}
四、拓扑结点、弧段和面之间的关系
for(int i=0；i<LineCountofArc-2；i++)
{
for(int j=i+2；j<LineCountofArc；j++)
{
if(直线段i和直线段j相交)
{
将交点Ｐ插如弧段，并在Ｐ处将弧段断为两截，分别存入first和second
中；
IsSelfCross=true；
Break；
6．多边形内部标识号的自动生成。
事实上，拓扑关系的生成过程中还涉及到许多工作，例如弧段两端角度的计算、悬挂结点和悬线的标识、多边形面积计算、点在多边性内外的判别等。
第一节拓扑关系的计算机表达
一、拓扑结点
结点用来描述如管线的交点、道路路口等现实世界的特征对象，结点可以用来检测弧段与弧段的连接关系和多边形特征是否能正确地完成。只与一条弧段相连接的起点或终点叫做悬挂结点。如图1所示P点就是悬挂结点：
L（悬挂弧段）
图2
图1
弧段一般包括：弧段号、弧段节点坐标串、弧段起始和终止结点、弧段左右多边形，弧段的数据结构可以表示如下：
class Arc
{
private:
long _ID; //弧段号
vector<Point > _Points; //弧段节点坐标串指针
Node * _start; //起始结点指针
弧段相交打断处理的算法描述如下：
//计算两直线段的交点情况
bool LineCross(Line * first, Line * second, Point & p1)
{
double dx1=first->end->x-first->head->x;
double dx2=second->end->x-second->head->x;
Point *head;
Point *end;
public:
OtherMethod();
}
//判断两个矩形是否相交或具有包含关系
bool RectangleCross(Rectangle one ,Rectangle two)
{
if(one.bottom>two.top || one.right<two.left || one.left>two.right || one.top >two.bottom)
图1
结点一般包括：结点号、结点坐标、与该结点连接的弧段集合，结点的数据结构可以表示如下：
class Node
{
private:
long _ID; //结点号
Point * _Point; //指向结点坐标的指针
vector<Arc<Point > *> ArcCollection ; //与该结点相联接的弧段集合指//针
double temp=dx1*dy2-dy1*dx2;
double u= (cx*dy2-cy*dx2)/temp;
double v=(cx*dy1-cy*dx1)/temp;
if(u>=0 && u<=1 && v>=0 && v<=1)
{
p1.x=first->head->x+dx1*u;
p1.y=first->head->y+dy1*u;
return true;
}
return false;
}
class Line
{
friend LineCross(Line * first,Line * second,Point & p1,Point & p2);
public:
Line(Point *,Point *);
~Line();
private:
class Polygon
{
public:
class Part
{
private:
typedef pair<Arc<Point> * ,bool > ArcElement;
vector<ArcElement> _Arcs;
public:
Part() {…};
~Part() {…};
OtherMethod() {…};
拓扑关系生成后，拓扑结点，拓扑弧段和拓扑面之间的关系如下几张表所示：
弧段号
结点号
A0
A1
.
An
N00N10
N10N11
.
NN0NN1
表1弧段——结点关系表
结点号
弧段号
Ni
AiAjAk……
表2结点——弧段关系表
弧段号
左多边形号
右多边形号
A0
A1
.
An
PL0
PL1
.
PLn
PR0
PR1
.
PRn
表它们的参数方程分别为：
和
判断两线段有无交点的关键变为判断 和 是否符合不等式 且 。令：
，
，
，
有： ，
如果 =0，说明两线段平行或者重合，没有交点，或者交点在两线段的头或尾上，否则如果满足不等式 且 ，两线段有交点，交点在两线段的中间。
二、自相交弧段处理
具有自相交特征的弧段至少具有四个（结）节点，由三个点或两个点组成的弧段不可能自相交。依次取出每一条弧段，如果弧段的（结）节点个数不少于四个，就利用直线段相交的方法，对组成弧段的各直线段进行判断，如果相交，将线段断开为两条，自相交的弧段可能不止有一处相交，可以通过递归的方法来将弧段分开，算法如下：
double dy1=first->end->y-first->head->y;
double dy2=second->end->y-second->head->y;
double cx=second->head->x-first->head->x;
double cy=second->head->y-first->head->y;
for(Position ite=Arcs.begin()；ite!=Arcs.end()；ite++)
{
Arc * arc =*ite；
BreakArc(arc)；
}
}
三、弧段相交打断处理
弧段与弧段相交关系的判断，可以通过取每一条弧段，与其他未判断过的所有弧段目标进行相交关系判断而得，从而要进行 次判断，具体方法为，取出第一条弧段，与其他 条弧段进行相交判断，求得交点后，将交点分别插入第一条弧段和与其相交弧段的对应位置上，并记录位置，将第一条弧段与所有其他弧段的相交关系判断完毕后，通过记录下的交点位置将第一条弧段分割，然后依次取出下一条弧段进行同样的处理，直到所有弧段处理完毕。