Nous utilisons des cookies pour vous garantir la meilleure expérience sur notre site. Si vous continuez à utiliser ce dernier, nous considèrerons que vous acceptez l'utilisation des cookies. J'ai compris ! ou En savoir plus !.
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

Printemps des cartes 2024

#1 Wed 09 July 2008 15:03

VIGNE Céline
Juste Inscrit !
Date d'inscription: 6 May 2008
Messages: 3

projection WGS84 et Lambert 2 etendu

Je débute dans  MapServer.
J’ai deux couches à projeter: l’une en WGS84 et l’autre en Lambert 2 étendu
Les couches se superposent parfaitement dans MapInfo mais dans MapServer, la couche en WGS84 se retrouve nettement plus basse que l’autre.
Les diverses recherches au sein du forum n’ont pas permis de corriger le décalage.
Je vous remercie d’avance pour toute aide.

Voici le mapfile

MAP
  NAME                PNG_
    EXTENT            0 1000000 1000000 3000000
  IMAGETYPE     png
  SIZE                     1840 2000
  SHAPEPATH     "../data"
  IMAGECOLOR 240 240 240
   
  DEBUG ON

  OUTPUTFORMAT
    NAME              png
    DRIVER            "GD/PNG"
    MIMETYPE     "image/png"
    IMAGEMODE                PC256
    EXTENSION    "png"
  END

  WEB
    log                     "../tmp/sig.log"
    TEMPLATE      "../sig.html"
    IMAGEPATH  "../tmp/"
    IMAGEURL     "http://localhost:8070/sig/"
  END

    PROJECTION
    "init=epsg:27562"
  END
 
  # Start of LAYER DEFINITIONS ---------------------------------------------

  LAYER # DEPARTEMENTS
   NAME         departements
   DATA         lim_adm_france_dom
   STATUS       DEFAULT
    TYPE         LINE

    CLASS
      COLOR        0 0 0
    END
  END #

  LAYER
    NAME                              "communes"
    STATUS                           DEFAULT
    TYPE                  POINT
    CONNECTIONTYPE     OGR 
    CONNECTION                               "../ressources/geoannuaire.ovf"

    PROJECTION
"proj=latlong"
                "ellps=WGS84"
                "datum=WGS84"
                #             "lat_0=0"
                #             "lon_0=0"
                #             "x_0=600000"
                #             "y_0=1200000"
            #  "init=epsg:4326"
                UNIT "Decimal_Degree"
    END
               
    CLASS
      NAME            "pointage_ville"
      STYLE

      SYMBOL        'circle'
        SIZE               15
        COLOR         255 200 0
      END
    END
  END

  # End of LAYER DEFINITIONS -------------------------------
END # end of map file/object

Hors ligne

 

#2 Wed 09 July 2008 20:38

VIGNE Céline
Juste Inscrit !
Date d'inscription: 6 May 2008
Messages: 3

Re: projection WGS84 et Lambert 2 etendu

Une simple erreur de epsg.
J'ai changé 27582 au lieu de 27562.
J’ai obtenu la superposition, à ceci près que les projections ne sont pas sur le même méridien (GREENWICH et PARIS).
Une petite addition de 2° 20’ 14.025" a fini de régler mon problème (source IGN).

Hors ligne

 

#3 Wed 09 July 2008 20:43

ThomasG
Membre
Lieu: Nantes
Date d'inscription: 9 Sep 2005
Messages: 943
Site web

Re: projection WGS84 et Lambert 2 etendu

Bonjour,

Sauf erreur le code EPSG lambert 2 étendu (bien que le libellé dise le contraire) est 27572 et pas 27562 qui lui correspond au lambert 2 centre.

Je me permet cette affirmation suite au sujet http://georezo.net/forum/viewtopic.php?id=44073

Cela semble se confirmer avec les définitions EPSG

Code:

PostGIS SRID:27562
+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

et

Code:

PostGIS SRID:27572
+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

Tu constates que la valeur y_0 a une valeur qui varie de 2 million qui vraisemblablement doit correspondre à ton décalage vertical et que les autres variables sont similaires.

La solution est donc de changer ton code EPSG dans le mapfile.

Cordialement

Thomas G

PS:Merci de prévenir si ça marche.

Hors ligne

 

#4 Wed 09 July 2008 20:47

ThomasG
Membre
Lieu: Nantes
Date d'inscription: 9 Sep 2005
Messages: 943
Site web

Re: projection WGS84 et Lambert 2 etendu

Rebonjour,

Tu dis utiliser 27582 comme code EPSG, il faut utiliser le 27572: l'autre est dit "deprecated" (obsolète)
Le lien dans mon post précédent explique aussi ça.

Cordialement

Thomas G

Hors ligne

 

#5 Fri 11 July 2008 16:02

VIGNE Céline
Juste Inscrit !
Date d'inscription: 6 May 2008
Messages: 3

Re: projection WGS84 et Lambert 2 etendu

Merci beaucoup pour l'explication
Je corrige 27582 par 27572

cordialement

Hors ligne

 

Pied de page des forums

Powered by FluxBB