Cross scatter correction method and device, computer equipment and readable storage medium

Abstract

The invention relates to a cross scatter correction method and device, computer equipment, a computer readable storage medium and a computer program product, relates to the technical field of medical images, and can improve the accuracy of cross scatter correction. The method is applied to a double-source or multi-source CT imaging system, the double-source or multi-source CT imaging system comprises at least two detectors, each detector comprises a plurality of channels, and the method comprises the following steps: acquiring a signal intensity hopping relationship corresponding to the double-source or multi-source CT imaging system; the signal intensity hopping relationship is determined according to the detection signal intensity of each channel of the plurality of detectors, and the detection signal intensity is obtained by scanning different motifs and / or the same motif at different positions by the detectors; obtaining scattering estimation information obtained by scanning the to-be-measured object by the dual-source or multi-source CT imaging system; and correcting the scattering estimation information according to the signal intensity hopping relationship.

Description

Technical Field

[0001] The present application relates to the field of medical imaging technology, and in particular to a cross-scatter correction method, apparatus, computer equipment, and computer-readable storage medium. Background Art

[0002] With the advancement of medical technology, computed tomography (CT) images can now be acquired using dual-source CT imaging systems. A dual-source CT imaging system consists of two sets of X-ray tubes and detectors. The two X-ray tubes can scan an object from different angles, and the two detectors can each acquire corresponding scan data. During the imaging process, cross-scattering is a physical factor that seriously affects image quality. Cross-scattering occurs when photons are scattered by an object and then transmitted to a side detector.

[0003] To suppress cross-scatter, relevant scatter suppression hardware, such as an anti-scatter grid (ASG), can be installed. However, the inventors have found in practice that even with scatter suppression hardware installed, the reconstructed image obtained after cross-scatter correction still has ring or banding artifacts, and the accuracy of cross-scatter correction still needs to be improved. Summary of the Invention

[0004] Based on this, it is necessary to provide a cross-scatter correction method, apparatus, computer equipment, computer-readable storage medium and computer program product to address the above technical problems.

[0005] In a first aspect, the present application provides a cross-scatter correction method, which is applied to a dual-source or multi-source CT imaging system, wherein the dual-source or multi-source CT imaging system includes at least two detectors, each of which includes multiple channels. The method includes:

[0006] Obtaining a signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system; the signal intensity jump relationship is determined based on detection signal intensities of respective channels of the plurality of detectors, the detection signal intensities being obtained by the detectors scanning different phantoms and / or the same phantom at different positions;

[0007] Acquiring scatter estimation information obtained by scanning the object to be measured by the dual-source or multi-source CT imaging system;

[0008] The scattering estimation information is corrected according to the signal strength jump relationship.

[0009] In one embodiment, the detection signal strength includes cross-scattering signal strength, and obtaining the signal strength jump relationship corresponding to the dual-source or multi-source CT imaging system includes:

[0010] Determine the jump amplitude and the target signal intensity before the jump when the signal intensity of the channel undergoes a jump based on the multiple cross-scattering signal intensities;

[0011] Perform function fitting based on the jump amplitude and the target signal intensity corresponding to each signal intensity jump, and obtain the signal intensity jump relationship according to the fitting result.

[0012] In one embodiment, determining the jump amplitude when the signal intensity of the channel undergoes a jump based on the multiple cross-scattering signal intensities includes:

[0013] Obtain cross-scattering signal curves corresponding to the multiple cross-scattering signal intensities;

[0014] Perform curve smoothing on the cross-scattering signal curves to obtain smoothed cross-scattering signal curves;

[0015] Determine the jump amplitude when the signal intensity of the channel undergoes a jump based on the difference between the cross-scattering signal curves and the smoothed cross-scattering signal curves.

[0016] In one embodiment, the determining the jump amplitude when the signal intensity of the channel undergoes a jump based on the difference between the cross-scattering signal curves and the smoothed cross-scattering signal curves includes:

[0017] Subtract the cross-scattering signal curves from the smoothed cross-scattering signal curves to obtain signal intensity differences corresponding to each channel;

[0018] For the signal intensity difference corresponding to each channel, when the signal intensity difference satisfies the signal intensity jump condition, determine the signal intensity difference as the jump amplitude corresponding to the channel.

[0019] In one embodiment, the obtaining the signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system includes:

[0020] Obtain a phantom reconstruction image based on the multiple detection signal intensities;

[0021] Perform pixel value smoothing on the abnormal regions with abnormal pixel values in the phantom reconstruction image to obtain a smoothed reconstruction image;

[0022] Determine the jump amplitude and the target signal intensity before the jump when the signal intensity of the channel undergoes a jump based on the difference between the phantom reconstruction image and the smoothed reconstruction image;

[0023] Function fitting is performed according to the jump amplitude corresponding to each signal strength jump and the target signal strength, and the signal strength jump relationship is obtained according to the fitting result.

[0024] In one embodiment, each pixel in the phantom reconstructed image corresponds to a corresponding channel;

[0025] The determining, based on the difference between the phantom reconstructed image and the smoothed reconstructed image, a jump amplitude when the signal strength of the channel jumps, includes:

[0026] Subtracting the phantom reconstructed image from the smoothed reconstructed image to obtain a pixel value difference corresponding to each pixel point in the phantom reconstructed image;

[0027] For the pixel value difference corresponding to each pixel point, when the pixel value difference meets the pixel value jump condition, the jump amplitude when the signal strength jump occurs in the channel corresponding to the pixel point is determined according to the pixel value difference.

[0028] In one embodiment, the correcting the scattering estimation information according to the signal strength jump relationship includes:

[0029] determining the cross-scatter signal strength of each of the channels according to the scatter estimation information;

[0030] Determining a jump amplitude of the cross-scattering signal intensity of the channel according to the signal intensity jump relationship;

[0031] The cross-scattering signal strengths in the scattering estimation information are adjusted according to the jump amplitudes to obtain corrected scattering estimation information.

[0032] In a second aspect, the present application further provides a cross-scatter correction device, which is applied to a dual-source or multi-source CT imaging system, wherein the dual-source or multi-source CT imaging system includes at least two detectors, each of which includes multiple channels, and the device includes:

[0033] a relationship acquisition module, configured to acquire a signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system; the signal intensity jump relationship is determined based on the detection signal strength of each channel of the plurality of detectors, the detection signal strength being obtained by scanning different phantoms and / or the same phantom at different positions by the detectors;

[0034] a scatter estimation acquisition module, configured to acquire scatter estimation information obtained by scanning the object to be measured by the dual-source or multi-source CT imaging system;

[0035] A calibration module for calibrating the scattering estimation information according to the signal strength jump relationship.

[0036] In a third aspect, the present application further provides a computer device, including a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the steps of the cross-scattering calibration method described in any one of the above are implemented.

[0037] In a fourth aspect, the present application further provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the cross-scattering calibration method described in any one of the above are implemented.

[0038] In a fifth aspect, the present application further provides a computer program product, including a computer program. When the computer program is executed by a processor, the steps of the cross-scattering calibration method described in any one of the above are implemented.

[0039] The above cross-scattering calibration method, device, computer device, computer-readable storage medium, and computer program product can be applied to a dual-source or multi-source CT imaging system. The dual-source or multi-source CT imaging system includes at least two detectors, and each detector includes multiple channels. In this method, on the one hand, the signal strength jump relationship corresponding to the dual-source or multi-source CT imaging system can be obtained. The signal strength jump relationship is determined according to the detection signal strengths of each channel of multiple detectors. The detection signal strength is obtained by scanning different phantoms and / or the same phantom at different positions by the detectors. On the other hand, the scattering estimation information obtained by scanning the object to be measured by the dual-source or multi-source CT imaging system can be obtained, and then the scattering estimation information can be calibrated according to the signal strength jump relationship. In this embodiment, by pre-obtaining the signal strength jump relationship, the signal strength jump situation of the cross-scattering signal strength detected by the detector channels under different signal strengths can be established, which can accurately reflect the influence on the detector channels under different scattering suppression efficiencies. Therefore, when the scattering estimation information of the object to be measured is obtained, by using the signal strength jump relationship to correct the scattering estimation information, the correction result can more accurately reflect the actual situation of the detector channels receiving the cross-scattering signal considering the difference in scattering suppression efficiency, and the abnormal high-frequency signals superimposed due to the difference in scattering suppression efficiency can be removed, improving the accuracy of the cross-scattering calibration effect. Description of the Drawings

[0040] To more clearly illustrate the technical solutions in the embodiments of the present application or the related art, the following will briefly introduce the accompanying drawings required for the description of the embodiments of the present application or the related art. Obviously, the accompanying drawings in the following description are only some embodiments of the present application. For those of ordinary skill in the art, without creative efforts, other related accompanying drawings can also be obtained based on these drawings.

[0041] Figure 1 It is a schematic flowchart of a cross-scattering correction method in an embodiment;

[0042] Figure 2 It is a phantom reconstruction image in an embodiment;

[0043] Figure 3 It is a schematic diagram of a group of cross-scattering signal curves in an embodiment;

[0044] Figure 4 It is a schematic flowchart of the steps of determining the signal intensity jump relationship in an embodiment;

[0045] Figure 5 It is a schematic diagram of a signal intensity jump relationship in an embodiment;

[0046] Figure 6 It is a schematic flowchart of another cross-scattering correction method in an embodiment;

[0047] Figure 7 It is a structural block diagram of a cross-scattering correction device in an embodiment;

[0048] Figure 8 It is an internal structure diagram of a computer device in an embodiment;

[0049] Figure 9 It is another internal structure diagram of a computer device in an embodiment. Detailed implementation manners

[0050] In order to make the objectives, technical solutions, and advantages of the present application clearer, the following further details the present application in combination with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0051] In one embodiment, as Figure 1As shown, a cross-scattering correction method is provided. In this embodiment, the method is described as being applied to a dual-source or multi-source CT imaging system. The dual-source or multi-source CT imaging system includes at least two detectors, and each detector may include multiple channels. It can be understood that this method can also be applied to a server, or to a system including a dual-source or multi-source CT imaging system and a server, and is implemented through the interaction between the CT imaging system and the server. In this embodiment, the method includes the following steps:

[0052] S101, obtain the signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system; the signal intensity jump relationship is determined according to the detection signal intensities of each channel of multiple detectors, and the detection signal intensities are obtained by scanning different phantoms and / or the same phantom at different positions by the detectors.

[0053] In practical applications, in order to suppress or reduce the cross-scattering of the dual-source CT imaging system, a widely used processing method is to install relevant scatter suppression hardware in the CT imaging system, such as an anti-scatter grid (ASG). However, it has been found in practice that due to the non-ideality of the relative pose between the scatter suppression hardware and the detectors during installation, such as slight changes in the installation angle or installation position, the scatter suppression efficiency of the scatter suppression hardware for each detector is different, and high-frequency signals will be superimposed on the basis of cross-scattering. Even if the traditional scatter correction method is used to correct the cross-scattering signal, obvious ring-shaped or strip-shaped artifacts will still exist in the reconstructed image obtained, for example Figure 2 in the reconstructed image shown, obvious ring-shaped artifacts can still be seen at the positions indicated by the black arrows.

[0054] In response to this, in this embodiment, the signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system can be obtained in advance.

[0055] Among them, the dual-source or multi-source CT imaging system may be a medical imaging device integrated with two or more sets of X-ray tube-detector units, that is, the dual-source or multi-source CT imaging system includes at least two detectors, and each detector has a supporting X-ray tube. In the detection signal detected by the detector, the signal from the supporting X-ray tube can be called the primary ray signal, and the signal from the non-supporting X-ray tube can be called the cross-scattering signal. Among them, each detector includes multiple channels. The channels of the detector can be understood as independent units or paths in the detector for receiving and measuring X-ray signals. Different channels correspond to different spatial positions. By the coordinated work of multiple channels, the incident position of the X-ray on the detector plane can be accurately determined, so as to provide accurate spatial information for subsequent image reconstruction and construct a tomographic image with high resolution.

[0056] It can be understood that due to the differences in the scatter suppression efficiency of the scatter suppression hardware for each detector, for different channels of different detectors in a dual-source or multi-source CT imaging system, the intensity of the detected cross-scatter signals may also vary. In some cases, the intensity of the cross-scatter signal received by a certain channel in the detector may change significantly relative to the intensity of the cross-scatter signals detected by other channels, that is, the intensity of the received cross-scatter signal undergoes a signal intensity jump.

[0057] In some embodiments, the signal intensity jump can be understood as follows: due to the differences in scatter suppression efficiency, the intensity of the cross-scatter signal detected by one channel in the detector jumps relative to the intensity of other cross-scatter signals used as a reference; in terms of the observation results, since the detected signal received by the channel may include cross-scatter signals from non-matching X-ray tubes, the differences in scatter suppression efficiency will result in different intensities of the cross-scatter signals received by different channels, so that when the cross-scatter signal intensity jumps, the intensity of the detected signal received by the detector channel changes significantly, for example, the change amplitude of the detected signal intensity is greater than a threshold value. Among them, the jump may include an upward jump, for example, rising from the detected signal intensity A to the detected signal intensity B. Of course, the jump may also include a downward jump, for example, dropping from the detected signal intensity C to the detected signal intensity B. In some examples, the signal intensity used as a reference may be the intensity of the detected signal detected by an adjacent channel. For example, for the intensity of the cross-scatter signal detected by a certain channel in the detector, the intensity of the cross-scatter signal detected by the adjacent channel connected to it in terms of spatial position can be used as a reference to determine whether a signal intensity jump occurs.

[0058] In this embodiment, a dual-source or multi-source CT imaging system can be used to scan different phantoms at the same position and / or the same phantom at different positions. In some examples, the phantom used during scanning may refer to a phantom whose chemical parameters and / or physical parameters are the same as or similar to those of the object to be measured during the actual scanning process. For example, if the object to be measured during actual scanning is a living body, the phantom can be a water film.

[0059] During the process of scanning the phantom, each channel of multiple detectors in the dual-source or multi-source CT imaging system can collect signals and determine the intensity of the detected signal corresponding to each channel, thereby obtaining multiple detected signal intensities. Furthermore, by analyzing the multiple detected signal intensities detected, the signal intensity jump relationship can be determined, where the signal intensity jump relationship can characterize the jump amplitude of the detector channel at each cross-scatter signal intensity.

[0060] In some embodiments, the detected signal intensity can be recorded in the form of a curve. For example, when the detected signal intensity is the cross-scatter signal intensity, it can be recorded asFigure 3 Recorded in the manner shown Figure 3 A set of cross-scattering signal curves is provided. In the Figure 3 coordinate system, the abscissa is the number of the detector channels, and the ordinate is the signal intensity. It can be understood that by analyzing one or more cross-scattering signal curves, the signal intensity jump relationship can be determined.

[0061] S102. Obtain the scattering estimation information obtained by scanning the object to be measured with a dual-source or multi-source CT imaging system.

[0062] In practical applications, the object to be measured can be scanned with a dual-source or multi-source CT imaging system, and the scattering estimation is performed and the scattering estimation information is obtained according to the detection signal intensities detected by the channels of multiple detectors in the system. Among them, the scattering estimation information can characterize the cross-scattering signal intensities detected by multiple channels, such as the scattering distribution; for example, the scattering photon energy, spatial distribution, intensity and other elements received by each channel in the dual-source or multi-source CT imaging system can be quantitatively estimated, and the obtained scattering distribution is used as the scattering estimation information.

[0063] S103. Correct the scattering estimation information according to the signal intensity jump relationship.

[0064] Since the signal intensity jump relationship can characterize the signal intensity jump situations of multiple channels caused by the difference in cross-scattering suppression efficiency, and at the same time, the scattering estimation information can characterize the cross-scattering signal intensities detected by multiple channels, accordingly, the scattering estimation information of the object to be measured can be adjusted and corrected according to the signal intensity jump relationship to obtain the updated scattering estimation information. Subsequently, signal correction and image reconstruction can be performed according to the scattering estimation information.

[0065] In the above cross-scattering correction method, it can be applied to a dual-source or multi-source CT imaging system. The dual-source or multi-source CT imaging system includes at least two detectors, and each detector contains multiple channels. In this method, on the one hand, the signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system can be obtained. This signal intensity jump relationship is determined based on the detection signal intensities of each channel of multiple detectors. The detection signal intensity is obtained by scanning different phantoms and / or the same phantom at different positions by the detectors; on the other hand, the scattering estimation information obtained by scanning the object to be measured by the dual-source or multi-source CT imaging system can be obtained. Furthermore, the scattering estimation information can be corrected according to the signal intensity jump relationship. In this embodiment, by pre-obtaining the signal intensity jump relationship, the signal intensity jump situation of the cross-scattering signal intensity detected by the detector channels under different signal intensities can be established, accurately reflecting the influence on the detector channels under different scattering suppression efficiencies. Therefore, when the scattering estimation information of the object to be measured is obtained, by using the signal intensity jump relationship to correct this scattering estimation information, the corrected result can more accurately reflect the actual situation of the detector channels receiving the cross-scattering signal considering the difference in scattering suppression efficiency, removing the abnormal high-frequency signals superimposed due to the difference in scattering suppression efficiency, and improving the accuracy of the cross-scattering correction effect.

[0066] In an exemplary embodiment, the detection signal intensity obtained when scanning a phantom using a dual-source or multi-source CT imaging system includes the cross-scattering signal intensity; correspondingly, as Figure 4 shown, in step S101, the following steps may be included:

[0067] S401, determine the jump amplitude and the target signal intensity before the jump when the signal intensity of the channel undergoes a signal intensity jump according to multiple cross-scattering signal intensities.

[0068] In a specific implementation, a dual-source or multi-source CT imaging system can be used to scan a phantom, and according to the detection signals obtained by each channel in the dual-source or multi-source CT imaging system, the cross-scattering signal intensity of the cross-scattering signal can be determined. In some embodiments, the signal type of the detection signal received by each detector can be identified to determine the cross-scattering signal intensity corresponding to the cross-scattering signal from a non-matching X-ray tube; in some other embodiments, the detection signal intensity detected by the detector can also be used as the cross-scattering signal intensity when the matching X-ray tube is turned off and the non-matching X-ray tube is turned on.

[0069] Then, by comparing and analyzing each cross-scattering signal intensity, the jump amplitude when the cross-scattering signal of the channel undergoes a signal intensity jump and the corresponding target signal intensity before the jump can be determined.

[0070] Among them, the jump amplitude can be understood as the change value of the cross-scattering signal intensity. The target signal intensity is the signal intensity before the signal intensity jump occurs, which can be understood as the starting point of the signal intensity jump. For example, if it significantly rises from signal intensity A to signal intensity B, then signal intensity A can be called the target signal intensity, and the change value between signal intensity B and signal intensity A is called the jump amplitude. It can be understood that for the same channel, there can be one or more target signal intensities, and the corresponding jump amplitudes can also be different under different target signal intensities.

[0071] Taking Figure 3 the cross-scattering signal curve shown as an example, it can be determined that there are relatively obvious jumps near channels 190 and 200. Figure 3 Among the three regions circled in

[0072] S402, perform function fitting based on the jump amplitude and target signal intensity corresponding to each signal intensity jump, and obtain the signal intensity jump relationship according to the fitting result.

[0073] Furthermore, after determining the jump amplitude and target signal intensity corresponding to each signal intensity jump, function fitting can be performed to determine the correlation relationship between the target signal intensity and the jump amplitude, so as to obtain the signal intensity jump relationship according to the fitting result of the function.

[0074] For example, multiple discrete points determined by "target signal intensity - jump amplitude" can be obtained according to step S401, and then, with the signal intensity as the abscissa and the jump amplitude as the ordinate, a linear function fitting is performed on the multiple discrete points to obtain the signal intensity jump relationship as shown in Figure 5

[0075] In some exemplary embodiments, the signal intensity jump relationships of different channels may be different. In this regard, for each channel where a signal intensity jump has occurred, function fitting can be performed according to the jump amplitude and target signal intensity corresponding to the signal intensity jump of this channel to obtain the signal intensity jump relationship corresponding to this channel.

[0076] In this embodiment, by pre-determining the jump amplitude and target signal intensity corresponding to each signal intensity jump according to the signal intensities of multiple cross-scattering signals and constructing the signal intensity jump relationship, the signal intensity jump differences under different cross-scattering signal intensities can be established. Furthermore, when obtaining the scattering estimation information of the actual object to be measured, the scattering estimation information can be quickly corrected through this signal intensity jump relationship, reducing the high-frequency scattering signals caused by the scattering suppression efficiency differences, and significantly improving the correction speed and correction effect of cross-scattering correction.

[0077] ​In an exemplary embodiment, in step S401, determining the jump amplitude when the signal intensity of the channel jumps according to multiple cross-scattering signal intensities may include the following steps:

[0078] Obtain cross-scattering signal curves corresponding to multiple cross-scattering signal intensities; perform curve smoothing on the cross-scattering signal curves to obtain the smoothed cross-scattering signal curves; determine the jump amplitude when the signal intensity of the channel jumps according to the difference between the cross-scattering signal curves and the smoothed cross-scattering signal curves.

[0079] In practical applications, cross-scattering signal curves can be generated according to multiple cross-scattering signal intensities, and then curve smoothing is performed on the cross-scattering signal curves. In some examples, it can be implemented by using methods for smoothing signal curves in related technologies, such as any one of the moving average method, Gaussian filtering method, median filtering method, and wavelet transform method.

[0080] By performing curve smoothing on the cross-scattering signal curves, one or more mutations existing in the curves can be eliminated. Furthermore, by comparing the cross-scattering signal curves with the smoothed cross-scattering signal curves to determine the difference between the curves, the jump amplitude corresponding to the signal intensity jump of the detector channel and the target signal intensity before the jump can be determined.

[0081] In some embodiments, when the computing resources are insufficient (for example, the amount of computing resources does not meet the preset conditions), the jump amplitude and the target signal intensity corresponding to the signal intensity jump of the detector channel can be determined by the method provided in this embodiment.

[0082] In this embodiment, by performing curve smoothing on the cross-scattering signal curves and then determining the jump amplitude corresponding to the signal intensity jump of the detector channel according to the difference between the cross-scattering signal curves and the smoothed cross-scattering signal curves, the jump amplitude corresponding to the signal intensity jump can be quickly determined on the basis of the originally obtained cross-scattering signals. The jump amplitude corresponding to the signal intensity jump can be quickly obtained without excessive complex processing, significantly improving the acquisition efficiency of the signal intensity jump relationship.

[0083] In an exemplary embodiment, determining the jump amplitude when the signal intensity of the channel jumps according to the difference between the cross-scattering signal curves and the smoothed cross-scattering signal curves includes:

[0084] Subtract the smoothed cross-scattering signal curves from the cross-scattering signal curves to obtain the signal intensity differences corresponding to each channel; for the signal intensity differences corresponding to each channel, when the signal intensity differences meet the signal intensity jump condition, determine the signal intensity differences as the jump amplitudes corresponding to the channels.

[0085] In a specific implementation, the cross-scattering signal curves before and after smoothing can be subtracted, that is, a difference operation is performed, and then the signal intensity difference corresponding to each channel can be obtained according to the difference result. During the difference operation, for each channel, the signal intensity on the original cross-scattering signal curve of this channel can be subtracted from the signal intensity on the smoothed cross-scattering signal curve to obtain the signal intensity difference of this channel.

[0086] Then, for the signal intensity difference corresponding to each channel, it can be determined whether the signal intensity difference satisfies the signal intensity jump condition. For example, whether the signal intensity difference is greater than a threshold, or whether the difference between the signal intensity difference and the signal intensity difference of an adjacent signal detection unit is greater than a threshold. If so, it can be determined that the signal intensity jump condition is satisfied. Furthermore, when the signal intensity difference satisfies the signal intensity jump condition, the signal intensity difference can be determined as the jump amplitude corresponding to the signal intensity jump.

[0087] In this embodiment, by subtracting the cross-scattering signal curve from the smoothed cross-scattering signal curve, and according to the obtained signal intensity difference and the signal intensity jump condition, the jump amplitude corresponding to the signal intensity jump is determined, and the signal intensity jump and its jump amplitude existing in the curve can be quickly identified through simple subtraction processing, while saving the computational resource overhead and improving the construction speed of the scattering jump relationship.

[0088] In an exemplary embodiment, in step S401, the following steps may be included:

[0089] Obtain a phantom reconstruction image according to multiple detection signal intensities; perform pixel value smoothing processing on the abnormal region with abnormal pixel values in the phantom reconstruction image to obtain a smoothed reconstruction image; determine the jump amplitude and the target signal intensity before the jump when the signal intensity of the channel jumps according to the difference between the phantom reconstruction image and the smoothed reconstruction image; perform function fitting according to the jump amplitude and the target signal intensity corresponding to each signal intensity jump, and obtain the signal intensity jump relationship according to the fitting result.

[0090] In practical applications, image reconstruction can be performed according to the detection signals collected by each channel in a dual-source or multi-source CT imaging system to obtain a reconstruction image of the phantom, that is, the phantom reconstruction image. At this time, the phantom reconstruction image may include an abnormal region with abnormal pixel values, such as image artifacts, and this abnormal pixel value can be understood as an abnormality caused by the difference in scattering suppression efficiency.

[0091] Then, pixel value smoothing processing can be performed on the abnormal region with abnormal pixel values in the phantom reconstruction image to obtain a smoothed reconstruction image. For example, it can be reconstructed according to the detection signal intensity to obtain Figure 2For the reconstructed phantom image shown, by analyzing the image content, the area where the annular artifact is located can be visually determined. This area can be used as an abnormal area with abnormal pixel values, and pixel value smoothing processing is performed to obtain the smoothed reconstructed image.

[0092] Furthermore, the reconstructed phantom image and the smoothed reconstructed image can be compared to determine the difference between them, and based on this difference, the jump amplitude and the target signal intensity corresponding to the signal intensity jump occurring in the detector channel can be determined.

[0093] Since it is necessary to determine the jump amplitude corresponding to the signal intensity jump through image reconstruction, in some embodiments, when the computing resources are sufficient (for example, the resource amount of the computing resources meets the preset conditions) or the requirement for the recognition accuracy of the signal intensity jump is higher than the preset conditions, the method provided in this embodiment can be used to determine the jump amplitude corresponding to the signal intensity jump occurring in the detector channel.

[0094] Then, function fitting can be performed according to the jump amplitude and the target signal intensity corresponding to each signal intensity jump, and the signal intensity jump relationship can be obtained based on the fitting result. The fitting process of the signal intensity jump relationship can be the same as that in one or more of the foregoing embodiments and will not be elaborated here.

[0095] In this embodiment, by performing pixel value smoothing processing on the abnormal area with abnormal pixel values in the reconstructed phantom image to obtain the smoothed reconstructed image, and based on the difference between the reconstructed phantom image and the smoothed reconstructed image, the jump amplitude and the target signal intensity before the jump are determined, and then function fitting is performed. It is possible to accurately back-infer the signal intensity jumps existing in multiple cross-scattering signals through the reconstructed phantom image. It can be understood that by analyzing the reconstructed phantom image, the area with abnormal pixel values can be more accurately and visually identified through the image content, thereby effectively improving the recognition accuracy of the signal intensity jump.

[0096] In an exemplary embodiment, each pixel point in the reconstructed phantom image corresponds to a corresponding channel. Specifically, different channels correspond to different spatial positions. Through the collaborative work of multiple channels, the incident position of the X-ray on the detector plane can be accurately determined, and accurate spatial information can be provided for subsequent image reconstruction. Correspondingly, each pixel point in the reconstructed phantom image can correspond to a corresponding channel.

[0097] In this embodiment, to determine the jump amplitude when the signal intensity jump occurs in the channel according to the difference between the reconstructed phantom image and the smoothed reconstructed image, the following steps can be included:

[0098] Subtract the reconstructed phantom image from the smoothed reconstructed image to obtain the pixel value difference corresponding to each pixel in the reconstructed phantom image. For each pixel value difference, when the pixel value difference meets the pixel value jump condition, determine the jump amplitude when the signal intensity jumps in the channel corresponding to the pixel according to the pixel value difference.

[0099] In practical applications, the reconstructed phantom image can be subtracted from the smoothed reconstructed image. This process can be understood as follows: for each pixel in the reconstructed phantom image, a corresponding pixel can be determined in the smoothed reconstructed image, and then the pixel value difference between this pixel in the reconstructed image and the corresponding pixel in the smoothed reconstructed image is obtained, thereby obtaining the pixel value difference corresponding to each pixel in the reconstructed phantom image.

[0100] Then, it can be determined whether the pixel value difference meets the pixel value jump condition, such as whether the pixel value difference is greater than 0 or other preset thresholds. If so, it is determined that the pixel value jump condition is met. Furthermore, for each pixel value difference that meets the pixel value jump condition, the jump amplitude when the signal intensity jumps in the channel corresponding to the pixel can be determined according to each pixel value difference. For example, for each pixel that meets the pixel value jump condition, the artifact region to be corrected in the reconstructed phantom image can be obtained according to each pixel and its pixel value. By performing a forward projection simulation on this artifact region, the pixel value jump size corresponding to the artifact region can be seen. Since the pixel value size corresponds to the signal intensity size and the pixel corresponds to the detector channel, the jump amplitude corresponding to the signal intensity jump occurring in the detector channel can be determined according to the pixel value jump size.

[0101] In this embodiment, by subtracting the reconstructed phantom image from the smoothed reconstructed image, the signal detection unit where the signal intensity jumps and its jump amplitude can be accurately identified according to the pixel value difference.

[0102] In an exemplary embodiment, multiple detection signal intensities are obtained through the following steps: scanning different phantoms placed at the same position, and determining multiple detection signal intensities according to the signal intensities of the detection signals collected by each channel of multiple detectors during each scan; and / or scanning the same phantom placed at different positions, and determining multiple detection signal intensities according to the signal intensities of the detection signals collected by each channel of multiple detectors during each scan.

[0103] In practical applications, a dual-source or multi-source CT imaging system can scan multiple preset phantoms of different sizes placed at the same position, or a dual-source or multi-source CT imaging system can scan multiple preset phantoms of the same size placed at different positions. Of course, it can also be a combination of the above two methods.

[0104] For the first method, phantoms of different sizes can be placed at the same position and scanned separately. For example, in each scan, a preset of a certain size is placed at the isocenter (i.e., isocenter, also known as iso center) of the multi-source CT system. For the scan of each size phantom, detection signals can be collected by multiple channels in the system, and then the signal intensity of the detection signals can be determined.

[0105] For the second method, phantoms of the same size can be placed at different positions and scanned separately. For example, in this scan, a phantom of a specified size can be placed under eccentricity 1 of the multi-source CT system for scanning. In the next scan, the same specified-size phantom can be placed under eccentricity 2 for scanning, where eccentricity 1 is different from eccentricity 2. For the scan of each position, detection signals can be collected by multiple channels in the system, and the corresponding signal intensity can be determined.

[0106] For example, the cross-scattering signal intensities of phantoms of different sizes or the same phantom with different eccentricities can be recorded separately, resulting in Figure 3 3 curves as shown. Each curve corresponds to one measurement / scan, such as a phantom of a different size or a different eccentricity.

[0107] In this embodiment, the detection signal intensities received by multiple channels in different situations can be obtained. By analyzing the multiple detection signal intensities, a refined signal intensity jump relationship can be constructed for each channel, and the difference in the jump amplitude at different signal intensities can be established. For different signal intensities involved in the scattering estimation information of different phantoms, the corresponding jump amplitude can be quickly and accurately determined and corrected, providing reliable support for correcting the scattering estimation information according to the signal intensity jump relationship in the subsequent process.

[0108] In an exemplary embodiment, in step S103, correcting the scattering estimation information according to the signal intensity jump relationship may include the following steps:

[0109] Determine the cross-scattering signal intensity of each channel according to the scattering estimation information; determine the jump amplitude of the cross-scattering signal intensity of the channel according to the signal intensity jump relationship; adjust each cross-scattering signal intensity in the scattering estimation information according to each jump amplitude to obtain the corrected scattering estimation information.

[0110] In a specific implementation, the scattering estimation information may record the cross-scattering signal intensities of multiple channels, and the cross-scattering signal intensities of each channel can be determined based on the scattering estimation information. Since the pre-acquired signal intensity jump relationship can record the target signal intensity before the signal intensity jump occurs and the jump amplitude of the signal intensity when the signal intensity jump occurs, therefore, by querying the signal intensity jump relationship, the jump amplitude of the cross-scattering signal intensity of each channel can be determined. For example, if the cross-scattering signal intensity of a certain channel recorded in the scattering estimation information is x1, the jump amplitude corresponding to the signal intensity x1 can be determined according to the signal intensity jump relationship.

[0111] Then, for each cross-signal intensity in the scattering estimation information, it can be adjusted according to its corresponding jump amplitude, so that the corrected scattering estimation information can be obtained. For example, during image reconstruction, the scattering distribution estimated for the real object to be measured can be calculated according to the signal intensity jump relationship, and then the updated scattering distribution can be obtained. Subsequently, the detection signal received by the detector can be corrected according to the updated scattering distribution, so as to correct the ring artifact caused by high-frequency scattering.

[0112] In this embodiment, according to the signal intensity jump relationship, the cross-scattering signal intensities of different channels can be differentially adjusted accurately and efficiently, effectively improving the accuracy and reliability of the scattering estimation information, and providing a reliable data basis for subsequent image reconstruction.

[0113] In order to enable those skilled in the art to better understand the above steps, the following exemplarily illustrates the embodiments of the present application through an example, but it should be understood that the embodiments of the present application are not limited thereto.

[0114] As Figure 6 shown, the following steps may be included:

[0115] S601, respectively scan phantoms of different sizes placed at the same position, and obtain multiple detection signal intensities according to the signal intensities of the detection signals collected by multiple channels each time, and / or, respectively scan the same phantom placed at different positions, and obtain multiple detection signal intensities according to the signal intensities of the detection signals collected by multiple channels each time.

[0116] S602, determine the jump amplitude and the target signal intensity before the jump when the signal intensity of the channel jumps according to the cross-scattering signal intensity in the multiple detection signal intensities.

[0117] S603. Obtain a phantom reconstruction image based on multiple detection signal intensities, perform pixel value smoothing processing on the abnormal regions with abnormal pixel values in the phantom reconstruction image to obtain a smoothed reconstruction image, and determine the jump amplitude and the target signal intensity before the jump when the signal intensity of the channel jumps according to the difference between the phantom reconstruction image and the smoothed reconstruction image.

[0118] S604. Perform function fitting according to the jump amplitude and the target signal intensity corresponding to each signal intensity jump, and obtain a signal intensity jump relationship according to the fitting result.

[0119] S605. Correct the scatter estimation information of the object to be measured according to the signal intensity jump relationship.

[0120] It should be understood that although the steps in the flowcharts involved in the above-described embodiments are sequentially shown according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless there is a clear description in this article, the execution of these steps has no strict order limit, and these steps can be executed in other orders. Moreover, at least a part of the steps in the flowcharts involved in the above-described embodiments may include multiple steps or multiple stages. These steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be executed alternately or alternately with at least a part of other steps or steps or stages in other steps.

[0121] Based on the same inventive concept, an embodiment of the present application also provides a cross-scatter correction device for implementing the above-mentioned cross-scatter correction method. The solution provided by this device to solve the problem is similar to the solution described in the above method. Therefore, the specific limitations in one or more cross-scatter correction device embodiments provided below can refer to the limitations on the cross-scatter correction method in the above text, and will not be repeated here.

[0122] In an exemplary embodiment, as Figure 7 shown, a cross-scatter correction device is provided, which is applied to a dual-source or multi-source CT imaging system. The dual-source or multi-source CT imaging system includes at least two detectors, and each detector includes multiple channels. The device includes:

[0123] A relationship acquisition module 701, configured to acquire a signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system; the signal intensity jump relationship is determined according to the detection signal intensities of each channel of multiple detectors, and the detection signal intensity is obtained by scanning different phantoms and / or the same phantom at different positions by the detectors;

[0124] A scatter estimation acquisition module 702, configured to acquire scatter estimation information obtained by scanning a subject to be measured by the dual-source or multi-source CT imaging system;

[0125] A correction module 703, configured to correct the scatter estimation information according to the signal strength jump relationship.

[0126] In one embodiment, the detection signal strength includes cross-scatter signal strength, and the relationship acquisition module 701 is configured to:

[0127] Determine a jump amplitude and a target signal strength before the jump when the signal strength of the channel jumps according to a plurality of the cross-scatter signal strengths;

[0128] Perform function fitting according to the jump amplitude and the target signal strength corresponding to each signal strength jump, and obtain the signal strength jump relationship according to the fitting result.

[0129] In one embodiment, the relationship acquisition module 701 is configured to:

[0130] Obtain cross-scatter signal curves corresponding to a plurality of the cross-scatter signal strengths;

[0131] Perform curve smoothing processing on the cross-scatter signal curves to obtain smoothed cross-scatter signal curves;

[0132] Determine a jump amplitude when the signal strength of the channel jumps according to the difference between the cross-scatter signal curves and the smoothed cross-scatter signal curves.

[0133] In one embodiment, the relationship acquisition module 701 is configured to:

[0134] Subtract the cross-scatter signal curves from the smoothed cross-scatter signal curves to obtain signal strength differences corresponding to each channel;

[0135] For the signal strength difference corresponding to each channel, when the signal strength difference satisfies the signal strength jump condition, determine the signal strength difference as the jump amplitude corresponding to the channel.

[0136] In one embodiment, the relationship acquisition module 701 is configured to:

[0137] Obtain a phantom reconstruction image according to a plurality of the detection signal strengths;

[0138] Perform pixel value smoothing processing on an abnormal region with abnormal pixel values in the phantom reconstruction image to obtain a smoothed reconstruction image;

[0139] Determine the jump amplitude and the target signal intensity before the jump when the signal intensity of the channel jumps according to the difference between the image reconstructed from the phantom and the smoothed reconstructed image;

[0140] Perform function fitting according to the jump amplitude and the target signal intensity corresponding to each signal intensity jump, and obtain the signal intensity jump relationship according to the fitting result.

[0141] In one embodiment, each pixel point in the image reconstructed from the phantom corresponds to a corresponding channel;

[0142] The relationship acquisition module 701 is used for:

[0143] Subtract the smoothed reconstructed image from the image reconstructed from the phantom to obtain the pixel value difference corresponding to each pixel point in the image reconstructed from the phantom;

[0144] For the pixel value difference corresponding to each pixel point, when the pixel value difference satisfies the pixel value jump condition, determine the jump amplitude when the signal intensity of the channel corresponding to the pixel point jumps according to the pixel value difference.

[0145] In one embodiment, the correction module 703 is used for:

[0146] Determine the cross-scattering signal intensity of each channel according to the scattering estimation information;

[0147] Determine the jump amplitude of the cross-scattering signal intensity of the channel according to the signal intensity jump relationship;

[0148] Adjust each cross-scattering signal intensity in the scattering estimation information according to each jump amplitude to obtain the corrected scattering estimation information.

[0149] Each module in the above cross-scattering correction device can be implemented in whole or in part by software, hardware and their combination. Each of the above modules can be embedded in the processor of the computer device in hardware form or be independent of it, or can be stored in the memory of the computer device in software form, so that the processor can call and execute the operations corresponding to each of the above modules.

[0150] In an exemplary embodiment, a computer device is provided. The computer device can be a server, and its internal structure diagram can be as Figure 8As shown in the figure. The computer device includes a processor, a memory, an input / output interface (Input / Output, abbreviated as I / O), and a communication interface. Among them, the processor, the memory, and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store the signal data collected by the signal detection unit. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with external terminals through a network connection. When the computer program is executed by the processor, it implements a cross-scattering correction method.

[0151] In an exemplary embodiment, a computer device is provided. The computer device can be a terminal, and its internal structure diagram can be as Figure 9 shown in the figure. The computer device includes a processor, a memory, an input / output interface, a communication interface, a display unit, and an input device. Among them, the processor, the memory, and the input / output interface are connected through a system bus, and the communication interface, the display unit, and the input device are connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with external terminals in a wired or wireless manner. The wireless manner can be implemented through WIFI, a mobile cellular network, near field communication (Near Field Communication, NFC), or other technologies. When the computer program is executed by the processor, it implements a cross-scattering correction method. The display unit of the computer device is used to form a visually visible picture, which can be a display screen, a projection device, or a virtual reality imaging device. The display screen can be a liquid crystal display screen or an electronic ink display screen. The input device of the computer device can be a touch layer covering the display screen, or a button, a trackball, or a touchpad provided on the computer device housing, or an external keyboard, touchpad, or mouse, etc.

[0152] Those skilled in the art can understand that Figure 8 and Figure 9The structure shown is only a block diagram of some structures related to the solution of this application, and does not constitute a limitation on the computer device to which the solution of this application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine some components, or have different component arrangements.

[0153] In one embodiment, a computer device is provided, including a memory and a processor. A computer program is stored in the memory, and when the processor executes the computer program, the steps in the above method embodiments are implemented.

[0154] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0155] In one embodiment, a computer program product is provided, including a computer program. When the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0156] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use, and processing of relevant data need to comply with relevant regulations.

[0157] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, database, or other medium used in the embodiments provided in this application can include at least one of non-volatile memory and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc. The databases involved in the embodiments provided in this application can include at least one of relational databases and non-relational databases. Non-relational databases can include distributed databases based on blockchain, etc., without limitation. The processors involved in the embodiments provided in this application can be general-purpose processors, central processors, graphics processors, digital signal processors, programmable logic devices, data processing logics based on quantum computing, artificial intelligence (AI) processors, etc., without limitation.

[0158] The technical features of the above embodiments can be combined arbitrarily. For the sake of brevity of description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered to be within the scope recorded in this application.

[0159] The above-described embodiments merely represent several implementation manners of the present application. The description thereof is relatively specific and detailed, but it should not be construed as a limitation on the patent scope of the present application. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present application, several modifications and improvements can still be made, and these all fall within the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the appended claims.

Claims

1. A cross-scattering correction method, characterized in that, Applied to a dual-source or multi-source CT imaging system, the dual-source or multi-source CT imaging system includes at least two detectors, and each of the detectors includes a plurality of channels. The method includes: Obtaining a signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system; the signal intensity jump relationship is determined according to the detection signal intensities of each channel of a plurality of the detectors, and the detection signal intensity is obtained by scanning different phantoms and / or the same phantom at different positions by the detectors; Obtaining scatter estimation information obtained by scanning a subject to be detected by the dual-source or multi-source CT imaging system; Correcting the scatter estimation information according to the signal intensity jump relationship.

2. The method according to claim 1, wherein The detection signal intensity includes a cross-scatter signal intensity. The obtaining of the signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system includes: Determining a jump amplitude and a target signal intensity before the jump when the signal intensity of the channel jumps according to a plurality of the cross-scatter signal intensities; Performing function fitting according to the jump amplitude and the target signal intensity corresponding to each signal intensity jump, and obtaining the signal intensity jump relationship according to the fitting result.

3. The method according to claim 2, wherein Determining the jump amplitude when the signal intensity of the channel jumps according to a plurality of the cross-scatter signal intensities includes: Obtaining cross-scatter signal curves corresponding to a plurality of the cross-scatter signal intensities; Performing curve smoothing processing on the cross-scatter signal curves to obtain smoothed cross-scatter signal curves; Determining the jump amplitude when the signal intensity of the channel jumps according to the difference between the cross-scatter signal curves and the smoothed cross-scatter signal curves.

4. The method according to claim 3, characterized in that, The determining of the jump amplitude when the signal intensity of the channel jumps according to the difference between the cross-scatter signal curves and the smoothed cross-scatter signal curves includes: Subtracting the cross-scatter signal curves from the smoothed cross-scatter signal curves to obtain signal intensity differences corresponding to each of the channels; For each signal intensity difference corresponding to the channel, when the signal intensity difference satisfies a signal intensity jump condition, determining the signal intensity difference as the jump amplitude corresponding to the channel.

5. The method according to claim 1, characterized in that The obtaining of the signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system includes: Obtaining a phantom reconstruction image according to a plurality of the detection signal intensities; Performing pixel value smoothing processing on an abnormal region with abnormal pixel values in the phantom reconstruction image to obtain a smoothed reconstruction image; Determining a jump amplitude and a target signal intensity before the jump when the signal intensity of the channel jumps according to the difference between the phantom reconstruction image and the smoothed reconstruction image; Performing function fitting according to the jump amplitude and the target signal intensity corresponding to each signal intensity jump, and obtaining the signal intensity jump relationship according to the fitting result.

6. The method according to claim 5, wherein Each pixel point in the phantom reconstruction image corresponds to a corresponding one of the channels; The determining of the jump amplitude when the signal intensity of the channel jumps according to the difference between the phantom reconstruction image and the smoothed reconstruction image includes: Subtract the reconstructed phantom image from the smoothed reconstructed image to obtain the pixel value difference corresponding to each pixel point in the reconstructed phantom image; For each pixel value difference corresponding to each pixel point, when the pixel value difference satisfies the pixel value jump condition, determine the jump amplitude when the signal intensity of the corresponding channel of the pixel point jumps according to the pixel value difference.

7. The method according to any one of claims 1 to 6, characterized in that, The correcting the scatter estimation information according to the signal intensity jump relationship includes: Determine the cross-scatter signal intensity of each channel according to the scatter estimation information; Determine the jump amplitude of the cross-scatter signal intensity of the channel according to the signal intensity jump relationship; Adjust each cross-scatter signal intensity in the scatter estimation information according to each jump amplitude to obtain the corrected scatter estimation information.

8. A cross-scattering correction device, characterized in that Applied to a dual-source or multi-source CT imaging system, the dual-source or multi-source CT imaging system includes at least two detectors, and each detector includes a plurality of channels. The device includes: A relationship acquisition module, configured to acquire the signal intensity jump relationship corresponding to the dual-source or multi-source CT imaging system; the signal intensity jump relationship is determined according to the detection signal intensities of each channel of a plurality of the detectors, and the detection signal intensity is obtained by scanning different phantoms and / or the same phantom at different positions by the detector; A scatter estimation acquisition module, configured to acquire the scatter estimation information obtained by scanning a subject to be detected by the dual-source or multi-source CT imaging system; A correction module, configured to correct the scatter estimation information according to the signal intensity jump relationship.

9. A computer device, comprising a memory and a processor, the memory storing a computer program, characterized in that, When the processor executes the computer program, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images