def test_neurite_volumes(): # note: cannot use SIMPLE since it lies in a plane total_volumes = _nf.total_volume_per_neurite(NRN) assert_allclose(total_volumes, [ 271.94122143951864, 281.24754646913954, 274.98039928781355, 276.73860261723024 ])
def test_total_volume_per_neurite(): vol = _nf.total_volume_per_neurite(NRN) nt.eq_(len(vol), 4) # calculate the volumes by hand and compare vol2 = [sum(_sf.section_volume(s) for s in n.iter_sections()) for n in NRN.neurites] nt.eq_(vol, vol2) # regression test ref_vol = [271.94122143951864, 281.24754646913954, 274.98039928781355, 276.73860261723024] assert_allclose(vol, ref_vol)
def test_neurite_volume_density(): vol = np.array(_nf.total_volume_per_neurite(NRN)) hull_vol = np.array([convex_hull(n).volume for n in nm.iter_neurites(NRN)]) vol_density = _nf.neurite_volume_density(NRN) nt.eq_(len(vol_density), 4) nt.ok_(np.allclose(vol_density, vol / hull_vol)) ref_density = [ 0.43756606998299519, 0.52464681266899216, 0.24068543213643726, 0.26289304906104355 ] assert_allclose(vol_density, ref_density)
def test_total_volume_per_neurite(): vol = _nf.total_volume_per_neurite(NRN) assert len(vol) == 4 # calculate the volumes by hand and compare vol2 = [ sum(sectionfunc.section_volume(s) for s in n.iter_sections()) for n in NRN.neurites ] assert vol == vol2 # regression test ref_vol = [ 271.94122143951864, 281.24754646913954, 274.98039928781355, 276.73860261723024 ] assert np.allclose(vol, ref_vol)