R: Gauss Boaga, problema di proiezione (RISOLTO)

la pausa pranzo mi ha portato consiglio

in qgis era rimasta settata la
proiezione WGS84 (quella nella quale normalmente lavoro)
e la procedura
"qgis, proprietà progetto, coordinate di riferimento spaziale" non
c'entra nulla con il settare il sistema di coordinate. Abituato a
lavorare sempre sullo stesso sistema non mi ero mai accorto che tale
procedura (che davo per assodata) era totalmente inutile (vedi capitolo
8 del manuale).

chiedo scusa
:-p

----Messaggio originale----
Da: a.
[hidden email]
Data: 02/12/2010 13.39
A: <[hidden email]>

Cc: <[hidden email]>
Ogg: R: [Gfoss] Gauss Boaga, problema di
proiezione

Hai impostato la proiezione del progetto ad EPSG 3004?


>----Messaggio originale----
>Da: [hidden email]
>Data:
02/12/2010 11.50
>A: "Mailing List Gfoss (ita)"<[hidden email]>

>Ogg: [Gfoss] Gauss Boaga, problema di proiezione
>
>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

>[hidden email]
>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







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
[hidden email]
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