US20140236488A1 - Computed Tomography Calibration Systems and Methods - Google Patents
Computed Tomography Calibration Systems and Methods Download PDFInfo
- Publication number
- US20140236488A1 US20140236488A1 US13/770,617 US201313770617A US2014236488A1 US 20140236488 A1 US20140236488 A1 US 20140236488A1 US 201313770617 A US201313770617 A US 201313770617A US 2014236488 A1 US2014236488 A1 US 2014236488A1
- Authority
- US
- United States
- Prior art keywords
- tissue
- basis material
- basis
- unknown
- pixel value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims description 70
- 238000002591 computed tomography Methods 0.000 title description 173
- 239000000463 material Substances 0.000 claims abstract description 346
- 239000012925 reference material Substances 0.000 claims description 114
- 239000000203 mixture Substances 0.000 claims description 82
- 210000001519 tissue Anatomy 0.000 claims description 67
- 210000000988 bone and bone Anatomy 0.000 claims description 64
- 210000001185 bone marrow Anatomy 0.000 claims description 41
- 230000004044 response Effects 0.000 claims description 39
- 238000003384 imaging method Methods 0.000 claims description 34
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 30
- 208000000913 Kidney Calculi Diseases 0.000 claims description 29
- 206010029148 Nephrolithiasis Diseases 0.000 claims description 29
- 238000012512 characterization method Methods 0.000 claims description 23
- 210000005228 liver tissue Anatomy 0.000 claims description 16
- 210000003205 muscle Anatomy 0.000 claims description 16
- 210000005013 brain tissue Anatomy 0.000 claims description 15
- 210000005003 heart tissue Anatomy 0.000 claims description 15
- 229910052742 iron Inorganic materials 0.000 claims description 15
- 210000004072 lung Anatomy 0.000 claims description 15
- 210000004923 pancreatic tissue Anatomy 0.000 claims description 15
- 210000005084 renal tissue Anatomy 0.000 claims description 15
- 210000004872 soft tissue Anatomy 0.000 claims description 15
- 230000001419 dependent effect Effects 0.000 claims description 13
- 238000013507 mapping Methods 0.000 claims description 11
- 230000001338 necrotic effect Effects 0.000 claims description 9
- 230000009471 action Effects 0.000 claims description 8
- 210000004884 grey matter Anatomy 0.000 claims description 8
- 210000004885 white matter Anatomy 0.000 claims description 8
- 239000000126 substance Substances 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 abstract description 16
- 230000009977 dual effect Effects 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000005259 measurement Methods 0.000 description 61
- 229910000396 dipotassium phosphate Inorganic materials 0.000 description 34
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 33
- 229910052500 inorganic mineral Inorganic materials 0.000 description 31
- 239000011707 mineral Substances 0.000 description 31
- ZPWVASYFFYYZEW-UHFFFAOYSA-L dipotassium hydrogen phosphate Chemical compound [K+].[K+].OP([O-])([O-])=O ZPWVASYFFYYZEW-UHFFFAOYSA-L 0.000 description 28
- 230000008569 process Effects 0.000 description 28
- 230000006870 function Effects 0.000 description 19
- 230000008859 change Effects 0.000 description 14
- 229940116269 uric acid Drugs 0.000 description 13
- LEHOTFFKMJEONL-UHFFFAOYSA-N Uric Acid Chemical compound N1C(=O)NC(=O)C2=C1NC(=O)N2 LEHOTFFKMJEONL-UHFFFAOYSA-N 0.000 description 12
- TVWHNULVHGKJHS-UHFFFAOYSA-N Uric acid Natural products N1C(=O)NC(=O)C2NC(=O)NC21 TVWHNULVHGKJHS-UHFFFAOYSA-N 0.000 description 12
- 238000013459 approach Methods 0.000 description 11
- 239000011575 calcium Substances 0.000 description 10
- 210000004003 subcutaneous fat Anatomy 0.000 description 10
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 9
- 229910052791 calcium Inorganic materials 0.000 description 9
- 210000004185 liver Anatomy 0.000 description 8
- XYJRXVWERLGGKC-UHFFFAOYSA-D pentacalcium;hydroxide;triphosphate Chemical compound [OH-].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O XYJRXVWERLGGKC-UHFFFAOYSA-D 0.000 description 8
- 229960003067 cystine Drugs 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000011002 quantification Methods 0.000 description 6
- 102000008186 Collagen Human genes 0.000 description 5
- 108010035532 Collagen Proteins 0.000 description 5
- LEVWYRKDKASIDU-IMJSIDKUSA-N L-cystine Chemical compound [O-]C(=O)[C@@H]([NH3+])CSSC[C@H]([NH3+])C([O-])=O LEVWYRKDKASIDU-IMJSIDKUSA-N 0.000 description 5
- 229920001436 collagen Polymers 0.000 description 5
- 239000000470 constituent Substances 0.000 description 5
- 235000019797 dipotassium phosphate Nutrition 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 239000004575 stone Substances 0.000 description 5
- 239000007864 aqueous solution Substances 0.000 description 4
- 230000037182 bone density Effects 0.000 description 4
- QXDMQSPYEZFLGF-UHFFFAOYSA-L calcium oxalate Chemical compound [Ca+2].[O-]C(=O)C([O-])=O QXDMQSPYEZFLGF-UHFFFAOYSA-L 0.000 description 4
- 238000001727 in vivo Methods 0.000 description 4
- 241000894007 species Species 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000009547 dual-energy X-ray absorptiometry Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- LEVWYRKDKASIDU-QWWZWVQMSA-N D-cystine Chemical compound OC(=O)[C@H](N)CSSC[C@@H](N)C(O)=O LEVWYRKDKASIDU-QWWZWVQMSA-N 0.000 description 2
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 2
- 238000005284 basis set Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000013170 computed tomography imaging Methods 0.000 description 2
- 238000011109 contamination Methods 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 238000000326 densiometry Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 239000008241 heterogeneous mixture Substances 0.000 description 2
- 230000003278 mimic effect Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 229910052717 sulfur Inorganic materials 0.000 description 2
- 239000011593 sulfur Substances 0.000 description 2
- 201000005569 Gout Diseases 0.000 description 1
- 102000001554 Hemoglobins Human genes 0.000 description 1
- 108010054147 Hemoglobins Proteins 0.000 description 1
- 208000008558 Osteophyte Diseases 0.000 description 1
- 208000001132 Osteoporosis Diseases 0.000 description 1
- MUBZPKHOEPUJKR-UHFFFAOYSA-N Oxalic acid Chemical compound OC(=O)C(O)=O MUBZPKHOEPUJKR-UHFFFAOYSA-N 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 206010003246 arthritis Diseases 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 159000000007 calcium salts Chemical class 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- 230000001143 conditioned effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- XUJNEKJLAYXESH-UHFFFAOYSA-N cysteine Natural products SCC(N)C(O)=O XUJNEKJLAYXESH-UHFFFAOYSA-N 0.000 description 1
- 235000018417 cysteine Nutrition 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 229910052588 hydroxylapatite Inorganic materials 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000004962 physiological condition Effects 0.000 description 1
- 238000011045 prefiltration Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 210000003813 thumb Anatomy 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/505—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of bone
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
- G01N23/20075—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials by measuring interferences of X-rays, e.g. Borrmann effect
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4241—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/483—Diagnostic techniques involving scattered radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/419—Imaging computed tomograph
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/423—Imaging multispectral imaging-multiple energy imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/601—Specific applications or type of materials density profile
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/612—Specific applications or type of materials biological material
- G01N2223/6123—Specific applications or type of materials biological material bone mineral
Definitions
- This disclosure relates generally to computer hardware and software for calibrating single or multiple energy computed tomography (CT) devices.
- the hardware and software runs, displays, provides and otherwise uses electronic content; more particularly, such hardware and software operate on certain content that is independent of the CT scanner being used, in order (among other things) to calibrate the scanner and to present one dimensional information the CT scanner generates (such as representative pixel values) in two dimensional basis material space.
- QCT Quantitative Computed Tomography
- DEQCT Dual energy QCT
- DXA Dual energy x-ray absorptiometry
- the diagnostic x-ray energy regime is also interesting in that theoretical models and tabulated x-ray data can be used to predict mass attenuation coefficients for tissues of biological interest with accuracies on the order of 1% to 2%. Energy-dependent mass attenuation coefficients can be used to derive basis material decompositions that are analogous to those measured with dual-energy x-ray methods.
- computer hardware and software can perform single or dual-energy decomposition calculations and predict CT values for either single-energy or dual energy CT from basis material decomposition estimates derived using theoretical models. Using that information, the computer hardware and software can present one dimensional information from a CT scanner, such as representative pixel value for an unknown material, in two dimensions, such as a line in two dimensional basis material space. In some embodiments, for example, the hardware and software can predict CT values for heterogeneous mixtures of multiple tissue types, such as a model of bone including collagen, bone mineral described as calcium hydroxyapatite, and a mixture of red and yellow marrow, based on atomic compositions for the tissue types.
- tissue types such as a model of bone including collagen, bone mineral described as calcium hydroxyapatite, and a mixture of red and yellow marrow, based on atomic compositions for the tissue types.
- the hardware and software can interpret and present observed CT values from a single energy CT scan in a basis material space relative to the model tissue constituents; for example, a representative pixel value acquired from a CT scan can be presented in two dimensional basis material space with attendant new benefits for CT calibration. As another example, pixel values from a dual energy CT scan can be presented as lines in basis material space, whose intersection gives information about the basis material composition of the scanned material.
- a method is executed by a computing device that includes these steps:
- Characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets; 4.
- Some embodiments relate to a system that includes a processor device and a tangible medium embodiment code.
- the code can be executable by the processor to cause the system to perform certain actions.
- the actions can include:
- Characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets; 4.
- Some embodiments relate to a tangible computer readable medium that includes program code that is executable by a processor to perform actions.
- the program code includes code for causing the processor to carry out the following method:
- Characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets; 4.
- the computer devices, processor devices, tangible medium embodiment code, and tangible computer readable medium that includes program code can perform any or all of the following actions:
- CT pixel values to present an unknown material line in two dimensional basis material space, in order to characterize the unknown material.
- Receiving imaging information carrying out methods described above, and presenting information for the unknown material at a first energy level and a second energy level.
- Receiving imaging information for the unknown material can include receiving imaging information for a tissue, mixture, or substance.
- Characterizing unknown material selected from the group consisting of: pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue.
- tissue selected from the group consisting of: fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
- steps 3 in the embodiments stated above can occur at a time (and/or on a CT scanner) different than step 4 in those embodiments; thus, some or all of step 3, that is, characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets can occur at a time (or on a CT scanner) different from some or all of step 4, that is from imaging information related to an unknown material received from the CT scanner at the first energy level; extracting at least one unknown material CT pixel value representative of a pre
- FIG. 1 is a functional block diagram of a system that is capable of carrying out methods according to certain embodiments of the invention.
- FIG. 2 is a flow diagram showing process flow for one method according to an embodiment of the invention.
- FIG. 3 is a diagram that shows a first presentation according to an embodiment of the invention.
- FIG. 4 is a diagram that shows a second presentation according to an embodiment of the invention.
- FIG. 5 is a diagram that shows a third presentation according to an embodiment of the invention.
- FIG. 6 is a diagram that shows a fourth presentation according to an embodiment of the invention.
- FIG. 7 is a diagram that shows a fifth presentation according to an embodiment of the invention.
- FIG. 8 is a graph showing pairs of MRS and DEQCT fat measurements.
- FIG. 9 is a graph showing changes in a fat signal as measured by both MRS and DEQCT methods.
- FIG. 1 depicts a system that is configured, according to certain embodiments of the invention, to do any or all of at least the following: Perform single or dual-energy decomposition calculations and predict CT values for either single-energy or dual energy CT from basis material decomposition estimates derived using theoretical models; present one dimensional information from a CT scanner, such as representative pixel value for an unknown material, in two dimensions, such as a line in two dimensional basis material space; predict CT values for heterogeneous mixtures of multiple tissue types, such as a model of bone including collagen, bone mineral described as calcium hydroxyapatite, and a mixture of red and yellow marrow, based on atomic compositions for the tissue types; present observed CT values from a single energy CT scan in a basis material space relative to the model tissue constituents, such as, for example, presenting a representative pixel value acquired from a CT scan in two dimensional basis material space; present, based on pixel values from a dual energy CT scan, lines in basis material space, whose intersection gives information about the basis material composition of
- the system generally includes a computing device 102 having a processor 104 that can execute code stored on a tangible computer-readable medium, such as a memory 106 , to cause the computing device 102 to perform processes as described more fully herein.
- the computing device 102 may be any device that can process data and execute code that is a set of instructions to perform actions. Examples of the computing device 102 include a database server, a web server, desktop personal computer, a laptop personal computer, a server device, a handheld computing device, and a mobile device.
- Examples of the processor 104 include a microprocessor, an application-specific integrated circuit (ASIC), a state machine, or other suitable processor.
- the processor 104 may include one processor or any number of processors.
- the processor 104 can access code stored in the memory 106 via a bus 108 .
- the memory 106 may be any non-transitory computer-readable medium configured for tangibly embodying code and can include electronic, magnetic, or optical devices. Examples of the memory 106 include random access memory (RAM), read-only memory (ROM), a floppy disk, compact disc, digital video device, magnetic disk, an ASIC, a configured processor, flash drive, thumb drive, cloud memory capacity, or other storage device or capacity.
- the bus 108 may be any device capable of transferring data between components of the computing device 102 .
- the bus 108 can include one device or multiple devices.
- Instructions can be stored in the memory 106 as executable code.
- the instructions can include processor-specific instructions generated by a compiler and/or an interpreter from code written in any suitable computer-programming language, such as C, C++, C#, Visual Basic, Java, Python, Perl, JavaScript, and ActionScript.
- the instructions can include a database server application, such as engine 112 .
- the executable code when executed by the processor 104 , can cause the computing device 102 to perform processes according to embodiments as explained in more detail in this disclosure.
- Memory 106 can also include a datastore 114 .
- the datastore 114 may be a relational database, a flat-file database, triplestore, cloud storage, or other data storage device or instrumentality. In other embodiments, the datastore 114 is separate from the computing device 102 but in communication with the computing device 102 through an input/output (I/O) interface 110 .
- I/O input/output
- the computing device 102 can share data with additional components through the I/O interface 110 .
- the I/O interface 110 can include a USB port, eSATA port, Ethernet port, serial bus interface, parallel bus interface, wireless connection interface, or any suitable interface capable of allowing data transfers between the computing device and another component.
- the additional components can communicate with I/O interface 110 over a network 116 .
- the network 116 can include the internet, an intranet, wide area network (WAN), local area network (LAN), virtual private network (VPN), or any suitable communications network that allows computing device 102 to communicate with other components.
- Network 116 may include one or more networks, and can operate using any desired protocols.
- the additional components can include a database 118 and a user device 120 .
- the database 118 may be remote from the computing device 102 .
- the database 118 can store various types of information as explained more fully below.
- the user device 120 can include a second computing device, such as a laptop or personal computer that is configured for processing commands to output a user interface to a display.
- the user device 120 is a display device that is coupled to the computing device 102 directly instead of over the network 116 .
- any desired portion of the system can be connected to a CT scanner or other device capable of performing computed tomography, in order to interoperate with or form a portion of the scanner or device, or function as otherwise desired to carry out processes as described in this disclosure.
- This exemplary system configuration is provided to illustrate configurations of certain embodiments. Other configurations and embodiments may of course be utilized.
- the basis material decomposition theory which may be leveraged in certain embodiments of the invention, leads to the expression of the energy-dependent (E), mass attenuation coefficient for Material A, ⁇ A (E), as a linear combination of the mass attenuation coefficients for two basis materials, Material 1 and Material 2.
- E energy-dependent
- ⁇ A mass attenuation coefficient for Material A
- ⁇ A mass attenuation coefficient for two basis materials
- ⁇ A is the density of material A and ⁇ 1A and ⁇ 2A are abstract densities physically characteristic of Material A.
- ⁇ A (E) is physically constrained to be non-negative for all energies E in the energy range where the basis material decomposition theory applies.
- ⁇ 1A or ⁇ 2A may be negative, but not both.
- ⁇ 1A or ⁇ 2A have units of density, they do not necessarily represent physically realizable mixtures of basis materials 1 and 2.
- Each of the densities, ⁇ A , ⁇ 1A and ⁇ 2A are effectively energy-independent at energies commonly used for whole-body CT imaging and so do not depend upon factors such as kVp.
- the basis material composition for material A with density ⁇ A is sometimes referred to using the tuple ( ⁇ 1A , ⁇ 2A ).
- CT A The CT value of Material A, CT A , as measured on a particular CT scanner, can be written as a linear combination of the basis material composition densities, ( ⁇ 1A , ⁇ 2A ) if the response of the CT scanner to a change in density for each of the basis materials is accurately modeled as being proportional to the density of each basis material.
- CT A can be predicted using an equation of the form:
- CT A ⁇ 1 ⁇ 1A + ⁇ 2 ⁇ 2A +C (Eq 2)
- Equation 2 shows that calibration of a CT scanner involves the determination of three parameters for a given CT scan protocol.
- the CT protocol parameters of particular interest in certain embodiments of the invention are parameters that affect the effective energy of the x-ray source.
- kVp is the primary protocol parameter under user-control that influences the effective energy of the x-ray source.
- Effective energy can also be influenced by changes in filtration of the x-ray source as well as by filtration of the x-ray source by attenuation due to the patient, phantoms and other materials within the field-of-view of the CT scanner.
- CT scanners used for medical applications are generally calibrated in Hounsfield Units (HUs) where air is defined to have a CT value of ⁇ 1000 HU and water is defined to have a CT value of 0 HU.
- HU Hounsfield Unit
- the constant in Equation 2 will typically have a value on the order of 1000.
- basis Material 1 is assumed to be water with a physical density of 1000 mg/cm3, then it can be seen that ⁇ 1 for a properly calibrated CT scanner should have a value of 1 HU/(mg/cm3).
- the set of calibration parameters ( ⁇ 1 , ⁇ 2 ,C) can be estimated by imaging a phantom or other object containing a set of standards with known basis material compositions.
- CT values for a set of three standards must be obtained for a particular CT scan protocol to calibrate a particular scanner; although, the number of standards may be reduced if one or more of the three calibration parameters are known a priori by convention, as is the case when using Hounsfield Units, or other means.
- Aqueous K 2 HPO 4 was selected by Cann and Genant as the bone mineral standard for their pioneering work in the development of QCT for vertebral bone mineral density assessment.
- C. Cann, H. Genant “ Precise measurement of vertebral mineral content using computed tomography ,” J Comput Assist Tomogr 4, 493-500 (1980), incorporated herein by this reference).
- Their CT calibration scheme involves measuring CT values for a set of aqueous solutions containing varying and known amounts of K 2 HPO 4 .
- a linear fit to CT value versus K 2 HPO 4 concentration then provides a mechanism for expressing average CT value in terms of an equivalent aqueous K 2 HPO 4 concentration that would yield the same average CT value.
- the calibration slope derived with this approach neither defines the response of the CT scanner to a change in water density nor a change in K 2 HPO 4 density, but rather describes calibration relative to a controlled variation in simultaneous changes in densities of both of these materials.
- ⁇ H2O (0) is the density of pure water, 998.2 mg/cm3, at 20° C.
- the basis material composition relative to the water and K 2 HPO 4 basis set for an aqueous K 2 HPO 4 solution with a K 2 HPO 4 density of ⁇ K2HPO4 is ( ⁇ H2O ( ⁇ K2HPO4 ), ⁇ K2HPO4 ).
- This information along with CT value measurements for the calibration standards can be used in Equation 2 to determine ( ⁇ 1 , ⁇ 2 ,C). It can be shown in the case of a CT scanner calibrated in HU's that ⁇ 1 has a value of 1.002 HU/(mg/cm3) for pure changes in water density at 20° C.
- ⁇ 2 has a value that is equal to the slope derived using the conventional Cann and Genant calibration approach plus 0.1842 where ⁇ 2 characterizes of the CT scanner to pure changes in K 2 HPO 4 density.
- the value of C in this case is ⁇ 1000 HU.
- subcutaneous fat from a patient could be used as an internal patient control (assuming, for example that the DEQCT signal for subcutaneous fat variation is sufficiently small across subjects).
- Imaging methods such as DXA and DEQCT can benefit from use of documented expressions derived from the theory of basis material decomposition to estimate equivalent densities of basis materials for unknown materials from images acquired at two different x-ray energies.
- a method for accomplishing this follows from Equation 2 when it is recognized that the basis material decomposition ( ⁇ 1A , ⁇ 2A ) is independent of kVp while the calibration of a CT scanner, ( ⁇ 1 , ⁇ 2 ), depends upon kVp.
- CT 80 and CT 140 Two such measurements as denoted as CT 80 and CT 140 , respectively. It is sought to derive estimates of the densities for two basis materials that accurately represent the mass attenuation coefficient for the unknown material. For simplicity, use is made of H 2 O and K 2 HPO 4 as the basis set, so that it is sought to derive an estimate for ( ⁇ H2O , ⁇ K2HPO4 ) for the unknown material. To accomplish this requires knowledge of the CT-scanner-dependent response for each of the basis materials at both 80 kVp and 140 kVp.
- CT 80 - C 80 CT 140 - C 140 ( ⁇ H ⁇ ⁇ 2 ⁇ ⁇ O , 80 ⁇ K ⁇ ⁇ 2 ⁇ ⁇ HPO ⁇ ⁇ 4 , 80 ⁇ H ⁇ ⁇ 2 ⁇ ⁇ O , 140 ⁇ K ⁇ ⁇ 2 ⁇ ⁇ HPO ⁇ ⁇ 4 , 140 ) ⁇ ( ⁇ H ⁇ ⁇ 2 ⁇ ⁇ O ⁇ K ⁇ ⁇ 2 ⁇ ⁇ HPO ⁇ ⁇ 4 ) ( Eq ⁇ ⁇ 4 )
- Equation 1 provides an alternative method for deriving basis material decompositions if the energy-dependent, mass attenuation coefficients for a pair of suitable basis materials and materials that are to be characterized using the basis material decomposition method are known or can be approximated to sufficient accuracy.
- material calibration can include calibrating for one or more reference materials and two or more basis materials.
- One or more reference material CT pixel values are acquired using the CT scanner at a first, and if desired a second (and perhaps other) energy levels, which may be received by a system according to certain embodiments of the invention.
- One or more reference material CT pixel values may be extracted for one or more reference material regions of interest. Equations that incorporate, for one or more energy levels, this value, and expected reference material pixel value in terms of a linear combination of basis material densities and an offset, may be solved for the linear weights of the basis materials and an offset. These calibrated linear weights and offset can be used with CT pixel values acquired using the scanner operating at the relevant energy level to image an unknown material, in order to characterize basis material compositions for the unknown material and thereby characterize the unknown material.
- Multiple Energy CT image data is acquired at the relevant energy levels and segmented into one or more regions of interest (ROI's) within one or more anatomically, physiologically or functional areas of the patient's body. This step can occur at times different from the times at which basis material calibration occurs.
- ROI regions of interest
- average pixel value may be measured, received by systems according to certain embodiments of the invention, and used in a set of equations that incorporate expected ROI pixel values in terms of unknown basis material composition for the ROI and the calibrated linear weights and offset. These equations may be solved for the ROI basis material composition, which can be used to characterize the area of the patient's body that corresponds to the ROI.
- the characterization is a two dimensional result (in terms of weights for the two basis materials) and can be presented as a point in two dimensional space. It may also be presented in two dimensional space relative to other materials of interest characterized in the space. It can also be presented as a point along a line or a curve defined relative to other materials of interest characterized in the space or a line or a curve otherwise defined in the space. It can also be presented as a projection onto a line or curve defined relative to other materials of interest characterized in the space or a line or curve otherwise defined in the space.
- Projection methods include projections onto a line or a curve from a particular direction or projection by finding the shortest distance from the point to the desired reference line or curve or other algorithm that uniquely maps a point in the 2D basis material space to a unique point on a reference line or curve.
- Such processes and systems are useful for characterizing unknown material such as (among others) pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue. They are also useful for distinguishing tissue, or portions of tissue from other portions of tissue, including, for example, fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
- Single Energy CT image data is acquired at the relevant energy level and segmented into one or more regions of interest (ROI's) within one or more anatomically, physiologically or functional areas of the patient's body. This step can occur at times different from the times at which basis material calibration occurs.
- ROI regions of interest
- average pixel value may be measured, received by systems according to certain embodiments of the invention, and used in an equation that incorporates an expected ROI pixel value in terms of unknown basis material composition for the ROI and the calibrated linear weights and offset. This equation defines a line in two dimensional space (in terms of weights for the two basis materials).
- Another line or curve in this two dimensional space is defined by a physiological condition such as an observed or modeled bone mineral density or liver fat content.
- the intersection of the line and the reference line or curve may be reported to characterize the area of the patient's body corresponding to the ROI. If the line and a reference curve do not intersect, then the minimum distance between the line and curve may be reported to characterize the area of the patient's body corresponding to the ROI.
- the concept of characterizing a patient-specific measurement line to a reference line or curve can be generalized to comparison of the patient-specific measurement to two or more such reference lines resulting in characterizing the patient measurement by two or more intersection points or other unique measures relative to the multiple reference lines or curves.
- Such processes and systems are useful for characterizing unknown material such as (among others) pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue. They are also useful for distinguishing tissue, or portions of tissue from other portions of tissue, including, for example, fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
- the computer devices, processor devices, tangible medium embodiment code, and tangible computer readable media that includes program code can perform the following actions, as shown in FIG. 2 :
- Characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets; 4.
- the presentation can be in the form of graphical output such as on a monitor, or data that can be interpreted and used by other computing devices, processors, CT scanners, or other equipment.
- the computer devices, processor devices, tangible medium embodiment code, and tangible computer readable media that includes program code can also perform any or all of the following actions:
- CT pixel values to present an unknown material line in two dimensional basis material space, in order to characterize the unknown material.
- Receiving imaging information carrying out methods described above, and presenting information for the unknown material at a first energy level and a second energy level.
- Receiving imaging information for the unknown material can include receiving imaging information for a tissue, mixture, or substance.
- Such processes and systems are useful for characterizing unknown material such as (among others) pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue. They are also useful for distinguishing tissue, or portions of tissue from other portions of tissue, including, for example, fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
- steps 3 in the embodiment stated above can occur at a time (and/or using a CT scanner) different than step 4 in that embodiment; thus, some or all of step 3, that is, characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets can occur at a time (and/or using a CT scanner) different from some or all of step 4, that is from imaging information related to an unknown material received from the CT scanner at the first energy level; extracting at least one unknown material CT pixel value representative of a pre
- a reference material is provided in the form of CT calibration phantom marketed by Mindways Software, Austin, Tex., with known density and mass attenuation coefficient.
- Water and K 2 HPO 4 are used as basis materials, and the XCOM database of National Institute of Standards and Technology is used for obtaining mass attenuation coefficients for the two basis materials (http//physics.nist.gov/PhysRefData/Xcom/html/xcom1.html).
- a linear response model is used as follows:
- CT Pixel Value (Water Equiv Density)* A +(K 2 HPO 4 Equiv Density)* B+C (Eq 8)
- Water Equivalent Density and “K2HPO4 Equivalent Density” are the two numbers that describe the basis material composition for some material.
- a linear response model is not necessary, but appears to be sufficient for many practical applications of interest relative to CT imaging of human subjects. It is not necessary, but it appears to be sufficient to assume the coefficient “A” has a value of (1 HU/(mg/cm3)). This is essentially what would be expected for a CT scanner properly calibrated relative to the Hounsfield Unit scale. It is not necessary, but in many circumstances it is reasonable to assume the offset coefficient “C” has a value of ( ⁇ 1000 HU).
- certain embodiments of the invention can be reduced to imaging a single reference standard to gather sufficient information to solve a single equation in a single unknown (“B”).
- Multiple reference materials can be used leading to a set of equations in one or more unknowns that can be solved for the unknowns using various methods as long as the problem is not underdetermined (e.g., trying to use one reference material and one equation to extract 3 unknowns (e.g., “A”, “B” and “C”) is not going to lead to a useful representation of the CT scanner response in the context of this invention).
- CT pixel value measurement may be represented in a pixel value space defined by a tool used for deriving basis material composition estimates from what amounts to physical principles rather than from, say, x-ray attenuation measurements.
- a linear response model is suitable and that A, B and C in Eq 8 are known.
- CT Pixel Value now refers to a measurement from a CT image
- Water Equivalent Density and “K 2 HPO 4 Equivalent Density” are now unknowns. Solving Eq. 8 for “K 2 HPO 4 Equivalent Density”:
- FIG. 3 shows that the same information can be expressed or presented as a line ( 310 ). Additionally, the “Fat (80 kVp)” line can be presented if desired with one or more lines relating to reference materials.
- FIG. 4 shows “Fat (80 kVp)” line 410 with line 412 representing the reference material (here, a bone mineral density calibration standard or phantom marketed by Mindways Software, Austin, Tex.).
- This bone density calibration standard is based upon attenuation properties of aqueous solutions of K 2 HPO 4 .
- FIG. 4 shows the line 410 in the basis material space relative to the bone mineral density standard line 412 —both are calculated at 80 kVp based upon the measured response of the CT scanner at that energy level.
- the intersection point for these two lines is the equivalent bone mineral density of the fat measurement.
- the bone mineral density at the intersection with the K 2 HPO 4 line is ⁇ 69.4 mg/cm3 which is the same as the K 2 HPO 4 density at the intersection point.
- a second CT scan is obtained of the subcutaneous fat at a second energy level of 120 kVp.
- the CT scanner response is different at 120 kVp, so the fat measurement line at 120 kVp is different.
- the bone mineral density standard is, however, a fixed standard in this space.
- FIG. 5 which can be generated and presented using systems such as shown in FIG. 1 and according to a process such as shown in FIG. 2 , shows lines for the fat measurements at both 80 kVp and 120 kVp, together with the bone mineral density standard.
- the water equivalent density scale has been expanded from that scale as shown in FIG. 4 .
- the fat measurement line at 80 kVp is shown as line 510 ; the bone mineral density standard line as 512 , and the 120 kVp fat measurement line as 514 . Note that the intersection of the Fat (120 kVp) measurement line 514 and the bone density standard line 512 is not the same as the point at which the Fat (80 kVp) measurement line 510 intersects the bone mineral density standard line 512 .
- the equivalent bone mineral density of the 120 kVp subcutaneous fat measurement is about ⁇ 87.5 mg/cm3.
- the bone mineral density standard is not particularly convenient for characterizing fat because the bone mineral density estimate will change depending upon kVp and beam hardening (which affects the average energy of the x-rays hitting the CT scanner detector—beam hardening can be intentional, such as pre-filtration of the x-ray source or coincidental, such as degree of hardening due to patient size).
- the bone mineral density standard line did not change as a result of changing how the CT scanner is used.
- the density for example, of an aqueous solution of 100 mg/cm3 of K 2 HPO 4 does not change depending upon whether measurement occurs at 80 kVp or 120 kVp. It will always be that composition that is at the end of the black BMD calibration standard line which intersects the horizontal line corresponding to a K 2 HPO 4 equivalent density of 100 mg/cm3.
- the unknown subcutaneous fat has the same physical attribute. Its density is not changed according to whether it is measured at 120 kVp or 80 kVp.
- the basis material composition for the subcutaneous fat that is consistent across the two measurements is at the intersection point of the 80 kVp and 120 kVp fat lines.
- estimates of the bone mineral density calibration standard for two hypothetical scenarios can also be considered: (1) bone with a marrow space filled with pure yellow marrow and (2) bone with the marrow space filled with pure red marrow.
- the fat 80 kVp line bears numeral 610 ; the bone mineral density line 612 ; the fat 120 kVp line 614 ; the yellow marrow line 616 and the red marrow line 618 .
- These calculations are based upon atomic composition and density estimates from the literature, and are scanner independent and do not rely on a material that mimics bone. That is, there is no attempt to construct and measure with x-rays some yellow marrow mimic or some red marrow mimic. Instead, the lines for these two scenarios are constructed and presented based on mass attenuation coefficients obtained from the literature.
- the line generated and/or presented in the basis material space which can be generated and presented using systems such as shown in FIG. 1 and according to a process such as shown in FIG. 2 , can be compared to a standard reference line such as a bone mineral density line.
- the comparison can take the form of an intersection of the two lines, or shortest distance between the two lines.
- the point could be estimated where “normal” liver sits in the basis material space. That point can be determined by experiment, such as the subcutaneous fat basis material composition point addressed above, or by calculation using known mass attenuation coefficients. Then, the closest distance between the normal liver reference point and the experimental measurement line can be considered.
- a line can be generated from the normal liver point based upon replacing liver with fat to calculate amount of fat based on where the experimental measurement line crosses the liver/fat line.
- the CT measurement can be dual energy measurements, in which case a basis material composition can be estimated for the unknown tissue. That estimate can allow measurements relative to the basis material composition rather than relative to the shortest distance between lines/points/curves.
- two (or more) components or constituent materials of a region of interest or measurement volume may be characterized, such as according to the process shown in FIG.
- fat in pancreatic tissue in (for example and without limitation) fat in pancreatic tissue, yellow/red bone marrow tissue, white-matter/grey-matter brain tissue, fat in muscle tissue, healthy/necrotic heart tissue, any healthy/necrotic soft tissue, any healthy or normal/diseased or pathological tissue, normal/diseased lung tissue, calcified stones/kidney tissue, and iron in liver tissue.
- DEQCT quantification methods is characterization of “bright” deposits.
- Kidney stones may be predominantly composed of uric acid or calcium.
- the origins of these different classes of kidney stones are different and may have implications regarding treatment.
- Deposition of uric acid crystals in joint spaces (gout) is also characteristically different than other forms of arthritis that lead to growth of calcium-based osteophytes.
- Uric acid and calcium have distinctly different basis material compositions.
- kidney stones and certain other structures
- the kidney stones themselves will have appreciable amounts of calcium in the case of hydroxyapatite and oxalate stones, while cystine has sulfur and uric acid is nitrogen rich.
- Each of these materials has fairly different characteristics when viewed in the basis material space due to the significant difference in atomic numbers of these species.
- Uric acid to some extent looks like water (in the sense of being parallel to the water axis), but is actually slightly below the water axis reflecting the effective atomic number for uric acid being slightly less than water.
- Cystine has a pronounced positive “K2HPO4” character due to sulfur (atomic number 16—double that of oxygen).
- the calcium containing materials (Ca atomic number 20) are still higher, in the range of 2000 mg/cm3 to 3000 mg/cm3 range for K2HPO4 component.
- FIG. 7 which can be generated and presented using systems such as shown in FIG. 1 and according to a process such as shown in FIG. 2 , and which shows lines that relate to various potential reference materials relevant to kidney stones, the intersection point of all of the reference lines is the “water” point.
- Some moving to the right (and downward) along the “uric acid” reference line 710 would represent a mixture that looks mostly like water with some of that water-like material replaced with uric acid crystals.
- a measurement classifier that could include factors such as estimated stone size, signal strength (basically how far up the K2HPO4 axis a particular measurement falls for the size of the stone relative to the sampled volume), as well as proximity to one of the reference lines above to identify the nature of the kidney stone (calcium/cystine/uric acid).
- some embodiments relate to the projection of a single value received from a CT scan back into a two dimensional space, which single value defines a line in the two dimensional (basis material) space.
- the two dimensional space is very rich because, among other things, various theoretical and experimental results in the two dimensional space can be generated and explored in a controlled (calibrated) manner without significant consideration for the source of the experimental data (i.e., what CT scanner was used and in exactly what mode it was used).
- the one-dimensional to two dimensional mapping can be reversed: Once something interesting is apparent in the device-independent two dimensional space, it can be transformed back into one dimension by applying the reverse transform derived for a specific CT scanner used in a particular manner.
- Such a transformation can be valuable, for instance, where there are various calibration problems that occur because by pixel values observed on a particular CT scanner depend upon the CT scanner and how it is used.
- One such example is above, where the subcutaneous fat example used above resulted in estimated CT values of ⁇ 108.4 HU at 80 kVp and ⁇ 83.7 HU at 120 kVp.
- MRS Magnetic Resonance Spectroscopy
- MRS results are presented as a fat ratio.
- the ratio of (fat signal) to (fat signal+water signal) is reported with MRS when quantifying marrow fat.
- This ratio is expected to approach unity if the marrow space is dominated by the presence of yellow marrow (fatty marrow) and it is expect to approach zero if there is little fat in the marrow space (as might be the case for mostly red marrow). So this is neither an absolute quantification of the amount of water present nor an absolute quantification of the amount of fat present, and the interpretation of the measured ratio is confounded by the presence of other materials in a tissue environment as highly diverse as cancellous bone. That is, there may be molecular species present in the bone (and bone marrow) environment that are neither accurately described as “fat” or “water” that might contaminate either or both “signals” and their ratio.
- One approach to this marrow fat problem in the DEQCT space is to consider projecting the DEQCT measurements onto a line that passes through the basis material composition expected for water (e.g., as derived using the XCOM data) as well as the basis material composition expected for pure yellow marrow (e.g., as derived by using the XCOM data). That projection of DEQCT results onto such a line can provide a simple means for mapping DEQCT measurements to the MRS measurement space is not obvious. It is possible, for example, to consider other lines such as projection onto a line perpendicular to the bone density reference standard which nominally characterizes the expected change in signal due to the displacement of water from a hydrated cadaveric bone sample by the presence of mineralized bone.
- Another option may be to project the DEQCT measurements onto lines of constant BMD for a single energy QCT measurement. These lines in the basis material space represent mixtures of atomic constituents that can lead to the same BMD estimate for single energy QCT at a particular energy even though they had different relative amounts of bone and marrow and marrow composition. It is possible, for example, to contemplate the replacement of some red marrow with yellow marrow (which by itself would decrease the estimate bone mineral density) while adding some additional mineralized bone to the characterized volume such that the overall BMD remained constant.
- FIGS. 8 and 9 show pair MR and DEQCT fat measurements.
- the square of the correlation coefficient (R) in this case is about 0.7 which implies R is about 0.84, which is a relatively strong correlation.
- FIG. 9 shows the change in the fat signal as measured by both methods.
- both the DEQCT and MRS measurements are potentially contaminated by materials unrelated to those that are really the subject of quantification. Subtracting the baseline measurement from a follow-up measurement might subtract out a substantial portion of the signal contamination if the degree of contamination does not change appreciably between baseline and follow-up. When one considers the difference in signals, it is apparent that R in this case for this data is about 0.93 (R2 of about 0.86). As such, it may be that DEQCT predicts the change in MRS fat signal as accurately as or more accurately than it does in predicting the absolute MRS fat signal.
- a computing device can include any suitable arrangement of components that provide a result conditioned on one or more inputs.
- Suitable computing devices include multipurpose microprocessor-based computer systems accessing stored software that programs or configures the computing system from a general-purpose computing apparatus to a specialized computing apparatus implementing one or more embodiments of the present subject matter. Any suitable programming, scripting, or other type of language or combinations of languages may be used to implement the teachings contained herein in software to be used in programming or configuring a computing device.
- Embodiments of the methods disclosed herein may be performed in the operation of such computing devices.
- the order of the blocks presented in the examples above can be varied—for example, blocks can be re-ordered, combined, and/or broken into sub-blocks. Certain blocks or processes can be performed in parallel.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Optics & Photonics (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Toxicology (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Crystallography & Structural Chemistry (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
According to certain embodiments of the invention, computer hardware and software can perform single or dual-energy decomposition calculations and predict CT values for either single-energy or dual energy CT from basis material decomposition estimates derived using theoretical models. Using that information, the computer hardware and software can (among other things) present one dimensional information from a CT scanner, such as representative pixel value, in two dimensions, such as a line in two dimensional basis material space.
Description
- NOT APPLICABLE
- NOT APPLICABLE
- This disclosure relates generally to computer hardware and software for calibrating single or multiple energy computed tomography (CT) devices. The hardware and software runs, displays, provides and otherwise uses electronic content; more particularly, such hardware and software operate on certain content that is independent of the CT scanner being used, in order (among other things) to calibrate the scanner and to present one dimensional information the CT scanner generates (such as representative pixel values) in two dimensional basis material space.
- Quantitative Computed Tomography (QCT) is an established method for assessing bone, for example, bone mineral density. In principle other tissue types or traits are amenable to characterization with QCT; although, presently non-bone applications of QCT have not evolved to common clinical use.
- Historically QCT bone densitometry has relied upon calibration of a CT image relative to a calibration standard with bone mineral density expressed in density units relative to the selected standard. While internal standards are sometimes used, such as calibration relative to fat and muscle attenuation as measured within a subject, the use of a standard external to a subject with well-known and stable x-ray attenuation properties generally provides advantages in terms of long-term measurement accuracy and precision. Most published descriptions of external calibration standards for bone densitometry involve the use of bone-mimicking standards. Regardless of the details, most calibration methods used for QCT rely upon standardizing QCT measurements to a one-dimensional scale.
- Limitations of using a one-dimensional scale for reporting bone mineral density are known. Those include effect on such use of variations in marrow composition, as to which the impact of marrow composition on bone mineral density estimates have been explored. Dual energy QCT (DEQCT) has been used to characterize marrow-fat content and to quantify bias in bone mineral density estimates induced by marrow fat.
- In this vein, the physical density of fat/yellow marrow (YM) and red marrow (RM) are different. The apparent bone marrow density of red marrow is approximately 50 mg/cm3 higher than the bone marrow density of yellow marrow. Since cancellous bone is composed mostly of marrow (>80%), change in marrow composition may confound true change in trabecular bone density if they are not distinguished. Additionally, chemotherapy and radiation cause marrow damage in vitro and in vivo. Although the concept of DEQCT to estimate bone mineral density corrected for marrow change was proposed earlier, DEQCT has not generally been translated to the study of osteoporosis in clinical practice.
- The basis material decomposition theory published by Alvarez and Macovski provides an elegant framework for describing and understanding DEQCT. (R. Alvarez, A. Macovski, “Energy-selective reconstructions in x-ray computerised tomography,” Physics in Medicine and Biology 21, 733-744 (1976); L. Lehmann, R. Alvarez, A. Macovski, W. Brody, N. Pelc, S. Riederer, A. Hall, “Generalized image combinations in dual KVP digital radiography,” Medical Physics 8, 659 (1981), both of which are incorporated herein by this reference.) Applicability of this theory is predicated upon a set of assumptions regarding x-ray attenuation properties of material being imaged at two different x-ray energies. Of particular importance is that the materials being imaged must not have K-edges in or slightly below the x-ray energy windows used for imaging. This assumption is generally reasonably satisfied in the diagnostic x-ray energy regime for a broad range of substances including biological tissues. Dual energy x-ray absorptiometry, DXA, is perhaps the most successful conventional demonstration of the applicability of the basis material decomposition to the interpretation of x-ray images acquired at different energies.
- The diagnostic x-ray energy regime is also interesting in that theoretical models and tabulated x-ray data can be used to predict mass attenuation coefficients for tissues of biological interest with accuracies on the order of 1% to 2%. Energy-dependent mass attenuation coefficients can be used to derive basis material decompositions that are analogous to those measured with dual-energy x-ray methods.
- According to certain embodiments of the invention, computer hardware and software can perform single or dual-energy decomposition calculations and predict CT values for either single-energy or dual energy CT from basis material decomposition estimates derived using theoretical models. Using that information, the computer hardware and software can present one dimensional information from a CT scanner, such as representative pixel value for an unknown material, in two dimensions, such as a line in two dimensional basis material space. In some embodiments, for example, the hardware and software can predict CT values for heterogeneous mixtures of multiple tissue types, such as a model of bone including collagen, bone mineral described as calcium hydroxyapatite, and a mixture of red and yellow marrow, based on atomic compositions for the tissue types. In some embodiments, the hardware and software can interpret and present observed CT values from a single energy CT scan in a basis material space relative to the model tissue constituents; for example, a representative pixel value acquired from a CT scan can be presented in two dimensional basis material space with attendant new benefits for CT calibration. As another example, pixel values from a dual energy CT scan can be presented as lines in basis material space, whose intersection gives information about the basis material composition of the scanned material.
- According to some embodiments, a method is executed by a computing device that includes these steps:
- 1. For a reference material and for each of a first basis material and a second basis material, receiving an energy-dependent, mass attenuation coefficient for the material as a function of its elemental composition, molecular composition or elemental and molecular composition;
2. Deriving a basis material composition of the reference material relative to the first basis material and the second basis material;
3. Characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets;
4. From imaging information related to an unknown material received from the CT scanner at the first energy level; extracting at least one unknown material CT pixel value representative of a predetermined region of the unknown material; and presenting the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the reference material relative to the first basis material and the second basis material. - Some embodiments relate to a system that includes a processor device and a tangible medium embodiment code. The code can be executable by the processor to cause the system to perform certain actions. The actions can include:
- 1. For a reference material and for each of a first basis material and a second basis material, receiving an energy-dependent, mass attenuation coefficient for the material as a function of its elemental composition, molecular composition or elemental and molecular composition;
2. Deriving a basis material composition of the reference material relative to the first basis material and the second basis material;
3. Characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets;
4. From imaging information related to an unknown material received from the CT scanner at the first energy level; extracting at least one unknown material CT pixel value representative of a predetermined region of the unknown material; and presenting the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the reference material relative to the first basis material and the second basis material. - Some embodiments relate to a tangible computer readable medium that includes program code that is executable by a processor to perform actions. The program code includes code for causing the processor to carry out the following method:
- 1. For a reference material and for each of a first basis material and a second basis material, receiving an energy-dependent, mass attenuation coefficient for the material as a function of its elemental composition, molecular composition or elemental and molecular composition;
2. Deriving a basis material composition of the reference material relative to the first basis material and the second basis material;
3. Characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets;
4. From imaging information related to an unknown material received from the CT scanner at the first energy level; extracting at least one unknown material CT pixel value representative of a predetermined region of the unknown material; and presenting the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the reference material relative to the first basis material and the second basis material. - In some embodiments, the computer devices, processor devices, tangible medium embodiment code, and tangible computer readable medium that includes program code can perform any or all of the following actions:
- Using one or more CT pixel values to present an unknown material line in two dimensional basis material space, in order to characterize the unknown material.
- Mapping the characterization of the response of the CT scanner to the first basis material and the second basis material to a single dimensional value relating to the reference material.
- Mapping the representation of the unknown material CT pixel values as a function of the characterization of the response of the CT scanner to the first basis material and the second basis material to a single dimensional value relating to the reference material.
- Comparing an unknown material line in two dimensional basis material space to a line for the reference material in two dimensional basis material space in order further to characterize the unknown material.
- Receiving imaging information, carrying out methods described above, and presenting information for the unknown material at a first energy level and a second energy level.
- Mapping the at least one unknown material CT pixel value back to a family of possible basis material compositions that are predicted to exhibit the same CT pixel value as the at least one unknown material CT pixel value.
- Receiving imaging information for the unknown material can include receiving imaging information for a tissue, mixture, or substance.
- Characterizing unknown material selected from the group consisting of: pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue.
- Distinguishing tissue, or portions of tissue from other portions of tissue, the tissue selected from the group consisting of: fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
- In some embodiments, processes in step 3 in the embodiments stated above can occur at a time (and/or on a CT scanner) different than step 4 in those embodiments; thus, some or all of step 3, that is, characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets can occur at a time (or on a CT scanner) different from some or all of step 4, that is from imaging information related to an unknown material received from the CT scanner at the first energy level; extracting at least one unknown material CT pixel value representative of a predetermined region of the unknown material; and presenting the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the reference material relative to the first basis material and the second basis material.
-
FIG. 1 is a functional block diagram of a system that is capable of carrying out methods according to certain embodiments of the invention. -
FIG. 2 is a flow diagram showing process flow for one method according to an embodiment of the invention. -
FIG. 3 is a diagram that shows a first presentation according to an embodiment of the invention. -
FIG. 4 is a diagram that shows a second presentation according to an embodiment of the invention. -
FIG. 5 is a diagram that shows a third presentation according to an embodiment of the invention. -
FIG. 6 is a diagram that shows a fourth presentation according to an embodiment of the invention. -
FIG. 7 is a diagram that shows a fifth presentation according to an embodiment of the invention. -
FIG. 8 is a graph showing pairs of MRS and DEQCT fat measurements. -
FIG. 9 is a graph showing changes in a fat signal as measured by both MRS and DEQCT methods. - a. Illustrative System Implementation
-
FIG. 1 depicts a system that is configured, according to certain embodiments of the invention, to do any or all of at least the following: Perform single or dual-energy decomposition calculations and predict CT values for either single-energy or dual energy CT from basis material decomposition estimates derived using theoretical models; present one dimensional information from a CT scanner, such as representative pixel value for an unknown material, in two dimensions, such as a line in two dimensional basis material space; predict CT values for heterogeneous mixtures of multiple tissue types, such as a model of bone including collagen, bone mineral described as calcium hydroxyapatite, and a mixture of red and yellow marrow, based on atomic compositions for the tissue types; present observed CT values from a single energy CT scan in a basis material space relative to the model tissue constituents, such as, for example, presenting a representative pixel value acquired from a CT scan in two dimensional basis material space; present, based on pixel values from a dual energy CT scan, lines in basis material space, whose intersection gives information about the basis material composition of the scanned material; and/or carry out any of the methods mentioned in this disclosure. Other system embodiments may be utilized. - The system generally includes a
computing device 102 having aprocessor 104 that can execute code stored on a tangible computer-readable medium, such as amemory 106, to cause thecomputing device 102 to perform processes as described more fully herein. Thecomputing device 102 may be any device that can process data and execute code that is a set of instructions to perform actions. Examples of thecomputing device 102 include a database server, a web server, desktop personal computer, a laptop personal computer, a server device, a handheld computing device, and a mobile device. - Examples of the
processor 104 include a microprocessor, an application-specific integrated circuit (ASIC), a state machine, or other suitable processor. Theprocessor 104 may include one processor or any number of processors. Theprocessor 104 can access code stored in thememory 106 via abus 108. Thememory 106 may be any non-transitory computer-readable medium configured for tangibly embodying code and can include electronic, magnetic, or optical devices. Examples of thememory 106 include random access memory (RAM), read-only memory (ROM), a floppy disk, compact disc, digital video device, magnetic disk, an ASIC, a configured processor, flash drive, thumb drive, cloud memory capacity, or other storage device or capacity. Thebus 108 may be any device capable of transferring data between components of thecomputing device 102. Thebus 108 can include one device or multiple devices. - Instructions can be stored in the
memory 106 as executable code. The instructions can include processor-specific instructions generated by a compiler and/or an interpreter from code written in any suitable computer-programming language, such as C, C++, C#, Visual Basic, Java, Python, Perl, JavaScript, and ActionScript. The instructions can include a database server application, such asengine 112. The executable code, when executed by theprocessor 104, can cause thecomputing device 102 to perform processes according to embodiments as explained in more detail in this disclosure. -
Memory 106 can also include adatastore 114. Thedatastore 114 may be a relational database, a flat-file database, triplestore, cloud storage, or other data storage device or instrumentality. In other embodiments, thedatastore 114 is separate from thecomputing device 102 but in communication with thecomputing device 102 through an input/output (I/O)interface 110. - The
computing device 102 can share data with additional components through the I/O interface 110. The I/O interface 110 can include a USB port, eSATA port, Ethernet port, serial bus interface, parallel bus interface, wireless connection interface, or any suitable interface capable of allowing data transfers between the computing device and another component. The additional components can communicate with I/O interface 110 over anetwork 116. Thenetwork 116 can include the internet, an intranet, wide area network (WAN), local area network (LAN), virtual private network (VPN), or any suitable communications network that allowscomputing device 102 to communicate with other components.Network 116 may include one or more networks, and can operate using any desired protocols. - The additional components can include a
database 118 and auser device 120. Thedatabase 118 may be remote from thecomputing device 102. In some embodiments, thedatabase 118 can store various types of information as explained more fully below. Theuser device 120 can include a second computing device, such as a laptop or personal computer that is configured for processing commands to output a user interface to a display. In some embodiments, theuser device 120 is a display device that is coupled to thecomputing device 102 directly instead of over thenetwork 116. In some embodiments, any desired portion of the system can be connected to a CT scanner or other device capable of performing computed tomography, in order to interoperate with or form a portion of the scanner or device, or function as otherwise desired to carry out processes as described in this disclosure. - This exemplary system configuration is provided to illustrate configurations of certain embodiments. Other configurations and embodiments may of course be utilized.
- b. Theory
- The basis material decomposition theory which may be leveraged in certain embodiments of the invention, leads to the expression of the energy-dependent (E), mass attenuation coefficient for Material A, σA(E), as a linear combination of the mass attenuation coefficients for two basis materials, Material 1 and Material 2. In particular
-
ρA·σA(E)=ρ1A·σ1(E)+ρ2A·σ2(E) (Eq 1) - where ρA is the density of material A and ρ1A and ρ2A are abstract densities physically characteristic of Material A. σA(E) is physically constrained to be non-negative for all energies E in the energy range where the basis material decomposition theory applies. Thus it is apparent that either ρ1A or ρ2A may be negative, but not both. So while ρ1A or ρ2A have units of density, they do not necessarily represent physically realizable mixtures of basis materials 1 and 2. Each of the densities, ρA, ρ1A and ρ2A are effectively energy-independent at energies commonly used for whole-body CT imaging and so do not depend upon factors such as kVp. For convenience, the basis material composition for material A with density ρA is sometimes referred to using the tuple (ρ1A, ρ2A).
- The CT value of Material A, CTA, as measured on a particular CT scanner, can be written as a linear combination of the basis material composition densities, (ρ1A, ρ2A) if the response of the CT scanner to a change in density for each of the basis materials is accurately modeled as being proportional to the density of each basis material. In this common case, CTA can be predicted using an equation of the form:
-
CTA=β1·ρ1A+β2·ρ2A +C (Eq 2) - where β1 and β2 are characteristic of the CT scanner and the CT scan protocol, and C is a constant generally related to convention. Equation 2 shows that calibration of a CT scanner involves the determination of three parameters for a given CT scan protocol. The CT protocol parameters of particular interest in certain embodiments of the invention are parameters that affect the effective energy of the x-ray source. kVp is the primary protocol parameter under user-control that influences the effective energy of the x-ray source. Effective energy can also be influenced by changes in filtration of the x-ray source as well as by filtration of the x-ray source by attenuation due to the patient, phantoms and other materials within the field-of-view of the CT scanner.
- CT scanners used for medical applications are generally calibrated in Hounsfield Units (HUs) where air is defined to have a CT value of −1000 HU and water is defined to have a CT value of 0 HU. When the HU scale is used, the constant in Equation 2 will typically have a value on the order of 1000. Additionally, if basis Material 1 is assumed to be water with a physical density of 1000 mg/cm3, then it can be seen that β1 for a properly calibrated CT scanner should have a value of 1 HU/(mg/cm3).
- Regardless of convention with respect to scaling CT scanner measurements, the set of calibration parameters (β1,β2,C) can be estimated by imaging a phantom or other object containing a set of standards with known basis material compositions. In general CT values for a set of three standards must be obtained for a particular CT scan protocol to calibrate a particular scanner; although, the number of standards may be reduced if one or more of the three calibration parameters are known a priori by convention, as is the case when using Hounsfield Units, or other means.
- Aqueous K2HPO4 was selected by Cann and Genant as the bone mineral standard for their pioneering work in the development of QCT for vertebral bone mineral density assessment. (C. Cann, H. Genant, “Precise measurement of vertebral mineral content using computed tomography,” J Comput Assist Tomogr 4, 493-500 (1980), incorporated herein by this reference). Their CT calibration scheme involves measuring CT values for a set of aqueous solutions containing varying and known amounts of K2HPO4. A linear fit to CT value versus K2HPO4 concentration then provides a mechanism for expressing average CT value in terms of an equivalent aqueous K2HPO4 concentration that would yield the same average CT value. The calibration slope derived with this approach neither defines the response of the CT scanner to a change in water density nor a change in K2HPO4 density, but rather describes calibration relative to a controlled variation in simultaneous changes in densities of both of these materials.
- Supplementing the information derived from the Cann and Genant calibration approach with information regarding the physical density of aqueous K2HPO4 solutions of known concentration leads to an approach that does yield (β1,β2,C) defined above. Using published solution density information, it can be shown that at 20° C. the density of water in a sample of an aqueous solution of containing K2HPO4 density is:
-
ρH2O(ρK2HPO4)=ρH2O(0)−0.1842·ρK2HPO4 (Eq 3) - where ρH2O(0) is the density of pure water, 998.2 mg/cm3, at 20° C. Thus the basis material composition relative to the water and K2HPO4 basis set for an aqueous K2HPO4 solution with a K2HPO4 density of ρK2HPO4 is (ρH2O(ρK2HPO4), ρK2HPO4). This information along with CT value measurements for the calibration standards can be used in Equation 2 to determine (β1,β2,C). It can be shown in the case of a CT scanner calibrated in HU's that β1 has a value of 1.002 HU/(mg/cm3) for pure changes in water density at 20° C. and β2 has a value that is equal to the slope derived using the conventional Cann and Genant calibration approach plus 0.1842 where β2 characterizes of the CT scanner to pure changes in K2HPO4 density. The value of C in this case is −1000 HU. Given uncertainties in the purity of water and the temperature of the water used when calibrating a CT scanner, and the precision of typical CT scan measurements, it is convenient to assume β1 has a value of 1 HU/(mg/cm3) in practical applications.
- It is also within the ambit of various embodiments of the invention to use human tissue as an internal control or basis material. For example, subcutaneous fat from a patient could be used as an internal patient control (assuming, for example that the DEQCT signal for subcutaneous fat variation is sufficiently small across subjects).
- (ii) Basis Material Decomposition Methods
- Imaging methods such as DXA and DEQCT can benefit from use of documented expressions derived from the theory of basis material decomposition to estimate equivalent densities of basis materials for unknown materials from images acquired at two different x-ray energies. A method for accomplishing this follows from Equation 2 when it is recognized that the basis material decomposition (ρ1A, ρ2A) is independent of kVp while the calibration of a CT scanner, (β1,β2), depends upon kVp.
- According to this method, consider the measurement of the CT value for an unknown material at two different energies. For simplicity, assume the energies correspond to imaging at 80 kVp and 140 kVp. Two such measurements as denoted as CT80 and CT140, respectively. It is sought to derive estimates of the densities for two basis materials that accurately represent the mass attenuation coefficient for the unknown material. For simplicity, use is made of H2O and K2HPO4 as the basis set, so that it is sought to derive an estimate for (ρH2O, ρK2HPO4) for the unknown material. To accomplish this requires knowledge of the CT-scanner-dependent response for each of the basis materials at both 80 kVp and 140 kVp. The symbol β with subscripts is used for material and energy to represent the scanner response terms. Finally, allowance is made for a scanner-dependent offset, C, that, in general, may depend upon kVp. Accordingly, the following set of equations is arrived at that must be solved to generate the basis material estimates:
-
- Direct matrix inversion leads to the following solution for the basis material densities:
-
- Thus measuring the CT value for a particular sample at two different kVp settings for a calibrated CT scanner yields a set of two linear equations in two unknowns, (ρ1A, ρ2A), that can be solved using standard methods.
- Equation 1 provides an alternative method for deriving basis material decompositions if the energy-dependent, mass attenuation coefficients for a pair of suitable basis materials and materials that are to be characterized using the basis material decomposition method are known or can be approximated to sufficient accuracy.
- c. Processes
- Certain processes according to certain embodiments of the invention can be carried out using certain embodiments of systems of the invention, with reference to some or all of the following:
- (i) Basis Material Calibration
- In accordance with theory addressed above, material calibration can include calibrating for one or more reference materials and two or more basis materials. First, expected energy dependent mass attenuation coefficients are calculated or provided. Second, basis material composition for the reference material is determined by finding a linear combination of basis material mass attenuation coefficients that is the best fit to the mass attenuation coefficient of the reference material.
- (ii) Two Dimensional Scanner Calibration
- One or more reference material CT pixel values are acquired using the CT scanner at a first, and if desired a second (and perhaps other) energy levels, which may be received by a system according to certain embodiments of the invention. One or more reference material CT pixel values may be extracted for one or more reference material regions of interest. Equations that incorporate, for one or more energy levels, this value, and expected reference material pixel value in terms of a linear combination of basis material densities and an offset, may be solved for the linear weights of the basis materials and an offset. These calibrated linear weights and offset can be used with CT pixel values acquired using the scanner operating at the relevant energy level to image an unknown material, in order to characterize basis material compositions for the unknown material and thereby characterize the unknown material.
- (iii) Multiple Energy CT Patient Measurement
- In the context of a patient, Multiple Energy CT image data is acquired at the relevant energy levels and segmented into one or more regions of interest (ROI's) within one or more anatomically, physiologically or functional areas of the patient's body. This step can occur at times different from the times at which basis material calibration occurs. For the relevant ROI, average pixel value may be measured, received by systems according to certain embodiments of the invention, and used in a set of equations that incorporate expected ROI pixel values in terms of unknown basis material composition for the ROI and the calibrated linear weights and offset. These equations may be solved for the ROI basis material composition, which can be used to characterize the area of the patient's body that corresponds to the ROI. The characterization is a two dimensional result (in terms of weights for the two basis materials) and can be presented as a point in two dimensional space. It may also be presented in two dimensional space relative to other materials of interest characterized in the space. It can also be presented as a point along a line or a curve defined relative to other materials of interest characterized in the space or a line or a curve otherwise defined in the space. It can also be presented as a projection onto a line or curve defined relative to other materials of interest characterized in the space or a line or curve otherwise defined in the space. Projection methods include projections onto a line or a curve from a particular direction or projection by finding the shortest distance from the point to the desired reference line or curve or other algorithm that uniquely maps a point in the 2D basis material space to a unique point on a reference line or curve.
- Such processes and systems are useful for characterizing unknown material such as (among others) pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue. They are also useful for distinguishing tissue, or portions of tissue from other portions of tissue, including, for example, fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
- (iv) Single Energy CT Patient Measurement
- In the context of a patient, Single Energy CT image data is acquired at the relevant energy level and segmented into one or more regions of interest (ROI's) within one or more anatomically, physiologically or functional areas of the patient's body. This step can occur at times different from the times at which basis material calibration occurs. For the relevant ROI, average pixel value may be measured, received by systems according to certain embodiments of the invention, and used in an equation that incorporates an expected ROI pixel value in terms of unknown basis material composition for the ROI and the calibrated linear weights and offset. This equation defines a line in two dimensional space (in terms of weights for the two basis materials). Another line or curve in this two dimensional space, which may be known, is defined by a physiological condition such as an observed or modeled bone mineral density or liver fat content. The intersection of the line and the reference line or curve may be reported to characterize the area of the patient's body corresponding to the ROI. If the line and a reference curve do not intersect, then the minimum distance between the line and curve may be reported to characterize the area of the patient's body corresponding to the ROI. The concept of characterizing a patient-specific measurement line to a reference line or curve can be generalized to comparison of the patient-specific measurement to two or more such reference lines resulting in characterizing the patient measurement by two or more intersection points or other unique measures relative to the multiple reference lines or curves. One possible application of such an approach would be cross-calibration of measurements to different standards. In clinical practice today, for example, bone mineral density is sometimes reported relative to the aqueous K2HPO4 density standard while in other instances it is reported relative to a calcium hydroxyapatite standard. Both standards may be represented simultaneously in the 2D basis material space and a patient measurement may be expressed relative to either or both standards.
- Such processes and systems are useful for characterizing unknown material such as (among others) pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue. They are also useful for distinguishing tissue, or portions of tissue from other portions of tissue, including, for example, fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
- d. Exemplary Processes
- In some embodiments, the computer devices, processor devices, tangible medium embodiment code, and tangible computer readable media that includes program code can perform the following actions, as shown in
FIG. 2 : - 1. For a reference material and for each of a first basis material and a second basis material, receiving an energy-dependent, mass attenuation coefficient for the material as a function of its elemental composition, molecular composition or elemental and molecular composition;
2. Deriving a basis material composition of the reference material relative to the first basis material and the second basis material;
3. Characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets;
4. From imaging information related to an unknown material received from the CT scanner at the first energy level; extracting at least one unknown material CT pixel value representative of a predetermined region of the unknown material; and presenting the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the reference material relative to the first basis material and the second basis material. - The presentation can be in the form of graphical output such as on a monitor, or data that can be interpreted and used by other computing devices, processors, CT scanners, or other equipment.
- In some embodiments, the computer devices, processor devices, tangible medium embodiment code, and tangible computer readable media that includes program code can also perform any or all of the following actions:
- Using one or more CT pixel values to present an unknown material line in two dimensional basis material space, in order to characterize the unknown material.
- Mapping the characterization of the response of the CT scanner to the first basis material and the second basis material to a single dimensional value relating to the reference material.
- Mapping the representation of the unknown material CT pixel values as a function of the characterization of the response of the CT scanner to the first basis material and the second basis material to a single dimensional value relating to the reference material.
- Comparing an unknown material line in two dimensional basis material space to a line for the reference material in two dimensional basis material space in order further to characterize the unknown material.
- Receiving imaging information, carrying out methods described above, and presenting information for the unknown material at a first energy level and a second energy level.
- Mapping the at least one unknown material CT pixel value back to a family of possible basis material compositions that are predicted to exhibit the same CT pixel value as the at least one unknown material CT pixel value.
- Receiving imaging information for the unknown material can include receiving imaging information for a tissue, mixture, or substance.
- Such processes and systems are useful for characterizing unknown material such as (among others) pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue. They are also useful for distinguishing tissue, or portions of tissue from other portions of tissue, including, for example, fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
- In some embodiments, processes in step 3 in the embodiment stated above can occur at a time (and/or using a CT scanner) different than step 4 in that embodiment; thus, some or all of step 3, that is, characterizing the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets can occur at a time (and/or using a CT scanner) different from some or all of step 4, that is from imaging information related to an unknown material received from the CT scanner at the first energy level; extracting at least one unknown material CT pixel value representative of a predetermined region of the unknown material; and presenting the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the reference material relative to the first basis material and the second basis material.
- According to one embodiment, which may be carried out using a system such as shown in
FIG. 1 described above and according to a process such as shown inFIG. 2 , a reference material is provided in the form of CT calibration phantom marketed by Mindways Software, Austin, Tex., with known density and mass attenuation coefficient. Water and K2HPO4 are used as basis materials, and the XCOM database of National Institute of Standards and Technology is used for obtaining mass attenuation coefficients for the two basis materials (http//physics.nist.gov/PhysRefData/Xcom/html/xcom1.html). For purposes of characterizing the response of the CT scanner for the reference material (here, human fat where the scanner is operated at an energy level of 80 kVp) relative to the first basis material and the second basis material, a linear response model is used as follows: -
CT Pixel Value=(Water Equiv Density)*A+(K2HPO4 Equiv Density)*B+C (Eq 8) - “Water Equivalent Density” and “K2HPO4 Equivalent Density” are the two numbers that describe the basis material composition for some material. A linear response model is not necessary, but appears to be sufficient for many practical applications of interest relative to CT imaging of human subjects. It is not necessary, but it appears to be sufficient to assume the coefficient “A” has a value of (1 HU/(mg/cm3)). This is essentially what would be expected for a CT scanner properly calibrated relative to the Hounsfield Unit scale. It is not necessary, but in many circumstances it is reasonable to assume the offset coefficient “C” has a value of (−1000 HU). As such, in their simplest form, certain embodiments of the invention can be reduced to imaging a single reference standard to gather sufficient information to solve a single equation in a single unknown (“B”). Multiple reference materials can be used leading to a set of equations in one or more unknowns that can be solved for the unknowns using various methods as long as the problem is not underdetermined (e.g., trying to use one reference material and one equation to extract 3 unknowns (e.g., “A”, “B” and “C”) is not going to lead to a useful representation of the CT scanner response in the context of this invention).
- Thus, according to this embodiment, CT pixel value measurement may be represented in a pixel value space defined by a tool used for deriving basis material composition estimates from what amounts to physical principles rather than from, say, x-ray attenuation measurements. For example, assume a linear response model is suitable and that A, B and C in Eq 8 are known. In Eq 8, therefore, “CT Pixel Value” now refers to a measurement from a CT image, and “Water Equivalent Density” and “K2HPO4 Equivalent Density” are now unknowns. Solving Eq. 8 for “K2HPO4 Equivalent Density”:
-
(K2HPO4 Equiv Density)=((Water Equiv Density)*A+C−(CT Pixel Value))/B (Eq 9) - That is, given a single pixel value measurement from a CT scanner, a line can be defined in the basis material space. The “Fat (80 kVp)” graph shown in
FIG. 3 , which can be generated and presented by the system such as shown inFIG. 1 , provides an example. It relates to a measurement from a CT scan from a region of subcutaneous fat at an energy level of 80 kVp. Whether or not the single value by itself is remarkable,FIG. 3 shows that the same information can be expressed or presented as a line (310). Additionally, the “Fat (80 kVp)” line can be presented if desired with one or more lines relating to reference materials.FIG. 4 shows “Fat (80 kVp)”line 410 withline 412 representing the reference material (here, a bone mineral density calibration standard or phantom marketed by Mindways Software, Austin, Tex.). This bone density calibration standard is based upon attenuation properties of aqueous solutions of K2HPO4. -
FIG. 4 shows theline 410 in the basis material space relative to the bone mineral densitystandard line 412—both are calculated at 80 kVp based upon the measured response of the CT scanner at that energy level. The intersection point for these two lines is the equivalent bone mineral density of the fat measurement. Thus, the space of possible basis material compositions consistent with the measured CT value for the subcutaneous fat measurement can be compared with a composition consistent with the bone mineral density calibration standard. The bone mineral density at the intersection with the K2HPO4 line is −69.4 mg/cm3 which is the same as the K2HPO4 density at the intersection point. - In another embodiment, a second CT scan is obtained of the subcutaneous fat at a second energy level of 120 kVp. The CT scanner response is different at 120 kVp, so the fat measurement line at 120 kVp is different. The bone mineral density standard is, however, a fixed standard in this space.
FIG. 5 , which can be generated and presented using systems such as shown inFIG. 1 and according to a process such as shown inFIG. 2 , shows lines for the fat measurements at both 80 kVp and 120 kVp, together with the bone mineral density standard. The water equivalent density scale has been expanded from that scale as shown inFIG. 4 . The fat measurement line at 80 kVp is shown asline 510; the bone mineral density standard line as 512, and the 120 kVp fat measurement line as 514. Note that the intersection of the Fat (120 kVp)measurement line 514 and the bone densitystandard line 512 is not the same as the point at which the Fat (80 kVp)measurement line 510 intersects the bone mineral densitystandard line 512. The equivalent bone mineral density of the 120 kVp subcutaneous fat measurement is about −87.5 mg/cm3. In fact the bone mineral density standard is not particularly convenient for characterizing fat because the bone mineral density estimate will change depending upon kVp and beam hardening (which affects the average energy of the x-rays hitting the CT scanner detector—beam hardening can be intentional, such as pre-filtration of the x-ray source or coincidental, such as degree of hardening due to patient size). - Note that the bone mineral density standard line did not change as a result of changing how the CT scanner is used. The density, for example, of an aqueous solution of 100 mg/cm3 of K2HPO4 does not change depending upon whether measurement occurs at 80 kVp or 120 kVp. It will always be that composition that is at the end of the black BMD calibration standard line which intersects the horizontal line corresponding to a K2HPO4 equivalent density of 100 mg/cm3. The unknown subcutaneous fat has the same physical attribute. Its density is not changed according to whether it is measured at 120 kVp or 80 kVp. The basis material composition for the subcutaneous fat that is consistent across the two measurements is at the intersection point of the 80 kVp and 120 kVp fat lines.
- In another embodiment, presentation of results for which are reflected in
FIG. 6 , which can be generated and presented using systems such as shown inFIG. 1 and according to a process such as shown inFIG. 2 , estimates of the bone mineral density calibration standard for two hypothetical scenarios can also be considered: (1) bone with a marrow space filled with pure yellow marrow and (2) bone with the marrow space filled with pure red marrow. Here, the fat 80 kVp line bearsnumeral 610; the bonemineral density line 612; the fat 120 kVpline 614; theyellow marrow line 616 and thered marrow line 618. These calculations are based upon atomic composition and density estimates from the literature, and are scanner independent and do not rely on a material that mimics bone. That is, there is no attempt to construct and measure with x-rays some yellow marrow mimic or some red marrow mimic. Instead, the lines for these two scenarios are constructed and presented based on mass attenuation coefficients obtained from the literature. - In another embodiment, the line generated and/or presented in the basis material space, which can be generated and presented using systems such as shown in
FIG. 1 and according to a process such as shown inFIG. 2 , can be compared to a standard reference line such as a bone mineral density line. The comparison can take the form of an intersection of the two lines, or shortest distance between the two lines. For example, the point could be estimated where “normal” liver sits in the basis material space. That point can be determined by experiment, such as the subcutaneous fat basis material composition point addressed above, or by calculation using known mass attenuation coefficients. Then, the closest distance between the normal liver reference point and the experimental measurement line can be considered. As another option, a line can be generated from the normal liver point based upon replacing liver with fat to calculate amount of fat based on where the experimental measurement line crosses the liver/fat line. As yet another option, the CT measurement can be dual energy measurements, in which case a basis material composition can be estimated for the unknown tissue. That estimate can allow measurements relative to the basis material composition rather than relative to the shortest distance between lines/points/curves. In similar or analogous fashion, two (or more) components or constituent materials of a region of interest or measurement volume may be characterized, such as according to the process shown inFIG. 2 , in (for example and without limitation) fat in pancreatic tissue, yellow/red bone marrow tissue, white-matter/grey-matter brain tissue, fat in muscle tissue, healthy/necrotic heart tissue, any healthy/necrotic soft tissue, any healthy or normal/diseased or pathological tissue, normal/diseased lung tissue, calcified stones/kidney tissue, and iron in liver tissue. Another potential application of DEQCT quantification methods according to certain embodiments of the invention is characterization of “bright” deposits. - For example, there is provided according to one embodiment a process for characterization of kidney stones. The process can be conducted using systems such as shown in
FIG. 1 and according to a process such as shown inFIG. 2 . Kidney stones may be predominantly composed of uric acid or calcium. The origins of these different classes of kidney stones are different and may have implications regarding treatment. Deposition of uric acid crystals in joint spaces (gout) is also characteristically different than other forms of arthritis that lead to growth of calcium-based osteophytes. Uric acid and calcium have distinctly different basis material compositions. - While there are conventional proposed methods to characterize kidney stones, such methods are CT-scanner specific methods—the results are specific to a particular machine. By contrast, processes according to these embodiments of the invention are not CT-scanner specific, because they make use of the basis material space to reduce or eliminate scanner dependence. These embodiments thus help overcome problems arising from, among other things, the fact that HU scales for non-water-like materials can vary significantly across CT scanners.
- The processes of these embodiments make use of the fact that in vivo kidney stones (and certain other structures) will be surrounded by water-like tissues. The kidney stones themselves will have appreciable amounts of calcium in the case of hydroxyapatite and oxalate stones, while cystine has sulfur and uric acid is nitrogen rich. Each of these materials has fairly different characteristics when viewed in the basis material space due to the significant difference in atomic numbers of these species. Uric acid to some extent looks like water (in the sense of being parallel to the water axis), but is actually slightly below the water axis reflecting the effective atomic number for uric acid being slightly less than water. Cystine has a pronounced positive “K2HPO4” character due to sulfur (atomic number 16—double that of oxygen). The calcium containing materials (Ca atomic number 20) are still higher, in the range of 2000 mg/cm3 to 3000 mg/cm3 range for K2HPO4 component.
- It can be difficult to obtain a CT image of just a kidney stone in vivo due to the spatial resolution limits of CT for reasonable clinical radiation doses. As such, what would be measured is some water-like tissue background along with a component characteristic of the composition of the kidney stone. In
FIG. 7 , which can be generated and presented using systems such as shown inFIG. 1 and according to a process such as shown inFIG. 2 , and which shows lines that relate to various potential reference materials relevant to kidney stones, the intersection point of all of the reference lines is the “water” point. Some moving to the right (and downward) along the “uric acid”reference line 710 would represent a mixture that looks mostly like water with some of that water-like material replaced with uric acid crystals. Similarly, with a kidney stone rich in calcium, there will be displacement somewhere along the calcium hydroxyapatite and calcium oxalate lines (712 and 714 respectively) (and perhaps well above the scale shown inFIG. 7 depending upon how much of the sampled volume is filled with stone versus water-like tissue). Because of the proximity of the calcium hydroxyapatite andcalcium oxalate lines line 716, is not nearly as strong an absorber as calcium hydroxyapatite or calcium oxalate, so perhaps it can be distinguished both based upon intensity and proximity to the cystine reference line. In any event, it is possible to readily distinguish, using processes according to these embodiments, uric acid kidney stones from kidney stones rich in calcium salts and/or cysteine. - It is possible to prepare a measurement classifier that could include factors such as estimated stone size, signal strength (basically how far up the K2HPO4 axis a particular measurement falls for the size of the stone relative to the sampled volume), as well as proximity to one of the reference lines above to identify the nature of the kidney stone (calcium/cystine/uric acid). The method of using bisectors as described in A.N. Primak, et al., Noninvasive Differentiation of Uric Acid versus Non-Uric Acid Kidney Stones Using Dual-Energy CT, Academic Radiology, p. 1441, Vol. 14, No. 12, December 2007 (which is incorporated herein by this reference), is one logical starting point. Bisectors for calcium/cystine as well as cystine/uric acid should be examined as well.
- In general, some embodiments relate to the projection of a single value received from a CT scan back into a two dimensional space, which single value defines a line in the two dimensional (basis material) space. The two dimensional space is very rich because, among other things, various theoretical and experimental results in the two dimensional space can be generated and explored in a controlled (calibrated) manner without significant consideration for the source of the experimental data (i.e., what CT scanner was used and in exactly what mode it was used). Moreover, the one-dimensional to two dimensional mapping can be reversed: Once something interesting is apparent in the device-independent two dimensional space, it can be transformed back into one dimension by applying the reverse transform derived for a specific CT scanner used in a particular manner. Such a transformation can be valuable, for instance, where there are various calibration problems that occur because by pixel values observed on a particular CT scanner depend upon the CT scanner and how it is used. One such example is above, where the subcutaneous fat example used above resulted in estimated CT values of −108.4 HU at 80 kVp and −83.7 HU at 120 kVp.
- Another embodiment of the invention involves comparing fat measurements made with Magnetic Resonance Spectroscopy (MRS) and DEQCT. MRS is an accepted method for in vivo quantification of marrow fat. The inventor has found that marrow fat as measured by DEQCT can accurately predict MRS fat measurements.
- Several issues immediately arise when considering the correlation of DEQCT and MRS derived estimates of marrow fat. MRS results are presented as a fat ratio. In particular, the ratio of (fat signal) to (fat signal+water signal) is reported with MRS when quantifying marrow fat. This ratio is expected to approach unity if the marrow space is dominated by the presence of yellow marrow (fatty marrow) and it is expect to approach zero if there is little fat in the marrow space (as might be the case for mostly red marrow). So this is neither an absolute quantification of the amount of water present nor an absolute quantification of the amount of fat present, and the interpretation of the measured ratio is confounded by the presence of other materials in a tissue environment as highly diverse as cancellous bone. That is, there may be molecular species present in the bone (and bone marrow) environment that are neither accurately described as “fat” or “water” that might contaminate either or both “signals” and their ratio.
- Compare the above MRS scenario to what is measured by DEQCT. There is no molecular level selectivity with DEQCT. Any atomic species present that can scatter x-rays (which means all atomic species present) in the sample volume will contribute to the measured DEQCT signal. Looking at a ratio such as the amount of fat contributing to the DEQCT signal versus the amount of (fat+water) contributing to the DEQCT signal is an ill-defined concept. How is the signal from calcium, collagen, hemoglobin and other materials present in cancellous bone removed from the signal so that the desired MRS-equivalent ratio can be predicted? This is a different scenario than, say, contemplating measuring the amount of fat present in liver tissue in cases where one can reasonably assume that only “fat” and “normal liver” tissue are present in a sample that is being quantified. Again, bone presents a very diverse environment that is simply not easy to describe in its substantial entirety by the amounts of two constituents.
- One approach to this marrow fat problem in the DEQCT space is to consider projecting the DEQCT measurements onto a line that passes through the basis material composition expected for water (e.g., as derived using the XCOM data) as well as the basis material composition expected for pure yellow marrow (e.g., as derived by using the XCOM data). That projection of DEQCT results onto such a line can provide a simple means for mapping DEQCT measurements to the MRS measurement space is not obvious. It is possible, for example, to consider other lines such as projection onto a line perpendicular to the bone density reference standard which nominally characterizes the expected change in signal due to the displacement of water from a hydrated cadaveric bone sample by the presence of mineralized bone. By projecting measurements onto this orthogonal line, it can be possible to subtract out the “mineralized collagen” component of the DEQCT bone signal. Another option may be to project the DEQCT measurements onto lines of constant BMD for a single energy QCT measurement. These lines in the basis material space represent mixtures of atomic constituents that can lead to the same BMD estimate for single energy QCT at a particular energy even though they had different relative amounts of bone and marrow and marrow composition. It is possible, for example, to contemplate the replacement of some red marrow with yellow marrow (which by itself would decrease the estimate bone mineral density) while adding some additional mineralized bone to the characterized volume such that the overall BMD remained constant.
- What can separate consideration of these various lines that may be interesting projection axes in the basis material space is if interest in them is not driven by predicting how a signal might change as a result of, for example, replacing some “normal liver” tissue with some “fat”, but rather by trying to find a way to quantify a DEQCT measurement in a way that excludes or significantly reduces the impact of some “uninteresting” part of a measurement (e.g., mineralized collagen) from the part that might be “interesting” (e.g., some measure of the relative amounts of water and fat in the marrow space). In this latter scenario, one does not have a detailed model of what's going on in bone, such that use of XCOM would allow ready prediction of an outcome and quantification of some signal relative to this predicted outcome. Rather, there is instead a somewhat ill-defined MRS measurement with some complex DEQCT measurements. It is not obvious that one can accurately predict the MRS results from the DEQCT measurements.
- One approach to this issue is to pair MRS and DEQCT measurements. One can ask the general question about the line that DEQCT measurements could be projected on that would maximize the correlation between the projected DEQCT measurements and the MRS measurements. One can restrict this search to a search of all possible lines that go through the expected basis material composition for water since “water” is an important component of the MRS signal. An extension of this approach is to generalize the search to all possible lines in the basis material plane. However, that is actually a more complex problem (because of degenerate solutions associated with the fact that whatever that “minimum” line is (and it may actually be more than one unique line) could be described relative to any point along the lines—all such descriptions would give the same value for the maximum correlation coefficient).
- The results of this exercise, which can be conducted using systems such as shown in
FIG. 1 and according to a process such as shown inFIG. 2 , are summarized inFIGS. 8 and 9 . They provide a potential example of possible utility of this approach to quantifying the results.FIG. 8 shows pair MR and DEQCT fat measurements. The square of the correlation coefficient (R) in this case is about 0.7 which implies R is about 0.84, which is a relatively strong correlation. -
FIG. 9 shows the change in the fat signal as measured by both methods. As mentioned in the comments above, both the DEQCT and MRS measurements are potentially contaminated by materials unrelated to those that are really the subject of quantification. Subtracting the baseline measurement from a follow-up measurement might subtract out a substantial portion of the signal contamination if the degree of contamination does not change appreciably between baseline and follow-up. When one considers the difference in signals, it is apparent that R in this case for this data is about 0.93 (R2 of about 0.86). As such, it may be that DEQCT predicts the change in MRS fat signal as accurately as or more accurately than it does in predicting the absolute MRS fat signal. - Coincidentally, it so happens that the projection line that maximizes R goes through the predicted water basis material composition (by construction) and it comes very close to going through the pure yellow marrow line. That is, projecting the DEQCT results on the line passing through yellow marrow and water (one of the options discussed above) seems to do a good job of isolating the “fat” and “water” components measured by MRS.
- e. General
- Numerous specific details are set forth herein to provide a thorough understanding of the claimed subject matter. However, those skilled in the art will understand that the claimed subject matter may be practiced without these specific details. In other instances, methods, apparatuses or systems that would be known by one of ordinary skill have not been described in detail so as not to obscure subject matter that may be claimed.
- The system or systems discussed herein are not limited to any particular hardware architecture or configuration. A computing device can include any suitable arrangement of components that provide a result conditioned on one or more inputs. Suitable computing devices include multipurpose microprocessor-based computer systems accessing stored software that programs or configures the computing system from a general-purpose computing apparatus to a specialized computing apparatus implementing one or more embodiments of the present subject matter. Any suitable programming, scripting, or other type of language or combinations of languages may be used to implement the teachings contained herein in software to be used in programming or configuring a computing device.
- Embodiments of the methods disclosed herein may be performed in the operation of such computing devices. The order of the blocks presented in the examples above can be varied—for example, blocks can be re-ordered, combined, and/or broken into sub-blocks. Certain blocks or processes can be performed in parallel.
- While the present subject matter has been described in detail with respect to specific embodiments thereof, it will be appreciated that those skilled in the art, upon attaining an understanding of the foregoing may readily produce alterations to, variations of, and equivalents to such embodiments. Accordingly, it should be understood that the present disclosure has been presented for purposes of example rather than limitation, and does not preclude inclusion of such modifications, variations and/or additions to the present subject matter as would be readily apparent to one of ordinary skill in the art.
Claims (26)
1. A method for execution by a computing device that comprises:
a. For a reference material and for each of a first basis material and a second basis material, receiving in a computing device an energy-dependent, mass attenuation coefficient for the material as a function of its elemental composition, molecular composition or elemental and molecular composition;
b. Deriving by the computing device a basis material composition of the reference material relative to the first basis material and the second basis material;
c. Characterizing by the computing device the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets;
d. From imaging information related to an unknown material received from the CT scanner at the first energy level; extracting by the computing device at least one unknown material CT pixel value representative of a predetermined region of the unknown material; and presenting by the computing device the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the reference material relative to the first basis material and the second basis material.
2. A method according to claim 1 further comprising using one or more CT pixel values to present an unknown material line in two dimensional basis material space, in order to characterize the unknown material.
3. A method according to claim 1 further comprising mapping a characterization of the response of the CT scanner to the first basis material and the second basis material to a single dimensional value relating to the reference material.
4. A method according to claim 1 further comprising mapping the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the first basis material and the second basis material to a single dimensional value relating to the reference material.
5. A method according to claim 1 further comprising comparing an unknown material line in two dimensional basis material space to a line for the reference material in two dimensional basis material space in order further to characterize the unknown material.
6. A method according to claim 1 further comprising receiving imaging information and presenting information for the unknown material at a first energy level and a second energy level.
7. A method according to claim 1 further comprising mapping the at least one unknown material CT pixel value back to a predetermined family of possible basis material compositions that exhibit the same CT pixel value as the at least one unknown material CT pixel value.
8. A method according to claim 1 wherein receiving imaging information for the unknown material includes receiving imaging information for a tissue, a mixture, or a substance.
9. A method according to claim 1 wherein presenting by the computing device the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the reference material relative to the first basis material and the second basis material includes presenting graphically.
10. A method according to claim 1 wherein the unknown material is selected from the group consisting of: pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue.
11. A method according to claim 1 further comprising distinguishing tissue selected from the group consisting of: fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
12. A system, comprising:
A processor device; and
A tangible medium embodying code that is executable by the processor device to cause the system to:
a. Receive, for a reference material and for each of a first basis material and a second basis material, an energy-dependent, mass attenuation coefficient for the material as a function of its elemental composition, molecular composition or elemental and molecular composition;
b. Derive a basis material composition of the reference material relative to the first basis material and the second basis material;
c. Characterize the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets;
d. From imaging information related to an unknown material received from the CT scanner at the first energy level, extract at least one unknown material CT pixel value representative of a predetermined region of the unknown material; and present the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the reference material relative to the first basis material and the second basis material.
13. A system according to claim 12 wherein the code is capable of causing the system to use one or more CT pixel values to present an unknown material line in two dimensional basis material space, in order to characterize the unknown material.
14. A system according to claim 12 wherein the code is capable of causing the system to map a characterization of the response of the CT scanner to the first basis material and the second basis material to a single dimensional value relating to the reference material.
15. A system according to claim 12 wherein the code is capable of causing the system to map the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the first basis material and the second basis material to a single dimensional value relating to the reference material.
16. A system according to claim 12 wherein the code is capable of causing the system to compare an unknown material line in two dimensional basis material space to a line for the reference material in two dimensional basis material space in order further to characterize the unknown material.
17. A system according to claim 12 wherein the code is capable of permitting the system to receive imaging information and present information for the unknown material at a first energy level and a second energy level.
18. A system according to claim 12 wherein the code is capable of causing the system to map the at least one unknown material CT pixel value back to a predetermined family of possible basis material compositions that exhibit the same CT pixel value as the at least one unknown material CT pixel value.
19. A method according to claim 12 wherein the code is configured to cause the system to characterize an unknown material selected from the group consisting of: pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue.
20. A method according to claim 12 further comprising code that is executable by the processor device to cause the system to distinguish tissue selected from the group consisting of: fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
21. A tangible computer readable medium that includes program code that is executable by a processor to perform actions, the program code comprising:
a. Code for causing the processor to receive, for a reference material and for each of a first basis material and a second basis material, an energy-dependent, mass attenuation coefficient for the material as a function of its elemental composition, molecular composition or elemental and molecular composition;
b. Code for causing the processor to derive a basis material composition of the reference material relative to the first basis material and the second basis material;
c. Code for causing the processor to characterize the response of a CT scanner to the reference material relative to the first basis material and the second basis material at a first energy level by receiving imaging information for the reference material from the CT scanner at the first energy level; extracting at least one reference material CT pixel value representative of the imaged reference material; describing the at least one reference material CT pixel value by a set of one or more equations predicting the at least one reference material CT pixel value as a function of the basis material composition of the reference material and one or more basis material weights and offsets; and solving the set of one or more equations for the basis material weights and offsets;
d. Code for causing the processor to do the following: From imaging information related to an unknown material received from the CT scanner at the first energy level; extract at least one unknown material CT pixel value representative of a predetermined region of the unknown material; and present the at least one unknown material CT pixel value as a function of the characterization of the response of the CT scanner to the reference material relative to the first basis material and the second basis material.
22. A tangible computer readable medium according to claim 21 wherein the program code includes code for causing the system to use one or more CT pixel values to present an unknown material line in two dimensional basis material space, in order to characterize the unknown material.
23. A tangible computer readable medium according to claim 21 wherein the program code includes code for causing the system to map a characterization of the response of the CT scanner to the first basis material and the second basis material to a single dimensional value relating to the reference material.
24. A tangible computer readable medium according to claim 21 wherein the program code includes code for permitting the system to receive imaging information and present information for the unknown material at a first energy level and a second energy level.
25. A method according to claim 21 wherein the code is configured to cause the system to characterize an unknown material selected from the group consisting of: pancreatic tissue, bone marrow tissue, brain tissue, fat in muscle tissue, heart tissue, soft tissue, non-bone tissue, lung tissue, kidney stones, kidney tissue, and iron in liver tissue.
26. A method according to claim 21 further comprising code that is executable by the processor to cause the system to distinguish tissue selected from the group consisting of: fat in pancreatic tissue, red and yellow bone marrow tissue, white and grey matter brain tissue, fat in muscle tissue, healthy and necrotic heart tissue, healthy and nectrotic soft tissue, healthy and diseased non-bone tissue, normal and diseased lung tissue, kidney stones in kidney tissue, and iron in liver tissue.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/770,617 US20140236488A1 (en) | 2013-02-19 | 2013-02-19 | Computed Tomography Calibration Systems and Methods |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/770,617 US20140236488A1 (en) | 2013-02-19 | 2013-02-19 | Computed Tomography Calibration Systems and Methods |
Publications (1)
Publication Number | Publication Date |
---|---|
US20140236488A1 true US20140236488A1 (en) | 2014-08-21 |
Family
ID=51351851
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/770,617 Abandoned US20140236488A1 (en) | 2013-02-19 | 2013-02-19 | Computed Tomography Calibration Systems and Methods |
Country Status (1)
Country | Link |
---|---|
US (1) | US20140236488A1 (en) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150031992A1 (en) * | 2013-07-29 | 2015-01-29 | The Cleveland Clinic Foundation | Identifying kidney stone composition from medical imaging |
US20170039735A1 (en) * | 2015-08-06 | 2017-02-09 | General Electric Company | Computed tomography self-calibration without calibration targets |
CN107072627A (en) * | 2014-10-20 | 2017-08-18 | 皇家飞利浦有限公司 | For the cardiac reconstruction for heart and the photon counting CT of lung image |
US20180082447A1 (en) * | 2014-08-16 | 2018-03-22 | Fei Company | Tomographic reconstruction for material characterization |
CN109959669A (en) * | 2017-12-26 | 2019-07-02 | 清华大学 | Sill decomposition method and device and computer readable storage medium |
EP3552547A1 (en) * | 2018-04-13 | 2019-10-16 | Siemens Healthcare GmbH | Method for providing conversion information about an imaging data set, x-ray device, computer program and electronically readable data carrier |
US10736601B2 (en) * | 2013-06-21 | 2020-08-11 | O.N.Diagnostics, LLC | Quantitative phantomless calibration of computed tomography scans |
DE102019110648A1 (en) * | 2019-04-25 | 2020-10-29 | China Medical University | Method for converting scan information from computed tomography scanners into bone parameters |
CN111881910A (en) * | 2020-07-31 | 2020-11-03 | 杭州依图医疗技术有限公司 | Information processing method based on vertebra image and computer readable storage medium |
CN112244884A (en) * | 2020-10-27 | 2021-01-22 | 沈阳先进医疗设备技术孵化中心有限公司 | Bone image acquisition method and device, console equipment and CT system |
CN113100803A (en) * | 2021-04-20 | 2021-07-13 | 西门子数字医疗科技(上海)有限公司 | Method, apparatus, computer device and medium for displaying venous thrombosis |
US11291416B2 (en) * | 2017-08-10 | 2022-04-05 | Fujifilm Healthcare Corporation | Parameter estimation method and X-ray CT system |
CN118383800A (en) * | 2024-06-26 | 2024-07-26 | 大庆君安医疗器械有限公司 | Intelligent ultrasonic inspection method and system for same-direction operation |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090207967A1 (en) * | 2008-02-15 | 2009-08-20 | Xin Liu | System and method for quantitative imaging of chemical composition to decompose multiple materials |
US20100223016A1 (en) * | 2007-08-17 | 2010-09-02 | Gary Gibson | Method and apparatus for identification and detection of liquids |
-
2013
- 2013-02-19 US US13/770,617 patent/US20140236488A1/en not_active Abandoned
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100223016A1 (en) * | 2007-08-17 | 2010-09-02 | Gary Gibson | Method and apparatus for identification and detection of liquids |
US20090207967A1 (en) * | 2008-02-15 | 2009-08-20 | Xin Liu | System and method for quantitative imaging of chemical composition to decompose multiple materials |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10736601B2 (en) * | 2013-06-21 | 2020-08-11 | O.N.Diagnostics, LLC | Quantitative phantomless calibration of computed tomography scans |
US9706970B2 (en) * | 2013-07-29 | 2017-07-18 | The Cleveland Clinic Foundation | Identifying kidney stone composition from medical imaging |
US20150031992A1 (en) * | 2013-07-29 | 2015-01-29 | The Cleveland Clinic Foundation | Identifying kidney stone composition from medical imaging |
US20180082447A1 (en) * | 2014-08-16 | 2018-03-22 | Fei Company | Tomographic reconstruction for material characterization |
US10311606B2 (en) | 2014-08-16 | 2019-06-04 | Fei Company | Correction of beam hardening artifacts in microtomography for samples imaged in containers |
US10354418B2 (en) * | 2014-08-16 | 2019-07-16 | Fei Company | Tomographic reconstruction for material characterization |
CN107072627A (en) * | 2014-10-20 | 2017-08-18 | 皇家飞利浦有限公司 | For the cardiac reconstruction for heart and the photon counting CT of lung image |
US20170039735A1 (en) * | 2015-08-06 | 2017-02-09 | General Electric Company | Computed tomography self-calibration without calibration targets |
US11291416B2 (en) * | 2017-08-10 | 2022-04-05 | Fujifilm Healthcare Corporation | Parameter estimation method and X-ray CT system |
CN109959669A (en) * | 2017-12-26 | 2019-07-02 | 清华大学 | Sill decomposition method and device and computer readable storage medium |
EP3552547A1 (en) * | 2018-04-13 | 2019-10-16 | Siemens Healthcare GmbH | Method for providing conversion information about an imaging data set, x-ray device, computer program and electronically readable data carrier |
CN110368006A (en) * | 2018-04-13 | 2019-10-25 | 西门子医疗有限公司 | For providing the method and X-ray equipment of the transitional information about image data set |
US11340364B2 (en) | 2018-04-13 | 2022-05-24 | Siemens Healthcare Gmbh | Providing an item of conversion information relating to an image dataset |
DE102019110648A1 (en) * | 2019-04-25 | 2020-10-29 | China Medical University | Method for converting scan information from computed tomography scanners into bone parameters |
CN111881910A (en) * | 2020-07-31 | 2020-11-03 | 杭州依图医疗技术有限公司 | Information processing method based on vertebra image and computer readable storage medium |
CN112244884A (en) * | 2020-10-27 | 2021-01-22 | 沈阳先进医疗设备技术孵化中心有限公司 | Bone image acquisition method and device, console equipment and CT system |
CN113100803A (en) * | 2021-04-20 | 2021-07-13 | 西门子数字医疗科技(上海)有限公司 | Method, apparatus, computer device and medium for displaying venous thrombosis |
CN118383800A (en) * | 2024-06-26 | 2024-07-26 | 大庆君安医疗器械有限公司 | Intelligent ultrasonic inspection method and system for same-direction operation |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20140236488A1 (en) | Computed Tomography Calibration Systems and Methods | |
Fuggle et al. | Fracture prediction, imaging and screening in osteoporosis | |
Adams et al. | Quantitative computer tomography in children and adolescents: the 2013 ISCD Pediatric Official Positions | |
Morgan et al. | Micro-computed tomography assessment of fracture healing: relationships among callus structure, composition, and mechanical function | |
Burghardt et al. | High-resolution computed tomography for clinical imaging of bone microarchitecture | |
Adams | Advances in bone imaging for osteoporosis | |
Treece et al. | Independent measurement of femoral cortical thickness and cortical bone density using clinical CT | |
MacNeil et al. | Improved reproducibility of high-resolution peripheral quantitative computed tomography for measurement of bone quality | |
MacNeil et al. | Load distribution and the predictive power of morphological indices in the distal radius and tibia by high resolution peripheral quantitative computed tomography | |
CN102048552B (en) | Beam hardening correction for ct perfusion measurements | |
Barnard et al. | Machine learning for automatic paraspinous muscle area and attenuation measures on low-dose chest CT scans | |
JP2009294209A (en) | Detection of material composition from calculation of effective atomic number | |
Saha et al. | Characterization of trabecular bone plate‐rod microarchitecture using multirow detector CT and the tensor scale: algorithms, validation, and applications to pilot human studies | |
JP6526669B2 (en) | Segmentation of bones from image data | |
Kim et al. | Reliability of skeletal muscle area measurement on CT with different parameters: a phantom study | |
Crabtree et al. | Dual-energy X-ray absorptiometry | |
US10360681B2 (en) | Method, a system and a computer program for determining data defining an estimate of the thickness and density of a cortical bone tissue structure of interest from imaging data | |
US11517278B2 (en) | System and method for basis material decomposition with general physical constraint for multi-energy computed tomography | |
Wang et al. | Characterisation of urinary stones in the presence of iodinated contrast medium using dual-energy CT: a phantom study | |
Li | Dual-energy CT with fast-kVp switching and its applications in orthopedics | |
McCollough et al. | Dual-energy algorithms and postprocessing techniques | |
Li et al. | Assessment of hepatic fatty infiltration using dual-energy computed tomography: a phantom study | |
Lala et al. | Measuring apparent trabecular structure with pQCT: a comparison with HR-pQCT | |
Lewiecki | Imaging technologies for assessment of skeletal health in men | |
Li et al. | Update on the clinical application of quantitative computed tomography (QCT) in osteoporosis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: MINDWAYS SOFTWARE, INC., TEXAS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:BROWN, J. KEENAN;REEL/FRAME:030146/0236 Effective date: 20130220 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |