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

Generate the report for OCULTest.

Pebble is a stat's visualization tool, this will provide a skeleton to develop a monitoring tool.

Gallery of applications built using bqplot and widget libraries like ipywidgets, ipydatagrid etc.

Visualize your pandas data with one-line code

Process dataframe in a easily way.

Here are my graphs for hw_02

Script to create an animated data visualisation for categorical timeseries data - GIF choropleth map with annotations.

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

Python histogram library - histograms as updateable, fully semantic objects with visualization tools. [P]ython [HYST]ograms.

Make scripted visualizations in blender

Rick and Morty Data Visualization with python

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

Flexitext is a Python library that makes it easier to draw text with multiple styles in Matplotlib

Bcc2telegraf: An integration that sends ebpf-based bcc histogram metrics to telegraf daemon

Insert SVGs into matplotlib

Application for viewing pokemon regional variants.

Fast 1D and 2D histogram functions in Python

Curvipy - The Python package for visualizing curves and linear transformations in a super simple way

A program that analyzes data from inertia measurement units installed in aircraft and generates g-exceedance curves.

Drag’n’drop Pivot Tables and Charts for Jupyter/IPython Notebook, care of PivotTable.js