Можно ли использовать проекцию гомодозина Гуда в QGIS?

11

Можно ли использовать проекцию гомодозина Гуда в QGIS? Как QGIS, так и GRASS, похоже, испытывают проблемы с перерывами, объединяя более высокие широты. Я пробовал + proj = goode + ellps =phere + datum = wgs84 как пользовательский CRS.

mattwigway
источник

Ответы:

0

Ранее задавался вопрос о прерывистых проекциях: прерывистые или прерывистые картографические проекции?

Был сделан вывод, что proj4 (который используется для обработки проекций в QGIS и GRASS) не поддерживает прерванные проекции.

Подземье
источник
Я должен был упомянуть билет на трекере Proj.4 для поддержки Goode Homolosine (который был создан после того, как был задан этот вопрос): trac.osgeo.org/proj/ticket/106 . Похоже, у меня не такой новый проект, как я думал; Я попробую обновить мой proj и посмотрю, как оно пойдет. Извините за двойной вопрос.
Mattwigway
1
Все проекции прерываются ... то есть встречная область является конечным пространством, а не бесконечным;)
Луис де Соуза
6

С QGIS Lisboa теперь есть базовая поддержка проекции гомодозина Гуда.

Я загрузил данные IGH с http://edc2.usgs.gov/glcc/fao/data_specifications.php

и настроить этот пользователь CRS:

+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

для Северной Америки.

lon_0 и x_0 взяты из приложения 03 к этому документу: http://proceedings.esri.com/library/userconf/proc98/proceed/TO850/PAP844/P844.HTM

Наложая границы естественного мира Земли, я получаю эту картину:

введите описание изображения здесь

Хорошо вписывается в Северной Америке (около 20 км не на месте) и искажается в других регионах.


РЕДАКТИРОВАТЬ

Только что заметил, что proj знает прерванный гомолозин Гуда:

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

дает мне эту картину границ мира Естественных Земель и векторную сетку с 10 степенями:

введите описание изображения здесь

Я не уверен, какие параметры кроме +proj=ighнеобходимы больше.

Andrej
источник
Используя приведенное выше определение, я получаю сдвиг пикселей. Для меня определение +proj=igh +towgs84=0,0,0работало очень хорошо
Томас
3

QGIS 2.18 почти полностью поддерживает гомолозин для растров. Вам просто нужно знать о проблеме # 959 : выполнить деформацию растра напрямую с помощью GDAL, используя CHECK_WITH_INVERT_PROJпараметры, и все в порядке.

Относительно векторов вещи все еще сложны. Некоторые из ошибок, которые AndreJ показывает в своем ответе преобладают. Причина проста: QGis не знает, что это картографическая программа, и рассматривает совместную область любой проекции как конечное пространство.

Это легко исправить: пересекайте любой векторный слой, импортированный в QGis, с гомодозиновым доменом. Некоторое время назад я сделал доступным вектор с раздельными областями, определяемый с интервалами в один градус долготы и широты.

DOI

введите описание изображения здесь

Луис де Соуза
источник