Articles of gis

Encuentre la línea más cercana a cada punto en el conjunto de datos grande, posiblemente utilizando shapely y rtree

Tengo un mapa simplificado de una ciudad que tiene calles como líneas y direcciones como puntos. Necesito encontrar el camino más cercano desde cada punto a cualquier línea de la calle. Tengo un script de trabajo que hace esto, pero se ejecuta en tiempo polinomial ya que se ha nested para bucle. Para 150 000 […]

Encuentra las coordenadas de destino dadas las coordenadas iniciales, el rumbo y la distancia.

Estoy viendo la fórmula que aparece aquí: http://www.movable-type.co.uk/scripts/latlong.html Parece que estoy teniendo problemas ya que las coordenadas resultantes no son lo que yo esperaría que fueran. Dada la siguiente información: Lat de inicio: 28.455556 Inicio lon: -80.527778 Rodamiento: 317.662819 (Grados) Distancia: 130.224835 (millas náuticas) def getEndpoint(lat1,lon1,bearing,d): R = 6378.1 #Radius of the Earth brng = […]

¿Cómo puedo obtener la versión instalada de GDAL / OGR de python?

¿Cómo puedo obtener la versión instalada de GDAL / OGR de python? Soy consciente del progtwig gdal-config y actualmente gdal-config usando lo siguiente: In [3]: import commands In [4]: commands.getoutput(‘gdal-config –version’) Out[4]: ‘1.7.2’ Sin embargo, sospecho que hay una forma de hacerlo usando la API de python. Cualquier dados?

Cómo utilizar la ruta de acceso recortada para el polígono de mapa base

Quiero usar imshow (por ejemplo) para mostrar algunos datos dentro de los límites de un país (para los fines de ejemplo, elegí los EE. UU.) El siguiente ejemplo ilustra lo que quiero: import numpy as np import matplotlib.pyplot as plt from matplotlib.patches import RegularPolygon data = np.arange(100).reshape(10, 10) fig = plt.figure() ax = fig.add_subplot(111) im […]

Amortiguamiento geodésico en python

Dados los polígonos terrestres como un MultiPolygon MultiPolygon , quiero encontrar el Polígono (Multi-) que representa, por ejemplo, el amortiguador de 12 millas náuticas alrededor de las líneas costeras. El uso del método de buffer Shapely no funciona, ya que utiliza cálculos euclidianos. ¿Alguien puede decirme cómo calcular los búferes geodésicos en python?

Interpolación local de orden superior de curvas implícitas en Python

Dado un conjunto de puntos que describen alguna trayectoria en el plano 2D, me gustaría proporcionar una representación fluida de esta trayectoria con interpolación local de alto orden. Por ejemplo, digamos que definimos un círculo en 2D con 11 puntos en la siguiente figura. Me gustaría agregar puntos entre cada par de puntos consecutivos en […]

Cálculos en 3D con GeoDjango.

estoy usando python 2.7.12 django 1.10.6 postgreSQL 9.5.6 postGIS 2.2.2 Primera pregunta Necesito usar GeoDjango para calcular la distancia entre dos puntos. Cuando verifiqué la documentación , dice que GeoQuerySet.distance () está en desuso y en su lugar usa Distancia () de django.contrib.gis.db.models.functions . El siguiente código funciona bien: from django.contrib.gis.db.models.functions import Distance p1 = […]

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 […]