CN111986778A - Intensity modulated plan optimization system, device and storage medium - Google Patents
Intensity modulated plan optimization system, device and storage medium Download PDFInfo
- Publication number
- CN111986778A CN111986778A CN202010763340.5A CN202010763340A CN111986778A CN 111986778 A CN111986778 A CN 111986778A CN 202010763340 A CN202010763340 A CN 202010763340A CN 111986778 A CN111986778 A CN 111986778A
- Authority
- CN
- China
- Prior art keywords
- field
- angle
- sub
- portal
- alternative
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 178
- 230000005855 radiation Effects 0.000 claims description 91
- 230000006870 function Effects 0.000 claims description 39
- 238000000034 method Methods 0.000 claims description 31
- 230000004907 flux Effects 0.000 claims description 23
- 238000012216 screening Methods 0.000 claims description 20
- 230000003068 static effect Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 4
- 241000218228 Humulus Species 0.000 description 124
- 238000001959 radiotherapy Methods 0.000 description 9
- 244000025221 Humulus lupulus Species 0.000 description 7
- 235000008694 Humulus lupulus Nutrition 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000002721 intensity-modulated radiation therapy Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 244000169660 Flemingia strobilifera Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Urology & Nephrology (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Surgery (AREA)
- General Engineering & Computer Science (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
The embodiment of the invention discloses an intensity modulated plan optimization system, an intensity modulated plan optimization device and a storage medium. The system comprises: one or more processors, a storage device for storing instructions, the processors for executing the instructions to implement the steps of: obtaining a plurality of alternative field angles, determining the shape and the hop count of the sub-field under each alternative field angle, determining an optimization model based on the shape and the hop count of the sub-field under each alternative field angle, performing iterative optimization on the optimization model to obtain a target field angle and the shape and the hop count of the sub-field under the target field angle, optimizing the shape and the hop count of the sub-field under the target field angle, and determining an optimized intensity modulation plan. The optimal field angle and the shape and hop count of the sub-field under the optimal field angle can be automatically screened out from a plurality of alternative field angles, the optimized intensity modulation plan is obtained, the labor cost is reduced, and the reliability of the intensity modulation plan is improved.
Description
Technical Field
The embodiment of the invention relates to an intensity modulated radiotherapy technology, in particular to an intensity modulated plan optimization system, an intensity modulated plan optimization device and a storage medium.
Background
Intensity Modulated Radiation Therapy (IMRT) is one type of three-dimensional conformal radiation therapy. Intensity modulated radiation therapy is based on three-dimensional conformal radiation, and adjusts the output dose of points in the radiation field section according to the required mode, and makes the radiation dose in the space distribution in vivo consistent with the pathological changes through rotary radiation, so as to form a high dose region.
In the current intensity modulation planning and optimization process, the user needs to set the irradiation angle in advance, which means that the selection of the number of the radiation fields and the angles is generally based on the experience of the user. However, the conditions of patients are different, and the optimal number of the radiation fields and the radiation field angle are not completely the same. Therefore, the current intensity modulation plan is easy to have the situation that the setting of the field angle is not good, and needs to be improved.
Disclosure of Invention
The embodiment of the invention provides an intensity modulation plan optimization system, an intensity modulation plan optimization device and a storage medium, which are used for selecting an optimal field angle from alternative angles according to different conditions so as to obtain an intensity modulation plan conforming to the actual condition.
In a first aspect, an embodiment of the present invention provides an intensity modulated plan optimization system, including: a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the storage device is for storing instructions, the processor is configured to perform the steps of:
acquiring a plurality of alternative field angles, and determining the shape and hop count of the sub-field under each alternative field angle;
determining an optimization model based on the shapes of the segments and the hop counts at each of the alternative portal angles, wherein the optimization model includes optimizing an objective function and a regular term, the regular term being determined according to the portals and/or the segments at each of the alternative portal angles, wherein the portals include at least one of the segments;
performing iterative optimization on the optimization model, and screening a target field angle from the alternative field angles according to an iterative optimization result;
and optimizing the shape and hop count of the sub-field under the target field angle, and determining an optimized intensity modulation plan.
In a second aspect, an embodiment of the present invention further provides an emphasis plan optimization apparatus configured in the processor for emphasis plan optimization, where the apparatus includes:
the device comprises a sub-field and hop count determining module, a shape determining module and a hop count determining module, wherein the sub-field and hop count determining module is used for acquiring a plurality of alternative radiation field angles and determining the shape and hop count of the sub-field under each alternative radiation field angle;
an optimization model determining module, configured to determine an optimization model based on the shapes of the segments and the hop counts at each of the candidate portal angles, where the optimization model includes optimizing an objective function and a regular term, the regular term is determined according to the portal and/or the segments at each of the candidate portal angles, and the portal includes at least one segment;
the target radiation field angle screening module is used for performing iterative optimization on the optimization model and screening a target radiation field angle from the alternative radiation field angles according to an iterative optimization result;
and the intensity modulation plan optimization module is used for optimizing the shape and the hop count of the sub-field under the target radiation field angle and determining an optimized intensity modulation plan.
In a third aspect, an embodiment of the present invention further provides a storage medium containing computer-executable instructions, which are configured by a computer processor to perform the following method:
acquiring a plurality of alternative field angles, and determining the shape and hop count of the sub-field under each alternative field angle;
determining an optimization model based on the shapes of the sub-fields and the hop counts under the alternative radiation field angles, wherein the optimization model comprises optimizing an objective function and a regular term, the regular term is determined according to the radiation field and/or the sub-fields under the alternative radiation field angles, and the radiation field comprises at least one sub-field;
performing iterative optimization on the optimization model, and screening a target field angle from the alternative field angles according to an iterative optimization result;
and optimizing the shape and hop count of the sub-field under the target field angle, and determining an optimized intensity modulation plan.
The technical solution provided by this embodiment is to design an emphasis plan optimization system, where the system includes one or more processors and a storage device, the storage device is used to store instructions, and the processors are used to execute the instructions to implement: obtaining a plurality of alternative field angles, determining the shape and the hop count of the sub-field under each alternative field angle, determining an optimization model based on the shape and the hop count of the sub-field under each alternative field angle, performing iterative optimization on the optimization model to screen a target field angle from the alternative field angles and determine the sub-field and the hop count under the target field angle, further optimizing the shape and the hop count of the sub-field under the target field angle, and determining an optimized intensity modulation plan. Compared with the prior art, the problem that the field angle is not matched with the actual condition or does not meet the actual application requirement due to the fact that the field angle is set according to the experience of the user is needed, the optimal field angle can be automatically screened out from a plurality of alternative field angles, the shape and the hop count of the sub-field under the optimal field angle are optimized, the optimized intensity modulation plan is obtained, the labor cost is reduced, and the reliability of the intensity modulation plan is improved.
Drawings
Fig. 1 is a schematic structural diagram of an intensity modulated plan optimization system according to an embodiment of the present invention;
fig. 2 is a schematic flowchart of an intensity modulated plan optimization method according to an embodiment of the present invention;
fig. 3 is a schematic flowchart of an intensity modulated plan optimization method according to a second embodiment of the present invention;
fig. 4 is a schematic structural diagram of an intensity modulated plan optimization apparatus according to a third 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
Fig. 1 is a schematic structural diagram of an intensity modulated plan optimization system provided in an embodiment of the present invention, where this embodiment is applicable to a situation of optimizing an intensity modulated plan, as shown in fig. 1, the system includes: a processor 10 and a memory device 20.
Wherein, the processor 10 may be one or more, the storage device 20 is configured to store instructions, and the processor 10 is configured to execute the instructions to implement the steps of the emphasis plan optimization method as described in fig. 2:
s110, acquiring a plurality of alternative field angles, and determining the shapes and hop counts of the sub-fields under the alternative field angles.
The candidate radiation field angle may be obtained by equally dividing one isocenter plane into a plurality of angles, or may be obtained by equally dividing a plurality of isocenter planes into a plurality of angles, that is, the candidate radiation field angle is coplanar or non-coplanar, and the isocenter plane may be understood as a plane passing through the isocenter and the radiation emission plane. Alternatively, one or more isocenter planes may be equally divided into 36 candidate portal angles, or one or more isocenter planes may be divided into other angles according to an intensity modulation plan, and the shape and hop count of the segments at each candidate portal angle may be determined.
Optionally, the method for determining the shape and the hop count of the sub-field at each candidate radiation field angle is a flux map optimization method, for example: determining a flux optimization variable in the initial intensity modulation plan; at each alternative radiation field angle, performing minimization calculation on a flux objective function determined according to the flux optimization variable; and (4) carrying out blade serialization processing on the minimization calculation result of the flux objective function, and determining the sub-field and the hop count under the alternative field angle. Wherein, the flux optimization variable can be ray energy and is determined according to the intensity modulation plan optimization. Specifically, the flux optimization variable may be understood as a radiation dose distribution, the flux objective function may be determined according to a target dose distribution and an actual dose distribution, and an expression of the flux objective function is as follows: obj ═ g (d, d)0) Wherein d is the actual dose distribution and d0The flux target function is iteratively adjusted and minimized calculation is carried out on the flux target function by adjusting the flux under each alternative radiation field angle to obtain a minimized calculation result, and the minimized calculation result can be expressed in a vector form; further, the minimization calculation result in the form of vector is converted into a flux data matrix through a sequence conversion program, and each alternative emission is determined according to the flux data matrixThe number of segments and hops under the field angle.
In other embodiments, the method for determining the shape and the hop count of the sub-field under each alternative radiation field angle may further include: and determining the shape and hop count of the sub-field under each alternative radiation field angle according to a direct machine parameter optimization method. The specific method comprises the following steps: determining the shape and the hop count (MU for short) of the sub-field as variables in the optimization process, establishing an objective function, and adjusting the shape and the hop count of the sub-field under each alternative radiation field angle to adjust the corresponding actual dose distribution so as to perform iterative adjustment on the objective function; and when the objective function reaches a minimum value, the shape and the hop count of the corresponding subfield are the optimization result. It should be noted that both direct machine parameter optimization and flux map optimization are prior art in the field, and are not described herein again.
Optionally, the sub-fields include static sub-fields and/or dynamic sub-fields, and the method for determining the sub-fields under each alternative radiation field angle may further be: and acquiring a dynamic plan and/or a static plan in the initial intensity modulation plan, and determining the sub-fields under the alternative field angles according to the dynamic plan and/or the static plan. Dynamic planning means that the position of mechanical parts of the radiotherapy device changes during the beam-out process; static planning means that the position of the mechanical parts of the radiotherapy apparatus does not change during the beam-out process. By way of example only and not to limit the scope of the present invention, the mechanical components of the radiation therapy device may be, for example, the collimator, the leaves of a multi-leaf collimator, or the gantry, and the change in position may include movement of the leaves, rotation of the collimator or rotation of the gantry, and the like. Wherein the dynamic plan and the static plan may be preset. For example, if the static plan is included in the initial intensity modulated plan, 10 to 20 sub-fields may be determined at each of the alternative field angles, if the dynamic plan is included in the initial intensity modulated plan, 30 to 180 sub-fields may be determined at each of the alternative field angles, and if the dynamic plan and the static plan are included in the initial intensity modulated plan, an appropriate number of sub-fields may be generated at the alternative field angles according to the static plan and the dynamic plan. The purpose of optimizing the mixed intensity modulation plan with the non-fixed angle is achieved by setting a static plan and a dynamic plan, determining the number of the sub-fields of each alternative radiation field angle under the static plan and the dynamic plan, and deleting and selecting the sub-fields under the static plan and the dynamic plan.
And S120, determining an optimization model based on the shapes and the hop counts of the sub-fields under the alternative radiation field angles.
The optimization model comprises the step of optimizing an objective function and a regular term, wherein the objective function is determined by the leaf position, the portal jump number, the sub-portal weight and the target dose distribution under each alternative portal angle, the regular term is determined according to the portal and/or the sub-portal, and the portal comprises at least one sub-portal under the corresponding portal angle. The number of the wild hops and the weight of the sub-field are related to the number of the sub-field hops. The regularization term may be specifically determined according to the hop counts of all the sub-fields under any alternative portal angle and/or the weight of any sub-field within any alternative portal angle, where the weight of any sub-field within any alternative portal angle may be understood as a ratio of the hop count of any sub-field within any alternative portal angle to the sum of the hop counts of all the sub-fields within the alternative portal angle (i.e., the hop count of the portal under the corresponding portal angle).
And S130, performing iterative optimization on the optimization model under each alternative radiation field angle, and screening a target radiation field angle from the alternative radiation field angles according to an iterative optimization result.
Optionally, the shape and the hop count of the sub-field at each alternative radiation field angle may be adjusted, and the optimization model is iteratively optimized until the optimization model reaches a set optimization threshold, and the iterative operation is ended. And adjusting the shape of the sub-field under each alternative radiation field angle, namely adjusting the position of the blade under each alternative radiation field angle. Optionally, the target field angles are coplanar or non-coplanar, and the alternative field angle and the target field angle are coplanar or non-coplanar. In one embodiment, the alternative field angles are coplanar such that the preferred target field angles determined therefrom are also coplanar, enabling coplanar treatment; in one embodiment, the alternative field angles are not coplanar but the target field angles determined therefrom may also be coplanar, thereby also achieving coplanar treatment; in one embodiment, the alternative radiation field angles are non-coplanar and the target radiation field angles determined therefrom are also non-coplanar, enabling non-coplanar radiation therapy.
Specifically, the actual dose distribution of the objective function is adjusted according to the blade position in the initial intensity modulation plan, the hop count in the candidate portal angle and the weight of any sub-portal in any candidate portal angle, the regular term is adjusted until the value of the optimization model reaches the set optimization threshold, the iterative optimization operation is ended, and the candidate portal angle when the value of the optimization model reaches the set optimization threshold is used as the target portal angle. For example, in the optimization process, when the value of the optimization model obtained in the current iteration is smaller than the value of the optimization model obtained in the previous iteration by less than the threshold, the optimization model is considered to reach the minimum value.
S140, optimizing the shape and the hop count of the sub-field under the target field angle, and determining the optimized intensity modulation plan.
It can be understood that, after the target portal angle is determined, the shape and the hop count of the sub-portal under the target portal angle can be obtained, the actual dose distribution corresponding to the objective function in the optimization model is adjusted again according to the leaf position and the hop count of the sub-portal under the target portal angle until the value of the optimization model reaches the set optimization threshold, the iterative optimization operation is ended, the shape and the target hop count of the target sub-portal under the target portal angle when the set optimization threshold is reached are determined, and the optimized intensity modulation plan is determined according to the shape and the target hop count of the target sub-portal.
The technical solution provided by this embodiment is to design an emphasis plan optimization system, where the system includes one or more processors and a storage device, the storage device is used to store instructions, and the processors are used to execute the instructions to implement: obtaining a plurality of alternative field angles, determining the shape and the hop count of the sub-field under each alternative field angle, determining an optimization model based on the shape and the hop count of the sub-field under each alternative field angle, performing iterative optimization on the optimization model to screen a target field angle from the alternative field angles and determine the sub-field and the hop count under the target field angle, further optimizing the shape and the hop count of the sub-field under the target field angle, and determining an optimized intensity modulation plan. Compared with the problem that the field angle is not matched with the actual condition due to the fact that the field angle needs to be set according to the experience of the user in the prior art, the optimal field angle and the shape and the hop count of the sub-field under the optimal field angle can be automatically screened out from a plurality of alternative field angles, the optimized intensity modulation plan is obtained, the labor cost is reduced, and the reliability of the intensity modulation plan is improved; moreover, the number of the sub-fields can be reduced as much as possible by the embodiment, and the complexity of the radiation therapy planning can be reduced.
Example two
Fig. 3 is a flowchart illustrating a method for optimizing a emphasis plan executed by the processor 10 according to a second embodiment of the present invention, which may be a variation of the second embodiment, and details steps in the second embodiment are described. Optionally, the screening a target radiation field angle from the candidate radiation field angles according to an iterative optimization result includes: and taking the alternative portal angle with the hop count exceeding a set threshold value as the target portal angle, and deleting the alternative portal angle with the hop count not exceeding the set threshold value. In the method, reference is made to the above-described embodiments for those parts which are not described in detail. Referring specifically to fig. 3, the method may include the steps of:
s210, acquiring a plurality of alternative field angles, and determining the shapes and hop counts of the sub-fields under the alternative field angles.
And S220, determining an optimization model based on the shapes and the hop counts of the sub-fields under the alternative radiation field angles.
As described in the previous embodiments, the optimization model includes optimizing an objective function and a regularization term, the objective function is determined by the leaf position, the number of field hops, the subfield weight and the target dose distribution at each alternative field angle, and the field hops number and the subfield weight are related to the subfield hops number. The regular term is determined according to the hop count of the sub-fields, the regular term comprises a first regular term and a second regular term, the first regular term is determined according to the hop count of the fields within the candidate field angle, and the second regular term is determined according to the weight of the sub-fields within the candidate field angle. Correspondingly, the hop count of the portal within the candidate portal angle is determined according to the hop counts of all the sub-portals within the candidate portal angle, and the weight of the sub-portals is determined according to the hop count of any sub-portal within the candidate portal angle and the hop counts of all the sub-portals within the candidate portal angle (i.e. the hop counts of the portal).
Optionally, the determining method of the optimization model includes: calculating the sum of the hop counts of each sub-field in any optional field angle (namely the hop count of the field), calculating the ratio of the hop count of any sub-field in any optional field angle to the sum of the hop counts, multiplying the sum of the hop counts of all optional field angles by a corresponding first weight value to obtain a first regular term, multiplying the sum of the ratio of all optional field angles by a corresponding second weight value to obtain a second regular term, adding the first regular term and the second regular term to obtain a regular term, and determining an optimization model according to the regular term and a target function.
And the sum of the ratios of the hop counts of all the sub-fields in any one alternative portal to the sum of the hop counts is 1, namely the sum of the ratios of the hop counts of all the sub-fields to the hop counts of the corresponding portal is 1 under each alternative portal angle. Hop count b in ith candidate field angle in this embodimenti(hop count of the portal) is the sum of the hop counts of all the sub-portals within the candidate portal angle, i.e. bi＝∑j mijWhere j is the subfield at the ith candidate field angle, mijThe hop count of the jth sub-field under the ith alternative field angle; the weight of the sub-field at the ith candidate portal angle is the ratio of the hop count of any sub-field at the ith candidate portal angle to the sum of the hop counts of all sub-fields at the candidate portal angle, that is, the expression of the weight of the sub-field at the ith candidate portal angle is as follows: w is aij＝mij/biAnd isj wij1. (ii) a The expression for the first regular term is α1∑i biThe expression of the second regular term is α2∑ij wijWherein α is1、α2Respectively representing a first weight value and a second weight valueAnd (4) weighting values.
Further, the first regular term and the second regular term are added to obtain a regular term, that is, the expression of the regular term is alpha1∑i bi+α2∑ij wij. In this embodiment, the regular term and the objective function may be added to obtain an expression of the optimization model, where the expression of the optimization model is: g (d (x, b, w), d0)+α1∑i bi+α2∑ij wijWherein g (d (x, b, w), d0) Is the objective function, α1∑i biIs a first regularization term, α2∑ij wijIs a second regularization term, x is the leaf position, b is the hop count of the portal within each alternative portal angle, w is the weight of any sub-portal within any alternative portal angle, b is the weight of any sub-portal within any alternative portal angleiIs the sum of the hops of all the sub-fields within the ith candidate field angle, wijIs the weight of the jth sub-field within the ith candidate field angle, d (x, b, w) is the actual dose distribution when the leaf position is x, the hop count of the field is b and the sub-field weight is w, d0Is the target dose distribution.
And S230, performing iterative optimization on the optimization model under each alternative portal angle, taking the alternative portal angle with the hop number exceeding the set threshold as the target portal angle, and deleting the alternative portal angle with the hop number not exceeding the set threshold.
Specifically, the leaf position (i.e., the shape of the sub-field), the number of hops within the alternative field angle, and the weight of any sub-field within any alternative field angle in the initial emphasis plan may be continually optimized, in combination with ∑j wijDynamically determining an optimization model g (d (x, b, w), d ═ 10)+α1∑i bi+α2∑ij wijUntil the value of the optimization model is minimum, obtaining a minimum optimization model: minx，b，w g(d(x，b，w)，d0)+α1∑i bi+α2∑ij wijEnding the iterative adjustment operation and optimizing the alternative radiation field angle corresponding to the model with the minimum valueAs the target field angle. It should be noted that, in this embodiment, by determining the optimization model according to the regular term and the objective function and adjusting the regular term and the objective function, the optimal candidate portal angle can be accurately screened out, and the optimal candidate portal angle is used as the target portal angle, which is beneficial to determining a reliable intensity modulation plan, and meanwhile, the number of portals can be effectively reduced, and the complexity of the radiotherapy plan is reduced.
Wherein the set threshold may be a minimum portal weight. If the hop count of the alternative portal angle exceeds the minimum portal weight, the alternative portal angle is used as the target portal angle, and if the hop count of the alternative portal angle is smaller than the minimum portal weight, the alternative portal angle is used as an invalid portal angle and is removed, so that the target portal angle is accurately screened out from the alternative portal angle.
Different from S230, this embodiment may also set a constraint condition of the maximum number of radiation fields, and if the constraint condition of the maximum number of radiation fields is set, the method for performing iterative optimization on the optimization model under each alternative radiation field angle includes: and carrying out iterative optimization on the optimization model based on the preset constraint condition of the maximum number of the radiation fields until the optimization model reaches a set optimization threshold value, and finishing iterative adjustment operation.
Specifically, the expression for the constraint is:
combining the constraint condition
Hop count and sigma in alternative field anglesj wijDynamically determining an optimization model g (d (x, b, w), d ═ 10)+α1∑i bi+α2∑ij wijDynamically adjusting the optimization model until the optimization model g (d (x, b, w), d0)+α1∑i bi+α2∑ij wijThe minimum value is obtained, and a minimum value optimization model is obtained: minx，b，w g(d(x，b，w)，d0)+α1∑ibi+α2∑ij wij. Wherein the content of the first and second substances,
is an exponential function of e, i.e.
E is the minimum portal weight.
It should be noted that the constraint condition may be understood as a linear constraint condition, and the constraint condition of the maximum number of fields in this embodiment may also be nonlinear. For example, the nonlinear constraint: sigmai sign(bi) Less than or equal to L, wherein,
in this embodiment, a constraint optimization method may be used to apply the nonlinear constraint condition: sigmaisign(bi) And (3) solving by less than or equal to L, and adopting an external penalty function method, an internal penalty function method and a Lagrange penalty function method to carry out nonlinear constraint conditions: sigmai sign(bi) And solving the solution with the L being less than or equal to the L.
Further, if the constraint condition of the maximum number of the radiation fields is set, the method for screening the target radiation field angle comprises the following steps: and taking the candidate portal angle with the maximum hop count corresponding to the maximum portal number as a target portal angle, and/or taking the candidate portal angle with the hop count exceeding a set threshold value as the target portal angle, and deleting the candidate portal angle with the hop count not exceeding the set threshold value. If the hop count of the alternative portal angle exceeds the minimum portal weight, taking the alternative portal angle as a target portal angle; if the hop count of the alternative portal angle is smaller than the minimum portal weight, taking the alternative portal angle as an invalid portal angle and removing the invalid portal angle; and/or taking the candidate portal angle with the maximum hop count corresponding to the maximum portal number as the target portal angle. Through setting up the constraint condition of the biggest radiation field number, can follow the target radiation field angle of selecting more accurately in the alternative radiation field angle, be favorable to reaching the purpose of accurate optimization accent plan.
S240, optimizing the shape and the hop count of the sub-field under the target field angle, and determining the optimized intensity modulation plan.
According to the technical scheme provided by the embodiment, a first regular term and a second regular term are calculated according to the sum of the hop count of each sub-field in any candidate portal and the ratio of the hop count of any sub-field in any candidate portal to the sum of the hop count, the regular term is determined according to the first regular term and the second regular term, the expression of an optimization model is further determined, the optimization model under each candidate portal angle is iteratively adjusted based on the leaf position in the initial intensity modulation plan, the hop count in the candidate portal angle and the weight of any sub-field in any candidate portal angle, and the sum of all sub-fields under any candidate portal angle is 1, or the optimization model under each candidate portal angle is iteratively optimized according to the constraint condition of the preset maximum portal number, and the candidate portal angle with the hop count exceeding a set threshold is used as a target portal angle, and deleting the alternative portal angle with the hop number not exceeding the set threshold, or taking the alternative portal angle with the hop number exceeding the set threshold as the target portal angle, and/or taking the alternative portal angle with the hop number exceeding the set threshold as the target portal angle, and deleting the alternative portal angle with the hop number not exceeding the set threshold, so that the optimal alternative portal angle can be accurately screened out, and the maximum alternative portal angle is taken as the target portal angle, which is favorable for achieving the purpose of accurately optimizing the intensity modulation plan.
EXAMPLE III
Fig. 4 is a schematic structural diagram of an emphasis plan optimization apparatus according to a third embodiment of the present invention, where the apparatus is configured in a processor of the emphasis plan optimization system. Referring to fig. 4, the apparatus includes: a subfield and hop count determination module 310, an optimization model determination module 320, a target portal angle screening module 330, and a emphasis plan optimization module 340.
The sub-field and hop count determining module 310 is configured to obtain a plurality of candidate portal angles, and determine the shape and hop count of the sub-field at each candidate portal angle;
an optimization model determining module 320, configured to determine an optimization model based on the shapes of the segments and the hop counts at each of the candidate portal angles, where the optimization model includes optimizing an objective function and a regular term, the regular term is determined according to the portal and/or the segments at each of the candidate portal angles, and the portal includes at least one of the segments;
the target radiation field angle screening module 330 is configured to perform iterative optimization on the optimization model, and screen a target radiation field angle from the alternative radiation field angles according to an iterative optimization result;
and an intensity modulation plan optimization module 340, configured to optimize the shape and the hop count of the sub-field at the target portal angle, and determine an optimized intensity modulation plan.
On the basis of the technical schemes, the alternative radiation field angles are coplanar, and the target radiation field angles are coplanar; alternatively, the first and second electrodes may be,
each alternative radiation field angle is not coplanar, and each target radiation field angle is coplanar; alternatively, the first and second electrodes may be,
each of the alternative field angles is non-coplanar and each of the target field angles is non-coplanar.
On the basis of the above technical solutions, the subfield and hop count determining module 320 is further configured to determine a flux optimization variable in the initial intensity modulation plan;
performing minimization calculation on a flux objective function determined according to the flux optimization variable at each alternative radiation field angle;
and carrying out blade serialization processing on the minimization calculation result of the flux objective function, and determining the sub-field and the hop count under the alternative radiation field angle.
On the basis of the above technical solutions, the sub-field and hop count determining module 320 is further configured to determine the shape and hop count of the sub-field under each candidate portal angle according to a direct machine parameter optimization method.
On the basis of the technical schemes, the sub-fields comprise static sub-fields and/or dynamic sub-fields.
On the basis of the technical schemes, the regular term is determined according to the hop count of the sub-fields.
On the basis of the above technical solutions, the regular term includes a first regular term and a second regular term, the first regular term is determined according to the hop count of the portal within the candidate portal angle, and the second regular term is determined according to the weight of the sub-portal within the candidate portal angle.
On the basis of the above technical solutions, the hop count of the portal within the candidate portal angle is determined according to the hop counts of all the sub-portals within the candidate portal angle, and the weight of the sub-portals is determined according to the hop count of any sub-portal within the candidate portal angle and the hop count of the portal within the candidate portal angle.
On the basis of the above technical solutions, the target portal angle screening module 330 is further configured to adjust the shape and the hop count of the sub-portal at each alternative portal angle, perform iterative optimization on the optimization model until the optimization model reaches a set optimization threshold, and end the iterative operation.
On the basis of the above technical solutions, the target portal angle screening module 330 is further configured to use the alternative portal angle with the hop count exceeding the set threshold as the target portal angle, and delete the alternative portal angle with the hop count not exceeding the set threshold.
On the basis of the above technical solutions, the target radiation field angle screening module 330 is further configured to perform iterative optimization on the optimization model based on a preset constraint condition of the maximum number of radiation fields until the optimization model reaches a set optimization threshold, and end the iterative adjustment operation.
On the basis of the foregoing technical solutions, the target portal angle screening module 330 is further configured to use the candidate portal angle with the largest hop count corresponding to the largest number of portals as the target portal angle, and/or use the candidate portal angle with the hop count exceeding a set threshold as the target portal angle, and delete the candidate portal angle with the hop count not exceeding the set threshold.
According to the technical scheme provided by the embodiment, a plurality of alternative field angles are obtained through a strength modulation plan optimization device in a strength design optimization system, the shape and the hop count of the sub-field under each alternative field angle are determined, an optimization model is determined based on the shape and the hop count of the sub-field under each alternative field angle, iterative optimization is performed on the optimization model, a target field angle is screened from the alternative field angles, the sub-field and the hop count under the target field angle are determined, the shape and the hop count of the sub-field under the target field angle are further optimized, and an optimized strength modulation plan is determined. Compared with the prior art that the field angle is not matched with the actual condition due to the fact that the field angle is set according to the experience of the user, the optimal field angle and the shape and the hop count of the sub-field under the optimal field angle can be automatically selected from a plurality of alternative field angles, the optimized intensity modulation plan is obtained, labor cost is reduced, the reliability of the intensity modulation plan is improved, meanwhile, the number of the sub-fields can be reduced as far as possible, and the complexity of the radiotherapy plan is reduced.
Example four
An embodiment of the present invention further provides a computer-readable storage medium, where the computer-executable instructions are configured by a computer processor to perform the following method:
acquiring a plurality of alternative field angles, and determining the shape and hop count of the sub-field under each alternative field angle;
determining an optimization model based on the shapes of the sub-fields and the hop counts under the alternative radiation field angles, wherein the optimization model comprises optimizing an objective function and a regular term, the regular term is determined according to the radiation field and/or the sub-fields under the alternative radiation field angles, and the radiation field comprises at least one sub-field;
performing iterative optimization on the optimization model, and screening a target field angle from the alternative field angles according to an iterative optimization result;
and optimizing the shape and hop count of the sub-field under the target field angle, and determining an optimized intensity modulation plan.
Of course, the computer program stored on the computer-readable storage medium provided by the embodiments of the present invention is not limited to the above method operations, and may also perform related operations in a method for enhancing plan optimization provided by any embodiments of the present invention.
Computer storage media for embodiments of the invention may employ any combination of one or more computer-readable media. The computer readable medium may be a computer readable signal medium or a computer readable storage medium. A computer readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, or device, or any combination of the foregoing. More specific examples (a non-exhaustive list) of the computer readable storage medium would include the following: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this document, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, or device.
The computer-readable signal medium may include, among other things, a candidate field angle, a shape and hop count of the sub-field at the candidate field angle, an objective function, a regularization term, a target field angle, and a shape and hop count of the sub-field at the target field angle, and may carry computer-readable program code. The alternative radiation field angle of the propagation, the shape and the hop count of the sub-field under the alternative radiation field angle, the objective function, the regular term, the target radiation field angle, the shape and the hop count of the sub-field under the target radiation field angle, and the like. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, or device.
Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
Computer program code for carrying out operations for aspects of the present invention may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C + +, or the like, as well as conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet service provider).
It should be noted that, in the embodiment of the emphasis plan optimization apparatus, the modules included in the apparatus are only divided according to functional logic, but are not limited to the above division as long as the corresponding functions can be implemented; 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.
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 described 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 (14)
1. An emphasis plan optimization system comprising: a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the storage device is for storing instructions, the processor is configured to perform the steps of:
acquiring a plurality of alternative field angles, and determining the shape and hop count of the sub-field under each alternative field angle;
determining an optimization model based on the shapes of the sub-fields and the hop counts under the alternative radiation field angles, wherein the optimization model comprises optimizing an objective function and a regular term, the regular term is determined according to the radiation field and/or the sub-fields under the alternative radiation field angles, and the radiation field comprises at least one sub-field;
performing iterative optimization on the optimization model, and screening a target field angle from the alternative field angles according to an iterative optimization result;
and optimizing the shape and hop count of the sub-field under the target field angle, and determining an optimized intensity modulation plan.
2. The system of claim 1, wherein each of the alternative field angles are coplanar and each of the target field angles are coplanar; alternatively, the first and second electrodes may be,
each alternative radiation field angle is not coplanar, and each target radiation field angle is coplanar; alternatively, the first and second electrodes may be,
each of the alternative field angles is non-coplanar and each of the target field angles is non-coplanar.
3. The system of claim 1, wherein said determining the shape and hop count of the segments at each of said candidate field angles comprises:
determining a flux optimization variable in the initial intensity modulation plan;
performing minimization calculation on a flux objective function determined according to the flux optimization variable at each alternative radiation field angle;
and carrying out blade serialization processing on the minimization calculation result of the flux objective function, and determining the sub-field and the hop count under the alternative radiation field angle.
4. The system of claim 1, wherein said determining the shape and hop count of the segments at each of said candidate field angles comprises:
and determining the shape and hop count of the sub-field under each alternative radiation field angle according to a direct machine parameter optimization method.
5. The system of claim 1, wherein the sub-fields comprise static sub-fields and/or dynamic sub-fields.
6. The system of claim 1, wherein the regularization term is determined according to a hop count of the sub-field.
7. The system of claim 1, wherein the regularization term comprises a first regularization term determined according to the hop count of the portal at the candidate portal angle and a second regularization term determined according to the weight of the segments within the candidate portal angle.
8. The system of claim 7, wherein the hop count of the portal under the alternative portal angle is determined according to the hop counts of all the sub-portals within the alternative portal angle, and the weights of the sub-portals are determined according to the hop count of any sub-portal within the alternative portal angle and the corresponding hop count of the portal within the alternative portal angle.
9. The system of claim 1, wherein the iterative optimization of the optimization model comprises:
and adjusting the shape and hop count of the sub-field under each alternative radiation field angle, performing iterative optimization on the optimization model until the optimization model reaches a set optimization threshold value, and ending the iterative operation.
10. The system of claim 1, wherein the screening of the candidate portal angles for a target portal angle according to the iterative optimization result comprises:
and taking the alternative portal angle with the hop count exceeding a set threshold value as the target portal angle, and deleting the alternative portal angle with the hop count not exceeding the set threshold value.
11. The system of claim 1, wherein the iterative optimization of the optimization model comprises:
and carrying out iterative optimization on the optimization model based on the preset constraint condition of the maximum number of the radiation fields until the optimization model reaches a set optimization threshold value, and finishing iterative adjustment operation.
12. The system of claim 11, wherein the screening of the candidate portal angles for a target portal angle according to the iterative optimization result comprises:
and taking the candidate portal angle with the maximum hop count corresponding to the maximum portal number as the target portal angle, and/or taking the candidate portal angle with the hop count exceeding a set threshold value as the target portal angle, and deleting the candidate portal angle with the hop count not exceeding the set threshold value.
13. An emphasis plan optimization device configured in a processor of the emphasis plan optimization system, the device comprising:
the device comprises a sub-field and hop count determining module, a shape determining module and a hop count determining module, wherein the sub-field and hop count determining module is used for acquiring a plurality of alternative radiation field angles and determining the shape and hop count of the sub-field under each alternative radiation field angle;
an optimization model determining module, configured to determine an optimization model based on the shapes of the segments and the hop counts at each of the candidate portal angles, where the optimization model includes optimizing an objective function and a regular term, the regular term is determined according to the portal and/or the segments at each of the candidate portal angles, and the portal includes at least one segment;
the target radiation field angle screening module is used for performing iterative optimization on the optimization model and screening a target radiation field angle from the alternative radiation field angles according to an iterative optimization result;
and the intensity modulation plan optimization module is used for optimizing the shape and the hop count of the sub-field under the target radiation field angle and determining an optimized intensity modulation plan.
14. A storage medium containing computer-executable instructions, wherein the computer-executable instructions are configured by a computer processor to perform a method comprising:
acquiring a plurality of alternative field angles, and determining the shape and hop count of the sub-field under each alternative field angle;
determining an optimization model based on the shapes of the sub-fields and the hop counts under the alternative radiation field angles, wherein the optimization model comprises optimizing an objective function and a regular term, the regular term is determined according to the radiation field and/or the sub-fields under the alternative radiation field angles, and the radiation field comprises at least one sub-field;
performing iterative optimization on the optimization model, and screening a target field angle from the alternative field angles according to an iterative optimization result;
and optimizing the shape and hop count of the sub-field under the target field angle, and determining an optimized intensity modulation plan.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010763340.5A CN111986778B (en) | 2020-07-31 | 2020-07-31 | Intensity-modulated plan optimization system, device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010763340.5A CN111986778B (en) | 2020-07-31 | 2020-07-31 | Intensity-modulated plan optimization system, device and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111986778A true CN111986778A (en) | 2020-11-24 |
| CN111986778B CN111986778B (en) | 2024-02-20 |
Family
ID=73445608
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010763340.5A Active CN111986778B (en) | 2020-07-31 | 2020-07-31 | Intensity-modulated plan optimization system, device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111986778B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113223664A (en) * | 2021-05-14 | 2021-08-06 | 山东第一医科大学附属肿瘤医院(山东省肿瘤防治研究院、山东省肿瘤医院) | Method, device and medium for calculating complexity of volume rotation intensity modulation plan |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103083821A (en) * | 2013-01-27 | 2013-05-08 | 合肥超安医疗科技有限公司 | System and a method of directly optimizing static-state intensity-modulated sub-field shapes and machine hop count thereof |
| CN103706044A (en) * | 2014-01-01 | 2014-04-09 | 杨金磊 | Design method for radiotherapeutic plan |
| CN105031820A (en) * | 2015-09-01 | 2015-11-11 | 上海联影医疗科技有限公司 | Intensity modulated radiation therapy reverse optimization method and device |
| CN105617536A (en) * | 2015-12-24 | 2016-06-01 | 上海联影医疗科技有限公司 | Inverse intensity modulated arc radioterapy optimization method and device |
| CN107469239A (en) * | 2017-08-11 | 2017-12-15 | 上海联影医疗科技有限公司 | Computer-readable recording medium and radiotherapy treatment planning system |
| 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 |
| CN107998519A (en) * | 2017-12-26 | 2018-05-08 | 荣科科技股份有限公司 | A kind of molecular dynamics re-optimization algorithm for IMRT |
| US20180243586A1 (en) * | 2017-02-28 | 2018-08-30 | Sun Nuclear Corporation | Radiation therapy treatment verification with electronic portal imaging device transit images |
| CN110211664A (en) * | 2019-04-25 | 2019-09-06 | 安徽大学 | A kind of system based on machine learning Automated Design radiation treatment plan |
| CN110246562A (en) * | 2018-03-07 | 2019-09-17 | 有方（合肥）医疗科技有限公司 | Determine the method, apparatus and computer system of the sub-beam intensity in radiotherapy system |
| CN110289075A (en) * | 2019-04-29 | 2019-09-27 | 中北大学 | A kind of direct Ziye optimization method and system based on fuzzy entropy |
| US20200101320A1 (en) * | 2018-09-28 | 2020-04-02 | Varian Medical Systems International Ag | Beam angle optimization for external beam radiation therapy using sectioning |
-
2020
- 2020-07-31 CN CN202010763340.5A patent/CN111986778B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103083821A (en) * | 2013-01-27 | 2013-05-08 | 合肥超安医疗科技有限公司 | System and a method of directly optimizing static-state intensity-modulated sub-field shapes and machine hop count thereof |
| CN103706044A (en) * | 2014-01-01 | 2014-04-09 | 杨金磊 | Design method for radiotherapeutic plan |
| CN105031820A (en) * | 2015-09-01 | 2015-11-11 | 上海联影医疗科技有限公司 | Intensity modulated radiation therapy reverse optimization method and device |
| CN105617536A (en) * | 2015-12-24 | 2016-06-01 | 上海联影医疗科技有限公司 | Inverse intensity modulated arc radioterapy optimization method and device |
| 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 |
| US20180243586A1 (en) * | 2017-02-28 | 2018-08-30 | Sun Nuclear Corporation | Radiation therapy treatment verification with electronic portal imaging device transit images |
| CN107469239A (en) * | 2017-08-11 | 2017-12-15 | 上海联影医疗科技有限公司 | Computer-readable recording medium and radiotherapy treatment planning system |
| CN107998519A (en) * | 2017-12-26 | 2018-05-08 | 荣科科技股份有限公司 | A kind of molecular dynamics re-optimization algorithm for IMRT |
| CN110246562A (en) * | 2018-03-07 | 2019-09-17 | 有方（合肥）医疗科技有限公司 | Determine the method, apparatus and computer system of the sub-beam intensity in radiotherapy system |
| US20200101320A1 (en) * | 2018-09-28 | 2020-04-02 | Varian Medical Systems International Ag | Beam angle optimization for external beam radiation therapy using sectioning |
| CN110211664A (en) * | 2019-04-25 | 2019-09-06 | 安徽大学 | A kind of system based on machine learning Automated Design radiation treatment plan |
| CN110289075A (en) * | 2019-04-29 | 2019-09-27 | 中北大学 | A kind of direct Ziye optimization method and system based on fuzzy entropy |
Non-Patent Citations (3)
| Title |
|---|
| AARON BABIER等: "Knowledge-based automated planning with three-dimensional generative adversarial networks", 《ARXIV》, pages 1 - 15 * |
| D. BERTSIMAS等: "A hybridapproachtobeamangleoptimizationinintensity-modulated radiation therapy", 《COMPUTERS &OPERATIONSRESEARCH》, pages 2187 - 2197 * |
| 石翔翔等: "宫颈癌术后调强计划的最优子野数目选择", 《中国医学装备》, vol. 16, no. 1, pages 32 - 36 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113223664A (en) * | 2021-05-14 | 2021-08-06 | 山东第一医科大学附属肿瘤医院(山东省肿瘤防治研究院、山东省肿瘤医院) | Method, device and medium for calculating complexity of volume rotation intensity modulation plan |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111986778B (en) | 2024-02-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109621228B (en) | Radiation dose calculation device, radiation dose calculation apparatus, and storage medium | |
| US11103727B2 (en) | Model based PBS optimization for flash therapy treatment planning and oncology information system | |
| US20230302296A1 (en) | Method, apparatus, and system for simulating a particle transport and determining human dose in a radiotherapy | |
| EP3388108B1 (en) | Systems and methods for dose calculation in generating radiation treatment plans | |
| CN110420396B (en) | Radiotherapy plan optimization system, radiotherapy plan optimization device and storage medium | |
| US20180178036A1 (en) | Selecting a dose prediction model based on clinical goals | |
| US20190275352A1 (en) | Treatment planning based on multiple modalities | |
| US10512790B2 (en) | Systems and methods for generating radiation treatment plans | |
| CN111001097B (en) | Radiotherapy dose evaluation system, device and storage medium | |
| EP3302699A1 (en) | Method of selecting beam geometries | |
| Yue et al. | Dose prediction via distance-guided deep learning: Initial development for nasopharyngeal carcinoma radiotherapy | |
| CN106408543A (en) | Blocking optical grating optimization method and device for scattering correction of cone-beam CT (computed tomography) image | |
| US10881875B2 (en) | Pre-optimization method for quick prediction of achievability of clinical goals in intensity modulated radiation therapy | |
| CN111986778A (en) | Intensity modulated plan optimization system, device and storage medium | |
| EP4252840A1 (en) | Monitor unit optimization constraint in radiation treatment planning | |
| US11173322B2 (en) | Systems and methods for treatment planning using continuous quality metrics | |
| US10737114B2 (en) | Translating different clinical protocols for particle therapy into a set of constraints | |
| CN115966281B (en) | Method, device, equipment and storage medium for generating arc radiotherapy plan | |
| WO2020047537A1 (en) | Dosevolume histogram and dose distribution based autoplanning | |
| US20230390579A1 (en) | Virtual Particle Based Monte Carlo Dose Calculation for Charged Particle Therapy Treatment Planning | |
| US9987501B1 (en) | Extracting protobeams for cancer radiation therapy | |
| CN117138260A (en) | Treatment planning system, method, equipment and medium for arc treatment of tumor by ion beam | |
| CN118160045A (en) | Method, device, equipment and storage medium for generating arc radiotherapy plan | |
| WO2021259690A1 (en) | Dynamic estimation of a biological effect of a variable composition of non-photon radiation |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |