CN112700841A - Dose calculation modeling method, model, device and storage medium in non-uniform mode - Google Patents
Dose calculation modeling method, model, device and storage medium in non-uniform mode Download PDFInfo
- Publication number
- CN112700841A CN112700841A CN202011644390.8A CN202011644390A CN112700841A CN 112700841 A CN112700841 A CN 112700841A CN 202011644390 A CN202011644390 A CN 202011644390A CN 112700841 A CN112700841 A CN 112700841A
- Authority
- CN
- China
- Prior art keywords
- dose
- distribution
- kernel
- calculation
- data
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000003860 storage Methods 0.000 title claims abstract description 18
- 238000009826 distribution Methods 0.000 claims abstract description 69
- 230000005855 radiation Effects 0.000 claims abstract description 24
- 230000006870 function Effects 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 10
- 238000011156 evaluation Methods 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000000329 molecular dynamics simulation Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 3
- 230000009878 intermolecular interaction Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000001959 radiotherapy Methods 0.000 description 6
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000005623 Carcinogenesis Diseases 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000036952 cancer formation Effects 0.000 description 1
- 231100000504 carcinogenesis Toxicity 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000010977 jade Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Surgery (AREA)
- Urology & Nephrology (AREA)
- Engineering & Computer Science (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
The embodiment of the application discloses a dose calculation modeling method, a model, equipment and a storage medium in a non-uniform f (FFF) mode. The dose calculation modeling method under the non-uniform mode mainly comprises the following steps: acquiring two-dimensional incident intensity distribution; acquiring an irradiation field dose nucleus database comprising a plurality of dose nuclei under different irradiation fields; calculating dose distributions under different radiation fields based on the two-dimensional incident intensity distribution and the radiation field dose nuclear database; and respectively fitting the data of the dose distribution under each different field, recalculating the dose distribution after continuously changing the dose kernel distribution, and comparing the recalculated dose distribution with the actually measured dose distribution to obtain an optimal dose kernel, thereby forming final dose kernel model data. The dose calculation model in the non-uniform (FFF) mode of the embodiment of the invention can quickly realize dose calculation and simultaneously improve the calculation precision.
Description
Technical Field
The invention relates to the technical field related to medical tumor radiotherapy, in particular to a dose calculation modeling method, a model, equipment and a storage medium under a non-uniform f (FFF) mode.
Background
Malignant tumors are the first killer of human beings, and the current tumor treatment mode is important in the aspects of radiotherapy, chemotherapy and surgery, and radiation therapy is one of the most common methods in tumor treatment at present.
The conventional radiotherapy requires the dose distribution in the radiation field to be equalized, and the mode is favorable for dose calculation and control in the treatment process, and has the defects that the dose rate is greatly reduced after the equalization of an equalizing block compared with that before the equalization, and the direct result is that the treatment time is too long. The treatment time is an important factor for evaluating the radiotherapy surgical technology, and the excessive treatment time can increase the physical burden of the patient, damage the normal tissue cells of the patient and generate the possibility of secondary carcinogenesis.
Dose calculation is the core content of a radiation treatment planning system, and the accuracy of dose distribution mainly depends on a dose calculation model adopted in the planning system. At present, the dose algorithm built in the radiotherapy planning system is mostly based on the pencil beam convolution superposition technology, and the monte carlo algorithm obtains the distribution of the deposition energy of particles in human tissues by randomly simulating the interaction between examples and substances, and is considered as the most accurate calculation method in all the current calculation methods. At present, the Monte Carlo algorithm is difficult to be applied to the reverse TMRT plan as long as the problem of how to solve the contradiction between the calculation precision and the calculation speed is high. On the premise of keeping the characteristic of high precision, how to accelerate the calculation speed is the main subject of the Monte Carlo dose calculation method, and the development of the dose calculation method with higher speed and higher precision is a future research hotspot.
Disclosure of Invention
In view of the deficiencies of the prior art, one aspect of the present invention provides a dose calculation modeling method in a non-uniform mode.
The dose calculation modeling method under the non-uniform mode mainly comprises the following steps:
acquiring two-dimensional incident intensity distribution;
acquiring an irradiation field dose nucleus database comprising a plurality of dose nuclei under different irradiation fields;
calculating dose distributions under different radiation fields based on the two-dimensional incident intensity distribution and the radiation field dose nuclear database;
and respectively fitting the data of the dose distribution under each different field, recalculating the dose distribution after continuously changing the dose kernel distribution, and comparing the recalculated dose distribution with the actually measured dose distribution to obtain an optimal dose kernel, thereby forming final dose kernel model data.
According to a preferred embodiment of the present invention, the two-dimensional incident intensity distribution is obtained by measuring the dose distribution of the maximum field XOY plane of the accelerator when in use according to the intensity of different positions when the rays enter the body surface.
According to a preferred embodiment of the invention, the radiation field dose kernel database comprising several different under-field dose kernels is composed of dose kernels K (F, z) of different depths, the dose kernel being composed of a dose of n volume units, i.e. K0、K1、K2......KnThe distance between the volume unit and the center of the dose core is r0、r1、r2......rn。
According to a preferred embodiment of the present invention, calculating the dose distribution under different radiation fields based on the two-dimensional incident intensity distribution and the radiation field dose nuclear database comprises the following steps:
according to the setting of the irradiation field, finding out a dose kernel parameter corresponding to the irradiation field in the database model, and using a dose calculation formula of the parameter in the homogeneous phantom as follows:
in the formula (1), the first and second groups,
Ki(F, z) is the dose nuclear parameter of the irradiation field F at depth z in the database model, i is the dose nuclear product unit index, r isiRepresenting the distance of the ith volume element from the center of the dose kernel;
i (x, y) represents the intensity of a point with coordinates (x, y) in the plane XOY;
order to
I.e. all distances from point (x, y) are riThe sum of the intensities of the volume units of (1), current point(x+ricosθ，y+risin θ) in the irradiated field, I (x + r)icosθ，y+risin θ) is 0, the formula (1) is simplified to the following formula:
according to a preferred embodiment of the present invention, fitting the data of the dose distribution under each of the different fields, recalculating the dose distribution after continuously changing the dose kernel distribution and comparing the recalculated dose distribution with the actually measured dose distribution to obtain the optimal dose kernel, so as to form the final dose kernel model data includes the following steps:
s401: simulating the dose distribution of the pencil beams by using a Monte Carlo algorithm, taking the distribution as the initial value of dose cores of all fields, and respectively fitting according to off-axis ratio data under different fields;
s402: calculating the dosage d (x) of the corresponding point on the coordinate axis by using a formula (2) according to off-axis ratio data of a certain depth in the measurement radiation field, wherein,
s403: comparing the calculated data with the measured off-axis ratio data using the following evaluation function:
in the formula (4), the first and second groups,
m represents the number of actually measured off-axis ratio data points;
wjrepresenting the weight of the jth point of the off-axis ratio data;
djrepresenting the dose value calculated by formula (3) at the point corresponding to the measured off-axis ratio data;
xjrepresenting the off-axis distance from the jth point of the axial ratio data;
s404: dose kernel parameter K according to equation (4)jCalculating partial derivatives to obtain formula (5):
s405: equation (5) is analogous to a molecular dynamics equation and uses fjExpressing the intermolecular interaction force yields equation (6):
at this time, K is addedk(F, z) is analogous to the position of the kth atom in molecular dynamics, the atomic weight of this atom being in mkExpressed, equation (7) is obtained:
s406: adopting a formula (8) and a formula (9) to carry out iteration, introducing a temporary time variable t in the iteration process, obtaining different dose kernel parameters by continuously changing the value of t,
s407: recalculating the evaluation function of the dose kernel parameter obtained in the step S406 by using a formula (4), and quitting the iteration process when the evaluation function is smaller than a set value or the iteration frequency reaches the set value to obtain a final dose kernel parameter;
s408: obtaining dose kernel data of each measured depth through steps S402-S407, and obtaining the dose kernel data of the depth which is not measured through interpolation to form a dose kernel model;
s409: and finally, multiplying the dose kernel data by the correction coefficient of each depth to obtain final dose kernel model data.
Another aspect of the invention provides a model for dose calculation in non-uniform mode.
The dose calculation model in the non-uniform mode is formed by any one of the dose calculation modeling methods in the non-uniform mode.
Another aspect of the invention provides a dose calculation method in non-uniform mode.
The dose calculation method in the non-uniform mode is realized through the dose calculation model in the non-uniform mode.
Another aspect of the invention provides a dose calculation device in non-leveling mode.
The dose calculation device in the non-leveling mode comprises a processor and a memory; the memory is used for storing instructions, and the instructions, when executed by the processor, cause the dose calculation device in the non-uniform mode to implement the dose calculation modeling method in the non-uniform mode as described in any one of the above or the corresponding operation of the dose calculation method in the non-uniform mode as described above.
Yet another aspect of the invention provides a computer-readable storage medium.
The storage medium stores computer instructions, and after the computer reads the computer instructions in the storage medium, the computer executes the dose calculation modeling method in the non-uniform mode or the dose calculation method in the non-uniform mode.
Compared with the prior art, the dose calculation modeling method, the model, the device and the storage medium in the non-uniform mode have the following beneficial effects:
the dose calculation model in the non-uniform (FFF) mode of the embodiment of the invention can quickly realize dose calculation and simultaneously improve the calculation precision.
Additional features of the invention will be set forth in part in the description which follows. Additional features of some aspects of the invention will become apparent to those of ordinary skill in the art upon examination of the following description and accompanying drawings or may be learned by the manufacture or operation of the embodiments. The features of the present disclosure may be realized and attained by practice or use of various methods, instrumentalities and combinations of the specific embodiments described below.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention without limiting the invention. Like reference symbols in the various drawings indicate like elements. Wherein,
FIG. 1 is a schematic illustration of a dose kernel distribution in an embodiment of the present invention.
Detailed Description
The following description is presented to enable one of ordinary skill in the art to make and use the application and is provided in the context of a particular application and its requirements. Various modifications to the disclosed embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be applied to other embodiments and applications without departing from the spirit and scope of the present application. Thus, the present application is not limited to the embodiments shown, but is to be accorded the widest scope consistent with the claims.
The terminology used herein is for the purpose of describing particular example embodiments only and is not intended to limit the scope of the present application. As used herein, the singular forms "a", "an" and "the" may include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
The features and characteristics of the present application, as well as the methods of operation and functions of the related elements of structure and the combination of parts and economies of manufacture, will become more apparent upon consideration of the following description of the drawings, which form a part hereof. It is to be understood, however, that the drawings are designed solely for the purposes of illustration and description and are not intended as a definition of the limits of the application. It should be understood that the drawings are not to scale.
Flow charts are used herein to illustrate operations performed by systems according to embodiments of the present application. It should be understood that the operations in the flow diagrams are not necessarily performed exactly in order. Rather, various steps may be processed in reverse order or simultaneously. Further, one or more other operations may be added to the flowchart. One or more operations may also be deleted from the flowcharts.
One aspect of the embodiments of the present invention discloses a dose calculation modeling method in a non-uniform mode.
The dose calculation modeling method under the non-uniform mode mainly comprises the following steps:
a two-dimensional incident intensity distribution is obtained.
The two-dimensional incident intensity distribution is obtained by measuring the dose distribution of the maximum field XOY plane of the accelerator when the accelerator is used according to the intensity of different positions when rays enter a body surface.
An irradiation field dose kernel database containing a plurality of dose kernels under different irradiation fields is obtained.
Wherein the radiation field dose kernel database comprising several different radiation field dose kernels consists of dose kernels K (F, z) of different depths, the dose kernel consisting of a dose of n volume units, i.e. K0、K1、K2......KnThe distance between the volume unit and the center of the dose core is r0、r1、r2......rnAs shown in fig. 1.
And calculating the dose distribution under different radiation fields based on the two-dimensional incident intensity distribution and the radiation field dose nuclear database.
Calculating the dose distribution under different radiation fields based on the two-dimensional incident intensity distribution and the radiation field dose nuclear database comprises the following steps:
according to the setting of the irradiation field, finding out a dose kernel parameter corresponding to the irradiation field in the database model, and using a dose calculation formula of the parameter in the homogeneous phantom as follows:
in the formula (1), the first and second groups,
Ki(F, z) is the dose nuclear parameter of the irradiation field F at depth z in the database model, i is the dose nuclear product unit index, r isiRepresenting the distance of the ith volume element from the center of the dose kernel;
i (x, y) represents the intensity of a point with coordinates (x, y) in the plane XOY;
order to
I.e. all distances from point (x, y) are riThe sum of intensities of the volume units of (c), when point (x + r)icosθ，y+risin θ) in the irradiated field, I (x + r)icosθ，y+risin θ) is 0, the formula (1) is simplified to the following formula:
and respectively fitting the data of the dose distribution under each different field, recalculating the dose distribution after continuously changing the dose kernel distribution, and comparing the recalculated dose distribution with the actually measured dose distribution to obtain an optimal dose kernel, thereby forming final dose kernel model data.
The method comprises the following steps of respectively fitting the data of dose distribution under each different field, recalculating the dose distribution after continuously changing the dose kernel distribution, and comparing the recalculated dose distribution with the actually measured dose distribution to obtain an optimal dose kernel, so as to form final dose kernel model data, wherein the method comprises the following steps:
s401: simulating the dose distribution of the pencil beams by using a Monte Carlo algorithm, taking the distribution as the initial value of dose cores of all fields, and respectively fitting according to off-axis ratio data under different fields;
s402: calculating the dosage d (x) of the corresponding point on the coordinate axis by using a formula (2) according to off-axis ratio data of a certain depth in the measurement radiation field, wherein,
s403: comparing the calculated data with the measured off-axis ratio data using the following evaluation function:
in the formula (4), the first and second groups,
m represents the number of actually measured off-axis ratio data points;
wjrepresenting the weight of the jth point of the off-axis ratio data;
djrepresenting the dose value calculated by formula (3) at the point corresponding to the measured off-axis ratio data;
xjrepresenting the off-axis distance from the jth point of the axial ratio data.
S404: dose kernel parameter K according to equation (4)jCalculating partial derivatives to obtain formula (5):
s405: equation (5) is analogous to a molecular dynamics equation and uses fjExpressing the intermolecular interaction force yields equation (6):
at this time, K is addedk(F, z) is analogous to the position of the kth atom in molecular dynamics, the atomic weight of this atom being in mkExpressed, equation (7) is obtained:
s406: adopting a formula (8) and a formula (9) to carry out iteration, introducing a temporary time variable t in the iteration process, obtaining different dose kernel parameters by continuously changing the value of t,
s407: recalculating the evaluation function of the dose kernel parameter obtained in the step S406 by using a formula (4), and quitting the iteration process when the evaluation function is smaller than a set value or the iteration frequency reaches the set value to obtain a final dose kernel parameter;
s408: obtaining dose kernel data of each measured depth through steps S402-S407, and obtaining the dose kernel data of the depth which is not measured through interpolation to form a dose kernel model;
s409: and finally, multiplying the dose kernel data by the correction coefficient of each depth to obtain final dose kernel model data.
Another aspect of an embodiment of the present invention discloses a dose calculation model in a non-averaging mode.
The dose calculation model in the non-uniform mode is formed by the dose calculation modeling method in the non-uniform mode.
Another aspect of the embodiments of the present invention discloses a dose calculation method in a non-uniform mode.
The dose calculation method in the non-uniform mode is realized through the dose calculation model in the non-uniform mode.
Another aspect of an embodiment of the invention discloses a dose calculation device in a non-leveling mode.
The dose calculation device in the non-leveling mode comprises a processor and a memory; the memory is configured to store instructions, and the instructions, when executed by the processor, cause the dose calculation apparatus in the non-uniform mode to implement any one of the dose calculation modeling method in the non-uniform mode or an operation corresponding to the dose calculation method in the non-uniform mode.
Another aspect of an embodiment of the present invention discloses a computer-readable storage medium. The storage medium stores computer instructions, and after the computer reads the computer instructions in the storage medium, the computer runs the dose calculation modeling method in the non-uniform mode or the dose calculation method in the non-uniform mode.
The dose calculation modeling method, model, device and computer-readable storage medium in the non-uniform mode according to the embodiments of the present invention can bring beneficial effects including, but not limited to, the following:
the dose calculation model in the non-uniform (FFF) mode of the embodiment of the invention can quickly realize dose calculation and simultaneously improve the calculation precision.
Having thus described the basic concept, it will be apparent to those skilled in the art that the foregoing detailed disclosure is to be considered merely illustrative and not restrictive of the broad application. Various modifications, improvements and adaptations to the present application may occur to those skilled in the art, although not explicitly described herein. Such modifications, improvements and adaptations are proposed in the present application and thus fall within the spirit and scope of the exemplary embodiments of the present application.
Also, this application uses specific language to describe embodiments of the application. Reference throughout this specification to "one embodiment," "an embodiment," and/or "some embodiments" means that a particular feature, structure, or characteristic described in connection with at least one embodiment of the present application is included in at least one embodiment of the present application. Therefore, it is emphasized and should be appreciated that two or more references to "an embodiment" or "one embodiment" or "an alternative embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment. Furthermore, some features, structures, or characteristics of one or more embodiments of the present application may be combined as appropriate.
Moreover, those skilled in the art will appreciate that aspects of the present application may be illustrated and described in terms of several patentable species or situations, including any new and useful combination of processes, machines, manufacture, or materials, or any new and useful improvement thereon. Accordingly, various aspects of the present application may be embodied entirely in hardware, entirely in software (including firmware, resident software, micro-code, etc.) or in a combination of hardware and software. The above hardware or software may be referred to as "data block," module, "" engine, "" unit, "" component, "or" system. Furthermore, aspects of the present application may be represented as a computer product, including computer readable program code, embodied in one or more computer readable media.
The computer storage medium may comprise a propagated data signal with the computer program code embodied therewith, for example, on baseband or as part of a carrier wave. The propagated signal may take any of a variety of forms, including electromagnetic, optical, etc., or any suitable combination. A computer storage medium may be any computer-readable medium that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code located on a computer storage medium may be propagated over any suitable medium, including radio, cable, fiber optic cable, RF, or the like, or any combination of the preceding.
Computer program code required for the operation of various portions of the present application may be written in any one or more programming languages, including an object oriented programming language such as Java, Scala, Smalltalk, Eiffel, JADE, Emerald, C + +, C #, VB.NET, Python, and the like, a conventional programming language such as C, Visualbasic, Fortran2003, Perl, COBOL2002, PHP, ABAP, a dynamic programming language such as Python, Ruby, and Groovy, or other programming languages, and the like. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user's computer through any network format, such as a Local Area Network (LAN) or a Wide Area Network (WAN), or the connection may be made to an external computer (for example, through the Internet), or in a cloud computing environment, or as a service, such as a software as a service (SaaS).
Additionally, the order in which elements and sequences of the processes described herein are processed, the use of alphanumeric characters, or the use of other designations, is not intended to limit the order of the processes and methods described herein, unless explicitly claimed. While various presently contemplated embodiments have been discussed in the foregoing disclosure by way of example, it should be understood that such detail is solely for that purpose and that the appended claims are not limited to the disclosed embodiments, but, on the contrary, are intended to cover all modifications and equivalent arrangements that are within the spirit and scope of the embodiments herein. For example, although the system components described above may be implemented by hardware devices, they may also be implemented by software-only solutions, such as installing the described system on an existing server or mobile device.
Similarly, it should be noted that in the preceding description of embodiments of the application, various features are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the embodiments. This method of disclosure, however, is not intended to require more features than are expressly recited in the claims. Indeed, the embodiments may be characterized as having less than all of the features of a single embodiment disclosed above.
Claims (9)
1. A dose calculation modeling method in a non-uniform mode is characterized by comprising the following steps:
acquiring two-dimensional incident intensity distribution;
acquiring an irradiation field dose nucleus database comprising a plurality of dose nuclei under different irradiation fields;
calculating dose distributions under different radiation fields based on the two-dimensional incident intensity distribution and the radiation field dose nuclear database;
and respectively fitting the data of the dose distribution under each different field, recalculating the dose distribution after continuously changing the dose kernel distribution, and comparing the recalculated dose distribution with the actually measured dose distribution to obtain an optimal dose kernel, thereby forming final dose kernel model data.
2. The dose calculation modeling method in non-uniform mode according to claim 1,
the two-dimensional incident intensity distribution is obtained by measuring the dose distribution of the maximum field XOY plane of the accelerator when the accelerator is used according to the intensity of different positions when rays enter a body surface.
3. The dose calculation modeling method in non-uniform mode according to claim 1,
the radiation field dose kernel database comprising several dose kernels under different radiation fields consists of dose kernels K (F, z) of different depths, the dose kernel consisting of a dose of n volume units, i.e. K0、K1、K2......KnThe distance between the volume unit and the center of the dose core is r0、r1、r2......rn。
4. The dose calculation modeling method in non-uniform mode according to claim 1,
calculating the dose distribution under different radiation fields based on the two-dimensional incident intensity distribution and the radiation field dose nuclear database comprises the following steps:
according to the setting of the irradiation field, finding out a dose kernel parameter corresponding to the irradiation field in the database model, and using a dose calculation formula of the parameter in the homogeneous phantom as follows:
in the formula (1), the first and second groups,
Ki(F, z) is the dose nuclear parameter of the irradiation field F at depth z in the database model, i is the dose nuclear product unit index, r isiRepresenting the distance of the ith volume element from the center of the dose kernel;
i (x, y) represents the intensity of a point with coordinates (x, y) in the plane XOY;
order to
I.e. all distances from point (x, y) are riThe sum of intensities of the volume units of (c), when point (x + r)icosθ，y+risin θ) in the irradiated field, I (x + r)icosθ，y+risin θ) is 0, the formula (1) is simplified to the following formula:
5. the dose calculation modeling method in non-uniform mode according to claim 4,
fitting the data of the dose distribution under each different field respectively, recalculating the dose distribution after continuously changing the dose kernel distribution and comparing the recalculated dose distribution with the actually measured dose distribution to obtain the optimal dose kernel, thereby forming the final dose kernel model data, comprising the following steps:
s401: simulating the dose distribution of the pencil beams by using a Monte Carlo algorithm, taking the distribution as the initial value of dose cores of all fields, and respectively fitting according to off-axis ratio data under different fields;
s402: calculating the dosage d (x) of the corresponding point on the coordinate axis by using a formula (2) according to off-axis ratio data of a certain depth in the measurement radiation field, wherein,
s403: comparing the calculated data with the measured off-axis ratio data using the following evaluation function:
in the formula (4), the first and second groups,
m represents the number of actually measured off-axis ratio data points;
wjrepresenting the weight of the jth point of the off-axis ratio data;
djrepresenting the dose value calculated by formula (3) at the point corresponding to the measured off-axis ratio data;
xjrepresenting the off-axis distance from the jth point of the axial ratio data;
s404: dose kernel parameter K according to equation (4)jCalculating partial derivatives to obtain formula (5):
s405: equation (5) is analogous to a molecular dynamics equation and usedfjExpressing the intermolecular interaction force yields equation (6):
at this time, K is addedk(F, z) is analogous to the position of the kth atom in molecular dynamics, the atomic weight of this atom being in mkExpressed, equation (7) is obtained:
s406: adopting a formula (8) and a formula (9) to carry out iteration, introducing a temporary time variable t in the iteration process, obtaining different dose kernel parameters by continuously changing the value of t,
s407: recalculating the evaluation function of the dose kernel parameter obtained in the step S406 by using a formula (4), and quitting the iteration process when the evaluation function is smaller than a set value or the iteration frequency reaches the set value to obtain a final dose kernel parameter;
s408: obtaining dose kernel data of each measured depth through steps S402-S407, and obtaining the dose kernel data of the depth which is not measured through interpolation to form a dose kernel model;
s409: and finally, multiplying the dose kernel data by the correction coefficient of each depth to obtain final dose kernel model data.
6. A non-uniform mode dose calculation model formed by the non-uniform mode dose calculation modeling method according to any one of claims 1 to 5.
7. A dose calculation method in non-uniform mode, which is implemented by the dose calculation model in non-uniform mode according to claim 6.
8. A dose calculation device in non-uniform mode is characterized by comprising a processor and a memory; the memory is used for storing instructions, and the instructions, when executed by the processor, cause the dose calculation device in the non-leveling mode to implement the dose calculation modeling method in the non-leveling mode according to any one of claims 1 to 5 or the corresponding operation of the dose calculation method in the non-leveling mode according to claim 7.
9. A computer-readable storage medium, wherein the storage medium stores computer instructions, and when the computer instructions in the storage medium are read by a computer, the computer executes the modeling method for dose calculation in non-uniform mode according to any one of claims 1 to 5 or the dose calculation method in non-uniform mode according to claim 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011644390.8A CN112700841B (en) | 2020-12-24 | 2020-12-24 | Dose calculation modeling method, model, device and storage medium in non-uniform mode |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011644390.8A CN112700841B (en) | 2020-12-24 | 2020-12-24 | Dose calculation modeling method, model, device and storage medium in non-uniform mode |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112700841A true CN112700841A (en) | 2021-04-23 |
| CN112700841B CN112700841B (en) | 2023-03-28 |
Family
ID=75514312
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011644390.8A Active CN112700841B (en) | 2020-12-24 | 2020-12-24 | Dose calculation modeling method, model, device and storage medium in non-uniform mode |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112700841B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160059040A1 (en) * | 2014-08-29 | 2016-03-03 | Wisconsin Alumni Research Foundation | System and method for control of external beam radiation |
| CN110310720A (en) * | 2018-03-20 | 2019-10-08 | 北京连心医疗科技有限公司 | A kind of Monte Carto dosage computing method, equipment and storage medium |
| CN110504016A (en) * | 2018-05-18 | 2019-11-26 | 北京连心医疗科技有限公司 | A kind of Monte Carlo grid parallel computation dose calculation methodology, equipment and storage medium |
| CN111408073A (en) * | 2020-03-24 | 2020-07-14 | 宁波同调医学科技有限公司 | Method and system for calculating dose of radiotherapy ray planar detector |
-
2020
- 2020-12-24 CN CN202011644390.8A patent/CN112700841B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160059040A1 (en) * | 2014-08-29 | 2016-03-03 | Wisconsin Alumni Research Foundation | System and method for control of external beam radiation |
| CN110310720A (en) * | 2018-03-20 | 2019-10-08 | 北京连心医疗科技有限公司 | A kind of Monte Carto dosage computing method, equipment and storage medium |
| CN110504016A (en) * | 2018-05-18 | 2019-11-26 | 北京连心医疗科技有限公司 | A kind of Monte Carlo grid parallel computation dose calculation methodology, equipment and storage medium |
| CN111408073A (en) * | 2020-03-24 | 2020-07-14 | 宁波同调医学科技有限公司 | Method and system for calculating dose of radiotherapy ray planar detector |
Non-Patent Citations (3)
| Title |
|---|
| 王敏杰: "非均整（FFF）模式下脑转移癌分次立体定向放射治疗的剂量学研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 * |
| 陈昭等: "Geant4中不同物理列表对硼中子俘获治疗剂量分布的影响", 《中国医学物理学杂志》 * |
| 龙凤翔等: "GZP3钴-60近距离治疗源剂量学参数的蒙特卡洛计算", 《中国医学物理学杂志》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112700841B (en) | 2023-03-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3730514B2 (en) | System and method for radiation dose calculation | |
| US11020615B2 (en) | Computing radiotherapy dose distribution | |
| US20210339049A1 (en) | Radiotherapy treatment plans using differentiable dose functions | |
| Baltas et al. | The physics of modern brachytherapy for oncology | |
| EP2760543B1 (en) | Dose computation for radiation therapy using heterogeneity compensated superposition | |
| EP3043862B1 (en) | Method for patient-specific radiotherapy treatment verification and quality assurance | |
| US20060259282A1 (en) | Deterministic computation of radiation transport for radiotherapy dose calculations and scatter correction for image reconstruction | |
| CN109464756A (en) | Verify the method, apparatus and radiotherapy equipment of radiotherapy dosage | |
| CN108415058A (en) | The dose calculation methodology and system of radioactive ray | |
| WO2005052721A2 (en) | Deterministic computation of radiation doses delivered to tissues and organs of a living organism | |
| CN104888364A (en) | Ion beam radiotherapy dosage verification method based on water equivalent coefficients | |
| CN110237445B (en) | EPID-based in-vivo three-dimensional dose monitoring and verifying method | |
| CN107551411B (en) | Proton heavy ion intensity modulated radiotherapy robust optimization method aiming at range uncertainty | |
| Toscano et al. | Impact of machine log-files uncertainties on the quality assurance of proton pencil beam scanning treatment delivery | |
| Kim et al. | Correlation analysis between 2D and quasi-3D gamma evaluations for both intensity-modulated radiation therapy and volumetric modulated arc therapy | |
| Gopishankar et al. | MAGAT gel and EBT2 film‐based dosimetry for evaluating source plugging‐based treatment plan in Gamma Knife stereotactic radiosurgery | |
| CN110310720A (en) | A kind of Monte Carto dosage computing method, equipment and storage medium | |
| CN113192554A (en) | Dose optimization method and medium for radiation particle internal radiotherapy interventional operation robot | |
| Li et al. | Extension and validation of a GPU‐Monte Carlo dose engine gDPM for 1.5 T MR‐LINAC online independent dose verification | |
| CN111494815B (en) | Three-dimensional dose calculation method, device and medium based on mixed variable-scale model | |
| US20160213948A1 (en) | System and method to create test objects for radiation therapy | |
| CN112700841B (en) | Dose calculation modeling method, model, device and storage medium in non-uniform mode | |
| Nakaguchi et al. | Validation of fluence‐based 3D IMRT dose reconstruction on a heterogeneous anthropomorphic phantom using Monte Carlo simulation | |
| CN110310743B (en) | Monte Carlo dose calculation method, equipment and storage medium | |
| Briscoe et al. | Establishing a threshold for rotational patient setup errors in linear accelerator-based stereotactic radiosurgery |
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 |