Comment ajouter une projection à un Shapefile à l'aide d'outils open-source?

On m'a donné des fichiers de formes SANS projection. J'ai donc besoin d'ajouter la projection afin de l'afficher sur google map. J'ai trouvé les informations de projection dont j'ai besoin sur http://spatialreference.org/ref/epsg/2066 .

Question:
Comment puis-je ajouter la projection en utilisant une méthode open source (je n'ai pas accès au logiciel ESRI)?

Merci
Chris

Utilisez ogr2ogr de GDAL :

ogr2ogr -a_srs EPSG:2066 output.shp input.shp

Ce qui est bien avec ogr, c'est qu'il peut également utiliser un fichier PRJ existant en entrée:

ogr2ogr -a_srs input.prj output.shp input.shp

La projection dans le fichier .shp se trouve dans shapefilename.prj. Il s'agit essentiellement d'un fichier texte, avec des informations que vous avez trouvées. Jetez un œil à certains de vos fichiers de formes et découvrez comment les informations sont écrites. Ou moins compliqué, si vous avez un autre fichier .shp dans la même projection, copiez-le et changez le nom en conséquence.

Créez un fichier .prj avec le même nom de base que votre fichier de formes (c'est-à-dire foo.prj), ouvrez le fichier dans un éditeur de texte et ajoutez le contenu de, dans ce cas, http://spatialreference.org/ref/epsg/2066/ prj /

Téléchargez le fichier .prj sur le lien que vous avez publié, puis renommez-le pour qu'il corresponde à votre fichier de forme.

Trouvez ci-joint un petit script qui écrit la référence de projection de spatialreference.org dans un fichier .prj. Il ajoute un fichier de projection à tous les fichiers spécifiés dans un répertoire. Par exemple, tous les fichiers de formes dans le répertoire E: \ .Je vous inquiétez seulement du code EPSG de la projection que vous souhaitez intégrer, de l'extension des fichiers auxquels vous souhaitez ajouter un fichier de projection et du répertoire où ils se trouvent. Il parcourra récursivement tous les sous-répertoires, donc utilisez-le avec précaution.

    import os

    def getWKT_PRJ (epsg_code):
     import urllib.request, urllib.parse, urllib.error

    # Access projection information
     wkt = urllib.request.urlopen("http://spatialreference.org/ref/epsg/{0}/prettywkt/".format(epsg_code))
     decoded = (wkt.read().decode('utf-8')) 

    # Remove spaces between charachters
     remove_spaces = decoded.replace(" ","")        
    # Place all the text on one line
     output = remove_spaces.replace("\n","")
     return output

    def referencer(folder_path, extension):
        for path, subdirs, files in os.walk(folder_path):
            for name in files:
                file_extension = os.path.splitext(name)[-1]
                if(extension in file_extension):
                    file_path = os.path.join(path,name)
                    file_name = os.path.splitext(file_path)[0]
                    prj = file_name + ".prj"
                    projection = open(prj,"w")
                    projection.write(epsg)
                    projection.close()

    epsg = getWKT_PRJ("25831")              
    referencer('E:\Testfolder', '.shp')

Essayez d'utiliser Geokettle http://www.spatialytics.org/ il lit les fichiers de formes et a une option pour ajouter un SRS dans ses transformations. Sinon, utilisez ogr2ogr gui (windows uniquement) pour voir le srid