Cardiovascular imaging phase determination method and apparatus, electronic device, and storage medium

Abstract

The application relates to a cardiac blood vessel imaging phase determination method, which comprises the following steps: acquiring multiple phase images of a cardiac blood vessel; extracting target coronary arteries based on the multiple phase images respectively to obtain corresponding multiple evaluation images; calculating image quality scores of the target coronary arteries in each evaluation image; when the target coronary arteries include at least two, acquiring weighting parameters corresponding to each target coronary artery; performing weighted calculation according to the image quality scores of the target coronary arteries and the corresponding weighting parameters to obtain the image quality score of each evaluation image; and determining the imaging phase of the cardiac blood vessel based on the quality scores of the multiple evaluation images. According to the application, the user can determine the optimal phase for reconstruction according to the use purpose, and the imaging quality of the target coronary artery is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of medical imaging, and in particular to a cardiac blood vessel imaging phase determination method and device, an electronic device and a storage medium. BACKGROUND

[0002] In medical image reconstruction, multi-phase data can be acquired by a CT scanner, and image reconstruction is performed according to the multi-phase data to obtain multiple images of a patient's part. Automatic determination of a suitable phase point for image reconstruction can improve the quality of the target reconstructed image.

[0003] In cardiac scanning, the quality of the coronary vessels determines the quality of the cardiac image. Due to the movement of the heart during CT scanning, the optimal imaging phase of the cardiac blood vessels needs to be selected for image reconstruction. However, the motion patterns of different coronary vessels are different. According to clinical experience, the best image quality of different coronary vessels (such as the left main stem, the left anterior descending branch, the left circumflex branch, and the right coronary artery) is usually at different phases. Therefore, the phase obtained by using a global image best phase algorithm may not be the best phase for the target coronary vessel, which affects the imaging quality of the target coronary vessel. SUMMARY

[0004] Embodiments of the present application provide a cardiac blood vessel imaging phase determination method, device, electronic device and storage medium to at least solve the problem of inappropriate imaging phase of a target coronary vessel in related technologies, which affects the imaging quality.

[0005] In a first aspect, embodiments of the present application provide a cardiac blood vessel imaging phase determination method, comprising:

[0006] obtaining multiple phase images of a cardiac blood vessel;

[0007] performing target coronary vessel extraction based on the multiple phase images respectively to obtain corresponding multiple evaluation images;

[0008] calculating an image quality score of the corresponding target coronary vessel in each evaluation image;

[0009] when the target coronary vessel includes at least two, obtaining a weighting parameter corresponding to each target coronary vessel;

[0010] performing weighted calculation according to the image quality score of each target coronary vessel and the corresponding weighting parameter to obtain an image quality score of each evaluation image;

[0011] determining an imaging phase of the cardiac blood vessel based on the quality scores of the multiple evaluation images.

[0012] In some embodiments, the target coronary extraction based on the plurality of phase images respectively to obtain the corresponding plurality of evaluation images comprises:

[0013] Target coronary positioning is performed on the phase images;

[0014] An image segmentation threshold is determined;

[0015] Target coronary extraction is performed on the phase images after target coronary positioning according to the image segmentation threshold to obtain the corresponding evaluation images;

[0016] The above steps are repeated to obtain a plurality of evaluation images.

[0017] In some embodiments, the determination of the image segmentation threshold comprises:

[0018] The pixel value or CT value of the phase image is obtained;

[0019] A preset subdivision parameter corresponding to the target coronary is obtained;

[0020] The image segmentation threshold of the phase image is calculated according to the pixel value or CT value and the subdivision parameter.

[0021] In some embodiments, before the target coronary extraction on the phase images after target coronary positioning according to the image segmentation threshold to obtain the corresponding evaluation images, at least one of the following processing steps is further included:

[0022] A segmentation center is marked according to the position of the target coronary in the phase image, and a segmentation region is demarcated based on the segmentation center to segment the phase image based on the segmentation region;

[0023] The phase image is reconstructed through image interpolation operation;

[0024] Morphological operation is performed on the phase image to weaken the image background.

[0025] In some embodiments, the calculation of the image quality score of the target coronary in each evaluation image comprises the calculation of the image quality score of the target coronary in a single evaluation image:

[0026] The weight parameters corresponding to a plurality of different quality evaluation indicators are determined;

[0027] A plurality of quality evaluation indicators of the target coronary in the evaluation image are calculated;

[0028] The quality score of the target coronary artery in the to-be-evaluated image is obtained based on the plurality of quality evaluation indexes corresponding to the target coronary artery and the weight parameters corresponding to each quality evaluation index.

[0029] In some embodiments, determining the imaging phase of the heart blood vessel based on the quality scores of the plurality of to-be-evaluated images comprises:

[0030] A plurality of to-be-evaluated images in a preset single-phase sliding window are obtained, and common images and non-common images corresponding to each phase are screened out;

[0031] The image quality score of the corresponding phase is calculated based on the image quality scores of the common images and the non-common images in the single phase and the number of image layers of the single phase;

[0032] The image quality scores of the plurality of to-be-evaluated images in each phase in the single-phase sliding window are calculated repeatedly by the above steps, and the imaging phase of the heart blood vessel is determined according to the image quality scores.

[0033] In some embodiments, the calculation of the image quality score of the corresponding phase based on the image quality scores of the common images and the non-common images in the single phase and the number of image layers of the single phase comprises:

[0034] The average score of the corresponding phase is calculated based on the quality score of the common image in the single phase in the single-phase sliding window and the number of layers of each phase with the common image;

[0035] The inter-phase deviation score is calculated according to the number of layers of the common image in the single phase and the average number of layers of the phases;

[0036] The intra-phase deviation score is determined based on the quality score of the non-common image in the single phase;

[0037] The image quality score of the corresponding phase is calculated by weighted calculation based on the average score, the inter-phase deviation score and the intra-phase deviation score.

[0038] In a second aspect, the embodiments of the present application provide a heart blood vessel imaging phase determination device, comprising:

[0039] A phase image acquisition unit is configured to acquire a plurality of phase images of a heart blood vessel;

[0040] A to-be-evaluated image acquisition unit is configured to extract a target coronary artery based on the plurality of phase images respectively, and obtain a plurality of corresponding to-be-evaluated images;

[0041] The first quality score calculation unit is configured to calculate the image quality score of each of the target coronary arteries in each of the images to be evaluated.

[0042] The weighting parameter acquisition unit is configured to acquire the weighting parameter corresponding to each of the target coronary arteries when the target coronary arteries include at least two.

[0043] The second quality score calculation unit is configured to calculate the image quality score of each of the images to be evaluated by weighting the image quality scores of the target coronary arteries and the corresponding weighting parameters.

[0044] The imaging phase determination unit is configured to determine the imaging phase of the heart blood vessels based on the quality scores of the images to be evaluated.

[0045] In a third aspect, an electronic device is provided, which includes a memory and a processor. The memory stores a computer program. The processor is configured to execute the computer program to perform the method for determining the imaging phase of the heart blood vessels as described above.

[0046] In a fourth aspect, a computer readable storage medium is provided, which stores a computer program. The computer program is executed by a processor to implement the method for determining the imaging phase of the heart blood vessels as described above.

[0047] Compared with the related art, the method for determining the imaging phase of the heart blood vessels provided by the embodiments of the present application extracts the target coronary arteries based on the plurality of phase images respectively to obtain the corresponding plurality of images to be evaluated, calculates the image quality score of each of the images to be evaluated by weighting the image quality scores of the target coronary arteries and the corresponding weighting parameters, and determines the imaging phase of the heart blood vessels based on the quality scores of the images to be evaluated. Therefore, the user can customize the weighting parameters of different target coronary arteries as needed, so that the image quality score ranking corresponding to the global coronary quality that is most suitable for the user can be obtained after the image quality score of each of the images to be evaluated is determined based on the weighting parameters of the target coronary arteries. The imaging phase of the heart blood vessels is determined based on the quality scores of the images to be evaluated, so that the user can determine the optimal phase for reconstruction according to the use purpose, and the imaging quality of the target coronary arteries is improved.

[0048] The details of one or more embodiments of the present application are presented in the following drawings and description to make other features, objects and advantages of the present application more apparent. BRIEF DESCRIPTION OF DRAWINGS

[0049] The accompanying drawings illustrated herein are used to provide further understanding of the present application, and constitute a part of the present application. The illustrative embodiments of the present application and their description serve to explain the present application, and do not constitute improper limitations on the present application. In the drawings:

[0050] Figure 1 is a flowchart of a method for determining an imaging phase of a heart blood vessel according to an embodiment of the present application;

[0051] Figure 2 is a flowchart of a method for extracting a target coronary artery based on a plurality of phase images according to an embodiment of the present application;

[0052] Figure 3 is a diagram of a result of positioning a target coronary artery according to an embodiment of the present application;

[0053] Figure 4 is a flowchart of a method for calculating an image quality score of a target coronary artery in a single image to be evaluated according to an embodiment of the present application;

[0054] Figure 5 is a block diagram of a device for determining an imaging phase of a heart blood vessel according to an embodiment of the present application;

[0055] Figure 6 is a diagram of an electronic device according to an embodiment of the present application.

[0056] BRIEF DESCRIPTION OF DRAWINGS: 11, left coronary; 12, right coronary; 201, phase image acquisition unit; 202, image to be evaluated acquisition unit; 203, first quality score calculation unit; 204, weighting parameter acquisition unit; 205, second quality score calculation unit; 206, imaging phase determination unit; 30, bus; 31, processor; 32, memory; 33, communication interface. DETAILED DESCRIPTION

[0057] In order to make the objectives, technical solutions, and superiorities of the present application clearer, the following will describe and explain the present application in combination with the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and should not be used to limit the present application. Based on the embodiments provided by the present application, all other embodiments obtained by those of ordinary skill in the art without any creative effort should fall within the scope of protection of the present application.

[0058] Obviously, the drawings in the following description are only some examples or embodiments of the present application, and for those of ordinary skill in the art, the present application can be applied to other similar scenarios without any creative effort based on these drawings. In addition, it should be understood that although the efforts made in this development process can be complex and lengthy, for those of ordinary skill in the art related to the content disclosed in the present application, some design, manufacture or production changes based on the technical content disclosed in the present application are only routine technical means and should not be understood as insufficient disclosure of the present application.

[0059] Reference to an "embodiment" in this application means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the application. The appearances of the phrase in various places in the specification are not necessarily all referring to the same embodiment, nor are they necessarily mutually exclusive of other embodiments. It is expressly understood that the described embodiments of the application are merely example structures selected for the descriptiveness and that other embodiments can be equally selected.

[0060] Unless otherwise defined, technical terms and scientific terms used in the present application shall have the ordinary meanings as understood by one of ordinary skill in the art to which this application pertains. The terms "a", "an", "one", "this", and the like, as used in the application, do not denote a limitation of quantity and can be used to denote one or more than one. The terms "including", "containing", "having", and the like, as used in the application, are meant to be inclusive, unless otherwise noted. For example, a process, method, system, product, or apparatus that comprises a list of steps or modules (units) is not limited to the listed steps or units, but can also include other steps or units not listed, or can also include other steps or units inherent to the process, method, product, or apparatus. The terms "connected", "coupled", and the like, as used in the application, do not necessarily mean physically or mechanically connected, but can also include electrically connected, whether directly or indirectly. The term "plurality" means two or more. The term "and / or" describes an association relationship between associated objects, which means that there can be three relationships, for example, "A and / or B" can mean that A exists alone, A and B exist together, and B exists alone. The character " / " generally means that the associated objects before and after it have an "or" relationship. The terms "first", "second", "third", and the like, as used in the application, are merely to distinguish similar objects, and do not represent a specific order of the objects.

[0061] Computed Tomography (CT) is a scanning device that uses X-rays to scan a certain thickness from many directions of the human body, converts the attenuated X-rays into visible light by a detector, and then converts the visible light into an electrical signal. Finally, after analog-to-digital conversion of the electrical signal, a computer device is used to reconstruct the image to obtain the final CT image. When using CT to reconstruct the image of the heart, the clarity of the coronary blood vessel visualization is the key to determining the quality of the reconstructed image of the heart. Due to the physiological characteristics of the movement of the heart, the optimal phase data needs to be selected for image reconstruction when reconstructing the image.

[0062] With the rapid development of computer technology, the acquisition, processing, display and storage of medical images have been digitized, and the image data processed by doctors and the workload of reading films have increased exponentially. Controlling human factors and improving the quality of image acquisition and processing are the key to correct disease diagnosis.

[0063] The embodiment also provides a cardiac vessel imaging phase determination method. Figure 1 The flowchart of the cardiac vessel imaging phase determination method according to the embodiment of the application is shown in FIG. 1, which includes the following steps: Figure 1

[0064] In step S101, a plurality of phase images of a cardiac vessel are acquired.

[0065] In the embodiment, the CT continuously scans the scanned cardiac vessel in a period of time and obtains corresponding scanning data when performing normal scanning work. The plurality of to-be-evaluated images are acquired according to the scanning data.

[0066] Specifically, the multi-phase reconstruction can be performed according to preset reconstruction parameters, and the reconstruction parameters include a preset reconstruction center and a preset reconstruction range. The cardiac protocol general parameters can be used. Since the position of the coronary artery in the chest cavity is not fixed, the coronary artery has a curved shape. Therefore, the reconstruction range of the multi-phase image needs to contain the entire region to be scanned, that is, all the phases that can be reconstructed. Meanwhile, considering the balance between the resolution of the image and the calculation efficiency, the size of the reconstruction matrix and the field of view of the reconstruction should not be too small or too large.

[0067] In step S102, target coronary extraction is performed based on the plurality of phase images respectively, and a plurality of corresponding to-be-evaluated images are obtained.

[0068] In the embodiment, the target coronary extraction can be performed on the plurality of phase images based on a deep learning neural network model, which is not described herein. Optionally, when there are multiple target coronaries, each to-be-evaluated image contains multiple extracted target coronaries. Specifically, by optimizing the network structure of the deep learning neural type, the training feature parameters and the loss function, different extraction scales can be determined, such as extracting only the left coronary artery (hereinafter referred to as the left coronary) and the right coronary artery (hereinafter referred to as the right coronary), or extracting more subdivided 16-level coronaries. After the quality score evaluation of the to-be-evaluated images, the quality scores of different target coronaries can be obtained, so that the corresponding best imaging phase for imaging can be selected.

[0069] Of course, in other embodiments, the target coronary extraction can also be performed based on image processing, enhancement filtering and region growing, which is not limited herein.

[0070] ​Step S103, calculate the image quality score of each image corresponding to the target coronary artery.

[0071] Specifically, in the present embodiment, any one of the quality evaluation indexes such as inter-region contrast, intra-region uniformity, shape smoothness measure, and region shape area difference can be used to evaluate the quality of each image to be evaluated, or multiple evaluation indexes can be combined to make the evaluation result more comprehensive.

[0072] Step S104, when the target coronary artery includes at least two, obtain the weighting parameter corresponding to each target coronary artery.

[0073] In the present embodiment, the target coronary artery can be a large branch such as the left anterior descending branch, the circumflex branch, the right coronary, or at least two of the left anterior descending branch, the circumflex branch, and the right coronary. The present application is not limited thereto. When the target coronary artery includes at least two, the user can customize the weighting parameter of each different target coronary artery according to the need. For example, the importance of the quality of different coronary arteries in the overall cardiac vascular image quality can be determined according to the clinical need, and the weighting parameter of each different target coronary artery can be configured accordingly.

[0074] Step S105, according to the image quality score of each target coronary artery and the corresponding weighting parameter, the image quality score of each image to be evaluated is obtained by weighted calculation.

[0075] Specifically, after obtaining the weighting parameter corresponding to each target coronary artery, the image quality score of each target coronary artery is weighted and summed with the corresponding weighting parameter to obtain the image quality score of each image to be evaluated. The image quality score of the image to be evaluated is the global coronary image quality corresponding to the phase sorting that best meets the actual clinical needs of the user, so that the user can determine the best phase for reconstruction according to the use purpose, and the imaging quality of the target coronary artery is improved.

[0076] For example, when the image quality scores of multiple target coronary arteries in a certain image to be evaluated are X1, X2, and X3, and the corresponding weighting parameters are w1, w2, and w3, the image quality score of the image to be evaluated can be determined as Y=X1*w1+X2*w2+X3*w3.

[0077] Step S106, based on the quality scores of multiple images to be evaluated, determine the imaging phase of the cardiac blood vessels.

[0078] In the present embodiment, determining the imaging phase of the cardiac blood vessels based on the quality scores of multiple images to be evaluated means that the user can select different phase positions for reconstruction according to the use purpose of the phase image.

[0079] Specifically, the quality score ranking result can be determined based on the quality scores of the plurality of images to be evaluated, and the index corresponding to the ranking result represents the ranking of the coronary quality in the phase. The best phase is the phase with the first ranking, and the user can select the phase for reconstruction according to the user's needs. Specifically: (1) the best phase can be selected for reconstruction to diagnose the coronary artery; (2) the best phase can be selected for motion correction, and a non-optimal phase (the worst or an intermediate level phase) can be selected for motion correction, and the quality of the two corrected images is compared to evaluate the correction effect; (3) the best phase in the systole or the best phase in the diastole can be selected for reconstruction to observe the myocardium or the coronary artery in different phases; (4) the quality score of the image to be evaluated in the best phase without correction can be retrieved, and the information of the patient (age, gender, heart rate, disease, etc.) is combined to compare different patients to explore the correlation between different conditions and the coronary quality. In addition, the method can also be used for reference for the quality control of scanning and imaging of the heart blood vessels, and the present application is not limited herein.

[0080] In summary, the method for determining the imaging phase of the heart blood vessels provided by the embodiments of the present application can extract the target coronary arteries based on the plurality of phase images respectively to obtain a plurality of images to be evaluated, perform weighted calculation according to the image quality scores of the target coronary arteries and the corresponding weighted parameters to obtain the image quality score of each image to be evaluated, so that the user can customize the weighted parameters of different target coronary arteries as needed, and then determine the image quality score of each image to be evaluated based on the weighted parameters of the target coronary arteries to obtain an image quality score ranking corresponding to the global coronary quality most suitable for the user. The imaging phase of the heart blood vessels is determined based on the quality scores of the plurality of images to be evaluated, so that the user can determine the best phase for reconstruction according to the use purpose, and the imaging quality of the target coronary artery is improved.

[0081] The embodiments of the present application will be described and explained below through preferred embodiments.

[0082] As shown in Figure 2 In some embodiments, the target coronary extraction based on the plurality of phase images respectively to obtain a plurality of images to be evaluated includes:

[0083] In step S1021, the target coronary is positioned on the phase image.

[0084] In the embodiment, a network model such as Vnet can be used for target coronary positioning. Specifically, there are four layers of up-sampling and four layers of down-sampling, and the probability of whether a pixel point in a segmentation result judgment phase image or feature map is foreground or background is predicted. The difference between the gold standard and the residual error is obtained, and the residual error is optimized to train the entire model. The trained model is used for target coronary positioning. It can be understood that other deep learning models can also be used for target coronary positioning, and the model is not specifically limited in the application. For example, as shown in Figure 3 The Vnet network model is used for target coronary positioning, and left coronary 11 and right coronary 12 are obtained.

[0085] In step S1022, the image segmentation threshold is determined.

[0086] In the embodiment, the image segmentation threshold can be a preset value, or can be obtained by calculation, such as a preset multiple of the maximum value of each pixel point in each phase image. The image segmentation threshold can also be configured by using image processing methods, etc.

[0087] In some embodiments, the image segmentation threshold can be obtained by calculation, and the determination of the image segmentation threshold includes the following steps: obtaining the pixel value or CT value of the phase image, and obtaining the preset subdivision parameter corresponding to the target coronary; according to the pixel value or CT value and the subdivision parameter, the image segmentation threshold of the phase image is calculated, and the calculation formula is as follows:

[0088] TL = T x Q

[0089] Wherein, TL is the image segmentation threshold, T is the image segmentation reference, the pixel value of the phase image includes the pixel value of all pixel points of the phase image, and the CT value of the phase image includes the CT value of all pixel points of the phase image. The preset subdivision parameter Q of the target coronary is the reference basis of the image segmentation subdivision degree. The preset subdivision parameter Q can be configured as 0-1. Specifically, the artifact morphology extracted by different subdivision parameters Q is different. The lower the subdivision parameter Q, the more motion artifacts contained in the segmented image. When the subdivision program of the image is required to be high, the preset subdivision parameter Q is large, and vice versa, the preset subdivision parameter Q is small.

[0090] In some embodiments, the image segmentation reference T can be determined according to the pixel value, and the image segmentation threshold TL of the phase image can be calculated according to the image segmentation reference T and the subdivision parameter Q. Optionally, the image segmentation reference T can be the maximum value of the pixel value of all pixel points of the phase image.

[0091] In some embodiments, the image segmentation reference T can be determined according to the CT value, and the image segmentation threshold TL of the phase image can be calculated according to the image segmentation reference T and the subdivision parameter Q. Optionally, the image segmentation reference T can be the median of the CT value of all pixel points of the phase image.

[0092] It can be understood that the determination method of the image segmentation reference T is not limited to the above, and other image processing methods such as Otsu method or fixed parameter can also be used to obtain the image segmentation threshold.

[0093] It should be noted that before determining the image segmentation threshold, if the phase image is pre-processed (such as image enhancement), the pixel value of the pixel point in the phase image is different from the CT value of the pixel point in the phase image, and at this time, the image segmentation reference T can be determined by using the pixel value of the phase image after the image pre-processing or the CT value of the original phase image.

[0094] Through the above steps, the image segmentation threshold can be flexibly determined, and the to-be-evaluated image containing the target coronary artery and the motion artifact carried by the target coronary artery can be extracted, so that the image quality of the coronary artery can be more accurately evaluated on the coronary artery containing the artifact.

[0095] Step S1023: performing target coronary extraction on the phase image positioned by the target coronary according to the image segmentation threshold, to obtain a corresponding to-be-evaluated image.

[0096] Step S1024: repeating the above steps to obtain multiple to-be-evaluated images.

[0097] Specifically, in the embodiment, the step of obtaining a single to-be-evaluated image includes: taking the image with a gray value greater than the image segmentation threshold in the phase image as the corresponding to-be-evaluated image. Using multiple image segmentation thresholds to segment the phase image will obtain multiple image regions in the to-be-evaluated image.

[0098] Through the above steps, the to-be-evaluated image is obtained by performing target coronary extraction on the phase image positioned by the target coronary based on the image segmentation threshold, and the to-be-evaluated image can be accurately determined according to the image segmentation threshold, so that the determination of the optimal imaging phase based on the to-be-evaluated image after the target coronary extraction will be more accurate.

[0099] On the basis of the above embodiment, in some embodiments, before the step of performing target coronary extraction on the phase image positioned by the target coronary according to the image segmentation threshold to obtain a corresponding to-be-evaluated image, at least one of the following processing steps is further included:

[0100] Step S1022A, a segmentation center is marked according to the position of the target coronary artery in the phase image, and a segmentation region is demarcated based on the segmentation center, so as to segment the phase image based on the segmentation region.

[0101] Specifically, the phase image can be pre-segmented to obtain a sub-image based on the segmentation center and the demarcated segmentation region, and the segmentation center is taken as a center point, and N*N pixels are taken as the segmented sub-image. The calculation process is as follows:

[0102] I sub I (R1:R2, R3:R4), R2-R1=R4-R3=N-1

[0103] Wherein, I (R1:R2, R3:R4) represents the pixel index range; R1 and R2 are the starting point and the ending point of the row in the segmentation region respectively; R3 and R4 are the starting point and the ending point of the column in the segmentation region respectively; N is the number of pixels in the row or column; I sub is the size of the pixel block; and sub is the pixel index.

[0104] It can be understood that the size of N can be self-defined and configured, and the complete coronary artery in the phase image can be covered. The segmentation center can be determined by taking a single marked coronary artery as a reference, or can be determined by taking multiple coronary arteries as a reference. The present application does not limit this.

[0105] Through the above steps, the sub-image is obtained by segmenting the phase image based on the segmentation region, and the corresponding to-be-evaluated image is obtained by extracting the target coronary artery from the sub-image, which can reduce the calculation amount of the target coronary artery extraction process.

[0106] Step S1022B, the phase image is reconstructed by image interpolation operation. Specifically, before determining the image segmentation threshold, the phase image can be reconstructed by interpolation.

[0107] Through the above steps, the phase image is reconstructed by image interpolation operation, which can improve the resolution of the image, thereby improving the accuracy of calculating the blood vessel shape and the blood vessel edge. Of course, in other embodiments, if the resolution of the phase image is sufficient, the phase image does not need to be reconstructed by image interpolation operation.

[0108] Step S1022C, morphological operation is performed on the phase image to weaken the image background.

[0109] Specifically, tophat transformation can be performed on the phase image.

[0110] It should be noted that before the target coronary artery is extracted from the phase image positioned by the target coronary artery according to the image segmentation threshold to obtain the corresponding to-be-evaluated image, at least one of S1022A-S1022C in the above processing steps is further included, and when multiple processing steps in S1022A-S1022C are adopted, the number and order of each processing step adopted by the present application are not limited. Through the above steps, the image background is weakened, and the region where the target coronary artery is located is highlighted.

[0111] As shown in the above embodiments, in some of the embodiments, the calculating the image quality score corresponding to the target coronary artery in each to-be-evaluated image includes calculating the image quality score corresponding to the target coronary artery in a single to-be-evaluated image. Specifically as follows: Figure 4

[0112] For different coronary arteries or different positions of the same coronary artery, the morphological information corresponding to the target coronary artery is not the same. For example, the front end and the rear end of the target coronary artery are long strip-shaped, the roundness is small, the sharpness is large, and the trend is different from the middle end. The morphological information of the target coronary artery can be determined by the blood vessel area in the 2D image, the trend of the coronary artery, the shape, or the topological structure in 3D, and the present application is not limited here.

[0113] For the target coronary artery with different morphological information, the importance of different evaluation indexes is different. For example, for a long strip-shaped coronary artery (hereinafter referred to as a long coronary artery, Long CA), it is required to have a little smaller roundness and a little larger sharpness. Therefore, on the cross section of the long coronary artery, the weight of the regularity is lower than that of the edge sharpness. In this embodiment, when there are multiple evaluation indexes, the weight parameters of different evaluation indexes can be self-defined and configured, so that the evaluation result can reflect the importance of different evaluation indexes, and the evaluation result is more accurate.

[0114] Step S1032, calculating multiple quality evaluation indexes corresponding to the target coronary artery in the to-be-evaluated image.

[0115] In the evaluation of cardiac vascular images, there is a clinical evaluation standard: the ideal coronary artery should have a clear edge, the second is that the edge is slightly blurred without obvious artifacts, the passing grade is that the coronary artery contour is visible and some slight artifacts are allowed, and the blood vessel edge blur motion artifact is serious and the blood vessel contour disappears, which all belong to the situation that cannot be diagnosed. Combined with the above clinical evaluation standard, multiple evaluation indexes can be proposed.

[0116] In the evaluation of cardiac vascular images, there is a clinical evaluation standard: the ideal coronary artery should have a clear edge, the second is that the edge is slightly blurred without obvious artifacts, the passing grade is that the coronary artery contour is visible and some slight artifacts are allowed, and the blood vessel edge blur motion artifact is serious and the blood vessel contour disappears, which all belong to the situation that cannot be diagnosed. Combined with the above clinical evaluation standard, multiple evaluation indexes can be proposed.

[0117] ​In the embodiment, the quality evaluation indexes of the target coronary artery can be shape regularity and boundary sharpness, which are respectively used to measure whether the boundary of the region of interest is blurred and the strength of motion artifacts. The two standards are quantified as shape regularity (strength of artifacts) and boundary sharpness (degree of boundary sharpness) respectively. The shape regularity and the boundary sharpness can basically cover coronary artery images of various quality levels (including vessels with stents, calcified vessels, etc.), and have stronger versatility. Of course, in other embodiments, other evaluation indexes can be used to quantify the morphological features of the heart vessels, such as the proportion of low CT values in the blood vessels (the CT values of artifacts are generally lower than the CT values of angiographic contrast agents) or entropy. The present application does not limit this.

[0118] In step S1033, the quality score of the target coronary artery in each of the images to be evaluated is obtained by weighted calculation based on the plurality of quality evaluation indexes corresponding to the target coronary artery and the weight parameters corresponding to each quality evaluation index.

[0119] Specifically, after obtaining the weight parameters corresponding to each quality evaluation index, the numerical values of the quality evaluation indexes of the target coronary artery are weighted and summed with the corresponding weight parameters to obtain the image quality score of the target coronary artery in each of the images to be evaluated.

[0120] Exemplarily, the image quality score of the target coronary artery can be obtained by the following formula:

[0121]

[0122] Wherein, regularity is the shape regularity, sharpness is the boundary sharpness, factorS is the weight parameter corresponding to the shape regularity, factorSL is the weight parameter corresponding to the boundary sharpness, QuaIdx is the image quality score of the target coronary artery, if LongCA is false indicates when the target coronary artery is not a long coronary artery, if LongCA is true indicates when the target coronary artery is a long coronary artery, + is addition operation, and × is multiplication operation.

[0123] In addition, since the magnitudes of the shape regularity and the boundary sharpness are not consistent, it is necessary to pull the two metrics to a baseline, which can be achieved by weighting or normalization, and the present application does not limit this.

[0124] By the above steps, the weight parameters corresponding to different quality evaluation indexes are determined. Thus, the weight parameters of different evaluation indexes are self-defined, so that the evaluation results can reflect the importance of different evaluation indexes, and the evaluation results of the image quality score of the heart blood vessels are more accurate and reliable.

[0125] On the basis of the above embodiments, in some of the embodiments, the determining the imaging phase of the heart blood vessels based on the quality scores of the plurality of images to be evaluated comprises:

[0126] In step S1061, a plurality of images to be evaluated in a preset single-phase sliding window are obtained, and common images and non-common images corresponding to each phase are screened.

[0127] In the scanning process of the heart blood vessels, the phases collected contain multiple cardiac cycles. For example, in a retrospective heart scan, the collected phases contain the entire cardiac cycle (0% to 100%), and in a prospective scan, the cardiac cycle usually includes the systole and diastole. However, the motion of the coronary artery in the z direction is not consistent in different cardiac cycles such as systole and diastole. In a real clinical scenario, the evaluation of the images to be evaluated based on 2D cross-section needs to consider the inconsistency of the coronary artery in the z direction at different phases. Therefore, the calculation of the image quality score in a single image to be evaluated cannot be simply compared at the same z direction coordinate. The z direction refers to the body length direction.

[0128] In this embodiment, the phase sliding window refers to the phase range within a preset window width in the entire cardiac cycle. It is considered that the motion pattern of the coronary artery in the z direction is similar within a single phase sliding window. The length of the phase sliding window can be adaptively configured. For example, the length of the phase sliding window can be adaptively adjusted according to the heart rate. The length of the phase sliding window should not be too wide. Optionally, the range of the phase sliding window is between 10% and 20%.

[0129] In this embodiment, based on the plurality of images to be evaluated within the phase sliding window, common images and non-common images corresponding to each phase can be screened, wherein the common images corresponding to each phase are the images to be evaluated having the same image reconstruction layer in each phase, and the non-common images (extraslice) are other images to be evaluated except the common images.

[0130] In step S1062, the image quality score of the corresponding phase is calculated based on the image quality scores of the common images and the non-common images in a single phase, and the number of image layers of the single phase.

[0131] In the present embodiment, firstly, the average score of the corresponding phase is calculated based on the quality score of the common image in the single phase within the single phase sliding window and the number of layers of each phase with the common image. Specifically, the quality score average VQSAvg is determined based on the image quality score of each phase with the common image in the single phase within the single phase sliding window; the number of layers average AvgRange is determined based on the number of layers of each phase with the common image in the single phase within the single phase sliding window; and the product of the quality score average VQSAvg and the number of layers average AvgRange is determined as the average score of the corresponding phase.

[0132] Then, the inter-phase deviation score VQOffset is calculated according to the number of layers of the common image in the single phase and the number of layers average of each phase. Specifically, the difference between the number of layers of the common image in the single phase and the number of layers average of each phase is determined as the inter-phase deviation score VQOffset. It can be understood that the inter-phase deviation score can be used to represent the single phase deviation situation, and the calculation method is not limited thereto.

[0133] Next, the intra-phase deviation score VQExt is determined based on the quality score of the non-common image in the single phase. In some embodiments, the intra-phase deviation score VQExt can be the quality score of the non-common image. In other embodiments, the intra-phase deviation score VQExt can be the difference between the quality score of the non-common image and the quality score average, which is not limited herein.

[0134] Finally, the image quality score VQ of the corresponding phase is calculated by weighted calculation according to the average score, the inter-phase deviation score VQOffset and the intra-phase deviation score VQExt, and the specific calculation method is as follows:

[0135] VQ = VQSAvg x AvgRange + VQOffset x w1 + VQExt x w2

[0136] Wherein, w1 and w2 are weight coefficients, and the range is 0-1 respectively.

[0137] It can be understood that in other embodiments, the calculation method of the image quality score of the corresponding phase is not limited thereto, for example, when the image quality scores of the common image and the non-common image in the single phase are obtained, and the number of image layers of the single phase, the image quality score of the corresponding phase can be directly weighted and summed and normalized, which is not limited herein.

[0138] Step S1063, the image quality scores of the multiple images to be evaluated of each phase within the single phase sliding window are calculated by repeating the above steps, and the imaging phase of the cardiac blood vessel is determined according to the image quality scores.

[0139] In the embodiment, the image quality scores of the multiple to-be-evaluated images of each phase in the single phase sliding window are calculated by repeating the above steps, the image quality scores of the multiple to-be-evaluated images of each phase are sorted, and different phase positions in the sorting are selected for reconstruction according to the use purpose of the phase image. Through the above steps, the image quality scores of the multiple to-be-evaluated images of each phase are calculated by introducing the phase sliding window, the motion mode of the coronary artery in the z direction is considered to be similar in the single phase sliding window, so as to better match the variability of the coronary artery in the z direction at different phases, so that the calculation result of the image quality score is more accurate, and the reliability of the optimal imaging phase and the imaging quality of the target coronary artery are ensured.

[0140] It should be noted that the steps shown in the above flow or the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from here.

[0141] The embodiment also provides a cardiac vascular imaging phase determination apparatus, which is used to implement the above embodiments and preferred embodiments, and details are not repeated. As used below, the terms "module", "unit", "sub-unit" and the like can be a combination of software and / or hardware that implements a predetermined function. Although the apparatus described in the following embodiments is preferably implemented in software, hardware or a combination of software and hardware is also possible and contemplated.

[0142] Figure 5 is a structural block diagram of a cardiac vascular imaging phase determination apparatus according to an embodiment of the present application, as shown in Figure 5 The apparatus includes a phase image acquisition unit 201, a to-be-evaluated image acquisition unit 202, a first quality score calculation unit 203, a weighting parameter acquisition unit 204, a second quality score calculation unit 205, and an imaging phase determination unit 206.

[0143] The phase image acquisition unit 201 is configured to acquire multiple phase images of a cardiac blood vessel.

[0144] The to-be-evaluated image acquisition unit 202 is configured to extract a target coronary artery based on the multiple phase images respectively to obtain corresponding multiple to-be-evaluated images.

[0145] The first quality score calculation unit 203 is configured to calculate an image quality score of the target coronary artery in each to-be-evaluated image.

[0146] The weighting parameter acquisition unit 204 is configured to acquire a weighting parameter corresponding to each target coronary artery when the target coronary artery includes at least two.

[0147] The second quality score calculation unit 205 is configured to calculate the image quality score of each of the to-be-evaluated images according to the image quality score of each of the target coronary arteries and the corresponding weighting parameter;

[0148] The imaging phase determination unit 206 is configured to determine the imaging phase of the cardiac blood vessel based on the quality scores of the to-be-evaluated images.

[0149] In some embodiments, the to-be-evaluated image acquisition unit 202 comprises a target coronary artery positioning module, a threshold determination module, a target coronary artery extraction module, and a loop module.

[0150] The target coronary artery positioning module is configured to position the target coronary artery in the phase image.

[0151] The threshold determination module is configured to determine an image segmentation threshold.

[0152] The target coronary artery extraction module is configured to extract the target coronary artery from the phase image positioned with the target coronary artery according to the image segmentation threshold to obtain the corresponding to-be-evaluated image.

[0153] The loop module is configured to repeat the above steps to obtain multiple to-be-evaluated images.

[0154] In some embodiments, the threshold determination module comprises a pixel parameter acquisition module, a subdivision parameter acquisition module, and a segmentation threshold calculation module.

[0155] The pixel parameter acquisition module is configured to acquire a pixel value or a CT value of the phase image.

[0156] The subdivision parameter acquisition module is configured to acquire a preset subdivision parameter corresponding to the target coronary artery.

[0157] The segmentation threshold calculation module is configured to calculate the image segmentation threshold of the phase image according to the pixel value or the CT value and the subdivision parameter.

[0158] In some embodiments, the cardiac blood vessel imaging phase determination apparatus further comprises at least one of a first preprocessing module, a second preprocessing module, and a second preprocessing module.

[0159] The first preprocessing module is configured to mark a segmentation center according to the position of the target coronary artery in the phase image, and to demarcate a segmentation region based on the segmentation center, so as to segment the phase image based on the segmentation region.

[0160] The second preprocessing module is configured to reconstruct the phase image through image interpolation operation.

[0161] A third pre-processing module is configured to perform a morphological operation on the phase image to weaken the background of the image.

[0162] In some embodiments, the first quality score calculation unit 203 comprises a weight parameter determination module, an evaluation index calculation module, and a first quality score calculation module.

[0163] The weight parameter determination module is configured to determine weight parameters corresponding to a plurality of different quality evaluation indexes;

[0164] The evaluation index calculation module is configured to calculate a plurality of quality evaluation indexes corresponding to the target coronary artery in the image to be evaluated;

[0165] The first quality score calculation module is configured to perform weighted calculation based on the plurality of quality evaluation indexes corresponding to the target coronary artery and the weight parameters corresponding to each quality evaluation index, to obtain a quality score of the image to be evaluated corresponding to the target coronary artery.

[0166] In some embodiments, the imaging phase determination unit 206 comprises:

[0167] The image acquisition module is configured to acquire a plurality of images to be evaluated within a preset single phase sliding window, and screen to obtain common images and non-common images corresponding to each phase;

[0168] The image quality score calculation module is configured to calculate an image quality score of a corresponding phase based on image quality scores of the common images and the non-common images within a single phase, and the number of image layers of the single phase;

[0169] The phase determination module is configured to repeatedly calculate the image quality scores of the plurality of images to be evaluated of each phase within the single phase sliding window, and determine the imaging phase of the heart blood vessel according to the image quality scores.

[0170] In some embodiments, the image quality score calculation module comprises:

[0171] The average score calculation module is configured to calculate an average score of a corresponding phase based on the quality scores of the common images in a single phase within a single phase sliding window and the number of layers of each phase having common images;

[0172] The inter-phase deviation score calculation module is configured to calculate an inter-phase deviation score based on the number of layers of the common images in the single phase and the average number of layers of the phases;

[0173] The intra-phase deviation score calculation module is configured to determine an intra-phase deviation score based on the quality scores of the non-common images in a single phase;

[0174] The weighting calculation module is configured to perform weighting calculation according to the average score, the inter-phase deviation score and the intra-phase deviation score to obtain an image quality score corresponding to the phase.

[0175] It should be noted that each of the above modules can be a functional module or a program module, and can be implemented by software or hardware. For the module implemented by hardware, each of the above modules can be located in the same processor; or each of the above modules can also be located in different processors in any combination.

[0176] In addition, in combination with Figure 6 The cardiac vascular imaging phase determination method described in the embodiments of the present application can be implemented by an electronic device. Figure 6 The hardware structure of the electronic device according to the embodiments of the present application is shown in the figure.

[0177] The electronic device can include a processor 31 and a memory 32 storing computer program instructions.

[0178] Specifically, the processor 31 can include a central processing unit (CPU), or an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application.

[0179] The memory 32 can include a mass storage for data or instructions. By way of example and not limitation, the memory 32 can include a hard disk drive (HDD), floppy disk drive, solid state drive (SSD), flash memory, compact disk read only memory (CD-ROM), digital versatile disk (DVD), optical disk, tape, universal serial bus (USB) drive, or a combination of two or more of these. The memory 32 can be removable and / or non-removable (or fixed) as appropriate. The memory 32 can be internal or external to the data processing device as appropriate. In particular embodiments, the memory 32 is a non-volatile memory. In particular embodiments, the memory 32 includes read only memory (ROM) and random access memory (RAM). The ROM can be a mask programmed ROM, programmable ROM (PROM), erasable PROM (EPROM), electrically erasable PROM (EEPROM), electrically alterable ROM (EAROM), or FLASH, or a combination of two or more of these, as appropriate. The RAM can be static random access memory (SRAM) or dynamic random access memory (DRAM), which can be Fast Page Mode Dynamic Random Access Memory (FPMDRAM), Extended Data Output Dynamic Random Access Memory (EDODRAM), Synchronous Dynamic Random Access Memory (SDRAM), etc., as appropriate.

[0180] The memory 32 can be used to store or buffer various data files required for processing and / or communication, and possible computer program instructions executed by the processor 31.

[0181] The processor 31 reads and executes the computer program instructions stored in the memory 32 to implement any one of the above-mentioned embodiments of the cardiac vascular imaging phase determination method.

[0182] In some embodiments, the electronic device can further include a communication interface 33 and a bus 30. In which, as shown in the figure, the processor 31, the memory 32, the communication interface 33 are connected through the bus 30 and complete the communication between each other. Figure 6

[0183] The communication interface 33 is used to realize the communication between each module, device, unit and / or equipment in the embodiments of the present application. The communication interface 33 can also realize data communication with other components, such as: external devices, image / data acquisition devices, databases, external storage, and image / data processing workstations, etc.

[0184] ​Bus 30 includes hardware, software, or both, to couple components of the electronic device to each other and to couple components to other components in communication. Although bus 30 is shown in Figure 1 as a single bus, alternative embodiments include any number of buses. Bus 30 includes, for example, one or more of a data bus, an address bus, a control bus, an expansion bus, a local bus, a graphics bus, an Accelerated Graphics Port (AGP), an Extended Industry Standard Architecture (EISA) bus, a Front Side Bus (FSB), a Hyper Transport (HT) interconnect, an Industry Standard Architecture (ISA) bus, an InfiniBand (IB) interconnect, a Low Pin Count (LPC) bus, a memory bus, a Micro Channel Architecture (MCA) bus, a Peripheral Component Interconnect (PCI) bus, a PCI-Express (PCI-X) bus, a Serial Advanced Technology Attachment (SATA) bus, a Video Electronics Standards Association Local (VLB) bus, or any other suitable bus or interconnect, or a combination of two or more of these. Where appropriate, bus 30 can include one or more buses. Although the present embodiments describe and show a particular bus, the present embodiments contemplate any suitable bus or interconnect.

[0185] The electronic device can execute the cardiac vessel imaging phase determination method in the embodiments of the present application based on the obtained program instructions, thereby realizing the cardiac vessel imaging phase determination method in the embodiments of the present application in combination with Figure 1 The cardiac vessel imaging phase determination method described.

[0186] In addition, in combination with the cardiac vessel imaging phase determination method in the above embodiments, the embodiments of the present application can provide a computer readable storage medium to realize. The computer readable storage medium has computer program instructions stored thereon; the computer program instructions are executed by the processor to realize any one of the cardiac vessel imaging phase determination methods in the above embodiments.

[0187] Any combination of the technical features in the above-described embodiments can be made. For the sake of brevity, the foregoing description is not intended to be exhaustive or to limit the scope of the application to the precise form disclosed. Modifications and alterations can occur to others upon reading and understanding the disclosure. Accordingly, the scope of the application is intended to be limited only by the appended claims.

[0188] The above-described embodiments only express several implementation manners of the application, and the description is relatively specific and detailed, but it should not be understood as a limitation on the patent scope of the application. It should be pointed out that, for ordinary skilled persons in the art, several modifications and improvements can be made without departing from the concept of the application, and these all belong to the protection scope of the application. Therefore, the patent protection scope of the application should be subject to the appended claims.

Claims

1. A method for determining phase of cardiovascular imaging, characterized in that: include: Acquire multiple phase images of cardiac vessels; Extracting target coronary arteries based on the multiple phase images to obtain corresponding multiple images to be evaluated; Calculating an image quality score corresponding to the target coronary artery in each image to be evaluated; When the target coronary arteries include at least two, obtaining a weighted parameter corresponding to each of the target coronary arteries; Performing weighted calculation according to the image quality score of each target coronary artery and the corresponding weighting parameter to obtain an image quality score of each image to be evaluated; determining the imaging phase of the cardiovascular system based on the quality scores of the plurality of images to be evaluated; The quality scores of the plurality of images to be evaluated are used to determine a quality score ranking result, and an index corresponding to the ranking result represents a ranking of the coronary artery quality at a corresponding phase.

2. The cardiovascular imaging phase determination method according to claim 1, wherein: The target coronary artery is extracted based on the multiple phase images to obtain the corresponding multiple images to be evaluated, including: performing target coronary artery positioning on the phase image; Determine the image segmentation threshold; Extracting the target coronary artery from the phase image after the target coronary artery is located according to the image segmentation threshold to obtain a corresponding image to be evaluated; Repeat the above steps to obtain multiple images to be evaluated.

3. The cardiovascular imaging phase determination method according to claim 2, wherein: Determining the image segmentation threshold comprises: Acquiring pixel values ​​or CT values ​​of the phase image; Acquiring preset subdivision parameters corresponding to the target coronary artery; An image segmentation threshold of the phase image is calculated according to the pixel value or CT value and the subdivision parameter.

4. The cardiovascular imaging phase determination method according to claim 2, wherein: Before extracting the target coronary artery from the phase image after the target coronary artery is located according to the image segmentation threshold to obtain the corresponding image to be evaluated, at least one of the following processing steps is also included: marking a segmentation center according to a position of a target coronary artery in the phase image, and defining a segmentation region based on the segmentation center, so as to segment the phase image based on the segmentation region; reconstructing the phase image through image interpolation operation; A morphological operation is performed on the phase image to weaken the image background.

5. The cardiovascular imaging phase determination method according to claim 1, wherein: Calculating the image quality score corresponding to the target coronary artery in each image to be evaluated includes calculating the image quality score corresponding to the target coronary artery in a single image to be evaluated: Determine weight parameters corresponding to multiple different quality evaluation indicators; Calculating a plurality of quality evaluation indicators corresponding to the target coronary artery in the image to be evaluated; A weighted calculation is performed based on the multiple quality evaluation indicators corresponding to the target coronary artery and the weight parameters corresponding to each quality evaluation indicator to obtain a quality score corresponding to the target coronary artery in the image to be evaluated.

6. The cardiovascular imaging phase determination method according to claim 1, wherein: Determining the imaging phase of the cardiovascular system based on the quality scores of the plurality of images to be evaluated includes: Acquire multiple images to be evaluated within a preset single phase sliding window, and screen out the common images and non-common images corresponding to each time; Calculating an image quality score of a corresponding phase based on the image quality scores of the shared image and the non-shared image in a single phase and the number of image layers in the single phase; Repeat the above steps to calculate and obtain image quality scores of multiple images to be evaluated at each time phase within the single phase sliding window, and determine the imaging phase of the cardiovascular system based on the image quality scores.

7. The cardiovascular imaging phase determination method according to claim 6, wherein: The calculating the image quality score of the corresponding phase based on the image quality scores of the shared images and the non-shared images in the single phase and the number of image layers of the single phase includes: Based on the quality score of the shared image in a single phase within a single phase sliding window and the number of layers of each phase having the shared image, an average score of the corresponding phase is calculated; Calculating an inter-phase deviation score according to the number of layers of the common image in the single phase and the average number of layers of each phase; determining an intra-phase deviation score based on the quality scores of the non-shared images in a single phase; A weighted calculation is performed based on the average score, the inter-phase deviation score, and the intra-phase deviation score to obtain an image quality score for the corresponding phase.

8. A cardiovascular imaging phase determination device, characterized in that: include: A phase image acquisition unit, used for acquiring multiple phase images of cardiovascular vessels; an image acquisition unit for evaluation, configured to extract target coronary arteries based on the plurality of phase images to obtain a corresponding plurality of images for evaluation; a first quality score calculation unit, configured to calculate an image quality score corresponding to the target coronary artery in each image to be evaluated; a weighted parameter acquisition unit, configured to acquire a weighted parameter corresponding to each of the target coronary arteries when the target coronary arteries include at least two; a second quality score calculation unit, configured to perform weighted calculation based on the image quality score of each target coronary artery and a corresponding weighting parameter to obtain an image quality score of each image to be evaluated; an imaging phase determining unit, configured to determine the imaging phase of the cardiovascular system based on the quality scores of the plurality of images to be evaluated; The quality scores of the plurality of images to be evaluated are used to determine a quality score ranking result, and an index corresponding to the ranking result represents a ranking of the coronary artery quality at a corresponding phase.

9. An electronic device comprising a memory and a processor, characterized in that: A computer program is stored in the memory, and the processor is configured to run the computer program to perform the cardiovascular imaging phase determination method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the cardiovascular imaging phase determination method according to any one of claims 1 to 7 is implemented.

Images