Brevet FR3079617A1 — publié le 20191004

METHODE DE DETECTION DE CELLULES PRESENTANT AU MOINS UNE ANOMALIE DANS UN ECHANTILLLON CYTOLOGIQUE La présente invention concerne une méthode de détection de cellules présentant au moins une anomalie dans un échantillon cytologique, à partir d'une image numérisée ou de plusieurs images numérisées du même échantillon. L'invention se rapporte également à un dispositif mettant en oeuvre cette méthode. L'invention se rapporte au domaine de la cytologie qui est l'étude par microscopie des cellules d'un échantillon cytologique après fixation sur une lame. Plus précisément, il s'agit de l'étude des aspects morphologiques, radiométriques ou biochimiques de cellules normales ou atypiques, c'est-à-dire présentant au moins une anomalie. Le dépistage du cancer est un des domaines d'application de la cytologie, et consiste à rechercher et à mettre en évidence un cancer à l'aide de prélèvements d'échantillons cytologiques sur un sujet. Le dépistage est réalisé en partie dans un laboratoire d'anatomopathologie par des cytotechniciens et cytopathologistes formés aux méthodes de détection de cellules atypiques. La cytologie est très utilisée pour la détection de certains cancers car elle permet de les diagnostiquer à différents stades, certains cancers ayant par exemple un stade précoce dit « bas grade » et un stade plus avancé dit « haut grade » (le grade peut aussi correspondre à un type de cancer). Il est évident que le diagnostic précoce d'un cancer est un grand enjeu avec d'importantes répercussions liées au taux de survie des sujets atteints et également au coût correspondant pour soigner ces derniers. Par exemple, chaque année dans le cadre du cancer de la vessie, 10 millions de cytologies urinaires sont effectuées dans le monde pour aider à le diagnostiquer, l'analyse de lames de cytologie d'échantillons provenant de mictions spontanées pour diagnostiquer un cancer de la vessie étant une pratique courante. 2 Cependant, cette méthode de diagnostic par l'analyse d'échantillon cytologiques est limitée d'une part par l'aspect humain, c'est-à-dire l'expérience et le nombre décroissants de cytopathologistes et de cytotechniciens, et d'autre part par une faible sensibilité pour certains 5 cancers. Typiquement, il existe une très faible sensibilité pour le diagnostic du cancer bas grade de la vessie car les cellules atypiques ou anormales n'ont pas encore subi de déformations morphologiques prononcées et sont alors difficilement détectables par l'oeil humain. D'où le besoin de rechercher des procédés utilisables en routine et peu coûteux qui permettraient 10 d'augmenter la performance et également la fiabilité de l'ensemble de cette pratique. A partir de cette situation, un des buts de l'invention consiste à fournir une méthode de détection de cellules présentant au moins une anomalie cellulaire dans un échantillon cytologique permettant d'aider par la 15 suite à la détection de cancer, notamment du cancer de la vessie, à partir d'image de lames numérisées. L'invention permet avantageusement, à l'aide d'un module de traitement mis en oeuvre par ordinateur, une analyse automatisée et exhaustive des cellules de l'échantillon présent sur la lame. Ainsi, l'analyse de l'échantillon cytologique et de toutes les cellules 20 potentiellement anormales est complète, rapide, plus précise et permet de palier les erreurs d'analyse courantes des cytopathologistes et des cytotechniciens. De plus, cette méthode assure une standardisation du résultat de l'analyse cytologique de l'échantillon. Mais encore, cette méthode permet en outre d'aider au diagnostic précoce de cancer dit bas grade, 25 notamment du cancer bas grade de la vessie, en détectant des cellules présentant de faibles anomalies. L'invention vise à remédier aux inconvénients de l'état de la technique cités ci-dessus et est motivée par l'absence de solutions alternatives. Pour atteindre ce but ou d'autres, un premier aspect de 30 l'invention propose une méthode de détection de cellules présentant au 3 moins une anomalie dans un échantillon cytologique à partir d'au moins une première image de microscopie électronique numérisée de l'échantillon, caractérisée en ce que la méthode est exécutée à l'aide d'un module de traitement mis en oeuvre par ordinateur, notamment entièrement exécutée 5 par le module de traitement mis en oeuvre par ordinateur. La méthode comprend une étape de traitement de la première image, comprenant les sous-étapes suivantes : - à partir de la première image, détecter chaque élément présent dans l'échantillon, choisi parmi au moins une cellule isolée ou un groupement 10 de cellules, à l'aide d'une détection colorimétrique, - filtrer chaque cellule isolée et/ou chaque groupement de cellules à l'aide de la détection colorimétrique, - à partir de chaque groupement de cellules détecté, individualiser au moins une cellule avec son noyau, à l'aide d'au moins une méthode de 15 segmentation fondée sur une détection des contours dans le groupement de cellules, De plus, la méthode comprend en outre une étape de décèlement d'anomalie cellulaire, comprenant les sous-étapes suivantes : - pour chaque cellule isolée ou individualisée dans l'étape précédente, 20 calculer des données d'au moins deux caractéristiques d'anomalies distinctes dites « caractéristiques expertes d'anomalie», et déterminer si ladite cellule présente au moins une anomalie, en comparant ces données calculées avec des données de références, - puis compter le nombre total de cellules dans l'échantillon présentant 25 au moins une anomalie, Et enfin, elle comprend en outre une première étape de caractérisation de l'échantillon suivant un niveau d'anomalie de l'échantillon, comprenant les sous-étapes suivantes : 4 - comparer le nombre total de cellules présentant au moins une anomalie à une valeur seuil du niveau d'anomalie, de préférence une valeur seuil égale à au moins trois cellules présentant chacune au moins une anomalie, 5 - si le nombre total de cellules présentant au moins une anomalie est inférieur à la valeur seuil du niveau d'anomalie, classer l'échantillon dans une classe représentative d'un échantillon avec un niveau d'anomalie nul, - si le nombre total de cellules présentant au moins une anomalie est 10 supérieur à la valeur seuil du niveau d'anomalie, classer l'échantillon dans une classe représentative d'un échantillon avec un niveau d'anomalie avéré, c'est-à-dire un niveau d'anomalie faible ou fort. Selon cet aspect de l'invention, la méthode de détection classe, in fine, l'échantillon cytologique selon un niveau d'anomalie cellulaire à partir de 15 l'analyse d'au moins une première image, de préférence avec une seconde image. Le module de traitement utilisé dans l'invention est apte à combiner des données calculées à partir de la première image de microscopie électronique à transmission numérisée et/ou d'une seconde image de microscopie à fluorescence numérisée de l'échantillon avec des données de 20 références connues en cytologie. On comprendra que lors de l'étape de traitement de la première image citée précédemment, plusieurs éléments constitutifs de la première image sont détectés, notamment des cellules isolées, des groupements de cellules, des éléments à rejeter ou des amas de cellules non analysables et qui sont donc 25 écartés pour la bonne poursuite de la méthode. Cette solution permet de bénéficier du classement de l'échantillon selon deux classes et de diagnostiquer si l'échantillon a un niveau d'anomalie cellulaire nul ou avéré. Par la suite, ce résultat permettra d'aider significativement le diagnostic d'un cancer ou non par un spécialiste. 5 Par ailleurs, on notera que selon l'invention une caractéristique experte d'anomalie d'une cellule peut être considérée comme un bon indicateur de la conformité d'une cellule. De préférence, une des caractéristiques expertes d'anomalie est un 5 rapport surfacique du noyau dans la cellule ou une irrégularité du noyau ou une colorimétrie mottée ou une hyper colorimétrie du noyau. On entend par colorimétrie mottée, la mesure de l'intensité de coloration de certains points ou de tâches. De préférence encore, la caractéristique experte d'anomalie est une 10 caractéristique du noyau, en particulier, un rapport de taille du noyau au regard de la taille de la cellule, un positionnement du noyau dans la cellule, ou une caractéristique de colorimétrie du noyau. Plus précisément la caractéristique experte d'anomalie est un rapport surfacique du noyau dans la cellule ou une irrégularité du noyau ou une colorimétrie mottée ou une 15 hyper colorimétrie du noyau. Avantageusement, au moins une des caractéristiques expertes d'anomalie est une irrégularité du noyau, le calcul de ladite irrégularité du noyau est un calcul d'écart correspondant à la comparaison entre chaque noyau de cellule analysée et un noyau type dit de référence d'une cellule dite 20 de « référence » en cytologie. Ce calcul d'écart est choisi parmi la liste suivante: une différence de surface, une différence de périmètre, une différence de la convexité du noyau, une différence de barycentre, une différence entre une longueur du contour de chaque noyau détecté, une différence d'irrégularité d'un contour du noyau, une différence du nombre de 25 pixels dans chaque noyau et une différence d'excentricité du noyau dans la cellule. Selon un mode de réalisation particulier, au moins une caractéristique experte d'anomalie est une longueur du contour du noyau, le calcul de cette caractéristique comprend en outre les sous-étapes suivantes : 6 - calculer un contour réel du noyau en utilisant une méthode d'extraction de contour du noyau, - déterminer un diamètre du noyau, de préférence le plus grand diamètre, en traçant une enveloppe englobante du noyau, 5 - tracer une courbe convexe fermée de rang C2 dont l'enveloppe englobante est celle d'un noyau de référence, - calculer une différence entre le contour réel du noyau et la courbe convexe tracée, représentant le contour d'un noyau de référence. Selon un mode de réalisation particulier, au moins une 10 caractéristique experte d'anomalie est une irrégularité d'un contour du noyau, le calcul de cette caractéristique comprend en outre les sous-étapes suivantes : - détecter au moins un point anguleux dans le contour du noyau détecté, 15 - calculer un nombre total de points anguleux dans le contour du noyau, - mesurer un angle pour chaque point anguleux correspondant, puis - déterminer un degré d'irrégularité du contour du noyau en combinant les données du nombre de points anguleux détecté par noyau et leurs angles. 20 Selon un mode de réalisation particulier, la méthode selon l'invention comprend une étape de sélection typologique, comprenant les sous-étapes suivantes: - à partir de la première image numérisée de l'échantillon, mesurer la taille des cellules, puis classer les cellules selon leur taille, 25 - associer chaque classe de cellules à un type cellulaire correspondant, puis sélectionner au moins un type cellulaire choisi correspondant à une taille de cellule, tel que le type cellulaire urothélial. 7 Selon ce mode de réalisation particulier, l'invention peut comprendre en outre une seconde étape de caractérisation de l'échantillon suivant le niveau d'anomalie des cellules de l'échantillon. Cette seconde étape comprend en outre les sous-étapes suivantes: 5 - à partir d'une seconde image de fluorescence numérisée de l'échantillon, sélectionner toutes les cellules d'un type cellulaire choisi, tel qu'urothélial, selon leur localisation dans l'échantillon prédéfinie au préalable dans la première image numérisée, - détecter la présence d'un halo de fluorescence autour de chaque 10 cellule sélectionnée précédemment, de préférence mesurer un taux de fluorescence, de préférence encore mesurer un taux d'étalement de la fluorescence, - si un halo de fluorescence est détecté, de préférence si un taux minimum de fluorescence est mesuré, de préférence encore un taux 15 minimum d'étalement de la fluorescence est mesuré, classer l'échantillon dans une classe représentative d'un échantillon avec un niveau d'anomalie faible. Selon un autre mode de réalisation particulier, l'invention comprend une étape de recommandation d'exploitabilité de l'échantillon selon un 20 nombre de cellules comprenant un noyau et selon une valeur seuil minimum d'exploitabilité, de préférence la valeur seuil minimum d'exploitabilité est égale à au moins sept cellules comprenant un noyau dans la première image. Selon l'invention, la méthode, à l'aide du module de traitement, 25 combine l'ensemble des données obtenues et classe l'échantillon cytologique analysé à partir de la première image et de la seconde image, dans l'une des trois classes associée au niveau d'anomalie cellulaire, nul, faible ou fort. 8 Un deuxième aspect de l'invention propose une méthode de détection des cellules présentant au moins anomalie dans un échantillon cytologique à partir d'une première image et/ou d'une deuxième image à l'aide d'une méthode d'apprentissage profond, à l'aide d'un module de 5 traitement mis en oeuvre par ordinateur, pour affecter l'échantillon à l'une des trois classes associée à un niveau d'anomalie cellulaire de l'échantillon, nul, faible ou fort. Cette méthode est avantageuse en ce qu'elle utilise les données fournies par au moins une des variantes du premier aspect de l'invention. 10 De préférence, l'échantillon cytologique est un échantillon urinaire de miction spontanée. Un troisième aspect de l'invention propose un dispositif mettant en oeuvre une méthode telle que décrite ci-dessus. Un quatrième aspect de l'invention propose une utilisation de la 15 méthode telle que décrite ci-dessous pour aider au diagnostic du cancer de la vessie. D'autres particularités et avantages de la présente invention apparaîtront dans la description ci-après d'exemples de réalisation non limitatifs, en référence aux dessins annexés dans lesquels : 20 - la figure 1 est un schéma illustrant les différentes étapes de la méthode mise en oeuvre selon un mode préféré de l'invention, - la figure 2 représente un premier exemple d'une première image de microscopie électronique numérisée d'un échantillon avant et après une étape de traitement. La première image comprend ici 25 notamment des cellules isolées, des amas cellulaires, des groupements de cellules et des éléments à rejeter, - la figure 3 représente un deuxième exemple d'une première image de microscopie électronique numérisée d'un deuxième échantillon avant et après une étape de traitement de la première 9 image. La première image comprend ici notamment des cellules isolées, des amas cellulaires, des groupements de cellules et des éléments à rejeter. - la figure 4 représente deux images de microscopie électronique 5 d'une cellule isolée et deux images de microscopie électronique de leurs noyaux extraits à l'aide d'un module de traitement mis en oeuvre par ordinateur, - la figure 5 représente des images de microscopie électronique d'une cellule isolée, de son noyau extrait et d'un exemple de tracé 10 pour mesurer le rapport taille noyau cellule, à l'aide du module de traitement mis en oeuvre par ordinateur, - la figure 6 est un schéma illustrant les différentes étapes de la méthode mis en oeuvre selon un autre mode de réalisation selon lequel elle comprend une deuxième étape de caractérisation de 15 l'échantillon à partir d'une seconde image numérisée, - la figure 7 est un schéma illustrant un deuxième aspect de l'invention, plus précisément une méthodologie d'implémentation pour un apprentissage profond appliqué à la classification d'au moins une image d'échantillon cytologique correspondant à un 20 ensemble d'imagettes. Il faut tout d'abord noter que les figures exposent l'invention de manière schématique. Ces figures représentent des exemples de réalisation donnés à titre non limitatifs et ne respectent pas nécessairement une chronologie particulière. 25 La présente invention s'inscrit dans le cadre général de l'application de la vision par ordinateur d'images numérisées d'objets complexes et nombreux. L'invention traite d'une méthode originale pour détecter les cellules présentant au moins une anomalie dans un échantillon cytologique à partir d'une première image numérisée, et de préférence à l'aide d'une 10 seconde image numérisée du même échantillon. La première image est obtenue à l'aide d'un microscope électronique en transmission et la seconde image est obtenue à l'aide d'une microscopie en fluorescence. Les images utilisées sont issues de la numérisation particulière d'un échantillon 5 cytologique ayant subi une préparation adaptée pour la méthode, qui a ensuite été étalée sur une lame de cytologie. La méthode selon l'invention peut utiliser un échantillon cytologique provenant indifféremment de fluide ou de prélèvement d'éléments du corps humain selon l'analyse souhaitée. Pour le mode de réalisation présenté ici, on prend un échantillon d'urine de miction 10 spontanée. Dans ce contexte, le résultat de cette méthode de détection de cellules présentant au moins une anomalie dans un échantillon cytologique permettra d'aider au diagnostic d'un cancer de la vessie, notamment à un stade précoce ou dit bas grade du cancer de la vessie. Dans la suite de la description et tel qu'illustré sur la figure 1, on va 15 décrire un mode de réalisation préféré comprenant plusieurs étapes pouvant être préliminaire, intermédiaire ou successive. On comprendra que selon le mode de réalisation choisi on pourra ajouter ou enlever une ou plusieurs de ces étapes à la méthode selon l'invention. Avantageusement, la méthode de détection selon l'invention 20 comprend des étapes exécutées à l'aide d'un module de traitement mis en oeuvre par un ordinateur apte à combiner des données calculées à partir de l'échantillon avec des données de références connues en cytologie. Les données dites calculées, sont calculées et recueillies à partir d'au moins une image numérisée d'une lame d'au moins un échantillon par le module de 25 traitement. A l'aide d'une méthode d'apprentissage le module de traitement peut ainsi combiner toutes les données calculées et recueillies à partir de l'échantillon ou de plusieurs échantillons, et avec des données de références connues en cytologie. Selon le mode de réalisation préféré présenté ici, la méthode de 30 détection selon l'invention comprend une étape préliminaire de traitement 11 d'image 110. Cette étape préliminaire de traitement d'image 110, illustrée dans les figures 1 à 3, a pour but d'affiner au maximum la sélection de départ des cellules pour la détection de cellules, notamment en écartant des éléments 4 constitutifs d'une première image numérisée 1 qui pourraient 5 interférer et diminuer la finesse du résultat de la méthode de détection. Pour ce faire, cette étape préliminaire de traitement 110 comprend une sous-étape de détection 112 des éléments 4 présents dans la première image 1. Cette sous-étape de détection 112 est au moins en partie basée sur une détection colorimétrique pouvant consister au moins en partie à une segmentation en 10 plusieurs groupes de l'ensemble des couleurs de la première image 1, le groupe le plus clair et le plus répandu correspondant à la couleur du fond de l'image 6 et le groupe le plus foncé correspondant à la couleur des noyaux cellulaires. On considéra que tout pixel n'étant pas un pixel du fond de l'image 6 est un pixel d'un élément 4, c'est-à-dire d'un objet biologique. Ainsi, 15 le fond de l'image 6 est déterminé et tout ce qui ne fait pas partie du fond de l'image 6 correspond à un élément 4, c'est-à-dire, une cellule, un groupement de cellules, un amas de cellules ou un élément à rejeter. Puis, tous les éléments 4 sont extraits du fond de l'image 6 en créant plusieurs imagettes 8 à l'aide d'une méthode de segmentation. Pour catégoriser les 20 différents éléments 4 détectés, on utilise en outre cette détection colorimétrique, détectant les noyaux, ici cette sous-étape est dite « grossière ». De plus, on pourra combiner cette première sous-étape de détection des noyaux, c'est-à-dire de détection colorimétrique à des mesures des 25 dimensions de chacun des éléments 4 détectés. Par exemple, si plusieurs noyaux sont détectés par détection colorimétrique la mesure de cet élément aidera à identifier le type d'élément 4 détecté. A l'issue de la sous-étape de détection 112, il est ainsi possible d'identifier différentes pré-classes des éléments 4 constitutifs de la première image 1 à 12 l'aide de la détection colorimétrique des noyaux, et éventuellement de la mesure de ces éléments 4, tels que : - Eléments à rejeter : C'est tout ce qui ne constitue pas une cellule à noyau. En effet, tout élément sans noyau est à rejeter, car non 5 interprétable, cela peut être soit un déchet, soit une fausse détection. - Cellules isolées : Ce sont les imagettes qui ne contiennent qu'un seul noyau et qui sont par conséquent des cellules isolées. - Groupement de cellules : Ce sont des cellules regroupés (qui se touchent ou qui se chevauchent) et qui contiennent de ce fait plusieurs 10 noyaux. - Amas : Ce sont des regroupements particuliers de cellules très denses (très foncés) et où les cellules se chevauchent d'une façon bien particulière. Ce qui différencie les amas des groupements classiques de cellules est le fait d'avoir dans les amas une superficie 15 de la classe "noyau" plus importante. De façon plus détaillée, la figure 3 illustre un exemple de première image 1 subissant l'étape préliminaire de traitement 110. Tel que décrit plus tôt, lors de cette étape 110, la première image 1 subit une segmentation pour obtenir les imagettes 8 entourant les éléments 4 présentant une couleur foncé, ces 20 imagettes 8 sont ensuite pré classées en quatre classes : éléments à rejeter 12, cellules isolées 10, groupement de cellules 14 et amas de cellules 16 (qui seront rejetées également). Pour le bon déroulement de l'invention, la méthode de détection, selon l'invention doit se faire à partir de cellules nécessairement isolées ou individualisées à partir d'un groupement de 25 cellules et comprenant un noyau. Les cellules isolées 10 et les groupements de cellules 14 sont filtrés lors d'une sous-étape de filtrage 114 utilisant la détection colorimétrique. En conséquence, la pré classe groupement de cellules 14 subit une sous-étape d'individualisation 116 après la sous-étape de filtrage 114. Cette 13 sous-étape 116 consiste à individualiser les cellules du groupement de cellules 14. Par individualiser, on entend extraire du groupement de cellules 14 au moins une cellule comprenant un noyau. On notera que cette sous-étape d'individualisation 116 n'est pas exécutée pour la pré-classe des 5 « amas de cellules» 16. En effet, les amas de cellules 16 sont des structures trop complexes à individualiser et leur prise en compte engendrerait beaucoup plus d'erreurs que d'informations pertinentes à utiliser. C'est pour cela que les amas 16 ne sont pas utilisés selon l'invention. La sous-étape d'individualisation 116 des cellules se fait à l'aide d'une nouvelle 10 segmentation qui est fondée, sur la détection de contours via le calcul du gradient de la première image ou via l'utilisation d'une transformée mathématique qui a un aspect multidirectionnel telles que Radon, Gabor ou les transformées multi échelles multidirectionnelles. Après avoir déterminé les contours et éventuellement à l'aide d'une approche de morphologie 15 mathématique, on extrait toute structure comprise dans un contour fermé, correspondant donc à une cellule, et on la définit comme étant une cellule individualisée, intéressante pour la suite de la méthode selon l'invention. Selon le mode de réalisation préféré ici, une étape optionnelle de recommandation d'exploitabilité 120 de l'échantillon pour détecter des 20 cellules potentiellement cancéreuses peut être exécutée. Selon les modes de réalisation, cette étape de recommandation d'exploitabilité 120 peut intervenir une ou plusieurs fois à plusieurs niveaux durant le déroulement de la méthode selon l'invention, de préférence à la fin de l'étape préliminaire de traitement 110. Cette étape de recommandation d'exploitabilité 120 revient à 25 estimer si la première image 1 contient un nombre de cellules comprenant un noyau, suffisant pour que la détection dans l'échantillon ne soit pas biaisée. Pour ce faire, le nombre total de cellules avec noyau présentes dans la première image 1 est déterminé. Ce nombre total de cellules sera communiqué à un utilisateur qui validera ou non la poursuite de la méthode 30 selon l'invention. De préférence, ce nombre total de cellules avec noyau est comparé à une valeur d'un seuil minimum d'exploitabilité qui sera déterminé 14 au préalable par l'utilisateur. En variante, l'utilisateur pourra rentrer la valeur du seuil minimum d'exploitabilité de l'échantillon souhaitée dans le module de traitement pour permettre la poursuite automatique ou non de la méthode inventive. Par conséquent, la méthode de détection peut se poursuivre 5 même s'il n'y a que très peu de cellules avec noyau dans la première image 1, dans la mesure où le nombre total de ces cellules est supérieur ou égal à la valeur du seuil minimum d'exploitabilité. Par exemple, la valeur du seuil minimum d'exploitabilité est égale à au moins sept cellules avec noyau détectées dans l'échantillon, ainsi la méthode se poursuit si le nombre total 10 de cellules avec noyau est supérieur ou égal à sept. Inversement, si le nombre de cellules avec noyau est insuffisant l'étape de recommandation d'exploitabilité 120 peut interrompre la méthode de détection. Selon le mode de réalisation préféré ici, une étape optionnelle préliminaire de sélection typologique 130 pour détecter des cellules d'un type 15 cellulaire choisi dans l'échantillon avec au moins une anomalie cellulaire, plus précisément, potentiellement cancéreuses, peut être exécutée. Pour ce faire, à partir de la première image 1, les cellules de l'échantillon sont mesurées puis classées selon leurs tailles, étape 132. A l'aide du module de traitement, les cellules classées selon leurs tailles sont affiliées à des types 20 cellulaires, puis sélectionnées selon le type cellulaire choisi, tel qu'urothélial ou tout autre type cellulaire dont un niveau d'anomalie cellulaire peut-être corrélé à un type de cancer, étape 134. Selon une variante, l'utilisateur pourra définir lui-même le type cellulaire préféré. Par exemple, pour une méthode de détection appliquée à des échantillons de cytologie urinaire pour 25 aider au diagnostic du cancer de la vessie, le type cellulaire préféré sera le type urothélial. De façon générale, l'étude des cellules urothéliales est un bon indicateur pour aider au diagnostic du cancer de la vessie, notamment pour déterminer le stade du cancer en fonction des anomalies rencontrées. Suite à l'étape préliminaire de traitement d'image 110, à l'étape 30 optionnelle de recommandation d'exploitabilité de l'échantillon 120, et à 15 l'étape optionnelle préliminaire de classification typologique 130 on obtient une sélection fine de cellules isolées ou individualisées, comprenant un noyau et d'un type cellulaire choisi. Pour rappel, la méthode selon l'invention aboutit à la détection de 5 cellules présentant au moins une anomalie dans l'échantillon cytologique à partir d'au moins une première image. Le niveau d'anomalie cellulaire de l'échantillon est corrélé au nombre et au type d'anomalie cellulaire rencontré dans l'échantillon cytologique. En d'autres termes, la méthode selon l'invention peut permettre de classer l'échantillon soit dans une classe de 10 niveau d'anomalie cellulaire dit avéré, soit dans une classe de niveau d'anomalie cellulaire dit nul. Pour ce faire, le mode de réalisation préféré illustré sur la figure 1 comprend en outre une étape de décèlement d'anomalie cellulaire 140 dans l'échantillon. Selon l'invention, une cellule présentant une anomalie est 15 définie comme une cellule présentant une caractéristique experte d'anomalie divergente d'une donnée dite de référence, c'est-à-dire divergente de la norme. En d'autres termes, il s'agit d'une cellule qui présente au moins une anomalie c'est-à-dire une caractéristique dite experte d'anomalie qui est anormale. Ainsi le calcul d'une caractéristique experte d'anomalie permet 20 d'obtenir une donnée qui permet d'aider à déterminer un écart avec la norme et donc la présence ou non d'une anomalie. On entend par le calcul de caractéristiques expertes d'anomalie, la mesure ou tout autre type de méthodes mathématiques. Selon l'invention, la caractéristique experte d'anomalie qui est calculée est une caractéristique du noyau, en particulier 25 elle correspond à un rapport de taille entre le noyau et le cytoplasme de la cellule ou à une mesure du positionnement du noyau ou à une mesure colorimétrique du noyau. Plus précisément, la caractéristique experte d'anomalie correspond à un rapport surfacique du noyau dans la cellule, ou à une irrégularité du noyau, ou à une colorimétrie mottée ou une hyper 30 colorimétrie du noyau. Les différentes caractéristiques expertes d'anomalie 16 mentionnées ci-dessus, à titre non limitatif, peuvent révéler des anomalies selon les valeurs calculées dans l'échantillon. On comprendra que d'autres types de caractéristiques expertes d'anomalies pourront être indifféremment calculées et utilisées comme paramètres pour classer les cellules. Les anomalies recherchées et détectées dans chaque cellule vont, plus particulièrement, concerner son noyau. Pour ce faire, il faut donc une détection fine du noyau pour pouvoir détecter les anomalies qui lui sont liées. La détection fine du noyau se fait par une segmentation. Cette segmentation peut notamment se faire par affinement à partir de la segmentation grossière 10 utilisant uniquement la colorimétrie. Cette technique décrite plus tôt dans la description consiste à partir d'une enveloppe englobante, de préférence formant un rectangle autour du noyau. Cette enveloppe englobante délimitera une zone autour du noyau, toujours incluse dans la cellule, et débordant du noyau de 10 à 40% (de la surface totale du noyau), de 15 préférence entre 20 à 30 %. Les figures 4 et 5 illustrent chacune le résultat d'une détection fine de noyau. Selon le mode de réalisation, le noyau peut déjà avoir été détecté à l'étape préliminaire de traitement de l'image 110 lors de la détection grossière des noyaux à l'aide d'une première segmentation. Un exemple de réalisation 20 consiste à prendre en compte l'enveloppe englobante utilisée pour la détection grossière du noyau et d'estimer dans cette enveloppe la distribution de couleur du noyau ainsi que celle de l'arrière-plan en utilisant un modèle de mélange de gaussiennes. Par la suite, on l'utilise pour construire un champ aléatoire de Markov sur les étiquettes des pixels, avec une fonction 25 d'énergie qui préfère les régions connectées ayant la même étiquette, et exécutant une optimisation basée sur le Graph cut pour interférer leurs valeurs. Comme cette estimation est susceptible d'être plus précise que l'original extraite de la boîte englobante, ces deux étapes sont répétées jusqu'à la convergence. 17 L'étape de décèlement d'anomalie cellulaire 140 dans les cellules de l'échantillon se fait en partie par le calcul d'au moins deux caractéristiques expertes d'anomalie distinctes, étape 142. Dans la suite de la description, on décrit le calcul de quatre types de 5 caractéristiques expertes d'anomalies mettant chacune en évidence la présence ou non d'une anomalie cellulaire morphologique. Premièrement, on va décrire le calcul de la caractéristique experte d'anomalie du rapport noyau/cytoplasme dans la cellule. On se propose de calculer à l'aide du module de traitement automatisé le 10 rapport RNC (Rapport Noyau Cellule) qu'on définit par les formules suivantes: RNC = longueur du plus grand diamètre du noyau / longueur du plus grand diamètre de la cellule correspondant. La figure 5 illustre un exemple de réalisation du calcul du rapport RNC. 15 Les étapes du calcul sont les suivantes à l'aide du module de traitement : 1. Tracer le plus grand axe correspondant à la cellule, 2. Calculer la longueur de l'axe délimitée par la cellule, 3. Calculer la longueur de ce même axe délimité par le noyau, On propose également une autre façon de calculer ce rapport en se basant 20 sur le calcul de surfaces. RNC=(Surface du noyau/surface de la cellule)112 Avec : Surface de la cellule = Nombre de pixels de la cellule. Surface du noyau = Nombre de pixels du noyau. 25 Ce rapport est un chiffre inférieur à 1 puisque le noyau est à l'intérieur de la cellule. Si le RNC est très proche de 1, la cellule est considérée comme étant « un noyau nu » et sera donc rejetée lors de la méthode par le module de 18 traitement. On se propose alors de donner la main à l'utilisateur pour qu'il saisisse un seuil à partir duquel la cellule est considérée comme « un noyau nu ». De même pour les seuils indiquant un RNC normal, critique ou anormal. Par défaut, on initialisera ces seuils avec les valeurs suivantes : 5 a) RNC 0,7 correspondant à un RNC anormal d) RNC > 0,9 correspondant à un noyau nu non pris en compte « à rejeter » 10 Plus précisément, le rapport noyau/cellule est calculé à partir du rapport entre le plus grand diamètre de la cellule et le plus grand diamètre du noyau. Dans le cas où la cellule n'est pas très sphérique, on ne prend pas le diamètre réel mais la distance ou longueur de bout en bout de la cellule, correspondant au plus grand diamètre de la cellule. Tous rapports compris 15 entre 0,9 et 1 ne seront pas interprétables puisque cela correspondra à des noyaux nus. Les rapports supérieurs à 0 ,7 correspondent à une forte anomalie et les rapports compris entre 0,5 et 0,7 correspondent à une faible anomalie. Deuxièmement, on va décrire le calcul d'une autre caractéristique 20 experte d'anomalie qui est l'irrégularité du noyau. On entend par noyau irrégulier, tout noyau dont le contour dévie d'une forme classique de contour de noyau dit de référence bien lisse et régulière. On peut le voir aussi comme un contour ayant des points anguleux ou coins. Pour ce faire, on se propose de calculer l'irrégularité du noyau via trois 25 méthodes dont la troisième méthode décrite est la fusion des deux autres méthodes. 19 Dans un premier temps, on va décrire la première méthode basée sur le calcul de différence avec un noyau de référence d'une cellule de référence. La différence calculée peut être une différence surfacique ou une différence de périmètres ou n'importe quelle autre différence telle que la différence de 5 barycentre, la différence de pixel dans chaque région du noyau, la différence entre la longueur des contours de chaque noyau (de référence et réel), la différence linéaire entre le centre de gravité de chaque noyau, la différence entre le nombre de pixels des régions convexes liées à chaque noyau, la différence entre l'excentricité du noyau réel et celle du noyau de référence. 10 Pour rappel, l'excentricité e=0 pour un cercle, 0 Moyenne (Zone entourant (x,y)) + 3* EcartType(Zone entourant (x,y)) alors (x,y) appartient à une zone non homogène, 10 Calculer une somme des points (ou blocs) non homogènes détectés par la cohérence (Nb coh), 4. Fusionner ces caractéristiques [V, Nb pics et Nb coh] et les utiliser pour détecter la colorimétrie mottée du noyau. On comprendra que l'on peut utiliser ces caractéristiques à l'entrée d'un 15 classifieur type SVM ou réseaux de neurones. Dans un deuxième temps, on va décrire la deuxième méthode basée sur l'estimation dans le domaine de la transformée: On se propose de définir une colorimétrie mottée par la présence de trois mottes au minimum dans la zone concernée, c'est-à-dire le noyau. Ici, l'idée 20 est de détecter un minimum de trois contours fermés à l'intérieure de la zone. Pour ce faire, on se propose de détecter les contours par une transformée appropriée comme la transformée en contourlette. Cette transformée a la particularité d'avoir des coefficients de faible énergie (proche de 0) dans les zones homogènes et des coefficients de très forte énergie dans les zones 25 anisotropes (non homogènes) et qui correspondent généralement aux contours de l'image. L'estimation de la colorimétrie mottée se fera en suivant les étapes les suivantes : 27 1. Faire la décomposition de l'image par une transformée multi échelle multidirectionnelle appropriée comme la transformée en contourlettes ou wave atoms, 2. Détecter les zones à fortes composante énergétique en appliquant un 5 seuillage pour ne garder que les coefficients correspondant aux contours de l'image, 3. Appliquer une fonction non linéaire de rehaussement de contraste aux coefficients en contourlettes retenus (seuillés), 4. Reconstruire l'image contour à partir des coefficients retenus en utilisant la 10 transformée inverse, 5. Détecter dans l'image contour les contours fermés et les énumérer, 6. S'il y a plus de trois contours fermés, on conclura que le noyau à une colorimétrie mottée. Dans un troisième temps, on va décrire la troisième méthode consistant à 15 fusionner la première et la seconde méthode de calcul de la colorimétrie mottée. De même que pour le calcul de l'irrégularité du noyau, on propose, ici, de combiner les résultats obtenus par chacune des deux méthodes précédentes pour obtenir un résultat fin pour la caractéristique experte d'anomalie de 20 colorimétrie mottée du noyau de chaque cellule. On comprendra qu'à l'aide du module de traitement, l'ensemble des caractéristiques expertes d'anomalies calculées (rapport noyau/cellule, irrégularité du noyau, hypercolorimétrie et colorimétrie mottée) sont combinées avec des données de référence consistant en des données de 25 références connues en cytologie. En particulier, ces données calculées peuvent être utilisées soit telles qu'elle, soit par un classifieur, soit fusionnées à d'autres caractéristiques calculées et générées par apprentissage profond pour déterminer si la cellule est normale ou présente au moins une anomalie, étape 144. Si l'échantillon ne comprend pas au moins une cellule présentant 28 au moins une anomalie, alors la méthode de détection s'arrête ici et si l'échantillon présente au moins une cellule avec au moins une anomalie, la méthode de détection continue. A la fin de l'étape de décèlement d'anomalie cellulaire 140, le module 5 de traitement détermine un nombre total de cellules présentant au moins une anomalie dans l'échantillon cytologique, étape 146. Puis, il peut intervenir, ici, une première étape de caractérisation de l'échantillon 150 suivant le niveau d'anomalie de l'échantillon dans sa globalité. Pour ce faire, le nombre total de cellules détectées présentant au 10 moins une anomalie est comparé à une valeur seuil d'un niveau d'anomalie de l'échantillon. Cette valeur seuil du niveau d'anomalie de l'échantillon est définie au préalable par l'utilisateur pour la méthode. Elle correspond à au moins trois cellules présentant chacune au moins une anomalie. Cette valeur seuil du niveau d'anomalie peut indifféremment être définie au préalable par 15 l'utilisateur selon la finesse du résultat de classification souhaitée. Selon l'invention, si le nombre total de cellule avec au moins une anomalie est inférieur à la valeur seuil du niveau d'anomalie, l'échantillon cytologique est classé dans une classe représentative d'un échantillon avec un niveau d'anomalie nul 152, si le nombre total de cellule avec au moins une anomalie 20 est supérieur ou égal à la valeur seuil du niveau d'anomalie, l'échantillon cytologique est classé dans une classe représentative d'un échantillon avec une classe avec un niveau d'anomalie non-négligeable 154, c'est-à-dire un faible niveau d'anomalie et un fort niveau d'anomalie 154. En bref, selon un autre aspect, l'invention à ce stade permet ainsi de classer 25 l'échantillon cytologique selon deux classes : une classe avec un niveau d'anomalie cellulaire nul 152 et une classe avec un niveau d'anomalie cellulaire non-négligeable 154. Ce résultat permet de mettre en évidence des anomalies cellulaires qui sont caractéristiques de cellules cancéreuses. Ainsi, il permet d'aider les cytotechniciens et cytopathologistes pour le 30 diagnostic du cancer, en apportant un résultat standardisé fiable et rapide. 29 Tel qu'illustré sur la figure 6, la méthode de détection peut comprendre une seconde étape de caractérisation de l'échantillon 160 suivant le niveau d'anomalie cellulaire de l'échantillon, mettant notamment en 5 évidence spécifiquement un niveau d'anomalie cellulaire faible dans l'échantillon. Ce niveau d'anomalie correspond à des cellules n'ayant pas encore subies de déformations morphologiques prononcées. De préférence, selon ce mode de réalisation particulier, l'invention permet de classer plus précisément l'échantillon soit dans une classe de niveau d'anomalie cellulaire 10 dit « faible », soit dans une classe de niveau d'anomalie cellulaire dit « fort », ou soit dans une classe de niveau d'anomalie cellulaire dit « nul ». Suite à la détection de cellule anormale, le classement de l'échantillon selon le niveau d'anomalie morphologique ou radiométrique de ces cellules de l'échantillon procure une réelle aide pour le médecin concernant le diagnostic 15 du cancer. En effet le niveau d'anomalie cellulaire dit « faible » peut être associé au stade précoce dit « bas-grade » du cancer et le niveau d'anomalie dit « fort» peut être associé au stade avancé du cancer dit « haut grade ». A cet effet, pour détecter les anomalies d'un niveau dit faible, l'invention 20 utilise une seconde image 2 numérisée d'imagerie de fluorescence. A partir de cette seconde image 2, on sélectionne les cellules d'au moins un type cellulaire choisi, tel qu'urothélial, dans une étape de sélection typologique 162, grâce à leurs localisations prédéfinies dans l'échantillon dans la première image 1 numérisée. La première image 1 et la seconde image sont 25 alors recalées, éventuellement à l'aide d'une méthode classique de recalage mis en oeuvre par ordinateur. On comprendra que pour exécuter cette seconde étape de caractérisation 160, l'étape préliminaire de sélection typologique 130 à partir de la première image 1 doit nécessairement avoir été exécutée plus tôt dans la méthode 30 selon l'invention pour pouvoir localiser les cellules du type cellulaire choisi. 30 Cette sous-étape de sélection typologique 162 d'un type cellulaire est suivie d'une sous-étape de détection d'un halo de fluorescence 164 autour de chaque cellule sélectionnée précédemment. On définit le terme « halo » comme un phénomène se traduisant par la présence d'une auréole 5 fluorescente qui entoure la cellule et qui est bien contrastée. Dans une variante de réalisation, la force du halo, le taux de fluorescence ou le taux d'étalement de la fluorescence est mesuré. Ainsi, si un halo de fluorescence est détecté, de préférence si un taux minimum de fluorescence, ou de préférence encore, un taux minimum d'étalement de la fluorescence est 10 détecté, on classe l'échantillon dans une classe représentative d'un échantillon avec un faible niveau d'anomalie 168. Inversement, les cellules ne présentant pas de fluorescence sont classées dans la classe de niveau d'anomalie nul et fort 166. Selon la présente invention, on propose trois méthodes pour effectuer l'étape 15 de seconde caractérisation de l'échantillon 160 à l'aide du module de traitement. On détaille ci-dessous, les sous-étapes de l'étape de seconde caractérisation de l'échantillon 160 dans le cadre d'une première méthode : - Utiliser la première image 1 et détecter le centre de chaque cellule d'un 20 type cellulaire choisi, à l'aide du calcul du barycentre de la cellule, - tracer un profil d'intensité le long de chaque segment de droite commençant du centre de la cellule jusqu'à la limite de chaque cellule selon différents angles. Par exemple, selon 8 angles : 0°, 45°, 90°, 135°,180°, 225°, 270° et 315°, 25 - détecter le pic Pmax (point maximum) pour chaque profil ainsi que le creux Pmin qui précède le maximum (le point minium), - tracer sur la courbe du profil d'intensité une ligne horizontale partant du creux Pmin situé à droite du pic Pmax jusqu'à l'intersection Pint avec la courbe du profil d'intensité à gauche du pic Pmax, 31 - calculer la distance linéaire horizontale dx (en abscisse) entre Xmax l'abscisse du pic et Xint l'abscisse du point Pint, - calculer la distance verticale dy (en ordonnée) entre Ymax l'ordonnée du pic et Yint=Ymin l'ordonnée de Pint et de Pmin, 5 - calculer pour chaque cellule détectée Dx et Dy qui sont la moyenne des dx et les dy selon les 8 angles. Dy représentera une première estimation de la force du halo de fluorescence du halo et Dx une première estimation de l'étalement du halo. On détaille ci-dessous les sous-étapes de l'étape de seconde caractérisation 10 de l'échantillon 160 dans le cadre d'une deuxième méthode : On propose ici pour mesurer et calculer l'étalement d'un halo de fluorescence de partir des profils Radon des imagettes des cellules. La transformée de radon représente une image comme une collection de projections dans diverses directions. Chaque projection selon une ligne 15 représente un profil Radon. Les étapes à suivre pour la deuxième méthode sont les suivantes : - repérer la position de chaque cellule d'un type cellulaire choisi dans la première image 1, - faire la correspondance dans la seconde image 2 pour extraire les 20 imagettes correspondantes à chaque cellule du type cellulaire choisi en dépassant de 10 à 30% en longueur et en largeur autour de chaque noyau, - calculer les profils Radon correspondants aux imagettes pour différents angles, de préférence les mêmes angles utilisés que dans la première méthode, 25 - pour chaque profil, détecter les deux pics Pmin et Pmax de chaque profil radon de chaque imagette dépassant de 10 à 30% en longueur et en largeur autour du noyau et suivre ensuite les étapes de la première méthode pour calculer une seconde estimation de l'intensité du halo et de son étalement. 32 Alternativement, dans le cadre d'une troisième méthode, l'étape de seconde caractérisation de l'échantillon 160 exécutée par le module de traitement, consiste en une méthode se basant sur le calcul de l'irrégularité du contour 5 du « halo » de fluorescence présent. En d'autres termes si une cellule présente un fort niveau d'irrégularité du contour du halo, cela signifie qu'elle présente un halo de fluorescence. Pour ce faire cette troisième méthode consiste en l'une des trois méthodes du calcul d'irrégularité du noyau présentée précédemment. 10 Pour l'ensemble des trois méthodes d'exécution de la seconde caractérisation de l'échantillon 160, décrites ci-dessus, à titre non limitatif, les résultats obtenus pour déterminer la classe de l'échantillon sont comparés à des valeurs seuils du niveau d'intensité de fluorescence et/ou d'étalement qui sont prédéterminées ou entrées au préalable par un utilisateur. 15 On comprendra que ces trois méthodes décrites peuvent être combinées et traité par méthode d'apprentissage pour obtenir un résultat fin. Cette étape de seconde caractérisation de l'échantillon 160 permet ainsi d'obtenir, in fine, la classification de la lame de l'échantillon, donc de l'échantillon dans trois classes possibles, c'est-à-dire: soit dans la classe de 20 niveau d'anomalie nul, pouvant être associé à un échantillon provenant d'un sujet sain ; soit dans la classe de niveau d'anomalie cellulaire faible, pouvant être associé à un échantillon présentant des cellules cancéreuses à un stade précoce ou dit cancer « bas-grade » ; soit dans la classe de niveau d'anomalie cellulaire fort, pouvant être associé à un échantillon présentant 25 des cellules cancéreuses à un stade avancé ou dit cancer « haut grade ». Un deuxième aspect de l'invention propose d'utiliser une deuxième méthode de détection d'anomalie cellulaire dans l'échantillon cytologique en utilisant une méthode d'apprentissage profond pour faire la classification directe d'une image de lame d'échantillon cytologique selon un niveau 30 d'anomalie cellulaire rencontré. De manière générale, une méthode 33 d'apprentissage profond apprend à un modèle informatique comment réaliser des tâches telles que la classification directement à partir d'images. Contrairement à ce qui se fait habituellement dans le domaine de l'apprentissage profond et du traitement d'image, l'originalité de cette 5 deuxième méthode est qu'elle classifie un ensemble d'imagettes sans faire de classification préalable de chacune des imagettes constituant l'ensemble de l'image. Plus précisément, cette méthode permet directement la classification de l'image de l'échantillon cytologique de la lame, c'est-à-dire de l'ensemble des cellules comprises dans l'image de l'échantillon sans 10 passer par une classification de chaque imagette ou cellule de l'échantillon. On notera que cette méthode de classification selon la présente invention est particulièrement utile lorsque la détection d'anomalies ou de pathologies peut seulement se faire en regardant l'ensemble d'une lame ou l'image de la lame, et non en regardant chaque cellule indépendamment les unes des 15 l'autres. Pour ce faire, cette deuxième méthode de détection de cellule présentant au moins une anomalie dans l'échantillon utilise une modèle informatique d'apprentissage profond, prenant un ensemble d'observations en entrée et produisant une classification de cet ensemble en sortie. Ces observations peuvent directement être des imagettes brutes de l'échantillon 20 cytologique et/ou une liste de valeurs caractérisant l'imagette (ou données de descripteurs), par exemple des valeurs calculées à l'aide de la première méthode (on comprendra que ces dernières données de descripteurs auront pour utilité d'apprendre au modèle informatique). Pour ce faire, cette deuxième méthode utilise une méthode d'apprentissage profond utilisant une 25 architecture de réseau de neurones profond (NN). Quand ce réseau de neurones profond (NN) est utilisé sur les imagettes brutes, il s'agit d'un réseau de neurones profond à convolution (CNN). Le CNN convolue les caractéristiques apprises avec les données d'entrées du descripteur, et utilise des couches à convolution 2D, parfaitement adaptée au traitement des 30 données 2D telles que les images. Le réseau de neurones à convolution peut extraire des représentations expertes des images directement via un 34 apprentissage. Néanmoins, la méthode d'apprentissage profond nécessite une puissance de calcul considérable. Les calculs étant hautement parrallélisables, ils sont généralement faits sur un processeur graphique ou GPU. Le GPU a de grandes performances et est doté d'une architecture 5 parallèle, qui est efficace pour la mise en oeuvre de cette méthode d'apprentissage. Le GPU est intégré dans une unité de traitement ou CPU. Toutefois, l'utilisation de GPU n'est pas possible en traitant chacune des entrées indépendamment car l'une des opérations les plus chronophages reste le transfert des 10 données de la mémoire au CPU et au GPU. De manière générale, pour pallier à ce problème, les méthodes d'apprentissage profond utilisent des tenseurs 4D (des tableaux à 4 dimensions) pour traiter de façon générale, les entrées par paquet ou batch. Grâce à une utilisation originale du tenseur 4D, cette deuxième méthode de 15 détection d'anomalie cellulaire classe l'image c'est-à-dire l'ensemble d'imagettes sans détection préalable de chacune d'entre elles. En d'autres termes, cette deuxième méthode de détection propose de détecter les cellules présentant au moins une anomalie cellulaire et de classer directement non pas une imagette mais un ensemble d'imagettes (l'imagette 20 brute et/ou caractérisée par des valeurs calculées d'une autre manière) représentant l'image d'une lame d'un échantillon cytologique (lame). Pour effectuer cette deuxième méthode, l'entrainement du modèle informatique d'apprentissage profond s'effectue à l'aide d'un vaste ensemble de lames labellisées. Sans cet entrainement, la deuxième méthode de classification ne 25 peut fonctionner. Dans ce cas, l'entrée est un ensemble d'imagettes (brutes et/ou caractérisées par la première méthode) et la sortie est une distribution de probabilités selon des classes pour affecter au finale une classe à chaque image de l'ensemble des imagettes de l'échantillon analysée. Selon l'invention, les classes seront choisies parmi une classe de niveau 30 d'anomalie cellulaire faible, une classe de niveau d'anomalie cellulaire fort et/ou une classe de niveau d'anomalie cellulaire nul. En outre, cette méthode effectue un apprentissage « de bout en bout » : c'est-à-dire connectant directement l'entrée à la sortie sans passer par une phase intermédiaire d'analyse de chaque imagette d'une image d'un 5 échantillon cytologique. Plus précisément à partir d'un ensemble d'imagettes de cellules de l'échantillon cytologique, un réseau se voit assigner des tâches à accomplir (une classification, par exemple) et apprend comment les automatiser. Un autre avantage majeur de cette méthode est le fait que les algorithmes de la méthode d'apprentissage évoluent avec les données et 10 continuent à s'améliorer en même temps que le volume des données augmente. De préférence, ces données proviennent au moins en partie des données obtenues grâce à la méthode de détection selon le premier aspect décrit précédemment. De préférence encore, ce deuxième aspect de la méthode de détection pourra avantageusement coupler les données d'une 15 image en transmission et d'une image en fluorescence, c'est-à-dire d'une première et d'une deuxième image. Le deuxième aspect de la méthode de détection propose ici la mis en oeuvre d'une méthode d'apprentissage profond, détournant la fonctionnalité de batch pour pouvoir traiter un ensemble. 20 Classiquement le tenseur 4Dimensions (4D) est composé de 2D pour la composante spatiale de l'image (1D pour la hauteur de l'image et 1D pour la largeur), 1D pour les canaux (RVB si on utilise Rouge Vert Bleu), et 1D pour le batch « B » (B :correspondant au nombre d'images pour construire un batch ou un paquet). 25 De manière générale, l'utilisation du tenseur 4D a comme unique utilité d'accélérer les calculs. L'originalité de notre méthode est d'utiliser ce tenseur pour une autre finalité qui est celle de classifier des populations, ici des ensembles d'imagettes. Contrairement à ce que font classiquement les autres méthodes d'apprentissage qui prennent un lot de B d'images et 30 produisent B décisions : 1 décision pour chaque image (le fait de prendre 36 B>1 sert juste à gagner du temps), nous, on prend un lot de B imagettes (B : correspondant au nombre d'imagettes constituant l'ensemble de l'image) pour donner 1 seule décision basée sur l'ensemble d'imagettes. L'originalité proposée ici est justement d'utiliser cette fonctionnalité de faire des lots qui 5 existe initialement pour des raisons de temps de calcul et d'en détourner sa fonctionnalité en l'utilisant pour un autre but et une autre finalité. Considérons un ensemble d'imagettes correspondant à un nombre B (B devrait être assez petit pour que le batch puisse être transféré au GPU qui a une mémoire limitée). Classiquement, les méthodes d'apprentissages 10 utilisent un batch de 32 à 64 imagettes à choisir pour accélérer le temps de calcul. Dans le contexte de l'invention, on ne choisit pas un nombre B prédéfini puisqu'on utilise ce nombre comme le nombre d'imagettes d'intérêts, par exemple les imagettes de cellules urothéliales par image de lame, ce nombre est différent d'une image de lame à une autre (selon le 15 nombre de cellules urothélilales présente dans l'échantillon cytologique sur la lame). Il est ainsi possible de former un tenseur 4D avec ces B imagettes : B x RVB x H x W (B : nombre d'imagettes, RVB : pour les canaux trois couleurs Rouge Vert Bleu, H et W respectivement pour la hauteur et la largeur des imagettes). Il est également possible de considérer un nombre 20 quelconque de canal classique RVB pour des images de microscopie électronique en transmission mais aussi un nombre quelconque de canal RVB associé à F (F : Fluorescence) pour un couple d'image de microscopie électronique en transmission et de microscopie à fluorescence, tel que RVB (Rouge vert Fluo). Ensuite, si nous transférons ce batch dans le réseau de 25 neurones jusqu'à une couche entièrement connectée (c'est à dire une couche de neurones qui ne tient pas compte de l'aspect spatial de la couche précédente et qui de ce fait produit comme sortie uniquement un vecteur 1D par imagette), la sortie aura donc la taille B x F où F est la taille de cette couche entièrement connectée (ainsi le tenseur 4D a été transformé en un 30 tenseur 2D). Maintenant, il est possible de transposer ce tenseur 2D B x F comme un 37 tenseur 1 xFxB et d'appliquer une opération symétrique sur B par exemple un pool 1D (max ou moyenne). Cette transposition est utile dans la mesure où les opérations symétriques sont nativement programmées pour avoir lieu sur les dimensions 3 et 4 et non sur la dimension 1 - la nécessité d'utiliser 5 une couche entièrement connectée étant justement de libérer ces dimensions 3 et 4 pour la transposition. Ce résultat donne ainsi un tenseur 1 x F, où les caractéristiques F devraient contenir des informations concernant l'ensemble des imagettes. On comprendra que, ce 1 résulte du fait que cette ensemble doit conduire à 1 seule décision de classification affiliée au batch 10 des B imagettes. De cette façon, il est possible d'utiliser une méthode d'apprentissage profond classique pour former un réseau de bout en bout pour une problématique de classification d'une image d'échantillon cytologique c'est-à-dire d'un ensemble d'imagettes. La figure 7 résume la méthodologie d'implémentation pour un apprentissage 15 profond appliqué à la classification d'une image d'échantillon cytologique correspondant à un ensemble d'imagettes. La classe finale est assignée à un ensemble d'imagettes. On comprendra que cette méthode d'implémentation est appliquée ici à la classification d'image d'échantillon cytologique mais peut indifféremment être appliquée à 20 d'autres domaines de façon non limitative. L'invention propose selon un premier aspect une méthode de classification d'échantillon cytologique permettant d'extraire et de calculer des caractéristiques expertes d'anomalie cellulaire présente dans un échantillon pour contribuer au diagnostic du cancer, à l'aide de la vision par 25 ordinateur. L'invention propose, ici, également selon un deuxième aspect, une deuxième méthode de classification d'échantillon cytologique selon des niveaux d'anomalie cellulaire rencontrés dans l'échantillon, en utilisant une méthode d'apprentissage profond, pouvant contribuer au diagnostic du cancer, à l'aide de la vision par ordinateur. L'invention se rapporte également 30 à un dispositif de classification d'échantillon cytologique, à partir d'au moins 38 une image numérisée, adapté à la mise en oeuvre d'une telle méthode. Chaque utilisateur peut définir les seuils qu'il prend en compte pour effectuer la méthode de classification. En outre, on propose une approche générale qui s'adapte à chaque utilisateur selon le besoin c'est-à-dire selon le 5 domaine d'application REVENDICATIONS1. Méthode de détection de cellules présentant au moins une anomalie dans un échantillon cytologique à partir d'au moins une première image (1) de microscopie électronique numérisée de l'échantillon, caractérisée en ce que la méthode est exécutée à l'aide d'un module de traitement mis en oeuvre par ordinateur, la méthode comprend une étape de traitement de la première image (110), comprenant les sous-étapes suivantes : - à partir de la première image (1), détecter chaque élément (4) présent dans l'échantillon (étape 112), puis filtrer chaque élément (4) choisi parmi au moins une cellule isolée (10) ou un groupement de cellules (14), à l'aide d'une détection colorimétrique, - à partir de chaque groupement de cellules (14) détecté, individualiser au moins une cellule avec son noyau, à l'aide d'au moins une méthode de segmentation fondée sur une détection des contours dans le groupement de cellules (14) (étape 116), la méthode comprend en outre une étape de décèlement d'anomalie cellulaire (étape 140), comprenant les sous-étapes suivantes : -pour chaque cellule isolée (10) ou individualisée dans l'étape précédente, calculer des données d'au moins deux caractéristiques d'anomalies distinctes dites « caractéristiques expertes d'anomalie», et déterminer si ladite cellule présente au moins une anomalie (étape 144), en comparant ces données calculées avec des données de références (étape 142), - puis compter le nombre total de cellules dans l'échantillon présentant au moins une anomalie (étape 146), 40 la méthode comprend en outre une première étape de caractérisation de l'échantillon (étape 150) suivant un niveau d'anomalie de l'échantillon, comprenant les sous-étapes suivantes: - comparer le nombre total de cellules présentant au moins une 5 anomalie à une valeur seuil du niveau d'anomalie, de préférence une valeur seuil égale à au moins trois cellules présentant chacune au moins une anomalie, - si le nombre total de cellules présentant au moins une anomalie est inférieur à la valeur seuil du niveau d'anomalie, classer l'échantillon 10 dans une classe représentative d'un échantillon avec un niveau d'anomalie nul (152), - si le nombre total de cellules présentant au moins une anomalie est supérieur à la valeur seuil du niveau d'anomalie, classer l'échantillon dans une classe représentative d'un échantillon avec un niveau 15 d'anomalie avéré (154). 2. Méthode selon la revendication 1, caractérisée en ce qu'au moins une des caractéristiques expertes d'anomalie est un rapport surfacique du noyau dans la cellule ou une irrégularité du noyau ou une colorimétrie mottée ou une hyper colorimétrie du noyau. 20 3. Méthode selon la revendication 1 ou 2, caractérisée en ce qu'au moins une des caractéristiques expertes d'anomalie est une irrégularité du noyau, le calcul de ladite irrégularité du noyau est un calcul d'écart correspondant à la comparaison entre chaque noyau de cellule analysée et un noyau type d'une cellule dite de « référence », ce calcul d'écart est 25 choisi parmi la liste suivante: une différence de surface, une différence de périmètre, une différence de la convexité du noyau, une différence de barycentre, une différence entre une longueur du contour de chaque noyau détecté, une différence d'irrégularité d'un contour du noyau, une différence 41 du nombre de pixels dans chaque noyau et une différence d'excentricité du noyau dans la cellule. 4. Méthode selon l'une quelconque des revendications précédentes, caractérisée en ce qu'elle comprend une étape de sélection typologique 5 (étape 130), comprenant les sous-étapes suivantes: - à partir de la première image (1) numérisée de l'échantillon, mesurer la taille des cellules, puis classer les cellules selon leur taille (étape 132), - associer chaque classe de cellules mesurées à un type cellulaire 10 correspondant, puis sélectionner au moins un type cellulaire choisi, tel que le type cellulaire urothélial (étape 134). 5. Méthode selon la revendication 4, caractérisée en ce qu'elle comprend une seconde étape de caractérisation de l'échantillon (étape 160) suivant un niveau d'anomalie des cellules de l'échantillon, comprenant 15 en outre les sous-étapes suivantes: - à partir d'une seconde image (2) de fluorescence numérisée de l'échantillon, sélectionner toutes les cellules d'un type cellulaire choisi, tel qu'urothélial, selon leur localisation dans l'échantillon prédéfinie (étape 134) au préalable dans la première image (1) numérisée, (étape 20 162), - détecter la présence d'un halo de fluorescence autour de chaque cellule sélectionnée précédemment, de préférence mesurer un taux de fluorescence (étape 164), de préférence encore mesurer un taux d'étalement de la fluorescence, 25 - si un halo de fluorescence est détecté, de préférence si un taux de fluorescence mesuré est supérieur à un taux minimum de fluorescence de référence, de préférence encore si un taux d'étalement de la fluorescence mesuré est supérieur à un taux minimum d'étalement de 42 la fluorescence, classer l'échantillon dans une classe représentative d'un échantillon avec un niveau d'anomalie cellulaire faible (168). 6. Méthode selon l'une quelconque des revendications précédentes, caractérisée en ce qu'elle comprend une étape de recommandation 5 d'exploitabilité (étape 120) de l'échantillon cytologique selon un nombre de cellules comprenant un noyau et selon une valeur seuil minimum d'exploitabilité, de préférence la valeur seuil minimum d'exploitabilité est égale à au moins sept cellules comprenant un noyau détectées dans la première image (1). 10 7. Méthode selon la revendication 5 et selon l'une quelconque des revendications précédentes, caractérisée en ce que la méthode à l'aide du module de traitement combine l'ensemble des données obtenues et classe l'échantillon cytologique analysé à partir de la première image (1) et de la seconde image (2), dans l'une des trois classes associée à un niveau 15 d'anomalie cellulaire de l'échantillon, nul, faible ou fort. 8. Méthode selon l'une quelconque des revendications précédentes, caractérisée en ce que ledit échantillon cytologique est un échantillon urinaire de miction spontanée.