CN104298871B - 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
- CN104298871B CN104298871B CN201410528520.XA CN201410528520A CN104298871B CN 104298871 B CN104298871 B CN 104298871B CN 201410528520 A CN201410528520 A CN 201410528520A CN 104298871 B CN104298871 B CN 104298871B
- Authority
- CN
- China
- Prior art keywords
- voxel
- dose
- organ
- gpu
- flux
- 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.)
- Active
Links
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 to the voxel manikin dose assessment accelerated method being accelerated based on GPU in a kind of radiation protection, it is
In calculating in voxel manikin radiation protection dosimeter, quickly and accurately determine radiation in irradiated voxel model
3-dimensional dose distribution method.
Background technology
Manikin is radiation proof important tool.Manikin is roughly divided into three by its development：Mass Model, number
Learn model and voxel model.Voxel model is made up of millions of corpusculum units, and voxel model more can reflect human dissection
Structure is more truly, therefore, more accurate based on voxel model Rapid Dose Calculation.
In dose assessment, quickly and accurately determine 3-dimensional dose distribution in body mould for the radiating particle be one basic
Problem.The method of the determination dosage using at present is that the radiation field being calculated based on Monte Carlo method or analytic method, using logical
The method of amount dosage conversion calculates the dosage of each little volume elements.
GPU (Graphic Processing Unit) initially applies the acceleration with figure shows.The single instruction multiple of GPU
According to stream (SIMD：Single Instruction Multiple Data) processing mode can concurrently large-scale data be entered
Row operation, is substantially shorter the calculating time.
GPU is applied in the Rapid Dose Calculation of voxel manikin the present invention.By carrying out the dosage of a large amount of voxels simultaneously
Calculate, its computational accuracy and using CPU (Central Process Unit) computational accuracy quite, its calculating time relatively uses
CPU platform calculates and greatly reduces.
Realize the parallel computation based on GPU, understand support GPU computing first becomes instrument.CUDA(Computer
Unified Device Architecture) it is that NVIDIA company becomes a brand-new software and hardware architecture of proposition for GPU,
It writes out the program of execution on GPU based on C language.May refer to CUDA Programming Guide with regard to CUDA detailed description
And CUDA reference manual.The present invention is formally accelerated to voxel manikin dosage using CUDA technology in GPU programming.
Content of the invention
The technology solve problem of the present invention：Overcome the deficiencies in the prior art, provide and in a kind of radiation protection, be based on GPU acceleration
Voxel manikin dose assessment accelerated method, can quickly and accurately determine three-dimensional agent in voxel model for the radiating particle
Amount distribution.
The technical scheme is that：Added based on the voxel manikin dose assessment that GPU accelerates in a kind of radiation protection
Fast method it is characterised in that：Employ CPU and combine with GPU and carry out Rapid Dose Calculation, for some, there is the process of parallel condition
Journey carries out acceleration process using GPU, and the serial computing such as complex logic and issued transaction still adopts CPU.Realize step to include：
(1) voxel manikin parameter is formed on CPU platform；
(2) select the GPU for determining voxel dose, be GPU distribution memory space, the body in GPU transfer step (1)
The calculating data of plain manikin parameter, need transmission data include voxel sequence, organ name and ID and U corresponding relation and
Voxel precision；
(3) utilize the dosage of GPU parallel computation voxel each voxel of manikin, dose assessment is in voxel model parameter
On the basis of, using flux dosage conversion method, the i.e. flux of known free space, by conversion factor, calculate often individual respectively
The absorbed dose of radiation of element, using the dosage of GPU each voxel of parallel computation, each GPU thread calculates the dosage of a voxel, finally
All voxel doses are averagely Absorbed Doses.
Present invention advantage compared with prior art is：The present invention combines the serial computing such as CPU complex logic and things,
Voxel equivalent dose in being calculated using graphic process unit (Graphic Processor Unit, GPU) treatment dosage, this process is deposited
Data volume big, data dependence is low, data has identical configuration processor, degree of parallelism and calculate density high the features such as.Thus it is logical
Cross GPU, can allow more computing unit executed in parallel, shorten the calculating time, improve computational efficiency, realize voxel manikin
Dosage calculates effect in real time.
Brief description
Fig. 1 is the flowchart of the present invention.
Specific embodiment
As shown in figure 1, the present invention to implement step as follows：
(1) voxel manikin parameter is formed on CPU platform, voxel manikin parameter includes voxel sequence, organ
Name and ID and U corresponding relation and voxel precision；
(2) select the GPU for determining voxel dose, be GPU distribution memory space, the body in GPU transfer step (1)
The calculating data of plain manikin parameter, need transmission data include voxel sequence, organ name and ID and U corresponding relation and
Voxel precision；
(3) utilize the dosage of GPU parallel computation voxel each voxel of manikin.For the voxel of organ discretization, allow
Voxel mark and thread identification correspond, the product of the flux of voxel dose position and flux dose conversion factor.Dosage tires out
Plus during by voxel mark with thread identification correspond, search the corresponding dosage of organ voxel, by each voxel of organ
Dosage carries out cumulative average again.
The Rapid Dose Calculation accuracy of voxel model is, the human organ of voxel manikin is voxelization, by a large amount of
Little cuboid voxel composition, each voxel comprises organ code name and material composition.Dose assessment is on the basis of this model, point
Do not calculate the absorbed dose of radiation of each voxel, using the dosage of GPU each voxel of parallel computation, finally will be average for all voxel doses
It is Absorbed Doses.Can be expressed from the next.
Wherein,For tissue T absorbed dose of radiation, DT,jFor the absorbed dose of radiation of this organ voxel j, N is the body that this organ comprises
Prime number mesh.
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, the flux of this grid is expressed as the flux of this voxel.
The absorbed dose of radiation of organ is the average of the flux forming each voxel location, with conversion under this energy for this organ
The result of factor product.Therefore, Absorbed Doses calculate and translate on the flux obtaining this each voxel of organ.
The voxel of voxel model is according to certain aligned transfer, and voxel fill order is first X-axis, then Y-axis, last Z
Axle.Fill order according to voxel to travel through to all voxel number consecutivelies, the corresponding organ of each voxel U, represents it
For forming the voxel of this organ.
For the parallel data processing function using GPU, accelerate calculating speed, distribute one to each voxel of voxel model
Thread is calculated.When calculating voxel dose, obtain voxel location first, according to position enquiring flux, using flux dosage
Conversion factor, obtains voxel dose.
DT,j=φj(x,y,z)*CF (2)
Wherein, φj(x, y, z) is the flux of (x, y, z) at voxel j, and CF is flux dose conversion factor.
In dosage cumulative process, voxel mark and thread identification are corresponded, to the voxel in same organ, parallel
Change out voxel dose, then all voxel doses add up to organ.
Wherein, DT, andFor tissue T dosage value preset.
Claims (2)
1. in a kind of radiation protection based on GPU accelerate voxel manikin dose assessment accelerated method it is characterised in that：Use
CPU is combined with GPU and carries out Rapid Dose Calculation, and the processing procedure having parallel condition for some carries out acceleration process using GPU,
And complex logic and issued transaction serial computing still adopt CPU；Realize step to include：
(1) voxel manikin parameter is formed on CPU platform；Voxel manikin parameter include voxel sequence, organ name and
ID and U corresponding relation and voxel precision；(2) select the GPU for determining voxel dose, be GPU distribution memory space, to
The calculating data of the voxel manikin parameter in GPU transfer step (1), needs the data of transmission to include voxel sequence, organ
Name and ID and U corresponding relation and voxel precision；
(3) utilize the dosage of GPU parallel computation voxel each voxel of manikin, for the voxel of organ discretization, allow voxel
Mark is corresponded with thread identification, the product of the flux of voxel dose position and flux dose conversion factor；Dosage added up
In journey, voxel mark and thread identification are corresponded, search the corresponding dosage of organ voxel, by the dosage of each voxel of organ
Carry out cumulative average again；
The Rapid Dose Calculation accuracy of voxel model is, the human organ of voxel manikin is voxelization, by substantial amounts of little
Cuboid voxel forms, and each voxel comprises organ code name and material composition；Dose assessment is on the basis of this model, counts respectively
Calculate the absorbed dose of radiation of each voxel, using the dosage of GPU each voxel of parallel computation, finally all voxel doses are averagely
Absorbed Doses, are expressed from the next：
Wherein,For tissue T absorbed dose of radiation, DT,jFor the absorbed dose of radiation of this organ voxel j, N is the number of voxels that this organ comprises
Mesh；
The acquisition of the voxel flux of voxel manikin is the centre coordinate by calculating voxel, then searches from radiation field grid
Go out the grid of arest neighbors, the flux of this grid is expressed as the flux of this voxel；
The absorbed dose of radiation of organ is the average of the flux forming each voxel location, takes advantage of with conversion factor under flux for this organ
Long-pending result；Absorbed Doses calculate and translate on the flux obtaining this each voxel of organ；
The voxel of voxel model is according to certain aligned transfer, and voxel fill order is first X-axis, then Y-axis, last Z axis, presses
Fill order according to voxel to travel through to all voxel number consecutivelies, the corresponding organ of each voxel U, represents that it is composition
The voxel of this organ；
For the parallel data processing function using GPU, accelerate calculating speed, distribute a thread to each voxel of voxel model
Calculated, when calculating voxel dose, obtained voxel location first, according to position enquiring flux, using the conversion of flux dosage
The factor, obtains voxel dose：
DT,j=φj(x,y,z)*CF (2)
Wherein, φj(x, y, z) is the flux of (x, y, z) at voxel j, and CF is flux dose conversion factor；
In dosage cumulative process, voxel mark and thread identification are corresponded, to the voxel in same organ, Parallel transformation
Go out voxel dose, then all voxel doses add up to organ,
Wherein, DT, andFor tissue T dosage value preset.
2. accelerated based on the voxel manikin dose assessment that GPU accelerates in a kind of radiation protection according to claim 1
Method it is characterised in that：Dosage includes photon doses, also includes neutron dose, and photon doses are used for evaluating photon irradiation to people
The harm of body, neutron dose is used for evaluating the harm to human body for the neutron exposure.
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 CN104298871A (en) | 2015-01-21 |
| CN104298871B true 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) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107290774A (en) * | 2017-07-21 | 2017-10-24 | 四川瑶天纳米科技有限责任公司 | Neutron dose fast determination method |
| 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 |
| CN113420491B (en) * | 2021-08-24 | 2021-12-17 | 中国人民解放军军事科学院军事医学研究院 | Method for evaluating organ radiation dose of experimental animal irradiated outside particles |
Family Cites Families (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 |
| US9468776B2 (en) * | 2012-06-01 | 2016-10-18 | Raysearch Laboratories Ab | Method and a system for optimizing a radiation treatment plan based on a reference dose distribution |
| CN103405852A (en) * | 2013-08-29 | 2013-11-27 | 曲桂红 | Rapid calculation method of stereotactic radiotherapy dosage field distribution |
-
2014
- 2014-10-09 CN CN201410528520.XA patent/CN104298871B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104298871A (en) | 2015-01-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jia et al. | Development of a GPU-based Monte Carlo dose calculation code for coupled electron–photon transport | |
| CN104298871B (en) | Method for accelerating voxel human body model dose evaluation based on GPU acceleration in radiation protection | |
| Yepes et al. | A GPU implementation of a track-repeating algorithm for proton radiotherapy dose calculations | |
| Douglass et al. | Development of a randomized 3D cell model for Monte Carlo microdosimetry simulations | |
| Kanzaki | Monte Carlo integration on GPU | |
| Karonis et al. | Distributed and hardware accelerated computing for clinical medical imaging using proton computed tomography (pCT) | |
| CN104182571B (en) | Kriging interpolation methods based on Delaunay and GPU | |
| 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 | |
| CN102201036A (en) | Realization of convolution superposition algorithm in graphics processing unit (GPU) acceleration dosage calculation | |
| Arter et al. | The SMARDDA approach to ray tracing and particle tracking | |
| CN103065358B (en) | Organ geometry reconstruction method based on image volume element operation | |
| Yegin | A new approach to geometry modeling for Monte Carlo particle transport: an application to the EGS code system | |
| CN107290774A (en) | Neutron dose fast determination method | |
| Henderson et al. | A CUDA Monte Carlo simulator for radiation therapy dosimetry based on Geant4 | |
| Pasmann et al. | iQMC: Iterative Quasi-Monte Carlo for k-Eigenvalue Neutron Transport Simulations | |
| 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. | |
| Bush et al. | A technique for generating phase-space-based Monte Carlo beamlets in radiotherapy applications | |
| Kapusta et al. | Application of GPU parallel computing for acceleration of finite element method based 3D reconstruction algorithms in electrical capacitance tomography | |
| Orok | Implementation of a tetrahedral mesh phantom geometry library for EGSnrc | |
| CN116612853B (en) | Radiotherapy verification plan dose generation method, radiotherapy verification plan dose generation system, electronic equipment and storage medium | |
| Zhukovskiy et al. | Hybrid parallelization of computing the electron fluxes produced by photon radiation |
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 |