Methods for acquiring multi spectral data of an object

• US 7,116,749 B2
• Filed: 05/19/2005
• Issued: 10/03/2006
• Est. Priority Date: 06/25/2003
• Status: Expired due to Fees

• Current Assignee: Forevision Imaging Technologies LLC
• Original Assignee: Guy M. Besson
• Inventors: Besson, Guy M.
• Primary Examiner(s): BRUCE, DAVID VERNON
• Application Number:
  US11/133,114
• Publication Number:
  US 20050220265A1
• Time in Patent Office: 502 Days
• Field of Search: 378/4, 378/8, 378/15, 378/16, 378/98.9, 378/901
• US Class Current: 378/16
• CPC Class Codes: A61B 6/027 characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
• CPC Class Codes: A61B 6/032 Transmission computed tomography [CT]
• CPC Class Codes: A61B 6/0487 Motor-assisted positioning
• CPC Class Codes: A61B 6/06 Diaphragms
• CPC Class Codes: A61B 6/4028 resulting in acquisition of views from substantially different positions, e.g. EBCT
• CPC Class Codes: A61B 6/4035 the source being combined with a filter or grating
• CPC Class Codes: A61B 6/4042 K-edge filters
• CPC Class Codes: A61B 6/405 Source units specially adapted to modify characteristics of the beam during the data acquisition process A61B6/4021, A61B6/4035 take precedence
• CPC Class Codes: A61B 6/4085 Cone-beams
• CPC Class Codes: A61B 6/4241 using energy resolving detectors, e.g. photon counting
• CPC Class Codes: A61B 6/4275 using a detector unit almost surrounding the patient, e.g. more than 180°
• CPC Class Codes: A61B 6/4441 the rigid structure being a C-arm or U-arm
• CPC Class Codes: A61B 6/4488 Means for cooling A61B6/045 takes precedence
• CPC Class Codes: A61B 6/482 involving multiple energy imaging
• CPC Class Codes: A61B 6/484 involving phase contrast X-ray imaging
• CPC Class Codes: A61B 6/488 involving pre-scan acquisition
• CPC Class Codes: A61B 6/508 for non-human patients
• CPC Class Codes: A61B 6/5288 involving retrospective matching to a physiological signal
• CPC Class Codes: A61B 6/541 involving acquisition triggered by a physiological signal
• CPC Class Codes: A61B 6/542 involving control of exposure
• CPC Class Codes: A61B 6/563 : involving image data transmission via a network
• CPC Class Codes: A61B 6/583 : using calibration phantoms
• CPC Class Codes: G21K 1/10 : Scattering devices; Absorbing devices; Ionising radiation filters
• CPC Class Codes: Y10S 378/901 : Computer tomography program or processor

Abstract

A multispectral X-ray imaging system uses a wideband source and filtration assembly to select for M sets of spectral data. Spectral characteristics may be dynamically adjusted in synchrony with scan excursions where an X-ray source, detector array, or body may be moved relative to one another in acquiring T sets of measurement data. The system may be used in projection imaging and/or CT imaging. Processed image data, such as a CT reconstructed image, may be decomposed onto basis functions for analytical processing of multispectral image data to facilitate computer assisted diagnostics. The system may perform this diagnostic function in medical applications and/or security applications.

Claims

1. A method of acquiring multispectral x-ray projection data of an object with an imaging apparatus, comprising:

(a) Selecting an x-ray source x-ray radiation spectral characteristics by adjusting at least two parameter values from the group consisting of (i) an x-ray tube target track material, (ii) a generator peak kilo-voltage, (iii) an x-ray tube current, (iv) a filter composition, and (v) an x-ray tube focal spot size;

(b) projecting, in a first relative orientation between the object and x-ray imaging apparatus, x-ray radiation according to selecting step (a) through the object and acquiring a first set of measurements defining a first set of projection lines through the object;

(c) for at least one additional relative orientation between the object and the x-ray imaging apparatus, projecting x-ray radiation according to selecting step (a) through the object and acquiring an additional set of measurements defining an additional set of projection lines through the object, the additional set of projection lines overlapping, at least in part, with the first set of projection lines; and

(d) repeating steps (a) through (c) for at least two different parameter values;

wherein multispectral sets of projection data are acquired for the object at a multiplicity of x-ray spectra in a single scan for the object through the imaging apparatus.

2. The method of claim 1, further comprising filtering projected x-ray radiation through the object using a double filter wheel.

3. The method of claim 2, further comprising independently controlling each filter of the double filter wheel.

4. The method of claim 1, further comprising digitally processing the acquired projection data to increase projection data signal-to-noise ratio.

5. The method of claim 4, the digital processing comprising sorting and averaging the acquired projection data.

6. The method of claim 4, the digital processing comprising utilizing statistical methods.

7. The method of claim 1, wherein selecting step (a) further comprises analyzing data collected during the scan and selecting the parameter values as a function of the object x-ray attenuation currently present along the x-ray projection lines.

8. The method of claim 1, further comprising setting the position of an x-ray beam with respect to a detector based upon the acquired projection data.

9. The method of claim 8, wherein setting the position of the x-ray beam comprises determining current position of the x-ray beam with respect to the detector.

10. The method of claim 1, further comprising recording signal information for a patient electrocardiogram or respiratory trace simultaneously to x-ray data acquisition during the scan.

11. The method of claim 10, further comprising sorting acquired x-ray projection data according to the recorded electrocardiogram or respiratory signal information.

12. The method of claim 1, further comprising measuring air signal and performing calibration and correction steps to provide estimates of line integrals of the object being imaged.

13. The method of claim 1, wherein projecting step (b) comprises illuminating part of a detector directly by a primary x-ray beam projection and measuring scattered radiation with the detector outside the primary beam projection.

14. The method of claim 13, further comprising estimating a scatter signal collected together with a primary signal from the scatter data measured outside the primary beam projection.

15. The method of claim 14, further comprising subtracting the estimated scatter signal from the recorded primary and scatter data to obtain an estimate of the primary signal.

16. The method of claim 1, wherein step (b) of projecting x-ray radiation comprises operating the x-ray source in a pulsed mode or a continuous mode.

17. A method of acquiring multispectral computed tomography data of an object with a computed tomography apparatus, comprising:

(a) selecting an x-ray source x-ray radiation spectral characteristics by adjusting at least two parameter values from the group consisting of (i) an x-ray tube target track material, (ii) a generator peak kilo-voltage, (iii) an x-ray tube current, (iv) a filter composition, and (v) an x-ray tube focal spot size;

(b) projecting, for at least one relative orientation between the x-ray source and the object, x-ray radiation according to the selecting step (a) through the object and acquiring a first set of measurements defining a first set of projection lines through the object;

(c) repeating steps (a) and (b) in a different relative orientation between the x-ray source and the object, such that at least two sets of projection lines through the object are acquired at at least two spectra;

wherein projection data sets acquired in a single scan for the object through the computed tomography apparatus are sufficient for the generation of at least two computed tomography image estimates at at least two spectra for at least part of one cross section of the object.

18. The method of claim 17, further comprising filtering projected x-ray radiation through the object using a double filter wheel.

19. The method of claim 18, further comprising independently controlling each filter of the double filter wheel.

20. The method of claim 17, comprising acquiring the data for each projection line spectral data set corresponding to one object cross section with a source excursion greater than 720 degrees.

21. The method of claim 20, further comprising weighting the acquired data prior to reconstruction to account for data redundancy.

22. The method of claim 17, comprising acquiring the data for each projection line spectral data set corresponding to one object cross section with a source excursion greater than 360 degrees.

23. The method of claim 22, further comprising weighting the acquired data prior to reconstruction to account for data redundancy.

24. The method of claim 17, comprising acquiring the data for each projection line spectral data set corresponding to one object cross section with a source excursion less than 360 degrees.

25. The method of claim 24, further comprising weighting the acquired data prior to reconstruction to account for data redundancy.

26. The method of claim 17, wherein projecting step (b) comprises illuminating part of a detector directly by a primary x-ray beam projection and measuring scattered radiation with the detector outside the primary beam projection.

27. The method of claim 26, further comprising estimating a scatter signal collected together with a primary signal from the scatter data measured outside the primary beam projection.

28. The method of claim 27, further comprising subtracting the estimated scatter signal from the recorded primary and scatter data and obtaining an estimate of the primary signal.

29. The method of claim 17, comprising changing the x-ray spectral parameters from one view to the next such that at least two sets of projection data each sufficient for image reconstruction are acquired for at least part of one object cross section for each source rotation of at least 360 degrees.

30. The method of claim 17, wherein selecting step (a) comprises analyzing data collected during the scan and selecting the parameter values as a function of the object x-ray attenuation currently present along the x-ray projection lines.

31. The method of claim 17, further comprising recording signal information for a patient electrocardiogram or respiratory trace simultaneously to x-ray data acquisition during the scan.

32. The method of claim 31, further comprising sorting the acquired x-ray projection data according to the recorded electrocardiogram or respiratory signal information.

33. The method of claim 17, wherein step (b) of projecting x-ray radiation comprises operating the x-ray source in a pulsed mode or a continuous mode.

34. A method of acquiring multispectral limited-angle computed tomography or tomosynthesis data of an object with an x-ray imaging system, comprising:

(a) selecting an x-ray source x-ray radiation spectral characteristics by adjusting at least two parameter values for the group consisting of (i) an x-ray tube target track material, (ii) a generator peak kilo-voltage, (iii) an x-ray tube current, (iv) a filter composition, and (v) an x-ray tube focal spot size;

(b) projecting, for at least one relative orientation of the x-ray source with respect to the object, x-ray radiation according to the selecting step (a) through the object and acquiring a first set of measurements defining a first set of projection lines through the object;

(c) repeating steps (a) and (b) in a different relative orientation of the x-ray source with respect to the object, such that at least two sets of projection lines through the object are acquired at at least two spectra;

wherein projection data sets acquired of the object in a single scan for the object through the x-ray imaging system are sufficient for the generation of at least two limited-angle computed tomography or tomosynthesis image estimates at at least two spectra for at least part of one cross section of the object.

35. The method of claim 34, wherein step (c) of repeating steps (a) and (b) in a different relative orientation of the x-ray source with respect to the object comprises moving an x-ray tube along a trajectory.

36. The method of claim 34, wherein step (c) of repeating steps (a) and (b) in a different relative orientation of the x-ray source with respect to the object comprises deflecting a focused electron beam along an x-ray target.

37. The method of claim 34, further comprising filtering projected x-ray radiation through the object by mechanically advancing at least one x-ray filter sheet in front of the source.

38. The method of claim 34, further comprising digitally processing the acquired projection data to increase the projection data signal-to-noise ratio.

39. The method of claim 38, the digital processing comprising sorting and averaging the data.

40. The method of claim 38, the digital processing comprising using statistical methods.

41. The method of claim 34, wherein selecting step (a) comprises analyzing data collected during the scan and selecting the parameter values as a function of the object x-ray attenuation currently present along the x-ray projection lines.

42. The method of claim 34, further comprising setting the position of an x-ray beam with respect to a detector based upon the acquired projection data.

43. The method of claim 42, wherein setting the position of the x-ray beam comprises determining current position of the x-ray beam with respect to the detector.

44. The method of claim 34, further comprising recording signal information for a patient electrocardiogram or respiratory trace simultaneously to x-ray data acquisition during the scan.

45. The method of claim 44, further comprising sorting acquired x-ray projection data according to the recorded electrocardiogram or respiratory signal information.

46. The method of claim 34, further comprising measuring air signal and performing calibration and correction steps to provide estimates of line integrals of the object being imaged.

47. The method of claim 34, wherein projecting step (b) comprises illuminating part of a detector directly by a primary x-ray beam projection and measuring with the detector scattered radiation outside the primary beam projection.

48. The method of claim 47, further comprising estimating a scatter signal collected together with a primary signal from the scatter data measured outside the primary beam projection.

49. The method of claim 48, further comprising subtracting the estimated scatter signal from the recorded primary and scatter data to obtain an estimate of the primary signal.

50. The method of claim 34, wherein step (b) of projecting x-ray radiation comprises operating the x-ray source in a pulsed mode or a continuous mode.

51. The method of claim 34, further comprising weighting the acquired projection data prior to reconstruction to account for the geometry of the x-ray acquisition system.

52. The method of claim 34, wherein the generation of limited-angle computed tomography or tomosynthesis reconstructed images comprises performing a filtered-backprojection method.

53. The method of claim 34, wherein the generation of limited-angle computed tomography or tomosynthesis reconstructed images comprises performing a matrix inversion method.

54. The method of claim 34, wherein the generation of limited-angle computed tomography or tomosynthesis reconstructed images comprises performing an iterative method.

55. The method of claim 54, wherein performing an iterative reconstruction method further comprises using a-priori information.