CN110548231A - Radiotherapy plan generation system, radiotherapy plan generation device and storage medium - Google Patents

Radiotherapy plan generation system, radiotherapy plan generation device and storage medium Download PDF

Info

Publication number
CN110548231A
CN110548231A CN201910865915.1A CN201910865915A CN110548231A CN 110548231 A CN110548231 A CN 110548231A CN 201910865915 A CN201910865915 A CN 201910865915A CN 110548231 A CN110548231 A CN 110548231A
Authority
CN
China
Prior art keywords
plan
current
value
sampling point
evaluation index
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910865915.1A
Other languages
Chinese (zh)
Other versions
CN110548231B (en
Inventor
彭雄峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201910865915.1A priority Critical patent/CN110548231B/en
Publication of CN110548231A publication Critical patent/CN110548231A/en
Application granted granted Critical
Publication of CN110548231B publication Critical patent/CN110548231B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/103—Treatment planning systems
    • 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/103—Treatment planning systems
    • A61N5/1031—Treatment planning systems using a specific method of dose optimization

Abstract

The embodiment of the invention discloses a radiotherapy plan generating system, a radiotherapy plan generating device and a storage medium. The system includes a processor to implement the steps of: acquiring an initial radiotherapy plan corresponding to a target area, quantitative evaluation index values of a reference plan of the initial radiotherapy plan, current constraint weights and current dose target values of all sampling points in a dose limiting range corresponding to the target area; optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point to generate an adjusted radiotherapy plan; and if the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is greater than a preset deviation threshold value, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value, and optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value. By the technical scheme, the optimization efficiency of the radiotherapy plan is improved.

Description

Radiotherapy plan generation system, radiotherapy plan generation device and storage medium
Technical Field
The embodiment of the invention relates to medical technology, in particular to a radiotherapy plan generating system, a radiotherapy plan generating device and a storage medium.
Background
In a radiation therapy (radiotherapy) system, a physicist needs to generate an executable radiotherapy plan, such as a radiation field shape and a hop count, by a radiotherapy planning system according to a radiotherapy prescription provided by a doctor, such as a required dose of a target region of a tumor, a maximum dose of an organ at risk, and the like.
In the existing radiotherapy planning system software, there is a complete optimization workflow, and a physicist needs to determine a reasonable constraint condition for generating a radiotherapy plan according to a radiotherapy prescription, then inputs the reasonable constraint condition into the radiotherapy planning system software to generate the radiotherapy plan, and continuously adjusts the constraint condition in the optimization process of the radiotherapy plan until an executable radiotherapy plan meeting the prescription requirement is generated.
In the process of optimizing the radiotherapy plan, a physical engineer is required to continuously intervene, so that the labor is consumed, and the efficiency of generating the executable radiotherapy plan is low.
disclosure of Invention
The embodiment of the invention provides a radiotherapy plan generation system, a radiotherapy plan generation device and a storage medium, which are used for realizing automatic optimization of a radiotherapy plan and improving optimization efficiency of the radiotherapy plan.
In a first aspect, an embodiment of the present invention provides a radiotherapy plan generating system, including a memory, a processor, and a computer program stored on the memory and executable on the processor, where the processor executes the computer program to implement the following steps:
Acquiring an initial radiotherapy plan corresponding to a target area, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limited range corresponding to the target area, wherein the constraint weight is used for representing the proportion of the dose target value corresponding to the sampling point in an objective function, and the dose target value is used for constraint conditions of radiotherapy plan optimization;
Optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generating an adjusted radiotherapy plan;
if the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is larger than a preset deviation threshold value, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so as to enable the deviation to be smaller than or equal to the preset deviation threshold value.
In a second aspect, an embodiment of the present invention further provides a radiotherapy plan generating apparatus, configured in a processor, the apparatus including:
An initial radiotherapy plan acquisition module, configured to acquire an initial radiotherapy plan corresponding to a target region, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limiting range corresponding to the target region, where the constraint weight is used to represent a specific gravity of a dose target value corresponding to the sampling point in an objective function, and the dose target value is used as a constraint condition for radiotherapy plan optimization;
The adjusting radiotherapy plan generating module is used for optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point and generating an adjusting radiotherapy plan;
and the radiotherapy plan optimization module is used for updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point and optimally adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so as to enable the deviation to be less than or equal to the preset deviation threshold value if the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is greater than a preset deviation threshold value.
in a third aspect, an embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program is used to, when executed by a processor, perform a method including:
acquiring an initial radiotherapy plan corresponding to a target area, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limited range corresponding to the target area, wherein the constraint weight is used for representing the proportion of the dose target value corresponding to the sampling point in an objective function, and the dose target value is used for constraint conditions of radiotherapy plan optimization;
optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generating an adjusted radiotherapy plan;
If the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is larger than a preset deviation threshold value, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so as to enable the deviation to be smaller than or equal to the preset deviation threshold value.
the embodiment of the invention obtains an initial radiotherapy plan corresponding to a target area, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limited range corresponding to the target area; optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generating an adjusted radiotherapy plan; if the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is larger than the preset deviation threshold, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value, and optimally adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the deviation is smaller than or equal to the preset deviation threshold. The method and the device realize the driving and adjustment of the current constraint weight and the current dose target value of each sampling point in the dose limiting range according to the deviation between the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value, further automatically optimize and adjust the radiotherapy plan, avoid the process that a physicist repeatedly and manually adjusts constraint conditions in the optimization process of the radiotherapy plan, achieve the effect of customizing the radiotherapy plan by one key, and improve the optimization efficiency of the radiotherapy plan.
Drawings
fig. 1 is a schematic structural diagram of a radiotherapy plan generating system according to a first embodiment of the present invention;
Fig. 2 is a schematic diagram illustrating an implementation process of a processor in a radiotherapy plan generating system according to a first embodiment of the present invention;
fig. 3 is a schematic diagram illustrating an implementation process of a processor in a radiotherapy plan generating system according to a second embodiment of the present invention;
Fig. 4 is a schematic diagram of an implementation process of a processor in a radiotherapy plan generating system according to a third embodiment of the present invention;
fig. 5 is a schematic structural diagram of a radiotherapy plan generating apparatus according to a fourth embodiment of the present invention.
Detailed Description
the present invention will be described in further detail with reference to the accompanying drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the invention and are not limiting of the invention. It should be further noted that, for the convenience of description, only some of the structures related to the present invention are shown in the drawings, not all of the structures.
Example one
the radiotherapy plan generating system provided by the embodiment can be suitable for automatically generating and optimizing a radiotherapy plan. Referring to fig. 1, the radiotherapy plan generation system comprises a memory 101, a processor 102 and a computer program stored on the memory 101 and executable on the processor 102. When the processor 102 executes the computer program, the steps shown in fig. 2 are implemented:
s110, obtaining an initial radiotherapy plan corresponding to the target area, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limiting range corresponding to the target area.
The initial radiotherapy plan is an initial radiotherapy plan prepared according to a radiotherapy prescription prescribed by a doctor, and the initial radiotherapy plan is not an executable radiotherapy plan which is most suitable for clinical requirements but a radiotherapy plan which needs to be optimized. The initial radiotherapy plan may be obtained from an external storage medium or a network, or may be generated by a radiotherapy plan generating system according to a radiotherapy prescription. The plan quantitative evaluation index value is a value of an index capable of quantitatively evaluating the radiotherapy effect of the radiotherapy plan, and may be, for example, a Dose-volume Histogram (DVH) corresponding to the radiotherapy plan, or a Dose value having a certain characteristic (for example, a maximum value or a minimum value) corresponding to the radiotherapy plan. The reference plan quantitative evaluation index value is a referential plan quantitative evaluation index set for automatically optimizing the radiotherapy plan. In the embodiment of the invention, the reference plan quantitative evaluation index value is set to be determined according to the initial radiotherapy plan instead of the high-quality tumor patient plan database, so that the association with the high-quality tumor patient plan database in the plan optimization process is avoided, and the time consumption for establishing the database is saved.
The dose limiting range refers to a pre-defined range in the planning image from which the constraint weights and dose target values of selected sampling points are used to adjust the dose distribution of the planned radiotherapy of the target region and the organs at risk, i.e., the related information of the sampling points in the dose limiting range can be used to optimize the radiotherapy plan. The dose limiting range may be disposed within the target volume, may be disposed outside the target volume (e.g., an annular region surrounding the target volume), or may be disposed in a connected region (e.g., a circular region covering the target volume) that encompasses the target volume, without limitation as to the shape of the dose limiting range.
the constraint weight is a preset weight value corresponding to the sampling point and is used for representing the proportion of a dose target value corresponding to the sampling point in an objective function. The dose target value is a target dose value of radiotherapy corresponding to the sampling point and is used as a constraint condition for optimizing a radiotherapy plan. For example, if the sampling point is within the target region, the dose target value is the dose lower limit value; and if the sampling point is outside the target area, the dose target value is the upper limit value of the dose. The current constraint weight and the current dose target value respectively refer to the constraint weight and the dose target value corresponding to the sampling point in the current operation.
Illustratively, the number of dose limit ranges and the width of the dose limit ranges are determined by at least one of a target volume, a distance between the target volume and each organ at risk, and a preset organ priority for each organ at risk. The preset organ priority weight refers to a preset weight representing the importance of organs at risk, and is set according to clinical requirements. For example, if the clinical requirement is to protect organs at risk preferentially during radiotherapy, the weight is higher. Since the dose limiting range is used to optimize the radiotherapy plan, the setting of the dose limiting range is directly related to the optimization of the radiotherapy plan. The smaller the number and wider the width of dose limiting ranges, the greater the number of points within which the radiotherapy plan can be optimized, and the better the optimization. However, in clinical applications, both effectiveness and efficiency are required, and therefore, the number of dose limit ranges and the width of each dose limit range can be determined comprehensively according to at least one of target region information such as the position, size and shape of the target region, preset organ priority of each organ at risk, and the distance between the target region and each organ at risk. For example, for more accurate planning results, the number of dose limiting ranges may be set to 1 and the width is the entire circular region including the target region, so the number of sampling points for optimizing the planning is larger. In order to improve the calculation efficiency, dose limiting rings with a certain width and distributed at intervals can be arranged outside the target area, for example, one dose limiting ring is distributed at intervals of 2 mm.
the radiotherapy plan optimization in the embodiment of the invention is realized by constraining the whole radiotherapy plan based on the dose target value of each sampling point in the dose limiting range and the corresponding constraint weight on the basis of the dose distribution of the target area and the critical organ corresponding to the initial radiotherapy plan. Specifically, the radiotherapy plan is adjusted according to the constraint weight and the dose target value of each sampling point, so as to determine a deviation between an adjusted plan quantitative evaluation index value corresponding to the adjusted radiotherapy plan and a reference plan quantitative evaluation index value, such as a ratio or a difference between two plan quantitative evaluation values, and then drive the update of the constraint weight and/or the dose target value of each sampling point by the deviation, so as to further iteratively adjust the radiotherapy plan, thereby forming a cyclic automatic optimization of the radiotherapy plan until a convergence condition of the radiotherapy plan optimization is satisfied, such as that the deviation satisfies a set deviation value (i.e., a preset deviation threshold value), or the adjusted iteration number reaches a preset number (i.e., a preset number upper limit), and the like.
To achieve the above radiotherapy plan optimization, an initial radiotherapy plan to be optimized needs to be determined first, and a reference plan quantitative evaluation index value is determined according to the initial radiotherapy plan. For example, when the plan quantitative evaluation index is a dose value, the dose value in the dose distribution corresponding to the target area and/or the organs at risk in the initial radiotherapy plan can be counted; and when the plan quantitative evaluation index is DVH, counting the DVH of each target area and the organs at risk corresponding to the initial radiotherapy plan, and further determining the DVH corresponding to the whole initial radiotherapy plan. In addition, sampling of points within the dose limits is required to determine sampling points, which are not changed during the subsequent planning optimization. For the first optimization, the current constraint weight (which may be referred to as the initial constraint weight) of each sample point may be set to a small weight value; the current dose target value (which may be referred to as the initial dose target value) for each sampling point may be assigned with a dose drop template for a tumor at a particular location, which is a template of attenuation of the dose at different distances from the tumor, which is related to the specific tumor location and tumor type; in addition, the initial dose target value may also be determined according to the dose distribution of the initial radiotherapy plan, that is, the dose value of the dose distribution corresponding to the initial radiotherapy plan at each sampling point is used as the initial dose target value of the sampling point. For subsequent iterations in the optimization process, the current constraint weights and the current dose target values are updated through operation S130.
illustratively, acquiring an initial radiotherapy plan corresponding to the target volume includes: and generating an initial radiotherapy plan corresponding to the target area by taking a standard plan quantitative evaluation index value corresponding to the radiotherapy prescription as a constraint condition and based on a plan optimization algorithm of the region of interest.
The standard planned quantitative evaluation index value refers to an ideal planned quantitative evaluation index value corresponding to a radiotherapy prescription, such as an ideal DVH. The standard quantitative evaluation index value of the plan can be obtained by matching or fitting from a high-quality tumor patient plan database, or can be obtained by generating a radiotherapy plan based on a general organ distribution state and a radiotherapy prescription obtained by clinical statistics and obtaining the radiotherapy plan.
the initial radiotherapy plan may be generated according to optimal constraints of the radiotherapy prescription and then needs to be optimized to be an executable radiotherapy plan. In this embodiment, an initial radiotherapy plan corresponding to a target region is generated by using a plan optimization algorithm of a volume of interest (VOI) with a standard plan quantitative evaluation index value corresponding to a radiotherapy prescription as a constraint condition. The advantage of this setting is that the initial radiotherapy plan with higher precision can be obtained, and then the quantitative evaluation index value of the reference plan with higher precision can be obtained, so as to enhance the optimization effect of the radiotherapy plan.
And S120, optimizing the initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generating an adjusted radiotherapy plan.
The process of optimizing the radiotherapy plan is an iterative process, each iteration adjusts the radiotherapy plan, and the resulting radiotherapy plan is called an adjusted radiotherapy plan. Each iteration is operated as follows: and (3) constructing an objective function by taking the current constraint weight and the current dose target value of each sampling point as constraint conditions, and adjusting the adjusted radiotherapy plan obtained last time (the initial radiotherapy plan in the first iteration) to obtain the adjusted radiotherapy plan.
S130, if the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is larger than a preset deviation threshold value, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the deviation is smaller than or equal to the preset deviation threshold value.
In the embodiment of the invention, the preset deviation threshold is used as one of convergence conditions for optimizing the radiotherapy plan, and after the adjusted radiotherapy plan is obtained in each iteration, the quantitative evaluation index value of the adjusted plan is obtained according to the adjusted radiotherapy plan and is compared with the quantitative evaluation index value of the reference plan to obtain the deviation between the index values. The deviation is then compared to a preset deviation threshold. And if the deviation is less than or equal to the preset deviation threshold value, the iteration is converged, and the adjusted radiotherapy plan at the moment is used as a final radiotherapy plan optimization result. If the deviation is larger than the preset deviation threshold value, the optimization effect of the plan is not up to the standard, and the radiotherapy plan needs to be continuously adjusted. At this time, the current constraint weight and the current dose target value of each sampling point need to be updated according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value. During specific updating, the updating relationship between the constraint weight and the dose target value can be constructed according to the influence of the deviation between the two index values corresponding to each target area and each crisis organ on each sampling point. The form of the update relationship may not be limited as long as the updated data can narrow the gap between the deviation and the preset deviation threshold. When updating the current constraint weights and the current dose target value, only one of the amounts, preferably the current dose target value, may be adjusted. Of course, the current constraint weight and the current dose target value may also be adjusted simultaneously. And then, continuously adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value and the operation of S120 until the obtained deviation is less than or equal to the preset deviation threshold value.
according to the technical scheme of the embodiment, an initial radiotherapy plan corresponding to a target area, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limited range corresponding to the target area are obtained; optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generating an adjusted radiotherapy plan; if the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is larger than the preset deviation threshold, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value, and optimally adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the deviation is smaller than or equal to the preset deviation threshold. The method and the device realize the driving and adjustment of the current constraint weight and the current dose target value of each sampling point in the dose limiting range according to the deviation between the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value, further automatically optimize and adjust the radiotherapy plan, avoid the process that a physicist repeatedly and manually adjusts constraint conditions in the optimization process of the radiotherapy plan, achieve the effect of customizing the radiotherapy plan by one key, and improve the optimization efficiency of the radiotherapy plan.
Example two
in this embodiment, a step of "adjusting the reference plan quantitative evaluation index value" is added on the basis of the first embodiment. On the basis, the step of removing high-dose points in the target area can be further added. Wherein explanations of the same or corresponding terms as those of the above embodiments are omitted. Referring to fig. 3, the implementation procedure of the processor in the radiotherapy plan generating system provided by this embodiment includes:
S210, obtaining an initial radiotherapy plan corresponding to the target area, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limiting range corresponding to the target area.
And S220, optimizing the initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generating an adjusted radiotherapy plan.
And S230, if the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is greater than a preset deviation threshold, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and optimally adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so as to enable the deviation to be less than or equal to the preset deviation threshold.
And S240, if the deviation is greater than a preset deviation threshold value and the current iteration number is greater than or equal to a preset number upper limit, adjusting the reference plan quantitative evaluation index value to reduce the constraint of the adjusted reference plan quantitative evaluation index value on the current constraint weight and the current dose target value of each sampling point.
In this embodiment, in addition to the iterative convergence condition of the preset deviation threshold, a convergence condition of the number of iterations is also set. When both convergence conditions cannot be satisfied, that is, the deviation between the adjustment plan quantitative evaluation index value obtained after each iteration operation and the reference plan quantitative evaluation index value is greater than the preset deviation threshold, and the iteration number is greater than or equal to the preset upper limit, in this embodiment, the reference plan quantitative evaluation index value is dynamically adjusted to relax the iteration constraint, that is, the constraint of the adjusted reference plan quantitative evaluation index value on the current constraint weight and the current dose target value of each sampling point is reduced. According to the above process, as long as both iterative convergence conditions cannot be satisfied, the once-referenced plan quantitative evaluation index value is adjusted until an executable radiotherapy plan satisfying the above convergence conditions is obtained.
Illustratively, adjusting the reference plan quantitative evaluation index value to reduce the constraint of the adjusted reference plan quantitative evaluation index value on the current constraint weight and the current dose target value of each sampling point comprises: and reducing the deviation between the adjusted reference plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value so as to reduce the constraint of the adjusted reference plan quantitative evaluation index value on the current constraint weight and the current dose target value of each sampling point.
the purpose of adjusting the reference plan quantitative evaluation index value is to reduce the deviation between the adjustment plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value. Taking the planned quantitative evaluation index value as the dose volume histogram as an example, in specific implementation, the effect of adjusting the whole histogram can be achieved by adjusting the dose value of each point in the histogram. The adjustment mode is selected according to the fact that the adjustment amplitude of the reference plan quantitative evaluation index value is smaller and smaller, for example, an adjustment step length can be preset, but the adjustment step length needs to meet the requirement of gradual reduction; algorithms such as dichotomy may also be used. The advantage of setting up like this is, can adjust reference plan quantization assessment index value more flexibly dynamically, and then further improve the optimization efficiency of radiotherapy plan.
And S250, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value, and optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the deviation between the adjustment plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value is smaller than or equal to a preset deviation threshold value.
The description of this operation can be referred to in S130, but the reference plan quantitative evaluation index value is changed to the adjusted reference plan quantitative evaluation index value.
and S260, setting the dose limiting range in the target area, taking the maximum dose value as a plan quantitative evaluation index value, referring to the plan quantitative evaluation index value as the maximum dose value corresponding to the adjusted radiotherapy plan optimized in the dose limiting range, and adjusting the plan quantitative evaluation index value as the maximum dose value corresponding to the adjusted radiotherapy plan continuously optimized in the dose limiting range.
If the dose limit range in the first embodiment and the optimization process of the radiotherapy plan in S210 to S250 is a dose limit ring outside the target region, and the dose target value is a dose upper limit value, the result of adjusting the radiotherapy plan by using the dose upper limit value of the sampling point in the dose limit ring as a constraint is to remove a high dose region in the dose limit ring, and the dose distribution of the radiotherapy plan obtained in this way in the target region may have some high dose points with an excessively high dose value. Therefore, in this embodiment, after the adjusted radiotherapy plan is obtained in S250, the dose limiting range is further set in the target volume, and the maximum dose value is used as the plan quantification assessment index value, so that the obtained adjusted radiotherapy plan is continuously optimized through the operations in S210 to S250, thereby removing high dose points in the target volume. In this case, the reference plan quantitative evaluation index value may be determined as a maximum dose value corresponding to the optimized adjusted radiotherapy plan within the dose limiting range, and may be obtained by counting the dose values of the sampling points selected within the dose limiting range, for example; the adjustment plan quantitative evaluation index value may be set to a maximum dose value corresponding to the adjustment radiotherapy plan which is continuously optimized within the dose limiting range, and may be obtained by, for example, counting the dose values of the respective sampling points within the dose limiting range in the dose distribution corresponding to the adjustment radiotherapy plan which is continuously optimized after the adjustment radiotherapy plan which is continuously optimized is obtained by continuously optimizing the adjustment radiotherapy plan. Thereafter, the process proceeds to S270 to S280.
S270, determining the current constraint weight and the current dose target value of each sampling point in the dose limiting range, and continuously optimizing and adjusting the radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point.
As the dose allowance range is altered, the current constraint weight and the current dose target value for each sampling point within the dose allowance range need to be re-determined. Like S110, the initial constraint weight may be set to be a small weight value, and the initial dose target value may be assigned by a dose drop template of the tumor at a specific location, or determined according to the dose distribution of the optimized adjusted radiotherapy plan obtained in S250, that is, the dose value of the dose distribution corresponding to the optimized adjusted radiotherapy plan at each sampling point is used as the initial dose target value of the sampling point. The current constraint weights and current dose target values in subsequent iterations are similarly updated by subsequent operations S280. It should be noted that the dose limiting range in this operation is within the target range, so the dose target value is the lower dose limiting value.
after the initial constraint weight and the initial dose target value of each sampling point are determined, an objective function is constructed by taking the initial constraint weight and the initial dose target value as constraint conditions, and the optimized adjusted radiotherapy plan obtained in the step S250 is adjusted to obtain an adjusted radiotherapy plan which is continuously optimized.
s280, if the adjustment plan quantitative evaluation index value corresponding to the continuously optimized adjustment radiotherapy plan is greater than or equal to the reference plan quantitative evaluation index value, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and continuously optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the adjustment plan quantitative evaluation index value is smaller than the reference plan quantitative evaluation index value.
and counting the dose value corresponding to the continuously optimized adjusted radiotherapy plan corresponding to each sampling point in the dose limiting range to determine an adjusted plan quantitative evaluation index value, and comparing the adjusted plan quantitative evaluation index value with the reference plan quantitative evaluation index value determined in the step S260. If the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan to be continuously optimized is greater than or equal to the reference plan quantitative evaluation index value, which indicates that a point with an excessively high dose value exists in the dose limit range at this time, in step S130, the current constraint weight and the current dose target value of each sampling point are updated by using the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value, and the iterative operation of the objective function is continuously performed by using the updated current constraint weight and the current dose target value of each sampling point, so as to continuously optimize and adjust the radiotherapy plan, and the above process is repeated until the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan to be continuously optimized is smaller than the reference plan quantitative evaluation index value. Thus, an executable radiotherapy plan with more uniform dose distribution in the target area can be obtained.
According to the technical scheme of the embodiment, if the deviation is greater than the preset deviation threshold value and the current iteration number is greater than or equal to the preset number upper limit, the reference plan quantitative evaluation index value is adjusted; and updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value, and optimally adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the deviation between the adjustment plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value is smaller than or equal to a preset deviation threshold value. The automatic adjustment of the quantitative evaluation index value of the reference plan is realized, so that the automatic optimization efficiency of the radiotherapy plan is further improved. Setting a dose limited range in a target area, quantitatively evaluating an index value by taking a maximum dose value as a plan, determining a current constraint weight and a current dose target value of each sampling point in the dose limited range, and continuously optimizing and adjusting a radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point; if the adjustment plan quantitative evaluation index value corresponding to the continuously optimized adjustment radiotherapy plan is greater than or equal to the reference plan quantitative evaluation index value, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and continuously optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the adjustment plan quantitative evaluation index value is smaller than the reference plan quantitative evaluation index value. The secondary optimization of the radiotherapy plan is realized, and the uniformity of the dose distribution in the target area is improved on the basis that the dose distribution of the target area in the obtained radiotherapy plan meets the radiotherapy prescription, so that the optimization effect of the radiotherapy plan is further improved.
EXAMPLE III
In the present embodiment, based on the first embodiment, further optimization is performed on "update of the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value". Wherein explanations of the same or corresponding terms as those of the above embodiments are omitted. In this embodiment, a dose volume histogram is taken as an example of a plan quantitative evaluation index. Referring to fig. 4, the implementation process of the processor in the radiotherapy plan generating system provided by this embodiment includes:
s310, obtaining an initial radiotherapy plan corresponding to the target area, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limiting range corresponding to the target area.
And S320, optimizing the initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generating an adjusted radiotherapy plan.
s330, if the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is larger than a preset deviation threshold, determining a comprehensive action result when the plan quantitative index deviation of each organ at risk acts on the sampling point for each sampling point, and taking the comprehensive action result as an organ deviation influence value of the sampling point.
The radiotherapy plan contains the dose distribution for each target volume and its surrounding organs at risk. The reference plan quantitative assessment index value corresponding to the initial radiotherapy plan also includes the reference plan quantitative assessment index value of each target region and each organ at risk. Similarly, after the generation of the adjusted radiotherapy plan, the dose distribution of each target and the organs at risk can be determined, and then the quantitative evaluation index value of the adjusted plan of the corresponding target and the corresponding organs at risk can be determined according to the dose distribution of each target and organ at risk. Therefore, the deviation of the plan quantitative evaluation index corresponding to the corresponding organs at risk can be calculated according to the reference plan quantitative evaluation index value and the adjustment plan quantitative evaluation index value corresponding to each organ at risk; similarly, the deviation of the quantitative assessment plan corresponding to each target region can be calculated according to the quantitative assessment reference plan value and the quantitative assessment adjustment plan value corresponding to each target region.
considering that the dose setting purposes of the target region and the organs at risk in radiotherapy planning are different, that is, the irradiation dose of the target region is higher and the irradiation dose of the organs at risk is lower, in the embodiment, when the current constraint weight and the current dose target value of the sampling point are updated, the deviation of the planning quantization index corresponding to the target region and the deviation of the planning quantization index corresponding to the organs at risk are distinguished. In addition, according to the description of the foregoing embodiment, the update of the constraint weight and the dose target value of the sampling point is realized by the above deviation, so that in this embodiment, it is necessary to determine the result of the deviation of the planned quantization index corresponding to each target area and each organ at risk acting on the sampling point, and further update the constraint weight and the dose target value.
In specific implementation, the results of the plan quantitative index deviation corresponding to each organ at risk acting on the sampling point are accumulated, so that the action results of the plan quantitative index deviations of all organs at risk acting on the sampling point can be determined and used as the organ deviation influence value of the sampling point. According to the process, the organ deviation influence value of each sampling point can be determined.
illustratively, determining the combined effect of the planned quantitative index deviation of each organ at risk on the sampling points includes: and determining a comprehensive action result when the deviation of the planned quantization index of each organ at risk acts on the sampling point according to the organ distance weight between the sampling point and each organ at risk, the deviation of the planned quantization index corresponding to each organ at risk and the priority weight of each preset organ.
Wherein, the organ distance weight refers to a weight value corresponding to the organ at risk, which is used for calculating the influence of the deviation of the planned quantitative index of the organ at risk on the sampling point, and which is a monotonically decreasing function of the distance between the sampling point and the organ at risk, and the organ distance weight α (r) is determined by, for example, the following formula 5ij)：
equation 5:
Wherein b is a constant coefficient for adjusting the distance weight with the distance rijThe decay rate of (c).
In particular, given the different distances between organs at risk and the target area, and thus the different organs at risk are protected to different degrees, the pre-set organ priority (which may also be the same) for each organ at risk may be set according to clinical needs. Meanwhile, considering that the effect result of the plan quantitative index deviation of the organs at risk different from the distance between the sampling points on the sampling points is different, an organ distance weight is set for each organ at risk in the embodiment. Therefore, the comprehensive action result of the plan quantitative index deviations of all organs at risk on the same sampling point can be determined according to the organ distance weight, the preset organ priority weight and the plan quantitative index deviation of each organ at risk.
for example, the comprehensive action result of the deviation of the planned quantitative index of each organ at risk when the deviation acts on the sampling point is determined according to the following formula 1:
Equation 1:
Wherein, Oirepresenting organ deviation influence values of the ith sampling point; sorganrepresenting a set of organs at risk in a constraint; r isijrepresenting the distance between the ith sample point and the jth organ at risk; alpha (r)ij) Is a distance rijrepresents the distance weight between the ith sampling point and the jth interested area, and alpha (r) when the interested area is an organ at riskij) Representing an organ distance weight;Indicating a deviation of the adjusted planned quantitative assessment index value corresponding to the jth organ-at-risk from the reference planned quantitative assessment index value,The ratio of the quantitative evaluation index value of the adjustment plan to the quantitative evaluation index value of the reference plan;representing the degree of deviation of the histogram deviation, in the range of 0-1,a is a constant coefficient and is used for adjusting the change sensitivity of the deviation degree;to the extent of deviationand a distance rijIs monotonousA decreasing function representing that the histogram deviation of the jth organ-at-risk acts on the single organ deviation contribution of the ith sample point; phi (j) represents a preset organ priority weight.
and S340, aiming at each sampling point, determining a comprehensive action result when the planned quantization index deviation of each target area acts on the sampling point, and taking the comprehensive action result as a target area deviation influence value of the sampling point.
In S330, the results of the deviation of the planned quantization index corresponding to each target area acting on the sampling point may be accumulated, and the action results of the deviation of the planned quantization index of all target areas at the sampling point may be determined as the target area deviation influence value of the sampling point. According to the process, the target deviation influence value of each sampling point can be determined.
Illustratively, determining a combined effect of the planned quantitative index deviation for each target zone acting on the sampling points comprises: and determining a comprehensive action result when the planned quantization index deviation of each target area acts on the sampling point according to the target area distance weight between the sampling point and the target area, the planned quantization index deviation corresponding to the target area and the priority weight of each preset target area.
the target distance weight is a concept corresponding to the organ distance weight, and refers to a weighted value corresponding to the target, which is used to calculate the influence of the deviation of the planned quantization index of the target on the sampling point, and is a monotonically decreasing function of the distance between the sampling point and the target, and the target distance weight is determined, for example, by the above formula 5. The preset target area priority weight refers to a preset weight representing the importance of the target area, and is set according to clinical requirements. For example, if the clinical requirement is to preferentially ensure the radiotherapy effect of the target region during radiotherapy, the weight value is larger. According to the description of the preset organ priority weight and the preset target area priority weight, the values of the two weights can be set according to the clinical requirement for preferentially protecting organs at risk or preferentially ensuring the radiotherapy effect of the target area. Illustratively, the preset target region priority weight is greater than the preset organ priority weight. This has the advantage that the amplitude of the repetitive oscillations between the target and the organs at risk during optimization of the radiotherapy plan can be reduced.
In specific implementation, considering that the action results of the planned quantization index deviations of the target regions with different distances from the sampling points on the sampling points are different, a target region distance weight is set for each target region in the embodiment. Therefore, the comprehensive action result of the plan quantization index deviations of all the target areas on the same sampling point can be determined according to the target area distance weight of each target area, the preset target area priority weight and the plan quantization index deviations.
for example, the comprehensive action result when the deviation of the planned quantization index of each target region acts on the sampling point is determined according to the following formula 2:
Equation 2:
wherein, TiRepresenting the target deviation influence value of the ith sampling point; stargeta set of target areas in the constraint representation; alpha (r)ij) Representing a target range weight; r isijrepresenting the distance between the ith sampling point and the jth target area;A deviation of an adjusted planned quantitative assessment index value corresponding to the jth target zone from a reference planned quantitative assessment index value,A ratio of the reference plan quantitative evaluation index value to the adjustment plan quantitative evaluation index value;Representing the degree of deviation of the histogram deviation, in the range of 0-1,a is a constant coefficient and is used for adjusting the change sensitivity of the deviation degree;to the extent of deviationAnd a distance rijThe monotone decreasing function of (a) represents that the histogram deviation of the jth target region acts on the single target region deviation influence value of the ith sampling point; phi (j) represents the preset target area priority weight.
The histogram deviation for the jth organ-at-risk or target as described above acts on the deviation contribution of the individual organ-at-risk or target at the ith sample pointCan be determined by the following equation 6:
Equation 6:
wherein c is a constant coefficient for adjusting the deviation influence value with the distance rijThe decay rate of (c).
And S350, adjusting the current dose target value of the corresponding sampling point according to the current dose target value, the preset dose adjusting amplitude, the organ deviation influence value and the target region deviation influence value of each sampling point.
The preset dose adjustment range refers to an adjustment range of the dose target value, and is an adjustment coefficient preset according to clinical requirements.
For each sampling point, the process of adjusting its current dose target value is: determining a comprehensive deviation influence value of the sampling point according to the organ deviation influence value and the target area deviation influence value corresponding to the sampling point; then determining the dose adjustment value of the sampling point according to the preset dose adjustment amplitude and the comprehensive deviation influence value; and finally, calculating to obtain a new current dose target value according to the current dose target value and the determined dose adjusting value.
For example, when the dose limit range is a dose limit loop outside the target region, the current dose target value of each sampling point can be adjusted by the following equation 3:
Equation 3: d'i＝Di·[1+γ1·(-Oi+Ti)]
Wherein, D'iIndicating the adjusted current dose target value, DiRepresenting the current dose target value before adjustment; gamma ray1Indicating a preset dose adjustment amplitude, by default gamma may be set1Is 0.5.
And S360, adjusting the current constraint weight of the corresponding sampling point according to the current constraint weight, the preset weight adjustment amplitude, the organ deviation influence value and the target area deviation influence value of each sampling point.
The preset weight adjustment range refers to an adjustment range of the constraint weight, and is an adjustment coefficient preset according to clinical requirements.
For each sample point, the process of adjusting the current constraint weight is as follows: determining a comprehensive deviation influence value of the sampling point according to the organ deviation influence value and the target area deviation influence value corresponding to the sampling point; then determining a constraint weight adjustment value of the sampling point according to a preset weight adjustment amplitude and the comprehensive deviation influence value; and finally, calculating to obtain a new current constraint weight according to the current constraint weight and the determined constraint weight adjusting value.
For example, when the dose limit range is a dose limit loop outside the target zone, the current constraint weight of each sampling point can be adjusted by the following equation 4:
Equation 4: w'i＝Wi+γ2·(-Oi+Ti)
Wherein, W'irepresenting the adjusted current constraint weight, WiRepresenting the current constraint weight before adjustment; gamma ray2Indicating a preset weight adjustment amplitude, by default gamma may be set2Is 500.
And S370, optimally adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so as to enable the deviation to be smaller than or equal to a preset deviation threshold value.
It should be noted that S360 and S370 may be executed alternatively or simultaneously.
According to the technical scheme, a comprehensive action result when the plan quantization index deviation of each organ at risk acts on the sampling point is determined and is used as an organ deviation influence value of the sampling point; determining a comprehensive action result when the plan quantization index deviation of each target area acts on the sampling point, and taking the comprehensive action result as a target area deviation influence value of the sampling point; adjusting the current dose target value of the corresponding sampling point according to the current dose target value, the preset dose adjustment amplitude, the organ deviation influence value and the target region deviation influence value of each sampling point; and/or adjusting the current constraint weight of the corresponding sampling point according to the current constraint weight, the preset weight adjustment amplitude, the organ deviation influence value and the target deviation influence value of each sampling point. More accurate adjustment of the current constraint weight and/or the current dose target value is realized, so that the optimization efficiency of the radiotherapy plan is further improved.
example four
The present embodiment provides a radiotherapy plan generating apparatus, configured in a processor, and referring to fig. 5, the radiotherapy plan generating apparatus specifically includes:
an initial radiotherapy plan obtaining module 510, configured to obtain an initial radiotherapy plan corresponding to the target region, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limiting range corresponding to the target region, where the constraint weight is used to represent a proportion of a dose target value corresponding to the sampling point in an objective function, and the dose target value is used as a constraint condition for radiotherapy plan optimization;
An adjusted radiotherapy plan generating module 520, configured to optimize an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generate an adjusted radiotherapy plan;
A radiotherapy plan optimizing module 530, configured to update the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value if a deviation between the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan is greater than a preset deviation threshold, so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and optimize and adjust the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point, so as to make the deviation less than or equal to the preset deviation threshold.
Optionally, on the basis of the foregoing apparatus, the apparatus further includes a reference plan quantitative evaluation index value adjustment module, configured to:
Optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and after generating an adjusted radiotherapy plan, if the deviation is greater than a preset deviation threshold value and the current iteration number is greater than or equal to a preset number upper limit, adjusting a reference plan quantitative evaluation index value to reduce the constraint of the adjusted reference plan quantitative evaluation index value on the current constraint weight and the current dose target value of each sampling point;
And updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value, and optimally adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the deviation between the adjustment plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value is smaller than or equal to a preset deviation threshold value.
Further, the reference plan quantitative evaluation index value adjustment module is specifically configured to:
And reducing the deviation between the adjusted reference plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value so as to reduce the constraint of the adjusted reference plan quantitative evaluation index value on the current constraint weight and the current dose target value of each sampling point.
Optionally, the radiotherapy plan optimization module 530 comprises a numerical update sub-module for:
aiming at each sampling point, determining a comprehensive action result of the plan quantization index deviation of each organ at risk when the plan quantization index deviation acts on the sampling point, and taking the comprehensive action result as an organ deviation influence value of the sampling point;
Aiming at each sampling point, determining a comprehensive action result when the planned quantitative index deviation of each target area acts on the sampling point, and taking the comprehensive action result as a target area deviation influence value of the sampling point;
Adjusting the current dose target value of the corresponding sampling point according to the current dose target value, the preset dose adjustment amplitude, the organ deviation influence value and the target region deviation influence value of each sampling point; and/or the presence of a gas in the gas,
And adjusting the current constraint weight of the corresponding sampling point according to the current constraint weight, the preset weight adjustment amplitude, the organ deviation influence value and the target region deviation influence value of each sampling point.
further, the value update submodule is specifically configured to:
determining a comprehensive action result when the deviation of the planned quantization index of each organ at risk acts on the sampling point according to the organ distance weight between the sampling point and each organ at risk, the deviation of the planned quantization index corresponding to each organ at risk and the priority weight of each preset organ;
Determining a comprehensive action result when the deviation of the planned quantization index of each target area acts on the sampling point comprises the following steps:
and determining a comprehensive action result when the planned quantization index deviation of each target area acts on the sampling point according to the target area distance weight between the sampling point and the target area, the planned quantization index deviation corresponding to the target area and the priority weight of each preset target area.
wherein the preset target area priority weight is greater than the preset organ priority weight.
Optionally, the initial radiotherapy plan acquisition module 510 is specifically configured to:
And generating an initial radiotherapy plan corresponding to the target area by taking a standard plan quantitative evaluation index value corresponding to the radiotherapy prescription as a constraint condition and based on a plan optimization algorithm of the region of interest.
optionally, the number of dose limit ranges and the width of the dose limit ranges are determined by at least one of a target volume, a distance between the target volume and each organ at risk, and a preset organ priority for each organ at risk.
optionally, on the basis of the above apparatus, the apparatus further includes a radiotherapy plan continuous optimization module, configured to:
Optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the deviation is less than or equal to a preset deviation threshold value, setting a dose limit range in a target area, and taking the maximum dose value as a plan quantitative evaluation index value, then referring to the plan quantitative evaluation index value as the maximum dose value corresponding to the optimized adjusted radiotherapy plan in the dose limit range, and adjusting the plan quantitative evaluation index value as the maximum dose value corresponding to the adjusted radiotherapy plan which is continuously optimized in the dose limit range;
determining the current constraint weight and the current dose target value of each sampling point in the dose limiting range, and continuously optimizing and adjusting the radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point;
If the adjustment plan quantitative evaluation index value corresponding to the continuously optimized adjustment radiotherapy plan is greater than or equal to the reference plan quantitative evaluation index value, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and continuously optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the adjustment plan quantitative evaluation index value is smaller than the reference plan quantitative evaluation index value.
By the radiotherapy plan generating device, the current constraint weight and the current dose target value of each sampling point in the dose limit range are driven and adjusted according to the deviation between the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value, so that the radiotherapy plan is automatically optimized and adjusted, the process that a physicist repeatedly and manually adjusts constraint conditions in the optimization process of the radiotherapy plan is avoided, the effect of customizing the radiotherapy plan by one key is achieved, and the optimization efficiency of the radiotherapy plan is improved.
the radiotherapy plan generating device provided by the embodiment of the invention can execute the radiotherapy plan generating method provided by any embodiment of the invention, and has corresponding functional modules and beneficial effects of the executing method.
It should be noted that, in the embodiment of the radiotherapy plan generating apparatus, the included units and modules are only divided according to functional logic, but are not limited to the above division as long as the corresponding functions can be realized; in addition, specific names of the functional units are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the present invention.
EXAMPLE five
The present embodiments provide a storage medium containing computer executable instructions which, when executed by a computer processor, are operable to perform a radiotherapy plan generation method, the method comprising:
Acquiring an initial radiotherapy plan corresponding to a target area, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limited range corresponding to the target area, wherein the constraint weight is used for representing the proportion of the dose target value corresponding to the sampling point in a target function, and the dose target value is used for constraint conditions of radiotherapy plan optimization;
optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generating an adjusted radiotherapy plan;
if the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is larger than the preset deviation threshold, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and optimally adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the deviation is smaller than or equal to the preset deviation threshold.
Of course, the storage medium provided by the embodiment of the present invention contains computer executable instructions, and the computer executable instructions are not limited to the above method operations, and may also perform related operations in the radiotherapy plan generation method provided by any embodiment of the present invention.
from the above description of the embodiments, it is obvious for those skilled in the art that the present invention can be implemented by software and necessary general hardware, and certainly, can also be implemented by hardware, but the former is a better embodiment in many cases. Based on such understanding, the technical solutions of the present invention or portions thereof that contribute to the prior art may be embodied in the form of a software product, where the computer software product may be stored in a computer-readable storage medium, such as a floppy disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a FLASH Memory (FLASH), a hard disk, or an optical disk of a computer, and includes several instructions to enable a computer device (which may be a personal computer, a server, or a network device) to execute the radiotherapy plan generating method provided in the embodiments of the present invention.
it is to be noted that the foregoing is only illustrative of the preferred embodiments of the present invention and the technical principles employed. It will be understood by those skilled in the art that the present invention is not limited to the particular embodiments illustrated herein, but is capable of various obvious changes, rearrangements and substitutions as will now become apparent to those skilled in the art without departing from the scope of the invention. Therefore, although the present invention has been described in greater detail by the above embodiments, the present invention is not limited to the above embodiments, and may include other equivalent embodiments without departing from the spirit of the present invention, and the scope of the present invention is determined by the scope of the appended claims.

Claims (11)

1. A radiotherapy plan generation system comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, comprising:
the processor, when executing the computer program, implements the steps of:
Acquiring an initial radiotherapy plan corresponding to a target area, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limited range corresponding to the target area, wherein the constraint weight is used for representing the proportion of the dose target value corresponding to the sampling point in an objective function, and the dose target value is used for constraint conditions of radiotherapy plan optimization;
Optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generating an adjusted radiotherapy plan;
If the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is larger than a preset deviation threshold value, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so as to enable the deviation to be smaller than or equal to the preset deviation threshold value.
2. The system of claim 1, further comprising, after optimizing an initial radiotherapy plan and generating an adjusted radiotherapy plan based on the current constraint weights and the current dose target values for each of the sampling points:
if the deviation is greater than the preset deviation threshold value and the current iteration number is greater than or equal to a preset number upper limit, adjusting the reference plan quantitative evaluation index value to reduce the constraint of the adjusted reference plan quantitative evaluation index value on the current constraint weight and the current dose target value of each sampling point;
Updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value, and optimally adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so that the deviation between the adjustment plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value is less than or equal to the preset deviation threshold value.
3. the system of claim 2, wherein adjusting the reference plan quantitative evaluation index value to reduce the constraint of the adjusted reference plan quantitative evaluation index value on the current constraint weight and the current dose target value of each of the sampling points comprises:
And reducing the deviation between the adjusted reference plan quantitative evaluation index value and the adjusted reference plan quantitative evaluation index value so as to reduce the constraint of the adjusted reference plan quantitative evaluation index value on the current constraint weight and the current dose target value of each sampling point.
4. The system of claim 1, wherein updating the current constraint weight and the current dose target value for each of the sampling points as a function of the adjusted planned quantitative assessment index value and the reference planned quantitative assessment index value comprises:
aiming at each sampling point, determining a comprehensive action result when the plan quantization index deviation of each organ at risk acts on the sampling point, and taking the comprehensive action result as an organ deviation influence value of the sampling point;
Determining a comprehensive action result of the planned quantization index deviation of each target area acting on each sampling point as a target area deviation influence value of the sampling point;
adjusting the current dose target value of the corresponding sampling point according to the current dose target value, the preset dose adjustment amplitude, the organ deviation influence value and the target region deviation influence value of each sampling point; and/or the presence of a gas in the gas,
And adjusting the current constraint weight of the corresponding sampling point according to the current constraint weight, the preset weight adjustment amplitude, the organ deviation influence value and the target deviation influence value of each sampling point.
5. the system of claim 4, wherein determining a composite effect of planned quantifier deviations for each of the organs-at-risk on the sample points comprises:
Determining a comprehensive action result when the plan quantization index deviation of each organ at risk acts on the sampling point according to the organ distance weight between the sampling point and each organ at risk, the plan quantization index deviation corresponding to each organ at risk and the priority weight of each preset organ;
determining a combined effect result when the planned quantitative index deviation of each target region acts on the sampling point, the combined effect result comprising:
and determining a comprehensive action result when the planned quantization index deviation of each target area acts on the sampling point according to the target area distance weight between the sampling point and each target area, the planned quantization index deviation corresponding to the target area and the priority weight of each preset target area.
6. The system of claim 4 or 5, wherein the preset target zone priority weight is greater than the preset organ priority weight.
7. The system of claim 1, wherein acquiring an initial radiotherapy plan corresponding to the target volume comprises:
And generating an initial radiotherapy plan corresponding to the target area based on a plan optimization algorithm of the region of interest by taking a standard plan quantitative evaluation index value corresponding to a radiotherapy prescription as a constraint condition.
8. the system of claim 1, wherein the number of dose limit ranges and the width of the dose limit ranges are determined by at least one of the target volume, a distance between the target volume and each organ-at-risk, and a preset organ priority weight for each organ-at-risk.
9. the system of claim 1, further comprising, after optimally adjusting the radiotherapy plan to make the deviation less than or equal to the preset deviation threshold according to the updated current constraint weights and current dose target values of each of the sampling points:
Setting a dose limit range in the target area, and taking a maximum dose value as a plan quantitative evaluation index value, wherein the reference plan quantitative evaluation index value is a maximum dose value corresponding to an adjusted radiotherapy plan optimized in the dose limit range, and the adjusted plan quantitative evaluation index value is a maximum dose value corresponding to an adjusted radiotherapy plan continuously optimized in the dose limit range;
Determining the current constraint weight and the current dose target value of each sampling point in the dose limiting range, and continuously optimizing the adjusted radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point;
if the adjustment plan quantitative evaluation index value corresponding to the continuously optimized adjustment radiotherapy plan is greater than or equal to the reference plan quantitative evaluation index value, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and continuously optimizing the adjustment radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so as to enable the adjustment plan quantitative evaluation index value to be smaller than the reference plan quantitative evaluation index value.
10. A radiotherapy plan generating apparatus, the apparatus being configured in a processor, the apparatus comprising:
An initial radiotherapy plan acquisition module, configured to acquire an initial radiotherapy plan corresponding to a target region, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limiting range corresponding to the target region, where the constraint weight is used to represent a specific gravity of a dose target value corresponding to the sampling point in an objective function, and the dose target value is used as a constraint condition for radiotherapy plan optimization;
The adjusting radiotherapy plan generating module is used for optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point and generating an adjusting radiotherapy plan;
And the radiotherapy plan optimization module is used for updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point and optimally adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so as to enable the deviation to be less than or equal to the preset deviation threshold value if the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is greater than a preset deviation threshold value.
11. A computer-readable storage medium having a computer program stored thereon, wherein the computer-executable instructions, when executed by a computer processor, are for performing a method comprising:
Acquiring an initial radiotherapy plan corresponding to a target area, a reference plan quantitative evaluation index value corresponding to the initial radiotherapy plan, and a current constraint weight and a current dose target value of each sampling point in a dose limited range corresponding to the target area, wherein the constraint weight is used for representing the proportion of the dose target value corresponding to the sampling point in an objective function, and the dose target value is used for constraint conditions of radiotherapy plan optimization;
Optimizing an initial radiotherapy plan according to the current constraint weight and the current dose target value of each sampling point, and generating an adjusted radiotherapy plan;
If the deviation between the adjustment plan quantitative evaluation index value corresponding to the adjustment radiotherapy plan and the reference plan quantitative evaluation index value is larger than a preset deviation threshold value, updating the current constraint weight and the current dose target value of each sampling point according to the adjustment plan quantitative evaluation index value and the reference plan quantitative evaluation index value so as to adjust the current constraint weight and/or the current dose target value of each sampling point, and optimizing and adjusting the radiotherapy plan according to the updated current constraint weight and the current dose target value of each sampling point so as to enable the deviation to be smaller than or equal to the preset deviation threshold value.
CN201910865915.1A 2019-09-12 2019-09-12 Radiotherapy plan generation system, radiotherapy plan generation device and storage medium Active CN110548231B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910865915.1A CN110548231B (en) 2019-09-12 2019-09-12 Radiotherapy plan generation system, radiotherapy plan generation device and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910865915.1A CN110548231B (en) 2019-09-12 2019-09-12 Radiotherapy plan generation system, radiotherapy plan generation device and storage medium

Publications (2)

Publication Number Publication Date
CN110548231A true CN110548231A (en) 2019-12-10
CN110548231B CN110548231B (en) 2021-11-05

Family

ID=68740317

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910865915.1A Active CN110548231B (en) 2019-09-12 2019-09-12 Radiotherapy plan generation system, radiotherapy plan generation device and storage medium

Country Status (1)

Country Link
CN (1) CN110548231B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110993059A (en) * 2019-12-11 2020-04-10 上海联影医疗科技有限公司 Maximum dose point optimization method and device, electronic equipment and storage medium
CN111388879A (en) * 2020-03-19 2020-07-10 上海联影医疗科技有限公司 Radiation dose determination system, radiation dose determination device and storage medium
CN112057753A (en) * 2020-09-23 2020-12-11 上海联影医疗科技股份有限公司 Radiotherapy plan adjusting system and device
CN112149669A (en) * 2020-09-23 2020-12-29 上海联影医疗科技股份有限公司 Radiotherapy plan generation system, radiotherapy plan generation device and storage medium
CN116130056A (en) * 2023-04-17 2023-05-16 福建自贸试验区厦门片区Manteia数据科技有限公司 Determination device for radiotherapy plan and electronic equipment

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090154644A1 (en) * 2007-12-12 2009-06-18 Varian Medical Systems International Ag Treatment planning system and method for radiotherapy
CN104519956A (en) * 2012-06-01 2015-04-15 光线搜索实验室公司 A method and a system for optimizing a radiation therapy treatment plan
CN105031819A (en) * 2015-08-25 2015-11-11 上海联影医疗科技有限公司 Dosage optimization method and system
CN105447330A (en) * 2015-12-30 2016-03-30 上海联影医疗科技有限公司 Weight adjustment method and weight adjustment device of intensity-modulated radiation therapy
CN106730411A (en) * 2016-12-21 2017-05-31 上海联影医疗科技有限公司 Radiotherapy treatment planning optimization method and system
WO2018031365A1 (en) * 2016-08-09 2018-02-15 University Of Maryland, Baltimore System and method for optimizing a treatment plan for irradiation therapy
WO2018048575A1 (en) * 2016-09-07 2018-03-15 Elekta, Inc. System and method for learning models of radiotherapy treatment plans to predict radiotherapy dose distributions
EP3384962A1 (en) * 2017-04-05 2018-10-10 RaySearch Laboratories AB System and method for automatic radiotherapy treatment planning
CN109453473A (en) * 2018-12-30 2019-03-12 上海联影医疗科技有限公司 Radiotherapy treatment planning system, determining device and storage medium
CN109499012A (en) * 2018-12-27 2019-03-22 苏州雷泰医疗科技有限公司 A kind of TPS automatic Iterative optimization algorithm of optimization dosage guidance

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090154644A1 (en) * 2007-12-12 2009-06-18 Varian Medical Systems International Ag Treatment planning system and method for radiotherapy
WO2009074288A1 (en) * 2007-12-12 2009-06-18 Varian Medical Systems International Ag Treatment planning system and method for radiotherapy
CN104519956A (en) * 2012-06-01 2015-04-15 光线搜索实验室公司 A method and a system for optimizing a radiation therapy treatment plan
CN105031819A (en) * 2015-08-25 2015-11-11 上海联影医疗科技有限公司 Dosage optimization method and system
CN105447330A (en) * 2015-12-30 2016-03-30 上海联影医疗科技有限公司 Weight adjustment method and weight adjustment device of intensity-modulated radiation therapy
WO2018031365A1 (en) * 2016-08-09 2018-02-15 University Of Maryland, Baltimore System and method for optimizing a treatment plan for irradiation therapy
WO2018048575A1 (en) * 2016-09-07 2018-03-15 Elekta, Inc. System and method for learning models of radiotherapy treatment plans to predict radiotherapy dose distributions
CN106730411A (en) * 2016-12-21 2017-05-31 上海联影医疗科技有限公司 Radiotherapy treatment planning optimization method and system
EP3384962A1 (en) * 2017-04-05 2018-10-10 RaySearch Laboratories AB System and method for automatic radiotherapy treatment planning
CN109499012A (en) * 2018-12-27 2019-03-22 苏州雷泰医疗科技有限公司 A kind of TPS automatic Iterative optimization algorithm of optimization dosage guidance
CN109453473A (en) * 2018-12-30 2019-03-12 上海联影医疗科技有限公司 Radiotherapy treatment planning system, determining device and storage medium

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110993059A (en) * 2019-12-11 2020-04-10 上海联影医疗科技有限公司 Maximum dose point optimization method and device, electronic equipment and storage medium
CN111388879A (en) * 2020-03-19 2020-07-10 上海联影医疗科技有限公司 Radiation dose determination system, radiation dose determination device and storage medium
CN111388879B (en) * 2020-03-19 2022-06-14 上海联影医疗科技股份有限公司 Radiation dose determination system, radiation dose determination device and storage medium
CN112057753A (en) * 2020-09-23 2020-12-11 上海联影医疗科技股份有限公司 Radiotherapy plan adjusting system and device
CN112149669A (en) * 2020-09-23 2020-12-29 上海联影医疗科技股份有限公司 Radiotherapy plan generation system, radiotherapy plan generation device and storage medium
CN112057753B (en) * 2020-09-23 2022-08-16 上海联影医疗科技股份有限公司 Radiotherapy plan adjusting system and device
CN112149669B (en) * 2020-09-23 2023-05-23 上海联影医疗科技股份有限公司 Radiotherapy plan generation system, radiotherapy plan generation device and storage medium
CN116130056A (en) * 2023-04-17 2023-05-16 福建自贸试验区厦门片区Manteia数据科技有限公司 Determination device for radiotherapy plan and electronic equipment
CN116130056B (en) * 2023-04-17 2023-08-22 福建自贸试验区厦门片区Manteia数据科技有限公司 Determination device for radiotherapy plan and electronic equipment

Also Published As

Publication number Publication date
CN110548231B (en) 2021-11-05

Similar Documents

Publication Publication Date Title
CN110548231B (en) Radiotherapy plan generation system, radiotherapy plan generation device and storage medium
US6222905B1 (en) Irradiation dose calculation unit, irradiation dose calculation method and recording medium
US11020613B2 (en) Physician directed radiation treatment planning
CN110420396B (en) Radiotherapy plan optimization system, radiotherapy plan optimization device and storage medium
EP3347094B1 (en) Knowledge based multi-criteria optimization for radiotherapy treatment planning
US9889317B2 (en) System and method for automated radiation treatment planning using physical objectives
US11224762B2 (en) System and method for automatic radiotherapy treatment planning
CN111128340A (en) Radiotherapy plan generation device, radiotherapy plan generation apparatus, and storage medium
CN110931107B (en) Radiotherapy plan generation system, radiotherapy plan generation device and storage medium
CN110706780B (en) Radiotherapy plan generation system and storage medium
US20240115883A1 (en) Neural network calibration for radiotherapy
Zhang et al. The minimum knowledge base for predicting organ-at-risk dose–volume levels and plan-related complications in IMRT planning
CN116341182A (en) Proton intensity-modulated multi-objective optimization method and system based on cooperation of multi-objective evolution and traditional optimization method
Liu et al. Isodose feature‐preserving voxelization (IFPV) for radiation therapy treatment planning
CN112149669B (en) Radiotherapy plan generation system, radiotherapy plan generation device and storage medium
CN112151146A (en) Flux map optimization system, flux map optimization device, and storage medium
CN114420252A (en) Method, device and medium for determining intensity modulated radiotherapy plan evaluation parameter matrix
Mathews et al. Enhanced optimization of volumetric modulated arc therapy plans using Monte Carlo generated beamlets
CN111986777A (en) Radiation therapy plan optimization system, device and storage medium
CN115966281B (en) Method, device, equipment and storage medium for generating arc radiotherapy plan
US10967199B2 (en) Knowledge-based brachytherapy planning based on primitive features computed from anatomical contour data
CN115223683B (en) Radiation treatment plan evaluation method, radiation treatment plan evaluation device, computer equipment and medium
CN117116421B (en) Method and device for determining radiotherapy plan
WO2023184178A1 (en) Treatment plan generation method and apparatus, and device, and medium
US20230117544A1 (en) Radiation therapy treatment planning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

GR01 Patent grant
GR01 Patent grant