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 > [email protected] > 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. _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione [email protected] 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.
