1、原理
从目标点引一条射线,判断这条线与曲线的交点个数为奇数还是偶数,为奇数表示点在曲线外,偶数表示点在曲线内。实际计算的时候,只需要判断一侧即可。
C 代码:
//xFlat为目标点,points为线段中的点
int nInsideFlag = 0;
for (int pIdx = 0; pIdx < nLines; ++pIdx)
{
points->GetPoint(pIdx, p1);
points->GetPoint(pIdx + 1, p2);
float fSlope = 0.f;
if (fabs(p2[0] - p1[0]) > VTK_SIZEOF_FLOAT)
{
fSlope = (p2[1] - p1[1]) / (p2[0] - p1[0]);
}
int nBetweenFlag = (p1[0] < xFlat[0] && xFlat[0] < p2[0]) || (p2[0] < xFlat[0] && xFlat[0] < p1[0]);
int nAboveFlag = xFlat[1] < (fSlope*(xFlat[0] - p1[0]) + p1[1]);
if (nBetweenFlag&&nAboveFlag)
{
nInsideFlag = !nInsideFlag;
}
}
if (nInsideFlag)
{
sign = -1;
}
else
{
sign = 1;
}
Js代码:
function pointInPolygon(point, lines) {
var inside = false;
for (var i=0; i<lines.length-1; i++) {
var slope = (lines[i+1].y - lines[i].y)/(lines[i+1].x - lines[i].x);
var xBetween = (lines[i].x < point.x && point.x <lines[i+1].x)
|| (lines[i+1].x < point.x && point.x <lines[i].x);
var yAbove = point.y < (slope * (point.x - lines[i].x) + lines[i].y);
if (xBetween && yAbove) {
inside = !inside;
}
}
return inside;
}