[Gfoss] proj parametri da UTM a GaussBoaga - in mapserver

Local Chapters GFOSS.it
1

ciao a tutti,
sto facendo delle prove con mapserver
sto visualizzando un layer WMS cliente del PCN minambiente
(http://wms.pcn.minambiente.it/cgi-bin/mapserv.exe?map=/ms_ogc/service/ortofoto_colore_06_f32.map)

però c'è un errore di posizione di circa 30-70 metri!!
il layer viene pubblicato sia in WGS84 e sia in UTM/WGS84 (fuso 32-33)

ho letto la risposta di antonio falciano
http://www.nabble.com/ancora%3A-QGIS-UTM-ED50-Emilia-Romagna-td19589110.html#a19623878,
che mi ha fatto capire molte cose
però mi rimangono ancora molti dubbi:
1) con mapserver e un layer WMS client, è fattibile forzare la proiezione
sostituendo il codice EPSG con i parametri ?
2) perchè se sostituisco il codice 32632 con la retiva stringa presa dal
proj "+proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs",
mapserver va in errore?
3) entrando nel merito di come è composto il proj, esiste una guida che
illustra come vanno combinati i parametri: elissoide, datum, primo
meridiano, ecc,
4) da wgs84 (coordinate geografiche) come faccio a correggerlo per
proiettarlo in gauss boaga?
5) da UTM (coordinate piane)?

ho fatto qualche esame di topografia e cartografia... ma le domande sono
rimaste senza risposta... :frowning:

mi rendo conto che questo tipo di domande riguarda la cartografia e i
sistemi di riferimento
però vi chiedo se avete qualche riferimento dettagliato da consultare

qui sotto il layer di esempio pe ril collegamento WMS al PCN del MinAmbiente

grazie mille a chi mi darà qualche spunto per chiarirmi definitivamente le
idee!!
emanuele masiero

LAYER
  NAME "wms_pcn_ortofoto_06"
  TYPE RASTER
  STATUS OFF
  CONNECTION
"http://wms.pcn.minambiente.it/cgi-bin/mapserv.exe?map=/ms_ogc/service/ortofoto_colore_06_f32.map"
  CONNECTIONTYPE WMS
  METADATA
    "DESCRIPTION" "ortofoto 2006 PNC MinAmbiente"
    "wms_srs" "EPSG:32632"
    "wms_name" "ortofoto_colore_06_f32"
    "wms_server_version" "1.1.1"
    "wms_format" "image/jpeg"
  END
  
  PROJECTION
    "init=epsg:32632"
    #"+proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"
#<32632>
  END
  
END
--
View this message in context: http://www.nabble.com/proj-parametri-da-UTM-a-GaussBoaga---in-mapserver-tp22331258p22331258.html
Sent from the Gfoss mailing list archive at Nabble.com.

2

Emanuele Masiero ha scritto:

ciao a tutti,
sto facendo delle prove con mapserver
sto visualizzando un layer WMS cliente del PCN minambiente
(http://wms.pcn.minambiente.it/cgi-bin/mapserv.exe?map=/ms_ogc/service/ortofoto_colore_06_f32.map)

però c'è un errore di posizione di circa 30-70 metri!!

rispetto a cosa?

il layer viene pubblicato sia in WGS84 e sia in UTM/WGS84 (fuso 32-33)

ho letto la risposta di antonio falciano
http://www.nabble.com/ancora%3A-QGIS-UTM-ED50-Emilia-Romagna-td19589110.html#a19623878,
che mi ha fatto capire molte cose
però mi rimangono ancora molti dubbi:
1) con mapserver e un layer WMS client, è fattibile forzare la proiezione
sostituendo il codice EPSG con i parametri ?

certo, il codice EPSG o la stringa proj che ci sta dietro sono equivalenti

2) perchè se sostituisco il codice 32632 con la retiva stringa presa dal
proj "+proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs",
mapserver va in errore?

forse si tratta di un banale problema di formattazione:
http://mapserver.org/mapfile/projection.html

3) entrando nel merito di come è composto il proj, esiste una guida che
illustra come vanno combinati i parametri: elissoide, datum, primo
meridiano, ecc.

No, esiste la documentazione di proj.4:
http://trac.osgeo.org/proj/

4) da wgs84 (coordinate geografiche) come faccio a correggerlo per
proiettarlo in gauss boaga?

se cerchi in questa lista troverai decine di risposte a questa domanda

5) da UTM (coordinate piane)?

come sopra

ciao
Antonio

3

grazie antonio!!
sono riuscito a risolvere il problema, ma mi restano ancora dei dubbi e ,
approffittando ancora di te antonio, vorrei riformulare il quesito anche per
chiarezza di altri:

ho realizzato un web-gis in pmapper/mapserver con dati shape in gauss boaga
ovest (3003)
ho aggiunto il layer WMS del minambiente che è in UTM32/WGS84 (32632)

nel mapfile ho indicato la PROJECTION :
'init=epsg:3003' per tutti i LAYER e per la MAP
'init=epsg:32632' per il layer WMS

con questa configurazione ho uno scostamento geometrico di circa x=-20m
y=x80

nella PROJECTION della MAP ho sostituito 'init=epsg:3003' con i parametri
corrispondenti aggiungendo in coda il +towgs84=.... questa è la stringa:
'+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +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'

in questo modo si è eliminato lo scostamento tra dati locali e WMS!!!

ora vorrei fare alcune osservazioni/domande per approfondire la questione:
1
inizialmente la mia idea era quella di forzare la projection del layer WMS
definendo dei parametri personalizzati, perchè dal mio punto di vista era
lui (il layer WMS) che non si proiettava correttamente sul gauss boaga ovest
(3003)
questa cosa però non mi è mai riuscita perchè mapserver mi ha sempre dato
errore (Layer must have an EPSG or AUTO projection code (in its PROJECTION
object or wms_srs metadata)
ho provato a mettere un codice epsg diverso per vedere il comportamento, in
questo caso mapserver visualizzava la mappa però al posto di visualizzare il
layer WMS compariva un messaggio (WMS server error. invalid SRS given : SRS
must be valid for all requested layers.)

la mia conclusione quindi è: non è possibile modificare il codice epsg lato
client di un layer WMS (qualcuno confema?)

2
a questo punto ho fatto un altro tipo di prova: ho aggiunto il parametro
+towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 nella projection di un
layer 3003 (gauss-boaga ovest)
mi aspettavo che si verificasse uno scostamento tra questo layer e qualunque
altro layer 3003 (non modificato) e che ci fosse corrispondenza con il layer
WMS di prima (min. Ambiente)
e invece non è successo nulla: le coordinate sono rimaste sempre le stesse!
pensavo che cambianto il datum di un layer (anche se le coordinate piane
rimanevano invariate) ci fosse qualche "movimento" rispetto ad un altro
layer.
c'è qualcuno che sa darmi una spiegazione "facile"?? o che ha qualche spunto
di riflessione?
esistono corsi di 1-2 giornate specifici sui sistemi di riferimento
geografici?

ringrazio per la pazienza (soprattutto antonio f.)
emanuele masiero

Antonio Falciano wrote:

Emanuele Masiero ha scritto:

ciao a tutti,
sto facendo delle prove con mapserver
sto visualizzando un layer WMS cliente del PCN minambiente
(http://wms.pcn.minambiente.it/cgi-bin/mapserv.exe?map=/ms_ogc/service/ortofoto_colore_06_f32.map)

però c'è un errore di posizione di circa 30-70 metri!!

rispetto a cosa?

il layer viene pubblicato sia in WGS84 e sia in UTM/WGS84 (fuso 32-33)

ho letto la risposta di antonio falciano
http://www.nabble.com/ancora%3A-QGIS-UTM-ED50-Emilia-Romagna-td19589110.html#a19623878,
che mi ha fatto capire molte cose
però mi rimangono ancora molti dubbi:
1) con mapserver e un layer WMS client, è fattibile forzare la proiezione
sostituendo il codice EPSG con i parametri ?

certo, il codice EPSG o la stringa proj che ci sta dietro sono equivalenti

2) perchè se sostituisco il codice 32632 con la retiva stringa presa dal
proj "+proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs",
mapserver va in errore?

forse si tratta di un banale problema di formattazione:
http://mapserver.org/mapfile/projection.html

3) entrando nel merito di come è composto il proj, esiste una guida che
illustra come vanno combinati i parametri: elissoide, datum, primo
meridiano, ecc.

No, esiste la documentazione di proj.4:
http://trac.osgeo.org/proj/

4) da wgs84 (coordinate geografiche) come faccio a correggerlo per
proiettarlo in gauss boaga?

se cerchi in questa lista troverai decine di risposte a questa domanda

5) da UTM (coordinate piane)?

come sopra

ciao
Antonio
_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
Gfoss@faunalia.com
http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
I messaggi di questa lista non rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.

--
View this message in context: http://www.nabble.com/proj-parametri-da-UTM-a-GaussBoaga---in-mapserver-tp22331258p22387033.html
Sent from the Gfoss mailing list archive at Nabble.com.