經度關閉180度與cartopy Orthographic和RotatedPole

Question

我試圖從北極沿着素數子午線（經度= 0）放置藍色點，而是看到沿着日期線（經度= 180）的點。 代碼：

#!/usr/bin/env python 
import matplotlib.pyplot as plt 
import cartopy.crs as ccrs 
Lon = 0 
ax = plt.axes(projection=ccrs.Orthographic(central_latitude=70, 
              central_longitude=Lon)) 
ax.set_global() 
vector_crs = ccrs.RotatedPole(pole_latitude=90, pole_longitude=Lon) 
ax.plot([Lon, Lon, Lon, Lon, Lon, Lon], # longitude 
     [ 90, 80, 70, 60, 50, 40], # latitude 
     'bo', markersize=5, transform=vector_crs) 
ax.stock_img() 
plt.show() 

可能是一些相關的變換，但我還沒有想出什麼。 Cartopy版本0.14.2，Python 3.6。

Answer 1

我認爲這個問題是來自變換，您已經定義的投影：

ax = plt.axes(projection=vector_crs) 
ax.coastlines() 
plt.show() 

注意，這個簡單的海岸線在變換投影看起來像一個長方情節繪製中心經度爲180 °。考慮到這一點，讓我們來看看爲了也儘量簡化你繪製到地圖上用長方投影的陰謀策劃你的樣本數據點：

ax = plt.axes(projection=ccrs.PlateCarree()) 
ax.plot([Lon, Lon, Lon, Lon, Lon, Lon], 
     [ 90, 80, 70, 60, 50, 40], 
     'bo', markersize=5, transform=vector_crs) 
ax.stock_img() 
plt.show() 

由於用你的例子，這些點不會出現在我們預期的地方。最後，讓我們嘗試使用長方投影爲他們繪製到一個正射時，變換你的觀點：

ax = plt.axes(projection=ccrs.Orthographic(central_latitude=70, 
              central_longitude=Lon)) 
ax.set_global() 
ax.plot([Lon, Lon, Lon, Lon, Lon, Lon], 
     [ 90, 80, 70, 60, 50, 40], 
     'bo', markersize=5, transform=ccrs.PlateCarree()) 
ax.stock_img() 
plt.show() 

這似乎提供更多您所尋找的情節。