1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:

function PtInPolyDouble (const Points: array of TPoint; Count: integer; X, Y: double): Boolean;

 var bInside : boolean;
     i,j : integer;
     rkU0,rkU1 : ^TPoint;
     fRHS, fLHS : double;
begin
    bInside := false;
    j:= Count-1;
    for i := 0 to Count-1 do
    begin
        rkU0 := @Points[i];
        rkU1 := @Points[j];

        if ( y < rkU1.y ) then  // U1 above ray
        begin
            if ( rkU0.y <= y ) then  // U0 on or below ray
            begin
                fLHS := (y-rkU0.y)*(rkU1.x-rkU0.x);
                fRHS := (x-rkU0.x)*(rkU1.y-rkU0.y);
                if ( fLHS > fRHS ) then
                bInside := not bInside;
            end;
        end
        else if ( y < rkU0.y ) then  // U1 on or below ray, U0 above ray
        begin
            fLHS := (y-rkU0.y)*(rkU1.x-rkU0.x);
            fRHS := (x-rkU0.x)*(rkU1.y-rkU0.y);
            if ( fLHS < fRHS ) then
                bInside := not bInside;
        end;

     j:=i;
    end;
  Result:=bInside;
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:

function PointInPolygon(var xp, yp : array of Single; x, y: Single) : Boolean;
// The code below is from Wm. Randolph Franklin <wrf@ecse.rpi.edu>
// with some minor modifications for speed.  It returns 1 for strictly
// interior points, 0 for strictly exterior, and 0 or 1 for points on
// the boundary.
var
   I, J: Integer;
begin
   Result:=False;
   if High(XP)=High(YP) then begin
      J:=High(XP);
      for I:=0 to High(XP) do begin
         if ((((YP[I]<=Y) and (Y<YP[J])) or ((YP[J]<=Y) and (Y<YP[I])) )
             and (X<(XP[J]-XP[I])*(Y-YP[I])/(YP[J]-YP[I])+XP[I])) then
            Result:=not Result;
         J:=I;
      end;
   end;
end;