Preguntas etiquetadas con gdalwarp

Tengo una imagen ráster geotiff que tiene un sistema de coordenadas con longitudes de 0 a 360. El centro horizontal de la imagen es 180 de longitud. Ver imagen a continuación: Quiero transformarlo a EPSG: 4326 SRS con -180 180 rango de longitud. Y quiero que el centro de la imagen esté en el...

Estoy tratando de recortar mi archivo Raster de acuerdo con una geometría de shapefile. Estoy usando el siguiente código gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif Pero está resultando un color negro fuera de la extensión de la geometría del archivo de forma. Estoy dando el ejemplo aquí. En...

Estoy tratando de volver a proyectar / volver a muestrear con los enlaces de Python GDAL, pero obtengo resultados ligeramente diferentes en comparación con los de la utilidad de línea de comandos gdalwarp. Vea la actualización a continuación para ver un ejemplo más corto Este script ilustra el...

Después de usar gdalwarppara proyectar y alinear a la cuadrícula (a través de -tap) una serie de rásteres, noté que los rásteres de salida eran significativamente más grandes que los rásteres originales. Una búsqueda en la web bastante exhaustiva reveló este problema de Trac : Frank Warmerdam...

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 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...

En QGIS 2.XX, las herramientas de la barra de menú desplegable tenían un cuadro de llamada editable en la parte inferior, donde los comandos no disponibles en la interfaz gráfica de usuario podían ingresarse manualmente (haciendo clic en un botón de lápiz). En 3.0, las herramientas de la barra de...

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...

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 =...

Estoy usando gdalwarp para manipular mosaicos SRTM cerca de la línea de fecha (es decir, 180 °, también conocido como el antimeridiano). Las baldosas SRTM tienen una superposición muy leve (1/2 píxel) con el meridiano. Puedes ver esto usando gdalinfo: gdalinfo S16W180.hgt Driver: SRTMHGT/SRTMHGT...

Estoy tratando de reasignar de geoestacionario a Lambert conforme usando gdalwarp. Mis datos de entrada están en netcdf, y están en coordenadas geográficas (grados) y me gustaría enviar los datos reasignados a netcdf. He creado un archivo vrt correspondiente para los datos de entrada netcdf....

Quiero georreferenciar un ráster usando pythony GDAL. Mi enfoque actual es llamar gdal_translatey gdalwarpusar os.systemy una lista fea de puntos de control en tierra. Realmente me gustaría una forma de hacer esto de forma nativa python. Este es el proceso actual que estoy usando: import...

He estado luchando un poco para generar mosaicos para una imagen de alta resolución que tenemos. La imagen actual que tenemos es una imagen muy grande (+ 20GB), guardada como un archivo GeoTiff. Me gustaría generar los mosaicos utilizando la utilidad de línea de comando gdal2tiles y luego abrirlos...

Estoy usando python y QGIS 2.0. Estoy tratando de recortar rásteres en una carpeta por una función de polígono. Es la primera vez que uso (digamos) "PyQGIS", estaba acostumbrado a hacer arcpy antes. De todos modos, no logro que mi script simple funcione, ¡cualquier sugerencia sería muy...

Estoy tratando de proyectar una trama. En R existe la projectRaster()función para esto (debajo de un ejemplo completamente reproducible): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project...

Estoy tratando de agregar algunos datos ráster (todos .tif) hasta una resolución más gruesa (de .05 grados a .25 grados) usando gdalwarp en python, pero el comando no funciona. En lugar de obtener una salida con una amplia gama de valores, todos los valores de la salida son 0. La resolución y la...