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 PDF

Info

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
Application number
CN201511031048.XA
Other languages
Chinese (zh)
Other versions
CN105664376A (en
Inventor
勾磐杰
Original Assignee
上海联影医疗科技有限公司
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 上海联影医疗科技有限公司 filed Critical 上海联影医疗科技有限公司
Priority to CN201511031048.XA priority Critical patent/CN105664376B/en
Publication of CN105664376A publication Critical patent/CN105664376A/en
Application granted granted Critical
Publication of CN105664376B publication Critical patent/CN105664376B/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1042X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head
    • A61N5/1045X-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/1047X-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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • A61N5/1036Leaf sequencing algorithms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1071Monitoring, verifying, controlling systems and methods for verifying the dose delivered by the treatment plan

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

A kind of device and equipment of dosage distribution determination and radiotherapy treatment planning optimization

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.
CN201511031048.XA 2015-12-31 2015-12-31 A kind of device and equipment of dosage distribution determination and radiotherapy treatment planning optimization CN105664376B (en)

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 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 (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107041999A (en) * 2017-02-16 2017-08-15 西南医科大学附属医院 There is 3-dimensional dose contribution distributed model construction method and its application of barrier material

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103038669A (en) * 2010-06-22 2013-04-10 卡尔·奥托 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

Also Published As

Publication number Publication date
CN105664376A (en) 2016-06-15

Similar Documents

Publication Publication Date Title
Unkelbach et al. Reducing the sensitivity of IMPT treatment plans to setup errors and range uncertainties via probabilistic treatment planning
US9630025B2 (en) Methods and apparatus for the planning and delivery of radiation treatments
Webb The physical basis of IMRT and inverse planning
EP1909904B1 (en) Methods and apparatus for the planning of radiation treatments
Jia et al. GPU-based fast Monte Carlo dose calculation for proton therapy
US7734010B2 (en) Method and apparatus for planning and delivering radiation treatment
Holmes et al. A comparison of three inverse treatment planning algorithms
Tillikainen et al. A 3D pencil-beam-based superposition algorithm for photon dose calculation in heterogeneous media
Wu et al. Re-optimization in adaptive radiotherapy
Miften et al. Implementation of FFT convolution and multigrid superposition models in the FOCUS RTP system
Shepard et al. An arc‐sequencing algorithm for intensity modulated arc therapy
Baum et al. Robust treatment planning for intensity modulated radiotherapy of prostate cancer based on coverage probabilities
Holmes et al. A filtered backprojection dose calculation method for inverse treatment planning
Van Santvoort et al. Dynamic multileaf collimation withouttongue-and-groove'underdosage effects
de la Zerda et al. Formulating adaptive radiation therapy (ART) treatment planning into a closed-loop control framework
Cabrera et al. An evolutionary algorithm for path synthesis of mechanisms
Jia et al. GPU-based fast Monte Carlo simulation for radiotherapy dose calculation
Oliver et al. Clinical significance of multi-leaf collimator positional errors for volumetric modulated arc therapy
Alber et al. Optimization of intensity modulated radiotherapy under constraints for static and dynamic MLC delivery
Cameron Sweeping-window arc therapy: an implementation of rotational IMRT with automatic beam-weight calculation
Altay et al. SPHRAY: a smoothed particle hydrodynamics ray tracer for radiative transfer
US7369645B2 (en) Information theoretic inverse planning technique for radiation treatment
Kamath et al. Leaf sequencing algorithms for segmented multileaf collimation
Bush et al. Monte Carlo simulation of RapidArc radiotherapy delivery
Tacke et al. Real‐time tumor tracking: Automatic compensation of target motion using the Siemens 160 MLC

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