2015-12-15 58 views
1

我有限定斜墨卡託投影以下proj4投影字符串:Proj4斜墨卡託給出奇怪的結果當α90和270之間包容

+凸出= 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變化到在相同緯度地球繞地球東進的點,直到α> 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

回答

1

它看起來像有Proj4版本之間的迴歸介紹了一個bug在此方案中4.7和4.8,我發現了以下在郵件列表歸檔:

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

的Mikael Rittri寫道:

不是我所期望的。據我所知,這些奇怪的結果出現在阿爾法 是南方(在閉區間90至270)。 Proj 4.7.0的操作不像這樣。

這與我的觀察一致。我的解決方案是顛倒alpha,並添加一個前/後處理層,它顛倒東向和北向值的符號。

+0

看起來這種行爲持續存在於我的環境中可用的最新版本「4.9.2」中。 – geotheory