CN114098777A - Method and device for acquiring cardiac phase, storage medium and computer equipment - Google Patents
Method and device for acquiring cardiac phase, storage medium and computer equipment Download PDFInfo
- Publication number
- CN114098777A CN114098777A CN202111149020.1A CN202111149020A CN114098777A CN 114098777 A CN114098777 A CN 114098777A CN 202111149020 A CN202111149020 A CN 202111149020A CN 114098777 A CN114098777 A CN 114098777A
- Authority
- CN
- China
- Prior art keywords
- coronary artery
- phase
- candidate
- region image
- 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.)
- Pending
Links
- 230000000747 cardiac effect Effects 0.000 title claims abstract description 168
- 238000000034 method Methods 0.000 title claims abstract description 45
- 210000004351 coronary vessel Anatomy 0.000 claims abstract description 273
- 230000008859 change Effects 0.000 claims abstract description 73
- 230000033001 locomotion Effects 0.000 claims abstract description 71
- 238000003709 image segmentation Methods 0.000 claims abstract description 15
- 230000011218 segmentation Effects 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 10
- 238000009499 grossing Methods 0.000 claims description 6
- 230000003205 diastolic effect Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 238000002583 angiography Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- 208000031226 Hyperlipidaemia Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 208000031481 Pathologic Constriction Diseases 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 235000005911 diet Nutrition 0.000 description 1
- 230000000378 dietary effect Effects 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 210000002837 heart atrium Anatomy 0.000 description 1
- 201000001421 hyperglycemia Diseases 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 208000037804 stenosis Diseases 0.000 description 1
- 230000036262 stenosis Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Vascular Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Cardiology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The application discloses a method and a device for acquiring a cardiac phase, a storage medium and computer equipment, wherein the method comprises the following steps: acquiring heart images to be segmented corresponding to a plurality of candidate periods; respectively carrying out image segmentation on the heart image to be segmented corresponding to each candidate period to obtain a coronary artery region image corresponding to each candidate period; and determining coronary artery motion change data among the candidate phase phases according to the coronary artery region image, and acquiring a target cardiac phase based on the coronary artery motion change data. According to the method and the device, the optimal target cardiac phase can be automatically obtained according to the motion change condition of the coronary artery, so that a clear and low-artifact coronary artery image can be obtained by reconstructing the image of the target cardiac phase, the workload of a doctor and the dependence on the experience of the doctor are reduced, and the efficiency and the accuracy of image reconstruction are improved.
Description
Technical Field
The present application relates to the field of image reconstruction technologies, and in particular, to a method and an apparatus for acquiring a cardiac phase, a storage medium, and a computer device.
Background
As is well known, cardiovascular and cerebrovascular diseases are diseases which have high morbidity and are easy to cause death in China, and particularly, with the aging of population, the living standard of people is continuously improved and the dietary structure is unbalanced, so that the problems of hyperlipidemia, hypertension, hyperglycemia and the like are easily caused, and coronary heart disease is one of the common diseases. With the continuous development of the CT technology, the rotating speed of a CT frame and the row number of detectors are continuously increased, the CT can bear more scanning tasks, the cardiac coronary artery angiography CTA becomes an important examination means for eliminating the coronary heart disease, and the characteristics of no wound, low cost and no need of hospitalization are easy to be accepted by patients.
However, when image construction for clinical diagnosis is actually performed, a doctor usually selects several target phases to perform image construction according to experience, and the image effect is sometimes not ideal; if the image is reconstructed in the whole scanning phase range, and then the doctor manually distinguishes the cardiac coronary artery angiography image with high quality and low artifact in the optimal phase in all the reconstructed images, the image reconstruction time is long, the CT use efficiency is low, and more time and energy of the doctor are occupied. If the optimal cardiac phase can be automatically obtained, effective assistance is provided for reducing the image reconstruction time, improving the optimal cardiac phase image obtaining efficiency, improving the cardiac image effect and reducing the workload of doctors.
Disclosure of Invention
In view of the above, the present application provides a method and an apparatus for acquiring a cardiac phase, a storage medium, and a computer device.
According to one aspect of the present application, there is provided a method of acquiring cardiac phase, comprising:
acquiring heart images to be segmented corresponding to a plurality of candidate periods;
respectively carrying out image segmentation on the heart image to be segmented corresponding to each candidate period to obtain a coronary artery region image corresponding to each candidate period;
and determining coronary artery motion change data among the candidate phase phases according to the coronary artery region image, and acquiring a target cardiac phase based on the coronary artery motion change data.
Optionally, the acquiring the to-be-segmented cardiac images corresponding to the multiple phases specifically includes:
acquiring cardiac scanning data, wherein a scanning phase range of the cardiac scanning data is determined based on a preset reconstruction phase range;
and according to the interval of the preset period, reconstructing the cardiac scanning data corresponding to the phase range of the preset reconstruction period, and determining the cardiac images to be segmented corresponding to a plurality of candidate periods.
Optionally, the determining, according to the coronary artery region image, coronary artery motion change data between candidate phases specifically includes:
acquiring a coronary artery region image i and a coronary artery region image i +1 corresponding to any two continuous candidate periods, wherein the coronary artery region image i represents a coronary artery region image corresponding to a candidate period phase i, the coronary artery region image i +1 represents a coronary artery region image corresponding to the candidate period phase i +1, and i is smaller than the number of the candidate period phases;
taking the intersection of the coronary artery region image i and the coronary artery region image i +1 as a comparison template i, determining coronary artery contour data i based on the coronary artery region image i and the comparison template i, and determining coronary artery contour data i +1 based on the coronary artery region image i +1 and the comparison template i;
and determining the coronary artery motion change data corresponding to the candidate phase i +1 according to the coronary artery contour data i +1 and the coronary artery contour data i.
Optionally, the acquiring a target cardiac phase based on the coronary artery motion change data specifically includes:
drawing a coronary artery change curve corresponding to the preset reconstruction period phase range according to the coronary artery motion change data;
and acquiring a phase corresponding to the minimum value in the coronary artery change curve as the target cardiac phase.
Optionally, the acquiring a phase corresponding to a minimum value in the coronary artery variation curve as the target cardiac phase specifically includes:
taking a preset cardiac systole as a first starting search point, searching a first minimum extreme point in a first search phase range corresponding to the preset cardiac systole, and acquiring a phase corresponding to the first minimum extreme point as a first target cardiac phase, wherein the target cardiac phase comprises the first target cardiac phase; and/or the presence of a gas in the gas,
and taking a preset diastole as a second initial search point, searching a second minimum extreme point in a second search phase range corresponding to the preset diastole, and acquiring a phase corresponding to the second minimum extreme point as a second target cardiac phase, wherein the target cardiac phase comprises the second target cardiac phase.
Optionally, the image segmentation is performed on the to-be-segmented cardiac image corresponding to each candidate period to obtain a coronary artery region image corresponding to each candidate period, and specifically includes:
respectively carrying out image segmentation on the multilayer heart images to be segmented corresponding to each candidate period to obtain multilayer coronary artery region images corresponding to each candidate period;
and performing smoothing treatment on the segmentation edges of the multilayer coronary artery region image corresponding to each candidate period to determine the coronary artery region image corresponding to each candidate period.
Optionally, after obtaining the image of the coronary artery region corresponding to each candidate period, the method further includes:
respectively determining the heart centroid corresponding to each coronary artery region image according to the coronary artery region image corresponding to each candidate period;
acquiring at least one target coronary artery region image in each coronary artery region image according to the position relation between the center of mass of the heart and at least one target coronary artery;
correspondingly, the determining the coronary artery motion change data between the candidate phases according to the coronary artery region image specifically includes:
and determining coronary artery motion change data among the candidate phases according to the target coronary artery region image.
According to another aspect of the present application, there is provided an acquisition apparatus of a cardiac phase, comprising:
the image acquisition module is used for acquiring cardiac images to be segmented corresponding to a plurality of candidate periods;
the coronary artery segmentation module is used for respectively carrying out image segmentation on the heart image to be segmented corresponding to each candidate period to obtain a coronary artery region image corresponding to each candidate period;
and the phase acquisition module is used for determining coronary artery motion change data among the candidate phases according to the coronary artery region image and acquiring a target cardiac phase based on the coronary artery motion change data.
Optionally, the image obtaining module is specifically configured to:
acquiring cardiac scanning data, wherein a scanning phase range of the cardiac scanning data is determined based on a preset reconstruction phase range;
and according to the interval of the preset period, reconstructing the cardiac scanning data corresponding to the phase range of the preset reconstruction period, and determining the cardiac images to be segmented corresponding to a plurality of candidate periods.
Optionally, the phase acquisition module is specifically configured to:
acquiring a coronary artery region image i and a coronary artery region image i +1 corresponding to any two continuous candidate periods, wherein the coronary artery region image i represents a coronary artery region image corresponding to a candidate period phase i, the coronary artery region image i +1 represents a coronary artery region image corresponding to the candidate period phase i +1, and i is smaller than the number of the candidate period phases;
taking the intersection of the coronary artery region image i and the coronary artery region image i +1 as a comparison template i, determining coronary artery contour data i based on the coronary artery region image i and the comparison template i, and determining coronary artery contour data i +1 based on the coronary artery region image i +1 and the comparison template i;
and determining the coronary artery motion change data corresponding to the candidate phase i +1 according to the coronary artery contour data i +1 and the coronary artery contour data i.
Optionally, the phase acquisition module is specifically configured to:
drawing a coronary artery change curve corresponding to the preset reconstruction period phase range according to the coronary artery motion change data;
and acquiring a phase corresponding to the minimum value in the coronary artery change curve as the target cardiac phase.
Optionally, the phase obtaining module is further configured to:
taking a preset cardiac systole as a first starting search point, searching a first minimum extreme point in a first search phase range corresponding to the preset cardiac systole, and acquiring a phase corresponding to the first minimum extreme point as a first target cardiac phase, wherein the target cardiac phase comprises the first target cardiac phase; and/or the presence of a gas in the gas,
and taking a preset diastole as a second initial search point, searching a second minimum extreme point in a second search phase range corresponding to the preset diastole, and acquiring a phase corresponding to the second minimum extreme point as a second target cardiac phase, wherein the target cardiac phase comprises the second target cardiac phase.
Optionally, the coronary artery segmentation module is specifically configured to:
respectively carrying out image segmentation on the multilayer heart images to be segmented corresponding to each candidate period to obtain multilayer coronary artery region images corresponding to each candidate period;
and performing smoothing treatment on the segmentation edges of the multilayer coronary artery region image corresponding to each candidate period to determine the coronary artery region image corresponding to each candidate period.
Optionally, the coronary segmentation module is further configured to:
after obtaining the coronary artery region image corresponding to each candidate period, respectively determining the heart centroid corresponding to each coronary artery region image according to the coronary artery region image corresponding to each candidate period; acquiring at least one target coronary artery region image in each coronary artery region image according to the position relation between the center of mass of the heart and at least one target coronary artery;
accordingly, the phase acquisition module is further configured to: and determining coronary artery motion change data among the candidate phases according to the target coronary artery region image.
According to yet another aspect of the present application, a storage medium is provided, on which a computer program is stored which, when being executed by a processor, carries out the above-mentioned method of acquiring a cardiac phase.
According to yet another aspect of the present application, there is provided a computer device comprising a storage medium, a processor and a computer program stored on the storage medium and executable on the processor, the processor implementing the above method for acquiring cardiac phases when executing the program.
By means of the technical scheme, the method and the device for acquiring the cardiac phase, the storage medium and the computer device provided by the application perform image reconstruction on the scanning data, acquire the to-be-segmented cardiac images corresponding to a plurality of candidate phases, segment the coronary artery region image in the image, and further count the motion data of the coronary artery between the candidate phases according to the coronary artery region image, so that the optimal target cardiac phase is selected according to the change data. Compared with the prior art, the method has the advantages that the doctor can automatically acquire the optimal target cardiac phase according to the motion change condition of the coronary artery in a mode of manually selecting the full-phase reconstruction image according to experience or the motion change condition of the coronary artery, so that the clear and low-artifact coronary artery image can be acquired by reconstructing the image of the target cardiac phase, the workload of the doctor and the dependence on the doctor experience are reduced, and the efficiency and the accuracy of image reconstruction are improved.
The foregoing description is only an overview of the technical solutions of the present application, and the present application can be implemented according to the content of the description in order to make the technical means of the present application more clearly understood, and the following detailed description of the present application is given in order to make the above and other objects, features, and advantages of the present application more clearly understandable.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the application and together with the description serve to explain the application and not to limit the application. In the drawings:
fig. 1 is a schematic flow chart illustrating a method for acquiring a cardiac phase according to an embodiment of the present disclosure;
fig. 2 is a schematic flow chart illustrating a coronary artery region segmentation provided by an embodiment of the present application;
fig. 3 is a schematic flow chart illustrating another method for acquiring cardiac phases according to an embodiment of the present disclosure;
FIG. 4 is a schematic diagram illustrating a coronary artery variation curve provided by an embodiment of the present application;
fig. 5 is a schematic structural diagram illustrating an apparatus for acquiring a cardiac phase according to an embodiment of the present disclosure.
Detailed Description
The present application will be described in detail below with reference to the accompanying drawings in conjunction with embodiments. It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict.
In this embodiment, a method for acquiring a cardiac phase is provided, as shown in fig. 1, the method including:
in the embodiment of the present application, the cardiac image to be segmented may be obtained by reconstructing CT scan data. In a specific application scenario, a prospective cardiac gating mode may be used to perform any phase scanning, such as a full-phase mode scanning covering the whole heartbeat cycle in one scanning time, or a retrospective cardiac gating spiral acquisition mode may be used to perform the full-phase mode scanning.
The candidate phase may be a plurality of all scan phases, for example, when the cardiac image to be segmented is obtained according to the scan data of the full phase, the full phase is 0% to 100%, and 20 of the full phase may be uniformly selected for image reconstruction to obtain the cardiac image to be segmented. The relative phase of the cardiac scan refers to the corresponding time range of the image reconstructed by the heart on the ECG signal, which accounts for the percentage of the range of two adjacent R peaks of the ECG signal, 0% corresponds to the position of the first R peak, 100% corresponds to the position of the second R peak, and the full phase covers 0% -100%.
Optionally, step 101 may specifically include: acquiring cardiac scanning data, wherein a scanning phase range of the cardiac scanning data is determined based on a preset reconstruction phase range; and according to the interval of the preset period, reconstructing the cardiac scanning data corresponding to the phase range of the preset reconstruction period, and determining the cardiac images to be segmented corresponding to a plurality of candidate periods.
In this embodiment, the cardiac scanning data may be data obtained by performing CT scanning in any mode, and since scanning data in a certain range before and after the phase of the image reconstruction period is required to be used in image reconstruction, in an actual application scenario, the phase range of the cardiac scanning period may be determined according to a preset reconstruction period phase range, for example, the preset reconstruction period phase range is 5% to 95%, and then the scanning period phase range of the heart may be 0% to 100%. And then, according to a preset phase interval, determining a plurality of candidate phase within a preset reconstruction phase range, and reconstructing the heart image to be segmented corresponding to each candidate phase, so as to reduce the image reconstruction workload and improve the efficiency. For example, the preset period interval may be selected to be 5%, and the candidate period phase may include 5%, 10%, 15% … …
In addition, in order to improve the image reconstruction efficiency, the range of the CT scanning is possibly large, the visual field range of the image reconstruction can be adjusted, and only images in a certain range around the heart need to be reconstructed, so that the image reconstruction efficiency is improved.
102, respectively carrying out image segmentation on the heart image to be segmented corresponding to each candidate period to obtain a coronary artery region image corresponding to each candidate period;
the main focus of cardiac image reconstruction is whether coronary arteries of a heart have lesions, plaques, stenosis and blockage, in the contraction and relaxation processes of the heart, the motion states of atria and ventricles are different, and the actual motion state of the coronary arteries cannot necessarily represent only by the size change condition of the whole heart contour, so that the obtained optimal period is inaccurate, and the reconstructed coronary arteries of the image have motion artifacts, which affects the evaluation and diagnosis of a doctor. Therefore, in the embodiment of the application, after the to-be-segmented cardiac image corresponding to each candidate phase is determined, each image can be segmented respectively to segment the coronary artery region image, so that the optimal cardiac phase is analyzed by using the coronary artery region image in the following process, and the accuracy of determining the optimal cardiac phase and the accuracy of image reconstruction are improved. In an actual application scenario, the segmentation of the coronary artery region may be realized by a coronary artery segmentation model, or may be realized by a coronary artery feature recognition mode, which is not limited herein.
And finally, determining the change data of the coronary artery between the candidate phase phases according to the segmented coronary artery region image, searching the phase with the most gentle motion state of the coronary artery according to the change data, and taking the phase with the most gentle motion state as the final target cardiac phase so as to reconstruct the image of the cardiac image of the target cardiac phase and obtain a high-quality and low-artifact cardiac reconstruction image.
In this embodiment of the present application, optionally, after step 103, the method may further include:
and 104, carrying out image reconstruction according to the cardiac scanning data and the target cardiac phase, and determining a cardiac reconstruction image corresponding to the target cardiac phase.
In addition, in this embodiment of the application, optionally, the step 103 and the step 104 may further include: the target cardiac phase is output and step 104 is performed in response to the image reconstruction instructions. The image reconstruction instruction may be a confirmation instruction for the target cardiac phase, or may be a data modification instruction for the target cardiac phase.
By applying the technical scheme of the embodiment, image reconstruction is performed on the scanning data, to-be-segmented cardiac images corresponding to a plurality of candidate phases are obtained, coronary artery region images in the images are segmented, and further, the motion data of coronary arteries among the candidate phases is counted according to the coronary artery region images, so that the optimal target cardiac phase is selected according to the change data. Compared with the mode that a doctor manually selects the full-phase reconstruction image in the prior art, the method and the device can automatically acquire the optimal target cardiac phase according to the motion change condition of the coronary artery, so that the clear and low-artifact coronary artery image can be acquired by reconstructing the image of the target cardiac phase, the workload of the doctor and the dependence on the experience of the doctor are reduced, and the efficiency and the accuracy of image reconstruction are improved.
In this embodiment of the application, optionally, the "determining coronary artery motion variation data between candidate phases according to the coronary artery region image" in step 103 may specifically include:
103-1, acquiring a coronary artery region image i and a coronary artery region image i +1 corresponding to any two continuous candidate periods, wherein the coronary artery region image i represents a coronary artery region image corresponding to the candidate period i, and the coronary artery region image i +1 represents a coronary artery region image corresponding to the candidate period i + 1;
step 103-2, taking an intersection of the coronary artery region image i and the coronary artery region image i +1 as a comparison template i, determining coronary artery contour data i based on the coronary artery region image i and the comparison template i, and determining coronary artery contour data i +1 based on the coronary artery region image i +1 and the comparison template i; and determining the coronary artery motion change data corresponding to the candidate phase i +1 according to the coronary artery contour data i +1 and the coronary artery contour data i.
In this embodiment, the coronary artery variation between the candidate phase i and the candidate phase i +1 can be determined by: firstly, obtaining a coronary artery region image i and a coronary artery region image i +1 corresponding to a candidate phase i, then, taking intersection of the coronary artery region image i and the coronary artery region image i +1, taking the intersection part as a comparison template i, performing coronary artery identification on an overlapped region of the coronary artery region image i and the comparison template i to obtain corresponding coronary artery contour data i, performing coronary artery identification on an overlapped region of the coronary artery region image i +1 and the comparison template i in the same way to obtain corresponding coronary artery contour data i +1, finally, subtracting the coronary artery contour data i +1 from the coronary artery contour data i, and taking the difference value as change data of the coronary artery from the candidate phase i to the candidate phase i + 1. Specifically, the coronary motion variation data may be used as the coronary motion variation corresponding to the phase [ i + (i +1) ]/2, so as to facilitate the subsequent screening of the optimal cardiac phase, for example, the phase with the smallest motion variation may be selected as the optimal cardiac phase.
In this embodiment of the present application, optionally, the "acquiring the target cardiac phase based on the coronary artery motion variation data" in step 103 may specifically include:
103-3, drawing a coronary artery change curve corresponding to the preset reconstruction period phase range according to the coronary artery motion change data;
and 103-4, acquiring a phase corresponding to the minimum value in the coronary artery change curve as the target cardiac phase.
In the above embodiment, in order to analyze the motion change conditions of all phases within the preset reconstruction phase range, the coronary artery motion change data corresponding to a plurality of candidate phases may be used to fit the corresponding coronary artery change curve within the preset reconstruction phase range, and specifically, the coronary artery motion change data may be subjected to interpolation processing to obtain the coronary artery change values corresponding to more phase points, improve the smoothness of the curve, and eliminate the influence of interference factors on the curve. Further, a minimum value corresponding to the coronary artery change curve is determined, a phase corresponding to the minimum value is used as an optimal target cardiac phase, and image reconstruction can be subsequently performed according to the target cardiac phase to obtain a high-quality and low-artifact cardiac coronary artery image.
In this embodiment of the application, in order to further improve the reliability of the optimal phase, optionally, step 103-4 may specifically include:
s1, taking a preset cardiac systole as a first initial search point, searching a first minimum extreme point in a first search phase range corresponding to the preset cardiac systole, and acquiring a phase corresponding to the first minimum extreme point as a first target cardiac phase, wherein the target cardiac phase comprises the first target cardiac phase; and/or the presence of a gas in the gas,
s2, taking a preset diastole as a second start search point, searching a second minimum extreme point in a second search phase range corresponding to the preset diastole, and obtaining a phase corresponding to the second minimum extreme point as a second target cardiac phase, where the target cardiac phase includes the second target cardiac phase.
In this embodiment, based on the heartbeat rule, the optimal cardiac phase is generally the systolic phase or the diastolic phase, and the approximate range of the systolic phase and the diastolic phase is determined, and generally the systolic phase is about 40% and the diastolic phase is about 75%, so that the target cardiac phase can be searched using the systolic phase and the diastolic phase as the starting search points. In the schematic diagram of the coronary artery variation curve shown in fig. 4, the horizontal axis represents the phase, and the vertical axis represents the amount of coronary artery motion variation, which may be in Hounsfield Unit (HU), such that the preset systolic period fsystole location (e.g. 40%) and the preset diastolic period fdistallelocation (e.g. 75%) estimated according to the empirical values are taken as the first start search point and the preset diastolic period fdistallelocation is taken as the second start search point. The first minimum extreme point is searched near the first search starting point (for example, 30% -50%), the second minimum extreme point is searched near the second search starting point (for example, 60% -90%), and the first minimum extreme point and the second minimum extreme point obtained through searching are used as the most gentle phase of the coronary artery movement, namely the optimal systolic phase OptimalSyle and the optimal diastolic phase OptimalDiastole, namely the first target cardiac phase and the second target cardiac phase. By the method, the phase deviation and the theoretical range of the found target cardiac phase caused by factors such as interference can be avoided, and the reliability of the optimal phase is improved.
In this embodiment of the present application, optionally, step 102 may specifically include:
102-1, respectively carrying out image segmentation on a plurality of layers of heart images to be segmented corresponding to each candidate period to obtain a plurality of layers of coronary artery region images corresponding to each candidate period;
and 102-2, performing smoothing treatment on the segmentation edges of the multilayer coronary artery region image corresponding to each candidate period, and determining the coronary artery region image corresponding to each candidate period.
In this embodiment, the coronary artery region image may be obtained through a pre-trained model, as shown in fig. 2, a model is trained by using a sample set to obtain a shape model of a target region, a classifier is constructed by using the model, and the classifier is used to classify a coronary artery region in an input image, specifically, after the image is input to the classifier, the classifier may classify whether each pixel point in the image belongs to a coronary artery, so as to generate a segmentation template (i.e., a coronary artery region). In an actual application scene, layered scanning is generally performed when the heart is scanned, a reconstructed image corresponding to each candidate period comprises a plurality of layers of heart images to be segmented, then coronary artery segmentation is performed on each layer of heart images to be segmented respectively, smoothing processing is performed on edge feature points of each layer of coronary artery obtained by segmentation, the edge feature points of the coronary artery regions of the adjacent layers are coherent and smooth, and therefore a three-dimensional coronary artery region image can be constructed according to the plurality of layers of coronary artery region images corresponding to each candidate period, and then coronary artery motion change data calculation and target heart period phase acquisition are performed.
Further, another method for acquiring a cardiac phase is provided in an embodiment of the present application, as shown in fig. 3, the method includes:
In this embodiment, the coronary arteries mainly include a right coronary artery, a left circumflex branch and a left anterior descending branch, and since the motion change conditions of different coronary arteries may be different, in an actual application scenario, the target cardiac phase may be found according to the overall motion change condition of the coronary artery, or according to the motion change condition of one or more coronary arteries. After obtaining the coronary artery region images corresponding to a plurality of candidate periods, the heart centroid position can be determined firstly, the target coronary artery can be determined according to actual needs, and then different target coronary artery region images can be obtained by combining the position relation of the heart centroid position and various target coronary arteries, wherein the right coronary artery is generally positioned at the right side of the heart centroid, the left anterior descending branch is positioned at the upper part of the left side of the heart centroid, and the left circumflex branch is positioned at the lower part of the left side of the heart centroid. Further, in the phase of searching the target cardiac phase, if the target coronary comprises one coronary, the target cardiac phase is determined directly according to the motion change condition of the target coronary. If the target coronary artery includes a plurality of kinds, a plurality of target cardiac phases may be determined according to the motion change condition of each target coronary artery, and one of the target cardiac phases may be selected as a final optimal phase. For example, the target coronary artery comprises a left circumflex and a left anterior descending branch, the target cardiac phase is determined to comprise a phase A corresponding to the left circumflex and a phase B corresponding to the left anterior descending branch by searching a minimum value point near the preset cardiac systole, and the phase closest to the preset cardiac systole is selected as a final optimal phase. The determination method of the motion change condition of the target coronary artery is similar to the determination method of the coronary artery motion change data, and is not described herein again. By the method, the optimal phase search is carried out, the problem that the motion change situation description is inaccurate due to the fact that the motion changes of various coronary arteries are mutually offset when the motion change situations of the various coronary arteries are calculated as a whole can be avoided, and the acquisition accuracy of the target cardiac phase is further improved.
Further, as a specific implementation of the method in fig. 1, an embodiment of the present application provides an apparatus for acquiring cardiac phase, as shown in fig. 5, the apparatus includes:
the image acquisition module is used for acquiring cardiac images to be segmented corresponding to a plurality of candidate periods;
the coronary artery segmentation module is used for respectively carrying out image segmentation on the heart image to be segmented corresponding to each candidate period to obtain a coronary artery region image corresponding to each candidate period;
and the phase acquisition module is used for determining coronary artery motion change data among the candidate phases according to the coronary artery region image and acquiring a target cardiac phase based on the coronary artery motion change data.
Optionally, the image obtaining module is specifically configured to:
acquiring cardiac scanning data, wherein a scanning phase range of the cardiac scanning data is determined based on a preset reconstruction phase range;
and according to the interval of the preset period, reconstructing the cardiac scanning data corresponding to the phase range of the preset reconstruction period, and determining the cardiac images to be segmented corresponding to a plurality of candidate periods.
Optionally, the phase acquisition module is specifically configured to:
acquiring a coronary artery region image i and a coronary artery region image i +1 corresponding to any two continuous candidate periods, wherein the coronary artery region image i represents a coronary artery region image corresponding to a candidate period phase i, the coronary artery region image i +1 represents a coronary artery region image corresponding to the candidate period phase i +1, and i is smaller than the number of the candidate period phases;
and taking the intersection of the coronary artery region image i and the coronary artery region image i +1 as a comparison image i +1, and determining the coronary artery motion change data corresponding to the candidate phase i + 1.
Optionally, the phase acquisition module is specifically configured to:
drawing a coronary artery change curve corresponding to the preset reconstruction period phase range according to the coronary artery motion change data;
and acquiring a phase corresponding to the minimum value in the coronary artery change curve as the target cardiac phase.
Optionally, the phase obtaining module is further configured to:
taking a preset cardiac systole as a first starting search point, searching a first minimum extreme point in a first search phase range corresponding to the preset cardiac systole, and acquiring a phase corresponding to the first minimum extreme point as a first target cardiac phase, wherein the target cardiac phase comprises the first target cardiac phase; and/or the presence of a gas in the gas,
and taking a preset diastole as a second initial search point, searching a second minimum extreme point in a second search phase range corresponding to the preset diastole, and acquiring a phase corresponding to the second minimum extreme point as a second target cardiac phase, wherein the target cardiac phase comprises the second target cardiac phase.
Optionally, the coronary artery segmentation module is specifically configured to:
respectively carrying out image segmentation on the multilayer heart images to be segmented corresponding to each candidate period to obtain multilayer coronary artery region images corresponding to each candidate period;
and performing smoothing treatment on the segmentation edges of the multilayer coronary artery region image corresponding to each candidate period to determine the coronary artery region image corresponding to each candidate period.
Optionally, the coronary segmentation module is further configured to:
after obtaining the coronary artery region image corresponding to each candidate period, respectively determining the heart centroid corresponding to each coronary artery region image according to the coronary artery region image corresponding to each candidate period; acquiring at least one target coronary artery region image in each coronary artery region image according to the position relation between the center of mass of the heart and at least one target coronary artery;
accordingly, the phase acquisition module is further configured to: and determining coronary artery motion change data among the candidate phases according to the target coronary artery region image.
It should be noted that, other corresponding descriptions of the functional units involved in the apparatus for acquiring cardiac phase provided in the embodiment of the present application may refer to the corresponding descriptions in the methods of fig. 1 to fig. 3, and are not repeated herein.
Based on the above-mentioned method shown in fig. 1 to 3, correspondingly, the present application further provides a storage medium, on which a computer program is stored, and the computer program, when executed by a processor, implements the above-mentioned method for acquiring cardiac phases shown in fig. 1 to 3.
Based on such understanding, the technical solution of the present application may be embodied in the form of a software product, which may be stored in a non-volatile storage medium (which may be a CD-ROM, a usb disk, a removable hard disk, etc.), and includes several instructions for enabling a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the method according to the implementation scenarios of the present application.
Based on the above methods shown in fig. 1 to fig. 3 and the virtual device embodiment shown in fig. 5, in order to achieve the above object, an embodiment of the present application further provides a computer device, which may specifically be a personal computer, a server, a network device, and the like, where the computer device includes a storage medium and a processor; a storage medium for storing a computer program; a processor for executing a computer program for implementing the above-described method for acquiring a cardiac phase as shown in fig. 1 to 3.
Optionally, the computer device may also include a user interface, a network interface, a camera, Radio Frequency (RF) circuitry, sensors, audio circuitry, a WI-FI module, and so forth. The user interface may include a Display screen (Display), an input unit such as a keypad (Keyboard), etc., and the optional user interface may also include a USB interface, a card reader interface, etc. The network interface may optionally include a standard wired interface, a wireless interface (e.g., a bluetooth interface, WI-FI interface), etc.
It will be appreciated by those skilled in the art that the present embodiment provides a computer device architecture that is not limiting of the computer device, and that may include more or fewer components, or some components in combination, or a different arrangement of components.
The storage medium may further include an operating system and a network communication module. An operating system is a program that manages and maintains the hardware and software resources of a computer device, supporting the operation of information handling programs, as well as other software and/or programs. The network communication module is used for realizing communication among components in the storage medium and other hardware and software in the entity device.
Through the above description of the embodiments, those skilled in the art can clearly understand that the present application may be implemented by software plus a necessary general hardware platform, or may implement image reconstruction on scan data by hardware, obtain cardiac images to be segmented corresponding to a plurality of candidate phases, segment coronary artery region images in the images, and further calculate motion data of coronary arteries between the candidate phases according to the coronary artery region images, thereby selecting an optimal target cardiac phase according to the change data. Compared with the prior art, the method has the advantages that the doctor can automatically acquire the optimal target cardiac phase according to the motion change condition of the coronary artery in a mode of manually selecting the full-phase reconstruction image according to experience or the motion change condition of the coronary artery, so that the clear and low-artifact coronary artery image can be acquired by reconstructing the image of the target cardiac phase, the workload of the doctor and the dependence on the doctor experience are reduced, and the efficiency and the accuracy of image reconstruction are improved.
Those skilled in the art will appreciate that the figures are merely schematic representations of one preferred implementation scenario and that the blocks or flow diagrams in the figures are not necessarily required to practice the present application. Those skilled in the art will appreciate that the modules in the devices in the implementation scenario may be distributed in the devices in the implementation scenario according to the description of the implementation scenario, or may be located in one or more devices different from the present implementation scenario with corresponding changes. The modules of the implementation scenario may be combined into one module, or may be further split into a plurality of sub-modules.
The above application serial numbers are for description purposes only and do not represent the superiority or inferiority of the implementation scenarios. The above disclosure is only a few specific implementation scenarios of the present application, but the present application is not limited thereto, and any variations that can be made by those skilled in the art are intended to fall within the scope of the present application.
Claims (10)
1. A method for acquiring a cardiac phase, comprising:
acquiring heart images to be segmented corresponding to a plurality of candidate periods;
respectively carrying out image segmentation on the heart image to be segmented corresponding to each candidate period to obtain a coronary artery region image corresponding to each candidate period;
and determining coronary artery motion change data among the candidate phase phases according to the coronary artery region image, and acquiring a target cardiac phase based on the coronary artery motion change data.
2. The method according to claim 1, wherein the acquiring of the cardiac images to be segmented corresponding to a plurality of phases specifically comprises:
acquiring cardiac scanning data, wherein a scanning phase range of the cardiac scanning data is determined based on a preset reconstruction phase range;
and according to the interval of the preset period, reconstructing the cardiac scanning data corresponding to the phase range of the preset reconstruction period, and determining the cardiac images to be segmented corresponding to a plurality of candidate periods.
3. The method according to claim 1, wherein determining coronary motion variation data between candidate phases based on the coronary artery region image comprises:
acquiring a coronary artery region image i and a coronary artery region image i +1 corresponding to any two continuous candidate periods, wherein the coronary artery region image i represents a coronary artery region image corresponding to a candidate period phase i, and the coronary artery region image i +1 represents a coronary artery region image corresponding to the candidate period phase i + 1;
taking the intersection of the coronary artery region image i and the coronary artery region image i +1 as a comparison template i, determining coronary artery contour data i based on the coronary artery region image i and the comparison template i, and determining coronary artery contour data i +1 based on the coronary artery region image i +1 and the comparison template i;
and determining the coronary artery motion change data corresponding to the candidate phase i +1 according to the coronary artery contour data i +1 and the coronary artery contour data i.
4. The method according to claim 3, wherein the obtaining a target cardiac phase based on the coronary motion variation data comprises:
drawing a coronary artery change curve corresponding to the preset reconstruction period phase range according to the coronary artery motion change data;
and acquiring a phase corresponding to the minimum value in the coronary artery change curve as the target cardiac phase.
5. The method according to claim 4, wherein the obtaining a phase corresponding to a minimum value in the coronary variation curve as the target cardiac phase comprises:
taking a preset cardiac systole as a first starting search point, searching a first minimum extreme point in a first search phase range corresponding to the preset cardiac systole, and acquiring a phase corresponding to the first minimum extreme point as a first target cardiac phase, wherein the target cardiac phase comprises the first target cardiac phase; and/or the presence of a gas in the gas,
and taking a preset diastole as a second initial search point, searching a second minimum extreme point in a second search phase range corresponding to the preset diastole, and acquiring a phase corresponding to the second minimum extreme point as a second target cardiac phase, wherein the target cardiac phase comprises the second target cardiac phase.
6. The method according to claim 1, wherein the image segmentation is performed on the cardiac image to be segmented corresponding to each candidate period to obtain a coronary artery region image corresponding to each candidate period, specifically including:
respectively carrying out image segmentation on the multilayer heart images to be segmented corresponding to each candidate period to obtain multilayer coronary artery region images corresponding to each candidate period;
and performing smoothing treatment on the segmentation edges of the multilayer coronary artery region image corresponding to each candidate period to determine the coronary artery region image corresponding to each candidate period.
7. The method of claim 1, wherein after obtaining the image of the coronary artery region corresponding to each candidate period, the method further comprises:
respectively determining the heart centroid corresponding to each coronary artery region image according to the coronary artery region image corresponding to each candidate period;
acquiring at least one target coronary artery region image in each coronary artery region image according to the position relation between the center of mass of the heart and at least one target coronary artery;
correspondingly, the determining the coronary artery motion change data between the candidate phases according to the coronary artery region image specifically includes:
and determining coronary artery motion change data among the candidate phases according to the target coronary artery region image.
8. An apparatus for acquiring a cardiac phase, comprising:
the image acquisition module is used for acquiring cardiac images to be segmented corresponding to a plurality of candidate periods;
the coronary artery segmentation module is used for respectively carrying out image segmentation on the heart image to be segmented corresponding to each candidate period to obtain a coronary artery region image corresponding to each candidate period;
and the phase acquisition module is used for determining coronary artery motion change data among the candidate phases according to the coronary artery region image and acquiring a target cardiac phase based on the coronary artery motion change data.
9. A storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the method of any of claims 1 to 7.
10. A computer device comprising a storage medium, a processor and a computer program stored on the storage medium and executable on the processor, characterized in that the processor implements the method of any one of claims 1 to 7 when executing the computer program.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111149020.1A CN114098777A (en) | 2021-09-29 | 2021-09-29 | Method and device for acquiring cardiac phase, storage medium and computer equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111149020.1A CN114098777A (en) | 2021-09-29 | 2021-09-29 | Method and device for acquiring cardiac phase, storage medium and computer equipment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114098777A true CN114098777A (en) | 2022-03-01 |
Family
ID=80441668
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111149020.1A Pending CN114098777A (en) | 2021-09-29 | 2021-09-29 | Method and device for acquiring cardiac phase, storage medium and computer equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114098777A (en) |
-
2021
- 2021-09-29 CN CN202111149020.1A patent/CN114098777A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10980502B2 (en) | Method, storage medium, and system for analyzing image sequences of periodic physiological activities | |
| US20070053482A1 (en) | Reconstruction of an image of a moving object from volumetric data | |
| JP5060719B2 (en) | Biological organ image processing apparatus, method, and computer program | |
| US20220296205A1 (en) | Ultrasound image-based diagnosis system for coronary artery lesion using machine learning and diagnosis method of same | |
| JP2004509686A (en) | Selection of medical images based on image data | |
| CA3104074A1 (en) | Methods and systems to quantify and remove asynchronous noise in biophysical signals | |
| CN110969633B (en) | Automatic optimal phase identification method for cardiac CT imaging | |
| CN112184690B (en) | Coronary vessel trend prediction method, prediction model training method and device | |
| CN113205880B (en) | LogitBoost-based heart disease prognosis prediction method and device | |
| CN110728730A (en) | Image reconstruction method and device, CT (computed tomography) equipment and CT system | |
| KR20200113336A (en) | Learning method for generating multiphase collateral image and multiphase collateral imaging method using machine learning | |
| CN111369636B (en) | Image reconstruction method, system, device and storage medium | |
| JP2011043489A (en) | Image processing apparatus, program and method for evaluating local wall motion state of cardiac muscle | |
| KR102257295B1 (en) | Diagnostic system for diagnosing vulnerable atheromatous plaque through ultrasound image-based machine learning and the diagnostic method thereof | |
| CN114098777A (en) | Method and device for acquiring cardiac phase, storage medium and computer equipment | |
| CN115564817A (en) | Three-dimensional blood vessel image determination method, device, equipment and storage medium | |
| CN113538419B (en) | Image processing method and system | |
| CN106910182A (en) | In cardiac function MRI in end diastole image blood pool dividing method | |
| CN112862827B (en) | Method, device, terminal and storage medium for determining opening parameters of left atrial appendage | |
| CN113116385B (en) | Arterial intracavity ultrasonic imaging method, device, equipment and system | |
| CN103006218A (en) | Data processing method and system for cardiac magnetic resonance real-time film imaging | |
| Wang et al. | Deep learning for automated echocardiogram analysis | |
| KR20180008130A (en) | A method for correcting a motion in a coronary image using a convolutional neural network | |
| CN111755104A (en) | Heart state monitoring method and system, electronic equipment and storage medium | |
| CN116831623B (en) | Ultrasound image correction method and device, electronic equipment and storage medium |
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 |