点在多边形内的判断

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下测试通过，能正确判断点与多边形的位置关系。