Trajectory optimization using dose estimation and conflict detection
Systems and methods for radiation treatment planning can include a computing system determining an estimate of a radiation dose distribution within an anatomical region of a patient, and determining a cost matrix representing an objective function, using the estimate of the radiation dose distribution. The computing system can project the cost matrix on each of a plurality of fluence planes. Each of the plurality of fluence planes can be associated with a corresponding gantry-couch orientation of a plurality of gantry-couch orientations of a medical linear accelerator. The computing system can determine, using projections of the cost matrix on each of the plurality of fluence planes, a sequence of gantry-couch orientations among the plurality of gantry-couch orientations representing a treatment path.
1. A method of radiation treatment planning comprising:
determining, by one or more processors, an estimate of a radiation dose distribution within an anatomical region of a patient;
determining, by the one or more processors, using the estimate of the radiation dose distribution, an objective function defined within the anatomical region of the patient in terms of the estimate of the radiation dose distribution and patient specific data;
determining, by the one or more processors, for each radiation beam orientation of a plurality of radiation beam orientations, a corresponding conflict severity value using the objection function, the corresponding conflict severity value indicative of a severity of violating one or more clinical goals if the anatomical region of the patient is radiated according to each radiation beam orientation of the plurality of radiation beam orientations; and
determining, by the one or more processors, using conflict severity values corresponding to the plurality of radiation beam orientations, a sequence of radiation beam orientations among the plurality of radiation beam orientations representing a treatment path.
2. The method of claim 1 , wherein determining the estimate of the radiation dose distribution within the anatomical region of the patient includes determining the estimate of the radiation dose distribution as a function of a distance from a planning target volume (PTV) of the anatomical region of the patient.
3. The method of claim 1 , wherein the objective function reflects one or more radiation constraints for the patient.
4. The method of claim 1 , wherein the objective function is defined to optimize an intensity-modulated radiation therapy (IMRT) based radiation plan.
5. The method of claim 1 , wherein the objective function is defined to optimize a volumetric-modulated arc therapy (VMAT) based radiation plan.
6. The method of claim 1 , further comprising:
determining the plurality of radiation beam orientations by discretizing a space of possible radiation beam orientations, each point of the discretized space of possible radiation beam orientations represents a corresponding gantry-couch radiation beam orientation of the plurality of radiation beam orientations.
7. The method of claim 1 , wherein determining, for each radiation beam orientation of the plurality of radiation beam orientations, the corresponding conflict severity value includes applying, for each radiation beam orientation, a weighting function to values of the objective function at voxels of the anatomical region of the patient.
8. The method of claim 7 , wherein applying the weighting function includes weighting the values of the objective function at the voxels of the anatomical region of the patient according to a depth relative to a planning target volume (PTV) inside the anatomical region of the patient in a direction along each radiation beam orientation of the plurality of radiation beam orientations.
9. The method of claim 1 , wherein determining the sequence of radiation beam orientations includes:
generating, by the one or more processors, a matrix of the conflict severity values corresponding to the plurality of radiation beam orientations among the plurality of radiation beam orientations; and
determining the sequence of radiation beam orientations among the plurality of radiation beam orientations using the matrix of the conflict severity values corresponding to the plurality of radiation beam orientations.
10. The method of claim 9 , wherein determining the sequence of radiation beam orientations among the plurality of radiation beam orientations includes determining a sequence of conflict severity values within the matrix having a minimum total sum, each conflict severity value of the sequence of conflict severity values corresponding to a radiation beam orientation of the sequence of radiation beam-orientations among the plurality of radiation beam orientations.
11. A radiation treatment planning system comprising:
one or more processors; and
a non-transitory memory to store computer code instructions, the computer code instructions when executed cause the one or more processors to:
determine an estimate of a radiation dose distribution within an anatomical region of a patient;
determine, using the estimate of the radiation dose distribution, an objective function defined in terms of the estimate of the radiation dose distribution and patient specific data;
determine, for each radiation beam orientation of a plurality of radiation beam orientations, a corresponding conflict severity value using the objection function, the corresponding conflict severity value indicative of a severity of violating one or more clinical goals if the anatomical region of the patient is radiated according to each radiation beam orientation of the plurality of radiation beam orientations; and
determine, using conflict severity values corresponding to the plurality of radiation beam orientations, a sequence of radiation beam orientations among the plurality of radiation beam orientations representing a treatment path.
12. The radiation treatment planning system of claim 11 , wherein in determining the estimate of radiation dose distribution within the anatomical region of the patient, the one or more processors are configured to determine the estimate of the radiation dose distribution as a function of a distance from a planning target volume (PTV) of the anatomical region of the patient.
13. The radiation treatment planning system of claim 11 , wherein the objective function reflects one or more radiation constraints for the patient.
14. The radiation treatment planning system of claim 11 , wherein the objective function is defined to optimize an intensity modulated radiation therapy (IMRT) based radiation plan or to optimize a volumetric-modulated arc therapy (VMAT) based radiation plan.
15. The radiation treatment planning system of claim 11 , wherein the one or more processors are further configured to:
determine the plurality of radiation beam orientations by discretizing a space of possible radiation beam orientations, each point of the discretized space of possible radiation beam orientations represents a corresponding radiation beam orientation of the plurality of radiation beam orientations.
16. The radiation treatment planning system of claim 11 , wherein in determining, for each radiation beam orientation of the plurality of radiation beam orientations, the corresponding conflict severity value, the one or more processors are configured to apply, for each radiation beam orientation, a weighting function to values of the objective function at voxels of the anatomical region of the patient.
17. The radiation treatment planning system of claim 16 , wherein in applying the weighting function, the one or more processors are configured to weight the values of the objective function at the voxels of the anatomical region of the patient according to a depth relative to a planning target volume (PTV) inside the anatomical region of the patient in a direction along each radiation beam orientation of the plurality of radiation beam orientations.
18. The radiation treatment planning system of claim 11 , wherein in determining the sequence of radiation beam orientations among the plurality of radiation beam orientations, the one or more processors are configured to:
generate a matrix of the conflict severity values corresponding to the plurality of radiation beam orientations; and
determine the sequence of radiation beam orientations among the plurality of radiation beam orientations using the matrix of the conflict severity values corresponding to the plurality of radiation beam orientations.
19. The radiation treatment planning system of claim 18 , wherein in determining the sequence of radiation beam orientations among the plurality of radiation beam orientations, the one or more processors are configured to determine within the matrix a sequence of conflict severity values having a minimum total sum, each conflict severity value of the sequence of conflict severity values corresponding to a radiation beam orientation of the sequence of radiation beam orientations among the plurality of radiation beam orientations.
20. A non-transitory computer-readable medium including computer code instructions stored thereon, the computer code instructions when executed cause one or more processors to:
determine an estimate of a radiation dose distribution within an anatomical region of a patient;
determine, using the estimate of the radiation dose distribution, an objective function defined within the anatomical region of the patient in terms of the estimate of the radiation dose distribution and patient specific data;
determine, for each radiation beam orientation of a plurality of radiation beam orientations, a corresponding conflict severity value using the objection function, the corresponding conflict severity value indicative of a severity of violating one or more clinical goals if the anatomical region of the patient is radiated according to each radiation beam orientation of the plurality of radiation beam orientations; and
determine, using conflict severity values corresponding to the plurality of radiation beam orientations, a sequence of radiation beam orientations among the plurality of radiation beam orientations representing a treatment path.