In [1]:
import pickle
import btmorph
In [2]:
d = pickle.load(open("res_data1L_21051982.pkl","r"))
In [3]:
res=[] # peak of the transfer impedance at each location
amp=[] # Q-factor at each location
for key in d.keys():
amp.append(d[key][0])
res.append(d[key][1])
for i in range(2): # some misnumbering in my code with 3 point soma
amp.append(amp[-1])
res.append(res[-1])
In [4]:
# this is figure 4A in the manuscript
btmorph.plot_2D_SWC("P20-DEV139.CNG.swc",cs=res,new_fig=False)
In [5]:
# this image shows the amplitude of the transfer impedance
btmorph.plot_2D_SWC("P20-DEV139.CNG.swc",cs=amp,new_fig=False)
# the color bar does not work well here because of the small range of amplitudes...
In [ ]: