Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tracking errors”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Particle tracking during Ostwald ripening using time-resolved laboratory X-ray microtomography

Laboratory X-ray microtomography is investigated as a method for obtaining time-resolved images of microstructural coarsening of the semisolid state of Al–5 wt.% Cu samples during Ostwald ripening. Owing to the 3D imaging capability of tomography, this technique uniquely provides access to the growth rates of individual particles, thereby not only allowing a statistical characterization of coarsening—as has long been possible by conventional metallography—but also enabling quantification of the influence of local environment on particle boundary migration. The latter information is crucial to understanding growth kinetics during Ostwald ripening at high volume fractions of the coarsening phase. Automated image processing and segmentation routines were developed to close gaps in the network of particle boundaries and to track individual particles from one annealing step to the next. The particle tracking success rate places an upper bound of only a few percent on the likelihood of segmentation errors for any given particle. The accuracy of particle size trajectories extracted from the time-resolved tomographic reconstructions is correspondingly high. Statistically averaged coarsening data and individual particle growth rates are in excellent agreement with the results of prior experimental studies and with computer simulations of Ostwald ripening. - Highlights: • Ostwald ripening in Al–5 wt.% Cu measured by laboratory X-ray microtomography • Time-resolved measurement of individual particle growth • Automated segmentation routines developed to close gaps in particle boundary network • Particle growth/shrinkage rates deviate from LSW model prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2016

An integrated model-driven method for in-treatment upper airway motion tracking using cine MRI in head and neck radiation therapy

Purpose: For the first time, MRI-guided radiation therapy systems can acquire cine images to dynamically monitor in-treatment internal organ motion. However, the complex head and neck (H&N) structures and low-contrast/resolution of on-board cine MRI images make automatic motion tracking a very challenging task. In this study, the authors proposed an integrated model-driven method to automatically track the in-treatment motion of the H&N upper airway, a complex and highly deformable region wherein internal motion often occurs in an either voluntary or involuntary manner, from cine MRI images for the analysis of H&N motion patterns. Methods: Considering the complex H&N structures and ensuring automatic and robust upper airway motion tracking, the authors firstly built a set of linked statistical shapes (including face, face-jaw, and face-jaw-palate) using principal component analysis from clinically approved contours delineated on a set of training data. The linked statistical shapes integrate explicit landmarks and implicit shape representation. Then, a hierarchical model-fitting algorithm was developed to align the linked shapes on the first image frame of a to-be-tracked cine sequence and to localize the upper airway region. Finally, a multifeature level set contour propagation scheme was performed to identify the upper airway shape change, frame-by-frame, on the entire image sequence. The multifeature fitting energy, including the information of intensity variations, edge saliency, curve geometry, and temporal shape continuity, was minimized to capture the details of moving airway boundaries. Sagittal cine MR image sequences acquired from three H&N cancer patients were utilized to demonstrate the performance of the proposed motion tracking method. Results: The tracking accuracy was validated by comparing the results to the average of two manual delineations in 50 randomly selected cine image frames from each patient. The resulting average dice similarity coefficient (93.28% ± 1.46%) and margin error (0.49 ± 0.12 mm) showed good agreement between the automatic and manual results. The comparison with three other deformable model-based segmentation methods illustrated the superior shape tracking performance of the proposed method. Large interpatient variations of swallowing frequency, swallowing duration, and upper airway cross-sectional area were observed from the testing cine image sequences. Conclusions: The proposed motion tracking method can provide accurate upper airway motion tracking results, and enable automatic and quantitative identification and analysis of in-treatment H&N upper airway motion. By integrating explicit and implicit linked-shape representations within a hierarchical model-fitting process, the proposed tracking method can process complex H&N structures and low-contrast/resolution cine MRI images. Future research will focus on the improvement of method reliability, patient motion pattern analysis for providing more information on patient-specific prediction of structure displacements, and motion effects on dosimetry for better H&N motion management in radiation therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Improving superficial target delineation in radiation therapy with endoscopic tracking and registration

Purpose: Target delineation within volumetric imaging is a critical step in the planning process of intensity modulated radiation therapy. In endoluminal cancers, endoscopy often reveals superficial areas of visible disease beyond what is seen on volumetric imaging. Quantitatively relating these findings to the volumetric imaging is prone to human error during the recall and contouring of the target. We have developed a method to improve target delineation in the radiation therapy planning process by quantitatively registering endoscopic findings contours traced on endoscopic images to volumetric imaging. Methods: Using electromagnetic sensors embedded in an endoscope, 2D endoscopic images were registered to computed tomography (CT) volumetric images by tracking the position and orientation of the endoscope relative to a CT image set. Regions-of-interest (ROI) in the 2D endoscopic view were delineated. A mesh created within the boundary of the ROI was projected onto the 3D image data, registering the ROI with the volumetric image. This 3D ROI was exported to clinical radiation treatment planning software. The precision and accuracy of the procedure was tested on two solid phantoms with superficial markings visible on both endoscopy and CT images. The first phantom was T-shaped tube with X-marks etched on the interior. The second phantom was an anatomically correct skull phantom with a phantom superficial lesion placed on the pharyngeal surface. Markings were contoured on the endoscope images and compared with contours delineated in the treatment planning system based on the CT images. Clinical feasibility was tested on three patients with early stage glottic cancer. Image-based rendering using manually identified landmarks was used to improve the registration. Results: Using the T-shaped phantom with X-markings, the 2D to 3D registration accuracy was 1.5-3.5 mm, depending on the endoscope position relative to the markings. Intraobserver standard variation was 0.5 mm. Rotational accuracy was within 2 deg. Using the skull phantom, registration accuracy was assessed by calculating the average surface minimum distance between the endoscopy and treatment planning contours. The average surface distance was 0.92 mm with 93% of all points in the 2D-endoscopy ROI within 1.5 mm of any point within the ROI contoured in the treatment planning software. This accuracy is limited by the CT imaging resolution and the electromagnetic (EM) sensor accuracy. The clinical testing demonstrated that endoscopic contouring is feasible. With registration based on em tracking only, accuracy was 5.6-8.4 mm. Image-based registration reduced this error to less than 3.5 mm and enabled endoscopic contouring in all cases. Conclusions: Registration of contours generated on 2D endoscopic images to 3D planning space is feasible, with accuracy smaller than typical set-up margins. Used in addition to standard 3D contouring methods in radiation planning, the technology may improve gross tumour volume (GTV) delineation for superficial tumors in luminal sites that are only visible in endoscopy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Three-dimensional liver motion tracking using real-time two-dimensional MRI

Purpose: Combined magnetic resonance imaging (MRI) systems and linear accelerators for radiotherapy (MR-Linacs) are currently under development. MRI is noninvasive and nonionizing and can produce images with high soft tissue contrast. However, new tracking methods are required to obtain fast real-time spatial target localization. This study develops and evaluates a method for tracking three-dimensional (3D) respiratory liver motion in two-dimensional (2D) real-time MRI image series with high temporal and spatial resolution. Methods: The proposed method for 3D tracking in 2D real-time MRI series has three steps: (1) Recording of a 3D MRI scan and selection of a blood vessel (or tumor) structure to be tracked in subsequent 2D MRI series. (2) Generation of a library of 2D image templates oriented parallel to the 2D MRI image series by reslicing and resampling the 3D MRI scan. (3) 3D tracking of the selected structure in each real-time 2D image by finding the template and template position that yield the highest normalized cross correlation coefficient with the image. Since the tracked structure has a known 3D position relative to each template, the selection and 2D localization of a specific template translates into quantification of both the through-plane and in-plane position of the structure. As a proof of principle, 3D tracking of liver blood vessel structures was performed in five healthy volunteers in two 5.4 Hz axial, sagittal, and coronal real-time 2D MRI series of 30 s duration. In each 2D MRI series, the 3D localization was carried out twice, using nonoverlapping template libraries, which resulted in a total of 12 estimated 3D trajectories per volunteer. Validation tests carried out to support the tracking algorithm included quantification of the breathing induced 3D liver motion and liver motion directionality for the volunteers, and comparison of 2D MRI estimated positions of a structure in a watermelon with the actual positions. Results: Axial, sagittal, and coronal 2D MRI series yielded 3D respiratory motion curves for all volunteers. The motion directionality and amplitude were very similar when measured directly as in-plane motion or estimated indirectly as through-plane motion. The mean peak-to-peak breathing amplitude was 1.6 mm (left-right), 11.0 mm (craniocaudal), and 2.5 mm (anterior-posterior). The position of the watermelon structure was estimated in 2D MRI images with a root-mean-square error of 0.52 mm (in-plane) and 0.87 mm (through-plane). Conclusions: A method for 3D tracking in 2D MRI series was developed and demonstrated for liver tracking in volunteers. The method would allow real-time 3D localization with integrated MR-Linac systems.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Technical Note: A novel leaf sequencing optimization algorithm which considers previous underdose and overdose events for MLC tracking radiotherapy

Purpose: Multileaf collimator (MLC) tracking radiotherapy is complex as the beam pattern needs to be modified due to the planned intensity modulation as well as the real-time target motion. The target motion cannot be planned; therefore, the modified beam pattern differs from the original plan and the MLC sequence needs to be recomputed online. Current MLC tracking algorithms use a greedy heuristic in that they optimize for a given time, but ignore past errors. To overcome this problem, the authors have developed and improved an algorithm that minimizes large underdose and overdose regions. Additionally, previous underdose and overdose events are taken into account to avoid regions with high quantity of dose events. Methods: The authors improved the existing MLC motion control algorithm by introducing a cumulative underdose/overdose map. This map represents the actual projection of the planned tumor shape and logs occurring dose events at each specific regions. These events have an impact on the dose cost calculation and reduce recurrence of dose events at each region. The authors studied the improvement of the new temporal optimization algorithm in terms of the L1-norm minimization of the sum of overdose and underdose compared to not accounting for previous dose events. For evaluation, the authors simulated the delivery of 5 conformal and 14 intensity-modulated radiotherapy (IMRT)-plans with 7 3D patient measured tumor motion traces. Results: Simulations with conformal shapes showed an improvement of L1-norm up to 8.5% after 100 MLC modification steps. Experiments showed comparable improvements with the same type of treatment plans. Conclusions: A novel leaf sequencing optimization algorithm which considers previous dose events for MLC tracking radiotherapy has been developed and investigated. Reductions in underdose/overdose are observed for conformal and IMRT delivery.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2012

A Measurement of the Two Track Charged Current Quasi-Elastic Cross Section with the MINER$\nu$A Detector

The MINER$\nu$A detector collected 1.2 × 10 20 POT of $\nu$ μ data with the NuMI beam at Fermi National Laboratory from 03/2010 to 07/2010. I will describe an analysis to identify two track charged current muon neutrino quasi-elastic interactions using this data set. The final state muon was identified by curvature using the MINOS detector. The final state contained proton was identified by the energy deposited per unit length at the end of the track. Two additional cuts were applied to separate signal from background: a cut on the extra energy in the event and a cut on a proton energy dependent variable. The dσ/dQ2 using this selection is presented. An initial list of systematic errors is introduced and calculated. In addition to this analysis work, I also developed both pattern recognition schemes used by the experiment, wrote MINOS muon track matching algorithms, implemented a Kalman Filter track fitting algorithm and contributed to various other tasks to the MINER$\nu$A experiment; I will provide an explanation of these. I will also give a brief overview of neutrino physics, the NuMI beam, and the MINER$\nu$A detector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Sci—Fri PM: Topics — 08: The Role and Benefits of Electromagnetic Needle-Tracking Technologies in Brachytherapy

In modern brachytherapy, application of large doses of ionizing radiation in a limited number of fractions is frequent. Furthermore, as with any surgical procedures, brachytherapy is subject to learning curve effects. In this context, there could be advantages of integrating real-time tracking of needles/catheters to existing protocols given the recent prominent advances in tracking technologies. In this work, we review the use of an electromagnetic tracking system (EMTS) based on the second generation Aurora® Planar Field Generator (Northern Digital Inc) and custom design needles (Philips Healthcare) for brachytherapy applications. The position and orientation information is obtained from 5 degrees of freedom sensors. Basic system performance characterization is performed in well-controlled conditions to establish accuracy and reproducibility as well as potential interference from standard brachytherapy equipment. The results show that sensor locations can be tracked to within 0.04mm (la) when located within 26cm of the generator. Orientation accuracy of the needle remained within ±1° in the same region, but rose quickly at larger distances. The errors on position and orientation strongly dependent the sensor position in the characterization volume (500×500×500mm{sup 3}). The presence of an ultrasound probe was shown to have negligible effects on tracking accuracy. The use of EMTS for automatic catheter/applicator reconstruction was also explored. Reconstruction time was less than 10 sec/channel and tips identification was within 0.69±0.29mm of the reference values. Finally, we demonstrate that hollow needle designs with special EM adaptation also allow for real-time seed drop position estimation. In phantom experiments showed that drop positions were on average within 1.6±0.9mm of the reference position measured from μCT. Altogether, EMTS offer promising benefits in a wide range of brachytherapy applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Gimballed Tracking Mount Pointing Angle Qualification

Tonopah Test Range (TTR), in support of its testing mission and modernization effort acquired a fleet of new gimballed tracking mounts (GTMs) manufactured by BAE Systems. The new GTMs can be operated remotely during flight tests and provide near real-time target tracking data. Furthermore, test vehicle Time-Space-Position-Information (TSPI) is evaluated using post-test synchronized imagery and pointing angle measurements acquired from each tracking mount. To comply with the Nuclear Enterprise Assurance Program (NEAP), all measurements devices must be certified. In keeping with the NEAP program, qualification of the new GTMs have been assessed to confirm that their pointing angle measurements produce acceptable TSPI results. This study only evaluated the four GTMs as a stand-alone solution and found that the GTMs meet their performance requirement of 0.006 degrees RMS error (or less) for post-processed pointing angles and produced TSPI solution with error volumes on the order of one meter or less. The new GTMs will be utilized in combination with existing optical tracking mounts, which will only improve the accuracy of the resulting TSPI data product. Details regarding the approach, analysis, summary results, and conclusions are presented.

42 ENGINEERING↗
DOE OSTI2014

A neural network clustering algorithm for the ATLAS silicon pixel detector

A novel technique to identify and split clusters created by multiple charged particles in the ATLAS pixel detector using a set of artificial neural networks is presented. Such merged clusters are a common feature of tracks originating from highly energetic objects, such as jets. Neural networks are trained using Monte Carlo samples produced with a detailed detector simulation. This technique replaces the former clustering approach based on a connected component analysis and charge interpolation. The performance of the neural network splitting technique is quantified using data from proton-proton collisions at the LHC collected by the ATLAS detector in 2011 and from Monte Carlo simulations. This technique reduces the number of clusters shared between tracks in highly energetic jets by up to a factor of three. It also provides more precise position and error estimates of the clusters in both the transverse and longitudinal impact parameter resolution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Multilevel Conditional Disturbance Rejection Control for Satellite Attitude Tracking

Recently, the conditional disturbance rejection controller (CDRC) was proposed to improve control performance by leveraging disturbances that have beneficial effects on the system. However, it only considers disturbances acting on the state variable directly influenced by them. Although fast convergence of this state can be achieved with the CDRC, it may unintentionally affect the convergence of the output (i.e., the primary state). Here, in this article, a multilevel CDRC is proposed to enhance satellite attitude control performance by accounting for the effect of disturbances on both attitude (output) and angular velocity. The extended state observer is employed to estimate the lumped disturbance, including the modeling errors and external disturbances. Then, a backstepping-based controller with the multilevel disturbance rejection law (ML-DRL) is designed for attitude tracking. The ML-DRL is developed to improve the control performance by using a disturbance with a damping effect on both attitude and velocity. Faster convergence of attitude and velocity can be achieved by conditionally compensating for the disturbance. The stability of the proposed control method is analyzed by demonstrating that the errors are bounded as time tends to infinity. The attitude control performance of the proposed method is evaluated through numerical examples conducted using the MATLAB/Simulink Multibody tool.

Active disturbance rejection control (ADRC)↗
DOE OSTI2013

Monte Carlo mesh tallies based on a Kernel Density Estimator approach using integrated particle tracks

A new Monte Carlo mesh tally based on a Kernel Density Estimator (KDE) approach using integrated particle tracks is presented. We first derive the KDE integral-track estimator and present a brief overview of its implementation as an alternative to the MCNP fmesh tally. To facilitate a valid quantitative comparison between these two tallies for verification purposes, there are two key issues that must be addressed. The first of these issues involves selecting a good data transfer method to convert the nodal-based KDE results into their cell-averaged equivalents (or vice versa with the cell-averaged MCNP results). The second involves choosing an appropriate resolution of the mesh, since if it is too coarse this can introduce significant errors into the reference MCNP solution. After discussing both of these issues in some detail, we present the results of a convergence analysis that shows the KDE integral-track and MCNP fmesh tallies are indeed capable of producing equivalent results for some simple 3D transport problems. In all cases considered, there was clear convergence from the KDE results to the reference MCNP results as the number of particle histories was increased. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A method for verification of treatment delivery in HDR prostate brachytherapy using a flat panel detector for both imaging and source tracking

Purpose: Verification of high dose rate (HDR) brachytherapy treatment delivery is an important step, but is generally difficult to achieve. A technique is required to monitor the treatment as it is delivered, allowing comparison with the treatment plan and error detection. In this work, we demonstrate a method for monitoring the treatment as it is delivered and directly comparing the delivered treatment with the treatment plan in the clinical workspace. This treatment verification system is based on a flat panel detector (FPD) used for both pre-treatment imaging and source tracking. Methods: A phantom study was conducted to establish the resolution and precision of the system. A pretreatment radiograph of a phantom containing brachytherapy catheters is acquired and registration between the measurement and treatment planning system (TPS) is performed using implanted fiducial markers. The measured catheter paths immediately prior to treatment were then compared with the plan. During treatment delivery, the position of the {sup 192}Ir source is determined at each dwell position by measuring the exit radiation with the FPD and directly compared to the planned source dwell positions. Results: The registration between the two corresponding sets of fiducial markers in the TPS and radiograph yielded a registration error (residual) of 1.0 mm. The measured catheter paths agreed with the planned catheter paths on average to within 0.5 mm. The source positions measured with the FPD matched the planned source positions for all dwells on average within 0.6 mm (s.d. 0.3, min. 0.1, max. 1.4 mm). Conclusions: We have demonstrated a method for directly comparing the treatment plan with the delivered treatment that can be easily implemented in the clinical workspace. Pretreatment imaging was performed, enabling visualization of the implant before treatment delivery and identification of possible catheter displacement. Treatment delivery verification was performed by measuring the source position as each dwell was delivered. This approach using a FPD for imaging and source tracking provides a noninvasive method of acquiring extensive information for verification in HDR prostate brachytherapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Reconstruction of Charged Particle Tracks in Realistic Detector Geometry Using a Vectorized and Parallelized Kalman Filter Algorithm

One of the most computationally challenging problems expected for the High-Luminosity Large Hadron Collider (HL-LHC) is finding and fitting particle tracks during event reconstruction. Algorithms used at the LHC today rely on Kalman filtering, which builds physical trajectories incrementally while incorporating material e ects and error estimation. Recognizing the need for faster computational throughput, we have adapted Kalman-filterbased methods for highly parallel, many-core SIMD and SIMT architectures that are now prevalent in high-performance hardware. Previously we observed significant parallel speedups, with physics performance comparable to CMS standard tracking, on Intel Xeon, Intel Xeon Phi, and (to a limited extent) NVIDIA GPUs. While early tests were based on artificial events occurring inside an idealized barrel detector, we showed subsequently that our mkFit software builds tracks successfully from complex simulated events (including detector pileup) occurring inside a geometrically accurate representation of the CMS-2017 tracker. Here, we report on advances in both the computational and physics performance of mkFit, as well as progress toward integration with CMS production software. Recently we have improved the overall eciency of the algorithm by preserving short track candidates at a relatively early stage rather than attempting to extend them over many layers. Moreover, mkFit formerly produced an excess of duplicate tracks; these are now explicitly removed in an additional processing step. We demonstrate that with these enhancements, mkFit becomes a suitable choice for the first iteration of CMS tracking, and eventually for later iterations as well. We plan to test this capability in the CMS High Level Trigger during Run 3 of the LHC, with an ultimate goal of using it in both the CMS HLT and oine reconstruction for the HL-LHC CMS tracker.

Cerati, Giuseppe↗
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

TH-AB-202-03: A Novel Tool for Computing Deliverable Doses in Dynamic MLC Tracking Treatments

Purpose: In tracked dynamic multi-leaf collimator (MLC) treatments, segments are continuously adapted to the target centroid motion in beams-eye-view. On-the-fly segment adaptation, however, potentially induces dosimetric errors due to the finite MLC leaf width and non-rigid target motion. In this study, we outline a novel tool for computing the 4d dose of lung SBRT plans delivered with MLC tracking. Methods: The following automated workflow was developed: A) centroid tracking, where the initial segments are morphed to each 4dCT phase based on the beams-eye-view GTV shift (followed by a dose calculation on each phase); B) re-optimized tracking, in which all morphed initial plans from (A) are further optimised (“warm-started”) in each 4dCT phase using the initial optimisation parameters but phase-specific volume definitions. Finally, both dose sets are accumulated to the reference phase using deformable image registration. Initial plans were generated according to the RTOG-1021 guideline (54Gy, 3-Fx, equidistant 9-beam IMRT) on the peak-exhale (reference) phase of a phase-binned 4dCT. Treatment planning and delivery simulations were performed in RayStation (research v4.6) using our in-house segment-morphing algorithm, which directly links to RayStation through a native C++ interface. Results: Computing the tracking plans and 4d dose distributions via the in-house interface takes 5 and 8 minutes respectively for centroid and re-optimized tracking. For a sample lung SBRT patient with 14mm peak-to-peak motion in sup-inf direction, mainly perpendicular leaf motion (0-collimator) resulted in small dose changes for PTV-D95 (−13cGy) and GTV-D98 (+18cGy) for the centroid tracking case compared to the initial plan. Modest reductions of OAR doses (e.g. spinal cord D2: −11cGy) were achieved in the idealized tracking case. Conclusion: This study presents an automated “1-click” workflow for computing deliverable MLC tracking doses in RayStation. Adding a non-deliverable re-optimized tracking scenario is expected to help quantify plan robustness for more challenging patients with anatomy deformations. We acknowledge support of the MLC tracking research from Elekta AB. MFF is supported by Cancer Research UK under Programme C33589/A19908. Research at ICR is also supported by Cancer Research UK under Programme C33589/A19727 and NHS funding to the NIHR Biomedical Research Centre at RMH and ICR.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-G-BRD-04: BEST IN PHYSICS (JOINT IMAGING-THERAPY): An Integrated Model-Based Intrafractional Organ Motion Tracking Approach with Dynamic MRI in Head and Neck Radiotherapy

Purpose: In-treatment dynamic cine images, provided by the first commercially available MRI-guided radiotherapy system, allow physicians to observe intrafractional motion of head and neck (H&N) internal structures. Nevertheless, high anatomical complexity and relatively poor cine image contrast/resolution have complicated automatic intrafractional motion evaluation. We proposed an integrated model-based approach to automatically delineate and analyze moving structures from on-board cine images. Methods: The H&N upper airway, a complex and highly deformable region wherein severe internal motion often occurs, was selected as the target-to-be-tracked. To reliably capture its motion, a hierarchical structure model containing three statistical shapes (face, face-jaw, and face-jaw-palate) was first built from a set of manually delineated shapes using principal component analysis. An integrated model-fitting algorithm was then employed to align the statistical shapes to the first to-be-detected cine frame, and multi-feature level-set contour propagation was performed to identify the airway shape change in the remaining frames. Ninety sagittal cine MR image sets, acquired from three H&N cancer patients, were utilized to demonstrate this approach. Results: The tracking accuracy was validated by comparing the results to the average of two manual delineations in 20 randomly selected images from each patient. The resulting dice similarity coefficient (93.28+/−1.46 %) and margin error (0.49+/−0.12 mm) showed good agreement with the manual results. Intrafractional displacements of anterior, posterior, inferior, and superior airway boundaries were observed, with values of 2.62+/−2.92, 1.78+/−1.43, 3.51+/−3.99, and 0.68+/−0.89 mm, respectively. The H&N airway motion was found to vary across directions, fractions, and patients, and highly correlated with patients’ respiratory frequency. Conclusion: We proposed the integrated computational approach, which for the first time allows to automatically identify the H&N upper airway and quantify in-treatment H&N internal motion in real-time. This approach can be applied to track other structures’ motion, and provide guidance on patient-specific prediction of intra-/inter-fractional structure displacements.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

The feasibility of head motion tracking in helical CT: A step toward motion correction

Purpose: To establish a practical and accurate motion tracking method for the development of rigid motion correction methods in helical x-ray computed tomography (CT). Methods: A commercially available optical motion tracking system provided 6 degrees of freedom pose measurements at 60 Hz. A 4 Multiplication-Sign 4 calibration matrix was determined to convert raw pose data acquired in tracker coordinates to a fixed CT coordinate system with origin at the isocenter of the scanner. Two calibration methods, absolute orientation (AO), and a new method based on image registration (IR), were compared by means of landmark analysis and correlation coefficient in phantom images coregistered using the derived motion transformations. Results: Transformations calculated using the IR-derived calibration matrix were found to be more accurate, with positional errors less than 0.5 mm (mean RMS), and highly correlated image voxel intensities. The AO-derived calibration matrix yielded larger mean RMS positional errors ( Asymptotically-Equal-To 1.0 mm), and poorer correlation coefficients. Conclusions: The authors have demonstrated the feasibility of accurate motion tracking for retrospective motion correction in helical CT. Their new IR-based calibration method based on image registration and function minimization was simpler to perform and delivered more accurate calibration matrices. This technique is a useful tool for future work on rigid motion correction in helical CT and potentially also other imaging modalities.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

FFA@CEBAF beam transport error and tolerance simulation studies

The Continuous Electron Beam Accelerator Facility (CEBAF) is a 12 GeV recirculating electron accelerator at the Thomas Jefferson National Accelerator Facility (JLAB). Major upgrades to the accelerator are being investigated which include a new 650 MeV injection beamline and state-of-the-art fixed-field alternating (FFA) gradient recirculation arcs. The upgrade will extend the energy of the electron beam to over 20 GeV. In this paper, we provide an error and tolerance simulation study of the amended beam optics transport of the existing accelerator tuned for 22 GeV operation. The study is conducted with the particle tracking codes elegant and Bmad in two parts. In the first part, we treat each section of the accelerator (electromagnetic arcs and linacs) modularly with ideal conditions at the beginning. The second part is a pseudo start-to-end (S2E) simulation with accumulated errors propagating from one beamline to the next.

Accelerator Physics↗