Predicting respiratory motion for four-dimensional radiotherapy.

Adapting radiation delivery to respiratory motion is made possible through corrective action based on real-time feedback of target position during respiration. The advantage of this approach lies with its ability to allow tighter margins around the target while simultaneously following its motion. A significant hurdle to the successful implementation of real-time target-tracking-based radiation delivery is the existence of a finite time delay between the acquisition of target position and the mechanical response of the system to the change in position. Target motion during the time delay leads to a resultant lag in the system's response to a change in tumor position. Predicting target position in advance is one approach to ensure accurate delivery. The aim of this manuscript is to estimate the predictive ability of sinusoidal and adaptive filter-based prediction algorithms on multiple sessions of patient respiratory patterns. Respiratory motion information was obtained from recordings of diaphragm motion for five patients over 60 sessions. A prediction algorithm that employed both prediction models-the sinusoidal model and the adaptive filter model-was developed to estimate prediction accuracy over all the sessions. For each session, prediction error was computed for several time instants (response time) in the future (0-1.8 seconds at 0.2-second intervals), based on position data collected over several signal-history lengths (1-7 seconds at 1-second intervals). Based on patient data included in this study, the following observations are made. Qualitative comparison of predicted and actual position indicated a progressive increase in prediction error with an increase in response time. A signal-history length of 5 seconds was found to be the optimal signal history length for prediction using the sinusoidal model for all breathing training modalities. In terms of overall error in predicting respiratory motion, the adaptive filter model performed better than the sinusoidal model. With the adaptive filter, average prediction errors of less than 0.2 cm (1sigma) are possible for response times less than 0.4 seconds. In comparing prediction error with system latency error (no prediction), the adaptive filter model exhibited lesser prediction errors as compared to the sinusoidal model, especially for longer response time values (>0.4 seconds). At smaller response time values (<0.4 seconds), improvements in prediction error reduction are required for both predictive models in order to maximize gains in position accuracy due to prediction. Respiratory motion patterns are inherently complex in nature. While linear prediction-based prediction models perform satisfactorily for shorter response times, their prediction accuracy significantly deteriorates for longer response times. Successful implementation of real-time target-tracking-based radiotherapy requires response times less than 0.4 seconds or improved prediction algorithms.

[1]  Monson H. Hayes,et al.  Statistical Digital Signal Processing and Modeling , 1996 .

[2]  D A Jaffray,et al.  The effects of intra-fraction organ motion on the delivery of dynamic intensity modulation. , 1998, Physics in medicine and biology.

[3]  R. Mohan,et al.  Quantifying the predictability of diaphragm motion during respiration with a noninvasive external marker. , 2003, Medical physics.

[4]  J. Adler,et al.  Robotic Motion Compensation for Respiratory Movement during Radiosurgery , 2000, Computer aided surgery : official journal of the International Society for Computer Aided Surgery.

[5]  P Zygmanski,et al.  A method of calculating a lung clinical target volume DVH for IMRT with intrafractional motion. , 2003, Medical physics.

[6]  Steve B. Jiang,et al.  An experimental investigation on intra-fractional organ motion effects in lung IMRT treatments. , 2003, Physics in medicine and biology.

[7]  L. Kochian Author to whom correspondence should be addressed , 2006 .

[8]  R. Mohan,et al.  Quantifying the effect of intrafraction motion during breast IMRT planning and dose delivery. , 2003, Medical physics.

[9]  Linda Hong,et al.  The effects of intra-fraction organ motion on the delivery of intensity-modulated field with a multileaf collimator. , 2003, Medical physics.

[10]  L. Papie,et al.  The leaf sweep algorithm for an immobile and moving target as an optimal control problem in radiotherapy delivery , 2003 .

[11]  Steve B. Jiang,et al.  Effects of intra-fraction motion on IMRT dose delivery: statistical analysis and simulation. , 2002, Physics in medicine and biology.

[12]  D. Menicucci,et al.  Deriving the respiratory sinus arrhythmia from the heartbeat time series using empirical mode decomposition , 2003, q-bio/0310002.

[13]  Dorin Todor,et al.  A singular value decomposition based algorithm for multicomponent exponential fitting of NMR relaxation signals , 1995 .

[14]  Steve B Jiang,et al.  Synchronized moving aperture radiation therapy (SMART): average tumour trajectory for lung patients. , 2003, Physics in medicine and biology.

[15]  R. Mohan,et al.  Motion adaptive x-ray therapy: a feasibility study , 2001, Physics in medicine and biology.

[16]  Gregory C Sharp,et al.  Prediction of respiratory tumour motion for real-time image-guided radiotherapy. , 2004, Physics in medicine and biology.