WO2021004157A1 - Procédé et appareil d'imagerie par balayage médicale, support de stockage et dispositif informatique - Google Patents
Procédé et appareil d'imagerie par balayage médicale, support de stockage et dispositif informatique Download PDFInfo
- Publication number
- WO2021004157A1 WO2021004157A1 PCT/CN2020/090865 CN2020090865W WO2021004157A1 WO 2021004157 A1 WO2021004157 A1 WO 2021004157A1 CN 2020090865 W CN2020090865 W CN 2020090865W WO 2021004157 A1 WO2021004157 A1 WO 2021004157A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- dynamic
- interest
- region
- attenuation map
- target object
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Definitions
- 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.
- 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.
- a medical scanning imaging method including:
- the acquiring the attenuation map obtained according to the medical scanning process of the target object includes any one of the following items:
- 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;
- 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.
- the determining the scattering response function of the target pixel in the region of interest includes any one of the following items:
- 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.
- the dynamic information includes: drug metabolism dynamic information or time dimension information generated by the movement of the target object.
- 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.
- 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:
- 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.
- FIG. 1 is a schematic flowchart of a medical scanning imaging method in an embodiment
- FIG. 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.
- 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.
- the processor 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.
- 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.
- 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.
- 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*
- 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.
- the processor further includes the step of obtaining the corresponding original chord diagram according to the dynamic information of the target object.
- the execution order of the steps of step S300 and the previous joint step 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.
- the processor 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.
- the processor 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.
- 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.
- 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.
- 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.
- 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
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- PET scans can be divided into TOF (Time of Flight)-PET and non-TOF-PET.
- TOF Time of Flight
- a radioactive tracer such as fluoroglucose
- 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.
- This embodiment uses the TOF-PET scanning mode.
- 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.
- 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.
- the non-TOF-PET scanning mode is adopted.
- 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.
- 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.
- the dynamic information 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.
- 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.
- the dynamic image model when the dynamic equation corresponding to the original chord diagram is obtained according to the dynamic image model and the scattering response function, 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.
- 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.
- the dynamic PET image when the dynamic information is drug metabolism dynamic information, the dynamic PET image can be expressed as:
- x(t) represents a dynamic PET image
- v b represents a plasma ratio image
- C p (t) represents a plasma input function
- K 1 and K 2 represent a forward parameter image and a reverse parameter image
- t represents time.
- the above formula can use 3 images to describe a group of dynamic images.
- images with different motion phases can be expressed as:
- x(t) represents a dynamic PET image
- T represents a conversion matrix
- X represents a reference image without motion
- t represents a motion phase.
- T can be obtained by roughly reconstructing the image and by image registration.
- 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).
- dynamic images can be expressed as:
- x(t) represents a dynamic PET image
- p 1 , p 2 ...p n represents an unknown image, such as v b , K 1 , K 2 in formula (1)
- q 2 ...q m represents known parameters, such as C p (t) in formula (1) and T in formula (2).
- y j (t) a j H i,j f(p 1 ,p 2 ...p n ,q 1 ,q 2 ...q m ,t)+s j (t)+r j
- 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
- m represents the number of nested iterations.
- the iteration method of the k-th parameter image is:
- G k is an iterative equation related to f.
- a medical scanning imaging device 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.
- 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.
- a computer device including a memory and a processor, and a computer program is stored in the memory.
- 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.
- 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.
- 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.
- 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.
- 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).
- the computer equipment includes a processor, a memory, a network interface, an input device, and a display screen connected through a system bus.
- 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.
- 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.
- 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.
- 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.
- a computer-readable storage medium on which a computer program is stored.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- SRAM static RAM
- DRAM dynamic RAM
- SDRAM synchronous DRAM
- DDRSDRAM double data rate SDRAM
- ESDRAM enhanced SDRAM
- SLDRAM synchronous chain Channel
- memory bus Radbus direct RAM
- RDRAM direct memory bus dynamic RAM
- RDRAM memory bus dynamic RAM
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Nuclear Medicine (AREA)
Abstract
L'invention concerne un procédé et un appareil d'imagerie par balayage médicale, un support de stockage et un dispositif informatique, consistant à : obtenir un graphe d'atténuation obtenu conformément à un processus de balayage médical d'un objet cible (S100) ; déterminer une région d'intérêt en fonction du graphe d'atténuation, et déterminer une fonction de réponse de diffusion d'un point de pixel cible dans la région d'intérêt (S200) ; obtenir un graphe de corde d'origine correspondant en fonction d'informations dynamiques de l'objet cible (S300) ; obtenir une équation dynamique correspondant au graphe de corde d'origine selon un modèle d'image dynamique et la fonction de réponse de diffusion (S400) ; et effectuer un traitement de reconstruction dynamique selon le graphe d'atténuation, la fonction de réponse de diffusion, le graphe de corde d'origine et l'équation dynamique correspondant au graphe de corde d'origine pour obtenir une image dynamique correspondant à l'objet cible (S500). Lorsqu'une correction de diffusion est effectuée, une estimation de diffusion est effectuée dans le processus de reconstruction dynamique sans utiliser de graphe d'activité, de telle sorte que l'efficacité de correction de diffusion peut être améliorée. En outre, la fonction de réponse de diffusion et le modèle d'image dynamique sont combinés de telle sorte qu'une estimation de diffusion à bas bruit peut être obtenue dans la reconstruction dynamique, ce qui permet d'améliorer la précision de reconstruction d'image et d'assurer la qualité d'image.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910616741.5A CN110415310B (zh) | 2019-07-09 | 2019-07-09 | 医学扫描成像方法、装置、存储介质及计算机设备 |
CN201910616741.5 | 2019-07-09 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2021004157A1 true WO2021004157A1 (fr) | 2021-01-14 |
Family
ID=68360790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2020/090865 WO2021004157A1 (fr) | 2019-07-09 | 2020-05-18 | Procédé et appareil d'imagerie par balayage médicale, support de stockage et dispositif informatique |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN110415310B (fr) |
WO (1) | WO2021004157A1 (fr) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113674254A (zh) * | 2021-08-25 | 2021-11-19 | 上海联影医疗科技股份有限公司 | 医学图像异常点识别方法、设备、电子装置和存储介质 |
US12008689B2 (en) | 2021-12-03 | 2024-06-11 | Canon Medical Systems Corporation | Devices, systems, and methods for deep-learning kernel-based scatter estimation and correction |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110415310B (zh) * | 2019-07-09 | 2022-12-20 | 上海联影医疗科技股份有限公司 | 医学扫描成像方法、装置、存储介质及计算机设备 |
CN111904379B (zh) * | 2020-07-13 | 2024-04-12 | 上海联影医疗科技股份有限公司 | 多模态医学设备的扫描方法和装置 |
CN112017258B (zh) * | 2020-09-16 | 2024-04-30 | 上海联影医疗科技股份有限公司 | Pet图像重建方法、装置、计算机设备以及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102496175A (zh) * | 2011-12-22 | 2012-06-13 | 中国科学院高能物理研究所 | 基于计算机断层成像ct创建被测量物衰减图的方法及装置 |
US20130141098A1 (en) * | 2011-07-19 | 2013-06-06 | Siemens Aktiengesellschaft | Alignment Phantom for MR/PET System |
CN103876772A (zh) * | 2014-03-20 | 2014-06-25 | 中北大学 | 一种多谱成像方法和装置 |
CN105678711A (zh) * | 2016-01-29 | 2016-06-15 | 中国科学院高能物理研究所 | 一种基于图像分割的衰减校正方法 |
CN110415310A (zh) * | 2019-07-09 | 2019-11-05 | 上海联影医疗科技有限公司 | 医学扫描成像方法、装置、存储介质及计算机设备 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7813783B2 (en) * | 2006-11-02 | 2010-10-12 | General Electric Company | Methods and systems for attenuation correction in medical imaging |
CN106415317B (zh) * | 2014-06-13 | 2019-11-19 | 美国西门子医疗解决公司 | 单光子发射计算机化断层摄影术中的多个发射能量 |
CN106491153B (zh) * | 2016-12-29 | 2017-10-27 | 上海联影医疗科技有限公司 | 一种pet散射校正方法、pet成像方法及pet成像系统 |
-
2019
- 2019-07-09 CN CN201910616741.5A patent/CN110415310B/zh active Active
-
2020
- 2020-05-18 WO PCT/CN2020/090865 patent/WO2021004157A1/fr active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130141098A1 (en) * | 2011-07-19 | 2013-06-06 | Siemens Aktiengesellschaft | Alignment Phantom for MR/PET System |
CN102496175A (zh) * | 2011-12-22 | 2012-06-13 | 中国科学院高能物理研究所 | 基于计算机断层成像ct创建被测量物衰减图的方法及装置 |
CN103876772A (zh) * | 2014-03-20 | 2014-06-25 | 中北大学 | 一种多谱成像方法和装置 |
CN105678711A (zh) * | 2016-01-29 | 2016-06-15 | 中国科学院高能物理研究所 | 一种基于图像分割的衰减校正方法 |
CN110415310A (zh) * | 2019-07-09 | 2019-11-05 | 上海联影医疗科技有限公司 | 医学扫描成像方法、装置、存储介质及计算机设备 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113674254A (zh) * | 2021-08-25 | 2021-11-19 | 上海联影医疗科技股份有限公司 | 医学图像异常点识别方法、设备、电子装置和存储介质 |
CN113674254B (zh) * | 2021-08-25 | 2024-05-14 | 上海联影医疗科技股份有限公司 | 医学图像异常点识别方法、设备、电子装置和存储介质 |
US12008689B2 (en) | 2021-12-03 | 2024-06-11 | Canon Medical Systems Corporation | Devices, systems, and methods for deep-learning kernel-based scatter estimation and correction |
Also Published As
Publication number | Publication date |
---|---|
CN110415310A (zh) | 2019-11-05 |
CN110415310B (zh) | 2022-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021004157A1 (fr) | Procédé et appareil d'imagerie par balayage médicale, support de stockage et dispositif informatique | |
US9754390B2 (en) | Reconstruction of time-varying data | |
US10395353B2 (en) | Model-based scatter in multi-modality multi-energy SPECT reconstruction | |
Novosad et al. | MR-guided dynamic PET reconstruction with the kernel method and spectral temporal basis functions | |
US9155514B2 (en) | Reconstruction with partially known attenuation information in time of flight positron emission tomography | |
US8457380B2 (en) | PET local tomography | |
US8068650B2 (en) | Lesion quantification and tracking using multiple modalities | |
EP2399238B1 (fr) | Imagerie fonctionnelle | |
US10863946B2 (en) | Respiratory motion estimation in projection domain in nuclear medical imaging | |
EP2814395B1 (fr) | Reconstruction d'image nucléaire spatialement corrigée | |
CN110136076B (zh) | 医学扫描成像方法、装置、存储介质及计算机设备 | |
US10064593B2 (en) | Image reconstruction for a volume based on projection data sets | |
US20140133707A1 (en) | Motion information estimation method and image generation apparatus using the same | |
EP3631762A1 (fr) | Systèmes et procédés pour fournir des valeurs de confiance en tant que mesure d'assurance quantitative pour des images reconstruites de manière itérative en tomographie d'émission | |
US11354831B2 (en) | Image reconstruction using tracer and various parameters | |
CN112365479A (zh) | Pet参数图像处理方法、装置、计算机设备及存储介质 | |
CN105631908B (zh) | 一种pet图像重建方法和装置 | |
Reader | The promise of new PET image reconstruction | |
Feng et al. | Non-rigid respiratory motion correction for 4D gated PET sinogram data | |
CN115774281A (zh) | 用于断层摄影不一致性检测的基于衰减校正的加权 | |
JP7459243B2 (ja) | 1以上のニューラルネットワークとしての画像形成のモデル化による画像再構成 | |
US11311207B2 (en) | Systems and methods for pulmonary ventilation from image processing | |
Tsui et al. | Advances in 4D gated cardiac PET imaging for image quality improvement and cardiac motion and contractility estimation | |
US20230069017A1 (en) | Systems and methods for image reconstruction | |
CN113674377B (zh) | 正电子发射断层扫描图像重建的系统和方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20837124 Country of ref document: EP Kind code of ref document: A1 |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 20837124 Country of ref document: EP Kind code of ref document: A1 |