<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Il tuo quadro di unione si trova in coordinate GEOGRAFICGE WGS84, non
sono PROIETTATE EPSG: 4326<br>
<br>
I dati istat sono in ED50 UTM32N, EPSG 23032<br>
<br>
Ora ne approfitto per fare un appunto perchè ho notato che sono in
molti a fanno fatica a distinguere i due tipi di coordinate:
geografiche e proiettate<br>
<br>
Le coordinate geografiche sono solitamente espresse in gradi LAT e LONG
costruite su di un elissoide di riferimento (per esempio WGS84), e  in
quanto tali non hanno una zona!!<br>
<br>
Ora questo elissoide non è altro che la forma matematica che va ad
approssimare la "forma" del pianeta. E' ovvio che a seconda delle zone
l'errore che si va a creare può essere maggiore o minore.<br>
<br>
Dato che un elissoide non può essere rappresentato sul piano, bisogna
proiettarlo su di un cilindro o un cono.<br>
<br>
La proiezione dunque permette di trasferire le coordinate da lat-long a
piane metriche. Per esempio UTM32N  WGS84 rappresenta coordinate
appartenenti all'elissoide di riferimento WGS84 riproiettate su di un
cilindro che utilizza il fuso 32 nord come centro, o "punto di
contatto" tra l'elissoide e il cilindro. Questo punto (che in questo
caso è una linea) ha distorsione 0 e mano mano che ci si allontana
questa distorsione aumenta.<br>
Ecco perchè esistono molti tipi di proiezioni  e diverse zone (32N,
33N, 32S, 33S, ecc)<br>
<br>
Aggiungerei una nota:<br>
<br>
Proiettare e riproiettare dati utilizzando lo stesso elissoide di
riferimento (per es.: WGS84  -> UTM32N WGS84 -> WGS84) non crea
errori,  e qualsiasi desktop gis è in grado di farlo. <br>
Nel caso in cui la riproiezione venga effettuata tra diversi elissoidi
(per es. WGS84 UTM32N -> ED50 UTM32N) questa può generare errori
anche di diversi metri. Per applicare una trasformazione corretta
bisognerebbe procurarsi l'elenco dei famosi 7 parametri, che vanno a
sostituire il processo di riproiezione. I sette parametri sono in
ordine: <br>
<ul>
  <li>X (spostamento di coordinate lungo l'asse X)</li>
  <li>Y (spostamento di coordinate lungo l'asse Y)</li>
  <li>Z (spostamento di coordinate lungo l'asse Z)</li>
  <li>dX (rotazione intorno all'asse X)</li>
  <li>dY (rotazione intorno all'asse Y)</li>
  <li>dZ (rotazione intorno all'asse Z)</li>
  <li>K fattore di scala</li>
</ul>
<br>
<div class="moz-signature">
<div class="moz-signature">
<table style="display: block;" align="left" border="0" cellpadding="0"
 cellspacing="0" width="100%">
  <tbody>
    <tr>
      <td colspan="2"> <font
 face="DejaVu Sans,Verdana,Arial,Helvetica,sans-serif" size="2"> <strong>Dr.
Pietro d'Orio</strong><br>
      <strong><font color="#ffa64c">R3 </font><font color="#9c9c9c">GIS</font></strong>
Srl - GmbH<br>
Via Johann Kravogl-Str. 2<br>
I-39012 Meran/Merano (BZ)<br>
web: <a href="http://www.r3-gis.com/" target="_blank">www.r3-gis.com</a><br>
Email: <a href="mailto:pietro.dorio@r3-gis.com" target="_blank">pietro.dorio@r3-gis.com</a><br>
Tel. : +39 0473 494949<br>
Fax : +39 0473 069902<br>
      </font></td>
    </tr>
    <tr>
      <td colspan="2"> <font color="#888888"
 face="DejaVu Sans,Verdana,Arial,Helvetica,sans-serif" size="1"><br>
      <strong>ATTENZIONE!</strong> Le informazioni contenute nella
presente e-mail e nei documenti eventualmente allegati sono
confidenziali. La loro diffusione, distribuzione e/o riproduzione da
parte di terzi, senza autorizzazione del mittente è vietata e può
violare il D. Lgs. 196/2003. In caso di ricezione per errore, Vogliate
immediatamente informare il mittente del messaggio e distruggere la
e-mail.<br>
      <br>
      <strong>ACHTUNG!</strong> Die in dieser Nachricht oder in den
beigelegten Dokumenten beinhalteten Informationen sind streng
vertraulich. Ihre Verbreitung und/oder ihre Wiedergabe durch Dritte ist
ohne Erlaubnis des Absenders verboten und verstößt gegen das
Legislativdekret 196/2003. Sollten Sie diese Mitteilung
irrtümlicherweise erhalten haben, bitten wir Sie uns umgehend zu
informieren und anschließend die Mitteilung zu vernichten.<br>
      <br>
      <strong>WARNING!</strong> This e-mail may contain confidential
and/or privileged information. If you are not the intended recipient
(or have received this e-mail in error) please notify the sender
immediately and destroy this e-mail. Any unauthorised copying,
disclousure or distribution of the material in this e-mail is strictly
forbidden and could be against the law (D. Lgs. 196/2003)<br>
      <br>
      <br>
      <br>
      <br>
      <br>
      </font></td>
    </tr>
  </tbody>
</table>
</div>
</div>
<br>
<br>
e-paul ha scritto:
<blockquote cite="mid:4D0F06E6.5030809@tiscali.it" type="cite">
  <pre wrap="">A questo punto... non ci sto più capendo nulla!
Scusate ma, come ho premesso, sono autodidatta e alle prime armi.

Dunque, ricapitolando;
- come faccio a sapere quale è il CRS associato al mio quadro di unione 
della CTR? Il file qu_10u33.prj contiene quanto segue:
GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]
quindi siamo in WGS 1984, ma in  che zona?
Il fatto che nel nome del file ci sia un "33" mi fa pensare che siamo 
nella zona 33, ma questa informazione è contenuta nel file?

- Il file ISTAT dei Comuni effettivamente fa riferimento alla zona 32N. 
Il file "com2010_s.prj" contiene quanto segue:
PROJCS["ED_1950_UTM_Zone_32N",GEOGCS["GCS_European_1950",DATUM["D_European_1950",SPHEROID["International_1924",6378388.0,297.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",9.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]

- riproiezione al volo vuol dire che i layer vengono correttamente 
rappresentati indipendentemente dal CRS adottato per la vista che sto 
utilizzando? Devo essere io comunque a specificare la trasformazione 
necessaria o il software individua i CRS di ciascun layer autonomamente?

Il 17/12/2010 10.22, Antonio Falciano ha scritto:
  </pre>
  <blockquote type="cite">
    <pre wrap="">Il 17/12/2010 9.20, Pietro d'Orio ha scritto:
    </pre>
    <blockquote type="cite">
      <pre wrap="">Io opererei in questo modo:

1. Apri una vista con codice EPSG 23030 e ci carichi i dati originali Istat
2. Apri il processo di riproiezioni a partire da 23030 al codice EPSG
della tua CTR
3. Chiudi la vista, ne fai una nuova con il codice EPSG della tua CTR e
carichi i nuovi confini riproiettati.
      </pre>
    </blockquote>
  </blockquote>
  <pre wrap=""><!---->
Alla fine ho aeguito un po' una via di mezzo tra le varie proposte che 
ho letto qui nella lista:
- ho creato una nuova vista con CRS 32633
- ho caricato il quadro di unione CTR (non dovrebbero essere necessarie 
riproiezioni, se ho capito bene, quindi non ne applico alcuna)
- ho caricato lo shapefile dei comuni ISTAT riproiettando al volo (cioè 
quando carico il file con trasformazione da WGS_84 a EPSG 32633

I layer sono ancora distanti 500 km :-(

Se invece riproietto (tirando a indovinare) il layer dei comuni ISTAT da 
WGS84 a EPSG 32634 si sovrappongono correttamente.

Qualcuno di voi potrebbe essere disposto a fare la stessa operazione di 
sovrapposizione per verificare i risultati? in caso posso inviare i 
files necessari in email privata.

Grazie ancora a tutti e scusate per il protrarsi di questa discussione!
Paolo Caporro
_______________________________________________
Gvsig_italian mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Gvsig_italian@lists.gvsig.org">Gvsig_italian@lists.gvsig.org</a>
<a class="moz-txt-link-freetext" href="https://gvsig.org/lists/mailman/listinfo/gvsig_italian">https://gvsig.org/lists/mailman/listinfo/gvsig_italian</a>
  </pre>
</blockquote>
</body>
</html>