CN110310743A - A kind of Monte Carto dosage computing method, equipment and storage medium - Google Patents
A kind of Monte Carto dosage computing method, equipment and storage medium Download PDFInfo
- Publication number
- CN110310743A CN110310743A CN201810227911.6A CN201810227911A CN110310743A CN 110310743 A CN110310743 A CN 110310743A CN 201810227911 A CN201810227911 A CN 201810227911A CN 110310743 A CN110310743 A CN 110310743A
- Authority
- CN
- China
- Prior art keywords
- grid
- weight
- particle
- monte
- 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.)
- Granted
Links
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
-
- 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
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
Abstract
The invention belongs to dose of radiations to calculate field, be related to a kind of Monte Carto dosage computing method, equipment and storage medium.The method comprising the steps of: (1) pre-processing, including input patient's image, delineate information, field size, direction of illumination and source parameter;(2) model treatment: rebuilding patient's image, is threedimensional model by two-dimensional patient's image reconstruction, and threedimensional model is carried out uniform grid;(3) determine weight grid: i. determines three dimensional weight grids;Ii. two dimension modulus grid is determined;(4) particle input simulation: when incoming particle passes through two dimension modulus grid, particle carries out weight changes according to importance, calls Monte Carlo database, carries out PARTICLE TRANSPORT FROM simulation using Monte Carlo PARTICLE TRANSPORT FROM principle;(5) result exports.The present invention creatively proposes weight grid method, can greatly increase the efficiency of PARTICLE TRANSPORT FROM;Precision is improved to the area-of-interest of user's concern.
Description
Technical field
The invention belongs to dose of radiations to calculate field, be related to a kind of Monte Carto dosage computing method, equipment and storage and be situated between
Matter.
Background technique
Rapid Dose Calculation is one of the core content of radiotherapy in the treatment planning system, is rapidly and accurately provided suffered in region of interest
The data of exposure dose, it is most important to the formulation of radiotherapy planning, how under the premise of guaranteeing Rapid Dose Calculation precision to reduce agent
Meter evaluation time is the main bottleneck for formulating radiotherapy treatment planning.
There are mainly two types of the methods for improving Rapid Dose Calculation speed, first is that using different dose calculation methodologies, second is that by
In the stronger hardware of computing capability.The Dose calculation algorithm studied and used in radiotherapy at present is divided into three kinds substantially, according to
Computational accuracy calculates the time from being more to successively Monte Carlo EGS4 method less from high to low, and differential convolution superposition algorithm and pencil beam are calculated
Method.Monte Carlo EGS4 method simulates the overall process of particle and matter interaction usually as the standard of Rapid Dose Calculation, can calculate
The dosage of various complex conditions is distributed, therefore Monte Carlo EGS4 method is referred to as the goldstandard of industry, is all Rapid Dose Calculation essences
Spend highest method.But the error that provides of existing Monte Carlo Calculation be it is non-uniform, cause error none very well
Standard, the big some places error of some places error is small, can not judging result quality.
If can provide a user Region Of Interest application condition it is uniform as a result, if greatly facilitate assessment errors, and
Computational accuracy can be controlled very well, and reach accurate calculation as a result, so that substantially reduce risk.
Summary of the invention
It is an object of the invention to provide a kind of Monte Carto dosage calculation side to overcome the defect of the above-mentioned prior art
Method, equipment and storage medium.
To achieve the above object, the invention adopts the following technical scheme:
A kind of Monte Carto dosage computing method, suitable for being executed in calculating equipment, comprising steps of
(1) it pre-processes, including inputs patient's image, delineates information, field size, direction of illumination and source parameter;
(2) model treatment: rebuilding patient's image, is threedimensional model by two-dimensional patient's image reconstruction, and by three
Dimension module carries out uniform grid;
(3) weight grid is determined:
I. three dimensional weight grids are determined;
Ii. two dimension modulus grid is determined;
(4) particle input simulation: when incoming particle passes through two dimension modulus grid, particle is weighed according to importance
Change again, call Monte Carlo database, carries out PARTICLE TRANSPORT FROM simulation using Monte Carlo PARTICLE TRANSPORT FROM principle;
(5) result exports: output includes that dosage distribution is distributed with standard deviation.
The parameter in the source includes the energy in source, position, direction, particle types, weight;The weights initialisation is
1;
The determination method of the three dimensional weight grids is customized by the user or by using based on physical agent, biology
One of medicine factor or combination, which calculate, to be obtained.
The physical agent is to reflect that the material of patient or die body constitutes, irradiates physical condition;The material structure
At including but not limited to: density, mass number, the atomicity of die body;The irradiation physical condition includes: launched field distribution, source point
Cloth.
The biomedical factor includes but is not limited to: organ-tissue exposure threshold, biological susceptibility, secondary cancer hair
Raw probability distribution, tumour kill probability distribution, damage probability.
The determination method of the two dimension modulus grid is that three dimensional weight mesh mappings are become two to source plane of incidence
Dimensional weight grid;Include the following steps:
(a) shape of two dimension modulus grid is set, can be arbitrary polygon, it is therefore preferable to regular polygon, Huo Zheneng
Enough merge the grid of rectangular polygon;
(b) weight of three-dimensional grid be mapped on two dimension modulus grid according to following formula:
Wherein, WijFor the weight of two dimension modulus grid (i, j);
For two dimension modulus grid (i, j) central point to the distance of three dimensional weight grid (i ', j ', k ') central points
wi′j′k′For the weight distribution of three dimensional weight grids (i ', j ', k ').
The effect of two dimension modulus grid is source particle after two dimension modulus grid, is changed to the weight of source particle.
Preferably, a weight threshold or threshold interval can be set, the population for being more than or equal to the weight threshold (section) is carried out
Division, the weighted value of population before weight adduction=division of each particle after division;To the particle for being less than the threshold value (section)
Russian roulette is carried out, continues to inject after gambling the threshold weights of the particle acquisition setting of win, the particle being out of the money will be killed.
In step (4), the Monte Carlo database includes cross-section library, material depot;The particle includes but unlimited
In photon, electronics, proton, heavy ion, neutron.
The present invention also provides a kind of calculating equipment, comprising:
One or more processors;
Memory;And
One or more programs, wherein one or more of programs are stored in the memory and are configured as by one
A or multiple processors execute, and one or more programs include the finger for above-mentioned Monte Carto dosage computing method
It enables.
The present invention also provides a kind of computer readable storage medium for storing one or more programs, described one or more
A program includes instruction, and described instruction is suitable for being loaded by memory and executing above-mentioned Monte Carto dosage computing method.
The invention has the following advantages:
1, the present invention creatively proposes weight grid method, by three-dimensional interested area maps to two-dimensional grid
In, the efficiency of PARTICLE TRANSPORT FROM can be considerably increased;
2, it is also avoided that the calculating of inactive area simultaneously, improves the computational accuracy of the area-of-interest of user's concern.
Detailed description of the invention
Fig. 1 is the flow chart of Monte Carto dosage computing method in a preferred embodiment of the invention.
Specific embodiment
The present invention is further illustrated with attached drawing with reference to embodiments.
Embodiment 1
A kind of Monte Carto dosage computing method, suitable for being executed in calculating equipment, comprising steps of
(1) 110 are pre-processed, including inputs patient's image, delineate information, field size, direction of illumination and source parameter;Wherein
The parameter in source includes the energy in source, position, direction, particle types, weight;It is 1 by weights initialisation;
(2) model treatment 120: rebuilding patient's image, is threedimensional model by two-dimensional patient's image reconstruction, and will
Threedimensional model carries out uniform grid;For sizing grid according to the customized determination of user, the smaller computational accuracy of grid is higher, calculates and appoints
It is engaged in more complicated;
(3) weight grid 130 is determined:
I. three dimensional weight grids 131 are determined;The wherein determination method of three dimensional weight grids, is customized by the user or by adopting
It is calculated and is obtained with based on one of physical agent, the biomedical factor or combination.Physical agent is reflection patient or mould
The material of body constitutes, irradiates physical condition;Material composition includes but is not limited to: density, mass number, the atomicity of die body;Irradiation
Physical condition includes: launched field distribution, source distribution.The biomedical factor includes but is not limited to: organ-tissue exposure threshold, and biology is quick
Perception, secondary cancer probability of happening distribution, tumour kill probability distribution, damage probability;
In the present embodiment, three-dimensional is calculated using die body density and the weight of biological susceptibility (grade organ of endangering in launched field)
Weight grid determines, the steps include:
A) density of die body is normalized, up to 1, minimum 0;Weight, danger in launched field are set to different zones
And organ is set as 1, tumour is set as 0.5, other are set as 0.3, and jeopardizing organ setting weight outside launched field is 0.5, other areas
Domain is set as 0;
B) weight summation and normalization: the weight for normalizing density and region weight are subjected to arithmetic, and returned again
One changes;Obtained normalized weight is three dimensional weight grids;The arithmetic formula is determined by user preset.This implementation
The weight parameter of unspecified biotic factor or physical agent is defaulted as 1 in example.
Ii. determine two dimension modulus grid 132: the determination method of two dimension modulus grid is to arrive three dimensional weight mesh mappings
Source plane of incidence becomes two dimension modulus grid;Include the following steps:
(a) shape of two dimension modulus grid is arranged in, can be arbitrary polygon, it is therefore preferable to regular polygon, Huo Zheneng
Enough merge the grid of rectangular polygon;
(b) weight of three-dimensional grid is mapped on two dimension modulus grid by according to formula (1):
Wherein, WijFor the weight of two dimension modulus grid (i, j);
For two dimension modulus grid (i, j) central point to the distance of three dimensional weight grid (i ', j ', k ') central points
wi′j′k′For the weight distribution of three dimensional weight grids (i ', j ', k ').
The effect of two dimension modulus grid is source particle after two dimension modulus grid, is changed to the weight of source particle.
In the present embodiment preferably, a weight threshold or threshold interval can be set, to more than or equal to the weight threshold (section)
Population is divided, the weighted value of population before weight adduction=division of each particle after division;To less than the threshold value (area
Between) particle carry out Russian roulette, gamble win particle obtain the threshold weights of setting after continue to inject, the particle that is out of the money will
It is killed.
(4) particle input simulation 140: when incoming particle pass through two dimension modulus grid when, particle according to importance into
Row weight changes;Monte Carlo database is called, carries out PARTICLE TRANSPORT FROM simulation using Monte Carlo PARTICLE TRANSPORT FROM principle;Obtain agent
Amount distribution (according to formula (2)) and standard deviation distribution (according to formula (3));Wherein, Monte Carlo database includes cross-section library, material
Expect library;Particle includes photon, electronics, proton, heavy ion, neutron;
D(x,y,z)=Edep/(ρ(x,y,z)V(x,y,z)) (2)
Wherein, D(x,y,z)For the dose value of grid locating for grid centre coordinate (x, y, z) in the body mould or human body of gridding;
EdepFor the sedimentary energy of grid locating for certain grid element center coordinate (x, y, z);
V(x,y,z)For the volume of grid locating for certain grid element center coordinate (x, y, z);
ρ(x,y,z)For the averag density of grid locating for certain grid element center coordinate (x, y, z);
Wherein,For the standard deviation for reaching the mean dose that all particles of volume element generate;
xiTo reach the dose value that i-th of particle of volume element generates;
For the mean dose of all particles of arrival volume element;
N is the particle number for reaching volume element;
(5) result output 150: output includes that dosage distribution is distributed with standard deviation.
Embodiment 2
A kind of calculating equipment, comprising:
One or more processors;
Memory;And
One or more programs, wherein one or more programs are stored in memory and are configured as by one or more
Processor executes, wherein one or more programs include the instruction for above-mentioned Monte Carto dosage computing method, and above-mentioned illiteracy is special
Caro dose calculation methodology comprising steps of
(1) it pre-processes: including input patient's image, delineating information, field size, direction of illumination and source parameter;
(2) model treatment: rebuilding patient's image, is threedimensional model by two-dimensional patient's image reconstruction, and by three
Dimension module carries out uniform grid;
(3) weight grid is determined:
I. three dimensional weight grids are determined;
Ii. two dimension modulus grid is determined;
(4) particle input simulation: when incoming particle passes through two dimension modulus grid, particle is weighed according to importance
Change again, call Monte Carlo database, carries out PARTICLE TRANSPORT FROM simulation using Monte Carlo PARTICLE TRANSPORT FROM principle;
(5) result exports: output includes that dosage distribution is distributed with standard deviation.
Embodiment 3
A kind of computer readable storage medium storing one or more programs, wherein one or more programs include referring to
It enables, which is suitable for being loaded by memory and being executed above-mentioned Monte Carto dosage computing method, and the method comprising the steps of:
(1) it pre-processes: including input patient's image, delineating information, field size, direction of illumination and source parameter;
(2) model treatment: rebuilding patient's image, is threedimensional model by two-dimensional patient's image reconstruction, and by three
Dimension module carries out uniform grid;
(3) weight grid is determined:
I. three dimensional weight grids are determined;
Ii. two dimension modulus grid is determined;
(4) particle input simulation: when incoming particle passes through two dimension modulus grid, particle is weighed according to importance
Change again, call Monte Carlo database, carries out PARTICLE TRANSPORT FROM simulation using Monte Carlo PARTICLE TRANSPORT FROM principle;
(5) result exports: output includes that dosage distribution is distributed with standard deviation.
The Monte Carto dosage that the above embodiment of the present invention provides calculates precision methods and creatively proposes weight grid
Method can considerably increase the efficiency of PARTICLE TRANSPORT FROM;The computational accuracy of the area-of-interest of user's concern can also be improved.
It should be appreciated that various technologies described herein are realized together in combination with hardware or software or their combination.From
And some aspects or part of the process and apparatus of the present invention or the process and apparatus of the present invention can take the tangible matchmaker of insertion
It is situated between, such as the program code in floppy disk, CD-ROM, hard disk drive or other any machine readable storage mediums (refers to
Enable) form, wherein when program is loaded into the machine of such as computer etc, and when being executed by the machine, which becomes real
Trample equipment of the invention.
By way of example and not limitation, computer-readable medium includes computer storage media and communication media.It calculates
Machine storage medium stores the information such as computer readable instructions, data structure, program module or other data.Communication media one
As with the modulated message signals such as carrier wave or other transmission mechanisms embody computer readable instructions, data structure, program
Module or other data, and including any information transmitting medium.Above any combination is also included within computer-readable
Within the scope of medium.
This hair can be understood and applied the above description of the embodiments is intended to facilitate those skilled in the art
It is bright.Person skilled in the art obviously easily can make various modifications to these embodiments, and described herein
General Principle is applied in other embodiments without having to go through creative labor.Therefore, the present invention is not limited to implementations here
Example, those skilled in the art's announcement according to the present invention, improvement and modification made without departing from the scope of the present invention all should be
Within protection scope of the present invention.
Claims (10)
1. a kind of Monte Carto dosage computing method, suitable for being executed in calculating equipment, it is characterised in that: comprising steps of
(1) it pre-processes: including input patient's image, delineating information, field size, direction of illumination and source parameter;
(2) model treatment: rebuilding patient's image, is threedimensional model by two-dimensional patient's image reconstruction, and by three-dimensional mould
Type carries out uniform grid;
(3) weight grid is determined:
I. three dimensional weight grids are determined;
Ii. two dimension modulus grid is determined;
(4) particle input simulation: when incoming particle passes through two dimension modulus grid, particle carries out weight according to importance and changes
Become, call Monte Carlo database, carries out PARTICLE TRANSPORT FROM simulation using Monte Carlo PARTICLE TRANSPORT FROM principle;
(5) result exports: output includes that dosage distribution is distributed with standard deviation.
2. Monte Carto dosage computing method according to claim 1, it is characterised in that: the source parameter includes source
One of energy, position, direction, particle types and weight or more than one.
3. Monte Carto dosage computing method according to claim 2, it is characterised in that: the weights initialisation is 1.
4. Monte Carto dosage computing method according to claim 1, it is characterised in that: the three dimensional weight grids
It determines method, be customized by the user or calculated by using based on one of physical agent, the biomedical factor or combination
It obtains.
5. Monte Carto dosage computing method according to claim 4, it is characterised in that: the physical agent is reflection
The material of patient or die body constitutes, irradiates physical condition;
Wherein material constitutes density, mass number or atomicity including die body;Irradiation physical condition includes: launched field distribution or source point
Cloth;
The biomedical factor includes organ-tissue exposure threshold, biological susceptibility, the distribution of secondary cancer probability of happening, swells
Tumor kills probability distribution or damage probability.
6. Monte Carto dosage computing method according to claim 1, it is characterised in that: the two dimension modulus grid
It determines method, is that three dimensional weight mesh mappings to source plane of incidence are become into two dimension modulus grid;Include the following steps:
(a) shape of two dimension modulus grid is set, and the shape of the grid is arbitrary polygon;
(b) weight of three-dimensional grid be mapped on two dimension modulus grid according to formula (1):
Wherein, WijFor the weight of two dimension modulus grid (i, j);
For two dimension modulus grid (i, j) central point to the distance of three dimensional weight grid (i ', j ', k ') central points
wi′j′k′For the weight distribution of three dimensional weight grids (i ', j ', k ').
7. Monte Carto dosage computing method according to claim 6, it is characterised in that: the two dimension modulus grid
Shape is regular polygon or can merge rectangular polygon.
8. Monte Carto dosage computing method according to claim 1, it is characterised in that: in step (4), the illiteracy is special
Caro database includes cross-section library, material depot;The particle includes photon, electronics, proton, heavy ion or neutron.
9. a kind of calculating equipment, comprising:
One or more processors;
Memory;And
One or more programs, wherein the storage of one or more of programs in the memory and be configured as by one or
Multiple processors execute, and one or more programs include for the Meng Teka any in the claims 1-8
The instruction of sieve dose calculation methodology.
10. a kind of computer readable storage medium for storing one or more programs, one or more programs include referring to
It enables, described instruction is suitable for being loaded by memory and being executed any Monte Carto dosage in the claims 1-8 and calculates
Method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810227911.6A CN110310743B (en) | 2018-03-20 | 2018-03-20 | Monte Carlo dose calculation method, equipment and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810227911.6A CN110310743B (en) | 2018-03-20 | 2018-03-20 | Monte Carlo dose calculation method, equipment and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110310743A true CN110310743A (en) | 2019-10-08 |
CN110310743B CN110310743B (en) | 2022-05-06 |
Family
ID=68073826
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810227911.6A Active CN110310743B (en) | 2018-03-20 | 2018-03-20 | Monte Carlo dose calculation method, equipment and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110310743B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110782952A (en) * | 2019-10-25 | 2020-02-11 | 深圳大学 | Quantitative Monte Carlo simulation method and system for absorbed dose and storage medium |
WO2022143333A1 (en) * | 2020-12-31 | 2022-07-07 | 中硼(厦门)医疗器械有限公司 | Radiation irradiation system and control method therefor |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020027971A1 (en) * | 2000-05-05 | 2002-03-07 | Deasy Joseph O. | Method and apparatus for radiotherapy treatment planning |
US20130030762A1 (en) * | 2009-07-15 | 2013-01-31 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Method for calculating doses deposited by ionizing radiation |
CN104407374A (en) * | 2014-12-10 | 2015-03-11 | 上海联影医疗科技有限公司 | Dose distribution determination method and device |
CN104750992A (en) * | 2015-04-01 | 2015-07-01 | 上海联影医疗科技有限公司 | Method, device and system for simulating particle transport and calculating body dosage in radiotherapy |
CN105930637A (en) * | 2016-04-02 | 2016-09-07 | 中北大学 | Radiotherapy scheme optimization method capable of automatically determining weight of target function |
-
2018
- 2018-03-20 CN CN201810227911.6A patent/CN110310743B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020027971A1 (en) * | 2000-05-05 | 2002-03-07 | Deasy Joseph O. | Method and apparatus for radiotherapy treatment planning |
US20130030762A1 (en) * | 2009-07-15 | 2013-01-31 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Method for calculating doses deposited by ionizing radiation |
CN104407374A (en) * | 2014-12-10 | 2015-03-11 | 上海联影医疗科技有限公司 | Dose distribution determination method and device |
CN104750992A (en) * | 2015-04-01 | 2015-07-01 | 上海联影医疗科技有限公司 | Method, device and system for simulating particle transport and calculating body dosage in radiotherapy |
CN105930637A (en) * | 2016-04-02 | 2016-09-07 | 中北大学 | Radiotherapy scheme optimization method capable of automatically determining weight of target function |
Non-Patent Citations (4)
Title |
---|
HELIO YORIYAZ 等: "Monte Carlo MCNP-4B-Based Absorbed Dose Distribution Estimates for Patient-Specific Dosimetry", 《THE JOURNAL OF NUCLEAR MEDICINE》 * |
ZAKARIAN, C 等: "Beamlet dose distribution compression and reconstruction using wavelets for intensity modulated treatment planning", 《MEDICAL PHYSICS》 * |
吴宜灿 等: "超级蒙特卡罗核计算仿真软件系统SuperMC", 《核科学与工程》 * |
李贵 等: "基于多算法的电子剂量场反演方法", 《核技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110782952A (en) * | 2019-10-25 | 2020-02-11 | 深圳大学 | Quantitative Monte Carlo simulation method and system for absorbed dose and storage medium |
CN110782952B (en) * | 2019-10-25 | 2023-06-27 | 深圳大学 | Quantitative Monte Carlo simulation method, system and storage medium for absorbed dose |
WO2022143333A1 (en) * | 2020-12-31 | 2022-07-07 | 中硼(厦门)医疗器械有限公司 | Radiation irradiation system and control method therefor |
Also Published As
Publication number | Publication date |
---|---|
CN110310743B (en) | 2022-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3695882A1 (en) | Computing radiotherapy dose distribution | |
Yepes et al. | A GPU implementation of a track-repeating algorithm for proton radiotherapy dose calculations | |
Jahnke et al. | GMC: a GPU implementation of a Monte Carlo dose calculation based on Geant4 | |
CN110310720A (en) | A kind of Monte Carto dosage computing method, equipment and storage medium | |
CN104407374B (en) | The determination method and device of dosage distribution | |
US20070071169A1 (en) | Method for verification of intensity modulated radiation therapy | |
CN107072624A (en) | system and method for automatic treatment plan | |
CN108766563A (en) | Radiotherapy prediction of result method and system based on dosage group | |
CN108415058A (en) | The dose calculation methodology and system of radioactive ray | |
CN110302475B (en) | Cloud Monte Carlo dose verification analysis method, equipment and storage medium | |
CN107551411B (en) | Proton heavy ion intensity modulated radiotherapy robust optimization method aiming at range uncertainty | |
CN110504016A (en) | A kind of Monte Carlo grid parallel computation dose calculation methodology, equipment and storage medium | |
US9251302B2 (en) | System and method for planning a radiation therapy treatment | |
US20140032185A1 (en) | Gpu-based fast dose calculator for cancer therapy | |
Penfold | Image reconstruction and Monte Carlo simulations in the development of proton computed tomography for applications in proton radiation therapy | |
Jones et al. | Characterization of Compton-scatter imaging with an analytical simulation method | |
Pastor-Serrano et al. | Millisecond speed deep learning based proton dose calculation with Monte Carlo accuracy | |
Da Silva et al. | Fast pencil beam dose calculation for proton therapy using a double-Gaussian beam model | |
Toscano et al. | Impact of machine log-files uncertainties on the quality assurance of proton pencil beam scanning treatment delivery | |
CN110310743A (en) | A kind of Monte Carto dosage computing method, equipment and storage medium | |
Wahl et al. | Efficiency of analytical and sampling-based uncertainty propagation in intensity-modulated proton therapy | |
Maneval et al. | pGPUMCD: an efficient GPU-based Monte Carlo code for accurate proton dose calculations | |
CN109407134B (en) | Dose distribution calculation method and system | |
Zhong et al. | Generation of a novel phase‐space‐based cylindrical dose kernel for IMRT optimization | |
CN110826207A (en) | OCT imaging rapid Monte Carlo model of dental tissue in any shape |
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 |