Полярный стереографический сюжет

Стереографическая проекция вызывает несколько головных болей и, вероятно, является той проекцией, которая вызвала больше всего проблем для кода преобразования многоугольника в 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()

Надеюсь, это покажет вам, как строить полярный стереографический контур в картографии. Если у вас возникли проблемы с данными, просмотрите открытые проблемы с меткой " Преобразование геометрии " и посмотрите, получаете ли вы что-то подобное, а если нет, продолжайте и откройте проблему, и мы сможем ее решить.

Примечание. Этот ответ относится к cartopy v0.5.x (т.е. непосредственно перед выпуском v0.5), и мы надеемся, что многие из упомянутых здесь ошибок будут исправлены в будущих выпусках.

Надеюсь, это поможет,