Brevet FR2639740A1 — publié le 19900601

La présente invention concerne un traitement d'image pour diagnostic médical et plus particulièrement un traitement d'image qui conduit à traiter un zoom d'images. Il est souvent désirable d'agrandir la dimension d'une partie d'une image affichée dans les systèmes de production d'images pour diagnostic médical. Une image zoom affichée est constituée d'un certain nombre de pixels qui sont contigus dans une image initiale maintenant située dans l'image zoom avec de nouveaux pixels entre les pixels initiaux. Par exemple, une région d'intérêt dite ROI est choisie constituée d'un nombre donné de pixels contigus. La région d'intéret est agrandie pour grandir l'écran d'affichage tout entier. Pour réaliser ceci, les pixels contigus initiaux de la région d'intérit sont séparés par de nouveaux pixels places entre les pixels initiaux. Les nouveaux pixels places entre les pixels initiaux sont prévus avec des données d'intensité obtenues par interpolation en utilisant les valeurs d'intensité des pixels initiaux. En général, l'interpolation est réalisée au moyen de techniques telles que l'extraction des valeurs d'intensité de pixels connus, disons deux pixels entourant un pixel inconnu, en ajoutant arithmétiquement les valeurs d'intensité et en les divisant par deux, dans le cas o il existe seulement un pixel inconnu entre deux pixels connus. Lorsqu'il n'existe plus qu'un pixel inconnu entre le pixel de valeur connue, alors des facteurs de pondération sont utilisés pour rendre compte des distances entre les pixels inconnus et les pixels connus. Un zoom est généralement réalisé lorsque l'on désire regarder une ROI particulière et afficher celle-ci sur l'écran entier. Dans le passé, un zoom a été réalisé principalement par des processus d'interpolation linéaire qui comprenne une interpolation arithmétique décrite précédemment ou des processus d'interpolation analytique. Dans une interpolation analytique, la valeur interpolée I, doit Itre trouvée en résolvant l'équation linéaire: f(I) = mX + b, o m est égal à la pente de la fonction linéaire, b est la valeur de la fonction lorsque X est égal à zéro, et o X est l'emplacement du pixel dans la direction horizontale. Ainsi, la connaissance des valeurs d'intensité de deux pixels quelconques fournit deux points sur un graphique d'intensité par rapport & l'emplacement. Les deux points définissent une-ligne droite qui présente une pente m et coupe l'axe d'intensité au point b. L'interpolation linéaire du cas bidimensionnel (production d'images) assure l'ajustement linéaire à la fois dans les directions X et Y. Un problème spécifique se produisant par l'emploi de l'interpolation linéaire est que dans la vie réelle il existe des emplacements o la pente m se modifie de manière critique - la pente peut, par exemples s'inverser. Ainsi, les deux pixels connus suivants doivent définir une ligne ayant une pente & l'oppose de la pente de la fonction linéaire initiale. La modification obtenue dans la pente de la ligne entraîne des artefacts tels que des artefacts en escalier dans l'image affichée. Il existe à présent des algorithmes qui lissent les résultats des changements profonds dans la direction de la pente qui provoque des angles aigus dans la fonction linéaire. Toutefois, ces algorithmes tendent en général à tâcher l'image et il est souvent douteux de savoir si le remède de l'artefact n'est pas fatal à la qualité de l'image. Un procédé conseillé pour améliorer les images pour lesquelles un zoom a été réalise consiste à effectuer l'interpolation des pixels en utilisant des fonctions polynomiales. ou, plus généralement, des fonctions binomiales au lieu de la fonction linéaire. Toutefois, le traitement pour une interpolation polynomiale devient très complexe et en général requiert beaucoup trop de temnps et d'espace mémoire du fait que dans une interpolation binomiale et polynomiale il existe de plus nombreuses étapes arithmétiques que dans une interpolation linéaire. Par suite, les spécialistes de la technique sont toujours à la recherche de procédés d'interpolations utiles pour des images zoom. Cette recherche est particulièrement intense en ce qui concerne la découverte d'un procédé d'interpolation qui permettra d'effectuer un "zoom" "en direct", c'est-à-dire la réalisation suffisamment rapide des interpolations pour traiter les données d'image entrantes durant l'acquisition, ou peut permettre à l'opérateur de faire varier interactivement le zoom à des vitesses dépassant 10 cadres par seconde. En conformité avec la présente invention, un procédé de traitement d'une image normale en entrée pour fournir une image zoom en sortie est énonce. L'image en entrée comprend une matrice d'entrée de rangées et de colonnes de "régions" ayant des valeurs d'intensité en entrée correspondant aux rangées et colonnes de pixels dans une image affichée normale. L'image zoom en sortie comprend une matrice de rangées et de colonnes de régions ayant des valeurs d'intensité correspondantes aux rangées et aux colonnes de pixels dans une image zoom affichée (ici les "régions" sont souvent également désignées comme "pixels") ledit procédé comprenant les opérations consistant à: spécifier un facteur de zoom et une région d'intért contenant des valeurs d'intensité en entrée par lesquelles un zoom doit Etre réalisé, transférer les valeurs d'intensité en entrée des régions spécifiées pour lesquelles un zoom doit Etre réalise aux régions désignées de la matrice zoom. lesdites régions désignées avec les valeurs dintensité transférées connues séparées par des régions de valeurs d'intensité connue, convoluer les valeurs d'intensité connue incluant les valeurs d'intensité voisines de celles des régions désignées de la matrice zoom pour obtenir des valeurs d'intensité pour les régions de valeurs d'intensité connue, et ladite opération de convolution incluant l'opération consistant à choisir un kernel mémorise individuel à chaque région de la matrice zoom, ledit kernel étant sélectionné comme une fonction de l'emplacement de la rangée et de la colonne de la région de valeurs d'intensité inconnue et étant alimenté avec des valeurs d'intensité connue pour la convolution comme une fonction du zoom, permettant que les valeurs d'intensité soient appliquées à chaque région de la matrice zoom. Une caractéristique de l'invention comprend une disposition d'une pluralité de kernels mémorises sous la forme de tables de consultation individuelles à chaque région de la matrice zoom en fonction de la rangée et de la colonne de la région. Une caractéristique supplémentaire de l'invention comprend l'étape consistant à diviser la matrice zoom en sous-divisions entre des régions d'intensité connue, chaque sous-division fournissant un pointeur aux tables de consultation, les tables de consultation ainsi désignées étant convoluées avec des régions voisines d'intensité connue et dans laquelle les régions d'intensité connues utilisées dans la convolution sont incrémentées comme une fonction du nombre de sous- divisions divise par le facteur de zoom. Une caractéristique encore supplémentaire de l'invention comprend l'opération consistant à déterminer une valeur d'incrément basée sur le nombre de sous- divisions divisé par le facteur de zoom chaque fois que ladite valeur d'incrément est utilisée pour 'avancer pas & pas" à travers les tables de consultation de kernels à convolution mémorisée. Une caractéristique supplémentaire de l'invention comprend l'utilisation des kernels sur la base de polynSmes. Une autre caractéristique de l'invention utilise les kernels sur la base de binômes. Une autre caractéristique de l'invention est qu'il est nécessaire de trouver les valeurs pour les tables de consultation seulement une fois pour chaque système de production d'image du fait que les kernels utilisés dans le processus de convolution sont indépendantes du format d'image, du facteur de zoom et du format de la matrice. Ainsi, les kernels sont calculés une fois et sont mémorisés sur un disque pour Etre charges dans la mémoire chaque fois que la fonction de zoom est appelée. Par suite, un zoom d'image interactif est pratique avec ce système. Une caractéristique supplémentaire en relation avec l'invention est l'emploi de l'opération de convolution caractérisant les kernels obtenus à partir d'un polynôme qui sont fonction seulement de l'emplacement de la région nécessitant une valeur interpolée en liaison avec tout traitement d'image quelconque qui pourrait nécessiter une interpolation, au lieu de l'interpolation linéaire bien connue, afin d'ainsi améliorer la qualité d'image en minimisant les artefacts tels que les artefacts d'étape. Encore une autre caractéristiques de l'invention a utilisé le procédé d'interpolation décrit pour translater une image de format matriciel (c'est-à-dire 340 x 340) en une autre image de format matriciel (c'est-àdire 512 x 512) avec un zoom égal à un (aucun zoom n'est effectué). Les caractéristiques et buts mentionnés ci-dessus et autres de l'invention seront le mieux compris lorsqu'ils seront considérés à la lumière de la description suivante d'un aspect principal de la présente invention prise en liaison avec les dessins suivants, dans lesquels: Les fig. la, lb sont des schémas synoptiques illustratifs de deux systèmes d'interpolation linéaire de la technique antérieure. La fig. lc est un graphique représentatif d'un procédé d'interpolation linéaire analytique. La fig. 2 est un schéma synoptique représentatif du présent système d'interpolation inventif utilisé pour un zoom. La fig. Z est une représentation graphique d'une fonction de polynômes définie par quatre valeurs d'intensité connue en adaptant au mieux une courbe aux: valeurs connues pour permettre la détermination de valeurs inconnues aux emplacements sélectionnés. La fig. 4 représente une table de consultation (:ernel) utilisée avec'l'interpolation polynomiale à convolutive de la présente invention, La fig. 5 est une représentation de matrice zoom représentant des grilles entre des régions de valeurs d'intensité connue recouvrant les régions devant Etre pourvues de valeurs d'intensité en conformité avec l'invention, et Les fig. 6a,6b,6c et 6d représentent un seul et mme pro- gramme d'interpolation effective zoom pr ordinateur en langage Fortram et se font suite les unes aux autres. Le zoom est une procédure de traitement d'image bien connue dans laquelle une région d'intérêt est agrandie, telle que par exemple une région plus petite que l'image affichée entière répartie sur l'image affichée complète. Pour réaliser l'agrandissement, des valeurs d'intensité connue des régions pour lesquelles un zoom a été réalisé sont interpolées pour remplir les valeurs d'intensité des nouvelles régions insérées entre les régions de valeurs d'intensité connue. Plus particulièrement, la donnée d'image en entrée acquise est disposée sous forme matricielle, c'est-àa- dire en rangées et en colonnes de S régions contenant des valeurs d'intensité mesurées avant le traitement de zoom. Pour effectuer un zoom, une région d'intérgt de la matrice d'image en entrée est traitée en plaçant symétriquement les valeurs d'intensité connue de la région matricielle dans la région d'intért sur une matrice agrandie comportant beaucoup plus de régions matricielles qu'il n'en est inclus dans la région deintérQt. Ceci laisse des régions vierges dans les rangées et colonnes entre les régions d'intensité connues prises à partir de la région d'intérêt de l'image en entrée. Dans la technique antérieure, les régions ayant des intensités inconnues sont alimentées avec des valeurs d'intensité par interpolation "linéaire" entre les régions de valeurs d'intensité connue. Il existe de nombreuses manières de réaliser une interpolation linéaire. Une manière est représentée à la fig. la qui peut Etre appelée un procédé analytique. Plus particulièrement, à la fig. la il est représenté une disposition de zoom en 11. La disposition de zoom comprend une image en entrée 12 inclut des rangées il, i2, i3,... in et des colonnes jl, j2, j3. in. En général, de telles matrices peuvent comprendre 512 rangées par 512 colonnes. Pour effectuer un zoom de l'image, l'opérateur sélectionne une région d'intért sur un dispositif d'entrée par opérateur tel qu'une unité de clavier à boule roulante 13 par exemple. En variante, l'opérateur peut sélectionner un point central tel que i3, j3 et un facteur de zoom, disons 8. Le processeur central du système détermine la dimension de la région d'intérêt (ROI), c'est-à-dire quelle quantité de la région matricielle doit doit donner lieu à un zoom. Les valeurs d'intensité des pixels dans la ROI sont extraites de la donnée d'image en entrée par l'unité d'extraction..DTD: 14 et écrites dans un opérateur d'interpolation 16. L'opérateur d'interpolation effectue ensuite une interpolation en utilisant une interpolation linéaire entre les pixels sélectionnés d'une matrice 12 pour fournir des valeurs d'intensité pour les pixels entre les pixels de la ROI qui sont extraits de l'image en entrée. Une unité d'écriture 17 transfère les valeurs ici d'intensité dans les pixels d'une image en sortie 18 qui sont extraits de la ROI de l'image en entrée 12. L'interpolation est effectuée mathématiquement. Ainsi, par exemple, s'il existe une région vierge entre un pixel ayant une valeur d'intensité de 5 et un pixel ayant une valeur d'intensité de 3, alors la valeur interpolée de la région entre celle-ci sera de 5+3 divisé par 2 ou par 4. S'il existe deux régions vides entre les régions ayant des intensités connues, alors un facteur de pondération peut Etre utilisé pour rendre compte de la distance d'une région vierge depuis la région ayant l'intensité connue. Si la région vierge tombe entre à la fois des rangées et des colonnes de pixels connus, l'interpolation linéaire est effectuée par pondération en conformité avec les distances entre quatre pixels voisins connus. De cette manière, l'image zoom en 18 est constituée de régions qui ont toutes des valeurs d'intensité pour fournir une image complète. Dans la disposition de la technique antérieure de la fig. lb, une solution analytique est utilisée. Ici de 350 nouveau, il existe une image en entrée 12a. Un dispositif d'entrée par opérateur tel qu'un clavier à boule roulante 13a permet la sélection d'un point central de région d'intért et d'un facteur de zoom. L'unité d'extraction 14a fonctionne en réponse aux informations reçues depuis la région d'intérît sélectionnée et extrait les données de la région d'intérêt sélectionnée de l'image en entrée 12a. Seules deux régions ayant des intensités connues doivent tre utilisées pour développer l'équation linéaire f(I) = mX + b (représenté à la fig. lc o X est la distance à laquelle se trouve la région du coté gauche de la matrice, m est la pente de la ligne droite définie par les deux points (il, jl) et (i3, j3) et o b est la valeur d'intensité de la ligne lorsque X est égal à zéro, en d'autres termes, lorsque la ligne coupe l'axe d'intensité I. A partir des deux points de donnée, la pente m peut Etre déterminée et b peut également Etre déterminé. Le graphique f(I) de la fig. lc représente l'intensité I par rapport à la distance horizontale X à partir de l'axe d'intensité I. Toute valeur quelconque sur la fonction peut être déterminée à partir de la fonction. Une approche similaire est utilisée pour déterminer l'intensité I par rapport à la distance verticale Y. Un problème présenté par l'interpolation linéaire est qu'elle conduit souvent à des artefacts connus comme artefacts "mixtes". Les artefacts mixtes apparaissent lorsque se produisent des discontinuités dans les pentes. Par exemple, lorsque sur le c8té droit d'un pixel connu la pente est positive et que sur le c8té gauche la pente est négative. Par suite, un procédé d'interpolation utilisant une analyse d'un ordre supérieur qui fournit une pente variant en continu est présenté. Le procédé d'interpolation est représenté à titre d'exemple dans le schéma synoptique de la fig. 2. Dans le passe, des interpolations d'un ordre supérieur n'ont pas été pratiquées du fait du grand nombre de calculs nécessaires exerçant une contrainte sur l'ordinateur et sur la mémoire du système. Toutefois, le système décrit ici qui ) effectue une convolution afin d'interpoler en utilisant des kernels indépendants de la dimension d'image, du facteur de zoom et des intensités rend possible des interpolations d'ordre supérieur. Le système 20 comprend une image en entrée d'une donnée acquise représentée en 21. Une donnée peut tre acquise par l'un quelconque des modes d'acquisition utilisé par le champ médical de diagnostic, par exemple. L'image en entrée est de nouveau constituée d'une matrice de rangées et de colonnes qui ne sont pas représentées. Des moyens tels qu'un clavier à boule roulante 22 sont prévus pour une interaction de l'opérateur avec l'image en entrée que l'opérateur utilise pour sélectionner la région d'intérêt ou le centre de la région d'intérêt plus un facteur de zoom. En mode interactif, l'image en sortie est présentée presque instantanément sur l'écran et l'opérateur peut modifier le zoom, faire un panoramique et un défilement direct. L'unité 22 amène une unité d'extraction 23 à extraire des valeurs d'intensité des régions appropriées de l'image en entrée. Les régions appropriées sont les pixels voisins nécessaires pour le processus de convolution. Dans un exemple polynomial, quatre valeurs horizontales et quatre valeurs régionales verticales voisines (c'est-à-dire un minimum de 16 pixels) sont extraites. L'information extraite est fournie & une mémoire 25 dans un processeur de convolution de zoom 24. C'est la fonction du processeur de convolution 24 de fournir des valeurs d'intensité aux régions matricielles entre les régions X:0 extraites qui font partie de l'image zoom en sortie. Ceci est effectué en utilisant les kernels sélectionnés d'une pluralité de kernels mémorisés dans le processeur de convolution de zoom à la mémoire de kernels 26. La mémoire de kernels est représentée en plus grand -;5 détail en 27 bien que dans un mode de réalisation préféré la mémoire de kernels constitue une partie intégrante du processeur de convolution de zoom. La pluralité de kernels dans un mode de réalisation préféré est établie en une pluralité de tables de consultation telles que la table de consultation 28. Des moyens sont prévus pour produire des tables de consultation qui contiennent les "poids" devant tre utilisés dans la convolution avec les valeurs de pixels connues pour fournir les valeurs de pixels manquantes. Dans un mode de réalisation préféré, les tables de consultation (kernel) sont indépendantes des données d'image. du facteur de zoom et du format matriciel. Des moyens sont prévus pour produire de telles tables de consultation indépendantes et pour sélectionner la table de consultation appropriée devant Etre utilisée dans l'opération de convolution. Les moyens impliquent l'établissement d'un système de grilles indépendantes s'étendant entre les régions d'intensité connues pour lesquelles un zoom a été réalisé. Dans un mode de réalisation préféré, la grille comprend 16 lignes horizontales et 16 lignes verticales entre chacune des régions de la ROI qui sont pour lesquelles un zoom a été réalisé. Le nombre de lignes de grille est généralement indiqué ici par isubX et isubY. Les poids pour les tables de consultation sont déterminés à l'aide d'une table de poids prévue pour chaque intersection de ligne de grille. Ainsi, la matrice en entrée est développée pour fournir un nombre de lignes sélectionnées entre chaque zoné d'intensité connue transférée vers la matrice zoom en sortie. Dans un mode de réalisation, 16 lignes sont utilisées s'étendant depuis le milieu de l'un des pixels d'intensité connue vers le milieu du pixel suivant d'intensité connue de la matrice zoom en sortie. En fonction du facteur de zoom, il existe un nombre fixe de régions qui sont insérées entre les régions d'intensité connues. L'intersection des lignes de grille les plus proches du centre de région insérée est le pointeur vers (ou'sélectionne) le kernel désiré devant Utre utilisé, le kernel est une table de consultation de 16 poids (dans notre exemple) qui est convoluée avec 16 pixels voisins d'intensité connue pour fournir la valeur d'intensité de la région insérée pour la matrice en sortie. Les poids de chaque intersection de la grille sont IC) précalculés comme indiqué à la fig. 3. On suppose qu'il existe 4 points de valeur connue tels que f1, f2, f3 et f4 entre lesquels les valeurs interpolées telles que la valeur d'intensité du point X doit Etre déterminée située entre les points f2 et f3 (une interpolation dimensionnelle est caractérisée). Une fonction f(X) reliant les points fi, f2, f3 et f4 peut Etre approchée par un polyn8me du troisième ordre f (X) p(X) = f1 + rjfl + r(r-1l)/2_t2f, + [r'r-, simplement un kernel à convolution dimensionnelle. Du fait que le zoom graphique implique une interpolation bidimensionnelle, des facteurs de pondération sont obtenus de manière identique dans la direction verticale. Le résultat est un kernel de 2) convolution 4 x 4 représenté à la fig. 4 construit en utilisant: w(isubx, isuby) = wij(isubx, isuby) = wxi(isubx) * wyi(isuby) Il conviendra de noter qu'il existe ISUB*ISUB de tels kernels de convolution, dans lequel chacun est utilise pour trouver la valeur à une intersection de grille spécifique par exemple, si la valeur du point interpolé (intersection de grille) se produisant en isubx = 3 et isuby = 10 est nécessaire pour que l'image en sortie du kernel approprié soit trouvée wij(3,10) et soit convolué avec les 16 zones matricielles voisines de la matrice d'entrée pour produire la matrice zoom en sortie. Avec ISUB égal à 16, il existe 4K mots nécessaires pour mémoriser les kernels. De cette manière, les tables de consultation (kernels de convolution) fournissent des poids pour une interpolation par convolution de la valeur de zone connue afin d'obtenir des valeurs pour le nombre de zoom situées entre elles comme imposées par le facteur de zoom. Une convolution avec les tables de consultation sur la donnée provenant de l'image en entrée fournit des données pour l'image zoom en sortie représentée en 31. Les tables de consultation sont choisies par l'emplacement de l'intersection de ligne de grille la plus proche du centre de la région pour la région recherchant une valeur d'intensité. Afin de faciliter la localisation de l'intersection de ligne de grille la plus proche, une valeur d'incrément est sélectionnée par une unité de sous-incrémentation 29. Cette valeur est ajoutée au pointeur de la table de consultation pour chaque nouvelle région ayant une valeur générée. Une ligne partant de l'unité 22 est représentée se dirigeant vers l'unité de lecture 23 pour extraire les y valeurs de zone ROI sélectionnées à partir de la matrice d'image en entrée 21. Les zones lues sont exploitées par les kernels qui sont sélectionnés en conformité avec le facteur de zoom déterminant les 16 zones voisines de valeur connue qui sont utilisées dans la convolution. Ainsi, une ligne provenant de l'unité de sous- incrémentation 29 est tracée jusqu'au processeur de convolution de zoom pour déterminer quelles zones de valeurs connues doivent Stre convoluées avec le kernel. Par exemple, à la fig. 5, les régions entre les zones O0 initiales de valeurs connues (f2,2, f2,3, f3,2 et f3,3) sont sous-divisées en une grille fixe. La grille fixe dans un mode de réalisation préféré est une grille 16 x 16 commençant au début d'une zone connue et s'étendant vers le commencement de la zone connue suivante. Il convient de comprendre que d'autres schémas de grilles cycliques peuvent également etre utilises dans le cadre de l'invention. Les valeurs des intersections de grille sont fournies en utilisant les kernels pour chaque intersection. Le processus de convolution commence aux positions de départ X, Y définies comme: Départ X = centre X - (format matriciel/2)(1/zoom) Départ Y = centre Y (format matriciel/2)(1/zoom) le centre X, le centre Y et le facteur de zoom sont définis par l'utilisateur. La dimension de matrice est généralement celle de l'image en entrée, mais peut tre d'autres formats. Les pointeurs de sous-divisions (intersection de grille) commence en subx = suby = 1 et sont incrémentés par: subinc = ISUB/zoom dans les directions X et Y à mesure que progresse la convolution. Ainsi, lorsque les points de sous-divisions deviennent supérieurs à ISUB), alors le kernel est décalé et le reste (ISUB-SUBx ou ISUB-SUBy est utilisé comme le nouveau pointeur de sous-divisions). Le processus de convolution peut Etre effectué soit rangée par rangée soit colonne par colonne. Toutefois, le processus devra Etre cyclique. En d'autres termes, si la convolution est effectuée rangée par rangée, quatre rangées doivent tre mémorisées dans la mémoire du processeur 25 pour fournir la valeur d'intensité des quatre zones de rangées voisines nécessaire pour fournir la valeur de chaque rangée unique insérée dans l'image zoom en sortie. Si l'une des rangées n'est plus désormais nécessaire, i.e. lorsque isubY est supérieur à ISUB, seule une rangée a besoin d'Ptre remplacée par la rangée de l'image en entrée suivante dans la mémoire 25 et le processus se poursuit sans nécessiter qu'un réseau entièrement nouveau soit entre en mémoire 25. Par exemple, à la fig. 5 f2, 25 f2,3, f32 et f3,3 sont des zones ayant des valeurs d'intensité connues transférées en provenance de l'image en entrée. Le facteur de zoom est 4 et ISUB est 16. Il existe trois zones nécessitant des valeurs d'interpolation entre f2,2 et f2,3 et entre f3,2 et f3,3 de mime que les trois ensembles de cinq zones ci-après f2,2, f2,3 et au-dessus f3,2 et f3s,53- En supposant que l'on effectue une intersection t rangée par rangée, alors la mémoire 25 aura quatre colonnes et quatre rangées provenant de l'imageinitiale qui y est mémorisée pour fournir les 16 valeurs d'intensité nécessaires des zones voisines requises pour le processus de convolution. Ainsi, afin d'obtenir la i5 valeur de zone A, le kernel pour l'intersection de grille 19, 24 est utilisé en même temps que les valeurs de la région à quatre rangées voisines et à quatre colonnes voisines: f1, 1, f1,2, fl,35 fl,4 - {f2,1' f2,2, f2,35, f2,4 f3,15 f3,25 f3,35 f3,4 f4,1, f4,2, f4,3, et f4,4. De manière similaire, pour obtenir la valeur de la zone B, le kernel pour l'intersection de grille 27, 34 est utilisé en mime temps que les valeurs des mimes 16 zones voisines. Afin d'obtenir la valeur d'intensité de zone i5, j3, le kernel pour l'intersection de grille la plus proche du centre de la région est utilisé en mime temps que les valeurs des régions connues nouvellement incrémentées: f2,1- f2,2, f2,35, f2,4 f3,1, f3,2, f 3,3, f3,4 f4,1, f4,25 f4,3, f4,4 f5,1, f5,2, f5,3, et f5,4. Ainsi, la cinquième rangée a été prise et la première a été rejetée. Le programme d'interpolation effective de zoom par ordinateur en langage Fortram est illustré aux figures 6a, 6b, 6c et 6d. Ainsi, un nouveau et unique système d'interpolation de zoom est créé lequel, grâce aux kernels qui sont seulement fonction d'un emplacement de l'interpolation, permet une convolution destinée à obtenir des valeurs d'intensité sur une zone insérée en utilisant des tables de consultation. Ceci conduit à une interpolation très rapide qui assure une qualité d'image comparable à la reconstruction avec un zoom et fournit efficacement une production d'image zoom en "temps réel". Tandis que l'invention a été décrite en utilisant un mode de réalisation à titre d'exemple, il convient de comprendre que ce mode de réalisation est décrit à titre d'exemple seulement et non comme limitation de la portée de l'invention. REVENDICATIONS: 1. Procéde de traitement d'une image normale en entrée (21) pour fournir une image zoomée en sortie (31), l'image en entrée (21) comprenant une matrice d'entrée (28) de rangées et de colonnes de régions ayant des valeurs d'intensité en entrée correspondant aux rangées et colonnes de pixels dans une image affichée normale. ladite image zoom en sortie (31) comprenant une matrice zoom (28) de rangées et de colonnes de régions ayant des valeurs d'intensité correspondant aux rangées et colonnes de pixels dans une image zoom affichée (31), ledit procédé comprenant: la spécification d'un facteur de zoom (15) et des régions de valeurs d'intensité en entrée pour lesquelles un zoom doit être réalisé, le transfert des valeurs d'intensité en entrée aux régions spécifiées pour lesquelles un zoom doit #tre réalisé vers les régions désignées de la matrice zoom (28), lesdites régions désignées étant séparées par des t:0 régions de valeurs d'intensité connue, la sélection des kernels mémorises individuels pour chacune desdites régions d'intensité connue de la matrice zoom (28), lesdits kernels mémorisés étant sélectionnés en fonction de l'emplacement de la région d'intensité inconnue de la matrice zoom (28), la sélection de régions voisines & partir d'emplacements déterminés par l'emplacement de la région d'intensité.inconnue et le facteur de zoom (15), et la convolution desdits kernels mémorisés sélectionnés avec des valeurs d'intensité desdites régions voisines sélectionnées de valeurs d'intensité connue incluant les régions désignées de la matrice zoom d'intensité inconnue. 2. Procédé selon la revendication 1, comprenant un moyen pour sélectionner 16 régions voisines afin de satisfaire une fonction polynomiale. Procédé selon la revendication 1, dans laquelle ladite étape de sélection comprend la sélection de 9 régions voisines pour satisfaire une fonction binomiale. 4. Procédé selon la revendication 1 dans laquelle l'étape de sélection desdits kernels mémorisés comprend les opérations consistant à: recouvrir une grille sur la matrice zoom (28), ladite grille comportant une pluralité de lignes horizontales entre lesdites régions d'intensité connue dans lesdites rangées et une pluralité de lignes verticales entre lesdites régions d'intensité connue dans lesdites colonnes, un kernel se trouvant associé à chaque intersection desdites lignes horizontales et verticales, et sélectionner ledit kernel le plus près du centre de ladite région d'intensité inconnue pour Etre utilisé dans l'étape de convolution (24) afin d'obtenir une valeur d'intensité pour ladite région d'intensité inconnue. 5. Procédé selon la revendication 4, comprenant l'étape consistant à utiliser une grille 16 x 16. 6. Procédé selon la revendication 1, dans lequel lesdits kernels sont indépendants des données d'image, du facteur de zoom (15) et du format matriciel. 7. Procédé selon la revendication 4, comprenant l'étape consistant à incrémenter les kernels sur la base du nombre de lignes de grille dans lesdites régions spécifiées de valeurs d'intensité connue divisées par le 0 facteur de zoom (15) de façon que lorsqu'il existe 16 lignes de grille et un facteur de zoom de 8 seul chaque second kernel soit utilisé. 8. Procédé selon la revendication 1, dans lequel lesdits kernel sont basés sur une fonction polynomiale. S5 9. Procédé selon la revendication 1, dans lequel lesdits kernels sont basés sur une fonction binomiale. li. Procédé selon la revendication 1. comprenant les étapes consistant à produire un kernel, lesdites étapes de production deun kernel comprenant les opérations consistant à: a) établir une fonction polynomiale qui réunit les régions spécifiée d'intensité connue en une rangée comme une fonction des lignes de grille entre lesdites régions spécifiées, ladite convolution étant: P(isubX) = fi (i-r + A/2 - B/6) + f2 (r-A + B/2) + fE (A/2-B/2) + f4 (B/6)..... (3) o fl f2, f3s f4 sont des régions voisines d'intensité connue dans une rangée, les multiplicateurs fi, f2, f3 et f4 sont les kernels dans la direction X, r = isubX-1)/ISUB +1 A = r(r-1) B = A(r-2> isubX est le nombre de lignes de grille verticales, et 201 ISUB est le nombre total de lignes de grille entre les régions spécifiées d'intensité connue; b) obtenir des facteurs de pondération utilisant la fonction polynomiale P(isubY) afin d'obtenir des kernels dans la direction X conduisant à un kernel de convolution 4 x 4 pour chaque intersection de grille. 11. Procédé d'interpolation dans une matrice pour obtenir des valeurs interpolées & partir de valeurs connues, ledit procédé comprenant les opérations consistant à produire des kernels qui sont indépendants des valeurs connues ou la dimension de matrice, mais sont fonction de l'emplacement séparant les pixels de valeur connue avec des pixels nécessitant des valeurs interpolées, à effectuer une convolution de pixels voisins de valeur connue avec un kernel déterminé par ô;5 emplacement du pixel ayant besoin d'une valeur interpolée pour obtenir des valeurs pour les pixels nécessitant les valeurs interpolées. 12. Système de traitement d'une image normale en entrée pour fournir une image zoom en sortie (31), l'image en entrée (21) comprenant une matrice d'entrée (27) de rangée et de colonne de région ayant des valeurs d'intensité en entrée correspondant aux rangée et colonne de pixels dans une image d'affichage normal, ladite image zoom en sortie (31) comprenant une matrice zoom (28) de rangées et de colonnes de régions ayant des valeurs d'intensité correspondant aux rangées et colonnes de pixels dans l'image zoom affichée (31), ledit procédé comprenant: un moyen pour spécifier un facteur de zoom et des régions de valeurs d'intensité en entrée pour lesquelles un zoom a été réalise, un moyen pour transférer les valeurs d'intensité en entrée aux régions spécifiées pour lesquelles un zoom doit Etre réalise afin de désigner des régions de la matrice zoom (28), lesdites régions désignées étant séparées par des régions de valeurs d'intensité inconnue, un moyen pour sélectionner des kernels mémorisés individuels pour chacune desdites régions d'intensité inconnues de la matrice zoom (28), lesdits kernels mémorisés étant sélectionnés en fonction de l'emplacement de la région d'intensité inconnue dans la matrice zoom (28>, un moyen servant à sélectionner des régions voisines à partir d'emplacements déterminés par l'emplacement de 0C} la région d'intensité inconnue et le facteur de zoom (15), et un moyen servant à effectuer une convolution desdits kernels mémorises sélectionnés avec des valeurs d'intensité desdites régions voisines sélectionnées avec -_. des valeurs d'intensité connue comprenant les régions désignées de la matrice zoom (28) afin d'obtenir des valeurs d'intensité pour les régions de valeurs d'intensité inconnue. 135. Système selon la revendication 12. comprenant un S moyen (24) servant à sélectionner 16 régions voisines pour satisfaire una fonction polynomiale. 14. Système selon la revendication 12, dans lequel ledit moyen de sélection (24) comprend un moyen (26) servant à sélectionner 9 régions voisines pour satisfaire une fonction binomiale. 15. Procédé selon la revendication 12, dans lequel le moyen de sélection (24) dudit kernel mémorisé comprend: un moyen servant à recouvrir une grille sur la matrice zoom (28), ladite grille comportant une pluralité 1' de lignes horizontales entre lesdites régions d'intensité connue dans lesdites rangées et une pluralité de lignes verticales entre lesdites régions d'intensité connue dans lesdites colonnes, un kernel associé à chaque intersection desdites lignes horizontales et verticales -- étant présent, et un moyen servant à sélectionner ledit kernel le plus voisin du centre de ladite région d'intensité inconnue pour utilisation par le moyen pour effectuer une convolution afin d'obtenir les valeurs d'intensité pour lesdites régions d'intensité inconnue. 16. Système selon la revendication 15 comprenant une grille 16 x 16. 17. Système selon la revendication 12, dans lequel lesdits kernels sont indépendants de la donnée d'image 3o *_ i(15). 19. Système selon la revendication 12, dans lequel lesdits kernels sont basés sur un polynôme. -. Système selon la revendication 12, dans lequel lesdits kernels sont basés sur une fonction binomiale. 21. Système pour interpolation permettant d'obtenir une valeur interpolée à partir de valeurs connues dans une matrice de pixels, ledit système comprenant un moyen (24) servant à produire les kernels qui sont indépendants de la donnée matricielle ou du format de la matrice, un moyen servant à séparer les pixels de valeur connue à l'aide de pixels nécessitant une valeur inconnue, un moyen servant à effectuer une convolution de pixels de valeur connue avec lesdits kernels produits afin d'obtenir des valeurs pour des pixels nécessitant des valeurs interpolées.