Domaine technique

L'invention concerne le domaine du traitement d'images produites par ordinateur, et, plus précisément, un procédé, un système et un produit-programme d'ordinateur pour déterminer la croissance de micro-organismes afin de détecter et d'identifier des micro-organismes au sein d'un objet tel qu'une boîte de Pétri, après l'incubation dudit objet.

Etat de la technique

Dans le domaine de la microbiologie, il est connu d'utiliser des images d'une boîte de Pétri acquises à différents temps d'incubation afin de détecter et d'identifier d'éventuels micro-organismes. La boîte de Pétri contient un milieu de culture tel qu'une gélose spécifique dans laquelle un échantillon adapté à ladite gélose est ensemencé. De manière connue, une première image de la boîte de Pétri est acquise à un instant TO, préalablement à l'incubation de la boîte de Pétri et au moyen d'un système d'imagerie adapté. Ainsi, la première image acquise ne comprend aucune colonie de micro-organismes. La boîte de Pétri est ensuite incubée pendant une période de temps déterminée, jusqu'à un instant Tl. Une deuxième image de la boîte de Pétri est alors acquise à l'instant Tl. Ensuite, la première image acquise au temps TO est comparée avec l'image acquise au temps Tl, afin de déterminer si des colonies de micro-organismes sont apparues lors de la période d'incubation, c'est-à-dire si l'échantillon ensemencé a généré la croissance de colonies de micro-organismes au sein de la boîte de Pétri. Cependant, la comparaison des images de la boîte de Pétri ne permet pas de détecter avec fiabilité une éventuelle présence de colonies de micro-organismes. En effet, de nombreux paramètres peuvent engendrer un résultat erroné de la comparaison des images de la boîte de Pétri.

Concernant la boîte de Pétri elle-même, différents paramètres peuvent engendrer une détection erronée des colonies de micro-organismes.

Un premier paramètre concerne le matériau de la boîte de Pétri. Ainsi, la boîte de Pétri est réalisée dans un matériau, tel que le verre ou le plastique, qui peut comporter des défauts tels que des rayures. Selon les conditions d'acquisition des images de la boîte de Pétri, les rayures peuvent apparaître sur l'image de la boîte de Pétri en tant qu'élément distinct. Ainsi une rayure peut être interprétée, de manière erronée, comme représentative d'une croissance de micro-organismes.

Un deuxième paramètre concerne des éléments associés à la boîte de Pétri. Ainsi, la boîte de Pétri comprend généralement une inscription sérigraphiée, par exemple un numéro permettant d'identifier l'origine de la production de la boîte de Pétri. De manière similaire aux rayures, selon les conditions d'acquisition des images de la boîte de Pétri, l'inscription sérigraphiée peut apparaître sur l'image de la boîte de Pétri en tant qu'élément distinct. Ainsi une inscription sérigraphiée peut être interprétée, de manière erronée, comme représentative d'une croissance de microorganismes.

Un troisième paramètre concerne le milieu de culture disposé au sein de la boîte de Pétri, tel que la gélose. En effet, la gélose nécessaire à la mise en culture de l'échantillon peut contenir des défauts structurels. Selon les conditions d'acquisition des images de la boîte de Pétri, ces défauts structurels tels que des bulles ou une rétractation de la gélose peuvent également apparaître en tant qu'éléments distincts sur l'image de la boîte de Pétri. Ainsi, les défauts structurels peuvent également être interprétés, de manière erronée, comme représentatifs d'une croissance de micro - organismes. Un quatrième paramètre concerne la présence de buée. En effet, lors de la période d'incubation, des particules d'eau formant de la buée peuvent apparaître sur les parois de la boîte de Pétri. Or, lors de l'analyse de la boîte de Pétri, celle-ci subit un retournement qui peut provoquer la chute des particules d'eau sur la gélose. Ainsi, selon les conditions d'acquisition d'une image de la boîte de Pétri, lesdites particules d'eau peuvent apparaître sous forme d'éléments distincts sur ladite image de la boîte de Pétri. Ainsi, la présence de buée peut également engendrer une interprétation erronée d'une croissance de micro-organismes. Enfin, un cinquième paramètre concerne d'éventuels déplacements de la boîte de Pétri en translation et/ou en rotation entre l'acquisition de la première image à T0 et l'acquisition de la deuxième image à Tl. Ces déplacements génèrent un décalage de la position de l'objet au sein de la première image acquise par rapport à la position dudit objet au sein de la deuxième image acquise. Or, un décalage entre les images acquises aux instant T0 et Tl ne permet pas de les comparer par superposition, en vue de détecter la présence de colonies de micro-organismes.

Dans l'art antérieur, le brevet EP0796319 décrit un procédé d'identification de colonies de micro-organismes comprenant l'acquisition d'images d'un dispositif de mise en culture tel qu'une boîte de Pétri ou un film dit Petrilfilm™. Le procédé comprend des étapes de recalage des images acquises en utilisant comme position de référence le centre du dispositif de mise en culture. Le recalage concerne uniquement un mouvement de translation du dispositif de mise en culture entre une première image acquise à un temps T0 et une deuxième image acquise à un temps Tl. Par conséquent, le brevet EP0796319 ne décrit pas d'étape de recalage permettant de considérer un éventuel mouvement de rotation du dispositif de mise en culture entre les deux temps d'acquisition T0 et Tl des images. De plus, les conditions d'acquisition des images sont telles que les conditions d'illumination d'une image acquise au temps T0 sont identiques à celles d'une image acquise au temps Tl. Ainsi, le brevet EP0796319 décrit des étapes de comparaison des images comprenant la soustraction des pixels de l'image acquise à T0 avec les pixels de l'image acquise à Tl afin d'obtenir une image de comparaison comprenant des pixels résultants dont l'intensité est éventuellement représentative d'une croissance de micro-organismes. Par conséquent, le brevet EP0796319 ne décrit pas d'étapes de comparaison d'images adaptées à des images acquises selon différentes conditions d'illumination.

Toujours dans l'art antérieur, le brevet US 5,510,246 divulgue un procédé pour déterminer la croissance de micro-organismes au sein d'un/sur un dispositif de culture substantiellement plan tel qu'une plaque de type Petrifilm™. Ce procédé comprend notamment des étapes préliminaires de filtrage d'images, plus précisément de filtrage de pixels, afin d'obtenir une image d'arrière-plan dépourvue de pixels « parasites ». Cette image d'arrière-plan est ensuite utilisée pour identifier les différences observées entre une image du dispositif de culture prise après une certaine durée d'incubation du milieu de culture (préalablement inoculé) et ladite image d'arrière-plan. Même si le brevet US 5,510,246 traite, de manière générale, d'un procédé de traitement d'images en vue de dénombrer les colonies microbiennes ayant poussé sur un milieu de culture présent au sein/sur un dispositif de culture, il est important de noter que US 5,510,246 enseigne que le procédé divulgué dans ce brevet est automatisé et considéré comme fiable en termes de positionnement du dispositif de culture, ce qui rend totalement inutile tout étape de recalage d'image, et qui plus est tout étape de recalage nécessitée par un mouvement de type rotation du dispositif de culture au cours du procédé d'acquisition d'images.

Dans l'art antérieur, la demande de brevet internationale WO 2014/059357 décrit également un procédé pour déterminer des paramètres de croissance de cellules ensemencées sur un support d'analyse tel qu'une plaque de verre. Le procédé comprend une étape de recalage des images du support d'analyse, laquelle comprend notamment l'utilisation :

d'un algorithme connu dans l'état de l'art, l'algorithme « Starmatch », qui utilise des matrices de transformation affine ainsi que le regroupement des objets recalés dans un espace euclidien à deux dimensions, afin de déterminer la liste des points d'intérêt dans chaque image puis d'en déduire la transformation optimale, et d'une fonction également connue, la fonction "Gompertz", appliquée aux zones de données des objets regroupés pour déterminer les zones de croissances de micro-organismes.

Un inconvénient majeur du procédé objet de WO 2014/059357 réside dans la présence obligatoire de points d'intérêt dans une image acquise à un temps initial afin d'être en mesure de retrouver ces points d'intérêt dans une image acquise à un temps ultérieur. Par conséquent, le procédé objet de WO 2014/059357 ne permet pas de recaler des images en l'absence de colonies de micro-organismes à un temps initial. La demande de brevet européen EP-A-1061127 divulgue, quant à elle, un procédé pour identifier des micro-organismes présents au sein d'un récipient d'analyse tel qu'une boîte de Pétri. Le procédé comprend la prise de deux images, l'une après l'autre, dans l'optique de mieux déterminer les contours et les détails d'au moins une colonie de micro-organismes. Plus précisément, la première image est acquise avec une première position de la caméra et la deuxième image est acquise avec une deuxième position de la caméra. Le déplacement entre la première et la deuxième image est donc un déplacement connu (et donc prévisible), aisément identifiable sur les images. EP-A-1061127 ne divulgue, pas plus que ne suggère, d'étape de recalage d'images lors d'un déplacement imprévu/inopiné du récipient d'analyse.

En outre, l'article de DEN HERTOG, A. et al., « Simplified Automated Image Analysis for Détection and Phenotyping of Mycobacte um tubercolisis on Porous Supports by Monitoring Growing Microcolonies ». PLos ONE, (2010). Volume 5, Issue 6, ell008, décrit une méthode d'analyse de la croissance de micro-organismes au sein d'un support d'oxyde d'aluminium poreux sous forme de bandes. Le procédé permet une détection automatique des colonies de micro-organismes tout en détectant de manière spécifique les micro-colonies individuelles afin de les surveiller au cours du temps et de les identifier en fonction de leur croissance. L'article ne décrit aucune étape de recalage relative à un déplacement du support d'analyse.

Une autre publication, à savoir CHEN W-B. « An automated bacterial colony counting and classification System ». Sp nger Science + Business Media, (2009). Published online, Inf Syst Front, 11:349-368, décrit un procédé de comptage automatique de colonies de micro-organismes au sein d'un support d'analyse qui peut être une boîte de Pétri. Le procédé comprend la détection du support d'analyse, l'identification des colonies, la séparation des colonies agrégées et la détermination du nombre de colonies. Cette publication ne divulgue aucune étape de recalage basée sur un déplacement du support d'analyse.

Par conséquent, il s'avère nécessaire d'améliorer la détection de croissance de micro-organismes au sein d'un objet tel qu'une boîte de Pétri en considérant le recalage d'image lié à une translation et/ou à une rotation - pouvant se produire de manière inopinée - de la boîte de Pétri lors de l'acquisition des images de la boîte de Pétri à des temps différents.

Objet de l'invention

La présente invention vise à surmonter au moins partiellement les problèmes mentionnés ci-dessus.

Ainsi, un premier objectif de l'invention consiste à fournir un procédé pour déterminer la croissance de micro-organismes dans un échantillon biologique susceptible de contenir des micro-organismes, ledit échantillon biologique étant contenu dans un récipient d'analyse tel qu'une boîte de Pétri, ledit récipient d'analyse étant soumis à une incubation d'une durée déterminée, ledit procédé comprenant les étapes suivantes :

- acquérir une première pluralité d'images initiales du récipient d'analyse à un premier temps d'acquisition Tl, avant ou pendant l'incubation ;

- acquérir une deuxième pluralité d'images finales du récipient d'analyse à un deuxième temps d'acquisition T2, pendant ou après l'incubation ;

- recaler chaque image initiale de la première pluralité d'images initiales acquises, avec chaque image finale correspondante de la deuxième pluralité d'images finales acquises ;

- localiser au moins une zone potentielle de croissance de micro-organismes au sein d'au moins une image de la deuxième pluralité d'images acquises ; - évaluer le contenu de la zone potentielle de croissance de microorganismes identifiée afin de déterminer la présence de micro-organismes.

De manière avantageuse, l'étape de recalage comprend une étape de recalage primaire (également dénommé « recalage grossier ») associée à un identifiant localisé sur le récipient d'analyse.

De manière avantageuse, l'étape de recalage comprend une étape de recalage secondaire associée au contenu du récipient d'analyse.

Comme cela est connu de l'homme du métier, le recalage de deux images consiste à mettre en correspondance le ou les objets réels visualisés sur chacune de ces deux images. Le recalage de deux images s'avère particulièrement utile lorsqu'un déplacement du dispositif de prise d'images et/ou du ou des objets d'intérêt se produit entre l'acquisition des deux images.

De manière avantageuse, l'étape de localisation d'au moins une zone potentielle de croissance de micro-organismes comprend une étape de détection d'une densité élevée de colonies de micro-organismes, une étape de détection de colonies de micro-organismes non-circulaires et une étape d'attribution d'un élément de localisation à ladite zone potentielle de croissance de micro-organismes.

De manière avantageuse, l'étape d'évaluation du contenu de la zone potentielle de croissance identifiée comprend une étape de détermination des valeurs de deux paramètres d'évolution.

De manière avantageuse, les deux paramètres d'évolution comprennent un paramètre de corrélation et un paramètre de contraste. De manière avantageuse, selon l'invention l'étape de recalage secondaire (également dénommé « recalage fin ») comprend la création d'un maillage de sous- images pour chaque image initiale et chaque image finale. Comme cela est connu de l'homme du métier, une « sous-image » s'entend d'une partie d'une image, de taille inférieure à celle de ladite image. Une « sous-image » permet de focaliser/centrer l'image sur un point d'intérêt ou un aspect local de l'image. La taille et/ou la position de la sous-image (au sein de l'image) est/sont variable(s) en fonction des souhaits de l'opérateur.

Le recalage secondaire permet de compenser, de manière précise, à savoir à l'échelle du pixel, le décalage (de type translation et/ou rotation) créé entre les images acquises à Tl et les images acquises à T2. Ainsi, les étapes ultérieures de localisation d'au moins une zone potentielle de croissance de micro-organismes et d'évaluation du contenu de cette zone potentielle de croissance de micro - organismes peuvent être réalisées de manière optimale, tout en limitant - voire supprimant - l'impact des artefacts (tels que la présence d'une rayure ou d'une inscription sérigraphiée sur le récipient d'analyse, par exemple sur une boîte de Pétri) sur l'interprétation des résultats.

Ainsi, un deuxième objectif de l'invention consiste à fournir un système pour déterminer la croissance de micro-organismes dans un échantillon biologique susceptible de contenir des micro-organismes, ledit échantillon biologique étant contenu dans un récipient d'analyse telle qu'une boîte de Pétri, ledit récipient d'analyse étant soumis à une incubation d'une durée déterminée, ledit système comprenant :

- un dispositif de capture d'images afin d'acquérir une pluralité d'images de l'objet à analyser à un premier temps d'acquisition Tl, avant ou pendant l'incubation, et à un deuxième temps d'acquisition T2, pendant ou après l'incubation ;

- un dispositif d'illumination comprenant une ou plusieurs sources lumineuses afin d'illuminer le récipient d'analyse ;

- un dispositif de contrôle pour contrôler l'application du procédé de détermination de la croissance de micro-organismes selon l'invention afin de déterminer la présence de micro-organismes. Ainsi, un troisième objectif de l'invention consiste à fournir un produit-programme d'ordinateur comprenant des instructions logicielles pour la mise en œuvre d'un procédé selon l'invention lorsque ledit programme est exécuté par un processeur de données.

Brève description des dessins

L'invention, sa fonctionnalité, ses applications ainsi que ses avantages seront mieux appréhendés à la lecture de la présente description, faite en référence aux figures, dans lesquelles :

- la figure 1 montre une vue schématique d'un système d'analyse d'images selon un mode de réalisation de l'invention ;

- la figure 2 montre un diagramme des étapes du procédé global de détermination de la croissance de micro-organismes selon un mode de réalisation de l'invention ;

- la figure 3 montre un diagramme des étapes du procédé de recalage primaire selon un mode de réalisation de l'invention ;

- la figure 4 montre une image d'une boîte de Pétri acquise au temps Tl ;

- la figure 5 montre une image de la boîte de Pétri selon la figure 4 avec une représentation ciblée sur les pixels distants du centre de la boîte de Pétri selon un rayon R-10 pixels ;

- la figure 6 montre une image de la boîte de Pétri selon la figure 4 avec une représentation ciblée sur les pixels distants du centre de la boîte de Pétri selon un Rayon R-50 pixels ;

- la figure 7 montre une image de la boîte de Pétri selon la figure 4 en application d'une opération de segmentation ;

- les figures 8 et 9 montrent une image de la boîte de Pétri selon la figure 4 après l'application d'une opération d'ouverture morphologique ;

- la figure 10 montre une image de la boîte de Pétri sur la figure 4 après application d'une fonction d'érosion ;

- la figure 11 montre en détail le contenu du rectangle représenté sur la figure 10 ; - la figure 12 montre une image de la boîte de Pétri selon la figure 9 après application d'une fonction d'ouverture ;

- la figure 13 montre en détail le contenu du rectangle représenté sur la figure 12 ;

- la figure 14 montre une représentation graphique de l'étiquette latérale théorique ;

- la figure 15 montre une représentation graphique de l'étiquette latérale théorique en traits discontinus, de l'angle d'orientation Θ en trait plein et de l'ellipsoïde équivalent sous forme de points ;

- la figure 16 montre une représentation graphique comparative des images de l'étiquette latérale acquise à Tl et à T2 ;

- la figure 17 montre une représentation graphique selon la figure 16 basée sur les moitiés supérieures desdites étiquettes latérales théoriques ;

- la figure 18 montre une représentation graphique des étiquettes selon la figure 16 basée sur les moitiés inférieures desdites étiquettes latérales théoriques ;

- la figure 19 montre une image d'une boîte de Pétri acquise à un temps Tl ;

- la figure 20 montre une image de la boîte de Pétri acquise à un temps T2 ;

- la figure 21 montre un diagramme des étapes concernant le procédé de recalage secondaire selon le mode de réalisation de l'invention ;

- la figure 22 montre la représentation d'un pixel au sein d'une cellule initiale

- la figure 23 montre en détail le contenu du carré représenté sur la figure 22

- la figure 24 montre la représentation d'un pixel au sein d'une cellule finale ; - la figure 25 montre en détail le contenu du carré représenté sur la figure 24

- la figure 26 montre un diagramme des étapes relatives au procédé de recalage local selon la figure 21 ;

- la figure 27 montre une image de la boîte de Pétri acquise à un temps Tl selon la condition d'acquisition « backlight » ;

- la figure 28 montre une image de la boîte de Pétri selon la figure 27 acquise à un temps T2 selon la condition d'acquisition « backlight » ; - la figure 29 montre en détail le contenu du carré représenté sur la figure

28 ;

- la figure 30 montre un diagramme des étapes relatives au procédé de détermination de croissance préalable de micro-organismes selon la figure 2 ;

- la figure 31 montre un diagramme des étapes du procédé d'analyse temporel ;

- la figure 32 montre un diagramme des étapes relatives au procédé de prétraitement d'images selon la figure 31 ;

- la figure 33 montre des exemples de sous-images acquise à Tl et à T2 selon les conditions d'acquisition « backlight », « bottom » et « médian » après l'application d'un procédé d'ajustement de contraste ;

- la figure 34 montre des exemples de sous-images acquise à Tl et à T2 selon les conditions d'acquisition « backlight », « bottom » et « médian » après l'application d'un filtre médian ;

- la figure 35 montre un diagramme des étapes relatives à l'évaluation des premier et deuxième paramètres d'évolution selon la figure 31 ;

- la figure 36 montre un diagramme des étapes relatives à l'évaluation du paramètre de contraste selon un mode de réalisation de l'invention ;

- la figure 37 montre un diagramme des étapes relatives à l'évaluation du paramètre de corrélation selon un mode de réalisation de l'invention ;

- la figure 38 montre un diagramme des étapes relatives à la détermination de la croissance effective selon la figure 31 selon un mode de réalisation de l'invention ;

- la figure 39 montre une représentation graphique de la croissance de micro-organismes telle que détectée, en fonction des valeurs des paramètres de contraste et de corrélation.

Description détaillée de l'invention

La description détaillée ci-après a pour but d'exposer l'invention de manière suffisamment claire et complète, notamment à l'aide d'exemples, mais ne doit en aucun cas être regardée comme limitant l'étendue de la protection aux modes de réalisation particuliers et aux exemples présentés ci-après.

La présente invention concerne l'analyse d'un échantillon. Selon la présente invention, l'échantillon peut être de diverses origines, par exemple d'origine alimentaire, environnementale, vétérinaire, clinique, pharmaceutique ou cosmétique.

Parmi les échantillons d'origine alimentaire, l'on peut citer, de façon non-exhaustive, un échantillon de produits lactés (yaourts, fromages, etc.), de viande, de poisson, d'œuf, de fruit, de légume, d'eau, de boisson (lait, jus de fruits, soda, etc.). Bien évidemment, ces échantillons d'origine alimentaire peuvent aussi provenir de sauces ou de plats plus élaborés ou des matières premières non transformées ou partiellement transformées. Un échantillon alimentaire peut également être issu d'une alimentation destinée aux animaux, telle que des tourteaux, des farines animales.

Tel qu'indiqué précédemment, l'échantillon peut être d'origine environnementale et peut consister, par exemple, en un prélèvement de surface, d'eau, d'air, etc.

L'échantillon peut également consister en un échantillon d'origine clinique, pouvant correspondre à des prélèvements de fluide biologique (urine, sang total, ou dérivés tel que sérum, salive, pus, liquide céphalo-rachidien, etc.), de selles (par exemple diarrhées cholériques), de prélèvements de nez, de gorge, de peau, de plaies, d'organes, de tissus ou de cellules isolées. Cette liste n'est évidemment pas exhaustive.

D'une manière générale, le terme « échantillon » se réfère à une partie ou à une quantité, plus particulièrement une petite partie ou une petite quantité, prélevée à partir d'une ou plusieurs entités aux fins d'analyse. Cet échantillon peut éventuellement avoir subi un traitement préalable, impliquant par exemple des étapes de mélange, de dilution ou encore de broyage, en particulier si l'entité de départ est à l'état solide.

L'échantillon prélevé est, en général, susceptible de - ou suspecté de - contenir au moins un micro-organisme externe cible tel qu'une bactérie, une levure, un champignon.

Selon la présente invention, l'analyse est réalisée à l'aide d'un récipient d'analyse tel qu'une boîte de Pétri comprenant un réceptacle et un couvercle. Le réceptacle et/ou le couvercle peuvent comprendre un identifiant. Le réceptacle comprend un milieu de culture tel qu'un milieu de culture gélosé. Un fluide est ensemencé sur le milieu de culture avant le début de l'analyse. Le but de l'analyse est de déterminer la présence de micro-organismes spécifiques et la concentration en micro-organismes au sein du fluide sur le milieu de culture, après incubation de la boîte de Pétri.

Selon un mode de réalisation préféré de l'invention, le récipient d'analyse est une boîte de Pétri qui comprend comme identifiant une étiquette latérale disposée sur la paroi latérale du réceptacle de la boîte de Pétri. Selon un mode de réalisation de l'invention, les conditions d'acquisition des images d'une boîte de Pétri sont identiques. Dans ce mode de réalisation, l'éclairage et l'angle de vue de la boîte de Pétri, lors de la capture d'images correspondantes à deux temps d'acquisition distincts, sont identiques pour chaque image acquise. Selon la présente invention, l'objectif de l'analyse d'images d'une boîte de Pétri est de recaler une image acquise avant incubation de la boîte de Pétri avec une image correspondante acquise après incubation de la boîte de Pétri, ou deux images de la boîte de Pétri acquises à des temps d'incubation différents. En effet, lors de la manipulation de la boîte de Pétri, à l'aide d'un dispositif automatisé ou manuellement, la position de la boîte de Pétri peut varier entre l'instant d'acquisition de l'image avant l'incubation et l'instant d'acquisition de l'image correspondante après incubation. Ainsi, la comparaison des deux images acquises à des instants différents permet de déterminer la présence ou l'absence de colonies de micro-organismes localisée au sein de la boîte de Pétri. La boîte de Pétri étant circulaire, le recalage peut nécessiter une rotation et/ou une translation de l'image de la boîte de Pétri afin d'obtenir la correspondance entre les deux images. La figure 1 montre un système d'analyse d'images 10 comprenant un dispositif de capture d'image 12, un dispositif d'illumination 14, un dispositif d'analyse d'images 16 et un dispositif de contrôle 18.

Le dispositif de capture d'images 12 comprend un dispositif numérique tel qu'une caméra numérique afin d'acquérir des images de l'objet à analyser. Le dispositif de capture d'images 12 acquiert une première série d'images à un temps Tl et une deuxième série d'images à un temps T2. Le temps Tl correspond à un instant initial c'est-à-dire avant l'ensemencement de micro-organismes et avant translation, ou à un instant initial en présence de micro-organismes et après une première période d'incubation au sein de la boîte de Pétri. Le temps T2 correspond à un instant ultérieur à l'instant initial Tl, après incubation de la boîte de Pétri, ou après une deuxième période d'incubation, c'est-à-dire, potentiellement, en présence de microorganismes au sein de la boîte de Pétri. A titre d'exemple, le dispositif de capture d'images 12 comprend un capteur d'images monochrome de type CCD consommant entre 12 et 24 Volts et utilisant la technologie de balayage avec une vitesse maximale d'acquisition de 17 images par seconde.

Le dispositif d'illumination 14 comprend par exemple une ou plusieurs sources lumineuses telles que des diodes électroluminescentes (DELs) afin d'illuminer l'objet à analyser selon des conditions déterminées par l'utilisateur. Les sources lumineuses sont disposées soit au-dessus, soit au-dessous de la boîte de Pétri à analyser. Le dispositif d'illumination 14 permet de générer des conditions d'acquisition relatives à différentes possibilités d'illumination de la boîte de Pétri telles qu'une illumination rouge, verte et/ou bleue (RVB) combinées à la présence ou à l'absence d'un cache situé au-dessous de la boîte de Pétri, c'est-à-dire au- dessous du réceptacle ou au-dessous du couvercle de la boîte de Pétri. Ainsi, l'illumination de type « backlight » correspond à l'acquisition d'images en présence d'une illumination arrière de la boîte de Pétri, en l'absence d'un cache situé au- dessous de la boîte de Pétri. L'illumination de type « bottom annular» ou «bottom» correspond à l'acquisition d'images en présence d'une illumination située au-dessus de la boîte de Pétri, et selon quatre directions par rapport à la boîte de Pétri, c'est-à- dire vers la gauche, vers la droite, vers la partie supérieure, et vers la partie inférieure de ladite boîte de Pétri. L'illumination de type « left bottom » correspond à l'obtention d'images basée sur une illumination située au-dessus de la boîte de Pétri de type « bottom » et du côté gauche de ladite boîte de Pétri. L'illumination de type « médian » correspond à l'obtention d'une image médiane basée sur les quatre images distinctes relatives à chaque direction, gauche, droite, partie inférieure, et partie supérieure de la boîte de Pétri, obtenues au moyen d'une illumination « bottom». D'autres conditions d'acquisition peuvent être envisagées en utilisant une combinaison de sources lumineuses disposées selon une direction spécifique par rapport à la boîte de Pétri en présence ou en absence d'un cache situé au- dessous de la boîte de Pétri afin d'obtenir un objet visible et sans artefact en fonction de la nature de l'échantillon et des micro-organismes à identifier.

Le dispositif d'illumination 14 peut également produire des illuminations verticales ou annulaires par rapport à la boîte de Pétri à analyser.

Le dispositif d'analyse d'images 16 comprend notamment un support d'objet et des moyens de maintien de l'objet tel qu'une pince afin de maintenir un objet tel qu'une boîte de Pétri. Le dispositif de contrôle 18 comprend notamment un micro-processeur (non montré), un écran d'affichage (non montré), une mémoire de stockage (non montrée), afin d'exécuter des algorithmes de traitement d'images. Le dispositif de contrôle 18 permet de définir les paramètres de fonctionnement du dispositif de capture d'images 12, du dispositif d'illumination 14 et du dispositif d'analyse d'images 16 afin de contrôler le fonctionnement du dispositif de capture d'images 12, du dispositif d'illumination 14 et du dispositif d'analyse d'images 16. Le système d'analyse d'images 10 fonctionne en association avec un dispositif d'incubation 20. Le dispositif d'incubation 20 permet d'incuber au moins une boîte de Pétri selon des conditions d'incubation spécifiques concernant la durée et la température, afin de favoriser la croissance de micro-organismes au sein des boîtes de Pétri.

Ainsi, selon un fonctionnement classique, et en considérant une seule boîte de Pétri à titre d'exemple, la boîte de Pétri ensemencée avec un échantillon est acheminée via un moyen de transport connu selon l'état de l'art, tel qu'un tapis roulant, au sein du dispositif d'incubation 20.

A l'issue de la période d'incubation, la boîte de Pétri est à nouveau acheminée au sein du système d'analyse d'images 10. Après le traitement des images obtenues au moyen du système d'analyse d'images 10, un système d'identification (non montré) permet de déterminer, selon un procédé manuel ou automatisé appliqué aux images acquises, la nature et le nombre de colonies de micro-organismes qui se sont éventuellement développées au sein de la boîte de Pétri.

La présente invention concerne l'analyse d'images disponibles au sein du système d'analyse d'images 10, après incubation d'une boîte de Pétri et acquisition des images de ladite boîte de Pétri. Au sein de la présente description, à titre d'exemple, l'objet considéré est une boîte de Pétri. La présente invention s'applique à une pluralité d'images ou séries d'images concernant une ou plusieurs boîtes de Pétri. Une première série d'images Il est acquise à un temps T=T1 et une deuxième série d'images 12 est acquise à un temps T=T2. La présente description, à titre d'exemple, fait référence à une image II issue de la première série d'images II acquise à Tl et à une image 12 issue de la deuxième série d'images 12 acquise à T2. Comme montré sur la figure 2, le procédé global selon la présente invention comprend une première étape d'acquisition 22 de la première série d'images II, à un temps Tl, une deuxième étape d'acquisition 24 de la deuxième série d'images 12, à un temps T2, ainsi qu'une étape 26 relative à un premier procédé de recalage primaire, une étape 28 relative à un deuxième procédé de recalage secondaire, une étape 30 relative à un troisième procédé de détection préalable de microorganismes. Dans l'étape 22, les images sont acquises selon un procédé décrit au sein de la demande internationale publiée sous le numéro WO2012/152769. Les premier, deuxième et troisième procédés sont distincts et leurs algorithmes respectifs sont stockés au sein de la mémoire de stockage du dispositif de contrôle 18.

Le premier procédé concerne le recalage primaire ou recalage grossier montré sur la figure 3 et relatif au positionnement de la boîte de Pétri entre la première image II acquise à Tl et la deuxième image 12 acquise à T2.

Le procédé de recalage primaire est basé sur la détection de l'étiquette latérale de la boîte de Pétri. Ainsi, le procédé de recalage primaire comprend, comme montré sur la figure 3, un procédé de détection de l'étiquette latérale 1001 et un procédé de détermination des paramètres de recalage primaire 1002 associés aux caractéristiques de l'étiquette latérale.

Le procédé de détection de l'étiquette latérale 1001 comprend l'utilisation d'au moins deux images de la boîte de Pétri, avec au moins une image acquise en utilisant la condition d'illumination « left bottom » et au moins une image acquise en utilisant la condition d'illumination « backlight ».

Le procédé de détection de l'étiquette latérale 1001 comprend une première étape 100 concernant la détection de la boîte de Pétri. La première étape 100 comprend une détection des bords de la boîte de Pétri, selon une méthode connue dans l'état de l'art. Les images utilisées dans la première étape correspondent aux images acquises au temps Tl et au temps T2 avec la condition d'illumination « backlight ». Ainsi, la première étape permet de déterminer les coordonnées du centre C et la longueur du rayon R de la boîte de Pétri, et ce pour chaque image acquise.

Dans les étapes suivantes, seuls les pixels situés à l'intérieur d'un disque caractérisé par les valeurs ainsi déterminées de R et C sont considérés au sein des images acquises. Le procédé de détection de l'étiquette latérale comprend une deuxième étape 110 concernant la localisation de l'étiquette latérale.

Les images utilisées dans la deuxième étape 110 correspondent aux images acquises au temps Tl comme montré sur la figure 4, et au temps T2 avec la condition d'illumination « left bottom». En effet l'étiquette latérale localisée sur le côté droit de la boîte de Pétri est visible, de manière optimale, selon cette condition d'illumination.

La deuxième étape 110 comprend la création d'un premier disque de rayon RI = R- dl et la création d'un deuxième disque de rayon R2=R-d2, dl et d2 étant des valeurs déterminées au préalable par l'utilisateur, avec dl<d2, telles que dl=10 pixels et d2=50 pixels. Pour une boîte de Pétri de rayon R=950 pixels, le premier disque, comme montré sur la figure 5, correspond à un rayon Rl=940 pixels et le deuxième disque, correspond à un rayon R2=900 pixels. La deuxième étape comprend la création d'un anneau résultant contenant les pixels situés à une distance d comprise entre RI et R2, soit entre 900 et 940 pixels, comme montré sur la figure 6.

Dans les étapes suivantes, seuls les pixels situés à l'intérieur de l'anneau résultant sont considérés au sein des images acquises.

La deuxième étape 110 comprend ensuite l'application d'une opération de segmentation ou binarisation sur l'image acquise de la boîte de Pétri. En effet, dans la mesure où l'étiquette latérale se distingue par une portion claire d'anneau du reste de l'anneau résultant, l'opération de segmentation permet d'isoler efficacement l'étiquette latérale. L'opération de segmentation peut être utilisée en appliquant un algorithme tel que l'algorithme des k-moyennes ou « k-means » du logiciel Matlab™ basé sur une recherche de deux groupes ou « deux clusters » en considérant, pour un cluster, un centroïde initialisé à la valeur minimale du triplet RVB pour approcher la valeur du blanc et regrouper ainsi les pixels clairs et, pour un autre cluster, un centroïde initialisé à la valeur maximale du triplet RVB pour approcher la valeur du noir, et regrouper ainsi les pixels foncés.

L'opération de segmentation permet de repérer différents objets au sein des images acquises, comme montré sur la figure 7. Ces objets comprennent par exemple l'étiquette latérale et les défauts de réflexion. Les défauts de réflexion sont générés par la réflexion des faisceaux lumineux du dispositif d'illumination 14 sur les bords de la boîte de Pétri. Les défauts de lumière ou « barbules » ont une forme généralement étroite et incurvée.

L'étiquette latérale comprend éventuellement des défauts d'illumination relatifs à une surexposition lors de l'illumination de la boîte de Pétri.

De manière générale, l'étiquette latérale est formée d'une seule pièce. Cependant, l'étiquette latérale peut être parfois constituée de différents fragments, de manière volontaire ou de manière involontaire. Ainsi, selon le type de conditionnement des boîtes de Pétri ou en fonction d'incidents liés par exemple à la qualité du transport des boîtes de Pétri, un fragment de l'étiquette latérale peut s'avérer manquant. Une étape de suppression des défauts d'illumination ne permet pas toujours de distinguer l'étiquette latérale des défauts d'illumination lorsque la taille de l'étiquette latérale est du même ordre de grandeur que la taille des défauts d'illumination. Ainsi, le procédé de détection de l'étiquette latérale comprend une troisième étape 120 concernant la suppression des défauts d'illumination de l'étiquette latérale. La troisième étape 120 comprend l'application d'une fonction d'ouverture morphologique, selon l'art antérieur, sur l'image binarisée, montrées sur les figures 8 et 9. La fonction d'ouverture correspond à la combinaison d'une fonction d'érosion et d'une fonction de dilatation. Pour l'application de la fonction d'ouverture, la fonction STREL du logiciel Matlab™ est utilisée afin de définir un élément structurant spécifique comprenant un segment de longueur déterminée, par exemple 10 pixels, concernant un nombre défini de pixels voisins du pixel considéré, par exemple 3 pixels, associé à un angle de rotation dudit élément structurant compris entre -90 et 90 degrés. Comme montré sur les figures 10 et 11, la valeur de l'angle de rotation est 0°. Comme montré sur les figures 12 et 13, la valeur de l'angle de rotation est 45°.

La fonction d'érosion, avec les valeurs définies ci-dessus, permet d'éliminer les objets ayant une épaisseur inférieure à 10 pixels quelle que soit la courbure de ces objets.

La fonction de dilatation permet d'obtenir une correction de la courbure de l'étiquette latérale, afin d'obtenir une courbure globale de forme convexe.

La fonction d'ouverture est appliquée pour chaque canal de couleur du triplet RVB.

Ainsi, à l'issue de l'application de la fonction d'ouverture, les défauts d'illumination de l'étiquette latérale sont supprimés sur les images acquises au temps Tl et au temps T2. Ensuite, le procédé de détection de l'étiquette latérale comprend la quatrième étape 130 relative à la détermination d'une étiquette latérale théorique formée d'un seul fragment. La zone de la boîte de Pétri associée à l'étiquette latérale théorique correspond à un arc de cercle de sommet situé au centre de la boîte de Pétri, de largeur 20 pixels et avec un angle d'étalement Θ.

L'angle d'étalement Θ de l'étiquette latérale théorique est déterminé par la valeur maximale 9max et la valeur minimale 9min de Θ afin de définir la surface d'étalement recouvrant l'ensemble des fragments de l'étiquette latérale détectés à l'étape précédente.

La quatrième étape 130 permet d'obtenir une étiquette latérale théorique, comme montré sur la figure 14. L'étiquette latérale théorique est située sur le pourtour de la boîte de Pétri, comme montré sur la figure 5, et contient tous les fragments de l'étiquette latérale.

Ainsi, le procédé de détermination des paramètres de recalage primaire 1002 peut être appliqué comme décrit ci-dessous.

Le procédé de détermination des paramètres de recalage primaire 1002 comprend une première étape 140 concernant la détermination des orientations des étiquettes latérales et la détermination d'un intervalle de valeurs des angles de rotation.

Pour définir une plage de recherche pour l'angle de rotation, on utilise l'orientation des étiquettes latérales théoriques au temps Tl et au temps T2.

De manière connue, la fonction « regionprops » associée à la propriété « Orientation » du logiciel Matlab™ permet de déterminer l'orientation d'objets au sein d'une image binarisée. L'orientation de l'étiquette théorique latérale est déterminée en obtenant les angles Θ1 et Θ2. Cette orientation correspond à l'orientation de l'ellipse ayant le même moment d'inertie que celui de l'étiquette latérale théorique, comme montré sur la figure 15.

Pour définir la plage de recherche pour l'angle de rotation, l'étiquette latérale théorique est utilisée selon trois aspects différents.

En effet, selon un premier aspect, montré sur la figure 16, le calcul de l'orientation de l'étiquette latérale théorique est basé sur la totalité de l'étiquette latérale théorique afin d'obtenir les angles correspondants 91tot et 92tot et l'écart angulaire A9tot. Selon un deuxième aspect, montré sur la figure 17, le calcul de l'orientation de l'étiquette latérale théorique est basé sur la moitié supérieure de l'étiquette latérale théorique afin d'obtenir les angles correspondants 91sup and 92sup et l'écart angulaire correspondant A9sup.

Selon un troisième aspect, montré sur la figure 18, le calcul de l'orientation de l'étiquette latérale théorique est basé sur la moitié inférieure de l'étiquette latérale théorique afin d'obtenir les angles correspondants 91inf and 92inf et l'écart angulaire correspondant Δθίηί.

Pour chacun des trois aspects, les écarts entre les angles d'orientation à Tl et à T2 sont calculés comme indiqués ci-dessous :

A0tot = |91tot - 92tot|,

A9sup = |91sup - 92sup|,

A9inf = |91inf - 92inf|,

Un intervalle des valeurs possibles de l'angle d'orientation optimal 9opt est défini ci- dessous :

Intervalle optimal : [min (A9tot, A9sup, Δθίηί) - 0,2°; max (Δθίοί, A9sup, Δθίηί) + 0,2°].

Le procédé de détermination des paramètres de recalage primaire 1002 comprend une deuxième étape 150 de recalage optimal concernant la rotation et la translation de la boîte de Pétri entre les temps Tl et T2, le centre de la rotation étant assimilé au centre de la boîte de Pétri.

L'étape de recalage optimal permet de déterminer, à partir de l'image acquise au temps T2, des paramètres optimaux de recalage primaire comprenant un angle de rotation et un vecteur de translation afin d'appliquer ces paramètres à l'image acquise au temps Tl pour que les deux images acquises soient superposables. La détermination des paramètres optimaux est obtenue en utilisant une fonction de corrélation croisée 2D, telle que disponible dans le logiciel Matlab™.

Pour estimer les paramètres optimaux, on test toutes les valeurs d'angle de rotation dans l'intervalle optimal avec un pas de 0,1° de l'image à T2 et on calcul la corrélation croisée entre l'image à T2 tournée et l'image à Tl. Chaque valeur résultante représente la translation optimale pour chaque valeur d'angle de rotation dans l'intervalle optimal. Le maximum de ces valeurs résultantes permet de déterminer la rotation optimale et sa translation optimale associée. Le deuxième procédé du procédé global de l'invention concerne le recalage secondaire ou recalage fin 28. Le deuxième procédé peut éventuellement être réalisé en premier lieu, sans que le premier procédé relatif au recalage primaire n'ait été réalisé au préalable. La figure 19 montre une première image ou image initiale d'une boîte de Pétri, acquise au moyen du système d'analyse d'images 10 à un temps Tl, avant l'incubation de la boîte de Pétri par exemple.

La figure 20 montre une deuxième image ou image finale de la boîte de Pétri montrée sur la figure 19 et acquise à un temps T2 après l'incubation de la boîte de Pétri au sein du dispositif d'incubation 20. Comme montré sur la figure 20, des colonies de micro-organismes sont présentes au sein de la boîte de Pétri au temps T2. Comme montré sur la figure 21, le procédé de recalage secondaire comprend une étape 260 relative à la création d'un maillage de sous-images au sein de l'image initiale II acquise au temps Tl et de l'image finale 12 acquise au temps T2, une étape 270 relative à un procédé de recalage local et une étape 280 relative à un procédé de recalage global.

Pour l'étape 260 de création d'un maillage de sous-images, une image initiale II est quadrillée, comme montré sur la figure 22, selon un maillage initial qui comprend des cellules initiales. La figure 23 montre un exemple d'une cellule initiale. L'image finale 12 est quadrillée comme montré sur la figure 24 selon un maillage final qui comprend des cellules finales. La figure 25 montre un exemple d'une cellule finale. La dimension d'une cellule initiale est supérieure à la dimension d'une cellule finale.

Comme montré sur la figure 21, le procédé de recalage secondaire comprend une étape 270 de recalage local concernant chaque couple de sous-images obtenu à l'issue de l'étape 260 décrite ci-dessus. L'étape 270 de recalage local, montré sur la figure 26, comprend une sous-étape 272 pour calculer la valeur d'un déplacement optimal entre les deux sous-images considérées au sein du couple de sous-images. Le déplacement optimal est estimé en considérant le déplacement de la cellule finale sur la cellule initiale correspondante. Dans la mesure où la dimension de la cellule finale est inférieure à la dimension de la cellule initiale, le déplacement de la cellule finale est possible selon une marge de déplacement paramétrée par l'utilisateur. Pour chaque déplacement une valeur d'un paramètre de similarité est calculée dans une étape 273. Ce paramètre de similarité peut être, par exemple, la mesure de distance ou d'information mutuelle. Dans la présente invention, le paramètre de similarité correspond au coefficient de corrélation. Pour chaque déplacement de la sous-image finale sur la sous-image initiale, une valeur du coefficient de corrélation est obtenue. La valeur maximale de ce coefficient de corrélation définit le déplacement optimal pour un couple de sous- images. La sous-étape 274 permet de calculer la valeur d'un paramètre de contraste entre les deux sous-images considérées au sein du couple de sous-images recalées selon le déplacement optimal estimé dans la sous-étape 272. Le second paramètre de contraste permet d'évaluer la quantité d'information présente au sein de chaque sous-image du couple de sous-images. Dans la présente invention, la valeur du paramètre de contraste est calculée comme l'écart-type de la distance entre les couleurs calculée pixel à pixel selon la formule ci-dessous : J E \ (X ~~ E [X] iEixn E[Xf

Le paramètre de contraste permet de déterminer les sous-images qui sont pertinentes pour le calcul ultérieur dans le procédé de recalage global.

Ainsi l'étape 270 permet d'obtenir un ensemble de transformations locales en translation entre l'image initiale acquise à Tl montrée sur la figure 27 et l'image finale acquise à T2 montrée sur la figure 28 par des vecteurs pour chaque point du maillage. La figure 29 montre en détail le contenu du carré représenté sur la figure 38. Comme montré sur les figures 28 et 29, ces transformations locales sont représentées par des vecteurs.

Comme montré sur la figure 21, le procédé de recalage d'images selon l'invention comprend une étape 280 de recalage global. L'étape 280 permet d'utiliser les résultats de l'étape 270 afin d'en déduire, pour chaque couple d'images II et 12, des paramètres de recalage global.

La détermination d'une transformation globale est réalisée sur la base, d'une part, de l'ensemble des transformations locales c'est-à-dire des translations et, d'autre part, des paramètres de similarité et de contraste pour tous les couples de sous- images.

Les valeurs des paramètres de similarité et de contraste permettent de sélectionner les couples de sous-images les plus pertinents. Ainsi, les transformations locales les plus pertinentes correspondantes peuvent être déduites, par exemple, en appliquant une étape de comparaison avec une valeur seuil.

Ensuite, une transformation globale est déterminée à partir des transformations locales les plus pertinentes par exemple en considérant la moyenne, la médiane ou la transformation locale la plus fréquemment représentée parmi les transformations locales. Comme montré sur la figure 2, le procédé global selon l'invention comprend une étape 30 de détermination préalable de croissance de micro-organismes sur l'image 12 acquise au temps T=T2. L'étape 30 de la figure 2 comprend l'application d'un procédé de détermination préalable de croissance afin de détecter l'éventuelle croissance de microorganismes au sein de la boîte de Pétri après incubation de ladite boîte de Pétri. Le procédé de détermination préalable de croissance comprend un procédé d'identification de zones potentielles de croissance et un procédé d'analyse temporelle desdites zones de croissance. Comme montré sur la figure 30, le procédé de détermination préalable de croissance comprend une étape 300 relative à l'obtention d'une ou plusieurs images de la boîte de Pétri, à l'issu du procédé de recalage secondaire, selon différentes conditions d'acquisition relatives notamment au dispositif d'illumination 14.

Le procédé de détermination préalable de croissance comprend également une étape 302 pour créer une image numérisée à partir des images obtenues lors de l'étape 300. Le procédé de détermination préalable de croissance comprend une étape 304 pour détecter la présence d'une densité élevée de colonies de micro-organismes ainsi que les colonies de micro-organismes de forme non-circulaire. Ainsi, l'étape 304 permet de détecter une zone de croissance de micro-organismes, la forme de ladite zone de croissance étant variable. Afin de faciliter l'application ultérieure des calculs, ladite zone de croissance est définie au moyen d'un rectangle dit « rectangle englobant ». Le « rectangle englobant » est centré sur la zone de croissance et permet d'agrandir la surface de la zone de croissance selon une dimension déterminée à l'avance, par exemple 10 pixels autour de la zone de croissance. Le procédé de détermination préalable de croissance comprend ensuite une étape 306 pour associer un élément de localisation tel qu'un rectangle englobant, à chaque élément de croissance déterminé par le procédé de détermination préalable de croissance.

Ainsi, le procédé de détermination préalable de croissance de micro-organismes permet de déterminer si des colonies de micro-organismes sont présentes au sein de la boîte de Pétri.

Afin d'améliorer les résultats issus du procédé de détermination préalable, un procédé d'analyse temporelle 40 de la croissance de micro-organismes est ensuite appliqué afin de valider chaque rectangle englobant une zone de croissance obtenue lors de l'application du procédé de détermination préalable de croissance de microorganismes et obtenir ainsi l'indication de la croissance effective de micro - organismes au sein de la boîte de Pétri. Le procédé d'analyse temporelle considère donc les images obtenues après l'application du procédé de recalage primaire optionnel et du procédé de recalage secondaire et, pour lesquelles, l'application du procédé de détermination préalable de croissance a généré des rectangles englobant des zones potentielles de croissance au sein des images correspondantes.

Selon le mode de réalisation de la présente invention, les images initiales et les images finales considérées correspondent à des images acquises avec les trois conditions d'acquisition suivantes : « backlight », « bottom» et « médian ». En effet, chaque condition d'acquisition du dispositif d'illumination 14 permet une observation spécifique de la boîte de Pétri. Ainsi, l'illumination de type « backlight » permet d'observer par transparence les éléments contenus au sein de la boîte de Pétri. Les conditions d'acquisition de type «bottom» et «médian» permettent d'observer les éléments localisés à la surface du milieu de culture gélosé au sein de la boîte de Pétri. Comme montré sur la figure 31, le procédé d'analyse temporelle comprend une étape 420 pour appliquer un procédé de prétraitement d'images aux images initiales et aux images finales, une étape 430 pour évaluer des critères d'évolution basés sur les images initiales et les images finales et une étape 440 pour déterminer la croissance effective de micro-organismes.

L'étape 420, relative au procédé de prétraitement d'images et montré sur la figure 32, comprend une sous-étape 424 pour améliorer le contraste au sein de chaque zone relative à rectangle englobant de ladite zone au sein des sous-images. Ainsi, la sous-étape 424 concerne l'application d'un procédé d'ajustement de contraste d'une image tel que l'égalisation d'histogramme connue dans l'art antérieur. Ainsi, les images obtenues après l'étape 424 contiennent différents niveaux de gris comme montré sur la figure 33, en fonction des différentes conditions d'illumination, c'est- à-dire « backlight », « bottom » et « médian ».

L'étape 420 comprend ensuite une sous-étape 425 pour réduire, au sein de l'image, le bruit généré par l'augmentation de contraste. Ainsi, la sous-étape 425 concerne l'application d'un filtre médian tel que connu dans l'art antérieur, par exemple un filtre pour une fenêtre 3x3.

Ainsi, à l'issue de la sous-étape 425, les zones relatives aux rectangles englobant les zones potentielles de croissance contiennent des transitions montrées sur la figure 34, qui sont plus lisses et plus homogènes que les transitions desdites zones montrées sur la figure 33, à l'issue de l'étape 424.

Comme montré sur la figure 31, le procédé d'analyse temporelle comprend une étape 430 pour déterminer, au sein des sous-images, la valeur d'au moins un paramètre d'évolution.

Selon un mode de réalisation préféré de l'invention, la valeur de deux paramètres d'évolution doit être déterminée pour chaque ensemble de sous-images. L'étape 430 comprend une sous-étape 461 montrée sur la figure 35 et relative au premier paramètre d'évolution qui permet de caractériser l'évolution de la couleur au sein des sous-images entre les temps d'acquisition Tl et T2. Selon un mode de réalisation de l'invention, le premier paramètre d'évolution correspond à un paramètre de contraste qui permet d'évaluer la stabilité globale de la sous-image entre les temps d'acquisition Tl et T2. Selon un mode de réalisation de l'invention, le premier paramètre d'évolution correspond à un paramètre de corrélation. L'étape 430 comprend une sous-étape 462 montrée sur la figure 35 et relative au deuxième paramètre d'évolution

Selon un mode de réalisation de l'invention, la méthode de mesure de contraste est basée sur le calcul de la déviation standard, par exemple au moyen de la fonction « std » du logiciel Matlab™, concernant les distances mesurées, par exemple au moyen de la fonction « pdist » du logiciel Matlab™, entre deux pixels, et ce pour tous les pixels localisés à l'intérieur du rectangle englobant :

C = std(pdtst{array) }

Pour la mesure du contraste entre des voxels dans une sous-image en couleur, la variable "array" correspond à un vecteur colonne n x 3.

Pour la mesure du contraste entre des pixels dans une sous-image en niveaux de gris, la variable "array" correspond à un vecteur colonne n. Comme montré sur la figure 36, la sous-étape 461 relative au calcul du paramètre de contraste comprend une pluralité d'étapes basées sur la mesure du contraste définie avec l'équation C ci-dessus, selon un mode de réalisation de l'invention.

La première étape 463 correspond à la mesure du contraste pour chaque sous- image d'une image initiale et pour chaque condition d'illumination « backlight », « bottom» et « médian ». La deuxième étape 464 correspond à la mesure du contraste pour chaque sous- image d'une image finale et pour chaque condition d'illumination « backlight », « bottom» et « médian ».

La troisième étape 465 correspond au calcul des ratios des mesures de contraste en considérant deux à deux les mesures de contraste d'une sous-image aux temps Tl et T2, pour chaque condition d'illumination.

La quatrième étape 466 correspond à la détermination de la valeur minimale des ratios issus de l'étape 465 pour déterminer la valeur du paramètre de contraste.

La sous-étape 461 correspond à l'application de l'équation ContC ci-dessous :

_^ . , c {Ti backlight} c {T, bottom} e {T, médian}

ContC= min(— ¹ , ¹ , ¹ )

Y ₂ backlight} T ₂ bottom} {T ₂ médian}

Selon un autre mode de réalisation de l'invention, une autre équation ContC_alt peut être appliquée en utilisant les mêmes mesures de contraste :

ContC ait = min(

L'équation ContC_alt présente des avantages par rapport à l'équation ContC. En effet, les valeurs obtenues avec l'équation ContC_alt sont toujours comprises dans l'intervalle [0 ; 1]. En effet, pour un élément stable, c'est-à-dire un élément qui n'a pas évolué au cours du temps, la valeur obtenue avec l'équation ContC_alt est proche de 1, selon une marge d'appréciation que l'utilisateur peut fixer. Ainsi, pour une marge d'appréciation de 0,2, la valeur obtenue pour un élément stable est comprise dans l'intervalle [0,8 ; 1]. Pour un élément de croissance, c'est-à-dire un élément ayant évolué au cours du temps, la valeur obtenue avec l'équation ContC_alt est proche de 0, selon une marge d'appréciation que l'utilisateur peut fixer. Ainsi, pour une marge d'appréciation de 0,2, la valeur obtenue pour un élément ayant évolué est comprise dans l'intervalle [0 ; 0,2].

Comme montré sur la figure 37, la sous-étape 462 comprend une pluralité d'étapes. Ainsi, la première étape 450 correspond à la sélection de la sous-image initiale en couleur la plus contrastée, parmi les sous-images identiques acquises selon les trois conditions d'acquisition « backlight », « bottom» et « médian ».

La deuxième étape 451 correspond à la sélection d'une sous-image finale prétraitée correspondante et telle qu'obtenue à l'issue des étapes 424 et 425.

La sélection de l'image la plus contrastée au sein des étapes 450 et 451 est basée sur la méthode de mesure de contraste décrite au sein de l'étape 461. La troisième étape 452 comprend la transformation de la sous-image colorée en sous-image en niveaux de gris.

La quatrième étape 453 comprend le calcul du coefficient de corrélation entre la sous-image acquise à Tl et la sous-image acquise à T2, toutes deux en niveaux de gris, en utilisant la fonction « corr2 » du logiciel Matlab™.

La cinquième étape 454 correspond au calcul de la valeur du paramètre de corrélation qui correspond à la valeur maximale entre les deux valeurs de coefficients de corrélation calculées à l'étape 453, au moyen de l'équation ci- dessous :

CorrP = max[corr2(Tlgris, T2gris), corr2 (Tlaméliorée, T2améliorée)]

Comme montré sur la figure 31, le procédé d'analyse temporelle comprend une étape 440 pour déterminer la croissance effective de micro-organismes au sein de la boîte de Pétri, à l'issue d'une période déterminée d'incubation. L'objectif de l'étape 440 consiste à vérifier si l'élément repéré par une zone de croissance déterminée lors de l'application du procédé de détermination préalable de croissance au sein d'une sous-image, lors de l'étape 30, correspond effectivement à un élément de croissance tel qu'une colonie ou un ensemble de colonies de micro- organismes.

Comme montré sur la figure 38, l'étape 440 comprend une étape 470 afin de lister les différentes valeurs obtenues pour les deux paramètres d'évolution lors de l'étape 430. Ces valeurs sont associées à chaque rectangle englobant des zones potentielles de croissance de la boîte de Pétri pour le temps d'acquisition T2

L'étape 440 comprend ensuite une étape 472 relative à l'application d'une méthode de régression logistique. Cette méthode permet d'obtenir une séparation entre deux groupes de données à discriminer sur la base de plusieurs paramètres au moyen d'une droite de régression définie. Au sein de la présente invention, la méthode de régression logistique utilisée permet de distinguer deux groupes de données. Le premier groupe de données concerne les zones potentielles de croissance qui sont des zones de croissance avérées, c'est-à-dire des zones positives. Le deuxième groupe de données concerne les zones potentielles de croissance qui ne correspondent pas à des zones de croissance. Il s'agit de zones négatives. La méthode de régression logistique au sein de la présente invention est basée sur deux paramètres, CorrP et ContC, selon l'équation ci-dessous :

CL = -Î0..3 · CorrP - 9 £5 · ContC_v2 + 12.

De plus, pour chaque valeur des paramètres de corrélation CorrP et de contraste ContC, l'équation P ci-dessous permet de déterminer une probabilité de présence d'un élément de croissance tel que des colonies de micro-organismes à partir de la valeur de CL définie ci-dessus.

î

P = : r

î + sxîî ( -CL } Afin de classifier les zones potentielles de croissance l'utilisateur peut fixer un critère de décision sur les valeurs de P obtenues.

Selon un mode de réalisation de la présente invention montré sur la figure 39, le critère de décision P est fixé avec les conditions ci-dessous :

- si P<0,25 alors la zone potentielle de croissance correspond à un élément stable dans le temps donc la zone potentielle de croissance n'est pas une zone de croissance, il s'agit d'une zone dite négative ;

- si P>0,6 alors la zone potentielle de croissance correspond à un élément qui évolue au cours du temps et donc la zone potentielle de croissance est une zone de croissance avérée, il s'agit d'une zone dite positive ;

- si 0,25<P<0.6 alors il n'est pas possible de déterminer si la zone potentielle de croissance correspond à un élément stable ou à un élément qui évolue au cours du temps, c'est-à-dire que la distinction entre une zone dite positive et une zone dite négative n'est pas possible.

Le procédé d'analyse temporelle permet ensuite de comparer dans une étape 474 la valeur de l'annotation locale avec le critère de décision P, donc de rechercher, pour chaque zone potentielle de croissance d'une image finale, la correspondance avec l'image initiale correspondante. Ainsi dans une étape 476, l'utilisateur peut déterminer l'absence ou la présence de colonies de micro-organismes au sein de la boîte de Pétri analysée. Après le calcul des paramètres de corrélation et de contraste, le procédé d'analyse temporelle permet de déterminer si la zone de croissance correspond à un élément ayant évolué ou n'ayant pas évolué au cours du temps.

Lorsque la zone de croissance correspond à un élément ayant évolué au cours du temps, le procédé d'analyse temporelle permet de confirmer la validité de la zone de croissance telle que déterminée par le procédé de détermination préalable de croissance. Ainsi, la zone de croissance correspond effectivement à une croissance de micro-organismes.

Lorsque la zone de croissance correspond à un élément n'ayant pas évolué au cours du temps, tel qu'un défaut lié au matériau de la boîte de Pétri, le procédé d'analyse temporelle permet d'obtenir une correction de la zone de croissance. Ainsi, la zone de croissance corrigée indique que l'élément observé au sein de ladite zone de croissance ne correspond pas à une croissance de micro-organismes.

Après la détermination de la croissance, l'utilisateur a la possibilité d'identifier la nature des colonies de micro -organismes au moyen d'un système d'identification (non montré).