Tengo un geopandas que GeoDataFrame contiene cientos de formas Polygony MultiPolygongeometrías. Los polígonos se superponen en muchos lugares. Me gustaría hacer una nueva geometría que contenga un recuento de cuántas de ellas se superponen. Algo como esto: ¿Alguien tiene alguna idea sobre cómo...
Cerrada . Esta pregunta necesita estar más centrada . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que se centre en un problema solo editando esta publicación . Cerrado hace 4 años . Estoy buscando...
Me gustaría llamar a la función del complemento de interpolación QGIS (método TIN) (Raster-> Interpolate) desde la consola de Python. No puedo encontrar la función correspondiente en QGIS API o en la lista de algoritmos de procesamiento. Encontré el algoritmo de triangulación SAGA, que...
Estoy sirviendo mosaicos vectoriales usando TileStache , tengo todo configurado como quiero. Mis datos se almacenan en Postgres y estoy usando el proveedor VecTiles para servir mosaicos GeoJSON . Quiero almacenar en caché todos mis mosaicos para que los mosaicos se sirvan más rápido. Estoy usando...
En la documentación de gdal_calc se indica Calculadora ráster de línea de comando con sintaxis numpy . Más adelante hay algunos ejemplos donde en uno de ellos: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - significa establecer valores de cero e...
Instalé GDAL y verifiqué en la línea de comandos, funciona y cuando ejecuto gdal2tiles.py en la línea de comandos funciona, pero cuando ejecuto el mismo archivo gdal2tiles.py con Python IDLE, obtengo "Ningún módulo llamado 'osgeo '"error. ¿Puedes
Logré agregar QgsVectorLayera QGIStravés de un complemento de Python, por ejemplo: vl = self.iface.addVectorLayer(uri.uri(), layerName, self.dbConn.getProviderName()) Pero, ¿cómo puedo
Vea este enlace para más detalles. El problema: Quiero recorrer un ráster continuo (uno que no tiene tabla de atributos), celda por celda, y obtener el valor de la celda. Quiero tomar esos valores y ejecutar condicionales sobre ellos, emulando los pasos de álgebra de mapas detallados a...
usando Gdal / OGR, ¿hay alguna forma de convertir una capa de línea shp, kml o PostGIS en una red de nodos y segmentos para que pueda usarse en paquetes como
[EDITAR: la solución a esto fue simplemente usar OGR para leer archivos de forma. Ver el ejemplo de geographika.] En un archivo de forma ESRI, no hay distinción entre Polígonos y MultiPolygons. Además, no existe una distinción explícita entre los agujeros interiores y los anillos exteriores...
Tengo un dataset ráster en formato ASCII. Usando Python, necesito determinar los valores miny maxdentro del conjunto de datos. Me han dicho que la información del encabezado es clave, que contiene elementos como el número de filas / columnas, el tamaño de la celda, etc. ¿No puede simplemente...
Me gustaría construir un shapefile a partir de un Marco de Datos Pandas usando las filas lon y lat. Tengo un archivo csv y lo proceso con pandas para crear un marco de datos que sea más fácil de manejar ¿Es posible hacer eso sin hacer un bucle línea por
¿Existe un paquete de Python que proporcione una implementación de un algoritmo de esqueleto recto? Soy consciente de que el proyecto CGAL de código abierto (C ++) contiene una implementación, pero parece que cgal-bindings no incluye este paquete CGAL . En cualquier caso, preferiría una...
Estoy interesado en aprender métodos para utilizar toda la potencia de procesamiento multinúcleo disponible en una computadora de escritorio. Arc afirma que el geoprocesamiento en segundo plano permite al usuario utilizar múltiples núcleos, sin embargo, las tareas esencialmente tienen que esperar...
Estoy escribiendo un código de Python para convertir los puntos X, Y en Excel a shapefile. En este proceso, tengo que leer los puntos xy de shhet1. ¿Cuál es el paso que debo incluir en mi proceso para leer datos de la hoja 1 de un libro de trabajo de Excel (97-2003)? Escribí un código de la...
Estoy trabajando con GDAL en Python y he querido que GDAL genere excepciones para poder hacer algo útil con ellos. Por ejemplo, si un archivo no se puede abrir con gdal.Open (), me gustaría lanzar una excepción en lugar de que gdal devuelva None. Solo me gustaría una forma de capturar los errores...
Tengo problemas para usar la función Polygonize en python. El ejemplo del libro de cocina para esto se puede encontrar aquí . La parte relevante de mi código es: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile =...
Estoy tratando de grabar un shapefile en un ráster usando RasterizeLayer de GDAL. Precreo un ráster de área de interés desde un archivo de formas diferente, dado un tamaño de píxel específico. Este AOI sirve entonces como base para todas las siguientes rasterizaciones (mismo número de columnas y...
Me gustaría actualizar un atributo de una característica. Sin embargo, no sé actualizarlo usando el objeto de función. Tengo que usar el proveedor de datos para actualizarlo. layers = QgsMapLayerRegistry.instance().mapLayersByName('my_line') layer = layers[0] dp = layer.dataProvider() it =...
¿Cuál es el mejor método para calcular la distancia desde la carretera más cercana para una serie de puntos de dirección en QGIS, con un conjunto de datos muy grande? (QGIS 2.0.1, Win XP) Leí las respuestas a preguntas similares anteriores: Medición de la distancia entre líneas y puntos en...