将立体图限制为底图中的四分之一或一部分

时间:2017-07-31 17:22:24

标签: matplotlib matplotlib-basemap map-projections

我一直在搜索stackoverflow和整个互联网,但只是找不到我的问题的解决方案。

我想将一个立体图缩小到具有ge绘制的实际有效点的区域。这意味着"切断"我附着的图像上的白色多余。 Basemap的文档并没有真正帮助。

使用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获取此结果的方法。谢谢!

0 个答案:

没有答案