Brevet FR2620544A1 — publié le 19890317

PR1CED D'INTEPLATION La présente invention à pour objet un procédé d'interpolation de la valeur d'un point dans une maille à N sommets. L'interpolation étudiée est de préférenoe de type bilinéaire encore que l'invention s'applique à des cas o elle ne serait ni linéaire ni même double mais plut8t multiple. Dans l'exemple décrit la maille comporte quatre sommnets mais dans certaines applications, en particulier dans des applications d'imagerie en représentation d'objets à trois dimensions, la maille pourrait comporter un nombre différent de sommets: par exemple trois, six ou huit sommets. Le procédé d'interpolation de l'invention trouve une application dans le domaine ndical co il peut être utilisé pour corriger des distorsions d'images acquises avec les gamma-caméras. Il peut aussi être utilisé dans le domaine vidéo o il permet de représenter des images de faible définition, par exemple à 64 X 64 éléments d'image, comme des images haute définition, par exemple à 256 X 256 éléments d'image. Le principal intérêt du procédé d'interpolation de l'invention réside dans la rapidité avec laquelle l'interpolation peut être faite. Cette rapidité autorise le traitement en temps réel en respectant pour visualiser les images les normes habituelles de représentation vidéo télévision. Les gamma-caméras sont utilisées en médecine nucléaire pour visualiser dans un organe la répartition de molécules marquées par un isotope radioactif que l'on a injecté à un patient. Une gamma-caméra camprend généralement un collimateur pour focaliser les photons gammas émis par le patient, un cristal scintillateur pour transformer les photons gamma en photons lumineux ou scintillations, et un réseau de tubes photomultiplicateurs qui transforment - 2 - chacune des scintillations en des impulsions électriques dites contributions électriques de tubes. Elles comportent en outre des circuits électroniques pour produire à partir des cortributions électriques fournies par les tubes photo- multiplicateurs, des signaux de coordonnées X et Y du lieu o s'est produite une scintillation, ainsi qu'un signal de validation Z quand l'énergie W de la scintillation appartient à une bande d'énergie préd&éterminée. Cette chaîne de détection est généralement suivie d'un ensemble de visualisation qui peut comnporter un oscilloscope cathodique commandé par les signaux de coordonnés X, Y, et par Z pour visualiser par un point lumineux sur l'écran le point d'impact du photon gamma sur le cristal. Cet impact est aussi appelé évènement d'image. L'ensemble de visualisation peut éventuellement conporter un dispositif Iphotographique pour former une image de l'organe observé en intégrant un grand nombre de points lumineux produits sur l'écran cathodique. Il peut par ailleurs comprendre un traitement numérique des images. En particulier l'ensemble de visualisation peut être adapté à la représentation de tomographies de l'organe observé. Pour atteindre oe but on acquiert plusieurs images de oet organe selon une pluralité d'orientations d'observation de la gamma- caméra par rapport a cet organe. Par des traitements du signal, analogues à oeux rencontrés dans les tomodensitometres, on peut reconstituer des images de coupes des organes examinés. En autres qualités une gamma-caméra doit posséder une bonne résolution spatiale, c'est à dire la capacité de distinguer des petites sources radioactives rapprochées, une bonne réponse en taux de comptage, c'est à dire la capacité de traiter un grand nombre d'événements par unité de temps, -3- et une qualité d'image le plus possible indépedante de l'énergie de l'isotope considéré. La résolution spatiale dépend de la précision des calculs des oordck-.nes X et Y de chacun des évènerents d'image. La qualité de l'élaboration de ces coordonnées dépend de la précision des mesures et des lois physiques régissant le fonctionnemnt des différentes parties de la gammna-caméra. Ainsi l'interaction d'un photon gamma avec le cristal donne naissance à une scintillation lumineuse dont l'intensité décroît expotentiellement avec le temps. Cette scintillation est vue par plusieurs tubes photo- multiplicateurs simultan&nent. Les photons lumineux caposant cette scintillation arrachent des photo-électrons aux photo-cathodes des tubes photo-multiplicateurs. Le nanbre des photo-électrons arrachés obéit, pour une scintillation donnée, à la loi statistique de Poisson. Par ailleurs, à énergie constante, la contribution électrique est une fonction sensiblement gaussienne de la distance séparant le centre du tube photomultiplicateur de la projection du lieu o s'est produite la scintillation. Si la scintillation se produit à l'aplanb du centre du tube, la - contribution électrique est maximum. Plus le lieu de la scintillation est éloigné du centre du tube plus la contribution électrique est faible. A titre d'exnemple, si une scintillation se produit à l'aplomb de la paroi du tube la contribution électrique de celui-ci est environ réduite de moitié par rapport à la contribution électrique maximum. Une scintillation est vue par plusieurs tubes pihoto-multiplicateurs simultanément, généralement 6 à 10 tubes. Aussi, la détermination de l'emplacement de cette scintillation sur le cristal, elle-même représentative du lieu d'émission du photon gammna d'excitation (et donc de l'évènement d'image), peut être obtenue en calculant 4- l'emplacement du barycentre des contributions électriques délivrés par 1' ensemble des tubes phto-multiplicateurs excites par cette scintillation. Ce calcul s'effectue simplement, selon le brevet américain ANG N 3011057, en injectant les contributions électriques au travers d'un jeu de matrices de résistanoes dont les valeurs des résistances sont fonction des positions des tubes photo-nultiplicateurs auxquels elles sont raccordées. Les positions de oes tubes sont repérées par rapport à des axes cartésiens de référence dont le point d'intersection est généralement situé au centre du réseau de tubes. Dans chaque matrioe il y a autant de résistances qu'il y a de tubes photo-multiplicateurs dans le réseau de tubes. Chacune des résistances est reliée d'une part à la sortie d'un tube photo-multiplicateur différent et d'autre part à un point commun qui constitue la sortie de la matrioe. Ces résistances effectuent- ainsi une pondération des contributions électriques de chacun des tubes de photamultiplicateurs qui les alimentent. Un des problèmes présentés par les gamna-caméra du type ANG est qu'elles présentent des distorsions géométriques de l'image liées à la structure de captation de la lumière: cristal scintillateur-tube photanultiplicateur- matrice de barycentration. Les progrès de la médecine nucléaire, et en particulier les améliorations faites pour recueillir davantage d'informations, et des informations meilleures, à partir des gamma-camnéras, par exemple pour la détection de petites tumeurs, onduisent à optimiser la résolution spatiale au détriment de la linéarité caractéristique inhérente à la coneption et à la réalisation des caméras. Il en résulte une distorsion spatiale génante des images. Ces distorsions sont à l'origine de défauts d'uniformité de densité qui peuvent être importants. A titre d'exemple, une --5 -- contraction de 0,4 mn du rayon de l' image d'une surface circulaire d'un centimétre de rayon conduit à un défaut d'uniformité de 8 %. Ces défauts d'uniforrr té peuvent être particulièrement gênants en tomographie o l'effet d'amplification des défauts d'uniformité dûs aux procédés de reconstruction peut être supérieur à un facteur 10. Dans l'état de la technique on a tenté de corriger les effets de oette distorsion. Le principe général de la correction est le suivant. On effectue la mesure d'une image d'une mire régulière placée entre une source radioactive d'némission uniforme et la gamna-caméra. Si la gamracaméra était parfaite, l'image représentée devrait correspondre à la répartition régulière des trous ou des fentes de la mire. La mesure effectuée montre en fait des irrégularités provenant des effets de distorsion spatiales. On peut cependant utiliser dans 1' image ainsi acquise la connaissance de la distorsion, mesurée en conparant l'image acquise à la répartition théorique à laquelle on aurait dû aboutir, pour effectuer les corrections des images ultérieurement acquises avec la mêne gamma-caméra. Dans un état de la technique, constitué par le brevet européen N 0021. 366, un calculateur élabore une matrice de correction de 64 X 64 nouvelles coordonnées U et V d'une image non distordue correspondant à un même narmbre de coordonnées X et Y de l'image distordue de la mire.En fonctionnement de correction, les coordonnées acquises X et Y des évneents d'images sont codées sur 12 bits. Les 6 bits de poids fort de X et Y adressent la mémoire des U et V. Puis une interpolation linéaire est effectuée en utilisant les 6 bits de poids faible de X et Y. On utilise ainsi le décalage d'un évènent d'image à placer par rapport aux quatre angles dé la maille élémentaire U, V, de l'image corrigée dans laquelle l'élément d'image doit finalement se trouver. Ccmnpte tenu d'un champ image d'environ 400 mm de -6- diamètre, une codification sur 12 bits des coordonnées des évènements d'image devrait normalement conduire à une définition d'image de l'ordre de 0,1 mm. On a pu montrer que cette définition ne pouvait être obtenue réellement que si la matrice de correction était en fait une matrice de 128 X 128 mailles et si donc on utilisait 7 bits de poids fort pour déterminer dans quelle maille de l'image corrigée doit se trouver l'évènement d'image à placer. Cette détermination est équivalente à attribuer des corrections de coordonnées aux sommets des 128 X 128 mailles de l'image acquise. Ceci étant, cette différence n'influe pas dans le principe de la présente invention. En effet la représentation des images avec 128 X 128 évènements d'image est aussi trop grossière. Il est nécessaire d'atteindre une résolution plus fine. Compte tenu du fait que les coordonnées acquises des évènements d'image (de l'image distordue) sont données sur 12 bits, il reste 5 bits pour faire une interpolation de la position vraie d'un évènement d'image dans la maille corrigée à laquelle on vient de déterminer qu'il doit appartenir. On calcule par interpolation le déplacement à appliquer à 1'évènement, à partir des corrections de coordonnées à appliquer aux sommets des mailles. Si on appelle x et y les coordonnées, relatives aux bits de poids faible, de l'élément d'image dans la maille acquise, l'interpolation consiste à déterminer DX DY les déplacements à appliquer à la position de l'évènement situé à l'intérieur d'une maille afin de la corriger. Cette correction dépend de x et de y et de la manière dont la maille acquise a été transformée pour devenir la maille corrigée. Si on appelle A,B, C, et D, les samnets d'une maille carrée dans l'image acquise, ces sammets avec la correction de distorsion doivent respectivement être transformés c'est à dire subir des corrections de coordonnées DXa-DYa, DXb-DYb, DXc-DYc et DXd- -7 - DYd selon les deux axes de référence de l'image. Ces corrections de coordonnées sont stockées dans une mémoire dite de correction. La correction de distorsion de la position de l'évrnlent d'image à replacer est alors donnée par les formules suivantes: DX = (1-x-y+xy) DXa + (x-xy) DXb + (y-xy) DXc + (xy) DXd DY = (1-x-y+xy) DYa + (x-xy) DYb + (y-xy) DYc + (xy) DYd En fait la méthode classique de calcul d'interpolation passe, pour le calcul de chaque déplacement DX ou DY, par deux niveaux de calculs internmédiaires. On calcule d'abord une première correction CI à partir des coordonnées x (ou respectivement y) et des corrections de coordonnées des deux sommets de la maille DXa (ou DYa) et DXb (ou DYb). Avec deux accès à la mànoire o sont enregistrées les coordonnées des sam:ets des mailles corrigées on obtient ainsi un résultat intermédiaire CIx (ou CIy). On calcule ensuite une autre correction CJ par interpolation entre les deux autres sar.ets C et D de la maille. On obtient de la même façon CJx (ou CJy) avec.deux autres accès à la mémnoire. Le calcul final consiste, à partir de CI et de CJ à obtenir le déplacement DX (ou DY). On a donc, pour obtenir chacun des déplacmeents d'image d'un l'évènement d'image dans l'image corrigée, procédé à quatre acos mémoire (en A, B, C et D), à quatre multiplications en tenant compte des coordonnées x (ou y), à deux additions pour déterminer CI et CJ, à deux accès aux valeurs intermédiaires CI et CJ, à deux autres multiplications en tenant compte des coordonnées y cette fois (ou x) et enfin à une addition pour obtenir le déplaoement final DX (ou DY). Ces mrnbreuses opérations nuisent à la rapidité de la présentation des images corrigées. -8- Lorsque le nombre de bits de poids faible à interpoler est peu important, par exemple trois bits, il est possible de précalculer pour chaque maille corrigée les corrections définitives qu'il faut attribuer à des évènements d'image placés à l'intérieur des mailles acquises. On procède ensuite à un adressage direct pour obtenir les valeurs préinterpolées. Comme la détermination de la maille corrigée a elle-même été obtenue par adressage, oette technique revient à effectuer un adressage direct de la totalité de la correction. Cette technique présente oependant un inconvénient: elle est inexploitable si le nombre de bits d'adressage est trop important. Ceci signifie qu'elle est inexploitable si la précision de la résolution attendue de l'image est trop grande ou, ce qui revient au mnme, si l'interpolation doit être pratiquée sur un nombre de bits plus grand. En effet lorsque les coordonnées sont données sur 12 bits tout procédé de correction d'image basé en définitive sur une telle technique d'adressage direct, nécessite une capacité mémoire exorbitante. En particulier pour une image à 4096 X 4096 points d'image, il faudrait une mamoire de 16 Méga mots de 24 bits (12 bits pour déterminer X corrigé, et 12 bits pour déterminer Y corrigé). Cette mémoire serait adressée par les X et Y acquis et codés chacun sur 12 bits. Cette capacité dépasse les limites technologiques actuellement accessibles. L'invention a pour objet de remédier à ces inconvénients en proposant un procédé d'interpolation ou la capacité mémoire n'est pas particulièrement augmentée mais o par contre le calcul d'interpolation peut être fait en temps réel. Le fait que le calcul d'interpolation puisse être fait en temps réel peut être utilisé dans le domaine vidéo. En effet la représentation des muvmeents nécessite l'acquisition et le stockage de nombreuses images. De manière à ne pas trop augmenter les capacités mémoire au -9- momnent de oe stockage, il peut être utile, étant donné qu'on n'est intéressé qu'au mouvement, de ne stocker que des images à définition dégradée. Dans certains cas, en particulier dans l'étude des mouvements cardiaques par gamma- caméra, les images obtenues ont dès le dcpart une faible définition. Cependant, au moment de la visualisation, il peut être nécessaire d'amortir les différences de niveau de luninosité et d'augmenter la résolution de l'image en recréant, a partir des images stockées dans la mémoire en faible nombre (faible résolution) des images à haute résolution (introduction d'évènements d'image intermédiaires dans chaque maille d'événements d'image stockée) et avec un nombre de niveaux de gris enrichi. On ne peut cependant créer des évènements d'images intermédiaires que dans la mesure o on peut leur donner, dans la même ganmme de niveaux de gris, des niveaux de gris intermédiaires. L'amortissement des transitions lumineuses ne se produit ainsi que si l'écart de niveau de gris entre éléments d'images voisins dans l'image faiblement définie est important. Le calcul en temps réel autorisé par l'invention permet alors la création "au vol" de telles images haute définition à partir d'images basse définition stockées en grand nrombre. L'invention a donc pour objet un procédé d'interpolation multiple de la valeur d'un point, dans une maille à N sotnets, caractérisé en ce qu'il comporte les étapes suivantes: - on détermine une finesse d'interpolation à obtenir, - on précalcule des contributions d'interpolation relatives à tous les lieux possibles de la maille, pour la finesse déterminée, et en fonction des coordonnées de ces lieux dans cette maille, - et on attribue au point conme valeur interpolée le résultat d'une composition algébrique des contributions, relatives aux coordonnées du point dans la mille, et de - 10-- coefficients pondérateurs des sommets de la maille. L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen des figures qui l'acccmpagnent. Celles-ci ne sont données qu'à titre indicatif et nullement limitatif de l'invention. Les figures montrent - figure 1: une représentation schématique fonctionnelle du pocédé d'interpolation selon l'invention; - figure 2: un schéma fonctionnel des moyens mis en oeuvre dans le procédé de l'invention; - figure 3: une variante préférée du schéma précédent dans une application de type vidéo. La figure 1 représente schlûatiquement le procédé d'interpolation selon 1' invention. Un point P, ou événement d'image, est défini, après une acquisition, par ses coordonnées X et Y. Dans un exemple préféré, l'acquisition est une acquisition obtenue à la suite d'une expérimentation avec une gamma-camera et, pour fixer les idées, les coordon- nées X et Y sont données sous forme binaire avec 12 bits de définition. Le point P appartient à une image distordue ID qui est inexploitable. Pour pouvoir être utilisée oette image doit être corrigée par un procédé de correction de distorsion d'image. En définitive, il faut calculer les coordonnées d'un point P' dans une image corrigée IC. Dans les procédés de correction de distorsion on attribue à chaque maille élémentaire M de 1' image distordue ID acquise une maille corrigée M' dont les sommets A, B, C et D sont définis par des corrections de coordonnées, dits aussi ici coefficients pondérateurs, DXa - DYa, DXb... DYd. En déterminant la résolution du maillage de correction il est possible d'effectuer une correction de distorsion en deux temps. Un premier temps consiste à attribuer à une maille M une maille M'. Le second temps consiste à interpoler la position du point P' dans la maille M' en fonction de la _-11 - position du point P dans la maille M et en fonction des coefficients pondérateurs des sommrets de la nmaille. La fléche 1 représente la première opération. Pour la mener à bien on utilise les bits de poids fort des coordonnées X et Y. Dans une réalisation particulière on utilise les sept bits de poids fort de manière à définir un maillage de 128 X 128 mailles dans les images distordues ID et corrigées IC. Dans l'invention on détermine la finesse d'interpolation à obtenir. Dans l'exemple présenté,la finesse d'interpolation à obtenir concerne les cinq bits de poids faible, notés x et y, des coordonnées X et Y. Si on appelle Ka la contribution (1-x-y+xy), Kb la contribution (x- xy), Kc la contribution (y-xy), et Kd la contribution (xy), l'invention a pour but de précalculer les contributions de Ka, Kb, Kc, et Kd, pour toutes les valeurs possibles de x et y correspondant à la finesse d'interpolation à obtenir. On stocke ensuite ces valeurs dans des mémoires mortes appelées matrices d'interpolation. Dans l'exemple cette finesse est mesurée sur cinq bits: chaque coordonnée x ou y peut donc prendre 32 valeurs. En adressant des mémoires préprogramikes contenant les différentes contributions Ka, Kb, Kc et Kd par les coordonnées x et y, on peut obtenir directement les valeur des fonctions représentatives de ces contributions. On détermine ultérieurement DX et DY en appliquant les formules suivantes: DX = Ka.DXa + Kb.DXb + Kc.DXc + Kd.DXd DY = Ka.DYa + Kb.DYb + Kc.DYb + Kd.DYd Puisque x et y sont ici codés sur cinq bits chacun, chaque mémoire morte Ka, Kb, Kc et Kd comnporte 1024 cases mémoire. Dans la pratique chaque case mémoire comporte une contribution codées sur huit bits. Il y a quatre matrices d'interpolation. Autrement dit l'interpolation nécessite une -12 - augmentation mémoire de 4 K octets pour stocker l'ensemble des contributions relatives. Ccrmpte tenu du fait que le maillage de correction comporte huit zones mémoires dites de correction de coordonnées et relatives aux déplacements des sommets des mailles: DXa, DYa, DXb,... DYd, du fait que ces corrections de coordonnées sont codées chacune sur 11 bits (un bit de signe, huit bits significatifs et deux bits fractionnaires), et du fait que le maillage de correction est relatif à une résolution 128 X 128, la taille mémoire globale est de 2 fois 16 K mots de 11 bits de mémoire vive, et de 4 K octets de mémoire morte. En effet les corrections de coordonnées doivent pouvoir changer, en fonction du vieillissement de la gamma-caméra par exemple, alors que les contributions relatives précalculées sont immuables. Elles peuvent donc être stockées dans une mémoire morte. La taille mémoire supplémentaire pour les matrices d'interpolation est par conséquent négligeable. Dans l'application il est utile de remarquer que le précalcul des contributions d'interpolation peut être simplifié en utilisant, pour la programmation des mémoires Ka, Kb et Kc les résultats obtenus pour le précalcul de la mémoire Kd. La multiplication x.y n'a pas être refaite. La figure 2 représente le mode préféré de détermination des coordonnées corrigées des point P' dans la maille M'. Un bus 2 achemine les coordonnées X et Y jusqu'aux mémoires de corrections de coordonnées 3 à 6 d'une part, et jusqu'aux mémoires d'interpolation 7 à 10 d'autre part. Un multiplicateur-accumulateur 15 relié aux sorties de ces memoires reçoit ensuite des mémoires de corrections, les corrections de coordonnées DX (ou DY) d'une part, et des mémoires d'interpolation 7 à 10 les contributions relatives d'interpolation K d'autre part. Il effectue les multiplications correspondantes et 1' accumulation des - 13 - résultats de oes multiplications. IL délivre selon le cas le déplacaement DX puis le déplacement DY. Avec le multiplicateur accumulateur 15 le traitement est effectué en série. Un circuit de correction 150 déli.-;e ensuite les coordonnées corrigées des événements d'image en additionnant aux coordonnées distordues les déplacements calculés. Ajouté au fait que la taille mémoire a ainsi peu augmenté le procédé de l'invention peut aussi permettre de plus de gagner du temps. En effet l'adressage des corrections de coordonnées des différents saomnets de la maille de correction M' peut être simultané. En effet X et Y, ou du moins les 7 bits de poids fort de X et Y, représentent le sacnet A de la maille. Les sanets B, C et D s'en déduisent par une addition ou une soustraction d'une unité de bit de poids le plus faible parmi les bits de poids fort. En conséquence on peut avoir accès simultanément aux 8 coefficients pondérateurs DXa, DYa,etc. Dans le mnême temps il est également possible d'accéder simultanément aux contributions Ka, Kb, Kc et Kd. En remplaçant alors le multiplicateur-accumulateur 15 par un jeu de quatre multiplicateurs affectés respectivement aux mémoires 3 et 7, 4 et 8, 5 et 9 et 6 et 10, et en les reliant à un additionneur on peut accélérer les calculs des déplacements. En effet lors d'un cycle d'horloge suivant celui ou les mémoires ont été adressées il est possible au moyen des multiplicateurs d'effectuer chacune des multiplications K.DX (ou K.DY). Au troisinème temps de cycle on dispose des valeurs de déplacement DX (ou DY). Trois temps de cycles encore plus tard on dispose de l'autre valeur de déplacement DY (ou DX). Il est même possible de disposer de ces deux valeurs de déplacement simultanément en doublant encore le nanombre des multiplicateurs. L'application vidéo représentée en partie sur la figure 3, est semblable. Dans celle-ci une trame d'une image -14 - disponible a une faible résolution. Par exemple elle comporte 64 lignes et 64 colonnes: il y a 4 K positions possibles des évènements d'image et dans un exemple chaque élément d'images peut être codé sur 8 bits soit 256 niveaux de gris,de luminanoe ou de chraninance. Le stockage des images est donc peu exigeant. La reconstruction, au vol, d'une image plus finement définie, par exemple avec un pas de 256 X 256 rend nécessaire le calcul du niveau de gris en 16 lieux d'une maille du maillage 64 X 64. Autrement dit une maille repérée par quatre saimets adjacents appartenant à deux lignes adjacentes de l'image stockée doit être morcelée en 4 x 4 = 16 parties. La coordonnée x codée sur 2 bits doit pouvoir valoir de 1 à 4 et la coordonnée y aussi codée sur 2 bits doit pouvoir prendre également les valeurs de 1 à 4. Connaissant le niveau de gris, le coefficient pordérateur, Ca, Cb, Cc et Cd en chacun des sarmets A, B, C et D de la maille, il est possible, par un calcul similaire à celui présenté plus haut de calculer les niveaux degris intermédiaires des points intermédiaires contenus dans la maille. Dans un perfectionnement de l'invention on s'est oepedant rendu compte que, plutôt que de n'effectuer que le calcul des contributions Ka, Kb,Kc et Kd, il était ici préférable et possible d'effectuer le précalcul de toutes les multiplications nécessaires. En effet, étant donné que les niveaux de gris ne sont donnés pour chaque point que sur huit bits (256 niveaux), il est possible de remplaoer les mémoires d'interpolation Ka, Kb, Kc et Kd à 16 cases mémoires (16 positions possibles dans la maille) par des mémoires d'interpolation comnposites à 16x256 cases mémoires. Les précalculs effectués dans ces mémoires composites sont effectués sur 8 bits. - 15 - Les transitions brutales de niveau de gris entre des éléments d'image voisins de 1' image faiblement définie sont alors adoucies par l'introduction d'éléments d'image intermédiaires, dont les niveaux de gris so:nt fixés à des niveaux intermédiaires dans la Mmne gamme de niveaux de gris, à ceux de ces éléments d'image voisins. Dans une première phase de ce traitement rapide des images sources de faible définition sont lues à chaque temps de cycle pour déterminer les niveaux de gris, ou coefficients pondérateurs, Ca à Cd des quatre saumets A à D d'une maille. La mémoire source peut être quadruplée de manière à être adressée simultanément en X et Y pour délivrer en nMe temps les coefficients pondérateurs Ca à Cd. De manière à éviter le coût imposé par le quadruplement de la mnmoire on peut l'organiser différemment en séparant les évènements d'image pairs en x et y. Cette organisation différente permet alors de lire simultanent, mais en des zones différentes de cette ménoire, les coefficients pondérateurs Ca à Cd des quatres sommnets de la maille. Au besoin un multiplexage de lecture de cette mémoire peut aussi être organisé. La mémoire de l'image source, ou pour simplifier l'explication les m&noires sources 16 à 19, délivrent alors les quatre coefficents pondérateurs sur 8 bits représentatifs chacun du niveau de gris d'un des samets de la maille. Ces niveaux de gris peuvent être utilisés comme entrée d'adresse dechacune des métoires 20 à 23 de précalcul composite ayant effectué le produit K.C,et recevant comme autre adresse des codes correspondant successivement à chacune des quatre positions intermédiaires des points interpolés 41 à 44 situés sur une même ligne dans la maille. Dans une ligne y est fixé, seul x varie (de 1à 4). Les mémoires conposites 20 à 23 délivrent alors chacune successivement quatre résultats corresporndant au calcul de -16 - la luminosité pondérée en quatre lieux dans la maille en fonction d'un des sanets de oette maille. Pour chaque lieu un jeu d'additionneurs 24 à 26 permet d'additionner les résultats délivrés par les quatre mémoires 20 à 23 et de délivrer la luminosité pondérée corresporndant à ce point. Avec un temps de cycle de lecture de 166 rnaosecondes, les quatres luminosités peuvent être interpolées en 664 rxnoiosecondes. Les luminosités en 256 points d'une ligne vidéo peuvent donc être ainsi calculées en 43 microsecondes. La visualisation de l'image peut alors être obtenue en mode vidéo classique. En effet, en lecture vidéo classique à 25 images par seconde,la lecture d'une ligne d'une image dure 64 microsecondes. La omnpatibilité de ces deux durées (64 microsecondes étant supérieur à 43 microsecondes) qualifie le procédé d'interpolation ainsi présenté pour la visualisation "au vol" d'images faiblement définies en des images hautement définies. Le précalcul des mémoires coapozsites présente de plus une particularité. Si on accepte de coter x= O et y = O le centre d'une maille, on se rend compte que la contribution composite attribuée à un point peut être simplifiée. Il suffit de ne calculer qu'une seule mémoire camposite et de la réaliser en quatre exemplaires. Seules les valeurs d'adressage utilisées pour accéder à ces mémoires canposites changent. Selon la contribution canposite envisagée, selon le sawnet, on fera subir aux adresses x et y des inversions de valeurs -x et -y. Cette manière d'agir, présente un intérêt. Lors de la fabrication il n'y a pas lieu de se préoccuper d'une différentiation quelconque entre les mnoires 20 à 23. Seul des décodeurs respectifs d'accès 28 à 31 sont ainsi particularisés en fonction de leur cbjet. Le mode d'interpolation préconisé dans ces deux - 17- applications est de préférence un mode linéaire. Canpte tenu du fait qu'on cherche à représenter des images à deux dimensions le mode d'interpolation sera nm..e bilinéaire. Cependant l'invention peut être transposée à une interpolation d'un point dans un espace gécnétrique à trois dimensions, ou même dans un espaoe théorique à M dimensions. En outre l'interpolation n'est pas rnécessairement du type linéaire. D'autres fonctions que la fonction linéaire sont aussi envisageables. S'agissant d'une maille à N saoets, il peut seulement être nécessaire que la saure des contributions relatives à chacun de oes scomnets soit égale à 1. En outre la définition d'une maille à N sacnets ne doit pas être comprise uniquement dans le cas o N vaut quatre one ce qui vient d'être décrit. Dans le cas o la maille à une forme différente, par exemple triangulaire ou hexagonale,le nombre de saommnets peut être différent. Dans ces derniers cas, plutôt que de calculer les contributions selon une définition cartésienne, il peut s'avérer préférable d'utiliser une définition polaire. - 18 - REVENDICATIONS 1 - Procédé d'interpolation multiple de la valeur d'un point dans une maille à N somnets, caractérisé en ce qu'il Homporte les étapes suivantes: - on détermine une finesse d'interpolation à cbtenir, - on prcalcule des contributions (Ka, Kb, Kc, Kd) d'interpolation relatives à tous les lieux de la mille (A,B, C,D), pour la finesse déterminée et en fonction des coordonnées de ces lieux dans cette maille, - on attribue au point comme valeur interpolée le résultat d'une composition algébrique des contributions relatives aux coordonnées (x,y) du point, et de coefficients pdérateurs des sommets de la maille. 2 - Procédé selon la revendication 1 caractérisé en oe que l'interpolation multiple est une interpolation double. 3 - Procédé selon la revendication 2 caractérisé en oe que l'interpolation double est une interpolation bilineaire. 4 - Procédé selon l'une quelconque des revendications i à 3 caractérisé en ce que N vaut quatre. - Procédé selon l'une quelconque des revendications 1 à 4 caractérisé en ce que la valeur représente une absisse (DX) d'un point à interpoler mesurée selon un axe. 6 - Procédé selon l'une quelconque des revendications 1 à 4 caract6érisé en ce que la valeur représente la luminance ou la chraminance d'un élément d'image. 7 - Procédé selon la revendlication 6 caractérisé en ce qu'on précalcule (20,23) en partie les valeurs interpolées en fonction de valeurs prédéterminées des coefficients pondérateurs. - 8 - Procédé selon la revendication 7 caractérisé en ce que le précalcul des contributions d'interpolation porte - - 19 - un précalcul omposite d'un jeu de contributions et un multiplexage (x, -x, y, -y) de l'adressage de oes contributions.