CN110215226B - Image attenuation correction method, image attenuation correction device, computer equipment and storage medium - Google Patents
Image attenuation correction method, image attenuation correction device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN110215226B CN110215226B CN201910452083.0A CN201910452083A CN110215226B CN 110215226 B CN110215226 B CN 110215226B CN 201910452083 A CN201910452083 A CN 201910452083A CN 110215226 B CN110215226 B CN 110215226B
- Authority
- CN
- China
- Prior art keywords
- pet
- data
- curve
- respiratory motion
- image
- 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.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 118
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000033001 locomotion Effects 0.000 claims abstract description 233
- 230000000241 respiratory effect Effects 0.000 claims abstract description 180
- 238000002591 computed tomography Methods 0.000 claims description 29
- 238000004590 computer program Methods 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 8
- 230000001360 synchronised effect Effects 0.000 claims description 7
- 238000003702 image correction Methods 0.000 claims description 4
- 230000029058 respiratory gaseous exchange Effects 0.000 abstract description 62
- 201000010099 disease Diseases 0.000 abstract description 7
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000013507 mapping Methods 0.000 description 5
- 230000004060 metabolic process Effects 0.000 description 4
- 210000000056 organ Anatomy 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 210000002216 heart Anatomy 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 210000004185 liver Anatomy 0.000 description 3
- 210000004072 lung Anatomy 0.000 description 3
- 210000000496 pancreas Anatomy 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 210000003484 anatomy Anatomy 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000035565 breathing frequency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 235000014113 dietary fatty acids Nutrition 0.000 description 1
- 229930195729 fatty acid Natural products 0.000 description 1
- 239000000194 fatty acid Substances 0.000 description 1
- 150000004665 fatty acids Chemical class 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 239000008103 glucose Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 102000039446 nucleic acids Human genes 0.000 description 1
- 108020004707 nucleic acids Proteins 0.000 description 1
- 150000007523 nucleic acids Chemical class 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 210000002345 respiratory system Anatomy 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- 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
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- 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
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
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)
- Apparatus For Radiation Diagnosis (AREA)
- Nuclear Medicine (AREA)
Abstract
The application relates to an image attenuation correction method, an image attenuation correction device, computer equipment and a storage medium. The method comprises the following steps: a corrected image, PET scan data, a first respiratory motion profile, and a second respiratory motion profile are acquired. Fitting the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve. And searching matching frame PET data corresponding to the phase information of the corrected image in the PET scanning data according to the corrected image, the PET scanning data and the respiratory motion fitting curve. And carrying out attenuation correction on the PET data of the matching frame according to the correction image. By adopting the method, the PET attenuation correction reconstructed image without breathing artifact can be obtained. A doctor can comprehensively and accurately judge the diseases through PET attenuation correction reconstructed images without breathing artifacts.
Description
Technical Field
The present application relates to the field of medical image processing technologies, and in particular, to an image attenuation correction method, an image attenuation correction device, a computer device, and a storage medium.
Background
With the continuous development of medical imaging technology, in order to better inspect the human body, a plurality of technology fusion modes are adopted to inspect the human body. For example, PET-CT uses PET (Positron Emission Computed Tomography, positron emission tomography) for organ and soft tissue detection, and CT (Computed Tomography, electronic computer tomography) for body layer detection of the human body. By simultaneously obtaining the CT image and the PET image, the advantages of the two images are complementary, so that doctors can obtain accurate anatomical positioning while knowing biological metabolism information, and thus, comprehensive and accurate judgment can be made on diseases.
However, at present, when a patient is detected by a PET-CT system, CT images are imaged corresponding to a single frame or adjacent frames of human respiratory motion due to the rapid CT scanning speed. Whereas PET images are typically scanned for a relatively long period of time, so that the PET images correspond to average imaging of human breath. Therefore, when detecting a patient, the patient respiratory motion is large, which causes a large difference between the CT image and the PET image. When the CT image is used for carrying out attenuation correction on the PET image, respiratory artifacts exist in the reconstructed PET image, so that the judgment of diseases is affected.
Disclosure of Invention
In view of the foregoing, it is desirable to provide an image attenuation correction method, apparatus, computer device, and storage medium that can reduce PET image breathing artifacts after PET image reconstruction using CT images.
A method of image attenuation correction, the method comprising:
acquiring a correction image and PET scanning data;
acquiring a first respiratory motion curve and a second respiratory motion curve;
fitting the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve;
searching matching frame PET data corresponding to the phase information of the corrected image in the PET scanning data according to the corrected image, the PET scanning data and the respiratory motion fitting curve;
and carrying out attenuation correction on the PET data of the matching frame according to the correction image.
In one embodiment, the acquiring a first respiratory motion profile comprises:
acquiring video data when scanning the corrected image and PET scanning data;
searching chest position video data in the video data;
performing motion recognition on the chest position video data to obtain a first motion curve;
and filtering the first motion curve to obtain the first respiratory motion curve.
In one embodiment, the acquiring the second respiratory motion profile comprises:
acquiring PET scanning data;
performing motion recognition on the PET scanning data to obtain a second motion curve;
and filtering the second motion curve to obtain the second respiratory motion curve.
In one embodiment, the performing motion recognition on the PET scan data to obtain a second motion profile includes: the PET scan data includes time-of-flight information;
determining photon annihilation point positions according to the flight time information, wherein the photon annihilation point positions are used for representing three-dimensional space coordinates of annihilation points;
and performing time dimension cutting on the photon annihilation point position to obtain a second motion curve.
In one of the embodiments of the present application,
the respiratory motion fitting curve is as follows: the first breathing motion fitting curve and the second breathing motion fitting curve which are identical in phase and synchronous in time sequence.
In one embodiment, the searching the PET scan data corresponding to the corrected image phase information in the PET scan data according to the corrected image, the PET scan data and the respiratory motion fitting curve, and using the PET scan data as matching frame PET data includes:
acquiring correction data of a correction image;
according to the correction data and the respiratory motion fitting curve, obtaining phase information of the correction data in the respiratory motion fitting curve;
and searching PET scanning data matched with the phase information in the PET scanning data according to the phase information, the respiratory motion fitting curve and the PET scanning data, and taking the PET scanning data as matched frame PET data.
In one embodiment, the obtaining the phase information of the correction data in the respiratory motion fitting curve according to the correction data and the respiratory motion fitting curve includes:
and obtaining phase information of the correction data in the first respiratory motion fitting curve according to the correction data and the first respiratory motion fitting curve.
In one embodiment, the searching the PET scan data matched with the phase information in the PET scan data according to the phase information, the respiratory motion fitting curve and the PET scan data, and using the PET scan data as the matched frame PET data includes:
searching a plurality of time point information corresponding to the phase information in a first respiratory motion fitting curve according to the phase information and the first respiratory motion fitting curve;
and according to the time point information and a second breathing fitting curve, searching PET scanning data matched with the time point information in the second breathing motion fitting curve, and taking the PET scanning data as matched frame PET data.
An image attenuation correction device, the device comprising:
the first acquisition module is used for acquiring a correction image and PET scanning data;
the second acquisition module is used for acquiring the first respiratory motion curve and the second respiratory motion curve;
the curve fitting module is used for fitting the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve;
the matching frame PET data obtaining module is used for searching matching frame PET data corresponding to the phase information of the correction image in the PET scanning data according to the correction image, the PET scanning data and the respiratory motion fitting curve;
and the image correction module is used for carrying out attenuation correction on the PET data of the matching frame according to the correction image.
A computer device comprising a memory storing a computer program and a processor which when executing the computer program performs the steps of:
acquiring a correction image and PET scanning data;
acquiring a first respiratory motion curve and a second respiratory motion curve;
fitting the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve;
searching matching frame PET data corresponding to the phase information of the corrected image in the PET scanning data according to the corrected image, the PET scanning data and the respiratory motion fitting curve;
and carrying out attenuation correction on the PET data of the matching frame according to the correction image.
A computer readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of:
acquiring a correction image and PET scanning data;
acquiring a first respiratory motion curve and a second respiratory motion curve;
fitting the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve;
searching matching frame PET data corresponding to the phase information of the corrected image in the PET scanning data according to the corrected image, the PET scanning data and the respiratory motion fitting curve;
and carrying out attenuation correction on the PET data of the matching frame according to the correction image.
The image attenuation correction method, the image attenuation correction device, the computer equipment and the storage medium are characterized in that a respiratory motion fitting curve is obtained by acquiring a correction image, PET scanning data, a first respiratory motion curve and a second respiratory motion curve and fitting the first respiratory motion curve and the second respiratory motion curve, and then the correction image, the PET scanning data and the respiratory motion fitting curve are used for searching matching frame PET data corresponding to the phase information of the correction image in the PET scanning data. And finally, carrying out attenuation correction on the PET data of the matching frame according to the correction image. And mapping the corrected image into a respiratory motion fitting curve, and searching matching frame PET data corresponding to the corrected image phase information in the PET scanning data. And carrying out attenuation correction on the PET scanning data of the matching frame according to the correction image, so as to obtain a PET attenuation correction reconstructed image without breathing artifacts.
Drawings
FIG. 1 is a diagram of an application environment for an image attenuation correction method in one embodiment;
FIG. 2 is a schematic diagram of a respiratory motion profile in one embodiment;
FIG. 3 is a block diagram of an image attenuation correction device in one embodiment;
fig. 4 is an internal structural diagram of a computer device in one embodiment.
Detailed Description
The present application will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present application more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the application.
Positron emission computed tomography (Positron Emission Computed Tomography, PET) is a relatively advanced clinical examination imaging technique in the field of nuclear medicine. Is to mix a substance, which is generally necessary for the metabolism of living things, such as: glucose, protein, nucleic acid, fatty acid, short-lived radionuclides (e.g., 18F,11C, etc.) are labeled, and after injection into the human body, the radionuclides release positrons during decay, one positron encounters an electron after travelling a few tenths of a millimeter to a few millimeters and annihilates, thereby generating a pair of photons with energy of 511KeV in opposite directions. This pair of photons is captured by a highly sensitive camera and corrected for scatter and random information by a computer. By carrying out the same analysis processing on different positrons, we can obtain a three-dimensional image of the aggregation situation in the organism, thereby achieving the purpose of diagnosis.
An electronic computed tomography (Computed Tomography, CT) is a cross-sectional scan of a part of the human body, one after the other, with a detector of extremely high sensitivity, using precisely collimated X-ray beams, gamma rays, ultrasound, etc. According to the difference of the absorption and transmittance of different tissues of the human body to the X-ray beam, the data obtained by measuring the human body through the detector is input into an electronic computer, and after the electronic computer processes the data, the cross section or the three-dimensional image of the inspected part of the human body can be photographed, so that the tiny lesions of any part in the human body can be found. The method has the characteristics of quick scanning time, clear images and the like, and can be used for checking various diseases. The rays used can be classified differently according to the type: x-ray CT (X-CT), ultrasonic CT (UCT), gamma-ray CT (gamma-CT), and the like.
PET-CT integrates PET and CT on an instrument to form a complete imaging system, which is called as a PET-CT system (integrated PET-CT system), and a patient can obtain CT human anatomy images and PET human tissue functional metabolism images simultaneously through rapid whole-body scanning during examination, and the advantages of the two images are complementary, so that a doctor can obtain accurate anatomy positioning while knowing biological metabolism information, and thus, the disease can be comprehensively and accurately judged.
In the prior art, due to the difference of PET and CT scanning speeds, the PET-CT equipment has obvious limitation in the application of chest and abdomen tumor diagnosis and treatment scheme formulation, which mainly means that the PET image appears artifact caused by the respiration of a patient. In the PET-CT system, CT scanning is fast, and only a fraction of a second is needed at the fastest, and the whole body scanning of popular CT is also only a few seconds, and the obtained image is almost a snapshot at a certain moment, so that the image is not influenced by respiration basically; whereas PET scanning is very slow, typically 2D PET acquisition takes 10 to 15 minutes, and 3D PET acquisition takes as much as 2 minutes, obtaining an average image over several tens of respiratory cycles, during which the patient's breathing attracts different degrees of movements of the heart, lung, liver, pancreas etc. organs, causing the radioactive radiation source injected therein to follow the movements as well, so that the PET image inevitably generates movement artefacts. Another source of PET image respiratory motion artifacts is CT attenuation correction artifacts, which are accomplished in PET-CT systems by tissue density images obtained by CT scanning, where the PET and CT images of the same patient are not perfectly matched in spatial position and phase, and an instantaneous CT image is used to attenuate the average PET image to account for errors that necessarily occur in the PET image.
In one embodiment, as shown in fig. 1, there is provided an image attenuation correction method including the steps of:
and step 101, acquiring corrected images and PET scanning data.
In this embodiment, a PET-CT system acquires a plurality of frames of original images and PET scan data, where the PET scan data is original PET data obtained by gating and scanning a PET object to be measured, and the original corrected image is an original corrected CT image obtained by scanning a CT object to be scanned, and both images are obtained by scanning the same position of the PET object to be scanned by the PET-CT system. When a PET-CT system is scanning a human body, the scanning is not terminated instantaneously, but continues for a period of time during which the human body is accompanied by respiratory motion. At this time, the CT scan and the PET scan are not performed simultaneously, and the CT scan is performed rapidly to obtain an original corrected CT image containing no respiratory motion information, and then the PET scan is performed for a period of time to obtain PET scan data containing respiratory motion information of several tens of respiratory cycles.
Step 102, a first respiratory motion profile and a second respiratory motion profile are acquired.
In this embodiment, the first respiratory motion profile is generated from a video signal acquired by a camera. The camera is used for shooting the whole course in the PET-CT scanning process of the patient, and the video signal of the camera contains the breathing motion information of the whole course of the patient in the PET-CT scanning process. The first respiratory motion profile is thus the complete respiratory motion profile of the patient during the PET-CT scan. The second respiratory motion profile is the respiratory motion profile of the patient during the PET scan.
Specifically, in step 102, acquiring a first respiratory motion profile includes: and acquiring video data during scanning the correction image and PET scanning data, searching chest position video data in the video data, performing motion recognition on the chest position video data to obtain a first motion curve, and filtering the first motion curve to obtain the first respiratory motion curve.
In this embodiment, during PET-CT scanning of a patient, the patient is imaged by a camera at the same time, respiratory motion information of the patient during the whole scanning process is recorded, and relevant video data is acquired. The respiratory information in the video data generally includes motion information that causes organs such as heart, lung, liver, pancreas, etc. to move to different extents when the patient breathes. Generally, the respiratory amplitude is obvious at the chest position near the diaphragm, so that chest position video data related to the chest position is extracted from the acquired video data, and motion recognition is performed on the chest position video data to obtain a first motion curve. The acquired first motion profile is a motion profile of the patient throughout the process of performing a PET-CT scan. At this time, the first motion profile also includes noise or a motion signal not generated by respiratory motion of the patient. It is necessary to obtain the first respiratory motion profile by respiratory band filtering the first motion profile, preserving the signal in the respiratory band.
As shown in fig. 2, the first respiratory motion profile and the second respiratory motion profile each include an X-axis representing time and a Y-axis representing respiratory amplitude. In performing a PET-CT scan, the patient is required to undergo a plurality of respiratory motion cycles, one motion respiratory cycle comprising an inhalation phase and an exhalation phase, the highest point of the respiratory motion curve corresponding to the extreme end of the inhalation phase. The respiratory motion curve can be divided into a plurality of respiratory phases in the X-axis direction.
Specifically, the first motion profile may be obtained by motion recognition of chest position video data and AI techniques.
Specifically, in step 102, acquiring a second respiratory motion profile includes: and acquiring PET scanning data, and performing motion recognition on the PET scanning data to obtain a second motion curve. And filtering the second motion curve to obtain the second respiratory motion curve.
In this embodiment, the PET scan data is raw data acquired by a PET scan. The raw data includes time-of-flight data. The time-of-flight data is the time at which a positron emitted by a radionuclide during decay, upon encountering an electron in the patient, annihilates, producing a pair of oppositely directed photons, each captured by a highly sensitive camera. The positions of the pair of photon annihilation points can be determined from the flight data. Because the respiration of the patient during the scanning causes different degrees of movement of organs such as heart, lung, liver, pancreas and the like, the radionuclide injected into the respiratory tract also follows the movement, and the position of the photon annihilation point correspondingly changes in each respiratory phase. And then, the photon annihilation point positions are subjected to time dimension cutting, namely, the annihilation point positions corresponding to each breathing phase are counted, so that a second motion curve is obtained. And similarly, carrying out breathing frequency band filtering on the second motion curve to obtain a second breathing motion curve.
In this embodiment, the position of the photon annihilation point is used to represent three-dimensional space coordinates of the annihilation point in three-dimensional space with the focus of the patient as the origin, that is, the position of the center of gravity of the photon annihilation point.
And step 103, fitting the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve.
In this embodiment, the respiratory motion fitting curve is obtained by fitting the first respiratory motion curve and the second respiratory motion curve. From the breathing motion fitting curve, the corresponding relation of each breathing amplitude or breathing phase of the first breathing motion curve and the second breathing motion curve can be clearly obtained.
In this embodiment, the respiratory motion fitting curve is: the first breathing motion fitting curve and the second breathing motion fitting curve which are identical in phase and synchronous in time sequence.
In other embodiments, the respiratory motion fitting curve is: the first breathing motion fitting curve and the second breathing motion fitting curve which are identical in phase, synchronous in time sequence and consistent in breathing amplitude.
Step 104, searching matching frame PET data corresponding to the phase information of the corrected image in the PET scanning data according to the corrected image, the PET scanning data and the respiratory motion fitting curve.
In this embodiment, by mapping the corrected image to the respiratory motion fitting curve, the respiratory amplitude or the respiratory phase corresponding to the corrected image is searched in the respiratory motion fitting curve, and then the matching frame PET data equivalent to the corrected image information can be found correspondingly through the respiratory amplitude or the respiratory phase.
In step 104, correction data of a correction image is obtained, and phase information of the correction data in the respiratory motion fitting curve is obtained according to the correction data and the respiratory motion fitting curve. And searching matching frame PET data matched with the phase information in the PET scanning data according to the phase information, the respiratory motion fitting curve and the PET scanning data.
Wherein, according to the correction data and the respiratory motion fitting curve, obtaining the phase information of the correction data in the respiratory motion fitting curve includes: according to the correction data and the first respiratory motion fitting curve, phase information of the correction data in the first respiratory motion fitting curve is obtained, and then according to the phase information and the first respiratory motion fitting curve, a plurality of time point information corresponding to the phase information is searched in the first respiratory motion fitting curve. And searching matching frame PET data matched with the time point information in the second respiratory motion fitting curve according to the time point information and the second respiratory motion fitting curve.
In this embodiment, the correction data of the corrected image is data acquired by CT scanning, and image reconstruction is performed using the correction data. And acquiring the starting time of CT scanning through correction data, and searching a respiratory phase corresponding to the corresponding starting time on a respiratory motion fitting curve according to the starting time. The respiratory motion fitting curve is: the first breathing motion fitting curve and the second breathing motion fitting curve which are identical in phase and synchronous in time sequence.
Specifically, according to the known starting time of the CT scanning, searching a corresponding breathing time phase in a CT scanning period in the first breathing motion fitting curve, and searching a plurality of corresponding time points corresponding to the breathing time phase in a PET scanning period in the first breathing motion fitting curve through the breathing time phase. And obtaining a plurality of matching frame PET data matched with the breathing time phase of the correction image according to the obtained mapping relation between the breathing time points and the second breathing motion fitting curve.
In other embodiments, the correction data of the corrected image is data acquired by a CT scan, and the image reconstruction is performed using the correction data. And acquiring the starting time of CT scanning through correction data, and acquiring corresponding respiratory amplitude according to the starting time. And searching corresponding breathing amplitude on a breathing motion fitting curve according to the breathing amplitude. The respiratory motion fitting curve at this time is: the first breathing motion fitting curve and the second breathing motion fitting curve which are identical in phase and consistent in breathing amplitude.
Specifically, according to the respiration amplitude corresponding to the known starting time of the CT scanning, searching the corresponding respiration amplitude in the CT scanning period in the first respiration motion fitting curve, and searching a plurality of corresponding respiration amplitudes in the PET scanning period in the first respiration motion fitting curve through the respiration amplitude. And obtaining a plurality of matching frame PET data matched with the breathing amplitude of the correction image according to the obtained mapping relation between the breathing amplitudes and the second breathing motion fitting curve.
And 105, performing attenuation correction on the PET data of the matching frame according to the correction image.
In this embodiment, the matching frame PET data is an original image in the same respiratory phase or the same respiratory amplitude as the corrected image, and the PET attenuation correction reconstructed image without respiratory artifacts is obtained by performing attenuation correction on the matching frame PET data and the corrected image. A doctor can comprehensively and accurately judge the diseases through PET attenuation correction reconstructed images without breathing artifacts.
In the image attenuation correction method, the correction image, the PET scanning data, the first breathing motion curve and the second breathing motion curve are obtained, the first breathing motion curve and the second breathing motion curve are fitted to obtain a breathing motion fitting curve, and then matching frame PET data corresponding to the phase information of the correction image is searched in the PET scanning data according to the correction image, the PET scanning data and the breathing motion fitting curve. And finally, carrying out attenuation correction on the matched frame image according to the correction image. Matching frame PET data consistent with the breathing phase of the corrected image is found by mapping the corrected image into a breathing motion fitting curve and using the first breathing motion fitting curve. And carrying out attenuation correction on the PET data of the matching frame according to the corrected image, thereby obtaining a PET attenuation correction reconstructed image without breathing artifacts. A doctor can comprehensively and accurately judge the diseases through PET attenuation correction reconstructed images without breathing artifacts.
It should be understood that, although the steps in the flowchart of fig. 1 are shown in sequence as indicated by the arrows, the steps are not necessarily performed in sequence as indicated by the arrows. The steps are not strictly limited to the order of execution unless explicitly recited herein, and the steps may be executed in other orders. Moreover, at least some of the steps in fig. 1 may include multiple sub-steps or stages that are not necessarily performed at the same time, but may be performed at different times, nor do the order in which the sub-steps or stages are performed necessarily performed in sequence, but may be performed alternately or alternately with at least a portion of other steps or sub-steps of other steps.
In one embodiment, as shown in fig. 3, there is provided an image attenuation correction apparatus including: a first acquisition module 201, a second acquisition module 202, a curve fitting module 203, a matching frame PET data obtaining module 204, and an image correction module 205, wherein:
a first acquisition module 201 is used for acquiring corrected images and PET scan data.
A second acquisition module 202 is configured to acquire the first respiratory motion curve and the second respiratory motion curve.
And the curve fitting module 203 is configured to fit the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve.
And the matching frame PET data obtaining module 204 is configured to search matching frame PET data corresponding to the phase information of the corrected image in the multi-frame original image according to the corrected image, PET scan data and the respiratory motion fitting curve.
An image correction module 205, configured to perform attenuation correction on the matching frame PET data according to the correction image.
For specific limitations of the image attenuation correction device, reference may be made to the above limitations of the image attenuation correction method, and no further description is given here. The respective modules in the image attenuation correction apparatus described above may be implemented in whole or in part by software, hardware, and combinations thereof. The above modules may be embedded in hardware or may be independent of a processor in the computer device, or may be stored in software in a memory in the computer device, so that the processor may call and execute operations corresponding to the above modules.
In one embodiment, a computer device is provided, which may be a terminal, and the internal structure of which may be as shown in fig. 4. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to implement a method of image attenuation correction. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, can also be keys, a track ball or a touch pad arranged on the shell of the computer equipment, and can also be an external keyboard, a touch pad or a mouse and the like.
It will be appreciated by persons skilled in the art that the architecture shown in fig. 4 is merely a block diagram of some of the architecture relevant to the present inventive arrangements and is not limiting as to the computer device to which the present inventive arrangements are applicable, and that a particular computer device may include more or fewer components than shown, or may combine some of the components, or have a different arrangement of components.
In one embodiment, a computer device is provided comprising a memory and a processor, the memory having stored therein a computer program, the processor when executing the computer program performing the steps of:
acquiring a correction image and PET scanning data;
acquiring a first respiratory motion curve and a second respiratory motion curve;
fitting the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve;
searching matching frame PET data corresponding to the phase information of the corrected image in the PET scanning data according to the corrected image, the PET scanning data and the respiratory motion fitting curve;
and carrying out attenuation correction on the PET data of the matching frame according to the correction image. In one embodiment, a computer readable storage medium is provided having a computer program stored thereon, which when executed by a processor, performs the steps of:
acquiring a correction image and PET scanning data;
acquiring a first respiratory motion curve and a second respiratory motion curve;
fitting the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve;
searching matching frame PET data corresponding to the phase information of the corrected image in the PET scanning data according to the corrected image, the PET scanning data and the respiratory motion fitting curve;
and carrying out attenuation correction on the PET data of the matching frame according to the correction image. Those skilled in the art will appreciate that implementing all or part of the above described methods may be accomplished by way of a computer program stored on a non-transitory computer readable storage medium, which when executed, may comprise the steps of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in embodiments provided herein may include non-volatile and/or volatile memory. The nonvolatile memory can include Read Only Memory (ROM), programmable ROM (PROM), electrically Programmable ROM (EPROM), electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous Link DRAM (SLDRAM), memory bus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), among others.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The above examples illustrate only a few embodiments of the application, which are described in detail and are not to be construed as limiting the scope of the application. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the application, which are all within the scope of the application. Accordingly, the scope of protection of the present application is to be determined by the appended claims.
Claims (11)
1. A method of image attenuation correction, the method comprising:
acquiring a correction image and PET scanning data; the correction image is an original correction CT image obtained after the CT scans the object to be scanned;
acquiring a first respiratory motion curve and a second respiratory motion curve; the first respiratory motion curve is generated from video signals acquired by a camera, and the second respiratory motion curve is a respiratory motion curve in the PET scanning process;
fitting the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve;
searching matching frame PET data corresponding to the phase information of the corrected image in the PET scanning data according to the corrected image, the PET scanning data and the respiratory motion fitting curve;
and carrying out attenuation correction on the PET data of the matching frame according to the correction image.
2. The method of claim 1, wherein the acquiring a first respiratory motion profile comprises:
acquiring video data when scanning the corrected image and the PET scanning data;
searching chest position video data in the video data;
performing motion recognition on the chest position video data to obtain a first motion curve;
and filtering the first motion curve to obtain the first respiratory motion curve.
3. The method of claim 1, wherein the acquiring a second respiratory motion profile comprises:
acquiring PET scanning data;
performing motion recognition on the PET scanning data to obtain a second motion curve;
and filtering the second motion curve to obtain the second respiratory motion curve.
4. A method according to claim 3, wherein the motion recognition of the PET scan data to obtain a second motion profile comprises: the PET scan data includes time-of-flight information;
determining photon annihilation point positions according to the flight time information, wherein the photon annihilation point positions are used for representing three-dimensional space coordinates of annihilation points;
and performing time dimension cutting on the photon annihilation point position to obtain a second motion curve.
5. The method of claim 4, wherein the step of determining the position of the first electrode is performed,
the respiratory motion fitting curve is as follows: the first respiratory motion fitting curve and the second respiratory motion fitting curve are identical in phase and synchronous in time sequence.
6. The method of claim 5, wherein said locating PET scan data corresponding to the corrected image phase information in the PET scan data as matching frame PET data based on the corrected image, the PET scan data, and the respiratory motion fitting curve comprises:
acquiring correction data of the correction image;
according to the correction data and the respiratory motion fitting curve, obtaining phase information of the correction data in the respiratory motion fitting curve;
and searching PET scanning data matched with the phase information in the PET scanning data according to the phase information, the respiratory motion fitting curve and the PET scanning data, and taking the PET scanning data as matched frame PET data.
7. The method of claim 6, wherein deriving phase information of the correction data in the respiratory motion fitting curve from the correction data and the respiratory motion fitting curve comprises:
and obtaining phase information of the correction data in the first respiratory motion fitting curve according to the correction data and the first respiratory motion fitting curve.
8. The method of claim 6, wherein the searching for PET scan data matching the phase information from the phase information, the respiratory motion fitting curve, and the PET scan data, and using the PET scan data as matching frame PET data comprises:
searching a plurality of time point information corresponding to the phase information in a first respiratory motion fitting curve according to the phase information and the first respiratory motion fitting curve;
and according to the time point information and a second respiratory motion fitting curve, searching PET scanning data matched with the time point information in the second respiratory motion fitting curve, and taking the PET scanning data as matched frame PET data.
9. An image attenuation correction device, the device comprising:
the first acquisition module is used for acquiring a correction image and PET scanning data; the correction image is an original correction CT image obtained after the CT scans the object to be scanned;
the second acquisition module is used for acquiring the first respiratory motion curve and the second respiratory motion curve; the first respiratory motion curve is generated from video signals acquired by a camera, and the second respiratory motion curve is a respiratory motion curve in the PET scanning process;
the curve fitting module is used for fitting the first respiratory motion curve and the second respiratory motion curve to obtain a respiratory motion fitting curve;
the matching frame PET data obtaining module is used for searching matching frame PET data corresponding to the phase information of the correction image in the PET scanning data according to the correction image, the PET scanning data and the respiratory motion fitting curve;
and the image correction module is used for carrying out attenuation correction on the PET data of the matching frame according to the correction image.
10. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor implements the steps of the method of any of claims 1 to 8 when the computer program is executed.
11. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method of any of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910452083.0A CN110215226B (en) | 2019-05-28 | 2019-05-28 | Image attenuation correction method, image attenuation correction device, computer equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910452083.0A CN110215226B (en) | 2019-05-28 | 2019-05-28 | Image attenuation correction method, image attenuation correction device, computer equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110215226A CN110215226A (en) | 2019-09-10 |
| CN110215226B true CN110215226B (en) | 2023-10-03 |
Family
ID=67818198
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910452083.0A Active CN110215226B (en) | 2019-05-28 | 2019-05-28 | Image attenuation correction method, image attenuation correction device, computer equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110215226B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111768463A (en) * | 2020-06-30 | 2020-10-13 | 上海联影医疗科技有限公司 | PET image attenuation correction method and device and computer equipment |
| CN114159085A (en) * | 2021-12-06 | 2022-03-11 | 深圳市联影高端医疗装备创新研究院 | PET image attenuation correction method and device, electronic equipment and storage medium |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102781331A (en) * | 2009-12-10 | 2012-11-14 | 皇家飞利浦电子股份有限公司 | Method and apparatus for using time of flight information to detect and correct for motion in imaging scans |
| CN103282941A (en) * | 2011-01-05 | 2013-09-04 | 皇家飞利浦电子股份有限公司 | Method and apparatus to detect and correct motion in list-ode pet data with a gated signal |
| WO2015040120A1 (en) * | 2013-09-23 | 2015-03-26 | Koninklijke Philips N.V. | Mr-based attenuation correction in pet/mr imaging with dixon pulse sequence |
| CN107095691A (en) * | 2017-07-11 | 2017-08-29 | 上海联影医疗科技有限公司 | A kind of PET imaging methods and system |
| CN107468267A (en) * | 2017-08-11 | 2017-12-15 | 上海联影医疗科技有限公司 | Data processing method and medical imaging devices |
| CN109035355A (en) * | 2017-07-10 | 2018-12-18 | 上海联影医疗科技有限公司 | System and method for PET image reconstruction |
| CN109741284A (en) * | 2018-01-27 | 2019-05-10 | 上海联影医疗科技有限公司 | System and method for correcting the mismatch as caused by respiratory movement in PET imaging |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7787675B2 (en) * | 2005-03-17 | 2010-08-31 | General Electric Company | System, program product, and methods for attenuation correction of emission data on PET/CT and SPEC/CT |
| EP2163201A1 (en) * | 2008-09-15 | 2010-03-17 | Westfälische Wilhelms-Universität Münster | List mode-based respiratory and cardiac gating in positron emission tomography |
| US10690782B2 (en) * | 2017-09-30 | 2020-06-23 | Uih America, Inc. | Systems and methods for positron emission tomography image reconstruction |
-
2019
- 2019-05-28 CN CN201910452083.0A patent/CN110215226B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102781331A (en) * | 2009-12-10 | 2012-11-14 | 皇家飞利浦电子股份有限公司 | Method and apparatus for using time of flight information to detect and correct for motion in imaging scans |
| CN103282941A (en) * | 2011-01-05 | 2013-09-04 | 皇家飞利浦电子股份有限公司 | Method and apparatus to detect and correct motion in list-ode pet data with a gated signal |
| WO2015040120A1 (en) * | 2013-09-23 | 2015-03-26 | Koninklijke Philips N.V. | Mr-based attenuation correction in pet/mr imaging with dixon pulse sequence |
| CN109035355A (en) * | 2017-07-10 | 2018-12-18 | 上海联影医疗科技有限公司 | System and method for PET image reconstruction |
| CN107095691A (en) * | 2017-07-11 | 2017-08-29 | 上海联影医疗科技有限公司 | A kind of PET imaging methods and system |
| CN107468267A (en) * | 2017-08-11 | 2017-12-15 | 上海联影医疗科技有限公司 | Data processing method and medical imaging devices |
| CN109741284A (en) * | 2018-01-27 | 2019-05-10 | 上海联影医疗科技有限公司 | System and method for correcting the mismatch as caused by respiratory movement in PET imaging |
Non-Patent Citations (1)
| Title |
|---|
| PET/CT图像呼吸运动伪影校正研究进展;许全盛等;《中国生物医学工程学报》;20090820(第04期);第573-580页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110215226A (en) | 2019-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2398390B1 (en) | Model-based extension of field-of-view in nuclear imaging | |
| US8107695B2 (en) | Methods and systems for assessing patient movement in diagnostic imaging | |
| US7813783B2 (en) | Methods and systems for attenuation correction in medical imaging | |
| US11633166B2 (en) | Spatial registration of positron emission tomography and computed tomography acquired during respiration | |
| CN102292743B (en) | Nuclear image reconstruction | |
| US20140010428A1 (en) | Method for extraction of a dataset from a medical image dataset and also medical imaging device | |
| CN110960241A (en) | Method and device for determining scanning parameters of medical image scanning and computer equipment | |
| US20080187094A1 (en) | Method and system for performing local tomography | |
| Hu et al. | Design and implementation of automated clinical whole body parametric PET with continuous bed motion | |
| CN110223247B (en) | Image attenuation correction method, device, computer equipment and storage medium | |
| CN110215226B (en) | Image attenuation correction method, image attenuation correction device, computer equipment and storage medium | |
| JP4997575B2 (en) | Attenuation coefficient map creation apparatus, attenuation coefficient map creation method and program | |
| US20030048937A1 (en) | Method of processing visual imagery from a medical imaging device | |
| WO2006006096A1 (en) | Image processing system for the processing of morphological and functional images | |
| CN111544022A (en) | PET attenuation correction method and device and computer equipment | |
| CN110811665A (en) | PET image attenuation correction method, apparatus, computer device and storage medium | |
| CN110507343B (en) | Playground acquisition method, PET image correction device and computer equipment | |
| EP2711738A1 (en) | A method and a device to generate virtual X-ray computed tomographic image data | |
| CN110215203B (en) | Electrocardiosignal acquisition method and device, computer equipment and storage medium | |
| CN111316327B (en) | Attenuation correction of PET data for moving objects | |
| Wang et al. | Motion correction strategies for enhancing whole-body PET imaging | |
| CN110992280B (en) | PET image attenuation correction method, device and computer equipment | |
| CN110443868B (en) | PET image acquisition method, PET image acquisition device, computer equipment and storage medium | |
| Scott-Jackson | Marker-Less Respiratory Gating for PET Imaging with Intelligent Gate Optimisation | |
| US20100074497A1 (en) | Stabilization of imaging methods in medical diagnostics |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co.,Ltd. Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |