Space-time microwave imaging for cancer detection

• US 7,570,063 B2
• Filed: 07/03/2002
• Issued: 08/04/2009
• Est. Priority Date: 07/06/2001
• Status: Active Grant

• Current Assignee: Wisconsin Alumni Research Foundation (University of Wisconsin System)
• Original Assignee: Wisconsin Alumni Research Foundation (University of Wisconsin System)
• Inventors: Li, Xu, Van Veen, Barry D., Hagness, Susan C., Bond, Essex Julian
• Primary Examiner(s): Casler; Brian L.
• Assistant Examiner(s): Cattungal; Sanjay
• Application Number:
  US10/190,352
• Publication Number:
  US 20030088180A1
• Time in Patent Office: 2,589 Days
• Field of Search: 600/407, 600/408, 600/430, 606 27- 52, 607154-156, 607/101, 324637-648
• US Class Current: 324/637
• CPC Class Codes: A61B 18/1815 using microwaves
• CPC Class Codes: A61B 5/05 Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves impedance plethysmography A61B5/0295; measuring movement of the entire body or parts thereof A61B5/11; detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof A61B5/24
• CPC Class Codes: A61B 5/0507 using microwaves or terahertz waves
• CPC Class Codes: A61B 5/7257 using Fourier transforms

Abstract

Microwave imaging via space-time beamforming is carried out by transmitting microwave signals from multiple antenna locations into an individual to be examined and receiving the backscattered microwave signals at multiple antenna locations to provide received signals from the antennas. The received signals are processed in a computer to remove the skin interface reflection component of the signal at each antenna to provide corrected signal data. The corrected signal data is provided to a beamformer process that time shifts the received signals to align the returns from a scatterer at a candidate location, and then passes the time aligned signals through a bank of filters, the outputs of which are summed, time-gated and the power therein calculated to produce the beamformer output signal at a candidate location. The beamformer is then scanned to a plurality of different locations in the individual by changing the time shifts, filter weights and time-gating of the beamformer process. The output power may be displayed as a function of scan location, with regions of large output power corresponding to significant microwave scatterers such as malignant lesions.

Claims

1. A space-time microwave imaging system for imaging an object comprising:

(a) an array of antennas for radiating and receiving microwaves;

(b) a microwave pulse source connected to the array of antennas to provide pulse signals to the antennas;

(c) a receiver connected to the antennas to detect the microwave signals received by the antennas and provide signal data corresponding thereto; and

(d) a computer connected to receive the signal data, the computer programmed to estimate an artifact reflection component of the signal at each antenna as a filtered combination of the signals at all other antennas and to subtract the estimated artifact reflection component from the signal data to provide corrected signal data, with weights of the filters chosen to minimize the residual signal over that portion of the received data dominated by the reflection, and the computer further programmed to process the corrected signal data in a beamformer process to time shift the corrected signal data to approximately align the returns from a scatterer at a candidate location, to pass the time aligned signals through a bank of filters with a filter for each antenna, to sum the outputs of the filters to form a summed signal, and to calculate the power in the summed signal to produce a beamformer output signal, wherein the beamformer filters are designed to present maximum gain to scattered signals originating from the candidate location, and to scan the beamformer process to a plurality of different candidate locations in the object to be imaged by changing the time shifts and filter weights to generate multi-dimensional output power data.

2. The imaging system of claim 1 including an output device connected to the computer to display the multi-dimensional output power data as a function of scanned locations.

3. The imaging system of claim 1 wherein the computer is further programmed to time gate the summed signal to form a time gated summed signal and to calculate the power in the time gated summed signal, and to scan the beamformer process to a plurality of different candidate locations in the object to be imaged by changing the time shifts, filter weights and time gates to generate multi-dimensional output power data.

4. The imaging system of claim 3 wherein the computer is further programmed to process the beamformer output signal from the filters prior to time gating using a parametric signal processing model to compensate frequency dependent scattering effects.

5. The imaging system of claim 1 wherein the microwave pulse source provides pulses having pulse widths on the order of 100 picoseconds or less in duration.

6. The imaging system of claim 1 wherein the microwave pulse source is connected to the antennas to provide pulse signals to one antenna at a time.

7. The imaging system of claim 1 wherein the microwave pulse source is connected to the antennas to provide pulse signals to all of the antennas simultaneously.

8. The imaging system of claim 1 wherein the computer is programmed in the beamformer process to apply a selected window to the time aligned signals before passing the time aligned signals through the bank of filters, and further to apply a selected window to the summed signal before the power in the summed signal is calculated.

9. The imaging system of claim 1 wherein the computer is programmed to calculate the power at each candidate location as the sum of the squares of each sample of the summed signal.

10. The imaging system of claim 1 wherein the beamformer filters are finite-impulse response filters designed to satisfy a penalized least squares condition to present unit gain to scattered signals originating from a candidate location.

11. The imaging system of claim 1 further including signal processing circuitry that receives the pulses from the source and passes the pulses through a delay and a filter for each antenna before providing the delayed and filtered pulses to the antennas, the delays and filters selected to focus the radiated microwave energy from the array of antennas at a selected candidate location in the object.

12. A method of carrying out space-time microwave imaging of an individual comprising:

(a) transmitting microwave signals from a plurality of antenna locations into an individual to be examined;

(b) receiving backscattered microwave signals at a plurality of antenna locations to provide received signals from the plurality of antenna locations;

(c) processing the received signals in a computer in a beamformer process by time shifting the received signal data to approximately align the returns from a scatterer at a candidate location, passing the time aligned signals through a bank of filters with a filter for each antenna location, summing the outputs of the filters to provide a summed signal, and equalizing the summed signal to reduce a sensitivity to a time-gate duration; and

time-gating the equalized signal to provide a time-gated signal, and calculating the power in the time-gated signal to produce a beamformer output signal, wherein the filters are designed to present maximum gain to scattered signals originating from the candidate location; and

(d) then scanning the beamformer process to a plurality of different candidate locations in the individual and repeating steps (a), (b) and (c) at each candidate location by changing the time shifts, filter weights and time gates for each candidate location to generate multi-dimensional output power data.

13. The method of claim 12 wherein the step of transmitting microwave signals comprises transmitting microwave pulses having pulse widths on the order of 100 picoseconds or less in duration.

14. The method of claim 12 wherein the step of transmitting microwave signals comprises transmitting microwave pulses having frequency content at 10 GHz or higher.

15. The method of claim 12 wherein the step of filtering is carried out in the time domain.

16. The method of claim 12 wherein the step of filtering is carried out in the frequency domain.

17. The method of claim 16 wherein the step of filtering in the frequency domain is carried out for each signal from an antenna location in a channel by performing a fast Fourier transform on the signal data to provide frequency domain data, applying weights to the frequency domain data in each channel and adding the weighted data;

and then performing an inverse fast Fourier transform on the weighted data to provide a filtered time domain data signal.

18. The method of claim 12 including compensating the received signal data for frequency dependent scattering before carrying out beamformer processing.

19. The method of claim 12 including transmitting the microwave signals from an array of antennas so as to focus the microwave power on a candidate location.

20. The method of claim 12 further including, before the step of processing the received signals in a beamformer process, the step of:

carrying out an artifact response subtraction process on the received signals in the computer by estimating the interface reflection component of the signal at each antenna location as a combination of the received signals at the other antenna locations passed through filters to provide corrected signal data, the filters having weights chosen to minimize the received signal over that portion of the received signal dominated by the artifact response, and providing the corrected signal data to the beamformer process.

21. The method of claim 12 wherein microwave signals are provided to one antenna at a time and backscattered microwave signals are received from one antenna at a time for each of the antenna locations.

22. The method of claim 12 wherein microwave signals are transmitted from all of the antennas simultaneously and backscattered microwave signals are received from all of the antennas simultaneously.

23. The method of claim 12 including in the beamformer process the steps of applying a selected window to the time aligned signals before filtering the signals and further applying a selected window to the summed signal before calculating the power in the summed signal.

24. The method of claim 12 wherein the power is calculated at each candidate location as the sum of the squares of each sample of the summed signal.

25. The method of claim 12 wherein filtering is carried out utilizing finite-impulse response filters designed to satisfy a penalized least squares condition to present unit gain to scattered signals originating from a candidate location.

26. The method of claim 12 wherein the step of transmitting the microwave signals is carried out simultaneously from all of the antenna locations by passing microwave pulses for each antenna at an antenna location through a delay and a filter for each antenna, the delays and filters selected to focus the radiated microwave energy from the antennas at a selected candidate location in the object.

27. A method of compensating for an artifact reflection in microwave imaging of an individual comprising:

(a) transmitting microwave signals from a plurality of antenna locations into an individual to be examined;

(b) receiving backscattered microwave signals at a plurality of antenna locations to provide received signals from the plurality of antenna locations;

(c) carrying out an artifact response subtraction process on the received signals in a computer by estimating the artifact reflection component of the signal at each antenna location as a combination of the received signals at the other antenna locations passed through filters to provide corrected signal data, the filters having weights chosen to minimize the received signal over that portion of the received signal dominated by the artifact response; and

(d) compensating the corrected signal data for frequency dependent scattering.

28. The method of claim 27 wherein the step of transmitting microwave signals comprises transmitting microwave pulses having pulse widths on the order of 100 picoseconds or less in duration.

29. The method of claim 27 wherein the step of transmitting microwave signals comprises transmitting microwave pulses having frequency content at 10 GHz or higher.

30. The method of claim 27 including transmitting the microwave signals from an array of antennas so as to focus the microwave power on a candidate location.

31. A method of carrying out space-time microwave imaging of an individual comprising:

(a) transmitting microwave signals from a plurality of antenna locations into an individual to be examined;

(b) receiving backscattered microwave signals at a plurality of antenna locations to provide received signals from the plurality of antenna locations;

(c) carrying out an artifact response subtraction process on the received signals in the computer by estimating the artifact reflection component of the signal at each antenna location as a combination of the received signals at the other antenna locations passed through filters to provide corrected signal data, the filter weights chosen to minimize the received signal over that portion of the received signal dominated by the artifact response;

(d) processing the corrected signal data in a computer in a beamformer process by time shifting the received signal data to approximately align the returns from a scatterer at a candidate location, passing the time aligned signals through a bank of filters with one filter in each channel, summing the outputs of the filters to provide a summed signal, and time-gating the summed signal to provide a time-gated signal, and calculating the power in the time-gated signal to produce a beamformer output signal, wherein the filters are designed to present maximum gain to scattered signals originating from the candidate location; and

(e) then scanning the beamformer process to a plurality of different candidate locations in the individual and repeating steps (a), (b) (c) and (d) at each candidate location by changing the time shifts, filter weights and time gates for each candidate location to generate multi-dimensional output power data.

32. The method of claim 31 wherein the step of transmitting microwave signals comprises transmitting microwave pulses having pulse widths on the order of 100 picoseconds or less in duration.

33. The method of claim 31 wherein the step of transmitting microwave signals comprises transmitting microwave pulses having frequency content at 10 GHz or higher.

34. The method of claim 31 wherein the step of filtering is carried out in the time domain.

35. The method of claim 31 wherein the step of filtering is carried out in the frequency domain.

36. The method of claim 35 wherein the step of filtering in the frequency domain is carried out for each signal from an antenna location in a channel by performing a fast Fourier transform on the signal data to provide frequency domain data, applying weights to the frequency domain data in each channel and adding the weighted data;

and then performing an inverse fast Fourier transform on the weighted data to provide a filtered time domain data signal.

37. The method of claim 31 including compensating the received signal data for frequency dependent scattering before carrying out beamformer processing.

38. The method of claim 31 including transmitting the microwave signals from an array of antennas so as to focus the microwave power on a candidate location.

39. The method of claim 31 wherein microwave signals are provided to one antenna at a time and backscattered microwave signals are received from one antenna at a time for each of the antenna locations.

40. The method of claim 31 wherein microwave signals are transmitted from all of the antennas simultaneously and backscattered microwave signals arc received from all of the antennas simultaneously.

41. The method of claim 31 including in the beamformer process the steps of applying a selected window to the time aligned signals before filtering the signals and further applying a selected window to the summed signal before calculating the power in the summed signal.

42. The method of claim 31 wherein the power is calculated at each candidate location as the sum of the squares of each sample of the summed signal.

43. The method of claim 31 wherein filtering is carried out utilizing finite-impulse response filters designed to satisfy a penalized least squares condition to present unit gain to scattered signals originating from a candidate location.

44. The method of claim 31 wherein the step of transmitting the microwave signals is carried out simultaneously from all of the antenna locations by passing microwave pulses for each antenna at an antenna location through a delay and a Filter for each antenna, the delays and filters selected to focus the radiated microwave energy from the antennas at a selected candidate location in the object.

45. An artifact reflection compensation system for use in space-time microwave imaging comprising:

(a) an array of antennas for radiating and receiving microwaves;

(b) a microwave pulse source connected to the array of antennas to provide pulse signals to the antennas;

(c) a receiver connected to the antennas to detect the microwave signals received by the antennas and provide signal data corresponding thereto; and

(d) a computer connected to receive the signal data, the computer programmed to estimate an artifact reflection component of the signal at each antenna as a filtered combination of the signals at all other antennas, to subtract the estimated artifact reflection component from the signal data to provide corrected signal data, and to compensate the corrected signal data for frequency dependent scattering.

46. The system of claim 45 wherein the microwave pulse source provides pulses having pulse widths on the order of 100 picoseconds or less.