É possível usar a projeção Goode Homolosine no QGIS?

11

É possível usar a projeção Goode Homolosine no QGIS? Tanto o QGIS quanto o GRASS parecem ter problemas com as interrupções, esmagando as latitudes mais altas. Eu tentei + proj = goode + ellps = esfera + dado = wgs84 como um CRS personalizado.

mattwigway
fonte

Respostas:

0

Houve uma pergunta sobre projeções interrompidas antes: Projeções de mapas descontínuos ou interrompidos?

A conclusão foi que o proj4 (que é usado para manipulação de projeção no QGIS e GRASS) não suporta projeções interrompidas.

underdark
fonte
Eu deveria ter mencionado um ticket no rastreador do Proj.4 para suporte ao Goode Homolosine (que foi criado após a pergunta): trac.osgeo.org/proj/ticket/106 . Parece que eu não tenho um projeto tão novo quanto eu pensava; Vou tentar atualizar meu proj e ver como vai. Desculpe pela pergunta duplicada.
mattwigway
11
Todas as projeções são interrompidos ... ou seja, o domínio contador é um espaço finito, não infinito;)
Luís de Sousa
6

Com o QGIS Lisboa, agora existe um suporte básico para a projeção de Goode Homolosine.

Eu baixei os dados do IGH em http://edc2.usgs.gov/glcc/fao/data_specifications.php

e configure este CRS de usuário:

+proj=goode +lon_0=-100 +x_0=-11119487.43 +y_0=0 +a=6370997 +b=6370997 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

para a América do Norte.

lon_0 e x_0 são do apêndice 03 deste documento: http://proceedings.esri.com/library/userconf/proc98/proceed/TO850/PAP844/P844.HTM

Sobrepondo os limites naturais do mundo terrestre, recebo esta imagem:

insira a descrição da imagem aqui

Ajuste adequado na América do Norte (cerca de 20 km extraviado) e distorcido em outras regiões.


EDITAR

Só notei que o proj sabe interrompeu goode homolosine agora:

+proj=igh +ellps=sphere +towgs84=0,0,0 +lon_0=100w +x_0=-11119487.43

dá-me esta imagem dos limites mundiais da Terra Natural e uma grade vetorial de 10 graus:

insira a descrição da imagem aqui

Não tenho mais certeza de quais parâmetros +proj=ighsão necessários.

AndreJ
fonte
Usando a definição acima, recebo uma mudança de pixel. Para mim, a definição +proj=igh +towgs84=0,0,0funcionou bem
Thomas
3

O QGIS 2.18 está quase apoiando totalmente o Homolosine para rasters. Você só precisa estar ciente do problema # 959 : executar qualquer distorção de varredura diretamente com o GDAL usando os CHECK_WITH_INVERT_PROJparâmetros e você é bom.

Em relação aos vetores, as coisas ainda são complexas. Alguns dos erros que AndreJ mostra em sua resposta prevalecem. A causa é simples: o QGis não sabe que é um programa cartográfico e considera o co-domínio de qualquer projeção um espaço finito.

Isso tem uma solução simples: intercepte qualquer camada vetorial importada para o QGis com o co-domínio Homolosine. Algum tempo atrás, disponibilizei um vetor de co-domínio discreto, definido em intervalos de um grau de longitude e latitude.

DOI

insira a descrição da imagem aqui

Luís de Sousa
fonte