[Gfoss] problema con g.proj

Bud_P_Bruegger · January 10, 2008, 1:50pm

Provo di convertire la descrizione della projezione in formato proj.4
in wkt per creare un SRID in postgis. Ma mi da un segmentation fault.

postgres@gis:~$ cat cs-gr-trans2.proj4
+proj=cass +lat_0=43.3174975 +lon_0=11.33244638 +x_0=-2.0 +y_0=-16.0
+ellps=intl +units=m +no_defs
postgres@gis:~$ cat cs-gr-trans2.proj4 |/usr/lib/grass/bin/g.proj -w
-proj4=- Segmentation fault

C'e' qualcuno che mi sa dire che spaglio?

mille grazie
-b

--
Bud P. Bruegger, Ph.D. +39-0564-488577 (voice), -21139 (fax)
   European Chair, Global Collaboration Forum on eID
   Chair, Porvoo Subgroup on collab. govs/operating systems
   Leader of the Permanent eID Status Observatory (PESO) project
Servizio Elaborazione Dati e-mail: bud@comune.grosseto.it
Comune di Grosseto jabber: bud@jabber.no
Via Ginori, 43 http://www.comune.grosseto.it/
58100 Grosseto (Tuscany, Italy)
http://www.comune.grosseto.it/interopEID/

neteler · January 10, 2008, 5:26pm

On Jan 10, 2008 2:50 PM, Bud P. Bruegger <bud@comune.grosseto.it> wrote:

Provo di convertire la descrizione della projezione in formato proj.4
in wkt per creare un SRID in postgis. Ma mi da un segmentation fault.

postgres@gis:~$ cat cs-gr-trans2.proj4
+proj=cass +lat_0=43.3174975 +lon_0=11.33244638 +x_0=-2.0 +y_0=-16.0
+ellps=intl +units=m +no_defs
postgres@gis:~$ cat cs-gr-trans2.proj4 |/usr/lib/grass/bin/g.proj -w
-proj4=- Segmentation fault

C'e' qualcuno che mi sa dire che spaglio?

Dovresti usare m.proj
http://grass.osgeo.org/grass63/manuals/html63_user/m.proj.html

ciao
markus

neteler · January 11, 2008, 5:20pm

2008/1/11 Bud P. Bruegger <bud@comune.grosseto.it>:

grazie mille!

On Thu, 10 Jan 2008 18:26:32 +0100
"Markus Neteler" <neteler@osgeo.org> wrote:

> On Jan 10, 2008 2:50 PM, Bud P. Bruegger <bud@comune.grosseto.it> wrote:
> > Provo di convertire la descrizione della projezione in formato proj.4
> > in wkt per creare un SRID in postgis. Ma mi da un segmentation fault.
> >
> > postgres@gis:~$ cat cs-gr-trans2.proj4
> > +proj=cass +lat_0=43.3174975 +lon_0=11.33244638 +x_0=-2.0 +y_0=-16.0
> > +ellps=intl +units=m +no_defs
> > postgres@gis:~$ cat cs-gr-trans2.proj4 |/usr/lib/grass/bin/g.proj -w
> > -proj4=- Segmentation fault
> >
> > C'e' qualcuno che mi sa dire che spaglio?
>
> Dovresti usare m.proj
> http://grass.osgeo.org/grass63/manuals/html63_user/m.proj.html

Scusa, ovviamente la mia risposta era nonsense.

cat cs-gr-trans2.proj4 | g.proj -w
PROJCS["Lambert Conformal Conic",
    GEOGCS["grs80",
        DATUM["North_American_Datum_1983",
            SPHEROID["Geodetic_Reference_System_1980",6378137,298.257222101]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",36.16666666666666],
    PARAMETER["standard_parallel_2",34.33333333333334],
    PARAMETER["latitude_of_origin",33.75],
    PARAMETER["central_meridian",-79],
    PARAMETER["false_easting",609601.22],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]

Funziona bene...

Il trucco è:
- non chiamare GRASS moduli direttamente con il path, ma
- chiamare g.proj nella GRASS shell senza path.

ciao
Markus

neteler · January 11, 2008, 5:46pm

On Jan 11, 2008 6:20 PM, Markus Neteler <neteler@osgeo.org> wrote:

2008/1/11 Bud P. Bruegger <bud@comune.grosseto.it>:
> grazie mille!
>
> On Thu, 10 Jan 2008 18:26:32 +0100
> "Markus Neteler" <neteler@osgeo.org> wrote:
>
> > On Jan 10, 2008 2:50 PM, Bud P. Bruegger <bud@comune.grosseto.it> wrote:
> > > Provo di convertire la descrizione della projezione in formato proj.4
> > > in wkt per creare un SRID in postgis. Ma mi da un segmentation fault.
> > >
> > > postgres@gis:~$ cat cs-gr-trans2.proj4
> > > +proj=cass +lat_0=43.3174975 +lon_0=11.33244638 +x_0=-2.0 +y_0=-16.0
> > > +ellps=intl +units=m +no_defs
> > > postgres@gis:~$ cat cs-gr-trans2.proj4 |/usr/lib/grass/bin/g.proj -w
> > > -proj4=- Segmentation fault
> > >
> > > C'e' qualcuno che mi sa dire che spaglio?
> >
> > Dovresti usare m.proj
> > http://grass.osgeo.org/grass63/manuals/html63_user/m.proj.html

Scusa, ovviamente la mia risposta era nonsense.

cat cs-gr-trans2.proj4 | g.proj -w
PROJCS["Lambert Conformal Conic",
    GEOGCS["grs80",
        DATUM["North_American_Datum_1983",
            SPHEROID["Geodetic_Reference_System_1980",6378137,298.257222101]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",36.16666666666666],
    PARAMETER["standard_parallel_2",34.33333333333334],
    PARAMETER["latitude_of_origin",33.75],
    PARAMETER["central_meridian",-79],
    PARAMETER["false_easting",609601.22],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]

Funziona bene...

... da punto di vista tecnico. Solo che il tuo script ha
dato una risposta diversa:

postgres@gis:~/py$ ./proj2wkt.py
PROJCS["unnamed",
   GEOGCS["International 1909 (Hayford)",
       DATUM["unknown",
           SPHEROID["intl",6378388,297]],
       PRIMEM["Greenwich",0],
       UNIT["degree",0.0174532925199433]],
   PROJECTION["Cassini_Soldner"],
   PARAMETER["latitude_of_origin",43.3174975],
   PARAMETER["central_meridian",11.33244638],
   PARAMETER["false_easting",-2],
   PARAMETER["false_northing",-16],
   UNIT["Meter",1]]

Cerco a controllarlo.

Markus

neteler · January 13, 2008, 1:13pm

On Jan 11, 2008 6:46 PM, Markus Neteler <neteler@osgeo.org> wrote:

On Jan 11, 2008 6:20 PM, Markus Neteler <neteler@osgeo.org> wrote:
> 2008/1/11 Bud P. Bruegger <bud@comune.grosseto.it>:
> > On Thu, 10 Jan 2008 18:26:32 +0100
> > "Markus Neteler" <neteler@osgeo.org> wrote:
> > > On Jan 10, 2008 2:50 PM, Bud P. Bruegger <bud@comune.grosseto.it> wrote:
> > > > Provo di convertire la descrizione della projezione in formato proj.4
> > > > in wkt per creare un SRID in postgis. Ma mi da un segmentation fault.
> > > >
> > > > postgres@gis:~$ cat cs-gr-trans2.proj4
> > > > +proj=cass +lat_0=43.3174975 +lon_0=11.33244638 +x_0=-2.0 +y_0=-16.0
> > > > +ellps=intl +units=m +no_defs
> > > > postgres@gis:~$ cat cs-gr-trans2.proj4 |/usr/lib/grass/bin/g.proj -w
> > > > -proj4=- Segmentation fault
> > > >
> > > > C'e' qualcuno che mi sa dire che spaglio?

...

postgres@gis:~/py$ ./proj2wkt.py
PROJCS["unnamed",
   GEOGCS["International 1909 (Hayford)",
       DATUM["unknown",
           SPHEROID["intl",6378388,297]],
       PRIMEM["Greenwich",0],
       UNIT["degree",0.0174532925199433]],
   PROJECTION["Cassini_Soldner"],
   PARAMETER["latitude_of_origin",43.3174975],
   PARAMETER["central_meridian",11.33244638],
   PARAMETER["false_easting",-2],
   PARAMETER["false_northing",-16],
   UNIT["Meter",1]]

OK, risolto (mancava 'proj4=-' come parametro):

cat cs-gr-trans2.proj4 | g.proj -w proj4=-
WARNING: Datum <unknown> not recognised by GRASS and no parameters found
PROJCS["Cassini",
    GEOGCS["international",
        DATUM["unknown",
            SPHEROID["International_1924",6378388,297]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Cassini_Soldner"],
    PARAMETER["latitude_of_origin",43.3174975],
    PARAMETER["central_meridian",11.33244638],
    PARAMETER["false_easting",-2],
    PARAMETER["false_northing",-16],
    UNIT["Meter",1]]

ciao
Markus