sandbox/nlemoine/groundwater/Meuvaines/colormaps.h

    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
    
    int NCCLASS;
    
    void discrete_blues (double cmap[NCMAP][3])
    {
      double startColor[3] = {0.,0.,0.5};
      double endColor[3] = {0.95,0.95,1.};
      double classWidth = (-1.+(double)NCMAP)/((double)NCCLASS);
    
      for (int i = 0; i < NCMAP; i++) {
        double f = i<(NCMAP-1) ? floor( ((double) i) / classWidth ) / (-1+(double)NCCLASS) : 1.0;
        cmap[i][0] = f * endColor[0] + (1.-f) * startColor[0];
        cmap[i][1] = f * endColor[1] + (1.-f) * startColor[1];
        cmap[i][2] = f * endColor[2] + (1.-f) * startColor[2];
      }
    }
    
    void discrete_reds (double cmap[NCMAP][3])
    {
      double startColor[3] = {0.45,0.,0.};
      double endColor[3] = {1.,0.9,0.9};
      double classWidth = (-1.+(double)NCMAP)/((double)NCCLASS);
    
      for (int i = 0; i < NCMAP; i++) {
        double f = i<(NCMAP-1) ? floor( ((double) i) / classWidth ) / (-1+(double)NCCLASS) : 1.0;
        cmap[i][0] = f * endColor[0] + (1.-f) * startColor[0];
        cmap[i][1] = f * endColor[1] + (1.-f) * startColor[1];
        cmap[i][2] = f * endColor[2] + (1.-f) * startColor[2];
      }
    }
    
    void discrete_grays (double cmap[NCMAP][3])
    {
      double startColor[3] = {0.2,0.2,0.2};
      double endColor[3] = {0.9,0.9,0.9};
      double classWidth = (-1.+(double)NCMAP)/((double)NCCLASS);
    
      for (int i = 0; i < NCMAP; i++) {
        double f = i<(NCMAP-1) ? floor( ((double) i) / classWidth ) / (-1+(double)NCCLASS) : 1.0;
        cmap[i][0] = f * endColor[0] + (1.-f) * startColor[0];
        cmap[i][1] = f * endColor[1] + (1.-f) * startColor[1];
        cmap[i][2] = f * endColor[2] + (1.-f) * startColor[2];
      }
    }
    
    void continuous_blues (double cmap[NCMAP][3])
    {
      double startColor[3] = {0.,0.,0.5};
      double endColor[3] = {0.95,0.95,1.};
    
      for (int i = 0; i < NCMAP; i++) {
        double f = ((double) i) / (-1.+(double)NCMAP);
        cmap[i][0] = f * endColor[0] + (1.-f) * startColor[0];
        cmap[i][1] = f * endColor[1] + (1.-f) * startColor[1];
        cmap[i][2] = f * endColor[2] + (1.-f) * startColor[2];
      }
    }
    
    void continuous_reds (double cmap[NCMAP][3])
    {
      double startColor[3] = {0.45,0.,0.};
      double endColor[3] = {1.,0.9,0.9};
    
      for (int i = 0; i < NCMAP; i++) {
        double f = ((double) i) / (-1.+(double)NCMAP);
        cmap[i][0] = f * endColor[0] + (1.-f) * startColor[0];
        cmap[i][1] = f * endColor[1] + (1.-f) * startColor[1];
        cmap[i][2] = f * endColor[2] + (1.-f) * startColor[2];
      }
    }
    
    void continuous_grays (double cmap[NCMAP][3])
    {
      double startColor[3] = {0.2,0.2,0.2};
      double endColor[3] = {0.9,0.9,0.9};
    
      for (int i = 0; i < NCMAP; i++) {
        double f = ((double) i) / (-1.+(double)NCMAP);
        cmap[i][0] = f * endColor[0] + (1.-f) * startColor[0];
        cmap[i][1] = f * endColor[1] + (1.-f) * startColor[1];
        cmap[i][2] = f * endColor[2] + (1.-f) * startColor[2];
      }
    }
    
    void topocolormap (double cmap[NCMAP][3])
    {
       double classColors[5][3] = {{0.,0.,0.6},
    			     {0.6,0.6,1.},
                                 {1.,1.,1.},
                                 {1.,0.6,0.6},
                                 {0.6,0.,0.}};
    
      double classWidth = (-1.+(double)NCMAP)/5.;
    
      for (int i = 0; i < NCMAP; i++) {
        int nc = i<(NCMAP-1) ? (int) floor( ((double) i) / classWidth ) : 5;
        cmap[i][0] = classColors[nc][0];
        cmap[i][1] = classColors[nc][1];
        cmap[i][2] = classColors[nc][2];
      }
    
    }
    
    void uniform_gray(double cmap[NCMAP][3])
    {
      for (int i = 0; i < NCMAP; i++) {
        cmap[i][0] = 0.7;
        cmap[i][1] = 0.7;
        cmap[i][2] = 0.7;
      }
    }