CN108671419A - A kind of calculating of quick afterloading radiotherapy dosage, verification method - Google Patents
A kind of calculating of quick afterloading radiotherapy dosage, verification method Download PDFInfo
- Publication number
- CN108671419A CN108671419A CN201810591778.2A CN201810591778A CN108671419A CN 108671419 A CN108671419 A CN 108671419A CN 201810591778 A CN201810591778 A CN 201810591778A CN 108671419 A CN108671419 A CN 108671419A
- Authority
- CN
- China
- Prior art keywords
- dosage
- radioactive source
- calculating
- dose
- rapid
- 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/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
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)
- Radiation-Therapy Devices (AREA)
Abstract
The invention discloses a kind of calculating of quick afterloading radiotherapy dosage, verification methods, and the corresponding DICOM file of patient is exported from treatment planning systems, are saved in a glue file folder；Glue file folder is scanned, required Rapid Dose Calculation information is read from DICOM file；Radioactive source parameter is obtained, Distribution of dose rate table is established, determines radioactive source direction, and carry out Rapid Dose Calculation；Dose Results in the Dose Results and treatment planning systems of calculating are compared, deviation Dev and γ verification result is obtained；After the completion of comparison, by deviation Dev and γ verification result preserve to glue file press from both sides in, and be automatically deleted glue file folder in DICOM file.The present invention can carry out dosage verifying before the treatment, simply, quickly, improve accuracy and efficiency, entire verification process only needs a few minutes, too many human-computer interaction is not needed, and can be used for various Afterloading radiotherapy planning systems, quality control and guarantee to Afterloading radiotherapy are improved, suitable for promoting.
Description
Technical field
The present invention relates to Afterloading radiotherapy technical fields, and in particular to a kind of calculating of quick afterloading radiotherapy dosage is tested
Card method.
Background technology
Afterloading radiotherapy is a kind of effective cancer treatment method, and the accuracy for assessing Rapid Dose Calculation is accurate to ensureing
Afterloading radiotherapy is most important；The flow of radiotherapy is after treatment plan approval, to treatment plan before implementing to start
Carry out dosage verifying.
External physics teacher or dosage teacher (should by treatment planning systems (Treatment Planning System, TPS)
Treatment planning systems are the existing software that physics teacher or dosage teacher are used to design patient radiation's treatment plan) it completes treatment plan and sets
After meter, it will usually there is a dosage verifying of other one veteran physics Shi Jinhang treatment plan, but China physics teacher
Clinical task is heavier, and it is difficult to realize to carry out dosage verifying to treatment plan by an other physics teacher.European radiotherapy association
It can just recommend each Afterloading radiotherapy center that should have a set of dosage independently of other than treatment planning systems early in 2004
Verification method, but China yet there are no the research of Afterloading radiotherapy dosage verifying method.
Invention content
Based on this, in view of the above-mentioned problems, it is necessary to propose a kind of simplified Afterloading radiotherapy flow, improve Rapid Dose Calculation,
Verifying speed, and ensure the calculating of the quick afterloading radiotherapy dosage of Afterloading radiotherapy accuracy, verification method.
The technical scheme is that：
A kind of calculating of quick afterloading radiotherapy dosage, verification method, include the following steps：
A, the corresponding DICOM file of patient is exported from treatment planning systems, is saved in a glue file folder；
B, glue file folder is scanned, required Rapid Dose Calculation information is read from DICOM file；
C, radioactive source parameter is obtained, Distribution of dose rate table is established；
D, it determines radioactive source direction, and Rapid Dose Calculation is carried out according to Rapid Dose Calculation information, obtain Rapid Dose Calculation result；
E, the Dose Results in the Rapid Dose Calculation result and treatment planning systems of acquisition are compared, obtain deviation
Dev and γ verification results；
F, after the completion of comparison, deviation Dev and γ verification result is preserved to glue file in pressing from both sides, and be automatically deleted fixation
DICOM file in file.
In the technical scheme, from the point of view of physics teacher or dosage teacher, automatic business processing as far as possible reduces people
The intervention of work, saves the time；By the result of calculation for the treatment of planning systems (Treatment Planning System, TPS) with
The relatively independent result of calculation of the technical program, which is compared, obtains deviation Dev and γ verification result, then according to deviation
Dev and γ verification results judge whether dosage is reasonable, avoid that malpractice occurs；The technical program is importing treatment planning systems
The parked position of middle radioactive source and after the time can calculate three-dimensional dosage distribution, be read by comparing in treatment planning systems automatically
The calculated dosage distribution of the dosage distribution taken and the technical program, which can detect, is usually less susceptible to the mistake being found, such as
The accidental change of radioactive source parameter database or damage, radioactive source calibrate date and activity input error etc., Jin Erbao after replacing
Demonstrate,prove safety；Moreover, after the completion of comparing, the technical program can be output to result in the file of user's name automatically, in case will
Come DICOM (the Digital Imaging and Communications in for verifying, and being automatically deleted in glue file folder
Medicine) file (i.e. the standard format files of medical image), avoiding data redundancy from causing, memory burden is big and data go out
It is wrong.
Preferably, the data acquisition radioactive source parameter recommended according to AAPM and ESTRO in the step c.In order to ensure this
The independence of software Rapid Dose Calculation, radioactive source parameter are not obtained from radiotherapy treatment planning system, and use AAPM (American
Association of Physicists in Medicine) and ESTRO (European Society for
Radiotherapy and Oncology) recommend data.
Preferably, Distribution of dose rate table is established in the step c to include the following steps：
Calculate the dosage rate at a certain Rapid Dose Calculation pointWherein,For the dosage rate at Rapid Dose Calculation point, r is Rapid Dose Calculation point to the distance of radiation source center, r0=1cm, θ are polar coordinate system
Middle dosage calculates the angle between point and the directions radioactive source long axis V, θ0=pi/2, SkFor air kerma intensity, Λ is dosage
Rate constant, G are geometrical factor, and g is radial dose function, and F is anisotropy function；
According to the symmetry that dosage around radioactive source is distributed, calculate along the radioactive source directions long axis V and perpendicular to radioactive source
The two-dimensional Distribution of dose rate table T (m, n) in the directions long axis U, and deposit in calculator memory.
In the technical scheme, it is one that the formula recommended according to AAPM TG-43 (Task Group 43) reports calculates certain
Amount calculates the dosage rate at pointThe formula that the AAPM TG-43 (Task Group 43) reports are recommended, which is this field, routinely to be made
Formula, according to the calculated dosage rateIt can obtain the dosage for the treatment of planning systems, dose ratio is carried out convenient for follow-up
It is right；Moreover, in order to improve Rapid Dose Calculation speed under the premise of ensureing Rapid Dose Calculation precision, it is distributed in conjunction with dosage around radioactive source
Symmetry, we first calculate one along radioactive source long axis direction (directions V) and perpendicular to putting before Rapid Dose Calculation starts
A two-dimensional Distribution of dose rate table T (m, n) for penetrating source long axis direction (directions U) is stored in calculator memory, wherein dosage
Resolution ratio of the rate distribution table in the directions U and the directions V all takes 0.1cm, the ranging from 20cm in the directions U, and the range in the directions V is from -20cm
To 20cm.
Preferably, determine that radioactive source direction includes the following steps in the step d：
I-th of parked position of radioactive source is set as Si(xi, yi, zi), and pass through parked position Si(xi, yi, zi) and under
One parked position Si+1(xi+1, yi+1, zi+1) composition vectorDetermine direction of the radioactive source in human body coordinate system,
And vectorDue to the use of radioactive source be line source, line source need
Consider the difference of radioactive source coordinate system and human body coordinate system, so, direction of the radioactive source in human body coordinate system is calculated, is convenient for
Carry out subsequent dose calculating.
Preferably, before carrying out Rapid Dose Calculation in the step d, it is P (x, y, z), applicating device pipe that setting dosage, which calculates point,
Road sum is NA, radioactive source parked position sum is N in every applicating device pipelineS, i-th in jth root applicating device pipeline be resident
Position is d to the dosage rate contribution at P (x, y, z)J, i, the residence time of i-th of parked position in jth root applicating device pipeline
For tJ, i.The preparation for carrying out Rapid Dose Calculation, according to the Rapid Dose Calculation information read from derived DICOM file so that with treatment
The data of planning system are identical, calculate same dose calculate point Dose Results, make the accuracy higher of calculating, convenient for subsequently into
Row compares.
Preferably, the Rapid Dose Calculation at P (x, y, z) is included the following steps in the step d：
D101, P (x, y, z) i-th away from radioactive source parked position S is calculatedi(xi, yi, zi) distance r ', then
D102, according to vector
It is calculatedWithBetween angle
D103, i-th away from radioactive source parked position S of P (x, y, z) is calculated separatelyi(xi, yi, zi) the directions V distance v and
The distance u in the directions U, wherein v=r ' cos θ,
D104, Distribution of dose rate table T (m, n) is checked, if m × 0.1cm≤u ＜ (m+1) × 0.1cm and n × 0.1cm≤v
＜ (n+1) × 0.1cm, then to dJ, i(u, v) corresponding Distribution of dose rate carries out bilinear interpolation calculating, obtains dJ, i(u, v)=
w1T (m, n)+w2T (m, n+1)+w3T (m+1, n)+w4T (m+1, n+1)
, wherein w1、w2、w3And w4T (m, n) respectively in Distribution of dose rate table, T (m, n+1), T (m+1, n) and T (m+1,
N+1 the weight at)；
Dosage D (x, y, z) at d105, calculating P (x, y, z),
In the technical scheme, each radioactive source contributes d to the dosage rate at P (x, y, z)J, i(u, v) is depended in memory
Obtained after the data simple interpolations of extraction, without using AAPM TG-43 (Task Group 43) report recommend formula into
Row computes repeatedly；This programme Rapid Dose Calculation sizing grid and treatment planning systems it is consistent, be all 0.1 × 0.1 × 0.1cm3, examine
The dosage for considering 0.5cm ranges around radioactive source is all very high, this subregion nor clinical dosage focus of attention, so
It is assumed that being equal to away from the dosage at radioactive source 0.5cm away from the dosage in radioactive source 0.5cm；And then ensureing the same of counting accuracy
When, improve calculating speed.
Preferably, deviation Dev is calculated in the step e to include the following steps：
E101, dosimetric parameters D is obtained from the RT Dose for the treatment of planning systemsTPS；
Dosimetric parameters D is obtained after e102, Rapid Dose CalculationQA；
E103, according to dosimetric parameters DTPSWith dosimetric parameters DQADeviation Dev is calculated, is obtained
In the technical scheme, the dosage distribution that treatment planning systems calculate is read from RT Dose, by itself and this programme
The dosage distribution of calculating is compared, we define DX%For the exposure dose that the volume of organ x% is subject to, DyccFor organ y cm3
The exposure dose that volume is subject to, recommends according to ESTRO, counts target area D100%, D90%, normal organ D0.1cc、D1ccAnd D2cc；From
And dosimetric parameters D is counted respectivelyQAWith dosimetric parameters DQA, then deviated value Dev according to contrast conting；Which energy
The accuracy that Rapid Dose Calculation is filled after rapid evaluation, ensures the accuracy of Afterloading radiotherapy.
Preferably, γ verifications include the following steps in the step e：：
E101, dosage is read from the RT Dose for the treatment of planning systems；
E102, using the dosage after Rapid Dose Calculation as standard, calculate γ values；
If e103, γ value are more than 1, obstructed overtreatment verification passes through dosage verifying if γ values are less than or equal to 1.
The γ is verified as another verification method, and for conventional verification method, the computational methods of wherein γ values are shown in document
(Low D A,Harms W B,Mutic S,et al.A technique for the quantitative evaluation
of dose distributions[J].Medical Physics,1998,25(5):656-661.), main is exactly with we
The calculated dosage of case is standard, is read with treatment planning systems RT Dose (file for storing dosage in DICOM file)
Dosage be compared, judge whether to be verified；Can be not only according to deviation Dev, but also dosage judged according to γ verifications
Accuracy ensures the accuracy of Afterloading radiotherapy, further increases safety, avoids that malpractice occurs.
The beneficial effects of the invention are as follows：
The present invention provides a kind of rear dress dosage verifying tools that is simple and quick, can carrying out dosage verifying before the treatment, carry
High verification accuracy and efficiency, entire verification process only need a few minutes, reduce the stand-by period of patient, need not be too many
Human-computer interaction has been saved time and cost of labor, and can be used for various Afterloading radiotherapy planning systems, is put to penetrate for after and be controlled
The quality control and quality assurance for the treatment of play the role of reinforcement, promotion, greatly reduce malpractice incidence, suitable for pushing away
Extensively, very valuable.
Description of the drawings
Fig. 1 is the flow chart of the calculating of quick afterloading radiotherapy dosage, verification method described in the embodiment of the present invention；
Fig. 2 is Rapid Dose Calculation coordinate system schematic diagram described in the embodiment of the present invention；
Fig. 3 is for CTV described in the embodiment of the present invention in D100%In the case for the treatment of planning systems and this programme dosimeter
Parameter differences schematic diagram；
Fig. 4 is for CTV described in the embodiment of the present invention in D90%In the case of the dosimeter for the treatment of planning systems and this programme join
Number difference schematic diagram；
Fig. 5 is for bladder described in the embodiment of the present invention in D0.1ccIn the case for the treatment of planning systems and this programme dosimeter
Parameter differences schematic diagram；
Fig. 6 is for bladder described in the embodiment of the present invention in D1ccIn the case for the treatment of planning systems and this programme dosimeter
Parameter differences schematic diagram；
Fig. 7 is for bladder described in the embodiment of the present invention in D2ccIn the case for the treatment of planning systems and this programme dosimeter
Parameter differences schematic diagram；
Fig. 8 is for rectum described in the embodiment of the present invention in D0.1ccIn the case for the treatment of planning systems and this programme dosimeter
Parameter differences schematic diagram；
Fig. 9 is for rectum described in the embodiment of the present invention in D1ccIn the case for the treatment of planning systems and this programme dosimeter
Parameter differences schematic diagram；
Figure 10 is for rectum described in the embodiment of the present invention in D2ccIn the case for the treatment of planning systems and this programme dosimeter
Parameter differences schematic diagram；
Figure 11 is for small intestine described in the embodiment of the present invention in D0.1ccIn the case for the treatment of planning systems and this programme dosage
Learn parameter differences schematic diagram；
Figure 12 is for small intestine described in the embodiment of the present invention in D1ccIn the case for the treatment of planning systems and this programme dosimeter
Parameter differences schematic diagram；
Figure 13 is for small intestine described in the embodiment of the present invention in D2ccIn the case for the treatment of planning systems and this programme dosimeter
Parameter differences schematic diagram；
Figure 14 is the flow chart that γ verifications are carried out described in the embodiment of the present invention.
Specific implementation mode
The embodiment of the present invention is described in detail below in conjunction with the accompanying drawings.
Embodiment 1
As shown in Figure 1, a kind of calculating of quick afterloading radiotherapy dosage, verification method, include the following steps：
A, the corresponding DICOM file of patient is exported from treatment planning systems, is saved in a glue file folder；
B, glue file folder is scanned, required Rapid Dose Calculation information is read from DICOM file；
C, radioactive source parameter is obtained, Distribution of dose rate table is established；
D, it determines radioactive source direction, and Rapid Dose Calculation is carried out according to Rapid Dose Calculation information, obtain Rapid Dose Calculation result；
E, the Dose Results in the Rapid Dose Calculation result and treatment planning systems of acquisition are compared, obtain deviation
Dev and γ verification results；
F, after the completion of comparison, deviation Dev and γ verification result is preserved to glue file in pressing from both sides, and be automatically deleted fixation
DICOM file in file.
In the present embodiment, from the point of view of physics teacher or dosage teacher, automatic business processing as far as possible reduces artificial
Intervention, save the time；Pass through the result of calculation for the treatment of planning systems (Treatment Planning System, TPS) and this
The relatively independent result of calculation of technical solution, which is compared, obtains deviation Dev and γ verification result, then according to deviation Dev
Judge whether dosage is reasonable with γ verification results, avoids that malpractice occurs；The present embodiment is put in importing treatment planning systems
It penetrates the parked position in source and after the time, three-dimensional dosage distribution can be calculated automatically, pass through what is read in comparison treatment planning systems
Dosage is distributed and the calculated dosage distribution of the technical program can detect and usually be less susceptible to the mistake being found, for example is radiated
The accidental change of source parameter database or damage, radioactive source is calibrated date and activity input error etc. after replacing, and then ensures peace
Quan Xing；Moreover, after the completion of comparing, the present embodiment can be output to result in the file of user's name automatically, in case verify in the future,
And it is automatically deleted the DICOM (Digital Imaging and Communications in Medicine) in glue file folder
File avoids data redundancy from causing memory burden big and corrupt data.
Embodiment 2
The present embodiment is put in the step c according to the AAPM and ESTRO data acquisitions recommended on the basis of embodiment 1
Penetrate source parameter.In order to ensure that the independence of the software Rapid Dose Calculation, radioactive source parameter are not obtained from radiotherapy treatment planning system,
And use AAPM (American Association of Physicists in Medicine) and ESTRO (European
Society for Radiotherapy and Oncology) recommend data.
Embodiment 3
The present embodiment is established Distribution of dose rate table on the basis of embodiment 2, in the step c and is included the following steps：
Calculate the dosage rate at a certain Rapid Dose Calculation pointWherein,For the dosage rate at Rapid Dose Calculation point, r is Rapid Dose Calculation point to the distance of radiation source center, r0=1cm, θ are polar coordinate system
Middle dosage calculates the angle between point and the directions radioactive source long axis V, θ0=pi/2, SkFor air kerma intensity, Λ is dosage
Rate constant, G are geometrical factor, and g is radial dose function, and F is anisotropy function；
According to the symmetry that dosage around radioactive source is distributed, calculate along the radioactive source directions long axis V and perpendicular to radioactive source
The two-dimensional Distribution of dose rate table T (m, n) in the directions long axis U, and deposit in calculator memory.
In the present embodiment, the formula recommended according to AAPM TG-43 (Task Group 43) reports calculates certain dose
Calculate the dosage rate at pointThe formula that the AAPM TG-43 (Task Group 43) reports are recommended is commonly used in the art
Formula, according to the calculated dosage rateIt can obtain the dosage for the treatment of planning systems, dose comparing is carried out convenient for follow-up；And
And in order to improve Rapid Dose Calculation speed under the premise of ensureing Rapid Dose Calculation precision, in conjunction with pair of dosage distribution around radioactive source
Title property, we first calculate one along radioactive source long axis direction (directions V) and perpendicular to radioactive source before Rapid Dose Calculation starts
One two-dimensional Distribution of dose rate table T (m, n) of long axis direction (directions U) is stored in calculator memory, wherein dosage rate point
Resolution ratio of the cloth table in the directions U and the directions V all takes 0.1cm, the ranging from 20cm in the directions U, the range in the directions V from -20cm to
20cm。
Embodiment 4
The present embodiment is on the basis of embodiment 2, as shown in Fig. 2, determining that radioactive source direction includes following in the step d
Step：
I-th of parked position of radioactive source is set as Si(xi, yi, zi), and pass through parked position Si(xi, yi, zi) and under
One parked position Si+1(xi+1, yi+1, zi+1) composition vectorDetermine direction of the radioactive source in human body coordinate system,
And vectorDue to the use of radioactive source be line source, line source need
Consider the difference of radioactive source coordinate system and human body coordinate system, so, direction of the radioactive source in human body coordinate system is calculated, is convenient for
Carry out subsequent dose calculating.
Embodiment 5
The present embodiment is on the basis of embodiment 1, and before carrying out Rapid Dose Calculation in the step d, setting dosage calculates point P
(x, y, z), applicating device pipeline sum are NA, radioactive source parked position sum is N in every applicating device pipelineS, jth root applicating device
I-th of parked position in pipeline is d to the dosage rate contribution at P (x, y, z)J, i, stay for i-th in jth root applicating device pipeline
It is t to stay the residence time of positionJ, i。
Embodiment 6
The present embodiment is on the basis of embodiment 4, and before carrying out Rapid Dose Calculation in the step d, setting dosage calculates point P
(x, y, z), applicating device pipeline sum are NA, radioactive source parked position sum is N in every applicating device pipelineS, jth root applicating device
I-th of parked position in pipeline is d to the dosage rate contribution at P (x, y, z)J, i, stay for i-th in jth root applicating device pipeline
It is t to stay the residence time of positionJ, i。
Carry out Rapid Dose Calculation preparation, according to the Rapid Dose Calculation information read from derived DICOM file so that with control
The data for treating planning system are identical, calculate the Dose Results that same dose calculates point, make the accuracy higher of calculating, convenient for follow-up
It is compared.
Embodiment 7
The present embodiment includes following to the Rapid Dose Calculation at P (x, y, z) in the step d on the basis of embodiment 6
Step：
D101, P (x, y, z) i-th away from radioactive source parked position S is calculatedi(xi, yi, zi) distance r ', then
D102, according to vector
It is calculatedWithBetween angle
D103, i-th away from radioactive source parked position S of P (x, y, z) is calculated separatelyi(xi, yi, zi) the directions V distance v and
The distance u in the directions U, wherein v=r ' cos θ,
D104, Distribution of dose rate table T (m, n) is checked, if m × 0.1cm≤u ＜ (m+1) × 0.1cm and n × 0.1cm≤v
＜ (n+1) × 0.1cm, then to dJ, i(u, v) corresponding dosage distribution carries out bilinearity transplanting calculating, obtains dJ, i(u, v)=w1T
(m, n)+w2T (m, n+1)+w3T (m+1, n)+w4T (m+1, n+1)
, wherein w1、w2、w3And w4T (m, n) respectively in Distribution of dose rate table, T (m, n+1), T (m+1, n) and T (m+1,
N+1 the weight at)；
Dosage D (x, y, z) at d105, calculating P (x, y, z),
In the present embodiment, each radioactive source contributes d to the dosage rate at P (x, y, z)J, i(u, v) is depended in memory and is carried
It is obtained after the data simple interpolations taken, without being computed repeatedly using the formula in treatment planning systems；This programme dosage
The consistent of sizing grid and treatment planning systems is calculated, is all 0.1 × 0.1 × 0.1cm3, it is contemplated that 0.5cm models around radioactive source
The dosage enclosed is all very high, this subregion is nor clinical dosage focus of attention, it is assumed that away from radioactive source 0.5cm
Dosage is equal to away from the dosage at radioactive source 0.5cm；And then while ensureing counting accuracy, calculating speed is improved.
Embodiment 8
The present embodiment calculates deviation Dev on the basis of embodiment 7, in the step e and includes the following steps：
E101, dosimetric parameters D is obtained from the RT Dose for the treatment of planning systemsTPS；
Dosimetric parameters D is obtained after e102, Rapid Dose CalculationQA；
E103, according to dosimetric parameters DTPSWith dosimetric parameters DQADeviation Dev is calculated, is obtained
In the present embodiment, the dosage distribution that treatment planning systems calculate is read from RT Dose, by itself and this programme meter
The dosage distribution of calculation is compared, we define DX%For the exposure dose that the volume of organ x% is subject to, DyccFor organ y cm3Body
The exposure dose that product is subject to, recommends according to ESTRO, counts target area D100%, D90%, normal organ D0.1cc、D1ccAnd D2cc；To,
Dosimetric parameters D is counted respectivelyQAWith dosimetric parameters DQA, then deviated value Dev according to contrast conting；Which can be fast
The accuracy that Rapid Dose Calculation is filled after speed assessment, ensures the accuracy of Afterloading radiotherapy；
For CTV (clinical tumor target area), the dosimetric parameters difference for the treatment of planning systems and this programme of the present invention is shown in Fig. 3
And Fig. 4, Fig. 2 D100%Difference value, wherein maximum value is 1.76%, minimum value 0.08%, and the difference of 20 patients is average
Value is 0.85%；Fig. 3 is D90%Difference value, wherein maximum value be 1.37%, minimum value 0.15%, the difference of 20 patients
Average value is 0.70%；
For jeopardizing organ, treatment planning systems and this programme dosimetric parameters difference of the present invention are shown in Fig. 5-Figure 13, choose
20 patients, bladder and small intestine D0.1ccMaximum deviation be respectively less than 1.10%, rectum D0.1ccMaximum deviation be 1.41%,
Bladder, rectum and small intestine D1ccDifference average value be respectively 0.49%, 0.66% and 0.52%, D2ccDifference be less than D1ccWith
D0.1cc。
Embodiment 9
The present embodiment is on the basis of embodiment 7, and as shown in figure 14, γ is verified and included the following steps in the step e：：
E101, dosage is read from the RT Dose for the treatment of planning systems；
E102, using the dosage after Rapid Dose Calculation as standard, calculate γ values；
If e103, γ value are more than 1, obstructed overtreatment verification passes through dosage verifying if γ values are less than or equal to 1.
The γ is verified as another verification method, and for conventional verification method, the computational methods of wherein γ values are shown in document
(Low D A,Harms W B,Mutic S,et al.A technique for the quantitative evaluation
of dose distributions[J].Medical Physics,1998,25(5):656-661.), main is exactly with we
The calculated dosage of case is standard, is read with treatment planning systems RT Dose (file for storing dosage in DICOM file)
Dosage be compared, judge whether to be verified；Can be not only according to deviation Dev, but also dosage judged according to γ verifications
Accuracy ensures the accuracy of Afterloading radiotherapy, further increases safety, avoids that malpractice occurs.
Specific implementation mode of the invention above described embodiment only expresses, the description thereof is more specific and detailed, but simultaneously
Cannot the limitation to the scope of the claims of the present invention therefore be interpreted as.It should be pointed out that for those of ordinary skill in the art
For, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to the guarantor of the present invention
Protect range.
Claims (8)
1. a kind of calculating of quick afterloading radiotherapy dosage, verification method, which is characterized in that include the following steps：
A, the corresponding DICOM file of patient is exported from treatment planning systems, is saved in a glue file folder；
B, glue file folder is scanned, required Rapid Dose Calculation information is read from DICOM file；
C, radioactive source parameter is obtained, Distribution of dose rate table is established；
D, it determines radioactive source direction, and Rapid Dose Calculation is carried out according to Rapid Dose Calculation information, obtain Rapid Dose Calculation result；
E, the Dose Results in the Rapid Dose Calculation result and treatment planning systems of acquisition are compared, obtain deviation Dev and
γ verification results；
F, after the completion of comparison, deviation Dev and γ verification result is preserved to glue file in pressing from both sides, and be automatically deleted glue file
DICOM file in folder.
2. the calculating of quick afterloading radiotherapy dosage according to claim 1, verification method, which is characterized in that the step c
The middle data acquisition radioactive source parameter recommended according to AAPM and ESTRO.
3. the calculating of quick afterloading radiotherapy dosage according to claim 1 or 2, verification method, which is characterized in that the step
Distribution of dose rate table is established in rapid c to include the following steps：
Calculate the dosage rate at a certain Rapid Dose Calculation pointWherein,For
Dosage rate at Rapid Dose Calculation point, r are Rapid Dose Calculation point to the distance of radiation source center, r0=1cm, θ are agent in polar coordinate system
Amount calculates the angle between point and the directions radioactive source long axis V, θ 0=pi/2s, SkFor air kerma intensity, Λ is that dosage rate is normal
Amount, G is geometrical factor, and g is radial dose function, and F is anisotropy function；
According to the symmetry that dosage around radioactive source is distributed, calculate along the radioactive source directions long axis V and perpendicular to radioactive source long axis U
The two-dimensional Distribution of dose rate table T (m, n) in direction, and deposit in calculator memory.
4. the calculating of quick afterloading radiotherapy dosage according to claim 1, verification method, which is characterized in that the step d
Middle determining radioactive source direction includes the following steps：
Set i-th of parked position S of radioactive sourcei(xi, yi, zi), and pass through parked position Si(xi, yi, zi) and next stay
Stay position Si+1(xi+1, yi+1, zi+1) composition vectorDetermine direction of the radioactive source in human body coordinate system, and to
Amount
5. the calculating of quick afterloading radiotherapy dosage according to claim 1 or 4, verification method, which is characterized in that described
Before carrying out Rapid Dose Calculation in step d, it is P (x, y, z) that setting dosage, which calculates point, and applicating device pipeline sum is NA, every applicating device
Radioactive source parked position sum is N in pipelineS, i-th of parked position in jth root applicating device pipeline is to the agent at P (x, y, z)
Dose rate contribution is dJ, i, the residence time of i-th of parked position in jth root applicating device pipeline is tJ, i。
6. the calculating of quick afterloading radiotherapy dosage according to claim 5, verification method, which is characterized in that in the step
The Rapid Dose Calculation at P (x, y, z) is included the following steps in rapid d：
D101, P (x, y, z) i-th away from radioactive source parked position S is calculatedi(xi, yi, zi) distance r ', then
D102, according to vector
It is calculatedWithBetween angle
D103, i-th away from radioactive source parked position S of P (x, y, z) is calculated separatelyi(xi, yi, zi) directions distance v and U in the directions V
Distance u, wherein v=r ' cos θ,
D104, Distribution of dose rate table T (m, n) is checked, if m × 0.1cm≤u ＜ (m+1) × 0.1cm and n × 01cm≤v ＜ (n+
1) × 0.1cm, then to dJ, i(u, v) corresponding Distribution of dose rate carries out bilinear interpolation calculating, obtains dJ, i(u, v)=w1T (m,
n)+w2T (m, n+1)+w3T (m+1, n)+w4T (m+1, n+1),
Wherein, w1、w2、w3And w4Respectively in Distribution of dose rate table at T (m, n), T (m, n+1), T (m+1, n) and T (m+1, n+1)
Weight；
Dosage D (x, y, z) at d105, calculating P (x, y, z),
7. the calculating of quick afterloading radiotherapy dosage according to claim 6, verification method, which is characterized in that the step e
Middle calculating deviation Dev includes the following steps：
E101, dosimetric parameters D is obtained from the RT Dose for the treatment of planning systemsTPS；
Dosimetric parameters D is obtained after e102, Rapid Dose CalculationQA；
E103, according to dosimetric parameters DTPSWith dosimetric parameters DQADeviation Dev is calculated, is obtained
8. the calculating of quick afterloading radiotherapy dosage according to claim 6, verification method, which is characterized in that the step e
Middle γ, which is verified, to be included the following steps,：
E101, dosage is read from the RT Dose for the treatment of planning systems；
E102, using the dosage after Rapid Dose Calculation as standard, calculate γ values；
If e103, γ value are more than 1, obstructed overtreatment verification passes through dosage verifying if γ values are less than or equal to 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810591778.2A CN108671419B (en) | 2018-06-11 | 2018-06-11 | Calculation and verification method for radiotherapy dose after rapid loading |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810591778.2A CN108671419B (en) | 2018-06-11 | 2018-06-11 | Calculation and verification method for radiotherapy dose after rapid loading |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108671419A true CN108671419A (en) | 2018-10-19 |
| CN108671419B CN108671419B (en) | 2021-05-11 |
Family
ID=63810693
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810591778.2A Active CN108671419B (en) | 2018-06-11 | 2018-06-11 | Calculation and verification method for radiotherapy dose after rapid loading |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108671419B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113117253A (en) * | 2021-04-20 | 2021-07-16 | 中北大学 | Dose calculation system based on anisotropic kernel |
| CN115607861A (en) * | 2022-12-19 | 2023-01-17 | 安徽慧软科技有限公司 | Nuclear magnetic guided three-dimensional Monte Carlo dose independent verification method and system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060293583A1 (en) * | 2005-06-27 | 2006-12-28 | Saracen Michael J | Method for automatic anatomy-specific treatment planning protocols based on historical integration of previously accepted plans |
| CN105854191A (en) * | 2016-04-26 | 2016-08-17 | 中国科学院合肥物质科学研究院 | System and method for three-dimensional dose verification in radiosurgery |
-
2018
- 2018-06-11 CN CN201810591778.2A patent/CN108671419B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060293583A1 (en) * | 2005-06-27 | 2006-12-28 | Saracen Michael J | Method for automatic anatomy-specific treatment planning protocols based on historical integration of previously accepted plans |
| CN105854191A (en) * | 2016-04-26 | 2016-08-17 | 中国科学院合肥物质科学研究院 | System and method for three-dimensional dose verification in radiosurgery |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113117253A (en) * | 2021-04-20 | 2021-07-16 | 中北大学 | Dose calculation system based on anisotropic kernel |
| CN115607861A (en) * | 2022-12-19 | 2023-01-17 | 安徽慧软科技有限公司 | Nuclear magnetic guided three-dimensional Monte Carlo dose independent verification method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108671419B (en) | 2021-05-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104338240B (en) | Automatic optimization device for on-line self-adaption radiotherapy plan | |
| CN100431642C (en) | Method for determining radiating field output dose accurately in conformalradiotherapy | |
| US10238892B2 (en) | System and method for patient-specific radiotherapy treatment verification and quality assurance system | |
| CN109562277B (en) | Automated method and implementation system for radiation therapy dose calibration, reconstruction and verification | |
| CN105854191A (en) | System and method for three-dimensional dose verification in radiosurgery | |
| Zourari et al. | Dosimetric accuracy of a deterministic radiation transport based 192Ir brachytherapy treatment planning system. Part III. Comparison to Monte Carlo simulation in voxelized anatomical computational models | |
| CN109464756A (en) | Verify the method, apparatus and radiotherapy equipment of radiotherapy dosage | |
| Renner et al. | A dose delivery verification method for conventional and intensity‐modulated radiation therapy using measured field fluence distributions | |
| CN104888364A (en) | Ion beam radiotherapy dosage verification method based on water equivalent coefficients | |
| CN109147952A (en) | A kind of Radiotherapy dosimetry verification method, processor and system | |
| Kim et al. | Correlation analysis between 2D and quasi-3D gamma evaluations for both intensity-modulated radiation therapy and volumetric modulated arc therapy | |
| CN108671419A (en) | A kind of calculating of quick afterloading radiotherapy dosage, verification method | |
| CN106456993A (en) | Radiation-therapy simulation apparatus | |
| Esposito et al. | Improving dose delivery accuracy with EPID in vivo dosimetry: results from a multicenter study | |
| CN109125952A (en) | Convolution based on nuclear model is superimposed dose calculation methodology | |
| Ecclestone et al. | Experimental validation of the van Herk margin formula for lung radiation therapy | |
| CN115120891B (en) | Dose transmission evaluation device, computer-readable storage medium and system | |
| CN110368006A (en) | For providing the method and X-ray equipment of the transitional information about image data set | |
| CN116785601A (en) | Method and system for robust radiotherapy treatment planning with dose mapping uncertainty | |
| CN105797284A (en) | Radioactive particle radiation dose verification method and system | |
| CN114344735A (en) | Radiation dose monitoring and verifying method based on scattered ray characteristics | |
| CN109999365B (en) | Computing device and storage medium for implementing radiotherapy reverse optimization method | |
| CN110975174A (en) | Method, device and system for manufacturing dose verification die body model | |
| Zhong et al. | Generation of a novel phase‐space‐based cylindrical dose kernel for IMRT optimization | |
| CN111921098B (en) | Robustness optimization system based on radiotherapy biological effect |
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 |