Articles of gis

La forma más sencilla de guardar una matriz en un archivo raster en Python

Con una matriz 2-d en forma de (100, 100), quiero guardarla en un archivo raster en formato .tiff . Puedo usar el paquete gdal para leer archivos tiff que ya existen. Pero aún no puedo encontrar una forma sencilla de transformar la matriz 2-d en un archivo tiff . Usando plt.imsave(“xx.tif”,array) o def array_to_raster(array): “””Array […]

Pandas: calculan la distancia del tallo dentro de cada grupo de filas

El CSV de muestra es así: user_id lat lon 1 19.111841 72.910729 1 19.111342 72.908387 2 19.111542 72.907387 2 19.137815 72.914085 2 19.119677 72.905081 2 19.129677 72.905081 3 19.319677 72.905081 3 19.120217 72.907121 4 19.420217 72.807121 4 19.520217 73.307121 5 19.319677 72.905081 5 19.419677 72.805081 5 19.629677 72.705081 5 19.111860 72.911347 5 19.111860 72.931346 5 […]

Calcular la distancia a la costa o la costa de un barco

Para un conjunto de datos de 200M GPS (lon, lat) coordenadas de embarcaciones, quiero calcular una distancia aproximada a la tierra o costa más cercana, como una función llamada distance_to_shore, que devolverá la distancia y el país de esa costa. Estoy usando un archivo de formas de los límites de los países y las líneas […]

¿Estructura de datos para realizar búsquedas rápidas de GPS?

Tengo un archivo de texto (UTF-8, ~ 50K líneas) con nombres de ciudades y coordenadas GPS. Líneas de ejemplo: San Pedro locality -3367 -5968 Argentina Buenos Aires San Pedro Talagante locality -3366 -7093 Chile Metropolitana Talagante Peñaflor locality -3362 -7092 Chile Metropolitana Talagante La tercera y cuarta columnas son las coordenadas GPS de las ciudades […]

Rectangularjs desde puntos usando Python

Tengo un archivo de texto lleno de puntos. Están separados en cada línea por un par limitado por comas (x, y). p.ej. -43.1234,40.1234\n -43.1244,40.1244\n etc. Ahora necesito crear un polígono alrededor de cada uno de estos puntos. El polígono debe tener un búfer de 15 kilómetros desde el punto. No tengo acceso a ArcGIS ni […]

Python, GEOS y Shapely en Windows 64

Al intentar instalar Shapely en mi computadora con Windows de 64 bits, no puedo hacer que la biblioteca GEOS funcione. Hasta ahora, he ejecutado el instalador OSGeo4W desde el cual instalé GDAL (creo que la biblioteca geos está incluida en ese paquete). Después de eso, verifiqué y tengo geos_c.dll en mi directorio C:\OSGeo4W\bin , pero […]

La consulta a distancia de Geodjango no recupera resultados correctos

Estoy tratando de recuperar algunos mensajes en función de su proximidad geográfica. Como puede ver en el código, estoy usando GeoDjango y el código se ejecuta dentro de una vista. El problema es que el filtro de distancia parece estar completamente ignorado. Cuando verifico las distancias en el queryset, obtengo las distancias esperadas (1m y […]

Cómo cambiar la resolución del archivo raster .tif sin perder datos

Por ejemplo, en un archivo .tiff de tamaño 400 * 200, puedo leerlo como una matriz 2D (400 x 200) en python. Quiero cambiar el tamaño de tiff a 200 x 100 u otra relación. Cómo darse cuenta de eso en Python o software GIS (QGIS, ArcGIS et.al). Sé que la herramienta de parecido en […]

Mejor biblioteca de Python GIS?

Estoy buscando recomendaciones sobre la mejor biblioteca de Python GIS actualmente disponible según los siguientes criterios (subjetivos): Facilidad de uso / Interfaz Pythonic / Documentación Potencia en términos de funciones disponibles en la actualidad, y potencialmente en el futuro (aún se está desarrollando o está muerto). Hasta ahora he mirado, con algún detalle: GEOS / […]

Encuentra todas las coordenadas dentro de un círculo en datos geográficos en python

Tengo millones de puntos geográficos. Para cada uno de estos, quiero encontrar todos los “puntos vecinos”, es decir, todos los otros puntos dentro de algún radio, digamos unos pocos cientos de metros. Hay una solución ingenua de O (N ^ 2) a este problema: simplemente calcule la distancia de todos los pares de puntos. Sin […]