Electrical property enhanced tomography (EPET) apparatus and method

• US 6,522,910 B1
• Filed: 03/08/2000
• Issued: 02/18/2003
• Est. Priority Date: 09/11/1997
• Status: Expired due to Term

• Current Assignee: West Virginia University Research Corporation (West Virginia University)
• Original Assignee: Wisys Technology Foundation Incorporated
• Inventors: Gregory, William D.
• Primary Examiner(s): Smith, Ruth S.
• Application Number:
  US09/508,250
• Time in Patent Office: 1,077 Days
• Field of Search: 600/425, 600/427, 600/547, 600/508, 600/428, 600/509, 600/521, 600/529
• US Class Current: 600/427
• CPC Class Codes: A61B 5/0536 Impedance imaging, e.g. by tomography
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: G01N 27/24 Investigating the presence of flaws

Abstract

Enhanced information on the electrical properties of the subregions of a sample is generated by an enhanced tomography (EPET) apparatus connected to a tomographic device, wherein the EPET apparatus holds a voltage constant and measures the electrical current of the sample inserted within the EPET apparatus. A matching medium is applied to the sample such that when the sample is placed within an electromagnetic field generated by the EPET apparatus, the matching medium contacts the segmented sensor plates of the EPET apparatus and allows the detection of the net total charges Q on the surface of the sample. The EPET apparatus uses these net total charges Q in the calculation of the additional information on the electrical properties of the sample by using the scale factor method or the iterative correction method. The EPET apparatus and method has application in the medical imaging field, such as in the obtaining of enhanced dielectric constant and conductivity electrical properties values related to the imaging of the interior of the human body, e.g., in the determination of whether tissue is cancerous, dead or healthy.

Claims

1. A method for providing enhanced information of the electrical properties of a sample measured with a tomographic device and an electrical property enhanced tomography (EPET) apparatus, comprising the steps of:

(a) generating an electromagnetic field in a sample holder by applying frequency dependent voltages of predetermined amplitude to an array of sensors that form a surface to be placed in electrical connection with the sample, such that time dependent net total electric charges are produced on the surface of the sample holder;

(b) determining the time dependent net total electric charges on the surface of the sample holder, wherein the sample holder is filled with a matching medium and the sample is not in the sample holder, thereby generating empty data values;

(c) inserting the sample into the sample holder wherein the matching medium is disposed between the sample and the sample holder to provide electrical connection between the sample and the sensors;

(d) inputting initial image data of the sample from the tomographic device;

(e) determining the time dependent net total electric charges on the surface of the sample holder when the sample is within the sample holder, thereby generating full data values;

(f) generating from said empty data values, said initial image data and said full data values enhanced information on the electrical properties of the sample; and

(g) outputting said enhanced information.

2. The method according to claim 1, wherein said step (f) uses the scale factor method, comprising the steps of:

(f)(1) generating first Fourier transform values 1004 for said empty data values;

(f)(2) generating second Fourier transform values 1004 for said full data values;

(f)(3) computing the difference 1004 between said second Fourier transform values and said first Fourier transform values, thereby generating third Fourier transform values;

(f)(4) inverting said third Fourier transform values 1004, thereby generating inverted Fourier transform values;

(f)(5) determining scale factors 1006 for said inverted Fourier transform values; and

(f)(6) generating enhanced information 1008 of the electrical properties of the sample 204 based on said scale factors.

3. The method according to claim 2, wherein said step (f)(4) uses singular value decomposition.

4. The method according to claim 2, wherein said step (f)(6) uses neural network analysis.

5. The method according to claim 2, wherein said step (f)(6) uses a look up table.

6. The method according to claim 1, wherein said step (f) uses the iterative correction method, comprising the steps of:

(f)(1) inputting 1104 image data on the sample 204 and a guess of the electrical properties values of the sample 204;

(f)(2) generating experimental Fourier transform values 1106;

(f)(3) generating a Fourier transform derivative matrix 1108;

(f)(4) generating theoretical Fourier transform values 1110;

(f)(5) calculating corrections 1112 to said guess of the electrical properties values of the sample 204,

(f)(6) determining whether said corrections are within an acceptable limit 1114; and

(f)(7) generating a new guess of the electrical properties values 116 of the sample 204 by adding said corrections to said guess electrical properties and returning to step (f)(4), if it is determined in step (f)(6) that said corrections are not within said acceptable limit.

7. The method according to claim 1, further comprising the step of:

(h) triggering said step (e) upon the completion of said step (d).

8. The method according to claim 1, further comprising the step of:

(h) triggering said step (e) when noise associated with the sample 204 is minimized.

9. The method according to claim 8, wherein the sample 204 is a human subject and said triggering of said step (h) occurs at an appropriate cardiac cycle of the human subject.

10. The method according to claim 8, wherein the sample 204 is a human subject and said triggering of said step (h) occurs at an appropriate pulmonary cycle of the human subject.

11. A computer program product for use with a computer system of an Electrical Property Enhanced Tomography (EPET) apparatus for providing information of the electrical properties of a sample being measured with a tomographic device, which comprises:

a computer usable medium having computer readable program code means embodied in said medium for generating object code for a computer, said computer readable program code means having;

a first computer readable program code means for operating the computer system to control an electromagnetic field in a sample holder by applying frequency dependent voltages of predetermined amplitude to an array of sensors that form a surface to be placed in electrical connection with the sample, such that time dependent net total electric charges are produced on the surface of the sample holder;

a second computer readable program code means for operating the computer system to determine the time dependent net total electric charges on the surface of the sample holder while the sample is not present in the sample holder to thereby generate empty data values;

a third computer readable program code means for operating the computer system to input initial Image data from the tomographic device when the sample is placed in the sample holder;

a fourth computer readable program code means for operating the computer system to determine the time dependent net total electrical charges on the surface of the sample holder when the sample is within the sample holder, thereby generating full data values;

a fifth computer readable program code means for operating the computer system in response to the generated empty data values, the received initial image data, and the generated full data values to produce enhanced information on the electrical properties of the sample; and

a sixth computer readable program code means for operating the computer system to output said enhanced information on the electrical properties of the sample.

12. The computer program product according to claim 11, wherein said fifth computer readable program code means, comprises:

a means for generating first Fourier transform values 1004 for said empty data values;

a means for generating second Fourier transform values 1004 for said full data values;

a means for computing the difference 1004 between said second Fourier transform values and said first Fourier transform values, thereby generating third Fourier transform values;

a means for inverting said third Fourier transform values 1004, thereby generating inverted Fourier transform values;

a means for determining scale factors 1006 for said inverted Fourier transform values; and

a means for generating said enhanced information 1008 of the electrical properties of the sample 204 based on said scale factors.

13. The computer program product according to claim 12, wherein said means for inverting uses singular value decomposition.

14. The computer program product according to claim 12, wherein said means for generating said enhanced information 1008 uses neural network analysis.

15. The computer program product according to claim 12, wherein said means for generating said enhanced information 1008 uses a look up table.

16. The computer program product according to claim 11, wherein said means for inverting uses the iterative correction method, comprising the steps of:

a means for inputting 1104 image data on the sample 204 and a guess of the electrical properties values of the sample 204;

a means for generating experimental Fourier transform values 1106;

a means for generating a Fourier transform derivative matrix 1108;

a means for generating theoretical Fourier transform values 1110;

a means for calculating corrections 1112 to said guess of the electrical properties values of the sample 204;

a means for determining whether said corrections are within an acceptable limit 1114; and

a means for generating a new guess of the electrical properties values 1116 of the sample 204 by adding said corrections to said guess electrical properties and returning to said means for generating theoretical Fourier transform values 1110, if it is determined by said means for determining whether said corrections are within said acceptable limit 1114 that said corrections are not within said acceptable limit.

17. The computer program product according to claim 11, further comprising:

a means for triggering said fourth computer readable program code means upon the completion of said third computer readable program code means.

18. The computer program product according to claim 11, further comprising:

a means for triggering said fourth computer readable program code means when noise associated with the sample 202 is minimized.

19. The computer program product according to claim 11, wherein said second computer readable program code means and said fourth computer readable program code means detects an instantaneous value of a time dependent net total electric charge.

20. The computer program product according to claim 11, wherein said second computer readable program code means and said fourth computer readable program code means detects a net total electrical current produced by the flow of a time dependent net total electric charge.

21. A method for producing an image indicative of an electrical characteristic of an object, the steps comprising:

a) applying a voltage to the surface of the object with an array of sensor elements that make electrical connection with the surface;

b) measuring the surface charge at each sensor element that results from the applied voltage;

c) transforming the surface charge measurements;

d) calculating from the transformed surface charge measurements and image data input from a tomographic device electrical characteristic values throughout the object; and

e) producing an image of the electrical characteristic of the object from said calculated electrical characteristic values.

22. The method as recited in claim 21 in which the array of sensor elements includes a first set of sensor elements disposed on one side of the object and a second set of sensor elements disposed on the other side of the object, and step a) is performed by applying the voltage across said first and second set of sensor elements.

23. The method as recited in claim 22 in which the surface charge is measured at each sensor element by measuring the current flowing therethrough.

24. The method as recited in claim 22 which includes:

f) repeating steps a), b), and c) without the object being disposed between the first and second set of sensor elements; and

step d) includes;

computing difference transformed values by subtracting the transformed surface charge measurements made with the object disposed between the first and second set of sensor elements and the transformed surface charge measurements made without the object disposed between the first and second set of sensor elements.

25. The method as recited in claim 24 in which step d) further includes:

inverting the difference transformed values to thereby produce inverted transformed values; and

determining scale factors for said inverted transformed values.

26. The method as recited in claim 25 in which the electrical characteristic values are calculated in step d) using a neural network analysis.

27. The method as recited in claim 25 in which the electrical characteristic values are calculated in step d) using a look up table.

28. The method as recited in claim 21 which includes inserting an impedance matching media between the sensor elements and the surface of the object to make electrical connection therebetween.

29. The method as recited in claim 21 in which the image produced in step e) indicates the electrical conductivity throughout the object.

30. The method as recited in claim 21 in which the image produced in step e) indicates the dielectric constant throughout the object.

31. The method as recited in claim 21 in which the image produced in step e) is produced by setting the intensity of each pixel in the image to a level corresponding to the value of the electrical property at the corresponding location in the object.

32. The method as recited in claim 21 in which the transformation performed in step c) is a Fourier transformation.

33. An imaging system for producing an image indicative of an electrical property of an object which comprises:

an array of sensor elements for placement in electrical connection with the surface of the object;

means for applying a predetermined voltage to the array of sensor elements;

means for measuring the electrical charge at each of the sensor elements while the voltage is applied;

means for transforming the measured electrical charges to a set of transformed surface charges;

means responsive to image data input from a tomographic device for calculating electrical characteristic values at locations throughout the object from the set of transformed surface charges; and

means for producing an image depicting the electrical characteristic throughout the object from the calculated electrical characteristic values.

34. The imaging system as recited in claim 33 in which the array of sensor elements includes a top array disposed on one side of the object and a bottom array disposed on the opposite side of the object.

35. The imaging system as recited in claim 34 which includes a matching medium disposed between each of the top and bottom arrays and the object to make electrical connection therebetween.

36. The imaging system as recited in claim 35 in which the impedance of the matching medium is substantially the same as that of the object being imaged.

37. The imaging system as recited in claim 34 in which each of the top and bottom arrays includes a segmented sensor plate, the means for applying a voltage includes means for controlling the voltage amplitude and frequency to each segment, and the means for measuring the electrical charge includes means for detecting the time dependent net total electric charges at each of said segments.

38. The imaging system as recited in claim 37 in which the means for detecting the time dependent net total electric charge includes means for detecting a net total electric current flowing through its associated segment.

39. The imaging system as recited in claim 34 in which the array of sensor elements also includes a pair of spaced side arrays disposed on opposite sides of the object and positioned with respect to the top and bottom arrays to substantially surround the object.

40. The imaging system as recited in claim 33 in which the means for applying a voltage includes a frequency source for applying an alternating voltage at a selected frequency to the array of sensor elements.

41. The imaging system as recited in claim 33 in which the object is a living animal and means for producing a trigger signal responsive to physiological motion in the living animal is coupled to the means for measuring the electrical charge to initiate said measurement.

42. The imaging system as recited in claim 33 in which the array of sensor elements encircles the object.

43. The imaging system as recited in claim 42 in which the array of sensor elements has a fixed geometric shape and a matching medium is disposed between the array of sensor elements and the object to make electric connection therebetween.

44. The imaging system as recited in claim 43 in which the means for measuring the electrical charge includes a plurality of controller units which enable the electrical charge to be measured on all of the sensor elements at substantially the same time.

45. The imaging system as recited in claim 33 in which the means for measuring the electrical charge includes a plurality of controller units which enable the electrical charge to be measured on all of the sensor elements at substantially the same time.