Brevet FR2844899A1 — publié le 20040326

PROCEDE ET DISPOSITIF DE MODELISATION ET DE RECONSTRUCTION DES VOIES AERIENNES D'UN PATIENT La présente invention se rapporte au domaine de la modélisation et de la reconstruction des voies aériennes de patients. La présente invention se rapporte plus particulièrement à un procédé et un dispositif de 10 modélisation et de reconstruction des voies aériennes d'un patient à partir de données tomodensitométriques acquises en mode spiralé. Dans un contexte clinique général d'investigation des voies respiratoires, le médecin dispose aujourd'hui de trois grandes approches: la bronchographie, les explorations fonctionnelles respiratoires et la tomodensitométrie. La bronchographie est la plus ancienne des 20 techniques de diagnostic des maladies bronchiques. Elle consiste en l'injection d'un produit de contraste à l'intérieur des voies aériennes étudiées et en l'acquisition de clichés radiologiques sous une incidence appropriée. En dépit de sa capacité à visualiser et à 25 analyser la morphologie globale de la partie étudiée de l'arbre bronchique, cette technique a été pratiquement abandonnée en raison de son caractère hautement invasif. Les explorations fonctionnelles respiratoires sont les méthodes traditionnelles employées dans le 30 diagnostic des maladies respiratoires. Elles reposent sur une estimation globale des modifications de calibre bronchique suite à l'administration de différents stimuli, en mesurant des paramètres tels que le débit d'air et la pression. A partir de ces mesures, des indices globaux sur 35 les calibres bronchiques peuvent être obtenus. Toutefois, une analyse morphologique globale et une analyse locale des voies aériennes n'est pas possible avec ces techniques. La tomodensitométrie (TDM) offre une représentation bidimensionnelle précise des structures anatomiques, sous forme d'images en coupe, reconstruites selon une orientation choisie par le médecin (généralement dans le plan axial). L'acquisition des coupes s'effectue soit en mode conventionnel (acquisition coupe par coupe) 5 soit en mode spiralé (acquisition continue suivie par une reconstruction des coupes aux endroits désignés). L'analyse des données TDM fournit les meilleures indications sur les calibres bronchiques et leur variation sous l'action de divers stimuli. C'est la raison principale pour laquelle la 10 TDM est devenue l'outil de référence dans le diagnostic des maladies pulmonaires. Toutefois, en routine clinique, plusieurs problèmes sont à noter: a Le radiologue ne dispose pas d'un aperçu 15 global de la morphologie 3D des bronches. Pour cela il doit construire mentalement une représentation 3D à partir des sections 2D disponibles sur les coupes axiales. Ce processus demande de l'expérience et un long entraînement préalable. * La visualisation de certaines configurations (distorsions,...) est parfois très difficile sur les seules coupes axiales. * Une analyse complète de la morphologie des bronches nécessite l'exploration d'un nombre élevé de 25 coupes axiales (pouvant aller jusqu'à 300 voire 400 images par patient). a La communication entre le radiologue et le chirurgien amené à pratiquer une intervention sur les voies respiratoires souffre de l'absence d'une modalité 30 unificatrice de repérage spatial 3D. * L'évaluation des calibres bronchiques reste fastidieuse pour une application en routine clinique et elle est limitée aux seules sections bronchiques dont l'axe est orthogonal au plan de coupe. 35 Pour résoudre ces différents problèmes, il est nécessaire de disposer d'une reconstruction 3D performante de l'arbre bronchique, ainsi que d'un mécanisme adéquat de visualisation, navigation et interaction avec les données 3D. Les techniques actuelles de reconstruction 3D reposent principalement sur deux approches: - l'approche 2D/3D consistant en une segmentation. coupe par coupe des sections bronchiques suivie par une restauration de l'ensemble 3D, et - l'approche 3D utilisant des données de volume et reposant en général sur des méthodes de croissance de 10 régions à partir d'un noyau sélectionné de manière automatique ou semi automatique. L'inconvénient majeur de ces techniques réside dans leur incapacité à permettre de reconstruire le réseau bronchique au-delà du 4ème ordre d'embranchement. Une technique de reconstruction 3D de l'arbre bronchique a été déjà développée et publiée dans [C. Fetita, F. Prêteux, "Bronchial tree modeling and 3D reconstruction", Proceedings SPIE Conference on Mathematical Modeling, Estimation and Imaging, San Diego, 20 CA, Vol. 4121, August 2000] et [C. Fetita, "Analyse morphofonctionnelle des voies aériennes en TDM spiralée volumique", Thèse de Doctorat, Université Paris V - René Descartes, France, Décembre 20001. Cette technique permet d'obtenir une reconstruction de l'arbre bronchique jusqu'au 25 6ème ordre d'embranchement, mais elle se limite au seul type de données spécifié dans lesdits ouvrages, et échoue sur des données fournies par autre type de scanner. La présente invention entend remédier aux 30 inconvénients de l'art antérieur en proposant un procédé plus général et plus robuste permettant, à partir de données TDM acquises en mode spiralé avec un protocole spécifique, de produire de manière automatique une représentation 3D interactive de l'arbre bronchique 35 jusqu'au 6 e voire jusqu'au 7èm ordre d'embranchement, permettant ainsi une analyse morphologique globale et locale à partir de l'axe central indexé de l'arbre. Pour ce faire, le procédé selon l'invention est du type décrit ci-dessus et il est remarquable, dans son acception la plus large, en ce qu'il comprend successivement: * une étape de mise en forme desdites données afin d'obtenir un volume pulmonaire en niveaux de gris Vg, * une étape d'extraction de l'arbre bronchique à partir dudit volume pulmonaire en niveaux de gris Vg par croissance morphologique conditionnelle à une 10 topologie colorée, afin d'obtenir un arbre binaire en trois dimensions, et * une étape finale de représentation de l'arbre bronchique, afin de permettre de naviguer dans lesdites voies aériennes. L'étape de mise en forme des données comporte de préférence un fenêtrage permettant la visualisation des tissus parenchyme pulmonaire, ainsi qu'éventuellement, une étape de filtrage. L'étape d'extraction de l'arbre bronchique par 20 croissance morphologique conditionnelle à une topologie colorée comporte de préférence une étape de filtrage morphologique en trois dimensions du volume pulmonaire en niveaux de gris Vg, puis une étape de reconstruction en trois dimensions par une technique de croissance 25 conditionnelle à une topologie colorée. Le filtrage est de préférence opéré de manière itérative selon la fonction suivante 0255, xe X Fn(x) =,OnSkVg(X),XESUPPVg-(YUX) LSk Vg (x), sinon pour l'étape d'initialisation et F,,'(x) = max(Fn' - (x), min F,'(y) ),Vxe supp Vg o X est l'ensemble des parois horizontales du volume VI, Y est l'ensemble des parois verticales du volume Vg, Sk un opérateur d'homogénéisation de la paroi bronchique 35 et n est un opérateur extensif tel que (OVg - Vg) (x) F N(x)= FN+I(x),Vx Esupp Vg. Sk est donné de préférence par: Vf:93 -_> 9, Vx e supp f, Skf(x) = min max {f(y) Iy-x X| o di représente une des 13 directions principales de l'espace 3D, |.- la norme L2 et k la taille 15 du filtre prenant valeurs dans l'intervalle [1, 10]. O, est de préférence donné par: v Of = min(f B, f D Bj), o B1 est un élément structurant en deux dimensions défini dans le plan axial par: Bn =t (iO)}t)u (1O,)u tu. (iu)} [On [O,n] [,n] [O,n] v et o B. est son symétrique par rapport à l'origine et Bl'addition Minkowski. Le filtrage morphologique 3D s'effectue de préférence de la manière suivante: 1)- Détection automatique de la section axiale de la trachée et de son diamètre DT; 2)- Filtrage itératif de Vq avec le filtre Fp, o p noté par A'. La deuxième étape de l'extraction de l'arbre bronchique A. implique une reconstruction 3D par une technique de croissance en topologie colorée, à partir du sous-ensemble A,. De manière générale, l'ensemble A, peut croître en occupant un point libre x de son voisinage, si le changement d'état du point x remplit une condition de croissance. La condition de croissance est spécifiée par une combinaison de fonctions logiques booléennes validant ou non l'appartenance au réseau bronchique d'un site quelconque x, exprimée par 2 0 Vkiv vL(X)çI i, j > 1, o v et A représentent respectivement la somme et le produit logique, et L,, une fonction logique booléenne définie à partir de l'information de densité Vg relative à une topologie donnée. L'implémentation du schéma de reconstruction peut envisager plusieurs modules, chacun définissant sa propre fonction booléenne (cf. figure 5). Ces modules peuvent être permutés, rajoutés ou enlevés (à condition de 30 conserver au minimum un d'entre eux), avec des répercussions plus ou moins importantes sur le résultat de la reconstruction. Les fonctions logiques booléennes L1j sont de 35 préférence données par: L11, (x) = L(x, x), avec L(x, p) = 3r c 9, (p) u 02 (p), tel que Vg(x) direction locale de croissance et F un chemin fermé. L,2(x)= min Vg(y) à partir duquel s'effectue la croissance locale. y(Vg(y)-Vg(x)) X (Vg(y)-Vg(x))2 L22(X)= -l y_ AnZ926(X) yEA-026(6(x) YE 26(X) E 1{,Al} 1{)yA} Y. 926 (x) YE 926 (x) (Vg(y)- Vg(X)) > (y -A Vg26(x) -A +{ Vg(x) Vg(x)> k} >YE 926(X) y 1 le,4 YE Z26 (x) o 826(x) désigne le voisinage 26-connexe de x, A csuppVg l'ensemble courant des points de l'arbre 15 reconstruit, et k une valeur de niveau de gris dans l'intervalle [50, 120]. Dans le procédé selon l'invention, on opère en 20 outre de préférence une indexation de l'arbre bronchique, ainsi qu'une extraction de l'axe central. Avantageusement, le procédé selon l'invention permet de modéliser et de reconstruire les voies aériennes 25 d'un patient à partir de données tomodensitométriques acquises en mode spiralé, quel que soit le scanner utilisé. Avantageusement également, plusieurs types de représentations sont possibles: * représentation en mode bronchographie TDM 30 virtuelle, * représentation surfacique ou en fil de fer, ou encore * représentation indexée (avec délimitation des segments anatomiques). L'ensemble de ces représentations forme ce qu'on appelle un arbre bronchique interactif qui offre la possibilité d'une navigation complexe endo- et exoluminale, une analyse tant globale que locale avec des 10 mesures de calibre bronchique. En outre, à partir d'informations morphométriques stockées dans l'axe central indexé, la synthèse de modèles de voies aériennes adaptés au patient est également possible. La présente invention se rapporte également à un dispositif de modélisation et de reconstruction des voies aériennes permettant la mise en oeuvre du procédé selon l'invention, ledit dispositif comportant au moins: * un module de mise en forme desdites données afin d'obtenir un volume pulmonaire en niveaux de gris Vg, * un module d'extraction de l'arbre bronchique à partir dudit volume pulmonaire en niveaux de gris Vg par croissance morphologique conditionnelle à une 25 topologie colorée, afin d'obtenir un arbre binaire en trois dimensions, et * un module de représentation de l'arbre bronchique, afin de permettre de naviguer dans lesdites voies aériennes. On comprendra mieux l'invention à l'aide de la description, faite ciaprès à titre purement explicatif, d'un mode de réalisation de l'invention, en référence aux figures annexées: * la figure 1 illustre un schéma général de l'invention * la figure 2 illustre une étape de mise en forme des données; * la figure 3 illustre une étape d'extraction de l'arbre bronchique; * la figure 4 illustre un filtrage morphologique 3D; a la figure 5 illustre une reconstruction 3D par la technique de croissance; et * la figure 6 illustre un exemple de représentation et d'investigation 3D interactive grâce au procédé selon l'invention. Le procédé selon l'invention, illustré figure 1, est un procédé de modélisation et de reconstruction des voies aériennes d'un patient à partir de données 15 tomodensitométriques acquises en mode spiralé. Le procédé selon l'invention repose sur un système logiciel qui reçoit en entrée des données tomodensitométriques (TDM) volumiques acquises en mode 20 spiralé avec une collimation fine (de 1 à 2 mm) et un pas d'hélice de 1 à 2, et fournit à la sortie une représentation 3D interactive de l'arbre bronchique proximal (jusqu'au 7ème ordre d'embranchement). Le système comporte trois modules principaux, 25 comme illustré figure 1, qui seront détaillés par la suite: * un module de lecture et mise en forme des données pour la mise en oeuvre de l'étape de lecture et de mise en forme des données; * un module d'extraction de l'arbre bronchique 3D pour la mise en oeuvre de l'étape d'extraction de l'arbre bronchique 3D; et * un module de représentation et investigation 3D interactive pour la mise en oeuvre de 35 l'étape de représentation et investigation 3D interactive. Le module de lecture et de mise en forme des données récupère les données volumiques à la sortie du scanner (se présentant soit sous forme d'une série de coupes 2D, soit sous forme d'un volume 3D) et les transforme en un volume en niveaux de gris, noté Vg, incluant le champ pulmonaire complet sur toute la longueur 5 de l'acquisition, comme illustré figure 2. Cette conversion inclut une étape de filtrage (pour les scanners fonctionnant en mode de reconstruction rétroprojection filtrée) et un fenêtrage permettant la visualisation des tissus de la parenchyme pulmonaire. Le filtre, réglable 10 directement sur la console du scanner, est choisi parmi ceux proposés par le constructeur de la console d'acquisition, dans une plage basse-moyenne fréquence. La fenêtre de visualisation présente, par exemple des paramètres compris entre: - -600 à -200 Unités Hounsfield (UH) pour le centre de la fenêtre, et - 1200 à 2000 UH pour la largeur de la fenêtre. Le module d'extraction de l'arbre bronchique 3D part du volume pulmonaire V, pour permettre d'isoler le réseau des lumières bronchiques (les régions de circulation de l'air à l'intérieur de l'arbre bronchique), noté Aw sur la figure 3. Le module est composé de deux sous-modules qui 25 effectuent successivement un filtrage morphologique 3D et une reconstruction 3D par une technique de croissance conditionnelle à une topologie colorée. Le filtrage morphologique 3D a pour objectif 30 l'isolement d'un sousensemble maximal des lumières bronchiques capable de reconstituer par croissance 3D l'ensemble complet des voies aériennes. Le filtrage s'effectue de manière itérative, en utilisant un opérateur définit ciaprès: Soit YcVg c9V l'ensemble des parois verticales (les surfaces de frontière en plan coronal et sagittal) du volume Vg. Soit XCVg c9{3 l'ensemble des parois horizontales (les surfaces de frontière en plan axial) du volume Vg. Soit On: 9 -- 9 un opérateur extensif, Onf 2 f. 5 f:3 _9, tel que (OnVg Vg) (x) appartenant à une zone spatiale de texture quasi-homogène et d'étendue spatiale supérieure à n dans le plan axial. Ici h est un seuil de faible valeur, comprise entre O et 25. Un choix possible de 0n (non-unique) est donné par: v O.f =min(f B.,f B.), o B1 est un élément structurant 2D (non-) causal défini dans le plan axial par: B= {U (i,O)}U {EU(O, i)}U tu, (i, i)} v V B, son symétrique par rapport à l'origine et 15 L1'addition Minkowski. Soit Sk un opérateur morphologique donné de préférence par: Vf:93 -_9, VxE supp f, Skf(x)= min max {f(y) Il -x xd du filtre prenant valeurs dans l'intervalle [1,10]. Le filtre itératif Fn est défini de manière récursive comme suit: 255, xe X F. (x) = OnSkVg(x), XESUPPVg-(YUX) [SkVg (x), sinon pour l'étape d'initialisation et 30 F' (x) = max(F,''(x), min F'. (y)), Vx e supp Vg, F. yt(x)maF.- o i désigne l'itération courante et &(x)le voisinage causal 3D de x par rapport à la direction de balayage. Le nombre d'itérations nécessaires N est donné par la condition d'idempotence: F"N(x)=F"N+l(x),VxEsuppVg. Le filtrage morphologique 3D s'effectue de la manière suivante (cf. figure 4): 1) Détection automatique de la section axiale de la trachée et de son diamètre DT; 2) Filtrage itératif de Vg avec le filtre FP, o p " DT. Généralement, p=1 ou p=2. Le résultat est stocké dans V1; 3) Filtrage itératif de V, ou de V, avec le filtre Fq, o q=D,. Le résultat est stocké dans V2; 4) Seuillage adaptatif de la différence V2-V1 en fonction de l'histogramme des valeurs de 0 à 50; ) Dilatation morphologique 3D sous contrainte de positivité du gradient local. Le résultat du filtrage est un sous-ensemble de la lumière bronchique, noté par A'. 20 Le deuxième sous-module de l'extraction de l'arbre bronchique implique une reconstruction 3D par une technique de croissance en topologie colorée, à partir du sous-ensemble A,. De manière générale, l'ensemble A. peut croître en occupant un point libre x de son voisinage, si le changement d'état du point x remplit une condition de croissance. La condition de croissance est spécifiée par une combinaison de fonctions logiques booléennes validant ou non l'appartenance au réseau bronchique d'un site quelconque x, exprimée par 3 5 4;lv{L,(x)}} i, j > 1, o v et A représentent respectivement la somme et le produit logique, et L,, une fonction logique booléenne définie à partir de l'information de densité Vg relative à une topologie donnée. L'implémentation du schéma de reconstruction 5 peut envisager plusieurs modules, chacun définissant sa propre fonction booléenne (cf. figure 5). Ces modules peuvent être permutés, rajoutés ou enlevés (à condition de conserver au minimum un d'entre eux), avec des répercussions plus ou moins importantes sur le résultat de 10 la reconstruction. Les fonctions logiques booléennes Lij sont de préférence données par: L, (x) = L(x,x) avec L(x, p) = 3F c i(p) u O(p), telqueVg(x) à partir duquel s'effectue la croissance locale. ZE(Vg (y) - Vg (x)) (Vg (y) - Vg (X))2 r yL22(X)= Y lQA} -,4A26(x) jEAr926(x) L22Yx) 26= lY1{A}A yE 926 (X) yE 926 (X) (yVg (y) - Vg (x)) > YlIA} YZ26(x)-A +'{Vg (x) Vg (x) > k} ye926 (X) x) E 926 (X) o Z926(x) désigne le voisinage 26-connexe de x, A c: supp V9 l'ensemble courant des points de l'arbre reconstruit, et k une valeur de niveau de gris dans l'intervalle [50, 1201. Le module de représentation et investigation interactive offre une représentation et une manipulation interactive de l'arbre bronchique reconstruit. Les outils fournis et les différents types de sorties sont illustrés figure 6. L'arbre 3D binaire A. et le volume pulmonaire en niveaux de gris V1 constituent les entrées de ce module. Le procédé selon l'invention permet ainsi de réaliser plusieurs types de représentations: * représentation en mode bronchographie TDM virtuelle à travers une technique de rendu de volume appliquée à l'ensemble d'intersection entre la dilatation 3D unitaire de l'arbre 3D binaire A et du volume pulmonaire en niveaux de gris Vg, * représentation surfacique ou en fil de fer à travers une technique classique de rendu de surface, * représentation indexée (avec délimitation des segments anatomiques), soit sous forme d'axe central soit sous forme d'arbre 3D indexé. 25 L'invention est décrite dans ce qui précède à titre d'exemple. Il est entendu que l'homme du métier est à même de réaliser différentes variantes de l'invention sans 30 pour autant sortir du cadre du brevet. REVENDICATIONS 1. Procédé de modélisation et de reconstruction des voies aériennes d'un patient à partir de 5 données tomodensitométriques acquises en mode spiralé, caractérisé en ce qu'il comprend successivement: * une étape de mise en forme desdites données afin d'obtenir un volume pulmonaire en niveaux de gris Vg, * une étape d'extraction de l'arbre 10 bronchique à partir dudit volume pulmonaire en niveaux de gris V., par croissance morphologique conditionnelle à une topologie colorée, afin d'obtenir un arbre binaire en trois dimensions, et a une étape finale de représentation de 15 l'arbre bronchique, afin de permettre de naviguer dans lesdites voies aériennes. 2. Procédé de modélisation et de reconstruction des voies aériennes selon la revendication 1, caractérisé en ce que ladite étape de mise en forme des 20 données comporte un fenêtrage, ainsi qu'éventuellement, une étape de filtrage. 3. Procédé de modélisation et de reconstruction des voies aériennes selon la revendication 1 ou la revendication 2, caractérisé en ce que ladite étape 25 d'extraction de l'arbre bronchique par croissance morphologique conditionnelle à une topologie colorée comporte une étape de filtrage morphologique en trois dimensions du volume pulmonaire en niveaux de gris Vg, puis une étape de reconstruction en trois dimensions par une 30 technique de croissance conditionnelle à une topologie colorée. 4. Procédé de modélisation et de reconstruction des voies aériennes selon la revendication 3, caractérisé en ce que ledit filtrage est opéré de 35 manière itérative selon la fonction suivante 0255, xeX F0 (x) = OnSkVg(X), xE SUPp Vg -(YUX) LSk Vg (x), sinon pour l'étape d'initialisation et Fni (x) = max(F. (x), min F. (y)),Vxe suppVg yet(x) o X est l'ensemble des parois horizontales du volume V, Y est l'ensemble des parois verticales du volume 5 Vg, Sk est un opérateur d'homogénéisation de la paroi bronchique et On, est un opérateur extensif tel que (O.f f) (x) nécessaires N étant donné par la condition d'idempotence: 15 FN(x) = FN+l (x),Vx EsuppVg. 5. Procédé de modélisation et de reconstruction des voies aériennes selon la revendication 4, caractérisé en ce que Sk est donné par: Vf:3 _> 9, Vxe supp f, Skf(x) = m[Oin,] max {f(y) Iy - x1 du filtre prenant valeurs dans l'intervalle [1,10]. 6. Procédé de modélisation et de reconstruction des voies aériennes selon la revendication 4, caractérisé en ce que 0n est donné par: v O, f = min(f À B,, f G B,), o B, est un élément structurant en deux dimensions défini dans le plan axial par: B,, = U (iO)} t{U, (0,i)}U U, (i,i)} v et o B. est son symétrique par rapport à l'origine et sl'addition Minkowski. 7. Procédé de modélisation et de reconstruction des voies aériennes selon l'une quelconque des revendications 4 à 6, caractérisé en ce que le filtrage morphologique 3D s'effectue de la manière suivante: 1)- Détection automatique de la section axiale de la trachée et de son diamètre DT; 2)Filtrage itératif de V; avec le filtre Fp, o p " DT( le résultat est stocké dans V1); 3)- Filtrage itératif de V avec le filtre Fqi o q=D,(le résultat est stocké dans V2); 4)- Seuillage adaptatif de la différence V2V1 en fonction de l'histogramme des valeurs de 0 à 50; et )- Dilatation morphologique en trois dimensions sous contrainte de positivité du gradient local. Le résultat du 15 filtrage est un sous-ensemble de la lumière bronchique, noté par A,. 8. Procédé de modélisation et de reconstruction des voies aériennes selon la revendication 20 3, caractérisé en ce que ladite croissance conditionnelle à une topologie colorée est opérée à partir d'une combinaison de fonctions logiques booléennes validant ou non l'appartenance au réseau bronchique d'un site quelconque x, exprimée par v v{LU w(x)s i,j> 1, o v et A représentent respectivement la somme et le produit logique, et Lî, une fonction logique booléenne définie à partir de l'information de densité V, relative à une topologie donnée. 9. Procédé de modélisation et de reconstruction des voies aériennes selon la revendication 8, caractérisé en ce que LU1 est donné par LI (x) = L(x, x), avec L(x, p) = 3F c Z (p) u z (p), tel que Vg (x) direction locale de croissance et F un chemin fermé. 10. Procédé de modélisation et de reconstruction des voies aériennes selon l'une quelconque des revendications 8 et 9, caractérisé en ce que L12 est donné par: L2 (x)= min Vg(y) reconstruction des voies aériennes selon l'une quelconque des revendications de 8 à 10, caractérisé en ce que L21 est donné par: L2,(x) = L(x,x) A L(p,x) A L(x, p) A L(p, p) A L(v,x), o p et v sont les voisins d'ordre 1 de x sur la direction locale de croissance, v étant le site (occupé) à partir duquel s'effectue la croissance locale. 12. Procédé de modélisation et de reconstruction des voies aériennes selon l'une quelconque des revendications de 8 à 11, caractérisé en ce que L22 est donné par: Z(Vg(y)-V(x)) (V (y) - V(x))2 L22(X) = 1{A} An 6(x) EA926(x) > Ye 26 (X) Y 1{>y.A Z1{y,.A} yE 926(x) Y 926(x) E(Vg (y) - Vg (X)) > 1-iA} -926(x)-A +{Vg(x) Vg(x)> k} >,5L26(X) {NA} YE 26(x) o 26(x) désigne le voisinage 26-connexe de x, 25 A c suppVg l'ensemble courant des points de l'arbre reconstruit, et k une valeur de niveau de gris dans l'intervalle [50, 120]. 13. Procédé de modélisation et de reconstruction des voies aériennes selon l'une quelconque 30 des revendications 1 à 12, caractérisé en ce que l'on opère en outre une indexation de l'arbre bronchique. 14. Procédé de modélisation et de reconstruction des voies aériennes selon l'une quelconque des revendications 1 à 13, caractérisé en ce que l'on opère en outre une extraction de l'axe central. 15. Dispositif de modélisation et de reconstruction des voies aériennes permettant la mise en oeuvre du procédé de modélisation et de reconstruction des voies aériennes selon l'une quelconque des revendications 1 à 8, caractérisé en ce qu'il comporte au moins: a un module de mise en forme desdites données afin d'obtenir un volume pulmonaire en niveaux de gris Vg, * un module d'extraction de l'arbre bronchique à partir dudit volume pulmonaire en niveaux de gris Vg, par croissance morphologique conditionnelle à une 15 topologie colorée, afin d'obtenir un arbre binaire en trois dimensions, et * un module de représentation de l'arbre bronchique, afin de permettre de naviguer dans lesdites voies aériennes.