Le portail francophone de la géomatique

Bonjour,

Merci pour ta réponse Lionel!!
Pour les coordonnées géographiques suivantes :
latitude = 48.816667 et longitude = 2.266667
J'obtiens des coordonnées en x et y mais en lambert 93 je crois.
x = 5715525.74450435 et y = 168079.951955856
et avec la variable init :
x = 2434149.98545833 et y = 938103.167756713

Alors que je les souhaiterais en Lambert 2 étendu.
Donc je rappelle à l'aide.


GH

Re

Je l'ai déja utilisé, voici mon code :

# To Convert from Lat/Long to LAMBERT:
  #!/usr/bin/perl

  use strict;
  use Geo::Proj4;

  my $proj = Geo::Proj4->new( proj => "lcc", ellps => "clrk80", init => "epsg:27572");
# , datum => NTF (Paris));
 
  my $lat = 48.816667;
  my $long = 2.266667;
 
  my ($x, $y) = $proj->forward($lat, $long);
  print "conversion to LAMBERT II : y is  $y\n";
  print "conversion to LAMBERT II : x is  $x\n";

  my ($lat, $long) = $proj->inverse($x, $y);
  print "inverse conversion: lat is $lat \n" ;
  print "inverse conversion: long is $long \n" ;

J'obtiens ce résultat avec :
x = 2434149.98545833 et y = 938103.167756713

C'est sur que je me rapproche du résultat voulu mais le "y" est un peu énorme.
A moins que je l'écris mal.
Des propositions??

GH

Merci pour le raccourci, c'est vrai quand on regarde bien le fichier epsg ça donne toutes les infos.
Et concernant, le résultat obtenu, qu'en pensez vous?
Y a i il une erreur quelque part???

Re,

Le point correspond à la ville "issy les moulineaux".
Je pense que c'est  la derniere : proj4  version 4.4.9.

Entre temps j'ai mis la derniere version 4.5.0.     
Les résultats que j'obtiens sont differents:
y is  2426646.42005569
x is  766532.73567244

G.H.

Dernière modification par gheloise (Mon 11 June 2007 16:27)

Re,

Ca c'est sur!!
Vous pensez que l'erreur vient de mon code ou d'ailleurs?
C'est frustrant d'être si pres du but.


G.H.

Re,

j'ai trouvé ça mais je n'arrive pas à l'intégrer à mon code

my $from  = Geo::Proj4->new("+proj=latlong +datum=NAD83");
my $to    = Geo::Proj4->new("+proj=utm +zone=10 +datum=WGS84");

my $point = [ 1.12, 3.25 ];  # See Geo::Point
my $pr_point = $from->transform($to, $point);

G.H.

my $pr    = $from->transform($to, [ $point1, $point2 ]);
my $pr_point1 = $pr->[0];
my $pr_point2 = $pr->[1];

Bonjour,

2.266667°E et 48.816667°N (WGS84)
Pour contrôle, les coordonnées Lambert II étendu retournées par Circé 2000 (IGN) pour ces coordonnées géographiques :
N(m) = 2424254 et E(m) = 594868
le long de la D130 à Issy-les-Moulineaux.

Ceci dit, dans "Proj4text" de "spatial_ref_sys" (srid=27572, voir ci-dessous), au passage, lon_0 = 0 et pm=paris???

Les coordonnées géographiques sont tjrs exprimées en référence au méridien de Greenwich, non?

proj4txt pour srid =27572 dans spatial_ref_sys :
+proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs


Cordialement

Dernière modification par ph_de_fr (Thu 14 June 2007 14:45)

Bonjour,

concernant proj4, vous n'obtiendrez  pas les mêmes coordonnées que CIRCE étant donné que CIRCE utilise une grille "made in france" différente des grilles utilisées par proj4 (l'une est basée sur un décalage des coordonnées cartésiennes, l'autre sur un décalage des coordonnées géographiques) en conséquence vous ne pourrez et ceci en récupérant les bons facteurs (ellipsoide et projection) qu'avoir la précision d'environ 2 mètres.

Pour information, vous devriez pouvoir trouver les bons paramètres sur le WEB : il faut redefinir l'elipsoide car en principe la lat et long sont en rapport à l'ellipsoide WGS84 et le lambert2 est en rapport avec l'eelipsoide Clarke80 (qui est un Clarke80 modifié par l'IGN donc vous ne trouverez les bons paramètres que chez eux -l'IGN). Donc vous ne pouvez utiliser proj mais devez utiliser cs2cs

Je sais également qu'un simili de grille existe et qu'il est possible de constituer une grille type NAD à partir de la grille IGN mais le système des grilles sur proj pose un gros problème en cela qu'il existe une multitude de version et que la grille n'est jamais adapté à la version que l'on a. (Bref, il est nécessaire de retoucher le code si on veut vraiment travailler avec des grilles)

Si vous ne trouvez pas les bons paramètres, "email me" !

(proj4 est tres bien et tres puissant mais comme beaucoup de libre, il n'est pas mis en valeur par une bonne documentation francaise)

PROJ4 est une librairie C faisant de la transfo de SRS, il se peut qu'on module perl existe (http://devpit.org/wiki/Perl/Geo-Proj4)

Sinon, GDAL/OGR; la lib SIG de reference, utilise proj4 et peut se piloter en python, c, java (sur) et peut etre en perl:
http://fedora.kiewel-online.de/extras/6 … perl.html.

Nicolas

Bonjour,

  je suis intéressé moi aussi par scripter cette transformation, avez vous pu avancer sur cela ?

Emmanuel