CN105664376B - A kind of device and equipment of dosage distribution determination and radiotherapy treatment planning optimization - Google Patents
A kind of device and equipment of dosage distribution determination and radiotherapy treatment planning optimization Download PDFInfo
- Publication number
- CN105664376B CN105664376B CN201511031048.XA CN201511031048A CN105664376B CN 105664376 B CN105664376 B CN 105664376B CN 201511031048 A CN201511031048 A CN 201511031048A CN 105664376 B CN105664376 B CN 105664376B
- Authority
- CN
- China
- Prior art keywords
- pen
- voxel
- bundle element
- dose
- continuously distributed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1042—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head
- A61N5/1045—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head using a multi-leaf collimator, e.g. for intensity modulated radiation therapy or IMRT
- A61N5/1047—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head using a multi-leaf collimator, e.g. for intensity modulated radiation therapy or IMRT with movement of the radiation head during application of radiation, e.g. for intensity modulated arc therapy or IMAT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
- A61N5/1036—Leaf sequencing algorithms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1071—Monitoring, verifying, controlling systems and methods for verifying the dose delivered by the treatment plan
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
This disclosure relates to the determination method and apparatus of radiotherapy treatment planning middle dosage distribution, comprising: turn to multiple pen-bundle elements and multiple voxels for beam and area-of-interest are discrete respectively, determine each pen-bundle element to the dose contribution value of each voxel；Determine the integral dose on each voxel with the continuously distributed of pen-bundle element according to the dose contribution value；According to the continuously distributed dosage distribution calculated at current vane position.Present disclosure also relates to radiotherapy treatment planning optimization method and devices.By means of the present invention or device, the blade of multi-leaf optical grating can be moved to the position of any optimal solution, be not only restricted to the resolution ratio of pen-bundle element, and method or apparatus of the invention, calculation amount is small, and speed is fast.
Description
Technical field
Present application relates generally to field of radiation therapy more particularly to radiotherapy treatment planning optimization method and equipment, Yi Jifang
Penetrate the determination method and apparatus for the treatment of plan middle dosage distribution.
Background technique
Multi-diaphragm collimator (MLC) is the equipment for generating the radiation field size of suitable shape, is widely used in medical domain.
Multi-diaphragm collimator usually reaches launched field dynamic or static forming by the movement of each blade.It is real in radiotherapy treatment planning optimization
The key of the function of existing multi-diaphragm collimator includes the determination of Ziye shape and intensity.
In the prior art, when determining Ziye shape and intensity, there are following challenges.
The continuous moving of blade is realized difficult in radiotherapy optimization problem.When formulating reverse optimization plan, need
The shift position of MLC blade in each Ziye shape is provided, usually by piecing together some pen-bundle elements (beamlet) dose matrix
To realize.At this moment blade can only navigate to the edge of pen-bundle element, and the mobile minimum range of blade is a pen-bundle element, therefore leaf
The mobile continuity of piece is limited by the resolution ratio of pen-bundle element；Resolution ratio is too low, and the minimum range for causing blade mobile is excessive, from
And the controllability of blade shift position is reduced, expected effect of optimization is not achieved；And resolution ratio is excessively high, then can make beam
The dose matrix data volume of member increased dramatically, and extend the calculating time significantly.
In addition, the data volume of pen-bundle element dose matrix is excessive.In VMAT planning optimization process, first have to calculate pen-bundle element
Dose matrix is used to store the pen-bundle element pen-bundle element dosage deposition (cGy/ of sampled point or voxel (voxel) nearby
MU).In order to preferably control the distribution of the dosage on each organ in some multiple organ cases, the number selection of voxel is relatively more,
This allows for that mass data can be generated in the lower situation of the resolution ratio in pen-bundle element, causes to consume a large amount of memory money
Source, and occupy and calculate the time.In addition, calculating speed is also one of problem of the prior art slowly.In the optimization of radiotherapy treatment planning
In the process, need to calculate the dosage distribution situation on current voxel, and according to present dose distribution situation, to determine in next step
Blade moving direction；Each more new dosage distribution situation or change leaf position require to do accumulative place to Ziye dosage
Reason reduces optimal speed so that calculation amount be made to increase.
For example, Jiang, Men and Jia are in 2011 in the U.S. Patent application that number of patent application is US 13/695,015
In mention, using column-generation (column generation) method, solve Ziye shape and strong respectively in two problems
Degree.This method has certain effect for solving Large-scale Optimization Problems, but this method variable space is big, consumes compared with multicomputer
Resource, and Segment weight and shape are separately optimized in optimization process, it is not that local optimum is found in global scope
Solution.
For another example, it is added by building network model, and by the mobile constraint of blade in United States Patent (USP) US 2013/0077751A1
In, which more focuses on the mobile continuity of blade, and ignores the effect of optimization relatively, and its Ziye shape and Segment weight
And be separately optimized.
In above two patents, the moveable position of blade be all it is discrete, moveable position is still limited by beam
Member, therefore its intense adjustment ability will will receive limitation largely.
Summary of the invention
According to the one side of the disclosure, a kind of determination method of radiotherapy treatment planning middle dosage distribution includes by beam and sense
Interest region is discrete respectively to turn to multiple pen-bundle elements and multiple voxels, determines each pen-bundle element to the dose contribution of each voxel
Value；Determine the integral dose on each voxel with the continuously distributed of pen-bundle element according to the dose contribution value；And according to described
The continuously distributed dosage distribution calculated at current vane position.
According to an embodiment of this aspect, determine the integral dose on each voxel with beam according to the dose contribution value
The continuously distributed of member further comprises depositing dose contribution value of each pen-bundle element to the voxel by row to each voxel
Storage, the row are parallel with the moving direction of blade；To the pen-bundle element of every row, calculate after each pen-bundle element is sequentially overlapped, the body
Integral dose on element with pen-bundle element distribution；And the integral dose is carried out with the distribution of pen-bundle element using fitting function
Fitting is to obtain the integral dose on the voxel with the continuously distributed of pen-bundle element.The fitting function can be one-dimensional functions,
The integral dose is integral dose with the curve of change in location with the continuously distributed of pen-bundle element.Alternatively, the fitting letter
Number can be two-dimensional function, and the integral dose is integral dose with the song of position and line number variation with the continuously distributed of pen-bundle element
Face.
According to this aspect further embodiment, the integral dose on all voxels is obtained with the continuous of pen-bundle element in fitting
It further include being modified to the integral dose being fitted on the obtained voxel with the continuously distributed of pen-bundle element after distribution.Root
According to another embodiment, the integral dose on the voxel obtained using interpolation method to fitting is with the continuously distributed of pen-bundle element
It is modified.According to further embodiments, the integral dose on the voxel obtained to fitting divides with the continuous of pen-bundle element
Cloth, which is modified, further comprises: the Optimized model of radiotherapy treatment planning is established, it is tired on each voxel obtained using fitting
Dosage is counted with the initial solution of the continuously distributed calculating Optimized model of pen-bundle element, according to the initial solution location to each
Integral dose on voxel is modified with the continuously distributed of pen-bundle element.Further, according to the initial solution location to every
Integral dose on a voxel with pen-bundle element it is continuously distributed be modified may include to each voxel, it is true according to the initial solution
Surely the position range for needing modified pen-bundle element, using interpolation method to accumulative in the position range for needing modified pen-bundle element
Dosage is modified.
According to another aspect of the present disclosure, a kind of radiotherapy treatment planning optimization method includes by beam and area-of-interest point
Taking leave of dispersion is multiple pen-bundle elements and multiple voxels, determines each pen-bundle element to the dose contribution value of each voxel；According to described
Dose contribution value determines the integral dose on each voxel with the continuously distributed of pen-bundle element；Establish the optimization mould of radiotherapy treatment planning
Type；And using the integral dose on each voxel with the optimal solution of the continuously distributed calculating Optimized model of pen-bundle element, and
Optimum results as radiotherapy treatment planning.
According to an embodiment of this aspect, determine the integral dose on each voxel with beam according to the dose contribution value
The continuously distributed of member further comprises depositing dose contribution value of each pen-bundle element to the voxel by row to each voxel
Storage, the row are parallel with the moving direction of blade；To the pen-bundle element of every row, calculate after each pen-bundle element is sequentially overlapped, the body
Integral dose on element with pen-bundle element distribution；And the integral dose is carried out with the distribution of pen-bundle element using fitting function
Fitting is to obtain the integral dose on the voxel with the continuously distributed of pen-bundle element.The fitting function can be one-dimensional functions,
The integral dose is integral dose with the curve of change in location with the continuously distributed of pen-bundle element.Alternatively, the fitting letter
Number can be two-dimensional function, and the integral dose is integral dose with the song of position and line number variation with the continuously distributed of pen-bundle element
Face.
According to the another embodiment of this aspect, the integral dose on all voxels is obtained in fitting and is divided with the continuous of pen-bundle element
It further include being modified to the integral dose being fitted on the obtained voxel with the continuously distributed of pen-bundle element after cloth.For example,
The integral dose on the voxel obtained using interpolation method to fitting is modified with the continuously distributed of pen-bundle element.
According to further embodiments, the integral dose on the voxel obtained to fitting is with the continuously distributed of pen-bundle element
The integral dose on each voxel obtained including the use of fitting is modified with the continuously distributed calculating optimization of pen-bundle element
The initial solution of model, according to the initial solution location to the integral dose on each voxel with pen-bundle element it is continuously distributed into
Row amendment.
According to the another embodiment of this aspect, according to the initial solution location to the integral dose on each voxel with
The continuously distributed of pen-bundle element is modified including determining the position for needing modified pen-bundle element according to the initial solution to each voxel
Range is set, the integral dose in the position range for needing modified pen-bundle element is modified using interpolation method.Further
Embodiment may also include the continuously distributed calculating optimal solution feelings according to the integral dose on the revised voxel with pen-bundle element
Dosage distribution under condition.
The disclosure further includes equipment corresponding with above-mentioned various aspects.
Detailed description of the invention
Fig. 1 is the example being distributed according to the pen-bundle element dosage of one embodiment of the disclosure.
Fig. 2 is the example being distributed according to the pen-bundle element dose accumulation of one embodiment of the disclosure.
Fig. 3 is the example according to the beam variable matrix uniline fitting result of one embodiment of the disclosure.
Fig. 4 is the result example that piecewise interpolation is carried out to accumulative pen-bundle element section according to one embodiment of the disclosure.
Specific embodiment
Various aspects will now be described with reference to the drawings.In the following description, numerous details are elaborated for explanatory purposes
To provide the thorough understanding to one or more aspects.It will, however, be evident that can also be practiced without these details such
In terms of (class).
This disclosure relates to determination method and apparatus, the radiotherapy treatment planning optimization method of the distribution of radiotherapy treatment planning middle dosage
And equipment.For example, all aspects of this disclosure can relate to, but are not limited to it is following one or more: by the method for curve matching come
Construct blade can continuous moving Optimized model；Beam variable matrix is replaced with resulting parameter is fitted, to greatly reduce
Data scale；The continuously distributed of integral dose is determined with the method for fitting, so that the flatness of optimization problem is enhanced, so that
Ideal solution can be converged to, the resolution ratio of pen-bundle element is not only restricted to；Optimal solution is quickly found using analytic method, for quasi-
It closes error to correct using such as interpolation method, and calculating speed and precision can have been taken into account based on amendment come re-optimization.
The various specific embodiments of the disclosure are given below.It will be appreciated that the specific embodiment of the disclosure is to show
Example and what infinite mode provided, the protection scope of the disclosure is not limited in any specific embodiment.Specifically, this public affairs
Claimed range is opened to be limited only by the claims.For example, this disclosure relates to but being not limited to multi-diaphragm collimator blade
Position Filled function.As one of ordinary skill understandable, the disclosure is equally adaptable to other related applications.
According to the aspect of the disclosure, multiple pen-bundle elements can be turned to by beam and area-of-interest are discrete respectively
(beamlet) and multiple voxels (voxel), determine each pen-bundle element to the dose contribution value of each voxel；And according to described
Dose contribution value determines the integral dose on each voxel with the continuously distributed of pen-bundle element.
For example, determining that each pen-bundle element may include calculating pen-bundle element dose matrix to the dose contribution value of each voxel.It can
Further by calculated pen-bundle element dose matrix by row storage.Skilled person will appreciate that, it is possible to use other storage sides
Formula, by column storage etc., without departing from the scope of the present disclosure.According to an embodiment, the capable moving direction with blade can be
Parallel.
Fig. 1 shows the dosage distribution situation of each pen-bundle element of a voxel in a row according to an exemplary embodiment,
Wherein x-axis can be the position of the pen-bundle element at the row from left to right, and y-axis can be pen-bundle element for the dose contribution value of the voxel
(cGy/MU), and it can remember that the corresponding dose contribution value of i-th of pen-bundle element is bi。
According to further exemplary and non-limiting embodiment, the optimization method of the disclosure may include from left to right right
The dose contribution value of pen-bundle element adds up, to obtain the distribution of the integral dose on the voxel.Such as those of ordinary skill in the art
It will be appreciated, the disclosure is not limited to add up from left to right, but can have other cumulative directions and/or side
Formula.For example, to the pen-bundle element of every row, can calculate each pen-bundle element dose contribution value be sequentially overlapped after, it is tired on the voxel
Dosage is counted with the distribution of pen-bundle element position.
According to an example, which can be as shown in Fig. 2, wherein x-axis be pen-bundle element position from left to right, y-axis
For corresponding position pen-bundle element dose contribution value it is cumulative after dosage aggregate-value.For example, yiIt can be the agent of i-th of pen-bundle element
The dosage aggregate-value on the voxel after contribution margin adds up is measured, can for example be calculated as shown in following formula (1):
yi=yi-1+bi (1)
According to further exemplary and non-limiting embodiment, according to the continuously distributed calculating current vane position
The dosage distribution at place may further include to be fitted to obtain to integral dose using fitting function with the distribution of pen-bundle element
Integral dose on to the voxel is with the continuously distributed of pen-bundle element, such as may include choosing fitting function to tired on the voxel
Meter dosage distribution is fitted.According to non-limiting example, S type fitting function can be chosen, the integral dose on the voxel is distributed
It is fitted.The choosing method of the fitting function has very much, the non-limiting example as exemplary, than if any Logistics
Curve, Gompertz curve, Bertalanffy curve, bent tangent function (tanh (x)), polynomial function etc..
In the following, for convenient for explaining, in conjunction with being illustrated for Logistics function to embodiment of the disclosure.Using it
His curve can be ready-made for those of ordinary skill in the art on the basis of the disclosure come the example being fitted, so by including
Within the scope of this invention.
For example, typical Logistics function can be written as:
Wherein a, b and K can be undetermined parameter.
For example, in order to be fitted using such Logistics function to the integral dose distribution on such as voxel, it can
By being adjusted such that integral dose distribution and Logistics function to any one or more of undetermined parameter a, b and K
Reach best fit.Target it is, for example, possible to use least mean-square error as fitting function.For example, mean square error minimizes
It can be written as:
Wherein N represents the number of pen-bundle element, by Optimum Theory it is found that then existing for example, by seeking local derviation to a, b, K respectively
Following equation group (4) a, b, K value corresponding when setting up is the optimal solution of the optimization problem.
Those of ordinary skill in the art can correspondingly adjust its parameter it is found that if having selected other fitting functions
It saves so that integral dose distribution function and selected fitting function reach best fit.
It is as shown in Figure 3 to the fitting result of beam variable matrix at uniline according to an exemplary embodiment.For each voxel
For, N number of data in beam variable matrix on original every row can be replaced with 3 parameters K, a, b after fitting, thus
The data volume having shared by beam variable matrix can be greatly reduced, reduce storage and/or transmit occupied resource.Not only such as
This, can also further accelerate the calculating of objective function and/or variable derivative using analytic method.
In upper example, fitting function can be one-dimensional functions, and integral dose is accumulative with the continuously distributed of pen-bundle element
Dosage with change in location curve.In another embodiment, which can be two-dimensional function, and integral dose is with pen-bundle element
The continuously distributed curved surface changed for integral dose with position and line number.
For example, in an exemplary and non-limiting example, for example, still by taking Logistics function as an example, it is two-dimensional quasi-
Closing function can be written as:
Wherein x represents leaf position, and r indicates the row where blade.The available further compression of data volume in this way,
Cost is to might have corresponding precision to reduce.
Determined the integral dose on each voxel with pen-bundle element it is continuously distributed after, on the one hand, can be according to described
The continuously distributed dosage distribution calculated at current vane position.
According to another aspect of the present disclosure, determined the integral dose on each voxel with pen-bundle element according to dose contribution value
It is continuously distributed after, can further establish the Optimized model of radiotherapy treatment planning, and utilize the integral dose on each voxel
With the optimal solution of the continuously distributed calculating Optimized model of pen-bundle element, and the optimum results as radiotherapy treatment planning.
According to the one side of the disclosure, a kind of analytic method of Optimized model that establishing radiotherapy treatment planning may include set it is excellent
The catalogue scalar functions of change are F, then F is the function about Ziye leaf position x and Segment weight w:
F=F (x, w)=F (d (x, w)) (5)
Wherein d is djThe vector of composition, djFor the dose value at j-th of voxel (or sampled point), expression formula can be
Such as:
K is the number of Ziye；wkFor the weight of k-th of Ziye, the range of k is from 1 to K；R is the line number of MLC blade；R is
Row, in the range of R；WithThe position of the left and right blade of respectively k-th Ziye r row；For k-th of Ziye r
The fitting function of j-th of voxel of row.
Then above formula need to only be substituted by calculating the dosage distribution situation at current vane position, without again to beam
The value of first dose matrix adds up.As a result, objective function can be calculated with accelerate.Accumulative agent is determined with the method for fitting
The continuously distributed flatness for enhancing optimization problem of amount, allows to converge to ideal solution.
On the other hand, using based on gradient optimization method carry out radiotherapy treatment planning when it needs to be determined that each variable is led
Number, can also very easily be acquired using above formula:
(7)
Wherein
Wherein flag is mark, works as xk,rIt is flag=-1, x for lobus sinister piecek,rFlag=1 when for lobus dexter piece.
Thus, it is possible to the derivative at current variable is readily available, for the such as, but not limited to optimization based on gradient
Method etc..
There may be certain errors for the dosage summation curve and truth obtained for example, by above-mentioned fitting.Therefore
After obtaining preliminary solution using analytic method as escribed above, need to be modified using accurate pen-bundle element dosage.As example
Property and not limiting example, using various interpolation methods to the integral dose on the obtained voxel of fitting with the continuous of pen-bundle element
Distribution is modified.It is, for example, possible to use polynomial interopolation, such as Lagrangian (Lagrange) interpolation, newton (Newton)
Interpolation；Hermit interpolation；Piecewise interpolation；Spline interpolation etc..
Accoding to exemplary embodiment, it after establishing the Optimized model of radiotherapy treatment planning, can use for example, by aforementioned
Integral dose on obtained each voxel is fitted with the initial solution of the continuously distributed calculating Optimized model of pen-bundle element, then root
The integral dose on each voxel is modified with the continuously distributed of pen-bundle element according to the initial solution location.
According to an exemplary embodiment, for example, accurately modified penetrate can be needed to determine first according to current vane position
Beam member position range.Such as to each voxel, modified beam can be needed to determine according to based on the obtained initial solution of fitting
First position range.Then, the integral dose in the modified pen-bundle element position range of needs is modified using interpolation method.
According to an exemplary embodiment, interpolation can be carried out to the section, so as to obtain one using each pen-bundle element position as node
By the exact dose change curve of each pen-bundle element integral dose value.
In an exemplary and non-limiting example, it is modified using hermit interpolation.It is exemplary in Fig. 4
Ground is shown centered on 300, to section 260-360 using arriving after 3 hermit interpolation of segmentation as a result, wherein x-axis
Unit can be millimeter (mm).
As from institute in Fig. 4 as it can be seen that being wherein true pen-bundle element position at ' * ', solid line is the result after interpolation.Utilize the party
After method amendment, exact dose of the available blade in specific location.
According to a further exemplary embodiment, the optimization process of aforementioned analytic method can be repeated based on this, from
And preliminary result can be corrected.For example, can be according to the integral dose on the revised voxel with the company of pen-bundle element
Continuous distribution, to calculate the dosage distribution in the case of optimal solution.Optimal solution quickly is found by using fitting process, fitting is missed
Difference is corrected using such as interpolation method, and has taken into account calculating speed and precision based on amendment come re-optimization.
Summarized according to the process of the exemplary and non-limiting embodiment of the disclosure one can include:
(1) pen-bundle element is calculated, dosage is deposited by each voxel and is for example added up by row, accumulative agent discrete on the row is obtained
Saltation point is measured, the point on these positions is fitted using such as ' S ' type curve, the parameter after record fitting is based on subsequent
It calculates；
(2) using for example based on the optimal method of gradient, the optimal solution of solving optimization model, solve dosage it is accumulative and
It can use existing fitting parameter when derivative quickly to calculate；
(3) point centered on leaf position obtained in (2) is chosen a small section as correcting region, is utilized
Interpolation method carries out local correction to the optimum results in (2), obtains exact dose distribution situation, and correct locally optimal solution.
It similarly, equally can be with class after being optimized using for example foregoing two-dimentional fitting function with analytic method
As be modified using interpolation method.
Realize various advantages according to disclosed method and equipment, including but not limited to for example it is following one or more:
1. reducing the data volume of pen-bundle element.Pen-bundle element agent is replaced with fitting parameter according to disclosed method and equipment
Moment matrix is substantially reduced data volume, to reduce the required stock number occupied of storage/transmission and greatly improve the efficiency.
2. utilizing fitting parameter, calculated locating for dosage distribution situation and each blade under current shape by analytic method
The derivative of position greatly accelerates the iteration speed of optimization compared with dosage stacking method.
3. the flatness of increased Optimized model, accelerates the convergence rate of optimization problem and to approach problem more optimal
Solution.
4. using two step optimisation strategies, the first step quickly orients the rough position of locally optimal solution using approximating method
It sets；Second step makes up the error in approximating method using interpolation method, accurately corrects locally optimal solution.Do not losing optimization precision
Under the premise of, increase substantially optimal speed.
It will be recognized by one of ordinary skill in the art that the beneficial effect of the disclosure is not by any single embodiment Lai all real
It is existing.Various combinations, modification and replacement are that those of ordinary skill in the art are illustrated on the basis of the disclosure.
In addition, term "or" is intended to indicate that inclusive "or" and nonexcludability "or".That is, unless otherwise specified or from upper and lower
Text can be clearly seen, otherwise phrase " X uses A or B " be intended to indicate that it is any naturally can and arrangement.That is, phrase " X using A or
B " is met by example any in following instance: X uses A；X uses B；Or X uses both A and B.In addition, the application and
The article " one " used in the appended claims and " certain " are generally to be understood as indicating " one or more ", unless in addition
It states or can be apparent from from the context to refer to singular.
Various aspects or feature by by may include several equipment, component, module, and the like system in the form of be in
It is existing.It should be understood that and understand, various systems may include optional equipment, component, module etc., and/or can not include in conjunction with attached drawing
Armamentarium, component, module for being discussed etc..Also the combination of these methods can be used.
It can be with general in conjunction with various illustrative logicals, logical block, module and the circuit that presently disclosed embodiment describes
Processor, digital signal processor (DSP), specific integrated circuit (ASIC), field programmable gate array (FPGA) or it is other can
Programmed logic device, discrete door or transistor logic, discrete hardware component or its be designed to carry out function described herein
Any combination realize or execute.General processor can be microprocessor, but in alternative, and processor, which can be, appoints
What conventional processor, controller, microcontroller or state machine.Processor is also implemented as calculating the combination of equipment, example
Such as DSP and the combination of microprocessor, multi-microprocessor, one or more microprocessors cooperateed with DSP core or any
Other such configurations.In addition, at least one processor may include may act on execute one or more steps described above and/
Or one or more modules of movement.
In addition, the method that describes in conjunction with aspect disclosed herein or the step of algorithm and/or movement can be directly hard
Implement in part, in the software module executed by processor or in combination of the two.
The element of the various aspects described in the whole text in the disclosure is that those of ordinary skill in the art are currently or hereafter known
Equivalent scheme in all structures and functionally is clearly included in this by citation, and is intended to be intended to be encompassed by the claims.
In addition, any content disclosed herein is all not intended to contribute to the public --- it is no matter such open whether in claim
It is explicitly recited in book.
Claims (19)
1. a kind of for determining the device of radiotherapy treatment planning middle dosage distribution characterized by comprising
Memory；And
Processor is coupled to the memory and is configured to:
Multiple pen-bundle elements and multiple voxels are turned to by beam and area-of-interest are discrete respectively, determine each pen-bundle element to every individual
The dose contribution value of element；
Determine the integral dose on each voxel with the continuously distributed of pen-bundle element according to the dose contribution value；And
According to the continuously distributed dosage distribution calculated at current vane position, wherein determined according to the dose contribution value
Integral dose on each voxel further comprises to integral dose using fitting function with beam with the continuously distributed of pen-bundle element
The distribution of member is fitted to obtain the integral dose on the voxel with the continuously distributed of pen-bundle element.
2. device as described in claim 1, which is characterized in that determined according to the dose contribution value accumulative on each voxel
Dosage further comprises with the continuously distributed of pen-bundle element:
To each voxel, by each pen-bundle element to the dose contribution value of the voxel by row storage, the row and blade
Moving direction is parallel；
To the pen-bundle element of every row, calculate after each pen-bundle element is sequentially overlapped, the integral dose on the voxel with pen-bundle element point
Cloth.
3. device as claimed in claim 2, which is characterized in that the fitting function be one-dimensional functions, the integral dose with
The continuously distributed of pen-bundle element is integral dose with the curve of change in location.
4. device as claimed in claim 2, which is characterized in that the fitting function be two-dimensional function, the integral dose with
The continuously distributed curved surface changed for integral dose with position and line number of pen-bundle element.
5. device as claimed in claim 2, which is characterized in that obtain the integral dose on all voxels with pen-bundle element in fitting
It is continuously distributed after, the processor is further configured to,
The integral dose on the voxel obtained to fitting is modified with the continuously distributed of pen-bundle element.
6. device as claimed in claim 5, which is characterized in that tired on the voxel obtained using interpolation method to fitting
Meter dosage is modified with the continuously distributed of pen-bundle element.
7. device as claimed in claim 5, which is characterized in that the described pair of integral dose being fitted on the obtained voxel with
Continuously distributed be modified of pen-bundle element further comprises:
The Optimized model of radiotherapy treatment planning is established,
Using the integral dose on the obtained each voxel of fitting with pen-bundle element the continuously distributed calculating Optimized model just
Begin solution,
The integral dose on each voxel is modified with the continuously distributed of pen-bundle element according to the initial solution location.
8. device as claimed in claim 7, which is characterized in that according to the initial solution location to tired on each voxel
Meter dosage is modified with the continuously distributed of pen-bundle element, including
To each voxel, the position range for needing modified pen-bundle element is determined according to the initial solution,
The integral dose in the modified pen-bundle element position range of needs is modified using interpolation method.
9. a kind of device for radiotherapy treatment planning optimization characterized by comprising
Memory；And
Processor is coupled to the memory and is configured to:
Multiple pen-bundle elements and multiple voxels are turned to by beam and area-of-interest are discrete respectively, determine each pen-bundle element to every individual
The dose contribution value of element；
Determine the integral dose on each voxel with the continuously distributed of pen-bundle element according to the dose contribution value；
Establish the Optimized model of radiotherapy treatment planning；And
Using the integral dose on each voxel with the optimal solution of the continuously distributed calculating Optimized model of pen-bundle element, and conduct
The optimum results of radiotherapy treatment planning, wherein determine the integral dose on each voxel with beam according to the dose contribution value
The continuously distributed of member further comprises being fitted to obtain institute with the distribution of pen-bundle element to integral dose using fitting function
The integral dose on voxel is stated with the continuously distributed of pen-bundle element.
10. device as claimed in claim 9, which is characterized in that determined according to the dose contribution value tired on each voxel
Dosage, which is counted, with the continuously distributed of pen-bundle element further comprises:
To each voxel, by each pen-bundle element to the dose contribution value of the voxel by row storage, the row and blade
Moving direction is parallel；
To the pen-bundle element of every row, calculate after each pen-bundle element is sequentially overlapped, the integral dose on the voxel with pen-bundle element point
Cloth.
11. device as claimed in claim 10, which is characterized in that the fitting function is one-dimensional functions, the integral dose
Continuously distributed with pen-bundle element is integral dose with the curve of change in location.
12. device as claimed in claim 10, which is characterized in that the fitting function is two-dimensional function, the integral dose
With the continuously distributed curved surface changed for integral dose with position and line number of pen-bundle element.
13. device as claimed in claim 10, which is characterized in that obtain the integral dose on all voxels with beam in fitting
Member it is continuously distributed after, further include,
The integral dose on the voxel obtained to fitting is modified with the continuously distributed of pen-bundle element.
14. device as claimed in claim 13, which is characterized in that on the voxel obtained using interpolation method to fitting
Integral dose is modified with the continuously distributed of pen-bundle element.
15. device as claimed in claim 13, which is characterized in that the described pair of integral dose being fitted on the obtained voxel
Continuously distributed with pen-bundle element is modified, including,
Using the integral dose on the obtained each voxel of fitting with pen-bundle element the continuously distributed calculating Optimized model just
Begin solution,
The integral dose on each voxel is modified with the continuously distributed of pen-bundle element according to the initial solution location.
16. device as claimed in claim 15, which is characterized in that according to the initial solution location on each voxel
Integral dose is modified with the continuously distributed of pen-bundle element, including
To each voxel, the position range for needing modified pen-bundle element is determined according to the initial solution,
The integral dose in the position range for needing modified pen-bundle element is modified using interpolation method.
17. device as claimed in claim 16, which is characterized in that the processor is further configured to, according to revised institute
The integral dose stated on voxel is distributed with the dosage in the case of the continuously distributed calculating optimal solution of pen-bundle element.
18. a kind of radiotherapy treatment planning middle dosage is distributed locking equipment really characterized by comprising
For turning to multiple pen-bundle elements and multiple voxels for beam and area-of-interest are discrete respectively, determine each pen-bundle element to every
The device of the dose contribution value of a voxel；
For determining the integral dose on each voxel with the continuously distributed device of pen-bundle element according to the dose contribution value；With
And
For the device according to the continuously distributed dosage distribution calculated at current vane position, wherein for according to
Dose contribution value determines that the integral dose on each voxel with the continuously distributed device of pen-bundle element further comprises for using
Fitting function is fitted to obtain the integral dose on the voxel with pen-bundle element with the distribution of pen-bundle element to integral dose
Continuously distributed device.
19. a kind of radiotherapy treatment planning optimizes equipment characterized by comprising
For turning to multiple pen-bundle elements and multiple voxels for beam and area-of-interest are discrete respectively, determine each pen-bundle element to every
The device of the dose contribution value of a voxel；
For determining the integral dose on each voxel with the continuously distributed device of pen-bundle element according to the dose contribution value；
For establishing the device of the Optimized model of radiotherapy treatment planning；And
For the optimal solution using the integral dose on each voxel with the continuously distributed calculating Optimized model of pen-bundle element, and
The device of optimum results as radiotherapy treatment planning, wherein for being determined on each voxel according to the dose contribution value
Integral dose further comprises for using fitting function to integral dose with pen-bundle element with the continuously distributed device of pen-bundle element
Distribution be fitted to obtaining the integral dose on the voxel with the continuously distributed device of pen-bundle element.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511031048.XA CN105664376B (en) | 2015-12-31 | 2015-12-31 | A kind of device and equipment of dosage distribution determination and radiotherapy treatment planning optimization |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511031048.XA CN105664376B (en) | 2015-12-31 | 2015-12-31 | A kind of device and equipment of dosage distribution determination and radiotherapy treatment planning optimization |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105664376A CN105664376A (en) | 2016-06-15 |
| CN105664376B true CN105664376B (en) | 2019-11-19 |
Family
ID=56189933
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511031048.XA Active CN105664376B (en) | 2015-12-31 | 2015-12-31 | A kind of device and equipment of dosage distribution determination and radiotherapy treatment planning optimization |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105664376B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018077709A1 (en) * | 2016-10-28 | 2018-05-03 | Koninklijke Philips N.V. | Graphical user interface for iterative treatment planning |
| EP3338859B1 (en) * | 2016-12-22 | 2019-06-26 | RaySearch Laboratories AB | System and method for modelling of scattering in ion radiotherapy treatment planning |
| CN107041999A (en) * | 2017-02-16 | 2017-08-15 | 西南医科大学附属医院 | There is 3-dimensional dose contribution distributed model construction method and its application of barrier material |
| EP3412340A1 (en) * | 2017-06-08 | 2018-12-12 | Koninklijke Philips N.V. | Treatment plan generation for radiation therapy treatment |
| CN108074636B (en) * | 2017-12-26 | 2021-10-26 | 北京纳米维景科技有限公司 | Surface incident dose calculation method, device and storage medium |
| JP7187972B2 (en) * | 2018-10-24 | 2022-12-13 | 富士通株式会社 | Apparatus, method and program |
| US10940330B2 (en) * | 2019-03-20 | 2021-03-09 | Varian Medical Systems International Ag | Dosimetric projection of collimator geometry |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011160235A1 (en) * | 2010-06-22 | 2011-12-29 | Karl Otto | System and method for estimating and manipulating estimated radiation dose |
| US9597529B2 (en) * | 2012-09-05 | 2017-03-21 | Best Medical International, Inc. | Rapid range stacking (RRS) for particle beam therapy |
| CN104645500B (en) * | 2015-02-12 | 2018-08-03 | 上海联影医疗科技有限公司 | A kind of intensity modulated radiation therapy optimization system |
-
2015
- 2015-12-31 CN CN201511031048.XA patent/CN105664376B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105664376A (en) | 2016-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105664376B (en) | A kind of device and equipment of dosage distribution determination and radiotherapy treatment planning optimization | |
| CN104318122B (en) | Disturbance-based intensity modulated arc therapy optimization method | |
| McQuaid et al. | IMRT delivery to a moving target by dynamic MLC tracking: delivery for targets moving in two dimensions in the beam's eye view | |
| US20180267399A1 (en) | Method of applying vertex based corrections to a semiconductor design | |
| CN105413068B (en) | Solve the optimization method and radiotherapy equipment of flux pattern | |
| Kamath et al. | Leaf sequencing algorithms for segmented multileaf collimation | |
| Ahunbay et al. | An on‐line replanning method for head and neck adaptive radiotherapy a | |
| CN110215623B (en) | Method and device for realizing orthogonal double-layer grating rotation intensity modulation based on subfield optimization | |
| CN106573152A (en) | Supervised 4-d dose map deformation for adaptive radiotherapy planning | |
| Jabbari et al. | Monte carlo simulation of siemens oncor linear accelerator with beamnrc and dosxyZnrc code | |
| Craft et al. | Deliverable navigation for multicriteria step and shoot IMRT treatment planning | |
| CN110570923B (en) | Mixed Monte Carlo radiotherapy reverse optimization method, equipment and storage medium | |
| FR2948030A1 (en) | METHOD FOR CALCULATING DOSES DEPOSITED BY IONIZING RADIATION | |
| CN110404185A (en) | A kind of strong dividing method of rotation tune and device of orthogonal double-layer grating | |
| CN105617536A (en) | Inverse intensity modulated arc radioterapy optimization method and device | |
| CN109499012A (en) | A kind of TPS automatic Iterative optimization algorithm of optimization dosage guidance | |
| Fuchs et al. | A pencil beam algorithm for helium ion beam therapy | |
| Childress et al. | Mobius3D white paper: dose calculation algorithm | |
| CN105617535A (en) | Dose distribution estimation method and sub-field optimization method | |
| CN105561485B (en) | The optimization device of radiotherapy treatment planning | |
| CN110248701A (en) | Thermal starting initialization to external beam radiation therapy planning optimization | |
| US20230214542A1 (en) | Static device for use in radiotherapy treatment and design method for such a device | |
| Yang et al. | Dose calculation for hypofractionated volumetric‐modulated arc therapy: approximating continuous arc delivery and tongue‐and‐groove modeling | |
| Saw et al. | Leaf sequencing techniques for MLC-based IMRT | |
| Moore et al. | Comparisons of treatment optimization directly incorporating random patient setup uncertainty with a margin‐based approach |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |