Insert point processes in a Starburst Cell model
Posted: Wed Sep 03, 2014 1:16 pm
I'm trying to model a Starburst Amacrine Cell receiving inputs from bipolar cells. I've done this before in a simplified model of a neuron with 4 perpendicular dendrites. I create a circular grid of bipolar inputs and then, based on the distance and angle, I determine over which dendrites the inputs falls. However, now I'm trying to use a morphological detailed model like this http://neuromorpho.org/neuroMorpho/Keyw ... =starburst.
The model consist on a detailed morphology coded by a long list of pt3dadd calls.
Does anybody have an idea about how to connect the bipolar inputs to the corresponding dendrite?
For clarity, I'll provide the code I use in the simplified model.
After this, I just insert the inputs on dendrite(raux).
However, this approach doesn't work with this more sophisticated model.
I've thought about getting the pt3 coordinates somehow, but the documentation hasn't been helpful.
Any ideas will be really appreciated.
Thanks!
The model consist on a detailed morphology coded by a long list of pt3dadd calls.
Does anybody have an idea about how to connect the bipolar inputs to the corresponding dendrite?
For clarity, I'll provide the code I use in the simplified model.
Code: Select all
//Initial values for the bipolar inputs
xmin = -ldend
xmax = ldend
ymin = -ldend
ymax = ldend
y = ymin
fila = 0
n_is = 0
//Synaptic inputs location
while (y<=ymax) {
x = xmin + (1+(-1)^fila)*d_is/4
while (x<=xmax) {
a = r1.repick()
b = r2.repick()
if (sqrt((x+a*cos(b))^2+(y+a*sin(b))^2)<ldend) { //Circular array
xaux.append(x+a*cos(b))
yaux.append(y+a*sin(b))
n_is += 1
}
x += d_is
}
y += d_is*sin(PI/3)
fila += 1
}
//Synaptic input to dendrite connection
objref numdend, raux
numdend = new Vector(n_is) //Dendrite number
raux = new Vector(n_is) //Distance from the soma
for k = 0, n_is-1 {
angle = atan2(yaux.x[k], xaux.x[k]) //Determines the angle
if (angle < 0) {
angle += 2*PI
}
numdend.x[k] = int(angle*ndend/(2*PI)+0.5) //Determines the corresponding dendrite
if (numdend.x[k]==ndend) {
numdend.x[k]=0
}
raux.x[k]=sqrt(xaux.x[k]^2+yaux.x[k]^2)/ldend //Normalised distance to the soma
}
However, this approach doesn't work with this more sophisticated model.
I've thought about getting the pt3 coordinates somehow, but the documentation hasn't been helpful.
Any ideas will be really appreciated.
Thanks!