WO2000024466A1 - Method and system for predictive physiological gating of radiation therapy - Google Patents

Method and system for predictive physiological gating of radiation therapy Download PDF

Info

Publication number
WO2000024466A1
WO2000024466A1 PCT/US1999/024946 US9924946W WO0024466A1 WO 2000024466 A1 WO2000024466 A1 WO 2000024466A1 US 9924946 W US9924946 W US 9924946W WO 0024466 A1 WO0024466 A1 WO 0024466A1
Authority
WO
WIPO (PCT)
Prior art keywords
period
physiological movement
data
physiological
patient
Prior art date
Application number
PCT/US1999/024946
Other languages
French (fr)
Inventor
Hassan Mostafavi
Original Assignee
Varian Medical Systems, Inc.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Varian Medical Systems, Inc. filed Critical Varian Medical Systems, Inc.
Priority to EP99970937A priority Critical patent/EP1123137B1/en
Priority to CA002348091A priority patent/CA2348091A1/en
Priority to JP2000578067A priority patent/JP4391023B2/en
Priority to DE69929628T priority patent/DE69929628T2/en
Priority to KR1020017005049A priority patent/KR20010083921A/en
Priority to AU12286/00A priority patent/AU771038B2/en
Publication of WO2000024466A1 publication Critical patent/WO2000024466A1/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/54Control of apparatus or devices for radiation diagnosis
    • A61B6/541Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • A61N5/1037Treatment planning systems taking into account the movement of the target, e.g. 4D-image based planning
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1064Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring

Definitions

  • the present invention relates to medical methods and systems. More particularly, the invention relates to a method and system for physiological gating of radiation therapy.
  • Radiation therapy involves medical procedures that selectively expose certain areas of a human body, such as cancerous tumors, to high doses of radiation.
  • the intent of the radiation therapy is to irradiate the targeted biological tissue such that the harmful tissue is destroyed.
  • the irradiation volume can be restricted to the size and shape of the tumor or targeted tissue region to avoid inflicting unnecessary radiation damage to healthy tissue.
  • conformal therapy is a radiotherapy technique that is often employed to optimize dose distribution by conforming the treatment volume more closely to the targeted tumor.
  • Normal physiological movement represents a limitation in the clinical planning and delivery of conventional radiotherapy and conformal therapy.
  • Normal physiological movement such as respiration or heart movement, can cause a positional movement of the tumor or tissue region undergoing irradiation. If the radiation beam has been shaped to conform the treatment volume to the exact dimensions of a tumor, then movement of that tumor during treatment could result in the radiation beam not being sufficiently sized or shaped to fully cover the targeted tumoral tissue.
  • One approach to this problem involves physiological gating of the radiation beam during treatment, with the gating signal synchronized to the movement of the patient's body.
  • instruments are utilized to measure the physiological state and/or movement of the patient. For example, respiration has been shown to cause movements in the position of a lung tumor in a patient's body.
  • a temperature sensor, strain gauge, preumotactrograph, or optical imaging system can be utilized to measure the patient's respiration cycle. These instruments can produce a signal indicative of the movement of the patient during the respiratory cycle.
  • the radiation beam can be gated based upon certain threshold amplitude levels of the measured respiratory signal, such that the radiation beam is disengaged or stopped during particular time points in the respiration signal that correspond to excessive movement of the lung tumor.
  • the gating system may include a switch for gating the radiation treatment, in which the switch requires a given time period ⁇ t to fully engage. If the switching time period ⁇ t is relatively slow compared to the measured physiological motion cycle, then a system employing such a switch in a reactive manner may not be able to effectively gate the application of radiation at appropriate time points during the radiation therapy.
  • the present invention provides an improved method and system for physiological gating for radiation therapy.
  • the invention comprises a method and system for detecting and predictably estimating regular cycles of physiological activity or movements.
  • Another aspect of the invention is directed to predictive actuation of gating system components.
  • Yet another aspect of the invention is directed to physiological gating of radiation treatment based upon the phase of the physiological activity.
  • FIG. 1 depicts the components of a system for physiological gating according to an embodiment of the invention.
  • Fig. 2 depicts an example of a respiratory motion signal chart.
  • Fig. 3 depicts a motion signal chart and a gating signal chart.
  • Fig. 4 is a flowchart showing process actions performed in an embodiment of the invention.
  • Fig. 5 is a flowchart showing process actions for detecting and predicting estimation of regular physiological movements.
  • Fig. 6a depicts a side view an embodiment of a camera that can be utilized in the invention.
  • Fig. 6b depicts a front view of the camera of Fig. 6a.
  • Fig. 7a depicts a retro-reflective marker according to an embodiment of the invention.
  • Fig. 7b depicts a cross-sectional view of the retro-reflective marker of Fig. 7a.
  • Fig. 8 depicts an apparatus for making a retro-reflective marker.
  • Fig. 9 depicts a phase chart synchronized with a gating signal chart.
  • Fig. 10 depicts an embodiment of a cylindrical marker block.
  • Fig. 11 depicts an embodiment of a hemispherical marker block.
  • Fig. 12 is a diagram of a computer hardware system with which the present invention can be implemented.
  • An aspect of the present invention comprises a method for detecting and predictively estimating regular cycles of physiological activity or movement.
  • the method of the invention can be employed for any regular physiological activity, including for example, the respiratory or cardiac cycles.
  • one or more sets of data representative of the physiological activity of interest are collected for the patient.
  • an electrocardiograph can be employed to generate data representative of the cardiac cycle.
  • a temperature sensor, strain gauge or preumotactrograph can be employed.
  • Other instruments or mechanisms can be employed within the scope of the invention to obtain sets of data representative of physiological activity or movements.
  • Fig. 1 depicts the components of an embodiment of a system 100 for physiological gating of radiation therapy, in which data representative of physiological activity is collected with an optical imaging apparatus.
  • System 100 comprises a radiation beam source 102 (such as a conventional linear accelerator) which is positionally configured to direct a radiation beam at a patient 106 located on treatment table 104.
  • a switch 116 is operatively coupled to the radiation beam source 102. Switch 116 can be operated to suspend the application of the radiation beam at patient 106.
  • switch 116 is part of the mechanical and electrical structure of radiation beam source 102.
  • switch 116 comprises an external apparatus that is connected to the control electronics of radiation beam source 102.
  • An optical or video image apparatus such as video camera 108, is aimed such that as least part of the patient 106 is within the camera's field of view.
  • Camera 108 monitors patient 106 for motion relating to the particular physiological activity being measured. For example, if respiration movements of the patient is being monitored, then camera 108 can be configured to monitor the motion of the patient's chest.
  • camera 108 is placed with its axis at approximately 45 degrees to the longitudinal axis of the patient 106.
  • the video image field of view is preferably set to view an approximately 20 cm by 20 cm area of the patient's chest.
  • a single camera 108 is shown in Fig. 1. However, the number of cameras 108 employed in the present invention can exceed that number, and the exact number to be used in the invention depends upon the particular application to which it is directed.
  • one or more illumination sources project light at the patient 106 on treatment table 104.
  • the generated light is reflected from one or more landmarks on the patient's body.
  • the camera 108 which is directed at patient 106, captures and detects the reflected light from the one or more landmarks.
  • the landmarks are selected based upon the physiological activity being studied. For example, for respiration measurements, landmarks are selected from one or more locations on the patient's chest.
  • the output signals of camera 108 are sent to a computer 110 or other type of processing unit having the capability to receive video images.
  • computer 110 comprises an Intel Pentium-based processor running Microsoft Windows NT and includes a video frame grabber card having a separate channel for each video source utilized in the system.
  • the images recorded by camera 108 are sent to computer 110 for processing. If camera 108 produces an analog output, the frame grabber converts the camera signals to a digital signal prior to processing by computer 110.
  • control signals can be sent from computer 110 to operate switch 116.
  • one or more passive markers 114 are located on the patient in the area to be detected for movement.
  • Each marker 114 preferably comprises a reflective or retro-reflective material that can reflect light, whether in the visible or invisible wavelengths. If the illumination source is co-located with camera 108, then marker 114 preferably comprises a retro-reflective material that reflects light mostly in the direction of the illumination source. Alternatively, each marker 114 comprises its own light source.
  • the marker 114 is used in place of or in conjunction with physical landmarks on the patient's body that is imaged by the camera 108 to detect patient movement. Markers 114 are preferably used instead of body landmarks because such markers 114 are easier to detect and track via the video image generated by camera 108. Because of the reflective or retro-reflective qualities of the preferred markers 114, the markers 114 inherently provide greater contrast in a video image to a light detecting apparatus such as camera 108, particularly when the camera 108 and illumination source are co-located.
  • a video or optical based system provides a reliable mechanism for repeating measurement results between uses on a given patient.
  • the method of the invention is noninvasive, and even if markers are used, no cables or connections must be made to the patient. Moreover, if the use of markers is impractical, the system can still be utilized without markers by performing measurements of physiological activity keyed to selected body landmarks.
  • the method of the invention is more accurate because it is based upon absolute measurement of external anatomical physical movement.
  • a possible inefficiency in tracking the markers 114 is that the marker may appear anywhere on the video frame, and all of the image elements of the video frame may have to be examined to determine the location of the marker 114.
  • the initial determination of locations for the marker 114 involves an examination of all of the image elements in the video frame. If the video frame comprise 640 by 480 image elements, then all 307200 (640*480) image elements are initially examined to find the location of the markers 114. For real-time tracking of the marker 114, examining every image element for every video frame to determine the location of the marker 114 in real-time could consume a significant amount of system resources. Thus, in an embodiment, the real-time tracking of marker 114 can be facilitated by processing a small region of the video frame, referred to herein as a "tracking gate", that is placed based on estimation of the location of the already-identified marker 114 in the video frame.
  • a small region of the video frame referred to herein as a "tracking gate"
  • the previously determined location of a marker 114 defined in the previous video frame is used to define an initial search range (i.e., the tracking gate) for that same marker in real-time.
  • the tracking gate is a relatively small portion of the video frame that is centered at the previous location of the marker 114.
  • the tracking gate is expanded only if it does not contain the new location of the marker 114.
  • the tracking gate is limited to a 50 by 50 area of the video frame, then the tracking gate for this example would comprise the image elements bound within the area defined by the coordinates (25,50), (75,50), (50,25), and (50,75).
  • the other portions of the video frame are searched only if the marker 106 is not found within this tracking gate.
  • Fig. 2 depicts an example of a motion signal chart 200 for respiratory movement that contains information regarding the movement of marker 114 during a given measurement period.
  • the horizontal axis represents points in time and the vertical axis represents the relative location or movement of the marker 114.
  • the illustrated signal in Fig. 2 comprises a plurality of discrete data points plotted along the motion signal chart 200.
  • threshold points can be defined over the amplitude range of the motion signal to determine the boundaries of the treatment intervals.
  • Motion of the patient that fall outside the boundaries of the treatment intervals correspond to movement that is predicted to cause unacceptable levels of movement to the tumor or tissue targeted for irradiation.
  • the treatment intervals correspond to the portion of the physiological cycle in which motion of the clinical target volume is minimized.
  • Other factors for determining the boundaries of the treatment intervals include identifying the portion of the motion signals involving the least movement of the target volume or the portion of the motion signal involving the largest separation of the target volume from organs at risk.
  • the radiation beam pattern can be shaped with the minimum possible margin to account for patient movement.
  • Fig. 3 depicted are examples of treatment intervals, indicated by signal range 302, that has been defined over the motion data shown in motion signal chart 200.
  • any movement of the measured body location that exceeds the value of 0.8 (shown by upper boundary line 304) or which moves below the value of 0.0 (shown by lower boundary line 306) falls outside the boundaries of the treatment intervals.
  • FIG. 3 Shown in Fig. 3 is an example of a gating signal chart 300 that is aligned with motion signal chart 200. Any motion signal that falls outside the treatment interval signal range 302 results in a "beam hold" gating signal threshold 310 that stops the application of radiation to the patient. Any motion signal that is within the treatment interval signal range 302 results in a "beam on" gating signal threshold 312 that allows radiation to be applied to the patient.
  • digital signals that represent the information shown in motion signal chart 200 are processed by computer 110 and compared to the threshold levels of the treatment interval signal range 302 to generate gating signal thresholds 310 and 312.
  • gating signal thresholds 310 and 312 can be obtained by feeding analog motion signals to a comparator to be compared with analog threshold signals that correspond to treatment interval signal range 302.
  • gating signal thresholds 310 and 312 are generated by computer 110 and are applied to the switch 116 that controls the operation of radiation beam source 102 (Fig. 1) to stop or start the application of a radiation beam at patient 106.
  • Fig. 4 is a flowchart of the process actions performed in an embodiment of the invention.
  • the first process action is to define boundaries for the treatment intervals over the range of motion signals to be detected by a camera (402). As indicated above, any motion that fall outside the boundaries of the treatment intervals correspond to motion that is predicted to result in unacceptable levels of movement of the tumor or tissue targeted for irradiation.
  • An optical or video imaging system such as a video camera, is used to measure the physiological motion of the patient (404), and the output signals of the optical or video imaging system are processed to compare the measured motion signals with the threshold boundaries of the treatment intervals (406).
  • a "beam off gating signal threshold is applied to a switch that is operatively coupled to the radiation beam source (408). If the radiation beam source is presently irradiating the patient (410), then the switch setting is operated to hold or stop the radiation beam (411). The process then returns back to process action 406.
  • a "beam on" gating signal threshold is produced (412) and is applied to a switch that is operatively coupled to the radiation beam source. If the radiation beam source is presently not being applied to the patient (413), then the switch setting is operated to turn on or apply the radiation beam source to irradiate the patient (414). The process then returns back to process action 406.
  • the radiation beam source can be disengaged if a significant deviation is detected in the regular physiological movements of the patient. Such deviations can result from sudden movement or coughing by the patient. The position and/or orientation of the targeted tissue may unacceptably shift as a result of this deviation, even though the amplitude range of the motion signal still falls within the boundaries of the treatment intervals during this deviation. Thus, detection of such deviations helps define the appropriate time periods to gate the radiation treatment.
  • the present invention provides a method for detecting and predictively estimating a period of a physiological activity. In effect, the present invention can "phase lock" to the physiological movement of the patient. Since the gating system phase locks to the movement period, deviations from that period can be identified and appropriately addressed.
  • the radiation treatment can be gated during these deviations from the regular period.
  • the present invention also provides a method for predictively estimating the period of the subsequent physiological movement to follow.
  • Fig. 5 is a process flowchart of an embodiment of the invention to perform predictive estimation and detection of regular physiological movement cycles.
  • an instrument or system (such as system 100 from Fig. 1) is employed to generate data signals representative of the physiological activity of interest.
  • the data signals comprises a stream of digital data samples that collectively form a signal wave pattern representative of the physiological movement under examination.
  • a number of discrete measurement samples are taken for the physiological activity during a given time period. For example, in an embodiment of the invention directed towards respiratory measurement, approximately 200-210 data samples are measured for each approximately 7 second time interval.
  • pattern matching analysis is performed against the measured data samples.
  • the most recent set of data samples for the physiological signal is correlated against an immediately preceding set of data samples to determine the period and repetitiveness of the signal.
  • An autocorrelation function can be employed to perform this pattern matching. For each new sample point of the physiological motion or physiological monitoring sensor signal, the process computes the autocorrelation function of the last n samples of the signal, where n corresponds to approximately 1.5 to 2 signal breathing periods. The secondary peak of the autocorrelation function is then identified to determine the period and repetitiveness of the signal.
  • an absolute difference function is used instead of an autocorrelation function. Instead of secondary peak, a secondary minimum in the absolute difference is searched for. For each new sample point of the physiological motion or physiological monitoring sensor signal, the process computes the minimum absolute difference between the two sets of data over a range of overlapping data samples. The secondary minimum corresponds to the data position that best matches the recent set of data samples with the preceding set of data samples. Yet another alternate embodiment performs a pattern matching based upon a model of the physiological activity being measured. The model is a dynamic representation of the physiological motion or physiological monitoring sensor signal for that physiological activity. The latest set of data samples is matched against the model to estimate parameters of the repetitive process.
  • Pattern matching using the measured physiological signal provides information regarding the degree of match, as well as a location of best match for the repetitive process. If an autocorrelation function is employed in process action 504, then the relative strength of secondary peak provides a measure of how repetitive the signal is.
  • a threshold range value is defined to provide indication of the degree of match between the two sets of data samples. If the strength of the secondary peak is within the defined threshold range (process action 508), then the degree of match indicates that the signal is repetitive, and the secondary peak location provides an estimate of the signal period.
  • an absolute difference function is used in process action 504, then the relative value of the secondary minimum provides a measure of how repetitive the signal is. If the value of the secondary minimum meets a defined threshold range (508), then the degree of match indicates that the signal is repetitive, and the secondary minimum location provides an estimate of the signal period.
  • the correlation value of the secondary peak or secondary minimum does not meet the defined threshold range, then a deviation from the regular physiological activity is detected, thereby indicating an irregularity in the regular physiological movement of the patient (510).
  • This irregularity could result, for example, from sudden movement or coughing of the patient.
  • this detected irregularity results in the generation of a "beam hold" signal that suspends the application of radiation at the patient.
  • the point of best match is tested to determine if the period is within a reasonable range.
  • the location of the secondary peak or secondary minimum provides an estimate of the period of the physiological activity.
  • the point of best match is compared to a threshold range (509). If the point of best match does not fall within the threshold range, than a deviation from regular physiological activity is detected (510). If the point of best match falls within the threshold range, then the signal is accepted as being repetitive (512).
  • process actions 504, 508, and 509 test for repetitiveness based upon a plurality of data samples over a range of such samples. However, in some circumstances, a significant deviation from normal physiological movement may actually occur within the new or most recent data sample(s) being analyzed, but because the overall set of data samples indicates repetitiveness (e.g., because of averaging of absolute differences over the range of data samples being compared), process actions 504, 508, and 509 may not detect the deviation. To perform a test for rapid deviation, the predicted value from process action 514 is compared with the next corresponding data sample (515).
  • the pattern matching process action (504) is performed over the entire range of data samples. Thereafter, the pattern matching process action can be performed over a limited search interval, which is defined by the results of the prior immediate execution of the process. For example, the predicted value from process action 514 can be used to define the location of the search interval for the next set of data samples. However, if process action 508, 509, and 514 detect deviation based upon analysis of the initial search interval, then the search interval can be expanded to ensure that a deviation has actually occurred. The process of Fig. 5 can be repeated with the increased search interval to attempt to find a point of best match outside of the initial search interval. In an embodiment, this increased search interval comprises the entire range of data samples. Alternatively, the increased search interval comprises only an expanded portion of the entire range of data samples.
  • physiological gating can be performed based upon the phase of the physiological activity, rather than its amplitude. This is in contrast to the example shown in Fig. 3, in which the amplitude of the physiological movement signal defines the boundaries of the treatment intervals for the application of radiation.
  • FIG. 9 depicted is an example of a chart 900 showing the phase progression of a physiological movement signal.
  • Treatment interval range 902 have been defined over the phase chart 900.
  • the boundaries of the treatment interval range 902 are defined by the phase of the detected signal. Radiation is applied to the patient only when the phase of the physiological movement signal falls within the boundaries of the treatment interval range 902.
  • Fig. 9 depicts examples of treatment interval range 902 having boundaries that span from 30 degrees to 300 degrees. Thus, the applied radiation to the patient is suspended or stopped when the phase of the physiological movement signal is between 301 degrees to 29 degrees.
  • Shown in Fig. 9 is a gating signal chart 906 that is aligned with phase chart 900.
  • “beam hold” signal threshold 910 results if the phase of the physiological movement signal falls outside the treatment interval range 902.
  • a “beam on” signal threshold 912 results if the phase of the physiological movement signal falls within the boundaries of the treatment interval range 902.
  • the "beam on” and “beam hold” signal thresholds 910 and 912 are applied to a switch 116 that operatively controls the operation of a radiation beam source 102 (Fig. 1). If radiation is being applied to the patient, application of the "beam hold” signal threshold 910 triggers switch 116 to suspend or stop the application of radiation. If radiation to the patient is not being applied, application of the "beam on” signal threshold 912 triggers the application of radiation to the patient.
  • the predictive qualities of the present invention permits operation of a gating system even if the measured physiological movement involves motion that that is relatively fast when compared to the effectively operating speeds of gating system components.
  • a gating system that includes a switch for gating the radiation treatment, in which the switch requires a known time period ⁇ t to fully engage. If the switching time period ⁇ t is relatively slow compared to the measured physiological motion cycle, then a system employing such a switch in a reactive manner may not be able to effectively gate the application of radiation at the patient.
  • the present invention allows predictive triggering of switch 116 to compensate for the amount of time ⁇ t required to fully engage the switch.
  • a predicted period for a physiological activity can be obtained by employing the process of Fig. 5.
  • a treatment interval range is defined over a portion of the period of the physiological activity. Based upon the time ⁇ t required to fully actuate the switch 116, the switch 116 can be pre-actuated by this time period ⁇ t prior to the time of the boundary of the treatment interval, so that the time for full actuation of the switch 116 coincides with the boundary of the treatment interval.
  • the radiation can be effectively gated at the boundaries of the treatment interval, regardless of the operating speeds of the switch 116.
  • the same procedure can also be employed to compensate for the operating speeds of other gating system components.
  • the following is an embodiment of the invention coded in the Visual Basic programming language. The following program code is directed to a process for detecting and predictively estimating the respiration cycle period using the absolute difference function:
  • StartJ TimeRefIdxBuf((i - 201 + BufLength) Mod BufLength)
  • CurrJ TimeRefldxBuf((i - 1 + BufLength) Mod BufLength)
  • Do j StartJ
  • variable "i” represents a counter or index to the data sample being processed.
  • the variable “Range” represents the search range that is to be analyzed. If the period of the physiological cycle has already been determined (i.e., from a prior execution of this program code), then the variable “Period” comprises the detected period. If the period has not yet been determined, then the variable "Period” is set to a default value representative of a normal respiration cycle (e.g., the number of data points in a normal breathing cycle, which is approximately 95 data samples in an embodiment of the invention where approximately 200-210 data samples are obtained over an approximately 7 second time period).
  • a default value representative of a normal respiration cycle e.g., the number of data points in a normal breathing cycle, which is approximately 95 data samples in an embodiment of the invention where approximately 200-210 data samples are obtained over an approximately 7 second time period.
  • MinAbsDiff is the minimum absolute difference value over the search range.
  • Diff represents a comparison between the actual value of a next data sample and the expected value of that next data sample.
  • the variables "j", “StartJ”, and “CurrJ” are counters or indexes into the data samples being processed.
  • the variable “k” is a counter to the search range.
  • the variable “MaxK” represents the position in the search range having the minimum absolute difference value.
  • the variable “AbsDiff maintains the sum of the absolute difference values for overlapping data samples.
  • the variable “NormaAbsDiff is the average absolute difference value for a particular position in the search range, which is represented as a percentage value.
  • the variable “n” is used to track the position of the data samples relative to the search range which is represented as a percentage value.
  • "Predict” is the predicted value that is returned by this program code.
  • the variable "MinAbsDiff is initialized to a high value so that so that any subsequent absolute difference value will be smaller than the initialized value.
  • the set of data samples being processed comprises 200 data points.
  • the variable "StartJ” is initialized back 201 data samples.
  • the variable “CurrJ” is initialized back one data sample. Because a circular array is being used, the "BufLength” variable is referenced during the initialization of both "StartJ" and "CurrJ".
  • the outer Do loop moves the current and preceding sets of data samples relative to each other. The outer Do loop is active while the variable "k" indicates that the program code is processing within the search range.
  • the search range is initially set at three positions to either side of a predicted position.
  • the predicted position is based upon the period obtained for an immediately prior execution of the program code. If the program code has not been executed immediately prior, then a default period value is used. If an acceptable minimum absolute difference value is not found within this search range, then the search range can be expanded to, for example, 50 positions to either side of the predicted position.
  • the variable "j" is initialized to the "StartJ” value.
  • the variables "AbsDiff and "n” are also initialized prior to execution of the inner Do loop.
  • the inner Do loop performs the computation of the absolute difference values between the present set and prior set of data samples.
  • the variable "AbsDiff maintains the sum of the absolute difference of values for overlapping data samples being compared. Note that the number of data samples being analyzed to determine the absolute difference values varies based upon the position in the search range being processed. This results because different positions in the search range has different numbers of data samples that overlap with the previous set of data samples being compared.
  • the absolute difference function is computed using every 10 th signal sample point, i.e., a subsampled subtraction is used.
  • the "Chartwidth” variable is referenced during the calculation of "AbsDiff .
  • the variables “MaxSignal” and “MinSignal” indicate a maximum and minimum range for signal values that have previously been sampled. These values can be established, for example, during a learning period for the system in which data samples are obtained for a plurality of respiratory cycles.
  • the "NormAbsDiff variable holds the average absolute difference value represented as a percentage value based upon the "MaxSignal” and "MinSignal” values.
  • the "NormAbsDiff value is less than or equal to a previously established "MinAbsDiff value, then the "MinAbsDiff variable is set to the "NormAbsDiff value.
  • the "MaxK” variable is set to the value of "k” if the "MinAbsDiff value is reset. The variable “k” is then incremented, and if the "k" value is still within the search range, then the program code returns back to the beginning of the outer Do loop.
  • the result of this program code is a candidate minimum absolute difference value and a candidate position for the minimum absolute difference.
  • the MaxK value is compared to predefined threshold values to ensure that it falls within a correct range of values for the physiological activity being processed. Thus, in an embodiment, the MaxK value is tested to make sure that it is greater than or equal to 40 and less than or equal to 150. If the mark value meets the threshold range, then the variable "Period” is set to the "MaxK” value. The variable “Predict” returns the predicted value for the next set of data samples to be processed.
  • the variable "Diff indicates the difference value between the current data sample value and the predicted data sample value, and is represented as a percentage to the "MaxSignal" and "MixSignal” values.
  • an image of a progress bar can be displayed to visually indicate the periodicity of the signal samples.
  • the visual progress bar is updated with the computed "MinAbsDiff value. Otherwise, the visual progress bar displays all other "MinAbsDiff values that exceed a 20% difference as a default value of "20".
  • Figs. 6a and 6b depict an embodiment of a camera 108 that can used in the present invention to optically or visually collect data representative of physiological movement.
  • Camera 108 is a charge-couple device ("CCD") camera having one or more photoelectric cathodes and one or more CCD devices.
  • a CCD device is a semiconductor device that can store charge in local areas, and upon appropriate control signals, transfers that charge to a readout point.
  • CCD device is a semiconductor device that can store charge in local areas, and upon appropriate control signals, transfers that charge to a readout point.
  • the electrons are captured in charge buckets located within the CCD device. The distribution of captured electrons in the charge buckets represents the image received at the camera.
  • camera 108 is a monochrome CCD camera having RS-170 output and 640x480 pixel resolution.
  • camera 408 can comprise a CCD camera having CCIR output and 756x567 pixel resolution.
  • an infra-red illuminator 602 (“IR illuminator”) is co-located with camera 108.
  • IR illuminator 602 produces one or more beams of infrared light that is directed in the same direction as camera 108.
  • IR illuminator 602 comprises a surface that is ringed around the lens 606 of camera body 608.
  • the surface of IR illuminator 602 contains a plurality of individual LED elements 604 for producing infrared light.
  • the LED elements 604 are arranged in a spiral pattern on the IR illuminator 602.
  • Infrared filters that may be part of the camera 108 are removed or disabled to increase the camera's sensitivity to infrared light.
  • Retro-reflective marker 700 comprises a raised reflective surface 702 for reflecting light.
  • Raised reflective surface 702 comprises a semi- spherical shape such that light can be reflected regardless of the input angle of the light source.
  • a flat surface 704 surrounds the raised reflective surface 702.
  • retro-reflective marker 700 is comprised of a retro-reflective material 3M#7610WS available from 3M Corporation.
  • marker 700 has a diameter of approximately .5 cm and a height of the highest point of raised reflective surface 702 of approximately .1 cm.
  • Fig. 8 depicts an apparatus 802 that can be employed to manufacture retro-reflective markers 700.
  • Apparatus 802 comprises a base portion 804 having an elastic ring 806 affixed thereto.
  • Elastic ring 806 is attached to bottom mold piece 808 having a bulge protruding from its center.
  • a control lever 810 can be operated to move top portion 812 along support rods 814.
  • Top portion 812 comprises a spring-loaded top mold piece 814.
  • Top mold piece 814 is formed with a semi-spherical cavity on its underside.
  • Control lever 810 is operated to move top portion 812 towards base portion 804.
  • the retro-reflective material is compressed and shaped between the bottom mold piece 808 and the top mold piece 814.
  • the top mold piece 814 forms the upper exterior of the retro-reflective material into a semi-spherical shape.
  • marker 114 comprises a marker block having one or more reference locations on its surface.
  • Each reference location on the marker block preferably comprises a retro-reflective or reflective material that is detectable by an optical imaging apparatus, such as camera 108.
  • Fig. 11 depicts an embodiment of a marker block 1100 having a cylindrical shape with multiple reference locations comprised of retro-reflective elements 1102 located on its surface.
  • Marker block 1100 can be formed as a rigid block (e.g., from Styrofoam). Blocks made in this fashion can be reused a plurality of times, even with multiple patients.
  • the retro-reflective elements 1102 can be formed from the same material used to construct retro-reflective markers 114 of Figs. 7a and 7b.
  • the marker block is preferably formed from a material that is lightweight enough not to interfere with normal breathing by the patient.
  • a marker block can be formed into any shape or size, as long as the size, spacing, and positioning of the reference locations are configured such that a camera or other optical imaging apparatus can view and generate an image that accurately shows the positioning of the marker block.
  • Fig. 10 depicts an alternate marker block 1000 having a hemispherical shape comprised of a plurality of retro-reflective elements 1002 attached to its surface.
  • the marker block can be formed with shapes to fit particular body parts. For example, molds or casts that match to specific locations on the body can be employed as marker blocks. Marker blocks shaped to fit certain areas of the body facilitate the repeatable placement of the marker blocks at particular locations on the patient.
  • the marker blocks can be formed to fit certain fixtures that are attached to a patient's body.
  • a marker block can be formed within indentations and grooves that allow it to be attached to eyeglasses.
  • the fixtures are formed with integral marker block(s) having reflective or retro-reflective markers on them.
  • An alternate embodiment of the marker block comprises only a single reference location/reflective element on its surface.
  • This embodiment of the marker block is used in place of the retro-reflective marker 406 to detect particular locations on a patient's body with an optical imaging apparatus. Further details regarding camera 108, markers 114, marker blocks, or procedures to optically measure physiological movements are described in copending U.S. Patent Application Ser. No. [Not Yet Assigned] (Varian Docket No. 98-34, Attorney Docket No. 236/223) and U.S. Patent Application Ser. No. [Not Yet Assigned] (Varian Docket No. 98-33, Attorney Docket No. 236/224), filed concurrently herewith, both of which are hereby incorporated by reference in their entirety.
  • FIG 12 is a block diagram that illustrates an embodiment of a computer system 1900 upon which an embodiment of the invention may be implemented.
  • Computer system 1900 includes a bus 1902 or other communication mechanism for communicating information, and a processor 1904 coupled with bus 1902 for processing information.
  • Computer system 1900 also includes a main memory 1906, such as a random access memory (RAM) or other dynamic storage device, coupled to bus 1902 for storing information and instructions to be executed by processor 1904.
  • Main memory 1906 also may be used for storing temporary variables or other intermediate information during execution of instructions to be executed by processor 1904.
  • Computer system 1900 further includes a read only memory (ROM) 1908 or other static storage device coupled to bus 1902 for storing static information and instructions for processor 1904.
  • ROM read only memory
  • a data storage device 1910 such as a magnetic disk or optical disk, is provided and coupled to bus 1902 for storing information and instructions.
  • Computer system 1900 may be coupled via bus 1902 to a display 1912, such as a cathode ray tube (CRT), for displaying information to a user.
  • a display 1912 such as a cathode ray tube (CRT)
  • An input device 1914 is coupled to bus 1902 for communicating information and command selections to processor 1904.
  • cursor control 1916 is Another type of user input device
  • cursor control 1916 such as a mouse, a trackball, or cursor direction keys for communicating direction information and command selections to processor 1904 and for controlling cursor movement on display 1912.
  • This input device typically has two degrees of freedom in two axes, a first axis (e.g., x) and a second axis (e.g., y), that allows the device to specify positions in a plane.
  • the invention is related to the use of computer system 1900 for detecting and predictively estimating physiological cycles. According to one embodiment of the invention, such use is provided by computer system 1900 in response to processor 1904 executing one or more sequences of one or more instructions contained in main memory 1906. Such instructions may be read into main memory 1906 from another computer-readable medium, such as storage device 1910. Execution of the sequences of instructions contained in main memory 1906 causes processor 1904 to perform the process steps described herein. One or more processors in a multi-processing arrangement may also be employed to execute the sequences of instructions contained in main memory 1906. In alternative embodiments, hard- wired circuitry may be used in place of or in combination with software instructions to implement the invention. Thus, embodiments of the invention are not limited to any specific combination of hardware circuitry and software.
  • Non- volatile media includes, for example, optical or magnetic disks, such as storage device 1910.
  • Volatile media includes dynamic memory, such as main memory 1906.
  • Transmission media includes coaxial cables, copper wire and fiber optics, including the wires that comprise bus 1902. Transmission media can also take the form of acoustic or light waves, such as those generated during radio wave and infrared data communications.
  • Computer-readable media include, for example, a floppy disk, a flexible disk, hard disk, magnetic tape, or any other magnetic medium, a CD-ROM, any other optical medium, punch cards, paper tape, any other physical medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, a carrier wave as described hereinafter, or any other medium from which a computer can read.
  • Various forms of computer readable media may be involved in carrying one or more sequences of one or more instructions to processor 1904 for execution.
  • the instructions may initially be carried on a magnetic disk of a remote computer.
  • the remote computer can load the instructions into its dynamic memory and send the instructions over a telephone line using a modem.
  • a modem local to computer system 1900 can receive the data on the telephone line and use an infrared transmitter to convert the data to an infrared signal.
  • An infrared detector coupled to bus 1902 can receive the data carried in the infrared signal and place the data on bus 1902.
  • Bus 1902 carries the data to main memory 1906, from which processor 1904 retrieves and executes the instructions.
  • the instructions received by main memory 1906 may optionally be stored on storage device 1910 either before or after execution by processor 1904.
  • Computer system 1900 also includes a communication interface 1918 coupled to bus 1902.
  • Communication interface 1918 provides a two-way data communication coupling to a network link 1920 that is connected to a local network 1922.
  • communication interface 1918 may be an integrated services digital network (ISDN) card or a modem to provide a data communication connection to a corresponding type of telephone line.
  • ISDN integrated services digital network
  • communication interface 1918 may be a local area network (LAN) card to provide a data communication connection to a compatible LAN.
  • LAN local area network
  • Wireless links may also be implemented.
  • communication interface 1918 sends and receives electrical, electromagnetic or optical signals that carry data streams representing various types of information.
  • Network link 1920 typically provides data communication through one or more networks to other devices.
  • network link 1920 may provide a connection through local network 1922 to a host computer 1924 or to medical equipment 1926 such as a radiation beam source or a switch operatively coupled to a radiation beam source.
  • the data streams transported over network link 1920 can comprise electrical, electromagnetic or optical signals.
  • the signals through the various networks and the signals on network link 1920 and through communication interface 1918, which carry data to and from computer system 1900, are exemplary forms of carrier waves transporting the information.
  • Computer system 1900 can send messages and receive data, including program code, through the network(s), network link 1920 and communication interface 1918.
  • computer 110 can be performed by any combination of hardware and software within the scope of the invention, and should not be limited to particular embodiments comprising just a particular definition of "computer".
  • the specification and drawings are, accordingly, to be regarded in an illustrative rather than restrictive sense.

Abstract

A method and system for physiological gating for radiation therapy is disclosed. A method and system for detecting and predictably estimating regular cycles of physiological activity or movements is disclosed. Another disclosed aspect of the invention is directed to predictive actuation of gating system components. Yet another disclosed aspect of the invention is directed to physiological gating of radiation treatment based upon the phase of the physiological activity.

Description

METHOD AND SYSTEM FOR PREDICTIVE PHYSIOLOGICAL GATING OF
RADIATION THERAPY
BACKGROUND OF THE INVENTION Field Of The Invention
The present invention relates to medical methods and systems. More particularly, the invention relates to a method and system for physiological gating of radiation therapy.
Background Radiation therapy involves medical procedures that selectively expose certain areas of a human body, such as cancerous tumors, to high doses of radiation. The intent of the radiation therapy is to irradiate the targeted biological tissue such that the harmful tissue is destroyed. In certain types of radiotherapy, the irradiation volume can be restricted to the size and shape of the tumor or targeted tissue region to avoid inflicting unnecessary radiation damage to healthy tissue. For example, conformal therapy is a radiotherapy technique that is often employed to optimize dose distribution by conforming the treatment volume more closely to the targeted tumor.
Normal physiological movement represents a limitation in the clinical planning and delivery of conventional radiotherapy and conformal therapy. Normal physiological movement, such as respiration or heart movement, can cause a positional movement of the tumor or tissue region undergoing irradiation. If the radiation beam has been shaped to conform the treatment volume to the exact dimensions of a tumor, then movement of that tumor during treatment could result in the radiation beam not being sufficiently sized or shaped to fully cover the targeted tumoral tissue. One approach to this problem involves physiological gating of the radiation beam during treatment, with the gating signal synchronized to the movement of the patient's body. In this approach, instruments are utilized to measure the physiological state and/or movement of the patient. For example, respiration has been shown to cause movements in the position of a lung tumor in a patient's body. If radiotherapy is being applied to the lung tumor, then a temperature sensor, strain gauge, preumotactrograph, or optical imaging system can be utilized to measure the patient's respiration cycle. These instruments can produce a signal indicative of the movement of the patient during the respiratory cycle. The radiation beam can be gated based upon certain threshold amplitude levels of the measured respiratory signal, such that the radiation beam is disengaged or stopped during particular time points in the respiration signal that correspond to excessive movement of the lung tumor.
Known approaches to physiological gating of radiation therapy are reactive, that is, known approaches utilize gating methods that slavishly react to measured levels of physiological movements. One drawback to reactive gating systems is that the measured physiological movement may involve motion that that is relatively fast when compared to the effectively operating speeds of gating system components. Thus, a purely reactive gating system may not be able to react fast enough to effectively gate the applied radiation. For example, the gating system may include a switch for gating the radiation treatment, in which the switch requires a given time period Δt to fully engage. If the switching time period Δt is relatively slow compared to the measured physiological motion cycle, then a system employing such a switch in a reactive manner may not be able to effectively gate the application of radiation at appropriate time points during the radiation therapy.
Therefore, there is a need for a system and method to address these and other problems of the related art. There is a need for a method and system for physiological gating which is not purely reactive to measure physiological movement signals.
SUMMARY OF THE INVENTION
The present invention provides an improved method and system for physiological gating for radiation therapy. According to an aspect, the invention comprises a method and system for detecting and predictably estimating regular cycles of physiological activity or movements. Another aspect of the invention is directed to predictive actuation of gating system components. Yet another aspect of the invention is directed to physiological gating of radiation treatment based upon the phase of the physiological activity.
These and other aspects, objects, and advantages of the invention are described below in the detailed description, drawings, and claims.
BRIEF DESCRIPTION OF THE DRAWINGS
The accompanying drawings are included to provide a further understanding of the invention and, together with the Detailed Description, serve to explain the principles of the invention. Fig. 1 depicts the components of a system for physiological gating according to an embodiment of the invention.
Fig. 2 depicts an example of a respiratory motion signal chart.
Fig. 3 depicts a motion signal chart and a gating signal chart.
Fig. 4 is a flowchart showing process actions performed in an embodiment of the invention.
Fig. 5 is a flowchart showing process actions for detecting and predicting estimation of regular physiological movements.
Fig. 6a depicts a side view an embodiment of a camera that can be utilized in the invention. Fig. 6b depicts a front view of the camera of Fig. 6a.
Fig. 7a depicts a retro-reflective marker according to an embodiment of the invention.
Fig. 7b depicts a cross-sectional view of the retro-reflective marker of Fig. 7a.
Fig. 8 depicts an apparatus for making a retro-reflective marker.
Fig. 9 depicts a phase chart synchronized with a gating signal chart. Fig. 10 depicts an embodiment of a cylindrical marker block.
Fig. 11 depicts an embodiment of a hemispherical marker block.
Fig. 12 is a diagram of a computer hardware system with which the present invention can be implemented. DETAILED DESCRIPTION
An aspect of the present invention comprises a method for detecting and predictively estimating regular cycles of physiological activity or movement. The method of the invention can be employed for any regular physiological activity, including for example, the respiratory or cardiac cycles.
In operation, one or more sets of data representative of the physiological activity of interest are collected for the patient. For example, an electrocardiograph can be employed to generate data representative of the cardiac cycle. To generate data representative of the respiration cycle, a temperature sensor, strain gauge or preumotactrograph can be employed. Other instruments or mechanisms can be employed within the scope of the invention to obtain sets of data representative of physiological activity or movements.
Fig. 1 depicts the components of an embodiment of a system 100 for physiological gating of radiation therapy, in which data representative of physiological activity is collected with an optical imaging apparatus. System 100 comprises a radiation beam source 102 (such as a conventional linear accelerator) which is positionally configured to direct a radiation beam at a patient 106 located on treatment table 104. A switch 116 is operatively coupled to the radiation beam source 102. Switch 116 can be operated to suspend the application of the radiation beam at patient 106. In an embodiment, switch 116 is part of the mechanical and electrical structure of radiation beam source 102. Alternatively, switch 116 comprises an external apparatus that is connected to the control electronics of radiation beam source 102.
An optical or video image apparatus, such as video camera 108, is aimed such that as least part of the patient 106 is within the camera's field of view. Camera 108 monitors patient 106 for motion relating to the particular physiological activity being measured. For example, if respiration movements of the patient is being monitored, then camera 108 can be configured to monitor the motion of the patient's chest. According to an embodiment, camera 108 is placed with its axis at approximately 45 degrees to the longitudinal axis of the patient 106. For measurement of respiration activity that could result in about 3-5 mm of chest motion, the video image field of view is preferably set to view an approximately 20 cm by 20 cm area of the patient's chest. For purposes of illustration only, a single camera 108 is shown in Fig. 1. However, the number of cameras 108 employed in the present invention can exceed that number, and the exact number to be used in the invention depends upon the particular application to which it is directed.
In an embodiment, one or more illumination sources (which are infrared sources in the preferred embodiment) project light at the patient 106 on treatment table 104. The generated light is reflected from one or more landmarks on the patient's body. The camera 108, which is directed at patient 106, captures and detects the reflected light from the one or more landmarks. The landmarks are selected based upon the physiological activity being studied. For example, for respiration measurements, landmarks are selected from one or more locations on the patient's chest.
The output signals of camera 108 are sent to a computer 110 or other type of processing unit having the capability to receive video images. According to a particular embodiment, computer 110 comprises an Intel Pentium-based processor running Microsoft Windows NT and includes a video frame grabber card having a separate channel for each video source utilized in the system. The images recorded by camera 108 are sent to computer 110 for processing. If camera 108 produces an analog output, the frame grabber converts the camera signals to a digital signal prior to processing by computer 110. Based upon the video signals received by computer 110, control signals can be sent from computer 110 to operate switch 116.
According to one embodiment, one or more passive markers 114 are located on the patient in the area to be detected for movement. Each marker 114 preferably comprises a reflective or retro-reflective material that can reflect light, whether in the visible or invisible wavelengths. If the illumination source is co-located with camera 108, then marker 114 preferably comprises a retro-reflective material that reflects light mostly in the direction of the illumination source. Alternatively, each marker 114 comprises its own light source. The marker 114 is used in place of or in conjunction with physical landmarks on the patient's body that is imaged by the camera 108 to detect patient movement. Markers 114 are preferably used instead of body landmarks because such markers 114 are easier to detect and track via the video image generated by camera 108. Because of the reflective or retro-reflective qualities of the preferred markers 114, the markers 114 inherently provide greater contrast in a video image to a light detecting apparatus such as camera 108, particularly when the camera 108 and illumination source are co-located.
Utilizing a video or optical based system to track patient movement provides several advantages. First, a video or optical based system provides a reliable mechanism for repeating measurement results between uses on a given patient. Second, the method of the invention is noninvasive, and even if markers are used, no cables or connections must be made to the patient. Moreover, if the use of markers is impractical, the system can still be utilized without markers by performing measurements of physiological activity keyed to selected body landmarks. Finally, the method of the invention is more accurate because it is based upon absolute measurement of external anatomical physical movement. A possible inefficiency in tracking the markers 114 is that the marker may appear anywhere on the video frame, and all of the image elements of the video frame may have to be examined to determine the location of the marker 114. Thus, in an embodiment, the initial determination of locations for the marker 114 involves an examination of all of the image elements in the video frame. If the video frame comprise 640 by 480 image elements, then all 307200 (640*480) image elements are initially examined to find the location of the markers 114. For real-time tracking of the marker 114, examining every image element for every video frame to determine the location of the marker 114 in real-time could consume a significant amount of system resources. Thus, in an embodiment, the real-time tracking of marker 114 can be facilitated by processing a small region of the video frame, referred to herein as a "tracking gate", that is placed based on estimation of the location of the already-identified marker 114 in the video frame. The previously determined location of a marker 114 defined in the previous video frame is used to define an initial search range (i.e., the tracking gate) for that same marker in real-time. The tracking gate is a relatively small portion of the video frame that is centered at the previous location of the marker 114. The tracking gate is expanded only if it does not contain the new location of the marker 114. As an example, consider the situation when the previously determined location of a particular marker is image element (50,50) in a video frame. If the tracking gate is limited to a 50 by 50 area of the video frame, then the tracking gate for this example would comprise the image elements bound within the area defined by the coordinates (25,50), (75,50), (50,25), and (50,75). The other portions of the video frame are searched only if the marker 106 is not found within this tracking gate.
The video image signals sent from camera 108 to computer 110 are used to generate and track motion signals representative of the movement of marker 114 and/or landmark structures on the patient's body. Fig. 2 depicts an example of a motion signal chart 200 for respiratory movement that contains information regarding the movement of marker 114 during a given measurement period. The horizontal axis represents points in time and the vertical axis represents the relative location or movement of the marker 114. According to an embodiment, the illustrated signal in Fig. 2 comprises a plurality of discrete data points plotted along the motion signal chart 200. An important aspect of physiological gating of radiotherapy is the determination of the boundaries of the "treatment intervals" for applying radiation. For gating purposes, threshold points can be defined over the amplitude range of the motion signal to determine the boundaries of the treatment intervals. Motion of the patient that fall outside the boundaries of the treatment intervals correspond to movement that is predicted to cause unacceptable levels of movement to the tumor or tissue targeted for irradiation. According to an embodiment, the treatment intervals correspond to the portion of the physiological cycle in which motion of the clinical target volume is minimized. Other factors for determining the boundaries of the treatment intervals include identifying the portion of the motion signals involving the least movement of the target volume or the portion of the motion signal involving the largest separation of the target volume from organs at risk. Thus, the radiation beam pattern can be shaped with the minimum possible margin to account for patient movement.
Radiation is applied to the patient only when the motion signal is within the designated treatment intervals. Referring to Fig. 3, depicted are examples of treatment intervals, indicated by signal range 302, that has been defined over the motion data shown in motion signal chart 200. In the example of Fig. 3, any movement of the measured body location that exceeds the value of 0.8 (shown by upper boundary line 304) or which moves below the value of 0.0 (shown by lower boundary line 306) falls outside the boundaries of the treatment intervals.
Shown in Fig. 3 is an example of a gating signal chart 300 that is aligned with motion signal chart 200. Any motion signal that falls outside the treatment interval signal range 302 results in a "beam hold" gating signal threshold 310 that stops the application of radiation to the patient. Any motion signal that is within the treatment interval signal range 302 results in a "beam on" gating signal threshold 312 that allows radiation to be applied to the patient. In an embodiment, digital signals that represent the information shown in motion signal chart 200 are processed by computer 110 and compared to the threshold levels of the treatment interval signal range 302 to generate gating signal thresholds 310 and 312. Alternatively, gating signal thresholds 310 and 312 can be obtained by feeding analog motion signals to a comparator to be compared with analog threshold signals that correspond to treatment interval signal range 302. In any case, gating signal thresholds 310 and 312 are generated by computer 110 and are applied to the switch 116 that controls the operation of radiation beam source 102 (Fig. 1) to stop or start the application of a radiation beam at patient 106.
Fig. 4 is a flowchart of the process actions performed in an embodiment of the invention. The first process action is to define boundaries for the treatment intervals over the range of motion signals to be detected by a camera (402). As indicated above, any motion that fall outside the boundaries of the treatment intervals correspond to motion that is predicted to result in unacceptable levels of movement of the tumor or tissue targeted for irradiation. An optical or video imaging system, such as a video camera, is used to measure the physiological motion of the patient (404), and the output signals of the optical or video imaging system are processed to compare the measured motion signals with the threshold boundaries of the treatment intervals (406).
If the motion signal is outside the boundaries of the treatment intervals, then a "beam off gating signal threshold is applied to a switch that is operatively coupled to the radiation beam source (408). If the radiation beam source is presently irradiating the patient (410), then the switch setting is operated to hold or stop the radiation beam (411). The process then returns back to process action 406.
If the motion signal is within the boundaries of the treatment intervals, then a "beam on" gating signal threshold is produced (412) and is applied to a switch that is operatively coupled to the radiation beam source. If the radiation beam source is presently not being applied to the patient (413), then the switch setting is operated to turn on or apply the radiation beam source to irradiate the patient (414). The process then returns back to process action 406.
According to one embodiment, the radiation beam source can be disengaged if a significant deviation is detected in the regular physiological movements of the patient. Such deviations can result from sudden movement or coughing by the patient. The position and/or orientation of the targeted tissue may unacceptably shift as a result of this deviation, even though the amplitude range of the motion signal still falls within the boundaries of the treatment intervals during this deviation. Thus, detection of such deviations helps define the appropriate time periods to gate the radiation treatment. The present invention provides a method for detecting and predictively estimating a period of a physiological activity. In effect, the present invention can "phase lock" to the physiological movement of the patient. Since the gating system phase locks to the movement period, deviations from that period can be identified and appropriately addressed. For example, when gating to the respiratory cycle, sudden movement or coughing by the patient can result in deviation from the detected period of the respiration cycle. The radiation treatment can be gated during these deviations from the regular period. The present invention also provides a method for predictively estimating the period of the subsequent physiological movement to follow.
Fig. 5 is a process flowchart of an embodiment of the invention to perform predictive estimation and detection of regular physiological movement cycles. In process action 502, an instrument or system (such as system 100 from Fig. 1) is employed to generate data signals representative of the physiological activity of interest. In an embodiment, the data signals comprises a stream of digital data samples that collectively form a signal wave pattern representative of the physiological movement under examination. A number of discrete measurement samples are taken for the physiological activity during a given time period. For example, in an embodiment of the invention directed towards respiratory measurement, approximately 200-210 data samples are measured for each approximately 7 second time interval.
In process action 504, pattern matching analysis is performed against the measured data samples. In an embodiment, the most recent set of data samples for the physiological signal is correlated against an immediately preceding set of data samples to determine the period and repetitiveness of the signal. An autocorrelation function can be employed to perform this pattern matching. For each new sample point of the physiological motion or physiological monitoring sensor signal, the process computes the autocorrelation function of the last n samples of the signal, where n corresponds to approximately 1.5 to 2 signal breathing periods. The secondary peak of the autocorrelation function is then identified to determine the period and repetitiveness of the signal.
In an alternate embodiment, an absolute difference function is used instead of an autocorrelation function. Instead of secondary peak, a secondary minimum in the absolute difference is searched for. For each new sample point of the physiological motion or physiological monitoring sensor signal, the process computes the minimum absolute difference between the two sets of data over a range of overlapping data samples. The secondary minimum corresponds to the data position that best matches the recent set of data samples with the preceding set of data samples. Yet another alternate embodiment performs a pattern matching based upon a model of the physiological activity being measured. The model is a dynamic representation of the physiological motion or physiological monitoring sensor signal for that physiological activity. The latest set of data samples is matched against the model to estimate parameters of the repetitive process. Pattern matching using the measured physiological signal (504) provides information regarding the degree of match, as well as a location of best match for the repetitive process. If an autocorrelation function is employed in process action 504, then the relative strength of secondary peak provides a measure of how repetitive the signal is. A threshold range value is defined to provide indication of the degree of match between the two sets of data samples. If the strength of the secondary peak is within the defined threshold range (process action 508), then the degree of match indicates that the signal is repetitive, and the secondary peak location provides an estimate of the signal period. If an absolute difference function is used in process action 504, then the relative value of the secondary minimum provides a measure of how repetitive the signal is. If the value of the secondary minimum meets a defined threshold range (508), then the degree of match indicates that the signal is repetitive, and the secondary minimum location provides an estimate of the signal period.
If the correlation value of the secondary peak or secondary minimum does not meet the defined threshold range, then a deviation from the regular physiological activity is detected, thereby indicating an irregularity in the regular physiological movement of the patient (510). This irregularity could result, for example, from sudden movement or coughing of the patient. In an embodiment, this detected irregularity results in the generation of a "beam hold" signal that suspends the application of radiation at the patient.
If the degree of match indicates repetitiveness, the point of best match is tested to determine if the period is within a reasonable range. The location of the secondary peak or secondary minimum provides an estimate of the period of the physiological activity. In an embodiment, the point of best match is compared to a threshold range (509). If the point of best match does not fall within the threshold range, than a deviation from regular physiological activity is detected (510). If the point of best match falls within the threshold range, then the signal is accepted as being repetitive (512).
The estimate of the period based on the point of best match can be used to predict the period and waveform parameters of the next set of data samples for the signal (514). Note that process actions 504, 508, and 509 test for repetitiveness based upon a plurality of data samples over a range of such samples. However, in some circumstances, a significant deviation from normal physiological movement may actually occur within the new or most recent data sample(s) being analyzed, but because the overall set of data samples indicates repetitiveness (e.g., because of averaging of absolute differences over the range of data samples being compared), process actions 504, 508, and 509 may not detect the deviation. To perform a test for rapid deviation, the predicted value from process action 514 is compared with the next corresponding data sample (515). If the predicted value does not match the actual data sample value within a defined threshold range, then deviation is detected (510). If a comparison of the predicted and actual data sample values fall within the defined threshold range, then repetitiveness is confirmed, and deviation is not detected for that data sample range (516).
In an embodiment, the first time the process of Fig. 5 is performed, the pattern matching process action (504) is performed over the entire range of data samples. Thereafter, the pattern matching process action can be performed over a limited search interval, which is defined by the results of the prior immediate execution of the process. For example, the predicted value from process action 514 can be used to define the location of the search interval for the next set of data samples. However, if process action 508, 509, and 514 detect deviation based upon analysis of the initial search interval, then the search interval can be expanded to ensure that a deviation has actually occurred. The process of Fig. 5 can be repeated with the increased search interval to attempt to find a point of best match outside of the initial search interval. In an embodiment, this increased search interval comprises the entire range of data samples. Alternatively, the increased search interval comprises only an expanded portion of the entire range of data samples.
According to an embodiment of the invention, physiological gating can be performed based upon the phase of the physiological activity, rather than its amplitude. This is in contrast to the example shown in Fig. 3, in which the amplitude of the physiological movement signal defines the boundaries of the treatment intervals for the application of radiation.
Referring to Fig. 9, depicted is an example of a chart 900 showing the phase progression of a physiological movement signal. Treatment interval range 902 have been defined over the phase chart 900. In the example of Fig. 9, the boundaries of the treatment interval range 902 are defined by the phase of the detected signal. Radiation is applied to the patient only when the phase of the physiological movement signal falls within the boundaries of the treatment interval range 902. Fig. 9 depicts examples of treatment interval range 902 having boundaries that span from 30 degrees to 300 degrees. Thus, the applied radiation to the patient is suspended or stopped when the phase of the physiological movement signal is between 301 degrees to 29 degrees. Shown in Fig. 9 is a gating signal chart 906 that is aligned with phase chart 900. A
"beam hold" signal threshold 910 results if the phase of the physiological movement signal falls outside the treatment interval range 902. A "beam on" signal threshold 912 results if the phase of the physiological movement signal falls within the boundaries of the treatment interval range 902. The "beam on" and "beam hold" signal thresholds 910 and 912 are applied to a switch 116 that operatively controls the operation of a radiation beam source 102 (Fig. 1). If radiation is being applied to the patient, application of the "beam hold" signal threshold 910 triggers switch 116 to suspend or stop the application of radiation. If radiation to the patient is not being applied, application of the "beam on" signal threshold 912 triggers the application of radiation to the patient. The predictive qualities of the present invention permits operation of a gating system even if the measured physiological movement involves motion that that is relatively fast when compared to the effectively operating speeds of gating system components. As just one example, consider a gating system that includes a switch for gating the radiation treatment, in which the switch requires a known time period Δt to fully engage. If the switching time period Δt is relatively slow compared to the measured physiological motion cycle, then a system employing such a switch in a reactive manner may not be able to effectively gate the application of radiation at the patient. The present invention allows predictive triggering of switch 116 to compensate for the amount of time Δt required to fully engage the switch. A predicted period for a physiological activity can be obtained by employing the process of Fig. 5. A treatment interval range is defined over a portion of the period of the physiological activity. Based upon the time Δt required to fully actuate the switch 116, the switch 116 can be pre-actuated by this time period Δt prior to the time of the boundary of the treatment interval, so that the time for full actuation of the switch 116 coincides with the boundary of the treatment interval. Thus, the radiation can be effectively gated at the boundaries of the treatment interval, regardless of the operating speeds of the switch 116. The same procedure can also be employed to compensate for the operating speeds of other gating system components. The following is an embodiment of the invention coded in the Visual Basic programming language. The following program code is directed to a process for detecting and predictively estimating the respiration cycle period using the absolute difference function:
Public Function Predict(ByVal i As Long, ByVal Range As Long, Period As Double, MinAbsDiff As Double, Diff As Double) As Double
Dim j As Long, Start! As Long, CurrJ As Long Dim k As Long, MaxK As Long Dim AbsDiff As Double Dim NormAbsDiff As Double, n As Long k = Period - Range MinAbsDiff = 10000000#
StartJ = TimeRefIdxBuf((i - 201 + BufLength) Mod BufLength) CurrJ = TimeRefldxBuf((i - 1 + BufLength) Mod BufLength) Do j = StartJ
AbsDiff = 0# n = 0 Do
AbsDiff = AbsDiff + Abs(SigBuf(SigRefIdxBufG)) - SigBuf(SigRefIdxBuf(G + k + ChartWidth) Mod ChartWidth))) n = n + 1 j = (j + 10) Mod ChartWidth
Loop While n <= (200 - k) / 10
NormAbsDiff = 100 * AbsDiff/ (n * Abs(MaxSignal - MinSignal)) If NormAbsDiff <= MinAbsDiff Then MinAbsDiff = NormAbsDiff MaxK = k
End If k = k + l Loop While k <= Period + Range If MaxK >= 40 And MaxK <= 150 Then Period = MaxK Predict = SigBuf(SigRefldxBuf((CurrJ - Period + ChartWidth) Mod ChartWidth))
Diff = 100 * Abs(SigBuf(SigRefIdxBuf(CurrJ)) - Predict) / Abs(MaxSignal - MinSignal) If MinAbsDiff <= 20 Then
ProgressBarl .Value = MinAbsDiff Else ProgressBarl .Value = 20
End If End Function
In this program code, the variable "i" represents a counter or index to the data sample being processed. The variable "Range" represents the search range that is to be analyzed. If the period of the physiological cycle has already been determined (i.e., from a prior execution of this program code), then the variable "Period" comprises the detected period. If the period has not yet been determined, then the variable "Period" is set to a default value representative of a normal respiration cycle (e.g., the number of data points in a normal breathing cycle, which is approximately 95 data samples in an embodiment of the invention where approximately 200-210 data samples are obtained over an approximately 7 second time period). The variable
"MinAbsDiff is the minimum absolute difference value over the search range. The variable
"Diff represents a comparison between the actual value of a next data sample and the expected value of that next data sample. The variables "j", "StartJ", and "CurrJ" are counters or indexes into the data samples being processed. The variable "k" is a counter to the search range. The variable "MaxK" represents the position in the search range having the minimum absolute difference value. The variable "AbsDiff maintains the sum of the absolute difference values for overlapping data samples. The variable "NormaAbsDiff is the average absolute difference value for a particular position in the search range, which is represented as a percentage value. The variable "n" is used to track the position of the data samples relative to the search range which is represented as a percentage value. "Predict" is the predicted value that is returned by this program code.
The variable "MinAbsDiff is initialized to a high value so that so that any subsequent absolute difference value will be smaller than the initialized value. In an embodiment, the set of data samples being processed comprises 200 data points. Thus, in this program code, the variable "StartJ" is initialized back 201 data samples. The variable "CurrJ" is initialized back one data sample. Because a circular array is being used, the "BufLength" variable is referenced during the initialization of both "StartJ" and "CurrJ". The outer Do loop moves the current and preceding sets of data samples relative to each other. The outer Do loop is active while the variable "k" indicates that the program code is processing within the search range. In an embodiment, the search range is initially set at three positions to either side of a predicted position. The predicted position is based upon the period obtained for an immediately prior execution of the program code. If the program code has not been executed immediately prior, then a default period value is used. If an acceptable minimum absolute difference value is not found within this search range, then the search range can be expanded to, for example, 50 positions to either side of the predicted position.
The variable "j" is initialized to the "StartJ" value. The variables "AbsDiff and "n" are also initialized prior to execution of the inner Do loop. The inner Do loop performs the computation of the absolute difference values between the present set and prior set of data samples. The variable "AbsDiff maintains the sum of the absolute difference of values for overlapping data samples being compared. Note that the number of data samples being analyzed to determine the absolute difference values varies based upon the position in the search range being processed. This results because different positions in the search range has different numbers of data samples that overlap with the previous set of data samples being compared. In the embodiment of this program code, the absolute difference function is computed using every 10th signal sample point, i.e., a subsampled subtraction is used.
Because a circular array is being used, the "Chartwidth" variable is referenced during the calculation of "AbsDiff . The variables "MaxSignal" and "MinSignal" indicate a maximum and minimum range for signal values that have previously been sampled. These values can be established, for example, during a learning period for the system in which data samples are obtained for a plurality of respiratory cycles. The "NormAbsDiff variable holds the average absolute difference value represented as a percentage value based upon the "MaxSignal" and "MinSignal" values.
If the "NormAbsDiff value is less than or equal to a previously established "MinAbsDiff value, then the "MinAbsDiff variable is set to the "NormAbsDiff value. The "MaxK" variable is set to the value of "k" if the "MinAbsDiff value is reset. The variable "k" is then incremented, and if the "k" value is still within the search range, then the program code returns back to the beginning of the outer Do loop.
The result of this program code is a candidate minimum absolute difference value and a candidate position for the minimum absolute difference. The MaxK value is compared to predefined threshold values to ensure that it falls within a correct range of values for the physiological activity being processed. Thus, in an embodiment, the MaxK value is tested to make sure that it is greater than or equal to 40 and less than or equal to 150. If the mark value meets the threshold range, then the variable "Period" is set to the "MaxK" value. The variable "Predict" returns the predicted value for the next set of data samples to be processed. The variable "Diff indicates the difference value between the current data sample value and the predicted data sample value, and is represented as a percentage to the "MaxSignal" and "MixSignal" values.
In an embodiment, an image of a progress bar can be displayed to visually indicate the periodicity of the signal samples. According to the program code, if the "MinAbsDiff value is less than or equal to a 20% difference, then the visual progress bar is updated with the computed "MinAbsDiff value. Otherwise, the visual progress bar displays all other "MinAbsDiff values that exceed a 20% difference as a default value of "20".
Figs. 6a and 6b depict an embodiment of a camera 108 that can used in the present invention to optically or visually collect data representative of physiological movement. Camera 108 is a charge-couple device ("CCD") camera having one or more photoelectric cathodes and one or more CCD devices. A CCD device is a semiconductor device that can store charge in local areas, and upon appropriate control signals, transfers that charge to a readout point. When light photons from the scene to be images are focussed on the photoelectric cathodes, electrons are liberated in proportion to light intensity received at the camera. The electrons are captured in charge buckets located within the CCD device. The distribution of captured electrons in the charge buckets represents the image received at the camera. The CCD transfers these electrons to an analog-to-digital converter. The output of the analog-to-digital converter is sent to computer 410 to process the video image and to calculate the positions of the retro-reflective markers 406. According to an embodiment of the invention, camera 108 is a monochrome CCD camera having RS-170 output and 640x480 pixel resolution. Alternatively, camera 408 can comprise a CCD camera having CCIR output and 756x567 pixel resolution.
In a particular embodiment of the invention, an infra-red illuminator 602 ("IR illuminator") is co-located with camera 108. IR illuminator 602 produces one or more beams of infrared light that is directed in the same direction as camera 108. IR illuminator 602 comprises a surface that is ringed around the lens 606 of camera body 608. The surface of IR illuminator 602 contains a plurality of individual LED elements 604 for producing infrared light. The LED elements 604 are arranged in a spiral pattern on the IR illuminator 602. Infrared filters that may be part of the camera 108 are removed or disabled to increase the camera's sensitivity to infrared light. According to an embodiment, digital video recordings of the patient in a session can be recorded via camera 108. The same camera 108 used for tracking patient movement can be used to record video images of the patient for future reference. A normal ambient light image sequence of the patient can be obtained in synchronization with the measured movement signals of markers 114. Figs. 7a and 7b depict an embodiment of a retro-reflective marker 700 that can be employed within the present invention. Retro-reflective marker 700 comprises a raised reflective surface 702 for reflecting light. Raised reflective surface 702 comprises a semi- spherical shape such that light can be reflected regardless of the input angle of the light source. A flat surface 704 surrounds the raised reflective surface 702. The underside of flat surface 704 provides a mounting area to attach retro-reflective marker 700to particular locations on a patient's body. According to an embodiment, retro-reflective marker 700 is comprised of a retro-reflective material 3M#7610WS available from 3M Corporation. In an embodiment, marker 700 has a diameter of approximately .5 cm and a height of the highest point of raised reflective surface 702 of approximately .1 cm. Fig. 8 depicts an apparatus 802 that can be employed to manufacture retro-reflective markers 700. Apparatus 802 comprises a base portion 804 having an elastic ring 806 affixed thereto. Elastic ring 806 is attached to bottom mold piece 808 having a bulge protruding from its center. A control lever 810 can be operated to move top portion 812 along support rods 814.
Top portion 812 comprises a spring-loaded top mold piece 814. Top mold piece 814 is formed with a semi-spherical cavity on its underside. In operation, a piece of retro-reflective material is placed on bottom mold piece 808. Control lever 810 is operated to move top portion 812 towards base portion 804. The retro-reflective material is compressed and shaped between the bottom mold piece 808 and the top mold piece 814. The top mold piece 814 forms the upper exterior of the retro-reflective material into a semi-spherical shape.
In an alternate embodiment, marker 114 comprises a marker block having one or more reference locations on its surface. Each reference location on the marker block preferably comprises a retro-reflective or reflective material that is detectable by an optical imaging apparatus, such as camera 108. Fig. 11 depicts an embodiment of a marker block 1100 having a cylindrical shape with multiple reference locations comprised of retro-reflective elements 1102 located on its surface. Marker block 1100 can be formed as a rigid block (e.g., from Styrofoam). Blocks made in this fashion can be reused a plurality of times, even with multiple patients. The retro-reflective elements 1102 can be formed from the same material used to construct retro-reflective markers 114 of Figs. 7a and 7b. The marker block is preferably formed from a material that is lightweight enough not to interfere with normal breathing by the patient.
A marker block can be formed into any shape or size, as long as the size, spacing, and positioning of the reference locations are configured such that a camera or other optical imaging apparatus can view and generate an image that accurately shows the positioning of the marker block. For example, Fig. 10 depicts an alternate marker block 1000 having a hemispherical shape comprised of a plurality of retro-reflective elements 1002 attached to its surface.
The marker block can be formed with shapes to fit particular body parts. For example, molds or casts that match to specific locations on the body can be employed as marker blocks. Marker blocks shaped to fit certain areas of the body facilitate the repeatable placement of the marker blocks at particular locations on the patient. Alternatively, the marker blocks can be formed to fit certain fixtures that are attached to a patient's body. For example, a marker block can be formed within indentations and grooves that allow it to be attached to eyeglasses. In yet another embodiment, the fixtures are formed with integral marker block(s) having reflective or retro-reflective markers on them. An alternate embodiment of the marker block comprises only a single reference location/reflective element on its surface. This embodiment of the marker block is used in place of the retro-reflective marker 406 to detect particular locations on a patient's body with an optical imaging apparatus. Further details regarding camera 108, markers 114, marker blocks, or procedures to optically measure physiological movements are described in copending U.S. Patent Application Ser. No. [Not Yet Assigned] (Varian Docket No. 98-34, Attorney Docket No. 236/223) and U.S. Patent Application Ser. No. [Not Yet Assigned] (Varian Docket No. 98-33, Attorney Docket No. 236/224), filed concurrently herewith, both of which are hereby incorporated by reference in their entirety.
Figure 12 is a block diagram that illustrates an embodiment of a computer system 1900 upon which an embodiment of the invention may be implemented. Computer system 1900 includes a bus 1902 or other communication mechanism for communicating information, and a processor 1904 coupled with bus 1902 for processing information. Computer system 1900 also includes a main memory 1906, such as a random access memory (RAM) or other dynamic storage device, coupled to bus 1902 for storing information and instructions to be executed by processor 1904. Main memory 1906 also may be used for storing temporary variables or other intermediate information during execution of instructions to be executed by processor 1904. Computer system 1900 further includes a read only memory (ROM) 1908 or other static storage device coupled to bus 1902 for storing static information and instructions for processor 1904. A data storage device 1910, such as a magnetic disk or optical disk, is provided and coupled to bus 1902 for storing information and instructions.
Computer system 1900 may be coupled via bus 1902 to a display 1912, such as a cathode ray tube (CRT), for displaying information to a user. An input device 1914, including alphanumeric and other keys, is coupled to bus 1902 for communicating information and command selections to processor 1904. Another type of user input device is cursor control 1916, such as a mouse, a trackball, or cursor direction keys for communicating direction information and command selections to processor 1904 and for controlling cursor movement on display 1912. This input device typically has two degrees of freedom in two axes, a first axis (e.g., x) and a second axis (e.g., y), that allows the device to specify positions in a plane.
The invention is related to the use of computer system 1900 for detecting and predictively estimating physiological cycles. According to one embodiment of the invention, such use is provided by computer system 1900 in response to processor 1904 executing one or more sequences of one or more instructions contained in main memory 1906. Such instructions may be read into main memory 1906 from another computer-readable medium, such as storage device 1910. Execution of the sequences of instructions contained in main memory 1906 causes processor 1904 to perform the process steps described herein. One or more processors in a multi-processing arrangement may also be employed to execute the sequences of instructions contained in main memory 1906. In alternative embodiments, hard- wired circuitry may be used in place of or in combination with software instructions to implement the invention. Thus, embodiments of the invention are not limited to any specific combination of hardware circuitry and software. The term "computer-readable medium" as used herein refers to any medium that participates in providing instructions to processor 1904 for execution. Such a medium may take many forms, including but not limited to, non-volatile media, volatile media, and transmission media. Non- volatile media includes, for example, optical or magnetic disks, such as storage device 1910. Volatile media includes dynamic memory, such as main memory 1906. Transmission media includes coaxial cables, copper wire and fiber optics, including the wires that comprise bus 1902. Transmission media can also take the form of acoustic or light waves, such as those generated during radio wave and infrared data communications.
Common forms of computer-readable media include, for example, a floppy disk, a flexible disk, hard disk, magnetic tape, or any other magnetic medium, a CD-ROM, any other optical medium, punch cards, paper tape, any other physical medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, a carrier wave as described hereinafter, or any other medium from which a computer can read.
Various forms of computer readable media may be involved in carrying one or more sequences of one or more instructions to processor 1904 for execution. For example, the instructions may initially be carried on a magnetic disk of a remote computer. The remote computer can load the instructions into its dynamic memory and send the instructions over a telephone line using a modem. A modem local to computer system 1900 can receive the data on the telephone line and use an infrared transmitter to convert the data to an infrared signal. An infrared detector coupled to bus 1902 can receive the data carried in the infrared signal and place the data on bus 1902. Bus 1902 carries the data to main memory 1906, from which processor 1904 retrieves and executes the instructions. The instructions received by main memory 1906 may optionally be stored on storage device 1910 either before or after execution by processor 1904.
Computer system 1900 also includes a communication interface 1918 coupled to bus 1902. Communication interface 1918 provides a two-way data communication coupling to a network link 1920 that is connected to a local network 1922. For example, communication interface 1918 may be an integrated services digital network (ISDN) card or a modem to provide a data communication connection to a corresponding type of telephone line. As another example, communication interface 1918 may be a local area network (LAN) card to provide a data communication connection to a compatible LAN. Wireless links may also be implemented. In any such implementation, communication interface 1918 sends and receives electrical, electromagnetic or optical signals that carry data streams representing various types of information. Network link 1920 typically provides data communication through one or more networks to other devices. For example, network link 1920 may provide a connection through local network 1922 to a host computer 1924 or to medical equipment 1926 such as a radiation beam source or a switch operatively coupled to a radiation beam source. The data streams transported over network link 1920 can comprise electrical, electromagnetic or optical signals. The signals through the various networks and the signals on network link 1920 and through communication interface 1918, which carry data to and from computer system 1900, are exemplary forms of carrier waves transporting the information. Computer system 1900 can send messages and receive data, including program code, through the network(s), network link 1920 and communication interface 1918. In the foregoing specification, the invention has been described with reference to specific embodiments thereof. It will, however, be evident that various modifications and changes may be made thereto without departing from the broader spirit and scope of the invention. For example, the operations performed by computer 110 can be performed by any combination of hardware and software within the scope of the invention, and should not be limited to particular embodiments comprising just a particular definition of "computer". The specification and drawings are, accordingly, to be regarded in an illustrative rather than restrictive sense.

Claims

CLAIMSWhat is claimed is:
1. A method of applying radiation comprising: measuring a first set of signal data representative of a physiological movement of a patient during a first time period; pattern matching the first set of signal data with a second set of signal data related to measured physiological movement of a patient during a second time period; and applying radiation to the patient based upon results of the pattern matching.
2. The method of claim 1 in which the first set of signal data and the second set of signal data are pattern matched using an autocorrelation function.
3. The method of claim 1 in which the first set of signal data and the second set of signal data are pattern matched using an absolute difference function.
4. The method of claim 1 further comprising: determining a degree of match between the first set of signal data and the second set of signal data.
5. The method of claim 4 in which the degree of match is determined by a secondary peak value of an autocorrelation function.
6. The method of claim 4 in which the degree of match is determined by a secondary minimum value of an absolute difference function.
7. The method of claim 4 further comprising: comparing the degree of match to a threshold range.
8. The method of claim 7 in which the degree of match outside the threshold range indicates deviation from a normal physiological movement.
9. The method of claim 7 in which the degree of match within the threshold range indicates a repetitive physiological movement.
10. The method of claim 9 in which a point of best match indicates a period of the physiological movement.
11. The method of claim 1 further comprising: predicting a period of the physiological movement during a third time period.
12. The method of claim 11 further comprising: predictively actuating a gating system component based upon the predicted period.
13. The method of claim 1 further comprising: determining a period of the physiological movement.
14. The method of claim 13 further comprising; defining a treatment interval to apply radiation to a patient.
15. The method of claim 14 in which the treatment interval is defined based upon a phase range of the period of the physiological movement.
16. The method of claim 1 in which the second set of signal data is a data model of the physiological movement of the patient.
17. A method for applying radiation comprising: determining a period of physiological movement; applying radiation based upon a phase range of the period of the physiological movement.
18. The method of claim 17 in which the period of physiological movement is determined by pattern matching a first set of data representative of the physiological movement during a first time period with a second set of data related to the physiological movement during a second time period.
19. The method of claim 18 in which the first set of data and the second set of data are pattern matched using an autocorrelation function.
20. The method of claim 18 in which the first set of data and the second set of data are pattern matched using an absolute difference function.
21. A method of applying radiation comprising: determining an amount of time required to actuate a gating component; predictively actuating the gating component to compensate for the amount of time required to actuate the gating component.
22. The method of claim 21 in which said gating component is a switch operatively coupled to a radiation source.
23. The method of claim 21 further comprising: predicting a next period of physiological movement.
24. The method of claim 23 in which a treatment interval is defined over the physiological movement.
25. The method of claim 24 in which the gating component is predictively actuated to coincide the full actuation of the gating component with a boundary of the treatment interval.
PCT/US1999/024946 1998-10-23 1999-10-22 Method and system for predictive physiological gating of radiation therapy WO2000024466A1 (en)

Priority Applications (6)

Application Number Priority Date Filing Date Title
EP99970937A EP1123137B1 (en) 1998-10-23 1999-10-22 Method and system for predictive physiological gating of radiation therapy
CA002348091A CA2348091A1 (en) 1998-10-23 1999-10-22 Method and system for predictive physiological gating of radiation therapy
JP2000578067A JP4391023B2 (en) 1998-10-23 1999-10-22 System for predictive physiological gating of radiation therapy
DE69929628T DE69929628T2 (en) 1998-10-23 1999-10-22 METHOD AND SYSTEM FOR THE PREDICTIVE CONTROL OF RADIOTHERAPY
KR1020017005049A KR20010083921A (en) 1998-10-23 1999-10-22 Method and system for predictive physiological gating of radiation therapy
AU12286/00A AU771038B2 (en) 1998-10-23 1999-10-22 Method and system for predictive physiological gating of radiation therapy

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US09/178,383 1998-10-23
US09/178,383 US6621889B1 (en) 1998-10-23 1998-10-23 Method and system for predictive physiological gating of radiation therapy

Publications (1)

Publication Number Publication Date
WO2000024466A1 true WO2000024466A1 (en) 2000-05-04

Family

ID=22652337

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US1999/024946 WO2000024466A1 (en) 1998-10-23 1999-10-22 Method and system for predictive physiological gating of radiation therapy

Country Status (9)

Country Link
US (3) US6621889B1 (en)
EP (1) EP1123137B1 (en)
JP (1) JP4391023B2 (en)
KR (1) KR20010083921A (en)
AT (1) ATE316403T1 (en)
AU (1) AU771038B2 (en)
CA (1) CA2348091A1 (en)
DE (1) DE69929628T2 (en)
WO (1) WO2000024466A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2820045A1 (en) 2001-01-29 2002-08-02 Chabunda Christophe Mwanza Stereo-radiotherapy apparatus for bicephalic imaging uses acceleration and collision of particles from two synchronous radiation beams
JP2004533889A (en) * 2001-06-26 2004-11-11 バリアン・メディカル・システムズ・インコーポレイテッド Method and system for predictive physiological gate drive
US10646188B2 (en) 1998-10-23 2020-05-12 Varian Medical Systems, Inc. Method and system for radiation application
US10667727B2 (en) 2008-09-05 2020-06-02 Varian Medical Systems, Inc. Systems and methods for determining a state of a patient

Families Citing this family (175)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6279579B1 (en) * 1998-10-23 2001-08-28 Varian Medical Systems, Inc. Method and system for positioning patients for medical treatment procedures
US7158610B2 (en) * 2003-09-05 2007-01-02 Varian Medical Systems Technologies, Inc. Systems and methods for processing x-ray images
US6973202B2 (en) * 1998-10-23 2005-12-06 Varian Medical Systems Technologies, Inc. Single-camera tracking of an object
US6980679B2 (en) * 1998-10-23 2005-12-27 Varian Medical System Technologies, Inc. Method and system for monitoring breathing activity of a subject
US6621889B1 (en) * 1998-10-23 2003-09-16 Varian Medical Systems, Inc. Method and system for predictive physiological gating of radiation therapy
FR2823679B1 (en) * 2001-04-23 2004-01-02 Dyn R METHOD AND DEVICE FOR DRIVING A RADIATION THERAPY UNIT FOR THE TREATMENT OF THORACIC TUMORS
US7769430B2 (en) * 2001-06-26 2010-08-03 Varian Medical Systems, Inc. Patient visual instruction techniques for synchronizing breathing with a medical procedure
WO2003039212A1 (en) 2001-10-30 2003-05-08 Loma Linda University Medical Center Method and device for delivering radiotherapy
US7221733B1 (en) * 2002-01-02 2007-05-22 Varian Medical Systems Technologies, Inc. Method and apparatus for irradiating a target
US20040002641A1 (en) * 2002-06-24 2004-01-01 Bo Sjogren Patient representation in medical machines
US7227925B1 (en) 2002-10-02 2007-06-05 Varian Medical Systems Technologies, Inc. Gantry mounted stereoscopic imaging system
WO2004066211A1 (en) * 2003-01-21 2004-08-05 Elekta Ab (Publ) Imaging internal structures
US20040199073A1 (en) * 2003-04-03 2004-10-07 Agency For Science, Technology And Research Method and apparatus for measuring motion of a body in a number of dimensions
US7367955B2 (en) * 2003-06-13 2008-05-06 Wisconsin Alumni Research Foundation Combined laser spirometer motion tracking system for radiotherapy
US7199382B2 (en) 2003-08-12 2007-04-03 Loma Linda University Medical Center Patient alignment system with external measurement and object coordination for radiation therapy system
WO2005018735A2 (en) 2003-08-12 2005-03-03 Loma Linda University Medical Center Modular patient support system
US8571639B2 (en) 2003-09-05 2013-10-29 Varian Medical Systems, Inc. Systems and methods for gating medical procedures
US7217913B2 (en) * 2003-12-18 2007-05-15 Micron Technology, Inc. Method and system for wavelength-dependent imaging and detection using a hybrid filter
US7542544B2 (en) * 2004-01-06 2009-06-02 The Regents Of The University Of Michigan Ultrasound gating of cardiac CT scans
WO2005081842A2 (en) 2004-02-20 2005-09-09 University Of Florida Research Foundation, Inc. System for delivering conformal radiation therapy while simultaneously imaging soft tissue
US7536219B2 (en) * 2004-02-20 2009-05-19 Siemens Aktiengesellschaft 4D imaging with a C-arm X-ray system
US8457717B2 (en) * 2004-04-08 2013-06-04 Stanford University Method and system of adaptive control for reducing motion artifacts and patient dose in four dimensional computed tomography
US8214012B2 (en) * 2004-06-17 2012-07-03 Psychology Software Tools, Inc. Magnetic resonance imaging having patient video, microphone and motion tracking
US7073508B2 (en) 2004-06-25 2006-07-11 Loma Linda University Medical Center Method and device for registration and immobilization
ES2558978T3 (en) 2004-07-21 2016-02-09 Mevion Medical Systems, Inc. Programmable radiofrequency waveform generator for a synchro-cyclotron
US20060074305A1 (en) * 2004-09-30 2006-04-06 Varian Medical Systems Technologies, Inc. Patient multimedia display
US7957507B2 (en) 2005-02-28 2011-06-07 Cadman Patrick F Method and apparatus for modulating a radiation beam
US7983380B2 (en) 2005-04-29 2011-07-19 Varian Medical Systems, Inc. Radiation systems
US8232535B2 (en) 2005-05-10 2012-07-31 Tomotherapy Incorporated System and method of treating a patient with radiation therapy
US7349522B2 (en) * 2005-06-22 2008-03-25 Board Of Trustees Of The University Of Arkansas Dynamic radiation therapy simulation system
US8442287B2 (en) 2005-07-22 2013-05-14 Tomotherapy Incorporated Method and system for evaluating quality assurance criteria in delivery of a treatment plan
CA2616304A1 (en) 2005-07-22 2007-02-01 Tomotherapy Incorporated System and method of delivering radiation therapy to a moving region of interest
WO2007014108A2 (en) 2005-07-22 2007-02-01 Tomotherapy Incorporated Method and system for evaluating quality assurance criteria in delivery of a treament plan
WO2007014104A2 (en) 2005-07-22 2007-02-01 Tomotherapy Incorporated System and method of evaluating dose delivered by a radiation therapy system
KR20080039916A (en) 2005-07-22 2008-05-07 토모테라피 인코포레이티드 System and method of delivering radiation therapy to a moving region of interest
JP5390855B2 (en) 2005-07-23 2014-01-15 トモセラピー・インコーポレーテッド Imaging and delivery of radiation therapy using coordinated movement of gantry and treatment table
US9119541B2 (en) * 2005-08-30 2015-09-01 Varian Medical Systems, Inc. Eyewear for patient prompting
CN101361156B (en) 2005-11-18 2012-12-12 梅维昂医疗系统股份有限公司 Charged particle radiation therapy
US7656258B1 (en) 2006-01-19 2010-02-02 Massachusetts Institute Of Technology Magnet structure for particle acceleration
EP1977632A2 (en) * 2006-01-19 2008-10-08 Massachusetts Institute Of Technology High-field superconducting synchrocyclotron
US8060177B2 (en) * 2006-03-06 2011-11-15 Siemens Medical Solutions Usa, Inc. Registration of computed tomography (CT) and positron emission tomography (PET) image scans with automatic patient motion correction
US7587022B1 (en) * 2006-03-23 2009-09-08 General Electric Company Correlation-based motion estimation of object to be imaged
US8195278B2 (en) * 2006-05-15 2012-06-05 Siemens Medical Solutions Usa, Inc. Method for automatically determining an image plane having a biopsy device therein
PL2023812T3 (en) 2006-05-19 2017-07-31 The Queen's Medical Center Motion tracking system for real time adaptive imaging and spectroscopy
DE502006003264D1 (en) * 2006-06-27 2009-05-07 Brainlab Ag Medical-technical marker tracking with determination of marker properties
GB2441550A (en) * 2006-09-05 2008-03-12 Vision Rt Ltd Surface-imaging breathing monitor
US9451928B2 (en) 2006-09-13 2016-09-27 Elekta Ltd. Incorporating internal anatomy in clinical radiotherapy setups
US7620147B2 (en) * 2006-12-13 2009-11-17 Oraya Therapeutics, Inc. Orthovoltage radiotherapy
US7535991B2 (en) 2006-10-16 2009-05-19 Oraya Therapeutics, Inc. Portable orthovoltage radiotherapy
US7894649B2 (en) * 2006-11-02 2011-02-22 Accuray Incorporated Target tracking using direct target registration
EP2088925B8 (en) 2006-11-17 2015-06-17 Varian Medical Systems, Inc. Dynamic patient positioning system
EP2095374A4 (en) 2006-11-21 2012-05-30 Univ Loma Linda Med Device and method for immobilizing patients for breast radiation therapy
JP5348870B2 (en) * 2006-11-24 2013-11-20 株式会社東芝 MRI equipment
US7609810B2 (en) * 2006-12-14 2009-10-27 Byong Yong Yi Treatment-speed regulated tumor-tracking
WO2008086434A2 (en) * 2007-01-09 2008-07-17 Cyberheart, Inc. Depositing radiation in heart muscle under ultrasound guidance
WO2008086430A1 (en) * 2007-01-09 2008-07-17 Cyberheart, Inc. Method for depositing radiation in heart muscle
US10716528B2 (en) 2007-03-08 2020-07-21 Sync-Rx, Ltd. Automatic display of previously-acquired endoluminal images
US11064964B2 (en) 2007-03-08 2021-07-20 Sync-Rx, Ltd Determining a characteristic of a lumen by measuring velocity of a contrast agent
US9968256B2 (en) 2007-03-08 2018-05-15 Sync-Rx Ltd. Automatic identification of a tool
US9375164B2 (en) 2007-03-08 2016-06-28 Sync-Rx, Ltd. Co-use of endoluminal data and extraluminal imaging
US9629571B2 (en) 2007-03-08 2017-04-25 Sync-Rx, Ltd. Co-use of endoluminal data and extraluminal imaging
JP5639764B2 (en) * 2007-03-08 2014-12-10 シンク−アールエックス,リミティド Imaging and tools for use with moving organs
US8542900B2 (en) * 2007-03-08 2013-09-24 Sync-Rx Ltd. Automatic reduction of interfering elements from an image stream of a moving organ
US11197651B2 (en) 2007-03-08 2021-12-14 Sync-Rx, Ltd. Identification and presentation of device-to-vessel relative motion
US20080317204A1 (en) * 2007-03-16 2008-12-25 Cyberheart, Inc. Radiation treatment planning and delivery for moving targets in the heart
US10974075B2 (en) 2007-03-16 2021-04-13 Varian Medical Systems, Inc. Radiation treatment planning and delivery for moving targets in the heart
US8512236B2 (en) 2008-01-11 2013-08-20 Oraya Therapeutics, Inc. System and method for positioning and stabilizing an eye
US8363783B2 (en) 2007-06-04 2013-01-29 Oraya Therapeutics, Inc. Method and device for ocular alignment and coupling of ocular structures
WO2009012577A1 (en) 2007-07-20 2009-01-29 Resonant Medical Inc. Methods and systems for compensating for changes in anatomy of radiotherapy patients
WO2009012576A1 (en) 2007-07-20 2009-01-29 Resonant Medical Inc. Methods and systems for guiding the acquisition of ultrasound images
US8135198B2 (en) 2007-08-08 2012-03-13 Resonant Medical, Inc. Systems and methods for constructing images
US7986227B2 (en) * 2007-09-20 2011-07-26 Cornell Research Foundation, Inc. System and method for position matching of a patient for medical imaging
US8175676B2 (en) 2007-10-01 2012-05-08 National Health Research Institutes Generating gating signals for thermal therapy
US8003964B2 (en) 2007-10-11 2011-08-23 Still River Systems Incorporated Applying a particle beam to a patient
US8002465B2 (en) * 2007-11-19 2011-08-23 Pyronia Medical Technologies, Inc. Patient positioning system and methods for diagnostic radiology and radiotherapy
US8933650B2 (en) 2007-11-30 2015-01-13 Mevion Medical Systems, Inc. Matching a resonant frequency of a resonant cavity to a frequency of an input voltage
US8581523B2 (en) 2007-11-30 2013-11-12 Mevion Medical Systems, Inc. Interrupted particle source
US7801271B2 (en) 2007-12-23 2010-09-21 Oraya Therapeutics, Inc. Methods and devices for orthovoltage ocular radiotherapy and treatment planning
US7792249B2 (en) 2007-12-23 2010-09-07 Oraya Therapeutics, Inc. Methods and devices for detecting, controlling, and predicting radiation delivery
DE102008019114A1 (en) * 2008-04-16 2009-10-22 Siemens Aktiengesellschaft Medical device and method for docking a positioning device to a shuttle
EP2110160A1 (en) 2008-04-16 2009-10-21 Ion Beam Applications S.A. Device for selecting a beam triggering apparatus
US8189738B2 (en) 2008-06-02 2012-05-29 Elekta Ltd. Methods and systems for guiding clinical radiotherapy setups
ES2450391T3 (en) 2008-06-19 2014-03-24 Sync-Rx, Ltd. Progressive progress of a medical instrument
JP2010069099A (en) * 2008-09-19 2010-04-02 Toshiba Corp Image processing apparatus and x-ray computed tomography apparatus
US9974509B2 (en) 2008-11-18 2018-05-22 Sync-Rx Ltd. Image super enhancement
US11064903B2 (en) 2008-11-18 2021-07-20 Sync-Rx, Ltd Apparatus and methods for mapping a sequence of images to a roadmap image
US9095313B2 (en) 2008-11-18 2015-08-04 Sync-Rx, Ltd. Accounting for non-uniform longitudinal motion during movement of an endoluminal imaging probe
US10362962B2 (en) 2008-11-18 2019-07-30 Synx-Rx, Ltd. Accounting for skipped imaging locations during movement of an endoluminal imaging probe
US9101286B2 (en) 2008-11-18 2015-08-11 Sync-Rx, Ltd. Apparatus and methods for determining a dimension of a portion of a stack of endoluminal data points
US9144394B2 (en) 2008-11-18 2015-09-29 Sync-Rx, Ltd. Apparatus and methods for determining a plurality of local calibration factors for an image
US8855744B2 (en) 2008-11-18 2014-10-07 Sync-Rx, Ltd. Displaying a device within an endoluminal image stack
CA2744005C (en) * 2008-11-21 2017-04-25 Tohoku University Signal processing apparatus, signal processing method, signal processing program, computer-readable recording medium storing signal processing program, and radiotherapy apparatus
US8337512B2 (en) * 2008-12-08 2012-12-25 Siemens Aktiengesellschaft Device and workflow for minimally-invasive therapy, in particular needle guidance
US8632448B1 (en) 2009-02-05 2014-01-21 Loma Linda University Medical Center Proton scattering analysis system
EP2384142B1 (en) 2009-02-12 2012-10-31 Vision RT Limited Patient monitor and method
JP2010216968A (en) * 2009-03-16 2010-09-30 Fuji Xerox Co Ltd System and program for measuring position
US10542962B2 (en) 2009-07-10 2020-01-28 Elekta, LTD Adaptive radiotherapy treatment using ultrasound
AU2010274004A1 (en) 2009-07-17 2012-03-08 Cyberheart, Inc. Heart treatment kit, system, and method for radiosurgically alleviating arrhythmia
WO2011030460A1 (en) * 2009-09-14 2011-03-17 三菱重工業株式会社 Tomographic method
US8400149B2 (en) * 2009-09-25 2013-03-19 Nellcor Puritan Bennett Ireland Systems and methods for gating an imaging device
US20110110570A1 (en) * 2009-11-10 2011-05-12 Avi Bar-Shalev Apparatus and methods for generating a planar image
US9146289B2 (en) * 2009-12-23 2015-09-29 General Electric Company Targeted thermal treatment of human tissue through respiratory cycles using ARMA modeling
US20110172526A1 (en) * 2010-01-12 2011-07-14 Martin Lachaine Feature Tracking Using Ultrasound
US9248316B2 (en) 2010-01-12 2016-02-02 Elekta Ltd. Feature tracking using ultrasound
US9000914B2 (en) * 2010-03-15 2015-04-07 Welch Allyn, Inc. Personal area network pairing
BR112013012005A2 (en) * 2010-11-17 2016-08-30 Koninkl Philips Electronics Nv image projection system for use in a patient examination scan room and method for enhancing patient comfort in a patient examination scan room through projected images on a room surface
AU2012259403B2 (en) 2011-03-07 2016-08-04 Loma Linda University Medical Center Systems, devices and methods related to calibration of a proton computed tomography scanner
JP6099640B2 (en) 2011-06-23 2017-03-22 シンク−アールエックス,リミティド Lumen background sharpening
US9606209B2 (en) 2011-08-26 2017-03-28 Kineticor, Inc. Methods, systems, and devices for intra-scan motion correction
US10561861B2 (en) 2012-05-02 2020-02-18 Viewray Technologies, Inc. Videographic display of real-time medical treatment
JP6134789B2 (en) 2012-06-26 2017-05-24 シンク−アールエックス,リミティド Image processing related to flow in luminal organs
CN104813747B (en) 2012-09-28 2018-02-02 梅维昂医疗系统股份有限公司 Use magnetic field flutter focused particle beam
WO2014052734A1 (en) 2012-09-28 2014-04-03 Mevion Medical Systems, Inc. Controlling particle therapy
JP6121546B2 (en) 2012-09-28 2017-04-26 メビオン・メディカル・システムズ・インコーポレーテッド Control system for particle accelerator
JP6121544B2 (en) 2012-09-28 2017-04-26 メビオン・メディカル・システムズ・インコーポレーテッド Particle beam focusing
CN108770178B (en) 2012-09-28 2021-04-16 迈胜医疗设备有限公司 Magnetic field regenerator
EP3581242B1 (en) 2012-09-28 2022-04-06 Mevion Medical Systems, Inc. Adjusting energy of a particle beam
JP6367201B2 (en) 2012-09-28 2018-08-01 メビオン・メディカル・システムズ・インコーポレーテッド Control of particle beam intensity
US10254739B2 (en) 2012-09-28 2019-04-09 Mevion Medical Systems, Inc. Coil positioning system
EP2901824B1 (en) 2012-09-28 2020-04-15 Mevion Medical Systems, Inc. Magnetic shims to adjust a position of a main coil and corresponding method
EP2911745B1 (en) 2012-10-26 2019-08-07 ViewRay Technologies, Inc. Assessment and improvement of treatment using imaging of physiological responses to radiation therapy
US9305365B2 (en) 2013-01-24 2016-04-05 Kineticor, Inc. Systems, devices, and methods for tracking moving targets
US9717461B2 (en) 2013-01-24 2017-08-01 Kineticor, Inc. Systems, devices, and methods for tracking and compensating for patient motion during a medical imaging scan
US10327708B2 (en) 2013-01-24 2019-06-25 Kineticor, Inc. Systems, devices, and methods for tracking and compensating for patient motion during a medical imaging scan
WO2014120734A1 (en) 2013-02-01 2014-08-07 Kineticor, Inc. Motion tracking system for real time adaptive motion compensation in biomedical imaging
JP6412020B2 (en) 2013-02-26 2018-10-24 アキュレイ インコーポレイテッド Electromagnetic multi-leaf collimator
US9351698B2 (en) 2013-03-12 2016-05-31 Lightlab Imaging, Inc. Vascular data processing and image registration systems, methods, and apparatuses
US9446263B2 (en) 2013-03-15 2016-09-20 Viewray Technologies, Inc. Systems and methods for linear accelerator radiotherapy with magnetic resonance imaging
US8791656B1 (en) 2013-05-31 2014-07-29 Mevion Medical Systems, Inc. Active return system
US9730308B2 (en) 2013-06-12 2017-08-08 Mevion Medical Systems, Inc. Particle accelerator that produces charged particles having variable energies
CN110237447B (en) 2013-09-27 2021-11-02 梅维昂医疗系统股份有限公司 Particle therapy system
US10675487B2 (en) 2013-12-20 2020-06-09 Mevion Medical Systems, Inc. Energy degrader enabling high-speed energy switching
US9962560B2 (en) 2013-12-20 2018-05-08 Mevion Medical Systems, Inc. Collimator and energy degrader
US9661736B2 (en) 2014-02-20 2017-05-23 Mevion Medical Systems, Inc. Scanning system for a particle therapy system
US9959631B2 (en) 2014-02-21 2018-05-01 Samsung Electronics Co., Ltd. Tomography apparatus and method for reconstructing tomography image thereof
WO2015126189A1 (en) * 2014-02-21 2015-08-27 Samsung Electronics Co., Ltd. Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus
CN106061376B (en) * 2014-03-03 2020-07-28 瓦里安医疗系统公司 System and method for patient position monitoring
WO2015148391A1 (en) 2014-03-24 2015-10-01 Thomas Michael Ernst Systems, methods, and devices for removing prospective motion correction from medical imaging scans
WO2016014718A1 (en) 2014-07-23 2016-01-28 Kineticor, Inc. Systems, devices, and methods for tracking and compensating for patient motion during a medical imaging scan
US9629565B2 (en) 2014-08-18 2017-04-25 Cameron Health, Inc. Peak selection for self correlation analysis of cardiac rate in an implantable medical devices
US9950194B2 (en) 2014-09-09 2018-04-24 Mevion Medical Systems, Inc. Patient positioning system
US10499813B2 (en) 2014-09-12 2019-12-10 Lightlab Imaging, Inc. Methods, systems and apparatus for temporal calibration of an intravascular imaging system
KR101816886B1 (en) 2014-10-22 2018-01-09 사회복지법인 삼성생명공익재단 System and method for inducing respiration
US10105107B2 (en) 2015-01-08 2018-10-23 St. Jude Medical International Holding S.À R.L. Medical system having combined and synergized data output from multiple independent inputs
KR20160117907A (en) * 2015-04-01 2016-10-11 가톨릭대학교 산학협력단 Monitoring system for thermoplastic mask using pressure sensor
US10222956B2 (en) 2015-05-17 2019-03-05 Lightlab Imaging, Inc. Intravascular imaging user interface systems and methods
US10109058B2 (en) 2015-05-17 2018-10-23 Lightlab Imaging, Inc. Intravascular imaging system interfaces and stent detection methods
US10646198B2 (en) 2015-05-17 2020-05-12 Lightlab Imaging, Inc. Intravascular imaging and guide catheter detection methods and systems
US9996921B2 (en) 2015-05-17 2018-06-12 LIGHTLAB IMAGING, lNC. Detection of metal stent struts
JP6867385B2 (en) 2015-07-25 2021-04-28 ライトラボ・イメージング・インコーポレーテッド Intravascular data visualization method
US9943247B2 (en) 2015-07-28 2018-04-17 The University Of Hawai'i Systems, devices, and methods for detecting false movements for motion correction during a medical imaging scan
US10786689B2 (en) 2015-11-10 2020-09-29 Mevion Medical Systems, Inc. Adaptive aperture
US9847210B2 (en) * 2015-11-11 2017-12-19 Mitsubishi Electric Corporation Particle beam irradiation apparatus for irradiating a subject with an arbitrary number of particles
WO2017087821A2 (en) 2015-11-18 2017-05-26 Lightlab Imaging, Inc. X-ray image feature detection and registration systems and methods
CN115998310A (en) 2015-11-23 2023-04-25 光学实验室成像公司 Detection and verification of shadows in intravascular images
WO2017091479A1 (en) 2015-11-23 2017-06-01 Kineticor, Inc. Systems, devices, and methods for tracking and compensating for patient motion during a medical imaging scan
US9966160B2 (en) 2015-11-24 2018-05-08 Viewray Technologies, Inc. Radiation beam collimating systems and methods
US10300303B2 (en) * 2016-01-29 2019-05-28 Elekta Ltd. Therapy control using motion prediction based on cyclic motion model
KR20180120705A (en) 2016-03-02 2018-11-06 뷰레이 테크놀로지스 인크. Particle therapy using magnetic resonance imaging
ES2908571T3 (en) 2016-04-14 2022-05-03 Lightlab Imaging Inc Identification of branches of a blood vessel
WO2017201026A1 (en) 2016-05-16 2017-11-23 Lightlab Imaging, Inc. Intravascular absorbable stent detection and diagnostic methods and systems
AU2017281519A1 (en) 2016-06-22 2019-01-24 Viewray Technologies, Inc. Magnetic resonance imaging at low field strength
US10027410B2 (en) * 2016-06-23 2018-07-17 Abl Ip Holding Llc System and method using a gated retro-reflector for visible light uplink communication
JP7059245B2 (en) 2016-07-08 2022-04-25 メビオン・メディカル・システムズ・インコーポレーテッド Decide on a treatment plan
US10806409B2 (en) 2016-09-23 2020-10-20 Varian Medical Systems International Ag Medical systems with patient supports
CN110382049A (en) 2016-12-13 2019-10-25 优瑞技术公司 Radiotherapy system and method
US11103730B2 (en) 2017-02-23 2021-08-31 Mevion Medical Systems, Inc. Automated treatment in particle therapy
EP3586309B1 (en) * 2017-02-24 2023-04-05 Brainlab AG Deep inspiration breath-hold setup using x-ray imaging
US11273326B2 (en) * 2017-06-29 2022-03-15 Canon Medical Systems Corporation Radiotherapy system and treatment support apparatus
JP6940676B2 (en) 2017-06-30 2021-09-29 メビオン・メディカル・システムズ・インコーポレーテッド Configurable collimator controlled using a linear motor
US11033758B2 (en) 2017-12-06 2021-06-15 Viewray Technologies, Inc. Radiotherapy systems, methods and software
RU182872U1 (en) * 2018-04-19 2018-09-05 Максим Валерьевич Ткачев DEVICE FOR LABELING TUMOR BEDS
US11209509B2 (en) 2018-05-16 2021-12-28 Viewray Technologies, Inc. Resistive electromagnet systems and methods
WO2020124204A1 (en) * 2018-12-21 2020-06-25 Dalhousie University Apparatus and method for gating delivery of radiation based on capacitive monitoring of respiratory motion
EP3934752A1 (en) 2019-03-08 2022-01-12 Mevion Medical Systems, Inc. Delivery of radiation by column and generating a treatment plan therefor

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4463425A (en) * 1980-07-17 1984-07-31 Terumo Corporation Period measurement system
US5446548A (en) * 1993-10-08 1995-08-29 Siemens Medical Systems, Inc. Patient positioning and monitoring system
US5538494A (en) * 1994-03-17 1996-07-23 Hitachi, Ltd. Radioactive beam irradiation method and apparatus taking movement of the irradiation area into consideration
US5727554A (en) * 1996-09-19 1998-03-17 University Of Pittsburgh Of The Commonwealth System Of Higher Education Apparatus responsive to movement of a patient during treatment/diagnosis
WO1998016151A1 (en) * 1996-10-16 1998-04-23 The Trustees Of Columbia University In The City Of New York Apparatus and method to gate a source for radiation therapy

Family Cites Families (70)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US262306A (en) * 1882-08-08 Animal-tie
US912656A (en) * 1908-10-26 1909-02-16 Charles F Engstrom Album and sheets therefor.
US3861807A (en) 1972-08-17 1975-01-21 Charles Lescrenier Position locating and maintaining method and means
GB1475308A (en) * 1973-07-21 1977-06-01 Emi Ltd Radiography
US3871360A (en) 1973-07-30 1975-03-18 Brattle Instr Corp Timing biological imaging, measuring, and therapeutic timing systems
FR2273505A1 (en) 1974-06-07 1976-01-02 Inst Nat Sante Rech Med RESPIRATORY AND HEART CYCLE CORRELATION DEVICE, AND APPLICATION TO HEART RATE MEASUREMENT
DE2718804C3 (en) 1977-04-27 1979-10-31 Karlheinz Prof. Dr. 3000 Hannover Renner Device for positioning control of patients and / or radiation sources
US4971065A (en) 1985-02-11 1990-11-20 Pearce Stephen D Transducer for detecting apnea
ATE58462T1 (en) 1985-05-23 1990-12-15 Heinrich Prof Dr Ing Reents DEVICE FOR MEASURING THE VITAL FUNCTIONS OF A HUMAN, ESPECIALLY AN INFANT.
FI79458C (en) 1986-04-16 1990-01-10 Seppo Peltola FOERFARANDE OCH ANORDNING FOER STYRNING AV EN BEHANDLINGSANORDNING.
US4853771A (en) 1986-07-09 1989-08-01 The United States Of America As Represented By The Secretary Of The Navy Robotic vision system
US4710717A (en) * 1986-12-29 1987-12-01 General Electric Company Method for fast scan cine NMR imaging
FR2637189A1 (en) 1988-10-04 1990-04-06 Cgr Mev SYSTEM AND METHOD FOR MEASURING AND / OR VERIFYING THE POSITION OF A PATIENT IN RADIOTHERAPY EQUIPMENT
US4994965A (en) 1988-11-23 1991-02-19 General Electric Company Method for reducing motion induced image artifacts in projection imaging
US5295483A (en) 1990-05-11 1994-03-22 Christopher Nowacki Locating target in human body
US6405072B1 (en) 1991-01-28 2002-06-11 Sherwood Services Ag Apparatus and method for determining a location of an anatomical target with reference to a medical apparatus
US5662111A (en) 1991-01-28 1997-09-02 Cosman; Eric R. Process of stereotactic optical navigation
US5279309A (en) 1991-06-13 1994-01-18 International Business Machines Corporation Signaling device and method for monitoring positions in a surgical operation
US5262945A (en) 1991-08-09 1993-11-16 The United States Of America As Represented By The Department Of Health And Human Services Method for quantification of brain volume from magnetic resonance images
US5619995A (en) 1991-11-12 1997-04-15 Lobodzinski; Suave M. Motion video transformation system and method
DE4207632C2 (en) 1992-03-11 1995-07-20 Bodenseewerk Geraetetech Device and method for positioning a body part for treatment purposes
US5603318A (en) 1992-04-21 1997-02-18 University Of Utah Research Foundation Apparatus and method for photogrammetric surgical localization
US5389101A (en) 1992-04-21 1995-02-14 University Of Utah Apparatus and method for photogrammetric surgical localization
US5265142A (en) 1992-05-08 1993-11-23 General Electric Company Image reconstruction technique for a computer tomography system
US5271055A (en) * 1992-08-19 1993-12-14 General Electric Company Methods for reducing motion induced artifacts in a projection imaging system
US5513646A (en) 1992-11-09 1996-05-07 I Am Fine, Inc. Personal security monitoring system and method
JPH0723945A (en) 1993-07-07 1995-01-27 Toshiba Corp Medical image photographing device
JP3393895B2 (en) * 1993-09-13 2003-04-07 株式会社東芝 Magnetic resonance imaging
US5394875A (en) 1993-10-21 1995-03-07 Lewis; Judith T. Automatic ultrasonic localization of targets implanted in a portion of the anatomy
DE4341324A1 (en) 1993-12-03 1995-06-08 Siemens Ag Medical equipment for diagnosis and/or therapy
US5912656A (en) 1994-07-01 1999-06-15 Ohmeda Inc. Device for producing a display from monitored data
US5622187A (en) 1994-09-30 1997-04-22 Nomos Corporation Method and apparatus for patient positioning for radiation therapy
US5582182A (en) 1994-10-03 1996-12-10 Sierra Biotechnology Company, Lc Abnormal dyspnea perception detection system and method
US5588430A (en) 1995-02-14 1996-12-31 University Of Florida Research Foundation, Inc. Repeat fixation for frameless stereotactic procedure
GB9515311D0 (en) 1995-07-26 1995-09-20 3D Scanners Ltd Stripe scanners and methods of scanning
DE19529639C2 (en) * 1995-08-11 1997-06-19 Siemens Ag Process for the temporal and location-resolved representation of functional brain activities of a patient and arrangement for carrying out the process
US5638819A (en) 1995-08-29 1997-06-17 Manwaring; Kim H. Method and apparatus for guiding an instrument to a target
JP3597918B2 (en) 1995-09-11 2004-12-08 株式会社日立メディコ X-ray CT system
US5771310A (en) 1996-12-30 1998-06-23 Shriners Hospitals For Children Method and apparatus for recording three-dimensional topographies
US5823192A (en) 1996-07-31 1998-10-20 University Of Pittsburgh Of The Commonwealth System Of Higher Education Apparatus for automatically positioning a patient for treatment/diagnoses
US5820553A (en) 1996-08-16 1998-10-13 Siemens Medical Systems, Inc. Identification system and method for radiation therapy
US6296613B1 (en) 1997-08-22 2001-10-02 Synthes (U.S.A.) 3D ultrasound recording device
US5784431A (en) 1996-10-29 1998-07-21 University Of Pittsburgh Of The Commonwealth System Of Higher Education Apparatus for matching X-ray images with reference images
FR2760277B1 (en) 1997-02-28 1999-03-26 Commissariat Energie Atomique METHOD AND DEVICE FOR LOCATING AN OBJECT IN SPACE
EP1011424A1 (en) 1997-03-03 2000-06-28 Schneider Medical Technologies, Inc. Imaging device and method
DE19732784C1 (en) 1997-07-30 1999-03-04 Bruker Medizintech Positioning system and method for exact position determination of a manually operated manipulator in an MR tomograph
US6434507B1 (en) 1997-09-05 2002-08-13 Surgical Navigation Technologies, Inc. Medical instrument and method for use with computer-assisted image guided surgery
DE19743500A1 (en) 1997-10-01 1999-04-29 Siemens Ag Medical apparatus with device for detecting position of object
US6061644A (en) 1997-12-05 2000-05-09 Northern Digital Incorporated System for determining the spatial position and orientation of a body
US6348058B1 (en) 1997-12-12 2002-02-19 Surgical Navigation Technologies, Inc. Image guided spinal surgery guide, system, and method for use thereof
US5993397A (en) 1998-01-23 1999-11-30 Branson; Krista Lynn Infant respiratory monitor
US6076005A (en) 1998-02-25 2000-06-13 St. Jude Children's Research Hospital Respiration responsive gating means and apparatus and methods using the same
US6198959B1 (en) 1998-03-27 2001-03-06 Cornell Research Foundation Inc. Coronary magnetic resonance angiography using motion matched acquisition
US6185446B1 (en) 1998-08-21 2001-02-06 William F. Carlsen, Jr. Method and apparatus for monitoring the breathing of a patient during magnetic resonance imaging
US6144874A (en) * 1998-10-15 2000-11-07 General Electric Company Respiratory gating method for MR imaging
US6973202B2 (en) 1998-10-23 2005-12-06 Varian Medical Systems Technologies, Inc. Single-camera tracking of an object
US6279579B1 (en) 1998-10-23 2001-08-28 Varian Medical Systems, Inc. Method and system for positioning patients for medical treatment procedures
US6980679B2 (en) 1998-10-23 2005-12-27 Varian Medical System Technologies, Inc. Method and system for monitoring breathing activity of a subject
US6621889B1 (en) * 1998-10-23 2003-09-16 Varian Medical Systems, Inc. Method and system for predictive physiological gating of radiation therapy
US6937696B1 (en) * 1998-10-23 2005-08-30 Varian Medical Systems Technologies, Inc. Method and system for predictive physiological gating
US6138302A (en) 1998-11-10 2000-10-31 University Of Pittsburgh Of The Commonwealth System Of Higher Education Apparatus and method for positioning patient
US6501981B1 (en) 1999-03-16 2002-12-31 Accuray, Inc. Apparatus and method for compensating for respiratory and patient motions during treatment
US6144875A (en) 1999-03-16 2000-11-07 Accuray Incorporated Apparatus and method for compensating for respiratory and patient motion during treatment
US6370217B1 (en) 1999-05-07 2002-04-09 General Electric Company Volumetric computed tomography system for cardiac imaging
US6661617B1 (en) 1999-12-14 2003-12-09 Seagate Technology Llc Structure and fabrication process for integrated moving-coil magnetic micro-actuator
US7058440B2 (en) 2001-06-28 2006-06-06 Koninklijke Philips Electronics N.V. Dynamic computed tomography imaging using positional state modeling
US7120679B2 (en) 2001-06-29 2006-10-10 Intel Corporation Configuration of headless devices using configuration service
DE10133237B4 (en) 2001-07-09 2007-04-19 Siemens Ag Method for computed tomography and computed tomography (CT) device
US7620444B2 (en) 2002-10-05 2009-11-17 General Electric Company Systems and methods for improving usability of images for medical applications
WO2004066211A1 (en) 2003-01-21 2004-08-05 Elekta Ab (Publ) Imaging internal structures

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4463425A (en) * 1980-07-17 1984-07-31 Terumo Corporation Period measurement system
US5446548A (en) * 1993-10-08 1995-08-29 Siemens Medical Systems, Inc. Patient positioning and monitoring system
US5538494A (en) * 1994-03-17 1996-07-23 Hitachi, Ltd. Radioactive beam irradiation method and apparatus taking movement of the irradiation area into consideration
US5727554A (en) * 1996-09-19 1998-03-17 University Of Pittsburgh Of The Commonwealth System Of Higher Education Apparatus responsive to movement of a patient during treatment/diagnosis
WO1998016151A1 (en) * 1996-10-16 1998-04-23 The Trustees Of Columbia University In The City Of New York Apparatus and method to gate a source for radiation therapy

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KUBO H D ET AL: "RESPIRATION GATED RADIOTHERAPY TREATMENT: A TECHNICAL STUDY", PHYSICS IN MEDICINE AND BIOLOGY,GB,TAYLOR AND FRANCIS LTD. LONDON, vol. 41, no. 1, 1 January 1996 (1996-01-01), pages 83 - 91, XP000548771, ISSN: 0031-9155 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9232928B2 (en) 1998-10-23 2016-01-12 Varian Medical Systems, Inc. Method and system for predictive physiological gating
US10646188B2 (en) 1998-10-23 2020-05-12 Varian Medical Systems, Inc. Method and system for radiation application
FR2820045A1 (en) 2001-01-29 2002-08-02 Chabunda Christophe Mwanza Stereo-radiotherapy apparatus for bicephalic imaging uses acceleration and collision of particles from two synchronous radiation beams
JP2004533889A (en) * 2001-06-26 2004-11-11 バリアン・メディカル・システムズ・インコーポレイテッド Method and system for predictive physiological gate drive
US10667727B2 (en) 2008-09-05 2020-06-02 Varian Medical Systems, Inc. Systems and methods for determining a state of a patient

Also Published As

Publication number Publication date
US6621889B1 (en) 2003-09-16
US6959266B1 (en) 2005-10-25
US7191100B2 (en) 2007-03-13
JP2002528193A (en) 2002-09-03
EP1123137A1 (en) 2001-08-16
JP4391023B2 (en) 2009-12-24
DE69929628D1 (en) 2006-04-13
CA2348091A1 (en) 2000-05-04
KR20010083921A (en) 2001-09-03
US20060004547A1 (en) 2006-01-05
ATE316403T1 (en) 2006-02-15
AU771038B2 (en) 2004-03-11
EP1123137B1 (en) 2006-01-25
DE69929628T2 (en) 2006-09-21
AU1228600A (en) 2000-05-15

Similar Documents

Publication Publication Date Title
US6621889B1 (en) Method and system for predictive physiological gating of radiation therapy
EP1402761B1 (en) Method and system for predictive physiological gating
US10646188B2 (en) Method and system for radiation application
EP1123138B1 (en) Method and system for physiological gating of radiation therapy
EP2155039B1 (en) Marker systems
US20040092815A1 (en) Method and apparatus for tracking an internal target region without an implanted fiducial
WO2005089194A2 (en) Breathing synchronized computed tomography image acquisition
AU2003252856A1 (en) Predictive Physiological Gating of Radiation Therapy

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref country code: AU

Ref document number: 2000 12286

Kind code of ref document: A

Format of ref document f/p: F

AK Designated states

Kind code of ref document: A1

Designated state(s): AU CA IL JP KR MX

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
WWE Wipo information: entry into national phase

Ref document number: 1999970937

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2348091

Country of ref document: CA

Ref country code: CA

Ref document number: 2348091

Kind code of ref document: A

Format of ref document f/p: F

Ref country code: JP

Ref document number: 2000 578067

Kind code of ref document: A

Format of ref document f/p: F

WWE Wipo information: entry into national phase

Ref document number: 1020017005049

Country of ref document: KR

WWE Wipo information: entry into national phase

Ref document number: 12286/00

Country of ref document: AU

WWP Wipo information: published in national office

Ref document number: 1999970937

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 1020017005049

Country of ref document: KR

WWG Wipo information: grant in national office

Ref document number: 12286/00

Country of ref document: AU

WWW Wipo information: withdrawn in national office

Ref document number: 1020017005049

Country of ref document: KR

WWG Wipo information: grant in national office

Ref document number: 1999970937

Country of ref document: EP