Articles of gis

Equivalente de PostGIS ST_MakeValid en Django GEOS

Al verificar la validez de un objeto de polígono utilizando Objects.polygon.valid , se GEOS_NOTICE: Self-intersection un GEOS_NOTICE: Self-intersection . Sé que esto puede solucionarse utilizando el método ST_MakeValid de PostGIS. Estoy usando Django 1.11 con el soporte de GEOS y no puedo encontrar su equivalente en los documentos de Django. ¿Hay alguna función equivalente para […]

insertar valores de variables en una cadena en python

Como introduce una variable [i] en una cadena en python. Por ejemplo, mire la siguiente secuencia de comandos, solo quiero poder darle un nombre a la imagen, por ejemplo geo [0]. Tiff … a geo [i]. tiff, o si utiliza un contador, ya que puedo reemplazar una parte de la cadena de valor para generar […]

¿Cómo convierto kilómetros a grados en Geodjango / GEOS?

Estoy usando el método “buffer” de la API de GEOS para crear un círculo basado en un punto y un radio en GeoDjango, de acuerdo con esta respuesta: GeoDjango: Cómo crear un círculo basado en el punto y el radio Como señaló un comentarista, tendrás que convertir el radio en km a grados … ¿pero […]

Cómo calcular la distancia 3D (incluida la altitud) entre dos puntos en GeoDjango

Prólogo: Esta es una pregunta que surge a menudo en SO: Cálculos en 3D con GeoDjango. Cálculo de la distancia entre dos puntos utilizando la latitud, la longitud y la altitud (elevación) Distancia entre dos puntos 3D en geodjango (postgis) Quería redactar un ejemplo en la Documentación SO, pero el capítulo de geodjango nunca despegó […]

Cómo convertir de UTM a LatLng en python o Javascript

Tengo un montón de archivos con coordenadas en forma UTM. Por cada coordenada tengo este, norte y zona. Necesito convertir esto a LatLng para usar con la API de Google Map para mostrar la información en un mapa. He encontrado algunas calculadoras en línea que hacen esto, pero no código o bibliotecas reales. http://trac.osgeo.org/proj4js/ es […]

¿Cómo mostrar la etiqueta de shapefile en la leyenda del mapa base de python?

El mapa base de python tiene dos contenidos trazados: un shapefile multilínea (IL_State_ln) y un diagtwig de dispersión de algunos puntos aleatorios dentro de la extensión del mapa base. Mi interés está en generar una leyenda que brinde información sobre el shapefile y los puntos de dispersión. Por ahora, solo puedo incluir los puntos en […]

¿Matplotlib imshow offset to match axis?

Estoy trazando un montón de coordenadas UTM utilizando un matplotlib.pyplot.scatter. También tengo una foto aérea de fondo que sé que coincide exactamente con la extensión de la figura. Cuando grafico mis datos y configuro el eje, puedo mostrar la dispersión correctamente. Si grafico la foto aérea utilizando imshow, usa el número de píxeles como la […]

Obtener lat / long dado el punto actual, la distancia y el rodamiento

Dado un punto existente en lat / long, distancia en (en KM) y rumbo (en grados convertidos a radianes), me gustaría calcular el nuevo lat / long. Este sitio aparece una y otra vez, pero no puedo conseguir que la fórmula funcione para mí. Las fórmulas que se han tomado del enlace anterior son: lat2 […]

La forma más rápida de obtener todos los puntos entre dos coordenadas (X, Y) en python

Así que tengo un LineString shapely LineString : print np.round(shapely_intersecting_lines.coords).astype(np.int) >>> array([[ 1520, -1140], [ 1412, -973]]) Esto se puede interpretar como una matriz numpy como se ve arriba. Quiero obtener todos los puntos intermedios, es decir, obtener los puntos de la línea intermedios como valores enteros. La salida debería ser algo como esto: array([[ […]