我正在努力研究 matplotlib 的 basemap API。我想做的是在圆柱投影上绘制一个非常小的区域(任一方向几公里)的纬度和经度数据。
问题是,我无法理解llcrnrlon、llcrnrlat、urcrnrlon 和urcrnrlat 参数Basemap 构造函数工作。
据我了解,llcrnrlon 是最西经度,llcrnrlat 是最南纬度,urcrnrlon 和urcrnrlat分别是最东端和最北端的经度和纬度。
在所有这些情况下,给定一组坐标,(数值上)最小的经度是最西的,最小的纬度是最南的,反之亦然。这种理解是否正确?
我可以通过在底层 Axes 对象上设置 xlim 和 ylim 来使绘图工作,但在 basemap 构造函数中使用相同的值似乎会将我的数据推离绘图。
编辑:请参阅下面的代码以重现问题:
from matplotlib import pyplot
from mpl_toolkits import basemap
import random
lat_bounds = 52.063443, 52.072587
long_bounds = 1.010408, 1.024502
ax1 = pyplot.subplot(121)
ax2 = pyplot.subplot(122)
ax1.set_title('With ll* rr*')
ax2.set_title('With default values')
my_map1 = basemap.Basemap(projection='cyl', llcrnrlat=lat_bounds[0], llcrnrlon=long_bounds[0],
urcrnrlat=lat_bounds[1], urcrnrlon=long_bounds[1], ax=ax1)
my_map2 = basemap.Basemap(projection='cyl', ax=ax2)
data_lats = [random.uniform(*lat_bounds) for i in xrange(50)]
data_lons = [random.uniform(*long_bounds) for i in xrange(50)]
my_map1.plot(data_lats, data_lons)
my_map2.plot(data_lats, data_lons)
pyplot.show()
在下图中,右侧图像是使用 Basemap(projection='cyl') 制作的,左侧图像是使用 Basemap(projection=' cyl',llcrnrlat=lat_bounds[0],llcrnrlon=long_bounds[0],urcrnrlat=lat_bounds[1],urcrnrlon=long_bounds[1],ax=ax1)
请注意右侧图像中的圆点,当使用 matplotlib 工具栏缩放时,它会变成第二张图像。
请您参考如下方法:
示例代码中的问题是您以错误的方式将参数传递给 Basemap.plot()。 Basemap.plot 的参数与 matplotlib.pyplot.plot 的参数完全相同,即:
plot(x,y,*args,**kwargs)
在圆柱坐标系中,经度是 x 坐标,纬度是 y 坐标,因此您应该执行 mymap1.plot(data_lons, data_lats)。它在你的第二个例子中似乎起作用的原因是经度 ~52 和纬度 ~1 是有意义的。这些点已绘制,但远离您的域。如果您将 ax1 的窗口平移得足够远,您就会看到它们(这与执行 ax.set_xlim(lat_bounds) 和 ax.set_ylim( lon_bounds)).
