如何使用 geopandas 通过使用底图在更详细的地图上绘制纬度和经度？

Question

我正在尝试在德里地图上绘制一些纬度和经度，我可以通过使用 geopandas 在 python3.8 中使用形状文件来完成这是形状文件的 link：

https://drive.google.com/file/d/1CEScjlcsKFCgdlME21buexHxjCbkb3WE/view?usp=sharing

以下是我在地图上绘制点的代码：

lo=[list of longitudes]
la=[list of latitudes]

delhi_map = gpd.read_file(r'C:\Users\Desktop\Delhi_Wards.shp')
fig,ax = plt.subplots(figsize = (15,15))
delhi_map.plot(ax = ax)
geometry = [Point(xy) for xy in zip(lo,la)]
geo_df = gpd.GeoDataFrame(geometry = geometry)
print(geo_df)
g = geo_df.plot(ax = ax, markersize = 20, color = 'red',marker = '*',label = 'Delhi')
plt.show()

结果如下：

现在这个地图不是很清楚，任何人都无法识别标记的地方所以我尝试通过以下代码使用底图获得更详细的地图：

df = gpd.read_file(r'C:\Users\Jojo\Desktop\Delhi_Wards.shp')
new_df = df.to_crs(epsg=3857)
print(df.crs)
print(new_df.crs)
ax = new_df.plot()
ctx.add_basemap(ax)
plt.show()

结果如下：

我正在获取底图，但我的 shapefile 与它重叠。我可以得到一张地图来绘制我的纬度和经度吗？地图上有更详细的地名或道路名称或类似 google 地图中的任何东西，甚至像被蓝色重叠的地图shapefile 地图？

这样的地图能画出来吗？？

https://www.researchgate.net/profile/P_Jops/publication/324715366/figure/fig3/AS:618748771835906@1524532611545/Map-of-Delhi-reproduced-from-Google-Maps-12.png

Answer 1

我不知道 geopandas。我建议的想法仅使用基本的 python 和 matplotlib。希望您能根据自己的需要进行调整。

背景是下面这张地图。我使用 google-maps.

计算出它角的 GPS 坐标

代码遵循我所说的三点。请注意，使用 imread 和 imshow 会反转 y 坐标。这就是函数 coordinatesOnFigur 在 x 和 y.

中看起来 non-symmetrical 的原因

运行代码在 Montijo 附近生成带有红色子弹的地图（最后有一个小测试）。

import numpy as np
import matplotlib as mpl
import matplotlib.pyplot as plt
from matplotlib import patches
from matplotlib.widgets import Button


NE = (-8.9551, 38.8799)
SE = (-8.9551, 38.6149)
SW = (-9.4068, 38.6149)
NW = (-9.4068, 38.8799)

fig = plt.figure(figsize=(8, 6))
axes = fig.add_subplot(1,1,1, aspect='equal')
img_array = plt.imread("lisbon_2.jpg")
axes.imshow(img_array)

xmax = axes.get_xlim()[1]
ymin = axes.get_ylim()[0]  # the y coordinates are reversed, ymax=0

# print(axes.get_xlim(), xmax)
# print(axes.get_ylim(), ymin)


def coordinatesOnFigure(long, lat, SW=SW, NE=NE, xmax=xmax, ymin=ymin):
    px = xmax/(NE[0]-SW[0])
    qx = -SW[0]*xmax/(NE[0]-SW[0])
    py = -ymin/(NE[1]-SW[1])
    qy = NE[1]*ymin/(NE[1]-SW[1])
    return px*long + qx, py*lat + qy


# plotting a red bullet that corresponds to a GPS location on the map
x, y = coordinatesOnFigure(-9, 38.7)
print("test: on -9, 38.7 we get", x, y)
axes.scatter(x, y, s=40, c='red', alpha=0.9)

plt.show()

Answer 2

使用zorder参数调整层的顺序（lower zorder表示下层），alpha调整多边形。无论如何，我猜，你正在绘制 df 两次，这就是它重叠的原因。

这是我的脚本和结果

import geopandas as gpd
import matplotlib.pyplot as plt
import contextily as ctx
from shapely.geometry import Point

long =[77.2885437011719, 77.231931, 77.198767, 77.2750396728516]
lat = [28.6877899169922, 28.663863, 28.648287, 28.5429172515869]
geometry = [Point(xy) for xy in zip(long,lat)]


wardlink = "New Folder/wards delimited.shp"

ward = gpd.read_file(wardlink, bbox=None, mask=None, rows=None)
geo_df = gpd.GeoDataFrame(geometry = geometry)

ward.crs = {'init':"epsg:4326"}
geo_df.crs = {'init':"epsg:4326"}

# plot the polygon
ax = ward.plot(alpha=0.35, color='#d66058', zorder=1)
# plot the boundary only (without fill), just uncomment
#ax = gpd.GeoSeries(ward.to_crs(epsg=3857)['geometry'].unary_union).boundary.plot(ax=ax, alpha=0.5, color="#ed2518",zorder=2)
ax = gpd.GeoSeries(ward['geometry'].unary_union).boundary.plot(ax=ax, alpha=0.5, color="#ed2518",zorder=2)

# plot the marker
ax = geo_df.plot(ax = ax, markersize = 20, color = 'red',marker = '*',label = 'Delhi', zorder=3)

ctx.add_basemap(ax, crs=geo_df.crs.to_string(), source=ctx.providers.OpenStreetMap.Mapnik)
plt.show()

我不知道 google 地图在上下文中，我认为它不可用。或者，您可以使用显示完全相同地名的 OpenStreetMap 底图，或您可以探索的任何其他底图。在参数中使用 `source` 关键字，例如 `ctx.add_basemap(ax, source=ctx.providers.OpenStreetMap.Mapnik)` 。以下是检查可用提供商和每个提供商提供的地图的方法：

>>> ctx.providers.keys()
dict_keys(['OpenStreetMap', 'OpenSeaMap', 'OpenPtMap', 'OpenTopoMap', 'OpenRailwayMap', 'OpenFireMap', 'SafeCast', 'Thunderforest', 'OpenMapSurfer', 'Hydda', 'MapBox', 'Stamen', 'Esri', 'OpenWeatherMap', 'HERE', 'FreeMapSK', 'MtbMap', 'CartoDB', 'HikeBike', 'BasemapAT', 'nlmaps', 'NASAGIBS', 'NLS', 'JusticeMap', 'Wikimedia', 'GeoportailFrance', 'OneMapSG'])
>>> ctx.providers.OpenStreetMap.keys()
dict_keys(['Mapnik', 'DE', 'CH', 'France', 'HOT', 'BZH'])

如何使用 geopandas 通过使用底图在更详细的地图上绘制纬度和经度？

How to use geopandas to plot latitude and longitude on a more detailed map with by using basemaps?

python

matplotlib

matplotlib-basemap

geopandas

contextily