A time-causal and time-recursive scale-covariant scale-space representation of temporal signals and past time

This article presents an overview of a theory for performing temporal smoothing on temporal signals in such a way that: (i) temporally smoothed signals at coarser temporal scales are guaranteed to constitute simplifications of corresponding temporally smoothed signals at any finer temporal scale (including the original signal) and (ii) the temporal smoothing process is both time-causal and time-recursive, in the sense that it does not require access to future information and can be performed with no other temporal memory buffer of the past than the resulting smoothed temporal scale-space representations themselves. For specific subsets of parameter settings for the classes of linear and shift-invariant temporal smoothing operators that obey this property, it is shown how temporal scale covariance can be additionally obtained, guaranteeing that if the temporal input signal is rescaled by a uniform temporal scaling factor, then also the resulting temporal scale-space representations of the rescaled temporal signal will constitute mere rescalings of the temporal scale-space representations of the original input signal, complemented by a shift along the temporal scale dimension. The resulting time-causal limit kernel that obeys this property constitutes a canonical temporal kernel for processing temporal signals in real-time scenarios when the regular Gaussian kernel cannot be used, because of its non-causal access to information from the future, and we cannot additionally require the temporal smoothing process to comprise a complementary memory of the past beyond the information contained in the temporal smoothing process itself, which in this way also serves as a multi-scale temporal memory of the past. We describe how the time-causal limit kernel relates to previously used temporal models, such as Koenderink’s scale-time kernels and the ex-Gaussian kernel. We do also give an overview of how the time-causal limit kernel can be used for modelling the temporal processing in models for spatio-temporal and spectro-temporal receptive fields, and how it more generally has a high potential for modelling neural temporal response functions in a purely time-causal and time-recursive way, that can also handle phenomena at multiple temporal scales in a theoretically well-founded manner. We detail how this theory can be efficiently implemented for discrete data, in terms of a set of recursive filters coupled in cascade. Hence, the theory is generally applicable for both: (i) modelling continuous temporal phenomena over multiple temporal scales and (ii) digital processing of measured temporal signals in real time. We conclude by stating implications of the theory for modelling temporal phenomena in biological, perceptual, neural and memory processes by mathematical models, as well as implications regarding the philosophy of time and perceptual agents. Specifically, we propose that for A-type theories of time, as well as for perceptual agents, the notion of a non-infinitesimal inner temporal scale of the temporal receptive fields has to be included in representations of the present, where the inherent nonzero temporal delay of such time-causal receptive fields implies a need for incorporating predictions from the actual time-delayed present in the layers of a perceptual hierarchy, to make it possible for a representation of the perceptual present to constitute a representation of the environment with timing properties closer to the actual present.

For specific subsets of parameter settings for the classes of linear and shift-invariant temporal smoothing operators that obey this property, it is shown how temporal scale covariance can be additionally obtained, guaranteeing that if the temporal input signal is rescaled by a uniform temporal scaling factor, then also the resulting temporal scale-space representations of the rescaled temporal signal will constitute mere rescalings of the temporal scale-space representations of the original input signal, complemented by a shift along the temporal scale dimension.The resulting time-causal limit kernel that obeys this property constitutes a canonical temporal kernel for processing temporal signals in realtime scenarios when the regular Gaussian kernel cannot be used, because of its non-causal access to information from the future, and we cannot additionally require the temporal smoothing process to comprise a complementary memory of the past beyond the information contained in the tempo- The support from the Swedish Research Council (Contract 2018-03586) is gratefully acknowledged.
A shortened version of this article is published Open Access in Biological Cybernetics under the Creative Commons Attribution 4.0 International License, see https://doi.org/10.1007/s00422-022-00953-6.Tony Lindeberg, Computational Brain Science Lab, Division of Computational Science and Technology, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden.E-mail: tony@kth.seORCID: 0000-0002-9081-2170 ral smoothing process itself, which in this way also serves as a multi-scale temporal memory of the past.
We describe how the time-causal limit kernel relates to previously used temporal models, such as Koenderink's scaletime kernels and the ex-Gaussian kernel.We do also give an overview of how the time-causal limit kernel can be used for modelling the temporal processing in models for spatiotemporal and spectro-temporal receptive fields, and how it more generally has a high potential for modelling neural temporal response functions in a purely time-causal and timerecursive way, that can also handle phenomena at multiple temporal scales in a theoretically well-founded manner.
We detail how this theory can be efficiently implemented for discrete data, in terms of a set of recursive filters coupled in cascade.Hence, the theory is generally applicable for both: (i) modelling continuous temporal phenomena over multiple temporal scales and (ii) digital processing of measured temporal signals in real time.
We conclude by stating implications of the theory for modelling temporal phenomena in biological, perceptual, neural and memory processes by mathematical models, as well as implications regarding the philosophy of time and perceptual agents.Specifically, we propose that for A-type theories of time, as well as for perceptual agents, the notion of a non-infinitesimal inner temporal scale of the temporal receptive fields has to be included in representations of the present, where the inherent non-zero temporal delay of such time-causal receptive fields implies a need for incorporating predictions from the actual time-delayed present in the layers of a perceptual hierarchy, to make it possible for a representation of the perceptual present to constitute a representation of the environment with timing properties closer to the actual present.

Introduction
When processing time-dependent measurement signals, there is often a need to perform temporal smoothing prior to more refined data analysis.A commonly stated general motivation for this need is to suppress measurement noise, often based on the assumption that there is a well-defined underlying noise free signal that has been corrupted with some amount of measurement noise.
A more fundamental approach to take on the need for performing temporal smoothing of temporal signals is to follow a multi-scale approach, based on the observation that measurements performed on real-world data may reflect different types of temporal structures at different temporal scales.In other words, even for the underlying noise free signal in the above signal+noise model, it may hold that the data reflect different types of underlying physical or biological processes at different temporal scales.The measurement process itself, by which a non-infinitesimal amount of energy needs to be integrated over some non-infinitesimal temporal duration on the physical sensor, does in this respect define an inner temporal scale of the measurements, beyond which there is no way to resolve temporal phenomena that occur faster than this inner temporal scale.Any real-world physical measurement does in this respect involve an inherent notion of temporal scale. 1pecifically, in the areas of image processing, computer vision, machine listening2 and computational modelling of visual and auditory perception, this need is well understood, and has lead to multi-scale approaches for spatial, spatiotemporal and spectro-temporal receptive fields expressed in terms of multi-scale representations over the spatial, spectral and temporal domains, where specifically the theoretical framework known as scale-space theory is based upon solid theory in terms of axiomatic derivations concerning how the multi-scale processing operations should be performed (Iijima 1962;Witkin 1983;Koenderink 1984;Koenderink andvan Doorn 1987, 1992;Lindeberg 1993bLindeberg , 1994Lindeberg , 2011Lindeberg , 2013b;;Florack 1997;Sporring et al. 1997;Weickert et al. 1999;ter Haar Romeny 2003).It has also been found that biological perception, memory and cognition has developed biological processes at multiple temporal scales (DeAngelis et al. 1995(DeAngelis et al. , 2004;;Gütig and Sompolinsky 2006;Gentner 2008); Holcombe 2009; Goldman 2009;Gauthier et al. 2012;Atencio and Schreiner 2012;Chait et al. 2015;Teng et al. 2016;Buzsáki and Llinás 2017;Tsao et al. 2018;Osman et al. 2018;Latimer et al. 2019;Bright et al. 2020;Cavanagh et al. 2020;Monsa et al. 2020;Spitmaan et al. 2020;Howard and Hasselmo 2020;Howard 2021;Guo et al. 2021;Miri et al. 2022); see Section 7.3 for a more detailed retrospective review.
The subject of this article is to describe a theoretical framework for representing temporal signals at multiple temporal scales, intended for a more general audience without background in these areas and with the focus on the temporal domain only, thus without the complementary spatial or spectral domains that this theory has previously been combined with for expressing spatio-temporal and spectrotemporal receptive fields (Lindeberg and Fagerström 1996;Lindeberg 1997aLindeberg , 2016Lindeberg , 2017Lindeberg , 2018aLindeberg , 2018bLindeberg , 2021b;;Lindeberg andFriberg 2015b, 2015a).This theoretical framework, referred to as temporal scale-space theory, guarantees noncreation of the temporal structures with increasing temporal scales, in the sense that it ensures that a temporal representation at any coarser temporal scale constitutes a simplification of a temporal representation at any finer temporal scale, in the respect that the number of local temporal extrema, alternatively the number of temporal zero-crossings, is guaranteed to not increase from finer to coarser temporal scales.
Additionally, these temporal scale-space representations are time-causal, in the sense that they do not require access to future data, and are time-recursive, in the respect that the temporal representation at the next temporal moment can be computed with no other additional memory of the past than the temporal scale-space representation itself.For a specific choice of temporal scale-space kernel, referred to as the time-causal limit kernel, the temporal scale-space representations are also scale covariant, meaning that the set of temporal scale-space representations is closed under temporal rescalings of the input.A rescaling of the input signal by a uniform scaling factor merely corresponds to a rescaling of the temporal scale-space representations complemented by a shift of the temporal scale levels in the temporal scale-space representation.In this way, the temporal scale-space representation ensures an internally consistent way of processing temporal signals that may be subject to temporal scaling transformations, by phenomena or events that may occur faster or slower in the world.
A main purpose of this article is to describe this theory in a self-contained manner, without need for the reader to digest the original references, where the information is distributed over several papers, and may require a substantial effort for a reader not previously familiar with this framework, to get an updated view of the latest version of this theory. 3Furthermore, we will describe explicit relations to other previously used temporal models, such as Koenderink's scale-time kernels (Koenderink 1988) and the ex-Gaussian model (Grushka 1972;Bright et al. 2020), making it possible to transfer modelling results from those temporal models to the time-causal limit kernel described in this article.
We will also relate the presented temporal scale-space theory to other approaches for processing signals at multiple temporal scales, such as wavelet analysis and timefrequency analysis.Specifically, we will outline how the temporal derivatives of the proposed time-causal limit kernel described and analyzed in this article allow for fully timecausal and time-recursive wavelet analysis methods, without need for additional temporal buffering, and thus enabling minimal temporal response times in a time-critical context.We will also outline how a complex-valued extension of the proposed time-causal limit kernel can be seen as a timecausal analogue of Gabor functions, thus allowing for capturing essentially similar transformations of temporal signals as for the family of Gabor functions, and thereby providing a way to define a scale-covariant time-frequency representation over a time-causal temporal domain, which by a slight modification can also be extended to additionally being implemented in terms strictly time-recursive operations.
Additionally, we will describe implications of using this theory for modelling perceptual, neural and memory processes in biological systems by mathematical models, as well as implications of the theory with regard to the philosophy of time and perceptual agents.Specifically, we will argue that when modelling a perceptual representation of the present, it is essential to include the inner temporal scales of the perceptual processes that lead to any percept, where the inherent temporal delays of such time-causal operations imply that a representation of the present will de facto constitute a representation of some temporal intervals in the past, unless complemented by prediction processes to enable better timing properties of a perceptual agent that interacts with a dynamic world.

Structure of this article
This paper is organized as follows: Section 2 introduces the problem of constructing a temporal scale-space representation, as constituting a multi-scale representation of temporal signals, with the property that a measure of the amount of structure in the signal, quantified as the number of local extrema over time, must not increase from any finer to any coarser temporal scale.A complete classification of the time-causal convolution kernels that enable this property is given, and it is shown that the only possible time-causal scale-space kernels over a continuous temporal domain consist of truncated exponential kernels coupled in cascade.
Section 3 then adds a complementary condition on this structure, in terms of temporal scale covariance, and meaning that if the temporal input signal is rescaled by a uniform temporal scaling factor, then the result of temporal scalespace filtering of this kernel should also be a mere rescaling of the result of performing temporal scale-space filtering on the input signal, complemented by a shift in the temporal scale channels and a possibly complementary shift in the magnitude of the signal.It is shown that a specific kernel, the time-causal limit kernel, defined from an infinite convolution of truncated exponential kernels in cascade, with specially chosen time constants, obeys temporal scale covariance.We do also show how this time-causal limit kernel relates to previously used temporal models, such as Koenderink's scale-time kernels and the ex-Gaussian kernel.
In Section 4, we complement the above treatment for continuous signals with a corresponding discrete theory, ensuring that the number of local extrema in a discrete signal is also guaranteed to not increase from any finer to any coarser temporal scale.The discrete analogue of the truncated exponential kernels are first-order recursive filters, which are then again coupled in cascade.Section 5 furthermore generalizes the above theory from temporal smoothing of a raw temporal signal, to the computation of temporal scale-space derivatives, which measure the amount of change in the signal with respect to any level of temporal scale.Section 6 outlines how the proposed temporal scale-space representation is related to other approaches for handling temporal signals at multiple temporal scales, specifically wavelet analysis and time-frequency analysis, with conceptual extensions of these notions with respect to strictly time-causal and timerecursive operations for real-time applications.
Section 7 describes how this general theory can be used for modelling time-dependent processes and mechanisms in perceptual and neural systems, with emphasis on spatio-temporal and spectro-temporal receptive fields as well as temporal memory processes.Section 8 outlines more general implications of the theory with regard to the philosophy of time and how time is handled by a perceptual agent.Specifically, we develop how the inner temporal scale associated with any biophysical measurement of time-dependent phenomena implies that a non-infinitesimal inner temporal scale needs to be included in a representation of the perceptual present, and also that the non-zero temporal delay of such time-causal kernels implies that a biophysical representation of the present will de facto constitute a representation of what has occurred over some temporal intervals in the past, in turn implying a need for prediction mechanisms to extrapolate the de facto time-delayed representation of the present into a better predicted representation of the actual present.
Section 9 gives a retrospective historic overview of the different parts of temporal scale-space theory that this paper is based on, follows and extends, as well as a conceptual overview of some of the main contributions to temporal scale-space theory made in this article.Finally, Section 10 summarizes some of the main results. 1 The main idea of a scale-space representation is to, given any input signal f (t), create a set of derived signals L(t; τ ) intended to represent the information in the original signal at a set of coarser levels of scale τ , with L(t; 0) = f (t).These derived signals should preferably constitute true simplifications of each other, in such a way that the signal at a coarser level of scale does not contain more structures or information than any signal at any finer level of scale.Over spatial image domains, the notion of scale-space representation has been extensively studied, where several axiomatic derivations have shown that the Gaussian kernel and its corresponding Gaussian derivatives constitute a canonical class of convolution kernels for generating a spatial scale-space representation and have also been demonstrated to constitute a suitable basis of image primitives for computing different types of features from spatial image data.In this paper, we develop the associated notion of temporal scale-space theory, based on the additional constraints that (i) the temporal scale-space kernels are not allowed to access information from the future in relation to any time moment and that (ii) the computations should be possible to perform in a purely time-recursive manner, implying no other need for a temporal memory of the past than the temporal scale-space representation itself.Furthermore, we add a complementary requirement of (iii) temporal scale covariance, meaning that under temporal scaling variations of the input, the temporal scale-space representations should also constitute mere temporal rescalings of the temporal scale-space representation computed from the original temporal signal before the temporal rescaling operation, complemented by a shift along the temporal scale axis.

Time-causal and time-recursive scale-space model for temporal signals
The problem that we consider is that we are given a temporal signal f (t) and want to define a set of successively smoothed temporal scale-space representations L(t; τ ) for different values of a temporal scale parameter τ ≥ 0, as schematically illustrated in Figure 1.We will throughout this treatment assume linearity and translational shift covariance, implying that the transformation from the original signal f : R → R to the temporal scale-space representation L : R × R + → R is given by convolution with some one-parameter family of scale-dependent convolution kernels h (1) A crucial condition on this family of temporal scale-space representations is that the temporal scale-space representa-tion L(t; τ 2 ) at any coarser temporal scale t 2 should correspond to a simplification of the temporal scale-space representation L(t; τ 1 ) at any finer temporal scale t 1 .
Following Lindeberg (1990), we shall measure this simplification property in terms of the number of local extrema in the signal at any temporal scale, and define a scale-space kernel as a kernel that obeys the property that the number of local extrema in the signal after convolution is guaranteed to not exceed the number of local extrema prior to the convolution operation, with the important qualifier that this property should hold for any input signal.Equivalently, this property can also be expressed by measuring the number of zero-crossings before and after the convolution operation.A scale-space kernel h(t; τ ) is referred to as a temporal scalespace kernel (Lindeberg and Fagerström 1996) if it additionally satisfies h(t; τ ) = 0 for t < 0, meaning that it does not require access to the future relative to any time moment.
To make the scale simplification property from finer to coarser temporal scales hold, we will assume that the family of temporal smoothing kernels h(u; τ ) should obey the following cascade smoothing property4 for any pair of temporal scales (τ 1 , τ 2 ) with τ 2 > τ 1 and for some family of transformation kernels (∆h)(t; τ 1 → τ 2 ).
We can then obtain a temporal scale-space representation if and only if the transformation kernel (∆h)(t; τ 1 → τ 2 ) between adjacent temporal scale levels t 1 and t 2 is always a temporal scale-space kernel.

Classification of scale-space kernels for continuous signals
A fundamental question with regard to smoothing of temporal signals concerns what convolution kernels satisfy the conditions of being scale-space kernels.

Complete classification of continuous scale-space kernels
Interestingly, the class of one-dimensional scale-space kernels can be completely classified based on classical results by Schoenberg (1930Schoenberg ( , 1946Schoenberg ( , 1947Schoenberg ( , 1948Schoenberg ( , 1950Schoenberg ( , 1953Schoenberg ( , 1988)), see also the excellent monograph by Karlin (1968).Summarizing the treatment in (Lindeberg 1993b, Section 3.5;2016, Section 3.2), a continuous smoothing kernel is a scale-space kernel if and only if it has a bilateral Laplace-Stieltjes transform of the form (Schoenberg 1950) for −c < Re(s) < c and some c > 0, where C = 0, γ ≥ 0, δ and a i are real and ∞ i=1 a 2 i is convergent.
2.1.2Basic classes of primitive scale-space kernels over a continuous signal domain Interpreted over the temporal domain, 5 this result means that there, beyond trivial rescaling and translation, are two main classes of one-dimensional scale-space kernels: convolution with Gaussian kernels convolution with truncated exponential functions for some strictly positive |λ|.
Moreover, the result means that a continuous smoothing kernel is a scale-space kernel if and only if it can be decomposed into a cascaded convolution of these primitives.
2.2 Time-causal temporal scale-space kernels over continuous temporal domain Among the above primitive smoothing kernels, we recognize the Gaussian kernel, which is a good and natural temporal smoothing kernel to use when analysing pre-recorded signals in offline scenarios.When analysing temporal signals in a real-time situation, or when modelling biological processes that operate in real time, we cannot, however, use a temporal smoothing kernel that requires access to information in the future relative to any time moment.
For building a time-causal temporal scale-space representation, the truncated exponential kernels are therefore the only possible primitive time-causal temporal smoothing kernels (Lindeberg and Fagerström 1996) 5 In the general expression (3) for the bilateral Laplace-Stieltjes transform of a continuous scale-space kernel, the factor e γs 2 is the Laplace-Stieltjes transform of the Gaussian kernel e −γξ 2 , the factor 1/(1 + a i s) is the Laplace-Stieltjes transform of a truncated exponential function e −a i ξ /a i with time constant a i , whereas the factors e δs and e a i s correspond to translations in the temporal domain.Furthermore, the general product form of this expression in the Laplace-Stieltjes domain corresponds to a convolution of the corresponding primitives over the original temporal domain.
where we will throughout this treatment adopt the convention of normalizing these kernels to unit L 1 -norm.The Laplace transform of such a kernel is given by Coupling K such kernels in cascade leads to a composed kernel having a Laplace transform of the form The temporal mean and variance of the composed kernel is The temporal mean m K is a coarse measure of the temporal delay of the time-causal temporal scale-space kernel, and the temporal variance τ K is a measure of the temporal duration, also referred to as the temporal scale.
In terms of physical models, repeated convolution with this class of temporal scale-space kernels corresponds to coupling a series of first-order integrators with time constants µ k in cascade with L(t; 0) = f (t), where the temporal scale-space representations for larger values of the scale parameter t k constitute successively temporally smoothed representations of each other.An important property of this type of temporal scale-space representation is that it is also time-recursive.
The temporal scale-space representations L(t; τ k ) constitute a sufficient temporal memory of the past to compute the temporal scale-space representation and the next temporal moment, given a new input in the input signal f (t).
An important consequence of the above necessity result, is that this type of scale-space representation constitutes the only way to compute a time-causal temporal scale-space representation, given the requirement that the number of local extrema, or equivalently the number of zero-crossings, in the signal must not increase from finer to coarser temporal scales.In this respect, the temporal scale-space representations can be seen as gradual simplifications of each other from finer to coarser temporal scales.
Figure 2 shows an illustration of this model in terms of an electric wiring diagram for transforming an input signal f in to an output signal f out using a set of first-order integrators coupled in cascade.Such first-order temporal integration can be used as a straightforward computational model for temporal processing in biological neurons; see also Koch (1999, Chapters 11-12) regarding physical modelling of the information transfer in the dendrites of neurons.

Logarithmic distribution of the temporal scale levels
When implementing this temporal scale-space concept in practice, a set of intermediate temporal scale levels τ k has to be distributed between some minimum and maximum temporal scale levels τ min = τ 1 and τ max = τ K .Then, it is natural to choose these temporal scale levels according to a geometric series, corresponding to a uniform distribution in units of effective temporal scale τ eff = log τ (Lindeberg 1993a).
If we have a free choice of what minimum temporal scale level τ min to use, a natural way of parameterizing these temporal scale levels is by using a distribution parameter c > 1 such that which by equation ( 10) implies that the time constants of the individual first-order integrators should be given by (Lindeberg 2016, Equations ( 19)-( 20)) If the temporal signal is on the other hand given at some minimum temporal scale τ min , corresponding to an a priori given inner temporal scale of the measurement device, we can instead determine 12) such that τ 1 = τ min and add K − 1 temporal scales with µ k according to (14).
Temporal smoothing kernels of this form, combined with temporal differentiation for different orders of differentiation, to obtain ripples of opposite contrast in the resulting temporal receptive fields, have been used for modelling the temporal part of the processing in models for spatio-temporal receptive fields (Lindeberg and Fagerström 1996;Lindeberg 2015Lindeberg , 2016Lindeberg , 2021b) ) and spectro-temporal receptive fields (Lindeberg andFriberg 2015b, 2015a).

Logarithmic memory of the past
When using a logarithmic distribution of the temporal scale levels according to either of these methods, the different levels in the temporal scale-space representation at increasing temporal scales will serve as a logarithmic memory of the past, with qualitative similarity to the mapping of the past onto a logarithmic time axis in the scale-time model by Koenderink (1988).Such a logarithmic memory of the past can also be extended to later stages in a visual, auditory or other form of neural hierarchy.
An alternative type of temporal memory structure can be obtained if the different truncated exponential kernels are applied, not in a cascade as above, but instead in parallel with a single temporal time constant for each temporal memory channel, for µ k = √ τ k , again with a logarithmic distribution of the temporal scale levels τ k .Such a model for temporal memory has been studied by Howard and his co-workers (Howard 2021;Bright et al. 2020).Then, each temporal memory channel is also a simplification of the input signal f (t), and a record of the past with a given temporal delay and temporal duration.Inversion from the temporal memory channels to the input signal is also more straightforward, from the conceptual similarity to a real-valued Laplace transform (Howard et al. 2018;Howard and Hasselmo 2020).The different temporal memory channels are, however, not guaranteed to constitute formal simplifications of each other, as they are for the cascade model.The theoretical framework for time-causal and time-recursive temporal scale-space representations presented earlier in (Lindeberg and Fagerström 1996;Lindeberg 2016) and here can be seen as providing a theoretical foundation for such time-recursive temporal memory models.

Uniform distribution of the temporal scale levels
An alternative approach to distributing the temporal scale levels is to use a uniform distribution of the intermediate temporal scales implying that the time constants in the individual smoothing steps are given by Then, a compact expression can be easily obtained for the composed convolution kernel corresponding to a cascade of K such kernels Such kernels have also been used in memory models (Goldman 2009).The temporal Poisson model studied in more detail in (Lindeberg 1997a) can be seen as the limit case of such a uniform distribution of the temporal scale levels in the time-discrete case, when the difference between adjacent temporal scales tends to zero, a limit case that, however, only exists for discrete temporal signals (Lindeberg and Fagerström 1996), and which also serves as a multi-scale temporal memory of the past (see the illustrations of how the temporal scale-space representation evolves over time and temporal scales in the time-scale diagrams in Figures 3-5 in (Lindeberg 1997a), which demonstrate the temporal memory properties of such a temporal scale-space representation -specifically observe the property that an event that occurs at a certain temporal moment first appears in the temporal scale-space representation at the finest temporal scale, and then moves to gradually coarser temporal scales as time passes by, and is thus also after some short times gradually forgotten at the finer temporal scales, being taken over temporal structures that appear after the initial temporal event).
For constructing temporal memory processes that are to operate over wide ranges of temporal scales, such models based on a uniform sampling of the temporal scale levels do, however, require a larger number of primitive temporal integrators, and thus more hardware or wetware, compared to a temporal memory model based on a logarithmic distribution of the temporal scale levels.
Combined with temporal differentiation of the smoothing kernel, such temporal kernels have been used for modelling the temporal response properties of neurons in the visual system (den Brinker and Roufs 1992) and for computing spatio-temporal image features in computer vision (Rivero-Moreno and Bres 2004;Berg et al. 2014).
For a given value of the temporal scale (the temporal variance) of such time-causal kernels, the temporal delay for a temporal kernel based on a uniform distribution of the temporal scale levels will, however, also be longer than for a temporal kernel constructed from a logarithmic distribution of the intermediate temporal scale levels.Thus, for formulating computational algorithms for expressing time-critical decision processes in computer vision or machine listening, as well as for modelling time-critical decision processes in biological perception or cognition, we argue that a logarithmic distribution of the temporal scale levels should be a much better choice.
For these reasons, we will henceforth in this treatment focus solely on models based on a logarithmic distribution of the temporal scale levels.
3 Time-causal temporal scale-space representations that also obey temporal scale covariance Beyond the task of representing temporal signals at multiple temporal scales, a main requirement on a temporal scale-space representation should also be the notion of temporal scale covariance,6 so as to be able to consistently handle temporal phenomena and events that occur faster or slower in the world.Temporal scale covariance means that if a signal f (t) is subject to a temporal scaling transformation and then processed, here with a temporal convolution kernel T (t ; τ ) that depends on a temporal scale parameter τ , the result should be essentially similar to the result of applying the same type of processing to the original signal and then rescaling the processed original signal (for other types of processes possibly also complemented with some minor modification, such as a correction of the magnitude of the response).For the task of temporal filtering in a temporal scale-space representation, this implies that the temporal scale-space kernel should commute with temporal scaling transformations, as illustrated in the commutative diagram in Figure 3.
Fig. 3 Commutative diagram for temporal receptive field responses under temporal scaling transformations of the temporal domain.Such transformations describe the effect of events occurring slower or faster in the world.(The commutative diagram should be read from the lower left corner to the upper right corner, and means that irrespective of whether the image is first convolved with a temporal smoothing kernel and then subject to temporal scaling transformation, or whether the temporal signal is first subject to a temporal scaling transformation and then convolved with a temporal smoothing kernel, we should get the same result provided that the temporal scale parameters τ and τ are properly matched to the relative temporal scaling factor S between the two temporal patterns.) This algebraic closedness property under temporal scaling transformations will imply that similar temporal phenomena that occur faster or slower in the world will be treated in a conceptually similar manner.Under variations caused by scaling transformations in the input, the output of applying scale-covariant processing to such temporally rescaled data will be mere temporal rescalings of each other, thus without bias to any particular scales, which would otherwise be a severe shortcoming, if the computational model is not well-behaved under temporal scaling transformations.
In this section, we will describe a theory for how to obtain time-causal temporal scale-space representations that also obey such temporal scale covariance, which in turn makes it possible to construct provably scale-invariant temporal representations at higher levels in a temporal processing hierarchy.The way that we will reach this goal is by constructing a limit kernel that is the convolution of an infinite number of truncated exponential kernels in cascade, with specially chosen time constants that correspond to a geometric distribution of the intermediate temporal scale levels.
Unfortunately, there is no known simple compact explicit expression for this limit kernel in the temporal domain, implying that some of the closed-form calculations using the limit kernel may be interpreted as somewhat technical at the first encounter with this function.Once these algebraic transformation properties have been established for the limit kernel, however, this function can be handled and used in a similar way as other standard functions in mathematics.
For practical implementations, the limit kernel can furthermore for the purpose of computing the representation at a single temporal scale often be very well approximated by a moderate finite number of truncated exponential kernels coupled in cascade, usually between 4 and 8 in our implementations of this concept, because of its rapid convergence properties for suitable values of its internal distribution parameter.In turn, for the purpose of computing another temporal scale-space representation at the next coarser temporal scale, applying a single truncated exponential kernel to the nearest finer temporal scale is sufficient.
In this section, we will first define the limit kernel and derive its transformation properties.Then, we will turn to relating and comparing the limit kernel to two other models used for expressing temporal variations over time.

The time-causal limit kernel
Consider the Fourier transform of the composed convolution kernel that we obtain by coupling K truncated exponential kernels in cascade with a logarithmic distribution of the temporal scale levels and thus time constants according to (13) and ( 14) for some c > 1: By formally letting the number of primitive smoothing steps K tend to infinity and renumbering the indices by a shift in terms of one unit, we obtain a limit object of the form (Lindeberg 2016, Equation (38)) By treating this limit kernel as an object by itself, which will be well-defined because of the rapid convergence by the summation of variances according to a geometric series, interesting relations can be expressed between the temporal scale-space representations obtained by convolution with this limit kernel.
3.1.1Self-similar recurrence relation for the time-causal limit kernel over temporal scales Using the limit kernel, an infinite number of discrete temporal scale levels is implicitly defined given the specific choice of one temporal scale τ = τ 0 : . . .
Directly from the definition of the limit kernel, we obtain the following recurrence relation between adjacent temporal scales: and in terms of the Fourier transform:

Behaviour under temporal rescaling transformations
From the Fourier transform of the limit kernel (25), we can observe that for any temporal scaling factor S it holds that Thus, the limit kernel transforms as follows under a scaling transformation of the temporal domain: If we, for a given choice of distribution parameter c, rescale the input signal f by a temporal scaling factor S = 1/c such that t = t/c, it then follows that the scale-space representation of f at temporal scale τ = τ /c  Hence, under a rescaling of the original signal by a temporal scaling factor c, a rescaled copy of the temporal scale-space representation of the original signal can be found at the next lower discrete temporal scale, relative to the temporal scalespace representation of the original signal.

Provable temporal scale covariance
Applied recursively, the above result implies that the temporal scale-space representation obtained by convolution with the limit kernel obeys a closedness property over all temporal scaling transformations t = c j t with temporal rescaling factors S = c j (j ∈ Z) that are integer powers of the distribution parameter c (Lindeberg 2016, Equation ( 47)), thus allowing for perfect scale covariance over the restricted subset of scaling factors S = c j that precisely matches the specific set of discrete temporal scale levels that is defined by a specific choice of the distribution parameter c.Based on this desirable and highly useful property, it is natural to refer to the limit kernel as the scale-covariant time-causal limit kernel (Lindeberg 2016, Section 5).

Qualitative properties
Figure 4 shows graphs of this time-causal limit kernel as well its first-and second-order temporal derivatives for a few values of the distribution parameter c.As can be seen from the graphs, the raw smoothing kernels have a skewed shape, where the temporal delay increases with decreasing values of the distribution parameter c, and with the explicit measures of the skewness γ 1 and kurtosis γ 2 of these kernels increasing as function of the distribution parameter c according to (Lindeberg 2016, Equations ( 130) and ( 131))

Experimental results
Figure 5 shows the result of smoothing two synthetic temporal signals with the time-causal limit kernel for different values of the temporal scale parameter τ .As can be seen from the graphs, the signal is gradually smoothed from finer to coarser temporal scales, here clearly seen in the way that finer-scale structures are suppressed before coarser-scale structures in the left column and that higher frequencies are suppressed before lower frequencies in the right column.In addition, the temporal delay increases from finer to coarser temporal scales, here seen in terms of different temporal offsets regarding the temporal moments at which the temporal peaks occur.When using a comparably large value of the distribution parameter c, as used in this figure, the temporal delay will be comparably low, which is a preferable property when needing to respond fast in a time-critical context.When using lower values of the distribution parameter, the temporal delay at a given temporal scale will be longer, which may be a preferable property if you want to use the temporal scale-space representations as temporal memory buffers, Fig. 6 Illustration of the temporal scale covariance property of the temporal scale-space representation defined from convolutions with the timecausal limit kernel.In the bottom row, the signal in the right column is a rescaling of the signal in the left column by a temporal scaling factor S = 2 (with the temporal rescaling performed relative to the center of the temporal interval).In the top row, the temporal scale-space representations at the matching temporal scale levels Fig. 8 Illustration of the temporal memory property of the time-causal temporal scale-space representation generated by convolutions with the time-causal limit kernel (and as always for this kernel using a logarithmic distribution of the temporal scale levels).The input signal is here a time-causal peak corresponding to a discrete approximation of the time-causal limit kernel for √ τ = 4.The time-causal temporal scalespace representation of this signal has then been computed by convolution with discrete approximations of the time-causal limit kernel for different temporal scales, using a very low value of the distribution parameter c = 8 √ 2 to enable a more clear visualization of the transition properties between different temporal scales.As can be seen from the illustration, the response at the onset of the temporal peak mainly occurs at fine temporal scales, whereas the trace of the peak moves to coarser temporal scales with increasing time, reflecting the temporal memory property of the time-causal scale-space representation.(Horizontal axis: time, Vertical axis: temporal scale index for temporal scale levels between Figure 6 gives an experimental illustration of the temporal scale covariant property of the time-causal limit kernel.Here, a synthetic signal generated from a simulated Wiener process has been rescaled by a temporal rescaling factor S = 2. From these two input signals, temporal scale-space representations have then been computed at the matching temporal scale levels √ τ = 128 and √ τ = 256.Due to the temporal scale covariance property, these temporal scale-space representations are then also related by the same temporal scaling factor S = 2. Figure 7 gives an illustration of the equality between the two different ways of computing the representation in the upper right corner from the signal in the lower left corner in Figure 6, using either a clockwise orientation or a counterclockwise orientation in the corresponding commutative diagram in Figure 3.As can be seen from the visualization, the results are essentially indistinguishable, showing that a good numerical approximation of temporal scale covariance can also be achieved in a discrete implementation (to be described further in Section 4).
Figure 8 shows an illustration of the temporal memory property of this type of temporal scale-space representation.For a very low value of the distribution parameter, c = 8 √ 2 to enable a dense sampling of the temporal scale levels, and thus a more clear visualization of the transitions from finer to coarser temporal scales, we show the response to a temporal peak as function of time and temporal scales in the temporal scale-space representation.As can be seen from the figure, the trace in temporal scale space moves from finer scales at the appearance of the temporal peak, and then to successively coarser levels of temporal scales as time flows.After a certain amount of time, the only trace of the temporal peak, which originated just after t = 0, has moved to only being present at coarser temporal scales.In this way, the temporal scale channels at successively coarser levels of temporal scales serve as a temporal memory of what has happened during different time intervals in the past.

Applications of the time-causal limit kernel
The time-causal limit kernel and its temporal derivatives has been used for modelling the temporal component in spatiotemporal receptive fields in the retina, the LGN and the primary visual cortex (V1) (Lindeberg 2021b), for modelling the temporal component in methods for spatio-temporal feature detection in video data (Lindeberg 2016), for expressing methods for temporal scale selection in temporal signals (Lindeberg 2017(Lindeberg , 2018b)), for modelling the temporal component of spatio-temporal smoothing in methods for spatiotemporal scale selection (Lindeberg 2018a(Lindeberg , 2018b) ) and for modelling the temporal component of smoothing in computer vision methods for video analysis (Jansson and Lindeberg 2018).
In Section 7.3, we do additionally propose to use the time-causal limit kernel for modelling temporal phenomena at multiple temporal scales in neural signals, and in Section 7.2 specifically to use this kernel for modelling the temporal variability in auditory receptive fields.
In Section 6.1 we outline how the time-causal limit kernel can be used for defining time-causal and time-recursive wavelet representations, and in Section 6.2 how the timecausal limit kernel makes it possible to define time-causal and time-recursive time-frequency representations (spectrograms) that additionally obey temporal scale covariance.

Alternative scale-covariant temporal models
An alternative type of temporal model that one could also consider from the general classification of temporal scalespace kernels is to use a set of parallel temporal channels formed by convolution of the input signal, with a single truncated exponential function in each channel, and with a geometric distribution of the their time constants, of the form (16).As previously explained in Section 2.4, such temporal models have been previously used as models of temporal memory in neuroscience (Howard 2021;Bright et al. 2020).
Because of the geometric distribution of the time constants in these temporal channels, they will obey temporal Fig. 9 Comparison between (blue curves) the time-causal limit kernel according to (25) and approximated using the first K = 7 components of the infinite convolution of truncated exponential kernels in cascade with its first-and second-order temporal derivatives and (brown curves) the temporal kernels in Koenderink's scale-time model ( 37) and their first-and second-order temporal derivatives.All kernels correspond to temporal scale (variance) τ = 1 with the additional parameters determined such that the temporal mean values (the first-order temporal moments) become equal in the limit case when the number of temporal scale levels K tends to infinity (equation ( 38  scale covariance.Temporal scale covariance will also apply to different types of generalizations of such a model, e.g. by having the same small number of truncated exponential kernels in cascade in each temporal channel, with the time constants between the different temporal channels coupled according to a geometric distribution.
A fundamental difference between such temporal models and the temporal scale-space model based on the timecausal limit kernel, however, is that in the first class of models the temporal channels for larger values of the scale parameter are not guaranteed to constitute simplifications of the temporal channels for smaller values of the scale parameter.By the temporal smoothing kernels being scale-space kernels, each temporal channel is guaranteed to be a simplification of the input signal.When relating different temporal scale channels to each other, however, the number of local extrema in a temporal channel for a larger value of the temporal scale parameter is not guaranteed to not exceed the number of local extrema in a temporal channel for a finer smaller value of the temporal scale parameter.
Because of the scale-recursive property (28) of the timecausal limit kernel, it is on the other hand formally guaranteed that the temporal scale-space representation at the next coarser temporal scale corresponds to the result of applying temporal smoothing with a truncated exponential kernel to the temporal scale-space representation at the nearest finer temporal scale.Applied recursively, the temporal scale-space representation at any coarser temporal scale corresponds to the result of applying a set of truncated exponential kernels in cascade to the representation at any finer temporal scale.In this way, for the temporal scale-space representation generated by convolution with the time-causal limit kernel for different values of the temporal scale parameter, every temporal scale-space representation at a given temporal scale is guaranteed to constitute a formal simplification of any other temporal scale-space representation at any finer temporal scale.
The time-causal limit kernel is special in that it both obeys temporal scale covariance and guarantees non-creation of new local extrema with increasing temporal scales with regard to convolutions over a time-causal temporal domain.

Relation to Koenderink's scale-time model
In his scale-time model, Koenderink (1988) proposed to perform a logarithmic mapping of the past via a temporal delay δ and then applying Gaussian smoothing with standard deviation σ in the transformed domain.If we additionally normalize these kernels to unit L 1 -norm, we obtain a timecausal kernel of the form (Lindeberg 2016, Equation ( 151)) In (Lindeberg 2016, Appendix 2) a formal mapping between this scale-time kernel and the time-causal limit kernel is derived, by requiring the first-and second-order moments of these two classes of kernels to be equal: which hold as long as c > 1 and σ < √ log 2 ≈ 0.832.Figure 9 shows a comparison between the time-causal limit kernel and Koenderink's scale-time kernels regarding the zero-order convolution kernels as well as their first-and second-order derivatives.As can be seen from the graphs, these two classes of kernels have qualitatively rather similar shapes.The time-causal limit kernel does, however, have the conceptual advantage that it can be computed in a timerecursive manner, whereas the scale-time kernel does not have any known time-recursive implementation, implying that it formally requires an infinite memory of the past (or some substantially extended temporal buffer, if the infinite temporal convolution integral is truncated at the tail).
While we do not have any compact explicit expression for the time-causal limit kernel over the temporal domain, if we approximate the time-causal limit kernel by a scaletime kernel according to the mapping (38), we obtain the following estimate for the location of the maximum point of the time-causal limit kernel: This estimate can be expected to be an overestimate, and is a better estimate of the temporal delay of the time-causal limit kernel than the temporal mean according to (10).
3.4 Relation to the ex-Gaussian model used by Bright et al.
In (Bright et al. 2020), a so-called ex-Gaussian model (Grushka 1972), that is the convolution of an unnormalized Gaussian function with an unnormalized truncated exponential kernel is used for fitting temporal response functions of neurons to an analytical temporal model.In Appendix A.1, a relation between this ex-Gaussian model and a corresponding model based on the time-causal limit kernel is derived by requiring the zero-, first-and second-order temporal moments of these kernels to be equal, if the DC-offsets a 0 and b 0 are disregarded and assumed to be equal.This leads to the following mapping between the parameters of the two models where δ and V denote the temporal mean and the temporal variance of the ex-Gaussian model for a 0 = 0 and M 0 , M 1 and M 2 denote the explicit expressions for the zero-, first-and second-order moments of the ex-Gaussian model for a 0 = 0, according to (88), ( 89) and (90).
Figure 10 shows a few examples of ex-Gaussian temporal models approximated by models based on the timecausal limit kernel in this way.As can be seen from the graphs, the two classes of kernels can capture qualitative similar temporal shapes in time-causal temporal data, 7 with the conceptual differences that: (i) the model based on the time-causal limit kernel always tends to zero at the temporal origin t = 0 when the DC-offset is zero, whereas the ex-Gaussian model may take non-zero values for t = 0, (ii) the time-causal limit kernel does not contain any internal non-causal temporal component as the time-shifted Gaussian kernel in (40) constitutes, and (iii) the time-causal limit kernel has a completely time-recursive implementation, which is essential when modelling temporal phenomena in real time as they, for example, occur in biological neurons.The model based on the time-causal limit kernel is also specifically possible to implement based on a cascade of first-order integrators in cascade, which is a natural model for the information transfer in the dendrites of neurons (Koch 1999, Chapters 11-12).

Extension to third-order moment-based model fitting involving also a flexible temporal offset
In Appendix A.2, an extension of the above second-order moment-based model to a third-order moment-based model is performed, which makes it possible to also determine a temporal offset t 0 and which may be relevant in situations when the temporal origin of the signal cannot be accurately determined in an experimental situation.Since the closed-form expressions for the solutions become more complex in this case (they are determined from the solutions of a fourth-order algebraic equation), we restrict ourselves to a conceptual and algorithmic description in this treatment, see Appendix A.2 for further theoretical details and experimental results.

Extension to model fitting for other signals or functions
The above general procedures, whereby the parameters in the model based on the time-causal limit kernel are determined from the lower-order temporal moments of the data, can also be more generally used for fitting models based on the time-causal limit kernel to other signals and functions that: (i) are defined for non-negative values of time, (ii) assume non-negative values only, (iii) have a roughly unimodal shape of first increasing and then decreasing and (iv) decay towards zero towards infinity.The approach for fitting basically implies replacing the temporal moments M 0 , M 1 , M 2 and optionally M 3 of the ex-Gaussian model by the temporal moments of the signal or function to be fit with a model based on the time-causal limit kernel, see Appendix A.3 for additional details.
4 Computational implementation of convolutions with the time-causal limit kernel on discrete temporal data In the theory presented so far, we have throughout assumed that the signal is continuous over time.When implementing this model on sampled temporal data, the theory must be transferred to a discrete temporal domain.
In this section, we will describe how the temporal receptive fields can be implemented in terms of corresponding discrete temporal scale-space kernels that possess scalespace properties over a discrete temporal domain, and in addition are both time-causal and fully time-recursive.
Following Lindeberg (1990) and in a corresponding way as the treatment in Section 2, let us define a discrete kernel as a discrete scale-space kernel if for any input signal it is guaranteed that the number of local extrema, alternatively the number of zero-crossings, cannot increase under convolution with the discrete scale-space kernel.

Classification of scale-space kernels for discrete signals
To characterize the class of discrete scale-space kernels, we can, in a corresponding way as for the continuous case, also build upon classical results by Schoenberg (1930Schoenberg ( , 1946Schoenberg ( , 1947Schoenberg ( , 1948Schoenberg ( , 1950Schoenberg ( , 1953Schoenberg ( , 1988)), and as further developed in the monograph by Karlin (1968).
Making a summary of the treatment in Lindeberg (1990, Section IV) (2016, Section 6.1), a discrete smoothing kernel is a discrete scale-space kernel if and only if it has its generating function of the sequence of filter coefficients ϕ(z) = ∞ n=−∞ c n z n of the form (Schoenberg 1948)

Basic classes of primitive scale-space kernels over a discrete signal domain
With regard to the original temporal domain,8 this characterization means that, besides trivial rescalings and translations, there are three basic classes of discrete smoothing transformations: two-point weighted average or generalized binomial smoothing moving average or first-order recursive filtering infinitesimal smoothing 9 or diffusion as arising from the continuous semi-groups made possible by the factor e (q−1z −1 +q1z) .
To transfer the continuous first-order integrators derived in Section 2.2 to a discrete implementation, we shall in this treatment focus on the first-order recursive filters (50), which by additional l 1 -normalization constitute both the discrete correspondence and a numerical approximation of time-causal and time-recursive first-order temporal integration (11).

Discrete temporal scale-space kernels based on recursive filters
Given a signal that has been sampled by some temporal frame rate r, the temporal scale σ t in the continuous model in units of seconds is first transformed to a temporal variance τ relative to a unit time sampling Then, a discrete set of intermediate temporal scale levels τ k is defined by ( 12) or ( 17), with the difference between successive scale levels according to with τ 0 = 0.
For implementing the temporal smoothing operation between two such adjacent scale levels (with the lower level in each pair of adjacent scales referred to as f in and the upper level as f out ), we make use of a first-order recursive filter normalized to the form and having a generating function of the form 9 These kernels correspond to infinitely divisible distributions as can be described with the theory of Lévy processes (Sato 1999), where specifically the case q −1 = q 1 corresponds to convolution with the non-causal discrete analogue of the Gaussian kernel (Lindeberg 1990) and the case q −1 = 0 corresponds to convolution with time-causal temporal Poisson kernel (Lindeberg and Fagerström 1996;Lindeberg 1997a).
which is a time-causal kernel and satisfies discrete scalespace properties of guaranteeing that the number of local extrema or zero-crossings in the signal will not increase with increasing scale (Lindeberg 1990;Lindeberg and Fagerström 1996).These recursive filters are the discrete analogue of the continuous first-order integrators (11).
Each primitive recursive filter (53) has temporal mean value m k = µ k and temporal variance ∆τ k = µ 2 k + µ k , and we compute µ k from ∆τ k in (52) according to By the additive property of variances under convolution, the discrete variances of the discrete temporal scale-space kernels will perfectly match those of the continuous model, whereas the temporal mean values and the temporal delays may differ somewhat.If the temporal scale τ k is large relative to the temporal sampling distance, the discrete model should be a good approximation in this respect.
By the time-recursive formulation of this temporal scalespace concept, the computations can be performed based on a compact temporal buffer over time, which contains the temporal scale-space representations at temporal scales τ k , and with no need for storing any additional temporal buffer of what has occurred in the past, to perform the corresponding temporal smoothing operations.
For practical implementations, we often approximate the time-causal limit kernel using 4 to 8 layers of recursive filters coupled in cascade using either c = √ 2 or c = 2.A summarizing algorithmic description of how to implement these temporal filtering operations in practice is given in Appendix B.

Computation of temporal scale-space derivatives
So far, we have been concerned with the problem of how to smooth a temporal signal in such a way that the smoothing transformation is guaranteed to not increase the number of local extrema in the signal, or equivalently the number of zero-crossings.In many applications, one is, however, more interested in studying the change in the signal over time, as can be modelled by temporal derivatives.
For a purely time-dependent signal, the first-order temporal derivative will lead to strong responses in the signal when the temporal slope is high, corresponding to e.g.onsets or offsets of a sound in auditory processing, or motion in the world, alternatively changes in the illumination, for video processing.Regarding visual processing over a purely spatial domain, first-order spatial derivatives will respond to edges in the image domain, which in turn may correspond to discontinuities in either depth, surface orientation, reflectance or illumination in the world.
For a purely time-dependent signal, the second-order derivatives may on the other hand often lead to strong responses near local maxima or minima over time, if the sign of the first-order temporal derivative changes rapidly at those points.Concerning audio processing, a second-order temporal derivative applied to a spectrogram representation may give a strong response to e.g. a beep or some other brief temporal sound, provided that the temporal scale is sufficiently near the temporal duration of the sound.Applying second-order derivatives with respect to logarithmic frequencies to a spectrogram will, in turn, enhance spectral bands and formants, provided that the logspectral scales are appropriately selected.Regarding visual processing, a second-order temporal derivative applied to a video stream may give a strong response to a flashing light, again assuming that the temporal scale is sufficiently near the temporal duration of the flash.Assuming that the visual observer does not fixate a moving object, second-order temporal derivatives may also give strong responses to image patterns that move relative to the viewing direction.For visual processing on a purely spatial domain, second-order spatial derivative operators can be specially designed to give strong responses to blob-like or corner-like image structures, which can be detected by interest point detectors.
Beyond such pointwise or regionwise responses over time, as described above, temporal derivatives can also be interpreted and used densely, for every time moment, and, for example, be combined according to a local Taylor expansion around any temporal moment t 0 : to characterize the local temporal structures in the temporal signal at any scale τ .Such a representation involving temporal derivatives up to order N is referred to as a temporal N -jet representation.10A practical complication that, however, arises, when computing temporal derivatives at multiple scales concerns how to compare the responses between different levels of scale.Due to the temporal smoothing operation, the amplitude of the temporal derivatives can be expected to decrease monotonically with increasing amount of temporal smoothing, provided that the temporal smoothing operation is sufficiently well-designed.This does, for example, hold for temporal smoothing with the truncated exponential kernels, which arise as the only possible temporal smoothing primitives in the time-causal scale-space kernels, including the time-causal limit kernel.
In this section, we will describe a way to reduce the problem of decreasing amplitude of temporal derivatives with increasing values of the temporal scale parameter, by instead using scale-normalized temporal derivatives.The intention is that by using appropriately designed scale-normalized derivative operators, it should be possible to judge if a temporal derivative response of a certain order at a certain temporal scale should be regarded as stronger or weaker than a corresponding temporal derivative response at some other temporal scale.We will also describe how temporal scale covariance can be obtained for temporal derivative operators that are combined with the time-causal limit kernel.

The scale-normalized derivative concept
For the non-causal Gaussian scale-space concept defined over a purely spatial domain, and corresponding to Gaussian smoothing at all scales, it can be shown that the canonical way of defining scale-normalized derivatives at different spatial scales s is according to (Lindeberg 1998a(Lindeberg , 1998b(Lindeberg , 2021a) where γ is a free parameter.Specifically, it can be shown (Lindeberg 1998a, Section 9.1) that this notion of γ-normalized derivatives corresponds to normalizing the m:th order Gaussian derivatives g ξ m over N -dimensional image space to constant L p -norms over scale with the power p in the L p -norm depending on the scale normalization power γ, the order of differentiation m and the spatial dimensionality N of the signal according to where the perfectly scale-invariant case γ = 1 corresponds to L 1 -normalization for all orders m.

Scale normalization for time-causal temporal derivatives
For temporal derivatives11 defined from the time-causal scalespace concept corresponding to convolution with truncated exponential kernels coupled in cascade, it can be shown to be meaningful to define time-causal scale-space derivatives in a corresponding manner (Lindeberg 2016(Lindeberg , 2017)): -By variance-based scale normalization, we define scalenormalized temporal derivatives according to where τ denotes the variance of the temporal smoothing kernel.
-By L p -norm-based scale normalization, we determine a temporal scale normalization factor α n,γ (τ ) such that the L p -norm (with p determined as function of γ according to ( 59)) of the corresponding composed scale-normalized temporal derivative computation kernel α n,γ (τ ) h t n equals the L p -norm of some other reference kernel, where we may initially take the L p -norm of the corresponding Gaussian derivative kernels (Lindeberg 2016, Section 7.3)

Scale covariance property of scale-normalized temporal derivatives
In the special case when the temporal scale-space representation is defined by convolution with the scale-covariant time-causal limit kernel according to ( 26) and ( 25), it is shown in (Lindeberg 2016, Appendix 3) that the corresponding scale-normalized derivatives become truly scale covariant under temporal scaling transformations t = c j t with scaling factors S = c j that are integer powers of the distribution parameter c between matching temporal scale levels τ = c 2j τ .Specifically, for γ = 1 corresponding to p = 1 the magnitude values of the scale-normalized temporal derivatives at matching scales become fully scale invariant allowing for well-defined comparisons between the magnitude values of different types of temporal structures in a signal at different temporal scales.

A canonical class of time-causal, time-recursive and scale-covariant temporal basis functions
The above scale covariance property implies that the scalenormalized temporal derivatives of the time-causal limit kernel constitute a canonical class of temporal basis functions over a time-causal temporal domain.
These kernels have been used as temporal basis functions for spatio-temporal receptive fields (Lindeberg 2016(Lindeberg , 2021b;;Jansson and Lindeberg 2018) and for expressing methods for temporal scale selection (Lindeberg 2017(Lindeberg , 2018b) ) and spatio-temporal scale selection (Lindeberg 2018a(Lindeberg , 2018b) that detect and compare temporal structures at different temporal scales in a completely scale-invariant manner.
In this treatment, we additionally propose to use this family of temporal basis functions to model the temporal variability of neurons over multiple scales (Section 7.3) and specifically the temporal variability in computational models of auditory receptive fields (Section 7.2).

Discrete approximations of scale-normalized temporal scale-space derivatives
For the discrete temporal scale-space concept over discrete time described in Section 4.2, discrete approximations of temporal derivatives are obtained by applying temporal difference operators to the discrete temporal scale-space representation at any temporal scale, which in turn is constructed from a cascade of first-order recursive filters of the form (53), with the time constants µ k given by (55) from the differences in temporal scale levels ∆τ k = τ k − τ k−1 with τ k according to (12).Scale normalization factors for discrete l p -normalization are then defined in an analogous way as for continuous signals, (60) or (61), with the only difference that the continuous L p -norm is replaced by a discrete l p -norm.

Experimental results
Figure 11 shows an illustration of computing discrete approximations of second-order scale-normalized temporal derivatives in this way, 12 for a synthetic input signal consisting of two temporal peaks generated from discrete approximations of the time-causal limit kernel for τ = 16 and τ = 256, respectively, and with some amount of relative temporal delay to separate the responses as well as a small amount of added white Gaussian noise.
Observe how the dominant responses to the finer-scale structures in the input signal are obtained at finer levels of scale in the temporal scale-space representation, whereas the dominant responses to the coarser-scale structures in the input signal are obtained at coarser levels of scale in the temporal scale-space. 12Here, using distribution parameter c = 2 for the time-causal limit kernel and scale normalization power γ = 1 for the scale-normalized temporal derivative operator.Do also observe how the responses at coarser temporal scales are associated with longer temporal delays, manifesting themselves as temporal peaks corresponding to the underlying signal structures appearing at later time moments at coarser levels of scale.
Do furthermore note that the range of values on the vertical axis in these graphs is the same for all the scale values, demonstrating the ability to make relative comparisons between the magnitudes of the derivative responses at different scales, due to the notion of scale normalization of the temporal derivatives, here with regard to the l 1 -norm.

Relations to wavelet analysis and time-frequency analysis
For analyzing temporal signals at multiple temporal scales, wavelet analysis (Grossmann and Morlet 1984;Mallat 1989Mallat , 1999;;Heil and Walnut 1989;Meyer 1992;Daubechies 1992;Chui 1992;Rioul and Duhamel 1992;Graps 1995;Debnath and Shah 2002) and time-frequency analysis (Gabor 1946;Cohen 1995;Feichtinger and Strohmer 1998;Qian and Chen 1999;Gröchenig 2001;Flandrin 2018) constitute two other main classes of conceptual tools.In this treatment, we do, however, not follow those notions as prototype models, instead adhering to the scale-space paradigm because of its special properties.Nevertheless, the presented temporal scale-space theory can be related to wavelet analysis and time-frequency analysis in the following ways:

Relations to wavelet analysis
By construction, the temporal derivatives of the time-causal limit kernel Ψ (t; τ, c) defined from (25) have integral equal to zero In this respect, the temporal derivatives of the time-causal kernel, complemented by normalization with respect to a suitably chosen norm, can serve 13 as a mother wavelet over a continuous time-causal temporal domain, in a similar way as Gaussian derivative kernels of a certain order 13 Additionally, one usually states a requirement that the wavelet function should decrease sufficiently fast at the tails, such that for some α > 0. As will be shown later, for the temporal derivatives of the time-causal limit kernel, fulfilment of this condition follows from an exponential decrease towards zero at the infinity, see Equation ( 75).
(68) such as the Mexican hat wavelet (Marr 1976(Marr , 1982)), also known as a Ricker wavelet (Ricker 1944;Hosken 1988), and corresponding to the second-order derivative of the Gaussian, can serve as a mother wavelet over a continuous noncausal temporal domain.
In wavelet analysis, one usually normalizes both the mother wavelet and the child wavelets to unit L 2 -norm, leading to translated and rescaled child wavelets of the form In scale-space theory, the most common way of normalizing the Gaussian derivative kernels as well as temporal derivatives of the time-causal limit kernel is to constant L 1norm over scales (and corresponding to scale-normalized derivatives for γ = 1 according to Section 5.1), although other scale normalizations, including L 2 -normalization, are also possible, as further described in Section 5.1.Such L 1normalization then leads to translated and rescaled child wavelets of the form In the following, we will describe how the corresponding wavelet representations obtained my mapping a signal f onto the child wavelets can be computed if the mother wavelet is chosen as a temporal derivative of the time-causal limit kernel.

Handling the transformation properties of the child wavelets within the algebra of the time-causal temporal scale-space representation
By using the transformation properties of scale-normalized derivatives of the time-causal scale-space representation of the time-causal limit kernel (63), it follows that under a scaling transformation of time t = c j t for some integer j with c being the distribution parameter of the time-causal limit kernel, and with a corresponding transformation of the temporal scale parameter τ = c 2j τ , similar transformation properties hold for the scale-normalized temporal derivatives of the time-causal limit kernel (let the input signal be the continuous delta function f (t) = δ(t) in ( 63)) where γ is the power in the temporal scale-normalized derivative concept and p is the power in the corresponding L pnorm that is kept constant over scale by the scale-normalized derivatives.This implies that if we choose the mother wavelet as a temporal derivative of the time-causal limit kernel according to (67), then the temporal scaling and translation operations of the child wavelets in ( 69) and ( 70) can be expressed fully within the algebra of the time-causal scale-space representation, provided that the temporal scaling factors a are chosen as integer powers of the distribution parameter c in the time-causal limit kernel according to a = c j .This does in turn imply that the result of expanding a temporal test signal onto the child wavelets can be directly extracted as the corresponding temporal derivatives of the time-causal temporal scale-space representation of the temporal test signal at the different temporal scales, possibly complemented by a scale-dependent scaling of the magnitude values, depending on the choice of L p -norm in the wavelet representation and the choice of scale normalization power γ in the scalenormalized derivative concept.

Finite L p -norms for the temporal derivatives of the time-causal limit kernel
A regularity requirement that one usually imposes on wavelet functions is that they should be in both L 1 (R) and L 2 (R).This property can be easily shown for the temporal derivatives of the time-causal limit kernel, as follows: Consider a partial fraction decomposition of the Laplace transform (9) of the infinite convolution of truncated exponential kernels that defines the time-causal limit kernel according to (25): with µ k as functions of τ and c according to ( 13) and ( 14), and where the coefficients A k can be determined by first multiplying both sides of the equation by (1 + µ k q) and then setting q = −1/µ k , leading to Interpreted over the original temporal domain, this means that the time-causal limit kernel can be written in terms of the following decomposition as a sum of truncated exponential functions: Thus, the n:th order temporal derivative of the time-causal limit kernel will have the following series representation: When time t tends to infinity, this function will in the limit tend towards zero, and as fast as exponentially with respect to he slowest time constant µ 1 .Since (∂ t n Ψ )(t; τ, c) is additionally finite for finite values of t, it follows that both the L 1 -and the L 2 -norms of ∂ t n Ψ will be finite, implying that , thus proving the result.

Time-causal and time-recursive wavelets for real-time and time-critical applications
These resulting wavelets described in this section, consisting of temporal derivatives of the time-causal limit kernel, will be completely time-causal.The convolutions14 between these wavelet kernels and a temporal measurement function can also be computed in a completely time-recursive way, thus eliminating the need for additional temporal buffering and in turn allowing for minimal temporal response times in a time-critical context.In these respects, the temporal derivatives of the time-causal limit kernel may thus have interesting potential use for wavelet analysis with regard to applications that are to be performed over time-causal and time-recursive temporal domains, such as for real-time signal analysis systems, or when modelling physical or biological systems for which access to the relative future in relation to any time moment is not possible.Another type of time-causal wavelet representation has been proposed and studied by Szu et al. (1992), based on linear combinations of sine and cosine waves multiplied by a truncated exponential function.In this context, the wavelets based on temporal derivatives of the time-causal limit kernel have the conceptual advantage that they are solely based on truncated exponential kernels coupled in cascade, and can therefore be implemented in a fully time-recursive manner. 15dditionally, with regard to the discrete implementation of such temporal receptive fields in terms of recursive filters coupled in cascade (according to Section 4.2), the computation of wavelets based on temporal derivatives of the timecausal limit kernel, an additional temporal scale level can be computed with just the addition of a single recursive filter, complemented with a discrete temporal difference operator (according to Section 5.5).

Relations to time-frequency analysis
If we combine the time-causal limit kernel Ψ (t; τ, c) defined according to (25) with pointwise multiplication by a complex exponential function e iωt , then we obtain a straightforward way of defining a time-causal time-frequency representation of a temporal signal f (t) according to where the complex-valued extension of the time-causal limit kernel χ(t, ω; τ, c) = Ψ (t; τ, c) e iωt (77) can be seen as a time-causal analogue16 of the Gabor function (Gabor 1946), with the role of the Gaussian kernel g(t; σ) in the Gabor function now replaced the by the time-causal limit kernel Ψ (t; τ, c) for τ = σ 2 .Figure 12 shows graphs of a few examples of such complex-valued extensions of the time-causal limit kernel for different values of the angular frequency ω in relation to a given temporal scale τ .In this context, the time-causal limit kernel serves as a temporal window function for computing a windowed Fourier transform, to give better localization properties in the temporal domain compared to a regular Fourier transform, and where the window function in this case, in contrast to the more common choice of a Gaussian window function, is fully time-causal, to allow for real-time processing as well as realistic modelling of real-world physical and biological processes, where access to the relative future in relation to any time moment is simply not possible.

Relations to the Gammatone filter
The complex-valued extension of the time-causal limit kernel in (77) is specifically closely related to the Gammatone filter (Johannesma 1972;Patterson et al. 1987Patterson et al. , 1995;;Hewitt and Meddis 1994) in auditory processing with the main difference being that the truncated exponential kernels used in this auditory filter have equal time constants, and can thus under a convolution operation be composed into a single monomial multiplied by the complex exponential, in analogy with Equation ( 19), and thereby corresponding to a uniform distribution of the temporal scale levels according to Section 2.5, whereas the temporal scale levels in the complex-valued extension of the time-causal limit kernel are constructed according to a geometric distribution of the temporal scale levels according to Section 2.3, thus, in turn, allowing for different and more rapid temporal dynamics.Another minor difference is that the phase of the Gammatone filter is represented as a phase angle α of a cosine function, whereas the phase of the complex-valued extension of the time-causal limit kernel is represented as the phase value of a complex exponential.

Relations to the Heisenberg group
The time-frequency representation defined according to (76) has the theoretically attractive property that it is closed under (i) translations over time, (ii) multiplicative shifts in the frequency of periodic or repetitive temporal signals and (iii) uniform scaling transformations of the temporal axis with discrete scaling factors S that are integer powers of the distribution parameter c.Hence, except for the necessary discretization of the temporal scale parameter according to a geometric distribution, which implies closedness over a discrete set of scaling factors as opposed to as over a continuum, this time-frequency representation has the ability to capture similar types of transformations of the signal as the Gabor family, and as can be modelled by the Heisenberg group, see (Feichtinger and Gröchenig 1992).In this way, the complex-valued time-causal limit kernel provides a way to define a scale-covariant time-frequency representation also over a time-causal temporal domain.

Extension to an additionally time-recursive time-frequency transform
If one additionally wants these time-frequency representations to also be time recursive, then it is possible to modify this construction slightly, by instead multiplying the input signal by a set of complex exponentials and then filtering the resulting complex-valued signal with the time-causal limit kernel (according to Equation ( 81)), thus implying that this time-frequency transform can be implemented discretely in terms of a set of recursive filters that operate over time on the pointwise multiplication of the input signal with a set of complex exponential functions.The difference will then be that the phase values will have to be compensated a posteriori, whereas the magnitude values of the corresponding spectrogram will be preserved.An earlier version17 of this type of theoretical model has been successfully used for computing auditory receptive fields (Lindeberg andFriberg 2015a, 2015b), as will be further described in Section 7.2.

Applications to modelling temporal variations in biological systems
In this section, we will describe different application domains of using the theory for temporal scale-space repre-sentation, specifically the time-causal limit kernel, to model temporal variations in biological signals.

Temporal basis functions in spatio-temporal receptive field models
In (Lindeberg 2011(Lindeberg , 2013a), a general model for spatio-temporal receptive fields is derived of the form where x = (x 1 , x 2 ) T denotes the image coordinates, t denotes time, s denotes the spatial scale, τ denotes the temporal scale, v = (v 1 , v 2 ) T denotes a local image velocity, -Σ denotes a spatial covariance matrix determining the spatial shape of an affine Gaussian kernel g(x; s, Σ) = in space-time and h(t; τ ) is a temporal smoothing kernel over time.
This model for zero-order spatio-temporal receptive fields should, in turn, be complemented by spatial and temporal differentiation to lead to spatio-temporal receptive fields with positive and negative lobes that are balanced in the sense of the integral of the filter weights being equal to zero.
In (Lindeberg 2016(Lindeberg , 2018a(Lindeberg , 2021b)), it is described how the time-causal limit kernel can be successfully be used as the temporal smoothing kernel in this context, i.e., h(t; τ ) = Ψ (t; τ ) with Ψ defined from its Fourier transform according to (25), and allowing for truly time-causal and timerecursive model of spatio-temporal receptive fields, which in turn enable provable scale covariance and scale invariance properties over the temporal domain.
By comparisons with biological visual receptive fields measured by electrophysiological cell recordings by DeAngelis et al. (1995DeAngelis et al. ( , 2004)), it is shown in (Lindeberg 2016(Lindeberg , 2021b) that this spatio-temporal receptive field model very well captures the qualitative shape of lagged and non-lagged LGN neurons as well as simple cells in the primary visual cortex (V1).

Temporal basis functions in spectro-temporal receptive field models
In (Lindeberg andFriberg 2015a, 2015b), a theoretical framework for idealized models of auditory receptive fields is presented, based on a two-stage model consisting of time-causal spectrograms followed by time-causal spectro-temporal receptive fields applied on these, and which comprises covariance and invariance properties under natural sound transformations, such as frequency shifts and glissando transformations.
The time-causal spectrograms in this model are defined according to where the temporal integration kernel h composed is from theoretical arguments constrained to be the convolution of a set of truncated exponential kernels coupled in cascade.Following the arguments in this paper, and further restricting this kernel to be a time-causal limit kernel Ψ , we can extend the previous theoretical framework for multi-scale spectrograms to also comprise temporal scale covariance.
In the second-stage model of spectro-temporal receptive fields in this theory, the idealized form of auditory receptive fields are from theoretical arguments constrained to be of the form where -∂ t α represents a temporal derivative operator of order α with respect to time t which could alternatively be replaced by a glissando-adapted temporal derivative of the form ∂ t = ∂ t + v ∂ ν , -∂ ν β represents a logspectral derivative operator of order β with respect to logarithmic frequency ν, -T (t; τ a ) represents a temporal smoothing kernel with temporal scale parameter τ a , which should in the timecausal case be a set of truncated exponential kernels coupled in cascade, g(ν − vt; s) represents a Gaussian spectral smoothing kernel over logarithmic frequencies v with logspectral scale parameter s and v representing a glissando parameter making it possible to adapt the receptive fields to variations in frequency ν = ν + vt over time.By comparison with biological auditory receptive fields measured by electrophysiological cell recordings by Qiu et al. (2003), Andoni et al. (2007), Machens et al. (2004), Elhilali et al. (2007) and Atencio and Schreiner (2012), it is shown in (Lindeberg and Friberg 2015a) that the idealized receptive fields from this model agree qualitatively very well with biological auditory receptive fields measured in the inferior colliculus (ICC) and primary auditory cortex (A1) of mammals.
By following the arguments regarding temporal smoothing in this paper, and constraining the temporal kernel in the above model to be a time-causal limit kernel, T (t; τ a ) = Ψ (t; τ a ), it follows that the auditory covariance properties in the spectro-temporal receptive field model can be extended to also comprise temporal scale covariance.

Temporal scales in neural signals
In this section, we describe previous evidence and use of multiple temporal scales in neural signals, with relations to the theory for processing temporal signals at multiple scales presented in this paper.
Concerning the use of multiple temporal scales for processing neural signals, Goldman (2009) shows how neural responses can be maintained by a purely feedforward mechanism, which thus implements a temporal memory.In his model, a set of first-order integrators with equal time constants is used.By instead using different time constants of the first-order integrators, as used for the implementation of the time-causal limit kernel, we can get a more compact model for the memory buffers, requiring less wetware or computational modules, with the additional benefit that the time constants obey a self-similar logarithmic distribution.Tsao et al. (2018) show how temporal information in the lateral entorhinal cortex is robustly encoded over a wide range of temporal scales, from time scales of seconds to hours, where specifically the brain handles multiple scales in parallel, consistent with the underlying construction of a multi-scale representation over the temporal domain, and specifically using a multi-scale temporal representation as a temporal memory.In a further study of the primate entorhinal cortex, Bright et al. (2020) experimentally model time cells in this brain area as single truncated exponentials, in line with theoretical model in Equation ( 16), although also complemented with a Gaussian smoothing step that leads to the ex-Gaussian model, and conclude that the time cells in the entorhinal cortex use a spectrum of time constants to construct a temporal record of the past in support of episodic memory.In a study of cerebellar unipolar brush cells, Guo et al. (2021) show that the population of neurons generates a continuum of multi-scale temporal representations, with essentially a logarithmic distribution of the temporal scale levels, consistent with the distribution of temporal scale levels used for the temporal scale-space representation and its associated temporal memory model based on the time-causal limit kernel.
In their computational model, of temporal memory, Howard and Hasselmo (2020) propose that time cells in the hippocampus can be understood as a compressed estimate of events as a function of the past, and that temporal context cells in the entorhinal cortex can be understood as the (real-valued) Laplace transform of that function, respectively, where the Laplace transform in turn arises from the integration with truncated exponential kernels with different time constants, as are used as the unique primitive time-causal temporal smoothing kernel that are guaranteed to not increase the number of local extrema or zero-crossings in the signal.Howard (2021) gives a more general overview of mechanisms for temporal memory, including the use of multiple first-order temporal integrators as arising from this theory.
In an fMRI study of memory recall in human subjects over large variations in the time elapsed after the event, Monsa et al. (2020) conclude that scale-selective activity characterizes autobiographical memory processing and may provide a basis for understanding how the human brain processes and integrates experiences across temporal scales in a hierarchical manner.
Holcombe ( 2009) gives a general overview of different temporal scale limits in visual perception, in particular describing a distinction into slow and fast temporal processes, which are hypothesized to originate from neural processes over different ranges of temporal scales.In an fMRI study of the human ventral stream, Gauthier et al. (2012) show that the widths of temporal integration windows increase at higher hierarchical levels in the visual hierarchy.
Regarding the use of multiple temporal scales in auditory perception, Atencio and Schreiner (2012) show examples of spectro-temporal receptive fields in the primary auditory cortex (A1) with different spectro-temporal scale characteristics; broadly tuned receptive fields with short temporal duration and narrowly tuned receptive fields with longer temporal duration.Chait et al. (2015) investigate how different temporal scales interact in speech perception and suggest that human speech perception uses multi-time resolution processing.Teng et al. (2016) provide evidence that the auditory system extracts fine-detail acoustic information using short temporal windows and uses long temporal windows to abstract global acoustic patterns.Concerning the specific area of birdsong, Gentner (2008) shows how the use of multiple temporal scales within the acoustic pattern hierarchy conveys information about the individual identity of the singer.Osman et al. (2018) also propose a hierarchy of temporal scales for discriminating and classifying the temporal shapes of sound in different auditory cortical areas.
In a wider study regarding the visual, somatosensory and auditory cortices, Latimer et al. (2019) found that the behaviour of the adaptive responses that they observe can be accounted for by fixed filters that operate over multiple time scales.By developing a method for estimating temporal scales in neuronal dynamics, Spitmaan et al. (2020) found that most neurons exhibited multiple temporal scales in their response, which consistently increased from parietal to prefrontal and cingulate cortex.Miri et al. (2022) in turn suggest that gaze control requires integration over distributed temporal scales.
We propose that if the aim is to build mathematical models of such neural, perceptual or memory processes, then the mathematical theory for time-causal scale-space kernels presented in this paper should be ideally suited for building such models that are both time-causal and time-recursive.Specifically, if the aim is to build such temporal models that can handle multiple temporal scales in a way that respects temporal scale covariance, and under a architectural setting that corresponds to multiple primitive temporal smoothing stages coupled in cascade, then the time-causal limit ker-nel (described in Section 3.1) with its temporal derivatives (described in Section 5) constitutes a canonical class of temporal basis functions to be used in such models.
As a consequence of the temporal delay of such timecausal kernels (Equations ( 10) and ( 39)), any time-causal perceptual process will be associated with an inherent temporal delay (complemented with the processing time of the neural processes that implement the corresponding computations), implying that the representation of the present (White 2020) will, in practice, be a representation of some (temporally extended) temporal moment(s) 18 in the past, unless complemented with extrapolation/prediction (White 2018) over a time period corresponding to the temporal delay(s) of the perceptual process that lead to that percept.Still, however, a representation of the present, with or without temporal prediction implying without or with an inherent temporal delay, will by necessity be a representation of a temporally "fuzzy" present.
In their review of the use of multiple temporal scales in the brain, Cavanagh et al. (2020) state that short temporal windows facilitate adaptive responding in dynamic environments, whereas longer temporal windows promote the gradual integration of information across time, and specifically concerning the notion of multiple temporal scales they conclude a heterogeneity of temporal receptive fields at the level of single neurons within a cortical region, consistent with the aims behind the theory for temporal scale-space representation described in this article.

Implications of the presented theory with regard to the philosophy of time and perceptual agents
The subject of this paper has been to describe a theoretical framework for handling the notions of time and temporal scales for a perceptual system or a neural system, in a both principled and theoretically well-founded manner.Since this subject has implications regarding how we consider the notion of time for a perceptual agent, we will in this section describe relations to the philosophy of time (Mölder et al. 2016;Callender (2017), which is still an open topic in the area of philosophy.
The notion of time is something that we usually take for granted.Still there is no fully established definition for this concept.Already St. Augustine (354-430) stated (Outler, transl. 1955, Book 11, page 193): 18 For a time-causal temporal filtering process that operates over multiple temporal scales, there will, in general, be a different temporal delay for each temporal scale, in the sense that the temporal delay will be shorter for temporal filtering over a short temporal scale and longer for temporal filtering over a long temporal scale.This raises an interesting theoretical problem concerning how to maintain an internally consistent representation of the time-delayed present, given that different components in such a representation may be associated with different temporal delays.
"What, then, is time?If no one asks me, I know what it is.If I wish to explain it to him who asks me, I do not know.Yet I say with confidence that I know that if nothing passed away, there would be no past time; and if nothing were still coming, there would be no future time; and if there were nothing at all, there would be no present time." According to Newtonian or Galilean space-time, we can treat time as flowing continuously and define a universally valid notion of global time.According to Einstein's relativity theory (1905,1916), different observers can measure time differently, being affected by the relative velocity between the observers.Thus, measurement of time is a local property (attached to the path that an observer or a clock follows in space-time), and (at very high relative velocities) different observers may not even be able to agree on the temporal ordering between different temporal events in the world. 192021 This treatment deals with the handling of time for a single perceptual agent that observes a dynamic world using time- 19 To understand how relative temporal ordering can be different for two observers with different relative velocities, consider two observers Alice and Bob in relation to a moving train with constant velocity.Alice is in a waggon of the train, whereas Bob is standing on the ground outside.Let us assume that Alice is positioned in the middle of the waggon, and uses a special flashlight to emit two photons simultaneously from exactly the middle of the waggon, one photon is emitted in the forward direction of the train and the other one is emitted in the backward direction.Since light always travels at the speed of light in relation to any Galilean frame according to Einstein, the two photons in the forward direction and the backward direction will hit the walls in the forward and the backward directions of the train simultaneously, from the viewpoint of Alice.From the viewpoint of Bob, the situation will, however, be different.Since the train moves in the forward direction, the photon emitted in forward direction will have to travel a longer distance from the temporal moment of emission to the temporal moment of arrival than the photon emitted in the backward direction, since the train is moving and changing the positions of the walls during the time it takes for the light to travel between the two positions.Therefore, the photon emitted in the forward direction will hit the wall in the forward direction after the photon emitted in the backward direction hits the wall in the backward direction.Thus, because of the way that space-time is transformed by high velocities, Alice and Bob will arrive at different conclusions regarding the relative temporal ordering of the two events. 20A minor note concerning this thought experiment: If you find the situation artificial in the respect that the two photons arrive at the two walls exactly simultaneously from the viewpoint of Alice, you could modify the thought experiment slightly: Move the flashlight just a tiny bit in the forward direction, so that the photon emitted in the forward direction arrives at the forward wall slightly before the photon emitted in the backward wall from the viewpoint of Alice, but not too much so that that the photon emitted in the forward direction arrives at the forward wall before the photon emitted in the backward direction arrives at the backward wall from the viewpoint of Bob.Then, we have a complete reversal of the temporal ordering of the two temporal events.
21 For a biological perceptual agent, the relativistic time corrections that he or she may encounter due to relative velocities between two observers who observe everyday phenomena in the world with their own perceptual systems only will, however, be much shorter than the inner time scales of their perceptual systems, implying that relativistic time effects can be ignored in a treatment of how to handle the notion of time for a perceptual agent that observes everyday phenomena.causal receptive fields as temporal primitives in its perceptual system.
Originating from a paper by McTaggart (1908), there are two main theories regarding time in the area of philosophy: According to the A-theory, A-series events are ordered by which are present, which are past, and which are future (tensed propositions), whereas according to B-theory, B-series events are ordered by which come before and which come after (tenseless propositions) (Zalta (ed.), Stanford Encyclopedia of Philosophy 2020).Thus, A-theory is closer to how we perceive time as humans (and similar to St. Augustine's view above), whereas B-theory is closer to how we describe temporal phenomena in physical theories of the world.
In a treatment about the notion of temporal presence, Power (2016) discusses how we are able to maintain a perception of changes in the world in our representation of the present.Essentially using the argument that the temporal present is a instantaneous property (valid at a single time moment only), while arguing that the perception of changes requires access to properties of the world over an extended temporal interval, he concludes that A-theory is false, since extended temporal properties cannot exist in a representation of the temporal presence at a single time moment. 22rom the viewpoint of a temporal multi-scale analysis as developed in this paper, where each measurement of properties in the world requires integration over a non-infinitesimal temporal interval, it does, however, follow that any perceptual measurement of the world will have to be performed at some non-infinitesimal inner temporal scale, and thus correspond to integration over a non-infinitesimal duration over time.From such a viewpoint there is no contradiction relative to a perceptual representation of the present, since a multi-scale representation of the present will always occur over multiple temporal scales, and will thus have the possibility to collect information about how properties in the world change over time over extended temporal intervals.
Additionally, in human perception, there are dedicated perceptual mechanisms for registering changes or motion over time;23 compare, for example, with the illusion of the motion after effect (Wohlgemuth 1911), implying that if you look out of a window of a moving train for a long time, and if the train suddenly stops, you may for a while perceive a (physically non-existent) motion in the opposite direction.Alternatively, you may encounter a similar illusion if look-Fig.13 Illustration of non-infinitesimal temporal duration of any physical measurement that arises as a consequence of a non-infinitesimal inner temporal scale in a physical temporal measurement device, specifically for any biological sensory or perceptual system, as well as the non-zero temporal delay of any time-causal temporal receptive field, which implies that the representation at any present moment will de facto instead be a representation of what happened some amount of time ago in the past.For the scale-covariant time-causal limit kernel proposed as the most idealized model of a temporal receptive field in this article, the temporal delay will specifically be proportional to the temporal scale measured in units of [time], thus implying longer temporal delays at coarser temporal scales.(For a physical or biological implementation of these notions, there will also be another complementary temporal delay, not treated further here, caused by the time it takes to carry out the actual computations in the perceptual system.)(The vertical arrow in this illustration is intended to represent the present moment.The blue curve, in turn, reflects how different information from different temporal moments in the past contribute to the representation of the present at that present moment.To represent the temporal duration of the time-causal temporal smoothing kernel, we have in this illustration drawn the "full width half maximum" (FWHM), which is proportional to the temporal standard deviation of the temporal scale-space kernel, in other words proportional to the square root of the temporal scale parameter τ .)ing at the motion of streaming water for a sufficiently long time, and then perceive motion in the opposite direction if you change your viewing direction to focus on a static object.There are also static stimuli that give rise to perception of motion (see e.g.Conway et al. (2005)).
The model for temporal multi-scale processing developed in this paper does thus make the following assumptions concerning the handling of the notion of time for a perceiving agent: The perceptual system of the perceiving agent has a lowest layer of biophysical sensors, which performs temporal integration of the underlying physical signal with some shortest time constant corresponding to the smallest possible inner temporal scale of the perceiving agent.Then, successive layers of such operations are coupled in cascade in a hierarchical manner over that first layer, leading to a layered architecture in the perception system, with successively longer effective time constants at higher layers corresponding to coarser temporal scales.Each such representation in any layer of the hierarchy operates on input information acquired in the present, possibly complemented with access to memory buffers of the past.Thus, from the perspective of the perceiving agent, he or she cannot have any access to the actual physical present in the external world ("das Ding an sich"; Kant 1783Kant , 1902)), but instead just access to a temporally blurred representation of the present, which from the perspective of the perceiving agent is the only available representation of the present24 (see Figure 13).
From the representation of the (temporally blurred fuzzy) present, the internal perceiving system of the agent may also compute representations at coarser temporal scales, which by the temporal delays inherent to the time-causal temporal processes will also serve as a temporal memories of the past.The perceiving agent has no access to a video or audio recording of the past.Instead, the only possible representation of the past is what is stored in the temporal memories of the perceiving agent.25Some of these memories may be of a short term nature and soon be overwritten by more recent information, while other memories may be stored for further longer term access.
A more technical problem in relation to temporal memory concerns making estimates of the duration of a temporal event.According to the standard methodology in physics, one would use a clock, register the times of the beginning and the end of the temporal event and compute the duration from the difference between these temporal moments (a Bseries type of measurement).A biological perceiving agent does, however, not have access to any explicit clock, and there is no evidence for an accurate inner clock in the human brain that a human perceiving agent could relate to for directly measuring the duration of temporal events (Wittman 2009).
From the viewpoint of a temporal multi-scale analysis, it is, however, in principle possible to estimate the duration of a temporal event by operating on representations at multiple temporal scales and comparing the relative strengths of their responses, thus using A-type measurements in the (timedelayed) present as opposed to quantitative B-type temporal relations for estimating temporal duration.In (Lindeberg 2018a), it is shown how it is possible to define multi-scale spatio-temporal visual operations that respond by their strongest Fig. 14 The temporal delays of the time-causal receptive fields resulting from the presented theory call for a mechanism for performing temporal prediction to extrapolate the de facto time-delayed representation of the present (here represented as the temporal peak of the temporal receptive field marked in blue) to a better representation of the actual present (here represented by the vertical line on the time axis), to enable better temporal dynamics for a perceptual agent that interacts with a dynamic world.(Additionally, it is, of course, for other purposes also preferably to also attempt to perform predictions into the actual future in relation to any time moment, to enable temporal planning and to compensate for the time it will take to execute the actions called for by the perceptual agent.The latter types of temporal predictions are, however, not assumed to influence the representation of the present in this treatment.)response over temporal scales at a temporal scale corresponding to the temporal duration of the temporal event, thus estimating the duration of a temporal event based on measurements at a single temporal moment only, although a very special temporal moment at which the response assumes extrema over both time and temporal scales.This is an extension of spatial scale selection (Lindeberg 1998a(Lindeberg , 2021a)), which makes it possible to estimate spatial scales without need for explicitly laying out a ruler. 26ue to the temporal delays of the time-causal receptive fields that drive this perceptual engine over time, any representation of the present will not be a representation of the actual present moment, but instead of what had occurred at some temporal moments (or rather temporal intervals) in the past.Furthermore, representations at coarser temporal scales will harbour the traces of events that occurred further in the past compared to representations at finer scales, thus providing basic mechanisms for temporal memory buffers.
To make it possible for the perceiving agent to handle fast occurring temporal events in a dynamic world, it is therefore extremely valuable for a perceiving agent to be able to perform predictions from the time-delayed perceptual present to at least the actual physical present, so as to be able to coordinate his or her actions with fast occurring temporal phenomena (Figure 14).Given that it will additionally take time to plan and execute an action in practice, it is in a similar way essential that the perceptual agent can perform predictions into the actual future in relation to the actual present moment when planning and executing an event.Even further predictions to the future may of course also be valuable for longer term planning, and to be able to make such longer term predictions, it is very valuable to have an explicit memory of the past over longer temporal scales.Thus, the notion of multiple temporal scales is also important for making predictions into the future, for different time scales into the future.
For the brain of a perceiving agent, its ability to predict what will happen in the future may therefore be one of the most critical factors that determine its ability to survive and reproduce in a competition between individuals and species in the survival of the fittest (Darwin 1859(Darwin , 2004;;Spencer 1864Spencer , 2020)).Minimizing the prediction error, has been proposed as main principle underlying brain function (Friston 2010;McCrone 2022).It has also been argued that the sensory cortex is optimized for prediction of future input (Singer et al. 2018), and furthermore been demonstrated that it is possible to learn the receptive fields of deep neural networks by training the networks to predict the relative future from pre-recorded video sequences of natural scenes (Singer et al. 2018;Kwan and Park 2019;Lotter et al. 2020).Lowlevel neuronal learning mechanisms have also been proposed in terms of predicting future activity (Luczak et al. 2022).
To conclude, we argue that in a A-theory type treatment of time for a perceptual agent, it is essential to complement previous such treatments with explicit notions of (i) noninfinitesimal temporal scales for any representations of the present, and also to incorporate (ii) the unavoidable temporal delays of time-causal receptive fields that determine the functional properties of perceptual systems.In a corresponding manner, given the extended temporal delays of even the fastest temporal processes in e.g.human vision, it is essential to complement the perceptual process with (iii) mechanisms for temporal predictions, since otherwise the actions of the perceiving agent will be too slow to be able to handle and cope with rapid temporal phenomena in the environ-ment. 2728These three notions are immediate consequences of treating temporal perception as a consequence of a temporal measurement problem, where information in physical stimuli has to be integrated over non-infinitesimal durations over time (a main assumption underlying the formulation of the presented temporal scale-space theory), and making a notion such as the instantaneous present de facto impossible for a perceptual agent.
Given the working hypothesis that perception has to involve some mechanisms for temporal prediction to compensate for the non-avoidable temporal delays associated with time-causal temporal integration over non-infinitesimal neighbourhoods over time, our conscious experience of the present in the world, thus has to synthesize a view of the world, created by our brain, and truly corresponding to "controlled hallucination" (Koenderink 2011;Clark 2016;Paolucci 2021;Seth 2021).It is a "hallucination" in the sense that the view of the present is not actually a view of how the world is or was at the moment that it was first registered and then passed on to further processing.It is on the other hand "controlled" in the sense that it is grounded on biophysical measurements of properties in the world, and processed by a biological system that has been refined over evolution over a very large number of generations.
Let us finally emphasize that this treatment does not make any claim of being able to judge about the properties of time itself, which can only be made by physical experiments, possibly complemented by theoretical modelling and analysis, as done in the area of theoretical physics.Instead, the treatment in this section concerns how the notion of time is handled by a perceptual agent, specifically how the notion of multiple temporal scales with their associated temporal de-27 Considering, for example, the sport of playing tennis.For a professional tennis player, the speed of the tennis ball after a serve may be up to the order 200 km/h, corresponding to more than 50 m/s.For among the faster temporal processes in human vision, the temporal limit of visual judgement is of the order of 20 ms, while for among some slower processes the limit of visual judgement is of the order of 100 ms (Holcombe 2009).Regarding spatio-temporal receptive fields in the primary visual cortex (V1), explicit modelling of examples of such receptive fields measured by DeAngelis et al. (1995DeAngelis et al. ( , 2004) ) lead to temporal scale values in the range from 50 to 80 ms (Lindeberg 2016).The composed temporal scale level for the entire visual hierarchy has to be longer, whereas the fastest possible visual reaction times are of the order of 200 ms (Jain et al. 2015).During temporal intervals in the range 20 to 100 ms, the tennis ball will be able to move by the order of 1 m to 5 m, or during 200 ms even as far as 10 m (if we neglect the loss of speed because of the air resistance).As a tennis player, you quickly learn that you have to fixate on the ball in order to be able to hit it properly, which compensates for some amount of the relative motion.Nevertheless, it seems unlikely that we would be able to judge the position and the timing of the ball properly, unless our conscious perception of it involves at least some component of temporal extrapolation or prediction.Specifically, the process of fixating on a rapidly moving target also needs an explicit temporal prediction mechanism. 28For further support of the working hypothesis that our perception of the present likely involves essential components of prediction or extrapolation in the forward direction of time, see Nijhawan (1994Nijhawan ( , 2008)), Grush (2007Grush ( , 2008)), Changizi et al. (2008) and White (2018).lays have to be considered in such a context, with a set of immediate implications thereof.
Let us also stress that the model used as basis for this treatment is continuous in time, whereas for a biological neural system that communicates with spikes between its neurons, the underlying communication channels are in reality discrete, however, here assumed to be operating at a temporal scale below the inner temporal scale of the functional processes in the perceptual system.9 Historical developments of temporal scale-space theory For the reader interested in a historical overview of previous developments of temporal scale-space theory, this section gives an overview of some the main previous contributions in this area that this paper is based on, follows and extends.Koenderink (1988) pioneered the area of the temporal scale-space representation by proposing his scale-time model based on applying Gaussian smoothing over a logarithmically transformed temporal domain.
A complete classification of the general class of continuous scale-space kernels was given in (Lindeberg 1993b).While this classification also included the truncated exponential kernels used as main temporal primitives in this paper, the main topic of that book was spatial computer vision, and the specific detailed structure of time-causal scalespace kernels was at first developed further in the more dedicated treatment in (Lindeberg and Fagerström 1996) aimed at video processing, specifically including the logarithmic distribution of the temporal scale parameter in the set of temporal scale channels.
The topic of temporal scale selection was first addressed in (Lindeberg 1997a), including detailed investigations of the response properties of time-causal receptive fields over temporal scales and time, and illustrating how a closely related temporal model based on the time-causal Poisson kernel, in turn assuming a semi-group property over temporal scales, can also serve as a temporal memory of the past.
In (Lindeberg 1997b(Lindeberg , 2001) ) the time-causal model based on the temporal Poisson kernel, specifically the temporal derivatives of this kernel, was used for modelling the temporal variability in biological spatio-temporal receptive fields.In ter Haar Romeny et al. (2001) the temporal variability in biological spatio-temporal receptive fields was modelled using temporal derivatives of Koenderink's scale-time kernels.
Other temporal scale-space models based on a semi-group property over temporal scales were then studied in (Fagerström 2005(Fagerström , 2007) ) and (Lindeberg 2011).
In (Lindeberg 2016) a substantial theoretical extension was made of the temporal model based on truncated exponential kernels coupled in cascade, by deriving the timecausal limit kernel, which allows for temporal scale covariance.In (Lindeberg 2017) this model was extended to tem-poral scale selection, including detailed studies of the temporal response properties and scale selection properties for the cases of a uniform sampling vs. a logarithmic sampling of the temporal scale parameter.A general proof was also presented, explaining how previous temporal models based on the assumption of a semi-group property over temporal scales lead to poor temporal dynamics, specifically undesirably long temporal delays.
In (Lindeberg 2016) the developments of the time-causal limit kernel were performed in the context of video processing, and were used for deriving theoretical models of spatiotemporal receptive fields with close relations to biological receptive fields in the lateral geniculate nucleus (LGN) and the primary visual cortex (V1).In (Lindeberg 2018a) this theoretical framework for spatio-temporal receptive fields was extended to scale-covariant spatio-temporal feature detection with integrated spatio-temporal scale selection.In (Lindeberg 2018b) corresponding extensions were made for dense temporal scale selection as well as dense spatio-temporal scale selection.In (Jansson and Lindeberg 2018) a specific application to video analysis was developed to analyze dynamic textures in a temporally scale-covariant manner.In (Lindeberg 2021b) the same theoretical model for spatiotemporal receptive fields based on using the time-causal limit kernel and its temporal derivatives as temporal basis functions was used for modelling biological vision in an axiomatic normative theory of visual receptive fields In (Lindeberg andFriberg 2015b, 2015a) parallel developments were made for auditory signals, showing how main classes of time-frequency transforms (spectrograms) can be derived in an axiomatic manner, as well as how auditory receptive fields at a higher level can also be axiomatically derived with very close similarities to biological auditory receptive fields.
Most of the previous developments of the temporal scalespace theory relevant for the treatment in this paper have, however, been performed with regard to visual processing, and in the context of models for spatio-temporal receptive fields.Some parallel developments have on the other hand been performed with regard to auditory processing.
Anticipating that this could be a cause to problems for a reader from a background in biology or signal processing, who is interested in analysing or modelling purely temporal phenomena using a corresponding theory, and wanting to get reasonably quickly into the associated concepts, a first main purpose of this article has therefore been to give a dedicated and self-contained treatment that develops the relevant temporal scale-space theory for the specific domain of purely temporal signals, without having the theory intertwined with concepts regarding spatial or frequency domains, as is the case in the previously available literature, dealing with visual or auditory processing.
We do additionally outline extensions of this temporal scale-space theory to forming time-causal and time-recursive wavelet representations as well as time-causal and scale-covariant time-frequency representations, which do both provide novel contributions with regard to these areas.
With regard to modelling of temporal phenomena in biology, we develop detailed comparisons to other purely temporal models that can be used for such purposes, including ways of translating results from those models to models based on the time-causal limit kernel studied in this paper.With regard to such purposes, we do also extensively relate to previous work on modelling temporal scales in neural signals, for which we proposed that the presented temporal scale-space model could provide a both theoretically and practically valuable tool.Specifically, we present a general procedure for fitting the time-causal limit kernel to nonnegative data, without any need for making use of an explicit expression of the time-causal limit kernel over the temporal domain.
We do finally present implications of the presented theory to fundamental concept formation in the area of the philosophy of time and regarding non-infinitesimal inner temporal scales for any temporal sensor measurement in a perceptual agent, including the resulting inevitable non-zero temporal delays implied by that, in turn implying a need for making predictions into the real present moment, to be able to handle rapid temporal phenomena in the environment.

Summary and conclusions
We have presented a theory for how temporal smoothing of temporal signals can be performed in such a way that it guarantees that the smoothing process does not create new artificial structure in the signal, in the sense that the number of local extrema in the signal, or equivalently the number of zero-crossings, is guaranteed to not increase from finer to coarser temporal scales.Additional critical components of this theory are temporal causality, implying that we are not allowed to access information from the future in relation to any time moment, and temporal recursivity, implying that the temporal smoothing process should not require any other temporal memory of the past than the resulting temporal scale-space representations themselves.
A complete classification of the linear and shift-invariant convolution kernels that obey these properties has been given, based on an earlier treatment in (Lindeberg and Fagerström 1996), in turn based on earlier classical results by Schoenberg (1948Schoenberg ( , 1950)).For continuous signals, the corresponding temporal scale-space kernels consist of truncated exponential kernels coupled in cascade, corresponding to firstorder integrators coupled in cascade, and for discrete signals, first-order recursive filters coupled in cascade (Section 2.2).
As a conceptual extension of this general approach, we have described a specific subset of choosing these kernels in such a way that temporal scale covariance is obtained.The corresponding time-causal limit kernel that permits scale co-variance, which is a novel construction in (Lindeberg 2016), is the limit case of an infinite number of truncated exponential kernels coupled in cascade, with specific choices of the temporal time constants (Section 3.1).
Temporal scale covariance in this context means that if the input signal is rescaled by some uniform temporal scaling factor S = c i , where c is the distribution parameter of the time-causal limit kernel and i is some integer, then the result of performing temporal smoothing on the rescaled temporal signal is the same as performing temporal smoothing on the input signal, followed by a corresponding rescaling of the processed original signal, and complemented by a shift of i units along the scale dimension (Section 3.1.3).
These temporal kernels, optionally combined with their temporal derivatives, do in this way constitute a canonical class of temporal basis functions for numerous purposes of temporal modelling, in situations when the temporal operations have to be time-causal and time-recursive, and in addition have the ability to handle temporal information over multiple temporal scales in a theoretically well-founded manner.With appropriate scale normalization of the temporal derivatives, the temporal derivatives of the time-causal limit kernel are also truly scale covariant, with preserved magnitude values of temporal derivatives at matching temporal scale levels under scaling transformations, in turn allowing for truly scale-invariant processing under temporal scaling transformations of the input signal (Section 5.3).
We have given an explicit expression for the time-causal limit kernel in the Fourier domain (25) and although the kernel lacks a compact closed-form expression over the temporal domain, we have shown how it can be related to other temporal models, such as Koenderink's scale-time kernels (Section 3.3) and the ex-Gaussian model, which is the convolution with an exponential kernel with a single truncated exponential function (Section 3.4).We have also presented a general methodology for how the parameters in a model based on a (temporally either unshifted or time-shifted) timecausal limit kernel can be determined from lower-order temporal moments of some other temporal function or temporal signal (Section 3.4.2and Appendix A.3).
We have described how these kernels can be implemented on discrete data, based on a set of first-order recursive filters coupled in cascade, where also the discrete implementation guarantees that new local extrema, or equivalently new zero-crossings, cannot be created from finer to coarser levels of scale (Section 4).The discrete implementation of temporal derivatives is straightforward, in terms of small support finite difference operators applied to the discrete temporal scale-space representation (Section 5.5).Thus, the discrete implementation is highly efficient and lends itself to realtime applications.
We propose that the presented theory, serving as a normative theory of purely temporal receptive fields, provides a canonical way of defining multi-scale representations of temporal signals in situations where the signal operations have to be truly time-causal, because of lack to access of future information in real-time scenarios, and time-recursive, because of a need to keep memory buffers of the past to a minimum in terms of memory requirements.Specifically, we propose that the time-causal limit kernel with its temporal derivates constitutes a canonical class of temporal basis functions in situations when the temporal scales may vary, especially when temporal scale covariance and temporal scale invariance are desirable properties.
We have also related the theory to other approaches for processing temporal signals at multiple temporal scales, specifically wavelet analysis and time-frequency analysis.We have outlined how the temporal derivatives of the time-causal limit kernel can serve as time-causal and time-recursive wavelet bases (Section 6.1) and how a complex-valued extension of the time-causal limit kernel can be seen as time-causal analogue of Gabor functions, in turn enabling truly scalecovariant time-frequency analysis also over time-causal and time-recursive temporal domains (Section 6.2).
Concerning applications of the presented theory, we have described how these time-causal kernels constitute a canonical class of temporal kernels for modelling spatio-temporal and spectro-temporal receptive fields in biological perception (Sections 7.1-7.2).We have also given a more general overview of the applicability of multiple temporal scale levels in perceptual, memory and cognitive processes in biological nervous systems, as well as given arguments proposing that the time-causal kernels treated in this paper should constitute a corresponding canonical class of temporal kernels when modelling neural signals as well as more general perceptual and temporal memory processes by explicit mathematical models (Section 7.3).
Finally, we have presented general arguments for the need for incorporating the notion of non-infinitesimal temporal scales with their associated non-zero temporal delays when considering a perceptual representation of the present (not the same concept as the instantaneous actual present, which a perceptual agent has no possible access to), which then also leads to a direct need for temporal extrapolation or prediction in order to compensate for the temporal delays associated with the time-causal temporal filtering operations in a time-causal perceptual system (Section 8).We propose that these arguments should have essential implications for the logical reasoning in A-type theories of time in the philosophy of time, as well as when modelling perceptual agents.
In (Bright et al. 2020, see Equations ( 2) and (3)), the authors fit a socalled ex-Gaussian model, which is the convolution of an unnormalized Gaussian function with an unnormalized truncated exponential kernel, to the temporal response functions of neurons.With slightly different naming of the variables to avoid notational clashes with the notation used elsewhere in this article, let us consider a temporal response function of the form A.1 Second-order moment-based method without flexible temporal offset parameter In this appendix, we will derive a relation between the above ex-Gaussian model and a corresponding model based on the time-causal limit kernel h limit-kern,gen (t) = b 0 + b 1 Ψ (t; τ, c), (85) with the time-causal limit kernel Ψ (t; τ, c) in (26) defined from its Fourier transform according to (25).
For simplicity, let us first assume that we are in range of the parameter space of the ex-Gaussian model where the temporal delay is small relative to the standard deviation and the time constant µ, such that we do not need to introduce an additional temporal delay in the model (85) based on the time-causal limit kernel.Let us also assume that we can assume that the DC levels in the two models should be equal, such that we can throughout assume that b 0 = a 0 .Then, our task is to derive a mapping to compute the parameters b 1 , τ and c in the model based on the time-causal limit kernel from the parameters a 1 , m, σ and µ in the ex-Gaussian model.
The approach that we shall follow is to compute the zero-, firstand second-order temporal moments of the two models with the DCoffsets a 0 and b 0 suppressed and determine the mapping between the parameters of the two models from the requirement that the integral, the temporal mean and the temporal variance should be equal.
Computing the (uncentered) temporal moments up to order two of the ex-Gaussian model (86) in Mathematica gives Using the fact that the temporal mean and the temporal variance of the time-causal limit kernel are given by (Lindeberg 2016, Equations (34) and ( 35)) identifying these expressions and solving for b 1 , c and τ in the model (87) based on the time-causal limit kernel gives which with δ and V according to ( 91) and (92) as well as M 0 , M 1 and M 2 according to (88), ( 89) and ( 90) gives the desired mapping between the ex-Gaussian model ( 83) and the model ( 85) based on the time-causal limit kernel.

A.1.2 Experimental results
Figure 10 shows examples of ex-Gaussian temporal models approximated by time-causal limit kernels in this way.A conceptual advantage of the time-causal limit kernel in this context, is that we do not need to use or modify a Gaussian kernel to model the initial transient phenomena in a time-causal temporal response function that decays towards zero in an exponential manner towards the tail.In this way, a neural response modelled by the model based time-causal limit kernel would also correspond to a biologically plausible implementation corresponding to temporal integration of the form illustrated in Figure 2.
A.2 Extension to a third-order moment-based method involving an additional temporal offset parameter As a remark concerning extensions, if the ex-Gaussian model is in a range of the parameter space where the temporal delay is large relative relative to temporal duration of temporal onset of the composed kernel, then an additional temporal offset t 0 can be added to the model ( 85) based on the time-causal limit kernel and an additional computation and identification of the third-order central moments be performed to determine also this parameter in the mapping between the two types of temporal models.
The explicit expression for the unnormalized and uncentered thirdorder temporal moment of the ex-Gaussian model for a 0 = 0 is whereas the expression for the normalized and centered third-order moment of the time-causal limit kernel is (Lindeberg 2016, Equation (36)) Using the relationship between the centered and uncentered third-order moments we obtain the following expression for the normalized and centered third-order moment of the ex-Gaussian model (102)

A.2.1 Method for third-order moment-based model fitting
To determine the parameters in the model based on the time-shifted time-causal limit kernel (98) with the DC-offset disregarded (b 0 = 0), we can hence proceed as follows: 1. Compute the unnormalized and uncentered moments M 0 , M 1 , M 2 and M 3 of the ex-Gaussian model according to (88), ( 89), ( 90) and ( 101). 2. Compute the variance V of the ex-Gaussian model according to (92) and let the variance τ of the time-causal limit kernel be equal to this value according to (94). 3. Identify the normalized and centered third-order moments of the ex-Gaussian model and the model based on the time-causal limit kernel according to (102) and (100).4. With the third-order moment κ 3 of the ex-Gaussian model computed according to (102) and the variance τ of the time-causal limit kernel according to (94), square the expression (100) and solve the resulting fourth-order algebraic equation in terms of the distribution parameter c of the time-causal limit kernel.This will give four roots for c, out of which only two of the roots can be expected to satisfy the original unsquared equation, because of the squaring operation that may introduce new false roots. 5. Select 29 the real root of the original equation (100) that additionally satisfies c > 1.Then, determine the temporal offset t 0 of the time-shifted time-causal limit kernel in (98) from the normalized and centered first-order moment of the time-causal limit kernel with δ identified with the normalized first-order moment of the ex-Gaussian model according to (91).6. Compute the amplitude b 1 of the time-shifted time-causal limit kernel in (98) according to (95).This procedure can either be carried out purely numerically or in a package for symbolic computation, such as Mathematica. 29The skewness measure γ 1 = κ 3 /V 3/2 = 2(c+1) √ c 2 − 1/(c 2 + c + 1), which is used for determining the distribution parameter c in the model based on the time-causal limit kernel with a flexible temporal offset t 0 , increases monotonically with c for c > 1 and assumes values in the range ]0, 2[.Hence, provided that the skewness measure determined from κ 3 and V is in this range, there will a unique real root for c that satisfies c > 1.

A.2.2 Experimental results
Figure 15 shows the result of applying this procedure for fitting a timeshifted time-causal limit kernel to an ex-Gaussian model that does not obey the assumptions for fitting a model based on the non-shifted time-causal limit kernel to it according to the previous second-order moment-based method.In the left figure, the result of the second-order moment-based method is shown, demonstrating a substantial difference because of the fixed zero offset of the original time-causal limit kernel.The right figure shows corresponding results for the third-order moment-based method, demonstrating a much better agreement between the two models, when an additional degree of flexibility is introduced into the model based on the time-causal limit kernel by adding the temporal offset parameter.
A.3 Fitting models with the time-causal limit kernel to other functions or signals Note that with replacement of the moments M 0 , M 1 , M 2 and optionally M 3 with the moments of some other non-negative function or signal, the same overall procedures can more generally be used for fitting models based on the time-causal limit kernel to other one-dimensional signals or functions that: (i) are defined for positive values of time, (ii) assume non-negative values only, (iii) have a roughly unimodal shape of first increasing and then decreasing and (iv) tend to zero towards infinity.The second-order moment-based fitting approach is in this context intended for situations when the temporal origin of the signal or function is known in advance and in some sense intended to be minimal, whereas the third-order moment-based fitting approach is intended for situations when the temporal origin in the data is unknown and hence needs to be adapted to each situation.
B Appendix B: Implementing temporal filtering with a discrete approximation of the time-causal limit kernel This appendix gives a brief explicit description about how to implement temporal filtering of a sampled discrete signal with a discrete approximation of the time-causal limit kernel.
For simplicity, assume 30 that the input signal has been sampled with a unit time increment ∆t = 1.Then, given a temporal standard deviation of the kernel σ in such units of time, compute its variance τ = σ 2 and choose a suitable value of the distribution parameter c > 1 that determines the sampling density in the temporal scale direction.1. Compute a set of temporal scale levels τ k according to a geometric distribution (12): 2. Compute a corresponding set of scale increments: with the additional definition τ 0 = 0. 3. Compute the time constants µ k for a set of temporal recursive filters with generating functions of the form (54) according to (55): 4. Couple the following sets of first-order recursive filters in cascade (53): Note that, in a real-time scenario or an offline scenario where memory efficiency is important, if the task is to compute a single temporal scale level only, such as the first temporal scale level in a cascade, this operation can be performed without explicitly storing the representations at the intermediate temporal scale levels, except for at the current and the previous temporal frames.Furthermore, when computing multiple temporal scale levels in parallel, the temporal scale-space representation at the next coarser temporal scale is most efficiently computed by applying a single recursive filter to the temporal scale-space representation at the nearest finer temporal scale (if we assume a dense representation over temporal scale levels, where all the temporal scale levels are assumed to be used in the later processing stages). 5. Optionally, compute discrete approximations of scale-normalized temporal derivatives for some γ > 0 (where γ = 1 is a standard default value) by applying the following discrete derivative approximation operators (according to Equations ( 60) and ( 65)) to the temporally smoothed signal, alternatively instead using Lpnormalization according to (61) as opposed to variance-based normalization according to (60).

Fig. 2
Fig.2Electric wiring diagram consisting of a set of resistors and capacitors that emulate a series of first-order integrators coupled in cascade, if we regard the time-varying voltage f in as representing the time varying input signal and the resulting output voltage and f out as representing the time varying output signal at a coarser temporal scale.Such first-order temporal integration can be used as a straightforward computational model for temporal processing in biological neurons; see alsoKoch (1999, Chapters 11-12)  regarding physical modelling of the information transfer in the dendrites of neurons.

Fig. 4
Fig. 4 Approximations of the time-causal limit kernel for τ = 1 using K = 7 truncated exponential kernels in cascade and their first-and second-order derivatives.(top row) Logarithmic distribution of the scale levels for c = √ 2. (bottom row) Logarithmic distribution for c = 2. (Horizontal axes: time.Vertical axes: function values.)

Fig. 5
Fig.5Illustration of temporal smoothing of two signals: (left) a Brownian noise signal generated from a simulated Wiener process and (right) a synthetic sine wave signal f (t) = sin(exp((b − t)/a)) for a = 200 b = 1000 with temporally varying frequency so that the wavelength increases with time t, computed using a discrete approximation of the time-causal limit kernel for c = 2 in terms of a set of recursive filters coupled in cascade.Observe how fine-scale structures corresponding to higher frequencies are successively suppressed when going from finer to coarser temporal scales, and also that the temporal scale-space representations at coarser temporal scales are associated with longer temporal delays, in this figure seen as different offsets in the positions of the peaks in the temporal signal at different temporal scales.(Horizontal axes: time.Vertical axes: signal values.)

Fig. 7
Fig.7Comparison between the two different ways of computing the representation in the upper right corner in Figure6from the corresponding representation in the lower left corner, using either the clockwise direction (marked in green) or the counterclockwise direction (marked in red).When generating this illustration, we have first essentially performed a rescaling of the scale-space representation of the signal in the left column and marked the result as solid green curve, and then overlayed the scale-space representation of the signal in the right column with a dashed red curve.(Technically, in the discrete implementation, we have, however, instead visualized the equivalent result of such a computation at a lower resolution, to avoid the formally ill-defined operation of interpolating the discrete signal in the left column to a higher resolution, and instead subsampled the signal in the right column, which explains the change in the labelling of the temporal axis.)(left) The result for the entire temporal interval used in the right column in Figure6.(right) Enlargement of a central region of the temporal interval.As can be seen from the visualization, the results computed in the clockwise or counterclockwise directions are basically indistinguishable, demonstrating the scale covariance property of the temporal scale-space representation defined by convolution with the time-causal limit kernel.(The result is best viewed by zooming in to a digital copy of the article.)(Horizontal axes: time.Vertical axes: signal values.) temporal scale representations then constituting memories of what has happened further in the past.
Fig.9Comparison between (blue curves) the time-causal limit kernel according to (25) and approximated using the first K = 7 components of the infinite convolution of truncated exponential kernels in cascade with its first-and second-order temporal derivatives and (brown curves) the temporal kernels in Koenderink's scale-time model (37) and their first-and second-order temporal derivatives.All kernels correspond to temporal scale (variance) τ = 1 with the additional parameters determined such that the temporal mean values (the first-order temporal moments) become equal in the limit case when the number of temporal scale levels K tends to infinity (equation (38)).(top row) Logarithmic distribution of the temporal scale levels for c = √ 2 (bottom row) Corresponding results for c = 2. (Horizontal axes: time.Vertical axes: function values.)

Fig. 11
Fig.11Illustration of the result of computing discrete approximations of second-order scale-normalized temporal derivatives L ζζ from the time-causal temporal scale-space representation L at different scales (using distribution parameter c = 2 and scale normalization power γ = 1), here for a synthetic input signal f consisting of two temporal peaks generated as discrete approximations to time-causal limit kernels for temporal scales τ = 16 and τ = 256 with a certain amount of relative temporal delay to separate the responses as well as a small amount of added white Gaussian noise.(Horizontal axes: time.Vertical axes: Signal values.)

Fig. 12
Fig. 12 Graphs of the complex-valued extension χ(t, ω; τ, c) = Ψ (t; τ, c) e iωt of the time-causal limit kernel Ψ (t; τ, c) for temporal scale τ = 1 with distribution parameter c = 2 and different values of the angular frequency ω. (left column) The real component, corresponding to the time-causal limit kernel multiplied by a cosine wave.(right column) The imaginary component, corresponding to the time-causal limit kernel multiplied by a sine wave.(top row) Angular frequency ω = 5. (bottom row) Angular frequency ω = 20.(Horizontal axes: time.Vertical axes: kernel values.)