2

我有以下定义斜墨卡托投影的 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

4

1 回答 1

1

看起来 Proj4 版本 4.7 和 4.8 之间存在回归,这为这种情况引入了一个错误,我在邮件列表存档中找到了以下内容:

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

米凯尔·里特里写道:

不是我所期望的。据我所知,当 alpha 偏南时(在 90 到 270 的封闭区间内),就会出现这些奇怪的结果。Proj 4.7.0 的 omerc 并没有这样表现。

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

于 2015-12-24T01:49:39.290 回答