Le portail francophone de la géomatique

Bonjour,

Merci pour votre solution du coup je pense utiliser la commande pour faire un VRT.

Par contre mon besoin de fusionner vient de que j'ai 100 dalles dans une projection EPSG: 54004 et si je change la projection cela engendre des bandes noir du coup c'est problématique pour l'utilisation sous Mapinfo.

Je souhaite aprés la fusion reprojeter le VRT puis redécouper le VRT en dalles pour ne pas avoir de bandes noirs
Pouvez-vous m'indiquer la démarche à suivre  et la commande s'il vous plait ?

Je pense que limite un tuto serait très utile à tout le monde.

Cordialement,

Dernière modification par Nosxo (Tue 04 July 2017 12:51)

Bonjour,

Le VRT est un catalogue, c'est un fichier texte qui répertorie le chemin des dalles.
Le passage par le VRT me semble du coup inutile!
En adaptant un peu le script précédent on peut arriver au résultat souhaité:

cd repertoire_origine/
mkdir dalles_L93/
for i in *.jp2;
do gdalwarp -s_srs EPSG:54004 -t_srs EPSG:4326 $i $i_reprojete.jp2;
do mv $i_reprojete.jp2 dalles_L93/;
done;

L'ordi travaillera tout seul pendant un moment et le tour est joué.
Les bandes devraient disparaître aussi, car les dalles seront reprojetées.

Sinon l'usage d'ignmap est probablement + simple, mais je ne connais pas ce programme.
Tiens nous au courant.

Bonne journée.

Bonjour,

Concernant la solution IGNMAP elle ne peut pas fonctionner car je suis de base sur des fichiers en 54004 et IGNMAP ne propose pas ce genre de projection.

A+

Bonjour,
Il semblerait que les vieilles versions de GDAL confondent l'EPSG 54004  et 3857:
https://gis.stackexchange.com/questions … -epsg54004
Peux tu ouvrir tes dalles (dans Qgis par exemple) avec la projection 3857?
Les bandes sont elles encore présentes?

Si elles ont disparues, ça vaut le coup de reprendre la commande gdalwrap avec a_srs 3857 pour tester.
Sinon je sèche...

Re,
Avec ton fichier j'obtiens une erreur due à la projection:

gdalwarp -of GTiff -s_srs EPSG:54004 -t_srs EPSG:2154 Aveyron_1\ 15.map aveyron_repro_pouet.tiffERROR 6: EPSG PCS/GCS code 54004 not found in EPSG support files.  Is this a valid
EPSG coordinate system?
ERROR 1: Translating source or target SRS failed

En fournissant la définition précise du système de coordonnées ça passe!

gdalwarp -of GTiff -s_srs '+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs' -t_srs EPSG:2154 Aveyron_1\ 15.map aveyron_repro_2154.tiff
Creating output file that is 4126P x 4112L.
Processing input file Aveyron_1 15.map.
0...10...20...30...40...50...60...70...80...90...100 - done.

Par contre comme toi j'ai "les bandes noires", voir capture d'écran en PJ.
Une solution (c'est peut-être un peu sale, mais c'est la facilité, car je connais déjà la commande! :p ), c'est de créer un vrt dans lequel tu spécifie que tu veux de la transparence s'il n'y a pas de donnée en entrée!

gdalbuildvrt -srcnodata 0 averyon.vrt *.tiff
0...10...20...30...40...50...60...70...80...90...100 - done.

C'est un peu crade et peut-être que quelqu'un proposera mieux.

À tester en 1er lieu avec plusieurs dalles pour vérifier qu'il n'y ai pas de chevauchement ni d'écart important entre chacune d'elles.
Bonne continuation!

Oui, mais si ignmap ne connait pas la projection 54004 tu ne pourras pas faire de reprojection.

Il semblerait qu'un petit script s'impose du coup...
À adapter:

cd repertoire_origine/
mkdir dalles_L93/
for i in *.map;
do gdalwarp -of GTiff -s_srs '+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs' -t_srs EPSG:2154 $i $i_reprojete.tif;
done;
mv *reprojete* dalles_L93/
mv dalles_L93 ../
gdalbuildvrt -srcnodata 0 mon_super_vrt.vrt ./dalles_L93.tiff

Je t'invite vraiment à lire la doc de GDAL (gdalbuildvrt et gdalwarp dans ton cas), car tu peux faire un tas de trucs sympa.
Dans mon exemple je passe par du .tif mais tu dois pouvoir créer un autre format + avantageux.

Bref à creuser encore je pense!
Tu n'es pas bien loin du but.
À +

Bonjour,

En effet pour palier au problème des bandes noires il faut passer par une fusion des raster au préalable de la projection.

Pour récapituler avec GDAL :
1. Combiner les raster gdalbuildvrt http://www.gdal.org/gdalbuildvrt.html
2. Re-projeter gdalwarp http://www.gdal.org/gdalwarp.html
3. Re-découper gdal_retile.py http://www.gdal.org/gdal_retile.html (sur la base d'un shapefile pour les emprises de découpage)

Martin

Re,

Je suis curieux d'avoir les commandes entrées, c'est toujours utile pour les suivants.
Peux-tu nous les transmettre, stp?

Si gdal te renvoie qu'il ne connaît pas la commande c'est peut-être qu'elle est mal orthographiée
Si tu bosses sous linux, le système est sensible à la casse, du coup Gdal_retile.py n'est pas la même chose que gdal_retile.py.
Dans mon cas c'est la seconde commande qui passe.

Pour ton shapefile, il faut créer un index.
Tu peux utiliser gdal_grid: http://www.gdal.org/gdal_grid.html
Ou encore la fonction grille de Qgis : vecteur > outils de recherche > grille vecteur.

Une fois que ton fichier est créé (dans la même projection que tes rasters, alors tu pourras l'appeler dans ta commande gdal_retile.py.

À +

EDIT: visiblement gdal_grid ne sert pas à ça, il faut peut-être plutôt chercher du côté de gdaltindex. Je suis en limite de compétence, du coup je vais laisser des gens + au courant te répondre! :p

Dernière modification par Mathieu Denat (Wed 05 July 2017 18:38)

Si tu as réalisé la projection avec gdal en EPSG : 2154 tu peux faire ton découpage avec ignmap maintenant

Personnellement je suis sous linux mint et gdal_retile.py fonctionne.

https://gdal.gloobe.org/gdal/gdal_retile.html