CN104338240A - Automatic optimization method for on-line self-adaption radiotherapy plan and device - Google Patents
Automatic optimization method for on-line self-adaption radiotherapy plan and device Download PDFInfo
- Publication number
- CN104338240A CN104338240A CN201410607268.1A CN201410607268A CN104338240A CN 104338240 A CN104338240 A CN 104338240A CN 201410607268 A CN201410607268 A CN 201410607268A CN 104338240 A CN104338240 A CN 104338240A
- Authority
- CN
- China
- Prior art keywords
- ziye
- radiotherapy
- original
- radiotherapy planning
- dose distribution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001959 radiotherapy Methods 0.000 title claims abstract description 224
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000005457 optimization Methods 0.000 title claims abstract description 49
- 230000005540 biological transmission Effects 0.000 claims abstract description 11
- 230000009191 jumping Effects 0.000 claims description 60
- 230000003044 adaptive effect Effects 0.000 claims description 35
- 230000003068 static effect Effects 0.000 claims description 35
- 238000012795 verification Methods 0.000 claims description 30
- 238000011404 fractionated radiotherapy Methods 0.000 claims description 23
- 210000000056 organ Anatomy 0.000 claims description 21
- 238000002560 therapeutic procedure Methods 0.000 claims description 10
- 238000004321 preservation Methods 0.000 claims description 9
- 238000007689 inspection Methods 0.000 claims description 7
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 231100000628 reference dose Toxicity 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 238000001514 detection method Methods 0.000 abstract 2
- 238000010606 normalization Methods 0.000 abstract 1
- 238000002721 intensity-modulated radiation therapy Methods 0.000 description 28
- 210000003484 anatomy Anatomy 0.000 description 5
- 238000000275 quality assurance Methods 0.000 description 5
- 238000002591 computed tomography Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000001225 therapeutic effect Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000002600 positron emission tomography Methods 0.000 description 3
- 206010028980 Neoplasm Diseases 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Landscapes
- Radiation-Therapy Devices (AREA)
Abstract
The invention discloses an automatic optimization method for an on-line self-adaption radiotherapy plan and a device. The method comprises the following steps: by reading an original radiotherapy plan and the information, such as acquired present gradation guide image and sketching target section, calculating a beam direction view of each frame angle, thereby acquiring the optimized newly generated sub-view and corresponding dose distribution, machine hop count and the like; finishing the on-line self-adaption automatic optimization of the radiotherapy plan; and meanwhile, performing secondary check on the dose distribution of the automatic radiotherapy plan through the calculation and assessment for a three-dimensional gamma index value; if failing to pass, performing normalization operation and updating the machine hop count of the new sub-view, thereby finishing the quality ensuring work. Furthermore, through the radiotherapy plan parameter transmission process check, the possibly present problem in a parameter transmitting and writing executing system is effectively avoided, and the safety monitoring is performed through real-time dosage detection during a plan executing process. The on-line self-adaption optimization of the radiotherapy plan, the quality ensuring and the full automation of parameter transmission check and real-time dosage detection are realized.
Description
Technical field
The present invention relates to optimization method and the device of radiotherapy planning, be specifically related to a kind of automatic optimization method and device of radiotherapy planning.
Background technology
At present, the process of tumour radiotherapy (Radiation Therapy) generates treatment plan based on the location CT of patient first before the treatment, then in therapeutic process subsequently, keep treatment plan constant, patient is implemented to the treatment of some gradation (Fraction).This Therapeutic mode does not consider the anatomical structure change of patient in therapeutic process, when great changes will take place for differences in patient, off-line self adaptation treatment (Offline Adaptive Radiation Therapy, Offline ART) is usually used to revise treatment plan.But off-line adaptive radiation therapy method uses traditional commercial therapeutic planning system at present, and its efficiency is low, takes time and effort, thus seldom adopt by hospital.
Online adaptive radiotherapy (Online Adaptive Radiation Therapy, Online ART) is before certain interval procedure of patient, obtains the anatomical structure image of patient, then generates the new treatment plan of patient fast.Current online adaptive radiotherapy, cannot use widely in clinical actual therapeutic, and main exist following problems and shortcoming:
1, online adaptive radiotheraping method is based on common commercial planning system, and when generating new treatment plan, need calculate and solve large-scale Rapid Dose Calculation and optimization problem, computation time is very long, cannot complete within the clinical acceptable time.
2, the method needs the optimizing process of radiation supervisor's plan of participating in the overall process, and utilizes setting and the adjustment of manually carrying out parameter, consuming time very long, and the quality for the treatment of plan is also directly confined to experience and the judgement of radiation supervisor.
3, finally new treatment plan needs on die body, again carry out the quality assurance (Quality Assurance, QA), cannot realize the optimization of online adaptive radiotherapy planning and the integration of automated quality guarantee.
Summary of the invention
Goal of the invention: in order to overcome the deficiencies in the prior art, the invention provides a kind of online adaptive radiotherapy planning automatic optimization method, a kind of online adaptive radiotherapy planning automatic optimizing equipment is provided simultaneously, the online adaptive optimization of radiotherapy planning can be completed rapidly, and complete the quality assurance of dose distribution quadratic search simultaneously, to adapt to generate the plan being more suitable for each interval procedure in conjunction with the change of anatomical structure, for you to choose, meet clinical needs well.
Technical scheme: for solving the problems of the technologies described above, online adaptive radiotherapy planning automatic optimization method provided by the invention, comprises the following steps:
Step 1), read original radiotherapy planning and comprise launched field information, dosage information and clinical information in interior raw information;
Step 2), obtain current gradation and comprise gradation navigational figure, current drawing target outline, the current profile jeopardizing organ in interior current gradation information;
Step 3), calculate the beam direction view of each frame angle;
Step 4), revise the shape of each frame angle Ziye based on beam direction view, the newly-generated Ziye after being optimized;
Step 5), adopt convolution superposition algorithm to calculate the dose distribution of each newly-generated Ziye;
Step 6), use the machine jumping figure of conjugate gradient rapid solving and each Ziye of optimization;
Step 7), carry out the quadratic search of automatic radiotherapy planning dose distribution: adopt Monte Carlo EGS4 method again to calculate the dose distribution of each newly-generated Ziye, and with step 5) dose distribution that obtains carries out three-dimensional gamma index value and calculates and assess, judge whether gamma percent of pass is less than default gamma percent of pass threshold value, in this way, then enter step 8), as no, then enter step 9);
Step 8), to step 7) adopt the dose distribution that calculates of Monte Carlo EGS4 method be normalized computing and upgrade the machine jumping figure of each Ziye;
Step 9), complete Automatic Optimal: the newly-generated Ziye after preservation optimization, the dose distribution of each newly-generated Ziye and machine jumping figure, as new radiotherapy planning, terminate.
As preferably, described step 1) in, described original radiotherapy planning is that static intensity modulating radiotherapy plan (Static IMRT) or arc adjust strong radiotherapy plan (Arc IMRT), wherein
The launched field information of original radiotherapy planning, comprising: frame angle, multi-blade collimator (Multi-Leaf Collimator, the MLC) sequence of each Ziye (Segment), machine jumping figure (Machine Unit, MU);
The dosage information of original radiotherapy planning, comprising: dose distribution (Dose Distribution) distribution and dose volume histogram (Dose Volume Histogram, DVH);
The clinical information of original radiotherapy planning, comprising: the original drawing target outline of patient and the original profile jeopardizing organ;
Described step 2) in, described gradation navigational figure comprises at least one in the images such as CT (Computed Tomography), Cone-Beam CT (Cone-Beam CT), ultrasonic (Ultrasound), PET (Positron Emission Tomography) or magnetic resonance (Magnetic Resonance, MR);
Described step 3) in, comprise the steps:
Step 3.1), by Siddon Ray-Tracing projection algorithm, calculate the beam direction view (BEV) of original drawing target outline at each frame angle fast;
Step 3.2), by Siddon Ray-Tracing projection algorithm, calculate the beam direction view (BEV) of current drawing target outline at each frame angle fast;
Described step 4) in, comprise the steps:
Step 4.1), when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), directly from original scheme file, read original scheme Ziye shape; When original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, from original radiotherapy planning, read the multi-blade collimator position at each control point, obtained the original scheme Ziye shape of two continuous control point intermediate stand angles by linear interpolation;
Step 4.2), by the beam direction view (BEV) of original drawing target outline and the geometrical relationship of original scheme Ziye shape, based on the beam direction view (BEV) of current drawing target outline, according to the criterion that current Ziye shape is consistent with the geometrical relationship of beam direction view (BEV) of current drawing target outline and the geometrical relationship of the beam direction view BEV of original drawing target outline and original Ziye, obtain current Ziye shape;
Step 4.3), when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), preserve current Ziye shape as the newly-generated Ziye after optimization; When original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, to be no more than the maximum movement speed of multi-blade collimator between adjacent Ziye and frame for constraint principle, the newly-generated Ziye after current Ziye shape amendment is optimized;
Described step 5) in, adopt convolution superposition algorithm to calculate the dose distribution D of each newly-generated Ziye
ij, D
ijfor a jth Ziye is to the dose contribution of i-th volume elements;
Described step 6) in, through type 1.1, uses conjugate gradient solve and optimize the machine jumping figure set of Ziye:
Make
Wherein, i is the index value of volume elements, and j is the index value of Ziye, and T is the set of target area, and S is the set of target area and crisis organ, N
tfor the volume elements quantity of each target area, N
sfor the volume elements quantity of each organ;
D
ijfor a jth Ziye is to the dose contribution of i-th volume elements;
X
jfor the machine jumping figure of each Ziye of optimization;
D
iit is the dose distribution of i-th volume elements;
P
iit is the reference dose distribution of i-th volume elements;
with
for kth time iteration organ weight factor;
Described step 7) in, described dose distribution quadratic search comprises the steps:
Step 7.1), adopt Monte Carlo EGS4 method again to calculate the dose distribution of each newly-generated Ziye;
Step 7.2), by step 7.1) dose distribution that calculates and step 5) dose distribution that obtains carries out three-dimensional gamma index value and calculates and assess, judge whether gamma percent of pass is less than default gamma percent of pass threshold value, in this way, then enter step 8), as no, then enter step 9);
Described step 8) in, according to the principle making the Gross Target Volume of 95% accept the prescribed dose of 100%, to step 7.1) adopt the dose distribution that calculates of Monte Carlo EGS4 method to be normalized computing to obtain proportionality factors lambda, by λ and step 6) the machine jumping figure linear multiplication of each Ziye that obtains, upgrade the machine jumping figure of each Ziye;
Described step 9) in, the step completing Automatic Optimal comprises: when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), newly-generated Ziye after preservation optimization, the dose distribution of each newly-generated Ziye and machine jumping figure generate new radiotherapy planning, terminate; When original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, newly-generated Ziye shape linear interpolation is gone back to the position at each control point in original radiotherapy planning, jointly preserve with the dose distribution of each newly-generated Ziye and machine jumping figure and generate new radiotherapy planning, terminate.
As further preferably, the step 9 of above-mentioned online adaptive radiotherapy planning automatic optimization method) the new radiotherapy planning that obtains comprises the following steps in transmitting procedure:
Step 10), radiotherapy planning parameter transmission process check: when new radiotherapy planning is transferred to verification system (the Record and Verify System of dosage execution, R & V) after, read (R & V) verification system and comprise multi-blade collimator leaf position, machine jumping figure, frame angle is in interior treatment parameter, compare with the parameter in new radiotherapy planning, judge that whether the treatment parameter in verification system is consistent with the parameter in new radiotherapy planning, as otherwise whether inspection record verification system (R & V) abnormal, in this way, then preserve new radiotherapy planning as this fractionated radiotherapy plan.
As further preferred, this fractionated radiotherapy plan described comprises the following steps in the process of implementation:
Step 11), real-time dose measurement: use electronic portal image device (Electronic Portal Imaging Device in the implementation of this fractionated radiotherapy plan described, EPID) multi-blade collimator position in therapeutic process is recorded, or use daily record (Log) the file real-time verification multi-blade collimator position of accelerator record, and compare with the multi-blade collimator in radiotherapy planning, real-time reconstruction 3-dimensional dose dose distribution, if the difference of the dose distribution of real-time reconstruction dose distribution and this fractionated radiotherapy plan, more than 3%, exports prompting message, and stop to perform.
Preferably, described step 7.2) in the calculating parameter of three-dimensional gamma index value be 3mm, 3%}, presetting gamma percent of pass threshold value is 95% or 96% or 97% or 98% or 99%; Certain above-mentioned parameter and threshold value can be arranged as required flexibly.
A kind of online adaptive radiotherapy planning automatic optimizing equipment that the present invention provides simultaneously, comprising:
Device 1), comprise launched field information, dosage information and the clinical information device in interior raw information for reading original radiotherapy planning;
Device 2), comprise gradation navigational figure, current drawing target outline, the current device of profile in interior current gradation information jeopardizing organ for obtaining current gradation;
Device 3), for calculating the device of the beam direction view of each frame angle;
Device 4), for revising the shape of each frame angle Ziye based on beam direction view, the device of the newly-generated Ziye after being optimized;
Device 5), for the device adopting convolution superposition algorithm to calculate the dose distribution of each newly-generated Ziye;
Device 6), for using the device of the machine jumping figure of conjugate gradient rapid solving and each Ziye of optimization;
Device 7), for carrying out the device of automatic radiotherapy planning dose distribution quadratic search, comprising:
For the device adopting Monte Carlo EGS4 method again to calculate the dose distribution of each newly-generated Ziye;
For will the dose distribution that again calculates of Monte Carlo EGS4 method and step 5 be adopted) dose distribution that obtains carries out three-dimensional gamma index value and calculates and assess, and judge whether gamma percent of pass is less than default gamma percent of pass threshold value and must assesses judgment means;
For according to assessment judgment means output valve, select jump to step 8) or step 9) selection redirect device;
Device 8), for step 7) dose distribution that adopts Monte Carlo EGS4 method to calculate is normalized computing and upgrades the device of the machine jumping figure of each Ziye;
Device 9), for preserving the newly-generated Ziye after optimization, the dose distribution of each newly-generated Ziye and the machine jumping figure device as new radiotherapy planning.
Preferably, described device 3) in, comprise as lower device:
Device 3.1), by Siddon Ray-Tracing projection algorithm, calculate the device of original drawing target outline at the beam direction view (BEV) of each frame angle fast;
Device 3.2) by Siddon Ray-Tracing projection algorithm, calculate the device of current drawing target outline at the beam direction view (BEV) of each frame angle fast;
Described device 4) in, comprise as lower device:
Device 4.1), when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), from original scheme file, directly read the device of original scheme Ziye shape; And when original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, from original radiotherapy planning file, read the multi-blade collimator position at each control point, obtained the device of the original scheme Ziye shape of two continuous control point intermediate stand angles by linear interpolation;
Device 4.2), by the beam direction view (BEV) of original drawing target outline and the geometrical relationship of original scheme Ziye shape, based on the beam direction view (BEV) of current drawing target outline, according to the criterion that the geometrical relationship of current Ziye shape and the geometrical relationship of the beam direction view (BEV) of current drawing target outline and the beam direction view (BEV) of original drawing target outline and original Ziye is consistent, obtain the device of current Ziye shape;
Device 4.3), when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), preserve the device of current Ziye shape as the newly-generated Ziye after optimization; When original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, to be no more than the maximum movement speed of multi-blade collimator between adjacent Ziye and frame for constraint principle, the device of the newly-generated Ziye after current Ziye shape amendment is optimized;
Described device 8) be, for according to the principle of prescribed dose making the Gross Target Volume of 95% accept 100%, to step 7.1) adopt the dose distribution that calculates of Monte Carlo EGS4 method to be normalized the device that computing obtains proportionality factors lambda, and for by λ and step 6) the machine jumping figure linear multiplication of each Ziye that obtains, upgrade the device of the machine jumping figure of each Ziye;
Described device 9) in, comprise when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), the newly-generated Ziye after preservation optimization, the dose distribution of each newly-generated Ziye and machine jumping figure generate the device of new radiotherapy planning; And when original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, newly-generated Ziye shape linear interpolation is gone back to the position at each control point in original radiotherapy planning, jointly preserve with the dose distribution of each newly-generated Ziye and machine jumping figure the device generating new radiotherapy planning.
As further advantageous embodiment, described online adaptive radiotherapy planning automatic optimizing equipment also comprises device 10), described device 10) be radiotherapy planning parameter transmission process check device, comprise when new radiotherapy planning be transferred to dosage perform verification system (Record and Verify System, R & V) after, read verification system and comprise multi-blade collimator leaf position, machine jumping figure, frame angle in interior treatment parameter, the device compared with the parameter in new radiotherapy planning; And judge that whether the treatment parameter in verification system is consistent with the parameter in new radiotherapy planning, as otherwise whether inspection record verification system (R & V) abnormal, in this way, then the judgement testing fixture of new radiotherapy planning as this fractionated radiotherapy plan is preserved.
As further advantageous embodiment, described online adaptive radiotherapy planning automatic optimizing equipment also comprises device 11), described device 11) be real time agent amount detecting device, be included in the implementation of this fractionated radiotherapy plan described and use electronic portal image device (Electronic Portal Imaging Device, EPID) multi-blade collimator position in therapeutic process is recorded, or use daily record (Log) the file real-time verification multi-blade collimator position of accelerator record, and compare with the multi-blade collimator in radiotherapy planning, the device of real-time reconstruction 3-dimensional dose dose distribution, and when real-time reconstruction dose distribution and the difference of the dose distribution of this fractionated radiotherapy plan export prompting message more than 3%, and stop the device of execution.
Beneficial effect: the invention provides a kind of online adaptive radiotherapy planning automatic optimization method, a kind of online adaptive radiotherapy planning automatic optimizing equipment is provided simultaneously, by reading the current gradation navigational figure of original radiotherapy planning and acquisition, the information such as drawing target outline, calculate the beam direction view of each frame angle, and then obtain the dose distribution of the newly-generated Ziye after optimizing and correspondence thereof, the contents such as machine jumping figure, while completing the online adaptive Automatic Optimal of radiotherapy planning, the dose distribution quadratic search that automatic radiotherapy planning is carried out in assessment is calculated by three-dimensional gamma index value, as quadratic search assessment is not passed through, then computing be normalized to dose distribution and upgrade the machine jumping figure of Ziye, difficulty action accomplishment guarantee work.
This method invention is in conjunction with the change of anatomical structure in patient process, calculating scheduling algorithm by fast dose makes whole radiotherapy planning optimizing process can complete within tens minutes even a few minutes, before each interval procedure starts, online adaptive Automatic Optimal and the Quality Assurance of whole radiotherapy planning can be completed fast, use for therapeutic choice.Further, the present invention is by the process check of radiotherapy planning parameter transmission, effectively avoid possibility produced problem in the process of parameter transmission and write executive system, further, by real-time dose measurement in plan implementation, relatively the dose distribution of real-time reconstruction dose distribution and this fractionated radiotherapy plan carries out risk and security monitoring, once occur that security risk is reminded in time and termination plan performs.Achieve the full-automation of the online adaptive Automatic Optimal of radiotherapy planning, quality assurance, parameter transmission inspection and real-time dose measurement.
Generally speaking, compared to based on the off-line self adaptation radiotheraping method of business radiotherapy planning system and existing online adaptive radiotheraping method, efficiency of the present invention is high, save valuable time and human cost, meet clinical needed for, can be applicable clinically, there is significant social meaning.
Accompanying drawing explanation
Fig. 1 is the flow chart of the method invention that embodiment 1 provides.
Detailed description of the invention
Below in conjunction with embodiment, the present invention is described in further detail, and this implementation column does not form restriction to the present invention.
Original radiotherapy planning is before certain course for the treatment of, first or before certain interval procedure performed, based on the image information of patient tumors and peripheral organs tissue thereof, through Target delineations, and the radiotherapy planning of confirmation.
First or before certain interval procedure perform after therapeutic process in, also will carry out interval procedure to patient, before such interval procedure performs, the present invention considers the change of patient anatomy in therapeutic process, be optimized to generate new radiotherapy planning to original radiotherapy planning, as this fractionated radiotherapy plan.
Embodiment 1: the online adaptive radiotherapy planning automatic optimization method that the present embodiment 1 provides, as shown in Figure 1, comprises the following steps:
Step 1), read original radiotherapy planning and comprise launched field information, dosage information and clinical information in interior raw information;
Wherein, described original radiotherapy planning is that static intensity modulating radiotherapy plan (Static IMRT) or arc adjust strong radiotherapy plan (Arc IMRT);
The launched field information of original radiotherapy planning comprises: frame angle, the multi-blade collimator sequence of each Ziye, machine jumping figure;
The dosage information of original radiotherapy planning comprises: dose distribution and dose volume histogram;
The clinical information of original radiotherapy planning comprises: the original drawing target outline of patient and the original profile jeopardizing organ;
Step 2), obtain current gradation and comprise gradation navigational figure, current drawing target outline, the current profile jeopardizing organ in interior current gradation information; Wherein said gradation navigational figure comprises at least one in CT, Cone-Beam CT, ultrasonic, PET or magnetic resonance;
Step 3), calculate the beam direction view of each frame angle, comprising:
Step 3.1), by Siddon Ray-Tracing projection algorithm, calculate the beam direction view (BEV) of original drawing target outline at each frame angle fast;
Step 3.2), by Siddon Ray-Tracing projection algorithm, calculate the beam direction view (BEV) of current drawing target outline at each frame angle fast;
Step 4), revise the shape of each frame angle Ziye based on beam direction view, the newly-generated Ziye after being optimized, comprising:
Step 4.1), when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), directly from original scheme file, read original scheme Ziye shape; When original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, from original radiotherapy planning, read the multi-blade collimator position at each control point, obtained the original scheme Ziye shape of two continuous control point intermediate stand angles by linear interpolation;
Step 4.2), by the beam direction view (BEV) of original drawing target outline and the geometrical relationship of original scheme Ziye shape, based on the beam direction view (BEV) of current drawing target outline, according to the criterion that the geometrical relationship of current Ziye shape and the geometrical relationship of the beam direction view (BEV) of current drawing target outline and the beam direction view (BEV) of original drawing target outline and original Ziye is consistent, obtain current Ziye shape;
Step 4.3), when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), preserve current Ziye shape as the newly-generated Ziye after optimization; When original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, to be no more than the maximum movement speed of multi-blade collimator between adjacent Ziye and frame for constraint principle, the newly-generated Ziye after current Ziye shape amendment is optimized;
Step 5), adopt convolution superposition algorithm to calculate the dose distribution of each newly-generated Ziye, be specially the dose distribution D adopting convolution superposition algorithm to calculate each newly-generated Ziye
ij, D
ijfor a jth Ziye is to the dose contribution of i-th volume elements; ;
Step 6), use the machine jumping figure of conjugate gradient rapid solving and each Ziye of optimization, be specially:
Through type 1.1, uses conjugate gradient solve and optimize the machine jumping figure set of Ziye:
Make
Wherein, i is the index value of volume elements, and j is the index value of Ziye, and T is the set of target area, and S is target area and crisis device
The set of official, N
tfor the volume elements quantity of each target area, N
sfor the volume elements quantity of each organ;
D
ijfor a jth Ziye is to the dose contribution of i-th volume elements;
X
jfor the machine jumping figure of each Ziye of optimization;
D
iit is the dose distribution of i-th volume elements;
P
iit is the reference dose distribution of i-th volume elements;
with
for kth time iteration organ weight factor;
Above-mentioned solving when calculating first with optimizing process (is 1.0 based on given one group of organ weight factor, also other values can be set to, as 0.5 or 1.5 or 2 or 5 or 10), in solving-optimizing iterative process subsequently, compare each and optimize organ original doses volume histogram DVH
pwith kth time iteration volume histogram DVH
kdifference in areas
through type (1.2) constantly updates organ weight factor
with
The above-mentioned iteration solved and optimize calculating terminates principle and is: as kth iteration dose volume histogram DVH
kwith original doses volume histogram DVH
pwhen difference in areas is less than preset area difference limen value (be 0.05, also can be set to other values as required, as 0.5 or 1.5 or 2 or 5 or 10), terminates the optimization of this step and calculate, the machine jumping figure set { x of the Ziye after the solving-optimizing obtained
j;
Step 7), carry out the quadratic search of automatic radiotherapy planning dose distribution, comprising:
Step 7.1), adopt Monte Carlo EGS4 method again to calculate the dose distribution of each newly-generated Ziye;
Step 7.2), by step 7.1) dose distribution that calculates and step 5) dose distribution that obtains carries out three-dimensional gamma index value and calculates and assess, judge whether gamma percent of pass is less than default gamma percent of pass threshold value, in this way, then enter step 8), as no, then enter step 9);
Step 8), to step 7) adopt the dose distribution that calculates of Monte Carlo EGS4 method be normalized computing and upgrade the machine jumping figure of each Ziye, that is: according to the principle making the Gross Target Volume of 95% accept the prescribed dose of 100%, to step 7.1) adopt the dose distribution that calculates of Monte Carlo EGS4 method to be normalized computing to obtain proportionality factors lambda, by step 6) the machine jumping figure of each Ziye that obtains and λ linear multiplication, obtain the machine jumping figure of each Ziye after upgrading;
Step 9), newly-generated Ziye after preservation optimization, the dose distribution of each newly-generated Ziye and machine jumping figure are as new radiotherapy planning, complete Automatic Optimal, comprise: when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), newly-generated Ziye after preservation optimization, the dose distribution of each newly-generated Ziye and machine jumping figure generate new radiotherapy planning, terminate; When original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, newly-generated Ziye shape linear interpolation is gone back to the position at each control point in original radiotherapy planning, jointly preserve with the dose distribution of each newly-generated Ziye and machine jumping figure and generate new radiotherapy planning, terminate.
In the present embodiment 1, above-mentioned steps 7.2) in the calculating parameter of three-dimensional gamma index value be that { 3mm, 3%}, presetting gamma percent of pass threshold value is 95%.
The online adaptive radiotherapy planning automatic optimizing equipment that embodiment 1 provides, comprising:
Device 1), comprise launched field information, dosage information and the clinical information device in interior raw information for reading original radiotherapy planning;
Device 2), comprise gradation navigational figure, current drawing target outline, the current device of profile in interior current gradation information jeopardizing organ for obtaining current gradation;
Device 3), for calculating the device of the beam direction view of each frame angle, comprising:
Device 3.1), by Siddon Ray-Tracing projection algorithm, calculate the device of original drawing target outline at the beam direction view (BEV) of each frame angle fast;
Device 3.2) by Siddon Ray-Tracing projection algorithm, calculate the device of current drawing target outline at the beam direction view (BEV) of each frame angle fast;
Device 4), for revising the shape of each frame angle Ziye based on beam direction view, the device of the newly-generated Ziye after being optimized, comprising:
Device 4.1), when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), from original scheme file, directly read the device of original scheme Ziye shape; And when original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, from original radiotherapy planning file, read the multi-blade collimator position at each control point, obtained the device of the original scheme Ziye shape of two continuous control point intermediate stand angles by linear interpolation;
Device 4.2), by the beam direction view (BEV) of original drawing target outline and the geometrical relationship of original scheme Ziye shape, based on the beam direction view (BEV) of current drawing target outline, according to the criterion that current Ziye shape is consistent with the geometrical relationship of beam direction view (BEV) of current drawing target outline and the geometrical relationship of the beam direction view BEV of original drawing target outline and original Ziye, obtain the device of current Ziye shape;
Device 4.3), when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), preserve the device of current Ziye shape as the newly-generated Ziye after optimization; When original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, to be no more than the maximum movement speed of multi-blade collimator between adjacent Ziye and frame for constraint principle, the device of the newly-generated Ziye after current Ziye shape amendment is optimized;
Device 5), for the device adopting convolution superposition algorithm to calculate the dose distribution of each newly-generated Ziye;
Device 6), for using the device of the machine jumping figure of conjugate gradient rapid solving and each Ziye of optimization;
Device 7), for carrying out the device of automatic radiotherapy planning dose distribution quadratic search, comprising:
For the device adopting Monte Carlo EGS4 method again to calculate the dose distribution of each newly-generated Ziye;
For will the dose distribution that again calculates of Monte Carlo EGS4 method and step 5 be adopted) dose distribution that obtains carries out three-dimensional gamma index value and calculates and assess, and judge whether gamma percent of pass is less than default gamma percent of pass threshold value and must assesses judgment means;
For according to assessment judgment means output valve, select jump to step 8) or step 9) selection redirect device;
Device 8), for step 7) dose distribution that adopts Monte Carlo EGS4 method to calculate is normalized computing and upgrades the device of the machine jumping figure of each Ziye, comprising:
For according to the principle of prescribed dose making the Gross Target Volume of 95% accept 100%, to step 7.1) adopt the dose distribution that calculates of Monte Carlo EGS4 method to be normalized the device that computing obtains proportionality factors lambda, and for by λ and step 6) the machine jumping figure linear multiplication of each Ziye that obtains, upgrade the device of the machine jumping figure of each Ziye;
Device 9), for preserving the newly-generated Ziye after optimization, the dose distribution of each newly-generated Ziye and the machine jumping figure device as new radiotherapy planning, comprise: comprise when original radiotherapy planning is the plan of static intensity modulating radiotherapy (Static IMRT), the newly-generated Ziye after preservation optimization, the dose distribution of each newly-generated Ziye and machine jumping figure generate the device of new radiotherapy planning; And when original radiotherapy planning be arc adjust strong radiotherapy plan (Arc IMRT) time, newly-generated Ziye shape linear interpolation is gone back to the position at each control point in original radiotherapy planning, jointly preserve with the dose distribution of each newly-generated Ziye and machine jumping figure the device generating new radiotherapy planning.
Embodiment 2: the online adaptive radiotherapy planning automatic optimization method that embodiment 2 provides and the method provided with embodiment 1 basically identical, something in common is not repeated, and difference is:
Described step 9) the new radiotherapy planning that obtains is further comprising the steps of in transmitting procedure:
Step 10), radiotherapy planning parameter transmission process check: when new radiotherapy planning is transferred to verification system (the Record and Verify System of dosage execution, R & V) after, read verification system and comprise multi-blade collimator leaf position, machine jumping figure, frame angle is in interior treatment parameter, compare with the parameter in new radiotherapy planning, judge that whether the treatment parameter in verification system is consistent with the parameter in new radiotherapy planning, as otherwise whether inspection record verification system (R & V) abnormal, in this way, then preserve new radiotherapy planning as this fractionated radiotherapy plan.
This fractionated radiotherapy plan described in commission comprises the following steps:
Step 11), real-time dose measurement: use electronic portal image device (Electronic Portal Imaging Device in the implementation of this fractionated radiotherapy plan described, EPID) multi-blade collimator position in therapeutic process is recorded, or use daily record (Log) the file real-time verification multi-blade collimator position of accelerator record, and compare with the multi-blade collimator in radiotherapy planning, real-time reconstruction 3-dimensional dose dose distribution, if the difference of the dose distribution of real-time reconstruction dose distribution and this fractionated radiotherapy plan, more than 3%, exports prompting message, and stop to perform.
Described step 7.2) in default gamma percent of pass threshold value be 98%.
The online adaptive radiotherapy planning automatic optimizing equipment that embodiment 2 provides and the device provided with embodiment 1 basically identical, something in common is not repeated, and difference is:
Also comprise device 10), described device 10) be radiotherapy planning parameter transmission process check device, comprise when new radiotherapy planning be transferred to dosage perform verification system (Record and Verify System, R & V) after, read verification system and comprise multi-blade collimator leaf position, machine jumping figure, frame angle in interior treatment parameter, the device compared with the parameter in new radiotherapy planning; And judge that whether the treatment parameter in verification system is consistent with the parameter in new radiotherapy planning, as otherwise whether inspection record verification system (R & V) abnormal, in this way, then the judgement testing fixture of new radiotherapy planning as this fractionated radiotherapy plan is preserved.
Also comprise device 11), described device 11) be real time agent amount detecting device, be included in the implementation of this fractionated radiotherapy plan described and use electronic portal image device (Electronic Portal Imaging Device, EPID) multi-blade collimator position in therapeutic process is recorded, or use daily record (Log) the file real-time verification multi-blade collimator position of accelerator record, and compare with the multi-blade collimator in radiotherapy planning, the device of real-time reconstruction 3-dimensional dose dose distribution; And when real-time reconstruction dose distribution and the difference of the dose distribution of this fractionated radiotherapy plan export prompting message more than 3%, and stop the device of execution.
Above implementation column does not form restriction to the present invention, and relevant staff is in the scope not departing from the technology of the present invention thought, and the various change carried out and amendment, all drop in protection scope of the present invention.
Claims (9)
1. an online adaptive radiotherapy planning automatic optimization method, is characterized in that comprising the following steps:
Step 1), read original radiotherapy planning and comprise launched field information, dosage information and clinical information in interior raw information;
Step 2), obtain current gradation and comprise gradation navigational figure, current drawing target outline, the current profile jeopardizing organ in interior current gradation information;
Step 3), calculate the beam direction view of each frame angle;
Step 4), revise the shape of each frame angle Ziye based on beam direction view, the newly-generated Ziye after being optimized;
Step 5), adopt convolution superposition algorithm to calculate the dose distribution of each newly-generated Ziye;
Step 6), use the machine jumping figure of conjugate gradient rapid solving and each Ziye of optimization;
Step 7), carry out the quadratic search of automatic radiotherapy planning dose distribution: adopt Monte Carlo EGS4 method again to calculate the dose distribution of each newly-generated Ziye, and with step 5) dose distribution that obtains carries out three-dimensional gamma index value and calculates and assess, judge whether gamma percent of pass is less than default gamma percent of pass threshold value, in this way, then enter step 8), as no, then enter step 9);
Step 8), to step 7) adopt the dose distribution that calculates of Monte Carlo EGS4 method be normalized computing and upgrade the machine jumping figure of each Ziye;
Step 9), complete Automatic Optimal: the newly-generated Ziye after preservation optimization, the dose distribution of each newly-generated Ziye and machine jumping figure, as new radiotherapy planning, terminate.
2. online adaptive radiotherapy planning automatic optimization method according to claim 1, is characterized in that:
Described step 1) in, described original radiotherapy planning is that the plan of static intensity modulating radiotherapy or arc adjust strong radiotherapy plan, wherein
The launched field information of original radiotherapy planning comprises: frame angle, the multi-blade collimator sequence of each Ziye, machine jumping figure;
The dosage information of original radiotherapy planning comprises: dose distribution and dose volume histogram;
The clinical information of original radiotherapy planning comprises: the original drawing target outline of patient and the original profile jeopardizing organ;
Described step 2) in, described gradation navigational figure comprises at least one in CT, Cone-Beam CT, ultrasonic, PET or magnetic resonance;
Described step 3) in, comprise the steps:
Step 3.1), by Siddon Ray-Tracing projection algorithm, calculate the beam direction view of original drawing target outline at each frame angle fast;
Step 3.2), by Siddon Ray-Tracing projection algorithm, calculate the beam direction view of current drawing target outline at each frame angle fast;
Described step 4) in, comprise the steps:
Step 4.1), when original radiotherapy planning is the plan of static intensity modulating radiotherapy, directly from original scheme file, read original scheme Ziye shape; When original radiotherapy planning be arc adjust strong radiotherapy plan time, from original radiotherapy planning, read the multi-blade collimator position at each control point, obtained the original scheme Ziye shape of two continuous control point intermediate stand angles by linear interpolation;
Step 4.2), by the beam direction view of original drawing target outline and the geometrical relationship of original scheme Ziye shape, based on the beam direction view of current drawing target outline, according to the criterion that the geometrical relationship of current Ziye shape and the geometrical relationship of the beam direction view of current drawing target outline and the beam direction view of original drawing target outline and original Ziye is consistent, obtain current Ziye shape;
Step 4.3), when original radiotherapy planning is the plan of static intensity modulating radiotherapy, preserve current Ziye shape as the newly-generated Ziye after optimization; When original radiotherapy planning be arc adjust strong radiotherapy plan time, to be no more than the maximum movement speed of multi-blade collimator between adjacent Ziye and frame for constraint principle, the newly-generated Ziye after being optimized to current Ziye shape amendment;
Described step 5) in, adopt convolution superposition algorithm to calculate the dose distribution D of each newly-generated Ziye
ij, D
ijfor a jth Ziye is to the dose contribution of i-th volume elements;
Described step 6) in, through type 1.1, uses conjugate gradient solve and optimize the machine jumping figure set of Ziye:
Make
Wherein, i is the index value of volume elements, and j is the index value of Ziye, and T is the set of target area, and S is the set of target area and crisis organ, N
tfor the volume elements quantity of each target area, N
sfor the volume elements quantity of each organ;
D
ijfor a jth Ziye is to the dose contribution of i-th volume elements;
X
jfor the machine jumping figure of each Ziye of optimization;
D
iit is the dose distribution of i-th volume elements;
P
iit is the reference dose distribution of i-th volume elements;
with
for kth time iteration organ weight factor;
Described step 7) in, described dose distribution quadratic search comprises the steps:
Step 7.1), adopt Monte Carlo EGS4 method again to calculate the dose distribution of each newly-generated Ziye;
Step 7.2), by step 7.1) dose distribution that calculates and step 5) dose distribution that obtains carries out three-dimensional gamma index value and calculates and assess, judge whether gamma percent of pass is less than default gamma percent of pass threshold value, in this way, then enter step 8), as no, then enter step 9);
Described step 8) in, according to the principle making the Gross Target Volume of 95% accept the prescribed dose of 100%, to step 7.1) adopt the dose distribution that calculates of Monte Carlo EGS4 method to be normalized computing to obtain proportionality factors lambda, by λ and step 6) the machine jumping figure linear multiplication of each Ziye that obtains, upgrade the machine jumping figure of each Ziye;
Described step 9) in, the step completing Automatic Optimal comprises: when original radiotherapy planning is the plan of static intensity modulating radiotherapy, newly-generated Ziye after preservation optimization, the dose distribution of each newly-generated Ziye and machine jumping figure generate new radiotherapy planning, terminate; When original radiotherapy planning be arc adjust strong radiotherapy plan time, newly-generated Ziye shape linear interpolation is gone back to the position at each control point in original radiotherapy planning, jointly preserve with the dose distribution of each newly-generated Ziye and machine jumping figure and generate new radiotherapy planning, terminate.
3. online adaptive radiotherapy planning automatic optimization method according to claim 2, is characterized in that: described step 9) the new radiotherapy planning that obtains comprises the following steps in transmitting procedure:
Step 10), radiotherapy planning parameter transmission process check: after new radiotherapy planning is transferred to the verification system of dosage execution, read verification system and comprise multi-blade collimator leaf position, machine jumping figure, frame angle in interior treatment parameter, compare with the parameter in new radiotherapy planning, judge that whether the treatment parameter in verification system is consistent with the parameter in new radiotherapy planning, as otherwise whether inspection record verification system abnormal, in this way, then new radiotherapy planning is preserved as this fractionated radiotherapy plan.
4. online adaptive radiotherapy planning automatic optimization method according to claim 3, is characterized in that: this fractionated radiotherapy plan described in commission comprises the following steps:
Step 11), real-time dose measurement: use multi-blade collimator position in electronic portal image device record therapeutic process in the implementation of this fractionated radiotherapy plan described, or use the journal file real-time verification multi-blade collimator position of accelerator record, and compare with the multi-blade collimator in radiotherapy planning, real-time reconstruction 3-dimensional dose dose distribution, if the difference of the dose distribution of real-time reconstruction dose distribution and this fractionated radiotherapy plan, more than 3%, exports prompting message, and stop to perform.
5. online adaptive radiotherapy planning automatic optimization method according to claim 2, it is characterized in that: described step 7.2) in the calculating parameter of three-dimensional gamma index value be { 3mm, 3%}, presetting gamma percent of pass threshold value is 95% or 96% or 97% or 98% or 99%.
6. an online adaptive radiotherapy planning automatic optimizing equipment, is characterized in that comprising:
Device 1), comprise launched field information, dosage information and the clinical information device in interior raw information for reading original radiotherapy planning;
Device 2), comprise gradation navigational figure, current drawing target outline, the current device of profile in interior current gradation information jeopardizing organ for obtaining current gradation;
Device 3), for calculating the device of the beam direction view of each frame angle;
Device 4), for revising the shape of each frame angle Ziye based on beam direction view, the device of the newly-generated Ziye after being optimized;
Device 5), for the device adopting convolution superposition algorithm to calculate the dose distribution of each newly-generated Ziye;
Device 6), for using the device of the machine jumping figure of conjugate gradient rapid solving and each Ziye of optimization;
Device 7), for carrying out the device of automatic radiotherapy planning dose distribution quadratic search, comprising:
For the device adopting Monte Carlo EGS4 method again to calculate the dose distribution of each newly-generated Ziye;
For will the dose distribution that again calculates of Monte Carlo EGS4 method and step 5 be adopted) dose distribution that obtains carries out three-dimensional gamma index value and calculates and assess, and judge whether gamma percent of pass is less than default gamma percent of pass threshold value and must assesses judgment means;
For according to assessment judgment means output valve, select jump to step 8) or step 9) selection redirect device;
Device 8), for step 7) dose distribution that adopts Monte Carlo EGS4 method to calculate is normalized computing and upgrades the device of the machine jumping figure of each Ziye;
Device 9), for preserving the newly-generated Ziye after optimization, the dose distribution of each newly-generated Ziye and the machine jumping figure device as new radiotherapy planning.
7. online adaptive radiotherapy planning automatic optimizing equipment according to claim 6, is characterized in that:
Described device 3) in, comprise as lower device:
Device 3.1), by Siddon Ray-Tracing projection algorithm, calculate the device of original drawing target outline at the beam direction view of each frame angle fast;
Device 3.2) by Siddon Ray-Tracing projection algorithm, calculate the device of current drawing target outline at the beam direction view of each frame angle fast;
Described device 4) in, comprise as lower device:
Device 4.1), when original radiotherapy planning is the plan of static intensity modulating radiotherapy, from original scheme file, directly read the device of original scheme Ziye shape; And when original radiotherapy planning be arc adjust strong radiotherapy plan time, from original radiotherapy planning file, read the multi-blade collimator position at each control point, obtained the device of the original scheme Ziye shape of two continuous control point intermediate stand angles by linear interpolation;
Device 4.2), by the beam direction view of original drawing target outline and the geometrical relationship of original scheme Ziye shape, based on the beam direction view of current drawing target outline, according to the criterion that the geometrical relationship of current Ziye shape and the geometrical relationship of the beam direction view of current drawing target outline and the beam direction view of original drawing target outline and original Ziye is consistent, obtain the device of current Ziye shape;
Device 4.3), when original radiotherapy planning is the plan of static intensity modulating radiotherapy, preserve the device of current Ziye shape as the newly-generated Ziye after optimization; When original radiotherapy planning be arc adjust strong radiotherapy plan time, to be no more than the maximum movement speed of multi-blade collimator between adjacent Ziye and frame for constraint principle, the device of the newly-generated Ziye after being optimized to current Ziye shape amendment;
Described device 8) be, for according to the principle of prescribed dose making the Gross Target Volume of 95% accept 100%, to step 7.1) adopt the dose distribution that calculates of Monte Carlo EGS4 method to be normalized the device that computing obtains proportionality factors lambda, and for by λ and step 6) the machine jumping figure linear multiplication of each Ziye that obtains, upgrade the device of the machine jumping figure of each Ziye;
Described device 9) in, comprise when original radiotherapy planning is the plan of static intensity modulating radiotherapy, the newly-generated Ziye after preservation optimization, the dose distribution of each newly-generated Ziye and machine jumping figure generate the device of new radiotherapy planning; And when original radiotherapy planning be arc adjust strong radiotherapy plan time, newly-generated Ziye shape linear interpolation is gone back to the position at each control point in original radiotherapy planning, jointly preserve with the dose distribution of each newly-generated Ziye and machine jumping figure the device generating new radiotherapy planning.
8. online adaptive radiotherapy planning automatic optimizing equipment according to claim 7, it is characterized in that: also comprise device 10), described device 10) be radiotherapy planning parameter transmission process check device, comprise after new radiotherapy planning is transferred to the verification system of dosage execution, read verification system and comprise multi-blade collimator leaf position, machine jumping figure, frame angle in interior treatment parameter, the device compared with the parameter in new radiotherapy planning; And judge that whether the treatment parameter in verification system consistent with the parameter in new radiotherapy planning, as otherwise whether inspection record verification system abnormal, in this way, then preserve the judgement testing fixture of new radiotherapy planning as this fractionated radiotherapy plan.
9. online adaptive radiotherapy planning automatic optimizing equipment according to claim 8, it is characterized in that: also comprise device 11), described device 11) be real time agent amount detecting device, be included in the implementation of this fractionated radiotherapy plan described and use multi-blade collimator position in electronic portal image device record therapeutic process, or use the journal file real-time verification multi-blade collimator position of accelerator record, and compare with the multi-blade collimator in radiotherapy planning, the device of real-time reconstruction 3-dimensional dose dose distribution; And when real-time reconstruction dose distribution and the difference of the dose distribution of this fractionated radiotherapy plan export prompting message more than 3%, and stop the device of execution.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410607268.1A CN104338240B (en) | 2014-10-31 | 2014-10-31 | Automatic optimization device for on-line self-adaption radiotherapy plan |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410607268.1A CN104338240B (en) | 2014-10-31 | 2014-10-31 | Automatic optimization device for on-line self-adaption radiotherapy plan |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104338240A true CN104338240A (en) | 2015-02-11 |
| CN104338240B CN104338240B (en) | 2017-02-22 |
Family
ID=52495213
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410607268.1A Active CN104338240B (en) | 2014-10-31 | 2014-10-31 | Automatic optimization device for on-line self-adaption radiotherapy plan |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104338240B (en) |
Cited By (33)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104636832A (en) * | 2015-02-27 | 2015-05-20 | 上海联影医疗科技有限公司 | Method and device for estimating dose distribution change, and method and system for direct machine parameter optimizing |
| CN104815392A (en) * | 2015-04-10 | 2015-08-05 | 石峰 | Interactive radiotherapy plan system optimization method and system thereof |
| CN104866928A (en) * | 2015-06-04 | 2015-08-26 | 上海联影医疗科技有限公司 | Optimization method and optimization system for adaptive radiotherapy |
| CN104933652A (en) * | 2015-04-27 | 2015-09-23 | 苏州敏宇医疗科技有限公司 | Cloud-computing based dose verification system and method of tumor radiotherapy |
| CN105031819A (en) * | 2015-08-25 | 2015-11-11 | 上海联影医疗科技有限公司 | Dosage optimization method and system |
| CN105999567A (en) * | 2016-06-22 | 2016-10-12 | 沈阳东软医疗系统有限公司 | Electric multi-leaf collimator leaf position control method and device |
| CN106846317A (en) * | 2017-02-27 | 2017-06-13 | 北京连心医疗科技有限公司 | A kind of feature based extracts the method for retrieving medicine image with Similarity matching |
| CN106920234A (en) * | 2017-02-27 | 2017-07-04 | 北京连心医疗科技有限公司 | A kind of method of the automatic radiotherapy planning of combined type |
| CN106932810A (en) * | 2017-04-01 | 2017-07-07 | 西安体医疗科技有限公司 | A kind of convolutional calculation method of gamma rays dosage |
| CN107510898A (en) * | 2017-10-11 | 2017-12-26 | 上海联影医疗科技有限公司 | Light field, Beam's eye view device and the verification method of Medical Devices |
| CN107731298A (en) * | 2017-09-28 | 2018-02-23 | 北京全域医疗技术有限公司 | Launched field method to set up and device based on radiotherapy planning system |
| CN107767947A (en) * | 2017-09-28 | 2018-03-06 | 北京全域医疗技术有限公司 | Radiotherapy planning appraisal procedure and device |
| CN108175953A (en) * | 2017-12-29 | 2018-06-19 | 苏州雷泰医疗科技有限公司 | A kind of strong optimization method of rotation tune and radiotherapy equipment based on dynamic raster |
| CN108434612A (en) * | 2018-01-05 | 2018-08-24 | 沈阳东软医疗系统有限公司 | A kind of QAT quality assurance test method and device |
| CN104866928B (en) * | 2015-06-04 | 2018-08-31 | 上海联影医疗科技有限公司 | Optimization method and system for adaptive radiation therapy |
| CN108648811A (en) * | 2018-05-11 | 2018-10-12 | 南方医科大学 | A kind of automatic radiotherapy planning quality evaluating method based on the prediction of indication item |
| CN109432611A (en) * | 2018-10-18 | 2019-03-08 | 北京连心医疗科技有限公司 | Generation method, equipment and the storage medium at extential form Dynamic wedge control point |
| CN109513119A (en) * | 2018-11-12 | 2019-03-26 | 武汉大学 | A kind of volume rotation intensity modulated radiation therapy plan collimator orientation optimization method |
| CN109513121A (en) * | 2018-12-28 | 2019-03-26 | 安徽大学 | Dose-guided adaptive radiotherapy plan re-optimization system and method |
| CN109999365A (en) * | 2018-01-05 | 2019-07-12 | 北京连心医疗科技有限公司 | A kind of radiotherapy reverse optimization method, equipment and storage medium mixing Monte Carlo |
| CN110246562A (en) * | 2018-03-07 | 2019-09-17 | 有方（合肥）医疗科技有限公司 | Determine the method, apparatus and computer system of the sub-beam intensity in radiotherapy system |
| CN110302475A (en) * | 2018-03-20 | 2019-10-08 | 北京连心医疗科技有限公司 | A kind of cloud Monte Carlo dose verifying analysis method, equipment and storage medium |
| CN111388879A (en) * | 2020-03-19 | 2020-07-10 | 上海联影医疗科技有限公司 | Radiation dose determination system, radiation dose determination device and storage medium |
| CN111589000A (en) * | 2020-05-27 | 2020-08-28 | 浙江省肿瘤医院 | Method for verifying parameters of medical linear accelerator |
| CN112354086A (en) * | 2020-10-16 | 2021-02-12 | 北京全域医疗技术集团有限公司 | Tumor radiotherapy positioning adjustment method and device and computer equipment |
| CN112618967A (en) * | 2020-12-17 | 2021-04-09 | 程明霞 | Device for adjusting radiotherapy dose distribution and generating corresponding radiotherapy plan |
| CN112926156A (en) * | 2021-02-25 | 2021-06-08 | 中科超精(南京)科技有限公司 | Direct subfield optimization device based on position-variable penumbra model and optimization method thereof |
| CN113827877A (en) * | 2021-09-15 | 2021-12-24 | 上海市胸科医院 | Method for realizing automatic dose verification based on accelerator log file |
| CN114010963A (en) * | 2021-11-08 | 2022-02-08 | 上海联影医疗科技股份有限公司 | System, method and computer-readable storage medium for dose determination |
| CN116328214A (en) * | 2023-05-30 | 2023-06-27 | 福建自贸试验区厦门片区Manteia数据科技有限公司 | Detection device for execution status of radiotherapy plan, electronic device, and storage medium |
| WO2023115261A1 (en) * | 2021-12-20 | 2023-06-29 | 西安大医集团股份有限公司 | Radiotherapy device control method and apparatus, device, and storage medium |
| CN116525064A (en) * | 2023-07-03 | 2023-08-01 | 福建自贸试验区厦门片区Manteia数据科技有限公司 | Radiation therapy plan detection device, electronic equipment and computer readable storage medium |
| WO2024152913A1 (en) * | 2023-01-16 | 2024-07-25 | 迈胜医疗设备有限公司 | Method and apparatus for generating arc therapy plan, and device and storage medium |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102306239B (en) * | 2011-07-22 | 2014-11-05 | 上海赛提菲克医疗器械有限公司 | Method for evaluating and optimizing radiotherapy dose based on cone beam CT (Computer Tomography) image CT value correction technology |
| CN104117151B (en) * | 2014-08-12 | 2017-01-25 | 章桦 | Optimization method of online self-adaption radiotherapy plan |
-
2014
- 2014-10-31 CN CN201410607268.1A patent/CN104338240B/en active Active
Non-Patent Citations (1)
| Title |
|---|
| AHUNBAY.E.E, ET AL: "An on-line replanning scheme for interfractional variations", 《MEDICAL PHYSICS》 * |
Cited By (53)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104636832B (en) * | 2015-02-27 | 2018-09-21 | 上海联影医疗科技有限公司 | Dosage changes in distribution evaluation method and device, direct Ziye optimization method and system |
| CN104636832A (en) * | 2015-02-27 | 2015-05-20 | 上海联影医疗科技有限公司 | Method and device for estimating dose distribution change, and method and system for direct machine parameter optimizing |
| CN104815392B (en) * | 2015-04-10 | 2018-01-26 | 石峰 | A kind of interactive radiotherapy treatment planning system optimization system |
| CN104815392A (en) * | 2015-04-10 | 2015-08-05 | 石峰 | Interactive radiotherapy plan system optimization method and system thereof |
| CN104933652A (en) * | 2015-04-27 | 2015-09-23 | 苏州敏宇医疗科技有限公司 | Cloud-computing based dose verification system and method of tumor radiotherapy |
| CN104866928B (en) * | 2015-06-04 | 2018-08-31 | 上海联影医疗科技有限公司 | Optimization method and system for adaptive radiation therapy |
| CN104866928A (en) * | 2015-06-04 | 2015-08-26 | 上海联影医疗科技有限公司 | Optimization method and optimization system for adaptive radiotherapy |
| CN105031819A (en) * | 2015-08-25 | 2015-11-11 | 上海联影医疗科技有限公司 | Dosage optimization method and system |
| CN105031819B (en) * | 2015-08-25 | 2018-11-06 | 上海联影医疗科技有限公司 | A kind of injectivity optimizing system |
| CN105999567A (en) * | 2016-06-22 | 2016-10-12 | 沈阳东软医疗系统有限公司 | Electric multi-leaf collimator leaf position control method and device |
| CN105999567B (en) * | 2016-06-22 | 2018-09-18 | 沈阳东软医疗系统有限公司 | A kind of the leaf position control method and device of electric multiblade collimator |
| CN106920234B (en) * | 2017-02-27 | 2021-08-27 | 北京连心医疗科技有限公司 | Combined automatic radiotherapy planning method |
| CN106920234A (en) * | 2017-02-27 | 2017-07-04 | 北京连心医疗科技有限公司 | A kind of method of the automatic radiotherapy planning of combined type |
| CN106846317A (en) * | 2017-02-27 | 2017-06-13 | 北京连心医疗科技有限公司 | A kind of feature based extracts the method for retrieving medicine image with Similarity matching |
| CN106846317B (en) * | 2017-02-27 | 2021-09-17 | 北京连心医疗科技有限公司 | Medical image retrieval method based on feature extraction and similarity matching |
| CN106932810A (en) * | 2017-04-01 | 2017-07-07 | 西安体医疗科技有限公司 | A kind of convolutional calculation method of gamma rays dosage |
| CN107731298A (en) * | 2017-09-28 | 2018-02-23 | 北京全域医疗技术有限公司 | Launched field method to set up and device based on radiotherapy planning system |
| CN107767947A (en) * | 2017-09-28 | 2018-03-06 | 北京全域医疗技术有限公司 | Radiotherapy planning appraisal procedure and device |
| CN107510898A (en) * | 2017-10-11 | 2017-12-26 | 上海联影医疗科技有限公司 | Light field, Beam's eye view device and the verification method of Medical Devices |
| CN107510898B (en) * | 2017-10-11 | 2021-04-16 | 上海联影医疗科技股份有限公司 | Light field and radiation field verification device and verification method for medical equipment |
| CN108175953A (en) * | 2017-12-29 | 2018-06-19 | 苏州雷泰医疗科技有限公司 | A kind of strong optimization method of rotation tune and radiotherapy equipment based on dynamic raster |
| CN108175953B (en) * | 2017-12-29 | 2020-07-21 | 苏州雷泰医疗科技有限公司 | Dynamic grating-based intensity modulation optimization method and radiotherapy equipment |
| CN108434612A (en) * | 2018-01-05 | 2018-08-24 | 沈阳东软医疗系统有限公司 | A kind of QAT quality assurance test method and device |
| CN109999365B (en) * | 2018-01-05 | 2021-04-20 | 北京连心医疗科技有限公司 | Computing device and storage medium for implementing radiotherapy reverse optimization method |
| CN109999365A (en) * | 2018-01-05 | 2019-07-12 | 北京连心医疗科技有限公司 | A kind of radiotherapy reverse optimization method, equipment and storage medium mixing Monte Carlo |
| CN108434612B (en) * | 2018-01-05 | 2020-07-10 | 东软医疗系统股份有限公司 | Quality assurance testing 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 |
| CN110302475A (en) * | 2018-03-20 | 2019-10-08 | 北京连心医疗科技有限公司 | A kind of cloud Monte Carlo dose verifying analysis method, equipment and storage medium |
| CN110302475B (en) * | 2018-03-20 | 2021-02-19 | 北京连心医疗科技有限公司 | Cloud Monte Carlo dose verification analysis method, equipment and storage medium |
| CN108648811A (en) * | 2018-05-11 | 2018-10-12 | 南方医科大学 | A kind of automatic radiotherapy planning quality evaluating method based on the prediction of indication item |
| CN109432611A (en) * | 2018-10-18 | 2019-03-08 | 北京连心医疗科技有限公司 | Generation method, equipment and the storage medium at extential form Dynamic wedge control point |
| CN109432611B (en) * | 2018-10-18 | 2021-02-12 | 北京连心医疗科技有限公司 | Method, equipment and storage medium for generating control point of extended dynamic wedge-shaped board |
| CN109513119A (en) * | 2018-11-12 | 2019-03-26 | 武汉大学 | A kind of volume rotation intensity modulated radiation therapy plan collimator orientation optimization method |
| CN109513121B (en) * | 2018-12-28 | 2021-01-01 | 安徽大学 | Dose-guided adaptive radiotherapy plan re-optimization system and method |
| CN109513121A (en) * | 2018-12-28 | 2019-03-26 | 安徽大学 | Dose-guided adaptive radiotherapy plan re-optimization system and method |
| CN111388879A (en) * | 2020-03-19 | 2020-07-10 | 上海联影医疗科技有限公司 | Radiation dose determination system, radiation dose determination device and storage medium |
| CN111388879B (en) * | 2020-03-19 | 2022-06-14 | 上海联影医疗科技股份有限公司 | Radiation dose determination system, radiation dose determination device and storage medium |
| CN111589000A (en) * | 2020-05-27 | 2020-08-28 | 浙江省肿瘤医院 | Method for verifying parameters of medical linear accelerator |
| CN112354086A (en) * | 2020-10-16 | 2021-02-12 | 北京全域医疗技术集团有限公司 | Tumor radiotherapy positioning adjustment method and device and computer equipment |
| CN112618967A (en) * | 2020-12-17 | 2021-04-09 | 程明霞 | Device for adjusting radiotherapy dose distribution and generating corresponding radiotherapy plan |
| CN112618967B (en) * | 2020-12-17 | 2024-03-15 | 程明霞 | Device for adjusting radiotherapy dosage distribution and generating corresponding radiotherapy plan |
| CN112926156A (en) * | 2021-02-25 | 2021-06-08 | 中科超精(南京)科技有限公司 | Direct subfield optimization device based on position-variable penumbra model and optimization method thereof |
| CN112926156B (en) * | 2021-02-25 | 2024-04-26 | 中科超精(南京)科技有限公司 | Direct sub-field optimization device and method based on position variable penumbra model |
| CN113827877A (en) * | 2021-09-15 | 2021-12-24 | 上海市胸科医院 | Method for realizing automatic dose verification based on accelerator log file |
| CN113827877B (en) * | 2021-09-15 | 2024-05-28 | 上海市胸科医院 | Method for realizing automatic dose verification based on accelerator log file |
| CN114010963A (en) * | 2021-11-08 | 2022-02-08 | 上海联影医疗科技股份有限公司 | System, method and computer-readable storage medium for dose determination |
| CN114010963B (en) * | 2021-11-08 | 2023-10-20 | 上海联影医疗科技股份有限公司 | System, method and computer readable storage medium for dose determination |
| WO2023115261A1 (en) * | 2021-12-20 | 2023-06-29 | 西安大医集团股份有限公司 | Radiotherapy device control method and apparatus, device, and storage medium |
| WO2024152913A1 (en) * | 2023-01-16 | 2024-07-25 | 迈胜医疗设备有限公司 | Method and apparatus for generating arc therapy plan, and device and storage medium |
| CN116328214A (en) * | 2023-05-30 | 2023-06-27 | 福建自贸试验区厦门片区Manteia数据科技有限公司 | Detection device for execution status of radiotherapy plan, electronic device, and storage medium |
| CN116328214B (en) * | 2023-05-30 | 2023-08-22 | 福建自贸试验区厦门片区Manteia数据科技有限公司 | Detection device for execution status of radiotherapy plan, electronic device, and storage medium |
| CN116525064A (en) * | 2023-07-03 | 2023-08-01 | 福建自贸试验区厦门片区Manteia数据科技有限公司 | Radiation therapy plan detection device, electronic equipment and computer readable storage medium |
| CN116525064B (en) * | 2023-07-03 | 2023-11-10 | 福建自贸试验区厦门片区Manteia数据科技有限公司 | Radiation therapy plan detection device, electronic equipment and computer readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104338240B (en) | 2017-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104338240B (en) | Automatic optimization device for on-line self-adaption radiotherapy plan | |
| Glide-Hurst et al. | Adaptive radiation therapy (ART) strategies and technical considerations: a state of the ART review from NRG oncology | |
| CN104117151B (en) | Optimization method of online self-adaption radiotherapy plan | |
| Meijers et al. | Evaluation of interplay and organ motion effects by means of 4D dose reconstruction and accumulation | |
| CN104548372B (en) | The dosage determining device of radiotherapy | |
| Beltran et al. | Planning target margin calculations for prostate radiotherapy based on intrafraction and interfraction motion using four localization methods | |
| US8222616B2 (en) | Method for adapting fractionation of a radiation therapy dose | |
| CN102671309B (en) | For determining the method and apparatus that irradiation is planned | |
| US8442287B2 (en) | Method and system for evaluating quality assurance criteria in delivery of a treatment plan | |
| US9511241B2 (en) | Irradiation planning for particle therapy | |
| US10417390B2 (en) | Methods and systems for radiotherapy treatment planning | |
| CA2616306A1 (en) | Method and system for processing data relating to a radiation therapy treatment plan | |
| CN103394167B (en) | A kind of Forecasting Methodology of the tumour radiotherapy complication based on biological effect dosage | |
| Nystrom et al. | Treatment planning for proton therapy: what is needed in the next 10 years? | |
| Bäumer et al. | Comprehensive clinical commissioning and validation of the RayStation treatment planning system for proton therapy with active scanning and passive treatment techniques | |
| McCowan et al. | Clinical implementation of a model-based in vivo dose verification system for stereotactic body radiation therapy–volumetric modulated arc therapy treatments using the electronic portal imaging device | |
| Boggula et al. | A new strategy for online adaptive prostate radiotherapy based on cone-beam CT | |
| Kejda et al. | Evaluation of the clinical feasibility of cone-beam computed tomography guided online adaption for simulation-free palliative radiotherapy | |
| Chen et al. | A clinically relevant online patient QA solution with daily CT scans and EPID-based in vivo dosimetry: A feasibility study on rectal cancer | |
| Varnava et al. | Dosimetric comparison between volumetric modulated arc therapy planning techniques for prostate cancer in the presence of intrafractional organ deformation | |
| Gueorguiev et al. | Clinical evaluation of a novel transmission detector for 3D quality assurance of IMRT and SBRT | |
| Kostovski | Artificial Intelligence in Radiotherapy | |
| Satherley | The impact of plan complexity on the accuracy of VMAT for the treatment of head and neck cancer | |
| Suvira | Dosimetric differences between scheduled and adapted plans generated from ethos adaptive radiotherapy for patients with prostate cancer | |
| Potter | Feasibility Study of the Flattening Filter Removal Design for Conventional Treatments Using a Linear Accelerator |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| GR01 | Patent grant | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20170126 Address after: 100084 Beijing Zhongguancun East Road, No. 1, building No. 8, ground floor, No. CB102-023, No. Applicant after: Beijing Lianxin Medical Technology Co Ltd Address before: Daming Road in Qinhuai District of Nanjing City, Jiangsu province 210022 No. 210, building 22, room 702, Nga Court 7 Applicant before: Zhang Hua |