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

BrowZen correlates your emotional states with the web sites you visit to give you actionable insights about how you spend your time browsing the web.

Productivity Tools for Plotly + Pandas

Create 3d loss surface visualizations, with optimizer path. Issues welcome!

Minimalistic tool to visualize how the routes to a given target domain change over time, feat. Python 3.10 & mermaid.js

A shimmer pre-load component for Plotly Dash

📊 Charts with pure python

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

Set of matplotlib operations that are not trivial

Pretty Confusion Matrix

A Python package that provides evaluation and visualization tools for the DexYCB dataset

Automatically generate GitHub activity!

A simple python tool for explore your object detection dataset

This is a small program that prints a user friendly, visual representation, of your current bsp tree

Data Visualizer for Super Mario Kart (SNES)

Statistics and Visualization of acceptance rate, main keyword of CVPR 2021 accepted papers for the main Computer Vision conference (CVPR)

MPL Plotter is a Matplotlib based Python plotting library built with the goal of delivering publication-quality plots concisely.

Create a table with row explanations, column headers, using matplotlib

The Python ensemble sampling toolkit for affine-invariant MCMC

An interactive dashboard built with python that enables you to visualise how rent prices differ across Sweden.

nptsne is a numpy compatible python binary package that offers a number of APIs for fast tSNE calculation.