CN105617536A - Inverse intensity modulated arc radioterapy optimization method and device - Google Patents

Inverse intensity modulated arc radioterapy optimization method and device Download PDF

Info

Publication number
CN105617536A
CN105617536A CN201510990902.9A CN201510990902A CN105617536A CN 105617536 A CN105617536 A CN 105617536A CN 201510990902 A CN201510990902 A CN 201510990902A CN 105617536 A CN105617536 A CN 105617536A
Authority
CN
China
Prior art keywords
launched field
ziye
alternative
radiation therapy
intensity modulated
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.)
Pending
Application number
CN201510990902.9A
Other languages
Chinese (zh)
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 CN201510990902.9A priority Critical patent/CN105617536A/en
Publication of CN105617536A publication Critical patent/CN105617536A/en
Pending legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

The invention provides an inverse intensity modulated arc radioterapy optimization method and device. The method comprises a step of obtaining the alternative radiation fields of multiple directions, dividing an organ into multiple organ vocals, dividing the alternative radiation field of each direction into multiple sub beams, and calculating the unit dose contribution of each sub beam of the alternative radiation field of each direction, a step of iterating sequentially and taking one or more alternative radiation fields which allows the planning quality of the intensity modulated radioterapy to be improved to the maximum in the alternative radiation fields of multiple directions as reference radiation fields, and determining the segment shape corresponding to each reference radiation field and the needed intensity, a step of optimizing the segment corresponding to each reference radiation field, a step of outputting each reference radiation field which satisfies a setting requirement and a corresponding segment shape and the intensity needed by the corresponding segment shape. According to the device and the device, a large-scale nonlinear problem is divided into two smaller problems to be solved separately, in each time of iteration, the global optimal solutions of the two smaller problems can be obtained in a polynomial time, the algorithm converges fast, and the workload of a user is reduced.

Description

Rotate reverse intensity modulated radiation therapy optimization method and device
Technical field
The invention mainly relates to radiotherapy technology field, particularly relate to the reverse intensity modulated radiation therapy optimization method of a kind of rotation and device.
Background technology
Radiotherapy is the treatment utilizing one or more ionizing radiation that malignant tumor and some benigns are carried out. Intensity modulated radiation therapy (intensitymodulatedradiationtherapy, IMRT) is namely the one of radiotherapy, three dimensional conformal radiation therapy develop. Intensity modulated radiation therapy requires that radiation field size (hereinafter launched field) interior intensity of beam is adjusted by certain requirement. When intensity modulated radiation therapy is that launched field is consistent with target area profile throughout, intensity of beam is adjusted by the concrete anatomy relationship for target area 3D shape and critical organ and target area, in single launched field, intensity of beam is unevenly distributed, but in whole Gross Target Volume, dose distribution is than three dimensional conformal radiation therapy evenly.
Intensity modulated radiation therapy is divided into multiple tiny Ziye each launched field. When formulating radiotherapy planning, according to target area 3D shape and to the relevant anatomy relationship jeopardized between organ, these Ziyes are distributed with different weights, make to produce in same launched field that optimize, uneven intensity distributions, to make the beam flux by jeopardizing organ reduce, and the beam flux of other parts of target area increases.
Inverse algorithm is adopted to formulate intensity modulated radiation therapy plan at present clinically, namely specified the target dose in the direction of launched field and number, tumor and the dose constraint jeopardizing on organ by user, optimized algorithm iteration provide Ziye number, shape and the weight that each launched field comprises. In the reverse intensity modulated therapy mode of rotation that United States Patent (USP) 12/834067 and 8229071B2 propose, in the continuous rotary course of radiotherapy head discontinuous go out bundle, generally two go out to restraint to exist between interval bundle interval, a pass, and Ziye shape remains unchanged when going out and restrainting, the consecutive variations when closing bundle, but, this patent does not provide concrete optimized algorithm. In addition, Kainz etc. propose to adopt simulated annealing to be optimized, but it is very slow that simulated annealing solves Large-scale Optimization Problems speed, and optimum results depends critically upon parameter and selects, referring specifically to open source literature " AplanninganddeliverystudyofarotationalIMRTtechniquewithb urstdelivery ".
Summary of the invention
It is an object of the invention to, it is provided that the reverse intensity modulated radiation therapy optimization method of a kind of rotation and device, it is possible to optimize launched field direction, solve in prior art computationally intensive, the problem that optimal speed is slow.
For solving above-mentioned technical problem, the present invention provides rotation reverse intensity modulated radiation therapy optimization method, including:
Obtaining the alternative launched field of multiple directions, organ is divided into multiple organ volume elements, the alternative launched field of all directions is divided into many sub-beams, the unit dose of each organ volume elements is contributed by each sub-beam of the alternative launched field calculating all directions;
The plan quality making intensity modulated radiation therapy in the alternative launched field of multiple directions is promoted maximum one or more alternative launched field as reference launched field by progressive alternate successively, and determines each reference launched field corresponding Ziye shape and required intensity thereof;
It is optimized with reference to the Ziye corresponding to launched field each;
Export to meet and set each reference launched field and corresponding Ziye shape thereof required, and the intensity needed for corresponding Ziye shape.
Optionally, trying to achieve with reference to launched field according to the first model and determine corresponding Ziye shape, choosing object function isWherein, wjFor the weight of jth organ volume elements, DijFor i-th sub-beam, the unit dose of jth organ volume elements is contributed, ykFor the dose contribution that the alternative launched field of kth is current, tjFor jth organ volume elements target dose, object function is to ykPartial derivative beDjAll with reference to the launched field actual doses to jth organ volume elements for current optimization, draw gkThe set being worth one group of minimum sub-beam composition is with reference to launched field.
Optionally, when carrying out optimizing for the first time, DjEqual to zero.
Optionally, according to the intensity needed for this reference launched field of one second model solution, according toObtain Dkj, DkjFor the kth of this current optimization with reference to the launched field dose contribution to jth organ volume elements, according to f (Yk)=�� wj(��YkDkj-tj)2To ykIt is optimized.
Optionally, number or number scope with reference to launched field are scheduled.
Optionally, the Ziye intensity of generation has bound to retrain, and Ziye shape must is fulfilled for machine constraint and the other influences constraint by dosage verifying.
Optionally, the step each each Ziye comprised with reference to launched field being optimized includes: carry out Ziye shape and weight optimization.
Accordingly, the present invention also provides for the reverse intensity modulated radiation therapy of a kind of rotation and optimizes device, including:
Launched field angle and Shape optimization module, it is configured to obtain the alternative launched field of multiple directions, organ is divided into multiple organ volume elements, the alternative launched field of all directions is divided into many sub-beams, the unit dose of each organ volume elements is contributed by each sub-beam of the alternative launched field calculating all directions, the plan quality making intensity modulated radiation therapy in the alternative launched field of multiple directions is promoted maximum one or more alternative launched field as reference launched field, and determines each reference launched field corresponding Ziye shape and required dosage thereof;
Direct Ziye optimizes module, connects this launched field angle and Shape optimization module, is configured to each each Ziye comprised with reference to launched field is optimized, and output meets the Ziye optimum results that setting requires.
Optionally, trying to achieve with reference to launched field according to the first model and determine corresponding Ziye shape, choosing object function isWherein, wjFor the weight of jth organ volume elements, DijFor i-th sub-beam, the unit dose of jth organ volume elements is contributed, ykFor the dose contribution that the alternative launched field of kth is current, tjFor jth organ volume elements target dose, object function is to ykPartial derivative beDjAll with reference to the launched field actual doses to jth organ volume elements for current optimization, draw gkThe set being worth one group of minimum sub-beam composition is with reference to launched field.
Optionally, when carrying out optimizing for the first time, DjEqual to zero.
Optionally, according to the intensity needed for this reference launched field of one second model solution, according toObtain Dkj, DkjFor the kth of this current optimization with reference to the launched field dose contribution to jth organ volume elements, according to f (Yk)=�� wj(��YkDkj-tj)2To ykIt is optimized.
Optionally, also include dose calculation module, connect this launched field angle and Shape optimization module and this direct Ziye optimizes module, be configured to calculate the dose distribution of launched field.
Optionally, each each Ziye comprised with reference to launched field is optimized by this direct Ziye optimization module by carrying out Ziye shape and weight adjustment.
Optionally, also include dose volume histogram computing module, connect this launched field angle and Shape optimization module and this direct Ziye optimizes module, be configured to calculate the dose volume histogram of launched field.
Compared with prior art, in the present invention, large-scale nonlinear PROBLEM DECOMPOSITION becomes two small-scale problems solve respectively, each iteration can respectively obtain the globally optimal solution of two small-scale problems in polynomial time, algorithmic statement is fast, reduces amount of user effort.
Accompanying drawing explanation
Fig. 1 is the flow chart rotating reverse intensity modulated radiation therapy optimization method in one embodiment of the invention;
Fig. 2 is the schematic diagram of multiple directions launched field in one embodiment of the invention;
Fig. 3 is the structural representation that organ is divided in one embodiment of the invention multiple volume elements;
Fig. 4 rotates reverse intensity modulated radiation therapy to optimize the flow chart of device in one embodiment of the invention.
Detailed description of the invention
It is described in more detail below in conjunction with the schematic diagram reverse intensity modulated radiation therapy optimization method of rotation to the present invention and device, which show the preferred embodiments of the present invention, should be appreciated that those skilled in the art can revise invention described herein, and still realize the advantageous effects of the present invention. Therefore, it is widely known that description below is appreciated that for those skilled in the art, and is not intended as limitation of the present invention.
The rotation reverse intensity modulated radiation therapy optimization method of the present invention includes: obtain the alternative launched field of multiple directions, organ is divided into multiple organ volume elements, the alternative launched field of all directions is divided into many sub-beams, and the unit dose of each organ volume elements is contributed by each sub-beam of the alternative launched field calculating all directions; The plan quality making intensity modulated radiation therapy in the alternative launched field of multiple directions is promoted maximum one or more alternative launched field as reference launched field by progressive alternate successively, and determines each reference launched field corresponding Ziye shape and required intensity thereof; It is optimized with reference to the Ziye corresponding to launched field each; Export to meet and set each reference launched field and corresponding Ziye shape thereof required, and the intensity needed for corresponding Ziye shape. In the present invention, large-scale nonlinear PROBLEM DECOMPOSITION becoming two small-scale problems solve respectively, can respectively obtain the globally optimal solution of two small-scale problems in each iteration in polynomial time, algorithmic statement is fast, reduces amount of user effort.
Being described in detail below in conjunction with the accompanying drawing 1��Fig. 4 reverse intensity modulated radiation therapy optimization method of rotation to the present invention and device, Fig. 1 is the flow chart rotating reverse intensity modulated radiation therapy optimization method, specifically includes following steps:
Perform step S1, with reference to, shown in Fig. 2, obtaining the alternative launched field of multiple directions, op1 in Fig. 2, op2, op3 ... wait direction, in the present embodiment, all directions of alternative launched field can be that user is set in advance, it is also possible to is drawn by the optimization method of the present invention. In Patients During Radiotherapy, at interval [CP0, CP1], [CP2, CP3], [CP4, CP5], [CP6, CP7] close bundle, namely do not carry out radiotherapy, and in closing bundle and be interval from once open bundle time Ziye shape vary continuously to once open bundle time required Ziye shape. And at the interval [-a of residue1,+a1], [-a2,+a2], [-a3,+a3] open bundle, in opening bundle and be interval, keep Ziye shape invariance and go out Shu Jinhang radiotherapy. In the present embodiment, it is necessary to find out the suitable launched field shape in all directions and intensity.
With reference to shown in Fig. 3, organ is separated into the set of multiple organ volume elements, such as, organ is divided into j volume elements, simultaneously, the alternative launched field of all directions is divided into many sub-beams, such as, the alternative launched field of all directions includes i sub-beam altogether, and the alternative launched field in each direction includes parton beam therein, afterwards, the unit dose of each organ volume elements is contributed D by each sub-beam of the alternative launched field calculating all directionsij, DijBeing the i-th sub-beam unit dose to jth organ volume elements to contribute, unit is cGy/MU. In an embodiment of the present invention, the method such as pencil beam (PencilBeam), convolution (Convolution) or Monte Carlo (MonteCarlo) can be used to calculate and to obtain unit dose contribution Dij��
Performing step S2, the plan quality making intensity modulated radiation therapy in the alternative launched field of multiple directions is promoted maximum one or more alternative launched field as reference launched field by progressive alternate successively, and determines each reference launched field corresponding Ziye shape and required intensity thereof.
Concrete, to try to achieve with reference to launched field according to the first model and determine corresponding Ziye shape, the first model is chosen object function and is
Wherein, wjFor the weight of jth organ volume elements, weight wjIt is set previously according to needs for user, DijFor i-th sub-beam, the unit dose of jth organ volume elements is contributed, ykFor the dose contribution that the alternative launched field of kth is current, tjFor jth organ volume elements target dose, target dose tjIt is set previously according to needs for user. Then, by object function to ykSeek local derviation, then object function is to ykPartial derivative be
g k = ∂ f ∂ y k = 2 Σ j w j { ( Σ Y k Σ i ∈ k D i j - t j ) Σ i ∈ k D i j } = 2 Σ j { w j ( D j - t j ) Σ i ∈ k D i j } = 2 Σ i ∈ k { Σ j w j ( D j - t j ) D i j }
Wherein, DjAll reference launched field actual doses to jth organ volume elements for current optimization, it is to be understood that when carrying out optimizing for the first time, DjEqual to zero, when carrying out optimizing for the second time, DjEqual to the actual dose of the reference launched field of first time optimization, when carrying out optimizing for the third time, DjEqual to the sum of first time and the actual dose of the reference launched field of second time optimization, and so on. In the present embodiment, by calculating gkThe set being worth one group of minimum sub-beam composition is with reference to launched field, say, that each with reference to the corresponding Ziye shape of launched field. With reference to the number of launched field less than or equal to the sum of alternative launched field, and, number or the number scope of the reference launched field finally chosen are scheduled, for instance can be preset according to tumor type by algorithm, or be specified by user.
Then, according to the intensity needed for this reference launched field of one second model solution, according toObtain Dkj, DkjFor the kth of this current optimization with reference to the launched field dose contribution to jth organ volume elements, further according to object function f (Yk)=�� wj(��YkDkj-tj)2To ykIt is optimized, wherein, Yk �� [MUmin, MUmax]. In the present embodiment, the Ziye intensity of generation has bound to retrain, and Ziye shape must is fulfilled for machine constraint and the other influences constraint by dosage verifying, it is to avoid produce the Ziye of poor efficiency, extended treatment time.
Afterwards, carry out progressive alternate again and the plan quality making intensity modulated radiation therapy in the alternative launched field of the multiple directions maximum alternative launched field of lifting is selected as reference launched field successively, form the set with reference to launched field, and determine that each is with reference to launched field corresponding Ziye shape and required intensity thereof. In the present invention, if the number of required reference launched field is less than alternative launched field total, the optimization method of the present invention can select from multiple alternative launched fields preferably launched field as reference launched field.
In the present invention, large-scale nonlinear PROBLEM DECOMPOSITION become two small-scale problems solve respectively, namely first passing through the first model solution and go out the shape with reference to launched field, go out the intensity with reference to launched field again through the second model solution, progressive alternate solves out with reference to launched field by all of. Further, can respectively obtaining the globally optimal solution of two small-scale problems in each iteration in polynomial time, algorithmic statement is fast, reduces amount of user effort.
Performing step S3, be optimized with reference to the Ziye corresponding to launched field each, the step that each each Ziye comprised with reference to launched field is optimized includes: carry out Ziye shape and weight optimization. For example, following methods can be adopted to adjust Ziye shape: the blade to each multi-leaf optical grating of each Ziye, randomly choose blade moving direction and step-length; Calculating the dose distribution after blade moves and the value of corresponding object function, if target function value declines, then accept this blade and move, otherwise refusal blade moves.
Perform step S4, export and meet each reference launched field and corresponding Ziye shape thereof setting requirement, and the intensity needed for corresponding Ziye shape. In the present embodiment, it is necessary to judge whether Ziye optimum results meets the requirement that user sets, if it is export result. If it is not, illustrate that this case is more complicated, it is possible to the calculating repeating step 304 supplements new alternative launched field, and repeats step S2��S3. Certainly, if there being bigger assurance expection Ziye optimum results to disclosure satisfy that setting requirement, then extra judgement need not be carried out, directly output Ziye optimum results.
Accordingly, the present invention also provides for the reverse intensity modulated radiation therapy of a kind of rotation and optimizes device, including:
Launched field angle and Shape optimization module 10, launched field angle and Shape optimization module 10 are configured to obtain the alternative launched field of multiple directions, organ is divided into multiple organ volume elements, the alternative launched field of all directions is divided into many sub-beams, the unit dose of each organ volume elements is contributed by each sub-beam of the alternative launched field calculating all directions, the plan quality making intensity modulated radiation therapy in the alternative launched field of multiple directions is promoted maximum one or more alternative launched field as reference launched field, and determines each reference launched field corresponding Ziye shape and required intensity thereof. In the present invention, launched field angle and Shape optimization module 10 are tried to achieve with reference to launched field according to the first model and determine corresponding Ziye shape, and choosing object function isWherein, wjFor the weight of jth organ volume elements, DijFor i-th sub-beam, the unit dose of jth organ volume elements is contributed, ykFor the dose contribution that the alternative launched field of kth is current, tjFor jth organ volume elements target dose, object function is to ykPartial derivative be g k = ∂ f ∂ y k = 2 Σ j w j { ( Σ Y k Σ i ∈ k D i j - t j ) Σ i ∈ k D i j } = 2 Σ j { w j ( D j - t j ) Σ i ∈ k D i j } = 2 Σ i ∈ k { Σ j w j ( D j - t j ) D i j } , DjAll with reference to the launched field actual doses to jth organ volume elements for current optimization, draw gkThe set being worth one group of minimum sub-beam composition is with reference to launched field. It should be noted that when carrying out optimizing for the first time, DjEqual to zero. Then, launched field angle and Shape optimization module 10 according to the intensity needed for this reference launched field of one second model solution, according toObtain Dkj, DkjFor the kth of this current optimization with reference to the launched field dose contribution to jth organ volume elements, according to f (Yk)=�� wj(��YkDkj-tj)2To ykIt is optimized. In the present embodiment, the Ziye intensity of generation has bound to retrain, and Ziye shape must is fulfilled for machine constraint and the other influences constraint by dosage verifying. Number or number scope with reference to launched field are scheduled.
Direct Ziye optimizes module 20, connects this launched field angle and Shape optimization module 10, is configured to each each Ziye comprised with reference to launched field is optimized, and output meets the Ziye optimum results that setting requires. The step that each each Ziye comprised with reference to launched field is optimized includes: carry out Ziye shape and weight optimization.
Data input/output module 30, connects launched field angle and Shape optimization module 10, is configured to accept the setting of maximum Ziye number under the bound constraint of each Ziye intensity, minimum Ziye area and each launched field. It is appreciated that data input/output module 30 can also accept out of Memory when needed, for instance obtain patient data, machine parameter and other user setup parameter.
Dose calculation module 40, connects this launched field angle and Shape optimization module 10 and this direct Ziye optimizes module 20, is configured to calculate the dose distribution of launched field. Launched field angle and Shape optimization module 10 and direct Ziye optimize module 20 can provide, to dose calculation module 40, the launched field being currently needed for optimizing, the latter return present dose distribution.
Dose volume histogram computing module 50, connects this launched field angle and Shape optimization module 10 and this direct Ziye optimizes module 20, is configured to calculate the dose volume histogram of launched field. Launched field angle and Shape optimization module 10 and direct Ziye optimize module 20 can provide present dose to be distributed to dose volume histogram computing module 50, the latter return dose volume histogram.
In sum, in the present invention, large-scale nonlinear PROBLEM DECOMPOSITION becoming two small-scale problems solve respectively, can respectively obtain the globally optimal solution of two small-scale problems in each iteration in polynomial time, algorithmic statement is fast, reduces amount of user effort.
Obviously, the present invention can be carried out various change and modification without deviating from the spirit and scope of the present invention by those skilled in the art. So, if these amendments of the present invention and modification belong within the scope of the claims in the present invention and equivalent technologies thereof, then the present invention is also intended to comprise these change and modification.

Claims (14)

1. one kind rotates reverse intensity modulated radiation therapy optimization method, it is characterised in that including:
Obtaining the alternative launched field of multiple directions, organ is divided into multiple organ volume elements, the alternative launched field of all directions is divided into many sub-beams, the unit dose of each organ volume elements is contributed by each sub-beam of the alternative launched field calculating all directions;
The plan quality making intensity modulated radiation therapy in the alternative launched field of multiple directions is promoted maximum one or more alternative launched field as reference launched field by progressive alternate successively, and determines each reference launched field corresponding Ziye shape and required intensity thereof;
It is optimized with reference to the Ziye corresponding to launched field each;
Export to meet and set each reference launched field and corresponding Ziye shape thereof required, and the intensity needed for corresponding Ziye shape.
2. the reverse intensity modulated radiation therapy optimization method of rotation as claimed in claim 1, it is characterised in that trying to achieve with reference to launched field according to the first model and determine corresponding Ziye shape, choosing object function isWherein, wjFor the weight of jth organ volume elements, DijFor i-th sub-beam, the unit dose of jth organ volume elements is contributed, ykFor the dose contribution that the alternative launched field of kth is current, tjFor jth organ volume elements target dose, object function is to ykPartial derivative beDjAll with reference to the launched field actual doses to jth organ volume elements for current optimization, draw gkThe set being worth one group of minimum sub-beam composition is with reference to launched field.
3. the reverse intensity modulated radiation therapy optimization method of rotation as claimed in claim 2, it is characterised in that when carrying out optimizing for the first time, DjEqual to zero.
4. the reverse intensity modulated radiation therapy optimization method of rotation as claimed in claim 2, it is characterised in that according to the intensity needed for this reference launched field of one second model solution, according toObtain Dkj, DkjFor the kth of this current optimization with reference to the launched field dose contribution to jth organ volume elements, according to f (Yk)=�� wj(��YkDkj-tj)2To ykIt is optimized.
5. the reverse intensity modulated radiation therapy optimization method of rotation as claimed in claim 1, it is characterised in that number or number scope with reference to launched field are scheduled.
6. the reverse intensity modulated radiation therapy optimization method of rotation as claimed in claim 1, it is characterised in that the Ziye intensity of generation has bound to retrain, Ziye shape must is fulfilled for machine constraint and the other influences constraint by dosage verifying.
7. the reverse intensity modulated radiation therapy optimization method of rotation as claimed in claim 1, it is characterised in that the step that each each Ziye comprised with reference to launched field is optimized includes: carry out Ziye shape and weight optimization.
8. one kind rotates reverse intensity modulated radiation therapy optimization device, it is characterised in that including:
Launched field angle and Shape optimization module, it is configured to obtain the alternative launched field of multiple directions, organ is divided into multiple organ volume elements, the alternative launched field of all directions is divided into many sub-beams, the unit dose of each organ volume elements is contributed by each sub-beam of the alternative launched field calculating all directions, the plan quality making intensity modulated radiation therapy in the alternative launched field of multiple directions is promoted maximum one or more alternative launched field as reference launched field, and determines each reference launched field corresponding Ziye shape and required dosage thereof;
Direct Ziye optimizes module, connects this launched field angle and Shape optimization module, is configured to each each Ziye comprised with reference to launched field is optimized, and output meets the Ziye optimum results that setting requires.
9. the reverse intensity modulated radiation therapy of rotation as claimed in claim 8 optimizes device, it is characterised in that trying to achieve with reference to launched field according to the first model and determine corresponding Ziye shape, choosing object function isWherein, wjFor the weight of jth organ volume elements, DijFor i-th sub-beam, the unit dose of jth organ volume elements is contributed, ykFor the dose contribution that the alternative launched field of kth is current, tjFor jth organ volume elements target dose, object function is to ykPartial derivative beDjAll with reference to the launched field actual doses to jth organ volume elements for current optimization, draw gkThe set being worth one group of minimum sub-beam composition is with reference to launched field.
10. the reverse intensity modulated radiation therapy of rotation as claimed in claim 9 optimizes device, it is characterised in that when carrying out optimizing for the first time, DjEqual to zero.
11. rotate reverse intensity modulated radiation therapy as claimed in claim 9 to optimize device, it is characterised in that according to the intensity needed for this reference launched field of one second model solution, according toObtain Dkj, DkjFor the kth of this current optimization with reference to the launched field dose contribution to jth organ volume elements, according to f (Yk)=�� wj(��YkDkj-tj)2To ykIt is optimized.
12. rotate reverse intensity modulated radiation therapy as claimed in claim 8 to optimize device, it is characterised in that also include dose calculation module, connect this launched field angle and Shape optimization module and this direct Ziye optimizes module, be configured to calculate the dose distribution of launched field.
13. rotate reverse intensity modulated radiation therapy as claimed in claim 8 to optimize device, it is characterised in that this direct Ziye optimizes module and by carrying out Ziye shape and weight adjustment, each each Ziye comprised with reference to launched field is optimized.
14. rotate reverse intensity modulated radiation therapy as claimed in claim 8 to optimize device, it is characterized in that, also include dose volume histogram computing module, connect this launched field angle and Shape optimization module and this direct Ziye optimizes module, be configured to calculate the dose volume histogram of launched field.
CN201510990902.9A 2015-12-24 2015-12-24 Inverse intensity modulated arc radioterapy optimization method and device Pending CN105617536A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510990902.9A CN105617536A (en) 2015-12-24 2015-12-24 Inverse intensity modulated arc radioterapy optimization method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510990902.9A CN105617536A (en) 2015-12-24 2015-12-24 Inverse intensity modulated arc radioterapy optimization method and device

Publications (1)

Publication Number Publication Date
CN105617536A true CN105617536A (en) 2016-06-01

Family

ID=56033121

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510990902.9A Pending CN105617536A (en) 2015-12-24 2015-12-24 Inverse intensity modulated arc radioterapy optimization method and device

Country Status (1)

Country Link
CN (1) CN105617536A (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106902478A (en) * 2017-02-22 2017-06-30 天津大学 A kind of systematization radiotherapy biological effect appraisal procedure
CN107349531A (en) * 2017-07-04 2017-11-17 上海联影医疗科技有限公司 A kind of radiotherapy apparatus, storage medium, Ziye shape generation method and device
CN110246562A (en) * 2018-03-07 2019-09-17 有方（合肥）医疗科技有限公司 Determine the method, apparatus and computer system of the sub-beam intensity in radiotherapy system
CN110582325A (en) * 2017-05-03 2019-12-17 光线搜索实验室公司 System and method for ion-based radiotherapy treatment planning
CN111437521A (en) * 2020-04-23 2020-07-24 华东医院 Non-uniform volume arc intensity modulation method
CN111785378A (en) * 2020-06-19 2020-10-16 新里程医用加速器(无锡)有限公司 Radiotherapy planning system and storage medium
CN111986778A (en) * 2020-07-31 2020-11-24 上海联影医疗科技股份有限公司 Intensity modulated plan optimization system, device and storage medium
CN112043974A (en) * 2020-09-08 2020-12-08 苏州雷泰医疗科技有限公司 Dynamic intensity modulation method and device based on orthogonal double-layer grating rotary sweep
CN112972912A (en) * 2021-04-22 2021-06-18 北京大学第三医院（北京大学第三临床医学院） Non-coplanar radiotherapy beam incident path optimization method
CN113521561A (en) * 2021-06-22 2021-10-22 苏州雷泰医疗科技有限公司 Method, storage medium and device for realizing reverse rotation intensity modulation of orthogonal double-layer grating
CN115966281A (en) * 2023-01-16 2023-04-14 迈胜医疗设备有限公司 Method, device, equipment and storage medium for generating arc radiotherapy plan

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11142600A (en) * 1997-11-12 1999-05-28 Mitsubishi Electric Corp Charged particle beam irradiation device and irradiation method
CN101422640A (en) * 2008-11-25 2009-05-06 中国科学院等离子体物理研究所 Multiple-objective optimization method and system capable of optimizing radiotherapy beam intensity distribution
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
CN105031820A (en) * 2015-09-01 2015-11-11 上海联影医疗科技有限公司 Intensity modulated radiation therapy reverse optimization method and device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11142600A (en) * 1997-11-12 1999-05-28 Mitsubishi Electric Corp Charged particle beam irradiation device and irradiation method
CN101422640A (en) * 2008-11-25 2009-05-06 中国科学院等离子体物理研究所 Multiple-objective optimization method and system capable of optimizing radiotherapy beam intensity distribution
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
CN105031820A (en) * 2015-09-01 2015-11-11 上海联影医疗科技有限公司 Intensity modulated radiation therapy reverse optimization method and device

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106902478A (en) * 2017-02-22 2017-06-30 天津大学 A kind of systematization radiotherapy biological effect appraisal procedure
CN110582325A (en) * 2017-05-03 2019-12-17 光线搜索实验室公司 System and method for ion-based radiotherapy treatment planning
CN110582325B (en) * 2017-05-03 2022-05-27 光线搜索实验室公司 System for ion-based radiotherapy treatment planning
CN107349531A (en) * 2017-07-04 2017-11-17 上海联影医疗科技有限公司 A kind of radiotherapy apparatus, storage medium, Ziye shape generation method and device
CN107349531B (en) * 2017-07-04 2019-12-13 上海联影医疗科技有限公司 radiotherapy equipment, storage medium, and method and device for generating subfield shape
CN110246562B (en) * 2018-03-07 2021-08-06 有方（合肥）医疗科技有限公司 Method, device and computer system for determining sub-beam intensity in radiotherapy system
CN110246562A (en) * 2018-03-07 2019-09-17 有方（合肥）医疗科技有限公司 Determine the method, apparatus and computer system of the sub-beam intensity in radiotherapy system
CN111437521A (en) * 2020-04-23 2020-07-24 华东医院 Non-uniform volume arc intensity modulation method
CN111785378B (en) * 2020-06-19 2022-04-26 新里程医用加速器(无锡)有限公司 Radiotherapy planning system and storage medium
CN111785378A (en) * 2020-06-19 2020-10-16 新里程医用加速器(无锡)有限公司 Radiotherapy planning system and storage medium
CN111986778A (en) * 2020-07-31 2020-11-24 上海联影医疗科技股份有限公司 Intensity modulated plan optimization system, device and storage medium
CN111986778B (en) * 2020-07-31 2024-02-20 上海联影医疗科技股份有限公司 Intensity-modulated plan optimization system, device and storage medium
CN112043974B (en) * 2020-09-08 2021-07-06 苏州雷泰医疗科技有限公司 Dynamic intensity modulation method and device based on orthogonal double-layer grating rotary sweep
CN112043974A (en) * 2020-09-08 2020-12-08 苏州雷泰医疗科技有限公司 Dynamic intensity modulation method and device based on orthogonal double-layer grating rotary sweep
WO2022052324A1 (en) * 2020-09-08 2022-03-17 苏州雷泰医疗科技有限公司 Dynamic intensity modulation method and apparatus based on orthogonal double-layer grating rotary sweeping
CN112972912A (en) * 2021-04-22 2021-06-18 北京大学第三医院（北京大学第三临床医学院） Non-coplanar radiotherapy beam incident path optimization method
CN112972912B (en) * 2021-04-22 2021-08-17 北京大学第三医院（北京大学第三临床医学院） Non-coplanar radiotherapy beam incident path optimization method
CN113521561A (en) * 2021-06-22 2021-10-22 苏州雷泰医疗科技有限公司 Method, storage medium and device for realizing reverse rotation intensity modulation of orthogonal double-layer grating
CN115966281A (en) * 2023-01-16 2023-04-14 迈胜医疗设备有限公司 Method, device, equipment and storage medium for generating arc radiotherapy plan
CN115966281B (en) * 2023-01-16 2024-02-09 迈胜医疗设备有限公司 Method, device, equipment and storage medium for generating arc radiotherapy plan

Similar Documents

Publication Publication Date Title
CN105617536A (en) Inverse intensity modulated arc radioterapy optimization method and device
EP3515557B1 (en) Controlling and shaping the dose distribution outside treatment targets in external-beam radiation treatments
CN105031820B (en) Intensity modulated radiation therapy reverse optimization method and device
Younge et al. Penalization of aperture complexity in inversely planned volumetric modulated arc therapy
CN102671309B (en) For determining the method and apparatus that irradiation is planned
US11103726B2 (en) Creating treatment field using initial field and patient specific geometry and achievable dose
US10076673B2 (en) Interactive dose gradient based optimization technique to control IMRT delivery complexity
US11291858B2 (en) Radiation treatment planning with multiple target subset optimization
Cho et al. Multisource modeling of flattening filter free (FFF) beam and the optimization of model parameters
Fuchs et al. A pencil beam algorithm for helium ion beam therapy
Wilson et al. A simple and robust trajectory‐based stereotactic radiosurgery treatment
CN109125952B (en) Convolution superposition energy deposition calculation method based on kernel model
JP2014054346A (en) Radiotherapy planning apparatus
Earl et al. Jaws‐only IMRT using direct aperture optimization
EP3524321A1 (en) Methods to optimize coverage for multiple targets simultaneously for radiation treatments
Lee et al. Monte Carlo based treatment planning for modulated electron beam radiation therapy
JP2018515274A (en) How to select beam geometry
Carlsson Combining segment generation with direct step‐and‐shoot optimization in intensity‐modulated radiation therapy
Kearney et al. A continuous arc delivery optimization algorithm for CyberKnife m6
CN105866821A (en) Method and device for obtaining beam energy spectrum, and dose distribution calculating method
CN107998519A (en) A kind of molecular dynamics re-optimization algorithm for IMRT
WO2022088341A1 (en) Radiation therapy control method and apparatus
CN105477789A (en) Dynamic intensity-modulated radiotherapy method based on quadratic programming model suppressing total beam-out time
Bedford et al. A phase space model of a Versa HD linear accelerator for application to Monte Carlo dose calculation in a real‐time adaptive workflow
Orton et al. The role of medical physicists and the AAPM in the development of treatment planning and optimization

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20160601