Python ts2vg package provides high-performance algorithm implementations to build visibility graphs from time series data.

Last update: Dec 17, 2022

Overview

ts2vg: Time series to visibility graphs

The Python ts2vg package provides high-performance algorithm implementations to build visibility graphs from time series data.

The visibility graphs and some of their properties (e.g. degree distributions) are computed quickly and efficiently, even for time series with millions of observations thanks to the use of NumPy and a custom C backend (via Cython) developed for the visibility algorithms.

The visibility graphs are provided according to the mathematical definitions described in:

Lucas Lacasa et al., "From time series to complex networks: The visibility graph", 2008.
Lucas Lacasa et al., "Horizontal visibility graphs: exact results for random time series", 2009.

An efficient divide-and-conquer algorithm is used to compute the graphs, as described in:

Xin Lan et al., "Fast transformation from time series to visibility graphs", 2015.

Installation

The latest released ts2vg version is available at the Python Package Index (PyPI) and can be easily installed by running:

pip install ts2vg

For other advanced uses, to build ts2vg from source Cython is required.

Basic usage

Visibility graph

Building visibility graphs from time series is very simple:

from ts2vg import NaturalVG

ts = [1.0, 0.5, 0.3, 0.7, 1.0, 0.5, 0.3, 0.8]

g = NaturalVG()
g.build(ts)

edges = g.edges

The time series passed can be a list, a tuple, or a numpy 1D array.

Horizontal visibility graph

We can also obtain horizontal visibility graphs in a very similar way:

from ts2vg import HorizontalVG

ts = [1.0, 0.5, 0.3, 0.7, 1.0, 0.5, 0.3, 0.8]

g = HorizontalVG()
g.build(ts)

edges = g.edges

Degree distribution

If we are only interested in the degree distribution of the visibility graph we can pass only_degrees=True to the build method. This will be more efficient in time and memory than computing the whole graph.

g = NaturalVG()
g.build(ts, only_degrees=True)

ks, ps = g.degree_distribution

Directed visibility graph

g = NaturalVG(directed='left_to_right')
g.build(ts)

Weighted visibility graph

g = NaturalVG(weighted='distance')
g.build(ts)

For more information and options see: Examples and API Reference.

Interoperability with other libraries

The graphs obtained can be easily converted to graph objects from other common Python graph libraries such as igraph, NetworkX and SNAP for further analysis.

The following methods are provided:

as_igraph()
as_networkx()
as_snap()

For example:

g = NaturalVG()
g.build(ts)

nx_g = g.as_networkx()

Command line interface

ts2vg can also be used as a command line program directly from the console:

ts2vg ./timeseries.txt -o out.edg

For more help and a list of options run:

ts2vg --help

Contributing

ts2vg can be found on GitHub. Pull requests and issue reports are welcome.

License

ts2vg is licensed under the terms of the MIT License.

You might also like...

The Spectral Diagram (SD) is a new tool for the comparison of time series in the frequency domain

The Spectral Diagram (SD) is a new tool for the comparison of time series in the frequency domain. The SD provides a novel way to display the coherence function, power, amplitude, phase, and skill score of discrete frequencies of two time series. Each SD summarises these quantities in a single plot for multiple targeted frequencies.

3 Oct 10, 2022

The windML framework provides an easy-to-use access to wind data sources within the Python world, building upon numpy, scipy, sklearn, and matplotlib. Renewable Wind Energy, Forecasting, Prediction

windml Build status : The importance of wind in smart grids with a large number of renewable energy resources is increasing. With the growing infrastr

125 Dec 24, 2022

Comments

help getting started

I am playing around with ts2vg and I am having a hard time with the plotting using igraph. I try to compute the natural vg for a short time series, but when trying to plot it I get this error:

Traceback (most recent call last):
  File "\anaconda3\envs\DK_01\lib\site-packages\IPython\core\interactiveshell.py", line 3398, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-1-9a1fdcf342e8>", line 1, in <cell line: 1>
    ig.plot(nx_g, target='graph.pdf')
  File "\anaconda3\envs\DK_01\lib\site-packages\igraph\drawing\__init__.py", line 512, in plot
    result.save()
  File "\anaconda3\envs\DK_01\lib\site-packages\igraph\drawing\__init__.py", line 309, in save
    self._ctx.show_page()
igraph.drawing.cairo.MemoryError: out of memory

The file created is corrupted.

Here is my code:

import numpy as np
from ts2vg import NaturalVG
import igraph as ig

import matplotlib.pyplot as plt

# time domain
t = np.linspace(1, 40)
dt = np.diff(t)

# build series
x1 = np.sin(2*np.pi/10*t)
x2 = np.sin(2*np.pi/15*t)

y = x1 + x2

plt.plot(t, y, '.-')
plt.show()

# build HVG
g = NaturalVG()
g.build(y)

nx_g = g.as_igraph()

# plotting
ig.plot(nx_g, target='graph.pdf')

I am using ts2vg 1.0.0, igraph 0.9.11, and pycairo 1.21.0

opened by ACatAC 1

Python ts2vg package provides high-performance algorithm implementations to build visibility graphs from time series data.

Related tags

Overview

ts2vg: Time series to visibility graphs

Installation

Basic usage

Visibility graph

Horizontal visibility graph

Degree distribution

Directed visibility graph

Weighted visibility graph

Interoperability with other libraries

Command line interface

Contributing

License

You might also like...

The Spectral Diagram (SD) is a new tool for the comparison of time series in the frequency domain

The windML framework provides an easy-to-use access to wind data sources within the Python world, building upon numpy, scipy, sklearn, and matplotlib. Renewable Wind Energy, Forecasting, Prediction

Kglab - an abstraction layer in Python for building knowledge graphs

Extensible, parallel implementations of t-SNE

Extensible, parallel implementations of t-SNE

Graphical display tools, to help students debug their class implementations in the Carcassonne family of projects

Draw interactive NetworkX graphs with Altair

Draw interactive NetworkX graphs with Altair

Generate graphs with NetworkX, natively visualize with D3.js and pywebview

Comments

help getting started

Releases(v1.0.0)

v1.0.0(Aug 19, 2021)

Owner

Carlos Bergillos

A GUI for Pandas DataFrames

LinkedIn connections analyzer

A declarative (epi)genomics visualization library for Python

A research of IT labor market based especially on hh.ru. Salaries, rate of technologies and etc.

Boltzmann visualization - Visualize the Boltzmann distribution for simple quantum models of molecular motion

Some problems of SSLC ( High School ) before outputs and after outputs

Epagneul is a tool to visualize and investigate windows event logs

HW 02 for CS40 - matplotlib practice

Minimal Ethereum fee data viewer for the terminal, contained in a single python script.

Tools for writing, submitting, debugging, and monitoring Storm topologies in pure Python

Example scripts for generating plots of Bohemian matrices

Uniform Manifold Approximation and Projection

Generate SVG (dark/light) images visualizing (private/public) GitHub repo statistics for profile/website.

visualize_ML is a python package made to visualize some of the steps involved while dealing with a Machine Learning problem

Print matplotlib colors

A streamlit component for bi-directional communication with bokeh plots.

Create SVG drawings from vector geodata files (SHP, geojson, etc).

IPython/Jupyter notebook module for Vega and Vega-Lite

Simple function to plot multiple barplots in the same figure.

Data Visualizations for the #30DayChartChallenge