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é ?

#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,


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

Hors ligne

 

#2 Wed 18 July 2012 10:43

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

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

Frago a écrit:

et avec GDAL (4.7.1) :


La dernière version de GDAL est la 1.9.1 smile

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

MathieuR
Membre
Lieu: aix-en-provence
Date d'inscription: 16 Feb 2009
Messages: 1690
Site web

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

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.

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

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 ?

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

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 ?

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

 

Pied de page des forums

Powered by FluxBB