color.C

Go to the documentation of this file.

#ifdef HAVE_CONFIG_H
#include <tumble-conf.h>
#endif /* HAVE_CONFIG_H */

#include "color.h"

using namespace std;

Color RED( 1.0, 0.0, 0.0 );
Color BLUE( 0.0, 0.0, 1.0 );
Color GREEN( 0.0, 1.0, 0.0 );
Color YELLOW(1.0, 1.0, 0.0);
Color CYAN(0.0,1.0,1.0);
Color MAGENTA(1.0,0.0,1.0);
Color WHITE(1.0, 1.0, 1.0);
Color BLACK(0.0,0.0,0.0);
Color GRAY1(0.2, 0.2, 0.2);
Color GRAY2(0.4, 0.4, 0.4);
Color GRAY3(0.6, 0.6, 0.6);
Color GRAY4(0.8, 0.8, 0.8);

Color DIMYELLOW(0.2, 0.2, 0.0);
Color DIMCYAN(0.0, 0.2, 0.2);
Color DIMMAGENTA(0.2, 0.0, 0.2);
Color LIGHTYELLOW(1.0,1.0,0.5);
Color PINK1(1.0,0.6,0.6);
Color PINK2(1.0,0.3,0.3);
Color DARKRED(0.8,0.0,0.0);
Color DARKRED2(0.4, 0.0, 0.0);

map<int, Color> opengl_color_map;
//map<int, Color> eps_rgb_map;
//map<int, float> eps_grayscale_map;

bool Color::is_valid() const
{
    if( v[0] < 0.0 || v[0] > 1.0 ) return false;
    if( v[1] < 0.0 || v[1] > 1.0 ) return false;
    if( v[2] < 0.0 || v[2] > 1.0 ) return false;
    return true;
}

Color spectrum_color(double data, double data_min, double data_max)
{
    double norm;
    double variation = data_max - data_min;

    if(variation == 0.0) norm = 0.0;
    else norm= 1.0/variation;

    data = (data - data_min) * norm;

    Color LBLUE(0.5, 0.5, 1.0);

    int num_breaks = 7;
    Color colors[7] = {BLACK, MAGENTA, BLUE, CYAN, GREEN, YELLOW, RED};
    double breaks[7]= {0.0, 0.16666, 0.333333, 0.5, 0.6666666, 0.83333333, 1.0};

    /*
      Multiscale colors
      There should probably be a parameter to set these
    int num_breaks = 10;
    Color colors[10] =  {BLACK, MAGENTA,BLUE,   CYAN,   GREEN,  YELLOW,RED,   MAGENTA,LBLUE, WHITE};
    double breaks[10] = {0.0,   1.0E-12,5.0E-12,4.0E-11,2.0E-10,2.0E-9,9.0E-9,1.0E-4, 9.0E-4,1.0E-2};
    */

    Color rval = BLACK;

    if (data <= breaks[0]) return colors[0];
    if (data >= breaks[num_breaks-1]) return colors[num_breaks-1];

    int glb = 0; // greatest lower bound
    for (int i = 1; i<num_breaks-1; i++)
      {
    if (breaks[i] <= data) glb = i;
      }

    double mag = (data - breaks[glb]) / (breaks[glb+1] - breaks[glb]);
    rval = (colors[glb+1] * mag) + (colors[glb] * (1.0-mag));

    return rval;
}

/*
  old color map :  purple,blue,cyan,green,yellow,red
    if (data >= 1.0){
      //data is above max
      r = 1.0;
      g = 1.0;
      b = 1.0;
    }else if (data >= breaks[5]){
      //red to white
      r = 1.0;
      g = (data-breaks[5])/(1.0-breaks[5]);
      b = (data-breaks[5])/(1.0-breaks[5]);
    }else if (data >= breaks[4]){
      //yellow to red
      r = 1.0;
      g = (breaks[5] - data)/(breaks[5] - breaks[4]);
    } else if (data >= breaks[3]){
      //green to yellow
      r = (data - breaks[3])/(breaks[4]-breaks[3]);
      g = 1.0;
    } else if (data >= breaks[2]){
      //cyan to green
      g = 1.0;
      b = (breaks[3]-data)/(breaks[3]-breaks[2]);
    } else if (data >= breaks[1]){
      //blue to cyan
      b= 1.0;
      g= (data - breaks[1])/(breaks[2]-breaks[1]);
    } else if (data >= breaks[0]) {
      //purple to blue
      b = 1.0;
      r = (breaks[1] - data)/(breaks[1]-breaks[0]);
    } else if (data > 0.0){
      //black to purple
      b = (1.0 * data)/(breaks[0]);
      r = (1.0 * data)/(breaks[0]);
    }
*/


//#ifndef SWIG
void default_opengl_color_map( ColorMap &opengl_color_map )
{
    opengl_color_map[0]  = LIGHTYELLOW;
    opengl_color_map[1]  = DARKRED;
    opengl_color_map[2]  = PINK1;
    opengl_color_map[3]  = PINK2;
    opengl_color_map[4]  = CYAN;
    opengl_color_map[5]  = MAGENTA;
    opengl_color_map[6]  = BLACK;
    opengl_color_map[7]  = WHITE;

    opengl_color_map[8]  = DARKRED;
    opengl_color_map[9]  = DARKRED;
    opengl_color_map[10] = DARKRED;
    opengl_color_map[11] = DARKRED;
    opengl_color_map[12] = DARKRED;
    opengl_color_map[13] = DARKRED;
    opengl_color_map[14] = DARKRED;
    opengl_color_map[15] = DARKRED;
}

void default_eps_rgb_map( ColorMap &eps_rgb_map)
{
    eps_rgb_map[0]  = LIGHTYELLOW;
    eps_rgb_map[1]  = DARKRED;
    eps_rgb_map[2]  = PINK1;
    eps_rgb_map[3]  = PINK2;
    eps_rgb_map[4]  = CYAN;
    eps_rgb_map[5]  = MAGENTA;
    eps_rgb_map[6]  = BLACK;
    eps_rgb_map[7]  = WHITE;

    eps_rgb_map[8]  = DARKRED;
    eps_rgb_map[9]  = DARKRED;
    eps_rgb_map[10] = DARKRED;
    eps_rgb_map[11] = DARKRED;
    eps_rgb_map[12] = DARKRED;
    eps_rgb_map[13] = DARKRED;
    eps_rgb_map[14] = DARKRED;
    eps_rgb_map[15] = DARKRED;
}

void default_eps_grayscale_map( GrayscaleMap &eps_grayscale_map)
{
    eps_grayscale_map[0]  = 0.8;
    eps_grayscale_map[1]  = 0.6;
    eps_grayscale_map[2]  = 0.4;
    eps_grayscale_map[3]  = 0.2;
    eps_grayscale_map[4]  = 0.0;
    eps_grayscale_map[5]  = 0.1;
    eps_grayscale_map[6]  = 1.0;
    eps_grayscale_map[7]  = 1.0;

    eps_grayscale_map[8]  = 1.0;
    eps_grayscale_map[9]  = 1.0;
    eps_grayscale_map[10] = 1.0;
    eps_grayscale_map[11] = 1.0;
    eps_grayscale_map[12] = 1.0;
    eps_grayscale_map[13] = 1.0;
    eps_grayscale_map[14] = 1.0;
    eps_grayscale_map[15] = 1.0;
}

#if 0
int load_map( ColorMap &map)
{
    /* fix me */
    return 0;
}

int load_map( GrayscaleMap &map)
{
    /* fix me */
    return 0;
}
#endif
//#endif /* SWIG */

---