FR3083358A1 - Systeme et procede pour produire un indicateur graphique multiparametrique a partir d'une image d'une coupe histologique - Google Patents

Systeme et procede pour produire un indicateur graphique multiparametrique a partir d'une image d'une coupe histologique Download PDF

Info

Publication number
FR3083358A1
FR3083358A1 FR1855953A FR1855953A FR3083358A1 FR 3083358 A1 FR3083358 A1 FR 3083358A1 FR 1855953 A FR1855953 A FR 1855953A FR 1855953 A FR1855953 A FR 1855953A FR 3083358 A1 FR3083358 A1 FR 3083358A1
Authority
FR
France
Prior art keywords
interest
quantities
component
estimated
area
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
FR1855953A
Other languages
English (en)
Inventor
Yvon Jule
Florent Tomi
Jean-Claude Gilhodes
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Biocellvia
Original Assignee
Biocellvia
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Biocellvia filed Critical Biocellvia
Priority to FR1855953A priority Critical patent/FR3083358A1/fr
Priority to US17/256,316 priority patent/US20210264596A1/en
Priority to PCT/FR2019/051614 priority patent/WO2020002859A1/fr
Priority to EP19752216.2A priority patent/EP3815037A1/fr
Publication of FR3083358A1 publication Critical patent/FR3083358A1/fr
Ceased legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—Two-dimensional [2D] image generation
    • G06T11/20—Drawing from basic elements
    • G06T11/26—Drawing of charts or graphs
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H15/00—ICT specially adapted for medical reports, e.g. generation or transmission thereof
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30061—Lung
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/40—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for data related to laboratory analysis, e.g. patient specimen analysis
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H70/00—ICT specially adapted for the handling or processing of medical references
    • G16H70/60—ICT specially adapted for the handling or processing of medical references relating to pathologies

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Databases & Information Systems (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Quality & Reliability (AREA)
  • Geometry (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

L'invention concerne un procédé (100) pour produire un indicateur graphique multiparamétrique (I) relatif à un remodelage de l'épithélium bronchique humain ou animal à partir d'une représentation numérique (RDI, MRI) d'une coupe histologique d'un poumon comportant un ou plusieurs composants (Ci) dont chacun décrit une paroi (CWi) enceignant une lumière (CLi). Un tel procédé permet notamment d'estimer une ou plusieurs quantités d'intérêt (QI) relatives à un ou plusieurs composants (Ci) pour produire une représentation graphique décrivant la ou lesdites quantités d'intérêt (QI) afin de fournir une aide au diagnostic d'une pathologie. Le procédé (100) est mis en œuvre par une unité de traitement d'un objet électronique d'un système d'aide au diagnostic d'une pathologie conforme à l'invention.

Description

Système et procédé pour produire un indicateur graphique multiparamétrique à partir d'une image d'une coupe histologique
L'invention concerne un système et un procédé pour produire un indicateur graphique multiparamétrique en lien avec le remodelage de l'épithélium d'un tissu d'un organe humain ou animal, à partir d'une image d'une coupe histologique et ainsi, délivrer une aide objective et reproductible à un personnel de santé pour que ce dernier puisse établir un diagnostic en lien avec une pathologie humaine ou animale éventuelle. L'invention prévoit en outre une aide objective et reproductible pour qu'un expérimentateur en laboratoire puisse estimer une pertinence curative d'un traitement donné au regard d'une telle pathologie.
L'imagerie biologique est à l'heure actuelle une des ressources majeures de l'exploration des organes et des différents tissus organiques. Elle intervient notamment de façon prépondérante dans les domaines de l'aide au diagnostic médical et de la recherche préclinique et clinique.
Différentes techniques sont actuellement mises en œuvre dans l'imagerie préclinique et clinique, telles que, de manière non limitative, l'imagerie à résonance magnétique, la microscopie optique, électronique et confocale, la micro-tomographie, l'échographie, le scanner. Ces techniques peuvent être mises en œuvre pour des observations in vivo ou ex vivo. Les images numériques ainsi obtenues permettent, dans le contexte des laboratoires de recherche institutionnels ou industriels, d'analyser plus particulièrement un état biologique de tissus organiques et d'évaluer certains effets bénéfiques et/ou toxiques de certaines substances en vue de leur sélection pour l'élaboration de futurs médicaments.
A l'ère de la transformation digitale, le développement de ces technologies d'imagerie numérique a ouvert de nouvelles perspectives pour l'analyse histologique dans son ensemble.
La possibilité d'accéder aux images numériques des coupes histologiques a permis de développer de nouvelles méthodes basées sur l'analyse descriptive et quantitative des images numériques desdites coupes histologiques, au moyen d'outils informatiques mettant en œuvre des algorithmes ou procédés innovants permettant une avancée en termes de précision, de fiabilité, de rapidité et de reproductibilité.
Cependant, la mise en œuvre des outils informatiques actuellement disponibles ne permet pas d'automatiser l'évaluation quantitative de certaines pathologies, comme à titre d'exemples non limitatifs, les afflictions des voies respiratoires. En effet, l'expérimentateur demeure encore bien trop présent dans le processus de mise en œuvre de cette évaluation. Son intervention manuelle et personnelle entraîne de grandes variabilités dans la caractérisation des composants des échantillons de lames histologiques expertisées.
Dans le cadre de l'aide au diagnostic de certaines pathologies touchant les voies respiratoires, plus particulièrement les petites voies aériennes, telles que par exemple le remodelage des petites voies aériennes ou « Small Airway Remodeling » selon une terminologie anglo saxonne, également connu sous l'acronyme SAR, il n'existe pas d'évaluation conventionnelle concernant ce genre d'afflictions respiratoires. Chaque laboratoire effectue ses propres mesures à partir de ses propres références, qui se rapportent le plus souvent à quelques mesures morphométriques, telles que des périmètres et/ou diamètres bronchiques. Ces mesures sont effectuées manuellement par un expérimentateur à partir d'images de microscopie optique.
On trouve également d'autres techniques permettant d'identifier et de quantifier un degré de sévérité d'une pathologie liée au SAR. Certaines exploitent des coupes histologiques, d'autres sont exploitées in vivo. Parmi celles-ci, nous pouvons notamment mentionner la spirométrie. Cette première technique connue et largement utilisée, à la fois dans le diagnostic et la
détermination d'un degré de gravité de maladies
pulmonaires, consiste en un test de la fonction
pulmonaire. Un diagnostic concernant une maladie
pulmonaire obstructive est posé lorsqu'un rapport entre un volume expiratoire maximal en une seconde (VEMS) et une capacité vitale forcée (CVF) est inférieur à soixante-dix pour cent. Bien qu'une réduction du VEMS peut refléter une obstruction du flux d'air, elle dépend toutefois et également des volumes pulmonaires, du phénomène de rétractation élastique pulmonaire, de la force musculaire respiratoire et/ou de l'effort du patient. Cette méthode trouve ainsi des limitations importantes quant à l'utilisation des paramètres précédemment évoqués dont la variation peut être induite par des patholoqies autres que le SAR. De telles limitations sont également présentes dans le cadre du diagnostic d'un SAR de faible gravité, et dont le diagnostic par une telle méthode peut s'avérer difficile.
Une deuxième technique consiste à évaluer par pléthysmographie des volumes pulmonaires. Elle délivre une mesure sensible des gaz piégés et de 1'hyperinflation pulmonaire, qui peut être définie comme une élévation anormale des volumes pulmonaires à la fin de l'expiration. Cette mesure donne une fonction de la limitation du débit d'air, du phénomène de rétractation élastique des poumons et de la compliance de la paroi thoracique d'un patient. Un rétrécissement des voies aériennes entraîne en effet une augmentation du temps d'expiration nécessaire pour évacuer l'intégralité de l'air contenu dans les poumons. Il se peut ainsi que les voies respiratoires se ferment, piégeant alors le gaz restant. Le Volume Résiduel (VR) de gaz restant consiste également en une mesure indiquant un dysfonctionnement des petites voies respiratoires. Ledit volume peut être directement corrélé avec le degré de changements morphologiques des voies respiratoires dus aux inflammations présentes dans les petites voies aériennes. Un rapport dudit Volume Résiduel sur la capacité pulmonaire totale (TLC), que nous pouvons noter RV/TLC, peut être facilement établi par cette méthode, au même titre qu'une résistance des voies aériennes qui s'accroît lors de maladies pulmonaires obstructives. Cependant, cet indicateur n'est pas spécifique aux petites voies respiratoires, ce qui limite son application dans le diagnostic et la surveillance des maladies de type SAR.
Également une troisième technique peut être employée, connue sous la dénomination de « tomographie d'émission monophotonique » assistée par ordinateur, ou « SPECT » selon un acronyme anglo-saxon. Une telle technique d'imagerie en trois dimensions consiste à exploiter plusieurs détecteurs de rayons gamma qui évoluent autour d'un patient en position allongée. Une reconstruction d'images obtenues permet alors de visualiser une distribution des radionucléides en trois dimensions, offrant ainsi une évaluation de la ventilation de différentes régions pulmonaires.
Bien que très largement utilisées à l'heure actuelle, ces trois méthodes de mesure présentent bon nombre d'inconvénients. Tout d'abord, elles nécessitent un temps de mise en œuvre long, c'est-à-dire parfois plusieurs heures, et sont dépendantes de l'expérimentateur qui les met en œuvre. Elles impliquent ainsi généralement une analyse complémentaire par un pathologiste spécialisé pour corroborer ou infirmer les premiers résultats. C'est notamment le cas pour la spirométrie. L'implication de plusieurs expérimentateurs ou personnels humains spécialisés induit donc une variabilité importante des résultats concernant l'évaluation d'une pathologie et ainsi des diagnostics parfois aléatoires ou contradictoires.
En outre, la technique d'analyse SPECT peut induire un inconvénient majeur et un risque sanitaire pour un patient lors de la recherche de pathologies, eu égard aux rayonnements appliqués audit patient. Des limitations quant à l'intensité desdites radiations et/ou du temps d'acquisition des signaux expérimentaux peuvent impacter l'évaluation de l'élimination des radionucléides après dépôt de ces derniers dans les différentes régions pulmonaires d'intérêt, obérant ainsi une évaluation pertinente de la ventilation desdites régions pulmonaires du patient.
Ainsi, l'invention permet de répondre à tout ou partie des inconvénients précédemment soulevés et permet d'offrir une aide précieuse à tout expérimentateur désireux d'estimer des quantités d'intérêt en vue de produire un indicateur graphique, pour faciliter l'établissement d'un diagnostic en lien avec une pathologie humaine ou animale, voire pour estimer la pertinence d'un traitement au regard de ladite pathologie.
Certaines maladies pulmonaires chroniques sont caractérisées par le remodelage des tissus pulmonaires (SAR) affectant le parenchyme alvéolaire, les bronches et les vaisseaux. L'invention permet une analyse morphométrique quantitative d'une grande précision des composants observés dans une coupe histologique, permettant ainsi de caractériser, par exemple, un remodelage des voies respiratoires, et permet ainsi de délivrer un indicateur graphique multiparamétrique en lien avec la morphologie de composants tubulaires d'un organe, notamment des poumons, et plus particulièrement des petites voies respiratoires. L'invention exploite une représentation numérique d'une coupe histologique dudit organe, décrivant une pluralité de composants annulaires résultant de la section desdits composants tubulaires pour pratiquer une coupe histologique. Lesdits composants tubulaires se manifestent par des composants annulaires sur une représentation en deux dimensions, à l'issue d'une coupe histologique. La notion de composant tubulaire sera décrite plus en détail par la suite en lien avec la figure 3.
Parmi les nombreux avantages procurés par l'invention, nous pouvons mentionner que celle-ci permet de :
- réduire considérablement le temps d'analyse nécessaire à l'établissement d'un diagnostic d'une pathologie par un expérimentateur, réduisant ledit temps à moins d'une minute suivant la puissance de calcul du dispositif ou du système électronique mettant en œuvre un procédé conforme à l'invention ;
- décupler la précision et la fiabilité des mesures d'un échantillon analysé ;
- s'affranchir de la variabilité de résultats entre différents expérimentateurs, délivrant des mesures objectives et reproductibles ;
- ne pas requérir de matériels coûteux et/ou de procédés nocifs ou invasifs pour les patients examinés.
A cette fin, l'invention concerne un procédé pour produire un indicateur graphique multiparamétrique relatif à un remodelage de l'épithélium bronchique humain ou animal à partir d'une représentation numérique d'une coupe histologique d'un poumon comportant un ou plusieurs composants de forme annulaire dont chacun décrit une paroi enceignant une lumière, le procédé étant mis en œuvre par une unité de traitement d'un système d'aide au diagnostic d'une pathologie, ledit système comportant en outre une interface homme-machine de sortie et une mémoire de données.
Pour pouvoir produire un indicateur graphique multiparamétrique, ledit procédé comporte :
- une étape pour estimer des quantités d'intérêts relatives aux morphologies respectives d'un ou plusieurs composants présents dans la coupe histologique, lesdites quantités d'intérêt appartenant à un ensemble de quantités d'intérêt comprenant :
i. le diamètre de Féret du contour extérieur de la paroi d'un composant ;
ii. une distance moyenne séparant ledit contour extérieur du contour intérieur, traduisant une épaisseur moyenne de la paroi d'un composant ;
iii. l'aire décrite par ledit contour intérieur d'un composant traduisant la lumière de ce dernier ;
iv. l'aire décrite par ledit contour extérieur d'un composant traduisant l'aire totale couverte par ce dernier ;
v. l'aire de la paroi d'un composant définie par la soustraction de l'aire décrite par ledit contour intérieur à l'aire décrite par ledit contour extérieur dudit composant ;
- une étape pour produire, par quantité d'intérêt estimée, une représentation graphique de celle-ci relative à une quantité d'intérêt étalon ;
- une étape pour provoquer la restitution graphique conjointe des représentations graphiques desdites quantités d'intérêt relatives aux quantités d'intérêt étalon respectives préalablement produites, par l'interface homme-machine de sortie du système.
Avantageusement, l'étape pour estimer des quantités d'intérêts relatives aux morphologies respectives d'un ou plusieurs composants peut comporter une étape pour inscrire dans la mémoire de données une structure de données associée à chaque composant, ladite structure de données comportant un champ pour enregistrer la valeur de chaque quantité d'intérêt estimée.
En variante ou en complément, afin de produire un indicateur graphique multiparamétrique relatif à une ou plusieurs quantités d'intérêt d'un ou plusieurs types de composants, le procédé peut comporter une étape pour caractériser un type de composant à partir de la valeur d'une des quantités d'intérêt estimées, l'étape pour inscrire dans la mémoire de données une structure de données associée à chaque quantité d'intérêt estimée d'un composant peut consister à inscrire dans un champ de ladite structure de données une valeur caractérisant un type de composant déterminé.
Préférentiellement mais non limitativement, afin de faciliter le typage des composants présents dans une représentation numérique d'une coupe histologique, lorsque l'une des quantités d'intérêt estimées consiste en le diamètre de Féret du contour extérieur de la paroi du composant, l'étape pour caractériser un type de composant d'un procédé conforme à l'invention peut comporter une opération de comparaison de la valeur dudit diamètre de Féret estimée à un seuil de diamètre haut et/ou un seuil de diamètre bas.
En variante ou en complément, afin de faciliter le typage des composants présents dans une représentation numérique d'une coupe histologique, lorsque l'une des quantités d'intérêt estimées consiste en l'épaisseur moyenne de la paroi du composant, l'étape pour caractériser un type de composant d'un procédé conforme à l'invention peut comporter une opération de comparaison de la valeur de ladite épaisseur moyenne à un seuil d'épaisseur haut et/ou un seuil d'épaisseur bas prédéterminés.
En variante ou en complément, pour faciliter le typage des composants présents dans une représentation numérique d'une coupe histologique, lorsque l'une des quantités d'intérêt estimées consiste en l'aire décrite par le contour intérieur du composant, l'étape pour caractériser un type de composant d'un procédé conforme à l'invention peut comporter une opération de comparaison de la valeur de ladite aire à un seuil d'aire haut et/ou un seuil d'aire bas prédéterminés.
Afin de fournir à un expérimentateur des indications relatives à la morphologie d'un ou plusieurs composants d'intérêt, l'étape pour produire, par quantité d'intérêt estimée, une représentation graphique de celle-ci relative à une quantité d'intérêt étalon et/ou l'étape pour provoquer la restitution graphique conjointe des représentations graphiques desdites quantités d'intérêt relatives aux quantités d'intérêt étalon ne peuvent être mises en œuvre que pour un type de composant caractérisé déterminé.
De manière avantageuse, afin de permettre à un expérimentateur d'identifier un remodelage de l'épithélium bronchique ou toute autre modification de la morphologie des petites voies aériennes, un procédé selon l'invention prévoit que le type de composant caractérisé déterminé peut être une bronchiole.
Préférentiellement mais non limitativement, afin de fournir à un expérimentateur une aide visuelle permettant de confronter instantanément des quantités d'intérêts estimées, traduisant un éventuel changement de morphologie des composants dont sont issues lesdites quantités d'intérêt estimées en vue d'orienter le diagnostic d'une pathologie, l'étape, pour provoquer la restitution graphique conjointe des représentations graphiques desdites quantités d'intérêt relatives aux quantités d'intérêt étalon respectives préalablement produites par l'interface homme-machine de sortie du système, peut consister en l'affichage par cette dernière d'un diagramme de Kiviat. Un tel diagramme de Kiviat peut décrire au moins trois représentations graphiques de quantités d'intérêt relatives aux quantités d'intérêt étalon respectives sur des axes normés.
En variante ou en complément, afin de fournir à un expérimentateur une aide visuelle permettant de comparer instantanément des quantités d'intérêts estimées traduisant un éventuel changement de morphologie des composants dont sont issues lesdites quantités d'intérêt estimées, en vue d'orienter le diagnostic d'une pathologie, l'étape pour provoquer la restitution graphique conjointe des représentations graphiques desdites quantités d'intérêt relatives aux quantités d'intérêt étalon respectives préalablement produites par l'interface homme-machine de sortie du système peut consister en l'affichage par cette dernière d'un diagramme en barres décrivant les représentations graphiques de quantités d'intérêt relatives aux quantités d'intérêt étalon respectives par des barres normées.
De manière à restituer à un expérimentateur un indicateur graphique multiparamétrique présentant des quantités d'intérêt à une même échelle et ainsi fournir une aide visuelle instantanée, la normalisation d'un axe ou d'une barre de l'étape pour provoquer la restitution graphique conjointe des représentations graphiques desdites quantités d'intérêt relatives aux quantités d'intérêt étalon respectives préalablement produites par l'interface homme-machine de sortie du système peut consister à exprimer la valeur d'une quantité d'intérêt estimée en pourcentage de la valeur de la quantité d'intérêt étalon associée.
En variante ou en complément, pour permettre à un expérimentateur de visualiser instantanément une différence significative entre des quantités d'intérêts estimées et des quantités d'intérêts étalon très proches, la normalisation d'un axe ou d'une barre de l'étape pour provoquer la restitution graphique conjointe des représentations graphiques desdites quantités d'intérêt relatives aux quantités d'intérêt étalon respectives préalablement produites par l'interface homme-machine de sortie du système peut consister à exprimer la valeur d'une quantité d'intérêt estimée relativement à la valeur de la quantité d'intérêt étalon associée sous la forme de trois valeurs prédéterminées décrivant respectivement des valeurs de quantité d'intérêt estimées sensiblement inférieures, proches ou supérieures aux valeurs des quantités d'intérêt étalon associées.
Selon un deuxième objet, l'invention concerne en outre un objet électronique d'un système d'aide au diagnostic d'une pathologie, ledit objet électronique comportant une unité de traitement et coopérant avec une interface homme-machine de sortie et avec une mémoire de données, ladite mémoire de données comportant :
une représentation numérique d'une coupe histologique d'un organe humain ou animal ; des instructions exécutables ou interprétables par les moyens de traitement dont l'interprétation ou l'exécution desdites instructions par lesdits moyens de traitement provoque la mise en œuvre d'un procédé conforme au premier objet de 1'invention.
Selon un troisième objet, l'invention concerne également un système d'aide au diagnostic d'une pathologie comportant un objet électronique selon le deuxième objet de l'invention et une interface hommemachine de sortie apte à restituer à un utilisateur un indicateur graphique multiparamétrique selon un procédé conforme au premier objet de l'invention et mis en œuvre par ledit objet électronique.
Selon un quatrième objet, l'invention concerne enfin un produit programme d'ordinateur comportant une ou plusieurs instructions interprétables ou exécutables par un objet électronique conforme au deuxième objet de l'invention, l'interprétation ou l'exécution desdites instructions par ladite unité de traitement provoque la mise en œuvre d'un procédé conforme au premier objet de 1'invention.
D'autres caractéristiques et avantages apparaîtront plus clairement à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent parmi lesquelles :
- la figure 1 présente une première représentation numérique d'une coupe histologique d'un poumon, d'un sujet atteint d'un SAR, ledit sujet étant en l'espèce une souris ;
- la figure 2 illustre une deuxième représentation numérique binaire, issue de celle présentée par la figure 1, ladite deuxième représentation numérique binaire mettant en exergue des pixels d'intérêt au regard des autres pixels ;
- la figure 3 schématise une section d'un composant tubulaire, à partir duquel sont estimées des quantités d'intérêt relatives à la morphologie d'un composant annulaire résultant de ladite section d'un composant tubulaire ;
- la figure 4A présente un premier exemple de représentation graphique sous la forme d'un diagramme en barres, d'un indicateur graphique multiparamétrique produit par un procédé conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu ;
- la figure 4B présente un deuxième exemple de représentation graphique sous la forme d'un diagramme de Kiviat, d'un indicateur graphique multiparamétrique produit par un procédé conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu ;
- la figure 5 présente un premier exemple de représentation graphique sous la forme d'un diagramme en barres, d'un indicateur graphique multiparamétrique produit par un procédé conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu ;
- la figure 6 présente un deuxième exemple de représentation graphique sous la forme d'un diagramme de Kiviat, d'un indicateur graphique multiparamétrique produit par un procédé conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu ;
- la figure 7A présente un premier exemple de représentation graphique sous la forme d'un diagramme en barres, d'un indicateur graphique multiparamétrique produit par un procédé conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu ;
- la figure 7B présente un deuxième exemple de représentation graphique sous la forme d'un diagramme de Kiviat, d'un indicateur graphique multiparamétrique produit par un procédé conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu ;
- la figure 8 présente un premier exemple de représentation graphique sous la forme d'un diagramme en barres, d'un indicateur graphique multiparamétrique produit par un procédé conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu ;
- la figure 9 présente un deuxième exemple de représentation graphique sous la forme d'un diagramme de Kiviat, d'un indicateur graphique multiparamétrique produit par un procédé conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu ;
- la figure 10 présente un organigramme simplifié illustrant un exemple non limitatif d'un procédé pour produire un indicateur graphique multiparamétrique relatif à un tissu d'un organe humain ou animal conforme à l'invention.
La figure 1 décrit une première représentation numérique RDI d'une coupe histologique d'un poumon d'un sujet atteint, par exemple, d'un remodelage des petites voies aériennes ou SAR. Ledit sujet atteint est en l'espèce une souris. Une telle première représentation RDI est généralement issue d'un processus de numérisation d'une coupe histologique. Une section histologique numérisée avec un agrandissement x20 délivre ladite première représentation numérique sous une forme matricielle d'environ deux cents millions de pixels, soit selon l'exemple de la figure 1, une représentation sous la forme d'un tableau de quinze mille lignes sur autant de colonnes, chaque élément dudit tableau encodant un triplet de valeurs entières comprises entre zéro et deux cents cinquante-cinq, selon le codage de couleurs RVB (acronyme pour « Rouge Vert Bleu »), également connu sous le sigle anglo-saxon RGB (acronyme pour « Red Green Blue »). Un tel codage informatique des couleurs est le plus proche des matériels disponibles. En général, les écrans d'ordinateurs reconstituent une couleur par synthèse additive à partir de trois couleurs primaires, un rouge, un vert et un bleu, formant sur l'écran une mosaïque généralement trop petite pour être discriminée par l'humain. Le codage RVB indique une valeur pour chacune de ces couleurs primaires. Une telle valeur est généralement codée sur un octet et appartient donc à un intervalle de valeurs entières comprises entre zéro et deux cents cinquante-cinq.
La représentation RDI présente, de manière distinguable, au centre de ladite représentation, le lobe L d'un poumon. Un tel organe comporte, en outre, de nombreux composants Ci sensiblement tubulaires dont les parois internes forment respectivement des lumières. A la suite de la coupe pratiquée au niveau du lobe pulmonaire, seules des sections desdits composants tubulaires Ci sont visibles en deux dimensions. Ces dernières peuvent être décrites, notamment en lien avec la figure 3, comme des structures annulaires, dont l'anneau est une section de la paroi du composant Ci, cerclant un trou associé à la lumière. De tels composants tubulaires Ci consistent principalement en des vaisseaux, des bronches, des bronchioles ou encore des alvéoles pouvant former des sacs alvéolaires. Le reste du tissu P dudit lobe est nommé ci-après « parenchyme ».
Lorsqu'un sujet est atteint, par exemple, d'une maladie respiratoire obstructive chronique, la lésion du poumon se traduit par des changements de la morphologie de certains des composants tubulaires formant le poumon, notamment par un remodelage des voies aériennes, plus particulièrement des bronchioles dans le cas du SAR. Les diamètres des bronchioles sont généralement compris, au sein du taxon des vertébrés, entre cent microns et mille microns. Un SAR implique ainsi, dans les premiers temps d'une atteinte des voies aériennes d'un sujet, une desquamation de l'épithélium bronchique, plus communément appelé « paroi bronchique », se caractérisant sur une représentation numérique d'une coupe histologique d'un poumon, par un anneau représentant une section d'une telle paroi bronchique. Une telle desquamation induit généralement une diminution de l'aire des bronchioles et de ladite paroi bronchique, se caractérisant sur une représentation numérique d'une coupe histologique d'un poumon, par une diminution de l'aire de l'anneau et éventuellement de l'aire cerclée par ledit anneau, c'est-à-dire l'aire de la section d'une lumière dudit composant tubulaire. Dans un second temps, une atteinte plus sévère induit un remodelage de telles bronchioles et une modification du métabolisme de celles-ci, se matérialisant par une inflammation chronique de la paroi bronchique entraînant des changements dans la morphologie desdites bronchioles. Il en résulte alors un accroissement de la lumière moyenne et de l'aire desdites bronchioles remodelées, ainsi qu'une augmentation de l'épaisseur moyenne de l'épithélium bronchique de ces dernières. Une telle atteinte se caractérise, sur une représentation numérique d'une coupe histologique d'un poumon, par une augmentation de l'aire de la structure annulaire, c'està-dire de l'aire de la section dudit composant tubulaire, et plus particulièrement de l'aire cerclée par ledit anneau, c'est-à-dire de l'aire de la section transversale d'une lumière dudit composant tubulaire. La figure 1 décrit ainsi une première représentation numérique RDI, d'un lobe L d'un poumon OG d'un sujet atteint d'un SAR.
Ladite première représentation numérique RDI est, selon l'état de l'art, difficilement exploitable pour qu'un expérimentateur puisse déterminer la présence de modifications morphométriques induites par un remodelage éventuel des petites voies aériennes, caractérisant une pathologie telle que le SAR. En effet, lesdites modifications morphométriques étant de l'ordre du micron, il est très difficile pour un expérimentateur de les constater à partir d'une coupe histologique. Ainsi, selon l'état de l'art, les techniques décrites précédemment sont préférées pour caractériser, chez des sujets, l'atteinte des petites voies aériennes.
Pour produire une aide objective, automatique et quasiment en temps réel à de tels expérimentateurs, l'invention prévoit de se concentrer sur des propriétés morphométriques des composants retrouvés dans une coupe histologique d'un poumon par l'exploitation d'une représentation numérique, obtenue par numérisation de ladite coupe histologique. De manière à simplifier la compréhension du procédé pour estimer une quantité d'intérêt selon l'invention, une seule représentation numérique d'une coupe histologique d'un sujet présentant une atteinte des petites voies aériennes est représentée, les différences avec une représentation numérique d'une coupe histologique d'un sujet sain étant imperceptibles par l'œil humain.
En outre, pour faciliter la recherche de composants d'intérêt, un procédé selon l'invention peut comporter une étape 10 préalable consistant en un traitement pour binariser ladite première représentation numérique RDI et produire une deuxième représentation numérique MRI. Une telle deuxième représentation numérique MRI peut être produite par tout type de traitement numérique connu visant à binariser une représentation numérique en couleur(s), telle que la représentation numérique RDI. Une telle représentation numérique ou image se traduit sous la forme d'un tableau comportant un même nombre d'éléments ou pixels que la première représentation numérique RDI d'une coupe histologique dont elle est issue, telle que la première représentation numérique RDI précédemment évoquée en lien avec la figure 1. Ladite deuxième représentation numérique MRI est dite binaire car chacun de ses éléments MRI(i,j), désigné par deux indexes i et j déterminant respectivement la ligne et la colonne dudit élément ou pixel dans le tableau MRI, comporte une valeur entière choisie parmi deux valeurs prédéterminées signifiant respectivement que le pixel RDI(i,j), c'est-à-dire d'une même colonne j et d'une même ligne i dans une première représentation numérique RDI, désigne ou non une partie du lobe L.
A titre d'exemple non limitatif, la figure 2 illustre un exemple de deuxième représentation binaire MRI, la deuxième représentation numérique binaire MRI étant en l'espèce issue de la première représentation numérique RDI décrite en lien avec la figure 1. Selon cet exemple, un élément MRI(i,j) du tableau MRI prend la valeur zéro si le pixel associé RDI (i,j), c'est-à-dire désigné par la ligne i et la colonne j dans la première représentation numérique RDI, ne correspond pas à un pixel du lobe L, c'est-à-dire que ledit pixel décrit soit une lumière formée par la coupe transversale d'un composant tubulaire, soit l'extérieur du lobe L du poumon. Dans le cas contraire, un tel élément MRI(i,j) prend la valeur deux cents cinquante-cinq. De cette façon, une telle deuxième représentation numérique binaire MRI peut être affichée en noir et blanc sur un écran d'ordinateur. L'invention ne saurait se limiter à l'emploi desdites valeurs zéro et deux cents cinquantecinq. En variante, d'autres valeurs prédéterminées auraient pu être choisies desdites valeurs zéro et deux cents cinquante-cinq pour caractériser l'absence d'intérêt ou l'intérêt d'un tel pixel.
De manière avantageuse mais non limitative, afin de faciliter l'estimation de quantités d'intérêt relatives à des composants annulaires présents dans une représentation numérique d'une coupe histologique d'un tissu d'un organe humain ou animal et in fine la production d'un indicateur graphique multiparamétrique, un traitement 10 pour produire une représentation binaire MRI d'une représentation numérique RDI, peut être mis en œuvre avant la mise en œuvre d'un procédé pour produire un indicateur multiparamétrique selon l'invention, tel qu'un procédé 100, dont un exemple non limitatif de réalisation est notamment décrit en lien avec la figure 10. Par mesure de concision et de simplification, nous parlerons de morphologie du composant Ci en lieu et place de la morphologie de la section de forme annulaire d'un composant tubulaire, telle que décrite en lien avec la figure 3. Ainsi, par application de la coupe histologique, un composant tubulaire d'intérêt Ci se traduit dans une représentation numérique RDI ou MRI par une structure annulaire en deux dimensions. Ledit traitement 10 visant à binariser une représentation numérique peut, en variante, constituer une première étape dudit procédé 100 évoqué précédemment. Un exemple non limitatif d'un traitement 10, visant à binariser une telle première représentation numérique RDI, peut comporter une première étape pour produire une première représentation numérique intermédiaire en nuances de gris, non illustrée par les figures à des fins de simplification. Ladite représentation numérique intermédiaire comporte un même nombre d'éléments ou pixels que la première représentation numérique RDI. Une telle première étape peut consister en la mise en œuvre de toute technique connue pour convertir, pour chaque pixel de la représentation RDI, le triplet de valeurs représentant les niveaux des couleurs primaires en une valeur entière représentant une luminosité ou une intensité lumineuse associée à un pixel de représentation ainsi produite. Ladite première étape peut en outre consister en l'application, sur la représentation numérique ainsi produite, d'un filtre médian ou bilatéral pour supprimer certaines aberrations.
Une deuxième étape d'un tel traitement 10 visant à binariser une représentation numérique selon l'invention peut consister à mettre en œuvre un seuillage automatique des pixels de la (première) représentation numérique intermédiaire, de manière à discriminer les pixels décrivant tout ou partie d'une lumière formée par la coupe transversale d'un composant tubulaire ou l'extérieur du lobe L du poumon. Les pixels associés à une lumière ou à l'extérieur du lobe prennent la valeur zéro, apparaissant ainsi en noir sur la figure 2. Les autres pixels prennent la valeur deux cents cinquantecinq et apparaissent en blanc. Ils sont associés à du tissu parenchymateux ou à certains composants tubulaires tels que des vaisseaux, alvéoles, sacs alvéolaires, bronches ou encore bronchioles. Cette deuxième étape produit ainsi une deuxième représentation numérique binaire MRI.
Un tel traitement 10 visant à binariser une représentation numérique peut, en complément, produire une troisième représentation numérique binaire, que nous nommerons « Masque de lobe », de mêmes dimensions que la deuxième représentation numérique binaire MRI, dont chaque élément comporte une première valeur spécifiant qu'un pixel associé au sein d'une représentation numérique RDI ou MRI appartient au lobe ou est extérieur à celui-ci. En effet, tenir compte notamment dans la deuxième représentation numérique binaire MRI de pixels associés à l'arrière-plan AP du lobe L peut altérer la pertinence des quantités d'intérêt produites et in fine la pertinence d'un indicateur graphique multiparamétrique produit par un procédé selon l'invention. Pour cela, dans un exemple non limitatif d'un procédé conforme à l'invention, un pixel de la deuxième représentation numérique binaire MRI ne sera pris en considération que si et seulement si, le pixel associé dans ledit masque de lobe, c'est-à-dire de mêmes index de ligne et de colonne, comporte une valeur caractérisant un pixel appartenant au lobe L examiné. Une telle troisième représentation numérique, non représentée par les figures à des fins de simplification, peut être produite en complément de la deuxième étape d'un traitement 10 visant à binariser une représentation numérique, par la recherche du plus grand contour par la mise en œuvre, par exemple, d'un algorithme « Flood Fill » selon une terminologie anglo-saxonne ou encore connu par l'appellation francophone « algorithme par remplissage par diffusion ».
La figure 10 décrit un exemple de réalisation non limitatif d'un procédé 100 pour produire un indicateur graphique multiparamétrique I relatif à un remodelage de l'épithélium d'un organe humain ou animal à partir d'une représentation numérique d'une coupe histologique dudit organe conforme à l'invention, qu'une telle représentation soit sous la forme d'une représentation numérique en couleur RDI ou binaire MRI.
Un tel procédé 100 et/ou un tel traitement préalable 10, visant à binariser une représentation numérique RDI et produire une représentation numérique MRI, peuvent être agencés pour être transcrits en un programme d'ordinateur, dont les instructions de programme peuvent être installées dans une mémoire de programmes d'un objet électronique, sous la forme par exemple d'un ordinateur disposant d'une puissance de calcul suffisante et/ou adaptée à l'analyse de représentations numériques ou d'images de tailles conséquentes, compte tenu de la précision nécessaire à l'analyse d'un lobe pulmonaire L.
Lesdites instructions de programme sont ainsi agencées pour provoquer la mise en œuvre desdits procédé 100 et traitement préalable 10 visant à binariser une représentation numérique RDI par l'unité de traitement d'un tel objet électronique. Au sens du présent document, on entend par « unité de traitement », un ou plusieurs microcontrôleurs ou microprocesseurs coopérant avec une mémoire de programmes hébergeant le programme d'ordinateur selon l'invention. Une telle unité de traitement peut en outre être agencée pour coopérer avec une mémoire de données pour héberger, c'est-à-dire enregistrer, les représentations numériques produites par la mise en œuvre d'un procédé 100 pour produire un indicateur graphique multiparamétrique I relatif à un remodelage de l'épithélium bronchique selon l'invention et/ou toutes autres données nécessaires à la mise en œuvre de ce dernier.
Une telle unité de traitement peut également être agencée pour coopérer avec une interface homme-machine ou périphérique de sortie, tel qu'un écran d'ordinateur, une imprimante ou toute autre interface pour délivrer la teneur dudit indicateur graphique multiparamétrique I à un humain, de manière perceptible par l'intermédiaire de l'un de ses sens.
En lien avec la figure 10, un tel procédé 100 pour produire un indicateur graphique multiparamétrique I peut comporter une séquence 110 éventuellement itérative, d'étapes visant à estimer, à partir d'une première représentation numérique RDI et/ou MRI d'une coupe histologique d'un tissu d'un organe humain ou animal, une ou plusieurs quantités d'intérêts QI relatives aux morphologies respectives des composants Ci présents dans ladite coupe histologique.
Ladite séquence itérative d'un procédé pour produire un indicateur graphique multiparamétrique I conforme à l'invention comporte une étape 111 pour estimer une ou plusieurs quantités d'intérêts QI. La mise en œuvre de l'étape 111 peut avantageusement comporter une succession de sous-étapes, non représentées sur la figure 10, consistant à « rechercher », à partir d'une représentation numérique RDI et/ou MRI d'une coupe histologique, une première lumière décrite par la coupe d'un composant sensiblement annulaire Ci. Par l'application d'une technique, telle que, celle décrite par Satochi Suzuki et al. « Topological structural analysis of digitized binary images by border following, Computer Vision, Graphics, and Image processing, 1985 » ou de toute autre technique équivalente, l'étape 111 consiste à déterminer, à partir de représentations numériques RDI et/ou MRI, des topologies de composants Ci présents dans ladite représentation numérique RDI et/ou MRI. Il est alors possible de déterminer le contour d'une aire associée à une lumière d'un composant Ci précédemment identifié, celle-ci s'apparentant sur une représentation numérique binaire MRI à un ensemble de pixels contigus décrivant par exemple un nombre de valeurs égales à zéro. L'étape 111 consiste alors à délimiter le contour intérieur de ladite lumière, c'està-dire par voie de conséquence, celui de la paroi interne dudit composant identifié. Le résultat de l'application d'une telle technique se traduit par une première polyligne dont les index, c'est-à-dire les colonnes et lignes, des pixels qui en constituent les points caractéristiques, décrivent ledit contour de la paroi interne, enceignant la lumière précédemment identifiée, d'un composant tubulaire Ci. L'étape 111 consiste à présent à déterminer, à partir d'une telle première polyligne, le contour de la paroi extérieure du composant Ci précédemment identifié. A titre d'exemple non limitatif, une telle étape 111 pour estimer une ou plusieurs quantités d'intérêt d'un tel procédé 100 pour produire un indicateur graphique multiparamétrique I peut mettre en œuvre une technique, telle que celle connue sous l'appellation « dilatation morphologique du contour interne » décrite dans l'ouvrage de Jean Serra, Image Analysis and Mathematical Morphologyr 1982, ou encore toute autre technique équivalente. La mise en œuvre d'une telle technique produit ainsi une deuxième polyligne dont les index des pixels qui en constituent les points caractéristiques permettent, conjointement avec ladite première polyligne précédemment produite, d'estimer la morphologie de la section d'un composant de forme annulaire définie par les première et deuxième polylignes ainsi déterminées. A partir desdites première et deuxième polylignes, l'étape 111 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention consiste à présent à estimer une quantité d'intérêt QI relative à ladite morphologie dudit composant Ci dont la lumière a précédemment été identifiée.
Un composant d'intérêt Ci présente ainsi, en deux dimensions, une forme sensiblement annulaire telle que décrite dans un exemple en lien avec la figure 3, avantageusement mais non limitativement pour un type CTi de composant Ci, en l'espèce une bronchiole. Ladite bronchiole, et plus généralement un type de composant Ci, comprend un contour extérieur CWOCi, un contour intérieur CWICi enceignant une lumière CLi, et une paroi CWi, dont les aires correspondantes sont désignées par une flèche à bout rond. Comme précisé précédemment, différentes quantités d'intérêt QI peuvent être estimées à l'étape 111 d'un procédé conforme à l'invention. Telles que décrites en lien avec la figure 3, de telles quantités d'intérêt peuvent consister en :
- le diamètre de Féret DFi d'un composant Ci préalablement identifié. Au sens de l'invention et dans tout le document, on entend par « diamètre de Féret », la plus grande distance entre deux tangentes à la deuxième polyligne décrivant le contour extérieur CWOCi apparent de la paroi CWi dudit composant Ci, lesdites tangentes étant parallèles entre elles et perpendiculaires à un vecteur de direction donnée. En faisant varier l'angle polaire dudit vecteur de 0 à 2 π, lors de la mise en œuvre de l'étape 111, un ensemble de diamètres peuvent être ainsi estimés, dont le plus grand correspond au diamètre de Féret DFi retenu pour caractériser ledit composant Ci ;
- La distance moyenne CWWi, que l'on nommera épaisseur moyenne CWWi par mesure de simplification, de la paroi CWi d'un composant Ci préalablement identifié. Au sens de l'invention et dans tout le document, on entend par « épaisseur moyenne CWWi », le rapport entre la somme des distances estimées séparant le contour extérieur CWOCi du contour intérieur CWICi dudit composant Ci préalablement identifié et le nombre de distances estimées. Dans le cadre de pathologies touchant les voies aériennes, telles que le SAR, la paroi CWi de certains composants Ci retrouvés dans les poumons peut subir une desquamation totale ou partielle. Ainsi, l'estimation d'une telle distance moyenne CWWi d'un composant Ci est particulièrement pertinente, puisqu'elle permet de décrire l'importance et l'étendue d'une telle desquamation de la paroi CWi d'un composant annulaire Ci et par voie de conséquence la desquamation du composant tubulaire Ci, dont le composant annulaire Ci en décrit une coupe ;
- l'aire ICAi décrite par le contour intérieur CWICi d'un composant Ci traduisant la lumière CLi de ce dernier. Au sens de l'invention et dans tout le document, on entend par « l'aire ICAi », la surface d'une section décrivant la lumière CLi d'un composant Ci, ladite surface étant délimitée par la première polyligne dudit composant Ci ;
- l'aire OCAi décrite par le contour extérieur CWOCi d'un composant Ci, traduisant l'aire totale couverte par ce dernier. Au sens de l'invention et dans tout le document, on entend par « l'aire OCAi », la surface décrivant la lumière CLi d'un composant Ci, ladite surface étant délimitée par la première polyligne dudit composant Ci, ainsi que la surface décrite par la paroi dudit composant Ci, ladite surface décrivant la paroi étant délimitée par la première polyligne et la deuxième polyligne dudit composant Ci ;
- l'aire CAi décrite par la paroi CWi d'un composant Ci, traduisant l'aire couverte par celle-ci. Au sens de l'invention et dans tout le document, on entend par « l'aire CAi » la surface décrivant la paroi CWi, ladite paroi CWi étant délimitée par les première et deuxième polyligne dudit composant Ci.
A titre d'exemple avantageux mais non limitatif, la mise en œuvre de l'étape 111 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention permet, à partir des première et deuxième polylignes précédemment produites, respectivement associées au contour intérieur CWICi et au contour extérieur CWOCi, d'estimer une ou plusieurs quantités d'intérêts QI, telles que le diamètre de Féret DFi d'un composant Ci. Un exemple non limitatif d'une telle étape 111 pour estimer un tel diamètre consiste à sommer les pixels de dimensions connues, afin d'en déduire une distance correspondante caractérisant ledit diamètre. La distance relative à chaque diamètre peut ainsi être estimée, et la plus grande distance correspondant au diamètre de Féret DFi pourra ainsi être inscrite lors de la mise en œuvre d'une étape 113 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention, dans une structure de données associée au composant Ci.
L'étape 111 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut, en variante ou en complément, consister à estimer une distance moyenne CWWi traduisant une épaisseur moyenne de la paroi CWi d'un composant Ci préalablement identifié. Un exemple d'une telle étape 111 pour estimer une telle distance moyenne CWWi peut consister à estimer, pour chaque pixel d'une première et/ou d'une deuxième polyligne décrivant respectivement le contour intérieur CWICi et le contour extérieur CWOCi apparent de la paroi CWi dudit composant Ci, la plus petite distance séparant ledit pixel desdites première et/ou deuxième polylignes. Une pluralité d'estimations de distances décrivant l'épaisseur de la paroi CWi du composant Ci correspondant peuvent être ainsi estimées, de sorte qu'en moyennant lesdites estimations, la mise en œuvre d'une telle étape 111 permet d'estimer une distance moyenne CWWi traduisant l'épaisseur moyenne de la paroi CWi dudit composant Ci. Ladite distance moyenne
pourra, elle aussi, être inscrite lors de la mise en
œuvre d'une étape 113 d'un procédé 100 pour produire un
indicateur graphique multiparamétrique I conforme à
1'invention, dans une structure de données associée au
composant Ci
L'étape 111 d'un procédé 100 pour produire un
indicateur graphique multiparamétrique I conforme à
l'invention peut, en variante ou en complément, consister à estimer une aire ICAi traduisant la surface de la lumière CLi d'un composant Ci préalablement identifié. Un exemple d'une telle étape 111 pour estimer une telle aire ICAi peut consister à sommer les pixels de dimensions connues, décrivant la surface délimitée par ladite première polyligne. La mise en œuvre d'une telle étape 111 permet ainsi d'estimer l'aire ICAi, correspondant à la somme des aires des pixels décrivant la lumière CLi, du composant Ci préalablement identifié. Ladite aire ICAi pourra, elle aussi, être inscrite lors de la mise en œuvre d'une étape 113 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention, dans une structure de données associée au composant Ci.
L'étape 111 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut en outre consister à estimer une aire OCAi traduisant la surface d'un composant Ci préalablement identifié. Un exemple d'une telle étape 111 pour estimer une telle aire OCAi peut consister à sommer les pixels de dimensions connues décrivant l'aire délimitée par ladite deuxième polyligne. La mise en œuvre d'une telle étape 111 permet d'estimer l'aire OCAi, correspondant à la somme des aires des pixels décrivant la totalité de l'aire du composant Ci préalablement identifié. Ladite aire OCAi pourra, elle aussi, être inscrite lors de la mise en œuvre d'une étape 113 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention, dans une structure de données associée au composant Ci.
De manière avantageuse, l'étape 111 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut consister à estimer une aire CAi traduisant la surface de la paroi CWi d'un composant Ci préalablement identifié. Un exemple d'une telle étape 111 pour estimer une telle aire CAi peut consister à soustraire l'aire ICAi décrite par le contour intérieur CWICi d'un composant Ci, à l'aire OCAi décrite par le contour extérieur CWOCi du même composant Ci, et précédemment estimées. Ladite aire CAi pourra, elle aussi, être inscrite lors de la mise en œuvre d'une étape 113 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention, dans une structure de données associée au composant Ci.
Comme d'ores et déjà mentionné, la séquence 110 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut également comprendre une étape 113 pour inscrire dans la mémoire de données, une structure de données associée à chaque composant Ci identifié et dont la morphologie a été caractérisée lors d'une étape 111 dudit procédé 100.
Ladite structure de données peut avantageusement comprendre un champ pour enregistrer la valeur de chaque quantité d'intérêt QI estimée correspondante en lien avec ledit composant Ci, telle que, par exemple, le diamètre de Féret DFi dudit composant Ci, une épaisseur moyenne CWWi de la paroi CWi dudit composant Ci, l'aire ICAi décrivant la lumière CLi dudit composant Ci, l'aire OCAi décrivant la surface totale dudit composant Ci et l'aire CAi décrivant la surface totale de la paroi CWi dudit composant Ci.
De manière avantageuse mais non limitative, tel que d'ores et déjà précisé, ladite séquence peut-être mise en œuvre de manière itérative pour un ou plusieurs contours d'une ou plusieurs lumières CLi caractéristiques et respectives d'un ou plusieurs composants Ci identifiés. A cet effet, conformément à un mode de réalisation non limitatif d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention décrit en lien avec la figure 10, une telle séquence peut également comprendre une étape de test 114, de sorte que lorsqu'il n'existe plus d'autre contour enceignant une lumière CLi caractéristique d'un composant Ci encore non identifié, situation illustrée par le lien 114n en figure 10, l'estimation d'une ou des quantités d'intérêts par itération, telle que représentée par le lien 114y entre l'étape 113 et l'étape 111, prend fin et les données relatives auxdites quantités d'intérêt sont prêtes à être exploitées par l'étape 120 du procédé 100 décrite à titre d'exemple non limitatif par la figure 10.
Afin de fournir une aide au diagnostic d'une pathologie pulmonaire à un expérimentateur, un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut comporter une étape 112, préalable à l'étape 113 pour inscrire dans la mémoire de données une structure de données associée à un composant Ci, pour caractériser un type CTi de composant Ci, à partir de la valeur d'une des quantités d'intérêt QI estimées lors de l'étape 111, présent dans une coupe histologique d'un poumon OG, tels qu'à titre d'exemples non limitatifs, des types de composants discriminant des alvéoles et/ou des sacs alvéolaires, des bronchioles, des bronches et des vaisseaux. En effet, lesdits types de composants tubulaires précédemment énoncés présentent des corps annulaires de morphologies très différentes, il peut être alors avantageux de caractériser le type de chaque composant Ci identifié, afin de pouvoir déterminer avec le plus de justesse ou de précision, si lesdits composants Ci présentent des modifications de leur morphologie, tel qu'il est souvent le cas dans le cadre de pathologies touchant les voies pulmonaires.
A titre d'exemple non limitatif, conformément à un mode de réalisation d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention, l'étape 112 pour caractériser un type CTi de composant Ci peut consister à comparer, en une sousétape 1121, la valeur du diamètre de Féret DFi d'un composant Ci, estimée lors de la mise en œuvre de l'étape 111, à un seuil de diamètre haut DFh et/ou un seuil de diamètre bas DFb, lesdits seuils DFh et DFb étant avantageusement et préalablement paramétrables. Ainsi, pour caractériser un composant Ci de type bronchiole, ladite sous-étape 1121 pourra consister à affecter une valeur prédéterminée, par exemple la valeur « 1 » à un champ dédié dans la structure de données associée audit composant Ci, si le diamètre de Féret DFi de ce dernier est supérieur à un seuil de diamètre bas DFb de cent micromètres et inférieur à un seuil de diamètre haut DFh de mille micromètres. En effet, les types CTi de composants Ci observés dans une coupe histologique d'un poumon OG présentent généralement un diamètre sensiblement compris entre dix micromètres et mille micromètres. L'utilisation du diamètre de Féret s'avère particulièrement avantageuse, puisque une telle utilisation facilite ainsi la classification des différents types CTi de composants Ci et permet à un expérimentateur de ne prendre en considération qu'un ou plusieurs types de composants d'intérêts parmi les différents types CTi de composants Ci identifiés, pour produire ultérieurement un indicateur graphique multiparamétrique I apportant une aide précieuse dans l'établissement d'un diagnostic d'une pathologie affectant les voies respiratoires.
En variante ou en complément, toujours en lien avec la figure 10, une étape 112 pour caractériser un type CTi d'un composant Ci d'un procédé d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut consister à comparer, en une sous-étape 1122, la valeur de l'épaisseur moyenne CWWi d'un composant Ci, estimée lors de la mise en œuvre de l'étape 111, à un seuil d'épaisseur moyenne haut CWWh et/ou un seuil d'épaisseur moyenne bas CWWb, lesdits seuils CWWh et CWWb étant avantageusement et préalablement paramétrables. Ainsi, pour caractériser un composant Ci de type bronchiole, ladite sous-étape 1122 pourra consister à affecter ladite valeur prédéterminée « 1 » dans un champ dédié dans la structure de données associée audit composant Ci, si l'épaisseur moyenne CWWi de ce dernier est supérieure à un seuil de diamètre haut CWWh de dix micromètres. En effet, les types CTI de composants Ci observés dans une coupe histologique d'un poumon OG présentent généralement une épaisseur moyenne CWWi supérieure à 10 micromètres. L'utilisation de l'épaisseur moyenne CWWi de la paroi CWi peut faciliter la classification des différents types CTI de composants Ci et permettre à un expérimentateur de ne prendre en considération que le ou les types CTI de composants Ci d'intérêt, pour produire ultérieurement, un indicateur graphique multiparamétrique I apportant une aide précieuse dans l'établissement d'un diagnostic d'une pathologie affectant les voies respiratoires.
De manière avantageuse, en variante ou en complément, selon la figure 10, une étape 112 pour caractériser un type CTI de composant Ci d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut consister à comparer, en une sous-étape 1123, la valeur de l'aire ICAi décrite par le contour intérieur CWICi d'un composant Ci, estimée lors de la mise en œuvre de l'étape 111, à un seuil d'aire haut ICAh et/ou un seuil d'aire bas ICAb, lesdits seuils ICAh et ICAb étant avantageusement et préalablement paramétrables. Ainsi, pour caractériser un composant Ci de type bronchiole, ladite sous-étape 1123 pourra consister à affecter la valeur prédéterminée « 1 » dans un champ dédié au type de composant dans la structure de données associée audit composant Ci, si l'aire ICAi de ce dernier est inférieur à un seuil d'aire bas ICAb de douze mille micromètres carrés et supérieur à un seuil d'aire haut ICAh de quinze mille micromètres carrés. En effet, les types CTi de composants Ci observés dans une coupe histologique d'un poumon OG, plus particulièrement les bronchioles, présentent généralement une aire ICAi sensiblement égale à treize mille micromètres carrés. L'utilisation de l'aire ICAi facilite ainsi la classification des différents types CTi de composants Ci et/ou la différenciation des différents stades d'atteinte tels que les stades aiqus et chroniques, dans le cadre d'une pathologie telle que le SAR, et permet à un expérimentateur de ne prendre en considération que le ou les types CTi de composants Ci jugés d'intérêt, pour produire ultérieurement un indicateur graphique multiparamétrique I apportant une aide précieuse dans l'établissement d'un diagnostic d'une pathologie affectant les voies respiratoires.
Trois exemples de caractérisation de types CTi de composants Ci viennent d'être précédemment décrits, en fonction notamment d'une quantité d'intérêt estimée. Toutefois, dans certains cas, l'emploi d'une seule quantité d'intérêt n'est pas suffisant pour discriminer certains types CTi de composants, par exemple les alvéoles et bronchioles. Pour pallier cet inconvénient, l'invention prévoit que l'étape 112 pour caractériser un type CTi de composant Ci d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention puisse être agencé pour exploiter conjointement différentes quantités d'intérêts.
Conformément à un mode de réalisation d'un tel procédé 100 décrit en liaison avec la figure 10, ladite étape 112 peut consister tout d'abord à comparer, en une sousétape 1124, la valeur du diamètre de Féret DFi d'un composant Ci estimée, lors de la mise en œuvre de l'étape 111 à un seuil de diamètre haut DFh et/ou un seuil de diamètre bas DFb, lesdits seuils DFh et DFb étant avantageusement et préalablement paramétrables. Si ladite valeur du diamètre de Féret DFi d'un composant Ci est comprise entre lesdits seuils DFb et DFh, ladite sous-étape 1124 peut ensuite consister à comparer la valeur de l'épaisseur moyenne CWWi dudit composant Ci, estimée lors de la mise en œuvre de l'étape 111, à un seuil d'épaisseur moyenne haut CWWh de dix micromètres, avantageusement à un seuil d'épaisseur moyenne haut CWWh de cinq micromètres. En effet, dans le cadre d'une atteinte des voies aériennes par une pathologie telle que le SAR, il est connu que les bronchioles subissent un remodelage, se caractérisant notamment par une augmentation de l'aire de la paroi bronchique et/ou de l'aire de la lumière CLi desdites bronchioles. De manière tout à fait ingénieuse, un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention prévoit que ladite sous-étape 1124 soit agencée, dans un premier temps, pour comparer la valeur du diamètre de Féret DFi d'un composant Ci à un seuil de diamètre bas DFb de cent micromètres et un seuil de diamètre haut DFh de mille micromètres. Afin de caractériser le type de différents composants Ci dont le diamètre de Féret DFi serait compris dans un tel intervalle, de cent micromètres à mille micromètres, en l'espèce les composants Ci de type bronchiole, alvéole et/ou sac alvéolaire, ladite sous-étape 1124 peut consister, dans un deuxième temps, à comparer la valeur CWWi dudit composant Ci, estimée lors de la mise en œuvre de l'étape 111, à un seuil d'épaisseur moyenne haut CWWh de cinq micromètres. Ainsi, si une telle épaisseur moyenne CWWi est supérieure audit seuil d'épaisseur CWWh, alors ladite sous-étape 1124 pourra consister à affecter ladite valeur prédéterminée « 1 » dans un champ dédié dans la structure de données associée audit composant Ci. Dans le cas contraire, si une telle épaisseur moyenne CWWi est inférieure audit seuil d'épaisseur CWWh, alors ladite sous-étape 1124 pourra consister à affecter ladite valeur prédéterminée « 2 » dans un champ dédié dans la structure de données associée audit composant Ci. Une telle valeur prédéterminée « 1 » pourra à titre d'exemple avantageux mais non limitatif, caractériser un composant Ci de type bronchiole. La mise en œuvre d'une telle sous-étape 1124, exploitant conjointement une pluralité de quantités d'intérêt, peut ainsi faciliter la classification des différents types CTI de composants Ci dont les morphologies sont proches, tels que, par exemple, les composants Ci de type bronchiole et de type alvéole et/ou sac alvéolaire. En effet, il est notamment avantageux de caractériser de telles bronchioles afin de pouvoir estimer le degré de sévérité de l'atteinte d'un sujet par une pathologie telle que le SAR, relatif au remodelage des composants Ci de type bronchiole.
Ainsi, selon un mode de réalisation préféré mais non limitatif d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention, l'étape 113 pour inscrire dans la mémoire de données une structure de données associée à chaque composant Ci identifié, peut consister en outre à inscrire dans un champ dédié de ladite structure de données une valeur prédéterminée pour caractériser, à partir d'une ou plusieurs quantités d'intérêt, un type CTi particulier d'un composant Ci, tel que par exemple le type bronchiole.
Par ailleurs, en complément, selon la figure 10, un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut comprendre une étape 115 pour filtrer les données préalablement estimées et inscrites dans la mémoire de données, permettant de ne tenir compte que des données relatives à un ou plusieurs types CTx déterminés parmi l'ensemble des valeurs de types CTi de composants Ci identifiés, afin de produire finalement et ultérieurement un indicateur graphique multiparamétrique I relatif auxdites données. A titre d'exemple non limitatif, une telle étape 115 peut consister à lire, dans la structure de données associée à un composant Ci, la valeur présente dans le champ caractérisant le type CTi dudit composant Ci. Pour ce faire, ladite étape 115 peut avantageusement être paramétrée de sorte que seules les données, plus particulièrement la ou les quantités d'intérêt préalablement estimées relatives aux composants Ci de type CTx « bronchiole », associées à un tel type de composant Ci comportant, dans le champ caractérisant le type CTi dudit composant Ci, une valeur prédéterminée CTx par exemple égale à « 1 », soient utilisées pour produire ultérieurement un indicateur graphique multiparamétrique I, en vue de fournir une aide au diagnostic d'une pathologie telle que à titre d'exemple non limitatif, le SAR.
En outre, afin de comparer les valeurs des quantités d'intérêt précédemment estimées à des valeurs de quantités d'intérêt étalon et finalement de faciliter la production d'un indicateur graphique multiparamétrique I, il peut être nécessaire de calculer ou estimer une moyenne des quantités d'intérêt. Ainsi, de manière avantageuse mais non limitative, un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut comprendre une étape 116 pour estimer une moyenne pour chaque ensemble de quantités d'intérêt QI par type CTi de composant Ci. L'étape 116 peut ainsi produire, pour un type CTi de composant Ci donné, éventuellement préalablement filtré lors de la mise en œuvre de l'étape 115, une ou plusieurs quantités d'intérêt QI moyennes, à partir des quantités d'intérêt QI estimées respectivement pour tous les composants Ci. On entend par « ensemble de quantité d'intérêt QI », toutes les valeurs précédemment estimées, par exemple relatives :
- au Diamètre de Féret DFi du contour extérieur CWOCi de la paroi CWi pour un type CTi de composant Ci donné ;
- à la distance moyenne CWWi séparant un contour extérieur CWOCi d'un contour intérieur CWICi pour un type CTi de composant Ci donné ;
- à l'aire ICAi décrite par un contour intérieur CWICi pour un type CTi de composant Ci donné ;
- à l'aire OCAi décrite par un contour extérieur CWOCi pour un type CTi de composant Ci donné,
- à l'aire CAi de la paroi CWi définie par la soustraction de l'aire ICAi décrite par un contour intérieur CWICi à l'aire OCAi décrite par un contour extérieur CWOCi pour un type CTI de composant Ci donné.
A titre d'exemples non limitatifs, pour produire un indicateur graphique multiparamétrique I permettant de caractériser une pathologie telle que le SAR, l'étape 116 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut consister à estimer, pour un ensemble de composants Ci de type bronchiole :
le diamètre de Féret moyen DFm correspondant au rapport entre la somme des diamètres de Féret DFi des composants Ci de type bronchiole et le nombre de composants Ci dont un diamètre de Féret DFi a été préalablement estimé ;
l'épaisseur moyenne CWWm correspondant au rapport entre la somme des épaisseurs moyennes CWWi des composants Ci de type bronchiole et le nombre desdits composants Ci dont ladite épaisseur CWWi a été préalablement estimée ; l'aire moyenne ICArn des lumières CLi correspondant au rapport entre la somme des aires ICAi estimées des composants Ci de type bronchiole et le nombre de composants Ci dont l'aire ICAi a été préalablement estimée ;
l'aire moyenne OCArn correspondant au rapport entre la somme des aires OCAi des composants Ci de type bronchiole et le nombre de composants dont l'aire OCAi a été préalablement estimée ;
l'aire moyenne CArn correspondant au rapport entre la somme des aires CAi des composants Ci de type bronchiole et le nombre de composants dont l'aire CAi a été préalablement estimée.
De manière avantageuse, un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut comporter une étape 120 consistant à produire une ou plusieurs représentations graphiques, destinées à transcrire une ou plusieurs quantités d'intérêt QI, telles que, à titre d'exemples non limitatifs, les quantités d'intérêt DFi, CWWi, ICAi, OCAi, CAi précédemment estimées, sous une forme graphique. Ledit procédé 100 peut alors comporter une étape 130 pour provoquer la restitution graphique conjointe desdites représentations graphiques préalablement produites à l'étape 120, et délivrer un indicateur graphique multiparamétrique I.
Préférentiellement mais non limitativement, une étape 120 d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut consister à produire une représentation graphique, pour un type CTI de composant Ci donné et/ou préalablement sélectionné au moyen d'une interface homme-machine d'entrée, d'une quantité d'intérêt QI relativement à une quantité d'intérêt étalon. On entend par « quantité d'intérêt étalon », une quantité d'intérêt de référence, éventuellement paramétrable, associée à chaque quantité d'intérêt QI estimée. Conformément aux exemples précédemment décrits, une telle quantité d'intérêt étalon peut consister non exhaustivement en :
le diamètre de Féret moyen DFe correspondant au rapport entre la somme des diamètres de Féret DFi des composants Ci de type bronchiole et le nombre de composants Ci dont un diamètre de Féret DFi a été préalablement estimé chez un sujet sain ;
l'épaisseur moyenne CWWe correspondant au rapport entre la somme des épaisseurs moyennes CWWi des composants Ci de type bronchiole et le nombre desdits composants Ci dont ladite épaisseur CWWi a été préalablement estimée chez un sujet sain ;
l'aire moyenne ICAe des lumières CLi correspondant au rapport entre la somme des aires ICAi estimées des composants Ci de type bronchiole et le nombre de composants Ci dont l'aire ICAi a été préalablement estimée chez un sujet sain ;
l'aire moyenne OCAe correspondant au rapport entre la somme des aires OCAi des composants Ci de type bronchiole et le nombre de composants dont l'aire OCAi a été préalablement estimée chez un sujet sain ;
l'aire moyenne CAe correspondant au rapport entre la somme des aires CAi des composants Ci de type bronchiole et le nombre de composants dont l'aire CAi a été préalablement estimée chez un sujet sain.
Un sujet sain est défini comme un sujet ne présentant aucune pathologie affectant les voies aériennes, telle que, à titre d'exemple non limitatif un SAR. Dans ce cas, les représentations graphiques des quantités d'intérêt QI d'un type CTi de composant Ci produites lors de l'étape 120 pour produire une représentation graphique et/ou l'étape 130 pour produire un indicateur graphique multiparamétrique I peuvent être exprimées en valeur relative ou normalisées, permettant ainsi auxdites quantités d'intérêt d'être exprimées selon une échelle commune et par voie de conséquence comparables de manière conjointe, de telles quantités d'intérêt étant préférentiellement mais non limitativement exprimées en pourcentage de quantités d'intérêt étalon correspondantes.
Préférentiellement, conformément à un premier exemple de représentation graphique décrit en lien avec la figure 4A, une telle étape 120 peut consister à produire une représentation graphique sous la forme d'une barre d'un diagramme en barres représentant, par exemple, une des quantités d'intérêt QI préalablement estimées, dont la hauteur ou la longueur, selon que la barre soit respectivement verticale ou horizontale, décrit la valeur relative de ladite quantité d'intérêt estimée au regard de ladite quantité d'intérêt étalon. Ladite représentation graphique peut en outre encoder une ou plusieurs données déterminées pour caractériser, lors de l'affichage de ladite représentation graphique via une interface homme-machine de sortie lors de l'étape 130, une texture, un motif ou un contour particulier, voire toute autre donnée permettant d'impacter le rendu graphique. En variante, de manière avantageuse, conformément à un deuxième exemple de représentation graphique décrit en lien avec la figure 4B, une étape 120 d'un procédé 100 conforme à l'invention peut consister à produire une représentation graphique sous la forme d'un point sur un axe d'un diagramme de Kiviat représentant, par exemple, une des quantités d'intérêt Ql préalablement estimées, dont la position sur ledit axe décrit la valeur relative de ladite quantité d'intérêt estimée au regard de ladite quantité d'intérêt étalon. Ainsi, conformément aux premier et deuxième exemples de représentations graphiques respectivement décrits en lien avec les figures 4A et 4B, les quantités d'intérêt moyennes DFm, CArn, CWWm, ICArn estimées pour un type CTI de composant Ci, d'une représentation numérique d'une coupe histologique d'un organe d'un sujet, peuvent être comparées conjointement et respectivement aux quantités d'intérêts étalon DFe, CAe, CWWe, ICAe correspondant à un sujet sain.
Selon un mode de réalisation préféré mais non limitatif d'un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention, la production des représentations graphiques à l'étape 120 dudit procédé 100 sera décrite pour un type de composant Ci d'intérêt sous la forme d'une bronchiole, en vue de produire un indicateur graphique multiparamétrique I. Un tel indicateur graphique multiparamétrique I est notamment décrit en lien avec les figures 4A, 4B, 5 et 6, après mise en œuvre de l'étape 130 pour provoquer la restitution graphique conjointe des représentations graphiques des quantités d'intérêt, telles que, par exemple, les quantités d'intérêt moyennes DFm, CWWm, ICArn, CArn, relativement aux quantités d'intérêt moyennes DFe, CWWe, ICAe, CAe chez un sujet sain. Avantageusement, lesdites quantités d'intérêt DFm, CWWm, ICArn et CArn peuvent être respectivement exprimées en pourcentage ou normalisées par rapport à une quantité d'intérêt étalon, la valeur numérique du pourcentage étant juxtaposée à la représentation numérique concernée. En variante ou en complément, une grille symbolisée, par exemple par des axes « 50% CTL » et « 100% CTL », qualifiés de « barres normées », tels que présentés en lien avec la figure 4A, ou encore par des lignes graduées de vingt à vingt, qualifiées de « axes normés », tels que présentées en lien avec la fiqure 4B, pourra être superposée aux représentations qraphiques lors de l'étape 130. Un expérimentateur peut donc instantanément constater que les quantités d'intérêt DFm, CWWm, ICArn et CArn, caractérisant la morphologie moyenne des composants de type bronchiole présents dans une représentation numérique d'une coupe histologique d'un poumon, sont différentes ou non des quantités d'intérêt étalon observées chez un sujet sain.
La figure 4A illustre un premier exemple non limitatif de représentation graphique sous la forme d'un diagramme en barres d'un indicateur graphique multiparamétrique produit par un procédé conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu. Selon ce premier exemple, conformément à la figure 4A, une telle étape 120 peut consister à exprimer une première quantité d'intérêt DFm, sous la forme d'une barre, par rapport à une quantité d'intérêt étalon DFe, symbolisée par l'axe « 100% CTL » de la grille. Ainsi, la représentation graphique DFm encode une valeur numérique de soixantedix-sept pour cent, signifiant que le diamètre de Féret moyen DFm chez le sujet examiné est inférieur de vingt trois pour cent au regard du diamètre de Féret moyen chez un sujet sain. De la même manière, et indépendamment de la représentation graphique de la première quantité d'intérêt DFm, l'étape 120 peut avantageusement consister à exprimer une deuxième quantité d'intérêt CWWm, sous la forme d'une barre, par rapport à une quantité d'intérêt étalon CWWe, symbolisée par l'axe « 100% CTL » de la grille. Ainsi, la représentation graphique CWWm encode une valeur numérique de soixante-dix-neuf pour cent, signifiant que l'épaisseur moyenne CWWm chez le sujet examiné est inférieure de vingt-un pour cent au regard de l'épaisseur moyenne chez un sujet sain. De la même manière, et indépendamment des représentations graphiques respectives des première et deuxième quantités d'intérêt DFm et CWWm, l'étape 120 peut avantageusement consister à exprimer une troisième quantité d'intérêt ICArn, sous la forme d'une barre, par rapport à une quantité d'intérêt étalon ICAe, symbolisée par l'axe « 100% CTL » de la grille. Ainsi, la représentation graphique ICArn encode une valeur numérique de quatre-vingt-sept pour cent, signifiant que l'aire moyenne ICArn des lumières chez le sujet examiné est inférieure de treize pour cent au regard de l'aire moyenne des lumières chez un sujet sain. De la même manière, et indépendamment des représentations graphiques respectives des première, deuxième et troisième quantités d'intérêt DFm, CWWm et ICArn, l'étape 120 peut avantaqeusement consister à exprimer une quatrième quantité d'intérêt CArn, sous la forme d'une barre, par rapport à une quantité d'intérêt étalon CAe, symbolisée par l'axe « 100% CTL » de la grille. Ainsi, la représentation graphique CArn encode une valeur numérique de soixante-dix pour cent, signifiant que l'aire moyenne CAm des parois chez le sujet examiné est inférieure de vinqt-et-un pour cent au regard de l'aire moyenne des parois chez un sujet sain. Il sera ainsi possible, lors de la mise en œuvre de l'étape 130, de provoquer la restitution graphique conjointe desdites représentations qraphiques ainsi produites via une interface homme-machine de sortie d'un système mettant en œuvre un procédé 100 selon l'invention, sous la forme, par exemple, de diagrammes en barres, afin de fournir à un expérimentateur un indicateur graphique multiparamétrique I, présentant côte à côte les différentes barres.
Un tel premier exemple de restitution graphique conjointe, sous la forme d'un diagramme en barres, est présenté en figure 4A pour un sujet respectivement atteint d'un SAR dit aigu. Ainsi, ladite figure 4A exprime conjointement les représentations graphiques respectives des première, deuxième, troisième et quatrième quantités d'intérêt DFm, CWWm, ICAm et CAm sous la forme d'une représentation graphique de type diagramme en barres, et expriment les valeurs respectives desdites quantités d'intérêt DFm, CWWm, ICAm et CAm relativement à celles des mêmes quantités d'intérêt DFe, CWWe, ICAe et CAe chez un sujet sain. Un expérimentateur peut ainsi prendre aisément en considération par la simple visualisation de l'indicateur graphique multiparamétrique I que les quantités d'intérêt moyennes DFm, CWWm, ICAm et CAm chez un sujet examiné sont plus faibles que celles DFe, CWWe, ICAe et CAe chez un sujet sain. Un tel indicateur graphique multiparamétrique I est alors adapté pour fournir une aide précieuse à un expérimentateur ou à un personnel de santé, afin d'éventuellement diagnostiquer une atteinte des voies pulmonaires par une pathologie, telle que le SAR par exemple. En effet, un indicateur graphique multiparamétrique I, tel que présenté en lien avec la figure 4A, fournit des informations pertinentes concernant la morphologie générale de l'ensemble des composants de type bronchiole présents dans une coupe histologique d'un poumon d'un sujet atteint d'un SAR aigu. Une telle atteinte induit généralement une desquamation plus ou moins hétérogène de l'épithélium bronchique, c'est-à-dire une diminution de l'épaisseur moyenne CWWi de la paroi CWi des bronchioles. Cette desquamation de la paroi bronchique entraîne, par voie de conséquence, une diminution de l'aire totale CAi de la bronchiole, éventuellement de l'aire ICAi de la lumière CLi et du diamètre de Féret DFi.
En variante ou en complément, pour renforcer la discrimination visuelle d'un indicateur graphique multiparamétrique I relatif à un remodelage de l'épithélium bronchique humain ou animal, les étapes 120 et/ou 130 d'un procédé pour produire un tel indicateur graphique multiparamétrique I peuvent, à titre d'exemples non limitatifs, consister à exprimer puis afficher des quantités d'intérêt moyennées, sous la forme d'une représentation graphique de type diagramme de Kiviat, également connu sous la dénomination « radar », par rapport ou relativement à ces mêmes quantités d'intérêt étalon, lesdites quantités d'intérêt étalon étant indiquées comme correspondant à 100% par rapport à un axe gradué. Un tel axe peut être gradué, à titre d'exemple non limitatif, à partir de zéro pour cent de la valeur de la quantité étalon correspondante chez un sujet sain. Ainsi, chaque quantité d'intérêt QI peut être normalisée au regard de la quantité d'intérêt étalon correspondante et peuvent être visualisables après restitution graphique lors de l'étape 130 dans une échelle continue. La figure 4B illustre un deuxième exemple non limitatif de représentation graphique sous la forme d'un diagramme de Kiviat d'un indicateur graphique multiparamétrique I produit par un procédé 100 conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu.
Selon ce deuxième exemple, conformément à la figure 4B, une telle étape 120 peut consister à exprimer une première quantité d'intérêt DFm, sous la forme d'une position sur un axe gradué, par rapport à une quantité d'intérêt étalon DFe, symbolisée par la position 100 sur l'axe gradué. De la même manière, et indépendamment de la représentation graphique de la première quantité d'intérêt DFm, l'étape 120 peut avantageusement consister à exprimer une deuxième quantité d'intérêt CWWm, sous la forme d'une position sur un axe gradué, par rapport à une quantité d'intérêt étalon CWWe, symbolisée par la position 100 sur l'axe gradué. De la même manière, et indépendamment des représentations graphiques respectives des première et deuxième quantités d'intérêt DFm et CWWm, l'étape 120 peut avantageusement consister à exprimer une troisième quantité d'intérêt ICArn, sous la forme d'une position sur un axe gradué, par rapport à une quantité d'intérêt étalon ICAe, symbolisée par la position 100 sur l'axe gradué. De la même manière, et indépendamment des représentations graphiques respectives des première, deuxième et troisième quantités d'intérêt DFm, CWWm et ICAm, l'étape 120 peut avantageusement consister à exprimer une quatrième quantité d'intérêt CArn, sous la forme d'une position sur un axe gradué, par rapport à une quantité d'intérêt étalon CAe, symbolisée par la position 100 sur l'axe gradué.
Il sera ainsi possible, lors de l'étape 130 de provoquer la restitution graphique conjointe desdites représentations qraphiques ainsi produites via une interface homme-machine de sortie d'un système mettant en œuvre un procédé 100 selon l'invention, tel que, à titre d'exemple non limitatif, sous la forme, par exemple d'un diagramme de Kiviat, afin de fournir à un expérimentateur un indicateur graphique multiparamétrique I. Une telle restitution graphique conjointe peut en outre illustrer la matérialisation d'une aire définie par les différentes positions respectives des quantités d'intérêt normalisées sur les axes normés. Un tel premier exemple de restitution graphique conjointe, sous la forme d'un diagramme de Kiviat, est présenté en figure 4B pour un sujet atteint d'un SAR dit aigu. Ainsi, ladite figure 4B exprime conjointement les représentations graphiques des première, deuxième, troisième et quatrième quantités d'intérêt DFm, CWWm, ICAm et CArn sous la forme d'une représentation graphique de type diagramme de Kiviat. L'aire, décrivant les positions conjointes des quantités d'intérêt DFm, CWWm, ICAm et CArn, délimitée par un contour en pointillés, peut ainsi indiquer aisément l'état pathologique de la morphologie des composants Ci d'intérêt représentées sur chaque axe du diagramme au regard de l'aire, indiquée par le contour en trait plein, décrivant une morphologie type de tels composants chez un sujet sain. Il est ainsi plus aisé, pour l'expérimentateur, de percevoir la distance ou l'écart entre l'état pathologique de la morphologie, décrit par l'aire délimitée en pointillés, d'un sujet examiné au regard d'une morphologie chez un sujet sain. En l'espèce, selon l'exemple décrit par la figure 4B, les composants Ci d'intérêt de type bronchiole, présents dans une représentation numérique d'une coupe histologique d'un poumon, présentent des quantités d'intérêt qui ont diminué par rapport aux mêmes quantités d'intérêt étalon.
La figure 5 illustre un troisième exemple non limitatif de représentation graphique sous la forme d'un diagramme en barres d'un indicateur graphique multiparamétrique produit par un procédé conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu. A l'instar des figures 4A et 4B, la figure 5 exprime conjointement des représentations graphiques respectives des première, deuxième, troisième et quatrième quantités d'intérêt DFm, CWWm, ICArn et CArn sous la forme d'une représentation graphique de type diagramme en barres. Toutefois, la figure 5 présente des représentations graphiques respectives des première, deuxième, troisième et quatrième quantités d'intérêt DFm, CWWm, ICArn et Cam sous la forme d'écarts relatifs signés, de tels écarts consistant en des valeurs numériques positives ou négatives calculées à partir de la formule mathématique Æm = gim-gfe suivante pour chaque quantité d'intérêt QI:
, où Em consiste en la valeur de l'écart moyen estimé, consiste en la valeur de la quantité d'intérêt moyenne estimée chez un sujet examiné, consiste en la valeur de la quantité d'intérêt étalon chez un sujet sain.
Selon ce troisième exemple, conformément à la figure 5, une telle étape 120 peut consister à exprimer des écarts relatifs signés entre des première, deuxième, troisième et quatrième quantités d'intérêt DFm, CWWm, ICAm et CArn, respectivement sous la forme de barres, par rapport à un horizon vertical ou horizontal nul décrit par un axe « 100% CTL ». De tels écarts peuvent ainsi représentés une diminution ou une augmentation relative ou normalisée desdites des première, deuxième, troisième et quatrième quantités d'intérêt DFm, CWWm, ICAm et CArn, par rapport à des quantités d'intérêt respectives étalons. En l'espèce, pour un horizon vertical, une diminution de telle quantité d'intérêt est affichée sur la gauche tandis qu'une augmentation de telle quantité d'intérêt est affichée sur la droite.
Il sera ainsi possible, lors de la mise en œuvre de l'étape 130, de provoquer la restitution graphique conjointe desdites représentations graphiques ainsi produites via une interface homme-machine de sortie d'un système mettant en œuvre un procédé 100 selon l'invention, sous la forme, par exemple, de diagrammes en barres, afin de fournir à un expérimentateur un indicateur graphique multiparamétrique I, présentant côte à côte les différentes barres.
La figure 6 illustre un quatrième exemple non limitatif de représentation graphique sous la forme d'un diagramme de Kiviat d'un indicateur graphique multiparamétrique I produit par un procédé 100 conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR aigu. A l'instar des figures 4A, 4B et 5, la figure 6 exprime conjointement des représentations graphiques respectives des première, deuxième, troisième et quatrième quantités d'intérêt DFm, CWWm, ICArn et CArn sous la forme d'une représentation graphique de type diagramme de Kiviat. Toutefois, la figure 6 présente des représentations graphiques respectives des première, deuxième, troisième et quatrième quantités d'intérêt DFm, CWWm, ICArn et Cam sous la forme d'écarts absolus signés, de tels écarts consistant en des valeurs numériques positives ou négatives calculées à partir de la formule mathématique suivante pour chaque quantité d'intérêt QI:
Em < 0 —> Em = 0,1
Em = Q Im- Qle —> <
Em ~ 0 —> Em = 0,5 où Em consiste en la valeur de l'écart moyen estimé, consiste en la valeur de la quantité d'intérêt moyenne estimée chez un sujet examiné, consiste en la valeur de la quantité d'intérêt étalon chez un sujet sain. Un tel quatrième exemple de représentation graphique, sous la forme d'un diagramme de Kiviat, s'avère particulièrement avantageux, puisqu'une telle représentation graphique permet à un expérimentateur d'aisément visualiser graphiquement tout différence significative entre des quantités d'intérêt QI estimées et des quantités d'intérêt QI étalon, quand bien même celles-ci sont relativement proches.
Selon ce quatrième exemple, conformément à la figure 6, une telle étape 120 peut consister à exprimer des écarts absolus déterminés entre des première, deuxième, troisième et quatrième quantités d'intérêt moyennes DFm, CWWm, ICAm et CArn et des première, deuxième, troisième et quatrième quantités d'intérêt étalon respectives et correspondantes, respectivement sous la forme de positions sur des axes gradués. De tels écarts peuvent ainsi représentés une diminution ou une augmentation absolue et normalisée, en l'espèce selon trois valeurs croissantes prédéterminées « 0,1 », « 0,5 » et « 1 », desdites des première, deuxième, troisième et quatrième quantités d'intérêt DFm, CWWm, ICAm et CArn, par rapport à des quantités d'intérêt respectives étalons. L'invention ne saurait être limitée aux valeurs prédéterminées des écarts Em. Avantageusement, les écarts ou lesdites quantités d'intérêt DFm, CWWm, ICAm et CArn peuvent être respectivement exprimés par un indice caractérisant une diminution, une équivalence ou une augmentation par rapport à une quantité d'intérêt étalon, en l'espèce une telle équivalence à une quantité d'intérêt étalon est représentée sur chaque axe du diagramme, par la graduation « 0,5 » correspondante. Aussi, l'invention prévoit un seuil de tolérance paramétrable ou prédéterminé, par exemple plus ou moins dix pour cent autour de la valeur étalon pour déterminer ainsi une telle équivalence entre la quantité d'intérêt estimée et la quantité étalon correspondante. De la même manière, l'invention prévoit qu'en deçà d'un tel seuil, ladite quantité d'intérêt soit considérée comme inférieure, valeur « 0,1 », ou supérieur valeur « 1 », à ladite quantité d'intérêt étalon. Afin de visualiser graphiquement une différence jugée significative par l'expérimentateur, un tel quatrième exemple de représentation graphique, sous la forme d'un diagramme de Kiviat, peut avantageusement être agencée pour afficher une augmentation ou une diminution pour chaque quantité d'intérêt QI estimée par rapport à une quantité d'intérêt étalon si et seulement si ladite quantité d'intérêt QI est supérieure ou inférieure à un seuil paramétrable.
En l'espèce, conformément à la figure 6, nous pouvons constater une diminution généralisée des première, deuxième, troisième et quatrième quantités d'intérêt moyennes DFm, CWWm, ICArn et CArn, lesdits écarts en lien avec lesdites première, deuxième, troisième et quatrième quantités d'intérêt moyennes étant toutes égales à « 0,1 ». Il sera ainsi possible, lors de l'étape 130 de provoquer la restitution graphique conjointe desdites représentations qraphiques ainsi produites via une interface homme-machine de sortie d'un système mettant en œuvre un procédé 100 selon l'invention, tel que, à titre d'exemple non limitatif, sous la forme, par exemple d'un diagramme de Kiviat, afin de fournir à un expérimentateur un indicateur graphique multiparamétrique I. Le diagramme de Kiviat présenté, en lien avec la figure 6, décrit ainsi directement une diminution des quantités d'intérêt DFm, CWWm, ICArn et CArn une diminution d'une telle quantité d'intérêt est affichée sur la graduation « 0,1 » de l'axe correspondant, tandis qu'une augmentation d'une telle quantité d'intérêt est affichée sur la graduation « 1 » de l'axe correspondant.
Une telle restitution graphique conjointe peut en outre illustrer la matérialisation d'une aire définie par les différentes positions respectives des quantités d'intérêt normalisées sur les axes normés. Un tel premier exemple de restitution graphique conjointe, sous la forme d'un diagramme de Kiviat, est présenté en figure 4B pour un sujet atteint d'un SAR dit aigu. A l'instar de la figure 4B, l'aire, décrivant les positions conjointes des quantités d'intérêt DFm, CWWm, ICArn et CArn, délimitée par un contour en pointillés, peut ainsi indiquer aisément l'état pathologique de la morphologie des composants Ci d'intérêt représentées sur chaque axe du diagramme au regard de l'aire, indiquée par le contour en trait plein, décrivant une morphologie type de tels composants chez un sujet sain. Conformément à la figure 6, l'aire décrivant l'état pathologique est réduite ou nulle au regard de l'air décrivant la morphologie type.
A l'instar des figures 4A, 4B, 5 et 6, les figures 7A, 7B, 8 et 9 illustrent respectivement des premier, deuxième, troisième et quatrième exemples non limitatifs de représentations graphiques sous la forme de diagrammes en barres et de diagrammes de Kiviat, d'un indicateur graphique multiparamétrique I produit par un procédé 100 conforme à l'invention, ledit indicateur décrivant une pluralité de quantités d'intérêt estimées relatives à une pluralité de quantités d'intérêt étalons respectives sur des axes normés, chez un sujet atteint d'un SAR chronique. Un tel indicateur graphique multiparamétrique I est ainsi particulièrement adapté pour décrire la morphologie des bronchioles. En effet, il est connu qu'une pathologie, telle que le SAR, peut induire à long terme une inflammation chronique des bronchioles, provoquant un remodelage de ces dernières. A l'instar des figures 4A, 4B, 5 et 6, les figures 7A, 7B, 8 et 9 expriment conjointement des représentations graphiques respectives des première, deuxième, troisième et quatrième quantités d'intérêt DFm, CWWm, ICArn et CArn sous la forme de représentations graphiques de type diagramme en barres et diagrammes de Kiviat.
Un même procédé 100 pour produire un indicateur graphique multiparamétrique I est appliqué à une deuxième coupe histologique d'un poumon d'un sujet atteint cette fois-ci d'un SAR chronique. Une telle pathologie est aisément diagnosticable par un expérimentateur à la vue d'un ou plusieurs indicateurs graphiques multiparamétriques I. A cet effet, un indicateur graphique multiparamétrique I, tel que présenté en figure 7A, fournit des informations pertinentes concernant la morphologie générale de l'ensemble des composants de type bronchiole présents dans une coupe histologique d'un poumon d'un sujet atteint d'un SAR chronique. Une telle atteinte induit généralement une inflammation plus ou moins importante des bronchioles. La figure 7A présente ainsi une augmentation globale des quantités d'intérêt DFm, ICArn, caractérisant une inflammation desdites bronchioles ainsi qu'une « stabilisation » relative des quantités d'intérêt CAm, CWm, caractérisant une régénération de l'épithélium bronchique, chacune desdites quantités d'intérêt étant exprimée en pourcentage de la quantité d'intérêt étalon correspondante, c'est-à-dire en l'espèce chez un sujet sain.
Par ailleurs chez un sujet atteint d'un SAR chronique, en lien avec la figure 7B, l'expérimentateur peut aisément constater que les mêmes quantités d'intérêts DFm, ICAm et CAm sont manifestement plus élevées que chez un sujet sain. L'utilisation d'un tel diagramme de Kiviat permet, en outre, d'opposer directement des quantités d'intérêt QI censées être corrélées entre elles. La production d'un indicateur graphique multiparamétrique I sous la forme d'un diagramme de Kiviat permet de visuellement mettre en évidence une augmentation de la surface interne des lumières CLi et plus généralement de la taille des bronchioles respectivement décrites, en lien avec la figure 7B, par l'aire ICAm et le diamètre de Féret DFm.
Les figures 8 et 9 mettent d'autant plus en exergue des quantités d'intérêt moyennes CWWm et CAm sensiblement identiques aux quantités d'intérêt étalon, tandis que les quantités d'intérêt moyennes DFm et ICAm sont en nette augmentation, de telles quantités d'intérêt moyennes étant caractéristiques d'un sujet atteint d'un SAR chronique.
Par ailleurs, en variante ou en complément, un procédé 100 pour produire un indicateur graphique multiparamétrique I conforme à l'invention peut comporter une étape 132 pour provoquer la restitution ou sortie graphique d'un tel indicateur graphique multiparamétrique I conforme à l'invention par toute interface homme-machine adaptée, telle que, par exemple, non limitativement un écran d'ordinateur, coopérant avec l'objet électronique mettant en œuvre ledit procédé 100. En variante ou en complément, une telle restitution ou sortie peut être écrite ou sonore, respectivement via un périphérique de sortie de type imprimante ou encore un périphérique de sortie sonore via un haut-parleur.
En complément de l'étape 132, un procédé 100 pour produire un indicateur graphique multiparamétrique I selon l'invention peut comporter une ou plusieurs étapes 131, 133, 134, pour provoquer une restitution graphique respectivement de représentations numériques de type RDI, MRI et des quantités d'intérêt QI estimées sous la forme de cartes. De telles cartes peuvent être restituées graphiquement au moyen d'un périphérique de sortie identique ou distinct de celui délivrant l'indicateur graphique multiparamétrique précédemment élaboré, quantité d'intérêt par quantité d'intérêt à l'étape 120.
De cette manière, l'utilisateur d'un objet électronique adapté pour mettre en œuvre un procédé conforme à l'invention tel que le procédé 100, dispose d'une pluralité d'informations objectives, reproductibles et instantanées aidant à diagnostiquer une pathologie telle qu'un SAR. L'ensemble des étapes 131 à 134 constituent ainsi un traitement 130, destiné à restituer à l'utilisateur un indicateur graphique multiparamétrique I, une ou plusieurs quantités d'intérêt QI estimées, voire une ou plusieurs cartes, en l'espèce, une ou plusieurs représentations numériques parmi les représentations numériques RDI, MRI évoquées précédemment.
Alternativement mais non limitativement, l'invention prévoit que les quantités d'intérêt DFi, CAi, CWWi, ICAi estimées, ou plus généralement une quantité d'intérêt QI, pour chaque type CTI de composant Ci, d'une représentation numérique d'une coupe histologique d'un organe d'un sujet examiné, puissent être comparées non plus aux quantités d'intérêt étalon correspondantes chez un sujet sain, mais aux quantités d'intérêts moyennes DFm, CArn, CWWm, ICArn estimées chez ce même sujet. Selon cette variante, l'invention prévoit qu'une ou plusieurs étapes de restitution graphique 133 ou 134 d'un procédé 100 conforme à l'invention puissent consister en outre à affecter, aux pixels associés à un composant Ci d'intérêt, une couleur et/ou une valeur déterminée exprimant un accroissement ou une diminution significative au regard de la morphologie moyenne des autres composants.
L'invention a été notamment décrite en lien avec l'analyse d'un lobe pulmonaire d'une souris. Toutefois, elle ne saurait être limitée à ce seul exemple de réalisation et/ou d'application. D'autres modifications peuvent être envisagées sans sortir du cadre de la présente invention pour adapter, en variante ou en complément, le procédé pour produire un indicateur graphique multiparamétrique chez l'humain ou un autre animal, voire à un organe présentant des similarités anatomiques avec le poumon.
Par ailleurs, l'invention ne saurait se limiter aux seules représentations graphiques sous la forme de diagrammes en barres et/ou de diagrammes de Kiviat précédemment décrites, dont les quantités d'intérêt sont 5 notamment exprimées en fonction de quantité d'intérêt étalon. En variante, d'autres représentations adaptées pour produire un indicateur graphique multiparamétrique I auraient pu être employées.

Claims (15)

  1. REVENDICATIONS
    1. Procédé (100) pour produire un indicateur graphique multiparamétrique (I) relatif à un remodelage de l'épithélium bronchique humain ou animal à partir d'une représentation numérique (RDI, MRI) d'une coupe histologique d'un poumon (OG) comportant un ou plusieurs composants (Ci) de forme annulaire dont chacun décrit une paroi (CWi) enceignant une lumière (CLi), ledit procédé (100) étant mis en œuvre par une unité de traitement d'un système d'aide au diagnostic d'une pathologie, ledit système comportant en outre une interface hommemachine de sortie et une mémoire de données, caractérisé en ce que ledit procédé (100) comporte :
    - une étape (111) pour estimer des quantités d'intérêt (QI) relatives aux morphologies respectives d'un ou plusieurs composants (Ci) présents dans la coupe histologique, lesdites quantités d'intérêt appartenant à un ensemble de quantités d'intérêt comprenant :
    i. le diamètre de Féret (DFi) du contour extérieur (CWOCi) de la paroi d'un composant (Ci) ; ii. la distance moyenne (CWW i) séparant ledit contour extérieur (CWOCi) du contour intérieur (CWICi), traduisant une épaisseur moyenne de la paroi (CWi) d'un composant (Ci) ;
    iii. l'aire (ICAi) décrite par ledit contour intérieur (CWICi) d'un composant (Ci) traduisant la lumière (CLi) de ce dernier;
    iv. l'aire (OCAi) décrite par ledit contour extérieur (CWOCi) d'un composant (Ci) traduisant l'aire totale couverte par ce dernier ;
    v. l'aire (CAi) de la paroi (CWi) d'un composant (Ci) définie par la soustraction de l'aire (ICAi) décrite par ledit contour intérieur (CWICi) à l'aire (OCAi) décrite par ledit contour extérieur (CWOCi) dudit composant (Ci) ;
    - une étape (120) pour produire, par quantité d'intérêt estimée, une représentation graphique de celle-ci relative à une quantité d'intérêt étalon ;
    - une étape (130) pour provoquer la restitution graphique conjointe des représentations graphiques desdites quantités d'intérêt relatives aux quantités d'intérêt étalon respectives préalablement produites, par l'interface homme-machine de sortie du système.
  2. 2. Procédé (100) selon pour lequel l'étape quantités d'intérêts respectives d'un ou comporte une étape mémoire de données
    La revendication précédente, (111) pour estimer des relatives aux morphologies plusieurs composants (Ci) 113) pour inscrire dans la une structure de données associée à chaque composant (Ci), ladite structure de données comportant un champ pour enregistrer la valeur de chaque quantité d'intérêt estimée.
  3. 3. Procédé (100) selon la revendication précédente, comportant une étape (112) pour caractériser un type (CTi) de composant (Ci) à partir de la valeur d'une des quantités d'intérêt (QI) estimées, l'étape (113) pour inscrire dans la mémoire de données une structure de données associée à chaque quantité d'intérêt estimée d'un composant (Ci) consiste à inscrire dans un champ de ladite structure de données une valeur caractérisant un type (CTi) de composant (Ci) déterminé.
  4. 4. Procédé (100) selon la revendication précédente, pour lequel, lorsque l'une des quantités d'intérêt estimées consiste en le diamètre de Féret (DFi) du contour extérieur (CWOCi) de la paroi du composant (Ci), l'étape (112) pour caractériser un type (CTi) de composant (Ci) comporte une opération (1121) de comparaison de la valeur dudit diamètre de Féret (DFi) estimée à un seuil de diamètre haut (DFh) et/ou un seuil de diamètre bas (DFb).
  5. 5. Procédé (100) selon la revendication 3, pour lequel, lorsque l'une des quantités d'intérêt estimées consiste en l'épaisseur moyenne (CWWi) de la paroi (CWi) du composant (Ci), l'étape (112) pour caractériser un type de composant (Ci) comporte une opération (1122) de comparaison de la valeur de ladite épaisseur moyenne à un seuil d'épaisseur haut (CWWh) et/ou un seuil d'épaisseur bas (CWWb) prédéterminés.
  6. 6. Procédé (100) selon la revendication 3, pour lequel, lorsque l'une des quantités d'intérêt estimées consiste en l'aire (ICAi) décrite par le contour intérieur (CWICi) du composant (Ci),
    1'étape (112) pour caractériser un type de composant (Ci) comporte une opération (1123) de comparaison de la valeur de ladite aire à un seuil d'aire haut et/ou un seuil d' aire bas
    prédéterminés.
  7. 7. Procédé (100) selon l'une quelconque des revendications 3 à 6, pour lequel l'étape (120) pour produire, par quantité d'intérêt estimée, une représentation graphique de celle-ci relative à une quantité d'intérêt étalon et/ou l'étape (130) pour provoquer la restitution graphique conjointe des représentations graphiques desdites quantités d'intérêt relatives aux quantités d'intérêt étalon ne sont mises en œuvre (115) que pour un type de composant (Ci) caractérisé déterminé (CTx).
  8. 8. Procédé (100) selon la revendication précédente, pour lequel le type (CTi) de composant (Ci) caractérisé déterminé est une bronchiole (CTx).
  9. 9. Procédé (100) selon l'une quelconque des revendications précédentes, pour lequel l'étape (130) pour provoquer la restitution graphique conjointe des représentations graphiques desdites quantités d'intérêt relatives aux quantités d'intérêt étalon respectives préalablement produites par l'interface homme-machine de sortie du système, consiste en l'affichage par cette dernière d'un diagramme de Kiviat, décrivant au moins trois représentations graphiques de quantités d'intérêt relatives aux quantités d'intérêt étalon respectives sur des axes normés.
  10. 10. Procédé (100) selon l'une quelconque des revendications 1 à 9, pour lequel l'étape (130) pour provoquer la restitution graphique conjointe des représentations graphiques desdites quantités d'intérêt relatives aux quantités d'intérêt étalon respectives préalablement produites, par l'interface homme-machine de sortie du système, consiste en l'affichage par cette dernière d'un diagramme en barres décrivant les représentations graphiques de quantités d'intérêt relatives aux quantités d'intérêt étalon respectives par des barres normées.
  11. 11. Procédé (100) selon la revendication 9 ou la revendication 10, pour lequel la normalisation d'un axe ou d'une barre consiste à exprimer la valeur d'une quantité d'intérêt estimée en pourcentage de la valeur de la quantité d'intérêt étalon associée.
  12. 12. Procédé (100) selon la revendication 9 ou la revendication 10, pour lequel la normalisation d'un axe ou d'une barre consiste à exprimer la valeur d'une quantité d'intérêt estimée relativement à la valeur de la quantité d'intérêt étalon associée sous la forme de trois valeurs prédéterminées décrivant respectivement des valeurs de quantité d'intérêt estimées sensiblement inférieures, proches ou supérieures aux valeurs des quantités d'intérêt étalon associées.
  13. 13. Objet électronique d'un système d'aide au diagnostic d'une pathologie, ledit objet électronique comportant une unité de traitement et coopérant avec une interface homme-machine de sortie et avec une mémoire de données et caractérisée en ce que ladite mémoire de données comporte :
    une représentation numérique (RDI, MRI) d'une coupe histologique d'un organe (OG) humain ou animal ;
    des instructions exécutables ou interprétables par l'unité de traitement dont l'interprétation ou l'exécution desdites instructions par ladite unité de traitement provoque la mise en œuvre d'un procédé (100) selon l'une quelconque des revendications 1 à 12 .
  14. 14. Système d'aide au diagnostic d'une pathologie comportant un objet électronique selon la revendication précédente et une interface hommemachine de sortie apte à restituer à un utilisateur un indicateur graphique multiparamétrique (I) selon un procédé (100) conforme à l'une quelconque des revendications 1 à 12 et mis en œuvre par ledit objet électronique.
  15. 15. Produit programme d'ordinateur comportant une ou plusieurs instructions interprétables ou exécutables par une unité de traitement d'un objet électronique conforme à la revendications 13, caractérisé en ce que l'interprétation ou l'exécution desdites instructions par ladite unité de traitement provoque la mise en œuvre d'un procédé (100) selon l'une quelconque des revendications 1 à 12.
FR1855953A 2018-06-29 2018-06-29 Systeme et procede pour produire un indicateur graphique multiparametrique a partir d'une image d'une coupe histologique Ceased FR3083358A1 (fr)

Priority Applications (4)

Application Number Priority Date Filing Date Title
FR1855953A FR3083358A1 (fr) 2018-06-29 2018-06-29 Systeme et procede pour produire un indicateur graphique multiparametrique a partir d'une image d'une coupe histologique
US17/256,316 US20210264596A1 (en) 2018-06-29 2019-06-28 System and method for producing a multiparameter graphic indicator from an image of a histological section
PCT/FR2019/051614 WO2020002859A1 (fr) 2018-06-29 2019-06-28 Système et procédé pour produire un indicateur graphique multiparamétrique à partir d'une image d'une coupe histologique
EP19752216.2A EP3815037A1 (fr) 2018-06-29 2019-06-28 Système et procédé pour produire un indicateur graphique multiparamétrique à partir d'une image d'une coupe histologique

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1855953A FR3083358A1 (fr) 2018-06-29 2018-06-29 Systeme et procede pour produire un indicateur graphique multiparametrique a partir d'une image d'une coupe histologique
FR1855953 2018-06-29

Publications (1)

Publication Number Publication Date
FR3083358A1 true FR3083358A1 (fr) 2020-01-03

Family

ID=65685421

Family Applications (1)

Application Number Title Priority Date Filing Date
FR1855953A Ceased FR3083358A1 (fr) 2018-06-29 2018-06-29 Systeme et procede pour produire un indicateur graphique multiparametrique a partir d'une image d'une coupe histologique

Country Status (4)

Country Link
US (1) US20210264596A1 (fr)
EP (1) EP3815037A1 (fr)
FR (1) FR3083358A1 (fr)
WO (1) WO2020002859A1 (fr)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3656308A1 (fr) * 2018-11-20 2020-05-27 Koninklijke Philips N.V. Appareil de visualisation à base d'alvéoles en imagerie par rayons x à champ sombre ou à contraste de phase

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140064582A1 (en) * 2011-02-08 2014-03-06 Aalborg Universitet System for determining flow properties of a blood vessel
US20170046839A1 (en) * 2015-08-14 2017-02-16 Elucid Bioimaging Inc. Systems and methods for analyzing pathologies utilizing quantitative imaging

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140064582A1 (en) * 2011-02-08 2014-03-06 Aalborg Universitet System for determining flow properties of a blood vessel
US20170046839A1 (en) * 2015-08-14 2017-02-16 Elucid Bioimaging Inc. Systems and methods for analyzing pathologies utilizing quantitative imaging

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JEAN SERRA: "dilatation morphologique du contour interne", IMAGE ANALYSIS AND MATHEMATICAL MORPHOLOGY, 1982
SAARY ET AL: "Radar plots: a useful way for presenting multivariate health care data", JOURNAL OF CLINICAL EPIDEMIOLOGY, PERGAMON, GB, vol. 61, no. 4, 6 March 2008 (2008-03-06), pages 311 - 317, XP022503126, ISSN: 0895-4356, DOI: 10.1016/J.JCLINEPI.2007.04.021 *
SATOCHI SUZUKI ET AL., TOPOLOGICAL STRUCTURAL ANALYSIS OF DIGITIZED BINARY IMAGES BY BORDER FOLLOWING, COMPUTER VISION, GRAPHICS, AND IMAGE PROCESSING, 1985

Also Published As

Publication number Publication date
EP3815037A1 (fr) 2021-05-05
WO2020002859A1 (fr) 2020-01-02
US20210264596A1 (en) 2021-08-26

Similar Documents

Publication Publication Date Title
EP3755197B1 (fr) Dispositif et procédé de réalisation d&#39;un classifieur numérique d&#39;images, afin de déterminer la qualité de visualisation des images de vidéocapsule endoscopique d&#39;un segment du tube digestif
CN112075922B (zh) 2型糖尿病眼底图像指标的测量及与糖尿病肾病相关性的分析方法
EP4128140A1 (fr) Système et procédé pour produire un indicateur à partir d&#39;une image d&#39;une coupe histologique
FR2862788A1 (fr) Procede et systeme d&#39;analyse de changements temporels d&#39;image
CN108038875B (zh) 一种肺部超声图像识别方法和装置
CN102481117A (zh) 用于探测在三维重建中的不良质量的系统和方法
CN116703915B (zh) 基于多种图像融合的糖尿病视网膜病变会诊方法及系统
FR3062498B1 (fr) Systeme et procede pour l&#39;evaluation de risques vasculaires
Li et al. Automatic quantification of cone photoreceptors in adaptive optics scanning light ophthalmoscope images using multi-task learning
WO2020002859A1 (fr) Système et procédé pour produire un indicateur graphique multiparamétrique à partir d&#39;une image d&#39;une coupe histologique
US11776115B2 (en) System and method for estimating a quantity of interest based on an image of a histological section
WO2017009070A1 (fr) Procede d&#39;aide au diagnostic medical
CN109214433A (zh) 一种卷积神经网络区分肝癌分化等级的方法
US20240153288A1 (en) Systems and methods for quantitative polarization imaging
Razaghi et al. Correction of retinal nerve fiber layer thickness measurement on spectral-domain optical coherence tomographic images using U-net architecture
EP3815038A1 (fr) Systeme et procede pour produire un indicateur graphique multiparametrique a partir d&#39;une image d&#39;une coupe histologique
CA3017720A1 (fr) Procede de diagnostic de maladies inflammatoires chroniques de l&#39;intestin
TW201309264A (zh) 自動化甲襞微循環分析系統
CN111325631A (zh) 一种基于眼底图像的健康保险控费的方法和系统
FR3051950A1 (fr) Procede de caracterisation geometrique des voies aeriennes d&#39;un poumon par irm
US20200349700A1 (en) System and method for estimating a quantity of interest based on an image of a histological section
CN106370618A (zh) 一种双波长光谱差异度指数人体舌质舌苔分离方法
Qiang et al. Vessel segmentation and quantification in ultra-wide-field OCTA images
FR2997536A1 (fr) Procede et systeme d&#39;estimation d&#39;un indice de severite differentiel d&#39;une pathologie exprimant une alteration cutanee
Tanveer et al. Cancer image quantification with and without, expensive whole slide imaging scanners

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 2

PLSC Publication of the preliminary search report

Effective date: 20200103

RX Complete rejection

Effective date: 20200611