Использование базовой карты Python etopo() генерирует IndexError

Question

Я пытаюсь работать с базой. Я установил matplotlib, baseemap и Pillow успешно и могу строить основные сюжеты. Код ниже работает для меня.

from mpl_toolkits.basemap import Basemap 
from PIL import Image 
import matplotlib.pyplot as plt 

m = Basemap(projection = 'mill', 
     llcrnrlat = -90, 
     llcrnrlon = -180, 
     urcrnrlat = 90, 
     urcrnrlon = 180, 
     resolution = 'c') 

m.drawcoastlines() 
#m.etopo() 

plt.show() 

Однако раскомментирование команды 'etopo()' приводит к следующей ошибке.

Traceback (most recent call last): 

    File "C:\Users\xxx\Desktop\helper.py", line 13, in <module> m.etopo() 

    File "C:\Python35-32\lib\site-packages\mpl_toolkits\basemap\__init__.py", line 4061, in etopo return self.warpimage(image='etopo',scale=scale,**kwargs) 

    File "C:\Python35-32\lib\site-packages\mpl_toolkits\basemap\__init__.py", line 4167, in warpimage np.concatenate((self._bm_lons,self._bm_lons+360),1) 

IndexError: axis 1 out of bounds [0, 1) 

Я попытался переустановить и обновить все соответствующие пакеты, но это не сработало. Я также не могу найти ничего об этой ошибке для этой ситуации. Я на ОС Windows 8.1 с 32 битным Python 3.5.1

Answer 1

Если вы только тестирование Basemap и .etopo() функциональность, я обновил сценарий:

from mpl_toolkits.basemap import Basemap 
from PIL import Image 
import matplotlib.pyplot as plt 

# make a miller cylindrical projection with defaults 
m = Basemap('mill') 

m.drawcoastlines() 
# now displays topo relief image 
m.etopo() 

plt.show() 

Answer 2

Это, как представляется, incompatibilty между базовой картой и недавно NumPy версии. Предлагаемое решение возвращается к проекции по умолчанию, поскольку первое ключевое слово - не проекция, а llcrnrlat.

Я не понял, что именно вызывает IndexError, однако возвращается к более старой версии numpy. Кажется, что numpy 1.9.2 является последней рабочей версией. Возможно, что-то было изменено командой concatenate.

Edit: я изменил линию 4167 в BaseMap/INIT .py, чтобы получить тот же результат, чем это должно быть с более старыми версиями Numpy. Это, кажется, сделать трюк:

Old: np.concatenate((self._bm_lons,self._bm_lons+360),1) 
New: np.concatenate((self._bm_lons.T,self._bm_lons.T+360)).T 

конкатенации вдоль 1D массивов и ненулевой оси, теперь осуждается, именно поэтому исключение.