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