System and computer program product for radiation inverse treatment planning

1. A radiation inverse treatment planning system for a linear accelerator, comprising:

a radiation source, configured for delivering individual radiotherapeutic beams, wherein the beam is comprised of a plurality of beamlets, each individual radiotherapeutic beam having a predetermined location and incidence angle inside and/or outside a target area, a size and a shape,

at least a data bus system,

a memory coupled to the data bus system, wherein the memory comprises a computer usable program code, and

a processing unit coupled to the data bus system, wherein the processing unit is configured to execute the computer usable program code to

pre-compute a dictionary composed of a set of individual possible radiotherapeutic beams including, for each individual radiotherapeutic beam its location and its incident angle,

associate a weight to each individual radiotherapeutic beam, based on one or more constraints,

wherein the processing unit executes the computer usable program code to

find the sparsest subset of individual radiotherapeutic weighted beams, so as to satisfy said one or more constraints.

2. The system of claim 1 , the weight associated to each individual radiotherapeutic beam comprising the time of irradiation.

3. The system of claim 1 , the weight associated to each individual radiotherapeutic beam comprising the dose rate.

4. The system of claim 1 , the weight associated to each individual radiotherapeutic beam comprising dose profile.

5. The system of claim 1 ,

a first physical support for said radiation source,

a second physical support arranged for receiving a patient,

wherein the first physical support and the second physical support are arranged to be moved one relative to the other,

and wherein the processing unit executes the computer sable program code to

find the sparsest subset of individual radiotherapeutic beams, so as to satisfy said one or more constraints each time that the first physical support is moved relatively to the second physical support and/or each time that the second physical support is moved relatively to the first physical support.

6. The system of claim 1 ,

wherein the processing unit executes the computer usable program code to

find the sparsest subset of individual radiotherapeutic beams, so as to satisfy said one or more constraints each time that a patient and/or an organ of the patient moves.

7. The system of claim 1 , wherein the processing unit executes the computer usable program code to

find the minimum number of non-zero weights so as to satisfy said one or more constraints.

8. The system of claim 1 , wherein the number of non-zero weights is at least 1/100 of the number of pre-computed individual radiotherapeutic beams.

9. The system of claim 1 , wherein the processing unit executes the computer usable program code to

minimize a weighted L 1 norm of the vector of weights while satisfying said one or more constraints, so as to obtain an optimal subset of individual radiotherapeutic beams.

10. The system of claim 1 , wherein the processing unit executes the computer usable program code to

minimize a weighted L 0 norm of the vector of weights while satisfying said one or more constraints, so as to obtain an optimal subset of individual radiotherapeutic beams.

11. The system of claim 1 , wherein the processing unit executes the computer usable program code to

minimize a weighted L 2 norm of the vector of weights while satisfying said one or more constraints, so as to obtain an optimal subset of individual radiotherapeutic beams.

12. The system of claim 1 , wherein the processing unit executes the computer usable program code to

locate each individual radiotherapeutic beam in a location of a three-dimensional grid.

13. The system of claim 1 , wherein the processing unit is configured to execute the computer usable program code in real-time.

14. The system of claim 1 , wherein the constraint comprises dose constraints applied to the target region and/or to other areas such as sensitive structures to be protected against too high dose radiation.

15. The system of claim 1 , wherein the processing unit is configured to execute the computer usable program code to

take into account the physical properties of the patient's anatomy during the pre-computing of the set of individual radiotherapeutic beams.

16. The system of claim 1 , wherein the processing unit executes the computer usable program code to

apply a convex optimization criterion.

17. The system of the previous claim 16 , said optimization criterion comprising minimizing a treatment time.

18. The system of claim 1 , wherein the radiation source is a linear accelerator.

19. The system of claim 1 , wherein the radiation source is a cobalt source or a proton beam.

20. A computer program product, comprising:

a tangible non-transitory computer usable medium including computer usable program code for a radiation inverse treatment planning system comprising a radiation source configured for delivering individual radiotherapeutic beams, wherein the beam is comprised of a plurality of beamlets, each individual radiotherapeutic beam having a predetermined location and incidence angle inside and/or outside a target area, a size and a shape, the computer usable program code being used to

pre-compute a dictionary composed of a set of individual possible radiotherapeutic beams including for each individual radiotherapeutic beam its location and its incident angle,

associate a weight to each individual radiotherapeutic beam, based on one or more constraints,

wherein a processing unit executes the computer usable program code to

find the sparsest subset of individual radiotherapeutic weighted beams so as to satisfy said one or more constraints.

21. A non-transitory computer data carrier storing presentation content created with a radiation inverse treatment planning method, comprising the following steps:

pre-compute a dictionary composed of a set of individual radiotherapeutic beams including, for each individual radiotherapeutic beam its location and its incident angle, each individual radiotherapeutic beam being generated by a radiation source and having a predetermined location and incidence angle inside and/or outside a target area, a size and a shape, wherein the beam is comprised of a plurality of beamlets,

associate a weight to each individual radiotherapeutic beam, based on one or more constraints,

wherein a processing unit executes the computer usable program code to

find the sparsest subset of individual radiotherapeutic weighted beams so as to satisfy said one or more constraints.