CN105709341B - A kind of tune hadron open country optimization method based on gradient method and unrestrained water completion method - Google Patents
A kind of tune hadron open country optimization method based on gradient method and unrestrained water completion method Download PDFInfo
- Publication number
- CN105709341B CN105709341B CN201610029746.4A CN201610029746A CN105709341B CN 105709341 B CN105709341 B CN 105709341B CN 201610029746 A CN201610029746 A CN 201610029746A CN 105709341 B CN105709341 B CN 105709341B
- Authority
- CN
- China
- Prior art keywords
- ziye
- gradient
- launched field
- unrestrained water
- water completion
- 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.)
- Expired - Fee Related
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
- A61N5/1031—Treatment planning systems using a specific method of dose optimization
Abstract
The invention provides a kind of tune hadron open country optimization method based on gradient method and unrestrained water completion method, Optimized model is established with the conditions such as organ information, beam parameter and dosage (volume) constraint are jeopardized according to target area first, and initializes the initial Ziye of each launched field；Secondly, Segment weight is optimized using gradient method and obtains the intensity gradient on launched field；Again, the new Ziye shape of the launched field is calculated using unrestrained water completion method based on intensity gradient；Finally, second step is returned to until exiting optimization.Not only solving speed is fast by the present invention, and effect of optimization is good, and algorithmic stability, strong robustness.
Description
Technical field
The present invention relates to a kind of tune hadron open country optimization method based on gradient method and unrestrained water completion method, belong to intelligence computation side
Law technology field.
Background technology
Static intensity modulating radiotherapy is clinically the most frequently used treatment means at present, and it is each to be generally that optimization first obtains
Meet the launched field intensity of desired amount distribution on direction, Ziye segmentation then is carried out to the intensity, is converted to dynamic multi-leaf optical grating
The Ziye that can implement, this mode are referred to as " two step method ".Due to when intensity distribution optimizes, it is difficult to or do not account for moving
The hardware constraints limitation of state multi-leaf optical grating, and such constraint limitation only considers when Ziye is split, therefore cause Ziye to be split
After to produce Ziye number more, substantially increase treatment time, it is uncertain to add treatment.Long-term clinical practice shows,
The defects of " two step method " is present, it can be summarized as：(1) the Ziye number of each irradiation field is more, and total irradiation hop count (MU) is often
It is even more than ten times of the several times of conventional or conformal irradiation field, such case can not only increase the time of irradiation, also produce increase
The influence that the factors such as raw leakage is penetrated, head scattering are brought；(2) Ziye of some small shape very irregulars may be generated, these
Ziye caused error in implementation process is larger, and ray utilization rate is low, it is thus possible to influences the degree of accuracy of radiation dosage of patient；
(3) when preferable intensity distribution is converted into the intensity distribution that can implement, two intensity have difference, and this difference is often
Plan quality is set to be deteriorated, the plan of this clinical requirement met may no longer meet.
Adjust hadron open country optimization method also there are many researchs both at home and abroad at present, be mainly divided to two major classes at present：The first kind is to use
Intelligent algorithm optimizes, and this kind of algorithm is not easy to be absorbed in minimax solution, but it is non-convex problem to adjust hadron open country optimization problem, and the optimization time is past
It is past to be difficult to make us receiving；Second class is optimized using Deterministic Methods, and this kind of algorithm model is simple, and optimal speed is fast, but excellent
It is often unsatisfactory to change result.
The content of the invention
The technology of the present invention solves problem：For overcome the deficiencies in the prior art, it is an object of the invention to provide a kind of base
In the tune hadron open country optimization method of gradient method and unrestrained water completion method, fast and effective can solve to adjust hadron open country optimization problem, its is excellent
The requirement of radiotherapy planning formulation can be fully met by changing result, and enable the calculating time can be with those skilled in the art
In the range of receiving.
Technical scheme is as follows：A kind of tune hadron open country optimization method based on gradient method and unrestrained water completion method, it is real
Existing step is as follows：
(1) Ziye shape is initialized
Seen according to launched field direction and initialize corresponding Ziye shape on each launched field direction, the Ziye shape is being penetrated by target area
Outline extends out to obtain on wild direction；
(2) optimize Segment weight using gradient method and obtain the intensity gradient on launched field and optimize current Ziye using gradient method
Corresponding weight, while each launched field direction intensity gradient matrix after calculation optimization；
(3) the new Ziye shape of the launched field is calculated using unrestrained water completion method based on intensity gradient
Minimum value is found according to the intensity gradient matrix of step (2), if minimum value is more than threshold value, stopped, otherwise base
Unrestrained water completion method is carried out in minimum value position to calculate, finds out the Ziye for meeting Ziye shape need, wherein unrestrained water completion method
Boundary threshold be less than minimum value；
All launched fields are all done with unrestrained water completion method and calculates the new Ziye of generation, then obtains one group of new Ziye shape；
(4) step (2) and (3) alternate cycles are carried out, and until meeting stopping criterion, finally give one group of optimal Ziye shape
Shape and weight.
Unrestrained water completion method calculating is carried out based on minimum value position in the step (3) and realizes that step is as follows：
(21) extension gradient matrix minimum value is expert at continuous row, until border reaches threshold value；
(22) row based on current extensions, in extension columns, respective column is looked in the position of adjacent lines, finds out minimum value, minimum
Value is more than threshold value, then stops, and otherwise carries out step (21) extension；
(23) repeat step (22), until new Ziye shape generates.
The present invention compared with prior art the advantages of be：
(1) present invention fast and effective can solve to adjust hadron open country optimization problem, and its optimum results can fully meet radiotherapy
The requirement of plan, and enable and calculate the time in the range of those skilled in the art's acceptable.Therefore, this hair
It is bright that there is extremely important industrial application value and overcome still unsolved technical barrier.
(2) present invention can not only establish the model for adjusting hadron open country optimization problem exactly, and realizing can be rapidly and efficiently
Optimization adjusts hadron wild, meets actual user's request.
Brief description of the drawings
Fig. 1 is implementation process figure of the present invention.
Embodiment
As shown in figure 1, the present invention is embodied as：
(1) Ziye shape is initialized, it is specific as follows：
Each pencil beam is obtained in k-th of launched field to the dose contribution matrix D of target area sampled point according to Rapid Dose Calculationk, lead to
Cross formula (1) and obtain matrix Sk,1。
In formula, ε is a constant, and general value is 0.001, and (i, j) represents the i-th row jth row.Sk,1As k-th of launched field
First initial Ziye shape.
(2) Segment weight is optimized using gradient method and calculates the intensity gradient on launched field, it is specific as follows：
On each launched field direction is obtained after initial Ziye shape, using gradient method optimize corresponding to Segment weight.Make
Obtain each organ sampled point dose value after accelerator irradiation and try one's best and meet that the organ dose that doctor gives constrains.Gradient method optimization
The object function of wild weight is as follows：
Fobj(xk)=FTarget(xk)+FOAR(xk) (2)
Wherein, FTargetAnd FOARIt is target area and the object function for jeopardizing organ respectively；NTargetAnd NOARTarget area is represented respectively
With the number for jeopardizing organ；ωtarget,iAnd ωoar,iI-th of target area is represented respectively and jeopardizes the weight of organ, and weight represents each
The significance level of organ, general value is 0~1；ftarget,iAnd foar,iI-th of target area is represented respectively and jeopardizes the target of organ
Function；Ntarget,iAnd Noar,iI-th of target area is represented respectively and jeopardizes the sampled point number of organ；djRepresent j-th sampled point
Dosage；δ represents that if the dosage of sampled point in restriction range, takes 0, otherwise takes 1；Dtarget,iRepresent the agent of i-th of target area
Amount constraint, if sampling dose point is less than target dose lower limit, then Dtarget,iTarget dose lower limit is equal to, if sampled point
The dosage about target dose upper limit, then Dtarget,iThe target dose upper limit is equal to, if sampling dose point is somebody's turn to do in constraint
Item is 0, is not involved in calculating；Doar,iRepresent i-th of upper dosage limit for jeopardizing organ；NrayRepresent the pencil beam number in launched field；
aj,mRepresent Dose Effect of m-th of pencil beam to j-th of sampled point；XmRepresent the intensity of m-th of pencil beam；After M is discrete
Launched field number；xkThe weight of k-th of launched field is represented, is the parameter for needing to optimize；segkPressed for the matrix of Ziye 0/1 of k-th of launched field
According to the vector of row major arrangement.
The gradient matrix T in each launched field direction can be produced in last gradient method optimization processk, k k-th of launched field of expression.
(3) the new Ziye shape of the launched field is calculated using unrestrained water completion method based on intensity gradient
According to the intensity gradient matrix T of k-th of launched fieldk, find minimum value τ(i,j), wherein, i, j represent respectively the i-th row and
Jth arranges.If τ(i,j)More than -0.01, then illustrate the Ziye that k-th of launched field intensity gradient increases newly without needs, stop.Otherwise,
With τ(i,j)As seed, unrestrained water completion method extension is carried out by the way of four-way connection.Need satisfaction following several in expansion process
Point limitation：
There can only be a continuum per a line；
All launched fields are all done and once overflow the new Ziye of water completion method generation, then obtain one group of new Ziye shape, then pass through
Gradient method optimizes weight corresponding to each Ziye, while calculates the intensity gradient of launched field, and unrestrained water completion method and gradient method are alternately followed
Ring is carried out, until meeting stopping criterion.
(4) one group of optimal Ziye shape and weight are finally given.
It is described in detail again with reference to specific example.
The present invention can quickly obtain optimal tune hadron open country, after the parameter that user sets correlation, pass through gradient method
One group of Ziye shape and weight are quickly obtained with unrestrained water completion method, for peripheral type carcinoma of lung radiotherapy case, specific implementation step is such as
Under：
(1) target dose is contributed according to launched field pencil beam, initialized under Ziye shape：
Peripheral type carcinoma of lung example adds 5 launched fields, and each pencil beam is obtained in k-th of launched field to target area according to Rapid Dose Calculation
The dose contribution matrix D of sampled pointk, DkThe i-th row jth row represent the pencil beam to all sampled point dose contribution summations in target area,
Matrix S is obtained by formula (1)k,1。
In formula, ε is a constant, and general value is 0.001, and (i, j) represents the i-th row jth row.Sk,1As k-th of launched field
First initial Ziye shape.In general, Sk,1Meet Ziye requirement.
(2) Segment weight is optimized using gradient method and calculates the intensity gradient of launched field, it is specific as follows：
After the initial Ziye shape of 5 launched fields is obtained, using gradient method optimize corresponding to Segment weight.So that accelerator
Each organ sampled point dose value after irradiation, which is tried one's best, meets that the organ dose that doctor gives constrains.Gradient method optimization Segment weight
Object function is as follows：
Fobj(xk)=FTarget(xk)+FOAR(xk) (11)
Wherein, FTargetAnd FOARIt is target area and the object function for jeopardizing organ respectively；NTargetAnd NOARTarget area is represented respectively
With the number for jeopardizing organ；ωtarget,iAnd ωoar,iI-th of target area is represented respectively and jeopardizes the weight of organ, and weight represents each
The significance level of organ, general value is 0~1；ftarget,iAnd foar,iI-th of target area is represented respectively and jeopardizes the target of organ
Function；Ntarget,iAnd Noar,iI-th of target area is represented respectively and jeopardizes the sampled point number of organ；djRepresent j-th sampled point
Dosage；δ represents that if the dosage of sampled point in restriction range, takes 0, otherwise takes 1；Dtarget,iRepresent the agent of i-th of target area
Amount constraint, if sampling dose point is less than target dose lower limit, then Dtarget,iTarget dose lower limit is equal to, if sampled point
The dosage about target dose upper limit, then Dtarget,iThe target dose upper limit is equal to, if sampling dose point is somebody's turn to do in constraint
Item is 0, is not involved in calculating；Doar,iRepresent i-th of upper dosage limit for jeopardizing organ；NrayRepresent the pencil beam number in launched field；
aj,mRepresent Dose Effect of m-th of pencil beam to j-th of sampled point；XmRepresent the intensity of m-th of pencil beam；After M is discrete
Launched field number；xkThe weight of k-th of launched field is represented, is the parameter for needing to optimize；segkPressed for the matrix of Ziye 0/1 of k-th of launched field
According to the vector of row major arrangement.
The gradient matrix T in each launched field direction can be produced in last gradient method optimization processk, k k-th of launched field of expression.
(3) the new Ziye shape of the launched field is calculated using unrestrained water completion method based on intensity gradient
According to the intensity gradient matrix T of k-th of launched fieldk, find minimum value τ(i,j), wherein, i, j represent respectively the i-th row and
Jth arranges.If τ(i,j)More than -0.01, then illustrate the Ziye that k-th of launched field intensity gradient increases newly without needs, stop.Otherwise,
With τ(i,j)As seed, unrestrained water completion method extension is carried out by the way of four-way connection.Need satisfaction following several in expansion process
Point limitation：
Boundary threshold is τ(i,j)1/3；
There can only be a continuum per a line；
All launched fields are all done and once overflow the new Ziye of water completion method generation, then obtain one group of new Ziye shape, then pass through
Gradient method optimizes weight corresponding to each Ziye, while calculates the intensity gradient of launched field, and unrestrained water completion method and gradient method are alternately followed
Ring is carried out, until meeting stopping criterion.
(4) one group of optimal Ziye shape and weight are finally given.
Above example is provided just for the sake of the description purpose of the present invention, and is not intended to limit the scope of the present invention.This
The scope of invention is defined by the following claims.The various equivalent substitutions that do not depart from spirit and principles of the present invention and make and repair
Change, all should cover within the scope of the present invention.
Claims (2)
1. a kind of tune hadron open country optimization method based on gradient method and unrestrained water completion method, it is characterised in that realize that step is as follows：
(1) Ziye shape is initialized
Corresponding Ziye shape on each launched field direction is initialized according to launched field direction, the Ziye shape is by target area in launched field direction
Upper outline extends out to obtain；
(2) Segment weight is optimized using gradient method and obtains the intensity gradient on launched field
Optimize weight corresponding to current Ziye, while each launched field direction intensity gradient square after calculation optimization using gradient method
Battle array；
(3) the new Ziye shape of the launched field is calculated using unrestrained water completion method based on intensity gradient
Minimum value is found according to the intensity gradient matrix of step (2), if minimum value is more than threshold value, stopped, otherwise based on most
Small value position carries out unrestrained water completion method and calculated, and finds out the Ziye for meeting Ziye shape need, wherein the side of unrestrained water completion method
Boundary's threshold value is less than minimum value；
All launched fields are all done with unrestrained water completion method and calculates the new Ziye of generation, then obtains one group of new Ziye shape；
(4) step (2) and (3) alternate cycles are carried out, until meeting stopping criterion, finally give one group of optimal Ziye shape and
Weight.
2. according to tune hadron open country optimization method of the claim 1 based on gradient method and unrestrained water completion method, it is characterised in that：The step
Suddenly unrestrained water completion method calculating is carried out based on minimum value position in (3) and realizes that step is as follows：
(21) extension gradient matrix minimum value is expert at continuous row, until border reaches threshold value；
(22) row based on current extensions, in extension columns, respective column is looked in the position of adjacent lines, finds out minimum value, and minimum value is big
In threshold value, then stop, otherwise carrying out step (21) extension；
(23) repeat step (22), until new Ziye shape generates.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610029746.4A CN105709341B (en) | 2016-01-15 | 2016-01-15 | A kind of tune hadron open country optimization method based on gradient method and unrestrained water completion method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610029746.4A CN105709341B (en) | 2016-01-15 | 2016-01-15 | A kind of tune hadron open country optimization method based on gradient method and unrestrained water completion method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105709341A CN105709341A (en) | 2016-06-29 |
| CN105709341B true CN105709341B (en) | 2018-03-27 |
Family
ID=56147132
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610029746.4A Expired - Fee Related CN105709341B (en) | 2016-01-15 | 2016-01-15 | A kind of tune hadron open country optimization method based on gradient method and unrestrained water completion method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105709341B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106169028B (en) * | 2016-07-05 | 2020-11-03 | 西南医科大学附属医院 | System for constructing prediction of absorption dose of organs at risk in intensity modulated radiotherapy |
| CN106730411B (en) * | 2016-12-21 | 2019-11-19 | 上海联影医疗科技有限公司 | Radiotherapy treatment planning optimization method and system |
| CN107297033B (en) * | 2017-07-13 | 2019-02-26 | 中科超精（安徽）科技有限公司 | Freeze the radiotherapy pencil beam weight optimization acceleration method and system of pencil beam |
| CN107823806B (en) * | 2017-09-15 | 2019-11-19 | 中北大学 | A kind of method and system for the direct Ziye optimization of intensity-modulated radiation therapy |
| CN108211136B (en) * | 2017-12-29 | 2021-01-15 | 苏州雷泰医疗科技有限公司 | Dynamic grating-based rotary intensity modulation optimization method and radiotherapy equipment |
| CN108771795B (en) * | 2018-06-29 | 2021-03-09 | 上海联影医疗科技股份有限公司 | Device for determining radiation distribution of radiation therapy and radiation therapy system |
| CN110237439B (en) * | 2018-12-21 | 2020-12-18 | 苏州雷泰医疗科技有限公司 | Static intensity modulated subfield optimization method for orthogonal double-layer grating device |
| CN110289075B (en) * | 2019-04-29 | 2022-04-22 | 中北大学 | Direct subfield optimization method and system based on fuzzy entropy |
| CN110215623B (en) * | 2019-06-11 | 2020-11-10 | 苏州雷泰智能科技有限公司 | Method and device for realizing orthogonal double-layer grating rotation intensity modulation based on subfield optimization |
| CN112999531B (en) * | 2021-02-25 | 2022-06-24 | 中科超精(南京)科技有限公司 | Dynamic optimization system based on respiratory motion |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1926552A (en) * | 2003-12-12 | 2007-03-07 | 西安大略大学 | Method and system for optimizing dose delivery of radiation |
| 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 |
| CN104318122A (en) * | 2014-11-06 | 2015-01-28 | 中国科学院合肥物质科学研究院 | Disturbance-based intensity modulated arc therapy optimization method |
-
2016
- 2016-01-15 CN CN201610029746.4A patent/CN105709341B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1926552A (en) * | 2003-12-12 | 2007-03-07 | 西安大略大学 | Method and system for optimizing dose delivery of radiation |
| 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 |
| CN104318122A (en) * | 2014-11-06 | 2015-01-28 | 中国科学院合肥物质科学研究院 | Disturbance-based intensity modulated arc therapy optimization method |
Non-Patent Citations (1)
| Title |
|---|
| 基于改进区域生长算法的肝脏分割方法研究;刘航，汪东;《中国医学物理学杂志》;20140930;第31卷;第5204-5207页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105709341A (en) | 2016-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105709341B (en) | A kind of tune hadron open country optimization method based on gradient method and unrestrained water completion method | |
| US10850122B2 (en) | Optimization methods for radiation therapy planning | |
| CN110124214B (en) | Predictive dose distribution guided intensity modulated radiotherapy plan optimization systems, media and devices | |
| CN106029170B (en) | The establishment of increment treatment plan | |
| CN104318122B (en) | Disturbance-based intensity modulated arc therapy optimization method | |
| CN105413068B (en) | Solve the optimization method and radiotherapy equipment of flux pattern | |
| Witte et al. | IMRT optimization including random and systematic geometric errors based on the expectation of TCP and NTCP | |
| Fiege et al. | PARETO: A novel evolutionary optimization approach to multiobjective IMRT planning | |
| An et al. | Robust treatment planning with conditional value at risk chance constraints in intensity‐modulated proton therapy | |
| CN110327554B (en) | Intensity modulated radiotherapy plan optimization method based on predicted dose distribution guidance and application | |
| Bedford et al. | Constrained segment shapes in direct‐aperture optimization for step‐and‐shoot IMRT | |
| Witte et al. | Dealing with geometric uncertainties in dose painting by numbers: introducing the ΔVH | |
| Nguyen et al. | A comparison of Monte Carlo dropout and bootstrap aggregation on the performance and uncertainty estimation in radiation therapy dose prediction with deep learning neural networks | |
| CN109453473A (en) | Radiotherapy treatment planning system, determining device and storage medium | |
| CN109499012A (en) | A kind of TPS automatic Iterative optimization algorithm of optimization dosage guidance | |
| CN105617536A (en) | Inverse intensity modulated arc radioterapy optimization method and device | |
| US11103728B2 (en) | Methods and systems for irradiation based on a fluence map | |
| Barragán Montero et al. | Performance of a hybrid Monte Carlo‐Pencil Beam dose algorithm for proton therapy inverse planning | |
| Zhu et al. | Search for IMRT inverse plans with piecewise constant fluence maps using compressed sensing techniques | |
| Yue et al. | Dose prediction via distance-guided deep learning: Initial development for nasopharyngeal carcinoma radiotherapy | |
| US20140330064A1 (en) | Station parameter optimized radiation therapy (SPORT): a novel scheme for treatment planning and delivery in radiation therapy | |
| CN108618796A (en) | A kind of Monte Carlo scattered photon analogy method of task based access control driving | |
| CN105561485B (en) | The optimization device of radiotherapy treatment planning | |
| Zhuang et al. | Dose-volume histogram prediction in volumetric modulated arc therapy for nasopharyngeal carcinomas based on uniform-intensity radiation with equal angle intervals | |
| D D'Souza et al. | A nested partitions framework for beam angle optimization in intensity-modulated radiation therapy |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180327 Termination date: 20210115 |