Proj4 Oblique Mercator 在 alpha 介于 90 和 270 之间时给出奇怪的结果

Question

我有以下定义斜墨卡托投影的 proj4 投影字符串：

+proj=omerc +lat_0=45 +lonc=-80 +x_0=0 +y_0=0 +alpha=45 +gamma=0 +k_0=1

对于我的输入坐标，我使用的是 0、0。当我将 alpha 角从 45 更改为 90 到 270（含）范围时，投影点从预期的 80W、45N 变为向东移动的点地球在同一纬度，直到 alpha > 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 之间的 alpha？
或者，我需要对我的参数做些什么吗？

编辑：我在 Github 上为此创建了一个问题：https://github.com/OSGeo/proj.4/issues/331

Answer 1

看起来 Proj4 版本 4.7 和 4.8 之间存在回归，为这种情况引入了错误，我在邮件列表档案中发现了以下内容：

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

Mikael Rittri 写道：

Not what I would expect. As far as I can understand, these weird results appear when alpha is southish (in the closed interval 90 to 270). The omerc of Proj 4.7.0 did not behave this way.

这与我的观察不谋而合。我的解决方案是反转 alpha 并添加一个 pre/post 处理层来反转东向和北向值的符号。

Proj4 Oblique Mercator 在 alpha 介于 90 和 270 之间时给出奇怪的结果

Proj4 Oblique Mercator giving strange results when alpha is between 90 and 270 inclusive

proj

mercator