Os dados do Shapefiles e do OpenStreetMap não se sobrepõem no QGIS?


10

Então, eu sou muito novo no GIS e estou trabalhando em um projeto de mapeamento usando dados para Vancouver, Canadá. Quero combinar os dados do OpenStreetMap com os dados fornecidos pela cidade de Vancouver. Mas os dados da cidade não se alinham com os dados do OSM quando os mapeio no QGIS.

Seguindo o conselho de uma resposta aos mapas projetados, o ArcMap não deve se sobrepor? , Verifiquei as duas fontes de dados contra uma terceira (um mapa do mundo). Os dados do OSM estão corretos - a Colúmbia Britânica está no lugar certo - mas os dados da City colocam Vancouver no sul da França.

Minhas fontes de dados são:

  1. Ciclovias e bibliotecas da cidade de Vancouver . (Ambos são arquivos SHP.) Estes são os dados que não se alinham com os outros.
  2. Um mapa da Colúmbia Britânica no mapa da rua aberta. Coloquei os dados em um banco de dados usando osm2pgsql.
  3. Um mapa do mundo que faz parte do projeto Open Street Style.

Projeções:

  1. Os dados da cidade são projetados na zona UTM10, NAD83 (CSRS). Selecionei isso na tela Propriedades da camada do QGIS para essas camadas. Corda de projeção:+proj=utm +zone=5 +ellps=GRS80 +datum=NAD83 +units=m +no_defs
  2. Os mapas OSM estão na projeção Mercator. O QGIS tem muitas opções de Mercator, e eu não tinha certeza de qual usar, então adivinhei "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. Para o mapa mundial, o Quantum GIS carregou a projeção a partir do arquivo .prj. Isso é+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

Os mapas projetados devem, mas não se sobrepõem, no ArcMap? fala sobre problemas de zona. Parece o que eu posso estar experimentando. Tentei alterar as bibliotecas da zona UTM 10N para outras zonas (15N e 60N), mas nada mudou.

Alguém tem sugestões para o que tentar a seguir?


2
Se seus dados de Vancouver estão na zona 10TM da UTM, é a string do projeto que você deve ler +proj=utm +zone=10e não +zone=5porque eu acho que a zona 5 está muito perto de algum lugar no Reino Unido (começa com um G ...) e é por isso que seus dados estão chegando na França
SaultDon

Você pode acessar lrdw.ca e fazer o download do shapefile da UTM Zones com uma projeção BC Albers. Entre como convidado e vá para a página 'Principais produtos'; está perto do topo. Puxe esse arquivo para o seu projeto QGIS e veja se seus dois conjuntos de dados se enquadram na zona 10N - usando BC Albers como a transformação CRS "on the fly" (deve ser ativada primeiro Shift + Ctrl + P e marque a caixa necessária) .
precisa saber é o seguinte

Respostas:


11

O que realmente precisamos descobrir aqui são as projeções dos seus dados de entrada.

  1. Dados da cidade

    O Shapefile vem com um arquivo .prj que é perfeitamente interpretado pelo QGIS. Não altere as configurações da camada!

  2. Dados OSM

    Se você estiver usando dados vetoriais OSM, eles estão nas coordenadas geográficas WGS84 (EPSG: 4326). Somente os blocos de varredura estão no Web Mercator (EPSG: 900913).

  3. Abrir dados do Street Style

    Se o QGIS reconheceu o conteúdo .prj, deixe-o assim.

Por fim, escolha um sistema de coordenadas do projeto (sugiro o UTM) e ative a reprojeção on-the-fly como o @Mike_Toews já explicou.


6

Você precisa reprojetar as diferentes camadas para a mesma projeção 'on the fly'. Vamos para:

Configurações> Propriedades do projeto ...> Sistema de referência de coordenadas (CRS)

e verificar em Ativar 'on the fly' transformação CRS


1
Um pequeno esclarecimento: Usar 'on the fly' não será reprojetado no verdadeiro sentido, pois os arquivos reterão sua projeção de origem, mas serão 'desenhados' dinamicamente com uma projeção diferente que você especificar. Esta resposta ajudará você a começar a solucionar problemas com o QGIS.
precisa saber é o seguinte

1
ignorar a maioria do meu comentário acima porque sua formulação é realmente correto, seu meu cérebro que não estava = P
SaultDon

1

Acabei de ter esse problema ao adicionar dois rasters, um no UTM e outro no Geographic. Eu segui as instruções permitindo a projeção instantânea, mas nada aconteceu. Removai as camadas e as adicionei novamente, e o problema foi resolvido


Obrigado John. Eu tive problemas para sobrepor duas camadas de vetor e tentei remover e adicionar novamente as camadas (já que nada mais funcionava). Isso não funcionou por si só, mas quando abreviei o nome do shapefile no qual estava salvando, estava tudo bem.

Por favor, não adicione "obrigado" como resposta. Depois de ter reputação suficiente , você poderá votar nas perguntas e respostas que achou úteis. - Da avaliação
Midavalo
Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.