Nous utilisons des cookies pour vous garantir la meilleure expérience sur notre site. Si vous continuez à utiliser ce dernier, nous considèrerons que vous acceptez l'utilisation des cookies. J'ai compris ! ou En savoir plus !.
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 Mon 24 September 2012 11:28

Lathyrus
Participant occasionnel
Date d'inscription: 23 Jul 2010
Messages: 14

GDAL Coordonées absente fichier ECW

Bonjour et merci d'avance si vous pouvez m'aider.

Voila mon problème je souhaite extraire une dalle ecw d'un fichier vrt issu d'un assemblage de dalles IGN en Lamb93 que j'ai reprojeter en LambE.

j'utilise la commande suivante pour extraire ma dalle.
gdal_translate -of ECW -co TARGET=90 -co LARGE_OK=YES -co DATUM=NTF -co PROJ=LAMBE -a_srs EPSG:27572 -a_nodata 255 assemblage_LAMBE.vrt -projwin 847996 2465180 899707 2416140 Ortho_2009_VdL_LAMBE.ecw
je n'ai pas eu de problème sinon de mémoire mais que j'ai réglé.
j'obtiens un fichier de 3 Go en sortie.
Là ou j'ai un souci c'est que je n'ai pas de coordonnée affiché aux coins de mon image.

Lorsque que je fait un teste sur une surface plus petite mon fichier est nickel.
autre question je pensai que le bord noir de mon image serait devenu blanc et ce n'est pas le cas !

Voici le résultat de gdalinfo:

Driver: ECW/ERMapper Compressed Wavelets
Files: Ortho_2009_VdL_LAMBE.ecw
       Ortho_2009_VdL_LAMBE.ecw.aux.xml
Size is 172370, 163467
Coordinate System is:
PROJCS["NTF (Paris) / Lambert zone II",
    GEOGCS["NTF (Paris)",
        DATUM["Nouvelle_Triangulation_Francaise_Paris",
            SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                AUTHORITY["EPSG","7011"]],
            TOWGS84[-168,-60,320,0,0,0,0],
            AUTHORITY["EPSG","6807"]],
        PRIMEM["Paris",2.33722917,
            AUTHORITY["EPSG","8903"]],
        UNIT["grad",0.01570796326794897,
            AUTHORITY["EPSG","9105"]],
        AUTHORITY["EPSG","4807"]],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    PROJECTION["Lambert_Conformal_Conic_1SP"],
    PARAMETER["latitude_of_origin",52],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",0.99987742],
    PARAMETER["false_easting",600000],
    PARAMETER["false_northing",2200000],
    AUTHORITY["EPSG","27572"],
    AXIS["X",EAST],
    AXIS["Y",NORTH]]
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,163467.0)
Upper Right (172370.0,    0.0)
Lower Right (172370.0,163467.0)
Center      (86185.0,81733.5)
Band 1 Block=172370x1 Type=Byte, ColorInterp=Red
  NoData Value=255
  Overviews: arbitrary
Band 2 Block=172370x1 Type=Byte, ColorInterp=Green
  NoData Value=255
  Overviews: arbitrary
Band 3 Block=172370x1 Type=Byte, ColorInterp=Blue
  NoData Value=255
  Overviews: arbitrary

Cordialement

Hors ligne

 

#2 Mon 24 September 2012 12:36

Nicolas Granier
Participant assidu
Date d'inscription: 19 Apr 2007
Messages: 271

Re: GDAL Coordonées absente fichier ECW

Salut,

-co PROJ=LAMBE


Sauf erreur de ma part LAMBE est le code de projection du registre IGNF, et ici il faut le code de projection d'ER-MAPPER qui est si je me souviens bien "LM2FRANC".

Essayez déjà en remplaçant ce code.
A+

Hors ligne

 

#3 Mon 24 September 2012 13:30

Lathyrus
Participant occasionnel
Date d'inscription: 23 Jul 2010
Messages: 14

Re: GDAL Coordonées absente fichier ECW

Ok
je relance l’opération.

Merci

Hors ligne

 

#4 Mon 24 September 2012 17:32

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

Re: GDAL Coordonées absente fichier ECW

Avec http://trac.osgeo.org/gdal/browser/trun … al_edit.py et GDAL >= 1.9, il n'y a pas besoin de relancer tout le processus de création pour modifier le géoréférencement. On peut modifier sur place pour un certain nombre de formats, en particulier l'ECW :

python gdal_edit.py ton.ecw -mo PROJ=LM2FRANC

et éventuellement -a_ullr avec les valeurs qui vont bien

Hors ligne

 

#5 Mon 24 September 2012 18:29

Lathyrus
Participant occasionnel
Date d'inscription: 23 Jul 2010
Messages: 14

Re: GDAL Coordonées absente fichier ECW

Cool !
Merci
C'est exactement ce que je cherchai

J’essaie demain

Hors ligne

 

#6 Wed 26 September 2012 11:25

Lathyrus
Participant occasionnel
Date d'inscription: 23 Jul 2010
Messages: 14

Re: GDAL Coordonées absente fichier ECW

Salut,

La manip avec le code de projection Er-mapper a marché.
Merci
Mais il faudrait que je fasse la même chose sur un fichier de 10 Go c'est pourquoi ta manip Rouault m'intéresse grandement.

je viens d'essayer ton patch (je ne sais pas si c'est comme ça que cela s’appel ?)
Voici le message que j'ai :

Z:\couches_sig\donnees_de_reference\orthophotographie\2009>python gdal_edit.py O
rtho_2009_SCOT_LAMBE.ecw -mo PROJ=LM2FRANC
Traceback (most recent call last):
  File "gdal_edit.py", line 32, in ?
    from osgeo import gdal
ImportError: No module named osgeo

J'utilise Fwtools 2.4.7
As-tu une idée de ce qui ne marche pas.


Merci d'avance

Hors ligne

 

#7 Wed 26 September 2012 14:31

jcr83
Participant actif
Date d'inscription: 17 Mar 2008
Messages: 147

Re: GDAL Coordonées absente fichier ECW

Even a bien précisé :

GDAL >= 1.9


.
FwTools 1.7.2 est obsolète depuis plusieurs années. Il faut installer GDAL 1.9.1 .

Hors ligne

 

#8 Wed 26 September 2012 16:53

Lathyrus
Participant occasionnel
Date d'inscription: 23 Jul 2010
Messages: 14

Re: GDAL Coordonées absente fichier ECW

Ok !
Mais faut savoir que Fwtools 2.4.7 est compilé avec un Gdal 1.7.2 !
Pas facile a trouver de premier abord.

Bon j'ai installé
python 3.2
GDAL 1.9.1
Les 3 Plugin ECW MrSID et Oracle

Et ...

j'ai toujour le meme message

Z:\couches_sig\donnees_de_reference\orthophotographie\2009>gdal_edit.py Ortho_20
09_SCOT_LAMBE.ecw -mo PROJ=LM2FRANC
Traceback (most recent call last):
  File "Z:\couches_sig\donnees_de_reference\orthophotographie\2009\gdal_edit.py"
, line 32, in <module>
    from osgeo import gdal
ImportError: No module named osgeo

???

Hors ligne

 

#9 Wed 26 September 2012 17:03

Yves
Membre du bureau
Lieu: Aix-les-Bains
Date d'inscription: 22 Mar 2006
Messages: 9869
Site web

Re: GDAL Coordonées absente fichier ECW

Bonjour,

Ce script python doit se baser sur une version pas trop ancienne (genre gdal 1.8) par exemple car il y a eut des évolutions dans la manière d'importer les lib python.

Et ce message d'erreur m'y fait penser :

from osgeo import gdal
ImportError: No module named osgeo


Es tu sur sur que python-gdal a été mis à jour ?

Y.


Yves Jacolin, bénévole de l'association GeoRezo.net, agit au nom et pour le compte de l'association - Partageons ce qui nous départage !!  - GeoRezo vous aide ? Aidez GeoRezo !

Hors ligne

 

#10 Mon 01 October 2012 12:27

Lathyrus
Participant occasionnel
Date d'inscription: 23 Jul 2010
Messages: 14

Re: GDAL Coordonées absente fichier ECW

Bonjour,

Je ne sais pas si les mise a jour sont bien prise en compte.
Pour plus de sureté j'ai désinstallé Fwtool.
mais comment savoir si la liaison Python Gdal fonctionne bien ?

Je suis sous windows Xp.

Merci

Hors ligne

 

#11 Mon 01 October 2012 12:53

jcr83
Participant actif
Date d'inscription: 17 Mar 2008
Messages: 147

Re: GDAL Coordonées absente fichier ECW

Essaie de taper directement la commande d'importation dans la console Python :

Code:

$ python
Python 2.7.3 (default, Sep 25 2012, 12:34:17) 
[GCC 4.5.4] on linux2
Type "help", "copyright", "credits" or "license" for more information.
>>> from osgeo import gdal
>>>

S'il n'y a pas de message d'erreur, comme ci-dessus, c'est OK.

Dernière modification par jcr83 (Mon 01 October 2012 12:53)

Hors ligne

 

#12 Tue 02 October 2012 16:10

Lathyrus
Participant occasionnel
Date d'inscription: 23 Jul 2010
Messages: 14

Re: GDAL Coordonées absente fichier ECW

Merci à vous tous qui m'avez aidé.

la commande gdal_edit.py est génial !

J'ai réglé mon problème de liaison Python Gdal par une réinstalation des différents composants.

1-Python
2-Gdal core
3-Gdal 1.9.1 Py

Peut être que je n'avais pas respecter l'ordre.

Merci encore
Lathyrus

Hors ligne

 

Pied de page des forums

Powered by FluxBB