banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Rencontres QGIS 2025

L'appel à participation est ouvert jusqu'au 19 janvier 2025!

#1 Fri 19 February 2010 15:49

Sam_Dalembert
Participant actif
Lieu: Mérignac
Date d'inscription: 5 Aug 2009
Messages: 121

[GDAL] Problème gdal_translate en ECW avec wktext

Bonjour à tous,


Je viens vous faire part d'un problème assez bizarre. C'est assez simple :

gdalwarp -of VRT -s_srs "+init=IGNF:LAMB2 +wktext" -t_srs "+init=IGNF:LAMB93" fichier_de_depart.ecw fichier_virtuel_reproj.vrt
gdal_translate -of ECW -co TARGET=90 -co LARGE_OK=YES -co PROJ=LMFRAN93 -co DATUM=RGF93 -projwin X Y X Y fichier_virtuel_reproj.vrt Extraction.ecw

Le fichier VRT est bien créé, mais gdal_translate n'arrive pas à créer le fichier ECW, ou plutôt si, il le crée, mais il reste à 1ko, dans la fenêtre Shell ça ne dépassé pas 0, et le processus gdal_translate ne dépasse pas les 15 000 ko (alors que ça monte normalement à 200 000, 500 000, ou même bien plus).
J'ai essayé d'enlever les options PROJ et DATUM, ça ne marche pas, j'ai essayé d'enlever le +wktext, ça fonctionne, mais le résultat est complètement faux.

Par contre, quand j'essaye de sortir le fichier en TIFF au lieu de le sortir en ECW, ça marche. Et ensuite la conversion en ECW fonctionne également.


Quelqu'un a une idée ? sad

Hors ligne

 

#2 Fri 19 February 2010 16:09

sylpingus
Participant occasionnel
Lieu: Aix en Provence
Date d'inscription: 9 Jan 2006
Messages: 34

Re: [GDAL] Problème gdal_translate en ECW avec wktext

Bonjour,

je ne suis pas spécialiste dans le domaine mais j'ai déjà rencontré le problème sous MS Windows en utilisant FWTools. La même manipulation sur Linux compilé avait bien fonctionné. Après, de là à vous dire pourquoi... (mauvaise gestion du VRT sous Windows ?)

Par rapport à l'option +wktext, avez-vous la grille IGNF qui est installée correctement et qui est bien prise en compte ? (cf http://lambert93.ign.fr/index.php?id=28) ? Si le résultat est faux en enlevant ce paramètre, c'est tout à fait normal étant donné qu'il sert à prendre la grille de transformation officielle pour passer du RGF93 au NTF.

Bonne journée

Hors ligne

 

#3 Mon 22 February 2010 08:47

Sam_Dalembert
Participant actif
Lieu: Mérignac
Date d'inscription: 5 Aug 2009
Messages: 121

Re: [GDAL] Problème gdal_translate en ECW avec wktext

Oui désolé, je n'ai pas précisé, mais oui les résultats sont bons (j'ai la grille IGNF) avec wktext.

Je vais essayer de bidouiller le VRT, mais bon...

Hors ligne

 

#4 Mon 22 February 2010 23:18

rouault
Participant assidu
Date d'inscription: 26 Apr 2009
Messages: 169

Re: [GDAL] Problème gdal_translate en ECW avec wktext

Hum, s'il s'agit d'un tout dernier FWTools 2.4.7 qui utilise proj 4.7.0, alors je pense qu'il pourrait s'agir d'une manifestation du problème suivant :
http://trac.osgeo.org/proj/ticket/63

Tu devrais essayer avec une version antérieure pour voir ?

Hors ligne

 

#5 Wed 24 February 2010 09:07

Sam_Dalembert
Participant actif
Lieu: Mérignac
Date d'inscription: 5 Aug 2009
Messages: 121

Re: [GDAL] Problème gdal_translate en ECW avec wktext

Désolé pour le temps de réponse, j'attendais d'avoir un accès à un ordi où les anciennes versions étaient installées.

C'est toi qui a posté le ticket sur trac.osgeo.org/proj/ ? smile

En tous cas merci beaucoup et bien vu tongue C'est exactement ça, j'ai testé avec FWTools 2.4.6, ça ne marchait pas, par contre, avec la version 2.4.2, ça fonctionne. J'ai pris le proj.dll de cette version, je l'ai mis à la place du proj.dll de la version 2.4.7, et ça fonctionne.

Par contre, si ce n'est pas trop demander, tu pourrais m'expliquer ce qu'il faut faire pour résoudre le problème ? Apparemment il y a un patch, mais je n'y comprends rien...:s
Ou est-ce que remplacer le proj.dll suffit ? Ca change quelque chose au niveau de la précision ou pas ?
Je vais continuer les vérifications, mais entre l'ECW converti à partir du TIFF (qui marche avec wktext) et l'ECW que je viens de créer (en remplaçant le dll), je ne vois pas de différence pour l'instant alors que sans le wtext, on voit une différence).

Dernière modification par Sam_Dalembert (Wed 24 February 2010 09:10)

Hors ligne

 

#6 Wed 24 February 2010 09:36

Sam_Dalembert
Participant actif
Lieu: Mérignac
Date d'inscription: 5 Aug 2009
Messages: 121

Re: [GDAL] Problème gdal_translate en ECW avec wktext

Par contre, y a encore un souci là. Les tests marchent avec le Lambert 93, mais quand j'essaye avec une zone CC, il ne reconnait tout simplement le code de la zone !!! (ex: LAMBCC48).

Je viens de tomber sur ce ticket http://trac.osgeo.org/proj/changeset/1440
Ca veut dire qu'on peut rajouter manuellement dans le fichier IGNF ? J'ai l'impression que c'est ça mais bon, je n'y connais pas grand chose donc je préfèrerais avoir l'avis de quelqu'un qui s'y connait justement...


EDIT: bizarre, les codes apparemment ne sont pas sous la forme LAMBCCXX mais RGF93CCXX
http://lambert93.ign.fr/index.php?id=28
Je vais tester si les résultats semblent bons ou non.
Mais dans ce cas, à quel problème correspond le ticket dont je parle ?

Dernière modification par Sam_Dalembert (Wed 24 February 2010 09:52)

Hors ligne

 

#7 Wed 24 February 2010 13:13

rouault
Participant assidu
Date d'inscription: 26 Apr 2009
Messages: 169

Re: [GDAL] Problème gdal_translate en ECW avec wktext

Oui, c'est moi qui avait créé le ticket http://trac.osgeo.org/proj/ticket/63. Ca date d'environ 1 mois.

Plusieurs solutions :
1) Soit tu sais te débrouiller avec un compilateur et tu appliques le patch à proj 4.7.0, ou alors tu prends un snapshot du HEAD du subversion de proj.
2) Tu attends la prochaine release de proj. (Je n'en connais pas le planning)
3) Tu utilises effectivement une version antérieure de proj.dll (4.6.1 par ex). J'ai été un peu surpris que fwtools 2.4.6 ait aussi le problème mais c'est tout à fait possible qu'elle ait déjà proj 4.7.0. Faudrait que tu vérifies en lancant cs2cs et voir le numéro de version qu'il te dit (s'il te dit 4.7.1, faut comprendre 4.7.0 : il y a eu un petit typo apparemment à la génération de 4.7.0). En tout cas utiliser une version antérieure ne devrait pas poser de problème car l'interface de proj est stable (évidemment si tu utilises des fonctionnalités avec des anomalies corrigées dans des versions ultérieures, tu peux avoir des soucis...)

Quant au changeset http://trac.osgeo.org/proj/changeset/1440, il date de 2 ans, donc je pense qu'il ne faut pas se prendre la tête avec du moment que tu n'utilises pas de version proj trop vieille. Par contre, celui qui est intéressant, c'est celui qui suit http://trac.osgeo.org/proj/changeset/1444 où effectivement on voit que les codes LAMBCCXX ont été renommés en RGF93CCXX

Hors ligne

 

#8 Wed 24 February 2010 13:31

Sam_Dalembert
Participant actif
Lieu: Mérignac
Date d'inscription: 5 Aug 2009
Messages: 121

Re: [GDAL] Problème gdal_translate en ECW avec wktext

Merci beaucoup pour ta réponse !

Avec le proj.dll de la version 2.4.7: Rel. 4.7.1, 23 September 2009
Avec le proj.dll de la version 2.4.2: Rel. 4.6.1, 21 August 2008

Par contre, je dois t'avouer qu'appliquer le patch, je ne sais pas du tout comment faire. Et "le HEAD du subversion de proj" c'est du chinois pour moi big_smile

Concrètement, je reprojette des photos aériennes (en Lambert II Centre) en CC48 et en RGF93, je fais des assemblages, des extractions, mais ça s'arrête là. Donc utiliser la version 4.6.1 de proj ne changera rien non ?

Et bon, en faisant des tests visuels (je superpose le cadastre en vecteur), le résultat est impeccable (entre l'image d'origine en L2 et la transfo en CC48 ou L93). Et avec des tests de coordonnées, c'est pareil (je calcule avec Circé, puis avec gdaltransform, et les résultats sont identiques, contrairement aux codes EPSG où il y a une différence avec Circé).
Donc je suppose que c'est bon...

Dernière modification par Sam_Dalembert (Wed 24 February 2010 13:31)

Hors ligne

 

Pied de page des forums

Powered by FluxBB