CN116982995A - Device and method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT - Google Patents

Device and method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT Download PDF

Info

Publication number
CN116982995A
CN116982995A CN202310954800.6A CN202310954800A CN116982995A CN 116982995 A CN116982995 A CN 116982995A CN 202310954800 A CN202310954800 A CN 202310954800A CN 116982995 A CN116982995 A CN 116982995A
Authority
CN
China
Prior art keywords
transmission
fluorescence
signal intensity
reconstruction
scattering
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.)
Pending
Application number
CN202310954800.6A
Other languages
Chinese (zh)
Inventor
李亮
陈志强
张丽
高河伟
王振天
邢宇翔
邓智
宋佳丹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN202310954800.6A priority Critical patent/CN116982995A/en
Publication of CN116982995A publication Critical patent/CN116982995A/en
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/483Diagnostic techniques involving scattered radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/486Diagnostic techniques involving generating temporal series of image data
    • A61B6/487Diagnostic techniques involving generating temporal series of image data involving fluoroscopy
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • A61B6/5241Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT combining overlapping images of the same imaging modality, e.g. by stitching

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • Toxicology (AREA)
  • Physiology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The invention relates to the technical field of radiation imaging, in particular to a device and a method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT, wherein the device comprises: an X-ray source disposed on the rotary scanning electromechanical device; the transmission CT detector is arranged in the direction of an extension line of a connecting line of the X-ray source and the detected object so as to detect the signal intensity of an initial signal and a transmission photon; the fluorescent CT and scattering CT detector is arranged in the connecting direction perpendicular to the connecting line of the X-ray source and the detected object so as to detect scattered photons and fluorescent photon signal intensity; and the control and data processing equipment is respectively connected with the components to control the emission of the ray beam and the execution of the rotary CT scanning, and performs image reconstruction and attenuation correction according to the detected signal intensity to obtain an optimal transmission CT image and a fluorescence CT image. Therefore, the problems that the influence of scattered photons is ignored in a bimodal imaging scheme, the obtained fluorescent CT reconstruction result is inaccurate and the like are solved.

Description

Device and method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT
Technical Field
The invention relates to the technical field of radiation imaging, in particular to a device and a method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT.
Background
The X-ray transmission CT can realize structural imaging of human tissues by utilizing the characteristic that attenuation effects of different substances on X-rays are different, but because line attenuation coefficients of various soft tissues, blood, muscles and the like of a human body are very close, and under the influence of noise, dosage and the like in some cases, effective distinction cannot be carried out only by the X-ray transmission CT, so a bimodal imaging scheme of the X-ray transmission CT and the X-ray fluorescence CT (XFCT) is proposed in the related art.
The bimodal imaging scheme is that another modality XFCT is added on X-ray transmission CT, the XFCT can image specific tissues, nano particles prepared by adopting high Z elements (such as gold, gadolinium and the like) can be gathered at tumor tissues, the high Z element nano particles are injected into a body, and then the characteristic that characteristic X rays with specific energy can be generated by X-ray irradiation of the high Z element can be utilized, the high Z element distribution can be imaged, and then tumor tissue imaging is realized, so that the defect of X-ray transmission CT is overcome.
Further, in the bimodal imaging scheme, a filtered back projection algorithm is used for performing tomographic reconstruction of X-ray transmission CT, an MLEM algorithm based on a Poisson noise model is used for performing tomographic reconstruction of X-ray fluorescence CT, for the X-ray fluorescence CT, the influence of attenuation of incident photons in an object on a reconstruction result is required to be considered, and generally, the reconstruction result of transmission CT is used for performing attenuation correction on the fluorescence CT. However, in practice, the detector for X-ray fluorescence CT can collect not only fluorescent photons, but also scattered photons generated by compton scattering between X-rays and an object, where the scattered photons include electron density information, and the electron density information can be used for dose monitoring, and at the same time, the electron density information can also be used for image reconstruction with transmission CT, so as to obtain line attenuation coefficients under different energies for attenuation correction. Obviously, the influence of scattered photons is ignored in the bimodal imaging scheme, and the obtained fluorescent CT reconstruction result is inaccurate.
Disclosure of Invention
The invention provides a device and a method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT (computed tomography), which are used for solving the problems that the influence of scattered photons is ignored in a bimodal imaging scheme, the obtained fluorescence CT reconstruction result is inaccurate and the like.
An embodiment of a first aspect of the present invention provides an apparatus for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT, comprising: an X-ray source for emitting a ray beam to a detected object;
a transmission CT detector for detecting an initial signal intensity when the detected object is not placed and a transmission photon signal intensity emitted by the X-ray source and passing through the detected object;
at least one group of fluorescent CT and scattering CT detector for detecting the scattered photon signal intensity generated by Compton scattering and the fluorescent photon signal intensity generated by exciting the X-ray beam in the detected object.
Optionally, the transmission CT detector is disposed on a rotary scanning electromechanical device and is located in an extension line direction of a connection line between the X-ray source and the detected object.
Optionally, the at least one set of fluorescent CT and scattering CT detectors is disposed on a rotating scanner electromechanical device and is located in a connection direction perpendicular to a connection line between the X-ray source and the object to be detected.
Optionally, the method further comprises:
the at least one shielding box is arranged in the at least one shielding box in a one-to-one correspondence manner with the scattered CT detectors, and each shielding box is provided with a collimation hole.
Optionally, the method further comprises:
the control and data processing equipment is respectively connected with the rotary scanning electromechanical equipment, the X-ray source, the transmission CT detector and the at least one group of fluorescent CT and scattering CT detector so as to control the X-ray source to emit ray beams, control the rotary scanning electromechanical equipment to perform rotary CT scanning on the detected object, acquire the initial signal intensity, the transmission photon signal intensity, the scattering photon signal intensity and the fluorescent photon signal intensity, and perform image reconstruction and attenuation correction to obtain an optimal transmission CT image and a fluorescent CT image.
Optionally, the control and data processing device comprises:
the CT scanning control module is respectively connected with the X-ray source, the rotary scanning electromechanical device, the transmission CT detector and the at least one group of fluorescent CT and scattering CT detector to control the rotary scanning electromechanical device to rotate the X-ray source, the transmission CT detector and the at least one group of fluorescent CT and scattering CT detector so as to perform rotary CT scanning on the detected object;
the transmission/scattering/fluorescence CT data acquisition module is respectively connected with the transmission CT detector, the at least one group of fluorescence CT and the scattering CT detector to respectively acquire the initial signal intensity, the transmission photon signal intensity, the scattering photon signal intensity and the fluorescence photon signal intensity;
the transmission/scattering/fluorescence CT image reconstruction module is connected with the transmission/scattering/fluorescence CT data acquisition module so as to perform transmission reconstruction according to the initial signal intensity and the transmission photon signal intensity, perform fault reconstruction according to the scattering photon signal intensity and perform fluorescence reconstruction according to the fluorescence photon signal intensity;
the transmission/scattering/fluorescence CT image display module is connected with the transmission/scattering/fluorescence CT image reconstruction module so as to solve the accurate value of the line attenuation coefficient according to the transmission reconstruction result, the fault reconstruction result and the fluorescence reconstruction result, and carry out attenuation correction on the fault reconstruction result and the fluorescence reconstruction result according to the accurate value to obtain the optimal transmission CT image and the fluorescence CT image.
Optionally, the transmission/scattering/fluorescence CT image reconstruction module includes:
the transmission reconstruction unit is used for carrying out transmission reconstruction on the initial signal intensity and the transmission photon signal intensity by utilizing a filtering back projection algorithm to obtain a transmission reconstruction result;
the fault reconstruction unit is used for carrying out fault reconstruction on the scattered photon signal intensity by adopting an ML-EM algorithm based on poisson noise to obtain a fault reconstruction result;
and the fluorescence reconstruction unit is used for performing fluorescence reconstruction by utilizing the fluorescence photon signal intensity based on the ML-EM algorithm of poisson noise to obtain a fluorescence reconstruction result.
Optionally, the transmission reconstruction unit comprises:
a comparison subunit, configured to compare the initial signal intensity with the transmitted photon signal intensity, and solve a logarithm of the comparison result;
and the processing subunit is used for processing the logarithm by adopting a filtering back projection algorithm to obtain the transmission reconstruction result.
An embodiment of the second aspect of the present invention provides a method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT, comprising the steps of: controlling an X-ray source to emit a wire harness, and acquiring initial signal intensity when no detected object exists by utilizing a transmission CT detector; enabling the detected object to enter a rotary scanning electromechanical device, controlling an X-ray source to emit a wire harness, and enabling the rotary scanning electromechanical device to perform rotary CT scanning on the detected object; detecting transmitted photon signal intensity, scattered photon signal intensity and fluorescence photon signal intensity respectively by using the transmitted CT detector and at least one group of fluorescence CT and scattering CT detectors; performing image reconstruction according to the initial signal intensity, the transmission photon signal intensity, the scattered photon signal intensity and the fluorescence photon signal intensity to obtain a transmission reconstruction result, a fault reconstruction result and a fluorescence reconstruction result; and performing attenuation correction by using the transmission reconstruction result, the fault reconstruction result and the fluorescence reconstruction result to obtain an optimal transmission CT image and a fluorescence CT image.
Optionally, the performing image reconstruction according to the initial signal intensity, the transmission photon signal intensity, the scattered photon signal intensity and the fluorescence photon signal intensity to obtain a transmission reconstruction result, a tomographic reconstruction result and a fluorescence reconstruction result includes:
performing transmission reconstruction on the initial signal intensity and the transmission photon signal intensity by using a filtering back projection algorithm to obtain a transmission reconstruction result;
performing fault reconstruction on the scattered photon signal intensity by adopting an ML-EM algorithm based on poisson noise to obtain a fault reconstruction result;
and performing fluorescence reconstruction by using the fluorescence photon signal intensity based on the ML-EM algorithm of poisson noise to obtain a fluorescence reconstruction result.
Optionally, the performing transmission reconstruction on the initial signal intensity and the transmission photon signal intensity by using a filtered back projection algorithm to obtain a transmission reconstruction result includes:
comparing the initial signal intensity with the transmitted photon signal intensity to solve a logarithm of the comparison;
and processing the logarithm by adopting a filtering back projection algorithm to obtain the transmission reconstruction result.
Optionally, performing attenuation correction by using the transmission reconstruction result, the tomographic reconstruction result and the fluorescence reconstruction result to obtain an optimal transmission CT image and a fluorescence CT image, including:
solving an accurate value of a linear attenuation coefficient according to the transmission reconstruction result, the fault reconstruction result and the fluorescence reconstruction result;
and carrying out attenuation correction on the fault reconstruction result and the fluorescence reconstruction result according to the accurate value to obtain the optimal transmission CT image and the optimal fluorescence CT image.
Compared with an X-ray transmission and fluorescence bimodal imaging system, the device and the method for simultaneously imaging the X-ray transmission, fluorescence and scattering multimode CT have the advantages that a scattering modal imaging is additionally arranged, an electron density reconstruction result is obtained, the electron density result is a supplement to the transmission CT and fluorescence CT result, and dose estimation can be carried out through electron density; meanwhile, the linear attenuation coefficients of the object under different energies can be solved by using the three-mode imaging result, so that the fluorescent CT reconstruction result is accurately subjected to attenuation correction, richer object information can be obtained, and more accurate quantitative imaging is realized.
Additional aspects and advantages of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
The foregoing and/or additional aspects and advantages of the invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings, in which:
FIG. 1 is a schematic diagram of an apparatus for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT according to an embodiment of the present invention;
FIG. 2 is a block schematic diagram of a control and data processing apparatus provided in accordance with an embodiment of the present invention;
FIG. 3 is a block diagram of a transmission/scattering/fluorescence CT image reconstruction module according to an embodiment of the present invention;
FIG. 4 is a flow chart of a method for simultaneous X-ray transmission, fluorescence and scattering multi-mode CT imaging according to an embodiment of the present invention.
Reference numerals illustrate: the device for simultaneous imaging of 100-X-ray transmission, fluorescence and scattering multimode CT comprises a 1-X-ray source, a 2-rotary scanning electromechanical device, a 3-transmission CT detector, 4-at least one group of fluorescence CT and scattering CT detector, a 5-control and data processing device, a 51-CT scanning control module, a 52-transmission/scattering/fluorescence CT data acquisition module, a 53-transmission/scattering/fluorescence CT image reconstruction module, a 531-transmission reconstruction unit, a 532-tomographic reconstruction unit, a 533-fluorescence reconstruction unit, a 54-transmission/scattering/fluorescence CT image display module and 6-at least one shielding box.
Detailed Description
Embodiments of the present invention are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements or elements having like or similar functions throughout. The embodiments described below by referring to the drawings are illustrative and intended to explain the present invention and should not be construed as limiting the invention.
Fig. 1 is a schematic structural diagram of an apparatus for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT according to an embodiment of the present invention.
As shown in fig. 1, the apparatus 100 for simultaneous X-ray transmission, fluorescence and scattering multi-mode CT imaging includes: an X-ray source 1, a transmission CT detector 3, at least one set of fluorescent CT and scatter CT detectors 4, a control and data processing device 5.
Wherein the X-ray source 1 is arranged on a rotating scanner electromechanical device 2 for emitting a radiation beam towards the object to be examined. The transmission CT detector 3 is arranged on the rotary scanning electromechanical device 2 and is positioned in the direction of an extension line of a connecting line of the X-ray source and the detected object so as to detect the initial signal intensity when the detected object is not placed and the transmission photon signal intensity emitted by the X-ray source and passing through the detected object. At least one group of fluorescent CT and scattering CT detectors 4 are arranged on the rotary scanning electromechanical device and are positioned in the connecting direction perpendicular to the connecting line of the X-ray source and the detected object so as to detect the intensity of scattered photon signals generated by Compton scattering and the intensity of fluorescent photon signals generated by excitation of the ray beam emitted by the X-ray source 1 in the detected object. The control and data processing device 5 is respectively connected with the rotary scanning electromechanical device 2, the X-ray source 1, the transmission CT detector 3 and at least one group of fluorescence CT and scattering CT detector 4 to control the X-ray source 1 to emit ray beams, control the rotary scanning electromechanical device 2 to perform rotary CT scanning on a detected object, acquire initial signal intensity, transmission photon signal intensity, scattering photon signal intensity and fluorescence photon signal intensity, and perform image reconstruction and attenuation correction to obtain an optimal transmission CT image and a fluorescence CT image.
Further, the X-ray source 1 in the embodiment of the present invention may include: the common X-ray machine, accelerator and synchrotron radiation source of various types can also be a device such as a radioactive isotope which can emit X-rays or gamma-rays, and the device is arranged on the rotary scanning electromechanical device 2 to rotate and emit ray beams to a tested object according to the instruction of the control and data processing device 5, wherein the rotary scanning electromechanical device 2 can adopt a slip ring structure, can complete acquisition through the rotary X-ray machine and a detector, can also fix a scanned object on a turntable, and can realize scanning through the rotating object.
Further, the transmission CT detector 3 in the embodiment of the present invention may be a conventional energy integration detector, such as a multi-row detector array or a flat panel detector based on a scintillator+photoelectric conversion device, which is commonly used in the current medical spiral CT, and may also be a photon counting energy spectrum detector. The transmission CT detector 3 is arranged on the rotary scanning electromechanical device 2, the position of the transmission CT detector is positioned in the direction of an extension line of a connecting line of the X-ray source and the detected object, initial signal intensity when the detected object is not placed is collected firstly according to the instruction of the control and data processing device 5, and then transmission photon signal intensity generated by the X-ray source 1 after the ray beam passes through the detected object is collected.
Optionally, the apparatus 10 for simultaneous X-ray transmission, fluorescence and scatter multi-mode CT imaging further comprises
At least one shielding box 6, at least one group of fluorescent CT and scattering CT detectors 4 are arranged in the shielding box 6, and a collimation hole is arranged on the shielding box 6.
Specifically, each group of fluorescent CT and scattering CT detector 4 is respectively disposed in a shielding box 6, the shielding box 6 is disposed on a rotary scanning electromechanical device, the position is located in the direction perpendicular to the connection line between the X-ray source and the detected object, the fluorescent CT and scattering CT detector 4 detects the scattered photon signal intensity generated by compton scattering through a collimation hole on the shielding box 6 according to the instruction of the control and data processing device 5, and the fluorescent photon signal intensity generated by excitation of the radiation beam emitted by the X-ray source 1 in the detected object.
Optionally, as shown in fig. 2, the control and data processing device 5 comprises:
the CT scanning control module 51 is respectively connected with the X-ray source 1, the rotary scanning electromechanical device 2, at least one group of fluorescent CT and the scattering CT detector 4 to control the rotary scanning electromechanical device 2 to rotate the X-ray source 1, the transmission CT detector 3 and at least one group of fluorescent CT and the scattering CT detector 4 so as to perform rotary CT scanning on a detected object;
the transmission/scattering/fluorescence CT data acquisition module 52, the transmission/scattering/fluorescence CT data acquisition module 52 is respectively connected with the transmission CT detector 3 and at least one group of fluorescence CT and the scattering CT detector 4 to respectively acquire initial signal intensity, transmission photon signal intensity, scattered photon signal intensity and fluorescence photon signal intensity;
the transmission/scattering/fluorescence CT image reconstruction module 53, the transmission/scattering/fluorescence CT image reconstruction module 53 is connected with the transmission/scattering/fluorescence CT data acquisition module 52 to perform transmission reconstruction according to the initial signal intensity and the transmission photon signal intensity, perform tomographic reconstruction according to the scattering photon signal intensity, and perform fluorescence reconstruction according to the fluorescence photon signal intensity;
the transmission/scattering/fluorescence CT image display module 54, the transmission/scattering/fluorescence CT image display module 54 is connected with the transmission/scattering/fluorescence CT image reconstruction module 53, so as to solve the accurate value of the line attenuation coefficient according to the transmission reconstruction result, the tomographic reconstruction result and the fluorescence reconstruction result, and perform attenuation correction on the tomographic reconstruction result and the fluorescence reconstruction result according to the accurate value, so as to obtain an optimal transmission CT image and a fluorescence CT image.
Specifically, the embodiment of the invention uses the CT scanning control module 51 to control the rotary scanning electromechanical device 2 to rotate the X-ray source 1, the transmission CT detector 3 and at least one group of fluorescent CT and scattering CT detector 4, and controls the X-ray source 1 to emit a ray beam to the detected object so as to perform rotary CT scanning on the detected object;
the transmission/scattering/fluorescence CT data acquisition module 52 controls the transmission CT detector 3 to acquire initial signal intensity and transmission photon signal intensity, and controls at least one group of fluorescence CT and scattering CT detector 4 to acquire scattered photon signal intensity and fluorescence photon signal intensity;
the transmission/scattering/fluorescence CT image reconstruction module 53 performs transmission reconstruction according to the initial signal intensity and the transmission photon signal intensity to obtain a transmission reconstruction result, performs tomographic reconstruction according to the scattering photon signal intensity to obtain a tomographic reconstruction result (i.e., an electron density reconstruction result), and performs fluorescence reconstruction according to the fluorescence photon signal intensity to obtain a fluorescence reconstruction result;
the transmission/scattering/fluorescence CT image display module 54 solves the accurate values of the line attenuation coefficients according to the transmission reconstruction result, the tomographic reconstruction result and the fluorescence reconstruction result, and performs attenuation correction on the tomographic reconstruction result and the fluorescence reconstruction result according to the accurate values, so as to obtain an optimal transmission CT image and a fluorescence CT image.
Optionally, as shown in fig. 3, the transmission/scattering/fluorescence CT image reconstruction module 53 in the embodiment of the present invention includes:
the transmission reconstruction unit 531 is used for performing transmission reconstruction on the initial signal intensity and the transmission photon signal intensity by utilizing a filtering back projection algorithm to obtain a transmission reconstruction result;
a tomographic reconstruction unit 532, configured to perform tomographic reconstruction on the scattered photon signal intensity by using an ML-EM algorithm based on poisson noise, to obtain a tomographic reconstruction result;
and a fluorescence reconstruction unit 533 for performing fluorescence reconstruction by using the fluorescence photon signal intensity of the ML-EM algorithm based on poisson noise, so as to obtain a fluorescence reconstruction result.
Still further, the transmission reconstruction unit 531 in the embodiment of the present invention includes:
a comparison subunit for comparing the initial signal intensity with the transmitted photon signal intensity and solving the logarithm after comparison;
and the processing subunit is used for processing the logarithm by adopting a filtering back projection algorithm to obtain a transmission reconstruction result.
In summary, according to the device for simultaneously imaging the X-ray transmission, fluorescence and scattering multimode CT provided by the embodiment of the invention, a scattering mode imaging is newly added, an electron density reconstruction result is obtained, the electron density result is a supplement to the transmission CT and fluorescence CT result, and the dosage estimation can be carried out through the electron density; meanwhile, the linear attenuation coefficients of the object under different energies can be solved by using the three-mode imaging result, so that the fluorescent CT reconstruction result is accurately subjected to attenuation correction, richer object information can be obtained, and more accurate quantitative imaging is realized.
Next, a method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT according to an embodiment of the present invention will be described with reference to the accompanying drawings.
Fig. 4 is a flowchart of a method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT according to an embodiment of the present invention.
As shown in fig. 4, the method for simultaneous imaging of the X-ray transmission, fluorescence and scattering multimode CT comprises the following steps:
in step S101, the X-ray source is controlled to emit a beam, and the transmission CT detector is used to acquire an initial signal intensity without the detected object.
Specifically, signals of a transmission CT detector without a detected object are acquired:
I'=∫ E' I 0 (E')Q(E')dE' (1)
wherein I' is the initial signal intensity acquired by the transmission CT detector when no detected object exists, I 0 And (E ') the emergent energy spectrum of the X-ray machine, Q (E') is the energy response function of the transmission CT detector, and can be solved through experiments.
In step S102, the object is entered into the rotary scanning electro-mechanical device, the X-ray source is controlled to emit a beam, and the rotary scanning electro-mechanical device is caused to perform a rotary CT scan on the object.
In step S103, the transmitted photon signal intensity, the scattered photon signal intensity, and the fluorescence photon signal intensity are detected using the transmitted CT detector and at least one set of fluorescence CT and scattered CT detectors, respectively.
Specifically, after the X-ray generated by the conventional X-ray machine passes through an object, the X-ray is collected by a transmission CT (computed tomography) by using an array detector, and the signal intensity collected by the detector is as follows:
wherein I is 0 (E) The emergent energy spectrum of the X-ray machine can be calculated by simulation software,for scanned object at->The energy is the line attenuation coefficient under E, l is the path of X-ray passing through the object in the process from the optical machine to the detector, E is natural logarithm, Q (E) is the energy response function of the detector, and the energy response function can be solved through experiments.
Compton scattering can occur and scattered photons can be generated when X rays irradiate an object, and the scattered photon intensity (namely scattered photon signal intensity) of the energy Esca collected by the fluorescence CT and the scattered CT detector is as follows:
wherein I is ISCT (E sca ) To scatter photon intensity, E sca To scatter photon energy, l D Is the intersection line between the straight line connected with the fluorescent CT and scattering CT detector and the collimation hole and the object, E 0 To enter photon energy, I 0 (E 0 ) The energy emitted by the X-ray source is E 0 E is natural logarithm,for the detected object at->Energy at E 0 Lower line attenuation coefficient, ">For X-ray source S to position->Path traversed by the place, < >>Is the position/>Electron density at the site->For incident energy E 0 Compton scattering differential section with scattering angle θ, < >>For the detected object at->Energy at E sca Lower line attenuation coefficient, ">For the path of the fluorescent photons through the object to the detector, angle is the position +.>To a corresponding solid angle at the fluorescent CT and the scatter CT detector D.
Due to the limitation of the pinhole on the geometric relationship, the positionOnce given, the position of the detector unit D is determined, the solid Angle can be solved according to the geometrical relationship. The scattering angle θ is ST (i.e., the position S of the source of X-rays and the position of the scattering point of the object to be detected +.>Connecting line) and TD (position of scattering point of detected object +.>And the fluorescent CT and the scattered CT detector D), when the energy E of scattered photons sca And the energy E of the incident photon when the scattering angle theta is given 0 Can pass Compton scattering formulaAnd (5) solving. Solving E 0 After that, the emergent energy spectrum I of the X-ray machine 0 (E 0 ) Can be calculated by simulation software. Compton scattering differential section dSigma c,e (E 0 θ)/dΩ can be determined by Klein-Nishina formula.
According to the above analysis, I in equation (3) 0 (E 0 )·dσ c,e (E 0 θ)/dΩ·angle can be solved,unknown to be solved. Wait for the quantity->The solution can be performed using classical reconstruction algorithms in the CT field.
The object containing high Z element (gold or gadolinium, etc.) is irradiated by a conventional X-ray machine, and the intensity of fluorescence photon signals collected by the fluorescence CT and scattering CT detector is as follows:
wherein I is XFCT For fluorescence photon signal intensity, l D Is the intersection line between the straight line connected with the fluorescent CT and scattering CT detector and the collimation hole and the object, E K K-edge energy, E, of high Z element max For maximum energy of the X-ray source emergent energy spectrum, I 0 (E) The emergent energy spectrum of the X-ray source, e is natural logarithm,for the detected object at->The line attenuation coefficient at energy E, < ->For the beam to emerge from the X-ray source to the position +.>Path traversed by the place, < >>Is the photoelectric effect mass attenuation coefficient corresponding to the high Z element when the energy is E, +.>For position->At the concentration of high Z element, ω is the fluorescence yield of high Z element, +.>For the path of the fluorescence photons through the object to the light CT and scatter CT detector, +.>For the linear attenuation coefficient of the fluorescent photons in the object, angle is the position +.>To a corresponding solid angle at the fluorescent CT and the scatter CT detector D.
Emission spectrum I of X-ray machine 0 (E) The linear attenuation coefficients of the elements under different energies can be calculated through simulation softwareIt can be seen that the fluorescence yield ω of this element is a known constant and that the solid Angle can be found given the geometric relationship.
In step S104, image reconstruction is performed according to the initial signal intensity, the transmission photon signal intensity, the scattered photon signal intensity, and the fluorescence photon signal intensity, so as to obtain a transmission reconstruction result, a tomographic reconstruction result, and a fluorescence reconstruction result.
Further, step S104 specifically includes:
performing transmission reconstruction on the initial signal intensity and the transmission photon signal intensity by using a filtering back projection algorithm to obtain a transmission reconstruction result;
carrying out fault reconstruction on scattered photon signal intensity by adopting an ML-EM algorithm based on poisson noise to obtain a fault reconstruction result;
and performing fluorescence reconstruction by using the ML-EM algorithm fluorescence photon signal intensity based on poisson noise to obtain a fluorescence reconstruction result.
Specifically, the transmission reconstruction is performed on the initial signal intensity and the transmission photon signal intensity by using a filtering back projection algorithm, so as to obtain a transmission reconstruction result, which comprises the following steps:
comparing the initial signal intensity with the transmitted photon signal intensity to solve a logarithm of the comparison;
the logarithm after comparison is as follows:
wherein I' is the initial signal intensity, I is the transmitted photon signal intensity, I 0 (E) E is natural logarithm, l is the part of the X-ray and detector connecting line passing through the scanned object,for scanned object at->The energy is the linear attenuation coefficient under E, Q (E) is the energy response function of the transmission CT detector, I 0 (E ') is the emission spectrum of the X-ray source, Q (E') is the energy response function of the detector, w (E) is the energy spectrum weight function, generally according to I 0 (E) And Q (E) solving.
Adopts classical filtering back projection algorithm to make logarithmic pairsAnd processing to obtain a transmission reconstruction result.
Performing fault reconstruction on scattered photon signal intensity by using an ML-EM algorithm based on Poisson noise to obtain a fault reconstruction result, wherein the method comprises the following steps:
due to the number of scattered photons I in the scattered photon signal intensity ISCT (E sca ) Rarely, the scattered photon signal intensity is generally reconstructed by using an ML-EM algorithm based on poisson noise, and a fault reconstruction result is obtained.
Performing fluorescence reconstruction by using the ML-EM algorithm fluorescence photon signal intensity based on poisson noise to obtain a fluorescence reconstruction result, wherein the fluorescence reconstruction result comprises the following steps:
ignoring incident photon attenuation in fluorescence photon signal intensityAnd fluorescence photon attenuationPosition +.>Concentration of high Z element->And (5) reconstructing to obtain a fluorescence reconstruction result.
In step S105, attenuation correction is performed using the transmission reconstruction result, the tomographic reconstruction result, and the fluorescence reconstruction result, and an optimal transmission CT image and a fluorescence CT image are obtained.
Further, performing attenuation correction by using the transmission reconstruction result, the tomographic reconstruction result and the fluorescence reconstruction result to obtain an optimal transmission CT image and a fluorescence CT image, including:
solving the accurate value of the linear attenuation coefficient according to the transmission reconstruction result, the fault reconstruction result and the fluorescence reconstruction result;
and carrying out attenuation correction on the fault reconstruction result and the fluorescence reconstruction result according to the accurate value to obtain an optimal transmission CT image and a fluorescence CT image.
Specifically, the transmission reconstruction result, the fault reconstruction result and the fluorescence reconstruction result can be used for solving the line attenuation coefficientAfter the accurate value is obtained, the reconstruction results of the fluorescence CT and the scattering CT can be subjected to attenuation correction, and the specific calculation process is as follows:
the theoretical solution formula of the line attenuation coefficient is:
wherein, the liquid crystal display device comprises a liquid crystal display device,to be only in +.>Related unknowns, sigma c,e (E) For the total cross section of the single electron Compton effect, μ can be solved by the Klein-Nishina formula given E Z (E) For the line attenuation coefficient of the high Z element at energy E, mu when the high Z element and energy E are given Z (E) Can be checked.
The acquisition signals that bring the theory of the line attenuation coefficient into transmission CT can be obtained:
taking the values solved in (3) and (4) as valuesAnd->Then can be found by using the formula (7)For->Proceeding withFiltered back projection can be obtained +.>Then bring it into formula (6) to calculate +.>Then willBy taking in equation (5) and equation (3), a new +.>And->Of course, since the initial value is not accurate, the reconstructed new +.>And->The values of (2) are not completely accurate, so that the embodiment of the invention adopts an iterative mode to reconstruct new +.>And->Repeating the above steps as new initial value until the reconstructed +.>And->The values each converge to their true values.
It should be noted that the explanation of the embodiment of the apparatus for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT is also applicable to the method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT in this embodiment, and will not be repeated here.
According to the method for simultaneously imaging the X-ray transmission, fluorescence and scattering multimode CT provided by the embodiment of the invention, a scattering mode imaging is newly added, an electron density reconstruction result is obtained, the electron density result is a supplement to the transmission CT and fluorescence CT result, and the dosage estimation can be carried out through the electron density; meanwhile, the linear attenuation coefficients of the object under different energies can be solved by using the three-mode imaging result, so that the fluorescent CT reconstruction result is accurately subjected to attenuation correction, richer object information can be obtained, and more accurate quantitative imaging is realized.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In the description of the present invention, the meaning of "plurality" means at least two, for example, two, three, etc., unless specifically defined otherwise.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, the different embodiments or examples described in this specification and the features of the different embodiments or examples may be combined and combined by those skilled in the art without contradiction.
While embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the invention, and that variations, modifications, alternatives and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the invention.

Claims (12)

1. An apparatus for simultaneous X-ray transmission, fluorescence and scatter multi-mode CT imaging, comprising:
an X-ray source for emitting a ray beam to a detected object;
a transmission CT detector for detecting an initial signal intensity when the detected object is not placed and a transmission photon signal intensity emitted by the X-ray source and passing through the detected object;
at least one group of fluorescent CT and scattering CT detector for detecting the scattered photon signal intensity generated by Compton scattering and the fluorescent photon signal intensity generated by exciting the X-ray beam in the detected object.
2. The apparatus for simultaneous X-ray transmission, fluorescence and scattering multi-mode CT imaging of claim 1, wherein the transmission CT detector is disposed on a rotary scanner electromechanical device and is positioned in a direction of an extension line of a line connecting the X-ray source and the object under examination.
3. The apparatus for simultaneous X-ray transmission, fluorescence and scatter multi-mode CT imaging of claim 1, wherein the at least one set of fluorescence CT and scatter CT detectors is disposed on a rotating scanner electromechanical device and is positioned in a line direction perpendicular to a line connecting the X-ray source and the object under examination.
4. The X-ray transmission, fluorescence and scatter multi-mode CT simultaneous imaging apparatus of claim 1, further comprising:
the at least one shielding box is arranged in the at least one shielding box in a one-to-one correspondence manner with the scattered CT detectors, and each shielding box is provided with a collimation hole.
5. The X-ray transmission, fluorescence and scatter multi-mode CT simultaneous imaging apparatus of claim 1, further comprising:
the control and data processing equipment is respectively connected with the rotary scanning electromechanical equipment, the X-ray source, the transmission CT detector and the at least one group of fluorescent CT and scattering CT detector so as to control the X-ray source to emit ray beams, control the rotary scanning electromechanical equipment to perform rotary CT scanning on the detected object, acquire the initial signal intensity, the transmission photon signal intensity, the scattering photon signal intensity and the fluorescent photon signal intensity, and perform image reconstruction and attenuation correction to obtain an optimal transmission CT image and a fluorescent CT image.
6. The X-ray transmission, fluorescence and scattering multi-mode CT simultaneous imaging apparatus of claim 5, wherein said control and data processing device comprises:
the CT scanning control module is respectively connected with the X-ray source, the rotary scanning electromechanical device, the transmission CT detector and the at least one group of fluorescent CT and scattering CT detector to control the rotary scanning electromechanical device to rotate the X-ray source, the transmission CT detector and the at least one group of fluorescent CT and scattering CT detector so as to perform rotary CT scanning on the detected object;
the transmission/scattering/fluorescence CT data acquisition module is respectively connected with the transmission CT detector, the at least one group of fluorescence CT and the scattering CT detector to respectively acquire the initial signal intensity, the transmission photon signal intensity, the scattering photon signal intensity and the fluorescence photon signal intensity;
the transmission/scattering/fluorescence CT image reconstruction module is connected with the transmission/scattering/fluorescence CT data acquisition module so as to perform transmission reconstruction according to the initial signal intensity and the transmission photon signal intensity, perform fault reconstruction according to the scattering photon signal intensity and perform fluorescence reconstruction according to the fluorescence photon signal intensity;
the transmission/scattering/fluorescence CT image display module is connected with the transmission/scattering/fluorescence CT image reconstruction module so as to solve the accurate value of the line attenuation coefficient according to the transmission reconstruction result, the fault reconstruction result and the fluorescence reconstruction result, and carry out attenuation correction on the fault reconstruction result and the fluorescence reconstruction result according to the accurate value to obtain the optimal transmission CT image and the fluorescence CT image.
7. The X-ray transmission, fluorescence and scattering multi-mode CT simultaneous imaging apparatus of claim 6, wherein the transmission/scattering/fluorescence CT image reconstruction module comprises:
the transmission reconstruction unit is used for carrying out transmission reconstruction on the initial signal intensity and the transmission photon signal intensity by utilizing a filtering back projection algorithm to obtain a transmission reconstruction result;
the fault reconstruction unit is used for carrying out fault reconstruction on the scattered photon signal intensity by adopting an ML-EM algorithm based on poisson noise to obtain a fault reconstruction result;
and the fluorescence reconstruction unit is used for performing fluorescence reconstruction by utilizing the fluorescence photon signal intensity based on the ML-EM algorithm of poisson noise to obtain a fluorescence reconstruction result.
8. The apparatus for simultaneous X-ray transmission, fluorescence and scattering multi-mode CT imaging of claim 7, wherein said transmission reconstruction unit comprises:
a comparison subunit, configured to compare the initial signal intensity with the transmitted photon signal intensity, and solve a logarithm of the comparison result;
and the processing subunit is used for processing the logarithm by adopting a filtering back projection algorithm to obtain the transmission reconstruction result.
9. A method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT, characterized in that a device for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT according to any of claims 1-8 comprises the steps of:
controlling an X-ray source to emit a wire harness, and acquiring initial signal intensity when no detected object exists by utilizing a transmission CT detector;
enabling the detected object to enter a rotary scanning electromechanical device, controlling an X-ray source to emit a wire harness, and enabling the rotary scanning electromechanical device to perform rotary CT scanning on the detected object;
detecting transmitted photon signal intensity, scattered photon signal intensity and fluorescence photon signal intensity respectively by using the transmitted CT detector and at least one group of fluorescence CT and scattering CT detectors;
performing image reconstruction according to the initial signal intensity, the transmission photon signal intensity, the scattered photon signal intensity and the fluorescence photon signal intensity to obtain a transmission reconstruction result, a fault reconstruction result and a fluorescence reconstruction result;
and performing attenuation correction by using the transmission reconstruction result, the fault reconstruction result and the fluorescence reconstruction result to obtain an optimal transmission CT image and a fluorescence CT image.
10. The method of simultaneous X-ray transmission, fluorescence and scattering multi-mode CT of claim 9, wherein performing image reconstruction from the initial signal intensity, the transmission photon signal intensity, the scattering photon signal intensity, and the fluorescence photon signal intensity to obtain transmission reconstruction results, tomographic reconstruction results, and fluorescence reconstruction results comprises:
performing transmission reconstruction on the initial signal intensity and the transmission photon signal intensity by using a filtering back projection algorithm to obtain a transmission reconstruction result;
performing fault reconstruction on the scattered photon signal intensity by adopting an ML-EM algorithm based on poisson noise to obtain a fault reconstruction result;
and performing fluorescence reconstruction by using the fluorescence photon signal intensity based on the ML-EM algorithm of poisson noise to obtain a fluorescence reconstruction result.
11. The method of simultaneous X-ray transmission, fluorescence and scattering multi-mode CT imaging of claim 10, wherein said performing a transmission reconstruction of said initial signal intensity and said transmitted photon signal intensity using a filtered back-projection algorithm, obtaining said transmission reconstruction result, comprises:
comparing the initial signal intensity with the transmitted photon signal intensity to solve a logarithm of the comparison;
and processing the logarithm by adopting a filtering back projection algorithm to obtain the transmission reconstruction result.
12. The method of simultaneous X-ray transmission, fluorescence and scattering multi-mode CT imaging of claim 11, wherein performing attenuation correction using the transmission reconstruction result, the tomographic reconstruction result, and the fluorescence reconstruction result to obtain an optimal transmission CT image and a fluorescence CT image comprises:
solving an accurate value of a linear attenuation coefficient according to the transmission reconstruction result, the fault reconstruction result and the fluorescence reconstruction result;
and carrying out attenuation correction on the fault reconstruction result and the fluorescence reconstruction result according to the accurate value to obtain the optimal transmission CT image and the optimal fluorescence CT image.
CN202310954800.6A 2023-07-31 2023-07-31 Device and method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT Pending CN116982995A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310954800.6A CN116982995A (en) 2023-07-31 2023-07-31 Device and method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310954800.6A CN116982995A (en) 2023-07-31 2023-07-31 Device and method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT

Publications (1)

Publication Number Publication Date
CN116982995A true CN116982995A (en) 2023-11-03

Family

ID=88533288

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310954800.6A Pending CN116982995A (en) 2023-07-31 2023-07-31 Device and method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT

Country Status (1)

Country Link
CN (1) CN116982995A (en)

Similar Documents

Publication Publication Date Title
US8315352B2 (en) System and method of spectral calibration and basis material decomposition for X-ray CT systems
Jaffray et al. Cone‐beam computed tomography with a flat‐panel imager: initial performance characterization
EP2041606B1 (en) Energy spectrum reconstruction
EP2296550B1 (en) K-edge imaging
EP0124015B1 (en) Method and apparatus for reducing image artifacts due to projection measurement inconsistencies
JP5367574B2 (en) X-ray CT apparatus and method
CN107315019A (en) Transmission of radiation and fluorescence CT imaging systems and imaging method
US9320477B2 (en) Method and apparatus for adaptive scatter correction
WO2007074772A1 (en) X-ray ct device
JP5622487B2 (en) Radiation diagnostic apparatus and image reconstruction method
EP2313865B1 (en) System and method for spectral x-ray imaging
JP2009508616A (en) CT image forming system
IL144048A (en) Apparatus and methods and computer programs for estimating and correcting scatter in digital radiographic and tomographic imaging
JP6595154B2 (en) Diagnostic imaging apparatus using X-ray CT
JP6021347B2 (en) Medical image capturing apparatus and medical image capturing method
CN113237903B (en) Double-layer flat panel detector cone-beam CT-based scattering correction method and device
US20040210132A1 (en) Simulation of nuclear medical imaging
JP2010243395A (en) X ray-gamma ray imaging device
CN116982995A (en) Device and method for simultaneous imaging of X-ray transmission, fluorescence and scattering multimode CT
JP2011163966A (en) Medical imagery diagnostic apparatus and control program for radiation dose calculation
JPH08304309A (en) Method for forming distribution diagram for radiation dose of material body after correction for attenuation caused by material body thereof
JP3763165B2 (en) SPECT absorption correction method
JP2013215472A (en) Radiation tomographic image generating apparatus, radiation tomographic apparatus, and radiation tomographic image generating program
JP4371636B2 (en) Nuclear medicine diagnostic equipment
Li et al. Simulation and reconstruction for three-dimensional x-ray fluorescence computed tomography

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination