Correction parameter determination method and apparatus, image reconstruction method and apparatus, and device and medium

Abstract

A correction parameter determination method and apparatus, an image reconstruction method and apparatus, and a device and a medium. The correction parameter determination method comprises: acquiring a first relational expression between real signals corresponding to two different energy values and a second relational expression between scattering signals corresponding to the two different energy values (S110); according to the first relational expression, scanning data of a correction phantom that is obtained by means of narrow scanning and corresponds to the two different energy values, and scanning data of an air background that is obtained by means of narrow scanning and corresponds to the two different energy values, determining a first scattering correction parameter set (S120); according to the second relational expression, scanning data of the correction phantom that corresponds to different scanning modes and corresponds to the two different energy values, and scanning data of the air background that is obtained in the different scanning modes and corresponds to the two different energy values, determining a second scattering correction parameter set (S130); and using a combined set of the first scattering correction parameter set and the second scattering correction parameter set as a target scattering correction parameter set (S140).

Description

Correction parameter determination method, image reconstruction method, device, equipment and medium

[0001] This application claims priority to the Chinese patent application filed with the China Patent Office on July 4, 2023, with application number 202310807222.3, the entire contents of which are incorporated by reference into this application. Technical Field

[0002] The present application relates to the technical field of electronic computed tomography (CT), for example, to a correction parameter determination method, an image reconstruction method, an apparatus, a device and a medium. Background Art

[0003] Cone-beam spectral CT scanning is at the forefront of research in achieving quantitative cone-beam CT imaging. By combining CT scan data from an object at different energies, cone-beam spectral CT imaging can calculate the attenuation of X-rays of varying energies by different tissues and materials, thereby obtaining quantitative information. However, when performing multi-energy scans using cone-beam CT imaging, X-rays scatter within the object, resulting in scattering artifacts in the CT images.

[0004] Related technologies introduce periodic X-ray obstructions into cone-beam CT imaging systems and establish relationships between projection data and the effective energy of transmitted X-rays, as well as between the scattered signal and the effective energy. This artifact removal method requires high precision in the effective energy. Furthermore, due to the varying sizes of different body parts, the same obstruction cannot be used, resulting in poor generalization.

[0005] Summary of the Invention

[0006] The present application provides a correction parameter determination method, image reconstruction method, device, equipment and medium to solve the problem of poor generalization of cone-beam CT imaging methods in related technologies.

[0007] According to one aspect of the present application, a method for determining scatter correction parameters for cone-beam CT is provided, the method comprising: obtaining a first relationship between real signals corresponding to two different energy values, and a second relationship between scattered signals corresponding to the two different energy values; determining a first scatter correction parameter set based on the first relationship, scan data of a correction phantom acquired in a narrow scan mode corresponding to the two different energy values, and scan data of an air background acquired in the narrow scan mode corresponding to the two different energy values, wherein the scan data acquired in the narrow scan mode does not include scatter signals; determining a second scatter correction parameter set based on the second relationship, scan data of the correction phantom acquired in the narrow scan mode corresponding to the two different energy values, scan data of the correction phantom acquired in a full scan mode corresponding to the two different energy values, scan data of the air background acquired in the narrow scan mode corresponding to the two different energy values, and scan data of the air background acquired in the full scan mode corresponding to the two different energy values; and using a combination of the first scatter correction parameter set and the second scatter correction parameter set as a target scatter correction parameter set corresponding to the two different energy values.

[0008] According to another aspect of the present application, a method for cone-beam CT image reconstruction is provided, the method comprising: obtaining, based on an energy value combination corresponding to two sets of data to be imaged, a third relationship between a target scatter correction parameter set of any embodiment of the present application corresponding to the energy value combination and a total signal corresponding to the energy value combination, the third relationship being determined based on the first relationship and the second relationship; determining, based on the third relationship, the target scatter correction parameter set corresponding to the energy value combination, and the two sets of data to be imaged, a target attenuation parameter and a scatter signal in the two sets of data to be imaged; removing the scatter signal from the two sets of data to be imaged to obtain a real signal; and performing image reconstruction on the real signal based on the target attenuation parameter to obtain a target image.

[0009] According to another aspect of the present application, a device for determining scatter correction parameters for cone-beam CT is provided. The device includes: a relationship acquisition module configured to acquire a first relationship between real signals corresponding to two different energy values, and a second relationship between scattered signals corresponding to the two different energy values; a first parameter module configured to determine a first scatter correction parameter set based on the first relationship, scan data of a calibration phantom acquired in a narrow scan mode corresponding to the two different energy values, and scan data of an air background acquired in the narrow scan mode corresponding to the two different energy values, wherein the scan data acquired in the narrow scan mode does not include scatter signals; a second parameter module configured to determine a second scatter correction parameter set based on the second relationship, scan data of the calibration phantom acquired in the narrow scan mode corresponding to the two different energy values, scan data of the calibration phantom acquired in a full scan mode corresponding to the two different energy values, scan data of the air background acquired in the narrow scan mode corresponding to the two different energy values, and scan data of the air background acquired in the full scan mode corresponding to the two different energy values; and a result module configured to use a combination of the first scatter correction parameter set and the second scatter correction parameter set as a target scatter correction parameter set corresponding to the two different energy values.

[0010] According to another aspect of the present application, an image reconstruction device for cone-beam CT is provided, the device comprising: an acquisition module, configured to acquire, based on an energy value combination corresponding to two groups of data to be imaged, a target scatter correction parameter set of any embodiment of the present application corresponding to the energy value combination, and a third relationship between the total signals corresponding to the energy value combination, wherein the third relationship is determined based on the first relationship and the second relationship; a scattered signal determination module, configured to determine, based on the third relationship, the target scatter correction parameter set corresponding to the energy value combination, and the two groups of data to be imaged, a target attenuation parameter and the scattered signals in the two groups of data to be imaged; a real signal module, configured to remove the scattered signals from the two groups of data to be imaged to obtain a real signal; and an image reconstruction module, configured to reconstruct the real signal based on the target attenuation parameter to obtain a target image.

[0011] According to another aspect of the present application, an electronic device is provided, comprising: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the scatter correction parameter determination method for cone-beam CT of any embodiment of the present application in response to a correction instruction, or execute the image reconstruction method for cone-beam CT of any embodiment of the present application in response to an image reconstruction instruction.

[0012] According to another aspect of the present application, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the method for determining scatter correction parameters of cone-beam CT of any embodiment of the present application or the image reconstruction method of cone-beam CT of any embodiment of the present application when executed. BRIEF DESCRIPTION OF THE DRAWINGS

[0013] The following is a brief introduction to the drawings required for the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0014] FIG1 is a flowchart of a method for determining scatter correction parameters for cone-beam CT according to an embodiment of the present application;

[0015] FIG2 is a schematic diagram of scanning data of a calibration phantom acquired by cone-beam spectral CT in narrow field of view mode according to an embodiment of the present application;

[0016] FIG3 is a schematic diagram of a cone-beam spectral CT scan data acquisition of a calibration phantom in full-field mode according to an embodiment of the present application;

[0017] FIG4 is a flowchart of a cone-beam CT image reconstruction method provided according to an embodiment of the present application;

[0018] FIG5 is a structural block diagram of a device for determining scatter correction parameters for cone-beam CT according to an embodiment of the present application;

[0019] FIG6 is a structural block diagram of another apparatus for determining scatter correction parameters for cone-beam CT according to an embodiment of the present application;

[0020] FIG7 is a structural block diagram of another apparatus for determining scatter correction parameters for cone-beam CT according to an embodiment of the present application;

[0021] FIG8 is a structural block diagram of an image reconstruction device for cone-beam CT according to an embodiment of the present application;

[0022] FIG9 is a structural block diagram of another cone-beam CT image reconstruction device provided according to an embodiment of the present application;

[0023] FIG10 is a structural block diagram of a cone-beam CT scatter correction device provided according to an embodiment of the present application;

[0024] FIG11 is a structural block diagram of a cone-beam CT device provided according to an embodiment of the present application;

[0025] FIG12 is a structural block diagram of an electronic device provided according to an embodiment of the present application. DETAILED DESCRIPTION

[0026] The following will describe the technical solutions in the embodiments of this application in conjunction with the drawings in the embodiments of this application. Obviously, the embodiments described are only part of the embodiments of this application. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of this application.

[0027] It should be noted that the terms "first", "second" and "third" in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects and are not necessarily used to describe a specific order or sequential order. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0028] Figure 1 is a flowchart of a method for determining scatter correction parameters for cone-beam CT provided according to an embodiment of the present application. This embodiment is applicable to scenarios in which scatter correction parameters for cone-beam CT are determined based on a correction phantom. The method can be configured in a device for determining scatter correction parameters for cone-beam CT. The device for determining scatter correction parameters for cone-beam CT can be implemented in the form of hardware and / or software, and can also be configured in a processor of an electronic device.

[0029] As shown in FIG1 , a method for determining scatter correction parameters for cone-beam CT includes the following steps.

[0030] S110 , obtaining a first relational expression between real signals corresponding to two different energy values, and a second relational expression between scattered signals corresponding to two different energy values.

[0031] In one embodiment, the first relationship between the real signals corresponding to the two different energy values ​​is determined based on a first-order linear fitting equation between the equivalent attenuation coefficients corresponding to the calibration phantom at the two different energy values. The second relationship is determined based on a first-order linear expression between the scattered signals corresponding to the calibration phantom at the two different energy values.

[0032] During CT scanning, some X-ray photons are absorbed when they pass through an object, causing the intensity of the X-rays to decay exponentially. This decay is related to the density inside the object, and the degree of decay is proportional to the density. Ideally, the attenuation process of X-rays can be expressed as formula (1): I = I0e-∫μ(l)dl

[0033] (1);

[0034] I0 is the number of photons before X-ray attenuation; I is the number of photons remaining after X-ray attenuation, which is the signal ultimately received by the CT detector; μ(l) is a function of the ray path l, which represents the linear attenuation coefficient of the object's X-rays, and ∫μ(l)dl represents the total attenuation of the X-rays along that direction.

[0035] In actual CT scanning of an object, X-rays scatter inside the object, resulting in scattering artifacts in the CT image. In other words, the signal I ultimately received by the CT detector includes the object's true signal P and the scattered signal S. I=P+S (2)

[0036] According to formula (2), the present application converts different energies (E1, E2, ..., E n ) is equivalent to multiple independent data acquisition processes, so the original signals I with different energies can be expressed by the following equations:

[0037] At the first energy (E1), the original signal received by the detector is The original signal Including real signals and scattered signals and They represent the signals received by the detector under the second energy (E2) The real signal and the scattered signal in and They represent the original signal received by the detector at the nth energy The real signal and the scattered signal in .

[0038] In one embodiment, the cone-beam CT scatter correction parameter determination method of the present invention is used in dual-energy CT scanning scenarios, such as dual-source dual-detector scanning, dual-layer flat-panel detector scanning, rapid tube voltage switching scanning, and single-source double scanning. It is also applicable to cone-beam CT scanning that obtains object energy spectrum information, such as energy-resolved photon counting detectors. For scanning methods involving only two energies, simply set n in formula (3) to 2. For photon counting detector scanning and other multi-energy scanning, the value of n is determined based on the number of energy partitions configured.

[0039] In an embodiment of a dual-energy CT scan, in order to establish a first relationship between the real signals at the first energy (E1) and the second energy (E2), first, the real signals at the first energy (E1) and the second energy (E2) are determined. and

[0040] and Respectively represent the scanning data of the air background collected in the narrow scanning mode at the first energy (E1) and the second energy (E2), that is, the number of photons emitted by the light source of the air background collected in the narrow scanning mode at the first energy (E1) and the second energy (E2) or the light intensity measured by the detector; is the equivalent attenuation coefficient of the scanned object during the equivalent scanning process at the first energy (E1), is the equivalent attenuation coefficient of the scanned object during the equivalent scanning process at the second energy (E2); L is the thickness of the scanned object corresponding to the corresponding projection angle. It can be understood that the equivalent attenuation coefficient is related to the material of the scanned object.

[0041] Then, the connection between the real signals corresponding to different energy values ​​is established based on the equivalent attenuation coefficients under different energies. and The first-order linear equation of the real signal under the first energy (E1) is established and the real signal at the second energy (E2) The first-order linear function is used to represent the equivalent attenuation coefficient of the first energy (E1) and the equivalent attenuation coefficient of the second energy (E2) Relationship:

[0042] α0 and α1 are the constant term fitting coefficient and the first-order term fitting coefficient of the first-order linear equation, respectively.

[0043] Formula (4) and formula (5) are taken as the first relationship.

[0044] In order to describe the relationship between the scattered signals at different energies, a relationship between the scattered signals at two different energies is established based on the first-order linear mathematical relationship, that is, the scattered signal at the first energy (E1) is used as The first-order linear function represents the scattered signal at the second energy (E2)

[0045] β 1,2 is the constant fitting coefficient of the first-order linear equation, γ1,2 is the first-order fitting coefficient of the first-order linear equation. Formula (6) is used as the scattered signal under the first energy (E1) and the scattered signal at the second energy (E2) The second relationship between .

[0046] S120. Determine a first scatter correction parameter set based on the first relationship, scan data of the calibration phantom acquired in the narrow scan mode corresponding to two different energy values, and scan data of the air background acquired in the narrow scan mode corresponding to two different energy values, where the scan data acquired in the narrow scan mode does not include a scatter signal.

[0047] In one embodiment, the calibration phantom is a cylinder made of a single material, and the size of the calibration phantom is determined based on the size of the scanned object. For example, the cylinder can be made of acrylic or pure water to serve as the calibration phantom.

[0048] Optionally, correction phantoms of other shapes can also be used, as long as the parameters such as the shape and size of the correction phantom are the same or approximately the same as the parameters such as the shape and size of the scanned part. For example, for CT scanning of the head, the correction phantom can use a sphere made of a single material, and the diameter of the sphere is close to the size of the head to be scanned.

[0049] In one embodiment, FIG2 is a schematic diagram of a cone-beam spectral CT system acquiring scan data from a calibration phantom in narrow field of view mode, according to an embodiment of the present application. As shown in FIG2 , an X-ray source 65 emits X-rays, which pass through a filter 66, narrowing the X-ray detection field of view 69. A narrow scan is performed vertically over a calibration phantom 67. The X-rays passing through the calibration phantom 67 are attenuated, and the attenuated X-rays are detected by an energy-resolving photon counting detector 68, generating scan data for the calibration phantom 67 in narrow scan mode. The width of the beam corresponding to the narrow scan mode is less than 10 mm in a set direction of the detector. The set direction is parallel to the central axis of the cylindrical calibration phantom, which is parallel to the direction of motion of the scanning bed during scan data acquisition. It is understood that the vertical width of the narrow field of view needs to be less than a first set width. It should be noted that to prevent the number of emitted photons from being reduced by the collimator at the light source, the field of view should be no less than a second set width, which is smaller than the first set width. For example, the first set width is 10 mm, and the second set width is 5 mm. This embodiment makes the scattering signal in the scanning data corresponding to the narrow scanning mode zero by making the ray beam corresponding to the narrow scanning mode smaller than 10 mm in the set direction of the detector. That is, the scanning data acquired based on the narrow scanning mode only includes real signals.

[0050] For two different energies, air background scan data were collected in narrow scan mode. The air background scan data corresponding to each energy was used as the photon quantity before X-ray attenuation at the corresponding energy, that is, as I0 at the corresponding energy (Formula (1)).

[0051] The scan data of the calibration phantom acquired in the narrow scan mode at the first energy is used as the real signal at the first energy (E1) The scan data of the calibration phantom collected in the narrow scan mode at the second energy is used as the real signal at the second energy (E2) The above data and the air background scan data collected in the narrow scan mode under the first energy (E1) are combined Scan data of air background collected in narrow scan mode at the second energy (E2) Substituting into the first relation, we can get the equivalent attenuation coefficient of the calibration phantom at the first energy (E1): and the equivalent attenuation coefficient of the calibration phantom at the second energy (E2) Substituting into the first relationship formed by formula (4) and formula (5), α0 and α1 are obtained as the first scatter correction parameter set.

[0052] S130. Determine a second scatter correction parameter set based on the second relationship, the scan data of the calibration phantom collected in the narrow scanning mode corresponding to two different energy values, the scan data of the calibration phantom collected in the full scanning mode corresponding to two different energy values, the scan data of the air background collected in the narrow scanning mode corresponding to two different energy values, and the scan data of the air background collected in the full scanning mode corresponding to two different energy values.

[0053] Full field of view means that the X-rays emitted by the X-ray tube can completely cover the entire detector in a set direction, where the set direction is the direction of movement of the scanning bed during the scanning data acquisition process; the scanning data collected in the full field of view scanning mode includes real signals and scattered signals. In one embodiment, FIG3 is a schematic diagram of a cone-beam spectral CT according to an embodiment of the present application, which collects scanning data of a calibration phantom in full field of view mode. As shown in FIG3 , an X-ray source 65 emits X-rays, which scan the calibration phantom 67 within a full field of view 70. The X-rays passing through the calibration phantom 67 are attenuated, and the attenuated X-rays are detected by the energy-resolving photon counting detector 68 to obtain the scanning data of the calibration phantom 67 in full scanning mode.

[0054] The scan data of the calibration phantom collected in the full scan mode corresponding to the first energy is used as the original signal under the first energy (E1) The original signal Including real signals and scattered signals The scan data of the air background collected in the full scan mode corresponding to the second energy is used as the original signal under the second energy (E2) The original signal Including real signals and scattered signals The scanning data of the calibration phantom acquired in the narrow scanning mode corresponding to the first energy are subtracted from the scanning data of the calibration phantom acquired in the full scanning mode corresponding to the first energy to obtain the scattering signal of the calibration phantom corresponding to the first energy; the scanning data of the calibration phantom acquired in the narrow scanning mode corresponding to the second energy are subtracted from the scanning data of the calibration phantom acquired in the full scanning mode corresponding to the second energy to obtain the scattering signal of the calibration phantom corresponding to the second energy; the scattering signal of the calibration phantom at the first energy and the scattering signal at the second energy are combined as a first scattering signal; the scanning data of the air background acquired in the narrow scanning mode corresponding to the first energy are subtracted from the scanning data of the air background acquired in the full scanning mode corresponding to the first energy to obtain the scattering signal of the air background at the first energy; the scanning data of the air background acquired in the narrow scanning mode corresponding to the second energy are subtracted from the scanning data of the air background acquired in the full scanning mode corresponding to the second energy to obtain the scattering signal of the air background at the second energy; the scattering signal of the air background at the first energy and the scattering signal of the air background at the second energy are combined as a second scattering signal; the scattering signal at the first energy and the scattering signal at the second energy are respectively determined based on the first scattering signal combination and the second scattering signal combination; the scattering signal at the first energy and the scattering signal at the second energy are substituted into the second relational expression to obtain β 1,2 and γ 1,2 , as the second scatter correction parameter set.

[0055] The second scatter correction parameter set (β 1,2 and γ 1,2 ) and the first scatter correction parameter set (α0 and α1) are both related to the two different energy values, that is, different energy value combinations correspond to different first scatter correction parameter sets and second scatter correction parameter sets.

[0056] It should be noted that the scan data of the calibration phantom and the air background acquired in the narrow scan mode at two different energy values, as well as the scan data of the calibration phantom and the air background acquired in the full scan mode at two different energy values, correspond to the same scanning parameters, including the tube voltage and total filter thickness. For any energy value, the tube voltage and total filter thickness are related to the X-ray signal intensity. If the tube voltage and total filter thickness change, the true signal and scattered signal corresponding to the different energy values ​​will also change, resulting in changes in both the first scatter correction parameter set and the second scatter correction parameter set. Therefore, when determining the scatter correction parameters for cone-beam CT, it is necessary to maintain the tube voltage and total filter thickness constant.

[0057] S140: Use the combination of the first scatter correction parameter set and the second scatter correction parameter set as a target scatter correction parameter set corresponding to two different energy values.

[0058] In one embodiment, the target scatter correction parameter sets corresponding to the two different energy values ​​include parameters in the first scatter correction parameter set and parameters in the second scatter correction parameter set. Exemplarily, the target scatter correction parameter set includes α0 and α1 in the first scatter correction parameter set, and β in the second scatter correction parameter set. 1,2 and γ 1,2 .

[0059] About the scanning data collection of the correction phantom. Under the current energy and scanning mode, the ray source is controlled to rotate 360 ​​degrees relative to the correction phantom and the corresponding scanning data is collected during the rotation. In one embodiment, the ray source and the detector are controlled to remain stationary, and the correction phantom is controlled to rotate 360 ​​degrees. During the rotation of the correction phantom, the scanning data of the correction phantom at multiple projection angles are collected through the cooperation of the ray source and the detector. In one embodiment, the correction phantom is kept stationary, and the X-ray source and the detector are controlled to rotate 360 ​​degrees around the correction phantom, and during the rotation of the ray source and the detector, the scanning data of the correction phantom are collected at multiple projection angles. It should be noted that if the shape of the correction phantom is regular (such as a cylinder) and the material is single, the scanning data of the scanning phantom can be obtained at only one projection angle, and then the scanning data of multiple projection angles covering 360 degrees of the correction phantom are determined based on the projection data of the one projection angle.

[0060] The technical solution of the embodiment of the present application is based on the premise that the scanning data collected in the narrow scanning mode only includes real signals, and the scanning data collected in the full scanning mode includes scattered signals and real signals; based on a first relationship between the real signals corresponding to two different energy values ​​and the scanning data of the correction phantom corresponding to the two different energy values ​​in the narrow scanning mode and the scanning data of the air background, a first scattering correction parameter set is determined; based on a second relationship between the scattering signals corresponding to the two different energy values, the scanning data of the correction phantom corresponding to the two different energy values ​​in the narrow scanning mode and the scanning data of the air background, and the scanning data of the correction phantom corresponding to the two different energy values ​​in the full scanning mode and the scanning data of the air background, a second scattering correction parameter set is determined; and the combination of the first scattering correction parameter set and the second scattering correction parameter set is used as the target scattering correction parameter set to improve the accuracy and generalization of the scattering correction parameter determination method.

[0061] Figure 4 is a flowchart of a cone-beam CT image reconstruction method provided according to an embodiment of the present application. This embodiment is applicable to scenarios of cone-beam CT image reconstruction based on scatter correction. The method can be configured in a cone-beam CT image reconstruction method device. The cone-beam CT image reconstruction method device can be implemented in the form of hardware and / or software, and can also be configured in a processor of an electronic device.

[0062] As shown in FIG4 , the method includes the following steps.

[0063] S210 , according to the energy value combination corresponding to the two sets of to-be-imaged data, obtaining a third relationship between the target scatter correction parameter set corresponding to the energy value combination and the total signal corresponding to the energy value combination, wherein the third relationship is determined based on the first relationship and the second relationship.

[0064] A first energy (E1) and a second energy (E2) respectively corresponding to two groups of data to be imaged are determined, and a combination of the first energy and the second energy is used as an energy value combination.

[0065] Substituting formula (5) into the second equation of formula (4), we get:

[0066] In formula (7), Denoted as variable χ, formula (4) can be:

[0067] Substituting formula (8) into formula (3), we obtain the third relationship between the target scatter correction parameter set corresponding to the energy value combination and the total signal corresponding to the energy value combination:

[0068] Acquiring a target scatter correction parameter set of any embodiment of the present application corresponding to the energy value combination according to the energy value combination corresponding to two groups of to-be-imaged data includes the following steps.

[0069] Step a1: Determine the energy value combination and scanning parameters corresponding to the two sets of data to be imaged.

[0070] The two sets of data to be scanned are analyzed to determine the energy values ​​and scanning parameters corresponding to the two sets of imaging data. The combination of the energy values ​​corresponding to the two sets of imaging scan data is used as the energy value combination. The scanning parameters include parameters such as the scanning mode, tube voltage, and total filter thickness. The scanning mode includes a target parameter, which is the width of the cone beam in the vertical direction of the detector, or other parameters used to determine the width.

[0071] Step a2: Obtain a target scatter correction parameter set of any embodiment of the present application corresponding to the energy value combination and the scanning parameter.

[0072] A correspondence between an energy value combination, a scanning parameter, and a target scatter correction parameter set is pre-created. After the energy value combination and the scanning parameter are determined, the target scatter correction parameter set corresponding to the energy value combination and the scanning parameter is determined based on the correspondence.

[0073] S220 , determining the target attenuation parameter and the scattering signals in the two sets of data to be imaged according to the third relationship, the target scattering correction parameter set corresponding to the energy value combination, and the two sets of data to be imaged.

[0074] Since χ is related to the attenuation coefficient of the scanned object (the material), the value of this parameter is different for different scanned objects. Therefore, χ is a variable here, and the third relationship shown in formula (9) contains χ and Two variables. According to the third relation, the two sets of imaging data ( and ), the four parameters α0, α1, β in the target scattering correction parameter set corresponding to the energy value combination 1,2 and γ 1,2 Substituting into the third relation and iteratively solving the equations in the third relation, we can obtain the target attenuation parameter (χ) and the scattering signals in the two sets of imaging data. Further obtain the scattered signal of the object at the second energy (E2)

[0075] S230 , removing scattered signals from the two sets of data to be imaged to obtain real signals.

[0076] The scattered signal corresponding to the first energy Scattered signal corresponding to the second energy After obtaining, the imaged data corresponding to the first energy are respectively Remove scattered signals To get the corresponding real signal From the data to be imaged corresponding to the second energy Remove the corresponding scattered signal To get the corresponding real signal

[0077] In one embodiment, high-frequency structural information in the scattered signals corresponding to the two energy values ​​in the energy value combination is removed to update the scattered signals corresponding to the two energy values; the updated corresponding scattered signals are removed from the two sets of imaging data to obtain the real signal.

[0078] It can be understood that the scattered signal includes high-frequency structural information. The obtained scattered signal can be filtered to remove the high-frequency structural information in the scattered signal to obtain an updated scattered signal; the updated scattered signal is removed from the two sets of imaging data to obtain the real signal.

[0079] In one embodiment, an iterative fuzzy method is used to remove high-frequency structural information from the scattered signals corresponding to the two energy values ​​in the energy value combination. For example, first, according to formula (10), the high-frequency structural information of the scattered signal is smoothed to obtain a smoothed scattered signal, which is the updated scattered signal.

[0080] S is the scattered signal, i and j are used to represent the position of the detector pixel (for example, the horizontal and vertical coordinates), k is the number of iterations, and ω and γ are weight factors, and their values ​​can be set according to actual conditions.

[0081] After the scattering signal is updated, the updated scattering signal is subtracted from the two sets of data to be imaged to obtain the real signal.

[0082] S240 , performing image reconstruction on the real signal based on the target attenuation parameter to obtain a target image.

[0083] After the target attenuation parameters and the real signal corresponding to the two sets of imaging data are determined, the real signal is reconstructed based on the target attenuation parameters to obtain the target image.

[0084] The technical solution of this embodiment removes high-frequency structural information from the scattering signals corresponding to the two energy values ​​in the energy value combination through smoothing filtering to obtain an updated scattering signal; removes the updated scattering signal from the two sets of imaging data to obtain a real signal; and reconstructs the real signal based on the target attenuation parameter to obtain the target image, thereby solving the problem of scattering artifacts in cone-beam CT images. Based on scatter correction, a formal signal is obtained for image reconstruction, thereby improving the accuracy and generalization of the method for determining the scatter correction parameters, and further improving the accuracy of the cone-beam CT image reconstruction method.

[0085] Figure 5 is a block diagram of a device for determining scatter correction parameters for cone-beam CT according to an embodiment of the present application. The device can be implemented in the form of hardware and / or software and integrated into a processor of an electronic device.

[0086] As shown in FIG5 , the apparatus for determining scatter correction parameters for cone-beam CT includes: a relationship acquisition module 501 configured to acquire a first relationship between real signals corresponding to two different energy values, and a second relationship between scattered signals corresponding to the two different energy values; a first parameter module 502 configured to determine a first scatter correction parameter set based on the first relationship, scan data of a calibration phantom acquired in a narrow scan mode corresponding to the two different energy values, and scan data of an air background acquired in the narrow scan mode corresponding to the two different energy values, wherein the scan data acquired in the narrow scan mode does not include scatter signals; a second parameter module 503 configured to determine a second scatter correction parameter set based on the second relationship, scan data of the calibration phantom acquired in the narrow scan mode corresponding to the two different energy values, scan data of the calibration phantom acquired in a full scan mode corresponding to the two different energy values, scan data of the air background acquired in the narrow scan mode corresponding to the two different energy values, and scan data of the air background acquired in the full scan mode corresponding to the two different energy values; and a result module 504 configured to use the combination of the first scatter correction parameter set and the second scatter correction parameter set as a target scatter correction parameter set corresponding to the two different energy values.

[0087] Optionally, as shown in Figure 6, the device also includes a first relationship determination module 505, which is configured to: determine the first relationship between the real signals corresponding to the two different energy values ​​based on the first-order linear fitting equation between the equivalent attenuation coefficients corresponding to the correction model at two different energy values.

[0088] Optionally, as shown in Figure 7, the device also includes a second relationship determination module 506, which is configured to: determine the second relationship between the scattering signals corresponding to two different energy values ​​based on the first-order linear expression between the scattering signals corresponding to the correction model at two different energy values.

[0089] The technical solution of this embodiment, through the mutual cooperation of multiple modules, uses the scan data collected in the narrow scan mode as the real signal, and the scan data collected in the full scan mode includes the scattered signal and the real signal. Based on a first relationship between the real signals corresponding to two different energy values, a first scatter correction parameter set is determined, and based on a second relationship between the scatter signals corresponding to the two different energy values, a second scatter correction parameter set is determined. In this way, the target scatter correction parameter sets corresponding to the two different energy values ​​are obtained, thereby improving the accuracy and generalization of the scatter correction parameter determination method.

[0090] The apparatus for determining scatter correction parameters for cone-beam CT provided in the embodiments of the present application can execute the method for determining scatter correction parameters for cone-beam CT provided in any embodiment of the present application, and has functional modules corresponding to the execution method.

[0091] FIG8 is a structural block diagram of an image reconstruction device for cone-beam CT according to an embodiment of the present application. The device may be implemented in the form of hardware and / or software and integrated into a processor of an electronic device.

[0092] As shown in FIG8 , the image reconstruction device of the cone-beam CT includes: an acquisition module 601, configured to acquire, based on an energy value combination corresponding to two sets of data to be imaged, a target scatter correction parameter set of any embodiment of the present application corresponding to the energy value combination, and a third relationship between the total signals corresponding to the energy value combination, wherein the third relationship is determined based on the first relationship and the second relationship; a scattered signal determination module 602, configured to determine, based on the third relationship, the target scatter correction parameter set corresponding to the energy value combination, and the two sets of data to be imaged, a target attenuation parameter and the scattered signals in the two sets of data to be imaged; a real signal module 603, configured to remove the scattered signals from the two sets of data to be imaged to obtain a real signal; and an image reconstruction module 604, configured to reconstruct the real signal based on the target attenuation parameter to obtain a target image.

[0093] Optionally, as shown in Figure 9, the device also includes a correction parameter set module 605, which is configured to: determine the energy value combination and scanning parameters corresponding to the two sets of data to be imaged; and obtain the target scattering correction parameter set of any embodiment of the present application corresponding to the energy value combination and scanning parameters.

[0094] Optionally, the real signal module 603 is configured to: remove high-frequency structural information in the scattered signals corresponding to the two energy values ​​in the energy value combination to update the scattered signals corresponding to the two energy values; and remove the updated scattered signals corresponding to the two energy values ​​from the two sets of data to be imaged to obtain the real signal.

[0095] The technical solution of this embodiment, through the interaction of multiple modules, obtains a third relationship between the target scatter correction parameter set of any embodiment of the present application corresponding to the energy value combination and the total signal corresponding to the energy value combination based on the energy value combination corresponding to the two sets of data to be imaged, wherein the third relationship is determined based on the first relationship and the second relationship; determines the target attenuation parameter and the scatter signals in the two sets of data to be imaged based on the third relationship, the target scatter correction parameter set corresponding to the energy value combination, and the two sets of data to be imaged; removes the scatter signals from the two sets of data to be imaged to obtain the true signal; and reconstructs the true signal based on the target attenuation parameter to obtain the target image. This solves the problem of scatter artifacts in cone-beam CT images and improves the accuracy of cone-beam CT image reconstruction.

[0096] The cone-beam CT image reconstruction device provided in the embodiment of the present application can execute the cone-beam CT image reconstruction device method provided in any embodiment of the present application, and has functional modules corresponding to the execution method.

[0097] FIG10 is a block diagram of a cone-beam CT scatter correction device according to an embodiment of the present application. As shown in FIG10 , the cone-beam CT scatter correction device 800 includes a cone-beam CT device 801 configured to collect rays of two set energy values ​​to scan a scanned object to obtain scan data corresponding to the two set energy values; a scatter correction device 802 configured to acquire scan data from the cone-beam CT device 801 to obtain scan data of a correction phantom collected in a narrow scan mode corresponding to the two set energy values, scan data of an air background collected in the narrow scan mode corresponding to two different energy values, scan data of the correction phantom collected in a full scan mode corresponding to two different energy values, and scan data of an air background collected in a full scan mode corresponding to two different energy values; at least one processor 803; and a memory communicatively connected to the at least one processor; wherein the memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for determining scatter correction parameters for cone-beam CT according to any embodiment of the present application in response to an image correction instruction.

[0098] At a first energy value, the cone-beam CT device 801 uses a narrow scanning mode and a full scanning mode to scan the correction phantom, respectively, to obtain scanning data of the correction phantom acquired in the narrow scanning mode corresponding to the first energy value and scanning data of the correction phantom acquired in the full scanning mode; at the first energy value, the narrow scanning mode and the full scanning mode are used to scan the air background, respectively, to obtain scanning data of the air background in the narrow scanning mode corresponding to the first energy value and scanning data of the air background acquired in the full scanning mode; at a second energy value, the full scanning mode and the narrow scanning mode of the cone-beam CT device 801 are used to scan the correction phantom, respectively, to obtain scanning data of the correction phantom acquired in the full scanning mode corresponding to the second energy value and scanning data of the correction phantom acquired in the narrow scanning mode corresponding to the second energy value; at the second energy value, the full scanning mode and the narrow scanning mode of the cone-beam CT device 801 are used to scan the correction phantom, respectively. Scan the air background to obtain scanning data of the air background in the full scanning mode corresponding to the second energy value and scanning data of the air background in the narrow scanning mode corresponding to the second energy value; the cone-beam CT device 801 sends the above scanning data to the processor 803, the processor 803 receives the above scanning data, and sends the above scanning data to the memory 804, and the memory 804 stores the above scanning data; the processor 803 obtains a first relationship between the real signals of the first energy value and the second energy value, and a second relationship between the scattered signals of the first energy value and the second energy value, and obtains a first scatter correction parameter set and a second scatter correction parameter set based on the first relationship, the second relationship and the above scanning data; the processor 803 combines the above first scatter correction parameter set and the second scatter correction parameter set to obtain a target scatter correction parameter set corresponding to the first energy and the second energy.

[0099] The technical solution of the cone-beam CT scatter correction device provided in this embodiment uses the scan data collected in the narrow scan mode as the real signal, and the scan data collected in the full scan mode includes the scattered signal and the real signal. Based on a first relationship between the real signals corresponding to two different energy values, a first scatter correction parameter set is determined; based on a second relationship between the scatter signals corresponding to the two different energy values, a second scatter correction parameter set is determined, thereby obtaining a target scatter correction parameter set corresponding to the two different energy values, thereby improving the accuracy and generalization of the scatter correction parameter determination method.

[0100] The cone-beam CT scatter correction device provided in the embodiment of the present application can implement the cone-beam CT scatter correction parameter determination method provided in any embodiment of the present application, and has a functional module corresponding to the execution method.

[0101] Figure 11 is a structural block diagram of a cone-beam CT device provided according to an embodiment of the present application. As shown in Figure 11, the cone-beam CT device 700 includes a scanning table 701, which is configured to support a scanned object; a scanning device 702, which is configured to scan the scanned object using rays of two set energy values ​​to determine two sets of imaging data corresponding to the scanned object; at least one processor 703; and a memory 704 communicatively connected to the at least one processor 703; wherein the memory 704 stores a computer program executable by the at least one processor 703, and the computer program is executed by the at least one processor 703 so that the at least one processor 703 can execute the image reconstruction method of the cone-beam CT of any embodiment of the present application in response to an image execution instruction.

[0102] The scanning device 702 uses two different energy values ​​to scan the object on the scanning platform 701. Two sets of data to be imaged are obtained. The scanning device 702 sends the two sets of data to be imaged to the processor 703. The processor 703 receives the two sets of data to be imaged and sends the two sets of data to be imaged to the memory 704. The memory 704 stores the two sets of data to be imaged. The memory 704 stores multiple sets of energy value combinations of two set energy values, corresponding sets of target scatter correction parameters, and a third relationship. Based on the energy value combinations of two different energy values ​​of the two sets of data to be imaged, the processor 703 obtains a third relationship between the set of target scatter correction parameters corresponding to the energy value combination and the total signal corresponding to the energy value combination from the memory 704. The processor 703 determines the target attenuation parameter and the scattering signal in the data to be imaged based on the third relationship, the set of target scatter correction parameters corresponding to the energy value combination, and the two sets of data to be imaged. The processor 703 then removes high-frequency structure information from the scattering signals corresponding to the two energy values ​​in the energy value combination to update the scattering signals corresponding to the two energy values. The scattering signals corresponding to the two updated energy values ​​are removed from the two sets of data to be imaged to obtain a true signal. The true signal is reconstructed based on the target attenuation parameter to obtain a target image.

[0103] The technical solution of the cone-beam CT device provided in this embodiment removes high-frequency structural information from the scattering signals corresponding to the two energy values ​​in the energy value combination through smoothing filtering to obtain an updated scattering signal; removes the updated scattering signal from the data to be imaged to obtain a real signal; and reconstructs the real signal based on the target attenuation parameter to obtain the target image, thereby solving the problem of scattering artifacts in cone-beam CT images, improving the accuracy and generalization of the method for determining the scattering correction parameters, and improving the accuracy of the cone-beam CT image reconstruction method.

[0104] The cone-beam CT device provided in the embodiments of the present application can implement the cone-beam CT image reconstruction method provided in any embodiment of the present application, and has functional modules corresponding to the execution method.

[0105] Figure 12 is a block diagram of an electronic device provided according to an embodiment of the present application. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present application described and / or required herein.

[0106] As shown in FIG12 , the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12 and a random access memory (RAM) 13, that is communicatively connected to the at least one processor 11. The memory stores a computer program that can be executed by the at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the ROM 12 or the computer program loaded from the storage unit 18 into the RAM 13. The RAM 13 can also store various programs and data required for the operation of the electronic device 10. The processor 11, ROM 12, and RAM 13 are connected to each other via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0107] Multiple components in the electronic device 10 are connected to the I / O interface 15, including an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a magnetic disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks.

[0108] The processor 11 can be a variety of general-purpose and / or specialized processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various specialized artificial intelligence (AI) computing chips, various processors that run machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the method for determining scatter correction parameters for cone-beam CT and / or the method for image reconstruction for cone-beam CT.

[0109] In some embodiments, the method for determining scatter correction parameters for cone-beam CT and the method for reconstructing images for cone-beam CT can be implemented as a computer program tangibly embodied in a computer-readable storage medium, such as the storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the method for determining scatter correction parameters for cone-beam CT and the method for reconstructing images for cone-beam CT described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to execute the method for determining scatter correction parameters for cone-beam CT and / or the method for reconstructing images for cone-beam CT by any other appropriate means (e.g., via firmware).

[0110] Various embodiments of the systems and techniques described herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard parts (ASSPs), system on chips (SOCs), complex programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include being implemented in one or more computer programs that are executable and / or interpreted on a programmable system that includes at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0111] Computer programs for implementing the methods of the present application may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when the computer program is executed by the processor, the functions / operations specified in the flowcharts and / or block diagrams are implemented. The computer program may be executed entirely on the machine, partially on the machine, as a stand-alone software package, partially on the machine and partially on a remote machine, or entirely on a remote machine or server.

[0112] In the context of the present application, a computer-readable storage medium can be a tangible medium that can contain or store a computer program for use by an instruction execution system, device or equipment or used in combination with an instruction execution system, device or equipment. A computer-readable storage medium can include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium can be a machine-readable signal medium. Examples of machine-readable storage media can include electrical connections based on one or more lines, portable computer disks, hard disks, RAM, ROM, erasable programmable read-only memory (EPROM) or flash memory, optical fibers, portable compact disc read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.

[0113] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a cathode ray tube (CRT) or a liquid crystal display (LCD) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0114] The systems and techniques described herein can be implemented in a computing system that includes back-end components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes front-end components (e.g., a user computer with a graphical user interface or web browser through which a user can interact with embodiments of the systems and techniques described herein), or a computing system that includes any combination of such back-end components, middleware components, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0115] A computing system may include a client and a server. The client and server are generally remote from each other and typically interact via a communication network. The client-server relationship arises through computer programs running on the respective computers and establishing a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, a host product within a cloud computing service ecosystem that addresses the management difficulties and limited scalability of traditional physical hosts and virtual private server (VPS) services.

[0116] It should be understood that the various forms of the processes shown above can be used to reorder, add, or delete steps. For example, the multiple steps described in this application can be performed in parallel, sequentially, or in a different order, as long as the desired results of the technical solution of this application can be achieved. This is not limited herein.

Claims

1. A method for determining scatter correction parameters for cone-beam computed tomography (CT), comprising: Acquire a first relational expression between real signals corresponding to two different energy values, and a second relational expression between scattered signals corresponding to the two different energy values; determining a first scatter correction parameter set according to the first relationship, scan data of the calibration phantom collected in the narrow scanning mode corresponding to the two different energy values, and scan data of the air background collected in the narrow scanning mode corresponding to the two different energy values, wherein the scan data collected in the narrow scanning mode does not include a scatter signal; determining a second scatter correction parameter set according to the second relationship, the scan data of the calibration phantom collected in the narrow scanning mode corresponding to the two different energy values, the scan data of the calibration phantom collected in the full scanning mode corresponding to the two different energy values, the scan data of the air background collected in the narrow scanning mode corresponding to the two different energy values, and the scan data of the air background collected in the full scanning mode corresponding to the two different energy values; A set of the first scatter correction parameter set and the second scatter correction parameter set is used as a target scatter correction parameter set corresponding to the two different energy values.

2. The method according to claim 1, wherein: Determining the first relationship includes: Based on the first-order linear fitting equation between the equivalent attenuation coefficients respectively corresponding to the calibration phantom at the two different energy values, a first relationship between the real signals corresponding to the two different energy values ​​is determined.

3. The method according to claim 1, wherein: Determining the second relationship includes: Based on a first-order linear expression between the scattered signals corresponding to the calibration phantom at the two different energy values, a second relationship between the scattered signals corresponding to the two different energy values ​​is determined.

4. The method according to claim 1, wherein: The calibration phantom is a cylinder made of a single material, and the size of the calibration phantom is determined based on the size of the scanned object.

5. The method according to claim 1, wherein: The scanning data of the calibration model and the scanning data of the air background collected in the narrow scanning mode corresponding to the two different energy values, and the scanning data of the calibration model and the scanning data of the air background collected in the full scanning mode corresponding to the two different energy values ​​correspond to the same scanning parameters, and the scanning parameters include tube voltage and total filter thickness.

6. The method according to claim 1, wherein: The width of the ray beam corresponding to the narrow scanning mode in the set direction of the detector is less than 10 mm, and the set direction is a direction parallel to the central axis of the cylindrical correction phantom, and the central axis of the cylindrical correction phantom is parallel to the movement direction of the scanning bed during the scanning data acquisition process.

7. A method for image reconstruction of cone-beam computed tomography (CT), comprising: According to the energy value combination corresponding to the two groups of data to be imaged, obtaining a third relationship between the target scatter correction parameter set according to any one of claims 1 to 6 corresponding to the energy value combination and the total signal corresponding to the energy value combination, wherein the third relationship is determined based on the first relationship and the second relationship; Determine a target attenuation parameter and scattering signals in the two sets of data to be imaged according to the third relationship, the set of target scatter correction parameters corresponding to the energy value combination, and the two sets of data to be imaged; removing the scattered signal from the two sets of data to be imaged to obtain a real signal; The real signal is reconstructed based on the target attenuation parameter to obtain a target image.

8. The method according to claim 7, wherein: The step of obtaining, according to the energy value combination corresponding to the two groups of data to be imaged, a target scatter correction parameter set according to any one of claims 1 to 6 corresponding to the energy value combination, comprises: Determining energy value combinations and scanning parameters corresponding to the two groups of data to be imaged; Acquire the target scatter correction parameter set according to any one of claims 1 to 6 corresponding to the energy value combination and the scanning parameter.

9. The method according to claim 7, wherein: The step of removing the scattered signal from the two sets of data to be imaged to obtain a real signal comprises: Removing high-frequency structural information in the scattered signals corresponding to two energy values ​​in the energy value combination to update the scattered signals corresponding to the two energy values; The scattering signals corresponding to the two updated energy values ​​are respectively removed from the two groups of data to be imaged to obtain real signals.

10. A device for determining scatter correction parameters for cone-beam computed tomography (CT), comprising: A relational expression acquisition module, configured to acquire a first relational expression between real signals corresponding to two different energy values, and a second relational expression between scattered signals corresponding to two different energy values; a first parameter module configured to determine a first scatter correction parameter set according to the first relationship, scan data of the calibration phantom collected in the narrow scan mode corresponding to the two different energy values, and scan data of the air background collected in the narrow scan mode corresponding to the two different energy values, wherein the scan data collected in the narrow scan mode does not include a scatter signal; a second parameter module, configured to determine a second scatter correction parameter set according to the second relationship, the scan data of the calibration phantom collected in the narrow scanning mode corresponding to the two different energy values, the scan data of the calibration phantom collected in the full scanning mode corresponding to the two different energy values, the scan data of the air background collected in the narrow scanning mode corresponding to the two different energy values, and the scan data of the air background collected in the full scanning mode corresponding to the two different energy values; The result module is configured to use a set of the first scatter correction parameter set and the second scatter correction parameter set as a target scatter correction parameter set corresponding to the two different energy values.

11. An image reconstruction device for cone-beam computed tomography (CT), comprising: an acquisition module, configured to acquire, according to an energy value combination corresponding to two groups of data to be imaged, a target scatter correction parameter set according to any one of claims 1 to 6 corresponding to the energy value combination, and a third relationship between a total signal corresponding to the energy value combination, wherein the third relationship is determined based on the first relationship and the second relationship; a scattered signal determination module, configured to determine a target attenuation parameter and scattered signals in the two sets of data to be imaged according to the third relationship, the correction parameter set corresponding to the energy value combination, and the two sets of data to be imaged; A real signal module, configured to remove the scattered signal from the two sets of data to be imaged to obtain a real signal; The image reconstruction module is configured to perform image reconstruction on the real signal based on the target attenuation parameter to obtain a target image.

12. An electronic device comprising: at least one processor; as well as a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the scatter correction parameter determination method of the cone-beam electronic computer tomography (CT) described in any one of claims 1-6 in response to a correction instruction, or execute the image reconstruction method of the cone-beam electronic computer tomography (CT) described in any one of claims 7-9 in response to an image reconstruction instruction.

13. A computer-readable storage medium storing computer instructions, wherein the computer instructions are used to enable a processor to implement the scatter correction parameter determination of the cone-beam electronic computer tomography (CT) according to any one of claims 1 to 6 or the image reconstruction method of the cone-beam electronic computer tomography (CT) according to any one of claims 7 to 9 when executed.

Images