CN118071853A

CN118071853A - Medical image processing method, medical image processing device, computer equipment, medium and imaging system

Info

Publication number: CN118071853A
Application number: CN202211458320.2A
Authority: CN
Inventors: 卢一奂; 赵一璋; 叶青; 董筠
Original assignee: Shanghai United Imaging Healthcare Co Ltd
Current assignee: Shanghai United Imaging Healthcare Co Ltd
Priority date: 2022-11-17
Filing date: 2022-11-17
Publication date: 2024-05-24

Abstract

The application relates to a medical image processing method, a medical image processing device, a medical image processing computer device, a medical image processing medium and an imaging system. The method comprises the following steps: acquiring a motion signal related to an object to be detected from an original scanning image set of the object to be detected; separating the respiratory motion signal and the nuclide motion signal in the motion signal to obtain the respiratory motion signal of the object to be detected; carrying out respiratory motion correction on an original scanning image set of the object to be detected according to respiratory motion signals of the object to be detected, and generating corrected scanning images; by the method, respiratory motion signals can be effectively extracted from the original scanning image in the early stage of radionuclide injection, respiratory motion correction is carried out on the original scanning image in the early stage of radionuclide injection, a high-quality scanning image of the radionuclide in unsteady state is obtained, and the scanning image quality in the whole scanning process is improved.

Description

Medical image processing method, medical image processing device, computer equipment, medium and imaging system

Technical Field

The present application relates to the field of nuclear medicine imaging technology, and in particular, to a medical image processing method, apparatus, computer device, medium, and imaging system.

Background

Positron emission computed tomography (Positron Emission Computed Tomography, PET) is a medical imaging technique that reflects the condition of a patient's vital metabolic activity by injecting a radionuclide into the patient and by imaging the metabolic accumulation of the radionuclide in the body. In the process of PET imaging, respiratory motion of a patient can cause motion artifacts in PET images, thereby causing the problem of PET image blurring. In order to solve the problem of motion artifact in the PET image, a respiratory motion correction method for the PET image may be used to eliminate motion artifact to improve the sharpness of the PET image.

In the conventional method, only when the radionuclide in the patient is in a stable state, the respiratory motion condition of the patient can be detected from the PET image, and then the respiratory motion correction is performed on the PET image of the patient based on the detected respiratory motion condition.

Thus, conventional methods fail to correct for respiratory motion in PET images of a patient with the radionuclide in the patient in an unsteady state, i.e., in a motion state.

Disclosure of Invention

In view of the foregoing, it is desirable to provide a medical image processing method, apparatus, computer device, computer readable storage medium, computer program product and PET imaging system that enable respiratory motion correction of PET images at an early stage of radionuclide injection, i.e., in the case where the radionuclide is in motion in the patient.

In a first aspect, the present application provides a medical image processing method. The method comprises the following steps:

Acquiring a motion signal related to an object to be detected from an original scanning image set of the object to be detected; the motion signals comprise respiratory motion signals of the object to be detected and nuclide motion signals of radionuclides in the body of the object to be detected; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments through the scanning equipment when the radionuclide in the object to be detected is in a motion state;

separating the respiratory motion signal and the nuclide motion signal in the motion signal to obtain the respiratory motion signal of the object to be detected;

and carrying out respiratory motion correction on the original scanning image set of the object to be detected according to the respiratory motion signal of the object to be detected, and generating a corrected scanning image.

In one embodiment, acquiring motion signals related to an object to be measured from an original scan image set of the object to be measured includes:

Determining a region of interest from an original scan image set of an object to be measured;

Fitting the time activity curve according to the pixel value of the region of interest to generate a time activity curve corresponding to the original scanning image set; the time activity curve carries a motion signal related to the object to be measured.

In one embodiment, the separation of the respiratory motion signal and the nuclide motion signal in the motion signal to obtain the respiratory motion signal of the object to be detected includes:

Converting the time activity curve from a time domain signal to a frequency domain signal, and generating a frequency curve corresponding to the time activity curve;

Filtering the frequency curve to generate a target frequency curve corresponding to the respiratory motion signal of the object to be detected;

And converting the target frequency curve from the frequency domain signal to the time domain signal to obtain the respiratory motion signal of the object to be detected.

In one embodiment, respiratory motion correction is performed on an original scan image set of a subject to be measured according to respiratory motion signals of the subject to be measured, and a corrected scan image is generated, including:

Determining a gating phase according to a respiratory motion signal of the object to be detected, and performing gating reconstruction on an original scanning image set of the object to be detected based on the gating phase to obtain multi-frame intermediate scanning images;

Acquiring a CT reconstruction image of an object to be detected;

And carrying out respiratory motion correction on the original scanning image set of the object to be detected according to the multi-frame intermediate scanning image and the CT reconstruction image of the object to be detected, and generating a corrected scanning image.

In one embodiment, respiratory motion correction is performed on an original scan image set of a subject to be measured according to a multi-frame intermediate scan image and a CT reconstructed image of the subject to be measured, and a corrected scan image is generated, including:

determining a reference scan image from the multi-frame intermediate scan image;

Determining deformation fields between the reference scan image and each intermediate scan image for each of the plurality of frames of intermediate scan images;

Respectively acting deformation fields between the reference scanning image and each intermediate scanning image on the CT reconstructed images to obtain a plurality of new CT reconstructed images;

And carrying out respiratory motion correction on the original scanning image set of the object to be detected according to the multi-frame new CT reconstructed image, and generating a corrected scanning image.

In one embodiment, the original scan image set includes a plurality of frames of original scan images arranged in time sequence, a gating phase is determined according to a respiratory motion signal of the object to be measured, and gating reconstruction is performed on the original scan image set of the object to be measured based on the gating phase, so as to obtain a plurality of frames of intermediate scan images, including:

Determining a plurality of gating phases according to respiratory motion signals of the object to be detected;

For each gating phase, determining a multi-frame original scanning image corresponding to the gating phase from the original scanning image set of the object to be detected;

And carrying out fusion processing on the multi-frame original scanning images corresponding to the gating phase to obtain an intermediate scanning image corresponding to the gating phase.

In one embodiment, before acquiring the motion signal of the object to be measured from the original scan image set of the object to be measured, the method further includes:

acquiring original scanning data of an object to be detected;

dividing original scanning data of an object to be detected according to preset time intervals to obtain the original scanning data of a plurality of time periods;

Performing image reconstruction on the original scanning data of each time period to obtain an original scanning image corresponding to the time period;

And generating an original scanning image set of the object to be detected based on the original scanning image corresponding to each time period.

In a second aspect, the application also provides a medical image processing device. The device comprises:

The first acquisition module is used for acquiring motion signals related to the object to be detected from the original scanning image set of the object to be detected; the motion signals comprise respiratory motion signals of the object to be detected and nuclide motion signals of radionuclides in the body of the object to be detected; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments through the scanning equipment when the radionuclide in the object to be detected is in a motion state;

The second acquisition module is used for separating the respiratory motion signal and the nuclide motion signal in the motion signal to acquire the respiratory motion signal of the object to be detected;

The first generation module is used for carrying out respiratory motion correction on the original scanning image set of the object to be detected according to the respiratory motion signal of the object to be detected, and generating corrected scanning images.

In a third aspect, the present application also provides a computer device. The computer device comprises a memory storing a computer program and a processor implementing the steps of the medical image processing method in the first aspect when said computer program is executed.

In a fourth aspect, the present application also provides a computer-readable storage medium. The computer readable storage medium has stored thereon a computer program which, when executed by a processor, implements the steps of the medical image processing method in the first aspect.

In a fifth aspect, the application also provides a computer program product comprising a computer program which, when executed by a processor, implements the steps of the medical image processing method in the first aspect.

In a sixth aspect, the present application also provides a PET imaging system comprising a scanning device and a medical image processing apparatus, wherein,

The scanning equipment is used for acquiring an original scanning image set of the object to be detected and sending the original scanning image set to the medical image processing device; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments when the radionuclide in the object to be detected is in a motion state;

Medical image processing means for acquiring motion signals related to the object to be measured from the original scan image set; the motion signals comprise respiratory motion signals of the object to be detected and nuclide motion signals of radionuclides in the body of the object to be detected;

the medical image processing device is also used for separating the respiratory motion signal and the nuclide motion signal in the motion signal to acquire the respiratory motion signal of the object to be detected; and carrying out respiratory motion correction on the original scanning image set of the object to be detected according to the respiratory motion signal of the object to be detected, and generating a corrected scanning image.

The medical image processing method, the medical image processing device, the computer equipment, the storage medium, the computer program product and the PET imaging system acquire a motion signal related to an object to be detected from an original scanning image set of the object to be detected; the motion signals comprise respiratory motion signals of the object to be detected and nuclide motion signals of radionuclides in the body of the object to be detected; then, separating the respiratory motion signal and the nuclide motion signal in the motion signal to obtain the respiratory motion signal of the object to be detected; carrying out respiratory motion correction on an original scanning image set of the object to be detected according to respiratory motion signals of the object to be detected, and generating corrected scanning images; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments through the scanning equipment when the radionuclide in the object to be detected is in a motion state; that is, in the present application, for an original scan image in which a radionuclide acquired at an early stage of radionuclide injection is in a motion state, a respiratory motion signal of a subject to be measured in the motion signal is obtained by acquiring a motion signal related to the subject to be measured from the original scan image, and separating the respiratory motion signal of the subject to be measured from a nuclide motion signal of the radionuclide in the subject to be measured, thereby obtaining a respiratory motion signal of the subject to be measured, and respiratory motion correction is performed on the original scan image based on the respiratory motion signal, to obtain a corrected scan image; the method can effectively extract respiratory motion signals from the original scanning image in the early stage of radionuclide injection, and realize respiratory motion correction on the original scanning image in the early stage of radionuclide injection, so that a high-quality scanning image of radionuclide in unsteady state is obtained, and the scanning image quality in the whole scanning process is improved.

Drawings

FIG. 1 is a diagram of an application environment for a medical image processing method according to one embodiment;

FIG. 2 is a flow chart of a method of medical image processing according to one embodiment;

FIG. 3 is a flow chart of a method of medical image processing according to another embodiment;

FIG. 4 is a flow chart of a method of medical image processing according to another embodiment;

FIG. 5 is a flow chart of a method of medical image processing according to another embodiment;

FIG. 6 is a block diagram showing the structure of a medical image processing apparatus according to one embodiment;

FIG. 7 is an internal block diagram of a computer device in one embodiment;

fig. 8 is a schematic diagram of the structure of a PET imaging system 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.

In general, during PET image scanning, radionuclides such as fluorodeoxyglucose ¹⁸F-FDG、¹³NH₃、⁸² Rb and the like are injected into a patient, and after the radionuclides are injected for a period of time, that is, when the radionuclides are in a stable state, PET images of the patient are acquired; however, when acquiring PET images of the heart, problems arise in that the PET images are blurred due to the presence of respiratory motion, the reconstructed image artifacts and the quantitative results are inaccurate; therefore, there is a need to detect respiratory motion of a patient and eliminate respiratory motion in PET images, resulting in high quality PET reconstructed images.

Conventionally, respiratory motion detection may be performed by an external device, such as ANZAI belt devices, an Electrocardiogram (ECG), etc., and a data-driven (data-driven) method may be used, such as COD (Centroid Of Distribution), intrinsic orthogonal decomposition (Proper orthogonal decomposition, POD), or principal component analysis (PRINCIPAL COMPONENTS ANALYSIS, PCA), etc., without relying on an external device.

When the data driving method is adopted for PET image reconstruction, a breathing change curve is determined through the data driving method, so that breathing phase change is obtained, and then PET image reconstruction is carried out based on the breathing phase change, so that influence caused by breathing motion is eliminated. However, the precondition of the traditional data driving method is that the radionuclide is not changed, namely the radionuclide is in a stable state, so that the respiratory motion correction can be accurately realized.

Therefore, by adopting the data driving method, respiratory motion information cannot be accurately detected in the early stage of radionuclide injection, namely when the radionuclide is still in a motion state, so that accurate respiratory motion correction and PET image reconstruction cannot be realized in the early stage of radionuclide injection.

Based on the above, the embodiment of the application provides an image reconstruction method, which is used for identifying respiratory motion in early stage of radionuclide injection by a data driving method and correcting the respiratory motion in early stage of injection to obtain a clear PET image in early stage of injection.

The following describes a technical scheme related to an embodiment of the present disclosure in conjunction with a scenario in which the embodiment of the present disclosure is applied.

The medical image processing method provided by the embodiment of the application can be applied to an application environment shown in figure 1. The medical image scanning equipment can be equipment for carrying out PET image scanning, and comprises PET single-mode imaging equipment, PET-MR, PET-CT and other mixed-mode imaging equipment and the like; it should be noted that, the medical image processing method may be applied to a medical image scanning apparatus, or may be applied to a processing apparatus or a server that is communicatively connected to the medical image scanning apparatus, which is not particularly limited in the embodiment of the present application.

In one embodiment, as shown in fig. 2, a medical image processing method is provided, and the method is applied to the medical image scanning apparatus in fig. 1 for illustration, and includes the following steps:

step 201, acquiring a motion signal related to the object to be measured from the original scan image set of the object to be measured.

The motion signals comprise respiratory motion signals of the object to be detected and nuclide motion signals of radionuclides in the body of the object to be detected; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments through scanning equipment when the radionuclide in the object to be detected is in a motion state, and the original scanning image set can generate an original dynamic image sequence; in this embodiment, the scanning device may be a medical image scanning device, such as a PET device, a PET-MR device, or a PET-CT device.

Optionally, after the scanning device collects original scanning data in early detection of PET, image reconstruction can be performed on the original scanning data according to a preset duration, so as to obtain a plurality of original scanning images corresponding to different moments, and an original scanning image set of the object to be detected is formed. Illustratively, the medical image scanning apparatus may acquire raw scan data of an object to be measured; dividing original scanning data of an object to be detected according to preset time intervals to obtain the original scanning data of a plurality of time periods; then, carrying out image reconstruction on the original scanning data of each time period aiming at the original scanning data of each time period to obtain an original scanning image corresponding to each time period; and generating an original scanning image set of the object to be detected based on the original scanning image corresponding to each time period.

For example, the medical image scanning apparatus may acquire raw scan data (which may also be referred to as PET raw data) within 3 minutes after the radionuclide is injected into the subject to be detected, and reconstruct an image every 100 milliseconds in the raw scan data to obtain a set of raw scan images with a period of 100 milliseconds, that is, the set of raw scan images includes a raw scan image corresponding to 100 milliseconds, a raw scan image corresponding to 200 milliseconds, a raw scan image corresponding to 300 milliseconds, and a raw scan image corresponding to up to 3 minutes. It should be noted that this example is only for illustrating one example of the generation of the original scan image set, and is not used to limit the original scan image set, for example: in practical use, the reconstructed data of each frame of the original scanned image may be different from 10 ms to 500 ms, and of course, may be a time interval outside this range, for example, 1000 ms, that is, the original scanned image reconstruction may be based on a second level, a sub-second level reconstruction, etc., which is not particularly limited in the embodiment of the present application.

Alternatively, each original scan image in the original scan image set may be a reconstructed image without attenuation correction based on original scan data, and the reconstructed image may be a reconstructed image obtained based on an ordered subset maximum expected value (Ord-ered Subsets Expectation Maximization, abbreviated as OSEM) method; of course, other reconstruction algorithms may be employed, such as: the image reconstruction is performed on the original scan data by a filtered back projection reconstruction algorithm (Filtered Back Projection, FBP), a reconstruction algorithm with a priori, regularization, AI, etc., to obtain an original scan image set, which is not particularly limited in the embodiment of the present application.

Optionally, the medical image scanning device may extract a time activity curve according to an original scan image set of the object to be detected, and further obtain a motion signal related to the object to be detected from the time activity curve; wherein the time activity curve is used to characterize the concentration change of the radionuclide in the subject. Optionally, when the medical image scanning device acquires the time activity curve according to the original scanning image set, the original scanning image set may be input into a preset time activity curve model to output and obtain the time activity curve corresponding to the original scanning image set.

Step 202, separating the respiratory motion signal and the nuclide motion signal in the motion signal to obtain the respiratory motion signal of the object to be detected.

Wherein the motion signal is a mixed motion signal of various different types of motion signals related to the object to be measured, such as: the mixed motion signal may be a mixed motion signal including a respiratory motion signal of the subject and a radionuclide motion signal of the radionuclide in the subject. Of course, the hybrid motion signal may also comprise other types of motion signals, such as: and (3) twisting, overturning and other autonomous motion signals of the object to be measured in a non-stationary state.

Optionally, after the medical image scanning device acquires the motion signal related to the object to be detected, the motion signal can be analyzed according to the characteristics of the respiratory motion signal, and the respiratory motion signal of the object to be detected is separated from the motion signal; optionally, the medical image scanning device may further analyze the motion to separate each different type of motion signal from the motion signals, and then acquire a respiratory motion signal of the object to be detected from each different type of motion signal according to the characteristics of the respiratory motion signal.

Step 203, respiratory motion correction is performed on the original scan image set of the object to be detected according to the respiratory motion signal of the object to be detected, and a corrected scan image is generated.

Optionally, after the medical image scanning device acquires the respiratory motion signal of the object to be detected, a preset respiratory motion correction algorithm may be adopted, and respiratory motion correction processing is performed on the original scanning image in the original scanning image set based on the respiratory motion signal of the object to be detected, so as to generate a corrected scanning image; optionally, the preset respiratory motion correction algorithm may be a respiratory correction algorithm based on data driving, a respiratory correction algorithm based on respiratory gating, a respiratory correction algorithm based on other algorithm principles or algorithm models, and the like, for example: respiratory correction algorithms based on neural network models, deep learning models, machine learning models, etc.; the embodiment of the present application is not particularly limited thereto.

In addition, for the corrected scanned image, it may include corrected scanned images respectively corresponding to the original scanned images in the original scanned image set, forming a corrected dynamic image sequence; the method can also comprise a new dynamic image sequence obtained by carrying out image fusion processing on each corrected scanning image according to a preset fusion strategy based on the corrected scanning image corresponding to each original scanning image; for example: based on the above example, in the case where the original scan image set includes the original scan image corresponding to every 100 milliseconds, after the respiratory motion correction is performed, the corrected scan image corresponding to every 100 milliseconds can be obtained, and then, further, the corrected scan images corresponding to every 100 milliseconds can be subjected to the merging process every 5 seconds, that is, the merging process is performed on the plurality of corrected scan images every 5 seconds, so that the new moving image sequence formed in the time sequence of one frame of 5 seconds can be obtained.

It should be noted that, for the preset fusion strategy, it may be determined according to a specific quantitative analysis performed on the PET image in practical application, for example: based on the new dynamic image sequence of one frame every 5 seconds, calculating the Myocardial Blood Flow (MBF) index of myocardial perfusion through parameter analysis; the embodiment of the present application is not particularly limited thereto. Optionally, corresponding fusion strategies can be preset for different types of quantitative analysis operations to form a preset corresponding relation, when image fusion is performed, a target fusion strategy corresponding to the target quantitative analysis operation can be determined according to the preset corresponding relation according to the target quantitative analysis operation input by a user, and then image fusion processing operation is performed on each corrected original scanning image corresponding to an original scanning image set according to the target fusion strategy to obtain a target dynamic image sequence; the target dynamic image sequence is used for carrying out target quantitative analysis operation to obtain a target quantitative analysis result.

In the medical image processing method, a motion signal related to an object to be detected is acquired from an original scanning image set of the object to be detected; the motion signals comprise respiratory motion signals of the object to be detected and nuclide motion signals of radionuclides in the body of the object to be detected; then, separating the respiratory motion signal and the nuclide motion signal in the motion signal to obtain the respiratory motion signal of the object to be detected; carrying out respiratory motion correction on an original scanning image set of the object to be detected according to respiratory motion signals of the object to be detected, and generating corrected scanning images; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments through the scanning equipment when the radionuclide in the object to be detected is in a motion state; that is, in the present application, for an original scan image in which a radionuclide acquired at an early stage of radionuclide injection is in a motion state, a respiratory motion signal of a subject to be measured in the motion signal is obtained by acquiring a motion signal related to the subject to be measured from the original scan image, and separating the respiratory motion signal of the subject to be measured from a nuclide motion signal of the radionuclide in the subject to be measured, thereby obtaining a respiratory motion signal of the subject to be measured, and respiratory motion correction is performed on the original scan image based on the respiratory motion signal, to obtain a corrected scan image; the method can effectively extract respiratory motion signals from the original scanning image in the early stage of radionuclide injection, and realize respiratory motion correction on the original scanning image in the early stage of radionuclide injection, so that a high-quality scanning image of radionuclide in unsteady state is obtained, and the scanning image quality in the whole scanning process is improved.

Fig. 3 is a flow chart of a method of processing a medical image according to another embodiment. This embodiment relates to an optional implementation process of the medical image scanning apparatus for acquiring a motion signal related to an object to be measured from an original scan image set of the object to be measured, where, based on the above embodiment, as shown in fig. 3, the step 201 includes:

Step 301, determining a region of interest from an original scan image set of an object to be measured.

The region of interest is a region where respiratory motion and nuclide motion occur in the original scanning image set; such as: the region of interest may be the region in which the heart is located, or a region of interest within a certain range around the heart; alternatively, the region of interest may be two-dimensional or three-dimensional, but may also be the entire image.

Optionally, the medical image scanning device may extract a region of interest of at least one original scan image in the original scan image set, and perform fusion processing on the region of interest of the at least one original scan image to obtain a region of interest corresponding to the original scan image set; when extracting the region of interest, the medical image scanning device can adopt a preset region of interest extraction algorithm to determine the region of interest of the original scanned image; optionally, the preset region of interest extraction algorithm may be a preset region of interest automatic sketching algorithm, and the automatic sketching algorithm may be a region of interest sketching algorithm based on a preset model such as neural network, deep learning, etc.; and automatically sketching the region with motion change in the original scanning image by a preset region-of-interest automatic sketching algorithm, so as to determine the region-of-interest corresponding to the original scanning image.

Step 302, fitting a time activity curve according to the pixel values of the region of interest, and generating a time activity curve corresponding to the original scanned image set.

Wherein the time activity curve carries a motion signal related to the object to be measured.

Optionally, for each original scanning image in the original scanning image set, obtaining average pixels of the region of interest of each original scanning image, and fitting a time activity curve according to the average pixels of the region of interest of each original scanning image to generate a time activity curve corresponding to the original scanning image set; the horizontal axis of the time activity curve is a time frame corresponding to each original scanned image in the original scanned image set, for example: 0. 100, 200, 300, etc., in milliseconds, the vertical axis of which is the average pixel of the region of interest, which may also be referred to as the pixel mean, that is, the time activity curve may also be used to characterize the change over time of the pixel mean of the region of interest for each frame of the original scan image in the original scan image set.

In the embodiment, the medical image scanning device determines a region of interest in which respiratory motion and nuclide motion occur from an original scan image set of an object to be detected; fitting the time activity curve according to the pixel value of the region of interest to generate a time activity curve corresponding to the original scanning image set; wherein, the time activity curve carries a motion signal related to the object to be detected; in this embodiment, the motion signal related to the object to be measured is reflected by the pixels of the motion change area, so as to form a time activity curve capable of representing the respiratory motion and the nuclide motion, so that the respiratory motion signal of the object to be measured and the nuclide motion signal of the radionuclide in the body of the object to be measured can be extracted based on the time activity curve, support and basis are provided for the extraction of the respiratory motion signal, and feasibility and accuracy of the respiratory motion signal extraction are improved.

Fig. 4 is a flow chart of a method of processing a medical image according to another embodiment. The embodiment relates to an optional implementation process of the medical image scanning apparatus for separating a respiratory motion signal and a nuclide motion signal in a motion signal to obtain a respiratory motion signal of a subject to be detected, where on the basis of the above embodiment, as shown in fig. 4, the step 202 includes:

step 401, converting the time activity curve from the time domain signal to the frequency domain signal, and generating a frequency curve corresponding to the time activity curve.

The time activity curve obtained above is a change curve of the pixel mean value of the region of interest over time, i.e. a signal change curve in the time domain.

Alternatively, when extracting the motion signal from the time activity curve, the frequency characteristic of the motion signal may be used, i.e. a motion signal having a certain frequency characteristic may be extracted from the frequency signal; that is, the time activity curve can be converted from a time domain signal to a frequency domain signal, so as to obtain a frequency curve corresponding to the time activity curve; different types of motion signals have different frequency characteristics, corresponding to different frequency bands in a frequency curve, i.e. appear as different peaks in the frequency curve. In other words, the respiratory motion signal of the object to be measured and the nuclide motion signal of the radionuclide in the body of the object to be measured correspond to two different frequency bands in the frequency curve corresponding to the time activity curve, and the two different frequency bands correspond to two peaks.

Alternatively, fourier transform processing may be performed on the time activity curve, to convert the time activity curve in the time domain into the time activity curve in the frequency domain, so as to obtain a frequency curve corresponding to the time activity curve.

Step 402, performing filtering processing on the frequency curve to generate a target frequency curve corresponding to the respiratory motion signal of the object to be tested.

Optionally, the frequency band corresponding to the respiratory motion signal may be determined according to the frequency characteristic of the respiratory motion signal, and then a preset band-pass filter is generated based on the frequency band corresponding to the respiratory motion signal, where the preset band-pass filter is used to filter from the frequency curve corresponding to the time activity curve to obtain the frequency band corresponding to the respiratory motion signal, so as to obtain the target frequency curve corresponding to the respiratory motion signal.

Similarly, for the nuclide movement signal, the same mode can be adopted, namely, according to the frequency characteristic of the nuclide movement signal, the frequency curve corresponding to the nuclide movement signal is obtained by filtering from the frequency curve corresponding to the time activity curve.

Step 403, converting the target frequency curve from the frequency domain signal to the time domain signal, to obtain the respiratory motion signal of the object to be detected.

After the target frequency curve corresponding to the respiratory motion signal is obtained, the target frequency curve can be subjected to Fourier inverse transformation, namely, the frequency domain signal is converted into the time domain signal, so that the respiratory motion signal of the object to be detected can be obtained through decomposition from the time activity curve.

Similarly, for a frequency curve corresponding to the nuclide movement signal, after Fourier inverse transformation processing is performed, the nuclide movement signal can be obtained, namely, the nuclide movement signal of the radionuclide in the body of the object to be detected is obtained by decomposing the time activity curve; the nuclear species motion signal may be used for other image processing operations.

In this embodiment, a time activity curve is converted from a time domain signal to a frequency domain signal, and a frequency curve corresponding to the time activity curve is generated; filtering the frequency curve to generate a target frequency curve corresponding to the respiratory motion signal of the object to be detected; converting the target frequency curve from the frequency domain signal to the time domain signal to obtain a respiratory motion signal of the object to be detected; in this embodiment, based on the frequency characteristic of the respiratory motion signal, the time activity curve is converted from the time domain to the frequency domain, so as to obtain a frequency curve corresponding to the time activity curve, then a target frequency curve corresponding to the respiratory motion is obtained by filtering from the frequency curve, and finally the target frequency curve is converted into the time domain curve, so that the respiratory motion signal can be obtained, and the respiratory motion signal of the object to be detected is decomposed from the time activity curve; the functional effect of accurately extracting the respiratory motion signals from various fused motion signals is realized, and the extraction efficiency and the extraction accuracy of the respiratory motion signals can be realized.

Fig. 5 is a flow chart of a method of processing a medical image according to another embodiment. The present embodiment relates to an optional implementation process of the medical image scanning apparatus performing respiratory motion correction on an original scan image set of a to-be-detected object according to a respiratory motion signal of the to-be-detected object to generate a corrected scan image, where on the basis of the above embodiment, as shown in fig. 5, the step 203 includes:

Step 501, determining a gating phase according to a respiratory motion signal of the object to be detected, and performing gating reconstruction on an original scanning image set of the object to be detected based on the gating phase to obtain multi-frame intermediate scanning images.

Wherein the original scanned image set comprises a plurality of frames of original scanned images arranged in time sequence; illustratively, the original scanned image corresponds every 100 milliseconds as in the example above.

Optionally, the medical image scanning device may determine a plurality of gating phases according to respiratory motion signals of the object to be detected; next, for each gating phase, determining a multi-frame original scanning image corresponding to the gating phase from the original scanning image set of the object to be detected; performing fusion processing on a plurality of frames of original scanning images corresponding to the gating phase to obtain an intermediate scanning image corresponding to the gating phase; by adopting the mode, the intermediate scanning images corresponding to the gating phases respectively can be obtained, namely, a multi-frame intermediate scanning image is obtained, wherein the number of the intermediate scanning images is the number of the gating phases.

For example, according to the respiratory motion signal of the object to be detected, the respiratory motion signal in one respiratory cycle may be divided into 6 gating phases, and assuming that one respiratory cycle corresponds to 12 frames of original scan images, then, for the first gating phase, each frame of original scan images of the 2 nd, 14 th, 26 th, 38 th, etc. in the set of original scan images with 12 frames as cycles may be extracted, and fusion processing may be performed on the frames of original scan images, to obtain an intermediate scan image corresponding to the first gating phase. Similarly, for the second gating phase, the original scanning images of frames 4, 16, 28, 40, etc. taking 12 frames as the period in the original scanning image set can be extracted, and the frame original scanning images are fused to obtain the intermediate scanning image corresponding to the second gating phase. By analogy, for each gating phase, an intermediate scan image corresponding to each gating phase may be obtained, which will not be discussed in detail.

Step 502, a CT reconstructed image of an object to be measured is acquired.

The CT reconstruction image is a CT scanning image obtained after CT scanning is carried out on the same scanning part of the object to be detected, and is used for representing the attenuation distribution of the object to be detected under X rays; when the PET image under the gamma rays is subjected to attenuation correction, the attenuation distribution of the object to be measured under the gamma rays can be calculated based on the attenuation distribution of the object to be measured under the X rays by the aid of the CT reconstructed image under the X rays, and further, the PET image of the object to be measured is subjected to attenuation correction according to the attenuation distribution of the object to be measured under the gamma rays, so that the PET scanning image subjected to attenuation correction is obtained.

Optionally, the medical image scanning device may acquire a CT reconstructed image of the object to be measured from the CT scanning device, or may acquire a CT reconstructed image of the object to be measured from the server; of course, when the medical image scanning device is a PET-CT device, the CT reconstruction image of the object to be detected can be obtained according to CT scanning data after the PET-CT device finishes CT scanning of the object to be detected; the embodiment of the application does not limit the acquisition mode of the CT reconstruction image in particular.

Step 503, performing respiratory motion correction on the original scan image set of the object to be tested according to the multi-frame intermediate scan image and the CT reconstructed image of the object to be tested, and generating a corrected scan image.

Optionally, the medical image scanning device may perform deformation processing on CT reconstructed images of the object to be detected according to multiple frames of intermediate scanned images with different gating phases, so as to obtain deformed new CT reconstructed images corresponding to each frame of intermediate scanned image, where the new CT reconstructed images are CT reconstructed images after respiratory motion correction; then, based on the multi-frame new CT reconstructed image, performing image reconstruction with attenuation correction on each original scanning image in the original scanning image set of the object to be detected, and generating a scanning image after respiration correction and attenuation correction.

In an alternative implementation, the medical image scanning apparatus may determine the reference scan image from a plurality of frames of intermediate scan images; determining deformation fields between the reference scan image and each intermediate scan image for each of the plurality of frames of intermediate scan images; then, respectively acting deformation fields between the reference scanning image and each intermediate scanning image on the CT reconstructed images to obtain a plurality of new CT reconstructed images; and further, carrying out respiratory motion correction on the original scanning image set of the object to be detected according to the multi-frame new CT reconstructed image, and generating a corrected scanning image.

Alternatively, the reference scan image may be any one of the intermediate scan images corresponding to one respiratory cycle, for example: an intermediate scan image corresponding to the end of breath, i.e., the last gating phase in a breathing cycle, may be used as the reference scan image; then, for the intermediate scanning images corresponding to the first 5 gating phases respectively, determining deformation fields between each intermediate scanning image and the reference scanning image to obtain deformation fields corresponding to the first 5 gating phases respectively; then, the deformation fields corresponding to the 5 gating phases are respectively acted on the CT reconstructed images to obtain 5 deformed new CT reconstructed images; and combining the 5 new CT reconstructed images and the CT reconstructed images of the undeformed object to be detected, and carrying out attenuation correction on the original scanning image set to obtain scanning images after attenuation correction in every 100 milliseconds, namely an attenuation corrected dynamic image sequence.

It should be noted that, for each frame of original scanned image in a gating phase, a new CT reconstructed image corresponding to the gating phase is adopted to perform attenuation correction on each original scanned image; for example: based on the example in the above step 501, when performing attenuation correction on the original scan image set, for a plurality of frames of original scan images in the 1 st gating phase, including the original scan images in the 1 st gating phase in each period of 1 st, 2, 13, 14, 25, 26, 37, 38 frames, etc., the new CT reconstructed image after respiratory motion correction corresponding to the 1 st gating phase is used for attenuation correction, so as to obtain a corrected scan image; furthermore, for the multi-frame original scanning image in the 6 th gating phase, including the original scanning image in the 6 th gating phase in each period of 11 th, 12 th, 23 th, 24 th, 35 th, 36 th, 47 th, 48 th and the like, the original CT reconstructed image corresponding to the 6 th gating phase is adopted for attenuation correction, so as to obtain a corrected scanning image. In addition, for the multiple frame original scan image in the 2 nd gating phase, the multiple frame original scan image in the 3 rd gating phase, the multiple frame original scan image in the 4 th gating phase, and the multiple frame original scan image in the 5 th gating phase, the attenuation correction method of the multiple frame original scan image in the 1 st gating phase may be referred to, so that attenuation correction is performed on each original scan image to obtain a corresponding corrected scan image, which is not discussed one by one here.

Further, after respiratory motion correction and attenuation correction are performed to obtain a corrected dynamic image sequence (denoted as a first dynamic image sequence), respiratory motion image registration may also be performed on the corrected dynamic image sequence (i.e., the first dynamic image sequence) to obtain a final respiratory motion corrected dynamic image sequence (denoted as a second dynamic image sequence); for example, for a first dynamic image sequence, each frame of scanned image in the first dynamic image sequence may be registered to a reference image for motion correction to obtain a second dynamic image sequence; the reference image may be any frame of scanned image in the first dynamic image sequence, or may be a reference image obtained by performing other processing operations such as fusion and merging based on at least one frame of scanned image in the first dynamic image sequence, which is not particularly limited in the embodiment of the present application.

In the embodiment, a gating phase is determined according to a respiratory motion signal of an object to be detected, and gating reconstruction is performed on an original scanning image set of the object to be detected based on the gating phase to obtain multi-frame intermediate scanning images; acquiring a CT reconstruction image of an object to be detected; according to the multi-frame intermediate scanning image and the CT reconstruction image of the object to be detected, respiratory motion correction is carried out on an original scanning image set of the object to be detected, and a corrected scanning image is generated; that is, in this embodiment, gating framing is performed based on the acquired respiratory motion signal of the object to be detected, so that the original scan images with the same phase are fused, that is, the original scan images with consistent respiratory motion are fused, and further, based on the fused intermediate scan image, the CT reconstructed image is subjected to motion correction, so that based on the CT reconstructed image after motion correction, each original scan image is subjected to attenuation correction, and finally, a corrected scan image is obtained; by adopting the method in the embodiment, the respiratory motion correction and the attenuation correction can be effectively carried out on the original scanning image acquired in the early stage of radionuclide injection, namely the radionuclide is still in a motion state, so that the scanning reconstruction image is obtained, and the quality of the scanning image acquired in the early stage of radionuclide injection is improved.

In one embodiment, a complete process flow of a medical image processing method is provided, which may include:

step 1, acquiring original scanning data of an object to be detected in a period of time after radionuclide injection.

And step 2, dividing the original scanning data of the object to be detected according to a preset time interval to obtain the original scanning data of a plurality of time periods.

And 3, reconstructing an OSEM image without attenuation correction on the original scanning data of each time period according to the original scanning data of each time period, and obtaining an original scanning image corresponding to the time period.

Step 4, generating an original scanning image set of the object to be detected based on the original scanning image corresponding to each time period; the original scan image set may also be used as an original dynamic image.

And 5, automatically sketching a region of interest containing respiratory motion information from an original scanning image set of the object to be tested.

Step 6, fitting a time activity curve according to the pixel value of the region of interest, and generating a time activity curve corresponding to the original scanning image set; the time activity curve carries a motion signal related to the object to be measured.

And 7, converting the time activity curve from the time domain signal to the frequency domain signal, and generating a frequency curve corresponding to the time activity curve.

Step 8, filtering the frequency curve by adopting a preset band-pass filter to generate a target frequency curve corresponding to the respiratory motion signal of the object to be detected; wherein the predetermined band-pass filter is related to a frequency characteristic of the respiratory motion signal.

And 9, converting the target frequency curve from the frequency domain signal to the time domain signal to obtain a respiratory motion signal of the object to be detected.

Step 10, determining a plurality of gating phases according to respiratory motion signals of the object to be tested.

And 11, determining a multi-frame original scanning image corresponding to the gating phase from the original scanning image set of the object to be detected aiming at each gating phase.

And step 12, performing fusion processing on the multi-frame original scanning images corresponding to the gating phase to obtain an intermediate scanning image corresponding to the gating phase.

And step 13, determining an intermediate scanning image of the end stage of respiration from the multi-frame intermediate scanning images as a reference scanning image. It should be noted that, the reference scan image may be any one of multiple frames of intermediate scan images.

Step 14, for each intermediate scan image in the multi-frame intermediate scan image, determining a deformation field between the reference scan image and each intermediate scan image.

And step 15, acquiring CT reconstructed images of the object to be detected, and respectively acting deformation fields between the reference scanning image and each intermediate scanning image on the CT reconstructed images to obtain a plurality of new CT reconstructed images.

Step 16, respectively carrying out image reconstruction with attenuation correction on each original scanning image in the original scanning image set of the object to be detected according to a plurality of frames of new CT reconstructed images, and generating scanning images with attenuation correction corresponding to each original scanning image; wherein each of the attenuation corrected scan images may form a first dynamic image sequence.

Step 17, carrying out respiratory motion image registration on the first dynamic image sequence to obtain a second dynamic image sequence; and registering each frame of scanning image except the reference image in the first dynamic image sequence to the reference image by taking any scanning image with attenuation correction in the first dynamic image sequence as the reference image, and generating the second dynamic image sequence according to each registered scanning image and the reference image.

Further, based on the second dynamic image sequence, a subsequent quantitative result analysis can be performed; for example, in the case that the original dynamic image sequence is a cardiac dynamic scan image, the second dynamic image sequence may be combined according to one frame of 5 seconds to obtain a third dynamic image sequence, and an index Myocardial Blood Flow (MBF) of myocardial perfusion may be calculated for the third dynamic image sequence through parameter analysis. Of course, when the third dynamic image sequence is obtained based on the second dynamic image sequence, image merging may be performed at other time intervals, such as a frame of 1 to 10 seconds.

Illustratively, the second dynamic image sequence may be further subjected to gating framing according to an ECG heartbeat signal, and the obtained gating image may be used for cardiac function analysis; such as end-systolic volume ESV, end-diastolic volume EDV, left ventricular ejection fraction LVEF, etc.

Alternatively, the deformation field between the reference scan image and each intermediate scan image obtained in step 14 may be applied to the partial quantitative analysis of the original dynamic image sequence, such as blood pool volume. In addition, besides correcting the image, the deformation field may also be used for gating and framing for subsequent analysis, which is not particularly limited in the embodiment of the present application.

It should be understood that, although the steps in the flowcharts related to the embodiments described above are sequentially shown as indicated by arrows, these steps are not necessarily sequentially performed in the order 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 the flowcharts described in the above embodiments may include a plurality of steps or a plurality of stages, which are not necessarily performed at the same time, but may be performed at different times, and the order of the steps or stages is not necessarily performed sequentially, but may be performed alternately or alternately with at least some of the other steps or stages.

Based on the same inventive concept, the embodiment of the application also provides a medical image processing device for realizing the medical image processing method. The implementation of the solution provided by the device is similar to the implementation described in the above method, so the specific limitations in one or more embodiments of the medical image processing device provided below may be referred to above as limitations of the medical image processing method, and will not be described herein.

In one embodiment, as shown in fig. 6, there is provided a medical image processing apparatus including: a first acquisition module 601, a second acquisition module 602, and a first generation module 603, wherein:

A first obtaining module 601, configured to obtain a motion signal related to an object to be measured from an original scan image set of the object to be measured; the motion signals comprise respiratory motion signals of the object to be detected and nuclide motion signals of radionuclides in the body of the object to be detected; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments through the scanning equipment when the radionuclide in the object to be detected is in a motion state;

the second obtaining module 602 is configured to separate a respiratory motion signal and a nuclide motion signal in the motion signal, and obtain a respiratory motion signal of the object to be detected;

The first generating module 603 is configured to perform respiratory motion correction on an original scan image set of the object to be tested according to the respiratory motion signal of the object to be tested, and generate a corrected scan image.

In one embodiment, the first obtaining module 601 includes a determining unit and a first generating unit; the determining unit is used for determining an interested region from an original scanning image set of the object to be detected; the first generation unit is used for fitting the time activity curve according to the pixel value of the region of interest and generating a time activity curve corresponding to the original scanning image set; the time activity curve carries a motion signal related to the object to be measured.

In one embodiment, the second acquiring module 602 includes a second generating unit, a third generating unit, and a first acquiring unit; the second generation unit is used for converting the time activity curve from a time domain signal to a frequency domain signal and generating a frequency curve corresponding to the time activity curve; the third generating unit is used for carrying out filtering processing on the frequency curve and generating a target frequency curve corresponding to the respiratory motion signal of the object to be detected; the first acquisition unit is used for converting the target frequency curve from a frequency domain signal to a time domain signal to obtain a respiratory motion signal of the object to be detected.

In one embodiment, the generating module includes a reconstructing unit, a second obtaining unit, and a fourth generating unit; the reconstruction unit is used for determining a gating phase according to the respiratory motion signal of the object to be detected, and performing gating reconstruction on an original scanning image set of the object to be detected based on the gating phase to obtain multi-frame intermediate scanning images; the second acquisition unit is used for acquiring CT reconstructed images of the object to be detected; and the fourth generation unit is used for carrying out respiratory motion correction on the original scanning image set of the object to be detected according to the multi-frame intermediate scanning image and the CT reconstruction image of the object to be detected, and generating a corrected scanning image.

In one embodiment, the fourth generating unit is configured to determine a reference scan image from a plurality of frames of intermediate scan images; determining deformation fields between the reference scan image and each intermediate scan image for each of the plurality of frames of intermediate scan images; respectively acting deformation fields between the reference scanning image and each intermediate scanning image on the CT reconstructed images to obtain a plurality of new CT reconstructed images; and carrying out respiratory motion correction on the original scanning image set of the object to be detected according to the multi-frame new CT reconstructed image, and generating a corrected scanning image.

In one embodiment, the original scan image set includes a plurality of frames of original scan images arranged in time sequence, and a reconstruction unit for determining a plurality of gating phases according to respiratory motion signals of the object to be detected; for each gating phase, determining a multi-frame original scanning image corresponding to the gating phase from the original scanning image set of the object to be detected; and carrying out fusion processing on the multi-frame original scanning images corresponding to the gating phase to obtain an intermediate scanning image corresponding to the gating phase.

In one embodiment, the apparatus further includes a third acquisition module, a division module, a reconstruction module, and a second generation module; the third acquisition module is used for acquiring original scanning data of the object to be detected; the dividing module is used for dividing the original scanning data of the object to be detected according to a preset time interval to obtain the original scanning data of a plurality of time periods; the reconstruction module is used for reconstructing images of the original scanning data of each time period according to the original scanning data of each time period to obtain an original scanning image corresponding to the time period; and the second generation module is used for generating an original scanning image set of the object to be detected based on the original scanning image corresponding to each time period.

The respective modules in the above-described medical image processing apparatus may be implemented in whole or in part by software, hardware, and a combination 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, and the computer device may be the medical image scanning device, or may be a terminal device or a server which is in communication connection with the medical image scanning device, and the internal structure diagram may be as shown in fig. 7. The computer device includes a processor, a memory, and a network interface 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, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The database of the computer device is used to store raw scan data. 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 medical image processing method.

It will be appreciated by those skilled in the art that the structure shown in FIG. 7 is merely a block diagram of some of the structures associated with the present inventive arrangements and is not limiting of the computer device to which the present inventive arrangements may be applied, 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 an embodiment a computer device is provided comprising a memory and a processor, the memory having stored therein a computer program, the processor implementing the steps of the medical image processing method in each of the embodiments described above when the computer program is executed.

In an embodiment, a computer-readable storage medium is provided, on which a computer program is stored which, when being executed by a processor, implements the steps of the medical image processing method in the respective embodiments described above.

In an embodiment a computer program product is provided comprising a computer program which, when executed by a processor, implements the steps of the medical image processing method in the respective embodiments described above.

In one embodiment, a PET imaging system is provided, as shown in FIG. 8. The PET imaging system 100 includes a scanning device 10 and a medical image processing apparatus 20, wherein,

A scanning device 10 for acquiring an original scan image set of an object to be measured and transmitting the original scan image set to a medical image processing apparatus; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments when the radionuclide in the object to be detected is in a motion state.

A medical image processing means 20 for acquiring a motion signal related to the object to be measured from the original scan image set; the motion signal comprises a respiratory motion signal of the object to be detected and a nuclide motion signal of the radionuclide in the body of the object to be detected.

The medical image processing device 20 is further configured to separate a respiratory motion signal and a nuclide motion signal in the motion signal, and obtain a respiratory motion signal of the object to be detected; and carrying out respiratory motion correction on the original scanning image set of the object to be detected according to the respiratory motion signal of the object to be detected, and generating a corrected scanning image.

The specific implementation may refer to the embodiment shown in fig. 2, and will not be described herein.

The PET imaging system in the embodiment comprises a scanning device and a medical image processing device, wherein an original scanning image set of an object to be detected in the early stage of radionuclide injection can be acquired through the scanning device, a mixed motion signal which contains a respiratory motion signal of the object to be detected and a nuclide motion signal of the radionuclide in the body of the object to be detected in the early stage of radionuclide injection can be acquired based on the original scanning image set through the medical image processing device, and different motion signals are split on the mixed motion signal in the early stage of radionuclide injection, so that the respiratory motion signal of the object to be detected can be accurately obtained; further, respiratory motion correction is carried out on the original scanning image set based on respiratory motion signals of the object to be detected, and a corrected early PET scanning image is generated; the PET imaging system can effectively extract respiratory motion signals from the original scanning images in the early stage of radionuclide injection, and realize respiratory motion correction on the original scanning images in the early stage of radionuclide injection, so that high-quality scanning images of radionuclides in unsteady state are obtained, and the scanning image quality in the whole scanning process is improved.

In addition, in the PET imaging system of the present embodiment, the medical image processing apparatus 20 is further configured to implement the steps of the medical image processing method of the foregoing embodiments, and specific implementation manner may refer to the foregoing method embodiments, which are not discussed in detail herein.

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, database, or other medium used in embodiments provided herein may include at least one of non-volatile and volatile memory. The nonvolatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical Memory, high density embedded nonvolatile Memory, resistive random access Memory (ReRAM), magneto-resistive random access Memory (Magnetoresistive Random Access Memory, MRAM), ferroelectric Memory (Ferroelectric Random Access Memory, FRAM), phase change Memory (PHASE CHANGE Memory, PCM), graphene Memory, and the like. Volatile memory can include random access memory (Random Access Memory, RAM) or external cache memory, and the like. By way of illustration, and not limitation, RAM can be in various forms such as static random access memory (Static Random Access Memory, SRAM) or dynamic random access memory (Dynamic Random Access Memory, DRAM), etc. The databases referred to in the embodiments provided herein may include at least one of a relational database and a non-relational database. The non-relational database may include, but is not limited to, a blockchain-based distributed database, and the like. The processor referred to in the embodiments provided in the present application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic unit, a data processing logic unit based on quantum computing, or the like, but is not limited thereto.

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 foregoing examples illustrate only a few embodiments of the application and are described in detail herein without thereby 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 the application should be assessed as that of the appended claims.

Claims

1. A medical image processing method, the method comprising:

Acquiring a motion signal related to an object to be detected from an original scanning image set of the object to be detected; the motion signals comprise respiratory motion signals of the object to be detected and nuclide motion signals of radionuclides in the body of the object to be detected; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments through scanning equipment when the radionuclide in the object to be detected is in a motion state;

Separating the respiratory motion signal and the nuclide motion signal in the motion signal to obtain a respiratory motion signal of the object to be detected;

2. The method of claim 1, wherein the acquiring motion signals related to the object from the original set of scanned images of the object comprises:

Fitting a time activity curve according to the pixel value of the region of interest, and generating a time activity curve corresponding to the original scanning image set; the time activity curve carries a motion signal related to the object to be detected.

3. The method of claim 2, wherein separating the respiratory motion signal and the nuclear species motion signal from the motion signal to obtain a respiratory motion signal of the subject comprises:

And converting the target frequency curve from a frequency domain signal to a time domain signal to obtain a respiratory motion signal of the object to be detected.

4. A method according to any one of claims 1 to 3, wherein said respiratory motion correction of the original set of scan images of the subject from the respiratory motion signal of the subject, generating corrected scan images, comprises:

determining a gating phase according to the respiratory motion signal of the object to be detected, and performing gating reconstruction on an original scanning image set of the object to be detected based on the gating phase to obtain multi-frame intermediate scanning images;

acquiring a CT reconstructed image of the object to be detected;

5. The method of claim 4, wherein said performing respiratory motion correction on the original scan image set of the subject from the multi-frame intermediate scan image and the CT reconstructed image of the subject to generate a corrected scan image comprises:

determining, for each of the plurality of intermediate scan images, a deformation field between the reference scan image and each of the intermediate scan images;

6. The method of claim 4, wherein the set of original scan images comprises a plurality of frames of original scan images arranged in a time sequence, wherein the determining a gating phase according to the respiratory motion signal of the object to be measured, and wherein the performing gating reconstruction on the set of original scan images of the object to be measured based on the gating phase, to obtain a plurality of frames of intermediate scan images, comprises:

determining a plurality of gating phases according to the respiratory motion signals of the object to be detected;

For each gating phase, determining a plurality of frames of original scanning images corresponding to the gating phase from the original scanning image set of the object to be detected;

And carrying out fusion processing on the multi-frame original scanning image corresponding to the gating phase to obtain an intermediate scanning image corresponding to the gating phase.

7. The method of claim 6, wherein prior to the acquiring the motion signal of the object from the original scan image set of the object, the method further comprises:

Acquiring original scanning data of the object to be detected;

Dividing the original scanning data of the object to be detected according to a preset time interval to obtain the original scanning data of a plurality of time periods;

Performing image reconstruction on the original scanning data of each time period according to the original scanning data of each time period to obtain an original scanning image corresponding to the time period;

8. A medical image processing apparatus, the apparatus comprising:

The first acquisition module is used for acquiring motion signals related to the object to be detected from the original scanning image set of the object to be detected; the motion signals comprise respiratory motion signals of the object to be detected and nuclide motion signals of radionuclides in the body of the object to be detected; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments through scanning equipment when the radionuclide in the object to be detected is in a motion state;

the second acquisition module is used for separating the respiratory motion signal and the nuclide motion signal in the motion signals to acquire the respiratory motion signal of the object to be detected;

9. 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 7 when the computer program is executed.

10. 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 7.

11. An imaging system comprising a scanning device and medical image processing means, wherein,

The scanning equipment is used for acquiring an original scanning image set of an object to be detected and sending the original scanning image set to the medical image processing device; the original scanning image set is a set of original scanning images obtained by scanning the object to be detected at different moments when the radionuclide in the object to be detected is in a motion state;

The medical image processing device is used for acquiring a motion signal related to the object to be detected from the original scanning image set; the motion signals comprise respiratory motion signals of the object to be detected and nuclide motion signals of radionuclides in the body of the object to be detected;

the medical image processing device is further used for separating the respiratory motion signal and the nuclide motion signal in the motion signals to obtain the respiratory motion signal of the object to be detected; and carrying out respiratory motion correction on the original scanning image set of the object to be detected according to the respiratory motion signal of the object to be detected, and generating a corrected scanning image.