CN110246562A - Determine the method, apparatus and computer system of the sub-beam intensity in radiotherapy system - Google Patents
Determine the method, apparatus and computer system of the sub-beam intensity in radiotherapy system Download PDFInfo
- Publication number
- CN110246562A CN110246562A CN201810185373.9A CN201810185373A CN110246562A CN 110246562 A CN110246562 A CN 110246562A CN 201810185373 A CN201810185373 A CN 201810185373A CN 110246562 A CN110246562 A CN 110246562A
- Authority
- CN
- China
- Prior art keywords
- launched field
- sub
- intensity
- launched
- spare
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
Landscapes
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Surgery (AREA)
- Urology & Nephrology (AREA)
- Engineering & Computer Science (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
The application provides the method, apparatus and computer system for determining the sub-beam intensity in radiotherapy system.According to the present processes based on compressive sensing theory by the way that the position of launched field and sub-beam to be introduced into traditional objective function to obtain objective function.Objective function is optimized to select spare launched field from multiple launched fields again.Based on spare launched field, sub-beam intensity is obtained by optimizing traditional intensity modulated radiation therapy plan.
Description
Technical field
The application relates generally to the method and device and computer system for determining the sub-beam intensity in radiotherapy system.
Background technique
Technical way one of of the radiotherapy as treatment malignant tumour, the purpose is to the maximum extent by radioactive ray
Cover makes the normal tissue of surrounding receive the smallest dosage as far as possible to tumour (needing radiotherapy tissue or target area).By
The IMRT (intensity modulated radiation therapy, intensity modulated radiation therapy) that three dimensional conformal radiation therapy develops is
One of radiotherapy.IMRT requires sub-beam (also referred to as tuftlet) intensity in launched field to be adjusted according to certain requirement
Section.It is throughout launched field and target area shape it is consistent under the conditions of, for the tool of target area 3D shape and vital organs and target area
Sub- intensity of beam is adjusted in body anatomy relationship, and dose contribution is in non-uniform but entire Gross Target Volume in single launched field
Dose contribution is treated more evenly than three-dimensional conformal.
Each launched field is divided into multiple sub-beams.When making a plan, according to target area 3D shape and to it is related
Jeopardize the anatomy relationship between organ, these sub-beams are distributed with different weights, makes to generate optimization in the same launched field
, non-uniform intensity contribution, to reduce the beam flux by jeopardizing organ, and the beam flux of target area other parts
Increase.
Intensity modulated radiation therapy plan is clinically formulated using inverse algorithm at present.It needs according to user to tumour uniform irradiation
The requirement of dosage and the normal tissue Organoprotective dose requirements of surrounding, need user rule of thumb specify launched field position and
Number and jeopardizes the dose constraint on organ at the target dose in tumour, obtained by optimization algorithm iteration include in each launched field
Sub-beam number, shape and weight (intensity corresponding to sub-beam).However, due to the position of launched field and number be user according to
Micro-judgment, the intensity of the sub-beam in each launched field can not be made to be optimal.Therefore, final dose distribution is being ensured compliance with
While reduce sub-beam number it is actually particularly important.
Summary of the invention
In view of the above-mentioned problems, the application offer can select spare launched field to more accurately determine radiotherapy from multiple launched fields
The method, apparatus and computer system of sub-beam intensity in system.
According to the one side of the application, a kind of method of the sub-beam intensity in determining radiotherapy system, radiotherapy system includes
Multiple launched fields, each launched field have multiple sub-beams, and method includes: that 1) will need radiotherapy tissue and needing around radiotherapy tissue
Surrounding tissue is divided into multiple voxels and determines radiological dose in each of multiple voxels;2) it is penetrated according to each launched field and son
Beam is relative to the position for needing radiotherapy tissue and surrounding tissue and to need radiotherapy tissue and the organizational information of surrounding tissue to simulate more
Each sub-beam in a launched field contributes the first unit dose of each voxel under unit strength;3) according to radiological dose and
First object function is established in the contribution of first unit dose;4) spare penetrate is selected from multiple launched fields by optimizing first object function
It is wild;5) determine each sub-beam in spare launched field under unit strength to each according to selected spare launched field and organizational information
Second unit dose of a voxel is contributed;6) the second objective function is constructed according to the contribution of the second unit dose and radiological dose;7)
The intensity of the sub-beam in spare launched field is determined by optimizing the second objective function.
It can be gone out respectively by Monte-carlo Simulation in step 2) and step 5) according to presently filed embodiment
The contribution of first unit dose and the contribution of the second unit dose.
According to presently filed embodiment, first object function can be constructed are as follows:
Wherein, p meets the restrictive condition of 0≤p≤1,Each of
Element corresponds to radiological dose, A1Indicate that the first unit dose is contributed,In each element correspond to multiple launched fields in one
Sub-beam intensity, xijIndicate the intensity of j-th of sub-beam of i-th of launched field of multiple launched fields,It indicates in multiple launched fields
I-th of launched field in maximum sub-beam intensity, β is weight coefficient.
According to presently filed embodiment, in the case where p=1, can be optimized by Newton method or gradient descent method
One objective function.
It,, can be by the way that first object function be turned in the case where p=1 in step 4) according to presently filed embodiment
It turns to:To be optimized to first object function, wherein yiMeet yi≥xijRestrictive condition,
T isThe quadratic form of item.
In transformed first object function, p meets the restrictive condition of 0≤p≤1,In each element correspond to put
Penetrate dosage, A1Indicate that the first unit dose is contributed,In each element correspond to multiple launched fields in a sub- intensity of beam,
xijIndicate the intensity of j-th of sub-beam of i-th of launched field of multiple launched fields, yiIt indicates in i-th of launched field in multiple launched fields most
Big sub-beam intensity, β are weight coefficients.
According to presently filed embodiment, in step 4), the second objective function can be constructed are as follows:
Wherein, A2Indicate that the second unit dose is contributed,In each element correspond to spare launched field in a sub-beam
Intensity,In each element correspond to radiological dose.
According to presently filed embodiment, in step 4), first object function can further be constructed are as follows:
Wherein, p meets the restrictive condition of 0≤p≤1,In each element be multiple voxels in need combination radiotherapy group
Knit the radiological dose of a corresponding voxel, APTVIndicate each sub-beam under unit strength to corresponding with radiotherapy tissue is needed
The unit dose of each voxel is contributed, and N indicates the quantity of the surrounding tissue, AOARkIndicate each sub-beam under unit strength
Unit dose contribution to each voxel corresponding with k-th of surrounding tissue,In each element correspond to multiple launched fields in
One sub- intensity of beam, xijIndicate the intensity of j-th of sub-beam of i-th of launched field in multiple launched fields,Indicate multiple
Maximum sub-beam intensity, β and λ in i-th of launched field in launched fieldkFor weight coefficient.
According to presently filed embodiment, in step 4), the first object letter can be optimized by global optimization approach
Number.
According to the another aspect of the application, a kind of device of the sub-beam intensity in determining radiotherapy system, radiotherapy system are provided
System include multiple launched fields, each launched field have multiple sub-beams, device include dosage determination unit, spare launched field selecting unit and
Strength determining unit.Dosage determination unit can be configured to that radiotherapy tissue will be needed and the surrounding tissue around radiotherapy tissue needed to draw
It is divided into multiple voxels and determines radiological dose in each of multiple voxels.Spare launched field selecting unit can be configured to: according to
Each launched field and sub-beam are relative to the position for needing radiotherapy tissue and surrounding tissue and need the group of radiotherapy tissue and surrounding tissue
Each sub-beam that information simulation goes out in multiple launched fields is knitted to contribute the first unit dose of each voxel under unit strength;Root
First object function is established according to radiological dose and the contribution of the first unit dose;And it is penetrated by optimization first object function from multiple
Select spare launched field in Yezhong.Strength determining unit can be configured to: be determined according to selected spare launched field and organizational information spare
Each sub-beam in launched field contributes the second unit dose of each voxel under unit strength;According to the second unit dose tribute
It offers and constructs the second objective function with radiological dose;And the sub-beam in spare launched field is determined by optimizing the second objective function
Intensity.
According to the another aspect of the application, a kind of computer system of the sub-beam intensity in determining radiotherapy system is provided,
Radiotherapy system includes multiple launched fields, and each launched field has multiple sub-beams, the system for determining the sub-beam intensity in radiotherapy system
Include: memory, is stored with computer instruction;Processor executes the computer instruction of memory storage to execute following operation:
It will need radiotherapy tissue and the surrounding tissue around radiotherapy tissue needed to be divided into multiple voxels and determine in multiple voxels
Each radiological dose;Relative to the position for needing radiotherapy tissue and surrounding tissue and radiotherapy is needed according to each launched field and sub-beam
Tissue and the organizational information of surrounding tissue simulate each sub-beam in multiple launched fields under unit strength to each voxel
The contribution of first unit dose;First object function is established according to radiological dose and the contribution of the first unit dose;Pass through optimization first
Objective function selects spare launched field from multiple launched fields;It is determined in spare launched field according to selected spare launched field and organizational information
Each sub-beam the second unit dose of each voxel is contributed under unit strength;It contributes and puts according to the second unit dose
Penetrate the second objective function of dose construction;The intensity of the sub-beam in spare launched field is determined by optimizing the second objective function.
Detailed description of the invention
The illustrative embodiments made in conjunction with the accompanying drawings are described below, these and or other aspects of the application
It will be apparent and it is more readily appreciated that in the accompanying drawings:
Fig. 1 is the process that the method for the sub-beam intensity in radiotherapy system is determined according to the application illustrative embodiments
Figure;
Fig. 2 is the schematic plan view according to the launched field including multiple sub-beams of the illustrative embodiments of the application;
Fig. 3 is to need showing for radiotherapy tissue and surrounding tissue by the irradiation of multiple launched fields according to the application illustrative embodiments
Meaning property perspective view;
Fig. 4 is to be selected in intensity modulated radiation therapy method by the method for the sub-beam intensity in the determination radiotherapy system of the application
Schematic diagram of the spare launched field to the dose contribution for needing radiotherapy tissue and surrounding tissue;
Fig. 5 is the spare launched field rule of thumb selected in intensity modulated radiation therapy method by user to needing radiotherapy tissue and surrounding group
The schematic diagram for the dose contribution knitted;
Fig. 6 is the block diagram for showing the device that the sub-beam intensity in radiotherapy system is determined according to the embodiment of the present application;And
Fig. 7 is to show the structural schematic diagram for the computer system for being suitable for carrying out the embodiment of the present application.
Specific embodiment
Since the disclosure allows a variety of changes and many embodiments, so will be illustrated in the drawings and write
Embodiment is described in detail in specification.When being referred to the embodiment described with reference to the accompanying drawings, the aspect of the disclosure and
Feature and realize that these method or a variety of methods will be apparent.However, such disclosure can be with many different shapes
Formula is realized, and should not be interpreted as limited to embodiment set forth herein.
It is as used herein like that, term "and/or" include any of one or more of related listed item and
All combinations.In addition, the use of "available" is related to " one or more implementations of the invention when describing embodiments of the present invention
Mode ".Such as " at least one of " when being expressed in after a column element, modification permutation element without modify column in
Individual component.In addition, term " exemplary " is intended to indicate that example or illustration.Cover the statement of plural number with the statement that odd number uses,
Unless it has significantly different meaning within a context.
It is as used herein like that, unless context clearly otherwise indicates, otherwise singular " one " and "one"
It is intended to also include plural form.It will also be understood that terms used herein " packet ", " comprising ", " comprising " and/or " comprising " point out institute
The feature of elaboration or the presence of component, but it is not excluded for the presence or addition of one or more other features or component.
As used herein " substantially ", " about ", " approximation " and similar wording is used as approximate wording and is not used as
The wording of degree, and Ei is intended to illustrate the inherence on the measured value or calculated value recognized by those of ordinary skill in the art partially
Difference.In addition, as used herein these wording include as by those of ordinary skill in the art in view of the measurement discussed and with spy
The relevant error of quantitative measurement (that is, limitation of measuring system) and determine particular value deviation tolerance interval in
Setting value and average value.For example, " about " can indicate in one or more standard deviations, or the value described ±
30%, in ± 20%, ± 10%, ± 5%.
When some embodiment can differently to realize when, particular procedure sequence can with described sequence differently by
It executes.For example, two processes continuously described can be executed with the roughly the same time, or with opposite with described sequence
Sequence executes two processes continuously described.
It is carried out below in conjunction with method and device of the Fig. 1 to Fig. 6 to the small beam intensity in the determination radiotherapy system of the application
Detailed description.
In the method for the sub-beam intensity in the determination radiotherapy system according to the application, for example, radiotherapy system may include
Launched field F shown in multiple Fig. 3 (also referred to as wild), and each launched field has multiple sub-beam B.It is shown in Fig. 2 to penetrate source F packet
Include 9 sub-beam B.For example, each launched field F by MLC (Multi-Leaf Collimato, multi-diaphragm collimator) be embodied as it is multiple
Sub-beam B.Launched field F shown in Fig. 2 can be by including penetrating source and MCL realization in radiotherapy system.The source of penetrating can emit X-ray
Or gamma-rays, i.e. launched field F are X-ray fields.As another example, positively charged ion (i.e. proton) can be emitted by penetrating source.But this
Be merely exemplary, according to the application to penetrate source without being limited thereto.In one example, the number of the sub-beam B in each launched field F
Can be same or different, and each sub-beam B can have different shapes.
Referring to Figure 1 and Figure 3, in step sl, radiotherapy tissue PTV will be needed and needing the surrounding tissue around radiotherapy tissue
OAR is divided into multiple voxels and determines radiological dose in each of multiple voxels.In the example depicted in fig. 3, it need to put
It treats tissue PTV and is divided into multiple voxels, and surrounding tissue OAR is divided into around the multiple voxels for needing radiotherapy tissue PTV.
That is, the three-D volumes that the voxel for needing radiotherapy tissue PTV and surrounding tissue OAR to be divided can be equal to each other.Voxel is numerical data
Minimum unit in three-dimensional space segmentation.20 launched field F are shown in Fig. 3, but this is merely exemplary.In multiple voxels
Each radiological dose indicates the radiological dose that each voxel needs to receive.For example, the radiation for needing the need of radiotherapy tissue PTV to receive
Dosage is greater than a certain predetermined value, and the radiological dose that the need of surrounding tissue OAR receive is less than a certain predetermined value.It is every in multiple voxels
A radiological dose can be determined according to the CT image of radiotherapy tissue PTV and surrounding tissue OAR is wanted.
In step s 2, according to each launched field F and sub-beam B relative to the position for needing radiotherapy tissue PTV and surrounding tissue OAR
Set and need the organizational information of radiotherapy tissue PTV and surrounding tissue OAR to simulate each sub-beam B in multiple launched field F to each
First unit dose of a voxel is contributed.
Need the organizational information of radiotherapy tissue PTV and surrounding tissue OAR can be by needing radiotherapy tissue PTV and surrounding tissue OAR
Computed tomography images (that is, CT image) obtain.For example, organizational information can be by based on reverse intensity modulated radiation therapy
CT (planing computed tomography, computed tomography) device is drawn to obtain.Organizational information may include needing radiotherapy
Organize volume, shape, position and the electron density information of the tissue of PTV and surrounding tissue OAR.Multiple launched field F and sub-beam B phase
For needing the position of radiotherapy tissue PTV and surrounding tissue OAR usually can artificially be configured.For example, the position of settable launched field F
It sets so that the ray in launched field F passes through and is radiated at from the gap between voxel corresponding with surrounding tissue OAR as much as possible
It needs on radiotherapy tissue PTV.But anyway be arranged launched field F position, surrounding tissue OAR receive dosage be impossible to for
Zero.
First unit dose contribution be each sub-beam B under unit strength to the dose contribution of each voxel.
In step s3, first object function is established according to identified radiological dose and the contribution of the first unit dose.
In step s 4, spare launched field P is selected from multiple launched field F by optimizing first object function.
In reverse intensity modulated radiation therapy, in the case where radiological dose and the contribution of the first unit dose determine, main purpose is
The intensity of all sub-beam B in multiple launched field F is determined, so that the dosage that each voxel receives meets it as far as possible and should receive
Radiological dose.
Therefore, objective function used in inverse intensity-modulated are as follows:
Wherein,In each element correspond to radiological dose determining in step sl, A expression determines in step s 2
The first unit dose contribution, vectorIn each element correspond to step S1 and step S2 in multiple launched field F a son
The intensity of beam B, vectorThe intensity of all sub-beam B be can be by the concatenated vector of row.The solution of traditional objective functionIt is to make
Obtain vectorSecond order norm obtain its minimum value all sub-beam B intensity by the concatenated vector of row.
Therefore, the position of launched field F and sub-beam B are introduced by conventional target letter according to the first object function of the application
In number.
It, can be by be introduced into traditional objective function in multiple launched field F i-th according to the application embodiment
The intensity of j-th of sub-beam B of launched field, to introduce the position of launched field F and sub-beam B.In other words, based on i-th launched field
The intensity of j-th of sub-beam B and the mapping relations of corresponding launched field F and sub-beam B introduce the position of launched field F and sub-beam B indirectly
It sets, and then further determines that F.
Further, ready for use penetrate is determined from multiple launched field F by experience different from doctor in traditional intensity modulated radiation therapy
Compressed sensing is applied to reverse intensity modulated radiation therapy to select from multiple launched field F by wild (hereinafter referred to spare launched field), the application
More accurate spare launched field P, so that it is determined that spare launched field P and including corresponding sub-beam B relative to needing combination radiotherapy group
Knit the position of PTV and surrounding tissue OAR.
Basic goal due to establishing first object function is that several spare launched field P are selected from multiple launched field F, that is,
While meeting the radiological dose to be implemented, the number of launched field F is reduced.
Therefore objective function (i.e. formula 1) base in the present processes, used in inverse intensity-modulated
On plinth, it is based on compressive sensing theory, passes through the strong of the sub-beam B that is arranged so that sub-beam intensity set in each launched field F
The sum of degree maximum value and the product of weight coefficient obtain its minimum value (minimum value can be preset value) and are used as restrictive condition, will be every
The maximum of intensity of sub-beam B in a launched field F seeks the sparse of sampled signal as sampled signal, and then by compressed sensing
Spare launched field P is chosen in Xie Laicong launched field F, wherein enable to need radiotherapy when sub-beam strength set contract in spare launched field P
The radiological dose that tissue PTV and surrounding tissue OAR finally receives meets or closest to intensity modulated radiation therapy target (that is, pre-set
Radiological dose).Therefore, compression sensing method is also considered as at signal according to the method for the application embodiment
Application in reason.
On the one hand, compressed sensing is the technology for rebuilding sparse or compressible signal.The mathematics of sparse signal shows
Comprising many close to zero or null coefficient when it is expressed in some domains.That is, if signal only limited non-zero is adopted
Sampling point, and other sampled points be zero (or close to zero), then it is sparse for claiming signal.Efficiently use compressed sensing technology (i.e.
The sparsity of signal to be processed) it can be used for reducing required measurement number of samples.Mathematically, intensity modulated radiation therapy Inverse Planning is similar
In being signal processing problems wait ask in the case where " signal " in above-mentioned multiple launched field F (i.e. all existing launched field F).As above
Described, Inverse Planning is an ill-conditioning problem, and usually there is the set of the spare launched field P of multiple groups (that is, from the multiple launched field F
The set of the launched field of middle selection, such as rule of thumb selected by doctor) radiation for needing radiotherapy tissue PTV to need to receive can be generated
Dosage, but in this case, surrounding tissue OAR may not can be made to receive least radiological dose.Therefore, in this application, lead to
It crosses using the maximum of intensity of the sub-beam B in each launched field F as sampled signal, the sparse solution for seeking sampled signal comes from multiple
Launched field F chooses spare launched field P, while enabling to the sub-beam intensity set in spare launched field P to need radiotherapy tissue PTV and week
The radiological dose that tissue OAR finally receives is enclosed to meet or (such as meet formula as restrictive condition closest to intensity modulated radiation therapy target
(1)).The number of nonzero value corresponds to the number of spare launched field P in striked sparse solution.
On the other hand, using the maximum of intensity of the sub-beam B in each launched field F as in the case where sampled signal, lead to
Each single item in the obtained sparse solution of overcompression cognitive method can correspond to the maximum sub-beam intensity in a launched field F, so
It is that zero item indicates that the intensity of all sub-beam B in corresponding launched field F is zero in sparse solution, i.e., it will be corresponding to this
Launched field F be selected as spare launched field P.In this way, it can be achieved that the position of each single item and launched field F in sparse solution maps pass correspondingly
System.
In addition, the number of settable alternative launched field P, such as weight coefficient are bigger, alternative launched field by setting weight coefficient
The number of P is fewer.
It, can be by first object function component according to the application embodiment are as follows:
Wherein, p meets the restrictive condition of 0≤p≤1,In each element correspond to voxel in multiple voxels
Radiological dose, A indicate the first unit dose contribution described above,In each element correspond to sub-beam intensity,It can be with
It is the intensity of all sub-beam B by the concatenated vector of row.xijIndicate j-th of sub-beam of i-th of launched field in multiple launched field F
Intensity,Indicate that maximum sub-beam intensity, β are weight coefficients in i-th of launched field in multiple launched field F.
The number of the alternative launched field P selected from multiple launched field F can be set by the way that the value of weight coefficient β is arranged.β value is got over
Greatly, the number of alternative launched field P is fewer.The value of β is greater than zero.
According to presently filed embodiment, first object function (for example, formula (2)) can be optimized by global optimization approach,
To obtain vectorVectorIn an element correspond to corresponding launched field F in corresponding sub-beam B intensity.That is,
The intensity of each sub-beam B corresponds to unique sub-beam B in unique launched field F.Therefore, vector can be passed throughIn each element
With the mapping relations of the sub-beam B in launched field F and launched field F, spare launched field P is selected from multiple launched field F.
According to the application embodiment, can by will be needed in first object function radiotherapy tissue PTV and around
The voxel of tissue OAR, which distinguishes and weight coefficient is arranged, needs radiotherapy tissue PTV to determine that the radiological dose finally realized is closer
Radiotherapy dosimetry or surrounding tissue OAR Radiotherapy dosimetry.
Therefore, according to the another exemplary embodiment of the application, first object function can further be constructed are as follows:
Wherein, p meets the restrictive condition of 0≤p≤1,In each element be multiple voxels in need combination radiotherapy group
Knit the radiological dose of the corresponding voxel of PTV, APTVIndicate each sub-beam B under unit strength to need radiotherapy tissue PTV
The unit dose of corresponding each voxel is contributed, and N indicates the quantity of the surrounding tissue, AOARkIndicate each sub-beam B in list
The unit dose of each voxel corresponding with k-th of surrounding tissue OAR is contributed under the intensity of position,In each element correspond to
The intensity of sub-beam B,The intensity of all sub-beam B be can be by the concatenated vector of row, xijIndicate i-th in multiple launched field F
The intensity of j-th of sub-beam of a launched field,Indicate that maximum sub-beam is strong in i-th of launched field in multiple launched field F
Degree, β is weight coefficient, λkFor weight coefficient
For example, the quantity of N is 6 in Fig. 3 into example shown in fig. 5, i.e., the quantity of surrounding tissue OAR is 6.It can pass through
λ is setkNumerical value put to determine that the spare launched field P finally selected meets radiological dose more needed for needing radiotherapy tissue PTV
Penetrate Radiotherapy dosimetry needed for dosage still more meets surrounding tissue OAR.
According to presently filed embodiment, first object function (for example, formula (3)) can be optimized by global optimization approach,
To obtain vectorAnd then multiple spare launched field P are selected from multiple launched field F.
In the case that p in first object function (for example, formula (2) or formula (3)) is equal to 1, can by such as Newton method,
The conventional iterative algorithm of gradient descent method, conjugate gradient method or quasi-Newton method optimizes first object function (for example, formula (2)
Or formula (3)).
That is, in the case where 0≤p≤1, can by global optimization approach to first object function (for example, formula
(2) or formula (3)) optimize, to obtain the globally optimal solution of objective function.Global optimization research is multivariable nonlinearity letter
The characteristic and construction of globally optimal solution of the number on some constraint seek the calculation method of globally optimal solution, and solution side
The theory property of method and calculating performance.For example, such as simulated annealing, genetic algorithm, TABU search, population can be used to calculate
Method, the global optimization approach of ant group algorithm optimize first object function.But the application is without being limited thereto.
First object function (for example, in the case that p in formula (2) or formula (3) is equal to 1, can by such as Newton method,
The iterative algorithm of gradient descent method, conjugate gradient method or quasi-Newton method optimizes first object function (for example, formula (2) or formula
(3))。
In the method according to the embodiment of the application, in first object function (for example, formula (2) or formula (3))
In p be equal to 1 in the case where, will be in 1 rank norm in formula (6) or formula (7)Item is used as constraint condition, will
Nonlinear problem is converted into linear problem, to be further reduced calculation amount.
By taking formula (2) as an example, that is, convert first object function are as follows:
Wherein, yiMeet yi≥xijRestrictive condition, T isThe quadratic form of item.
In formula (4),In each element correspond to the radiological dose of a voxel in multiple voxels, A indicates above
The first unit dose contribution,In each element correspond to sub-beam intensity,It can be the strong of all sub-beam B
Degree is by the concatenated vector of row, xijIndicate the intensity of j-th of sub-beam of i-th of launched field in multiple launched field F,It indicates
Maximum sub-beam intensity, β are weight coefficients in i-th of launched field in multiple launched field F.
In this way, the iterative algorithms such as Newton method, gradient descent method can be passed through in the case where reducing calculation amount to target
Function optimizes, to obtain vectorFurther, vector can be passed throughIn each element and launched field F and launched field F in
Sub-beam B mapping relations, spare launched field P is selected from multiple launched field F.
Multiple launched field F (examples can be selected from 20 launched field F shown in Fig. 3 by step S3 and step S4 referring to Fig. 3
Such as, 3 launched field F) it is used as spare launched field P.
In step s 5, each sub-beam in spare launched field P is determined according to selected spare launched field P and organizational information
B contributes the second unit dose of each voxel under unit strength.For example, according to each launched field F and sub-beam B relative to need
The position of radiotherapy tissue PTV and surrounding tissue OAR and radiotherapy tissue PTV and the organizational information of surrounding tissue OAR is needed to simulate
The contribution of second unit dose.For example, being contributed by the second unit dose of Monte-carlo Simulation.Due to establishing first object
Function (step S3) consideration enables to the sub-beam intensity set in spare launched field P to need radiotherapy tissue PTV and surrounding tissue
The radiological dose that OAR finally receives meets or closest to intensity modulated radiation therapy target as restrictive condition, therefore, passes through spare launched field P
The the second unit dose contribution simulated enables to the result of Inverse Planning closer to ideal, i.e., the sub-beam finally acquired is strong
Degree can make surrounding tissue OAR receive least radiological dose while making to need radiotherapy tissue PTV acceptance presupposition radiological dose.
In step s 6, the second objective function is constructed according to the contribution of the second unit dose and radiological dose.
In presently filed embodiment, the second objective function can correspond to the conventional target letter with intensity modulated radiation therapy in the works
Number, it may be assumed that
Wherein, A2Indicate the second unit dose contribution,In each element correspond to alternative launched field P in one
Sub-beam intensity,All sub-beam intensity in alternative launched field P be can be by the concatenated vector of row,In each element pair
Radiological dose described in Ying Yu.
In the step s 7, the intensity of the sub-beam B in spare launched field P is determined by optimizing the second objective function.For
The case where two objective functions are formula (5), can optimize formula (5) by Newton method or gradient descent method.As the application's
Global optimization approach can be used to optimize formula (5) for another example.But the application is without being limited thereto.
The spare launched field pair that the method that Fig. 4 shows the sub-beam intensity in the determination radiotherapy system by the application is selected
Need the schematic diagram of the dose contribution of radiotherapy tissue PTV and surrounding tissue OAR.Fig. 5 be in intensity modulated radiation therapy method by user according to
Schematic diagram of the spare launched field P that experience is selected to the dose contribution for needing radiotherapy tissue PTV and surrounding tissue OAR.
In figures 4 and 5, the radiological agent of voxel receiving is schematically shown with the quantity of lines in each tissue
Amount, that is to say, that the radiological dose that the quantity of lines more multilist shows that how corresponding tissue receives is more.Pass through step S3 and step
S4 selects 3 as spare launched field P and again by step S5 to step S7 based on spare from 20 launched field F shown in Fig. 3
The intensity of launched field P acquisition sub-beam B.Surrounding tissue compared to radiological dose suffered by surrounding tissue OAR in Fig. 5, in Fig. 4
Radiological dose suffered by OAR significantly reduces.
In formula (2) into (9), of the alternative launched field P selected from multiple launched field F can be set by the way that the value of β is arranged
Number.β value is bigger, and the number of alternative launched field P is fewer.The value of β can be greater than zero.
It can be implemented by computer system referring to figs. 1 to the method for Fig. 5 description.The computer system can including being stored with
The memory and processor executed instruction.Processor is communicated with memory to execute executable instruction to implement referring to Fig.1
The method described to Fig. 5.Alternatively, or in addition, it can be deposited by non-transitory computer referring to figs. 1 to the method for Fig. 5 description
Storage media is implemented.The medium storing computer readable instruction, when these instruction be performed make processor execute referring to figs. 1 to
The method of Fig. 5 description.
Fig. 6 is the block diagram for showing the device 100 that the sub-beam intensity in radiotherapy system is determined according to the embodiment of the present application.Really
The device 100 for determining the sub-beam intensity in radiotherapy system includes: dosage determination unit 110, be configured to will to need radiotherapy tissue PTV with
Each of and be divided into multiple voxels in the surrounding tissue OAR needed around radiotherapy tissue and determine the multiple voxel
Radiological dose;Spare launched field selecting unit 120, being configured that needs combination radiotherapy group for described according to each launched field F and sub-beam B
It knits and the position of the surrounding tissue and the organizational information of radiotherapy tissue PTV and surrounding tissue OAR is needed to simulate multiple launched fields
In each sub-beam the first unit dose of each voxel is contributed under unit strength;According to radiological dose and first
First object function is established in unit dose contribution;And spare penetrate is selected from multiple launched field F by optimizing first object function
Wild P;Strength determining unit is configured that and is determined in spare launched field P according to selected spare launched field P and the organizational information
Each sub-beam B contributes the second unit dose of each voxel under unit strength;It contributes and puts according to the second unit dose
Penetrate the second objective function of dose construction;The intensity of the sub-beam B in spare launched field P is determined by second objective function.
According to the application embodiment, spare launched field selecting unit 120 and strength determining unit 130 can lead to respectively
It crosses Monte-carlo Simulation and goes out the contribution of the first unit dose and the contribution of the second unit dose.
According to the application embodiment, spare launched field selecting unit 120 by first object function is constructed are as follows:
In formula (6), p meets the restrictive condition of 0≤p≤1,In each element correspond to multiple voxels in one
The radiological dose of voxel, A indicate the first unit dose contribution described above,In each element correspond to sub-beam intensity,The intensity of all sub-beam B be can be by the concatenated vector of row.xijIndicate j-th of son of i-th of launched field in multiple launched field F
The intensity of beam,Indicate that maximum sub-beam intensity, β are weight coefficients in i-th of launched field in multiple launched field F.
According to presently filed embodiment, first object function (for example, formula (6)) can be optimized by global optimization approach,
To obtain vectorVectorIn an element correspond to corresponding launched field F in corresponding sub-beam B intensity.That is,
The intensity of each sub-beam B corresponds to unique sub-beam B in unique launched field F.Therefore, vector can be passed throughIn each element
With the mapping relations of the sub-beam B in launched field F and launched field F, spare launched field P is selected from multiple launched field F.
According to the application embodiment, spare launched field selecting unit 120 can also by the first object function into
The building of one step are as follows:
Wherein, p meets the restrictive condition of 0≤p≤1,In each element be multiple voxels in need radiotherapy tissue
The radiological dose of the corresponding voxel of PTV, APTVIndicate each sub-beam B under unit strength to need PTV pairs of radiotherapy tissue
The unit dose for each voxel answered is contributed, and N indicates the quantity of the surrounding tissue, AOARkIndicate each sub-beam B in unit
The unit dose of each voxel corresponding with k-th of surrounding tissue OAR is contributed under intensity,In each element correspond to son
The intensity of beam B,The intensity of all sub-beam B be can be by the concatenated vector of row, xijIndicate i-th in multiple launched field F
The intensity of j-th of sub-beam of launched field,Indicate maximum sub-beam intensity, β in i-th of launched field in multiple launched field F
It is weight coefficient, λkFor weight coefficient.
According to presently filed embodiment, first object function (for example, formula (7)) can be optimized by global optimization approach,
To obtain vectorAnd then multiple spare launched field P are selected from multiple launched field F.
That is, spare launched field selecting unit 120 can be by global optimization approach in the case where 0≤p≤1
One objective function (for example, formula (6) or formula (7)) optimizes, to obtain the globally optimal solution of objective function.Global optimization research
The characteristic for being globally optimal solution of the multivariable nonlinearity function on some constraint and construction seek globally optimal solution
The theory property of calculation method and method for solving and calculating performance.For example, such as simulated annealing, heredity can be used to calculate
Method, TABU search, particle swarm algorithm, the global optimization approach of ant group algorithm optimize first object function.But the application
It is without being limited thereto.
In the case that p in first object function (for example, formula (6) or formula (7)) is equal to 1, can by such as Newton method,
The iterative algorithm of gradient descent method, conjugate gradient method or quasi-Newton method optimizes first object function (for example, formula (6) or formula
(7)).According to the application embodiment, the p in first object function (for example, formula (6) or formula (7)) is equal to 1 feelings
Under condition, spare launched field selecting unit 120 can pass through such as Newton method, gradient descent method, conjugate gradient method or quasi-Newton method
Conventional iterative algorithm optimizes first object function (for example, formula (6) or formula (7)).
It will be in 1 rank norm in formula (6) or formula (7) in the case where p is equal to 1 according to the application embodiment
'sItem is used as constraint condition, linear problem is converted by nonlinear problem, to be further reduced calculation amount.
By taking formula (6) as an example, spare launched field selecting unit 120 is by converting first object function are as follows:
In transformed first object function (i.e. formula (8)), T isThe quadratic form of item, yiMeet yi≥xij's
Restrictive condition, that is,
In transformed first object function,In voxel corresponding in multiple voxels of each element putting
Dosage is penetrated, A indicates the first unit dose contribution described above,In each element correspond to sub-beam intensity,It can be
The intensity of all sub-beam B is by the concatenated vector of row, xijIndicate the strong of j-th of sub-beam of i-th of launched field in multiple launched field F
Degree,Indicate that maximum sub-beam intensity, β are weight coefficients in i-th of launched field in multiple launched field F.
According to the application embodiment, strength determining unit 130 can construct the second objective function are as follows:
Wherein, A2Indicate the second unit dose contribution,In each element correspond to alternative launched field P in one
Sub-beam intensity,All sub-beam intensity in alternative launched field P be can be by the concatenated vector of row,In each element pair
Radiological dose described in Ying Yu.
In the embodiment of the application, strength determining unit 130 can be determined standby by the second objective function of optimization
With the intensity of the sub-beam B in launched field P.The case where being formula (9) for the second objective function, can be declined by Newton method or gradient
Method optimizes formula (9).As another example of the application, global optimization approach can be used to optimize formula (9).But this
Apply without being limited thereto.
The device of the sub-beam intensity in determination radiotherapy system described referring to Fig. 6 can be implemented by computer system.
The computer system may include the memory and processor for being stored with executable instruction.Processor is communicated with memory to execute
Executable instruction is to implement the device referring to Fig. 6 description.Alternatively, or in addition, can pass through referring to Fig. 6 device described
Non-transitory computer storage medium is implemented.The medium storing computer readable instruction can when these instructions are by computer
Execute the function of the device referring to Fig. 6 description.
It is to show the structural representation for the computer system 7000 for being suitable for carrying out the embodiment of the present application referring now to Fig. 7, Fig. 7
Figure.
As shown in fig. 7, computer system 7000 may include processor (such as central processing unit (CPU) 7001, image processor
(GPU) etc.), it can be loaded at random according to the program being stored in read-only memory (ROM) 7002 or from storage section 7008
It accesses the program in memory (RAM) 7003 and executes various movements appropriate and process.In RAM 7003, can also it be stored with
Various programs required for system 7000 operates and data.CPU 7001, ROM 7002 and RAM 7003 by bus 7004 that
This connection.Input/output I/O interface 7005 is also connect with bus 7004.
The following are the components that can be connect with I/O interface 7005: the importation 7006 including keyboard, mouse etc.;Including yin
The output par, c 7007 of extreme ray pipe CRT, liquid crystal display LCD and loudspeaker etc.;Storage section 7008 including hard disk etc.;
And the communications portion 7009 including network interface card (such as LAN card and modem).Communications portion 7009 can be by all
Such as internet network executes communication process.As needed, driver 7010 can also be connect with I/O interface 7005.Such as disk, light
The detachable media 7011 of disk, magneto-optic disk, semiconductor memory etc. is mountable on driver 7010, in order to read from thereon
Computer program be mounted into storage section 7008 as needed.
Specifically, in accordance with an embodiment of the present disclosure, the method described above by reference to Fig. 1 to Fig. 5 it is implementable for computer it is soft
Part program.For example, embodiment of the disclosure may include computer program product, the product is machine readable including being tangibly embodied in
Computer program in medium.The computer program includes for executing the method program code referring to figs. 1 to Fig. 5 description.?
In this embodiment, computer program can be downloaded and be installed from network by communications portion 709, and/or can be from removable
Unload the installation of medium 7011.
Flow chart and block diagram in attached drawing are illustrated according to the system of various embodiments of the invention, method and computer journey
The architecture, function and operation in the cards of sequence product.In this regard, each box in flowchart or block diagram can generation
A part of one module, program segment or code of table, a part of the module, program segment or code include one or more
Executable instruction for implementing the specified logical function.It should also be noted that in some implementations as replacements, institute in box
The function of mark can also occur in a different order than that indicated in the drawings.For example, two boxes succeedingly indicated are practical
On can be basically executed in parallel, they can also be executed in the opposite order sometimes, and this depends on the function involved.Also it wants
It is noted that the combination of each box in block diagram and or flow chart and the box in block diagram and or flow chart, Ke Yiyong
The dedicated hardware based system of defined functions or operations is executed to realize, or can be referred to specialized hardware and computer
The combination of order is realized.
Unit or module involved in embodiments herein can be implemented by software or hardware.Described unit or mould
Block may also set up in the processor.The title of these units or module is not construed as limiting these units or module.
According to presently filed embodiment, the method, apparatus of the sub-beam intensity in provided determining radiotherapy system and
Computer system can be used in any kind of radiotherapy system.For example, realizing radiotherapy by X-ray, gamma-rays or proton
Radiotherapy system, but the application is without being limited thereto.
Above description is only the exemplary embodiment of the application and the explanation to institute's application technology principle.Art technology
Personnel should be appreciated that range involved in the application, however it is not limited to technical side made of the specific combination of above-mentioned technical characteristic
Case, while should also cover without departing substantially from the inventive concept, appointed by above-mentioned technical characteristic or its equivalent feature
Other technical solutions of meaning combination and formation.Such as features described above and the technical characteristic with similar functions disclosed herein
The technical solution replaced mutually and formed.
Claims (17)
1. a kind of method of the sub-beam intensity in determining radiotherapy system, the radiotherapy system includes multiple launched fields, each described
Launched field has multiple sub-beams, which comprises
1) radiotherapy tissue will be needed and be divided into multiple voxels in the surrounding tissue needed around radiotherapy tissue and determine institute
State radiological dose in each of multiple voxels;
2) according to each launched field and the sub-beam relative to it is described need the position of radiotherapy tissue and the surrounding tissue with
And each sub-beam for needing radiotherapy tissue and the organizational information of the surrounding tissue to simulate in the multiple launched field
The first unit dose of each voxel is contributed under unit strength;
3) first object function is established according to the radiological dose and first unit dose contribution;
4) spare launched field is selected from the multiple launched field by optimizing the first object function;
5) determine that each sub-beam in the spare launched field exists according to selected spare launched field and the organizational information
The second unit dose of each voxel is contributed under unit strength;
6) the second objective function is constructed according to second unit dose contribution and the radiological dose;
7) intensity of the sub-beam in the spare launched field is determined by optimizing second objective function.
2. the method for claim 1, wherein passing through Monte-carlo Simulation respectively in step 2) and step 5)
The first unit dose contribution and second unit dose contribution out.
3. method according to claim 1 or 2, wherein the first object function building are as follows:
Wherein, p meets the restrictive condition of 0≤p≤1,In each element correspond to the radiological dose, A1Indicate described
The contribution of one unit dose,In each element correspond to the multiple launched field in a sub- intensity of beam, xijDescribed in expression
The intensity of j-th of sub-beam of i-th of launched field of multiple launched fields,It indicates in i-th of launched field in the multiple launched field
Maximum sub-beam intensity, β are weight coefficients.
4. method as claimed in claim 3, wherein in the case where p=1, optimized by Newton method or gradient descent method
The first object function.
5. method as claimed in claim 4, wherein in step 4), by converting the first object function are as follows:The first object function optimized,
Wherein, yiMeet yi≥xijRestrictive condition, T isThe quadratic form of item, and
Wherein,In each element correspond to the radiological dose, A1Indicate the first unit dose contribution,Each of
Element corresponds to a sub- intensity of beam in the multiple launched field, xijIndicate j-th of i-th of launched field of the multiple launched field
The intensity of sub-beam, yiIndicate that maximum sub-beam intensity, β are weight coefficients in i-th of launched field in the multiple launched field.
6. method as claimed in claim 5, wherein construct second objective function are as follows:
Wherein, A2Indicate the second unit dose contribution,In each element correspond to the spare launched field in a son
Intensity of beam,In each element correspond to the radiological dose.
7. method as claimed in claim 3, wherein further construct the first object function are as follows:
Wherein, p meets the restrictive condition of 0≤p≤1,In each element be multiple voxels in it is corresponding with radiotherapy tissue is needed
A voxel radiological dose, APTVIndicate that each sub-beam needs radiotherapy tissue corresponding under unit strength to described
Each voxel unit dose contribution, N indicates the quantity of the surrounding tissue, AOARkIndicate each sub-beam in unit
The unit dose of each voxel corresponding with k-th of surrounding tissue is contributed under intensity,In each element correspond to institute
State a sub- intensity of beam in multiple launched fields, xijIndicate the strong of j-th of sub-beam of i-th of launched field in the multiple launched field
Degree,Indicate maximum sub-beam intensity, β and λ in i-th of launched field in the multiple launched fieldkFor weight coefficient.
8. method as claimed in claim 3, wherein in step 4), optimize the first object by global optimization approach
Function.
9. a kind of device of the sub-beam intensity in determining radiotherapy system, the radiotherapy system includes multiple launched fields, each described
Launched field has multiple sub-beams, and described device includes:
Dosage determination unit is configured to that radiotherapy tissue will be needed and is divided into the surrounding tissue needed around radiotherapy tissue more
Radiological dose in each of a voxel and determining the multiple voxel;
Spare launched field selecting unit, is configured that
According to each launched field and the sub-beam relative to the position for needing radiotherapy tissue and the surrounding tissue and
The each sub-beam for needing radiotherapy tissue and the organizational information of the surrounding tissue to simulate in the multiple launched field exists
The first unit dose of each voxel is contributed under unit strength;
First object function is established according to the radiological dose and first unit dose contribution;And
Spare launched field is selected from the multiple launched field by optimizing the first object function;And
Strength determining unit is configured that
Determine each sub-beam in the spare launched field in list according to selected spare launched field and the organizational information
The second unit dose of each voxel is contributed under the intensity of position;
The second objective function is constructed according to second unit dose contribution and the radiological dose;And
The intensity of the sub-beam in the spare launched field is determined by optimizing second objective function.
10. device as claimed in claim 9, wherein the spare launched field selecting unit and strength determining unit difference
Go out the first unit dose contribution and second unit dose contribution by Monte-carlo Simulation.
11. the device as described in claim 9 or 10, wherein the spare launched field selecting unit is by the first object function
Building are as follows:
Wherein, p meets the restrictive condition of 0≤p≤1,In each element correspond to the radiological dose, A1Indicate described
The contribution of one unit dose,In each element correspond to the multiple launched field in a sub- intensity of beam, xijDescribed in expression
The intensity of j-th of sub-beam of i-th of launched field of multiple launched fields,It indicates in i-th of launched field in the multiple launched field
Maximum sub-beam intensity, β are weight coefficients.
12. method as claimed in claim 11, wherein in the case where p=1, the spare launched field selecting unit passes through ox
Method or gradient descent method optimize the first object function.
13. device as claimed in claim 12, wherein the spare launched field selecting unit is by by the first object function
Conversion are as follows:The first object function optimized,
Wherein, wherein yiMeet yi≥xijRestrictive condition, T isThe quadratic form of item, and
Wherein,In each element correspond to the radiological dose, A1Indicate the first unit dose contribution,Each of
Element corresponds to a sub- intensity of beam in the multiple launched field, xijIndicate j-th of i-th of launched field of the multiple launched field
The intensity of sub-beam, yiIndicate that maximum sub-beam intensity, β are weight coefficients in i-th of launched field in the multiple launched field.
14. device as claimed in claim 9, wherein the strength determining unit constructs second objective function are as follows:
Wherein, A2Indicate the second unit dose contribution,In each element correspond to the spare launched field in a son
Intensity of beam,In each element correspond to the radiological dose.
15. method as claimed in claim 11, wherein the spare launched field selecting unit is by the first object function into one
Step building are as follows:
Wherein, p meets the restrictive condition of 0≤p≤1,In each element be multiple voxels in it is corresponding with radiotherapy tissue is needed
A voxel radiological dose, APTVIndicate that each sub-beam needs radiotherapy tissue corresponding under unit strength to described
Each voxel unit dose contribution, N indicates the quantity of the surrounding tissue, AOARkIndicate each sub-beam in unit
The unit dose of each voxel corresponding with k-th of surrounding tissue is contributed under intensity,In each element correspond to institute
State a sub- intensity of beam in multiple launched fields, xijIndicate the strong of j-th of sub-beam of i-th of launched field in the multiple launched field
Degree,Indicate maximum sub-beam intensity, β and λ in i-th of launched field in the multiple launched fieldkFor weight coefficient.
16. device as claimed in claim 11, wherein the spare launched field selecting unit optimizes institute by global optimization approach
State first object function.
17. a kind of computer system of the sub-beam intensity in determining radiotherapy system, the radiotherapy system includes multiple launched fields, often
A launched field has multiple sub-beams, and the system of the sub-beam intensity in the determining radiotherapy system includes:
Memory is stored with computer instruction;
Processor executes the computer instruction of the memory storage to execute as described in any one of claim 1 to 8
Method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810185373.9A CN110246562B (en) | 2018-03-07 | 2018-03-07 | Method, device and computer system for determining sub-beam intensity in radiotherapy system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810185373.9A CN110246562B (en) | 2018-03-07 | 2018-03-07 | Method, device and computer system for determining sub-beam intensity in radiotherapy system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110246562A true CN110246562A (en) | 2019-09-17 |
CN110246562B CN110246562B (en) | 2021-08-06 |
Family
ID=67876215
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810185373.9A Active CN110246562B (en) | 2018-03-07 | 2018-03-07 | Method, device and computer system for determining sub-beam intensity in radiotherapy system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110246562B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111986778A (en) * | 2020-07-31 | 2020-11-24 | 上海联影医疗科技股份有限公司 | Intensity modulated plan optimization system, device and storage medium |
CN115920257A (en) * | 2023-01-17 | 2023-04-07 | 安徽大学 | Beam angle optimization system for intensity modulated radiation therapy |
CN115966281A (en) * | 2023-01-16 | 2023-04-14 | 迈胜医疗设备有限公司 | Method, device, equipment and storage medium for generating arc radiotherapy plan |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1604134A (en) * | 2003-11-26 | 2005-04-06 | 戴建荣 | Method and system for optimizing radiotherapeutic radiation field orientation and strength distribution |
CN101120871A (en) * | 2006-12-29 | 2008-02-13 | 成都川大奇林科技有限责任公司 | Precise radiotherapy planning system |
CN101422640A (en) * | 2008-11-25 | 2009-05-06 | 中国科学院等离子体物理研究所 | Multiple-objective optimization method and system capable of optimizing radiotherapy beam intensity distribution |
US20110153547A1 (en) * | 2009-06-23 | 2011-06-23 | The Johns Hopkins University | System and method for shape based retrieval of prior patients for automation and quality control of radiation therapy treatment plans |
CN103083821A (en) * | 2013-01-27 | 2013-05-08 | 合肥超安医疗科技有限公司 | System and a method of directly optimizing static-state intensity-modulated sub-field shapes and machine hop count thereof |
CN103282967A (en) * | 2010-08-17 | 2013-09-04 | 德克萨斯州立大学董事会 | Automated treatment planning for radiation therapy |
CN104338240A (en) * | 2014-10-31 | 2015-02-11 | 章桦 | Automatic optimization method for on-line self-adaption radiotherapy plan and device |
CN104645500A (en) * | 2015-02-12 | 2015-05-27 | 上海联影医疗科技有限公司 | Method and system for optimizing intensity modulated radiation therapy |
CN105031819A (en) * | 2015-08-25 | 2015-11-11 | 上海联影医疗科技有限公司 | Dosage optimization method and system |
CN105031820A (en) * | 2015-09-01 | 2015-11-11 | 上海联影医疗科技有限公司 | Intensity modulated radiation therapy reverse optimization method and device |
CN105413068A (en) * | 2015-12-31 | 2016-03-23 | 上海联影医疗科技有限公司 | Optimization method for solving flux map and radiotherapy equipment |
CN105617536A (en) * | 2015-12-24 | 2016-06-01 | 上海联影医疗科技有限公司 | Inverse intensity modulated arc radioterapy optimization method and device |
US20170128029A1 (en) * | 2010-04-16 | 2017-05-11 | Scott Penfold | Simultaneous / single patient position x-ray and proton imaging apparatus and method of use thereof |
CN107073285A (en) * | 2014-06-17 | 2017-08-18 | 直觉疗法股份公司 | System and computer program product for reverse treatment plan |
-
2018
- 2018-03-07 CN CN201810185373.9A patent/CN110246562B/en active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1604134A (en) * | 2003-11-26 | 2005-04-06 | 戴建荣 | Method and system for optimizing radiotherapeutic radiation field orientation and strength distribution |
CN101120871A (en) * | 2006-12-29 | 2008-02-13 | 成都川大奇林科技有限责任公司 | Precise radiotherapy planning system |
CN101422640A (en) * | 2008-11-25 | 2009-05-06 | 中国科学院等离子体物理研究所 | Multiple-objective optimization method and system capable of optimizing radiotherapy beam intensity distribution |
US20110153547A1 (en) * | 2009-06-23 | 2011-06-23 | The Johns Hopkins University | System and method for shape based retrieval of prior patients for automation and quality control of radiation therapy treatment plans |
US20170128029A1 (en) * | 2010-04-16 | 2017-05-11 | Scott Penfold | Simultaneous / single patient position x-ray and proton imaging apparatus and method of use thereof |
CN103282967A (en) * | 2010-08-17 | 2013-09-04 | 德克萨斯州立大学董事会 | Automated treatment planning for radiation therapy |
CN103083821A (en) * | 2013-01-27 | 2013-05-08 | 合肥超安医疗科技有限公司 | System and a method of directly optimizing static-state intensity-modulated sub-field shapes and machine hop count thereof |
CN107073285A (en) * | 2014-06-17 | 2017-08-18 | 直觉疗法股份公司 | System and computer program product for reverse treatment plan |
CN104338240A (en) * | 2014-10-31 | 2015-02-11 | 章桦 | Automatic optimization method for on-line self-adaption radiotherapy plan and device |
CN104645500A (en) * | 2015-02-12 | 2015-05-27 | 上海联影医疗科技有限公司 | Method and system for optimizing intensity modulated radiation therapy |
CN105031819A (en) * | 2015-08-25 | 2015-11-11 | 上海联影医疗科技有限公司 | Dosage optimization method and system |
CN105031820A (en) * | 2015-09-01 | 2015-11-11 | 上海联影医疗科技有限公司 | Intensity modulated radiation therapy reverse optimization method and device |
CN105617536A (en) * | 2015-12-24 | 2016-06-01 | 上海联影医疗科技有限公司 | Inverse intensity modulated arc radioterapy optimization method and device |
CN105413068A (en) * | 2015-12-31 | 2016-03-23 | 上海联影医疗科技有限公司 | Optimization method for solving flux map and radiotherapy equipment |
Non-Patent Citations (3)
Title |
---|
GEORGIOS KALANTZIS AND ADITYA APTE: "A Novel Reduced-Order Prioritized Optimization", 《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 * |
朱力: "基于L_(1/2)范数X射线发光断层成像的稀疏算法重建研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
闵志方: "调强放疗中的数学规划问题研究", 《中国博士学位论文全文数据库·基础科学辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111986778A (en) * | 2020-07-31 | 2020-11-24 | 上海联影医疗科技股份有限公司 | Intensity modulated plan optimization system, device and storage medium |
CN111986778B (en) * | 2020-07-31 | 2024-02-20 | 上海联影医疗科技股份有限公司 | Intensity-modulated plan optimization system, device and storage medium |
CN115966281A (en) * | 2023-01-16 | 2023-04-14 | 迈胜医疗设备有限公司 | Method, device, equipment and storage medium for generating arc radiotherapy plan |
CN115966281B (en) * | 2023-01-16 | 2024-02-09 | 迈胜医疗设备有限公司 | Method, device, equipment and storage medium for generating arc radiotherapy plan |
CN115920257A (en) * | 2023-01-17 | 2023-04-07 | 安徽大学 | Beam angle optimization system for intensity modulated radiation therapy |
Also Published As
Publication number | Publication date |
---|---|
CN110246562B (en) | 2021-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kontaxis et al. | DeepDose: Towards a fast dose calculation engine for radiation therapy using deep learning | |
JP6437812B2 (en) | Method and radiation therapy system for optimizing radiation dose for radiation therapy treatment | |
Wieser et al. | Development of the open‐source dose calculation and optimization toolkit matRad | |
Gupta et al. | Generation of synthetic CT images from MRI for treatment planning and patient positioning using a 3-channel U-net trained on sagittal images | |
Kontaxis et al. | A new methodology for inter-and intrafraction plan adaptation for the MR-linac | |
US9507886B2 (en) | Multi-objective radiation therapy optimization method | |
Yepes et al. | Comparison of Monte Carlo and analytical dose computations for intensity modulated proton therapy | |
Chen et al. | A fast optimization algorithm for multicriteria intensity modulated proton therapy planning | |
US20080317317A1 (en) | Method and Apparatus For Accelerated Elastic Registration of Multiple Scans of Internal Properties of a Body | |
JP2022534324A (en) | Radiotherapy planning using a differentiable dose function | |
Tsekas et al. | DeepDose: a robust deep learning-based dose engine for abdominal tumours in a 1.5 T MRI radiotherapy system | |
CN110246562A (en) | Determine the method, apparatus and computer system of the sub-beam intensity in radiotherapy system | |
Kusch et al. | A robust collision source method for rank adaptive dynamical low-rank approximation in radiation therapy | |
Pastor-Serrano et al. | Millisecond speed deep learning based proton dose calculation with Monte Carlo accuracy | |
CN108771794A (en) | System and method for generating the Rapid Dose Calculation of radiation treatment in the works | |
Zhu et al. | Search for IMRT inverse plans with piecewise constant fluence maps using compressed sensing techniques | |
EP3445449B1 (en) | Fractionation selection tool in radiotherapy planning | |
Jhanwar et al. | Domain knowledge driven 3D dose prediction using moment-based loss function | |
Neylon et al. | A nonvoxel‐based dose convolution/superposition algorithm optimized for scalable GPU architectures | |
CN111494815B (en) | Three-dimensional dose calculation method, device and medium based on mixed variable-scale model | |
Matsuura et al. | Predictive gamma passing rate of 3D detector array-based volumetric modulated arc therapy quality assurance for prostate cancer via deep learning | |
Li et al. | A new Monte Carlo-based treatment plan optimization approach for intensity modulated radiation therapy | |
Vazquez et al. | A deep learning-based approach for statistical robustness evaluation in proton therapy treatment planning: a feasibility study | |
Sauget et al. | Efficient domain decomposition for a neural network learning algorithm, used for the dose evaluation in external radiotherapy | |
Lee et al. | Automated gamma knife radiosurgery treatment planning with image registration, data‐mining, and Nelder‐Mead simplex optimization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |