[Gfoss] proj4text per 3003 parte peninsulare

moovida · July 24, 2007, 12:01pm

Ciao, ho un dubbio riguardo a come fare ad avere in postgis una gauss
boaga con l'informazione del towgs per le trasformazioni. Non sono
riuscito a trovare tale info in una delle proiezioni.
Quindi parto dalla proiezione identificata da epsg 3003.

Mi basta aggiungere a srtext e proj4text i towgs ad entrambi, vero?

PROJCS["Monte Mario / Italy zone 1 / Pensinsular optimized",
   GEOGCS["Monte Mario",
     DATUM["Monte Mario",
       SPHEROID["International 1924", 6378388.0, 297.0,
AUTHORITY["EPSG","7022"]],
      TOWGS84[104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68]
       AUTHORITY["EPSG","6265"]],
     PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
     UNIT["degree", 0.017453292519943295],
     AXIS["Geodetic longitude", EAST],
     AXIS["Geodetic latitude", NORTH],
     AUTHORITY["EPSG","4265"]],
   PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]],
   PARAMETER["central_meridian", 9.0],
   PARAMETER["latitude_of_origin", 0.0],
   PARAMETER["scale_factor", 0.9996],
   PARAMETER["false_easting", 1500000.0],
   PARAMETER["false_northing", 0.0],
   UNIT["m", 1.0],
   AXIS["Easting", EAST],
   AXIS["Northing", NORTH],
   AUTHORITY["EPSG","3003"]]

e prj4text:

+proj=tmerc +lat_0=0 +lon_0=9 +k=0.999600 +x_0=1500000 +y_0=0
+ellps=intl +units=m +no_defs
+towgs84=104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68

A questo punto vi chiedo anche una cosa che magari qualcuno sa come
ricostruire in modo programmatico. GRASS nel PROJ_INFO esattamente cosa
contiene? Qualcuno ha idea di come si possa ricostruire senza l'ausilio
delle proj4?

Ciao,
Andrea

neteler · July 24, 2007, 12:38pm

On Tue, Jul 24, 2007 at 02:01:28PM +0200, Andrea Antonello wrote:

Ciao, ho un dubbio riguardo a come fare ad avere in postgis una gauss
boaga con l'informazione del towgs per le trasformazioni. Non sono
riuscito a trovare tale info in una delle proiezioni.
Quindi parto dalla proiezione identificata da epsg 3003.

Info:
GRASS contiene parametri TOWGS84 che mancano in QGIS, PostGIS etc.

Mi basta aggiungere a srtext e proj4text i towgs ad entrambi, vero?

PROJCS["Monte Mario / Italy zone 1 / Pensinsular optimized",
   GEOGCS["Monte Mario",
     DATUM["Monte Mario",
       SPHEROID["International 1924", 6378388.0, 297.0,
AUTHORITY["EPSG","7022"]],
      TOWGS84[104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68]
       AUTHORITY["EPSG","6265"]],
     PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
     UNIT["degree", 0.017453292519943295],
     AXIS["Geodetic longitude", EAST],
     AXIS["Geodetic latitude", NORTH],
     AUTHORITY["EPSG","4265"]],
   PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]],
   PARAMETER["central_meridian", 9.0],
   PARAMETER["latitude_of_origin", 0.0],
   PARAMETER["scale_factor", 0.9996],
   PARAMETER["false_easting", 1500000.0],
   PARAMETER["false_northing", 0.0],
   UNIT["m", 1.0],
   AXIS["Easting", EAST],
   AXIS["Northing", NORTH],
   AUTHORITY["EPSG","3003"]]

e prj4text:

+proj=tmerc +lat_0=0 +lon_0=9 +k=0.999600 +x_0=1500000 +y_0=0
+ellps=intl +units=m +no_defs
+towgs84=104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68

... un problema di QGIS, PostGIS etc...

A questo punto vi chiedo anche una cosa che magari qualcuno sa come
ricostruire in modo programmatico. GRASS nel PROJ_INFO esattamente cosa
contiene?

GRASS 6.3.cvs (pat):~ > g.proj -w
PROJCS["Transverse Mercator",
    GEOGCS["international",
        DATUM["Monte_Mario",
            SPHEROID["International_1924",6378388,297],
            TOWGS84[-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",9],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",1500000],
    PARAMETER["false_northing",0],
    UNIT["meter",1]]

GRASS 6.3.cvs (pat):~ > g.proj -p
-PROJ_INFO-------------------------------------------------
name : Transverse Mercator
datum : rome40
datumparams: towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68
proj : tmerc
ellps : international
a : 6378388.0000000000
es : 0.0067226700
f : 297.0000000000
lat_0 : 0.0000000000
lon_0 : 9.0
k_0 : 0.9996000000
x_0 : 1500000.0000000000
-PROJ_UNITS------------------------------------------------
unit : meter
units : meters
meters : 1.0

I dati stanno li':
grep rome grass63/lib/gis/datumtransform.table
rome40 "towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68" "Italy (Peninsular Part)" "Accuracy 3-4m"
rome40 "towgs84=-168.6,-34.0,38.6,-0.374,-0.679,-1.379,-9.48" "Italy (Sardinia)" "Accuracy 3-4m"
rome40 "towgs84=-50.2,-50.4,84.8,-0.690,-2.012,0.459,-28.08" "Italy (Sicily)" "Accuracy 3-4m"

Qualcuno ha idea di come si possa ricostruire senza l'ausilio
delle proj4?

Potresti usare GRASS con SWIG Utilizzando lib/proj/ che contiene
il codice aggiunto.

Magari potresti anche usare GDAL/OGR con SWIG.

ciao
Markus

PS: Ieri ho aggiornato GRASS 6.3-CVS su EPSG 6.13, seguendo GDAL.

------------------
ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler
ITC -> since 1 March 2007 Fondazione Bruno Kessler
------------------

moovida · July 24, 2007, 1:15pm

Ciao, ho un dubbio riguardo a come fare ad avere in postgis una gauss
boaga con l'informazione del towgs per le trasformazioni. Non sono
riuscito a trovare tale info in una delle proiezioni.
Quindi parto dalla proiezione identificata da epsg 3003.

Info:
GRASS contiene parametri TOWGS84 che mancano in QGIS, PostGIS etc.

Si', e' da li' che li ho presi ricordando una nostra chiaccherata di
qualche anno fa.

[...]

Qualcuno ha idea di come si possa ricostruire senza l'ausilio
delle proj4?

Potresti usare GRASS con SWIG Utilizzando lib/proj/ che contiene
il codice aggiunto.

Magari potresti anche usare GDAL/OGR con SWIG.

Il problema e' che chi usa metodi diversi per le proiezioni, fa
fatica ad interfacciarsi al workspace GRASS.
Udig si basa sul WKT con i parametri towgs per le riproiezioni.
Ora supporta le location di GRASS, ma non riesco a fargli capire la
proiezione senza fare uso delle proj. Come mai GRASS non usa qualcosa di
piu' accessibile? C'e' un motivo tecnico?
Al momento all'import di una location, viene fatta scegliere la
proiezione e scritta in un PROJ_INFO.WKT per uso JGrass-udig. Non e'
carino, ma non saprei proprio come fare altrimenti...

Qualche idea?
Andrea

ciao
Markus

PS: Ieri ho aggiornato GRASS 6.3-CVS su EPSG 6.13, seguendo GDAL.

------------------
ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler
ITC -> since 1 March 2007 Fondazione Bruno Kessler
------------------

_______________________________________________
Gfoss mailing list: 235 iscritti (23-07-2007)
Gfoss@faunalia.com
http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss

Leonardo_Lami · July 24, 2007, 2:26pm

Ciao a tutti,
mi inserisco nella conservazione con un problemino aggiuntivo:

stavo provando ad inserire un nuovo SRID nella tabella spatial_ref_sys
di un database di PostGis che contiene tutti gli SRID standard EPSG, il
nuovo SRID corrisponde al Gauss Boaga parte peninsulare sviluppato da
Markus che presenta la stringa:

+proj=tmerc +lat_0=0 +lon_0=9 +k=0.999600 +x_0=1500000 +y_0=0
+ellps=intl +units=m +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68

Ho provato poi a convertire uno strato di PostGIS in latlong wgs84 (SRID
4326) in questa nuova proiezione tramite il comando "transform",
visualizzando il tutto in QGIS ottengo un errore di circa 178 metri
verso sud.

Altre conversioni tra sistemi di proiezioni diverse danno invece dei
risultati validi. Il tutto è stato testato su QGIS.
Sembrerebbe che i fattori di correzzione inseriti non vengano letti in
maniera corretta da PostGis o comunque dal suo sistema di riproiezione,
mentre funzionano in modo corretto su QGIS e GRASS.

Qualcuno ha mai affrontato questo problema?

Ciao e grazie

Leonardo

Andrea Antonello ha scritto:

Ciao, ho un dubbio riguardo a come fare ad avere in postgis una gauss
boaga con l'informazione del towgs per le trasformazioni. Non sono
riuscito a trovare tale info in una delle proiezioni.
Quindi parto dalla proiezione identificata da epsg 3003.


Info:
GRASS contiene parametri TOWGS84 che mancano in QGIS, PostGIS etc.

Si', e' da li' che li ho presi ricordando una nostra chiaccherata di
qualche anno fa.

[...]


Qualcuno ha idea di come si possa ricostruire senza l'ausilio
delle proj4?


Potresti usare GRASS con SWIG Utilizzando lib/proj/ che contiene
il codice aggiunto.

Magari potresti anche usare GDAL/OGR con SWIG.

Il problema e' che chi usa metodi diversi per le proiezioni, fa
fatica ad interfacciarsi al workspace GRASS.
Udig si basa sul WKT con i parametri towgs per le riproiezioni.
Ora supporta le location di GRASS, ma non riesco a fargli capire la
proiezione senza fare uso delle proj. Come mai GRASS non usa qualcosa di
piu' accessibile? C'e' un motivo tecnico?
Al momento all'import di una location, viene fatta scegliere la
proiezione e scritta in un PROJ_INFO.WKT per uso JGrass-udig. Non e'
carino, ma non saprei proprio come fare altrimenti...

Qualche idea?
Andrea

ciao
Markus

PS: Ieri ho aggiornato GRASS 6.3-CVS su EPSG 6.13, seguendo GDAL.

------------------
ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler
ITC -> since 1 March 2007 Fondazione Bruno Kessler
------------------

_______________________________________________
Gfoss mailing list: 235 iscritti (23-07-2007)
Gfoss@faunalia.com
http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss

_______________________________________________
Gfoss mailing list: 235 iscritti (23-07-2007)
Gfoss@faunalia.com
http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss

neteler · July 24, 2007, 3:40pm

On Tue, Jul 24, 2007 at 03:15:31PM +0200, Andrea Antonello wrote:
...

Il problema e' che chi usa metodi diversi per le proiezioni, fa
fatica ad interfacciarsi al workspace GRASS.
Udig si basa sul WKT con i parametri towgs per le riproiezioni.
Ora supporta le location di GRASS, ma non riesco a fargli capire la
proiezione senza fare uso delle proj. Come mai GRASS non usa qualcosa di
piu' accessibile? C'e' un motivo tecnico?

Mhh, non ti capisco...

g.manual g.proj
-> WKT (anche per creare una nuova location)

Ciao
Markus

------------------
ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler
ITC -> since 1 March 2007 Fondazione Bruno Kessler
------------------

moovida · July 24, 2007, 3:48pm

Il mio problema e' che potrei decidere di installare jgrass-udig su
windows senza avere i vari grass, proj, etc. A questo punto se importo
una location di GRASS, non sono in grado di riconoscere esattamente la
proiezione dal file PROJ_INFO.

Dai dati che ci sono in PROJ_INFO non sono in grado di "costruire" un
WKT della proiezione. Giusto?

Ciao
Andrea

Markus Neteler probaly wrote:

On Tue, Jul 24, 2007 at 03:15:31PM +0200, Andrea Antonello wrote:
...

Il problema e' che chi usa metodi diversi per le proiezioni, fa
fatica ad interfacciarsi al workspace GRASS.
Udig si basa sul WKT con i parametri towgs per le riproiezioni.
Ora supporta le location di GRASS, ma non riesco a fargli capire la
proiezione senza fare uso delle proj. Come mai GRASS non usa qualcosa di
piu' accessibile? C'e' un motivo tecnico?

Mhh, non ti capisco...

g.manual g.proj
-> WKT (anche per creare una nuova location)

Ciao
Markus

------------------
ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler
ITC -> since 1 March 2007 Fondazione Bruno Kessler
------------------

neteler · July 25, 2007, 4:47pm

Andrea Antonello wrote on 07/24/2007 05:48 PM:

Il mio problema e' che potrei decidere di installare jgrass-udig su
windows senza avere i vari grass, proj, etc. A questo punto se importo
una location di GRASS, non sono in grado di riconoscere esattamente la
proiezione dal file PROJ_INFO.

Dai dati che ci sono in PROJ_INFO non sono in grado di "costruire" un
WKT della proiezione. Giusto?

Secondo me potresti essere in grado di estrare la funzionalita` da
g.proj che chiama OGR.
Perche leggia sia PROJ_INFO, sia WKT. Tutto GPL ovviamente...
PROJ/OGR funzionano anche su MS-Windows. Ma magari non ho capito bene il
problema... ti propongo un caffe` a Trento/Povo

ciao
markus

------------------
ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler
ITC -> since 1 March 2007 Fondazione Bruno Kessler
------------------

Leonardo_Lami · July 26, 2007, 4:13pm

Ciao a tutti,
ripropongo il problema per sapere se qualcuno (in particolare chi lavora
con PostGis abbia mai affrontato il problema dell'inseriemento di nuovi
SRID in Postgis e se questi funzionino nella riproiezione.

Ho provato ad inserire un nuovo SRID nella tabella spatial_ref_sys
di un database di PostGis che contiene tutti gli SRID standard EPSG, il
nuovo SRID corrisponde al Gauss Boaga parte peninsulare sviluppato da
Markus che presenta la stringa:

+proj=tmerc +lat_0=0 +lon_0=9 +k=0.999600 +x_0=1500000 +y_0=0
+ellps=intl +units=m +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68

Ho provato poi a convertire uno strato di PostGIS in latlong wgs84 (SRID
4326) in questa nuova proiezione tramite il comando "transform",
visualizzando il tutto in QGIS ottengo un errore di circa 178 metri
verso sud.

Altre conversioni tra sistemi di proiezioni diverse danno invece dei
risultati validi. Il tutto è stato testato su QGIS.
Sembrerebbe che i fattori di correzzione inseriti non vengano letti in
maniera corretta da PostGis o comunque dal suo sistema di riproiezione,
mentre funzionano in modo corretto su QGIS e GRASS.

Sapete se la cosa può funzioanre? o per caso avete affrontato un problema
simile??

Ciao e grazie

Leonardo

Andrea Antonello ha scritto:

Ciao, ho un dubbio riguardo a come fare ad avere in postgis una gauss
boaga con l'informazione del towgs per le trasformazioni. Non sono
riuscito a trovare tale info in una delle proiezioni.
Quindi parto dalla proiezione identificata da epsg 3003.

Info:
GRASS contiene parametri TOWGS84 che mancano in QGIS, PostGIS etc.

Si', e' da li' che li ho presi ricordando una nostra chiaccherata di
qualche anno fa.

[...]

Qualcuno ha idea di come si possa ricostruire senza l'ausilio
delle proj4?

Potresti usare GRASS con SWIG Utilizzando lib/proj/ che contiene
il codice aggiunto.

Magari potresti anche usare GDAL/OGR con SWIG.

Il problema e' che chi usa metodi diversi per le proiezioni, fa
fatica ad interfacciarsi al workspace GRASS.
Udig si basa sul WKT con i parametri towgs per le riproiezioni.
Ora supporta le location di GRASS, ma non riesco a fargli capire la
proiezione senza fare uso delle proj. Come mai GRASS non usa qualcosa di
piu' accessibile? C'e' un motivo tecnico?
Al momento all'import di una location, viene fatta scegliere la
proiezione e scritta in un PROJ_INFO.WKT per uso JGrass-udig. Non e'
carino, ma non saprei proprio come fare altrimenti...

Qualche idea?
Andrea

ciao
Markus

PS: Ieri ho aggiornato GRASS 6.3-CVS su EPSG 6.13, seguendo GDAL.

------------------
ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler
ITC -> since 1 March 2007 Fondazione Bruno Kessler
------------------

_______________________________________________
Gfoss mailing list: 235 iscritti (23-07-2007)
Gfoss@faunalia.com
http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss

_______________________________________________
Gfoss mailing list: 235 iscritti (23-07-2007)
Gfoss@faunalia.com
http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss