CN105866821A - Method and device for obtaining beam energy spectrum, and dose distribution calculating method - Google Patents
Method and device for obtaining beam energy spectrum, and dose distribution calculating method Download PDFInfo
- Publication number
- CN105866821A CN105866821A CN201610176331.XA CN201610176331A CN105866821A CN 105866821 A CN105866821 A CN 105866821A CN 201610176331 A CN201610176331 A CN 201610176331A CN 105866821 A CN105866821 A CN 105866821A
- Authority
- CN
- China
- Prior art keywords
- dose
- distribution
- value
- depth
- spectral distribution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Measurement Of Radiation (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
The invention discloses a method for obtaining a beam energy spectrum, and the method comprises the steps: setting the initial energy spectrum distribution of a beam; calculating the dose distribution on a central axis of the beam based on the initial energy spectrum distribution; judging whether a calculated depth value of a maximum dose place and the dose value of the maximum dose place are matched with the measurement value of dose distribution on the central axis of the beam or not: enabling the current energy spectrum distribution to serve as the energy spectrum distribution of the beam if the calculated depth value and the dose value are matched with the measurement value, or else, adjusting the current energy spectrum distribution, and calculating the dose distribution on the central axis of the beam based on the adjusted energy spectrum distribution till the calculated depth value and the dose value are matched with the measurement value. The method is simple and effective.
Description
Technical field
The present invention relates to field of medical device, a kind of method particularly relating to beam power spectrum that can obtain medical electronic linear accelerator.
Background technology
It is known that radiotherapy is one of most important, most basic means of oncotherapy.In order to improve the cure rate of tumor, it is desirable to the beam of high dose accurately to be rendered to target area (tumor), reduce the injury to normal surrounding tissue as far as possible, therefore, need accurately to calculate the dose distribution of beam when formulating radiotherapy treatment planning.
Photon beam dose calculation methodology (also referred to as Collapsed Cone dose calculation methodology based on convolution/Additive Model, it is called for short CC algorithm) dose distribution of degree of precision can be obtained, but it needs the beam power spectrum produced according to medical electronic linear accelerator to calculate, the only spectral information of beam is the most accurate, the most calculated dose distribution is only possible to enough accurate, and the beam power spectrum therefore obtaining medical electronic linear accelerator is the key of collapsed cone dose calculation methodology.
But, the beam energy that accelerator treatment head produces is high, and its power spectrum is difficult to directly measure.Prior art 1:CN101071172A discloses a kind of method measuring medical accelerator power spectrum, solves system of linear equations by relaxative iteration and obtains the power spectrum of accelerator beam.But, its reconstruction model does not comprise the constraints of power spectrum, and the solution that the method is obtained is the most unstable, the biggest with actual power spectrum deviation.Prior art 2:CN103336295A discloses a kind of photon beam energy spectrum of medical electric linear accelerator acquisition methods, utilize the method for Monte Carlo to calculate these Single photons dose distribution in central depths by the model of a series of Single photons, and use optimization method to be fitted obtaining a spectral distribution as initial condition.The measurement data but this power spectrum obtained by covering card simulation can not be coincide in CC algorithm well, such as intersect with experiment curv certain position in declining district when value of calculation curve, at this moment, in the case of uncomfortable energy-conservation spectrum, regulating other parameter always can not make the two ends data coincidence of intersection point all improve, it is common that one end improves and the other end degenerates.Therefore, the power spectrum obtained by the method is the most accurate.
Accordingly, it would be desirable to a kind of method seeking new acquisition medical electronic linear accelerator power spectrum.
Summary of the invention
In order to solve above-mentioned technical problem, the invention provides a kind of new method obtaining beam power spectrum, including:
The initial spectral distribution of beam is set；
Based on initial spectral distribution, calculate the dose distribution on beam centre axle；
Judge that the depth value at calculated maximal dose and the dose value at the depth capacity measured value the most all with the dose distribution on beam centre axle coincide；
If coincideing, the most current spectral distribution is the spectral distribution of beam；
If misfitting, then regulate current spectral distribution, and calculate the dose distribution on described beam centre axle based on the spectral distribution after regulation, till the depth value at calculated maximal dose and the dose value at depth capacity all coincide with measured value.
Optionally, described initial spectral distribution is Unimodal Distribution.
Optionally, the basic attenuation quotient of beam is set and with the attenuation quotient of change in depth, calculates the dose distribution on described beam centre axle by CC algorithm.
Optionally, the measured value of the dose distribution on described beam centre axle is obtained by 3 d water tank system or film measurement.
Optionally, before judging that the measured value whether with the dose distribution on beam centre axle of the dose value at calculated depth capacity coincide, to the dose distribution on calculated beam centre axle and the dose distribution normalization respectively measured on the beam centre axle obtained.
Optionally, regulation energy time spectrum, power spectrum meets described Unimodal Distribution all the time.
Optionally, regulation energy time spectrum, the multiple adjacent energy ingredient of adjusted in concert.
Optionally, the spectral distribution that described regulation is current, including,
If the depth value at calculated maximal dose is more than the depth value at maximal dose in measured value, then increases low-energy component and/or reduce high-energy component, otherwise, then reduce low-energy component and/or increase high-energy component；
If the dose value at calculated depth capacity is more than the dose value at depth capacity in measured value, then increase low-energy component and/or reduce high-energy component, otherwise, then reduce low-energy component and/or increase high-energy component, wherein, described high-energy component represents at least one energy more than peak energy, and described low-energy component represents at least one energy less than peak energy.
Optionally, also include regulating described basis attenuation quotient and the attenuation quotient with change in depth so that in the dose value at calculated depth capacity and measured value, the dose value at depth capacity coincide.
Accordingly, the invention provides a kind of device obtaining beam power spectrum, including,
Input block, for arranging the initial spectral distribution of beam, described initial spectral distribution is Unimodal Distribution；
Computing unit, for spectral distribution based on beam, calculates the dose distribution on beam centre axle by CC algorithm；
Judging unit, for judging that the depth value at calculated maximal dose and the dose value at the depth capacity measured value the most all with the dose distribution on beam centre axle coincide；
Regulation unit, for regulating current spectral distribution；
If the spectral distribution that judging unit output result is identical, the most current is the spectral distribution of beam；
If judging unit output result is for misfitting, the spectral distribution that then regulation of regulation unit is current, and the spectral distribution after regulation is input to computing unit, computing unit calculates the dose distribution on described beam centre axle based on the spectral distribution after regulation, till the depth value at calculated maximal dose and the dose value at depth capacity all coincide with measured value.
Optionally, described regulation unit also regulates the basic attenuation quotient in CC algorithm and the attenuation quotient with change in depth so that in the dose value at calculated depth capacity and measured value, the dose value at depth capacity coincide.
Present invention also offers the another kind of method obtaining beam power spectrum, including:
Arranging the initial spectral distribution of beam, described initial spectral distribution is Unimodal Distribution；
Based on described initial spectral distribution, calculate the dose distribution on beam centre axle by CC algorithm；
Judge whether the depth value at calculated maximal dose coincide with the depth value measured at the maximal dose obtained；
If coincideing, the most current spectral distribution is the spectral distribution of beam；
If misfitting, then regulate current spectral distribution, and calculate the dose distribution on described beam centre axle based on the spectral distribution after regulation, till the depth value at calculated maximal dose coincide with the depth value measured at the maximal dose obtained.
Present invention also offers a kind of dose distributions computation method, including,
The method utilizing aforementioned acquisition beam power spectrum obtains the spectral distribution of beam；
The off-axis ratio of regulation beam, updates the dose distribution of described beam by CC algorithm, until the dose distribution after Geng Xining is coincide with measuring the dose distribution obtained, thus the dose distribution after updating is the dose distribution of calculated beam.
Optionally, the regulation output factor is also included so that the dose distribution of beam and measured value coincide on absolute value.
Optionally, the dose distribution that described measurement obtains is obtained by 3 d water tank system or film measurement.
The present invention also provides for a kind of dose distributions computation method, including:
The method utilizing aforementioned acquisition beam power spectrum obtains the spectral distribution of beam；
The basic attenuation quotient of regulation beam and with the attenuation quotient of change in depth, updates described beam dose distribution on center shaft by CC algorithm, until the dose distribution after Geng Xining is identical with the measured value of the dose distribution on beam centre axle；
The off-axis ratio of regulation beam, by CC algorithm update described beam with the dose distribution in beam centre axle vertical direction, until the measured value of the dose distribution in the dose distribution after Geng Xining and beam centre axle vertical direction is identical；
Dose distribution after the renewal that three-dimensional all coincide with measured value is the calculated beam dose distribution at three-dimensional.
Compared to prior art, the method obtaining beam power spectrum of the present invention is based on initial spectral distribution, it is calculated dose distribution by CC algorithm, by the dose distribution of calculating is compared with measuring the dose distribution obtained, spectral distribution is regulated according to result of the comparison, therefore the measurement of power spectrum is converted into the regulation of parameter by the present invention, conveniently carries out automation mechanized operation, and computational efficiency is higher；Further, by reducing the step-length of parameter, the precision of power spectrum can be improved.
Further, first the present invention regulates spectral distribution and the depth value at maximal dose and the dose value at depth capacity is all matched with measured value, regulating other parameter again makes the dose distribution on the dosage of other depth point and other direction all match with measured value, and the method for this calculating power spectrum is simply, easily realize.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, the accompanying drawing used required in embodiment or description of the prior art will be briefly described below, apparently, accompanying drawing in describing below is only some embodiments of the present invention, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to obtain other accompanying drawing according to these accompanying drawings.
Fig. 1 is the method flow diagram obtaining beam power spectrum of the embodiment of the present invention；
Fig. 2 is the initial spectral distribution of the embodiment of the present invention；
Fig. 3 is the regulation schematic diagram reducing Low Energy Region energy weight of the embodiment of the present invention；
Fig. 4 is reducing Low Energy Region energy weight and increasing the regulation schematic diagram of high energy region energy weight of the embodiment of the present invention；
Fig. 5 is the regulation schematic diagram increasing high energy region energy weight of the embodiment of the present invention；
Fig. 6 is the method flow diagram obtaining beam power spectrum of another embodiment of the present invention；
Fig. 7 is the initial spectral distribution of 6MV accelerator beam；
Fig. 8 is the spectral distribution after the regulation of 6MV accelerator beam；
Fig. 9 is the PDD curve comparison diagram with experiment curv of 6MV accelerator beam.
Figure 10 is the device schematic diagram obtaining beam power spectrum of the embodiment of the present invention.
Detailed description of the invention
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments.Based on the embodiment in the present invention, the every other embodiment that those of ordinary skill in the art are obtained on the premise of not making creative work, broadly fall into the scope of protection of the invention.
Beam for medical electronic linear accelerator outgoing, its be by the electronics of certain energy beat treatment head in target on produce bremsstrahlung light molecular, and the power spectrum of bremsstrahlung photon is the widest, have from mental retardation to high energy, if it is not known that the spectral distribution of photon, cannot accurately determine beam dose distribution in human body.
The invention provides a kind of method obtaining beam power spectrum, it is calculated dose distribution based on initial spectral distribution, by the dose distribution of calculating is compared with measuring the dose distribution obtained, spectral distribution is regulated according to result of the comparison, thus obtain the power spectrum being consistent with actual power spectrum, the most calculated dose distribution precision is higher, meets the dose distribution of reality.
Fig. 1 is the method flow diagram obtaining beam power spectrum of the embodiment of the present invention；Fig. 2 is the initial spectral distribution of the embodiment of the present invention；Fig. 3 is the regulation schematic diagram reducing Low Energy Region energy weight of the embodiment of the present invention；Fig. 4 is reducing Low Energy Region energy weight and increasing the regulation schematic diagram of high energy region energy weight of the embodiment of the present invention；Fig. 5 is the regulation schematic diagram increasing high energy region energy weight of the embodiment of the present invention.
Described in reference diagram 1, the method for the present embodiment comprises the steps:
In step 101, the initial spectral distribution of beam is set.
The beam that electronics occurs bremsstrahlung to produce in target has continuum, and is Unimodal Distribution, and the initial power spectrum being therefore preferably provided with beam is Unimodal Distribution, such as Gauss distribution.
Calculating for convenience, be finely divided by the energy of beam, can divide equally, it is also possible to do not divide equally, the average energy taking each interval represents the photon energy in this interval, and the mono-energetic beams of each energy weight shared by beam forms the power spectrum of this beam.Such as, the energy of beam is divided into n interval, mental retardation to high energy is designated as E successively0, E1..., En.With reference to shown in Fig. 2, having 15 energy, in figure, the straight line along the vertical direction of transverse axis distribution represents each ENERGY E0, E1..., E14, the longitudinal axis represents the weight that each energy is corresponding, i.e. from left to right, and the y value that the straight line of the 1st article of vertical direction is corresponding with the intersection point of Gaussian function curve is ENERGY E0Weight, the y value that the straight line of the 2nd article of vertical direction is corresponding with the intersection point of Gaussian function curve is ENERGY E1Weight, by that analogy, the y value that the straight line of the 15th article of vertical direction is corresponding with the intersection point of Gaussian function curve is ENERGY E14Weight.
Initial spectral distribution can be rule of thumb set, if actual beam is the hardest, the high energy region that the peak energy of initial power spectrum is positioned in power spectrum is then set, if actual beam is the softest, the Low Energy Region that the peak energy of initial power spectrum is positioned in power spectrum is then set, such that it is able to quickening convergence rate, comparatively fast obtain accurate power spectrum.In the present embodiment, the peak energy of initial power spectrum is positioned at the ENERGY E in high energy region more to the right in power spectrum, i.e. Fig. 210Place.
In the present embodiment, if the energy of beam being split the thinnest, when i.e. n is bigger, it is believed that photon energy approximately equal in each interval, in this case, the information that the initial spectral distribution of setting comprises is the most, can improve the precision of final calculated power spectrum, and with actual power spectrum closer to.
In step 102, based on initial spectral distribution, calculate the dose distribution on beam centre axle.
In the present embodiment, CC algorithm is used to calculate the dose distribution on beam centre axle.The basic attenuation quotient of beam is set and with the attenuation quotient of change in depth, based on initial spectral distribution, calculates the dose distribution on described beam centre axle by CC algorithm.
Two kinds of convolution kernels can be used in CC algorithm to carry out Rapid Dose Calculation, two kinds of the most corresponding maximum launched fields of convolution kernel and the situation of minimum launched field.Described basis attenuation quotient is used for determining under current launched field and uses which kind of convolution kernel to carry out Rapid Dose Calculation, in the case of other parameter constants, basis attenuation quotient is the biggest, the convolution kernel used under current launched field is closer to the convolution kernel of maximum launched field, then ray is relatively soft, basis attenuation quotient is the least, and the convolution kernel used under current launched field is closer to the convolution kernel of minimum launched field, then ray is relatively hard.The described attenuation quotient with change in depth represents when ray transmits in die body, the decay of its energy with its in die body the relation of unit distance of process, the attenuation quotient with change in depth that softer ray is corresponding is bigger, characterize this ray in die body in unit length decay energy more, otherwise the attenuation quotient with change in depth that harder ray is corresponding is less, characterize this ray in die body in unit length decay energy less.
In step 103, it is judged that the depth value at calculated maximal dose and the dose value at depth capacity are all the most identical with the measured value of the dose distribution on beam centre axle.
In the present embodiment, the dose distribution on 3 d water tank systematic survey beam centre axle can be utilized, obtain the measured value of dose distribution, and the dose distribution on calculated beam centre axle is compared therewith, thus judge whether current spectral distribution meets the spectral distribution of reality.
Specifically, first pass through the radiacmeter being arranged in water tank, distinguish the dose value at different depth on measurement beam central shaft, thus obtain the measured value of the dose distribution on beam centre axle；
Then, whether the relatively depth value at the depth value at calculated maximal dose and maximal dose in measured value coincide, if the error of two depth values is in the range of first threshold, then it is assumed that both are identical, if the error of two depth values exceedes first threshold scope, then it is assumed that both misfit；Equally, whether the relatively dose value at the dose value at calculated depth capacity and depth capacity in measured value coincide, if the error of two dose values is in the range of Second Threshold, then it is assumed that both are identical, if the error of two dose values exceedes Second Threshold scope, then both misfit.First threshold and Second Threshold can be set according to required precision, such as first threshold is set as 2mm, in depth value at the most calculated maximal dose and measured value, the absolute error of depth value at maximal dose is less than 2mm, then think that both coincide, Second Threshold is set as 2%, in dose value at the most calculated depth capacity and measured value, the relative error between dose value at depth capacity is less than 2%, then it is assumed that both coincide.
Compare for convenience, before the dose value at depth capacity is compared, the dose distribution normalization respectively that calculated dose distribution and measurement are obtained, preferably before the depth value at maximal dose and the dose value at depth capacity are compared, the dose distribution normalization respectively that calculated dose distribution and measurement are obtained.
If the depth value at calculated maximal dose and the dose value at depth capacity all coincide with measured value, the most current spectral distribution is the spectral distribution of beam, and calculating terminates.
If the depth value at maximal dose and at least one in the dose value at depth capacity and measured value are misfitted, then enter step 104.
In step 104, regulate current spectral distribution, and recalculate the dose distribution on described beam centre axle based on the spectral distribution after regulation.
If the depth value at calculated maximal dose is more than the depth value at maximal dose in measured value, then increases low-energy component and/or reduce high-energy component, otherwise, then reduce low-energy component and/or increase high-energy component.In the present embodiment, described low-energy component refers to that energy is less than at least one energy of peak energy, and described high-energy component refers to that energy is higher than at least one energy of peak energy.During power spectrum regulates, low-energy component and high-energy component are followed power spectrum and are changed.
With reference to shown in Fig. 3, depth value at calculated maximal dose is less than the depth value at maximal dose in measured value, need the regulation schematic diagram during energy weight reducing Low Energy Region, the weight keeping peak energy and the energy higher than peak energy is constant, by the portion of energy lower than peak energy is moved to the left, the weight of this portion of energy is reduced；
Power spectrum entirety can also be moved to the left so that the peak value of power spectrum moves towards high energy region, not only the weight of Low Energy Region be reduced, and the weight of high energy region be raised, as shown in Figure 4, all power spectrums are all moved to the left, it can be seen that the peak value of power spectrum moves towards high energy region；
It is, of course, also possible to the weight of holding peak energy and the energy lower than peak energy is constant, by the portion of energy higher than peak energy is moved to the left, the weight of this portion of energy is raised, as shown in Figure 5, it can be seen that the energy weight of high energy region adds；
When only by any of the above-described kind of mode regulate power spectrum weight two depth values can not be made to coincide time, combination in any by the way can regulate power spectrum, such as, when two depth values can not be made after only reducing Low Energy Region weight to coincide, again power spectrum entirety is moved to the left, while reducing Low Energy Region weight, raise the weight of high energy region.
If the dose value at calculated depth capacity is more than the dose value at depth capacity in measured value, then increases low-energy component and/or reduce high-energy component, otherwise, then reduce low-energy component and/or increase high-energy component.
Regulating step refer to as above regulate the description of the depth value at maximal dose, regulates power spectrum by any of the above mode or its combination in any.
During above-mentioned power spectrum regulates, the power spectrum of beam meets the Unimodal Distribution arranged in step 101 all the time；When regulating every time, need the multiple adjacent energy ingredient of adjusted in concert, i.e., the energy ingredient that every time will be above peak energy moves an identical step-length simultaneously, or the energy ingredient that will be less than peak energy moves an identical step-length, then or all energy ingredients moves an identical step-length simultaneously simultaneously.
Every time after regulation spectral distribution, it is based on the spectral distribution after regulation and recalculates the dose distribution on described beam centre axle, again return to step 103, judge that the depth value at calculated maximal dose and the dose value at the depth capacity measured value the most all with the dose distribution on beam centre axle coincide, if misfitting, then continue regulation power spectrum, until the depth value at maximal dose and the dose value at depth capacity are all identical with the measured value of the dose distribution on beam centre axle.
When regulating spectral distribution, energy ingredient moves a step-length every time, then dose distribution is calculated, and compare with measured value, if misfitting, then energy ingredient is continued to move to a step-length, again calculated dose distribution comparing, it is iterated circulation in this way, until meeting the condition of convergence.Therefore, the precision of power spectrum can be controlled by the size controlling step-length, such as, reduce step-length, the precision of power spectrum can be improved.
In other embodiments, energy ingredient can not also be moved, the weighted value of energy is directly reduced or increased, such as, when needing the weight reducing mental retardation part, the weight keeping peak energy and energetic portions is constant, the weighted value of mental retardation part is reduced, and recalculate dose distribution based on the spectral distribution after change, it may be judged whether it coincide with measured value, be iterated circulation by that analogy；
It is similar to, during regulation spectral distribution, can only regulate mental retardation or the energy weight of energetic portions, it is also possible to regulation mental retardation and the energy weight of energetic portions simultaneously, or utilize the combination in any of these three mode to be adjusted.
Should be noted, during above-mentioned regulation power spectrum, the depth value at depth value at calculated maximal dose and maximal dose in measured value is preferentially made to coincide, the dose value at dose value at calculated depth capacity and depth capacity in measured value is made to coincide and then, in addition, in addition it is also necessary to the dose value of other position on beam centre axle is all higher than or the measured value of respectively less than correspondence position.Such goal condition is easily achieved in the present embodiment.
After in depth value at calculated maximal dose and measured value, depth value at maximal dose coincide, it is adjusted merely by power spectrum when the dose value at dose value at calculated depth capacity and depth capacity in measured value can not be made to coincide, coordinate regulation basis attenuation quotient and the attenuation quotient with change in depth so that in the dose value at calculated depth capacity and measured value, the dose value at depth capacity coincide.
In other embodiments, it is also possible to utilize the dose distribution on film or other measurement device beam centre axle, the measured value of certain dose distribution can also be measured in advance by measurement apparatus and is stored in memorizer.
The method obtaining beam power spectrum of the present embodiment is based on initial spectral distribution, it is calculated dose distribution by CC algorithm, by the dose distribution of calculating is compared with measuring the dose distribution obtained, spectral distribution is regulated according to result of the comparison, therefore the measurement of power spectrum is converted into the regulation of parameter by the present invention, conveniently carrying out automation mechanized operation, computational efficiency is higher；Further, by reducing the step-length of parameter, the precision of power spectrum can be improved.
In the method for above-mentioned acquisition beam power spectrum, can only compare the depth value at maximal dose in the depth value at calculated maximal dose and measured value, both are coincide, and the dose value of other position on beam centre axle is all higher than or the measured value of respectively less than correspondence position, specifically refer to Fig. 6.
Fig. 6 is the method flow diagram obtaining beam power spectrum of another embodiment of the present invention, comprises the steps:
In step 601, the initial spectral distribution of beam is set；
Being preferably provided with described initial spectral distribution is Unimodal Distribution.
In step 602, based on initial spectral distribution, calculate the dose distribution on beam centre axle；
The basic attenuation quotient of beam is set and with the attenuation quotient of change in depth, based on initial spectral distribution, calculates the dose distribution on described beam centre axle by CC algorithm.
In step 603, it is judged that whether the depth value at calculated maximal dose coincide with the depth value measured at the maximal dose obtained；
If coincideing, the most current spectral distribution is the spectral distribution of beam, and calculating terminates, if misfitting, then enters step 604；
In step 604, regulate current spectral distribution, and recalculate the dose distribution on described beam centre axle based on the spectral distribution after regulation, until the depth value at calculated maximal dose coincide with the depth value measured at the maximal dose obtained.
The concrete ins and outs of this embodiment refer to the embodiment of Fig. 1, does not repeats them here.
Obtained the spectral distribution of beam by the method for above-mentioned acquisition beam power spectrum, utilize this spectral distribution can be calculated the dose distribution of beam.
The basic attenuation quotient of regulation beam, with change in depth attenuation quotient with off axis than, the spectral distribution obtained based on the method in Fig. 1 or Fig. 6, the described beam dose distribution at three-dimensional is calculated by CC algorithm, until the dose distribution of calculated three-dimensional is coincide with measuring the 3-dimensional dose distribution obtained, the dose distribution of the most calculated three-dimensional is final dose distribution.In the present embodiment, it is possible to use threshold value judges whether the dose distribution of calculated three-dimensional and the dose distribution measuring the three-dimensional obtained coincide, and the setting of threshold value is referred to the method in Fig. 1.
From the method in Fig. 1, the calculated dose distribution of power spectrum utilizing the method for Fig. 1 to obtain, on beam centre axle, depth value at maximal dose and the dose value at depth capacity all coincide with measured value, now it is believed that the dose distribution on calculated beam centre axle is all coincide with measured value；From the method for Fig. 6, the calculated dose distribution of power spectrum utilizing the method for Fig. 6 to obtain, on beam centre axle, in depth value at maximal dose and measured value, the depth value at maximal dose coincide, and the dose value of other position is all higher than or respectively less than measured value, basic attenuation quotient by regulation beam and the attenuation quotient with change in depth, it is easy to the dose distribution on beam centre axle is all coincide with measured value.
Such as, on beam centre axle, in depth value at maximal dose and measured value, the depth value at maximal dose coincide, and the dose value of other position is all higher than measured value, the basic attenuation quotient of regulation beam and with the attenuation quotient of change in depth, described beam dose distribution on center shaft is updated so that it is identical with the measured value of the dose distribution on beam centre axle that the dose value of other positions is progressively smaller until the dose distribution after renewal by CC algorithm.
Ratio represents on the direction vertical with beam centre axle off axis, with the increase with central shaft distance, the attenuation of dosage.The off-axis ratio of regulation beam, by CC algorithm update described beam with the dose distribution in beam centre axle vertical direction, until the measured value of the dose distribution in the dose distribution after Geng Xining and beam centre axle vertical direction is identical.
Owing to field size is different, the output factor of medical electronic linear accelerator is the most different, and regulation suitably exports the factor so that beam coincide with measured value in the dose distribution of three-dimensional on absolute value.
Such as, the regulation output factor, make on beam centre axle dose value at calculated depth capacity and the dose value measured at the depth capacity obtained completely the same on absolute value, the dose value of other depth can also be used as reference, i.e. so that on beam centre axle, the dose value of other certain position calculated and the dosage measurement value of correspondence position are completely the same on absolute value.
The dose distributions computation method of the present embodiment, first the depth value at maximal dose and the dose value at depth capacity is made all to match with measured value, and keeping the dose value of other position on beam centre axle to be all higher than or the dosage measurement value of respectively less than correspondence position, such goal condition easily realizes；Regulating other parameter again makes the dose distribution on other direction all match with measured value, and the method for this calculating power spectrum is simple, and utilizes this power spectrum calculated dose distribution accuracy higher.
Power spectrum computational methods and dosage distribution calculation method to the present invention are verified the most as one example.Fig. 7 is the initial spectral distribution of 6MV accelerator beam；Fig. 8 is the spectral distribution after the regulation of 6MV accelerator beam；Fig. 9 is the PDD curve comparison diagram with experiment curv of 6MV accelerator beam.
As a example by the accelerator of 6MV, it is respectively 100kev, 200kev, 300kev, 400kev, 500kev with energy ingredient, 600kev, 800kev, 1000kev, 1250kev, 1500kev, the ray collection of 2000kev, 3000kev, 4000kev, 5000kev, 6000kev simulates power spectrum.
nullInitial spectral distribution is set as shown in Figure 7，In figure, the straight line along the vertical direction of transverse axis distribution represents each energy ingredient，The longitudinal axis represents the weight that each energy ingredient is corresponding，Then according to the method for above-mentioned introduction is iterated circulation，Dose value at depth value at calculated maximal dose and depth capacity respectively with the maximal dose in measurement data at depth value and depth capacity at dose value coincide，Power spectrum now is as shown in Figure 8，In figure, the straight line along the vertical direction of transverse axis distribution represents each energy ingredient，The longitudinal axis represents the weight that each energy ingredient is corresponding，As can be seen from Figure 8，The energy weight of Low Energy Region increases，Process is regulated according to aforesaid power spectrum，The power spectrum variation of Fig. 7 to Fig. 8 shows corresponding PDD (percentage depth dose，It is called for short PDD，Chinese is for percentage depth dose) it is the degree of depth at corresponding maximal dose on curve and diminishes，Dosage at depth capacity have also been smaller.New power spectrum calculated PDD curve is used to refer to Fig. 9 with the contrast of experiment curv, wherein, dotted line represents calculated curve, solid line represents experiment curv, in figure, transverse axis represents the distance to die body surface of the point on beam centre axle, the i.e. degree of depth, unit cm, the longitudinal axis represents dosage and the ratio (%) of maximal dose on beam centre axle of corresponding depth, and on central shaft, the ratio at maximal dose is 100%.
As can be seen from Figure 9, utilize calculated power spectrum that dose distribution is calculated, the PDD curve obtained essentially coincides with experiment curv, concordance is fine, and, do not occur that the Rapid Dose Calculation value curve described in background technology intersects with dosage measurement curve certain position in declining district, in the case of uncomfortable energy-conservation spectrum, regulate the situation that other parameter always can not make the two ends data coincidence of intersection point all improve, in the method for the present embodiment, Rapid Dose Calculation value curve is positioned at same one side of dosage measurement curve declining district, need not regulate power spectrum, calculated dose distribution can be easily achieved match with measuring the dose distribution that obtains by regulating other parameter, measurement data simultaneously also by water tank demonstrates the accuracy of this method calculating power spectrum and dose distribution.
Accordingly, present invention also offers a kind of device 1000 obtaining beam power spectrum, please examine Figure 10, this device 1000 includes input block 1001, and for arranging the initial spectral distribution of beam, described initial spectral distribution is Unimodal Distribution；Computing unit 1002, for spectral distribution based on beam, calculates the dose distribution on beam centre axle by CC algorithm；Judging unit 1003, for judging that the depth value at calculated maximal dose and the dose value at the depth capacity measured value the most all with the dose distribution on beam centre axle coincide；Regulation unit 1004, for regulating current spectral distribution；If judging unit 1003 exports the spectral distribution that result is identical, the most current and is the spectral distribution of beam, can be exported by output unit 1005；If judging unit 1003 exports result for misfitting, then regulation unit 1004 regulates current spectral distribution, and the spectral distribution after regulation is input to computing unit 1002, computing unit 1002 based on the spectral distribution after regulation again calculate the dose distribution on described beam centre axle, the depth value at calculated maximal dose and the dose value at depth capacity all identical with measured value till.
Described regulation unit 1004 also regulates the basic attenuation quotient in CC algorithm and the attenuation quotient with change in depth so that in the dose value at calculated depth capacity and measured value, the dose value at depth capacity coincide.
Wherein, the initial spectral distribution of beam directly can be inputted by input block 1001 by user, it is also possible to directly reads from memorizer, or is directly generated according to function by computer.The detail of this device refer to the embodiment of the method for foregoing description.
The above is the preferred embodiment of the present invention; it should be pointed out that, for those skilled in the art, under the premise without departing from the principles of the invention; can also make some improvements and modifications, these improvements and modifications are also considered as protection scope of the present invention.
Claims (16)
1. the method obtaining beam power spectrum, including:
The initial spectral distribution of beam is set；
Based on initial spectral distribution, calculate the dose distribution on beam centre axle；
Judge the depth value at calculated maximal dose and the dose value at depth capacity the most all with
The measured value of the dose distribution on beam centre axle coincide；
If coincideing, the most current spectral distribution is the spectral distribution of beam；
If misfitting, then regulate current spectral distribution, and calculate institute based on the spectral distribution after regulation
State the dose distribution on beam centre axle, until the depth value at calculated maximal dose and maximum
Till the dose value of depth all coincide with measured value.
The method of acquisition beam power spectrum the most according to claim 1, it is characterised in that at the beginning of described
Beginning spectral distribution is Unimodal Distribution.
The method of acquisition beam power spectrum the most according to claim 2, it is characterised in that arrange and penetrate
The basic attenuation quotient of bundle and with the attenuation quotient of change in depth, by Collapsed Cone Rapid Dose Calculation
Method calculates the dose distribution on described beam centre axle.
The method of acquisition beam power spectrum the most according to claim 1, it is characterised in that described in penetrate
The measured value of the dose distribution on beam center axle is obtained by 3 d water tank system or film measurement.
The method of acquisition beam power spectrum the most according to claim 3, it is characterised in that judging
Dose value at calculated depth capacity whether with the measured value of the dose distribution on beam centre axle
Before Wen Heing, the beam centre that the dose distribution on calculated beam centre axle and measurement are obtained
The normalization respectively of dose distribution on axle.
The method of acquisition beam power spectrum the most according to claim 3, it is characterised in that regulation energy
Time spectrum, power spectrum meets described Unimodal Distribution all the time.
The method of acquisition beam power spectrum the most according to claim 6, it is characterised in that regulation energy
Time spectrum, the multiple adjacent energy ingredient of adjusted in concert.
The method of acquisition beam power spectrum the most according to claim 7, it is characterised in that described tune
Save current spectral distribution, including,
If the depth value at calculated maximal dose is more than the degree of depth at maximal dose in measured value
Value, then increase low-energy component and/or reduce high-energy component, otherwise, then reduce low-energy component and/or increase
High-energy component；
If the dose value at calculated depth capacity is more than the dosage at depth capacity in measured value
Value, then increase low-energy component and/or reduce high-energy component, otherwise, then reduce low-energy component and/or increase
High-energy component, wherein, described high-energy component represents at least one energy more than peak energy, described
Low-energy component represents at least one energy less than peak energy.
The method of acquisition beam power spectrum the most according to claim 8, it is characterised in that also include
Regulate described basis attenuation quotient and the attenuation quotient with change in depth so that calculated maximum is deep
In dose value at degree and measured value, the dose value at depth capacity coincide.
10. obtain a device for beam power spectrum, including,
Input block, for arranging the initial spectral distribution of beam, described initial spectral distribution is unimodal
Distribution；
Computing unit, for spectral distribution based on beam, by Collapsed Cone Rapid Dose Calculation side
Method calculates the dose distribution on beam centre axle；
Judging unit, for judging at depth value at calculated maximal dose and depth capacity
Dose value is all the most identical with the measured value of the dose distribution on beam centre axle；
Regulation unit, for regulating current spectral distribution；
If the spectral distribution that judging unit output result is identical, the most current is the spectral distribution of beam；
If judging unit output result is for misfitting, then the spectral distribution that the regulation of regulation unit is current, and
Spectral distribution after regulation is input to computing unit, and computing unit is based on the spectral distribution meter after regulation
Calculate the dose distribution on described beam centre axle, until the depth value at calculated maximal dose and
Till dose value at depth capacity all coincide with measured value.
The device of 11. acquisition beam power spectrums according to claim 10, it is characterised in that described
Regulation unit also regulates the basic attenuation quotient in Collapsed Cone dose calculation methodology and with the degree of depth
The attenuation quotient of change so that the dose value at calculated depth capacity is maximum deep with measured value
Dose value at degree coincide.
12. 1 kinds of methods obtaining beam power spectrum, including:
Arranging the initial spectral distribution of beam, described initial spectral distribution is Unimodal Distribution；
Based on described initial spectral distribution, calculated in beam by Collapsed Cone dose calculation methodology
Dose distribution in mandrel；
Judge the depth value at calculated maximal dose whether with measure at the maximal dose obtained
Depth value coincide；
If coincideing, the most current spectral distribution is the spectral distribution of beam；
If misfitting, then regulate current spectral distribution, and calculate institute based on the spectral distribution after regulation
State the dose distribution on beam centre axle, until the depth value at calculated maximal dose and measurement
Till the depth value at maximal dose obtained coincide.
13. 1 kinds of dose distributions computation methods, including,
The method described in any one of claim 1-9 is utilized to obtain the spectral distribution of beam；
The off-axis ratio of regulation beam, updates described beam by Collapsed Cone dose calculation methodology
Dose distribution, until the dose distribution after Geng Xining is coincide with measuring the dose distribution obtained, thus updates
After dose distribution be the dose distribution of calculated beam.
14. a kind of dose distributions computation methods according to claim 13, it is characterised in that also
Including the regulation output factor so that the dose distribution of beam and measured value coincide on absolute value.
15. a kind of dose distributions computation methods according to claim 13, it is characterised in that institute
State and measure the dose distribution that obtains and obtained by 3 d water tank system or film measurement.
16. 1 kinds of dose distributions computation methods, including:
The method described in claim 12 is utilized to obtain the spectral distribution of beam；
The basic attenuation quotient of regulation beam and with the attenuation quotient of change in depth, passes through Collapsed
Cone dose calculation methodology updates described beam dose distribution on center shaft, until the agent after Geng Xining
Amount distribution coincide with the measured value of the dose distribution on beam centre axle；
The off-axis ratio of regulation beam, updates described beam by Collapsed Cone dose calculation methodology and exists
With the dose distribution in beam centre axle vertical direction, until the dose distribution after Geng Xining and beam centre
The measured value of the dose distribution in axle vertical direction coincide；
Dose distribution after the renewal that three-dimensional all coincide with measured value is calculated beam and exists
The dose distribution of three-dimensional.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610176331.XA CN105866821B (en) | 2016-03-25 | 2016-03-25 | A kind of method, apparatus and dose distributions computation method obtaining beam power spectrum |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610176331.XA CN105866821B (en) | 2016-03-25 | 2016-03-25 | A kind of method, apparatus and dose distributions computation method obtaining beam power spectrum |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105866821A true CN105866821A (en) | 2016-08-17 |
| CN105866821B CN105866821B (en) | 2018-11-30 |
Family
ID=56624781
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610176331.XA Active CN105866821B (en) | 2016-03-25 | 2016-03-25 | A kind of method, apparatus and dose distributions computation method obtaining beam power spectrum |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105866821B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111060955A (en) * | 2019-12-10 | 2020-04-24 | 中国人民解放军96901部队23分队 | Neutron energy spectrum measuring device based on multilayer neutron moderated fluorescence and optical fiber conduction |
| CN112231899A (en) * | 2020-09-30 | 2021-01-15 | 天津医科大学 | Method and device for determining X-ray continuous energy spectrum, electronic equipment and storage medium |
| CN113856069A (en) * | 2021-09-30 | 2021-12-31 | 南方医科大学 | Method for fusing dose distribution by high-low energy beams |
| WO2022032591A1 (en) * | 2020-08-13 | 2022-02-17 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for modeling radiation source |
| WO2024020874A1 (en) * | 2022-07-27 | 2024-02-01 | 上海联影医疗科技股份有限公司 | Energy spectrum acquisition method and system, and storage medium |
| CN118131302A (en) * | 2024-05-07 | 2024-06-04 | 迈胜医疗设备有限公司 | Computing method and related equipment for primary energy spectrum of medical particle accelerator |
| CN118486473A (en) * | 2024-07-16 | 2024-08-13 | 华硼中子科技(杭州)有限公司 | Dose distribution displacement device design method for BNCT, storage medium, terminal and computer program product |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100998497A (en) * | 2006-12-29 | 2007-07-18 | 四川大学 | Method for determining distribution of electron beam dosage |
| CN101120871A (en) * | 2006-12-29 | 2008-02-13 | 成都川大奇林科技有限责任公司 | Precise radiotherapy planning system |
| CN101571596A (en) * | 2008-04-29 | 2009-11-04 | 清华大学 | System and method for measuring pulse type ray energy spectrum |
| US20110224963A1 (en) * | 2010-03-12 | 2011-09-15 | Synopsys, Inc. | Fast Photolithography Process Simulation to Predict Remaining Resist Thickness |
| US20120001096A1 (en) * | 2006-03-28 | 2012-01-05 | Hampton University | Hadron treatment planning with adequate biological weighting |
| CN103083820A (en) * | 2013-01-13 | 2013-05-08 | 中国科学院合肥物质科学研究院 | Dosage leading and accurately emitting treatment system |
| CN104407374A (en) * | 2014-12-10 | 2015-03-11 | 上海联影医疗科技有限公司 | Dose distribution determination method and device |
| CN105204056A (en) * | 2015-09-29 | 2015-12-30 | 广州瑞多思医疗科技有限公司 | On-line dose monitor and verification method in radiotherapy |
-
2016
- 2016-03-25 CN CN201610176331.XA patent/CN105866821B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120001096A1 (en) * | 2006-03-28 | 2012-01-05 | Hampton University | Hadron treatment planning with adequate biological weighting |
| CN100998497A (en) * | 2006-12-29 | 2007-07-18 | 四川大学 | Method for determining distribution of electron beam dosage |
| CN101120871A (en) * | 2006-12-29 | 2008-02-13 | 成都川大奇林科技有限责任公司 | Precise radiotherapy planning system |
| CN101571596A (en) * | 2008-04-29 | 2009-11-04 | 清华大学 | System and method for measuring pulse type ray energy spectrum |
| US20110224963A1 (en) * | 2010-03-12 | 2011-09-15 | Synopsys, Inc. | Fast Photolithography Process Simulation to Predict Remaining Resist Thickness |
| CN103083820A (en) * | 2013-01-13 | 2013-05-08 | 中国科学院合肥物质科学研究院 | Dosage leading and accurately emitting treatment system |
| CN104407374A (en) * | 2014-12-10 | 2015-03-11 | 上海联影医疗科技有限公司 | Dose distribution determination method and device |
| CN105204056A (en) * | 2015-09-29 | 2015-12-30 | 广州瑞多思医疗科技有限公司 | On-line dose monitor and verification method in radiotherapy |
Non-Patent Citations (5)
| Title |
|---|
| 周正东等: "Collapsed Cone 光子束剂量计算方法研究", 《中国生物医学工程学报》 * |
| 康念坎等: "多能谱电子束能量分布的蒙特卡罗模拟和空间电荷效应", 《电子科学学刊》 * |
| 康盛伟等: "电子束能谱宽度及角分布对蒙特卡罗方法模拟计算剂量分布影响", 《中华放射肿瘤学杂志》 * |
| 张从华等: "关于百分深度剂量计算组织最大剂量比的处理方法及数据对比", 《中国医学物理学杂志》 * |
| 李诗豪等: "百分深度剂量的计算方法", 《中国医疗器械杂志》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111060955A (en) * | 2019-12-10 | 2020-04-24 | 中国人民解放军96901部队23分队 | Neutron energy spectrum measuring device based on multilayer neutron moderated fluorescence and optical fiber conduction |
| WO2022032591A1 (en) * | 2020-08-13 | 2022-02-17 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for modeling radiation source |
| CN112231899A (en) * | 2020-09-30 | 2021-01-15 | 天津医科大学 | Method and device for determining X-ray continuous energy spectrum, electronic equipment and storage medium |
| CN112231899B (en) * | 2020-09-30 | 2023-02-03 | 天津医科大学 | Method and device for determining X-ray continuous energy spectrum, electronic equipment and storage medium |
| CN113856069A (en) * | 2021-09-30 | 2021-12-31 | 南方医科大学 | Method for fusing dose distribution by high-low energy beams |
| CN113856069B (en) * | 2021-09-30 | 2023-11-03 | 南方医科大学 | Method for fusing high-low energy beam with dose distribution |
| WO2024020874A1 (en) * | 2022-07-27 | 2024-02-01 | 上海联影医疗科技股份有限公司 | Energy spectrum acquisition method and system, and storage medium |
| CN118131302A (en) * | 2024-05-07 | 2024-06-04 | 迈胜医疗设备有限公司 | Computing method and related equipment for primary energy spectrum of medical particle accelerator |
| CN118131302B (en) * | 2024-05-07 | 2024-08-30 | 迈胜医疗设备有限公司 | Computing method and related equipment for primary energy spectrum of medical particle accelerator |
| CN118486473A (en) * | 2024-07-16 | 2024-08-13 | 华硼中子科技(杭州)有限公司 | Dose distribution displacement device design method for BNCT, storage medium, terminal and computer program product |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105866821B (en) | 2018-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105866821A (en) | Method and device for obtaining beam energy spectrum, and dose distribution calculating method | |
| Ojala et al. | Performance of dose calculation algorithms from three generations in lung SBRT: comparison with full Monte Carlo‐based dose distributions | |
| Fragoso et al. | Dosimetric verification and clinical evaluation of a new commercially available Monte Carlo-based dose algorithm for application in stereotactic body radiation therapy (SBRT) treatment planning | |
| Fracchiolla et al. | Characterization and validation of a Monte Carlo code for independent dose calculation in proton therapy treatments with pencil beam scanning | |
| Bush et al. | Monte Carlo simulation of RapidArc radiotherapy delivery | |
| Yan et al. | Clinical implementation and evaluation of the Acuros dose calculation algorithm | |
| CN104117151A (en) | Optimization method of online self-adaption radiotherapy plan | |
| CN105617536A (en) | Inverse intensity modulated arc radioterapy optimization method and device | |
| Mackin et al. | Improving spot‐scanning proton therapy patient specific quality assurance with HPlusQA, a second‐check dose calculation engine | |
| Gordon et al. | Evaluation of dosimetric margins in prostate IMRT treatment plans | |
| Li et al. | Impact of dose calculation accuracy during optimization on lung IMRT plan quality | |
| JP6588158B2 (en) | Dose error distribution calculation device and particle beam therapy device equipped with dose error distribution calculation device | |
| Caron et al. | Deterministic model for the transport of energetic particles: Application in the electron radiotherapy | |
| Verhaegen et al. | Incorporating dynamic collimator motion in Monte Carlo simulations: an application in modelling a dynamic wedge | |
| Li et al. | Extension and validation of a GPU‐Monte Carlo dose engine gDPM for 1.5 T MR‐LINAC online independent dose verification | |
| US9878181B2 (en) | Method and device for computing physical and/or radiobiological quantities induced by an ion beam in a volume | |
| Cho et al. | Practical implementation of a collapsed cone convolution algorithm for a radiation treatment planning system | |
| JP2021097997A (en) | Obtaining energy spectrum of focused ion beam | |
| US12053645B2 (en) | System and method for robust treatment planning in radiation therapy | |
| Boylan et al. | The use of a realistic VMAT delivery emulator to optimize dynamic machine parameters for improved treatment efficiency | |
| Kalet et al. | Validation of the photon convolution‐superposition algorithm applied to fast neutron beams | |
| Yuan et al. | A virtual source model for Monte Carlo simulation of helical tomotherapy | |
| Lorenz et al. | An independent dose calculation algorithm for MLC‐based stereotactic radiotherapy | |
| Chen et al. | TOPAS Monte Carlo simulation for a scanning proton therapy system in SPHIC | |
| Craig et al. | Commissioning a fast Monte Carlo dose calculation algorithm for lung cancer treatment planning |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |