CN110675932A - Monte Carlo-based point dose calculation method, equipment and storage medium - Google Patents
Monte Carlo-based point dose calculation method, equipment and storage medium Download PDFInfo
- Publication number
- CN110675932A CN110675932A CN201810710921.5A CN201810710921A CN110675932A CN 110675932 A CN110675932 A CN 110675932A CN 201810710921 A CN201810710921 A CN 201810710921A CN 110675932 A CN110675932 A CN 110675932A
- Authority
- CN
- China
- Prior art keywords
- monte carlo
- dose calculation
- calculation
- dose
- dimensional
- 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.)
- Pending
Links
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/10—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
- G16H20/17—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients delivered via infusion or injection
-
- 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
Abstract
The invention belongs to the field of dose calculation, and relates to a dot dose calculation method based on Monte Carlo, a device and a storage medium. The method comprises the following steps: selecting a dose calculation point in a three-dimensional image of a patient or phantom; determining a calculation area; carrying out three-dimensional gridding on the calculation area, wherein each grid is a voxel; dividing any section in the incident direction of the beam into a two-dimensional flux grid; dose calculation is performed using a two-dimensional flux grid based on a monte carlo particle transport model. The invention adopts a point dose calculation method based on Monte Carlo, determines a calculation area of point dose by an effective electronic range, can fully consider the influence of particle scattering on the dose, reduces unnecessary dose calculation tasks while ensuring accurate calculation, improves the calculation speed and is beneficial to clinical application.
Description
Technical Field
The invention belongs to the field of dose calculation, and relates to a dot dose calculation method based on Monte Carlo, a device and a storage medium.
Background
Radiotherapy, chemotherapy and surgical operation are combined as three major means for treating tumors. In the existing method for treating cancer by using a linear accelerator, a treatment plan is generally made for a tumor region of each patient by using the TPS, and the quality of the plan is influenced by various factors, such as the level of planning made by a physicist, the performance of a TPS optimization algorithm, whether dose calculation can be accurately simulated, and the like. The dose calculation precision directly influences the iteration direction of the optimization algorithm and the judgment of a physicist, and becomes a basic index for measuring the quality of the TPS.
The Monte Carlo (MC) dose calculation method is the most accurate dose calculation method widely accepted in the industry, and is a random sampling simulation method which can be used for any absorption medium and is suitable for any ray, and the effect of a large number of single photons on a substance in a transportation process is simulated by a statistical method. After a photon enters a certain medium (e.g., human tissue), some or all of the energy of the ionizing radiation is transferred by interaction with atoms in the medium. The process of the action is simulated by a Monte Carlo method, and each action is predicted by applying the law of radiation physics and the result of the action is counted. The Monte Carlo simulation method is a three-dimensional dose calculation method with the strongest applicability, can be used for calculating dose distribution under various complex conditions, and is called as a 'gold standard' of dose calculation. However, at present, due to the statistical averaging of a large number of identical processes in Monte Care simulation, a large amount of time cost must be sacrificed while ensuring accuracy. Therefore, the Monte Carlo method has long time for calculating the dose (generally several to dozens of hours), and cannot meet the requirement of clinical real-time performance, so that the Monte Carlo method is not widely applied to a clinical treatment planning system.
In the prior art, the calculation of the deposited dose at the point dose position is large in calculation area, so that the dose calculation task based on the Monte Carlo particle transport model is increased, the calculation speed is reduced, and the clinical practical application is not facilitated.
Disclosure of Invention
It is an object of the present invention to provide a monte carlo-based spot dose calculation method, apparatus and storage medium that overcome the disadvantages of the prior art.
In order to achieve the purpose, the invention adopts the following technical scheme:
a monte carlo-based spot dose calculation method adapted to be executed in a computing device, characterized by: the method comprises the following steps:
selecting a dose calculation point in a three-dimensional image of a patient or phantom;
determining a calculation area;
performing three-dimensional gridding on a three-dimensional image of a patient or a phantom, wherein each grid is a voxel;
dividing any section in the incident direction of the beam into a two-dimensional flux grid; dose calculation is performed using a two-dimensional flux grid based on a monte carlo particle transport model.
Further preferably, the dose calculation point is located in a region of interest.
The method for determining the region of interest of the user comprises the following steps: obtained by user-defined or by calculation based on one or a combination of physical factors and biomedical factors.
The physical factors reflect the material composition and the irradiation physical conditions of the patient or the phantom; wherein the material composition of the patient or the phantom comprises: density, CT value, mass number and atomic number of the die body; the irradiation physical conditions include: field distribution and source distribution;
the biomedical factors comprise: organ tissue irradiation threshold, biological sensitivity, secondary cancer occurrence probability distribution, tumor killing probability distribution, and injury probability.
The method for determining the calculation area comprises the following steps:
(1) determining a sphere by taking the position of the point dose as a sphere center and taking the distance of an effective electron range as a radius;
(2) the portion of the sphere below the surface of the phantom is the calculated area.
The two-dimensional flux grid is a uniform two-dimensional grid.
Preferably, the three-dimensional voxels are uniform in size.
The present invention also provides a computing device comprising:
one or more processors;
a memory; and
one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for a Monte Carlo-based spot dose calculation method.
A computer readable storage medium storing one or more programs, the one or more programs comprising instructions adapted to be loaded from a memory and to carry out the monte carlo-based spot dose calculation method described above.
The invention has the following beneficial effects:
the invention adopts a point dose calculation method based on Monte Carlo, determines a calculation area of point dose by an effective electronic range, can fully consider the influence of particle scattering on the dose, reduces unnecessary dose calculation tasks while ensuring accurate calculation, improves the calculation speed and is beneficial to clinical application.
Drawings
FIG. 1 is a schematic diagram of a calculation region in a preferred embodiment of the present invention.
FIG. 1(a) the calculation region is a sphere;
the calculated area in fig. 1(b) is the spherical cap.
FIG. 2 is a flow chart of a Monte Carlo-based spot dose calculation method in a preferred embodiment of the present invention.
Detailed Description
The invention is further illustrated below with reference to the figures and examples.
Example 1
A monte carlo-based spot dose calculation method, adapted to be executed in a computing device, comprising the steps of (as shown in fig. 2):
step 210, selecting a dose calculation point in a three-dimensional image of a patient or a phantom;
wherein the dose calculation point is located in the region of interest. In this embodiment, preferably, the method for determining the region of interest of the user includes: obtained by user-defined or by calculation based on one or a combination of physical factors and biomedical factors.
Further preferably, the physical factor reflects the material composition of the patient or the phantom and the irradiation physical condition; wherein the material composition of the patient or the phantom comprises: density, CT value, mass number and atomic number of the die body; the irradiation physical conditions include: field distribution, source distribution, etc.
Further preferably, the biomedical factors include organ tissue irradiation threshold, biological sensitivity, secondary cancer occurrence probability distribution, tumor killing probability distribution, and damage probability.
fig. 1 is a schematic diagram of a calculation region in a preferred embodiment of the present invention, wherein a determination method of the calculation region includes the following steps:
(1) determining a sphere by taking the position of the point dose as a sphere center and taking the distance of an effective electron range as a radius;
(2) the portion of the sphere below the surface of the phantom is the calculated area.
In the embodiment shown in FIG. 1(a), the shortest distance between the spot dose and the mold body surface is greater than the effective electron range δ, and thus the calculation region is determined as a sphere with the spot dose as the center and δ as the radius.
In the embodiment shown in FIG. 1(b), where the shortest distance between the spot dose and the phantom surface is less than an effective electron range δ, the calculated area is determined as the portion of the spherical cap of the sphere below the phantom surface, such as the gray area in FIG. 2 (b).
it is further preferred in an embodiment that the three-dimensional voxels are uniform in size.
the present embodiment further preferably wherein the two-dimensional flux grid is a uniform two-dimensional grid.
The formula for calculating the voxel deposition dose at the point dose is as follows:
wherein the content of the first and second substances,
i is the two-dimensional flux grid designation,
n is the total number of flux grids,
j is the three-dimensional voxel label in the calculation region,
m is the total number of voxels within the computed region,
ωifor each flux grid weight in the monte carlo algorithm,
Dijfor the dose contributed by the ith flux grid to the jth voxel,
Djis the total dose deposited by the jth voxel.
Example 2
The present invention also provides a computing device comprising:
one or more processors;
a memory; and
one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, the one or more programs comprising instructions for a Monte Carlo-based spot dose calculation method, the method comprising the steps of:
selecting a dose calculation point in a three-dimensional image of a patient or phantom;
determining a calculation area;
carrying out three-dimensional gridding on the calculation area, wherein each grid is a voxel;
dividing any section in the incident direction of the beam into a two-dimensional flux grid; dose calculation is performed using a two-dimensional flux grid based on a monte carlo particle transport model.
Example 3
A computer readable storage medium storing one or more programs, the one or more programs comprising instructions adapted to be loaded from a memory and to perform the monte carlo-based spot dose calculation method described above, the method comprising the steps of:
selecting a dose calculation point in a three-dimensional image of a patient or phantom;
determining a calculation area;
carrying out three-dimensional gridding on the calculation area, wherein each grid is a voxel;
dividing any section in the incident direction of the beam into a two-dimensional flux grid; dose calculation is performed using a two-dimensional flux grid based on a monte carlo particle transport model.
In the embodiment 1 of the present invention, a monte carlo-based point dose calculation method is adopted, and the calculation region set by the present invention can fully consider the influence of particle scattering on the dose, thereby reducing unnecessary dose calculation tasks while ensuring accurate calculation, increasing the calculation speed, and facilitating clinical application.
Those skilled in the art will appreciate that the modules in the device in an embodiment may be adaptively changed and disposed in one or more devices different from the embodiment. The modules or units or components of the embodiments may be combined into one module or unit or component, and furthermore they may be divided into a plurality of sub-modules or sub-units or sub-components. All of the features disclosed in this specification (including any accompanying claims, abstract and drawings), and all of the processes or elements of any method or apparatus so disclosed, may be combined in any combination, except combinations where at least some of such features and/or processes or elements are mutually exclusive. Each feature disclosed in this specification (including any accompanying claims, abstract and drawings) may be replaced by alternative features serving the same, equivalent or similar purpose, unless expressly stated otherwise.
As used herein, unless otherwise specified the use of the ordinal adjectives "first", "second", "third", etc., to describe a common object, merely indicate that different instances of like objects are being referred to, and are not intended to imply that the objects so described must be in a given sequence, either temporally, spatially, in ranking, or in any other manner.
Furthermore, those skilled in the art will appreciate that while some embodiments described herein include some features included in other embodiments, rather than other features, combinations of features of different embodiments are meant to be within the scope of the invention and form different embodiments. For example, in the following claims, any of the claimed embodiments may be used in any combination.
It should be understood that the various techniques described herein may be implemented in connection with hardware or software or, alternatively, with a combination of both. Thus, the methods and apparatus of the present invention, or certain aspects or portions thereof, may take the form of program code (i.e., instructions) embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other machine-readable storage medium, wherein, when the program is loaded into and executed by a machine, such as a computer, the machine becomes an apparatus for practicing the invention.
By way of example, and not limitation, computer readable media may comprise computer storage media and communication media. Computer storage media store information such as computer readable instructions, data structures, program modules or other data. Communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media. Combinations of any of the above are also included within the scope of computer readable media.
The embodiments described above are intended to facilitate one of ordinary skill in the art in understanding and using the present invention. It will be readily apparent to those skilled in the art that various modifications to these embodiments may be made, and the generic principles described herein may be applied to other embodiments without the use of the inventive faculty. Therefore, the present invention is not limited to the embodiments described herein, and those skilled in the art should make improvements and modifications within the scope of the present invention based on the disclosure of the present invention.
Claims (9)
1. A monte carlo-based spot dose calculation method adapted to be executed in a computing device, characterized by: the method comprises the following steps:
selecting a dose calculation point in a three-dimensional image of a patient or phantom;
determining a calculation area;
carrying out three-dimensional gridding on the calculation area, wherein each grid is a voxel;
dividing any section in the incident direction of the beam into a two-dimensional flux grid; dose calculation is performed using a two-dimensional flux grid based on a monte carlo particle transport model.
2. The monte carlo-based spot dose calculation method according to claim 1, wherein: the dose calculation point is located in the region of interest.
3. The monte carlo-based spot dose calculation method according to claim 2, wherein: the method for determining the region of interest of the user comprises the following steps: obtained by user-defined or by calculation based on one or a combination of physical factors and biomedical factors.
4. The monte carlo-based spot dose calculation method according to claim 3, wherein: the physical factors reflect the material composition and the irradiation physical conditions of the patient or the phantom; wherein the material composition of the patient or the phantom comprises: density, CT value, mass number and atomic number of the die body; the irradiation physical conditions include: field distribution and source distribution;
the biomedical factors comprise: organ tissue irradiation threshold, biological sensitivity, secondary cancer occurrence probability distribution, tumor killing probability distribution, and injury probability.
5. The monte carlo-based spot dose calculation method according to claim 1, wherein: the method for determining the calculation area comprises the following steps:
(1) determining a sphere by taking the position of the point dose as a sphere center and taking the distance of an effective electron range as a radius;
(2) the portion of the sphere below the surface of the phantom is the calculated area.
6. The monte carlo-based spot dose calculation method according to claim 1, wherein: the two-dimensional flux grid is a uniform two-dimensional grid.
7. The monte carlo-based spot dose calculation method according to claim 1, wherein: the three-dimensional voxels are uniform in size.
8. A computing device, comprising:
one or more processors;
a memory; and
one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for the monte carlo-based spot dose calculation method of any of the preceding claims 1-7.
9. A computer readable storage medium storing one or more programs, the one or more programs comprising instructions adapted to be loaded from a memory and to perform the monte carlo based spot dose calculation method of any of the preceding claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810710921.5A CN110675932A (en) | 2018-07-03 | 2018-07-03 | Monte Carlo-based point dose calculation method, equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810710921.5A CN110675932A (en) | 2018-07-03 | 2018-07-03 | Monte Carlo-based point dose calculation method, equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110675932A true CN110675932A (en) | 2020-01-10 |
Family
ID=69065623
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810710921.5A Pending CN110675932A (en) | 2018-07-03 | 2018-07-03 | Monte Carlo-based point dose calculation method, equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110675932A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113539420A (en) * | 2021-06-30 | 2021-10-22 | 上海联影医疗科技股份有限公司 | Dose control system, device and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6029079A (en) * | 1997-05-22 | 2000-02-22 | Regents Of The University Of California | Evaluated teletherapy source library |
| US6285969B1 (en) * | 1998-05-22 | 2001-09-04 | The Regents Of The University Of California | Use of single scatter electron monte carlo transport for medical radiation sciences |
| US20020106054A1 (en) * | 2000-09-22 | 2002-08-08 | Numerix, Llc. | Radiation therapy treatment method |
| CN101477203A (en) * | 2009-01-22 | 2009-07-08 | 中国科学技术大学 | Resolution Monte Carto dosage computing method |
| CN104407374A (en) * | 2014-12-10 | 2015-03-11 | 上海联影医疗科技有限公司 | Dose distribution determination method and device |
-
2018
- 2018-07-03 CN CN201810710921.5A patent/CN110675932A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6029079A (en) * | 1997-05-22 | 2000-02-22 | Regents Of The University Of California | Evaluated teletherapy source library |
| US6285969B1 (en) * | 1998-05-22 | 2001-09-04 | The Regents Of The University Of California | Use of single scatter electron monte carlo transport for medical radiation sciences |
| US20020106054A1 (en) * | 2000-09-22 | 2002-08-08 | Numerix, Llc. | Radiation therapy treatment method |
| CN101477203A (en) * | 2009-01-22 | 2009-07-08 | 中国科学技术大学 | Resolution Monte Carto dosage computing method |
| CN104407374A (en) * | 2014-12-10 | 2015-03-11 | 上海联影医疗科技有限公司 | Dose distribution determination method and device |
Non-Patent Citations (1)
| Title |
|---|
| 范佳锦等: "蒙特卡罗法内照射剂量计算中不同物理模型的比较", 《核电子学与探测技术》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113539420A (en) * | 2021-06-30 | 2021-10-22 | 上海联影医疗科技股份有限公司 | Dose control system, device and storage medium |
| CN113539420B (en) * | 2021-06-30 | 2023-07-25 | 上海联影医疗科技股份有限公司 | Dose control system, device and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2760543B1 (en) | Dose computation for radiation therapy using heterogeneity compensated superposition | |
| CN102036712B (en) | System for radiation therapy and method for determining radiation therapy parameters | |
| US10417390B2 (en) | Methods and systems for radiotherapy treatment planning | |
| US20060259282A1 (en) | Deterministic computation of radiation transport for radiotherapy dose calculations and scatter correction for image reconstruction | |
| CN104548372A (en) | Radiotherapy planning method and device, radiotherapy dose determining method and device and radiotherapy quality guaranteeing method and device | |
| CN110504016B (en) | Monte Carlo grid parallel dose calculation method, equipment and storage medium | |
| EP3742702B1 (en) | Standardized cloud radiotherapy planning method and storage medium | |
| CN110556176B (en) | Monte Carlo-based dose optimization method, equipment and storage medium | |
| CN110310720B (en) | Monte Carlo dose calculation method, equipment and storage medium | |
| CN107050667B (en) | Proton and heavy ion dose calculation methodology under magnetic field based on GPU Monte carlo algorithm | |
| CN107073286A (en) | Radiation therapy planning device, radiation therapy planning method and radiation treatment systems | |
| Penfold | Image reconstruction and Monte Carlo simulations in the development of proton computed tomography for applications in proton radiation therapy | |
| US20140032185A1 (en) | Gpu-based fast dose calculator for cancer therapy | |
| US9251302B2 (en) | System and method for planning a radiation therapy treatment | |
| CN110570923A (en) | mixed Monte Carlo radiotherapy reverse optimization method, equipment and storage medium | |
| Jones et al. | Characterization of Compton-scatter imaging with an analytical simulation method | |
| CN113677396A (en) | Method and system for evaluating radiation therapy treatment plans | |
| Sikora | Virtual source modelling of photon beams for Monte Carlo based radiation therapy treatment planning | |
| CN110675932A (en) | Monte Carlo-based point dose calculation method, equipment and storage medium | |
| CN110310743B (en) | Monte Carlo dose calculation method, equipment and storage medium | |
| CN110554423B (en) | Radiation dose calculation system | |
| Cranmer-Sargison et al. | EUD-based radiotherapy treatment plan evaluation: incorporating physical and Monte Carlo statistical dose uncertainties | |
| CN116350959A (en) | Radiation therapy dose analysis method and system | |
| Scholz | Development and evaluation of advanced dose calculations for modern radiation therapy techniques | |
| CN114984464A (en) | Dose monitoring method, system and storage medium |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200110 |