Le portail francophone de la géomatique

Bonjour René-Luc,
c'est assez bien décrit dans le ticket. Tu ouvres une couche déclarée en 2154 ou IGNF:LAMB93, et tu essaie de changer la projection du projet vers du LAMB2C par exemple. Ou encore tu exporte la couche vers du LAMB2C.
Apparemment c'est proj4 qui n'arrive pas à exploiter la grille ign spécifiée dans les déclarations de projection. (paramètre +nadgrids=ntf_r93.gsb,null). la grille est bien présente dans la distribution dans C:\OSGeo4W\share\proj\ntf_r93.gsb.

Je soupçonne un problème de compilation de proj4 version OSGEO4W. Un test comparatif avec une version compilée hors OSGEO4W serait une bonne aide au diagnostic.
En espérant que ce n'est que ça...

Bonjour,

J'avais vu ce bug il y a quelques temps et il y a eu... les vacances.

En regardant le track de proj.4, il y a eu des modifications il y a 2-4-9 mois dans le moteur géodésique( fichiers geod_* ). Je vais regarder ce soir de chez moi (mais sous linux). Peut-être que temporairement, il faudrait revenir à une version antérieure de proj.4 ?

Cordialement.

Hien.

Bonjour,

J'ai également regardé de mon côté.

C'est un problème dans QGIS : en ligne de commande, la transformation d'un point par cs2cs fonctionne sans problème.

J'ai également fait les test suivant :
- dans les options, en demandant un CRS par défaut pour un nouveau projet de type IGNF:xxxxxxx, la création d'un nouveau projet ne fait pas apparaître de message d'erreur.
- puis la création d'une couche de même projection ne provoque pas d'erreur.

Par contre, chaque action suivante génère une erreur :
- sur un projet, changer son CRS vers un CRS IGNF:xxxxx
- activer la reprojection à la volée

Il semble donc que c'est liée à la transformation de CRS.
En regardant le dépôt Github, le fichier 'src/core/qgscoordinatetransform.cpp' a été modifié, pour permettre la transformation en une fois des points d'un objet. Je ne sais pas encore si c'est lié ou non. Je continue mes recherches.

Bonjour
Pourriez vous nous dire ce que fait le bug exactement ?
De mon côté j'ai récupéré les departements geofla (en lambert93), exporter en lamb2c, reprojeter en lamb2c ou lambert 93 ou lat/long wgs84, il me semble que je n'ai aucun bug (sauf si je suis passé à côté)
Je suis sous debian, avec proj4 4.7.0

Voici les versions que j'ai testé :
- sous windows 7, QGIS 1.8 et QGIS master installés via OSGEO4W sont concernés par l'erreur.
- le message d'erreur à la transformation est le suivant (CRS du projet : IGNF:LAMB1)

Transform error caught: forward transform of
(0.027226, 0.846246)
(0.030423, 0.846282)
(0.033620, 0.846314)
(0.036816, 0.846340)
(0.040014, 0.846362)
(0.043211, 0.846378)
(0.046408, 0.846390)
(0.049605, 0.846397)
(0.027186, 0.847730)
(0.030387, 0.847766)
(0.033589, 0.847797)
(0.036791, 0.847824)
(0.039994, 0.847845)
(0.043196, 0.847862)
(0.046398, 0.847874)
(0.049601, 0.847880)
(0.027145, 0.849213)
(0.030352, 0.849250)
(0.033559, 0.849281)
(0.036766, 0.849308)
(0.039973, 0.849329)
(0.043181, 0.849346)
(0.046389, 0.849357)
(0.049596, 0.849364)
(0.027104, 0.850697)
(0.030316, 0.850733)
(0.033528, 0.850765)
(0.036741, 0.850791)
(0.039953, 0.850813)
(0.043166, 0.850829)
(0.046379, 0.850841)
(0.049592, 0.850848)
(0.027063, 0.852180)
(0.030280, 0.852216)
(0.033497, 0.852248)
(0.036715, 0.852275)
(0.039933, 0.852296)
(0.043151, 0.852313)
(0.046369, 0.852325)
(0.049587, 0.852331)
(0.027021, 0.853663)
(0.030244, 0.853700)
(0.033467, 0.853731)
(0.036690, 0.853758)
(0.039913, 0.853780)
(0.043136, 0.853796)
(0.046359, 0.853808)
(0.049583, 0.853815)
(0.026980, 0.855146)
(0.030208, 0.855183)
(0.033436, 0.855214)
(0.036664, 0.855241)
(0.039892, 0.855263)
(0.043121, 0.855279)
(0.046350, 0.855291)
(0.049578, 0.855298)
(0.026938, 0.856629)
(0.030171, 0.856666)
(0.033405, 0.856697)
(0.036638, 0.856724)
(0.039872, 0.856746)
(0.043106, 0.856763)
(0.046340, 0.856774)
(0.049574, 0.856781)
PROJ.4: +proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs +to +proj=lcc +lat_1=49.5 +lat_0=49.5 +lon_0=0 +k_0=0.99987734 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356514.999978254 +nadgrids=ntf_r93.gsb,null +pm=2.337229167 +units=m +no_defs
Error: failed to load datum shift file

- ce message apparaît, que le projet soit vide ou non (64 points non transformés, quelque soit la taille de la couche) : ce qui fait penser que c'est la grille qui n'arrive pas à être transformée/utilisée ? Les valeurs diffèrent selon l'emprise de la couche à reprojeter.

- la couche est bien chargée (nombre d'entités de la couche cohérent), mais aucune entité n'est affichée.

- ce problème apparaît uniquement sous Windows. Je n'ai pas ce problème sous Linux. (Ubuntu 12.04, QGIS 1.8 et QGIS Master)

Oui, il est bien là, dans %OSGEO4W_ROOT%\share\proj

Confirmé chez moi! je passe l'info aux développeurs natifs. peut-être pas la peine de faire un pull request pour ça.
Merci mille fois!

26 Août 16h28 : appel aux contributeurs
27 Août 14h02 : bug résolu

Pas sûr qu'un SAV "éditeur" soit aussi efficace

Bravo à tous !

Bizarrement il semblerait que le problème soit revenu dans la 2.0 sortie récemment (2.0.1-2, téléchargée hier soir). Rencontrez vous le même problème ?

notes : problème également dans la 1.9 35 mais pas dans la 2.1 beta

w7 64b

Dernière modification par Didier Vanden Berghe (Fri 27 September 2013 10:59)

Est-ce que ce problème vient d'un QGIS packager ou de celui installer via OSGeo4W ?

Un bug avait été ouvert à ce propos :
https://hub.qgis.org/issues/7941

Le problème ne venait pas de QGIS mais d'un script d'OSGeo4W. GRASS modifiait des variables d'environnement ne permettant plus à QGIS de trouver "ntf_r93.gsb".

Salut
installé via l'installateur indépendant 64bits, sur w7

pour ma part le bug était bien présent sur la version installeur indépendant 64 bit.
Le problème vient de ce renouveler avec un collègue de bureau: impossible d'afficher une ortho en lambert 2 étendu, l'image n'apparait pas.
Le bug est vraiment handicapant quand on doit traiter un grand nombre de données pas toujours récente (donc pas en lambert 93).
J'en suis pour le moment revenu à la version 1.8, plus stable à mon goût, en attendant une mise à jour.

Je n'ai pas de windows 64bits sous la main, mais est-il possible d'avoir les variables d'environnements de vos configurations ? Pour cela, allez dans le menu Préférences -> Options -> Système et veuillez m'indiquer les valeurs correspondantes si présentes de OSGEO4W_ROOT, PATH, PROJ_LIB, QGIS_PREFIX_PATH. Je soupçonne une mauvaise config sur une de ces variables.

Ce qui est étrange, c'est que les noms de répertoire soient tous en format court (...~1) alors que même en 32bits, mes noms de répertoires sont au format long (C:\Program Files\...)

Pouvez-vous forcer la variable PROJ_LIB à "C:\Program Files\QGIS Dufour\share\proj" dans la partie haute de la même fenêtre (option écraser) puis redémarrer ?