Plotting a grid to a map with a dynamic plot using qmesh. Xarray dataset is lats/lons and heights. Not sure if there is an easy fix to have a bigger sampling size for each raster cube.
shoal_grid = data_xr.hvplot.quadmesh(x='X_lon', y='Y_lat', rasterize=True,
geo=True,dynamic=True)
geo_bg = hv.element.tiles.EsriImagery().opts(width=600, bgcolor='black')
tap = hv.streams.Tap(source=shoal_grid,x=-93.34010, y=29.92420)
def drilldown(x,y):
test_point = np.argmin((x-lons)**2+(y-lats)**2)
temp_ts = elvs.sel(points=test_point)
temp_ts['time'] = time_series
return temp_ts.hvplot.line(x='time',y='elevations',width=500,dynamic=False)
tap_dmap = hv.DynamicMap(drilldown, streams=[tap])
geo_bg*shoal_grid.opts(height=500, tools=['hover','box_select', 'lasso_select'],active_tools=
['box_select'], **geo_opts) + tap_dmap
Zoomed in: