Actualmente estoy experimentando un error de 'Acceso denegado' después de usar rasterio para abrir los archivos de imagen Sentinel-2 AWS .jp2 (actualmente alojado por Sinergise y se puede encontrar aquí ). Mi código se ve así: import rasterio access_key = '*****************' secret_access_key =...
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...
He estado aprendiendo geoprocesamiento con Python usando el curso GEOG 825 (geoprocesamiento de python) de MSU . Lamentablemente, el enlace está muerto ahora. Sé que esto es una posibilidad remota, pero ¿alguien tiene una copia de los materiales de ese curso? Realmente me gustaría terminarlo. Ya...
Mi gerente me ha pedido que presente algún tipo de medida de competencia para que los empleados comprendan las secuencias de comandos / programación. Debido a que soy un programador autodidacta, realmente no tengo una comprensión formal de cómo medir el nivel de comprensión cuando se trata de...
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...
¿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 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...
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...
Instalé algunas bibliotecas de Python usando pip (pandas, scipy, simplekml) y desde entonces cada vez que intento cargar QGIS me sale este error: Warning: loading of qgis translation failed [/usr/share/qgis/i18n//qgis_en_US] Warning: loading of qt translation
Me gustaría obtener la intersección de múltiples polígonos. Usando el shapelypaquete de Python , puedo encontrar la intersección de dos polígonos usando la intersectionfunción. ¿Existe una función eficiente similar para obtener la intersección de múltiples polígonos? Aquí hay un fragmento de...
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
¿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...
Tengo un Arc / Info Binary Grid --- específicamente, un ráster de acumulación de flujo ArcGIS --- y me gustaría identificar todas las celdas que tienen un valor específico (o en un rango de valores). En última instancia, me gustaría un archivo de formas de puntos que representen estas...
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 buscando un método para calcular la coordenada X, Y de un punto interno (también conocido como labelPoint en arcmap) de un polígono (irregular). Estoy tratando de encontrar una solución sin el uso del módulo ArcGIS La diferencia entre un "centroide" y un labelPoint es que el centroide...
Estoy realmente sorprendido de los enormes pasos que ha dado el desarrollo de QGIS en los últimos meses. ¡El nuevo conjunto de iconos, el mejorado compositor de impresión y muchos más cambios que se verán en QGIS 2.0 son realmente increíbles! Entre muchas otras mejoras interesantes también está...
¿Hay alguna manera de realizar la misma tarea que la utilidad gdalbuildvrt usando enlaces GDAL Python? Hasta ahora no he encontrado ninguna manera de hacer esto que no sea crear un vrt de un único conjunto de datos y editar manualmente el xml. Me gustaría crear un vrt a partir de múltiples rásteres...
Estoy leyendo un archivo de forma usando Fiona: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] esto me da un dict que contiene una lista de listas de coordenadas y escribe 'MultiPolygon'. Sin embargo, no estoy seguro de...
Pasando de la lista de correo de gdal-dev: El lunes 2 de septiembre de 2013 a las 7:09 p.m., David Shean escribió: Hola lista, estoy tratando de empaquetar una serie temporal de rásteres GTiff con idéntica proyección / extensión / resolución como un único archivo NetCDF para su distribución. Pasé...
Tengo una capa intermedia (polígono verde) que quiero dividir en dos polígonos cada vez que cruza una barrera (línea azul). He estado tratando de usar el método "splitGeometry", pero no puedo hacerlo funcionar. Mi código hasta ahora es este: while ldbuffprovider.nextFeature(feat): while...