(c) 2020 by Sophie Laturnus, Adam von Daranyi, Ziwei Huang and Philipp Berens.
MorphoPy is a Python3 package that uses networkX to compute and show information about neurites. The input can be passed with single swc-files or it can handle whole directories with multiple files at once. You can use MorphoPy imported in Python or from command line as batch-tool as well.
The current working build:
version 0.6
All builds are tested on Linux (Debian and CentOS) and Windows 10.
You can find all working builds at pypi.org.
In the following, all commands written in boxes need to be entered into your terminal.
- homebrew (to install latest version of python3)
- python >3.4: without homebrew you can find python here
- pip
1) If you want the latest python3 version you need to install homebrew, for that just run this command in your terminal:
ruby -e "$(curl -fsSL https://raw.githubusercontent.com/Homebrew/install/master/install)"
2) Then you can easily install the newest version of python3:
brew install python
3) Now you should have python3 and pip3 installed. You can check the version like this:
python3 --version
pip3 -V
If both are enough you can continue with installing MorphoPy, see below.
- python >3.4
- pip 3
Install Python3 on Debian/Ubuntu systems:
apt-get install python3 python3-pip
Install Python3 on RedHat/CentOS/Fedora systems (perhaps python version differs):
yum install python36 python3-pip
That's it. All requirements are met and you can continue with the MorphoPy installation. See below.
- python >3.5 (x64): download from here
- pip : be sure you selected it during installation with the base python package
- Microsoft Build Tools >14.0 you can download them here
1) Install python with pip by executing the downloaded installation file and be sure to check the option to add python paths to enviroment variable at the first step!
2) Install Microsoft Build Tools with the automatic installation tool from the website above. Select C/C++ Compiler Tools at the selection page
All requirements are met now, see below for continue with MorphoPy installation.
Install the python package with pip3 and get the latest build:
pip3 install morphopy
You can simply uninstall the package with pip3:
pip3 uninstall morphopy
After the installation you can simply call:
MorphoPy.py <options>
Help:
Usage: MorphoPy.py -c <compute_feature> [--wide | --func <persistence_function> | --conf <config_file>]
[-f <swc_file> | -d <directory>] [-o <output directory>]
Options:
-c, --compute parameter for selecting the computing feature:
persistence: compute persistence data
stats : compute morphometric statistics
density : create density maps
statistics options:
--long you can change your output format, in long
format you get all values in a separate row.
(default: all values in one row)
persistence options:
--func if persistence is selected as feature, you can
specify with this option a method function you
want to use at computing the persistence.
(default: radial distance function)
density map options:
--conf if density map is selected, you can pass a
config file with more parameters for creating
the density maps. (optional)
-f, --file specifies a swc-file as input for Morphopy,
if no file is selected, directory is used
-d, --directory specifies a directory as input for swc-files.
(default: working directory)
-o, --output specifies the output directory for saving the
results in. (default: same as source)
Available functions for persistence at the moment are:
- radial_distance (default function)
- height
- path_length
- branch_order
A sample config file for density maps looks like this (stored in a text file):
[global]
# specific distance for resampling nodes:
distance: 1
# width of each bin in microns across all dimensions
#bin_size: 20
# number of bins for each dimension (only if you don't use bin_size)
n_bins_x: 20
n_bins_y: 20
n_bins_z: 20
# if true: probabilty density is returned, count histogram otherwise
density: True
# smoothing the density data
smooth: True
# sigma used at smoothing
sigma: 2
# normalization bounds for density map:
[norm_bound]
r_max_x: 238.85
r_max_y: 140.95
r_max_z: 285.97
r_min_x: -236.17
r_min_y: -24.2
r_min_z: -173.72
We tested MorphoPy to the best of our knowledge and abilities in the scope of several projects. If you still find a bug or you are missing a feature, please do not hesitate to contact us via GitHub issues. Please try to provide a minimal example that reproduces the bug you want to be fixed. If you want to develop the code base further, you can work with git pull requests. Please make sure that you document the code and add tests and examples of how to use your code.
If you are using MorphoPy for your research or your work project please make sure to cite us and this repository:
@misc{morphopy,
author = {Laturnus, Sophie and von Daranyi, Adam and Huang, Ziwei and Berens, Philipp},
title = {MorphoPy},
publisher = {GitHub},
journal = {GitHub repository},
url = {https://github.com/berenslab/MorphoPy}
}
Important: MorphoPy requires the soma to be one single point. If more than 3 somatic points are present in the reconstruction file they will be automatically collapsed to the centroid of their convex hull on file loading. If the soma is described by 2 to 3 points they will be automatically collapsed to their mean (also see utils.get_standardized_swc).
MorphoPy currently only supports neurites that connect back to the soma. This means, axons that emerge from dendritic structures can not be handled.
A neuron is represented as a directed acyclic graph with node attributes id, x-, y-, z- position, radius and type_id (soma: 1, axon: 2, dendrite: 3, apical dendrite: 4), and with edge attributes path_length and euclidean_dist. Positions, radius and length mesaures are assumed to be given in microns.
Fig. 1: Node and edge attributes associated with each neuron graph.
All data is stored in the tidy data format.
Please also refer to our tutorial.
Density maps are marginal histograms over the neural mass. MorphoPy allows you to create density maps of different projections through the function compute_denisty_maps(). Per default it computes x, y, z, xy, xz and yz density maps from the point cloud of the original reconstruction. The point cloud is constructed through resampling along all neurites with a default distance of 1 micron. The resulting point cloud is then binned into bins of 20 microns and smoothed using Gaussian smoothing with std of 1.
However, you can customize all these parameters by passing a config file to the function (see above).
MorphoPy offers a default selection of 28 single-valued morphometric statistics, namely:
- number of branch points
- width (x-extend), depth (y-extend), height (z-extend)
- number of tips
- number of neurites extending from the soma directly (stems)
- the total path length (in microns)
- average and maximal radius thickness (with the soma excluded)
- total surface and volume
- maximal path distance to the soma
- maximal branch order
- maximal, min and median path angle
- average soma exit angle
- maximal path length of a segment
- median intermediate and median terminal segment length
- log of max, min and median tortuosity across all edges (= path length/euclidean length)
- max, min and average branch angle
- maximal branching degree (with soma excluded)
- weighted proportional sum of absolute deviations as a measure of tree asymmetry (for more details see https://www.sciencedirect.com/science/article/pii/0165027086901196)
Fig. 2: Explanatory schematic of the morphometric statistics that can be computed on all nodes. Left: distance measures, Right: angles.
Morphometric distributions are not (yet) available via the command line tool.
Frequency histograms or Gaussian kernel density estimates can be queried for all listed key statistics using the
methods get_histogram(key) or get_kde_distribution(key). If you provide a distance measure (e.g. branch order,
path distance from soma or radial distance) the returned distribution will be two-dimensional and allows to investigate
a features' development across space.
Additionally, it is possible to compute Sholl intersection profiles using the function get_sholl_intersection_profile().
Key statistics are
- branch orders
- Strahler order
- branch angles
- path angles
- root angles
- thickness
- segment lengths
- path length to soma
- radial distance
Persistence diagrams are a concept from topology. They have been introduced as descriptors of neural morphologies by Kanari et al. and Li et al.
Fig. 3: Schematic of how a persistence diagram is generated. The longest branch with the longest 'lifetime' is marked in red. Taken from Kanari et al. 2018.
The recorded birth and death times in the figure above are based on a certain distance (or lifetime) function. MorphoPy implements four different distance functions to choose from: radial distance (default), height, path length and branch order.
They all compute the distance of a point with respect to the soma. In the command line tool you can switch between them using the --func keyword (see above).
To provide your own distance function, add its code and its keyword to the persistence_functions.py file, but make sure
that the distance functions interface fits the specification custom_distance(networkx.DiGraph,node_id_end, node_id_start)
(see Fig. 4 and the tutorial for an example).
Fig. 4: How to add a custom persistence distance function. To be able to call it from the command line you need to add it to the functions list.
If you are using the API you can simply pass a function to the get_persistence()-method (see the tutorial for an example).
You want to compute your own features? Go for it! We recommend you to check out networkx and shapely for more options.