pcolor/cmap.h at main · ggdoe/pcolor · GitHub

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
#pragma once
#include <math.h>
#include <stdint.h>

uint32_t cmap_cool_warm(double v)
{
    static float data[4][3] =
            {
                {59.000115, 75.999945, 191.999955},
                {220.57500, 220.57500, 220.575000},
                {179.99991, 4.0000070, 38.0001000},
                {179.99991, 4.0000070, 38.0001000},
            };

    v = (v < 0.0) ? 0.0 : (v > 1.0) ? 1.0 : v;
    double index = v * 2.;
    double frac = index - floor(index);
    size_t id = index;

    uint32_t r = ((1-frac) * data[id][0] + frac * data[id+1][0]);
    uint32_t g = ((1-frac) * data[id][1] + frac * data[id+1][1]);
    uint32_t b = ((1-frac) * data[id][2] + frac * data[id+1][2]);

    return (r<<16) | (g<<8) | b;
}

uint32_t cmap_nipy_spectral(double v)
{
    static float data[22][3] =
            {
                {  0.,   0.,   0.},
                {119.,   0., 136.},
                {136.,   0., 153.},
                {  0.,   0., 170.},
                {  0.,   0., 221.},
                {  0., 119., 221.},
                {  0., 153., 221.},
                {  0., 170., 170.},
                {  0., 170., 136.},
                {  0., 153.,   0.},
                {  0., 187.,   0.},
                {  0., 221.,   0.},
                {  0., 255.,   0.},
                {187., 255.,   0.},
                {238., 238.,   0.},
                {255., 204.,   0.},
                {255., 153.,   0.},
                {255.,   0.,   0.},
                {221.,   0.,   0.},
                {204.,   0.,   0.},
                {204., 204., 204.},
                {204., 204., 204.}
            };

    v = (v < 0.0) ? 0.0 : (v > 1.0) ? 1.0 : v;
    double index = v * 20.;
    double frac = index - floor(index);
    size_t id = index;

    uint32_t r = ((1-frac) * data[id][0] + frac * data[id+1][0]);
    uint32_t g = ((1-frac) * data[id][1] + frac * data[id+1][1]);
    uint32_t b = ((1-frac) * data[id][2] + frac * data[id+1][2]);

    return (r<<16) | (g<<8) | b;
}