Questions marquées «gdalwarp»

J'ai une image raster géotiff qui a un système de coordonnées avec des longitudes de 0 à 360. Le centre horizontal de l'image est 180 longtitude. Voir l'image ci-dessous: Je veux le transformer en EPSG: 4326 SRS avec une plage de -180 180 longtitude. Et je veux que le centre de l'image soit au...

J'essaie de découper mon fichier raster selon une géométrie de fichier de formes. J'utilise le code suivant gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif Mais il en résulte une couleur noire en dehors de l'étendue de la géométrie du fichier de formes. Je donne l'exemple ici. Dans la première...

J'essaie de re-projeter / rééchantillonner avec les liaisons Python GDAL, mais j'obtiens des résultats légèrement différents par rapport à ceux de l'utilitaire de ligne de commande gdalwarp. Voir la mise à jour ci-dessous pour un exemple plus court Ce script illustre l'approche Python: from osgeo...

Après avoir utilisé gdalwarppour projeter et aligner sur la grille (via -tap) un certain nombre de rasters, j'ai remarqué que les rasters en sortie étaient beaucoup plus volumineux que les rasters d'origine. Une recherche Web assez approfondie a révélé ce problème Trac : Frank Warmerdam a expliqué...

Je voulais fusionner environ 20 fichiers de démonstration chacun autour de 100 Mo, et j'ai utilisé gdalwarp (je connais l'option de fichier virtuel): gdalwarp $(list_of_tiffs) merged.tiff Cependant, cette commande a pris très longtemps, et après environ 15 minutes, elle traitait la 5ème image de...

J'écris un utilitaire simple pour recadrer des lots de fichiers raster géotiff multi-bandes dans la même zone (plus petite). En utilisant gdalwarp, je peux facilement recadrer un fichier à l'aide d'un fichier de formes de découpage à un seul polygone: gdalwarp -cutline clipper.shp -crop_to_cutline...

Dans QGIS 2.XX, les outils de la barre de menu déroulante avaient une zone d'appel modifiable en bas, où les commandes non disponibles dans l'interface graphique pouvaient être saisies manuellement (en cliquant sur un bouton de crayon). Dans 3.0, les outils de la barre de menu sont les mêmes que...

Je mosaïque certaines images avec gdal et j'aimerais améliorer le résultat final en utilisant une couche d'altération progressive / graduelle vers le bord de chaque image pour supprimer les bords nets au milieu de la mosaïque. Le problème que j'ai, c'est que la partie de chaque image individuelle...

J'ai essayé d'utiliser le code de QGis Save Raster as Rendered Image . iest une couche raster: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs =...

J'utilise gdalwarp pour manipuler des tuiles SRTM près de la ligne de données (c'est-à-dire 180 °, alias l'antiméridien). Les carreaux SRTM ont un très léger chevauchement (1/2 pixel) avec le méridien. Vous pouvez le voir en utilisant gdalinfo: gdalinfo S16W180.hgt Driver: SRTMHGT/SRTMHGT File...

J'essaie de remapper de Geostationary à Lambert conforme en utilisant gdalwarp. Mes données d'entrée sont en netcdf, et sont en coordonnées géographiques (degrés) et je voudrais sortir les données remappées vers netcdf. J'ai créé un fichier vrt correspondant pour les données d'entrée netcdf....

Je souhaite géoréférencer un raster à l'aide de pythonet GDAL. Mon approche actuelle est d'appeler gdal_translateet d' gdalwarputiliser os.systemet une vilaine liste de points de contrôle au sol. J'aimerais vraiment un moyen de le faire nativement à l'intérieur python. Voici le processus que...

J'ai du mal à générer des tuiles pour une image haute résolution que nous avons. L'image actuelle que nous avons est une très grande image (+ 20 Go), enregistrée en tant que fichier GeoTiff. Je voudrais générer les tuiles à l'aide de l'utilitaire de ligne de commande gdal2tiles, puis l'ouvrir et...

J'utilise python et QGIS 2.0. J'essaie de découper les rasters dans un dossier par une entité polygonale. C'est la première fois que j'utilise (disons) "PyQGIS", j'étais habitué à arcpy avant. Quoi qu'il en soit, je ne fais pas fonctionner mon script simple, toute suggestion serait très appréciée!...

J'essaie de projeter un raster. Dans R, il y a la projectRaster()fonction à cela (ci-dessous un exemple entièrement reproductible): # 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 to...

J'essaie d'agréger certaines données raster (toutes .tif) jusqu'à une résolution plus grossière (de 0,05 degré à 0,25 degré) en utilisant gdalwarp en python, mais la commande ne fonctionne pas. Au lieu d'obtenir une sortie avec une large plage de valeurs, toutes les valeurs de la sortie sont 0. La...