prj2epsg

[Filippo]

Penso che già conoscevate

http://prj2epsg.org/search

non sempre infallibile....ma a volte e a qualcuno potrebbe essere utile...:-)

[Paolo Corti]

> Penso che già conoscevate
>
> http://prj2epsg.org/search
>
> non sempre infallibile....ma a volte e a qualcuno potrebbe essere
> utile...:-)
>

Usando GDAL (testato sulla 1.8) e senza usare un servizio on line.
Crea il file esriprj2standards ed inserisci le seguenti linee:

import sys
from osgeo import osr

def esriprj2standards(shapeprj_path):
    prj_file = open(shapeprj_path, 'r')
    prj_txt = prj_file.read()
    srs = osr.SpatialReference()
    srs.ImportFromESRI([prj_txt])
    print 'Shape prj is: %s' % prj_txt
    print 'WKT is: %s' % srs.ExportToWkt()
    print 'Proj4 is: %s' % srs.ExportToProj4()
    srs.AutoIdentifyEPSG()
    print 'EPSG is: %s' % srs.GetAuthorityCode(None)

esriprj2standards(sys.argv[1])

a questo punto puoi eseguire il programmino, dandogli in pasto un prj:

$ python esriprj2standards.py /home/pcorti/data/shapefile/country.prj
Shape prj is: GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]
WKT is: GEOGCS["GCS_WGS_1984",DATUM["WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]
Proj4 is: +proj=longlat +datum=WGS84 +no_defs
EPSG is: 4326

Divertiti a modificarlo come meglio credi, buon divertimento

ciao
P

--
Paolo Corti
Geospatial software developer
web: http://www.paolocorti.net
twitter: @capooti
_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
Non inviate messaggi commerciali.
I messaggi di questa lista non rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.
502 iscritti all'11.2.2011

[Filippo]

Grazie Paolo,
proverò...;-)

[pcav]

Grazie Paolo.
Puoi metterlo sul wiki, cosi' che non rimanga seppellito nelgi archivi
della ML?
Grazie.

Il giorno mer, 06/04/2011 alle 16.37 +0200, Paolo Corti ha scritto:

> Usando GDAL (testato sulla 1.8) e senza usare un servizio on line.
> Crea il file esriprj2standards ed inserisci le seguenti linee:
>
> import sys

--
http://www.faunalia.it/pc

_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
Non inviate messaggi commerciali.
I messaggi di questa lista non rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.
502 iscritti all'11.2.2011