E Teftef

M.-J Escobar

A Astudillo

C Carva Jal

B Cessac

A G Palacios

T Viéville

F Alexandre

Modeling non-standard retinal in/out function using computer vision variational methods

Keywords:

We propose a computational approach using a variational specication of the visual front-end, where ganglion cells with properties of retinal Konio cells (K-cells), are considered as a network, yielding a mesoscopic view of the retinal process. The variational framework is implemented as a simple mechanism of diusion in a two-layered non-linear ltering mechanism with feedback, as observed in synaptic layers of the retina, while its biological plausibility, and capture functionalities as (i) stimulus adapted response; (ii) non-local noise reduction (i.e. segmentation); (iii) visual event detection, taking several visual cues into account: contrast and local texture, color or edge channels, and motion base in natural images. Those functionalities could be implemented in the biological tissues We use computer vision methods to propose an eective link between the observed functions and their possible implementation in the retinal network base on a two-layers network with nonseparable local spatio-temporal convolution as input, and recurrent connections performing nonlinear diusion before prototype based visual event detection. The numerical robustness of the proposed model has been experimentally checked on real natural images. Finally, we discuss in base of experimental biological and computational results the generality of our description.

Introduction

Recently it has been proposed that the retina, a accessible part of the brain, sustain more complex behaviors than expected before [START_REF] Masland | The unsolved mystery of vision[END_REF], including spatial, temporal and motion recognition [START_REF] Schwartz | Sophisticated Temporal Pattern Recognition in Retinal Ganglion Cells[END_REF]) [START_REF] Olveczky | Segregation of object and background motion in the retina[END_REF]. The retina correspond to a multi-stream device including from phototransduction to early visual processing and neural coding, at dierent spatial and temporal scales. At the application level, a computing parallel information streams system for (e.g.) visual prosthesis is likely to need an important post-processing level before visual signals feed the nervous system (Barriga-Rivera & Suaning, 2011). For the later, an adequate computational architecture should have at least two-layers network with non-separable local spatio-temporal convolution as input, and recurrent connections performing non-linear diusion before prototype based visual event detection. The present article is organized as follow: In the section 2, we review key biological facts at the early dorsal and ventral K-cells visual streams, computing sophisticated spatial and temporal pattern recognition as review in the next section, 3, base on [START_REF] Gollisch | Eye Smarter than Scientists Believed: Neural Computations in Circuits of the Retina[END_REF][START_REF] Litke | What Does the Eye Tell the Brain?: Development of a System for the Large-Scale Recording of Retinal Output Activity[END_REF][START_REF] Schwartz | Sophisticated Temporal Pattern Recognition in Retinal Ganglion Cells[END_REF][START_REF] Olveczky | Segregation of object and background motion in the retina[END_REF][START_REF] Sterling | How retinal circuits optimize the transfer of visual information[END_REF]). In the subsequent section 4, we implement the computational principles raised by the study on the retinal Kcells [START_REF] Hendry | The Koniocellular Pathway in Primate Vision[END_REF][START_REF] Yoonessi | Functional assessment of magno, parvo and koniocellular pathways; current state and future clinical applications[END_REF] using a variational specication of the visual front-end, base on a network of retinal ganglion cells. In the section 5, we implemented a numerical model to study real image sequences and nally model predictions are presented to veried biological validity.

2

From standard to non-standard early-vision front-end

Standard front-end retinal streams

The retinal output correspond to several dierent spatio-temporal processing, base on a diversity of ganglion cells (Gc) types, the output of the retina to the nervous system, of the incoming image sequence. Such streams are embodied as separate strata that span the retina at the Gc surface.

It is generaly accepted, that the (i) parvo (P-cells) and (ii) magno (M-cells) streams, compute respectively (i) color image details and contrast in central vision and (ii) monochrome intensity temporal variation, including at very low contrast, in peripheral vision and (iii) konio streams (K-cells) projecting to the LGN (middle layer) receive information from blue cones, sending information to the V1 color blobs [START_REF] Hendry | The Koniocellular Pathway in Primate Vision[END_REF] and from small bistratied Gc having a surround which is sensitive to yellow (G. Field et al., 2007). Here, we are going to focus on a subset of K-cells, beyond their color property, as detailed in the sequel. For other types of Gc see, e.g., [START_REF] Callaway | Local circuits in primary visual cortex of the macaque monkey[END_REF].

At the modeling level, information streams (i) and (ii) are well represented by LN models, i.e., a one layer spatio-temporal ltering followed by a static non-linearity, as represented in the right part of Fig. 2, this unit being tuned by a gain control mechanism, as reviewed e.g., in [START_REF] Wohrer | Model and large-scale simulator of a biological retina with contrast gain control[END_REF]. Qualitatively, such non-linear ltering tends to remove spatial correlations in the visual world, producing a less redundant output, as required to transmit information in the optic nerve bers of limited spatial capacity (see [START_REF] Pitkow | Decorrelation and ecient coding by retinal ganglion cells[END_REF] for a discussion), resulting in a sparse coding for retinal spike trains in response to the statistics of natural images [START_REF] Simoncelli | Vision and the statistics of the visual environment[END_REF].

Figure 1: Schematic representation of projections of the K-cells from [START_REF] Masland | The unsolved mystery of vision[END_REF] and [START_REF] Nassi | Parallel processing strategies of the primate visual system[END_REF].

Considering retinal output K-streams While standard magno/parvo Gc downstream correspond to the standard visual dorsal/ventral streams (leftward drawing) after a relay in the ventral/dorsal LNG layers (rightward drawing); K-cells projections are more miscellaneous (leftward drawing) and interact with the parvo/magno streams in the LGN and via the V1 2/3 (rightward drawing). Connections from the retina are feed-forward, while all other brain connections include feed-backs. In a purposive view of visual perception, the parietal cortex performs from such input an unconscious, ecient, specialized and rapid processing of the whole visual eld, and prepares actions adapted to the characteristics of the environment. In addition K-cells connect via the LGN to important multi-sensorial areas like the amygdala lateral nucleus (ALN) involved in the control of emotion [START_REF] Ciocchi | Encoding of conditioned fear in central amygdala inhibitory circuits[END_REF], thus in direct link with survival actions.

Importance of the K-streams Quantitatively, 80% of the Gc are midget towards the P-cells in the thalamus, because of the need of ne spatial resolution to perceive visual details. Only 10% are parasol Gc towards the M-cells in the thalamus, and 10% are bistratied Gc (K-cells) towards the superior colliculus and the LGN layers (G. D. Field et al., 2007). Though such non-standard retinal cells are in proportion few in primates (whereas up to 75% in phylogenetically earlier mammalians) they constitute robust and well-identied cell mapping of the visual eld [START_REF] Gauthier | Receptive elds in primate retina are coordinated to sample visual space more uniformly[END_REF]. This lower proportion is easily explained: Having large receptive eld (RF) about 10 deg, the cover of the whole retinal eld required less units.

In the human eye, about 130 millions of photoreceptors, where 120 millions are rods and 10 Inria millions are cones, concentrate onto about 1.5 millions of Gc, including the 10 5 K-cells considered here.

The key computational aspect of K-streams

The current computational assumption here consider that K-cells streams provide rough but fast event or object -detection of visual events [START_REF] Hendry | The Koniocellular Pathway in Primate Vision[END_REF]) [START_REF] Masland | The unsolved mystery of vision[END_REF] in order to induce survival goal-directed actions on time (V. [START_REF] Lamme | The distinct modes of vision oered by feedforward and recurrent processing[END_REF], and to drive higher-level iterative processes with prior information [START_REF] Callaway | Local circuits in primary visual cortex of the macaque monkey[END_REF][START_REF] Koivisto | Recurrent Processing in V1/V2 Contributes to Categorization of Natural Scenes[END_REF]. Given natural image sequences, we interpret these two functionalities as image segmentation for visual objects detection and natural statistical recognition, including temporal pattern recognition. In the other hand, K-cells participate to blind-sight (V. A. [START_REF] Lamme | Blindsight: the role of feedforward and feedback corticocortical connections[END_REF][START_REF] Cowey | Visual System: How Does Blindsight Arise?[END_REF], i.e., or not consciously seen, including, slowly moving targets detection, rough localization of stimulus. This is typically information not for seeing but for visually guided behavior. There are some evidences, that such computation starts in the retina:

(i) Some retinal Gc are able to produce sophisticated spatial and temporal pattern recognition on their own [START_REF] Gollisch | Eye Smarter than Scientists Believed: Neural Computations in Circuits of the Retina[END_REF][START_REF] Werblin | The retinal hypercircuit: a repeating synaptic interactive motif underlying visual function[END_REF], including motion detection, directional selectivity, local edge detection, object motion and looming detection.

(ii) During phylogenetical evolution, the eye was not a simple feed-forward system: Feed-backs from the remainder of the nervous system had been able to produce adaptive learning of sophisticated visual functions optimizing the transfer of information [START_REF] Sterling | How retinal circuits optimize the transfer of visual information[END_REF].

(iii) Fast stimulus categorization develops roughly 150 ms after stimulus onset [START_REF] Thorpe | Speed of processing in the human visual system[END_REF]. In that respect the visual processing need to perform this using only few computation steps [START_REF] Viéville | Using an hebbian learning rule for multi-class svm classiers[END_REF].

(iv) From an information processing point of view, regarding downstream computation, it is optimal to take into-account the whole spatial and temporal resolution of the photo-receptors before the necessary compression occurring in the optical nerve to implement such complex non-linear ltering [START_REF] Simoncelli | Vision and the statistics of the visual environment[END_REF].

These are the key elements that dierentiate the K-cells from parvo/magno streams. In nutschell, the K-stream provides fast, a-priori event detection assumptions, to the remainder of the visual system (V. [START_REF] Lamme | The distinct modes of vision oered by feedforward and recurrent processing[END_REF].

Retinal architecture and biological constraints

What are the biological constraints regarding what can be computed in the retina (i.e., the computational characteristics)?.

In order to derive a mesoscopic model of K-cells input/output function let us collect known facts about the computational properties of the retina, as schematized in Fig. 2.

Here are the three major features we propose to take into account:

1. Not one but two layers: It is clear that, by no means, a simple feed-forward ltering layer, even non-linear, can compute complex visual cues such as those reviewed above. On the other hand, it is known that a two layers non-linear feed-forward network, i.e., with a hidden layer, is a universal approximator of any function. For a static input/output relationship, i.e., a unique image as input with a unique related output, a neural network

RR n°8217

Figure 2: left: Schematic representation of the retinal architecture, after [START_REF] Wohrer | Model and large-scale simulator of a biological retina with contrast gain control[END_REF] and [START_REF] Gollisch | Eye Smarter than Scientists Believed: Neural Computations in Circuits of the Retina[END_REF]. Excitatory connections are in red, inhibitory in blue. Gap junctions are in red with a black contour, neighboring cells of the same type being generally linked through gap junctions. The direct pathway from light receptors to Gc is a two-layers process dening two successive non-linear ltering stages with dierent excitatory-inhibitory patterns.

The present mesoscopic model is compatible with such an architecture, but requires a subset of this architecture (highlighted in yellow). This pathway is modulated by the interstitial cells in each layer. Here, we have emphasized that non-linearities are present on each link. right:

The basic view of cell unit processing: (i) spatial ltering (gray arrow) specied by the wired connectivity, combined by a (ii) band-pass temporal lter, thus including delays, (blue rectangle) and (iii) followed by a thresholding operation.

(e.g. retina) with two layers is a sucient architecture for generic computations [START_REF] Hornik | Multilayer feedforward networks are universal approximators[END_REF]. For a dynamic input, i.e., an image sequence, this basic idea has to be generalized, and we propose the following.

2. Non separable spatio-temporal ltering: At a given time, we consider that the retina input state is a function of the short-term visual information and that the output is a causal function of this 2D+T image volume.

Quantitatively, we may consider the last 150 ms interval as temporal depth [START_REF] Li | Dierent retinal ganglion cells have dierent functional goals[END_REF], with a visual minimal event-time of about 10 ms: This thus correspond to a dozen of frames when time is discretized.

-This means that the spatio-temporal ltering is local, i.e., only spatial connections from neighborhood cells are taken into account, while in the temporal domain the ltering is only due to delays in the biological substrate and connections. This seems to be fairly true for the outerplexiform layer (OPL) (i.e. for horizontal cells) which corresponds to the ltering mechanism considered here. This is less obvious for amacrine cells in the IPL [START_REF] Werblin | The retinal hypercircuit: a repeating synaptic interactive motif underlying visual function[END_REF], likely involved in mechanisms not addressed here.

-This means that the spatio-temporal lters are constrained: not all weight values can be used, but only values with signs compatible with the excitatory/inhibitory connections, and value ranges taking into account the nature of the connection (e.g., gap junction versus synaptic junctions). All synapses have a delay, increasing as a function of the connection length.

-This means that the output at a given time is a static function a 2D+T volume, thus computing local motion in the image sequence, but no long-term motion cue. Therefore the previous argument of a two-layers architecture is sucient for complex function approximation still stands for dynamic stimuli in this restrained framework.

-This also means that the spatio-temporal ltering is non-separable: In other words, we do not consider a sequence of 2D image (thus spatial ltering each image and then temporal ltering the result) but the ltering of the 2D+T volume.

In fact, as made explicit before [START_REF] Gollisch | Eye Smarter than Scientists Believed: Neural Computations in Circuits of the Retina[END_REF][START_REF] Dong | Spatiotemporal Inseparability of Natural Images and Visual Sensitivities Inria Motion Vision[END_REF], the retinal structure induces a separable spatio-temporal ltering if we consider a cell alone in a standard parvo-stream, but dierent if we consider interactions between cells, and also consider that the retina is adapted to natural image stimuli, as suggested by [START_REF] Dong | Spatiotemporal Inseparability of Natural Images and Visual Sensitivities Inria Motion Vision[END_REF].

3. Recurrent connections: The general architecture denitely shows that the hard-wiring provides feed-forward connections and recurrent (horizontal and feed-back) connections, as made explicit in Fig. 2.

-Recurrent connections as local diusion. These connections are local, continuous and with tiny delays (i.e. based on gap-junctions on local non-spiking very short dendrite/axon connections) [START_REF] Sterling | How retinal circuits optimize the transfer of visual information[END_REF]. Qualitatively, this corresponds to local diusion of the state value, here mainly corresponding to the membrane voltage value.

-Static non-linearities everywhere. Another key aspect is the fact that each connection is subject to a static non-linearity with a main characteristic: Threshold corresponding to a rectication of the signal. The rectication is not related to action potential spiking (except for the Gc output), but to bio-chemical membrane mechanisms.

RR n°8217

Functional characteristics of the retinal network

The retina is able to compute local isotropic contrast and intensity temporal variation detection, but requires specic functionalities.

-A-stimulus adapted response

From raw light intensity event detection (including dim light ashes), to RF responses corresponding to the basis element of natural images decomposition [START_REF] Hyvärinen | Natural Image Statistics[END_REF], the primary aspect of the retinal input-output function is to be adapted not only to the general statistics of natural images, as observed in the LGN projections [START_REF] Dan | Ecient coding of natural scenes in the lateral geniculate nucleus: experimental test of a computational theory[END_REF], but also to dynamic changes in the statistics of a given environment [START_REF] Hosoya | Dynamic predictive coding by the retina[END_REF].

Such stimulus adaptation is implemented by a suitable choice of the local 2D+T convolution kernel induced by the local network connectivity [START_REF] Wohrer | Model and large-scale simulator of a biological retina with contrast gain control[END_REF]. Such lter elements span the space of possible relevant local input [START_REF] Hyvärinen | Natural Image Statistics[END_REF].

However, beyond the hard-wiring of such connections, a local adaptation must occur to optimize the response. This also means that the retinal ltering is not only adapted to any natural image sequence, but likely also adapted to a the present visual environment, as observed at the LGN level [START_REF] Truccolo | Dynamic temporal decorrelation: An information-theoretic and Inria biophysical model of the functional role of the lateral geniculate nucleus[END_REF].

A way to interpret this, is to generalize the contrast adaptation mechanism proposed in, e.g. [START_REF] Wohrer | Model and large-scale simulator of a biological retina with contrast gain control[END_REF], to second-order local statistics adaptation. This means to assume that the retinal processing is able to tune not only the gain but also the local spatio-temporal correlations to a given environment and to react (i.e., to detect) a specic second-order spatio-temporal statistics in the signal. This is the rst main computational assumption of our work, with two consequences:

1. On one hand, as illustrated in Fig. 3, second-order statistics is in one to one correspondence with the magnitude of the (here: 2D+T) spectrum 1 . Detecting second-order spatio- temporal statistics means detecting a given spectrum magnitude signature. On the reverse, the spectrum phase is a function of the higher-order and second-order statistics, so that assuming second-order processing means that the retinal adaptation and detection capability is only related to the spectrum magnitude signature. This is due to the Wiener-Khintchine theorem that relates the signal auto-correlation, thus the second-order moments, the power spectral density. It has been shown that, in fact, this is a relevant cue to characterize natural image categories [START_REF] Torralba | Statistics of natural image categories[END_REF], e.g., discriminate between human-made environment or not, categorize a given landscape, etc. We would like to hypothesize that what is true for large scale natural image categories is also true for small scale image categories, e.g., detect stone vs vegetation, etc. This is plausible since the natural image statistics power spectra is 1/f , meaning that the second order statistics is scale invariant, which suggests (but it is not established) that the same cue may be relevant at a local scale.

This is going to be veried numerically in section 5. This is not obvious and this does not means that the whole statistics is scale invariant. It is true, in average, for natural images [START_REF] Simoncelli | Natural Image Statistics and Neural Representation[END_REF]. It is not the case for a given subset of images and there are evidence at least in V1 that even if the statistic is preserved our visual system treats images dierently [START_REF] Simoncini | Nature Neuroscience[END_REF].

2. On the other hand, tuning the second-order but not the higher-order statistics has a very The semantic of the spatial spectrum magnitude versus the spectrum phase: When (on the left, from [START_REF] Hyvärinen | Natural Image Statistics[END_REF]) corrupting the magnitude while preserving the phase, the image (here of a bull and an insect) is still recognizable, whereas (on the right, for the well-known Lena image) when canceling the phase and preserving the magnitude the image has still natural image statistics but the shape has vanished. Here, we propose as a minimal model, to consider that K-cells processing is based on second-order information, thus on spectrum amplitude, since in one-to-one correspondence 1 .

RR n°8217

a signal with respect to its correlations, thus second-order statistics, than an independent component analysis (ICA) which decomposes the signal with respect to higher-order moments.

The fact that early vision optimized for either only a second-order statistics or for higherorder statistical dependencies is an open question and recent results tends to show that V1 RF function not only make their outputs as independent as possible, but also facilitate the processing of higher-order statistical regularities [START_REF] Karklin | Is Early Vision Optimized for Extracting Higher-order Dependencies[END_REF]. Similarly, it is assumed that high-order statistics are taken into account in redundancy reduction [START_REF] Barlow | Redundancy reduction revisited[END_REF].

However, results suggest that adaptation in the retina system is not sensitive to changes in commonly measured higher-order statistics, retinal cells exhibit remarkably invariant behavior to changes in higher-order statistics [START_REF] Tka£ik | Retinal adaptation and invariance to changes in higher-order stimulus statistics[END_REF].

This invariance is coherent with what has been observed in V1 with pink-noise, i.e., the fact it is sucient to imitate second-order statistics of natural images to obtain responses similar to natural image statistics [START_REF] Kayser | Responses to natural scenes in cat V1[END_REF]. To this debate, we bring the following argument: If natural images categories are essentially characterized by their spectrum magnitude [START_REF] Torralba | Statistics of natural image categories[END_REF], thus second-order momenta, it is sucient to detect visual events related to such second-order cues.

No mistake, this statement only concerns the K-cell processing considered here, whereas analyzing natural images using ICA at the retina level is another important issue considered elsewhere [START_REF] Keskes | Modélisation des premières étapes de la vision avec des réseaux de neurones[END_REF]. And we are going to numerically verify that it is sucient to use such cues to perform the complex early-vision processing that is attributed to K-cells.

We also wonder, if, since the eye is expected to produce visual input from natural images anyway, it is relevant not to transmit the obvious information that the image is natural, but the information remainder. How could it be, in that case, that though invariant with respect to higher-order momenta (such as kurtosis) the retina still computes higher-order momenta? In our framework, this is due to another processing step : Non-local ltering, now described and which is a non-linear transformation of the image, so that second-order and higher-order momenta are mixed, thus higher-order momenta indirectly taken into account, even if the basic detection mechanism is only due to second-order cues.

-B-non-local noise reduction

Since the main biological constraint of the retina in terms of information transmission is the 10 2 compression factor for P-cells up to 10 3 for other Gc cells in the optical nerve, reducing the noise with respect to the signal is a major issue [START_REF] Simoncelli | Natural Image Statistics and Neural Representation[END_REF]. The noise added by the photo-transduction mechanism is obviously to be ltered [START_REF] Wohrer | Model and large-scale simulator of a biological retina with contrast gain control[END_REF]. However, what is to be considered as pertinent signal vs distracting information or noise in the incoming visual ow is task dependent. For instance, the background motion induced by ego-motion is a noise with respect to object motion recognition but a relevant signal with respect to gaze stabilization (the opto-kinetic reex for instance has the background motion as main input cue).

Furthermore, there is a local notion of noise (i.e., a random spurious additive signal, locally correlated in space and time, or not), but regarding the visual perception, there are other notions:

A notion of outliers (i.e., spurious piece of signal not to be taken into account to avoid bias, e.g., a spurious light reection) and a notion of inliers (i.e., piece of signal coherent with the ground base data but inducing a mis-perception, e.g., a shadow.

The concept of noise reduction is no more related to a local lter, but to non-local operators, since the context matters. In other words, ltering means estimating a given zone in the image, Inria thus to perform image segmentation. This non-local idea of ltering corresponds to the grouping of homogeneous information, reducing small or spurious non-signicant variations, and then produce a reduced representation in which outliers are eliminated because they're not homogeneous with respect to the information neighborhood, while inliers are grouped in another region to be further treated as signal or noise depending on the perceptual task. Such mechanism of synchronization is well-known in the retina [START_REF] Gollisch | Rapid neural coding in the retina with relative spike latencies[END_REF][START_REF] Shlens | Synchronized ring in the retina[END_REF] where the shared noise from photoreceptors is a major cause of correlated ring (Greschner et al., n.d.), so that synchronized Gc outputs correspond to activity with the same stochastic signal characteristics.

Furthermore, it is important to point out that we do not group regions of homogeneous intensity but regions of homogeneous response to the 2D+T ltering discussed above, e.g., of the intensity, intensity variation, contrast, local edge or motion pattern, etc. It is important to remark that time, here, is considered only as a local feature. As a consequence, such mechanism should be able to segment on complex local spatio-temporal patterns. In particular, as explained in [START_REF] Gollisch | Eye Smarter than Scientists Believed: Neural Computations in Circuits of the Retina[END_REF] after [START_REF] Hochstein | Linear and nonlinear spatial subunits in Y cat retinal ganglion cells[END_REF], texture motion, even if the average local intensity remains constant, is locally detected: During texture motion, dierent intensityvariation cells respond at dierent time to dierent spatial patterns at their RF locations, showing the same activity prole, inducing a cumulative eect related to strong rectication (summing only positive responses in the neighborhood with both excitatory and inhibitory responses). It is a perspective of this work to numerically verify that texture motion is also compatible with such variational approach [START_REF] Olveczky | Segregation of object and background motion in the retina[END_REF]) [START_REF] Aubert | Can the nonlocal characterization of sobolev spaces by bourgain etal. be useful to solve variational problems?[END_REF].

To further understand this aspect of retinal processing, we are going to formalize in details to which extents the retinal architecture is able to not only smooth but also segment the visual signal in order to provide a multi-scale representation of the input in terms of regions of homogeneous signal values (upper-scale), with a coding at ner-scale of the region indexing for a given unit. In this framework we thus only consider the spiking rate and phase, and are going to code them by two analog signals, the former being a positive value, the latter being a relative (i.e., only dened by the dierence with another value). Since phases synchronize, we obtain only a discrete set of distinguishable values and this is equivalent to an indexing.

-C-visual event detection

The retina is able to detect sophisticated visual events such as:

(i) Fine spatial information, as observed qualitatively with natural image sequences, a Gc being often either silent or sparsely by deterministic ring for a given image category [START_REF] Soo | Fine Spatial Information Represented in a Population of Retinal Ganglion Cells[END_REF];

(ii) local approaching motion, which is known as being computable by a local divergence lter [START_REF] Subbarao | Bounds on time-to-collision and rotational component from rstorder derivatives of image ow[END_REF], i.e., the weighted sum of dedicated local masks.

Such functionality includes anticipation, either motion extrapolation or omitted stimulus response [START_REF] Schwartz | Sophisticated Temporal Pattern Recognition in Retinal Ganglion Cells[END_REF], these complex computations coming about through the interplay of rather generic features of retinal circuitry [START_REF] Gollisch | Eye Smarter than Scientists Believed: Neural Computations in Circuits of the Retina[END_REF]. In both cases, as expected, spatio-temporal pattern yields adaptation (delay reduction in the former case, response habituation in the later) and an unexpected event induces a strong reaction. The time window order of magnitude is the second. This is thus not directly derived from a 2D+T ltering mechanism (its temporal window being too short) but by a local fading memory of previous input, thus using regressive feed-backs, and the simpler model we propose here is a simple auto-regressive mechanism. The key point is that this induces a switching mechanism, the detection of a visual event modifying the local parameters on the retina. With (Gollisch

RR n°8217

& Meister, 2010) we assume that the switching may be initiated by much more subtle image features, precisely by the detection of a given visual event.

Regarding the detection of such fast visual event, it has already been investigated in detail how Gc and subsequent layers can produce early responses to a visual event [START_REF] Rullen | Rate coding versus temporal order coding: What the retina ganglion cells tell the visual cortex[END_REF][START_REF] Thorpe | Seeking categories in the brain[END_REF]. We assume here, that if such a mechanism is able to produce a detection at early stages in the visual system, it is likely produced by a large contribution of the retinal processing. This is plausible because the corresponding bio-plausible most ecient computational method is a one layer method implementable a simple non-linear lter with threshold eect [START_REF] Viéville | Using an hebbian learning rule for multi-class svm classiers[END_REF]. This is related to the so-called Support

Vector Machine (SVM) algorithm, and the bio-plausible implementation is a simple competition between units coding similarities with respect to prototypes. This may be implemented as a two steps process: correlation with prototyping responses and competition between the outputs, as implemented in a standard one-layer architecture with lateral inhibitory connections, which is the case of the OPL.

The statistical learning theory and the related SVM algorithm [START_REF] Viéville | Using an hebbian learning rule for multi-class svm classiers[END_REF] allows us to assume that only a little set of prototypes (the support vectors) has to be recruited for the detection of an event. Moreover, approximate prototypes are sucient, since the prototype of a category can be wrong as far as it is less wrong than a prototype of another category. It is thus reasonable to assume that this is a biologically plausible functional model of the retina, since it is able to detect rather complex spatio-temporal visual events. As already pointed-out, the major argument is that it is a pertinent to use the full retinal resolution to perform it before the optical nerve spatial compression.

This view is also coherent with the observation that in the presence of natural stimulus Gc have a sparse and deterministic behavior. Sparse because they react in terms of detection of not a given visual event. Deterministic in the sense that they re a spike response at reproducible times. This is also coherent with observations already reviewed here, that the response is not limited to a local visual eld, since there is no signicant visual events at such small case, but at higher sizes of the visual eld. When considering the objects size in a standard natural scene (?, ?), the fact that the K-cells response corresponds to events in a 10deg of the visual eld seems properly tuned to natural visual events [START_REF] Hendry | The Koniocellular Pathway in Primate Vision[END_REF].

Variational specication of the visual front-end

Let us now translate into equations the previous elements, as developed elsewhere [START_REF] Viéville | Using an hebbian learning rule for multi-class svm classiers[END_REF][START_REF] Viéville | How do high-level specications of the brain relate to variational approaches?[END_REF] and adapted to the retinal processing.

Non-local IPL ltering

The main mechanism is illustrated in Fig. 4, where non-local ltering is precisely dened from the well-known work of [START_REF] Mumford | Boundary detection by minimizing functionals[END_REF], and following here the presentation of [START_REF] Aubert | Can the nonlocal characterization of sobolev spaces by bourgain etal. be useful to solve variational problems?[END_REF]. A step further, following [START_REF] Gérard | From variational to spiking network imagesegmentation techniques[END_REF], we consider the well-dened discrete approximation of the Mumford-Shah criterion proposed by [START_REF] Chambolle | Image segmentation by variational methods: Mumford and Shah functional and the discrete approximation[END_REF][START_REF] Chambolle | Discrete approximations of the MumfordShah functional in dimension two[END_REF]. This allows us to implement this mechanism, taking the retinal IPL architecture into account, as follows. Over a 2D array of pixels, at a resolution h, Inria it writes:

min O ij | W * I(ω ij ) -O(ω ij )| 2 Input ltering + ij,dω∈{(0,1),(1,0)} min α |O(ω ij ) -O(ω ij + dω)| 2 , β h
Region & border optimization , with I(ω ij ) = ω∈P ixelij I(ω)/h 2 being the integrated value of a pixel. It has been shown by [START_REF] Chambolle | Image segmentation by variational methods: Mumford and Shah functional and the discrete approximation[END_REF] that this discrete criterion is a well-dened discretization of such a function.

Contrary to other numerical schemes [START_REF] Aubert | Can the nonlocal characterization of sobolev spaces by bourgain etal. be useful to solve variational problems?[END_REF], its distributed implementation corresponds to a simple-layer network computation with feed-backs. Another technical aspect is the fact that here we consider not a scalar but a vectorial image, and in, e.g., [START_REF] Viéville | How do high-level specications of the brain relate to variational approaches?[END_REF], it has been derived how such mechanisms trivially generalize to the multi-channel case. More precisely, the minimization is implemented by the following non-linear iterative lter:

O(ω ij ) t+1 ← O(ω ij ) t + δ ∇O(ω ij ) with: ∇O(ω ij ) ≡ W * I(ω ij ) -O(ω ij ) + dω∈{(0,1),(1,0)} 0 if |O(ω) -O(ω + dω)| 2 > β h α α (O(ω + dω) -O(ω)) otherwise .
Making explicit this recurrent equations allows us to relate it to the IPL retinal architecture discussed previously. It is thus an iterative lter, implemented as a local feed-back with recurrent connections and is based on local diusion, as made explicit in Fig. 2 and discussed in section ??, regarding recurrent connections. The input gain is positive (excitatory connections) and the local diusion (inuence of the neighborhood activity) is positive too. On the reverse, the local feedback is negative, again as predicted by our knowledge of the IPL architecture. The proposed iterative lter is based on a piece-wise linear feedback, with a very precise static non-linearity:

if the local contrast magnitude |∇O(ω)| 2 is higher than a threshold, diusion is canceled. This is easily implemented by local inhibition, triggered by a contrast threshold, again compatible with elements of Fig. 2. This is the reason why we propose that the IPL ltering stage may be considered as implementing such non-local ltering.

Visual event detection in the OPL

As far as the choice of the non separable 2D+T convolution kernels are concerned, beyond usual contrast and image intensity variation kernels, we thus assume in this modeling that natural image local categories can be related to a spectral magnitude signature, i.e., normalized spectral magnitude prole. We thus hypothesize that retinal computational units are tuned to a given spatio-temporal spectral response Wn (f ). This means that a given input I corresponds to the pattern Wn , if its spectrum I(f ) Between regions, we want the border length to be as small as possible to draw smooth and minimal edges between regions. Minimizing the weighted sum of these three criterion allows us to specify an image segmentation, tuning the ne/coarse grained segmentation with β and tuning the region scale with α. The success of the method comes from the fact that, under realistic conditions, the optimization converges towards a piece-wise a constant intensity prole, i.e., perfect homogeneous regions. Implementing this minimization, from the derivation of the criterion gradient using the related Euler-Lagrange equation, leads to a distributed algorithm with precise local weights adaptation rules [START_REF] Viéville | How do high-level specications of the brain relate to variational approaches?[END_REF]. Region indexing (pixels of the same region have similar output values) and edge detection (edge pixels have non negligible local contrast) is also obviously provided by such process.

a negligible measure for the given dot-product. Applying usual equalities we can relate this mechanism to a simple convolution with normalization and static non-linearity as follows:

Wn (f ) ≃ λ I(f ) ⇔ ǫ ≃ f | Wn(f ) | W| -I(f ) | Ī| | 2 = ω | Wn(ω)) | W| -I(ω) | Ī| | 2 ⇔    Wn(ω) * I(ω)| ω=0 | Ī| ≃ f Wn (f ) Wn (f ) = 1 W⊥ n (ω) * I(ω)| ω=0 | Ī| = f W⊥ n (f ) W(f ) = 0
This obvious textbook derivations make explicit two facts: if spectral magnitude signatures are similar, this simply means that the two normalized signals are similar, while using the reference signal as a normalized convolution (convolution after normalization of the input signal) yields a unitary response.

Local spectral signature detectors are thus easily implementable in the biological substrate as a convolution with a normalized input signal, as it is the case for a simple neural cell with gain control, a static non-linearity providing the thresholding mechanism after the convolution, to cancel negative or negligible parts of the signal. A step further, if two coupled convolution kernels are coupled with their negation, i.e. with kernel representing the orthogonal response, the visual event detection can be based on more robust mechanisms of margin comparison, detecting the event if the positive channel response is signicantly higher than the negative one.

An eective implementation of visual events detection

In order to address the question How is it possible for the retina `alone' to eciently compute, on natural image sequences, such sophisticated visual responses, we have developed a software [START_REF] Teftef | Formalisation de la transformation analogique / événementielle des mécanismes non-standards des cellules ganglionnaires de la rétine[END_REF][START_REF] Teftef | Formalization of the input/output retinal transformation regarding non-standard ganglion cells behavior[END_REF]) that implements the previous described functionalities.

This piece of code is based on the CImg 2 image processing open-source library, and is itself a piece of open-source code available under the K-mrs 3 nick-name. For the generic segmentation routine, we reuse the Léonard Gérard implementation [START_REF] Gérard | From variational to spiking network imagesegmentation techniques[END_REF] of the related Chambolle algorithm [START_REF] Chambolle | Discrete approximations of the MumfordShah functional in dimension two[END_REF], and have contributed to CImg on this aspect. The event detection SVM mechanism has been implemented using the open-source libsvm 4 of Chih- Chung Chang and Chih-Jen Lin.

As far as software engineering is concerned, the key point is that we work on vectorial data volumes, i.e. an image is a data structure with pixels indexed by their 2D location (x, y) and temporal depth t, each pixel being a vector c of channels, starting with the red, green, blue channels, as input, extended with computed channels (e.g., spectrum signature). As a consequence, all functions developed here are usable for both static and dynamic event detection, and all related processes.

While the numerical verications proposed here only consider static images, spatio-temporal convolutions and other operators are available for the software to process 2D+T images sequence slices.

The SVM learning of the spatial event detection prototype was a complete standard process:

We have selected images of variable sizes of each category, and run the calculation of the intensity channels and spectral channel, as shown and explained in Fig. 7. Given the red (R), green (G)

and blue (B), original channels we have chosen the redundant coding : G-R, G-B, G in addition to the intensity R+G+B. The only reason was to obtain numbers rather stable, given the data The key point, for the present process, is the choice of the prototypes. As pointed out in [START_REF] Viéville | Using an hebbian learning rule for multi-class svm classiers[END_REF], the very powerful statistical learning theory, leading to the SVM algorithm, can be revisited as a competitive comparison between prototypes. These prototypes, are the support vectors, close to the hyper-surface that separates two categories. The key parameter is the margin between these prototypes and the frontier between two categories: The higher the margin, the more robust the detection on the subsequent data set (generalization property). In the present numerical experimentation, a standard SVM supervised mechanism have been used.

Inria

set intensity distribution, with for instance, a predominance of green, while red or blue coloration of the image is often relative to green zones. We recall that the ltering is not linear, so that even if the SVM algorithm is invariant by an ane transformation of the features space, the segmentation step depends on such coding choice.

Obviously, the representation is easily adaptable, for instance to dichromate visual systems, like in rodents, using B+G or UV+G, so somehow dierent from trichromats vision [START_REF] Peichl | Eye and vision in the subterranean rodent cururo (Spalacopus cyanus, Octodontidae)[END_REF]. As soon as these two colors allow to discriminate the natural scene objects, the previous framework applies. This is in coherence withj the fact that the retinal visual system is adapted to the natural environment of the animal [START_REF] Sterling | How retinal circuits optimize the transfer of visual information[END_REF].

This event detection is performed on the non-local ltered images (in fact it has been qualitatively experimented that event detection does not work on the raw images in this context).

The choice of considering the sum of the spectrum amplitude responses at 50% of the maximal spectrum response is a pure experimental choice, as detailed in Fig. 7, but guided by two considerations. On one hand, this can be implemented by a simple linear / non-linear lter, as follows: each point of the spatial spectrum corresponds to the correlation of the signal by the corresponding trigonometric functions. The maximal spectrum response always corresponds to the continuous (or DC) component of the spectrum, i.e., the average intensity value. Tuning, in a given 2D direction, the spatial frequency in order the response to correspond to 50% of the DC component, is a simple gain control feedback. Finally, computing the signature reduces to a simple summation. By no means we assume here that this is what happens in the retina, nor that this is the only solution. In fact, as developed in the previous section, the biological substrate oers a much more generic mechanism, computing the correlation between the observed and expected spectrum amplitude.

Though, as discussed previously, all functions can be implemented as distributed computations on a biologically plausible neural network, for obvious computation time performances issues, we have used optimized code in this numerical implementation. Previous studies on these aspects [START_REF] Viéville | Using an hebbian learning rule for multi-class svm classiers[END_REF][START_REF] Viéville | How do high-level specications of the brain relate to variational approaches?[END_REF] have extensively addressed the issue of both the segmentation and detection implementation in a biological network.

The originality of the numerical experimentation is the fact that we have considered a rather specic set of visual stimuli, as made explicit in Fig. 6, and involving studying a specic animal model (see [START_REF] Delgado | The GABAergic system in the retina of neonate and adult Octodon degus, studied by immunohistochemistry and electroretinography[END_REF][START_REF] Ardiles | Postsynaptic dysfunction is associated with spatial and object recognition memory loss in a natural model of Alzheimer's disease[END_REF] for more details on these aspects) 5 . Restraining the issue to computer vision process only, it means that we have considered degraded visual conditions (with respect to state of the art image sequence data base), the issue being to check whether and to which extents such processes still work in such restrained conditions.

In order to experiment in a rather tricky context, we have applied the proposed framework to the ability to recognize shapes in the image [START_REF] Masland | The unsolved mystery of vision[END_REF], i.e. restraining to 2D+T volumes with a thickness of one. We thus have implemented the fact that we consider that the retina is able to compute sophisticated second-order visual cues, including spectrum amplitude signature, since this is implementable as a linear/non-linear lter, see Fig. Interesting enough is the fact that zones 3. and 5., despite their very dierent visual aspect have a rather similar spectrum signature, dierent from 1. Though zones 2. and 4. have spurious spectral responses (for 2. this is the artifact in the sky that yields the spectrum and for 4. the spectrum is not relevant in such a dark zone), since the system is also using the local color as a cue, this yields no mistake. We numerically observed that not only the whole image second-order statistics [START_REF] Torralba | Statistics of natural image categories[END_REF], but also very local second-order statistics seem to be relevant for natural image categorization. More precisely, in the present numerical implementation limited to static images, the color and the sum of the spectrum amplitude responses at 50% of the maximal amplitude are taken as cues to characterize the local texture, in order to perform region detection. The fact that such very simple scalar parameter calculated on the spectrum magnitude leads to relevant results (see Fig. 9) is a good argument in favor of the proposed choice. The generalization to spatio-temporal volumes, thus 2D+T spectra, is straightforward and is a direct perspective of this work. Interesting enough, several key points can be observed here : The output is morally an articial image corresponding to the natural input, where the forms in the image have been preserved.

More precisely, the edges (e.g., the owers boundaries) are preserved even when the image is smoothed : The fact the mechanism lters the noise and not the edges, is simply due to the non-linearity of the diusion operator, since small variations diuse and are thus ltered, whereas higher intensity variations related to edges do not. Furthermore, diusion is anisotropic, more precisely performed along the edge, but not across it (in a more mathematical language : Tangential to the local edge orientation, but not orthogonal to it), thus preserving it. In the present computer implementation, region synchronization is implemented by the propagation of a phase index as expected in a biological neural network. In the sky, the artifact (likely a tree limb) has disappeared, because one non-local eect of this process is small enough so such regions are absorbed though the non-linear diusion. In a nutshell, we have here a full resolution, but simplied image. to directly detect zones with potential pray/predators hidden in the green. The categorization algorithm has been trained with a very small learning set (a few dozen of samples) and the recognition test is no more than a comparison of the region synchronized parameters (color, spectral signature) with respect to prototypes in competition. Such functionality is easy to generalize to another visual events (e.g. snake crawling, predator approaching motion, etc..) and thus provides an eective early-vision generation of a-priory information. We have numerically veried that such detection is robust along the image sequence. However, we also experimented that the non-linear ltering (i.e., segmentation) step is mandatory: Without this sophisticated early-vision step, the present categorization has very poor performances. This output represents sparse responses of the retina in the presence of specic visual events.

Inria

RR n°8217

this article.

6 Discussion: What is predicted by such a modeling ?

With respect to usual models regarding the retina processing, the present work simply proposes to instantiate usual ideas about the so-called non-standard Gc behaviors, but with a computational framework that make them compatible with realistic image processing. Furthermore, this framework raises some new predictions about retinal processing that are not yet tested, but could be worth tested to better understand the biological behavior of the retina.

Here are experimental some facts that will falsify the present framework, regarding the K-cells behavior considered here.

About adaptation

We consider K-cells stream and speculate that the related processing mainly involve second-order statistics, because this is sucient to characterize natural image categories. It may thus be a second-order issue [START_REF] Simoncelli | Natural Image Statistics and Neural Representation[END_REF]) [START_REF] Torralba | Statistics of natural image categories[END_REF]. A simple way to falsify the present statement would be to study to which kind of stimulus statistics a K-cell in the LGN is sensitive to. We expect such cell to be invariant to global pink noise whitening (e.g., to response to an articial image with the same features a natural one except its spatial spectra magnitude).

About segmentation

We claim that the IPL ltering not only reduces local noise but performs an implicit clustering of visual zones, likely based on non-linear diusion mechanisms yielding synchronization between cell responses of the same class. As a consequence, the compression process is not only based on very local spectral or statistical properties, but tuned to non-local cues. This means that the context must strongly inuence the local response. More precisely, a Gc response must be averaged and synchronized with respect to Gc with similar responses (i.e., cells of the same region) but not with cells of another one. If we can observe, for instance, un-synchronized cells corresponding clearly to the same value of the corresponding visual cue in the same image area, the present model is lapsed.

A key falsiable consequence of the present assumption is the fact that the Gc channels output image (i.e., 2D map) must not correspond anymore to the statistics of natural images with a scale invariance of a spectrum in 1/f . This is simply incompatible with the non-local region segmentation mechanism. If such statistics are preserved, the proposal is defeated. In other words, the decoding of the LGN transmitted image to the brain must have the look of an articial image more than a natural one.

A less falsiable consequence on this topic is the fact that high spatial frequencies are ltered, but not in a homogeneous way: inside a region high-frequency are expected to be canceled, information diusion leading to a uniform responses of cells averaging the response value, whereas on the border between regions, no ltering occur. This is a weaker argument because, anyway, as the spatial frequencies are higher on edges than between edges. Here the inhomogeneity in terms of spatial high-frequency ltering must have been amplied by the IPL.

Inria

About detection

We nally have considered the visual event detection to be based on a prototype competition.

This means that a given visual event must be represented by a few spatio-temporal patterns and if the input is close to one of these patterns, the Gc detection must re. If, on the contrary, the representation of a given visual event (e.g. a crawling motion) is fully distributed, as when represented by additive information or dynamic attraction, the proposal is to be deprecated. Another key aspect is that such prototype must represent event and no-event cases, since the underlying algorithm compares two alternatives to trigger an event. In other words, the mechanism is a comparison with respect to a ground state, not a threshold. Furthermore, we expect border-line prototypes rather than typical prototypes, to be coded for the detection. This means that we expect the local spatio-temporal activity to be correlated with responses that, for dierent categories, are close together. This is again due to the fact that the system does not have to detect an event in absolute way but relatively to an alternative. Therefore, it is more important to be able to characterize the data with respect to prototypes close to the border of such alternatives (notion of support prototypes), than to conrm that the similarity with a caricaturistic prototype. If such characteristic is detected at the experimental level, the proposed framework becomes plausible, otherwise it remains questionable.

At the mesoscopic level, i.e. considering the response of a given retinal stream at the network level, it implies that a single cell type may serve quite dierent roles (e.g., approaching motion and structure motion), because it is recruited for a dierent event detection mechanism, depending on the current state of the network. This has already been partially observed, but if this observation is marginal, whereas the majority of cells are dedicated to a unique function, the model fails.

A less falsiable consequence is that we have to assume that the observed switching between local scale functions must be initiated by much more subtle image features, i.e., the visual event context. If such switching is coherent with the related event detection (i.e., switch does not occur when stimuli correspond to detecting similar events, but does otherwise), the proposed modeling is not a fake.

Conclusion

We propose here the use a variational framework to model and simulate, given natural image sequences, the mesoscopic collective non-standard behavior of some retinal input/output functions that correspond to the output of a subclass of the so-called K-cells. We hypothesize that from sophisticated temporal pattern recognition, to image segmentation, or specic natural statistical recognition, a unique generic two-layers non-linear ltering mechanism with feedback is implemented in the biological tissues, while not the individual but the collective behavior of the retinal cells answer for such input/output functions. Taking the retinal architecture and related biological constraints into account and considering the wider class of early-vision non-standard sensory-motor functionality known as non-standard behaviors, we use computer vision methods to propose an eective link between the observed functions and their possible implementation in the retinal network.

At the application level, we simply propose a generic visual front-end to be used for either the simulation of the sensory input of a larger simulator of the brain behavior (e.g., considering sensory-motor loops, embedded vision, ..) or for industrial applications considering such similar visual input, which concretely means degraded visual input.

Such framework raises falsiable predictions to be veried at the biological level and points out how the retina can computationally deliver visual event candidates driving subsequent visual processes. With such framework, it is not possible to assume that the brain visual system stands on the reception of an homogeneous visual pipe of ltered images coming from the retina,

RR n°8217

whereas it is connected to several heterogeneous sources of information at dierent spatial and visual scales, and dierent integration levels, while tuned to the natural image statistics. At the application level, on one hand, this completely changes the way articial systems could be inspired by the biological retina. Roughly speaking, this encourages to compute in parallel several information streams depending on the sensory-motor task, instead of thinking of a large information pipe. On the other hand, this makes very questionable the idea of designing relevant visual prostheses by simply connecting a visual sensor array to the brain, as it is. In a nutshell, visual prostheses are likely to be though as a plural of sensors with a non-negligible data processing, before feeding the nervous system.

  Figure 3:Top view : Relevant representation of natural image categories from the spectrum magnitude signature (here the 60, 80 (in blue) and 90% (in red) level curves), from[START_REF] Torralba | Statistics of natural image categories[END_REF]: Each environment class has a relatively specic signature. In a bio-plausible framework such signature correlation can be computed by a linear/non-linear network subset, since the distance to a given spectral signature is equivalent to the energy response of a lter with spectral property corresponds the given signature, as made explicit in the sequel. Bottom view:

Figure 4 :

 4 Figure4: The computational principle of the Mumford-Shah image segmentation variational approach for a multi-channel image, after[START_REF] Aubert | Can the nonlocal characterization of sobolev spaces by bourgain etal. be useful to solve variational problems?[END_REF]. Given an area of pixels Ω and an input image I : Ω → R 3 which associate to each pixel its color (here a red-green-blue (RGB) value), the goal is to compute an output image O : Ω → R C , segmented in regions, writing K the border between regions. The specication is three-fold. (i) Input ltering. The output image is closed to the local ltering of the incoming image by 2D+T convolution kernels parameterized by the weights W and computing C channels such as contrast, intensity variation, or local spectral property. (ii) Region homogeneity. Inside a region we want the output value variations to be minimal (here ∇ stand for the gradient operator, thus the value variation), yielding homogeneous regions. (iii) Border parsimony. Between regions, we want the border length to be as small as possible to draw smooth and minimal edges between regions. Minimizing the weighted sum of these three criterion allows us to specify an image segmentation, tuning the ne/coarse grained segmentation with β and tuning the region scale with α. The success of the method comes from the fact that, under realistic conditions, the optimization converges towards a piece-wise a constant intensity prole, i.e., perfect homogeneous regions. Implementing this minimization, from the derivation of the criterion gradient using the related Euler-Lagrange equation, leads to a distributed algorithm with precise local weights adaptation rules[START_REF] Viéville | How do high-level specications of the brain relate to variational approaches?[END_REF]. Region indexing (pixels of the same region have similar output values) and edge detection (edge pixels have non negligible local contrast) is also obviously provided by such process.

Figure 5 :

 5 Figure 5: Schematic representation of middle-scale visual events as a support-vector-machine (SVM) implemented in terms of a competitive winner-takes-all connectivity. Here, areas of 3 to 10 deg of the visual eld are considered with, as examples, owers (upper-left), gravel-oor (middle) and vegetation (upper-right). The discrimination between such regions is based on the trivial competitive response of the dierent detectors. Such winner-takes-all mechanism is obviously implemented with spike-coding (Rullen & Thorpe, 2001) since the higher the response, the sooner the spike generation, which can inhibit other responses and thus win the competition.

Figure 6 :Figure 7 :

 67 Figure 6: The La Campana database: a natural habitat of Octodon degus rodent (40 sequences of 10 2 images of about 800 × 800 pixels). In order to be able to perform both biological and numerical experiments on degus visual system, image sequences of about 10 2 images each have been taken, in the natural environment of the animal, at the height, lighting conditions, image resolution and with displacements roughly corresponding to the rodent usual displacements. Itis important to note that this species use to live in a specic visual environment: In the bottomleft view, the red-green-glue (RGB) intensity histogram is not Gaussian as it is usually the case[START_REF] Simoncelli | Natural Image Statistics and Neural Representation[END_REF], but biased, especially in the green channel (oscillation in the distribution), due to the statistics of the vegetation colors. In the bottom-central view, two spectrum magnitude proles in log coordinates, taken over hundred of frames, show that the 1/f d , d ≃ 1 expected prole (thus a line with a negative slope in the gure) does not occur here, whereas the spectrum is a function of such visual context.

Figure 8 :

 8 Figure 8: Qualitative examples of non-local ltering using diusion mechanisms as discussed in (Viéville et al., 2007) regarding their biological implementation. Lower images are the non-local ltering output of one channel (here intensity) computed on the corresponding upper image.

Figure 9 :

 9 Figure 9: Two examples of detection of visual events (here, spatial events only). Top views correspond to the input image, bottom views to the detection output. The most salience regions are shown in both cases, i.e. those with maximal size and/or channel average values (here, for the present gures, size only is considered). Left view A sky zone and four dark/light tree zone have been detected. Right view The main tree zone have been detected, while the ground is less salient since not structured in terms of regions. Clearly this kind of information is relevant

RR n°8217

Inria

http://cimg.sourceforge.net

M-mrs, for KEOpS mesoscopic retina simulator http://project.inria.fr/keops/K-mr

http://www.csie.ntu.edu.tw/∼cjlin/libsvm RR n°8217

As an example, an interesting issue regarding the will to work with such rodent specic visual data, is to wonder whether the rodent visual system is optimized for such visual environment, or if on the contrary such rodent visual system is phylogenetically conserve with respect to such ecological variation.RR n°8217

351, Cours de la LibérationBâtiment A 29

Publisher Inria Domaine de Voluceau -Rocquencourt BP 105 -78153 Le Chesnay Cedex inria.fr ISSN 0249-6399

Acknowledgment A big thank you to Stéphane Deny for recent discussions and Pierre Kornprobst for older discussions, both at the origin of some important aspects of this work. We also acknowledge FONDECYT Nro.1120570 (Chile), CONICYT PIA Nro. 79100014 (Chile), ANR-47 CONICYT (Chile) and DGIP USM 231117.

project and CORTINA associated team.