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_neurite_volume_density_failed_convex_hull():
with patch('neurom.features.neuritefunc.convex_hull',
side_effect=scipy.spatial.qhull.QhullError('boom')):
vol_density = _nf.neurite_volume_density(NRN)
nt.ok_(vol_density, np.nan)