Brevet FR2761552A1 — publié le 19981002

La présente invention concerne, d'une manière générale, un procédé et un dispositif pour le traitement de données, notamment la représentation, la génération, la restitution, la régénération ou le traitement itératif de données. o1 Plus particulièrement, la présente invention vise la transmission (c'est-à-dire la représentation à distance) et/ou le stockage (c'est-à-dire la représentation dans une mémoire) de données, avec ou sans compression de celles-ci. Ces données peuvent avantageusement être constituées par les données numériques ou analogiques d'une image ou séquence d'images, et/ou les données numériques ou analogiques relatives à du son (musique, parole, etc.),. et/ou les données numériques ou analogiques relatives à un signal mono ou multidimensionnel quelconque. Plus généralement, la présente invention vise la représentation de données mono ou multidimensionnelles quelconques. Avant d'exposer les objectifs et moyens de l'invention, on se propose ci-après de donner les définitions suivantes: "Ensemble de données": ensemble de données quelconques représentant des grandeurs physiques (le plus souvent des niveaux de tensions) pouvant représenter elles-mêmes d'autres grandeurs physiquement perceptibles ou mesurables. Dans des applications privilégiées, les ensembles de données concernés sont des images, du son ou des signaux mono ou multidimensionnels, etc... Dans l'application concernant le traitement d'images, on visera parfois l"'image" en lieu et place de l"'ensemble de données" relatives à l'image. Ce qui sera exposé ci-après à propos des "images" ou "sous- images" est respectivement applicable aux "ensembles de données" ou "sous- 3o ensembles de données" et vice versa. "Représentation de données": tout "traitement" d'un ensemble de données d'un type donné, résultant dans la transformation parfaite ou imparfaite dudit ensemble de données, dans un autre type. Exemple: les données pourront être constituées par les 256 niveaux de gris des pixels d'une image, leur représentation consistant en une série de mots binaires aptes à être stockés ou transmis; réciproquement, la représentation des données constituées par lesdits mots binaires pourra consister en leur transformation pour retrouver des données du même type que les données initiales. "Représentation primaire": par convention, tout traitement o10 résultant dans la transformation des données d'un premier type vers des données d'un second type. On parlera dans ce cas d'un "traitement primaire" des données "Représentation secondaire": par convention, la transformation de données du second type résultant d'un traitement primaire. On parlera dans ce cas de "traitement secondaire". "Restitution de données": cas particulier d'une représentation secondaire dans lequel les données du second type sont transformées dans des données du premier type. Cette restitution peut être parfaite ou imparfaite. Un "espace métrique" est un ensemble d'éléments muni d'une fonctionnelle distance qui est symétrique, positive et satisfait l'inégalité triangulaire. Cet espace est "complet" lorsqu'il contient tous les points limites de toutes les suites convergentes. Une "application Lipschitzienne": une application qui transforme les points d'un espace métrique dans le même espace et pour laquelle I'ensemble des rapports de la distance de deux éléments transformés par ladite application, sur la distance des deux dits éléments, est borné. Une "application contractive" ou "contraction": une application Lipschitzienne pour laquelle le plus petit des majorants (facteur de contraction) dudit ensemble de rapports est inférieur a l'unité. Toutefois, au sens de la présente invention, sont contractives toutes les applications convergentes, c'est-à-dire possédant un point fixe (permettant ainsi la mise en oeuvre d'approximations successives) dans un sous-ensemble de l'espace métrique. Une "similarité" est une application Lipschitzienne pour laquelle le rapport de la distance de deux éléments transformés sur la distance desdits éléments est une quantité fixe. Une similarité linéaire est une "similitude". Une "similarité contractive" est une application Lipschitzienne pour laquelle le rapport de la distance de deux éléments transformés sur la distance desdits éléments est une quantité fixe strictement inférieure à l'unité. Le "point fixe" d'une contraction d'un espace métrique complet l0 dans le même espace (ou d'un sous-ensemble de cet espace dans lui-même) est l'unique élément qui est laissé invariant par ladite application contractive. La "construction" d'une application contractive sur un ensemble de données consiste à constituer une famille de contractions susceptibles de transformer les données et à sélectionner les paramètres de l'une desdites contractions afin de satisfaire une ou plusieurs conditions prédéterminées. La "méthode des approximations successives" permet d'approcher de façon itérative aussi près que l'on le souhaite, le point fixe d'une contraction. En partant d'un élément arbitraire, on applique sur ce dernier ladite contraction. On applique ensuite la même contraction sur la transformée préalablement obtenue. En réitérant ce processus on s'approche successivement et inéluctablement du point fixe de la contraction. "Une meilleure approximation" d'un élément dans un espace métrique est un point d'un sous-ensemble de candidats, qui sont eux-mêmes des points dudit espace, qui minimise la distance audit élément. "Une bonne approximation" d'un élément dans un espace métrique est un point d'un sous-ensemble dudit espace qui est proche avec une erreur prédéterminée d'une meilleure approximation prédéterminée. On connaît, dans l'état de la technique, divers procédés de représentation d'image, avec compression, utilisant la technique dite "fractale". Par le document US-A-5.065.447, on connaît un procédé et un dispositif permettant le traitement o le stockage d'image avec compression des données relatives à l'image initiale. Dans ce procédé, les données relatives à l'image initiale sont divisées, c'est-à-dire que l'image est découpée en une pluralité de sous images élémentaires (appelées "domain blocks"). Ce procédé consiste ensuite à générer un dictionnaire ordonné d'un ensemble de sous images de référence (appelées "range blocks") constituées à partir de portions d'image de taille prédéterminée ayant subi un certain nombre d'opérations prédéterminées telles que rotation et retournement selon divers axes de symétrie. Ensuite, pour chaque sousimage élémentaire, on procède à une comparaison avec l'ensemble des sous images de référence du dictionnaire et o10 on sélectionne dans le dictionnaire la sous image de référence la plus similaire de la sous image élémentaire considérée. Enfin, le procédé consiste à traiter ou à stocker des informations paramétriques relatives aux adresses des sous images de référence sélectionnées dans le dictionnaire, pour représenter chacune des sous images élémentaires d'origine. C'est par cet ensemble d'opérations que le procédé décrit dans ce document permet d'obtenir une première représentation de l'image avec compression des données. Pour réaliser, à partir desdites informations paramétriques, la restitution de l'image initiale, le procédé et le dispositif décrits dans ce document, effectuent les opérations suivantes: à partir d'une image arbitraire de départ, on procède à des étapes similaires à celles ci-dessus, c'est-à-dire que l'on partitionne l'image arbitraire de départ et que l'on constitue un dictionnaire à partir des sous-images élémentaires de celle-ci. Toutefois, le dictionnaire n'est constitué que de façon partielle en n'effectuant pour chaque sous image que l'opération prédéterminée correspondant a l'adresse concernée du dictionnaire pour la position de la sous-image considérée. Les données ainsi obtenues sont utilisées pour réitérer ce processus jusqu'à ce que la différence entre deux images consécutivement obtenues soit inférieure à un seuil prédéterminé. L'image obtenue en dernier lieu est alors une approximation de l'image initiale. Cette méthode, intéressante au plan théorique, présente pour inconvénient majeur qu'elle est parfois difficilement praticable au plan industriel avec les moyens connus à ce jour. En effet, chaque image doit donner lieu à une longue analyse avec constitution d'un dictionnaire particulièrement important, les sous images élémentaires de l'image à stocker ou à transmettre étant toutes comparées à chacune des sous images de référence présentes dans le dictionnaire. Les inventeurs, qui ont procédé à des simulations, ont ainsi constaté que pour certaines images de taille et de résolution moyennes (512 x 512 pixels à 3 x 256 niveaux de couleur) le temps de traitement pour la 1o compression était de l'ordre de 1000 à 2000 secondes sur une station de travail, soit de l'ordre d'une demi-heure. Un tel temps de traitement est évidemment prohibitif pour quasiment toute application industrielle. On connaît par le document WO-A-93/17519, un autre procédé et dispositif permettant le traitement d'images avec compression des données relatives à l'image initiale. Dans ce procédé les données relatives à l'image initiale sont également partitionnées, c'est-à-dire que l'image est divisée en une pluralité de sous images élémentaires (appelées "domains"). Dans ce procédé on détermine ensuite un ensemble de transformées fractales liant les sous images à l'ensemble des données initiales, c'est-à-dire que l'on recherche dans l'image initiale une portion d'image semblable à l'image élémentaire considérée, de façon à minimiser l'erreur entre les données initiales et la donnée approximée ainsi obtenue. Les coefficients de l'ensemble de transformations fractales ainsi réalisées constituent la première représentation de l'image initiale. Pour réaliser la restitution de l'image, le procédé et le dispositif décrits dans ce document procèdent à la recherche des points fixes des transformations fractales élémentaires dont les coefficients ont été reçus sans toutefois expliquer comment cette recherche est effectuée. L'image initiale est ainsi restituée en procédant à une réunion des sous-images restituées. Cette méthode est également intéressante au plan théorique. Toutefois, les inventeurs ont constaté des dégradations importantes de la qualité de l'image restituée. Les inventeurs ont cherché à développer un procédé de traitement de données du même genre que les procédés de l'état de la technique analysé ci-dessus, mais présentant de meilleures performances tant au plan de la qualité des données restituées qu'au plan de la vitesse de traitement. Au cours de leur recherche, les inventeurs ont découvert que la l0 raison pour laquelle les procédés décrits dans US-A-5.065.447 et WO- A- 93/17519 ont les défauts susmentionnés tient au fait que les deux procédés font intervenir des transformations géométriques de portions d'images constituées essentiellement par des similarités contractives (construction du dictionnaire dans US-A-5.065.447 et transformation fractale dans WO-A- 93/17519). Ils ont en fait pu mettre en évidence qu'il s'agit dans les deux cas de la mise en oeuvre d'une idée largement partagée par l'homme du métier concerné par la compression fractale au moment du dépôt de la présente demande, idée selon laquelle les images seraient, quelle que soit leur origine, constituées essentiellement d'éléments ou de parties présentant entre-eux d'importants degrés de ressemblance. Or en fait, les performances de WO-A-93/17519 sont très faibles sur le plan de la qualité de l'image restituée car les inventeurs ont découvert que les capacités de restitution de similitudes contractives sont très limitées contrairement à la croyance desdits hommes de l'art à la date du dépôt de la présente demande. De même les performances de US-A-5.065.447 sont très limitées sur le plan du temps de traitement. Comme les capacités de restitution des similitudes contractives sont intrinsèquement limitées, I'auteur de US-A- 5.065.447 est amené à la réalisation d'un important dictionnaire qui est construit sur la base de nombreuses transformations de sous-images dans le procédé selon US-A-5.065.447: il s'agit en effet de multiplier artificiellement les éléments de ressemblance, ce qui permet certes d'améliorer sensiblement la qualité de l'image restituée, mais aux dépens du temps de calcul. Les inventeurs ont donc cherché à réaliser un procédé et un dispositif pour le traitement de données, notamment d'images, permettant d'obtenir de meilleures performances que l'état de la technique précité tant sur le plan de la qualité des données ou de l'image restituée que sur celui de la rapidité de calcul, et susceptible d'aboutir à une compression significative des données à transmettre ou à stocker. Dans le cas o une reconstruction d'une image avec un nombre de données différent de celui de l'image ayant donné lieu à la compression, est nécessaire, il est connu: - soit d'effectuer les trois opérations successives suivantes: décompression d'image, duplication de pixels pour augmenter la résolution de l'image puis lissage des valeurs de chaque composante, - soit d'effectuer une décompression d'image puis une interpolation pour ajouter ou retirer des pixels. La valeur de chaque composante des pixels de l'image finale est une interpolation de la valeur de la même composante de pixels de l'image initiale. Dans les deux cas, le traitement est compliqué puisqu'il possède plusieurs étapes et sa complexité varie avec la méthode d'interpolation choisie et avec les dimensions de l'image reconstruite. Les inventeurs ont donc cherché une méthode de décompression d'image qui soit simple et qui le reste quelles que soient les dimensions de l'image reconstruite. 2 5 Les procédés et dispositifs selon l'invention résolvent les problèmes susmentionnes. D'une manière générale, la présente invention vise tout d'abord un procédé de traitement primaire d'un premier ensemble de données représentatives de grandeurs physiques, caractérisé en ce qu'il comporte une étape de "construction" d'une application contractive globale dite "de premier type", pour ledit premier ensemble de données, dont le point fixe constitue une approximation de tout ou partie de cet ensemble, en mettant en oeuvre, pour ce faire, au moins une application appartenant au groupe d'applications constitué des applications multidimensionnelles et des applications non linéaires et en déterminant des paramètres de l'application contractive globale de façon à permettre la mise en oeuvre, sur un deuxième ensemble de données comportant un nombre de données différent du nombre de données du premier ensemble, d'une méthode d'approximations successives convergeant vers le point fixe d'une application contractive globale dite "de second type"; l'ensemble des paramètres ainsi déterminés constituant une représentation primaire dudit ensemble de données. D'une manière générale, la présente invention vise ensuite un dispositif de traitement primaire d'un ensemble de données représentatives de grandeurs physiques, caractérisé en ce qu'il comporte: - des moyens d'entrée de données à traiter; - des moyens de sortie de paramètres; - des premiers moyens logiques, de transformation des données à traiter en paramètres; lesdits premiers moyens logiques comportant: - des moyens de "construction" adaptés: à construire pour chaque ensemble de données à traiter, une application contractive globale, dite "de premier type" dont le point fixe constitue une approximation de tout ou partie de cet ensemble de données, en mettant en oeuvre, pour ce faire, au moins une application appartenant au groupe d'applications constitué des applications multidimensionnelles et des applications non linéaires, de façon: à permettre la mise en oeuvre d'une méthode d'approximations successives convergeant vers le point fixe de ladite application contractive globale, sur un deuxième ensemble de données comportant un nombre de données différent du nombre de données du premier ensemble de données; à délivrer sur les moyens de sortie les paramètres ainsi déterminés, ceux-ci constituant une représentation primaire dudit ensemble de données. Grâce à ces dispositions, la présente invention atteint ces objectifs. Il y a lieu de noter tout d'abord qu'il a été découvert que, d'une manière générale, conformément à l'invention, pour traiter, et notamment représenter, un ensemble de données initiales représentatives de grandeurs lo physiques, notamment avec compression desdites données initiales, il suffit de construire, par des moyens appropriés, une application contractive globale, agissant sur lesdites données. Une autre caractéristique importante de l'invention réside dans le fait que pour construire ladite application contractive globale il y a lieu de mettre en oeuvre au moins une application appartenant au groupe des applications constitué des applications multidimensionnelles et des applications non linéaires. Ce faisant, on exclut tout recours à des transformations mathématiques appartenant à la classe des similarités contractives, ce qui permet d'éviter les inconvénients relevés plus haut. Toutefois, pour construire ladite application contractive globale agissant sur les données, I'homme du métier a à sa disposition un très large éventail d'outils mathématiques dans lequel il peut puiser. Il peut par conséquent choisir parmi ces outils ceux qui sont le plus favorables aux performances, par exemple en conjonction avec le type de données à représenter ou l'application particulière visée. Une telle possibilité n'existe pas dans l'état de la technique analysé plus haut et s'avère très avantageuse. Ainsi, par exemple, dans le cas d'un traitement local de l'image, par partition de cette dernière en sous- images élémentaires, I'homme du métier pourra choisir, pour le traitement local, des applications qui ne soient pas forcément contractives, ce qui lui donne un important choix dans les possibilités de traitement local de l'image et lui laisse ainsi une grande liberté. En particulier, il pourra choisir au niveau local des applications Lipschitziennes qui lui permettront de faciliter la mise en oeuvre du caractère contractif de l'application globale. Par ailleurs, la mise en oeuvre d'une méthode d'approximations successives, laisse à l'homme du métier, la liberté de choisir l'écart final entre ledit point fixe et les données ainsi restituées ou générées. Enfin, les paramètres étant indépendants de la dimension du sous-ensemble élémentaire auquel ils sont associés, une reconstruction peut être effectuée avec la même simplicité et avec la même efficacité quelles que soient les dimensions de l'image reconstruite. En évitant la constitution d'un dictionnaire et la longue analyse décrite dans le document US-A-5.065.447, et grâce aux caractéristiques de la présente invention, le traitement d'une image est notablement accéléré par rapport au document US précité. En mettant en oeuvre leur mode préféré de réalisation, les inventeurs ont pu constater que pour une image de 512 x 512 pixels à 3 x 256 niveaux de couleur, le temps de traitement, était de l'ordre de moins de 100 secondes sur une station de travail; ils ont en outre tout lieu de croire que des perfectionnements apportés au procédé et dispositif succinctement définis plus haut, permettront d'obtenir des performances en matière de vitesse. On notera en outre que comme les applications élémentaires construites appartiennent au groupe d'applications constitué des applications multidimensionnelles et des applications non linéaires, on évite les applications appartenant à la classe des similitudes, avec les avantages observés plus haut. Dans un autre de ses aspects, la présente invention vise un procédé de traitement secondaire d'un ensemble de paramètres résultant d'un traitement primaire d'un premier ensemble de données, procédé de traitement secondaire caractérisé en ce qu'il comporte une étape de "calcul itératif" au cours de laquelle à partir dudit ensemble de paramètres, on met en oeuvre, sur un deuxième ensemble de données possédant un nombre de données différent du premier ensemble de données, une méthode d'approximations successives convergeant vers le point fixe d'une application contractive globale dite "de second type", en mettant en oeuvre, pour ce faire, au moins une application appartenant au groupe d'applications constitué des applications multidimensionnelles et des applications non linéaires, ce point fixe constituant une représentation secondaire dudit premier ensemble de données. Dans un autre de ses aspects, la présente invention vise un dispositif de traitement secondaire de données convenant pour le traitement secondaire d'un ensemble de paramètres résultant du traitement primaire d'un premier ensemble de données initiales, caractérisé en ce qu'il comporte: - des moyens d'entrée des paramètres; - des moyens de sortie de données; - des seconds moyens logiques de transformation des paramètres en données; lesdits seconds moyens logiques comportant: - des moyens de "calcul itératif' adaptés à mettre en oeuvre, à partir dudit ensemble de paramètres et d'un deuxième ensemble de données possédant un nombre de données différent du nombre de données du premier ensemble de données, une méthode d'approximations successives convergeant vers le point fixe d'une application contractive globale dite de "second type", ce dernier point fixe constituant une représentation secondaire dudit ensemble des données initiales, et à délivrer les données ainsi calculées sur les moyens de sortie. Ces aspects de l'invention présentent les mêmes avantages que ceux exposés ci-dessus. Ils ne sont donc pas rappelés ici. Dans l'hypothèse mise en oeuvre selon la présente invention, o I'application contractive globale mise en oeuvre lors du traitement primaire et l'application contractive globale mise en oeuvre lors du traitement secondaire sont de types distincts, en particulier parce que leur nombre de données est différent, la représentation secondaire ne constitue plus une reconstruction des données initiales mais constitue en réalité une transformation de celles-ci. C'est aussi le cas par exemple lorsqu'une image initiale en couleur a fait l'objet d'un traitement primaire et qu'elle est ensuite représentée par une image possédant le même nombre de données, mais en niveau de gris: dans les deux cas il s'agit bien de la même image, toutefois les données relatives à l'image obtenue à l'issue du traitement secondaire ne sont plus du même type que les données initiales. Grâce à ces dispositions, la méthode d'approximations successives utilisée pour le traitement secondaire est mise en oeuvre de façon particulièrement simple à partir des paramètres obtenus à l'issue du traitement primaire, parce que ces paramètres sont indépendants des dimensions du sous-ensemble élémentaire auquel ils sont associés. Les inventeurs ont pu constater, dans le domaine de traitement d'images avec compression, de bonnes performances tant sur le plan de la qualité de l'image restituée que sur celui de la vitesse de traitement rappelée supra. Un aspect de l'invention combine, au sein d'un même procédé de traitement, une phase de traitement primaire et une phase de traitement secondaire, sur des ensembles de données possédant un nombre de données différent. Il s'agit dans ce cas d'un traitement dit "complexe". L'invention vise aussi une imprimante, un photocopieur, un appareil photographique et un dispositif de télécopie comportant un dispositif de traitement tel que succinctement exposé ci-dessus. L'invention vise encore une mémoire qui conserve: - des paramètres et des informations de partition issues d'un traitement primaire tel que succinctement exposé ci-dessus et - des instructions pour un logiciel permettant la mise en oeuvre d'un procédé de traitement secondaire tel que succinctement exposé ci-dessus. Les avantages de ces dispositifs et mémoire étant identiques à ceux qui sont exposés ci-dessus, ils ne sont pas rappelés ici. Les caractéristiques, avantages et autres objets de l'invention ressortiront d'ailleurs lors de la description des dessins annexés sur lesquels: - la figure 1 illustre de façon schématique une partition d'une image de 512 x 512 pixels en sous ensembles élémentaires, constitués de carrés de 32 x 32 pixels, - la figure 2 illustre de façon schématique une division d'une image de 512 x 512 pixels en un recouvrement par ardoises, constituées par des carrés de 32 x 32 pixels, décalés horizontalement et verticalement de 16 pixels, - la figure 3 illustre de façon schématique l'agencement d'un sous- ensemble élémentaire K, et de ses quatre sous-ensembles de travail L., L2. L3 et L4, qui, dans un premier mode de réalisation de l'invention, sont de taille identique à celle du sous-ensemble élémentaire K, et décalés dans les deux directions, par rapport audit sous-ensemble élémentaire, de manière prédéterminée, - la figure 3A est similaire à la figure 3 et illustre de façon schématique, un deuxième mode, préféré, de réalisation de l'invention, dans lequel l'agencement d'un sous-ensemble élémentaire K, et de ses quatre sous- ensembles de travail L1, L2, L3 et L4, dont chaque dimension est double de celle du sous-ensemble élémentaire Ki associé, et qui est décalée dans les deux directions, par rapport audit sous-ensemble élémentaire, de manière prédéterminée, - la figure 4 est un schéma bloc simplifié d'un mode particulier de réalisation du dispositif de traitement primaire d'un ensemble de données, - la figure 5 est un schéma de principe d'un dispositif mettant en oeuvre le dispositif de la figure 4, - la figure 6 est un schéma bloc simplifié d'un mode particulier de réalisation du dispositif de traitement secondaire d'un ensemble de données, - la figure 7 est un schéma de principe d'un dispositif mettant en oeuvre le dispositif de la figure 6, - la figure 8 est un schéma de principe d'un dispositif de traitement primaire et de traitement secondaire, - la figure 9 est un schéma de principe d'un autre dispositif de traitement primaire et de traitement secondaire, - la figure 9A est un schéma bloc simplifié d'un mode particulier de réalisation du dispositif de traitement de la figure 9, - la figure 10 est un organigramme simplifié de compression d'une image couleur de 512 x 512 pixels mis en oeuvre dans le dispositif des figures 4, 9 et9A, - les figures 11, 12, 12A, 12B, 13, 14 et 15 sont des organigrammes illustrant certaines opérations de l'organigramme de la figure , - la figure 16 est un organigramme simplifié de décompression lo d'une image couleur mis en oeuvre dans le dispositif des figures 6, 9 et 9A, - la figure 16A est un organigramme illustrant une opération de l'organigramme de la figure 16 - les figures 17, 17A et 17B sont des organigrammes illustrant une opération de l'organigramme de la figure 16, - la figure 18 est un organigramme d'un programme de configuration mis en oeuvre dans le dispositif de la figure 9A. - les figures 19 et 20 sont des organigrammes illustrant certainesopérations de l'organigramme de la figure 10 et d'une de ses variantes, et - la figure 21A est une représentation d'une succession de partitions de sous-ensembles élémentaires, - la figure 21B est une représentation de structure d'informations de partition et de paramètres représentant cette succession de partitions et - la figure 21C est une représentation de l'organisation des paramètres et des informations de partition résultant du traitement primaire. La présente description est divisée en cinq parties: - la première partie a trait à la description générale du procédé. Cette partie est divisée en trois sections: À la première section a trait à la théorie générale à la base de l'invention. Il y est exposé les concepts mathématiques qui sous-tendent, d'une manière générale, l'invention telle qu'elle peut être comprise dans ses termes les plus généraux; la seconde section, expose, notamment à l'appui des figures 1 à 3, certaines caractéristiques particulières du premier mode de mise en oeuvre du procédé permettant le traitement primaire des données d'une image de 512 x 512 pixels aboutissant à la compression de celle-ci. Cette section explique plus particulièrement la structure de l'application multidimensionnelle choisie dans ce mode particulier de mise en oeuvre et les moyens permettant de calculer les coefficients de cette application. la troisième section expose les caractéristiques particulières du premier mode de mise en oeuvre du procédé de représentation secondaire des données de l'image, et notamment la restitution des données de celle-ci à partir des données compressées obtenues à l'issue d'un traitement primaire. Cette section expose plus particulièrement la méthode d'approximations successives mise en oeuvre ici pour ce faire. - la seconde partie a trait à la description d'un premier mode de réalisation d'un dispositif de traitement et elle est divisée en six sections: la première section a trait à la description d'un dispositif de traitement primaire d'un ensemble de données, ò la seconde section a trait à un dispositif de traitement secondaire d'un ensemble de paramètres résultant du traitement primaire de données 2 0 initiales. Un tel dispositif permet, par exemple, de restituer l'image compressée par le dispositif de traitement primaire, À la troisième section a trait à la description d'un dispositif incorporant à la fois un dispositif de traitement primaire et un dispositif de traitement secondaire, À la quatrième section a trait à la description d'un organigramme particulier de compression d'une image en couleurs de 512 x 512 pixels mis en oeuvre dans les moyens logiques que comporte le dispositif de traitement primaire, À la cinquième section a trait à la description de deux modes de réalisation d'un organigramme de décompression mis en oeuvre dans les moyens logiques que comporte le dispositif de traitement secondaire, et la sixième section a trait à la description de l'organigramme du programme de configuration du dispositif incorporant un dispositif de traitement primaire et un dispositif de traitement secondaire. - la troisième partie a trait à la description du mode de réalisation préféré, mettant en oeuvre des sous-ensembles de travail présentant une taille différente de celle du sous-ensemble élémentaire correspondant. Dans cette partie, il est uniquement décrit les modifications qui sont apportées au mode de réalisation décrit dans la deuxième partie; - dans la quatrième partie il est décrit une variante de réalisation dont l'une met en oeuvre une application contractive globale composée d'applications non linéaires. Dans cette quatrième partie seuls sont exposés les éléments modifiés par rapport à ceux du mode de réalisation décrit dans la seconde partie; et - dans la cinquième partie, il est décrit une variante de réalisation mettant en oeuvre un critère de qualité mesuré à l'aide de l'application contractive globale. 1. Description générale du procédé On va maintenant décrire à l'appui des figures 1 à 3, le procédé tel qu'il est mis en oeuvre dans le premier mode de réalisation de la présente invention. Cette description va être faite à l'appui des figures précitées et d'un certain nombre de concepts mathématiques appartenant d'une manière générale à la théorie de l'analyse et plus particulièrement à celles de la topologie et de l'analyse fonctionnelle. Dans ce mode particulier de réalisation, on se propose de compresser les données relatives à une image, pour aboutir à une "représentation primaire" de celle-ci, puis de restituer les données initiales, c'est-àdire de réaliser une "représentation secondaire" des données compressées pour obtenir un ensemble de données du même type que les données initiales avec un écart maximum prédéterminé entre les données ainsi restituées et les données initiales. Pour fixer les idées, il est ici considéré une image couleur de 512 x 512 pixels à trois composantes, I'une de luminance et deux de chrominance, possédant chacune 256 niveaux. En général, une telle image couleur est captée sous forme de trois composantes relatives aux trois couleurs fondamentales, rouge, verte et bleue, puis convertie, en mettant en oeuvre des matrices de conversion bien connues de l'homme du métier, en composantes de luminance et de chrominance. Cette image numérique nécessite 512 x 512 x 24 bits à stocker soit, sans compression de données, 6 Mbit. Afin de mesurer l'écart entre I'image initiale et l'image restituée, il est considéré la racine carrée de l'écart quadratique moyen, moyenné sur tous les pixels et tous les composants (noté classiquement RMSE). L'efficacité d'un algorithme de compression est exprimé par le lien entre le taux de compression obtenu pour un écart (RMSE) donné. Le mode particulier de réalisation décrit ici permet d'obtenir l'efficacité suivante Tableau 1 Nombre itérations à RMSE bits/pixel la décompression - 7 7 - 15 0.75 - 8 5 - 10 1.50 6 - 8 4 - 8 2.25 7-9 3 -7 3.00 7 - 10 2 - 5 4.50 2 0 Théorie générale à la base de l'invention La représentation d'un ensemble de données, représentant elles- mêmes des grandeurs physiques, telles que les données afférentes à l'image décrite ci-dessus, surtout lorsque l'on souhaite comprimer ces données à l'occasion de ladite représentation, peut s'analyser comme une approximation de l'ensemble des données. L'invention procède de la mise en évidence, par les inventeurs, du fait que la théorie du point fixe peut être mise en oeuvre pour réaliser une telle approximation. On considère que l'approximation d'un ensemble de données initiales est une opération consistant à rechercher dans une famille de candidats à la représentation, le (ou les) candidat(s) qui minimise(nt) un critère prédéterminé. Le plus souvent, le critère s'exprime par une distance entre deux 1o éléments d'un espace métrique, ce qui peut s'exprimer par l'équation suivante: d: (g,h) e E x E -> d(g, h) E + (1) o: - d est la fonction distance, - Eest un espace métrique, -R+ est l'ensemble des réel positifs ou zéro. - g, h sont des éléments de l'espace métrique E, c'est-à-dire de l'espace qui contient toutes les applications d'un domaine de définition donné Q c R' dans un domaine d'application Rd. Les données à traiter sont des applications de Q c R' dans R': g. (x, y) e 2 c RC --> g(x, y) E R' O Roq et Rd sont des espaces produits de R d'un nombre (n ou d) donné de dimensions. L'ensemble de données à traiter est un espace métrique E donné. On notera que par exemple dans le cadre d'une image couleur 512 x 512 à 3 x 256 niveaux de chacune des trois composantes, les données de l'image initiale sont constituées par un tableau de correspondance qui à chaque position (x,y) de pixel dans le domaine de définition [0,511] x [0,511] associe trois valeurs positives de niveaux des composantes de luminance, Y, de chrominance rouge, Cr, et de chrominance bleue, Cb, comprises entre 0 et 255. On observe que ces données de l'image initiale constituent en réalité une application f EE d'un domaine de l'espace R2 dans un domaine de l'espace R3 g:(x,y) e [O,511] x [O, 511]cR2 -> [0,255]3 c R3 Le problème de l'approximation d'un ensemble de données à traiter /f E consiste donc à rechercher dans une famille de candidats à la représentation G c E le ou les candidats g* dans G qui minimisent la distance (1). Ceci s'exprime par l'équation suivante: ?g' e G:d(f,g)= inf {(f,g)l g eG}. Lesdits candidats appartiennent à l'ensemble: PG (f)={ g*e GI d(f, g*)=inf (f,g)lgeG, (2) o: Ou: - P; est l'opérateur de projection métrique, - inf {...} est le plus grand des minorants de l'ensemble considéré '.... Dans la suite, on suppose que l'espace métrique est complet, ce qui garantit l'existence d'un point fixe unique pour chaque contraction de l'espace métrique dans lui-même. Les inventeurs ont mis en évidence que des candidats à la représentation peuvent être définis comme étant des points fixes relatifs à des contractions, une famille de contractions étant conjointement déterminée. Dans ces conditions, la recherche de la meilleure approximation revient à déterminer le point fixe relatif à la contraction minimisant ledit critère prédéterminé d'erreur dans cette famille de contractions. Ces conditions peuvent être exprimées ainsi ?a a e((X):d(f,a) _ t E E (4) o t est le point fixe de la contraction T (T (t) = t), CON(E) I'ensemble des contractions de E. Le problème d'approximation peut, dans ces conditions, être reformulé ainsi qu'il suit: il s'agit de rechercher la contraction dont le point fixe est la meilleure approximation des données initiales. On observera que ce faisant les inventeurs ont déplacé le problème à résoudre puisqu'il ne s'agit plus maintenant de chercher un candidat constituant en lui-même une meilleure approximation mais une contraction dont le point fixe constitue la meilleure approximation. Ceci peut être exprimé par l'inéquation: ?T* eX: d(f, (T*)) o X est un sous-ensemble de contractions de E (X c CON(E)). Les inventeurs ont ainsi mis en évidence que, conformément au procédé général découvert par eux, I'étape essentielle pour aboutir à une représentation de données à traiter, notamment dans le cadre d'une compression des données, consiste: - à constituer une famille de contractions susceptibles de transformer les données, et - à sélectionner les paramètres de l'une desdites contractions afin de satisfaire une ou plusieurs conditions prédéterminées, en l'espèce minimiser la distance entre les données initiales et le point fixe constituant la meilleure approximation desdites données initiales. Il s'agit là, conformément à l'invention, d'une étape de "construction" de l'application contractive globale. Ensuite, conformément à l'invention, les paramètres ainsi déterminés doivent permettre, pour la restitution de l'ensemble initial de données, la mise en oeuvre d'une méthode d'approximations successives convergeant vers le point fixe de ladite transformée contractive globale représentatif dudit ensemble de données initiales, ce qui permet pour les lo raisons exposées plus haut de restituer avec la meilleure approximation ledit ensemble de données à traiter. On observera que pour la mise en oeuvre du procédé général selon l'invention, et notamment dans l'étape de construction de l'application contractive globale, I'homme du métier a, pour constituer la famille de contractions susceptibles de transformer les données, un large éventail de choix d'applications mathématiques dans le groupe des applications constitué des applications multidimensionnelles et des applications non linéaires. Dans le mode de réalisation décrit ci-après, la contraction globale sera construite par morceaux, et pour chaque morceau il va être utilisé une application linéaire multidimensionnelle spécifique. Il pourrait utiliser d'autres applications, par exemple des applications non linéaires. Le procédé général, dont un mode de mise en oeuvre particulier va être décrit en détail ci-après, permet ainsi de très nombreuses possibilités de traitement des données, en vue de les représenter. On va maintenant décrire les moyens utilisés dans le mode particulier de réalisation pour construire l'application contractive globale. Les inventeurs ont mis en évidence le fait que pour rechercher une contraction dont le point fixe constitue la meilleure approximation des données initiales, il fallait tout d'abord chercher un candidat constituant en lui- même une bonne approximation, ce qui constitue une approche pratique pour 3o résoudre un tel problème. Les inventeurs ont mis en évidence que ladite approche pratique peut consister à trouver la contraction qui minimise la distance entre les données initiales et transformées, en prenant en compte un facteur de contraction estimé. Cette reformulation du problème à résoudre peut être exprimée par le problème d'optimisation principal suivant: ?T E X: d f =min)j = T e X (6) o con (T) est le facteur de contraction de l'application contractive T. Comme, en réalité, il est parfois difficile de trouver ledit facteur de contraction de l'application globale, les inventeurs ont trouvé que, dans ces hypothèses, on peut se contenter de résoudre un problème d'optimisation simplifié aux termes duquel il faut trouver une contraction globale qui minimise la distance entre les données initiales et leur transformée. On observera, ici, que le facteur 1/1-con(T) n'est plus pris en compte. Ce problème d'optimisation simplifié s'exprime comme suit: ? T e X:d(f, T(f))= min.{(f, T(f)" T e X, (7) Premier mode de traitement primaire des données de l'image. Selon une caractéristique importante d'un des aspects de l'invention telle que mise en oeuvre ici, la contraction globale est construite par morceaux. Pour ce faire, il est procédé à une étape de "partition" de l'ensemble des données relatives à l'image, au cours de laquelle on partitionne ledit ensemble en m "sous-ensembles élémentaires". Dans le mode de réalisation particulier qui concerne le traitement d'image, d'une manière générale, on partitionne l'image initiale, en fi sous- images élémentaires, elles aussi appelées "sous-ensembles élémentaires" pour des raisons de clarté de la description. On rappellera ici que d'une manière générale une "partition" de l'espace métrique Q, revient à définir l'ensemble T = {KiJKi c Qi = I.....in. o Ki,..., K,,, sont des sous-ensembles définis dans ledit espace métrique Q, devant respecter les quatre conditions suivantes: 1) les sous- ensembles K, doivent être fermés, ce qui s'exprime par K, = adh Ki, i = l,...,m. 2) aucun des sous-ensembles Ki ne peut avoir un intérieur vide, ce qui s'exprime par: int Ki = 0, i = I,....m. 3) leur union doit produire la fermeture de l'espace métrique Q considéré, ce qui s'exprime par: adhQ = U K, 1=1 4) les intérieurs de chaque sous-ensemble K, doivent être disjoints deux à deux, ce qui s'exprime par: int K, nr int K, = 0 i,. = 1.....m, i j. On notera ici également que la partition d'une donnée initiale d'un type donné (g e E) revient, d'une manière générale à définir une restriction g K, pour chaque sous-ensemble Ki de la partition T. Cette caractéristique permet, avantageusement, de faciliter le travail de construction de la contraction globale, en ramenant ce travail à celui de trouver des transformations mathématiques élémentaires ayant pour objet chacune de permettre la restitution du sous-ensemble élémentaire K/. En effet, par définition, les sous-ensembles de données comportent moins de données que l'ensemble initial, et par conséquent il est sensiblement plus facile de trouver des transformations mathématiques permettant de restituer ces données, étant précisé à nouveau que l'ensemble des transformations mathématiques doit, au niveau de l'ensemble de l'image, constituer une application contractive globale. Par conséquent, l'application contractive globale T est ainsi construite par morceau, chaque morceau étant, conformément à cet aspect de l'invention, une application Lipstchizienne. On observera ici d'ailleurs que ces applications Lipstchiziennes lo élémentaires ne sont pas nécessairement contractives, et qu'à ce niveau également, grâce à l'invention, l'homme de l'art a à sa disposition un grand éventail d'applications mathématiques lui permettant de traiter les données relatives aux sous-ensembles. On notera toutefois que dans certaines hypothèses, il n'y aura pas is lieu de procéder à cette partition, par exemple lorsque l'image donne lieu à une structure géométrique simple. Selon l'aspect de l'invention décrit ici, pour chacun des m sous- ensembles élémentaires (KI,....Km), on considère n sous-ensembles de travail (n = 4 dans ce mode de réalisation). Pour chacun des m sousensembles élémentaires réalisés au cours de la partition on construit l'application Lipstchizienne élémentaire considérée à partir de ces n sous-ensembles de travail. De manière à pouvoir considérer les n sousensembles de travail liés à chacun des m sous-ensembles élémentaires, on va réaliser une "division" de l'ensemble de données initiales. Cette division r = {L L, c Qj =.... de l'ensemble de données, qui est définie dans l'espace métrique Q, doit respecter les trois conditions suivantes: 1) les sous-ensembles L doivent être fermés, ce qui s'exprime par: L/ = adh L/, = 1...I 2) Aucun des sous-ensembles Li ne peut avoir un intérieur vide ce qui s'exprime: int L/#,j =- 1...l. 3) Leur union est incluse dans l'espace métrique Q, ce qui / s'exprime par: U Lj c Q j=l On observe à ce stade que la division ainsi réalisée englobe d'une part, les partitions, telles qu'elles ont été définies supra, mais englobe aussi, d'autre part, ce qu'on appellera un "recouvrement par ardoises" (notion connue o en anglais par "tiling") pour lesquels les sous-ensembles L, sont susceptibles de se recouvrir, ce qui s'exprime par l'équation int Lin int L., ó 0, i,j = 1,...,1,i ó j. Dans le premier mode de réalisation, I'image de 512 x 512 pixels est, selon une partition initiale, partitionnée en 256 sous-ensembles élémentaires carrés de 32 x 32 pixels (voir figure 1). Conformément à une habitude dans le domaine de traitement de l'image, les données sont organisées comme suit: - le point de référence de coordonnées (0,0) est constitué par le point en bas à gauche de l'image; 2 0 - les colonnes (abscisse x) sont numérotées de gauche à droite; - les lignes (ordonnée y) sont numérotées de bas en haut. Sur la figure 1 on a illustré de façon schématique en ordonnées le rang de certaines lignes de l'image et en abscisses le rang de certaines colonnes. Pour des raisons de clarté, l'échelle utilisée sur l'axe des abscisses et l'axe des ordonnées n'est pas uniforme. Sur la partie en bas à gauche de la figure 1 les 9 premiers carrés de 32 x 32 pixels sont illustrés, de façon à faire ressortir les rangs et colonnes constituant la frontière de ces derniers. Les carrés de 32 x 32 pixels constituent un exemple de partition du domaine Q = [0,511] x [0,511] o l'image considérée est définie. Dans le premier mode de réalisation, la division de l'image de 512 x 512 pixels en sous-ensembles de travail est illustré en figure 2. Elle est similaire à la partition de la figure 1, à ceci près que les sousensembles de travail de 32 x 32 pixels sont décalés de 16 pixels dans les deux directions, et dans chaque sens de déplacement, par rapport aux sous-ensembles élémentaires. Sur cette figure, qui illustre le coin inférieur gauche de l'image de 512 x 512 pixels, les sous-ensembles élémentaires résultant de la partition sont illustrés en traits gras tandis que les sous-ensembles de travail résultant de la division sont illustrés en traits pointillés et sont hachurés. En figure 2, n'ont pas été représentés les sous-ensembles de travail touchant les bords de l'image. De nombreuses solutions existent pour fournir à chaque sous-ensemble élémentaire qui n'est pas associé à des sous- ensembles de travail intégralement à l'intérieur de l'image initiale, des sous- ensembles de travail de substitution. On note ici, que dans toute la suite de la description, le choix * d'une telle solution est appliqué, aussi bien aux sous-ensembles de travail pour le traitement primaire, pour le traitement secondaire et pour le traitement complexe: les points des sous-ensembles de travail définis ci-dessus qui n'appartiennent pas à l'image, c'est-à-dire ceux dont l'une au moins des coordonnées est négative, se voient fictivement associer une valeur moyenne, par exemple 127 pour la composante de luminance et les éventuelles composantes de chrominance. La figure 3 illustre comment, dans ce premier mode de réalisation, chaque sous-ensemble de travail Lú est lié à un sous-ensemble élémentaire Ki donné. Le sous-ensemble élémentaire K, ayant pour coordonnées (x,,y,) est illustré en traits gras, tandis que les sous- ensembles de travail correspondants sont illustrés en traits pointillés, les autres sous-ensembles étant ébauchés en traits fins. On observe que, dans ce premier mode de réalisation, à chacun des sous-ensembles Ki sont liés quatre sous-ensembles de travail L, (numérotés L,, 30..., L4) décalés de 16 pixels dans les deux directions. Il s'agit en réalité des quatre sous-ensembles issus de la division qui recouvrent le sous- ensemble élémentaire considéré (issu de la partition). Les expressions ci-dessous montrent les coordonnées d'un sous-ensemble K, et des sous- ensembles de travail L/,.... L4. Ki: [xi, xi+ 31] x [yi, Yi + 31] LI: [xi-16, x i + 15] X [yi +16, yi + 47] L2: [xi+16, xi+ 47] X [yi + 16, yi + 47] L3: [xE-16, Xi + 15] X [y1 - 16, y, + 15] L4: [xi+16, xi+ 47] x [Yi- 16, yi + 15] On notera que, dans ce premier mode de réalisation, les sous- ensembles résultant de la division présentent les mêmes dimensions que ceux résultant de la partition; en outre, dans ce premier mode de réalisation, la division de l'image en sous-ensembles de travail ne constitue pas un recouvrement par ardoises mais une simple partition en raison du fait que les sous-ensembles de travail présentent une taille identique aux sous-ensembles et que le décalage, entre un sous-ensemble de travail et le sous-ensemble élémentaire correspondant, est égal à un demi côté des sous-ensembles carrés 2 0 considérés. Infra il sera décrit en particulier à l'appui de la figure 3A un deuxième mode, préféré, de réalisation dans lequel, conformément à un autre aspect de l'invention, les sous-ensembles de travail présentent une taille différente de celle du sous-ensemble élémentaire auquel ils correspondent, la division de l'image résultant ainsi en un recouvrement par ardoises. On va maintenant définir l'application globale utilisée dans le mode particulier de réalisation décrit à l'appui des figures 2 et 3. On rappelle tout d'abord que l'application contractive globale consiste à prendre un élément g de l'espace métrique E pour le transformer en un élément T(g) du même espace métrique E ce qui s'exprime par T: g E E-> T(g) E E (8) L'appartenance de g a l'espace E (g E E) peut être exprimée également de la façon suivante: g: (x,y) E Q c Rq -> g(x,y) E IF. (9) On observe en effet que l'équation (9) ne fait qu'expliciter le type o10 de données spéciales considérées: à chaque couple de valeurs (x.^) pris dans le domaine de définition (en l'espèce l'image plane Q cR2) correspond une valeur dans le domaine d'application Rd, c'est-à-dire trois valeurs parmi 256 niveaux des trois composantes, Y, Cr et Cb. On observera d'ailleurs que s'agissant du traitement d'images, le domaine de définition est R2, tandis que le domaine d'application est R3 puisque l'image couleur a trois composantes, I'une de luminance et deux de chrominance. Dans le mode particulier de réalisation des figures 2 et 3 décrit ci- après et conformément à un autre aspect de l'invention, chaque "morceau" de l'application contractive globale, c'est-à-dire une application élémentaire Lipschitzienne, est une application multilinéaire du type: 1' T(g)|K,:(x,y)Ki -> [T(g)](x,y) = a,-.g,, (x.y)+/b (10) /=l o: T(g)l Ki la restriction de T(g) au sous-ensemble élémentaire K, n > I la dimension de l'application multilinéaire a,... an les coefficients multiplicatifs b le facteur de translation Li,... L, les sous-ensembles de travail, à position prédéterminée par rapport à K, et lié à ce sous-ensemble. Ainsi dans l'exemple de la figure 3, les valeurs g (xy) des pixels de coordonnées (x,y), sont transformées comme suit: [T(g)](x,y) = a,. g(x16, y+16)+ a2. g(x+ 16, y+16)+ a3. g(x- 16, y-16)+ a4. g(x + 16, y-16)+h (1 Obis) o: g (x,y) est la valeur du pixel (x,y) dans l'image g (x, y) est un pixel du sous-ensemble élémentaire K,. 1 0 g(x - 16, y + 16) est la valeur du pixel correspondant dans le sous- ensemble de travail L, g(x + 16, y+ 16) est la valeur du pixel correspondant dans le sous- ensemble de travail L2 g(x - 16, y- 16) est la valeur du pixel correspondant dans le sous- ensemble de travail L3 g (x + 16, y- 16) est la valeur du pixel correspondant dans le sous- ensemble de travail L4 a.... a4 les coefficients multiplicatifs b le facteur de translation. On voit que les équations (10) et (lObis) ci-dessus sont de nature linéaire et multidimensionnelle c'est-à-dire multilinéaire, les coefficients a, pour.j = 1,..., n et b étant respectivement les coefficients multiplicatifs et le facteur de translation. Les coefficients (aI.... an, b) sont ici les "paramètres", au sens de I'invention, de l'application élémentaire relative au sous-ensemble Ki. L'ensemble des coefficients (ai.. an,b)K' relatifs aux rn sous-ensembles élémentaires (K,.... Km) constitue, dans ce mode de réalisation, "l'ensemble des paramètres" de l'application contractive globale, au sens de l'invention. On rappelle que dans le mode particulier de réalisation, n = 4, ce qui signifie que seulement quatre sous-ensembles de travail L,, j = I n tels que ceux définis à l'appui de la figure 3, sont liés à chaque sousensemble élémentaire Ki. Ainsi, à chaque couple de valeurs (x,y) d'un sous-ensemble élémentaire Ki déterminé, on associe une valeur calculée à partir de l'équation l0 multilinéaire (10) (voir l'exemple de la figure 3 et l'équation 10 bis). Ceci est fait pour l'ensemble des m sous-ensembles élémentaires (K,/..., K,,,). On va maintenant expliquer comment les coefficients de l'application multilinéaire aj, j = 1,..., n et b sont déterminés pour chacun des sous-ensembles élémentaires de travail Ki. Pour résoudre le problème d'optimisation simplifié (7) il faut, d'une manière générale, pour chaque sous-ensemble élémentaire K,. minimiser la distance entre la restriction au sous-ensemble élémentaire K, et satransformée par l'application contractive globale. En d'autres termes, il faut minimiser les distances: d (flKi, T(f) Ki), i= I,...m (11) Dans le mode particulier de réalisation, on met en oeuvre une méthode de minimisation des moindres carrés en présence de contraintes, o la quantité à minimiser est la distance définie dans l'équation (11) et s'exprime sous la forme: [d (fK, T(f)|K) = i (lf(xy) - [T(f)] (X))2 (12) (x,y) e Ki ai] o I En minimisant l'expression (12) sous les contraintes (13), on détermine chaque restriction de T(f) à Ki et en même temps les coefficients ac et h de l'équation (10). Les coefficients aj, b associés à chacun des sous-ensembles K constituent la représentation primaire de l'ensemble de données initiales. Le problème de minimisation selon les équations (10) et (12) sous les contraintes (13), afin de déterminer les coefficients (a.. a,,. h) pour chaque Ki: i = i,..., m sera détaillé ci-après. Pour chaque sous-ensemble élémentaire, nous considérons la fonction objectif 4 K: Rn+l -> R définie comme suit: K, (al,...anb)= A f(x, y) aj. /'I (X'y)+. (14) (x,y)K, L = I Le problème de minimisation sous contraintes revient donc à rechercher les (n+l) valeurs (al*,..... an, b*) e C c R." qui satisfont: d K, (al,...,a,,", b')= min# K, (al,....,a,,,b)l(a,..... a,.h) ' CI}, (15) o C =a,...a,=,..., n b)R- Rappelons tout d'abord que la résolution d'un problème de minimisation sans contraintes est bien connu de l'homme de l'art. Il s'agit dans ce cas de minimiser en général une fonction objectif du type: Y(cI, cm)= ó (x, y) - m-l ci.g(x,y)+c, r,(16) (x, y) K = o m est le nombre de paramètres de la fonction objectif et lorsque (c...., c,) parcourt tout l'espace Rm. Le vecteur (cl*. Cm..,*) recherché satisfait: YI,....cm)= min (....Cm cm) (cl.Cm)E Rm). Ce vecteur (c*...c.,*) est caractérisé par l'équation d' Euler: Cl *.,cm)=O, (17) o I'(ci*... cm*) est la dérivée de Y (au sens de Fréchet). C'est une application linéaire et continue de R' a R. 2 o Cette dernière équation est remplacée ici par les dérivées partielles de Y au point (cl.... cm*), par rapport a chaque variable. On obtient donc que (cl c...,Cm) est caractérisé par: c (Cl... Cm=0, E 5,"P(cî*,... cm*)=O. i. (18) IF ud J I.....Cm*=0 Tenant compte de la forme particulière de la fonctionnelle (16), les équations (18) sont équivalentes à: m ci Gji = Si; i = I..m, (19) j=l o G. g1 (xy).-gi(xj4 {J.mi (xy)EK Gmi Gim = Y gi(xy) I (X,y)EK Gmm = nombre de points (xy) dans Ki, Si = gj(xy) F(xy) i (x,y)EK Sm = Z F(xy). (x,y) eK Les équations (19) forment un système linéaire d'équations facile à résoudre par diverses méthodes bien connues dans les mathématiques. La résolution du problème de minimisation (15) de la fonctionnelle (14) sous les contraintes (13) comprendra les étapes suivantes: a) Résolution du problème de minimisation sans contraintes de la fonctionnelle (14) sur R+' par la méthode décrite supra par les équations (16) à (19). La solution sera notée (l,... a, b) e R" b) Si (a- ,..., à,b) e C, donc si -3 -13 - UVest l'ensemble des indices j e {i,...,n} pour lesquels [3 La solution du dernier problème (24) et (25) sera (c*... c,,). Finalement la solution recherchée du problème avec contraintes sera a, =c 1j À j = l...,m-1 o i E NV, (26) ak = À k E UV, (27) a, = - IE LV, (28) b = Cm. (29) si toutefois (26) respecte les contraintes. Dans le cas contraire, il faudra réitérer le processus. Un organigramme permettant la partition de l'image de 512 x 512 pixels, sa division et le calcul, pour chacun des sous-espaces élémentaires K, de coefficients aj et b est décrit infra, à l'appui des figures 10 à 15. On observe ici que, lorsque tous les sous-ensembles élémentaires possèdent les mêmes dimensions, un compromis doit être effectué. En effet, lorsque ces dimensions sont trop réduites, le taux de compression est limité puisque celui-ci est directement fonction du nombre de pixels qui sont représentés par les cinq paramètres associés à chaque sous- ensemble élémentaire. En revanche, lorsque les dimensions des sous- ensembles élémentaires sont importantes, la qualité de l'image reconstruite après un traitement secondaire de ces paramètres, peut être trop limitée pour une exploitation confortable de l'image reconstruite, les plus petits détails pouvant avoir été mal traités au cours du traitement primaire. Dans une image, il est courant que certaines zones présentent très peu de détail et soient pratiquement uniformes et que d'autres zones de l'image présentent une grande complexité (complexité pouvant être o10 représentée par la variance des valeurs de chaque composante des pixels de la zone considérée). L'invention propose, en complément de la mise en oeuvre de la théorie exposée ci-dessus, une méthode permettant d'utiliser des sous- ensembles élémentaires de grande dimensions dans les zones de l'image les 1s plus uniformes et des sous-ensembles élémentaires de plus petite dimension dans les zones de l'image les moins uniformes. A cet effet, pour chaque sous- ensemble élémentaire en cours de traitement, un critère de complexité, dit de "qualité", est déterminé et une estimation d'opportunité de partition est effectuée. Lorsque cette estimation ne répond pas au critère de qualité, c'est-à- dire lorsque que le sous-ensemble élémentaire considéré est supposé présenter une certaine complexité, il est partitionné en sous-ensembles élémentaires de plus petite dimension qui sont traités à leur tour, jusqu'à ce qu'une dimension minimale ait été atteinte. Dans les modes de réalisation décrits et représentés, cette partition consiste à diviser le sous-ensemble élémentaire en quatre sous-ensembles élémentaires égaux, de forme homothétique à celle du sous-ensemble élémentaire initial et décallés entre eux dans chaque direction. Premier mode de traitement secondaire des données 3 o compressées. On va maintenant décrire un premier mode de mise en oeuvre du procédé de représentation secondaire d'un ensemble de données constitués par des paramètres résultant d'un traitement primaire de données, conformément au premier mode de mise en oeuvre du procédé de traitement primaire décrit ci-dessus. Dans ce premier mode de réalisation, les paramètres issus du traitement primaire font l'objet du traitement secondaire, en vue de restituer, avec le taux d'erreur mentionné supra, les données initiales. Dans d'autres modes de réalisation, dont l'un est décrit infra, le 1o traitement secondaire a pour objet de délivrer des données d'un autre type que celui des données initiales. D'une manière générale, dans le cas d'une restitution des données initiales, il est procédé a une étape de "calcul itératif" au cours de laquelle, à partir desdits paramètres (al,...,a",b) on met en oeuvre une méthode d'approximations successives convergeant vers le point fixe de l'application contractive globale (T) ce point fixe constituant une représentation dudit ensemble de données initiales (f E E). Pour mettre en oeuvre la méthode d'approximations successives dans le cadre d'une restitution, on applique conformément au premier mode de mise en oeuvre de l'invention, sur un ensemble de données arbitrairement choisi (x0 E E), I'application contractive globale (T*) telle que définie plus haut (équations (8) et (10)). On applique ensuite la même contraction (T*) sur la transformée ainsi obtenue et on réitère ce processus jusqu'à ce que la distance entre deux transformées consécutives soit inférieure à un seuil donné. Ceci s'exprime par l'équation suivante: T* (T(T(o... T*(T* (x())....(T (xo (30) En appliquant q fois l'application contractive globale sur le résultat de la transformée précédente, on obtient une bonne approximation de fi ce qui s'exprime par: (T* r(xo) f (31) le critère d'arrêt des itérations étant exprimé comme suit: d((T *)I(xo), ( T*)/+I (xo)) I'invention. Dans la description qui précède, I'application contractive globale définie aux équations (8) et (10) est utilisée dans le processus de restitution de données. En effet, comme le but du traitement secondaire décrit plus haut est, dans ce mode de réalisation, d'obtenir une restitution aussi fidèle que possible des données initiales, c'est la même application contractive globale qui est utilisée lors de l'étape de restitution. En revanche, dans la variante décrite infra, à l'appui des figures 9, 9A, 18 et 19 o l'on souhaite obtenir, à l'issue du traitement secondaire, des données d'un type différent de celui des données initiales, une autre application contractive globale, et/ou un autre ensemble de définition seront utilisés. D'une manière générale, I'application contractive globale utilisée lors du traitement primaire est dite de "premier type"; celle utilisée lors du traitement secondaire est dite de "second type". En cas de restitution des données initiales, I'application contractive globale de second type est identique à celle de premier type. Dans le cas o au cours du traitement secondaire les données d'un type différent de celui des données initiales sont obtenues, l'application contractive globale de second type sera différente de celle de premier type. 2. Description d'un premier mode de réalisation d'un dispositif de traitement Dispositif de traitement primaire d'un ensemble de données S On va maintenant décrire à l'appui de la figure 4 un premier mode de réalisation d'un dispositif de traitement d'un ensemble de données représentatives de grandeurs physiques, et destiné à permettre une représentation primaire dudit ensemble par la mise en oeuvre du procédé décrit supra. Dans ce mode de réalisation, qui concerne un dispositif de traitement d'images en couleur, les données de l'image en 512 x 512 pixels sont admises en entrée, tandis qu'en sortie le dispositif délivre les paramètres ai, i =,.... 4 et b et des informations de partition, relatifs à chacun des sous- ensembles constituant une partition de l'image. On notera ici que la présente invention n'est certainement pas limitée à un format d'image particulier mais qu'elle permet le traitement d'images de toutes tailles. Toutefois, dans un mode de réalisation décrit et représenté, l'image traitée et restituée est une image de 512 x 512 pixels à 256 niveaux sur chacune des trois composantes Y, luminance, Cr, chrominance rouge, et Cb, chrominance bleue. Sur le schéma-bloc de la figure 4, le dispositif de codage porte la référence générale 500. Il comporte, raccordés à un bus de données et d'adresses, désigné sous la référence 501: - un contrôleur (ALU + CTR) 550 incorporant une unité de logique arithmétique. Ce contrôleur est constitué, dans ce mode de réalisation, par un microprocesseur du type de celui commercialisé par la Société INTEL sous la référence i486, - une mémoire vive RAM, portant la référence 510. Cette mémoire o10 comporte un certain nombre de registres qui seront décrits ci-après et une partie adressable rapidement (cache) 514, - une mémoire programme 530, de type ROM, - un premier registre 560 de mémoire vive à accès direct (REG 1), - un second registre 570 de mémoire vive à accès direct (REG 2), - un tampon d'entrée (IBFR) 590, raccordé à un connecteur d'entrée 591 et - un tampon de sortie (OBFR) 580, raccordé à un connecteur de sortie 581. Dans la ROM 530 est enregistré un programme permettant le traitement des données en entrée 591, avec compression de celles-ci conformément au procédé décrit ci-dessus. Le programme de compression, qui porte sur la figure 4, la référence 531, est décrit infra à l'appui des figures 10 à , 19,20et21. L'ensemble des moyens 510, 530, 531, 550, 560 et 570 constitue ici un mode de réalisation de "premiers moyens logiques de transformation" au sens de l'invention. Dans un mode de réalisation particulier, I'ensemble de ces éléments est incorporé dans un seul et même circuit intégré spécifique (ASIC). Dans ce mode de réalisation, la mémoire vive RAM 510 comporte, entre autres les registres suivants destinés à la mise en oeuvre du programme 531: - un registre 511 (IM), destiné à mémoriser les données de l'image, un compteur 513 (CNT) utilisé ainsi qu'il est expliqué infra (figure ). la partie 514 (CACHE) de la mémoire vive. Cette partie est une partie de mémoire adressable très rapidement, dans laquelle sont, ainsi qu'il est connu, stockées les données les plus couramment utilisées à un instant donné. La partie (CACHE) 514 de la mémoire 510 comporte les registres suivants: - un registre 515 CSI destiné à mémoriser, à un instant donné, les coordonnées du sous-ensemble élémentaire Ki en cours de traitement, des registres f(L,)- f(L4) référencés 516 à 519, dans lesquels sont enregistrés les valeurs [f(Li)](x,y) correspondent à chaque point (x,y) (niveau parmi 256 de la composante en cours de traitement) des sous- ensembles de travail L, - L4 associés au sous-ensemble élémentaire Ki considéré, - un registre 520 destiné à mémoriser les coefficients ap, a2, a;, a4 et b, relatifs à un sous-ensemble élémentaire considéré Ki, et - un registre 522 destiné à conserver des facteurs de division A correspondants au nombre de partition complémentaires effectuées pour obtenir le sous-ensemble élémentaire considéré Ki, - les registres 521 et 523 à 526, destinés à conserver des variables intermédiaires d'un traitement générant les références de sous- ensembles élémentaires formant une partition initiale d'une image à compresser, - le registre 560 (REG 1), destiné à cumuler, pour chaque sous- ensemble élémentaire, les coefficients a,, a2, a3, a4 et h, ainsi que l'éventuelle information de partition, relatifs à chaque sous-ensemble de travail. - le registre 570 (REG 2), destiné à mémoriser la liste des adresses de tous les sous-ensembles élémentaires Ki. Sur la figure 5, on a représenté, raccordés respectivement aux moyens d'entrée et de sortie: - une source 592 de données non compressées, et - des moyens 582 utilisateurs de données compressées. La source 592 peut comporter de nombreux moyens à la portée de l'homme du métier. Par exemple, elle pourra consister en un lecteur de disque dur ou de disquette ou de compact disques (connus sous les initiales CD-ROM). Les données de l'image à compresser pourront donc avantageusement être enregistrées sur les supports correspondants: disque dur, disquette ou CD-ROM (compact disque). Cet enregistrement peut être fait de toute manière connue de l'homme de l'art. La source 592 pourra également être constituée par un interface vidéo, susceptible de délivrer les données de l'image à comprimer à l'entrée 591. La source 592 pourra également être constituée par des moyens de réception des données raccordées à un réseau de transmission, la source 592 étant susceptible de transformer les données reçues en un format adéquat pour qu'elles puissent être présentées en entrée du dispositif de traitement primaire. De même, les moyens utilisateurs 582 pourront consister en des moyens de stockage de l'image compressée (disque dur, disquette, CM - ROM etc...) ou un simple interface destiné à délivrer sous un format connu le signal vidéo, à d'autres moyens utilisateurs, tels que des moyens destinés par exemple à transmettre cette image comprimée à distance ou à la stocker. Ces moyens 582 pourront d'ailleurs consister dans un dispositif de transmission de données raccordé à un réseau de télécommunications. Dispositif de traitement secondaire d'un ensemble de paramètres résultant du traitement primaire de données initiales. Dans ce mode de réalisation, qui concerne le traitement d'images, le dispositif de traitement secondaire est apte à restituer, à partir des coefficients a, a2, a3, a4 et b et des informations de partition, relatifs à chacun des sous-ensembles élémentaires, I'image initiale, en couleur. Il est également apte à délivrer les données relatives à une image l0 à niveaux de gris correspondant à une image originale en couleurs. Il est, en outre, adapté à fournir une image possédant un nombre de points élémentaires différent du nombre de points élémentaires de l'image traitée par le dispositif de traitement primaire, ces images étant dites de définition ou de résolution différentes. Ces données en sortie du dispositif de traitement secondaire sont is d'un autre type que le type des données de l'image initiale ayant fait l'objet du traitement primaire. Sur le schéma-bloc de la figure 6, le dispositif de restitution ou de calcul itératif porte la référence générale 600. Il comporte, raccordés à un bus de données et d'adresses, désigné sous la référence 601: - un contrôleur (ALU + CTR) 650 incorporant une unité de logique arithmétique. Ce contrôleur est constitué, dans ce mode de réalisation, par un microprocesseur du type de celui commercialisé par la Société INTEL sous la référence i486. - une mémoire vive RAM, portant la référence 610, qui comporte un certain nombre de registres qui seront décrits ci-après et une partie adressable rapidement (dite "mémoire cache") 614. - une mémoire programme 630, de type ROM. Dans cette mémoire, sont enregistrés notamment: 3 0. un programme 631 de décompression des données de l'image présentée en entrée 681 en vue de la reconstruction d'une approximation de l'image initiale de résolution identique ou différente, un programme 632 de décompression des données de l'image en vue de la reconstruction d'une image à niveaux de gris en lieu et place d'une image couleur, et un facteur de zoom z dans un registre "z" 628. Les programmes 631 et 632 mettent en oeuvre le procédé décrit ci-dessus et leurs organigrammes sont détaillés ci-après à l'appui des figures 16 à 17. - un premier registre 560 à accès direct (REG 1), et - un second registre 570 à accès direct (REG 2) identiques à ceux décrits à l'appui de la figure 4, - un tampon de sortie (OBFR) 690, raccordé à un connecteur de sortie 691 et - un tampon d'entrée (IBFR) 680, raccordé à un connecteur d'entrée 681. L'ensemble des moyens 610, 630, 631, 650, 560, et 570 constitue ici un mode de réalisation de "seconds moyens logiques" au sens de l'invention. L'ensemble des moyens 610, 630, 632, 650, 560 et 570 constitue ici un autre mode de réalisation de "seconds moyens logiques" au sens de l'invention. Dans un mode de réalisation particulier, I'ensemble de ces éléments est incorporé dans un seul et même circuit intégré spécifique (ASIC) convenant pour le traitement secondaire. Dans ce mode de réalisation, la mémoire vive RAM 610 comporte, entre autres les registres suivants destinés à la mise en oeuvre des programmes 631 et 632: - un registre 511 (IM), décrit à l'appui de la figure 4, - un registre 621 (RE), destiné à mémoriser les données d'une image, - un registre (CHTY) 622, destiné à différencier le traitement secondaire général d'une restitution, - un registre (IT) 612 destiné à mémoriser la valeur courante du nombre d'itérations effectuées pour réaliser la représentation secondaire, - un registre (p) 623 destiné à mémoriser le nombre d'itérations pour mettre en oeuvre la méthode des approximations successives afin d'effectuer le traitement secondaire de l'image, - un compteur 513 (CNT) décrit à l'appui de la figure 4. Ce compteur est ici utilisé ainsi qu'il est expliqué infra (figure 16), - la partie 614 (CACHE) de la mémoire vive 610. Cette partie est une partie de mémoire adressable très rapidement, dans laquelle sont, ainsi qu'il est connu, stockées les données les plus couramment utilisées à un instant donné, la partie (CACHE) 614 de la mémoire 610 comporte des registres 515 à 520 et 627 respectivement identiques aux registres 515 à 520 et 522 décrits à l'appui de la figure 4, - le registre 560 (REG 1), destiné à stocker, pour chaque sous- ensemble élémentaire, les coefficients a. a2, a3, a4 et h et l'éventuelle information de partition, associés à ce sous-ensemble, et - le registre 570 (REG 2), destiné à mémoriser la liste des adresses de tous les sous-ensembles élémentaires K,. Sur la figure 7, on a représenté, raccordés respectivement aux moyens d'entrée et de sortie: - une source 682 de données compressées, et - des moyens 692 utilisateurs de données décompressées. La source 682 pourra consister en des moyens de stockage de l'image comprimée (disque dur, disquette, CD-ROM etc...) ou un simple interface destiné à délivrer celles-ci. Les moyens 682 pourront aussi consister dans un système de réception de données raccordé à un réseau de télécommunications. Les moyens 692 pourront consister par exemple en des moyens de stockage des données décompressées ou en un simple interface délivrant sur sa sortie les données décomprimées. Ils pourront également consister en un dispositif de transmission des données décompressées sur un réseau de télécommunications. Description d'un dispositif incorporant un dispositif de traitement primaire et un dispositif de traitement secondaire Sur le schéma bloc de la figure 8, un dispositif incorporant un dispositif de traitement primaire 500 et un dispositif de traitement secondaire 600 est représenté et porte la référence générale 700. Ce dispositif est capable de traiter parallèlement des compressions d'images et des décompressions. Il est dès lors composé des dispositifs décrits à l'appui des figures 4 et 6. Il s'agit en l'occurrence du dispositif de traitement primaire 500 et du dispositif de traitement secondaire 600. Ce type d'appareil est appelé "full-duplex". Dans un autre mode de réalisation (half-duplex), le dispositif ne peut exécuter qu'un traitement primaire ou un traitement secondaire mais jamais en même temps. Il peut en outre réaliser un traitement "complexe" comportant un traitement primaire et un traitement secondaire réalisés l'un à la suite de l'autre. Le schéma bloc d'un tel dispositif est représenté dans les figures 9 et 9A. Ce dispositif qui porte la référence 800 comporte conjointement les ressources des dispositifs 500 et 600 décrites supra, en évitant le double emploi. C'est ainsi que: - les registres des mémoires vives 510 et 610 décrits supra sont incorporés dans une même unité de mémoire vive 810, une même unité de mémoire morte 830 incorpore les programmes 531, 631 et 632 mentionnés supra et un programme de configuration 831, - un même microprocesseur 850 incorpore les moyens 550 et 650. Le programme de configuration 831 qui sera décrit infra à l'appui de la figure 18 permet au dispositif 800, de réaliser: - une compression des données d'une image couleur, - une décompression des données soit en vue de restituer l'image couleur initiale, soit en vue de la restituer avec une résolution différente, soit en vue de délivrer une image à niveaux de gris correspondantes, et - un traitement complexe permettant de transformer une image couleur en une image à niveaux de gris. Description d'un organigramme particulier de compression d'une image couleur 512 x 512 pixels i5 Cet organigramme est illustré en figure 10. Les figures 11 à 15 sont des organigrammes explicitant certaines opérations de l'organigramme de compression de la fig u re 10. Les organigrammes des figures 11 et 12 détaillent également par ailleurs des opérations de l'organigramme de décompression des figures 16 et 17 qui sera décrit infra. Les références entre les paramètres renvoient aux étapes de la figure 16. On rappelle tout d'abord que les données relatives à l'image matérialisent une application telle que définie dans l'équation (9) exposée supra. En l'espèce, à chaque couple de valeurs (xy) (coordonnées colonne, ligne) pris dans le domaine de définition (en l'espèce à l'image plane 512 x 512 pixels) correspond une valeur fjtx,y) dans le domaine d'application R3', c'est à dire une valeur parmi 256 niveaux relatifs à chacune des trois composantes Y, Cr et Cb. On organise trois suites de nombres correspondant chacun au niveau d'une composante (soit 512 x 512 = 262.144 valeurs /(x,y) par composante), en balayant les données de l'image comme suit: - les lignes sont balayées de gauche a droite; - les lignes sont incrémentées de bas en haut. Ces trois suites sont ici enregistrées dans la source de données non compressées 592. Elles peuvent être lues sur le port d'entrée 591 à la commande du microprocesseur 550. En 200, il est procédé à une initialisation des registres (IM) 511, (CNT) 513 et des deux registres (REG 1) 560 et (REG 2) 570 à la valeur o10 uniforme nulle. Les fichiers relatifs aux composantes Y, Cr et Cb sont accessibles séparément et sont fournis dans un ordre prédéterminé, à savoir: - d'abord les données de l'image considérée et relatives a la composante de luminance (Y), ensuite les données de l'image relatives à la composante de chrominance rouge (Cr) et finalement les données de l'image relatives à la composante de chrominance bleue (Cb). En 201, il est procédé à un premier test aux fins de déterminer si le compteur (CNT) 513 est ou non égal à 3. Dans la négative le programme se branche sur l'étape 202, au cours de laquelle: - le compteur 513 (CNT) est incrémenté d'une unité; - le microprocesseur 550 va lire sur le port d'entrée 591 la composante suivante à traiter. On rappelle que chaque composante, Y, Cr ou Cb, se présente comme une suite de 262.144 valeurs /(x,y) correspondant au niveau de la composante considérée pour chacun des pixels de l'image. Cette suite est chargée dans le registre (IM) 511. On rappelle que, dans ce mode particulier de réalisation, chaque sous-ensemble de travail a autant de données que les sous-ensembles élémentaires. Le programme se branche ensuite à l'étape 204 (explicitée en figure 11), au cours de laquelle il va être procédé à la partition initiale de l'image de 512 x 512 pixels en 256 sous-ensemble élémentaires Ki, i = (,. ,256) de 32 x 32 pixels. Cette partition se fait en déterminant les..DTD: coordonnées (x,y) des coins inférieurs gauche de chacun des ni sous- ensembles Ki considérés. C'est ainsi que le sous-ensemble K, illustré en figure 3, a pour coordonnées (xiyi), qui correspondent à celles de son coin inférieur gauche. Chacune des adresses des coins inférieurs gauches des blocs de 32 x 32 pixels correspond à une adresse particulière dans la liste enregistrée dans le registre (IM) 511 de la mémoire vive 510. i0 La liste de ces adresses est chargée dans le registre (REG 2). En référence à la figure 11, après une initialisation d'un compteur constitué par le registre i 521 relatif à l'indice des sous-ensembles K, en l'étape 211, le programme se branche au test 212 au cours duquel il est vérifié si la liste est complète ou non. Dans la négative, le programme se branche à l'étape 213 au cours de laquelle le facteur de division A, qui décompte chaque partition complémentaire qui est effectuée sur un sous-ensemble élémentaire, est initialisée à la valeur nulle et o l'on détermine une variable B comme le quotient d'une division entière de (i-1) par 16. Le reste de la division (variable C) est mémorisé dans le registre 524. Les valeurs x, = C * 32 et y, = B: 32 sont calculées au cours de l'opération 214 et temporairement mémorisées dans les registres 525 et 526. Ces coordonnées qui sont celles du coin inférieur gauche du sous-ensemble élémentaire Ki considéré, et le facteur de division A associé, sont mémorisées dans les variables (x,,y) dans le registre 570 (REG 2) au cours de l'étape 215. Pendant l'étape 216 le compteur i est incrémenté et le programme exécute ensuite à nouveau le test 212. Lorsque le test 212 est positif, le programme effectue le test 205 (figure 10) au cours duquel il est vérifié si la liste enregistrée dans le registre (REG 2) est vide ou non. Dans la négative, cela signifie qu'il y a encore des sous-ensembles élémentaires K, à traiter et le programme se branche à l'étape 206. A l'étape 206, il est supprimé le premier sous-ensemble élémentaire de la liste enregistrée dans le registre (REG 2). Ce sous-ensemble est enregistré dans un registre CSI 515 de la mémoire cache de la RAM 510. Le programme se branche alors à l'étape 207 détaillée en figure s 19 ou, selon une variante, en figure 20. Au cours de cette étape, un test 255 détermine si le sous-ensemble élémentaire considéré possède les dimensions minimales, c'est-à-dire dans le mode de réalisation décrit et représenté, 2 fois 2 pixels, ou non. Lorsque le résultat du test 255 est positif, l'opération 250 consiste à déterminer les quatre sous-images de travail Li associées au sousensemble de travail considéré. Cette opération 250 est détaillée en figure 12. Au cours de l'opération 250, à partir de l'adresse enregistrée dans le registre CSI et du facteur de division A, il est procédé à la recherche des quatre sous-ensembles de travail voisins L,, L2, L3, L4 (voir figure 3 pour les sous-ensembles de travail associés à la partition initiale). A cet effet, la règle rappelée supra est appliquée - chaque sous-ensemble de travail possède les mêmes dimensions que le sous-ensemble élémentaire auquel il est lié, soit dans chaque direction: 32 divisé par 2 porté à la puissance A; et - chaque sous-ensemble de travail est décalé d'une demi dimension dans chaque direction, par rapport au sous-ensemble élémentaire associé. L'opération 250 a donc pour but de déterminer les valeurs f/(x,y) de la composante considérée associée à chacun des pixels des sous-ensembles de travail L,, L2, L3 et L4 liés géométriquement et de manière prédéterminée à lI'ensemble élémentaire considéré Ki. L'ensemble des valeurs /(x,y) liées aux pixels du sous-ensemble LI (respectivement L,, L3 et L4) sera noté ci-après J(L,) (respectivementj(L2) J(L3) etj(L4)). En référence à la figure 12, pour chaque sous-ensemble élémentaire K, considéré, le contrôleur 550 lit dans le registre 570 (REG 2) les coordonnées (xiy) du pixel inférieur gauche de K,, à l'étape 217. Ensuite, à l'étape 218, l'unité arithmétique 550 lit les valeurs de la composante considérée pour chaque pixel des quatre sous-ensemble de travail Lt à L4, définis selon le tableau suivant: L,: [x,-16/2A, x,-1 + 16/2A)] x [y, + 16/2A, y -1 +(3x16)/2"] L2: [xi+16/2A, x i-1 + (3x16)/2A] x [Yi + 16/2A, Yi-1 +(3x16)/2A] L3: [x-16/2A, x,-1 + 16/2A] x [y,- 16/2, yi-] +16/2"] L4: [xi+16/2A, x i- + (3x16)/2A] x [y- 16/2A, yi- +16/2'] On observe ici que chacun des sous-ensembles de travail possède ici le même nombre de pixels que le sous-ensemble élémentaire considéré, soit (32/2)2. Ensuite, à l'étape 219, l'unité arithmétique 550 va chercher, pour la composante, Y, Cr ou Cb, en cours de traitement, la valeur correspondant à chacun des (32/24)2 pixels de L, (respectivement L,, L3 et L4) et stocke l'ensemble des (32/12)2 valeurs dans le registre f(L,) 516 (respectivement 1(L,) 517, f(L3) 518 etj(L4) 519) du cache 514). A la suite de l'opération 250 (figure 19), le programme effectue l'opération 251 au cours de laquelle, d'une manière générale, les paramètres a, (en l'espèce ap, a2, a3, a4) et b sont calculés en mettant en oeuvre la méthode des moindres carrés avec les contraintes décrites supra (voir équations 12 et 13). Une description simplifiée de l'étape 251 est donnée dans l'organigramme représenté dans la figure 13. On rappelle que pour déterminer les coefficients al, a2, a3, a4 et b il y a lieu de résoudre un problème de minimisation d'une fonctionnelle sous contraintes. Supra il a été vu, à l'appui des équations 14 à 29 que cela revient à résoudre itérativement deux types de problèmes de minimisation sans contraintes. Les explications données supra pour un cas général sont mises en pratique pour la détermination des quatre 3 0 paramètres aI, a2, a3, a4 et du paramètre b. A l'étape 220 on procède à la résolution du premier problème de minimisation sans contraintes. On minimise la fonctionnelle: +K (a,,a4, b)= Y f(x, Y)_ a,. +JL bk-" K(a(x.y)EK, ô LI()I= Ix pour des valeurs (a,,..., a4,b) parcourant tout R5. Pour ce faire on met en oeuvre la méthode des moindres carrés décrite supra par les équations (16) à (19). La solution unique à ce problème (a,,..., a4,b) dans R5 est la solution du 1o système linéaire (19). Comme illustré en figure 14, l'étape 220 est subdivisée en quatre étapes 230 - 233. Au cours de l'étape 230, on calcule des coefficients Gu i = 1,...,5, j = 1,...,5 en mettant en oeuvre les formules suivantes: Gji = E [f(Lj)](x, Y) [f(Li)](x, y) ij l.4 (x,y)EK G5s = Gi5 = [f(Li)] (x,y) i= (x.y) K G55 = 322 Les coefficients [Gj] constituent la matrice à inverser afin de résoudre le système linéaire d'équations (19). Les formules ci-dessus permettant de déterminer les coefficients Gu, sont comme on le voit des simples additions et multiplications mettant en oeuvre des valeurs de pixels correspondant aux niveaux de la composante considérée, et cela pour chacun des pixels du sous ensemble élémentaire considéré. De même à l'étape 231 on calcule les coefficients, %'.. S5 constituant le second membre du système linéaire. Les formules mises en oeuvre sont: =, z S1 E |[f (L, Y)])(xy) f(xy)i=...4 (x,y)EK S5= E f(x,y) (xy)EK L'étape 232 fait appel à une routine classique d'inversion d'un système d'équations linéaires, en l'espèce à cinq inconnues. Le système linéaire à inverser est le suivant: -G I G21 G31 G41 G51 al SI G12 G22 G32 G42 G52 a2 S2 G13 G23 G33 G43 G53 a3 S3 GI4 G24 G34 G44 G a4 5 4 G15 G25 G35 G45 G55 b _SJ La solution a,.... a4, b sera mémorisée à l'étape 233 dans le registre 520 de la partie 514 (CACHE) de la mémoire vive 510. Le programme passe ensuite à l'étape 221 (figure 13) o, conformément aux équations (20) à (23), on détermine une partition disjointe de l'ensemble des entiers {1,2,3,4} en NV, UVet LVdéfinis comme NV = { E {1, .,4}|- 1 Le programme se branche ensuite au test 222 au cours duquel il est vérifié si les ensembles UVet LV sont vides. Si ce test est positif les valeurs actuellement présentes dans le registre 520 constituent (étape 224) la solution du problème avec contraintes. Dans la négative, certaines contraintes sont actives pour la solution du problème avec contraintes. Dans ce cas, (étape 223 voir figure 15) il y a lieu de résoudre un problème restreint. Il s'agit en l'occurrence de minimiser sous contraintes la fonctionnelle modifiée: UJ(CIC,...,cm)= Z F(xy) - C Ji ()+X,Y) + Cm (xy)eK y Li 1l o F(x,y) = f(x,y)- V. (Lj)x, y)+ f.[f(L"3)](X.y), jeUV k ELI' et m La solution de ce problème est (C....,cni) et est obtenue en considérant les équations (16) à (18) et en résolvant le système linéaire (19). La solution du problème initial avec contraintes est: ai =cj: j = 1,...,m1 ou i e NV, ak = k e UV, a, = -3. le LV, b =Cm. L'étape 223 va être décrite à l'appui de l'organigramme de la figure 15 en supposant m = 3. Il convient de considérer le problème de minimisation sans contraintes, décrit dans l'étape 223 qui est subdivisée en quatre étapes 240 à 243. Au cours de l'étape 240, on calcule des coefficients Hij i = 1,... 4 j = 1.... 4 en mettant en oeuvre les formules suivantes: H-, = [=(L,)](XIy). ['(Li)] (xy). ... .3 (x,y)EK H4i = H4:= [fLj)](xy) (x,y)EK.. H44 = 322 Les coefficients [HI,] constituent la matrice à inverser afin de résoudre le système linéaire d'équations (19). Les formules ci-dessus permettant de déterminer les coefficients H, sont comme on le voit des simples additions et multiplications mettant en oeuvre des valeurs de pixels correspondant aux niveaux de la composante o considérée, et cela pour chacun des pixels du sous ensemble élémentaire considéré. De même à l'étape 241 on calcule les coefficients T,... T, constituant le second membre du système linéaire. Les formules mises en oeuvre sont: Ti = [f (L)](x Y)- F(x, y).=.3 (x,y)EK T4 = Z F(x,y) (X,y)EK L'étape 242 fait appel à une routine classique d'inversion d'un système d'équations linéaires, en l'espèce à 4 inconnues. Le système linéaire à inverser est le suivant: HI Hi, H31 H41 C Ti H12 H22 H 32 H42 c2 T7 H,3 H23 H33 H43 c3 T3 H,4 H24 H34 H44 T4r A l'issue de l'étape 223, il y a lieu de revoir si la solution ainsi obtenue satisfait encore les contraintes. Pour ce faire, l'étape 221 et le test 222 sont exécutés à nouveau. Ce processus est réitéré jusqu'à ce que le test 222 soit positif. * Au cours de l'opération 251, les coefficients ai, a,, a3, a4, b sont enregistrés dans le registre 560 (REG 1). Puis le test 205 est réitéré. On observe ici que dans la mesure o le sous-ensemble élémentaire possède des l0 dimensions minimales, aucune information de partition n'est nécessaire, puisqu'aucune partition complémentaire ne peut être effectuée. Lorsque le résultat du test 255 est négatif, l'opération 252 consiste à mesurer la variance des niveaux de la composante considérée, dans le sous- ensemble élémentaire considéré. Ensuite, l'opération 253 consiste à déterminer un critère dit "de qualité" applicable au sous-ensemble élémentaire considéré. Préférentiellement, la définition du critère de qualité consiste en la définition d'un seuil en deçà duquel le critère de qualité est considéré comme respecté. Préférentiellement aussi, la valeur du seuil prend en compte les dimensions du sous-ensemble élémentaire considéré. Dans le mode de réalisation décrit et o représenté ce seuil est la fonction correspondant au tableau suivant: dimensions du sous-ensemble en pixels seuil 2x2 200 4x4 160 8x8 100 16 x 16 60 32 x 32 20 Selon des variantes non représentées: - la valeur du seuil est constante quelle que soit les dimensions des sous-ensembles élémentaires ou - la valeur du seuil est une fonction inverse d'une dimension du sous-ensemble élémentaire considéré. Le test 254 détermine si le critère de qualité est respecté ou non. Dans le mode de réalisation décrit et représenté, il est respecté lorsque la mesure de la variance des niveaux de la composante considérée est inférieure à la valeur du seuil qui correspond, dans le tableau cidessus, aux dimensions l0 du sous-ensemble élémentaire considéré. Lorsque le résultat de l'un des tests 254 ou 255 est positif, l'opération 250 consiste à déterminer les quatre sous-images de travail comme exposé ci-dessus. Ensuite, l'opération 258 est identique à l'opération 251, si ce n'est que sont mémorisés, dans le registre 560 (REG 1), successivement I'information de partition valant "0" puisqu'aucune partition complémentaire n'a été effectuée et les coefficients ai, a2, a3, a4, b. Ensuite, le test 205 est réitéré. Lorsque le résultat du test 254 est négatif, l'opération 256 consiste à effectuer une partition complémentaire du sous-ensemble élémentaire considéré. Dans le mode de réalisation décrit et représenté, cette partition 2 0 complémentaire consiste à diviser chaque sous-ensemble élémentaire en quatre sous ensembles élémentaires carrés égaux. A cet effet, au cours de l'opération 256, il est ajouté en tête de la liste liste_sousimage les coordonnées des quatre points qui définissent les coins inférieurs gauches des quatre sous-ensembles élémentaires résultant de la partition complémentaire. A ces coordonnées sont associées un facteur de division A qui est égal au facteur de division du sous-ensemble élémentaire que l'on partitionne, incrémenté de la valeur 1. En outre une information de partition valant "1", puisqu'une partition complémentaire a été effectuée, est mémorisée dans le registre 560 (REG 1). A la suite de l'opération 256, le test 205 est réitéré. Les figures 21A, 21B et 21C illustrent l'effet des partitions complémentaires qui sont effectuées sur les sous-ensembles élémentaires issus de la partition initiale ou de partitions complémentaires précédentes. Dans un but explicatif, et par exception aux règles générales de la description, on suppose que, dans le mode de réalisation exposé en figures 21A à 21C, les dimensions minimales des sous-ensembles élémentaires sont quatre fois plus petites que leurs dimensions maximales. La figure 21A représente une partie d'une image, correspondant à douze sous-ensembles élémentaires référencés B0 à B1i résultant de la lo partition initiale (opération 204, figure 11), dans laquelle: - le sousensemble élémentaire B1 a subi une partition complémentaire en quatre sous-ensembles élémentaires B10o à B1,3, - le sous-ensemble élémentaire B1,3 a subi une partition complémentaire en quatre sous-ensembles élémentaires B1,3,0 à B1,33, - le sous-ensemble élémentaire B6 a subi une partition complémentaire en quatre sous- ensembles élémentaires B6,0 à B6,3, - le sous-ensemble élémentaire B6,1 a subi une partition complémentaire en quatre sous-ensembles élémentaires B6,1,o à B6,13, et le sous-ensemble élémentaire Bll a subi une partition complémentaire en quatre sous-ensembles élémentaires B.,0 à B1,3. La figure 21B illustre, pour les sous-ensembles élémentaires ainsi partitionnés, les informations de partition qui sont associées aux sous- ensembles élémentaires résultants de la partition initiale (sousensembles pour lesquels la valeur du facteur de division A est nulle) ou aux sous-ensembles élémentaires résultants de partitions complémentaires. Ainsi, les sous-ensembles élémentaire B0, B2, B3, B4, B5, B7, B8, B9 et B10 sont associés à une information de partition nulle, les sousensembles élémentaires B1, B6 et Bl, sont associés à une information de partition valant 1, leurs sous ensembles élémentaires B1,0, B.,1, B1,2, B6,0, B6,2, B6,3, B110, B,11, B11,2 et B1,3 sont associés à une information de partition nulle, les sous ensembles élémentaires B.,3 et B6,1, sont associés à une information de partition égale a "1" et les sous-ensembles élémentaires résultant de leur partition complémentaire ne sont associés a aucune information de partition puisque les dimensions de ces sous-ensembles élémentaires sont les dimensions minimales. Les dimensions de tous ces sous-ensembles élémentaires peuvent ainsi être définies sans ambiguïté par la suite d'informations de partition La figure 21C propose une représentation des paramètres a,, a2, a3, a4 et b, qui, conjointement, sont, référencés ici D et des informations de partition qui leur sont associées, en utilisant la suite exposée cidessus. On y observe que la séquence des données compressées correspondant à la partie d'image représentée en figure 21A comporte: - un en-tête, représentant en particulier les dimensions de l'image (valant ici 512 x 512 pixels) les dimensions initiales des sous-ensembles élémentaires (correspondant ici à 32 x 32 pixels) et la définition des dimensions minimales des sous-ensembles élémentaires (correspondant ici à 2 x 2 pixels) et une information de traitement complémentaire et l'éventuel ratio de dimensions entre la résolution de la composante de luminance et celle des composantes de chrominance; - la suite des informations de partition présentée ci-dessus: associée aux éventuelles données qui leur correspondent, au cours de la mémorisation effectuée pendant les étapes 251 (sans information de partition) et 258 (avec information de partition), sachant que les informations de partition valant I ne sont associées à aucune donnée. L'ensemble des opérations 205-258 ou 205-256 est réitérée tant que la mémoire (REG 2) n'est pas vide Lorsque le test 205 (figure 10) est positif, le programme se branche à nouveau au test 201 aux fins de vérifier si le compteur (CNT) 513 est ou non égal à 3. Dans l'affirmative cela signifie que les trois composantes ont été traitées. Le programme se branche alors sur l'étape de détermination de traitement complémentaire 209, au cours de laquelle il est procédé: - d'une part, à la détermination de l'utilité ou non d'effectuer un traitement complémentaire de lissage sur l'image résultant d'un traitement secondaire des paramètres et informations de partition et - d'autre part, à la mémorisation de l'en-tête du fichier des paramètres et informations de partition, en-tête comportant notamment les dimensions de l'image, les dimensions maximales et minimales des sous- ensembles élémentaires, I'éventuel ratio de dimensions entre la résolution de la composante de luminance et celle des composantes de chrominance, et l'information de traitement complémentaire représentant l'utilité dudit traitement complémentaire. Pour déterminer cette utilité, plusieurs critères peuvent être utilisés. Le premier consiste à comparer la dimension minimale des sous- ensembles élémentaires à une dimension qui sert de seuil. Lorsque la dimension minimale est supérieure à la dimensions seuil, cela signifie que des sous-ensembles élémentaires relativement grands (par exemple de 8 x 8 pixels) vont être utilisés même si le critère de qualité, défini à l'étape 253, n'est pas respecté et qu'il sera donc profitable de réduire la sensation de bord affectant de tels sous-ensembles élémentaires. Le second consiste à comparer les seuils utilisés dans le tableau définissant le critère de qualité avec des valeurs dites "critique", de seuils prédéterminées. Lorsque les seuils effectivement utilisés sont supérieurs à ces valeurs critiques, cela signifie que le critère de qualité est aisément respecté sur des sous-ensembles élémentaires de grandes dimensions, sous-ensemble dont les bords peuvent être visibles. Il est alors, de nouveau, intéressant de prévoir un traitement 3 o complémentaire. Ensuite est réalisée l'étape de fin 210. Les paramètres de l'application contractive globale se trouvent associées aux informations de partition, dans le registre (REG 1) 560. Ces paramètres constituent une représentation primaire des données initiales qui étaient enregistrées dans la source de données 592. Les performances suivantes ont pu être observées sur une image couleur 512 x 512 pixels à 256 niveaux par composante: Nombre de données au départ: 3x 512 x 512 x 8bits = 6 291456 Nombre de données à l'arrivée: 3x 256 x (4 x 64 + 10) = 204288 par conséquent: - taux de compression: 30 - nombre de bits par pixels: 0,78 Temps de traitement: 3 x 256 x 80 msec.= 61 sec. (80msec. = temps unitaire de calcul de a,...., a) Bien entendu, I'invention s'applique aussi bien au traitement primaire d'images en niveaux de gris (possédant une seule composante) qu'à celui d'images en couleurs (possédant plusieurs composantes). Description de deux modes de réalisation d'un organigramme de décompression On va décrire à l'appui des figures 16 et 17 I'organigramme de programme 631 permettant de reconstruire les données de l'image couleur d'origine. Au cours de l'opération 400, figure 16, il est procédé à une initialisation des registres (IM), (RE) et (IT) de la mémoire vive 610, des deux mémoires (REG 1) 560 et (REG 2) 570, et du compteur (CNT) 513, à la valeur uniforme nulle. Le test 401 détermine si le contenu du compteur 513 (CNT) est ou non égal à 3. Lorsque le résultat du test 401 est positif, les trois composantes ont été restituées et le programme de décompression est terminé. Le programme effectue alors le test 416 décrit infra. Dans la négative, le programme effectue l'opération 403 au cours de laquelle les données compressées de la mémoire d'entrée sont chargées dans la mémoire (REG 1) 560. Il est rappelé que les données compressees sont formées d'une suite de cinq ou six nombres, à savoir les coefficients ap, a2, a3, a4 et b et, éventuellement, d'une information de partition, relatifs à chaque sous-ensemble élémentaire Ki. Dans le cas d'une image couleur de 512 x 512 pixels, chaque composante a été compressée en suites au cours du traitement primaire tel que mis en oeuvre par le programme 531 décrit supra. Cette séquence de suites de nombres est chargée dans (REG 1) 560. Au cours de l'opération 420, le facteur de zoom z est lu dans la mémoire morte ROM 630. On note ici que la sélection du facteur de zoom z peut aussi se faire par tout moyen à la portée de l'homme du métier. En particulier, il pourrait être fait au moyen d'un commutateur (non représenté) 2 0 agencé sur le boîtier du dispositif de traitement secondaire objet de la présente invention, commutateur à plusieurs positions correspondant chacune à un facteur de zoom prédéterminé. Ce commutateur permettrait d'envoyer les signaux appropriés à la mémoire vive RAM 610, de manière à charger les valeurs indiquées dans un registre non représenté z. Cette initialisation pourra également être faite dans des conditions similaires au moyen d'un clavier (non représenté), moyennant des dispositions techniques non représentées mais bien connues de l'homme du métier. Le programme effectue ensuite l'opération 404, au cours de laquelle il est procédé à la création de la liste des références de chaque sous- ensemble élémentaire. Il s'agit en l'occurrence d'une liste des coordonnées (xi,y,) du pixel inférieur gauche de chaque sous-ensemble élémentaire K,, dont le nombre de pixels est de ((z x 32)/2X)2 pixels, le facteur de division A étant recalculé en fonction des informations de partition comme exposé avec le traitement primaire. Cette liste est déterminée en prenant en compte I'information de partition qui est associée à chaque sous-ensemble élémentaire ainsi que le facteur de zoom z qui est associé à toute l'image résultant du traitement secondaire. On note ici que le facteur de zoom z est un facteur multiplicatif qui s'applique aux dimensions de l'image reconstruite au cours du traitement 1o secondaire à partir des paramètres provenant d'une image traitée au cours d'un traitement primaire de compression. Ce facteur multiplicatif s'applique aussi aux dimensions de chaque sous-ensemble élémentaire, qu'il provienne de la partition initiale ou de partitions complémentaires, et à chaque dimension de sous-ensemble de travail associé à un de ces sous-ensembles élémentaires. Lorsqu'un ratio entre les résolutions des composantes de luminance a été transmis dans l'en-tête, le facteur de zoom z est appliqué à la luminance et le produit du facteur de zoom z par ledit ratio est appliqué aux composantes de chrominance, de telle manière que tous les ensembles de données sur lesquels on reconstruit une composante possèdent la même résolution. Conformément à l'invention, aucun facteur de zoom z ne provoque une duplication de pixel formant l'apparition, dans l'image reconstruite, de pavés uniformes. Le programme se branche ensuite sur l'étape 406 au cours de laquelle: conformément au mode de mise en oeuvre de l'invention décrit ici, une image arbitraire initiale est mémorisée dans le registre (IM) 511 de la mémoire vive 610. Dans le cas présent de décompression d'une image couleur o chaque composante possède 256 niveaux ayant respectivement les valeurs 0 à 255, on initialise le registre (IM) à la valeur numérique 128 (correspondant sensiblement à la valeur moyenne entre 0 et 255) pour chaque pixel pour la première composante traitée. Pour les composantes suivantes, on utilise comme initialisation le résultat précédemment obtenu. - le registre (IT) 612 est initialisé avec une valeur p lue dans le registre 623 de la RAM 610. Le nombre d'itérations est enregistré dans le registre (p), compte-tenu du RMSE souhaité (voir supra tableau 1 de la description générale du procédé). Les inventeurs ont en effet constaté au cours d'essais avec le prototype que, pour obtenir un RMSE donné à partir d'un taux de compression, un nombre limité d'itérations successives au cours du calcul d'approximation était nécessaire. Typiquement: 2 images très complexes p peut avoir une valeur supérieure. - le compteur (CNT) 513 est incrémenté d'une unité. Le programme exécute ensuite un nombre p d'itérations référencé sur la boucle 450 afin de permettre la restitution d'une composante de l'image. Le programme se branche tout d'abord au test 407 au cours duquel on vérifie si le contenu du registre (IT) 612 est positif et non nul. Si le test 407 est positif, le programme se branche au test 408 au cours duquel il est vérifié si la liste enregistrée dans (REG 2) 570 est vide ou non. Si le test 408 est négatif, le programme se branche sur l'étape 2 0 410 et exécute celle-ci jusqu'à ce que tousles sous-ensembles élémentaires ki aient été traités. Le détail de l'étape 410 sera maintenant décrit dans la figure 17. L'étape 410 comprend elle-même les étapes 411 à 415: A l'étape 411, on retire les coordonnées (x,,y,) correspondants au sous-ensemble élémentaire Ki considéré et on les copie dans le registre CSI 515. Parallèlement, lescoefficients (a,...,. a4,b),relatifs au sous-ensemble élémentaire considéré sont copiés du registre (REG 1) vers le registre 520 de la partie 614 (CACHE) de la mémoire vive 610. L'étape 412 a pour but de déterminer les sous-ensembles de travail liés géométriquement et de manière prédéterminée a l'ensemble élémentaire Ki considéré, prenant en compte le fait que chaque sous- ensemble de travail possède les mêmes dimensions que le sous-ensemble élémentaire considéré (tenant ainsi compte du facteur de zoom z) et est décalé, par rapport à ce sous-ensemble élémentaire, d'une demi dimension. Cette étape comporte les même opérations qu'à l'étape 250 décrite supra à l'appui de la figure 12. Les valeurs fiL/(), /(L2), /(L3) et 1(L4) des pixels relatifs respectivement à L,, L2, L3 et L4 sont respectivement copiées dans les registres 516, 517, 518 et 519 de la partie 614 (CACHE) de la mémoire vive 610. Le programme se branche ensuite aux étapes 413 et 414 au cours desquelles on calcule pour chaque pixel (x,y) de Ki la valeur E ai [f(Li)](x, y)+ b i-= à l'aide des coefficients stockés dans le registre 520 et des données f(L,)..., f(L4) stockées dans les registres 516 à 519, (étape 413) et on mémorise (étape 414) cette valeur dans le registre (RE) 621 de la RAM 610. Si le test 408 est positif, le contenu du registre (RE) 621 est copié dans le registre (IM) 511, et le contenu du registre (IT) 612 est diminué d'une unité. Il est à nouveau procédé au test 407. Tant que ce test est positif cela signifie que le nombre d'itérations prévu par le nombre enregistré dans le registre (p), n'est pas terminé et les opérations 408 à 410 sont réitérées. Lorsque le test 407 est négatif ou nul, le nombre d'itérations prévu est terminé, pour la composante en cours de traitement, et le programme sort de la boucle 450. Cette composante est donc décompressée. On procède à l'étape 405 qui consiste à copier le contenu du registre (IM) 511 dans le tampon de sortie (OBFR) 690. Le programme se branche ensuite au test 401. Tant que le test 401 est négatif, cela signifie qu'il reste au moins une composante à décompresser. On observe que pour la seconde et la troisième composante il n'est pas procédé à une nouvelle initialisation du registre (IM) 511. Le registre (IM) comporte en effet la liste des valeurs correspondant à chaque pixel pour la composante précédemment décompressée. Cette liste de valeur est l'image arbitraire initiale pour la composante suivante. La figure 16A illustre à cet égard le déroulement de l'étape 406. Cependant, selon une variante non représentée, il est procédé de manière parallèle aux traitements des trois composantes, Y, Cr et Cb, de l'image, en prenant comme ensemble de données initial, un ensemble de données arbitraires. On observe que, selon des caractéristiques générales de l'invention, le traitement primaire comporte une étape de "construction", opération 251 ou 258, d'une application contractive globale dite "de premier type", pour ledit premier ensemble de données, dont le point fixe constitue une approximation de tout ou partie de cet ensemble, en mettant en oeuvre, pour ce faire, au moins une application appartenant au groupe d'applications constitué des applications multidimensionnelles et des applications non linéaires. Les paramètres ai et b de l'application contractive globale permettent la mise en oeuvre, sur un deuxième ensemble de données comportant un nombre de données qui dépend du facteur de zoom z, et qui est généralement différent du nombre de données du premier ensemble, d'une méthode d'approximations successives convergeant vers le point fixe d'une application contractive globale. En figure 16A, au cours du test 430, il est procédé au test du compteur 513 (CNT). Si son contenu est égal à zéro, cela signifie qu'il s'agit de la première composante, celle de luminance, et le registre (IM) 511 est initialisé à 128 pour chaque pixel (étape 431). Le registre (IT) est initialisé a la valeur p (étape 432). Si le contenu du compteur 513 (CNT) est supérieur à zéro, cela signifie qu'une composante au moins a été décompressée et, à l'étape 433, on charge le registre (IM) 511 avec le contenu du registre (RE) 621, c'est- à-dire avec les valeurs correspondant à la composante précédemment décompressée. Ceci permet de décroître le nombre d'itérations p de deux unités (cf. étape 434 o l'on initialise (IT) à la valeur p-2) car l'image initiale est très proche de l'image finale. Le temps de calcul est diminué grâce à cette caractéristique. Le compteur 513 (CNT) est ensuite incrémenté d'une unité (étape l0 435). Lorsque les trois composantes ont été décompressées, le test 401 est positif. Le test 416 consiste alors à déterminer si: - soit l'information de traitement complémentaire présente dans l'en-tête des paramètres et informations de partition est représentatif d'une utilité d'effectuer un traitement complémentaire, - soit le facteur de zoom z possède une valeur qui rend utile un traitement complémentaire (z étant, par exemple, supérieur ou égal à 4) ou non. Lorsque le résultat du test 416 est négatif, l'étape de fin 402 est effectuée directement. Lorsque le résultat du test 416 est positif, l'opération 417 de traitement complémentaire est effectuée. Cette opération 417 consiste en une opération connue par l'homme du métier sous le nom de "post-processing", ce qui fait référence à un traitement effectué après la décompression d'image illustrée en figure 16. Ce traitement comporte, de manière connue, un filtrage spatial, destiné à diminuer l'éventuelle sensation de frontière des sous- ensembles élémentaires. Préférentiellement, les paramètres de ce filtrage sont fonctions des dimensions du sous-ensemble traité. Par exemple, les dimensions du filtre utilisé sont proportionnelles à celles du sous-ensemble considéré, d'une part, et le nombre de pixels sur lesquels s'applique ce filtre est proportionnel à chaque dimension du sous-ensemble élémentaire considéré, d'autre part. Selon chacun de ces exemples, le facteur de zoom intervient dans le paramétrage du traitement complémentaire. A la fin de l'opération 417, l'étape de fin 402 est effectuée. Au cours de l'étape de fin 402, les données relatives à chacune des trois composantes sont, tour à tour, copiées dans le tampon de sortie (OBFR) 690 d'o elles sont transférées vers les moyens utilisateurs de données décompressées 682 (figure 7). On observe ici que, pour obtenir une image reconstruite en niveaux de gris, sans couleur, il suffit de n'effectuer un traitement primaire que sur les données provenant du traitement primaire de la luminance de l'image couleur, cette limitation étant par ailleurs désignée sous le nom de programme 632. Description de l'organigramme du programme de : 15 configuration On va maintenant décrire à l'appui de la figure 20 l'organigramme du programme de configuration mis en oeuvre dans le dispositif 800 des figures 9 et 9A. Le programme de configuration porte la référence 831 aux figures 9A et 20. Le dispositif 800 illustré en figure 9A possède quatre modes de fonctionnement qui seront détaillés ci-après. Le mode de fonctionnement de l'appareil est tout d'abord sélectionné (étape 832). La sélection du mode se fait en rentrant dans le registre MOD 822, une valeur égale à 0, 1 ou 2 et dans le registre (CHTY) une valeur égale à 0 ou 1. Cet enregistrement peut se faire par tout moyen à la portée de l'homme du métier. En particulier, il pourrait être fait au moyen d'un commutateur agencé sur le boîtier de l'appareil 800 à quatre positions correspondantes chacune au mode de fonctionnement souhaité. Ce commutateur permettrait d'envoyer les signaux appropriés à la mémoire RAM 810, de manière à charger les valeurs indiquées dans les registres MOD ou (CHTY). Cette initialisation pourra également être faite dans des conditions similaires au moyen d'un clavier, non représenté sur la figure 9A. A l'étape 833, il est procédé à un premier test sur le registre MOD aux fins de déterminer si son contenu est ou non égal à zéro. Si le test 833 est positif, cela signifie que le mode de fonctionnement choisi est le premier de ceux mentionnés ci-dessus et le programme 531 est mis en oeuvre à partir de données alimentées depuis la source de données 810, les données compressées étant délivrées en sortie 881 aux moyens utilisateurs de données 0o compressées 860 (figure 9). Si le test 833 est négatif, il est procédé, à l'étape 834, à un test aux fins de déterminer si le contenu du registre MOD est ou non égal à 1. Si le test 834 est positif, il est alors procédé à un test 840 sur le contenu du registre (CHTY) aux fins de déterminer s'il est ou non égal à 1. Si le test 840 est négatif, cela signifie que le dispositif des figures 9 et 9A est dans le second mode de fonctionnement, c'est-à-dire dans celui o il est utilisé pour restituer les données d'origine d'une image couleur à partir de paramètres (données compressées) et d'informations de partition, fournies en entrée 881, à partir d'une source 861 (figure 9). Les données décompressées délivrées en sortie 892 du dispositif 800 sont destinées à des moyens utilisateurs de données décompressées 820. Si le test 840 est positif, cela signifie que le dispositif 800 est dans le troisième mode de fonctionnement, c'est-à-dire celui o il est utilisé pour réaliser un traitement secondaire de données compressées et délivrer des données d'un autre type que celui des données originales, en l'espèce des données relatives à une image à niveaux de gris correspondant à une image couleur originelle. Dans ce cas, le programme 632 est mis en oeuvre à partir de paramètres et informations de partition délivrés par une source de données compressées 861. Les données décompressées sont délivrées en sortie 892 et 3 o alimentées vers des moyens utilisateurs 820. Si le test 834 est négatif, cela signifie que le contenu du registre MOD est égal a 2. Le dispositif 800, qui est alors dans son quatrième mode de fonctionnement, est utilisé pour transformer une image couleur en image à niveaux de gris. Il s'agit ici d'un exemple de traitement complexe au sens de la présente invention, comportant une phase de traitement primaire et une phase de traitement secondaire. Ainsi, dans le mode de réalisation décrit ici, conformément à cet aspect de l'invention, les programmes de traitement primaire et de traitement secondaire sont mis en oeuvre séquentiellement. Le programme 531 permet, à partir de données délivrées par une source de données non compressées 810 sur l'entrée 891 du dispositif 800, de compresser ces données, pour obtenir des paramètres a,. hb, A relatifs à chacun des sous-ensembles élémentaires Ki de l'image d'origine. Ces paramètres et informations de partition sont, conformément à ce qui a été décrit supra à l'appui de la figure 10, tous présents dans le registre (REG 1) 560 à la fin du programme 531. Le programme principal de la figure 20 commande alors un branchement direct à l'issue de l'étape 404 du programme 631 décrit supra. Les données transformées sont alors présentes sur le registre tampon (OBFR) 893 du dispositif 800 et transférées vers les moyens 820 utilisateurs de données décompressées. Ces données sont celles relatives à une image à niveaux de gris correspondant à l'image couleur dont les données ont été transférées dans le dispositif 800 à partir de la source 810. On notera que le sous-programme comportant l'ensemble des étapes 404, 452, 406, 450, 405 et 402 du programme 632 constitue ici un mode particulier de mise en oeuvre d'un procédé de traitement secondaire au sens de l'invention. C'est ce sous-programme particulier qui est ici mis en oeuvre à la suite du programme 631 pour constituer un exemple d'un traitement "complexe" au sens de la présente invention. 3. Description d'une variante préférée mettant en oeuvre des sous-ensembles de travail présentant une taille différente de celle du sous-ensemble élémentaire correspondant Lorsque les sous-ensembles de travail possèdent la même dimension que les sous- ensembles élémentaires auxquels ils sont associés, et que l'application contractive est linéaire sur chaque sous-ensemble de travail, l'image reconstruite ne peut avoir dans le sous-ensemble élémentaire considéré une complexité plus grande que celles qu'elle a dans chaque sous- l0 ensemble de travail. En outre, la convergence de la méthode d'approximations successives peut être relativement lente. Pour ces deux raisons, et au moins dans le cas de la mise en oeuvre d'applications contractives multilinéaires, il est préférable d'utiliser des sous-ensembles de travail possédant des dimensions supérieures à celles des sous-ensembles élémentaires. Aussi, on va maintenant décrire à l'appui des figures 3A, 10, 11, 12A, 13 à 16 et 17A, une variante préférée de réalisation du procédé et dispositif décrit supra, dans laquelle au moins un sous-ensemble de travail 2 0 présente une taille différente de celle du sous-ensemble élémentaire correspondant. Dans le mode de réalisation particulier de cette variante, quatre sous-ensembles de travail sont associés à chaque sous-ensemble élémentaire, ces quatre sous-ensembles de travail présentant des dimensions supérieures à celles du sous-ensemble élémentaire correspondant. L'image à traiter est ici une image de 512 x 512 pixels partitionnée initialement en sous-ensembles carrés de 32 x 32 pixels, conformément à ce qui a été décrit supra à l'appui de la figure 1. On va maintenant décrire la division initiale de l'image en sous- ensembles de travail, sachant que chaque partition complémentaire d'un sous- ensemble élémentaire donne lieu à la définition de sous-ensembles de travail de dimensions proportionnelles à celle de ce sous-ensemble élémentaire. Les sous-ensembles de travail présentent ici une taille double de celle des sous- ensembles élémentaires, de telle sorte que l'image de 512 x 512 pixels est divisée en sous-ensembles carrés de 64 x 64 pixels. Toutefois, de manière à ce qu'à un sous-ensemble élémentaire donné correspondent quatre sous- ensembles de travail (voir figure 3A), le décalage utilisé pour la division de l'image est de 32 pixels, la division de l'image résultant alors en un recouvrement par ardoises (en anglais "tiling"). En figure 3A le sous-ensemble élémentaire K, ayant pour coordonnées (xiy) est illustré en traits gras, tandis que les sous-ensembles de travail correspondants sont illustrés en traits pointillés, les autres sous- ensembles étant ébauchés en traits fins. On observe que, dans ce mode de réalisation, à chacun des sous-ensembles Ki sont liés quatre sousensembles de travail Lú (numérotés L,,..., L4) décalés de 64 pixels dans les deux directions et disposés symétriquement par rapport au sousensemble Ki considéré. Il s'agit en réalité des quatre sous-ensembles issus de la division qui recouvrent le sous-ensemble élémentaire considéré (issu de la partition initiale). Les expressions ci-dessous montrent les coordonnées d'un sous-ensemble K, et des sous-ensembles de travail L,,.... L4, pour toutes valeurs du facteur de division A: Tableau 2 Ki: [xi, xi - 1 +(32/2)] x [y, yi-1 +(3212/)] LI [xi - (48/2A), x -1 + (16/2A] x [y, + (16/'2), yi -1 + (80/2")] L2: [xi + (16/2A), xi -1 + (8012A] x [y, + (16/2Y), yi -1 + (80/2A)] L3 [xi- (48/), xi -1 + (16/2A] x [yi - (48/2), yi -1 + (16/2)] L4 [xi + (16/2A), x -1 + (80/2A] x [yi (48/2A), yj-1 + (16/2A)] On rappelle ici que les points des sousensembles de travail définis ci-dessus qui n'appartiennent pas à l'image, c'est-à-dire ceux dont l'une au moins des coordonnées est négative, se voient fictivement associer une valeur moyenne, par exemple 127, pour chaque composante. Les dispositifs de compression et de décompression utilisés pour mettre en oeuvre cette variante de réalisation sont semblables à ceux décrits à l'appui des figures 4 et 6. Toutefois, les programmes 531 et 631 sont modifiés ainsi qu'il va être expliqué ci-après. S'agissant tout d'abord du programme de compression 531 décrit supra à l'appui de la figure 10, sa structure est conservée. Seule change l'opération 250. Cette opération modifiée 250A va être décrite ci-après à l'appui l0 de la figure 12A. A l'étape 250A, il est procédé à la recherche des quatre sous- ensembles de travail voisins LI, L2, L3, L4 (voir figure 3A) à partir de l'adresse enregistrée dans le registre CSI. A cet effet, la règle rappelée supra est appliquée. is L'étape 250A a pour but de déterminer les valeurs f/(x,y) de la composante considérée, associée à chacun des pixels des sous-ensembles de travail L1, L2, L3, L4 liés géométriquement et de manière prédéterminée à l'ensemble élémentaire considéré Ki. L'ensemble des valeurs f/(x,y) liées aux pixels du sous- ensemble LI (respectivement L2, L3 et L4) sera noté ci-après f(1LI) (respectivement fL2),./(L3) et.f/(L4)) conformément aux notations utilisées supra. L'élément essentiel qui différencie l'étape 250A (figure 12A) de l'étape 250 (figure 12), réside dans le fait que l'ensemble de valeurs représentées par./(Li), i = 1,.... 4 concerne des sous-ensembles de travail Li, i = 1,.... 4 qui ont une taille différente de celle du sous-ensemble élémentaire K,. Afin de pouvoir exécuter les opérations de comparaison avec le sousensemble K,, il y a lieu de ramener les ensembles f(L) à la même taille. Plusieurs méthodes sont à la disposition de l'homme du métier. Dans ce mode de réalisation o la taille des sous-ensembles de travail est de 64 x 64 pixels et celle des sous-ensembles élémentaires de 32 x 32 pixels, on retient une valeur moyenne par carré de 4 pixels. L'ensemble de ces valeurs moyennes sera noté pour un ensemble de travail Li considéré: f/(L). L'organigramme détaillé de cette étape est représenté dans la figure 12A. En référence à la figure 12A à l'étape 317, on copie dans l'unité arithmétique 550 les coordonnées (xi,y) du pixel inférieur gauche du sousensemble élémentaire Ki considéré ainsi que le facteur de division A, et on les stocke dans le registre (REG 2) 570. Ensuite, I'unité arithmétique 550 détermine les coordonnées de chaque pixel de chaque sous-ensemble de travail selon le tableau 2 comme indiqué à l'étape 318. A l'étape 319, I'unité arithmétique 550 va chercher, pour la composante en cours de traitement, la valeur correspondant à chacun des (64/2')2 pixels de L1 (respectivement L2, L3 et L4) et stocke l'ensemble des (64/2A)2 valeurs dans le registresj(L) 516 (respectivement /(L2) 517, f/(L,) 518 et /(L4) 519) du cache 514. 1.5 A l'étape 320, I'unité arithmétique 550 traite chaque carré de 2 x 2 valeurs deJ(L,) et y substitue une valeur correspondante de./(f L) au registre 516. Dans ce mode de réalisation la valeur de f1(L[) correspondante sera obtenue en calculant la moyenne des valeurs dudit carré de 2 x 2 valeurs de f(L1). Les autres ensembles f(L2).... j(L4) sont ensuite traités de la même façon. A l'étape 321, les valeurs de (L)......(,L4) sont stockées dans la partie 514 (CACHE) de la mémoire vive 510, dans les registres./(L,) 516 àf(L4) 519. Le programme se branche ensuite à l'étape 208 (figure 13) au cours de laquelle, comme indiqué supra, les paramètres a, (en l'espèce a1, a2, a3, a4) et b sont calculés en mettant en oeuvre la méthode des moindres carrés avec les contraintes décrites supra (voir équations 12 et 13). Toutefois, dans ce mode de réalisation, la valeur de a sera inférieure à -12. On observe que cette dernière valeur est plus élevée que celle utilisée dans le premier mode de réalisation. En, effet, dans la mesure o, au cours de l'étape 320, on a placé dans les registres 516 à 519 des sousensembles de travail de substitution dont la valeur de chaque pixel est calculée comme la moyenne des valeurs de quatre pixels de l'image initiale, on a déjà effectué une application contractive. Les contraintes appliquée à la contractivité de l'application globale qui succède à cette substitution peuvent, en conséquence, être réduites. Le fonctionnement du programme de compression est, pour le reste, identique au fonctionnement du programme de compression 531 décrit supra à l'appui de la figure 10. La seule différence étant que, dans les registres 516 à 519, sont stockées les valeurs dLú),..., f(L4). Dans cette variante de réalisation, le programme de décompression fonctionne de façon similaire au programme 631 décrit supra à l'appui des figures 16 et 17. Toutefois, l'opération 410 doit être modifiée comme exposé ci-après et sera référencée comme 41 OA, illustrée en figure 17A. Le détail de l'étape 410A est décrit à l'appui de la figure 17A. L'étape 410A comprend elle-même les étapes 459 à 465: Les étapes 459 et 460 sont identiques aux étapes 411 et 412 décrites supra (figure 17), mais elles prennent en compte, pour la définition des sous- ensemble élémentaires et des sous-ensembles de travail le facteur multiplicatif dit "de zoom" z qui multiplie chacune de leur dimension dans un ensemble de données qui possède lui-même des dimensions qui sont multipliées par ce facteur de zoom. Le programme se branche ensuite aux étapes 461, 462 et 463 au cours desquelles on détermine les valeurs f(L),...,f(L4). Ensuite à l'étape 462 on détermine les valeursf/LI)..(, iLs)j. Pour ce faire, l'unité arithmétique 550 traite chaque carré de 4 valeurs de /(Ll) et y substitue une valeur correspondante de.Lú) au registre 516. Dans ce mode de réalisation la valeur de I(LI/) correspondante est obtenue par calcul de la moyenne des valeurs dudit carré ou tout autre filtrage passe-bas bien connu de l'homme du métier. Les autres ensembles f(L2)..., f(L4) sont ensuite traités de la même façon. A l'étape de 463 on calcule pour chaque pixel (x,y) de K, la valeur a [f (L)](, y)+ h i=l à l'aide des coefficients de 520 et les données.(Lú)..... (L4) de 516 à 519, (étape 462) et on mémorise (étape 463) cette valeur au pixel correspondant dans le registre (RE) 621 de la RAM 610. Finalement, l'étape 464 est identique à l'étape 414 (figure 17). Le fonctionnement du programme de décompression modifié comme exposé cidessus à l'appui de la figure 17A est pour le reste identique au fonctionnement du programme de décompression 631 décrit supra à l'appui de la figure 16. Dans cette variante de réalisation, il peut également être réalisé un traitement complexe permettant de transformer des données, comme décrit supra à l'appui de la figure 20, en mettant en oeuvre les programmes de compression 531 et de décompression 631 modifiés tel qu'exposé ci-dessus à l'appui des figures 10, 11, 12A, 13,14 et 17A. 4. Description d'une variante de réalisation mettant en oeuvre une application contractive globale composée d'applications non linéaires Dans cette variante de réalisation des procédés et dispositifs conformes à l'invention, il est mis en oeuvre une application contractive globale composée d'applications non linéaires. Dans cette variante de réalisation, qui concerne également la compression et la décompression d'images, I'image à traiter est également une image 512 x 512 pixels. Cette image est initialement partitionnée en sous- ensembles élémentaires et divisée en sous-ensembles de travail de la même façon que ce qui est décrit supra à l'appui des figures 1 à 3. On va maintenant définir l'application globale utilisée dans un mode particulier de cette variante de réalisation. Conformément à cet aspect de l'invention, chaque "morceau" de l'application contractive globale, est une application élémentaire Lipschitzienne non linéaire du type: i T(f) Ki ' (x,y) = Ki --> [T(f)] (x,y) =a+ g[] i (x,y) +, 1=1 (10 ter) o: T)I Ki la restriction de T(f) au sous-ensemble élémentaire Ki g1 [f,., gj[l] sont des transformations non linéaires de f(g, E F (E. E)) tels que des puissances fractionnaires ou entières de f (par exemple g[f] = JXi) n > I le nombre de sous-ensembles de travail (n = 4 dans ce mode de réalisation) a/,... an les coefficients multiplicatifs h le facteur de translation Li, L, les sous- ensembles de travail, à position prédéterminée par rapport à Ki et lié à ce sous-ensemble. On voit que l'équation (10 ter) ci-dessus est de nature non linéaire essentiellement à cause des transformations non linéaires,P.....g.., les coefficients a1, j =...n et b étant respectivement les coefficients multiplicatifs et le facteur de translation. Les coefficients(a1,...,a,,, bh sont ici les "paramètres" au sens de l'invention de l'application élémentaire relative au sous-ensemble K,. L'ensemble des coefficients (aI,..., an, b)K,. relatifs aux m sous-ensembles élémentaires (KI,.... Km,,) constitue ici aussi l"'ensemble des paramètres" de l'application contractive globale au sens de l'invention. On rappelle que, dans ce mode de réalisation, n = 4, ce qui signifie que seulement quatre sous-ensembles de travail Lj. = 1.... n tels que ceux définis à l'appui de la figure 3, sont liés à chaque sous-ensemble élémentaire Ki. Ainsi, à chaque couple de valeurs (x,y) d'un sous-ensemble élémentaire Ki déterminé, on associe une valeur calculée à partir de l'équation non linéaire (10 ter). Ceci est fait pour l'ensemble des m sous- ensembles élémentaires (K,.., Km) résultants aussi bien de la partition initiale que de partition complémentaires, comme décrit supra. Pour résoudre le problème d'optimisation simplifié (voir supra, notamment équation (7)) qui permet de déterminer les coefficients de I'application non linéaire a.,,j =,.... n et b pour chacun des sous- ensembles élémentaires de travail Ki, il faut, d'une manière générale, pour chaque sous- ensemble élémentaire Ki, minimiser la distance entre la restriction au sous- ensemble élémentaire K, et sa transformée par l'application contractive globale. En d'autres termes, il faut minimiser les distances: d (.fTl,, r(f ^-,), i =....ni (11) A cet égard, on applique la théorie exposée à l'appui des équations (12) à (29) dans la partie relative à la théorie générale. Par ailleurs, on observera que, d'une manière générale dans cette variante de réalisation, le procédé de représentation secondaire des paramètres issus du traitement primaire, de données initiales, conforme à la variante décrite ci-dessus, met en oeuvre les principes généraux décrits supra à l'appui des équations (30) à (32). Les dispositifs de compression et de décompression utilisés pour mettre en oeuvre cette variante de réalisation sont semblables à ceux décrits supra à l'appui des figures 4 et 6. Toutefois, les programmes 531 et 631 sont modifiés ainsi qu'il va être expliqué ci-après. S'agissant tout d'abord du programme de compression 531 décrit supra à l'appui de la figure 10, sa structure est conservée. Seule change l'étape 251 relative aux calculs des coefficients. Cette variante de l'étape 251 (référencée comme 251B) va maintenant être décrite notamment à l'appui de la figure 12B. L'étape 251B comporte les étapes 322 et 323 respectivement identiques aux étapes 217 et 218 décrites supra à l'appui de la figure 12. A l'étape 324, on lit dans le registre (IM) 511, les valeursf(L,)...., /f(L4). A l'étape 325, on détermine les valeurs g1 [f (L,)J. ,g4[f(L4)] en calculant par application de la formule g,[f(L, ): (Li) Enfin, à l'étape 326, on mémorise dans les registres 516 à 519 de la partie 514 (CACHE) de la mémoire vive 510, les valeurs e. [f (Li Il...* 14 [f(L.4z Le fonctionnement du programme de compression avec l'opération 250 modifiée ainsi qu'il vient d'être décrit est, pour le reste, identique..DTD: au fonctionnement du programme 531 décrit supra à l'appui de la figure 10. * De même, dans cette variante de réalisation, le programme de décompression fonctionne de façon similaire au programme 631 décrit supra à l'appui des figures 16 et 17. Toutefois l'opération 410 doit être modifiée comme exposé ci-après. Cette nouvelle opération est référencée 410B et est illustrée dans la figure 17B. Le détail de l'étape 410B sera maintenant décrit à l'appui de la lo figure 17B. L'étape 410B comprend elle-même les étapes 471 à 475: A l'étape 471, on retire les coordonnées (x,,,y,) du sous-ensemble élémentaire Ki considéré et on les copie dans le registre CSI 515. Les coefficients (a,..., a4, b, relatifs aux sous-ensembles Ki considérés sont copiés du registre (REG 1) vers le registre 520 de la partie 614 (CACHE) de la mémoire vive 610. L'étape 472 a pour but de déterminer les sous-ensembles de travail liés géométriquement et de manière prédéterminée à l'ensemble élémentaire Ki considéré. Cette étape comporte les même opérations qu'à l'étape 250 décrite supra à l'appui de la figure 12. Les valeurs g1 [fi] i, g2 Lf] /2, g3 [f] l,3 et g4 [f], 4 des pixels relatifs respectivement à LI, L2, L3 et L4 sont respectivement copiées dans les registres 516, 517, 518 et 519 du (CACHE) 614. Le programme se branche ensuite aux étapes 473 et 474 au cours desquelles on calcule pour chaque pixel (x,y) de K, la valeur 4 a [ ' Y) + i=a,.[gi[/],,J](x, y) +/h à l'aide des coefficients conservés dans le registre 520 et les données 1'(L,)...,.f(L4) conservés dans le registre de 516 à 519, (étape 474) et on mémorise cette valeur au pixel correspondant dans le registre (RE) 621 de la mémoire vive RAM 610. Dans cette variante de réalisation, il peut également être réalisé un traitement complexe permettant de transformer des données comme décrit supra à l'appui de la figure 20, en mettant en oeuvre des programmes de compression 531 et de décompression 631 modifiés tel qu'exposé ci-dessus à l'appui des figures 12B et 17B. 5. Description d'une variante de réalisation mettant en oeuvre un critère de qualité mesuré à l'aide de l'application contractive globale Selon une variante, le critère de qualité utilisé pour déterminer la nécessité d'effectuer une partition complémentaire n'est pas, comme présenté en regard de la figure 19, la comparaison, d'une part, d'une mesure statistique effectuée sur les niveaux de la composante considérée dans un sous-ensemble élémentaire avec, d'autre part, un seuil, mais la comparaison d'une mesure de différence entre, d'une part, I'image obtenue après lui avoir appliquée une application contractive globale, avec cette image initiale. Cette variante est illustrée en figure 20, dans laquelle les tests et opérations communs avec la figure 19 portent les mêmes références. Dans cette variante, à la suite de l'opération 206, I'opération 250 est réalisée, en prenant en compte la valeur actualisée du facteur de division A même si cette valeur n'est pas la valeur qui sera, finalement, associée au sous-ensemble élémentaire considéré. L'opération 250 détermine, on le rappelle, les sous- ensembles de travail associés au sous-ensemble élémentaire considéré. Ensuite, l'opération 260 est identique à l'opération 251 (figure 19) à l'exception de la mémorisation dans le registre de mémoire 560 (REG 1). Elle fournit donc des paramètres de l'application contractive globale, dans le sous- ensemble élémentaire considéré. Le test 255 succède à l'opération 260 et détermine, comme exposé supra, si le sous-ensemble élémentaire possède des dimensions minimales ou non. Lorsque le résultat du test 255 est positif, l'opération 263 consiste à mémoriser dans le registre 560 (REG 1) les cinq paramètres déterminés au cours de l'opération 260, sans leur associer d'information de partition puisqu'aucune nouvelle partition complémentaire n'a eu lieu. Ensuite, le test 205 est réitéré. Lorsque le résultat du test 255 est négatif, l'opération consiste à déterminer un critère dit "de qualité" applicable au sous-ensemble élémentaire considéré. Préférentiellement, la définition du critère de qualité consiste en la définition d'un seuil en deçà duquel le critère de qualité est considéré comme respecté. Préférentiellement aussi, la valeur du seuil prend en compte les is dimensions du sous-ensemble élémentaire considéré. Dans le mode de réalisation décrit et représenté ce seuil est la fonction correspondant au tableau suivant: dimensions du sous-ensemble en pixels seuil 2 x 2 800 4x4 640 8x8 400 16 x 16 240 32 x 32 80 On observe ici qu'il est préférentiel d'utiliser un tableau correspondant à une fonction décroissante. Le test 254 détermine ensuite si le critère de qualité est respecté ou non. Dans le mode de réalisation décrit est représenté, il est respecté lorsque la mesure de l'écart quadratique moyen entre la composante considérée et son image par l'application contractive globale, mesure restreinte au sous-ensemble élémentaire considéré, est inférieur à la valeur du seuil qui correspond, dans le tableau ci-dessus aux dimensions du sous-ensemble élémentaire considéré. Lorsque le résultat du test 254 est négatif, l'opération 256 consiste à effectuer une partition complémentaire du sous-ensemble élémentaire considéré. Dans le mode de réalisation décrit et représenté, cette partition complémentaire consiste à diviser chaque sous-ensemble élémentaire en quatre sous ensembles élémentaires carrés égaux. A cet effet, au cours de l'opération 256, il est ajouté en tête de la liste liste_sousimage les coordonnées des quatre points qui définissent les coins inférieurs gauches des quatre sous-ensembles élémentaires résultant de la partition complémentaire. A ces coordonnées est associée un facteur de division A qui est égal au facteur de division du sous-ensemble élémentaire que l'on partitionne, incrémenté de la valeur 1. Enfin, une information de partition valant "1" est is mémorisée dans le registre 560 (REG 1) puisqu'une partition complémentaire a été effectuée. Lorsque le résultat du test 254 est positif, les cinq paramètres déterminés au cours de l'opération 260 sont mémorisés dans le registre 560 (REG 1) avec une information de partition valant "0" puisqu'aucune partition complémentaire n'a eu lieu. Dans la variante exposée en figure 20, il est intéressant de mesurer la moyenne des valeurs de l'écart quadratique moyen de tous les sousensembles élémentaires traités pour lesquels cet écart est mesuré au cours de l'opération 254, pour déterminer l'utilité d'un traitement complémentaire au cours du traitement secondaire: lorsque cette moyenne est supérieure à une valeur critique globale prédéterminée, il est utile de prévoir un traitement complémentaire. On observe ici que les variantes présentées en regard des figures 19 et 20 se combinent particulièrement bien, en faisant se succéder l'opération 262 de la figure 20 à la fin des opérations 258 de la figure 19. Selon des variantes non représentés, les partitions complémentaires sont effectuées en effectuant un découpage non homothétique des sous-ensembles élémentaires, comme par exemple un découpage en triangles, en deux rectangles égaux, successivement en formant deux rectangles côte-à-côte horizontalement ou verticalement,... Selon une autre variante, le facteur de zoom n'est pas uniforme. Il peut ainsi, par exemple, avoir une valeur horizontale, qui sert de coefficient multiplicatif aux dimensions horizontales des sous-ensembles élémentaires et une valeur verticale qui sert de coefficient multiplicatif aux dimensions verticales des sous-ensembles élémentaires. Le facteur de zoom peut aussi avoir une valeur qui varie en fonction de la distance au centre de l'image... Selon une autre variante, préférentielle, le facteur de zoom est différent pour les différentes composantes de l'image traitée. Par exemple, le facteur de zoom est le double pour chaque composante de chrominance de ce qu'elle est pour la composante de luminance. Dans ce cas, il est effectué un filtrage préliminaire au traitement primaire (c'est-à-dire ici à la compression), filtrage qui divise la résolution des composantes de chrominance d'un facteur de deux dans chaque direction, en donnant à chaque composante de chaque pixel de l'image de plus faible résolution la valeur moyenne de la composante correspondante de quatre pixels de l'image initiale. Ce facteur, est transmis dans l'en-tête sous la forme du ratio susmentionné. Conjointement, le nombre de données initiales utilisées pour la reconstruction de la composante de luminance est, dans chaque direction, le double de ce qu'il est pour chacune des composantes de chrominance. En effet, il est courant que la complexité des composantes de chrominance soit moindre que celle de la composante de luminance, d'une part et que cette complexité soit moins perçue par l'oeil humain. En combinaison avec cette dernière variante, le dispositif de traitement primaire est adapté à effectuer, avant les opérations illustrées en figure 10, un filtrage des composantes de chrominance afin de diviser par deux leur résolution dans chaque direction de l'image. Bien entendu la présente invention ne se limite nullement aux modes de réalisation choisis et représentés mais englobe, bien au contraire, toutes variantes à la portée de l'homme du métier. REVENDICATIONS 1. Procédé de traitement primaire d'un premier ensemble de données (IM,f représentatives de grandeurs physiques, caractérisé en ce qu'il comporte une étape de "construction" (205 à 209) d'une application contractive globale dite "de premier type", pour ledit premier ensemble de données, dont le point fixe constitue une approximation de tout ou partie de cet ensemble, en mettant en oeuvre, pour ce faire, au moins une application appartenant au groupe d'applications constitué des applications multidimensionnelles et des applications non linéaires et en déterminant des paramètres (ai, b) de l'application contractive globale de façon à permettre la mise en oeuvre, sur un deuxième ensemble de données comportant un nombre de données différent du nombre de données du premier ensemble, d'une méthode d'approximations successives (407 à 410) convergeant vers le point fixe d'une application contractive globale dite "de second type"; l'ensemble des paramètres (ai, b) ainsi déterminés constituant une représentation primaire dudit ensemble de données. 2. Procédé de traitement primaire selon la revendication 1, caractérise en ce que: - il comporte une étape de "partition initiale" (204) au cours de laquelle, on partitionne ledit premier ensemble de données en sous-ensembles dits "élémentaires initiaux" (Ki), - au cours de l'étape de construction (205 à 209), pour chacun des sous- ensembles élémentaires (Ki), on détermine les paramètres (ai, b) de l'application élémentaire qui constitue une restriction de l'application contractive globale audit sous-ensemble élémentaire. - on permet à la méthode d'approximations successives effectuée sur le deuxième ensemble de données, de prendre en compte des sous-ensembles élémentaires correspondant chacun à un sous-ensemble élémentaire du premier ensemble de données et comportant un nombre de données différent de celui du sous-ensemble élémentaire du premier ensemble de données correspondant. 3. Procédé de traitement primaire selon la revendication 2, caractérisé en ce que les nombres de données des sous- ensembles élémentaires du deuxième ensemble de données sont proportionnels aux nombres de données des sous-ensembles élémentaires du premier ensemble de données qui leur correspondent respectivement. 4. Procédé de traitement primaire selon l'une quelconque des revendication 2 ou 3, caractérisé en ce que ladite application contractive globale est constituée, partiellement au moins, d'applications Lipschitziennes élémentaires liées aux sous-ensembles élémentaires (Ki). 5. Procédé de traitement primaire selon l'une quelconque des revendications 2 à 4, caractérisé en ce que: - au cours de l'étape de construction (205 à 209), pour chacun des sous- ensembles élémentaires (Ki) résultant de la partition initiale, on considère n sous-ensembles de travail (Li), n étant non nul, et on construit une application élémentaire des n sous-ensembles de travail dans le sous-ensemble élémentaire considéré en déterminant les paramètres (ai, b) de l'application élémentaire, - on permet à la méthode d'approximations successives effectuée sur le deuxième ensemble de données, de prendre en compte des sous-ensembles de travail correspondant chacun à un sous-ensemble de travail du premier ensemble de données et comportant un nombre de données différent de celui du sous-ensemble de travail du premier ensemble de données correspondant. 3 o 6. Procédé de traitement primaire selon la revendication 5, caractérisé en ce que les nombres de données des sous-ensembles de travail du deuxième ensemble de données sont proportionnels aux nombres de données des sous-ensembles de travail du premier ensemble de données qui leur correspondent respectivement. 7. Procédé de traitement primaire selon l'une quelconque des revendications 5 ou 6, caractérisé en ce qu'au moins un sous-ensemble de travail (Li) parmi lesdits n sous-ensembles de travail (Li) présente une taille différente de celle du sous-ensemble élémentaire correspondant. 8. Procédé de traitement primaire selon l'une quelconque des revendications 5 à 7, caractérisé en ce que lesdits n sous-ensembles de travail (Li) sont liés géométriquement de façon prédéterminée au sous-ensemble élémentaire considéré. 9. Procédé de traitement primaire selon l'une quelconque des revendications 1 à 8, caractérisé en ce que, ledit premier ensemble de données constituant une image couleur à plusieurs composantes (Y, Cr, Cb), il comporte une opération de filtrage, au cours de laquelle on filtre au moins l'une desdites composantes et on modifie le nombre de données de l'ensemble de données liés à au moins une composante. 10. Procédé de traitement primaire selon la revendication 9, caractérisé en ce que, au cours de l'opération de filtrage, on réduit le nombre de données des ensembles de données liés aux composantes de chrominance, par rapport au nombre de données de l'ensemble de données lié à la composante de luminance. 11. Procédé de traitement secondaire d'un ensemble de paramètres (ai, b) résultant d'un traitement primaire d'un premier ensemble de données, procédé de traitement secondaire caractérisé en ce qu'il comporte une étape de "calcul itératif" (450) au cours de laquelle à partir dudit ensemble de paramètres (ai, b), on met en oeuvre, sur un deuxième ensemble de données possédant un nombre de données différent du premier ensemble de données, une méthode d'approximations successives (407 à 410) convergeant vers le point fixe d'une application contractive globale dite "de second type", en mettant en oeuvre, pour ce faire, au moins une application appartenant au groupe d'applications constitué des applications multidimensionnelles et des applications non linéaires, ce point fixe constituant une représentation secondaire (RE) dudit premier ensemble de données. 12. Procédé de traitement secondaire selon la revendication 11, caractérisé en ce que: - il comporte une étape de "partition initiale" (204) au cours de laquelle, on partitionne ledit deuxième ensemble de données en sous-ensembles dits "élémentaires" correspondants chacun à un sous-ensemble élémentaire du premier ensemble de données et comportant un nombre de données différent de celui du sous-ensemble élémentaire du premier ensemble de données correspondant. 13. Procédé de traitement secondaire selon la revendication 12, 2 0 caractérisé en ce que les nombres de données des sous-ensembles élémentaires du deuxième ensemble de données sont proportionnels aux nombres de données des sous-ensembles élémentaires du premier ensemble de données qui leur correspondent respectivement. 14. Procédé de traitement secondaire selon l'une quelconque des revendications 12 ou 13, caractérisé en ce que, le traitement primaire prenant en compte, pour chacun des sous-ensembles élémentaires (Ki) résultant de la dernière étape de partition, n sous-ensembles de travail (Li), n étant non nul, et construisant une application élémentaire des n sous- ensembles de travail dans le sous-ensemble élémentaire considéré en déterminant les paramètres (ai, b) de l'application élémentaire, - il prend en compte des sous-ensembles de travail correspondant chacun à un sous-ensemble de travail du deuxième ensemble de données et comportant un nombre de données différent de celui du sous-ensemble de travail du premier ensemble de données correspondant. 15. Procédé de traitement secondaire selon la revendication 14, caractérisé en ce que les nombres de données des sous-ensembles de travail du deuxième ensemble de données sont proportionnels aux nombres de données des sous-ensembles de travail du premier ensemble de données qui lo leur correspondent respectivement. 16. Procédé de traitement d'un ensemble de données (IM, f) représentatives de grandeurs physiques, comprenant une phase de traitement primaire mettant en oeuvre un procédé de traitement primaire selon l'une quelconque des revendications 1 à 10 et une phase de traitement secondaire mettant en oeuvre un procédé de traitement secondaire selon l'une quelconque des revendications 11 à 15, la phase de traitement secondaire traitant l'ensemble de paramètres (ai, b) résultant du traitement primaire. 17. Procédé selon l'une quelconque des revendications 11 à 16, caractérisé en ce qu'au cours de ladite étape de "calcul itératif" (450), on prend en compte d'une part, ledit ensemble de paramètres (ai, b), et, d'autre part, un deuxième ensemble de données constitué d'un ensemble de données arbitraires dit ensemble d"'initialisation", et pour réaliser ledit calcul itératif: c1) on partitionne l'ensemble d'initialisation en sous-ensembles élémentaires (Ki); c2) on considère, pour chacun des sous-ensembles élémentaires (Ki) ainsi déterminés, les paramètres (ai, b) des applications élémentaires correspondant aux sous-ensembles élémentaires (Ki) résultant de la partition (204) de I'ensemble des données initiales, et on met en oeuvre lesdites applications élémentaires sur les sous-ensembles élémentaires (K,) de l'ensemble d'initialisation; c3) on réunit l'ensemble des données obtenues par les mises en oeuvre d'applications élémentaires, ceci constituant une représentation intermédiaire desdites données initiales, affectée d'une erreur; c4) on réitère les opérations c1-c3, en prenant comme ensemble d'initialisation, ladite représentation intermédiaire, de façon à s'approcher, du point fixe de ladite application contractive globale de second type; l'ensemble de données finalement obtenu constituant une représentation lo secondaire (RE) de l'ensemble de données ayant été l'objet d'un traitement primaire. 18. Procédé de traitement selon l'une quelconque des revendications 11 à 17 dans lequel lesdits paramètres (ai, b) sont une représentation primaire d'une image couleur à plusieurs composantes (Y, Cr, Cb) caractérisé en ce qu'après une étape de calcul itératif (450) sur l'une des composantes, on utilise comme initialisation de l'étape de calcul itératif (450) relatif à la composante suivante la représentation secondaire (RE) obtenue à la fin de l'étape de calcul itératif (450) relative à l'une des composantes déjà 2 o traitées. 19. Procédé de traitement selon l'une quelconque des revendications 11 à 18 dans lequel le traitement primaire étant effectué avec des premiers ensembles de données possédant des dimensions différentes pour les différentes composantes (Y, Cr, Cb) d'une image couleur, au cours de l'étape de "calcul itératif" (450) du traitement secondaire, on met en oeuvre des deuxièmes ensembles de données possédant les mêmes dimensions pour les différentes composantes. 20. Procédé de traitement selon la revendication 19, caractérisé en ce que lesdits premiers ensembles de données mis en oeuvre au cours du traitement primaire possèdent, pour chaque composante de chrominance (Cr, Cb) des dimensions inférieures aux dimensions du premier ensemble de données utilisé pour la composante de luminance (Y). 21. Dispositif de traitement primaire (500) d'un ensemble de données (IM, f) représentatives de grandeurs physiques, caractérisé en ce qu'il comporte: - des moyens d'entrée (591) de données à traiter; - des moyens de sortie (581) de paramètres (ai, b); - des premiers moyens logiques (510, 530, 531, 532, 550, 560, 570, 200 à 209), de transformation des données à traiter en paramètres (ai, b); lesdits premiers moyens logiques comportant: - des moyens de "construction" (510, 530, 531, 550, 560, 570, 205 à 209) adaptés: 1.5 à construire pour chaque ensemble de données à traiter, une application contractive globale, dite "de premier type" dont le point fixe constitue une approximation de tout ou partie de cet ensemble de données, en mettant en oeuvre, pour ce faire, au moins une application appartenant au groupe d'applications constitué des applications multidimensionnelles et des applications non linéaires, de façon: à permettre la mise en oeuvre d'une méthode d'approximations successives (407 à 410) convergeant vers le point fixe de ladite application contractive globale, sur un deuxième ensemble de données comportant un nombre de données différent du nombre de données du premier ensemble de données; 25. à délivrer sur les moyens de sortie (581) les paramètres (ai, b) ainsi déterminés, ceux-ci constituant une représentation primaire dudit ensemble de données. 22. Dispositif de traitement primaire selon la revendication 21, caractérisé en ce que lesdits premiers moyens logiques comportent des moyens de "partition" (204) adaptés à partitionner ledit ensemble de données en "sous-ensembles élémentaires" (Ki), et en ce que les moyens de construction sont adaptés, pour chacun des sous-ensembles élémentaires (Ki), à déterminer les paramètres (ai, b) de l'application élémentaire qui constitue une restriction de l'application contractive globale audit sous-ensemble élémentaire, pour permettre à la méthode d'approximations successives effectuée sur le deuxième ensemble de données, de prendre en compte des sous-ensembles élémentaires correspondant chacun à un sous-ensemble élémentaire du premier ensemble de données et comportant un nombre de données différent de celui du sous- ensemble élémentaire du premier ensemble de données correspondant. 23. Dispositif de traitement primaire selon la revendication 22, caractérisé en ce que les nombres de données des sous-ensembles élémentaires du deuxième ensemble de données sont proportionnels aux nombres de données des sous-ensembles élémentaires du premier ensemble de données qui leur correspondent respectivement. 24. Dispositif de traitement primaire selon l'une quelconque des revendications 22 ou 23, caractérisé en ce que les moyens de construction (205 à 209) sont adaptés à constituer ladite application contractive globale, partiellement au moins, d'applications Lipschitziennes élémentaires liées aux sous-ensembles élémentaires (Ki). 25. Dispositif de traitement primaire selon l'une quelconque des revendications 22 à 24, caractérisé en ce que les moyens de construction sont adaptés, pour chacun des sous-ensembles élémentaires (Ki) résultant de l'étape de partition initiale, à prendre en compte n sous-ensembles de travail (Li), n étant non nul, et à construire une application élémentaire des n sous- ensembles de travail dans le sous-ensemble élémentaire considéré en déterminant les paramètres (ai, b) de l'application élémentaire, pour permettre à la méthode d'approximations successives effectuée sur le deuxième ensemble de données, de prendre en compte des sous-ensembles de travail correspondant chacun à un sous-ensemble de travail du premier ensemble de données et comportant un nombre de données différent de celui du sous- ensemble de travail du premier ensemble de données correspondant. 26. Dispositif de traitement primaire selon la revendication 25, caractérisé en ce que les nombres de données des sous-ensembles de travail du deuxième ensemble de données sont proportionnels aux nombres de données des sous-ensembles de travail du premier ensemble de données qui lo leur correspondent respectivement. 27. Dispositif de traitement primaire selon l'une quelconque des revendications 25 ou 26, caractérisé en ce qu'au moins un sous-ensemble de travail (Li) parmi lesdits n sous-ensembles de travail (Li) présente une taille différente de celle du sous-ensemble élémentaire correspondant. 28. Dispositif de traitement primaire selon l'une quelconque des revendications 25 à 27, caractérisé en ce que les moyens de construction (510, 530, 531, 550, 570, 205 à 209) sont adaptés à prendre en compte pour chacun desdits sous-ensembles élémentaires (Ki), n sous-ensembles de travail (Li) géométriquement liés de manière prédéterminée au sousensemble élémentaire considéré. 29. Dispositif de traitement primaire selon l'une quelconque des revendications 21 à 28, caractérisé en ce que les premiers moyens logiques comportent: -un microprocesseur (550) - une mémoire morte (530) comportant un programme de représentation primaire de données, et - une mémoire vive (510) comportant des registres adaptés à enregistrer des variables modifiées au cours de l'exécution dudit programme de représentation primaire. 30. Dispositif de traitement primaire selon l'une quelconque des revendications 21 à 29, caractérisé en ce que, ledit premier ensemble de données constituant une image couleur à plusieurs composantes (Y, Cr, Cb), il comporte un moyen de filtrage, adapté à filtrer au moins l'une desdites composantes et à modifier le nombre de données de l'ensemble de données liés à au moins une composante. 31. Dispositif de traitement primaire selon la revendication 30, caractérisé en ce que, le moyen de filtrage est adapté à réduire le nombre de données des ensembles de données liés aux composantes de chrominance, par rapport au nombre de données de l'ensemble de données lié à la composante de luminance. 32. Dispositif de traitement secondaire (600) de données convenant pour le traitement secondaire d'un ensemble de paramètres (ai, b) résultant du traitement primaire d'un premier ensemble de données initiales (IM, f), caractérisé en ce qu'il comporte: - des moyens d'entrée (681) des paramètres (ai, b); - des moyens de sortie (691) de données; - des seconds moyens logiques (560, 570, 610, 630, 631, 632, 650, 400 à 410) de transformation des paramètres (ai, b) en données; lesdits seconds moyens logiques comportant: - des moyens de "calcul itératif" (650, 450) adaptés a mettre en oeuvre, à partir dudit ensemble de paramètres (ai, b) et d'un deuxième ensemble de données possédant un nombre de données différent du nombre de données du premier ensemble de données, une méthode d'approximations successives (407 à 410) convergeant vers le point fixe d'une application contractive globale dite de "second type", ce dernier point fixe constituant une représentation secondaire (RE) dudit ensemble des données initiales, et à délivrer les données ainsi calculées sur les moyens de sortie. 33. Dispositif de traitement secondaire selon la revendication 32, caractérisé en ce qu'il comporte des moyens de partition adaptés à partitionner ledit deuxième ensemble de données en sous-ensembles dits "élémentaires" correspondants chacun à un sous-ensemble élémentaire du premier ensemble de données et comportant un nombre de données différent de celui du sous- 1o ensemble élémentaire du premier ensemble de données correspondant. 34. Dispositif de traitement secondaire selon la revendication 33, caractérisé en ce que les nombres de données des sous-ensembles élémentaires du deuxième ensemble de données sont proportionnels aux nombres de données des sous-ensembles élémentaires du premier ensemble de données qui leur correspondent respectivement. 35. Dispositif de traitement secondaire selon l'une quelconque des revendications 33 ou 34, caractérisé en ce que, le traitement primaire prenant en compte, pour chacun des sous-ensembles élémentaires (Ki) résultant de la dernière étape de partition, n sous-ensembles de travail (Li), n étant non nul, et construisant une application élémentaire des n sous- ensembles de travail dans le sous-ensemble élémentaire considéré en déterminant les paramètres (a,, b) de l'application élémentaire, - il prend en compte des sous-ensembles de travail correspondant chacun à un sous-ensemble de travail du deuxième ensemble de données et comportant un nombre de données différent de celui du sous-ensemble de travail du premier ensemble de données correspondant. 36. Dispositif de traitement secondaire selon la revendication 35, caractérisé en ce que les nombres de données des sous-ensembles de travail du deuxième ensemble de données sont proportionnels aux nombres de données des sous-ensembles de travail du premier ensemble de données qui leur correspondent respectivement. 37. Dispositif de traitement d'un ensemble de données (IM, f) représentatives de grandeurs physiques, comprenant un dispositif de traitement primaire selon l'une quelconque des revendications 20 à 30 et un dispositif de traitement secondaire selon l'une quelconque des revendications 31 à 35, le dispositif de traitement secondaire étant adapté à traiter l'ensemble 0o de paramètres (ai, b) sortant du dispositif de traitement primaire. 38. Dispositif de traitement selon l'une quelconque des revendications 32 à 37, caractérisé en ce que lesdits seconds moyens logiques comportent: - un microprocesseur (650), - une mémoire morte (630) comportant un programme de représentation secondaire (RE) de données, et une mémoire vive (610) comportant des registres adaptés à enregistrer des variables modifiées au cours de l'exécution dudit programme de représentation secondaire (RE). 39. Dispositif de traitement selon l'une quelconque des revendications 32 à 38, dans lequel lesdits paramètres (ai, b) sont une représentation primaire d'une image couleur à plusieurs composantes (Y, Cb, Cr) caractérisé en ce que lesdits seconds moyens logiques sont adaptés à fournir au moyen de calcul itératif (450) un ensemble d'initialisation arbitraire pour le traitement des paramètres (ai, b) relatifs à une première composante et, comme ensemble d'initialisation pour le traitement des composantes suivantes, la représentation secondaire (RE) obtenue à la fin de l'étape de calcul itératif (450) relative à l'une des composantes déjà traitées. 40. Dispositif de traitement secondaire selon l'une quelconque des revendications 32 à 38, dans lequel lesdits paramètres (a,, b) sont une représentation primaire d'une image couleur à plusieurs composantes (Y, Cr, Cb) caractérisé en ce que lesdits seconds moyens logiques sont adaptés à fournir aux moyens de calcul itératif (450) un ensemble d'initialisation constitué de données arbitraires pour le traitement des paramètres (a,, b) relatifs à chaque composante et en ce que les moyens de calcul itératif sont adaptés à mettre en oeuvre en parallèle une méthode d'approximations successives (407 à 410) convergeant vers le point fixe de chaque application contractive globale l0 correspondant à une composante. 41. Dispositif de traitement selon l'une quelconque des revendications 32 à 40 dans lequel le traitement primaire étant effectué avec des premiers ensembles de données possédant des dimensions différentes pour les différentes composantes (Y, Cr, Cb) d'une image couleur, les moyens de calcul itératif mettent en oeuvre des deuxièmes ensembles de données possédant les mêmes dimensions pour les différentes composantes. dans lequel lesdits paramètres (ai, b) sont une représentation primaire d'une image couleur à plusieurs composantes (Y, Cr, Cb) lesdits premiers ensembles de données étant de dimensions différentes pour lesdites composantes, caractérisé en ce qu'il met en oeuvre desdits deuxièmes ensembles de données de dimensions égales pour lesdites composantes. 42. Dispositif de traitement selon la revendication 41, caractérisé en ce que lesdits premiers ensembles de données mis en oeuvre au cours du traitement primaire possèdent, pour chaque composante de chrominance (Cr, Cb) des dimensions inférieures aux dimensions du premier ensemble de données utilisé pour la composante de luminance (Y). 43. Imprimante, caractérisée en ce qu'elle comporte un dispositif de traitement selon l'une quelconque des revendications 32 à 42. 44. Photocopieur, caractérisée en ce qu'elle comporte un dispositif de traitement selon l'une quelconque des revendications 32 à 42. 45. Appareil photographique, caractérisé en ce qu'il comporte un dispositif de traitement selon l'une quelconque des revendications 21 à 42. 46. Dispositif de télécopie, caractérisé en ce qu'il comporte un dispositif de traitement selon l'une quelconque des revendications 21 à 42.