Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Tracking”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Computerized analysis of coronary artery disease: Performance evaluation of segmentation and tracking of coronary arteries in CT angiograms

Purpose: The authors are developing a computer-aided detection system to assist radiologists in analysis of coronary artery disease in coronary CT angiograms (cCTA). This study evaluated the accuracy of the authors’ coronary artery segmentation and tracking method which are the essential steps to define the search space for the detection of atherosclerotic plaques. Methods: The heart region in cCTA is segmented and the vascular structures are enhanced using the authors’ multiscale coronary artery response (MSCAR) method that performed 3D multiscale filtering and analysis of the eigenvalues of Hessian matrices. Starting from seed points at the origins of the left and right coronary arteries, a 3D rolling balloon region growing (RBG) method that adapts to the local vessel size segmented and tracked each of the coronary arteries and identifies the branches along the tracked vessels. The branches are queued and subsequently tracked until the queue is exhausted. With Institutional Review Board approval, 62 cCTA were collected retrospectively from the authors’ patient files. Three experienced cardiothoracic radiologists manually tracked and marked center points of the coronary arteries as reference standard following the 17-segment model that includes clinically significant coronary arteries. Two radiologists visually examined the computer-segmented vessels and marked the mistakenly tracked veins and noisy structures as false positives (FPs). For the 62 cases, the radiologists marked a total of 10191 center points on 865 visible coronary artery segments. Results: The computer-segmented vessels overlapped with 83.6% (8520/10191) of the center points. Relative to the 865 radiologist-marked segments, the sensitivity reached 91.9% (795/865) if a true positive is defined as a computer-segmented vessel that overlapped with at least 10% of the reference center points marked on the segment. When the overlap threshold is increased to 50% and 100%, the sensitivities were 86.2% and 53.4%, respectively. For the 62 test cases, a total of 55 FPs were identified by radiologist in 23 of the cases. Conclusions: The authors’ MSCAR-RBG method achieved high sensitivity for coronary artery segmentation and tracking. Studies are underway to further improve the accuracy for the arterial segments affected by motion artifacts, severe calcified and noncalcified soft plaques, and to reduce the false tracking of the veins and other noisy structures. Methods are also being developed to detect coronary artery disease along the tracked vessels.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

DMLC tracking and gating can improve dose coverage for prostate VMAT

Purpose: To assess and compare the dosimetric impact of dynamic multileaf collimator (DMLC) tracking and gating as motion correction strategies to account for intrafraction motion during conventionally fractionated prostate radiotherapy. Methods: A dose reconstruction method was used to retrospectively assess the dose distributions delivered without motion correction during volumetric modulated arc therapy fractions for 20 fractions of five prostate cancer patients who received conventionally fractionated radiotherapy. These delivered dose distributions were compared with the dose distributions which would have been delivered had DMLC tracking or gating motion correction strategies been implemented. The delivered dose distributions were constructed by incorporating the observed prostate motion with the patient's original treatment plan to simulate the treatment delivery. The DMLC tracking dose distributions were constructed using the same dose reconstruction method with the addition of MLC positions from Linac log files obtained during DMLC tracking simulations with the observed prostate motions input to the DMLC tracking software. The gating dose distributions were constructed by altering the prostate motion to simulate the application of a gating threshold of 3 mm for 5 s. Results: The delivered dose distributions showed that dosimetric effects of intrafraction prostate motion could be substantial for some fractions, with an estimated dose decrease of more than 19% and 34% from the planned CTVD{sub 99%} and PTV D{sub 95%} values, respectively, for one fraction. Evaluation of dose distributions for DMLC tracking and gating deliveries showed that both interventions were effective in improving the CTV D{sub 99%} for all of the selected fractions to within 4% of planned value for all fractions. For the delivered dose distributions the difference in rectum V{sub 65%} for the individual fractions from planned ranged from −44% to 101% and for the bladder V{sub 65%} the range was −61% to 26% from planned. The application of tracking decreased the maximum rectum and bladder V{sub 65%} difference to 6% and 4%, respectively. Conclusions: For the first time, the dosimetric impact of DMLC tracking and gating to account for intrafraction motion during prostate radiotherapy has been assessed and compared with no motion correction. Without motion correction intrafraction prostate motion can result in a significant decrease in target dose coverage for a small number of individual fractions. This is unlikely to effect the overall treatment for most patients undergoing conventionally fractionated treatments. Both DMLC tracking and gating demonstrate dose distributions for all assessed fractions that are robust to intrafraction motion.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Markerless tracking of small lung tumors for stereotactic radiotherapy

Purpose: (1) To validate retrospective markerless tracking software for small lung tumors by comparing tracked motion in 4-dimensional planning computed tomography (4DCT) derived kV projection images and known tumor motion in the same 4DCT. (2) To evaluate variability of tumor motion using kV projection images from cone-beam computed tomography (CBCT) scans acquired on different days. Methods: Nonclinical tumor tracking software (TTS) used a normalized cross correlation algorithm to track the tumor on enhanced kV projection images (e.g., from a CBCT scan). The reference dataset consisted of digitally reconstructed radiographs (DRRs) from one phase of a planning 4DCT. TTS matches two in-plane coordinates and obtains the out-of-plane coordinate by triangulating with match results from other projections. (1) To validate TTS, tracking results were compared with known 4DCT tumor motion for two patients (A and B). Projection images (1 image/1°) were digitally reconstructed for each 4DCT phase. From these, kV projection series were composed simulating full breathing cycles every 20° of gantry rotation [breathing period = 20°/(6°/s) = 3.33 s]. Reference templates were 360 “tumor enhanced” DRRs from the 4DCT expiration phase. TTS-derived tumor motion was compared to known tumor motion on 4DCT. (2) For five patients, TTS-assessed motion during clinical CBCT acquisition was compared with motion on the planning 4DCT, and the motion component in the Y (cranio–caudal)-direction was compared with the motion of an external marker box (RPM, real-time position management). Results: (1) Validation results: TTS for case A (tumor 6.2 cm{sup 3}, 32 mm axial diameter) over 360° showed mean motion X (medial–lateral) = 3.4, Y = 11.5, and Z (ventral–dorsal) = 4.9 mm (1 SD < 1.0 mm). Corresponding 4DCT motion was X = 3.1, Y = 11.3, and Z = 5.1 mm. Correlation coefficients between TTS tumor motion and displacement of the tumor’s center of mass (CoM) on 4DCT were 0.64, 0.96, and 0.82 (X, Y, and Z, respectively). For case B (4.1 cm{sup 3}, 20 mm diameter), due to temporarily decreased tumor visibility preventing TTS from resolving the tumor, robust tracking data were only available between angles 300°–40° and 120°–220°. Mean motion according to TTS was X = 2.0, Y = 7.7, and Z = 8.2 mm (1 SD < 0.9 mm). Tumor motion on 4DCT was X = 1.8, Y = 7.6, and Z = 9.5 mm and correlation coefficients between TTS motion and CoM displacement were 0.59, 0.95, and 0.93 (X, Y, and Z, respectively). (2) Clinical results: TTS revealed a mean intrafraction variation in tumor motion in Y-direction of >2.0 mm (1 SD) in four of five patients. In addition, clinical tumor motion amplitude differed from that seen on planning 4DCT. Internal and external structures that create abrupt density change (e.g., table-top edge, interface between lung/mediastinum and lung/heart) were observed to prevent 360° tracking of the tumor. Correlation coefficients between TTS motion in the Y-direction and the RPM signal (22 observations) ranged from 0.78 to 0.96. In 2D, 241 TTS matches at end-inspiration and end-expiration were visually validated: mean difference was 0.8 mm (SD = 0.7) for both. Conclusions: TTS can track small lung tumors if these are visible in kV projections. A 4DCT dataset can be used to validate kV tracking of moving targets. TTS and 4DCT displacement agreed to within 2 mm. TTS and RPM motion were closely associated but tumor motion during CBCT can vary from the planning 4DCT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-JeP1-12: Head-To-Head Performance Characterization of Two Multileaf Collimator Tracking Algorithms for Radiotherapy

Purpose: Multi-leaf collimator (MLC) tracking is being clinically pioneered to continuously compensate for thoracic and abdominal motion during radiotherapy. The purpose of this work is to characterize the performance of two MLC tracking algorithms for cancer radiotherapy, based on a direct optimization and a piecewise leaf fitting approach respectively. Methods: To test the algorithms, both physical and in silico experiments were performed. Previously published high and low modulation VMAT plans for lung and prostate cancer cases were used along with eight patient-measured organ-specific trajectories. For both MLC tracking algorithm, the plans were run with their corresponding patient trajectories. The physical experiments were performed on a Trilogy Varian linac and a programmable phantom (HexaMotion platform). For each MLC tracking algorithm, plan and patient trajectory, the tracking accuracy was quantified as the difference in aperture area between ideal and fitted MLC. To compare algorithms, the average cumulative tracking error area for each experiment was calculated. The two-sample Kolmogorov-Smirnov (KS) test was used to evaluate the cumulative tracking errors between algorithms. Results: Comparison of tracking errors for the physical and in silico experiments showed minor differences between the two algorithms. The KS D-statistics for the physical experiments were below 0.05 denoting no significant differences between the two distributions pattern and the average error area (direct optimization/piecewise leaf-fitting) were comparable (66.64 cm2/65.65 cm2). For the in silico experiments, the KS D-statistics were below 0.05 and the average errors area were also equivalent (49.38 cm2/48.98 cm2). Conclusion: The comparison between the two leaf fittings algorithms demonstrated no significant differences in tracking errors, neither in a clinically realistic environment nor in silico. The similarities in the two independent algorithms give confidence in the use of either algorithm for clinical implementation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Improving ICARUS Track Reconstruction Algorithms

The ICARUS experiment is part of the Short-Baseline Neutrino (SBN) program at Fermilab. The main goal of the experiment is to investigate the possibility of sterile neutrinos in the O(1 eV) mass region and provide clarification of the anomaly detected from the Liquid Scintillator Neutrino Detector (LSND) and MiniBooNE experiments. The ICARUS-T600 detector is a Liquid Argon Time Projection Chamber (LAr-TPC), that can provide excellent 3D imaging and calorimetric reconstruction of any ionizing particles. This detection technique allows a detailed study of neutrino interactions, spanning a wide energy spectrum (from a few keV to several hundreds of GeV). The detector consists of two identical adjacent modules, filled with a total of 760 tons of ultra-pure liquid argon. Each module houses two LAr-TPCs separated by a common cathode with a maximum drift distance of 1.5 m, equivalent to about 1 ms drift time for the nominal $500$ V/m electric drift field. The anode is made of three parallel wire planes positioned 3 mm apart, where the stainless-steel wires are oriented on each plane at a different angle with respect to the horizontal direction ($+60^\degree$,$-60^\degree$,$0^\degree$). The first two planes (Induction 1 and Induction 2) provide a non-destructive charge measurement, whereas the ionization charge is fully collected by the last collection plane. In total, 53248 wires with a 3 mm pitch and length up to 9 m are installed in the detector. In the first stage of the reconstruction, segments of waveforms corresponding to physical signals (hits) are searched for in the deconvolved wire waveform with a threshold-based hit-finding algorithm. Each hit is then fitted with a Gaussian, whose area is proportional to the number of drift electrons generating the signal. In the second stage of the reconstruction, hits are passed as input to Pandora, a framework software composed of different pattern recognition algorithms, that performs a 3D reconstruction of the full image recorded in the collected event, including the identification of interaction vertices and tracks and showers inside the TPC. These are organized into a hierarchical structure (called slice) of particles generated starting from a primary interaction vertex. In some cases, related to the inefficiencies in the hit detection or excessive deflection of the particle trajectory, Pandora breaks the particle's track into two or more smaller pieces and considers each piece as an independent track. We studied this phenomenon focusing on primary muons from ν_μ CC interactions contained in a single module with a track at least 20 cm long, to exclude delta rays. The study determined that about $7-8\%$ of the muon tracks are broken. Approximately $80\%$ of the times, Pandora assigns all segments of the track to the same slice (intra-slice track split), while in the remaining $20\%$ of the cases, one of the segments is associated with another slice (extra-slice track split). To mitigate this phenomenon, we designed an algorithm that detects and stitches the tracks broken by Pandora for the intra-slice split. In Monte Carlo simulations, the algorithm showed an efficiency exceeding $80\%$ and a purity exceeding $93\%$.

Ricci, Alessandro Maria [Pisa U.; INFN, Pisa] (ORC↗
NASA NTRS1972

Satellite-satellite tracking for estimating geopotential coefficients.

Results of simulation studies of tracking of low-altitude satellites by a high-altitude satellite using range measurements between satellites. The high-altitude satellite is tracked by ground-based range trackers. The low-altitude satellites are tracked only by the high-altitude satellites. Orbital elements for one tracking satellite and four tracked low-altitude satellites are estimated, along with a set of low degree and order geopotential coefficients. Estimation of geopotential coefficients of low degree and order is shown to be considerably more accurate using satellite-to-satellite tracking than using ground-based tracking. The use of data from one high-altitude satellite tracking four low-altitude satellites in a range of inclinations is shown to produce geopotential coefficient estimates with smaller errors than presently exist in individual coefficients or would be obtained using ground-based tracking of the same satellites.

Martin, C. F.↗
NASA NTRS1982

Fission track studies of xenolithic chondrites - Implications regarding brecciation and metamorphism

Fission tracks in phosphates from one gas-poor chondrite and three gas-rich ones are studied to determine their thermal history and brecciation time scales. Model fission track ages are calculated for given track densities due to Pu-244 and U-238; track densities and uranium measurements for whitlockite are shown, and possible sources of the tracks are mentioned. Details of track density and uranium measurements are discussed for each meteorite separately. Whitlockites from all the meteorites give model fission track ages of 4.4 Gyr assuming a Pu/U ratio at 4.55 Gyr of 0.045. The final brecciation event definitely did not reset the track clock in phosphates of one meteor and probably not in another two. It is concluded that the observed fission track ages date the end of metamorphic cooling in the meteorite parent bodies and support the planetesimal model for the formation of xenolithic chondrites.

Kothari, B. K.↗
NASA NTRS1990

The extended tracking network and indications of baseline precision and accuracy in the North Andes

The CASA Uno Global Positioning System (GPS) experiment (January-February 1988) included an extended tracking network which covered three continents in addition to the network of scientific interest in Central and South America. The repeatability of long baselines (400-1000 km) in South America is improved by up to a factor of two in the horizontal vector baseline components by using tracking stations in the Pacific and Europe to supplement stations in North America. In every case but one, the differences between the mean solutions obtained using different tracking networks was equal to or smaller than day-to-day rms repeatabilities for the same baselines. The mean solutions obtained by using tracking stations in North America and the Pacific agreed at the 2-3 millimeter level with those using tracking stations in North America and Europe. The agreement of the extended tracking network solutions suggests that a broad distribution of tracking stations provides better geometric constraints on the satellite orbits and that solutions are not sensitive to changes in tracking network configuration when an extended network is use. A comparison of the results from the North Andes and a baseline in North America suggests that the use of a geometrically strong extended tracking network is most important when the network of interest is far from North America.

Freymueller, Jeffrey T.↗
NASA NTRS2012

Changes in Extratropical Storm Track Cloudiness 1983-2008: Observational Support for a Poleward Shift

Climate model simulations suggest that the extratropical storm tracks will shift poleward as a consequence of global warming. In this study the northern and southern hemisphere storm tracks over the Pacific and Atlantic ocean basins are studied using observational data, primarily from the International Satellite Cloud Climatology Project, ISCCP. Potential shifts in the storm tracks are examined using the observed cloud structures as proxies for cyclone activity. Different data analysis methods are employed, with the objective to address difficulties and uncertainties in using ISCCP data for regional trend analysis. In particular, three data filtering techniques are explored; excluding specific problematic regions from the analysis, regressing out a spurious viewing geometry effect, and excluding specific cloud types from the analysis. These adjustments all, to varying degree, moderate the cloud trends in the original data but leave the qualitative aspects of those trends largely unaffected. Therefore, our analysis suggests that ISCCP data can be used to interpret regional trends in cloudiness, provided that data and instrumental artefacts are recognized and accounted for. The variation in magnitude between trends emerging from application of different data correction methods, allows us to estimate possible ranges for the observational changes. It is found that the storm tracks, here represented by the extent of the midlatitude-centered band of maximum cloud cover over the studied ocean basins, experience a poleward shift as well as a narrowing over the 25 year period covered by ISCCP. The observed magnitudes of these effects are larger than in current generation climate models (CMIP3). The magnitude of the shift is particularly large in the northern hemisphere Atlantic. This is also the one of the four regions in which imperfect data primarily prevents us from drawing firm conclusions. The shifted path and reduced extent of the storm track cloudiness is accompanied by a regional reduction in total cloud cover. This decrease in cloudiness can primarily be ascribed to low level clouds, whereas the upper level cloud fraction actually increases, according to ISCCP. Independent satellite observations of radiative fluxes at the top of the atmosphere are consistent with the changes in total cloud cover. The shift in cloudiness is also supported by a shift in central position of the mid-troposphere meridional temperature gradient. We do not find support for aerosols playing a significant role in the satellite observed changes in cloudiness. The observed changes in storm track cloudiness can be related to local cloud-induced changes in radiative forcing, using ERBE and CERES radiative fluxes. The shortwave and the longwave components are found to act together, leading to a positive (warming) net radiative effect in response to the cloud changes in the storm track regions, indicative of positive cloud feedback. Among the CMIP3 models that simulate poleward shifts in all four storm track areas, all but one show decreasing cloud amount on a global mean scale in response to increased CO2 forcing, further consistent with positive cloud feedback. Models with low equilibrium climate sensitivity to a lesser extent than higher-sensitivity models simulate a poleward shift of the storm tracks.

Bender, Frida A-M.↗
NASA NTRS2015

Solar Flare Track Exposure Ages in Regolith Particles: A Calibration for Transmission Electron Microscope Measurements

Mineral grains in lunar and asteroidal regolith samples provide a unique record of their interaction with the space environment. Space weathering effects result from multiple processes including: exposure to the solar wind, which results in ion damage and implantation effects that are preserved in the rims of grains (typically the outermost 100 nm); cosmic ray and solar flare activity, which result in track formation; and impact processes that result in the accumulation of vapor-deposited elements, impact melts and adhering grains on particle surfaces. Determining the rate at which these effects accumulate in the grains during their space exposure is critical to studies of the surface evolution of airless bodies. Solar flare energetic particles (mainly Fe-group nuclei) have a penetration depth of a few millimeters and leave a trail of ionization damage in insulating materials that is readily observable by transmission electron microscope (TEM) imaging. The density of solar flare particle tracks is used to infer the length of time an object was at or near the regolith surface (i.e., its exposure age). Track measurements by TEM methods are routine, yet track production rate calibrations have only been determined using chemical etching techniques [e.g., 1, and references therein]. We used focused ion beam-scanning electron microscope (FIB-SEM) sample preparation techniques combined with TEM imaging to determine the track density/exposure age relations for lunar rock 64455. The 64455 sample was used earlier by [2] to determine a track production rate by chemical etching of tracks in anorthite. Here, we show that combined FIB/TEM techniques provide a more accurate determination of a track production rate and also allow us to extend the calibration to solar flare tracks in olivine.

Berger, Eve L.↗
NASA NTRS2019

A Kuiper Belt Source for Solar Flare Track-Rich Interplanetary Dust Particles

The presence of solar flare particle tracks in mineral grains within interplanetary dust particles (IDPs) has long been accepted as proof of their extraterrestrial origin [e.g. 1]. The 10-20 micrometers diameter IDPs released by dust producing objects in the solar system (mainly comets and asteroids) spiral in towards the Sun under the influence of Poynting-Robertson (PR) drag forces [2] and accumulate solar flare energetic particle tracks during their journey. The number of IDPs with well-constrained track density measurements is small, owing to the difficulty in the measurements and the lack of appropriatelysized crystals in which to image them. In order to use track densities as a chronometer of space exposure, the track production rate must be known. All previous work relied on track production rates determined by chemical etching techniques [e.g. 3], but tracks in IDPs are measured using TEM imaging. Here we report measurements of track densities in IDPs from both the anhydrous and hydrated IDP groups. Using the track production calibration determined from TEM observations of anorthite and olivine in lunar rock 64455 [4] we estimate space exposure times for these IDPs to constrain their parent body sources.

Keller, L. P.↗
NASA NTRS2020

Automatic Feature Tracking on Small Bodies for Autonomous Approach

Abstract—The autonomous approach of a spacecraft to an asteroid or comet (a small body) relies heavily on visual feature tracking to aid in estimating relative trajectories and the properties of the small body. Feature tracking for small bodies brings several challenges, including changing lighting, poor visual texture, and a concentration of features in a small part of an image. Six existing, open-source algorithms for feature tracking were tested on a simulated dataset and compared to the ground truth in the path of features. The main finding is that none of the algorithms provide all of the desired characteristics of long feature tracks with low errors and few outliers. Instead, there is a trade-off between long feature tracks and low error. The feature-matching algorithms SIFT, and BRISK provide good error characteristics, but short feature tracks, whereas the optical flow algorithm KLT provides long feature tracks, but with many features of large error. Given the challenges in feature tracking, it is recommended to focus development on each component of a feature tracking system: detection, description, and outlier rejection.

Morrell, Benjamin J↗
DOE OSTI2022

Reconstruction of Large Radius Tracks with the Exa.TrkX pipeline

Particle tracking is a challenging pattern recognition task at the Large Hadron Collider (LHC) and the High Luminosity-LHC. Conventional algorithms, such as those based on the Kalman Filter, achieve excellent performance in reconstructing the prompt tracks from the collision points. However, they require dedicated configuration and additional computing time to efficiently reconstruct the large radius tracks created away from the collision points. We developed an end-to-end machine learning-based track finding algorithm for the HL-LHC, the Exa.TrkX pipeline. The pipeline is designed so as to be agnostic about global track positions. In this work, we study the performance of the Exa.TrkX pipeline for finding large radius tracks. Trained with all tracks in the event, the pipeline simultaneously reconstructs prompt tracks and large radius tracks with high efficiencies. This new capability offered by the Exa.TrkX pipeline may enable us to search for new physics in real time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

A Confidence-Guided Technique for Tracking Time-Varying Features

Application scientists often employ feature tracking algorithms to capture the temporal evolution of various features in their simulation data. However, as the complexity of the scientific features is increasing with the advanced simulation modeling techniques, quantification of reliability of the feature tracking algorithms is becoming important. One of the desired requirements for any robust feature tracking algorithm is to estimate its confidence during each tracking step so that the results obtained can be interpreted without any ambiguity. To address this, we develop a confidence-guided feature tracking algorithm that allows reliable tracking of user-selected features and presents the tracking dynamics using a graph-based visualization along with the spatial visualization of the tracked feature. Here, the efficacy of the proposed method is demonstrated by applying it to two scientific datasets containing different types of time-varying features.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Data association algorithm for large-scale multi-object tracking with complex interactions

We present an online multi-object tracking algorithm to track multiple objects across a large number of image frames. Our work is motivated by the need to study evolution of nanoscale objects by transmission electron microscopy. The proposed approach is based on the existing multi-way data association tracking algorithm that is capable of tracking interacting objects with complex behaviors (i.e., merge, split, overlap, and appearance or disappearance). The multi-way data association is an offline algorithm to associate objects across all image frames at one step with a global optimization, which does not scale very well for large number of image frames. The proposed online tracking algorithm processes image frames as they arrive by detecting all objects in the newly arrived image frame and making the associations of the objects to those detected from the previous frame by the multi-way data association. This frameby-frame association scheme can cause fragmented traces of the objects that are occasionally misdetected for some image frames. We overcome this issue by allowing previously unassociated objects to be associated when the objects reappear within a fixed number of future image frames, namely the frame-delayed association. We combine the multi-way data association with the frame-delayed association to be able to track interacting objects with accurate handling of object disappearance events. The proposed method is validated through applications to simulated multi-object tracking problem and a real multi-object tracking problem. Here, the outcome of the proposed method is compared with four state-of-the-art algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Analysis of the track- and dose-averaged LET and LET spectra in proton therapy using the GEANT4 Monte Carlo code

Purpose: The motivation of this study was to find and eliminate the cause of errors in dose-averaged linear energy transfer (LET) calculations from therapeutic protons in small targets, such as biological cell layers, calculated using the GEANT 4 Monte Carlo code. Furthermore, the purpose was also to provide a recommendation to select an appropriate LET quantity from GEANT 4 simulations to correlate with biological effectiveness of therapeutic protons. Methods: The authors developed a particle tracking step based strategy to calculate the average LET quantities (track-averaged LET, LET{sub t} and dose-averaged LET, LET{sub d}) using GEANT 4 for different tracking step size limits. A step size limit refers to the maximally allowable tracking step length. The authors investigated how the tracking step size limit influenced the calculated LET{sub t} and LET{sub d} of protons with six different step limits ranging from 1 to 500 μm in a water phantom irradiated by a 79.7-MeV clinical proton beam. In addition, the authors analyzed the detailed stochastic energy deposition information including fluence spectra and dose spectra of the energy-deposition-per-step of protons. As a reference, the authors also calculated the averaged LET and analyzed the LET spectra combining the Monte Carlo method and the deterministic method. Relative biological effectiveness (RBE) calculations were performed to illustrate the impact of different LET calculation methods on the RBE-weighted dose. Results: Simulation results showed that the step limit effect was small for LET{sub t} but significant for LET{sub d}. This resulted from differences in the energy-deposition-per-step between the fluence spectra and dose spectra at different depths in the phantom. Using the Monte Carlo particle tracking method in GEANT 4 can result in incorrect LET{sub d} calculation results in the dose plateau region for small step limits. The erroneous LET{sub d} results can be attributed to the algorithm to determine fluctuations in energy deposition along the tracking step in GEANT 4. The incorrect LET{sub d} values lead to substantial differences in the calculated RBE. Conclusions: When the GEANT 4 particle tracking method is used to calculate the average LET values within targets with a small step limit, such as smaller than 500 μm, the authors recommend the use of LET{sub t} in the dose plateau region and LET{sub d} around the Bragg peak. For a large step limit, i.e., 500 μm, LET{sub d} is recommended along the whole Bragg curve. The transition point depends on beam parameters and can be found by determining the location where the gradient of the ratio of LET{sub d} and LET{sub t} becomes positive.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Technical aspects of real time positron emission tracking for gated radiotherapy

Purpose: Respiratory motion can lead to treatment errors in the delivery of radiotherapy treatments. Respiratory gating can assist in better conforming the beam delivery to the target volume. We present a study of the technical aspects of a real time positron emission tracking system for potential use in gated radiotherapy. Methods: The tracking system, called PeTrack, uses implanted positron emission markers and position sensitive gamma ray detectors to track breathing motion in real time. PeTrack uses an expectation–maximization algorithm to track the motion of fiducial markers. A normalized least mean squares adaptive filter predicts the location of the markers a short time ahead to account for system response latency. The precision and data collection efficiency of a prototype PeTrack system were measured under conditions simulating gated radiotherapy. The lung insert of a thorax phantom was translated in the inferior–superior direction with regular sinusoidal motion and simulated patient breathing motion (maximum amplitude of motion ±10 mm, period 4 s). The system tracked the motion of a {sup 22}Na fiducial marker (0.34 MBq) embedded in the lung insert every 0.2 s. The position of the was marker was predicted 0.2 s ahead. For sinusoidal motion, the equation used to model the motion was fitted to the data. The precision of the tracking was estimated as the standard deviation of the residuals. Software was also developed to communicate with a Linac and toggle beam delivery. In a separate experiment involving a Linac, 500 monitor units of radiation were delivered to the phantom with a 3 × 3 cm photon beam and with 6 and 10 MV accelerating potential. Radiochromic films were inserted in the phantom to measure spatial dose distribution. In this experiment, the period of motion was set to 60 s to account for beam turn-on latency. The beam was turned off when the marker moved outside of a 5-mm gating window. Results: The precision of the tracking in the IS direction was 0.53 mm for a sinusoidally moving target, with an average count rate ∼250 cps. The average prediction error was 1.1 ± 0.6 mm when the marker moved according to irregular patient breathing motion. Across all beam deliveries during the radiochromic film measurements, the average prediction error was 0.8 ± 0.5 mm. The maximum error was 2.5 mm and the 95th percentile error was 1.5 mm. Clear improvement of the dose distribution was observed between gated and nongated deliveries. The full-width at halfmaximum of the dose profiles of gated deliveries differed by 3 mm or less than the static reference dose distribution. Monitoring of the beam on/off times showed synchronization with the location of the marker within the latency of the system. Conclusions: PeTrack can track the motion of internal fiducial positron emission markers with submillimeter precision. The system can be used to gate the delivery of a Linac beam based on the position of a moving fiducial marker. This highlights the potential of the system for use in respiratory-gated radiotherapy.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

SU-G-BRA-06: Quantification of Tracking Performance of a Multi-Layer Electronic Portal Imaging Device

Purpose: The purpose of this study was to quantify the improvement in tumor tracking, with and without fiducial markers, afforded by employing a multi-layer (MLI) electronic portal imaging device (EPID) over the current state-of-the-art, single-layer, digital megavolt imager (DMI) architecture. Methods: An ideal observer signal-to-noise ratio (d’) approach was used to quantify the ability of an MLI EPID and a current, state-of-the-art DMI EPID to track lung tumors from the treatment beam’s-eye-view. Using each detector modulation transfer function (MTF) and noise power spectrum (NPS) as inputs, a detection task was employed with object functions describing simple three-dimensional Cartesian shapes (spheres and cylinders). Marker-less tumor tracking algorithms often use texture discrimination to differentiate benign and malignant tissue. The performance of such algorithms is simulated by employing a discrimination task for the ideal observer, which measures the ability of a system to differentiate two image quantities. These were defined as the measured textures for benign and malignant lung tissue. Results: The NNPS of the MLI ∼25% of that of the DMI at the expense of decreased MTF at intermediate frequencies (0.25≤<f≤<0.5 cycles/mm). Lung tissue textures were found to follow a power-law trend in the frequency domain. It was found that for both fiducial marker tracking and lung tissue texture differentiation, the signal power of each task is retained at low frequencies. Thus, little difference may be observed in signal power when DMI and MLI are compared. In all cases, improvements in tracking were greater than a factor of 2. Conclusion: MLI performance in tumor tracking is greatly improved by the additional imager layers. This implies that further improvements in tracking may be gained through increasing the thickness of each MLI layer. For tracking, the MLI performance is limited by noise response. Losses in MTF result in negligible differences in d’. The project was partially supported by a grant from Varian Medical Systems, Inc. and grant No. R01CA188446-01 from the National Cancer Institute.

60 APPLIED LIFE SCIENCES↗