Le portail francophone de la géomatique

Frago · Wed 18 July 2012 10:18

Bonjour,

Je rencontre un petit soucis de reprojection de fichier .VRT (du Lambert 93 au Lambert II ou Lambert II Etendue)... Le résultat donc un décalage de plusieurs dizaines mètres.

Merci de trouver ci-joints, le fichier VRT L93 original et le résultat de reprojection pas Qgis et Gdal.

J'ai testé la reprojection avec QGis 1.8 :

Code:

gdalwarp -s_srs "+proj=lcc +towgs84=0.0000,0.0000,0.0000 +a=6378137.0000 +rf=298.2572221010000 +lat_0=46.500000000 +lon_0=3.000000000 +lat_1=44.000000000 +lat_2=49.000000000 +x_0=700000.000 +y_0=6600000.000 +units=m +no_defs" -t_srs "+proj=lcc +nadgrids=ntf_r93.gsb,null +a=6378249.2000 +rf=293.4660210000000 +pm=2.337229167 +lat_0=46.800000000 +lon_0=0.000000000 +k_0=0.99987742 +lat_1=46.800000000 +x_0=600000.000 +y_0=200000.000 +units=m +no_defs" -of VRT "l93.vrt" "l2.vrt"

et avec GDAL (4.7.1) :

Code:

gdalwarp -of VRT -rc -s_srs "+init=IGNF:LAMB93 +wktext" -t_srs "+init=IGNF:LAMB2" l93.vrt l2.vrt

Aussi, le gdalinfo ne me donne pas le même résultat avec les deux fichiers :

Code:

Driver: VRT/Virtual Raster
Files: [b]l2_qgis.vrt[/b]
       l93.vrt
Size is 22871, 22871
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["unnamed ellipse",
        DATUM["unknown",
            SPHEROID["unnamed",6378249.2,293.466021],
            EXTENSION["PROJ4_GRIDS","ntf_r93.gsb,null"]],
        PRIMEM["unnamed",2.337229167],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_1SP"],
    PARAMETER["latitude_of_origin",46.8],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",0.99987742],
    PARAMETER["false_easting",600000],
    PARAMETER["false_northing",200000],
    UNIT["Meter",1]]
Origin = (532318.881982434310000,346534.555956897500000)
Pixel Size = (2.647762553580680,-2.647762553580680)
Corner Coordinates:
Upper Left  (  532318.882,  346534.556) (  0d54'31.70"W, 48d 6'52.30"N)
Lower Left  (  532318.882,  285977.579) (  0d53'58.20"W, 47d34'11.81"N)
Upper Right (  592875.859,  346534.556) (  0d 5'44.40"W, 48d 7'4.82"N)
Lower Right (  592875.859,  285977.579) (  0d 5'40.87"W, 47d34'24.21"N)
Center      (  562597.371,  316256.067) (  0d29'58.79"W, 47d50'40.88"N)

Code:

Driver: VRT/Virtual Raster
Files: [b]l2_gdal.vrt[/b]
       l93.vrt
Size is 22871, 22870
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["unnamed ellipse",
        DATUM["unknown",
            SPHEROID["unnamed",6378249.2,293.466021]],
        PRIMEM["unnamed",2.337229167],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_1SP"],
    PARAMETER["latitude_of_origin",46.8],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",0.99987742],
    PARAMETER["false_easting",600000],
    PARAMETER["false_northing",200000],
    UNIT["Meter",1]]
Origin = (532263.659772898190000,346527.532997859990000)
Pixel Size = (2.647813999941495,-2.647813999941495)
Corner Coordinates:
Upper Left  (  532263.660,  346527.533) (  0d54'34.37"W, 48d 6'52.05"N)
Lower Left  (  532263.660,  285972.027) (  0d54'0.84"W, 47d34'11.61"N)
Upper Right (  592821.814,  346527.533) (  0d 5'47.01"W, 48d 7'4.59"N)
Lower Right (  592821.814,  285972.027) (  0d 5'43.45"W, 47d34'24.03"N)
Center      (  562542.737,  316249.780) (  0d30'1.42"W, 47d50'40.66"N)

Auriez-vous déjà eu le même problème ?
Pensez vous que cette anomalie provient de mes fichiers de configuration des proj ?

-> A noter également que je viens de reprojecter en CC48 sans décalage..

Merci d'avance pour votre aide,

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.

Yves · Wed 18 July 2012 10:43

Bonjour,

Je pense déjà que le "+init=IGNF:LAMB93 +wktext" devrait être "+init=IGNF:LAMB93".

Y.

Frago · Wed 18 July 2012 10:49

Bonjour Yves,

Désolé, j'ai oublié de mentionner que j'avais déjà essayer sans le "wktext", mais le résultat est exactement le même.

jcr83 · Wed 18 July 2012 10:56

Frago a écrit:

et avec GDAL (4.7.1) :

La dernière version de GDAL est la 1.9.1

J'ai testé les trois chaînes avec PROJ4, et je n'ai pas vu de différence :

$ cs2cs +init=IGNF:LAMB93 +to +init=IGNF:LAMB2
583000 6780000
532324.94 346025.95 0.00

$ cs2cs +init=IGNF:LAMB93 +wktext +to +init=IGNF:LAMB2
583000 6780000
532324.94 346025.95 0.00

$ cs2cs +proj=lcc +towgs84=0.0000,0.0000,0.0000 +a=6378137.0000 +rf=298.2572221010000 +lat_0=46.500000000 +lon_0=3.000000000 +lat_1=44.000000000 +lat_2=49.000000000 +x_0=700000.000 +y_0=6600000.000 +units=m +no_defs +to +proj=lcc +nadgrids=ntf_r93.gsb,null +a=6378249.2000 +rf=293.4660210000000 +pm=2.337229167 +lat_0=46.800000000 +lon_0=0.000000000 +k_0=0.99987742 +lat_1=46.800000000 +x_0=600000.000 +y_0=200000.000 +units=m +no_defs
583000 6780000
532324.94 346025.95 -0.00

Dernière modification par jcr83 (Wed 18 July 2012 11:07)

MathieuR · Wed 18 July 2012 11:07

entre le lambert 2 étendu epsg, standard code 27572 et l'ignf lamb2E, il y a un gros décalage en effet

dominique.lys · Wed 18 July 2012 11:15

Bonjour,

Quelle est la données de référence qui te permet de constater le décalage; est-elle en Lambert 2 ou en Lambert 93? En tout cas, si elle est en RGF93 c'est peut-être simplement lié à la reprojection à la volée.

JRM · Wed 18 July 2012 11:19

Je ne sais plus si j'avais eu cette conversation ici ou sur forumsig, je résume : le décalage métrique depuis une projection RGF93/Lambert93 vers une projection NTF en utilisant les codes IGNF -> détail connu de l'IGN et considéré comme normal

Qu'il n'y ai pas de problème entre un LAMB93 et un CC48 est normal puisque le second est juste une excroissance du premier et également basé sur le RGF93. Pas besoin d'utiliser l'appel +wktext lorsqu'on convertit entre ces 2 systèmes.

Frago · Wed 18 July 2012 11:26

Je constate une différente en fonction de versions N- x années du fichier à reprojeter et en fonction d'autres référentiels (ortho et fichiers vecteurs). La reprojection à la volée est désactivée.

Baobazz : je passe souhaite pas passer par EPSG mais reprojeter exclusivement via IGNF

JCR83 : j'obtiens 4.7.1 en tapant cs2cs dans ma fenêtre shell sous Windows. J'ai Qgis 1.8 Lisboa d'installer. Comment voir si j'obtiens les mêmes résultats que toi en fonction de ma configuration ?

jcr83 · Wed 18 July 2012 14:10

Frago a écrit:

JCR83 : j'obtiens 4.7.1 en tapant cs2cs dans ma fenêtre shell sous Windows. J'ai Qgis 1.8 Lisboa d'installer. Comment voir si j'obtiens les mêmes résultats que toi en fonction de ma configuration ?

cs2cs ne fait pas partie de GDAL, mais de PROJ.4. J'ai la version 4.8.0.
Je ne suis pas sûr de comprendre ta question : il suffit de taper les commandes citées plus haut.

Frago · Wed 18 July 2012 17:11

En fait, je ne sais pas comment faire pour connaître ma version de GDAL...

Quant aux trois chaines à tester, j'ai -0.00 aux deux premières. Est ce que cela a une incidence ?

Complément d'information :
si avec PROJ4 je reprojecte la coordonnée du centre de mon raster, je n'obtiens pas les mêmes coordonnées que celles des VRT reprojetés avec GDAL (avec ou sans wktext)

Dernière modification par Frago (Wed 18 July 2012 17:19)

jcr83 · Wed 18 July 2012 17:54

1 / gdalinfo --version
2 / Il me semble que -0 est égal à 0, non ?
3/ C'est normal, lorsqu'on reprojette un raster, l'étendue du résultat n'est pas exactement la même que l'étendue du raster initial (à cause de la rotation induite par la reprojection). Il faut utiliser un point de référence (une borne IGN par exemple) pour effectuer la comparaison.

Frago · Thu 19 July 2012 13:52

1/ je suis sous 1.7.0b2 via FwTools.

--> question : étant donné que QGIS intégre GDAL 1.9.1, est-il possible de "faire pointer" FwTools vers cette version à jour ?

2/ il m'est déjà arrivé d'avoir le nord au sud et inversement. Je me demandais si ça pouvait venir de là.

3/ effectivement, je n'y avais pas pensé.

... je continu mes tests de reprojection, avec ou sans QGIS, avec ou sans wktext. Si vous voyez un moyen de régler cette anomalie, je suis preneur !

Frago · Fri 20 July 2012 09:08

Bonjour,

J'ai installé QGIS via OSGeo4W et passe par le shell où j'ai bien la version 1.9.1. Par contre, je ne parviens pas à rajouter le format ECW (malgré avoir copié les dll ERDAS dans le bin).

Je reprojete mes vrt avec OsGeo4W.
Pour générer l'ECW je reviens donc dans Fwtools toujours sous 1.7.0b2 et là décalage ! Par contre, si je génère un Gtiff pas de décalage (!! c'est à devenir fou).

Une idée sur le problème ?

jcr83 · Fri 20 July 2012 09:31

Frago a écrit:

J'ai installé QGIS via OSGeo4W et passe par le shell où j'ai bien la version 1.9.1. Par contre, je ne parviens pas à rajouter le format ECW (malgré avoir copié les dll ERDAS dans le bin).

Les variables d'environnement GDAL_DATA et GDAL_DRIVER_PATH pointent-elles sur les bons dossiers ? Si oui, la commande "gdalinfo --formats" liste-t-elle le format ECW ?

Frago · Fri 20 July 2012 10:35

Dans mon --formats, j'ai bien :

Code:

  ECW (rov): ERDAS Compressed Wavelets (SDK 4.2)
  JP2ECW (rov): ERDAS JPEG2000 (SDK 4.2)

gdal_translate -of ECW me renvoi

Code:

Input file size is 7627, 7627
0ERROR 6: GDALDriver::Create() ... no create method implemented for this format.

... par contre dans les variables d'environnement je n'ai aucun ligne concernant GDAL. Dois-je les rajouter et pointer vers mon C:\OSGeo4W\bin ?

jcr83 · Fri 20 July 2012 11:15

Frago a écrit:

gdal_translate -of ECW me renvoie
Code:
ERROR 6: GDALDriver::Create() ... no create method implemented for this format.

Vous êtes certain que c'est gdal_translate et pas gdalwarp qui renvoie cette erreur ?

Frago · Fri 20 July 2012 11:38

Oui certain sur gdal_translate (conversion raster vrt -> ecw)

jcr83 · Fri 20 July 2012 12:05

Pouvez-vous donner la commande complète utilisée ?

Frago · Fri 20 July 2012 12:08

J'ai toujours utilisé la commande suivante :

Code:

gdal_translate -of ECW -co TARGET=90 -co LARGE_OK=YES -co PROJ=LM1FRE2D -co DATUM=NTF l2.vrt l2.ecw --config GDAL_ECW_CACHE_MAXMEM 130000000

jcr83 · Fri 20 July 2012 13:01

Essayez de passer par un fichier TIFF intermédiaire (ce format possède une méthode Create).

Frago · Fri 20 July 2012 14:08

Oui, j'arrive à générer toute sorte d'image en Gtiff... mais il me faut du ECW j'ai pas des Go à allouer à ces données !

jcr83 · Fri 20 July 2012 18:33

J'ai bien écrit "intermédiaire". Dès que vous avez le Gtiff, convertissez-le en ECW. Je pense que votre problème est dû au fait que vous essayez de convertir directement un VRT en ECW.

michel56 · Fri 20 July 2012 18:45

Bonjour,

http://ignmap.ign.fr/

Il permet de reprojeter les tuiles et les redécoupage également.

Fonctionne sous Windows. Sous linux Linux?

Voir:
http://georezo.net/forum/viewtopic.php?id=80233

Dernière modification par michel56 (Fri 20 July 2012 18:45)

rouault · Sat 21 July 2012 13:58

A noter que OSGeo4W utilise le SDK ECW 4.2 RO, lecture seule. Il n'est donc pas possible de créer des ECW avec. On peut utiliser un des binaires disponibles à http://vbkto.dyndns.org/sdk/ qui utilise l'ancien SDK 3.3 qui permet la création d'ECW

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

#1 Wed 18 July 2012 10:18

QGIS 1.8 / Gdal : Problème de reprojection IGNF

Code:

Code:

Code:

Code:

#2 Wed 18 July 2012 10:43

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#3 Wed 18 July 2012 10:49

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#4 Wed 18 July 2012 10:56

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#5 Wed 18 July 2012 11:07

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#6 Wed 18 July 2012 11:15

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#7 Wed 18 July 2012 11:19

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#8 Wed 18 July 2012 11:26

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#9 Wed 18 July 2012 14:10

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#10 Wed 18 July 2012 17:11

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#11 Wed 18 July 2012 17:54

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#12 Thu 19 July 2012 13:52

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#13 Fri 20 July 2012 09:08

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#14 Fri 20 July 2012 09:31

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#15 Fri 20 July 2012 10:35

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

Code:

Code:

#16 Fri 20 July 2012 11:15

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

Code:

#17 Fri 20 July 2012 11:38

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#18 Fri 20 July 2012 12:05

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#19 Fri 20 July 2012 12:08

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

Code:

#20 Fri 20 July 2012 13:01

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#21 Fri 20 July 2012 14:08

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#22 Fri 20 July 2012 18:33

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#23 Fri 20 July 2012 18:45

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

#24 Sat 21 July 2012 13:58

Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF

Pied de page des forums