Ситуация: я рисую географические границы США (в частности, озеро Мичиган) с помощью Cartopy, а затем рисую множество полилиний с помощью Matplotlib LineCollection. должным образом. Если я использую projection=ccrs.Miller(), геометрия LineCollection смещается на юго-восток (+ долгота, -широта).
Мой вопрос, тогда:
Правильный ли следующий код?
axes = pyplot.gca() # this is a GeoAxesSubplot
transform = axes.projection
my_lines = matplotlib.collections.LineCollection(
numpy.vstack(paths),
transform=projection
)
axes.add_collection(my_lines)
Мета-вопрос: если это - это правильный способ получить LineCollection в GeoAxes, как мне go отладить это?
Приложение: если я использую Робинсона или Меркатора, мой набор данных вообще не отображается. Ясно, что я не могу провести трансформацию в каком-нибудь важном месте.