python - 極平射図法

Question

cartopy の最新バージョンを使用して極平射図上に等高線をプロットすることは可能ですか? 私は自分でそれを解決するのに苦労しているので、これがどのように行われるかの例を見たいです!

Answer 1

立体射影はいくつかの頭痛の種を引き起こしています。おそらく、cartopy のポリゴン変換コードで最も問題を引き起こしたのは射影です。

次の例は、カートピーを使用して極平射図を作成する方法を示しています。注: このコードを使用しても、サンプル データの解像度を微調整すると、プロットが実際にレンダリングされるまでに約 30 分かかることがわかります (これはバグであり、遅かれ早かれ分類する必要があります)。

import cartopy.crs as ccrs
import matplotlib.pyplot as plt

from cartopy.examples.waves import sample_data

ax = plt.axes(projection=ccrs.NorthPolarStereo())

x, y, z = sample_data((100, 200))
cs = ax.contourf(x, y, z, 50,
                 transform=ccrs.PlateCarree(),
                 cmap='gist_ncar')
ax.coastlines()

# without the set_global, currently, the plot is tiny because the limits
# are being erroneously being set (opened issue for that)
ax.set_global()

plt.show()

うまくいけば、カートピーで極平射図等高線図を作成する方法がわかります。データに問題がある場合は、「Geometry transforms」というタグが付けられた未解決の問題を見て、類似した問題が発生しているかどうかを確認してください。そうでない場合は、問題を開いて調査してください。

注: この回答は cartopy v0.5.x (つまり、v0.5 リリースの直前) に関するものであり、ここで言及されているバグの多くは将来のリリースで解消されることを願っています。

それが役立つことを願って、