Geometry in Pascal/de

//  Die Funktion gibt den Wert True zurück, wenn der Punkt x,y innerhalb eines Polygons liegt.
//  Die Funktion gibt den Wert False zurück, wenn der Punkt x,y nicht in einem Polygon liegt.
//
//  Der Original C-Code: http://www.visibone.com/inpoly/inpoly.c.txt
//
//  Übersetzng von C nach Pascal: Felipe Monteiro de Carvalho
//
//  Lizenz: Public Domain

function IsPointInPolygon(AX, AY: Integer; APolygon: array of TPoint): Boolean;
var
  xnew, ynew,
  xold, yold,
  x1, y1,
  x2, y2:     Cardinal;
  i, npoints: Integer;
  inside:     Integer = 0;

begin
  Result := False;

  npoints := Length(APolygon);

  if (npoints < 3) then Exit;

  xold := APolygon[npoints-1].X;
  yold := APolygon[npoints-1].Y;

  for i := 0 to npoints - 1 do begin
    xnew := APolygon[i].X;
    ynew := APolygon[i].Y;

    if (xnew > xold) then begin
      x1 := xold;
      x2 := xnew;
      y1 := yold;
      y2 := ynew;
    end else begin
      x1 := xnew;
      x2 := xold;
      y1 := ynew;
      y2 := yold;
    end;

    if (((xnew < AX) = (AX <= xold)) and    // Die Ecke ist am linken Ende "offen"
       ((AY - y1) * (x2 - x1) < (y2 - y1) * (AX - x1))) then
    begin
      inside := not inside;
    end;

    xold := xnew;
    yold := ynew;

  end;

  Result := inside <> 0;
end;