我一直在搜索stackoverflow和整个互联网,但只是找不到我的问题的解决方案。
我想将一个立体图缩小到具有ge绘制的实际有效点的区域。这意味着"切断"我附着的图像上的白色多余。 Basemap的文档并没有真正帮助。
我用来绘制图形的代码如下:
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap, cm
field = np.ones(shape=(20,80))
plt.figure(figsize=(12, 12))
lonsize=[-50,30]
latsize=[60,80]
m = Basemap(width=4500000,height=3000000,
resolution='i',projection='stere',boundinglat=45,
lat_0=70,lon_0=-10.)
m.drawmapboundary(fill_color='white')
m.drawparallels(np.arange(-90., 91., 30.),labels=[1,1,1,1])
m.drawmeridians(np.arange(-180., 180., 30.),labels=[1,1,0,1],size=10)
m.fillcontinents(color='grey',lake_color='grey')
lonnew = np.linspace(-50, 30, 80)
latnew = np.linspace(60, 80, 20)
lons, lats = np.meshgrid(lonnew, latnew)
xx, yy = m(lons, lats)
m.pcolormesh(xx, yy, field, cmap='Spectral_r')
plt.show()
我希望有人知道用Basemap获取此结果的方法。谢谢!