Brevet FR3004834A1 — publié le 20141024

Procédé de production de données représentatives de contraintes de traitement de radiothérapie, dispositifs et programme correspondant 1. Domaine de l'invention L'invention se rapporte au traitement de données d'imagerie médicale. L'invention se rapporte plus particulièrement à une méthode de traitement de données médicales en vue d'une planification, par une console de planification, d'une ou plusieurs actions à mener, comme par exemple des actions d'application d'un rayonnement sur une ou plusieurs zones données dans le cadre d'une radiothérapie. On rappelle que la radiothérapie est une méthode de traitement de cancers qui utilise des radiations pour détruire les cellules cancéreuses. L'irradiation a pour but de détruire des cellules tumorales tout en épargnant des tissus sains périphériques. Pour permettre une irradiation, un accélérateur linéaire d'électron est utilisé, lequel produit un faisceau d'irradiation. 2. Art antérieur L'histoire de la radiothérapie est relativement récente. La radiothérapie gagne en efficacité en fonction de l'évolution des progrès techniques. Ainsi, grâce à la mise en oeuvre d'un scanner, la planification des traitements de radiothérapie en trois-dimensions devint possible, ce qui représentait une avancée majeure par rapport aux traitements en deux dimensions. Les traitements basés sur l'utilisation du scanner permettent aux radiooncologues de déterminer plus précisément la distribution de la dose de radiation en utilisant des images tomodensitométriques de l'anatomie du patient. Des techniques d'imagerie ont permis l'apparition de la radiothérapie guidée par l'image (IGRT) qui permet de mieux contrôler la position de la zone à traiter au fur et à mesure du traitement. Des innovations au niveau des appareils de traitement comme l'apparition des collimateurs multi-lames ont permis l'apparition de la radiothérapie conformationnelle avec modulation d'intensité (RCMI) qui permet d'adapter plus précisément l'irradiation à la forme des organes à traiter. Il a été ainsi possible de visualiser et de traiter plus efficacement les tumeurs, ce qui s'est traduit par un meilleur pronostic pour les patients, une meilleure préservation des organes sains et moins d'effets secondaires. Il subsiste cependant des problématiques. En effet, bien qu'il soit désormais possible de moduler l'intensité de l'irradiation en fonction de zones de traitement déterminées, il n'en reste pas moins que certaines régions anatomiques sont plus sensibles que d'autres aux irradiations. Ces zones de traitement peuvent être des organes, mais il s'agit également de portions d'organes. Il a été démontré que certaines portions d'organes, qui sont surirradiées, peuvent être responsables à elles-seules, d'une part conséquente des effets secondaires rencontrés lors du traitement par radiothérapie. Ainsi, par exemple, lors d'un traitement d'un cancer des voies aéro-digestives supérieures (ou cancer ORL), il est particulièrement important de ne pas faire subir d'irradiation aux glandes salivaires. En effet, il a été démontré que l'irradiation de ces glandes génère des effets secondaires désastreux et contre-productifs comme par exemple la survenance d'un cancer des glandes salivaire elles-mêmes. Il est donc nécessaire de préserver ces glandes salivaires. Il en est de même pour le cancer de la prostate. Il a été démontré qu'une sous-région spécifique du rectum, lorsqu'elle est surirradiée, est corrélée à de possibles saignements rectaux qui surviennent postérieurement au traitement. Il est donc important de ne pas faire subir d'irradiation excessive à cette zone qui peut être plus sensible lors du traitement radiothérapique. D'autres zones et/ou d'autres organes sont également concernés. Ces contraintes étant posées, il est nécessaire de connaître précisément la localisation de ces zones à éviter lors du traitement. Or ces zones ne sont la plupart du temps pas connues et dépendent des caractéristiques individuelles, de la forme du rectum, de son volume, de la forme de la vessie, du poids du patient, etc. Ainsi, à l'heure actuelle, il n'est pas possible de localiser, pour un patient spécifique, avec une précision suffisante, les sous-régions des organes qui doivent être évitées lors des traitements. En d'autres termes, des méthodes d'utilisation de l'information spatiale dans des modèles de prédiction, et leur inclusion dans une étape de planning en radiothérapie n'ont pas encore été rapportés. Il y a une énorme variabilité anatomique et clinique qui entrave le transfert des informations d'un modèle de population vers un patient spécifique. En conclusion, il existe des matériels qui peuvent permettre d'éviter ces zones à risque et offrir une meilleure qualité de soins aux patients, mais il n'existe pas de méthodes permettant exploiter efficacement ces matériels. 3. Résumé de l'invention La technique proposée ne présente pas ces inconvénients de l'art antérieur. Plus particulièrement, la technique proposée se rapporte à un procédé de production de données représentatives de contraintes de traitement de radiothérapie associées à un patient à traiter. Un tel procédé comprend : une phase de mise en correspondance d'au moins un atlas courant dudit ensemble d'atlas en fonction d'au moins une donnée représentative d'un patient à traiter, délivrant la carte de probabilité de région à risque ; une phase de génération d'au moins une donnée représentative d'une contrainte à appliquer à un traitement destiné audit patient en fonction de ladite carte de probabilité de région à risque. Ainsi, la technique proposée permet d'identifier des régions à risque au sein des organes de patient en fonction de données préexistantes, qui sont traitées et normalisées pour correspondre à un patient donné. Ainsi, la méthode proposée permet d'identifier des zones sensibles à l'intérieur des organes à risque, associées à un effet secondaire. La méthode permet ainsi la définition de contraintes spatiales dans des sous-régions afin de diminuer le risque de toxicité dans un contexte de planning personnalisé Selon un mode de réalisation particulier, ledit procédé comprend en outre une phase de planification d'au moins un traitement de radiothérapie en fonction de ladite au moins une donnée représentative d'une contrainte (C). Selon un mode de réalisation particulier, ledit procédé comprend en outre une phase de prédiction d'une toxicité pour ledit patient en fonction de ladite au moins une donnée représentative dudit patient à traiter de ladite au moins une donnée représentative d'une contrainte (C) générée. Selon un mode de réalisation particulier, ladite phase de mise en correspondance d'au moins un atlas courant dudit ensemble d'atlas comprend : une étape de recalage, à partir d'au moins un paramètre de recalage déterminé, d'au moins une partie desdites données dudit atlas courant sur ladite au moins une donnée représentative dudit patient, délivrant un atlas courant recalé contenant des données probabilistes brutes ; une étape de calcul, à l'aide de ladite au moins une donnée représentative dudit patient, de ladite carte de probabilité de région à risque en fonction desdites données probabilistes brutes. Selon un mode de réalisation particulier, ladite phase de mise en correspondance d'au moins un atlas courant dudit ensemble d'atlas comprend en outre : une étape de calcul d'une similarité entre ledit atlas courant et ladite au moins une donnée représentative dudit patient, délivrant une valeur de similarité ; et une étape de pondération de données probabilistes contenues dans ledit atlas courant recalé en fonction de ladite valeur de similarité ; et en ce que ladite une étape de calcul, à l'aide de ladite au moins une donnée représentative dudit patient, de ladite carte de probabilité de région à risque est réalisée en fonction des données probabilistes pondérées. Selon un mode de réalisation particulier, ladite phase de génération d'au moins une donnée représentative d'une contrainte comprend : une étape de seuillage de la carte de probabilité de région à risque en fonction d'un paramètre de seuillage prédéterminé ; une étape d'obtention en fonction de ladite carte de probabilité de région à risque seuillée, d'au moins une position tridimensionnelle de ladite zone à risque pour ledit patient. une étape de génération ladite au moins une donnée représentative d'une contrainte de traitement en fonction de ladite position. Dans un autre mode de réalisation, la technique se rapporte également à un dispositif de production de données représentatives de contraintes de traitement de radiothérapie associées à un patient à traiter. Un tel dispositif comprend : des moyens de mise en correspondance d'au moins un atlas courant dudit ensemble d'atlas en fonction d'au moins une donnée représentative d'un patient à traiter, délivrant la carte de probabilité de région à risque ; des moyens de génération d'au moins une donnée représentative d'une contrainte à appliquer à un traitement destiné audit patient en fonction de ladite carte de probabilité de région à risque. La technique décrite se rapporte également à un système de production de données représentatives de contraintes de traitement de radiothérapie associées à un patient à traiter. Un tel système comprend : des moyens de mise en correspondance d'au moins un atlas courant dudit ensemble d'atlas en fonction d'au moins une donnée représentative d'un patient à traiter, délivrant la carte de probabilité de région à risque ; des moyens de génération d'au moins une donnée représentative d'une contrainte à appliquer à un traitement destiné audit patient en fonction de ladite carte de probabilité de région à risque. Selon une implémentation préférée, les différentes étapes des procédés selon l'invention sont mises en oeuvre par un ou plusieurs logiciels ou programmes d'ordinateur, comprenant des instructions logicielles destinées à être exécutées par un processeur de données d'un module relais selon l'invention et étant conçu pour commander l'exécution des différentes étapes des procédés. En conséquence, l'invention vise aussi un programme, susceptible d'être exécuté par un ordinateur ou par un processeur de données, ce programme comportant des instructions pour commander l'exécution des étapes d'un procédé tel que mentionné ci-dessus. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée, ou dans n'importe quelle autre forme souhaitable. L'invention vise aussi un support d'informations lisible par un processeur de données, et comportant des instructions d'un programme tel que mentionné ci-dessus. Le support d'informations peut être n'importe quelle entité ou dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen d'enregistrement magnétique, par exemple une disquette (floppy disc) ou un disque dur. D'autre part, le support d'informations peut être un support transmissible tel qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet. Alternativement, le support d'informations peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé en question. Selon un mode de réalisation, l'invention est mise en oeuvre au moyen de composants logiciels et/ou matériels. Dans cette optique, le terme "module" peut correspondre dans ce document aussi bien à un composant logiciel, qu'à un composant matériel ou à un ensemble de composants matériels et logiciels. Un composant logiciel correspond à un ou plusieurs programmes d'ordinateur, un ou plusieurs sous-programmes d'un programme, ou de manière plus générale à tout élément d'un programme ou d'un logiciel apte à mettre en oeuvre une fonction ou un ensemble de fonctions, selon ce qui est décrit ci-dessous pour le module concerné. Un tel composant logiciel est exécuté par un processeur de données d'une entité physique (terminal, serveur, passerelle, routeur, etc.) et est susceptible d'accéder aux ressources matérielles de cette entité physique (mémoires, supports d'enregistrement, bus de communication, cartes électroniques d'entrées/sorties, interfaces utilisateur, etc.). De la même manière, un composant matériel correspond à tout élément d'un ensemble matériel (ou hardware) apte à mettre en oeuvre une fonction ou un ensemble de fonctions, selon ce qui est décrit ci-dessous pour le module concerné. Il peut s'agir d'un composant matériel programmable ou avec processeur intégré pour l'exécution de logiciel, par exemple un circuit intégré, une carte à puce, une carte à mémoire, une carte électronique pour l'exécution d'un micrologiciel (firmware), etc. Chaque composante du système précédemment décrit met bien entendu en oeuvre ses propres modules logiciels. Les différents modes de réalisation mentionnés ci-dessus sont combinables entre eux pour la mise en oeuvre de l'invention. 4. Figures D'autres caractéristiques et avantages de l'invention apparaîtront plus clairement à la lecture de la description suivante d'un mode de réalisation préférentiel, donné à titre de simple exemple illustratif et non limitatif, et des dessins annexés, parmi lesquels : la figure 1 présente un synoptique de la technique proposée ; la figure 2 illustre la première phase de la technique proposée ; la figure 3 illustre la deuxième phase de la technique proposée ; la figure 4 décrit un dispositif pour la mise en oeuvre de la technique proposée. 5. Description d'un mode de réalisation 5.1. Rappel du principe de la technique La technique proposée repose sur plusieurs éléments distincts. En premier lieu, la technique repose sur la création d'un ensemble de données appelé ensemble d'atlas. Un atlas est un ensemble de données représentatives de certaines caractéristiques. Plus particulièrement, un atlas comprend une image (une image médicale, de préférence en trois dimensions, comprenant des pixels en trois dimensions, nommés voxels) à laquelle plusieurs caractéristiques sont associées. Plus particulièrement, un atlas intègre les emplacements et les formes des structures anatomiques et les relations spatiales entre ces formes et ces structures. En règle générale, on distingue deux catégories d'atlas. La première catégorie est directement obtenue à partir de données d'un seul patient. Il s'agit d'un atlas associé à un patient. La seconde catégorie est obtenue à partir de données correspondant à plusieurs patients. Dans ce cas, il s'agit d'atlas comprenant des données combinées, ou moyennes, dont l'objectif est de servir de référence. Un problème lié à ces atlas de l'art antérieur est qu'ils correspondent à un ensemble d'organes donnés. Plus particulièrement, dans le cas de données correspondant au cancer de la prostate, les atlas se rapportent à l'ensemble de la région prostatique, laquelle comprend outre la prostate, le rectum, la vessie, l'urètre, etc. Ces atlas ne se rapportent donc pas à un élément précis. Comme l'objet de la technique décrite est d'éviter ou de diminuer les effets secondaires et que ces effets secondaires sont attribués à une irradiation, à une trop forte dose, d'une ou plusieurs zones sensibles à l'intérieur des organes et que ces zones sensibles ne sont pas a priori connues pour un patient donné, il est nécessaire d'utiliser un autre type d'atlas. Ce nouveau type d'atlas est de type combinatoire, comprenant une image en trois dimensions issue de la combinaison de plusieurs images en trois dimensions d'origine. Ces images en trois dimensions d'origine sont soit issues d'un même patient (dans ce cas les images d'origine ont été obtenues à des moments différents), soit issues de plusieurs patients. On dispose donc, dans cet atlas, de données anatomiques combinées. Ce nouveau type d'atlas présente la particularité de comprendre, en plus des données anatomiques combinées, des données cliniques obtenues postérieurement. Parmi les données cliniques obtenues postérieurement, on a notamment la survenance d'effets secondaires. Dans le cas du cancer de la prostate, il s'agit par exemple de saignements du rectum. Dans le cas d'un cancer ORL, il s'agit par exemple de la perte du gout ou de l'absence de salive. Dans ce nouvel atlas, on dispose également d'une distribution de doses combinée. Cette distribution de dose de l'atlas résulte d'un calcul réalisé sur des distributions de doses unitaires. Ces distributions de doses unitaires sont celles qui sont associées aux images unitaires qui ont été utilisées pour produire l'image combinée comprise dans l'atlas. En fonction des données cliniques (effets secondaires), des distributions de doses, on produit un atlas à partir duquel une analyse par voxel peut être mise en oeuvre pour déterminer la dose reçue pour un ou un ensemble de voxels donné et corréler cette information à la survenance des effets secondaires. Un atlas, dans la technique proposée, porte donc une information associée à une zone probable liée aux effets secondaires mesurés. Ainsi, à partir de données de plusieurs patients dans une population, une méthode automatique permet de construire au moins un atlas qui définit des sous-régions à l'intérieur des organes qui sont en relation avec la dose et la toxicité. L'idée, à la base de la technique proposée, est de disposer de plusieurs atlas de ce type. Ces atlas sont regroupés dans un ensemble d'atlas. Les atlas de cet ensemble d'atlas sont disjoints par au moins une caractéristique (au moins une donnée clinique et/ou au moins une donnée anatomique). Compte tenu de la disponibilité des effets secondaires associés aux patients, on dispose, pour un atlas courant de l'ensemble d'atlas, en fonction de la distribution de dose combinée, d'une localisation probable d'une zone ayant contribuée à la survenance de l'effet secondaire que l'on souhaite éviter (il s'agit de données probabilistes de localisation de zones associées audit atlas courant). Ce qui vient d'être exposé constitue le principe de la technique proposée. On décrit, en relation avec la figure 1, les différentes phases de la technique proposée. En premier lieu, on dispose d'un ensemble d'atlas (E(A)), comprenant x atlas (A1 à Ax). Cet ensemble d'atlas est utilisé, dans une première phase (Ph1), pour être mis en correspondance avec un patient à planifier pour lequel on dispose d'une part d'une image (le) et d'autre part de données cliniques initiales (Da). Ces données cliniques initiales se rapportent par exemple à l'age, au poids, à la taille, aux antécédents médicaux et à la pathologie cancéreuse. Ainsi, à l'aide de l'ensemble d'atlas (E(A)) et d'au moins une image en trois dimensions (le) correspondant à un patient à planifier (et à traiter), on réalise au moins une mise en correspondance des atlas de l'ensemble des atlas et du patient à planifier. Ainsi, pour un nouveau patient à planifier, le(s) atlas sont mis en correspondance (spatiale, clinique, temporelle). De manière complémentaire, lorsque plusieurs atlas sont successivement mis en correspondance avec l'image de base, ces atlas peuvent être pondérés en tenant compte de la similarité entre chaque atlas et le patient lors de la mise en correspondance. Cette phase de mise en correspondance permet d'obtenir, pour le patient à planifier, une carte de probabilité de région à risque (Pee). Cette carte de probabilité permet de situer, sur le patient à planifier, la zone ou les zones qui seront, pour ce patient, les plus probablement à l'origine d'effets secondaires. Cette carte est un ensemble de données comprenant, pour un voxel donné, une probabilité de se trouver dans une région associée à une radio-sensibilité plus importante. La deuxième phase (Ph2) consiste à générer des contraintes spécifiques au patient, en fonction de la carte de probabilité de région à risque (PRR). Ces contraintes spécifiques au patient peuvent être générées par la console de traitement, de manière automatique. On appelle contrainte (C) une structure de données représentative d'une zone ou sous région que l'on suppose être sensible aux radiations. Cette structure de données peut avoir plusieurs formes. Une première forme est une image tridimensionnelle binaire. Cette image comprend des 0 aux voxels qui ne font pas partie de la zone à risque et des 1 aux voxels qui font partie de la zone à risque. Une deuxième forme est une image tridimensionnelle probabiliste dans laquelle les voxel portent une information probabiliste. Une autre forme est une structure de données descriptive, identifiant les coordonnées de la zone à risque dans l'espace de coordonnées de la console de planification. Les zones à l'intérieur des organes à risques étant identifiées, la troisième phase (Ph3) consiste à effectuer un planning de traitement du patient. Cette tâche est également mise en oeuvre de manière automatique par l'intermédiaire de la console de traitement. Les contraintes permettent à la console de planifier le traitement de ce patient par exemple en évitant les sous-régions des organes à risques ou en réduisant la dose planifié sur ces mêmes sous-régions. En fonction de la forme que prennent les contraintes, l'intégration de celles-ci dans la console de planification est modifiée. La quatrième phase (Ph4), qui se tient en parallèle de la troisième, consiste à mettre en oeuvre une prédiction de la toxicité pour le patient en fonction des caractéristiques individuelles et des nouvelles contraintes générées au regard d'un traitement en IMRT, ARC, Cyber-Knife, qui permettent le ciblage de petites régions de manière précise. Cette planification a pour but de diminuer le risque de toxicité lié à la zone délimité. Plus particulièrement, cette phase de prédiction comprend au moins une étape d'évaluation d'une dose reçue par ledit patient lors dudit traitement, comme par exemple au moins une étape de calcul d'une dose volumique reçue au sein de ladite zone à risque en fonction de la carte de probabilité de région à risque (PRR) et au moins une étape de calcul d'une probabilité de toxicité résultante de l'application d'une dose sur ladite zone à risque. Ce calcul est itéré autant de fois qu'il y a de séances de traitement afin d'obtenir une probabilité de survenance d'effet secondaire du traitement. Dans un contexte de contrôle de la tumeur et de diminution de la toxicité, la technique proposée fournit donc un ensemble d'étapes nécessaires à la personnalisation du traitement avec une quantification possible du bénéfice en termes de toxicité du traitement. Pour finir, la dernière phase (Ph5) consiste à réintroduire, sous une forme anonyme, dans l'ensemble d'atlas, les données cliniques ultérieures (Dcu) à l'issue de la troisième phase dans l'ensemble d'atlas (E(A)). Plus particulièrement, à partir de l'introduction de ces nouvelles données, il est possible de recalculer un nouvel ensemble d'atlas, lequel prend en considération notamment les effets secondaires du patient qui vient d'être traité. 5.2. Description d'un mode de réalisation Dans ce mode de réalisation, la technique précédemment décrite est appliquée au traitement du cancer de la prostate. Dans un contexte d'individualisation du traitement de radiothérapie par modulation d'intensité (IMRT), la technique présentée dans ce mode de réalisation génère des contraintes spatiales spécifiques à un patient donné dans le but de réduire la toxicité des doses et/ou de prédire la toxicité des doses pour des organes ou des zones à risque. Plus particulièrement, un objet de cette technique, dans une application particulière est de réduire les saignements rectaux, voire de supprimer leur apparition. On considère, dans ce mode de réalisation, que l'on dispose d'un ensemble d'atlas (E(A)). Cet ensemble d'atlas est construit par une méthode de construction d'atlas qui est décrite par la suite. Cet ensemble d'atlas est construit préalablement à l'application de la présente technique. Selon une première variante, il est construit immédiatement avant la mise en oeuvre de la présente technique. Selon une deuxième variante, l'ensemble d'atlas est construit en amont, par exemple lors d'une phase d'installation ou de paramétrage d'un dispositif de traitement (une console de traitement). Selon une troisième variante l'atlas est disponible par l'intermédiaire d'un réseau de communication et est par exemple mutualisé entre différents équipements (différentes consoles) et la technique proposée est mise en oeuvre sur ce réseau de communication. Dans ce mode de réalisation la technique comprend : une phase (Phi) de mise en correspondance d'au moins un atlas courant (A,) dudit ensemble d'atlas (E(A)) en fonction d'au moins une donnée représentative d'un patient à traiter, délivrant la carte de probabilité de région à risque (PRR) ; une phase de génération (Ph2) d'au moins une donnée représentative d'une contrainte (C) à appliquer à un traitement destiné audit patient en fonction de ladite carte de probabilité de région à risque (PRR). Plus particulièrement, en relation avec la figure 2, la phase de mise en correspondance (Ph1) comprend, pour ledit atlas courant (A,) dudit ensemble d'atlas (E(A)) : une étape de recalage (Ph10), à partir d'au moins un paramètre de recalage déterminé (PR), d'au moins une partie desdites données dudit atlas courant (A,) sur ladite au moins une donnée représentative dudit patient (le, Dc,), délivrant un atlas courant recalé (AIR) contenant des données probabilistes brutes ; Cette phase (Phl) de mise en correspondance peut être mise en oeuvre pour chaque atlas dudit ensemble d'atlas (E(A)). En fonction d'un paramètre de configuration, cette phase comprend en outre : une étape de calcul d'une similarité (Ph11) entre ledit atlas courant recalé (AIR) et ladite au moins une donnée représentative dudit patient (le, Dc,), délivrant une valeur de similarité (SimAIR) ; et une étape de pondération (Ph12) de données probabilistes contenues dans ledit atlas courant recalé (AIR) en fonction de ladite valeur de similarité délivrant une valeur de pondération (PAIR) ; Par ailleurs, la phase de mise en correspondance comprend également : une étape de calcul (Ph13), à l'aide de ladite au moins une donnée représentative dudit patient (le, Dc,), de ladite carte de probabilité de région à risque (PRR) en fonction soit desdites données probabilistes brutes, soit en fonction des données probabilistes pondérées. Lorsque l'ensemble d'atlas comprend plusieurs atlas, l'étape de calcul Ph13 est réalisée à l'aide des données des atlas qui ont fait l'objet d'un recalage. Ainsi, les étapes Ph10, Ph11 et Ph12 peuvent être répétées autant de fois qu'il y a d'atlas dans l'ensemble d'atlas (E(A)) ou sur un sous ensemble des atlas de l'ensemble d'atlas (E(A)). Ceci est représenté sur la figure par les Atlas Al2 à A1+2. Dès lors, l'étape Ph13 est soit mise en oeuvre au coup par coup en fonction des résultats des étapes précédentes, soit directement en compilant en une seule étape les recalages et les pondérations précédemment réalisées. On dispose ainsi, à l'issue de cette phase, de données probabilistes (PRO, en quelque sorte calquées sur les données du patient, et plus particulièrement sur les données anatomiques du patient, sous la forme d'une carte de probabilité de région à risque (PRR) en trois dimensions. Cette carte comprend, pour chaque voxel, en sus des données usuelles, une probabilité d'une présence d'une zone à risque. Les précédentes étapes sont mises en oeuvre de manière automatique, soit par la console de traitement elle-même, soit par l'intermédiaire d'une ressource de calcul accessible par la voie d'un réseau de communication. Plus particulièrement, la phase de génération (Ph2) de ladite au moins une donnée représentative d'une contrainte (C) comprend : une étape de seuillage (Ph20) de la carte de probabilité de région à risque (PRR) en fonction d'un paramètre de seuillage prédéterminé (PS). Dans un mode de réalisation spécifique, ledit paramètre de seuillage est égal à 0,5. On considère ainsi que toute probabilité d'une sensibilité inférieure à 0,5 équivaut à une absence de sensibilité et qu'une probabilité d'une sensibilité supérieure ou égale à 0,5 équivaut à une sensibilité. Dans un mode de réalisation, cette étape de seuillage est mise en oeuvre sur la carte de probabilité de région à risque (PRR) indépendamment voxel par voxel. Dans un autre mode de réalisation, cette étape de seuillage est mise en oeuvre de manière ensembliste, en considérant par exemple qu'un voxel de faible probabilité (inférieure à 0,5) qui comprend dans son voisinage un certain nombre de voxel (par exemple plus de treize) de forte probabilité (supérieure à 0,5) est un anachronisme et est considéré comme un voxel à forte probabilité. En d'autres termes, le procédé comprend une étape de seuillage des probabilités de présence d'une zone sensible en fonction d'un paramètre de seuillage prédéterminé, délivrant au moins une donnée représentative d'une présence d'une zone sensible. - une étape d'obtention (Ph21) en fonction de ladite carte de probabilité de région à risque seuillée (PRRs), d'au moins une position tridimensionnelle de ladite zone à risque pour ledit patient (P3DzR). L'obtention de cette position peut soit être réalisée par la transmission, par l'intermédiaire d'un réseau de communication, d'une position calculée au sein du réseau ou encore par l'intermédiaire d'un calcul réalisé au sein de la console de traitement elle-même. Dans au moins un mode de réalisation particulier, la zone à risque est assimilée à un organe au sein de ladite console. Ceci permet de contourner une impossibilité, pour certaines consoles, d'identifier des zones particulière de certains organes. En d'autres termes, on a une étape d'intégration, au sein d'une application de planification de traitement, des données représentatives d'une présence d'une zone sensible. une étape de génération (Ph22) ladite au moins une donnée représentative d'une contrainte (C) de traitement en fonction de ladite position (P3DzR). On dispose ainsi, à l'issue de cette phase, de contraintes associées au patient dans le cadre du traitement qui lui est destiné. Ces contraintes peuvent soit prendre la forme de doses spécifiques à appliquer à la ou les zones à risques identifiées. Ces contraintes peuvent également se traduire par un ensemble de coordonnées spatiales définies dans la console et marquées de manière spécifiques pour l'opérateur de traitement de radiothérapie afin qu'il détermine lui-même les doses qui vont être appliquées. 5.3. Construction de l'ensemble d'atlas L'ensemble d'atlas (E(A)) est un modèle multi-varié de prédiction de toxicité. C'est un résultat d'une analyse locale (par voxel) de doses individuelles sur une population donnée. En d'autres termes, à la différence des atlas usuels, l'ensemble d'atlas (E(A)) de la présente technique comprend des données cliniques corrélées aux doses planifiées ou effectivement reçues par les patients. Ces données sont intégrées à l'atlas. Ainsi, cet ensemble d'atlas est un modèle multi-varié (comprenant plusieurs grades de toxicité, plusieurs espaces de référence, plusieurs types de variables cliniques) résultat d'une analyse de population. Il permet d'établir des relations entre des distributions de doses en trois dimensions, obtenues lors de la planification ou par mesure et la toxicité de ces doses reçues. Plusieurs espaces de référence (ou atlas) sont utilisés en fonction de critères de similarité inter-individus. On connaît, dans les techniques antérieures, des méthodes permettant d'obtenir un atlas à partir d'une population donnée. À la différence de ces techniques antérieure, la technique décrite ici permet d'obtenir plusieurs atlas à partir d'une population donnée. Ces différents atlas peuvent être obtenus à partir d'une méthode de base en revanche, les différents atlas correspondent chacun à un paramètre préalablement déterminé. En effet, dans les techniques antérieures, un atlas correspond à un ensemble de données correspondant à une synthèse de paramètres. Dans la technique objet de la présente demande, un atlas correspond à une synthèse orientée à partir d'un seul paramètre de départ. Plus précisément, la technique ici présentée comprend : une phase de reconnaissance et de marquage, par un ou plusieurs méthodes experts, de la position d'un ou plusieurs sous-régions des organes dans les différentes images destinées à l'atlas ; Les experts utilisés sont des méthodes automatique (basées sur des analyses de population). une phase de détermination, parmi les images précédentes, d'une image de référence ; une phase d'application, en fonction d'une pluralité de doses préalablement connues, d'une distribution de doses sur les images marquées. Cette phase d'application peut être mise en oeuvre par l'intermédiaire d'une déformation appliquée sur les images en fonction d'une image de référence ; une phase de normalisation des images en fonction de l'image de référence ; une phase de détermination d'une distribution de dose normalisée à l'aide des informations précédentes. Ainsi, on dispose d'une distribution de dose cohérente, appliquée sur un patient particulier. 5.4. Dispositif pour la mise en oeuvre de la technique proposée On présente, en relation avec la figure 4, une architecture simplifiée d'un dispositif apte à mettre en oeuvre la technique décrite. Un tel dispositif comprend une mémoire 41, une unité de traitement 42 équipée par exemple d'un microprocesseur, et pilotée par le programme d'ordinateur 43, mettant en oeuvre au moins une partie du procédé tel que décrit. Dans au moins un mode de réalisation, l'invention est mise en oeuvre sous la forme d'une application installée sur un dispositif de planification. Un tel dispositif comprend : des moyens de mise en correspondance d'au moins un atlas courant (A,) dudit ensemble d'atlas (E(A)) en fonction d'au moins une donnée représentative d'un patient à traiter, délivrant la carte de probabilité de région à risque (PRR). Ces moyens peuvent se matérialiser sous la forme d'une interface de connexion (I) à un ou plusieurs réseaux de communication. Il peut s'agir d'interfaces logicielles ou d'interfaces matérielles (de type carte réseau ou modules matériels de communication réseau). Ces moyens peuvent se présenter sous la forme de ressources de calcul dédiées permettant de réaliser les traitements au sein même du dispositif; des moyens de génération (Ph2) d'au moins une donnée représentative d'une contrainte (C) à appliquer à un traitement destiné audit patient en fonction de ladite carte de probabilité de région à risque (PRR) des moyens de transmission (T), à un dispositif de traitement, des données précédemment générées. Dans au moins un mode de réalisation, la technique décrite peut être mis en oeuvre par l'intermédiaire d'un réseau de communication auquel un dispositif est connecté. Dans au moins un mode de réalisation, la technique décrite est mise en oeuvre au sein d'une console de planification de traitement de radiothérapie.5 REVENDICATIONS1. Procédé de production de données représentatives de contraintes de traitement de radiothérapie associées à un patient à traiter, procédé caractérisé en ce qu'il comprend : une phase (Phi) de mise en correspondance d'au moins un atlas courant (AI) d'un ensemble d'atlas (E{A)) en fonction d'au moins une donnée représentative d'un patient à traiter, délivrant la carte de probabilité de région à risque (PRR) ; 10 une phase de génération (Ph2) d'au moins une donnée représentative d'une contrainte (C) à appliquer à un traitement destiné audit patient en fonction de ladite carte de probabilité de région à risque (PRR). 2. Procédé selon la revendication 1, caractérisé en ce qu'il comprend en outre une phase de planification d'au moins un traitement de radiothérapie en fonction de 15 ladite au moins une donnée représentative d'une contrainte (C). 3. Procédé selon la revendication 1, caractérisé en ce que comprend en outre une phase de prédiction d'une toxicité pour ledit patient en fonction de ladite au moins une donnée représentative dudit patient à traiter de ladite au moins une donnée représentative d'une contrainte (C) générée. 20 4. Procédé selon la revendication 1, caractérisé en ce que ladite phase (Phi) de mise en correspondance d'au moins un atlas courant (A1) dudit ensemble d'atlas (E(A)) comprend : une étape de recalage, à partir d'au moins un paramètre de recalage déterminé, d'au moins une partie desdites données dudit atlas courant sur ladite au moins 25 une donnée représentative dudit patient, délivrant un atlas courant recalé contenant des données probabilistes brutes ; une étape de calcul, à l'aide de ladite au moins une donnée représentative dudit patient, de ladite carte de probabilité de région à risque (PRR) en fonction desdites données probabilistes brutes. 30Procédé selon la revendication 4, caractérisé en ce que ladite phase (Phi) de .mise en correspondance d'au moins un atlas courant (A1) dudit ensemble d'atlas (E(t)) comprend en outre : une étape de calcul d'une similarité entre ledit atlas courant et ladite au moins une donnée représentative dudit patient, délivrant une valeur de similarité ; et une étape de pondération de données probabilistes contenues dans ledit atlas courant recalé en fonction de ladite valeur de similarité ; et en ce que ladite une étape de calcul, à l'aide de ladite au moins une donnée représentative dudit patient, de ladite carte de probabilité de région à risque 10 (PRR) est réalisée en fonction des données probabilistes pondérées. 6. Procédé selon la revendication 1, caractérisé en ce que une ladite phase de génération (Ph2) d'au moins une donnée représentative d'une contrainte (C) comprend : 15 une étape de seuillage (Ph20) de la carte de probabilité de région à risque (PRR) en fonction d'un paramètre de seuillage prédéterminé (PS) ; une étape d'obtention en fonction de ladite carte de probabilité de région à risque seuillée, d'au moins une position tridimensionnelle de ladite zone à risque pour ledit patient. 20 une étape de génération ladite au moins une donnée représentative d'une contrainte (C) de traitement en fonction de ladite position. 7. Dispositif de production de données représentatives de contraintes de traitement de radiothérapie associées à un patient à traiter, dispositif 25 caractérisé en ce qu'il comprend : des moyens de mise en correspondance d'au moins un atlas courant (Al) d'un ensemble d'atlas (E(A)) en fonction d'au moins une donnée représentative d'un patient à traiter, délivrant la carte de probabilité de région à risque (PRR);des moyens de génération (Ph2) d'au moins une donnée représentative d'une contrainte (C) à appliquer à un traitement destiné audit patient en fonction de ladite carte de probabilité de région à risque (PRR). 8. Système de production de données représentatives de contraintes de traitement de radiothérapie associées à un patient à traiter, système caractérisé en ce qu'il comprend : des moyens de mise en correspondance d'au moins un atlas courant (A) d'un ensemble d'atlas (E(A)) en fonction d'au moins une donnée représentative d'un 10 patient à traiter, délivrant la carte de probabilité de région à risque (PRO ; des moyens de génération (Ph2) d'au moins une donnée représentative d'une contrainte (C) à appliquer à un traitement destiné audit patient en fonction de ladite carte de probabilité de région à risque (PRO; 15 9. Système de production de données selon la revendication 8, caractérisé en ce qu'il comprend des moyens de prédiction d'une toxicité pour ledit patient en fonction de ladite au moins une donnée représentative dudit patient à traiter et de ladite au moins une donnée représentative d'une contrainte (C) générée. 20 10. Produit programme d'ordinateur téléchargeable depuis un réseau de communication et/ou stocké sur un support lisible par ordinateur et/ou exécutable par un microprocesseur, caractérisé en ce qu'il comprend des instructions de code de programme pour l'exécution d'un procédé de fourniture selon la revendication 1, lorsqu'il est exécuté sur un ordinateur. 25