Rééchantillonner les images GeoTIFF à la même résolution?

21

J'ai deux images géotiff que je voudrais rééchantillonner la même résolution. Si j'utilise ArcGIS, ce serait facile à faire. Mais en utilisant QGIS, je ne sais pas comment faire. Une aide sera très appréciée. Merci.

Gil Gonçalves
la source

Réponses:

29

C'est aussi facile dans QGIS, mais un peu moins évident. Vous pouvez le faire de plusieurs manières:

  • Calculatrice raster - utilisez simplement la calculatrice raster et vous pouvez y définir la résolution et l'étendue et les faire correspondre à un autre raster en sélectionnant la bande raster que vous souhaitez faire correspondre dans la liste des bandes raster, puis en cliquant sur le bouton "Etendue de la couche actuelle". Les champs de colonnes et de lignes vous permettront de définir la résolution. Cependant, cette méthode ne vous donne aucun contrôle sur la méthode de rééchantillonnage.
  • Utilisation de GDAL_Warp - cet outil vous permet de définir la résolution de sortie soit en spécifiant la largeur et la hauteur du raster en sortie, soit en spécifiant le commutateur -tr (voir la documentation). Vous pouvez accéder à l'outil GDAL_warp en allant Raster-> Projections-> Warp (j'ai dit que ce n'était pas évident du point de vue du rééchantillonnage!).

    • (v2.x) Si vous souhaitez utiliser le commutateur -tr, remplissez toutes les cases pour le raster en entrée et la sortie, etc. (vos valeurs SRS source et cible seront probablement les mêmes dans ce cas - mais cela ne doit pas être vous reprojetez aussi). Cliquez ensuite sur la petite icône en forme de crayon en bas et modifiez la ligne de commande gdal-warp générée automatiquement pour inclure votre commutateur -tr. Gdal_wrap vous permet de spécifier l'algorithme que vous souhaitez utiliser pour le rééchantillonnage et est donc un peu moins un instrument contondant que d'utiliser la calculatrice raster.
    • (v3.x) Le commutateur -tr est activé à l'aide de la zone Résolution du fichier de sortie dans les unités cibles géoréférencées . Par exemple, pour sous-échantillonner un DEM de 1 m à un DEM de 2 m, vous pouvez entrer 2 dans ce champ. Cependant, il n'y a pas d'option pour passer deux arguments différents pour les pixels non carrés. Supposons que votre taille de pixel cible soit 0.3125,0.25, c'est-à-dire l' xresest 0.3125et l' yresest 0.25. Si vous passez maintenant la valeur 0.3125dans cette case, elle sera définie -tr 0.3125 0.3125dans la commande. Pour contrer cette limitation, copiez simplement le code, collez-le sur la ligne de commande, modifiez l'indicateur -tr et exécutez. Par exemple:

      gdalwarp -t_srs EPSG:4326 -tr 0.3125 0.25 -r near -te 71.40625 24.875 84.21875 34.375 -te_srs EPSG:4326 -of GTiff foo.tiff bar.tiff

      (en fonction de vos variables d'installation et d'environnement, vous devrez peut-être également indiquer explicitement le chemin d'accès à gdalwarp).

MappaGnosis
la source
1
La 2ème option ne fonctionne que dans QGIS v. 2.x, l'édition des commandes GDAL n'est pas disponible dans v. 3x issues.qgis.org/issues/15090
reima
C'est vrai, malheureusement! Cependant, vous pouvez utiliser gdal_warp à partir de la ligne de commande. Avec QGIS installé, vous aurez accès à cela, donc, bien que gênant, la deuxième option est toujours disponible.
MappaGnosis
3

J'utilise normalement l'outil d'alignement raster.

Il est localisé dans l'onglet raster et peut également être utilisé pour définir les rasters dans la même mesure. Il peut également effectuer certaines opérations lors du rééchantillonnage à une nouvelle résolution (comme la valeur moyenne, min et max, etc.).

Voici un lien avec des informations à ce sujet: https://docs.qgis.org/2.18/en/docs/user_manual/working_with_raster/raster_analysis.html?highlight=raster%20alignment#id3

J'utilise la version 2.18.13. Je ne sais pas si l'outil est disponible dans les anciennes versions.

Gabriel Guariglia Perez
la source
1
Cela ne semble pas fonctionner pour le rééchantillonnage - ce qui réduit la résolution d'un raster.
Jakub Sisak GeoGraphics
étrange, cela a fonctionné pour moi la dernière fois que j'ai utilisé. Cela peut être un petit buggy.
Gabriel Guariglia Perez
J'ai utilisé Warp (reprojet) et cela a fonctionné comme un charme avec des agrégats et d'autres options ...
Jakub Sisak GeoGraphics