FR 3042058 A1 20170407 FR 1563122 A 20151222 1 PROCÉDÉ DE CLASSIFICATION D'UN PRODUIT PAR ANALYSE DE LA DISTRIBUTION DES COULEURS DU PRODUIT DESCRIPTION Domaine technique [01] La présente invention se rapporte au domaine de la classification d'un produit parmi une pluralité de produits en se basant sur l'analyse de distribution des couleurs. État de la technique [02] Il est connu de l'art antérieur qu'on peut connaître certaines caractéristiques d'un produit par analyse colorimétrique. [03] Les mesures colorimétriques sont utilisées dans de nombreux domaines. En particulier : - Le domaine de l'agroalimentaire, où la couleur est utilisée pour contrôler la qualité d'une production ou le degré de maturité des fruits et légumes par exemple. Ces contrôles sont effectués en illuminant les produits avec une source lumineuse calibrée et avec des instruments étalonnés (colorimètre). Le contrôle d'un fruit dans un verger se fait par contre à l'oeil. A noter que les vignerons utilisent un réfractomètre afin d'avoir une indication du degré potentiel à partir du jus d'un raisin écrasé traversé par la lumière du soleil ; - Le domaine de l'imprimerie et de la peinture, où la couleur est bien évidement très importante. Dans ce domaine il existe plusieurs niveaux de mesure de la couleur. Pour retrouver rapidement une couleur parmi un choix de peintures on utilise des nuanciers composés des différents coloris disponibles (nuancier RAL K5 par exemple). L'identification se fait généralement à l'oeil mais il existe maintenant des applications pour téléphone intelligent permettant de retrouver la couleur la plus proche dans le nuancier à partir d'une photographie (ColorBlindC). Évidemment, en fonction des 3042058 2 conditions de prises de vue et de l'observateur, l'identification de la couleur est plus ou moins précise. Pour la restauration de fresque ou de tableau, la mesure de la couleur doit être extrêmement précise afin de pouvoir doser les pigments. Celle-ci ne peut alors 5 plus se faire à l'oeil du fait notamment du phénomène de métamérisme qui fait que deux couleurs peuvent apparaitre similaires sous un éclairage donné mais différentes sous un autre éclairage. C'est particulièrement le cas pour les bleus. Il est alors nécessaire d'utiliser un spectromètre qui décompose la lumière 10 issue de l'objet, et donc la couleur. Les pigments sont alors choisis de manière à reproduire le spectre de la peinture et non sa couleur sous une lumière donnée. - Dans le domaine de la santé, on citera les nuanciers utilisés par les dermatologues pour définir les peaux à risques ou encore ceux 15 fournis avec certains dentifrices pour contrôler, à l'oeil, l'effet du produit sur la blancheur des dents. - Enfin, dans le domaine de la recherche scientifique, la colorimétrie se fait généralement à l'aide de méthodes spectroscopiques. Elle est très utilisée en biophysique pour l'étude des agents de 20 contrastes optiques et des fluorophores. D'autres applications utilisent la couleur à des fins de détermination comme la microscopie optique sur des préparations de roches ou de glaces par exemple. [4] Pour des besoins sophistiqués, par exemple l'analyse du 25 vieillissement d'un matériau par l'évolution de sa couleur, on procède généralement à un travail dans un environnement contrôlé. [5] En général, ce type d'analyse commence par un calibrage/étalonnage de la chaine de prise de vue de façon à obtenir un référentiel absolu en couleur. Pour cela, après avoir mis en place 30 l'éclairage qui sera utilisé dans toute la suite, une mire colorimétrique est photographiée et chaque appareil de la chaine de prise de vue (l'appareil 3042058 3 photographique, l'écran de visualisation, l'imprimante, etc.) est étalonné de sorte que les couleurs reproduites à chaque étape soient le plus proches possibles de celles de la mire colorimétrique. [6] Puis le produit à analyser/classifier est photographié dans les mêmes 5 conditions d'éclairage pour garantir que les couleurs acquises correspondent à celles du référentiel. [7] Ensuite, la distribution de couleurs du produit est extraite, éventuellement normalisée et sert de base à une classification / comparaison permettant de déduire certaines informations physico- 10 chimiques du produit. [8] Toutefois, ce processus est lourd dans sa mise en place et demande énormément d'attention dans son exécution. Par exemple, si une lampe d'éclairage défaille et doit être remplacée, toute l'opération de calibrage/étalonnage doit être recommencée. 15 [09] Pour éviter cela, il a été proposé de laisser la mire colorimétrique visible sur chaque photographie du produit à analyser/classifier de sorte à avoir sur chaque photographie une référence colorimétrique permettant de calibrer/étalonner la photographie avant de faire l'analyse du produit. Par exemple, la demande de brevet US 2012/288195 décrit une telle solution. 20 [10] Malgré cela, les procédés d'analyse colorimétrique actuellement utilisés présentent les inconvénients de nécessiter la connaissance de la couleur vraie de l'objet à caractériser. Par couleur vraie, l'on entend son positionnement dans un référentiel connu de l'état de l'art. Cet état de fait implique l'utilisation d'une mire colorimétrique standard constituée de 25 références couleurs positionnées dans le référentiel connu suivant leurs valeurs théoriques. Cette recherche de la valeur vraie des couleurs limite énormément le procédé de classification : Premièrement par la limitation des valeurs prises par les couleurs de références, deuxièmement, par la contrainte de rester dans un référentiel couleur connu. 30 [11] Il existe donc un réel besoin de mise au point d'un procédé de classification d'un produit palliant ces défauts, inconvénients et obstacles 3042058 4 de l'art antérieur, en particulier d'un procédé permettant de s'affranchir de la recherche de couleur vraie, réduire les coûts et encombrements, et d'améliorer la procédure de classification en libérant la contrainte de repérage dans un référentiel standard. 5 Description de l'invention [12] Pour résoudre un ou plusieurs des inconvénients cités précédemment, un procédé de classification d'un produit à classifier comprend : - un sous-procédé de création d'une distribution des couleurs du 10 produit à classifier comprenant : - prise d'une photo numérique du produit à classifier et de la référence mire dans des conditions d'éclairement identiques ; - extraction, à partir de la photo numérique, d'une distribution des couleurs du produit à classifier et d'une distribution des couleurs 15 de la référence mire ; puis - classification du produit à classifier par rapport à une pluralité de produits prédéterminés, chaque produit prédéterminé étant associé à une distribution des couleurs et à une distribution des couleurs de la référence mire, ces deux distributions ayant été obtenues par 20 utilisation du sous-procédé de création de distribution des couleurs, une mesure de similarité étant définie entre les distributions de couleurs du produit à classifier et chacun des produits prédéterminés, ladite mesure de similarité étant calculée après une calibration des distributions de couleurs des produits par une 25 fonction de calibrage minimisant la différence entre les distributions de la référence mire associées auxdits produits. [13] Des caractéristiques ou des modes de réalisation particuliers, utilisables seuls ou en combinaison, sont : - la référence mire est installée près du produit à classifier et une 30 photo numérique unique est prise des deux objets ensemble ; 3042058 5 - une distribution des couleurs étant définie comme une répartition du nombre de pixels associés à chaque couleur, avant la classification, une normalisation du nombre de pixels des distributions est réalisée ; 5 - les distributions des couleurs des produits sont transformées par une fonction de transformation avant ou après la calibration ; - une distribution des couleurs étant définie comme une pluralité de distributions de couleurs élémentaires dans un référentiel prédéterminé, notamment un référentiel RGB, la fonction de 10 transformation transforme la distribution des couleurs en un vecteur unique de distributions des couleurs élémentaires ; - la mesure de similarité est obtenue par l'application d'un algorithme de classification à l'ensemble constitué de la pluralité de produits prédéterminés plus le produit à classifier, notamment un algorithme 15 de minimisation des distances entre produits dans un référentiel prédéterminé. - la référence mire est un objet dont l'image présente une distribution de couleurs recouvrant la palette de couleurs des produits prédéterminés. 20 [14] Dans un second aspect de l'invention, le procédé précédent est utilisés lorsque les produits sont des roches broyées, les produits prédéterminés étant des roches dont des caractéristiques physico-chimiques ont été mesurées préalablement, de telle sorte que la classification du produit à classifier détermine ses caractéristiques physico- 25 chimiques [15] Un autre aspect de l'invention se rapporte à un programme d'ordinateur à installer dans un ordinateur, comprenant des instructions pour mettre en oeuvre les étapes d'un procédé de classification tel que défini ci-dessus lors d'une exécution du programme par une unité de calcul 30 dudit appareil. 3042058 6 Brève description des figures [16] L'invention sera mieux comprise à la lecture de la description qui suit, faite uniquement à titre d'exemple, et en référence aux figures en annexe dans lesquelles : 5 - La figure 1 représente un histogramme du nombre de pixels en fonction de la valeur du canal rouge de la photo ; - La figure 2 représente un ordinogramme du procédé selon un mode de réalisation de l'invention ; - La figure 3 représente une mire colorimétrique de calibration ; 10 - La figure 4 représente la structure d'un système argentique (a), d'un détecteur numérique CCD contenant une unique couche de pixels agencés selon la matrice de Bayer (b) et du détecteur numérique Foveon X3 (c) ; - La figure 5 représente les valeurs RGB de l'échelle de gris 15 théorique d'une mire en fonction des valeurs mesurées et l'approximation par une courbe de tendance de forme polynomiale d'ordre 2 ; [17] La figure 6 représente les histogrammes en boite permettant de comparer les différentes poudres dans le canal rouge, vert et bleu après correction par un polynôme d'ordre 5 et, sur la droite, le pourcentage 2 0 d'identification pour les 21 poudres ; [18] La figure 7 représente le vecteur de la poudre 14FR27 avec les histogrammes du canal R (a), du canal G (b) et du canal B (c) ; - La figure 8 représente les 5 premières lignes et dernière ligne de deux tableaux de correspondance pour une poudre non correctement 25 identifiée (14FR24) et une poudre correctement identifiée (14FR07) ; et - La figure 9 représente les diagrammes de taux d'identification en fonction des poudres de la base de données pour les poudres 14FR22 (a), 14FR08 (b), 14FR12 (c). Modes de réalisation 30 [19] En préliminaire, on rappellera quelques notions utiles à la compréhension de l'invention. 3042058 7 [20] Le système rouge, vert, bleu (RGB pour « Red, Green, Blue » en anglais) est le système d'encodage utilisé par tous les systèmes d'imagerie optiques dans le visible, depuis la pellicule argentique à l'oeil humain, en passant par les capteurs CCD ou les téléviseurs. À partir de la synthèse 5 additive des trois couleurs primaires il est possible de restituer l'ensemble des couleurs de chromaticité moindre. En informatique, chaque couleur est généralement encodée en 8 bits (256 niveaux de 0 à 255) ce qui permet d'avoir 2563 couleurs différentes, soit environ 16,78 millions de couleurs. [21] Une image couleur peut être représentée par une matrice de 10 dimension MxNx3 où M et N sont respectivement la longueur et la largeur de l'image en pixels et 3 représente le nombre de canaux : un rouge, un vert et un bleu. Les pixels sont donc repérés par leurs coordonnées (x,y,Z) ou Z = R, G ou B. La dynamique de l'image est représentée par trois histogrammes correspondant aux distributions des valeurs de R, G et B sur 15 l'ensemble des pixels. Le maximum d'une distribution est ainsi la valeur de R, G ou B prise par le plus grand nombre de pixels dans l'image. On peut voir par exemple sur la figure 1 que le maximum pour le canal rouge est de 88. [22] On définit le paramètre de luminance d'une couleur en RGB par : 20 p=(R+G+B)/3 (2.1) [23] Si le système RGB est le système le plus couramment utilisé, il est parfois utile de travailler avec d'autres formats d'encodage des couleurs obtenus à partir des valeurs RGB. Les principaux sont rgb, YCbCr, tsv, Yxy et L*a*b*. Ces différents formats sont bien connus de l'homme du métier. 25 [24] Dans la suite de cette description, nous utiliserons le format RGB comme référentiel. [25] Par ailleurs, nous prendrons à titre d'exemple d'utilisation du procédé la classification des roches volcaniques. [26] En référence à la figure 2, un procédé de classification d'un produit 30 comprend un sous-procédé 10 de création d'une distribution des couleurs du produit à classifier. 3042058 8 [27] Le produit est photographié, étape 12, ainsi qu'une mire colorimétrique, étape 14, dans des conditions d'éclairement identiques. [28] Il est possible de réaliser deux photographies successives : le produit puis la mire mais, pour assurer les conditions d'éclairement identiques, il 5 est souvent préférable de poser la mire à côté du produit et de prendre une photographie unique des deux éléments côte à côte. [29] Dans l'application des roches volcaniques, la roche à étudier est d'abord broyée puis étalée dans une coupelle. Une mire colorimétrique est posée à côté de la coupelle et une photographie est prise. 10 [30] La mire est, par exemple, une mire colorimétrique du marché, Figure 3 comme la mire de calibration ColorCheckee composée de carrés de différentes couleurs. Les valeurs théoriques Rth, Gth et i'-sth de chaque carré dans le système sRGB GMB (standard RGB GretagMacbeth) sont reportées sur la figure. On note en particulier la ligne du bas dont les 15 carrés représentent l'échelle de gris correspondant à une équiproportion des coefficients chromatiques. Pour des raisons pratique, ces carrés sont numérotés de 1 à 6 du blanc vers le noir, et leurs valeurs théoriques notés Rthi, Gthi et avec i = 1 à 6. [31] La mire peut être également une mire spécialement conçue pour 20 l'application et les produits à classer. Dans ce cas, la mire contiendra une gamme de couleurs recouvrant la palette de couleurs des produits. [32] Puis, à partir de la(es) photo(s), on extrait, étape 16, une distribution des couleurs du produit à classifier et une distribution des couleurs de la référence mire. 25 [33] Par exemple, comme expliqué ci-dessus, on extrait les histogrammes des canaux RGB. [34] Dans une étape optionnelle 18, les histogrammes sont normalisées c'est-à-dire que par application d'un ratio adapté, ils prennent en compte un nombre de pixels prédéterminés : par exemple, 10 000 pixels. Ou bien 30 le point haut de chaque histogramme est positionné sur une même valeur, par exemple 100. 3042058 9 [35] Le produit photographié est alors classifier, étape 20 par rapport à une pluralité de produits prédéterminés. Chaque produit prédéterminé est associé à une distribution des couleurs et à une distribution des couleurs de la référence mire, ces deux distributions ayant été obtenues par 5 utilisation du sous-procédé de création de distribution des couleurs décrit ci-dessus. [36] Dans un mode de réalisation particulier, chaque produit est mémorisé sous forme d'un enregistrement d'une base de données, chaque enregistrement comprenant les informations suivantes : 10 - Un identifiant qui peut être un nom ou une référence ; - Une distribution des couleurs du produit ; - Une distribution des couleurs de la mire référence. Ces deux distributions ont été obtenues via une prise de vue sous condition d'éclairement identique comme indiqué à propos des étapes 12, 15 14 ; - Des caractéristiques physico-chimiques associées au produit. [37] Une mesure de similarité est définie entre les distributions de couleurs du produit à classifier et chacun des produits prédéterminés. Cette mesure de similarité est calculée 22 après une calibration 24 des 2 0 distributions de couleurs des produits par une fonction de calibrage minimisant la différence entre les distributions de la référence mire associées auxdits produits. [38] Ainsi, par exemple, le traitement est le suivant, pour chaque couple produit à classifier / produit prédéterminé : 25 - Transformation des distributions, étape 26. Cette étape optionnelle vise à appliquer une transformation mathématique des distributions permettant de faciliter les traitements suivants. On notera que cette transformation mathématique est la même pour toutes les distributions de couleurs. Par ailleurs, cette transformation peut 30 s'appliquer au moment du calcul de similarité ou préalablement, par exemple après l'étape de normalisation. Un exemple de 3042058 10 transformation est donné dans l'exemple détaillé ci-dessous : les distributions des couleurs RVB sont transformées en un vecteur unique de distributions des couleurs élémentaires. Il est alors avantageux que les enregistrements de la base de données 5 contiennent les distributions transformées. Pour la suite, le terme « distribution » sera utilisé indifféremment pour nommer la distribution ou sa version transformée. - Calibration des distributions, étape 24. Cette calibration consiste à comparer les distributions de la mire référence associées au produit 10 à classifier et au produit de la base de données et à rechercher une fonction mathématique qui minimise, si ce n'est annule, les différences entre ces deux distributions, puis à appliquer cette fonction mathématique sur la distribution du produit à classifier, par exemple, de façon à ce que les deux distributions des produits 15 partagent le même référentiel. Comme décrit ci-dessous, une telle fonction mathématique peut être une fonction polynomiale passant par les points de mesure. Il est à noter que cette étape de calibration est une calibration relative : il n'y a pas de recherche d'une couleur « vraie » comme dans l'art antérieur mais seulement 20 d'un partage de référentiel relatif d'un produit par rapport à un autre. Il est envisageable de prendre un produit, par exemple le premier photographié, et de le labelliser « étalon » et alors de calibrer tous les autres produits par rapport à ce produit « étalon ». Mais même, dans ce mode de réalisation particulier, il n'y a pas de recherche de 25 « vraies » couleurs. - Calcul de la similarité, étape 22. On définit donc une mesure de similarité, ou une distance, dans l'espace des distributions des produits. On notera qu'il est habituel de définir une mesure de similarité comme une fonction qui croit en même temps que la 30 similarité augmente et de définir une distance comme une fonction qui diminue en même temps que les produits se rapprochent. On 3042058 11 utilise l'un ou l'autre terme indifféremment pour indiquer une fonction qui varie de façon monotone avec la similarité des distributions. [39] Le produit à classifier est alors classé, étape 30, par une recherche 5 du minimum de la distance, ou du maximum de la mesure de similarité, dans le référentiel prédéterminé. [40] Dans une étape optionnelle 35 les caractéristiques physico-chimiques du produit à classifier sont déduites des caractéristiques physico-chimiques du/des produit(s) le(s) plus similaire(s). 10 [41] On conçoit que la mise en oeuvre du procédé est réalisé par un ordinateur convenablement programmé par un produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur. Exemple 15 [42] Un exemple détaillé de mise en oeuvre du procédé ci-dessus va maintenant être décrit. [43] Il se rapporte à la classification des roches volcaniques. [44] Un appareil photo est équipé d'un détecteur Foveon X3 développé et fourni par la compagnie Sigmê. Ce détecteur a une architecture originale 20 qui lui permet de se rapprocher au mieux du principe des appareils photos argentiques (Figure 4). En effet, il est constitué de trois couches de pixels à savoir, une pour le rouge, une pour le vert et enfin une pour le bleu (Fig. 4 c). [45] Dans un premier temps, les images ont été réalisées avec un 25 appareil photographique numérique Sigma SD15 commercial. Cependant, lors de la prise de vue, ce système applique différentes corrections au signal. Parmi celles-ci, les principales sont : - La linéarisation Look Up Table (LUT) : la fonction représentant le signal d'entrée en fonction du signal de sortie n'est pas linéaire. Le 30 but de la linéarisation LUT est donc de remplacer cette fonction par une application affine ; 3042058 12 - Le Dark Frame Substract : il permet de retirer du bruit à une image linéarisée en soustrayant à cette dernière une image du noir (« dark ») elle-même linéarisée ; - Des corrections sont appliquées pour retirer du bruit ou de 5 «mauvais pixels» ceux-ci étant remplacés par la valeur moyenne de leurs voisins. Ces corrections peuvent poser problème dans le cadre du procédé car les « mauvais pixels » peuvent être en fait dû à une particularité de la poudre, - Enfin, différentes opérations matricielles sont ensuite effectuées 10 afin de modifier les couleurs sur chacun des trois canaux dans le but d'obtenir des couleurs proches de celles observées par l'oeil humain et de compenser une luminosité trop faible ou trop forte. [46] Ces algorithmes sont particulièrement gênants puisqu'ils vont induire des changements de couleurs difficiles à contrôler. Il est donc préférable 15 de travailler sur le signal brut en sortie du détecteur. Pour ce faire, nous avons utilisé un boîtier équipé d'un détecteur Foveon Xe. Il est fixé sur un statif vertical permettant de prendre les photos par le dessus. Des lampes réglables sont fixées de part et d'autre du champ d'observation. Le boîtier est équipé d'un objectif à focale fixe dont l'ouverture du diaphragme est 20 réglée manuellement. Enfin, il est relié à un ordinateur muni du logiciel EPIX XCAP permettant de traiter les informations brutes envoyées par le détecteur et de contrôler les paramètres d'acquisition dont le temps d'exposition (c'est-à-dire le temps durant lequel le détecteur compte les photons qui lui arrive), et le format d'enregistrement. Les images peuvent 25 être enregistrées au format texte ou au format tif. L'ensemble des données est alors traité avec le logiciel Matlab 2010. [47] Pour cet exemple, 21 poudres de roches volcaniques sont utilisées. Les poudres sont non tamisées mais les grains supérieurs à 2mm sont retirés. Les différentes poudres sont posées en tas de forme aléatoire à 30 côté de la mire. Différentes photos sont prises des différents tas de 3042058 13 poudres dans différentes conditions d'éclairage, avec différents temps d'exposition et différentes ouvertures du diaphragme de l'objectif. [48] Le format utilisé ne permettant pas de faire apparaître des paramètres particulièrement pertinents, le travail est réalisé en RGB, ce 5 format étant représentatif de la physique du détecteur. [49] L'objectif premier est d'obtenir une méthode d'étalonnage permettant d'obtenir une même image à partir de deux images d'un même objet réalisées dans des conditions de prises de vue différentes. On définit le paramètre de qualité d'étalonnage, Qe : 111-121 (3.A.1) 10 Qe = 100 x 100 12 [50] où 11 et 12 sont deux images couleurs. [51] Sur la mire, les valeurs minimales et maximales de R, G et B sont associées respectivement au carré noir 6 et au carré blanc 1 (Fig. 3). Il apparait pertinent de ré-étaler la dynamique de la mire en replaçant les 15 valeurs moyennes de ces carrés à leur valeurs théoriques en utilisant les formules suivantes : [52] RL th th R1 -R6 x (R - R6) + R6h (3.A.5) [53] GL R1 -R6 x (G - G6) + Ge (3.A.6) [54] BL th th G1 -G6 x (B - B6) + B6h (3.A.7) G1 -G6 th th B1 -B6 B1 -B6 20 [55] où : [56] R1 (resp. G1 et B1) est la valeur moyenne mesurée du carré blanc dans le canal R (resp. G et B), [57] R6 (resp. G6 et B6) est la valeur moyenne mesurée du carré noir dans le canal R (resp. G et B). 25 [58] Avec cette méthode, le facteur Qe passe de Qe = 85,71% avant étalonnage à Qe = 95,27% après étalonnage. [59] Les couleurs des poudres formant un nuancier de gris, l'échelle de gris de la mire est utilisée pour étalonner les images. La figure 5 montre les valeurs théoriques des carrés de l'échelle de gris en fonction des valeurs 3042058 14 moyennes mesurées sur une photo réalisée avec un temps d'exposition de 110 ms et une ouverture de 5,6. [60] On observe que sur aucun canal la réponse du détecteur n'est linéaire. Pour utiliser tous les carrés de l'échelle de gris, le polynôme 5 d'ordre 5 permettant de relier les valeurs moyennes mesurées pour chaque carré aux valeurs théoriques de la mire est calculé. [61] L'obtention des polynômes ce fait en résolvant les 3 systèmes d'équations suivants : (3.A.8) aC b -4 (.[C 1 +-eCL-1- C'," c1(.7+cC. aC:1)C. +cC, eCr+f=C,'" aC.;-bC +eC. bC; +eC_+f- [62] où C1,..., C6 représente les valeurs moyennes mesurées des carrés de l'échelle de gris, C = R, G ou B. La résolution de ce système a permis d'obtenir les coefficients a, b, c, d, e 20 et f b--45 -a° sous la forme : [63] (3.A.9) 25 d= 0 0-b e= )- -c e)-( f= t C2-eC, -aCT1-b( [64] avec : 30 - c(= - P=C15-c52 - y=c14_c24 - EI=C13-c32 - e=c12_c22 - =c1-c2 5 - 71 = Cr - Cr ce(C2 -C3) - 61 = cl - ci ie(C2-C3) - t = Cl - cl- Y(C2-C3) - K =0-c, 8(C -C3) - À = £(C2-C3) 3042058 15 10 - = - n(q -Cn ^ (ne a 4. (C3 C4) - V = - 0(q-C4)+ /Os - = -C4 fl) - C4) - 0 = - -C4) ( te y) (c 3 c 4) -C4)+ ke 8\ ) (C3 C ) REVENDICATIONS1. Procédé de classification d'un produit à classifier comprenant : - un sous-procédé de création (10) d'une distribution des couleurs du produit à classifier comprenant : « prise d'une photo numérique (12, 14) du produit à classifier et de la référence mire dans des conditions d'éclairement identiques ; - extraction (16), à partir de la photo numérique, d'une distribution des couleurs du produit à classifier et d'une distribution des couleurs de la référence mire ; puis - classification (20) du produit à classifier par rapport à une pluralité de produits prédéterminés, chaque produit prédéterminé étant associé à une distribution des couleurs et à une distribution des couleurs de la référence mire, ces deux distributions ayant été obtenues par utilisation du sous-procédé de création de distribution des couleurs, une mesure de similarité étant définie (22) entre les distributions de couleurs du produit à classifier et chacun des produits prédéterminés, ladite mesure de similarité étant calculée après une calibration (24) des distributions de couleurs des produits par une fonction de calibrage minimisant la différence entre les distributions de la référence mire associées auxdits produits. 2. Procédé selon la revendication 1, caractérisé en ce que la référence mire est installée près du produit à classifier et une photo numérique unique est prise des deux objets ensemble. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que, une distribution des couleurs étant définie comme une répartition du nombre de pixels associés à chaque couleur, avant la classification, une normalisation du nombre de pixels des distributions est réalisée. 3042058 21 4. Procédé selon la revendication 1, 2 ou 3, caractérisé en ce que les distributions des couleurs des produits sont transformées (26) par une fonction de transformation avant ou après la calibration. 5 5. Procédé selon la revendication 4, caractérisé en ce que, une distribution des couleurs étant définie comme une pluralité de distributions de couleurs élémentaires dans un référentiel prédéterminé, notamment un référenciel RGB, la fonction de transformation transforme la distribution des couleurs en un vecteur 10 unique de distributions des couleurs élémentaires. 6. Procédé selon la revendication 5, caractérisé en ce que la mesure de similarité est obtenue par l'application d'un algorithme de classification à l'ensemble constitué de la pluralité de produits 15 prédéterminés plus le produit à classifier, notamment un algorithme de minimisation des distances entre produits dans un référentiel prédéterminé. 7. Procédé selon l'une quelconque des revendications précédentes, 20 caractérisé en ce que la référence mire est un objet dont l'image présente une distribution de couleurs recouvrant la palette de couleurs des produits prédéterminés. 8. Utilisation du procédé selon l'une quelconque des revendications 25 précédentes, caractérisé en ce que les produits sont des roches broyées, les produits prédéterminés étant des roches dont des caractéristiques physico-chimiques ont été mesurées préalablement, de telle sorte que la classification du produit à classifier détermine ses caractéristiques physico-chimiques. 30 3042058 22 9. Produit programme d'ordinateur comprenant des portions de code de programme téléchargeables depuis un réseau de communication et/ou enregistré sur un support lisible par ordinateur et/ou exécutable par un processeur, caractérisé en ce qu'il comprend des instructions 5 de code de programme pour la mise en oeuvre du procédé selon l'une au moins des revendications 1 à 8 lorsque ledit programme fonctionne sur un ordinateur.