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 PDF

Info

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
Application number
CN201410607268.1A
Other languages
Chinese (zh)
Other versions
CN104338240B (en
Inventor
李永宝
章桦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Lianxin Medical Technology Co Ltd
Original Assignee
章桦
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 章桦 filed Critical 章桦
Priority to CN201410607268.1A priority Critical patent/CN104338240B/en
Publication of CN104338240A publication Critical patent/CN104338240A/en
Application granted granted Critical
Publication of CN104338240B publication Critical patent/CN104338240B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

A kind of online adaptive radiotherapy planning automatic optimization method and device
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:
{ x j } = arg min { Σ s ∈ T a s k , - Σ i = 1 , . . . , N T [ max ( 0 , P i - d i ) ] 2 + Σ s ∈ S a s k , + Σ i = 1 , . . . , Ns [ max ( 0 , d i - P i ) ] 2 }
Make d i = Σ j D ij x j , x j ≥ 0 - - - ( 1.1 ) ;
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:
{ x j } = arg min { Σ s ∈ T a s k , - Σ i = 1 , . . . , N T [ max ( 0 , P i - d i ) ] 2 + Σ s ∈ S a s k , + Σ i = 1 , . . . , Ns [ max ( 0 , d i - P i ) ] 2 }
Make d i = Σ j D ij x j , x j ≥ 0 - - - ( 1.1 ) ;
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
α s k = α s 0 · f ( Δ s k ) - - - ( 1.2 ) ;
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:
{ x j } = arg min { Σ s ∈ T a s k , - Σ i = 1 , . . . , N T [ max ( 0 , P i - d i ) ] 2 + Σ s ∈ S a s k , + Σ i = 1 , . . . , N S [ max ( 0 , d i - P i ) ] 2 }
Make d i = Σ j D ij x j , x j ≥ 0 - - - ( 1.1 ) ;
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.
CN201410607268.1A 2014-10-31 2014-10-31 Automatic optimization device for on-line self-adaption radiotherapy plan Active CN104338240B (en)

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 (32)

* Cited by examiner, † Cited by third party
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
CN109513121A (en) * 2018-12-28 2019-03-26 安徽大学 A kind of dosage guidance adaptive radiation therapy plan re-optimization system and method
CN109513119A (en) * 2018-11-12 2019-03-26 武汉大学 A kind of volume rotation intensity modulated radiation therapy plan collimator orientation optimization 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

Family Cites Families (2)

* Cited by examiner, † Cited by third party
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

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
AHUNBAY.E.E, ET AL: "An on-line replanning scheme for interfractional variations", 《MEDICAL PHYSICS》 *

Cited By (52)

* Cited by examiner, † Cited by third party
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
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
CN109999365B (en) * 2018-01-05 2021-04-20 北京连心医疗科技有限公司 Computing device and storage medium for implementing radiotherapy reverse optimization method
CN110246562A (en) * 2018-03-07 2019-09-17 有方（合肥）医疗科技有限公司 Determine the method, apparatus and computer system of the sub-beam intensity in radiotherapy system
CN110302475B (en) * 2018-03-20 2021-02-19 北京连心医疗科技有限公司 Cloud Monte Carlo dose verification analysis method, equipment and storage medium
CN110302475A (en) * 2018-03-20 2019-10-08 北京连心医疗科技有限公司 A kind of cloud Monte Carlo dose verifying 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 安徽大学 A kind of dosage guidance adaptive radiation therapy 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
CN116328214B (en) * 2023-05-30 2023-08-22 福建自贸试验区厦门片区Manteia数据科技有限公司 Detection device for execution status of radiotherapy plan, electronic 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
CN116525064B (en) * 2023-07-03 2023-11-10 福建自贸试验区厦门片区Manteia数据科技有限公司 Radiation therapy plan detection device, electronic equipment and computer readable storage medium
CN116525064A (en) * 2023-07-03 2023-08-01 福建自贸试验区厦门片区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
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
US9511241B2 (en) Irradiation planning for particle therapy
US10417390B2 (en) Methods and systems for radiotherapy treatment planning
US20110112351A1 (en) Method and system for evaluating quality assurance criteria in delivery of a treatment plan
Cai et al. A practical implementation of physics quality assurance for photon adaptive radiotherapy
CA2616306A1 (en) Method and system for processing data relating to a radiation therapy treatment plan
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
Osewski et al. Clinical examples of 3D dose distribution reconstruction, based on the actual MLC leaves movement, for dynamic treatment techniques
Kejda et al. Evaluation of the clinical feasibility of cone-beam computed tomography guided online adaption for simulation-free palliative radiotherapy
Shrimali et al. Pitfalls and challenges to consider before setting up a lung cancer intensity-modulated radiotherapy service: A review of the reported clinical experience
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
Mijnheer Patient-Specific Quality Assurance: In Vivo 3D Dose Verification
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

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