1,156 research outputs found
GPU-based ultra-fast direct aperture optimization for online adaptive radiation therapy
Online adaptive radiation therapy (ART) has great promise to significantly
reduce normal tissue toxicity and/or improve tumor control through real-time
treatment adaptations based on the current patient anatomy. However, the major
technical obstacle for clinical realization of online ART, namely the inability
to achieve real-time efficiency in treatment re-planning, has yet to be solved.
To overcome this challenge, this paper presents our work on the implementation
of an intensity modulated radiation therapy (IMRT) direct aperture optimization
(DAO) algorithm on graphics processing unit (GPU) based on our previous work on
CPU. We formulate the DAO problem as a large-scale convex programming problem,
and use an exact method called column generation approach to deal with its
extremely large dimensionality on GPU. Five 9-field prostate and five 5-field
head-and-neck IMRT clinical cases with 5\times5 mm2 beamlet size and
2.5\times2.5\times2.5 mm3 voxel size were used to evaluate our algorithm on
GPU. It takes only 0.7~2.5 seconds for our implementation to generate optimal
treatment plans using 50 MLC apertures on an NVIDIA Tesla C1060 GPU card. Our
work has therefore solved a major problem in developing ultra-fast
(re-)planning technologies for online ART
Fast Monte Carlo Simulation for Patient-specific CT/CBCT Imaging Dose Calculation
Recently, X-ray imaging dose from computed tomography (CT) or cone beam CT
(CBCT) scans has become a serious concern. Patient-specific imaging dose
calculation has been proposed for the purpose of dose management. While Monte
Carlo (MC) dose calculation can be quite accurate for this purpose, it suffers
from low computational efficiency. In response to this problem, we have
successfully developed a MC dose calculation package, gCTD, on GPU architecture
under the NVIDIA CUDA platform for fast and accurate estimation of the x-ray
imaging dose received by a patient during a CT or CBCT scan. Techniques have
been developed particularly for the GPU architecture to achieve high
computational efficiency. Dose calculations using CBCT scanning geometry in a
homogeneous water phantom and a heterogeneous Zubal head phantom have shown
good agreement between gCTD and EGSnrc, indicating the accuracy of our code. In
terms of improved efficiency, it is found that gCTD attains a speed-up of ~400
times in the homogeneous water phantom and ~76.6 times in the Zubal phantom
compared to EGSnrc. As for absolute computation time, imaging dose calculation
for the Zubal phantom can be accomplished in ~17 sec with the average relative
standard deviation of 0.4%. Though our gCTD code has been developed and tested
in the context of CBCT scans, with simple modification of geometry it can be
used for assessing imaging dose in CT scans as well.Comment: 18 pages, 7 figures, and 1 tabl
Beam Orientation Optimization for Intensity Modulated Radiation Therapy using Adaptive l1 Minimization
Beam orientation optimization (BOO) is a key component in the process of IMRT
treatment planning. It determines to what degree one can achieve a good
treatment plan quality in the subsequent plan optimization process. In this
paper, we have developed a BOO algorithm via adaptive l_1 minimization.
Specifically, we introduce a sparsity energy function term into our model which
contains weighting factors for each beam angle adaptively adjusted during the
optimization process. Such an energy term favors small number of beam angles.
By optimizing a total energy function containing a dosimetric term and the
sparsity term, we are able to identify the unimportant beam angles and
gradually remove them without largely sacrificing the dosimetric objective. In
one typical prostate case, the convergence property of our algorithm, as well
as the how the beam angles are selected during the optimization process, is
demonstrated. Fluence map optimization (FMO) is then performed based on the
optimized beam angles. The resulted plan quality is presented and found to be
better than that obtained from unoptimized (equiangular) beam orientations. We
have further systematically validated our algorithm in the contexts of 5-9
coplanar beams for 5 prostate cases and 1 head and neck case. For each case,
the final FMO objective function value is used to compare the optimized beam
orientations and the equiangular ones. It is found that, our BOO algorithm can
lead to beam configurations which attain lower FMO objective function values
than corresponding equiangular cases, indicating the effectiveness of our BOO
algorithm.Comment: 19 pages, 2 tables, and 5 figure
GPU-based Iterative Cone Beam CT Reconstruction Using Tight Frame Regularization
X-ray imaging dose from serial cone-beam CT (CBCT) scans raises a clinical
concern in most image guided radiation therapy procedures. It is the goal of
this paper to develop a fast GPU-based algorithm to reconstruct high quality
CBCT images from undersampled and noisy projection data so as to lower the
imaging dose. For this purpose, we have developed an iterative tight frame (TF)
based CBCT reconstruction algorithm. A condition that a real CBCT image has a
sparse representation under a TF basis is imposed in the iteration process as
regularization to the solution. To speed up the computation, a multi-grid
method is employed. Our GPU implementation has achieved high computational
efficiency and a CBCT image of resolution 512\times512\times70 can be
reconstructed in ~5 min. We have tested our algorithm on a digital NCAT phantom
and a physical Catphan phantom. It is found that our TF-based algorithm is able
to reconstrct CBCT in the context of undersampling and low mAs levels. We have
also quantitatively analyzed the reconstructed CBCT image quality in terms of
modulation-transfer-function and contrast-to-noise ratio under various scanning
conditions. The results confirm the high CBCT image quality obtained from our
TF algorithm. Moreover, our algorithm has also been validated in a real
clinical context using a head-and-neck patient case. Comparisons of the
developed TF algorithm and the current state-of-the-art TV algorithm have also
been made in various cases studied in terms of reconstructed image quality and
computation efficiency.Comment: 24 pages, 8 figures, accepted by Phys. Med. Bio
Effect of long-term fluorination on surface electrical performance of ethylene propylene rubber
To investigate the effect of fluorination on surface electrical performance of ethylene propylene rubber (EPR), four pieces of EPR specimens are prepared and fluorinated for different duration ranging from 120 min to 480 min. The surface morphology and element compositions of experimental specimens are tested. The surface potential decay and complex permittivity are measured. The tracking discharge property and the erosion properties after tracking test are investigated. The surface charge transport and electrical property tailoring mechanism are discussed. The results show that when the fluorination duration is from 120 min to 360 min, the surface morphology gradually becomes flat and compact with time, and the surface resistance to electrical tracking is improved. While an excessive fluorination with the duration of 480 min leads to a significant increase in surface conductivity and a weakening in the resistance to electrical tracking, which is mainly due to a change in surface morphology. The recommended fluorination duration for an optimised surface property is between 240 and 360 min for EPR specimen
GPU-based Fast Cone Beam CT Reconstruction from Undersampled and Noisy Projection Data via Total Variation
Purpose: Cone-beam CT (CBCT) plays an important role in image guided
radiation therapy (IGRT). However, the large radiation dose from serial CBCT
scans in most IGRT procedures raises a clinical concern, especially for
pediatric patients who are essentially excluded from receiving IGRT for this
reason. The goal of this work is to develop a fast GPU-based algorithm to
reconstruct CBCT from undersampled and noisy projection data so as to lower the
imaging dose. Methods: The CBCT is reconstructed by minimizing an energy
functional consisting of a data fidelity term and a total variation
regularization term. We developed a GPU-friendly version of the
forward-backward splitting algorithm to solve this model. A multi-grid
technique is also employed. Results: It is found that 20~40 x-ray projections
are sufficient to reconstruct images with satisfactory quality for IGRT. The
reconstruction time ranges from 77 to 130 sec on a NVIDIA Tesla C1060 GPU card,
depending on the number of projections used, which is estimated about 100 times
faster than similar iterative reconstruction approaches. Moreover, phantom
studies indicate that our algorithm enables the CBCT to be reconstructed under
a scanning protocol with as low as 0.1 mAs/projection. Comparing with currently
widely used full-fan head and neck scanning protocol of ~360 projections with
0.4 mAs/projection, it is estimated that an overall 36~72 times dose reduction
has been achieved in our fast CBCT reconstruction algorithm. Conclusions: This
work indicates that the developed GPU-based CBCT reconstruction algorithm is
capable of lowering imaging dose considerably. The high computation efficiency
in this algorithm makes the iterative CBCT reconstruction approach applicable
in real clinical environments.Comment: Accepted as a letter in Med. Phys., brief clarifying comments and
updated references. 6 pages and 2 figure
Development of a GPU-based Monte Carlo dose calculation code for coupled electron-photon transport
Monte Carlo simulation is the most accurate method for absorbed dose
calculations in radiotherapy. Its efficiency still requires improvement for
routine clinical applications, especially for online adaptive radiotherapy. In
this paper, we report our recent development on a GPU-based Monte Carlo dose
calculation code for coupled electron-photon transport. We have implemented the
Dose Planning Method (DPM) Monte Carlo dose calculation package (Sempau et al,
Phys. Med. Biol., 45(2000)2263-2291) on GPU architecture under CUDA platform.
The implementation has been tested with respect to the original sequential DPM
code on CPU in phantoms with water-lung-water or water-bone-water slab
geometry. A 20 MeV mono-energetic electron point source or a 6 MV photon point
source is used in our validation. The results demonstrate adequate accuracy of
our GPU implementation for both electron and photon beams in radiotherapy
energy range. Speed up factors of about 5.0 ~ 6.6 times have been observed,
using an NVIDIA Tesla C1060 GPU card against a 2.27GHz Intel Xeon CPU
processor.Comment: 13 pages, 3 figures, and 1 table. Paper revised. Figures update
An efficient imaging algorithm for GNSS-R bi-static SAR
Global Navigation Satellite System Reflectometry (GNSS-R) based Bi-static Synthetic Aperture Radar (BSAR) is becoming more and more important in remote sensing, given its low power, low mass, low cost, and real-time global coverage capability. Due to its complex configuration, the imaging for GNSS-R BSAR is usually based on the Back-Projection Algorithm (BPA), which is very time consuming. In this paper, an efficient and general imaging algorithm for GNSS-R BSAR is presented. A Two Step Range Cell Migration (TSRCM) correction is firstly applied. The first step roughly compensates the RCM and Doppler phase caused by the motion of the transmitter, which simplifies the SAR data into the quasi-mono-static case. The second step removes the residual RCM caused by the motion of the receiver using the modified frequency scaling algorithm. Then, a cubic phase perturbation operation is introduced to equalize the Doppler frequency modulation rate along the same range cell. Finally, azimuth phase compensation and geometric correction are completed to obtain the focused SAR image. A simulation and experiment are conducted to demonstrate the feasibility of the proposed algorithm, showing that the proposed algorithm is more efficient than the BPA, without causing significant degradation in imaging quality
Orbit optimization for ASTROD-GW and its time delay interferometry with two arms using CGC ephemeris
ASTROD-GW (ASTROD [Astrodynamical Space Test of Relativity using Optical
Devices] optimized for Gravitation Wave detection) is an optimization of ASTROD
to focus on the goal of detection of gravitation waves. The detection
sensitivity is shifted 52 times toward larger wavelength compared to that of
LISA. The mission orbits of the 3 spacecraft forming a nearly equilateral
triangular array are chosen to be near the Sun-Earth Lagrange points L3, L4 and
L5. The 3 spacecraft range interferometrically with one another with arm length
about 260 million kilometers. In order to attain the requisite sensitivity for
ASTROD-GW, laser frequency noise must be suppressed below the secondary noises
such as the optical path noise, acceleration noise etc. For suppressing laser
frequency noise, we need to use time delay interferometry (TDI) to match the
two different optical paths (times of travel). Since planets and other
solar-system bodies perturb the orbits of ASTROD-GW spacecraft and affect the
(TDI), we simulate the time delay numerically using CGC 2.7 ephemeris
framework. To conform to the ASTROD-GW planning, we work out a set of 20-year
optimized mission orbits of ASTROD-GW spacecraft starting at June 21, 2028, and
calculate the residual optical path differences in the first and second
generation TDI for one-detector case. In our optimized mission orbits for 20
years, changes of arm length are less than 0.0003 AU; the relative Doppler
velocities are less than 3m/s. All the second generation TDI for one-detector
case satisfies the ASTROD-GW requirement.Comment: 17 pages, 7 figures, 1 tabl
Continuous twin screw rheo-extrusion of an AZ91D magnesium alloy
© The Minerals, Metals & Materials Society and ASM International 2012The twin screw rheo-extrusion (TSRE) is designed to take advantage of the nondendritc microstructure and thixotropic characterization of semisolid-metal slurries and produce simple metal profiles directly from melts. The extrusion equipment consists of a rotor-stator high shear slurry maker, a twin screw extruder, and a die assembly. The process is continuous and has a potential for significantly saving energy, manufacturing cost, and enhancing efficiency. The present investigation was carried out to study the process performance for processing rods of an AZ91D magnesium alloy and the microstructure evolution during processing. The semisolid slurry prepared by the process was characterized by uniformly distributed nondendritic granular primary phase particles. AZ91D rods with uniform and fine microstructures and moderate mechanical properties were produced. For the given slurry making parameters, decreasing extrusion temperature was found to improve microstructures and properties. The mechanisms of particle granulation and refinement and the effect of processing parameters on process performance and thermal management are discussed. © 2012 The Minerals, Metals & Materials Society and ASM International.EPSRC (UK) and Rautomead Lt
- …