Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “motion”

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

Using an external surrogate for predictor model training in real-time motion management of lung tumors

Purpose: Precise prediction of respiratory motion is a prerequisite for real-time motion compensation techniques such as beam, dynamic couch, or dynamic multileaf collimator tracking. Collection of tumor motion data to train the prediction model is required for most algorithms. To avoid exposure of patients to additional dose from imaging during this procedure, the feasibility of training a linear respiratory motion prediction model with an external surrogate signal is investigated and its performance benchmarked against training the model with tumor positions directly. Methods: The authors implement a lung tumor motion prediction algorithm based on linear ridge regression that is suitable to overcome system latencies up to about 300 ms. Its performance is investigated on a data set of 91 patient breathing trajectories recorded from fiducial marker tracking during radiotherapy delivery to the lung of ten patients. The expected 3D geometric error is quantified as a function of predictor lookahead time, signal sampling frequency and history vector length. Additionally, adaptive model retraining is evaluated, i.e., repeatedly updating the prediction model after initial training. Training length for this is gradually increased with incoming (internal) data availability. To assess practical feasibility model calculation times as well as various minimum data lengths for retraining are evaluated. Relative performance of model training with external surrogate motion data versus tumor motion data is evaluated. However, an internal–external motion correlation model is not utilized, i.e., prediction is solely driven by internal motion in both cases. Results: Similar prediction performance was achieved for training the model with external surrogate data versus internal (tumor motion) data. Adaptive model retraining can substantially boost performance in the case of external surrogate training while it has little impact for training with internal motion data. A minimum adaptive retraining data length of 8 s and history vector length of 3 s achieve maximal performance. Sampling frequency appears to have little impact on performance confirming previously published work. By using the linear predictor, a relative geometric 3D error reduction of about 50% was achieved (using adaptive retraining, a history vector length of 3 s and with results averaged over all investigated lookahead times and signal sampling frequencies). The absolute mean error could be reduced from (2.0 ± 1.6) mm when using no prediction at all to (0.9 ± 0.8) mm and (1.0 ± 0.9) mm when using the predictor trained with internal tumor motion training data and external surrogate motion training data, respectively (for a typical lookahead time of 250 ms and sampling frequency of 15 Hz). Conclusions: A linear prediction model can reduce latency induced tracking errors by an average of about 50% in real-time image guided radiotherapy systems with system latencies of up to 300 ms. Training a linear model for lung tumor motion prediction with an external surrogate signal alone is feasible and results in similar performance as training with (internal) tumor motion. Particularly for scenarios where motion data are extracted from fluoroscopic imaging with ionizing radiation, this may alleviate the need for additional imaging dose during the collection of model training data.

60 APPLIED LIFE SCIENCES↗

SU-F-J-133: Adaptive Radiation Therapy with a Four-Dimensional Dose Calculation Algorithm That Optimizes Dose Distribution Considering Breathing Motion

Purpose: To model patient motion and produce four-dimensional (4D) optimized dose distributions that consider motion-artifacts in the dose calculation during the treatment planning process. Methods: An algorithm for dose calculation is developed where patient motion is considered in dose calculation at the stage of the treatment planning. First, optimal dose distributions are calculated for the stationary target volume where the dose distributions are optimized considering intensity-modulated radiation therapy (IMRT). Second, a convolution-kernel is produced from the best-fitting curve which matches the motion trajectory of the patient. Third, the motion kernel is deconvolved with the initial dose distribution optimized for the stationary target to produce a dose distribution that is optimized in four-dimensions. This algorithm is tested with measured doses using a mobile phantom that moves with controlled motion patterns. Results: A motion-optimized dose distribution is obtained from the initial dose distribution of the stationary target by deconvolution with the motion-kernel of the mobile target. This motion-optimized dose distribution is equivalent to that optimized for the stationary target using IMRT. The motion-optimized and measured dose distributions are tested with the gamma index with a passing rate of >95% considering 3% dose-difference and 3mm distance-to-agreement. If the dose delivery per beam takes place over several respiratory cycles, then the spread-out of the dose distributions is only dependent on the motion amplitude and not affected by motion frequency and phase. This algorithm is limited to motion amplitudes that are smaller than the length of the target along the direction of motion. Conclusion: An algorithm is developed to optimize dose in 4D. Besides IMRT that provides optimal dose coverage for a stationary target, it extends dose optimization to 4D considering target motion. This algorithm provides alternative to motion management techniques such as beam-gating or breath-holding and has potential applications in adaptive radiation therapy.

60 APPLIED LIFE SCIENCES↗

SU-G-JeP4-08: External and Internal Motion Assessments with Amplitude and Time Based Reconstruction and AveCT

Purpose: To investigate the effects of external surrogate and tumour motion by observing the reconstructed phases and AveCT in an Amplitude and Time based 4DCT. Methods: Based on patient motion studies, Cos6 and sinusoidal motions were simulated as external surrogate and tumour motions in a motion phantom. The diaphragm and tumour motions may or may not display the same waveform therefore the same and different waveforms were programmed into the phantom, scanned and reconstructed based on Amplitude and Time. The AveCT and phases were investigated with these different scenarios. The AveCT phantom images were also compared with CBCT phantom images programmed with the same motions. Results: For the same surrogate and tumour sin motions, the phases (Amplitude and Time) and AveCT indicated similar motions based on the position of the BB at the slice and displayed contrast values respectively. For cos6 motions, due to the varied time the tumour spends at each position, the Amplitude and Time based phases differed. The AveCT images represented the actual tumour motions and the Time and Amplitude based phases were represented by the surrogate with varied times. Conclusion: Different external surrogate and tumour motions may result in different displayed image motions when observing the AveCT and reconstructed phases. During the 4DCT, the surrogate motion is readily available for observation of the amplitude and time of the diaphragm position. Following image reconstruction, the user may need to observe the AveCT in addition to the reconstructed phases to comprehend the time weightings of the tumour motion during the scan. This may also apply to 3D CBCT images where the displayed tumour position in the images is influenced by the long duration of the CBCT. Knowledge of the tumour motion represented by the greyscale of the AveCT may also assist in CBCT treatment beam verification matching.

60 APPLIED LIFE SCIENCES↗

The use of vestibular models for design and evaluation of flight simulator motion

Quantitative models for the dynamics of the human vestibular system are applied to the design and evaluation of flight simulator platform motion. An optimal simulator motion control algorithm is generated to minimize the vector difference between perceived spatial orientation estimated in flight and in simulation. The motion controller has been implemented on the Vertical Motion Simulator at NASA Ames Research Center and evaluated experimentally through measurement of pilot performance and subjective rating during VTOL aircraft simulation. In general, pilot performance in a longitudinal tracking task (formation flight) did not appear to be sensitive to variations in platform motion condition as long as motion was present. However, pilot assessment of motion fidelity by means of a rating scale designed for this purpose, were sensitive to motion controller design. Platform motion generated with the optimal motion controller was found to be generally equivalent to that generated by conventional linear crossfeed washout. The vestibular models are used to evaluate the motion fidelity of transport category aircraft (Boeing 727) simulation in a pilot performance and simulator acceptability study at the Man-Vehicle Systems Research Facility at NASA Ames Research Center. Eighteen airline pilots, currently flying B-727, were given a series of flight scenarios in the simulator under various conditions of simulator motion. The scenarios were chosen to reflect the flight maneuvers that these pilots might expect to be given during a routine pilot proficiency check. Pilot performance and subjective rating of simulator fidelity was relatively insensitive to the motion condition, despite large differences in the amplitude of motion provided. This lack of sensitivity may be explained by means of the vestibular models, which predict little difference in the modeled motion sensations of the pilots when different motion conditions are imposed.

Bussolari, Steven R.↗

Transfer of Training on the Vertical Motion Simulator

This paper describes a transfer-of-training study in the NASA Ames Vertical Motion Simulator (VMS). Sixty-one general aviation pilots, divided in four groups, trained on four challenging commercial transport tasks with four motion conditions: no motion, small hexapod, large hexapod, and VMS motion. Then, every pilot repeated the tasks with VMS motion to determine if training with different motion conditions had an effect. New objective motion criteria guided the selection of the motion parameters for the small and large hexapod conditions. Considering results that were statistically significant, or marginally, the motion condition used in training affected 1) longitudinal and lateral touchdown location; 2) the number of secondary stall warnings in a stall recovery; 3) pilot ratings of motion utility and maximum load factor obtained in an overbanked upset recovery; and 4) pilot ratings of motion utility and pedal input reaction time in the engine-out-on-takeoff task. Since the motion condition revealed statistical differences on objective measures in all the tasks, even with some in the direction not predicted, trainers should be cautious not to oversimplify the effects of platform motion. Evidence suggests that the new objective motion criteria may offer valid standardization benefits, as instances arose when the higher-fidelity hexapod motion, as predicted by the criteria, provided better cues in training than the lower-fidelity hexapod motion.

Zaal, Peter M. T.↗

Effects of Motion Cues on the Training of Multi-Axis Manual Control Skills

The study described in this paper investigated the effects of two different hexapod motion configurations on the training and transfer of training of a simultaneous roll and pitch control task. Pilots were divided between two groups which trained either under a baseline hexapod motion condition, with motion typically provided by current training simulators, or an optimized hexapod motion condition, with increased fidelity of the motion cues most relevant for the task. All pilots transferred to the same full-motion condition, representing motion experienced in flight. A cybernetic approach was used that gave insights into the development of pilots use of visual and motion cues over the course of training and after transfer. Based on the current results, neither of the hexapod motion conditions can unambiguously be chosen as providing the best motion for training and transfer of training of the used multi-axis control task. However, the optimized hexapod motion condition did allow pilots to generate less visual lead, control with higher gains, and have better disturbance-rejection performance at the end of the training session compared to the baseline hexapod motion condition. Significant adaptations in control behavior still occurred in the transfer phase under the full-motion condition for both groups. Pilots behaved less linearly compared to previous single-axis control-task experiments; however, this did not result in smaller motion or learning effects. Motion and learning effects were more pronounced in pitch compared to roll. Finally, valuable lessons were learned that allow us to improve the adopted approach for future transfer-of-training studies.

flight simulators↗

Novel Assessment of Renal Motion in Children as Measured via Four-Dimensional Computed Tomography

Objectives: Abdominal intensity-modulated radiation therapy and proton therapy require quantification of target and organ motion to optimize localization and treatment. Although addressed in adults, there is no available literature on this issue in pediatric patients. We assessed physiologic renal motion in pediatric patients. Methods and Materials: Twenty free-breathing pediatric patients at a median age of 8 years (range, 2-18 years) with intra-abdominal tumors underwent computed tomography simulation and four-dimensional computed tomography acquisition (slice thickness, 3 mm). Kidneys and diaphragms were contoured during eight phases of respiration to estimate center-of-mass motion. We quantified center of kidney mass mobility vectors in three dimensions: anteroposterior (AP), mediolateral (ML), and superoinferior (SI). Results: Kidney motion decreases linearly with decreasing age and height. The 95% confidence interval for the averaged minima and maxima of renal motion in children younger than 9 years was 5-9 mm in the ML direction, 4-11 mm in the AP direction, and 12-25 mm in the SI dimension for both kidneys. In children older than 9 years, the same confidence interval reveals a widening range of motion that was 5-16 mm in the ML direction, 6-17 mm in the AP direction, and 21-52 mm in the SI direction. Although not statistically significant, renal motion correlated with diaphragm motion in older patients. The correlation between diaphragm motion and body mass index was borderline (r = 0.52, p = 0.0816) in younger patients. Conclusions: Renal motion is age and height dependent. Measuring diaphragmatic motion alone does not reliably quantify pediatric renal motion. Renal motion in young children ranges from 5 to 25 mm in orientation-specific directions. The vectors of motion range from 5 to 52 mm in older children. These preliminary data represent novel analyses of pediatric intra-abdominal organ motion.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-C-210-02: Impact of Intrafractional Motion On TomoTherapy Stereotactic Body Radiotherapy (SBRT) 4D Dosimetry

Purpose: TomoTherapy treatment has unique challenges in handling intrafractional motion compared to conventional LINAC. This study is aimed to gain a realistic and quantitative understanding of motion impact on TomoTherapy SBRT treatment of lung and prostate cancer patients. Methods: A 4D dose engine utilizing GPUs and including motion during treatment was developed for the efficient simulation of TomoTherapy delivered dosimetry. Two clinical CyberKnife lung cases with respiratory motion tracking and two prostate cases with a slower non-periodical organ motion treated by LINAC plus Calypso tracking were used in the study. For each disease site, one selected case has an average motion (6mm); the other has a large motion (10mm for lung and 15mm for prostate). SBRT of lung and prostate cases were re-planned on TomoTherapy with 12 Gyx4 fractions and 7Gyx5 fractions, respectively, all with 95% PTV coverage. Each case was planned with 4 jaw settings: 1) conventional 1cm static, 2) 2.5cm static, 3) 2.5cm dynamic, and 4) 5cm dynamic. The intrafractional rigid motion of the target was applied in the dose calculation of individual fractions of each plan and total dose was accumulated from multiple fractions. Results: For 1cm static jaw plans with motions applied, PTV coverage is related to motion type and amplitude. For SBRT patients with average motion (6mm), the PTV coverage remains > 95% for lung case and 74% for prostate case. For cases with large motion, PTV coverage drops to 61% for lung SBRT and 49% for prostate SBRT. Plans with other jaws improve uniformity of moving target, but still suffer from poor PTV coverage (< 70%). Conclusion: TomoTherapy lung SBRT is less motion-impacted when average amplitude of respiratory-induced intrafractional motion is present (6mm). When motion is large and/or non-periodic (prostate), all studied plans lead to significantly decreased target coverage in actual delivered dosimetry.

60 APPLIED LIFE SCIENCES↗

Validation of earthquake ground-motion models in southern California, USA, using precariously balanced rocks

Accurate estimates of earthquake ground shaking rely on uncertain ground-motion models derived from limited instrumental recordings of historical earthquakes. A critical issue is that there is currently no method to empirically validate the resultant ground-motion estimates of these models at the timescale of rare, large earthquakes; this lack of validation causes great uncertainty in ground-motion estimates. Here, we address this issue and validate ground-motion estimates for southern California utilizing the unexceeded ground motions recorded by 20 precariously balanced rocks. We used cosmogenic 10Be exposure dating to model the age of the precariously balanced rocks, which ranged from ca. 1 ka to ca. 50 ka, and calculated their probability of toppling at different ground-motion levels. With this rock data, we then validated the earthquake ground motions estimated by the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3) seismic-source characterization and the Next Generation Attenuation (NGA)-West2 ground-motion models. We found that no ground-motion model estimated levels of earthquake ground shaking consistent with the observed continued existence of all 20 precariously balanced rocks. The ground-motion model I14 estimated ground-motion levels that were inconsistent with the most rocks; therefore, I14 was invalidated and removed. At a 2475 year mean return period, the removal of this invalid ground-motion model resulted in a 2−7% reduction in the mean and a 10−36% reduction in the 5th−95th fractile uncertainty of the ground-motion estimates. Our findings demonstrate the value of empirical data from precariously balanced rocks as a validation tool for removing invalid ground-motion models and, in turn, reducing the uncertainty in earthquake ground-motion estimates.

Geology↗

Changes in plasma vasopressin during motion sickness in cats

Changes in levels of plasma vasopressin (AVP) and cortisol (C) have been shown to be correlated with motion sickness and nausea in man. As part of the research aimed at validation of the cat as an appropriate animal model for motion sickness research, levels of these hormones were investigated in the cat during motion sickness elicited by vertical linear acceleration of approximately 0.6 Hz and 1 +/- 0.6 G. In Study 1, 15 cats previously screened for susceptibility to motion sickness were prepared with indwelling jugular catheters to permit withdrawl of blood with minimal disruption of the stimulus and minimum stress to the animal. AVP and C were measured in blood samples obtained during exposure to vertical linear acceleration and during control sessions in which the animals were placed in the stationary apparatus. 10 min and 1 min prior to duration; 1, 5, 10, and 20 min after start of motion. Total duration of exposure to motion was 20 min. The data indicate that both AVP and C are elevated during exposure to motion if emesis occurs. AVP reaches maximum levels during or about the same time as emesis, while C increases gradually throughout the period of vertical acceleration. In Study 2, four cats were prepared with indwelling catheters and AVP was measured in blood withdrawn during exposure to the vertical linear acceleration. A single pre-motion sample consisting of three samples drawn 5 min prior to motion onset. Two series of samples consisting of three samples drawn at 3-min intervals were obtained during motion. The first series was initiated at emesis, and the second 25 min after emesis. Results show that levels of circulating AVP were elevated (2 to 27 times the control and pre-motion levels) in the samples taken during emesis and decreased, but remained 1 to 6 times above the pre-motion or control levels within 25 min. The results of these two studies indicate that AVP is elevated during motion-produced emesis than is C. These findings are in general agreement with those obtained from humans under motion sickness conditions, and indicate that it is appropriate to continue to use the cat in studies of hormone changes during motion sickness.

R Fox↗

Feasibility of Earthbound Motion in Lunar Gravity

BACKGROUND: Marginal stability of astronaut movement while performing lunar surface tasks has been well documented, and is clearly demonstrated in videos of falls, and near falls, during Apollo Lunar Extravehicular Activities. Referencing mission reports from Apollo 15 and 16 [1, 2], suspected causes for falls include: surface conditions, visibility, and gravitational effects (hypogravity). In this preliminary test, we employ the open-source biomechanical tool OpenSim [3, 4] to analyze the impact of lunar gravity (Lg) on two object-pickup motions performed by a single shirt-sleeved subject. Specifically, we attempt to answer the following questions based on an estimation of the Center of Mass Projection (CoMP) and ground reaction force Center of Pressure (COP) as it relates to the astronaut Base of Support (BOS) for 1g and Lg conditions: 1. Is the task motion, as performed in 1g, dynamically feasible in 1g and Lg? 2. Can we make the motion dynamically feasible in Lg by slowing it down? 3. Is the Lg COP equal to the 1g COP at a theoretically predicted reduction in motion speed? METHODS AND RESULTS: To answer the first question, the gravitational acceleration in the OpenSim model is modified from a nominal 9.81 m/s2 to 1.64 m/s2, and the 1g joint trajectory is input to an OpenSim based method [5] for estimating ground reaction forces and moments. From this method, the position of the CoMP and COP can be estimated and checked to see if they remain within a simulated BOS formed from the footprint of the OpenSim model to determine whether the motion is dynamically feasible. As expected, both of the motions were estimated to be feasible in a 1g environment, however, both motions had periods of infeasibility in Lg. It is well known that crew members make adjustments to motion trajectories in altered gravity fields to maintain balance. As a first step, we considered the simple adjustment of slowing the motion in Lg by a constant factor. This was accomplished by scaling the time stamps in the motion trajectory file by that factor. For the two motions considered, it was found that scale factors of 1.3 and 1.4 kept the COP just within the subject BOS. The CoMP is unchanged by the gravity field. Simple analysis of an inverted pendulum in the Lg environment, which generalizes to a general multibody system, leads to a theoretical prediction that a reduction in speed factor of √1g/Lg, or 2.445, will make the COP trajectory in Lg equivalent to that in 1g. When the above procedure was performed with a factor of 2.445, the estimated COP in Lg, was observed to be very close to that in 1g. In summary, we have developed a method for estimating the CoMP and COP in Lg, for subject motion collected in 1g. We believe this method can prove to be a valuable check and balance for simulated Lg training and testing by exposing potential simulator-induced artifacts that make the simulated task motion seem possible, when in fact, it would violate the above criteria. We also note that a reduction in task speed should tend the task motion towards stability, with a theoretical slowdown factor of √1g/Lg making the motion stability equal to that in 1g according to the CoMP and COP criteria.

R K Huffman↗

Motion coherence affects human perception and pursuit similarly

Pursuit and perception both require accurate information about the motion of objects. Recovering the motion of objects by integrating the motion of their components is a difficult visual task. Successful integration produces coherent global object motion, while a failure to integrate leaves the incoherent local motions of the components unlinked. We compared the ability of perception and pursuit to perform motion integration by measuring direction judgments and the concomitant eye-movement responses to line-figure parallelograms moving behind stationary rectangular apertures. The apertures were constructed such that only the line segments corresponding to the parallelogram's sides were visible; thus, recovering global motion required the integration of the local segment motion. We investigated several potential motion-integration rules by using stimuli with different object, vector-average, and line-segment terminator-motion directions. We used an oculometric decision rule to directly compare direction discrimination for pursuit and perception. For visible apertures, the percept was a coherent object, and both the pursuit and perceptual performance were close to the object-motion prediction. For invisible apertures, the percept was incoherently moving segments, and both the pursuit and perceptual performance were close to the terminator-motion prediction. Furthermore, both psychometric and oculometric direction thresholds were much higher for invisible apertures than for visible apertures. We constructed a model in which both perception and pursuit are driven by a shared motion-processing stage, with perception having an additional input from an independent static-processing stage. Model simulations were consistent with our perceptual and oculomotor data. Based on these results, we propose the use of pursuit as an objective and continuous measure of perceptual coherence. Our results support the view that pursuit and perception share a common motion-integration stage, perhaps within areas MT or MST.

NASA Center ARC↗

Motion Sickness and Concerns for Urban Air Mobility Vehicles: A Literature Review

Motion sickness is a general term for a constellation of signs and symptoms, generally due to exposure to abrupt, periodic, or unnatural accelerations, especially when traveling in a vehicle. Motion sickness results from a mismatch of the visual and nonvisual (vestibular and kinesthetic) information, the observed scene and the motion felt or lack of it. Motion sickness onset is associated with a pattern of physiological changes in heart rate, peripheral blood flow, respiration, and skin conductance and the pattern is repeatable for a particular subject but variable between subjects. Demographic factors such as gender and age that affect motion sickness are well known with children, women, and older adults more likely to be susceptible. Often motion sickness is assessed and quantified using variations of the motion sickness susceptibility questionnaires including the Pensacola Diagnostic Rating Scale and the Simulator Sickness Questionnaire. Even though symptoms are easily identified by such questionnaires, they commonly are subjective. Tools such as these questionnaires for screening individuals susceptible to motion sickness are useful, however, they are only mildly predictive. Moreover, models for predicting motion sickness, which have largely been developed for sea sickness, do not consider task characteristics. Predictions of motion sickness rates and prevalence for Urban Air Mobility (UAM) vehicles are not possible at present because data from actual flight or full-fidelity simulation are simply not yet available. Extrapolation from other modes of transportation (i.e., automobiles, buses, trains, boats, other types of aircraft) is difficult because of differences in the motion stimulus experienced, trip duration, and other factors. How UAM vehicles will change the social dynamics of passenger interaction and vehicle interior design changes (e.g., seat orientations) is unknown. Should motion sickness prove to be an issue, vehicle design modifications such as having passengers face forward, providing additional seat recline, giving each person their own climate control for airflow, perhaps ensuring the horizon is visible to all passengers (reducing visual occlusion by the headrest) and visually stabilizing displays on carry-on devices (smart phones, tablets, etc.) may benefit passengers. Several commercial companies provide wearable devices for physiological monitoring that have been validated and are suitable for use with passengers in UAM vehicles or high-fidelity simulators. Potential countermeasures for motion sickness include user-worn devices, anti-motion sickness medications, and non-pharmacological approaches such as biofeedback and Autogenic Feedback Training Exercise. Both simulator and in-vehicle UAM research is needed to evaluate the effectiveness of any potential countermeasure.

motion sickness↗

Effect of Soil Properties and Input Motion on Site Amplification Using Validated Nonlinear Soil Model

Here, three-dimensional (3-D) nonlinear site response analyses are conducted using finite element models of actual soil profiles from ten nuclear power plant (NPP) sites in the United States to investigate the effects of soil properties and input motions on site amplification. The modeling approach developed in this study combines several novel elements, such as 3-D analysis (including vertical motions), nonlinear inelastic behavior of soil (strain-dependent shear modulus reduction and hysteretic damping), formulation of nonreflecting boundary conditions at the base, and generation of realistic outcrop ground motions for specific sites. All these elements of the modeling approach are first validated using actual data from five earthquakes at three downhole array stations recorded in the Kiban-Kyoshin network (KiK-net), Japan. The same approach is then used to develop site models of ten NPP sites in the United States and corresponding ground motions that are spectrally matched to the site hazard spectra. Eight sets of three-component input motions are used in the study and are categorized on the basis of presence or absence of a near-field pulse in the seed ground motions used for spectral matching. It is found that all sites retain a definite site amplification function regardless of the input motion, provided that the seed motion is spectrally matched to the site hazard spectra. The magnitude of site amplification and frequencies at which they occur depend upon soil properties, particularly the shear wave velocity profile and the constitutive relationship (strain-dependent shear modulus reduction and hysteretic damping) of soil. Amplification of spectral acceleration in the vertical direction (up-down motion) is found to be just as much as, if not more than, the amplification in the horizontal direction. Peak shear strain is found to be about 20% larger for near-field motions compared to far-field motions whereas maximum horizontal site amplification for far-field motions is found to be consistently larger than that of near-field motions, even though the differences between the two remain within the scatter resulting from individual ground motions.

99 GENERAL AND MISCELLANEOUS↗

Evaluation of interpolation methods for surface-based motion compensated tomographic reconstruction for cardiac angiographic C-arm data

Purpose: For interventional cardiac procedures, anatomical and functional information about the cardiac chambers is of major interest. With the technology of angiographic C-arm systems it is possible to reconstruct intraprocedural three-dimensional (3D) images from 2D rotational angiographic projection data (C-arm CT). However, 3D reconstruction of a dynamic object is a fundamental problem in C-arm CT reconstruction. The 2D projections are acquired over a scan time of several seconds, thus the projection data show different states of the heart. A standard FDK reconstruction algorithm would use all acquired data for a filtered backprojection and result in a motion-blurred image. In this approach, a motion compensated reconstruction algorithm requiring knowledge of the 3D heart motion is used. The motion is estimated from a previously presented 3D dynamic surface model. This dynamic surface model results in a sparse motion vector field (MVF) defined at control points. In order to perform a motion compensated reconstruction, a dense motion vector field is required. The dense MVF is generated by interpolation of the sparse MVF. Therefore, the influence of different motion interpolation methods on the reconstructed image quality is evaluated. Methods: Four different interpolation methods, thin-plate splines (TPS), Shepard's method, a smoothed weighting function, and a simple averaging, were evaluated. The reconstruction quality was measured on phantom data, a porcine model as well as on in vivo clinical data sets. As a quality index, the 2D overlap of the forward projected motion compensated reconstructed ventricle and the segmented 2D ventricle blood pool was quantitatively measured with the Dice similarity coefficient and the mean deviation between extracted ventricle contours. For the phantom data set, the normalized root mean square error (nRMSE) and the universal quality index (UQI) were also evaluated in 3D image space. Results: The quantitative evaluation of all experiments showed that TPS interpolation provided the best results. The quantitative results in the phantom experiments showed comparable nRMSE of Almost-Equal-To 0.047 {+-} 0.004 for the TPS and Shepard's method. Only slightly inferior results for the smoothed weighting function and the linear approach were achieved. The UQI resulted in a value of Almost-Equal-To 99% for all four interpolation methods. On clinical human data sets, the best results were clearly obtained with the TPS interpolation. The mean contour deviation between the TPS reconstruction and the standard FDK reconstruction improved in the three human cases by 1.52, 1.34, and 1.55 mm. The Dice coefficient showed less sensitivity with respect to variations in the ventricle boundary. Conclusions: In this work, the influence of different motion interpolation methods on left ventricle motion compensated tomographic reconstructions was investigated. The best quantitative reconstruction results of a phantom, a porcine, and human clinical data sets were achieved with the TPS approach. In general, the framework of motion estimation using a surface model and motion interpolation to a dense MVF provides the ability for tomographic reconstruction using a motion compensation technique.

60 APPLIED LIFE SCIENCES↗

TU-F-17A-03: An Analytical Respiratory Perturbation Model for Lung Motion Prediction

Purpose: Breathing irregularity is common, causing unreliable prediction in tumor motion for correlation-based surrogates. Both tidal volume (TV) and breathing pattern (BP=ΔVthorax/TV, where TV=ΔVthorax+ΔVabdomen) affect lung motion in anterior-posterior and superior-inferior directions. We developed a novel respiratory motion perturbation (RMP) model in analytical form to account for changes in TV and BP in motion prediction from simulation to treatment. Methods: The RMP model is an analytical function of patient-specific anatomic and physiologic parameters. It contains a base-motion trajectory d(x,y,z) derived from a 4-dimensional computed tomography (4DCT) at simulation and a perturbation term Δd(ΔTV,ΔBP) accounting for deviation at treatment from simulation. The perturbation is dependent on tumor-specific location and patient-specific anatomy. Eleven patients with simulation and treatment 4DCT images were used to assess the RMP method in motion prediction from 4DCT1 to 4DCT2, and vice versa. For each patient, ten motion trajectories of corresponding points in the lower lobes were measured in both 4DCTs: one served as the base-motion trajectory and the other as the ground truth for comparison. In total, 220 motion trajectory predictions were assessed. The motion discrepancy between two 4DCTs for each patient served as a control. An established 5D motion model was used for comparison. Results: The average absolute error of RMP model prediction in superior-inferior direction is 1.6±1.8 mm, similar to 1.7±1.6 mm from the 5D model (p=0.98). Some uncertainty is associated with limited spatial resolution (2.5mm slice thickness) and temporal resolution (10-phases). Non-corrected motion discrepancy between two 4DCTs is 2.6±2.7mm, with the maximum of ±20mm, and correction is necessary (p=0.01). Conclusion: The analytical motion model predicts lung motion with accuracy similar to the 5D model. The analytical model is based on physical relationships, requires no training, and therefore is potentially more resilient to breathing irregularities. On-going investigation introduces airflow into the RMP model for improvement. This research is in part supported by NIH (U54CA137788/132378). AY would like to thank MSKCC summer medical student research program supported by National Cancer Institute and hosted by Department of Medical Physics at MSKCC.

60 APPLIED LIFE SCIENCES↗

SU-F-I-15: Evaluation of a New MR-Compatible Respiratory Motion Device at 3T

Purpose: Recent advances in MRI-guided radiotherapy has inspired the development of MRI-compatible motion devices that simulate patient periodic motion in the scanner, particularly respiratory motion. Most commercial devices rely on non MR-safe ferromagnetic stepper motors which are not practical for regular QA testing. This work evaluates the motion performance of a new fully MRI compatible respiratory motion device at 3T. Methods: The QUASAR™ MRI-compatible respiratory motion phantom has been recently developed by Modus QA Inc., London, ON, Canada. The prototype is constructed from diamagnetic materials with linear motion generated using MRI-compatible piezoelectric motors that can be safely inserted in the scanner bore. The tumor was represented by a fillable sphere and is attached to the linear motion generator. The spherical tumor-representative and its surroundings were filled with different concentrations of MnCl2 to produce realistic relaxation times. The motion was generated along the longitudinal (H/F) axis of the bore using sinusoidal reference waveform (amplitude = 15 mm, frequency 0.25 Hz). Imaging was then performed on 3T Philips Achieva using a 32-channel cardiac coil. Fast 2D spoiled gradient-echo was used with a spatial resolution of 1.8 × 1.8 mm{sup 2} and slice thickness of 4 mm. The motion waveform was then measured on the resultant image series by tracking the centroid of the sphere through the time series. This image-derived measured motion was compared to the software-generated reference waveform. Results: No visible distortions from the device were observed on the images. Excellent agreement between the measured and the reference waveforms were obtained. Negligible motion was observed in the lateral (R/L) direction. Conclusion: Our investigation demonstrates that this piezo-electric motor design is effective at simulating periodic motion and is a potential candidate for MRI-radiotherapy respiratory motion simulation. Future work should focus on evaluating non-sinusoidal waveforms, fast 3D pulse sequences, and perform dosimetric QA.

60 APPLIED LIFE SCIENCES↗

Preliminary results on the feasibility of using ultrasound to monitor intrafractional motion during radiation therapy for pancreatic cancer

Purpose: Substantial intrafraction organ motion during radiation therapy (RT) for pancreatic cancer is well recognized as a major limiting factor for accurate delivery of RT. The aim of this work is to determine the feasibility of monitoring the intrafractional motion of the pancreas or surrounding structures using ultrasound for RT delivery. Methods: Transabdominal ultrasound (TAUS) and 4DCT data were acquired on ten pancreatic cancer patients during radiation therapy process in a prospective study. In addition, TAUS and MRI were collected for five healthy volunteers. The portal vein (PV) and the head of the pancreas (HP) along with other structures were contoured on these images. Volume changes, distance between the HP and PV, and motion difference between the HP and PV were measured to examine whether PV can be used as a motion surrogate for HP. TAUS images were acquired and processed using a research version of the Clarity autoscan ultrasound system (CAUS). Motion monitoring was performed with the ultrasound probe mounted on an arm fixed to the couch. Video segments of the monitoring sessions were captured. Results: On TAUS, PV is better visualized than HP. The measured mean volume deviation for all patients for the HP and PV was 1.4 and 0.6 ml, respectively. The distance between the HP and PV was close to a constant with 0.22 mm mean deviation throughout the ten breathing phases. The mean of the absolute motion difference for all patients was 1.7 ± 0.8 mm in LR, 1.5 ± 0.5 mm in AP, and 2.3 ± 0.7 mm in SI, suggesting that the PV is a good surrogate for HP motion estimation. By using this surrogate, the HP motion tracking using TAUS was demonstrated. Conclusions: Large intrafractional organ motion due to respiratory and/or bowel motion is a limiting factor in administering curative radiation doses to pancreatic tumors. The authors investigate the use of real-time ultrasound to track pancreas motion. Due to the poor visibility of the pancreas head on an ultrasound image, the portal vein is identified as a surrogate. The authors have demonstrated the feasibility of tracking HP motion through the localization of the PV using TAUS. This will potentially allow real-time tracking of intrafractional motion to justify small PTV-margins and to account for unusual motions, thus, improving normal tissue sparing.

60 APPLIED LIFE SCIENCES↗