# Lieu planckien

En physique et en colorimétrie, le lieu plankien (ou lieu de Planck, ou encore lieu du corps noir) est la trajectoire (ou lieu géométrique) que suit la couleur d'un corps noir incandescent dans un diagramme de chromaticité donné. Cette trajectoire part du rouge profond, passe par l’orange, le jaune, le blanc, et s'achève dans le bleu à très hautes températures.
Un espace de couleur est un espace tridimensionnel, c'est-à-dire qu'une couleur y est spécifiée par trois nombres (les coordonnées X, Y et Z de l'espace CIE RGB de 1931, par exemple ; ou tout autre système comme « teinte, saturation, luminance »), qui permettent de spécifier la couleur et la brillance d'un stimulus visuel donné. Le diagramme de chromaticité est une projection à deux dimensions de cette représentation, obtenue en éliminant la luminance. L'espace CIE RGB transforme ainsi par projection ses trois coordonnées X, Y et Z en deux coordonnées de chromaticité notées x et y, conduisant au diagramme de chromaticité de la figure ci-contre. C'est généralement dans un tel diagramme de chromaticité qu'est représenté le lieu de Planck.
Inversement, la projection d'un point du diagramme de chromaticité sur ce lieu de Planck détermine la température de couleur de ce point.

## Calcul

La représentation du lieu de Planck dans un espace de couleur s'appuie sur les fonctions colorimétriques qui lui sont associées.
Dans l'espace CIE XYZ, par exemple, les trois coordonnées définissant la couleur d'une émission lumineuse sont définies par :
    $X_{T}=\int _{0}^{\infty }X(\lambda )M(\lambda ,T)\,d\lambda$
    $Y_{T}=\int _{0}^{\infty }Y(\lambda )M(\lambda ,T)\,d\lambda$
    $Z_{T}=\int _{0}^{\infty }Z(\lambda )M(\lambda ,T)\,d\lambda$
où M(λ,T) est l'exitance de l'émission lumineuse d'un corps noir de température T pour la longueur d'onde λ ; et X(λ), Y(λ) et Z(λ) sont les fonctions colorimétriques (représentées sur le diagramme ci-contre) représentant la sensibilité suivant chacun de ses trois récepteurs de l'observateur de référence en fonction de la longueur d'onde λ.
Le lieu de Planck est alors déterminé en substituant dans l'équation précédente la formule qui donne le spectre de la luminance énergétique spectrale du corps noir, laquelle suit la loi de Planck :
    $M(\lambda ,T)={\frac {2\mathrm {h} c^{2}}{\lambda ^{5}}}{\frac {1}{\mathrm {e} ^{\mathrm {h} c/(\lambda \mathrm {k} T)}-1}},$
où
* ν est la fréquence du rayonnement du corps noir, en s$^{−1}$ ;
* λ = λ$_0$/n, avec λ$_0$ la longueur d'onde du rayonnement du corps noir dans le vide et n l'indice de réfraction du milieu, est la longueur d'onde du rayonnement du corps noir dans le milieu, en m ;
* h = 6,626 070 15 × 10$^{−34}$ J s est la constante de Planck ;
* c = c$_0$/n, avec c$_0$ = 299 792 458 m s$^{−1}$ la vitesse de la lumière dans le vide et n l'indice de réfraction du milieu, est la vitesse de propagation du rayonnement du corps noir dans le milieu ;
* k = 1,380 649 × 10$^{−23}$ J K$^{−1}$ est la constante de Boltzmann ;
* T est la température de la surface du corps noir, en K.
L'intégration de ces équations donne alors les trois coordonnées du lieu de Planck dans l'espace CIE XYZ.

## Représentation dans l'espace de chromaticité

Si ces coordonnées sont notées X$_T$, Y$_T$ et Z$_T$, les coordonnées dans l'espace de chromaticité s'en déduiront par les formules :
    $x_{T}={\frac {X_{T}}{X_{T}+Y_{T}+Z_{T}}}$
    $y_{T}={\frac {Y_{T}}{X_{T}+Y_{T}+Z_{T}}}$
Le facteur $2\mathrm {h} c^{2}$ de la formule, qui dimensionne la luminance énergétique de la source, s'élimine dans l'espace de chromaticité, qui ne dépend que de la couleur.

### Paramétrage en mired

Le lieu de Planck dans l'espace xy correspond à la courbe représentée dans le diagramme de chromaticité ci-dessus. Il est bien entendu possible de calculer point par point les coordonnées de cette courbe en utilisant les formules ci-dessus, mais il est plus rapide dans les applications pratiques d'utiliser une formule approximative. Pour une telle approximation, il est usuel d'utiliser comme paramètre le mired, fonction inverse de la température, parce que le déplacement du point représentatif est plus régulier suivant ce paramètre qu'il ne l'est suivant la température absolue.

### Approximation dans l'espace CIE XYZ

Kim et al. utilisent comme approximation une spline cubique$^,$ :
$x_{c}={\begin{cases}-0.2661239{\frac {10^{9}}{T^{3}}}-0.2343580{\frac {10^{6}}{T^{2}}}+0.8776956{\frac {10^{3}}{T}}+0.179910&1667{\text{K}}\leq T\leq 4000{\text{K}}\\-3.0258469{\frac {10^{9}}{T^{3}}}+2.1070379{\frac {10^{6}}{T^{2}}}+0.2226347{\frac {10^{3}}{T}}+0.240390&4000{\text{K}}\leq T\leq 25000{\text{K}}\end{cases}}$
$y_{c}={\begin{cases}-1.1063814x_{c}^{3}-1.34811020x_{c}^{2}+2.18555832x_{c}-0.20219683&1667{\text{K}}\leq T\leq 2222{\text{K}}\\-0.9549476x_{c}^{3}-1.37418593x_{c}^{2}+2.09137015x_{c}-0.16748867&2222{\text{K}}\leq T\leq 4000{\text{K}}\\+3.0817580x_{c}^{3}-5.87338670x_{c}^{2}+3.75112997x_{c}-0.37001483&4000{\text{K}}\leq T\leq 25000{\text{K}}\end{cases}}$

### Approximation dans l'espace CIE UVW

Le lieu planckien peut également être approximé dans l'espace CIE UVW par l'expression suivante
${\bar {u}}(T)={\frac {0.860117757+1.54118254\times 10^{-4}T+1.28641212\times 10^{-7}T^{2}}{1+8.42420235\times 10^{-4}T+7.08145163\times 10^{-7}T^{2}}}$
${\bar {v}}(T)={\frac {0.317398726+4.22806245\times 10^{-5}T+4.20481691\times 10^{-8}T^{2}}{1-2.89741816\times 10^{-5}T+1.61456053\times 10^{-7}T^{2}}}$
Cette approximation donne de bons résultats, avec des écarts dans l'intervalle $\left|u-{\bar {u}}\right|<8\times 10^{-5}$ et $\left|v-{\bar {v}}\right|<9\times 10^{-5}$, pour des températures dans l'intervalle $1000K<T<15000K$.