ts2vg 1.2.1

Build visibility graphs from time series data.

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.

The visibility graphs are provided according to the mathematical definitions presented 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.

Supported graph types

Root graph types

• Natural Visibility Graphs (NVG) (ts2vg.NaturalVG)

• Horizontal Visibility Graphs (HVG) (ts2vg.HorizontalVG)

Available variations

• Weighted Visibility Graphs (via the weighted parameter)

• Directed Visibility Graphs (via the directed parameter)

• Parametric Visibility Graphs (via the min_weight and max_weight parameters)

• Limited Penetrable Visibility Graphs (LPVG) (via the penetrable_limit parameter)

Note that multiple graph variations can be combined and used simultaneously.

Documentation

Usage and reference documentation for ts2vg can be found at carlosbergillos.github.io/ts2vg.

Basic usage

Visibility graph

To build a visibility graph:

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 (ts) can be any one-dimensional iterable, such as a list or a numpy 1D array.

By default, the input observations are assumed to be equally spaced in time. Alternatively, a second 1D iterable (xs) can be provided for unevenly spaced time series.

Horizontal visibility graph

Horizontal visibility graphs can be obtained 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 storing the whole graph.

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

ks, ps = g.degree_distribution

Directed visibility graph

Directed graphs can be obtained by using the directed parameter. See the reference documentation for the different available options.

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

Weighted visibility graph

Weighted graphs can be obtained by using the weighted parameter. See the reference documentation for the different available options.

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.