29.03.04 · psychology / sensation-perception

Hubel and Wiesel's visual cortex architecture: orientation columns, ocular dominance, and hypercolumns

shipped3 tiersLean: none

Anchor (Master): Hubel & Wiesel primary literature 1959-1977; Bonhoeffer & Grivvald 1991 (Nature 353); Bonhoeffer & Grivenvald 1993 (Nature 363); Stryker & Harris 1986 (J. Neurosci. 6); Ohki, Chung, Ch'ng, Kara & Reid 2005 (Nature 433); Lien & Scanziani 2013 (Nature 493)

Intuition Beginner

When you look at a tree branch, your eye sends a raw, pixel-like image to the back of your brain. But the brain does not see pixels. It sees edges, lines, and motion. David Hubel and Torsten Wiesel, working at Harvard from 1958, discovered how. They inserted microelectrodes into the visual cortex of cats (and later macaques) and listened to individual neurons fire. Most neurons stayed silent until a bar or an edge appeared at one specific angle in one specific spot — and then the cell erupted.

These orientation-tuned neurons are called simple cells. The deeper finding was architectural. The cortex is organised as columns: every neuron stacked in one vertical column prefers the same angle, and neighbouring columns prefer slightly different angles. Walk an electrode sideways across the cortical surface and the preferred angle rotates smoothly, sweeping through the full every to micrometres. The cortex lays out orientation as a spatial map, with nearby neurons encoding similar features.

This work, recognised by the 1981 Nobel Prize in Physiology or Medicine (shared with Roger Sperry), established that sensory cortex is not a passive relay. It is an active feature extractor with hierarchical, columnar architecture — the founding insight of systems neuroscience. The brain builds its picture of the world one edge at a time, and the cortical surface is the canvas on which that building is laid out.

Visual Beginner

The picture shows the cortical surface as a tiled mosaic of hypercolumns, each about one square millimetre. Within each hypercolumn, two axes of organisation cross at right angles: orientation columns that sweep through to , and ocular-dominance stripes that alternate between the left and right eyes. The orientation map bends into pinwheels at the centre of each column, where every angle meets at a single point.

The "ice-cube model" of Stryker and Harris pictures orientation as one face of the cube and ocular dominance as the adjacent face: two independent feature axes laid out on the same patch of cortex. Real cortex is messier — orientation preference bends into pinwheels rather than lying in straight stripes — but the principle holds. Each hypercolumn handles one small patch of the visual field and represents that patch at every orientation and through both eyes.

Worked example Beginner

A bar of light moves across a cat's visual field at an angle of from vertical. A microelectrode, its tip inside the cat's primary visual cortex (V1, at the back of the brain), is positioned next to a single neuron. For most orientations the cell is quiet, firing at perhaps spikes per second (its baseline rate). The moment the bar enters the cell's small receptive field at exactly , the firing rate leaps to about spikes per second.

Step 1: rotate the bar to . The cell falls silent, dropping back to its baseline of spikes per second. The cell is selective for one specific orientation, not just any edge.

Step 2: lower the electrode micrometres straight down through the cortex. A new neuron comes into view. Show the bar at — it fires at spikes per second. The neuron immediately below prefers , not . Orientation preference is laid out as a function of position in the cortex.

Step 3: lower another micrometres. The next neuron prefers . Then one that prefers (vertical). Then one at . Walking the electrode through the cortex rotates the preferred angle through the full over about micrometres of travel.

What this tells us: the cortex uses spatial position to encode features. A vertical electrode finds neurons that all share one preference; a horizontal traversal walks through angle space. The cortical surface is a physical map of orientation.

Check your understanding Beginner

Formal definition Intermediate+

The primary visual cortex (V1), also called Brodmann area 17 or striate cortex, is the cortical region that receives input from the lateral geniculate nucleus (LGN) of the thalamus and is the first cortical stage of visual processing. Its architecture has eight interlocking features, established in cat by Hubel and Wiesel (1962) [HubelWiesel1962] and in primate by Hubel and Wiesel (1968) [HubelWiesel1968].

Definition (Retinotopic map). V1's surface maps the contralateral visual field. neighbouring points in the visual field project to neighbouring points on V1's surface, with strong foveal magnification: roughly half of primate V1 represents the central to of vision. The cortical magnification factor , defined as millimetres of cortex per degree of visual angle at eccentricity , obeys approximately with and in macaque.

Definition (Lamination). V1 has six layers (numbered through from pia to white matter). Afferents from the LGN terminate primarily in layer 4, with the magnocellular pathway arriving in layer and the parvocellular pathway in layer . Simple cells are concentrated in layers and ; complex cells in layers and ; output to higher cortical areas (V2, V4, MT) leaves from layers , , and .

Definition (Simple cell). A simple cell has an oriented receptive field with spatially distinct, antagonistic ON and OFF subregions that can be mapped by small spots of light [HubelWiesel1959]. It responds optimally to an elongated stimulus (a bar or edge) of the right orientation at the right position. Its response is approximately linear in the stimulus: the firing rate is well predicted by convolving the image with the cell's receptive-field kernel.

Definition (Complex cell). A complex cell is orientation-tuned and typically binocular, but lacks separable ON and OFF subregions. It responds to the right orientation anywhere within a larger receptive field — its response is approximately position-invariant within that field. The firing rate is approximately a nonlinear function (typically half-wave rectification) of a weighted sum of simple-cell responses at offset positions.

Definition (Hypercomplex / end-stopped cell). A hypercomplex cell (also called end-stopped) responds selectively to line endings, corners, or stopped edges. It is built by convergence from multiple complex cells with offset receptive fields of the same orientation, so that simultaneous activation by a long line produces mutual inhibition.

Definition (Orientation column). An orientation column is a vertical slab of cortical tissue, roughly deep, whose neurons share a common orientation preference. As one moves tangentially across the cortical surface, the preferred orientation rotates smoothly, completing the full cycle over a distance of to (the hypercolumn width).

Definition (Ocular-dominance column). An ocular-dominance (OD) column is a vertical slab whose neurons preferentially respond to input from one eye. OD columns for the left and right eye alternate as parallel stripes roughly wide across the cortical surface, prominent in layer 4 where they originate from the eye-segregated LGN inputs. In the Stryker-Harris "ice-cube model" [StrykerHarris1986], orientation and ocular dominance are orthogonal axes on the cortical surface.

Definition (Hypercolumn). A hypercolumn is the roughly block of cortical surface containing the full set of orientation columns (covering to ) plus both ocular-dominance stripes for a single visual-field location [HubelWiesel1977]. It is the basic repeating unit of V1: roughly hypercolumns tile primate V1, one per small patch of visual field. Every hypercolumn carries the complete machinery needed to represent one piece of the visual scene.

Counterexamples to common slips

  • Every V1 cell is orientation-tuned. No. Neurons in layer of primates retain the circular centre-surround receptive fields inherited from their parvocellular LGN inputs, and are not orientation-selective. Orientation selectivity emerges downstream of layer 4, in layers , , , and .

  • Pinwheels are an imaging artefact. No. The orientation pinwheel structure that Bonhoeffer and Grivenvald first visualised by optical imaging in 1991 (and 1993) was confirmed at single-cell resolution by Ohki, Reid and colleagues using two-photon calcium imaging in 2005 and 2007. The pinwheel is a real, repeatable spatial feature of V1.

  • Monocular deprivation always shifts ocular-dominance columns. No. The shift is permanent only if the deprivation occurs during the critical period (postnatal weeks to in cats). Deprivation in adulthood produces minimal reorganisation, because the cortical circuitry has stabilised.

  • The Hubel-Wiesel hierarchical model is the actual circuit. No. Lien and Scanziani (2013) [LienScanziani2013] showed that recurrent amplification within V1, not pure feed-forward LGN convergence, is required to generate the observed orientation tuning. The Hubel-Wiesel model was a profound simplification that captured the input-output structure but understated the role of intracortical recurrence.

Key model: the hierarchical convergence architecture Intermediate+

Model (Hubel-Wiesel hierarchical convergence). A simple cell's receptive field is the linear arrangement of LGN centre-surround inputs along an oriented bar, and a complex cell's receptive field is the rectified sum of many same-orientation simple cells at offset positions. The model predicts (i) that simple cells have separable ON and OFF subregions and respond linearly to the stimulus, and (ii) that complex cells are position-invariant within a larger receptive field and respond nonlinearly.

Derivation. (i) Simple cell from LGN inputs. A magnocellular or parvocellular LGN relay cell has a circular centre-surround receptive field whose response to a stimulus is

with the centre-surround kernel for surround scale to . Suppose a simple cell receives convergent excitatory input from a row of LGN cells whose centres lie along an axis at orientation through the position , with the -th LGN afferent at offset for spacing . The simple-cell firing rate is then approximately

with a threshold-linear rectifier. Pulling the sum inside the integral gives an effective simple-cell receptive-field kernel , which is an oriented bar of alternating ON-centre contributions flanked by antagonistic surrounds. This kernel is maximal when the stimulus is a bright bar at angle overlapping the ON row, darkened by the antagonistic flanks. The model therefore predicts:

(a) The cell fires only for a narrow range of orientations around — orientation selectivity.

(b) ON and OFF subregions are spatially separated and can be mapped by small spots — the defining property of a simple cell.

(c) The cell is monocular (it inherits its eye-of-origin from the LGN inputs, which themselves are eye-segregated).

(d) The response is approximately linear in the stimulus, with the predicted firing rate well-fit by .

These four predictions are the empirical signature by which simple cells are identified experimentally.

(ii) Complex cell from simple-cell convergence. Now suppose a target cell receives convergent input from simple cells sharing the same orientation preference but with receptive-field centres at offset positions covering a larger region of diameter . If each simple cell contributes a rectified response , the complex-cell firing rate is

Because the simple cells have offset receptive fields, an edge at orientation placed anywhere within the diameter- region will excite at least one of them, and the outer sum produces a response. The complex cell therefore responds to orientation at any position within its receptive field — position-invariant orientation tuning, the empirical signature of the complex cell. Because each is half-wave rectified (i.e. only positively correlated stimuli drive the cell), the complex cell responds to both ON and OFF edges of the right orientation, again matching experiment: complex cells fire for both light-on-dark and dark-on-light edges of the preferred orientation. Binocularity arises if some of the convergent simple cells carry left-eye input and others right-eye input; the pooled complex cell then responds to either eye, which is what is observed.

(iii) Hypercomplex cell from complex-cell convergence. Extending the construction: if a target cell receives convergent input from a row of complex cells sharing orientation and arranged along the same axis, mutual inhibition from the off-axis complex cells when a long line is presented produces end-stopping: the cell fires vigorously for a short bar that fits within its receptive field but is silenced when the bar extends past the inhibitory flanks. This is the empirical signature of a hypercomplex (end-stopped) cell.

The convergence model therefore predicts the full hierarchy: LGN simple complex hypercomplex, with each stage constructed by summing many same-orientation inputs from the previous stage at offset positions or with offset connectivity. The empirical match is excellent for the input-output behaviour of V1 cells.

Caveat. The model is a feed-forward simplification. Lien and Scanziani (2013) [LienScanziani2013] showed in mouse V1 that blocking intracortical recurrence abolishes orientation tuning even when the thalamocortical (LGN V1) input is preserved, demonstrating that the observed tuning requires recurrent amplification within V1. The Hubel-Wiesel convergence architecture correctly describes the input-output receptive field but understates the contribution of intracortical feedback to the sharpness and robustness of tuning.

Bridge. The hierarchical convergence model builds toward the critical-period plasticity work of Hubel and Wiesel (1970) [HubelWiesel1970] on ocular-dominance column formation, where this is exactly the substrate whose wiring is permanently reshaped by monocular deprivation in the postnatal window. The central insight — that cortex extracts features through a hierarchy of tuned neurons — generalises across every sensory modality and appears again in 18.13.02 for hair-cell mechanotransduction and the cochlear place-frequency map, where a similar labelled-line architecture encodes acoustic frequency. The bridge is between a feed-forward convergence circuit and a self-organising cortical sheet whose recurrent amplification was the structural correction Hubel and Wiesel did not need to make their architecture correct in the large.

Exercises Intermediate+

Advanced results Master

Result 1 (Hubel-Wiesel 1959: discovery of orientation selectivity). In their first V1 paper, Hubel and Wiesel [HubelWiesel1959] reported that neurons in cat striate cortex fire selectively for the orientation of an elongated stimulus within a small receptive field. The discovery was famously precipitated by an accidental observation: a slide slipped in the projector, producing an edge rather than the spot stimuli the experimenters had intended, and the cell they were recording from suddenly fired vigorously. The cells we now call simple cells were characterised as oriented, monocular, with separable ON and OFF subregions mappable by small spots.

Result 2 (Hubel-Wiesel 1962: the simple / complex / hypercomplex hierarchy). In a long J. Physiol. paper [HubelWiesel1962], Hubel and Wiesel introduced the complex cell (oriented, position-invariant, binocular) and the hypercomplex cell (end-stopped), proposed the hierarchical convergence architecture (simple LGN bar; complex many simple; hypercomplex many complex), and inferred the existence of orientation columns from oblique electrode penetrations in which orientation preference shifted smoothly with depth.

Result 3 (Hubel-Wiesel 1968: macaque V1 architecture). Extending to primate [HubelWiesel1968], Hubel and Wiesel characterised roughly cells in macaque striate cortex and established that the columnar organisation generalises to primates. The paper fixed the layering scheme (simple cells in layers and , complex in and ) and identified the eye-segregated termination of LGN inputs in layer as the structural substrate of ocular-dominance columns.

Result 4 (Hubel-Wiesel 1970: critical-period plasticity). In a series of deprivation experiments [HubelWiesel1970], Hubel and Wiesel showed that monocular eyelid closure during a postnatal critical period (weeks to in cats) permanently shifts ocular-dominance columns toward the open eye: the open eye expands its territory in layer 4, the closed eye contracts, and the animal behaves as functionally blind in the deprived eye. Deprivation after the critical period produces minimal reorganisation. This work founded the modern field of critical-period developmental neuroscience.

Result 5 (Hubel-Wiesel 1977: the hypercolumn hypothesis). The Ferrier Lecture to the Royal Society [HubelWiesel1977] crystallised the hypercolumn as the basic repeating unit of V1: a block of cortical surface containing a full set of orientation columns (covering to ) plus both ocular-dominance stripes for a single visual-field location. Roughly hypercolumns tile primate V1. The hypercolumn is the structural unit that makes V1 a complete representation: every visual-field location is analysed at every orientation through both eyes.

Result 6 (Stryker-Harris 1986: the ice-cube model and activity-dependence). Stryker and Harris [StrykerHarris1986] refined the Hubel-Wiesel picture with the ice-cube model, in which orientation preference varies along one axis of the cortical surface and ocular dominance along the orthogonal axis. They further showed that binocular impulse blockade (preventing retinal activity during development) abolishes the formation of ocular-dominance columns, establishing that the columns are shaped by activity-dependent competition rather than specified chemotactically.

Result 7 (Bonhoeffer-Grivenvald 1991, 1993: optical imaging of pinwheels). Bonhoeffer and Grivenvald, using voltage-sensitive-dye optical imaging [BonhoefferGrivenvald1991] in cat area 18 (and the 1993 follow-up in area 17), directly visualised the two-dimensional orientation map across the cortical surface. The map is not a set of parallel stripes but a lattice of pinwheels — singular points around which orientation preference rotates through , with typically four branches per pinwheel. The pinwheel structure could not have been inferred from single-electrode penetrations; optical imaging was the enabling technology.

Result 8 (Ohki-Reid 2005, 2007: two-photon confirmation at single-cell resolution). Ohki, Chung, Ch'ng, Kara and Reid used two-photon calcium imaging to map orientation preference of individual neurons within pinwheel centres in cat and rodent V1. They confirmed that even neurons at the pinwheel singularity — surrounded by cells of all orientations — have sharp orientation tuning, ruling out the hypothesis that pinwheel neurons are unselective. The hypercolumn architecture was thus confirmed at the resolution of single cells.

Result 9 (Lien-Scanziani 2013: the recurrent-circuit revision). Lien and Scanziani [LienScanziani2013], working in mouse V1, showed that optogenetic silencing of intracortical activity abolishes orientation tuning even when thalamocortical LGN inputs are preserved. The Hubel-Wiesel feed-forward convergence architecture captures the structural bias but understates the role of recurrent amplification: sharp, robust orientation tuning requires that the weak feed-forward bias provided by LGN convergence be amplified by intracortical recurrence. The convergence model is the skeleton; recurrence is the gain.

Synthesis. The foundational reason V1 is organised as a tiled mosaic of hypercolumns is that orientation and ocular dominance must be represented at every retinotopic location, and this is exactly the design constraint that produces the repeating module. The central insight of Hubel and Wiesel — that cortex uses hierarchical convergence to build progressively more sophisticated receptive fields — generalises across every sensory area: putting these together with the critical-period plasticity work, the same columnar architecture is both the substrate of feature extraction and the locus of experience-dependent wiring. The bridge is between the single-electrode discovery, the optical-imaging visualisation of pinwheels, and the recurrent-circuit revision, and the pattern recurs in the modern two-photon confirmation at single-cell resolution. The Lien-Scanziani correction identifies the feed-forward convergence skeleton with the structural bias and the recurrent amplification with the gain that makes the tuning sharp.

Full proof set Master

Proposition (Inverse-linear cortical magnification and the foveal over-representation). Assume the cortical magnification factor obeys with in macaque, and that visual eccentricity extends to . Then the fraction of V1 surface area devoted to the central degrees of vision is

For and , , i.e. roughly of primate V1 is devoted to the central of vision, consistent with empirical measurements of foveal magnification.

Proof. The cortical surface area devoted to eccentricity range is proportional to the integral of the local cortical-area density. A small annulus at eccentricity of angular width has visual-field area ; the cortical territory devoted to it is (the accounts for isotropic magnification in two dimensions). Hence

where the constant cancels in the ratio .

The indefinite integral evaluates as follows. Substitute , :

At : . For : .

Numerator (, ): .

Denominator (, ): .

Hence , i.e. about of V1 is devoted to the central of vision. Empirical measurements in macaque V1 (using electrophysiological retinotopic mapping and MRI-based surface reconstruction) report that the central occupies roughly of striate cortex, in close agreement with the formula. The inverse-linear magnification law is therefore not an arbitrary fit: it follows from the geometry of mapping a finite cortical sheet onto a hemispheric visual field with foveal-weighting, and it predicts the observed area budget to within the scatter of biological measurement.

Proposition (Position invariance by hierarchical convergence). Let be a family of oriented simple-cell receptive-field kernels of common orientation , with centres densely covering a region of diameter . Let the complex-cell response be with a non-negative rectifier. Then for any translation with , up to boundary effects at the edges of .

Proof. Let be the translated stimulus. The response of the -th simple cell is

Substituting , :

where is the kernel centred at . By hypothesis, the centres densely cover , so for the shifted centre lies within , and the index corresponds to a genuine member of the convergent family. Therefore , where the boundary terms vanish for translations small compared to . Hence , i.e. position invariance holds up to boundary corrections at the receptive-field edge.

Connections Master

  • Sensation and perception survey 29.03.01 is the immediate prerequisite and parent survey for this unit, which deepens the brief treatment of V1 architecture sketched there into the full Hubel-Wiesel columnar account. The general framework of 29.03.01 — receptor-potential to labelled-line, feature extraction in subcortical and cortical relays — is exactly the framework within which the Hubel-Wiesel hypercolumn sits as the canonical worked example of cortical feature representation.

  • Hair cell mechanotransduction and cochlear tuning 18.13.02 is the canonical sensory-systems depth peer and another Nobel-tier story of sensory encoding. The cochlea's place-frequency (tonotopic) map — a one-dimensional frequency axis laid out along the basilar membrane — is the auditory analogue of V1's two-dimensional orientation-plus-retinotopy map. Both are cases where a sensory epithelium lays out its feature space as a spatial map on the receptor surface, and both illustrate that the brain's feature vocabulary is set by the geometry of the periphery.

  • The cognitive revolution and the computational theory of mind 29.14.02 is the chapter-level narrative into which the Hubel-Wiesel story fits as the empirical pillar. The cognitive revolution's claim — that the mind is an information-processing system whose intermediate representations can be characterised independently of their substrate — receives its strongest single piece of empirical support from the discovery that V1 represents visual features by explicit, mappable cortical codes. The hypercolumn is the cleanest known instance of a representation that is simultaneously neural (it is made of neurons) and computational (it implements a feature decomposition).

  • Adaptive optics, interferometry, and high-angular-resolution astronomy 28.09.02 pending provides the optical-resolution framework that bounds what the eye's photoreceptor mosaic can resolve before the visual cortex ever sees an image. The diffraction limit of the eye's optics, the sampling density of the foveal cone mosaic, and the cortical magnification factor together determine the foveal acuity that V1 represents at such disproportionate cortical cost. The Hubel-Wiesel hypercolumn is the cortical end of a pipeline whose properties are set at every stage by physical limits — optical, photoreceptor, and cortical.

Historical & philosophical context Master

David Hubel and Torsten Wiesel began their collaboration at the Harvard Medical School in 1958, having been assigned adjacent laboratory space by Stephen Kuffler. Their first discovery of orientation selectivity in cat striate cortex [HubelWiesel1959] was famously precipitated by an accidental observation during an experiment intended to map receptive fields with spots of light: a slide slipped in the projector, producing an edge rather than a spot, and the cell they were recording from suddenly fired vigorously. The accident prompted a systematic re-mapping of V1 cells with oriented bars, which became the foundation of the simple-cell characterisation.

The long J. Physiol. paper of 1962 [HubelWiesel1962] introduced the complex cell and the hypercomplex cell, the hierarchical convergence architecture, and the columnar-organisation hypothesis inferred from oblique electrode penetrations. The 1968 paper extended the work to macaque [HubelWiesel1968], characterising roughly cells and fixing the layering scheme. The 1970 paper on the critical period [HubelWiesel1970] showed that monocular deprivation during postnatal weeks to in cats permanently shifts ocular-dominance columns, founding modern critical-period developmental neuroscience and motivating early clinical intervention for congenital cataracts and strabismus. The 1977 Ferrier Lecture [HubelWiesel1977] synthesised twenty years of work into the hypercolumn hypothesis. Hubel and Wiesel shared the 1981 Nobel Prize in Physiology or Medicine with Roger Sperry; Hubel's Nobel lecture remains the canonical first-person account.

The columnar organisation that Hubel and Wiesel inferred from single-electrode penetrations was confirmed and refined byBonhoeffer and Grivenvald, whose 1991 optical-imaging paper [BonhoefferGrivenvald1991] revealed the orientation map as a lattice of pinwheels rather than parallel stripes. The ice-cube model of Stryker and Harris [StrykerHarris1986] had already established that orientation and ocular dominance are orthogonal axes and that activity during development is required for OD column formation. The single-cell-resolution confirmation of pinwheel architecture came from Ohki, Reid and colleagues using two-photon calcium imaging in 2005 and 2007. The modern circuit revision by Lien and Scanziani in 2013 [LienScanziani2013] showed that the feed-forward convergence architecture understates the role of intracortical recurrent amplification, and is the principal correction needed to bring the Hubel-Wiesel model into line with contemporary circuit physiology.

Bibliography Master

@article{HubelWiesel1959,
  author = {Hubel, David H. and Wiesel, Torsten N.},
  title = {Receptive fields of single neurones in the cat's striate cortex},
  journal = {Journal of Physiology},
  volume = {148},
  year = {1959},
  pages = {574--591},
}

@article{HubelWiesel1962,
  author = {Hubel, David H. and Wiesel, Torsten N.},
  title = {Receptive fields, binocular interaction and functional architecture in the cat's visual cortex},
  journal = {Journal of Physiology},
  volume = {160},
  year = {1962},
  pages = {106--154},
}

@article{HubelWiesel1968,
  author = {Hubel, David H. and Wiesel, Torsten N.},
  title = {Receptive fields and functional architecture of monkey striate cortex},
  journal = {Journal of Physiology},
  volume = {195},
  year = {1968},
  pages = {215--243},
}

@article{HubelWiesel1970,
  author = {Hubel, David H. and Wiesel, Torsten N.},
  title = {The period of susceptibility to the physiological effects of unilateral eye closure in kittens},
  journal = {Journal of Physiology},
  volume = {206},
  year = {1970},
  pages = {419--436},
}

@article{HubelWiesel1977,
  author = {Hubel, David H. and Wiesel, Torsten N.},
  title = {Ferrier lecture: Functional architecture of macaque monkey visual cortex},
  journal = {Proceedings of the Royal Society of London B},
  volume = {198},
  year = {1977},
  pages = {1--59},
}

@article{StrykerHarris1986,
  author = {Stryker, Michael P. and Harris, William A.},
  title = {Binocular impulse blockade prevents the formation of ocular dominance columns in cat visual cortex},
  journal = {Journal of Neuroscience},
  volume = {6},
  year = {1986},
  pages = {2117--2133},
}

@article{BonhoefferGrivenvald1991,
  author = {Bonhoeffer, Tobias and Grivenvald, Amiram},
  title = {Iso-orientation domains in cat visual cortex are arranged in pinwheel-like patterns},
  journal = {Nature},
  volume = {353},
  year = {1991},
  pages = {429--431},
}

@article{BonhoefferGrivenvald1993,
  author = {Bonhoeffer, Tobias and Grivenvald, Amiram},
  title = {The layout of iso-orientation domains in area 18 of cat visual cortex: optical imaging reveals a pinwheel-like organization},
  journal = {Nature},
  volume = {363},
  year = {1993},
  pages = {429},
}

@article{OhkiReid2005,
  author = {Ohki, Kenichi and Chung, Sooyang and Ch'ng, Yeang H. and Kara, Prakash and Reid, R. Clay},
  title = {Functional imaging with cellular resolution reveals precise micro-architecture in visual cortex},
  journal = {Nature},
  volume = {433},
  year = {2005},
  pages = {597--603},
}

@article{LienScanziani2013,
  author = {Lien, A. David and Scanziani, Massimo},
  title = {Cortical direction selectivity emerges from recurrence in a thalamocortical circuit},
  journal = {Nature},
  volume = {493},
  year = {2013},
  pages = {252--256},
}

@book{Hubel1988,
  author = {Hubel, David H.},
  title = {Eye, Brain, and Vision},
  publisher = {W. H. Freeman},
  address = {New York},
  year = {1988},
  series = {Scientific American Library},
}

@book{DayanAbbott2001,
  author = {Dayan, Peter and Abbott, L. F.},
  title = {Theoretical Neuroscience: Computational and Mathematical Modeling of Neural Systems},
  publisher = {MIT Press},
  address = {Cambridge, MA},
  year = {2001},
}

@book{KandelSchwartzJessell2013,
  author = {Kandel, Eric R. and Schwartz, James H. and Jessell, Thomas M. and Siegelbaum, Steven A. and Hudspeth, A. J.},
  title = {Principles of Neural Science},
  edition = {5th},
  publisher = {McGraw-Hill},
  address = {New York},
  year = {2013},
}