我正在尝试将星星的B-V颜色指数转换为明显的RGB颜色。除了查找表和颜色渐变之外，似乎没有众所周知的算法可以做到这一点。
什么是B-V颜色指数？
这是天文学家分配给星星的数字，用于指示其明显的颜色。热星（低B-V）为蓝/紫色，而冷星（高B-V）为红色，中间为白/橙色的星星。
初始算法
B-V到开尔文的转换
更新
根据Ozan的评论，我使用了错误的矩阵将XYZ转换为RGB。由于sRGB是Web上的默认色彩空间（或者不是吗？），我现在使用正确的矩阵，然后跟着进行伽马校正函数（
什么是B-V颜色指数？
这是天文学家分配给星星的数字，用于指示其明显的颜色。热星（低B-V）为蓝/紫色，而冷星（高B-V）为红色，中间为白/橙色的星星。
初始算法
B-V到开尔文的转换
var t = 4600 * ((1 / ((0.92 * bv) + 1.7)) +(1 / ((0.92 * bv) + 0.62)) );
开尔文转xyY
如果您将星体建模为黑体，则可以使用普朗克辐射定律的数值近似来计算xy坐标（CIE色度）。
// t to xyY
var x, y = 0;
if (t>=1667 && t<=4000) {
x = ((-0.2661239 * Math.pow(10,9)) / Math.pow(t,3)) + ((-0.2343580 * Math.pow(10,6)) / Math.pow(t,2)) + ((0.8776956 * Math.pow(10,3)) / t) + 0.179910;
} else if (t > 4000 && t <= 25000) {
x = ((-3.0258469 * Math.pow(10,9)) / Math.pow(t,3)) + ((2.1070379 * Math.pow(10,6)) / Math.pow(t,2)) + ((0.2226347 * Math.pow(10,3)) / t) + 0.240390;
}
if (t >= 1667 && t <= 2222) {
y = -1.1063814 * Math.pow(x,3) - 1.34811020 * Math.pow(x,2) + 2.18555832 * x - 0.20219683;
} else if (t > 2222 && t <= 4000) {
y = -0.9549476 * Math.pow(x,3) - 1.37418593 * Math.pow(x,2) + 2.09137015 * x - 0.16748867;
} else if (t > 4000 && t <= 25000) {
y = 3.0817580 * Math.pow(x,3) - 5.87338670 * Math.pow(x,2) + 3.75112997 * x - 0.37001483;
}
xyY 转 XYZ (Y = 1)
// xyY to XYZ, Y = 1
var Y = (y == 0)? 0 : 1;
var X = (y == 0)? 0 : (x * Y) / y;
var Z = (y == 0)? 0 : ((1 - x - y) * Y) / y;
XYZ转RGB
var r = 0.41847 * X - 0.15866 * Y - 0.082835 * Z;
var g = -0.091169 * X + 0.25243 * Y + 0.015708 * Z;
var b = 0.00092090 * X - 0.0025498 * Y + 0.17860 * Z;
问题
我使用B-V颜色指数运行了这个算法：1.2、1.0、0.59、0.0、-0.29。以下是输出结果。
更新
根据Ozan的评论，我使用了错误的矩阵将XYZ转换为RGB。由于sRGB是Web上的默认色彩空间（或者不是吗？），我现在使用正确的矩阵，然后跟着进行伽马校正函数（
a = 0.055）。
我现在得到了这个漂亮的颜色渐变,
演示
现在还有一个fiddle，您可以尝试一下。
更新2
如果使用0.5的Gamma值，并将B-V颜色指数的范围扩展到4.7至-0.5，则可以在一个极端获得红色，但仍无法获得紫色。这是更新后的fiddle链接。
![BV=f(temp[K])](https://istack.dev59.com/NPuaq.webp)
temp <--> BV（包括开尔文和摄氏度）之间的转换例程，以及RGB -> BV的转换例程... - Spektre