Scattering correction methods, devices, electronic equipment and storage media

Abstract

This invention provides a scattering correction method, apparatus, electronic device, and storage medium. The method includes: obtaining at least two data subsets based on energy spectrum information of PET data; and determining preliminary scattering correction data based on the at least two data subsets and a preset scattering estimation model. This invention can determine the preliminary scattering correction data based on at least two data subsets and a preset scattering estimation model, reducing the iterative process required for model estimation of scattering correction, improving the speed of acquiring preliminary scattering correction data, and thus improving scattering correction efficiency.

Description

Technical Field

[0001] This invention relates to the field of scattering correction technology, and specifically to a scattering correction method, apparatus, electronic device, and storage medium. Background Technology

[0002] With the continuous development of medical imaging technology, positron emission tomography (PET) has emerged to better examine the human body. PET is a relatively advanced clinical imaging technique in the field of medical imaging. However, during PET scanning, due to the Compton effect, two gamma photons, while losing energy, deviate from their original direction of motion. Although the two photons detected by the detector originate from the same electron-positron annihilation event, at least one photon undergoes one or more scattering events with the medium; such events are called scattering coincidence events. Scattering coincidence events can cause severe image noise, poor contrast, and inaccurate quantitative analysis, seriously affecting image quality. Therefore, scattering correction is necessary in modern PET imaging.

[0003] Currently, most scattering correction algorithms in the industry employ model estimation-based methods, such as single-scattering simulation, two-scattering simulation, and Monte Carlo simulation. While model estimation-based methods are relatively accurate, they require precise distributions of radioactivity and attenuation for model calculation. Due to the lack of accurate radioactivity distribution information in the initial conditions, accurate estimates can only be obtained through multiple iterations from uncorrected PET images, resulting in a relatively time-consuming process. Furthermore, as the axial length of the PET system increases, the computational load required for model estimation algorithms increases, especially for two-scattering and Monte Carlo simulation algorithms, where multiple iterations necessitate even more computation time. Summary of the Invention

[0004] In view of this, it is necessary to provide a scattering correction method, apparatus, electronic device and storage medium to solve the technical problems of long scattering correction time and low efficiency in the prior art.

[0005] On one hand, the present invention provides a scattering correction method, comprising:

[0006] Based on the energy spectrum information of PET data, at least two data subsets are obtained;

[0007] Preliminary scattering correction data are determined based on the at least two subsets of data and a preset scattering estimation model.

[0008] In some possible implementations, the energy spectrum information based on PET data yields at least two data subsets, including:

[0009] Determine the upper limit of energy discrimination, the lower limit of energy discrimination, at least one first energy discrimination threshold, and at least one second energy discrimination threshold;

[0010] The at least two data subsets are obtained from the energy spectrum information of the PET data based on the energy discrimination upper limit, energy discrimination lower limit, at least one first energy discrimination threshold, and at least one second energy discrimination threshold.

[0011] In some possible implementations, the at least two data subsets include a first data subset and a second data subset; determining the preliminary scattering correction data based on the at least two data subsets includes:

[0012] Determine the number of events that match the events in the first data subset;

[0013] The preliminary scattering correction data is determined based on the number of events counted and the preset scattering estimation model.

[0014] In some possible implementations, after determining the preliminary scattering correction data, the method further includes:

[0015] Based on a preset scattering estimation model, estimated scattering correction data is obtained, and the preliminary scattering correction data is optimized based on the estimated scattering correction data to generate optimized scattering correction data.

[0016] In some possible implementations, after determining the preliminary scattering correction data, the method further includes:

[0017] Based on the preliminary scattering correction data and the preset image reconstruction algorithm, the radioactivity distribution map is obtained;

[0018] Obtain the attenuation information distribution map, and optimize the preliminary scattering correction data based on the attenuation information distribution map and the radioactivity distribution map to generate optimized scattering correction data.

[0019] In some possible implementations, after determining the preliminary scattering correction data, the method further includes:

[0020] Based on the preliminary scattering correction data and the preset image reconstruction algorithm, a radioactivity distribution map is obtained;

[0021] A deep learning network is constructed, and the radioactivity distribution map is optimized based on the deep learning network to generate an optimized radioactivity distribution map;

[0022] Obtain the attenuation information distribution map, and optimize the preliminary scattering correction data based on the attenuation information distribution map and the optimized radioactivity distribution map to generate optimized scattering correction data.

[0023] In some possible implementations, optimizing the preliminary scattering correction data based on the attenuation information distribution map and the radioactivity distribution map to generate optimized scattering correction data includes:

[0024] The preliminary scattering correction data is optimized based on the attenuation information distribution map and the radioactivity distribution map to generate transition scattering correction data;

[0025] Determine whether the number of optimization attempts for the initial scattering correction data exceeds the threshold number of optimization attempts;

[0026] If the number of optimization attempts is greater than the threshold number of optimization attempts, then the transition scattering correction data is the optimized scattering correction data;

[0027] If the number of optimization attempts is less than or equal to the threshold number of optimization attempts, then the transition scattering correction data is optimized based on the attenuation information distribution map and the radioactivity distribution map.

[0028] On the other hand, the present invention also provides a scattering correction device, comprising:

[0029] The data subset acquisition unit is used to obtain at least two data subsets based on the energy spectrum information of PET data;

[0030] A scattering correction unit is used to determine preliminary scattering correction data based on the at least two subsets of data.

[0031] On the other hand, the present invention also provides an electronic device, including a memory and a processor, wherein,

[0032] The memory is used to store programs;

[0033] The processor, coupled to the memory, is used to execute the program stored in the memory to implement the steps of the scattering correction method described in any of the above implementations.

[0034] On the other hand, the present invention also provides a computer-readable storage medium for storing a computer-readable program or instructions, which, when executed by a processor, can implement the steps of the scattering correction method described in any of the above implementations.

[0035] The beneficial effects of the above embodiments are as follows: The scattering correction method provided by the present invention obtains at least two data subsets based on the energy spectrum information of PET data. The preliminary scattering correction data can be determined according to the at least two data subsets and the preset scattering estimation model, which reduces the iterative process required for model estimation of scattering correction, improves the acquisition speed of preliminary scattering correction data, and thus improves the scattering correction efficiency. Attached Figure Description

[0036] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0037] Figure 1 A schematic diagram of an embodiment of a scattering coincidence event is provided for the present invention;

[0038] Figure 2 A schematic flowchart of an embodiment of the scattering correction method provided by the present invention;

[0039] Figure 3 A schematic diagram of an embodiment of the one-dimensional energy spectrum information provided by the present invention;

[0040] Figure 4 A schematic diagram of an embodiment of the two-dimensional energy spectrum information provided by the present invention;

[0041] Figure 5 For the present invention Figure 2 A schematic diagram of an embodiment of S201;

[0042] Figure 6 For the present invention Figure 2 A schematic diagram of an embodiment of S202;

[0043] Figure 7 A schematic flowchart of an embodiment of the present invention for optimizing preliminary scattering correction data;

[0044] Figure 8 This is a schematic flowchart of another embodiment of the present invention for optimizing preliminary scattering correction data.

[0045] Figure 9 For the present invention Figure 7 A schematic diagram of an embodiment of S702;

[0046] Figure 10 A schematic diagram of an embodiment of the scattering correction device provided by the present invention;

[0047] Figure 11 A schematic diagram of an embodiment of the electronic device provided by the present invention. Detailed Implementation

[0048] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of them. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative effort are within the scope of protection of the present invention.

[0049] In the description of the embodiments of the present invention, unless otherwise stated, "multiple" means two or more. "And / or" describes the relationship between related objects, indicating that there can be three relationships. For example, A and / or B can represent three situations: A exists alone, A and B exist simultaneously, and B exists alone.

[0050] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of the invention. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.

[0051] This invention provides a scattering correction method, apparatus, electronic device, and storage medium, which are described below.

[0052] Before demonstrating the embodiments, the scattering coincidence event in PET will be described in detail. For example... Figure 1 As shown, the PET scanner includes a detector ring 11 containing multiple detectors. An annihilation event occurs at point H, generating a pair of gamma photons, L and R. Photon L travels along path HA to detector A, where it is recorded. However, photon R encounters a scattering point S during its flight, undergoing Compton scattering at S, which changes its flight direction, causing it to reach detector B and be recorded by detector B. Therefore, if the scattering coincidence event is not considered, the PET scanning device would assume that the annihilation event occurred on the response lines of detectors A and B. If image reconstruction is performed directly using this data, the resulting image would be inaccurate.

[0053] To address the above problems, embodiments of the present invention provide a scattering correction method, such as... Figure 2 As shown, the scattering correction method includes:

[0054] S201. Based on the energy spectrum information of PET data, obtain at least two data subsets;

[0055] S202. Determine preliminary scattering correction data based on at least two data subsets and a pre-defined scattering estimation model.

[0056] It should be noted that in some embodiments of the present invention, the PET data is in the form of list-mode data, which records information such as the energy deposition location, energy, and detection event of photons in each annihilation event.

[0057] In step S202, the preliminary scattering correction data is generally presented as a scattering correction chord diagram.

[0058] In some embodiments of the present invention, detector detection energy spectrum information can be generated based on the energy information in the list pattern data, and the energy spectrum information conforming to the single end of the detector is as follows: Figure 3 The one-dimensional energy spectrum information shown corresponds to the energy spectrum information at both ends of the detector as follows: Figure 4 The diagram shows two-dimensional energy spectrum information. In one-dimensional energy spectrum information, the horizontal axis represents the gamma-ray energy value received by a single-ended detector, and the vertical axis represents the number of coincidence events. In two-dimensional energy spectrum information, the horizontal and vertical axes represent the gamma-ray energy values ​​received by two-ended detectors, where the horizontal axis represents the energy value received by the gamma-ray hitting one detector, and the vertical axis represents the energy value received by the gamma-ray hitting the other detector. Points represent coincidence events.

[0059] Compared with the prior art, the scattering correction method provided in this embodiment of the invention obtains at least two data subsets based on the energy spectrum information of PET data. The preliminary scattering correction data can be determined based on the at least two data subsets and the preset scattering estimation model. This reduces the iterative process required for model estimation of scattering correction, improves the acquisition speed of preliminary scattering correction data, and thus improves the scattering correction efficiency.

[0060] In some embodiments of the present invention, such as Figure 5 As shown, step S201 includes:

[0061] S501. Determine the upper limit of energy discrimination, the lower limit of energy discrimination, at least one first energy discrimination threshold, and at least one second energy discrimination threshold;

[0062] S502. Obtain at least two data subsets from the energy spectrum information of PET data based on the upper limit of energy discrimination, the lower limit of energy discrimination, at least one first energy discrimination threshold, and at least one second energy discrimination threshold.

[0063] By determining an upper limit for energy discrimination, a lower limit for energy discrimination, at least one first energy discrimination threshold, and at least one second energy discrimination threshold, the present invention can divide at least two data subsets from the energy spectrum information, thereby improving the speed of data subset division and further improving the scattering correction efficiency.

[0064] In some embodiments of the present invention, the upper and lower limits of energy discrimination in step S501 should be determined based on the original conformity window range in the PET data. Specifically, the upper and lower boundaries of the original conformity window range are the upper and lower limits of energy discrimination, respectively. Figure 3 and Figure 4 As shown, E LLD E is the lower limit for energy discrimination. ULD The upper limit for energy identification.

[0065] In some embodiments of the present invention, the first energy discrimination threshold and the second energy discrimination threshold in step S501 should be determined based on the energy of the PET detection of a true coincidence event and the energy resolution of the PET system.

[0066] In a specific embodiment of the present invention, the energy of a true coincidence event detected by PET is 511 keV, and the energy resolution of the PET system is 12%. Since the energy resolution follows a Gaussian distribution, its standard deviation δ is 26 keV. Therefore, the first energy discrimination threshold is 511 - 2*δ = 459 keV, and the second energy discrimination threshold is 511 + 2*δ = 563 keV. Figure 3 and Figure 4 As shown, E A The first energy discrimination threshold; E B This is the second energy discrimination threshold.

[0067] And such as Figure 3 and Figure 4 As shown, there are at least two data subsets, namely the first data subset and the second data subset, where the first data subset is E. LLD <E<E A ∪E B <E<E ULD The second data subset is E A <E<E B .

[0068] That is: the first subset of data includes a higher coincidence window E LLD <E<E A and a lower conformity window E B <E<E ULD By using this higher and lower coincidence energy windows, the coincidence scattering events within the original coincidence energy window are estimated together.

[0069] In some embodiments of the present invention, such as Figure 6 As shown, step S202 includes:

[0070] S601. Determine the number of events that match the events in the first data subset;

[0071] S602. Determine preliminary scattering correction data based on the number of events and the preset scattering estimation model.

[0072] In some embodiments of the present invention, the event count in step S601 includes E. LLD <E<E A The number of events and E in B <E<E ULD The number of events counted.

[0073] In some embodiments of the present invention, the preset scattering estimation model in step S602 includes a one-dimensional scattering estimation model and a two-dimensional scattering estimation model. Specifically, the one-dimensional scattering estimation model is:

[0074]

[0075] The two-dimensional scattering estimation model is as follows:

[0076]

[0077] In the formula, SC represents the scattering coincidence event estimate within the second data subset. Count the number of events that match the event in the first data subset.

[0078] Using the one-dimensional or two-dimensional scattering estimation model in the above embodiments, the scattering coincidence event estimates in the second data subset can be obtained. By superimposing the scattering coincidence event estimates in the second data subset with the coincidence events in the first data subset, preliminary scattering correction data can be obtained.

[0079] To further reduce computational steps and improve the efficiency of scattering correction, in some embodiments of the present invention, the scattering estimation model can directly determine the preliminary scattering correction data. Specifically, the one-dimensional scattering estimation model is as follows:

[0080]

[0081] The two-dimensional scattering estimation model is as follows:

[0082]

[0083] In the formula, SC1 represents the preliminary scattering correction data.

[0084] Since the above embodiments use all coincidence events in the first data subset as estimates of scattering coincidence events to obtain preliminary scattering correction data, and since all coincidence events in the first data subset include not only scattering coincidence events but also some true coincidence events, in order to further improve the accuracy of scattering correction, in some embodiments of the present invention, after determining the preliminary scattering correction data, the scattering correction method further includes:

[0085] The initial scattering correction data is optimized to generate optimized scattering correction data.

[0086] In some embodiments of the present invention, the preliminary scattering correction data is optimized to generate optimized scattering correction data, including:

[0087] Based on a preset scattering estimation model, estimated scattering correction data is obtained, and the preliminary scattering correction data is optimized based on the estimated scattering correction data to generate optimized scattering correction data.

[0088] Optimizing the preliminary scattering correction data using the estimated scattering correction data obtained from the scattering estimation model can improve the accuracy of the scattering correction data.

[0089] In a specific embodiment of the present invention, an empirical correction model can be determined based on the estimated scattering correction data, and the preliminary scattering correction data can be optimized using the empirical correction model.

[0090] For example, preliminary scattering correction data and estimated scattering correction data are estimated using the scattering correction method and scattering estimation model in the embodiments of the present invention, respectively, and the difference between the preliminary scattering correction data and the estimated scattering correction data is determined. An empirical correction model for the preliminary scattering correction data is determined by means of smooth fitting, etc., and the preliminary scattering correction data is optimized by the empirical correction model to generate optimized scattering correction data.

[0091] In addition to the method for optimizing preliminary scattering correction data in the above embodiments, in other embodiments of the present invention, such as... Figure 7 As shown, the method for optimizing preliminary scattering correction data to generate optimized scattering correction data may further include:

[0092] S701. Based on the preliminary scattering correction data and the preset image reconstruction algorithm, obtain the radioactivity distribution map;

[0093] S702. Obtain the attenuation information distribution map, and optimize the preliminary scattering correction data based on the attenuation information distribution map and the radioactivity distribution map to generate optimized scattering correction data.

[0094] In some embodiments of the present invention, the image reconstruction algorithm preset in step S701 can be any one of the following: Filtered Back-projection (FBP) reconstruction algorithm and Ordered Subset Expectation Maximum (OSEM) reconstruction algorithm.

[0095] In some embodiments of the present invention, the attenuation information distribution map in step S702 may be pre-stored in a storage medium and retrieved directly from the storage medium when step S702 is performed, or the attenuation information distribution map may be generated and obtained in real time when step S702 is performed.

[0096] In one specific embodiment of the present invention, the attenuation information distribution map can be obtained by computed tomography (CT) or magnetic resonance (MR) scanning.

[0097] In some embodiments of the present invention, step S702 optimizes the preliminary scattering correction data based on the attenuation information distribution map and the radioactivity distribution map, specifically as follows:

[0098] The preliminary scattering correction data are optimized based on the attenuation information distribution map, the radioactivity distribution map, and the preset model estimation algorithm.

[0099] The preset model estimation algorithms include, but are not limited to, Single Scatter Simulation (SSS), Double Scatter Simulation (DSS), and Monte Carlo Simulation (MCS).

[0100] To improve the accuracy of the radioactivity distribution map obtained in step S701, in some embodiments of the present invention, step S701 may further be specifically: obtaining the radioactivity distribution map by combining image reconstruction algorithms with Time-of-Flight (TOF) information, thereby improving the accuracy of the radioactivity distribution map.

[0101] In some embodiments of the present invention, in order to reduce the inaccuracy of optimized scattering correction data due to inaccuracies in the radioactivity distribution map, in some embodiments of the present invention, such as Figure 8 As shown, the method for optimizing preliminary scattering correction data to generate optimized scattering correction data may further include:

[0102] S801. Based on the preliminary scattering correction data and the preset image reconstruction algorithm, obtain the radioactivity distribution map;

[0103] S802. Construct a deep learning network, optimize the radioactivity distribution map based on the deep learning network, and generate an optimized radioactivity distribution map.

[0104] S803. Obtain the attenuation information distribution map, and optimize the preliminary scattering correction data based on the attenuation information distribution map and the optimized radioactivity distribution map to generate optimized scattering correction data.

[0105] The embodiments of the present invention can further improve the reliability of scattering correction data by optimizing the radioactivity distribution map.

[0106] It should be understood that before optimizing the radioactivity distribution map based on a deep learning network, the deep learning network should be trained with training samples, and the trained deep learning network should be used to optimize the radioactivity distribution map.

[0107] In some embodiments of the present invention, such as Figure 9 As shown, step S702 further includes:

[0108] S901. Optimize the preliminary scattering correction data based on the attenuation information distribution map and the radioactivity distribution map to generate transition scattering correction data;

[0109] S902. Determine whether the number of optimizations for the initial scattering correction data is greater than the threshold number of optimizations.

[0110] S903. If the number of optimization attempts exceeds the threshold number of optimization attempts, then the transition scattering correction data is the optimized scattering correction data.

[0111] S904. If the number of optimizations is less than or equal to the threshold number of optimizations, then the transition scattering correction data is optimized based on the attenuation information distribution map and the radioactivity distribution map.

[0112] In some embodiments of the present invention, the number of threshold optimization iterations in step S902 is a positive integer greater than 1. The specific value of the number of threshold optimization iterations should be adjusted according to the requirements of accuracy and speed.

[0113] Specifically: when the accuracy requirement for optimizing scattering data is high, the threshold optimization number can be set to a larger number; when the speed requirement for optimizing scattering data is high, the threshold optimization number can be set to a smaller number. The specific value of the threshold optimization number is not limited here.

[0114] To better implement the scattering correction method in the embodiments of the present invention, based on the scattering correction method, correspondingly, as follows: Figure 10As shown, this embodiment of the invention also provides a scattering correction device, the scattering correction device 1000 comprising:

[0115] The data subset acquisition unit 1001 is used to obtain at least two data subsets based on the energy spectrum information of the list pattern data.

[0116] Scattering correction unit 1002 is used to determine preliminary scattering correction data based on at least two subsets of data.

[0117] The scattering correction device 1000 provided in the above embodiments can realize the technical solutions described in the above scattering correction method embodiments. The specific implementation principles of each module or unit can be found in the corresponding content in the above scattering correction method embodiments, and will not be repeated here.

[0118] like Figure 11 As shown, the present invention also provides an electronic device 1100. The electronic device 1100 includes a processor 1101, a memory 1102, and a display 1103. Figure 11 Only some components of the electronic device 1100 are shown, but it should be understood that it is not required to implement all the components shown, and more or fewer components may be implemented instead.

[0119] In some embodiments, memory 1102 may be an internal storage unit of electronic device 1100, such as a hard disk or memory of electronic device 1100. In other embodiments, memory 1102 may also be an external storage device of electronic device 1100, such as a plug-in hard disk, smart media card (SMC), secure digital (SD) card, flash card, etc. equipped on electronic device 1100.

[0120] Furthermore, the memory 1102 may include both internal storage units of the electronic device 1100 and external storage devices. The memory 1102 is used to store application software and various types of data installed on the electronic device 1100.

[0121] In some embodiments, processor 1101 may be a central processing unit (CPU), microprocessor or other data processing chip, used to run program code stored in memory 1102 or process data, such as the scattering correction method of the present invention.

[0122] In some embodiments, display 1103 may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, or an OLED (Organic Light-Emitting Diode) touchscreen. Display 1103 is used to display information from electronic device 1100 and to display a visual user interface. Components 1101-1103 of electronic device 1100 communicate with each other via a system bus.

[0123] In one embodiment, when processor 1101 executes the scattering correction program in memory 1102, the following steps may be performed:

[0124] Based on the energy spectrum information of PET data, at least two data subsets are obtained;

[0125] Preliminary scattering correction data are determined based on at least two subsets of data.

[0126] It should be understood that when the processor 1101 executes the scattering correction program in the memory 1102, in addition to the functions mentioned above, it can also perform other functions, as can be found in the description of the corresponding method embodiments above.

[0127] Furthermore, the embodiments of the present invention do not specifically limit the type of electronic device 1100 mentioned. Electronic device 1100 can be a mobile phone, tablet computer, personal digital assistant (PDA), wearable device, laptop computer, or other portable electronic device. Exemplary embodiments of portable electronic devices include, but are not limited to, portable electronic devices running iOS, Android, Microsoft, or other operating systems. The aforementioned portable electronic device can also be other portable electronic devices, such as a laptop computer with a touch-sensitive surface (e.g., a touch panel). It should also be understood that in some other embodiments of the present invention, electronic device 1100 may not be a portable electronic device, but rather a desktop computer with a touch-sensitive surface (e.g., a touch panel).

[0128] Accordingly, this application also provides a computer-readable storage medium for storing a computer-readable program or instruction. When the program or instruction is executed by a processor, it can implement the scattering correction method steps or functions provided in the above-described method embodiments.

[0129] The scattering correction method, apparatus, electronic device, and storage medium provided in this invention obtain at least two data subsets based on the energy spectrum information of list pattern data. Preliminary scattering correction data can be determined based on these at least two data subsets, reducing the iterative process required for model estimation of scattering correction, improving the acquisition speed of preliminary scattering correction data, and thus improving scattering correction efficiency.

[0130] Furthermore, by setting various methods to optimize the preliminary scattering correction data, the embodiments of the present invention generate optimized scattering correction data, which can improve the accuracy of the scattering correction data.

[0131] Those skilled in the art will understand that all or part of the processes of the methods described in the above embodiments can be implemented by a computer program instructing related hardware (such as a processor, controller, etc.), and the computer program can be stored in a computer-readable storage medium. The computer-readable storage medium may be a disk, optical disk, read-only memory, or random access memory, etc.

[0132] The scattering correction method, apparatus, electronic device, and storage medium provided by the present invention have been described in detail above. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A scattering correction method, characterized in that, include: Based on the energy spectrum information of PET data, at least two data subsets are obtained; Preliminary scattering correction data are determined based on the at least two subsets of data and the preset scattering estimation model; The energy spectrum information based on PET data yields at least two data subsets, including: Determine the upper limit of energy discrimination, the lower limit of energy discrimination, at least one first energy discrimination threshold, and at least one second energy discrimination threshold; The data in the energy spectrum information whose energy is greater than the lower limit of energy discrimination and less than the first energy discrimination threshold, and whose energy is greater than the second energy discrimination threshold and less than the upper limit of energy discrimination, are taken as the first data subset, and the data in the energy spectrum information whose energy is greater than the first energy discrimination threshold and less than the second energy discrimination threshold are taken as the second data subset. The first data subset and the second data subset constitute the at least two data subsets. The step of determining preliminary scattering correction data based on the at least two data subsets and a preset scattering estimation model includes: Determine the number of events that match the events in the first data subset; Based on the number of coincidence events in the first data subset and the preset scattering estimation model, the scattering coincidence event estimates in the second data subset are determined. The scattering coincidence event estimates in the second data subset are then superimposed with the coincidence events in the first data subset to obtain the preliminary scattering correction data.

2. The scattering correction method according to claim 1, characterized in that, After determining the preliminary scattering correction data, the process also includes: Based on a preset scattering estimation model, estimated scattering correction data is obtained, and the preliminary scattering correction data is optimized based on the estimated scattering correction data to generate optimized scattering correction data.

3. The scattering correction method according to claim 1, characterized in that, After determining the preliminary scattering correction data, the process also includes: Based on the preliminary scattering correction data and the preset image reconstruction algorithm, a radioactivity distribution map is obtained; Obtain the attenuation information distribution map, and optimize the preliminary scattering correction data based on the attenuation information distribution map and the radioactivity distribution map to generate optimized scattering correction data.

4. The scattering correction method according to claim 1, characterized in that, After determining the preliminary scattering correction data, the process also includes: Based on the preliminary scattering correction data and the preset image reconstruction algorithm, a radioactivity distribution map is obtained; A deep learning network is constructed, and the radioactivity distribution map is optimized based on the deep learning network to generate an optimized radioactivity distribution map; Obtain the attenuation information distribution map, and optimize the preliminary scattering correction data based on the attenuation information distribution map and the optimized radioactivity distribution map to generate optimized scattering correction data.

5. The scattering correction method according to claim 3, characterized in that, The step of optimizing the preliminary scattering correction data based on the attenuation information distribution map and the radioactivity distribution map to generate optimized scattering correction data includes: The preliminary scattering correction data is optimized based on the attenuation information distribution map and the radioactivity distribution map to generate transition scattering correction data; Determine whether the number of optimization attempts for the initial scattering correction data exceeds the threshold number of optimization attempts; If the number of optimization attempts is greater than the threshold number of optimization attempts, then the transition scattering correction data is the optimized scattering correction data; If the number of optimization attempts is less than or equal to the threshold number of optimization attempts, then the transition scattering correction data is optimized based on the attenuation information distribution map and the radioactivity distribution map.

6. A scattering correction device, characterized in that, include: The data subset acquisition unit is used to obtain at least two data subsets based on the energy spectrum information of PET data; A scattering correction unit is used to determine preliminary scattering correction data based on the at least two data subsets and a preset scattering estimation model; The energy spectrum information based on PET data yields at least two data subsets, including: Determine the upper limit of energy discrimination, the lower limit of energy discrimination, at least one first energy discrimination threshold, and at least one second energy discrimination threshold; The data in the energy spectrum information whose energy is greater than the lower limit of energy discrimination and less than the first energy discrimination threshold, and whose energy is greater than the second energy discrimination threshold and less than the upper limit of energy discrimination, are taken as the first data subset, and the data in the energy spectrum information whose energy is greater than the first energy discrimination threshold and less than the second energy discrimination threshold are taken as the second data subset. The first data subset and the second data subset constitute the at least two data subsets. The step of determining preliminary scattering correction data based on the at least two data subsets and a preset scattering estimation model includes: Determine the number of events that match the events in the first data subset; Based on the number of coincidence events in the first data subset and the preset scattering estimation model, the scattering coincidence event estimates in the second data subset are determined. The scattering coincidence event estimates in the second data subset are then superimposed with the coincidence events in the first data subset to obtain the preliminary scattering correction data.

7. An electronic device, characterized in that, Including memory and processor, among which, The memory is used to store programs; The processor, coupled to the memory, is used to execute the program stored in the memory to implement the steps of the scattering correction method according to any one of claims 1 to 5.

8. A computer-readable storage medium, characterized in that, Used to store computer-readable programs or instructions, which, when executed by a processor, can implement the steps of the scattering correction method according to any one of claims 1 to 5.

Images