Skip to main content

How Uncertainty Bounds the Shape Index of Simple Cells

Abstract

We propose a theoretical motivation to quantify actual physiological features, such as the shape index distributions measured by Jones and Palmer in cats and by Ringach in macaque monkeys. We will adopt the uncertainty principle associated to the task of detection of position and orientation as the main tool to provide quantitative bounds on the family of simple cells concretely implemented in primary visual cortex.

Mathematics Subject Classification (2000)2010:62P10, 43A32, 81R15.

1 Introduction

One of the fundamental tasks performed by simple cells in primary visual cortex is that of detecting position and local orientation of a stimulus [1]. On the other hand, the functional behavior of simple cells as visual detectors is characterized in terms of standard linear filtering and with other so-called nonclassical behaviors [2]. We will concentrate on linear aspects, and consider classical receptive profiles modeled with a planar oscillation under a spatially localizing window. In [3], such receptive profiles were studied in terms of two dimensionless indexes of shape ( n x , n y ) corresponding to the product of the frequency of the oscillation and the sizes of the window in the direction of the oscillation and in the orthogonal one, showing that the distribution of such feature on V1 simple cells of macaque monkeys is confined to a specific region. This result is summarized in Fig. 1. Remarkably, the same confinement was found also in cats [4, 5], and this suggests that this pattern can be associated with some criteria of optimality with respect to perceptive tasks. Notable proposals of such criteria were stated in terms of sparse coding in [6], already discussed in [3], and more recently in [7], or in terms of Bayesian learning [8].

Fig. 1
figure 1

Distribution of receptive profiles in terms of their shapes. Figure extracted from [3]

In this paper, we will focus on the task of position and orientation detection, and propose theoretical motivations based on the uncertainty principle for the corresponding geometry to explain such confinement. In general, the uncertainty principle is indeed a tool that gives information on the possible localization of functions with respect to competing symmetries that in this case are those of the well-known group of translations and rotations of the Euclidean plane. The role of symmetries in the mechanisms of visual perception in V1 is a well recognized point [911], as well as the uncertainty principle was already invoked to explain relevant cortical morphologies [12, 13]. Here, we will use such concepts to characterize the resolution that can be obtained with joint spatial and angular measurements, based on the localization properties of receptive profiles. In terms of such characterizations, we will deduce the bounds observed in Fig. 1 as the result of intrinsic notions of balance between joint measurements resolutions.

2 Receptive Profiles and Relevant Symmetries

We will assume isotropic Gaussian Gabor filters as a model for standard V1 simple cells classical receptive profiles, defined on the Euclidean image plane:

ψ q , p σ (x)= 1 σ π e i p ( x q ) e | x q | 2 / ( 2 σ 2 ) ,x R 2
(1)

with parameters q R 2 , p R 2 , σ R + . Each V1 simple cell is assumed to perform a linear filtering with a function shaped as in (1), so that it can be characterized by these parameters. Their mapping on the two dimensional cortical layers are referred to as cortical maps [14]. In particular, the centers q of receptive fields are in a so-called retinotopic correspondence on the cortex [1], while the size σ is in average larger at the periphery and smaller close to the fovea [15]. The frequency parameters p are generally considered in polar coordinates p=|p|(cosθ,sinθ), where |p| is called spatial frequency and the angle θ up to a factor of π is called preferred orientation, and their cortical maps are also well studied [11, 13, 16].

The family of functions (1) were proposed in [12] due to their optimal localization in space and frequency with respect to the classical Heisenberg uncertainty principle, and their fitness to model the linear behavior of simple cells was thoroughly tested [3]. We note, however, that here we are dealing with a simplified model of isotropic receptive fields, since as we will see this provides enough information for the present study, with the advantage that the results can be stated in a clearer form. We also recall that the real and imaginary parts in (1) correspond to so-called even and odd cells

ψ q , p σ ( x ) = 1 σ π cos ( p ( x q ) ) e | x q | 2 / ( 2 σ 2 ) + i 1 σ π sin ( p ( x q ) ) e | x q | 2 / ( 2 σ 2 )

but it will be sufficient for our purposes to deal with the full complex function as a whole.

2.1 Groups of Transformations

Let us introduce the following unitary operators on L 2 ( R 2 ):

  1. (i)

    translations: T q f(x)=f(xq), q R 2 ,

  2. (ii)

    modulations: M p f(x)= e i p x f(x), p R 2 ,

  3. (iii)

    dilations: Σ σ f(x)= 1 σ f( x σ ), σ R + ,

  4. (iv)

    rotations: R α f(x)=f( r α x), α S 2 ,

where r α stands for the usual counterclockwise rotation of an angle θ on the Euclidean plane. In particular, we note that rotations commute with dilations, and it is easy to see that

R α M p = M r α p R α .
(2)

If we denote with g 1 , a L 2 ( R 2 ) normalized isotropic Gaussian with unit standard deviation

g 1 (x)= 1 π e | x | 2 / 2

then we can characterize the functions (1) in terms of the operators (i), (ii), and (iii) as

ψ q , p σ (x)= T q M p Σ σ g 1 (x).

Such a family is the prototype of a so-called wave packet systems [17], and much is known about these structures [18, 19].

In this work, we will deal with the localization properties of (1) with respect to translations and local rotations, i.e., making use of the symmetries (i) and (iv), since they constitute two fundamental symmetries related to the mechanisms of visual perception in V1 (see, e.g., [9] and references therein).

Local rotations are defined by

R α q f(x)f ( ( r α q ) 1 x ) = T q R α T q f(x),
(3)

where r α q x= r α (xq)+q is a rotation of the Euclidean plane around point q, and with respect to these transformations we have the following.

Lemma 2.1 Let ψ q , p σ be as in (1). Then

R α q ψ q , p σ (x)= ψ q , r α p σ (x).
(4)

Proof Using (2) and the definition of local rotations (3), we get

R α q ψ q , p σ (x)= T q R α T q T q M p Σ σ g 1 (x)= T q R α M p Σ σ g 1 (x)= T q M r θ p R α Σ σ g 1 (x)

so (4) follows since rotations commute with dilations and g 1 is isotropic, i.e., R α g 1 (x)= g 1 (x). □

Actually, the fact that g 1 is isotropic allows to write the whole family (1) in terms of all the operators (i) to (iv). Indeed, denoting with θ the polar angle of p, that means p=|p|(cosθ,sinθ), we can write (1) as

ψ q , p σ (x)= T q R θ M ( | p | 0 ) Σ σ g 1 (x),

where M ( | p | 0 ) f(x)= e i | p | x 1 f(x), so another way to characterize the system of functions (1) is to consider a family { g σ , | p | = M ( | p | 0 ) Σ σ g 1 (x),|p|,σ R + } and rotate and translate each of its members. The aim of next section is actually to deduce properties on the localization of ψ q , p σ with respect to the parameters q and θ, expressed in terms of the parameters |p| and σ.

3 Measures of Uncertainty

In this section, we characterize the uncertainty associated to joint measurements of positions and local orientations in terms of the properties of the measurement devices, expressed by L 2 functions, and quantify such uncertainties for the case of receptive profiles.

We recall that the generators P j of translations along the Cartesian axis are given by partial derivatives

P j d d q j | q = 0 T q = x j ,j=1,2
(5)

while the generator J q of a rotation around point q can be written in terms of the ordinary infinitesimal rotation operator J= x 2 x 1 x 1 x 2

J q d d α | α = 0 R θ q = T q J T q
(6)

and acts as the skew self-adjoint operator on L 2 ( R 2 )

J q f(x)= d d α | α = 0 f ( r α 1 ( x q ) + q ) = ( ( x 2 q 2 ) x 1 ( x 1 q 1 ) x 2 ) f(x).

We will measure averages and variances using the standard definitions for operators on L 2 , denoting with , L 2 ( R 2 ) the L 2 ( R 2 ) scalar product and with L 2 ( R 2 ) the associated norm.

Definition 3.1 Let L be a densely defined skew self-adjoint linear operator on D L 2 ( R 2 ). We define its mean value over fD as

E f (L) ( i L ) f , f L 2 ( R 2 ) R
(7)

and its variance over fD as

( Δ f L ) 2 E f ( ( L E f ( L ) ) 2 ) = ( L E f ( L ) ) f L 2 ( R 2 ) 2 .
(8)

Since skew self-adjoint operators are the infinitesimal generators of a one parameter group of unitary transformations, the meaning of the average (7) is that of measuring the deformation of f under such transformations

E f (L)=i lim t 0 ( exp t L ) f , f L 2 ( R 2 ) f , f L 2 ( R 2 ) t

and the imaginary constant is merely a convention to ensure the result to be real. With this averaging, the variance (8) has the usual meaning of strength of the fluctuations of f under the considered transformations that corresponds to the second moment of the distribution t ( exp t L ) f , f L 2 ( R 2 ) . This means then that the variance (8) provides a measure of the localization of f with respect to the symmetry exptL.

When applied to the operators (5) and (6) of linear and rotational derivatives, these variances correspond respectively to a measure of linear and rotational fluctuations of a function f. The more f is insensitive to translations (f smooth and close to a constant function), the smaller is its P variance, while a small J q variance means that f has little sensitivity to rotations around q.

The notion of localization in orientation that arises indicates that a function consisting of a set of parallel stripes, independently on their widths, is maximally localized in orientation, while a function that is circular symmetric around q is minimally localized.

If we are interested in the joint localization properties of a function with respect to a two parameters group of unitary transformations, generated by two skew self-adjoint operators L 1 and L 2 , we are led to consider the distribution

( t 1 , t 2 ) ( exp t 1 L 1 ) ( exp t 2 L 2 ) f , f L 2 ( R 2 ) .
(9)

In this case, if the operators L 1 and L 2 do not commute, then the second moments of the distribution (9) are influenced by their commutator. Such an effect of competing symmetries is quantified by the uncertainty principle.

3.1 The SE(2) Uncertainty Principle

The operators (5) and (6) satisfy the commutation relations of angular momentum [20]

[ J q , P 1 ] = P 2 ; [ J q , P 2 ] = P 1 .
(10)

These commutators define the algebra of the SE(2) group (see, e.g., [9] and references therein), and for them the following generalized uncertainty principle holds [13, 21], with respect to the quantities of Definition 3.1. Since we are dealing with densely defined operators, we will skip in what follows the technicalities related to operator domains, and refer the statements simply to L 2 ( R 2 ). For more details, see [21].

Theorem 3.2 (SE(2) uncertainty principle)

For any f L 2 ( R 2 ), it holds

{ ( Δ f J q ) ( Δ f P 1 ) 1 2 | E f ( P 2 ) | , ( Δ f J q ) ( Δ f P 2 ) 1 2 | E f ( P 1 ) | .
(11)

These inequalities play the same role for the noncommutative symmetries of rotations and translations as the one played by the ordinary uncertainty inequality for the noncommutativity of quantum mechanical operators. The main difference is that in this case if we consider separately each of the two inequalities, we cannot obtain a constant lower bound. Indeed for a function f the product of variances of an infinitesimal rotations and a translations along one axis can be arbitrarily small, provided that the average of translations along the other axis on f is small. This effect disappears when we consider translations on both axis, which is natural whenever we do not want to discriminate one direction over the other. In this case, we can actually recast the two inequalities (11) into one inequality with a constant lower bound.

The following definition is closely related to that of [22], and for this reason we use the same notation Angv.

Definition 3.3 Let us define the functional

Angv[f] Δ f P E f ( P ) ,

where

E f (P) ( E f ( P 1 ) ) 2 + ( E f ( P 2 ) ) 2 and Δ f P ( Δ f P 1 ) 2 + ( Δ f P 2 ) 2 .

We denote with ΔΘ[f] the corresponding measure of angular uncertainty

ΔΘ[f]arctan ( Angv [ f ] ) .
(12)

With this definition, a direct consequence of the SE(2) uncertainty principle is the following.

Theorem 3.4 For all f L 2 ( R 2 )

( Δ f J q ) Angv[f] 1 2 .
(13)

This inequality resembles the ordinary Heisenberg uncertainty inequality, since the presence of a constant lower bound provides a clear constraint on the joint localizations quantified by Δ f J q and Angv[f]. However, as first noted in [23], the SE(2) uncertainty inequalities (11) cannot be simultaneously minimized, so also (13) does not admit minimizers. This is related to the issue of nonexistence of a canonically conjugate observable for angular momentum [24, 25]. Indeed, if we had a well defined self-adjoint operator canonically commuting with angular momentum, we would end up with a well-known complex equation defining minimal uncertainty states [21], while in this case we have two such equations, whose solutions provide CR function functions on the R 2 × S 1 for two noncompatible almost complex structures [26].

3.2 SE(2) Autocorrelations

We pass now to the study of the properties of the distribution (9) applied to the symmetries under study, that we call autocorrelation since it has the form of the autocorrelation of a function with respect to the group of rotations and translations, and extends naturally the ordinary definition of autocorrelation with respect to translations. We will actually restrict the analysis to the square modulus of correlations, since as we will see it contains enough information for the present purposes. In particular, we will show that such correlations can be used to characterize the uncertainty in the detection of position and local preferred angle associated to a function.

Definition 3.5 Given f in L 2 ( R 2 ), we define its SE(2) autocorrelation centered at q as

C q [f](ξ,α)=| T ξ R α q f , f L 2 ( R 2 ) | 2 .
(14)

In general, C q [f](ξ,α) provides a natural way to study the joint localization properties of f with respect to position and local preferred angle. Indeed, when we specialize to translations we get the usual autocorrelation, and by Plancherel theorem

C q [f](ξ,0)=| R 2 f(xξ) f ( x ) ¯ dx | 2 =|F ( | F f | 2 ) (ξ) | 2
(15)

so we have that by Young inequality and the Riemann–Lebesgue lemma C q [f](ξ,0) is bounded and goes to 0 as ξ becomes large. Moreover, by the usual uncertainty principle, we have that when f is well localized in space, then Ff is broadly localized, hence passing under another Fourier transform C q [f](ξ,0) will decay rapidly, uniformly on q, and vice versa.

On the other hand, if we consider correlations only with respect to rotations, for simplicity centered at q=0

C 0 [f](0,α)=| R 2 f( r α x) f ( x ) ¯ dx | 2
(16)

essentially the same argument applies to the decay of correlations for functions that are localized with respect to rotations.

Remark 3.6 (What does “essentially the same argument” mean)

Since F R θ f= R θ Ff, we get

C 0 [f](0,α)=| R 2 f( r α x) f ( x ) ¯ dx | 2 =| R 2 R α (Ff)(k) F f ( k ) ¯ dk | 2

so setting polar coordinates, with the notation ϕ κ (φ)=(Ff)(κcosφ,κsinφ)

C 0 [f](0,α)= R + κdκ 0 2 π ϕ κ (φα) ϕ κ ( φ ) ¯ = R + κdκ n Z e 2 π i n α | ϕ ˆ κ (n) | 2 ,

where ϕ ˆ κ (n)= 0 2 π e 2 π i n φ ϕ κ (φ)dφ and the last transition is Parseval identity.

Since L 2 ( R 2 ) L 2 ( R + ,κdκ) L 2 ( S 1 ) as a tensor product of Hilbert spaces, and since f is localized with respect to rotations in the real plane if and only if it is localized with respect to rotations in the Fourier plane, then we can assume without loss of generality that ϕ κ (φ)=r(κ)Φ(φ), where Φ(φ) decays rapidly away from φ=0 and R + | r ( κ ) | 2 κdκ=c<. So,

C 0 [f](0,α)=c n Z e 2 π i n α | Φ ˆ (n) | 2

and now strictly the same argument used for (15) applies.

3.3 Uncertainty Associated to SE(2) Measurements with Receptive Profiles

When specialized to receptive profiles, the introduced uncertainties can be explicitly computed. In the proofs, we will use the shorthand notation g σ = Σ σ g 1 , and θ will be the polar angle of p.

Lemma 3.7 The variance of the operators (5) on receptive profiles (1) is

Δ ψ q p σ P j = 1 2 σ .

Proof Since x j ψ q p σ (x)=(i p j ( x j q j ) σ 2 ) ψ q p σ (x), we get E ψ q p σ ( P j )=i p j , and

E ψ q p σ ( P j 2 ) = P j ψ q p σ L 2 ( R 2 ) 2 = R 2 | p j + i ( x j q j ) σ 2 | 2 | g σ ( x q ) | 2 d x = p j 2 + 1 π σ 2 R 2 y j 2 e y 2 d y = p j 2 + 1 2 σ 2 .

 □

Lemma 3.8 The variance of the operator (6) on receptive profiles (1) is

ΔJ Δ ψ q p σ J q = | p | σ 2
(17)

and we will call it angular momentum variance.

Proof Since

x j ψ q p σ (x)= ( i p j ( x j q j ) σ 2 ) ψ q p σ (x)

then

J q ψ q p σ (x)=i ( ( x 2 q 2 ) p 1 ( x 1 q 1 ) p 2 ) ψ q p σ (x).

Its mean value vanishes on ψ q p σ , due to the isotropy of g σ :

E ψ q p σ ( J q ) = J q ψ q p σ , ψ q p σ L 2 ( R 2 ) = i R 2 ( x 2 p 1 x 1 p 2 ) | g σ ( x ) | 2 d x = i | p | R 2 ( r θ x ) 2 | g σ ( x ) | 2 d x = i | p | R 2 x 2 | g σ ( x ) | 2 d x = 0 .

To compute the variance, by analogous arguments

( Δ ψ q p σ J q ) 2 = | p | 2 R 2 ( ( r θ x ) 2 ) 2 | g σ ( x ) | 2 d x = | p | 2 ( σ π ) 2 ( R e x 1 2 / σ 2 d x 1 ) ( R x 2 2 e x 2 2 / σ 2 d x 2 ) = | p | 2 σ 2 2 π R y ( 2 y e y 2 ) d y = ( | p | σ 2 ) 2 .

 □

We have then obtained the following proposition, which shows that for receptive profiles the angular momentum variance ΔJ is inversely proportional to the angular uncertainty quantified in terms of Angv.

Proposition 3.9 Let ψ q p σ be as in (1). Then

ΔJAngv [ ψ q p σ ] = 1 2 .

Proof Using Definition 3.3, Lemma 3.7, and Lemma 3.8, we have

(Angv) [ ψ q p σ ] = 1 σ | p | = 1 2 Δ J .
(18)

 □

We will now consider SE(2) autocorrelations of receptive profiles, and see that they indeed contain precisely the desired joint information on localizations in space and local orientation associated to the uncertainties we computed.

Proposition 3.10 Let ψ q p σ be defined by (1). Then its SE(2)-autocorrelation reads

C q [ ψ q p σ ] (ξ,α)= e | ξ | 2 / ( 2 σ 2 ) e ( Δ J ) 2 ( 1 cos α ) .
(19)

Proof By Lemma 2.1, and computing the Fourier transform of a Gaussian

T ξ R α q ψ q , p σ , ψ q , p σ L 2 ( R 2 ) = T ξ ψ 0 , r α p σ , ψ 0 , p σ L 2 ( R 2 ) = 1 σ 2 π R 2 e i ( r α p ) ( x ξ ) e | x ξ | 2 / ( 2 σ 2 ) e i p x e | x | 2 / ( 2 σ 2 ) d x = e i ( r α p ) ξ σ 2 π e | ξ | 2 / ( 4 σ 2 ) R 2 e i ( p ( r α p ) ) x e | x ξ / 2 | 2 / ( σ 2 ) d x = e i ( ( p + ( r α p ) ) / 2 ) ξ e | ξ | 2 / ( 4 σ 2 ) e σ 2 | p ( r α p ) | 2 / 4

so the result follows since |p( r α p) | 2 =2|p | 2 (1cosα). □

This proposition shows that the decay of the autocorrelation in space is a Gaussian with the same width of the corresponding receptive profile that characterize spatial uncertainty. With respect to rotations, we have ended up with a Von Mises distribution in orientations. Such distributions appear naturally when discussing the SE(2) uncertainty principle [20, 27], but they also provide a good model for orientation tuning of simple cells [28], which is defined as the response curve of a cell to oriented stimuli [28, 29]. This confirms that the introduced notion of localization is compatible with the resolution of measurements performed with receptive profiles. Moreover, we note that the commonly used circular variance [30] of the Von Mises distribution in (19), up to a normalization constant, is

CircVar(ΔJ)=1 I 1 ( ( Δ J ) 2 ) I 0 ( ( Δ J ) 2 )

which results to be numerically close to what we have introduced as angular uncertainty (12) when applied to receptive profiles (18)

( Δ Θ [ ψ q , p σ ] ) 2 = ( arctan ( 1 2 Δ J ) ) 2 .

In particular, as we will see in next section, typical values of ΔJ in the filters encountered in V1 are around 1.7, where the difference between these two notions of variance is around 5 10 2 .

4 Bounds on the Shape Index Induced by Uncertainty

In this section, we will use the measures of uncertainty referred to receptive profiles (18) and (19) to deduce relevant features about the physiological data measured in [3] and depicted in Fig. 1. In particular, we will see how the information provided by the analysis of uncertainty relations of Sect. 3 are sufficient to establish bounds on the number of subregions observed in the family of filters implemented in V1, and permit to reobtain characteristic sampling rates commonly used in image analysis.

A receptive profile ψ q p σ consists of an oscillation of frequency ν= | p | 2 π under a Gaussian bell of width σ, so it appears natural to define a dimensionless index of shape [3]

n=νσ.
(20)

This quantity is related to the number of subregions defining a receptive profiles, since if we let N k be the number of half wavelength of receptive profile’s oscillation within k standard deviations σ, we obtain N k =4kn. As it is apparent from the data measured in [3], we see that approximately k=2 standard deviations are sufficient to represent the main content of the filters, so that we can relate the effective subregions N to n as N=8n.

In terms of n, the angular momentum variance (17) of a receptive profile ψ q p reads

ΔJ= 2 πn
(21)

while its angular variance (12), after (18), reads

ΔΘ [ ψ q p σ ] =arctan ( 1 2 π n ) .
(22)

4.1 Lower Bound for Orientation Measurements

By the discussion in Sect. 3, we have seen how we can quantify with ΔΘ the angle resolution allowed by a linear filtering. If we refer to the task of orientation detection, we can set as a reasonable bound that of angle uncertainty less than π/2, that is expressed by

ΔΘ [ ψ q p σ ] π 4 .

This condition can be stated in terms of the shape index using (22)

n 1 2 π 0.16 n min .

As we can see in Fig. 1 and by the discussions in [3], cells which show a selectivity in orientation all lie above this threshold. Moreover, we note that for indexes n< n min , it can be a hard task to distinguish an even cell from being represented only by a Gaussian, while odd cells under this threshold all appear identical up to a multiplicative factor, so the parametric fit of the Gabor model (1) is quite delicate in this region. We can then interpret the bunch of broadly tuned cells around the zero value of the shape index n as generally below the minimal uncertainty bound that allows a consistent detection of orientations.

4.2 Upper Bound

In order to discuss the upper bound, we introduce a notion of characteristic length associated to a specific level set of the correlations (1), intrinsically related to the task of detection of positions and local orientations. Its purpose is to quantify the minimum distance that one needs to cover in order to decorrelate a function f as much as f is decorrelated when compared at orthogonal directions.

Definition 4.1 The correlation length for f L 2 ( R 2 ) is the smallest distance λ for which

C[f](ξ,0)C[f] ( 0 , π 2 ) |ξ|=λ.
(23)

If we apply this notion to receptive profiles (1), we obtain the following.

Proposition 4.2 The shape index (20) is bounded from above by the ratio of the correlation length λ and the spatial uncertainty σ

2πn λ σ .
(24)

Proof Condition (23) on receptive profiles ψ q p σ , by (19) reads ΔJ 2 2 λ σ , since

e λ 2 / ( 2 σ 2 ) e Δ J 2 λ 2 2 σ 2 Δ J 2

so (24) follows by the relation (21) between ΔJ and the shape index. □

On the other hand, as discussed when dealing with the relation between the shape index and the number of subregions, we have also that the effective field of influence of a receptive profile can be set within two standard deviations σ. From this point of view, we can then assume that the distance d at which a receptive profile is effectively spatially uncorrelated corresponds to the distance that one has to cover in order to let its effective field of influence not intersect with its translation at a distance d, i.e., d=4σ.

In order to couple with both position and orientation measurements, we will then consider the hypothesis of balance of the two characteristic scales introduced that is the identification λ=d. By (24), this condition can be stated in terms of the shape index as

n 2 π 0.64 n max .

To compare this bound with Fig. 1, we recall that here we are dealing with the simplified model of isotropic receptive fields, while in [3] the analysis is performed considering two anisotropic indexes n x and n y . In terms of such indexes, we can see that the largest part of the population lies within two bounds n x 0.5 and n y 0.76, and n max looks in good accordance with their mean value.

The question of whether this identification of characteristic distances is truly implemented in the cortex cannot be answered at this point, but we note that a cortical scale related to the symmetries under study that is possibly compatible with the proposed relation is the mean correlation length of orientation preference maps (see, e.g., [13, 14] and references therein). Indeed, by the measurements performed in [31], we see that such scale is comparable with the size of a so-called cortical point image, that is, the cortical region that is activated after a highly spatially localized stimulus, and at least when we reduce to linear behavior of cells this notion corresponds to what we have indicated as effective field of influence.

4.3 Sampling on Orientations

Another intriguing consequence of the performed uncertainty analysis can be stated in terms of optimal sampling rates for orientation detection. Indeed, if we consider the mean value on shape index measured in [3], or equivalently, in terms of the deduced bounds, for n= n max + n min 2 0.4, we have that

ΔΘ [ ψ q p σ ] =arctan ( 1 0.8 π ) π 8 .

With respect to Gabor filters possessing such n, one way to use such result is to consider that the detection of orientations at angles that are closer than this uncertainty do not provide an actual improvement in the resolution of the local orientation present in the stimulus, so that it can be sufficient to cover the interval of orientations [0,π) with a sampling having a π/8 spacing. This actually compares well with the notions of optimal sampling adopted in image analysis tasks (see, e.g., in [32] and references therein), generally justified with independent arguments. Moreover, this uncertainty analysis permits to set clear sampling spacings depending on the shape index of the filter used.

5 Conclusions

In this paper, we have studied theoretical aspects of an analytic characterization of uncertainty that generalizes the well-known Heisenberg uncertainty principle to the symmetries associated with the task of joint measurements of position and local orientation. The implications of this analysis, together with an hypothesis of balance between characteristic correlation distances, allowed us to obtain bounds comparable with experimental data on the shape index of the V1 simple cells that are selective for orientation, and to separate them from broadly tuned cells, which lie below the uncertainty bound for consistent orientation detection.

We remark that this was possible even if our working assumptions on the functional behavior of simple cells were reduced to linear filtering with symmetric receptive fields, and the only considered task is the one associated to the sole symmetries of rotations and translations.

Whether such elementary principles could be directly responsible of the observed distribution of receptive profiles is a question that can hardly find an answer. Nevertheless, the present study shows that they are sufficient to describe many of the relevant features that concern the shape of simple cells.

References

  1. Hubel DH, Wiesel TN: Functional architecture of macaque monkey visual cortex. Proc R Soc Lond B 1977, 198: 1–59. 10.1098/rspb.1977.0085

    Article  Google Scholar 

  2. Graham NV: Beyond multiple pattern analyzers modeled as linear filters (as classical V1 simple cells): useful additions of the last 25 years. Vis Res 2011, 51: 1397–1430. 10.1016/j.visres.2011.02.007

    Article  Google Scholar 

  3. Ringach DL: Spatial structure and symmetry of simple cell receptive fields in macaque primary visual cortex. J Neurophysiol 2002, 88: 455–463.

    Google Scholar 

  4. Jones JP, Palmer LA: The two-dimensional spatial structure of simple receptive fields in cat striate cortex. J Neurophysiol 1987, 58: 1187–1211.

    Google Scholar 

  5. Jones JP, Palmer LA: An evaluation of the two-dimensional Gabor filter model of simple receptive fields in cat striate cortex. J Neurophysiol 1987, 58: 1233–1258.

    Google Scholar 

  6. Olshausen BA, Field DJ: Sparse coding with an overcomplete basis set: a strategy employed by V1? Vis Res 1997, 37: 3311–3325. 10.1016/S0042-6989(97)00169-7

    Article  Google Scholar 

  7. Shelton J, Sterne P, Bornschein J, Sheikh AS, Luecke J: Why MCA? Nonlinear sparse coding with spike-and-slab prior for neurally plausible image encoding. Adv Neural Inf Process Syst 2012, 25: 2285–2293.

    Google Scholar 

  8. Hosoya H: Multinomial Bayesian learning for modeling classical and nonclassical receptive field properties. Neural Comput 2012, 24: 2119–2150. 10.1162/NECO_a_00310

    Article  MathSciNet  Google Scholar 

  9. Citti G, Sarti A: A cortical based model of perceptual completion in the roto-translation space. J Math Imaging Vis 2006, 24: 307–326. 10.1007/s10851-005-3630-2

    Article  MathSciNet  Google Scholar 

  10. Chossat P, Faugeras O: Hyperbolic planforms in relation to visual edges and textures perception. PLoS Comput Biol 2009, 5: 1–16.

    Article  MathSciNet  Google Scholar 

  11. Keil W, Wolf F: Coverage, continuity and visual cortical architecture. Neural Syst Circuits 2011. 10.1186/2042-1001-1-17

    Google Scholar 

  12. Daugman JG: Uncertainty relation for resolution in space, spatial frequency, and orientation optimized by two dimensional visual cortical filters. J Opt Soc Am A 1985, 2: 1160–1169. 10.1364/JOSAA.2.001160

    Article  Google Scholar 

  13. Barbieri D, Citti G, Sanguinetti G, Sarti A: An uncertainty principle underlying the functional architecture of V1. J Physiol (Paris) 2012, 106: 183–193. 10.1016/j.jphysparis.2012.03.001

    Article  Google Scholar 

  14. Basole A, White LE, Fitzpatrick D: Mapping multiple features in the population response of visual cortex. Nature 2003, 423: 986–990. 10.1038/nature01721

    Article  Google Scholar 

  15. Hubel DH, Wiesel TN: Uniformity of monkey striate cortex: a parallel relationship between field size, scatter, and magnification factor. J Comp Neurol 1974, 158: 295–305. 10.1002/cne.901580305

    Article  Google Scholar 

  16. Ohki K, Chung S, Kara P, Hubener M, Bonhoeffer T, Reid RC: Highly ordered arrangement of single neurons in orientation pinwheels. Nature 2006, 442: 925–928. 10.1038/nature05019

    Article  Google Scholar 

  17. Córdoba A, Fefferman C: Wave packets and Fourier integral operators. Commun Partial Differ Equ 1978, 3: 979–1005. 10.1080/03605307808820083

    Article  MATH  Google Scholar 

  18. Kalisa C, Torrésani B: N-dimensional affine Weyl–Heisenberg wavelets. Ann Inst Henri Poincaré. Phys Théor 1993, 59: 201–236.

    MATH  Google Scholar 

  19. Labate D, Weiss G, Wilson E: An approach to the study of wave packet systems. Contemporary Mathematics 345. In Wavelets, Frames and Operator Theory. Birkhäuser, Basel; 2004:215–235.

    Chapter  Google Scholar 

  20. Carruthers P, Nieto MM: Phase and angle variables in quantum mechanics. Rev Mod Phys 1968, 40: 411–440. 10.1103/RevModPhys.40.411

    Article  Google Scholar 

  21. Folland GB, Sitaram A: The uncertainty principle: a mathematical survey. J Fourier Anal Appl 1997, 3: 207–238. 10.1007/BF02649110

    Article  MATH  MathSciNet  Google Scholar 

  22. Breitenberger E: Uncertainty measures and uncertainty relations for angle observables. Found Phys 1985, 15: 353–364. 10.1007/BF00737323

    Article  MathSciNet  Google Scholar 

  23. Jackiw R: Minimum uncertainty product, number phase uncertainty product and coherent states. J Math Phys 1968., 9: Article ID 339 Article ID 339

    Google Scholar 

  24. Dubin DA, Hennings MA, Smith TB: Quantization in polar coordinates and the phase operator. Publ Res Inst Math Sci 1994, 30: 479–532. 10.2977/prims/1195165908

    Article  MATH  MathSciNet  Google Scholar 

  25. Kastrup HA: Quantization of the canonically conjugate pair angle and orbital angular momentum. Phys Rev A 2006., 73: Article ID 052104 Article ID 052104

    Google Scholar 

  26. arXiv: http://arxiv.org/abs/arXiv:1301.3783 Barbieri D, Citti G: Reproducing kernel Hilbert spaces of CR functions for the Euclidean Motion group. arXiv:1301.3783.

  27. Hradil Z, Řehác̆ek J, Bouchal Z, C̆elechovský R, Sánchez-Coto LL: Minimum uncertainty measurement of angle and angular momentum. Phys Rev Lett 2006., 97: Article ID 243601 Article ID 243601

    Google Scholar 

  28. Swindale NV: Orientation tuning curves: empirical description and estimation of parameters. Biol Cybern 1998, 78: 45–56. 10.1007/s004220050411

    Article  MATH  Google Scholar 

  29. Mooser F, Bosking WH, Fitzpatrick D: A morphological basis for orientation tuning in primary visual cortex. Nat Neurosci 2004, 7: 872–879. 10.1038/nn1287

    Article  Google Scholar 

  30. Jammalamadaka SR, Sengupta A: Topics in Circular Statistics. World Scientific, Singapore; 2001.

    MATH  Google Scholar 

  31. Bosking WH, Crowley JC, Fitzpatrick D: Spatial coding of position and orientation in primary visual cortex. Nat Neurosci 2009, 5: 874–882.

    Article  Google Scholar 

  32. Lee TS: Image representation using 2D Gabor wavelets. IEEE Trans Pattern Anal Mach Intell 1996, 18: 959–971. 10.1109/34.541406

    Article  Google Scholar 

Download references

Acknowledgements

The research of the first author was supported by Grant DIM2011—Région Île de France. The research of the second author was supported by Project AGAPE.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to D Barbieri.

Additional information

Competing Interests

The authors declare that they have no competing interests.

Authors’ Contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Barbieri, D., Citti, G. & Sarti, A. How Uncertainty Bounds the Shape Index of Simple Cells. J. Math. Neurosc. 4, 5 (2014). https://doi.org/10.1186/2190-8567-4-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/2190-8567-4-5

Keywords