#1 Wed 18 July 2012 10:18
- Frago
- Participant actif
- Date d'inscription: 5 Oct 2011
- Messages: 53
QGIS 1.8 / Gdal : Problème de reprojection IGNF
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,
Hors ligne
#2 Wed 18 July 2012 10:43
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
Bonjour,
Je pense déjà que le "+init=IGNF:LAMB93 +wktext" devrait être "+init=IGNF:LAMB93".
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
#3 Wed 18 July 2012 10:49
- Frago
- Participant actif
- Date d'inscription: 5 Oct 2011
- Messages: 53
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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.
Hors ligne
#4 Wed 18 July 2012 10:56
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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)
Hors ligne
#5 Wed 18 July 2012 11:07
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
entre le lambert 2 étendu epsg, standard code 27572 et l'ignf lamb2E, il y a un gros décalage en effet
geodata au cerema et petits billets en géomatique
Hors ligne
#6 Wed 18 July 2012 11:15
- dominique.lys
- Participant assidu
- Date d'inscription: 5 Oct 2006
- Messages: 473
- Site web
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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.
Hors ligne
#7 Wed 18 July 2012 11:19
- JRM
- Participant assidu
- Lieu: Arras
- Date d'inscription: 15 Apr 2009
- Messages: 521
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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.
Hors ligne
#8 Wed 18 July 2012 11:26
- Frago
- Participant actif
- Date d'inscription: 5 Oct 2011
- Messages: 53
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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 ?
Hors ligne
#9 Wed 18 July 2012 14:10
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: QGIS 1.8 / Gdal : Problème de reprojection 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 ?
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.
Hors ligne
#10 Wed 18 July 2012 17:11
- Frago
- Participant actif
- Date d'inscription: 5 Oct 2011
- Messages: 53
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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)
Hors ligne
#11 Wed 18 July 2012 17:54
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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.
Hors ligne
#12 Thu 19 July 2012 13:52
- Frago
- Participant actif
- Date d'inscription: 5 Oct 2011
- Messages: 53
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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 !
Hors ligne
#13 Fri 20 July 2012 09:08
- Frago
- Participant actif
- Date d'inscription: 5 Oct 2011
- Messages: 53
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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 ?
Hors ligne
#14 Fri 20 July 2012 09:31
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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 ?
Hors ligne
#15 Fri 20 July 2012 10:35
- Frago
- Participant actif
- Date d'inscription: 5 Oct 2011
- Messages: 53
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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 ?
Hors ligne
#16 Fri 20 July 2012 11:15
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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 ?
Hors ligne
#17 Fri 20 July 2012 11:38
- Frago
- Participant actif
- Date d'inscription: 5 Oct 2011
- Messages: 53
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
Oui certain sur gdal_translate (conversion raster vrt -> ecw)
Hors ligne
#18 Fri 20 July 2012 12:05
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
Pouvez-vous donner la commande complète utilisée ?
Hors ligne
#19 Fri 20 July 2012 12:08
- Frago
- Participant actif
- Date d'inscription: 5 Oct 2011
- Messages: 53
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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
Hors ligne
#20 Fri 20 July 2012 13:01
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
Essayez de passer par un fichier TIFF intermédiaire (ce format possède une méthode Create).
Hors ligne
#21 Fri 20 July 2012 14:08
- Frago
- Participant actif
- Date d'inscription: 5 Oct 2011
- Messages: 53
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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 !
Hors ligne
#22 Fri 20 July 2012 18:33
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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.
Hors ligne
#23 Fri 20 July 2012 18:45
- michel56
- Participant assidu
- Lieu: Lorient
- Date d'inscription: 14 Jul 2012
- Messages: 474
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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)
Opérateur cartographique & DAO & PHP.
Hors ligne
#24 Sat 21 July 2012 13:58
- rouault
- Participant assidu
- Date d'inscription: 26 Apr 2009
- Messages: 168
Re: QGIS 1.8 / Gdal : Problème de reprojection IGNF
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
Hors ligne