OpenGL研究3.0 多边形区域填充
DionysosLai(906391500@qq.com)2014-06-22
所谓多边形区域填充。就是将多边形内部区域，所有已相同色块填充。注意：这里讨论的多边形是简单多边形（即不考虑诸如五角星这样的相交多边形）。简单多边形，分为凹多边形和凸多边形。
多边形区域填充有下面几种方法：
1. 逐点扫描方法：
原理：扫描多边形区域，逐点推断点是否在多边形内。
难点：在于怎样推断点是否在区域内;
经常使用怎样推断点是否在区域内方法：射线法、面积法。
面积法原理:取一个点。连接多边形各个点，依据三个点形成一个三角形原理，我们能够求得三角形面积。推断面积的大小，就能够推断该店是否在多边形内了。
射线法：这种方法。是我们这里要重点解说的一个方法。原理：取一个点。向左或者向右做一条射线过去，推断射线与多边形的交点。依据多边形交点熟练和本身多边形边情况，推断点是否在多边形内。
首先，射线从左向右，最左边点肯定在多边形区域为（我们这里假定射线方向水平向左）。那么与多边形相交第一个点，必定表明射线右部分在多变形内，与多边形相交第二个点。表明射线右部分在多边形外面。因此。通过推断射线与多变的交点奇偶性。推断点是否在多边形内。
这里。有几种特殊情况，例如以下图所看到的：
图a，射线与多变形顶点相交，顶点算一个；图b。射线与多边形顶点的交点。不被计算在内（注意图a和图b的差别->顶点纵坐标大小差别）；图c和图d，射线与多边形一条边重合，这条边被忽略不计。
因此。我们能够设计例如以下：取点向左做一条射线，1. 对于水平边不做考虑；2. 对于多边形顶点与射线交点情况。假设其纵坐标是所属边较大顶点，则计数（參考图a），否则不计数（图b）。3.对于点在多边形上情况，可直接推断点在多边形内。
伪代码例如以下：
count ← 0; 以P为端点。作从右向左的射线L; for 多边形的每条边s do if P在边s上 then return true; if s不是水平的 then if s的一个端点在L上 if 该端点是s两端点中纵坐标较大的端点 then count ← count+1 else if s和L相交 then count ← count+1; if count mod 2 = 1 then return true; else return false;
对应代码例如以下：注：我是在coco2dx2.3版本号内測试的。因此可能移植要改些类名。
///@brief 推断点是否在多边形
///@param[in] p0--要推断点， poly--多边形点集合， numberOfPoints--多边形点数量
///@return 2---点在多边形内， 1---点在多边边上，0---点不在多边形内
///@author DionysosLai,906391500@qq.com
///@retval
///@post
///@version 1.0
///@data 2014-04-11
int HelloWorld::pointIsInPolygon(const CCPoint& p0, const CCPoint* poly, const unsigned int numberOfPoints)
{
unsigned int count = 0; ///< 用来标记射线L与多边形的交点数；
CCSize winsize = CCDirector::sharedDirector()->getWinSize();
/// 已点p0向左向右做一条射线L；
CCPoint leftPoint = ccp(-100.0f, p0.y);
CCPoint rightPoint = ccp(winsize.width+100.0f, p0.y);
/// 推断每条边
for (unsigned int i = 0; i < numberOfPoints-1; i++)
{
/// 先推断点p0是否在边s上；
if (pointIsAtLine(p0, poly[i], poly[(i+1)%(numberOfPoints)]))
{
CCLOG("Point is at the %dth line", i);
return 1;
}
/// 推断边s是否是平行线；
if (poly[i].y != poly[(i+1)%(numberOfPoints)].y)
{
do
{
/// 推断边s的是否有端点在L上 同一时候 再推断该点是否是边s纵坐标较大的一个点
if (pointIsAtLine(poly[i], leftPoint, rightPoint))
{
if (poly[i].y > poly[(i+1)%(numberOfPoints)].y)
{
count += 1;
}
break;
}
if (pointIsAtLine(poly[(i+1)%(numberOfPoints)], leftPoint, rightPoint))
{
if (poly[i].y < poly[(i+1)%(numberOfPoints)].y)
{
count += 1;
}
break;
}
/// 假设边s没有端点在L上，则推断s与L是否相交
if (segmentLineIsIntersect(leftPoint, rightPoint, poly[i], poly[(i+1)%(numberOfPoints)]))
{
count += 1;
}
} while (0);
}
}
if (1 == count%2)
{
CCLOG("Point is not in polygon!");
return 0;
}
else
{
CCLOG("Point is in polygon!");
return 2;
}
}这里有个pointIsAtLine。是用来推断点是否在边上函数；segmentLineIsIntersect。是用来推断两条线段是否相交函数。可參考我的还有一边博文：http://blog.csdn.net/dionysos_lai/article/details/24418697计算几何文档一系列文章（眼下仅仅写了一篇，实际上是几乎相同写了经常使用几何算法，还没写成博文。怪楼主太懒了。）
Ok，逐点扫描推断方法就是差不都这样了。
2. 扫描线算法
逐点扫描算法，没有充分考虑到像素之间的连贯性。效率低。
扫描线算法。就是要利用像素之间的连贯性，提高算法效率。
所谓连贯性：有三个概念，1.边的连贯性，AB边与扫描线1相交，也可能与扫描线2相交；2.扫描线连贯性：当前扫描线与多边形边交点顺序，可能与下一条扫描线交点情况一致或者类似；3.区间连贯性：同一区间像素取同一颜色属性。
扫描线原理：将整个多边形区域扫描问题分解到一条条扫描线问题。仅仅要完毕每条扫描线的绘制，就实现了多边形区域填充问题。
一条扫描线与多边形有偶数个交点（0就不算了），按顺序每2个点形成一个区间，仅仅要绘制这个区间就可以。
难点这与扫描线与多边形边交点推断，这个是高中问题了，通过线段一般方程ax+by+c=0,两立方程求解。
只是这样的方法。要计算各种參数，比較费时，更好的方法是分成各种情况分开讨论（尽管比較麻烦）。能够关注我的《计算几何算法》系类文章。