VC++开发GIS系统（284）判断面元素之间的拓扑关系

本文介绍GIS系统中面元素之间拓扑关系的判断方法。

关于面元素之间的拓扑

通常情况下可以将面元素之间拓扑关系做如下划分：
两个面完全重合
两个面相交
两个面相离
一个面内含另外一个面

其中为了简化，可以将1与4看作同一种情况。

面元素之间的拓扑算法思路

将一个面拆分成一条折线，并判断其是否与另一面相交，有一条相交则两个面相交。
如果一个面所有的点在另外一个面内，则包含
非上述情况则相离

算法实现

在 _malaBase.h 中声明如下的函数：

int  polygonAndPolygon(vector<malaPoint>&Polygon1, vector<malaPoint>&Polygon2);//面与面的拓扑关系 相交返回2，内含返回3,相离返回4

@ 

Polygon1

 面要素1
@ 

Polygon2

 面要素2

并在 _malaBase.cpp中实现上述函数，代码如下：

/*
* 面与面的拓扑关系  相等返回1，相交返回2，内含返回3,相离返回4，覆盖返回5
*/
int malaLogic::polygonAndPolygon(vector<malaPoint>&Polygon1, vector<malaPoint>&Polygon2)
{
if (polyLineAndPolygon(Polygon1, Polygon2) == 2)
return 2;
int length = Polygon1.size();

int index = 0;
int i;
for (i = 0; i < length; i++)
{
if (isPointInPolygon(Polygon1[i], Polygon2))
index++;
}
if (index == length)
return 3;
int length2 = Polygon2.size();
index = 0;

for (i = 0; i < length2; i++)
{
if (isPointInPolygon(Polygon2[i], Polygon1))
index++;
}
if (index == length2)
return 5;
return 4;
}

其中函数 

polyLineAndPolygon

 用于判断折线与面的拓扑关系，可以参考文章《VC++开发GIS系统（283）判断折线与面的拓扑关系》 

isPointInPolygon

 用于判断点与面的拓扑关系，可以参考文章《VC++开发GIS系统（280）判断点与面的拓扑关系》

详细代码可以参考：GitHub
VC++开发GIS系统（283）判断折线与面的拓扑关系VC++开发GIS系统（285）判断拓扑关系类的设计及实现