点在多边形内的判断

c语言代码：

typedef struct tagVERTEX
{
double x;
double y;
}VERTEX;

//判断点与多边形关系，0-外，1-上，2-内
int PointInPolygon(VERTEX &pt,VERTEX *vtx,int N)
{
int i,j,count=0;
double cp;

for(i=0,j=N-1;i<N;j=i++)//J为后端点，I为前端点
{
if(pt.y>=vtx[j].y&&pt.y<=vtx[i].y||pt.y>=vtx[i].y&&pt.y<=vtx[j].y)//点在线段y范围内
{
cp = (vtx[i].x-vtx[j].x)*(pt.y-vtx[j].y)-(pt.x-vtx[j].x)*(vtx[i].y-vtx[j].y);//向量JI与JP叉积
if(cp==0&&(pt.x>=vtx[j].x&&pt.x<=vtx[i].x||pt.x>=vtx[i].x&&pt.x<=vtx[j].x))//点在线段上
{
return 1;
}
if(cp>0&&vtx[i].y>vtx[j].y&&pt.y!=vtx[i].y||cp<0&&vtx[i].y<vtx[j].y&&pt.y!=vtx[j].y)//假定以pt为起点的向右水平射线与线段有交点
{
count++;
}
}
}

return (count%2 ? 2 : 0);
}

上述函数在vc++6.0下测试通过，能正确判断点与多边形的位置关系。