Three-dimensional tracking of a target in a body
Disclosed is a method and system for three-dimensional tracking of a target located within a body, the method performed using at least one processing system. A two-dimensional scanned image of the body including the target is processed to obtain a two-dimensional image of the target. A first present dataset of the target is predicted using a previous dataset of the target and a state transition model, the first present dataset includes a three-dimensional present position value of the target. A second present dataset of the target is measured by template-matching of the two-dimensional image of the target with a model of the target. A third present dataset of the target is estimated by statistical inference using the first present dataset and the second present dataset. The previous dataset of the target is updated to match the third present dataset.
1. A method for three-dimensional tracking of a target located within a body, the method performed by at least one processing system, and comprising the steps of:
processing a two-dimensional scanned image of the body including the target to obtain a two-dimensional image of the target, wherein the two-dimensional scanned image is acquired from a scanning device;
predicting a first present dataset of the target by using a previous dataset of the target and a state transition model, wherein the first present dataset includes at least a three-dimensional present position value of the target, and wherein the previous dataset includes at least a three-dimensional previous position value of the target;
measuring a second present dataset of the target by template-matching of the two-dimensional image of the target with a model of the target, wherein the second present dataset includes at least a two-dimensional present position value of the target;
estimating a third present dataset of the target by statistical inference using the first present dataset and the second present dataset, wherein the third present dataset includes at least a three-dimensional present position value of the target; and
updating the previous dataset of the target to match the third present dataset.
2. The method of claim 1 , wherein the statistical inference utilizes an extended Kalman filter.
3. The method of claim 1 , wherein the model of the target provides information of three-dimensional spatial distribution of the target at one or more points in time.
4. The method of claim 1 , wherein the model of the target is derived from one or more prior images.
5. The method of claim 4 , wherein the model of the target is further derived from a surrogate signal indicative of a motion of the target.
6. The method of claim 1 , wherein each present dataset of the target further includes a measure of an uncertainty in the respective position value of the target.
7. The method of claim 6 , wherein the measure of the uncertainty of the first present dataset is a predicted covariance matrix of a distribution of the target position.
8. The method of claim 7 , wherein the predicted covariance matrix is a function of an uncertainty of the state transition model.
9. The method of claim 8 , wherein the uncertainty of the state transition model is calculated using a target position tracked in a past time frame.
10. The method of claim 8 , wherein the present position value of the target in the first present dataset equals a mean position of the target within a same motion phase of the target in the past 10 seconds to 140 seconds.
11. The method of claim 10 further wherein the uncertainty of the state transition model is estimated to equal a covariance matrix of the target positions within the same motion phase in the past 10 seconds to 140 seconds.
12. The method of claim 1 , wherein measuring the second present dataset further comprises processing the two-dimensional image of the target prior to template-matching.
13. The method of claim 1 , wherein measuring the second present dataset further includes calculating a template-matching metric.
14. The method of claim 1 , wherein the second present dataset further includes a measure of an uncertainty in the two-dimensional present position value of the target.
15. The method of claim 14 , wherein the measure of the uncertainty in the two-dimensional present position value of the target is a function of the template-matching metric.
16. The method of claim 1 , wherein the target is a tumor.
17. The method of claim 1 , wherein the target is a lung tumor.
18. The method of claim 1 , wherein the model of the target accounts for the periodic nature of tumor motion.
19. The method of claim 1 , wherein predicting the first present dataset of the target further comprises accounting for the periodic nature of tumor motion.
20. The method of claim 1 , wherein the body is a human body.