CN113397581A - Method and device for reconstructing medical dynamic image - Google Patents
Method and device for reconstructing medical dynamic image Download PDFInfo
- Publication number
- CN113397581A CN113397581A CN202110953441.3A CN202110953441A CN113397581A CN 113397581 A CN113397581 A CN 113397581A CN 202110953441 A CN202110953441 A CN 202110953441A CN 113397581 A CN113397581 A CN 113397581A
- Authority
- CN
- China
- Prior art keywords
- interest
- time
- activity
- region
- curve
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine (AREA)
Abstract
本发明提供了一种重建医学动态影像的方法，包括：提取医学动态影像的感兴趣区域，医学动态影像包括在多个时段拍摄的人体部位的影像；获取感兴趣区域所注入的放射性药物的活度随时间变化的曲线；根据曲线确定感兴趣时段，感兴趣时段包含活度的最大值所在的时刻；根据感兴趣区域的像素的活度变化，重建感兴趣区域在感兴趣时段内的影像。根据本发明的技术方案，通过对动态影像进行分析，重建出感兴趣区域最适合观测的时段的影像，以便阅片人能够直观地看到感兴趣区域的信号强度变化情况。
The present invention provides a method for reconstructing a medical dynamic image, which includes: extracting a region of interest of the medical dynamic image, where the medical dynamic image includes images of human body parts captured in multiple time periods; The curve of the degree of change with time; the period of interest is determined according to the curve, and the period of interest includes the moment when the maximum activity is located; the image of the region of interest within the period of interest is reconstructed according to the activity change of the pixels in the region of interest. According to the technical solution of the present invention, by analyzing the dynamic image, the image of the time period most suitable for observation of the region of interest is reconstructed, so that the reader can intuitively see the change of the signal intensity of the region of interest.
Description
Technical Field
The present invention relates to the field of medical imaging, and more particularly, to a method and apparatus for reconstructing a medical dynamic image.
Background
Nuclear medicine is an emerging discipline that employs nuclear technology to diagnose, treat, and study disease, including clinical and basic nuclear medicine. Nuclear medicine imaging technology has made breakthrough progress since the 70 s due to the development of single photon Emission Computed Tomography (PET) technology, and the innovation and development of radiopharmaceuticals. It complements and verifies with CT, nuclear magnetic resonance, ultrasonic technology and the like, and greatly improves the diagnosis and research level of diseases, so nuclear medicine imaging is a very active branch and an important component in the field of recent clinical medical image diagnosis.
In the course of clinical trials of new drugs, it is necessary to label them with isotopes in order to investigate various problems of drug metabolism, for example, in order to determine the active metabolites of drugs and to evaluate their pharmacological effects. For example, the metabolic status of a tracer drug in an organ of an animal can be followed by medical imaging.
However, the current image acquisition method has difficulty in ensuring that the acquisition time window is just the period of time when the tracer drug metabolism of the region of interest is most vigorous.
Disclosure of Invention
In order to solve the above problems, the present invention provides a method and an apparatus for reconstructing a medical dynamic image.
In a first aspect, an embodiment of the present invention provides a method for reconstructing a medical dynamic image, including: extracting an interested region of a medical dynamic image, wherein the medical dynamic image comprises images of a human body part shot in a plurality of time intervals; acquiring a curve of the activity of the radiopharmaceutical injected into the region of interest along with the change of time; determining an interesting time period according to the curve, wherein the interesting time period comprises the moment of the maximum value of the activity; and reconstructing an image of the region of interest in the interested period according to the activity change of the pixels of the region of interest.
In some embodiments of the invention, obtaining a time-dependent profile of activity of the radiopharmaceutical injected into the region of interest comprises: establishing a dynamic equation based on an atrioventricular model of the region of interest, wherein the atrioventricular model comprises a central chamber and peripheral chambers, and the dynamic equation comprises an equation of the change of the drug activity of the central chamber along with time and an equation of the change of the drug activity of the peripheral chambers along with time; determining a proportionality coefficient between an equation of time-varying drug activity in the central chamber and an equation of time-varying drug activity in the peripheral chambers; and determining a curve according to a kinetic equation and a proportionality coefficient.
In some embodiments of the invention, the curve is determined by the following equation: creal(T)=k1C1real(T)+k2C2real(T) wherein Creal(T) isActivity of the region of interest at time T, C1real(T) drug activity of the central compartment at time T, C2real(T) drug activity in the peripheral compartment at time T, k1Is the proportionality coefficient, k, of the equation relating the activity of the drug in the central compartment to the time2Proportionality coefficient, k, of equation of time-dependent change in drug activity in the peripheral compartment1And k2Satisfy k1+k2=1。
In some embodiments of the invention, determining the period of interest from the curve comprises: presetting an activity threshold according to the maximum value of the activity; the period of interest is determined from the maximum value of the activity and the activity threshold.
In some embodiments of the invention, determining the period of interest based on the maximum value of the activity and the activity threshold comprises: and determining a first time and a second time when the activity value on the curve is equal to the activity threshold value, wherein the first time is less than the second time, and the interested time period is from the first time to the second time.
In some embodiments of the invention, determining the period of interest from the curve comprises: determining the slope of a rising edge and the slope of a falling edge of the curve; the period of interest is determined based on the maximum value of the activity, the slope of the rising edge and the slope of the falling edge.
In a second aspect, an embodiment of the present invention provides an apparatus for reconstructing a medical dynamic image, including: the extraction module is used for extracting an interested region of a medical dynamic image, and the medical dynamic image comprises images of a human body part shot in a plurality of time intervals; the acquisition module is used for acquiring a curve of the change of the activity of the injected radiopharmaceutical in the region of interest along with time; a determining module, configured to determine, according to the curve, a time period of interest, where the time period of interest includes a time at which a maximum value of the activity is located; and the reconstruction module is used for reconstructing an image of the interested region in the interested time period according to the activity change of the pixels of the interested region.
In a third aspect, an embodiment of the present invention provides a computer-readable storage medium, which stores a computer program for executing the method for reconstructing a medical dynamic image according to the first aspect.
In a fourth aspect, an embodiment of the present invention provides an electronic device, including: a processor; a memory for storing processor executable instructions, wherein the processor is adapted to perform the method of reconstructing a medical dynamic image of the first aspect.
According to the technical scheme of the invention, the dynamic image is analyzed, and the image of the region of interest in the time period most suitable for observation is reconstructed, so that a reader can visually see the signal intensity change condition of the region of interest.
Drawings
Fig. 1 is a schematic diagram of an implementation environment for reconstructing a medical dynamic image according to an embodiment of the present invention.
Fig. 2 is a flowchart illustrating a method for reconstructing a medical dynamic image according to an embodiment of the present invention.
Fig. 3 is a schematic flowchart of calculating an activity curve of a region of interest according to another embodiment of the present invention.
Fig. 4 is a flowchart illustrating a process of determining a period of interest according to another embodiment of the present invention.
Fig. 5 is a flowchart illustrating a process of determining a period of interest according to another embodiment of the present invention.
Fig. 6a is a schematic diagram illustrating a comparison of single and multiple PET image acquisition processes according to an exemplary embodiment of the present invention.
FIG. 6b is a schematic view of a two component compartmental model according to an exemplary embodiment of the present invention.
Fig. 6c is a schematic diagram illustrating comparison of an activity curve of a region of interest with an observed value and reconstruction of the region of interest according to an exemplary embodiment of the present invention.
Fig. 7 is a schematic structural diagram of an apparatus for reconstructing a medical dynamic image according to another embodiment of the present invention.
Fig. 8 is a block diagram of an electronic device for executing a method of reconstructing a medical dynamic image according to an exemplary embodiment of the present invention.
Detailed Description
Embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While certain embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be construed as limited to the embodiments set forth herein, but rather are provided for a more thorough and complete understanding of the present invention. It should be understood that the drawings and the embodiments of the present invention are illustrative only and are not intended to limit the scope of the present invention.
The term "include" and variations thereof as used herein are intended to be open-ended, i.e., "including but not limited to". The term "according to" is "at least partially according to". The term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one additional embodiment". Relevant definitions for other terms will be given in the following description.
PET (Positron Emission Computed Tomography) is a relatively advanced clinical examination imaging technology in the field of nuclear medicine, and the method is as follows: after the substances (such as glucose, protein, nucleic acid and fatty acid) necessary for biological life metabolism are labeled with short-lived radioactive nuclide and injected into human body, the condition of life metabolic activity can be reflected by the aggregation degree of said substances in metabolism so as to attain the goal of diagnosis.
The radionuclide commonly used in PET is FDG-18 (fluorodeoxyglucose, the fluorine in the molecule being selected from positron-emitting radioisotopes)18F, FDG-18 for short), the mechanism is as follows: the metabolic states of different tissues of a human body are different, glucose metabolism is vigorous and more glucose is accumulated in high-metabolic malignant tumor tissues, and the characteristics can be reflected through images, so that the pathological changes can be diagnosed and analyzed.
For example, PET imaging can be used to assess the activity of tumors by the intensity of uptake of injected radionuclide drugs by the tumor. In this process, in order to reduce the damage to the body of the subject caused by injection, the shorter half-life FDG-18 is often used as a tracer drug for PET imaging.
In the acquisition process of PET images, on one hand, since the tracer generally enters the body of the subject by injection, it takes a certain time for the tracer to diffuse from the injection site to the whole body, and the signal intensity of the acquired images needs to be gradually strengthened after waiting for a certain time after the injection is completed. On the other hand, however, the signal intensity of the acquired image decreases with the decay of the tracer, requiring the subject to acquire the image as soon as possible after the injection is completed.
From the above two points, it can be seen that an optimal time period in which the tracer is stably distributed needs to be found for image signal acquisition. Due to the difference of human body metabolism and drug transport capacity and the difference of drug absorption capacity of the region of interest, it is difficult to ensure that the time window for collection is just the most suitable time window for observation of the region of interest in the actual collection process.
Fig. 6a is a schematic diagram illustrating a comparison of single and multiple PET image acquisition processes according to an exemplary embodiment of the present invention. As can be seen in connection with section a of fig. 6a, it is difficult for a single PET image acquisition to ensure that the time window of acquisition is just the period of most vigorous tracer drug metabolism in the region of interest, especially for patients without prior PET recordings. Therefore, multiple PET image acquisitions of the B-part (i.e. dynamic PET image acquisitions) are often employed to show as many PET images of the drug metabolic processes in the region of interest as possible. However, dynamic PET image acquisition may also suffer from the same problems, e.g. the periods of maximum metabolism of the tracer drug in the region of interest may intersect two consecutive time windows, or span more than two time windows, or only occupy a small part of a certain time window. These conditions all lead to poor visualization of dynamic PET images, i.e. it is difficult to ensure that the time window of acquisition is exactly the period of time in which the tracer drug metabolism in the region of interest is most vigorous.
In order to solve the above problems, the present invention provides a method of reconstructing a medical image.
FIG. 1 is a schematic diagram of an implementation environment provided by an embodiment of the invention. The implementation environment includes a PET acquisition device 110 and a computer device 120.
The computer device 120 may acquire medical motion images from the PET acquisition device 110. For example, the computer device 120 may communicate with the PET acquisition device 110 over a wired or wireless network.
The PET collecting device 110 is used for collecting the signal intensity of the injected radioactive tracer to obtain the PET medical dynamic image of the human body part. In one embodiment, the PET acquisition device 110 acquires the brain of the human body, and a PET medical dynamic image of the brain can be obtained.
The computer device 120 may be a general-purpose computer or a computer device composed of an application-specific integrated circuit, and the like, which is not limited in the embodiment of the present invention. For example, the Computer device 120 may be a mobile terminal device such as a tablet Computer, or may be a Personal Computer (PC), such as a laptop portable Computer and a desktop Computer.
One skilled in the art will appreciate that the number of computer devices 120 described above may be one or more, and that the types may be the same or different. For example, the number of the computer devices 120 may be one, or the number of the computer devices 120 may be several tens or hundreds, or more. The number and types of the computer devices 120 are not limited in the embodiments of the present invention.
In some optional embodiments, the computer device 120 extracts a region of interest of a PET medical motion image from the PET acquisition device 110, the PET medical motion image including images of a human body part taken at a plurality of time intervals, calculates a curve of activity change over time of the radiopharmaceutical injected into the region of interest, determines a time interval of interest from the curve, the time interval of interest including a time at which a maximum value of the activity is located, and reconstructs an image of the region of interest within the time interval of interest from activity change of pixels of the region of interest.
Fig. 2 schematically shows a flowchart of a method for reconstructing a medical dynamic image according to an embodiment of the present invention. The method of fig. 2 is performed by a computing device (e.g., a server or the computer device of fig. 1), but the embodiments of the invention are not limited thereto. The server may be one server, or may be composed of several servers, or may be a virtualization platform, or a cloud computing service center, which is not limited in the embodiment of the present invention. As shown in fig. 2, the method includes the following.
S210: a region of interest of a medical dynamic image is extracted, the medical dynamic image including images of a human body part taken at a plurality of time periods.
In particular, the medical motion image may be a PET medical motion image. For example, a radiotracer drug may be injected into a body part and a PET medical dynamic image of the body part may be taken. For example, embodiments of the present invention may employ FDG-18 as the radiotracer. Because the half-life of the FDG-18 radioactive tracer drug is 6586.2 seconds, and the distribution concentration of the drug is combined, the time period of each image acquisition can be 5min, 6min, 8min, 10min, 12min, 15min and 20min in sequence. The duration and the number of times of the image capturing period can be increased or decreased properly, and 5 to 8 medical images can be captured. It should be understood that the medical dynamic image of the embodiment of the present invention is not limited to the PET medical dynamic image, but may be other nuclear medicine images.
In an embodiment, a user may select a certain region as the interested area on the user interface by using a selection box, for example, the user may select a certain region of interest of one frame of image in the medical dynamic image, and the computer device may extract an image of the corresponding region of interest in each frame of medical image in the medical dynamic image accordingly. For example, the PET acquisition device 110 acquires a brain of a human body to obtain a PET medical dynamic image of the brain, and a certain region on the PET medical dynamic image of the brain is used as a region of interest.
S220: a curve of the activity of the radiopharmaceutical injected in the region of interest as a function of time is acquired.
In one embodiment, a two-component compartmental model (i.e., kinetic two-component model) may be used to calculate a time-dependent curve of the activity of the radiopharmaceutical injected into the region of interest. FIG. 6b is a schematic view of a two component compartmental model according to an exemplary embodiment of the present invention. As shown in fig. 6b, the two-component compartmental model includes a plasma chamber, a central chamber, and a peripheral chamber, which may also be referred to collectively as a peripheral chamber,the peripheral chamber is also referred to as an isolation chamber or a peripheral chamber, wherein A1Represents the central compartment, A2Represents an isolation chamber, AGIRepresenting the concentration of the tracer drug in the plasma compartment, the variance can be calculated from the above model as:
wherein k isaThe rate of uptake of the tracer drug into the central compartment, also known as the first order rate constant, k12And k21The transfer rates, k, between the central and peripheral chambers, respectively10Is the discharge rate of the central chamber, ka、k12、k21And k10Collectively referred to as pharmacokinetic parameters, C1(t) and C2(t) the concentration of the drug is tracked in the central and peripheral chambers, respectively, at time t. Integrating the above equation then has:
wherein, V1And V2The volumes of the central and peripheral chambers, respectively. In pharmacokinetics, the metabolism of a drug in the body proceeds according to first-order kinetic processes, i.e., the drug also has a relatively stable half-life in the body, and thus needs to be multiplied by a decay factor e-βtObtaining:
C1real(T)=C1(T) e-βt
C2real(T)=C2(T) e-βt
wherein β = (ln2)/T1/2，T1/2Half-life of the radiopharmaceutical, C1real(T) drug activity of the central compartment at time T, C2real(T) drug activity in the peripheral compartment at time T, C1(T) and C2(T) the concentration of the drug is tracked in the central and peripheral chambers, respectively, at time T.
In one embodiment, the time-dependent activity profile of the radiopharmaceutical injected into the region of interest is determined by the following equation:
Creal(T)=k1C1real(T)+ k2C2real(T)
wherein, Creal(T) is the activity of the radiopharmaceutical in the region of interest at time T (i.e. the curve of the activity of the radiopharmaceutical injected in the region of interest as a function of time), C1real(T) drug activity of the central compartment at time T, C2real(T) drug activity in the peripheral compartment at time T, k1Is the proportionality coefficient of the central chamber, k2Is the proportionality coefficient of the peripheral chamber, k1And k2Satisfy k1+k2And =1, namely, the proportion of the central chamber and the peripheral chamber in the region of interest is determined by adopting a linear superposition method.
In fact, any region is not an ideal central or peripheral chamber, and the region of interest can be considered as a superposition of the central and peripheral chambers.
Based on the above formula, the pharmacokinetic parameter k is set according to the initial condition (the first set of images acquired in 5 min) by using least squaresa、k12、k21And k10And setting the above-mentioned proportionality coefficient k1And k2Then, according to other images, gradually fitting the images in an iterative manner, and determining the offset step length of the iterative fitting according to the steepest descent method until the images derived according to the set pharmacokinetic parameter values and the proportionality coefficient values gradually approach to the real images of each group of images (for example, the error is less than 2%), then the pharmacokinetic parameter values and the proportionality coefficient values obtained by fitting are considered to be correct,so as to obtain a curve of the activity of the radiopharmaceutical injected in the region of interest as a function of time.
In the embodiment, the proportion of the central chamber and the peripheral chambers in the region of interest is determined by adopting a linear superposition method, so that the two-component chamber model has higher theoretical performance and accuracy.
S230: from the curve, a period of interest is determined, which contains the moment at which the maximum of the activity is located.
In one embodiment, determining the period of interest from the curve comprises: presetting an activity threshold according to the maximum value of the activity; the period of interest is determined from the maximum value of the activity and the activity threshold.
Specifically, after a curve of the activity of the radiopharmaceutical injected into the region of interest over time is fitted, 70% of the maximum value of the activity may be set as an activity threshold, and a time period in the curve where the activity is higher than the activity threshold may be set as a time period of interest, so as to further establish an image of a time period most suitable for observation in the region of interest based on the time period of interest. The specific value of the activity threshold may be other values, which is not limited in the present invention.
In an embodiment, determining the period of interest from the maximum of the activities and the activity threshold comprises: and determining a first time and a second time when the activity value on the curve is equal to the activity threshold value, wherein the first time is less than the second time, and the interested time period is from the first time to the second time.
Specifically, after the activity threshold is set, an equation that the activity value on the curve is equal to the activity threshold is established, a first time and a second time are obtained, and the activity value on the curve at the first time or the second time is the activity threshold. The first time may be less than the second time. When more than two moments satisfy the equation, the first moment may be the minimum value thereof, and the second moment may be the maximum value thereof, so as to ensure that the maximum activity value is within the selected period of interest.
In one embodiment, determining the period of interest from the curve comprises: determining the slope of a rising edge and the slope of a falling edge of the curve; the period of interest is determined based on the maximum value of the activity, the slope of the rising edge and the slope of the falling edge.
Specifically, after fitting a curve of the activity of the radiopharmaceutical injected in the region of interest over time, the time period of interest may be determined according to the morphology of the curve, so as to further establish an image of the time period in the region of interest that is most suitable for observation based on the time period of interest. For example, the period of interest on the curve may be determined by calculating the slope of the rising edge and the slope of the falling edge of the curve, and when the slope is greater than a certain value, a shorter period around the maximum value of the activity in the curve is taken as the period of interest, and when the slope is less than a certain value, a longer period around the maximum value of the activity in the curve is taken as the period of interest, so that the reconstructed image of the region of interest within the period of interest has a higher signal-to-noise ratio.
S240: and reconstructing an image of the region of interest in the interested period according to the activity change of the pixels of the region of interest.
Specifically, according to activity variation of pixels of the interested region in the acquired PET image, the image of the interested region in the interested period is reconstructed. For example, an image of the region of interest over the period of interest may be obtained by separately reconstructing the projection data for each time frame. The single-frame reconstruction method may employ an analytical reconstruction method (such as a conventional FBP algorithm) or a statistical iterative reconstruction method, such as a Maximum likelihood estimation-based expectation Maximization (ML-EM) algorithm. The time information can also be introduced into the image reconstruction process in the form of a time basis function, so that the time activity curve of the dynamic image is smoother, and the signal-to-noise ratio of the time activity curve can be improved. In addition, a direct parametric imaging method can also be adopted, and the image of the interested region in the interested period can be directly reconstructed from the projection data by combining the dynamic parameter estimation and the dynamic image reconstruction.
Fig. 6c is a schematic diagram illustrating comparison of an activity curve of a region of interest with an observed value and reconstruction of the region of interest according to an exemplary embodiment of the present invention. As shown in fig. 6c, the fitted curve is very close to the actual observed value, which provides a more intuitive and accurate analysis method for the reader to evaluate the distribution change rule and the transmission capability of the tracer drugs.
In the aspect of engineering implementation, java is used as a development language, for a model fitting part, mathmatic in wolfram is used for implementation, and a model fitting algorithm in wolfram is integrated into a java product by using J/link in workbench in wolfram.
According to the technical scheme of the invention, the dynamic image is analyzed, and the image of the region of interest in the time period most suitable for observation is reconstructed, so that a reader can visually see the signal intensity change condition of the region of interest.
Fig. 3 is a schematic flowchart of a process for calculating an activity curve of a region of interest according to another embodiment of the present invention, including:
s310: and establishing a kinetic equation based on the atrioventricular model, wherein the atrioventricular model comprises a plasma chamber, a central chamber and a peripheral chamber, and the kinetic equation is established based on the atrioventricular model of the region of interest, the atrioventricular model comprises the central chamber and the peripheral chamber, and the kinetic equation comprises an equation of the drug activity of the central chamber changing along with time and an equation of the drug activity of the peripheral chamber changing along with time.
S320: a proportionality coefficient is determined between the equation of time-varying drug activity in the central chamber and the equation of time-varying drug activity in the peripheral chambers.
S330: and determining a curve according to a kinetic equation and a proportionality coefficient.
For specific contents of S310 to S320, reference may be made to the description in the above embodiments, and details are not repeated here to avoid repetition.
Fig. 4 is a schematic flowchart of determining a period of interest according to another embodiment of the present invention, which includes:
s410: the activity threshold is preset according to the maximum value of the activity.
S420: the period of interest is determined from the maximum value of the activity and the activity threshold.
For specific contents of S410 to S420, reference may be made to the description in the above embodiments, and details are not repeated here to avoid repetition.
Fig. 5 is a schematic flowchart of determining a period of interest according to another embodiment of the present invention, which includes:
s510: the slope of the rising edge and the slope of the falling edge of the curve are determined.
S520: the period of interest is determined based on the maximum value of the activity, the slope of the rising edge and the slope of the falling edge.
For specific contents of S510 to S520, reference may be made to the description in the above embodiments, and details are not repeated here to avoid repetition.
Fig. 7 is a schematic structural diagram of an apparatus for reconstructing a medical dynamic image according to another embodiment of the present invention, including:
an extracting module 710, configured to extract a region of interest of a medical dynamic image, where the medical dynamic image includes images of a human body part captured at multiple time intervals;
an obtaining module 720, configured to obtain a time-varying curve of activity of the radiopharmaceutical injected into the region of interest;
the determining module 730 is configured to determine a time period of interest according to the curve, the time period of interest including a time instant at which the maximum value of the activity is located.
And a reconstructing module 740 configured to reconstruct an image of the region of interest in the interested time period according to the activity change of the pixels of the region of interest.
According to the device for reconstructing the medical dynamic image, the dynamic image is analyzed, and the image of the region of interest in the time period most suitable for observation is reconstructed, so that a viewer can visually see the signal intensity change condition of the region of interest.
According to an embodiment of the present invention, the obtaining module 720 establishes a kinetic equation based on an atrioventricular model of the region of interest, the atrioventricular model including a central chamber and a peripheral chamber, the kinetic equation including an equation of a change in drug activity of the central chamber with time and an equation of a change in drug activity of the peripheral chamber with time; determining a proportionality coefficient between an equation of time-varying drug activity in the central chamber and an equation of time-varying drug activity in the peripheral chambers; and determining a curve according to a kinetic equation and a proportionality coefficient.
According to an embodiment of the present invention, the curve of the obtaining module 720 is determined by the following formulaAnd (3) determining: creal(T)=k1C1real(T)+k2C2real(T) wherein Creal(T) Activity of the region of interest at time T, C1real(T) drug activity of the central compartment at time T, C2real(T) drug activity in the peripheral compartment at time T, k1Is the proportionality coefficient, k, of the equation relating the activity of the drug in the central compartment to the time2Proportionality coefficient, k, of equation of time-dependent change in drug activity in the peripheral compartment1And k2Satisfy k1+k2=1。
According to an embodiment of the present invention, the determining module 730 presets an activity threshold according to a maximum value of the activity; the period of interest is determined from the maximum value of the activity and the activity threshold.
According to an embodiment of the present invention, the determining module 730 determines a first time and a second time on the curve where the activity value is equal to the activity threshold, the first time is less than the second time, and the period of interest is from the first time to the second time.
The determination module 730 determines the slope of the rising edge and the slope of the falling edge of the curve, according to an embodiment of the invention; the period of interest is determined based on the maximum value of the activity, the slope of the rising edge and the slope of the falling edge.
For specific limitations of the apparatus for reconstructing a medical dynamic image, reference may be made to the above limitations of the method for reconstructing a medical dynamic image, which are not described herein again.
Fig. 8 is a block diagram of an electronic device for executing a method for reconstructing a medical dynamic image according to an exemplary embodiment of the present application, which includes a processor 810 and a memory 820.
The memory 820 is used to store processor executable instructions. The processor is used for executing the executable instructions to execute the method for reconstructing the medical dynamic image in any one of the above embodiments.
The present application further provides a computer-readable storage medium, which stores a computer program for executing the method for reconstructing a medical dynamic image according to any one of the above embodiments.
The method for reconstructing the medical dynamic image adopts an improved two-component chamber model method, analyzes the dynamic image which is decayed to a certain degree by injecting the tracer agent to the tracer agent, obtains the curve of the activity change of the radiopharmaceutical injected into the region of interest along with the time, determines the most suitable interesting time period for observation in the region of interest according to the maximum value of the activity of the curve, and reconstructs the image of the most suitable time period for observation in the region of interest, so that a reader can visually see the signal intensity change condition of the region of interest.
In the above embodiments, all or part of the implementation may be realized by software, hardware, firmware or any other combination. When implemented in software, may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. The procedures or functions according to the embodiments of the invention are brought about in whole or in part when the computer program instructions are loaded and executed on a computer. The computer may be a general purpose computer, a special purpose computer, a network of computers, or other programmable device. The computer instructions may be stored on a computer readable storage medium or transmitted from one computer readable storage medium to another computer readable storage medium, for example, the computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wire (e.g., coaxial cable, fiber optic, Digital Subscriber Line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.). The computer-readable storage medium can be any available medium that can be accessed by a computer or a data storage device, such as a server-side, data center, etc., that includes one or more of the available media. The usable medium may be a magnetic medium (e.g., a floppy disk, a hard disk, a magnetic tape), an optical medium (e.g., a Digital Video Disk (DVD)), or a semiconductor medium (e.g., a Solid State Disk (SSD)), among others.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
In the embodiments provided in the present invention, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, a division of a unit is merely a logical division, and an actual implementation may have another division, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
Units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present invention, and the changes or substitutions should be covered within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110953441.3A CN113397581B (en) | 2021-08-19 | 2021-08-19 | Method and device for reconstructing medical dynamic image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110953441.3A CN113397581B (en) | 2021-08-19 | 2021-08-19 | Method and device for reconstructing medical dynamic image |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113397581A true CN113397581A (en) | 2021-09-17 |
| CN113397581B CN113397581B (en) | 2021-11-30 |
Family
ID=77688889
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110953441.3A Active CN113397581B (en) | 2021-08-19 | 2021-08-19 | Method and device for reconstructing medical dynamic image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113397581B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113989231A (en) * | 2021-10-28 | 2022-01-28 | 上海联影医疗科技股份有限公司 | Determination method, apparatus, computer equipment and storage medium of kinetic parameters |
| CN115115826A (en) * | 2022-05-30 | 2022-09-27 | 深圳技术大学 | Feature selection and extraction method and device, anomaly detection model and construction method thereof |
| CN118452959A (en) * | 2024-05-07 | 2024-08-09 | 广州医科大学附属第一医院（广州呼吸中心） | Calculation method of radioactivity distribution in both lungs based on PET/CT dynamic imaging |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6745066B1 (en) * | 2001-11-21 | 2004-06-01 | Koninklijke Philips Electronics, N.V. | Measurements with CT perfusion |
| US20100046820A1 (en) * | 2008-07-22 | 2010-02-25 | Jerome Declerck | Framing of positron emission tomography data to assess activity peak |
| US20110060602A1 (en) * | 2009-09-09 | 2011-03-10 | Grudzinski Joseph J | Treatment Planning System For Radiopharmaceuticals |
| US20110311021A1 (en) * | 2010-06-16 | 2011-12-22 | Shinsuke Tsukagoshi | Medical image display apparatus and x-ray computed tomography apparatus |
| US20130109964A1 (en) * | 2011-10-14 | 2013-05-02 | Matthew David Kelly | Methods and apparatus for analyzing medical imaging data |
| US20160216246A1 (en) * | 2012-06-05 | 2016-07-28 | The University Of Utah Research Foundation | Reduced parameter space kinetic modeling systems and methods |
| US20160232330A1 (en) * | 2013-09-27 | 2016-08-11 | Commonwealth Scientific And Industrial Research Organisation | Manifold Diffusion of Solutions for Kinetic Analysis of Pharmacokinetic Data |
| US20160296178A1 (en) * | 2015-04-07 | 2016-10-13 | Siemens Aktiengesellschaft | Determining an initialization time point of imaging using a contrast medium |
| CN108320793A (en) * | 2018-01-17 | 2018-07-24 | 江苏赛诺格兰医疗科技有限公司 | The method and apparatus for determining scanning duration |
| CN108932741A (en) * | 2018-06-14 | 2018-12-04 | 上海联影医疗科技有限公司 | Dynamic PET parameter imaging method, device, system and computer readable storage medium |
| CN110269590A (en) * | 2019-06-20 | 2019-09-24 | 上海联影医疗科技有限公司 | Pharmacokinetic parameter determines method, apparatus, computer equipment and storage medium |
| CN110312475A (en) * | 2016-11-11 | 2019-10-08 | 奥胡斯大学 | Method and system for modeling the human heart and atrium |
| CN110996800A (en) * | 2018-08-01 | 2020-04-10 | 联影美国公司 | System, method for determining kinetic parameters of PET imaging |
| US20200367846A1 (en) * | 2018-01-08 | 2020-11-26 | The Regents Of The University Of California | Time-varying kinetic modeling of high temporal-resolution dynamic pet data for multiparametric imaging |
| US20210022697A1 (en) * | 2018-02-15 | 2021-01-28 | The Regents Of The University Of California | Hepatic inflammation analysis with dynamic pet |
-
2021
- 2021-08-19 CN CN202110953441.3A patent/CN113397581B/en active Active
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6745066B1 (en) * | 2001-11-21 | 2004-06-01 | Koninklijke Philips Electronics, N.V. | Measurements with CT perfusion |
| US20100046820A1 (en) * | 2008-07-22 | 2010-02-25 | Jerome Declerck | Framing of positron emission tomography data to assess activity peak |
| US20110060602A1 (en) * | 2009-09-09 | 2011-03-10 | Grudzinski Joseph J | Treatment Planning System For Radiopharmaceuticals |
| US20110311021A1 (en) * | 2010-06-16 | 2011-12-22 | Shinsuke Tsukagoshi | Medical image display apparatus and x-ray computed tomography apparatus |
| US20130109964A1 (en) * | 2011-10-14 | 2013-05-02 | Matthew David Kelly | Methods and apparatus for analyzing medical imaging data |
| US20160216246A1 (en) * | 2012-06-05 | 2016-07-28 | The University Of Utah Research Foundation | Reduced parameter space kinetic modeling systems and methods |
| US20160232330A1 (en) * | 2013-09-27 | 2016-08-11 | Commonwealth Scientific And Industrial Research Organisation | Manifold Diffusion of Solutions for Kinetic Analysis of Pharmacokinetic Data |
| US20160296178A1 (en) * | 2015-04-07 | 2016-10-13 | Siemens Aktiengesellschaft | Determining an initialization time point of imaging using a contrast medium |
| CN110312475A (en) * | 2016-11-11 | 2019-10-08 | 奥胡斯大学 | Method and system for modeling the human heart and atrium |
| US20200367846A1 (en) * | 2018-01-08 | 2020-11-26 | The Regents Of The University Of California | Time-varying kinetic modeling of high temporal-resolution dynamic pet data for multiparametric imaging |
| CN108320793A (en) * | 2018-01-17 | 2018-07-24 | 江苏赛诺格兰医疗科技有限公司 | The method and apparatus for determining scanning duration |
| US20210022697A1 (en) * | 2018-02-15 | 2021-01-28 | The Regents Of The University Of California | Hepatic inflammation analysis with dynamic pet |
| CN108932741A (en) * | 2018-06-14 | 2018-12-04 | 上海联影医疗科技有限公司 | Dynamic PET parameter imaging method, device, system and computer readable storage medium |
| CN110996800A (en) * | 2018-08-01 | 2020-04-10 | 联影美国公司 | System, method for determining kinetic parameters of PET imaging |
| CN110269590A (en) * | 2019-06-20 | 2019-09-24 | 上海联影医疗科技有限公司 | Pharmacokinetic parameter determines method, apparatus, computer equipment and storage medium |
Non-Patent Citations (1)
| Title |
|---|
| 张浩: "基于分数阶微积分的药物代谢动力学建模及其分析", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113989231A (en) * | 2021-10-28 | 2022-01-28 | 上海联影医疗科技股份有限公司 | Determination method, apparatus, computer equipment and storage medium of kinetic parameters |
| CN113989231B (en) * | 2021-10-28 | 2025-09-26 | 上海联影医疗科技股份有限公司 | Method, device, computer equipment and storage medium for determining kinetic parameters |
| CN115115826A (en) * | 2022-05-30 | 2022-09-27 | 深圳技术大学 | Feature selection and extraction method and device, anomaly detection model and construction method thereof |
| CN115115826B (en) * | 2022-05-30 | 2024-10-22 | 深圳技术大学 | Feature selection, extraction method and device, anomaly detection model and construction method thereof |
| CN118452959A (en) * | 2024-05-07 | 2024-08-09 | 广州医科大学附属第一医院（广州呼吸中心） | Calculation method of radioactivity distribution in both lungs based on PET/CT dynamic imaging |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113397581B (en) | 2021-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113397581B (en) | Method and device for reconstructing medical dynamic image | |
| CN112770838B (en) | Systems and methods for image enhancement using self-attention deep learning | |
| Zanotti-Fregonara et al. | Comparison of eight methods for the estimation of the image-derived input function in dynamic [18F]-FDG PET human brain studies | |
| US9275451B2 (en) | Method, a system, and an apparatus for using and processing multidimensional data | |
| CN110996800B (en) | System, method, and non-transitory computer readable medium for determining PET imaging kinetic parameters | |
| US20100260402A1 (en) | Image analysis | |
| JP5214624B2 (en) | Image generation based on limited data sets | |
| CN110415310B (en) | Medical scanning imaging method, device, storage medium and computer equipment | |
| US11896417B2 (en) | Time-varying kinetic modeling of high temporal-resolution dynamic pet data for multiparametric imaging | |
| US12100075B2 (en) | Image reconstruction by modeling image formation as one or more neural networks | |
| US20140016850A1 (en) | Voxel-Resolution Myocardial Blood Flow Analysis | |
| US20140121511A1 (en) | Rapid Stress-Rest Cardiac PET Imaging Systems and Methods | |
| Yao et al. | Simplified protocol for whole‐body Patlak parametric imaging with 18F‐FDG PET/CT: Feasibility and error analysis | |
| WO2024109762A1 (en) | Pet parameter determination method and apparatus, and device and storage medium | |
| Marin et al. | Numerical surrogates for human observers in myocardial motion evaluation from SPECT images | |
| US20140133707A1 (en) | Motion information estimation method and image generation apparatus using the same | |
| Lehnert et al. | Large-scale Bayesian spatial-temporal regression with application to cardiac MR-perfusion imaging | |
| US20070165926A1 (en) | Data processing system for compartmental analysis | |
| KR102723565B1 (en) | Method, device, computing device and storage medium for determining blood flow rate | |
| CN106548473B (en) | A kind of method and device constructing phase image | |
| CN103263275B (en) | PET (positron emission tomography) physiological parameter reconstruction method based on H-infinity filtering | |
| CN117425434A (en) | System and method for parametric imaging | |
| US20250322565A1 (en) | Pet parameter determination method and apparatus, and device and storage medium | |
| CN106548464B (en) | A kind of method and device constructing phase image | |
| US20230154067A1 (en) | Output Validation of an Image Reconstruction Algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |