FR 2841357 A1 20031226 FR 0207759 A 20020624 I'acces, vent des sites juridiques eVou economiques concernant l'immobilier. 2841 357 PROCEDE DE RESTAURATION SPATIALE DES COULEURS D'UNE IMAGE La presente invention appartient au domaine technologique de l'imagerie. La presente invention a pour objet un procede automatique de s restauration de la couleur d'une image, apres que ladite couleur ait ete alteree par rapport a la couleur d'origine de ['image. La renovation ou restauration des images en couleur enregistrees sur des supports argentiques, est generalement operee en faisant un traitement d'une version numerisee de ces images argentiques alterees. La restauration des images alterees permet qu'a liexamen visuel, ['image recouvre ses couleurs d'origine. Les supports traditionnels argentiques des images vent constitues generalement par du film photographique, ou du papier photographique. On designee dans la description qui suit, par image argentique, une image en couleur obtenue a partir d'emulsions photographiques aux halogenures d'argent. L'alteration ou degradation des images en couleur enregistrees sur ces supports argentiques est due par exemple au vieillissement, ou a ['exposition de ces supports argentiques a des sources lumineuses, a ['occasion de leur manipulation. Les procedes de restauration des couleurs d'une image, connus de l'homme du metier, consistent par exemple a numeriser une image prealablement enregistree sur un support argentique (film ou papier photographique), puis a traiter ['image numerisee afin de restaurer ses couleurs. Ces pro cedes de restauration consistent a transformer les couleurs degradees sur ['image numerisee, en mettant en ccuvre des algorithmes de logiciels de traitement dtimage. L'image numerique, transformee par le procede de traitement, 2s presente ainsi des couleurs restaurees qui correspondent approximativement a celles de ['image d'origine. Le brevet US 5 796 874 decrit la restauration d'images imprimees sur un support materiel, par exemple de type papier, decolorees ou fanees au cours du temps. Ce brevet decrit une methode et des moyens qui permettent a un operateur n'ayant pas un niveau de competence eleve, de restaurer et facilement une image dont les couleurs ont ete alterees. L'operateur n'a notamment pas a choisir et a ajouter des filtres dans le scanner utilise pour numeriser ['image a restaurer. Mais, l'operateur doit appliquer des model es de degradation de [ 'image au cours du temps a [ 'image a restaurer. C'est- a-dire que la restauration est op eree en fonction du choix d'une variable de temps: par exemple, on suppose que ['image s a une dizaine d'annees d'age et on applique le modele de restauration correspondent a ce temps (une dizaine d'annees) afin de restaurer les couleurs de ['image alteree. La methode de restauration permet d'obtenir et de selectionner automatiquement un modele de restauration des couleurs de ['image, en fonction d'une variable temporelle affectee a ['image. L'image ainsi restauree peut etre affichee sur un 0 ecran de type moniteur. De plus, l'operateur peut aussi choisir une autre modele, lie a un autre variable de temps, si le premier resultat de restauration ne le satisfait pas. La regle de correction utilisee est appliquee uniformement a ['ensemble de ['image. Le document de l'Universite de La Rochelle, France, de M. Chambah et B. Besserer, intitule Digital Color Restoration of Faded Motion Pictures, et presente a la premiere Conference Internationale CGIP (Color in Graphics and Image Processing) de Saint-Etienne, France, du 01 au 04 octobre 2000, decrit une methode de restauration des couleurs dtimages fanees sur de vieux films. Dans un premier temps, le film est numerise avec un scanner; puis l'effet du recouvrement des densites spectrales des differents colorants est soustrait, en utilisant une matrice d'ajustement. Dans un second temps, les canaux de couleurs de ['image vent equilibres, en utilisant une autre matrice de correction. Finalement, le contraste est augmente, pour ameliorer la qualite visuelle. Dans cette methode, la regle de correction utilisee est appliquee uniformement a ['ensemble de ['image. Dans un article (IEEE Transactions On Image Processing, volume 9, numero 5, de mai 2000), intitule "Adaptive Image Contrast Enhancement Using Generalizations of Histogram Equalization", J. Alex Stark decrit une methode d'augmentation du contraste dans une image. Cette methode utilise une fonction cumulative des niveaux de gris dans une zone ou une fenetre (par exemple de forme carree) autour d'un pixel, sans contrainte d'interdependance entre ces zones. Les corrections vent independantes, pour chaque zone de representation de ['image. Les moyens de restauration de la technique anterieure ne prennent pas en compte les variations spatiales de ['alteration des couleur dans une image. Les moyens de la technique anterieure utilisent des modeles ou fonctions mathematiques permettant de prendre en compte certains paramekes importants qui influencent la degradation des couleurs d'une image enregistree sur un film ou un papier photographique, pour ensuite restaurer les couleurs de ['ensemble de ['image d'une maniere uniforme. Le procede de la presente invention, pour la restauration spatiale des couleurs d'une image alteree, permet, a partir de la detection des alterations ou lo degradations de couleur qui varient spatialement d'un bord a l'auke de ['image, une restauration des couleurs conformement aux couleurs d'origine de ladite image. L'avantage de ce procede de ['invention est de restaurer automatiquement, rapidement, et de maniere fiable les couleurs d'origine de ['image, alors que lesdites couleurs ont ete alterees non uniformement. 1 S Plus particulierement, ['invention concerne un pro cede automatique de restauration de la couleur d'une image argentique dont la couleur est alteree d'une facon non uniforme par rapport a la couleur d'origine de ['image argentique, spatialement dans le plan de ['image argentique, et suivant une direction parallele a un axe principal de ['image; ce procede comprend les etapes suivantes: a) numeriser ['image argentique alteree en utilisant une echelle de mesure des densites optiques; b) filtrer ['ensemble des pixels de ['image numerisee; c) diviser ['image numerique en bandes de pixels disposees perpendiculairement a la direction de l'axe principal de ['image 2s suivant lequel se produit ['alteration de couleur, la reunion de ['ensemble des pixels de chaque bande representant ['ensemble des pixels de ['image numerique; d) calculer, pour chaque bande de pixels, et pour chacun des canaux de couleur formant ['image, une distribution des densites optiques des pixels formant ladite bande; e) calculer, pour chacune desdites distributions de l'etape d), une valeur de reference brute superieure de densite optique, et une valeur de reference brute inferieure de densite optique; f) calculer, a partir de ['ensemble des valeurs de reference brutes superieures et inferieures de l'etape e), des valeurs de reference filtrees superieures de densite optique et des valeurs de reference filtrees inferieures de densite optique correspondantes; g) determiner, a part* des valeurs de reference filtrees superieures et inferieures de densite optique de l' etape f), des courb es de o reference superieures de densite optique et des courbes de reference inferieures de den s it e opti que, l es dites courb es de reference etant representatives du profil de degradation de couleur de ['image; h) determiner, pour chacun des canaux de couleur formant ['image, s et a partir de la valeur maximale de densite optique de la courbe de reference superieure de densite optique, la bande de pixels de moindre degradation; i) appliquer a chaque bande de pixels de ['image numerique autre - que la bande de pixels de moindre degradation, et pour chacun des canaux de couleur formant ['image, une kansformation lineaire permettant de transformer les valeurs filtrees de densite optique placees sur la courbe de reference superieure et sur la courbe de reference inferieure de densite optique de ladite chaque bande de pixels, respectivement en valeurs superieure et 2s inferieure de densite optique de la bande de pixels de moindre degradation. Le procede de ['invention permet d'obtenir une image restauree. Une caracteristique de ['invention qui se traduit par un avantage majeur par rapport a la technique anterieure, est de prendre en compte les variations spatiales non uniformes de couleur ayant altere par exemple une image argentique. Ce procede permet de restaurer automatiquement les couleurs fanees de ['image, sans faire intervenir les competences photographiques d'un operateur. Un autre avantage du procede de ['invention pour executer des restaurations de couleur d'images alterees ou fanees est le temps gagne, et par consequent une productivite accrue, par rapport a un traitement manuel de la restauration spatiale des images par un s operateur. D'autres caracteristiques et avantages apparaitront a la lecture de la description, faite en reference aux dessins des differentes figures. La figure 1 represente schematiquement un film photographique en couleur et un moyen destine a son enroulement. lo La figure 2 represente schematiquement un exemple d'environnement materiel (hardware) dans lequel est par exemple mis en ccovre le procede de ['invention. La figure 3 represente schematiquement une image numerisee de ['image a restaurer. La figure 4 represente un exemple de distribution, dans un canal de couleur, des densites optiques d'une bande de pixels d'une image numerique. La figure 5 represente schematiquement, pour un canal de couleur, une evolution de la distribution des densites optiques dans ['image, selon un axe principal suivant lequel se produit ['alteration de couleur. La figure 6 represente schematiquement un exemple de diagramme, dans chaque canal de couleur, de representations spatiales de la distribution des densites optiques dans ['image filtree de son contenu. La figure 7 represente la ffinction de restauration des densites optiques. 2s La figure 8 represente un exemple de distribution des densites optiques cumulees de la figure 4. La figure 9 represente un principe de filtrage de ['image par suivi des maximum. La figure 10 represente schematiquement, pour un canal de couleur, la distribution des densites optiques apres filtrage, de la figure 5. La description suivante est une description detaillee de ['invention en reference aux dessins dans lesquels les memes references numeriques identifient les memes elements dans chacune des differentes figures. Le procede de ['invention est utilise pour restaurer automatiquement s les couleurs d'une image dont la couleur est alteree, c'est-a-dire fanee par rapport a la couleur d'origine de ['image. Les images a restaurer peuvent etre fixes ou animees. Les images fixes vent par exemple des photographies en couleur (portraits, paysages, etc.). Les images animees vent par exemple des sequences de film cinematographiques. Ces images vent enregistrees sur des supports lo photographiques argentiques. Le probleme d'alteration non homogene des couleurs des images enregistrees sur un support photographique argentique est rencontre notamment lorsque, suivant la figure 1, ces supports d'images argentiques 2 sQnt enroules sur des bobines 1. Le support dtimages argentiques est par exemple un film cinematographique enregistre en couleur. Au cours du temps ou suite a des variations de temperature et/ou dthygrometrie non maitrisees, des degradations de colorants se produisent. Pratiquement, les degradations des colorants de ['image se produisent plus fortement sur les bords 4F, SF qu'au centre du film 2. La ligne materialisant le milieu du film 2 et des images argentiques 3 est representee par un premier axe principal 6. La figure 2 represente un environnement materiel preferential dans lequel le procede de ['invention est mis en ccuvre de maniere interactive. Un operateur qui desire par exemple restaurer une ou plusieurs images d'origine 3 dont les couleurs ont ete alterees au cours du temps par rapport a la couleur d'origine, utilise un terminal 7. Le terminal 7 est par exemple un ordinateur personnel ou PC 2s (Personal Computer) comprenant un ecran d'affichage 11. Le terminal 7 est interface a une souris 9 comprenant au moins un bouton de commande 10. Le systeme d'exploitation du terminal 7 est par exemple Unix, Windows, ou Linux. Tous ces systemes d'exploitation permettent de reconna^tre des peripheriques comme un scanner ou numeriseur 13, une imprimante 14, ou tout autre dispositif 15 relic au terminal 7, et permettant la capture d'images ou de sequences d'images. Le procede de ['invention permet de numeriser automatiquement ['image argentique 3 dont les couleurs vent alterees, en utilisant par exemple le scanner 13. L'image 3 numerisee est codee en valeurs numeriques. Pour numeriser ['image, on utilise une echelle de mesure des densites optiques. Autrement dit, les valeurs numeriques associees a chaque pixel, et pour chacun des canaux de couleur s formant ['image, vent proportionnelles (fonction lineaire) aux densites optiques. L'image 3 est numerisee en une image 12 qui peut ainsi 8tre visualisee sur l'ecran 11. L'image numerique 12 peut 8tre memorisee dans une memoire interne au terminal 7. Dans un autre mode de realisation, la mise en ceuvre du procede de ['invention peut aussi ne pas 8tre interactive, et ntutilise pas forcement des moyens comprenant un ecran d'affichage 11. Dans ce mode de realisation, le procede de ['invention est mis en ccuvre en utilisant un processeur dedie, c'est-a-dire un processeur configure avec des portes et des fonctions logiques de base permettant de mettre en aeuvre ledit procede. On utilise par exemple un processeur FPGA S (Field Programmable Gate Array). Ce processeur est par exemple integre dans le scanner 13, et execute automatiquement la restauration juste apres ['operation de numerisation de ['image argentique 3. L' alteration, c 'est- a-dire l a deco l oration ou l e fanage des coul eurs d'une image argentique est causee en particulier par ['exposition du support de ['image argentique a la lumiere et/ou a l'humidite. Les films vent par exemple stockes enroules sur des bobines 1; dans ces conditions, ['alteration des couleurs de ['image argentique 3 enregistree sur le film 2 se produit d'abord au voisinage des bords 4F, 5F audit film 2. Au cours du temps, ['alteration de couleur est done plus importante au voisinage des bords 4F, SF du film 2 enrouIe sur la bobine 1, qutau 2s voisinage de l'axe principal 6 situe au centre du film 2. En consequence, une premiere hypothese fondee sur des resultats experimentaux, et prealable a la mise en aeuvre du pro cede de ['invention, consiste a admettre d'une part que ['alteration de couleur varie spatialement dans le plan de ['image 3, suivant une direction perpendiculaire au premier axe principal 6, ou parallele a un second axe principal 1 g de ['image 3, 12. Donc, ['alteration de couleur varie en fonction de la distance par rapport aux bords 4F, SF du film 2. D'autre part, cette premiere hypothese est liee a une seconde hypothese qui consiste a admettre egalement que la degradation de couleur est uniforme (c'est-a-dire ne varie pas) suivant le premier axe principal 6, parallele aux bords 4F, 5F. Cette hypothese d'uniformite de la degradation de couleur doit etre verifiee au moins sur l'etendue de ['image selon cet axe principal s 6. Selon la figure 3, et par analogie avec la figure 1, les bords 16, 17 de ['image numerique 12 (image 3 numerisee) correspondent respectivement aux bords 4I, SI de ['image argentique 3. Les bords 26, 27 de ['image numerique 12 correspondent aux deux autres bords de ['image argentique 3, lesdits autres bords o etant paralleles au second axe principal 19. Le premier axe principal 6 est perpendiculaire au second axe principal 19. Les deux axes principaux 6, 19 de ['image 3, 12 se coupent en un point 20 qui est place au centre de ['image. Pour determiner les variations des densites de couleur spatialement dans le plan de ['image numerique 12 par rapport a la couleur d'origine de ['image argentique 3 correspondante, l'algorithme du pro cede de ['invention permet de calculer automatiquement des distributions de pixels, en fonction des valeurs numeriques codees de leurs densites optiques respectives, dans une bande mobile 18 de ['image 12. Le procede de ['invention permet de calculer ces distributions dans des tranches ou bandes 18 representatives d'une partie de ['image 12. La largeur 21 de la bande 18 inclus par exemple cinq colonnes de pixels paralleles a l'axe 6. La longueur de la bande 18 correspond a ['extension totale de ['image 12 suivant la direction de l'axe 6, perpendiculaire a la direction de l'axe 19; la direction de l'axe 19 caracterise la direction de la variation de ['alteration des couleurs. L'extension totale de ['image 12 est la dimension de ['image: par exemple 2s sa hauteur ou sa largeur. Dans un mode de realisation particulier de ['invention, la bande 18 la moins large est reduite a une senle rangee de pixels parallele a l'axe 6. La figure 4 represente un exemple de distribution des densites optiques d'une bande de pixels, dans un canal de couleur donne, par exemple le rouge. Une courbe 22 represente la distribution des pixels (nombre de pixels N suivant l'axe des ordonnees) en fonction des valeurs codees de densite optique VC(DO) desdits pixels (suivant l'axe des abscissas). La courbe 22 represente, pour une position spatiale donnee de la bande 18 dans le plan de ['image 12, la distribution des valeurs codees de densite optique VC(DO) des pixels formant ladite bande 18, entre une valeur minimale de densite optique representee par le point 23 et une valeur maximale de densite optique representee par le point 24. Generalement, ['image est affectee par du bruit (signal parasite). Les causes de bruit peuvent etre multiples. Dans le contexte de ['invention, le bruit genant est du principalement a des poussieres. Les poussieres ont genere par exemple des points noirs dans ['image, ou bien des points blancs, et perturbent la distribution des valeurs codees de densite optique VC(DO). Selon la figure 4, les valeurs codees de lo densite optique des points noirs vent par exemple representees par un lobe 28, et les valeurs codees de densite optique des points blancs vent par exemple representees par un lobe 29. Le procede de ['invention permet d'eliminer automatiquement ces lobes aberrants 28, 29. La figure 8 represente un exemple de courbe 70, de la distribution des densites optiques cumulees de la figure 4. Selon un mode de realisation particulier de ['invention, la distribution cumulee de la figure 8 est par exemple divisee en class es d'effectifs egaux de pixel s ayant des valeurs voi sines de densites optiques. Preferentiellement, la distribution est divisee en centiles (100 classes). Mais la distribution peut etre divisee en un autre nombre arbikaire de classes. Pour eliminer les lobes aberrants eventuels 28, 29 dus au bruit, le procede de ['invention ne traite que la partie de la courbe 70 comprise entre les points d'ordonnees respectives X et (100 - X'). Dans un mode de realisation prefere: X = X' = 5%. Mais les valeurs de X et X' vent parametrables, c'est-a-dire peuvent prendre d'autres valeurs (pourcentage), et ne pas etre egales entre elles. Une troisieme hypothese de mise en ccuvre du procede de ['invention consiste a considerer que le changement de densite optique correspondent a ['alteration de couleur s'opere suivant un modele lineaire, caracterise par exemple par une premiere equation Eq(1) du type: DF = a DR + P. Le mode de realisation prefere de realisation de ['invention utilise le modele lineaire en l'appliquant aux densites optiques de pixels de bandes 18 de l'irnage numerique 12. Dans un autre mode de realisation, on peut egalement utiliser ce modele lineaire en l'appliquant aux intensites des couleurs. Dans ['equation DF = a DR +,B, DF represente la valeur codee de la densite optique dicn pixel altere, c'est-a-dire un pixel dont la valeur codee de s densite optique a par exemple diminue par rapport a la valeur codee de densite optique de reference DR audit pixel. La valeur codee de densite optique de reference DR correspond a la densite optique de la couleur d'origine du point correspondent spatialement dans le plan de ['image argentique 3, au pixel de ['image numerique 12. Les valeurs de a et varient selon l'axe principal 19 de lo ['image numerique 12. Une quatrieme hypothese prealable a la mise en ceuvre du procede de ['invention, admet que le profil de degradation reel de couleur de ['image varie de maniere continue (c'est- a-dire sans rupture), lisse (ctesta-dire avec une variation lente) et monotone (c'est-a-dire avec une degradation toujours decroissante du bord 1S vers le point de moindre degradation). Le point de moindre degradation de ['image 12 est represente par exemple par le point 20 place au centre de ['image 12. Mais ce point de moindre degradation peut etre place spatialement n'importe ou dans le plan de ['image 12. Le profil de degradation detecte a pour objectif de se rapprocher du profil de degradation reel. Donc, le resultat de restauration de couleur est d'autant meilleur que le profil de degradation detecte par le procede de ['invention est continu, li sse, et monotone. C'est- a-dire que le pro fil s'affranchit des variations chromatiques inherentes au contenu meme de ['image. Des moyens de filtrage connus de l'homme du metier, permettent de rendre le profil de degradation detecte, continu, lisse, et monotone. 2s Dans un mode de realisation prefere de ['invention, les etapes de detection de ['alteration des couleurs de ['image vent generalement precedees d'un filtrage initial automatique de ['ensemble des pixels de ['image numerique 12. Le filtrage initial permet de minimiser ['influence du contenu de ['image et des poussieres (artefacts dans ['image). Mais, un autre mode de realisation particulier permet de ne pas faire ce filtrage initial. Le filtrage de ['image numerique 12 est opere en utilisant des operateurs de morphologic mathematique connus de l'homme du metier. Pour chaque pixel de ['image 12, est operee une ouverture (erosion suivie d'une dilatation) puis une fermeture (dilatation suivie d'une erosion), en utilisant un element structurant ayant par exemple une forme carree et contenant 9 ou 25 pixels. Ces ouverture et fermeture peuvent etre operees consecutivement n s fois, n etant un nombre entier. Le procede de ['invention permet de diviser automatiquement ['image numerique 12 en bandes de pixels 18 disposees perpendiculairement a la direction de l'axe principal l 9 de ['image suivant lequel se produit ['alteration de couleur. La reunion de ltensemble des pixels de chaque bande 18 represente o ['ensemble des pixels de ['image 12. La division de ['image 12 en bandes de pixels 18 est operee suivant des bandes 18 paralleles adjacentes; selon un premier mode de realisation, ces bandes paralleles adjacentes n'ont aucun pixel commun entre elles, ou bien selon un second mode de realisation, la division de ['image est operee suivant des bandes paralleles se recouvrant partiellement deux a deux suivant au moins une rangee de pixels; un troisieme mode de realisation permet de diviser ['image a la fois en bandes adjacentes ou recouvrantes. Le mode de realisation choisi parmi ces trots variantes depend de l'algorithme utilise pour deplacer, dans ['image 12, la bande mobile 18 de largeur 21. Dans un mode de realisation preferential, les bandes 18 vent de meme largeur 21. Mais l'algorithme mettant en ceuvre le procede de ['invention permet egalement de choisir des bandes 18 de largeurs 21 inegales. A partir de la division de ['image numerique 12 en bandes de pixels 18, le procede de ['invention permet de calculer automatiquement, pour chaque bande de pixels 18, et selon la figure 6, pour chacun des plans des canaux colores 2s 40, 50, 60 formant ['image 12 (par exemple les canaux rouge, vert et bleu), une distribution des densites optiques de chaque pixel de ladite bande 18. Un exemple de distribution des densites optiques est represente a la figure 4. La determination des variations spatiales de ['alteration des couleurs dans ['ensemble de ['image 12, et pour chacun des canaux de couleur, par exemple rouge, vert, bleu, consiste a separer ou filtrer la contribution du contenu de ['image, de la contribution de ['alteration des couleurs. Ceci en tenant compte de la variation de la distribution des valeurs codees de densite optique, comme represente dans la figure 4. La figure 5 represente, pour un canal de couleur de ['image, une representation spatiale de la distribution reelle des densites optiques dans ['extension Xe de ladite image. L'extension Xe represente la dimension de ['image suivant laquelle se produit ['alteration ou degradation de couleur (axe principal 19). Les lignes 41 et 45 correspondent aux bords de ['image. La distribution reelle des densites optiques est representee par exemple par les courbes 30, 31, 32, 33, 34. Ces courbes 30, 31, 32, 33, 34 correspondent a un mixage du profil de degradation lo des couleurs et du contenu de ['image. Le procede de ['invention permet de separer automatiquement les signaux lies au contenu de ['image d'une part, et au profil de degradation des couleurs d'autre part. Selon la figure 5, le procede de ['invention permet de calculer automatiquement, pour chacune des distributions de densites optiques des pixels de la bande de pixels 44, une valeur de reference brute superieure de densite optique 36, et une valeur de reference brute inferieure de densite optique 37. Selon la figure 5, la distribution des densites optiques est representee par une famille de courbes 30, 31, 32, 33, 34. Les courbes 30 et 34 representent par exemple les courbes de densite optique respectivement du maximum et du minimum de la distribution. Dans un mode de realisation prefere, l'algorithme du procede de ['invention elimine automatiquement ces deux extremes de la distribution (les courbes 30 et 34) pour operer les calculs. Les valeurs de reference brute superieure 36 et inferieure 37 vent calculees a partir des courbes 31 et 32 qui correspondent par exemple respectivement aux courbes de densite optique de la famille de courb es placees dans la distribution, au voisinage respectivement des courb es extremes 30 et 34. Toutefois, dans un autre mode de realisation, les valeurs de reference brutes superieures de densite optique 39, et inferieures 38, vent calculees a partir des courbes de densite optique du maximum 30 et du minimum 34. A partir des valeurs de reference brutes superieure 36 et inferieure 37 de densite optique, le procede de ['invention permet de calculer des valeurs de reference filtrees superieure 46 et inferieure 47 de densite optique. Le pro cede de [ 'invention, a partir des valeurs de reference filtrees s 46, 47, determine automatiquement des courbes de reference de densite optique: la courbe de reference superieure de densite optique 42, et la courbe de reference inferieure de densite optique 43. Les courbes de reference de densite optique 42, 43 vent representatives du profil de degradation de couleur de ['image. Ces courbes de reference 42, 43 vent obtenues par filtrage, en tenant compte de la quatrieme hypothese liee a la mise en ceuvre de ['invention. La quatrieme hypothese est basee sur un profil de degradation de couleur de ['image qui varie de maniere continue, lisse, et monotone. Dans un mode de realisation prefere, ['operation de filtrage, dont le principe est represente sur la figure 9, est operee de maniere a appliquer cette quatrieme hypothese, eliminer les points extremes de la distribution correspondent, apres filtrage, aux courbes extremes de la distribution CO et C1OO, et separer le contenu de ['image du profil de degradation de couleur de ladite image. Les courbes CO et CIOO vent les courbes de reference inferieure et superieure de densite optique qui representent les courbes extremes de la distribution. Selon la figure 9 et dans le mode de realisation prefere de ['invention, ['operation de filtrage est operee en utilisant la methode cite du suivi des maximum. Cette methode de filtrage revient, a partir d'une courbe 31 representative d'une distribution, a obtenir la courbe resultante qui represente le minima des deux courbes 35C et 36C. (: eci permet d'eliminer les creux de la courbe initiale 31 et ainsi de respecter lthypothese de monotonic. Selon la figure 6, 2s pour chacun des canaux de couleur 40, 50, 60 formant ['image, les valeurs de reference filtrees superieures de densite optique 46, 56, 66, et les valeurs de reference filtrees inferieures de densite optique 47, 57, 67, vent calculees a partir de l'utilisation d'un filtrage de ['image 12 par la methode du suivi des maximum dans chaque bande de pixels 18. Selon les figures 6 et 10, le procede de ['invention permet ainsi diobtenir, apres filtrage, des courbes corrigees lissees 42, 43, 52, 53, 62, 63. Suivant un mode prefere de realisation, les courbes Co et CIOO, correspondent aux extremes de la distribution (0% et 100%), ont ete exclues. La figure 6 represente les courbes de reference superieures de densite optique 42, 52, 62 et respectivement inferieures 43, 53, 63. Les cites courbes de reference etant representees dans chaque plan 40, 50, 60 correspondent aux canaux de couleur. Par analogie avec la representation graphique de la figure 3, ['extension de ['image suivant laquelle se produit la variation de ['alteration de couleur est representee en abscisse sur la figure 6 par la distance comprise entre les lignes 41 et 45 pour le canal de couleur 40, entre les lignes 45 et 55 pour le canal de couleur 50, et entre lo les lignes 55 et 65 pour le canal de couleur 60. L'axe des ordonnees de la figure 6 represente les valeurs codees des densites optiques VC(DO). La couche d'emulsion photographique de ['image 3 correspondent a la couleur rouge est representee en abscisse par exemple par le plan 40 de ['image; la couche d'emulsion photographique correspondent a la couleur verse est representee en abscisse par exemple par le plan 50 de ['image; La couche d'emulsion photographique correspondent a la couleur bleue est representee en abscisse par exemple par le plan 60 de ['image. Selon la figure 6, le procede de ['invention permet de determiner automatiquement, pour chacun des canaux de couleur 40, 50, 60 formant ['image, et a partir de la valeur maximale de densite optique 46M, 56M, 66M de la courbe de reference superieure de densite optique, la bande de pixels de moindre degradation 44M, 54M, 64M. Les valeurs de reference filtrees superieures de densite optique des bandes de pixels de moindre degradation 44M, 54M, 64M vent respectivement 46M, 56M, 66M. Les valeurs de reference filtrees inferieures de densite optique des bandes de pixels de moindre degradation 44M, 54M, 64M vent respectivement 47M, 57M, 67M. Cette bande de pixels de moindre degradation correspond a la bande de pixels de ['image ayant subi une alteration ou degradation minimale de couleur. Selon un mode particulier de realisation du procede de ['invention, la bande de moindre degradation est la bande de pixels comprenant une senle rangee de pixels, et passant par le centre 20 de ['image 12. La fi gure 7 repres ente la foncti on de re st aurati on de s d ens it es optiques. La transformation lineaire est operee suivant une equation de kansformation Eq(2) du type: DR = a DF + b; equation dans laquelle a et b vent independents pour chaque canal colore rouge (ar, br), vert (av, by), bleu (ate, bb), et dans laquelle a et b vent determines par l'algorithme permettant de mettre en ccuvre le pro cede de ['invention. C ette seconde equation Eq(2) represente la transformation inverse de celle de la premiere equation Eq(l). L'application simultanee de ['equation Eq(2) aux trots canaux colores s'exprime sous la forme: o DRr = arDFr + br DRV = avDFv + bv DRb = abDFb + bb soit, sous la forme matricielle: DRr ar O O DFr br DRV = O av O DF | + bv DRb O O ab DFb bb Dans un mode de realisation particulier de ['invention, les termes non diagonaux de la matrice a trots lignes et trots colonnes, peuvent etre non nuls pour tenir compte de l'effet du recouvrement spectral des absorptions, comme decrit par exemple dans le brevet US 5 796 874 precise. Selon les figures 6, 7, et 10, le procede de ['invention permet ensuite d'appliquer a chaque bande de pixels 44, 54, 64 de ['image numerique autre que la bande de pixels de moindre degradation 44M, 54M, 64M, et pour chacun des canaux de couleur formant ['image, une transformation lineaire permettant de 2s transformer les valeurs filtrees de densite optique 46, 56, 66 placees sur la courbe de reference superieure de densite optique 42, 52, 62 et les valeurs filtrees de densite optique 47, 57, 67 placees sur la courbe de reference inferieure de densite optique 43, 53, 63 de ladite chaque bande de pixels 44, 54, 64, respectivement en valeurs superieures 46M, 56M, 66M et inferieures 47M, 57M, 67M de densite optique de la bande de pixels de moindre degradation 44M, 54M, 64M. Le procede de ['invention permet done d'obtenir automatiquement une image numerique 12 dont les couleurs vent restaurees, c'est-a-dire dont les couleurs ont ete transformees pour correspondre a celles de ['image argentique d'origine. L'operateur peut par exemple, apres avoir applique le procede de ['invention sur une image alteree, afficher sur l'ecran 11 ['image numerique 12 dont les couleurs vent restaurees; il peut egalement restituer ensuite cette image restauree, par exemple sur du papier en utilisant une imprimante laser ou a jet d'encre, ou bien sur un film photographique. L'inventi on a ete decrite en fai s ant reference en p articulier a des o modes de realisation preferes de celle-ci; mais il est evident que des variantes et modifications peuvent etre apportees a l'interieur de la portee des revendications. RE VEND IC ATI ON S 1. Procede automatique de restauration de la couleur d'une image argentique (3) dont la couleur est alteree par rapport a la couleur d'origine de ladite image argentique (3), ['alteration de la couleur se produisant non s uniformement, spatialement dans le plan de ['image, et suivant une direction parallele a un axe principal (19) de ['image, ledit procede comprenant les etapes suivantes: a) numeriser ['image argentique alteree en utilisant une echelle de mesure des densites optiques; b) filtrer ['ensemble des pixels de ['image numerisee; c) diviser ['image numerique (12) en bandes de pixels (18) disposees perpendiculairement a la direction de l'axe principal (19) de ['image suivant lequel se produit ['alteration de couleur, la reunion de ['ensemble des pixels de chaque bande representant ltensemble des pixels de ['image numerique; d) calculer, pour chaque bande de pixels, et pour chacun des canaux de couleur (40), (50), (60) formant ['image, une distribution des densites optiques des pixels formant ladite bande; e) calculer, pour chacune desdites distributions de l'etape d), une valeur de reference brute superieure de densite optique (36), et une valeur de reference brute inferieure de densite optique (37); f) calculer, a partir de ['ensemble des valeurs de reference brutes superieures et inferieures de l'etape e), respectivement des valeurs de reference filtrees superieures de densite optique (46) et des valeurs de 2s reference filtrees inferieures de densite optique (47) correspondantes; g) determiner, a partir des valeurs de reference filtrees superieures et inferieures de densite optique de l'etape f), respectivement des courbes de reference superieures de densite optique (42) et des courbes de reference inferieures de densite optique (43), lesdites courbes de reference etant representatives du profil de degradation de couleur de ['image; h) determiner, pour chacun des canaux de couleur formant ['image, et a partir de la valeur maximale de densite optique (46M), (56M), (66M) de la courbe de reference superieure de densite optique (42), la bande de pixels de moindre degradation (44M), (54M), (64M); i) appliquer a chaque bande de pixels (44), (54), (64) de ['image numerique autre que la bande de pixels de moindre degradation (44M), (54M), (64M), et pour chacun des canaux de couleur formant ['image, une transformation lineaire permettant de transformer les valeurs filtrees de densite optique (46), (56), (66), (47) , (57), (67) placees sur la courbe de o reference superieure (42), (52), (62) et sur la courbe de reference inferieure (43), (53), (63) de densite optique de ladite chaque bande de pixels (44), (54), (64), respectivement en valeurs superieure (46M), (56M), (66M) et inferieure (47M), (57M), (67M) de densite optique de la bande de pixels de moindre degradation. 2. Procede selon la revendication 1, dans lequel le filtrage initial de ['image numerique alteree est opere en utilisant des operateurs de morphologic mathematique du type ouverture d'ordre (n) suivie de fermeture d'ordre (n), (n) etant un nombre entier. 3. Procede selon l'une quelconque des revendications 1 ou 2, dans lequel les valeurs de reference filtrees superieures de densite optique (46), (56), (66), et les valeurs de reference filtrees inferieures de densite optique (47), (57), (67), vent calculees a partir de l'utilisation d'un filtrage de ['image (12) par la methode du suivi des maximum dans chaque bande de pixels (18). 4. Procede selon l'une quelconque des revendications 1 a 3, dans loquel la division de ['image numerique (12) en bandes de pixels (18) est operee suivant des bandes paralleles adjacentes ntayant aucun pixel commun entre elles. 2841 357 5. Procede selon l'une quelconque des revendications 1 a 3, dans lequel la division de ['image numerique (12) en bandes de pixels (18) est operee suivant des bandes paralleles se recouvrant partiellement deux a deux suivant au moins une rangee de pixels. 6. Procede selon l'une quelconque des revendications precedentes, dans lequel les valeurs de reference brutes superieures (36) et inferieures (37) de densite optique correspondent respectivement aux courbes de la distribution de o densite optique (31), (32) placees au voisinage des courbes de densite optique maximum (30) et minimum (34) de ladite distribution. .. 7. Procede selon l'une quelconque des revendications 1 a 5, dans lequel les valeurs de reference brutes superieures (39) et inferieures (38) de densite optique correspondent respectivement aux courbes de la distribution de densite optique maximum (30) et minimum (34) de ladite distribution. 8. Procede selon l'une quelconque des revendications precedentes, dans lequel la bande de moindre degradation (44M), (54M), (64M) est une bande de o pixels de ['image comprenant le pixel de valeur maximale de densite optique (46M), (56M), (66M) de la courbe de reference superieure de densite optique (42), (52),(62). 9. Procede selon l'une quelconque des revendications 1 a 7, dans lequel la bande de moindre degradation est la bande de pixels comprenant une senle rangee de pixels, ladite bande passant par le centre (20) de ['image (12). 10. Procede selon l'une quelconque des revendications precedentes, dans lequel les plans des canaux de couleur formant ['image vent le rouge, le vert et le