CN104268437B - It is a kind of to realize that rapid Optimum adjusts the wild method of hadron with point group - Google Patents
It is a kind of to realize that rapid Optimum adjusts the wild method of hadron with point group Download PDFInfo
- Publication number
- CN104268437B CN104268437B CN201410559541.8A CN201410559541A CN104268437B CN 104268437 B CN104268437 B CN 104268437B CN 201410559541 A CN201410559541 A CN 201410559541A CN 104268437 B CN104268437 B CN 104268437B
- Authority
- CN
- China
- Prior art keywords
- blade
- msub
- mrow
- group
- point group
- 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.)
- Active
Links
Landscapes
- Radiation-Therapy Devices (AREA)
Abstract
Realize that rapid Optimum adjusts the wild method of hadron with point group the invention discloses a kind of, solving prior art optimization needs successive ignition, the problem of time is longer.The present invention comprises the following steps：（1）, then will be for forming the blade point group of launched field it is determined that for forming the blade needed for launched field according to the 3D shape and the related anatomy relationship jeopardized between organ of target area；（2）The blade of each group is judged respectively using simulated annealing method.The present invention carries out a point group using before simulated annealing to blade, then optimizes calculating with a CPU core respectively to each group；Each a point group is judged with an independent object function, is changed the position of point group's Leaf, is calculated object function, judges whether the change of this blade can be received according to the rule that receives of solution.So as to reduce the optimization time on the basis of ensureing reduction Ziye number and machine hop count, improving efficiency of the practice；The operational capability of multi-core CPU is make use of simultaneously, and reduction is wasted.
Description
Technical field
Realize that rapid Optimum adjusts the wild method of hadron with point group the present invention relates to a kind of.
Background technology
Realize that preferable dosage is distributed by adjusting intensity of beam distribution, referred to as adjust strong technology.It is the present age to adjust strong technology
The most important developing direction of radiotherapy technology, is the core of radiotherapy technology, can effectively improve the cure rate of tumour and improve patient
Quality of life.Adjust at present mainly has following several ways by force：MLC static intensity modulatings, MLC dynamically adjust strong, arc to adjust strong, step-by-step movement
Tomotherapy.MLC static intensity modulatings method mainly has two kinds of substep reverse optimization method and direct Ziye optimization (DAO).In plan matter
In the case that amount is suitable, direct Ziye optimization method can obviously reduce Ziye number and machine hop count, improve efficiency of the practice.But calculate
Method optimization needs successive ignition, and the time is longer, and computational efficiency is low.
The CPU such as Intel, AMD productions commercial city employs multi-core technology to lift cpu performance, or even proposes crowd core CPU's
Concept.It means that to give full play to the performance of multi-core CPU, program just must be by the way of multi-thread concurrent calculating, tradition
Serial program will greatly waste the operational capability of multi-core CPU；So it can not preferably be applied to existing intensity modulated radiation therapy
In technology.
The content of the invention
Realize that rapid Optimum adjusts the wild method of hadron with point group it is an object of the invention to provide a kind of, solve prior art
Optimization needs successive ignition, the problem of time is longer.
To achieve these goals, the technical solution adopted by the present invention is as follows：
It is a kind of to realize that rapid Optimum adjusts the wild method of hadron with point group, comprise the following steps：
(1) according to the anatomy relationship between the direction of irradiation field, the 3D shape of target area and related crisis organ, it is determined that with
Blade needed for launched field is formed, then by the blade point group for forming launched field；
(2) blade of each group is judged respectively using simulated annealing method；
Simulated annealing method in the step (2), comprises the following steps：
A, initialization：Temperature T, initial solution state S are initialized, the iterations of each T values is L；
B, to K=1,2,3 ... L, perform C to E steps；
C, generation new explanation S '；
D, calculating increment △ F=F (S ')-F (S), wherein F (S) are evaluation function；
If E, △ F ＜ 0, receive S ' as new current solution, S ' is otherwise received as newly using probability exp (- △ F/T)
Current solution；
F, T are gradually reduced, and repeat B to E, try to achieve some new current solutions, judge whether all new current solutions are full
Sufficient end condition, if the current solution that then output meets end condition terminates program as optimal solution；Wherein, the T limit tend to 0；
The evaluation function is as follows：
Wherein, j is the numbering of histoorgan profile；ω is institutional framework weight；I is the intensity of a beam；dk dIt is to calculate
Dosage, dk pIt is prescribed dose, subscript k is volume elements ordinal number.
Further, the step (1) must is fulfilled for each group of blade Leafs in isocentric overall width minimum value:
R is the width of photon beam dose nuclear scattering in formula, and C is considers the revised width redundancy factor of pencil-beam, and i is leaf
The sequence number of piece, δiIt is the i-th blade in isocentric width, wminIt is blade in isocentric overall width minimum value.
Yet further, when the step C produces new explanation, each point group's Leaf quantity N is identical, and selection is each in point group
A blade carry out position change, be used for occurring number of the blade sequence number difference of position change for point group's blade in adjacent point of group
Measure N.
The present invention compared with prior art, with advantages below and beneficial effect：
The present invention carries out a point group using before simulated annealing to blade, and then each group is entered with a CPU core respectively
Row optimization is calculated；Each a point group is judged with an independent object function, is changed the position of point group's Leaf, is calculated mesh
Scalar functions, judge whether the change of this blade can be received according to the rule that receives of solution.So as to ensure reduce Ziye number and
The optimization time is reduced on the basis of machine hop count, raising efficiency of the practice；The operational capability of multi-core CPU, reduction are make use of simultaneously
Waste.
Brief description of the drawings
Fig. 1 is blade of the present invention point group's schematic diagram.
Fig. 2 is embodiments of the invention dose volume histogram.
Fig. 3 is schematic flow sheet of the invention.
Embodiment
The invention will be further described with reference to the accompanying drawings and examples, and embodiments of the present invention include but is not limited to
The following example.
Embodiment
As shown in figure 3, a kind of realize that rapid Optimum adjusts the wild method of hadron with point group, it is embodied as follows：
According to the direction of irradiation field, the 3D shape of target area and the related anatomy relationship jeopardized between organ, it is determined that being used for
The blade formed needed for launched field, then by the blade point group for forming launched field, launched field outer leafs are not considered, such as Fig. 1, figure
The blade numbering that launched field is formed in 1 is 5~24.Each group optimizes calculating with a CPU core respectively, and each group is only with one
Vertical object function is judged, is changed the position of group's Leaf, is calculated object function, this is judged according to the rule that receives of solution
Whether secondary blade changes and can be received, and is then specifically that the blade of each group is judged respectively using simulated annealing method.
Wherein, the quantity for dividing group is not to get The more the better, and according to the quantity and blade of CPU core, gains in depth of comprehension are wide in waiting
Spend to determine.The step (1) must is fulfilled for each group of blade Leafs in isocentric overall width minimum value:
R is the width of photon beam dose nuclear scattering in formula, and C is considers the revised width redundancy factor of pencil-beam, and i is leaf
The sequence number of piece, δiBe the i-th blade in isocentric width, W is blade in isocentric overall width minimum value.
According to actual conditions, it is exemplified below the blade in Fig. 1 being initialised and is divided into 4 points of groups, i.e. N value and is
5.First point of group's blade numbering is 5-9；Second point of group's blade numbering is 10-14；3rd point of group's blade numbering is 15-19；4th
It is 20-24 to divide group's blade numbering.
The simulated annealing method, comprises the following steps：
A, initialization：Initialize temperature T, initial solution state S, the iterations L of each T values；
B, to K=1,2,3 ... L, perform C to E steps；
C, generation new explanation S '；
D, calculating increment △ F=F (S ')-F (S), wherein F (S) are evaluation function；
If E, △ F ＜ 0, receive S ' as new current solution, S ' is otherwise received as newly using probability exp (- △ F/T)
Current solution；
F, T are gradually reduced, and repeat B to E, try to achieve some new current solutions, judge whether all new current solutions are full
Sufficient end condition, if then exporting current solution as optimal solution, terminates program；Wherein, the T limit tend to 0；
The evaluation function is as follows：
Wherein, j is the numbering of histoorgan profile；ω is institutional framework weight；I is the intensity of a beam；dk dIt is to calculate
Dosage, dk pIt is prescribed dose, subscript k is volume elements ordinal number.
In order that being independent of each other between each group, it is necessary to avoid choosing the situation that adjacent point of group's blade interacts, adopt
Mode is that a pair of blades in each group of selection carry out the leaf for being used to occur position change in position change, adjacent point of group
Piece sequence number difference is the quantity N of point group's blade, can so make the space requirement of score group minimum, for example, divide group 1 to select blade 6
During change, the blade numbering of other corresponding three points of groups is then 11,16,21.
End condition has：Circulating temperature T is less than setup parameter value, iteration number of success and reaches setup parameter value, target letter
Numerical value is less than setup parameter value etc..
The present invention is applied to personal computer (the core dominant frequency 3.0GHz of CPU 4, internal memory 4G, operating system win7), does one
Typical nasopharyngeal carcinoma adjusts strong case, solid with 7 directions (frame angle is respectively 0 °, 52 °, 104 °, 156 °, 208 °, 260 °, 312 °)
Determine open country and do the strong plan of tune to be tested, show that the dose volume histogram of plan is shown in Fig. 2, solid line is the agent regardless of each organ under group
Volume histogram is measured, the dose volume histogram of each organ when round dot is divides 4 groups.As can be seen from the figure group is divided and regardless of group
In the case of obtained effect of optimization it is basically identical, but optimization efficiency improves 1.5 times when launched field is divided into two groups, is divided into 3
Efficiency improves 2.2 times during group, and it is that efficiency improves 3.1 times to be divided into 4 groups.So as to illustrate the relatively existing method of the present invention its
Efficiency substantially increases, and is to increase substantially, and effectively shortens the tune strong time.
According to above-described embodiment, the present invention just can be realized well.What deserves to be explained is, before said structure design
Put, to solve same technical problem, some made in the present invention are used without substantial change or polishing
Technical scheme essence still as the present invention, therefore it should also be as within the scope of the present invention.
Claims (3)
1. a kind of realize that rapid Optimum adjusts the wild method of hadron with point group, it is characterised in that comprises the following steps：
(1) according to the anatomy relationship between the direction of irradiation field, the 3D shape of target area and related crisis organ, it is determined that for shape
Into the blade needed for launched field, then by the blade point group for forming launched field；
(2) blade of each group is judged respectively using simulated annealing method；
Simulated annealing method in the step (2), comprises the following steps：
A, initialization：Temperature T, initial solution state S are initialized, the iterations of each T values is L；
B, to K=1,2,3 ... L, perform C to E steps；
C, generation new explanation S '；
D, calculating increment △ F=F (S ')-F (S), wherein F (S) are evaluation function；
If E, △ F ＜ 0, receive S ' as new current solution, S ' is otherwise received as newly current using probability exp (- △ F/T)
Solution；
F, T are gradually reduced, and repeat B to E, try to achieve some new current solutions, judge whether all new current solutions meet end
Only condition, if then exporting the current solution for meeting end condition as optimal solution, terminates program；Wherein, the T limit tend to 0；
The evaluation function is as follows：
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<mi>j</mi>
</munder>
<msub>
<mi>&omega;</mi>
<mi>j</mi>
</msub>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>&Element;</mo>
<mi>&omega;</mi>
</mrow>
</munder>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>k</mi>
</msub>
<msup>
<msub>
<mi>d</mi>
<mi>k</mi>
</msub>
<mi>d</mi>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>d</mi>
<mi>k</mi>
</msub>
<mi>d</mi>
</msup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
Wherein, j is the numbering of histoorgan profile；ω is institutional framework weight；I is the intensity of a beam；dk dIt is the dosage calculated,
dk pIt is prescribed dose, subscript k is volume elements ordinal number.
2. a kind of method for realizing that rapid Optimum tune hadron is wild with point group according to claim 1, it is characterised in that described
Step (1) must is fulfilled for each group of blade Leafs in isocentric overall width minimum value:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>w</mi>
<mi>min</mi>
</msub>
<mo>=</mo>
<mi>C</mi>
<mo>&CenterDot;</mo>
<mi>r</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mstyle>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
</mstyle>
<msub>
<mi>&delta;</mi>
<mi>i</mi>
</msub>
<mo>&GreaterEqual;</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
R is the width of photon beam dose nuclear scattering in formula, and C is considers the revised width redundancy factor of pencil-beam, and i is blade
Sequence number, δiIt is the i-th blade in isocentric width, wminIt is blade in isocentric overall width minimum value.
3. a kind of method for realizing that rapid Optimum tune hadron is wild with point group according to claim 1 or 2, it is characterised in that
When the step C produces new explanation, each point group's Leaf quantity N is identical, and a blade in each point of group of selection carries out position
Change, be used for occurring quantity N of the blade sequence number difference of position change for point group's blade in adjacent point of group.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410559541.8A CN104268437B (en) | 2014-10-17 | 2014-10-17 | It is a kind of to realize that rapid Optimum adjusts the wild method of hadron with point group |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410559541.8A CN104268437B (en) | 2014-10-17 | 2014-10-17 | It is a kind of to realize that rapid Optimum adjusts the wild method of hadron with point group |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104268437A CN104268437A (en) | 2015-01-07 |
| CN104268437B true CN104268437B (en) | 2017-09-26 |
Family
ID=52159958
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410559541.8A Active CN104268437B (en) | 2014-10-17 | 2014-10-17 | It is a kind of to realize that rapid Optimum adjusts the wild method of hadron with point group |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104268437B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106902480B (en) * | 2017-03-07 | 2019-12-03 | 西安一体医疗科技有限公司 | A kind of parallel Quantum annealing target spot distribution calculation method |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6411675B1 (en) * | 2000-11-13 | 2002-06-25 | Jorge Llacer | Stochastic method for optimization of radiation therapy planning |
-
2014
- 2014-10-17 CN CN201410559541.8A patent/CN104268437B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6411675B1 (en) * | 2000-11-13 | 2002-06-25 | Jorge Llacer | Stochastic method for optimization of radiation therapy planning |
Non-Patent Citations (2)
| Title |
|---|
| 《基于多核的并行模拟退火放疗计划优化加速方法》;王慧亮;《北京大学学报（自然科学版）》;20081120;第44卷(第6期);第855页第1-2栏 * |
| 《调强放疗治疗计划制定中的改进模拟退火算法》;高磊;《医疗卫生装备》;20060815;第27卷(第8期);第23页第2栏 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104268437A (en) | 2015-01-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103083821B (en) | The system and method for a kind of direct optimization static intensity modulating Ziye shape and machine jumping figure thereof | |
| Xhaferllari et al. | Automated IMRT planning with regional optimization using planning scripts | |
| Gintz et al. | Initial evaluation of automated treatment planning software | |
| Togni et al. | Electromagnetic redesign of the HYPERcollar applicator: toward improved deep local head-and-neck hyperthermia | |
| US7987053B2 (en) | Monitor units calculation method for proton fields | |
| CN110124214A (en) | Intensity modulated radiation therapy plan optimization method and application based on predicted dose distribution guidance | |
| Bedford et al. | Constrained segment shapes in direct‐aperture optimization for step‐and‐shoot IMRT | |
| Cao et al. | Incorporating deliverable monitor unit constraints into spot intensity optimization in intensity-modulated proton therapy treatment planning | |
| Rowbottom et al. | Simultaneous optimization of beam orientations and beam weights in conformal radiotherapy | |
| CN105617536A (en) | Inverse intensity modulated arc radioterapy optimization method and device | |
| Widesott et al. | Is there a single spot size and grid for intensity modulated proton therapy? Simulation of head and neck, prostate and mesothelioma cases | |
| CN105031820A (en) | Intensity modulated radiation therapy reverse optimization method and device | |
| CN105709341A (en) | Intensity-modulated sub-field optimization method based on gradient method and floodfill algorithm | |
| Joy et al. | Dosimetric effects of jaw tracking in step‐and‐shoot intensity‐modulated radiation therapy | |
| CN104268437B (en) | It is a kind of to realize that rapid Optimum adjusts the wild method of hadron with point group | |
| Binny et al. | Effects of changing modulation and pitch parameters on tomotherapy delivery quality assurance plans | |
| Mahnam et al. | Simultaneous delivery time and aperture shape optimization for the volumetric-modulated arc therapy (VMAT) treatment planning problem | |
| CN110237439A (en) | Static intensity modulating Ziye optimization method for the orthogonal double-deck grating device | |
| CN107551411A (en) | For the probabilistic proton heavy ion intensity modulated radiation therapy robust Optimal methods of range | |
| Jensen et al. | Energy layer optimization strategies for intensity‐modulated proton therapy of lung cancer patients | |
| Wang et al. | An integrated solution of deep reinforcement learning for automatic IMRT treatment planning in non-small-cell lung cancer | |
| CN102247660A (en) | Reverse planning method for treatment plan and treatment plan system | |
| CN109248385B (en) | Radiation therapy plan optimization system based on Monte Carlo tree search | |
| Mortezazadeh et al. | Optimization of refueling cycle length by an enhanced PSO with novel mutation operator | |
| CN105561485A (en) | Radiotherapy treatment planning optimization method and device |
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 |