Sistemas de Información Geográfica

9
Leer desde vistas generales de TIF usando rasterio

Estoy leyendo una ventana de un tif usando rasterio. ¿Cómo utilizo las vistas generales disponibles para obtener una versión de muestra reducida de los datos en lugar de leer en la resolución verdadera? En este momento estoy leyendo una ventana como esta: rmin, cmin = source.index(xmin,...

9
¿Eliminar una banda del archivo ráster?

Estaba trabajando con imágenes de drones (cámara normal) y existe esta banda alfa que me impide hacer una clasificación de máxima probabilidad. Quiero eliminar la banda alfa. Hasta ahora utilicé la calculadora ráster, separé la imagen en 3 imágenes para cada banda (R, G, B) y las usé con la fusión...

9
¿Cómo particiono espacialmente muchos puntos?

Tengo datos de ubicación de clientes, con más de 130 millones de registros, repartidos por todo el país. Esta es una tabla espacial en PostGIS Ahora quiero dividir todo el país en 'áreas de servicio', de modo que cada área de servicio tenga: Al menos 30 clientes Tiene como máximo 1000 km2. He...

9
QGIS 2.16 Merge carece de la opción de carpeta

He instalado QGIS 2.16 y he notado que la combinación (menú de vectores> gestión de datos) ya no parece tener una opción de carpeta. Esta fue la opción que le permitió procesar / fusionar una carpeta completa sin tener que cargar capas en su proyecto / espacio de trabajo primero. He echado un...

9
ST_Intersecta con LINESTRING degenerado

Tengo, en PostGIS, un montón de objetos LINESTRING Z (que representan polos), y quiero encontrar cuáles están dentro de un área (representados por un POLÍGONO). Para el propósito de este ejercicio, podemos asumir con seguridad que un poste es bastante vertical, por lo que no se cruza con el límite...

9
¿Duplicar la capa en la memoria usando pyqgis?

Tengo una capa en QGIS, y quiero duplicarla a través de un complemento para poder usar la copia que desee, sin modificar el original. Por supuesto layer2 = layer1, no funcionará, porque todo lo que le sucede a la capa2 también le ocurrirá a la capa1, ya que es el mismo objeto detrás de todo...

9
ggmap: trazar polígono desde un archivo de forma

Usando ggmap, quiero incluir los bordes del municipio (polígono) de un archivo de forma en un mapa que contiene algunos puntos de ubicación. Este script hace todo excepto trazar el polígono: library(rgdal) library(ggmap) # Get shapefile with Drammen municipality