import datashader as ds
import holoviews as hv
import numpy as np
import pandas as pd
from holoviews.operation.datashader import rasterize

hv.extension("bokeh")

num = 10000
np.random.seed(1)

dists = {
    cat: pd.DataFrame({
        "x": np.random.normal(x, s, num),
        "y": np.random.normal(y, s, num),
        "s": s,
        "val": val,
        "cat": cat,
    })
    for x, y, s, val, cat in [
        (2, 2, 0.03, 0, "d1"),
        (2, -2, 0.10, 1, "d2"),
        (-2, -2, 0.50, 2, "d3"),
        (-2, 2, 1.00, 3, "d4"),
        (0, 0, 3.00, 4, "d5"),
    ]
}

df = pd.concat(dists, ignore_index=True)

points = hv.Points(df)
plot = rasterize(points, aggregator=ds.count_cat("cat"))

plot + plot.select(cat=["d1", "d3"])