CN107469239A - Computer-readable recording medium and radiotherapy treatment planning system - Google Patents
Computer-readable recording medium and radiotherapy treatment planning system Download PDFInfo
- Publication number
- CN107469239A CN107469239A CN201710686337.6A CN201710686337A CN107469239A CN 107469239 A CN107469239 A CN 107469239A CN 201710686337 A CN201710686337 A CN 201710686337A CN 107469239 A CN107469239 A CN 107469239A
- Authority
- CN
- China
- Prior art keywords
- organ
- discrete point
- point
- flux
- mark
- 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 21
- 210000000056 organ Anatomy 0.000 claims abstract description 171
- 230000004907 flux Effects 0.000 claims abstract description 116
- 238000013507 mapping Methods 0.000 claims abstract description 54
- 238000005457 optimization Methods 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 35
- 238000004590 computer program Methods 0.000 claims abstract description 9
- 230000015572 biosynthetic process Effects 0.000 claims abstract 2
- 230000006870 function Effects 0.000 claims description 8
- 238000009499 grossing Methods 0.000 claims description 3
- 230000001154 acute effect Effects 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000002721 intensity-modulated radiation therapy Methods 0.000 description 2
- 230000002285 radioactive effect Effects 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008365 organ synthesis Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 210000000952 spleen Anatomy 0.000 description 1
- 210000002784 stomach Anatomy 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
- A61N2005/1054—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using a portal imaging system
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N2005/1092—Details
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
The present invention relates to computer-readable recording medium, is stored thereon with computer program, can be used for performing a kind of flux pattern optimization method when the computer program is executed by processor, for producing the flux pattern in each launched field direction in radiotherapy treatment planning；Methods described includes：The discrete point of discretization is established in the plane perpendicular to the launched field direction；Each organ sampled point in launched field is mapped to the mapping point corresponding to formation in the plane in the launched field direction；The corresponding relation established between each mapping point and each discrete point；The distance of organ mark and telorism's irradiation source to each discrete point each mapping point according to corresponding to it determines initial value；Initial flux figure is formed according to the initial value of each discrete point；And the flux pattern after flux is optimized is carried out to the initial flux figure.Above-mentioned flux pattern optimization method can reduce the peak of acute variation and low ebb in the flux pattern after optimization.The present invention also provides a kind of radiotherapy treatment planning system.
Description
Technical field
The present invention relates to radiation therapy technology field, more particularly to a kind of computer-readable recording medium and radiotherapy
Planning system.
Background technology
Flux pattern refers to the distribution map for the flux (Fluence) that radioactive source obtains after multi-leaf optical grating.In radiotherapy
It is referred to as flux pattern optimization as intensity modulated radiation therapy is related to a kind of optimized algorithm in the works, its effect is to being needed to use in Rapid Dose Calculation
Flux pattern directly optimize, so as to be approached the plan for either meeting that physics teacher or doctor require, finally by near
Split to obtain executable Ziye shape and hop count like optimal flux pattern.Traditional flux pattern optimization method is to flux pattern
After optimizing, obtained flux pattern there may be the violent peak of change and low ebb.For such flux pattern, afterwards
Less opening can be produced in cutting procedure and is unfavorable for Rapid Dose Calculation, and more sons can be produced on the premise of precision is not reduced
Open country, the time is performed so as to add treatment.If be defined to Ziye number, during approximate segmentation is carried out to flux pattern
More errors can be produced.
The content of the invention
Based on this, it is necessary to provide the peak of acute variation and the meter of low ebb in a kind of flux pattern that can be reduced after optimizing
Calculation machine readable storage medium storing program for executing and radiotherapy treatment planning system.
A kind of computer-readable recording medium, is stored thereon with computer program, and the computer program is held by processor
It can be used for performing a kind of flux pattern optimization method during row；The flux pattern optimization method is used to produce in radiotherapy treatment planning respectively to penetrate
The flux pattern in wild direction；Methods described includes：
Discrete point is established in the plane perpendicular to the launched field direction；
By each organ sampled point in launched field be mapped to it is described formed in the plane in the launched field direction corresponding to reflect
Exit point；
The corresponding relation established between each mapping point and each discrete point；
The distance of organ mark and telorism's irradiation source to each discrete point each mapping point according to corresponding to it is true
Determine initial value；
Initial flux figure is formed according to the initial value of each discrete point；And
The flux pattern after flux is optimized is carried out to the initial flux figure.
Above computer readable storage medium storing program for executing, the flux pattern optimization method of its memory storage are right when forming initial flux figure
Each discrete point organ mark of each mapping point and the distance of telorism's irradiation source can determine initial value according to corresponding to it, from
And the initial Distribution value of rational flux can be provided.Pass through the initial flux figure that is initially worth to relatively ideal flux
Figure, therefore during subsequently flux optimization is carried out to flux pattern, can only less adjustment be carried out to initial flux figure, so as to
Reduce the peak of acute variation and low ebb in the flux pattern after optimization indirectly.
In one of the embodiments, the step of establishing discrete point in the plane perpendicular to the launched field direction is wrapped
Include：The grid of discretization is established in the plane perpendicular to the launched field direction；Described in grid node on the grid is used as
Discrete point.
In one of the embodiments, it is described establish between each mapping point and each discrete point corresponding relation the step of be,
Each mapping point is corresponded on discrete point according to the distance between each mapping point and each discrete point.
In one of the embodiments, the organ mark to each discrete point each mapping point according to corresponding to it and described
The step of distance of telorism's irradiation source determines initial value includes, to each discrete point：
The organ mark of each mapping point and the distance of telorism's irradiation source determine distance according to corresponding to discrete point
The nearest organ mark of the irradiation source；And
The initial value of the discrete point is determined according to the organ mark nearest apart from the irradiation source.
In one of the embodiments, it is described that the discrete point is determined according to the organ mark nearest apart from the irradiation source
Initial value the step of in：
When the organ nearest apart from the irradiation source is identified as target area organ mark, the initial value of the discrete point according to
The target dose of target area organ determines；
When the organ nearest apart from the irradiation source, which is identified as, jeopardizes organ mark, the initial value of the discrete point according to
The upper dosage limit for jeopardizing organ determines with weight.
In one of the embodiments, the organ mark to each discrete point each mapping point according to corresponding to it and described
The step of distance of telorism's irradiation source determines initial value includes, to each discrete point：
The organ mark of each mapping point and the distance of telorism's irradiation source determine distance according to corresponding to discrete point
The nearest organ mark of the irradiation source；
When the organ nearest apart from the irradiation source is identified as target area organ mark：
In the presence of when jeopardizing organ in the range of the pre-determined distance of target area organ, jeopardize organ and the target area device according to described
Official determines the initial value of the discrete point；
Be not present in the range of the pre-determined distance of target area organ when jeopardizing organ, according to determining the target area organ from
The initial value of scatterplot；
When the organ nearest apart from the irradiation source, which is identified as, jeopardizes organ mark：
When the organ mark for target area organ mark being not present in the discrete point, institute is determined according to the organ that jeopardizes
State the initial value of discrete point；
When the organ mark of target area organ mark is existed in the discrete point, according to the target area organ and the danger
And organ determines the initial value of the discrete point.
In one of the embodiments, it is described establish between each mapping point and each discrete point corresponding relation the step of after
Also include：Each discrete point is labeled；
Described the step of being labeled to each discrete point for by the organ of mapping point corresponding to each discrete point mark and it is described
The distance of telorism's irradiation source is labeled on discrete point；Or will be nearest apart from irradiation source in mapping point corresponding to each discrete point
The organ of an organ be labeled on discrete point.
In one of the embodiments, the step of methods described also includes being smoothed the flux pattern after optimization.
In one of the embodiments, the step of flux pattern after described pair of optimization is smoothed includes, according to pre-
If smooth optimization model is smoothed to the flux pattern after optimization；The object function of the default smooth optimization model includes
Variations of flux penalty and smooth penalty；The variations of flux penalty is used to provide flux new on each discrete point
The punishment that the variable quantity of value and former amount of flux is brought；The smooth penalty is used to provide amount of flux on adjacent multiple discrete points
Difference band come punishment.
A kind of radiotherapy treatment planning system, including memory and processor；The memory includes such as foregoing any implementation
Computer-readable recording medium described in example.
Brief description of the drawings
Fig. 1 is the flow chart of the flux pattern optimization method in an embodiment；
Fig. 2 is that each mapping point is corresponded to after discrete point to the macroscopical schematic diagram formed in an embodiment；
Fig. 3 is the partial enlarged drawing in Fig. 2；
Fig. 4 is the particular flow sheet of step S140 in an embodiment；
Fig. 5 is the particular flow sheet of step S140 in another embodiment.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, it is right below in conjunction with drawings and Examples
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
Radiotherapy treatment planning system in one embodiment, the pattern that can be used in its radiotherapy treatment planning, which includes adjusting, puts by force
Penetrate treatment (Intensity Modulated Radiation Therapy, MRT), rotation IMRT (Intensity
Modulated Arc Radiation Therapy, IMAT), Tomotherapy radiotherapy (Tomotherapy) etc. and it is mixed
Syntype.The radiotherapy treatment planning system includes memory and processor.Memory includes computer-readable recording medium.Meter
Computer program is stored with calculation machine readable storage medium storing program for executing, the computer program can be used for performing one kind when being executed by processor
Flux pattern optimization method.The flux pattern optimization method, which is used to produce, obtains the flux pattern in each launched field direction in radiotherapy treatment planning.
Caused flux pattern can be as the reference information of follow-up radiotheraphy device in above-mentioned radiotherapy treatment planning system.In an embodiment
In, radiotherapy treatment planning system and radiotheraphy device can integrate in a system or two independent systems.
Flux pattern optimization method can be used for the flux pattern for obtaining each launched field direction.In rotation IMRT
(Arc) in the works, if can turn to all arcs are discrete per one launched field direction of mass dryness fraction, by frame (gantry) angle, leafy
A launched field direction is merged into the attribute identical launched field such as grating (Collimator) angle direction, then to each launched field direction
Handled using the flux pattern optimization method in inventive embodiments.
Referring to Fig. 1, the flux pattern optimization method in an embodiment comprises the following steps：
Step S110, the discrete point of discretization is established in the plane perpendicular to launched field direction.
Beam axle visual angle (beam ' s eye veiw, abbreviation BEV) is to contemplate station in source positions, along beam centre axle side
To observation launched field.BEV planes are the planes vertical with launched field direction.If represent launched field side with the direction of beam central shaft
To then BEV planes are perpendicular to the plane in launched field direction, its plane substantially vertical with beam centre axle.In one embodiment,
BEV planes can be arranged at radioactive source and wait between central plane.In one embodiment, BEV planes can be arranged at leafy light
Plane where grid.The discrete point of discretization can arrange in any way, such as be rendered as the point of rule arrangement or random row
The point of cloth, rule arrangement can be uniformly arrangement, or uneven arrangement.
Step S120, each organ sampled point in launched field direction is mapped to this and formed in the plane in launched field direction pair
The mapping point answered.
Each organ sampled point in launched field direction can be mapped to perpendicular to launched field direction along the opposite direction of beam direction
In plane.Organ sampled point can carry out selection setting as needed, such as each organ is sampled using the method for sampling, from
And obtain each organ sampled point.Each organ sampled point can be to be uniformly distributed or the point of random distribution.Not to adopting in the application
Quadrat method is defined, and the method for sampling corresponding to each organ be able to can also be differed with identical, also without limit in the application
It is fixed.
Step S130, the corresponding relation established between each mapping point and each discrete point.
Each mapping point is corresponded on corresponding discrete point according to preset rules.
Step S140, to the organ of each discrete point each mapping point according to corresponding to it mark and telorism's irradiation source
Distance determines initial value.
Phase is determined according to the distance of the organ mark of corresponding each mapping point on each discrete point and telorism's irradiation source
The field size dosage of adaptation, so as to set the relatively reasonable initial Distribution value of flux for each discrete point.The dosage of sampled point with it is right
The amount of flux for the discrete point answered is the relation of direct ratio, and the dosage of sampled point is bigger, and the amount of flux of discrete point is bigger.Discrete point leads to
Value size can be several times or the part for the corresponding sampling dose point being calculated.
Step S150, initial flux figure is formed according to the initial value of each discrete point.
Step S160, the flux pattern after being optimized is optimized to initial flux figure.
The known Optimized model of those skilled in the art can be used to carry out flux optimization processing to initial flux figure, obtained
Flux pattern after optimization.
Above-mentioned flux pattern optimization method, when forming initial flux figure, each discrete point can respectively be mapped according to corresponding to it
The organ mark of point and the distance of telorism's irradiation source determine initial value, so as to provide rational flux initial value point
Cloth.For the mode that the initial value of traditional given each discrete point is zero, the initial value in above-mentioned optimization method is
Close to optimal dosage requirement, so that obtained initial flux figure relatively ideal flux figure.Therefore follow-up right
During flux pattern carries out flux optimization, can only less adjustment be carried out to initial flux figure, so as to reduce optimization indirectly
The peak of acute variation and low ebb in flux pattern afterwards.Therefore less opening will not be produced in cutting procedure afterwards, favorably
In Rapid Dose Calculation.Also, on the premise of precision is not reduced, it is possible to reduce the number of caused Ziye, so as to shorten treatment
Perform the time.If be defined to Ziye number, caused error is smaller during approximate segmentation is carried out to flux pattern, has
Beneficial to raising precision.
In one embodiment, step S110 includes：The step of the grid of discretization is established in the plane perpendicular to launched field direction
Suddenly.Grid node on grid is as discrete point.Grid can be uniform grid, non-uniform grid or Random Discrete grid
Deng so that its grid node shows uniform, non-homogeneous or stochastic regime.In one embodiment, perpendicular to launched field side
To plane namely BEV planes can be arranged at plane where multi-leaf optical grating.Now, with the blade bearing of trend of multi-leaf optical grating
For X-direction, using width of blade direction as Y direction.The grid interval of the grid of foundation in the X-axis direction sets for radiotherapy
The minimum blade of multi-leaf optical grating (MLC) in standby closes default times of interval or sampling interval, such as 1.1 times.Sampling interval is
Refer to the sampling interval sampled to each organ.The grid spacing of grid in the Y-axis direction is the width of blade of multi-leaf optical grating.
In other embodiments, interval of the grid in X-direction and Y direction can be arranged as required to.
In one embodiment, it is by each mapping point pair in step S130 according to the distance between each mapping point and each discrete point
It should arrive on discrete point (or grid node).For example, each mapping point is corresponded on neighbouring flux point according to nearby principle.Tool
Body, each mapping point can be corresponded on four neighbouring discrete points according to nearby principle, so that it is guaranteed that what is finally given is logical
Spirogram is as close possible to ideal flux figure.In the present embodiment, the distance between four discrete points and mapping point can not be entered
Row is limited, it is only necessary to which mapping point is corresponded into default 4 neighbouring discrete points according to distance.In other embodiments
In, each mapping point can also be corresponded to neighbouring one, on two or more discrete points according to nearby principle.Corresponding to
Cheng Zhong, the distance between mapping point and discrete point can be defined, such as mapping point is corresponded in the range of pre-determined distance
All discrete points on.In another embodiment, mapping point can also be corresponded on a nearest discrete point.
In one embodiment, after step s 130, the step of in addition to each discrete point being labeled.Specifically, will
The organ mark of mapping point corresponding to each discrete point and the distance of telorism's irradiation source are labeled on discrete point.Organ mark
Know the type for identifying the organ.Organ mark can directly be words identification or corresponding symbol, pattern identification etc..
In one embodiment, in labeling process, using the distance of each telorism's irradiation source corresponding to each mapping point as its Z axis coordinate.
Therefore, amount of flux can be determined according to the organ mark on each mapping point and its corresponding Z axis coordinate in step S140
Initial value.
In one embodiment, the schematic diagram after mark can be exported to display interface and checked to operator., can be with during output
The organ mark being marked only is shown, and the distance of telorism's irradiation source can't be shown, as shown in Figure 2.Fig. 2 is one
Each mapping point is corresponded to after discrete point to the macroscopical schematic diagram formed in embodiment.From figure 2 it can be seen that target detection position
Including target area organ (PTV) and jeopardize organ (OAR).Fig. 3 is in Fig. 2 20 partial enlarged drawing.In figure 3, solid line circle is
Target area organ mark, and dashed circle then identifies to jeopardize organ, triangle is discrete point., can be by respectively according to nearby principle
Mapping point is corresponded on four neighbouring discrete points, therefore target area organ mark has been only existed on some discrete point, such as a left side in Fig. 3
Four discrete points of side；Target area organ mark on some discrete point being present, there is also jeopardize organ mark such as upper right two in Fig. 3
Individual discrete point；And some discrete point then only exists and jeopardizes organ mark such as the discrete point of bottom right two in Fig. 3.Implement in others
In example, only the organ of an organ nearest apart from irradiation source in mapping point corresponding to each discrete point can also be labeled in
On discrete point.
In one embodiment, in step S140, its is initial can be determined according to the flow shown in Fig. 4 to each discrete point
Value, is specifically included：
Step S210, the organ mark of each mapping point and the distance of telorism's irradiation source determine according to corresponding to discrete point
The organ mark nearest apart from irradiation source.
Compare telorism's irradiation source in each mapping point corresponding to discrete point distance determine it is nearest apart from irradiation source
Organ identifies, and influences maximum organ so as to which organ corresponding to organ mark is defined as into beam.
Step S220, the initial value of discrete point is determined according to the organ mark nearest apart from irradiation source.
Specifically, when the organ determined is identified as target area organ mark, the initial value of discrete point is according to target area device
The target dose of official determines.Target dose can be an occurrence, or one has the numerical value model of higher limit and lower limit
Enclose, be set by the user.The ratio between such as the initial value of discrete point is equal to target dose and the launched field number of target area organ.
Herein, launched field number is the summation of the launched field quantity on all launched field directions.
When the organ determined, which is identified as, jeopardizes organ mark, the initial value of discrete point jeopardizes organ according to
Upper dosage limit determines with weight.Upper dosage limit represents this and jeopardizes the upper limit that organ can bear dosage, is set by the user；Weight generation
Table corresponds to the relative importance of organ.For example, the initial value of discrete point is set smaller than or jeopardizes the upper of organ equal to this
Limit the ratio between dosage and launched field number.By by the initial value of discrete point be set smaller than or equal to jeopardize the upper dosage limit of organ with
The ratio between launched field number, it can be protected in the allowed band for jeopardize organ.
In one embodiment, in step S140, its is initial can be determined according to the flow shown in Fig. 5 to each discrete point
Value, is specifically included：
Step S310, the organ mark of each mapping point and the distance of telorism's irradiation source determine according to corresponding to discrete point
The organ mark nearest apart from irradiation source.
When determining that the organ nearest apart from irradiation source is identified as target area organ mark, step S320 is performed；When true
Fix the organ nearest apart from irradiation source to be identified as when jeopardizing organ mark, perform step S330.
Step S320, judge to whether there is in the range of the pre-determined distance of target area organ and jeopardize organ.
Pre-determined distance scope can be set as needed, as determined that this is default according to the significance level for jeopardizing organ
Distance range etc..Exist in the range of the pre-determined distance of target area organ when jeopardizing organ, perform step S340, otherwise perform step
S350。
Step S340, according to jeopardizing organ and target area organ determines the initial value of discrete point.
Specifically, according to the initial value for jeopardizing organ and target area organ synthesis determination discrete point.In one embodiment, can be with
The significance level for respectively jeopardizing organ is defined.In one embodiment, the weight of corresponding organ can be determined according to organ mark
Degree is wanted, such as organ mark can be with liver, spleen, stomach, backbone etc., the significance level of these organs is height.Jeopardizing organ
Significance level it is higher when, the initial value of discrete point is set smaller than or jeopardizes the upper dosage limit of organ and launched field equal to this
The ratio between number, protected so as to jeopardize organ to this., then can be by the first of discrete point when the important procedure for jeopardizing organ is relatively low
Initial value is equal to the ratio between target dose and launched field number of target area organ.In other examples, when jeopardizing organ
, can also be by the protection that suitably reduces the initial value of discrete point to realize to jeopardizing organ when significance level is relatively low.Jeopardize device
The height of the significance level of official can predefine, so as to get jeopardize organ mark when i.e. can be known this jeopardize the weight of organ
Want the height situation of degree.
Step S350, the initial value of discrete point is determined according to target area organ.
The step of determining the initial value of discrete point according to target area organ by the agency of in preceding embodiment, does not go to live in the household of one's in-laws on getting married herein
State.
Step S330, judge to whether there is the organ mark for target area organ mark in the discrete point.
When existing for the organ mark of target area organ mark, step S360 is performed, otherwise perform step S370.
Step S360, according to target area organ and jeopardize the initial value that organ determines discrete point.
When existing for the organ mark of target area organ mark, the initial value of discrete point is equal to this and jeopardizes organ
Upper dosage limit and the ratio between launched field number.By the way that the initial value of discrete point to be equal to jeopardize to the upper dosage limit of organ with penetrating
The ratio between wild number, it can both be protected in the allowed band for jeopardize organ, target area organ can also be radiated
Treatment.In another embodiment, it is also contemplated that jeopardize the upper dosage limit of organ and weight determines the initial value of discrete point.
Step S370, according to the initial value for jeopardizing organ and determining discrete point.
When being identified in the absence of the organ identified for target area organ, the initial value of discrete point is set smaller than or equal to this
Jeopardize the ratio between upper dosage limit and launched field number of organ.In one embodiment, can be according to the upper dosage limit and power for jeopardizing organ
The initial value of discrete point is determined again.In one embodiment, the initial value of the discrete point can be arranged to zero, so as to avoid jeopardizing
Organ is by unnecessary irradiation.
By the above method, consider that beam influences maximum organ type in selected initial value, it is rational so as to provide
The initial Distribution value of flux, final obtained initial flux figure have been closer to preferable flux pattern, therefore subsequently to first
, only can be logical after optimization so as to reduce indirectly to the less adjustment of initial flux figure during beginning flux pattern optimizes
The peak of acute variation and low ebb in spirogram.Alternatively, it is also possible to accelerate the decrease speed of optimization object function, Optimized Iterative is reduced
Number, shorten the optimization time.
In one embodiment, above-mentioned flux pattern optimization method also includes step S170, and the flux pattern after optimization is put down
The step of sliding processing, as shown in Figure 1.It is in the present embodiment, independent by the way that the flux of smoothing processing and flux pattern is optimized,
The problem of dimension disunity can be avoided to cause effect of optimization difference, occurs.
In one embodiment, the flux pattern after optimization is smoothed using default smooth optimization model.It is default flat
The object function of sliding Optimized model includes variations of flux penalty and smooth penalty.Variations of flux penalty and smooth
Penalty sum forms the object function.Wherein, variations of flux penalty is used to provide flux new on each discrete point
The punishment that the variable quantity of value and former amount of flux is brought；Smooth penalty is then used to provide amount of flux on adjacent multiple discrete points
The punishment that difference band comes.
In one embodiment, the object function for presetting smooth optimization model is as follows：
Wherein, δ () is variations of flux penalty, to provide amount of flux f new on i-th of grid nodeiWith original
Amount of fluxThe punishment that brings of variable quantity, such as δ (x)=(x)2.K is change penalty coefficient；For smooth penalty,
To represent amount of flux on adjacent multiple discrete pointsDifference band come punishment, such as(a1, a2... ..., an) ∈ N represent the labels of adjacent multiple discrete points, N is the set of all adjacent discrete points.
In the present embodiment, the common optimum method of quadratic function can be used to optimize.For example, musical instruments used in a Buddhist or Taoist mass is solved to target using BFGS
Function carries out successive ignition optimization, and iterations can be determined according to the precision for finally needing to obtain, such as to target
Function carries out 20 iteration optimizations.It is appreciated that it can also be optimized by the way of other need not be iterated.
In traditional flux pattern optimization process, it can be carried out being added to the smoothing processing of flux pattern in flux optimization together.
Because variations of flux punishment is different with the dimension smoothly punished, effect of optimization is substantially influenceed by both coefficients.And use this
Flux pattern optimization method in embodiment, flux pattern is smoothly punished independently of flux optimization process, so as to the effectively amount of avoiding
The problem of guiding principle difference causes effect of optimization impacted.
One embodiment of the invention also provides a kind of computer-readable recording medium, is stored thereon with computer program.The meter
When calculation machine program is executed by processor, it can be used for performing the flux pattern optimization method as described in foregoing any one embodiment.
Each technical characteristic of embodiment described above can be combined arbitrarily, to make description succinct, not to above-mentioned reality
Apply all possible combination of each technical characteristic in example to be all described, as long as however, the combination of these technical characteristics is not deposited
In contradiction, the scope that this specification is recorded all is considered to be.
Embodiment described above only expresses the several embodiments of the present invention, and its description is more specific and detailed, but simultaneously
Can not therefore it be construed as limiting the scope of the patent.It should be pointed out that come for one of ordinary skill in the art
Say, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to the protection of the present invention
Scope.Therefore, the protection domain of patent of the present invention should be determined by the appended claims.
Claims (10)
1. a kind of computer-readable recording medium, is stored thereon with computer program, it is characterised in that the computer program quilt
It can be used for performing a kind of flux pattern optimization method during computing device；The flux pattern optimization method is based on producing radiotherapy
The flux pattern in each launched field direction in drawing；Methods described includes：
Discrete point is established in the plane perpendicular to the launched field direction；
Each organ sampled point in launched field is mapped to the mapping point corresponding to formation in the plane in the launched field direction；
The corresponding relation established between each mapping point and each discrete point；
The distance of organ mark and telorism's irradiation source to each discrete point each mapping point according to corresponding to it determines just
Initial value；
Initial flux figure is formed according to the initial value of each discrete point；And
The flux pattern after flux is optimized is carried out to the initial flux figure.
2. computer-readable recording medium according to claim 1, it is characterised in that described perpendicular to the launched field side
To plane on establish discrete point the step of include：The grid of discretization is established in the plane perpendicular to the launched field direction；
Grid node on the grid is as the discrete point.
3. computer-readable recording medium according to claim 1, it is characterised in that it is described establish each mapping point with respectively from
The step of corresponding relation between scatterplot is, according to the distance between each mapping point and each discrete point by each mapping point correspond to from
In scatterplot.
4. computer-readable recording medium according to claim 1, it is characterised in that described right according to its to each discrete point
The step of organ mark for each mapping point answered and the distance of telorism's irradiation source determine initial value includes, to each discrete
Point：
The organ mark of each mapping point and the distance of telorism's irradiation source are determined described in distance according to corresponding to discrete point
The nearest organ mark of irradiation source；And
The initial value of the discrete point is determined according to the organ mark nearest apart from the irradiation source.
5. computer-readable recording medium according to claim 4, it is characterised in that described according to apart from the irradiation source
In the step of nearest organ mark determines the initial value of the discrete point：
When the organ nearest apart from the irradiation source is identified as target area organ mark, the initial value of the discrete point is according to target area
The target dose of organ determines；
When the organ nearest apart from the irradiation source, which is identified as, jeopardizes organ mark, the initial value of the discrete point is according to
The upper dosage limit and weight for jeopardizing organ determine.
6. computer-readable recording medium according to claim 1, it is characterised in that described right according to its to each discrete point
The step of organ mark for each mapping point answered and the distance of telorism's irradiation source determine initial value includes, to each discrete
Point：
The organ mark of each mapping point and the distance of telorism's irradiation source are determined described in distance according to corresponding to discrete point
The nearest organ mark of irradiation source；
When the organ nearest apart from the irradiation source is identified as target area organ mark：
Exist in the range of the pre-determined distance of target area organ when jeopardizing organ, jeopardize organ and the target area organ is true according to described
The initial value of the fixed discrete point；
In the absence of when jeopardizing organ in the range of the pre-determined distance of target area organ, the discrete point is determined according to the target area organ
Initial value；
When the organ nearest apart from the irradiation source, which is identified as, jeopardizes organ mark：
In the discrete point be not present for target area organ mark organ mark when, according to it is described jeopardize organ determine described in from
The initial value of scatterplot；
When the organ mark of target area organ mark is existed in the discrete point, according to the target area organ and described jeopardize device
Official determines the initial value of the discrete point.
7. computer-readable recording medium according to claim 1, it is characterised in that it is described establish each mapping point with respectively from
Also include after the step of corresponding relation between scatterplot：Each discrete point is labeled；
Described the step of being labeled to each discrete point for by the organ of mapping point corresponding to each discrete point mark and the organ
Distance apart from irradiation source is labeled on discrete point；Or by mapping point corresponding to each discrete point apart from irradiation source it is nearest one
The organ of individual organ is labeled on discrete point.
8. computer-readable recording medium according to claim 1, it is characterised in that after methods described also includes to optimization
Flux pattern the step of being smoothed.
9. computer-readable recording medium according to claim 8, it is characterised in that the flux pattern after described pair of optimization enters
The step of row smoothing processing, includes, and the flux pattern after optimization is smoothed according to default smooth optimization model；It is described pre-
If the object function of smooth optimization model includes variations of flux penalty and smooth penalty；The variations of flux punishes letter
The punishment that number is brought for providing the variable quantity of amount of flux new on each discrete point and former amount of flux；The smooth penalty is used
In the punishment that the difference band for providing amount of flux on adjacent multiple discrete points comes.
10. a kind of radiotherapy treatment planning system, including memory and processor；Characterized in that, the memory is included as weighed
Profit requires the computer-readable recording medium described in 1~9 any one.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710686337.6A CN107469239B (en) | 2017-08-11 | 2017-08-11 | Computer-readable storage medium and radiation therapy planning system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710686337.6A CN107469239B (en) | 2017-08-11 | 2017-08-11 | Computer-readable storage medium and radiation therapy planning system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107469239A true CN107469239A (en) | 2017-12-15 |
| CN107469239B CN107469239B (en) | 2020-07-14 |
Family
ID=60600292
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710686337.6A Active CN107469239B (en) | 2017-08-11 | 2017-08-11 | Computer-readable storage medium and radiation therapy planning system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107469239B (en) |
Cited By (2)
| 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 |
| CN113577581A (en) * | 2021-08-30 | 2021-11-02 | 上海联影医疗科技股份有限公司 | Radiation therapy dose determination system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103038669A (en) * | 2010-06-22 | 2013-04-10 | 卡尔·奥托 | System and method for estimating and manipulating estimated radiation dose |
| CN105413068A (en) * | 2015-12-31 | 2016-03-23 | 上海联影医疗科技有限公司 | Optimization method for solving flux map and radiotherapy equipment |
| EP3075416A1 (en) * | 2015-04-02 | 2016-10-05 | RaySearch Laboratories AB | System and method for determining a radiation treatment plan and a radiation treatment machine |
| US20170021192A1 (en) * | 2011-09-28 | 2017-01-26 | Varian Medical Systems, Inc. | Radiation therapy treatment plan improvement through use of a knowledge base |
-
2017
- 2017-08-11 CN CN201710686337.6A patent/CN107469239B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103038669A (en) * | 2010-06-22 | 2013-04-10 | 卡尔·奥托 | System and method for estimating and manipulating estimated radiation dose |
| US20170021192A1 (en) * | 2011-09-28 | 2017-01-26 | Varian Medical Systems, Inc. | Radiation therapy treatment plan improvement through use of a knowledge base |
| EP3075416A1 (en) * | 2015-04-02 | 2016-10-05 | RaySearch Laboratories AB | System and method for determining a radiation treatment plan and a radiation treatment machine |
| CN105413068A (en) * | 2015-12-31 | 2016-03-23 | 上海联影医疗科技有限公司 | Optimization method for solving flux map and radiotherapy equipment |
Cited By (4)
| 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 |
| CN113577581A (en) * | 2021-08-30 | 2021-11-02 | 上海联影医疗科技股份有限公司 | Radiation therapy dose determination system |
| CN113577581B (en) * | 2021-08-30 | 2024-02-20 | 上海联影医疗科技股份有限公司 | Radiotherapy dosage determination system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107469239B (en) | 2020-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11986671B2 (en) | System and method for estimating and manipulating estimated radiation dose | |
| US10076673B2 (en) | Interactive dose gradient based optimization technique to control IMRT delivery complexity | |
| CN106730407A (en) | A kind of scanning illuminating method for particle therapy, device and treatment head | |
| Yue et al. | Dose prediction via distance-guided deep learning: Initial development for nasopharyngeal carcinoma radiotherapy | |
| Hwang et al. | A novel model to correlate hydrogel spacer placement, perirectal space creation, and rectum dosimetry in prostate stereotactic body radiotherapy | |
| CN107469239A (en) | Computer-readable recording medium and radiotherapy treatment planning system | |
| JP5401391B2 (en) | Particle beam therapy planning apparatus and therapy planning method | |
| CN115862812A (en) | Radiotherapy plan processing device, storage medium and electronic device | |
| Yuan et al. | A dose distribution evaluation technique using the k‐d tree for nearest neighbor searching | |
| CN105561485A (en) | Radiotherapy treatment planning optimization method and device | |
| van der Meer et al. | Bi‐objective optimization of catheter positions for high‐dose‐rate prostate brachytherapy | |
| CN114420252A (en) | Method, device and medium for determining intensity modulated radiotherapy plan evaluation parameter matrix | |
| CN111494815B (en) | Three-dimensional dose calculation method, device and medium based on mixed variable-scale model | |
| Hirayama et al. | Difference in LET‐based biological doses between IMPT optimization techniques: Robust and PTV‐based optimizations | |
| WO2022088341A1 (en) | Radiation therapy control method and apparatus | |
| CN112635023B (en) | Method for generating dose prediction model of nasopharyngeal carcinoma, dose prediction method and device | |
| Böck | On adaptation cost and tractability in robust adaptive radiation therapy optimization | |
| Gaddy et al. | Improving the VMERGE treatment planning algorithm for rotational radiotherapy | |
| Tian et al. | Moving GPU‐Open CL‐based Monte Carlo dose calculation toward clinical use: Automatic beam commissioning and source sampling for treatment plan dose calculation | |
| Hild et al. | Fast optimization and dose calculation in scanned ion beam therapy | |
| US20240091559A1 (en) | Enhanced dose rate/ultra-high dose rate radiation and spatially fractioned radiation therapy | |
| EP4427803A1 (en) | Radiation treatment planning systems and methods with interpolable spot placements and dose distributions | |
| CN117695533A (en) | Ion beam arc intensity modulated radiotherapy energy selection method and system | |
| CN117138260A (en) | Treatment planning system, method, equipment and medium for arc treatment of tumor by ion beam | |
| CN118593914A (en) | Treatment planning method and system for controlling dose distribution uniformity of a radiation treatment field |
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 | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |