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
- intensity
- fields
- objective function
- tissue
- field
- 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
- 238000001959 radiotherapy Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 65
- 230000006870 function Effects 0.000 claims description 101
- 230000005855 radiation Effects 0.000 claims description 75
- 238000004422 calculation algorithm Methods 0.000 claims description 28
- 238000005457 optimization Methods 0.000 claims description 16
- 238000011478 gradient descent method Methods 0.000 claims description 10
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 5
- 238000002721 intensity-modulated radiation therapy Methods 0.000 abstract description 8
- 230000008569 process Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 11
- 238000004590 computer program Methods 0.000 description 6
- 238000002591 computed tomography Methods 0.000 description 5
- 230000014509 gene expression Effects 0.000 description 5
- 238000013507 mapping Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 238000002939 conjugate gradient method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 210000000920 organ at risk Anatomy 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000002922 simulated annealing Methods 0.000 description 2
- 201000011510 cancer Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 239000011541 reaction mixture Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
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 present application relates generally to methods and apparatus and computer systems for determining beamlet intensities in a radiotherapy system.
Background
Radiotherapy is one of the main technical means for treating malignant tumors, and aims to irradiate a tumor (tissue to be treated with radiotherapy or a target region) with radiation at a maximum concentration and to receive a minimum dose of surrounding normal tissues as much as possible. IMRT (intensity modulated radiation therapy), developed from three-dimensional conformal radiation therapy, is one type of radiation therapy. IMRT requires that the intensity of the sub-beams (also called beamlets) within the field be adjusted to certain requirements. The method is characterized in that the intensity of the sub-beams is adjusted according to the three-dimensional shape of the target area and the specific anatomical relationship between the vital organs and the target area under the condition that the fields are consistent with the shape of the target area, and the dose contribution in a single field is uneven but is more uniform in the whole volume of the target area than three-dimensional conformal treatment.
Each field is divided into a plurality of beamlets. In planning, these beamlets are assigned different weights according to the three-dimensional shape of the target and the anatomical relationship with the associated organs-at-risk to produce an optimized, non-uniform intensity contribution within the same field of view such that the beam flux through the organs-at-risk is reduced while the beam flux in other portions of the target is increased.
Currently, the inverse algorithm is adopted clinically to make an intensity modulated radiotherapy plan. The number, shape and weight (corresponding to the intensity of the sub-beams) of each field are obtained iteratively by an optimization algorithm according to the requirement of the user on the dose for uniform irradiation of the tumor and the requirement on the dose for protecting the surrounding normal tissue organs, and the requirement of the user on the position and number of the fields, the target dose on the tumor and the dose constraint on the organs at risk, which are specified empirically. However, since the positions and the number of the fields are empirically determined by the user, the intensity of the sub-beams in each field cannot be optimized. Therefore, it is particularly important to reduce the number of beamlets while ensuring compliance with the final dose distribution.
Disclosure of Invention
In view of the above, the present application provides a method, apparatus and computer system capable of selecting a spare field from a plurality of fields to more accurately determine beamlet intensities in a radiotherapy system.
According to an aspect of the application, a method of determining beamlet intensities in a radiotherapy system, the radiotherapy system comprising a plurality of fields, each field having a plurality of beamlets, the method comprising: 1) dividing tissue to be radiotreated and surrounding tissue around the tissue to be radiotreated into a plurality of voxels and determining a radiation dose for each of the plurality of voxels; 2) simulating a first unit dose contribution of each beamlet of the plurality of beamlets to each voxel at unit intensity based on the position of each field and beamlet relative to the tissue to be treated and surrounding tissue and the tissue information of the tissue to be treated and surrounding tissue; 3) establishing a first objective function based on the radiation dose and the first unit dose contribution; 4) selecting a spare field from the plurality of fields by optimizing a first objective function; 5) determining a second unit dose contribution of each beamlet in the spare field to each voxel at unit intensity from the selected spare field and the tissue information; 6) constructing a second objective function based on the second unit dose contribution and the radiation dose; 7) the intensity of the beamlets in the standby field is determined by optimizing a second objective function.
According to an embodiment of the present application, in step 2) and step 5), the first unit dose contribution and the second unit dose contribution may be modeled by a monte carlo method, respectively.
According to an embodiment of the present application, the first objective function may be constructed as:
wherein p satisfies the constraint condition that p is more than or equal to 0 and less than or equal to 1,corresponds to the radiation dose, A1The contribution of the first unit dose is expressed,corresponds to a sub-beam intensity, x, in a plurality of fieldsijAn intensity of a jth sub-beam representing an ith field of the plurality of fields,representing the maximum beamlet intensity in the ith field of the plurality of fields, β is a weight coefficient.
According to an embodiment of the present application, in the case where p ═ 1, the first objective function may be optimized by a newton method or a gradient descent method.
According to an embodiment of the present application, in step 4), in the case where p ═ 1, the first objective function may be converted into:to optimize a first objective function, where yiSatisfy yi≥xijWith the proviso that T isQuadratic form of the term.
In the transformed first objective function, p satisfies the constraint of 0. ltoreq. p.ltoreq.1,corresponds to the radiation dose, A1The contribution of the first unit dose is expressed,corresponds to a sub-beam intensity, x, in a plurality of fieldsijIntensity, y, of a jth sub-beam representing an ith field of the plurality of fieldsiRepresenting the maximum beamlet intensity in the ith field of the plurality of fields, β is a weight coefficient.
According to an embodiment of the present application, in step 4), the second objective function may be constructed as:
wherein A is2The contribution of the second unit dose is expressed,each element in (a) corresponds to one of the spare fieldsThe intensity of the sub-beams,corresponds to the radiation dose.
According to an embodiment of the present application, in step 4), the first objective function may be further configured to:
wherein p satisfies the constraint condition that p is more than or equal to 0 and less than or equal to 1,is the radiation dose of one of the voxels corresponding to the tissue to be radiotreated, APTVRepresenting the unit dose contribution of each sub-beam at unit intensity to each voxel corresponding to the tissue to be radiotherapy treated, N representing the amount of said surrounding tissue, AOARkRepresenting the unit dose contribution of each beamlet at unit intensity to each voxel corresponding to the kth surrounding tissue,corresponds to a sub-beam intensity, x, in a plurality of fieldsijAn intensity of a jth sub-beam representing an ith field of the plurality of fields,represents the maximum beamlet intensity in the ith field of the plurality of fields, β and λkAre weight coefficients.
According to an embodiment of the present application, in step 4), the first objective function may be optimized by a global optimization algorithm.
According to another aspect of the present application, there is provided an apparatus for determining beamlet intensities in a radiotherapy system, the radiotherapy system comprising a plurality of fields, each field having a plurality of beamlets, the apparatus comprising a dose determination unit, a spare field selection unit and an intensity determination unit. The dose determination unit may be configured to divide the tissue to be radiotherapy and surrounding tissue surrounding the tissue to be radiotherapy into a plurality of voxels and to determine a radiation dose for each of the plurality of voxels. The spare field selection unit may be configured to: simulating a first unit dose contribution of each beamlet of the plurality of beamlets to each voxel at unit intensity based on the position of each field and beamlet relative to the tissue to be treated and surrounding tissue and the tissue information of the tissue to be treated and surrounding tissue; establishing a first objective function based on the radiation dose and the first unit dose contribution; and selecting a spare field from the plurality of fields by optimizing the first objective function. The intensity determination unit may be configured to: determining a second unit dose contribution of each beamlet in the spare field to each voxel at unit intensity from the selected spare field and the tissue information; constructing a second objective function based on the second unit dose contribution and the radiation dose; and determining the intensity of the beamlets in the standby field by optimizing the second objective function.
According to yet another aspect of the present application, there is provided a computer system for determining beamlet intensities in a radiotherapy system, the radiotherapy system comprising a plurality of fields, each field having a plurality of beamlets, the system for determining beamlet intensities in the radiotherapy system comprising: a memory storing computer instructions; a processor executing computer instructions stored by the memory to perform the following operations: dividing tissue to be radiotreated and surrounding tissue around the tissue to be radiotreated into a plurality of voxels and determining a radiation dose for each of the plurality of voxels; simulating a first unit dose contribution of each beamlet of the plurality of beamlets to each voxel at unit intensity based on the position of each field and beamlet relative to the tissue to be treated and surrounding tissue and the tissue information of the tissue to be treated and surrounding tissue; establishing a first objective function based on the radiation dose and the first unit dose contribution; selecting a spare field from the plurality of fields by optimizing a first objective function; determining a second unit dose contribution of each beamlet in the spare field to each voxel at unit intensity from the selected spare field and the tissue information; constructing a second objective function based on the second unit dose contribution and the radiation dose; the intensity of the beamlets in the standby field is determined by optimizing a second objective function.
Drawings
These and/or other aspects of the present application will become apparent and more readily appreciated from the following description of the exemplary embodiments, taken in conjunction with the accompanying drawings of which:
figure 1 is a flow chart of a method of determining beamlet intensities in a radiotherapy system according to an exemplary embodiment of the present application;
fig. 2 is a schematic plan view of a beam field including a plurality of beamlets according to an exemplary embodiment of the present application;
FIG. 3 is a schematic perspective view of irradiation of tissue to be radiotreated and surrounding tissue through multiple fields in accordance with an exemplary embodiment of the present application;
FIG. 4 is a graph of the dose contribution of a spare field selected by the method of determining beamlet intensities in a radiotherapy system of the present application to tissue to be treated and surrounding tissue in an intensity modulated radiotherapy method;
FIG. 5 is a schematic diagram of the dose contribution of a spare field empirically selected by a user to the tissue to be treated and surrounding tissue in an intensity modulated radiation therapy method;
figure 6 is a block diagram illustrating an apparatus for determining beamlet intensities in a radiotherapy system according to an embodiment of the present application; and
FIG. 7 is a block diagram illustrating a computer system suitable for implementing embodiments of the present application.
Detailed Description
While the disclosure is susceptible to various modifications and alternative embodiments, specific embodiments have been shown in the drawings and will be described in detail in this written description. Aspects and features of the present disclosure and methods or methods of accomplishing the same will become apparent when reference is made to the embodiments described with reference to the accompanying drawings. This disclosure may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein.
As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items. Furthermore, when describing embodiments of the present invention, the use of "may" refer to "one or more embodiments of the present invention. Expressions such as "at least one of" when following a column of elements modify the entire column of elements without modifying individual elements in the column. Additionally, the term "exemplary" is intended to mean exemplary or illustrative. The use of the singular forms "a", "an" and "the" encompass plural referents unless the context clearly dictates otherwise.
As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises," "comprising," "includes" and/or "including," when used herein, specify the presence of stated features or components, but do not preclude the presence or addition of one or more other features or components.
As used herein, "substantially," "about," "approximately," and similar terms are used as terms of approximation and not as terms of degree, and are intended to account for inherent deviations in measured or calculated values that would be recognized by one of ordinary skill in the art. Further, these terms, as used herein, include the set point and average value within an acceptable range of deviation for a particular value as determined by one of ordinary skill in the art in view of the measurement in question and the error associated with the measurement of the particular value (i.e., the limitations of the measurement system). For example, "about" can mean within one or more standard deviations, or within ± 30%, ± 20%, ± 10%, ± 5% of the stated values.
While certain embodiments may be implemented differently, the particular process sequence may be performed differently than described. For example, two consecutively described processes may be performed at approximately the same time, or in the reverse order to the described order.
The method and apparatus for determining the intensity of beamlets in a radiotherapy system of the present application will be described in detail with reference to fig. 1 to 6.
In a method of determining beamlet intensities in a radiotherapy system according to the present application, for example, the radiotherapy system may comprise a plurality of fields F (also referred to as fields) as shown in fig. 3, and each field has a plurality of beamlets B. The source F shown in fig. 2 comprises 9 beamlets B. For example, each field F is realized as a plurality of beamlets B by an MLC (Multi-Leaf collimator). The field F shown in FIG. 2 can be achieved by a radiation source and MCL included in the radiotherapy system. The source may emit X-rays or gamma-rays, i.e. the field F is a field of X-rays. As another example, the source may emit positively charged ions (i.e., protons). This is merely exemplary and the source according to the application is not limited thereto. In one example, the number of beamlets B in each field F may be the same or different, and each beamlet B may have a different shape.
Referring to fig. 1 and 3, in step S1, a tissue to be radiotherapy PTV and a surrounding tissue OAR around the tissue to be radiotherapy are divided into a plurality of voxels and a radiation dose of each of the plurality of voxels is determined. In the example shown in fig. 3, the tissue to be radiotherapy PTV is divided into a plurality of voxels, and the surrounding tissue OAR is divided into a plurality of voxels surrounding the tissue to be radiotherapy PTV. That is, the voxels divided by the tissue PTV to be treated with radiation and the surrounding tissue OAR may correspond to equal three-dimensional volumes. A voxel is the smallest unit of digital data on a three-dimensional partition. In fig. 3, 20 fields F are shown, but this is merely exemplary. The radiation dose for each of the plurality of voxels represents the radiation dose that each voxel needs to receive. For example, the radiation dose to be received by the tissue PTV to be treated is greater than a predetermined value, and the radiation dose to be received by the surrounding tissue OAR is less than a predetermined value. The radiation dose for each of the plurality of voxels may be determined from a CT image of the tissue PTV and surrounding tissue OAR to be radiotreated.
In step S2, the first unit dose contribution of each beamlet B in the plurality of fields F to each voxel is modeled based on the position of each field F and beamlet B relative to the tissue PTV and surrounding tissue OAR to be treated and the tissue information of the tissue PTV and surrounding tissue OAR to be treated.
The tissue information of the tissue PTV to be radiotherapy and the surrounding tissue OAR can be acquired by a computed tomography image (i.e., CT image) of the tissue PTV to be radiotherapy and the surrounding tissue OAR. For example, the tissue information may be acquired by a planning CT (computed tomography) device for inverse intensity modulated radiation therapy. The tissue information may include volume, shape, location and electron density information of the tissue to be treated for radiotherapy PTV and surrounding tissue OAR. The position of the plurality of fields F and beamlets B relative to the tissue PTV to be treated and the surrounding tissue OAR may typically be set manually. For example, the field F may be set so that the rays in the field F pass through the gaps between the voxels corresponding to the surrounding tissue OAR as much as possible and irradiate on the tissue PTV to be radiotreated. However, no matter how the position of the field F is set, the dose received by the OAR of the surrounding tissue cannot be zero.
The first unit dose contribution is the dose contribution of each beamlet B to each voxel at unit intensity.
In step S3, a first objective function is established based on the determined radiation dose and the first unit dose contribution. In step S4, a spare field P is selected from the plurality of fields F by optimizing the first objective function.
In inverse intensity modulated radiotherapy, with the radiation dose and first unit dose contribution determination, the main objective is to determine the intensities of all the beamlets B in the multiple fields F such that each voxel receives a dose that satisfies as much as possible the radiation dose it should receive.
Thus, the objective function used in the backprojection radiotherapy planning is:
wherein,corresponds to the radiation dose determined in step S1, a denotes the first unit dose contribution, vector, determined in step S2Corresponds to the intensity, vector, of one beamlet B of the plurality of fields F in steps S1 and S2It may be a vector in which the intensities of all the beamlets B are connected in series by a row. Solution of the conventional objective functionIs such that the vectorThe second-order norm of (a) takes the vector of all beamlets B whose minimum value is connected in series by rows.
Thus, the first objective function according to the present application introduces the position of the field F and the beamlet B into the conventional objective function.
According to an embodiment of the present application, the field F and the position of the beamlet B may be introduced by introducing the intensity of the jth beamlet B of the ith field of the plurality of fields F in a conventional objective function. In other words, F is further determined based on the mapping relationship of the intensity of the jth sub-beam B of the ith field to the corresponding field F and sub-beam B, indirectly introducing the positions of field F and sub-beam B.
Further, unlike the conventional intensity modulated radiotherapy in which a doctor empirically determines a field to be used (hereinafter referred to as a spare field) from a plurality of fields F, the present application applies compressive sensing to inverse intensity modulated radiotherapy to select a more accurate spare field P from the plurality of fields F, thereby determining the position of the spare field P and the corresponding beamlets B included therein with respect to the tissue PTV to be radiotherapy and the surrounding tissue OAR.
Since the fundamental purpose of establishing the first objective function is to select several spare fields P from the plurality of fields F, i.e., to reduce the number of fields F while satisfying the radiation dose to be delivered.
Therefore, in the method of the present application, based on the objective function (i.e. equation 1) used in the inverse intensity modulated radiotherapy plan, based on the compressive sensing theory, the intensity set of the sub-beams is set such that the product of the sum of the intensity maxima of the sub-beams B in each field F and the weighting coefficient takes its minimum value (which may be a preset value) as a limiting condition, the intensity maxima of the sub-beams B in each field F are taken as sampling signals, and then the spare field P is selected from the field F by solving the sparse solution of the sampling signals through the compressive sensing, wherein the intensity set of the sub-beams in the spare field P simultaneously enables the radiation dose finally received by the tissue PTV to be irradiated and the surrounding tissue OAR to meet or most closely meet the intensity modulated radiotherapy target (i.e. the preset radiation dose). Therefore, the method according to an embodiment of the present application may also be considered as an application of the compressed sensing method in signal processing.
In one aspect, compressed sensing is a technique for reconstructing sparse or compressible signals. The mathematical manifestation of a sparse signal is that it contains many coefficients close to or equal to zero when expressed in some domains. That is, a signal is said to be sparse if it has only a finite non-zero sample point, while other sample points are zero (or close to zero). Efficient use of compressed sensing techniques (i.e., sparsity of the signal to be processed) can be used to reduce the number of measurement samples required. Mathematically, intensity modulated radiation therapy inverse planning is similar to the signal processing problem described above in the case where multiple fields F (i.e., all existing fields F) are the "signals" to be sought. As mentioned above, inverse planning is a pathological problem and there are usually a plurality of sets of spare fields P (i.e. a set of fields selected from the plurality of fields F, for example selected empirically by a physician) that can generate the radiation dose that the tissue PTV to be treated with radiation needs to receive, but in this case it is not necessarily possible to receive the minimum radiation dose to the surrounding tissue OAR. Therefore, in the present application, the spare field P is selected from the plurality of fields F by solving for sparseness of the sampling signals by using the maximum intensity value of the beamlet B in each field F as the sampling signal, and the set of the beamlet intensities in the spare field P enables the radiation dose finally received by the tissue PTV to be radiotherapy and the surrounding tissue OAR to satisfy or most approximate to the target of intensity modulated radiotherapy as the limiting condition (for example, equation (1) is satisfied). The number of non-zero values in the sparse solution found corresponds to the number of spare fields P.
On the other hand, when the maximum intensity value of the sub-beam B in each field F is used as the sampling signal, each term in the sparse solution obtained by the compressed sensing method may correspond to the maximum sub-beam intensity in one field F, so that a term of zero in the sparse solution indicates that the intensities of all the sub-beams B in the corresponding field F are all zero, that is, the field F corresponding to the term is not selected as the spare field P. In this way, a one-to-one mapping relationship of each item in the sparse solution with the position of the portal F can be realized.
Further, the number of candidate fields P can be set by setting the weight coefficient, and for example, the larger the weight coefficient, the smaller the number of candidate fields P.
According to one embodiment of the present application, the first objective function component may be:
wherein p satisfies the constraint condition that p is more than or equal to 0 and less than or equal to 1,corresponds to the radiation dose of one of the plurality of voxels, a denotes the first unit dose contribution described above,each element in (a) corresponds to a beamlet intensity,it may be a vector in which the intensities of all the beamlets B are connected in series by a row. x is the number ofijAn intensity of a jth sub-beam representing an ith field of the plurality of fields F,representing the maximum beamlet intensity in the ith field of the plurality of fields F, β is a weight coefficient.
The number of candidate fields P selected from the plurality of fields F can be set by setting the value of the weight coefficient β the larger the value of β, the smaller the number of candidate fields P, the value of β is larger than zero.
According to an embodiment of the present application, the first objective function (e.g., equation (2)) may be optimized by a global optimization algorithm to obtain a vectorVector quantityCorresponds to the intensity of the respective beamlet B in the respective field F. That is, the intensity of each beamlet B corresponds to a unique beamlet B in a unique field F. Thus, the vector can be passedEach element in (1) is related to the field F andthe mapping relation of the sub-beam B in the field F selects a spare field P from the plurality of fields F.
According to one embodiment of the present application, it is possible to determine whether the radiation dose to be finally achieved is closer to the radiation dose of the tissue to be treated PTV or the radiation dose of the surrounding tissue OAR by differentiating the voxels of the tissue to be treated PTV and the surrounding tissue OAR in the first objective function and setting the weighting coefficients.
Thus, according to another exemplary embodiment of the present application, the first objective function may be further structured as:
wherein p satisfies the constraint condition that p is more than or equal to 0 and less than or equal to 1,is the radiation dose, A, of one of the voxels corresponding to the tissue to be radiotherapy PTVPTVRepresenting the unit dose contribution of each beamlet B at unit intensity to each voxel corresponding to the tissue PTV to be irradiated, N representing the amount of said surrounding tissue, AOARkRepresenting the unit dose contribution of each beamlet B at unit intensity to each voxel corresponding to the k-th surrounding tissue OAR,corresponds to the intensity of the beamlet B,may be a vector, x, in which the intensities of all beamlets B are connected in series by a lineijAn intensity of a jth sub-beam representing an ith field of the plurality of fields F,representing the ith ray in a plurality of fields FThe maximum beamlet intensity in the field, β, is the weight coefficient, λkIs a weight coefficient
For example, in the examples shown in fig. 3 to 5, the number of N is 6, i.e., the number of the surrounding tissues OAR is 6. Can be set by setting lambdakDetermines whether the last selected spare field P is such that the radiation dose is more than that required for the tissue PTV to be treated or that required for the surrounding tissue OAR.
According to an embodiment of the present application, the first objective function (e.g., equation (3)) may be optimized by a global optimization algorithm to obtain a vectorFurther, a plurality of spare fields P are selected from the plurality of fields F.
In the case where p in the first objective function (e.g., equation (2) or equation (3)) is equal to 1, the first objective function (e.g., equation (2) or equation (3)) may be optimized through a conventional iterative algorithm such as a newton method, a gradient descent method, a conjugate gradient method, or a quasi-newton method.
That is, in the case where 0 ≦ p ≦ 1, the first objective function (e.g., equation (2) or equation (3)) may be optimized by the global optimization algorithm to obtain a globally optimal solution for the objective function. The global optimization researches the characteristics and the construction of the global optimal solution of the multivariable nonlinear function in a certain constraint area to find the computing method of the global optimal solution, and the theoretical properties and the computing performance of the solving method. For example, the first objective function may be optimized using a global optimization algorithm such as simulated annealing algorithm, genetic algorithm, tabu search, particle swarm algorithm, ant colony algorithm. But the application is not limited thereto.
In the case where p is equal to 1 in the first objective function (e.g., equation (2) or equation (3)), the first objective function (e.g., equation (2) or equation (3)) may be optimized by an iterative algorithm such as a newton method, a gradient descent method, a conjugate gradient method, or a quasi-newton method.
At the rootIn the method according to an embodiment of the present application, in the case where p in the first objective function (e.g., expression (2) or expression (3)) is equal to 1, within a norm of order 1 in expression (6) or expression (7)The term is used as a constraint to convert the non-linear problem into a linear problem to further reduce the amount of computation.
Taking equation (2) as an example, the first objective function is converted into:
wherein, yiSatisfy yi≥xijWith the proviso that T isQuadratic form of the term.
In the formula (4), the reaction mixture is,corresponds to the radiation dose of one of the plurality of voxels, a denotes the first unit dose contribution described above,each element in (a) corresponds to a beamlet intensity,may be a vector, x, in which the intensities of all beamlets B are connected in series by a lineijAn intensity of a jth sub-beam representing an ith field of the plurality of fields F,representing the maximum beamlet intensity in the ith field of the plurality of fields F, β is a weight coefficient.
In this way, the objective function can be optimized by an iterative algorithm such as newton's method, gradient descent method, etc. with a reduced amount of computation to obtain a vectorFurther, the vector may be passedThe mapping relation between each element in (B) and the field F and the sub-beam B in the field F, and selects a spare field P from the plurality of fields F.
Referring to fig. 3, through steps S3 and S4, a plurality of fields F (e.g., 3 fields F) can be selected as spare fields P from the 20 fields F shown in fig. 3.
In step S5, a second unit dose contribution of each beamlet B in the spare field P to each voxel at unit intensity is determined from the selected spare field P and the tissue information. For example, the second unit dose contribution is modeled based on the position of each field F and beamlet B relative to the tissue to be treated PTV and surrounding tissue OAR and tissue information of the tissue to be treated PTV and surrounding tissue OAR. The second unit dose contribution is simulated, for example, by the monte carlo method. Since the first objective function (step S3) is established by considering the set of beamlet intensities in the spare field P as the limiting condition, so that the radiation dose finally received by the tissue PTV to be treated and the surrounding tissue OAR meets or is closest to the intensity modulated radiation treatment target, the second unit dose contribution modeled by the spare field P can make the inverse planning result closer to the ideal, i.e. the finally obtained beamlet intensities can make the surrounding tissue OAR receive the minimum radiation dose while making the tissue PTV to be treated receive the preset radiation dose. In step S6, a second objective function is constructed from the second unit dose contribution and the radiation dose.
In an embodiment of the present application, the second objective function may correspond to a conventional objective function in an intensity modulated radiotherapy plan, that is:
wherein A is2Represents the contribution of the second unit dose,corresponds to one sub-beam intensity in the alternative field P,it may be a vector in which all the beamlet intensities in the candidate field P are connected in series by a row,corresponds to the radiation dose.
In step S7, the intensity of the beamlet B in the spare field P is determined by optimizing the second objective function. In the case where the second objective function is equation (5), equation (5) can be optimized by newton's method or gradient descent method. As another example of the present application, equation (5) may be optimized using a global optimization algorithm. But the application is not limited thereto.
Fig. 4 shows a schematic diagram of the dose contribution of the spare field selected by the method of determining the intensity of the sub-beams in a radiotherapy system of the present application to the tissue to be radiotherapy PTV and the surrounding tissue OAR. Fig. 5 is a schematic diagram of the dose contribution of the spare field P empirically selected by the user to the tissue PTV to be treated and the surrounding tissue OAR in the intensity modulated radiation therapy method.
In fig. 4 and 5, the radiation dose received by the voxel is schematically represented by the number of lines in each tissue, i.e. a larger number of lines indicates a larger radiation dose received by the corresponding tissue. 3 from the 20 fields F shown in fig. 3 are selected as the spare fields P by steps S3 and S4 and the intensity of the sub-beam B is acquired based on the spare fields P again by steps S5 to S7. The radiation dose to the surrounding tissue OAR in fig. 4 is significantly reduced compared to the radiation dose to the surrounding tissue OAR in fig. 5.
In equations (2) to (9), the number of candidate fields P selected from the plurality of fields F can be set by setting the value of β the larger the value of β, the smaller the number of candidate fields P, the value of β can be larger than zero.
The methods described with reference to fig. 1-5 may be implemented by a computer system. The computer system includes a memory storing executable instructions and a processor. The processor communicates with the memory to execute the executable instructions to implement the methods described with reference to fig. 1-5. Alternatively or additionally, the methods described with reference to fig. 1-5 may be implemented by a non-transitory computer storage medium. The medium stores computer readable instructions that, when executed, cause a processor to perform the method described with reference to fig. 1-5.
Figure 6 is a block diagram illustrating an apparatus 100 for determining beamlet intensities in a radiotherapy system according to an embodiment of the present application. An apparatus 100 for determining beamlet intensities in a radiotherapy system comprising: a dose determination unit 110 configured to divide a tissue to be radiotherapy PTV and a surrounding tissue OAR around the tissue to be radiotherapy into a plurality of voxels and determine a radiation dose of each of the plurality of voxels; a spare field selection unit 120 configured to: simulating a first unit dose contribution of each of the beamlets in a plurality of fields to each voxel at a unit intensity based on the position of each field F and beamlet B with respect to the tissue to be treated and the surrounding tissue and tissue information of the tissue to be treated PTV and surrounding tissue OAR; establishing a first objective function based on the radiation dose and the first unit dose contribution; selecting a standby radiation field P from the plurality of radiation fields F by optimizing the first objective function; an intensity determination unit configured to: determining a second unit dose contribution of each beamlet B in the spare field P to each voxel at a unit intensity from the selected spare field P and the tissue information; constructing a second objective function based on the second unit dose contribution and the radiation dose; the intensity of the sub-beam B in the backup field P is determined by means of the second objective function.
According to an embodiment of the present application, the spare field selection unit 120 and the intensity determination unit 130 may simulate the first unit dose contribution and the second unit dose contribution by the monte carlo method, respectively.
According to an embodiment of the present application, the spare field selection unit 120 constructs the first objective function as:
in the formula (6), p satisfies the constraint of 0. ltoreq. p.ltoreq.1,corresponds to the radiation dose of one of the plurality of voxels, a denotes the first unit dose contribution described above,each element in (a) corresponds to a beamlet intensity,it may be a vector in which the intensities of all the beamlets B are connected in series by a row. x is the number ofijAn intensity of a jth sub-beam representing an ith field of the plurality of fields F,representing the maximum beamlet intensity in the ith field of the plurality of fields F, β is a weight coefficient.
According to an embodiment of the present application, the first objective function (e.g., equation (6)) may be optimized by a global optimization algorithm to obtain a vectorVector quantityCorresponds to the intensity of the respective beamlet B in the respective field F. That is, the intensity of each beamlet B corresponds to a unique beamlet B in a unique field F. Thus, the vector can be passedThe mapping relation between each element in (B) and the field F and the sub-beam B in the field F, and selects a spare field P from the plurality of fields F.
According to an embodiment of the present application, the spare field selection unit 120 may further configure the first objective function to:
wherein p satisfies the constraint condition that p is more than or equal to 0 and less than or equal to 1,is the radiation dose, A, of one of the voxels corresponding to the tissue to be radiotherapy PTVPTVRepresenting the unit dose contribution of each beamlet B at unit intensity to each voxel corresponding to the tissue PTV to be irradiated, N representing the amount of said surrounding tissue, AOARkRepresenting the unit dose contribution of each beamlet B at unit intensity to each voxel corresponding to the k-th surrounding tissue OAR,corresponds to the intensity of the beamlet B,may be a vector, x, in which the intensities of all beamlets B are connected in series by a lineijRepresenting in a plurality of fields FThe intensity of the jth beamlet of the ith field,denotes the maximum beamlet intensity in the ith field of the plurality of fields F, β is a weight coefficient, λkAre weight coefficients.
According to an embodiment of the present application, the first objective function (e.g., equation (7)) may be optimized by a global optimization algorithm to obtain a vectorFurther, a plurality of spare fields P are selected from the plurality of fields F.
That is, in the case where 0 ≦ p ≦ 1, the spare field selection unit 120 may optimize the first objective function (e.g., equation (6) or equation (7)) through a global optimization algorithm to obtain a global optimal solution for the objective function. The global optimization researches the characteristics and the construction of the global optimal solution of the multivariable nonlinear function in a certain constraint area to find the computing method of the global optimal solution, and the theoretical properties and the computing performance of the solving method. For example, the first objective function may be optimized using a global optimization algorithm such as simulated annealing algorithm, genetic algorithm, tabu search, particle swarm algorithm, ant colony algorithm. But the application is not limited thereto.
In the case where p in the first objective function (e.g., equation (6) or equation (7)) is equal to 1, the first objective function (e.g., equation (6) or equation (7)) may be optimized through an iterative algorithm such as a newton method, a gradient descent method, a conjugate gradient method, or a quasi-newton method. According to an embodiment of the present application, in the case where p is equal to 1 in the first objective function (e.g., equation (6) or equation (7)), the spare field selection unit 120 may optimize the first objective function (e.g., equation (6) or equation (7)) through a conventional iterative algorithm such as a newton method, a gradient descent method, a conjugate gradient method, or a quasi-newton method.
According to one embodiment of the present application, in the case where p is equal to 1, the 1 st order in formula (6) or formula (7) is normalizedWithin a number ofThe term is used as a constraint to convert the non-linear problem into a linear problem to further reduce the amount of computation.
Taking equation (6) as an example, the spare field selection unit 120 converts the first objective function into:
in the transformed first objective function (i.e., equation (8)), T isQuadratic form of the term, yiSatisfy yi≥xijThe limit conditions of (a), that is,
in the transformed first objective function,corresponds to the radiation dose of one of the plurality of voxels, a denotes the first unit dose contribution described above,each element in (a) corresponds to a beamlet intensity,may be a vector, x, in which the intensities of all beamlets B are connected in series by a lineijAn intensity of a jth sub-beam representing an ith field of the plurality of fields F,express a plurality ofThe largest beamlet intensity in the ith of the individual fields F, β, is the weighting factor.
According to an embodiment of the present application, the intensity determination unit 130 may construct the second objective function as:
wherein A is2Represents the contribution of the second unit dose,corresponds to one sub-beam intensity in the alternative field P,it may be a vector in which all the beamlet intensities in the candidate field P are connected in series by a row,corresponds to the radiation dose.
In one embodiment of the present application, the intensity determination unit 130 may determine the intensity of the beamlet B in the spare field P by optimizing the second objective function. In the case where the second objective function is equation (9), equation (9) can be optimized by newton's method or gradient descent method. As another example of the present application, equation (9) may be optimized using a global optimization algorithm. But the application is not limited thereto.
The apparatus for determining beamlet intensities in a radiotherapy system described with reference to figure 6 may be implemented by a computer system. The computer system may include a memory storing executable instructions and a processor. The processor communicates with the memory to execute executable instructions to implement the apparatus described with reference to fig. 6. Alternatively or additionally, the apparatus described with reference to fig. 6 may be embodied by a non-transitory computer storage medium. The medium stores computer readable instructions that, when executed by a computer, perform the functions of the apparatus described with reference to fig. 6.
Referring now to FIG. 7, FIG. 7 is a block diagram that illustrates a computer system 7000 that is suitable for implementing embodiments of the present application.
As shown in fig. 7, the computer system 7000 may include a processor (such as a Central Processing Unit (CPU)7001, a Graphic Processing Unit (GPU), or the like) that can perform various appropriate actions and processes in accordance with a program stored in a Read Only Memory (ROM)7002 or a program loaded from a storage section 7008 into a Random Access Memory (RAM) 7003. In the RAM 7003, various programs and data required for the operation of the system 7000 are also stored. The CPU 7001, ROM 7002, and RAM 7003 are connected to one another through a bus 7004. An input/output I/O interface 7005 is also connected to the bus 7004.
The following are components that may be connected to the I/O interface 7005: an input portion 7006 including a keyboard, a mouse, and the like; an output portion 7007 including a cathode ray tube CRT, a liquid crystal display device LCD, a speaker, and the like; a storage portion 7008 including a hard disk and the like; and a communication portion 7009 including a network interface card (such as a LAN card and a modem). The communication portion 7009 can perform communication processing over a network such as the internet. The driver 7010 may also be connected to the I/O interface 7005 as necessary. A removable medium 7011 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, or the like may be mounted on the drive 7010 so that a computer program read out therefrom is installed into the storage portion 7008 as needed.
In particular, according to embodiments of the present disclosure, the methods described above with reference to fig. 1 to 5 may be implemented as computer software programs. For example, embodiments of the disclosure may include a computer program product comprising a computer program tangibly embodied in a machine-readable medium. The computer program comprises program code for performing the method described with reference to fig. 1 to 5. In such an embodiment, the computer program can be downloaded from a network through the communication section 709 and installed, and/or can be installed from the removable medium 7011.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The units or modules referred to in the embodiments of the present application may be implemented by software or hardware. The described units or modules may also be provided in a processor. The names of these units or modules should not be construed as limiting these units or modules.
According to the embodiments of the present application, the provided method, apparatus and computer system for determining the intensity of sub-beams in a radiotherapy system can be used in any kind of radiotherapy system. For example, a radiotherapy system that effects radiotherapy by X-rays, gamma-rays, or protons, but the application is not limited thereto.
The above description is only exemplary of the present application and illustrative of the principles of the technology employed. It will be appreciated by a person skilled in the art that the scope of the present application is not limited to the embodiments with a specific combination of the above-mentioned features, but also covers other embodiments with any combination of the above-mentioned features or their equivalents without departing from the inventive concept. For example, the above features and the technical features having similar functions disclosed in the present application are mutually replaced to form the technical solution.
Claims (17)
1. A method of determining beamlet intensities in a radiotherapy system comprising a plurality of fields, each field having a plurality of said beamlets, the method comprising:
1) dividing tissue to be radiotreated and surrounding tissue surrounding the tissue to be radiotreated into a plurality of voxels and determining a radiation dose for each of the plurality of voxels;
2) modeling a first unit dose contribution of each of the beamlets in the plurality of fields to each of the voxels at a unit intensity as a function of the position of each of the fields and the beamlets relative to the tissue to be treated and the surrounding tissue and the tissue information for the tissue to be treated and the surrounding tissue;
3) establishing a first objective function based on the radiation dose and the first unit dose contribution;
4) selecting a spare field from the plurality of fields by optimizing the first objective function;
5) determining a second unit dose contribution of each of the beamlets in the back-up field to each of the voxels at a unit intensity from the selected back-up field and the tissue information;
6) constructing a second objective function from the second unit dose contribution and the radiation dose;
7) determining the intensity of the beamlets in the alternate field of radiation by optimizing the second objective function.
2. The method of claim 1, wherein in step 2) and step 5), the first unit dose contribution and the second unit dose contribution are modeled by a monte carlo method, respectively.
3. The method of claim 1 or 2, wherein the first objective function is constructed as:
wherein p satisfies the constraint condition that p is more than or equal to 0 and less than or equal to 1,corresponds to the radiation dose, A1Represents the first unit dose contribution and the second unit dose contribution,each element in (1)Corresponding to a sub-beam intensity, x, in the plurality of fieldsijAn intensity of a jth sub-beam representing an ith field of the plurality of fields,representing the maximum beamlet intensity in the ith of the plurality of fields, β is a weight coefficient.
4. The method of claim 3, wherein the first objective function is optimized by Newton's method or gradient descent method in case of p ═ 1.
5. The method of claim 4, wherein in step 4), the first objective function is transformed into:to optimize the first objective function,
wherein, yiSatisfy yi≥xijWith the proviso that T isQuadratic form of the term, and
wherein,corresponds to the radiation dose, A1Represents the first unit dose contribution and the second unit dose contribution,corresponds to a sub-beam intensity, x, of the plurality of fieldsijIntensity, y, of a jth sub-beam representing an ith field of the plurality of fieldsiRepresenting the maximum beamlet intensity in the ith of the plurality of fields, β is a weight coefficient.
6. The method of claim 5, wherein the second objective function is constructed as:
wherein A is2Represents the contribution of the second unit dose,corresponds to one sub-beam intensity in the spare field,corresponds to the radiation dose.
7. The method of claim 3, wherein the first objective function is further constructed as:
wherein p satisfies the constraint condition that p is more than or equal to 0 and less than or equal to 1,is the radiation dose of one of the voxels corresponding to the tissue to be radiotreated, APTVRepresenting the unit dose contribution of each of said sub-beams to each voxel corresponding to said tissue to be radiotherapy at a unit intensity, N representing the amount of said surrounding tissue, AOARkRepresenting a unit dose contribution of each of said sub-beams at a unit intensity to each voxel corresponding to a k-th of said surrounding tissue,corresponds to a sub-beam intensity, x, of the plurality of fieldsijAn intensity of a jth sub-beam representing an ith field of the plurality of fields,represents the maximum beamlet intensity in the ith of the plurality of fields, β and λkAre weight coefficients.
8. The method of claim 3, wherein in step 4), the first objective function is optimized by a global optimization algorithm.
9. An apparatus for determining beamlet intensities in a radiotherapy system comprising a plurality of fields, each field having a plurality of said beamlets, the apparatus comprising:
a dose determination unit configured to divide a tissue to be radiotherapy and a surrounding tissue around the tissue to be radiotherapy into a plurality of voxels and determine a radiation dose of each of the plurality of voxels;
a spare field selection unit configured to:
modeling a first unit dose contribution of each of the beamlets in the plurality of fields to each of the voxels at a unit intensity as a function of the position of each of the fields and the beamlets relative to the tissue to be treated and the surrounding tissue and the tissue information for the tissue to be treated and the surrounding tissue;
establishing a first objective function based on the radiation dose and the first unit dose contribution; and
selecting a spare field from the plurality of fields by optimizing the first objective function; and
an intensity determination unit configured to:
determining a second unit dose contribution of each of the beamlets in the back-up field to each of the voxels at a unit intensity from the selected back-up field and the tissue information;
constructing a second objective function from the second unit dose contribution and the radiation dose; and
determining the intensity of the beamlets in the alternate field of radiation by optimizing the second objective function.
10. The apparatus of claim 9, wherein the alternate field selection unit and the intensity determination unit simulate the first and second unit dose contributions, respectively, by a monte carlo method.
11. The apparatus of claim 9 or 10, wherein the alternate field selection unit constructs the first objective function as:
wherein p satisfies the constraint condition that p is more than or equal to 0 and less than or equal to 1,corresponds to the radiation dose, A1Represents the first unit dose contribution and the second unit dose contribution,corresponds to a sub-beam intensity, x, of the plurality of fieldsijAn intensity of a jth sub-beam representing an ith field of the plurality of fields,representing the maximum beamlet intensity in the ith of the plurality of fields, β is a weight coefficient.
12. The method of claim 11, wherein the alternate field selection unit optimizes the first objective function by newton's method or gradient descent method in case p ═ 1.
13. The apparatus of claim 12, wherein the alternate field selection unit converts the first objective function into:to optimize the first objective function,
wherein, yiSatisfy yi≥xijWith the proviso that T isQuadratic form of the term, and
wherein,corresponds to the radiation dose, A1Represents the first unit dose contribution and the second unit dose contribution,corresponds to a sub-beam intensity, x, of the plurality of fieldsijIntensity, y, of a jth sub-beam representing an ith field of the plurality of fieldsiRepresenting the maximum beamlet intensity in the ith of the plurality of fields, β is a weight coefficient.
14. The apparatus of claim 9, wherein the intensity determination unit constructs the second objective function as:
wherein A is2Represents the second unit agentThe contribution of the amount is such that,corresponds to one sub-beam intensity in the spare field,corresponds to the radiation dose.
15. The method of claim 11, wherein the alternate field selection unit further constructs the first objective function as:
wherein p satisfies the constraint condition that p is more than or equal to 0 and less than or equal to 1,is the radiation dose of one of the voxels corresponding to the tissue to be radiotreated, APTVRepresenting the unit dose contribution of each of said sub-beams to each voxel corresponding to said tissue to be radiotherapy at a unit intensity, N representing the amount of said surrounding tissue, AOARkRepresenting a unit dose contribution of each of said sub-beams at a unit intensity to each voxel corresponding to a k-th of said surrounding tissue,corresponds to a sub-beam intensity, x, of the plurality of fieldsijAn intensity of a jth sub-beam representing an ith field of the plurality of fields,represents the maximum beamlet intensity in the ith of the plurality of fields, β and λkAre weight coefficients.
16. The apparatus of claim 11, wherein the alternate field selection unit optimizes the first objective function by a global optimization algorithm.
17. A computer system for determining beamlet intensities in a radiotherapy system, the radiotherapy system comprising a plurality of fields, each field having a plurality of beamlets, the system for determining beamlet intensities in the radiotherapy system comprising:
a memory storing computer instructions;
a processor executing computer instructions stored by the memory to perform the method of any of claims 1 to 8.
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 (6)
| 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 |
| WO2024152913A1 (en) * | 2023-01-16 | 2024-07-25 | 迈胜医疗设备有限公司 | Method and apparatus for generating arc therapy plan, and device and storage medium |
| 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 | |
| Schuemann et al. | Site-specific range uncertainties caused by dose calculation algorithms for proton therapy | |
| Kontaxis et al. | Towards adaptive IMRT sequencing for the MR-linac | |
| Inaniwa et al. | Adaptation of stochastic microdosimetric kinetic model for charged-particle therapy treatment planning | |
| Marchant et al. | Accuracy of radiotherapy dose calculations based on cone-beam CT: comparison of deformable registration and image correction based methods | |
| Lu | A non-voxel-based broad-beam (NVBB) framework for IMRT treatment planning | |
| CN110246562B (en) | Method, device and computer system for determining sub-beam intensity in radiotherapy system | |
| Wu et al. | Improving proton dose calculation accuracy by using deep learning | |
| Pastor-Serrano et al. | Millisecond speed deep learning based proton dose calculation with Monte Carlo accuracy | |
| JP2015136625A (en) | Method for optimizing radiation dose for radiation therapy treatment and radiation therapy system | |
| CA2796159A1 (en) | Multi-objective radiation therapy optimization method | |
| Tsekas et al. | DeepDose: a robust deep learning-based dose engine for abdominal tumours in a 1.5 T MRI radiotherapy system | |
| Kusch et al. | A robust collision source method for rank adaptive dynamical low-rank approximation in radiation therapy | |
| Bergman et al. | Direct aperture optimization for IMRT using Monte Carlo generated beamlets | |
| Da Silva et al. | Sub-second pencil beam dose calculation on GPU for adaptive proton therapy | |
| US9251302B2 (en) | System and method for planning a radiation therapy treatment | |
| WO2012037472A2 (en) | Gpu-based fast dose calculator for cancer therapy | |
| Naqvi et al. | Convolution/superposition using the Monte Carlo method | |
| Li et al. | A new approach to integrate GPU-based Monte Carlo simulation into inverse treatment plan optimization for proton therapy | |
| Jhanwar et al. | Domain knowledge driven 3D dose prediction using moment-based loss function | |
| EP3445449A1 (en) | Fractionation selection tool in radiotherapy planning | |
| Vazquez et al. | A deep learning-based approach for statistical robustness evaluation in proton therapy treatment planning: a feasibility study | |
| Garrido-Hernandez et al. | Hypoxia adapted relative biological effectiveness models for proton therapy: a simulation study | |
| Li et al. | A new Monte Carlo-based treatment plan optimization approach for intensity modulated radiation therapy | |
| Tsekas et al. | Robust deep learning-based forward dose calculations for VMAT on the 1.5 T MR-linac |
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 |