1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
|
# -*- coding: utf-8 -*-
# -----------------------------------------------------------------------------
# Copyright (c) Vispy Development Team. All Rights Reserved.
# Distributed under the (new) BSD License. See LICENSE.txt for more info.
# -----------------------------------------------------------------------------
import numpy as np
from .mesh import MeshVisual
class HistogramVisual(MeshVisual):
"""Visual that calculates and displays a histogram of data
Parameters
----------
data : array-like
Data to histogram. Currently only 1D data is supported.
bins : int | array-like
Number of bins, or bin edges.
color : instance of Color
Color of the histogram.
orientation : {'h', 'v'}
Orientation of the histogram.
"""
def __init__(self, data, bins=10, color='w', orientation='h'):
# 4-5
# | |
# 1-2/7-8
# |/| | |
# 0-3-6-9
data = np.asarray(data)
if data.ndim != 1:
raise ValueError('Only 1D data currently supported')
if not isinstance(orientation, str) or \
orientation not in ('h', 'v'):
raise ValueError('orientation must be "h" or "v", not %s'
% (orientation,))
X, Y = (0, 1) if orientation == 'h' else (1, 0)
# do the histogramming
data, bin_edges = np.histogram(data, bins)
# construct our vertices
rr = np.zeros((3 * len(bin_edges) - 2, 3), np.float32)
rr[:, X] = np.repeat(bin_edges, 3)[1:-1]
rr[1::3, Y] = data
rr[2::3, Y] = data
bin_edges.astype(np.float32)
# and now our tris
tris = np.zeros((2 * len(bin_edges) - 2, 3), np.uint32)
offsets = 3 * np.arange(len(bin_edges) - 1,
dtype=np.uint32)[:, np.newaxis]
tri_1 = np.array([0, 2, 1])
tri_2 = np.array([2, 0, 3])
tris[::2] = tri_1 + offsets
tris[1::2] = tri_2 + offsets
MeshVisual.__init__(self, rr, tris, color=color)
|
