[Gfoss] stringa proj

Local Chapters GFOSS.it
1

Ciao a tutti.
ho uno shapefile con un prj contenente questa stringa.

PROJCS["ED_1950_Lambert_Conformal_Conic",GEOGCS["GCS_European_1950",DATUM["D_European_1950",SPHEROID["International_1924",6378388.0,297.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",13.0],PARAMETER["Standard_Parallel_1",41.91666666666666],PARAMETER["Standard_Parallel_2",42.08333333333334],PARAMETER["Latitude_Of_Origin",42.0],UNIT["Meter",1.0]]

il dato riguarda il territorio italiano.

qualcuno sa darmi qualche dritta su come tradurla in proj4?

si tratta di una lambert conforme conica (lcc) secante su due paralleli
praticamente attaccati.

in particolare non ho idea di come trattare i due paralleli in questione
(standard parallel1 e 2)...

magari qualcuno di voi ha già affrontato e risolto il problema...

grazie..

ciao

--
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

2

ivan marchesini ha scritto:

Ciao a tutti.
ho uno shapefile con un prj contenente questa stringa.

PROJCS["ED_1950_Lambert_Conformal_Conic",GEOGCS["GCS_European_1950",DATUM["D_European_1950",SPHEROID["International_1924",6378388.0,297.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",13.0],PARAMETER["Standard_Parallel_1",41.91666666666666],PARAMETER["Standard_Parallel_2",42.08333333333334],PARAMETER["Latitude_Of_Origin",42.0],UNIT["Meter",1.0]]

il dato riguarda il territorio italiano.

qualcuno sa darmi qualche dritta su come tradurla in proj4?

con testepsg la validazione fallisce...
Qui c'e' lo schema da seguire:
http://www.remotesensing.org/geotiff/proj_list/lambert_conic_conformal_2sp.html

si tratta di una lambert conforme conica (lcc) secante su due paralleli
praticamente attaccati.

in particolare non ho idea di come trattare i due paralleli in questione
(standard parallel1 e 2)...

+lat_0 Latitude of origin
+lat_1 Latitude of first standard parallel
+lat_2 Latitude of second standard parallel
Fonte: http://trac.osgeo.org/proj/wiki/GenParms

magari qualcuno di voi ha già affrontato e risolto il problema...

In definitiva, se non erro, in formato proj.4 dovrebbe essere:

+proj=lcc +ellps=intl +x_0=0 +y_0=0 +lat_0=42 +lon_0=13
+lat_1=41.91666666666666 +lat_2=42.08333333333334 +units=m +no_defs

Facci sapere se va a buon fine.

ciao
Antonio

--
Antonio Falciano
http://www.linkedin.com/in/antoniofalciano

3

Grazie mille Antonio...
oggi non ho fatto in tempo a provarci ma appena ci metto mano vi faccio
sapere...
una sola domanda: essendo il datum ed50, non servono i parametri
towgs84? per definirlo?

ciao grazie

Il giorno gio, 04/06/2009 alle 10.17 +0200, Antonio Falciano ha scritto:

ivan marchesini ha scritto:
> Ciao a tutti.
> ho uno shapefile con un prj contenente questa stringa.
>
> PROJCS["ED_1950_Lambert_Conformal_Conic",GEOGCS["GCS_European_1950",DATUM["D_European_1950",SPHEROID["International_1924",6378388.0,297.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",13.0],PARAMETER["Standard_Parallel_1",41.91666666666666],PARAMETER["Standard_Parallel_2",42.08333333333334],PARAMETER["Latitude_Of_Origin",42.0],UNIT["Meter",1.0]]
>
> il dato riguarda il territorio italiano.
>
> qualcuno sa darmi qualche dritta su come tradurla in proj4?

con testepsg la validazione fallisce...
Qui c'e' lo schema da seguire:
http://www.remotesensing.org/geotiff/proj_list/lambert_conic_conformal_2sp.html

> si tratta di una lambert conforme conica (lcc) secante su due paralleli
> praticamente attaccati.
>
> in particolare non ho idea di come trattare i due paralleli in questione
> (standard parallel1 e 2)...

+lat_0 Latitude of origin
+lat_1 Latitude of first standard parallel
+lat_2 Latitude of second standard parallel
Fonte: http://trac.osgeo.org/proj/wiki/GenParms

> magari qualcuno di voi ha già affrontato e risolto il problema...

In definitiva, se non erro, in formato proj.4 dovrebbe essere:

+proj=lcc +ellps=intl +x_0=0 +y_0=0 +lat_0=42 +lon_0=13
+lat_1=41.91666666666666 +lat_2=42.08333333333334 +units=m +no_defs

Facci sapere se va a buon fine.

ciao
Antonio

4

ivan marchesini ha scritto:

Grazie mille Antonio...
oggi non ho fatto in tempo a provarci ma appena ci metto mano vi faccio
sapere...
una sola domanda: essendo il datum ed50, non servono i parametri
towgs84? per definirlo?

Certo, la stringa proj.4 che ho abbozzato era di primo tentativo e non
avendo dati sotto mano ero curioso di sapere con che grado di
errore/approssimazione avveniva la trasformazione.
In alternativa al +towgs84, si potrebbe utilizzare anche un grigliato
per la trasformazione verso WGS84 (+nadgrids=...). Il tutto naturalmente
dipende dalla precisione che ti occorre.

ciao

--
Antonio Falciano
http://www.linkedin.com/in/antoniofalciano

5

Ottimo antonio...
grazie...
appena riesco a metterci le mani vi faccio sapere...
ciao

Il giorno gio, 04/06/2009 alle 19.22 +0200, Antonio Falciano ha scritto:

ivan marchesini ha scritto:
> Grazie mille Antonio...
> oggi non ho fatto in tempo a provarci ma appena ci metto mano vi faccio
> sapere...
> una sola domanda: essendo il datum ed50, non servono i parametri
> towgs84? per definirlo?

Certo, la stringa proj.4 che ho abbozzato era di primo tentativo e non
avendo dati sotto mano ero curioso di sapere con che grado di
errore/approssimazione avveniva la trasformazione.
In alternativa al +towgs84, si potrebbe utilizzare anche un grigliato
per la trasformazione verso WGS84 (+nadgrids=...). Il tutto naturalmente
dipende dalla precisione che ti occorre.

ciao

--
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a
06125
Perugia (Italy)
Socio fondatore GFOSS "Geospatial Free and Open Source Software" http://www.gfoss.it
e-mail: marchesini@unipg.it
        ivan.marchesini@gmail.com
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: geoivan73@jabber.org