CN104338240B - Automatic optimization device for on-line self-adaption radiotherapy plan - Google Patents
Automatic optimization device for on-line self-adaption radiotherapy plan Download PDFInfo
- Publication number
- CN104338240B CN104338240B CN201410607268.1A CN201410607268A CN104338240B CN 104338240 B CN104338240 B CN 104338240B CN 201410607268 A CN201410607268 A CN 201410607268A CN 104338240 B CN104338240 B CN 104338240B
- 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.)
- Active
Links
Landscapes
- Radiation-Therapy Devices (AREA)
Abstract
The invention discloses an automatic optimization device for an on-line self-adaption radiotherapy plan. 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 the optimization method of radiotherapy planning and device are and in particular to a kind of automatic optimization method of radiotherapy planning
And device.
Background technology
At present, the process of tumour radiotherapy (Radiation Therapy) is determining based on patient before the treatment first
Position CT, to generate treatment plan, then keeps treatment plan constant in subsequent therapeutic process, implements some gradation to patient
(Fraction) treatment.This Therapeutic mode does not account for the anatomical structure change of patient in therapeutic process, when patient's solution
When cuing open structure great changes will take place, offline self adaptation treatment (Offline Adaptive Radiation Therapy,
Offline ART) it is usually used to correction treatment plan.However, offline adaptive radiation therapy method is at present using traditional
Commercial therapeutic planning system, its efficiency is low, takes time and effort, so seldom being adopted by hospital.
Online adaptive radiotherapy (Online Adaptive Radiation Therapy, Online ART) is in patient
Certain interval procedure before, obtain patient anatomical structure image, then quickly generate the new treatment plan of patient.Online at present
Self adaptation radiotherapy, it is impossible to be widely used in clinical actual therapeutic, is primarily present following some problems and disadvantages:
1st, online adaptive radiotheraping method is based on common commercial planning system, when generating new treatment plan, needs to count
Calculate and solve large-scale Rapid Dose Calculation and optimization problem, the time that calculates very long it is impossible to complete within the clinical acceptable time
Become.
2nd, the method needs radiation supervisor to participate in the overall process the optimization process of plan, using the setting manually entering line parameter and tune
Whole, take very long, and the quality for the treatment of plan is also directly confined to experience and the judgement of radiation supervisor.
3rd, finally new treatment plan needs again to carry out the quality assurance (Quality Assurance, QA) on die body,
The integration that the optimization of online adaptive radiotherapy planning and automated quality ensure cannot be realized.
Content of the invention
Goal of the invention:In order to overcome the deficiencies in the prior art, the present invention provides a kind of online adaptive radiotherapy meter
Draw automatic optimization method, a kind of online adaptive radiotherapy planning automatic optimizing equipment is provided simultaneously, radiotherapy can be quickly completed
The online adaptive optimization of plan, and complete the quality assurance of dose distribution quadratic search, with the change with reference to anatomical structure simultaneously
Change and adapt to generate the plan being more suitable for each interval procedure, for selection, meet clinical needs well.
Technical scheme:For solving above-mentioned technical problem, the online adaptive radiotherapy planning Automatic Optimal side that the present invention provides
Method, comprises the following steps:
Step 1), read original radiotherapy planning and include launched field information, dosage information and clinical information in interior original letter
Breath;
Step 2), obtain and current include gradation navigational figure, current drawing target outline by several times, currently jeopardize the profile of organ
In interior current gradation information;
Step 3), calculate the beam direction view of each frame angle;
Step 4), change the shape of each frame angle Ziye based on beam direction view, the new life after being optimized
Become Ziye;
Step 5), the dose distribution of each newly-generated Ziye is calculated using convolution superposition algorithm;
Step 6), using the machine jumping figure of conjugate gradient rapid solving and each Ziye of optimization;
Step 7), carry out automatic radiotherapy planning dose distribution quadratic search:Each is calculated again using Monte Carlo EGS4 method
The dose distribution of newly-generated Ziye, and with step 5) dose distribution that obtains carries out three-dimensional gamma index value and calculates assessment, judge
Gamma percent of pass, whether less than default gamma percent of pass threshold value, in this way, then enters step 8), such as no, then enter step 9);
Step 8), to step 7) be normalized computing and update using the calculated dose distribution of Monte Carlo EGS4 method
The machine jumping figure of each Ziye;
Step 9), complete Automatic Optimal:Preserve optimize after newly-generated Ziye, the dose distribution of each newly-generated Ziye and
Machine jumping figure, as new radiotherapy planning, terminates.
As preferred, described step 1) in, described original radiotherapy planning is static intensity modulating radiotherapy plan (Static
IMRT) or arc adjust strong radiotherapy plan (Arc IMRT), wherein
The launched field information of original radiotherapy planning, including:Frame angle, the multi-blade collimator of each Ziye (Segment)
(Multi-Leaf Collimator, MLC) sequence, machine jumping figure (Machine Unit, MU);
The dosage information of original radiotherapy planning, including:Dose distribution (Dose Distribution) distribution and dose volume
Rectangular histogram (Dose Volume Histogram, DVH);
The clinical information of original radiotherapy planning, including:The original drawing target outline of patient and the original profile jeopardizing organ;
Described step 2) in, described gradation navigational figure includes CT (Computed Tomography), Cone-Beam CT (Cone-
Beam CT), ultrasonic (Ultrasound), PET (Positron Emission Tomography) or magnetic resonance (Magnetic
Resonance, MR) etc. at least one in image;
Described step 3) in, comprise the steps:
Step 3.1), by Siddon Ray-Tracing projection algorithm, quickly calculate original drawing target outline in each machine
The beam direction view (BEV) of frame angle;
Step 3.2), by Siddon Ray-Tracing projection algorithm, quickly calculate current drawing target outline in each machine
The beam direction view (BEV) of frame angle;
Described step 4) in, comprise the steps:
Step 4.1), when original radiotherapy planning is static intensity modulating radiotherapy plan (Static IMRT), directly from former
Original scheme Ziye shape is read in beginning project documentation;When original radiotherapy planning is that arc adjusts strong radiotherapy plan (Arc
From original radiotherapy planning, when IMRT), read the multi-blade collimator position at each control point, two companies are obtained by linear interpolation
The original scheme Ziye shape of continuous control point intermediate stand angle;
Step 4.2), by the beam direction view (BEV) of original drawing target outline and the geometry of original scheme Ziye shape
Relation, based on the beam direction view (BEV) of current drawing target outline, according to the beam of current Ziye shape and current drawing target outline
The geometrical relationship of direction view (BEV) and the beam direction view BEV of original drawing target outline and the geometrical relationship one of original Ziye
The criterion causing, obtains current Ziye shape;
Step 4.3), when original radiotherapy planning is static intensity modulating radiotherapy plan (Static IMRT), preserve current
Ziye shape is as the newly-generated Ziye after optimizing;When original radiotherapy planning is that arc adjusts strong radiotherapy plan (Arc IMRT)
When, with the maximum movement speed less than the multi-blade collimator between adjacent Ziye and frame for constraining principle, to current Ziye
Shape changes the newly-generated Ziye after being optimized;
Described step 5) in, dose distribution D of each newly-generated Ziye is calculated using convolution superposition algorithmij, DijFor jth
The dose contribution to i-th volume elements for the individual Ziye;
Described step 6) in, by formula 1.1, solve and optimize the machine jumping figure set of Ziye using conjugate gradient:
So that
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 organ
Set, NTFor the volume elements quantity of each target area, NSVolume elements quantity for each organ;
DijFor j-th Ziye dose contribution to i-th volume elements;
xjThe machine jumping figure of each Ziye for optimizing;
diDose distribution for i-th volume elements;
PiReference dose distribution for i-th volume elements;
WithFor kth time iteration organ weight factor;
Described step 7) in, described dose distribution quadratic search comprises the steps:
Step 7.1), the dose distribution of each newly-generated Ziye is calculated again using Monte Carlo EGS4 method;
Step 7.2), by step 7.1) calculated dose distribution and step 5) dose distribution that obtains carries out three-dimensional gal
Horse index value calculates assessment, judges that gamma percent of pass, whether less than default gamma percent of pass threshold value, in this way, then enters step 8),
As no, then entrance step 9);
Described step 8) in, accept the principle of 100% prescribed dose according to the Gross Target Volume making 95%, to step
7.1) computing is normalized using the calculated dose distribution of Monte Carlo EGS4 method and obtains proportionality factors lambda, by λ and step
6) the machine jumping figure linear multiplication of each Ziye obtaining, updates the machine jumping figure of each Ziye;
Described step 9) in, the step completing Automatic Optimal includes:When original radiotherapy planning is static intensity modulating radiotherapy
During plan (Static IMRT), preserve the newly-generated Ziye after optimizing, the dose distribution of each newly-generated Ziye and machine jumping figure
Generate new radiotherapy planning, terminate;When original radiotherapy planning is that arc adjusts strong radiotherapy plan (Arc IMRT), by new life
The Ziye shape linear interpolation becoming goes back to the position at each control point in original radiotherapy planning, divides with the dosage of each newly-generated Ziye
Cloth and machine jumping figure jointly preserve and generate new radiotherapy planning, terminate.
As further preferred, the step 9 of above-mentioned online adaptive radiotherapy planning automatic optimization method) obtain new
Radiotherapy planning comprises the following steps in transmitting procedure:
Step 10), radiotherapy planning parameter transmitting procedure checks:When new radiotherapy planning is transferred to the checking system of dosage execution
System (Record and Verify System, R&V) after, read (R&V) checking system include multi-blade collimator leaf position,
Machine jumping figure, frame angle, in interior treatment parameter, are compared with the parameter in new radiotherapy planning, judge in checking system
Treatment parameter whether consistent with the parameter in new radiotherapy planning, such as otherwise whether inspection record checking system (R&V) abnormal,
In this way, then preserve new radiotherapy planning as this fractionated radiotherapy plan.
As further, preferably, this fractionated radiotherapy plan described comprises the following steps in the process of implementation:
Step 11), real-time dose measurement:Using electronics portal imaging in the implementation procedure of this fractionated radiotherapy plan described
Multi-blade collimator position in device (Electronic Portal Imaging Device, EPID) record therapeutic process, or make
With daily record (Log) the file real-time verification multi-blade collimator position of accelerator record, and with radiotherapy planning in multi-blade collimator
It is compared, real-time reconstruction 3-dimensional dose dose distribution, if the dosage of real-time reconstruction dose distribution and this fractionated radiotherapy plan
The difference of the distribution prompting message of the output more than 3%, and stop to execute.
Preferably, described step 7.2) in three-dimensional gamma index value calculating parameter be { 3mm, 3% }, default gamma lead to
Crossing rate threshold value is 95% or 96% or 97% or 98% or 99%;Certainly above-mentioned parameter and threshold value can flexibly set as needed
Put.
Present invention simultaneously provides a kind of online adaptive radiotherapy planning automatic optimizing equipment, including:
Device 1), include launched field information, dosage information and clinical information interior former for reading original radiotherapy planning
The device of beginning information;
Device 2), include gradation navigational figure, current drawing target outline, currently jeopardize organ for obtaining current gradation
Profile is in the device of interior current gradation information;
Device 3), for calculating the device of the beam direction view of each frame angle;
Device 4), for changing the shape of each frame angle Ziye based on beam direction view, after being optimized
The device of newly-generated Ziye;
Device 5), for calculating the device of the dose distribution of each newly-generated Ziye using convolution superposition algorithm;
Device 6), for the device using conjugate gradient rapid solving and the machine jumping figure optimizing each Ziye;
Device 7), for carrying out the device of automatic radiotherapy planning dose distribution quadratic search, including:
For calculating the device of the dose distribution of each newly-generated Ziye again using Monte Carlo EGS4 method;
For by using Monte Carlo EGS4 method calculated dose distribution and step 5 again) dose distribution that obtains enters
Row three-dimensional gamma index value calculates assessment, and judges whether gamma percent of pass must assess judgement less than default gamma percent of pass threshold value
Device;
For according to assessment judgment means output valve, select jump to step 8) or step 9) selection redirect device;
Device 8), for step 7) computing is normalized simultaneously using the calculated dose distribution of Monte Carlo EGS4 method
Update the device of the machine jumping figure of each Ziye;
Device 9), for preserving the newly-generated Ziye after optimizing, the dose distribution of each newly-generated Ziye and machine jumping figure
Device as new radiotherapy planning.
Preferably, described device 3) in, including following device:
Device 3.1), by Siddon Ray-Tracing projection algorithm, quickly calculate original drawing target outline in each machine
The device of the beam direction view (BEV) of frame angle;
Device 3.2) pass through Siddon Ray-Tracing projection algorithm, quickly calculate current drawing target outline in each frame
The device of the beam direction view (BEV) of angle;
Described device 4) in, including following device:
Device 4.1), when original radiotherapy planning is static intensity modulating radiotherapy plan (Static IMRT), directly from former
The device of original scheme Ziye shape is read in beginning project documentation;And when original radiotherapy planning is that arc adjusts strong radiotherapy meter
When drawing (Arc IMRT), read the multi-blade collimator position at each control point from original radiotherapy planning file, by linearly inserting
It is worth to the device of the original scheme Ziye shape of two continuous control point intermediate stand angles;
Device 4.2), by the beam direction view (BEV) of original drawing target outline and the geometry of original scheme Ziye shape
Relation, based on the beam direction view (BEV) of current drawing target outline, according to the beam of current Ziye shape and current drawing target outline
The geometrical relationship of direction view (BEV) and the beam direction view (BEV) of original drawing target outline and the geometrical relationship of original Ziye
Consistent criterion, obtains the device of current Ziye shape;
Device 4.3), when original radiotherapy planning is static intensity modulating radiotherapy plan (Static IMRT), preserve current
Ziye shape is as the device of the newly-generated Ziye after optimizing;When original radiotherapy planning is that arc adjusts strong radiotherapy plan (Arc
When IMRT), with the maximum movement speed less than the multi-blade collimator between adjacent Ziye and frame for constraining principle, to current
Ziye shape changes the device of the newly-generated Ziye after being optimized;
Described device 8) be, for accepting the principle of 100% prescribed dose according to the Gross Target Volume making 95%, to step
Rapid 7.1) are normalized, using the calculated dose distribution of Monte Carlo EGS4 method, the device that computing obtains proportionality factors lambda, with
And for by λ and step 6) the machine jumping figure linear multiplication of each Ziye that obtains, update the dress of the machine jumping figure of each Ziye
Put;
Described device 9) in, including when original radiotherapy planning is static intensity modulating radiotherapy plan (Static IMRT),
Preserve the dress that the newly-generated Ziye after optimizing, the dose distribution of each newly-generated Ziye and machine jumping figure generate new radiotherapy planning
Put;And when original radiotherapy planning is that arc adjusts strong radiotherapy plan (Arc IMRT), by newly-generated Ziye shape line
Property interpolation go back to the position at each control point in original radiotherapy planning, with the dose distribution of each newly-generated Ziye and machine jumping figure altogether
With the device preserving the new radiotherapy planning of generation.
As further advantageous embodiment, described online adaptive radiotherapy planning automatic optimizing equipment also includes device 10),
Described device 10) it is radiotherapy planning parameter transmitting procedure check device, it is transferred to dosage execution including when new radiotherapy planning
After checking system (Record and Verify System, R&V), read checking system and include multi-blade collimator blade position
Put, machine jumping figure, frame angle be in interior treatment parameter, the device being compared with the parameter in new radiotherapy planning;And
Judge whether the treatment parameter in checking system is consistent with the parameter in new radiotherapy planning, such as otherwise inspection record checking system
(R&V) whether abnormal, in this way, then preserve new radiotherapy planning as the judgement check device of this fractionated radiotherapy plan.
As further advantageous embodiment, described online adaptive radiotherapy planning automatic optimizing equipment also includes device 11),
Described device 11) it is real-time dose measuring systems, penetrated using electronics including in the implementation procedure of this fractionated radiotherapy plan described
Multi-blade collimator position in wild device for image (Electronic Portal Imaging Device, EPID) record therapeutic process
Put, or using accelerator record daily record (Log) file real-time verification multi-blade collimator position, and with radiotherapy planning in multipage
Collimator is compared, the device of real-time reconstruction 3-dimensional dose dose distribution;And when real-time reconstruction dose distribution and this gradation
The difference of the dose distribution of the radiotherapy planning prompting message of the output more than 3%, and stop the device executing.
Beneficial effect:The present invention provides a kind of online adaptive radiotherapy planning automatic optimization method, provides one kind to exist simultaneously
Line self adaptation radiotherapy planning automatic optimizing equipment, by reading original radiotherapy planning and the current gradation navigational figure obtaining, hook
Draw the information such as target area, calculate the beam direction view of each frame angle, and then obtain the newly-generated Ziye after optimizing and its right
The contents such as the dose distribution answered, machine jumping figure, complete radiotherapy planning online adaptive Automatic Optimal while, by three-dimensional gal
Horse index value calculates the dose distribution quadratic search that assessment carries out automatic radiotherapy planning, and such as quadratic search assessment is not passed through, then right
Dose distribution is normalized computing and updates the machine jumping figure of Ziye, difficulty action accomplishment guarantee work.
Present method invention combines the change of anatomical structure during patient, calculates scheduling algorithm by fast dose and makes
Whole radiotherapy planning optimization process can complete within more than ten minutes or even a few minutes, can be fast before each interval procedure starts
Speed completes online adaptive Automatic Optimal and the Quality Assurance of whole radiotherapy planning, so that therapeutic choice uses.Further
, the present invention can be checked by radiotherapy planning parameter transmitting procedure, is effectively prevented from parameter transmission and the mistake of write execution system
The problem being likely to occur in journey, further, can pass through real-time dose measurement in plan implementation procedure, compare real-time reconstruction agent
Amount distribution and the dose distribution of this fractionated radiotherapy plan carry out risk and security monitoring, once occur security risk remind in time and
Stop plan execution.Achieve online adaptive Automatic Optimal, quality assurance, parameter transfer check and the real time agent of radiotherapy planning
The full-automation of amount detection.
Generally speaking, compared to the offline self adaptation radiotheraping method based on business radiotherapy planning system and existing adaptive online
Answer radiotheraping method, efficiency high of the present invention, save valuable time and human cost, meet clinical required, can clinically promote
Application, has significant social meaning.
Brief description
Fig. 1 is the flow chart of the method invention that embodiment 1 provides.
Specific embodiment
With reference to embodiment, the present invention is described in further detail, and this enforcement row do not constitute restriction to the present invention.
Before original radiotherapy planning is first or before certain interval procedure execution of certain course for the treatment of, based on patient tumors and about
The image information of organ-tissue, through Target delineations, the radiotherapy planning of confirmation.
First or before certain interval procedure execution after therapeutic process in, also will carry out interval procedure to patient,
Before such interval procedure execution, the present invention considers the change of patient anatomy in therapeutic process, to original radiotherapy
Plan is optimized to generate new radiotherapy planning, as this fractionated radiotherapy plan.
Embodiment 1:The online adaptive radiotherapy planning automatic optimization method that the present embodiment 1 is provided, as shown in figure 1, bag
Include following steps:
Step 1), read original radiotherapy planning and include launched field information, dosage information and clinical information in interior original letter
Breath;
Wherein, described original radiotherapy planning is static intensity modulating radiotherapy plan (Static IMRT) or arc is adjusted and put by force
Penetrate radiotherapy planning (Arc IMRT);
The launched field information of original radiotherapy planning includes:Frame angle, the multi-blade collimator sequence of each Ziye, machine is jumped
Number;
The dosage information of original radiotherapy planning includes:Dose distribution and dose volume histogram;
The clinical information of original radiotherapy planning includes:The original drawing target outline of patient and the original profile jeopardizing organ;
Step 2), obtain and current include gradation navigational figure, current drawing target outline by several times, currently jeopardize the profile of organ
In interior current gradation information;Wherein said gradation navigational figure is included in CT, Cone-Beam CT, ultrasonic, PET or magnetic resonance at least
A kind of;
Step 3), calculate the beam direction view of each frame angle, including:
Step 3.1), by Siddon Ray-Tracing projection algorithm, quickly calculate original drawing target outline in each machine
The beam direction view (BEV) of frame angle;
Step 3.2), by Siddon Ray-Tracing projection algorithm, quickly calculate current drawing target outline in each machine
The beam direction view (BEV) of frame angle;
Step 4), change the shape of each frame angle Ziye based on beam direction view, the new life after being optimized
Become Ziye, including:
Step 4.1), when original radiotherapy planning is static intensity modulating radiotherapy plan (Static IMRT), directly from former
Original scheme Ziye shape is read in beginning project documentation;When original radiotherapy planning is that arc adjusts strong radiotherapy plan (Arc
From original radiotherapy planning, when IMRT), read the multi-blade collimator position at each control point, two companies are obtained by linear interpolation
The original scheme Ziye shape of continuous control point intermediate stand angle;
Step 4.2), by the beam direction view (BEV) of original drawing target outline and the geometry of original scheme Ziye shape
Relation, based on the beam direction view (BEV) of current drawing target outline, according to the beam of current Ziye shape and current drawing target outline
The geometrical relationship of direction view (BEV) and the beam direction view (BEV) of original drawing target outline and the geometrical relationship of original Ziye
Consistent criterion, obtains current Ziye shape;
Step 4.3), when original radiotherapy planning is static intensity modulating radiotherapy plan (Static IMRT), preserve current
Ziye shape is as the newly-generated Ziye after optimizing;When original radiotherapy planning is that arc adjusts strong radiotherapy plan (Arc IMRT)
When, with the maximum movement speed less than the multi-blade collimator between adjacent Ziye and frame for constraining principle, to current Ziye
Shape changes the newly-generated Ziye after being optimized;
Step 5), calculate the dose distribution of each newly-generated Ziye using convolution superposition algorithm, specially adopt convolution to fold
Computation system calculates dose distribution D of each newly-generated Ziyeij, DijFor j-th Ziye dose contribution to i-th volume elements;;
Step 6), using conjugate gradient rapid solving and optimize each Ziye machine jumping figure, specially:
By formula 1.1, solve and optimize the machine jumping figure set of Ziye using conjugate gradient:
So that
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, NTFor the volume elements quantity of each target area, NSVolume elements quantity for each organ;
DijFor j-th Ziye dose contribution to i-th volume elements;
xjThe machine jumping figure of each Ziye for optimizing;
diDose distribution for i-th volume elements;
PiReference dose distribution for i-th volume elements;
WithFor kth time iteration organ weight factor;
Above-mentioned solution and optimization process when calculating first be based on given one group of organ weight factor (for 1.0 it is also possible to
It is set to other values, such as 0.5 or 1.5 or 2 or 5 or 10), in subsequent solving-optimizing iterative process, compare each optimizer
Official original doses volume histogram DVHpWith kth time iteration volume histogram DVHkDifference in areasBy formula (1.2) constantly more
Neologism weight factorWith
The iteration of above-mentioned solution and optimization calculating terminates principle and is:As kth iteration dose volume histogram DVHkWith original
Dose volume histogram DVHpDifference in areas is less than during preset area difference limen value (for 0.05 it is also possible to be arranged as required to as other
Value, such as 0.5 or 1.5 or 2 or 5 or 10), terminate this step and optimize calculating, the machine jumping figure collection of the Ziye after the solving-optimizing obtaining
Close { xj};
Step 7), carry out automatic radiotherapy planning dose distribution quadratic search, including:
Step 7.1), the dose distribution of each newly-generated Ziye is calculated again using Monte Carlo EGS4 method;
Step 7.2), by step 7.1) calculated dose distribution and step 5) dose distribution that obtains carries out three-dimensional gal
Horse index value calculates assessment, judges that gamma percent of pass, whether less than default gamma percent of pass threshold value, in this way, then enters step 8),
As no, then entrance step 9);
Step 8), to step 7) be normalized computing and update using the calculated dose distribution of Monte Carlo EGS4 method
The machine jumping figure of each Ziye, that is,:Accept the principle of 100% prescribed dose according to the Gross Target Volume making 95%, to step
7.1) computing is normalized using the calculated dose distribution of Monte Carlo EGS4 method and obtains proportionality factors lambda, with step 6)
The machine jumping figure of each Ziye arriving and λ linear multiplication, the machine jumping figure of each Ziye after being updated;
Step 9), preserve the newly-generated Ziye after optimizing, the dose distribution of each newly-generated Ziye and machine jumping figure conduct
New radiotherapy planning, completes Automatic Optimal, including:When original radiotherapy planning is static intensity modulating radiotherapy plan (Static
When IMRT), preserve the newly-generated Ziye after optimizing, the dose distribution of each newly-generated Ziye and machine jumping figure and generate new radiotherapy
Plan, terminates;When original radiotherapy planning is that arc adjusts strong radiotherapy plan (Arc IMRT), by newly-generated Ziye shape
Linear interpolation goes back to the position at each control point in original radiotherapy planning, dose distribution and the machine jumping figure with each newly-generated Ziye
Common preservation generates new radiotherapy planning, terminates.
In the present embodiment 1, above-mentioned steps 7.2) in the calculating parameter of three-dimensional gamma index value be { 3mm, 3% }, preset
Gamma percent of pass threshold value is 95%.
The online adaptive radiotherapy planning automatic optimizing equipment that embodiment 1 provides, including:
Device 1), include launched field information, dosage information and clinical information interior former for reading original radiotherapy planning
The device of beginning information;
Device 2), include gradation navigational figure, current drawing target outline, currently jeopardize organ for obtaining current gradation
Profile is in the device of interior current gradation information;
Device 3), for calculating the device of the beam direction view of each frame angle, including:
Device 3.1), by Siddon Ray-Tracing projection algorithm, quickly calculate original drawing target outline in each machine
The device of the beam direction view (BEV) of frame angle;
Device 3.2) pass through Siddon Ray-Tracing projection algorithm, quickly calculate current drawing target outline in each frame
The device of the beam direction view (BEV) of angle;
Device 4), for changing the shape of each frame angle Ziye based on beam direction view, after being optimized
The device of newly-generated Ziye, including:
Device 4.1), when original radiotherapy planning is static intensity modulating radiotherapy plan (Static IMRT), directly from former
The device of original scheme Ziye shape is read in beginning project documentation;And when original radiotherapy planning is that arc adjusts strong radiotherapy meter
When drawing (Arc IMRT), read the multi-blade collimator position at each control point from original radiotherapy planning file, by linearly inserting
It is worth to the device of the original scheme Ziye shape of two continuous control point intermediate stand angles;
Device 4.2), by the beam direction view (BEV) of original drawing target outline and the geometry of original scheme Ziye shape
Relation, based on the beam direction view (BEV) of current drawing target outline, according to the beam of current Ziye shape and current drawing target outline
The geometrical relationship of direction view (BEV) and the beam direction view BEV of original drawing target outline and the geometrical relationship one of original Ziye
The criterion causing, obtains the device of current Ziye shape;
Device 4.3), when original radiotherapy planning is static intensity modulating radiotherapy plan (Static IMRT), preserve current
Ziye shape is as the device of the newly-generated Ziye after optimizing;When original radiotherapy planning is that arc adjusts strong radiotherapy plan (Arc
When IMRT), with the maximum movement speed less than the multi-blade collimator between adjacent Ziye and frame for constraining principle, to current
Ziye shape changes the device of the newly-generated Ziye after being optimized;
Device 5), for calculating the device of the dose distribution of each newly-generated Ziye using convolution superposition algorithm;
Device 6), for the device using conjugate gradient rapid solving and the machine jumping figure optimizing each Ziye;
Device 7), for carrying out the device of automatic radiotherapy planning dose distribution quadratic search, including:
For calculating the device of the dose distribution of each newly-generated Ziye again using Monte Carlo EGS4 method;
For by using Monte Carlo EGS4 method calculated dose distribution and step 5 again) dose distribution that obtains enters
Row three-dimensional gamma index value calculates assessment, and judges whether gamma percent of pass must assess judgement less than default gamma percent of pass threshold value
Device;
For according to assessment judgment means output valve, select jump to step 8) or step 9) selection redirect device;
Device 8), for step 7) computing is normalized simultaneously using the calculated dose distribution of Monte Carlo EGS4 method
Update the device of the machine jumping figure of each Ziye, including:
For accepting the principle of 100% prescribed dose according to the Gross Target Volume making 95%, to step 7.1) using illiteracy
The calculated dose distribution of special Caro algorithm is normalized the device that computing obtains proportionality factors lambda, and for by λ and step
The machine jumping figure linear multiplication of each Ziye that rapid 6) obtain, updates the device of the machine jumping figure of each Ziye;
Device 9), for preserving the newly-generated Ziye after optimizing, the dose distribution of each newly-generated Ziye and machine jumping figure
As the device of new radiotherapy planning, including:Including being static intensity modulating radiotherapy plan (Static when original radiotherapy planning
When IMRT), preserve the newly-generated Ziye after optimizing, the dose distribution of each newly-generated Ziye and machine jumping figure and generate new radiotherapy
The device of plan;And when original radiotherapy planning is that arc adjusts strong radiotherapy plan (Arc IMRT), by newly-generated son
Wild shape linear interpolation goes back to the position at each control point in original radiotherapy planning, dose distribution and the machine with each newly-generated Ziye
Device jumping figure preserves the device generating new radiotherapy planning jointly.
Embodiment 2:Embodiment 2 provide online adaptive radiotherapy planning automatic optimization method and with embodiment 1 provide
Method is basically identical, and something in common is not repeated, except that:
Described step 9) the new radiotherapy planning that obtains is further comprising the steps of in transmitting procedure:
Step 10), radiotherapy planning parameter transmitting procedure checks:When new radiotherapy planning is transferred to the checking system of dosage execution
After system (Record and Verify System, R&V), read checking system and include multi-blade collimator leaf position, machine
Jumping figure, frame angle, in interior treatment parameter, are compared with the parameter in new radiotherapy planning, judge controlling in checking system
Whether consistent with the parameter in new radiotherapy planning treat parameter, such as otherwise whether inspection record checking system (R&V) is abnormal, such as
It is then to preserve new radiotherapy planning as this fractionated radiotherapy plan.
This fractionated radiotherapy plan described comprises the following steps in commission:
Step 11), real-time dose measurement:Using electronics portal imaging in the implementation procedure of this fractionated radiotherapy plan described
Multi-blade collimator position in device (Electronic Portal Imaging Device, EPID) record therapeutic process, or make
With daily record (Log) the file real-time verification multi-blade collimator position of accelerator record, and with radiotherapy planning in multi-blade collimator
It is compared, real-time reconstruction 3-dimensional dose dose distribution, if the dosage of real-time reconstruction dose distribution and this fractionated radiotherapy plan
The difference of the distribution prompting message of the output more than 3%, and stop to execute.
Described step 7.2) in default gamma percent of pass threshold value be 98%.
Embodiment 2 provide online adaptive radiotherapy planning automatic optimizing equipment and with embodiment 1 provide device basic
Unanimously, something in common is not repeated, except that:
Also include device 10), described device 10) it is radiotherapy planning parameter transmitting procedure check device, put including when new
After treating the checking system (Record and Verify System, R&V) that planned transmission executes to dosage, read in checking system
Treatment parameter including multi-blade collimator leaf position, machine jumping figure, frame angle, with the parameter in new radiotherapy planning
The device being compared;And judge whether the treatment parameter in checking system is consistent with the parameter in new radiotherapy planning, such as
Otherwise whether inspection record checking system (R&V) is abnormal, in this way, then preserves new radiotherapy planning as this fractionated radiotherapy plan
Judge check device.
Also include device 11), described device 11) it is real-time dose measuring systems, including in this fractionated radiotherapy plan described
Implementation procedure in using electronic portal image device (Electronic Portal Imaging Device, EPID) record control
Multi-blade collimator position during treatment, or daily record (Log) the file real-time verification multi-blade collimator position using accelerator record,
And be compared with the multi-blade collimator in radiotherapy planning, the device of real-time reconstruction 3-dimensional dose dose distribution;And when in real time
Rebuild the difference prompting message of the output more than 3% of the dose distribution of dose distribution and this fractionated radiotherapy plan, and stop to execute
Device.
Above row of implementing do not constitute restriction to the present invention, and relevant staff is in the scope without departing from the technology of the present invention thought
Interior, carried out various change and modifications, all fall within protection scope of the present invention.
Claims (3)
1. a kind of online adaptive radiotherapy planning automatic optimizing equipment is it is characterised in that include:
Device 1), include launched field information, dosage information and clinical information in interior original letter for reading original radiotherapy planning
The device of breath;
Device 2), current include gradation navigational figure, current drawing target outline by several times, currently jeopardize the profile of organ for obtaining
Device in interior current gradation information;
Device 3), for calculating the device of the beam direction view of each frame angle;
Device 4), for changing the shape of each frame angle Ziye, the new life after being optimized based on beam direction view
Become the device of Ziye;
Device 5), for calculating the device of the dose distribution of each newly-generated Ziye using convolution superposition algorithm;
Device 6), for the device using conjugate gradient rapid solving and the machine jumping figure optimizing each Ziye;
Device 7), for carrying out the device of automatic radiotherapy planning dose distribution quadratic search, including:
For calculating the device of the dose distribution of each newly-generated Ziye again using Monte Carlo EGS4 method;
For by using Monte Carlo EGS4 method calculated dose distribution and device 5 again) dose distribution that obtains carries out three
Dimension gamma index value calculates assessment, and judges whether gamma percent of pass must assess judgement dress less than default gamma percent of pass threshold value
Put;
For according to assessment judgment means output valve, select jump to device 8) or device 9) selection redirect device;
Device 8), for device 7) be normalized computing and update using the calculated dose distribution of Monte Carlo EGS4 method
The device of the machine jumping figure of each Ziye;
Device 9), for preserving the newly-generated Ziye after optimizing, the dose distribution of each newly-generated Ziye and machine jumping figure conduct
The device of new radiotherapy planning;
Wherein:Described device 3) in, including following device:
Device 3.1), by Siddon Ray-Tracing projection algorithm, quickly calculate original drawing target outline at each frame angle
The device of the beam direction view of degree;
Device 3.2) pass through Siddon Ray-Tracing projection algorithm, quickly calculate current drawing target outline in each frame angle
Beam direction view device;
Described device 4) in, including following device:
Device 4.1), when original radiotherapy planning is static intensity modulating radiotherapy plan, directly read from original scheme file
The device of original scheme Ziye shape;And when original radiotherapy planning is that arc adjusts strong radiotherapy plan, from original radiotherapy
Read the multi-blade collimator position at each control point in project documentation, machine in the middle of two continuous control points is obtained by linear interpolation
The device of the original scheme Ziye shape of frame angle;
Device 4.2), by the beam direction view of original drawing target outline and the geometrical relationship of original scheme Ziye shape, it is based on
The beam direction view of current drawing target outline, according to the geometry of current Ziye shape and the beam direction view of current drawing target outline
The relation criterion consistent with the beam direction view of original drawing target outline and the geometrical relationship of original Ziye, obtains current Ziye shape
The device of shape;
Device 4.3), when original radiotherapy planning is static intensity modulating radiotherapy plan, preserve current Ziye shape as optimization
The device of newly-generated Ziye afterwards;When original radiotherapy planning is that arc adjusts strong radiotherapy plan, with less than adjacent Ziye
Between multi-blade collimator and frame maximum movement speed be constraint principle, to current Ziye shape modification optimized after
The device of newly-generated Ziye;
Described device 8) be, for accepting the principle of 100% prescribed dose according to the Gross Target Volume making 95%, to step
7.1) device that computing obtains proportionality factors lambda is normalized using the calculated dose distribution of Monte Carlo EGS4 method, and
For by λ and step 6) the machine jumping figure linear multiplication of each Ziye that obtains, update the device of the machine jumping figure of each Ziye;
Described device 9) in, including when original radiotherapy planning is static intensity modulating radiotherapy plan, preserve the new life after optimizing
Ziye, the dose distribution of each newly-generated Ziye and machine jumping figure is become to generate the device of new radiotherapy planning;And ought original put
When treatment plan is that arc adjusts strong radiotherapy plan, newly-generated Ziye shape linear interpolation is returned each in original radiotherapy planning
The position at control point, preserves, with the dose distribution of each newly-generated Ziye and machine jumping figure, the dress generating new radiotherapy planning jointly
Put.
2. online adaptive radiotherapy planning automatic optimizing equipment according to claim 1 it is characterised in that:Also include device
10), described device 10) it is radiotherapy planning parameter transmitting procedure check device, hold including when new radiotherapy planning is transferred to dosage
After the checking system of row, read that checking system includes multi-blade collimator leaf position, machine jumping figure, frame angle are controlled in interior
Treat the device that the parameter in parameter, with new radiotherapy planning is compared;And judge treatment parameter in checking system whether
Consistent with the parameter in new radiotherapy planning, such as otherwise whether inspection record checking system is abnormal, in this way, then preserves new radiotherapy
It is intended to be the judgement check device of this fractionated radiotherapy plan.
3. online adaptive radiotherapy planning automatic optimizing equipment according to claim 2 it is characterised in that:Also include device
11), described device 11) be real-time dose measuring systems, including in the implementation procedure of this fractionated radiotherapy plan described using electricity
Multi-blade collimator position in sub- portal imaging device record therapeutic process, or the journal file real-time verification using accelerator record
Multi-blade collimator position, and be compared with the multi-blade collimator in radiotherapy planning, real-time reconstruction 3-dimensional dose dose distribution
Device;And when difference output more than 3% of real-time reconstruction dose distribution and the dose distribution of this fractionated radiotherapy plan is reminded
Information, and stop the device executing.
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 CN104338240A (en) | 2015-02-11 |
CN104338240B true 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) |
Families Citing this family (31)
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 |
CN104815392B (en) * | 2015-04-10 | 2018-01-26 | 石峰 | A kind of interactive radiotherapy treatment planning system optimization system |
CN104933652A (en) * | 2015-04-27 | 2015-09-23 | 苏州敏宇医疗科技有限公司 | Cloud-computing based dose verification system and method of tumor radiotherapy |
CN105031819B (en) * | 2015-08-25 | 2018-11-06 | 上海联影医疗科技有限公司 | A kind of injectivity optimizing system |
CN105999567B (en) * | 2016-06-22 | 2018-09-18 | 沈阳东软医疗系统有限公司 | A kind of the leaf position control method and device of electric multiblade collimator |
CN106846317B (en) * | 2017-02-27 | 2021-09-17 | 北京连心医疗科技有限公司 | Medical image retrieval method based on feature extraction and similarity matching |
CN106920234B (en) * | 2017-02-27 | 2021-08-27 | 北京连心医疗科技有限公司 | Combined automatic radiotherapy planning method |
CN106932810B (en) * | 2017-04-01 | 2018-02-23 | 西安一体医疗科技有限公司 | 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 |
CN107510898B (en) * | 2017-10-11 | 2021-04-16 | 上海联影医疗科技股份有限公司 | Light field and radiation field verification device and verification method for medical equipment |
CN108175953B (en) * | 2017-12-29 | 2020-07-21 | 苏州雷泰医疗科技有限公司 | Dynamic grating-based intensity modulation optimization method and radiotherapy equipment |
CN109999365B (en) * | 2018-01-05 | 2021-04-20 | 北京连心医疗科技有限公司 | Computing device and storage medium for implementing radiotherapy reverse optimization method |
CN108434612B (en) * | 2018-01-05 | 2020-07-10 | 东软医疗系统股份有限公司 | Quality assurance testing method and device |
CN110246562B (en) * | 2018-03-07 | 2021-08-06 | 有方(合肥)医疗科技有限公司 | Method, device and computer system for determining 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 |
CN108648811A (en) * | 2018-05-11 | 2018-10-12 | 南方医科大学 | A kind of automatic radiotherapy planning quality evaluating method based on the prediction of indication item |
CN109432611B (en) * | 2018-10-18 | 2021-02-12 | 北京连心医疗科技有限公司 | Method, equipment and storage medium for generating control point of extended dynamic wedge-shaped board |
CN109513119B (en) * | 2018-11-12 | 2020-07-17 | 武汉大学 | Volume rotation intensity modulation radiotherapy plan collimator angle optimization method |
CN109513121B (en) * | 2018-12-28 | 2021-01-01 | 安徽大学 | Dose-guided adaptive radiotherapy plan re-optimization system and method |
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 |
CN112618967B (en) * | 2020-12-17 | 2024-03-15 | 程明霞 | Device for adjusting radiotherapy dosage distribution and generating corresponding radiotherapy plan |
CN112926156B (en) * | 2021-02-25 | 2024-04-26 | 中科超精(南京)科技有限公司 | Direct sub-field optimization device and method based on position variable penumbra model |
CN113827877B (en) * | 2021-09-15 | 2024-05-28 | 上海市胸科医院 | Method for realizing automatic dose verification based on accelerator log file |
CN114010963B (en) * | 2021-11-08 | 2023-10-20 | 上海联影医疗科技股份有限公司 | System, method and computer readable storage medium for dose determination |
CN118382395A (en) * | 2021-12-20 | 2024-07-23 | 西安大医集团股份有限公司 | Radiotherapy equipment control method, radiotherapy equipment control device, radiotherapy equipment control equipment and storage medium |
CN115966281B (en) * | 2023-01-16 | 2024-02-09 | 迈胜医疗设备有限公司 | Method, device, equipment and storage medium for generating arc radiotherapy plan |
CN116328214B (en) * | 2023-05-30 | 2023-08-22 | 福建自贸试验区厦门片区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 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102306239A (en) * | 2011-07-22 | 2012-01-04 | 李宝生 | Method for evaluating and optimizing radiotherapy dose based on cone beam CT (Computer Tomography) image CT value correction technology |
CN104117151A (en) * | 2014-08-12 | 2014-10-29 | 章桦 | Optimization method of online self-adaption radiotherapy plan |
-
2014
- 2014-10-31 CN CN201410607268.1A patent/CN104338240B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102306239A (en) * | 2011-07-22 | 2012-01-04 | 李宝生 | Method for evaluating and optimizing radiotherapy dose based on cone beam CT (Computer Tomography) image CT value correction technology |
CN104117151A (en) * | 2014-08-12 | 2014-10-29 | 章桦 | Optimization method of online self-adaption radiotherapy plan |
Non-Patent Citations (1)
Title |
---|
An on-line replanning scheme for interfractional variations;Ahunbay.E.E, et al;《Medical Physics》;20080731;第35卷(第8期);第3608-3615页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104338240A (en) | 2015-02-11 |
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 | |
CN104548372B (en) | The dosage determining device of radiotherapy | |
US8222616B2 (en) | Method for adapting fractionation of a radiation therapy dose | |
Sun et al. | Evaluation of the efficiency and effectiveness of independent dose calculation followed by machine log file analysis against conventional measurement based IMRT QA | |
US7639854B2 (en) | Method and system for processing data relating to a radiation therapy treatment plan | |
Bossuyt et al. | Evaluation of automated pre-treatment and transit in-vivo dosimetry in radiotherapy using empirically determined parameters | |
CN112581475B (en) | Method for predicting gamma passing rate of radiotherapy plan | |
Cai et al. | A practical implementation of physics quality assurance for photon adaptive radiotherapy | |
US11043142B2 (en) | Virtual radiation oncology clinic and methods of use thereof | |
CN110415785A (en) | The method and system of artificial intelligence guidance radiotherapy planning | |
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 | |
US20220296925A1 (en) | Method and system for robust radiotherapy treatment planning for biological uncertainties | |
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 | |
Miften et al. | A dose‐volume‐based tool for evaluating and ranking IMRT treatment plans | |
Cilla et al. | Epid‐based in vivo dose verification for lung stereotactic treatments delivered with multiple breath‐hold segmented volumetric modulated arc therapy | |
Newpower et al. | Spot delivery error predictions for intensity modulated proton therapy using robustness analysis with machine learning | |
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 | |
Kisling | Development of automated radiotherapy treatment planning for cervical and breast cancer for resource-constrained clinics | |
Mijnheer | Patient-Specific Quality Assurance: In Vivo 3D Dose Verification | |
Kostovski | Artificial Intelligence in Radiotherapy | |
Gueorguiev et al. | Clinical evaluation of a novel transmission detector for 3D quality assurance of IMRT and SBRT | |
Satherley | The impact of plan complexity on the accuracy of VMAT for the treatment of head and neck cancer | |
Conibear | Assessment of Target Volume and Organ at Risk Contouring Variability within the Context of UK Head and Neck and Lung Cancer Radiotherapy Clinical Trials |
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 |