Brevet FR2930358A1 — publié le 20091023

Procédé et dispositif d'établissement de carte représentative de zones de risque de prolifération d'une maladie infectieuse, carte et support d'enregistrement associés La présente invention concerne un procédé et un dispositif d'établissement de cartes représentatives de zones de risque de prolifération d'une maladie. Il est connu d'établir des cartes de prévisions météorologiques à partir du traitement d'images satellite. Toutefois, il n'est pas connu de procédé ou de dispositif permettant d'obtenir des cartes représentatives de zones de risque de prolifération d'une maladie infectieuse à partir du traitement d'images satellite. La présente invention a pour but de proposer un procédé ou un dispositif permettant d'obtenir des cartes représentatives de zones de risque de prolifération d'une maladie infectieuse. A cet effet, l'invention a pour objet un procédé d'établissement d'une carte représentant un lieu géographique et délimitant sur ce lieu géographique des zones de risque de prolifération d'une maladie infectieuse (M) provoquée par un agent pathogène, caractérisé en ce que le procédé comporte les étapes suivantes : a) définir à partir d'informations environnementales et climatiques recueillies in situ en fonction de l'agent pathogène, au moins un facteur (F, F') de développement et de transmission de l'agent pathogène, appelé facteur ; b) acquérir au moins une image satellite du lieu géographique à une échelle prédéfinie ; c) traiter la ou chaque image acquise pour générer au moins une première et une deuxième images sur chacune desquelles des zones de présence du facteur (F, F') sont extraites ; d) déterminer les zones communes à la fois aux zones extraites de la première image et aux zones extraites de la deuxième image ; et e) établir une carte représentant lesdites zones communes, lesdites zones communes étant des zones de risque de prolifération de la maladie infectieuse (M). Suivant des modes particuliers de réalisation, le procédé comporte l'une ou plusieurs des caractéristiques suivantes prises isolément ou en combinaison : - à l'étape d) le procédé comporte en outre les étapes suivantes : - obtention d'une image représentant le lieu géographique à l'échelle prédéfinie et délimitant des zones de présence d'êtres vivants susceptibles d'être infectés par la maladie infectieuse (M), ladite image étant appelée image de peuplement ; - détermination de zones communes à la fois aux zones extraites de la première image, aux zones extraites de la deuxième image et aux zones de présence d'êtres vivants de l'image de peuplement, - à l'étape b) le procédé comporte les étapes suivantes : - définir des caractéristiques d'acquisition d'images satellite à partir des informations environnementales et climatiques recueillies in situ en fonction de l'agent pathogène ; - sélectionner dans une base de données au moins un satellite présentant les caractéristiques d'acquisition définies, ladite base de données comportant les caractéristiques d'acquisition d'images de télédétection de plusieurs satellites ; et - acquérir au moins une image satellite du lieu géographique à l'aide du ou de chaque satellite sélectionné, - à l'étape d) le procédé peut également comprendre en outre pour chaque facteur défini (F, F'), les étapes suivantes : dl) sélectionner au moins un traitement (P1, P2) faisant apparaître le facteur défini (F, F') ; d2) traiter la ou chaque image acquise par le ou chaque traitement sélectionné (P1, P2) pour obtenir au moins une première et une deuxième images traitées ; d3) déterminer au moins une valeur seuil (S) du facteur défini (F, F') à partir de laquelle l'agent pathogène est susceptible de se développer ; d4) extraire de ou de chaque image traitée les pixels ayant un compte numérique supérieur ou inférieur à la valeur seuil déterminée (S), lesdits pixels formant des zones extraites, - à l'étape d3) le procédé comprend une étape de détermination d'une première (Scmin Samin) et d'une deuxième (Scmax Sc') valeurs seuil du facteur de défini (F, F') ; dans lequel l'étape d4) comprend une étape d'extraction des pixels ayant un compte numérique supérieur à la première valeur seuil (Scmin Samin) lesdits pixels formant des premières zones de ou de chaque image traitée ; et une étape d'extraction des pixels ayant un compte numérique compris entre la première (Scmin Samin) et la deuxième (Scmax Samax) valeurs seuil, lesdits pixels formant des deuxièmes zones de ou de chaque image traitée ; et dans lequel l'étape e) comprend une étape d'établissement d'une carte représentant les pixels appartenant aux premières zones de la première et de la deuxième images traitées, lesdits pixels formant des zones de présence manifeste d'un agent pathogène, ladite carte représentant les pixels appartenant aux deuxièmes zones de la première et de la deuxième images traitées, lesdits pixels formant des zones de présence potentielle de l'agent pathogène, - le procédé comprend en outre les étapes suivantes : i) calculer le nombre (A) de pixels des premières zones et le nombre (B) de pixels des deuxièmes zones d'une image traitée définie ; j) calculer la moyennne (C) entre le nombre (A) de pixels des premières zones et le nombre (B) de pixels des deuxièmes zones ; k) calculer l'écart (D) entre le nombre (A) de pixels des premières zones et le nombre (B) de pixels des deuxièmes zones ; I) calculer le coefficient d'incertitude (I) de l'image traitée définie, ledit coefficient d'incertitude étant le rapport entre d'une part l'écart (D) divisé par deux et d'autre part la moyenne (C) ; m) répéter les étapes i) à I) pour chaque image traitée et combiner les coefficients d'incertitude (I) calculés pour chaque image traitée pour obtenir le coefficient d'incertitude de la carte représentant les zones de risque de prolifération de la maladie infectieuse (M) , - les caractéristiques d'acquisition d'images satellite comprennent au moins une caractéristique parmi la fréquence d'acquisition des images satellite du lieu géographique, l'échelle spatiale et la résolution spatiale et spectrale des images satellite acquises, - l'étape de définition d'au moins un facteur (F, F') comprend les étapes suivantes : - déterminer le type de transmission de la maladie infectieuse (M) parmi les types suivants : maladies à transmission par l'air, maladies à transmission hydrique et maladies à transmission par vecteurs ; - déterminer des informations environnementales et climatiques en fonction du type de la maladie infectieuse (M) déterminé ; et - définir le ou chaque facteur (F, F') en fonction desdites informations environnementales et climatiques, - une étape de correction de chaque image traitée à partir de prélèvements d'informations environnementales prélevées in situ en fonction de l'agent pathogène, - le procédé peut également comprendre les étapes suivantes : - rechercher la densité (D) d'agent pathogène présente dans les zones communes de la carte, - marquer les zones communes de la carte en fonction d'une échelle de marquage préétablie en fonction de la densité (D) d'agent pathogène, ledit marquage représentant la gravité du risque de prolifération de la maladie infectieuse (M), - à l'étape b) est réalisée à un instant courant t, le procédé comprenant les étapes suivantes, pour chaque facteur (F, F') : - définir un graphe d'évolution temporelle du facteur (F, F') ; le graphe comportant une première valeur (F1) du facteur à l'instant courant t et une deuxième valeur (F2) du facteur à un instant ultérieur t+1 - modifier la taille des zones extraites de la ou de chaque image proportionnellement au rapport de la différence entre la première (F1) et la deuxième (F2) valeurs du facteur, et de la variation de temps entre l'instant courant t et l'instant ultérieur t+1 ; et - établir au moins une image prédictive du facteur représentant les zones modifiées, lesdites zones modifiées représentant les zones de présence du facteur à l'instant ultérieur t+1, - l'image de peuplement a été établie à un instant courant t, le procédé comportant en outre les étapes suivantes : - définir un graphe d'évolution temporelle de la répartition des êtres vivants susceptibles d'être infectés par la maladie infectieuse (M) ; le graphe comportant une première (R1) et une deuxième (R2) valeurs de répartition à l'instant courant t et respectivement à un instant ultérieur t+1 ; - modifier la taille des zones de l'image de peuplement proportionnellement au rapport de la différence entre la première (R1) et la deuxième (R2) valeurs de répartition, et la variation de temps entre l'instant courant t et l'instant ultérieur t+1 ; et - établir au moins une image prédictive de peuplement représentant lesdites zones modifiées, lesdites zones modifiées représentant les zones de répartition des êtres vivants à l'instant ultérieur t+1, - une étape de détermination des zones communes entre les zones modifiées de l'image prédictive du facteur et les lesdites zones modifiées de l'image prédictive de peuplement, l'instant courant t et l'instant ultérieur t+1 considérés dans le graphe d'évolution temporelle du facteur (F, F') étant sensiblement identique à l'instant courant t et l'instant ultérieur t+1 considérés dans le graphe d'évolution temporelle de la répartition des êtres vivants, - la ou chaque image a été acquise au cours de l'étape b) à un instant courant t, le procédé comprenant les étapes suivantes : - définir un graphe d'évolution temporelle de la densité d'agent pathogène (D) ; le graphe (96) comportant une valeur (D1) de densité à un instant ultérieur t+1 ; - créer une image prédictive représentative de la gravité du risque de prolifération de la maladie infectieuse (M) à l'instant ultérieur t+1 en marquant les zones de présence du facteur (F, F') de la ou chaque image en fonction du marquage correspondant à la valeur (D1) de densité de l'échelle de marquage préétablie. L'invention concerne également un support d'enregistrement d'informations qui comporte des instructions pour l'exécution d'un procédé lorsque ces instructions sont exécutées par un calculateur électronique. L'invention concerne également une carte représentant sur un lieu géographique des zones de risque de prolifération d'une maladie infectieuse (M) provoquée par un agent pathogène, ladite carte étant obtenue à partir d'un procédé conforme au procédé de l'invention. L'invention concerne également un dispositif d'établissement d'une carte représentant un lieu géographique et délimitant sur ce lieu géographique des zones de risque de prolifération d'une maladie infectieuse (M) provoquée par un agent pathogène, le dispositif comporte : - une unité de définition d'au moins un facteur (F, F') de développement et de transmission de l'agent pathogène à partir d'informations environnementales et climatiques recueillies in situ en fonction de l'agent pathogène ; - une unité d'acquisition d'au moins une image satellite du lieu géographique à une échelle prédéfinie ; - une unité de traitement de la ou chaque image acquise pour générer au moins une première et une deuxième images sur lesquelles des zones de présence d'un facteur de développement (F, F') sont extraites ; - une unité de détermination de zones communes à la fois aux zones extraites de la première image et aux zones extraites de la deuxième image, lesdites zones communes étant des zones de risque de prolifération de la maladie infectieuse (M) ; l'unité de détermination étant propre à établir une carte représentant lesdites zones communes. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés sur lesquels : - la figure 1 est un schéma simplifié du dispositif d'établissement selon l'invention ; - les figures 2A, 2B sont des diagrammes du procédé d'établissement de cartes selon l'invention ; - les figures 3A et 3B sont des diagrammes d'une variante de mise en oeuvre d'une partie du procédé d'établissement selon l'invention ; - les figures 4A et 4b des diagrammes d'une variante de mise en oeuvre d'une partie du procédé d'établissement selon l'invention ; et - la figure 5 est un schéma simplifié illustrant le diagramme des figures 4A et 4B. En référence à la figure 1, le dispositif 2 d'établissement de cartes selon l'invention comprend des moyens de calcul 4 de type ordinateur, des bases de données 6, 7, 8, 9 reliées à l'ordinateur 4, et des moyens 10 d'acquisition d'images satellite propres à communiquer avec l'ordinateur 4. La base de données 6 comprend la liste des satellites actuellement en fonctionnement, ainsi que les caractéristiques d'acquisition des images satellite de chacun de ces satellites. Les satellites sont par exemple les satellites MODIS, MERIS, JASON, la série SPOT, NOAA-AVHRR, Terra-SARX, Quickbird, Ikonos. Les caractéristiques d'acquisition d'images satellite comprennent notamment la localisation géographique des lieux géographiques imagés par les différents capteurs embarqués à bord des satellites, la fréquence d'acquisition d'images satellite de ces lieux, l'échelle spatiale, la résolution spatiale et spectrale des images satellite acquises. La résolution spatiale correspond à la taille des plus petits objets terrestres pouvant être détectés par les capteurs embarqués. Elle peut donc être assimilée à la distance minimale devant exister entre deux objets au sol pour être correctement distingués. La résolution spectrale correspond à l'aptitude d'un capteur de détection à distinguer deux longueurs d'onde voisines. Elle se caractérise donc par le nombre et la largeur des bandes spectrales (c'est-à-dire les longueurs d'onde) dans lesquelles travaille le capteur embarqué à bord du satellite. L'échelle spatiale est la zone de couverture imagée par les capteurs embarqués. La base de données 7 comprend des informations environnementales et climatiques recueillies in situ en fonction du type de transmission de la maladie infectieuse M. La base de données 8 comprend des procédés de traitement des images satellite propres à mettre en évidence un facteur F, F' de développement et de transmission de l'agent pathogène. Ces procédés comprennent par exemple le procédé SST4, le procédé OC3 de l'Administration Nationale de l'Aéronautique et de l'Espace plus connue sous le nom de NASA, ou le procédé ALGAE1 et ALGAE2, le procédé faisant apparaître la matière suspendue et la � Yellow � substance de l'Agence Spatiale Européenne plus connue sous le nom de ESA, le procédé 005 de l'Institut Français de Recherche pour l'Exploitation de la MER plus connu sous le nom de l'IFREMER, le procédé de détection des mares temporaires (gîtes larvaires) du Centre National d'Etudes Spatiales plus connu sous le nom du CNES et du groupement d'intérêt public MEDIAS-France. Ces procédés de traitement sont fournis avec des correspondances entre les comptes numériques de l'image traitée et le facteur de développement mis en évidence par le traitement. Par exemple, un compte numérique de 50 correspond à une température de 5° C. Les procédés de traitement ne font pas uniquement apparaître des facteurs de développement représentatifs d'un état tels que la température ou la concentration en chlorophylle mais également des facteurs de développement tels qu'une couleur particulière permettant de mettre en évidence un milieu favorable au développement de l'agent pathogène. La base de données 9 comprend des images représentant des zones de présence d'êtres vivants à différentes échelles spatiales. Ces images sont appelées ci-après images de peuplement. Les images de cette base de données 9 ont été préalablement établies en fonction de la distribution des êtres vivants susceptibles d'attraper la maladie, de la caractérisation de leurs types d'habitat, de l'immunité de la population et du modèle épidémiologique de la maladie. Selon l'invention, un être vivant est un humain ou un animal susceptible d'être infecté par l'agent pathogène de la maladie infectieuse M. On parle alors d'hôtes. La base de données 9 est régulièrement mise à jour. En référence à la figure 1, les moyens d'acquisition d'images 10 comprennent des capteurs embarqués sur des satellites 24 et propres à imager un lieu géographique 26 et à transmettre l'image satellite résultante à l'ordinateur 4, via l'unité de communication 16. Pour des raisons de simplification, un seul satellite a été représenté sur la figure 1. L'ordinateur 4 comprend une interface homme-machine 12, une unité de sélection 14 reliée à l'interface 12 et une unité de communication 16 reliée à l'unité de sélection 14 et à l'interface 12. L'unité de communication 16 est apte à échanger des données avec les moyens 10 d'acquisition d'images satellite. L'ordinateur 4 comprend en outre une unité de traitement 18 reliée à l'unité de communication 16 et à l'unité de sélection 14, et une unité de détermination 20 reliée à l'unité de traitement 18. L'unité de détermination 20 est propre à générer une carte 22 représentative de zones 23 de risque de prolifération d'une maladie infectieuse M. L'ordinateur 4 est également propre à recevoir des informations environnementales et climatiques prélevées sur le lieu 26 à l'aide d'un dispositif de prélèvement 28. Le procédé d'établissement de cartes selon l'invention est décrit ci-dessous dans un exemple de réalisation de cartes représentatives de zones de risque de prolifération du choléra. Toutefois, le procédé peut être appliqué à d'autres maladies, telles que par exemple la fièvre de la vallée du Rift, le paludisme, la dengue, le chikungunya, les maladies respiratoires liées à la pollution, etc. En référence à la figure 2A, le procédé débute par une étape 30 de détermination de l'agent pathogène. Dans le cas de la maladie du choléra, cet agent est le Vibrio cholerae. Au cours d'une étape 31, le type de transmission de la maladie est déterminé. La maladie infectieuse M peut être une maladie à transmission par l'air, une maladie à transmission hydrique ou une maladie à transmission par vecteurs. Par exemple, le choléra est une maladie à transmission hydrique. Au cours d'une étape 32, des informations environnementales et climatiques sont déterminées en fonction du type de transmission de la maladie infectieuse M. Par exemple, pour une maladie infectieuse M transmise par l'air, les informations environnementales et climatiques concernent les données sur les concentrations en aérosols et en gaz des constituants atmosphériques, les données climatiques et des informations de transmissibilité entre un être vivant susceptible d'être infecté par l'agent pathogène de la maladie infectieuse M, aussi appelées données de contact hôte-agent. Lorsque la maladie infectieuse M est à transmission hydrique, les informations environnementales et climatiques comprennent des informations sur l'écologie des agents infectieux (vibrion, virus, bactérie, parasite...), des données microbiologiques, des données environnementales, des données climatiques et des données de contact hôte-agent. Lorsque la maladie infectieuse M est une maladie à transmission par vecteurs, les informations environnementales et climatiques concernent des données relatives aux vecteurs, c'est-à-dire aux moustiques, aux tiques, aux rongeurs, aux oiseaux, etc. Ces données comprennent des données de biologie animale, des données entomologiques, des données environnementales, des données climatiques, des données microbiologiques et des données de contact hôte-agent vecteurs. Dans l'exemple de la maladie du choléra, les données microbiologiques déterminées comprennent la mesure journalière in situ de la densité de Vibrio cholerae sur plusieurs sites expérimentaux, et la mesure journalière in situ de la concentration en chlorophylle-a, phytoplanctons et zooplanctons. Les données environnementales comprennent des mesures journalières in situ de la salinité, de l'acidité et de la composition chimique de l'eau. Les données climatiques concernent la mesure journalière in situ de la température de l'eau. Les données de contact hôte-agent comprennent des données de suivi de courantologie des océans. Les informations environnementales et climatiques sont recueillies in situ, puis stockées dans la base de données 7. Au cours d'une étape 33, des facteurs de développement et de transmission F, F' de l'agent pathogène sont déterminés à partir des informations environnementales et climatiques recueillies in situ et stockées dans la base de données 7. Selon l'invention, un facteur de développement et de transmission est une caractéristique du milieu terrestre, atmosphérique ou climatique favorisant l'apparition et la croissance d'un agent pathogène. Dans l'exemple du choléra, ces facteurs F, F' comprennent la température de la surface de la mer, la concentration en chlorophylle-a de la mer, la courantologie de la mer et la concentration en algues. Pour d'autres maladies, les facteurs de développement sont par exemple la surface d'eaux stagnantes, le battement des mares temporaires, la température de l'air, le taux d'humidité, la pluviométrie, la surface de zones en végétation, etc. Au cours d'une étape 34, les caractéristiques d'acquisition d'images satellite sont définies en fonction des informations environnementales et climatiques et des facteurs de développement et de transmission de l'agent pathogène. Ces informations permettent de définir l'échelle temporelle, l'échelle spatiale et la résolution spatiale adéquates. L'échelle temporelle permet de déterminer la fréquence d'acquisition des images satellite. Comme vu précédemment, ces caractéristiques comprennent la localisation des lieux géographiques imagés par les capteurs des satellites, la fréquence d'acquisition des images satellite et la résolution spatiale et spectrale. Dans l'exemple du choléra, les caractéristiques d'acquisition des images satellite comprennent une résolution spatiale de 1 Km sur une portion du lieu 26 ayant une superficie de 100 Km à 1000 Km, avec une fréquence d'acquisition des images de journalière à hebdomadaire, et des longueurs d'onde comprises entre le bleu et l'infrarouge thermique. Au cours d'une étape 35, l'unité de sélection 14 sélectionne dans la base de données 6 un ou plusieurs satellites répondant aux caractéristiques d'acquisition d'images satellite définies au cours de l'étape 34. Dans l'exemple du choléra, les satellites sélectionnés à partir de ces caractéristiques d'acquisition sont les satellites MODIS, MERIS et JASON. Le satellite MODIS est apte à transmettre de façon journalière des images satellite du lieu 26 avec une résolution spatiale de 1 Km. Le satellite MERIS est apte à transmettre des images satellite de façon journalière du lieu 26 avec une résolution spatiale de 300 m. Le satellite JASON est apte à fournir de façon journalière des images de courantologie du lieu 26 avec une résolution spatiale de 1 /8 de degré. L'unité de sélection 14 transmet alors via l'unité de communication 16 une requête d'images satellite aux moyens 10 d'acquisition d'images satellite. Au cours d'une étape 36, les satellites sélectionnés 24 prennent l'image satellite du lieu 26 et la transmettent à l'unité de traitement 18 via l'unité de communication 16. Au cours d'une étape 37, l'unité de sélection 14 sélectionne plusieurs procédés de traitement de chaque image satellite réceptionnée par l'unité de communication 16. Les procédés de traitement sont déterminés de façon à faire apparaître des zones de présence d'un facteur de développement de l'agent pathogène. Au cours d'une étape 38, les images satellite acquises sont traitées par les procédés de traitement sélectionnés. L'image satellite générée par le satellite MERIS est traitée par des algorithmes standards de l'ESA pour obtenir des images représentant les concentrations en algues et en matière suspendue ainsi que la présence des substances jaunes appelées � Yellow � substance. Les données acquises par le satellite JASON sont traitées pour obtenir des images de courantologie. L'image satellite 39 acquise par le satellite MODIS est traitée par l'algorithme SST4 de la NASA pour obtenir des images représentatives de températures de surface de la mer. Cette image satellite 39 est traitée par le procédé de traitement IFREMER 005 pour obtenir une image représentant des concentrations en chlorophylle-a. Cette image satellite 39 est également traitée par le procédé de traitement NASA OC3 pour obtenir une image 40 représentant des concentrations en chlorophylle-a, illustrée sur la figure 2B. L'image satellite 39 est également traitée par les procédés de traitement ESA ALGAE2 pour obtenir une image 42 faisant apparaître des concentrations en algues, illustrée sur la figure 2B. Pour des raisons de simplification, seul le post-traitement des images 40 et 42, est décrit ci- dessous, mais l'ensemble des images traitées est post-traité de cette façon. Au cours d'une étape 44, un utilisateur définit pour chaque facteur de développement et de transmission au moins une valeur seuil S d'une quantité de ce facteur. Cette valeur seuil S définit la quantité du facteur de développement à partir de laquelle l'agent pathogène est susceptible de se développer. Cette saisie est réalisée à l'aide de l'interface homme-machine 12. Un facteur de développement F, F' a été mis en évidence sur chaque image traitée 40, 42. Une valeur seuil de ce facteur de développement F, F' est saisie pour le traitement de cette image 40, 42. Dans l'exemple décrit ci-dessus, le facteur de développement F mis en évidence sur l'image traitée 40 concerne la concentration en chlorophylle-a et le facteur de développement F' mis en évidence sur l'image traitée 40 concerne la concentration en algues. Une valeur seuil Sc d'une concentration en chlorophylle-a F de 0.7 mg/m3 est déterminée et saisie, via l'interface homme-machine 12, pour le post traitement de l'image 40, et une valeur seuil Sa d'une concentration en algues F' de 0.9 mg/m3 est déterminée et saisie pour le post traitement de l'image 42. Au cours d'une étape 45, l'unité de traitement 18 extrait les zones 46 de concentrations en chlorophylle-a par recherche des pixels ayant un compte numérique supérieur à la valeur seuil Sc saisie pour le facteur chlorophylle-a. De même, elle extrait les zones 47 de présence d'algues par recherche des pixels ayant un compte numérique supérieur à la valeur seuil Sa saisie pour le facteur algues. Le compte numérique est la mesure physique faite par un capteur embarqué sur le satellite sur un pixel dans un domaine de longueurs d'onde donné. Au cours d'une étape 51, l'unité de détermination 20 extrait des zones communes 52 aux zones 46 de forte concentration en chlorophylle-a, mises en évidence sur l'image 40, et aux zones 47 de forte concentration en algues, mises en évidence par le traitement ALGAE2 sur l'image 42. A cet effet, l'unité de détermination 20 établit une image 53 représentant les pixels appartenant à la fois aux zones 46 extraites de l'image 40 et aux zones 47 extraites de l'image 42. Ces pixels forment des zones 52 de présence potentielle de l'agent pathogène. Au cours d'une étape 55, une image 56 de peuplement est sélectionnée par l'unité de sélection 14 dans la base de données 9 de manière à ce que cette image 56 de peuplement représente le même lieu géographique que le lieu géographique 26 imagé au cours de l'étape 36 et de manière à ce que l'image 56 soit à la même échelle et à la même résolution spatiale que celles définies au cours de l'étape 34. Des zones 58 de présence d'êtres vivants susceptibles d'attraper la maladie M sont délimitées sur cette image de peuplement 56. En variante, cette image de peuplement 56 est obtenue par traitement d'une image satellite acquise par le satellite 24 à l'échelle et à la résolution spatiales déterminées au cours de l'étape 34. Dans ce cas, le traitement choisi permet de mettre en évidence les êtres vivants susceptibles d'attraper la maladie. Comme visible sur la figure 2B, l'image de peuplement 56 du lieu 26 comprend deux zones 58 de présence d'êtres vivants susceptibles d'attraper le choléra. Au cours d'une étape 59, l'unité de détermination 20 détermine des zones communes 23 à la fois aux zones 58 de présence d'êtres vivants de l'image de peuplement 56 et aux zones 52 de présence potentielle de l'agent pathogène de l'image 53 générée au cours de l'étape 51 pour obtenir une carte 22 représentative de zones 23 de risque de prolifération du choléra, comme illustré sur la figure 2B. Selon l'invention, l'ordinateur 4 établit journalièrement de nouvelles cartes représentatives de zones de risque à partir des images satellite qu'il reçoit chaque jour du satellite 24. En variante, l'unité de traitement 18 est propre à corriger les images satellite reçues en fonction d'informations environnementales prélevées sur le lieu 26 et transmises par le dispositif de prélèvement 28. Ces informations environnementales comprennent par exemple des prélèvements d'eau. En variante, le procédé s'arrête à l'étape 51. Dans ce cas, l'image 53 constitue une carte sur laquelle les zones 52 sont des zones du risque de prolifération de la maladie infectieuse M. Selon un deuxième mode de réalisation de l'invention et dans le but d'établir une incertitude, le procédé selon l'invention comprend à la place des étapes 44 à 51, les étapes 62 à 77 illustrées sur les figures 3A et 3B. A cet effet, l'étape de traitement 38 est suivie par une étape 62 de détermination d'une valeur seuil minimale Scmin et d'une valeur seuil maximale Sc' du facteur de développement F mis en évidence sur l'image traitée 40. Cette étape 62 est suivie d'une étape 64 au cours de laquelle l'unité de traitement 18 détermine les pixels ayant un compte numérique supérieur à la valeur seuil minimale Scmin lesdits pixels formant des premières zones 66 de l'image satellite traitée 40, et les pixels ayant un compte numérique compris entre les valeurs seuil minimale Scmin et maximale Samax lesdits pixels formant des deuxième zones 67 de l'image satellite traitée. Les étapes 62 et 64 sont également mises en oeuvre pour la deuxième image 42. A cet effet, une valeur seuil minimale Samin et une valeur seuil maximale Samax sont déterminées pour le facteur de développement F' mis en évidence sur l'image 42. Les pixels ayant un compte numérique supérieur à la valeur seuil minimale Samin et les pixels ayant un compte numérique compris entre les valeurs seuil minimale Samin et maximale Samax sont déterminés. Ces pixels forment des premières 68 et respectivement des deuxièmes 69 zones de l'image 42. Au cours d'une étape 65, l'unité de détermination 20 établit une image 70 représentant les pixels appartenant à la fois aux premières zones 66 de la première image 40 et aux premières zones 68 de la deuxième image 42, lesdits pixels formant des zones 71 de présence très probable d'un agent pathogène et donc de risque élevé. Sur cette image 70, l'unité de détermination 20 détermine les pixels appartenant à la fois aux premières 66 et aux deuxièmes 67 zones de la première image 40, et aux deuxièmes zones 69 de la deuxième image 42, et les pixels appartenant à la fois aux premières 68 et aux deuxièmes 69 zones de la deuxième image 42 et aux deuxièmes zones 67 de la première image 40 ; lesdits pixels formant des zones 72 de présence moins probable de l'agent pathogène et donc de risque faible. L'image 70 représente donc des zones de risque élevé, matérialisées par les lignes 71, et des zones de risque plus faible matérialisées par les lignes 72. Il existe donc un risque élevé d'être infecté par la maladie M dans les zones 71, et un risque faible dans les zones 72. Puis, au cours d'une étape 73, le nombre A de pixels des premières zones 66 et le nombre B de pixels des deuxièmes zones 67 de la première image 40 sont comptabilisés. Au cours d'une étape 74, l'unité de détermination 20 détermine la moyenne C du nombre de pixels appartenant aux premières 66 et aux deuxièmes 67 zones selon la relation C = A+B 2 Au cours d'une étape 75, l'unité de détermination 20 calcule l'écart D entre le nombre A de pixels des premières zones 66 et le nombre B de pixels des deuxièmes zones 67, puis calcule au cours d'une étape 76 le rapport entre l'écart D calculé au cours de l'étape 75, divisé par 2 et le nombre moyen C calculé à l'étape 74. Ce rapport définit un coefficient d'incertitude I de la première image traitée 40. I = (D/2)/C Les étapes 73 à 76 sont mises en oeuvres pour la deuxième image traitée 42. Puis, l'unité de détermination 20 combine au cours d'une étape 77 l'ensemble des coefficients d'incertitude I estimés pour chaque image traitée 40, 42. La combinaison de ces coefficients d'incertitude de chaque image traitée correspond au coefficient d'incertitude global de l'image 70 des zones de risque de prolifération d'une maladie. Cette combinaison consiste par exemple en une somme de ces coefficients, une somme pondérée de ceux-ci, un rapport, etc. Les étapes 55 et 59 sont appliquées à l'image 70 pour obtenir une nouvelle carte présentant deux niveaux de risque et tenant compte des zones de présence des êtres vivants. En variante, le procédé s'arrête à l'étape 65. Dans ce cas, l'image 70 constitue une carte sur laquelle les zones 71, 72 sont des zones de risque de prolifération de la maladie infectieuse M. En variante, les zones 66 et 67 sont déterminées par photo interprétation, c'est-à-dire par délimitation visuelle de zones sans faire une discrimination automatique à partir des valeurs seuil S ou Scmin et Scmax En référence à la figure 4A, une variante du procédé d'établissement selon l'invention est illustrée. Cette variante permet d'obtenir des images prédictives 122, 222 des zones de risque de prolifération de la maladie infectieuse M. Selon l'invention, une image prédictive est une image de prévision de l'évolution des zones à un instant ultérieur à l'instant d'acquisition des images satellite. Au cours d'une étape 78, un graphe d'évolution temporelle 80 du facteur de développement est défini à partir des informations environnementales et climatiques recueillies in situ et d'un modèle d'évolution prédéfini et qui est par exemple stocké dans la base de données 7. Le modèle d'évolution permet d'appréhender l'évolution spatiale du facteur de développement au cours du temps. Au cours d'une étape 82, la taille des zones extraites 46 de la première image 40 au cours de l'étape 45, est modifiée en fonction du graphe d'évolution 80. L'étape 82 est réalisée pour une première image 40, acquise au cours de l'étape 36 à un instant courant t. La taille des zones extraites 46 est modifiée proportionnellement au rapport de la différence entre une première F1 et une deuxième F2 valeurs du facteur, et de la variation de temps entre l'instant courant t et un instant ultérieur t+1. Cette modification est réalisée un nombre de fois correspondant au nombre d'images prédictives que l'on souhaite obtenir. Par exemple, sur la figure 4A, deux images prédictives 140, 240 à des instants ultérieurs t+1 et t+2 ont été établies à partir de la première image 40 acquise à l'instant courant t. Le graphe d'évolution temporelle 80 choisi pour réaliser l'étape 82 correspond au graphe d'évolution du facteur qui a été mis en évidence lors du traitement de l'image à laquelle ce graphe sera appliqué. Par exemple, l'image satellite 39 a été traitée par le procédé de traitement P1 pour faire apparaître le facteur de développement F. Au cours de l'étape 82, le graphe d'évolution du facteur F est utilisé pour modifier la taille des zones extraites 46 de la première image 40. De même, l'image satellite 39 a été traitée par le procédé de traitement P2 pour faire apparaître le facteur de développement F'. Au cours de l'étape 82, le graphe d'évolution du facteur F' est utilisé pour modifier la taille des zones extraites 47 de la deuxième image 42. Au cours d'une étape 84, les étapes 78 et 82 sont répétées pour chaque image traitée. Au cours d'une étape 86, un graphe 88 d'évolution temporelle de la répartition des êtres vivants est défini. Ce graphe 88 est par exemple issu de la base de données 9. Au cours d'une étape 90, la taille des zones 58 de l'image de peuplement 56 est modifiée en fonction du graphe 88. La modification de la taille des zones 58 est proportionnelle au rapport de la différence entre une première valeur de répartition R1 définie à l'instant courant t et une deuxième valeur de répartition R2 définie à l'instant ultérieur t+1 divisée par la variation de temps entre l'instant courant t et l'instant ultérieur t+1. L'image de peuplement 56 a été établie à un instant courant t qui correspond sensiblement à l'instant t d'acquisition de l'image satellite 39. L'instant ultérieur t+1 du graphe d'évolution 88 correspond sensiblement à l'instant t+1 du graphe d'évolution 80. Les images de peuplement 156 et 256 créées à l'instant t+1 et respectivement à l'instant t+2 sont des images prédictives de peuplement. Au cours d'une étape 92, les zones communes 123, 223 à la fois aux zones 146 de l'image prédictive 140 à l'instant ultérieur t+1 et aux zones 158 de l'image prédictive 156 de peuplement au même instant ultérieur t+1 sont déterminées. Cette étape est réalisée pour chaque instant t+1, t+2 pour lesquelles des images prédictives ont été réalisées à des instants correspondants sensiblement. Au cours d'une étape 94, un graphe 96 d'évolution temporelle de la densité d'agent pathogène D est défini, par exemple à partir des informations climatiques et environnementales déterminées au cours de l'étape 82 ou à partir d'un modèle d'évolution de la densité de l'agent pathogène. Ce graphe 96 est par exemple stocké dans la base de données 7. Au cours d'une étape 98, les zones 23, 123, 223 de risque de prolifération de la maladie infectieuse M de la carte 22, 122, 222 établie à l'instant t, à l'instant t+1 et respectivement à l'instant t+2 sont marquées en fonction de la densité d'agent pathogène. Ce marquage est par exemple réalisé en fonction d'une échelle 100 de coloris préétablis. Ainsi, le procédé permet d'obtenir une carte 22 représentative de la gravité du risque à un instant courant t et des cartes prédictives 122, 222 de la gravité du risque aux instants ultérieurs t+1 et t+2. En variante, l'étape de marquage 98 est réalisée après l'étape de détermination 51 pour une seule image 53 acquise à l'instant t par marquage des zones communes 52 en fonction d'une information de densité d'agent pathogène présent dans les zones communes et de l'échelle 100 de marquage préétablie en fonction de la densité d'agent pathogène. En variante, l'étape de marquage 98 est réalisée après l'étape de modification 82 sur les images prédictives 140, 240 obtenues par cette étape de modification. L'invention concerne également une carte obtenue par ce procédé, ainsi qu'un support d'enregistrement d'information comportant des instructions pour l'exécution de ce procédé. REVENDICATIONS1.- Procédé d'établissement d'une carte (22, 53, 70) représentant un lieu géographique (26) et délimitant sur ce lieu géographique des zones (23, 52, 71, 72) de risque de prolifération d'une maladie infectieuse (M) provoquée par un agent pathogène, caractérisé en ce que le procédé comporte les étapes suivantes : a) définir (31, 32, 33) à partir d'informations environnementales et climatiques recueillies in situ en fonction de l'agent pathogène, au moins un facteur (F, F') de développement et de transmission de l'agent pathogène, appelé facteur ; b) acquérir (34, 35, 36) au moins une image satellite (39) du lieu géographique (26) à une échelle prédéfinie ; c) traiter (37, 38, 44, 45, 62, 64) la ou chaque image acquise (39) pour générer au moins une première (40) et une deuxième (42) images sur chacune desquelles des zones de présence du facteur (F, F') sont extraites ; d) déterminer (51, 59, 65) les zones communes (23, 52, 71, 72) à la fois aux zones (46, 66, 67) extraites de la première image (40) et aux zones (47, 68, 69) extraites de la deuxième image (42) ; et e) établir (51, 59, 65) une carte (22, 53, 70) représentant lesdites zones communes (23, 52, 71, 72), lesdites zones communes (23, 52, 71, 72) étant des zones de risque de prolifération de la maladie infectieuse (M). 2.- Procédé selon la revendication 1, dans lequel l'étape d) comporte en outre les étapes suivantes : - obtention d'une image (56) représentant le lieu géographique (26) à l'échelle prédéfinie et délimitant des zones (58) de présence d'êtres vivants susceptibles d'être infectés par la maladie infectieuse (M), ladite image (56) étant appelée image de peuplement ; - détermination (59) de zones communes (23, 52, 71, 72) à la fois aux zones extraites (46, 66, 67) de la première image (40), aux zones extraites (47, 68, 69) de la deuxième image (42) et aux zones (58) de présence d'êtres vivants de l'image de peuplement (56). 3.- Procédé selon l'une quelconque des revendications 1 et 2, dans lequel l'étape b) comporte les étapes suivantes : - définir (34) des caractéristiques d'acquisition d'images satellite à partir des informations environnementales et climatiques recueillies in situ en fonction de l'agent pathogène ; - sélectionner (35) dans une base de données (6) au moins un satellite (24) présentant les caractéristiques d'acquisition définies, ladite base de données (6) comportant les caractéristiques d'acquisition d'images de télédétection de plusieurs satellites ; et - acquérir (36) au moins une image satellite (39) du lieu géographique (26) à l'aide du ou de chaque satellite sélectionné (24). 4.- Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape d) comprend en outre pour chaque facteur défini (F, F'), les étapes suivantes : dl) sélectionner (37) au moins un traitement (P1, P2) faisant apparaître le facteur défini (F, F') ; d2) traiter (38) la ou chaque image acquise (39) par le ou chaque traitement sélectionné (P1, P2) pour obtenir au moins une première (40) et une deuxième (42) images traitées ; d3) déterminer (44, 62) au moins une valeur seuil (S) du facteur défini (F, F') à partir de laquelle l'agent pathogène est susceptible de se développer ; d4) extraire (45, 64) de ou de chaque image traitée (40, 42) les pixels ayant un compte numérique supérieur ou inférieur à la valeur seuil déterminée (S), lesdits pixels formant des zones extraites (46, 47). 5.- Procédé selon la revendication 4, dans lequel l'étape d3) comprend une étape de détermination (62) d'une première (Scmin Samin) et d'une deuxième (Scmax Scmax) valeurs seuil du facteur de défini (F, F') ; dans lequel l'étape d4) comprend une étape d'extraction (64) des pixels ayant un compte numérique supérieur à la première valeur seuil (Scmin Samin) lesdits pixels formant des premières zones (66, 68) de ou de chaque image traitée (40, 42) ; et une étape d'extraction des pixels ayant un compte numérique compris entre la première (Scmin Samin) et la deuxième (Scmax Scmax) valeurs seuil,lesdits pixels formant des deuxièmes zones (67, 69) de ou de chaque image traitée ; et dans lequel l'étape e) comprend une étape d'établissement (65) d'une carte (70) représentant les pixels appartenant aux premières zones (66, 68) de la première (40) et de la deuxième (42) images traitées, lesdits pixels formant des zones (71) de présence manifeste d'un agent pathogène, ladite carte (70) représentant les pixels appartenant aux deuxièmes zones (67, 69) de la première (40) et de la deuxième (42) images traitées, lesdits pixels formant des zones (72) de présence potentielle de l'agent pathogène. 6.- Procédé selon la revendication 5, qui comprend en outre les étapes suivantes : i) calculer (73) le nombre (A) de pixels des premières zones (66, 68) et le nombre (B) de pixels des deuxièmes zones (67, 69) d'une image traitée définie (40, 42) ; j) calculer (74) la moyennne (C) entre le nombre (A) de pixels des premières zones (66, 68) et le nombre (B) de pixels des deuxièmes zones (67, 69) ; k) calculer (75) l'écart (D) entre le nombre (A) de pixels des premières zones (66, 68) et le nombre (B) de pixels des deuxièmes zones (67, 69) ; I) calculer (76) le coefficient d'incertitude (I) de l'image traitée définie (40, 42), ledit coefficient d'incertitude étant le rapport entre d'une part l'écart (D) divisé par deux et d'autre part la moyenne (C) ; m) répéter les étapes i) à I) pour chaque image traitée (40, 42) et combiner (77) les coefficients d'incertitude (I) calculés pour chaque image traitée (40, 42) pour obtenir le coefficient d'incertitude de la carte (22, 53, 70) représentant les zones de risque de prolifération de la maladie infectieuse (M) . 7.- Procédé selon l'une quelconque des revendications 3 à 6, dans lequel lesdites caractéristiques d'acquisition d'images satellite comprennent au moins une caractéristique parmi la fréquence d'acquisition des images satellite du lieu géographique, l'échelle spatiale et la résolution spatiale et spectrale des images satellite acquises. 8.- Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape de définition (33) d'au moins un facteur (F, F') comprend les étapes suivantes :- déterminer (31) le type de transmission de la maladie infectieuse (M) parmi les types suivants : maladies à transmission par l'air, maladies à transmission hydrique et maladies à transmission par vecteurs ; - déterminer (32) des informations environnementales et climatiques en fonction du type de la maladie infectieuse (M) déterminé ; et - définir (33) le ou chaque facteur (F, F') en fonction desdites informations environnementales et climatiques. 9.- Procédé selon l'une quelconque des revendications précédentes, qui comporte en outre une étape de correction de chaque image traitée (40, 42) à partir de prélèvements d'informations environnementales prélevées in situ en fonction de l'agent pathogène. 10.- Procédé selon l'une quelconque des revendications précédentes, dans lequel le procédé comprend les étapes suivantes : - rechercher (94) la densité (D) d'agent pathogène présente dans les zones communes (23, 52, 71, 72) de la carte (22, 53, 70), - marquer (98) les zones communes (23, 52, 71, 72) de la carte (22, 53, 70) en fonction d'une échelle de marquage préétablie en fonction de la densité (D) d'agent pathogène, ledit marquage représentant la gravité du risque de prolifération de la maladie infectieuse (M). 11.- Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape b) est réalisée à un instant courant t, le procédé comprenant les étapes suivantes, pour chaque facteur (F, F') : - définir (78) un graphe (80) d'évolution temporelle du facteur (F, F') ; le graphe (80) comportant une première valeur (F1) du facteur à l'instant courant t et une deuxième valeur (F2) du facteur à un instant ultérieur t+1 - modifier (82) la taille des zones extraites (46, 47) de la ou de chaque image (40, 42) proportionnellement au rapport de la différence entre la première (F1) et la deuxième (F2) valeurs du facteur, et de la variation de temps entre l'instant courant t et l'instant ultérieur t+1 ; et - établir (82) au moins une image prédictive du facteur (140, 240) représentant les zones modifiées (146, 147), lesdites zones modifiées (146, 147) représentant les zones de présence du facteur à l'instant ultérieur t+1. 12.- Procédé selon l'une quelconque des revendications 2 à 11, dans lequel l'image de peuplement (56) a été établie à un instant courant t, le procédé comportant en outre les étapes suivantes : - définir (86) un graphe (88) d'évolution temporelle de la répartition des êtres vivants susceptibles d'être infectés par la maladie infectieuse (M) ; le graphe (88) comportant une première (R1) et une deuxième (R2) valeurs de répartition à l'instant courant t et respectivement à un instant ultérieur t+1 ; - modifier (90) la taille des zones (58) de l'image de peuplement (56) proportionnellement au rapport de la différence entre la première (R1) et la deuxième (R2) valeurs de répartition, et la variation de temps entre l'instant courant t et l'instant ultérieur t+1 ; et - établir (90) au moins une image prédictive de peuplement (156, 256) représentant lesdites zones modifiées (158, 258), lesdites zones modifiées représentant les zones de répartition des êtres vivants à l'instant ultérieur t+1. 13.- Procédé selon la revendication 11 et 12, qui comprend une étape de détermination des zones communes (123, 223) entre les zones modifiées (146, 147) de l'image prédictive du facteur (140, 240) et les lesdites zones modifiées (158, 258) de l'image prédictive de peuplement (156, 256), l'instant courant t et l'instant ultérieur t+1 considérés dans le graphe (80) d'évolution temporelle du facteur (F, F') étant sensiblement identique à l'instant courant t et l'instant ultérieur t+1 considérés dans le graphe (88) d'évolution temporelle de la répartition des êtres vivants. 14.- Procédé selon la revendication 10, dans lequel la ou chaque image a été acquise au cours de l'étape b) à un instant courant t, le procédé comprenant les étapes suivantes : - définir (94) un graphe (96) d'évolution temporelle de la densité d'agent pathogène (D) ; le graphe (96) comportant une valeur (D1) de densité à un instant ultérieur t+1 ; - créer (98) une image prédictive (122, 222) représentative de la gravité du risque de prolifération de la maladie infectieuse (M) à l'instant ultérieur t+1 en marquant les zones de présence du facteur (F, F') de la ou chaque image (40, 42) en fonction du marquage correspondant à la valeur (D1) de densité de l'échelle de marquage préétablie. 15. Support d'enregistrement d'informations, caractérisé en ce qu'il comporte des instructions pour l'exécution d'un procédé conforme à l'une quelconque des revendications 1 à 14, lorsque ces instructions sont exécutées par un calculateur électronique. 16- Carte (22, 53, 70) représentant sur un lieu géographique (26) des zones (23, 52, 71, 72) de risque de prolifération d'une maladie infectieuse (M) provoquée par un agent pathogène, caractérisé en ce que ladite carte (22, 53, 70) étant obtenue à partir d'un procédé conforme à l'une quelconque des revendications 1 à 14. 17.- Dispositif (4) d'établissement d'une carte (22, 53, 70) représentant un lieu géographique (26) et délimitant sur ce lieu géographique des zones (23, 52, 71, 72) de risque de prolifération d'une maladie infectieuse (M) provoquée par un agent pathogène, caractérisé en ce que le dispositif (2) comporte : - une unité de définition (7) d'au moins un facteur (F, F') de développement et de transmission de l'agent pathogène à partir d'informations environnementales et climatiques recueillies in situ en fonction de l'agent pathogène ; - une unité d'acquisition (16, 24) d'au moins une image satellite (39) du lieu géographique (26) à une échelle prédéfinie ; - une unité de traitement (18) de la ou chaque image acquise (39) pour générer au moins une première (40) et une deuxième (42) images sur lesquelles des zones de présence d'un facteur de développement (F, F') sont extraites ; - une unité de détermination (20) de zones communes (23, 52, 71, 72) à la fois aux zones (46, 66, 67) extraites de la première image (40) et aux zones (47, 68, 69) extraites de la deuxième image (42), lesdites zones communes (23, 52, 71, 72) étant des zones de risque de prolifération de la maladie infectieuse (M) ; l'unité de détermination (20) étant propre à établir une carte (22, 53, 70) représentant lesdites zones communes (23, 52, 71, 72).