Используя cartopy
, я пытаюсь построить заполненные контуры в проекции, центрированной по ненулевой долготе:
from pylab import *
from cartopy import crs
xi = linspace(-180, 180, 37)
yi = linspace(-90, 90, 19)
XI,YI = meshgrid(xi, yi)
ZI = cos(pi*XI/180)*cos(pi*YI/180)
ax = axes( projection = crs.PlateCarree(central_longitude = 90) )
contourf(
XI, YI, ZI, 10,
transform=crs.PlateCarree(),
)
ax.coastlines()
savefig('mwe.pdf')
Это дает следующий график, который (в некоторых масштабах) отображаеттонкая вертикальная линия (проходящая через центр более темной области), соответствующая разрыву в значениях долготы (180 градусов восточной долготы или западной долготы).Вероятно, это связано с тем, что формы, выводимые contourf
, неправильно соединены по всему разрыву.
Есть ли способ исправить это?Или же я могу произвольно редактировать базовые координаты, используемые coastlines()
для имитации проекции со смещением по долготе?