CN104298871A - Method for accelerating voxel human body model dose evaluation based on GPU acceleration in radiation protection - Google Patents
Method for accelerating voxel human body model dose evaluation based on GPU acceleration in radiation protection Download PDFInfo
- Publication number
- CN104298871A CN104298871A CN201410528520.XA CN201410528520A CN104298871A CN 104298871 A CN104298871 A CN 104298871A CN 201410528520 A CN201410528520 A CN 201410528520A CN 104298871 A CN104298871 A CN 104298871A
- Authority
- CN
- China
- Prior art keywords
- voxel
- gpu
- dosage
- dose
- manikin
- 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
Landscapes
- Radiation-Therapy Devices (AREA)
Abstract
The invention discloses a method for accelerating voxel human body model dose evaluation based on GPU acceleration in radiation protection. The method is different from a conventional method in that a voxel equivalent dose in dose calculation is processed by use of a Graphic Processor Unit (GPU) in combination with serial calculation of CPU complex logics, events and the like; the process has the characteristics of large data size, data dependency, same executive program for data, high parallelism and calculation density and the like. As a result, with the GPU, more calculation units can be put in parallel execution, and therefore, the calculation time is shortened, the calculation efficiency is improved and the real-time calculation of the voxel human body model dose is realized.
Description
Technical field
The present invention relates in a kind of radiation protection based on the voxel manikin dose assessment accelerated method that GPU accelerates, be in the calculation of voxel manikin radiation protection quantimeter, determine the method for the 3-dimensional dose distribution of radiation in irradiated voxel model quickly and accurately.
Background technology
Manikin is radiation proof important tool.Manikin is roughly divided into three by its development: Mass Model, mathematical model and voxel model.Voxel model is made up of millions of corpusculum units, and voxel model more can reflect human anatomic structure, more truly, therefore, more accurate based on voxel model Rapid Dose Calculation.
In dose assessment, determine that the 3-dimensional dose distribution of radiating particle in body mould is a basic problem quickly and accurately.The method of the determination dosage of current use is, based on the radiation field that Monte Carlo method or analytic method calculate, the method utilizing flux dosage to change calculates the dosage of each little volume elements.
GPU (Graphic Processing Unit) applies the acceleration shown with figure at first.The processing mode of the single-instruction multiple-data stream (SIMD) (SIMD:Single Instruction Multiple Data) of GPU can operate large-scale data concurrently, greatly can shorten computing time.
GPU is applied in the Rapid Dose Calculation of voxel manikin by the present invention.By carrying out the Rapid Dose Calculation of a large amount of voxel simultaneously, the computational accuracy of its computational accuracy and use CPU (Central Process Unit) is suitable, comparatively uses its computing time CPU platform to calculate and greatly reduces.
Realize the parallel computation based on GPU, that first understands support GPU computing becomes instrument.The software and hardware architecture that of CUDA (Computer Unified Device Architecture) to be NVIDIA company be GPU becomes proposition is brand-new, it is based on C language, writes out the program performed on GPU.Can referring to CUDA Programming Guide and CUDA reference manual about CUDA detailed description.During the present invention formally utilizes GPU to programme, CUDA technology is accelerated voxel manikin dosage.
Summary of the invention
Technology of the present invention is dealt with problems: overcome the deficiencies in the prior art, provides the voxel manikin dose assessment accelerated method based on GPU acceleration in a kind of radiation protection, can determine the 3-dimensional dose distribution of radiating particle in voxel model quickly and accurately.
Technical scheme of the present invention is: based on the voxel manikin dose assessment accelerated method that GPU accelerates in a kind of radiation protection, it is characterized in that: employ CPU and GPU and combine and carry out Rapid Dose Calculation, the processing procedure having parallel condition for some uses GPU to carry out acceleration process, and the serial computing such as complex logic and issued transaction still adopts CPU.Performing step comprises:
(1) on CPU platform, voxel manikin parameter is formed;
(2) GPU for determining voxel dosage is selected, for GPU memory allocated space, to the calculating data of the voxel manikin parameter in GPU transfer step (1), the data transmitted are needed to comprise voxel sequence, organ name and ID and U corresponding relation and voxel precision;
(3) dosage of each voxel of GPU parallel computation voxel manikin is utilized, dose assessment is on voxel model parameter basis, utilize flux dosage conversion method, the i.e. flux of known free space, by conversion factor, calculate the absorbed dose of each voxel respectively, utilize the dosage of each voxel of GPU parallel computation, the dosage of each GPU thread computes voxel, is finally on average Absorbed Doses by all voxel dosage.
The present invention's advantage is compared with prior art: the present invention is in conjunction with serial computing such as CPU complex logic and things, use graphic process unit (Graphic Processor Unit, GPU) voxel equivalent dose in treatment dosage calculating, this process exists that data volume is large, data dependence is low, data have identical executive routine, degree of parallelism and bulk density high.Thus by GPU, more computing unit executed in parallel can be allowed, shorten computing time, improve counting yield, realize voxel manikin dosage and calculate effect in real time.
Accompanying drawing explanation
Fig. 1 is realization flow figure of the present invention.
Embodiment
As shown in Figure 1, specific implementation step of the present invention is as follows:
(1) on CPU platform, form voxel manikin parameter, voxel manikin parameter comprises voxel sequence, organ name and ID and U corresponding relation and voxel precision;
(2) GPU for determining voxel dosage is selected, for GPU memory allocated space, to the calculating data of the voxel manikin parameter in GPU transfer step (1), the data transmitted are needed to comprise voxel sequence, organ name and ID and U corresponding relation and voxel precision;
(3) dosage of each voxel of GPU parallel computation voxel manikin is utilized.For the voxel of organ discretize, voxel is allowed to identify and thread identification one_to_one corresponding, the flux of voxel Dose location and the product of flux dose conversion factor.By voxel mark and thread identification one_to_one corresponding in dosage cumulative process, search the dosage that organ voxel is corresponding, the dosage of each for organ voxel is carried out cumulative average again.
The Rapid Dose Calculation accuracy of voxel model is, the human organ of voxel manikin is voxelization, is made up of a large amount of little rectangular parallelepiped voxels, and each voxel comprises organ code name and material composition.Dose assessment is on this model basis, calculates the absorbed dose of each voxel respectively, utilizes the dosage of each voxel of GPU parallel computation, finally all voxel dosage is on average Absorbed Doses.Can be expressed from the next.
Wherein,
for tissue T absorbed dose, D
t,jfor the absorbed dose of this organ voxel j, N is the voxel number that this organ comprises.
The acquisition of the voxel flux of voxel manikin is the centre coordinate by calculating voxel, then from radiation field grid, find out the grid of arest neighbors, and the flux of this grid is expressed as the flux of this voxel.
The absorbed dose of organ is the average of the flux forming each voxel location, with the result of the conversion factor product of this organ under this energy.Therefore, Absorbed Doses calculating is just converted on the flux of each voxel of this organ of acquisition.
The voxel of voxel model is that voxel fill order is first X-axis, then Y-axis, last Z axis according to certain aligned transfer.According to the fill order of voxel, all voxel number consecutivelies are traveled through, corresponding No. U, the organ of each voxel, represent that it is the voxel of this organ of composition.
In order to utilize the Synchronous data dispose function of GPU, accelerating computing velocity, distributing a thread to each voxel of voxel model and calculate.When calculating voxel dosage, first obtain voxel location, according to position enquiring flux, utilizing flux dose conversion factor, obtaining voxel dosage.
D
T,j＝φ
j(x,y,z)*CF (2)
Wherein, φ
j(x, y, z) for the flux at voxel j place (x, y, z), CF be flux dose conversion factor.
By voxel mark and thread identification one_to_one corresponding in dosage cumulative process, to the voxel at same organ, Parallel transformation goes out voxel dosage, then adds up to all voxel dosage of organ.
Wherein, D
t, andfor tissue T dosage and value.
Claims (2)
1. in a radiation protection based on GPU accelerate voxel manikin dose assessment accelerated method, it is characterized in that: employ CPU and GPU and combine and carry out Rapid Dose Calculation, the processing procedure having parallel condition for some uses GPU to carry out acceleration process, and the serial computing such as complex logic and issued transaction still adopts CPU; Performing step comprises:
(1) on CPU platform, voxel manikin parameter is formed;
(2) GPU for determining voxel dosage is selected, for GPU memory allocated space, to the calculating data of the voxel manikin parameter in GPU transfer step (1), the data transmitted are needed to comprise voxel sequence, organ name and ID and U corresponding relation and voxel precision;
(3) dosage of each voxel of GPU parallel computation voxel manikin is utilized, dose assessment is on voxel model parameter basis, utilize flux dosage conversion method, the i.e. flux of known free space, by conversion factor, calculate the absorbed dose of each voxel respectively, utilize the dosage of each voxel of GPU parallel computation, the dosage of each GPU thread computes voxel, is finally on average Absorbed Doses by all voxel dosage.
2. in a kind of radiation protection according to claim 1 based on GPU accelerate voxel manikin dose assessment accelerated method, it is characterized in that: dosage comprises photon doses, also neutron dose is comprised, photon doses is for evaluating the harm of photon irradiation to human body, and neutron dose is for evaluating the harm of neutron irradiation to human body.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410528520.XA CN104298871B (en) | 2014-10-09 | 2014-10-09 | Method for accelerating voxel human body model dose evaluation based on GPU acceleration in radiation protection |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410528520.XA CN104298871B (en) | 2014-10-09 | 2014-10-09 | Method for accelerating voxel human body model dose evaluation based on GPU acceleration in radiation protection |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104298871A true CN104298871A (en) | 2015-01-21 |
| CN104298871B CN104298871B (en) | 2017-02-08 |
Family
ID=52318594
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410528520.XA Active CN104298871B (en) | 2014-10-09 | 2014-10-09 | Method for accelerating voxel human body model dose evaluation based on GPU acceleration in radiation protection |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104298871B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107290774A (en) * | 2017-07-21 | 2017-10-24 | 四川瑶天纳米科技有限责任公司 | Neutron dose fast determination method |
| CN107832545A (en) * | 2017-11-23 | 2018-03-23 | 哈尔滨工程大学 | A kind of human body external irradiation dose emulation mode based on voxel |
| CN110404184A (en) * | 2019-06-13 | 2019-11-05 | 苏州同调医学科技有限公司 | A kind of method and system of measuring and calculating radiotherapy roentgen dose X distribution and dose objective function |
| CN113420491A (en) * | 2021-08-24 | 2021-09-21 | 中国人民解放军军事科学院军事医学研究院 | Method for evaluating organ radiation dose of experimental animal irradiated outside particles |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101954148A (en) * | 2010-09-15 | 2011-01-26 | 四川大学 | Method for accelerating dosage calculation in radiotherapy based on GPU (Graphics Processing Unit) |
| CN102122323A (en) * | 2011-02-21 | 2011-07-13 | 四川大学 | Method for quickly realizing Gamma analysis method based on GPU (graphic processing unit) |
| CN102201036A (en) * | 2011-05-16 | 2011-09-28 | 四川大学 | Realization of convolution superposition algorithm in graphics processing unit (GPU) acceleration dosage calculation |
| CN103405852A (en) * | 2013-08-29 | 2013-11-27 | 曲桂红 | Rapid calculation method of stereotactic radiotherapy dosage field distribution |
| US20130324784A1 (en) * | 2012-06-01 | 2013-12-05 | Raysearch Laboratories Ab | Method and a system for optimizing a radiation treatment plan |
-
2014
- 2014-10-09 CN CN201410528520.XA patent/CN104298871B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101954148A (en) * | 2010-09-15 | 2011-01-26 | 四川大学 | Method for accelerating dosage calculation in radiotherapy based on GPU (Graphics Processing Unit) |
| CN102122323A (en) * | 2011-02-21 | 2011-07-13 | 四川大学 | Method for quickly realizing Gamma analysis method based on GPU (graphic processing unit) |
| CN102201036A (en) * | 2011-05-16 | 2011-09-28 | 四川大学 | Realization of convolution superposition algorithm in graphics processing unit (GPU) acceleration dosage calculation |
| US20130324784A1 (en) * | 2012-06-01 | 2013-12-05 | Raysearch Laboratories Ab | Method and a system for optimizing a radiation treatment plan |
| CN103405852A (en) * | 2013-08-29 | 2013-11-27 | 曲桂红 | Rapid calculation method of stereotactic radiotherapy dosage field distribution |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107290774A (en) * | 2017-07-21 | 2017-10-24 | 四川瑶天纳米科技有限责任公司 | Neutron dose fast determination method |
| CN107832545A (en) * | 2017-11-23 | 2018-03-23 | 哈尔滨工程大学 | A kind of human body external irradiation dose emulation mode based on voxel |
| CN107832545B (en) * | 2017-11-23 | 2020-11-06 | 哈尔滨工程大学 | Voxel-based human body external irradiation dose simulation method |
| CN110404184A (en) * | 2019-06-13 | 2019-11-05 | 苏州同调医学科技有限公司 | A kind of method and system of measuring and calculating radiotherapy roentgen dose X distribution and dose objective function |
| CN113420491A (en) * | 2021-08-24 | 2021-09-21 | 中国人民解放军军事科学院军事医学研究院 | Method for evaluating organ radiation dose of experimental animal irradiated outside particles |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104298871B (en) | 2017-02-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Qin et al. | Recent developments and comprehensive evaluations of a GPU-based Monte Carlo package for proton therapy | |
| CN104298871B (en) | Method for accelerating voxel human body model dose evaluation based on GPU acceleration in radiation protection | |
| CN101954148A (en) | Method for accelerating dosage calculation in radiotherapy based on GPU (Graphics Processing Unit) | |
| Tse et al. | Efficient reconfigurable design for pricing asian options | |
| CN103116667B (en) | Method of obtaining fusion reactor experimental covering module neutronics parameters | |
| CN103405852A (en) | Rapid calculation method of stereotactic radiotherapy dosage field distribution | |
| Karonis et al. | Distributed and hardware accelerated computing for clinical medical imaging using proton computed tomography (pCT) | |
| CN102201036A (en) | Realization of convolution superposition algorithm in graphics processing unit (GPU) acceleration dosage calculation | |
| CN103065358B (en) | Organ geometry reconstruction method based on image volume element operation | |
| Gladkov et al. | Graphics processing unit based direct simulation Monte Carlo | |
| Aoyama et al. | First application of lattice QCD to Pezy-SC processor | |
| Henderson et al. | A CUDA Monte Carlo simulator for radiation therapy dosimetry based on Geant4 | |
| CN102163319B (en) | Method and system for realization of iterative reconstructed image | |
| Nelson | Improved Convergence Rate of Multi-Group Scattering Moment Tallies for Monte Carlo Neutron Transport Codes | |
| Canal et al. | High energy electromagnetic particle transportation on the GPU | |
| Vlachoudis et al. | Numerically robust geometry engine for compound solid geometries | |
| Murakami et al. | Geant4 based simulation of radiation dosimetry in CUDA | |
| Burke | Kernel Density Estimation Techniques for Monte Carlo Reactor Analysis. | |
| Shukla et al. | Multi-science applications with single codebase-GAMER-for massively parallel architectures | |
| Jevremovic et al. | AGENT code: Open-architecture analysis and configuration of research reactors-Neutron transport modeling with numerical examples | |
| Townson et al. | A hybrid phase-space and histogram source model for GPU-based Monte Carlo radiotherapy dose calculation3 | |
| Zieb | Novel Monte Carlo Methods for Proton Transport Simulations in the Heterogeneous CPU-GPU Environment | |
| Kummerländer | Symbolic generation of GPU-based Lattice Boltzmann implementations | |
| Burke | Development of MGMC: A proxy Multi-Group Monte Carlo Particle Transport Application | |
| Naidoo | HI intensity mapping and cross-correlation science with HIRAX. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |