Ciao sono l'unico ad avere problemi con la proiezione Gauss Boaga?
Ho un dataset (si tratta di un file xls) con una serie di punti per i quali sono riportate le coordinate. Ad occhio credo che siano in Gauss Boaga Fuso Est (sono in Italia Centrale ed uno dei miei dati ha coordinate 4740352 2275360). Ho inoltre una mappa *tif per la quale ho anche i file *rrd e *tfw. Tale mappa dovrebbe essere nel medesimo sistema di coordinate, in quanto i punti si dispongono coerentemente nella mappa. Il mio obiettivo è creare degli *shp dal dataset in *xls, cosa fra l'altro già fatta altre volte lavorando nel sistema WGS84. Vado su qgis, proprietà progetto, coordinate di riferimento spaziale Transverse Mercator - Monte Mario/Italy zone 2 (EPSG 3004). Poi, dopo aver trasformato il mio *xls in *csv, con il plugin "aggiungi layer testo delimitato", carico il mio dataset. Poi vado su layer -> salva con nome, formato ESRI Shapefile, seleziono il nome, CRS "Transverse Mercator - Monte Mario/Italy zone 2 (EPSG 3004)" e mi compare la seguente schermata: "Esportazione nel file vettoriale fallita. Errore: Failed to transform a point while drawing a feature of type ''. Writing stopped. (Exception: trasformazione diretta di (39712.5, 82734.8) fallito con errore: latitude or longitude exceeded limits)" Creandomi un shp vuoto ed un prj con la seguente scritta: "PROJCS ["Transverse_Mercator",GEOGCS["GCS_International 1909 (Hayford)",DATUM ["D_unknown",SPHEROID["intl",6378388,297]],PRIMEM["Greenwich",0],UNIT ["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",15], PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",2520000], PARAMETER["false_northing",0],UNIT["Meter",1]]" il che mi sembra coerente a quello che mi aspetterei (vedi link) http://spatialreference. org/ref/sr-org/4686/ se invece setto "CRS originale", mi crea un *shp "coincidente" con il mio *csv (se carico i due file su qgis, si sovrappongono perfettamente), ma il *prj creato non è coerente "GEOGCS ["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984", 6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree", 0.017453292519943295]]" Non contento, vado su Grass e creo una nuova location usando il codice EPSG 3004, settando "Italia Peninsulare" quando mi chiede di specificare (??) carico il *tiff (con "r.in.gdal") e sulla shell mi compare la seguente schermata: ERRORE: Projection of dataset does not appear to match current location. Location PROJ_INFO is: name: Transverse Mercator proj: tmerc datum: rome40 ellps: international lat_0: 0 lon_0: 15 k: 0.9996 x_0: 2520000 y_0: 0 no_defs: defined towgs84: -104.1,-49.1,-9.9,0.971, -2.917,0.714,-11.68 Import dataset PROJ_INFO is: cellhd.proj = 0 (unreferenced/unknown) You can use the -o flag to r.in.gdal to override this check and use the location definition for the dataset. Consider generating a new location from the input dataset using the 'location' parameter. Settando quindi "-o", che vorrebbe dire: sovrascrivi proiezione (usa proiezione della location)" sembrerebbe che la mappa è stata caricata "correttamente". Discorso analogo con il file *shp, quello che nel *prj sembra che sia in WGS84, tale file è stato caricato correttamente usando "-o". Per concludere se esporto il file *shp da GRASS, mi si crea un *shp con le coordinate ed un *prj giusti... In sostanza, con questo escamotage riesco a risolvere il problema, ma ho paura che "barando" non abbia considerato qualcosa ed alla fine ho paura di aver inevitabilmente commesso qualche errore... scusate la mail chilometrica, ma magari qualcuno si è trovato nella mia medesima situazione.... grazie mille!! marco Supera i limiti: raddoppia la velocità da 10 a 20 Mega! Risparmia con Tutto Incluso: telefono + adsl 20 mega a soli 29,95 € al mese per due anni! SCONTO DI 240 EURO! http://abbonati.tiscali.it/telefono-adsl/prodotti/tc/tuttoincluso/?WT.mc_id=01fw _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione Gfoss@lists.gfoss.it http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss Questa e' una lista di discussione pubblica aperta a tutti. Non inviate messaggi commerciali. I messaggi di questa lista non rispecchiano necessariamente le posizioni dell'Associazione GFOSS.it. 485 iscritti al 20.11.2010