CN102930109B - Dose model matching method and system - Google Patents
Dose model matching method and system Download PDFInfo
- Publication number
- CN102930109B CN102930109B CN201210450646.0A CN201210450646A CN102930109B CN 102930109 B CN102930109 B CN 102930109B CN 201210450646 A CN201210450646 A CN 201210450646A CN 102930109 B CN102930109 B CN 102930109B
- Authority
- CN
- China
- Prior art keywords
- parameter
- curve
- simulant
- tps
- dose
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 230000008569 process Effects 0.000 claims abstract description 31
- 238000005259 measurement Methods 0.000 claims description 37
- 238000004088 simulation Methods 0.000 claims description 17
- 230000008878 coupling Effects 0.000 claims description 10
- 238000010168 coupling process Methods 0.000 claims description 10
- 238000005859 coupling reaction Methods 0.000 claims description 10
- 230000001105 regulatory effect Effects 0.000 abstract description 6
- 238000012067 mathematical method Methods 0.000 description 13
- 238000001959 radiotherapy Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012913 prioritisation Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 238000002922 simulated annealing Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 238000011179 visual inspection Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a dose model matching method and a system. The method comprises following steps: a, a measure dose curve is obtained and a simulative dose curve is calculated according to a parameter of transactions per second (TPS); b, a matching degree between the measure dose curve and the simulative dose curve is calculated and the matching degree needs to be judged whether the matching degree is larger than a first threshold value; c, if the matching degree is larger than the first threshold value, the parameter of the TPS is regulated and the simulative dose curve is recalculated according to the regulated parameter of the TPS; d, the step b and the step c are carried out repeatedly until the matching degree is not larger than the first threshold value, and then the process is over.
Description
Technical field
The present invention relates to technical field of data processing, particularly a kind of dose modal matching process and system.
Background technology
Radiotherapy TPS (Treatment Planning System) software is a kind of medical software being applied to Patients During Radiotherapy at present, its role is to the dosage distribution produced by the irradiation human body of ray in simulation Patients During Radiotherapy, formulate more reasonably treatment plan.So for TPS, most important is exactly that its analog result will as much as possible with the process of actual therapeutic, and the dosage distribution that accelerator produces radiation exposure human body is consistent.The consistance of the two is higher then illustrates that the performance of TPS software is better.
In order to ensure the accuracy of TPS software, all can another process of debugging in the use procedure of TPS.Namely first can irradiate die body practically with accelerator, and the dosage distribution record of generation is got off, import in TPS software.The dosage distribution utilizing dose modal to calculate with TPS again distributes with the dosage imported and compares, and the difference of the two is less, illustrates that the degree of accuracy of TPS is higher.If degree of accuracy does not meet the demands, dose modal can be re-established by regulating accordingly the parameters in TPS dose modal simulation process, and calculated dose distribution, until its degree of accuracy reaches requirement.Generally speaking TPS generates according to the dosage distribution that reality produces and measures dose curve, then generates simulant discharge curve according to the dosage distribution of simulation.By directly contrasting the coincidence degree of two curves, namely have expressed the accuracy of TPS intuitively.
In fact, due to the inevitable otherness of prototype part, the dosage distribution that each accelerator for generation of radiotherapy ray produces is incomplete same.TPS, before any one accelerator of simulation, has to pass through debugging, and namely TPS needs, for a set of corresponding parameter of each accelerator setting, namely to set up a dose modal, mate the actual conditions of accelerator with this dose modal.If without debugging and the adjustment of parameter, TPS software can not the dosage of all accelerators of accurate simulation distribute.
But the parameter affecting the distribution of TPS simulation dosage reaches dozens or even hundreds of.And the general also dependence of the TPS debugging of present stage is artificial, do not realize automated debug.When each artificial adjustment, need rule of thumb first to change a parameter, then utilize TPS to carry out Rapid Dose Calculation, more whether observe the registration measuring dose curve and simulant discharge curve higher.If registration is higher, retain this value, if otherwise readjusted.But because parameter is numerous, so repeatedly revise needs and expend a lot of time, and the height of registration relies on the judgement of naked eyes completely, does not have the support in gross data.Artificial adjustment process for an accelerator generally reaches about two days, and efficiency is very low, not only bring hard work amount to staff, and accuracy is also stable not.
Summary of the invention
In view of this, the object of the present invention is to provide a kind of dose modal matching process and system, make TPS dose modal realize the dosage distribution of Auto-matching accelerator.
For achieving the above object, the present invention has following technical scheme:
A kind of dose modal matching process, said method comprising the steps of:
The measurement dose curve of a, acquisition accelerator; Simulant discharge curve is calculated according to the parameter of TPS;
The matching degree of b, computation and measurement dose curve and simulant discharge curve, and judge whether this matching degree is greater than first threshold;
If c matching degree is greater than first threshold, adjust TPS parameter, recalculate simulant discharge curve according to the TPS parameter after adjustment;
D, repeated execution of steps b-c, until matching degree is not more than first threshold, then set up dose modal according to the parameter after adjustment.
The matching degree of described computation and measurement dose curve and simulant discharge curve is specially:
Pre-establish correspondence function, the coordinate figure measuring dose curve and simulant discharge curve is substituted in correspondence function, using the value of correspondence function that the obtains matching degree as measurement dose curve and simulant discharge curve.
The described correspondence function that pre-establishes is specially:
Get the value of quadratic sum as correspondence function of measurement dose curve and simulant discharge curve corresponding point coordinate difference, then the expression formula of described correspondence function is;
K is variable, and represent abscissa value, on C (k) representative simulation dose curve, horizontal ordinate is the coordinate figure of the point of k, and it is the coordinate figure of the point of k that horizontal ordinate on dose curve is measured in M (k) representative; The value of f and correspondence function.
The described correspondence function that pre-establishes is specially:
Get the quadratic sum measuring part corresponding point coordinate difference in dose curve and simulant discharge curve open country; Be added the coupling measurement dose curve of weight and the quadratic sum of simulant discharge curve penumbra region part corresponding point slope differences again, as the value of correspondence function; Then the expression formula of described correspondence function is,
K is variable, represents abscissa value; C
ak in () representative simulation dose curve u'eno, part horizontal ordinate is the coordinate figure of the point of k; M
ak part horizontal ordinate in dose curve u'eno is measured in () representative is the coordinate figure of the point of k, and ω represents weight coefficient, T
con (k+ △ k) representative simulation dose curve, horizontal ordinate is the slope of the point of k+ △ k; T
mk horizontal ordinate on dose curve is measured in () representative is the slope of the point of k; △ k represents the side-play amount of two curve penumbra region parts.
Described adjustment TPS parameter specifically comprises the following steps:
Find and measure the dose curve point coordinate maximum relative to simulant discharge curve side-play amount;
According to the coordinate figure of this point coordinate, the TPS parameter affecting this point coordinate is selected to adjust.
The described TPS parameter affecting this point coordinate of selecting is specially:
For likely affecting the parametric distribution weight probability of this point coordinate, and according to a kind of parameter that may affect this point coordinate of weight probability Stochastic choice.
Described adjustment TPS parameter is specially:
Pre-set Second Threshold, and make TPS parameter increase or deduct Second Threshold.
A kind of dose modal matching system, described system specifically comprises:
Analog module, for calculating simulant discharge curve according to the parameter of TPS;
Matching module, for obtaining the measurement dose curve of accelerator and described simulant discharge curve, the matching degree of computation and measurement dose curve and simulant discharge curve;
Judge module, for judging whether the matching degree measuring dose curve and simulant discharge curve is greater than first threshold, then terminates to mate flow process when matching degree is not more than first threshold;
Adjusting module, for when described matching degree is greater than first threshold, adjustment TPS parameter, recalculates simulant discharge curve according to the TPS parameter after adjustment, and is sent to matching module re-starts coupling by recalculating the simulant discharge curve obtained;
MBM, for setting up dose modal according to the parameter after adjustment.
Described matching module specifically comprises:
Function unit, for pre-establishing and preserving correspondence function;
First computing unit, for obtaining the measurement dose curve of accelerator, correspondence function and simulant discharge curve; The coordinate figure measuring dose curve and simulant discharge curve is substituted in correspondence function, using the value of correspondence function that the obtains matching degree as measurement dose curve and simulant discharge curve.
Described adjusting module is specially:
Coordinate selection unit, measures the dose curve point coordinate maximum relative to simulant discharge curve side-play amount for finding;
Parameter selection unit, for the coordinate figure according to this point coordinate, selects the TPS parameter affecting this point coordinate;
Parameter adjustment unit, for adjusting by the TPS parameter selected;
Second computing unit, for recalculating simulant discharge curve according to the TPS parameter after adjustment;
Transmitting element, for being sent to matching module re-starting coupling by recalculating the simulant discharge curve obtained.
Described adjusting module also comprises:
Weight-assigning unit, for the parametric distribution weight probability for likely affecting this point coordinate;
Then described parameter selection unit is according to the coordinate figure of this point coordinate, selects the TPS parameter affecting this point coordinate to be specially, according to a kind of parameter that may affect this point coordinate of weight probability Stochastic choice.
Described adjusting module also comprises:
Second Threshold unit, for arranging Second Threshold;
Then described parameter adjustment unit adjustment is specially by the TPS parameter selected, and makes selected TPS parameter increase or deduct Second Threshold.
As seen through the above technical solutions, the automatic Debugging that present invention achieves TPS dose modal with mate, greatly improve the efficiency of TPS debugging, shorten the time of debugging, and improve the matching degree measuring dose curve and simulant discharge curve, the simulation degree of accuracy making TPS distribute for dosage is higher.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is dose profile form schematic diagram in prior art;
Fig. 2 is method flow diagram described in the embodiment of the present invention;
Fig. 3 is method flow diagram described in another embodiment of the present invention;
Fig. 4 is system architecture schematic diagram described in the embodiment of the present invention.
Embodiment
For making the object of the embodiment of the present invention, technical scheme and advantage clearly, below in conjunction with the accompanying drawing in the embodiment of the present invention, technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.
The Auto-matching that the dosage distribution that the object of the invention is to the dose modal simulation realizing TPS software distributes with the actual dosage produced of accelerator, improve the degree of accuracy of coupling, and it is consuming time to reduce coupling.First the description conveniently in follow-up specific embodiment and explanation, first following problem is made explanations herein:
Shown in Figure 1 is the concrete form of generally dose profile.In following examples of the present invention, dose profile is divided into measures dose curve and simulant discharge curve two kinds.The acquisition mode measuring dose curve is, transmitter generation radiation exposure die body produces dosage distribution, and utilizes the dosage distribution of this actual generation to form measurement dose curve; Simulant discharge curve is then that the dosage distribution that TPS software utilizes dose modal to simulate is formed.
Visible by Fig. 1, the geometric shape of routine dose profile is "Ji" type, and the horizontal ordinate of this distribution curve represents on die body that each point is to the distance of die body central point, and namely the initial point of horizontal ordinate represents the central point of die body; The ordinate of this curve represents dose value.Its top slope is less, almost parallel with transverse axis part be known as curve open country in part; Two sideleads are comparatively large, are almost called as the penumbra region part of curve with the part of longitudinal axis; The part that lowermost end two parts slope is less is called as field part.
Shown in Figure 2, be a specific embodiment of dose modal matching process of the present invention, in the present embodiment, said method comprising the steps of:
The measurement dose curve of step 201, acquisition accelerator; Simulant discharge curve is calculated according to the parameter of TPS;
Due to the inevitable otherness of prototype part, the dosage distribution that each accelerator for generation of radiotherapy ray produces is incomplete same.The dose modal of TPS software needs the adjustment according to being carried out parameter by the actual conditions of accelerator simulated.So in order to this accelerator of TPS dose modal adjustment parameter matching, then need to use this accelerator generation radiotherapy radiation exposure die body in advance, thus the distribution of the generation dosage of reality, and form measurement dose curve according to the dosage distribution that reality produces.Namely TPS measures the benchmark of dose curve as debugging using this in the process of debugging.
And then in TPS parameters, analogue accelerator dosage distribution, the dosage distribution situation simulated according to it calculates simulant discharge curve.The method of the invention in fact namely makes the method for simulant discharge curve Matched measurement dose curve by adjustment TPS parameter.
The matching degree of step 202, computation and measurement dose curve and simulant discharge curve, and judge whether this matching degree is greater than first threshold;
After obtaining measurement dose curve and simulant discharge curve, first need the matching degree both being calculated by mathematical method.In existing TPS debug process, general only by the coincidence degree of visual inspection two curves; But in the present embodiment, be incorporated herein the mathematical criterion that mathematical method tries to achieve the two matching degree.So, just achieve the matching degree of determining quantitative analysis two curves, be no longer subject to the impact of subjective judgement.
The unspecific concrete mathematical method limiting calculating matching degree in the present embodiment; But present stage, usual mathematical method generally can not exceed the scope of least square method or similar mathematical method.If for modal least square method, the matching degree index that it calculates is 0, then two curve perfect matchings are described, matching degree index value larger both explanations matching degree is lower; Generally also there is this feature in other similar mathematical methods.A first threshold is set according to this feature, if matching degree is not more than first threshold namely illustrate that both can meet the demands at matching degree in the present embodiment.The concrete numerical value of first threshold is arranged according to concrete user demand and mathematical method.In fact matching degree can think that TPS software has been debugged successfully after meeting the demands, and stops the adjustment of TPS parameter, process ends.
If step 203 matching degree is greater than first threshold, adjust TPS parameter, recalculate simulant discharge curve according to the TPS parameter after adjustment;
If the matching degree obtained in previous step is greater than first threshold, then illustrate that the matching degree of two curves is undesirable, regenerate simulant discharge curve after needing adjustment TPS parameter to generate, and then recalculate the matching degree of two curves.Specifically, in TPS software debugging, the parameter affecting simulant discharge curve is numerous, and concrete need adjust which or which parameter and matching degree index can be made to level off to 0, selects not unique; And so-called parameter adjustment, specifically a certain parameter increased or reduce, and increasing or the link such as the amplitude that reduces, also needing to decide as the case may be; The selection mode that the adjustment of parameter faces in practical operation is varied, cannot enumerate.Concrete restriction is not made in the present embodiment.
In the TPS debugging of present stage, those skilled in that art completely can according to mathematical theory and artificial the adjusting parameter of operating experience, thus improve simulant discharge curve and the matching degree measuring dose curve as much as possible, make matching degree index level off to 0.In the present embodiment when setting up dose modal, can, in conjunction with theory in the past and experience, program be utilized to realize adjustment TPS parameter automatically equally.
Step 204, repeated execution of steps 202-step 203, until matching degree is not more than first threshold, then set up dose modal according to the parameter after adjustment.
The simulant discharge curve regenerated in step 203 is returned step 202, recalculates simulant discharge curve and the matching degree measuring dose curve.And the judgement of matching degree is carried out according to identical flow process.If matching degree is still greater than first threshold, then again perform step 203, further adjustment TPS parameter.The execution step 202-step 203 of circulation like this, until the simulant discharge curve generated after repeatedly parameter adjustment is no longer greater than first threshold with the matching degree measuring dose curve, then think that parameter adjustment is successful, then set up dose modal according to the parameter after adjustment, namely newly-built dose modal can mate the actual conditions of this accelerator, can simulate the dosage distribution of this accelerator accurately.
In addition, also it should be noted that, if in certain parameter adjustment once, again the matching degree index obtained after adjustment does not more level off to 0, numerical value is larger on the contrary, then illustrate that mistake appears in last adjustment, the matching degree adjusting rear two curves reduces on the contrary.So in this case, last adjustment is no longer retained; Namely by before all Parameter reconstruction to last adjustment, then re-start adjustment.
Due to the basic embodiment that the present embodiment is described method, so do not make unnecessary restriction to concrete mathematical method and the process step relevant to mathematical method.Visible by above technical scheme, the beneficial effect that the present embodiment exists is, utilize the matching degree of mathematical method computation and measurement dose curve and simulant discharge curve, and founding mathematical models realizes program automatically adjusts TPS parameter, avoid the complicated processes of artificial adjustment TPS parameter, also eliminate the subjectivity of artificial judgment matching degree; Improve efficiency and the accuracy of TPS parameter adjustment.
Shown in Figure 3 is another specific embodiment of dose modal matching process of the present invention.The present embodiment, mainly on basis embodiment illustrated in fig. 2, has carried out concrete supplementary notes for the mathematical algorithm related to, and by description that a part of sport technique segment is specialized.The present embodiment is a preferred embodiment of the method for the invention.
Concrete, in the present embodiment, said method comprising the steps of:
The measurement dose curve of step 301, acquisition accelerator; Simulant discharge curve is calculated according to the parameter of TPS;
Step 301 described in the present embodiment with embodiment illustrated in fig. 2 in step 201 completely the same.
Step 302, pre-establishing correspondence function, substituting in correspondence function by measuring the coordinate figure of dose curve with simulant discharge curve, using the value of correspondence function that the obtains matching degree as measurement dose curve and simulant discharge curve;
After obtaining measurement dose curve and simulant discharge curve, first need the matching degree both being calculated by mathematical method.In the present embodiment, the mathematical method adopted sets up a correspondence function in advance, and substituted in the middle of function by coordinate figure corresponding for two curves.
Concrete, formulate correspondence function according to least square method in the present embodiment, the expression formula of correspondence function is as follows:
Wherein k is variable, and represent abscissa value, on C (k) representative simulation dose curve, horizontal ordinate is the coordinate figure of the point of k, and it is the coordinate figure of the point of k that horizontal ordinate on dose curve is measured in M (k) representative; The f then obtained represents the quadratic sum of corresponding point coordinate difference on two curves; In this, as the index weighing Curve Matching degree in the present embodiment; Namely the value of f measures the matching degree of dose curve and simulant discharge curve.
If two curve perfect matching, so f equals 0.But it is completely the same to be in fact difficult to accomplish to make two curves.The value of f more levels off to 0, then illustrate that the matching degree of two curves is higher.
Step 303, judge whether the matching degree measuring dose curve and simulant discharge curve is greater than first threshold; If matching degree is greater than first threshold, then enters step 304 matching degree be not more than first threshold if adjust TPS parameter, adjust TPS parameter then termination process;
According to the feature of correspondence function, set a first threshold.If the f value obtained is not more than first threshold, then illustrate that the matching degree of two curves meets the real needs in the present embodiment, stops the adjustment of TPS parameter, process ends; Otherwise then illustrate that matching degree is undesirable, must carry out the adjustment of TPS parameter.
Step 304, then adjust TPS parameter when matching degree is greater than first threshold, then find and measure the dose curve point coordinate maximum relative to simulant discharge curve side-play amount;
The result that correspondence function draws only has indicated the height of both matching degrees, if but matching degree is lower, then and correspondence function cannot indicate both unmatched reasons.In the middle of curve as shown in Figure 1, the optional position of whole piece curve all likely becomes matching degree the poorest a bit (namely measure the dose curve point maximum relative to simulant discharge curve side-play amount, side-play amount is larger, illustrates and does not mate).According to experience in the past, the form of curve diverse location is subject to different parameters impact; So the situation that side-play amount is maximum appears in different positions, the design parameter needing specific aim to carry out adjusting is not identical yet.In the process of adjustment parameter, must first determine to mate the poorest point coordinate, namely find and measure the dose curve point coordinate maximum relative to simulant discharge curve side-play amount.
The mode of the point coordinate that concrete discovery side-play amount is maximum, automatically can be realized by program completely in present stage this area, not repeat at this.
Step 305, according to the coordinate figure measuring the dose curve point maximum relative to simulant discharge curve side-play amount, the TPS parameter affecting this point coordinate is selected to adjust;
The coordinate figure of the point that side-play amount is maximum is different, so also not identical on its parameter producing main impact.In fact according to experience and the theoretic analysis of TPS debugging in the past, the scope of coordinate figure and the corresponding relation of major effect parameter have been well known in the art.Find this point coordinate value can determine the parameter needing to carry out regulating accordingly.
After determining this parameter, can adjust this parameter.Set up the mathematical model of adjustment parameter in the present embodiment, and a Second Threshold is set, apply this mathematical model when each adjustment parameter, the initial value of this parameter is increased or deducts Second Threshold.In general Second Threshold is larger, and the effect of adjustment is more obvious, but degree of accuracy is relatively low, and arranging of Second Threshold can be determined as the case may be.In addition, in the process of regulating parameter, concrete selection increases Second Threshold or deducts Second Threshold, then positive and negative depending on side-play amount; The positive and negative corresponding relation increased and decreased with parameter of side-play amount is similarly known in the art.
Step 306, according to adjustment after TPS parameter recalculate simulant discharge curve;
After adjustment parameter, TPS recalculates simulant discharge curve according to parameter after adjustment, follow-up again with new simulant discharge curve continuation Matched measurement dose curve.
Step 307, repeated execution of steps 302-step 306, until matching degree is not more than first threshold, then set up dose modal according to the parameter after adjustment.
The simulant discharge curve regenerated in step 306 is returned step 302, recalculate simulant discharge curve and the matching degree measuring dose curve, and judge whether to need continue adjustment TPS parameter and how to adjust TPS parameter according to the matching degree newly calculated, the description of said process and step 302-step 306 is completely the same.The execution step 302-step 306 of circulation like this, until the simulant discharge curve generated after repeatedly parameter adjustment is no longer greater than first threshold with the matching degree measuring dose curve, then think that parameter adjustment is successful, then set up dose modal according to the parameter after adjustment, namely newly-built dose modal can mate the actual conditions of this accelerator, can simulate the dosage distribution of this accelerator accurately.
In addition, also it should be noted that, if in certain parameter adjustment once, again the matching degree index obtained after adjustment does not more level off to 0, numerical value is larger on the contrary, then illustrate that mistake appears in last adjustment, the matching degree adjusting rear two curves reduces on the contrary.So in this case, last adjustment is no longer retained; Namely by before all Parameter reconstruction to last adjustment, then re-start adjustment.
The present embodiment is on basis embodiment illustrated in fig. 2, and the beneficial effect realized further is: in the present embodiment, technical scheme is more complete, open more abundant; Achieve the automatic adjustment of TPS parameter, avoid the complicated processes of artificial adjustment TPS parameter, also eliminate the subjectivity of artificial judgment matching degree; Improve efficiency and the accuracy of TPS parameter adjustment.
On basis embodiment illustrated in fig. 3, its correspondence function can also be changed and obtain another following specific embodiment.All the other steps of following examples and embodiment illustrated in fig. 3 completely the same, difference is only the difference of correspondence function.
The parameter related to due to TPS is of a great variety, and a kind of correspondence function is not enough to for all parameters.Time the correspondence function obtained by least square method in embodiment illustrated in fig. 3 is for spectrum parameters, electronic pollution parameter and wild intrinsic parameter, user demand can be met.But for other a part of parameters such as such as Gaussian width parameter, above-mentioned correspondence function then cannot be suitable for.
For the adjustment of Gaussian width parameter, measurement dose curve and simulant discharge curve interior part out of office can be made to overlap as far as possible, and penumbra region, both sides part cannot overlap, can only slope be equal as much as possible, namely parallel.For reflecting this characteristic, the correspondence function expression formula for Gaussian width parameter is as follows:
Wherein k is variable, represents abscissa value; C
ak in () representative simulation dose curve u'eno, part horizontal ordinate is the coordinate figure of the point of k; M
ak part horizontal ordinate in dose curve u'eno is measured in () representative is the coordinate figure of the point of k, so first of this formula sum term represents wild interior matching degree partly; ω represents weight coefficient, utilizes this weight coefficient to weigh second sum term for the impact of matching degree, T in this formula
con (k+ △ k) representative simulation dose curve, horizontal ordinate is the slope of the point of k+ △ k; T
mk horizontal ordinate on dose curve is measured in () representative is the slope of the point of k; △ k represents the side-play amount of two curve penumbra region parts; Second sum term of this formula represents the matching degree of penumbra region part.
It should be noted that, due to the parallel offset of two curve penumbra regions, so the impact of side-play amount and △ k must be considered when contrasting its slope, so by T
c(k+ △ k) and T
mk () is carried out wrong item and is subtracted each other.In fact △ k can be 0, and positive number or negative, be determined on a case-by-case basis.△ k can be determined by mathematical algorithm in the present embodiment, and the mathematical algorithm calculating △ k is not unique, is not construed as limiting, can be obtained the value of △ k those skilled in the art by mathematical algorithm in the present embodiment.
Can also comprise in the present embodiment simultaneously embodiment illustrated in fig. 2 described in above-mentioned correspondence function in correspondence function and the present embodiment.Select accordingly when adjustment different parameters.The beneficial effect that the present embodiment exists is that the correspondence function by optimizing makes overall technical architecture more complete, and the matching degree after dose profile is finally adjusted is higher.
Illustrate according to aforesaid explanation, known in curve, its top slope is less, almost parallel with transverse axis part be known as curve open country in part; Two sideleads are comparatively large, are almost called as the penumbra region part of curve with the part of longitudinal axis; The part that lowermost end two parts slope is less is called as field part.But above-mentioned saying is only for the judgement of vision.Be applied in the middle of mathematical computations, wild interior part and penumbra region part can be distinguished by slope, specifically can arrange a slope threshold value, some slope on curve being less than slope threshold value is classified as wild interior part, and the point that slope is greater than slope threshold value is classified as penumbra region part.
In addition, for above embodiment, also there is a kind of prioritization scheme selecting mechanism based on simulated annealing, the coordinate figure according to measuring the dose curve point maximum relative to simulant discharge curve side-play amount can being applied in, select to affect in the process of TPS parameter of this point coordinate.Simulated annealing selects mechanism to be present stage a kind of known method, and therefore not to repeat here, is applied in the method for the invention, specific as follows:
In previous embodiment, for the selection mode of parameter to be adjusted, disclose following content: according to experience in the past, the form of curve diverse location is subject to different parameters impact; So the situation that side-play amount is maximum appears in different positions, the design parameter needing specific aim to carry out adjusting is not identical yet.In fact according to experience and the theoretic analysis of TPS debugging in the past, the scope of coordinate figure and the corresponding relation of major effect parameter have been well known in the art.Find this point coordinate value can determine the parameter needing to carry out regulating accordingly.
In fact, also often just a kind of parameter of the tracing pattern within the scope of a certain coordinate figure of major effect.For example, if in the scope that coordinate figure is in 30-40 time, there are A, B, C tri-kinds of parameters and significantly impact produced for the form of curve; But Primary regulation process can not adjust three kinds of parameters again simultaneously, so therefrom can only get one adjustment.So in the process selected, first according to the size of three kinds of parameters to influence power within the scope of curve 30-40, be often kind of parametric distribution weight probability, such as, for the parameter A probability that assigns weight is 60%, the weight probability distributed for parameter B is 30%, and the weight probability distributed for parameter C is 10%; Namely illustrating that the parameter playing main impact is A, is secondly that the impact that parameter B, parameter C play is minimum.Next, based on the probability distributed, in A, B, C tri-parameters, Stochastic choice one regulates.Namely select the probability of A to be 60%, the probability selecting B is 30%, and the probability selecting C is 10%.The many kinds of parameters of the adjustment possibility influence curve balanced can be realized by the way.
With reference to shown in Fig. 4, it is a specific embodiment of dose modal matching system of the present invention.System described in the present embodiment is corresponding above-mentioned dose modal matching process and exists, and the technical scheme of the two is consistent in itself, the present embodiment and the embodiment of the method shown in corresponding diagram 1.
In the present embodiment, described system comprises:
Analog module, for calculating simulant discharge curve according to the parameter of TPS;
Matching module, for obtaining the measurement dose curve of accelerator and described simulant discharge curve, the matching degree of computation and measurement dose curve and simulant discharge curve;
Judge module, for judging whether the matching degree measuring dose curve and simulant discharge curve is greater than first threshold, then terminates to mate flow process when matching degree is not more than first threshold;
Adjusting module, for when described matching degree is greater than first threshold, adjustment TPS parameter, recalculates simulant discharge curve according to the TPS parameter after adjustment, and is sent to matching module re-starts coupling by recalculating the simulant discharge curve obtained;
MBM, for setting up dose modal according to the parameter after adjustment.
The beneficial effect that the present embodiment exists is, utilize the matching degree of mathematical method computation and measurement dose curve and simulant discharge curve, and founding mathematical models realizes program automatically adjusts TPS parameter, avoid the complicated processes of artificial adjustment TPS parameter, also eliminate the subjectivity of artificial judgment matching degree; Improve efficiency and the accuracy of TPS parameter adjustment.
Accordingly, for the various prioritization schemes in described embodiment of the method, the system embodiment shown in Fig. 4 can also be optimized, specifically can have following:
Described matching module specifically comprises:
Function unit, for pre-establishing and preserving correspondence function;
First computing unit, for obtaining the measurement dose curve of accelerator, correspondence function and simulant discharge curve; The coordinate figure measuring dose curve and simulant discharge curve is substituted in correspondence function, using the value of correspondence function that the obtains matching degree as measurement dose curve and simulant discharge curve.
Described adjusting module is specially:
Coordinate selection unit, measures the dose curve point coordinate maximum relative to simulant discharge curve side-play amount for finding;
Parameter selection unit, for the coordinate figure according to this point coordinate, selects the TPS parameter affecting this point coordinate;
Parameter adjustment unit, for adjusting by the TPS parameter selected;
Second computing unit, for recalculating simulant discharge curve according to the TPS parameter after adjustment;
Transmitting element, for being sent to matching module re-starting coupling by recalculating the simulant discharge curve obtained.
Described adjusting module also comprises:
Weight-assigning unit, for the parametric distribution weight probability for likely affecting this point coordinate;
Then described parameter selection unit is according to the coordinate figure of this point coordinate, selects the TPS parameter affecting this point coordinate to be specially, according to a kind of parameter that may affect this point coordinate of weight probability Stochastic choice.
Described adjusting module also comprises:
Second Threshold unit, for arranging Second Threshold;
Then described parameter adjustment unit adjustment is specially by the TPS parameter selected, and makes selected TPS parameter increase or deduct Second Threshold.
The above is only 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 also should be considered as protection scope of the present invention.
Claims (12)
1. a dose modal matching process, is characterized in that, said method comprising the steps of:
The measurement dose curve of a, acquisition accelerator; Simulant discharge curve is calculated according to the parameter of TPS;
The matching degree of b, computation and measurement dose curve and simulant discharge curve, and judge whether this matching degree is greater than first threshold;
If c matching degree is greater than first threshold, adjust TPS parameter, recalculate simulant discharge curve according to the TPS parameter after adjustment;
D, repeated execution of steps b-c, until matching degree is not more than first threshold, then set up dose modal according to the parameter after adjustment.
2. method according to claim 1, it is characterized in that, the matching degree of described computation and measurement dose curve and simulant discharge curve is specially:
Pre-establish correspondence function, the coordinate figure measuring dose curve and simulant discharge curve is substituted in correspondence function, using the value of correspondence function that the obtains matching degree as measurement dose curve and simulant discharge curve.
3. method according to claim 2, is characterized in that, described in pre-establish correspondence function and be specially:
Get the value of quadratic sum as correspondence function of measurement dose curve and simulant discharge curve corresponding point coordinate difference, then the expression formula of described correspondence function is;
K is variable, and represent abscissa value, on C (k) representative simulation dose curve, horizontal ordinate is the coordinate figure of the point of k, and it is the coordinate figure of the point of k that horizontal ordinate on dose curve is measured in M (k) representative; The value of f and correspondence function.
4. method according to claim 2, is characterized in that, described in pre-establish correspondence function and be specially:
Get the quadratic sum measuring part corresponding point coordinate difference in dose curve and simulant discharge curve open country; Be added the coupling measurement dose curve of weight and the quadratic sum of simulant discharge curve penumbra region part corresponding point slope differences again, as the value of correspondence function; Then the expression formula of described correspondence function is,
K is variable, represents abscissa value; C
ak in () representative simulation dose curve u'eno, part horizontal ordinate is the coordinate figure of the point of k; M
ak part horizontal ordinate in dose curve u'eno is measured in () representative is the coordinate figure of the point of k, and ω represents weight coefficient, T
con (k+ △ k) representative simulation dose curve, horizontal ordinate is the slope of the point of k+ △ k; T
mk horizontal ordinate on dose curve is measured in () representative is the slope of the point of k; △ k represents the side-play amount of two curve penumbra region parts.
5. method according to claim 1, it is characterized in that, described adjustment TPS parameter specifically comprises the following steps:
Find and measure the dose curve point coordinate maximum relative to simulant discharge curve side-play amount;
According to the coordinate figure of this point coordinate, the TPS parameter affecting this point coordinate is selected to adjust.
6. method according to claim 5, is characterized in that, the described TPS parameter affecting this point coordinate of selecting is specially:
For likely affecting the parametric distribution weight probability of this point coordinate, and according to a kind of parameter that may affect this point coordinate of weight probability Stochastic choice.
7. method according to claim 1-6 any one, is characterized in that, described adjustment TPS parameter is specially:
Pre-set Second Threshold, and make TPS parameter increase or deduct Second Threshold.
8. a dose modal matching system, is characterized in that, described system specifically comprises:
Analog module, for calculating simulant discharge curve according to the parameter of TPS;
Matching module, for obtaining the measurement dose curve of accelerator and described simulant discharge curve, the matching degree of computation and measurement dose curve and simulant discharge curve;
Judge module, for judging whether the matching degree measuring dose curve and simulant discharge curve is greater than first threshold, then terminates to mate flow process when matching degree is not more than first threshold;
Adjusting module, for when described matching degree is greater than first threshold, adjustment TPS parameter, recalculates simulant discharge curve according to the TPS parameter after adjustment, and is sent to matching module re-starts coupling by recalculating the simulant discharge curve obtained;
MBM, for setting up dose modal according to the parameter after adjustment.
9. system according to claim 8, it is characterized in that, described matching module specifically comprises:
Function unit, for pre-establishing and preserving correspondence function;
First computing unit, for obtaining the measurement dose curve of accelerator, correspondence function and simulant discharge curve; The coordinate figure measuring dose curve and simulant discharge curve is substituted in correspondence function, using the value of correspondence function that the obtains matching degree as measurement dose curve and simulant discharge curve.
10. system according to claim 8, it is characterized in that, described adjusting module is specially:
Coordinate selection unit, measures the dose curve point coordinate maximum relative to simulant discharge curve side-play amount for finding;
Parameter selection unit, for the coordinate figure according to this point coordinate, selects the TPS parameter affecting this point coordinate;
Parameter adjustment unit, for adjusting by the TPS parameter selected;
Second computing unit, for recalculating simulant discharge curve according to the TPS parameter after adjustment;
Transmitting element, for being sent to matching module re-starting coupling by recalculating the simulant discharge curve obtained.
11. systems according to claim 10, it is characterized in that, described adjusting module also comprises:
Weight-assigning unit, for the parametric distribution weight probability for likely affecting this point coordinate;
Then described parameter selection unit is according to the coordinate figure of this point coordinate, selects the TPS parameter affecting this point coordinate to be specially, according to a kind of parameter that may affect this point coordinate of weight probability Stochastic choice.
12. according to claim 10 or 11 system, it is characterized in that, described adjusting module also comprises:
Second Threshold unit, for arranging Second Threshold;
Then described parameter adjustment unit adjustment is specially by the TPS parameter selected, and makes selected TPS parameter increase or deduct Second Threshold.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210450646.0A CN102930109B (en) | 2012-11-12 | 2012-11-12 | Dose model matching method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210450646.0A CN102930109B (en) | 2012-11-12 | 2012-11-12 | Dose model matching method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102930109A CN102930109A (en) | 2013-02-13 |
| CN102930109B true CN102930109B (en) | 2015-07-01 |
Family
ID=47644906
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210450646.0A Expired - Fee Related CN102930109B (en) | 2012-11-12 | 2012-11-12 | Dose model matching method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102930109B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103565459B (en) * | 2013-09-25 | 2015-07-01 | 沈阳东软医疗系统有限公司 | Method and device for judging imaging quality |
| CN107016228B (en) * | 2016-11-25 | 2023-08-18 | 斯图尔特平李 | Medical imaging system based on HMDS |
| CN108573751A (en) * | 2017-03-08 | 2018-09-25 | 深圳大森智能科技有限公司 | The method of inspection and device of therapeutic effect |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1438602A (en) * | 2003-03-17 | 2003-08-27 | 吴大可 | Radiotherapy system of automatic computing prescription dosage |
-
2012
- 2012-11-12 CN CN201210450646.0A patent/CN102930109B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN102930109A (en) | 2013-02-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6783398B2 (en) | Systems and methods for modeling dose calculations in radiation therapy planning | |
| CN104338240B (en) | Automatic optimization device for on-line self-adaption radiotherapy plan | |
| US10076673B2 (en) | Interactive dose gradient based optimization technique to control IMRT delivery complexity | |
| US7542545B2 (en) | Method and device for calculating the radiation dose distribution for a radiation treatment system for the purpose of radiation therapy of an animal body | |
| CN112581475B (en) | Method for predicting gamma passing rate of radiotherapy plan | |
| CN110415785A (en) | The method and system of artificial intelligence guidance radiotherapy planning | |
| CN102930109B (en) | Dose model matching method and system | |
| CN105447330B (en) | The weight regulating method and device of intensity-modulated radiation therapy | |
| CN109166613A (en) | Radiotherapy treatment planning assessment system and method based on machine learning | |
| CN106029170A (en) | Incremental treatment planning | |
| US20220008748A1 (en) | Pareto Optimal Projection Search (POPS) for Automated Radiation Therapy Treatment Planning | |
| CN111784203B (en) | Electric power spot market risk simulation analysis method suitable for generator set participation | |
| CN112133440A (en) | Dose calculation system, dose quality assurance system and storage medium | |
| CN114420252A (en) | Method, device and medium for determining intensity modulated radiotherapy plan evaluation parameter matrix | |
| CN115938564A (en) | Tumor radiotherapy plan design method and device, electronic device and storage medium | |
| WO2022088341A1 (en) | Radiation therapy control method and apparatus | |
| JP2021524782A (en) | Systems, Computer Program Products, and Methods for Ion-Based Radiation Therapy Planning | |
| US20240157171A1 (en) | Tumor radiotherapy planning design method and apparatus, electronic device and storage medium | |
| CN109447574A (en) | Assets based on Fuzzy Optimum Neural Network turn solid project processing method | |
| EP4257181A1 (en) | Methods and system related to radiotherapy treatment planning | |
| Szekeres | Checking the Evidence for Declining Discount Rates | |
| CN113521559B (en) | Flux iteration optimization method for calculating dosage by adopting MC, TPS system and manufacturing method | |
| US20240189624A1 (en) | Dose coverage margin metric for evaluating radiotherapy treatment plans | |
| Panda et al. | Analysis of critical conditions in electric power systems by feed forward and layer recurrent neural networks | |
| CN114461695A (en) | Distribution network operation data analysis method |
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 | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180529 Address after: 110167 room 310, 177-2 innovation road, Hunnan District, Shenyang, Liaoning. Patentee after: SHENYANG NEUSOFT ZHIRUI RADIOTHERAPY TECHNOLOGY Co.,Ltd. Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16 Patentee before: SHENYANG NEUSOFT MEDICAL SYSTEMS Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150701 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |