FR 2909444 A1 20080606 FR 0610576 A 20061204 1 PROCEDE ET INSTALLATION D'IMAGERIE La présente invention est relative aux procédés et installations d'imagerie. Le document GB 2 350 187 décrit un exemple d'un procédé dans lequel, lors la détection d'un signal, on vérifie si la détection correspond plutôt à un mode de comptage, ou à un mode d'intégration, et dans lequel, en fonction de cette détermination, on obtient: une valeur estimée du signal correspondant soit au signal de comptage, soit au signal d'intégration, voire à une combinaison de ces deux signaux. L'intérêt d'un tel procédé est en particulier de pouvoir améliorer l'échelle dynamique de détection. Toutefois, la mesure décrite dans ce document est obtenue à partir d'un photomultiplicateur, de sorte qu'il s'agit d'une approche mono-pixel, pour laquelle on cherche à connaître le niveau de détection, mais pas l'emplacement d'origine du signal. Il existe donc un besoin d'un procédé d'imagerie à bonne échelle dynamique de détection, permettant de fournir en plus une image de la détection. A cet effet, l'invention se rapporte principalement à un procédé d'imagerie comprenant les étapes suivantes : (a)on fournit pour chacune d'un ensemble de zones (x ; y) à observer d'une région d'observation, et pour un instant 5t d'une période d'observation T, une mesure Fo(x, y, 5t) des signaux détectés correspondant à une émission de rayonnement en provenance de ladite zone pendant l'instant 5t, (b)on fournit d'une part une première valeur F1(x, y, ôt),dite d' intégration , estimative de l'intensité des signaux détectés en provenance de ladite zone pendant l'instant, et d'autre part une deuxième valeur F2(x, y, 5t) dite de comptage , estimative du nombre de signaux détectés en provenance de ladite zone pendant l'instant, 2909444 2 F2(x, y, ôt) étant obtenue en affectant, pour un ensemble de zones contiguës où la mesure Fo(x, y, 5t) est non nulle formant une tache de détection : une unité de signal au barycentre de ladite 5 tache, et zéro unité de signal aux autres zones de la tache, (c) en au moins une zone considérée (xo ; yo) dudit ensemble de zones, on estime une valeur Fe(xo, yo, 5t) du 10 nombre de signaux en la zone considérée (xo ; yo) pendant ledit instant à partir d'au moins une combinaison des première F1 (xo, yo, 5t) et deuxième F2(xo, yo, 5t) valeurs, dans lequel, au cours de l'étape (c), on met en oeuvre les étapes suivantes : 15 (cl) on définit un voisinage R(xo ; yo) formant une partie dudit ensemble de zones, ledit voisinage comprenant au moins la zone considérée (xo ; yo) et des zones voisines (xi ; yi) distinctes de la zone considérée, (c2) on estime un critère de détection dans le 20 voisinage à partir des mesures Fo(x, y, 6t) pour des zones voisines (x ; y) comprises dans ledit voisinage R(xo ; yo) au cours de la période d'observation T, et (c3) on estime la valeur Fe(xo, yo, 6t) du nombre de signaux en fonction en outre au moins dudit critère de 25 détection dans le voisinage. Grâce à ces dispositions, on obtient un procédé qui est bien adapté à la détection multi-pixels. En particulier, l'estimation de la valeur Fe du nombre de signaux est obtenue en prenant en compte la détection dans 30 le voisinage spatial et/ou temporel du pixel considéré. Dans des modes de réalisation préférés de l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes : - lorsque le critère de détection dans le 35 voisinage indique un niveau de détection bas, on estime la 2909444 3 valeur Fe(xo, yo, bt) du nombre de signaux comme la deuxième valeur F2(xo, yo, bt) ; - lorsque le critère de détection dans le voisinage indique un niveau de détection élevé, on estime 5 la valeur Fe(xo, yo, bt) du nombre de signaux comme la première valeur F1 (xo, yo, 5t) ; lorsque le critère de détection dans le voisinage indique un niveau de détection intermédiaire entre les niveaux bas et élevé, on estime la valeur Fe(xo, 10 yo, 5t) du nombre de signaux comme une combinaison de la première valeur F1 (xo, yo, bt) et la deuxième valeur F2(xo, yo, bt) chacune étant pondérée par un coefficient non nul; - lorsque le critère de détection dans le voisinage indique ledit niveau de détection intermédiaire, 15 on estime la valeur Fe(xo, yo, 5t) du nombre de signaux comme une combinaison linéaire de la première valeur F1(xo, yo, 5t) et la deuxième valeur F2(xo, yo, bt) selon Fe (xo, yo, 5t) =a F1. (xo, yo, 5t)+(1ù a) F2 (xo, yo, 5t), a étant un coefficient compris dans l'intervalle 20 ]0;1[ dépendant du critère de détection dans le voisinage ; - pour la zone (xo, yo) et l'instant bt considérés, a dépend de la mesure Fo(xo, yo, bt) en cette zone pour cet instant ; - au cours de l'étape (c2), on estime le critère 25 de détection dans le voisinage à partir au moins des mesures Fo(x, y, bt) pour des zones voisines (xi ;Yi) comprises dans ledit voisinage R(xo ; yo) au cours de l'instant bt ; - le critère de détection est relatif à un nombre 30 de zones (x1 ;yl) dudit voisinage dans lesquelles la mesure Fo(x, y, bt) est inférieure à un seuil prédéterminé ; - le critère est relatif à une moyenne de la mesure Fo(x, y, bt) dans ledit voisinage ; -les étapes (a), (b) et (c) sont mises en oeuvre 35 au moins pour un premier instant 5t1 et un deuxième instant 2909444 4 5t2 distinct du premier instant, et on estime pour ledit premier instant ledit critère de détection dans le voisinage au moins à partir des mesures Fo(x, y, 5t2) pour des zones voisines (x1 ; y1) comprises dans ledit voisinage 5 R(xo ; yo) au cours du deuxième instant ôte ; - l'étape (c) est mise en oeuvre en une pluralité de zones considérées ; - au cours de l'étape (b), on obtient ladite première valeur F1(x, y, 5t) en fonction de ladite mesure 10 Fo(x, y, 6t) ; - avant l'étape (a), on met en oeuvre une étape (z) au cours de laquelle on détecte pour chacune d'un ensemble de zones (x ; y) à observer d'une région d'observation, et pour un instant ôt d'une période d'observation, un signal 15 correspondant à une émission de rayonnement en provenance de ladite zone pendant l'instant 5t, ladite détection fournissant ladite mesure Fo(x, y, 6t) ; - au cours de l'étape (z), on détecte pour chacune dudit ensemble de zones (x ; y) à observer, et pour ledit 20 instant 5t, un signal correspondant à une émission de rayonnement lumineux en provenance de ladite zone pendant l'instant 5t ; - au cours de l'étape (c2), on estime ledit critère de détection dans le voisinage à partir en outre de 25 la mesure Fo (xo, yo, 5t) pour la zone considérée (xo, yo). Selon un autre aspect, l'invention se rapporte à un produit programme d'ordinateur comprenant des portions de code de programmes adaptées pour mettre en oeuvre un tel procédé lorsque ledit programme est exécuté sur une machine 30 programmable. Selon un autre aspect, l'invention se rapporte à une installation d'imagerie comprenant : (A)une mémoire stockant, pour chacune d'un ensemble de zones (x ; y) à observer d'une région d'observation, et 35 pour un instant 5t d'une période d'observation, une mesure 2909444 5 Fo(x, y, 5t) de signaux détectés correspondant à une émission de rayonnement en provenance de ladite zone pendant l'instant 6t, (B)un calculateur adapté pour fournir d'une part 5 une première valeur F1(x, y, ôt) dite d' intégration estimative de l'intensité des signaux détectés en provenance de ladite zone pendant l'instant, et d'autre part une deuxième valeur F2(x, y, ôt) dite de comptage estimative du nombre de signaux détectés en provenance de 10 ladite zone pendant l'instant, F2(x, y, ôt) étant obtenue en affectant, pour un ensemble de zones contiguës où la mesure Fo(x, y, 6t) est non nulle formant une tache de détection : une unité de signal au barycentre de ladite 15 tache, et zéro unité de signal aux autres zones de la tache, ledit calculateur étant adapté pour estimer, en au moins une zone considérée (xo ; yo) choisie dans ledit 20 ensemble de zones, une valeur Fe(xo, yo, 6t) du nombre de signaux en la zone considérée (xo; yo) pendant ledit instant à partir d'au moins une combinaison des première F1 (xo, yo, 6t) et deuxième F2 (xo, yo, ôt) valeurs, le calculateur étant adapté pour définir un 25 voisinage R(xo ; yo) formant une partie dudit ensemble de zones, ledit voisinage comprenant au moins la zone considérée (xo ; yo) et des zones voisines (xi ; yl) distinctes de la zone considérée, le calculateur étant adapté pour estimer un critère 30 de détection dans le voisinage, au moins à. partir des mesures Fo(x, y, 5t) pour des zones voisines (xi ; y1) comprises dans ledit voisinage R(xo ; yo) au cours de la période d'observation T, et le calculateur étant adapté pour estimer la valeur 35 Fe(xo, yo, 5t) du nombre de signaux en fonction en outre au 2909444 6 moins dudit critère de détection dans le voisinage. Dans certains modes de réalisation, on peut en outre prévoir : - une unité de détection adaptée pour détecter des 5 signaux correspondants à une émission de rayonnement en provenance de ladite zone pendant l'instant 5t, et pour stocker la mesure correspondante Fo(x, y, 6t) dans ladite mémoire ; - l'unité de détection est adaptée pour détecter 10 des signaux optiques. D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante d'une de ses formes de réalisation, donnée à titre d'exemple non limitatif, en regard des dessins joints. 15 Sur les dessins : - la figure 1 représente schématiquement une installation d'imagerie, - la figure 2 représente schématiquement un exemple d'obtention des valeurs F1(x, y, ôt) et F2(x, y, 5t), 20 - la figure 3 représente schématiquement une étape d'obtention de Fe(xo, yo, ôt), et - les figures 4a et 4b sont des graphes présentant deux variantes possibles pour a(C). Sur les différentes figures, les mêmes références 25 désignent des éléments identiques ou similaires. La figure 1 :représente schématiquement une installation 1 d'imagerie comprenant une enceinte dans laquelle est disposé un échantillon 3 à imager. L'échantillon 3 émet un rayonnement 4 détecté et 30 éventuellement amplifié par un détecteur o-i capteur 5 spécialement adapté pour détecter le rayonnement 4. Le détecteur 5 comprend une pluralité de cellules 6 arrangées classiquement en lignes et en colonnes (une seule ligne visible sur la figure 1) dans le plan perpendiculaire au 35 plan de la figure 1. Le détecteur 5 est relié en sortie à 2909444 7 un système informatique 7 comprenant une unité centrale 8 comprenant une mémoire 9 dans laquelle sont stockées les mesures obtenues par le détecteur 5, et une unité de traitement 10 programmée pour traiter les données stockées 5 dans la mémoire 9 de la manière explicitée par la suite. Le système informatique 7 peut également comporter un écran 11 permettant d'afficher à un utilisateur le signal détecté et/ou le résultat du traitement par l'unité de traitement 10, ou autre. 10 A titre d'exemple purement illustratif, on a représenté sur la figure 1 une expérience de détection en luminescence au cours de laquelle on détecte un rayonnement lumineux 4 émis depuis l'intérieur d'un animal 3, par exemple une souris modifiée génétiquement, disposée dans 15 une boite noire parfaitement opaque à la lumière extérieure. Ce signal lumineux est par exemple représentatif d'une réaction chimique se produisant dans l'animal et peut servir de traceur pour une telle réaction chimique. Le rayonnement 4 peut également servir à 20 quantifier l'expression d'un gène donné porté par l'animal 3 et qui, lorsqu'il s'exprime, occasionne la réaction chimique mentionnée précédemment. Pour mesurer un tel rayonnement, on utilise à titre de détecteur 5 une caméra CCD refroidie, un ICCD, un EMCCD, 25 (CCD à multiplication interne) ou autre présentant par exemple environ un (ou plusieurs) million de cellules de détection 6 disposées en mille colonnes de mille rangées. L'observation est effectuée pendant une période d'observation T, par exemple de quelques minutes, ou de 30 plusieurs heures, subdivisée en instants d'observation 5tl, 6t2,..., 6ti, , 6tn , pendant chacun desquels une mesure Fo(x, y, 6ti) est obtenue pour chacune des cellules 6, repérée par ses coordonnées (x, y) dans le plan de détection du détecteur 5 par rapport à une origine 35 prédéterminée, et stockée dans la mémoire 9 du système 2909444 8 informatique 7. La figure 2 représente à gauche, le plan du détecteur 5, et à droite en haut une portion agrandie du signal mesuré dans une région 12 de 20 cellules de 5 détection. En bas de la figure 2, sont représentés schématiquement des première FI et deuxième F2 valeurs obtenues à partir de la mesure Fo tel qu'explicité par la suite. Comme représenté sur la figure 2, la détection d'un 10 évènement, correspondant à l'émission d'un rayonnement 4 unique de l'échantillon 3 peut être détectée simultanément par plusieurs cellules de détection 6 au cours d'un instant ôt donné. Par exemple, une cellule donnée représentée au centre de la région agrandie 12 sur la figure 2 étant la 15 plus proche de l'emplacement effectif de l'émission du rayonnement 4 correspondra à une mesure maximale, schématisée par une assez forte densité de hachures sur la représentation de la mesure Fo de la figure 2. Les cellules avoisinantes pourront toutefois avoir détecté une 20 partie de cet unique rayonnement, de sorte qu'une mesure Fo non nulle est également mémorisée pour ces cellules, comme représenté par les hachures de faible densité sur la représentation de la mesure Fo sur la figure 2. Le détecteur 5 présentant par ailleurs un certain seuil de 25 détection, la mesure Fo pour les cellules (x, y) dans lesquelles ce seuil n'a pas été atteint est stockée égale à zéro dans la mémoire 9. Ces cellules sont représentées en blanc sur la représentation de la mesure Fo sur la figure 2. 30 Pour chaque pixel repéré par ses coordonnées (x, y), la mesure Fo stockée dans la mémoire 9 correspond donc à zéro si un signal trop faible a été détecté, pendant l'instant 5t, ou à une mesure non nulle pouvant correspondre à une partie d'un ou de plusieurs rayonnements 35 ayant eu lieu à proximité du pixel considéré. 2909444 9 Une fois les mesures Fo(x, y, ôt) stockées dans la mémoire 9 pour l'ensemble des instants 5t de la période d'observation T, on obtient pour chaque instant 5t une image dite de comptage F2 estimative du nombre 5 d'évènements de rayonnement émis en la zone considérée, pendant l'instant 5t. Par exemple, pour chaque tache de détection, correspondant à un ensemble de pixels contigus de mesures Fo non nulles, on identifie le barycentre de cette tache, et on lui assigne une valeur dite de 10 comptage , correspondant sensiblement à une unité de signal, par exemple normalisée. On assigne la valeur F2 = 0 aux autres pixels de la tache de détection. considérée. La valeur de comptage F2 peut également être corrigée si l'on connaît expérimentalement ou statistiquement la 15 relation pour le détecteur entre le comptage calculé comme explicité ci-dessus et le nombre de rayonnement effectivement émis dans la zone considérée donnant ce comptage calculé. Statistiquement, on peut en effet savoir à l'avance que pour un nombre No de rayonnement émis 20 pendant un instant ôt en une zone donnée, seul un nombre inférieur N1, dépendant de No, sera mesuré. Une image, dite image de comptage obtenue à partir des valeurs F2, fournit une très bonne résolution. De plus, il est connu que le comptage fournit des données 25 statistiquement optimum, indépendantes de toute source de fluctuations dues au capteur. Toutefois, une telle image peut s'avérer trop fausse pour la qualité des résultats recherchés lorsque de trop nombreux signaux proches sont émis pendant un même instant 5t, rendant difficile de 30 discerner les évènements individuels. En effet, dans les cas où de nombreux rayonnements sont détecr.és, pendant l'instant 6t, dans une même zone du détecteur 5, les taches associées à chaque rayonnement individuel se superposent de sorte qu'il n'est plus possible de distinguer, pendant 35 l'instant ôt, les évènements individuels. 2909444 10 Dans ce dernier cas, on peut préférer une image de valeur F1 dite d'intégration correspondant sensiblement, pour chaque pixel de coordonnée (x, y) à la quantité (intensité) de signal détecté pendant l'instant 6t 5 en ce pixel. Dans l'exemple présenté, F1 est. par exemple directement égal à Fo pour chaque pixel considéré. Une telle image présente une plus mauvaise résolution que l'image de comptage , mais présente l'avantage de fournir une valeur plus juste de la quantité incidente sur 10 le capteur, en particulier lors de fortes détections. Une fois le calcul des valeurs F1 et F2 effectué par l'unité de traitement 10 pour chaque pixel de coordonnées (x, y) et pour chaque instant 5t, et le résu=ltat mémorisé dans la mémoire 9, l'unité de traitement 10 applique le 15 traitement représenté schématiquement sur la figure 3. On considère un pixel considéré de coordonnée (xo ; yo). Par exemple, on considère tour à tour l'ensemble des pixels correspondant chacun à une cellule 6 du détecteur 5. En alternative, on considère uniquement les pixels pour 20 lesquels la valeur F2(xo, yo, bt) est non nulle. Pour le pixel considéré, de coordonnées (xo ; yo), on définit un voisinage R(xo, yo), constitué des quelques pixels, par exemple 25, ou 100 pixels, entourant le pixel considéré de coordonnées (xo ; yo). 25 On détermine un critère C(xo ; yo), relatif au niveau détection au voisinage R(xo ; yo) du pixel considéré de coordonnées (xo ; yo). Le critère C(xo ; yo) indique si, dans le voisinage R(xo ; yo) du pixel considéré, la détection est plutôt de haut niveau ou de bas niveau. 30 C(xo ; yo) varie évidemment selon le pixel considéré, et l'instant considéré. Pour déterminer le critère C(xo, yo), on utilise à titre d'exemple, la mesure Fo obtenue pour les pixels de coordonnées (x1 ; yl) autres que celui considéré, compris 35 dans le voisinage R(xo, yo). On utilise soit: directement 2909444 11 la mesure Fo, ou toute autre valeur obtenue d_Lrectement ou indirectement à partir de Fo, tel qu'en particulier F1 ou F2 en ces pixels. Par exemple, on utilise à titre de critère C(xo, 5 yo) la proportion de pixels du voisinage R(xo, yo) pour lesquels la mesure Fo (x1, y1) est nulle, ou inférieure à un seuil prédéterminé, pour l'instant ôt considéré. Si un grand nombre de pixels du voisinage présente une mesure Fo nulle, en d'autres termes si on détecte qu'on 10 est dans un cas de faible détection autour du pixel considéré, on décidera de donner à la valeur Fe(xo, yo, 5t) du nombre de signaux une valeur dans laquelle la valeur de comptage F2(xo, yo, 5t) du nombre de signaux sera prépondérante. En effet, comme expliqué précédemment, 15 F2 est préférable dans les régions de faible détection. Si on détecte que le nombre de pixeLs voisins à mesure nulle est inférieur à une valeur prédéterminée, à savoir qu'on est dans un cas de forte détection autour du pixel considéré (xo, yo), on choisira, pour la valeur Fe(xo, 20 yo, 6t) du nombre de signaux détectés une valeur dans laquelle la valeur d'intégration Fe(xo, yo, 5t) sera prépondérante. En effet, comme expliqué précédemment, F1 est préférable dans les régions de forte détection. En tout état de cause, on pourrait écrire Fe comme 25 une combinaison de F1 et de F2, par exemple une combinaison linéaire telle que : Fe(xo, yo, 5t)=a F1 (xo, yo, 5t)+(1- a) F2(xo, yo, 5t), où la fonction a du critère C(xo, yo) peut être comprise entre 0 et 1 en étant par exemple constante par morceau telle que 30 représentée sur la figure 4a, auquel cas la valeur Fe du nombre de signaux sera égale soit à F1 soit à. F2 selon la valeur du critère C par rapport à une valeur prédéterminée Co du critère C. En variante, comme représenté sur la figure 4b, la fonction a (C) peut comporter une portion 35 croissante, par exemple continue (en pointillés), par 2909444 12 exemple linéaire (en traits pleins), entre deux valeurs prédéterminées C1 et C2 du critère. Dans la portion intermédiaire comprise entre les valeurs C1 et C2 du critère, a est compris dans l'intervalle ]0 ;1[. On notera 5 que le paramètre a dépend du critère C en (xo, yo ) et dépend par conséquent directement ou indirectement de la mesure Fo en(xo, yo) ou/et au voisinage de (xo, yo) D'autres critères que le nombre de pixels présentant une mesure inférieure à un seuil prédéterminé 10 dans le voisinage R(xo, yo) peuvent être utilisés. A titre d'exemple, on peut par exemple prévoir que le critère correspondra à la moyenne de la mesure Fo(x, y, bt) sur l'ensemble des pixels de coordonnées (x ; y) du voisinage R(xo, yo) pendant l'instant bt. 15 Selon une autre variante, le critère C(xo, yo) peut également dépendre du temps. Dans ce cas, le critère, qui peut par exemple être l'un des deux critères précédemment décrit, est estimé non seulement pour l'instant bti pour lequel le calcul 20 Fe(xo, yo, bti) est en cours, mais également pour les instants bti_j précédents et bti+k suivants dans un voisinage temporel de l'instant considéré. L'utilisation d'un tel critère temporel peut permettre de détecter que, pour l'instant bti considéré, un rayonnement parasite de type 25 cosmique ou autre a été également capté par le détecteur 5. La comparaison de la mesure dans le voisinage du pixel considéré, de coordonnées (xo ; yo) pour l'instant bti avec les mesures dans ce même voisinage en des instants proches va permettre de déterminer si la mesure à :_'instant bti 30 inclut la mesure de ce rayonnement cosmique, et le calcul d'une valeur estimée Fe(xo, yo, bti) affranchie de ce rayonnement. Le procédé de traitement ici décrit en relation avec la figure 3 pourrait également s'appliquer à un signal 35 obtenu par détection de tout autre type de rayonnement que 2909444 13 des rayonnements lumineux, tel que présenté ic:_. REVENDICATIONS 1. Procédé d'imagerie comprenant les étapes suivantes : (a)on fournit pour chacune d'un ensemble de zones (x ; y) à observer d'une région d'observation, et pour un instant 5t d'une période d'observation T, une mesure Fo(x, y, 5t) des signaux (4) détectés correspondant à une émission de rayonnement en provenance de ladite zone pendant l'instant 5t, (b)on fournit d'une part une première valeur F1(x, y, 5t), dite d' intégration , estimative de l'intensité des signaux détectés en provenance de ladite zone pendant l'instant, et d'autre part une deuxième valeur F2(x, y, 5t) dite de comptage , estimative du nombre de signaux détectés en provenance de ladite zone pendant l'instant, F2(x, y, 5t) étant obtenue en affectant, pour un ensemble de zones contiguës où la mesure Fo(x, y, 5t) est non nulle formant une tache de détection : une unité de signal au barycentre de ladite tache, et zéro unité de signal aux autres zones de la tache, (c)en au moins une zone considérée (xo ; yo) dudit ensemble de zones, on estime une valeur Fe(xo, yo, 5t) du nombre de signaux en la zone considérée (xo; yo) pendant ledit instant à partir d'au moins une combinaison des première F1 (xo, yo, 6t) et deuxième F2 (xo, yo, ôt;) valeurs, dans lequel, au cours de l'étape (c), on met en oeuvre les étapes suivantes : (cl) on définit un voisinage R(xo ; yo) formant une partie dudit ensemble de zones, ledit voisinage comprenant au moins la zone considérée (xo ; yo) et des zones voisines (x1; y1) distinctes de la zone considérée, 2909444 15 (c2) on estime un critère de détection dans le voisinage à partir des mesures Fo(x, y, 5t) pour des zones voisines (xi ; y') comprises dans ledit voisinage R(xo ; yo) au cours de la période d'observation T, et 5 (c3) on estime la valeur Fe(xo, yo, ôt) du nombre de signaux en fonction en outre au moins dudit critère de détection dans le voisinage. 2. Procédé d'imagerie selon la revendication 1, dans lequel: 10 - lorsque le critère de détection dans le voisinage indique un niveau de détection bas, on estime la valeur Fe (xo, yo, 6t) du nombre de signaux comme la deuxième valeur F2 (xo, yo, 6t) , - lorsque le critère de détection dans le voisinage 15 indique un niveau de détection élevé, on estime la valeur Fe(xo, yo, ôt) du nombre de signaux comme la première valeur F1 (xo, yo, 5t). 3. Procédé d'imagerie selon la revendication 2, dans lequel, lorsque le critère de détection dans le 20 voisinage indique un niveau de détection intermédiaire entre les niveaux bas et élevé, on estime la valeur Fe (xo, yo, ôt) du nombre de signaux comme une combinaison de la première valeur F1 (xo, yo, 5t) et la deuxième valeur F2(xo, yo, ôt), chacune étant pondérée par un coefficient non nul. 25 4. Procédé d'imagerie selon la revendication 3 dans lequel, lorsque le critère de détection dans le voisinage indique ledit niveau de détection intermédiaire, on estime la valeur Fe(xo, yo, 6t) du nombre de signaux comme une combinaison linéaire de la première valeur F1(xD, yo, 6t) et 30 la deuxième valeur F2(xo, yo, 6t) selon : Fe(xo, yo, 5t) =a F1 (xo, yo, 5t) + (1- a) F2(xo, yo, 5t), a étant un coefficient compris dans l'intervalle ]0 ;1[ dépendant du critère de détection dans le voisinage. 5. Procédé d'imagerie selon la revendication 4, 35 dans lequel, pour la zone (xo, yo) et l'instant ôt 2909444 16 considéré, a dépend de la mesure Fo(xo, yo, 3t) en cette zone pour cet instant. 6. Procédé selon l'une des revendications précédentes dans lequel au cours de l'étape (c2), on estime 5 le critère de détection dans le voisinage à partir au moins de la mesure Fo(x, y, 5t) pour des zones voisines (xi ;yl) comprises dans ledit voisinage R(xo ; yo) au cours de l'instant 5t. 7. Procédé d'imagerie selon la revendication 6 dans 10 lequel le critère de détection est relatif à un nombre de zones (x1 ;y1) dudit voisinage dans lesquelles la mesure Fo(x, y, 5t) est inférieure à un seuil prédéterminé. 8. Procédé d'imagerie selon la revendication 6 dans lequel le critère de détection est relatif à une moyenne de 15 la mesure Fo(x, y, ôt) dans ledit voisinage. 9. Procédé d'imagerie selon l'une des revendications précédentes dans lequel les étapes (a), (b) et (c) sont mises en oeuvre au moins pour un premier instant (5t1) et un deuxième instant (5t2) distinct du premier 20 instant, et dans lequel on estime pour ledit premier instant ledit critère de détection dans le voisinage au moins à partir des mesures Fo(x, y, 5t2) pour des zones voisines (x1 ; y1) comprises dans ledit voisinage R(xo ; yo) au cours du deuxième instant (5t2). 25 10.Procédé d'imagerie selon l'une des revendications précédentes dans lequel l'étape (c) est mise en ouvre en une pluralité de zones considérées. 11.Procédé d'imagerie selon l'une des revendications précédentes dans lequel, au cours de l'étape 30 (b), on obtient ladite première valeur F1(x, y, ôt) en fonction de ladite mesure Fo(x, y, 5t). 12.Procédé d'imagerie selon l'une des revendications précédentes dans lequel, avant l'étape (a), on met en oeuvre une étape (z) au cours de laquelle on 35 détecte pour chacune d'un ensemble de zones (x ; y) à 2909444 17 observer d'une région d'observation, et pour un instant 5t d'une période d'observation, un signal (4) correspondant à une émission de rayonnement en provenance de ladite zone pendant l'instant 5t, ladite détection fournissant ladite 5 mesure Fo(x, y, 5t). 13.Procédé d'imagerie selon la revendication 12 dans lequel, au cours de l'étape (z), on détecte pour chacune dudit ensemble de zones (x ; y) à observer, et pour ledit instant 5t, un signal (4) correspondant à une 10 émission de rayonnement lumineux en provenance de ladite zone pendant l'instant 5t. 14.Procédé selon l'une quelconques des revendications précédentes, dans lequel, au cours de l'étape (c2), on estime ledit critère de détection dans le 15 voisinage à partir en outre de la mesure Fo (xo ;yo, 5t) pour la zone considérée (x:o ;yo) 15.Produit programme d'ordinateur comprenant des portions de code de programmes adaptées pour mettre en oeuvre un procédé selon l'une quelconque des revendications 20 précédentes lorsque ledit programme est exécuté sur une machine programmable. 16.Installation d'imagerie comprenant : (A)une mémoire (9) stockant, pour chacune d'un ensemble de zones (x ; y) à observer d'une région 25 d'observation, et pour un instant 5t d'une période d'observation, une mesure Fo(x, y, 6t) de signaux détectés correspondant à une émission de rayonnement en provenance de ladite zone pendant l'instant 5t, (B)un calculateur (10) adapté pour fournir d'une 30 part une première valeur F1(x, y, ôt) dite d' intégration estimative de l'intensité des signaux détectés en provenance de ladite zone pendant l'instant, et d'autre part une deuxième valeur F2(x, y, 5t) dite de comptage , estimative du nombre de signaux détectés en 35 provenance de ladite zone pendant l'instant, F2(x, y, 5t) 2909444 18 étant obtenue en affectant, pour un ensemble de zones contiguës où la mesure Fo(x, y, 6t) est non nulle formant une tache de détection : - une unité de signal au barycentre de ladite 5 tache, et zéro unité de signal aux autres zones de la tache, ledit calculateur (10) étant adapté pour estimer, en au moins une zone considérée (xo ; yo) choisie dans 10 ledit ensemble de zones, une valeur Fe(xo, yo, 5t) du nombre de signaux en la zone considérée (xo ; yo) pendant ledit instant à partir d'au moins une combinaison des première F1 (xo, yo, 6t) et deuxième F2(xo, yo, 5t) valeurs, le calculateur (10) étant adapté pour définir un 15 voisinage R(xo ; yo) formant une partie dudit ensemble de zones, ledit voisinage comprenant au moins la zone considérée (xo ; yo) et des zones voisines (xi ; yi) distinctes de la zone considérée, le calculateur (10) étant adapté pour estimer un 20 critère de détection dans le voisinage, au moins à partir des mesures Fo(x, y, 5t) pour des zones voisines (x1 ; Yi) comprises dans ledit voisinage R(xo ; yo) au cours de la période d'observation T, et le calculateur (10) étant adapté pour estimer la 25 valeur Fe(xo, yo, ôt) du nombre de signaux en fonction en outre au moins dudit critère de détection dans le voisinage. 17.Installation d'imagerie selon la revendication 16, comprenant en outre une unité de détection (5) adaptée 30 pour détecter des signaux correspondants à une émission de rayonnement en provenance de ladite zone pendant l'instant 6t, et pour stocker la mesure correspondante Fo(x, y, 6t) dans ladite mémoire (9). 18.Installation d'imagerie selon la revendication 35 17 dans laquelle l'unité de détection (5) est adaptée pour 2909444 19 détecter des signaux optiques.