////////////////////////////////////////////////////////////////
//计算点是否在多边形内
int msIntersectPointPolygon(pointObj *point, shapeObj *poly)
{//主函数1-在 0-不在
int i;
int status=MS_FALSE;
for(i=0; i<poly->numlines; i++)
{
if(msPointInPolygon(point, &poly->line[i]) == MS_TRUE) /* ok, the point is in a line */
status = !status;
}
return(status);
}
int msPointInPolygon(pointObj *p, lineObj *c)
{//算法非常奇怪!!!
int i, j, status = MS_FALSE;
for (i = 0, j = c->numpoints-1; i < c->numpoints; j = i++)
{ if ((((c->point[i].y<=p->y) && (p->y<c->point[j].y)) || ((c->point[j].y<=p->y) && (p->y<c->point[i].y))) && (p->x < (c->point[j].x - c->point[i].x) * (p->y - c->point[i].y) /(c->point[j].y - c->point[i].y) + c->point[i].x))
status = !status;
}
return status;
}