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 253 records · Page 14

Adaptation of the modified Bouc–Wen model to compensate for hysteresis in respiratory motion for the list-mode binning of cardiac SPECT and PET acquisitions: Testing using MRI

Purpose: Binning list-mode acquisitions as a function of a surrogate signal related to respiration has been employed to reduce the impact of respiratory motion on image quality in cardiac emission tomography (SPECT and PET). Inherent in amplitude binning is the assumption that there is a monotonic relationship between the amplitude of the surrogate signal and respiratory motion of the heart. This assumption is not valid in the presence of hysteresis when heart motion exhibits a different relationship with the surrogate during inspiration and expiration. The purpose of this study was to investigate the novel approach of using the Bouc–Wen (BW) model to provide a signal accounting for hysteresis when binning list-mode data with the goal of thereby improving motion correction. The study is based on the authors’ previous observations that hysteresis between chest and abdomen markers was indicative of hysteresis between abdomen markers and the internal motion of the heart. Methods: In 19 healthy volunteers, they determined the internal motion of the heart and diaphragm in the superior–inferior direction during free breathing using MRI navigators. A visual tracking system (VTS) synchronized with MRI acquisition tracked the anterior–posterior motions of external markers placed on the chest and abdomen. These data were employed to develop and test the Bouc–Wen model by inputting the VTS derived chest and abdomen motions into it and using the resulting output signals as surrogates for cardiac motion. The data of the volunteers were divided into training and testing sets. The training set was used to obtain initial values for the model parameters for all of the volunteers in the set, and for set members based on whether they were or were not classified as exhibiting hysteresis using a metric derived from the markers. These initial parameters were then employed with the testing set to estimate output signals. Pearson’s linear correlation coefficient between the abdomen, chest, average of chest and abdomen markers, and Bouc–Wen derived signals versus the true internal motion of the heart from MRI was used to judge the signals match to the heart motion. Results: The results show that the Bouc–Wen model generated signals demonstrated strong correlation with the heart motion. This correlation was slightly larger on average than that of the external surrogate signals derived from the abdomen marker, and average of the abdomen and chest markers, but was not statistically significantly different from them. Conclusions: The results suggest that the proposed model has the potential to be a unified framework for modeling hysteresis in respiratory motion in cardiac perfusion studies and beyond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accelerated acquisition of tagged MRI for cardiac motion correction in simultaneous PET-MR: Phantom and patient studies

Purpose: Degradation of image quality caused by cardiac and respiratory motions hampers the diagnostic quality of cardiac PET. It has been shown that improved diagnostic accuracy of myocardial defect can be achieved by tagged MR (tMR) based PET motion correction using simultaneous PET-MR. However, one major hurdle for the adoption of tMR-based PET motion correction in the PET-MR routine is the long acquisition time needed for the collection of fully sampled tMR data. In this work, the authors propose an accelerated tMR acquisition strategy using parallel imaging and/or compressed sensing and assess the impact on the tMR-based motion corrected PET using phantom and patient data. Methods: Fully sampled tMR data were acquired simultaneously with PET list-mode data on two simultaneous PET-MR scanners for a cardiac phantom and a patient. Parallel imaging and compressed sensing were retrospectively performed by GRAPPA and kt-FOCUSS algorithms with various acceleration factors. Motion fields were estimated using nonrigid B-spline image registration from both the accelerated and fully sampled tMR images. The motion fields were incorporated into a motion corrected ordered subset expectation maximization reconstruction algorithm with motion-dependent attenuation correction. Results: Although tMR acceleration introduced image artifacts into the tMR images for both phantom and patient data, motion corrected PET images yielded similar image quality as those obtained using the fully sampled tMR images for low to moderate acceleration factors (<4). Quantitative analysis of myocardial defect contrast over ten independent noise realizations showed similar results. It was further observed that although the image quality of the motion corrected PET images deteriorates for high acceleration factors, the images were still superior to the images reconstructed without motion correction. Conclusions: Accelerated tMR images obtained with more than 4 times acceleration can still provide relatively accurate motion fields and yield tMR-based motion corrected PET images with similar image quality as those reconstructed using fully sampled tMR data. The reduction of tMR acquisition time makes it more compatible with routine clinical cardiac PET-MR studies.

60 APPLIED LIFE SCIENCES↗

MagicPlate-512: A 2D silicon detector array for quality assurance of stereotactic motion adaptive radiotherapy

Purpose: Spatial and temporal resolutions are two of the most important features for quality assurance instrumentation of motion adaptive radiotherapy modalities. The goal of this work is to characterize the performance of the 2D high spatial resolution monolithic silicon diode array named “MagicPlate-512” for quality assurance of stereotactic body radiation therapy (SBRT) and stereotactic radiosurgery (SRS) combined with a dynamic multileaf collimator (MLC) tracking technique for motion compensation. Methods: MagicPlate-512 is used in combination with the movable platform HexaMotion and a research version of radiofrequency tracking system Calypso driving MLC tracking software. The authors reconstruct 2D dose distributions of small field square beams in three modalities: in static conditions, mimicking the temporal movement pattern of a lung tumor and tracking the moving target while the MLC compensates almost instantaneously for the tumor displacement. Use of Calypso in combination with MagicPlate-512 requires a proper radiofrequency interference shielding. Impact of the shielding on dosimetry has been simulated by GEANT4 and verified experimentally. Temporal and spatial resolutions of the dosimetry system allow also for accurate verification of segments of complex stereotactic radiotherapy plans with identification of the instant and location where a certain dose is delivered. This feature allows for retrospective temporal reconstruction of the delivery process and easy identification of error in the tracking or the multileaf collimator driving systems. A sliding MLC wedge combined with the lung motion pattern has been measured. The ability of the MagicPlate-512 (MP512) in 2D dose mapping in all three modes of operation was benchmarked by EBT3 film. Results: Full width at half maximum and penumbra of the moving and stationary dose profiles measured by EBT3 film and MagicPlate-512 confirm that motion has a significant impact on the dose distribution. Motion, no motion, and motion with MLC tracking profiles agreed within 1 and 0.4 mm, respectively, for all field sizes tested. Use of electromagnetic tracking system generates a fluctuation of the detector baseline up to 10% of the full scale signal requiring a proper shielding strategy. MagicPlate-512 is also able to reconstruct the dose variation pulse-by-pulse in each pixel of the detector. An analysis of the dose transients with motion and motion with tracking shows that the tracking feedback algorithm used for this experiment can compensate effectively only the effect of the slower transient components. The fast changing components of the organ motion can contribute only to discrepancy of the order of 15% in penumbral region while the slower components can change the dose profile up to 75% of the expected dose. Conclusions: MagicPlate-512 is shown to be, potentially, a valid alternative to film or 2D ionizing chambers for quality assurance dosimetry in SRS or SBRT. Its high spatial and temporal resolutions allow for accurate reconstruction of the profile in any conditions with motion and with tracking of the motion. It shows excellent performance to reconstruct the dose deposition in real time or retrospectively as a function of time for detailed analysis of the effect of motion in a specific pixel or area of interest.

60 APPLIED LIFE SCIENCES↗

SU-E-J-60: Evaluation of Temporal Lag in Radiotherapy Gating for Tumor Motion Trajectories

Purpose: Evaluating timing differences between LINAC beam ON/OFF and the estimation of tumor positioning using gating systems is essential for establishing confidence when treating with gating during radiotherapy, and is an annual requirement of TG-142. Temporal discrepancies between the trajectories of external marker surrogates and beam delivery may vary depending upon the type of external marker motion, which is quantified in this work for several trajectories. Methods: A precise robotic 3D motion stage performed several trajectories typically used for gating phantoms, including sinusoidal and Lujan-type motion; a commercial respiratory motion simulator was also employed. The true motions were monitored using variable resistors. The beam ON/OFF was controlled separately by two RPM (Varian) systems, an integrated version delivered by a Varian Truebeam LINAC and version 1.6 delivered by a Varian Trilogy, and measured using a diode. The resistor and diode signals were read by a multichannel digital oscilloscope, and timing differences between beam ON/OFF as measured by the diode and the phantom motion were determined using a peak detection algorithm. Results: Timing differences between beam ON/OFF and 3D stage motion peaks (diode—true motion timing) were computed to be 79.4 & 57.7ms for sinusoidal motion and 109.1 & 63.6ms for Lujan-type motion on the Truebeam LINAC, for beam ON and OFF, respectively. Timing differences for the Trilogy LINAC were 34.4 & 55.2ms for the sinusoidal motion and 29.0 & 26.3ms for the Lujan-type motion, for beam ON and OFF, respectively. With the commercial motion simulator, the timing differences were found to be −9ms and −78ms for beam ON/OFF, respectively, with the Truebeam, and −97.6ms and −60.9ms for beam ON/OFF, respectively, with the Trilogy. Conclusion: Setup-dependent temporal lags were found using this methodology. These discrepancies have the potential to influence quality assurance on gating systems and ultimately the treatment dose during radiotherapy. NIH National Institute of Biomedical Imaging and Bioengineering grant T32 EB002103-21.

60 APPLIED LIFE SCIENCES↗

Characterizing spatiotemporal information loss in sparse-sampling-based dynamic MRI for monitoring respiration-induced tumor motion in radiotherapy

Purpose: Sparse-sampling and reconstruction techniques represent an attractive strategy to achieve faster image acquisition speeds, while maintaining adequate spatial resolution and signal-to-noise ratio in rapid magnetic resonance imaging (MRI). The authors investigate the use of one such sequence, broad-use linear acquisition speed-up technique (k-t BLAST) in monitoring tumor motion for thoracic and abdominal radiotherapy and examine the potential trade-off between increased sparsification (to increase imaging speed) and the potential loss of “true” information due to greater reliance on a priori information. Methods: Lung tumor motion trajectories in the superior–inferior direction, previously recorded from ten lung cancer patients, were replayed using a motion phantom module driven by an MRI-compatible motion platform. Eppendorf test tubes filled with water which serve as fiducial markers were placed in the phantom. The modeled rigid and deformable motions were collected in a coronal image slice using balanced fast field echo in conjunction with k-t BLAST. Root mean square (RMS) error was used as a metric of spatial accuracy as measured trajectories were compared to input data. The loss of spatial information was characterized for progressively increasing acceleration factor from 1 to 16; the resultant sampling frequency was increased approximately from 2.5 to 19 Hz when the principal direction of the motion was set along frequency encoding direction. In addition to the phantom study, respiration-induced tumor motions were captured from two patients (kidney tumor and lung tumor) at 13 Hz over 49 s to demonstrate the impact of high speed motion monitoring over multiple breathing cycles. For each subject, the authors compared the tumor centroid trajectory as well as the deformable motion during free breathing. Results: In the rigid and deformable phantom studies, the RMS error of target tracking at the acquisition speed of 19 Hz was approximately 0.3–0.4 mm, which was smaller than the reconstructed pixel resolution of 0.67 mm. In the patient study, the dynamic 2D MRI enabled the monitoring of cycle-to-cycle respiratory variability present in the tumor position. It was seen that the range of centroid motion as well as the area covered due to target motion during each individual respiratory cycle was underestimated compared to the entire motion range observed over multiple breathing cycles. Conclusions: The authors’ initial results demonstrate that sparse-sampling- and reconstruction-based dynamic MRI can be used to achieve adequate image acquisition speeds without significant information loss for the task of radiotherapy guidance. Such monitoring can yield spatial and temporal information superior to conventional offline and online motion capture methods used in thoracic and abdominal radiotherapy.

60 APPLIED LIFE SCIENCES↗

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↗

SU-F-J-119: Pilot Study On the Location-Based Lung Motion Assessment

Purpose: In most of lung treatment cases with various radiotherapy beam modalities, 4DCT images are obtained in order to define ITV. ITV is defined with the signal from motion monitoring system, e.g. RPM. However, the signal is not consistent with tumor motion because it varies with location, its size, age, gender, etc. In the present study, the location-based motion assessment is presented. Methods: 4DCT images of 70 patients were reviewed: 28-left-lung and 42-right-lung patients; 36-female and 34-male patients; the age range of 51.2–89.9; tumor-size range of 0.75–9.50cm with 25% of these adherent to bony-anatomy. Philips Big-Bore Simulation CT and RPM systems were used. The study was performed as follows. First, RPM signal and tumor motion in superior-inferior direction was compared. Second, the tumor size and its motion amplitude in all directions were measured at multiple locations. Third, the average tumor motion was calculated to assess general motion amplitudes at various locations. Results: RPM amplitude is not consistent with lung tumor motion amplitude. The tumors of similar sizes at similar location present various motion amplitude up to 1.1cm difference, but in average, the standard deviation was <0.5cm. Almost regardless of tumor sizes, the tumor motion was greatest at lower lobe location (>=1.0cm), and the smallest at upper lobe location and when adherent to bony-anatomy (<=0.5cm). Conclusion: The tumor size affects the motion amplitude less than does the tumor location. However, as the study results indicate that tumor motion has noticeable variation and so further study with more patient cases is needed. Also, for the same patient, the RPM signal presents instability of breathing, and clinically the patient with the instability of RPM breathing of <=10% is selected for respiratory-gated radiotherapy and ∼25% of patients under current study was treated. Patient-specific motion-uncertainty margins are considered to be added following further study.

60 APPLIED LIFE SCIENCES↗

SU-F-J-128: Dosimetric Impact of Esophagus Motion in Spine Stereotactic Body Radiotherapy

Purpose: Acute esophageal toxicity is a common side effect in spine stereotactic body radiotherapy (SBRT). The respiratory motion may alter esophageal position from the planning scan resulting in excessive esophageal dose. Here we assessed the dosimetric impact resulting from the esophageal motion using 4DCT. Methods: Nine patients treated to their thoracic spines in one fraction of 24 Gy were identified for this study. The original plan on a free breathing CT was copied to each phase image of a 4DCT scan, recalculated, scaled, and accumulated to the free breathing CT using deformable image registration. A segment of esophagus was contoured in the vicinity of treatment target. Esophagus dose volume histogram (DVH) was generated for both the original planned dose and the accumulated 4D dose for comparison. In parallel, we performed a chained deformable registration of 4DCT phase images to estimate the motion magnitude of the esophagus in a breathing cycle. We examined the correlation between the motion magnitude and the dosimetric deviation. Results: The esophageal motion mostly exhibited in the superior-inferior direction. The cross-sectional motion was small. Esophagus motion at T1 vertebra level (0.7 mm) is much smaller than that at T11 vertebra level (6.5 mm). The difference of Dmax between the original and 4D dose distributions ranged from 9.1 cGy (esophagus motion: 5.6 mm) to 231.1 cGy (esophagus motion: 3.1 mm). The difference of D(5cc) ranged from 5 cGy (esophagus motion: 3.1 mm) to 85 cGy (esophagus motion: 3.3 mm). There was no correlation between the dosimetric deviation and the motion magnitude. The V(11.9Gy)<5cc constraint was met for each patient when examining the DVH calculated from the 4D dose. Conclusion: Respiratory motion did not result in substantial dose increase to esophagus in spine SBRT. 4DCT simulation may not be necessary with regards to esophageal dose assessment.

60 APPLIED LIFE SCIENCES↗

SU-G-BRA-13: An Advanced Deformable Lung Phantom for Analyzing the Dosimetric Impact of Respiratory Motion

Purpose: The difference between three-dimensional (3D) and four-dimensional (4D) dose is affected by factors such as tumor size and motion. To quantitatively analyze the effects of these factors, a phantom that can independently control for each factor is required. The purpose of this study is to develop a deformable lung phantom with the above attributes and evaluate characteristics. Methods: A phantom was designed to simulate diaphragm motion with amplitude in the range 1 to 7 cm and various periods of regular breathing. To simulate different size tumors, tumors were produced by pouring liquid silicone into custom molds created by a 3D printer. The accuracy of phantom diaphragm motion was assessed using calipers and protractor. To control tumor motion, tumor trajectories were evaluated using 4D computed tomography (CT), and diaphragm-tumor correlation curve was calculated by curve fitting method. Three-dimensional dose and 4D dose were calculated and compared according to tumor motion. Results: The accuracy of phantom diaphragm motion was less than 1 mm. Maximum tumor motion amplitudes in the left-right and anterior-posterior directions were 0.08 and 0.12 cm, respectively, in a 10 cm{sup 3} tumor, and 0.06 and 0.27 cm, respectively, in a 90 cm{sup 3} tumor. The diaphragm-tumor correlation curve showed that tumor motion in the superior-inferior direction was increased with increasing diaphragm motion. In the 10 cm{sup 3} tumor, the tumor motion was larger than the 90 cm{sup 3} tumor. According to tumor motion, variation of dose difference between 3D and 4D was identified. Conclusion: The developed phantom can independently control factors such as tumor size and motion. In potentially, this phantom can be used to quantitatively analyze the dosimetric impact of respiratory motion according to the factors that influence the difference between 3D and 4D dose. This research was supported by the Mid-career Researcher Program through NRF funded by the Ministry of Science, ICT & Future Planning of Korea (NRF-2014R1A2A1A10050270) and by the Radiation Technology R&D program through the National Research Foundation of Korea funded by the Ministry of Science, ICT & Future Planning (No. 2013M2A2A7038291)

60 APPLIED LIFE SCIENCES↗

MO-FG-BRA-09: Towards an Optimal Breath-Holding Procedure for Radiotherapy: Differences in Organ Motion During Inhalation and Exhalation Breath-Holds

Purpose: Breath-holding (BH) is often used to reduce organ motion during radiotherapy. The aim of this study was to determine the differences in pancreatic and diaphragmatic motion during BH between inhalation and exhalation BHs with variable lung volumes and to investigate whether motion increases/decreases during BH. Methods: Sixteen healthy volunteers were asked to perform four different 60-second BHs, from fully inflated to fully deflated lungs (i.e. lung volumes of: 100%, ∼70%, ∼30% and 0% of inspiratory capacity) three times (total of 192 BHs). During each BH, we obtained single-slice (coronal) magnetic-resonance scans with spatial resolution 0.93×0.93×8.0 mm3 and temporal resolution 0.6 s. We used 2-dimensional image correlation to obtain the motion of pancreatic head and diaphragm during BH. Motion magnitude in inferior-superior direction was obtained by determining the maximum displacement during BH. Results: Pancreatic and diaphragmatic drifts occurred during BH and were mostly in the superior direction. We observed significantly smaller pancreatic and diaphragmatic motion magnitudes in inferior-superior direction during exhalation BHs (BH{sub 30%} and BH{sub 0%}) compared to inhalation BHs (BH{sub 100%} and BH{sub 70%}). The mean motion magnitudes of the pancreatic head were 7.0, 6.5, 4.4 and 4.2 mm during BH{sub 100%}, BH{sub 70%}, BH{sub 30%} and BH{sub 0%}, respectively, and mean BH durations were 59.9, 59.1, 59.0 and 52.7 s. For the diaphragm, mean motion magnitudes were 9.8, 9.0, 5.6 and 4.3 mm, respectively. When considering 30-second BHs, as often used in the clinic, the motion was most pronounced during the first 10 s and excluding these from the analysis (yielding an effective BH period of 20 s) significantly reduced (P≤0.002) organ motion. Conclusion: Organ motion was significantly smaller during exhalation BHs compared to inhalation BHs. Also, motion was largest at the start of BH. Hence, waiting for 10 s may significantly decrease motion of the pancreas and diaphragm during treatment.

60 APPLIED LIFE SCIENCES↗

DISTANCE AND PROPER MOTION MEASUREMENT OF THE RED SUPERGIANT, PZ CAS, IN VERY LONG BASELINE INTERFEROMETRY H{sub 2}O MASER ASTROMETRY

We present the very long baseline interferometry H{sub 2}O maser monitoring observations of the red supergiant, PZ Cas, at 12 epochs from 2006 April to 2008 May. We fitted maser motions to a simple model composed of a common annual parallax and linear motions of the individual masers. The maser motions with the parallax subtracted were well modeled by a combination of a common stellar proper motion and a radial expansion motion of the circumstellar envelope. We obtained an annual parallax of 0.356 {+-} 0.026 mas and a stellar proper motion of {mu}{sub {alpha}}{sup *} cos {delta} = -3.7 {+-} 0.2 and {mu}{sup *}{sub {delta}}=-2.0{+-}0.3 mas yr{sup -1} eastward and northward, respectively. The annual parallax corresponds to a trigonometric parallax of 2.81{sup +0.22}{sub -0.19} kpc. By rescaling the luminosity of PZ Cas in any previous studies using our trigonometric parallax, we estimated the location of PZ Cas on a Hertzsprung-Russell diagram and found that it approaches a theoretically evolutionary track around an initial mass of {approx}25 M{sub Sun }. The sky position and the distance to PZ Cas are consistent with the OB association, Cas OB5, which is located in a molecular gas super shell. The proper motion of PZ Cas is close to that of the OB stars and other red supergiants in Cas OB5 measured by the Hipparcos satellite. We derived the peculiar motion of PZ Cas of U{sub s} = 22.8 {+-} 1.5, V{sub s} = 7.1 {+-} 4.4, and W{sub s} = -5.7 {+-} 4.4 km s{sup -1}. This peculiar motion has rather a large U{sub s} component, unlike those of near high-mass star-forming regions with negatively large V{sub s} motions. The uniform proper motions of the Cas OB5 member stars suggest random motions of giant molecular clouds moving into local potential minima in a time-dependent spiral arm, rather than a velocity field caused by the spiral arm density wave.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of time delays in gated radiotherapy for realistic respiratory motions

Purpose: Gated radiotherapy is used to reduce internal motion margins, escalate target dose, and limit normal tissue dose; however, its temporal accuracy is limited. Beam-on and beam-off time delays can lead to treatment inefficiencies and/or geographic misses; therefore, AAPM Task Group 142 recommends verifying the temporal accuracy of gating systems. Many groups use sinusoidal phantom motion for this, under the tacit assumption that use of sinusoidal motion for determining time delays produces negligible error. The authors test this assumption by measuring gating time delays for several realistic motion shapes with increasing degrees of irregularity. Methods: Time delays were measured on a linear accelerator with a real-time position management system (Varian TrueBeam with RPM system version 1.7.5) for seven motion shapes: regular sinusoidal; regular realistic-shape; large (40%) and small (10%) variations in amplitude; large (40%) variations in period; small (10%) variations in both amplitude and period; and baseline drift (30%). Film streaks of radiation exposure were generated for each motion shape using a programmable motion phantom. Beam-on and beam-off time delays were determined from the difference between the expected and observed streak length. Results: For the system investigated, all sine, regular realistic-shape, and slightly irregular amplitude variation motions had beam-off and beam-on time delays within the AAPM recommended limit of less than 100 ms. In phase-based gating, even small variations in period resulted in some time delays greater than 100 ms. Considerable time delays over 1 s were observed with highly irregular motion. Conclusions: Sinusoidal motion shapes can be considered a reasonable approximation to the more complex and slightly irregular shapes of realistic motion. When using phase-based gating with predictive filters even small variations in period can result in time delays over 100 ms. Clinical use of these systems for patients with highly irregular patterns of motion is not advised due to large beam-on and beam-off time delays.

60 APPLIED LIFE SCIENCES↗

Visual motion integration for perception and pursuit

To examine the relationship between visual motion processing for perception and pursuit, we measured the pursuit eye-movement and perceptual responses to the same complex-motion stimuli. We show that humans can both perceive and pursue the motion of line-figure objects, even when partial occlusion makes the resulting image motion vastly different from the underlying object motion. Our results show that both perception and pursuit can perform largely accurate motion integration, i.e. the selective combination of local motion signals across the visual field to derive global object motion. Furthermore, because we manipulated perceived motion while keeping image motion identical, the observed parallel changes in perception and pursuit show that the motion signals driving steady-state pursuit and perception are linked. These findings disprove current pursuit models whose control strategy is to minimize retinal image motion, and suggest a new framework for the interplay between visual cortex and cerebellum in visuomotor control.

NASA Center ARC↗

The Non-Gaussian Nature of Prostate Motion Based on Real-Time Intrafraction Tracking

Purpose: The objective of this work is to test the validity of the Gaussian approximation for prostate motion through characterization of its spatial distribution. Methods and Materials: Real-time intrafraction prostate motion was observed using Calypso 4-dimensional (4D) nonradioactive electromagnetic tracking system. We report the results from a total of 1024 fractions from 31 prostate cancer patients. First, the correlation of prostate motion in right/left (RL), anteroposterior (AP), and superoinferior (SI) direction were determined using Pearson's correlation of coefficient. Then the spatial distribution of prostate motion was analyzed for individual fraction, individual patient including all fractions, and all patients including all fractions. The displacement in RL, AP, SI, oblique, or total direction is fitted into a Gaussian distribution, and a Lilliefors test was used to evaluate the validity of the hypothesis that the displacement is normally distributed. Results: There is high correlation in AP/SI direction (61% of fractions with medium or strong correlation). This is consistent with the longitudinal oblique motion of the prostate, and likely the effect from respiration on an organ confined within the genitourinary diaphragm with the rectum sitting posteriorly and bladder sitting superiorly. In all directions, the non-Gaussian distribution is more common for individual fraction, individual patient including all fractions, and all patients including all fractions. The spatial distribution of prostate motion shows an elongated shape in oblique direction, indicating a higher range of motion in the AP and SI directions. Conclusions: Our results showed that the prostate motion is highly correlated in AP and SI direction, indicating an oblique motion preference. In addition, the spatial distribution of prostate motion is elongated in an oblique direction, indicating that the organ motion dosimetric modeling using Gaussian kernel may need to be modified to account for the particular organ motion character of prostate.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Dynamic volume vs respiratory correlated 4DCT for motion assessment in radiation therapy simulation

Purpose: Conventional (i.e., respiratory-correlated) 4DCT exploits the repetitive nature of breathing to provide an estimate of motion; however, it has limitations due to binning artifacts and irregular breathing in actual patient breathing patterns. The aim of this work was to evaluate the accuracy and image quality of a dynamic volume, CT approach (4D{sub vol}) using a 320-slice CT scanner to minimize these limitations, wherein entire image volumes are acquired dynamically without couch movement. This will be compared to the conventional respiratory-correlated 4DCT approach (RCCT). Methods: 4D{sub vol} CT was performed and characterized on an in-house, programmable respiratory motion phantom containing multiple geometric and morphological ''tumor'' objects over a range of regular and irregular patient breathing traces obtained from 3D fluoroscopy and compared to RCCT. The accuracy of volumetric capture and breathing displacement were evaluated and compared with the ground truth values and with the results reported using RCCT. A motion model was investigated to validate the number of motion samples needed to obtain accurate motion probability density functions (PDF). The impact of 4D image quality on this accuracy was then investigated. Dose measurements using volumetric and conventional scan techniques were also performed and compared. Results: Both conventional and dynamic volume 4DCT methods were capable of estimating the programmed displacement of sinusoidal motion, but patient breathing is known to not be regular, and obvious differences were seen for realistic, irregular motion. The mean RCCT amplitude error averaged at 4 mm (max. 7.8 mm) whereas the 4D{sub vol} CT error stayed below 0.5 mm. Similarly, the average absolute volume error was lower with 4D{sub vol} CT. Under irregular breathing, the 4D{sub vol} CT method provides a close description of the motion PDF (cross-correlation 0.99) and is able to track each object, whereas the RCCT method results in a significantly different PDF from the ground truth, especially for smaller tumors (cross-correlation ranging between 0.04 and 0.69). For the protocols studied, the dose measurements were higher in the 4D{sub vol} CT method (40%), but it was shown that significant mAs reductions can be achieved by a factor of 4-5 while maintaining image quality and accuracy. Conclusions: 4D{sub vol} CT using a scanner with a large cone-angle is a promising alternative for improving the accuracy with which respiration-induced motion can be characterized, particularly for patients with irregular breathing motion. This approach also generates 4DCT image data with a reduced total scan time compared to a RCCT scan, without the need for image binning or external respiration signals within the 16 cm scan length. Scan dose can be made comparable to RCCT by optimization of the scan parameters. In addition, it provides the possibility of measuring breathing motion for more than one breathing cycle to assess stability and obtain a more accurate motion PDF, which is currently not feasible with the conventional RCCT approach.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Effectiveness of external respiratory surrogates for in vivo liver motion estimation

Purpose: Due to low frame rate of MRI and high radiation damage from fluoroscopy and CT, liver motion estimation using external respiratory surrogate signals seems to be a better approach to track liver motion in real-time for liver tumor treatments in radiotherapy and thermotherapy. This work proposes a liver motion estimation method based on external respiratory surrogate signals. Animal experiments are also conducted to investigate related issues, such as the sensor arrangement, multisensor fusion, and the effective time period. Methods: Liver motion and abdominal motion are both induced by respiration and are proved to be highly correlated. Contrary to the difficult direct measurement of the liver motion, the abdominal motion can be easily accessed. Based on this idea, our study is split into the model-fitting stage and the motion estimation stage. In the first stage, the correlation between the surrogates and the liver motion is studied and established via linear regression method. In the second stage, the liver motion is estimated by the surrogate signals with the correlation model. Animal experiments on cases of single surrogate signal, multisurrogate signals, and long-term surrogate signals are conducted and discussed to verify the practical use of this approach. Results: The results show that the best single sensor location is at the middle of the upper abdomen, while multisurrogate models are generally better than the single ones. The estimation error is reduced from 0.6 mm for the single surrogate models to 0.4 mm for the multisurrogate models. The long-term validity of the estimation models is quite satisfactory within the period of 10 min with the estimation error less than 1.4 mm. Conclusions: External respiratory surrogate signals from the abdomen motion produces good performance for liver motion estimation in real-time. Multisurrogate signals enhance estimation accuracy, and the estimation model can maintain its accuracy for at least 10 min. This approach can be used in practical applications such as the liver tumor treatment in radiotherapy and thermotherapy.

60 APPLIED LIFE SCIENCES↗

An externally and internally deformable, programmable lung motion phantom

Purpose: Most clinically deployed strategies for respiratory motion management in lung radiotherapy (e.g., gating and tracking) use external markers that serve as surrogates for tumor motion. However, typical lung phantoms used to validate these strategies are based on a rigid exterior and a rigid or a deformable-interior. Such designs do not adequately represent respiration because the thoracic anatomy deforms internally as well as externally. In order to create a closer approximation of respiratory motion, the authors describe the construction and experimental testing of an externally as well as internally deformable, programmable lung phantom. Methods: The outer shell of a commercially available lung phantom (RS-1500, RSD, Inc.) was used. The shell consists of a chest cavity with a flexible anterior surface, and embedded vertebrae, rib-cage and sternum. A custom-made insert was designed using a piece of natural latex foam block. A motion platform was programmed with sinusoidal and ten patient-recorded lung tumor trajectories. The platform was used to drive a rigid foam “diaphragm” that compressed/decompressed the phantom interior. Experimental characterization comprised of determining the reproducibility and the external–internal correlation of external and internal marker trajectories extracted from kV x-ray fluoroscopy. Experiments were conducted to illustrate three example applications of the phantom—(i) validating the geometric accuracy of the VisionRT surface photogrammetry system; (ii) validating an image registration tool, NiftyReg; and (iii) quantifying the geometric error due to irregular motion in four-dimensional computed tomography (4DCT). Results: The phantom correctly reproduced sinusoidal and patient-derived motion, as well as realistic respiratory motion-related effects such as hysteresis. The reproducibility of marker trajectories over multiple runs for sinusoidal as well as patient traces, as characterized by fluoroscopy, was within 0.25 mm RMS error. The motion trajectories of internal and external radio-opaque markers as measured by fluoroscopy were found to be highly correlated (R > 0.95). Using the phantom, it was demonstrated that the motion trajectories of regions-of-interest on the surface as measured by VisionRT are highly consistent with corresponding fluoroscopically acquired surface marker trajectories, with RMS errors within 0.26 mm. Furthermore, it was shown that the trajectories of external and internal marker trajectories derived from NiftyReg deformation vector fields were within 1 mm root mean square errors comparing to trajectories obtained by segmenting markers from individual fluoro frames. Finally, it was shown that while 4DCT can be used to localize internal markers for sinusoidal motion with reasonable accuracy, the localization error increases significantly (by a factor of ∼2) in the presence of cycle-to-cycle variations that are observed in patient-derived respiratory motion. Conclusions: The authors have developed a realistic externally and internally deformable, programmable lung phantom that will serve as a valuable tool for clinical and investigational motion management studies in thoracic and abdominal radiation therapies.

60 APPLIED LIFE SCIENCES↗

Geometric validation of self-gating k-space-sorted 4D-MRI vs 4D-CT using a respiratory motion phantom

Purpose: MRI is increasingly being used for radiotherapy planning, simulation, and in-treatment-room motion monitoring. To provide more detailed temporal and spatial MR data for these tasks, we have recently developed a novel self-gated (SG) MRI technique with advantage of k-space phase sorting, high isotropic spatial resolution, and high temporal resolution. The current work describes the validation of this 4D-MRI technique using a MRI- and CT-compatible respiratory motion phantom and comparison to 4D-CT. Methods: The 4D-MRI sequence is based on a spoiled gradient echo-based 3D projection reconstruction sequence with self-gating for 4D-MRI at 3 T. Respiratory phase is resolved by using SG k-space lines as the motion surrogate. 4D-MRI images are reconstructed into ten temporal bins with spatial resolution 1.56 × 1.56 × 1.56 mm{sup 3}. A MRI-CT compatible phantom was designed to validate the performance of the 4D-MRI sequence and 4D-CT imaging. A spherical target (diameter 23 mm, volume 6.37 ml) filled with high-concentration gadolinium (Gd) gel is embedded into a plastic box (35 × 40 × 63 mm{sup 3}) and stabilized with low-concentration Gd gel. The phantom, driven by an air pump, is able to produce human-type breathing patterns between 4 and 30 respiratory cycles/min. 4D-CT of the phantom has been acquired in cine mode, and reconstructed into ten phases with slice thickness 1.25 mm. The 4D images sets were imported into a treatment planning software for target contouring. The geometrical accuracy of the 4D MRI and CT images has been quantified using target volume, flattening, and eccentricity. The target motion was measured by tracking the centroids of the spheres in each individual phase. Motion ground-truth was obtained from input signals and real-time video recordings. Results: The dynamic phantom has been operated in four respiratory rate (RR) settings, 6, 10, 15, and 20/min, and was scanned with 4D-MRI and 4D-CT. 4D-CT images have target-stretching, partial-missing, and other motion artifacts in various phases, whereas the 4D-MRI images are visually free of those artifacts. Volume percentage difference for the 6.37 ml target ranged from 5.3% ± 4.3% to 10.3% ± 5.9% for 4D-CT, and 1.47 ± 0.52 to 2.12 ± 1.60 for 4D-MRI. With an increase of respiratory rate, the target volumetric and geometric deviations increase for 4D-CT images while remaining stable for the 4D-MRI images. Target motion amplitude errors at different RRs were measured with a range of 0.66–1.25 mm for 4D-CT and 0.2–0.42 mm for 4D-MRI. The results of Mann–Whitney tests indicated that 4D-MRI significantly outperforms 4D-CT in phase-based target volumetric (p = 0.027) and geometric (p < 0.001) measures. Both modalities achieve equivalent accuracy in measuring motion amplitude (p = 0.828). Conclusions: The k-space self-gated 4D-MRI technique provides a robust method for accurately imaging phase-based target motion and geometry. Compared to 4D-CT, the current 4D-MRI technique demonstrates superior spatiotemporal resolution, and robust resistance to motion artifacts caused by fast target motion and irregular breathing patterns. The technique can be used extensively in abdominal targeting, motion gating, and toward implementing MRI-based adaptive radiotherapy.

60 APPLIED LIFE SCIENCES↗