Spatial networks: 4x3 grid with pyramidal cells and interneuronsΒΆ


Run this example as a Jupyter notebook:

See our guide for more information and troubleshooting.


Create a 4x3 grid with one pyramidal cell and one interneuron at each position.

BCCN Tutorial @ CNS*09 Hans Ekkehard Plesser, UMB

import matplotlib.pyplot as plt
import nest
import numpy as np

nest.ResetKernel()

nest.CopyModel("iaf_psc_alpha", "pyr")
nest.CopyModel("iaf_psc_alpha", "in")

pos = nest.spatial.grid(shape=[4, 3], extent=[2.0, 1.5])

ctx_pyr = nest.Create("pyr", positions=pos)
ctx_in = nest.Create("in", positions=pos)

nest.PrintNodes()

# extract position information
ppyr = nest.GetPosition(ctx_pyr)
pin = nest.GetPosition(ctx_in)

ppyr_x = np.array([x for x, y in ppyr])
ppyr_y = np.array([y for x, y in ppyr])

pin_x = np.array([x for x, y in pin])
pin_y = np.array([y for x, y in pin])

# plot
plt.clf()
plt.plot(pin_x - 0.05, ppyr_y - 0.05, "bo", markersize=20, label="Pyramidal", zorder=2)
plt.plot(pin_x + 0.05, pin_y + 0.05, "ro", markersize=20, label="Interneuron", zorder=2)
plt.plot(
    pin_x, ppyr_y, "o", markerfacecolor=(0.7, 0.7, 0.7), markersize=60, markeredgewidth=0, zorder=1, label="_nolegend_"
)

# beautify
plt.axis([-1.0, 1.0, -1.0, 1.0])
plt.axes().set_aspect("equal", "box")
plt.axes().set_xticks((-0.75, -0.25, 0.25, 0.75))
plt.axes().set_yticks((-0.5, 0, 0.5))
plt.grid(True)
plt.xlabel("4 Columns, Extent: 1.5")
plt.ylabel("3 Rows, Extent: 1.0")
plt.legend(numpoints=1)

plt.show()

# plt.savefig('ctx_2n.png')

Gallery generated by Sphinx-Gallery