A beneficio di altri:<br><br>Vista la situazione abbastanza caotica prodotta dalle nuove proj4.8.0.<br>Ho ritenuto di fare delle prove per vedere come si puo' compensare la cosa, facendo uso dei grigliati IGM (in formato NTv2)<br>
<br><a href="http://87.30.244.175/nota_tecnica_NTv2.pdf">http://87.30.244.175/nota_tecnica_NTv2.pdf</a><br><br>Il banco di prova per i tests ho usato un mapserver sotto Linux:<br><br>Il test si è sviluppato cosi':<br>
dati di partenza in epsg:25832 (realizzazione etrs89-etrf2000) alias: UTM-ETRF2000<br><br>dati attesi in output in epsg:3003 (alias GaussBoaga) e/o in epsg:4326 (alias WGS84)<br><br>Ipotesi di partenza:<br><br>nella installazione di proj si e' fatto uso delle cartelle di default e quindi la definizioni epsg sono contenute nella cartella:<br>
/usr/local/share/proj<br><br>Esecuzione:<br><br>Step 1)<br><br>supponendo che il file con i grigliati IGM si chiami "44101110_44301150_R40_S89.gsb", tale file va copiato nella cartella di default <br>delle definizioni epsg (/usr/local/share/proj).<br>
Attenzione che potrebbero essere piu' di uno ogni file racchiude i grigliati per passare da un sistema di riferimento a un altro.<br>In questo caso si parla di un file che contiene i grigliati per passare da Roma40 a ETRS89.<br>
<br>Ovviamente se ci sono piu' file andranno copiati tutti quanti e usato di volta in volta quello opportuno.<br><br>Step 2)<br><br>sempre nella cartella /usr/local/share/proj<br>si apre il file "epsg".<br><br>
In esso si <br>modifica la definizione di epsg:3003 cambiandola da <br><br><3003> +proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +units=m +no_defs  <><br>
<br>a<br><br><3003> +proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl +nadgrids=44101110_44301150_R40_S89.gsb +units=m +no_defs  <><br><br>Si lascia invece invariata la definizione di epsg:25832 cosi' come è proposta nel file epsg:<br>
<br><25832> +proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs  <><br><br>Si lascia altresi' invariata la definzione di epsg:4326<br>come proposta nel file:<br><br><4326> +proj=longlat +datum=WGS84 +no_defs  <><br>
<br>Step 3)<br><br>imposto un progettino MapServer per vedere dei risultati di riproiezione al volo:<br><br>Nel mapfile imposto il seguente parametro per i dati in input:<br><br>  PROJECTION<br>    "+init=epsg:25832 +towgs84=0,0,0,0,0,0,0"<br>
  END<br><br>inizio le prove:<br><br>Step 4)<br><br>Provando con Mapserver a farsi inviare una mappa riproiettata al volo da Mapserver su epsg:3003 usando dati nativi in epsg:25832 il risultato collima molto bene <br>con quello atteso per dei dati nativi in gaussBoaga (a essere pessimista ci potrà essere qualche centimetro di differenza, ma probabilmente è inferiore)<br>
<br>Step5)<br><br>Provando a farsi inviare da mapserver una mappa su epsg:25832 (il sistema di riferimento nativo dei dati in input) il dato arriva naturlamente corretto.<br><br>Step 6)<br><br>Provando a farsi inviare da MapServer una mappa su epsg:4326 (sistema geografico), il dato arriva riproietato correttamente con una buona collimazione nei confronti del<br>
medesimo dato nativo in lat-long calcolato sulla realizzazione etrs89-etrf2000.<br><br>Step 7)<br>giusto per completezza ho provato a inserire i grigliati (tramite nadgrid sulla definizione di epsg:4326,<br>ma il risultato viene molto errato.<br>
<br>[Considerazioni finali]<br><br>Nella ipotesi che i dati in ingresso siano epsressi in epsg:25832 (utm-etrf2000),<br>e ipotizzando che i sistemi ammessi in output siano epsg:4326, epsg:3003 e epsg:25832,<br>i grigliati vanno inseriti solamente nella definizione di epsg:3003.<br>
<br>Mentre in quella epsg:4326 non servono (idem a maggior ragione in quella epsg:25832).<br><br>Una ultima considerazione:<br><br>Ovviamente se si è in situazione opposta:<br><br>dati nativi in gaussBoaga, e in output si vuole supportare epsg:3003, epsg:25832 e epsg:4326, i grigliati andranno inseriti nella sola definizione di epsg:25832<br>
<br>Andrea.<br><br clear="all"><br>-- <br>-----------------<br>Andrea Peri<br>. . . . . . . . . <br>qwerty àèìòù<br>-----------------<br><br>