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

Annonce

Rencontres QGIS 2025

L'appel à participation est ouvert jusqu'au 19 janvier 2025!

#1 Wed 16 February 2011 18:29

Bertrand DORNER
Participant actif
Lieu: COLMAR
Date d'inscription: 5 Sep 2005
Messages: 90

Lambert 93 ou Lambert 93 ?

Bonsoir,

Derrière ce titre énigmatique, un téléchargement de données depuis le site de l'Inventaire Forestier National. Des fichiers shape, en Lambert 93. En principe, l'importation dans MapInfo doit être sans souci.
Le fichier .prj, associé, contient l'expression PROJCS["RGF93_Lambert_93",GEOGCS["GCS_RGF_1993",DATUM["D_RGF_1993",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",700000.0],PARAMETER["False_Northing",6600000.0],PARAMETER["Central_Meridian",3.0],PARAMETER["Standard_Parallel_1",44.0],PARAMETER["Standard_Parallel_2",49.0],PARAMETER["Latitude_Of_Origin",46.5],UNIT["Meter",1.0]]

L'importation dans MapInfo se fait, en utilisant le fichier.prj fourni. J'affiche les données, les coordonnées sont bien exprimées en Lambert 93. mais, par mauvais (?) esprit, je cherche à vérifier la projection exacte, et j'ai un message me signalant une projection inconnue de ma version de MapInfo (10.0). Dans le mif-mid créé dans la foulée, je lis
Version   300
Charset "WindowsLatin1"
Delimiter ","
CoordSys Earth Projection 3, 104, "m", 3, 46.5, 44, 49, 700000, 6600000 Bounds (-112761336.461, -100805724.411) (114161336.461, 126116948.511)
Columns 6
  AREA Float
  PERIMETER Float
  RNIFN250_ Decimal(11, 0)
  RNIFN250_I Decimal(11, 0)
  REGN Char(3)
  REGIONN Char(60)
Data

Region  1
  482
611685.4 7094947.6
613155 7095683
613909.4 7096011.6

J'arrive néanmoins à reprojeter la couche, proprement, en Lambert 93, bornes Europe, et le mif associé à cette reprojection me donne
Version   300
Charset "WindowsLatin1"
Delimiter ","
CoordSys Earth Projection 3, 33, "m", 3, 46.5, 44, 49.00000000001, 700000, 6600000 Bounds (-792421, 5278231) (3520778, 9741029)
Columns 6
  AREA Float
  PERIMETER Float
  RNIFN250_ Decimal(11, 0)
  RNIFN250_I Decimal(11, 0)
  REGN Char(3)
  REGIONN Char(60)
Data

Region  1
  482
611685.441 7094947.629
613154.991 7095682.971
Les mêmes coordonnées pour les objets, ou très proches, mais pas la même définition de couches

Lorsque, depuis MapInfo, je traduis en shape, et en Lambert 93 cette couche reprojetée, voici l'expression du .prj que j'obtiens
PROJCS["_MI_0",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["Geodetic_Reference_System_of_1980",6378137,298.2572221009113]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["standard_parallel_1",44],PARAMETER["standard_parallel_2",49.00000000001],PARAMETER["latitude_of_origin",46.5],PARAMETER["central_meridian",3],PARAMETER["false_easting",700000],PARAMETER["false_northing",6600000],UNIT["Meter",1]]

Ce n'est pas tout à fait la même chose que le premier .prj

Que s'est-il passé ?
Dans quelle projection, selon le cas, sont représentés les objets ?

Merci d'avance si vous avez des explications

Hors ligne

 

#2 Wed 16 February 2011 21:19

Daniel
Membre
Lieu: Victoriaville, Québec
Date d'inscription: 6 Sep 2005
Messages: 2028

Re: Lambert 93 ou Lambert 93 ?

Salut

Je repondrai ceci, malgré ma "non-connaissance" des projections Française.
Mon expérience personnelle me fait dire de ne jamais utiliser la projection du PRJ (j'enlève le petit crochet) et je fournis moi même la projetction.
Maintenant, est-ce que ça peut faire une différence, je ne sais pas.  Ce que je sais, c'est que parfois, le PRJ créé par le traducteur de Mapinfo n'est pas toujours le même que le PRJ venant directement de ArcGIS (SHP).

Est-ce qu'il y a une solution???  d'autres personnes plus connaissante que moi se feront probablement un plaisir d'y répondre!!!!!

Est-ce que ça vous a aidé???  Pas certain mais ça fait du bien d'en parler!!! smile

Daniel


De retour à la géomatique

Hors ligne

 

#3 Fri 18 February 2011 09:19

Bertrand DORNER
Participant actif
Lieu: COLMAR
Date d'inscription: 5 Sep 2005
Messages: 90

Re: Lambert 93 ou Lambert 93 ?

Bonjour, et merci pour votre réponse.
Le problème s'est posé pour des couches venant de l'Inventaire Forestier National français, mais concerne aussi nos échanges entre les services déconcentrés de l'Etat (dont mon ministère utilise MapInfo, et celui de la santé utilise Arc Gis) et ceux de la collectivité territoriale correspondante (département).
Je vais poser la même question sur le forum ArcGis, mais pour la réception des couches de ma structure, en shape, traduit depuis mapInfo, projetées en Lambert 93, cette réponse http://georezo.net/forum/viewtopic.php? … pe+mapinfo a semble-t-il été utile.
Cordialement,
Bertrand

Hors ligne

 

#4 Fri 18 February 2011 14:15

Robin
GeoRezo forever
Lieu: France
Date d'inscription: 31 Aug 2005
Messages: 13614
Site web

Re: Lambert 93 ou Lambert 93 ?

Bonjour,

Lorsque, depuis MapInfo, je traduis en shape, et en Lambert 93 cette couche reprojetée, voici l'expression du .prj que j'obtiens
PROJCS["_MI_0",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["Geodetic_Reference_System_of_1980",63781


Pour ma part, quand je vois cette projection dans Arcgis, qui est à mon avis un bug pur et dur, je la remplace sans ménagement par la projection théorique (ici Lambert 93) via Arcatalog, un écrasé de projection brut de pomme. Et en général, ça résout le problème.

Hors ligne

 

#5 Fri 18 February 2011 14:59

Bertrand DORNER
Participant actif
Lieu: COLMAR
Date d'inscription: 5 Sep 2005
Messages: 90

Re: Lambert 93 ou Lambert 93 ?

Bonjour,

Cela signifie bien que MapInfo et Arc Gis, sur la projection Lambert 93, ne sont pas vraiment interopérables. Je traduis pourtant mes tables MapInfo, (depuis MapInfo 10.0) en shape, la projection Lambert 93 étant lue dans le .map, via le traducteur universel. Du temps du Lambert zone ou du Lambert 2 Carto, je ne rencontrais pas ce genre de problèmes, pesants dans le cadre d'échanges avec des partenaires.

Ce qui m'intrigue est l'expression "CoordSys Earth Projection 3, 104, "m", 3, 46.5, 44, 49, 700000, 6600000 Bounds (-112761336.461, -100805724.411) (114161336.461, 126116948.511)". En effet, je ne retrouve pas l'expression de cette projection dans le fichier MapInfow.prj que j'utilise, à savoir celui de la version 10, enrichi de quelques descriptions de projections associées au RGF93, et fournies l'an dernier par ce qui était le MEEDDM*.
Quelle différence par rapport à "CoordSys Earth Projection 3, 33, "m", 3, 46.5, 44, 49.00000000001, 700000, 6600000 Bounds (-792421, 5278231) (3520778, 9741029)" ?

Les représentants de MapInfo et ESRI ont-ils été avertis de ces dysfonctionnements ?

* à savoir
"Longitude / Latitude (ETRS89)\p4937",1,115
"Longitude / Latitude (RGF93)\p4965",1,33
"Longitude / Latitude (RGF93 2D)\p4171", 1, 33
et
"--- Universal Transverse Mercator (RGF93) ---"
"UTM Zone 30 (RGF93)",8,33,7,-3,0,0.9996,500000,0
"UTM Zone 31 (RGF93)",8,33,7,3,0,0.9996,500000,0
"UTM Zone 32 (RGF93)",8,33,7,9,0,0.9996,500000,0

Dernière modification par Bertrand DORNER (Fri 18 February 2011 15:43)

Hors ligne

 

#6 Fri 18 February 2011 15:50

Robin
GeoRezo forever
Lieu: France
Date d'inscription: 31 Aug 2005
Messages: 13614
Site web

Re: Lambert 93 ou Lambert 93 ?

On va essayer de reconstituer les élements pour comprendre (voir p 404 du manuel de référence MI - http://reference.mapinfo.com/software/m … rGuide.pdf)

Alors, voyons :

CoordSys Earth Projection 3, 33, "m", 3, 46.5, 44, 49.00000000001, 700000, 6600000


> 3 : c'est le numéro de la projection : le Lambert Conformal Conic
> 33 : Geodetic Reference System 1980 (GRS 80) > ca correspondant à la ligne du prj : SPHEROID["GRS_1980"..
> Le reste sont les autres paramètres

Ensuite :

CoordSys Earth Projection 3, 104, "m", 3, 46.5, 44, 49, 700000, 6600000


> 3 : c'est le numéro de la projection : le Lambert Conformal Conic
> 104, ca correspond au World Geodetic System 1984(WGS 84)
> Le reste sont les autres paramètres

Pourquoi ces différences ? c'est expliqué sur le Wikipédia anglais :
http://en.wikipedia.org/wiki/World_Geod … em:_WGS_84

The WGS 84 originally used the GRS 80 reference ellipsoid, but has undergone some minor refinements in later editions since its initial publication. Most of these refinements are important for high-precision orbital calculations for satellites but have little practical effect on typical topographical uses. The following table lists the primary ellipsoid parameters.
Ellipsoid reference     Semi-major axis a     Semi-minor axis b     Inverse flattening (1/f)
GRS 80     6,378,137.0 m     ≈ 6,356,752.314 140 m     298.257 222 101
WGS 84     6,378,137.0 m     ≈ 6,356,752.314 245 m     298.257 223 563

The very small difference in the flattening thus results in a—very theoretical—difference of 105 µm in the semi polar axis. For most purposes, the differing polar axes can be merged to 6,356,752.3 m, with the inverse flattening rounded to 298.257.


Pour résumer (sauf erreur): WGS84 = ellipsoïde de référence GRS80 + des petites différences.

Encore une coquille qui traine dans les fonds de tiroir...

Dernière modification par Robin (Sat 19 February 2011 14:46)

Hors ligne

 

#7 Fri 18 February 2011 20:40

Maurice
Membre
Lieu: Montpellier
Date d'inscription: 5 Sep 2005
Messages: 5331

Re: Lambert 93 ou Lambert 93 ?

Bonjour
D'accord avec Robin, aux fautes de frappe près... une fois corrigé il faut lire
> 104, ca correspond au World Geodetic System 1984(WGS 84) 
> 33: Geodetic Reference System 1980 (GRS 80)

Si on est tous d'accord, je corrige le post de Robin et supprime celui-ci

Hors ligne

 

#8 Tue 09 July 2013 17:53

Bertrand DORNER
Participant actif
Lieu: COLMAR
Date d'inscription: 5 Sep 2005
Messages: 90

Re: Lambert 93 ou Lambert 93 ?

Bonjour,

Je reviens sur cette discussion, pour signaler une solution qui semble fort pertinente pour sécuriser la transformation de fichiers natifs MapInfo en fichiers shape : ogr2gui, téléchargeable sur http://www.ogr2gui.ca/fr/index.php

Les librairies de base de données (dll) dont il est question sont à placer dans le même répertoire que l'exécutable qui transforme les fichiers, et qui ne fait appel ni à MapInfo, ni à ArcGis.

N'ayant pas ArcGis là où je travaille je n'ai pas pu faire la vérification jusqu'au bout, mais le .prj de l'export en shape obtenu est visiblement correct

Cordialement

Bertrand Dorner

PS : la solution a été donnée par un membre de la liste labo.qgis-owner@texte-a-enlever.developpement-durable.gouv.fr, je ne fais que l'exploiter

Hors ligne

 

Pied de page des forums

Powered by FluxBB