This package was developed at the Wyss Institute at Harvard by A. Burykin, S. Mariani, T. Silva and T. Henriques. It is described in:
Burykin A, Mariani S, Henriques T, Silva T, Schnettler W, Costa MD, Goldberger AL. “Remembrance of time series past: simple chromatic method for visualizing trends in biomedical signals.” Physiol Meas 2015;36(7):N95.Please cite the above publication when referencing this material, and also include the standard citation for PhysioNet: Goldberger AL, Amaral LAN, Glass L, Hausdorff JM, Ivanov PCh, Mark RG, Mietus JE, Moody GB, Peng CK, Stanley HE. "PhysioBank, PhysioToolkit, and PhysioNet: Components of a New Research Resource for Complex Physiologic Signals," Circulation 101(23):e215e220 [Circulation Electronic Pages; http://circ.ahajournals.org/content/101/23/e215.full]; 2000 (June 13).
Description
Analysis of biomedical time series plays a key role in clinical management and basic investigation. However, most conventional monitors streaming data in realtime show only the most recent values, not referenced to past dynamics. The proposed visualization method (termed “data chromatix”) was developed to address this challenge by bringing memory of the system’s past behavior into the current display window.
The function DataChromatix.m (version 1.0) assigns a color to each data point of a time series. The color is determined by the values of a normalized histogram (estimated probability density function) computed from a preselected segment of the data. The algorithm receives the time series as input and generates a video of its colorized version as it would look on a typical monitor display, as well as a static graph of the entire colorized signal.
The algorithm has the following parameters: the memory and colorization window lengths, the shift (s), the histogram bin size, and the number of colors (c) in the chromatic map (Fig. 1).
Figure 1: Memory window and colorization window at a given time tC and at the following time instant tC+s. Both windows are shifted to the right by the shift s. Adapted from (1).
At each step, a normalized histogram of the data points in the memory window is computed. Then, the interval [0,1] is divided into c adjacent intervals, and each interval is assigned a color. If the jet colormap is used, the interval [0,1/c) corresponds to dark red and the interval ((c1)/c,1] corresponds to dark blue. Subsequently, each data point in the colorization window is assigned the color of the histogram bin into which it falls. Finally, the colorization window is advanced by the shift, s, and the memory window is either extended or advanced by the same amount (Fig. 2).
Figure 2: Schematic illustration of the colorization algorithm: the histogram of the time series (left panel) is colorized according to the bins height. Each point of the time series (right panel) is then colorized according to the bin it belongs to. Adapted from (1).
This colorization algorithm is intended to facilitate analysis of physiologic and nonphysiologic time series. Future studies will help assess its utility.
The function has been tested on MATLAB R2014a, R2014b and R2015a.
Sample Input
Mandatory inputs to the function are:
 The time series to be colorized
 The units to show on the x and y axis
 The length of the window to display in the video
 The shift
Optional inputs to the function are:
 The length of the memory window (default from the beginning of recording to the end of the colorization window)
 The length of the colorization window (default equal to the display window)
 The number of bins of the histogram (default determined by means of FriedmanDiaconis rule)
 The number of colors in the colormap (default = 64)
 The smoothing parameter for the histogram (default = 10)
 The data sampling frequency (default = 1 Hz)
 The name of the video (if = 0 does not create the video, default = ‘Myvideo’)
 A flag for data time format, when time units are seconds (if =1, the format is mm:ss, if =0 the format is s, default = 0)
 The desired colormap (default jet)
Along with the function, we provide two examples that employ a fetal heart rate time series from the CTUUHB Intrapartum Cardiotocography Database on PhysioNet, one using a memory window starting from the beginning of the recording, the other using a moving memory window of fixed length. Please note that for loading the time series, the MATLAB version of the wfdb library (http://physionet.org/physiotools/matlab/wfdbappmatlab/) must be installed.
Acknowledgments
This package was developed at the Wyss Institute at Harvard by A. Burykin, S. Mariani, T. Silva and T. Henriques.
Users of our software should cite: (1) Burykin A*, Mariani S*, Henriques T, Silva T, Schnettler W, Costa MD**, Goldberger AL**. “Remembrance of time series past: simple chromatic method for visualizing trends in biomedical signals.” Physiol Meas 2015;36(7):N95.
* Joint first authors
** Joint senior authors
Name Last modified Size Description
example2.m 20150612 17:56 485 example1.m 20150612 17:56 519 DataChromatix.m 20150612 17:56 9.1K
If you would like help understanding, using, or downloading content, please see our Frequently Asked Questions. If you have any comments, feedback, or particular questions regarding this page, please send them to the webmaster. Comments and issues can also be raised on PhysioNet's GitHub page. Updated Friday, 28 October 2016 at 16:58 EDT 
PhysioNet is supported by the National Institute of General Medical Sciences (NIGMS) and the National Institute of Biomedical Imaging and Bioengineering (NIBIB) under NIH grant number 2R01GM10498709.
