Блог

Я пытаюсь создавать отдельные стрелки за раз (т. Е. Просто Показывать вектор ветра в одном месте), используя plt.quiver на карте. Однако при этом используется (с помощью ccrs.Преобразование PlateCarree ()) Я получаю проблему, когда какая-то подфункция пытается получить количество затемнений массива 0-dim / int.

Вот несколько минимальных рабочих примеров:

 from matplotlib import pyplot as plt
import cartopy 
import cartopy.crs as ccrs
%matplotlib inline

#### Works: simple quiver plot
plt.quiver(0,0,1,-1)


#### Works: simple quiver plot on a map axis, without specifying transform (this obviously would get the location wrong, but just to show what works/doesn't)
plt.subplot(projection=ccrs.EckertIV())
plt.quiver(0,0,1,-1)

#### Works: simple quiver plot on a map axis with some other transform
plt.subplot(projection=ccrs.EckertIV())
plt.quiver(0,0,1,-1,transform=ccrs.EckertIV())

#### Doesn't work: simple quiver plot on a map axis with the PlateCarree transform
plt.subplot(projection=ccrs.EckertIV())
plt.quiver(0,0,1,-1,transform=ccrs.PlateCarree())
 

Этот последний вызов выдает следующую ошибку:

 ---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
<ipython-input-93-dec36b09994d> in <module>
      1 plt.subplot(projection=ccrs.EckertIV())
----> 2 plt.quiver(0,0,1,-1,transform=ccrs.PlateCarree())

~/.conda/envs/climate1/lib/python3.7/site-packages/matplotlib/pyplot.py in quiver(data, *args, **kw)
   2791 def quiver(*args, data=None, **kw):
   2792     __ret = gca().quiver(
-> 2793         *args, **({"data": data} if data is not None else {}), **kw)
   2794     sci(__ret)
   2795     return __ret

~/.conda/envs/climate1/lib/python3.7/site-packages/cartopy/mpl/geoaxes.py in wrapper(self, *args, **kwargs)
    308 
    309         kwargs['transform'] = transform
--> 310         return func(self, *args, **kwargs)
    311     return wrapper
    312 

~/.conda/envs/climate1/lib/python3.7/site-packages/cartopy/mpl/geoaxes.py in quiver(self, x, y, u, v, *args, **kwargs)
   1837             # Transform the vectors if the projection is not the same as the
   1838             # data transform.
-> 1839             if (x.ndim == 1 and y.ndim == 1) and (x.shape != u.shape):
   1840                 x, y = np.meshgrid(x, y)
   1841             u, v = self.projection.transform_vectors(t, x, y, u, v)

AttributeError: 'int' object has no attribute 'ndim'
 

Обходной путь прямо сейчас состоит в том, чтобы просто дважды нанести одну и ту же стрелку поверх самой себя, например.

 plt.subplot(projection=ccrs.EckertIV())
plt.quiver(np.array([0,0]),np.array([0,0]),np.array([1,1]),np.array([-1,-1]),transform=ccrs.PlateCarree())
 

но мне было интересно, не упустил ли я чего-то с точки зрения того, чтобы просто заставить это нормально работать.

Я использую cartopy 0.18.0 и matplotlib 3.2.1.

Заранее спасибо за любые советы!

Иногда необходимо строго придерживаться подписи вызова.

Подпись вызова:

колчан ([X, Y], U, V, [C], ** kw)

• X, Y: 1D или 2D-массив, подобный

Если все числа X, Y, U, V помещены в форму массива:

 plt.quiver(np.array([0]),np.array([0]), 
           np.array([1]),np.array([-1]), transform=ccrs.PlateCarree())
 

Это будет работать.