I file Shapefile e OpenStreetMap non si sovrappongono in QGIS?


10

Quindi sono abbastanza nuovo per GIS e sto lavorando a un progetto di mappatura usando dati per Vancouver, in Canada. Voglio combinare i dati OpenStreetMap con i dati forniti dalla città di Vancouver. Ma i dati della città non si allineano con i dati OSM quando li mappa in QGIS.

Seguire il consiglio di una risposta alle mappe proiettate dovrebbe ma non si sovrappongono in ArcMap? , Ho verificato le due origini dati rispetto a un terzo (una mappa del mondo). I dati OSM sono giusti - la Columbia Britannica è nel posto giusto - ma i dati della città mettono Vancouver nel sud della Francia.

Le mie fonti di dati sono:

  1. Piste ciclabili e biblioteche della città di Vancouver . (Entrambi sono file SHP.) Questi sono i dati che non si allineano con gli altri.
  2. Una mappa della Columbia Britannica da Open Street Map. Ho inserito i dati in un database usando osm2pgsql.
  3. Una mappa del mondo che fa parte del progetto Open Street Style.

proiezioni:

  1. I dati della città sono proiettati in UTM zone10, NAD83 (CSRS). L'ho selezionato nella schermata Proprietà layer di QGIS per quei layer. Stringa di proiezione:+proj=utm +zone=5 +ellps=GRS80 +datum=NAD83 +units=m +no_defs
  2. Le mappe OSM sono in proiezione Mercator. QGIS ha molte scelte di Mercator e non ero sicuro di quale usare, quindi ho indovinato "Popular Visualization CRS / Mercator".+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
  3. Per la mappa del mondo, Quantum GIS ha caricato la proiezione dal file .prj. Quello è+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +over +no_defs

Le mappe proiettate dovrebbero ma non si sovrappongono in ArcMap? parla di problemi di zona. Sembra quello che potrei sperimentare. Quindi ho provato a cambiare le librerie dalla zona UTM 10N ad altre zone (15N e 60N) ma non è cambiato nulla.

Qualcuno ha suggerimenti su cosa provare dopo?


2
Se i tuoi dati di Vancouver si trovano nella zona UTM 10N, si suppone di leggere la stringa del proj +proj=utm +zone=10e non +zone=5perché penso che la zona 5 sia molto vicina a qualche parte del Regno Unito (inizia con una G ...) ed è per questo che arrivano i tuoi dati in Francia
SaultDon

Puoi andare su lrdw.ca e scaricare il file di forma delle zone UTM con una proiezione BC Albers. Accedi come ospite e vai alla pagina "I migliori prodotti"; è vicino alla cima. Inserisci quel file nel tuo progetto QGIS e vedi se i tuoi due set di dati rientrano nella zona 10N - usando BC Albers come trasformazione CRS "al volo" (devi prima abilitare Maiusc + Ctrl + P, quindi metti un segno di spunta nella casella richiesta) .
SaultDon

Risposte:


11

Quello che dobbiamo davvero scoprire qui sono le proiezioni dei tuoi dati di input.

  1. Dati della città

    Lo Shapefile viene fornito con un file .prj che è perfettamente interpretato da QGIS. Non modificare le impostazioni del livello!

  2. Dati OSM

    Se si utilizzano i dati vettoriali OSM, si trovano nelle coordinate geografiche WGS84 (EPSG: 4326). Solo le tessere raster sono in Web Mercator (EPSG: 900913).

  3. Apri dati Street Style

    Se QGIS ha riconosciuto il contenuto .prj, lascialo così.

Infine, ho scelto un sistema di coordinate del progetto (suggerirei quello UTM) e attivare la riproiezione al volo come già spiegato da @Mike_Toews.


6

Devi riproiettare i diversi livelli nella stessa proiezione "al volo". Vai a:

Impostazioni> Proprietà progetto ...> Sistema di riferimento coordinate (CRS)

e verificare su Enable 'al volo' CRS trasformazione


1
Un piccolo chiarimento: l'uso di "al volo" non riproietterà nel vero senso poiché i file manterranno la loro proiezione di origine, ma saranno invece "disegnati" dinamicamente con una diversa proiezione specificata. Questa risposta ti aiuterà a iniziare a risolvere i problemi con QGIS.
SaultDon

1
ignoro la maggior parte del mio commento sopra perché le sue parole sono in realtà corrette, è il mio cervello che non era = P
SaultDon

1

Ho appena avuto questo problema aggiungendo due raster, uno in UTM e uno in Geographic. Ho seguito le indicazioni abilitandole al volo, ma non è successo nulla. Ho rimosso i livelli e li ho nuovamente aggiunti e il problema è stato risolto


Grazie John Ho avuto problemi a sovrapporre due livelli vettoriali e ho provato a rimuovere e aggiungere nuovamente i livelli (poiché nient'altro ha funzionato). Questo non ha funzionato da solo, MA quando ho abbreviato il nome dello shapefile che stavo salvando, andava tutto bene.

Ti preghiamo di non aggiungere "grazie" come risposta. Una volta che avrai una reputazione sufficiente , potrai votare le domande e le risposte che hai trovato utili. - Dalla recensione
Midavalo
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.