Matched filter in ionoacoustic signal processing for ion beam range determination and dosimetry

Patent metadata

Field	Value
Publication number	US-12329990-B2
Application number	US-202318176047-A
Country	US
Kind code	B2
Filing date	Feb 28, 2023
Priority date	Mar 3, 2022
Publication date	Jun 17, 2025
Grant date	Jun 17, 2025

Abstract

A method and an apparatus for determining information regarding the location of energy deposition of an ion beam, in particular a proton beam, in an absorptive medium, such as in the tissue of a patient undergoing radiation therapy.

First claim

What is claimed is:

1. A method of determining information regarding a location of energy deposition of an ion beam in an absorptive medium, comprising the following steps:
generating a pulsed ion beam,
detecting a time-resolved acoustic signal, said time-resolved acoustic signal comprising an energy-deposition-signal component attributable to the energy deposition of an individual pulse of said pulsed ion beam in said absorptive medium,
determining relative timing information of the energy deposition of said individual pulse with respect to said time-resolved acoustic signal,
providing a matched filter for processing the time-resolved acoustic signal, said matched filter being configured to facilitate detecting said energy-deposition-signal component within said time-resolved acoustic signal,
applying said matched filter to said time-resolved acoustic signal to thereby obtain a filtered time-resolved signal,
deriving, from said filtered time-resolved signal, occurrence timing information related to an occurrence of said energy-deposition signal component in said time-resolved acoustic signal, and
deriving information regarding the location of the energy deposition based, at least in part, on a delay between said energy deposition of said individual pulse and said occurrence timing information.

2. The method of claim 1 , wherein said absorptive medium is a tissue of a patient undergoing radiation therapy.

3. The method of claim 1 , wherein said step of providing the matched filter comprises providing a template, said template representing a predicted waveform of the energy-deposition-signal component.

4. The method of claim 3 , wherein applying said matched filter to said time-resolved acoustic signal comprises correlating said template with said time-resolved acoustic signal, or convolving the time-resolved acoustic signal with a conjugated time-reversed version of the template.

5. The method of claim 3 , wherein the step of providing said template comprises determining the template at least in part by computer simulation, wherein said simulation comprises simulating the template as an ionoacoustic signal caused by an estimated time-dependent energy deposition distribution of the ion beam in said absorptive medium.

6. The method of claim 5 , further comprising a step of obtaining said estimated time-dependent energy deposition by computer simulation as well.

7. The method of claim 5 , wherein said simulation of said template by said ionoacoustic signal is based on space-resolved information about mass density and/or speed of sound in a tissue of a patient, and the space-resolved information is obtained from medical images obtained during one or both of a treatment planning and during radiation therapy.

8. The method of claim 3 , wherein said step of providing said template further comprises accounting for transducer characteristics of a transducer used for detecting the time-resolved acoustic signal, wherein the transducer characteristics are accounted for using a transfer function associated with said transducer, or a bandpass filter approximating the transducer's transfer function.

9. The method of claim 3 , wherein said step of providing said template comprises optimizing the template in an iterative procedure, to thereby increase the signal-to-noise ratio (SNR) of the filtered time resolved signal.

10. The method of claim 3 , wherein providing said template comprises selecting one or more templates from a database of previously generated templates, wherein said previously generated templates comprise templates which are, at least in part, based on energy-deposition-signal components obtained by experiment using doses higher than 40 Gy.

11. The method of claim 10 , wherein said one or more selected templates are further modified to account for one or both of a shape of the individual pulse and the characteristics of a sensor used for detecting the time-resolved acoustic signal.

12. The method of claim 1 , wherein said relative timing information of the energy deposition of said individual pulse with respect to said time-resolved acoustic signal is determined based on a trigger signal indicating an arrival of a pulse of said pulsed ion beam in the absorptive medium, wherein said trigger signal is one of
an electric signal provided by electronic components involved with generating said pulsed ion beam,
a signal provided by a detector provided in front of the absorptive medium, and
a signal provided by a detector for detecting secondary radiation due to an interaction of the ion beam with the absorptive medium.

13. The method of claim 1 , wherein said time-resolved acoustic signal is detected using at least one detection apparatus,
wherein said at least one detection apparatus comprises an ionoacoustic transducer for detecting said time-resolved acoustic signal, and wherein
said ionoacoustic transducer is adapted to record ultrasonic images of said absorptive medium as well, or
wherein an ultrasound transducer is provided in a known spatial relationship to said ionoacoustic transducer.

14. The method of claim 13 , wherein said location of the energy deposition is derived with respect to a coordinate system associated with an ultrasound image recorded with said ionoacoustic transducer or with said ultrasound transducer in said known spatial relationship to said ionoacoustic transducer further comprising a step of providing a medical image in which said location of energy deposition is indicated, wherein said medical image is one of said ultrasound image, an image derived from said ultrasound image, and a medical image acquired with an imaging modality different from ultrasound imaging that is co-registered with said ultrasound image.

15. The method of claim 14 , wherein a target area for energy deposition according to a treatment plan is associated with said medical image, said method further comprising a step of determining whether the location of energy deposition with regard to said target area deviates from said treatment plan, and in case the deviation exceeds a predetermined threshold, deriving control or operating parameters such as to decrease the deviation, wherein the control or operating parameters relate to energy of the ion beam or to positioning parameters for a treatment table on which a patient is placed.

16. The method of claim 1 , wherein the ions are protons, and a pulse duration t d given as the full width at half maximum (FWHM) of a Gaussian fitted in a least square fit to an actual pulse shape is chosen as a function of proton energy E P , such that
ƒ( E P )≤ t d ≤g ( E P ),

wherein t d is measured in μs, E P is measured in MeV and chosen from a range between 20 MeV and 260 MeV, and ƒ(E P ), g(E P ) are second order polynomials defined as
ƒ( E P )= p 1( E P ) 2 +p 2 E P +p 3 and

g ( E P )= q 1( E P ) 2 +q 2 E P +q 3, respectively,

wherein the coefficients p i and q i with i=1 . . . 3 are chosen as follows:
p1=1.476*10 −5 , p2=0.002486, p3=−0.05185, and
q1=0.000118, q2=0.03548, q3=−0.9488.

17. The method of claim 16 , wherein for a given dose limit of protons of proton energy E P to be irradiated into tissue of a patient, the dose is split into two or more pulses, each of which, have a pulse duration t d .

18. An apparatus for determining information regarding the location of energy deposition of a pulsed ion beam in an absorptive medium, said apparatus comprising:
a detection apparatus for detecting a time-resolved acoustic signal, said time-resolved acoustic signal compris

Assignees

• Univ Der Bundeswehr Muenchen
• Univ Muenchen Tech
• Univ Muenchen Ludwig Maximilians

Inventors

• Schauer Jannis
• Dollinger Günther
• Ruser Heinrich
• Huang Yuanhui
• Ntziachristos Vasilis
• Wieser Hans Peter
• Parodi Katia

Classifications

• A61N2005/1041

  using a library of previously administered radiation treatment applied to other patients · CPC title

• A61N2005/1087

  Ions; Protons · CPC title

• A61B8/4488

  the transducer being a phased array · CPC title

• A61N5/1071 Primary

  for verifying the dose delivered by the treatment plan · CPC title