Как переключаться между EPSG в Geo pandas (прокладка и добавление фоновой карты) - PullRequest
Новая
       45

Как переключаться между EPSG в Geo pandas (прокладка и добавление фоновой карты)

0 голосов
/ 28 февраля 2020

Я пытаюсь нарисовать координаты "станции" на карте NY C. Я использую Geo pandas (gpd), чтобы создать GeoDataFrame и создать карту NY C из gpd.datasets с EPSG 4326.

Кроме того, я хотел добавить фоновое изображение на мою карту следующим учебник по Geo pandas (https://geopandas.org/gallery/plotting_basemap_background.html). Тем не менее, это работает на EPSG 3857. Есть ли способ сохранить информацию и изменить EPSG?

Вот моя база кода 

# Creating GeoDataFrame
gdf = gpd.GeoDataFrame(stations, geometry=gpd.points_from_xy(stations.Longitude, stations.Latitude))
# Find the NYC Borough map from contextily pkg
nyc = gpd.read_file(gpd.datasets.get_path('nybb')).to_crs(epsg=4326)
ax = nyc.plot(figsize=(10,10), alpha=0.3, edgecolor="k")

# Contextily adding background image
def add_basemap(ax, zoom, url='http://tile.stamen.com/terrain/tileZ/tileX/tileY.png'):
    xmin, xmax, ymin, ymax = ax.axis()
    basemap, extent = ctx.bounds2img(xmin, ymin, xmax, ymax, zoom=zoom, url=url)
    ax.imshow(basemap, extent=extent, interpolation='bilinear')
    # restore original x/y limits
    ax.axis((xmin, xmax, ymin, ymax))

add_basemap(ax, zoom=10, url='http://tile.stamen.com/toner-lite/tileZ/tileX/tileY.png')
# nyc.to_crs(epsg=4326)
gdf.plot(ax=ax, color="red")
ax.set_axis_off()
...