Картографическая тепловая карта на фоне OpenStreetMap

Question

У меня есть созданный участок карты открытых улиц с использованием картопы:

from __future__ import division
import numpy as np 
import matplotlib as mpl        
import matplotlib.pyplot as plt 
import cartopy.crs as ccrs
import cartopy.io.img_tiles as cimgt

request = cimgt.OSM()
extent = [-89, -88, 41, 42]

ax = plt.axes(projection=request.crs)
ax.set_extent(extent)

ax.add_image(request, 8)
plt.show()

Теперь у меня также есть список точек долготы и широты,Как я могу наложить тепловую карту этих точек долготы и широты на карту улиц?

Я пытался использовать hist2d, но это не работает.

lons = (-88 --89)*np.random.random(100)+-89
lats = (41 - 42)*np.random.random(100)+42
ax.hist2d(lons,lats)
plt.show()

Но это не такработа.

Полагаю, мне нужно где-нибудь добавить аргумент преобразования в команду построения?Но я не уверен, как это сделать.

Спасибо!

Answer 1

Вы правы насчет необходимости преобразования координат.Вот рабочий код с результирующим графиком.

import numpy as np 
import matplotlib as mpl        
import matplotlib.pyplot as plt 
import cartopy.crs as ccrs
import cartopy.io.img_tiles as cimgt

request = cimgt.OSM()
fig, ax = plt.subplots(figsize=(10,16),
                       subplot_kw=dict(projection=request.crs))
extent = [-89, -88, 41, 42]  # (xmin, xmax, ymin, ymax)
ax.set_extent(extent)
ax.add_image(request, 8)

# generate (x, y) centering at (extent[0], extent[2])
x = extent[0] + np.random.randn(1000)
y = extent[2] + np.random.randn(1000)

# do coordinate conversion of (x,y)
xynps = ax.projection.transform_points(ccrs.Geodetic(), x, y)

# make a 2D histogram
h = ax.hist2d(xynps[:,0], xynps[:,1], bins=40, zorder=10, alpha=0.5)
#h: (counts, xedges, yedges, image)

cbar = plt.colorbar(h[3], ax=ax, shrink=0.45, format='%.1f')  # h[3]: image

plt.show()

Конечный график:

Редактировать 1

Когда ax создается с plt.subplots(), для него определена определенная проекция.В этом случае проекция определяется ключевым словом projection=request.crs.Чтобы построить что-то на ax, вы должны использовать его систему координат.

Преобразование координат выполняется с помощью функции transform_points() в выражении

xynps=ax.projection.transform_points(ccrs.Geodetic(), x, y)

где

• (x, y) - список значений (долгота, широта),
• ccrs.Geodetic () означает, что значения (долгота, широта) в градусах

Возвращаемые значения xynps - это массив координат карты.Он имеет 2 столбца для x и y, которые находятся в соответствующей системе координат, используемой текущим ax.