  We present an algorithm for computing $\epsilon$-coresets for $(k,
\ell)$-median clustering of polygonal curves in $\mathbb{R}^d$ under the
Fr\'echet distance. This type of clustering is an adaption of Euclidean
$k$-median clustering: we are given a set of $n$ polygonal curves in
$\mathbb{R}^d$, each of complexity (number of vertices) at most $m$, and want
to compute $k$ median curves such that the sum of distances from the given
curves to their closest median curve is minimal. Additionally, we restrict the
complexity of the median curves to be at most $\ell$ each, to suppress
overfitting, a problem specific for sequential data. Our algorithm has running
time linear in $n$, sub-quartic in $m$ and quadratic in $\epsilon^{-1}$. With
high probability it returns $\epsilon$-coresets of size quadratic in
$\epsilon^{-1}$ and logarithmic in $n$ and $m$. We achieve this result by
applying the improved $\epsilon$-coreset framework by Langberg and Feldman to a
generalized $k$-median problem over an arbitrary metric space. Later we combine
this result with the recent result by Driemel et al. on the VC dimension of
metric balls under the Fr\'echet distance. Furthermore, our framework yields
$\epsilon$-coresets for any generalized $k$-median problem where the range
space induced by the open metric balls of the underlying space has bounded VC
dimension, which is of independent interest. Finally, we show that our
$\epsilon$-coresets can be used to improve the running time of an existing
approximation algorithm for $(1,\ell)$-median clustering.
