PyMassSpec¶
Python Toolkit for Mass Spectrometry
Docs |
|
---|---|
Tests |
|
PyPI |
|
Anaconda |
|
Activity |
|
QA |
|
Other |
PyMassSpec is a Python package for processing gas chromatography-mass spectrometry data. PyMassSpec provides a framework and a set of components for rapid development and testing of methods for processing of chromatography–mass spectrometry data. PyMassSpec can be used interactively through the Python shell, in a Jupyter Notebook, or the functions can be collected into scripts when it is preferable to perform data processing in the batch mode.
Forked from the original PyMS Repository: https://github.com/ma-bio21/pyms. Originally by Andrew Isaac, Sean O’Callaghan and Vladimir Likić. The original publication can be found here: https://bmcbioinformatics.biomedcentral.com/articles/10.1186/1471-2105-13-115
The original project seems to have been abandoned as there has been no activity since 2017.
The PyMassSpec project¶
The directory structure of PyMassSpec is as follows:
/
├── pyms: The PyMassSpec code
│
├── pyms-data: Example GC-MS data files
│
├── pyms-demo: Examples of how to use PyMassSpec
│
├── tests: pytest tests
│
└── doc-source: Sphinx source for documentation
Features¶
Installation¶
python3 -m pip install PyMassSpec --user
First add the required channels
conda config --add channels https://conda.anaconda.org/bioconda
conda config --add channels https://conda.anaconda.org/conda-forge
conda config --add channels https://conda.anaconda.org/domdfcoding
Then install
conda install PyMassSpec
python3 -m pip install git+https://github.com/PyMassSpec/PyMassSpec@master --user
Usage¶
A tutorial illustrating various PyMassSpec features in detail is provided in subsequent chapters of this User Guide. The commands executed interactively are grouped together by example, and can be found here.
The data used in the PyMassSpec documentation and examples is available here.
In the “Demos and Examples” section there is a page corresponding to each example, coded with the chapter number (ie. “pyms-demo/20a/” corresponds to the Example 20a, from Chapter 2).
Each example has a script named ‘proc.py’ which contains the commands given in the example. These scripts can be run with the following command:
$ python3 proc.py
Example processing GC-MS data¶
Download the file gc01_0812_066.jdx
and save it in the folder data
.
This file contains GC-MS data in the the JCAMP-DX format.
First the raw data is loaded:
>>> from pyms.GCMS.IO.JCAMP import JCAMP_reader
>>> jcamp_file = "data/gc01_0812_066.jdx"
>>> data = JCAMP_reader(jcamp_file)
-> Reading JCAMP file 'Data/gc01_0812_066.jdx'
>>> data
<pyms.GCMS.Class.GCMS_data at 0x7f3ec77da0b8>
The intensity matrix object is then built by binning the data:
>>> from pyms.IntensityMatrix import build_intensity_matrix_i
>>> im = build_intensity_matrix_i(data)
In this example, we show how to obtain the dimensions of the newly created intensity matrix, then loop over all ion chromatograms, and for each ion chromatogram apply Savitzky-Golay noise filter and tophat baseline correction:
>>> n_scan, n_mz = im.size
>>> from pyms.Noise.SavitzkyGolay import savitzky_golay
>>> from pyms.TopHat import tophat
>>> for ii in range(n_mz):
... print("working on IC", ii)
... ic = im.get_ic_at_index(ii)
... ic1 = savitzky_golay(ic)
... ic_smooth = savitzky_golay(ic1)
... ic_base = tophat(ic_smooth, struct="1.5m")
... im.set_ic_at_index(ii, ic_base)
The resulting noise and baseline corrected ion chromatogram is saved back into the intensity matrix.
Further examples can be found in the documentation
License¶
PyMassSpec is Free and Open Source software released under the GNU General Public License version 2.
Issues¶
If you encounter any problems, please file an issue along with a detailed description.
Installation¶
python3 -m pip install PyMassSpec --user
First add the required channels
conda config --add channels https://conda.anaconda.org/bioconda
conda config --add channels https://conda.anaconda.org/conda-forge
conda config --add channels https://conda.anaconda.org/domdfcoding
Then install
conda install PyMassSpec
python3 -m pip install git+https://github.com/PyMassSpec/PyMassSpec@master --user
User Guide
Documentation
Demos and Examples
Contributing
View the Function Index or browse the Source Code.