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 …
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 image que je …
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 import osr, gdal import numpy …
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é la raison: "Après …
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 l'ensemble. Je le résilie. J'ai ensuite …
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 input.tif output.tif Cependant, la zone réelle …
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 avec la …
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 Format Files: S16W180.hgt Size is 1201, 1201 …
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. Gdalwarp affichera le fichier netcdf, mais les …
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'utilise actuellement: import os os.system('gdal_translate -of GTiff …
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 …
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 newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- …
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! import qgis.core, qgis,utils …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.