Géoréférencement Raster ECW

[Sylvain M.]

Bonjour à tous,

J'essaye de me mettre petit à petit à QGis pour avoir une alternative à nos 2 licences ArcGis saturées dans notre équipe.
Comme tous nos raster (Bd Ortho) sont en ECW, je dois absolument trouver une solution pour les lire correctement dans QGis.
J'ai bien réussi à ouvrir ces raster grâce au post de Odurr, et la qualité d'affichage est parfaite.
Cependant, j'ai beau définir le SCR (Système de Référence de Coordonnées) de la couche (en l’occurrence Lambert 93), celle-di ne s'affiche pas correctement (son coin supérieur gauche prend les coordonnées X et Y = 0)
Je suppose que les fichiers de géoréférencement (tab par exemple) fournis par l'IGN ne sont pas lus par QGis.

Avez-vous rencontré cette difficulté, et comment puis la contourner ?

Merci d'avance.

Sylvain M.

[Barbapapa]

Salut Sylvain,
Pourtant si tu arrives à les ouvrir dans ArcGis, ils devraient s'ouvrir sans problèmes dans Qgis.
Tu n'as pas de fichier .eww ou .tab avec par hasard ?

[Sylvain M.]

J'arrive effectivement à les ouvrir correctement dans ArcGis.
Ces ecw sont accompagnés d'un .TAB, d'un .AUX.XML (généré par ArcGis), et d'un .GRF, donc pas de .EWW

[Barbapapa]

Essaies en laissant uniquement les fichiers ecw tous seuls (ts les autres fichiers dans un sous dossier provisoirement).

[Sylvain M.]

Effectivement,

J'ai vérifié, le pb ne vient pas des .TAB, mais des .AUX.XML
Ceux-ci sont générés par ArcGis au moment de la définition du système de coordonnées de mes raster.
En effet, ArcGis reconnait les ECW de l'IGN comme étant en "Lambert Conformal Conic" :

Lambert_Conformal_Conic
false_easting: 700000.000000
false_northing: 6600000.000000
central_meridian: 3.000000
standard_parallel_1: 44.000000
standard_parallel_2: 49.000000
latitude_of_origin: 46.500000
Unité linéaire: Meter

FRANCE 1993
Datum: D_RGF93

alors que mes données SIG en sont généralement en "RGF93 Lambert 93" :

RGF93_Lambert_93
Projection:  Lambert_Conformal_Conic
False_Easting: 700000.000000
False_Northing: 6600000.000000
Central_Meridian: 3.000000
Standard_Parallel_1: 44.000000
Standard_Parallel_2: 49.000000
Latitude_Of_Origin: 46.500000
Unité linéaire: Meter

GCS_RGF_1993
Datum: D_RGF_1993

Les paramètres semblent identiques, mais comme le nom est différent, je préfère redéfinir le système de projection des données IGN, ce qui me créé ces fichiers .AUX.XML avec les ECW...

Maintenant, si je les supprime, QGis lit bien mes raster, mais ArcGis les considère dans un système différent de celui de mes données...

Que me conseillez vous ?

<édition>
Pour compléter, voici le contenu d'un de mes fichiers .AUX.XML :

<PAMDataset>
  <SRS>PROJCS[&quot;RGF93_Lambert_93&quot;,GEOGCS[&quot;GCS_RGF_1993&quot;,DATUM[&quot;RGF_1993&quot;,SPHEROID[&quot;GRS_1980&quot;,6378137.0,298.257222101]],PRIMEM[&quot;Greenwich&quot;,0.0],UNIT[&quot;Degree&quot;,0.0174532925199433]],PROJECTION[&quot;Lambert_Conformal_Conic_2SP&quot;],PARAMETER[&quot;False_Easting&quot;,700000.0],PARAMETER[&quot;False_Northing&quot;,6600000.0],PARAMETER[&quot;Central_Meridian&quot;,3.0],PARAMETER[&quot;Standard_Parallel_1&quot;,44.0],PARAMETER[&quot;Standard_Parallel_2&quot;,49.0],PARAMETER[&quot;Latitude_Of_Origin&quot;,46.5],UNIT[&quot;Meter&quot;,1.0]]</SRS>
</PAMDataset>

</édition>

[Sylvain M.]

Bonjour à tous,

Utilisateur habituel d'ArcGis, j'essaie d'explorer un peu les potentialités de QGis, que je délaissais jusqu'à présent pour la complexité d'intégration des raster ECW (quasi la totalité de mes raster).
Effectivement, la dernière version (1.7.4) lit très bien les ECW.

Malheureusement, je rencontre un problème pour le géoréférencement de ceux-ci : mes rasters se chargent bien, mais leur positionnement est faux.
Leur coin supérieur gauche se place aux coordonnées [0;0], bien que le SCR des couches et du projet aient été définis (RGF93 / Lambert-93 ; EPSG 2154).

En fait, je ne sais pas quel fichier analyse QGis pour le géoréférencement des ECW : est-ce les données de l'entête ECW ? (j'ai essayé de regarder dedans avec ECW Header Editor, mais ça semble correct), ou est-ce dans un XML ou autre ?

Merci pour votre aide !

Sylvain M.

[Sylvain M.]

Je me permet de remonter la discussion... au cas où !

[jcr83]

Bonjour,
QGIS passe par GDAL pour lire les fichiers ECW. Que donne la commande gdalinfo sur ton fichier ?

[Sylvain M.]

Merci de te pencher sur mon problème.
Excuse moi, mais ne connaissant pas (encore) GDAL ni QGIS, comment puis-je saisir cette commande GDAL info ? Cela se passe-t-il dans QGis ou dois-je lancer une autre application ?
(désolé si la question est idiote !)

A+

Sylvain
*******
EDIT
*******
Bonjour,
J'ai trouvé comment utiliser "gdalinfo" avec FWTOOLS (peut-être était-ce possible directement dans QGis ?), et voila ce que ça donne avec un de mes ECW (BD-Ortho IGN) :

C:\Program Files\FWTools2.4.7>gdalinfo E:\SIG\DONNEES\REFERENTIELS\IGN\BDOrtho\D
EPT\61-2001-L2E\61-2001-0360-2395-LA2E-C10.ecw
Driver: ECW/ERMapper Compressed Wavelets
Files: E:\SIG\DONNEES\REFERENTIELS\IGN\BDOrtho\DEPT\61-2001-L2E\61-2001-0360-239
5-LA2E-C10.ecw
       E:\SIG\DONNEES\REFERENTIELS\IGN\BDOrtho\DEPT\61-2001-L2E\61-2001-0360-239
5-LA2E-C10.ecw.aux.xml
Size is 10000, 10000
Coordinate System is:
PROJCS["NTF_France_II_degrees",
    GEOGCS["GCS_NTF",
        DATUM["Nouvelle_Triangulation_Francaise",
            SPHEROID["Clarke_1880_IGN",6378249.2,293.46602]],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_1SP"],
    PARAMETER["False_Easting",600000.0],
    PARAMETER["False_Northing",2200000.0],
    PARAMETER["Central_Meridian",2.337229166666667],
    PARAMETER["Standard_Parallel_1",46.8],
    PARAMETER["Scale_Factor",0.99987742],
    PARAMETER["Latitude_Of_Origin",46.8],
    UNIT["Meter",1.0]]
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,10000.0)
Upper Right (10000.0,    0.0)
Lower Right (10000.0,10000.0)
Center      ( 5000.0, 5000.0)
Band 1 Block=10000x1 Type=Byte, ColorInterp=Red
  Overviews: arbitrary
Band 2 Block=10000x1 Type=Byte, ColorInterp=Green
  Overviews: arbitrary
Band 3 Block=10000x1 Type=Byte, ColorInterp=Blue
  Overviews: arbitrary

Et effectivement, le coin supérieur gauche a pour coordonnées 0;0... Maintenant, je ne comprends pas trop le reste et où se trouve l'erreur !

Quelqu'un pour m'aider ?

Sylvain M.

[Lud]

Peux tu rééssayer en copiant le fichier ecw (et juste celui là sans le ecw.aux.xml) ailleurs dans un répertoire ?
Je note par contre que tu dis avoir défini dans qgis la couche en lambert93 alors que là gdal semble t'indiquer qu'elle est en lambert 2 étendue, à voir donc. (mais ça n'explique pas le 0,0)

[Sylvain M.]

Effectivement, j'ai pour le test "gdalinfo" utilisé un raster en Lambert 2 étendu (trompé de dossier).
J'ai justement aussi essayer de charger mes raster dans QGis sans le fichier XML, et là, pas de souci. Le problème vient donc de mes XML (générés par ArcGis).

Voici par exemple le contenu d'un XML de raster en Lambert 93 :

<PAMDataset>
  <SRS>PROJCS[&quot;RGF93_Lambert_93&quot;,GEOGCS[&quot;GCS_RGF_1993&quot;,DATUM[&quot;RGF_1993&quot;,SPHEROID[&quot;GRS_1980&quot;,6378137.0,298.257222101]],PRIMEM[&quot;Greenwich&quot;,0.0],UNIT[&quot;Degree&quot;,0.0174532925199433]],PROJECTION[&quot;Lambert_Conformal_Conic_2SP&quot;],PARAMETER[&quot;False_Easting&quot;,700000.0],PARAMETER[&quot;False_Northing&quot;,6600000.0],PARAMETER[&quot;Central_Meridian&quot;,3.0],PARAMETER[&quot;Standard_Parallel_1&quot;,44.0],PARAMETER[&quot;Standard_Parallel_2&quot;,49.0],PARAMETER[&quot;Latitude_Of_Origin&quot;,46.5],UNIT[&quot;Meter&quot;,1.0]]</SRS>
</PAMDataset>

Est-ce que vous voyez quelque chose qui cloche là dedans ?

Sylvain M.

[Lud]

en fait il est possible qu'il n'y ait pas de pb dans le aux.xml mais que ça trouble gdal et qgis et qu'il cherche les coordonnées de l'image dedans alors qu'elles se trouvent dans l'en tête de l'ecw.

[Sylvain M.]

Je me rend compte que je perds la mémoire... J'avais déjà posté pour le meme sujet, avec une version précédente de QGis...
Mais pas de solution à l'époque non plus...

Sujet à cette adresse....

Finalement, le problème reste identique d'une version à l'autre de QGis...

[Lud]

Merci pour le lien, j'ai fusionné les discussions.

[jcr83]

Bonjour,
Le problème vient probablement de GDAL. La solution est de ne garder que les fichiers ECW et de supprimer les XML.
Tu peux aussi envoyer un rapport de bug au projet GDAL : http://trac.osgeo.org/gdal/