FR 3072806 A1 20190426 FR 1759814 A 20171019 Procédé de calcul d'un descripteur global d'une image La présente invention concerne le domaine de l'extraction et de l'analyse de motifs ou caractéristiques visuelles dans une image, une partie d'une image ou une séquence de plusieurs images. L'invention concerne notamment le domaine de la vision par ordinateur et porte plus précisément sur un procédé de calcul d'un descripteur global d'une image. Un descripteur global d'une image a pour fonction de caractériser l'information contenue dans l'image, par exemple les propriétés générales de 10 la scène de l'image. Il peut être constitué de descripteurs locaux correspondants chacun à un pixel de l'image parmi tous les pixels de l'image, d'une partie de l'image ou d'un ensemble de points d'intérêts de l'image. L'élaboration de descripteurs d'image est couramment utilisée dans le domaine de la vision par ordinateur, notamment pour des applications de 15 reconnaissance d'objets (de visages par exemple), de classification ou d'analyse de contenu. L'invention est applicable à tout type d'image numérique et en particulier aux images multi-spectrales dont chaque pixel comporte plusieurs composantes à différentes longueurs d'onde. L'invention est aussi applicable 20 à des images à trois composantes dites RGB. De nombreuses méthodes de calcul de descripteurs sont disponibles dans la littérature. On peut citer notamment le descripteur SURF « Speeded Up Robust Features » décrit dans la demande de brevet européenne EP 25 1850270, le descripteur SIFT « Scale Invariant Feature Transform » décrit dans le brevet américain U5671 1293 ou encore le descripteur LBP « Local Binary Patterns » décrit notamment dans la référence [1]. La figure 1 schématise le principe de la méthode de calcul d'un descripteur LBP selon l'état de l'art. Un descripteur local est calculé pour 30 chaque pixel de l'image ou d'une partie de l'image. La première étape de la méthode consiste à définir un voisinage V autour du pixel, par exemple un 2 bloc de 3x3 pixels autour d'un pixel central P. Les huit pixels entourant le pixel central Po forment une boucle qui est parcourue dans un sens choisi par convention. Chaque pixel de la boucle est comparé avec le pixel central et on lui affecte respectivement le signe « = », «-» ou « + » selon que sa 5 valeur est égale, inférieure ou supérieure à celle du pixel central. On encode ensuite ce résultat en affectant au signe «-» la valeur 0 et aux signes « = » et « + » la valeur 1. On multiplie ensuite ce résultat par une matrice de pondération M qui comporte, pour chaque pixel parcouru dans l'ordre initialement choisi, une puissance de deux incrémentée à chaque nouveau 10 pixel. Le résultat de cette multiplication donne le descripteur local du pixel central Po qui vaut 241 dans l'exemple de la figure 1. Le calcul de descripteur LBP est applicable à des images monochromatiques, c'est-à-dire pour lesquelles un pixel n'a qu'une seule 15 valeur. Cette méthode n'est pas adaptée à des images en couleurs ni à des images multi-spectrales comprenant plusieurs composantes pour différentes longueurs d'onde. De façon générale, la méthode de calcul de descripteur LBP souffre d'un manque de robustesse car de faibles changements d'intensité dans une 20 image peuvent engendrer des changements importants dans les motifs extraits à l'aide de cette méthode. En effet, cette méthode ne prend pas en compte les changements d'intensité des pixels d'un voisinage. Si tous les pixels d'un voisinage ont une intensité diminuée d'un même facteur, le descripteur local calculé pour ce voisinage sera toujours le même. 25 L'invention propose un procédé de calcul d'un descripteur global d'une image ou d'une partie d'image qui permet d'obtenir un descripteur spatial et spectral conjointement à la différence de la méthode LBP qui ne peut s'appliquer qu'à une seule longueur d'onde à la fois, de façon indépendante. 30 L'invention est applicable à des images multispectrales, notamment des images RGB, ou des images hyperspectrales qui comportent plusieurs 3 longueurs d'ondes allant de trois canaux à plusieurs centaines. Elle permet de calculer un seul descripteur pour tout le domaine spectral disponible de l'image et se démarque des méthodes de l'art antérieur qui appliquent en général un traitement séparé et indépendant pour chaque longueur d'onde. 5 Par ailleurs, l'invention permet de générer une base de descripteurs exhaustive qui comprend tous les motifs qui peuvent être extraits d'une image. L'invention a pour objet un procédé, mis en oeuvre par ordinateur, de calcul d'un descripteur global d'au moins une partie d'une image à décrire, le descripteur global comprenant un descripteur local associé à chaque pixel de l'au moins une partie d'image à décrire, le procédé comprenant le calcul de tous les descripteurs locaux d'un descripteur global en : - Définissant un voisinage du pixel associé au descripteur local, - Attribuant une étiquette à chaque pixel du voisinage, deux pixels de valeurs différentes ayant deux étiquettes différentes, deux pixels de valeurs identiques ayant la même étiquette, l'ensemble des étiquettes du voisinage formant un motif, - Affectant la valeur du descripteur local à un représentant d'une classe d'équivalence à laquelle appartient le motif. Selon un mode de réalisation particulier, le procédé selon l'invention comprend au préalable les étapes de : - Obtenir une première image acquise avec un premier temps d'exposition, - Simuler au moins une seconde image résultant d'une prise de vue de la première image avec au moins un second temps d'exposition inférieur au premier temps d'exposition, - Calculer un descripteur global pour au moins une partie d'une image à décrire prise parmi la première image acquise ou l'au moins une seconde image simulée. 4 Selon un aspect particulier de l'invention, l'étape de simuler une seconde image est réalisée en diminuant l'intensité des pixels de la première image d'un facteur égal au rapport entre le premier temps d'exposition et le second temps d'exposition. 5 Selon un aspect particulier de l'invention, l'au moins une partie de la première image acquise est composée de pixels dont les valeurs sont toutes différentes entre elles. Selon un mode de réalisation particulier, le procédé selon l'invention comprend une étape de génération de la première image acquise à partir 10 d'une image multi-spectrale en sélectionnant le nombre minimal de composantes spectrales de l'image nécessaire pour obtenir que l'au moins une partie de la première image acquise est composée de pixels dont les valeurs sont toutes différentes entre elles. Selon un mode de réalisation particulier, le procédé selon l'invention 15 comprend les étapes de faire varier le second temps d'exposition entre le premier temps d'exposition et zéro et de calculer un descripteur global pour au moins une partie de la première image acquise et de l'au moins une seconde image simulée. Selon un mode de réalisation particulier, le procédé selon l'invention 20 comprend une étape d'affichage des descripteurs globaux calculés en fonction du second temps d'exposition. Selon un mode de réalisation particulier, le procédé selon l'invention comprend un calcul d'entropie des descripteurs locaux pour chaque descripteur global calculé et la sélection d'au moins un descripteur global 25 d'au moins une image à décrire en fonction des entropies calculées. Selon un mode de réalisation particulier, le procédé selon l'invention comprend la sélection du descripteur global dont l'entropie calculée est la plus élevée. Selon un aspect particulier de l'invention, la classe d'équivalence est 30 une classe de chaines d'étiquettes invariantes par rotation. 5 Selon un aspect particulier de l'invention, la classe d'équivalence est une classe de chaines d'étiquettes qui est en outre invariante par permutation de valeurs d'étiquettes. Selon un aspect particulier de l'invention, la classe d'équivalence est 5 une classe de chaines d'étiquettes qui est en outre invariante par symétrie. Selon un aspect particulier de l'invention, une image à décrire est une image multi-spectrale dont chaque pixel comporte plusieurs composantes à différentes longueurs d'ondes. L'invention a encore pour objet un programme d'ordinateur 10 comportant des instructions pour l'exécution du procédé de calcul d'un descripteur global selon l'invention, lorsque le programme est exécuté par un processeur et un support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution du procédé de calcul d'un descripteur global selon l'invention, 15 lorsque le programme est exécuté par un processeur. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés qui représentent : 20 - La figure 1, un schéma illustrant un exemple de calcul de descripteur local selon une méthode de l'art antérieur, - La figure 2, un organigramme décrivant les étapes de mise en oeuvre d'un procédé selon un premier mode de réalisation de l'invention, - La figure 3, un schéma illustrant un exemple de calcul de descripteur 25 selon le premier mode de réalisation de l'invention, - La figure 4, un schéma représentant tous les colliers non-étiquetés de longueur égale à 4 et dans un alphabet à 4 éléments, - La figure 5, un organigramme décrivant les étapes de mise en oeuvre d'un procédé selon un deuxième mode de réalisation de l'invention, 30 - La figure 6, un exemple d'histogramme à trois dimensions produit par le procédé selon le deuxième mode de réalisation de l'invention, 6 - La figure 7, un exemple de descripteurs globaux représentés sous forme de séquence d'images pour différents temps d'exposition, - La figure 8, un organigramme décrivant les étapes de mise en oeuvre d'un procédé selon un troisième mode de réalisation de l'invention, 5 - La figure 9 représente, sur trois diagrammes, un exemple de résultat du procédé décrit à la figure 8, - La figure 10, plusieurs représentations d'une même image avec un calcul de descripteurs selon l'art antérieur d'une part et un calcul de descripteurs selon le deuxième mode de réalisation de l'invention 10 d'autre part. Les figures 2 et 3 illustrent, respectivement sur un organigramme et un schéma, les étapes de mise en oeuvre d'un procédé de calcul d'un descripteur global d'une image I ou partie d'image selon un premier mode de 15 réalisation de l'invention. Un descripteur global est composé d'une pluralité de descripteurs locaux associés chacun à un pixel de l'image I ou d'une zone de l'image que l'on souhaite décrire. La figure 3 illustre le calcul d'un descripteur local pour un pixel P. L'image I considérée peut être une image multi-spectrale, hyper- 20 spectrale ou une image à une seule composante spectrale (image niveau de gris). De façon générale, chaque pixel de l'image est représenté par un vecteur à N composantes spectrales correspondants à N longueurs d'ondes différentes, N étant un entier au moins égal à un. Dans une première étape 201 on sélectionne un voisinage V du pixel 25 P. Un voisinage est constitué de plusieurs pixels entourant le pixel P. Le pixel P peut être compris dans le voisinage ou en être exclu. Sur la figure 3, on a représenté un exemple particulier de voisinage V qui est égal à une boucle de 8 pixels entourant le pixel P. Alternativement, le voisinage V peut être égal à un motif carré de 2 pixels par 2 pixels, comprenant le pixel P, ou 30 encore à un motif rectangle de 3 pixels par 2 pixels, comprenant le pixel P ou 7 tout autre motif de pixels voisin du pixel P. Sur l'exemple de la figure 3, le voisinage V est composé des pixels Po à P7. Dans une deuxième étape 202 du procédé, on détermine une étiquette ou label pour chaque pixel du voisinage V en fonction de la valeur 5 du pixel. Lorsque l'image est monochromatique, la valeur d'un pixel est un nombre entier. Lorsque l'image est multi-spectrale, la valeur d'un pixel est un vecteur à plusieurs composantes. Une étiquette identique est attribuée aux pixels qui présentent exactement les mêmes valeurs. Une étiquette différente est attribuée aux pixels qui présentent des valeurs différentes. Une étiquette 10 est une valeur numérique, par exemple prise dans l'ensemble des nombres entiers {0,..,Mp} où Mp est le nombre de pixels du voisinage V. Lorsque l'image I est une image multi-spectrale, la comparaison des valeurs des pixels est effectuée sur toutes les composantes d'un pixel. Autrement dit, une étiquette identique est attribuée aux pixels qui présentent entre eux les 15 mêmes valeurs pour chaque composante spectrale. Sur l'exemple de la figure 3, les pixels Po et P1 ont la même valeur, on leur attribue l'étiquette 0. Les pixels P2 et P4 ont la même valeur, on leur attribue l'étiquette 1. Le pixel P3 a une valeur différente de tous les autres pixels, on lui attribue l'étiquette 2. Les pixels P5,136,P7 ont la même valeur, on 20 leur attribue l'étiquette 3. Les étiquettes 0,1,2,3 correspondent ainsi à quatre valeurs différentes de pixels. A l'issue de l'étape 202, on obtient donc une chaine d'étiquettes E correspondant à une suite de nombres. Dans l'exemple de la figure 3, la chaine E est égale à {0 ;0 ;1 ;2 ;1 ;3 ;3 ;3}, en considérant le pixel Po comme 25 le premier élément de la chaine, ce choix étant fait par convention. N'importe quel pixel du voisinage V peut être pris comme premier élément de la chaine E, dans la mesure où cette convention est conservée pour tous les pixels de l'image I. La dernière étape 203 du procédé consiste à attribuer un identifiant à 30 la chaine d'étiquettes E, cet identifiant étant le descripteur local du pixel P. 8 Autrement dit, il s'agit d'identifier la chaine d'étiquettes E avec l'un des motifs de la base de motifs. Pour cela, on utilise un outil de mathématique combinatoire appelé collier ou « necklace » en anglais. Le document [2] décrit le principe de cet 5 outil mathématique. Un collier de k perles de longueur n est une classe d'équivalence mathématique de suites de n symboles sur un alphabet de taille k, avec k inférieur ou égal à n, invariante par rotation. Autrement dit, un collier est une classe d'équivalence qui comprend plusieurs représentants qui sont tous obtenus par rotation les uns par rapport aux autres. 10 Un bracelet est un collier particulier, encore appelé collier libre ou collier réversible. Il correspond à une classe d'équivalence sous les trois opérations de rotation (ou décalage circulaire), de symétrie (ou retournement) et de permutation. Un collier non-étiqueté est un autre collier particulier. Il correspond à 15 une classe d'équivalence sous les deux opérations de rotation (ou décalage circulaire) et de permutation. Par la suite, on décrit l'invention en utilisant comme classe d'équivalence, les colliers non-étiquetés. Cependant, l'invention s'applique 20 de façon identique en utilisant comme classe d'équivalence les colliers simples, les bracelets ou toute autre classe d'équivalence comprenant plusieurs représentants obtenus par une transformation mathématique les uns par rapport aux autres. 25 La figure 4 représente, pour illustration, tous les colliers non-étiquetés possibles pour n=k=4. Un collier non-étiqueté est une classe d'équivalence qui comprend plusieurs représentants. Sur la figure 4, on a illustré, pour chaque collier non-étiqueté, l'un de ses représentants. On obtient 7 colliers non-étiquetés différents possibles invariants par rotation et par permutation 30 des étiquettes. Par exemple le collier non-étiqueté 0 comprend le représentant {0 ;0 ;0 ;0} illustré sur la figure 4 mais aussi les représentants 9 {1 ;1 ;1 ;1}, {2 ;2 ;2 ;2}, {3 ;3 ;3 ;3}. De même, le collier non-étiqueté 6 comprend le représentant {0 ;1 ;2 ;3} illustré sur la figure 4 mais aussi, de manière non exhaustive, les représentants {1 ;2 ;3 ;0}, {2 ;3 ;0 ;1}, {3 ;0 ;1 ; 2}. 5 On peut, de façon similaire, déterminer les colliers non-étiquetés possibles pour n=k=6 ou n=k=8 ou toute autre valeur de taille de voisinage V. Le descripteur local est pris égal à l'identifiant d'un représentant de la classe d'équivalence, qui est ici un collier non-étiqueté, à laquelle la chaine d'étiquettes E du voisinage V appartient. Par exemple, l'identifiant d'un 10 représentant est égale à un nombre obtenu à partir de la chaine d'étiquettes du représentant. Par exemple, le représentant choisi est celui qui donne l'identifiant le plus petit. Sur la figure 3, on a représenté un premier représentant E1 du collier non-étiqueté auquel appartient la chaine d'étiquettes E, ce premier 15 représentant étant obtenu par décalage circulaire de deux valeurs. On a ensuite représenté le représentant E0 choisi pour identifier le collier non-étiqueté et qui est obtenu par permutation des étiquettes du premier représentant El. L'identifiant du représentant E0 est Id=00011232. Ensuite, on transforme l'identifiant Id en un descripteur local DL=84 qui est le rang du 20 collier non-étiqueté parmi l'ensemble des colliers possibles. L'ensemble des colliers possibles est rangé dans un ordre prédéterminé, par exemple un ordre lexicographique ou tout autre ordre. Les étapes 201,202,203 sont itérées pour tous les pixels de l'image I 25 ou de la zone de l'image I à représenter. Lorsque tous les pixels sont parcourus, on obtient un descripteur global de l'image I ou zone de l'image I qui correspond à la concaténation de tous les descripteurs locaux calculés. Ce descripteur global peut être représenté sous la forme d'une image DG dont chaque pixel est égal au 30 descripteur local calculé. Alternativement, le descripteur global peut être un histogramme des valeurs calculées des descripteurs locaux. 10 Comme indiqué, le procédé décrit aux figures 2 et 3 s'applique à l'identique en remplaçant la classe d'équivalence de type collier non-étiqueté par une classe d'équivalence de type collier simple, bracelet ou autre classe 5 d'équivalence connue de l'Homme du métier. Une classe d'équivalence de type collier simple comprend toutes les chaines d'étiquettes invariantes par décalage circulaire. Par exemple, pour une telle classe, la chaine {0 ;0 ;1 ;2} est équivalente à la chaine {2 ;0, 0 ;1}. Une classe d'équivalence de type bracelet comprend toutes les 10 chaines d'étiquettes invariantes par décalage circulaire, symétrie et permutation des étiquettes. Par exemple, pour une telle classe, la chaine {0 ;0 ;1 ;2} est équivalente à la chaine {0 ;1 ;2 ;2} obtenue par permutation d'étiquettes de la chaine {2 ; 1 ; 0 ; 0}. Un autre exemple est donné par les chaines {0 ;0 ;1 ;0 ;1 ;2} et 15 10 ;0 ;1 ;2 ;0 ;21 qui appartiennent à deux colliers non-étiquetés différents mais au même bracelet. En effet, si l'on écrit la chaine {0 ;0 ;1 ;0 ;1 ;2} à l'envers (de droite à gauche), on obtient la chaine {2 ;1 ;0 ;1 ;0 ;0} à laquelle on peut appliquer une rotation pour obtenir la chaine {0 ;0 ;2 ;1 ;0 ;1} puis une permutation des étiquettes pour obtenir la chaine {0 ;0 ;1 ;2 ;0 ;2}. 20 Un objectif du calcul d'un descripteur, selon l'invention, est de pouvoir caractériser le contenu d'une image, en particulier la texture, les contours ou toute information contenue dans l'image dans le but de permettre une classification du contenu. 25 On peut remarquer que l'information apportée par le descripteur dépend de la répartition des valeurs des pixels de l'image. Par exemple, si tous les pixels d'une image ont la même valeur, tous les descripteurs locaux auront également la même valeur. A l'inverse, si tous les pixels d'une image ont tous des valeurs différentes, tous les descripteurs locaux auront 30 également la même valeur puisque le calcul d'un descripteur local dépend des différences entre pixels d'un voisinage. 11 On voit donc que la quantité d'information maximale apportée par le descripteur est obtenue pour une image intermédiaire entre une image dont tous les pixels sont identiques et une image dont tous les pixels sont différents. 5 Parce que les valeurs des pixels et de leurs différences dépendent du temps de pose, le descripteur local ou motif extrait en dépend également ; on s'intéresse alors ici à la variation locale du motif au cours du temps de pose. A partir de cette constatation, on propose un deuxième mode de 10 réalisation de l'invention illustré par l'organigramme de la figure 5. Dans ce deuxième mode, on simule une variation du temps de pose ou durée d'exposition, de l'image I, à différentes valeurs temporelles, puis on calcule un descripteur global de chacune des images obtenues, selon la méthode décrite aux figures 2 et 3. 15 L'analyse en motifs du contenu d'une image dépend fortement de la dynamique utilisée pour numériser les pixels, autrement dit le nombre de bits par pixels, mais aussi du nombre de canaux dans le cas d'une image multispectrale. Par ailleurs, elle dépend également de la durée d'exposition de l'image lors de sa prise de vue. En effet, les contours des objets et formes 20 présents dans la scène apparaissent progressivement en fonction de la durée d'exposition, lors de la prise de vue. De la même manière, les motifs qui peuvent être extraits d'une image à partir des descripteurs calculés selon l'invention, varient également en fonction de la durée d'exposition de l'image lors de son acquisition. 25 En partant de cette constatation, il est proposé, à partir d'une image I acquise avec un temps d'exposition donné T, de simuler la même image qui aurait été obtenue avec un temps d'exposition inférieur. Dans ce deuxième mode de réalisation de l'invention, on introduit ainsi une étape supplémentaire 204 qui consiste à générer une seconde image lt, 30 par simulation de l'acquisition de la même scène que l'image originale I, mais pour un temps d'exposition t inférieur au temps d'exposition T initial. 12 On calcule un descripteur global pour chacune des images It simulées, selon la méthode décrite à la figure 2. Le temps d'exposition t varie entre 0 et le temps d'exposition T initial. 5 Différentes méthodes peuvent être envisagées pour simuler une image It à un temps d'exposition inférieur au temps d'exposition T réel. Par exemple, un modèle linéaire de quantification peut être utilisé pour calculer la valeur Pt de chaque pixel de la nouvelle image lt à partir de la valeur P du pixel correspondant de l'image acquise I en utilisant la relation : 10 Pt= I(t/T).P1, où 1.1 désigne la fonction partie entière Plus généralement, le modèle de simulation peut être de tout type, il peut être linéaire ou non linéaire et est basé sur une modélisation de l'imageur utilisé pour la prise de vue de l'image I qui permet de restituer 15 l'image à différents temps d'exposition. Le nombre d'images simulées est un paramètre de l'invention. Les valeurs des temps d'exposition sont choisies pour suivre une distribution prédéterminée dans l'intervalle [ 0 ;T ]. Par exemple, la distribution peut être une distribution en T*(i/j) avec i et j deux nombres entiers premiers entre eux, 20 i étant strictement inférieur à j. Cette distribution présente l'avantage de comprendre toutes les valeurs possibles de temps d'exposition entre 0 et T et ainsi de permettre une détection de tous les changements d'intensité en fonction du temps de pose. Cette distribution fournit la version la plus détaillée possible de la variation dans le temps d'un descripteur local. 25 Pour limiter le nombre d'images à simuler, la distribution choisie peut être une distribution en T/i, avec i variant de 1 à une valeur maximum I. Une telle distribution est adaptée à la statistique des intensités des pixels dans une image, en particulier pour la distribution en intensité des scènes naturelles. 30 A la fin de l'exécution du procédé décrit à la figure 5, on obtient un descripteur global pour chaque temps d'exposition simulé. 13 Le résultat du procédé peut être produit 205 sous forme de séquence d'images ou séquence vidéo comprenant un descripteur global, représenté sous la forme d'une image, évoluant en fonction du temps entre 0 et le temps d'exposition T de l'image I. 5 Le résultat du procédé peut aussi être produit 205 sous la forme d'un histogramme à trois dimensions dont la première dimension correspond aux différentes classes d'équivalences possibles (selon le voisinage et le type de classe choisi pour calculer un descripteur local), la deuxième dimension est le temps d'exposition et la troisième dimension est l'intensité de la densité 10 des descripteurs locaux qui peut être exprimée en probabilité ou en occurrence d'apparition de chaque classe d'équivalence parmi les descripteurs locaux calculés. La figure 6 illustre un exemple d'histogramme à trois dimensions, pour un exemple d'image particulier et pour des descripteurs calculés à partir de 15 colliers non-étiquetés de taille n=k=6. La figure 7 illustre, pour une image particulière, un exemple de séquence de descripteurs globaux représentés sous forme d'images pour différents temps d'exposition. Une image comporte, pour chaque pixel, le descripteur local associé qui a une valeur numérique comprise entre 0 et le 20 nombre maximum de classes d'équivalences possibles. Le résultat produit par le procédé selon le deuxième mode de réalisation de l'invention peut être restitué en étant affiché sur un écran ou toute interface homme machine. 25 Pour exploiter au mieux la dynamique et le spectre de l'image, l'image correspondant au temps d'exposition maximum simulé doit être une image comprenant des pixels de valeurs toutes différentes avec un nombre de bits par pixels et un nombre de canaux minimum. Autrement dit, l'image correspondant au temps d'acquisition effectif Tpose est acquise de telle sorte 30 à ce que tous ses spectres différent significativement par rapport à un modèle standard de bruit. Cela conduit à choisir un nombre de canaux 14 spectraux et une sensibilité photométrique minimale. L'invention convient également pour des images présentant une moindre variabilité des pixels au temps d'acquisition Tpose mais conduit alors à une description moins détaillée de la scène. C'est par exemple le cas des images RGB qui présentent 5 typiquement plusieurs fois la même couleur à plusieurs endroits. Ainsi, dans une variante du deuxième mode de réalisation de l'invention, on détermine 206, à partir de l'image I acquise initialement, une autre image qui correspond au temps d'acquisition effectif ou temps de pose effectif T pose - 10 Pour cela, l'image I est réduite à une dynamique prédéterminée, par exemple 256 bits par pixel et à un seul canal. Ensuite, on augmente progressivement le nombre de canaux spectraux jusqu'à obtenir une image multispectrale pour laquelle tous les pixels sont différents. Alternativement, cette condition peut être limitée aux voisinages de chaque pixel de l'image. 15 Dans cette alternative, on augmente le nombre de canaux spectraux jusqu'à obtenir une image pour laquelle tous les pixels sont différents au sein de chaque voisinage d'un pixel. Dans une autre variante, la dynamique des pixels peut être différente pour différents canaux spectraux. Dans encore une autre variante, la 20 dynamique, c'est-à-dire le nombre de bits par pixel peut être fixé pour l'ensemble des canaux, le nombre de bits par canal étant ainsi diminué lorsque le nombre de canaux augmente. Une fois que l'image correspondant au temps de pose effectif Tpose est 25 acquise l'image pour un temps T inférieur à Tpose est simulé en prenant en compte les caractéristiques connues du capteur de l'image, suivant le procédé décrit à la figure 5. La figure 8 schématise, sur un organigramme, une troisième variante 30 de réalisation de l'invention. Dans cette variante, le procédé tel que décrit à la figure 5 est repris à l'identique puis une étape supplémentaire 207 de 15 calcul de l'entropie de chaque descripteur global est ajoutée. L'entropie de Shannon H(t) de la densité des intensités des descripteurs locaux peut être calculée à l'aide de la relation suivante : H (t) = ENki e i p(k, t)log [ p(k, 0], avec p(k,t) la densité d'état des 5 descripteurs locaux, k variant de 1 à Nice OÙ Nice est le nombre de classes d'équivalence différentes. La figure 9 représente, sur trois diagrammes, les entropies en fonction du temps d'exposition ramené au temps d'exposition maximal normalisé (t 10 varie entre 0 et 1). Les trois diagrammes correspondent respectivement, de haut en bas, à des classes d'équivalences de type colliers non-étiquetés de paramètres respectifs n=k=4 ; n=k=6 et n=k=8. Sur chaque diagramme, on a représenté quatre courbes d'entropie calculées respectivement pour la résolution originale de l'image, puis pour des résolutions diminuées 15 respectivement d'un facteur 2, 4 et 8. Les maxima de chaque courbe sont identifiés par un point. On remarque que l'entropie présente un maximum relativement stable en fonction de la résolution et de la taille du collier non-étiqueté. On voit donc qu'il existe un temps d'exposition optimal pour lequel l'entropie des densités 20 de descripteurs locaux est maximale. On peut en déduire que le descripteur global associé à ce temps d'exposition optimal est celui qui comporte le plus d'information relative au contenu de l'image initiale. On remarque aussi que la valeur de l'entropie maximale augmente avec la taille du collier non-étiqueté. 25 Dans une dernière étape 208 du procédé décrit à la figure 8, on sélectionne le descripteur global correspondant au temps d'exposition optimal pour lequel l'entropie des descripteurs locaux est maximale. La figure 10 illustre une comparaison de résultats obtenus avec d'une 30 part un calcul de descripteurs selon la méthode de l'art antérieur LBP (image 16 101) et d'autre part un calcul de descripteurs selon l'invention (images 110160), pour un exemple d'image. L'image 100 représentée sur la figure 10 est l'image originale. L'image 101 est l'image composée des descripteurs locaux calculés 5 pour chaque pixel de l'image 100 par application de la méthode LBP. Les images 110-160 sont composées des descripteurs locaux calculés pour chaque pixel de l'image 100 par application de l'invention et en faisant varier le temps de pose avec les valeurs indiquées en dessous des images 110160. 10 Une application possible du calcul de descripteurs de l'image 100 consiste à analyser les descripteurs pour segmenter et classifier certaines zones de l'image et ainsi permettre une reconnaissance automatique de certaines zones de l'image. Dans l'exemple de la figure 10, l'image d'origine 100 est une image d'un visage. Dans cet exemple, on cherche à segmenter 15 les éléments du visage, c'est-à-dire les yeux, les narines, les oreilles et les cheveux, notamment. Cette segmentation permet ensuite d'effectuer une reconnaissance automatique du visage de l'image 100, à partir des descripteurs calculés. Sur l'image 101, obtenue par application de la méthode LBP, on 20 remarque que les informations relatives aux différentes caractéristiques du visage sont mélangées. En outre, certains éléments du visage ne sont pas segmentables, par exemple les narines, car les descripteurs calculés pour ces éléments sont trop proches de descripteurs calculés pour d'autres parties du visage. Cela est du, notamment, au fait que la méthode LBP ne 25 permet pas de sélectionner un temps de pose optimal pour obtenir un niveau d'information maximal pour les descripteurs calculés. A l'inverse, la méthode selon l'invention permet de générer plusieurs descripteurs 110-160 pour différents temps de pose. Grâce à cet aspect de l'invention, un temps de pose optimal peut être sélectionné en fonction de 30 l'information que l'on cherche à classifier dans l'image. Par ailleurs, l'invention permet aussi de gérer de manière conjointe les différents canaux 17 spectraux de l'image ce qui favorise une meilleure segmentation des détails de l'image. Ainsi, on remarque par exemple que le descripteur représenté par l'image 140 permet une meilleure segmentation des narines que le descripteur 5 représenté par l'image 101. Sans sortir du cadre de l'invention, de nombreuses applications sont envisageables pour analyser des caractéristiques d'une image de façon automatique, à partir des descripteurs calculés. Par exemple, il est possible de comparer deux images en comparant les histogrammes des descripteurs 10 calculés pour chacune des deux images. Le procédé selon l'invention peut être mis en oeuvre en tant que programme d'ordinateur, le procédé étant appliqué à une image I préalablement acquise à l'aide d'un capteur d'image, par exemple un 15 appareil photo, une caméra ou une caméra multispectrale. L'invention peut être mise en oeuvre en tant que programme d'ordinateur comportant des instructions pour son exécution. Le programme d'ordinateur peut être enregistré sur un support d'enregistrement lisible par un processeur. La référence à un programme d'ordinateur qui, lorsqu'il est 20 exécuté, effectue l'une quelconque des fonctions décrites précédemment, ne se limite pas à un programme d'application s'exécutant sur un ordinateur hôte unique. Au contraire, les termes programme d'ordinateur et logiciel sont utilisés ici dans un sens général pour faire référence à tout type de code informatique (par exemple, un logiciel d'application, un micro logiciel, un 25 microcode, ou toute autre forme d'instruction d'ordinateur) qui peut être utilisé pour programmer un ou plusieurs processeurs pour mettre en oeuvre des aspects des techniques décrites ici. Les moyens ou ressources informatiques peuvent notamment être distribués ("Cloud computing"), éventuellement selon des technologies de pair-à-pair. Le code logiciel peut 30 être exécuté sur n'importe quel processeur approprié (par exemple, un microprocesseur) ou coeur de processeur ou un ensemble de processeurs, 18 qu'ils soient prévus dans un dispositif de calcul unique ou répartis entre plusieurs dispositifs de calcul (par exemple tels qu'éventuellement accessibles dans l'environnement du dispositif). Le code exécutable de chaque programme permettant au dispositif programmable de mettre en 5 oeuvre les processus selon l'invention, peut être stocké, par exemple, dans le disque dur ou en mémoire morte. De manière générale, le ou les programmes pourront être chargés dans un des moyens de stockage du dispositif avant d'être exécutés. L'unité centrale peut commander et diriger l'exécution des instructions ou portions de code logiciel du ou des 10 programmes selon l'invention, instructions qui sont stockées dans le disque dur ou dans la mémoire morte ou bien dans les autres éléments de stockage précités. Les résultats du procédé selon l'invention peuvent être affichés sur un écran ou toute autre interface homme machine. 15 19 Références [1] M.Pietikâinen, A. Hadid, G. Zhao, and T. Ahonen, Local binary patterns for still images, Springer, 2011. 5 [2] Jean-Philippe Gagnon, « Colliers et bracelets dont les perles importent peu », Faculté des sciences et de Génie, REVENDICATIONS1. Procédé, mis en oeuvre par ordinateur, de calcul d'un descripteur global (DG) d'au moins une partie d'une image à décrire (I), le descripteur global (DG) comprenant un descripteur local (DL) associé à chaque pixel de l'au moins une partie d'image à décrire, le procédé comprenant le calcul de tous les descripteurs locaux d'un descripteur global en : - Définissant (201) un voisinage du pixel associé au descripteur local, - Attribuant (202) une étiquette à chaque pixel du voisinage, deux pixels de valeurs différentes ayant deux étiquettes différentes, deux pixels de valeurs identiques ayant la même étiquette, l'ensemble des étiquettes du voisinage formant un motif (E), - Affectant (203) la valeur du descripteur local à un représentant d'une classe d'équivalence à laquelle appartient le motif (E). 2. Procédé de calcul d'un descripteur global selon la revendication 1 comprenant au préalable les étapes de : - Obtenir (206) une première image acquise avec un premier temps d'exposition, - Simuler (204) au moins une seconde image résultant d'une prise de vue de la première image avec au moins un second temps d'exposition inférieur au premier temps d'exposition, - Calculer (202,203) un descripteur global pour au moins une partie d'une image à décrire prise parmi la première image acquise ou l'au moins une seconde image simulée. 3. Procédé de calcul d'un descripteur global selon la revendication 2 dans lequel l'étape de simuler (204) une seconde image est réalisée en 21 diminuant l'intensité des pixels de la première image d'un facteur égal au rapport entre le premier temps d'exposition et le second temps d'exposition. 5 4. Procédé de calcul d'un descripteur global selon l'une des revendications 2 ou 3 dans lequel l'au moins une partie de la première image acquise est composée de pixels dont les valeurs sont toutes différentes entre elles. 5. Procédé de calcul d'un descripteur global selon la revendication 4 10 comprenant une étape de génération (206) de la première image acquise à partir d'une image multi-spectrale en sélectionnant le nombre minimal de composantes spectrales de l'image nécessaire pour obtenir que l'au moins une partie de la première image acquise est composée de pixels dont les valeurs sont toutes différentes entre elles. 15 6. Procédé de calcul d'un descripteur global selon l'une des revendications 2 à 5 comprenant les étapes de faire varier le second temps d'exposition entre le premier temps d'exposition et zéro et de calculer un descripteur global pour au moins une partie de la première image acquise et de l'au 20 moins une seconde image simulée. 7. Procédé de calcul d'un descripteur global selon la revendication 6 comprenant une étape d'affichage (205) des descripteurs globaux calculés en fonction du second temps d'exposition. 25 8. Procédé de calcul d'un descripteur global selon l'une des revendications 6 ou 7 comprenant un calcul d'entropie (207) des descripteurs locaux pour chaque descripteur global calculé et la sélection (208) d'au moins un descripteur global d'au moins une image à décrire en fonction des 30 entropies calculées. 9. Procédé de calcul d'un descripteur global selon la revendication 8 comprenant la sélection (208) du descripteur global dont l'entropie calculée est la plus élevée. 22 10. Procédé de calcul d'un descripteur global selon l'une des revendications précédentes dans lequel la classe d'équivalence est une classe de chaines d'étiquettes invariantes par rotation. 5 11. Procédé de calcul d'un descripteur global selon la revendication 10 dans lequel la classe d'équivalence est une classe de chaines d'étiquettes qui est en outre invariante par permutation de valeurs d'étiquettes. 10 12. Procédé de calcul d'un descripteur global selon l'une des revendications 10 ou 11 dans lequel la classe d'équivalence est une classe de chaines d'étiquettes qui est en outre invariante par symétrie. 13. Procédé de calcul d'un descripteur global selon l'une des revendications 15 précédentes dans lequel une image à décrire est une image multi- spectrale dont chaque pixel comporte plusieurs composantes à différentes longueurs d'ondes. 20 14. Programme d'ordinateur comportant des instructions pour l'exécution du procédé de calcul d'un descripteur global selon l'une quelconque des revendications 1 à 13, lorsque le programme est exécuté par un processeur. 25 15.Support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution du procédé de calcul d'un descripteur global selon l'une quelconque des revendications 1 à 13, lorsque le programme est exécuté par un processeur. 30