forked from widgetti/ipyvolume
-
Notifications
You must be signed in to change notification settings - Fork 0
/
transferfunction.py
138 lines (119 loc) · 6.16 KB
/
transferfunction.py
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
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
"""The transferfunction module of ipvyolume."""
from __future__ import absolute_import
__all__ = ['TransferFunction', 'TransferFunctionJsBumps', 'TransferFunctionWidgetJs3', 'TransferFunctionWidget3']
import numpy as np
import ipywidgets as widgets # we should not have widgets under two names
import traitlets
from traitlets import Unicode
from traittypes import Array
import ipyvolume._version
from ipyvolume import serialize
N = 1024
x = np.linspace(0, 1, N, endpoint=True)
semver_range_frontend = "~" + ipyvolume._version.__version_js__
@widgets.register
class TransferFunction(widgets.DOMWidget):
_model_name = Unicode('TransferFunctionModel').tag(sync=True)
_view_name = Unicode('TransferFunctionView').tag(sync=True)
_model_module = Unicode('ipyvolume').tag(sync=True)
_view_module = Unicode('ipyvolume').tag(sync=True)
style = Unicode("height: 32px; width: 100%;").tag(sync=True)
rgba = Array(default_value=None, allow_none=True).tag(sync=True, **serialize.ndarray_serialization)
_view_module_version = Unicode(semver_range_frontend).tag(sync=True)
_model_module_version = Unicode(semver_range_frontend).tag(sync=True)
class TransferFunctionJsBumps(TransferFunction):
_model_name = Unicode('TransferFunctionJsBumpsModel').tag(sync=True)
_model_module = Unicode('ipyvolume').tag(sync=True)
levels = traitlets.List(traitlets.CFloat(), default_value=[0.1, 0.5, 0.8]).tag(sync=True)
opacities = traitlets.List(traitlets.CFloat(), default_value=[0.01, 0.05, 0.1]).tag(sync=True)
widths = traitlets.List(traitlets.CFloat(), default_value=[0.1, 0.1, 0.1]).tag(sync=True)
def control(self, max_opacity=0.2):
return widgets.VBox()
class TransferFunctionWidgetJs3(TransferFunction):
_model_name = Unicode('TransferFunctionWidgetJs3Model').tag(sync=True)
_model_module = Unicode('ipyvolume').tag(sync=True)
level1 = traitlets.Float(0.1).tag(sync=True)
level2 = traitlets.Float(0.5).tag(sync=True)
level3 = traitlets.Float(0.8).tag(sync=True)
opacity1 = traitlets.Float(0.01).tag(sync=True)
opacity2 = traitlets.Float(0.05).tag(sync=True)
opacity3 = traitlets.Float(0.1).tag(sync=True)
width1 = traitlets.Float(0.1).tag(sync=True)
width2 = traitlets.Float(0.1).tag(sync=True)
width3 = traitlets.Float(0.1).tag(sync=True)
def control(self, max_opacity=0.2):
l1 = widgets.FloatSlider(min=0, max=1, step=0.001, value=self.level1)
l2 = widgets.FloatSlider(min=0, max=1, step=0.001, value=self.level2)
l3 = widgets.FloatSlider(min=0, max=1, step=0.001, value=self.level3)
o1 = widgets.FloatSlider(min=0, max=max_opacity, step=0.001, value=self.opacity1)
o2 = widgets.FloatSlider(min=0, max=max_opacity, step=0.001, value=self.opacity2)
o3 = widgets.FloatSlider(min=0, max=max_opacity, step=0.001, value=self.opacity2)
widgets.jslink((self, 'level1'), (l1, 'value'))
widgets.jslink((self, 'level2'), (l2, 'value'))
widgets.jslink((self, 'level3'), (l3, 'value'))
widgets.jslink((self, 'opacity1'), (o1, 'value'))
widgets.jslink((self, 'opacity2'), (o2, 'value'))
widgets.jslink((self, 'opacity3'), (o3, 'value'))
return widgets.VBox(
[
widgets.HBox([widgets.Label(value="levels:"), l1, l2, l3]),
widgets.HBox([widgets.Label(value="opacities:"), o1, o2, o3]),
]
)
class TransferFunctionWidget3(TransferFunction):
level1 = traitlets.Float(0.1).tag(sync=True)
level2 = traitlets.Float(0.5).tag(sync=True)
level3 = traitlets.Float(0.8).tag(sync=True)
opacity1 = traitlets.Float(0.4).tag(sync=True)
opacity2 = traitlets.Float(0.1).tag(sync=True)
opacity3 = traitlets.Float(0.1).tag(sync=True)
width1 = traitlets.Float(0.1).tag(sync=True)
width2 = traitlets.Float(0.1).tag(sync=True)
width3 = traitlets.Float(0.1).tag(sync=True)
def __init__(self, *args, **kwargs):
super(TransferFunctionWidget3, self).__init__(*args, **kwargs)
N = range(1, 4)
self.observe(
self.recompute_rgba, ["level%d" % k for k in N] + ["opacity%d" % k for k in N] + ["width%d" % k for k in N]
)
self.recompute_rgba()
def recompute_rgba(self, *_ignore):
import matplotlib
rgba = np.zeros((1024, 4))
N = range(1, 4)
levels = [getattr(self, "level%d" % k) for k in N]
opacities = [getattr(self, "opacity%d" % k) for k in N]
widths = [getattr(self, "width%d" % k) for k in N]
# print(levels, opacities, widths)
colors = [np.array(matplotlib.colors.colorConverter.to_rgb(name)) for name in ["red", "green", "blue"]]
# TODO: vectorize
for i in range(rgba.shape[0]):
position = i / (1023.0)
intensity = 0.0
# self.rgba[i, 0, :] = 0
for j in range(3):
intensity = np.exp(-((position - levels[j]) / widths[j]) ** 2)
rgba[i, 0:3] += colors[j] * opacities[j] * intensity
rgba[i, 3] += opacities[j] * intensity
rgba[i, 0:3] /= rgba[i, 0:3].max()
rgba = np.clip(rgba, 0, 1)
self.rgba = rgba
def control(self, max_opacity=0.2):
l1 = widgets.FloatSlider(min=0, max=1, value=self.level1)
l2 = widgets.FloatSlider(min=0, max=1, value=self.level2)
l3 = widgets.FloatSlider(min=0, max=1, value=self.level3)
o1 = widgets.FloatSlider(min=0, max=max_opacity, step=0.001, value=self.opacity1)
o2 = widgets.FloatSlider(min=0, max=max_opacity, step=0.001, value=self.opacity2)
o3 = widgets.FloatSlider(min=0, max=max_opacity, step=0.001, value=self.opacity2)
widgets.jslink((self, 'level1'), (l1, 'value'))
widgets.jslink((self, 'level2'), (l2, 'value'))
widgets.jslink((self, 'level3'), (l3, 'value'))
widgets.jslink((self, 'opacity1'), (o1, 'value'))
widgets.jslink((self, 'opacity2'), (o2, 'value'))
widgets.jslink((self, 'opacity3'), (o3, 'value'))
return widgets.VBox(
[
widgets.HBox([widgets.Label(value="levels:"), l1, l2, l3]),
widgets.HBox([widgets.Label(value="opacities:"), o1, o2, o3]),
]
)