Czy w QGIS można zastosować projekcję Goode Homolosine?

11

Czy w QGIS można zastosować projekcję Goode Homolosine? Zarówno QGIS, jak i GRASS wydają się mieć problemy z przerwami, miażdżąc razem wyższe szerokości geograficzne. Próbowałem + proj = goode + ellps = kula + dane = wgs84 jako niestandardowy CRS.

mattwigway
źródło

Odpowiedzi:

0

Wcześniej pojawiło się pytanie o przerywane projekcje: Nieciągłe lub przerywane projekcje map?

Wniosek był taki, że proj4 (który jest używany do obsługi projekcji w QGIS i GRASS) nie obsługuje przerwanych projekcji.

podmrok
źródło
Powinienem wspomnieć o biletach w module śledzącym Proj.4 na wsparcie Goode Homolosine (który został stworzony po zadaniu tego pytania): trac.osgeo.org/proj/ticket/106 . Wygląda jednak na to, że nie mam tak nowego proj, jak myślałem; Spróbuję zaktualizować mojego proj i zobaczę, jak idzie. Przepraszamy za duplikat pytania.
mattwigway
1
Wszystkie projekcje są przerywane ... tzn. Kontr-domena jest przestrzenią skończoną, a nie nieskończoną;)
Luís de Sousa
6

Dzięki QGIS Lisboa istnieje teraz podstawowa obsługa projekcji Goode Homolosine.

Pobrałem dane IGH z http://edc2.usgs.gov/glcc/fao/data_specifications.php

i skonfiguruj CRS tego użytkownika:

+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

dla Ameryki Północnej.

lon_0 i x_0 pochodzą z załącznika 03 do tego dokumentu: http://proceedings.esri.com/library/userconf/proc98/proceed/TO850/PAP844/P844.HTM

Nakładając naturalne granice świata Ziemi_m, otrzymuję ten obraz:

wprowadź opis zdjęcia tutaj

Pasuje dobrze w Ameryce Północnej (około 20 km źle umieszczony) i zniekształcony w innych regionach.


EDYTOWAĆ

Właśnie zauważyłem, że proj zna już przerwaną goode homolosine:

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

daje mi ten obraz granic świata Natural Earths i 10-stopniowej siatki wektorowej:

wprowadź opis zdjęcia tutaj

Nie jestem pewien, które parametry oprócz +proj=ighsą konieczne.

AndreJ
źródło
Korzystając z powyższej definicji, uzyskuję przesunięcie pikseli. Dla mnie definicja +proj=igh +towgs84=0,0,0działała dobrze
Thomas
3

QGIS 2.18 jest bliski pełnego wsparcia homolozyny dla rastrów. Musisz tylko zdawać sobie sprawę z problemu # 959 : wykonaj dowolne wypaczenie rastra bezpośrednio za pomocą GDAL, używając CHECK_WITH_INVERT_PROJparametrów i jesteś dobry.

Jeśli chodzi o wektory, rzeczy są nadal skomplikowane. Niektóre błędy, które AndreJ pokazuje w swojej odpowiedzi, przeważają. Przyczyna jest prosta: QGis nie zdaje sobie sprawy, że jest to program kartograficzny i uważa, że ​​domena każdej projekcji jest skończoną przestrzenią.

Ma to prostą poprawkę: przecina każdą warstwę wektorową zaimportowaną do QGis z domeną homolozynową. Jakiś czas temu udostępniłem dyskretny wektor w domenie wspólnej zdefiniowany w odstępach długości i szerokości geograficznej jednego stopnia.

DOI

wprowadź opis zdjęcia tutaj

Luís de Sousa
źródło