1

Buenas tardes, estoy con el siguiente problema. Implementé en Genexus usando Map Control el algoritmo para saber si un Punto dentro de un polìgono. Como estaba en c# lo traduje a genexus, este es el resultado:

parm(in:&GxMapPoint,in:&GxMapPolygon,out:&salida);

&salida='N'
&PuntoXY=&GxMapPoint.PointLat.Trim()+", "+&GxMapPoint.PointLong.Trim()
&y = &GeoLocationAPI1.GetLatitude(&PuntoXY)
&x= &GeoLocationAPI1.GetLongitude(&PuntoXY)
&j=&GxMapPolygon.Paths.Count
for &i= 1 to (&GxMapPolygon.Paths.Count)
    &PuntoI=&GxMapPolygon.Paths.Item(&i).PathLat.Trim()+", "+&GxMapPolygon.Paths.Item(&i).PathLong
    &PuntoJ=&GxMapPolygon.Paths.Item(&j).PathLat.Trim()+", "+&GxMapPolygon.Paths.Item(&j).PathLong
    &PointLatI=&GeoLocationAPI1.GetLatitude(&PuntoI)    
    &PointLongI=&GeoLocationAPI1.GetLongitude(&PuntoI)
    &PointLatJ=&GeoLocationAPI1.GetLatitude(&PuntoJ)
    &PointLongJ=&GeoLocationAPI1.GetLongitude(&PuntoJ)
    if (&PointLongI >= &x)<>(&PointLongJ>=&x)
        &valor= (&PointLatJ - &y) * (&x - &PointLongI) / ( &PointLongJ - &PointLongI) + &PointLatI      
        if (&y < &valor)    
            &salida='S'
            exit
        endif
    endif
    &j=&i
endfor

Pero tengo un problema, me está asumiendo algunos puntos que no están dentro del polígono,y me está haciendo dudar si no le pifié a algo en el algoritmo que no me estoy dando cuenta. Les muestro un ejemplo de lo que detecté con un par de imágenes.

Acá funciona correctamente determina que 2 puntos están en la zona:

introducir la descripción de la imagen aquí

Acá funciona mal los puntos marcados estan fuera de la zona y los toma como que están dentro:

introducir la descripción de la imagen aquí

La pregunta es, ¿Esta correcto el algoritmo que implemente? Parece como si mirando el vértice del triángulo que está mas a la izquierda, si quedan puntos por debajo a la derecha, los incluye como dentro de la zona. Y eso está mal.

Otro dato, las variables que usé son Numeric signed de 16 con 12 posiciones decimales.

Agradezco desde ya alguna ayuda, observación, corrección, etc.

Gracias, saludos Bruno

1

1 respuesta 1

1

Primero que nada, comentarte que en la versión 15 de genexus ya tienes integrada esta funcionalidad a partir del tipo de datos geográfico. (la cual utiliza las funcionalidades de la base de datos para resolverla).

Esto significa que si almacenas los valores en un atributo de tipo Geography, tienes una función Intersect que retorna un booleano, si el punto (u otro tipo geography) está incluido en el polígono. El código queda muy sencillo, sería algo así:

for each 
   where AttGeography.Intersect(&GeoPoint) 
   ...
endfor 

Acá está el detalle en la documentación : Geography data type

En el caso que no sea viable pasar a la última versión y usar el tipo de datos, hay que revisar con mas detalle la implementación de Ray casting que hiciste.

Basándome en la implementación que esta aquí, parecería que no esta iterando de la misma forma. Una forma rápida de probarlo es copiar ese código y embeberlo en un procedimiento, algo del tipo:

&lat = -34.908168
&long =  -56.152393
&vector(1) = "-34.902473,-56.164040"
&vector(2) = "-34.903248,-56.158421"
...
&nvert = len(&vector)
&c = False
csharp for([!&i!] = 0, [!&j!] = [!&nvert!] - 1; [!&i!] = &long )  (&longitude1 >= &long)) and
        (&lat 

Si eso sí te funciona, vemos de resolverlo sin embeber la líneas, en un código genérico multi lenguaje.

1
  • Gracias Javier! Creo que te quedó cortado el código. Pero no tengo c# , uso generador java. Debe estar x ahi el problema en la iteración, no pude entender el "for" cuando dice: for (i = 0, j = nvert-1; i < nvert; j = i++) no me queda claro como itera i y como itera j . Quizá se pueda resolver con Gx con 2 fors anidados xq no se como escribir esa linea en GX.
    – bcamargo75
    Commented el 4 nov. 2016 a las 11:29

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.