3D neuron models extracted from EM image segmentation

- skeletonization. extract the neuron skeletons based on image segmentation (colored labels).
- morphological analysis with a lot of features.
- neural networks including synapses. Most of current morphological analysis tools do not have synapses.
- cell type classification based on NBLAST. use NBLAST to compute similarity scores between neurons.

run this inside julia package REPL (hit `]`

to enter package mode):

```
add RealNeuralNetworks
```

`julia skeletonize.jl -h`

```
sudo docker build . -t realneuralnetworks
```

The secret files containing cloud authentication keys follows the cloudvolume secret format.

```
docker run -v /tmp:/tmp -v $HOME/.cloudvolume/secrets:/root/.cloudvolume/secrets --net=host realneuralnetworks julia skeletonize.jl -h
```

```
using RealNeuralNetworks.NodeNets
using RealNeuralNetworks.Neurons
using RealNeuralNetworks.SWCs
nodeNet = NodeNet(seg::Array{UInt32,3}; obj_id = convert(UInt32,77605))
neuron = Neuron( nodeNet )
swc = SWC(neuron)
SWCs.save(swc, tempname()*".swc")
```

- [x] arbor density
- [x] total path length
- [x] number of segment points
- [x] Median segment length is the median dendritic segment length of all the segments starting and ending at irreducible nodes (in μm). Irreducible nodes are the points of the dendritic arbor corresponding to soma, branching points or terminal points.
- [x] 3D sholl analysis.
- [ ] Hull area is the area of the tightest convex hull containing the z-projection of the dendritic arbor (in μm2).
- [ ] volume of the convex hull around all neurites
- [x] Average angle is the mean of the positive angle between (parent node, node) and (node, child node) vectors, where node, parent node and child node are irreducible nodes (in radians).
- [x] Average tortuosity is the average value of the ratio of the actual dendritic length to the Euclidean distance between irreducible nodes.
- [x] Asymmetry is the distance of the soma node to the dendritic arbor (skeleton) centre of mass (in nm).
- [x] Typical radius (λ) is the root-mean-square distance of dendritic arbor points to the centre of mass (in nm).
- [x] fractal dimension.
- [x] longest neurite length.
- [ ] Strahler number
- [ ] distribution of Euclidian distance of segmentes from soma (third principal component)
- [ ] distribution of Euclidian distance of segmentes from soma as a function of segment order (third principal component)
- [ ] number of segmentes per segment order (second principal component)
- [ ] distribution of morphological distance of segmentes from soma along the skeleton as a function of segment order (first principal component)

- [x] ratio of tail diameter to head. could be useful to identify spines.
- [x] segment order
- [x] segment length
- [x] branching angle. computation using dot product
- [x] tortuosity / curvature. caa4486b501f936743f781782e6561833da7e413
- [x] distance to root path length
- [ ] segment asymmetry
- [x] average radius. easy to compute with radius list.

The skeletonization was originally implemented in Matlab by Alexander Bae using TEASAR algorithm, which was translated to Julia by Nicholas Turner.

09/18/2017

6 months ago

266 commits