使用数据着色器、全息视图和散景进行地理映射时遇到问题

Question

我正在尝试使用全息视图、数据着色器和散景将谷歌手机历史位置映射到地图上。大多数情况下与 datashader 网站中给出的示例非常相似。但是当我做地图覆盖时，因为纬度/经度被弄乱了，所以地图覆盖不起作用。

import datashader as ds
import geoviews as gv
import holoviews as hv
from holoviews.operation.datashader import datashade, dynspread
from datashader import transfer_functions as tf
from colorcet import fire
hv.extension('bokeh')

> df2.head()

lat long
0   -37.7997515 144.9636466
1   -37.7997515 144.9636466
2   -37.7997369 144.9636036
3   -37.7997387 144.9636358
4   -37.7997515 144.9636466

这可以生成数据的图像，

ds_viz = ds.Canvas().points(df2,'lat','long')  
tf.set_background(tf.shade(ds_viz, cmap=fire),"black")

但是，当我尝试用地图覆盖它时，它不起作用，

从 bokeh.models 导入 WMTSTileSource

url = 'https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{Z}/{Y}/{X}.jpg'
tile_opts  = dict(width=1000,height=600,bgcolor='black',show_grid=False)
map_tiles  = gv.WMTS(url).opts(style=dict(alpha=0.5), plot=tile_opts)
points     = hv.Points(df2, kdims=['long','lat'])
trips = datashade(points, cmap=fire,width=1000, height=600)

map_tiles * trips

我究竟做错了什么？

Answer 1

看起来您的点在 lon,lat 中，但您的地图在 Web Mercator 坐标中，因此您需要在叠加之前将您的点投影到 Web Mercator 中。GeoViews 提供了对投影的全面支持，但对于这种特定情况，Datashader 提供了专用功能datashader.utils.lnglat_to_meters。这样的事情应该工作：

df2.loc[:, 'lon'], df.loc[:, 'lat'] = lnglat_to_meters(df2.lon,df2.lat)

投影可能很慢，因此您可能希望将生成的 df2 保存到 Parquet 文件中，这样您只需执行一次。