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


14

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

+proj=wintri

Но эта строка определения не работает. Это ошибка в моей версии 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 поддерживать это сейчас?
Ахмед Фасих
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.