Proj4 Oblique Mercator дает странные результаты, когда альфа составляет от 90 до 270 включительно

Question

У меня есть следующая строка проекции proj4, определяющая проекцию Oblique Mercator:

+ proj = omerc + lat_0 =45 + lonc = -80 + x_0 = 0 + y_0 = 0 + альфа =45 + гамма =0 +k_0=1

Для моей входной координаты я использовал 0, 0. Когда я изменяю угол альфа от 45 до диапазона от 90 до 270 включительно, проецируемая точка меняется с ожидаемых 80 Вт, 45N на точку, которая в восточном направлении распространяется вокруг Земли в то же время широта до альфа> 270, где я снова получаю ожидаемые результаты.

invproj +proj=omerc +lat_0=45 +lonc=-80 +x_0=0 +y_0=0 +alpha=45 +gamma=0 +k_0=1
0 0
80dW    45dN  (this is expected)

invproj +proj=omerc +lat_0=45 +lonc=-80 +x_0=0 +y_0=0 +alpha=89 +gamma=0 +k_0=1
0 0
80dW    45dN  (this is expected)

invproj +proj=omerc +lat_0=45 +lonc=-80 +x_0=0 +y_0=0 +alpha=91 +gamma=0 +k_0=1
0 0
77d10'18.699"W  45dN  (this is not expected)

Я не уверен, почему я должен получать эти разные результаты, в качестве обходного пути, я мог бы повернуть угол от 91 до 271, но я думаю, что мне может понадобиться также изменить U, V в проецируемой системе, было бы неплохо если Proj4 сделал это.

• Указание альфы от 90-270 не поддерживается?
• Или мне нужно что-то сделать с моими параметрами?

Изменить: я создал проблему в Github для этого по адресу: https://github.com/OSGeo/proj.4/issues/331

Answer 1

Похоже, что произошла регрессия между версиями 4.7 и 4.8 Proj4, которая привела к ошибке для этого сценария, я обнаружил следующее в архивах списков рассылки:

http://lists.maptools.org/pipermail/proj/2012-June/006331.html

Микаэль Риттри написал:

Не то, что я ожидал. Насколько я понимаю, эти странные результаты появляются, когда альфа южная (в закрытом интервале от 90 до 270). Omerc Proj 4.7.0 не вел себя таким образом.

Это совпадает с моими наблюдениями. Мое решение состояло в том, чтобы полностью изменить альфу и добавить слой предварительной / последующей обработки, который инвертирует знак значений восточного и северного направлений.