8 votos

QGIS: Crear archivo de formas poligonales de elipses a partir de puntos usando la herramienta de geometría por expresión

En QGIS 3.16.6 necesito crear un shapefile de polígonos de elipses a partir de puntos. Estoy haciendo lo siguiente:

Caja de herramientas de Procesamiento -> Geometría de Vectores -> Geometría por expresión

  • Capa de entrada - points_in.shp

  • Tipo de geometría de salida - polígono

  • Expresión de geometría - make_ellipse(transform(make_point("Lon-DD","Lat-DD"), 'EPSG:4269', 'EPSG:26913'), 635, 190, 90)

  • Geometría modificada - ellipses_out.shp

Donde "Lon-DD" y "Lat-DD" son campos en el shapefile de puntos, el SRC de entrada es GDD NAD83, el SRC de salida es NAD83/Zona UTM 13N en metros, y los ejes semi-mayor y semi-menor son 635 metros y 190 metros y el azimut es 90 grados.

La herramienta Geometría por expresión se ejecuta sin quejas pero el shapefile de salida está vacío.

¿Qué estoy haciendo mal y cómo puedo arreglarlo?

1 votos

¿Revisaste la tabla de atributos de la capa "ellipses_out"? ¿Hiciste zoom a la capa? ¿Significa "Vacio" que no ves elipses en el mapa o que la tabla de atributos de "ellipses_out" está vacía?

0 votos

Ampliar a la capa de salida es diferente que ampliar a la capa de entrada. La capa de salida está ahí, con más de 200 elipses. Entonces, ¿ya sé cuál es el problema, ahora cuál es la solución? Nuevamente, transformar ('EPSG:4269', 'EPSG:26913') parece que no está funcionando.

10voto

nitinsavant Puntos 6

El problema es que esperas que make_ellipse devuelva una elipse en coordenadas EPSG:26913. Sí, devuelve coordenadas EPSG:26913 (números grandes como 500000 que son coordenadas proyectadas), pero QGIS considera que están en el CRS de la capa fuente, EPSG:4269 (latitud, longitud). Si verificas el CRS de ellipses_out, verás que su CRS es EPSG:4269.

Para una solución, usa la herramienta Procesamiento > Vector general > Asignar proyección para la capa de salida. Elige EPSG:26913 como nuevo CRS.

introducir descripción de la imagen aquí

0 votos

Correcto. Gracias por identificar el problema, pero lo que pregunté fue por una solución: "¿Qué estoy haciendo mal y cómo lo arreglo?" Aparentemente transform('EPSG:4269', 'EPSG:26913') no está funcionando. ¿Supongo que si proyecto los puntos de entrada desde GDD a UTM primero, entonces esto debería funcionar?

0 votos

¡Gracias, problema resuelto! ¿Cómo marco esta pregunta como respondida y completada?

1 votos

Flecha hacia arriba y marca de verificación ¡Hecho!

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X