Sistemas de Informacion Geografico

[Mefistofeles]

Hola gente, alguien sabe de alguna funcion que verifique si un punto pertenece a un poligono? Estoy trabajando con SIG (Sistemas de Informacion Geografico), utilizo una tabla con un campo tipo POLYGON y quiero saber si un valor tipo POINT pertece a ese campo POLYGON. Espero que alguien este trabajando con estos tipos de datos. Desde ya, muchas gracias.

[JCarlosas]

Hola, de la gaveta empolvado saque esto algoritmos que alguna vez use.
Ahora no recuerdo la diferencia entre ellos, pero creo que radica en que si el ulitmo punto tiene o no que ser igual al primero.
En su epoca me funcionaron bien.
Espero que resuelvas.
Saludos
Juan Carlos

Código Delphi [-]
Function EstaElPuntoDentroDelPoligono(P:TPointDouble;Pol:Array of TPointDouble):Boolean;
Var
 Nc, //Número de cortes
 I : Integer;
 Alfa,
 XQ : Double;      //Era antes Real
 Db, Da: TCoordenada;   //Eran antes entero
Begin
  Result := False;
  Nc := -1;
  for I := 0 to length(Pol)-1 do
    begin
      If I = Length(Pol)-1 then
        Begin
          Da := Pol[i].Y - P.Y;
          Db := Pol[0].Y - P.Y;
          If ((Da >= 0) And (Db < 0)) Or ((Db >= 0) And (Da < 0)) then
            Begin
              Alfa := Da/(Da-Db);
              XQ := Pol[i].X + Alfa * (Pol[0].X-Pol[i].X);
              If XQ > P.X Then Inc(Nc);
            End;
        End
      Else
        Begin
          Da := Pol[i].Y - P.Y;  //dA = yA - yP y dB = yB - yP.
          Db := Pol[I+1].Y - P.Y;
          If ((Da >= 0) And (Db < 0)) Or ((Db >= 0) And (Da < 0)) then
            Begin
              Alfa := Da/(Da-Db);
              XQ := Pol[i].X + Alfa * (Pol[I+1].X-Pol[i].X); //xQ = xA + alfa(xB-xA).
              If XQ > P.X Then Inc(Nc);
            End;
        End;
    end;
  If (Not Odd(Nc)) And (Nc >= 0) Then Result := true;
End;

Function EstaElPuntoDentroDelPoligono2(P:TPointDouble;Pol:Array of TPointDouble):Boolean;
Var
 Nc, //Número de cortes
 I, LongPol : Integer;
 Alfa,
 XQ : Double;      //Era antes Real
 Db, Da: TCoordenada;   //Eran antes entero
Begin
  Result := False;
  Nc := -1;
  for I := 0 to length(Pol)-2 do
        Begin
          Da := Pol[i].Y - P.Y;  //dA = yA - yP y dB = yB - yP.
          Db := Pol[I+1].Y - P.Y;
          If ((Da >= 0) And (Db < 0)) Or ((Db >= 0) And (Da < 0)) then
            Begin
              Alfa := Da/(Da-Db);
              XQ := Pol[i].X + Alfa * (Pol[I+1].X-Pol[i].X); //xQ = xA + alfa(xB-xA).
              If XQ > P.X Then Inc(Nc);
            End;
        End;
  If (Not Odd(Nc)) And (Nc >= 0) Then Result := true;
End;
Function EstaElPuntoDentroDelPoligono3(P:TPointDouble;Pol:Array of TPointDouble):Boolean;
Var
 Nc, //Número de cortes
 I : Integer;
 Alfa,
 XQ : Double;      //Era antes Real
 Db, Da: TCoordenada;   //Eran antes entero
Begin
  Result := False;
  Nc := -1;
  for I := 0 to length(Pol)-1 do
    begin
      If I = Length(Pol)-1 then
        Begin
          Da := Pol[i].Y - P.Y;
          Db := Pol[0].Y - P.Y;
          If ((Da >= 0) And (Db < 0)) Or ((Db >= 0) And (Da < 0)) then
            Begin
              Alfa := Da/(Da-Db);
              XQ := Pol[i].X + Alfa * (Pol[0].X-Pol[i].X);
              If XQ > P.X Then Inc(Nc);
            End;
        End
      Else
        Begin
          Da := Pol[i].Y - P.Y;  //dA = yA - yP y dB = yB - yP.
          Db := Pol[I+1].Y - P.Y;
          If ((Da >= 0) And (Db < 0)) Or ((Db >= 0) And (Da < 0)) then
            Begin
              Alfa := Da/(Da-Db);
              XQ := Pol[i].X + Alfa * (Pol[I+1].X-Pol[i].X); //xQ = xA + alfa(xB-xA).
              If XQ > P.X Then Inc(Nc);
            End;
        End;
    end;
  If Not (Odd(Nc)) or (Nc < 0) Then Result := true;
End;