灰度到红-绿-蓝 (MATLAB Jet) 色标

Question

我得到了一个本质上是图像的数据集，但是图像中的每个像素都表示为从 -1 到 1 的值。我正在编写一个应用程序，需要将这些 -1 到 1 灰度值映射到 MATLAB“Jet”色标（红-绿-蓝颜色渐变）的关联 RGB 值。

我很好奇是否有人知道如何取一个线性值（如 -1 到 1）并将其映射到这个比例。请注意，我实际上并没有为此使用 MATLAB（我也不能），我只需要获取灰度值并将其放在 Jet 渐变上。

谢谢，亚当

原文由 Adam Shook 发布，翻译遵循 CC BY-SA 4.0 许可协议

Answer 1

我希望这是您正在寻找的：

 double interpolate( double val, double y0, double x0, double y1, double x1 ) {
  return (val-x0)*(y1-y0)/(x1-x0) + y0;
}
double blue( double grayscale ) {
  if ( grayscale < -0.33 ) return 1.0;
  else if ( grayscale < 0.33 ) return interpolate( grayscale, 1.0, -0.33, 0.0, 0.33 );
  else return 0.0;
}
double green( double grayscale ) {
  if ( grayscale < -1.0 ) return 0.0; // unexpected grayscale value
  if  ( grayscale < -0.33 ) return interpolate( grayscale, 0.0, -1.0, 1.0, -0.33 );
  else if ( grayscale < 0.33 ) return 1.0;
  else if ( grayscale <= 1.0 ) return interpolate( grayscale, 1.0, 0.33, 0.0, 1.0 );
  else return 1.0; // unexpected grayscale value
}
double red( double grayscale ) {
  if ( grayscale < -0.33 ) return 0.0;
  else if ( grayscale < 0.33 ) return interpolate( grayscale, 0.0, -0.33, 1.0, 0.33 );
  else return 1.0;
}

我不确定这个比例是否与您链接的图像 100% 相同，但它应该看起来非常相似。

更新 我已经根据 此处 找到的 MatLab 的 Jet 调色板的描述重写了代码

double interpolate( double val, double y0, double x0, double y1, double x1 ) {
    return (val-x0)*(y1-y0)/(x1-x0) + y0;
}

double base( double val ) {
    if ( val <= -0.75 ) return 0;
    else if ( val <= -0.25 ) return interpolate( val, 0.0, -0.75, 1.0, -0.25 );
    else if ( val <= 0.25 ) return 1.0;
    else if ( val <= 0.75 ) return interpolate( val, 1.0, 0.25, 0.0, 0.75 );
    else return 0.0;
}

double red( double gray ) {
    return base( gray - 0.5 );
}
double green( double gray ) {
    return base( gray );
}
double blue( double gray ) {
    return base( gray + 0.5 );
}

原文由 Ilya Denisov 发布，翻译遵循 CC BY-SA 4.0 许可协议

Answer 2

将生成 Jet 和 HotAndCold RGB 的 Java（处理）代码。我按照上面 Amro 帖子中的 RGB 分配方案创建了这段代码。

 color JetColor(float v,float vmin,float vmax){
       float r=0, g=0, b=0;
       float x = (v-vmin)/(vmax-vmin);
       r = 255*constrain(-4*abs(x-0.75) + 1.5,0,1);
       g = 255*constrain(-4*abs(x-0.50) + 1.5,0,1);
       b = 255*constrain(-4*abs(x-0.25) + 1.5,0,1);
       return color(r,g,b);
    }

color HeatColor(float v,float vmin,float vmax){
       float r=0, g=0, b=0;
       float x = (v-vmin)/(vmax-vmin);
       r = 255*constrain(-4*abs(x-0.75) + 2,0,1);
       g = 255*constrain(-4*abs(x-0.50) + 2,0,1);
       b = 255*constrain(-4*abs(x) + 2,0,1);
       return color(r,g,b);
    }
    //Values are calculated on trapezoid cutoff points in format y=constrain(a(x-t)+b,0,1)
    //Where a=((delta)y/(delta)x), t=x-offset value to symetric middle of trapezoid, and b=y-a(x-t) for the last peak point (x,y)

原文由 schornagel86 发布，翻译遵循 CC BY-SA 4.0 许可协议