CN109513118A - A kind of the photon energy synthetic method and system of clinac - Google Patents
A kind of the photon energy synthetic method and system of clinac Download PDFInfo
- Publication number
- CN109513118A CN109513118A CN201811312180.1A CN201811312180A CN109513118A CN 109513118 A CN109513118 A CN 109513118A CN 201811312180 A CN201811312180 A CN 201811312180A CN 109513118 A CN109513118 A CN 109513118A
- Authority
- CN
- China
- Prior art keywords
- photon
- energy
- dosimetric data
- data
- dosimetric
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1042—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
- A61N5/1031—Treatment planning systems using a specific method of dose optimization
-
- 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/1064—Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring
- A61N5/1065—Beam adjustment
-
- 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/1077—Beam delivery systems
- A61N5/1084—Beam delivery systems for delivering multiple intersecting beams at the same time, e.g. gamma knives
Abstract
The present invention relates to the photon energy synthetic methods and system of clinac, wherein, the photon energy synthetic method of clinac, comprising the following steps: obtain the first photon energy beam dosimetric data, the second photon energy beam dosimetric data and intermediate energy photon beam dosimetric data;Mathematical Fitting is carried out to the first photon energy beam dosimetric data and the second photon energy beam dosimetric data according to preset fitting coefficient, obtains synthesis photon energy beam dosimetric data;Verifying is compared with intermediate energy photon beam dosimetric data to synthesis photon energy beam dosimetric data, obtains dosimetric data fitness bias;Dosimetric data fitness bias is compared with preset threshold, obtains comparison result.It in the inventive solutions, being capable of the photon beam synthesis using two grades of energy of height and photon beam similar in photon energy beam beam behaviour any in both energy ranges.
Description
Technical field
The present invention relates to radiation therapy technology field more particularly to a kind of photon energy synthesis sides of clinac
The photon energy synthesis system of method and a kind of clinac.
Background technique
Clinac is a kind of radiotherapy equipment, is divided into electric wire and two kinds of photon line.One simplest
Clinac at least will include accelerating tube, high-power microwave source, Wave guide system, control system, ray expansion system and
The systems such as guard system.Clinac is the product that the high, precision and frontier technology in medical instrument combines, and is Medical treatment device
One of highest product of technology content in tool field also only has a few developed country to be able to produce in the world.
Conventional medical linear accelerator is divided into two class of traveling wave and standing wave, basic structure class according to the characteristics of microwave transmission
Seemingly, it generates high energy electron or photon beam mechanism is same.By utilizing the high energy electron and HIGH-POWERED MICROWAVES with certain energy
Microwave electromagnetic field interactions, to obtain higher energy electronic beam.Electron beam is directly led out, and can make electric wire treatment.
High energy electron hits heavy metal tungsten target, to generate bremstrahlen, emits sigmatron, resulting photon has continuous
Energy spectrum, energy range is from 0 to initiating electron beam energy.It is arrived for example, a 6MV electron beam can produce an energy 0
Photon beam between 6MV, this photon beam are commonly referred to as 6MV photon beam.
The penetration depth of different-energy ray is different, to the difference of tumor dose contribution, compared with low energy photon beam (≤
6MV) due to its percentage depth dose distribution within the organization, it is often used to treatment head and neck neoplasm or superficial tumor;It is higher
Stronger, the more tumour that be used to treat deep of the penetration capacity of the photon beam (> 6MV) of energy within the organization.High-energy photon
Line (> 6MV) has the characteristics that skin injury is small, penetrability is strong, the ray uniformity is high and guarantees that normal tissue effect is good, is suitable for
Treat deep tumor.Therefore it is directed to different parts difference tumour, the ray of different-energy may be selected when formulating treatment plan scheme
Carry out radiotherapy.
But the requirement of complex process and output dose rate due to accelerator engineering, the design of accelerator only allow a small amount of
Photon energy (usually 2 to 3 grades) is available.As the selectable energy gear number of clinac output light beamlet increases,
Technology complexity increases, so that manufacturing cost sharply increases, this also translates into higher purchasing price.Mainly pass through at present
Hardware controlling method adjusts clinac energy, including following two method: first is that by adjusting modulator high pressure,
Change magnetron or klystron output power, multiple kinds of energy selection is realized, the disadvantage is that energy-conditioning range is small;Second is that using energy
Method of switching carries out multiple kinds of energy selection, the disadvantage is that only there are two gears for the selectable energy of accelerator, can not continuously adjust.
To ensure that dosage is accurately and reliably irradiated by prescription, every grade of photon energy beam must be calibrated and be tested respectively,
And every grade of photon energy beam must individually carry out modeling test in treatment planning systems.Debugging need to be using three-dimensional water when checking and accepting
Case measures the dosimetric parameters under various geometrical conditions in water, such as percent depth-dose curves (Percentage Depth
Dose, PDD), off-axis dose profile (Off Axial Ratio, OAR) and the output factor (Output Factor, OF)
Deng importing treatment planning systems and generate dose computation model appropriate.Physics teacher is when planning design from treatment planning systems
Wherein one grade of available energy photon carries out plan design for middle selection.Therefore, the cost debugged and safeguarded is usually and medical linear
The quantity for the photon energy that accelerator can be used is directly proportional.It, will significantly when the quantity for the photon energy that can be used increases
Increase the operating cost of clinac.
Summary of the invention
The present invention is directed to solve at least one of the technical problems existing in the prior art or related technologies.
For this purpose, it is an object of the present invention to provide a kind of photon energy synthetic method of clinac,
It can be closed using photon beams characteristics such as the percent depth-dose curves of the photon of two grades of known energies and off-axis dose profiles
At with photon beam similar in the beam behaviour of energy any in both energy ranges.
It is another object of the present invention to provide a kind of photon energy synthesis systems of clinac, can
The step of realizing the photon energy synthesis of clinac involved in the above method.
To achieve the above object, the technical solution of first aspect present invention provides a kind of photon of clinac
Energy synthetic method, comprising the following steps: obtain the first photon energy beam dosimetric data, the second energy photon beam dosimetry number
Accordingly and corresponding to the first photon energy beam dosimetric data and the second photon energy beam dosimetric data intermediate energy photon
Beam dosimetry data;According to preset fitting coefficient to the first photon energy beam dosimetric data and the second energy photon beam dose
It learns data and carries out Mathematical Fitting, obtain synthesis photon energy beam dosimetric data;To synthesis photon energy beam dosimetric data with
Intermediate energy photon beam dosimetric data is compared verifying, obtains dosimetric data fitness bias;Dosimetric data is fitted
Deviation is compared with preset threshold, obtains comparison result;When comparison result is that dosimetric data fitness bias is less than or equal to
When preset threshold, synthesis photon energy beam dosimetric data is final dose data;When comparison result is quasi- for dosimetric data
When closing deviation greater than preset threshold, it is iterated according to dosimetric data fitness bias and optimizes and revises fitting coefficient.
In the technical scheme, the ray of the high and low two grades of energy exported from existing clinac is according to upper
The fitting coefficient that face determines mixes, and penetrates similar in the beam behaviour of any energy in the above two energy range of simulated implementation
Line, so that realizing can be identical or another medical with different photon energies by the radiation treatment plan of known photon beam
It realizes on linear accelerator, is distributed with reaching close to identical dosage;Meanwhile outfit two is not needed on clinac
A above photon energy can save the spending and later standing charges of equipment research and development;May further be used to offer arbitrary continuation can
Photon energy is adjusted, provides maximum freedom degree for user, this is not accomplished on current hardware.
In the above-mentioned technical solutions, it is preferable that the first photon energy beam dosimetric data, the second energy photon beam dosimetry
Data and intermediate energy photon beam dosimetric data are the dosimetric data measured on same clinac.
In the above-mentioned technical solutions, it is preferable that the first photon energy beam dosimetric data, the second energy photon beam dosimetry
Data and intermediate energy photon beam dosimetric data are the dosimetric data measured on different clinacs.
In any of the above-described technical solution, it is preferable that dosimetric data includes percent depth-dose curves and off-axis dosage
Distribution curve.
In any of the above-described technical solution, it is preferable that according to preset fitting coefficient to the first energy photon beam dosimetry
Data and the second photon energy beam dosimetric data carry out Mathematical Fitting, obtain synthesis photon energy beam dosimetric data, specifically
Formula are as follows:
PDDsyn=α PDDlow+β·PDDhigh;
OARsyn=α OARlow+β·OARhigh;
Wherein, syn indicates the corresponding ENERGY E of synthesis photon energy beam dosimetric datasyn, low the first energy photon of expression
The corresponding ENERGY E of beam dosimetry datalow, the corresponding ENERGY E of high the second photon energy beam dosimetric data of expressionhigh, α and β
Represent fitting coefficient.
In any of the above-described technical solution, it is preferable that synthesis photon energy beam dosimetric data and intermediate energy photon
Beam dosimetry data are compared verifying, obtain dosimetric data fitness bias, specific formula are as follows:
Wherein, mid indicates the corresponding ENERGY E of intermediate energy photon beam dosimetric datamid, RjFor the synthesis energy of jth step
Photon beam dosimetric data and intermediate energy photon beam dosimetric data fitness bias, WPDDAnd WOARIndicate the adjustable weight of user
Coefficient.
In any of the above-described technical solution, it is preferable that comparison result is that dosimetric data fitness bias is less than or equal in advance
If threshold value, specifically: Rj≤δ;
Comparison result is that dosimetric data fitness bias is greater than preset threshold, specifically: Rj> δ.
The technical solution of second aspect of the present invention provides a kind of photon energy synthesis system of clinac, packet
It includes: obtaining module, for obtaining the first photon energy beam dosimetric data, the second photon energy beam dosimetric data and correspondence
In the intermediate energy photon beam dosimetry number of the first photon energy beam dosimetric data and the second photon energy beam dosimetric data
According to;Mathematical Fitting module is used for according to preset fitting coefficient to the first photon energy beam dosimetric data and the second energy light
Beamlet dosimetric data carries out Mathematical Fitting and obtains synthesis photon energy beam dosimetric data;Authentication module, for synthesis energy
It measures photon beam dosimetric data and is compared verifying with intermediate energy photon beam dosimetric data, it is inclined to obtain dosimetric data fitting
Difference;Comparison module obtains comparison result for dosimetric data fitness bias to be compared with preset threshold;Wherein, when than
When relatively result is that dosimetric data fitness bias is less than or equal to preset threshold, synthesis photon energy beam dosimetric data is final
Dosimetric data;When comparison result is that dosimetric data fitness bias is greater than preset threshold, it is fitted according to dosimetric data inclined
Difference, which is iterated, optimizes and revises fitting coefficient.
In the above-mentioned technical solutions, it is preferable that the first photon energy beam dosimetric data, the second energy photon beam dosimetry
Data and intermediate energy photon beam dosimetric data are the dosimetric data of same clinac.
In the above-mentioned technical solutions, it is preferable that the first photon energy beam dosimetric data, the second energy photon beam dosimetry
Data and intermediate energy photon beam dosimetric data are the dosimetric data of different clinacs.
In any of the above-described technical solution, it is preferable that dosimetric data includes percent depth-dose curves and off-axis dosage
Distribution curve.
In any of the above-described technical solution, it is preferable that according to preset fitting coefficient to the first energy photon beam dosimetry
The dosimetric data of data and the second photon energy beam dosimetric data carries out Mathematical Fitting, obtains synthesis energy photon beam dose
Learn data, specific formula are as follows:
PDDsyn=α PDDlow+β·PDDhigh;
OARsyn=α OARlow+β·OARhigh;
Wherein, syn indicates the corresponding ENERGY E of synthesis photon energy beam dosimetric datasyn, low the first energy photon of expression
The corresponding ENERGY E of beam dosimetry datalow, the corresponding ENERGY E of high the second photon energy beam dosimetric data of expressionhigh, α and β
Represent default fitting coefficient;
Verifying is compared with intermediate energy photon beam dosimetric data to synthesis photon energy beam dosimetric data, is obtained
Dosimetric data fitness bias, specific formula are as follows:
Wherein, mid indicates the corresponding ENERGY E of intermediate energy photon beam dosimetric datamid, RjFor the synthesis energy of jth step
Photon beam dosimetric data and intermediate energy photon beam dosimetric data fitness bias, WPDDAnd WOARIndicate the adjustable weight of user
Coefficient;
Comparison result is that dosimetric data fitness bias is less than or equal to preset threshold, specifically: Rj≤δ;
Comparison result is that dosimetric data fitness bias is greater than preset threshold, specifically: Rj> δ.
Detailed description of the invention
Above-mentioned and/or additional aspect of the invention and advantage will become from the description of the embodiment in conjunction with the following figures
Obviously and it is readily appreciated that, in which:
Fig. 1 shows the flow chart element of the photon energy synthetic method of the clinac of some embodiments of the invention
Figure;
Fig. 2 shows the structural frames of the photon energy synthesis system of the clinac of some embodiments of the invention
Figure;
Fig. 3 is clinac 6 in the specific embodiment of the invention, 10, the PDD curve graph of 15MV photon energy beam;
Fig. 4 is the curve graph of fitting coefficient α in the specific embodiment of the invention, β value;
Fig. 5 is the PDD curve comparison figure of synthesis dosimetric data and intermediate dosimetric data in the embodiment of the present invention;
Fig. 6 is the OAR curve comparison figure of synthesis dosimetric data and intermediate dosimetric data in the embodiment of the present invention;
Fig. 7 is the dosage plane comparison diagram of synthesis dosimetric data and intermediate dosimetric data in the embodiment of the present invention;
Fig. 8 is pair of the dose volume histogram of synthesis dosimetric data and intermediate dosimetric data in the embodiment of the present invention
Than figure.
Specific embodiment
To better understand the objects, features and advantages of the present invention, with reference to the accompanying drawing and specific real
Applying mode, the present invention is further described in detail.It should be noted that in the absence of conflict, the implementation of the application
Feature in example and embodiment can be combined with each other.
In the following description, numerous specific details are set forth in order to facilitate a full understanding of the present invention, still, the present invention may be used also
To be implemented using other than the one described here other modes, therefore, protection scope of the present invention is not limited to following public affairs
The limitation for the specific embodiment opened.
The photon energy synthesis of clinac involved in the embodiment of the present invention is described referring to Fig. 1 to Fig. 8
Method and system.
As shown in Figure 1, the photon energy synthetic method of the clinac according to some embodiments of the invention, including
Following steps:
S100 obtains the corresponding first photon energy beam dosimetric data (percent depth-dose curves of low-energy light beamlet
With off-axis dose profile etc.), energy be higher than corresponding second energy photon of the high-energy light beamlet of the first photon energy beam
Beam dosimetry data and corresponding to centre corresponding to the intermediate energy photon beam of low-energy light beamlet and high-energy light beamlet
Photon energy beam dosimetric data;
S200, according to preset fitting coefficient to the first photon energy beam dosimetric data and the second energy photon beam dose
The dosimetric data for learning data carries out Mathematical Fitting, obtains the corresponding synthesis energy photon beam dosimetry number of synthesis photon energy beam
According to;
S300 is compared synthesis photon energy beam dosimetric data with intermediate energy photon beam dosimetric data and tests
Card, obtains dosimetric data fitness bias;
Dosimetric data fitness bias is compared with preset threshold, obtains comparison result by S400;
S500, when comparison result is that dosimetric data fitness bias is less than or equal to preset threshold, in step S200
Photon energy beam dosimetric data is synthesized as final dose data;
It, can be according to dosimetric data fitness bias when comparison result is that dosimetric data fitness bias is greater than preset threshold
It is iterated and optimizes and revises fitting coefficient, according to fitting coefficient adjusted to the first photon energy beam and the second photon energy beam
Mathematical Fitting is carried out, obtain and newly synthesizes photon energy beam similar in intermediate energy photon beam beam behaviour, repeats step S200-
S400, until the corresponding dosimetric data fitness bias of new synthesis photon energy beam is less than or equal to preset threshold.
In this embodiment, the ray of the high and low two grades of energy exported from existing clinac is according to above
Determining fitting coefficient mixes, ray similar in the beam behaviour of any energy in the above two energy range of simulated implementation,
It can be identical or another medical with different photon energies by the radiation treatment plan of known photon light beam to realize
It realizes on linear accelerator, is distributed with reaching close to identical dosage;Meanwhile outfit two is not needed on clinac
A above photon energy can save the spending and later standing charges of equipment research and development;May further be used to offer arbitrary continuation can
Photon energy is adjusted, provides maximum freedom degree for user, this is not accomplished on current hardware.
In the present embodiment, the first photon energy beam dosimetric data and the second photon energy beam dosimetric data are medical
The photon beam of the different-energy of linear accelerator, such as the first photon energy beam dosimetric data are one of clinac
Or multiple low energy ElowThe corresponding dosimetric data of photon beam, the second photon energy beam dosimetric data are one or more high
ENERGY EhighThe corresponding dosimetric data of photon beam;Intermediate energy photon beam dosimetric data can be clinac
The corresponding dosimetric data of one or more intermediate energy photon beams, and the energy value E of intermediate energy photon beammidPositioned at low energy
Measure the energy value E of photon beamlowWith the energy value E of high-energy light beamlethighBetween.
Wherein, dosimetric data be using 3 d water tank and corresponding measurement detector according under regulation geometrical condition to adding
What fast device relative dosage output measured, this set photon energy synthetic method is unrelated with planned treatment system, it is directly acted on
It is the foreground of planning system, it is possible to apply on all treatment planning systems in the data of water tank measurement.
In an embodiment of the invention, the first photon energy beam dosimetric data, the second energy photon beam dosimetry number
It is the dosimetric data measured on same clinac according to intermediate energy photon beam dosimetric data.
In another embodiment of the present invention, the first photon energy beam dosimetric data, the second energy photon beam dose
It learns data and intermediate energy photon beam dosimetric data is the dosimetric data measured on different clinacs.
In this embodiment, the first photon energy beam dosimetric data, the second photon energy beam dosimetric data and centre
Photon energy beam dosimetric data can be one or more dosimetric data measured on clinac;Or the
One photon energy beam dosimetric data, the second photon energy beam dosimetric data and intermediate energy photon beam dosimetric data can be with
For the dosimetric data measured on different clinacs;Or the first photon energy beam dosimetric data, the second energy
Photon beam dosimetric data and intermediate energy photon beam dosimetric data can be the agent measured on same clinac
Amount learns data;Or the first photon energy beam dosimetric data, the second photon energy beam dosimetric data and intermediate energy photon
At least two in beam dosimetry data can be the dosimetric data measured on same clinac.
In any of the above-described embodiment, it is preferable that dosimetric data includes percent depth-dose curves (percentage
Depth dose, PDD), off-axis dose profile (OffAxialRatio, OAR).
In this embodiment, percent depth-dose curves (percentage depth dose, PDD) are defined as in body film
Absorbed dose rate D on ray center axis at a certain depth ddWith reference depth d0Locate absorbed dose rate D0The ratio between percentage, be
The physical quantity of relative dosage distribution at ray center axis different depth is described.
In any of the above-described embodiment, it is preferable that according to preset fitting coefficient to the first energy photon beam dosimetry number
Mathematical Fitting is carried out according to the second photon energy beam dosimetric data, obtains synthesis photon energy beam dosimetric data, it is specific public
Formula are as follows:
PDDsyn=α PDDlow+β·PDDhigh;
OARsyn=α OARlow+β·OARhigh;
Wherein, syn indicates the corresponding ENERGY E of synthesis photon energy beam dosimetric datasyn, low the first energy photon of expression
The corresponding ENERGY E of beam dosimetry datalow, the corresponding ENERGY E of high the second photon energy beam dosimetric data of expressionhigh, α and β
Default fitting coefficient is represented, the initial value of default fitting coefficient can be set as (0.5,0.5).
In any of the above-described embodiment, it is preferable that synthesis photon energy beam dosimetric data and intermediate energy photon beam
Dosimetric data is compared verifying, obtains dosimetric data fitness bias, specific formula are as follows:
Wherein, mid indicates the corresponding ENERGY E of intermediate energy photon beam dosimetric datamid, RjFor the synthesis energy of jth step
Photon beam dosimetric data and intermediate energy photon beam dosimetric data fitness bias, WPDDAnd WOARIndicate the adjustable weight of user
Coefficient.
In any of the above-described embodiment, it is preferable that comparison result is dosimetric data fitness bias less than or equal to default
Threshold value, specifically: Rj≤δ;
Comparison result is that dosimetric data fitness bias is greater than preset threshold, specifically: Rj> δ.
As shown in Fig. 2, the photon energy synthesis system of the clinac of some embodiments according to the invention
1000, comprising: module 100 is obtained, for obtaining the corresponding first photon energy beam dosimetric data of low-energy light beamlet, energy
The second photon energy beam dosimetric data corresponding higher than the high-energy light beamlet of low-energy light beamlet and correspond to low energy
Intermediate energy photon beam dosimetric data corresponding to the intermediate energy photon beam of photon beam and high-energy light beamlet;Mathematical Fitting
Module 200 is used for according to preset fitting coefficient to the first photon energy beam dosimetric data and the second energy photon beam dose
It learns data and carries out Mathematical Fitting, obtain the corresponding synthesis photon energy beam dosimetric data of synthesis photon energy beam;Authentication module
300, for being compared verifying with intermediate energy photon beam dosimetric data to synthesis photon energy beam dosimetric data, obtain
Dosimetric data fitness bias;Comparison module 400 is obtained for dosimetric data fitness bias to be compared with preset threshold
To comparison result;Wherein, when comparison result is that dosimetric data fitness bias is less than or equal to preset threshold, energy photon is synthesized
Beam dosimetry data are final dose data;When comparison result is that dosimetric data fitness bias is greater than preset threshold, according to
Dosimetric data fitness bias, which is iterated, optimizes and revises fitting coefficient, according to fitting coefficient adjusted to the first energy photon
Beam and the second photon energy beam carry out Mathematical Fitting, obtain newly synthesizing energy light with similar in intermediate energy photon beam beam behaviour
Beamlet, until the corresponding dosimetric data fitness bias of new synthesis photon energy beam is less than or equal to preset threshold.
In this embodiment, which can be realized the light of the clinac as involved in above-mentioned any embodiment
The step of sub- energy synthetic method.
In the above embodiment, it is preferable that the first photon energy beam dosimetric data, the second energy photon beam dosimetry number
It is the dosimetric data measured on same clinac according to intermediate energy photon beam dosimetric data.
In the above embodiment, it is preferable that the first photon energy beam dosimetric data, the second energy photon beam dosimetry number
It is the dosimetric data measured on different clinacs according to intermediate energy photon beam dosimetric data.
In any of the above-described embodiment, it is preferable that dosimetric data includes percent depth-dose curves and off-axis dosage point
Cloth curve.
In any of the above-described embodiment, it is preferable that according to default fitting coefficient to the first photon energy beam dosimetric data
Mathematical Fitting is carried out with the dosimetric data of the second photon energy beam dosimetric data, obtains synthesis energy photon beam dosimetry number
According to specific formula are as follows:
PDDsyn=α PDDlow+β·PDDhigh;
OARsyn=α OARlow+β·OARhigh;
αj+1=αj+τ;
βj+1=βj+τ;
Wherein, syn indicates the corresponding ENERGY E of synthesis photon energy beam dosimetric datasyn, low the first energy photon of expression
The corresponding ENERGY E of beam dosimetry datalow, the corresponding ENERGY E of high the second photon energy beam dosimetric data of expressionhigh, α and β
Default fitting coefficient is represented, the initial value of default fitting coefficient can be set as (0.5,0.5), and iteration step length τ can be set as 0.05;
Verifying is compared with intermediate energy photon beam dosimetric data to synthesis photon energy beam dosimetric data, is obtained
Dosimetric data fitness bias, specific formula are as follows:
Wherein, mid indicates the corresponding ENERGY E of intermediate energy photon beam dosimetric datamid, RjFor the synthesis energy of jth step
Photon beam dosimetric data and intermediate energy photon beam dosimetric data fitness bias, WPDDAnd WOARIndicate the adjustable weight of user
Coefficient;
Comparison result is that dosimetric data fitness bias is less than or equal to preset threshold, specifically: Rj≤δ;
Comparison result is that dosimetric data fitness bias is greater than preset threshold, specifically: Rj> δ.
Specific embodiment
The photon beam of clinac can be in 4~18MV optional one to third gear, such as 6MV, 10MV and 15MV third gear.
The debugging of third gear energy measures the dosimetric data under each various geometrical conditions of shelves energy after checking and accepting, synthesize 10MV by 6MV and 15MV
Photon energy beam, specific implementation process are as follows:
1, the first photon energy beam dosimetric data, the second photon energy beam dosimetric data are obtained and corresponding to first
The intermediate energy photon beam dosimetric data of photon energy beam dosimetric data and the second photon energy beam dosimetric data;First
The corresponding low-energy light beamlet of photon energy beam dosimetric data and the corresponding high-energy of the second photon energy beam dosimetric data
Photon beam is respectively the photon beam of 6MV and 15MV, the photon that the corresponding intermediate energy photon beam of intermediate dosimetric data is 10MV
Beam;
2, correspond to low energy (6MV) photon beam, intermediate energy (10MV) photon beam and high-energy (15MV) photon beam
PDD curve, as shown in Figure 3.The initial value that α and β is arranged is (0.5,0.5), and iteration step length τ can be set as 0.05 conjunction being calculated
It is as shown in Figure 4 at the corresponding dosimetric data of photon energy beam.
3, PDD the and OAR data and medical linear for the synthesis 10_syn photon energy beam that each Mathematical Fitting optimizes
PDD the and OAR data of the known 10MV energy of accelerator measurement are fitted verifying, as shown in Figure 5 and Figure 6.
4, as shown in Figure 7 and Figure 8, in planning system, we are counted using 6MV and 15MV fitting 10MV energy-ray
Meter is delineated, and is intended to result and is compared with actual 10MV energy-ray plan.We have found that the dosage of two groups of plans point
The dose volume histogram result of cloth and each organ is almost the same.
In the present invention, term " first ", " second ", " third " are only used for the purpose of description, and should not be understood as indicating
Or imply relative importance;Term " multiple " then refers to two or more, unless otherwise restricted clearly.Term " installation ",
The terms such as " connected ", " connection ", " fixation " shall be understood in a broad sense, for example, " connection " may be a fixed connection, being also possible to can
Dismantling connection, or be integrally connected;" connected " can be directly connected, can also be indirectly connected through an intermediary.For this
For the those of ordinary skill in field, the specific meanings of the above terms in the present invention can be understood according to specific conditions.
In description of the invention, it is to be understood that the instructions such as term " on ", "lower", "left", "right", "front", "rear"
Orientation or positional relationship is to be based on the orientation or positional relationship shown in the drawings, and is merely for convenience of the description present invention and simplification is retouched
It states, rather than the device or unit of indication or suggestion meaning must have specific direction, be constructed and operated in a specific orientation,
It is thus impossible to be interpreted as limitation of the present invention.
In the description of this specification, the descriptions such as term " one embodiment ", " some embodiments ", " specific embodiment " are anticipated
Refer to that particular features, structures, materials, or characteristics described in conjunction with this embodiment or example are contained at least one implementation of the invention
In example or example.In the present specification, schematic expression of the above terms are not necessarily referring to identical embodiment or example.
Moreover, the particular features, structures, materials, or characteristics of description can be in any one or more of the embodiments or examples with suitable
Mode combine.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made any to repair
Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.
Claims (10)
1. a kind of photon energy synthetic method of clinac, which comprises the following steps:
It obtains the first photon energy beam dosimetric data, the second photon energy beam dosimetric data and corresponds to first energy
Measure the intermediate energy photon beam dosimetric data of photon beam dosimetric data and the second photon energy beam dosimetric data;
According to preset fitting coefficient to the first photon energy beam dosimetric data and the second energy photon beam dose
It learns data and carries out Mathematical Fitting, obtain synthesis photon energy beam dosimetric data;
Verifying is compared with the intermediate energy photon beam dosimetric data to the synthesis photon energy beam dosimetric data,
Obtain dosimetric data fitness bias;
The dosimetric data fitness bias is compared with preset threshold, obtains comparison result;
When the comparison result is that the dosimetric data fitness bias is less than or equal to the preset threshold, the synthesis energy
Amount photon beam dosimetric data is final dose data;
It is quasi- according to the dosimeter when the comparison result is that the dosimetric data fitness bias is greater than the preset threshold
Conjunction deviation, which is iterated, optimizes and revises the fitting parameter.
2. the photon energy synthetic method of clinac according to claim 1, it is characterised in that: described first
Photon energy beam dosimetric data, the second photon energy beam dosimetric data and the intermediate energy photon beam dosimetry number
According to the dosimetric data to be measured on same or different clinacs.
3. the photon energy synthetic method of clinac according to claim 2, it is characterised in that: the dosage
Learning data includes percent depth-dose curves and off-axis dose profile.
4. the photon energy synthetic method of clinac according to claim 3, it is characterised in that: according to default
Fitting coefficient to the first photon energy beam dosimetric data and the second photon energy beam dosimetric data mathematics into
Row fitting obtains synthesis photon energy beam dosimetric data, specific formula are as follows:
PDDsyn=α PDDlow+β·PDDhigh;
OARsyn=α OARlow+β·OARhigh;
Wherein, syn indicates the corresponding ENERGY E of the synthesis photon energy beam dosimetric datasyn, low expression first energy
The corresponding ENERGY E of photon beam dosimetric datalow, the corresponding energy of high expression the second photon energy beam dosimetric data
Ehigh, α and β represent the fitting coefficient.
5. the photon energy synthetic method of clinac according to claim 4, it is characterised in that: to the conjunction
It is compared verifying with the intermediate energy photon beam dosimetric data at photon energy beam dosimetric data, obtains dosimeter number
According to fitness bias, specific formula are as follows:
Wherein, mid indicates the corresponding ENERGY E of the intermediate energy photon beam dosimetric datamid, RjFor the synthesis of jth step
Photon energy beam dosimetric data and the intermediate energy photon beam dosimetric data fitness bias, WPDDAnd WOARIndicate that user can
The weight coefficient of tune.
6. the photon energy synthetic method of clinac according to claim 5, it is characterised in that: the comparison
As a result it is less than or equal to the preset threshold for the dosimetric data fitness bias, specifically: Rj≤δ;
The comparison result is that the dosimetric data fitness bias is greater than the preset threshold, specifically: Rj> δ.
7. a kind of photon energy synthesis system of clinac characterized by comprising
Module is obtained, for obtaining the first photon energy beam dosimetric data, the second photon energy beam dosimetric data and right
The intermediate energy photon of first photon energy beam dosimetric data and the second photon energy beam dosimetric data described in Ying Yu
Beam dosimetry data;
Mathematical Fitting module, for according to preset fitting coefficient to the first photon energy beam dosimetric data and described the
Two photon energy beam dosimetric datas carry out Mathematical Fitting, obtain synthesis photon energy beam dosimetric data;
Authentication module, for the synthesis photon energy beam dosimetric data and the intermediate energy photon beam dosimetric data
It is compared verifying, obtains dosimetric data fitness bias;
Comparison module obtains comparison result for the dosimetric data fitness bias to be compared with preset threshold;
Wherein, described when the comparison result is that the dosimetric data fitness bias is less than or equal to the preset threshold
Synthesis photon energy beam dosimetric data is final dose data;
When the comparison result is that the dosimetric data fitness bias is greater than the preset threshold, according to the dosimeter number
It is iterated according to fitness bias and optimizes and revises the fitting coefficient.
8. the photon energy synthesis system of clinac according to claim 7, it is characterised in that: described first
Photon energy beam dosimetric data, the second photon energy beam dosimetric data and the intermediate energy photon beam dosimetry number
According to the dosimetric data to be measured on same or different clinacs.
9. the photon energy synthesis system of clinac according to claim 8, it is characterised in that: the dosage
Learning data includes percent depth-dose curves and off-axis dose profile.
10. the photon energy synthesis system of clinac according to claim 9, it is characterised in that: according to pre-
If fitting coefficient the first photon energy beam dosimetric data and the second photon energy beam dosimetric data are carried out
Mathematical Fitting obtains synthesis photon energy beam dosimetric data, specific formula are as follows:
PDDsyn=α PDDlow+β·PDDhigh;
OARsyn=α OARlow+β·OARhigh;
Wherein, syn indicates the corresponding ENERGY E of the synthesis photon energy beam dosimetric datasyn, low expression first energy
The corresponding ENERGY E of photon beam dosimetric datalow, the corresponding energy of high expression the second photon energy beam dosimetric data
Ehigh, α and β represent the fitting coefficient;
Verifying is compared with the intermediate energy photon beam dosimetric data to the synthesis photon energy beam dosimetric data,
Obtain dosimetric data fitness bias, specific formula are as follows:
Wherein, mid indicates the corresponding ENERGY E of the intermediate energy photon beam dosimetric datamid, RjFor the synthesis of jth step
Photon energy beam dosimetric data and the intermediate energy photon beam dosimetric data fitness bias, WPDDAnd WOARIndicate that user can
The weight coefficient of tune;
The comparison result is that the dosimetric data fitness bias is less than or equal to the preset threshold, specifically: Rj≤δ;
The comparison result is that the dosimetric data fitness bias is greater than the preset threshold, specifically: Rj> δ.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811312180.1A CN109513118B (en) | 2018-11-06 | 2018-11-06 | Photon energy synthesis method and system of medical linear accelerator |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811312180.1A CN109513118B (en) | 2018-11-06 | 2018-11-06 | Photon energy synthesis method and system of medical linear accelerator |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109513118A true CN109513118A (en) | 2019-03-26 |
CN109513118B CN109513118B (en) | 2021-05-18 |
Family
ID=65774216
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811312180.1A Active CN109513118B (en) | 2018-11-06 | 2018-11-06 | Photon energy synthesis method and system of medical linear accelerator |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109513118B (en) |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1237079A (en) * | 1999-06-25 | 1999-12-01 | 清华大学 | Energy switch for axis-coupled standing wave accelerator tube |
US20020168048A1 (en) * | 2001-03-13 | 2002-11-14 | Wong James R. | Method for producing a range of therapeutic radiation energy levels |
CN1455635A (en) * | 2003-06-06 | 2003-11-12 | 南京大学 | Energy-adjustable electronic linear accelerator |
GB0505090D0 (en) * | 2005-03-12 | 2005-04-20 | Elekta Ab | Linear accelerator |
US20060033044A1 (en) * | 2004-08-04 | 2006-02-16 | Gentry John R | Treatment planning tool for multi-energy electron beam radiotherapy |
EP2248551A1 (en) * | 2009-05-05 | 2010-11-10 | 7Sigma N.V. | Method for the verification of a radiotherapy treatment apparatus |
CN202019491U (en) * | 2011-01-20 | 2011-10-26 | 孙安 | Standing wave irradiation accelerator with continuously variable energy |
CN102921115A (en) * | 2012-10-25 | 2013-02-13 | 合肥工业大学 | Method for establishing measurement data-based simple and convenient irradiation source model of medical linear accelerator |
CN102938969A (en) * | 2012-11-25 | 2013-02-20 | 中国原子能科学研究院 | Energy regulation method for traveling wave electron linear accelerator |
CN203233589U (en) * | 2012-12-28 | 2013-10-09 | 清华大学 | Standing wave electron linear accelerator apparatus |
CN104244562A (en) * | 2013-06-12 | 2014-12-24 | 梅维昂医疗系统股份有限公司 | Particle accelerator that produces charged particles having variable energies |
CN104470192A (en) * | 2013-09-22 | 2015-03-25 | 同方威视技术股份有限公司 | Electron linear accelerator system |
CN105288873A (en) * | 2015-12-04 | 2016-02-03 | 金乔 | Radiotherapy dosage rapid quantitative method |
CN105457170A (en) * | 2014-08-20 | 2016-04-06 | 上海联影医疗科技有限公司 | Electron linear accelerator source model building method and building device |
CN105916285A (en) * | 2016-06-30 | 2016-08-31 | 中国科学院上海应用物理研究所 | Dual-energy electron linear accelerator based on hybrid traveling wave accelerating structure |
CN106879158A (en) * | 2017-03-16 | 2017-06-20 | 东莞中子科学中心 | Medical proton linear accelerator |
US9910166B2 (en) * | 2008-05-22 | 2018-03-06 | Stephen L. Spotts | Redundant charged particle state determination apparatus and method of use thereof |
-
2018
- 2018-11-06 CN CN201811312180.1A patent/CN109513118B/en active Active
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1237079A (en) * | 1999-06-25 | 1999-12-01 | 清华大学 | Energy switch for axis-coupled standing wave accelerator tube |
US20020168048A1 (en) * | 2001-03-13 | 2002-11-14 | Wong James R. | Method for producing a range of therapeutic radiation energy levels |
CN1455635A (en) * | 2003-06-06 | 2003-11-12 | 南京大学 | Energy-adjustable electronic linear accelerator |
US20060033044A1 (en) * | 2004-08-04 | 2006-02-16 | Gentry John R | Treatment planning tool for multi-energy electron beam radiotherapy |
US7202486B2 (en) * | 2004-08-04 | 2007-04-10 | Standard Imaging, Inc. | Treatment planning tool for multi-energy electron beam radiotherapy |
GB0505090D0 (en) * | 2005-03-12 | 2005-04-20 | Elekta Ab | Linear accelerator |
US9910166B2 (en) * | 2008-05-22 | 2018-03-06 | Stephen L. Spotts | Redundant charged particle state determination apparatus and method of use thereof |
EP2248551A1 (en) * | 2009-05-05 | 2010-11-10 | 7Sigma N.V. | Method for the verification of a radiotherapy treatment apparatus |
CN202019491U (en) * | 2011-01-20 | 2011-10-26 | 孙安 | Standing wave irradiation accelerator with continuously variable energy |
CN102921115A (en) * | 2012-10-25 | 2013-02-13 | 合肥工业大学 | Method for establishing measurement data-based simple and convenient irradiation source model of medical linear accelerator |
CN102938969A (en) * | 2012-11-25 | 2013-02-20 | 中国原子能科学研究院 | Energy regulation method for traveling wave electron linear accelerator |
CN203233589U (en) * | 2012-12-28 | 2013-10-09 | 清华大学 | Standing wave electron linear accelerator apparatus |
CN104244562A (en) * | 2013-06-12 | 2014-12-24 | 梅维昂医疗系统股份有限公司 | Particle accelerator that produces charged particles having variable energies |
CN104470192A (en) * | 2013-09-22 | 2015-03-25 | 同方威视技术股份有限公司 | Electron linear accelerator system |
CN105457170A (en) * | 2014-08-20 | 2016-04-06 | 上海联影医疗科技有限公司 | Electron linear accelerator source model building method and building device |
CN105288873A (en) * | 2015-12-04 | 2016-02-03 | 金乔 | Radiotherapy dosage rapid quantitative method |
CN105916285A (en) * | 2016-06-30 | 2016-08-31 | 中国科学院上海应用物理研究所 | Dual-energy electron linear accelerator based on hybrid traveling wave accelerating structure |
CN106879158A (en) * | 2017-03-16 | 2017-06-20 | 东莞中子科学中心 | Medical proton linear accelerator |
Non-Patent Citations (1)
Title |
---|
夏祥等: "Varian RapidARC直线加速器6 MeV-X线虚拟源模型的建立", 《中国医学物理学杂志》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109513118B (en) | 2021-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ma et al. | Monte Carlo verification of IMRT dose distributions from a commercial treatment planning optimization system | |
US7987053B2 (en) | Monitor units calculation method for proton fields | |
AU2009201299A1 (en) | A programmable particle scatterer for radiation therapy beam formation | |
Arnfield et al. | Dosimetric validation for multileaf collimator-based intensity-modulated radiotherapy: a review | |
Hauri et al. | A general model for stray dose calculation of static and intensity‐modulated photon radiation | |
US20200094073A1 (en) | Methods to optimize coverage for multiple targets simultaneously for radiation treatments | |
CN105866821A (en) | Method and device for obtaining beam energy spectrum, and dose distribution calculating method | |
Tegaw et al. | Dosimetric characteristics of the INTRABEAM® system with spherical applicators in the presence of air gaps and tissue heterogeneities | |
Vidal et al. | Fast optimized Monte Carlo phase-space generation and dose prediction for low energy x-ray intra-operative radiation therapy | |
Asadi et al. | Development and validation of an optimal GATE model for proton pencil-beam scanning delivery | |
Fischer et al. | Commissioning of a double‐focused micro multileaf collimator (μMLC) | |
CN106061555A (en) | Particle beam treatment-planning apparatus and method for simulating particle beam irradiation | |
CN109513118A (en) | A kind of the photon energy synthetic method and system of clinac | |
Bedford et al. | A phase space model of a Versa HD linear accelerator for application to Monte Carlo dose calculation in a real‐time adaptive workflow | |
Fielding | Monte-Carlo techniques for radiotherapy applications II: equipment and source modelling, dose calculations and radiobiology | |
Badali et al. | Characterization of an x-ray source with a partitioned diamond-tungsten target for electronic brachytherapy with 3D beam directionality | |
Chen et al. | TOPAS Monte Carlo simulation for a scanning proton therapy system in SPHIC | |
Lazarus et al. | Validation of Monte Carlo-based calculations for megavolt electron beams for IORT and FLASH-IORT | |
O’Shea et al. | Monte Carlo simulation of medical accelerator electron treatment heads | |
Amour et al. | Verification of Depth Dose Curves Derived on Beeswax, Paraffin and Water Phantoms Using FLUKA Monte Carlo Code | |
Kirby | Radiation dosimetry of conventional and laser-driven particle beams | |
Zhuang et al. | Generating arbitrary one-dimensional dose profiles using rotational therapy | |
Neupane | Efficacy of the Virtual Cone Method Using Fixed Small Multi-Leaf Collimator Field for Stereotactic Radiosurgery | |
Galal et al. | Dosimetric and mechanical characteristics of a commercial dynamic μMLC used in SRS | |
Mavroidis et al. | Expected clinical impact of the differences between planned and delivered dose distributions in helical tomotherapy for treating head and neck cancer using helical megavoltage CT images |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20211215 Address after: North Carolina Patentee after: DUKE University Address before: 100191 jimenli hejingyuan, Haidian District, Beijing Patentee before: Wu Qiuwen Patentee before: Liu Bo Patentee before: Zhou Fugen |