<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Neurosci.</journal-id>
<journal-title>Frontiers in Neuroscience</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Neurosci.</abbrev-journal-title>
<issn pub-type="epub">1662-453X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fnins.2025.1610766</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Signal-to-event encoding parameter selection for multiple event classification with spiking neural networks</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Pabian</surname> <given-names>Mateusz</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2802918/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Rzepka</surname> <given-names>Dominik</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/3057337/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Pawlak</surname> <given-names>Miros&#x00142;aw</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Mi&#x0015B;kowicz</surname> <given-names>Marek</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/986178/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Sroka</surname> <given-names>Ryszard</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Measurement and Electronics, AGH University of Krakow</institution>, <addr-line>Krak&#x000F3;w</addr-line>, <country>Poland</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Electrical and Computer Engineering, University of Manitoba</institution>, <addr-line>Winnipeg, MB</addr-line>, <country>Canada</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Milutin Stanacevic, Stony Brook University, United States</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Mingguo Zhao, Tsinghua University, China</p>
<p>Siying Liu, Tsinghua University, China</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Miros&#x00142;aw Pawlak <email>Miroslaw.Pawlak&#x00040;umanitoba.ca</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>23</day>
<month>06</month>
<year>2025</year>
</pub-date>
<pub-date pub-type="collection">
<year>2025</year>
</pub-date>
<volume>19</volume>
<elocation-id>1610766</elocation-id>
<history>
<date date-type="received">
<day>12</day>
<month>04</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>06</day>
<month>06</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2025 Pabian, Rzepka, Pawlak, Mi&#x0015B;kowicz and Sroka.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Pabian, Rzepka, Pawlak, Mi&#x0015B;kowicz and Sroka</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>Event-driven systems can operate either on discrete-time event streams or on analog signals transformed into the event domain by a predefined encoding scheme. This paper studies the problem of optimal event-based signal encoding if data are to be processed by a machine learning model, such as the spiking neural network (SNN). We introduce a method of encoding parameter selection that evaluates a k-Nearest Neighbor (k-NN) classifier operating on a measure of the event stream distance in multiple trials of a Bayesian optimization process. The efficiency of the proposed method is assessed by relating the classification performance with the number of events produced by a signal-to-event encoding scheme. The proposed method is validated for vehicle monitoring sensor data with three event-based encoding schemes: level-crossing encoding, send-on-delta, and leaky integrate-and-fire encoder. The best-performing sets of encoding parameters give an average accuracy of up to 0.912 for the k-NN classification, while producing 97.8% fewer number of samples than for the classical periodic discrete-time signal representation. Additionally, we train the SNN classifiers on data encoded according to the selected sets of parameters, achieving an average classification accuracy of up to 0.946, improving upon the k-NN baseline. This shows that the proposed model-agnostic signal-to-event encoding parameter selection is promising for training sophisticated machine learning models.</p></abstract>
<kwd-group>
<kwd>event-based signal encoding</kwd>
<kwd>van Rossum distance</kwd>
<kwd>Bayesian optimization</kwd>
<kwd>multiple event classification</kwd>
<kwd>spiking neural networks</kwd>
<kwd>k-NN classifier</kwd>
</kwd-group>
<contract-num rid="cn001">327/447</contract-num>
<contract-num rid="cn002">DEC-2017/27/B/ST7/03082</contract-num>
<contract-num rid="cn003">PM/SP/0041/2021/1</contract-num>
<contract-sponsor id="cn001">Natural Sciences and Engineering Research Council of Canada<named-content content-type="fundref-id">https://doi.org/10.13039/501100000038</named-content></contract-sponsor>
<contract-sponsor id="cn002">Narodowe Centrum Nauki<named-content content-type="fundref-id">https://doi.org/10.13039/501100004281</named-content></contract-sponsor>
<contract-sponsor id="cn003">Ministerstwo Edukacji i Nauki<named-content content-type="fundref-id">https://doi.org/10.13039/501100004569</named-content></contract-sponsor>
<counts>
<fig-count count="9"/>
<table-count count="2"/>
<equation-count count="15"/>
<ref-count count="77"/>
<page-count count="17"/>
<word-count count="12553"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Neuromorphic Engineering</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1 Introduction</title>
<p>Event sequence data is a natural representation for modeling of discrete-event systems, in which the state changes occur at discrete points in time due to specific events (Cassandras and Lafortune, <xref ref-type="bibr" rid="B6">2008</xref>). The discrete-event systems (DESs) are commonly used to model and analyze processes that evolve in steps such as the arrival of a data packet, completion of a task, or a request for service. The examples of DESs are manufacturing and computer operating systems (banking, social media, etc.), as well as telecommunication or healthcare systems.</p>
<p>The event representation can also be applied to continuous dynamics present in event-based control and signal processing (Heemels et al., <xref ref-type="bibr" rid="B15">2012</xref>; Tsividis, <xref ref-type="bibr" rid="B68">2003</xref>; Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B40">2015</xref>), especially if implemented within the Industrial Internet of Things (IIoT) (Aranda-Escol&#x000E1;stico et al., <xref ref-type="bibr" rid="B2">2024</xref>). Event-based representations of continuous-time processes focus on capturing key moments in the system evolution, when significant changes&#x02014;considered &#x0201C;events&#x0201D;&#x02014;occur. The event rate for a given continuous input is usually much lower than when that input is represented at regular time intervals (Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B38">2003</xref>). This helps to reduce data transmission and energy consumption in applications such as environmental monitoring, industrial process automation, smart grids and power systems (Aranda-Escol&#x000E1;stico et al., <xref ref-type="bibr" rid="B2">2024</xref>). Event-based techniques enable resource-aware design and are especially effective in resource-constrained applications.</p>
<p>The classical event-based representations of continuous-time signals are level-crossing models where events are produced when a signal crosses a specific threshold or a set of levels (Mark and Todd, <xref ref-type="bibr" rid="B29">1981</xref>), or send-on-delta reporting when events are triggered by signal changes by a prespecified increment (Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B39">2006</xref>). The send-on-delta encoding is used in event-triggered control (Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B40">2015</xref>) and state estimation (Ge et al., <xref ref-type="bibr" rid="B12">2020</xref>), as well as in biologically-inspired event sensors that mimic the function of eyes, ears, nose, or touch (Tayarani-Najaran and Schmuker, <xref ref-type="bibr" rid="B64">2021</xref>; Cheng et al., <xref ref-type="bibr" rid="B7">2019</xref>). On the other hand, the level-crossing representations are applied in event-based signal processing (Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B40">2015</xref>) including signal reconstruction (Rzepka et al., <xref ref-type="bibr" rid="B56">2018</xref>), bandwidth estimation (Rzepka et al., <xref ref-type="bibr" rid="B57">2017</xref>), ECG signal analysis (Ravanshad et al., <xref ref-type="bibr" rid="B54">2013</xref>), design of analog-to-digital converters (Ravanshad et al., <xref ref-type="bibr" rid="B54">2013</xref>; Wang et al., <xref ref-type="bibr" rid="B70">2018</xref>), and digital filters (Huber and Liu, <xref ref-type="bibr" rid="B19">2019</xref>).</p>
<p>Some event-based representations are biologically inspired. One example is the integrate-and-fire encoder (Stein, <xref ref-type="bibr" rid="B60">1965</xref>; Thao et al., <xref ref-type="bibr" rid="B65">2023</xref>), that simplifies real neuron dynamics by mimicking two essential behaviors: integration as a sum of incoming electrical inputs over time, and firing as a generation of spike (event) when the input reaches a certain threshold. The integrate-and-fire model is frequently used in computational neuroscience and machine learning to simulate the behavior of biological neurons in spiking neural networks (SNNs) (Nunes et al., <xref ref-type="bibr" rid="B49">2022</xref>). This enables networks to process information based on the timing of events (spikes) rather than continuous values, offering advantages for temporal data analysis, especially for time-series data or sensory processing tasks.</p>
<p>This paper explores the problem of optimal event-based signal encoding for further processing of event streams by classification models. The proposed approach involves fitting a k-NN classifier to an event stream using the van Rossum distance measure through multiple trials in a Bayesian optimization process. Based on the results of Bayesian optimization, we determine a fixed set of signal-to-event encoding parameters suitable for further processing by more advanced classification techniques.</p>
<p>To validate the proposed methodology, we use sensor data from an Intelligent Transportation System (ITS) for vehicle monitoring (Marsza&#x00142;ek et al., <xref ref-type="bibr" rid="B31">2023</xref>). Our method makes no assumptions about the underlying signal, ensuring its adaptability to other data domains. To demonstrate how the encoding parameters selected by the k-NN classifier can support multiple event classifications, we train a time-to-first-spike spiking neural network (SNN) (Mostafa, <xref ref-type="bibr" rid="B45">2018</xref>) on event-data streams. The efficiency of the proposed method is evaluated by analyzing the relationship between classification performance and the number of events generated by the event-based encoding scheme. This allows us to assess the impact of signal- to-event encoding on classification accuracy for this dataset.</p>
<p>This paper is structured as follows. Section 2 provides a brief overview of event-based signal encoding schemes, spiking neural network training methods, and introduces the technical context of the vehicle classification problem used to evaluate the proposed methodology. Section 3.1 presents the paper&#x00027;s main contribution: a model-agnostic method for selecting signal-to-event encoding parameters. Section 3.2 focuses on the time-to-first-spike SNN model trained on event data encoded according to the proposed signal-to-event encoding schemes. Section 4 details an extensive simulation study assessing the methodology in the context of the vehicle classification problem. Finally, Section 5 summarizes the findings and provides recommendations for future research.</p></sec>
<sec id="s2">
<title>2 Background and related work</title>
<sec>
<title>2.1 Event-based signal encoding</title>
<p>The objective of event-based encoding is to create a framework for efficient representation of a continuous-time signal by capturing its values only at important instants (events), rather than registering it at regular time intervals as in classical methods of signal discretization. The event rate depends on the nature of the input signal, growing for example when the input changes substantially, and becoming smaller during periods of low activity. This allows for minimizing the number of events required to accurately represent the signal, reducing unnecessary computations and data transmissions (Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B38">2003</xref>).</p>
<p>The efficient event-based signal representation preserves signal characteristics by retaining its original key features, such as dynamics and frequency content, while avoiding artifacts or noise contribution that can cause generation of the events too frequently. With a flexible definition of events, event-based encoding of continuous-time signals can meet diverse system design objectives and application needs. Most definitions refer to the event as to a change of the state of an object occurring at an instant (Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B40">2015</xref>).</p>
<p>A class of event-based encoding that may be referred literally to the concept of the event as a significant change of the particular signal parameter are <italic>threshold-based criteria</italic>, known in the context of communication as the <italic>send-on-delta</italic> encoding (Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B39">2006</xref>). In the <italic>send-on-delta</italic> scheme, an input signal is encoded as an event when it deviates by a certain threshold (&#x00394;) (<xref ref-type="fig" rid="F1">Figure 1a</xref>). If <italic>x</italic>(<italic>t</italic>) is a continuous-time signal to be encoded, the send-on-delta encoding times <italic>t</italic><sub><italic>n</italic></sub>&#x02265;0 for <italic>n</italic> &#x0003D; 1, 2, &#x02026;&#x02009; with the threshold &#x00394; are determined such that the following condition holds:</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo class="qopname">min</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x0003E;</mml:mo><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>:</mml:mo><mml:mo>|</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>-</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>|</mml:mo><mml:mo>=</mml:mo><mml:mo>&#x00394;</mml:mo></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Examples of event-based encoding criteria. Black marker symbols denote different event types produced by the encoding. <bold>(a)</bold> Send-on-delta. <bold>(b)</bold> Leaky integrate-and-fire. <bold>(c)</bold> Level-crossing encoding.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-19-1610766-g0001.tif"/>
</fig>
<p>The size of the threshold &#x00394; defines a trade-off between encoding resolution and the event rate. The threshold can be fixed or modified during system operation depending on the safety and/or performance requirements of the system (Tang et al., <xref ref-type="bibr" rid="B63">2025</xref>), or as a response to variation of a signal slope (Trakimas and Sonkusale, <xref ref-type="bibr" rid="B67">2011</xref>; Wang et al., <xref ref-type="bibr" rid="B70">2018</xref>).</p>
<p>A different integration-based encoding scheme is the integrate-and-fire (IF) encoder that mimics the behavior of biological neurons by describing how a potential of a neuron membrane changes in response to input stimuli (Moser et al., <xref ref-type="bibr" rid="B44">2025</xref>). The IF encoder accumulates an input over time until the total input reaches a threshold <italic>V</italic><sub>thr</sub>, which causes emission of a spike and a reset of the accumulator to zero:</p>
<disp-formula id="E2"><label>(2)</label><mml:math id="M2"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo class="qopname">min</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x0003E;</mml:mo><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>:</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003C4;</mml:mi></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">int</mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo class="qopname">&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mi>x</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">thr</mml:mtext></mml:mrow></mml:msub></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003C4;<sub>int</sub> is the integration constant.</p>
<p>In a more biologically realistic neuron model, called leaky integrate-and-fire (LIF) encoder, the content of the accumulator decays exponentially with the time constant &#x003C4;<sub>leak</sub> (<xref ref-type="fig" rid="F1">Figure 1b</xref>) (Thao et al., <xref ref-type="bibr" rid="B65">2023</xref>):</p>
<disp-formula id="E3"><label>(3)</label><mml:math id="M3"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo class="qopname">min</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x0003E;</mml:mo><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>:</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003C4;</mml:mi></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">int</mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo class="qopname">&#x0222B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mo class="qopname">exp</mml:mo><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mfrac><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003C4;</mml:mi></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">leak</mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mi>x</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>V</mml:mi></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">thr</mml:mtext></mml:mrow></mml:msub></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>Another generic class of event-based encoding relevant to many applications in event-based signal processing is not related to detection of an explicit change of the input but focused on matching a certain reference. In the <italic>reference-crossing encoding</italic>, the input is registered when it crosses a prespecified function. The examples are the sine-wave crossing (Selva, <xref ref-type="bibr" rid="B58">2012</xref>) or level-crossing encoding (Mark and Todd, <xref ref-type="bibr" rid="B29">1981</xref>).</p>
<p>In the level-crossing encoding, the reference function can be defined as a single level (Abrahams, <xref ref-type="bibr" rid="B1">1986</xref>), or multiple reference levels (<italic>L</italic><sub>1</sub>, &#x02026;, <italic>L</italic><sub><italic>m</italic></sub>) disposed in the amplitude domain (Mark and Todd, <xref ref-type="bibr" rid="B29">1981</xref>) (<xref ref-type="fig" rid="F1">Figure 1c</xref>).</p>
<p>The level-crossing encoding with the set of levels <italic>L</italic><sub>1</sub>, &#x02026;, <italic>L</italic><sub><italic>m</italic></sub> is defined as follows:</p>
<disp-formula id="E4"><mml:math id="M4"><mml:msub><mml:mi>t</mml:mi><mml:mi>n</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>min</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x0003E;</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>&#x02212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>:</mml:mo></mml:mrow></mml:mrow></mml:math></disp-formula>
<disp-formula id="E5"><label>(4)</label><mml:math id="M"><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable columnalign='left'><mml:mtr columnalign='left'><mml:mtd columnalign='left'><mml:mrow><mml:mi>x</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mtext>or&#x000A0;</mml:mtext><mml:mi>x</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mtext>if&#x000A0;</mml:mtext><mml:msup><mml:mi>x</mml:mi><mml:mo>&#x02032;</mml:mo></mml:msup><mml:mo stretchy='false'>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>&#x02212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo stretchy='false'>)</mml:mo><mml:mo>&#x0003E;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr><mml:mtr columnalign='left'><mml:mtd columnalign='left'><mml:mrow><mml:mi>x</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mtext>or&#x000A0;</mml:mtext><mml:mi>x</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi>L</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x02212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mtext>if&#x000A0;</mml:mtext><mml:msup><mml:mi>x</mml:mi><mml:mo>&#x02032;</mml:mo></mml:msup><mml:mo stretchy='false'>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>&#x02212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo stretchy='false'>)</mml:mo><mml:mo>&#x0003C;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:math></disp-formula>
<p>where <italic>L</italic><sub><italic>k</italic></sub> and <italic>L</italic><sub><italic>k</italic>&#x0002B;1</sub> are reference levels, <italic>L</italic><sub><italic>k</italic>&#x0002B;1</sub>&#x0003E;<italic>L</italic><sub><italic>k</italic></sub>, 1 &#x0003C; <italic>k</italic>&#x0003C;<italic>m</italic>.</p>
<p>This condition specifies that the next level crossing after the actual crossing of <italic>L</italic><sub><italic>k</italic></sub> will occur for the same level (<italic>L</italic><sub><italic>k</italic></sub>) or the higher level (<italic>L</italic><sub><italic>k</italic>&#x0002B;1</sub>) in case of the up-crossing (<inline-formula><mml:math id="M6"><mml:msup><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>&#x0003E;</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula>) of <italic>L</italic><sub><italic>k</italic></sub> at <italic>t</italic><sub><italic>n</italic>&#x02212;1</sub>, or for the same level (<italic>L</italic><sub><italic>k</italic></sub>) or the lower level (<italic>L</italic><sub><italic>k</italic>&#x02212;1</sub>) in case of the down-crossing (<inline-formula><mml:math id="M7"><mml:msup><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>&#x0003C;</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula>) of <italic>L</italic><sub><italic>k</italic></sub> at <italic>t</italic><sub><italic>n</italic>&#x02212;1</sub>.</p>
<p>The level-crossing encodings provide not only time instants when the relevant reference levels are crossed but also imply that the signal remains between consecutive levels during the time between the level crossings (<xref ref-type="fig" rid="F1">Figure 1c</xref>). This means that for <italic>t</italic>&#x02208;(<italic>t</italic><sub><italic>n</italic>&#x02212;1</sub>, <italic>t</italic><sub><italic>n</italic></sub>):</p>
<disp-formula id="E6"><label>(5)</label><mml:math id="M8"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0003C;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>&#x0003C;</mml:mo><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mtext class="textrm" mathvariant="normal">if</mml:mtext><mml:msup><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>&#x0003E;</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x0003C;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>&#x0003C;</mml:mo><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mtext class="textrm" mathvariant="normal">if</mml:mtext><mml:msup><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>&#x0003C;</mml:mo><mml:mn>0</mml:mn><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>The extra knowledge defined by the above inequalities, called the implicit information, implies that no event occurs between consecutive encoding instants (Rzepka et al., <xref ref-type="bibr" rid="B56">2018</xref>). The implicit information is relevant not only to the level-crossing encodings but also to any event-based encoding scheme.</p>
<p>The reference levels in the level-crossing encoding are usually uniformly distributed in the amplitude domain although adaptive (Senay et al., <xref ref-type="bibr" rid="B59">2010</xref>) or optimal distributions of levels (Kozat et al., <xref ref-type="bibr" rid="B25">2013</xref>) have been also proposed.</p>
<p>The mean rate of level-crossing and send-on-delta encodings depends on the average slope of the signal (Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B39">2006</xref>), while the rate of spikes in IF encodings is defined by the average values of the input signal. If the hysteresis is adopted to triggering (multiple) level crossings (i.e., the repeated crossings of the same level are not encoded), then such level-crossing criterion coincides with the threshold-based encoding (send-on-delta scheme).</p>
<p>The general model of threshold-based encoding that provides a unified framework for send-on-delta, leaky integrate-and-fire, and other threshold-based schemes (e.g., leaky send-on-delta) has been introduced in Moser (<xref ref-type="bibr" rid="B43">2017</xref>).</p>
<p>The threshold-based event encoding is a fundamental paradigm in event-based control (Heemels et al., <xref ref-type="bibr" rid="B15">2012</xref>; Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B40">2015</xref>) and state estimation (Ge et al., <xref ref-type="bibr" rid="B12">2020</xref>) designed for efficient utilization of computation and communication resources, especially when implemented with wireless connectivity in IIoT (Aranda-Escol&#x000E1;stico et al., <xref ref-type="bibr" rid="B2">2024</xref>; Mi&#x0015B;kowicz, <xref ref-type="bibr" rid="B39">2006</xref>). The send-on-delta encoding is applied to event sensors that mimic the function of eyes (silicon retina), ears (silicon cochlea), nose (e-nose), or touch (e-skin) and emit events when a signal representing relevant modality (vision, sound, olfaction, touch) changes by a prespecified threshold (Tayarani-Najaran and Schmuker, <xref ref-type="bibr" rid="B64">2021</xref>; Cheng et al., <xref ref-type="bibr" rid="B7">2019</xref>). On the other hand, the level-crossing encoding is an effective signal representation in energy-efficient signal processing applications (e.g., ECG analysis with wearable devices) (Ravanshad et al., <xref ref-type="bibr" rid="B54">2013</xref>).</p>
<p>Integrate-and-fire encoders as simplified models of neuronal activity are foundational in computational neuroscience and extensively used in the development of SNNs (Thao et al., <xref ref-type="bibr" rid="B65">2023</xref>; Moser et al., <xref ref-type="bibr" rid="B44">2025</xref>). The SNNs are particularly well-suited for processing temporal data due to their inherent ability to handle time-dependent information (e.g., speech recognition, real-time sensory processing, and dynamic vision) (Nunes et al., <xref ref-type="bibr" rid="B49">2022</xref>). One of the active research areas focuses on improving the efficiency of encoding techniques to enhance the performance and applicability of SNNs in various domains (Auge et al., <xref ref-type="bibr" rid="B3">2021</xref>).</p>
<p>A recent work (Zanoli et al., <xref ref-type="bibr" rid="B74">2025</xref>) has shown that optimization-based level-crossing encoding parameter selection methodologies achieve better signal reconstruction in ECG signal monitoring. In this paper, we analyze signal-to-event encoding schemes by examining how signal representation affects the performance of a machine learning model that processes event data using an event distance measure within the Bayesian optimization process. The analysis of event-based encoding schemes is conducted using sensor data for vehicle classification in an intelligent transportation system.</p>
</sec>
<sec>
<title>2.2 Spiking neural networks (SNN)</title>
<p>Event stream data can be processed by specialized machine learning models such as the spiking neural networks (SNN). These models process data using impulses that propagate asynchronously through the network (Pfeiffer and Pfeil, <xref ref-type="bibr" rid="B52">2018</xref>). This approach allows SNNs to mimic biological networks more closely than traditional neural networks. While the development of the training rules for SNNs is an active area of research, several major training paradigms have been identified: plasticity-aware training (Falez et al., <xref ref-type="bibr" rid="B9">2019</xref>; Mozafari et al., <xref ref-type="bibr" rid="B46">2019</xref>); network conversion, which maps each component of the source network to its spiking equivalent (Rueckauer et al., <xref ref-type="bibr" rid="B55">2017</xref>; Midya et al., <xref ref-type="bibr" rid="B37">2019</xref>; St&#x000F6;ckl and Maass, <xref ref-type="bibr" rid="B61">2021</xref>); and training with backpropagation (Lee et al., <xref ref-type="bibr" rid="B27">2016</xref>; Wu et al., <xref ref-type="bibr" rid="B71">2018</xref>; Rasmussen, <xref ref-type="bibr" rid="B53">2019</xref>). The latter approach leverages existing algorithms and best practices developed for deep learning with nonspiking neural networks.</p>
<p>Backpropagation-based training of SNNs must address the challenge that the spike-generating function is non-differentiable. The methods to overcome this problem are of two types, depending on how much information from the forward pass is needed to compute the surrogate gradient signal. In event-driven learning the error is propagated only through spikes (Zhang and Li, <xref ref-type="bibr" rid="B75">2020</xref>; Wunderlich and Pehle, <xref ref-type="bibr" rid="B72">2021</xref>; Zhu et al., <xref ref-type="bibr" rid="B77">2022</xref>). In contrast, in a learning scheme similar to that used in the Recurrent Neural Networks (RNNs) the error information is propagated also through time steps which did not elicit a spike (Wu et al., <xref ref-type="bibr" rid="B71">2018</xref>; Neftci et al., <xref ref-type="bibr" rid="B48">2019</xref>; Xing et al., <xref ref-type="bibr" rid="B73">2020</xref>; Bauer et al., <xref ref-type="bibr" rid="B4">2023</xref>). Training models with these methods usually requires simulating the state of the entire network over a finite time window with a fixed time step. However, it is possible to design learning rules that do not require such extensive simulation. One example is the time-to-first-spike SNN (Mostafa, <xref ref-type="bibr" rid="B45">2018</xref>; Kheradpisheh and Masquelier, <xref ref-type="bibr" rid="B22">2020</xref>; Zhou et al., <xref ref-type="bibr" rid="B76">2021</xref>). This model is trained using locally exact gradients of the spike-generating function, making both the forward- and backward-pass through the network event-centric. Recently, the time-to-first-spike SNN model has been extended to the scenario where each postsynaptic neuron can generate multiple events in response to observed spike trains (Pabian et al., <xref ref-type="bibr" rid="B51">2024</xref>).</p>
</sec>
<sec>
<title>2.3 Inductive loop vehicle magnetic profiles (VMP)</title>
<p>The goal of Intelligent Transportation Systems (ITS) is to enhance the efficient use of existing transportation infrastructure by leveraging accurate traffic data acquisition technologies to monitor and manage traffic flow. Despite recent advancements in sensor technology, inductive loop (IL) sensors remain the most widely used in modern traffic control systems (Klein et al., <xref ref-type="bibr" rid="B24">2006</xref>). This technology features low installation costs, a long lifespan due to its simple design, robustness to weather conditions such as rain, fog or snow, and a flexible architecture that can be adapted to various applications. Systems based on IL sensors are capable of vehicle classification (Gajda et al., <xref ref-type="bibr" rid="B11">2001</xref>; Oh et al., <xref ref-type="bibr" rid="B50">2002</xref>; Ki and Baik, <xref ref-type="bibr" rid="B23">2006</xref>; Jeng and Ritchie, <xref ref-type="bibr" rid="B20">2008</xref>; Meta and Cinsdikici, <xref ref-type="bibr" rid="B36">2010</xref>), vehicle re-identification and tracking (Kwon, <xref ref-type="bibr" rid="B26">2005</xref>; Ndoye et al., <xref ref-type="bibr" rid="B47">2011</xref>; Guilbert et al., <xref ref-type="bibr" rid="B14">2013</xref>, <xref ref-type="bibr" rid="B13">2014</xref>), speed estimation (Sun and Ritchie, <xref ref-type="bibr" rid="B62">1999</xref>; Coifman et al., <xref ref-type="bibr" rid="B8">2003</xref>; Lu et al., <xref ref-type="bibr" rid="B28">2012</xref>), as well as wheel and axle detection (Gajda et al., <xref ref-type="bibr" rid="B10">2012</xref>; Marsza&#x00142;ek et al., <xref ref-type="bibr" rid="B33">2011</xref>, <xref ref-type="bibr" rid="B34">2015</xref>, <xref ref-type="bibr" rid="B35">2018</xref>).</p>
<p>At minimum, an IL detector consists of two components: a wire loop with one or more turns mounted on or embedded in the roadway pavement, and a controller cabinet that houses an electronics unit (Klein et al., <xref ref-type="bibr" rid="B24">2006</xref>). When a ferromagnetic or conducting metallic mass passes over the IL connected to a conditioning circuit, the currents induced in the object alter the magnetic field distribution affecting the loop inductance and resistance (Mochol&#x000ED;-Salcedo et al., <xref ref-type="bibr" rid="B42">2017</xref>). These dynamic interactions between the loop and the vehicle can be monitored and analyzed as impedance changes, generating a waveform signal known as the vehicle magnetic profile (VMP), where R-VMP represents its resistive component and X-VMP its reactive component.</p>
<p>Recently, a quad loop, four-channel VMP measurement system was proposed by Marsza&#x00142;ek and Duda (<xref ref-type="bibr" rid="B30">2020</xref>). It consists of two standard loops (IL1, IL3) and two slim loops (IL2, IL4) arranged in series. The standard loops are well-suited for general-purpose tasks, while the slim loops are preferable for axle identification and wheel rim detection (Gajda et al., <xref ref-type="bibr" rid="B10">2012</xref>; Marsza&#x00142;ek and Sroka, <xref ref-type="bibr" rid="B32">2016</xref>). Using two loops of each type instead of a single one introduces redundancy enhancing the system robustness against signal interference and simplifying the estimation of traffic parameters, such as vehicle speed (Mochol&#x000ED; Belenguer et al., <xref ref-type="bibr" rid="B41">2019</xref>). Additionally, three different excitation frequencies are applied to each channel. The simultaneous multi-frequency measurement increases the system resilience to noise, reducing the risk of poor signal quality affecting downstream processing tasks. <xref ref-type="fig" rid="F2">Figure 2</xref> presents exemplary VMP profiles acquired by this measurement system. We use the quad loop sensor VMP signals originally discretized at 1 kHz in our event-driven signal analysis. A detailed description of the VMP data used in this study is provided in Section 4.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Exemplary VMP profiles acquired by the quad loop VMP measurement system, divided into separate IL sensors, VMP components and loop excitation frequencies. Note that signals originating from the IL sensors of the same type (IL1-IL3; IL2-IL4) are slightly shifted in time. Furthermore, the loop excitation frequency impacts the signal magnitude, but does not significantly change its characteristic features.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-19-1610766-g0002.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="materials and methods" id="s3">
<title>3 Materials and methods</title>
<p>In this section we describe the proposed generic framework for encoding parameter selection, which is the main contribution of this paper. Furthermore, in order to show that the parameters selected by the proposed method can be reused when training more complex models, we also introduce the SNN model used in the experimental validation study.</p>
<sec>
<title>3.1 Model-agnostic signal-to-event encoding parameter selection</title>
<p>The main goal of this paper is to identify a set of signal-to-event encoding parameters that preserve information relevant to a multiple event classification task. We operate under the assumption that, for a simple classification model such as the k-NN classifier, the choice of spike encoding parameters has a measurable impact on classification performance and, intuitively, on the pairwise similarity of the generated event data. Consequently, event data that can be accurately classified by the k-NN model are expected to perform well as input to other event-based classifiers, assuming a positive correlation in their outcomes. This makes our method model-agnostic as it does not rely on any assumptions regarding the specific classifier to be used.</p>
<p><xref ref-type="fig" rid="F3">Figure 3</xref> illustrates the proposed methodology for selecting signal-to-event encoding parameters. The approach involves performing V-fold cross-validation using a k-NN classifier that leverages an event-based distance metric to classify signals from the validation subsets. Within each data fold, encoding parameters are sampled via a Bayesian optimization process, guided by the classification accuracy observed in prior iterations. The resulting scores are aggregated into ranked lists (one per data fold) along with their corresponding encoding parameters. The final set of encoding parameters is then selected based on these rankings. Detailed descriptions of the proposed methodology are provided in Sections 3.1.1-3.1.3. It is important to note that if the distance metric used by the k-NN classifier is itself parameterized, its impact on the classification accuracy must be marginalized, e.g., by sampling multiple distance parameter values for each fixed set of signal encoding parameters during the Bayesian optimization process.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Flowchart of the proposed model-agnostic signal-to-event encoding parameter selection. The k-NN classifier accuracy obtained at various optimization steps is aggregated into ranked lists, along with the corresponding encoding parameters used to encode signals into their event-based representation. The outer loop iterates over different data splits (V-fold cross-validation), while the inner loop performs sampling of encoding parameters for each split using Bayesian optimization. The dotted line representing classifier accuracy indicates that this value serves as the objective function, guiding the parameter sampling in the subsequent iteration of the Bayesian optimization process.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-19-1610766-g0003.tif"/>
</fig>
<sec>
<title>3.1.1 Van Rossum distance</title>
<p>A crucial component of the strategy for selection of signal-to-event encoding parameters based on the k-NN classifier is the choice of an appropriate distance metric to assess the similarity between event sequences. The choice of the distance measure heavily depends on the properties of the original signal. For relatively short signals (such as the VMP), typically spanning only a few seconds, the distance measures commonly used in computational neuroscience are well-suited for processing this type of data (Kass et al., <xref ref-type="bibr" rid="B21">2014</xref>).</p>
<p>Due to its prominence in neuroscience, we opt to use the multi-neuron van Rossum distance (van Rossum, <xref ref-type="bibr" rid="B69">2001</xref>; Houghton and Sen, <xref ref-type="bibr" rid="B18">2008</xref>) to assess the similarity of the event sequences produced by spike-encoding functions<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref>. It is a purely deterministic distance measure and does not take into account the stochastic nature of event data. Furthermore, it assumes that there is no functional dependence between events from different channels.</p>
<p>Let <inline-formula><mml:math id="M9"><mml:mrow><mml:mi mathvariant="script">U</mml:mi></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>u</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>u</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>u</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo>}</mml:mo></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M10"><mml:mrow><mml:mi mathvariant="script">V</mml:mi></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>v</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>v</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>v</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo>}</mml:mo></mml:mrow></mml:math></inline-formula> be two populations of <italic>P</italic> independent neurons. Furthermore, let</p>
<disp-formula id="E7"><label>(6)</label><mml:math id="M11"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>f</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>u</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munder class="msub"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:munder></mml:mstyle><mml:mi>h</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:mi>P</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>be a spike train with the event sequence <inline-formula><mml:math id="M12"><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo>}</mml:mo></mml:mrow></mml:math></inline-formula> in the <italic>p</italic>-th spike train of the population <inline-formula><mml:math id="M13"><mml:mrow><mml:mi mathvariant="script">U</mml:mi></mml:mrow></mml:math></inline-formula>. This spike train is smoothed by a causal exponential kernel</p>
<disp-formula id="E8"><label>(7)</label><mml:math id="M14"><mml:mi>h</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:mi>t</mml:mi><mml:mo>;</mml:mo><mml:mi>&#x003C4;</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable columnalign='left'><mml:mtr columnalign='left'><mml:mtd columnalign='left'><mml:mn>0</mml:mn></mml:mtd><mml:mtd columnalign='left'><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x0003C;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr><mml:mtr columnalign='left'><mml:mtd columnalign='left'><mml:mrow><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x02212;</mml:mo><mml:mi>t</mml:mi><mml:mo>/</mml:mo><mml:mi>&#x003C4;</mml:mi></mml:mrow></mml:msup></mml:mrow></mml:mtd><mml:mtd columnalign='left'><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x02265;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow><mml:mo>.</mml:mo></mml:math></disp-formula>
<p>The spike train <italic>f</italic>(<italic>t</italic>, <bold>v</bold><sup><italic>p</italic></sup>) with the event sequence <inline-formula><mml:math id="M15"><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo>}</mml:mo></mml:mrow></mml:math></inline-formula> is defined analogously. Note that the event sequences <inline-formula><mml:math id="M16"><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo>}</mml:mo></mml:mrow></mml:math></inline-formula> need not have the same number of events. Then, the multi-neuron van Rossum distance (Houghton and Kreuz, <xref ref-type="bibr" rid="B17">2012</xref>) between <inline-formula><mml:math id="M17"><mml:mrow><mml:mi mathvariant="script">U</mml:mi></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M18"><mml:mrow><mml:mi mathvariant="script">V</mml:mi></mml:mrow></mml:math></inline-formula> is</p>
<disp-formula id="E9"><label>(8)</label><mml:math id="M19"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>d</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">U</mml:mi></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:mi mathvariant="script">V</mml:mi></mml:mrow><mml:mo>;</mml:mo><mml:mi>&#x003C4;</mml:mi><mml:mo>,</mml:mo><mml:mi>c</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:mi>c</mml:mi><mml:mstyle displaystyle="true"><mml:munder class="msub"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>q</mml:mi><mml:mo>&#x02260;</mml:mo><mml:mi>p</mml:mi></mml:mrow></mml:munder></mml:mstyle><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mrow></mml:msqrt><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>with</p>
<disp-formula id="E10"><label>(9)</label><mml:math id="M20"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munder class="msub"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:munder></mml:mstyle><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mo>|</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mo>|</mml:mo><mml:mo>/</mml:mo><mml:mi>&#x003C4;</mml:mi></mml:mrow></mml:msup><mml:mo>&#x0002B;</mml:mo><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mo>|</mml:mo><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mo>|</mml:mo><mml:mo>/</mml:mo><mml:mi>&#x003C4;</mml:mi></mml:mrow></mml:msup><mml:mo>-</mml:mo><mml:mn>2</mml:mn><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mo>|</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mo>|</mml:mo><mml:mo>/</mml:mo><mml:mi>&#x003C4;</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>being the <italic>labeled line</italic> term representing the single-neuron van Rossum distance (i.e., a distance between neurons that directly correspond to one another), and</p>
<disp-formula id="E11"><label>(10)</label><mml:math id="M21"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munder class="msub"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:munder></mml:mstyle><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mo>|</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mi>q</mml:mi></mml:mrow></mml:msubsup><mml:mo>|</mml:mo><mml:mo>/</mml:mo><mml:mi>&#x003C4;</mml:mi></mml:mrow></mml:msup><mml:mo>&#x0002B;</mml:mo><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mo>|</mml:mo><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mi>q</mml:mi></mml:mrow></mml:msubsup><mml:mo>|</mml:mo><mml:mo>/</mml:mo><mml:mi>&#x003C4;</mml:mi></mml:mrow></mml:msup><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mo>|</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mi>q</mml:mi></mml:mrow></mml:msubsup><mml:mo>|</mml:mo><mml:mo>/</mml:mo><mml:mi>&#x003C4;</mml:mi></mml:mrow></mml:msup><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mo>|</mml:mo><mml:msubsup><mml:mrow><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mi>q</mml:mi></mml:mrow></mml:msubsup><mml:mo>|</mml:mo><mml:mo>/</mml:mo><mml:mi>&#x003C4;</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>is the <italic>summed population</italic> term representing the cross-neuron distance. Assuming that the sequences in populations <inline-formula><mml:math id="M22"><mml:mrow><mml:mi mathvariant="script">U</mml:mi></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M23"><mml:mrow><mml:mi mathvariant="script">V</mml:mi></mml:mrow></mml:math></inline-formula> are sorted, the computational complexity of <xref ref-type="disp-formula" rid="E9">Equation 8</xref> is <inline-formula><mml:math id="M24"><mml:mrow><mml:mi mathvariant="script">O</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msup><mml:mrow><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:math></inline-formula>, where &#x003BC; is the longest event sequence across populations <inline-formula><mml:math id="M25"><mml:mrow><mml:mi mathvariant="script">U</mml:mi></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:mi mathvariant="script">V</mml:mi></mml:mrow></mml:math></inline-formula>; and can be reduced to <inline-formula><mml:math id="M26"><mml:mrow><mml:mi mathvariant="script">O</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:math></inline-formula> using the so-called markage trick (Houghton and Kreuz, <xref ref-type="bibr" rid="B17">2012</xref>). Note that <xref ref-type="disp-formula" rid="E10">Equations 9</xref>, <xref ref-type="disp-formula" rid="E10">10</xref> can be computed even when one event sequence is empty while the other is not, as well as when both sequences are empty.</p>
<p>The multi-neuron van Rossum distance defined in <xref ref-type="disp-formula" rid="E9">Equation 8</xref> is parameterized by <italic>c</italic> and &#x003C4;, which model various phenomena observed in neuroscience. The mixing parameter 0 &#x02264; <italic>c</italic> &#x02264; 1 weighs the importance of treating each neuron separately (<italic>c</italic> &#x0003D; 0) vs. viewing the entire population as a single-unit (<italic>c</italic> &#x0003D; 1). The decay constant &#x003C4; influences the range of interspike dependencies. For &#x003C4; &#x02192; &#x0221E; each individual spike contributes to the distance computed at each subsequent spike. Conversely, for &#x003C4; &#x02192; 0 spikes impact only their direct neighborhood with &#x003C4; &#x0003D; 0 causing the measure to count the number of coinciding events in the spike trains. These properties show the versatility of this distance function. Unfortunately, this also means that the choice of <italic>c</italic>, &#x003C4; impacts the perceived (dis)similarity of spike trains produced by different spike encoding functions.</p></sec>
<sec>
<title>3.1.2 Encoding parameter selection with Bayesian optimization</title>
<p>To identify the optimal set of parameters for the encoding functions, we employ Bayesian optimization using the Tree-structured Parzen Estimator (TPE) approach, implemented via the Hyperopt software package (Bergstra et al., <xref ref-type="bibr" rid="B5">2013</xref>). Bayesian optimization is a strategy for locating the extrema of an objective function that is expensive to evaluate. It constructs a surrogate model of the objective function based on prior evaluations at sampled points within the parameter space. This model is iteratively refined with each new observation, and subsequent sampling is guided by an acquisition function derived from the surrogate model. The key distinction between Bayesian optimization and a purely random search lies in the sampling strategy: while random search samples the parameter space uniformly, Bayesian optimization leverages prior knowledge by preferentially sampling regions that are more likely to yield improvements, based on past observations.</p>
<p>For a predetermined training and validation data split, the procedure for selecting signal-to-event encoding parameters using Bayesian optimization can be summarized as follows:</p>
<list list-type="order">
<list-item><p>Apply min-max normalization to scale each signal individually to the range 0&#x02013;1. This normalization simplifies the process of selecting encoding parameters by making the resulting event stream independent of the original signal&#x00027;s amplitude range.</p></list-item>
<list-item><p>Define the search space for the parameters of the objective function to be optimized. The dimensionality of the optimization search space depends on the selected spike encoding scheme and the choice of the distance measure. In case of the van Rossum distance the search space must include the parameters &#x003C4; and <italic>c</italic>, which we constrain to the search space log(&#x003C4;)&#x0007E;<bold>U</bold>(&#x02212;1, 1) and <italic>c</italic>&#x0007E;<bold>U</bold>(0, 1), where <bold>U</bold>(<italic>a, b</italic>) denotes the uniform random variable on the interval [<italic>a, b</italic>]. It is important to choose the parameter bounds for the encoding functions such that they are able to produce informative event streams&#x02014;too few events in the sequence might not carry enough information, whereas too many events might encode contributions from artifacts or noise components in the signal that should be ignored.</p></list-item>
<list-item><p>Define the objective function used in the Bayesian optimization process. This function is the classification accuracy of the k-NN algorithm. For each individual record in the validation dataset, the classifier assigns a class label by identifying <italic>k</italic> &#x0003D; 7 most similar signals from the training dataset based on the multi-neuron van Rossum distance and applying majority voting. The number of the nearest neighbors <italic>k</italic> is fixed to avoid introducing an additional parameter into the optimization search space.</p></list-item>
<list-item><p>Run the Bayesian optimization process over 150 iterations. Each iteration involves sampling a new set of parameters (from the search space defined in step 2), extracting event sequences from the input signal, evaluating the objective function, and updating the surrogate model of the objective function conditioned on these parameters. Note that parameter sets which result in empty event sequences are considered invalid and are excluded from subsequent analysis.</p></list-item>
</list>
<p>Finally, to assess the robustness of the optimization procedure with respect to the data distribution, we apply the stratified 10-fold cross-validation, resulting in a total of 10 &#x000D7; 150 distinct optimization steps for each spike encoding scheme.</p>
<p>Note that we limit the scope of our optimization process to a single objective: maximizing the k-NN classification accuracy. One might consider using a different random search type algorithm to optimize a multi-objective function that also considers the event density to selectively promote solutions that achieve high classification performance while producing fewer events. We leave this topic for further research.</p></sec>
<sec>
<title>3.1.3 Weighted median parameter selection strategy</title>
<p>The results of the Bayesian optimization process conducted over V different data splits yield several independent ranked lists <italic>Q</italic> &#x0003D; {<italic>q</italic><sub>1</sub>, <italic>q</italic><sub>2</sub>, &#x02026;, <italic>q</italic><sub><italic>V</italic></sub>}. In each ranking <italic>q</italic>, higher position is assigned to parameter sets that achieve better k-NN classifier scores (ties are permitted). Our objective is to determine a single parameter set (per encoding type) from these rankings that can be reliably used by other multiple event classification algorithms operating on the same data. One possible approach would be to select the parameter set with the highest overall performance across all data splits. However, the absolute score depends not only on the selected parameters&#x02014;which may not have been selected in other splits during the Bayesian optimization<xref ref-type="fn" rid="fn0002"><sup>2</sup></xref>&#x02014;but also on the specific stratified data sample. Moreover, it is preferable to identify the parameter set that performs consistently well across all data splits. Therefore, a selection strategy based solely on top-performing scores in individual rankings is inadvisable.</p>
<p>To address these concerns, we propose the <italic>weighted median selection strategy</italic>&#x02014;a parameter selection method that accounts for the local (i.e., within-ranking) rank of each result across all data splits. The approach is outlined in <xref ref-type="table" rid="T3">Algorithm 1</xref>. Each set of parameters sampled during the Bayesian optimization process is assigned a rank reciprocal score <italic>r</italic> based on its position in the corresponding ranking <italic>q</italic>. The ranking lists is then aggregated, and the rank-reciprocal scores are normalized to produce a weight vector <bold>w</bold><sub>&#x003B1;</sub>, which assigns a weight to each sampled parameter set. This weight vector is used to compute a weighted median of the rank reciprocal scores, defining the final set of encoding parameters returned by the proposed method. If multiple parameter sets share the same rank-reciprocal score equal to the weighted median, the final selection is the median of all such sets. Overall, this strategy fulfills the objective of selecting a parameter set that performs consistently well across different data splits. It achieves this by mitigating the influence of outliers in individual rankings and by emphasizing configurations associated with higher k-NN classifier accuracy through the use of rank-reciprocal weighting.</p>
<table-wrap position="float" id="T3">
<label>Algorithm 1</label>
<caption><p>Weighted median selection strategy for the parameter sets found by the Bayesian optimization procedure.</p></caption>
<graphic xlink:href="fnins-19-1610766-i0001.tif"/>
</table-wrap>
</sec>
<sec>
<title>3.1.4 Computational complexity and problem-specific adaptations</title>
<p>It is important to note that the proposed methodology might require problem-specific adjustments as one cannot guarantee that the optimization hyperparameters related to the number of iterations, the number of data splits in the V-fold cross-validation, or the number of nearest neighbors in the k-NN algorithm suits all possible research problems. Fortunately, adjusting these settings does not meaningfully change the proposed methodology. Furthermore, the entire approach is modular enough that it is possible to make bigger changes to some parts of the algorithm without impacting the general framework summarized in <xref ref-type="fig" rid="F3">Figure 3</xref>. We have already mentioned that there might be some benefit to using a different distance measure than the van Rossum distance, one that is better suited to the data under analysis; or to using a different algorithm for the Bayesian optimization procedure.</p>
<p>The computational complexity of the proposed method depends on three main components: the distance measure computation (as established in Section 3.1.2 the van Rossum distance has the complexity <inline-formula><mml:math id="M30"><mml:mrow><mml:mi mathvariant="script">O</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:math></inline-formula>, where <italic>P</italic> is the number of channels of the signal and &#x003BC; is the number of events in longest event sequence), the iteration over data points in k-NN classification (with a brute force algorithm complexity of <inline-formula><mml:math id="M31"><mml:mrow><mml:mi mathvariant="script">O</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi></mml:mrow></mml:msub><mml:mi>k</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:math></inline-formula>, where <italic>N</italic><sub><italic>t</italic></sub> is the number of training samples and <italic>N</italic><sub><italic>v</italic></sub> is the number of validation samples in the given data split), and the choice of the signal-to-spike encoding function. This means that for some large-scale problems (e.g., when the dataset has many examples or the signal varies enough that the produced event sequences have numerous events) the proposed methodology might be too expensive to compute and an alternative approach should be devised. Note that the signal must be encoded as an event stream regardless of the chosen optimization method, hence its contribution to the computational complexity analysis should be excluded from comparison.</p>
</sec>
</sec>
<sec>
<title>3.2 Time-to-first-spike SNN classifier</title>
<p>To assess how the set of encoding parameters selected by our methodology can be reused to train a more sophisticated event classification model, we train a spiking neural network. Specifically, we use the multiple-input, multiple-output (MIMO) time-to-first-spike SNN introduced by Pabian et al. (<xref ref-type="bibr" rid="B51">2024</xref>). This model extends the work by Mostafa (<xref ref-type="bibr" rid="B45">2018</xref>) by relaxing the implicit assumption of an infinitely-long refractory period &#x003C4;<sub>ref</sub>.</p>
<p>We introduce a <italic>neuron response delay factor</italic> &#x003C4;<sub>delay</sub> in the first hidden layer of the MIMO SNN. This parameter defines the earliest time at which a given neuron is capable of responding to input events. Functionally, this is equivalent to initializing the neuron in a refractory state lasting &#x003C4;<sub>delay</sub>. We hypothesize that assigning a different &#x003C4;<sub>delay</sub> value to each neuron in the layer enables them to observe slightly different event sequences, thereby mitigating the issue of early events disproportionally influencing the training process. Let us consider the following scenario: assume that all neurons start in a resting state, i.e., &#x003C4;<sub>delay</sub> &#x0003D; 0 for all neurons. Additionally, assume that one neuron becomes specialized through training to detect a pattern occurring relatively late in the input event sequence. In this case, this neuron must either adopt smaller synaptic weights to avoid firing prematurely, or maintain comparable weights to other neurons, while producing spikes that are not informative (i.e., spikes triggered by early, irrelevant events). Such behavior arises when a neuron reacts to parts of the input sequence that are not aligned with its intended function.</p>
<p><xref ref-type="fig" rid="F4">Figure 4</xref> presents an example raster plot of a network trained with nonzero &#x003C4;<sub>delay</sub> values assigned to neurons in the first layer, increasing linearly across the layer. We do not apply &#x003C4;<sub>delay</sub> to neurons in deeper layers as the SNN-generated event sequences contain significantly fewer events than the network input, simplifying neuron specialization in those layers.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Spike raster plot of a MIMO SNN model trained with a nonzero neuron response delay factor &#x003C4;<sub>delay</sub> in the first hidden layer. The dashed line denotes the &#x003C4;<sub>delay</sub> value assigned to each neuron in the layer. It is evident that neurons with shorter response delays elicit their first spikes earlier than those with longer &#x003C4;<sub>delay</sub> values.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-19-1610766-g0004.tif"/>
</fig>
<p>The MIMO SNN models are trained by minimizing the following risk function</p>
<disp-formula id="E12"><label>(11)</label><mml:math id="M32"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">total</mml:mtext></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>y</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>y</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:mi>&#x003B3;</mml:mi><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">spiking</mml:mtext></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:mo>&#x003BB;</mml:mo><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where</p>
<disp-formula id="E13"><label>(12)</label><mml:math id="M33"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>y</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:msubsup><mml:mrow><mml:mi>y</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:msubsup><mml:mo class="qopname">ln</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mo class="qopname">exp</mml:mo><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mo class="qopname">exp</mml:mo><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>is the modified cross-entropy loss for a single example indexed by <italic>n</italic> with:</p>
<list list-type="bullet">
<list-item><p><italic>P</italic>- the number of output channels,</p></list-item>
<list-item><p><italic>y</italic><sub><italic>p</italic></sub> - a binary indicator (0 or 1) of the desired output channel <italic>p</italic> spiking first,</p></list-item>
<list-item><p><italic>z</italic><sub><italic>p</italic></sub> - the transformed spike time of the <italic>p</italic>-th output channel <italic>z</italic>(<italic>t</italic>) &#x0003D; exp(<italic>t</italic>).</p></list-item>
</list>
<p>The parameter &#x003B3; is the synaptic regularization parameter for the spike-firing penalty</p>
<disp-formula id="E14"><label>(13)</label><mml:math id="M34"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mtext class="textrm" mathvariant="normal">spiking</mml:mtext></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>h</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>H</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>h</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:munderover></mml:mstyle><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <inline-formula><mml:math id="M35"><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo class="qopname">max</mml:mo><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:munder class="msub"><mml:mrow><mml:mo class="qopname">&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x02208;</mml:mo><mml:msub><mml:mrow><mml:mi>B</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:munder><mml:msub><mml:mrow><mml:mi>w</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:math></inline-formula> with <italic>B</italic><sub><italic>mh</italic></sub> &#x0003D; {<italic>k</italic>:<italic>t</italic><sub><italic>kh</italic></sub> &#x0003C; &#x0221E;} being the set of valid inputs for the <italic>m</italic>-th output of the <italic>h</italic>-th postsynaptic neuron. Lastly, the parameter &#x003BB; controls the strength of the <italic>L</italic><sub>2</sub> regularization term of the network weights</p>
<disp-formula id="E15"><label>(14)</label><mml:math id="M36"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munder class="msub"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:munder></mml:mstyle><mml:msubsup><mml:mrow><mml:mi>w</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>For completeness <italic>R</italic><sub><italic>mh</italic></sub> &#x0003D; 0 if {<italic>k</italic>:<italic>t</italic><sub><italic>kh</italic></sub> &#x0003C; &#x0221E;} &#x0003D; &#x02205; . Models are trained with the RMSprop optimizer (Tieleman and Hinton, <xref ref-type="bibr" rid="B66">2012</xref>) with a learning rate of 0.001 over 2,500 iterations with a batch size of 50 examples. We set &#x003BB; &#x0003D; 10<sup>&#x02212;5</sup> in all experiments. Similarly to Pabian et al. (<xref ref-type="bibr" rid="B51">2024</xref>), the regularization parameter &#x003B3; is initially set to a large value of 10<sup>5</sup> in order to guide the model toward a solution that is capable of propagating event throughout the entire network. Then, after &#x003B7; iterations, the value of &#x003B3; is decreased to 10<sup>3</sup> in order to increase the relative importance of solving the classification task.</p>
<p>We use Hyperopt (Bergstra et al., <xref ref-type="bibr" rid="B5">2013</xref>) to optimize the hyperparameters that control how the information is processed by the network. Specifically, the following Bayesian optimization search space was defined:</p>
<list list-type="bullet">
<list-item><p>the refractory period: log(&#x003C4;<sub>ref</sub>)&#x0007E;<bold>U</bold>(&#x02212;0.6, 0),</p></list-item>
<list-item><p>the longest neuron response delay &#x003C4;<sub>delay</sub> for a linearly-spaced grid: &#x003C4;<sub>max delay</sub>&#x0007E;<bold>N</bold>(<italic>t</italic><sub>avg</sub>, 1) subject to &#x003C4;<sub>max delay</sub>&#x02265;0, where <bold>N</bold>(&#x003BC;, &#x003C3;<sup>2</sup>) denotes a Gaussian random variable with mean &#x003BC; and variance &#x003C3;<sup>2</sup>, and <italic>t</italic><sub>avg</sub> is the average event time of input sequences in the current training dataset,</p></list-item>
<list-item><p>the number of iterations to train with a larger spike-firing penalty: &#x003B7;&#x0007E;<bold>U</bold>{300, 800}.</p></list-item>
</list>
<p>This parameter selection process is executed over 30 optimization iterations on a single data split from the V-fold cross-validation. Subsequently, the five best-performing parameter sets are selected and used to train MIMO SNN models on the remaining data splits. In case of tied performance scores, the top-5 Hyperopt models are chosen according to the following heuristic rule: prioritize lower &#x003C4;<sub>max delay</sub>, higher &#x003C4;<sub>ref</sub>, and lower &#x003B7;, in that order of importance. This selection strategy favors models that generate fewer events and exhibit faster response time.</p>
</sec>
</sec>
<sec sec-type="results" id="s4">
<title>4 Results</title>
<sec>
<title>4.1 Experimental setup</title>
<p>For this study, we used the same VMP dataset as Marsza&#x00142;ek et al. (<xref ref-type="bibr" rid="B31">2023</xref>). The dataset consists of 3,328 records categorized into six classes as follows: motorcycle (24.32%), bicycle (20.42%), electric scooter (19.52%), car (15.02%), delivery van (14.11%), and truck (6.61%). Each record contains simultaneous multi-frequency measurements of the real and imaginary components of the VMP signal as described in Section 2.3. For our analysis, we select a subset of channels, retaining only four of them for further processing (the R-VMP signals recorded at the lowest excitation frequency). This choice streamlines our analysis while also encouraging future research that explores the full scope of the available data.</p>
<p>In our analysis of the vehicle classification dataset we consider three encoding types introduced in Section 2.1: level-crossing, send-on-delta and leaky integrate-and-fire (LIF). For the proposed encoding parameter selection with Bayesian optimization we constrain the search space to the following set of probability distributions for each encoding type:</p>
<list list-type="bullet">
<list-item><p>number of uniformly-distributed amplitude levels of the level-crossing encoding: <italic>L</italic>&#x0007E;<bold>U</bold>{4, 5, &#x02026;, 16},</p></list-item>
<list-item><p>send-on-delta encoding threshold: &#x00394;&#x0007E;<bold>U</bold>(0.02, 0.2),</p></list-item>
<list-item><p>LIF encoding parameters: &#x003C4;<sub>int</sub>&#x0007E;<bold>U</bold>(0.05, 0.4); &#x003C4;<sub>leak</sub>&#x0007E;<bold>U</bold>(0.2, 0.6); <italic>V</italic><sub>thr</sub>&#x0007E;<bold>U</bold>(0.2, 0.4).</p></list-item>
</list>
<p>These constraints were chosen based on preliminary experiments, which showed that certain parameter values resulted in suboptimal event stream densities (e.g., the send-on-delta encoding with &#x00394;&#x0003E;0.2 produced event streams that were too sparse).</p>
<p>Before the signal-to-event encoding, we upsample the digital VMP signals from 1 kHz to 10 kHz in order to more accurately assign the event occurrence time instants. Furthermore, as shown in <xref ref-type="fig" rid="F2">Figure 2</xref>, the individual signals exhibit either all-positive or all-negative amplitudes with only minor variations near the zero baseline. Hence, we take the absolute value of each signal and apply min-max normalization to scale them to the 0&#x02013;1 range. This normalization eliminates the need to account for variability in amplitude ranges across different VMP sensors. As a side effect of this signal normalization process, the LIF-negative event type cannot occur in the encoded sequences. Additionally, the spike sequences are time-shifted so that the first event within each sequence (across all event types) occurs at a relative time <italic>t</italic> &#x0003D; 0. This adjustment serves two purposes. Firstly, the sensors are positioned in a series, making it essential to preserve the relative time shifts between channels as the vehicle moves over the measurement system. Secondly, any sensor data recorded before the vehicle enters the measurement space should be discarded, enforcing a consistent effective signal start time <italic>t</italic> &#x0003D; 0.</p>
<p><xref ref-type="fig" rid="F5">Figure 5</xref> shows an example of event sequences obtained for a VMP time series. Each signal-to-event encoding scheme produces a different number of distinct event types for each VMP sensor: level-crossing encoding produces one event type per amplitude threshold; send-on-delta encoding generates two event types (one for &#x0201C;rising&#x0201D; and one for &#x0201C;falling&#x0201D; signal amplitude); LIF encoding produces one event type due to the amplitude normalization step described earlier.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>An example of a VMP signal and its different event domain representations. Curves and raster plot points are color-coded to signify the correspondence between the input VMP channel and the output event sequences. The vertical axis for the three encoding raster plots represents different event types: 48 for level-crossing (4 channels &#x000D7; 12 amplitude levels), 8 for send-on-delta, and 4 for the LIF encoding. The number of thresholds in the level-crossing encoding was chosen arbitrarily for illustrative purposes.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-19-1610766-g0005.tif"/>
</fig>
</sec>
<sec>
<title>4.2 Encoding parameter selection for the VMP dataset</title>
<p>The results of the Hyperopt optimization runs are summarized in <xref ref-type="fig" rid="F6">Figure 6</xref>, and show a wide range of the obtained classifier accuracy spanning approximately from 0.775 to 0.950. This variation highlights the significant impact of encoding parameter selection on the k-NN classifier performance. It is evident that some data splits are more challenging to classify than others, although none of the classifiers deviates significantly from their respective overall-average performance. Furthermore, the results indicate that suboptimal LIF encoding parameters can lead to significantly worse performance than compared to the level-crossing and send-on-delta encoding schemes. In particular, the van Rossum distance parameters were also optimized, influencing the final classification scores. However, determining whether one encoding scheme is definitely superior to the others is not possible based solely on these raw results, due to the factors discussed earlier in Section 3.1.3.</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p>Classification accuracy of k-NN models trained during the stratified 10-fold cross-validated Bayesian optimization procedure. The dashed lines denote the average score over all data splits.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-19-1610766-g0006.tif"/>
</fig>
<p>Additional insights into the encoding schemes were gained by analyzing the classification scores of models trained with Hyperopt in relation to the average number of events produced per VMP sensor<xref ref-type="fn" rid="fn0003"><sup>3</sup></xref>. The total number of events in the spike train was normalized by the number of VMP sensors to account for the fact that signals observed by the sensors tend to generate roughly the same number of events after encoding. As a result, the total number of events generated in the system scales linearly with the number of channels, regardless of the selected coding scheme. Additionally, this normalization allows for an assessment of the computational burden associated with adding an additional VMP sensor to the system. The results presented in <xref ref-type="fig" rid="F7">Figure 7</xref> indicate that the number of events alone is not the sole determinant of classification performance. Specifically, send-on-delta models maintained similar performance regardless of the number of events produced by the encoding scheme. In contrast, different LIF parameter sets, even when generating event sequences of similar length, led to significantly varied k-NN classifier performance.</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Kernel-smoothed density estimators of the k-NN classifier accuracy scores vs. the average number of events produced by a given encoding scheme per the VMP sensor. Heat maps representing the estimated density were min-max normalized separately.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-19-1610766-g0007.tif"/>
</fig>
<p>The selected parameter sets, determined using this strategy, are presented in <xref ref-type="table" rid="T1">Table 1</xref>, which also includes the stratified 10-fold cross-validated k-NN classifier performance for data splits preprocessed with the chosen encoding schemes. According to these results, level-crossing encoding performs slightly better than the alternatives for the VMP vehicle type classification task. Each signal prior to event encoding contains an average of 1,590 samples. This means that the most data-intensive encoding scheme&#x02014;send-on-delta&#x02014;produces a representation that uses approximately 97.8% fewer samples than the original signal.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Signal-to-event encoding parameters chosen according to the weighted median selection strategy.</p></caption>
<table frame="box" rules="all">
<thead>
<tr style="background-color:#8f9496;color:#ffffff">
<th valign="top" align="left"><bold>Encoding type</bold></th>
<th valign="top" align="center"><bold>Chosen parameters</bold></th>
<th valign="top" align="center"><bold>Number of events per VMP sensor</bold></th>
<th valign="top" align="center" colspan="2"><bold>Stratified 10-fold cross-validated k-NN performance</bold></th>
</tr>
</thead>
<tbody>
<tr style="background-color:#8f9496;color:#ffffff">
<td/>
<td/>
<td/>
<td valign="top" align="center"><bold>Accuracy</bold></td>
<td valign="top" align="center"><bold><italic>F</italic>1-score</bold></td>
</tr> <tr>
<td valign="top" align="left">level-crossing</td>
<td valign="top" align="center"><italic>L</italic> &#x0003D; 12</td>
<td valign="top" align="center">29.484 &#x000B1; 8.996</td>
<td valign="top" align="center">0.912 &#x000B1; 0.011</td>
<td valign="top" align="center">0.907 &#x000B1; 0.016</td>
</tr> <tr>
<td valign="top" align="left">send-on-delta</td>
<td valign="top" align="center">&#x00394; &#x0003D; 0.06</td>
<td valign="top" align="center">36.536 &#x000B1; 9.040</td>
<td valign="top" align="center">0.910 &#x000B1; 0.014</td>
<td valign="top" align="center">0.909 &#x000B1; 0.019</td>
</tr> <tr>
<td valign="top" align="left">LIF</td>
<td valign="top" align="center">&#x003C4;<sub>int</sub> &#x0003D; 0.1</td>
<td valign="top" align="center">8.317 &#x000B1; 5.673</td>
<td valign="top" align="center">0.905 &#x000B1; 0.011</td>
<td valign="top" align="center">0.900 &#x000B1; 0.015</td>
</tr>
 <tr>
<td/>
<td valign="top" align="center">&#x003C4;<sub>leak</sub> &#x0003D; 0.5</td>
<td/>
<td/>
<td/>
</tr>
 <tr>
<td/>
<td valign="top" align="center"><italic>V</italic><sub>thr</sub> &#x0003D; 0.2</td>
<td/>
<td/>
<td/>
</tr> <tr>
<td valign="top" align="left">LIF (<italic>post-hoc</italic> analysis)</td>
<td valign="top" align="center">&#x003C4;<sub>int</sub> &#x0003D; 0.1</td>
<td valign="top" align="center">20.655 &#x000B1; 13.768</td>
<td valign="top" align="center">0.910 &#x000B1; 0.010</td>
<td valign="top" align="center">0.906 &#x000B1; 0.015</td>
</tr>
 <tr>
<td/>
<td valign="top" align="center">&#x003C4;<sub>leak</sub> &#x0003D; 0.35</td>
<td/>
<td/>
<td/>
</tr>
 <tr>
<td/>
<td valign="top" align="center"><italic>V</italic><sub>thr</sub> &#x0003D; 0.1</td>
<td/>
<td/>
<td/>
</tr></tbody>
</table>
<table-wrap-foot>
<p>The multi-neuron van Rossum distance parameters used to evaluate the k-NN classifier are omitted. The last row shows the results of a <italic>post-hoc</italic> experiment that repeated the LIF parameter selection analysis with a wider sampling space for the <italic>V</italic><sub>thr</sub> parameter and a higher number of the Bayesian optimization iterations. The original result for the LIF encoding parameters was used in subsequent analysis of the SNN classifier training.</p>
</table-wrap-foot>
</table-wrap>
<p>Next, we analyzed the average number of events produced by the selected encoding schemes for each VMP sensor separately. Together with the classification results, this analysis provides insights into the encoding efficiency. <xref ref-type="fig" rid="F8">Figure 8</xref> summarizes our findings. To visualize these distributions, we use letter-value plots (Hofmann et al., <xref ref-type="bibr" rid="B16">2017</xref>)&#x02014;an extension of the classical boxplot that more accurately represents distribution tails for large datasets. The results reveal a clear distinction between the empirical distributions of event sequences produced by the two standard inductive loops (IL1 and IL3) and the two slim loops (IL2 and IL4), regardless of the signal-to-event encoding scheme. Additionally, the level-crossing and send-on-delta encoding schemes produce similar overall number of events, with send-on-delta sequences containing slightly more events on average. In contrast, the LIF encoding generates significantly fewer events with the median number of events being approximately three times lower for IL1&#x02013;IL3 sensors and about six times lower for IL2&#x02013;IL4 sensors. Moreover, the LIF encoding is capable of adapting to the original signal variability, as evidenced by the broad range of event counts across all data samples. Given that all three classifiers in <xref ref-type="table" rid="T1">Table 1</xref> achieve comparable classification accuracy, the LIF encoding appears to be more efficient in terms of the amount of information encoded per event.</p>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p>Letter-value plots of the total number of events produced by signal-to-event encoding schemes selected according to the weighted median strategy for each input VMP sensor.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-19-1610766-g0008.tif"/>
</fig>
<p>Finally, relating the set of selected encoding parameters summarized in <xref ref-type="table" rid="T1">Table 1</xref> with the original optimization search space reveals that the proposed method did not trivially select the parameter sets that would maximize the number of events in the event stream sequences. In fact, only the threshold voltage of the LIF encoding was on the sampling space boundary [selected <italic>V</italic><sub>thr</sub> &#x0003D; 0.2 for <italic>V</italic><sub>thr</sub>&#x0007E;<bold>U</bold>(0.2, 0.4)] with all other parameters being selected away from the boundary. This result suggests that the proposed method is not biased in favor of high event density encodings. However, the method itself does not have any safeguards and relies on carefully choosing the parameter bounds to avoid trivial solutions.</p>
<p>Selecting a parameter value for <italic>V</italic><sub>thr</sub> that lies on the search space boundary raises a question of whether the parameter space for the LIF encoding has been fully explored. In order to investigate this topic and challenge our initial choice of the search space, we repeated the LIF encoding parameter selection in a <italic>post-hoc</italic> experiment. In this analysis the LIF threshold voltage parameter sampling space was specified as <italic>V</italic><sub>thr</sub>&#x0007E;<bold>U</bold>(0.05, 0.4). Furthermore, due to an increase in the search space scope, we increased the number of iterations for the Bayesian optimization from 150 to 250. The new set of LIF encoding parameters selected by our method, as well as the average number of generated events and the classifier performance metrics, is presented in <xref ref-type="table" rid="T1">Table 1</xref>. Not only is the new choice of <italic>V</italic><sub>thr</sub> away from the sampling space boundary, but also the k-NN classifier performance for this encoding is on par with the one achieved by the send-on-delta scheme. Additionally, compared to the original result, this new set of LIF encoding parameters more than doubles the number of generated events (on average). Overall, this shows that the proposed method will achieve suboptimal results when the choice of the parameter sampling space for the Bayesian optimization is too constrained. As this result was obtained in a <italic>post-hoc</italic> experiment, we did not use this updated set of parameters of the LIF encoding in subsequent analysis.</p>
</sec>
<sec>
<title>4.3 Vehicle classification with the MIMO SNN</title>
<p>The set of parameters identified in the previous analysis was used to encode VMP signals into the spiking domain before training the MIMO SNN models for the vehicle type identification task. We applied the same stratified 10-fold data split for cross-validation. The network parameter settings and the hyperparameter optimization grid were described in Section 3.2. In all experiments we used the same <italic>C</italic>-64-128-128-6 architecture, where the number of input neurons <italic>C</italic> depends entirely on the encoding scheme. The specific values of <italic>C</italic> for the encoded VMP data are as follows: <italic>C</italic> &#x0003D; 48 for level-crossing encoding, <italic>C</italic> &#x0003D; 8 for send-on-delta encoding, and <italic>C</italic> &#x0003D; 4 for LIF encoding. This results in a small network with a total of 25, 344&#x0002B;64<italic>C</italic> parameters. However, it is important to note that the effective capacity of the MIMO SNN model is higher than what the raw parameter count suggests, due to the repeated firing of the IF neurons, which are influenced by the refractory period &#x003C4;<sub>ref</sub>. In total, we optimized the weights of 225 models across the three spike encoding schemes.</p>
<p><xref ref-type="fig" rid="F9">Figure 9a</xref> summarizes the classification accuracy of models trained during the Bayesian optimization procedure on a single data split from the stratified 10-fold cross-validation. On average, the LIF encoding performed worse than the alternatives. Recall that all three encoding schemes when parameterized using the weighted median selection strategy achieved similar k-NN classifier accuracy of approximately 0.908 (<xref ref-type="table" rid="T1">Table 1</xref>). This indicates that the Bayesian optimization procedure failed to find a LIF model that outperformed its own baseline. In contrast, the level-crossing and send-on-delta models performed better than their respective k-NN classifiers, with both achieving similar median accuracy. However, the presence of outliers in the lower accuracy range suggests that poor hyperparameter selection during training can significantly degrade the final model performance. Nonetheless, these results confirm that running the optimization procedure over multiple iterations can lead to a high-performing model, provided that the signal-to-event encoding is chosen correctly.</p>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption><p><bold>(a)</bold> Summary of the SNN model performance for models trained with Hyperopt (30 optimization iterations on a single data split of the stratified 10-fold cross-validation). <bold>(b)</bold> Classification accuracy of the SNN models trained with top-5 parameter sets found by the Bayesian optimization process on data split &#x00023;1. The dashed lines denote the average score over all data splits.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-19-1610766-g0009.tif"/>
</fig>
<p><xref ref-type="fig" rid="F9">Figure 9b</xref> presents the classification accuracy for models trained with the top-5 best-performing parameter sets selected by Hyperopt, separately for each encoding type and each cross-validation data split. This analysis evaluated the robustness of classifiers trained with these hyperparameters when applied to different input data. Several differences emerge between models operating on differently encoded data. Once again, LIF models exhibited significantly lower absolute performance compared to the other encoding schemes. Additionally, both level-crossing and LIF models evaluated on data splits &#x00023;2&#x02013;10 demonstrated worse classification accuracy than the reference split &#x00023;1, which was used to determine the top-5 hyperparameter sets. In contrast, send-on-delta model evaluated on data splits &#x00023;3, &#x00023;7, &#x00023;8, and &#x00023;10 achieved performance levels comparable to scenario &#x00023;1. This suggests that the hyperparameters selected for send-on-delta encoding are more robust across different training datasets. However, this encoding scheme also exhibited an extreme performance outlier in data split &#x00023;9, deviating from all other trained models. Overall, these findings indicate that send-on-delta encoding is less sensitive to variations in training data, whereas level-crossing and the LIF encoding require more precise tuning to achieve optimal performance.</p>
<p>Lastly, <xref ref-type="table" rid="T2">Table 2</xref> presents the top-5 parameter sets for the three encoding types along with the stratified 10-fold cross-validated performance metrics for models trained with these hyperparameters. Notably, the top-5 parameter sets selected for the send-of-delta and LIF models are significantly more similar to each other than those chosen for the level-crossing encoding. This may suggest that, for this specific task, 30 steps of the Bayesian optimization process are sufficient to converge to a good-enough (local minimum) solution for these two model types, but not enough for the level-crossing scheme. Interestingly, the selected values of &#x003C4;<sub>max delay</sub> were smaller for the send-on-delta and LIF models when compared to level-crossing scheme. This suggests that a nonzero &#x003C4;<sub>max delay</sub> is more beneficial when the number of neurons in the input layer is large (as in level-crossing encoder) rather than when they generate a large number of events (as in send-on-delta encoder). Finally, recall that the prior on the number of training iterations with a larger spike-firing penalty &#x003B7; was <inline-formula><mml:math id="M37"><mml:mrow><mml:mi mathvariant="script">U</mml:mi></mml:mrow><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mn>300</mml:mn><mml:mo>,</mml:mo><mml:mn>800</mml:mn></mml:mrow><mml:mo>}</mml:mo></mml:mrow></mml:math></inline-formula>. Almost all of the selected values of &#x003B7; ended up at the higher end of this range. This indicates that training the model for longer with a relatively smaller weight assigned to the task-specific loss leads to better final model performance.</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Summary of the top-5 parameter sets found by the Bayesian optimization process used to train the SNN models for the respective spike encoding type.</p></caption>
<table frame="box" rules="all">
<thead>
<tr style="background-color:#8f9496;color:#ffffff">
<th valign="top" align="left"><bold>Encoding type</bold></th>
<th valign="top" align="center"><bold>Number of events per VMP sensor</bold></th>
<th valign="top" align="center" colspan="3"><bold>Top-5 parameter sets</bold></th>
<th valign="top" align="center" colspan="2"><bold>Stratified 10-fold cross-validated SNN performance</bold></th>
</tr>
</thead>
<tbody>
<tr style="background-color:#8f9496;color:#ffffff">
<td/>
<td/>
<td valign="top" align="center">&#x003C4;<sub>max delay</sub></td>
<td valign="top" align="center">&#x003C4;<sub>ref</sub></td>
<td valign="top" align="center">&#x003B7;</td>
<td valign="top" align="center"><bold>Accuracy</bold></td>
<td valign="top" align="center"><bold><italic>F</italic>1-score</bold></td>
</tr> <tr>
<td valign="top" align="left">Level-crossing</td>
<td valign="top" align="center">29.484 &#x000B1; 8.996</td>
<td valign="top" align="center">0.9</td>
<td valign="top" align="center">1.0000</td>
<td valign="top" align="center">700</td>
<td valign="top" align="center">0.926 &#x000B1; 0.019</td>
<td valign="top" align="center">0.929 &#x000B1; 0.018</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">1.2</td>
<td valign="top" align="center">0.2512</td>
<td valign="top" align="center">400</td>
<td valign="top" align="center">0.917 &#x000B1; 0.024</td>
<td valign="top" align="center">0.920 &#x000B1; 0.027</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.9</td>
<td valign="top" align="center">0.6310</td>
<td valign="top" align="center">600</td>
<td valign="top" align="center">0.915 &#x000B1; 0.035</td>
<td valign="top" align="center">0.917 &#x000B1; 0.033</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.9</td>
<td valign="top" align="center">0.2512</td>
<td valign="top" align="center">700</td>
<td valign="top" align="center">0.919 &#x000B1; 0.024</td>
<td valign="top" align="center">0.922 &#x000B1; 0.024</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.7</td>
<td valign="top" align="center">0.7943</td>
<td valign="top" align="center">500</td>
<td valign="top" align="center">0.913 &#x000B1; 0.028</td>
<td valign="top" align="center">0.919 &#x000B1; 0.028</td>
</tr> <tr>
<td valign="top" align="left">Send-on-delta</td>
<td valign="top" align="center">36.536 &#x000B1; 9.040</td>
<td valign="top" align="center">0.3</td>
<td valign="top" align="center">0.3981</td>
<td valign="top" align="center">600</td>
<td valign="top" align="center">0.937 &#x000B1; 0.030</td>
<td valign="top" align="center">0.928 &#x000B1; 0.031</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.4</td>
<td valign="top" align="center">0.5012</td>
<td valign="top" align="center">600</td>
<td valign="top" align="center">0.926 &#x000B1; 0.036</td>
<td valign="top" align="center">0.913 &#x000B1; 0.044</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.4</td>
<td valign="top" align="center">0.5012</td>
<td valign="top" align="center">700</td>
<td valign="top" align="center"><bold>0.946 &#x000B1; 0.018</bold></td>
<td valign="top" align="center"><bold>0.937 &#x000B1; 0.021</bold></td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.4</td>
<td valign="top" align="center">0.5012</td>
<td valign="top" align="center">800</td>
<td valign="top" align="center">0.937 &#x000B1; 0.026</td>
<td valign="top" align="center">0.926 &#x000B1; 0.027</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.0</td>
<td valign="top" align="center">0.3162</td>
<td valign="top" align="center">500</td>
<td valign="top" align="center">0.909 &#x000B1; 0.073</td>
<td valign="top" align="center">0.907 &#x000B1; 0.068</td>
</tr> <tr>
<td valign="top" align="left">LIF</td>
<td valign="top" align="center">8.317 &#x000B1; 5.673</td>
<td valign="top" align="center">0.0</td>
<td valign="top" align="center">0.5012</td>
<td valign="top" align="center">300</td>
<td valign="top" align="center">0.861 &#x000B1; 0.024</td>
<td valign="top" align="center">0.859 &#x000B1; 0.027</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.2</td>
<td valign="top" align="center">0.5012</td>
<td valign="top" align="center">600</td>
<td valign="top" align="center">0.868 &#x000B1; 0.031</td>
<td valign="top" align="center">0.856 &#x000B1; 0.040</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.0</td>
<td valign="top" align="center">0.3981</td>
<td valign="top" align="center">700</td>
<td valign="top" align="center">0.854 &#x000B1; 0.032</td>
<td valign="top" align="center">0.852 &#x000B1; 0.038</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.0</td>
<td valign="top" align="center">0.5012</td>
<td valign="top" align="center">400</td>
<td valign="top" align="center">0.859 &#x000B1; 0.024</td>
<td valign="top" align="center">0.858 &#x000B1; 0.024</td>
</tr>
 <tr>
<td/>
<td/>
<td valign="top" align="center">0.0</td>
<td valign="top" align="center">0.5012</td>
<td valign="top" align="center">700</td>
<td valign="top" align="center">0.855 &#x000B1; 0.031</td>
<td valign="top" align="center">0.851 &#x000B1; 0.038</td>
</tr></tbody>
</table>
<table-wrap-foot>
<p>Bold values indicate the best classifier accuracy and F1-score across all experimental settings.</p>
</table-wrap-foot>
</table-wrap>
<p>When compared to the performance of the baseline k-NN classifiers (<xref ref-type="table" rid="T1">Table 1</xref>), training the SNN results in lower error rates for the vehicle type classification task, provided that the signal-to-event encoding scheme is properly selected. While processing an example through the SNN is undoubtedly more computationally demanding than computing the multi-neuron van Rossum distance between a pair of examples, a k-NN-based system prediction time scales poorly with the size of the reference database, making it less efficient as the number of training examples increases.</p>
<p>Considering these factors, the send-on-delta model emerged as the most effective event-based encoding scheme for the VMP-based vehicle classification. This model achieved the highest stratified 10-fold cross-validation performance, demonstrated greater robustness across different data splits, and allowed the Bayesian optimization procedure to converge to a locally optimal solution within the given number of optimization steps. In particular, this cannot be solely attributed to the large average number of events generated by the encoding, as it is comparable to that for the level-crossing event sequences (as shown in <xref ref-type="fig" rid="F8">Figure 8</xref>).</p>
</sec>
</sec>
<sec id="s5">
<title>5 Discussions and conclusions</title>
<p>We proposed a novel methodology for selecting signal-to-event encoding parameters that preserves critical signal information for multiple event classification tasks. This approach is model-agnostic and independent of the complexity of the classification model that will eventually process the data. By decoupling encoding parameter selection from the hyperparameter tuning phase of model development one can simplify the overall prototyping process. Additionally, this method enables a clearer assessment of how different input encoding parameters influence final model performance, ensuring a more consistent and meaningful comparison across models.</p>
<p>This study evaluated three event-based signal encoding schemes in a vehicle classification task: level-crossing, send-on-delta and leaky integrate-and-fire (LIF) models. The aim of the analysis was to assess the impact of encoding schemes on machine learning performance. Our findings reveal that the accuracy of the k-NN classifier varies significantly, ranging from approximately 0.775 to 0.950. This wide performance gap underscores the critical importance of selecting appropriate encoding parameters. To address this issue, we introduced a <italic>weighted median selection strategy</italic>, which constructs performance ranking lists for each data split and determines the median performance across different hyperparameter settings. This approach enabled us to establish a single optimal parameter set for each encoding scheme, yielding average classification accuracies ranging from 0.905 (for the LIF encoding) to 0.912 (for send-on-delta). These results provide a robust baseline for SNN models, given the comparable cross-validated performance of the three k-NN classifiers. Moreover, the selected encoding parameter sets significantly reduce data redundancy compared to the original digitized signal. The most data-intensive encoding&#x02014;send-on-delta&#x02014;still produces approximately 97.8% fewer samples than the original signal representation, demonstrating the efficiency of event-based encoding. All three encoding schemes achieved comparable levels of classification accuracy while producing encodings of different event density. This highlights the importance of testing multiple alternative encoding schemes in order to figure out the one that best satisfies the system design constraints.</p>
<p>For the SNN trained on data encoded by the event-based schemes, the send-on-delta models emerged as the best-performing group of networks. They achieved not only the highest overall cross-validated accuracy (0.946 on average), but also demonstrated greater robustness across different data splits. Additionally, the Bayesian optimization procedure successfully converged to a locally optimal solution within the limited number of optimization steps. However, this was achieved while generating the highest average number of events among the tested signal-to-event encoding schemes, slightly more events than the second-best approach, level-crossing encoding. In contrast, the SNNs trained on time-series data encoded using the LIF scheme performed significantly worse than the tested alternatives. We hypothesize that three key factors contributed to this outcome: the small number of distinct event types produced by the encoding (only one per input sensor), the relatively low number of events observed by the network (approximately three to six times fewer than in the alternative schemes), and the complexity of optimizing multiple encoding parameters (compared to a single parameter for send-on-delta and level-crossing schemes). Further investigation is needed to determine whether these factors universally impair the SNN classifier performance, regardless of the encoding type. Note that in a <italic>post-hoc</italic> experiment for the LIF parameter selection a slightly different set of parameters was selected, one that more than doubles the number of generated events (on average). While using this parameter set to train an SNN could have addressed the problem of low event stream density, it would not have helped with the small number of distinct event types produced by the LIF encoding. Nevertheless, the SNN training experiment should be considered as supplementary to the main goal of the paper (choosing an event-based signal encoding parameters) and is not intended to be a commentary on which encoding scheme&#x02014;if any&#x02014;is the best. Overall, our findings suggest the presence of a trade-off between classification performance and energy efficiency, measured by the number of generated events. These considerations must be carefully balanced when designing an end-to-end event-centric machine learning solution.</p>
<p>In this paper, we employed a simple multichannel spike train model as defined in <xref ref-type="disp-formula" rid="E7">Equation 6</xref>. This model assumes that there is no functional dependence between events from different channels. However, in certain scenarios, this assumption might not hold, as events from multiple channels can interact. To account for such dependencies, an alternative spike train model can be defined as follows:</p>
<disp-formula id="E16"><mml:math id="M38"><mml:mtable class="eqnarray" columnalign="center"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:mrow><mml:mi mathvariant="script">U</mml:mi></mml:mrow></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mo>=</mml:mo></mml:mtd><mml:mtd><mml:mstyle displaystyle="true"><mml:munder class="msub"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:munder></mml:mstyle><mml:mi>h</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:mi>&#x003B1;</mml:mi><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mtable style="text-align:axis;" equalrows="false" columnlines="none none none none none none none none none" equalcolumns="false" class="array"><mml:mtr><mml:mtd><mml:mi>s</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>s</mml:mi><mml:mo>&#x02260;</mml:mo><mml:mi>p</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mrow><mml:mi>P</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:mstyle displaystyle="true"><mml:munder class="msub"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:munder></mml:mstyle><mml:mi>h</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:msubsup><mml:mrow><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E17"><label>(15)</label><mml:math id="M39"><mml:mtable class="eqnarray" columnalign="center"><mml:mtr><mml:mtd><mml:mi>p</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:mi>P</mml:mi><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where the parameter &#x003B1;&#x0003E;0 controls the strength of the multichannel interactions. If &#x003B1; &#x0003D; 0, the model reduces to the one defined in <xref ref-type="disp-formula" rid="E7">Equation 6</xref>. Extending our learning strategy to this dependent multichannel spike train model presents an interesting option for further research. Additionally, this study could be extended to incorporate other event-based encoding schemes and to train more complex classifiers using the proposed methodology. Finally, our optimization objective considered only the classification accuracy and we made <italic>post hoc</italic> insights into the encoding event density efficiency for parameter sets selected according to such objective. It would be interesting to see a multi-objective optimization approach that also considers this event-based efficiency to allow better control over the trade-off between classification performance and the number of generated samples.</p></sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The data analyzed in this study is subject to the following licenses/restrictions. Data provided by a third-party&#x02014;the authors of this article do not hold the right to share the data. Requests to access these datasets should be directed to Zbigniew Marsza&#x00142;ek, <email>antic&#x00040;agh.edu.pl</email>.</p>
</sec>
<sec sec-type="author-contributions" id="s7">
<title>Author contributions</title>
<p>MPab: Conceptualization, Investigation, Methodology, Software, Writing &#x02013; original draft. DR: Conceptualization, Methodology, Writing &#x02013; original draft. MPaw: Supervision, Writing &#x02013; review &#x00026; editing. MM: Project administration, Writing &#x02013; review &#x00026; editing. RS: Project administration, Writing &#x02013; review &#x00026; editing.</p>
</sec>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the NSERC Grant 327447, the Polish National Center of Science under Grant DEC-2017/27/B/ST7/03082, and by the Ministry of Education and Science under the program called &#x0201C;Preparation of metrological security for the introduction of direct mass enforcement (e-WIM) WIM systems in Poland&#x0201D; (project number PM/SP/0041/2021/1).</p>
</sec>
<ack><p>We would like to thank Dr. Zbigniew Marsza&#x00142;ek for providing us with the sensor data to validate our methodology.</p>
</ack>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s9">
<title>Generative AI statement</title>
<p>The author(s) declare that no Gen AI was used in the creation of this manuscript.</p></sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x00027;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<fn-group>
<fn id="fn0001"><p><sup>1</sup>The van Rossum distance was originally defined for neural spike responses, therefore its formulation refers to various (populations of) neurons. We adopt this terminology for consistency with previous works in the field, noting that different &#x0201C;neurons&#x0201D; relate to &#x0201C;event types&#x0201D;, whereas &#x0201C;populations of neurons&#x0201D; intuitively correspond to &#x0201C;sets of event types.&#x0201D;</p></fn>
<fn id="fn0002"><p><sup>2</sup>In Bayesian optimization, the sampling process depends on the sequence of scores obtained in previous iterations. Since each data split contains a different subset of examples, it is not possible to guarantee that the same parameter sets will be sampled across splits even if the sampling is pseudo-random, i.e., initialized with a fixed random seed, unless the sequence of scores in successive iterations is identical.</p></fn>
<fn id="fn0003"><p><sup>3</sup>Recall that the term &#x0201C;VMP sensor&#x0201D; can refer to either the physical IL sensor or the excitation frequency of the loop. Throughout our analyses we used four &#x0201C;VMP sensors&#x0201D; representing the R-VMP signals registered by the quad-loop system at a single excitation frequency per loop.</p></fn>
</fn-group>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Abrahams</surname> <given-names>J.</given-names></name></person-group> (<year>1986</year>). <article-title>&#x0201C;A survey of recent progress on level-crossing problems for random processes,&#x0201D;</article-title> in <source>Communications and Networks: A Survey of Recent Advances</source> (<publisher-loc>New York, NY</publisher-loc>: <publisher-name>ACM</publisher-name>), <fpage>6</fpage>&#x02013;<lpage>25</lpage>. <pub-id pub-id-type="doi">10.1007/978-1-4612-4904-7_2</pub-id></citation>
</ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aranda-Escol&#x000E1;stico</surname> <given-names>E.</given-names></name> <name><surname>Guinaldo</surname> <given-names>M.</given-names></name> <name><surname>Mi&#x0015B;kowicz</surname> <given-names>M.</given-names></name> <name><surname>Dormido</surname> <given-names>S.</given-names></name></person-group> (<year>2024</year>). <article-title>Event-based control in industry practice: paving the way toward resource-efficient industrial internet of things</article-title>. <source>IEEE Ind. Electron. Mag</source>. <volume>18</volume>, <fpage>38</fpage>&#x02013;<lpage>47</lpage>. <pub-id pub-id-type="doi">10.1109/MIE.2023.3286984</pub-id></citation>
</ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Auge</surname> <given-names>D.</given-names></name> <name><surname>Hille</surname> <given-names>J.</given-names></name> <name><surname>Mueller</surname> <given-names>E.</given-names></name> <name><surname>Knoll</surname> <given-names>A.</given-names></name></person-group> (<year>2021</year>). <article-title>A survey of encoding techniques for signal processing in spiking neural networks</article-title>. <source>Neural Process. Lett</source>. <volume>53</volume>, <fpage>4693</fpage>&#x02013;<lpage>4710</lpage>. <pub-id pub-id-type="doi">10.1007/s11063-021-10562-2</pub-id></citation>
</ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bauer</surname> <given-names>F. C.</given-names></name> <name><surname>Lenz</surname> <given-names>G.</given-names></name> <name><surname>Haghighatshoar</surname> <given-names>S.</given-names></name> <name><surname>Sheik</surname> <given-names>S.</given-names></name></person-group> (<year>2023</year>). <article-title>EXODUS: stable and efficient training of spiking neural networks</article-title>. <source>Front. Neurosci</source>. <volume>17</volume>:<fpage>1110444</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2023.1110444</pub-id><pub-id pub-id-type="pmid">36845419</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Bergstra</surname> <given-names>J.</given-names></name> <name><surname>Yamins</surname> <given-names>D.</given-names></name> <name><surname>Cox</surname> <given-names>D.</given-names></name></person-group> (<year>2013</year>). <article-title>&#x0201C;Making a science of model search: hyperparameter optimization in hundreds of dimensions for vision architectures,&#x0201D;</article-title> in <source>Proceedings of the 30th International Conference on Machine Learning, volume 28 of Proceedings of Machine Learning Research</source>, eds. S. Dasgupta, and D. McAllester (<publisher-loc>Atlanta, GA</publisher-loc>: <publisher-name>PMLR</publisher-name>), <fpage>115</fpage>&#x02013;<lpage>123</lpage>.</citation>
</ref>
<ref id="B6">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Cassandras</surname> <given-names>C. G.</given-names></name> <name><surname>Lafortune</surname> <given-names>S.</given-names></name></person-group> (<year>2008</year>). <source>Introduction to Discrete Event Systems</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Springer</publisher-name>. <pub-id pub-id-type="doi">10.1007/978-0-387-68612-7</pub-id></citation>
</ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cheng</surname> <given-names>G.</given-names></name> <name><surname>Dean-Leon</surname> <given-names>E.</given-names></name> <name><surname>Bergner</surname> <given-names>F.</given-names></name> <name><surname>Olvera</surname> <given-names>J. R. G.</given-names></name> <name><surname>Leboutet</surname> <given-names>Q.</given-names></name> <name><surname>Mittendorfer</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>A comprehensive realization of robot skin: sensors, sensing, control, and applications</article-title>. <source>Proc. IEEE</source> <volume>107</volume>, <fpage>2034</fpage>&#x02013;<lpage>2051</lpage>. <pub-id pub-id-type="doi">10.1109/JPROC.2019.2933348</pub-id></citation>
</ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Coifman</surname> <given-names>B.</given-names></name> <name><surname>Dhoorjaty</surname> <given-names>S.</given-names></name> <name><surname>Lee</surname> <given-names>Z.-H.</given-names></name></person-group> (<year>2003</year>). <article-title>Estimating median velocity instead of mean velocity at single loop detectors</article-title>. <source>Transp. Res. C Emerg. Technol</source>. <volume>11</volume>, <fpage>211</fpage>&#x02013;<lpage>222</lpage>. <pub-id pub-id-type="doi">10.1016/S0968-090X(03)00025-1</pub-id></citation>
</ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Falez</surname> <given-names>P.</given-names></name> <name><surname>Tirilly</surname> <given-names>P.</given-names></name> <name><surname>Bilasco</surname> <given-names>I. M.</given-names></name> <name><surname>Devienne</surname> <given-names>P.</given-names></name> <name><surname>Boulet</surname> <given-names>P.</given-names></name></person-group> (<year>2019</year>). <article-title>Unsupervised visual feature learning with spike-timing-dependent plasticity: how far are we from traditional feature learning approaches?</article-title> <source>Pattern Recognit</source>. <volume>93</volume>, <fpage>418</fpage>&#x02013;<lpage>429</lpage>. <pub-id pub-id-type="doi">10.1016/j.patcog.2019.04.016</pub-id></citation>
</ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gajda</surname> <given-names>J.</given-names></name> <name><surname>Piwowar</surname> <given-names>P.</given-names></name> <name><surname>Sroka</surname> <given-names>R.</given-names></name> <name><surname>Stencel</surname> <given-names>M.</given-names></name> <name><surname>Zeglen</surname> <given-names>T.</given-names></name></person-group> (<year>2012</year>). <article-title>Application of inductive loops as wheel detectors</article-title>. <source>Transp. Res. C: Emerg. Technol</source>. <volume>21</volume>, <fpage>57</fpage>&#x02013;<lpage>66</lpage>. <pub-id pub-id-type="doi">10.1016/j.trc.2011.08.010</pub-id></citation>
</ref>
<ref id="B11">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Gajda</surname> <given-names>J.</given-names></name> <name><surname>Sroka</surname> <given-names>R.</given-names></name> <name><surname>Stencel</surname> <given-names>M.</given-names></name> <name><surname>Wajda</surname> <given-names>A.</given-names></name> <name><surname>Zeglen</surname> <given-names>T.</given-names></name></person-group> (<year>2001</year>). <article-title>&#x0201C;A vehicle classification based on inductive loop detectors,&#x0201D;</article-title> in <source>IMTC 2001. Proceedings of the 18th IEEE Instrumentation and Measurement Technology Conference. Rediscovering Measurement in the Age of Informatics (Cat. No.01CH 37188), Volume 1</source> (<publisher-loc>Budapest</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>460</fpage>&#x02013;<lpage>464</lpage>. <pub-id pub-id-type="doi">10.1109/IMTC.2001.928860</pub-id></citation>
</ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ge</surname> <given-names>X.</given-names></name> <name><surname>Han</surname> <given-names>Q.-L.</given-names></name> <name><surname>Zhang</surname> <given-names>X.-M.</given-names></name> <name><surname>Ding</surname> <given-names>L.</given-names></name> <name><surname>Yang</surname> <given-names>F.</given-names></name></person-group> (<year>2020</year>). <article-title>Distributed event-triggered estimation over sensor networks: a survey</article-title>. <source>IEEE Trans. Cybern</source>. <volume>50</volume>, <fpage>1306</fpage>&#x02013;<lpage>1320</lpage>. <pub-id pub-id-type="doi">10.1109/TCYB.2019.2917179</pub-id><pub-id pub-id-type="pmid">31199279</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guilbert</surname> <given-names>D.</given-names></name> <name><surname>Ieng</surname> <given-names>S.-S.</given-names></name> <name><surname>Bastard</surname> <given-names>C. L.</given-names></name> <name><surname>Wang</surname> <given-names>Y.</given-names></name></person-group> (<year>2014</year>). <article-title>Robust blind deconvolution process for vehicle reidentification by an inductive loop detector</article-title>. <source>IEEE Sens. J</source>. <volume>14</volume>, <fpage>4315</fpage>&#x02013;<lpage>4322</lpage>. <pub-id pub-id-type="doi">10.1109/JSEN.2014.2345755</pub-id></citation>
</ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guilbert</surname> <given-names>D.</given-names></name> <name><surname>Le Bastard</surname> <given-names>C.</given-names></name> <name><surname>Ieng</surname> <given-names>S.-S.</given-names></name> <name><surname>Wang</surname> <given-names>Y.</given-names></name></person-group> (<year>2013</year>). <article-title>&#x0201C;Re-identification by inductive loop detector: experimentation on target origin-destination matrix,&#x0201D;</article-title> in <source>2013 IEEE Intelligent Vehicles Symposium (IV)</source>, <fpage>1421</fpage>&#x02013;<lpage>1427</lpage>. <pub-id pub-id-type="doi">10.1109/IVS.2013.6629666</pub-id></citation>
</ref>
<ref id="B15">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Heemels</surname> <given-names>W. P.</given-names></name> <name><surname>Johansson</surname> <given-names>K. H.</given-names></name> <name><surname>Tabuada</surname> <given-names>P.</given-names></name></person-group> (<year>2012</year>). <article-title>&#x0201C;An introduction to event-triggered and self-triggered control,&#x0201D;</article-title> in <source>2012 IEEE 51st IEEE Conference on Decision and Control (CDC)</source> (<publisher-loc>Maui, HI</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>3270</fpage>&#x02013;<lpage>3285</lpage>. <pub-id pub-id-type="doi">10.1109/CDC.2012.6425820</pub-id></citation>
</ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hofmann</surname> <given-names>H.</given-names></name> <name><surname>Wickham</surname> <given-names>H.</given-names></name> <name><surname>Kafadar</surname> <given-names>K.</given-names></name></person-group> (<year>2017</year>). <article-title>Letter-value plots: boxplots for large data</article-title>. <source>J. Comput. Graphical Stat</source>. <volume>26</volume>, <fpage>469</fpage>&#x02013;<lpage>477</lpage>. <pub-id pub-id-type="doi">10.1080/10618600.2017.1305277</pub-id></citation>
</ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Houghton</surname> <given-names>C.</given-names></name> <name><surname>Kreuz</surname> <given-names>T.</given-names></name></person-group> (<year>2012</year>). <article-title>On the efficient calculation of van Rossum distances</article-title>. <source>Netw. Comput. Neural Syst</source>. <volume>23</volume>, <fpage>48</fpage>&#x02013;<lpage>58</lpage>. <pub-id pub-id-type="doi">10.3109/0954898X.2012.673048</pub-id><pub-id pub-id-type="pmid">22568695</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Houghton</surname> <given-names>C.</given-names></name> <name><surname>Sen</surname> <given-names>K.</given-names></name></person-group> (<year>2008</year>). <article-title>A new multineuron spike train metric</article-title>. <source>Neural Comput</source>. <volume>20</volume>, <fpage>1495</fpage>&#x02013;<lpage>1511</lpage>. <pub-id pub-id-type="doi">10.1162/neco.2007.10-06-350</pub-id><pub-id pub-id-type="pmid">18194108</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huber</surname> <given-names>A. E.</given-names></name> <name><surname>Liu</surname> <given-names>S.-C.</given-names></name></person-group> (<year>2019</year>). <article-title>Filtering of nonuniformly sampled bandlimited functions</article-title>. <source>IEEE Signal Process. Lett</source>. <volume>26</volume>, <fpage>1036</fpage>&#x02013;<lpage>1040</lpage>. <pub-id pub-id-type="doi">10.1109/LSP.2019.2918435</pub-id></citation>
</ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jeng</surname> <given-names>S.-T.</given-names></name> <name><surname>Ritchie</surname> <given-names>S. G.</given-names></name></person-group> (<year>2008</year>). <article-title>Real-time vehicle classification using inductive loop signature data</article-title>. <source>Transp. Res. Rec</source>. <volume>2086</volume>, <fpage>8</fpage>&#x02013;<lpage>22</lpage>. <pub-id pub-id-type="doi">10.3141/2086-02</pub-id></citation>
</ref>
<ref id="B21">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Kass</surname> <given-names>R. E.</given-names></name> <name><surname>Eden</surname> <given-names>U. T.</given-names></name> <name><surname>Brown</surname> <given-names>E. N.</given-names></name></person-group> (<year>2014</year>). <source>Analysis of Neural Data</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Springer New York</publisher-name>. <pub-id pub-id-type="doi">10.1007/978-1-4614-9602-1</pub-id></citation>
</ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kheradpisheh</surname> <given-names>S. R.</given-names></name> <name><surname>Masquelier</surname> <given-names>T.</given-names></name></person-group> (<year>2020</year>). <article-title>Temporal backpropagation for spiking neural networks with one spike per neuron</article-title>. <source>Int. J. Neural Syst</source>. <volume>30</volume>:<fpage>2050027</fpage>. <pub-id pub-id-type="doi">10.1142/S0129065720500276</pub-id><pub-id pub-id-type="pmid">32466691</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ki</surname> <given-names>Y.-K.</given-names></name> <name><surname>Baik</surname> <given-names>D.-K.</given-names></name></person-group> (<year>2006</year>). <article-title>Vehicle-classification algorithm for single-loop detectors using neural networks</article-title>. <source>IEEE Trans. Veh. Technol</source>. <volume>55</volume>, <fpage>1704</fpage>&#x02013;<lpage>1711</lpage>. <pub-id pub-id-type="doi">10.1109/TVT.2006.883726</pub-id></citation>
</ref>
<ref id="B24">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Klein</surname> <given-names>L. A.</given-names></name> <name><surname>Mills</surname> <given-names>M. K.</given-names></name> <name><surname>Gibson</surname> <given-names>D. R. P.</given-names></name></person-group> (<year>2006</year>). <source>Traffic Detector Handbook: Volume I. Technical Report, Federal Highway Administration</source>. <publisher-loc>McLean, VA</publisher-loc>: <publisher-name>Turner-Fairbank Highway Research Center</publisher-name>.</citation>
</ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kozat</surname> <given-names>S. S.</given-names></name> <name><surname>Guan</surname> <given-names>K. M.</given-names></name> <name><surname>Singer</surname> <given-names>A. C.</given-names></name></person-group> (<year>2013</year>). <article-title>Tracking the best level set in a level-crossing analog-to-digital converter</article-title>. <source>Digit. Signal Process</source>. <volume>23</volume>, <fpage>478</fpage>&#x02013;<lpage>487</lpage>. <pub-id pub-id-type="doi">10.1016/j.dsp.2012.07.018</pub-id></citation>
</ref>
<ref id="B26">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Kwon</surname> <given-names>T. M.</given-names></name></person-group> (<year>2005</year>). <article-title>&#x0201C;Route tracking of border crossing vehicles using inductance signatures of loop detectors,&#x0201D;</article-title> in <source>Proceedings of the 2005 IEEE International Workshop on Measurement Systems for Homeland Security, Contraband Detection and Personal Safety Workshop, 2005 (IMS 2005)</source> (<publisher-loc>Orlando, FL</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>103</fpage>&#x02013;<lpage>109</lpage>. <pub-id pub-id-type="doi">10.1109/MSHS.2005.1502566</pub-id></citation>
</ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>J. H.</given-names></name> <name><surname>Delbruck</surname> <given-names>T.</given-names></name> <name><surname>Pfeiffer</surname> <given-names>M.</given-names></name></person-group> (<year>2016</year>). <article-title>Training deep spiking neural networks using backpropagation</article-title>. <source>Front. Neurosci</source>. <volume>10</volume>:<fpage>508</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2016.00508</pub-id><pub-id pub-id-type="pmid">27877107</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lu</surname> <given-names>X.-Y.</given-names></name> <name><surname>Varaiya</surname> <given-names>P.</given-names></name> <name><surname>Horowitz</surname> <given-names>R.</given-names></name> <name><surname>Guo</surname> <given-names>Z.</given-names></name> <name><surname>Palen</surname> <given-names>J.</given-names></name></person-group> (<year>2012</year>). <article-title>Estimating traffic speed with single inductive loop event data</article-title>. <source>Transp. Res. Rec</source>. <volume>2308</volume>, <fpage>157</fpage>&#x02013;<lpage>166</lpage>. <pub-id pub-id-type="doi">10.3141/2308-17</pub-id></citation>
</ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mark</surname> <given-names>J.</given-names></name> <name><surname>Todd</surname> <given-names>T.</given-names></name></person-group> (<year>1981</year>). <article-title>A nonuniform sampling approach to data compression</article-title>. <source>IEEE Trans. Commun. COM</source> <volume>29</volume>, <fpage>24</fpage>&#x02013;<lpage>32</lpage>. <pub-id pub-id-type="doi">10.1109/TCOM.1981.1094872</pub-id></citation>
</ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marsza&#x00142;ek</surname> <given-names>Z.</given-names></name> <name><surname>Duda</surname> <given-names>K.</given-names></name></person-group> (<year>2020</year>). <article-title>Multifrequency vector measurement system for reliable vehicle magnetic profile assessment</article-title>. <source>Sensors</source> <volume>20</volume>:<fpage>4933</fpage>. <pub-id pub-id-type="doi">10.3390/s20174933</pub-id><pub-id pub-id-type="pmid">32878201</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marsza&#x00142;ek</surname> <given-names>Z.</given-names></name> <name><surname>Duda</surname> <given-names>K.</given-names></name> <name><surname>Piwowar</surname> <given-names>P.</given-names></name> <name><surname>Stencel</surname> <given-names>M.</given-names></name> <name><surname>Zeglen</surname> <given-names>T.</given-names></name> <name><surname>Izydorczyk</surname> <given-names>J.</given-names></name></person-group> (<year>2023</year>). <article-title>Load estimation of moving passenger cars using inductive-loop technology</article-title>. <source>Sensors</source> <volume>23</volume>:<fpage>2063</fpage>. <pub-id pub-id-type="doi">10.3390/s23042063</pub-id><pub-id pub-id-type="pmid">36850661</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Marsza&#x00142;ek</surname> <given-names>Z.</given-names></name> <name><surname>Sroka</surname> <given-names>R.</given-names></name></person-group> (<year>2016</year>). <article-title>&#x0201C;Signal fusion of changes in the inductive loop impedance components for vehicle axle detection,&#x0201D;</article-title> in <source>2016 21st International Conference on Methods and Models in Automation and Robotics (MMAR)</source> (<publisher-loc>Miedzyzdroje</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>25</fpage>&#x02013;<lpage>29</lpage>. <pub-id pub-id-type="doi">10.1109/MMAR.2016.7575082</pub-id></citation>
</ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marsza&#x00142;ek</surname> <given-names>Z.</given-names></name> <name><surname>Sroka</surname> <given-names>R.</given-names></name> <name><surname>Stencel</surname> <given-names>M.</given-names></name></person-group> (<year>2011</year>). <article-title>A new method of inductive sensors impedance measurement applied to the identification of vehicle parameters</article-title>. <source>Metrol. Meas. Syst</source>. <volume>18</volume>, <fpage>69</fpage>&#x02013;<lpage>76</lpage>. <pub-id pub-id-type="doi">10.2478/v10178-011-0007-2</pub-id></citation>
</ref>
<ref id="B34">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Marsza&#x00142;ek</surname> <given-names>Z.</given-names></name> <name><surname>Sroka</surname> <given-names>R.</given-names></name> <name><surname>Zeglen</surname> <given-names>T.</given-names></name></person-group> (<year>2015</year>). <article-title>&#x0201C;Inductive loop for vehicle axle detection from first concepts to the system based on changes in the sensor impedance components,&#x0201D;</article-title> in <source>2015 20th International Conference on Methods and Models in Automation and Robotics (MMAR)</source> (<publisher-loc>Miedzyzdroje</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>765</fpage>&#x02013;<lpage>769</lpage>. <pub-id pub-id-type="doi">10.1109/MMAR.2015.7283972</pub-id></citation>
</ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marsza&#x00142;ek</surname> <given-names>Z.</given-names></name> <name><surname>Zeglen</surname> <given-names>T.</given-names></name> <name><surname>Sroka</surname> <given-names>R.</given-names></name> <name><surname>Gajda</surname> <given-names>J.</given-names></name></person-group> (<year>2018</year>). <article-title>Inductive loop axle detector based on resistance and reactance vehicle magnetic profiles</article-title>. <source>Sensors</source> <volume>18</volume>:<fpage>2376</fpage>. <pub-id pub-id-type="doi">10.3390/s18072376</pub-id><pub-id pub-id-type="pmid">30037101</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Meta</surname> <given-names>S.</given-names></name> <name><surname>Cinsdikici</surname> <given-names>M. G.</given-names></name></person-group> (<year>2010</year>). <article-title>Vehicle-classification algorithm based on component analysis for single-loop inductive detector</article-title>. <source>IEEE Trans. Veh. Technol</source>. <volume>59</volume>, <fpage>2795</fpage>&#x02013;<lpage>2805</lpage>. <pub-id pub-id-type="doi">10.1109/TVT.2010.2049756</pub-id></citation>
</ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Midya</surname> <given-names>R.</given-names></name> <name><surname>Wang</surname> <given-names>Z.</given-names></name> <name><surname>Asapu</surname> <given-names>S.</given-names></name> <name><surname>Joshi</surname> <given-names>S.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name> <name><surname>Zhuo</surname> <given-names>Y.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Artificial neural network (ANN) to spiking neural network (SNN) converters based on diffusive memristors</article-title>. <source>Adv. Electron. Mater</source>. <volume>5</volume>:<fpage>1900060</fpage>. <pub-id pub-id-type="doi">10.1002/aelm.201900060</pub-id></citation>
</ref>
<ref id="B38">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Mi&#x0015B;kowicz</surname> <given-names>M.</given-names></name></person-group> (<year>2003</year>). <article-title>&#x0201C;The event-triggered sampling optimization criterion for distributed networked monitoring and control systems,&#x0201D;</article-title> in <source>IEEE International Conference on Industrial Technology, 2003, Volume 2</source> (<publisher-loc>Maribor</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>1083</fpage>&#x02013;<lpage>1088</lpage>. <pub-id pub-id-type="doi">10.1109/ICIT.2003.1290814</pub-id></citation>
</ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mi&#x0015B;kowicz</surname> <given-names>M.</given-names></name></person-group> (<year>2006</year>). <article-title>Send-on-delta concept: an event-based data reporting strategy</article-title>. <source>Sensors</source> <volume>6</volume>, <fpage>49</fpage>&#x02013;<lpage>63</lpage>. <pub-id pub-id-type="doi">10.3390/s6010049</pub-id></citation>
</ref>
<ref id="B40">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Mi&#x0015B;kowicz</surname> <given-names>M.</given-names></name></person-group> (<year>2015</year>). <source>Event-based Control and Signal Processing</source>. <publisher-loc>Boca Raton, FL</publisher-loc>: <publisher-name>CRC Press</publisher-name>.</citation>
</ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mochol&#x000ED; Belenguer</surname> <given-names>F.</given-names></name> <name><surname>Mochol&#x000ED; Salcedo</surname> <given-names>A.</given-names></name> <name><surname>Guill Iba&#x000F1;ez</surname> <given-names>A.</given-names></name> <name><surname>Mili&#x000E1;n S&#x000E1;nchez</surname> <given-names>V.</given-names></name></person-group> (<year>2019</year>). <article-title>Advantages offered by the double magnetic loops versus the conventional single ones</article-title>. <source>PLoS ONE</source> <volume>14</volume>:<fpage>e0211626</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0211626</pub-id><pub-id pub-id-type="pmid">30753200</pub-id></citation></ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mochol&#x000ED;-Salcedo</surname> <given-names>A.</given-names></name> <name><surname>Arroyo-N&#x000FA;&#x000F1;ez</surname> <given-names>J. H.</given-names></name> <name><surname>Mili&#x000E1;n-S&#x000E1;nchez</surname> <given-names>V. M.</given-names></name> <name><surname>Palomo-Anaya</surname> <given-names>M. J.</given-names></name> <name><surname>Arroyo-N&#x000FA;&#x000F1;ez</surname> <given-names>A.</given-names></name></person-group> (<year>2017</year>). <article-title>Magnetic field generated by the loops used in traffic control systems</article-title>. <source>IEEE Trans. Intell. Transp. Syst</source>. <volume>18</volume>, <fpage>2126</fpage>&#x02013;<lpage>2136</lpage>. <pub-id pub-id-type="doi">10.1109/TITS.2016.2632972</pub-id><pub-id pub-id-type="pmid">30753200</pub-id></citation></ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Moser</surname> <given-names>B. A.</given-names></name></person-group> (<year>2017</year>). <article-title>Similarity recovery from threshold-based sampling under general conditions</article-title>. <source>IEEE Trans. Signal Process</source>. <volume>65</volume>, <fpage>4645</fpage>&#x02013;<lpage>4654</lpage>. <pub-id pub-id-type="doi">10.1109/TSP.2017.2712121</pub-id></citation>
</ref>
<ref id="B44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Moser</surname> <given-names>B. A.</given-names></name> <name><surname>Werzi</surname> <given-names>A.</given-names></name> <name><surname>Lunglmayr</surname> <given-names>M.</given-names></name></person-group> (<year>2025</year>). <article-title>Integrate-and-fire from a mathematical and signal processing perspective</article-title>. <source>arXiv</source> [Preprint]. arXiv:2501.11453. <pub-id pub-id-type="doi">10.48550/arXiv.2501.11453</pub-id></citation>
</ref>
<ref id="B45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mostafa</surname> <given-names>H.</given-names></name></person-group> (<year>2018</year>). <article-title>Supervised learning based on temporal coding in spiking neural networks</article-title>. <source>IEEE Trans. Neural Netw. Learn. Syst</source>. <volume>29</volume>, <fpage>3227</fpage>&#x02013;<lpage>3235</lpage>. <pub-id pub-id-type="doi">10.1109/TNNLS.2017.2778987</pub-id></citation>
</ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mozafari</surname> <given-names>M.</given-names></name> <name><surname>Ganjtabesh</surname> <given-names>M.</given-names></name> <name><surname>Nowzari-Dalini</surname> <given-names>A.</given-names></name> <name><surname>Thorpe</surname> <given-names>S. J.</given-names></name> <name><surname>Masquelier</surname> <given-names>T.</given-names></name></person-group> (<year>2019</year>). <article-title>Bio-inspired digit recognition using reward-modulated spike-timing-dependent plasticity in deep convolutional networks</article-title>. <source>Pattern Recognit</source>. <volume>94</volume>, <fpage>87</fpage>&#x02013;<lpage>95</lpage>. <pub-id pub-id-type="doi">10.1016/j.patcog.2019.05.015</pub-id></citation>
</ref>
<ref id="B47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ndoye</surname> <given-names>M.</given-names></name> <name><surname>Totten</surname> <given-names>V. F.</given-names></name> <name><surname>Krogmeier</surname> <given-names>J. V.</given-names></name> <name><surname>Bullock</surname> <given-names>D. M.</given-names></name></person-group> (<year>2011</year>). <article-title>Sensing and signal processing for vehicle reidentification and travel time estimation</article-title>. <source>IEEE Trans. Intell. Transp. Syst</source>. <volume>12</volume>, <fpage>119</fpage>&#x02013;<lpage>131</lpage>. <pub-id pub-id-type="doi">10.1109/TITS.2010.2092769</pub-id></citation>
</ref>
<ref id="B48">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Neftci</surname> <given-names>E. O.</given-names></name> <name><surname>Mostafa</surname> <given-names>H.</given-names></name> <name><surname>Zenke</surname> <given-names>F.</given-names></name></person-group> (<year>2019</year>). <article-title>Surrogate gradient learning in spiking neural networks: bringing the power of gradient-based optimization to spiking neural networks</article-title>. <source>IEEE Signal Process. Mag</source>. <volume>36</volume>, <fpage>51</fpage>&#x02013;<lpage>63</lpage>. <pub-id pub-id-type="doi">10.1109/MSP.2019.2931595</pub-id></citation>
</ref>
<ref id="B49">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nunes</surname> <given-names>J. D.</given-names></name> <name><surname>Carvalho</surname> <given-names>M.</given-names></name> <name><surname>Carneiro</surname> <given-names>D.</given-names></name> <name><surname>Cardoso</surname> <given-names>J. S.</given-names></name></person-group> (<year>2022</year>). <article-title>Spiking neural networks: a survey</article-title>. <source>IEEE Access</source> <volume>10</volume>, <fpage>60738</fpage>&#x02013;<lpage>60764</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2022.3179968</pub-id></citation>
</ref>
<ref id="B50">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Oh</surname> <given-names>S.</given-names></name> <name><surname>Ritchie</surname> <given-names>S. G.</given-names></name> <name><surname>Oh</surname> <given-names>C.</given-names></name></person-group> (<year>2002</year>). <article-title>Real-time traffic measurement from single loop inductive signatures</article-title>. <source>Transp. Res. Rec</source>. <volume>1804</volume>, <fpage>98</fpage>&#x02013;<lpage>106</lpage>. <pub-id pub-id-type="doi">10.3141/1804-14</pub-id></citation>
</ref>
<ref id="B51">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pabian</surname> <given-names>M.</given-names></name> <name><surname>Rzepka</surname> <given-names>D.</given-names></name> <name><surname>Pawlak</surname> <given-names>M.</given-names></name></person-group> (<year>2024</year>). <article-title>Iteration over event space in time-to-first-spike spiking neural networks for Twitter bot classification</article-title>. <source>arXiv</source> [Preprint]. arXiv:2407.08746. <pub-id pub-id-type="doi">10.48550/arXiv.2407.08746</pub-id></citation>
</ref>
<ref id="B52">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pfeiffer</surname> <given-names>M.</given-names></name> <name><surname>Pfeil</surname> <given-names>T.</given-names></name></person-group> (<year>2018</year>). <article-title>Deep learning with spiking neurons: opportunities and challenges</article-title>. <source>Front. Neurosci</source>. <volume>12</volume>:<fpage>774</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2018.00774</pub-id><pub-id pub-id-type="pmid">30410432</pub-id></citation></ref>
<ref id="B53">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rasmussen</surname> <given-names>D.</given-names></name></person-group> (<year>2019</year>). <article-title>NengoDL: combining deep learning and neuromorphic modelling methods</article-title>. <source>Neuroinformatics</source> <volume>17</volume>, <fpage>611</fpage>&#x02013;<lpage>628</lpage>. <pub-id pub-id-type="doi">10.1007/s12021-019-09424-z</pub-id><pub-id pub-id-type="pmid">30972529</pub-id></citation></ref>
<ref id="B54">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ravanshad</surname> <given-names>N.</given-names></name> <name><surname>Rezaee-Dehsorkh</surname> <given-names>H.</given-names></name> <name><surname>Lotfi</surname> <given-names>R.</given-names></name> <name><surname>Lian</surname> <given-names>Y.</given-names></name></person-group> (<year>2013</year>). <article-title>A level-crossing based QRS-detection algorithm for wearable ECG sensors</article-title>. <source>IEEE J. Biomed. Health Inform</source>. <volume>18</volume>, <fpage>183</fpage>&#x02013;<lpage>192</lpage>. <pub-id pub-id-type="doi">10.1109/JBHI.2013.2274809</pub-id><pub-id pub-id-type="pmid">24403416</pub-id></citation></ref>
<ref id="B55">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rueckauer</surname> <given-names>B.</given-names></name> <name><surname>Lungu</surname> <given-names>I.-A.</given-names></name> <name><surname>Hu</surname> <given-names>Y.</given-names></name> <name><surname>Pfeiffer</surname> <given-names>M.</given-names></name> <name><surname>Liu</surname> <given-names>S.-C.</given-names></name></person-group> (<year>2017</year>). <article-title>Conversion of continuous-valued deep networks to efficient event-driven networks for image classification</article-title>. <source>Front. Neurosci</source>. <volume>11</volume>:<fpage>682</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2017.00682</pub-id><pub-id pub-id-type="pmid">29375284</pub-id></citation></ref>
<ref id="B56">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rzepka</surname> <given-names>D.</given-names></name> <name><surname>Mi&#x0015B;kowicz</surname> <given-names>M.</given-names></name> <name><surname>Ko&#x0015B;cielnik</surname> <given-names>D.</given-names></name> <name><surname>Thao</surname> <given-names>N. T.</given-names></name></person-group> (<year>2018</year>). <article-title>Reconstruction of signals from level-crossing samples using implicit information</article-title>. <source>IEEE Access</source> <volume>6</volume>, <fpage>35001</fpage>&#x02013;<lpage>35011</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2018.2839186</pub-id></citation>
</ref>
<ref id="B57">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rzepka</surname> <given-names>D.</given-names></name> <name><surname>Pawlak</surname> <given-names>M.</given-names></name> <name><surname>Ko&#x0015B;cielnik</surname> <given-names>D.</given-names></name> <name><surname>Mi&#x0015B;kowicz</surname> <given-names>M.</given-names></name></person-group> (<year>2017</year>). <article-title>Bandwidth estimation from multiple level-crossings of stochastic signals</article-title>. <source>IEEE Trans. Signal Process</source>. <volume>65</volume>, <fpage>2488</fpage>&#x02013;<lpage>2502</lpage>. <pub-id pub-id-type="doi">10.1109/TSP.2017.2664041</pub-id></citation>
</ref>
<ref id="B58">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Selva</surname> <given-names>J.</given-names></name></person-group> (<year>2012</year>). <article-title>Efficient sampling of band-limited signals from sine wave crossings</article-title>. <source>IEEE Trans. Signal Process</source>. <volume>60</volume>, <fpage>503</fpage>&#x02013;<lpage>508</lpage>. <pub-id pub-id-type="doi">10.1109/TSP.2011.2170171</pub-id></citation>
</ref>
<ref id="B59">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Senay</surname> <given-names>S.</given-names></name> <name><surname>Chaparro</surname> <given-names>L. F.</given-names></name> <name><surname>Sun</surname> <given-names>M.</given-names></name> <name><surname>Sclabassi</surname> <given-names>R. J.</given-names></name></person-group> (<year>2010</year>). <article-title>&#x0201C;Adaptive level-crossing sampling and reconstruction,&#x0201D;</article-title> in <source>Proc. of European Signal Processing Conference EUSIPCO 2010</source> (<publisher-loc>Aalborg</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>1296</fpage>&#x02013;<lpage>1300</lpage>.<pub-id pub-id-type="pmid">30761433</pub-id></citation></ref>
<ref id="B60">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stein</surname> <given-names>R. B.</given-names></name></person-group> (<year>1965</year>). <article-title>A theoretical analysis of neuronal variability</article-title>. <source>Biophys. J</source>. <volume>5</volume>, <fpage>173</fpage>&#x02013;<lpage>194</lpage>. <pub-id pub-id-type="doi">10.1016/S0006-3495(65)86709-1</pub-id><pub-id pub-id-type="pmid">14268952</pub-id></citation></ref>
<ref id="B61">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>St&#x000F6;ckl</surname> <given-names>C.</given-names></name> <name><surname>Maass</surname> <given-names>W.</given-names></name></person-group> (<year>2021</year>). <article-title>Optimized spiking neurons can classify images with high accuracy through temporal coding with two spikes</article-title>. <source>Nat. Mach. Intell</source>. <volume>3</volume>, <fpage>230</fpage>&#x02013;<lpage>238</lpage>. <pub-id pub-id-type="doi">10.1038/s42256-021-00311-4</pub-id></citation>
</ref>
<ref id="B62">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sun</surname> <given-names>C.</given-names></name> <name><surname>Ritchie</surname> <given-names>S. G.</given-names></name></person-group> (<year>1999</year>). <article-title>Individual vehicle speed estimation using single loop inductive waveforms</article-title>. <source>J. Transp. Eng</source>. <volume>125</volume>, <fpage>531</fpage>&#x02013;<lpage>538</lpage>. <pub-id pub-id-type="doi">10.1061/(ASCE)0733-947X(1999)125:6(531)</pub-id></citation>
</ref>
<ref id="B63">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tang</surname> <given-names>V.</given-names></name> <name><surname>Gomes</surname> <given-names>C.</given-names></name> <name><surname>Lucani</surname> <given-names>D. E.</given-names></name></person-group> (<year>2025</year>). <article-title>Precision on demand: propositional logic for event-trigger threshold regulation</article-title>. <source>IEEE Internet Things J</source>. <volume>12</volume>, <fpage>2674</fpage>&#x02013;<lpage>2689</lpage>. <pub-id pub-id-type="doi">10.1109/JIOT.2024.3476922</pub-id></citation>
</ref>
<ref id="B64">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tayarani-Najaran</surname> <given-names>M.-H.</given-names></name> <name><surname>Schmuker</surname> <given-names>M.</given-names></name></person-group> (<year>2021</year>). <article-title>Event-based sensing and signal processing in the visual, auditory, and olfactory domain: a review</article-title>. <source>Front. Neural Circuits</source> <volume>15</volume>:<fpage>610446</fpage>. <pub-id pub-id-type="doi">10.3389/fncir.2021.610446</pub-id><pub-id pub-id-type="pmid">34135736</pub-id></citation></ref>
<ref id="B65">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thao</surname> <given-names>N. T.</given-names></name> <name><surname>Rzepka</surname> <given-names>D.</given-names></name> <name><surname>Mi&#x0015B;kowicz</surname> <given-names>M.</given-names></name></person-group> (<year>2023</year>). <article-title>Bandlimited signal reconstruction from leaky integrate-and-fire encoding using POCS</article-title>. <source>IEEE Trans. Signal Process</source>. <volume>71</volume>, <fpage>1464</fpage>&#x02013;<lpage>1479</lpage>. <pub-id pub-id-type="doi">10.1109/TSP.2023.3256269</pub-id></citation>
</ref>
<ref id="B66">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tieleman</surname> <given-names>T.</given-names></name> <name><surname>Hinton</surname> <given-names>G.</given-names></name></person-group> (<year>2012</year>). <article-title>Lecture 6.5-RMSprop: divide the gradient by a running average of its recent magnitude</article-title>. <source>COURSERA</source> <volume>4</volume>, <fpage>26</fpage>&#x02013;<lpage>31</lpage>.</citation>
</ref>
<ref id="B67">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Trakimas</surname> <given-names>M.</given-names></name> <name><surname>Sonkusale</surname> <given-names>S. R.</given-names></name></person-group> (<year>2011</year>). <article-title>An adaptive resolution asynchronous ADC architecture for data compression in energy constrained sensing applications</article-title>. <source>IEEE Trans. Circuits Syst. I Regular Papers</source> <volume>58</volume>, <fpage>921</fpage>&#x02013;<lpage>934</lpage>. <pub-id pub-id-type="doi">10.1109/TCSI.2010.2092132</pub-id></citation>
</ref>
<ref id="B68">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tsividis</surname> <given-names>Y.</given-names></name></person-group> (<year>2003</year>). <article-title>Continuous-time digital signal processing</article-title>. <source>Electron. Lett</source>. <volume>39</volume>, <fpage>1551</fpage>&#x02013;<lpage>1552</lpage>. <pub-id pub-id-type="doi">10.1049/el:20031015</pub-id></citation>
</ref>
<ref id="B69">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>van Rossum</surname> <given-names>M. C. W.</given-names></name></person-group> (<year>2001</year>). <article-title>A novel spike distance</article-title>. <source>Neural Comput</source>. <volume>13</volume>, <fpage>751</fpage>&#x02013;<lpage>763</lpage>. <pub-id pub-id-type="doi">10.1162/089976601300014321</pub-id><pub-id pub-id-type="pmid">11255567</pub-id></citation></ref>
<ref id="B70">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>H.</given-names></name> <name><surname>Schembari</surname> <given-names>F.</given-names></name> <name><surname>Mi&#x0015B;kowicz</surname> <given-names>M.</given-names></name> <name><surname>Staszewski</surname> <given-names>R. B.</given-names></name></person-group> (<year>2018</year>). <article-title>An adaptive-resolution quasi-level-crossing-sampling ADC based on residue quantization in 28-nm CMOS</article-title>. <source>IEEE Solid-State Circuits Lett</source>. <volume>1</volume>, <fpage>178</fpage>&#x02013;<lpage>181</lpage>. <pub-id pub-id-type="doi">10.1109/LSSC.2019.2899723</pub-id></citation>
</ref>
<ref id="B71">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname> <given-names>Y.</given-names></name> <name><surname>Deng</surname> <given-names>L.</given-names></name> <name><surname>Li</surname> <given-names>G.</given-names></name> <name><surname>Zhu</surname> <given-names>J.</given-names></name> <name><surname>Shi</surname> <given-names>L.</given-names></name></person-group> (<year>2018</year>). <article-title>Spatio-temporal backpropagation for training high-performance spiking neural networks</article-title>. <source>Front. Neurosci</source>. <volume>12</volume>:<fpage>331</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2018.00331</pub-id><pub-id pub-id-type="pmid">29875621</pub-id></citation></ref>
<ref id="B72">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wunderlich</surname> <given-names>T. C.</given-names></name> <name><surname>Pehle</surname> <given-names>C.</given-names></name></person-group> (<year>2021</year>). <article-title>Event-based backpropagation can compute exact gradients for spiking neural networks</article-title>. <source>Sci. Rep</source>. <volume>11</volume>:<fpage>12829</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-021-91786-z</pub-id><pub-id pub-id-type="pmid">34145314</pub-id></citation></ref>
<ref id="B73">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xing</surname> <given-names>Y.</given-names></name> <name><surname>Di Caterina</surname> <given-names>G.</given-names></name> <name><surname>Soraghan</surname> <given-names>J.</given-names></name></person-group> (<year>2020</year>). <article-title>A new spiking convolutional recurrent neural network (SCRNN) with applications to event-based hand gesture recognition</article-title>. <source>Front. Neurosci</source>. <volume>14</volume>:<fpage>590164</fpage>. <pub-id pub-id-type="doi">10.3389/fnins.2020.590164</pub-id><pub-id pub-id-type="pmid">33324153</pub-id></citation></ref>
<ref id="B74">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zanoli</surname> <given-names>S.</given-names></name> <name><surname>Ansaloni</surname> <given-names>G.</given-names></name> <name><surname>Teijeiro</surname> <given-names>T.</given-names></name> <name><surname>Atiena</surname> <given-names>D.</given-names></name></person-group> (<year>2025</year>). <article-title>Mathematical model of parameters relevance in adaptive level-crossing sampling for electrocardiogram signals</article-title>. <source>arXiv</source> [Preprint]. arXiv:2501.10829. <pub-id pub-id-type="doi">10.48550/arXiv.2501.10829</pub-id></citation>
</ref>
<ref id="B75">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>W.</given-names></name> <name><surname>Li</surname> <given-names>P.</given-names></name></person-group> (<year>2020</year>). <article-title>Temporal spike sequence learning via backpropagation for deep spiking neural networks</article-title>. <source>Adv. Neural Inf. Process. Syst</source>. <volume>33</volume>, <fpage>12022</fpage>&#x02013;<lpage>12033</lpage>.<pub-id pub-id-type="pmid">37512665</pub-id></citation></ref>
<ref id="B76">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname> <given-names>S.</given-names></name> <name><surname>Li</surname> <given-names>X.</given-names></name> <name><surname>Chen</surname> <given-names>Y.</given-names></name> <name><surname>Chandrasekaran</surname> <given-names>S. T.</given-names></name> <name><surname>Sanyal</surname> <given-names>A.</given-names></name></person-group> (<year>2021</year>). <article-title>Temporal-coded deep spiking neural network with easy training and robust performance</article-title>. <source>Proc. AAAI Conf. Artif. Intell</source>. <volume>35</volume>, <fpage>11143</fpage>&#x02013;<lpage>11151</lpage>. <pub-id pub-id-type="doi">10.1609/aaai.v35i12.17329</pub-id></citation>
</ref>
<ref id="B77">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhu</surname> <given-names>Y.</given-names></name> <name><surname>Yu</surname> <given-names>Z.</given-names></name> <name><surname>Fang</surname> <given-names>W.</given-names></name> <name><surname>Xie</surname> <given-names>X.</given-names></name> <name><surname>Huang</surname> <given-names>T.</given-names></name> <name><surname>Masquelier</surname> <given-names>T.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Training spiking neural networks with event-driven backpropagation</article-title>. <source>Adv. Neural Inf. Process. Syst</source>. <volume>35</volume>, <fpage>30528</fpage>&#x02013;<lpage>30541</lpage>.</citation>
</ref>
</ref-list>
</back>
</article>