[Gfoss] projection WMS client in mapserver

[Emanuele Masiero]

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 a 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.