私は、例えば立体投影を使って私のデータをプロットしたいと思います。しかし、私がプロットしようとすると、大陸を除いて何もプロットされません。Python - pcolormeshとベースマップでプロットする
longrid_t
は経度である
plt.pcolormesh(longrid_t, latgrid_t,totvart_t)
、latgrid_t
は緯度とtotvart_t
私はプロットしたいデータである:
説明するが、これは、私のデータプロットの図です。
ここで図はplt.pcolormesh(longrid_t, latgrid_t,totvart_t)
として(ベースマップなし)のみpcolormeshプロットである:
ここで、Iは、立体射影を使用してこれらのデータをプロットしてみました:
map = Basemap(width=10000000,height=6000000, resolution='l',projection='stere',
lat_0=40.,lon_0=-40.)
map.fillcontinents(color='darkgrey',lake_color='darkgrey')
map.drawcoastlines()
map.drawparallels(arange(-90, 91., 10.), labels=[1, 0, 0, 1],
dashes=[1, 1], linewidth=0.25, color='0.5')
map.drawmeridians(arange(-180., 181., 10.), labels=[1, 0, 0, 1],
dashes=[1, 1], linewidth=0.25, color='0.5')
map.pcolormesh(longrid_t, latgrid_t, totvart_t, cmap='jet')
map.colorbar()
しかし、この立体投影の結果は私が期待したものではありません:
私は他の予測も試しました。また、pcolormeshの代わりにpcolorとcontourfを使用することもできます。しかし、それは私に同じ結果をもたらします。
アイデアはありますか?あなたは、投影のxy座標にあなたの座標データを変換し、変換後の座標にプロットする必要が
'map.pcolormesh'の呼び出しで' latlon = True' kwargを使ってみましたか? – farenorth