Проблема с несовпадением шейп-файлов GSHHS и Natural Earth.Python

Программы на Python
Anonymous
Проблема с несовпадением шейп-файлов GSHHS и Natural Earth.

Сообщение Anonymous »

Я пытаюсь построить данные по региону Американского Самоа, но у меня возникли проблемы с береговыми линиями. Обычно я использую данные с более высоким разрешением из Глобальной самосогласованной иерархической географической базы данных высокого разрешения, но береговая линия выглядит смещенной по сравнению с данными, которые я рисую, а также из других наборов данных, таких как NaturalEarth и карты Google.
Кто-нибудь знает источник этого смещения? Это проблема с тем, как я рисую береговые линии, проблема с картографической проекцией, которую я не учел, проблема с набором данных GSHHS или, может быть, проблема с его реализацией в Cartopy?
Ниже приведен пример кода, иллюстрирующий мою проблему со смещением острова Американское Самоа на запад в GSHHS (синий).

Код: Выделить всё

import matplotlib.pyplot as plt
import cartopy

proj = cartopy.crs.PlateCarree()
fig, ax = plt.subplots(subplot_kw=dict(projection=proj), figsize=(12,12))

ax.add_feature(cartopy.feature.GSHHSFeature(scale='f'), facecolor="blue")

ax.add_feature(cartopy.feature.NaturalEarthFeature("physical", "land", "10m"),
ec="red", fc="yellow", lw=2, alpha=0.4)

lon_min_ = -170.9
lon_max_ = -170.55
lat_min_ = -14.4
lat_max_ = -14.2

ax.set_extent([lon_min_, lon_max_, lat_min_, lat_max_],crs=proj)

ax.gridlines(draw_labels=True,alpha=.4,linewidth=2, color='black', linestyle='--')

plt.show()
Изображение


Подробнее здесь: https://stackoverflow.com/questions/789 ... shapefiles

Вернуться в «Python»