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