Python 底图：使用 shadedrelief、bluemarble 或 etopo 时出错（假经度格式？）

Question

我想绘制以太平洋为中心的南半球地图，并使用 python matplotlib 底图绘制一些内容。

除非我尝试使用底图例程 shadedrelief、bluemarble 或 etopo 绘制背景图像，否则一切正常。代码（没有我想在地图上绘制的东西）如下所示：

import numpy as np
from mpl_toolkits.basemap import Basemap
from matplotlib.backends.backend_pdf import PdfPages

latmin = -72.5
latmax = 40.
lonmin = 60.
lonmax = 370.

pp = PdfPages('datamap.pdf')

m = Basemap(projection='merc', llcrnrlat=latmin, urcrnrlat=latmax, llcrnrlon=lonmin, urcrnrlon=lonmax, resolution="c")

m.drawcoastlines(linewidth=0.25)
#m.shadedrelief()

pp.savefig()
pp.close()

当我取消注释 m.shadedrelief() 时，我得到以下信息：

---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
/xyz/datamap.py in <module>()
     32
---> 33 m.shadedrelief()
     34 

/usr/local/lib/python2.7/site-packages/mpl_toolkits/basemap/__init__.pyc in shadedrelief(self, ax, scale, **kwargs)
   3997             return self.warpimage(image='shadedrelief',ax=ax,scale=scale,**kwargs)
   3998         else:
-> 3999             return self.warpimage(image='shadedrelief',scale=scale,**kwargs)
   4000 
   4001     def etopo(self,ax=None,scale=None,**kwargs):

/usr/local/lib/python2.7/site-packages/mpl_toolkits/basemap/__init__.pyc in warpimage(self, image, scale, **kwargs)
   4115                 # any range of longitudes may be plotted on a world map.
   4116                 self._bm_lons = \
-> 4117                 np.concatenate((self._bm_lons,self._bm_lons+360),1)
   4118                 self._bm_rgba = \
   4119                 np.concatenate((self._bm_rgba,self._bm_rgba),1)
IndexError: axis 1 out of bounds [0, 1)

当我选择（出于测试目的）经度也不大于 180 度的较小地图时，一切（包括阴影浮雕）都可以正常工作。这让我假设经度格式的东西在这里不起作用。我尝试了一些东西，但我似乎无法找到一种方法来解决这个问题，同时仍然绘制相同的地图部分

你知道我如何在我的地图背景中绘制阴影浮雕吗？

最佳，

西联

Answer 1

底图/init 代码中似乎存在错误。事实上，如果您查看 github 上的最新版本 [请参阅 https://github.com/matplotlib/basemap/blob/master/lib/mpl_toolkits/basemap/init.py#L4139 以及它与您的错误的第 4117 行有何不同] 您会发现该错误已经解决。这意味着您可以：

# np.concatenate((self._bm_lons,self._bm_lons+360),1)
np.concatenate((self._bm_lons,self._bm_lons+360))

在 init.py 代码中或获取较新版本的底图。两者中的任何一个都可以解决您的问题。

Python 底图：使用 shadedrelief、bluemarble 或 etopo 时出错（假经度格式？）

Python Basemap: Error using shadedrelief, bluemarble or etopo (false longitude format?)

python

matplotlib

matplotlib-basemap