Parece que debería ser sencillo, pero no puedo entenderlo por mi vida. Dado un GeoTiff existente que contiene una sola banda, ¿cómo puede agregar otra banda a ese GeoTiff? ¿Es la única forma de hacer esto usando GDALCreate () para crear un nuevo archivo de salida con una banda adicional y luego...
Estoy tratando de instalar GDAL 2.X en CentOS7 sin compilar desde la fuente porque necesitaré configurarlo en otros servidores y la instalación por código fuente será una tarea muy agotadora y lenta. Encontré algunos repositorios como opengeo y ELGIS pero son para CentOS7. En EPEL repo no tiene...
Estoy tratando de optimizar mi código python usando pypy . Si bien puedo importar la mayoría de las bibliotecas y ejecutarlas con éxito, no puedo importar las bibliotecas GDAL. Las mismas bibliotecas cuando se importan en Python 2.7 funcionan bien. Intenté agregar la ubicación de GDAL a PYTHONPATH...
Tengo un ráster entero para el que me gustaría construir una tabla de atributos ráster usando Python y GDAL. Puedo crear una tabla de atributos ráster GDAL en Python de la siguiente manera: >>> rat = gdal.RasterAttributeTable() Esto funciona bien, como podemos ver: >>>...
Solo exporto una tabla Postgis a shp usando estos consejos, pero no puedo importar un shp a Postgis usando la misma biblioteca (ogr). ¿Alguna idea? Muchas gracias
Tengo un archivo de imagen tif y un archivo mundial .tfw (creado a partir de datos OSM con Mapertive). ¿Cómo puedo (usando la línea de comando de Linux) convertir eso a un archivo GeoTIFF? (es decir, quiero tener un archivo, no 2). Estoy seguro de que es posible con gdal o algo...
Soy extremadamente nuevo en SIG. Estoy usando gdalpara leer en un mapa de uso / cobertura del suelo y necesito seleccionar el lat / lng de ciertos tipos de cobertura del suelo para indexar en un conjunto de datos diferente que se expresa solo en lat / lng. Desafortunadamente, no entiendo la forma...
Estoy escribiendo una utilidad simple para recortar lotes de archivos de trama geotiff multibanda en la misma área (más pequeña). Con gdalwarp, puedo recortar fácilmente un archivo usando un archivo de forma de recorte de un solo polígono: gdalwarp -cutline clipper.shp -crop_to_cutline input.tif...
Quería fusionar alrededor de 20 dem tiffs cada uno alrededor de 100 MB, y usé gdalwarp (sé sobre la opción de archivo virtual): gdalwarp $(list_of_tiffs) merged.tiff Sin embargo, este comando tardó mucho tiempo y, después de unos 15 minutos, estaba procesando la quinta imagen del conjunto. Lo...
Estoy fusionando algunos archivos de forma y tuve algunos problemas para hacerlo dentro de QGIS, así que estoy usando ogr2ogr directamente. Estoy haciendo esto (en un lote): ogr2ogr -overwrite %destination% %n1% ogr2ogr -update -append %destination% %n2% -nln all_new ogr2ogr -update -append...
Quiero crear dos polígonos. Uno es de las extensiones rectangulares de mi trama. Sé que puedo obtener la extensión de la trama usando r <- raster("band5.tif") e <- extent(r) plot(e)me da esto: ¿Cómo puedo crear UN polígono del límite del ráster como se muestra a continuación? Así es...
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
Soy un tipo de vectores en su mayoría, pero para un proyecto actual necesito trabajar con un montón de viejos mapas de papel escaneados (¡mapas de daños de bombas ww2 para Londres, si estás interesado!) Escaneamos y georreferenciamos los mapas, y ahora queremos crear una capa compuesta en mosaico...
¿Hay alguna manera de establecer todos los píxeles de valor 0 en un ráster en nodata usando GDAL? Me gustaría usar fillnodata para rellenar algunos agujeros, pero los agujeros están establecidos en 0, por lo que primero debo
Estoy tratando de procesar algunos datos ráster usando ogr / gdal y parece que no puedo obtener la plena utilización de todos los núcleos de mi máquina. Cuando solo ejecuto el proceso en un solo núcleo, obtengo el 100% de utilización de ese núcleo. Cuando trato de dividirme en multinúcleo (en el...
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...
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...
Regularmente creo mis propios rásteres GeoTIFF con GDAL en Python, por ejemplo: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds = driver.Create('MyRaster.tif', 200, 300) band =
Tengo datos ráster en una tabla de postgres que quiero ingresar a python como una matriz numpy. Estoy usando psycopg2 para conectarme a la base de datos. Puedo descargar los datos pero vuelven como una cadena (probablemente un binario serializado). ¿Alguien sabe cómo tomar esta cadena y...
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...