Paper deep dive
Scaling Time Series Classification via XAI-Driven Data Reduction
Davide Italo Serramazza, Thach Le Nguyen, Georgiana Ifrim
Intelligence
Status: succeeded | Model: Gemma-4-26B-A4B | Prompt: intel-v1 | Confidence: 92%
Last extracted: 7/23/2026, 2:03:47 AM
Summary
The paper introduces drXAI, a methodology that leverages Explainable AI (XAI) attribution methods to perform data reduction for Time Series Classification (TSC). By using a fast GPU-accelerated classifier (Hydra) to generate local attributions and aggregating them into global feature importance scores, drXAI selects salient features via an automated elbow-cut heuristic. This approach achieves 80-90% data reduction while maintaining classification accuracy, enabling resource-intensive models like ConvTran to scale to large datasets previously inaccessible due to memory constraints.
Entities (8)
Relation Signals (7)
Hydra → hascomplexity → linear
confidence 95% · Since k is usually a fixed constant, the time complexity is effectively linear O(ndL).
drXAI → uses → Hydra
confidence 95% · drXAI addresses this by using a fast, GPU-accelerated classifier (Hydra) to generate local attributions.
drXAI → enables → ConvTran
confidence 92% · we show that drXAI allows resource-intensive models like ConvTran to scale to datasets that were previously inaccessible due to memory constraints.
drXAI → employs → Feature Ablation
confidence 90% · apply lightweight explainers such as Feature Ablation.
Transformers → hascomplexity → quadratic
confidence 90% · state-of-the-art models, such as Transformers, exhibit quadratic complexity relative to sequence length
drXAI → employs → Shapley Value Sampling
confidence 88% · We use two permutation-based attribution methods... Shapley Value Sampling (SVS)
drXAI → utilizes → Proto
confidence 85% · propose Proto, a single-instance background outperforming the standard zeros baseline
Cypher Suggestions (0)
No Cypher suggestions yet.
Abstract
Abstract:Explainable AI (XAI) for time series has seen significant algorithmic growth, but its utility in providing measurable performance gains for downstream tasks remains under-explored. This paper bridges this gap by introducing drXAI, a novel methodology that repurposes XAI attribution methods for effective data reduction in Time Series Classification (TSC). The core challenge in modern TSC is scalability; state-of-the-art models, such as Transformers, exhibit quadratic complexity relative to sequence length and linear complexity relative to the number of channels. This renders them computationally prohibitive for massive datasets. drXAI addresses this by using a fast, GPU-accelerated classifier (Hydra) to generate local attributions. We aggregate these into global feature importance scores and employ an automated elbow-cut heuristic to select the most salient features without requiring manual thresholds. We evaluate our approach on both synthetic and real-world univariate and multivariate datasets. On synthetic benchmarks, drXAI successfully recovers ground-truth features where traditional baselines fail. On real-world data, drXAI achieves between 80% and 90% data reduction while maintaining classification accuracy comparable to models trained on the full dataset. Most importantly, we show that drXAI allows resource-intensive models like ConvTran to scale to datasets that were previously inaccessible due to memory constraints. Our results show the benefits of using XAI not just for interpretability, but as a robust tool for feature selection and scalability in time series analysis. All our code and data are openly available.
Tags
Links
- Source: https://arxiv.org/abs/2607.15774v2
- Canonical: https://arxiv.org/abs/2607.15774v2
Trouble viewing inline? Open PDF directly →
Full Text
47,927 characters extracted from source content.
Expand or collapse full text
Scaling Time Series Classification via XAI-Driven Data Reduction Davide Italo Serramazza, Thach Le Nguyen, and Georgiana Ifrim School of Computer Science, University College Dublin, Ireland davide.serramazza@ucdconnect.ie thach.lenguyen,georgiana.ifrim@ucd.ie Abstract. Explainable AI (XAI) for time series has seen significant al- gorithmic growth, but its utility in providing measurable performance gains for downstream tasks remains under-explored. This paper bridges this gap by introducing drXAI, a novel methodology that repurposes XAI attribution methods for effective data reduction in Time Series Classi- fication (TSC). The core challenge in modern TSC is scalability; state- of-the-art models, such as Transformers, exhibit quadratic complexity relative to sequence length and linear complexity relative to the number of channels. This renders them computationally prohibitive for massive datasets. drXAI addresses this by using a fast, GPU-accelerated classifier (Hydra) to generate local attributions. We aggregate these into global feature importance scores and employ an automated elbow-cut heuristic to select the most salient features without requiring manual thresholds. We evaluate our approach on both synthetic and real-world univariate and multivariate datasets. On synthetic benchmarks, drXAI successfully recovers ground-truth features where traditional baselines fail. On real- world data, drXAI achieves 80–90% data reduction while maintaining classification accuracy comparable to models trained on the full dataset. Most importantly, we show that drXAI allows resource-intensive mod- els like ConvTran to scale to datasets that were previously inaccessible due to memory constraints. Our results show the benefits of using XAI not just for interpretability, but as a robust tool for feature selection and scalability in time series analysis. All our code and data are openly available. 1 Introduction The field of Explainable AI (XAI) has experienced significant growth in recent years, particularly within the Time Series Classification (TSC) domain. A major focus of XAI is feature attribution, which quantifies the importance of input features for the model predictions. This area has seen considerable advances in the efficiency and effectiveness of attribution methods for time series data [24]. Despite these algorithmic advancements, the utility of XAI as a tool for achieving measurable performance gains in other computational tasks remains underexplored. Our work offers a paradigm shift: viewing XAI not only as in- terpreting complex models but also as a practical tool for enhancing model effi- ciency and scalability. This paper bridges this gap by introducing drXAI (Data arXiv:2607.15774v2 [cs.LG] 22 Jul 2026 2Serramazza et al. Reduction with XAI), a novel attribution-agnostic methodology that repur- poses attribution methods for effective feature selection in TSC. For many State- of-the-Art (SOTA) TSC models, the time complexity grows linearly with the number of features or even quadratically with sequence length (e.g., transformers); additionally, these models also require a huge amount of memory. This computational burden is often a limitation for training on massive, high- dimensional datasets. A way to scale these models is to train them on a reduced, yet highly informative, set of features. Our work directly enables this by using XAI to identify the most important features, thus reducing data dimensionality without sacrificing critical information. To our knowledge, this work is the first to use XAI for feature selection to scale time series classification methods. drXAI is a wrapper feature selection method that operates in two stages: it first trains a TSC model, referred to as the explainer classifier, and then gener- ates explanations for its predictions. In this work, accounting for speed, we use the GPU implementation of Hydra [4] as the explainer classifier, restrict attri- bution to the explainer set, a subset of the training set, and apply lightweight explainers such as Feature Ablation. For Multivariate Time Series Classification (MTSC) datasets, we focus on channel selection. We propose a fast approach to compute global channel im- portance by aggregating the local attribution values, ranking the channels, and selecting a subset. We then retrain SOTA classifiers on the reduced data and measure the accuracy and computational gains. For UTSC datasets, we apply the same algorithm for computing global time point importance, which is then used for data reduction via consecutive time point selection. This is especially effective for very long time series, where SOTA TSC methods do not scale well due to the requirement of extensive computational resources. Our main contributions in this paper are: 1. We develop drXAI, a general framework that leverages XAI attribution to select features from MTSC and UTSC datasets. The algorithm is attribution- agnostic and is general enough to support both channel and time-point se- lection. As part of this, we also investigate the background data used for simulating data missingness in attribution methods and propose Proto, a single-instance background outperforming the standard zeros baseline com- mon in popular attribution methods [10]. 2. Using synthetic MTSC data, we show that drXAI selects only informative channels, unlike baseline methods. On synthetic UTSC data, it successfully selects over 90% of relevant features, whereas baselines are limited to at most 33%. On real-world datasets, drXAI provides the best trade-off between data reduction and accuracy, most of the time matching the models trained on all features, for both MTSC and UTSC experiments. 3. We demonstrate that by using a fast classifier (Hydra) and an attribution method (Feature Ablation) for data selection, we can effectively train more accurate but resource-intensive models (e.g., ConvTran, MultiRocket-Hydra [13]). This enables computationally expensive models to train successfully, overcoming the bottleneck in use cases where the model cannot be trained on the entire dataset due to memory limits. XAI-Driven Data Reduction3 4. To encourage further research and reproducibility in this area, we make all our data and code publicly available 1 . 2 Background 2.1 Problem Definition We represent a time series dataset D as a tensor, of dimensions n× d× L where n, d, and L respectively represent the number of samples, channels, and time points in the data. In the UTS case, d = 1. A (trained) classification model clf predicts the class c i of a time series instance D i (1≤ i≤ n). An explainer exp explains the clf predictions on D by producing a set of attribution maps A, a tensor of attributions with the same dimensions as D (n×d×L) whereA i,j,k indicates the attribution (i.e., relevance) of D i,j,k for the prediction of the model clf on the instance D i . We denote the set of features as F = f 1 , f 2 , . . . , f m where m is the total number of features. Each feature represents a portion of the time-series data, e.g., data within the same time segment, the same time step, or the same channel. The union of all features is the complete time series. We focus on the task of selecting a subset of features F sel ⊂ F. We informally refer to this new set as selected features, denoting as D F sel the new dataset containing only these features. In this work, we focus on data reduction for TSC, organised into two subtasks: the channel selection problem for MTSC datasets and time points selection for UTSC datasets. For the former, each feature represents a channel (m = d) while for the latter, each feature represents a time point (m = L). To avoid confusion, we specify which type of features we focus on when describing specific data reduction tasks. 2.2 Time Series Classification Recent work in TSC has significantly advanced the availability of extensive TSC benchmarks [2,5] as well as the accuracy and efficiency of TSC algorithms [7,4]. While many algorithms achieve top accuracy on benchmarks, they are resource- intensive, especially for large-scale datasets or very long time series [13]. We discuss a few relevant SOTA classifiers, as well as their computational complexity. Explainer Classifier. Our proposed method, drXAI, is a wrapper method for feature selection and requires a fast classifier to be explained. Hydra is a TS transformation algorithm combining convolution-based and dictionary-based aspects: g groups of k convolutional kernels slide over the TS, and for each group, the closest-matching kernels are counted at each time-point. Features are then fed to a Ridge Classifier. The time complexity of Hydra is dominated by the convolution and the competitive counting process. Since k is usually a fixed constant, the time complexity is effectively linear O(ndL). Hydra’s memory 1 https://github.com/mlgig/drXAI 4Serramazza et al. complexity is O(nk + kd). In our experiments, we used the fast Hydra GPU implementation from [4], which scales well to very large datasets. SOTA Time Series Classifiers. To assess the quality of data selection F sel , we trained the following 3 SOTA classifiers on the reduced datasets D F sel . MultiRocket-Hydra (MRH) [13] uses the Hydra features concatenated with those from MultiRocket [23]. MultiRocket applies kernels also to the first- order difference of the series and extracts features via 4 pooling operators (MPV, MIPV, LSPV, PPV). The MultiRocket pipeline is heavier computationally than Hydra. While still linear in the inputs, the constants for MultiRocket are larger, as well as the number of features extracted (by default around 50,000 features), which makes it more expensive both time and memory wise, especially for long time series. Moreover, since by default each kernel is applied to at most 8 random channels, both algorithms also benefit from selections shorter than 8 channels. Beyond efficiency, restricting the input to informative features can further improve model robustness and accuracy. ConvTran [7] is a recent transformer model, tailored for MTSC. It extracts features from the raw series using convolutional layers that are then fed into the transformer after tokenization. ConvTran comes with a significant jump in complexity as compared to Hydra. While Hydra uses random kernels and a linear classifier, ConvTran is a fully trainable Transformer-CNN hybrid that utilizes self-attention. ConvTran’s complexity is driven by the quadratic nature of the attention mechanism and the parameters of its convolutional embedding layers. Its cost is quadratic w.r.t the series length O(nL 2 d embedding ) and linear w.r.t the number of channels. For very long sequences where L > 1000, this becomes significantly slower than Hydra’s O(nL). InceptionTime [9] is another deep learning classifier, specifically, an ensem- ble method composed of 5 vanilla Inception networks (CNN). Its complexity is primarily determined by its deep Inception modules, which apply parallel con- volutions of varying kernel sizes through bottleneck layers, capturing both local and long-range temporal relationships. Unlike ConvTran, InceptionTime is linear with respect to sequence length and channels, but since it requires backpropa- gation, the constant factors and hardware requirements are much higher than those of Hydra. 2.3 XAI Methods for Time Series Classification Our methodology is attribution-method agnostic, i.e., it requires only an attri- bution map as input, regardless of which algorithm computed it. In this work, we select two permutation-based attribution methods that were adapted for TS and work efficiently in this domain [19,25]. These methods require a back- ground set to simulate data missingness when computing the attribution. No gradient-based attribution methods were considered for this study primarily because they can’t work with non-gradient-based algorithms, like Hydra and MultiRocket. For the rest of this paper, we use explainers to refer to the XAI attribution methods. XAI-Driven Data Reduction5 Shapley Value Sampling (SVS) is an approximation of the Shapley values as described in the original work [11]. It applies sampling to the SHAP formula, randomly permuting the features to be explained, adding them sequentially to each sample in the background set. The attribution values are the change in the model output resulting from these substitutions. Although shown to be effec- tive with regard to pointing out important features [21], this method requires a vast amount of computation time, due to the high number of feature permuta- tions. Some adaptations for time series, such as grouping features through TS segmentation, have been shown to drastically reduce computation time, while preserving accuracy [19]. Feature Ablation (FA) [10] is much simpler and faster when compared to SHAP. It sequentially replaces each feature with the corresponding values of the samples found in the background dataset. As for SVS, the attribution values are the differences in the model output after substituting feature values. We used the implementations provided in [20] for both explainers. 2.4 Feature Selection for TSC We first discuss recent work on channel selection for MTSC. In [6], the au- thors propose a supervised algorithm for selecting a subset of channels for MTSC datasets. The algorithm has two variants, ECS and ECP, where the channels are selected based on their discriminative power, estimated using the Euclidean dis- tance between class centroids: a higher distance implies higher discriminative power. According to the original paper, this filter method can reduce, on aver- age, 70% of the data without compromising the accuracy of classifiers. Another recent work that directly compares to ECS and ECP is TSelect [15], a wrapper method which trains for each channel a Logistic Regression classifier based on 5 computationally cheap features extracted channel-wise. These mod- els are then used to determine which of the relative channels to retrain: various filters based on these model accuracies, redundant predictions, etc., select which channels to discard. Reduction of time series data can also be performed by projecting the high- dimensional time series to a lower-dimensional latent space [1]. Known methods include Autoencoders, Principal Component Analysis, Singular Value Decompo- sition, Discrete Fourier/Wavelet Transform, and Down-sampling. However, these are dimensionality reduction techniques that project time series into a dif- ferent space, losing the original features. For many applications, it is important to keep the data in the original representation to be able to audit important features, for example, medical applications such as monitoring human health. Therefore, these transformation methods are outside the scope of this paper. For time point selection techniques for UTSC, a simple and effective approach, especially for long time series, is to use Random Forest importance (RFI), using the importance (impurity reduction) of each feature computed dur- ing training. Another alternative is the Mutual Information (MI) between each feature and targets for a TSC task [17]. 6Serramazza et al. 3 Proposed Methodology In this section, we describe drXAI, our algorithm for time series data reduction using XAI attribution methods to identify and select important features from TSC datasets. Algorithm 1 drXAI Algorithm to Select a Feature Subset from a Time Series Dataset Require: clf, exp, D, n samp ▷ explainer classifier, explainer, train set, n. samples per class 1: Train clf on D 2: D exp = sample n samp per class 3: A = Explain clf classifications for samples in D exp using exp 4: A = equation (1) if channel selection else equation (2) ▷ details in Section 3.1 5: fr avg = aggr_avg(A) ; fr abs = aggr_abs(A)▷ details in Section 3.2 6: F sel_avg = elbow_cut(fr avg ) ; F sel_abs = elbow_cut(fr abs ) ▷ details in Section 3.3 7: F sel = F sel_abs ∩ F sel_avg 8: return F sel Our methodology (Algorithm 1) relies on the explanation set D exp , which is a subset of the training set (assuring there is no information leakage), composed by randomly sampling n samp per class, the explanation classifier clf i.e., the model which is explained (in this work Hydra) using the explainer exp. The key functions are: – The aggr_avg and aggr_abs functions, which aggregate the attribution val- ues of different samples based on two complementary strategies. – Elbow cut which selects the top K features F sel . 3.1 Computing Feature Attributions After training clf using the training set D and instantiating the explanation set D exp , the result of applying the explainer exp to the classifier clf and the explanation set D exp are the local attributions A ∈ R n×d×L where n = |D exp |. The next step is to calculate the feature attribution of the feature set F for each time series in D exp . The attribution of a feature is simply the average attribution of all data points represented by the feature. In particular, if features represent channels (for channel selection), the attribution of a feature is: A i,j = 1 L L X t=1 A i,j,t (1) where A i,j is the attribution of channel j (or feature F j ) in time series D i . This is the row aggregation step in Figure 1. XAI-Driven Data Reduction7 Row Aggregation Time series data D Local Attribution D 1 D 2 D n Explainer Model Local Channel Attribution A ... aggr_abs aggr_avg Channel Scores Selection with Elbow Cut Selected Channels Intersection Selected Channels Fig. 1: drXAI: Channel selection for MTSC using XAI scores computed from time series attributions. Similarly, if each feature represents a time point (for time point selection), the feature attribution is: A i,t = 1 d d X j=1 A i,j,t (2) where A i,t is the attribution of time point t (or feature F t ) in time series D i . Stacking the attributions for each feature and each sample, we obtain the matrix A. To be noted that the next steps are executed regardless of the selection type (channel or time points). Attribution Background Set. As mentioned in Section 2.3, both explainers require a background dataset. Since the computational complexity of explaining linearly increases with the cardinality of this set, we study two different single- sample backgrounds b. Zeros background (Zeros). A default choice for most explanation frameworks, i.e., a time series full of zeros. Although this is conceptually very simple, it is an unrealistic sample, potentially leading to unreasonable explanations. Mathe- matically, this is defined as: b = 0 d×L Class prototypes average (Proto). The background we propose uses the prototype of each class c: p c = 1 |D c | X D i ∈D c D i (3) where D c ⊂ D is the set of all samples in the training set, belonging to class c. Let C be the set of all classes. The proto background is defined as the average of class prototypes: b = 1 |C| X c∈C p c (4) 8Serramazza et al. 3.2 Aggregating Feature Attribution over Samples The previous step computes the sample-wise attributions for each feature in the feature set F. The next step is to aggregate these local attributions over the explanation set D exp to obtain global attributions. This results in the feature importance scores fr avg ,fr abs ∈ R m + , which can be used for feature selection (e.g., this is the channel scores vector in Figure 1). We use two complementary approaches for this aggregation which proved effective in our experiments. aggr_avg: In this scenario, the feature importance score is averaged over all samples first, then the absolute value is computed. This strategy aims to de- emphasize uncertain features, i.e., features that have mixed negative and positive attribution signs (thus roles) across the dataset. f r avg =| 1 n n X i=1 A i |(5) where A i is the i-th a row of A, i.e., attribution vector of the sample D i . aggr_abs: On the other hand, this strategy averages the absolute value of the attributions for each feature: f r abs = 1 n n X i=1 |A i |(6) This strategy simply detects the most active features, regardless of the sign. The dimension of these vectors is m = d for channel selection and m = L for time point selection. 3.3 Feature Selection using Elbow Cut The elbow cut of a curve is a common heuristic to choose a point where intuitively the diminishing returns are no longer worth the additional cost; e.g., it is often used to select the number of clusters while running k-means algorithms, and in the channel selection context for the ECP and ECS methods as described in [6]. In drXAI, the elbow cut is used after all features are sorted by their computed score, to automatically select the number of top features to retain: F sel_abs and F sel_avg are respectively the results of the elbow cut applications on f r abs and f r avg . The set of final selected features F sel is the intersection of f r abs and f r avg (Line 7 Algorithm 1). Intuitively, this set contains only features that are both magnitude- and sign-wise important, thereby including only the essential ones. 3.4 Time and Space Complexity for drXAI The computational efficiency of the drXAI framework is a primary contribu- tion, specifically designed to mitigate the prohibitive costs of training classifiers on high-dimensional time series. The total complexity of the pipeline is the sum of three distinct phases: (i) Explainer Classifier Training, (i) Attribution Gen- eration, and (i) Feature Selection. XAI-Driven Data Reduction9 Explainer Classifier Training. We use Hydra as the core classifier for its linear-scaling properties. For a dataset D(n, d, L), Hydra extracts features using k random convolutional kernels Since k is a fixed hyperparameter, this phase remains O(n· d· L), which is asymptotically optimal for time series processing. Attribution Generation. The complexity of this phase is the product of the number of samples to be explained and the cost of the chosen explainer: – Feature Ablation (FA): For each sample, FA requires a forward pass for each feature. For channel selection, this is O(d·Cost Hydra ), and for time-point selection, O(L· Cost Hydra ). – SHAP: SHAP estimates values via sampling. While the number of samples s scales with the feature space, Hydra’s GPU-accelerated inference, coupled with TS segmentation, allows SHAP to remain tractable even for L > 1, 000. Feature Selection. This is a step that has a constant time involving sorting the features by importance and applying the elbow cut. Space Complexity. The memory footprint is dominated by the storage of attribution maps O(D exp · d· L). However, the subsequent elbow-cut selection enables a significant reduction in the memory required for training deep models. As demonstrated in our results (Section 4), this reduction is the critical factor enabling ConvTran to run on datasets where it would otherwise trigger Out- of-Memory (OOM) errors. 4 Experiments In our experiments, we consider 4 configurations of our methods: the combina- tions of previously listed explainers and backgrounds, SHAP Proto, FA Proto, SHAP zeros, and FA zeros. We hypothesize that the informative Proto back- ground gives a small boost to our method compared to the uninformative zeros. We also hypothesize that, because SHAP is a more complex algorithm, it yields better selection than FA. To assess the quality of the selected features, the following pipeline is applied to each SOTA classifier described in Section 2.2, and each dataset described in Sections 4.2 and 4.3. – Hydra is trained on the current dataset. This allows the application of the drXAI algorithm in the 4 configurations previously listed (e.g., drXAI- SHAP-Proto). Specifically, for each of those, we get the selected features F sel . – Each baseline, along with a random selection method, is evaluated on the dataset to obtain its F sel . – For each of the F sel , we trained each SOTA classifier 3 times using the reduced dataset D F sel , recording the mean accuracy and the mean time for training plus inference. Multiple training rounds (3) were done to account for the stability of the trained classifiers. – Our evaluation of each selection is based on mean accuracy over the 3 runs, and on the percentage of saved data. 10Serramazza et al. drXAI-FA-ProtodrXAI-SHAP-ProtodrXAI-FA-zerosdrXAI-SHAP-zerosECPECSTSelectrandom n. informative 201994212442 n. uninformative 20000017265 mean accuracy.735.658.586.555.634.578.5550.536 Table 1: Number of informative/uninformative channels selected for synthetic MTSC data and mean accuracy (across 3 runs of SOTA classifiers). Experiments were conducted using a machine with AMD EPYC 9654P CPU (96 cores, 192 threads), NVIDIA GeForce RTX 4090 GPU (24GB VRAM), and 1.5Tb of RAM. Considering the breadth of experiments, the main article reports only the most salient results, while the Appendix provides more details. 4.1 Classifier Training For MRH, we implemented our version based on MultiRocket and Hydra trans- formations from the aeon library [12] and the RidgeCV Classifier in sklearn [17]. For Hydra, we used the default hyperparameters, setting the batch size to 256. For ConvTran and InceptionTime training, we used the strategy of [5], i.e., re- serving 10% of the training set as a validation set. We allow up to 100 epochs with early stopping, using the validation loss as criterion. We set the batch size to 256. Finally, for computationally demanding UTSC datasets, we make two adjust- ments: we reduce the batch size for InceptionTime and ConvTran to fit the GPU memory and switch to an iterative solver for the Ridge Classifier in MultiRocket- Hydra, which is required when its input matrix exceeds a maximum size. 4.2 Channel Selection for MTSC Datasets The MTSC datasets have fewer samples than the UTSC ones (Appendix). Thus, we set n samp , the number of sampled instances per class composing D exp , to 50. We also set the maximum non-improving epochs before early stopping for ConvTran and InceptionTime to 20. drXAI is compared against three recent baselines, ECP, ECS, and TSelect, as well as a random selection baseline that samples both the number of channels to retain and the channels to keep from a uniform distribution. Synthetic Dataset. For the MTSC synthetic dataset, we used the data gener- ator code from [19]. Each channel is composed of a lower frequency sine wave. For each sample, two higher-frequency support waves, shorter than the previous ones and with frequencies varying within a specific range, are injected into two randomly selected informative channels. The binary classification task is whether the sum of these two frequencies exceeds a threshold. We generate 5,000 samples for each of the train and validation sets. The series length is set to 1,000 time points; there are 20 informative and 20 uninformative channels (40 total). Using this controlled dataset, we can evaluate how many of XAI-Driven Data Reduction11 Fig. 2: Mean accuracy of each selection (and All Features) for the 4 MTSC real- world datasets and the 3 SOTA classifiers, yielding 12 results in total. the first 20 informative and the last 20 uninformative channels are selected: the ideal selection is all of the informative, none of the uninformative channels. Table 1 shows that only our method exclusively selects informative chan- nels; using the Proto baseline, FA selects 19 out of 20 channels, and SHAP selects 9. This results in these two selections achieving the top two accuracies. Specifically, the drXAI-FA-Proto configuration outperforms the best baseline (ECP) by 10 percentage points in mean accuracy. Real-world Data. For real-world datasets, we focus on those with a large number of channels and time points, where the benefits of feature selection are more pronounced. In order to test drXAI on large-scale MTSC datasets, we selected Face Detection [16] (144 channels, 62 time points), Arc Loss [26] (96 channels, 1101 time points), Military Press and Rowing (mean-centered datasets, 50 channels and 161 time points each) [22]. We summarize the findings using Multi Comparison Matrices (MCM) [8] for mean accuracy and for percentage of data saved, respectively, in Figure 2 and 3. Figure 2 highlights that channel reduction is an important task, as 3 con- figurations of our method and 2 baselines have a higher average accuracy than using all features in the dataset. Our most accurate configuration, drXAI-SHAP- Proto, stands between the best 2 baselines accuracies, i.e., ECP and ECS. Using the faster FA explainer, our proposed background, Proto, is worse than the zero background, mainly due to poor performance when coupled with the MP dataset (see Appendix), although it is better than both ECP and ECS on 3 out of 12 experiments. TSelect has the lowest accuracy among the compared methods. Analyzing the percentage of saved data (i.e., how many channels are discarded), Figure 3 shows that each configuration of our method considerably saves more data than ECP and ECS, with at least 80% reduction. In this regard, the zero-background achieves a better reduction than Proto. Globally, our method is the best trade-off between accuracy and data saved, as different configurations achieve very good performance in both aspects, while ECP and ECS excel only in accuracy and TSelect in data reduction. 12Serramazza et al. Fig. 3: Mean percentage of data saved by each selection for the 4 MTSC datasets. drXAI-FAProtodrXAI-SHAPProtodrXAI-FAzerosdrXAI-SHAPzerosRFIMIrandom n. informative 10k9,00010,00060003000325217676,000 n. uninformative 10k0000002,000 mean accuracy.8500.882.801.608.694.7490.627 Table 2: Number of informative and uninformative time points for synthetic UTSC data and relative mean accuracy of compared feature selection methods (across the 3 runs of each of the 3 SOTA classifiers). 4.3 Time Point Selection on UTSC Datasets Since UTSC datasets have more samples than the MTSC ones (see Appendix), we empirically set n samp to max 100 samples per class, rather than 50. We also decrease the number of non-improving epochs before early stopping to 10. Lastly, since defining each time point as a feature in F sel would be extremely expensive for computing attribution, we instead group them in 20 consecutive, equal-length windows as done in [19]. This means that all time points within a window are either all selected or all discarded. The baselines considered in this section are: the filter method Mutual In- formation (MI) [17], the wrapper method Random Forest feature importance (RFI), and random selection. Similar to the channel selection procedure, ran- dom selection samples both the number and the specific time points from a uniform distribution. Other possible baselines, such as recursive feature elimi- nation, are not considered due to their high computational cost, especially for datasets with long time series. Synthetic Dataset. For the UTSC synthetic data, we used the data generator code included in [14], mirroring the MTSC case, but with d = 1 channels. Two support waves, injected in random places within the first 10,000 informative time points, define the same binary task as in the MTSC case. This informative area is followed by another 10,000 uninformative points. In this case, since using such a long series has a severe consequence on the running times of models, we used only 1,000 samples for training, and kept 5,000 for the test set. As in the former XAI-Driven Data Reduction13 Fig. 4: Mean accuracy of each selection (and All Features) for the 5 UTSC datasets used and the 3 classifiers, yielding 15 results in total. MTSC case, ideally, the first half of the features is selected, while the remaining second is discarded. Table 2 shows that each evaluated method only selects from relevant time points. Among our configurations, drXAI-SHAP-Proto has a perfect se- lection, selecting all 10,000 important features; drXAI-FA-Proto selected 9,000. Among the remaining methods, only drXAI-FA-Zeros, selecting 6,000 features, achieves a good accuracy, since the baselines RFI and MI, as well as SHAP-Zero, select at most 3,252 features, retaining insufficient signal for accurate classifica- tion. ConvTran cannot be trained using all features due to exceeding the GPU VRAM available. Nevertheless, it can run on reduced data of each selection: using drXAI- SHAP-Proto, it achieves a 0.921 average accuracy (Appendix). Real-world Data. For real-world datasets, we used large-scale datasets Cor- nell Whale Challenge, Mosquitos Sound and Whale Sounds from the MONSTER benchmark [5], Right Whale Calls [3] and Urban Sound [18]. These datasets were chosen due to the number of time points (2.5-44k length) and samples (2.7-84k). We note that using the Urban Sound dataset, which has 44k time points, Con- vTran can be trained only using the reduced data after feature selection, as it otherwise exceeds the GPU memory limit. Figure 4 and 5 show, respectively, the MCM for accuracy across all datasets and classifiers and the MCM for percentages of data saved of each selection. In this case, all configurations of our method achieve higher mean accuracy compared to the baselines. The most accurate configuration is drXAI-FA- Proto, having accuracy comparable to the original datasets using all features. Configurations using the informative Proto background have a big margin over those using zeros. Comparing the baselines, MI outperforms RFI. Focusing on data saved, RFI has the largest average save, followed by our method, and lastly MI. Among our configurations, the more expensive and accurate SHAP can save more data than the FA, and for both explainers, the Proto background saves 14Serramazza et al. Fig. 5: Mean percentage of data saved of each selection for the 5 UTSC datasets. more data than the zero background. We note that our configurations, exclud- ing FA Zeros, save at least 80% of the data. As with MTSC, drXAI is the only method achieving a good trade-off between accuracy and data reduction. RFI reduces data aggressively but at the cost of the worst accuracy, while MI shows no clear advantage in either aspect. 4.4 Accuracy-Time Trade-off Analysis An aspect worth analyzing is the total time required by drXAI compared to directly training the model using the original data. In our method, in addition to the cost of training the SOTA classifiers using the reduced data, the time to train Hydra and to compute the explanation (and thus the selection) must be considered. Figure 6 shows the accuracy vs total training time of 2 drXAI configurations vs All Features (training on non-reduced data) for ConvTran (the most expensive classifier), on 4 of the UTSC datasets, where the model can run using All Features. The plots show that the accuracy is comparable, while the total time is reduced by one order of magnitude. The Appendix shows a detailed analysis of this cost. Overall, based on our experiments, SHAP is preferable when aggressive data reduction is the priority, while FA is the best choice under time constraints, offering a faster yet effective selection. 5 Conclusion This paper introduces drXAI, a novel methodology that repurposes XAI at- tribution methods to drive effective data reduction in TSC. By leveraging the GPU-accelerated Hydra classifier and fast explainers, we successfully bridge the gap between XAI interpretability and practical model scalability. Our frame- work leverages two complementary aggregation strategies (absolute and aver- age), identifying the most salient features using an automated elbow-cut heuris- tic, avoiding the need for manual thresholds. We validated drXAI across large-scale synthetic and real-world benchmarks using three SOTA classifiers: MultiRocket-Hydra, InceptionTime, and ConvTran. XAI-Driven Data Reduction15 Fig. 6: Accuracy vs Time (minutes in log scale) for ConvTran on UTSC datasets. Total time compares training on All Features vs drXAI pipeline (Hydra training + explanation + training on reduced data). On 2 of 6 datasets, All Features cannot run due to exceeding GPU memory. Our results demonstrate that drXAI consistently achieves over 80-90% data re- duction while maintaining, and in some cases exceeding, classification accuracy on the full dataset. Specifically, drXAI-SHAP-Proto emerged as the most ac- curate configuration for channel selection in MTSC, while drXAI-FA-Proto led in performance for time-point selection in UTSC. The proposed Proto back- ground marginally outperforms the zero background, at no additional cost; SHAP achieves superior data reduction compared to Feature Ablation but at the expense of longer computation time. Crucially, we demonstrate that data reduction enabled by the lightweight Hydra model allows resource-intensive architectures like ConvTran to scale to massive datasets that were previously inaccessible due to GPU memory con- straints. While limitations exist regarding the computational overhead of SHAP on extremely long sequences, by grouping features drXAI shows that XAI can be a robust, flexible tool for high-performance feature selection. Future work will ex- tend this framework to regression tasks and explore the simultaneous reduction of both channels and time points. Acknowledgments This publication has emanated from research supported in part by a grant from Science Foundation Ireland under Grant number 18/CRT/6183. For the purpose of Open Access, the author has applied a C BY public copyright licence to any Author Accepted Manuscript version arising from this submission. References 1. M. Ashraf et al. A survey on dimensionality reduction techniques for time-series data. IEEE Access, 11, 2023. 2. A. J. Bagnall et al. The UEA multivariate time series classification archive, 2018. CoRR, abs/1811.00075, 2018. 3. T. M. Cox et al. Understanding the impacts of anthropogenic sound on beaked whales. 2006. 16Serramazza et al. 4. A. Dempster et al. Highly scalable time series classification for very large datasets. In International Workshop on Advanced Analytics and Learning on Temporal Data, pages 80–95. Springer, 2024. 5. A. Dempster et al. Monster: Monash scalable time series evaluation repository, 2025. 6. B. Dhariyal, T. Le Nguyen, and G. Ifrim. Scalable classifier-agnostic channel se- lection for multivariate time series classification. Data Mining and Knowledge Discovery, 37(2):1010–1054, 2023. 7. N. M. Foumani et al. Improving position encoding of transformers for multivariate time series classification. Data mining and knowledge discovery, 38(1):22–48, 2024. 8. A. Ismail-Fawaz et al. An approach to multiple comparison benchmark evalu- ations that is stable under manipulation of the comparate set. arXiv preprint arXiv:2305.11921, 2023. 9. H. Ismail Fawaz et al. Inceptiontime: Finding alexnet for time series classification. Data mining and knowledge discovery, 34(6):1936–1962, 2020. 10. N. Kokhlikyan et al. Captum: A unified and generic model interpretability library for pytorch. arXiv preprint arXiv:2009.07896, 2020. 11. S. M. Lundberg and S.-I. Lee. A unified approach to interpreting model predictions. Advances in neural information processing systems, 30, 2017. 12. M. Middlehurst, A. Ismail-Fawaz, A. Guillaume, C. Holder, D. Guijo-Rubio, G. Bu- latova, L. Tsaprounis, L. Mentel, M. Walter, P. Schäfer, and A. Bagnall. aeon: a python toolkit for learning from time series. Journal of Machine Learning Research, 25(289):1–10, 2024. 13. M. Middlehurst, P. Schäfer, and A. Bagnall. Bake off redux: a review and experi- mental evaluation of recent time series classification algorithms. Data Mining and Knowledge Discovery, 2024. 14. T. L. Nguyen and G. Ifrim. Tshap: Fast and exact shap for explaining time series classification and regression. In Joint European Conference on Machine Learning and Knowledge Discovery in Databases, pages 60–77. Springer, 2025. 15. L. Nuyts et al. Tselect: selecting relevant and non-redundant channels for mul- tivariate time series classification: L. nuyts et al. Data Mining and Knowledge Discovery, 39(6):76, 2025. 16. E. Olivetti et al. Meg decoding across subjects. In 2014 international workshop on pattern recognition in neuroimaging, pages 1–4. IEEE, 2014. 17. F. Pedregosa et al. Scikit-learn: Machine learning in Python. 12:2825–2830, 2011. 18. J. Salamon et al. A dataset and taxonomy for urban sound research. In Proceedings of the 22nd ACM international conference on Multimedia, pages 1041–1044, 2014. 19. D. I. Serramazza et al. Improving the evaluation and actionability of explanation methods for multivariate time series classification. In Joint European Confer- ence on Machine Learning and Knowledge Discovery in Databases, pages 177–195. Springer, 2024. 20. D. I. Serramazza et al. A short tutorial for multivariate time series explanation using tscaptum. Software Impacts, 22:100723, 2024. 21. D. I. Serramazza, T. T. Nguyen, T. Le Nguyen, and G. Ifrim. Evaluating explana- tion methods for multivariate time series classification. In International Workshop on Advanced Analytics and Learning on Temporal Data, pages 159–175. Springer, 2023. 22. A. Singh et al. An examination of wearable sensors and video data capture for human exercise classification. In Joint European Conference on Machine Learning and Knowledge Discovery in Databases, pages 312–329. Springer, 2023. XAI-Driven Data Reduction17 23. C. W. Tan et al. Multirocket: multiple pooling operators and transformations for fast and effective time series classification: Cw tan. Data Mining and Knowledge Discovery, 36(5):1623–1646, 2022. 24. A. Theissler et al. Explainable ai for time series classification: A review, taxonomy and research directions. IEEE Access, 10:100700–100724, 2022. 25. H. Turbé, M. Bjelogrlic, C. Lovis, and G. Mengaldo. Evaluation of post-hoc in- terpretability methods in time-series classification. Nature Machine Intelligence, 5(3):250–260, 2023. 26. I. Yousef et al. The arc loss dataset, Feb 2025. 1 Appendix drXAI Dataset Characteristics dataset namen. samples trainn. samples testn. channelsn. time pointsn. classes synthetic MTSC500050004010002 Arc Loss25816459611012 FaceDetection58903524144622 MP1426595501614 Rowing1838790501615 synthetic UTSC100050001200002 CornellWhaleChallenge240006000140002 MosquitoSound55913223653137506 RightWhaleCalls109341962140002 UrbanSound2717271814410010 WhaleSounds8413121032125008 Table 1. Characteristics of all datasets used in the paper. 2 MTSC Accuracy datasetclassifier All Features FA Proto SHAP Proto FA zeros SHAP zeros ECPECSTSelectrandom Arc LossConvT..748.749.754.740.744.745.753.671.725 Arc LossIncepT..726.716.719.729.709.724.736.654.678 Arc LossMRH.740.734.726.726.730.729.733.644.674 Face D.ConvT..560.570.593.596.559.603.621.564.566 Face D.IncepT..677.589.598.610.599.634.639.659.663 Face D.MRH.602.553.561.563.551.594.602.599.592 Military PressConvT..538.599.683.733.707.638.648.243.509 Military PressIncepT..661.527.641.604.575.665.555.261.545 Military PressMRH.7.768.82.807.787.833.78.248.648 RowingConvT..671.774.800.758.71.77.752.553.547 RowingIncepT..605.672.669.675.6.689.729.567.434 RowingMRH.755.729.744.725.719.737.705.568.688 Synthetic MTSCConvT..828.858.693.597.568.714.597.589.53 Synthetic MTSCIncepT..518.527.653.589.551.511.572.524.529 Synthetic MTSCMRH.830.820.628.571.547.677.566.552.53 Table 2. MTSC accuracies. Face D. stands for Face detection, ConvT. for ConvTran and IncepT. for InceptionTime. MTSC Data Selection Details datasetFA ProtoSHAP ProtoFA zerosSHAP zerosECPECSTSelectrandom n%n%n%n%n%n%n%n% 15 Arc Loss140.146120.125120.12580.083380.396590.61510.01150.16 FaceDetection100.06950.03570.049130.09120.083120.083870.6041170.81 Rowing110.22140.2870.1490.18220.44140.2810.02320.64 MP180.36140.2870.1450.1250.5140.2820.0440.08 Table 3. Selection length for real-world MTSC datasets. For each selection of our method and baselines, number of selected channels (n) and relative percentage of the total. 3 UTSC Accuracy datasetclassifier All Features FA Proto SHAP Proto FA zeros SHAP zeros RFIMIrandom CornellWhaleChallengeConvT..897.865.835.872.805.778.892.827 CornellWhaleChallengeIncepT..863.867.843.875.813.783.867.835 CornellWhaleChallengeMRH.838.807.756.825.759.727.832.757 MosquitoSoundConvT..703.726.741.714.726.684.5770.66 MosquitoSoundIncepT..930.819.843.783.843.756.658.892 MosquitoSoundMRH.885.771.794.727.786.691.537.838 RightWhaleCallsConvT..860.867.857.878.792.739.867.831 RightWhaleCallsIncepT..846.871.854.850.784.740.850.811 RightWhaleCallsMRH.810.814.781.819.711.688.799.753 WhaleSoundsConvT..726.666.666.666.666.644.672.694 WhaleSoundsIncepT..691.703.703.703.704.681.713.727 WhaleSoundsMRH.637.638.639.639.637.607.643.634 UrbanSoundConvT.NA.499.437NANA.459.444NA UrbanSoundIncepT..729.627.588.645.674.571.484NA UrbanSoundMRH.817.683.636.707.731.559.515.732 synthetic UTSCConvT.NA.838.921.710.608.598.605.647 synthetic UTSCIncepT..508.813.647.828.589.673.792.571 synthetic UTSCMRH.961.898.957.865.626.810.851.664 Table 4. UTSC accuracies. Experiments where the method exceeded the available memory are denote with NA. ConvT. stands for ConvTran, IncepT. stands for Incep- tionTime. 4 UTSC Data Selection Details datasetFA ProtoSHAP ProtoFA zerosSHAP zerosRFIMIrandom n%n%n%n%n%n%n% CornellWhaleChallenge11990.303990.1012000.308000.203410.0937610.946000.15 MosquitoSound3760.105640.151880.055640.152080.061040.0320650.55 RightWhaleCalls16000.406000.1530000.755990.153550.0936360.911000.25 WhaleSounds2500.102500.102500.102500.102030.084080.1616250.65 UrbanSound44100.1022050.0566150.1588200.2025260.0638040.0988200.2 Table 5. Selection length for real-world UTSC datasets. For each selection of our method and baselines, number of selected channels (n) and relative percentage of the total 5 Time analysis on UTSC data datasetclassifierAF TimedrXAI-FA ProtodrXAI-SHAP Proto hydra timeexp timetrain timeTOTALexp timetrain timeTOTAL RWCCON102.30.30.0411.511.841.042.33.64 RWCINC30.870.30.0411.611.941.044.45.74 RWCMRH5.60.30.042.42.741.041.32.64 CWCCON216.060.60.0412.112.741.025.36.92 CWCINC76.430.60.0425.125.741.026.27.82 CWCMRH15.390.60.047.37.941.025.16.72 MSCON278.563.30.139.713.1338.614.9 MSINC165.383.30.139.813.23317.623.9 MSMRH168.913.30.1376.479.8338187.3 WSCON378.891.20.119.220.52.519.222.9 WSINC201.491.20.111.412.72.511.214.9 WSMRH179.711.20.1140.8142.12.5145.1148.8 USCONNA2.56.5148.657.61150.6411.3164.44 USINC195.572.56.5115.424.41150.6410.4163.54 USMRH31.672.56.512.711.71150.641.3154.44 SYCONNA1.10.6593.194.8514.4988.4103.99 SYINC12.871.10.655.67.3514.495.521.09 SYMRH15.11.10.656.48.1514.497.222.79 Table 6. Detailed time comparison (in minutes) among all features (AF TIME), drXAI-FA Proto and drXAI-SHAP Proto, across the 3 SOTA classifiers Con- vTran (CON), IncepT.ionTime (INC) and MultiRocket-Hydra(MRH) using the UTSC dataset RightWhaleCalls (RWC), CornellWhaleChallenge (CWC), MosquitoSound (MS), WhaleSounds (WS), UrbanSound (US), and synthetic UTSC (SY). To have a lither table, we omit the Hydra time For drXAI-SHAP Proto , as it is identical to the one for drXAI-FA Proto. 6 datasetclassifierdrXAI-FA Proto time/AF TimedrXAI-SHAP Proto time/AF Time RightWhaleCallsConvT.0.1160.036 RightWhaleCallsIncepT.0.3870.186 RightWhaleCallsMRH0.4890.471 CornellWhaleChallengeConvT.0.0590.032 CornellWhaleChallengeIncepT.0.3370.102 CornellWhaleChallengeMRH0.5160.437 MosquitoSoundConvT.0.0470.053 MosquitoSoundIncepT.0.0800.145 MosquitoSoundMRH0.4730.517 WhaleSoundsConvT.0.0540.060 WhaleSoundsIncepT.0.0630.074 WhaleSoundsMRH0.7910.828 UrbanSoundConvT.NANA UrbanSoundIncepT.0.1250.836 UrbanSoundMRH0.3704.877 synthetic UTSCConvT.NANA synthetic UTSCIncepT.0.5711.639 synthetic UTSCMRH0.5401.509 Table 7. Proportion between the total time of drXAI-FA Proto and All Features time (AF Time) as well as the proportion between the total time of drXAI-FA SHAP and All Features time for each SOTA classifier and for each UTSC dataset. We denote with NA where the classifier can’t run using all features due to GPU exceeding the available memory.