Определение Winkel Tripel в proj 4.8.0?

14

У меня QGIS с установленным proj 4.8.0. Я хотел бы добавить проекцию Winkel Tripel, которая, по-видимому, определяется как

+proj=wintri

Но эта строка определения не работает. Это ошибка в моей версии proj, или какая строка определения правильная?

qgis coordinate-system proj Подземье
источник
1
У меня работает: эхо "34 -34" | proj "+ proj = wintri"
Пол Рэмси,
1
Да, это что-то эхо. Но если я определю пользовательский CRS с этой строкой, в QGIS ничего не будет перепроектировано. Угадай, мне нужны еще параметры?
Подземье
7
@Paul Ramsey: Кажется, что хотя проекция Winkel Tripel определена в projбиблиотеке и может вызываться из командной строки, ее нельзя использовать в качестве пользовательского CRS в QGIS, поскольку в projбиблиотеке нет обратного преобразования (см. Списки .maptools.org / pipermail / proj / 2012-May / 006302.html ). И кажется, что обратное преобразование не тривиально (и всегда будет приближением ...).
Джейк
8
Я думаю, что вы должны написать это в ответ, так как это ... ответ. Инструменты командной строки также отказываются выполнять обратную проекцию.
Пол Рэмси
1
Как упоминалось выше коллегами, если у вас также с proj + proj = aitoff есть те же проблемы, то это вопрос отсутствия обратной проекции, см .: (+ proj = aitoff; Rel. 4.7.1, 23 сентября 2009 г. <invproj>: обратная проекция недоступна) и (invproj + proj = wintri; Rel. 4.7.1, 23 сентября 2009 г. <invproj>: обратная проекция недоступна). Странно, но полезно знать
huckfinn

Ответы:

4

Подводя итог обсуждению выше:

Хотя проекция Winkel Tripel определена в библиотеке proj и может вызываться из командной строки, ее нельзя использовать в качестве пользовательского CRS в QGIS, поскольку в библиотеке proj нет обратного преобразования .

Запрос на расширение для добавления этой функциональности был закрыт, поскольку кажется, что обратное преобразование не тривиально.

Подземье
источник
Может ли это быть по-другому сейчас? Я с pyproj==1.9.5.1установленным, могу сделать: pp = Proj(proj='wintri', lon_0=47.); pp(*pp(0,0), inverse=True)что возвращает (0,0)(к точности машины), то есть я могу отменить прямое преобразование с inverse=Trueключевым словом. Поскольку этот pyproj, похоже, имеет обратную сторону, может ли QGIS поддерживать это сейчас?
Ахмед Фасих