北极立体投影不起作用

Question

我正在尝试使用从北极偏移的底图创建立体图，但东西方向显然是相反的。这是我的实现错误还是错误？

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np

m = Basemap(projection='stere',
        lat_0=90, lon_0=270,  lat_ts=(90.+35.)/2.,
        llcrnrlon=150,urcrnrlon=-60,llcrnrlat=50,urcrnrlat=50)

m.drawmeridians(np.arange(0,360,30),labels=[1,1,1,0])
m.drawparallels(np.arange(-90,90,5))
m.drawcoastlines()
m.shadedrelief()
plt.show()

结果如下： result from script

如何重现以下地图（以偏移为中心并旋转？）

Restricted map

Answer 1

使用方位角类型的地图投影总是需要一组适当的参数才能获得良好的结果。在这种情况下，Stereographic projection 以北极为中心，其适当的参数不是您经常使用的 PlateCaree 投影实现时通常使用的参数。这是您可以尝试的工作代码。

# Stereographic projection coverage
#   should be specified less than half of a hemisphere

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np

m = Basemap(projection='stere', resolution='c',
        lat_0=90, lon_0=270,  lat_ts=(90.+35.)/2., width=15000000, height=10000000)
# (width, height) is the plot extents in meters

m.drawmeridians(np.arange(0, 360, 30), labels=[1,1,1,0])
m.drawparallels(np.arange(0, 90, 10), labels=[0,0,0,1])
m.drawcoastlines()
m.shadedrelief()
plt.show()

结果图（地图 1）：

要让世界的其他部分进入绘图区域，可以通过重新居中地图来实现。

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np

# projection center point
lon0 = 180
lat0 = 60

m = Basemap(projection='stere', resolution='c',
        lat_0=lat0, lon_0=lon0,  lat_ts=lat0, width=15000000, height=10000000)

m.drawmeridians(np.arange(0, 360, 30), labels=[1,0,0,1])  # left, right, top, bottom
m.drawparallels(np.arange(0, 90, 10), labels=[0,1,1,0])
m.drawcoastlines()
m.shadedrelief()
plt.show()

输出图（地图 2）：

通过指定适当的llcrnrlon、urcrnrlon、llcrnrlat、urcrnrlat值，在Basemap()中可以根据需要得到地图范围。这是 OP 要求的另一个情节示例。

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np

fig = plt.figure(figsize=(8,8))
m = Basemap(projection='stere', resolution='c',
        lat_0=90, lon_0=-90,  lat_ts=(90.+35.)/2.,
        llcrnrlon=-142, urcrnrlon=78, llcrnrlat=19, urcrnrlat=45)

m.drawmeridians(np.arange(0, 360, 30), labels=[1,0,1,0])  # left, right, top, bottom
m.drawparallels(np.arange(0, 90, 10), labels=[0,1,0,1])
m.drawcoastlines()
m.shadedrelief()
plt.show()

结果图（地图 3）：

北极立体投影不起作用

North polar stereographic projection is not working

map-projections

python-3.x

matplotlib-basemap