Le portail francophone de la géomatique

sigdu80 · Thu 28 November 2013 09:54

Bonjour à tous,

je souhaitais convertir des shapefiles L93 en L2 étendu.

J'ai appliqué la commande :

Code:

ogr2ogr -s_srs EPSG:2154 -t_srs EPSG:27572 -a_srs EPSG:27572 -f "ESRI
Shapefile" C:\tests_gdal\Query4_region_l2e\Query4_region_l2e.shp C:\tests_gdal\Q
uery4_region\Query4_region.shp

Cela fonctionne correctement.
Par contre, quand je l'ouvre dans QGIS, et je clique sur Propriétés, concernant cette couche, il me dit : SRC généré

Je m'attendais à part l'identification de l'EPSG:27572 que j'avais désigné.
Est-ce normal ?
Est-ce que j'aurai oublié un paramètre dans l'exécution de ogr2ogr ?

Merci d'avance pour vos précisions.

Nicolas Granier · Fri 29 November 2013 08:30

Bonjour,

je ne réponds pas directement à la question mais vous pouvez attribuer une projection à un shape en passant directement le PRJ dans les paramètres d'OGR. Je m'explique :
- vous identifiez un shape pour lequel la projection est bien lu/reconnu dans Qgis.
- vous faites une copie du fichier PRJ, sous la forme par exemple "prj_lambert93_modele.prj"
- vous passez directement ce PRJ là dans votre transformation OGR dans la commande -a_srs

Code:

ogr2ogr -s_srs EPSG:2154 -t_srs EPSG:27572 -a_srs "c:/mon_dossier/prj_lambert93_modele.prj" -f "ESRI
Shapefile" C:\tests_gdal\Query4_region_l2e\Query4_region_l2e.shp C:\tests_gdal\Q
uery4_region\Query4_region.shp

et là normalement la projection sera reconnue.

A+
Nicolas GRANIER

sigdu80 · Fri 29 November 2013 17:26

Bonjour Nicolas,

si si tu réponds directement à ma question, merci pour ta piste.
je viens d'essayer et en fait, il me met un message d'erreur :

Code:

ERROR 6 : No translation for Lambert_Conformal_Conic to PROJ.4 format is known. ...

C'est dommage car cela aurait résolue le souci.

Alban NOIR · Mon 02 December 2013 15:42

Bonjour,
Ton code d'origine est juste même s'il n'est pas nécessaire de spécifier le -a_srs puisqu'il y a le -t_srs. Avec ogr2ogr il n'y a pas à ma connaissance d'option -dsco pour créer le Query4_region_l2e.prj de la table Query4_region_l2e.shp.
C'est probablement pourquoi si tu tentes d'ouvrir ce SHP qui n'est pas associé à PRJ dans qgis il faut indiquer une projection, ou il va chercher celle du projet ou une des dernières que tu as utilisé.
Alban

Nicolas Granier · Tue 03 December 2013 08:43

Bonjour,
Effectivement même erreur mais on peut la contourner de cette manière.
Tout dabord

Ton code d'origine est juste même s'il n'est pas nécessaire de spécifier le -a_srs puisqu'il y a le -t_srs.

le t_srs et le -a_srs sont deux paramètres différents, le premier s'utilise en combinaison avec -s_srs pour faire une reprojection, alors que le -a_srs permet de définir une projection (comme dans ArcGis avec l'outil "Définir une projection).

Effectivement si tu passes un fichier prj d'Arcgis le formatage ne plait pas à ogr et tu as bien une erreur. Par contre, tu peux contourner le problème de cette façon.
- tu choisis un shape pour lequel la projection est bien lue/reconnue dans Qgis.
- avec ogrinfo tu affiches la projection (option -al), et celui ci te renvoies quelque chose du genre (exemple en L93)

Code:

PROJCS["RGF_1993_Lambert_93",
    GEOGCS["GCS_RGF_1993",
        DATUM["Reseau_Geodesique_Francais_1993",
            SPHEROID["GRS_1980",6378137.0,298.2572221]],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    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]]

- tu copies/colles ceci dans un fichier texte que tu appelles L93_new.prj
- tu peux utiliser ce nouveau PRJ dans le paramètre -a_srs
Je viens de faire le test avec du Lambert 93 et ça fonctionne.

A+

Dernière modification par Nicolas Granier (Tue 03 December 2013 08:56)

sigdu80 · Fri 06 December 2013 17:34

Salut à vous,

Merci.
Petite question : le prj généré permettra de bien détecter la projection du shapefile uniquement dans QGIS ou dans tout logiciel SIG ?

Je pense que l'on tient le bon bout.

Nicolas Granier · Tue 10 December 2013 12:52

le prj généré permettra de bien détecter la projection du shapefile uniquement dans QGIS ou dans tout logiciel SIG ?

Détecté oui, après il faut savoir qu'un Lambert 93 d'Arcgis est différent d'un Lambert93 de Mapinfo et donc probablement de Qgis et que même entre certaines versions de logiciels SIG, on peut avoir des versions de systèmes de projections différents.
A+

Dernière modification par Nicolas Granier (Tue 10 December 2013 12:52)

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 Thu 28 November 2013 09:54

[OGR2OGR] Projection du shapefile de destination

Code:

#2 Fri 29 November 2013 08:30

Re: [OGR2OGR] Projection du shapefile de destination

Code:

#3 Fri 29 November 2013 17:26

Re: [OGR2OGR] Projection du shapefile de destination

Code:

#4 Mon 02 December 2013 15:42

Re: [OGR2OGR] Projection du shapefile de destination

#5 Tue 03 December 2013 08:43

Re: [OGR2OGR] Projection du shapefile de destination

Code:

#6 Fri 06 December 2013 17:34

Re: [OGR2OGR] Projection du shapefile de destination

#7 Tue 10 December 2013 12:52

Re: [OGR2OGR] Projection du shapefile de destination

Pied de page des forums