Berechnen Farbtemperatur in K

Wenn ich Sie richtig verstehe, sehen Sie sich bitte die pdf, wo eine Reihe von Methoden (kurz 🙁 ) beschrieben. Aus diesem PDF:

[CCT1] Berechnungen werden durchgeführt unter Verwendung der CIE 1960 uniform color space
die Koordinaten u und V Koordinaten u und v ergeben sich aus x und y mit der
Formeln:
u = 4x/(12y-2x+3) und v = 6y/(12y-2x+3)
Die Farbtemperatur ist definiert als die Temperatur eines schwarzen Körpers, die
liegt am nächsten an der u,v-Koordinaten der test-Quelle. Zwei separate Methoden
Ergebnisse: die eine ist eine iterative Methode basiert auf der definition, und die andere ist die
Häufig verwendete Robertson ' s Methode mit interpolation basierend auf einer Tabelle von 30
pre-berechnet u, v und inverse slope Parameter.
Iteration
Wenn uS und vS
sind die Werte für eine test-Quelle, sowie uT und vT
sind die schwarzen Werte
bei der Temperatur T, die korrelierte Farbtemperatur ist der Wert von T, wo:

sqrt( (uS - uT)^2 + (vS - vT)^2 )

minimiert wird. Einstellung der Temperatur T erhalten Sie das minimum dieser Funktion
wurde durchgeführt unter Verwendung von Tabellenkalkulation (Quattro Pro 8 und Excel 97). Beide Tabellen gab
identische Werte.

Weiß nicht, ob es tatsächlich hilft Ihnen.

Ich getan einige Graben in einige open-source-apps, und fand etwas, das in UFRaw. Ich habe noch nicht ganz herausgefunden, was Los ist genau.

Fand auch ein Papier das scheint so zu decken das Thema ganz gut.

Konvertiert php und das ist, was ich habe, so weit:

$temp = array(9500, 7000, 5500, 3750, 3000, 2700, 2250, 1800, 1500);
$hex = array('9DBEFF', 'E4EEFF', 'FFE4BE', 'FFA04C', 'FF7A26', 'FF6A19', 'FF500B', 'FF3403', 'FF2300');

echo '<h3>K -> RGB</h3>';
foreach ($temp as $k) {
    $rgb = ColourConverter::temperature2rgb($k);
    echo sprintf('<div style="background-color:rgb(%s); text-align: center; width: 100px; height: 25px; clear: both;">%s</div>', implode(', ', $rgb), $k);
}

echo '<h3>RGB -> K</h3>';
foreach ($hex as $v) {
    $rgb = array_values(ColourConverter::hex2rgb($v));
    $k = round(ColourConverter::rgb2temperature($rgb[0], $rgb[1], $rgb[2]));
    echo sprintf('<div style="background-color:rgb(%s); text-align: center; width: 100px; height: 25px; clear: both;">%s</div>', implode(', ', $rgb), $k);
}

Meine Ausgabe:

Ziemlich nah dran, aber nicht 100% noch. (Gefunden bug in meinem code, und es ist jetzt fast perfekt)

• Die Farben sind etwas abseits gehen von der k -> rgb -
• Funktioniert es nicht tun, k -> rgb -> k. Sie nicht wieder auf den gleichen Wert.

Code

UFRaw Linie 246-294

void Temperature_to_RGB(double T, double RGB[3])
{
    int c;
    double xD, yD, X, Y, Z, max;
    //Fit for CIE Daylight illuminant
    if (T <= 4000) {
        xD = 0.27475e9 / (T * T * T) - 0.98598e6 / (T * T) + 1.17444e3 / T + 0.145986;
    } else if (T <= 7000) {
        xD = -4.6070e9 / (T * T * T) + 2.9678e6 / (T * T) + 0.09911e3 / T + 0.244063;
    } else {
        xD = -2.0064e9 / (T * T * T) + 1.9018e6 / (T * T) + 0.24748e3 / T + 0.237040;
    }
    yD = -3 * xD * xD + 2.87 * xD - 0.275;

    //Fit for Blackbody using CIE standard observer function at 2 degrees
    //xD = -1.8596e9/(T*T*T) + 1.37686e6/(T*T) + 0.360496e3/T + 0.232632;
    //yD = -2.6046*xD*xD + 2.6106*xD - 0.239156;

    //Fit for Blackbody using CIE standard observer function at 10 degrees
    //xD = -1.98883e9/(T*T*T) + 1.45155e6/(T*T) + 0.364774e3/T + 0.231136;
    //yD = -2.35563*xD*xD + 2.39688*xD - 0.196035;

    X = xD / yD;
    Y = 1;
    Z = (1 - xD - yD) / yD;
    max = 0;
    for (c = 0; c < 3; c++) {
        RGB[c] = X * XYZ_to_RGB[0][c] + Y * XYZ_to_RGB[1][c] + Z * XYZ_to_RGB[2][c];
        if (RGB[c] > max) max = RGB[c];
    }
    for (c = 0; c < 3; c++) RGB[c] = RGB[c] / max;
}

void RGB_to_Temperature(double RGB[3], double *T, double *Green)
{
    double Tmax, Tmin, testRGB[3];
    Tmin = 2000;
    Tmax = 23000;
    for (*T = (Tmax + Tmin) / 2; Tmax - Tmin > 0.1; *T = (Tmax + Tmin) / 2) {
        Temperature_to_RGB(*T, testRGB);
        if (testRGB[2] / testRGB[0] > RGB[2] / RGB[0])
            Tmax = *T;
        else
            Tmin = *T;
    }
    *Green = (testRGB[1] / testRGB[0]) / (RGB[1] / RGB[0]);
    if (*Green < 0.2) *Green = 0.2;
    if (*Green > 2.5) *Green = 2.5;
}

Ich weiß, das ist eine alte Frage, aber ich war zu kämpfen, um eine Antwort zu finden, wie gut. Ich kam schließlich über diese nette Rechner, dass auch Beiträge, die die Formel wie folgt:

n = (x-0.3320)/(0.1858-y)

CCT = 437*n^3 + 3601*n^2 + 6861*n + 5517

Hoffe, dies hilft jemand, wer noch auf der Suche.