WO2021004157A1 - Medical scanning imaging method and apparatus, storage medium, and computer device - Google Patents

Abstract

A medical scanning imaging method and apparatus, a storage medium, and a computer device, comprising: obtaining an attenuation graph obtained according to a medical scanning process of a target object (S100); determining a region of interest according to the attenuation graph, and determining a scattering response function of a target pixel point in the region of interest (S200); obtaining a corresponding original chord graph according to dynamic information of the target object (S300); obtaining a dynamic equation corresponding to the original chord graph according to a dynamic image model and the scattering response function (S400); and performing dynamic reconstruction processing according to the attenuation graph, the scattering response function, the original chord graph, and the dynamic equation corresponding to the original chord graph to obtain a dynamic image corresponding to the target object (S500). When scattering correction is carried out, scattering estimation is carried out in the dynamic reconstruction process without using an activity graph, so that the scattering correction efficiency can be improved. Besides, the scattering response function and the dynamic image model are combined so that low-noise scattering estimation can be provided in dynamic reconstruction, thereby improving the image reconstruction accuracy and ensuring the image quality.

Description

Medical scanning imaging method, device, storage medium and computer equipment

cross reference

This application claims the priority of the Chinese application with the application number 201910616741.5 filed on July 9, 2019, and the entire content is incorporated herein by reference.

Technical field

This application relates to the field of image processing technology, and in particular to a medical scanning imaging method, device, storage medium and computer equipment.

Background technique

PET (Positron Emission Computed Tomography, Positron Emission Computed Tomography) is a relatively advanced clinical examination imaging technology in the field of nuclear medicine. This technology injects a substance labeled with a radionuclide into the human body, and detects the accumulation of the substance in the body's metabolism to reflect the body's life metabolism, thereby achieving the purpose of diagnosis.

When performing medical scanning of a target object based on PET scanning technology or multimodal scanning technology including PET scanning, it is necessary to perform scatter correction on the PET image obtained through image reconstruction processing to obtain correct quantitative results. The scatter correction methods in the prior art need to use activity maps and attenuation maps. However, in the dynamic reconstruction process, the activity maps at each time point are different. Therefore, each time point needs to be estimated separately. Thereby increasing the time of scattering correction and reducing the efficiency of scattering correction. In addition, a single high-noise image will bring high-noise scatter estimates, thereby reducing image quality.

Summary of the invention

Based on this, it is necessary to provide a medical scanning imaging method, device, storage medium and computer equipment with higher efficiency and accuracy in response to the problems in the prior art.

A medical scanning imaging method, including:

Obtain the attenuation map obtained according to the medical scanning process of the target object;

Determine the region of interest according to the attenuation map, and determine the scattering response function of the target pixel in the region of interest;

Obtaining the corresponding original chord diagram according to the dynamic information of the target object;

Obtaining the dynamic equation corresponding to the original chord diagram according to the dynamic image model and the scattering response function;

Perform dynamic reconstruction processing according to the attenuation map, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram, to obtain a dynamic image corresponding to the target object.

In an embodiment, the acquiring the attenuation map obtained according to the medical scanning process of the target object includes any one of the following items:

the first item:

Acquiring an attenuation map obtained according to a multi-modal medical scanning process of the target object, the multi-modal medical scanning process including PET scanning and other modal scanning, the attenuation map being obtained based on the scanning data of the other modal scanning;

second section:

Acquire PET scan data of the target object, and obtain a corresponding attenuation map according to the PET scan data.

In an embodiment, the determining the region of interest according to the attenuation map includes any one of the following items:

Item 1: Obtain the user's region of interest selection result, and determine the region of interest in the attenuation map according to the region of interest selection result;

Item 2: Determine the region of interest in the attenuation map by performing image segmentation processing on the attenuation map;

The third item: Define the area in the attenuation graph where the attenuation value is greater than the attenuation threshold as the region of interest.

In an embodiment, the determining the scattering response function of the target pixel in the region of interest includes any one of the following items:

the first item:

Determine the scattering response function of the target pixel in the region of interest at each chord diagram coordinate according to the attenuation map and the system geometric model corresponding to the medical scanning system;

second section:

According to the attenuation map and the system geometric model corresponding to the medical scanning system, the scattering response function of the target pixel in the region of interest at each chord diagram coordinate and each flight time interval is determined.

In an embodiment, the dynamic information includes: drug metabolism dynamic information or time dimension information generated by the movement of the target object.

In an embodiment, the dynamic image model includes at least one of a one-chamber model, a two-chamber model, a Patlak model, a rigid body motion model, and a non-rigid body motion model.

In an embodiment, performing dynamic reconstruction processing according to the attenuation map, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram to obtain the dynamic image corresponding to the target object includes:

According to the attenuation map, the scattering response function, the system geometric model corresponding to the medical scanning system, the original chord diagram and the dynamic equation corresponding to the original chord diagram, dynamic reconstruction processing is performed through the nested-maximum likelihood expectation algorithm To obtain the dynamic image corresponding to the target object.

A medical scanning imaging device includes: a first acquisition module, a first processing module, a second acquisition module, a second processing module, and a dynamic reconstruction module;

The first obtaining module is used to obtain an attenuation map obtained according to the medical scanning process of the target object;

The first processing module is configured to determine a region of interest according to the attenuation map, and determine a scattering response function of a target pixel in the region of interest;

The second acquiring module is used to acquire the corresponding original chord diagram according to the dynamic information of the target object;

The second processing module is used to obtain the dynamic equation corresponding to the original chord diagram according to the dynamic image model;

The dynamic reconstruction module is configured to perform dynamic reconstruction processing according to the attenuation map, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram, to obtain a dynamic image corresponding to the target object.

A computer device includes a memory and a processor, the memory stores a computer program, and the processor implements the steps of the foregoing method when the computer program is executed.

A computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the above method are realized.

The aforementioned medical scanning imaging method, device, storage medium and computer equipment are used to obtain the attenuation map obtained according to the medical scanning process of the target object; determine the region of interest according to the attenuation map, and determine the scattering response function of the target pixel in the region of interest; Obtain the corresponding original chord diagram according to the dynamic information of the target object; obtain the dynamic equation corresponding to the original chord diagram according to the dynamic image model and the scattering response function; perform according to the attenuation diagram, the scattering response function, the original chord diagram and the dynamic equation corresponding to the original chord diagram The dynamic reconstruction process obtains the dynamic image corresponding to the target object. When performing scatter correction, there is no need to use activity maps, but scatter estimation is performed during the dynamic reconstruction process, which can improve the efficiency of scatter correction. In addition, combining the scattering response function and the dynamic image model can provide low-noise scattering estimation in dynamic reconstruction, improve the accuracy of image reconstruction, and ensure image quality.

Description of the drawings

FIG. 1 is a schematic flowchart of a medical scanning imaging method in an embodiment;

2 is a schematic diagram of the structure of a medical scanning imaging device in an embodiment;

Figure 3 is an internal structure diagram of a computer device in an embodiment.

Detailed ways

In order to make the purpose, technical solutions, and advantages of this application clearer, the following further describes this application in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the application, and not used to limit the application.

In one embodiment, as shown in FIG. 1, a medical scanning imaging method is provided. The method is applied to a processor that can perform medical scanning imaging for explanation. The method includes the following steps:

Step S100: Obtain an attenuation map obtained according to the medical scanning process of the target object.

When the processor performs medical scan imaging, it first obtains the attenuation map of the target object. The attenuation map can be obtained by performing a medical scan on the target object. The medical scan can be a multi-modal scan, such as PET/CT (Positron Emission Computed Tomography/ Computed Tomography, Positron Emission Computed Tomography/Magnetic Resonance, PET/MR (Positron Emission Computed Tomography/Magnetic Resonance, Positron Emission Computed Tomography/Magnetic Resonance Scan), etc. Medical scanning can also be Single mode scan, such as PET scan, etc.

Step S200: Determine the region of interest according to the attenuation map, and determine the scattering response function of the target pixel in the region of interest.

The physical meaning of the scattering response function refers to the probability that an annihilation event on a pixel point in the image domain will be scattered through the influence of the attenuation map, and the probability of the scattered photons collected at the chord diagram point. The region of interest may be an image region containing the scanned part of the target object. In addition to the scanned part of the target object, the image obtained by the processor may also include other non-essential content. The non-essential content may have a certain interference effect on the medical analysis process of the target object. Therefore, you can determine the interest Area to remove other non-essential content.

After determining the region of interest in the attenuation map, the processor determines the scattering response function of the target pixel in the region of interest, and the scattering response function is used for scattering correction in the image reconstruction process. The target pixel may refer to all pixels in the region of interest. For example, if the image size of the region of interest is 255*255, then all pixels in the region of interest are selected as target pixels. The target pixel can also be a part of the pixel selected according to the actual situation. For example, if the image size of the region of interest is 255*255, then some pixels in the region of interest are extracted as the target pixel, and some pixels are The extraction process can be realized by the existing pixel extraction algorithm. The target pixel can also be the pixel selected after image processing of the original region of interest image. For example, the image size of the original region of interest image is 256*256, and the original region of interest image is down-sampled to obtain 128* For the 128 new image, select all or part of the pixels in the new image as target pixels; in addition, the image processing can also be other processing such as upsampling, which is not limited here.

Step S300: Obtain the corresponding original chord diagram according to the dynamic information of the target object.

During the process of medical scanning and imaging, the processor further includes the step of obtaining the corresponding original chord diagram according to the dynamic information of the target object. It should be noted that the execution order of the steps of step S300 and the previous joint step (the joint step refers to the step composed of step S100 and step S200) is not strictly limited. Step S300 and the joint step are independent steps, and can be performed at the same time. The execution may also be sequenced, and the specific execution order of the steps can be determined according to actual conditions.

Step S400: According to the dynamic image model and the scattering response function, a dynamic equation corresponding to the original chord diagram is obtained.

After obtaining the scattering response function of the target pixel in the region of interest of the attenuation map, the processor obtains the dynamic equation corresponding to the original chord diagram according to the dynamic image model.

Step S500: Perform dynamic reconstruction processing according to the attenuation map, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram, to obtain a dynamic image corresponding to the target object.

After the processor obtains the dynamic equation corresponding to the original chord diagram, it performs dynamic reconstruction processing according to the attenuation diagram, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram, so as to obtain the dynamic image corresponding to the target object and realize the dynamic image reconstruction.

It should be noted that the dynamic image reconstruction process in this embodiment refers to the process of obtaining corresponding dynamic images based on scan data. That is, in the method of this embodiment, only one dynamic image reconstruction process is required to complete medical Scatter correction of scanned images.

This embodiment provides a medical scanning imaging method. When performing scatter correction, no activity map is used, but scatter estimation is performed during the dynamic reconstruction process, thereby improving the efficiency of scattering correction. In addition, combining the scattering response function and the dynamic image model can provide low-noise scattering estimation in dynamic reconstruction, improve the accuracy of image reconstruction, and ensure image quality.

In one embodiment, obtaining the attenuation map obtained according to the medical scanning process of the target object includes: obtaining the attenuation map obtained according to the multi-modal medical scanning process of the target object. The multi-modal medical scanning process includes PET scanning and other modalities. Scan, the attenuation map is obtained based on the scan data of other modal scans.

Specifically, when the medical scan is a multi-modal scan including PET scans and other modal scans, the attenuation map is obtained based on the scan data of other modalities, for example: when using PET/CT to perform multi-modal medical scans on the target object In the process, the attenuation map can be obtained by CT; when using PET/MR to perform multi-modal medical scanning of the target object, the attenuation map can be obtained by MRI. Since the existing scanning protocol obtains the attenuation map before the patient undergoes the PET scan, the method can be performed simultaneously with the patient's scanning process, thereby improving efficiency.

In one embodiment, acquiring the attenuation map obtained according to the medical scanning process of the target object includes: acquiring the PET scan data of the target object, and obtaining the corresponding attenuation map according to the PET scan data. Since the existing scanning protocol obtains the attenuation map in the first half of the PET scan of the patient, the method can be performed simultaneously with the scanning process of the patient, thereby improving efficiency.

In one embodiment, determining the region of interest according to the attenuation map includes: obtaining a user's region of interest selection result, and determining the region of interest in the attenuation map according to the result of the region of interest selection. In the process of determining the region of interest, the region of interest can be manually guided, that is, the user selects the region of interest from the attenuation map through the interactive device, and the processor obtains the user's region of interest selection result, and determines according to the result of the region of interest selection The area of interest in the attenuation map.

In one embodiment, determining the region of interest according to the attenuation map includes: determining the region of interest in the attenuation map by performing image segmentation processing on the attenuation map. In the process of determining the region of interest, the processor may divide the attenuation map into different regions through image segmentation processing, and then select an appropriate region from the divided regions as the region of interest.

In one embodiment, determining the region of interest according to the attenuation map includes: defining a region in the attenuation map with an attenuation value greater than 0 as the region of interest. In the process of determining the region of interest, the processor may determine the region of interest according to the attenuation value. Since the attenuation value of the image of unnecessary content is usually less than 0, the attenuation value of the image of the scanned part of the target object is greater than 0, therefore, The area where the attenuation value is greater than 0 in the attenuation map can be defined as the region of interest. For example, the contour map of the patient can be generated based on the attenuation map, and the area within the contour is the region of interest.

In one embodiment, determining the scattering response function of the target pixel in the region of interest includes: determining the coordinates of the target pixel in the region of interest at each chord diagram and each coordinate according to the attenuation map and the system geometric model corresponding to the medical scanning system. The scattering response function over a flight time interval.

When performing medical scans on target objects, PET scans can be divided into TOF (Time of Flight)-PET and non-TOF-PET. Before positron emission tomography, it is necessary to inject a radioactive tracer (such as fluoroglucose) into the target object, and the tracer can be metabolized by human tissues. Compared with normal tissues, tumors have a higher level of metabolism. The principle of PET imaging is: the decay of the tracer produces positrons, and the annihilation of the positron and the negative electron emits two pairs of photons with opposite directions and equal energy. Each photon flies at the speed of light. After the detector detects the photon pair, it performs a series of signal processing to reconstruct an image with clinical diagnostic significance. If the time difference between the two photons' arrival at the detector can be measured, since the diameter of the detector and the speed of light are known, it is possible to determine the position of the photon, that is, the emission position of the positron, that is, the position of the tracer decay. Call this technology Time of Flight (TOF). The photon occurrence position can be calculated by the photon flight time difference, the detector diameter and the speed of light: Δx=Δt*C/2, Δx represents the distance between the annihilation position and the center of the detector, Δt represents the flight time difference of two photons, and C represents the speed of light.

This embodiment uses the TOF-PET scanning mode. When calculating the scattering response function, the processor determines the target pixel point in the region of interest on each chord diagram coordinate and according to the attenuation map and the system geometric model corresponding to the medical scanning system. The scattering response function in each flight time interval, the scattering response function is used for the scattering correction during the image dynamic reconstruction process.

In one embodiment, determining the scattering response function of the target pixel in the region of interest includes: determining the coordinates of the target pixel in the region of interest on each chord diagram coordinate according to the attenuation map and the system geometric model corresponding to the medical scanning system Scattering response function.

In this embodiment, the non-TOF-PET scanning mode is adopted. When calculating the scattering response function, the processor calculates the target pixel in the region of interest according to the attenuation map and the system geometric model corresponding to the medical scanning system. The scattering response function on the coordinate, the scattering response function is used for the scattering correction during the image dynamic reconstruction process.

In one embodiment, the method of calculating the scattering response function is not limited. The specific implementation method can also be determined according to the preset PET scan time. When a longer PET scan is scheduled, a slow calculation speed but high accuracy algorithm can be selected ( Such as Monte Carlo method), on the contrary, you can choose an algorithm with fast calculation speed but lower accuracy.

In one embodiment, when the corresponding original chord diagram is obtained according to the dynamic information of the target object, the dynamic information includes: drug metabolism dynamic information or time dimension information generated by the movement of the target object. Specifically, the drug metabolism dynamic information may refer to the dynamic information caused by the metabolic process of the drug in the target object; the time dimension information may refer to the time dimension information generated due to the movement of the target object (such as breathing exercise), which is not done here. Specific restrictions.

In one embodiment, when the dynamic equation corresponding to the original chord diagram is obtained according to the dynamic image model and the scattering response function, the dynamic image model includes: one-chamber model, two-chamber model, Patlak model, rigid body motion model, and non-rigid body motion model. At least one of. Specifically, one or more of the models may be used according to actual needs to obtain the dynamic equation corresponding to the original chord diagram, which is not specifically limited here.

In one embodiment, the dynamic reconstruction process is performed according to the attenuation map, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram to obtain the dynamic image corresponding to the target object, including: according to the attenuation diagram, the scattering response function, and the medical scan The system geometric model, the original chord diagram and the dynamic equations corresponding to the original chord diagram are dynamically reconstructed by the nested-ML-EM (nested-Maximum Likelihood Expectation Maximum) algorithm to obtain the target object The corresponding dynamic image.

In an embodiment, taking the dynamic image model as a one-room model as an example, when the dynamic information is drug metabolism dynamic information, the dynamic PET image can be expressed as:

Among them, x(t) represents a dynamic PET image, v _b represents a plasma ratio image, C _p (t) represents a plasma input function, K ₁ and K ₂ represent a forward parameter image and a reverse parameter image, and t represents time. The above formula can use 3 images to describe a group of dynamic images.

When the dynamic information is the time dimension information generated by the motion of the target object, images with different motion phases can be expressed as:

x(t)=TX (2)

Among them, x(t) represents a dynamic PET image, T represents a conversion matrix, X represents a reference image without motion, and t represents a motion phase. T can be obtained by roughly reconstructing the image and by image registration.

In addition, when the dynamic information includes both the drug metabolism dynamic information and the time dimension information generated by the movement of the target object, the images of different movement phases and different time points can be obtained by combining the above formulas (1) and (2). In different situations Below, dynamic images can be expressed as:

x(t)=f(p ₁ ,p ₂ …p _n ,q ₁ ,q ₂ …q _m ,t)

Among them, x(t) represents a dynamic PET image, p ₁ , p ₂ …p _n represents an unknown image, such as v _b , K ₁ , K ₂ in formula (1), and X, q ₁ , in formula (2) q ₂ …q _m represents known parameters, such as C _p (t) in formula (1) and T in formula (2).

Then the dynamic equation corresponding to the scattering chord diagram can be expressed as:

The projection relationship between the acquired original chord diagram y _j (t) and the dynamic PET image x _i (t) is:

y _j (t)=a _j H _i,j x _i (t)+s _j (t)+r _j

Combining the above projection relationship with the dynamic equation x(t) of the previously acquired dynamic PET image, the dynamic equation based on the original chord diagram is obtained:

y _j (t)=a _j H _i,j f(p ₁ ,p ₂ …p _n ,q ₁ ,q ₂ …q _m ,t)+s _j (t)+r _j

Further, in combination with the expression of the dynamic equation s _j (t) corresponding to the scattering chord diagram, we can get:

According to the dynamic equation between the original chord diagram y _j (t) and the dynamic parameters, an iterative function can be obtained.

In the dynamic reconstruction process, the nested-ML-EM method is used, combined with the dynamic equation corresponding to the scattered chord diagram and the dynamic equation of the original chord diagram, and the final dynamic image calculation formula is as follows:

Among them, i represents the coordinates of the pixels in the image domain, j represents the coordinates of the pixels in the projection domain, _Hi,j represents the system geometric model corresponding to the medical scanning system, S _i,j represents the scattering response function, a _j represents the attenuation chord diagram, r _j Represents the chord diagram of random events, y _j represents the original data, n represents the total number of iterations, and m represents the number of nested iterations.

The iteration method of the k-th parameter image is:

Among them, G _k is an iterative equation related to f.

Under reasonable conditions, it should be understood that, although the steps in the flowcharts involved in the foregoing embodiments are displayed in sequence as indicated by the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless specifically stated in this article, the execution of these steps is not strictly limited in order, and these steps can be executed in other orders. Moreover, at least part of the steps in each flowchart may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The order of execution is not necessarily performed sequentially, but may be performed alternately or alternately with other steps or at least a part of sub-steps or stages of other steps.

In one embodiment, as shown in FIG. 2, a medical scanning imaging device is provided, which includes: a first acquisition module 100, a first processing module 200, a second acquisition module 300, a second processing module 400, and dynamic reconstruction Module 500;

The first acquiring module 100 is configured to acquire the attenuation map obtained according to the medical scanning process of the target object;

The first processing module 200 is configured to determine the region of interest according to the attenuation map, and determine the scattering response function of the target pixel in the region of interest;

The second obtaining module 300 is configured to obtain the corresponding original chord diagram according to the dynamic information of the target object;

The second processing module 400 is used to obtain the dynamic equation corresponding to the original chord diagram according to the dynamic image model;

The dynamic reconstruction module 500 is used to perform dynamic reconstruction processing according to the attenuation diagram, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram, to obtain a dynamic image corresponding to the target object.

For the specific definition of the medical scanning imaging device, please refer to the above definition of the medical scanning imaging method, which will not be repeated here. Each module in the above medical scanning imaging device can be implemented in whole or in part by software, hardware, and a combination thereof. The foregoing modules may be embedded in the form of hardware or independent of the processor in the computer device, or may be stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to the foregoing modules.

In one embodiment, a computer device is provided, including a memory and a processor, and a computer program is stored in the memory. When the processor executes the computer program, the following steps are implemented: acquiring an attenuation map obtained according to a medical scanning process of a target object; Determine the region of interest according to the attenuation map, and determine the scattering response function of the target pixel in the region of interest; obtain the corresponding original chord diagram according to the dynamic information of the target object; obtain the corresponding original chord diagram according to the dynamic image model and the scattering response function Dynamic equation: Perform dynamic reconstruction processing according to the attenuation diagram, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram to obtain a dynamic image corresponding to the target object.

In an embodiment, the processor further implements any one of the following items when executing the computer program:

The first item: Obtain the attenuation map obtained according to the multi-modal medical scanning process of the target object. The multi-modal medical scanning process includes PET scanning and other modal scanning. The attenuation map is obtained based on the scanning data of other modal scanning;

The second item: Obtain the PET scan data of the target object, and obtain the corresponding attenuation map according to the PET scan data.

In an embodiment, the processor further implements any one of the following items when executing the computer program:

The first item: Obtain the user's region of interest selection result, and determine the region of interest in the attenuation map according to the result of the region of interest selection;

The second item: Determine the region of interest in the attenuation map by performing image segmentation processing on the attenuation map;

The third item: Define the area in the attenuation graph whose attenuation value is greater than the attenuation threshold as the region of interest.

In an embodiment, the processor further implements any one of the following items when executing the computer program:

The first item: Determine the scattering response function of the target pixel in the region of interest at each chord diagram coordinate according to the attenuation map and the corresponding system geometric model of the medical scanning system;

The second item: Determine the scattering response function of the target pixel in the region of interest at each chord diagram coordinate and each flight time interval according to the attenuation map and the system geometric model corresponding to the medical scanning system.

In one embodiment, the processor further implements the following steps when executing the computer program: according to the attenuation map, the scattering response function, the system geometric model corresponding to the medical scanning system, the original chord diagram, and the dynamic equations corresponding to the original chord diagram, through nesting- The maximum likelihood expectation algorithm performs dynamic reconstruction processing to obtain the dynamic image corresponding to the target object.

Fig. 3 shows an internal structure diagram of a computer device in an embodiment. The computer device may specifically be a terminal (or server). As shown in Figure 3, the computer equipment includes a processor, a memory, a network interface, an input device, and a display screen connected through a system bus. Among them, the memory includes a non-volatile storage medium and an internal memory. The non-volatile storage medium of the computer device stores an operating system and can also store a computer program. When the computer program is executed by the processor, the processor can realize the video bit rate control method and the video transcoding method. A computer program may also be stored in the internal memory. When the computer program is executed by the processor, the processor can execute the video rate control method and the video transcoding method. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen. The input device of the computer equipment can be a touch layer covered on the display screen, or a button, trackball or touch pad set on the housing of the computer equipment. It can be an external keyboard, touchpad, or mouse.

Those skilled in the art can understand that the structure shown in FIG. 3 is only a block diagram of part of the structure 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 Including more or fewer parts than shown in the figure, or combining some parts, or having a different arrangement of parts.

In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the following steps are implemented: acquiring an attenuation map obtained according to a medical scanning process of a target object; Determine the region of interest, and determine the scattering response function of the target pixel in the region of interest; obtain the corresponding original chord diagram according to the dynamic information of the target object; obtain the dynamic equation corresponding to the original chord diagram according to the dynamic image model and the scattering response function; According to the attenuation diagram, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram, the dynamic reconstruction process is performed to obtain the dynamic image corresponding to the target object.

In an embodiment, when the computer program is executed by the processor, any one of the following items is also implemented:

The first item: Obtain the attenuation map obtained according to the multi-modal medical scanning process of the target object. The multi-modal medical scanning process includes PET scanning and other modal scanning. The attenuation map is obtained based on the scanning data of other modal scanning;

The second item: Obtain the PET scan data of the target object, and obtain the corresponding attenuation map according to the PET scan data.

In an embodiment, when the computer program is executed by the processor, any one of the following items is also implemented:

The first item: Obtain the user's region of interest selection result, and determine the region of interest in the attenuation map according to the result of the region of interest selection;

The second item: Determine the region of interest in the attenuation map by performing image segmentation processing on the attenuation map;

The third item: Define the area in the attenuation graph where the attenuation value is greater than the attenuation threshold as the region of interest.

In an embodiment, when the computer program is executed by the processor, any one of the following items is also implemented:

The first item: Determine the scattering response function of the target pixel in the region of interest at each chord diagram coordinate according to the attenuation map and the corresponding system geometric model of the medical scanning system;

The second item: Determine the scattering response function of the target pixel in the region of interest at each chord diagram coordinate and each flight time interval according to the attenuation map and the system geometric model corresponding to the medical scanning system.

In one embodiment, when the computer program is executed by the processor, the following steps are also implemented: according to the attenuation map, the scattering response function, the system geometric model corresponding to the medical scanning system, the original chord diagram and the dynamic equations corresponding to the original chord diagram, through nesting -The maximum likelihood expectation algorithm performs dynamic reconstruction processing to obtain a dynamic image corresponding to the target object.

A person of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing relevant hardware through a computer program, which can be stored in a non-volatile computer readable storage medium. When the computer program is executed, it may include the processes of the above-mentioned method embodiments. Wherein, any reference to memory, storage, database or other media used in the embodiments provided in this application may include non-volatile and/or volatile memory. Non-volatile memory may include read only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. As an illustration and not a limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous chain Channel (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

The technical features of the above-mentioned embodiments can be combined arbitrarily. In order to make the description concise, all possible combinations of the various technical features in the above-mentioned embodiments are not described. However, as long as there is no contradiction in the combination of these technical features, All should be considered as the scope of this specification.

The above-mentioned embodiments only express several embodiments of the present invention, and the descriptions are more specific and detailed, but they should not be understood as limiting the scope of the invention patent. It should be pointed out that for those of ordinary skill in the art, without departing from the concept of the present invention, several modifications and improvements can be made, and these all fall within the protection scope of the present invention. Therefore, the protection scope of the patent of the present invention should be subject to the appended claims.

Claims (10)

1. A medical scanning imaging method is characterized in that it comprises:

   Obtain the attenuation map obtained according to the medical scanning process of the target object;

   Determine the region of interest according to the attenuation map, and determine the scattering response function of the target pixel in the region of interest;

   Obtaining the corresponding original chord diagram according to the dynamic information of the target object;

   Obtaining the dynamic equation corresponding to the original chord diagram according to the dynamic image model and the scattering response function;

   Perform dynamic reconstruction processing according to the attenuation map, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram, to obtain a dynamic image corresponding to the target object.
2. The medical scanning imaging method according to claim 1, wherein the acquiring the attenuation map obtained according to the medical scanning process of the target object includes any one of the following items:

   the first item:

   Acquiring an attenuation map obtained according to a multi-modal medical scanning process of the target object, the multi-modal medical scanning process including PET scanning and other modal scanning, the attenuation map being obtained based on the scanning data of the other modal scanning;

   second section:

   Acquire PET scan data of the target object, and obtain a corresponding attenuation map according to the PET scan data.
3. The medical scanning imaging method according to claim 1, wherein the determining the region of interest according to the attenuation map includes any one of the following items:

   Item 1: Obtain the user's region of interest selection result, and determine the region of interest in the attenuation map according to the region of interest selection result;

   Item 2: Determine the region of interest in the attenuation map by performing image segmentation processing on the attenuation map;

   The third item: Define the area in the attenuation graph where the attenuation value is greater than the attenuation threshold as the region of interest.
4. The medical scanning imaging method according to claim 1, wherein the determining the scattering response function of the target pixel in the region of interest includes any one of the following:

   the first item:

   Determine the scattering response function of the target pixel in the region of interest at each chord diagram coordinate according to the attenuation map and the system geometric model corresponding to the medical scanning system;

   second section:

   According to the attenuation map and the system geometric model corresponding to the medical scanning system, the scattering response function of the target pixel in the region of interest at each chord diagram coordinate and each flight time interval is determined.
5. The medical scanning imaging method according to claim 1, wherein the dynamic information comprises: drug metabolism dynamic information or time dimension information generated by the movement of the target object.
6. The medical scanning imaging method of claim 1, wherein the dynamic image model comprises at least one of a one-chamber model, a two-chamber model, a Patlak model, a rigid body motion model, and a non-rigid body motion model.
7. The medical scanning imaging method of claim 1, wherein the dynamic reconstruction process is performed according to the attenuation map, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram to obtain the The dynamic image corresponding to the target object includes:

   According to the attenuation map, the scattering response function, the system geometric model corresponding to the medical scanning system, the original chord diagram and the dynamic equation corresponding to the original chord diagram, dynamic reconstruction processing is performed through the nested-maximum likelihood expectation algorithm To obtain the dynamic image corresponding to the target object.
8. A medical scanning imaging device, characterized by comprising: a first acquisition module, a first processing module, a second acquisition module, a second processing module, and a dynamic reconstruction module;

   The first obtaining module is used to obtain an attenuation map obtained according to the medical scanning process of the target object;

   The first processing module is configured to determine a region of interest according to the attenuation map, and determine a scattering response function of a target pixel in the region of interest;

   The second obtaining module is configured to obtain the corresponding original chord diagram according to the dynamic information of the target object;

   The second processing module is used to obtain the dynamic equation corresponding to the original chord diagram according to the dynamic image model;

   The dynamic reconstruction module is configured to perform dynamic reconstruction processing according to the attenuation map, the scattering response function, the original chord diagram, and the dynamic equation corresponding to the original chord diagram, to obtain a dynamic image corresponding to the target object.
9. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method in any one of claims 1 to 7 when the computer program is executed by the processor.
10. A computer-readable storage medium with a computer program stored thereon, wherein the computer program implements the steps of the method according to any one of claims 1 to 7 when the computer program is executed by a processor.

Images