Estoy tratando de fusionar 60 archivos geotiff en un archivo enorme para poder crear mosaicos a partir de ese enorme archivo geotiff resultante usando gdal2tiles. El proceso de fusión funcionó, pero por alguna razón produce un archivo que tiene un tamaño de 78 GB en lugar de los esperados 19 GB....
Estoy usando gSimplify(paquete rgeos) para simplificar las geometrías de un shapefile. La función funciona bien, pero ahora no puedo escribir la salida en un nuevo archivo de forma. Intenté algunas formas: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') yo obtengo obj debe...
Estoy tratando de configurar un Python virtualenv aislado para trabajar en proyectos SIG en mi Mac OS X. Parece que necesitaré instalar GEOS, PROJ, GDAL / OGR desde el sitio de Kyngchaos aquí http://www.kyngchaos.com/ software / frameworks y no fácil_instalarlo en mi virtualenv. ¿Es eso...
Intenté usar el código de QGis Save Raster como imagen renderizada . ies una capa de trama: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs =...
Mosaico algunas imágenes con gdal y me gustaría mejorar el resultado final mediante el uso de una capa alfa desvanecimiento / gradual hacia el borde de cada imagen para eliminar los bordes afilados en el medio del mosaico. El problema que tengo es que la parte de cada imagen individual con la capa...
Estoy trabajando en una proyección de selección de sitio donde no podré realizar ninguna observación in situ. Un componente clave del proceso de selección es la robustez del terreno subyacente. Planeo usar GDALdem para crear los mapas de robustez derivados. He leído tanto el Wilson, et al....
Ha habido un montón de preguntas al respecto, pero ninguna de ellas resolvió mi problema. Tengo una máquina con: Windows 7 x64 Python 3.4.3 x64 gdal instalado: gdalinfo --version --> GDAL 1.11.4, released 2016/01/25; gdal se instaló desde una rueda GDAL-1.11.4-cp34-none-win_amd64.whldesde...
ArcGIS Desktop tiene un conjunto de datos de mosaico ráster , QGIS y sus amigos tienen tablas ráster virtuales . Una ojeada de la documentación revela mucha similitud. Ambos crean o habilitan una representación alternativa de datos que se encuentra en otro lugar sin tener que crear una nueva...
Estoy tratando de convertir datos ráster en datos poligonales. Anteriormente estaba usando shapefile.py para crear datos de puntos en lat / lon individuales en el ráster. Pero eso creará una gran cantidad de datos y su consulta será un gran problema. Por lo tanto, sería mejor si pudiera convertir...
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...
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 =...
Cuando intento 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 que instalé GDAL (creo que la biblioteca geos está incluida en ese paquete). Después de eso, verifiqué y tengo...
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...
¿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 usando Python 2.6.5 (32 bits) con Numpy 1.3 y Gdal 1.9.1 instalado en Windows 7 de 64 bits. Estoy tratando de leer un dataset ráster Imagine (.img) de 800 MB en una matriz Numpy para hacer algo de álgebra ráster, pero tan pronto como ejecuto el siguiente código, Python.exe se bloquea. from...
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é...
A veces escribo guiones cortos / instrucciones de trabajo para colegas que no son muy hábiles en GDAL / programación en general. Para evitar problemas con los scripts, al fusionar / convertir geotiffs con GDAL, trato de usar la opción BIGTIFF = YES independientemente de los tamaños de archivo que...
Estoy tratando de usar ogrinfo para obtener algunos detalles sobre un archivo de forma que descargué. Actualmente, la única forma en que sé cómo hacer esto es cargarlo en QGIS y hacer clic manualmente para encontrar cualquier información, como abrir la tabla de atributos. Solo quiero poder ver...
Sé que esa pregunta es bastante vaga, pero por favor tengan paciencia conmigo. Estoy tratando de tener una idea de qué tipo de rendimiento del producto, específicamente el tiempo, las personas han visto para diversas metodologías que han utilizado para crear mosaicos de mapas de google / bing....