Le portail francophone de la géomatique

dominique.lys · Mon 06 July 2009 12:44

hello,

Je sais pas si vous avez la même, mais chez moi lorsque je génère un shp avec avec Qgis il ne se cale jamais correctement que ce soit dans Qgis lui même ou dans ArcGIS. En consultant le fichier prj je m'aperçois que le datum n'est pas renseigné :

Code:

 PROJCS["Lambert_Conformal_Conic",GEOGCS["unnamed ellipse",DATUM["D_unknown",SPHEROID["Unknown",6378249.2,293.4660212936269]],PRIMEM["paris",2.3372291666985],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Conformal_Conic"],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]]

du coup je me demande si je ne fais qqchose de travers? une idée?

dominique.lys · Thu 09 July 2009 12:32

Bon je me répond à moi même

Problème: la projection lambert 2 carto des shp créés par QGIS (EPSG 27592) n'est reconnu dans ArcMap

**************************
QGIS : définition EPSG 27572 dans librairie proj4

+proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs

Les pramètres a et b définissent l'ellipsoide de clark1880. Ils peuvent être remplacés par l'écriture +ellps=clrk80

x_0 et y_0 correspondent au décalage Est et Nord (false easting and northing). Le y étant préfixé ou non du numéro de zone

pm indique le méridien d'origine : Greenwitch ou Paris

towgs84 n'est utilisé que pour faire de la reprojection à la volée dans les programme utilisant la librairie proj4

les paramètres lat et long définissent la projection :
long_0 est la longitude d'origine (=0 pour Greenwitch ou =2.33 pour Paris)
lat_0 est latitude d'origine (=46.8)

k correspond au facteur d'échelle

*************************Comparaison fichiers prj

ESRI

PROJCS["NTF_Lambert_II_étendu",
GEOGCS["GCS_NTF",DATUM["D_NTF",SPHEROID["Clarke_1880_IGN",6378249.2,293.46602]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic"],
PARAMETER["False_Easting",600000.0],
PARAMETER["False_Northing",2200000.0],
PARAMETER["Central_Meridian",2.3372291667],
PARAMETER["Standard_Parallel_1",45.8989188889],
PARAMETER["Standard_Parallel_2",47.6960144444],
PARAMETER["Scale_Factor",1.0],
PARAMETER["Latitude_Of_Origin",46.8],
UNIT["Meter",1.0]]

QGIS (27572)

PROJCS["Lambert_Conformal_Conic",
GEOGCS["unnamedellipse",DATUM["D_unknown",SPHEROID["Unknown",6378249.2,293.466021]],
PRIMEM["Unknown",2.337229167],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Lambert_Conformal_Conic"],
PARAMETER["latitude_of_origin",46.8],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",0.99987742],
PARAMETER["false_easting",600000],
PARAMETER["false_northing",2200000],UNIT["Meter",1]]

-------->

Dans le fichier prj définie par QGIS il manque deux paramètres : Standard_Parallel_1 et Standard_Parallel_2 ce sont ces paramètres qu'il faut définir dans la projection proj4 :

+proj=lcc +a=6378249.2 +b=6356515 +lat_0=46.8 +lon_0=2.337229167 +lat_1=45.8989188889 +lat_2=47.6960144444 +x_0=600000 +y_0=2200000, +k_0=1, +towgs84=-168,-60,+320 +units=m +no_defs

lat_1 et lat_2 sont les "Standard_Parallel"

Il suffit donc de créer une projection personnalisée dans Qgis avec ces paramètres pour que le calage soit correct dans ArcGIS

Jaan · Wed 15 July 2009 15:12

Bonjour à tous ceux qui font travailler encore leur neurones après le 15 juillet!

Voila, je suis encore un peu novice mais là je crois avoir tout essayer avant de poster et pas de solution!
J'ai reçu des données au format Shape et je travaille avec Mapinfo 8.0,donc j'utilise le fameux traducteur universel pour traduire mes couches en choisissant bien la bonne projection, ici du lambert II étendu.
Malheureusement, j'ai toujours un décalage d'à peu près 200 km à l'ouest au moment de l'affichage.
J'ai donc essayé toutes les possibilités logiques de projection au cas où mon interlocuteur ait joué avec moi sur la description des couches dans le dossier, mais rien n'y fait je suis dans une impasse!

Je vous met un extrait d'un des fichiers prj d'une des couches reçues:

PROJCS["Lambert II Carto, Méridien de Paris",
GEOGCS["unnamed",DATUM["D_NTF-PM",
SPHEROID["Clarke - 1880",6378249.145,293.465006079117]],
PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic"],
PARAMETER["standard_parallel_1",45.898918964419],
PARAMETER["standard_parallel_2",47.696014502038],
PARAMETER["latitude_of_origin",46.8],
PARAMETER["central_meridian",0],
PARAMETER["false_easting",600000],
PARAMETER["false_northing",2200000],UNIT["METER",1]]

Voilà, ça peut aider à comprendre...
Ce que je trouve pour ma part étrange c'est ce "unnamed" à la ligne 2
Serait-il possible que l'origine de mon problème soit une mauvaise définition de la projection des couches sous arcview ?

Par avance, merci à la communauté!

EDIT:
je viens de me rendre compte que l'ouverture en direct des couches résud mon problème et que quand un décalage se produit sur l'axe est ouest s'est souvent dut à une mauvaise définition du méridien...
Quand pensez-vous? et comment corriger cette erreur sur les fichier shape d'origine?
Merci pour vos réponses...

Dernière modification par Jaan (Wed 15 July 2009 15:47)

Robin · Wed 15 July 2009 16:14

Hello,

Ca c'est bien de la donnée qui vient de Mapinfo, le bug est récurrent. Notez bien l'incohérence :
PROJCS["Lambert II Carto, Méridien de Paris",....PRIMEM["Greenwich"...

Essayez avec ça plutôt :

PROJCS["NTF_Lambert_II_étendu",
GEOGCS["GCS_NTF",DATUM["D_NTF",
SPHEROID["Clarke_1880_IGN",6378249.2,293.46602]],
PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic"],
PARAMETER["False_Easting",600000.0],
PARAMETER["False_Northing",2200000.0],PARAMETER["Central_Meridian",2.3372291667],
PARAMETER["Standard_Parallel_1",45.8989188889],
PARAMETER["Standard_Parallel_2",47.6960144444],
PARAMETER["Scale_Factor",1.0],
PARAMETER["Latitude_Of_Origin",46.8],
UNIT["Meter",1.0]]

Robin · Wed 15 July 2009 16:19

Ps : j"'oubliais, il faut juste remplacer le contenu du fichier PRJ de vos couches

Jaan · Thu 16 July 2009 14:54

Bonjour!

Merci pour la réponse d'hier robin, mais désolé pour la solution!
je me retrouve non plus avec un décalage de 2OO km à l'ouest mais bien 2000 au sud!!!
Je vais essayer de me pencher un peu là dessus...à voir si je m'en sort!
Sinon, toujours preneur de conseils avisés...

Jaan · Thu 16 July 2009 16:06

Re,
bon voila je me suis pas laisser démonter par ce fichu fichier prj!
voila le code qui remet tout à sa place:

Code:

PROJCS["Lambert II Carto, Méridien de Paris",
GEOGCS[[b]"NTF-PM"[/b],DATUM["D_NTF-PM",
SPHEROID["Clarke - 1880",6378249.145,293.4650060791175]],
PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic"],
PARAMETER["standard_parallel_1",45.898918964419],
PARAMETER["standard_parallel_2",47.696014502038],
PARAMETER["latitude_of_origin",46.8],
PARAMETER["central_meridian",[b]2.337229104484[/b]],
PARAMETER["false_easting",600000],
PARAMETER["false_northing",2200000],UNIT["METER",1]]

Voili voilou, en gras les deux corrections apportées à mon fichier de départ!
Pour moi le sujet est résolu mais si quelqu'un veut apporter des précisions sur
la remarque de Robin:

Ca c'est bien de la donnée qui vient de Mapinfo, le bug est récurrent. Notez bien l'incohérence :
PROJCS["Lambert II Carto, Méridien de Paris",....PRIMEM["Greenwich"...

car mon code comporte toujours cette abération mais pourtant fonctionne!
Merci à toute la communauté!

Robin · Fri 17 July 2009 10:57

Hello,

Ca fait des années que je vois ce bug, et pourtant, il me perturbe toujours En fait, dans le principe, même si l'intitulé du début est faux mais les paramètres sont juste, je suppose que c'est pas gênant...

Désolé Jaan si je t'ai égaré... mais merci pour le suivi !
Robin.

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

Annonce

#1 Mon 06 July 2009 12:44

QGIS 1.0 : Projection Lambert 2 datum non défini dans le fichier *.prj

Code:

#2 Thu 09 July 2009 12:32

Re: QGIS 1.0 : Projection Lambert 2 datum non défini dans le fichier *.prj

#3 Wed 15 July 2009 15:12

Re: QGIS 1.0 : Projection Lambert 2 datum non défini dans le fichier *.prj

#4 Wed 15 July 2009 16:14

Re: QGIS 1.0 : Projection Lambert 2 datum non défini dans le fichier *.prj

#5 Wed 15 July 2009 16:19

Re: QGIS 1.0 : Projection Lambert 2 datum non défini dans le fichier *.prj

#6 Thu 16 July 2009 14:54

Re: QGIS 1.0 : Projection Lambert 2 datum non défini dans le fichier *.prj

#7 Thu 16 July 2009 16:06

Re: QGIS 1.0 : Projection Lambert 2 datum non défini dans le fichier *.prj

Code:

#8 Fri 17 July 2009 10:57

Re: QGIS 1.0 : Projection Lambert 2 datum non défini dans le fichier *.prj

Pied de page des forums