Brevet FR2591005A1 — publié le 19870605

Procédé d'identification de structures arborescentes dans des images numériques et son application à un dispositif de traitement d'images. L'invention se rapporte au domaine de la reconnaissance de formes dans des images, et plus particulièrement à l'identification de structures arborescentes. Elle s'applique au traitement d'images numéri- ques, notamment dans le domaine médical, par exemple pour l'identifi- cation des vaisseaux sanguins sur des images obtenues par angiographie numérique. De nombreux problèmes de reconnaissance de forme peuvent se ramener à l'étude de structures arborescentes, la reconnaissance étant obtenue lorsqu'il est possible de faire correspondre une structure dite "modèle" à une structure "objet" détectée dans l'image à traiter. L'invention a pour objet un procédé d'identification de structures arborescentes dans des images numériques. Pour cela, les données détec- tées dans l'image d'une part et celles du modèle d'autre part doiventêtre représentées conventionnellement de la même manière, par un ensemble de données caractérisant une structure arborescente. Il existe une méthode dite "méthode de Selkow" décrite par exemple par l'article intitulé "Selkow: tree to tree editing problem" tiré de "Information Processing letters" volume 6, numéro 6, décembre 1977, qui permet de comparer deux structures arborescentes en évaluant une "distance" entre les deux structures. Le procédé selon l'invention utilise cette méthode pour trouver, une distance minimale entre deux arbres, puis détermine le chemin conduisant à cette distance minimale, pour aboutir à l'identification entre la structure arborescente objet et le modèle auquel on cherche à l'identifier. Selon l'invention un procédé d'identification de structures arbores- centes détectées dans des images numériques, par lequel les branches d'un arbre objet peuvent être identifiées à celles d'un arbre modèle, ces deux arbres étant définis sous forme d'une racine et d'une liste de sousarbres disjoints, formés chacun de la même manière ou éventuellement réduits à 2591005 s leur racine, et chacune des branches étant affectée d'un vecteur étiquette dont les composantes sont fonction des caractéristiques de ces branches, est caractérisé en ce qu'il consiste: dans une première phase, à transformer l'arbre objet en l'arbre modèle par une suite d'opérations élémentaires sur l'ensemble des sous- arbres, chaque opération étant affectée d'un coût fonction des composan- tes des vecteurs étiquettes des branches constituant les sous-arbres, ces opérations étant soit des changements d'étiquette, soit des insertions de sous-arbres de l'arbre modèle, soit des destructions de sous-arbres de l'arbre objet, et à sélectionner la suite d'opérations la moins coûteuse, dont le coût global mesure la distance entre les arbres selon la méthode itérative dite de "Selkow" connue en soi, - puis, dans une seconde phase, à reprendre la suite cd'opérations en sens inverse à partir de la dernière opération effectuée, de façon à réaliser pour chaque étape du calcul itératif, selon l'opération élémentaire qui a permis d'obtenir le coût minimal: - soit un changement d'étiquette: une sous-arbre de l'arbre objet étant identifié à un sous-arbre de l'arbre de référence; - soit l'insertion d'un sous-arbre de l'arbre modèle, dans ce cas il n'y a pas de sous-arbre correspondant dans l'arbre objet; - soit la destruction d'un sous-arbre de Parbre objet; dans ce cas il n'y a pas de sous-arbre correspondant dans l'arbre de référence; le résultat de l'identification étant formé par l'ensemble des changements d'étiquette des sous-arbres réduits à leur racine obtenus dans cette seconde phase et le coût associé qui correspond à la distance minimale des sous-arbres ayant ces branches pour racines. L'invention a également pour objet l'utilisation de ce procédé dans un dispositif de traitement d'images numériques, par exemple pour la reconnaissance des vaisseaux sanguins dans des images obtenues par angiographie numérique. L'invention sera mieux comprise et d'autres caractéristiques appa- raîtront à l'aide de la description qui suit en référence aux figures annexées. La figure 1 représente un exemple d'arbre modèle; La figure 2 est le schéma des relations entre les branches de l'arbre représenté sur la figure 1; 2 5 9 1005 La figure 3 représente un exemple d'arbre objet détecté dans une image numérique; La figure 4 représente le schéma des relations entre les branches de l'arbre objet représenté sur la figure 3. Les figures 5 à 12 sont des tableaux illustrant le calcul de la distance entre deux arbres. Pour le procédé d'identification de structures arborescentes dans une image numérique suivant l'invention, les données caractéristiques de ces structures dans l'image et dans le modèle sont représentées conven- tionnellement de la même manière, selon le schéma suivant: Une structure en forme d'arbre est définie par une racine, qui est par exemple la partie entre deux lignes quasi-parallèles ayant la plus grande largeur, et des branches également formées par des couples de lignes quasi-parallèles reliées à la racine elle-même ou à d'autres bran- ches par des arcs, chaque branche étant la racine de sous-arbre. Chaque branche d'un arbre est caractérisée par un vecteur dit vecteur étiquette dont les composantes sont des paramètres caractéristi- ques; par exemple dans des images obtenues par angiographie numérique, chaque branche représente un vaisseau sanguin pour laquelle le vecteur étiquette peut comporter des composantes telles que la largeur, la longueur, la forme, du vaisseau considéré auquel on associe également le nom de vaisseau. Les données extraites de l'image numérique sont structurées de la même manière par détection puis paramètrisation des contours dans l'image, détection des branches par reconnaissance de lignes quasi- parallèles dans l'image, puis formation du schéma des relations entre les branches, chaque branche étant représentée par un vecteur étiquette dont les composantes ont la même structure que celles de l'arbre modèle. La figure I représente un exemple cParbre modèle T: chaque branche bi a un vecteur étiquette noté Ei, i =1 à 7; b est la racine, b2, b3, b6 et b7 sont des branches reliées à la racine bl, b4 est reliée à b3 et b5 est reliée à b6. L'étiquette de la racine de T, b1 dans l'exemple ci- desus est notée E(T). La figure 2 représente les relations entre ces branches telles qu'elles apparaissent sur la figure 1. Chaque branche est un noeud dans ce graphe des relations. La figure 3 représente un exemple d'arbre objet, T' détecté dans une image numérique, chaque branche bj ayant un vecteur étiquette Ej, j = a, b,..., h. Selon la notation ci-dessus E(T') = Ea. La figure 4 représente le schéma des relations entre ces branches telles qu'elles apparaissent sur la figure 3: bb, bc, bd et bg sont reliées à la racine ba, bfest reliée à be elle même reliée à bd, et bh est reliée à bg. Chaque branche est un noeud dans ce graphe des relations. Pour la description de cet arbre modèle T et de l'arbre objet T', chaque arbre est défini comme comportant une suite de sous-arbres disjoints, un ordre étant introduit entre ces sous-arbres: par exemple les sous-arbres sont décrits dans le sens trigonométrique, à partir de la racine à laquelle on donne un sens. Le modèle T peut alors être représenté par une liste ordonnée L(T) de sous-arbres: L(T) = (bl (b2Xb3(b4")) (b5(b6)) (b7)), cette liste précise les liaisons entre la racine T 0= b et les sous-arbres T1 = b2, T2 = b3 (b4) T3 = b5 (b6), et T4 = b7. Ti est le sous-arbres de T d'ordre i. De la même manière, larbre des données est représenté par la liste L(T') suivante: L(T') = (ba(bb)) (b) (bde(e(bf))) (bg(bh"), qui précise les liaisons entre la racine T0' = ba et les sous-arbres suivants T1' = bb, T2' = bc, T3' = bd(be(bf)) et T4' = bg(bh). Un sous-arbre tel que T2 par exemple comporte lui-même une racine b3 et un sous-arbre b4. Des sousarbres tels que b4 ou T1 = b2, sont réduits à leurs racines. Comme indiqué ci-dessus, pour calculer une "distance" entre arbres qui permettra ensuite de rechercher la distance "minimale" abou- tissant à ridentification de l'arbre objet à l'arbre modèle, il est nécessaire de définir cette "distance". Un arbre étant défini par la succession des sous-arbres qui le composent, ordonnée comme indiqué dans les listes ci- dessus, il est possible de doter l'ensemble des sous-arbres d'un certain nombre cd'opérations qui permettront ensuite de définir la distance entre deux arbres. On définit d'abord l'égalité entre deux arbres: deux arbres sont dits égaux si leurs racines ont la même étiquette, et si les listes ordonnées de leurs sous-arbres sont identiques. Trois opérations élémentaires sont alors définies dans l'ensemble des arbres par des opérations sur les étiquettes: 1. Changement d'étiquette d'un arbre: Un arbre T peut être transformé en un arbre T par un change- ment d'étiquette de la manière suivante: l'étiquette de sa racine E(T) est remplacée par l'étiquette de la racine de T, E(Td); la liste des sous- arbres de T est la même que la liste des sous-arbres de T. 2. L'insertion d'un sous-arbre dans la liste descriptive d'un arbre: un sous-arbre peut être inséré dans une liste ordonnée de sous-arbres descriptive d'un arbre T, à un rang donné, pour former un arbre modifié T^. L'étiquette de sa racine ne change pas: E(T") = E(T) et la liste initiale: L(T) = TO, Til... Ti, Ti+...Tm des sous-arbres de T devient la liste L(Tx) après insertion d'un sous-arbre A: L(Tt) = TO, TI,...Ti, A, Ti+ Tm. 3. Destruction d'un sous-arbre T. dans la liste ordonnée des sous- arbres constituant un arbre T pour former un arbre T^: L'étiquette de la racine ne change pas: E(TE) = E(T), et la liste initiale L(T) = TO Ti... Tiil Ti, T +i,...Tm devient après destruction de Ti: L(T") = T0, Tl... Ti_1, Ti+i...Tm 1+1 m A chacune de ces opérations élémentaires est associé un coût qui dépend des valeurs affectées à la racine ou aux différentes branches des sous-arbres insérés ou détruits. La valeur affectée à une branche peut être une combinaison de composantes du vecteur étiquette de cette branche. L'exemple le plus simple consiste à prendre l'une des compo- santes du vecteur étiquette soit par exemple la largeur des branches. Pour l'opération de changement d'étiquette, il est nécessaire de disposer d'une matrice des coûts de transformation d'une étiquette Ei en une étiquette Ej soit CT (Ei, E.) pour tous les i et j. Par exemple, ce coût peut être la différence entre les largeurs pour des sous-arbres réduits à des branches. Pour les autres opérations, destruction et insertion de sous-arbres dans un arbre, il est nécessaire de disposer des coûts élémentaires de destruction et d'insertion d'un arbre réduit à une branche étiquetée Ei: soient CD (Ei) et Ci (Ei) respectivement ces coûts élémentaires de des- truction et d'insertion. Le coût d'insertion Ci(A), d'un sous-arbre A dans une liste ordonnée de sous-arbres est égal à la somme des coûts d'insertion des sous-arbres constitutifs de A. Pour un sous-arbre réduit à une branche ce coût d'insertion est seulement le coût élémentaire indiqué ci-dessus, soit la largeur d'une branche par exemple. De la même manière, le coût de destruction d'un sous-arbre A est égal à la somme des coûts de destruction des sous-arbres constitutifs du sous- arbre A. Pour un sous-arbre réduit à une branche, ce coût peut être, de la même manière que pour le coût d'insertion, une valeur affectée à cette branche, soit la largeur par exemple. La mesure de la distance entre deux arbres A et B est alors obtenue au moyen d'un procédé récursif, tel que décrit par l'algorithme de Selkow indiqué ci-dessus qui affecte un coût global aux suites de transfor- mations possibles pour transformer l'arbre A en l'autre B: ces suites de transformations se font par une suite d'opérations élémentaires affectées chacune d'un coût élémentaire. La distance minimale d(A,B) entre les deux arbres étiquetés A et B sera le coût global de la façon la moins coûteuse de transformer l'arbre de départ en larbre cd'arrivée. Ce coût prend en compte à la fois les modifications de structures, insertions et suppressions de sous-arbres, et les modifications d'étiquettes. Soient deux arbres, A de sous-arbres, Ai,...Am, et B de sous- arbres, B1...Bn; A0 et B0 sont les racines, et T(i) est Parbre obtenu en enlevant d'un arbre T les sous-arbres de l'ordre i + 1 au dernier de la liste: si Tm est le dernier sous-arbre de la liste descriptive de T, on a en particulier T(m) = T et T(0) = T0' racine de T. Selkow a montré que la distance d(A, B) entre un premier arbre A et un second arbre B est telle que: _ d(A(0), B(j)) = CT (E(A), (E(B)) + i C (Bk) k=l I i - d(A(i), B(0)) = CT (E(A), E(B)) + CD(Ak) k=l - d(A(i), B(j)) = Min.d(A(i-1), B(j-I)) + d(Ai, Bj) 30.d(A(i), B(j-l)) + Cl (Bj) d(A(i-l")), B(j) + CD (Ai). L'étude détaillée montre en effet que les opérations, changement d'étiquette, insertion, ou destruction d'un sous-arbre, ne changent jamais l'ordre des sous-arbres. D'o, pour la mesure des distances un calcul de type programmation dynamique. Le calcul s'en déduit alors immédiatement: 1. Faire pour j = 0 àn: j d(A(0), B(j)) = CT (E(A), EB)) + T: CI (Bk) k=l 2. Faire pour i = 0 à m d(A(i), B(O")) = CT (E(A), E(B)) ± CD(Ak) k=l 3. Faire pour i = 1 à n Faire pour j = 1 à m d(A(i), B(j) = Min (d(A(i-l), B(j-l) + d(Ai, Bj); d(A(i), B(j-I) + C1 (Bj); d(A(i-l), B(j) + CD (Ai) Ce calcul est récursif puisqu'il est nécessaire de calculer d(Ai, Bj) pour calculer d(A(i), B(j>), Ai et Bj étant eux-mêmes des arbres. Les figures 5 et 6 représentent respectivement un autre exemple d'arbre modèle et d'arbre objet auquel on doit l'identifier. L'exemple de calcul qui suit appliqué aux arbres représentés par les schémas des figures 5 et 6 illustre le calcul de distance selon la méthode, connue en elle-même, décrite ci-dessus. La valeur d'une branche est notée entre parenthèses dans les noeuds représentés sur les figures 5 et 6 sous l'étiquette de la branche notée seulement par son indice a à g dans l'arbre modèle et 1 à 6 dans l'arbre des données. On définit les coûts suivants: Coût de changement d'étiquette CT (Ei, Ej) = cCi- Cj CoOt d'insertion = coût de destruction = Ck. La solution est représentée dans une suite de tableaux sur les figures 7 à 12. Le tableau TAB0 permet de calculer les distances entre les sousarbres des premières listes. L'arbre modèle est représenté par la liste L(T) = (a(b(cXd))(e)(f(g))) avec les sous-arbres suivants: -a racine;T0; b(cXd) sous-arbre T1, (c) et (d) sous-arbres de T1; - (e) sous-arbre T2; pas de sous-arbre; - f(g) sous-arbres T3; (g) sous-arbre de T3; De rmême rarbre des données est représenté par la liste L(T') = (1(2(3X4X)(5))(6)); - 1: racine rT'0; - 2(3X4X5) sous-arbre ' 1' (3), (4) et (5) étant des sous-arbres; - (6) sous-arbre T'2; pas de sous-arbre. - La figure 7 représente le premier tableau illustrant la solution: TAB0. Les entrées de ce tableau sont en ligne, la liste des sous-arbres de T et en colonne la liste des sous-arbres de T'. Pour chaque sous-arbre a été calculée une valeur égale à la somme des valeurs des branches et indiquée sous le sous-arbre correspondant; ces valeurs sont cerclées. Sur la première ligne apparaissent les distances calculées entre la racine lo et les arbres T(0), T(1), T(2), T(3) = T, ces distances sont calculées dans la première phase de calcul: d(T'(0); T(j)). De même, dans la première colonne apparaissent les distances calculées entre la racine TO et les arbres T'(0), T'(1), T'(2) = T', ces distances sont calculées dans la seconde phase de calcul d(T'i), T(0). Les autres calculs sont faits dans la troisième phase en utilisant éventuelle- ment les calculs intermédiaires nécessaires tels qu'il ressort des tableaux suivants. Le tableau TAB1 représenté sur la figure 8 permet de calculer la distance minimale entre T1 et Tr, soit 15; le tableau TAB 2 représenté sur la figure 9 permet de calculer la distance entre T2 et T'1 soit 42; le tableau TAB3 représenté sur la figure 10 illustre le calcul de la distance entre T1 et T'2, soit 41; le tableau TAB4 représenté sur la figure 11 illustre le calcul de la distance minimale entre T3 et T'l,soit 19; et le tableau TAB5 de la figure 12 illustre le calcul de la distance entre T3 et T'2, soit 13; la distance entre T2 et '2 est obtenue immédiatement du fait que les 2 sous-arbres sont réduits à leurs racines par la différence des valeurs (20 + 10) soit 10. Ces valeurs soulignées dans TAB0 sont alors utilisées pour le calcul final de la distance minimale entre les arbres T' et T tel qu'illustré par le tableau TAB0 de la figure 7. La distance minimale ainsi obtenue est 48 dans lexemple illustré par les figures 5 à 12. Le problème à résoudre alors est d'identifier les noeuds objet et modèle à partir de ces résultats, c'est-à-dire remonter le chemin qui a conduit à l'obtention de la distance minimale entre l'arbre objet et Parbre modèle. Pour cela chaque coût obtenu à une phase de calcul est affecté dans la case correspondante d'une flèche indiquant les choix effectués à chaque étape lors des tests sur les coûts minimaux, ou la provenance de la 2 5 9 1 0 0 5 valeur prise comme valeur de départ pour le calcul de coût lorsqu'il n'y a pas eu de choix. La méthode pour retrouver le chemin conduisant à la distance minimale consiste à remonter le chemin en suivant les flèches reliant les coûts, ce qui donne dans l'exemple ci-dessus les associations suivantes obtenues à la lecture du tableau TAB0 T'2 = T3; coût:13 (coût complexe) insérer T2; coût:10 (coût d'insertion) T'1 = T1; coût: 15 (coût complexe) 10.T'0 = To; coût 10 (par changement d'étiquette). Soit un coût total, égal à 48. L'identification n'est pas terminée car les changements d'étiquet- te T'1 = T1 et T'2 = T3 sont des opérations complexes qui doivent être décomposées. Les décompositions résultent de la lecture du tableau TAB1 pour l'identification T'1 = T1 et de la lecture du tableau 5 pour l'identifi- cation T'2 T3. Par le tableau TAB1 de la figure 8 on obtient: (5) = (d) changement d'étiquette: coût 1 (4) détruit: coût de destruction: 3 (3) = (c) changement d'étiquette: coût I (2) = (b) changement crd'étiquette: coût complexe 15. Par le tableau TAB5 de la figure 12 on obtient insérer (g) (non retrouvé dans l'arbre objet): coût 3 (6) =(f)changement d'étiquette coût complexe 13. L'identification est alors complète puisque réalisée jusqu'au niveau des branches. D'une manière plus générale, pour obtenir l'identification entre un arbre objet et un arbre modèle, représentés par des listes ayant comme premier élément la racine et comme éléments suivants les sous- arbres, si les indices i = I à n correspondent aux sous-arbres de l'arbre modèle et j-1 à n aux sous-arbres de rl'arbre objet, comme indiqué ci- dessus, le chemin ayant permis d'aboutir à la distance minimale est remonté à partir de (m, n). La description qui suit décrit l'application de ce procédé à un dispositif de traitement d'images, notamment pour des images obtenues par angiographie numérique. L'angiographie est une technique radiologique de visualisation des vaisseaux sanguins dont le principe est le suivant: un produit de contraste est injecté dans un vaisseau et on prend, à intervalles de temps réguliers, des radiographies de la région du corps à examiner. La séquence d'images obtenues par soustraction avec une image de référence, prise avant injection du produit de contraste, fournit l'image des vaisseaux et donne des indications sur la vitesse de propagation du produit de contraste dans ces vaisseaux. Cet examen est utilisé pour détecter des pathologies, principale- ment les sténoses (obstruction) ou les anévrismes (dilatation) de ces vaisseaux. L'examen se fait actuellement de manière entièrement manuelle. Dans le but d'automatiser un tel examen, une première étape consiste à reconnaître de manière automatique les vaisseaux par rapport à un modèle anatomique du réseau sanguin dans la région d'examen considérée. Les différentes phases du traitement utilisé sont les suivantes: 1. Choix dans la séquence d'images différence d'une image présen- tant un bon contraste. 2. Détection des contours, filtrage, pontage et paramétrisation de ces contours 3. Détection des vaisseaux. 4. Formation de rarbre des relations entre vaisseaux. 5. Mise en correspondance de l'arbre vasculaire détecté et du modèle anatomique mis sous la même forme. La méthode décrite précédemment est utilisée pour réaliser cette dernière phase de mise-en correspondance. Les phases préalables I à 3, nécessaires à la formation dans les images numériques des structures arborescentes sont en dehors du cadre de la présente invention. Pour Papplication du procédé d'identification selon l'invention les éléments suivants sont considérés pour la formation de la structure modèle: le réseau artériel a naturellement une structure d'arbre en trois dimensions. Si Pon prend comme racine la plus grosse artère, celle-ci se divise ensuite en vaisseaux plus petits et ainsi de suite. L'image angiographique corres- pond à la projection de cette structure sur un plan. Pour une même région, les images obtenues peuvent être très différentes suivant Pl'orientation de ce plan de visualisation. On considérera donc un modèle différent pour chaque type d'angiographie (région d'intérêt et angle de visualisation). Ces modèles sont fournis par les praticiens. Le modèle sous forme d'arbre pour lapplication du procédé est - constitué de la manière suivante: La plus grosse artère correspond à la racine de l'arbre. Les vaisseaux qui en partent sont les branches, reliées à la racine par un arc, et ils sont eux-mêmes racines des sous-arbres formés par leurs vaisseaux affluents. Un ordre est introduit sur ces sous-arbres. Par exemple, on parcourt les sous-arbres dans le sens trigonométrique à partir de la racine orientée dans le sens de la circulation du sang. Comme indiqué ci-dessus chaque branche de l'arbre est étiquetée par un ensemble de paramètres caractéristiques.du vaisseau qu'il représente. Cette étiquette est en fait un vecteur constitué d'un certain nombre de composantes telles que largeur, longueur, forme etc... avec également le nom du vaisseau considéré. Les données extraites de l'image angiographique doivent être naturellement structurées de la même manière. C'est l'étape 4 indiquée ci-dessus appelée "formation de Parbre des relations entre vaisseaux" au paragraphe précédent. Il faut déterminer les liaisons entre vaisseaux et le sens de ces liaisons et lever les différentes indéterminations (croisement ou jonction, sens). On arrive alors a une description de rl'image à étiqueter sous forme d'arbre o chaque branche est également représentée par un vecteur de composantes de même structure que celui caractérisant les branches du modèle. Pour utiliser la méthode de mise en correspondance décrite précédemment, il faut définir les différents coûts élémentaires: coût de changement détiquette, coût d'insertion et coût de destruction. Ces coûts sont calculés à partir du vecteur étiquette associé à chaque branche. A titre d'exemple, comme indiqué ci-dessus, les coûts d'insertion ou de destruction sont égaux à la largeur de vaisseau, et le coût de substitution à la valeur absolue de la différence des largeurs des deux vaisseaux considérés. Pour la simulation, une image test du type de celle représentée sur la figure 3 a été choisie, faisant apparaître un grand nombre de vaisseaux et donnant donc un arbre assez complet. De plus, des segments parasites et des trous dans les contours ont été introduits pour se rapprocher d'un cas réel. Le modèle choisi pour effectuer la mise en correspondance correspond à l'arbre représenté sur la figure 2. Une représentation possible est donnée figure 1. Pour que le test soit assez critique, le modèle choisi était assez différent des données: la branche b7 du modèle n'apparaît pas dans les données; les branches bc et bf des données n'apparaissent pas dans le modèle). Le résultat de 'association obtenue est: ba = b; bb= b2, b3= bd, b4 = be, b5 = bh, b6 = bg. Tous les vaisseaux de l'image des données existant sur le modèle ont bien été identifiés. Les autres ont été détectés comme devant être détruits ou insérés dans le processus de mise en correspondance. De même, l'application du procédé à des données réelles a permis de montrer la validité du procédé. L'invention n'est pas limitée au procédé tel que décrit dans ses détails, ni à son application à des images numériques angiographiques. REVENDICATIONS 1. Procédé didentification de structures arborescentes détectées dans des images numériques, par lequel les branches d'un arbre objet peuvent être identifiées à celles d'un arbre modèle, ces deux arbres étant définis sous forme d'une racine et d'une liste de sous-arbres disjoints, formés chacun de la même manière ou éventuellement réduits à une branche, et chacune des branches étant affectée d'un vecteur étiquette dont les composantes sont fonction des caractéristiques de ces branches, caractérisé en ce qu'il consiste: - dans une première phase, à transformer l'arbre objet en l'arbre modèle par une suite d'opérations élémentaires sur l'ensemble des sous- arbres, chaque opération étant affectée d'un coût fonction des compo- santes des vecteurs étiquettes des branches constituant les sous-arbres, ces opérations étant soit des changements d'étiquette, soit des insertions de sous-arbres de l'arbre modèle, soit des destructions de sous-arbres de l'arbre objet, et à sélectionner la suite d'opérations la moins coûteuse, dont le coût global mesure la distance entre les arbres selon la méthode itérative dite de "Selkow" connue en soi, - puis, dans une seconde phase, à reprendre la suite d'opérations en sens inverse à partir de la dernière opération effectuée, de façon à réaliser pour chaque étape du calcul itératif, selon l'opération élémentaire qui a permis d'obtenir le coût minimal: - soit un changement d'étiquette: un sous-arbre de l'arbre objet étant identifié à un sous-arbre de l'arbre de référence; - soit l'insertion d'un sous-arbre de l'arbre modèle, dans ce cas il n'y a pas de sous-arbre correspondant dans l'arbre objet; - soit la destruction d'un sous-arbre de l'arbre objet; dans ce cas il n'y a pas de sous-arbre correspondant dans l'arbre de référence; le résultat de l'identification étant formé par l'ensemble des changements d'étiquette des sous-arbres réduits à leur racine obtenus dans cette seconde phase et le coût associé qui correspond à la distance minimale des sous-arbres ayant ces branches pour racines. 2. Application du procédé selon la revendication 1 à un dispositif de traitement d'images numériques. 3. Application selon la revendication 2 à des images obtenues par angiographie numérique pour l'identification des vaisseaux sanguins.