CN102930109A - Dose model matching method and system - Google Patents

Dose model matching method and system Download PDF

Info

Publication number
CN102930109A
CN102930109A CN2012104506460A CN201210450646A CN102930109A CN 102930109 A CN102930109 A CN 102930109A CN 2012104506460 A CN2012104506460 A CN 2012104506460A CN 201210450646 A CN201210450646 A CN 201210450646A CN 102930109 A CN102930109 A CN 102930109A
Authority
CN
China
Prior art keywords
matching degree
parameter
curve
simulant
tps
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012104506460A
Other languages
Chinese (zh)
Other versions
CN102930109B (en
Inventor
姚鹏
王玉
马锐兵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenyang Neusoft Zhi Rui Radiotherapy Technology Co., Ltd.
Original Assignee
Neusoft Medical Systems Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Neusoft Medical Systems Co Ltd filed Critical Neusoft Medical Systems Co Ltd
Priority to CN201210450646.0A priority Critical patent/CN102930109B/en
Publication of CN102930109A publication Critical patent/CN102930109A/en
Application granted granted Critical
Publication of CN102930109B publication Critical patent/CN102930109B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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

A kind of dose modal matching process and system
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 that is applied at present Patients During Radiotherapy, its role is to the dosage distribution by the irradiation human body generation of ray in the simulation Patients During Radiotherapy, formulate more reasonably treatment plan.So for TPS, most important be exactly its analog result will be as much as possible with the process of actual therapeutic in, the dosage that accelerator produces the radiation exposure human body distributes consistent.The consistance of the two is higher to illustrate that then the performance of TPS software is better.
In order to guarantee the accuracy of TPS software, process that all can another debugging in the use procedure of TPS.Namely can shine practically die body with accelerator first, and the dosage distribution record that produces is got off, import in the TPS software.The dosage that utilizes dose modal to calculate with TPS again distributes and compares with the dosage distribution that imports, and the difference of the two is less to illustrate that then the degree of accuracy of TPS is higher.If not meeting the demands, degree of accuracy can rebulid dose modal by the parameters in the corresponding adjusting TPS dose modal simulation process, and calculated dose distribution, until its degree of accuracy reaches requirement.Generally speaking TPS generates the measurement dose curve according to the dosage distribution of reality generation, again according to the dosage distribution generation simulant discharge curve of simulating.By the coincidence degree of two curves of direct contrast, namely expressed intuitively the accuracy of TPS.
In fact, because the inevitable otherness of former device, it is incomplete same that each dosage that produces for generation of the accelerator of radiotherapy ray distributes.TPS must be through debugging before any one accelerator of simulation, and namely TPS need to set the corresponding parameter of a cover for each accelerator, namely sets up a dose modal, the actual conditions of mating accelerator with this dose modal.If without the adjustment of debugging and parameter, a TPS software can not be simulated the dosage distribution of all accelerators accurately.
But the parameter that affects the distribution of TPS simulation dosage reaches dozens or even hundreds of.And the TPS of present stage debugging generally also relies on manually, does not realize automated debug.In the time of each artificial adjustment, need rule of thumb to change first a parameter, then utilize TPS to carry out Rapid Dose Calculation, whether highlyer observe again the registration of measuring dose curve and simulant discharge curve.Keep this value if registration is higher, if otherwise readjust.But because parameter is numerous, so repeatedly revises and to expend a lot of times, and the height of registration does not have the support on the gross data fully by means of the judgement of naked eyes.Artificial adjustment process for an accelerator generally reached about two days, and efficient is very low, not only brings the hard work amount to the 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 the TPS dose modal realize that the dosage of Auto-matching accelerator distributes.
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; Calculation of parameter according to TPS goes out to simulate dose curve;
The matching degree of b, computation and measurement dose curve and simulant discharge curve, and judge that whether this matching degree is greater than first threshold;
If the c matching degree is then adjusted the TPS parameter greater than first threshold, recomputate the simulant discharge curve according to the TPS parameter after adjusting;
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 adjusting.
The matching degree of described computation and measurement dose curve and simulant discharge curve is specially:
Pre-establish the matching degree function, will measure in the coordinate figure substitution matching degree function of dose curve and simulant discharge curve, with the value of the matching degree function that the obtains matching degree as measurement dose curve and simulant discharge curve.
The described matching degree function that pre-establishes is specially:
Get the quadratic sum of measurement dose curve and simulant discharge curve corresponding point coordinate difference as the value of matching degree function, then the expression formula of described matching degree function is;
f = Σ k [ C ( k ) - M ( k ) ] 2
K is variable, the expression abscissa value, and horizontal ordinate is the coordinate figure of the point of k on C (k) the representative simulation dose curve, horizontal ordinate is the coordinate figure of the point of k on M (k) the representative measurement dose curve; F is the value of matching degree function.
The described matching degree function that pre-establishes is specially:
Get the quadratic sum of measuring part corresponding point coordinate difference in dose curve and the simulant discharge curve open country; The measurement dose curve of weight and the quadratic sum of simulant discharge curve penumbra region part corresponding point slope differences are mated in addition again, as the value of matching degree function; Then the expression formula of described matching degree function is,
f = Σ k [ C A ( k ) - M A ( k ) ] 2 + ω Σ k [ T C ( k + Δk ) - T M ( k ) ] 2
K is variable, the expression abscissa value; C A(k) the part horizontal ordinate is the coordinate figure of the point of k in the representative simulation dose curve u'eno; M A(k) the part horizontal ordinate is the coordinate figure of the point of k in the representative measurement dose curve u'eno, and ω represents weight coefficient, T CHorizontal ordinate is the slope of the point of k+ △ k on (k+ △ k) representative simulation dose curve; T M(k) horizontal ordinate is the slope of the point of k on the representative measurement dose curve; △ k represents the side-play amount of two curve penumbra region parts.
Described adjustment TPS parameter specifically may further comprise the steps:
Find and measure dose curve with respect to the point coordinate of simulant discharge curve side-play amount maximum;
According to the coordinate figure of this point coordinate, the TPS parameter of selecting to affect this point coordinate is adjusted.
The TPS parameter that described selection affects this point coordinate is specially:
For might affecting the parametric distribution weight probability of this point coordinate, and select a kind of parameter that may affect this point coordinate according to the weight probability is random.
Described adjustment TPS parameter is specially:
Set in advance Second Threshold, and make the increase of TPS parameter or deduct Second Threshold.
A kind of dose modal matching system, described system specifically comprises:
Analog module is used for going out to simulate dose curve according to the calculation of parameter of TPS;
Matching module is for the measurement dose curve and the described simulant discharge curve that obtain accelerator, the matching degree of computation and measurement dose curve and simulant discharge curve;
Judge module is used for judging that the matching degree of measuring dose curve and simulant discharge curve whether greater than first threshold, then finishes to mate flow process when matching degree is not more than first threshold;
Adjusting module is used for during greater than first threshold, adjusting the TPS parameter in described matching degree, recomputates the simulant discharge curve according to the TPS parameter after adjusting, and will recomputate the simulant discharge curve that obtains and be sent to matching module and re-start coupling;
MBM is used for setting up dose modal according to the parameter after adjusting.
Described matching module specifically comprises:
Function unit is used for pre-establishing and preserving the matching degree function;
The first computing unit is for the measurement dose curve, matching degree function and the simulant discharge curve that obtain accelerator; To measure in the coordinate figure substitution matching degree function of dose curve and simulant discharge curve, with the value of the matching degree function that the obtains matching degree as measurement dose curve and simulant discharge curve.
Described adjusting module is specially:
The coordinate selection unit is used for finding and measures dose curve with respect to the point coordinate of simulant discharge curve side-play amount maximum;
Parameter selection unit is used for the coordinate figure according to this point coordinate, selects to affect the TPS parameter of this point coordinate;
The parameter adjustment unit is used for adjusting selecteed TPS parameter;
The second computing unit is used for recomputating the simulant discharge curve according to the TPS parameter after adjusting;
Transmitting element is used for recomputating the simulant discharge curve that obtains and is sent to matching module and re-starts coupling.
Described adjusting module also comprises:
The weight allocation unit, being used to affect the parametric distribution of this point coordinate weight probability;
Then described parameter selection unit is according to the coordinate figure of this point coordinate, and the TPS parameter of selecting to affect this point coordinate is specially, according to a kind of parameter that may affect this point coordinate of the random selection of weight probability.
Described adjusting module also comprises:
The Second Threshold unit is used for arranging Second Threshold;
Then described parameter adjustment unit is adjusted selecteed TPS parameter and is specially, and selected TPS parameter is increased or deducts Second Threshold.
As seen through the above technical solutions, the present invention has realized automatic Debugging and the coupling of TPS dose modal, improved greatly the efficient of TPS debugging, shortened the time of debugging, and improved the matching degree of measuring dose curve and simulant discharge curve, so that TPS is higher for the simulation degree of accuracy that dosage distributes.
Description of drawings
In order to be illustrated more clearly in the embodiment of the invention or technical scheme of the prior art, the below will do to introduce simply to the accompanying drawing of required use in embodiment or the description of the Prior Art, apparently, accompanying drawing in the following describes is some embodiments of the present invention, for those of ordinary skills, under the prerequisite of not paying creative work, can also obtain according to these accompanying drawings other accompanying drawing.
Fig. 1 is dosage distribution curve form schematic diagram in the prior art;
Fig. 2 is the described method flow diagram of the embodiment of the invention;
Fig. 3 is the described method flow diagram of another embodiment of the present invention;
Fig. 4 is the described system architecture schematic diagram of the embodiment of the invention.
Embodiment
For the purpose, technical scheme and the advantage that make the embodiment of the invention clearer, below in conjunction with the accompanying drawing in the embodiment of the invention, technical scheme in the embodiment of the invention is clearly and completely described, obviously, described embodiment is the present invention's part embodiment, rather than whole embodiment.Based on the embodiment among the present invention, those of ordinary skills belong to the scope of protection of the invention not making the every other embodiment that obtains under the creative work prerequisite.
The dosage that the object of the invention is to realize the dose modal simulation of the TPS software Auto-matching that dosage with actual generations of accelerator distributes that distributes improve the degree of accuracy of coupling, and minimizing is mated consuming time.At first for description and explanation in the convenient follow-up specific embodiment, first following problem is made explanations herein:
Referring to Figure 1 shows that the generally speaking concrete form of dosage distribution curve.In following examples of the present invention, the dosage distribution curve is divided into measures two kinds of dose curve and simulant discharge curves.The acquisition mode of measuring dose curve is that transmitter generation radiation exposure die body produces dosage and distributes, and utilizes the dosage distribution of this actual generation to form the measurement dose curve; The simulant discharge curve then is that TPS software utilizes the dosage of dose modal simulation to distribute to form.
By Fig. 1 as seen, the geometric shape of routine dosage distribution curve is "Ji" type, and the horizontal ordinate of this distribution curve represents that each point is to the distance of die body central point on the die body, and the initial point of horizontal ordinate namely represents the central point of die body; The ordinate of this curve represents dose value.Its top slope is less, almost with the part of the transverse axis part in curve wild of being known as; Two sideleads are larger, and almost the part parallel with the longitudinal axis is called as the penumbra region part of curve; The part that lowermost end two parts slope is less is called as open-air part.
Referring to shown in Figure 2, a specific embodiment for 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; Calculation of parameter according to TPS goes out to simulate dose curve;
Because the inevitable otherness of former device, it is incomplete same that each dosage that produces for generation of the accelerator of radiotherapy ray distributes.The dose modal of TPS software needs to carry out the adjustment of parameter according to the actual conditions of the accelerator that simulated.So for this accelerator of TPS dose modal adjustment parameter matching, then need to use in advance this accelerator to produce radiotherapy radiation exposure die body, thereby the generation dosage of reality distributes, and forms the measurement dose curve according to the dosage distribution that reality produces.TPS namely measures dose curve as the benchmark of debugging with this in the process of debugging.
And then in TPS parameters, the dosage of analogue accelerator distributes, and calculates the simulant discharge curve according to its dosage distribution situation of simulating.The method of the invention in fact namely makes the method for simulant discharge curve Matched measurement dose curve by adjusting the TPS parameter.
The matching degree of step 202, computation and measurement dose curve and simulant discharge curve, and judge that whether this matching degree is greater than first threshold;
After having obtained measurement dose curve and simulant discharge curve, at first need to calculate by mathematical method the matching degree of the two.In existing TPS debug process, general the coincidence degree by two curves of visual inspection; But in the present embodiment, be incorporated herein the mathematics index that mathematical method is tried to achieve the two matching degree.So, just realize the matching degree of two curves of quantitative analysis, no longer be subject to the impact of subjective judgement.
The concrete mathematical method of matching degree is calculated in the present embodiment unspecific restriction; But present stage, habitual mathematical method generally can not exceed the scope of least square method or similar mathematical method.If take modal least square method as example, the matching degree index that it calculates is 0, and two curve perfect matchings then are described, the two matching degree of the larger explanation of matching degree index value is lower; Generally also there are these characteristics in other similar mathematical methods.Set a first threshold according to these characteristics in the present embodiment, illustrate that namely both matching degrees can meet the demands if matching degree is not more than first threshold.The concrete numerical value of first threshold arranges according to concrete user demand and mathematical method.In fact matching degree can think that TPS software debugged successfully after meeting the demands, and stops the adjustment of TPS parameter, process ends.
If step 203 matching degree is then adjusted the TPS parameter greater than first threshold, recomputate the simulant discharge curve according to the TPS parameter after adjusting;
If the matching degree that obtains in the previous step, illustrates then that the matching degree of two curves is undesirable greater than first threshold, need to regenerate the simulant discharge curve after the generation of adjustment TPS parameter, and then recomputate the matching degree of two curves.Specifically, in the TPS software debugging, the parameter that affects the simulant discharge curve is numerous, and concrete need are adjusted which or which parameter energy so that the matching degree index levels off to 0, select not unique; And so-called parameter adjustment specifically increases a certain parameter or dwindles, and the links such as amplitude that increase or dwindle, and also needs to decide as the case may be; The selection mode that the adjustment of parameter faces in practical operation is varied, can't enumerate one by one.Do not make concrete restriction in the present embodiment.
In the TPS debugging of present stage, those skilled in that art fully can be according to mathematical theory and artificial the adjusting parameter of operating experience, thereby improve as much as possible simulant discharge curve and the matching degree of measuring dose curve, so that the matching degree index levels off to 0.When setting up dose modal, can in conjunction with theory and experience in the past, utilize program to realize adjusting automatically the TPS parameter equally in the present embodiment.
Step 204, repeated execution of steps 202-step 203 until matching degree is not more than first threshold, are then set up dose modal according to the parameter after adjusting.
The simulant discharge curve that regenerates in the step 203 is returned step 202, recomputate simulant discharge curve and the matching degree of measuring dose curve.And the identical flow process of foundation is carried out the judgement of matching degree.If matching degree is still greater than first threshold, then execution in step 203 is again further adjusted the TPS parameter.The execution in step 202-step 203 that so circulates, until through the simulant discharge curve that generates after the parameter adjustment repeatedly and the matching degree of measuring dose curve no longer greater than first threshold, think that then parameter adjustment is successful, then set up dose modal according to the parameter after adjusting, newly-built dose modal namely can mate the actual conditions of this accelerator, and the dosage that can simulate accurately this accelerator distributes.
In addition, also need to prove, if in certain parameter adjustment once, again the matching degree index that obtains after the adjustment does not more level off to 0, numerical value is larger on the contrary, and then mistake appears in the last adjustment of explanation, and the matching degree of adjusting rear two curves has reduced on the contrary.So in this case, no longer keep last adjustment; Namely all parameters are returned to before last the adjustment, re-start again adjustment.
Because the present embodiment is a basic embodiment of described method, so concrete mathematical method and the process step relevant with mathematical method are not made unnecessary restriction.By above technical scheme as seen, 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 set up mathematical model realization program and automatically adjust the TPS parameter, avoid the complicated processes of artificial adjustment TPS parameter, also got rid of the subjectivity of artificial judgment matching degree; Efficient and the accuracy of TPS parameter adjustment have been improved.
Referring to another specific embodiment that Figure 3 shows that dose modal matching process of the present invention.The present embodiment mainly is on basis embodiment illustrated in fig. 2, has carried out concrete supplementary notes for the mathematical algorithm that relates to, and a part of sport technique segment has been carried out the description of specializing.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; Calculation of parameter according to TPS goes out to simulate dose curve;
Step 301 described in the present embodiment is in full accord with the step 201 in embodiment illustrated in fig. 2.
Step 302, pre-establish the matching degree function, will measure in the coordinate figure substitution matching degree function of dose curve and simulant discharge curve, with the value of the matching degree function that the obtains matching degree as measurement dose curve and simulant discharge curve;
After having obtained measurement dose curve and simulant discharge curve, at first need to calculate by mathematical method the matching degree of the two.In the present embodiment, the mathematical method that adopts is to set up in advance a matching degree function, and in the middle of the coordinate figure substitution function that two curves is corresponding.
Concrete, formulate the matching degree function according to least square method in the present embodiment, the expression formula of matching degree function is as follows:
f = Σ k [ C ( k ) - M ( k ) ] 2 ;
Wherein k is variable, the expression abscissa value, and horizontal ordinate is the coordinate figure of the point of k on C (k) the representative simulation dose curve, horizontal ordinate is the coordinate figure of the point of k on M (k) the representative measurement dose curve; The f that then obtains represents the quadratic sum of corresponding point coordinate difference on two curves; In the present embodiment with this as the index of weighing the Curve Matching degree; The value of f is namely measured the matching degree of dose curve and simulant discharge curve.
If curve perfect matching two, f equals 0 so.But in fact be difficult to accomplish to make two curves in full accord.The value of f more levels off to 0, illustrates that then the matching degree of two curves is higher.
Whether step 303, judgement measure the matching degree of dose curve and simulant discharge curve greater than first threshold; If then enter step 304 matching degree and be not more than first threshold and then adjust then termination process of TPS parameter if matching degree is then adjusted the TPS parameter greater than first threshold;
According to the characteristics of matching degree function, set a first threshold.If the f value that obtains is not more than first threshold, the matching degree that two curves then are described satisfies the real needs in the present embodiment, stops the adjustment of TPS parameter, process ends; Otherwise illustrate that then matching degree is undesirable, must carry out the adjustment of TPS parameter.
Step 304, then adjust the TPS parameter when matching degree greater than first threshold, then find and measure dose curve with respect to the point coordinate of simulant discharge curve side-play amount maximum;
The result that the matching degree function draws has only indicated the height of both matching degrees, if but matching degree is lower, and then the matching degree function can't indicate both unmatched reasons.In the middle of the curve as shown in Figure 1, the optional position of whole piece curve all might become matching degree the poorest a bit (namely measure dose curve with respect to the point of simulant discharge curve side-play amount maximum, the larger then explanation of side-play amount is not more mated).According to experience in the past, the form of curve diverse location is subject to the different parameters impact; So the situation of side-play amount maximum appears in different positions, the design parameter that needs specific aim to adjust is not identical yet.In the process of adjusting parameter, must at first determine the point coordinate that coupling is the poorest, namely find and measure dose curve with respect to the point coordinate of simulant discharge curve side-play amount maximum.
Can automatically realize by program fully in the mode of the point coordinate of concrete discovery side-play amount maximum, present stage this area, not do at this and give unnecessary details.
Step 305, according to measuring the coordinate figure of dose curve with respect to the point of simulant discharge curve side-play amount maximum, the TPS parameter of selecting to affect this point coordinate is adjusted;
The coordinate figure of the point of side-play amount maximum is different, and is also not identical on its parameter that produces main impact so.In fact according to experience and the theoretic analysis of in the past TPS debugging, the corresponding relation of the scope of coordinate figure and major effect parameter has been well known in the art.Find this point coordinate value can determine accordingly the parameter that need to regulate.
After determining this parameter, can adjust this parameter.Set up to adjust the mathematical model of parameter in the present embodiment, and a Second Threshold is set, use this mathematical model each when adjusting parameter, the initial value of this parameter is increased or deduct Second Threshold.In general Second Threshold is larger, and the effect of adjustment is more obvious, but degree of accuracy is relatively low, and the arranging to come as the case may be of Second Threshold determined.In addition, the concrete selection increases Second Threshold or deducts Second Threshold in the process of regulating parameter, then decides on the positive and negative of side-play amount; The positive and negative corresponding relation with the parameter increase and decrease of side-play amount is similarly known in the art.
Step 306, recomputate the simulant discharge curve according to the TPS parameter after adjusting;
Adjust after the parameter, TPS recomputates the simulant discharge curve according to parameter after adjusting, and is 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, are then set up dose modal according to the parameter after adjusting.
The simulant discharge curve that regenerates in the step 306 is returned step 302, recomputate simulant discharge curve and the matching degree of measuring dose curve, and judging whether that according to the matching degree that newly calculates needs continue to adjust the TPS parameter and how to adjust the TPS parameter, the description of said process and step 302-step 306 is in full accord.The execution in step 302-step 306 that so circulates, until through the simulant discharge curve that generates after the parameter adjustment repeatedly and the matching degree of measuring dose curve no longer greater than first threshold, think that then parameter adjustment is successful, then set up dose modal according to the parameter after adjusting, newly-built dose modal namely can mate the actual conditions of this accelerator, and the dosage that can simulate accurately this accelerator distributes.
In addition, also need to prove, if in certain parameter adjustment once, again the matching degree index that obtains after the adjustment does not more level off to 0, numerical value is larger on the contrary, and then mistake appears in the last adjustment of explanation, and the matching degree of adjusting rear two curves has reduced on the contrary.So in this case, no longer keep last adjustment; Namely all parameters are returned to before last the adjustment, re-start again adjustment.
The present embodiment is on basis embodiment illustrated in fig. 2, and the beneficial effect of further realizing is: technical scheme is more complete in the present embodiment, and is open more abundant; Realize the automatic adjustment of TPS parameter, avoided the complicated processes of artificial adjustment TPS parameter, also got rid of the subjectivity of artificial judgment matching degree; Efficient and the accuracy of TPS parameter adjustment have been improved.
On basis embodiment illustrated in fig. 3, can also change its matching degree function and obtain another following specific embodiment.All the other steps of following examples and embodiment illustrated in fig. 3 in full accord, difference only is the difference of matching degree function.
Because the parameter that TPS relates to is of a great variety, a kind of matching degree function is not enough to for all parameters.The matching degree function that obtains by least square method in embodiment illustrated in fig. 3 can satisfy user demand for spectrum parameters, electronic pollution parameter and wild intrinsic parameter the time.But for such as other a part of parameters such as Gaussian width parameters, above-mentioned matching degree function then can't be suitable for.
For the adjustment of Gaussian width parameter, the measurement dose curve is overlapped as far as possible with simulant discharge curve interior part out of office, and penumbra region, both sides part can't overlap, slope equates as much as possible, and is namely parallel.For reflecting this characteristic, as follows for the matching degree function expression of Gaussian width parameter:
f = Σ k [ C A ( k ) - M A ( k ) ] 2 + ω Σ k [ T C ( k + Δk ) - T M ( k ) ] 2
Wherein k is variable, the expression abscissa value; C A(k) the part horizontal ordinate is the coordinate figure of the point of k in the representative simulation dose curve u'eno; M A(k) the part horizontal ordinate is the coordinate figure of the point of k in the representative measurement dose curve u'eno, so first sum term of this formula represents the matching degree of wild interior part; ω represents weight coefficient, utilizes this weight coefficient to weigh second sum term for the impact of matching degree, T in this formula CHorizontal ordinate is the slope of the point of k+ △ k on (k+ △ k) representative simulation dose curve; T M(k) horizontal ordinate is the slope of the point of k on the representative measurement dose curve; △ 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.
Need to prove, because the parallel offset of two curve penumbra regions, so must consider when contrasting its slope that side-play amount is the impact of △ k, so with T C(k+ △ k) and T M(k) carrying out the mistake item subtracts each other.In fact △ k can be 0, and positive number or negative are determined on a case-by-case basis.△ k can determine by mathematical algorithm in the present embodiment, and the mathematical algorithm that calculates △ k is not unique, is not construed as limiting in the present embodiment, can obtain by mathematical algorithm the value of △ k those skilled in the art.
Can also comprise simultaneously in the present embodiment embodiment illustrated in fig. 2 described in above-mentioned matching degree function in matching degree function and the present embodiment.When adjusting different parameters, select accordingly.The beneficial effect that the present embodiment exists is so that overall technical architecture is more complete, so that the matching degree that the dosage distribution curve is final after adjusting is higher by the matching degree function optimized.
Explain according to aforesaid, in curve, its top slope is less as can be known, almost with the part of the transverse axis part in curve wild of being known as; Two sideleads are larger, and almost the part parallel with the longitudinal axis is called as the penumbra region part of curve; The part that lowermost end two parts slope is less is called as open-air part.But above-mentioned saying is only for the judgement of vision.Be applied in the middle of the mathematical computations, can distinguish wild interior part and penumbra region part by slope, a slope threshold value specifically can be set, with slope on the curve less than the point of slope threshold value be classified as wild in part, slope is classified as the penumbra region part greater than the point of slope threshold value.
In addition, for above embodiment, also there is a kind of prioritization scheme of selecting mechanism based on simulated annealing, can be applied in according to measuring the coordinate figure of dose curve with respect to the point of simulant discharge curve side-play amount maximum, select to affect in the process of TPS parameter of this point coordinate.It is known a kind of method of present stage that simulated annealing is selected mechanism, and therefore not to repeat here, and it is applied in the method for the invention, specific as follows:
In the 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 the different parameters impact; So the situation of side-play amount maximum appears in different positions, the design parameter that needs specific aim to adjust is not identical yet.In fact according to experience and the theoretic analysis of in the past TPS debugging, the corresponding relation of the scope of coordinate figure and major effect parameter has been well known in the art.Find this point coordinate value can determine accordingly the parameter that need to regulate.
In fact, the tracing pattern in a certain coordinate figure scope of major effect often is not a kind of parameter yet.For instance, if when coordinate figure is in the scope of 30-40, exist A, B, three kinds of parameters of C to produce significantly impact for the form of curve; But the Primary regulation process can not be adjusted three kinds of parameters again simultaneously, so can only therefrom get one adjustment.So in the process of selecting, at first according to the size of three kinds of parameters to influence power in the curve 30-40 scope, for every kind of parameter is distributed a weight probability, for example, for the parameter A probability that assigns weight is 60%, the weight probability that distributes for B parameter is 30%, and the weight probability that distributes for parameters C is 10%; It is A that the parameter of main impact is played in i.e. explanation, is B parameter secondly, and the impact that parameters C plays is minimum.Next take the probability that distributes as the basis, select one and regulate at random in A, B, three parameters of C.The probability of namely selecting A is 60%, and selecting the probability of B is 30%, and selecting the probability of C is 10%.Can realize by the way the many kinds of parameters of the adjustment possibility influence curve of balance.
With reference to shown in Figure 4, be 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 existing, and the technical scheme of the two is consistent in itself, and the present embodiment i.e. embodiment of the method corresponding shown in Figure 1.
In the present embodiment, described system comprises:
Analog module is used for going out to simulate dose curve according to the calculation of parameter of TPS;
Matching module is for the measurement dose curve and the described simulant discharge curve that obtain accelerator, the matching degree of computation and measurement dose curve and simulant discharge curve;
Judge module is used for judging that the matching degree of measuring dose curve and simulant discharge curve whether greater than first threshold, then finishes to mate flow process when matching degree is not more than first threshold;
Adjusting module is used for during greater than first threshold, adjusting the TPS parameter in described matching degree, recomputates the simulant discharge curve according to the TPS parameter after adjusting, and will recomputate the simulant discharge curve that obtains and be sent to matching module and re-start coupling;
MBM is used for setting up dose modal according to the parameter after adjusting.
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 set up mathematical model realization program and automatically adjust the TPS parameter, avoid the complicated processes of artificial adjustment TPS parameter, also got rid of the subjectivity of artificial judgment matching degree; Efficient and the accuracy of TPS parameter adjustment have been improved.
Accordingly, can also system embodiment shown in Figure 4 be optimized for the various prioritization schemes in the described embodiment of the method, specifically can have following:
Described matching module specifically comprises:
Function unit is used for pre-establishing and preserving the matching degree function;
The first computing unit is for the measurement dose curve, matching degree function and the simulant discharge curve that obtain accelerator; To measure in the coordinate figure substitution matching degree function of dose curve and simulant discharge curve, with the value of the matching degree function that the obtains matching degree as measurement dose curve and simulant discharge curve.
Described adjusting module is specially:
The coordinate selection unit is used for finding and measures dose curve with respect to the point coordinate of simulant discharge curve side-play amount maximum;
Parameter selection unit is used for the coordinate figure according to this point coordinate, selects to affect the TPS parameter of this point coordinate;
The parameter adjustment unit is used for adjusting selecteed TPS parameter;
The second computing unit is used for recomputating the simulant discharge curve according to the TPS parameter after adjusting;
Transmitting element is used for recomputating the simulant discharge curve that obtains and is sent to matching module and re-starts coupling.
Described adjusting module also comprises:
The weight allocation unit, being used to affect the parametric distribution of this point coordinate weight probability;
Then described parameter selection unit is according to the coordinate figure of this point coordinate, and the TPS parameter of selecting to affect this point coordinate is specially, according to a kind of parameter that may affect this point coordinate of the random selection of weight probability.
Described adjusting module also comprises:
The Second Threshold unit is used for arranging Second Threshold;
Then described parameter adjustment unit is adjusted selecteed TPS parameter and is specially, and selected TPS parameter is increased or deducts Second Threshold.
The above only is preferred implementation of the present invention; should be pointed out that for those skilled in the art, under the prerequisite that does not break away from the principle 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; Calculation of parameter according to TPS goes out to simulate dose curve;
The matching degree of b, computation and measurement dose curve and simulant discharge curve, and judge that whether this matching degree is greater than first threshold;
If the c matching degree is then adjusted the TPS parameter greater than first threshold, recomputate the simulant discharge curve according to the TPS parameter after adjusting;
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 adjusting.
2. described method according to claim 1 is characterized in that, the matching degree of described computation and measurement dose curve and simulant discharge curve is specially:
Pre-establish the matching degree function, will measure in the coordinate figure substitution matching degree function of dose curve and simulant discharge curve, with the value of the matching degree function that the obtains matching degree as measurement dose curve and simulant discharge curve.
3. described method according to claim 2 is characterized in that, the described matching degree function that pre-establishes is specially:
Get the quadratic sum of measurement dose curve and simulant discharge curve corresponding point coordinate difference as the value of matching degree function, then the expression formula of described matching degree function is;
f = Σ k [ C ( k ) - M ( k ) ] 2
K is variable, the expression abscissa value, and horizontal ordinate is the coordinate figure of the point of k on C (k) the representative simulation dose curve, horizontal ordinate is the coordinate figure of the point of k on M (k) the representative measurement dose curve; F is the value of matching degree function.
4. described method according to claim 2 is characterized in that, the described matching degree function that pre-establishes is specially:
Get the quadratic sum of measuring part corresponding point coordinate difference in dose curve and the simulant discharge curve open country; The measurement dose curve of weight and the quadratic sum of simulant discharge curve penumbra region part corresponding point slope differences are mated in addition again, as the value of matching degree function; Then the expression formula of described matching degree function is,
f = Σ k [ C A ( k ) - M A ( k ) ] 2 + ω Σ k [ T C ( k + Δk ) - T M ( k ) ] 2
K is variable, the expression abscissa value; C A(k) the part horizontal ordinate is the coordinate figure of the point of k in the representative simulation dose curve u'eno; M A(k) the part horizontal ordinate is the coordinate figure of the point of k in the representative measurement dose curve u'eno, and ω represents weight coefficient, T CHorizontal ordinate is the slope of the point of k+ △ k on (k+ △ k) representative simulation dose curve; T M(k) horizontal ordinate is the slope of the point of k on the representative measurement dose curve; △ k represents the side-play amount of two curve penumbra region parts.
5. described method according to claim 1 is characterized in that, described adjustment TPS parameter specifically may further comprise the steps:
Find and measure dose curve with respect to the point coordinate of simulant discharge curve side-play amount maximum;
According to the coordinate figure of this point coordinate, the TPS parameter of selecting to affect this point coordinate is adjusted.
6. described method according to claim 5 is characterized in that, the TPS parameter that described selection affects this point coordinate is specially:
For might affecting the parametric distribution weight probability of this point coordinate, and select a kind of parameter that may affect this point coordinate according to the weight probability is random.
7. the described method of any one is characterized in that according to claim 1-6, and described adjustment TPS parameter is specially:
Set in advance Second Threshold, and make the increase of TPS parameter or deduct Second Threshold.
8. a dose modal matching system is characterized in that, described system specifically comprises:
Analog module is used for going out to simulate dose curve according to the calculation of parameter of TPS;
Matching module is for the measurement dose curve and the described simulant discharge curve that obtain accelerator, the matching degree of computation and measurement dose curve and simulant discharge curve;
Judge module is used for judging that the matching degree of measuring dose curve and simulant discharge curve whether greater than first threshold, then finishes to mate flow process when matching degree is not more than first threshold;
Adjusting module is used for during greater than first threshold, adjusting the TPS parameter in described matching degree, recomputates the simulant discharge curve according to the TPS parameter after adjusting, and will recomputate the simulant discharge curve that obtains and be sent to matching module and re-start coupling;
MBM is used for setting up dose modal according to the parameter after adjusting.
9. described system according to claim 8 is characterized in that, described matching module specifically comprises:
Function unit is used for pre-establishing and preserving the matching degree function;
The first computing unit is for the measurement dose curve, matching degree function and the simulant discharge curve that obtain accelerator; To measure in the coordinate figure substitution matching degree function of dose curve and simulant discharge curve, with the value of the matching degree function that the obtains matching degree as measurement dose curve and simulant discharge curve.
10. described system according to claim 8 is characterized in that, described adjusting module is specially:
The coordinate selection unit is used for finding and measures dose curve with respect to the point coordinate of simulant discharge curve side-play amount maximum;
Parameter selection unit is used for the coordinate figure according to this point coordinate, selects to affect the TPS parameter of this point coordinate;
The parameter adjustment unit is used for adjusting selecteed TPS parameter;
The second computing unit is used for recomputating the simulant discharge curve according to the TPS parameter after adjusting;
Transmitting element is used for recomputating the simulant discharge curve that obtains and is sent to matching module and re-starts coupling.
11. described system is characterized in that according to claim 10, described adjusting module also comprises:
The weight allocation unit, being used to affect the parametric distribution of this point coordinate weight probability;
Then described parameter selection unit is according to the coordinate figure of this point coordinate, and the TPS parameter of selecting to affect this point coordinate is specially, according to a kind of parameter that may affect this point coordinate of the random selection of weight probability.
12. according to claim 10 or 11 described systems, it is characterized in that, described adjusting module also comprises:
The Second Threshold unit is used for arranging Second Threshold;
Then described parameter adjustment unit is adjusted selecteed TPS parameter and is specially, and selected TPS parameter is increased or deducts Second Threshold.
CN201210450646.0A 2012-11-12 2012-11-12 Dose model matching method and system Active CN102930109B (en)

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 true CN102930109A (en) 2013-02-13
CN102930109B CN102930109B (en) 2015-07-01

Family

ID=47644906

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210450646.0A Active CN102930109B (en) 2012-11-12 2012-11-12 Dose model matching method and system

Country Status (1)

Country Link
CN (1) CN102930109B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103565459A (en) * 2013-09-25 2014-02-12 沈阳东软医疗系统有限公司 Method and device for judging imaging quality
WO2018095346A1 (en) * 2016-11-25 2018-05-31 平李⋅斯图尔特 Medical imaging system based on hmds
CN108573751A (en) * 2017-03-08 2018-09-25 深圳大森智能科技有限公司 The method of inspection and device of therapeutic effect

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1438602A (en) * 2003-03-17 2003-08-27 吴大可 Radiotherapy system of automatic computing prescription dosage

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1438602A (en) * 2003-03-17 2003-08-27 吴大可 Radiotherapy system of automatic computing prescription dosage

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GEORGE STARKSCHALL: "A constrained least-squares optimization method for external beam radiation therapy treatment planning", 《MEDICAL PHYSICS》 *
吴新根等: "有约束模拟退火法优化3D立体定向放射治疗计划", 《生物医药学工程学杂志》 *
黄禹: "适形调强放射治疗系统若干关键技术的研究", 《中国优秀博硕士学位论文全文数据库（博士）工程科技Ⅱ辑》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103565459A (en) * 2013-09-25 2014-02-12 沈阳东软医疗系统有限公司 Method and device for judging imaging quality
CN103565459B (en) * 2013-09-25 2015-07-01 沈阳东软医疗系统有限公司 Method and device for judging imaging quality
WO2018095346A1 (en) * 2016-11-25 2018-05-31 平李⋅斯图尔特 Medical imaging system based on hmds
AU2017366502B2 (en) * 2016-11-25 2021-05-20 David Wei Lee Medical imaging system based on HMDS
CN108573751A (en) * 2017-03-08 2018-09-25 深圳大森智能科技有限公司 The method of inspection and device of therapeutic effect

Also Published As

Publication number Publication date
CN102930109B (en) 2015-07-01

Similar Documents

Publication Publication Date Title
CN104338240B (en) Automatic optimization device for on-line self-adaption radiotherapy plan
JP6783398B2 (en) Systems and methods for modeling dose calculations in radiation therapy planning
US10076673B2 (en) Interactive dose gradient based optimization technique to control IMRT delivery complexity
CN112581475B (en) Method for predicting gamma passing rate of radiotherapy plan
CN105447330B (en) The weight regulating method and device of intensity-modulated radiation therapy
CN102921115A (en) Method for establishing measurement data-based simple and convenient irradiation source model of medical linear accelerator
CN106029170A (en) Incremental treatment planning
CN107644370A (en) Price competing method and system are brought in a kind of self-reinforcing study together
CN107551411B (en) Proton heavy ion intensity modulated radiotherapy robust optimization method aiming at range uncertainty
CN106228403A (en) A kind of user based on step analysis algorithm is worth methods of marking and system
CN102930109B (en) Dose model matching method and system
CN101908172A (en) Power market hybrid simulation method adopting various intelligent agent algorithms
CN108122079A (en) Computational methods, system and the storage medium of thermal power plant's sharing of load
CN109979564A (en) A kind of intelligence radiotherapy planning method, equipment and storage medium
US20140275703A1 (en) Method and system for dose determination of radiation therapy
CN111784203B (en) Electric power spot market risk simulation analysis method suitable for generator set participation
US20230385490A1 (en) Hydrological model considering uncertainty of runoff production structure and method for quantifying its impact on surface-subsurface hydrological process
CN108671419A (en) A kind of calculating of quick afterloading radiotherapy dosage, verification method
WO2022088341A1 (en) Radiation therapy control method and apparatus
CN107944552A (en) A kind of industrial Internet of Things parameter prediction method based on Elman neutral nets
CN114996715A (en) Game vulnerability intelligent repairing method based on deep learning
Petrasek et al. A Monte Carlo methodology for solving the optimal timber harvest problem with stochastic timber and carbon prices
Mesgari et al. Modeling the spatial distribution of crop cultivated areas at a large regional scale combining system dynamics and a modified Dyna-CLUE: A case from Iran
CN109999365B (en) Computing device and storage medium for implementing radiotherapy reverse optimization method
CN109447574A (en) Assets based on Fuzzy Optimum Neural Network turn solid project processing 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
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 Zhi Rui Radiotherapy Technology Co., Ltd.

Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16

Patentee before: Dongruan Medical Systems Co., Ltd., Shenyang