WO2014170490A2 - Procédé de contrôle de qualité de positionnement radiothérapique - Google Patents
Procédé de contrôle de qualité de positionnement radiothérapique Download PDFInfo
- Publication number
- WO2014170490A2 WO2014170490A2 PCT/EP2014/058024 EP2014058024W WO2014170490A2 WO 2014170490 A2 WO2014170490 A2 WO 2014170490A2 EP 2014058024 W EP2014058024 W EP 2014058024W WO 2014170490 A2 WO2014170490 A2 WO 2014170490A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image
- repositioning
- planning
- displacement
- tumor structure
- Prior art date
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
- A61N2005/1061—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using an x-ray imaging system having a separate imaging source
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1064—Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring
- A61N5/1069—Target adjustment, e.g. moving the patient support
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1075—Monitoring, verifying, controlling systems and methods for testing, calibrating, or quality assurance of the radiation treatment apparatus
Definitions
- the invention relates to the processing of medical imaging data.
- the invention relates more particularly to a medical data processing method allowing control after radiotherapy treatment.
- radiotherapy is a method of treating cancers that uses radiation to destroy cancer cells.
- the purpose of irradiation is to destroy tumor cells while sparing healthy peripheral tissue.
- a linear electron accelerator is used which produces an irradiation beam.
- Image-guided radiation therapy is an irradiation technique that uses accelerator-based imaging, also known as follow-up imaging, (most commonly CBCT “cone beam computed tomography” for “Digital tomography based on a conical beam”) to locate the tumor target within the patient while the patient is lying under the accelerator, and thus ensure the correct positioning of the tumor target before the session of irradiation. If the tumor target is shifted under the accelerator relative to its theoretical position defined in the planning scanner (the planning scanner delivering one or more planning images for its part), the table on which the patient is lying is moved.
- accelerator-based imaging also known as follow-up imaging, (most commonly CBCT “cone beam computed tomography” for “Digital tomography based on a conical beam”
- This displacement of the patient is based on a "rigid" registration operation of the tumor target between imaging data from the planning scanner and the imaging data from the CBCT.
- This registration is called “rigid” because it is based on identical transformations applied in all points of the volume (rotation, translation).
- This registration must be both accurate and performed very quickly (a few seconds) during the irradiation session.
- the use of materials necessary for irradiation is expensive and the number of patients to be treated is large (it is therefore necessary that patients be treated as quickly as possible) and secondly prolonged irradiation patients (for example by multiplying imaging data from CBCT) is not desirable.
- the registration task which is essential because it allows the patient to be properly repositioned, is therefore performed as quickly as possible, which creates a potential risk of error.
- a quality control tool is crucial because its absence poses a high risk of recurrence if the registration tumoral performed is approximate. Indeed, at best an approximate registration leads to the partial treatment of the tumor target. At worst, an approximate registration results in irradiation of an unaffected organ to the detriment of the treatment of the tumor target.
- the first solution corresponds to the use of radio-opaque markers implanted in the tumor before the treatment and visualized by two orthogonal images.
- the mapping of the markers is more accurate than the two registration modalities described above. This mapping can be used for both real-time registration but can also be used to quantify the accuracy of the readjustments performed.
- This first solution requires the insertion of radiopaque markers. This solution is therefore invasive and limited because of the fact that the three-dimensional limits of the tumor are not taken into account when matching the orthogonal images.
- the second solution is the manual matching of several anatomical reference points visualized on the one hand on the imaging data from the scanner and on the other hand on the imaging data from the CBCT.
- This approach includes some approximation because of the difficulty of identifying very specific anatomical points.
- the proposed technique does not have these disadvantages of the prior art. More particularly, the proposed technique relates to a method for determining the efficiency of a displacement previously performed on a patient during a radiotherapy session, a method implemented by a data processing device from least one treatment planning image and at least one treatment tracking image. Such a method comprises: a phase of obtaining at least one position of a tumor structure within said tracking image, taking into account on the one hand said at least one planning image and at least one specific datum associated with said image planning;
- a first calculation phase starting from said at least one position of said tumor structure within said at least one tracking image, and at least one data representative of a displacement applied to said patient, of a displacement of said tumor structure, delivering at least one datum representative of a repositioning performed;
- a second calculation phase starting from at least one position of said tumor structure within said at least one planning image and said at least one position of said tumor structure within said at least one tracking image, at least one datum representative of optimal repositioning;
- this technique makes it possible to calculate, after the radiotherapy session (s), a value of the positioning achieved, for example by controlling the operations that have been carried out by the practitioners so that they can, during future operations, improve their techniques. .
- said at least one specific datum associated with said planning image comprises at least one datum representative of a contour of said tumor structure.
- said phase of obtaining at least one position of a tumor structure comprises:
- said first calculation phase comprises:
- said second calculation phase comprises:
- the proposed technique also relates, in at least one embodiment, to a device for determining the efficiency of a displacement previously performed on a patient during a radiotherapy session, from at least one planning image of treatment and at least one image of treatment monitoring.
- a device for determining the efficiency of a displacement previously performed on a patient during a radiotherapy session from at least one planning image of treatment and at least one image of treatment monitoring.
- Such a device comprises:
- calculating means from said at least one position of said tumor structure within said at least one tracking image, and at least one data representative of a displacement applied to said patient, of a displacement of said tumor structure, delivering at least one datum representative of a repositioning performed;
- calculating means from at least one position of said tumor structure within said at least one planning image and said at least one position of said tumor structure within said at least one tracking image, delivering at least one representative datum of optimal repositioning;
- the proposed technique also relates to a method for calculating from at least one position of a tumor structure within at least one treatment planning image and at least one position of a tumor structure within at least one image of treatment follow-up, of at least one data representative of optimal repositioning.
- Such a method comprises:
- At least one step of calculating at least one optimal recovery index is possible to calculate such a recovery index, for example using a suitable device, during clinical routines, to improve the effectiveness of the treatments performed.
- the proposed technique also relates to a device for implementing the previously proposed calculation method.
- the various steps of the methods according to the invention are implemented by one or more software or computer programs, comprising software instructions intended to be executed by a data processor of a relay module according to the invention. invention and being designed to control the execution of the various process steps.
- the invention is also directed to a program that can be executed by a computer or a data processor, which program includes instructions for controlling the execution of the steps of a method as mentioned above.
- This program can use any programming language, and be in the form of source code, object code, or intermediate code between source code and object code, such as in a partially compiled form, or in any other desirable form.
- the invention also provides a data carrier readable by a data processor, and including instructions of a program as mentioned above.
- the information carrier may be any entity or device capable of storing the program.
- the medium may comprise storage means, such as a ROM, for example a CD ROM or a microelectronic circuit ROM, or a magnetic recording medium, for example a floppy disk or a disk. hard.
- the information medium may be a transmissible medium such as an electrical or optical signal, which may be conveyed via an electrical or optical cable, by radio or by other means.
- the program according to the invention can be downloaded in particular on an Internet type network.
- the information carrier may be an integrated circuit in which the program is incorporated, the circuit being adapted to execute or to be used in the execution of the method in question.
- the invention is implemented by means of software and / or hardware components.
- module may correspond in this document as well to a software component, a hardware component or a set of hardware and software components.
- a software component corresponds to one or more computer programs, one or more subroutines of a program, or more generally to any element of a program or software capable of implementing a function or a program. set of functions, as described below for the module concerned.
- Such a software component is executed by a data processor of a physical entity (terminal, server, gateway, router, etc.) and is capable of accessing the hardware resources of this physical entity (memories, recording media, bus communication cards, input / output electronic cards, user interfaces, etc.).
- a material component corresponds to any element of a hardware set (or hardware) able to implement a function or a set of functions, as described below for the module concerned. It may be a hardware component that is programmable or has an integrated processor for executing software, for example an integrated circuit, a smart card, a memory card, an electronic card for executing a firmware ( firmware), etc.
- FIG. 3 illustrates the second phase of the proposed technique
- Figure 4 illustrates the third phase of the proposed technique
- Figure 5 describes a device for implementing the proposed technique
- Figure 6 generally describes the method of the proposed technique. 5. Description of an embodiment
- the proposed solution relies on the use of data obtained from tracking images to obtain optimal repositioning parameters.
- This optimal repositioning serves as a reference to evaluate the repositioning performed in the clinic.
- the proposed method makes it possible to compare an ideal repositioning with a repositioning actually performed.
- One of the problems is that it is necessary to simulate the repositioning performed after the processing itself (more precisely, it is necessary to recalculate the repositioning actually performed a posteriori, to calculate an ideal repositioning, to compare these two repositionings)
- the terms "displacement" and "repositioning" when referring to the patient, are synonymous.
- the advantage of the proposed approach is that it does not require any implantation or manual selection of markers, which makes it effective in a quality control context where measurements must be performed regularly.
- At least one planning image generally CT (computed tomography) or scanner (CT scan image), acquired prior to the treatment to define an irradiation ballistics that must be used during the various treatment sessions (for example five sessions of treatment per week for eight weeks, ie forty sessions in total); This image is called the planning image.
- CT computed tomography
- scanner CT scan image
- At least one tracking image generally cone-beam CT (or Cone Beam Computed Tomography (CBCT) image), acquired during the treatment, for example at the beginning of the session, so as to reposition the tumor before a treatment session. 'irradiation.
- This image is called a tracking image.
- Displacement data is also available, whether actual displacement data or data provided by software integrated into the treatment system.
- Quantification can be used to answer at least two problems: on the one hand to evaluate the quality of the trips made and on the other hand to evaluate the risks of toxicity of the treatment according to the displacements carried out or proposed.
- the general principle of the proposed method is to use specific data available on the tumor target.
- the specific data that can be used are, for example, tumor target contours, which can be either directly present on one or the other of the images or calculated automatically from a suitable device or module. .
- the specific data may also be differences in contrast, or even color, which also make it possible to identify a zone, an organ or a tumor target.
- the specific data can also be embodied as complementary data associated with the images, for example in the form of descriptive metadata of these images.
- the important thing is to have one or more reference images, these reference images can be used to perform one or more transformations on current images to evaluate the image. on the one hand, the location of a tumor target and, on the other hand, the efficiency of a displacement performed on a treatment device during one or more treatment sessions.
- the proposed technique presented in connection with FIG. 6, comprises the following phases:
- PH 2 a first calculation phase, starting from said at least one position (PSTI) of said tumor structure within said at least one image of su ivi, and at least one datum representative of a displacement ( D PT ) applied to said patient, a displacement of said tumor structure, delivering at least one data representative of a repositioning performed (DR RR );
- a second calculation phase starting from at least one position (P S T2) of said tumor structure (S T i) within said at least one planning image (1 P ) and said at least one a position (PSTI) of said tumor structure within said at least one tracking image (1 s ), of at least one representative representative of an optimal repositioning (DR R0 ); a comparison phase (PH4) of said at least one data representative of a repositioning performed (DR RR ) and said at least one data representative of an optimal repositioning (DR R0 ), delivering at least one efficiency indicator of displacement (l ED ).
- Obtaining the position (PSTI) of the tumor structure (S T i) within the tracking image (l s ) taking into account the planning image (l P ) and a specific piece of data (D s ) associated with the planning image (l P ) can be obtained by a follow-up image registration operation on the planning image or by manual delineation or by a combination of these two methods or other suitable methods.
- PSTI position of the tumor structure
- D s a specific piece of data associated with the planning image
- the solution proposed in this embodiment relies on the use of contou rs obtained on tracking images, for example by elastic propagation and correction by an expert, for example a human expert (a software expert can also be considered in some cases ).
- This obtaining of contours (delineation) is used to obtain the parameters of an optimal repositioning (resulting from an automatic registration step).
- This optimal repositioning serves as a benchmark for assessing repositioning performed in the clinic.
- One of the problems is that it is necessary to simulate the repositioning performed in the clinic after the treatment itself.
- repositioning of the patient under the accelerator during a fraction of treatment.
- This repositioning (of the patient) aims to position the die in the same position as that considered during the planning.
- registration of the tumor between its position during the treatment session and its position during the planning.
- This registration of the tumor can itself be based on a step of "image registration” which is a method of image processing allowing, from two images, automatically estimate the geometric transformation between what is represented in two images.
- a registration of images between the planning image and the image "per-treatment” thus allows a registration of the tumor which results in a repositioning of the patient.
- the general principle of this embodiment is to use the contours of the tumor target generated, on the one hand, on the planning image and, on the other hand, on the tracking image, in order to make a comparison repositioning actually performed under the processing device (or proposed by the accelerator registration algorithm) with an ideal theoretical repositioning based on contour information and calculated later.
- contours of the planning image are routinely generated (in a standard way) routinely generated during the treatment planning phase. This delineation is done manually, using a spread from an atlas or another automatic or semi-automatic method.
- the contours of the tracking image can be generated either manually or, in order to save time, thanks to an elastic propagation of the contours of the planning image, this propagation being followed by a validation and, if necessary, manual correction by a human operator.
- Figure 2 describes such an elastic propagation of the contours.
- An image registration method estimates the non-rigid transformation between the two images.
- the implemented method relies on a set of control points distributed along a regular grid in the planning image space.
- At each of the control points is associated a geometric transformation, corresponding to the displacement, in the three-dimensional space, of this control point.
- This grid makes it possible, by interpolation, to define a geometrical transformation on the whole of the planning image.
- the transformation associated with each control point is iteratively modified to optimize a metric, measuring the match quality between the transformed planning image according to the control points and the tracking image.
- the metric used in this example is mutual information.
- contour data planning image and tracking image
- table trip and isocenter position It corresponds to the center of rotation of the treatment unit
- the displacement of the table can correspond to that actually achieved under the treatment apparatus, or to that proposed by the system integrated into the accelerator (integrating or not rotations). This depends on the purpose of the quality control performed. When one wishes to carry out a quality control of the automatic procedures, then one uses the displacements proposed by the system. When one wishes to carry out a quality control of the actions actually carried out, one uses the actual displacements.
- the isocenter of treatment is considered as analogous point in the two images. It allows to position the two images in the same repository. It is materialized on the two images considered, generally by a specific structure for the planning image (located for example in the center of the target) and by the center of the image for the tracking image.
- the geometric transformation (translation and possibly rotation) corresponding to the repositioning is applied to the contours of the tracking image so as to obtain the position of the structures during the irradiation session.
- the calculation of the recovery indices is based on a representation of each structure by a binary 3D image (that is to say containing two values, for example 0 and 1: 0 is assigned to the voxels outside the structure and 1 is attributed to the voxels inside the structure). Different recovery indices are then calculated, with for example (here for two volumes X and Y, and N (X) the number of voxels of the set X):
- the "optimal" repositioning is also calculated. It is based on an automatic method of image registration. This is a rigid registration. The objective of this ideal repositioning is to determine what the operator or the software could have done ideally in terms of displacement (repositioning), given the tools at its disposal. We are not talking here about a re-planning during treatment, but about an estimate of what it would have been possible to do ideally.
- a rigid transformation (integrating or not the rotations) is estimated between the binary volumes representing the tumor targets.
- the registration is initialized by superimposing the centers of the structures.
- Iterative gradient-descent optimization or exhaustive search is used to minimize / maximize / optimize one or more metrics, for example the sum of squared differences between: (i) the target bit volume in the planning image after the estimated geometric transformation and (ii) the bit volume corresponding to the target in the tracking image.
- This metric can also correspond directly to a recovery measure, such as the Dice or Jaccard indices described above. It can also integrate not only geometric criteria, but also dosimetric criteria, related to the dose received by the organs. For example, it can incorporate the minimal dose received by the target.
- the planned dose matrix is positioned on the tracking image following the position of the isocenter of treatment. This makes it possible to better estimate the optimal repositioning that should have been done.
- Optimal registration may also consist in finding the best compromise (optimizing) on a set of metrics (for example: dice index on the tumor and minimal dose on a nearby risk organ).
- the recovery indices calculated previously can also be calculated for this optimal repositioning. These "optimal" recovery indices are recorded for a later phase of the control. These optimal recovery indices summarize in a way the optimal result that could have been achieved. Insofar as the data on which the optimal repositioning is calculated are of the same nature as the data on which the real repositioning is calculated, the calculated indices can therefore be the same (Dice index and Jaccard index).
- the interest volumes derived from the planning and tracking images can also be compared. They may be different because of the contouring uncertainties, but also of a tumor melting during treatment. A strong divergence corresponding to this second cause may also trigger an adaptive radiotherapy procedure (replanifications).
- the whole of this approach can be integrated in a global approach associating also an evaluation of the cumulated dose in the various structures by elastic registration.
- the method described is preferably implemented by means of a dedicated computing device described in relation with FIG. 5.
- a dedicated computing device described in relation with FIG. 5.
- Such equipment comprises a memory 41, a processing unit 42 equipped for example with a microprocessor, and driven by the computer program 43, implementing at least a part of the method as described.
- the invention is implemented in the form of an application installed on a planning device.
- a device comprises:
- simulation means from said at least one position of said tumor structure within said at least one tracking image, and at least one data representative of a displacement applied to said patient, of a displacement of said tumor structure, delivering at least one datum representative of a repositioning performed;
- calculating means from at least one position of said tumor structure within said at least one planning image and at least one position of said tumor structure within said at least one tracking image, delivering at least one datum representative of optimal repositioning;
- means for comparing said at least one data representative of a repositioning performed and said at least one data representative of an optimal repositioning, delivering at least one displacement efficiency indicator are implemented through dedicated software and / or hardware modules. These modules are able to implement the reception and data transmission interfaces present on said device (I and T).
- the technique described can be implemented via a communication network to which a device is connected. In at least one embodiment, the technique described is implemented within a radiotherapy treatment planning console.
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
L'invention se rapporte à un Procédé de détermination d'une efficience d'un déplacement précédemment effectué sur un patient lors d'une séance de radiothérapie, procédé mis en œuvre par un dispositif de traitement de données à partir d'au moins une image de planification de traitement et d'au moins une image de suivi de traitement. Selon l'invention un tel procédé comprend : - une phase d'obtention d'au moins une position d'une structure tumorale au sein de ladite image de suivi; - une première phase de calcul d'un déplacement de ladite structure tumorale, délivrant au moins une donnée représentative d'un repositionnement réalisé; - une deuxième phase de calcul d'une donnée représentative d'un repositionnement optimal; - une phase de comparaison de ladite au moins une donnée représentative d'un repositionnement réalisé et de ladite au moins une donnée représentative d'un repositionnement optimal, délivrant au moins un indicateur d'efficience de déplacement.
Description
Procédé de contrôle de qualité de positionnement radiothérapique
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 permettant un contrôle postérieur au traitement radiothérapique.
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'électrons 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 œuvre 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 radiothérapeutes 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 innovations au niveau des appareils de traitement comme l'apparition des collimateurs multi-lames ont permis le développement de la radiothérapie conformationnelle avec modulation d'intensité (RCM I) qui permet d'adapter plus précisément l'irradiation à la forme des organes à traiter. De plus, 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. Il a été ainsi possible de visualiser et de traiter plus efficacement les tumeurs, ce qui se traduit par un meilleur pronostic pour les patients, une meilleure préservation des organes sains et moins d'effets secondaires.
La radiothérapie guidée par l'image est une technique d'irradiation utilisant une imagerie embarquée avec l'accélérateur, également appelée imagerie de suivi par la suite, (le plus souvent le CBCT « cone beam computed tomography » pour
« tomographie numérisée à base d'un faisceau conique ») pour localiser la cible tumorale au sein du patient alors que celui-ci est allongé sous l'accélérateur, et ainsi s'assurer du bon positionnement de la cible tumorale avant la séance d'irradiation. Si la cible tumorale est décalée sous l'accélérateur par rapport à sa position théorique définie dans le scanner de planification (le scanner de planification délivrant quant à lui une ou plusieurs images de planification), la table sur laquelle est allongé le patient est déplacée.
Ce déplacement du patient est basé sur une opération de recalage « rigide » de la cible tumorale entre des données d'imagerie issues du scanner de planification et les données d'imagerie issues du CBCT. Ce recalage est appelé « rigide » car il est basé sur des transformations identiques appliquées en tous points du volume (rotation, translation). Ce recalage doit être à la fois précis et réalisé très rapidement (quelques secondes) lors de la séance d'irradiation. En effet, d'une part l'utilisation des matériels nécessaires à l'irradiation coûte cher et le nombre de patients à traiter est grand (il faut donc que les patients soient traités le plus rapidement possible) et d'autre part une irradiation prolongée des patients (par exemple en multipliant les données d'imagerie issues du CBCT) n'est pas souhaitable. La tâche de recalage, qui est essentielle car elle permet le bon repositionnement du patient, est donc effectuée le plus rapidement possible, ce qui engendre un risque d'erreur potentielle. Il y a principalement deux modalités de recalage :
soit par une visualisation de la structure tumorale dans différents plans de l'espace et une mise en correspondance manuelle de la structure tumorale visualisée avec celle obtenue par l'intermédiaire de l'image de planification ; soit à l'aide d'algorithmes de recalage d'images intégrés au système de pilotage des accélérateurs, le recalage automatiquement réalisé étant ensuite validé visuellement par un opérateur après une éventuelle correction manuelle.
Le problème tient à des défauts des outils de contrôle de qualité de la précision du recalage tumoral effectué sous l'accélérateur et correspondant au traitement effectué (quelles que soient par ailleurs les modalités de recalage). Un outil de contrôle de qualité est crucial car son absence expose à un risque élevé de récidive si le recalage
tumoral effectué est approximatif. En effet, au mieux un recalage approximatif entraine le traitement partiel de la cible tumorale. Au pire, un recalage approximatif entraine une irradiation d'un organe non atteint au détriment du traitement de la cible tumorale.
Les solutions préexistantes de contrôle qualité du recalage sont basées sur deux approches :
La première solution correspond à l'utilisation de marqueurs radio-opaques implantés dans la tumeur avant le traitement et visualisés par deux images orthogonales. La mise en correspondance des marqueurs est plus précise que les deux modalités de recalage décrites ci-dessus. Cette mise en correspondance peut servir à la fois au recalage en temps réel mais peut également être utilisée pour quantifier la précision des recalages effectués.
Cette première solution suppose cependant l'insertion de marqueurs radio- opaques. Cette solution est donc invasive et limitée du fait de l'absence de prise en compte des limites tridimensionnelles de la tumeur lors de la mise en correspondance des images orthogonales.
La deuxième solution est la mise en correspondance manuelle de plusieurs points de références anatomiques visualisés d'une part sur les données d'imagerie issues du scanner et d'autre part sur les données d'imagerie issues du CBCT. Cette approche comprend une certaine approximation du fait de la difficulté d'identifier des points anatomiques très spécifiques.
Il est donc nécessaire de proposer une méthode de contrôle de qualité des recalages effectués préalablement aux séances de radiothérapie qui ne comprennent pas ces inconvénients de l'art antérieur.
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 détermination d'une efficience d'un déplacement précédemment effectué sur un patient lors d'une séance de radiothérapie, procédé mis en œuvre par un dispositif de traitement de données à partir d'au moins une image de planification de traitement et d'au moins une image de suivi de traitement. Un tel procédé comprend :
une phase d'obtention d'au moins une position d'une structure tumorale au sein de ladite image de suivi en tenant compte d'une part de ladite au moins une image de planification et d'au moins une donnée spécifique associée à ladite image de planification ;
une première phase de calcul, à partir de ladite au moins une position de ladite structure tumorale au sein de ladite au moins une image de suivi, et d'au moins une donnée représentative d'un déplacement appliqué audit patient, d'un déplacement de ladite structure tumorale, délivrant au moins une donnée représentative d'un repositionnement réalisé ;
une deuxième phase de calcul, à partir d'au moins une position de ladite structure tumorale au sein de ladite au moins une image de planification et de ladite au moins une position de ladite structure tumorale au sein de ladite au moins une image de suivi, d'au moins une donnée représentative d'un repositionnement optimal ;
une phase de comparaison de ladite au moins une donnée représentative d'un repositionnement réalisé et de ladite au moins une donnée représentative d'un repositionnement optimal, délivrant au moins un indicateur d'efficience de déplacement.
Ainsi, cette technique permet de calculer, après la ou les séances de radiothérapie, une valeur du positionnement réalisé, à titre par exemple de contrôle des opérations qui ont été effectuées par les praticiens afin qu'ils puissent lors de futurs opérations, améliorer leurs techniques.
Selon un mode de réalisation particulier, ladite au moins une donnée spécifique associée à ladite image de planification comprend au moins une donnée représentative d'un contour de ladite structure tumorale.
Selon un mode de réalisation particulier, ladite phase d'obtention d'au moins une position d'une structure tumorale comprend :
une étape de calcul d'au moins une transformation géométrique à appliquer à ladite image de planification en fonction de ladite image de suivi ;
une étape d'application de ladite transformation géométrique à ladite au moins une donnée spécifique.
Selon un mode de réalisation particulier, ladite première phase de calcul comprend :
une étape de calcul d'au moins une transformation géométrique à appliquer à ladite image de suivi en fonction d'une donnée représentative d'un déplacement appliqué audit patient ;
au moins une étape de calcul d'au moins un indice de recouvrement réalisé. Selon un mode de réalisation particulier, ladite deuxième phase de calcul comprend :
une étape de recalage appliqué sur ladite image de suivi à partir de ladite image de planification ;
au moins une étape de calcul d'au moins un indice de recouvrement optimal. La technique proposée concerne également, dans au moins un mode de réalisation, un dispositif de détermination d'une efficience d'un déplacement précédemment effectué sur un patient lors d'une séance de radiothérapie, à partir d'au moins une image de planification de traitement et d'au moins une image de suivi de traitement. Un tel dispositif comprend :
des moyens d'obtention d'au moins une position d'une structure tumorale au sein de ladite image de suivi en tenant compte d'une part de ladite au moins une image de planification et d'au moins une donnée spécifique associée à ladite image de planification ;
des moyens de calcul, à partir de ladite au moins une position de ladite structure tumorale au sein de ladite au moins une image de suivi, et d'au moins une donnée représentative d'un déplacement appliqué audit patient, d'un déplacement de ladite structure tumorale, délivrant au moins une donnée représentative d'un repositionnement réalisé ;
des moyens de calcul, à partir d'au moins une position de ladite structure tumorale au sein de ladite au moins une image de planification et de ladite au moins une position de ladite structure tumorale au sein de ladite au moins une
image de suivi, délivrant au moins une donnée représentative d'un repositionnement optimal ;
des moyens de comparaison de ladite au moins une donnée représentative d'un repositionnement réalisé et de ladite au moins une donnée représentative d'un repositionnement optimal, délivrant au moins un indicateur d'efficience de déplacement.
Selon un autre aspect, la technique proposée se rapporte également à un procédé de calcul à partir d'au moins une position d'une structure tumorale au sein d'au moins une image de planification de traitement et d'au moins une position d'une structure tumorale au sein d'au moins une image de suivi de traitement, d'au moins une donnée représentative d'un repositionnement optimal . Un tel procédé comprend :
une étape de recalage appliqué sur ladite image de suivi à partir de ladite image de planification ;
au moins une étape de calcul d'au moins un indice de recouvrement optimal. Ainsi, il est possible de calculer un tel indice de recouvrement, par exemple à l'aide d'un dispositif idoine, lors de routines clinique, afin d'améliorer l'efficacité des traitements réalisés.
La technique proposée se rapporte également à un dispositif pour la mise en œuvre du procédé de calcul préalablement proposé.
Selon une implémentation préférée, les différentes étapes des procédés selon l'invention sont mises en œuvre 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 dise) 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 œuvre 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 œuvre 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 œuvre 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 œuvre ses propres modules logiciels.
Les différents modes de réalisation mentionnés ci-dessus sont combinables entre eux pour la mise en œuvre 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 illustre la troisième phase de la technique proposée ;
la figure 5 décrit un dispositif pour la mise en œuvre de la technique proposée ; la figure 6 décrit de manière générale le procédé de la technique proposée. 5. Description d'un mode de réalisation
5.1. Rappel du principe de la technique
La solution proposée repose sur l'utilisation des données obtenues à partir d'images de suivi pour obtenir des paramètres d'un repositionnement optimal Ce repositionnement optimal sert de référence pour évaluer le repositionnement réalisé en clinique. En d'autres termes, la méthode proposée permet de comparer un repositionnement idéal à un repositionnement effectivement réalisé. L'une des problématiques est qu'il est nécessaire de simuler le repositionnement réalisé postérieurement au traitement lui-même (plus précisément, il est nécessaire de recalculer le repositionnement effectivement réalisé a posteriori ; de calculer un repositionnement idéal ; de comparer ces deux repositionnements) Dans la présente
description, les termes « déplacement » et « repositionnement », lorsqu'ils se rapportent au patient, sont synonymes.
L'avantage de l'approche proposée est qu'elle ne nécessite aucune implantation ni sélection manuelle de marqueurs, ce qui la rend efficace dans un contexte de contrôle qualité où des mesures doivent être effectuées régulièrement.
De manière générale, on utilise, pour mettre en œuvre la technique proposée, au moins deux images tridimensionnelles :
(i) au moins une image de planification, généralement TDM (tomodensitométrique ou image scanner ou image CT (Computed Tomography)) de planification, acquise préalablement au traitement pour définir une balistique d'irradiation qui doit être utilisée lors des différentes séances de traitement (par exemple cinq séances de traitement par semaine pendant huit semaines, soit quarante séances au total) ; Cette image est appelée image de planification.
(ii) au moins une image de suivi, généralement TDM à faisceau conique (ou image CBCT (Cone Beam Computed Tomography)), acquise au cours du traitement, par exemple en début de séance, de façon à repositionner la tumeur avant une séance d'irradiation. Cette image est appelée image de suivi.
Pour poser le problème du contrôle de la qualité des déplacements qui sont mis en œuvre pour le patient lors de la séance de traitement, on dispose ainsi de données d'imagerie qui sont issues de deux phases différentes liées au traitement. La première phase est une phase préparatoire, la deuxième est une phase exécutoire. On dispose également de données de déplacement, que ce soit des données de déplacement réelles ou des données fournies par un logiciel intégré au système de traitement.
On utilise ces données pour quantifier l'efficacité (-ou l'efficience) du déplacement réellement effectué ou proposé. La quantification peut être utilisée pour répondre à au moins deux problématiques : d'une part évaluer la qualité des déplacements effectués et d'autre part évaluer les risques de toxicité du traitement en fonction des déplacements réalisés ou proposés.
Le principe général de la méthode proposée est d'utiliser des données spécifiques disponibles sur la cible tumorale,.
Les données spécifiques qui peuvent être utilisées sont par exemple des contours de cible tumorale, q ui peuvent être soit directement présents sur l'une ou l'autre des images ou bien calculés automatiquement à partir d'un d ispositif ou d'un module adéquat. Les données spécifiques peuvent également être des différences de contraste, voire de couleu r, lesquelles permettent également d'identifier une zone, un organe ou une cible tumorale. Les données spécifiq ues peuvent également être matérialisées sous la forme de données complémentaires, associées aux images, par exemple sous la forme de métadonnées descriptives de ces images. Dans la cadre de la mise en œuvre de la technique décrite, l'important est de disposer d'une ou plusieurs images de références, ces images de référence pouvant être utilisées pour réaliser une ou plusieurs transformations sur des images cou rantes pour évaluer d'une part la localisation d'u ne cible tumorale et d'autre part l'efficience d'un déplacement effectué sur un dispositif de traitement lors d'une ou plusieurs séances de traitement.
La technique proposée, présentée en relation avec la figure 6, comprend les phases suivantes :
une phase d'obtention (PH I) d'au moins une position (PSTI) d'une structure tumorale (STi) au sein de ladite image de suivi (ls) en tenant compte d'u ne part de ladite au moins une image de planification (lP) et d'au moins une donnée spécifique (Ds) associée à ladite image de planification (lP) ;
une première phase de calcul (PH 2), à partir de ladite au moins une position (PSTI) de ladite structure tumorale au sein de ladite au moins une image de su ivi, et d'au moins une donnée représentative d'un déplacement (DPT) appliqué audit patient, d'un déplacement de ladite structure tumorale, délivrant au moins une donnée représentative d'un repositionnement réalisé (DRRR) ;
une deuxième phase de calcul (PH3), à partir d'au moins une position (PST2) de ladite structure tumorale (STi) au sein de ladite au moins une image de planification (lP) et de ladite au moins une position (PSTI) de ladite structure tumorale au sein de ladite au moins une image de suivi (ls), d'au moins une donnée représentative d'un repositionnement optimal (DRR0) ;
une phase de comparaison (PH4) de ladite au moins une donnée représentative d'un repositionnement réalisé (DRRR) et de ladite au moins une donnée représentative d'un repositionnement optimal (DRR0), délivrant au moins u n indicateur d'efficience de déplacement (lED).
L'obtention de la position (PSTI) de la structu re tumorale (STi) au sein de l'image de suivi (ls) en tenant compte de l'image de planification (lP) et d'une donnée spécifique (Ds) associée à l'image de planification (lP) peut soit être obtenue par une opération de recalage de l'image de suivi sur l'image de planification ou bien par une délinéation manuelle ou par une combinaison de ces deux méthodes ou d'autres méthodes adaptées. Dans la suite, on présente une mise en œuvre de la technique proposée basée sur des contours de cibles tumorales. Toute autre mise en œuvre, notamment basée sur d'autres données de recalage peut bien entend ue être envisagée.
5.2. Description d'un mode de réalisation
On décrit, en relation avec la figure 1, dans ce mode de réalisation, une mise en œuvre de la technique proposée. La solution proposée dans ce mode de réalisation repose sur l'utilisation de contou rs obtenus sur les images de suivi, par exemple par propagation élastique puis correction par un expert, par exemple un expert humain (un expert logiciel peut également être envisagé dans certains cas). Cette obtention de contours (délinéation) est utilisée pour obtenir les paramètres d'un repositionnement optimal (issu d'une étape de recalage automatique). Ce repositionnement optimal sert de référence pour évaluer le repositionnement réalisé en clinique. L'une des problématiques est qu'il est nécessaire de simuler le repositionnement réalisé en clinique postérieurement au traitement lui-même.
On rappelle que l'on doit faire une distinction entre le repositionnement et le recalage. On parle de "repositionnement" d u patient sous l'accélérateur lors d'une fraction de traitement. Ce repositionnement (du patient) a pour objectif de positionner la tu meur à la même position que celle considérée lors de la planification. On peut alors parler de "recalage" de la tumeur entre sa position lors de la séance de traitement et sa position lors de la planification. Ce recalage de la tumeur peut lui-même reposer sur une étape de "recalage d'images" qui est une méthode de traitement d'images permettant,
à partir de deux images, d'estimer automatiquement la transformation géométrique entre ce qui est représenté dans deux images.
Un recalage d'images entre l'image de planification et l'image « per-traitement » permet donc un recalage de la tumeur qui se traduit par un repositionnement du patient.
Le principe général de ce mode de réalisation est d'utiliser les contours de la cible tumorale générés, d'une part, sur l'image de planification et, d'autre part, sur l'image de suivi, afin de réaliser une comparaison du repositionnement effectivement réalisé sous l'appareil de traitement (ou proposé par l'algorithme de recalage de l'accélérateur) avec un repositionnement théorique idéal basé sur l'information des contours et calculé ultérieurement.
Cette comparaison est réalisée à trois niveaux :
la comparaison des contours délinéés sur les deux images (planification et suivi), qui permet d'estimer une éventuelle modification de volume ou de forme de la tumeur (cette première comparaison ne prend pas le repositionnement du patient en compte) ;
après estimation du déplacement optimal du patient, comparaison entre ce déplacement optimal et le déplacement réalisé sous machine (ou celui proposé par le logiciel ou le système de traitement) ;
- comparaison : (i) d'un indice de recouvrement entre le volume tumoral délinéé à la planification et le volume tumoral correspondant à la séance de traitement suivant le repositionnement réalisé ; (ii) du même indice correspondant au repositionnement optimal.
Ces indices de recouvrement permettant de quantifier la correspondance géométrique entre la tumeur à la planification (à partir de laquelle le traitement a été planifié) et la tumeur lors du traitement, ils sont très informatifs sur la qualité du repositionnement, et donc sur la conformité de l'irradiation de ladite tumeur par rapport à la planification.
5.2.1. Génération des contours
Les contours que l'on cherche à générer sur les images de départ se rapportent à la tumeur et éventuellement à des organes ou des zones à risque.
Les contours de l'image de planification sont, en routine clinique (c'est-à-dire de manière standard), générés systématiquement lors de la phase de planification de traitement. Cette délinéation est réalisée manuellement, à l'aide d'une propagation à partir d'un atlas ou une autre méthode automatique ou semi- automatique.
Les contours de l'image de suivi peuvent être générés soit manuellement soit, dans un objectif de gain de temps, grâce à une propagation élastique des contours de l'image de planification, cette propagation étant suivie d'une validation et, si besoin, d'une correction manuelle par un opérateur humain.
La figure 2 décrit une telle propagation élastique des contours.
Une méthode de recalage d'images permet d'estimer la transformation non- rigide entre les deux images. Par exemple, la méthode implémentée repose sur un ensemble de points de contrôle répartis suivant une grille régulière dans l'espace de l'image de planification. À chacun des points de contrôle est associée une transformation géométrique, correspondant au déplacement, dans l'espace tridimensionnel, de ce point de contrôle. Cette grille permet, par interpolation, de définir une transformation géométrique sur l'ensemble de l'image de planification. La transformation associée à chaque point de contrôle est modifiée itérativement de façon à optimiser une métrique, mesurant la qualité de correspondance entre l'image de planification transformée suivant les points de contrôle et l'image de suivi. La métrique utilisée dans cet exemple est l'information mutuelle.
Une fois le contour propagé, il est validé et éventuellement corrigé par un expert.
5.2.2. Simulation de l'impact du repositionnement réalisé sous l'appareil de traitement sur la position des organes
À partir d'une part des données de contours (image de planification et image de suivi) et d'autre part du relevé des déplacements de table et de la position de l'isocentre
(Il correspond au centre de rotation de l'appareil de traitement) de traitement pour les deux images, l'impact du repositionnement sur la position des organes est simulé.
Le déplacement de table peut correspondre à celui effectivement réalisé sous l'appareil de traitement, ou à celui proposé par le système intégré à l'accélérateur (intégrant ou non les rotations). Ceci dépend de l'objet du contrôle qualité effectué. Lorsque l'on souhaite effectuer un contrôle qualité des procédures automatiques, alors on utilise les déplacements proposés par le système. Lorsque l'on souhaite effectuer un contrôle qualité des actions effectivement réalisées, on utilise les déplacements réels.
L'isocentre de traitement est considéré comme point analogue dans les deux images. Il permet de positionner les deux images dans le même référentiel. Il est matérialisé sur les deux images considérées, généralement par une structure spécifique pour l'image de planification (située par exemple au centre de la cible) et par le centre de l'image pour l'image de suivi. La transformation géométrique (translation et éventuellement rotation) correspondant au repositionnement est appliquée aux contours de l'image de suivi de façon à obtenir la position des structures lors de la séance d'irradiation.
Suivant les contours disponibles sur l'image de suivi (cible tumorale et éventuellement organes à risque), différents indices de recouvrement des structures sont calculés, de façon à évaluer la couverture de la cible et l'irradiation des organes à risque.
Par exemple, la couverture du volume cible anatomo-clinique (tumeur, CTV (Clinical Target Volume)) par le volume cible prévisionnel (ou PTV (Planning Target Volume) correspondant au volume cible anatomo-clinique auquel des marges de sécurité sont ajoutées) est calculée (cette couverture devant être la plus élevée possible). De même la couverture des organes à risque par le volume cible prévisionnel peut être calculée (cette couverture devant être la plus faible possible).
Le calcul des indices de recouvrement repose sur une représentation de chaque structure par une image 3D binaire (c'est-à-dire contenant deux valeurs, par exemple 0 et 1 : 0 est attribué aux voxels à l'extérieur de la structure et 1 est attribué aux voxels à l'intérieur de la structure).
Différents indices de recouvrement sont alors calculés, avec par exemple (ici pour deux volumes X et Y, et N(X) le nombre de voxels de l'ensemble X):
2N(X Π Y)
■ in Ad-ice de r Dvice : ° ~ N(X) J + N(Y) J
_ N(X n Y)
■ ind-ice de J ,accar Ad : 1 = N(X u J Ces indices de recouvrement sont alors enregistrés pour une phase ultérieure du contrôle.
5.2.3. Calcul du repositionnement optimal
Le repositionnement "optimal", basé sur les contours dans ce mode de réalisation, est également calculé. Il repose sur une méthode automatique de recalage d'images. Il s'agit d'un recalage rigide.. L'objectif de ce repositionnement idéal est de déterminer ce que l'opérateur ou le logiciel aurait pu faire idéalement en termes de déplacement (de repositionnement), compte tenu des outils à sa disposition. On ne parle pas ici d'une replanification en cours de traitement, mais bien d'une estimation de ce qu'il aurait été possible de faire dans l'idéal.
Par exemple, une transformation rigide (intégrant ou pas les rotations) est estimée entre les volumes binaires représentant les cibles tumorales. Le recalage est initialisé en superposant les centres des structures. Une optimisation itérative par descente de gradient ou une recherche exhaustive est utilisée pour minimiser/maximiser/optimiser une ou plusieurs métriques, par exemple la somme des différences au carré entre : (i) le volume binaire correspondant à la cible dans l'image de planification après la transformation géométrique estimée et (ii) le volume binaire correspondant à la cible dans l'image de suivi. Cette métrique peut aussi correspondre directement à une mesure de recouvrement, comme les indices de Dice ou de Jaccard décrits précédemment. Elle peut aussi intégrer non seulement des critères géométriques, mais aussi des critères dosimétriques, liés à la dose reçue par les organes. Par exemple, elle peut intégrer la dose minimale reçue par la cible. Dans ce cas, la matrice de dose planifiée est positionnée sur l'image de suivi suivant la position de
l'isocentre de traitement. Ceci permet de mieux estimer encore le repositionnement optimal qui aurait dû être réalisé. Le recalage optimal peut aussi consister à trouver le meilleur compromis (optimiser) sur un ensemble de métriques (par exemple : indice de dice sur la tumeur et dose minimale sur un organe à risque avoisinant).
Au résultat de ce recalage correspondent des déplacements spécifiques de la cible tumorale (intégrant ou non des rotations). C'est-à-dire que pour passer de l'image de planification à l'image de suivi, on dispose de déplacements pouvant être mis en œuvre par le dispositif à partir de la transformation rigide. Il peut s'agir par exemple de deux déplacements translatifs horizontaux, un déplacement vertical et une rotation. Ces déplacements forment une transformation rigide (une transformation qu'il n'implique pas une déformation des organes, mais uniquement leurs déplacements).
Les indices de recouvrement calculés précédemment peuvent également être calculés pour ce repositionnement optimal. Ces indices de recouvrement « optimaux » sont enregistrés pour une phase ultérieure du contrôle. Ces indices de recouvrement optimaux résument en quelque sorte le résultat optimal qui aurait pu être atteint. Dans la mesure où les données sur lesquelles le repositionnement optimal est calculé sont de même nature que les données sur lesquelles le repositionnement réel est calculé, les indices calculés peuvent donc être les mêmes (indice de Dice et Indice de Jaccard).
5.2.4. Comparaison des deux modalités de repositionnement (réalisé et idéal)
La précision du repositionnement réalisé sous machine est finalement quantifiée par comparaison des indices de repositionnement et des déplacements correspondant d'une part au repositionnement réalisé et d'autre part au repositionnement optimal reposant un recalage contours (Figure 1).
Les volumes d'intérêt dérivés des images de planification et de suivi peuvent également être comparés. Ils peuvent être différents du fait des incertitudes de contourage, mais aussi d'une fonte tumorale en cours de traitement. Une forte divergence correspondant à cette deuxième cause peut d'ailleurs déclencher éventuellement une procédure de radiothérapie adaptative (replanifications).
L'ensemble de cette approche peut s'intégrer dans une approche globale associant également une évaluation de la dose cumulée dans les différentes structures par recalage élastique.
5.3. Dispositif de mise en œuyre.
La méthode décrite est de préférence mise en œuvre par l'intermédiaire d'un équipement de calcul dédié décrit en relation avec la figure 5. Un tel équipement 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 œuvre au moins une partie du procédé tel que décrit. Dans au moins un mode de réalisation, l'invention est mise en œuvre sous la forme d'une application installée sur un dispositif de planification. Un tel dispositif comprend :
des moyens d'obtention d'au moins une position d'une structure tumorale au sein de ladite image de suivi en tenant compte d'une part de ladite au moins une image de planification et d'au moins une donnée spécifique associée à ladite image de planification ;
des moyens de simulation, à partir de ladite au moins une position de ladite structure tumorale au sein de ladite au moins une image de suivi, et d'au moins une donnée représentative d'un déplacement appliqué audit patient, d'un déplacement de ladite structure tumorale, délivrant au moins une donnée représentative d'un repositionnement réalisé ;
des moyens de calcul, à partir d'au moins une position de ladite structure tumorale au sein de ladite au moins une image de planification et d'au moins une position de ladite structure tumorale au sein de ladite au moins une image de suivi, délivrant au moins une donnée représentative d'un repositionnement optimal ;
des moyens de comparaison de ladite au moins une donnée représentative d'un repositionnement réalisé et de ladite au moins une donnée représentative d'un repositionnement optimal, délivrant au moins un indicateur d'efficience de déplacement.
Ces moyens sont mis en œuvre par l'intermédiaire de modules logiciels et/ou matériels dédiés. Ces modules sont aptes à mettre en œuvre les interfaces de réception et de transmission de données présentes sur ledit dispositif (I et T).
Dans au moins un mode de réalisation, la technique décrite peut être mis en œuvre 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 œuvre au sein d'une console de planification de traitement de radiothérapie.
Claims
REVENDICATIONS
Procédé de détermination d'une efficience d'un déplacement précédemment effectué sur un patient lors d'une séance de radiothérapie, procédé mis en œuvre par u n dispositif de traitement de données à partir d'au moins une image de planification de traitement (lP) et d'au moins une image de suivi de traitement (ls), procédé caractérisé en ce qu'il comprend :
une phase d'obtention (PH I) d'au moins une position (PSTI) d'une structure tumorale (STi) au sein de ladite image de suivi (ls) en tenant compte d'u ne part de ladite au moins une image de planification (lP) et d'au moins une donnée spécifique (Ds) associée à ladite image de planification (lP) ;
une première phase de calcul (PH 2), à partir de ladite au moins une position (PSTI) de ladite structure tumorale au sein de ladite au moins une image de su ivi, et d'au moins une donnée représentative d'un déplacement (DPT) appliqué audit patient, d'un déplacement de ladite structure tumorale, délivrant au moins une donnée représentative d'un repositionnement réalisé (DRRR) ;
une deuxième phase de calcul (PH3), à partir d'au moins une position (PST2) de ladite structure tumorale (STi) au sein de ladite au moins une image de planification (lP) et de ladite d'au moins une position (PSTI) de lad ite structure tumorale au sein de ladite au moins une image de suivi (ls), d'au moins une donnée représentative d'un repositionnement optimal (DRR0) ;
une phase de comparaison (PH4) de ladite au moins une donnée représentative d'un repositionnement réalisé (DRRR) et de ladite au moins une donnée représentative d'un repositionnement optimal (DRR0), délivrant au moins u n indicateur d'efficience de déplacement (lED).
Procédé selon la revendication 1, caractérisé en ce que ladite au moins une donnée spécifique associée à ladite image de planification comprend au moins une donnée représentative d'un contour de ladite structure tumorale.
Procédé selon la revendication 1, caractérisé en ce que ladite phase d'obtention d'au moins une position d'une structure tumorale comprend :
une étape de calcul d'au moins une transformation géométrique à appliquer à ladite image de planification en fonction de ladite image de suivi ;
une étape d'application de ladite transformation géométrique à ladite au moins une donnée spécifique.
Procédé selon la revendication 1, caractérisé en ce que ladite première phase de calcul comprend :
une étape de calcul d'au moins une transformation géométrique à appliquer à ladite image de suivi en fonction d'une donnée représentative d'un déplacement appliqué audit patient ;
au moins une étape de calcul d'au moins un indice de recouvrement réalisé.
Procédé selon la revendication 1, caractérisé en ce que ladite deuxième phase de calcul comprend :
une étape de recalage appliqué sur ladite image de suivi à partir de ladite image de planification ;
au moins une étape de calcul d'au moins un indice de recouvrement optimal.
Dispositif de détermination d'une efficience d'un déplacement précédemment effectué sur un patient lors d'une séance de radiothérapie, à partir d'au moins une image de planification de traitement et d'au moins une image de suivi de traitement, caractérisé en ce qu'il comprend :
des moyens d'obtention d'au moins une position d'une structure tumorale au sein de ladite image de suivi en tenant compte d'une part de ladite au moins une image de planification et d'au moins une donnée spécifique associée à ladite image de planification ;
des moyens de calcul, à partir de ladite au moins une position de ladite structure tumorale au sein de ladite au moins une image de suivi, et d'au moins une
donnée représentative d'un déplacement appliqué audit patient, d'un déplacement de ladite structure tumorale, délivrant au moins une donnée représentative d'un repositionnement réalisé ;
des moyens de calcul, à partir d'au moins une position de ladite structure tumorale au sein de ladite au moins une image de planification et de ladite au moins une position de ladite structure tumorale au sein de ladite au moins une image de suivi, délivrant au moins une donnée représentative d'un repositionnement optimal ;
des moyens de comparaison de ladite au moins une donnée représentative d'un repositionnement réalisé et de ladite au moins une donnée représentative d'un repositionnement optimal, délivrant au moins un indicateur d'efficience de déplacement.
7. 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 détermination d'une efficience d'un déplacement précédemment effectué sur un patient lors d'une séance de radiothérapie, procédé mis en œuvre par un dispositif de traitement de données à partir d'au moins une image de planification de traitement et d'au moins une image de suivi de traitement, procédé caractérisé en ce qu'il comprend :
une phase d'obtention d'au moins une position d'une structure tumorale au sein de ladite image de suivi en tenant compte d'une part de ladite au moins une image de planification et d'au moins une donnée spécifique associée à ladite image de planification ;
une première phase de calcul, à partir de ladite d'au moins une position de ladite structure tumorale au sein de ladite au moins une image de suivi, et d'au moins une donnée représentative d'un déplacement appliqué audit patient,
d'un déplacement de ladite structure tumorale, délivrant au moins une donnée représentative d'un repositionnement réalisé ;
une deuxième phase de calcu l, à partir de ladite d'au moins une position de ladite structure tumorale au sein de ladite au moins une image de planification et de ladite d'au moins une position de ladite structure tumorale au sein de ladite au moins une image de suivi, d'au moins une donnée représentative d'un repositionnement optimal ;
une phase de comparaison de ladite au moins u ne donnée représentative d'un repositionnement réalisé et de ladite au moins une donnée représentative d'un repositionnement optimal, délivrant au moins un indicateur d'efficience de déplacement..
Procédé de calcul à partir d'au moins une position (PST2) d'une structure tumorale (STi) au sein d'au moins une image de planification de traitement (lP) et d'au moins une position (PSTI) d'une structu re tumorale au sein d'au moins une image de suivi (ls) de traitement, d'au moins une donnée représentative d'un repositionnement optimal (DRR0), procédé caractérisé en ce qu'il comprend : une étape de recalage appliqué sur ladite image de suivi à partir de ladite image de planification ;
au moins une étape de calcul d'au moins un indice de recouvrement optimal.
Dispositif de calcul à partir d'au moins une position (PST2) d'u ne structure tumorale (STi) au sein d'au moins une image de planification de traitement (lP) et d'au moins une position (PSTI) d'une structu re tumorale au sein d'au moins une image de suivi (ls) de traitement, d'au moins une donnée représentative d'un repositionnement optimal (DRR0), procédé caractérisé en ce qu'il comprend : des moyens de recalage appliqué su r ladite image de suivi à partir de ladite image de planification ;
des moyens de calcul d'au moins un indice de recouvrement optimal.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1353554A FR3004653A1 (fr) | 2013-04-18 | 2013-04-18 | Procede de controle de qualite de positionnement radiotherapique |
FR1353554 | 2013-04-18 |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2014170490A2 true WO2014170490A2 (fr) | 2014-10-23 |
WO2014170490A3 WO2014170490A3 (fr) | 2014-12-11 |
Family
ID=49474504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/EP2014/058024 WO2014170490A2 (fr) | 2013-04-18 | 2014-04-18 | Procédé de contrôle de qualité de positionnement radiothérapique |
Country Status (2)
Country | Link |
---|---|
FR (1) | FR3004653A1 (fr) |
WO (1) | WO2014170490A2 (fr) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1342488A2 (fr) * | 2002-03-07 | 2003-09-10 | Siemens Aktiengesellschaft | Procédé et dispositif pour le repositionnement relatif répété d'un patient |
US20050096515A1 (en) * | 2003-10-23 | 2005-05-05 | Geng Z. J. | Three-dimensional surface image guided adaptive therapy system |
EP1759733A2 (fr) * | 2005-08-31 | 2007-03-07 | Hitachi, Ltd. | Appareil et méthode de positionnement pour la radiothérapie |
WO2012095172A1 (fr) * | 2011-01-12 | 2012-07-19 | Brainlab | Détermination d'une position de partie corporelle durant une radiothérapie |
-
2013
- 2013-04-18 FR FR1353554A patent/FR3004653A1/fr active Pending
-
2014
- 2014-04-18 WO PCT/EP2014/058024 patent/WO2014170490A2/fr active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1342488A2 (fr) * | 2002-03-07 | 2003-09-10 | Siemens Aktiengesellschaft | Procédé et dispositif pour le repositionnement relatif répété d'un patient |
US20050096515A1 (en) * | 2003-10-23 | 2005-05-05 | Geng Z. J. | Three-dimensional surface image guided adaptive therapy system |
EP1759733A2 (fr) * | 2005-08-31 | 2007-03-07 | Hitachi, Ltd. | Appareil et méthode de positionnement pour la radiothérapie |
WO2012095172A1 (fr) * | 2011-01-12 | 2012-07-19 | Brainlab | Détermination d'une position de partie corporelle durant une radiothérapie |
Non-Patent Citations (1)
Title |
---|
CAZOULAT G ET AL: "Benefice volumetrique de la radiotherapie guidee par l'image dans les cancers prostatiques : marges et cartographies de probabilite de presence", CANCER RADIOTHERAPIE, ELSEVIER, PARIS, FR, vol. 13, no. 5, 1 septembre 2009 (2009-09-01), pages 365-374, XP026497706, ISSN: 1278-3218, DOI: 10.1016/J.CANRAD.2009.06.002 [extrait le 2009-07-29] * |
Also Published As
Publication number | Publication date |
---|---|
FR3004653A1 (fr) | 2014-10-24 |
WO2014170490A3 (fr) | 2014-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6895019B2 (ja) | ディープ畳み込みニューラルネットワークを用いたビームモデルパラメータの決定 | |
JP6886565B2 (ja) | 表面の動きを追跡する方法及び装置 | |
CN104548372B (zh) | 放射治疗的剂量确定装置 | |
CN107206251B (zh) | 用于构建四维图像信息的方法和成像*** | |
JP6053792B2 (ja) | 可変画像レジストレーションのワークフローにおけるユーザー入力と変形ベクトル場の修正との統合 | |
Lawson et al. | Quantitative evaluation of a cone‐beam computed tomography–planning computed tomography deformable image registration method for adaptive radiation therapy | |
US11138768B2 (en) | System and method for artifact reduction in an image | |
WO2018223261A1 (fr) | Système de radiothérapie, procédé de traitement de données et support de stockage | |
CN103785112B (zh) | 用于不规则形状目标的基于图像的检测和跟踪的模板匹配方法 | |
WO2014170489A1 (fr) | Optimisation de la planification en radiothérapie adaptative | |
EP2900323B1 (fr) | Procédé d'évaluation d'une méthode de détermination de la dose cumulée reçue au cours d'un traitement de radiothérapie | |
EP3322347B1 (fr) | Systeme de guidage autonome d'un equipement porte-aiguille | |
WO2014170490A2 (fr) | Procédé de contrôle de qualité de positionnement radiothérapique | |
KR102024993B1 (ko) | 안구의 이미지 분석 방법 | |
BE1025557B1 (fr) | Enregistrement de données de radiographie par faisceau de particules | |
EP2987137B1 (fr) | Procédé de production de données représentatives de contraintes de traitement de radiothérapie, dispositifs et programme correspondant | |
Lee et al. | Prostate brachytherapy seed reconstruction with Gaussian blurring and optimal coverage cost | |
Lozares-Cordero et al. | An open-source development based on photogrammetry for a real-time IORT treatment planning system | |
WO2024127708A1 (fr) | Dispositif de support de radiothérapie et système de radiothérapie | |
US20240245932A1 (en) | Techniques for adaptive radiotherapy based on cbct projection correction and reconstruction | |
US20240249451A1 (en) | Techniques for removing scatter from cbct projections | |
Di Franco | Dosimetric impact of prostate intrafraction motion during hypofractionated prostate cancer radiotherapy | |
Li | Sensitivity of CT Images in Tracking Savi Based APBI HDR Brachytherapy | |
Riboldi et al. | Challenges and opportunities in image guided particle therapy | |
CN118414194A (zh) | 连续性放射治疗计划 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 14718441 Country of ref document: EP Kind code of ref document: A2 |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 14718441 Country of ref document: EP Kind code of ref document: A2 |