from tomogram_datasets import all_fm_tomograms
import numpy as np
from visualize_voxels import visualize as viz

tomos = all_fm_tomograms()
n_flagellar_motors = [len(tomo.annotation_points()) for tomo in tomos]
super_tomo = tomos[np.argsort(n_flagellar_motors)[-3]]

viz(
    super_tomo.get_data(),
    marks=super_tomo.annotation_points(),
    markalpha=0.5,
    axis=0,
    slices=np.linspace(80, 320, 100),
    fps=16
)

Rendering visualization:   0%|          | 0/100 [00:00<?, ?frame/s]