Cross scatter correction method and device and computer equipment

Abstract

The invention relates to a cross scatter correction method and device and computer equipment, relates to the technical field of medical images, and can improve the accuracy of a cross scatter correction result. 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, and the method comprises the steps that respective detection signal curves of the detectors are obtained, and the detection signal curves are determined based on detection signals collected by the detectors for a scanned object; for each detector, determining a first convolution region corresponding to a detection signal curve according to boundary position features of the detection signal curve of the detector, and determining scattering estimation information corresponding to the detector according to a convolution result of the first convolution region and a preset convolution kernel; and determining a cross scattering correction result of each detector according to the scattering estimation information corresponding to each detector.

Description

Technical Field

[0001] The present application relates to the field of medical imaging technology, and in particular to a cross-scatter correction method, device, and computer equipment. 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 obtain corresponding scan data. In dual-source CT imaging systems, cross-scattering is a physical factor that seriously affects image quality. Cross-scattering refers to the scattering of photons by an object and then transmission to the side detectors.

[0003] In related technologies, when estimating and correcting the scattering signal in cross scattering, preset phantoms of different sizes and at different eccentric positions are scanned in advance to obtain scattering signals under different conditions. During actual imaging, the size and position of the actual phantom are approximated with the various recorded preset phantom sizes and preset phantom positions to determine the closest preset phantom and select its corresponding scattering signal, and then scattering estimation is performed based on the scattering signal.

[0004] However, the above method has high roughness and is difficult to obtain accurate cross-scattering correction results. Summary of the Invention

[0005] 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.

[0006] 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, and the method includes:

[0007] Acquire a detection signal curve of each of the plurality of detectors, wherein the detection signal curve is determined based on a detection signal acquired by each detector with respect to the scanned object;

[0008] For each of the detectors, determining a first convolution region corresponding to the detection signal curve based on boundary position characteristics of the detection signal curve of the detector, and determining scattering estimation information corresponding to the detector based on a convolution result of the first convolution region with a preset convolution kernel;

[0009] A cross scatter correction result of each detector is determined according to the scatter estimation information corresponding to each detector.

[0010] In one embodiment, determining a first convolution region corresponding to the detection signal curve according to the boundary position feature of the detection signal curve of the detector includes:

[0011] According to the boundary position feature of the detection signal curve of the detector, determine a first starting position and a first ending position in the detector channel direction of the detection signal curve, and determine a first edge according to the first starting position and the first ending position;

[0012] Determine a second edge with a length less than a preset curve amplitude threshold in the signal intensity direction of the detection signal curve;

[0013] According to the first edge and the second edge, determine a first convolution region corresponding to the detection signal curve.

[0014] In one embodiment, determining a first convolution region corresponding to the detection signal curve according to the first edge and the second edge includes:

[0015] Determine a rectangular region according to the first edge and the second edge;

[0016] Obtain a first convolution region corresponding to the detection signal curve according to the rectangular region.

[0017] In one embodiment, before determining a first convolution region corresponding to the detection signal curve according to the boundary position feature of the detection signal curve of the detector, it further includes:

[0018] Obtain a preset curve amplitude threshold; the curve amplitude threshold is less than the maximum amplitude of the detection signal curve;

[0019] Determine curve points in the detection signal curve whose amplitudes match the curve amplitude threshold from the detection signal curve;

[0020] Obtain the boundary position feature of the detection signal curve according to the abscissa of the curve points.

[0021] In one embodiment, the preset convolution kernel is determined by the following steps:

[0022] Obtain a primary ray signal and a cross-scattering signal collected by a target detector for a phantom; the target detector is any one of the multiple detectors;

[0023] Determine a second convolution region corresponding to the main ray signal curve according to boundary position characteristics of the main ray signal curve corresponding to the main ray signal, and determine scattering estimation information corresponding to the target detector according to the cross-scattering signal and known scatter parameters of the phantom;

[0024] Perform deconvolution processing on the scattering estimation information corresponding to the target detector according to the second convolution region to obtain a preset convolution kernel.

[0025] In one embodiment, the determining a second convolution region corresponding to the main ray signal curve according to boundary position characteristics of the main ray signal curve corresponding to the main ray signal includes:

[0026] Determine a second starting position and a second ending position in the detector channel direction of the main ray signal curve according to boundary position characteristics of the main ray signal curve corresponding to the main ray signal, and determine a third edge according to the second starting position and the second ending position;

[0027] Determine a fourth edge with a length less than a preset curve amplitude threshold in the signal intensity direction of the main ray signal curve;

[0028] Determine a second convolution region corresponding to the main ray signal curve according to the third edge and the fourth edge.

[0029] In one embodiment, the first convolution region includes signal intensities of multiple detector channels in the detector;

[0030] The determining scattering estimation information corresponding to the detector according to a convolution result of the first convolution region and a preset convolution kernel includes:

[0031] For each detector channel in the first convolution region, obtain a signal intensity deviation of the detector channel according to a multiplication result of the signal intensity of the detector channel and the preset convolution kernel;

[0032] Determine scattering estimation information corresponding to the detector according to a summation result of signal intensity deviations of the multiple detector channels.

[0033] In one embodiment, the determining a cross-scattering correction result for each detector according to scattering estimation information corresponding to each detector includes:

[0034] Obtain a signal intensity of the detection signal collected by each detector;

[0035] For each detector, adjust the signal intensity according to the scattering estimation information corresponding to the detector to obtain an adjusted signal intensity;

[0036] Based on the adjusted signal intensity of each of the detectors, a cross-scattering correction result for each of the detectors is obtained.

[0037] In a second aspect, the present application further provides a cross-scattering correction device, 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 the device includes:

[0038] A curve acquisition module, configured to acquire the detection signal curves of the multiple detectors respectively, where the detection signal curves are determined based on the detection signals acquired by each of the detectors for the scanned object;

[0039] A convolution module, configured to, for each of the detectors, determine a first convolution region corresponding to the detection signal curve according to the boundary position characteristics of the detection signal curve of the detector, and determine the scattering estimation information corresponding to the detector according to the convolution result of the first convolution region and a preset convolution kernel;

[0040] A correction result acquisition module, configured to determine the cross-scattering correction result of each of the detectors according to the scattering estimation information corresponding to each of the detectors.

[0041] 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 cross-scattering correction method described in any one of the above is implemented.

[0042] 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 cross-scattering correction method described in any one of the above is implemented.

[0043] 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 cross-scattering correction method described in any one of the above is implemented.

[0044] The above cross-scattering correction method, device, computer device, computer-readable storage medium, and computer program product are 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. In this method, detection signal curves of multiple detectors can be obtained. The detection signal curves are determined based on the detection signals collected by each detector for the scanned object. Then, for each detector, according to the boundary position characteristics of the detection signal curve of the detector, a first convolution region corresponding to the detection signal curve is determined, and according to the convolution result of the first convolution region and a preset convolution kernel, scattering estimation information corresponding to the detector is determined. Furthermore, the cross-scattering correction result of each detector can be determined according to the scattering estimation information corresponding to each detector. In this application, on the one hand, according to the convolution result of the first convolution region and the preset convolution kernel, the scattering estimation information corresponding to the detector is determined, so that in the process of estimating cross-scattering, it is not necessary to approximate the actual phantom situation to a limited number of preset phantom situations, but the scattering can be directly estimated according to the detection signal image corresponding to the scanned object and the determined first convolution region, ensuring that the scattering estimation information accurately matches the actual situation of the scanned object. On the other hand, in this embodiment, by determining the first convolution region corresponding to the detection signal curve according to the boundary position characteristics of the detection signal curve, the interference of the intensity change amplitude of the detection signal curve on the scattering estimation can be reduced. Therefore, the accuracy of cross-scattering estimation can be effectively improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] To more clearly illustrate the technical solutions in the embodiments of the present application or related technologies, the following will briefly introduce the drawings required for use in the description of the embodiments of the present application or related technologies. Obviously, the 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 drawings can also be obtained based on these drawings.

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

[0047] Figure 2a It is a schematic diagram of cross-scattering of a dual-source CT imaging system in an embodiment;

[0048] Figure 2b It is a schematic diagram of cross-scattering of another dual-source CT imaging system in an embodiment;

[0049] Figure 3 It is a schematic comparison diagram of detection signal curves in an embodiment;

[0050] Figure 4 It is a schematic flowchart of the steps for determining a first convolution region in an embodiment;

[0051] Figure 5a Schematic diagram of a first convolution region in an embodiment;

[0052] Figure 5b Schematic diagram of another first convolution region in an embodiment;

[0053] Figure 6 Schematic diagram of a shoulder structure in an embodiment;

[0054] Figure 7a Schematic diagram of a convolution process in an embodiment;

[0055] Figure 7b Schematic diagram of another convolution process in an embodiment;

[0056] Figure 8 Schematic flow diagram of another cross-scattering correction method in an embodiment;

[0057] Figure 9 Structural block diagram of a cross-scattering correction device in an embodiment;

[0058] Figure 10 Internal structure diagram of a computer device in an embodiment;

[0059] Figure 11 Internal structure diagram of another computer device in an embodiment. Detailed implementation manners

[0060] In order to make the objectives, technical solutions and advantages of the present application clearer and more understandable, the present application will be further described in detail below with reference to 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.

[0061] In one embodiment, as Figure 1 shown, a cross-scattering correction method is provided. In this embodiment, the method is described by taking its application to a dual-source or multi-source CT imaging system as an example. Among them, the dual-source or multi-source CT imaging system includes at least two detectors. It can be understood that the method can also be applied to a server, and can also be applied 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:

[0062] S101. Obtain the detection signal curves of multiple detectors respectively, where the detection signal curves are determined based on the detection signals collected by each detector for the scanned object.

[0063] In specific implementation, the scanned object can be scanned by a dual-source or multi-source CT imaging system.

[0064] Among them, a dual-source or multi-source CT imaging system can be a medical imaging device integrated with two or more X-ray tube-detector sets. The object to be scanned can be any object that needs to obtain CT images. Exemplarily, the object to be scanned can be a living body, tissue organs, or one or more of them.

[0065] During the process of scanning the object to be scanned using a dual-source or multi-source CT imaging system, multiple detectors in the dual-source or multi-source CT imaging system can be used for signal acquisition, and each detector among the multiple detectors can obtain a corresponding detection signal. The detection signal received by the detector can simultaneously include a primary ray signal and a cross-scattered signal. Of course, the detection signal can also only include the primary ray signal or only include the cross-scattered signal. For example, by performing type recognition on the signal received by the detector, the primary ray signal therein can be determined and used as the detection signal to be analyzed subsequently. Another example is that by performing corresponding closing and opening processes on the paired and unpaired X-ray tubes, the detector can be made to receive only the primary ray signal or only the cross-scattered signal.

[0066] Among them, the primary ray signal refers to a signal whose signal source comes from the paired, opposite-side X-ray tube, and the cross-scattered signal refers to a signal scattered from an unpaired, side X-ray tube. For example, in the dual-source CT imaging system shown in Figure 2a Figure, X-ray tube 1 is paired with the opposite-side detector 1, and X-ray tube 2 is paired with the opposite-side detector 2. When detectors 1 and 2 are simultaneously turned on, if X-rays are emitted from X-ray tube 2, a part of the rays can be received by detector 2, and this part of the signal is called the primary ray signal. Another part is cross-scattered and the corresponding signal is received by detector 1, and this signal is the cross-scattered signal. Another example is that in the dual-source CT imaging system shown in Figure 2b Figure, if X-rays are emitted from X-ray tube 1, then detector 2 can receive the cross-scattered signal whose signal source is X-ray tube 1.

[0067] For the detection signal obtained by each detector, a corresponding detection signal curve can be generated. In some examples, the detector channel can be used as the abscissa and the signal intensity of the detection signal can be used as the ordinate to generate the detection signal curve. The detection signal curve can represent the signal intensity corresponding to the detection signals received on multiple channels in the detector.

[0068] S102. For each detector, according to the boundary position characteristics of the detection signal curve of the detector, determine the first convolution region corresponding to the detection signal curve, and according to the convolution result of the first convolution region and a preset convolution kernel, determine the scattering estimation information corresponding to the detector.

[0069] In the related art, there are mainly two schemes for cross-scattering correction. One is to add detectors outside the field of view to monitor the scattering signal. At the same time, the scattering signal within the field of view is considered as a linear interpolation of the scattering signal outside the field of view. Thus, through the scattering signal monitored by the detectors outside the field of view, this scheme still meets the requirements of feasibility when the collimation is narrow. However, as the collimation increases, the scattering distribution presents a two-dimensional complex performance, and linear interpolation or even high-order interpolation is not sufficient to give a better scattering estimation result. The other scheme is to scan phantoms of different sizes and at different eccentric positions in advance to obtain the scattering signals in different situations. During actual imaging, by approximating the size and positioning of the scanned object with the sizes and positions of various phantoms for which the corresponding scattering signals have been recorded, the closest phantom is determined and its corresponding scattering signal is selected, and then the scattering estimation is performed based on this scattering signal. This scheme has a high roughness and poor correction effect.

[0070] In this regard, in this embodiment, the convolution region corresponding to the detection signal curve can be determined. For the convenience of distinction, this convolution region is referred to as the first convolution region. After obtaining the first convolution region, a convolution operation can be performed on the first convolution region using a preset convolution kernel to obtain the corresponding convolution result. Then, the scattering estimation information corresponding to the detector can be determined based on this convolution result, where the scattering estimation information can be a scattering distribution. For example, one or more elements such as the energy, spatial distribution, and intensity of the scattered photons received by the detector can be quantitatively estimated, and the obtained scattering distribution is used as the scattering estimation information. The process of performing the convolution operation on the first convolution region using the preset convolution kernel can be understood as modeling the spatial propagation of scattering through the preset convolution kernel, converting the physical scattering process into operations such as neighborhood weighted summation in signal processing, so as to achieve efficient scattering estimation and correction. Since in this embodiment, the scattering estimation can be performed based on the convolution operation, this scattering estimation process does not need to approximate the actual situation of the scanned object with the situations of a limited number of phantoms, but directly performs the scattering estimation according to the actual projection value of the scanned object (that is, the first convolution region corresponding to the detection signal curve), so the correction effect is better.

[0071] In some exemplary embodiments, when determining the convolution region, it can be determined according to the region enclosed by the signal curve and the horizontal axis in the coordinate system. The region enclosed in this way is related to the boundary of the signal curve and the change amplitude of the signal. For example, as Figure 3The two detected signal curves shown can be collected during different scanning processes. There are significant differences in the variation amplitudes of detected signal curve 1 and detected signal curve 2. When determining the convolution region, the convolution region of detected signal curve 1 is determined according to the region formed by detected signal curve 1 and the X-axis, and the convolution region of detected signal curve 2 is determined according to the region formed by detected signal curve 2 and the X-axis. However, it is found in practice that when estimating cross-scattering, the contribution of the signal variation amplitude to the cross-scattering estimation is less. If the convolution region for subsequent processing is determined considering the variation amplitude, it may increase the computational amount and may also affect the estimation accuracy of cross-scattering.

[0072] Based on this, in the process of determining the first convolution region in this embodiment, the detected signal curve of the detector can also be first analyzed for its boundary to obtain the boundary position feature reflecting the boundary position of the detected signal curve. Among them, the boundary position feature can be understood as the feature that can characterize the intersection position of the detected signal curve with the horizontal or vertical axis in the coordinate system, and it can be information that accurately reflects or roughly reflects the starting position and ending position of the detected signal curve in the detected signal image. Furthermore, the first convolution region for subsequent convolution operation processing can be determined according to the boundary position feature. The first convolution region in this application can change correspondingly with the change of the boundary region feature. When the boundary region feature remains unchanged, the first convolution region can also remain unchanged. For example, as Figure 3 shown by the two detected signal curves, although there are significant differences in the variation amplitudes of detected signal curve 1 and detected signal curve 2, since the boundary positions of detected signal curve 1 and detected signal curve 2 are the same, the first convolution regions determined respectively for detected signal curve 1 and detected signal curve 2 can be the same or relatively similar.

[0073] S103. Determine the cross-scattering correction result of each detector according to the scattering estimation information corresponding to each detector.

[0074] After obtaining the scattering estimation information corresponding to each detector through step S102, the cross-scattering correction result of each detector can be obtained according to the scattering estimation information corresponding to each detector. Specifically, for each detector, the signal intensity of the detected signal collected by the detector can be corrected according to the scattering estimation information corresponding to the detector, and the obtained correction result is used as the cross-scattering correction result of the detector. For example, after obtaining the scattering estimation information of detector 1 according to the detected signal curve of detector 1, the detected signal collected by detector 1 is corrected using the scattering estimation information of detector 1. Similarly, after obtaining the scattering estimation information of detector 2 according to the detected signal curve of detector 2, the detected signal collected by detector 2 is corrected using the scattering estimation information of detector 2.

[0075] The above cross-scattering correction method 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. In this method, the detection signal curves of the respective detectors can be obtained, and the detection signal curve is determined based on the detection signals collected by the detector for the scanned object. Then, for each detector, according to the boundary position characteristics of the detection signal curve of the detector, a first convolution region corresponding to the detection signal curve is determined, and according to the convolution result of the first convolution region and a preset convolution kernel, the scattering estimation information corresponding to the detector is determined. Furthermore, the cross-scattering correction result of each detector can be determined according to the scattering estimation information corresponding to each detector. In this embodiment, on the one hand, according to the convolution result of the first convolution region and the preset convolution kernel, the scattering estimation information corresponding to the detector is determined, so that in the process of estimating cross-scattering, it is not necessary to approximate the situation of the scanned object with the situation of a finite phantom, but the scattering estimation can be directly performed according to the detection signal image corresponding to the scanned object and the determined first convolution region, ensuring that the scattering estimation information accurately matches the actual situation of the scanned object. On the other hand, in this embodiment, by determining the first convolution region corresponding to the detection signal curve according to the boundary position characteristics of the detection signal curve, the interference of the intensity change amplitude of the detection signal curve on the scattering estimation can be reduced. Thus, the accuracy of cross-scattering estimation can be effectively improved.

[0076] In addition, the cross-scattering correction method provided in this embodiment can directly perform scattering estimation and correction on the raw data (i.e., the detection signal curve determined according to the detection signal), without the need to pre-reconstruct the image and then perform scattering estimation, saving computing resources and significantly improving the correction speed, making the cross-scattering correction speed faster.

[0077] In an exemplary embodiment, as Figure 4 shown, in step S102, according to the boundary position characteristics of the detection signal curve of the detector, determining the first convolution region corresponding to the detection signal curve may include the following steps:

[0078] S401, according to the boundary position characteristics of the detection signal curve of the detector, determine a first starting position and a first ending position in the detector channel direction of the detection signal curve, and determine a first edge according to the first starting position and the first ending position.

[0079] Specifically, the detection signal curve can be recorded in a two-dimensional coordinate system, and the two-dimensional coordinate system can be composed of a detector channel direction and a signal intensity direction that are perpendicular to each other, so that the intensity change of the detection signals received by different channels in the detector can be described by this coordinate system.

[0080] Meanwhile, the boundary position features of the detection signal curve may include the starting position and / or the ending position of the curve of the detection signal curve, or may be other positions that are a preset distance away from the starting position and / or the ending position of the curve.

[0081] Correspondingly, in this step, the starting position and the ending position can be determined in the detector channel direction according to the boundary position features of the detection signal curve. For the sake of distinction, the starting position and the ending position are referred to as the first starting position and the first ending position. Then, the corresponding region edge is determined according to the first starting position and the first ending position. For the sake of distinction, the region edge determined according to the first starting position and the first ending position is referred to as the first edge, and the first edge can limit the span of the first convolution region in the detector channel direction.

[0082] S402. Determine a second edge with a length less than a preset curve amplitude threshold in the signal intensity direction of the detection signal curve.

[0083] As mentioned above, when performing cross-scattering estimation and correction, it is mainly related to the boundary (or span) of the detection signal curve, and is not sensitive to the change amplitude of the detection signal curve in the signal intensity direction. Directly determining the convolution region according to the change amplitude of the detection signal curve in the signal intensity direction may introduce redundant resource consumption and affect the accuracy of the cross-scattering result. For this reason, in this embodiment, a curve amplitude threshold can be determined in advance, and the region above the curve amplitude threshold can be removed or the weight of the region above the curve amplitude threshold can be reduced. At the same time, according to the region below the curve amplitude threshold, the first convolution region is determined.

[0084] In some exemplary embodiments, the curve amplitude threshold can be an empirical value, which can be determined according to multiple test results or experimental results. For example, in the case of a known preset convolution kernel and target scattering estimation information as a reference, multiple curve amplitude thresholds can be set and the corresponding scattering estimation information can be calculated respectively according to each curve amplitude threshold. By determining the difference between each scattering estimation information and the target scattering estimation information, the curve amplitude threshold corresponding to the smallest difference is determined.

[0085] Furthermore, after obtaining the detection signal curve, the edge of the detection signal curve in the signal intensity direction can be determined according to the preset curve amplitude threshold. For the sake of distinction, the edge determined according to the curve amplitude threshold is referred to as the second edge, and the length of the second edge is less than the curve amplitude threshold. For example, in the signal intensity direction, the distance corresponding to the curve amplitude threshold can be determined as the second edge.

[0086] S403. Determine the first convolution region corresponding to the detection signal curve according to the first edge and the second edge.

[0087] After obtaining the first edge and the second edge, since the first edge defines the span of the first convolution region in the detector channel direction and the second edge defines the amplitude of the first convolution region in the signal intensity direction, the first convolution region can be determined based on the first edge and the second edge.

[0088] In this embodiment, by determining the first edge according to the first starting position and the first ending position in the detector channel direction, and determining the second edge with a length less than a preset curve amplitude threshold in the signal intensity direction, and then determining the first convolution region according to the first edge and the second edge, on the one hand, it can accurately delimit the horizontal range matching the boundary of the detection signal curve in the detector channel direction, facilitating subsequent signal analysis focusing on specific spatial positions (such as specific detector channels). On the other hand, it can only focus on the signal part with amplitude changes within a certain range, exclude irrelevant signal interference, and reduce the influence of the change amplitude of the detection signal curve on the cross-scattering estimation. Compared with convolving the entire range enclosed by the detection signal curve and the abscissa, determining the first convolution region through this embodiment can reduce the calculation amount, while highlighting the signal features of interest, making the convolution result better reflect the key information of the signal, and contributing to the subsequent understanding, analysis, and processing of the signal, improving the accuracy and efficiency of subsequent analysis.

[0089] In an exemplary embodiment, in step S403, determining the first convolution region corresponding to the detection signal curve according to the first edge and the second edge may include the following steps:

[0090] Determine a rectangular region according to the first edge and the second edge; obtain the first convolution region corresponding to the detection signal curve based on the rectangular region.

[0091] Specifically, the first edge is the range defined by the first starting position and the first ending position in the detector channel direction, which can define the detector channel range interval. At the same time, the second edge is the range determined in the signal intensity direction according to the curve amplitude threshold, which can define the signal intensity amplitude interval. By combining the first edge and the second edge, a rectangular region can be formed.

[0092] Furthermore, the first convolution region corresponding to the detection signal curve can be determined based on this rectangular region. In some examples, the enclosed region of the detection signal curve and the horizontal axis in the coordinate system can be determined, and the first convolution region can be determined according to the intersection of this enclosed region and the rectangular region. For example, as Figure 5a and Figure 5b shown, the first convolution region 501 and the first convolution region 502 can be obtained in this way.

[0093] In this embodiment, by determining the first convolution region according to the rectangular region defined by the first edge and the second edge, irrelevant backgrounds or unimportant amplitude change parts that may exist in the signal can be quickly removed, enabling subsequent convolution operations to focus more on detecting the core information of the signal. Moreover, this method is simple and convenient, capable of quickly obtaining the first convolution region and improving the processing efficiency of cross-scattering correction.

[0094] In an exemplary embodiment, before determining the first convolution region corresponding to the detection signal curve according to the boundary position characteristics of the detection signal curve of the detector, the following steps may further be included:

[0095] Obtain a preset curve amplitude threshold; the curve amplitude threshold is less than the maximum amplitude of the detection signal curve; determine the curve points in the detection signal curve whose amplitudes match the curve amplitude threshold; obtain the boundary position characteristics of the detection signal curve according to the abscissas of the curve points.

[0096] In practical applications, the curve amplitude threshold can be set according to the processed detection signal curve and relevant prior indications. In some examples, the curve amplitude threshold is less than the maximum amplitude of any detection signal curve. Relevant personnel can determine a suitable curve amplitude threshold through multiple tests, which will not be elaborated here. By making the curve amplitude threshold less than the maximum amplitude of any detection signal curve, the threshold can be within a reasonable range. Specifically, if the curve amplitude threshold is greater than or equal to the maximum amplitude of the detection signal curve, it is difficult to effectively screen and analyze the signal. By ensuring that the threshold is less than the maximum amplitude, it helps to distinguish the roles of different amplitude parts in the detection signal.

[0097] Then, it is possible to determine the points in the detection signal curve whose amplitudes match the curve amplitude threshold, that is, the curve points. In some exemplary embodiments, the points in the detection signal curve whose signal amplitudes match the curve amplitude threshold can be used as the matching curves. For example, if the curve amplitude threshold is h, a straight line with a height of h and parallel to the horizontal axis of the coordinate system can be drawn, and the intersection points of this straight line and the detection signal curve can be determined as the matching curve points. Of course, in some other embodiments, the points (such as all points or some points) in the detection signal curve can also be compared with the preset curve amplitude threshold. When the amplitude of a certain point is equal to or close to (such as the difference between the two is less than the threshold) the curve amplitude threshold, it can be determined as the matching curve point.

[0098] Furthermore, the boundary position characteristics of the detection signal curve can be obtained according to the abscissas of the curve points. In some embodiments, if there are two matching curve points, these two points can be determined as the boundary position characteristics. If there are more than two matching curve points, the two curve points located at the beginning and end can be determined according to the abscissas, and then these two points can be determined as the boundary position characteristics.

[0099] In this embodiment, by setting a curve amplitude threshold smaller than the maximum amplitude of the detection signal curve and accordingly determining the corresponding curve points, it is possible to focus on the points in the detection signal whose amplitudes are below the curve amplitude threshold, thereby reducing the interference of the amplitude of the detection signal curve on the cross-scattering estimation and correction, and accurately identifying the boundary position features that can reflect the boundary position of the detection signal curve.

[0100] In an exemplary embodiment, the detection signal curve may include sub-signal curves corresponding to at least two tissue structures of the scanned object. For example, in Figure 6 the schematic diagram of the shoulder structure shown, it may include sub-signal curves corresponding to the neck and the two arms on both sides respectively; the multiple sub-signal curves are separated from each other. Correspondingly, in step S102, according to the boundary position features of the detection signal curve of the detector, determining the first convolution region corresponding to the detection signal curve may include:

[0101] For each sub-signal curve, according to the boundary position features of the sub-signal curve, determining the first convolution region corresponding to the sub-signal curve.

[0102] In practical applications, the detection signal curve may include sub-signal curves corresponding to at least two tissue structures of the scanned object. For example, Figure 6 shows a schematic diagram of a shoulder structure, where the central structure is the neck, the structures on the left and right sides are the raised arms, and the lower arc-shaped structure is the headrest. Combining Figure 2a with the cross-scattering scenario, when the tube 1 and the tube 2 are turned on simultaneously, as Figure 7a shown, the detection signal curve corresponding to the detection signal received by the detector 1 can obtain three sub-signal curves 701 corresponding to the arm, the neck, and the arm in sequence.

[0103] Furthermore, for each of the three sub-signal curves, the boundary position features of the sub-signal curve can be determined, and the first convolution region corresponding to the sub-signal curve can be determined according to the boundary position features. Among them, the determination method of the boundary position features of the word signal curve and the determination method of the first convolution region can refer to the introduction of determining the boundary position features and the first convolution region of the detection signal curve in the foregoing embodiments, which will not be elaborated here. For example, as Figure 7a shown, the first convolution regions of the three sub-signal curves can be obtained, that is, the rectangular regions enclosed by the dotted lines and the X-axis in the figure, and the scattered distribution signal shown in 702 can be obtained through the convolution operation.

[0104] On the other hand, combining Figure 2b with the cross-scattering scenario, when the tube 1 and the tube 2 are turned on simultaneously, the detection signal curve received by the detector 2 is as Figure 7bAs shown by 703 therein, at this time, since the detector 2 is in the same direction as the arm, neck, and tube, a detection signal curve is obtained at this time. As described above, after determining the first convolution region and performing the convolution operation, the scattered distribution signal shown by 704 can be obtained.

[0105] In this embodiment, by respectively determining the boundary position characteristics of each sub-signal curve and determining the first convolution region corresponding to the sub-signal curve according to the boundary position characteristics of each sub-signal curve, targeted analysis can be performed on the detection signal curve corresponding to each section of soft tissue, improving the recognition accuracy of the first convolution region.

[0106] In the convolution operation, a preset convolution kernel is involved. The preset convolution kernel can be determined by measurement or simulation. In an exemplary embodiment, the preset convolution kernel is determined through the following steps:

[0107] S1. Obtain the primary ray signal and cross-scattered signal collected by the target detector for the phantom; the target detector is any one of multiple detectors.

[0108] In specific implementation, the phantom can be scanned by a multi-source CT imaging system. Among them, the phantom can be any phantom with known phantom scattering parameters. In some examples, the phantom can be a water film. When scanning, the water film with known size and relevant phantom scattering parameters can be placed at the rotation center, or of course, it can also be placed at other positions. During the scanning process, the target detector in the multi-source CT imaging system can collect signals to obtain the primary ray signal and cross-scattered signal collected for the preset phantom.

[0109] Taking a dual-source CT imaging system as an example, tube 1 is paired with detector 1, and tube 2 is paired with detector 2. On the one hand, tube 1 can be turned on and the signal of detector 2 can be measured. At this time, since tube 2 is not turned on, the only signal received by detector 2 is the cross-scattered signal; on the other hand, tube 2 can be turned on and tube 1 can be turned off to measure the signal of detector 2. At this time, only the primary ray signal is received.

[0110] S2. According to the boundary position characteristics of the primary ray signal curve corresponding to the primary ray signal, determine the second convolution region corresponding to the primary ray signal curve, and, according to the cross-scattered signal and the known phantom scattering parameters of the phantom, determine the scattering estimation information corresponding to the target detector.

[0111] In this step, on the one hand, the main ray signal curve corresponding to the main ray signal can be determined, and then boundary analysis can be performed on the main ray signal curve to obtain the boundary position characteristics of the main ray signal curve. Then, based on the boundary position characteristics, the convolution region of the main ray signal curve can be determined. For the sake of distinction, this convolution region is called the second convolution region. On the other hand, since the scatter parameter of the phantom is known, the scatter estimation information corresponding to the target detector can be calculated based on the cross-scatter signal obtained by the target detector and the known scatter parameter of the phantom.

[0112] S3. Perform deconvolution processing on the scatter estimation information corresponding to the target detector according to the second convolution region to obtain a preset convolution kernel.

[0113] Furthermore, deconvolution processing can be performed on the scatter estimation information corresponding to the target detector according to the second convolution region. Here, deconvolution processing can be understood as the inverse operation of convolution processing. Thus, based on the calculation result obtained by the deconvolution processing, the preset convolution kernel for subsequent convolution processing can be determined.

[0114] In this embodiment, by scanning a preset phantom with known scatter parameters, the second convolution region and scatter estimation information are determined based on the obtained main ray signal and cross-scatter signal, and then deconvolution processing is performed, so that a preset convolution kernel that can accurately describe the actual scatter process can be obtained, improving the accuracy of subsequent cross-scatter correction.

[0115] In one embodiment, step S2 may include the following processing:

[0116] S21. Determine a second starting position and a second ending position in the detector channel direction of the main ray signal curve according to the boundary position characteristics of the main ray signal curve corresponding to the main ray signal, and determine a third edge according to the second starting position and the second ending position.

[0117] In practical applications, the main ray signal curve can be recorded in a two-dimensional coordinate system; the boundary position characteristics of the main ray signal curve may include the curve starting position and / or the curve ending position of the main ray signal curve, or may be other positions that are a preset distance away from the curve starting position and / or the curve ending position. In this step, the starting position and the ending position can be determined in the detector channel direction according to the boundary position characteristics of the main ray signal curve. For the sake of distinction, the starting position and the ending position are called the second starting position and the second ending position.

[0118] Then, the corresponding region edge is determined according to the second starting position and the second ending position. For the sake of distinction, the region edge determined according to the second starting position and the second ending position is called the third edge, and the third edge can limit the span of the second convolution region in the detector channel direction.

[0119] S22. Determine a fourth edge with a length less than a preset curve amplitude threshold in the signal intensity direction of the main ray signal curve.

[0120] In this step, the edge of the main ray signal curve in the signal intensity direction can be determined according to the preset curve amplitude threshold, and the curve amplitude threshold can be the same as the curve amplitude threshold used when determining the first convolution region. For the convenience of distinction, the edge determined in this step is called the fourth edge, and the length of the fourth edge is less than the curve amplitude threshold. For example, in the signal intensity direction, the distance corresponding to the curve amplitude threshold can be determined as the fourth edge.

[0121] S23. Determine a second convolution region corresponding to the main ray signal curve according to the third edge and the fourth edge.

[0122] After obtaining the third edge and the fourth edge, since the third edge defines the span of the second convolution region in the detector channel direction, and the fourth edge defines the amplitude of the second convolution region in the signal intensity direction, the second convolution region can be determined according to the third edge and the fourth edge.

[0123] In this embodiment, by determining the second convolution region in the same way as the convolution region determination method, on the one hand, it can reduce the interference of the amplitude of the signal curve on the cross-scattering correction result and reduce the calculation resource overhead. On the other hand, it can improve the matching degree between the preset convolution kernel obtained subsequently and the actual correction process with the first convolution region, and improve the accuracy of the obtained scattering estimation information.

[0124] It can be understood that the process of determining the second convolution region can be the same as the process of determining the first convolution region. The specific processing can refer to the foregoing embodiments and will not be elaborated here.

[0125] In an exemplary embodiment, the first convolution region includes the signal intensities of multiple detector channels in the detector. It can be understood that through the foregoing operation of determining the first convolution region, the signal intensity of each detector channel can be controlled within a certain amplitude range. In step S102, according to the convolution result of the first convolution region and the preset convolution kernel, determining the scattering estimation information corresponding to the detector may include the following steps:

[0126] For each detector channel in the first convolution region, obtain the signal intensity deviation of the detector channel according to the multiplication result of the signal intensity of the detector channel and the preset convolution kernel; determine the scattering estimation information corresponding to the detector according to the summation result of the signal intensity deviations of multiple detector channels.

[0127] In the implementation application, for the first convolution region, convolution operations can be performed along the detector channel direction. Specifically, for each detector channel in the first convolution region, the signal intensity of the detector channel can be multiplied by a preset convolution kernel, and the resulting multiplication result can be used as the signal intensity deviation of the detector channel. Furthermore, the signal intensity deviations of multiple detector channels can be summed, and the sum result can be used as the scattering estimation information corresponding to the detectors of the multiple detector channels.

[0128] In some examples, the scattering estimation information can be determined in the following manner :

[0129]

[0130] where is the preset convolution kernel, which can be a function related to the rotation angle and is the signal intensity corresponding to the detector channel that needs to be corrected at the corresponding rotation angle .

[0131] In this embodiment, on the one hand, the signal intensity deviation can be calculated independently for each detector channel, which can accurately capture the signal differences between different channels. Since the influence of scattering on each detector channel is random and individual, analyzing each channel separately can finely locate the signal anomalies of each channel and provide more accurate basic data for subsequent correction. On the other hand, by multiplying the signal intensity by the preset convolution kernel, the signal changes caused by scattering are converted into specific numerical deviations, making the originally abstract scattering influence measurable and analyzable, providing a clear numerical basis for subsequent image correction, and ensuring that the correction operation is targeted and conforms to the actual signal change law.

[0132] In an exemplary embodiment, in step S103, according to the scattering estimation information corresponding to each detector, determining the cross-scattering correction result of each detector may include the following steps:

[0133] Obtain the signal intensity of the detection signal collected by each detector; for each detector, adjust the signal intensity according to the scattering estimation information corresponding to the detector to obtain the adjusted signal intensity; based on the adjusted signal intensity of each detector, obtain the cross-scattering correction result of each detector.

[0134] In specific implementation, the signal intensity of the detection signal collected by each detector can be determined, and this signal intensity can be the signal intensity of all detection signals actually received by the detector (that is, including the primary ray signal and the cross-scattering signal).

[0135] Since the scattering estimation information corresponding to each detector has been obtained, for each detector, the signal intensity of the detector can be adjusted using the scattering estimation information. In some embodiments, the scattering estimation information can be superimposed on the signal intensity to obtain the adjusted signal intensity. For example, the signal superposition can be performed in the following manner:

[0136]

[0137] where is the adjusted signal intensity, is the signal intensity of the detection signal collected by the detector before adjustment, is the scattering estimation information.

[0138] Furthermore, the adjusted signal intensity of each detector can be used as the cross-scattering correction result of each detector. Of course, the adjusted signal intensity can also be converted into corrected projection data, and the corrected projection data can be used as the cross-scattering correction result.

[0139] In this embodiment, based on the scattering estimation information of each detector obtained previously, the signal intensity of each detector can be differentially and specifically corrected. Since the degree of influence of cross-scattering on different detectors is different, this adjustment method can specifically offset the interference of scattering on the signal, eliminate the interference added to the signal due to scattering, make the signal closer to the true physical attenuation value, and lay a solid data foundation for high-quality imaging.

[0140] To enable those skilled in the art to better understand the above steps, the following gives an exemplary illustration of the embodiments of the present application through an example. However, it should be understood that the embodiments of the present application are not limited thereto.

[0141] As Figure 8 shown, the following steps are included in this example:

[0142] S801, obtain the primary ray signal and the cross-scattering signal collected by the target detector in the dual-source or multi-source CT imaging system for the water film.

[0143] S802, determine the second convolution region corresponding to the primary ray signal curve according to the boundary position characteristics of the primary ray signal curve corresponding to the primary ray signal, and determine the scattering estimation information corresponding to the target detector according to the cross-scattering signal and the known phantom scattering parameters of the water film.

[0144] S803, perform deconvolution processing on the scattering estimation information corresponding to the target detector according to the second convolution region to obtain a preset convolution kernel.

[0145] S804. Obtain the detection signal curves of multiple detectors respectively. The detection signal curves are determined based on the detection signals collected by the detectors in a dual-source or multi-source CT imaging system for a scanned object.

[0146] S805. For each detector, determine a first starting position and a first ending position in the detector channel direction according to the boundary position characteristics of the detection signal curve of the detector, determine a first edge according to the first starting position and the first ending position, determine a second edge with a length less than a preset curve amplitude threshold in the signal intensity direction, determine a first convolution region corresponding to the detection signal curve according to the first edge and the second edge, and determine scattering estimation information according to the convolution result of the first convolution region and a preset convolution kernel.

[0147] S806. Determine the cross-scattering correction result of each detector according to the scattering estimation information corresponding to each detector.

[0148] It should be understood that although the steps in the flowcharts involved in the above-described embodiments are shown in sequence according to the arrows, these steps do not necessarily have to be executed in the order indicated by the arrows. Unless there is a clear indication in this article, there is no strict order limit for the execution of these steps, 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 do not necessarily have to be executed at the same time, but can be executed at different times. The execution order of these steps or stages does not necessarily have to be sequential, but can be executed alternately or in turn with at least a part of other steps or steps or stages in other steps.

[0149] Based on the same inventive concept, an embodiment of the present application further provides a cross-scattering correction device for implementing the above-mentioned cross-scattering 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 embodiments of the following cross-scattering correction device can refer to the limitations on the cross-scattering correction method in the above text, and will not be repeated here.

[0150] In an exemplary embodiment, as Figure 9 shown, a cross-scattering 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. The device includes:

[0151] A curve acquisition module 901, configured to obtain the detection signal curves of multiple detectors respectively. The detection signal curves are determined based on the detection signals collected by each detector for a scanned object;

[0152] The convolution module 902 is configured to, for each of the detectors, determine a first convolution region corresponding to the detection signal curve according to the boundary position feature of the detection signal curve of the detector, and determine the scatter estimation information corresponding to the detector according to the convolution result of the first convolution region and a preset convolution kernel;

[0153] The correction result acquisition module 903 is configured to determine the cross-scatter correction result of each detector according to the scatter estimation information corresponding to each detector.

[0154] In one embodiment, the convolution module 902 is configured to:

[0155] According to the boundary position feature of the detection signal curve of the detector, determine a first start position and a first end position in the detector channel direction of the detection signal curve, and determine a first edge according to the first start position and the first end position;

[0156] Determine a second edge with a length less than a preset curve amplitude threshold in the signal intensity direction of the detection signal curve;

[0157] Determine a first convolution region corresponding to the detection signal curve according to the first edge and the second edge.

[0158] In one embodiment, the convolution module 902 is configured to:

[0159] Determine a rectangular region according to the first edge and the second edge;

[0160] Obtain a first convolution region corresponding to the detection signal curve according to the rectangular region.

[0161] In one embodiment, the convolution module 902 is further configured to:

[0162] Obtain a preset curve amplitude threshold; the curve amplitude threshold is less than the maximum amplitude of the detection signal curve;

[0163] Determine curve points in the detection signal curve whose amplitudes match the curve amplitude threshold;

[0164] Obtain the boundary position feature of the detection signal curve according to the abscissa of the curve points.

[0165] In one embodiment, the detection signal curve includes sub-signal curves corresponding to at least two tissue structures of the scanned object, and the multiple sub-signal curves are separated from each other;

[0166] The step of determining a first convolution region corresponding to the detection signal curve according to the boundary position feature of the detection signal curve of the detector includes:

[0167] For each of the sub-signal curves, a first convolution region corresponding to the sub-signal curve is determined according to a boundary position feature of the sub-signal curve.

[0168] In one embodiment, the apparatus further includes a convolution kernel determination module, wherein the convolution kernel determination module is configured to:

[0169] Acquire the main ray signal and the cross-scattering signal collected by the target detector for the phantom; the target detector is any one of the multiple detectors;

[0170] determining a second convolution region corresponding to the main ray signal curve based on boundary position features of a main ray signal curve corresponding to the main ray signal, and determining scattering estimation information corresponding to the target detector based on the cross scattering signal and known phantom scattering parameters of the phantom;

[0171] Deconvolution processing is performed on the scattering estimation information corresponding to the target detector according to the second convolution region to obtain a preset convolution kernel.

[0172] In one embodiment, the convolution kernel determination module is used to:

[0173] determining a second starting position and a second ending position in the detector channel direction of the main ray signal curve according to boundary position characteristics of the main ray signal curve corresponding to the main ray signal, and determining a third edge according to the second starting position and the second ending position;

[0174] Determining a fourth edge having a length less than a preset curve amplitude threshold in the signal intensity direction of the main ray signal curve;

[0175] A second convolution region corresponding to the main ray signal curve is determined according to the third edge and the fourth edge.

[0176] In one embodiment, the first convolution region includes the signal intensity of each of a plurality of detector channels in the detector;

[0177] The convolution module is used to:

[0178] For each of the detector channels in the first convolution region, obtaining a signal intensity deviation of the detector channel according to a multiplication result of the signal intensity of the detector channel and the preset convolution kernel;

[0179] The scattering estimation information corresponding to the detector is determined according to a summation result of the signal intensity deviations of the plurality of detector channels.

[0180] In one embodiment, the correction result acquisition module is used to:

[0181] Obtain the signal intensity of the detection signal collected by each of the detectors;

[0182] For each of the detectors, adjust the signal intensity according to the scattering estimation information corresponding to the detector to obtain the adjusted signal intensity;

[0183] Based on the adjusted signal intensity of each of the detectors, obtain the cross-scattering correction result of each of the detectors.

[0184] 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 independent of it, or 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.

[0185] 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 10 shown. 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 detectors. 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, a cross-scattering correction method is implemented.

[0186] 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 11As 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 computer programs. The internal memory provides an environment for the operation of the operating system and computer programs 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, and the wireless manner can be implemented through WIFI, a mobile cellular network, 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 covered on the display screen, or a button, a trackball, or a touchpad set on the computer device housing, or an external keyboard, touchpad, or mouse, etc.

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

[0188] 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.

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

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

[0191] 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.

[0192] 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 processing units, graphics processing units, digital signal processors, programmable logic devices, data processing logics based on quantum computing, artificial intelligence (AI) processors, etc., without limitation.

[0193] The technical features of the above embodiments can be combined arbitrarily. For the sake of concise 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 as the scope recorded in this application.

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

Claims

1. A cross-scattering correction method applied to a dual-source or multi-source CT imaging system, the dual-source or multi-source CT imaging system including at least two detectors, characterized in that, The method includes: Obtaining the detection signal curves of each of the multiple detectors, where the detection signal curves are determined based on the detection signals collected by each detector for the scanned object; For each detector, determining a first convolution region corresponding to the detection signal curve according to the boundary position characteristics of the detection signal curve of the detector, and determining the scatter estimation information corresponding to the detector according to the convolution result of the first convolution region and a preset convolution kernel; Determining the cross-scatter correction result of each detector according to the scatter estimation information corresponding to each detector.

2. The method according to claim 1, wherein The determining the first convolution region corresponding to the detection signal curve according to the boundary position characteristics of the detection signal curve of the detector includes: According to the boundary position characteristics of the detection signal curve of the detector, determining a first start position and a first end position in the detector channel direction of the detection signal curve, and determining a first edge according to the first start position and the first end position; Determining a second edge with a length less than a preset curve amplitude threshold in the signal intensity direction of the detection signal curve; Determining the first convolution region corresponding to the detection signal curve according to the first edge and the second edge.

3. The method according to claim 2, wherein The determining the first convolution region corresponding to the detection signal curve according to the first edge and the second edge includes: Determining a rectangular region according to the first edge and the second edge; Obtaining the first convolution region corresponding to the detection signal curve according to the rectangular region.

4. The method according to claim 1, characterized in that, Before the determining the first convolution region corresponding to the detection signal curve according to the boundary position characteristics of the detection signal curve of the detector, it further includes: Obtaining a preset curve amplitude threshold; the curve amplitude threshold is less than the maximum amplitude of the detection signal curve; Determining curve points in the detection signal curve whose amplitudes match the curve amplitude threshold; Obtaining the boundary position characteristics of the detection signal curve according to the abscissa of the curve points.

5. The method according to claim 1, characterized in that The preset convolution kernel is determined through the following steps: Obtaining the primary ray signal and the cross-scatter signal collected by a target detector for a phantom; the target detector is any one of the multiple detectors; Determining a second convolution region corresponding to the primary ray signal curve according to the boundary position characteristics of the primary ray signal curve corresponding to the primary ray signal, and determining the scatter estimation information corresponding to the target detector according to the cross-scatter signal and the known phantom scatter parameters of the phantom; Performing deconvolution processing on the scatter estimation information corresponding to the target detector according to the second convolution region to obtain a preset convolution kernel.

6. The method according to claim 5, characterized in that, The determining the second convolution region corresponding to the primary ray signal curve according to the boundary position characteristics of the primary ray signal curve corresponding to the primary ray signal includes: According to the boundary position characteristics of the primary ray signal curve corresponding to the primary ray signal, determining a second start position and a second end position in the detector channel direction of the primary ray signal curve, and determining a third edge according to the second start position and the second end position; Determine a fourth edge with a length less than a preset curve amplitude threshold in the signal intensity direction of the main ray signal curve; Determine a second convolution region corresponding to the main ray signal curve according to the third edge and the fourth edge.

7. The method according to claim 1, characterized in that, The first convolution region includes the signal intensities of multiple detector channels in the detector; The determining the scatter estimation information corresponding to the detector according to the convolution result of the first convolution region and a preset convolution kernel includes: For each detector channel in the first convolution region, obtain the signal intensity deviation of the detector channel according to the multiplication result of the signal intensity of the detector channel and the preset convolution kernel; Determine the scatter estimation information corresponding to the detector according to the summation result of the signal intensity deviations of the multiple detector channels.

8. The method according to any one of claims 1 to 7, characterized in that, The determining the cross-scatter correction result of each detector according to the scatter estimation information corresponding to each detector includes: Obtain the signal intensity of the detection signal collected by each detector; For each detector, adjust the signal intensity according to the scatter estimation information corresponding to the detector to obtain an adjusted signal intensity; Based on the adjusted signal intensity of each detector, obtain the cross-scatter correction result of each detector.

9. A cross-scattering correction device 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 is characterized in that, The device includes: A curve acquisition module, configured to acquire detection signal curves of multiple detectors respectively, where the detection signal curves are determined based on detection signals collected by each detector for a scanned object; A convolution module, configured to, for each detector, determine a first convolution region corresponding to the detection signal curve according to the boundary position feature of the detection signal curve of the detector, and determine the scatter estimation information corresponding to the detector according to the convolution result of the first convolution region and a preset convolution kernel; A correction result acquisition module, configured to determine the cross-scatter correction result of each detector according to the scatter estimation information corresponding to each detector.

10. 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 8 are implemented.

Images