Suite

Créer GeoTIFF à partir de TIFF en utilisant GDAL ?


J'ai une image tiff et j'aimerais la géoréférencer, et ainsi créer un Geotiff. Pour ce faire, je dois utiliser GDAL. J'utilise notamment, l'exécutable store dans le fichier QGIS (gdalwarp.exe, gdaltranslate.exe… ).

Pour l'instant, j'utilise ce code dans une commande DOS (ligne de commande), et cela ne semble pas fonctionner.

gdal_translate -of GTiff -gcp 33.9983 669.668 -833589 -692.339 -gcp 124.661 1377.97 -833520 -1379,2 -gcp 6788.34 -10.2993 -826876 6.06375 -gcp 6839.33 748.997 -826813 -755.83 -gcp 3014.52 363.68239 -gc812. -1021.34 "D:/Users/(… )/Desktop/image_projetee_UTM.tiff" "C:/Users/(… )/AppData/Local/Temp/image_projetee_UTM.tiff" gdalwarp -r near -order 1 -co COMPRESS=NONE -dstalpha "C:/Users/(… )/AppData/Local/Temp/image_projetee_UTM.tiff" "D:/Users/(… )/Desktop/image_projetee_UTM_georef.tif"

1 - Après avoir exécuté ce code, l'image de sortie n'est pas un GTiff mais un Tiff comme l'image d'entrée.

2 - Si j'essaye d'ouvrir l'image de sortie dans Google Earth par exemple, elle ne semble pas du tout géoréférencée

Je pense que j'ai raté quelque chose dans la façon de géoréférencé un fichier Tiff…

Quelqu'un peut-il m'aider ?


J'ai donc réussi à géoréférencer un fichier tif avec GDAL avec l'utilisation de 4 gcps (points de contrôle au sol). Pour faire cette reprojection, j'utilise gdal en ligne de commande.

Tout d'abord, utilisez gdal translate comme ceci :

gdal_translate -of GTiff -gcp 0 0 -6.848326 45.501053 -gcp 6862 0 -6.490975 45.501503 -gcp 0 1379 -6.762872 45.377363 -gcp 6862 1379 -6.545354 45.382523 "Inputimage.tif" "OutputImage.tif"

Les points de contrôle au sol sont construits comme ceci : coordonnées en pixels dans l'image (x, y) puis localisation géographique (longitude, latitude).

Je suggère de vérifier avec gdalinfo si le fichier de sortie après gdaltranslate a été correctement alimenté avec le gcps. Avec gdalinfo, vous devez voir dans la commande dos quelque chose comme ça :

Le système de coordonnées est "Projection GCP = GCP[ 0] : Id=1, Info= (0,0) -> (45.501053,-6.848326,0) GCP[ 1] : Id=2, Info= (6862,0) - > (45.501503,-6.490975,0) GCP[ 2] : Id=3, Info= (0,1379) -> (45.377363,-6.762872,0) GCP[ 3] : Id=4, Info= (6862,1379 ) -> (45.382523,-6.545354,0)

À ce stade, l'image n'est pas encore géoréférencée. Vous devez utiliser gdalwarp comme ceci (il y a plusieurs façons de reprojeter, ici j'utilise "near" :

gdalwarp -r near -order 1 -co COMPRESS=NONE -dstalpha "OutputImageF_From_gdal_translate.tif" "FinalImage.tif"


Voir la vidéo: Convert between CSV and GeoTIFF with GDAL in Python (Octobre 2021).