Géoréférencement avec GDAL?

16

J'ai une image .tiff que j'essaie de superposer dans google earth. Pour cela, je dois d'abord le géoréférencer. La projection d'image est UTM avec NAD83 comme référence. Je connais le lat / long des quatre coins de l'image.

Je fais de la manière suivante:

  1. gdalinfo pour obtenir les quatre coordonnées.
  2. gdal_translate pour spécifier les informations lat / long
  3. gdal2tiles.py pour carreler l'image.

La seule chose qui me dérange, c'est que je n'ai jamais spécifié UTM et NAD83, alors comment se fait-il qu'il puisse superposer mon image? Quelqu'un sait-il comment spécifier ces informations de projection?

Akd
la source

Réponses:

23

suivez cette voie:

  1. Trouvez votre code de référence spatiale à partir d' ici

  2. Apprenez votre fichier tif coordonnée supérieure gauche et coordonnée inférieure droite.

  3. utilisez cette commande pour la coordonner:

    • gdal_translate -of GTiff -a_ullr ullon ullat lrlon lrlat -a_srs EPSG:4269 input.tif output.tif
  4. utilisez cette commande pour mercator: (epsg: 3857)

    • gdalwarp -of GTiff -t_srs EPSG:3857 input.tif output.tif
  5. puis utilisez gdal2tiles.py pour paver.

Aragon
la source
1
Merci pour la réponse. Mais j'ai besoin de superposer ma projection UTM sur Google Earth qui utilise WGS84 comme référence. Ne devrais-je donc pas définir le système de référence cible sur ESPG: 4326?
Akd
1
@Akd Aragon gdalwarp to EPSG: 3857 est bon pour Google. Voir wiki.openstreetmap.org/wiki/EPSG:3857 pour en savoir plus.
Dave X
1
Selon la documentation et le test: -a_ullr ulx uly lrx lryou avec des termes lat / lon -a_ullr ullon ullat lrlon lrlatest la bonne façon d'utiliser la commande.
zetah
2
J'ai édité le message pour l'ordre correct de x et y dans -a_ullr.
Nikos Alexandris
2
J'essaie également de géoréférencer un tiff, mais je me demande comment pouvons-nous connaître les coordonnées supérieure gauche et inférieure droite? Pourriez-vous m'aider @Aragon?
Manimalis