IP Library › Granted Patent US 12,376,812
Granted Patent B2
US 12,376,812 · App. 17/252,682 · Granted Aug 5, 2025

Methods and apparatus for dual energy X-ray imaging

Inventors: Mike Sattarivand (Halifax, CA); Michael Reno (Halifax, CA)
Assignee: DALHOUSIE UNIVERSITY
A61B6/482A61B6/032A61B6/40A61B6/505A61B6/5235G06T5/50G06T5/70G06T2207/10116G06T2207/30061G06T2207/30096
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,376,812
App. No.
17/252,682
Granted
Aug 5, 2025
Kind
B2
Abstract

Dual energy x-ray images may be generated by combining higher and lower energy images of a subject's tissue region with weighting factor values determined based on thicknesses of different tissue types corresponding to the imaged tissue region. Different weighting factor values may be applied to different pixels of the higher or lower energy images. Weighting factor values may, for example, be retrieved from memory, calculated, interpolated and/or extrapolated or the like. The higher and lower energy images and the weighting factor values may be linearly combined to generate a dual energy x-ray image. In some embodiments the weighting factor values may increase contrast of a tissue type.

Claims (55)

1. A method for dual energy x-ray imaging, the method comprising:

acquiring first and second images of a tissue region comprising a first tissue type having a first density and a second tissue type having a second density, the first image corresponding to a first radiation beam having a first energy spectrum and the second image corresponding to a second radiation beam having a second energy spectrum different from the first energy spectrum; and

combining the first and second images to generate a combined image by, separately for each of multiple parts of the first and second images, determining a weighting factor corresponding to the part based on determining a thickness of tissue of each of at least two classes corresponding to the part and pixelwise combining the first and second images using the weighting factors,

wherein each of the parts of the first and second images is one pixel.

2. The method according to claim 1 wherein combining the first and second images comprises computing a weighted difference of logarithms of pixels values in the first and second images by performing the computation:

ln( I DE )=ln( I HE )−ω ln( I LE ),

or a mathematical equivalent thereof, wherein I HE represents pixel values in the first image, I LE represents pixel values in the second image, ω represents the weighting factors and I DE represents pixel values in the combined image.

3. The method according to claim 1 wherein combining the first and second images comprises applying a non-rigid transformation to at least one of the first and second images.

4. The method according to claim 1 wherein determining the thickness of tissue of each of the at least two classes comprises processing volumetric data for each tissue region by ray tracing, the volumetric data obtained by CT scan.

5. The method according to claim 1 wherein determining the thickness of tissue of each of the at least two classes comprises processing CT data to classify voxels of the CT data into the at least two classes based on densities indicated by the CT data for the voxels, generating a volumetric mask for each of the at least two classes, each volumetric mask comprising data indicating for voxels of the volumetric mask whether the voxel does or does not correspond to tissue of the corresponding class, and processing the volumetric masks.

6. The method according to claim 5 wherein processing the volumetric masks comprises computing digitally reconstructed radiographs (DRRs) for each of the volumetric masks from the point of view of the first and second images.

7. The method according to claim 6 wherein processing the volumetric masks is based on a geometry of the first and second radiation beams and comprises determining the thickness of tissue of each of the at least two classes along rays of the first and second radiation beams corresponding to pixels of the first and second images.

8. The method according to claim 1 comprising applying the thickness of tissue of each of the at least two classes as inputs to a calculation to compute the weighting factor for each of the parts.

9. The method according to claim 1 comprising using the thickness of tissue of each of the at least two classes as keys to retrieve the weighting factor for each of the parts from a lookup table.

10. The method according to claim 1 comprising acquiring one or more dark images of the tissue model; and pixelwise subtracting the one or more dark images from each of the first and second images of the tissue model.

11. The method according to claim 1 comprising acquiring at least one of first and second flood images, the flood images respectively corresponding to the first and second radiation beams and pixelwise dividing at least the first image of the tissue model with the first flood image or the second image of the tissue model with the second flood image.

12. The method according to claim 1 comprising:

identifying a region of interest;

wherein determining the weighting factors comprises:

for pixels within the region of interest selecting the weighting factors that cause increased contrast for a first selected tissue type; and

for pixels outside of the region of interest selecting the weighting factors that cause de-emphasis for a second selected tissue type.

13. The method according to claim 12 wherein determining the weighting factors comprises:

for pixels within the region of interest using a first algorithm to determine the weighting factors; and

for pixels outside of the region of interest using a second algorithm different from the first algorithm to determine the weighting factors.

14. The method according to claim 12 wherein determining the weighting factors comprises:

for pixels within the region of interest selecting weighting factors from a first set of weighting factors; and

for pixels outside of the region of interest selecting the weighting factors from a second set of weighting factors.

15. An x-ray imaging method comprising:

acquiring a first image of a tissue region by controlling a radiation source to emit a first radiation beam having a first energy range;

acquiring a second image of the tissue region by controlling the radiation source to emit a second radiation beam having a second energy range wherein a maximum energy of the second energy range is lower than a maximum energy of the first energy range;

classifying tissues having a density greater than a threshold density in a first one of plural density ranges and classifying tissues having density less than the threshold density in a second one of the plural density ranges;

for each pixel of the first and second images, determining amounts of tissue in each of the first and second ones of the plural density ranges lying on a ray extending between the radiation source and a location on a detector used to acquire the pixel of the first and second images and determining a weighting factor for each pixel based on the amounts;

pixelwise combining the first and second images to generate a combined image by, separately for each pixel of the combined image, combining corresponding pixels of the first and second images using the corresponding one of the weighting factors.

16. The method according to claim 15 wherein combining the first and second images comprises computing a weighted difference of logarithms of pixels values in the first and second images by performing the computation:

ln( I DE )=ln( I HE )−ω ln( I LE ),

or a mathematical equivalent thereof, wherein I HE represents pixel values in the first image, I LE represents pixel values in the second image, ω represents the weighting factors and I DE represents pixel values in the combined image.

17. The method according to claim 15 wherein determining the amount of tissue in each of the plural density ranges comprises processing volumetric data for the tissue region.

18. The method according to claim 17 wherein the volumetric data comprises computed tomography (CT) data and determining the amount of tissue in each of the plural density ranges comprises processing the CT data to classify voxels of the CT data into the plural density ranges based on densities indicated by the CT data for the voxels, generating a volumetric mask for each of the plural density ranges, each volumetric mask comprising data indicating for voxels of the volumetric mask whether the voxel does or does not correspond to tissue of the corresponding density range, and processing the volumetric masks.

19. The method according to claim 18 wherein processing the volumetric masks comprises computing digitally reconstructed radiographs (DRRs) for each of the masks from the point of view of the first and second images.

20. The method according to claim 19 wherein processing the volumetric masks is based on a geometry of the first and second radiation beams and comprises determining the thickness of tissue of each of the at least two classes along rays of the first and second radiation beams corresponding to pixels of the first and second images.

21. The method according to claim 15 comprising:

identifying a region of interest;

wherein determining the weighting factors comprises:

for pixels within the region of interest selecting the weighting factors that cause increased contrast for a first selected tissue type; and

for pixels outside of the region of interest selecting the weighting factors that cause de-emphasis for a second selected tissue type.

22. The method according to claim 21 wherein determining the weighting factors comprises:

for pixels within the region of interest using a first algorithm to determine the weighting factors; and

for pixels outside of the region of interest using a second algorithm different from the first algorithm to determine the weighting factors.

23. The method according to claim 21 comprising, for pixels in a margin along a boundary of the region of interest for pixels within the region of interest selecting the weighting factors by blending the weighting factors that cause increased contrast for a first selected tissue type and the weighting factors that cause de-emphasis for a second selected tissue type.

24. A dual-energy x-ray imaging system, the system comprising:

an x-ray imaging device comprising a radiation source and a radiation detector;

an image acquisition module, the image acquisition module configured to generate instructions instructing the x-ray imaging device to acquire a first image of a tissue region using a first radiation beam having a first energy range and to acquire a second image of the tissue region using a second radiation beam having a second energy range, the second energy range lower than the first energy range;

an imaged tissue thickness computation module, the image tissue thickness module configured to compute a thickness of tissue of each of at least two classes corresponding to each of multiple pixels of the first and second images;

a weighting factor selection module, the weighting factor selection module configured to determine, for each pixel, a weighting factor based at least in part on the thicknesses of tissue of each of at least two classes corresponding to the pixel; and

a dual energy image generation module, the dual energy image generation module configured to combine the first and second images to generate a combined image by, separately for each of the multiple pixels of the first and second images, pixelwise combining the first and second images using the determined weighting factor.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Dec 15, 2020
From: SATTARIVAND, MIKE; RENO, MICHAEL
To: DALHOUSIE UNIVERSITY
Reel/Frame 054658/0942 →
Continuity (2)
Provisional Application 62685776 · Jun 15, 2018
Related Publication 20210267563A1 · Sep 2, 2021
References Cited (82)
US 4792900A · Sones et al. · 1988 [cited by applicant]
US 5841833A · Mazess et al. · 1998 [cited by applicant]
US 6661873B2 · Jabri et al. · 2003 [cited by applicant]
US 6683934B1 · Zhao et al. · 2004 [cited by applicant]
US 6816564B2 · Charles, Jr. et al. · 2004 [cited by applicant]
US 7010092B2 · Winsor · 2006 [cited by applicant]
US 7068826B2 · Jabri et al. · 2006 [cited by applicant]
US 7627084B2 · Jabri et al. · 2009 [cited by applicant]
US 7873141B2 · Imai et al. · 2011 [cited by applicant]
US 9351695B2 · Wang et al. · 2016 [cited by applicant]
US 9492134B2 · Takasaki · 2016 [cited by applicant]
US 9662078B2 · Berglund et al. · 2017 [cited by applicant]
US 9693742B2 · Grasruck et al. · 2017 [cited by applicant]
US 9754387B2 · Leng et al. · 2017 [cited by applicant]
US 9820712B2 · Takasaki · 2017 [cited by applicant]
US 9886765B2 · Naito · 2018 [cited by applicant]
US 9905003B2 · Maack · 2018 [cited by examiner]
US 20120082295A1 · Agano et al. · 2012 [cited by applicant]
US 20140219423A1 · Bertens · 2014 [cited by applicant]
US 20170116730A1 · Yamanaka · 2017 [cited by applicant]
US 20170000439A1 · Takasaki · 2017 [cited by applicant]
US 20180015306A1 · Maurer, Jr. et al. · 2018 [cited by applicant]
US 20180061088A1 · Bhagalia et al. · 2018 [cited by applicant]
US 20180068442A1 · Kawamura · 2018 [cited by examiner]
DE 102011004120B4 · 2017 [cited by applicant]
EP 2548510B1 · 2016 [cited by applicant]
EP 2243020B1 · 2017 [cited by applicant]
EP 3287774A1 · 2018 [cited by applicant]
JP 05236351A · 1993 [cited by applicant]
JP 0785247A · 1995 [cited by applicant]
JP 2002152593A · 2002 [cited by applicant]
JP 2007044520A · 2007 [cited by applicant]
JP 4526110B2 · 2010 [cited by applicant]
JP 2010194261A · 2010 [cited by applicant]
JP 201385967A · 2013 [cited by applicant]
JP 5680718B2 · 2015 [cited by applicant]
JP 6129125B2 · 2017 [cited by applicant]
JP 2017122705A · 2017 [cited by applicant]
JP 6242631B2 · 2017 [cited by applicant]
JP 6261915B2 · 2018 [cited by applicant]
JP 2018015664A · 2018 [cited by applicant]
JP 6275374B2 · 2018 [cited by applicant]
JP 2018506385A · 2018 [cited by applicant]
KR 101848177B1 · 2018 [cited by applicant]
WO 2013037659A1 · 2013 [cited by applicant]
WO 2016168194A1 · 2016 [cited by applicant]
WO 2018005721A1 · 2018 [cited by applicant]
Seidlitz, J., Sponheim, C., Glen, D., Frank, Q.Y., Saleem, K.S., Leopold, D.A., Ungerleider, L. and Messinger, A., 2018. A population MRI brain template and analysis tools for the macaque. Neuroimage, 170, pp. 121-131. [cited by examiner]
Wan Y, Hu H, Xu Y, Chen Q, Wang Y, Gao D. A Robust and Accurate Non-rigid Medical Image Registration Algorithm Based on Multi-level Deformable Model. Iran J Public Health. Dec. 2017;46(12):1679-1689. PMID: 29259943; PMC… [cited by examiner]
Alvarez, R. E. et al., “Energy-Selective Reconstructions in X-Ray Computerized Tomography”, Physics in Medicine and Biology 21(5):733-744, 1976. [cited by applicant]
Ay, M. R. et al., “Assessment of Different Computational Models for Generation of X-Ray Spectra in Diagnostic Radiology and Mommography”, Med Phys 32:1660-1675, 2006. [cited by applicant]
Bowman, W. A. et al., “Optimizing Dual-Energy X-Ray Parameters for the ExacTrac Clinical Stereoscopic Imaging System to Enhance Soft-Tissue Imaging”, Medical Physics 44(3):823-831, doi:10.1002/mp.12093, 2017. [cited by applicant]
Burton, C. S. et al., “Energy Subtraction Angiography is Comparable to Digital Subtraction Angiography in Terms of Iodine Rose SNR”, Medical Physics 43(11):5925, doi:10.1118/1.4962651, 2016. [cited by applicant]
Dhont, J., D. et al., “Feasibility of Markerless Tumor Tracking by Sequential Dual-Energy Fluoroscopy on a Clinical Tumor Tracking System”, Radiotherapy and Oncology: Journal of the European Society for Therapeutic Radi… [cited by applicant]
Ducote, J. L. et al., “Optimization of a Flat-Panel Based Real Time Dual-Energy System for Cardiac Imaging”, Medical Physics 33(6):1562-1568, doi:10.1118/1.2174131, 2006. [cited by applicant]
Ergun, D. L. et al., “Single-Exposure Dual-Energy Computed Radiography: Improved Detection and Processing”, Radiology 174(1):243-249, doi:10.1148/radiology.174.1.2294555, 1990. [cited by applicant]
Fischbach, F. et al., “Dual-Energy Chest Radiography with a Flat-Panel Digital Detector: Revealing Calcified Chest Abnormalities”, AJR, AM. J. Roentgenol. 181:1519-1524, 2003. [cited by applicant]
Ho, J. T. et al., “Comparison of Dual and Single Exposure Techniques in Dual-Energy Chest Radiography”, Medical Physics 16(2):202-208, doi:10.1118/1.596372, 1989. [cited by applicant]
Hoggarth, M. A. et al., “Dual Energy Imaging using a Clinical on-Board Imaging System”, Physics in Medicine and Biology 58(12):4331-4340, doi:10.1088/0031-9155/58/12/4331, 2013. [cited by applicant]
Kalender, W. A. et al., “An Algorithm for Noise Suppression in Dual Energy CT Material Density Images”, IEEE Transactions on Medical Imaging 7(3):218-224, doi:10.1109/42.7785, 1988. [cited by applicant]
Karunamuni, R. et al., “Search for Novel Contrast Materials in Dual-Energy X-Ray Breast Imaging using Theoretical Modeling of Contrast-to-Noise Ratio”, Physics in Medicine and Biology 59(15):4311-4324, doi:10.1088/0031-… [cited by applicant]
Lehmann, L. A. et al., “Generalized Image Combinations in Dual KVP Digital Radiography”, Medical Physics 8:659-667, 1981. [cited by applicant]
Li, H. et al., “Implementation of Dual-Energy Technique for Virtual Monochromatic and Linearly Mixed CBCTs”, Medical Physics 39(10):6056-6064, doi:10.1118/1.4752212, 2012. [cited by applicant]
McCollough, C. H. et al., “A Correlated Noise Reduction Algorithm for Dual-Energy Digital Subtraction Angiography”, Medical Physics 16(6):873-880, doi:10.1118/1.596436, 1989. [cited by applicant]
Menten, M. J. et al., “Using Dual-Energy X-Ray Imaging to Enhance Automated Lung Tumor Tracking during Real-Time Adaptive Radiotherapy”, Medical Physics 42(12):6987-6998, doi:10.1118/1.4935431, 2015. [cited by applicant]
Molloi, S. Y. et al., “Quantification Techniques for Dual-Energy Cardiac Imaging”, Medical Physics 16(2):209-217, doi:10.1118/1.596418, 1989. [cited by applicant]
Rebuffel, V. et al., “Dual-Energy X-Ray Imaging: Benefits and Limits”, Ecndt 49:1-12, 2006. [cited by applicant]
Richard, S. et al., “Cascaded Systems Analysis of Noise Reduction Algorithms in Dual-Energy Imaging”, Medical Physics 35(2):586-601, doi:10.1118/1.2826556, 2008. [cited by applicant]
Richard, S. et al., “Optimization of Dual-Energy Imaging Systems using Generalized NEQ and Imaging Task”, Medical Physics 34(1):127-139, doi:10.1118/1.2400620, 2007. [cited by applicant]
Sabol, J. et al., “Simulated and Experimental Technique Optimization of Dual-Energy Radiography: Abdominal Imaging Applications”, Medical Imaging 6142:545-556, 2006. [cited by applicant]
Saito, M. et al., “A Simple Formulation for Deriving Effective Atomic Numbers Via Electron Density Calibration from Dual-Energy CT Data in the Human Body”, Medical Physics (44):2293-2303, 2017. [cited by applicant]
Shkumat, N. A. et al., “Optimization of Image Acquisition Techniques for Dual-Energy Imaging of the Chest”, Medical Physics 34(10):3904-3915, doi:10.1118/1.2777278, 2007. [cited by applicant]
Siddon, R. L. “Fast Calculation of the Exact Radiological Path for a Three-Dimensional CT Array”, Medical Physics 12(2):252-255, doi:10.1118/1.595715, 1985. [cited by applicant]
Szucs-Farkas, Z. et al., “Single-Exposure Dual-Energy Subtraction Chest Radiography: Detection of Pulmonary Nodules and Masses in Clinical Practice”, European Radiology 18(1):24-31, doi:10.1007/s00330-007-0758-z, 2008. [cited by applicant]
Vock, P. et al., “Dual Energy Subtraction: Principles and Clinical Applications”, European Journal of Radiology 72(2):231-237, doi:10.1016/j.ejrad.2009.03.046, 2009. [cited by applicant]
Williams, D. B. et al., “Optimal kVp Selection for Dual-Energy Imaging of the Chest: Evaluation by Task-Specific Observer Preference Tests”, Medical Physics 34(10):3916-3925, doi:10.1118/1.2776239, 2007. [cited by applicant]
Xu, T. et al., “Dynamic Dual-Energy Chest Radiography: A Potential Tool for Lung Tissue Motion Monitoring and Kinetic Study”, Physics in Medicine and Biology 56(4):1191-1205, doi:10.1088/0031-9155/56/4/019, 2011. [cited by applicant]
Yu, L. et al., “Virtual Monochromatic Imaging in Dual-Source Dual-Energy CT: Radiation Dose and Image Quality”, Medical Physics 38:6371-6379, 2011. [cited by applicant]
“Dual Energy Radiography Acquisition and Processing”, SUNY Upstate Medical University, http://www.upstate.edu/radiology/education/rsna/radiography/dual.php, Mar. 14, 2018. [cited by applicant]
Schmidt, T., “Optimal ‘image-based’ weighting for energy-resolved CT”, Med. Phys. 36(7):3018-27, doi:10.1118/1.3148535, 2009. [cited by applicant]
Cruz-Bastida, J.-P. et al., Contrast Optimization in Clinical Contrast-Enhanced Digital Mammography Images. In: Maidment A.D.A., Bakic P.R., Gavenonis S. (eds) Breast Imaging. IWDM 2012. Lecture Notes in Computer Scienc… [cited by applicant]
Vijayakumar, S. et al., “Beam's Eye View-based Radiation Therapy: Description of Methods”, RadioGraphics 1992; 12:961-968. [cited by applicant]