A collection of CLIs and python function related to morphology handling
Project description
A toolbox for morphology editing. It aims to provide small helper programs that perform simple tasks.
Currently MorphTool provides:
A morphology diffing tool (via CLI or python)
A file converter: to convert morphology files from/to the following formats: SWC, ASC, H5
A soma area calculator (as computed by NEURON, requires the NEURON python module)
Installation
In a fresh virtualenv:
pip install morph-tool
Usage
In a shell, do:
morph-tool --help
Currently the three sub-command are:
morph-tool convert input_file output_file
morph-tool diff morph1 morph2
morph-tool soma-surface input_file
Morphology diffing
One can compare two morphologies with the CLI
morph-tool diff morph1 morph2
whose error code is 0 if morphologies are the same, else 1.
Morphologies with different formats can be compared.
Morphologies are considered different if one or more of the following properties differ:
number of root sections
sections type
sections point array
sections diameter array
sections perimeter array
sections number of children
The soma are not taken into consideration
The same functionality is also avalaible through the python API:
from morph_tool import diff
# The result can be used as a boolean:
if diff(filename1, filename2):
print('morphologies differ')
# And also contains information about how morphologies differ
result = diff(filename1, filename2)
print(result.info)
Converter
What is supported ?
The converter can be used to write morphology in different formats. Currently, supported formats are ASC, SWC and H5.
As each format has its own specificities, data specific to a given format will be discarded. This means the following will be lost during conversion:
spines (present in the ASC format)
all H5 metadata
the perimeter and mitochondrial data of the H5 format
Soma intricacies
Multiple formats are being used to represent somas (mainly) depending on the file format. For more information about file format, see the neuromorpho.org specification or MorphIO specification
Because different soma format represent soma in different planes, soma format conversion is not a bijective transformation. For example, it is not possible to have an accurate conversion from a soma contour in the XY plane from a H5 file to a SWC soma which is represented as a cylinder along Y.
As a result, it has been chosen to take the soma surface as an invariant. This means soma surfaces of the input and output morphologies, as computed by NEURON, should be preserved.
Here are the possible cases for the soma conversion:
SWC input file:
SWC output file -> no conversion
H5 or ASC output file:
Depending on the original soma type:
Soma stack of cylinders: The soma is converted to a contour in the XY plane. The points of the new contour are the outline of the soma stack projected in the XY plane.
Soma three point cylinder: The soma becomes a sphere of same surface. The contour made by the circle of biggest section in the XY plane is sampled in 20 points written to disk.
Soma sphere (soma represented by a single point representing the center of a sphere and its radius): the contour made by the circle of biggest section in the XY plane is sampled in 20 points written to disk.
other: Not in SWC spec -> not supported
H5 or ASC input file:
H5 output file -> no conversion needed
ASC output file.
Depending on soma type:
Soma single point sphere (soma represented by a single point representing the center of a sphere and its radius): the contour made by the circle of biggest section in the XY plane is sampled in 20 points written to disk.
Soma contour: no conversion needed
other: not in H5/ASC specs -> not supported
SWC:
Depending on soma format:
Soma single point sphere: no conversion needed
Soma contour: A soma stack of cylinder is generated. Each cylinder of the stack has its center and its axis along the principal direction of the contour. The radius of each stack is chosen such that it minimises the distance between the cylinder and the contour.
other: not in H5/ASC specs -> not supported
Example:
from morph_tool import convert
convert(inputfile, outputfile)
# Additionally the morphology can be recentered or written according to the NEURON neurite order during the conversion
convert(inputfile, outputfile, recenter=True, nrn_order=True)
The same for bash
morph-tool convert file ./inputfile ./outputfile
# with additional options
morph-tool convert file --recenter --nrn-order ./inputfile ./outputfile
# or an entire folder
morph-tool convert folder -ext SWC ./h5_input_folder ./swc_output_folder
# for more info use
morph-tool convert folder --help
NRN simulator compartment coordinates
The NRN simulator splits each section into chunks of equal length (equal only among a given section). These compartments do not really exist in the physical world but we can remap them to paths along the section. Each compartment can be associated to a path (a list of 3D points) such that the path and the compartment have the same path-length.
The following function can be used to access the mapping NeuroM section ID -> list of paths for the section:
morph_tool.nrnhines.NeuroM_section_to_NRN_compartment_paths
Example (in 2D) for one section:
(1, 2) ------ (2, 2)
|
|
|
|
|
|
|
|
|
(0, 0) ------- (1, 0)
Splitting this section into 3 compartments would results in the following paths:
[[0. , 0. ],
[1. , 0. ],
[1. , 0.33333333]]
[[1. , 0.33333333],
[1. , 1.66666667]]
[[1. , 1.66666667],
[1. , 2. ],
[2. , 2. ]]
Dendrogram with synapses
This functionality is available only when the package is installed with dendrogram extras:
pip install morph-tool
Draw NeuroM dendrogram with synapses on it. Synapses must be represented as a DataFrame. Required columns in this dataframe are:
from morph_tool import dendrogram
required_columns = [dendrogram.SOURCE_NODE_ID, dendrogram.TARGET_NODE_ID,
dendrogram.POST_SECTION_ID, dendrogram.POST_SECTION_POS,
dendrogram.PRE_SECTION_ID, dendrogram.PRE_SECTION_POS]
or equivalently
required_columns = ['@source_node', '@target_node',
'afferent_section_id', 'afferent_section_pos',
'efferent_section_id', 'efferent_section_pos']
For usage examples look at examples/dendrogram.py.
Contributing
If you want to improve the project or you see any issue, every contribution is welcome. Please check the contribution guidelines for more information.
Acknowledgements
This research was supported by the EBRAINS research infrastructure, funded from the European Union’s Horizon 2020 Framework Programme for Research and Innovation under the Specific Grant Agreement No. 945539 (Human Brain Project SGA3).
License
morph-tool is licensed under the terms of the GNU Lesser General Public License version 3. Refer to COPYING.LESSER and COPYING for details.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
File details
Details for the file morph-tool-2.3.0.tar.gz
.
File metadata
- Download URL: morph-tool-2.3.0.tar.gz
- Upload date:
- Size: 1.6 MB
- Tags: Source
- Uploaded using Trusted Publishing? No
- Uploaded via: twine/3.3.0 pkginfo/1.7.0 requests/2.25.1 setuptools/52.0.0 requests-toolbelt/0.9.1 tqdm/4.56.0 CPython/3.8.7
File hashes
Algorithm | Hash digest | |
---|---|---|
SHA256 | 08583995632633c67e585d8d465c491a3dc28fd77950d6735c94c460d5dfd46f |
|
MD5 | 97437cd800305eb2ced1013c2932f991 |
|
BLAKE2b-256 | a847ff57817cc1e72c58c4a8316e12bd0e76335c9f249b62df85a2eaf9d3d639 |