Name | LFPykit JSON |
Version |
0.5.1
JSON |
| download |
home_page | https://github.com/LFPy/LFPykit |
Summary | Electrostatic models for multicompartment neuron models |
upload_time | 2023-04-12 07:10:54 |
maintainer | |
docs_url | None |
author | LFPy-team |
requires_python | >=3.7 |
license | |
keywords |
|
VCS |
|
bugtrack_url |
|
requirements |
No requirements were recorded.
|
Travis-CI |
No Travis.
|
coveralls test coverage |
No coveralls.
|
# LFPykit
This Python module contain freestanding implementations of electrostatic
forward models incorporated in LFPy
(https://github.com/LFPy/LFPy, https://LFPy.readthedocs.io).
The aim of the `LFPykit` module is to provide electrostatic models
in a manner that facilitates forward-model predictions of extracellular
potentials and related measures from multicompartment neuron models, but
without explicit dependencies on neural simulation software such as
NEURON (https://neuron.yale.edu, https://github.com/neuronsimulator/nrn),
Arbor (https://arbor.readthedocs.io, https://github.com/arbor-sim/arbor),
or even LFPy.
The `LFPykit` module can then be more easily incorporated with
these simulators, or in various projects that utilize them such as
LFPy (https://LFPy.rtfd.io, https://github.com/LFPy/LFPy).
BMTK (https://alleninstitute.github.io/bmtk/, https://github.com/AllenInstitute/bmtk),
etc.
Its main functionality is providing class methods that return two-dimensional
linear transformation matrices **M**
between transmembrane currents
**I** of multicompartment neuron models and some
measurement **Y** given by **Y**=**MI**.
The presently incorporated volume conductor models have been incorporated in
LFPy (https://LFPy.rtfd.io, https://github.com/LFPy/LFPy),
as described in various papers and books:
1. Linden H, Hagen E, Leski S, Norheim ES, Pettersen KH, Einevoll GT
(2014) LFPy: a tool for biophysical simulation of extracellular
potentials generated by detailed model neurons. Front.
Neuroinform. 7:41. doi: 10.3389/fninf.2013.00041
2. Hagen E, Næss S, Ness TV and Einevoll GT (2018) Multimodal Modeling of Neural
Network Activity: Computing LFP, ECoG, EEG, and MEG
Signals With LFPy 2.0. Front. Neuroinform. 12:92.
doi: 10.3389/fninf.2018.00092
3. Ness, T. V., Chintaluri, C., Potworowski, J., Leski, S., Glabska,
H., Wójcik, D. K., et al. (2015). Modelling and analysis of electrical
potentials recorded in microelectrode arrays (MEAs).
Neuroinformatics 13:403–426. doi: 10.1007/s12021-015-9265-6
4. Nunez and Srinivasan, Oxford University Press, 2006
5. Næss S, Chintaluri C, Ness TV, Dale AM, Einevoll GT and Wójcik DK
(2017). Corrected Four-sphere Head Model for EEG Signals. Front. Hum.
Neurosci. 11:490. doi: 10.3389/fnhum.2017.00490
## Build Status
[![DOI](https://zenodo.org/badge/288660131.svg)](https://zenodo.org/badge/latestdoi/288660131)
[![Coverage Status](https://coveralls.io/repos/github/LFPy/LFPykit/badge.svg?branch=master)](https://coveralls.io/github/LFPy/LFPykit?branch=master)
[![Documentation Status](https://readthedocs.org/projects/lfpykit/badge/?version=latest)](https://lfpykit.readthedocs.io/en/latest/?badge=latest)
[![flake8 lint](https://github.com/LFPy/LFPykit/actions/workflows/flake8.yml/badge.svg)](https://github.com/LFPy/LFPykit/actions/workflows/flake8.yml)
[![Python application](https://github.com/LFPy/LFPykit/workflows/Python%20application/badge.svg)](https://github.com/LFPy/LFPykit/actions?query=workflow%3A%22Python+application%22)
[![Upload Python Package](https://github.com/LFPy/LFPykit/workflows/Upload%20Python%20Package/badge.svg)](https://pypi.org/project/LFPykit)
[![Conda Recipe](https://img.shields.io/badge/recipe-lfpykit-green.svg)](https://anaconda.org/conda-forge/lfpykit)
[![Conda Downloads](https://img.shields.io/conda/dn/conda-forge/lfpykit.svg)](https://anaconda.org/conda-forge/lfpykit)
[![Conda Version](https://img.shields.io/conda/vn/conda-forge/lfpykit.svg)](https://anaconda.org/conda-forge/lfpykit)
[![Conda Platforms](https://img.shields.io/conda/pn/conda-forge/lfpykit.svg)](https://anaconda.org/conda-forge/lfpykit)
[![Binder](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/LFPy/LFPykit.git/master)
[![License](http://img.shields.io/:license-GPLv3+-green.svg)](http://www.gnu.org/licenses/gpl-3.0.html)
## Features
`LFPykit` presently incorporates different electrostatic forward models for extracellular potentials
and magnetic signals that has been derived using volume conductor theory.
In volume-conductor theory the extracellular potentials can be calculated from a distance-weighted sum of contributions from transmembrane currents of neurons.
Given the same transmembrane currents, the contributions to the magnetic field recorded both inside and outside the brain can also be computed.
The module presently incorporates different classes.
To represent the geometry of a multicompartment neuron model we have:
* `CellGeometry`:
Base class representing a multicompartment neuron geometry in terms
of segment x-, y-, z-coordinates and diameter.
Different classes built to map transmembrane currents of `CellGeometry` like instances
to different measurement modalities:
* `LinearModel`:
Base class representing a generic forward model
for subclassing
* `CurrentDipoleMoment`:
Class for predicting current dipole moments
* `PointSourcePotential`:
Class for predicting extracellular potentials
assuming point sources and point contacts
* `LineSourcePotential`:
Class for predicting extracellular potentials assuming
line sourcers and point contacts
* `RecExtElectrode`:
Class for simulations of extracellular potentials
* `RecMEAElectrode`:
Class for simulations of in vitro (slice) extracellular
potentials
* `OneSphereVolumeConductor`:
For computing extracellular potentials within
sand outside a homogeneous sphere
* `LaminarCurrentSourceDensity`:
For computing the 'ground truth' current source density across
cylindrical volumes aligned with the z-axis
* `VolumetricCurrentSourceDensity`:
For computing the 'ground truth' current source density on regularly
spaced 3D grid
Different classes built to map current dipole moments (i.e., computed using `CurrentDipoleMoment`)
to extracellular measurements:
* `eegmegcalc.FourSphereVolumeConductor`:
For computing extracellular potentials in
4-sphere head model (brain, CSF, skull, scalp)
from current dipole moment
* `eegmegcalc.InfiniteVolumeConductor`:
To compute extracellular potentials in infinite volume conductor
from current dipole moment
* `eegmegcalc.InfiniteHomogeneousVolCondMEG`:
Class for computing magnetic field from current dipole moments under the assumption
of infinite homogeneous volume conductor model
* `eegmegcalc.SphericallySymmetricVolCondMEG`:
Class for computing magnetic field from current dipole moments under the assumption
of a spherically symmetric volume conductor model
* `eegmegcalc.NYHeadModel`:
Class for computing extracellular potentials in detailed head volume
conductor model (https://www.parralab.org/nyhead)
Each class (except `CellGeometry`) should have a public method `get_transformation_matrix()`
that returns the linear map between the transmembrane currents or current dipole moment
and corresponding measurements (see Usage below)
## Usage
A basic usage example using a mock 3-segment stick-like neuron,
treating each segment as a point source in a linear, isotropic and homogeneous volume conductor,
computing the extracellular potential in ten different locations
alongside the cell geometry:
>>> # imports
>>> import numpy as np
>>> from lfpykit import CellGeometry, PointSourcePotential
>>> n_seg = 3
>>> # instantiate class `CellGeometry`:
>>> cell = CellGeometry(x=np.array([[0.] * 2] * n_seg), # (µm)
y=np.array([[0.] * 2] * n_seg), # (µm)
z=np.array([[10. * x, 10. * (x + 1)]
for x in range(n_seg)]), # (µm)
d=np.array([1.] * n_seg)) # (µm)
>>> # instantiate class `PointSourcePotential`:
>>> psp = PointSourcePotential(cell,
x=np.ones(10) * 10,
y=np.zeros(10),
z=np.arange(10) * 10,
sigma=0.3)
>>> # get linear response matrix mapping currents to measurements
>>> M = psp.get_transformation_matrix()
>>> # transmembrane currents (nA):
>>> imem = np.array([[-1., 1.],
[0., 0.],
[1., -1.]])
>>> # compute extracellular potentials (mV)
>>> V_ex = M @ imem
>>> V_ex
array([[-0.01387397, 0.01387397],
[-0.00901154, 0.00901154],
[ 0.00901154, -0.00901154],
[ 0.01387397, -0.01387397],
[ 0.00742668, -0.00742668],
[ 0.00409718, -0.00409718],
[ 0.00254212, -0.00254212],
[ 0.00172082, -0.00172082],
[ 0.00123933, -0.00123933],
[ 0.00093413, -0.00093413]])
A basic usage example using a mock 3-segment stick-like neuron,
treating each segment as a point source,
computing the current dipole moment and computing the potential in ten different
remote locations away from the cell geometry:
>>> # imports
>>> import numpy as np
>>> from lfpykit import CellGeometry, CurrentDipoleMoment, \
>>> eegmegcalc
>>> n_seg = 3
>>> # instantiate class `CellGeometry`:
>>> cell = CellGeometry(x=np.array([[0.] * 2] * n_seg), # (µm)
y=np.array([[0.] * 2] * n_seg), # (µm)
z=np.array([[10. * x, 10. * (x + 1)]
for x in range(n_seg)]), # (µm)
d=np.array([1.] * n_seg)) # (µm)
>>> # instantiate class `CurrentDipoleMoment`:
>>> cdp = CurrentDipoleMoment(cell)
>>> M_I_to_P = cdp.get_transformation_matrix()
>>> # instantiate class `eegmegcalc.InfiniteVolumeConductor` and map dipole moment to
>>> # extracellular potential at measurement sites
>>> ivc = eegmegcalc.InfiniteVolumeConductor(sigma=0.3)
>>> # compute linear response matrix between dipole moment and
>>> # extracellular potential
>>> M_P_to_V = ivc.get_transformation_matrix(np.c_[np.ones(10) * 1000,
np.zeros(10),
np.arange(10) * 100])
>>> # transmembrane currents (nA):
>>> imem = np.array([[-1., 1.],
[0., 0.],
[1., -1.]])
>>> # compute extracellular potentials (mV)
>>> V_ex = M_P_to_V @ M_I_to_P @ imem
>>> V_ex
array([[ 0.00000000e+00, 0.00000000e+00],
[ 5.22657054e-07, -5.22657054e-07],
[ 1.00041193e-06, -1.00041193e-06],
[ 1.39855769e-06, -1.39855769e-06],
[ 1.69852477e-06, -1.69852477e-06],
[ 1.89803345e-06, -1.89803345e-06],
[ 2.00697409e-06, -2.00697409e-06],
[ 2.04182029e-06, -2.04182029e-06],
[ 2.02079888e-06, -2.02079888e-06],
[ 1.96075587e-06, -1.96075587e-06]])
## Physical units
Notes on physical units used in `LFPykit`:
- There are no explicit checks for physical units
- Transmembrane currents are assumed to be in units of (nA)
- Spatial information is assumed to be in units of (µm)
- Voltages are assumed to be in units of (mV)
- Extracellular conductivities are assumed to be in units of (S/m)
- current dipole moments are assumed to be in units of (nA µm)
- Magnetic fields are assumed to be in units of (nA/µm)
## Dimensionality
- Transmembrane currents are represented by arrays with shape `(n_seg, n_timesteps)`
where `n_seg` is the number of segments of the neuron model.
- Current dipole moments are represented by arrays with shape `(3, n_timesteps)`
- Response matrices **M** have shape `(n_points, input.shape[0])` where `n_points` is
for instance the number of extracellular recording sites and `input.shape[0]`
the first dimension of the input; that is, the number of segments in case of
transmembrane currents or 3 in case of current dipole moments.
- predicted signals (except magnetic fields using `eegmegcalc.InfiniteHomogeneousVolCondMEG` or
`eegmegcalc.SphericallySymmetricVolCondMEG`) have shape `(n_points, n_timesteps)`
## Documentation
The online Documentation of `LFPykit` can be found here:
https://lfpykit.readthedocs.io/en/latest
## Dependencies
`LFPykit` is implemented in Python and is written (and continuously tested) for `Python >= 3.7`.
The main `LFPykit` module depends on `numpy`, `scipy` and `MEAutility` (https://github.com/alejoe91/MEAutility, https://meautility.readthedocs.io/en/latest/).
Running all unit tests and example files may in addition require `py.test`, `matplotlib`,
`neuron` (https://www.neuron.yale.edu),
(`arbor` coming) and
`LFPy` (https://github.com/LFPy/LFPy, https://LFPy.readthedocs.io).
## Installation
### From development sources (https://github.com/LFPy/LFPykit)
Install the current development version on https://GitHub.com using `git` (https://git-scm.com):
$ git clone https://github.com/LFPy/LFPykit.git
$ cd LFPykit
$ python setup.py install # --user optional
or using `pip`:
$ pip install . # --user optional
For active development, link the repository location
$ python setup.py develop # --user optional
### Installation of stable releases on PyPI.org (https://www.pypi.org)
Installing from the Python Package Index (https://www.pypi.org/project/lfpykit):
$ pip install lfpykit # --user optional
To upgrade the installation using pip:
$ pip install --upgrade --no-deps lfpykit
### Installation of stable releases on conda-forge (https://conda-forge.org)
Installing `lfpykit` from the `conda-forge` channel can be achieved by adding `conda-forge` to your channels with:
$ conda config --add channels conda-forge
Once the `conda-forge` channel has been enabled, `lfpykit` can be installed with:
$ conda install lfpykit
It is possible to list all of the versions of `lfpykit` available on your platform with:
$ conda search lfpykit --channel conda-forge
Raw data
{
"_id": null,
"home_page": "https://github.com/LFPy/LFPykit",
"name": "LFPykit",
"maintainer": "",
"docs_url": null,
"requires_python": ">=3.7",
"maintainer_email": "",
"keywords": "",
"author": "LFPy-team",
"author_email": "lfpy@users.noreply.github.com",
"download_url": "https://files.pythonhosted.org/packages/b4/e0/b5935801a38839b177267f53cf6781941e36d8f832e06612773fdf89504e/LFPykit-0.5.1.tar.gz",
"platform": null,
"description": "# LFPykit\n\nThis Python module contain freestanding implementations of electrostatic\nforward models incorporated in LFPy\n(https://github.com/LFPy/LFPy, https://LFPy.readthedocs.io).\n\nThe aim of the `LFPykit` module is to provide electrostatic models\nin a manner that facilitates forward-model predictions of extracellular\npotentials and related measures from multicompartment neuron models, but\nwithout explicit dependencies on neural simulation software such as\nNEURON (https://neuron.yale.edu, https://github.com/neuronsimulator/nrn),\nArbor (https://arbor.readthedocs.io, https://github.com/arbor-sim/arbor),\nor even LFPy.\nThe `LFPykit` module can then be more easily incorporated with\nthese simulators, or in various projects that utilize them such as\nLFPy (https://LFPy.rtfd.io, https://github.com/LFPy/LFPy).\nBMTK (https://alleninstitute.github.io/bmtk/, https://github.com/AllenInstitute/bmtk),\netc.\n\nIts main functionality is providing class methods that return two-dimensional\nlinear transformation matrices **M**\nbetween transmembrane currents\n**I** of multicompartment neuron models and some\nmeasurement **Y** given by **Y**=**MI**.\n\nThe presently incorporated volume conductor models have been incorporated in\nLFPy (https://LFPy.rtfd.io, https://github.com/LFPy/LFPy),\nas described in various papers and books:\n\n1. Linden H, Hagen E, Leski S, Norheim ES, Pettersen KH, Einevoll GT\n (2014) LFPy: a tool for biophysical simulation of extracellular\n potentials generated by detailed model neurons. Front.\n Neuroinform. 7:41. doi: 10.3389/fninf.2013.00041\n\n2. Hagen E, N\u00e6ss S, Ness TV and Einevoll GT (2018) Multimodal Modeling of Neural\n Network Activity: Computing LFP, ECoG, EEG, and MEG\n Signals With LFPy 2.0. Front. Neuroinform. 12:92.\n doi: 10.3389/fninf.2018.00092\n\n3. Ness, T. V., Chintaluri, C., Potworowski, J., Leski, S., Glabska,\n H., W\u00f3jcik, D. K., et al. (2015). Modelling and analysis of electrical\n potentials recorded in microelectrode arrays (MEAs).\n Neuroinformatics 13:403\u2013426. doi: 10.1007/s12021-015-9265-6\n\n4. Nunez and Srinivasan, Oxford University Press, 2006\n\n5. N\u00e6ss S, Chintaluri C, Ness TV, Dale AM, Einevoll GT and W\u00f3jcik DK\n (2017). Corrected Four-sphere Head Model for EEG Signals. Front. Hum.\n Neurosci. 11:490. doi: 10.3389/fnhum.2017.00490\n\n\n## Build Status\n\n[![DOI](https://zenodo.org/badge/288660131.svg)](https://zenodo.org/badge/latestdoi/288660131)\n[![Coverage Status](https://coveralls.io/repos/github/LFPy/LFPykit/badge.svg?branch=master)](https://coveralls.io/github/LFPy/LFPykit?branch=master)\n[![Documentation Status](https://readthedocs.org/projects/lfpykit/badge/?version=latest)](https://lfpykit.readthedocs.io/en/latest/?badge=latest)\n[![flake8 lint](https://github.com/LFPy/LFPykit/actions/workflows/flake8.yml/badge.svg)](https://github.com/LFPy/LFPykit/actions/workflows/flake8.yml)\n[![Python application](https://github.com/LFPy/LFPykit/workflows/Python%20application/badge.svg)](https://github.com/LFPy/LFPykit/actions?query=workflow%3A%22Python+application%22)\n[![Upload Python Package](https://github.com/LFPy/LFPykit/workflows/Upload%20Python%20Package/badge.svg)](https://pypi.org/project/LFPykit)\n[![Conda Recipe](https://img.shields.io/badge/recipe-lfpykit-green.svg)](https://anaconda.org/conda-forge/lfpykit)\n[![Conda Downloads](https://img.shields.io/conda/dn/conda-forge/lfpykit.svg)](https://anaconda.org/conda-forge/lfpykit)\n[![Conda Version](https://img.shields.io/conda/vn/conda-forge/lfpykit.svg)](https://anaconda.org/conda-forge/lfpykit)\n[![Conda Platforms](https://img.shields.io/conda/pn/conda-forge/lfpykit.svg)](https://anaconda.org/conda-forge/lfpykit)\n[![Binder](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/LFPy/LFPykit.git/master)\n[![License](http://img.shields.io/:license-GPLv3+-green.svg)](http://www.gnu.org/licenses/gpl-3.0.html)\n\n\n## Features\n\n`LFPykit` presently incorporates different electrostatic forward models for extracellular potentials\nand magnetic signals that has been derived using volume conductor theory.\nIn volume-conductor theory the extracellular potentials can be calculated from a distance-weighted sum of contributions from transmembrane currents of neurons.\nGiven the same transmembrane currents, the contributions to the magnetic field recorded both inside and outside the brain can also be computed.\n\nThe module presently incorporates different classes.\nTo represent the geometry of a multicompartment neuron model we have:\n\n* `CellGeometry`:\n Base class representing a multicompartment neuron geometry in terms\n of segment x-, y-, z-coordinates and diameter.\n\nDifferent classes built to map transmembrane currents of `CellGeometry` like instances\nto different measurement modalities:\n\n* `LinearModel`:\n Base class representing a generic forward model\n for subclassing\n* `CurrentDipoleMoment`:\n Class for predicting current dipole moments\n* `PointSourcePotential`:\n Class for predicting extracellular potentials\n assuming point sources and point contacts\n* `LineSourcePotential`:\n Class for predicting extracellular potentials assuming\n line sourcers and point contacts\n* `RecExtElectrode`:\n Class for simulations of extracellular potentials\n* `RecMEAElectrode`:\n Class for simulations of in vitro (slice) extracellular\n potentials\n* `OneSphereVolumeConductor`:\n For computing extracellular potentials within\n sand outside a homogeneous sphere\n* `LaminarCurrentSourceDensity`:\n For computing the 'ground truth' current source density across\n cylindrical volumes aligned with the z-axis\n* `VolumetricCurrentSourceDensity`:\n For computing the 'ground truth' current source density on regularly\n spaced 3D grid\n\nDifferent classes built to map current dipole moments (i.e., computed using `CurrentDipoleMoment`)\nto extracellular measurements:\n\n* `eegmegcalc.FourSphereVolumeConductor`:\n For computing extracellular potentials in\n 4-sphere head model (brain, CSF, skull, scalp)\n from current dipole moment\n* `eegmegcalc.InfiniteVolumeConductor`:\n To compute extracellular potentials in infinite volume conductor\n from current dipole moment\n* `eegmegcalc.InfiniteHomogeneousVolCondMEG`:\n Class for computing magnetic field from current dipole moments under the assumption\n of infinite homogeneous volume conductor model\n* `eegmegcalc.SphericallySymmetricVolCondMEG`:\n Class for computing magnetic field from current dipole moments under the assumption\n of a spherically symmetric volume conductor model\n* `eegmegcalc.NYHeadModel`:\n Class for computing extracellular potentials in detailed head volume\n conductor model (https://www.parralab.org/nyhead)\n\nEach class (except `CellGeometry`) should have a public method `get_transformation_matrix()`\nthat returns the linear map between the transmembrane currents or current dipole moment\nand corresponding measurements (see Usage below)\n\n\n## Usage\n\nA basic usage example using a mock 3-segment stick-like neuron,\ntreating each segment as a point source in a linear, isotropic and homogeneous volume conductor,\ncomputing the extracellular potential in ten different locations\nalongside the cell geometry:\n\n >>> # imports\n >>> import numpy as np\n >>> from lfpykit import CellGeometry, PointSourcePotential\n >>> n_seg = 3\n >>> # instantiate class `CellGeometry`:\n >>> cell = CellGeometry(x=np.array([[0.] * 2] * n_seg), # (\u00b5m)\n y=np.array([[0.] * 2] * n_seg), # (\u00b5m)\n z=np.array([[10. * x, 10. * (x + 1)]\n for x in range(n_seg)]), # (\u00b5m)\n d=np.array([1.] * n_seg)) # (\u00b5m)\n >>> # instantiate class `PointSourcePotential`:\n >>> psp = PointSourcePotential(cell,\n x=np.ones(10) * 10,\n y=np.zeros(10),\n z=np.arange(10) * 10,\n sigma=0.3)\n >>> # get linear response matrix mapping currents to measurements\n >>> M = psp.get_transformation_matrix()\n >>> # transmembrane currents (nA):\n >>> imem = np.array([[-1., 1.],\n [0., 0.],\n [1., -1.]])\n >>> # compute extracellular potentials (mV)\n >>> V_ex = M @ imem\n >>> V_ex\n array([[-0.01387397, 0.01387397],\n [-0.00901154, 0.00901154],\n [ 0.00901154, -0.00901154],\n [ 0.01387397, -0.01387397],\n [ 0.00742668, -0.00742668],\n [ 0.00409718, -0.00409718],\n [ 0.00254212, -0.00254212],\n [ 0.00172082, -0.00172082],\n [ 0.00123933, -0.00123933],\n [ 0.00093413, -0.00093413]])\n\n\nA basic usage example using a mock 3-segment stick-like neuron,\ntreating each segment as a point source,\ncomputing the current dipole moment and computing the potential in ten different\nremote locations away from the cell geometry:\n\n >>> # imports\n >>> import numpy as np\n >>> from lfpykit import CellGeometry, CurrentDipoleMoment, \\\n >>> eegmegcalc\n >>> n_seg = 3\n >>> # instantiate class `CellGeometry`:\n >>> cell = CellGeometry(x=np.array([[0.] * 2] * n_seg), # (\u00b5m)\n y=np.array([[0.] * 2] * n_seg), # (\u00b5m)\n z=np.array([[10. * x, 10. * (x + 1)]\n for x in range(n_seg)]), # (\u00b5m)\n d=np.array([1.] * n_seg)) # (\u00b5m)\n >>> # instantiate class `CurrentDipoleMoment`:\n >>> cdp = CurrentDipoleMoment(cell)\n >>> M_I_to_P = cdp.get_transformation_matrix()\n >>> # instantiate class `eegmegcalc.InfiniteVolumeConductor` and map dipole moment to\n >>> # extracellular potential at measurement sites\n >>> ivc = eegmegcalc.InfiniteVolumeConductor(sigma=0.3)\n >>> # compute linear response matrix between dipole moment and\n >>> # extracellular potential\n >>> M_P_to_V = ivc.get_transformation_matrix(np.c_[np.ones(10) * 1000,\n np.zeros(10),\n np.arange(10) * 100])\n >>> # transmembrane currents (nA):\n >>> imem = np.array([[-1., 1.],\n [0., 0.],\n [1., -1.]])\n >>> # compute extracellular potentials (mV)\n >>> V_ex = M_P_to_V @ M_I_to_P @ imem\n >>> V_ex\n array([[ 0.00000000e+00, 0.00000000e+00],\n [ 5.22657054e-07, -5.22657054e-07],\n [ 1.00041193e-06, -1.00041193e-06],\n [ 1.39855769e-06, -1.39855769e-06],\n [ 1.69852477e-06, -1.69852477e-06],\n [ 1.89803345e-06, -1.89803345e-06],\n [ 2.00697409e-06, -2.00697409e-06],\n [ 2.04182029e-06, -2.04182029e-06],\n [ 2.02079888e-06, -2.02079888e-06],\n [ 1.96075587e-06, -1.96075587e-06]])\n\n\n## Physical units\n\nNotes on physical units used in `LFPykit`:\n\n- There are no explicit checks for physical units\n\n- Transmembrane currents are assumed to be in units of (nA)\n\n- Spatial information is assumed to be in units of (\u00b5m)\n\n- Voltages are assumed to be in units of (mV)\n\n- Extracellular conductivities are assumed to be in units of (S/m)\n\n- current dipole moments are assumed to be in units of (nA \u00b5m)\n\n- Magnetic fields are assumed to be in units of (nA/\u00b5m)\n\n\n## Dimensionality\n\n- Transmembrane currents are represented by arrays with shape `(n_seg, n_timesteps)`\n where `n_seg` is the number of segments of the neuron model.\n\n- Current dipole moments are represented by arrays with shape `(3, n_timesteps)`\n\n- Response matrices **M** have shape `(n_points, input.shape[0])` where `n_points` is\n for instance the number of extracellular recording sites and `input.shape[0]`\n the first dimension of the input; that is, the number of segments in case of\n transmembrane currents or 3 in case of current dipole moments.\n\n- predicted signals (except magnetic fields using `eegmegcalc.InfiniteHomogeneousVolCondMEG` or\n `eegmegcalc.SphericallySymmetricVolCondMEG`) have shape `(n_points, n_timesteps)`\n\n\n## Documentation\n\nThe online Documentation of `LFPykit` can be found here:\nhttps://lfpykit.readthedocs.io/en/latest\n\n\n## Dependencies\n\n`LFPykit` is implemented in Python and is written (and continuously tested) for `Python >= 3.7`.\nThe main `LFPykit` module depends on `numpy`, `scipy` and `MEAutility` (https://github.com/alejoe91/MEAutility, https://meautility.readthedocs.io/en/latest/).\n\nRunning all unit tests and example files may in addition require `py.test`, `matplotlib`,\n`neuron` (https://www.neuron.yale.edu),\n(`arbor` coming) and\n`LFPy` (https://github.com/LFPy/LFPy, https://LFPy.readthedocs.io).\n\n\n## Installation\n\n### From development sources (https://github.com/LFPy/LFPykit)\n\nInstall the current development version on https://GitHub.com using `git` (https://git-scm.com):\n\n $ git clone https://github.com/LFPy/LFPykit.git\n $ cd LFPykit\n $ python setup.py install # --user optional\n\nor using `pip`:\n\n $ pip install . # --user optional\n\nFor active development, link the repository location\n\n $ python setup.py develop # --user optional\n\n### Installation of stable releases on PyPI.org (https://www.pypi.org)\n\nInstalling from the Python Package Index (https://www.pypi.org/project/lfpykit):\n\n $ pip install lfpykit # --user optional\n\nTo upgrade the installation using pip:\n\n $ pip install --upgrade --no-deps lfpykit\n\n### Installation of stable releases on conda-forge (https://conda-forge.org)\n\nInstalling `lfpykit` from the `conda-forge` channel can be achieved by adding `conda-forge` to your channels with:\n\n $ conda config --add channels conda-forge\n\nOnce the `conda-forge` channel has been enabled, `lfpykit` can be installed with:\n\n $ conda install lfpykit\n\nIt is possible to list all of the versions of `lfpykit` available on your platform with:\n\n $ conda search lfpykit --channel conda-forge\n",
"bugtrack_url": null,
"license": "",
"summary": "Electrostatic models for multicompartment neuron models",
"version": "0.5.1",
"split_keywords": [],
"urls": [
{
"comment_text": "",
"digests": {
"blake2b_256": "6bbe5a72ee4deb4f1410e3266874447a38adde6f56cddf3e958506647f2e3ecc",
"md5": "7cf77f4e1455065591a788b0c20f0b8d",
"sha256": "25509ac8d3692c31231551d30c88720b49141fe6b8ed217efc54b0c59c68a076"
},
"downloads": -1,
"filename": "LFPykit-0.5.1-py3-none-any.whl",
"has_sig": false,
"md5_digest": "7cf77f4e1455065591a788b0c20f0b8d",
"packagetype": "bdist_wheel",
"python_version": "py3",
"requires_python": ">=3.7",
"size": 71217,
"upload_time": "2023-04-12T07:10:52",
"upload_time_iso_8601": "2023-04-12T07:10:52.625789Z",
"url": "https://files.pythonhosted.org/packages/6b/be/5a72ee4deb4f1410e3266874447a38adde6f56cddf3e958506647f2e3ecc/LFPykit-0.5.1-py3-none-any.whl",
"yanked": false,
"yanked_reason": null
},
{
"comment_text": "",
"digests": {
"blake2b_256": "b4e0b5935801a38839b177267f53cf6781941e36d8f832e06612773fdf89504e",
"md5": "71ba9f7b5232f335de2e2a75f438e703",
"sha256": "a962029460c8173c4fec3923204e04d64c370b05740b64c2d07efbe29cbe63a6"
},
"downloads": -1,
"filename": "LFPykit-0.5.1.tar.gz",
"has_sig": false,
"md5_digest": "71ba9f7b5232f335de2e2a75f438e703",
"packagetype": "sdist",
"python_version": "source",
"requires_python": ">=3.7",
"size": 2953825,
"upload_time": "2023-04-12T07:10:54",
"upload_time_iso_8601": "2023-04-12T07:10:54.970672Z",
"url": "https://files.pythonhosted.org/packages/b4/e0/b5935801a38839b177267f53cf6781941e36d8f832e06612773fdf89504e/LFPykit-0.5.1.tar.gz",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2023-04-12 07:10:54",
"github": true,
"gitlab": false,
"bitbucket": false,
"github_user": "LFPy",
"github_project": "LFPykit",
"travis_ci": false,
"coveralls": false,
"github_actions": true,
"requirements": [],
"lcname": "lfpykit"
}