A blood vessel contour generation method, a generation device, an equipment and a storage medium

Abstract

The application provides a blood vessel contour generation method, a blood vessel contour generation device, equipment and a storage medium. The method comprises the following steps: acquiring a blood vessel center line and an initial blood vessel contour of a target blood vessel in a coronary angiography image; determining a convex interval on the initial blood vessel contour based on the blood vessel center line and the initial blood vessel contour; for each convex interval, correcting a convex blood vessel contour corresponding to the convex interval based on blood vessel contours before and after the convex interval on the initial blood vessel contour to obtain a target blood vessel contour of the target blood vessel. In this way, the convex blood vessel contour on the initial blood vessel contour can be accurately located and corrected and optimized, so that a more accurate blood vessel contour can be quickly generated, which is helpful for subsequent blood vessel analysis.

Description

Technical Field

[0001] This application relates to the field of medical image processing technology, and in particular to a method, apparatus, device, and storage medium for generating blood vessel contours. Background Technology

[0002] Coronary angiography is a commonly used and effective method for diagnosing coronary atherosclerotic heart disease (CAD). It is a relatively safe and reliable invasive diagnostic technique, widely used in clinical practice and considered the "gold standard" for diagnosing CAD. The procedure involves injecting a contrast agent into the coronary arteries via a catheter, followed by X-ray imaging to obtain 2D coronary angiographic images. Accurate acquisition of the vessel contour from these images is a crucial step in subsequent analysis, affecting the assessment of vascular stenosis and the creation of a three-dimensional model, ultimately influencing the diagnostic results and treatment plan. Therefore, rapidly obtaining accurate vessel contours is a vital step in the vascular analysis of angiographic images.

[0003] Existing methods for generating blood vessel contours primarily rely on the gradient of the image, i.e., the boundary between the blood vessel and the image background; or they use deep learning to learn the contour size of each point on the centerline. The former is often affected by factors such as noise, while the latter often requires extensive data annotation and long-term model training and optimization. Furthermore, both methods share the drawback of difficulty in defining blood vessel branches or overlapping areas, resulting in bulging contours and inaccurate blood vessel outlines, which can negatively impact disease diagnosis and treatment. Summary of the Invention

[0004] In view of this, the purpose of this application is to provide a method, apparatus, device and storage medium for generating blood vessel contours, which can accurately locate the protruding blood vessel contours on the initial blood vessel contour and correct and optimize them, thereby quickly generating more accurate blood vessel contours, which is helpful for subsequent blood vessel analysis.

[0005] This application provides a method for generating blood vessel contours, the method comprising:

[0006] Obtain the centerline and initial contour of the target vessel in the coronary angiography image;

[0007] Based on the vessel centerline and the initial vessel contour, determine the convex regions on the initial vessel contour;

[0008] For each convex region, based on the vascular contours located before and after the convex region on the initial vascular contour, the convex vascular contour corresponding to the convex region is modified to obtain the target vascular contour of the target vascular.

[0009] Furthermore, the vessel centerline includes multiple vessel center points; the initial vessel contour includes an initial left vessel contour and an initial right vessel contour, the initial left vessel contour including multiple left vessel contour points, and the initial right vessel contour including multiple right vessel contour points; each vessel center point corresponds to one left vessel contour point and one right vessel contour point; after acquiring the vessel centerline and initial vessel contour of the target vessel in the coronary angiography image, and before determining the convex region on the initial vessel contour based on the vessel centerline and the initial vessel contour, the generation method further includes:

[0010] For each vessel center point, based on the coordinates of the left and right contour points of the vessel corresponding to that center point, determine whether the vessel center point has experienced contour disorder.

[0011] For each disordered vessel center point with a disordered contour, the coordinates of the vessel contour point on the side with the disordered contour corresponding to the disordered vessel center point are corrected based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel.

[0012] Furthermore, for each disordered vessel center point with a disordered contour, based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel, the coordinates of the vessel contour points on the side with the disordered contour corresponding to the disordered vessel center point are corrected, including:

[0013] For each disordered vessel center point with a distorted contour, based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel, the left contour direction, right contour direction, initial left contour length, and initial right contour length of the vessel center point are determined. Here, the contour direction on each side of the vessel refers to the direction from the vessel center point to the contour point on that side; the initial contour length on each side of the vessel refers to the distance from the vessel center point to the contour point on that side.

[0014] Compare the initial left contour length and the initial right contour length of the disordered blood vessel center point to determine whether the side where the contour is disordered at the center point of the disordered blood vessel is the left or right contour.

[0015] The direction of the blood vessel contour on the side where the center point of the disordered blood vessel is located is reversed, and the coordinates of the blood vessel contour point on the side where the contour is disordered are corrected based on the reversed contour direction.

[0016] Furthermore, the vessel centerline includes multiple vessel center points; the initial vessel contour includes multiple vessel contour points; each vessel center point corresponds to one vessel contour point; determining the convex region on the initial vessel contour based on the vessel centerline and the initial vessel contour includes:

[0017] For each vessel center point, the initial vessel contour length corresponding to the vessel center point is determined based on the vessel contour point corresponding to the vessel center point; wherein, the initial vessel contour length refers to the distance from the vessel center point to the corresponding vessel contour point.

[0018] Based on the initial blood vessel contour length corresponding to each blood vessel center point, convolution calculation is used to determine the minimum blood vessel contour length corresponding to each blood vessel center point.

[0019] Based on the minimum blood vessel contour length corresponding to each blood vessel center point and the initial blood vessel contour length, the convex interval on the initial blood vessel contour is determined.

[0020] Furthermore, the step of determining the minimum blood vessel contour length corresponding to each blood vessel center point using convolution calculation based on the initial blood vessel contour length corresponding to each blood vessel center point includes:

[0021] The initial blood vessel contour length sequence is composed of the initial blood vessel contour length corresponding to the center point of each blood vessel.

[0022] On the initial blood vessel contour length sequence, convolution calculation is performed using a preset convolution kernel through sliding multiplication. The minimum blood vessel contour length corresponding to each blood vessel center point is determined by the following formula:

[0023]

[0024] In the formula, i represents the i-th blood vessel center point on the blood vessel center line, j represents the j-th number in the preset convolution kernel, M is the length of the preset convolution kernel; f(i) represents the initial blood vessel contour length corresponding to the subscript i in the initial blood vessel contour length sequence; k(i) represents the value corresponding to the subscript i in the preset convolution kernel.

[0025] Furthermore, determining the convex region on the initial blood vessel contour based on the minimum blood vessel contour length corresponding to each blood vessel center point and the initial blood vessel contour length includes:

[0026] Determine the length ratio between the minimum vessel contour length and the initial vessel contour length corresponding to each vessel center point;

[0027] Based on the length ratio corresponding to each blood vessel center point and a preset ratio threshold, at least one protruding blood vessel center point is selected from the plurality of blood vessel center points.

[0028] At least one protruding vessel center point is divided into at least one group of protruding vessel center points, and the protruding interval on the initial vessel contour is determined according to the vessel center point number of each group of protruding vessel center points; wherein, the vessel center point number of each group of protruding vessel center points is consecutive.

[0029] Furthermore, for each protruding interval, based on the vessel contours located before and after the protruding interval on the initial vessel contour, the protruding vessel contour corresponding to the protruding interval is corrected to obtain the target vessel contour of the target vessel, including:

[0030] Determine the first initial blood vessel contour length corresponding to the blood vessel center point of the previous sequence number in the protrusion interval with the smallest sequence number, and determine the second initial blood vessel contour length corresponding to the blood vessel center point of the next sequence number in the protrusion interval with the largest sequence number.

[0031] Based on the first initial blood vessel contour length and the second initial blood vessel contour length, the initial blood vessel contour length corresponding to the center point of each protruding blood vessel in the protruding region is corrected using the following formula to obtain the target blood vessel contour length corresponding to the center point of each protruding blood vessel:

[0032]

[0033] In the formula, the protruding interval is represented as [i,j], where i and j represent the sequence number of the vessel center point, and n represents the sequence number of the protruding vessel center point in the protruding interval; r n This represents the target vessel contour length corresponding to the center point of the protruding vessel with sequence number n; r i-1 r represents the length of the first initial blood vessel contour. j+1 This represents the length of the second initial blood vessel contour;

[0034] For each protruding blood vessel center point, the corrected coordinates of the protruding blood vessel contour point corresponding to the center point of the protruding blood vessel are determined based on the target blood vessel contour length corresponding to the center point of the protruding blood vessel.

[0035] The target blood vessel contour is determined based on the corrected coordinates of the protruding blood vessel contour points and the coordinates of other blood vessel contour points.

[0036] This application embodiment also provides a blood vessel contour generation device, the generation device comprising:

[0037] The acquisition module is used to acquire the centerline and initial contour of the target blood vessel in the coronary angiography image;

[0038] The determination module is used to determine the protrusion region on the initial blood vessel contour based on the blood vessel centerline and the initial blood vessel contour;

[0039] The correction module is used to correct the protruding blood vessel contour corresponding to each protruding interval based on the blood vessel contours located before and after the protruding interval on the initial blood vessel contour, so as to obtain the target blood vessel contour of the target blood vessel.

[0040] This application also provides an electronic device, including: a processor, a memory, and a bus. The memory stores machine-readable instructions executable by the processor. When the electronic device is running, the processor communicates with the memory via the bus. When the machine-readable instructions are executed by the processor, the steps of the method for generating a blood vessel contour as described above are performed.

[0041] This application also provides a computer-readable storage medium storing a computer program that, when executed by a processor, performs the steps of the method for generating a blood vessel contour as described above.

[0042] This application provides a method, apparatus, device, and storage medium for generating a blood vessel contour, comprising: acquiring the vessel centerline and initial blood vessel contour of a target blood vessel in a coronary angiography image; determining a raised section on the initial blood vessel contour based on the vessel centerline and the initial blood vessel contour; and for each raised section, modifying the raised blood vessel contour corresponding to the raised section based on the blood vessel contours on the initial blood vessel contour located before and after the raised section, thereby obtaining the target blood vessel contour of the target blood vessel.

[0043] This allows for the accurate localization and correction of protruding vascular contours on the initial vascular profile, leading to the rapid generation of more accurate vascular contours, which is beneficial for subsequent vascular analysis.

[0044] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, preferred embodiments are described below in detail with reference to the accompanying drawings. Attached Figure Description

[0045] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0046] Figure 1 A flowchart illustrating a method for generating blood vessel contours according to an embodiment of this application is shown;

[0047] Figures 2(a) and 2(b) show schematic diagrams of a method for eliminating and correcting protruding blood vessel contours according to an embodiment of this application;

[0048] Figures 3(a) to 3(c) This illustration shows a schematic diagram of a method for correcting disordered blood vessel contours according to an embodiment of this application;

[0049] Figure 4 A schematic diagram of a blood vessel contour generation device provided in an embodiment of this application is shown.

[0050] Figure 5 A schematic diagram of the structure of an electronic device provided in an embodiment of this application is shown. Detailed Implementation

[0051] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. The components of the embodiments of this application described and shown in the accompanying drawings can generally be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of this application provided in the accompanying drawings is not intended to limit the scope of the claimed application, but merely represents selected embodiments of this application. Based on the embodiments of this application, every other embodiment obtained by those skilled in the art without inventive effort falls within the scope of protection of this application.

[0052] Research has found that existing methods for generating blood vessel contours mainly rely on the gradient of the reference image, i.e., the boundary between the blood vessel and the image background; or they use deep learning methods to learn the contour size of each point on the centerline. The former is often affected by factors such as noise, while the latter often requires a large amount of data annotation and a long period of model training and optimization. In addition, both of these methods share the drawback of difficulty in defining the branches or overlapping areas of blood vessels, resulting in bulging of the generated blood vessel contours, making the contours inaccurate and affecting the diagnosis and treatment of diseases.

[0053] Based on this, embodiments of this application provide a method, apparatus, device, and storage medium for generating blood vessel contours, which can accurately locate and correct the protruding blood vessel contours on the initial blood vessel contour, thereby quickly generating more accurate blood vessel contours, which is helpful for subsequent blood vessel analysis.

[0054] Please see Figure 1 , Figure 1 This is a flowchart illustrating a method for generating blood vessel contours provided in an embodiment of this application. Figure 1 As shown in the embodiments of this application, the generation method includes:

[0055] S101. Obtain the centerline and initial contour of the target vessel in the coronary angiography image.

[0056] S102. Based on the blood vessel centerline and the initial blood vessel contour, determine the convex region on the initial blood vessel contour.

[0057] The vessel centerline includes multiple vessel center points; the initial vessel contour includes multiple vessel contour points; each vessel center point corresponds to one vessel contour point.

[0058] S103. For each protruding interval, based on the blood vessel contours located before and after the protruding interval on the initial blood vessel contour, the protruding blood vessel contour corresponding to the protruding interval is modified to obtain the target blood vessel contour of the target blood vessel.

[0059] In this way, the protruding vascular contours on the initial vascular contour can be accurately located and corrected and optimized, thereby quickly generating more accurate vascular contours, which is helpful for subsequent vascular analysis.

[0060] In one possible implementation, after step S101 and before step S102, the generation method further includes:

[0061] S104. For each blood vessel center point, based on the coordinates of the left contour point and the right contour point of the blood vessel corresponding to the center point, determine whether the contour of the blood vessel center point has been disordered.

[0062] S105. For each disordered vessel center point with a disordered contour, based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel, the coordinates of the vessel contour points on the side with the disordered contour corresponding to the disordered vessel center point are corrected.

[0063] In this way, it is also possible to accurately locate disordered vascular contours on the initial vascular contour and correct and optimize the coordinates of the disordered side; furthermore, after correcting the disordered vascular contours, combined with the subsequent protruding vascular contour correction steps, a more accurate vascular contour can be quickly generated, which is more helpful for subsequent vascular analysis.

[0064] The following section will describe the implementation process of each step in detail with specific examples.

[0065] In one possible implementation, for step S101, the vessel centerline is a line used to identify the center of the vessel, usually represented by a line connecting spaced points and short line segments. The vessel contour is a line that forms the outer edge of the vessel along its natural extension direction. Exemplarily, both the vessel centerline and the initial vessel contour obtained in this step can be represented as a set of points in coordinate form; in specific implementations, the vessel centerline and initial vessel contour of the target vessel in the coronary angiography image can be obtained based on any method in the prior art, and this application is not limited thereto.

[0066] In one possible implementation, the vessel centerline includes multiple vessel center points; the initial vessel contour includes multiple vessel contour points; each vessel center point corresponds to one vessel contour point. Then step S102 may include the following steps:

[0067] S1021. For each blood vessel center point, determine the initial blood vessel contour length corresponding to the blood vessel center point based on the blood vessel contour point corresponding to the blood vessel center point; wherein, the initial blood vessel contour length refers to the distance from the blood vessel center point to the corresponding blood vessel contour point.

[0068] In this step, for each blood vessel center point, the initial blood vessel contour length corresponding to the blood vessel center point can be determined using the distance formula between the two points based on the coordinates of the blood vessel center point and the coordinates of the corresponding blood vessel contour point.

[0069] S1022. Based on the initial blood vessel contour length corresponding to each blood vessel center point, the minimum blood vessel contour length corresponding to each blood vessel center point is determined by convolution calculation.

[0070] In specific implementation, step S1022 may include:

[0071] Step 1: Use the initial vessel contour length corresponding to the center point of each vessel to form an initial vessel contour length sequence.

[0072] Step 2: On the initial blood vessel contour length sequence, perform convolution calculation using a preset convolution kernel through sliding multiplication, and determine the minimum blood vessel contour length corresponding to each blood vessel center point using the following formula:

[0073]

[0074] In the formula, i represents the i-th blood vessel center point on the blood vessel center line, j represents the j-th number in the preset convolution kernel, M is the length of the preset convolution kernel; f(i) represents the initial blood vessel contour length corresponding to the subscript i in the initial blood vessel contour length sequence; k(i) represents the value corresponding to the subscript i in the preset convolution kernel.

[0075] In one experiment, this application sets the convolution kernel to [0, 1 / 3, 1 / 3, 1 / 3, 0] with a length of 5. By using sliding multiplication, the minimum value near each point can be determined, thus obtaining the minimum vessel contour length corresponding to each vessel center point.

[0076] S1023. Based on the minimum blood vessel contour length corresponding to each blood vessel center point and the initial blood vessel contour length, determine the convex interval on the initial blood vessel contour.

[0077] In specific implementation, step S1023 may include:

[0078] Step 1: Determine the length ratio between the minimum blood vessel profile length and the initial blood vessel profile length corresponding to each blood vessel center point.

[0079] In this step, after determining the minimum blood vessel contour length corresponding to any blood vessel center point, the length ratio can be obtained by dividing the initial blood vessel contour length corresponding to the blood vessel center point by the minimum blood vessel contour length.

[0080] Step 2: Based on the length ratio corresponding to each blood vessel center point and a preset ratio threshold, at least one protruding blood vessel center point is selected from the multiple blood vessel center points.

[0081] The theoretical basis for this step is that, under normal circumstances, blood vessels should gradually narrow from beginning to end, and generally do not suddenly widen. Therefore, based on the length ratio and a preset ratio threshold, the center points of blood vessels with length ratios exceeding the preset ratio threshold can be selected as the center points of protruding blood vessels. The preset ratio threshold can be set according to the actual condition of the blood vessels; for example, the preset ratio threshold can be 1.1.

[0082] Step 3: Divide at least one protruding blood vessel center point into at least one group of protruding blood vessel center points, and determine the protruding interval on the initial blood vessel contour according to the blood vessel center point sequence number of each group of protruding blood vessel center points; wherein, the blood vessel center point sequence number of each group of protruding blood vessel center points is consecutive.

[0083] In this step, consecutive protruding vessel center points can be merged into a group based on their serial numbers to obtain a set of protruding vessel center points. For example, assuming the serial numbers of the protruding vessel center points are [1,2,4,5,6,9,10], they can be divided into three groups: [[1,2],[4,5,6],[9,10]]. Then, protrusion elimination correction is performed on each group.

[0084] In one possible implementation, step S103 may include:

[0085] S1031. Determine the first initial blood vessel contour length corresponding to the blood vessel center point of the preceding number of the smallest number in the protrusion interval, and determine the second initial blood vessel contour length corresponding to the blood vessel center point of the following number of the largest number in the protrusion interval.

[0086] S1032. Based on the first initial blood vessel contour length and the second initial blood vessel contour length, the initial blood vessel contour length corresponding to the center point of each protruding blood vessel in the protruding interval is corrected using the following formula to obtain the target blood vessel contour length corresponding to the center point of each protruding blood vessel:

[0087]

[0088] In the formula, the protruding interval is represented as [i,j], where i and j represent the sequence number of the vessel center point, and n represents the sequence number of the protruding vessel center point in the protruding interval; r n This represents the target vessel contour length corresponding to the center point of the protruding vessel with sequence number n; r i-1 r represents the length of the first initial blood vessel contour. j+1 This represents the length of the second initial blood vessel contour.

[0089] Here, considering that blood vessels should gradually taper from beginning to end under normal circumstances, the contour difference between the center point of the normal blood vessel in the previous sequence and the center point of the normal blood vessel in the next sequence is evenly distributed to the center point of each protruding blood vessel in the protruding interval. This makes the corrected blood vessel contour appear as a slow and gradual taper process under normal conditions, eliminating the abnormal protruding contour of the protruding interval.

[0090] S1033. For each protruding blood vessel center point, determine the corrected coordinates of the protruding blood vessel contour point corresponding to the center point of the protruding blood vessel based on the target blood vessel contour length corresponding to the center point of the protruding blood vessel.

[0091] In this step, after determining the target vessel contour length corresponding to the center point of the protruding vessel, the corrected coordinates of the protruding vessel contour point can be obtained by calculating based on the coordinates of the center point of the protruding vessel and the target vessel contour length. For example, it can be assumed that the corrected protruding vessel contour point still lies on the straight line connecting the center point of the protruding vessel to the original protruding vessel contour point; that is, it can be assumed that the direction between the center point of the protruding vessel and the corrected protruding vessel contour point remains unchanged. The coordinates of the point on this line, at a distance equal to the target vessel contour length from the center point of the protruding vessel, are calculated, and this coordinate point is determined as the corrected protruding vessel contour point, thus obtaining the corrected coordinates of the protruding vessel contour point.

[0092] S1034. Determine the target blood vessel contour based on the corrected coordinates of the protruding blood vessel contour points and the coordinates of other blood vessel contour points.

[0093] In this step, the target blood vessel contour after protrusion elimination can be determined based on the corrected coordinates of the protruding blood vessel contour points and the initial coordinates of other normal blood vessel contour points.

[0094] Please refer to Figures 2(a) and 2(b), which are schematic diagrams illustrating the elimination and correction of protruding blood vessel contours according to embodiments of this application. The process of eliminating and correcting protruding blood vessel contours will be described in detail below with reference to Figures 2(a) and 2(b).

[0095] Studies have shown that, as illustrated in Figure 2(a), the initial vascular contour is prone to bulges at overlapping or branching points. It should be noted that in planar coronary angiography images, the vascular contour can be divided into a left and a right contour. Both the left and right contours can be eliminated using the generation method provided in this application. Taking the right contour in Figure 2(a) as an example, a pre-defined convolution kernel is used to perform convolution calculations via sliding multiplication on the right contour. The minimum vascular contour length corresponding to each vascular center point is determined, resulting in the vascular contour shown by the dashed line in Figure 2(a).

[0096] Then, the convex interval can be determined by the length ratio between the minimum blood vessel contour length and the initial blood vessel contour length. For example, the convex blood vessel contour corresponding to a certain convex interval is shown in Figure 2(b), that is, [2,3,4] is a convex interval. Find the previous normal point and the next normal point of this convex interval, which are points 1 and 5 respectively. Then, the target blood vessel contour length from point 2 to point 4 can be determined by the above formula. Then, the corrected blood vessel contour is obtained according to the coordinates of the blood vessel center point 2 to the blood vessel center point 4, as shown by the dashed line in Figure 2(b). Among them, the optimized target blood vessel contour length of point 2 can be expressed as r2=r1-(r1+r5) / 4; the optimized target blood vessel contour length of point 3 can be expressed as r3=r1-2*(r1+r5) / 4.

[0097] Furthermore, this embodiment of the application also takes into account the problem of disordered blood vessel contours. Therefore, after step S101 and before step S102, the generation method further includes steps S104 and S105 for identifying and correcting disordered contours.

[0098] Please see Figures 3(a) to 3(c) , Figures 3(a) to 3(c) This is a schematic diagram illustrating the correction of disordered blood vessel contours provided in an embodiment of this application. The following will be combined with... Figures 3(a) to 3(c) This section details the process of correcting disordered blood vessel contours.

[0099] In one possible implementation, the vessel centerline includes multiple vessel center points; the initial vessel contour includes an initial left vessel contour and an initial right vessel contour, the initial left vessel contour including multiple left vessel contour points, and the initial right vessel contour including multiple right vessel contour points; each vessel center point corresponds to one left vessel contour point and one right vessel contour point. Then, in step S104, for each vessel center point, based on the coordinates of the corresponding left and right vessel contour points, it is determined whether the vessel center point has experienced contour distortion.

[0100] In practice, firstly, the coordinates of the left and right contour points corresponding to the center point of each blood vessel can be subtracted from the coordinates of the center point of the blood vessel to obtain the direction of the left and right contours of the blood vessel; the calculation formula for this process can be expressed as:

[0101]

[0102]

[0103] In the formula, i = 1, 2, ..., N, where N is the number of blood vessel center points; c i Indicates the coordinates of the center point of the i-th blood vessel; This represents the coordinates of the left contour point of the blood vessel corresponding to the center point of the i-th blood vessel. This represents the coordinates of the right contour point of the blood vessel corresponding to the center point of the i-th blood vessel. Indicates the direction of the left contour of the blood vessel at the center point of the i-th blood vessel; This indicates the direction of the right contour of the blood vessel at the center point of the i-th blood vessel.

[0104] Secondly, the angle between the directions of the left and right contours of the blood vessel can be used to determine whether contour misalignment has occurred at the center point of the vessel. The formula for calculating the angle between the left and right contours of the blood vessel can be expressed as:

[0105]

[0106] In the formula, ∝ i This represents the angle between the left and right contour directions of the blood vessel at the center point of the i-th blood vessel.

[0107] Then, according to ∝ i The relationship with a preset angle threshold can determine whether contour disorder has occurred at the center point of the i-th blood vessel. For example, if the included angle is less than 90°, it indicates that contour disorder has occurred. As shown in Figure 3(a), the left arrow in Figure 3(a) represents the normal left and right contour direction, and the right arrow in Figure 3(a) represents the left and right contour direction when contour disorder has occurred.

[0108] It should be noted that, in the specific calculation, the embodiment of this application adopts the batch matrix calculation method shown in Figure 3(b), which can calculate all blood vessel center points and blood vessel contour points at one time without the need for iterative calculation, thus speeding up the algorithm.

[0109] In one possible implementation, step S105 may include:

[0110] S1051. For each disordered vessel center point with a disordered contour, based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel, determine the left contour direction, right contour direction, initial left contour length, and initial right contour length of the vessel center point; wherein, the contour direction of each side of the vessel refers to the direction from the vessel center point to the contour point of that side; the initial contour length of each side of the vessel refers to the distance from the vessel center point to the contour point of that side.

[0111] In this step, the directions of the left and right contours of the blood vessel can be calculated using the same method as in step S104, and will not be repeated here. The initial lengths of the left and right contours of the blood vessel can also be determined using the formula for the distance from the contour point to the center point of the blood vessel. The formula can be expressed as:

[0112]

[0113] In the formula, l i v represents the initial contour length corresponding to the center point of the i-th blood vessel (applicable to both left and right sides); i (x) and v i (y) represents the x and y coordinates corresponding to the contour direction vector.

[0114] S1052. Compare the initial left contour length and the initial right contour length of the disordered blood vessel center point to determine whether the side where the contour of the disordered blood vessel center point is disordered is the left or right contour.

[0115] In this step, the initial left contour length and the initial right contour length of the disordered blood vessel center point are compared. If the initial left contour length is greater than the initial right contour length, it means that the left contour of the blood vessel crosses the right contour, that is, the left side of the blood vessel center point is disordered.

[0116] S1053. Invert the direction of the blood vessel contour on the side where the center point of the disordered blood vessel is located, and correct the coordinates of the blood vessel contour point on the side where the contour is disordered based on the inverted contour direction.

[0117] In this step, the direction of the blood vessel contour on the side with the disordered contour can be reversed (i.e., rotated 180°), and then the coordinates of the reversed contour direction and the center point of the disordered blood vessel are summed to obtain the coordinates of the corrected blood vessel contour point on the side with the disordered contour. This process can be expressed by the formula:

[0118]

[0119] In the formula, This represents the corrected coordinates of the left contour point of the i-th blood vessel center point after the disordered contour has been corrected; v i This indicates the direction of the blood vessel contour on the side where the contour is disordered at the center point of the i-th blood vessel.

[0120] Figure 3(c) shows a schematic diagram of the vascular contour after correction of the disordered contour. As can be seen from Figure 3(c), using the disordered contour correction method provided in this embodiment, the corrected vascular contour may exhibit abnormal protrusions. This protrusion may be due to overlapping or branching points of the vessels, or it may be caused by the correction of the original disordered contour. Therefore, after correcting the disordered contour, the corrected vascular contour can be redefined as the initial vascular contour, and then the protruding vascular contour can be eliminated and optimized using the methods in S102 and S103. The combination of these two correction methods—first correcting the disordered contour locally, then eliminating the protrusion globally—optimizes the vascular contour from both local and global perspectives, thereby quickly generating a more accurate vascular contour, which is beneficial for subsequent vascular analysis.

[0121] Please see Figure 4 , Figure 4 This is a schematic diagram of a blood vessel contour generation device provided in an embodiment of this application. Figure 4 As shown, the generating apparatus 400 includes:

[0122] The acquisition module 410 is used to acquire the center line and initial contour of the target blood vessel in the coronary angiography image.

[0123] The determining module 420 is used to determine the protrusion region on the initial blood vessel contour based on the blood vessel centerline and the initial blood vessel contour;

[0124] The correction module 430 is used to correct the protruding blood vessel contour corresponding to each protruding interval based on the blood vessel contours located before and after the protruding interval on the initial blood vessel contour, so as to obtain the target blood vessel contour of the target blood vessel.

[0125] Furthermore, the vessel centerline includes multiple vessel center points; the initial vessel contour includes an initial left vessel contour and an initial right vessel contour, the initial left vessel contour including multiple left vessel contour points, and the initial right vessel contour including multiple right vessel contour points; each vessel center point corresponds to one left vessel contour point and one right vessel contour point; after acquiring the vessel centerline and initial vessel contour of the target vessel in the coronary angiography image, and before determining the convex region on the initial vessel contour based on the vessel centerline and the initial vessel contour, the correction module 430 is further configured to:

[0126] For each vessel center point, based on the coordinates of the left and right contour points of the vessel corresponding to that center point, determine whether the vessel center point has experienced contour disorder.

[0127] For each disordered vessel center point with a disordered contour, the coordinates of the vessel contour point on the side with the disordered contour corresponding to the disordered vessel center point are corrected based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel.

[0128] Furthermore, when the correction module 430 corrects the coordinates of the vessel contour points on the side of the vessel with contour disorder corresponding to each disordered vessel center point based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right vessel contour points, the correction module 430 is used to:

[0129] For each disordered vessel center point with a distorted contour, based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel, the left contour direction, right contour direction, initial left contour length, and initial right contour length of the vessel center point are determined. Here, the contour direction on each side of the vessel refers to the direction from the vessel center point to the contour point on that side; the initial contour length on each side of the vessel refers to the distance from the vessel center point to the contour point on that side.

[0130] Compare the initial left contour length and the initial right contour length of the disordered blood vessel center point to determine whether the side where the contour is disordered at the center point of the disordered blood vessel is the left or right contour.

[0131] The direction of the blood vessel contour on the side where the center point of the disordered blood vessel is located is reversed, and the coordinates of the blood vessel contour point on the side where the contour is disordered are corrected based on the reversed contour direction.

[0132] Furthermore, the vessel centerline includes multiple vessel center points; the initial vessel contour includes multiple vessel contour points; each vessel center point corresponds to one vessel contour point; when the determining module 420 determines the convex region on the initial vessel contour based on the vessel centerline and the initial vessel contour, the determining module 420 is used to:

[0133] For each vessel center point, the initial vessel contour length corresponding to the vessel center point is determined based on the vessel contour point corresponding to the vessel center point; wherein, the initial vessel contour length refers to the distance from the vessel center point to the corresponding vessel contour point.

[0134] Based on the initial blood vessel contour length corresponding to each blood vessel center point, convolution calculation is used to determine the minimum blood vessel contour length corresponding to each blood vessel center point.

[0135] Based on the minimum blood vessel contour length corresponding to each blood vessel center point and the initial blood vessel contour length, the convex interval on the initial blood vessel contour is determined.

[0136] Furthermore, when determining the minimum blood vessel contour length corresponding to each blood vessel center point using convolution calculation based on the initial blood vessel contour length corresponding to each blood vessel center point, the determining module 420 is used to:

[0137] The initial blood vessel contour length sequence is composed of the initial blood vessel contour length corresponding to the center point of each blood vessel.

[0138] On the initial blood vessel contour length sequence, convolution calculation is performed using a preset convolution kernel through sliding multiplication. The minimum blood vessel contour length corresponding to each blood vessel center point is determined by the following formula:

[0139]

[0140] In the formula, i represents the i-th blood vessel center point on the blood vessel center line, j represents the j-th number in the preset convolution kernel, M is the length of the preset convolution kernel; f(i) represents the initial blood vessel contour length corresponding to the subscript i in the initial blood vessel contour length sequence; k(i) represents the value corresponding to the subscript i in the preset convolution kernel.

[0141] Furthermore, when determining the convex region on the initial blood vessel contour based on the minimum blood vessel contour length corresponding to each blood vessel center point and the initial blood vessel contour length, the determining module 420 is used to:

[0142] Determine the length ratio between the minimum vessel contour length and the initial vessel contour length corresponding to each vessel center point;

[0143] Based on the length ratio corresponding to each blood vessel center point and a preset ratio threshold, at least one protruding blood vessel center point is selected from the plurality of blood vessel center points.

[0144] At least one protruding vessel center point is divided into at least one group of protruding vessel center points, and the protruding interval on the initial vessel contour is determined according to the vessel center point number of each group of protruding vessel center points; wherein, the vessel center point number of each group of protruding vessel center points is consecutive.

[0145] Furthermore, when the correction module 430 corrects the protruding vessel contour corresponding to each protruding interval based on the vessel contours located before and after the protruding interval on the initial vessel contour to obtain the target vessel contour, the correction module 430 is used to:

[0146] Determine the first initial blood vessel contour length corresponding to the blood vessel center point of the previous sequence number in the protrusion interval with the smallest sequence number, and determine the second initial blood vessel contour length corresponding to the blood vessel center point of the next sequence number in the protrusion interval with the largest sequence number.

[0147] Based on the first initial blood vessel contour length and the second initial blood vessel contour length, the initial blood vessel contour length corresponding to the center point of each protruding blood vessel in the protruding region is corrected using the following formula to obtain the target blood vessel contour length corresponding to the center point of each protruding blood vessel:

[0148]

[0149] In the formula, the protruding interval is represented as [i,j], where i and j represent the sequence number of the vessel center point, and n represents the sequence number of the protruding vessel center point in the protruding interval; r n This represents the target vessel contour length corresponding to the center point of the protruding vessel with sequence number n; r i-1 r represents the length of the first initial blood vessel contour. j+1 This represents the length of the second initial blood vessel contour;

[0150] For each protruding blood vessel center point, the corrected coordinates of the protruding blood vessel contour point corresponding to the center point of the protruding blood vessel are determined based on the target blood vessel contour length corresponding to the center point of the protruding blood vessel.

[0151] The target blood vessel contour is determined based on the corrected coordinates of the protruding blood vessel contour points and the coordinates of other blood vessel contour points.

[0152] Please see Figure 5 , Figure 5 This is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. Figure 5As shown, the electronic device 500 includes a processor 510, a memory 520, and a bus 530.

[0153] The memory 520 stores machine-readable instructions executable by the processor 510. When the electronic device 500 is running, the processor 510 and the memory 520 communicate via the bus 530. When the machine-readable instructions are executed by the processor 510, they can perform the operations described above. Figure 1 The steps of the method for generating a blood vessel contour in the method embodiment shown in Figure 3 can be found in the method embodiment for specific implementation, and will not be repeated here.

[0154] This application also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, can perform the above-described actions. Figure 1 The steps of the method for generating a blood vessel contour as shown in Figure 3 can be found in the method embodiment, and will not be repeated here.

[0155] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0156] In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods can be implemented in other ways. The apparatus embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. Furthermore, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the shown or discussed mutual couplings, direct couplings, or communication connections may be through some communication interfaces; indirect couplings or communication connections between devices or units may be electrical, mechanical, or other forms.

[0157] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0158] In addition, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0159] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a processor-executable, non-volatile, computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0160] Finally, it should be noted that the above-described embodiments are merely specific implementations of this application, used to illustrate the technical solutions of this application, and not to limit them. The scope of protection of this application is not limited thereto. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments, or make equivalent substitutions for some of the technical features, within the scope of the technology disclosed in this application. Such modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application, and should all be covered within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for generating blood vessel contours, characterized in that, The generation method includes: The centerline and initial contour of the target vessel in the coronary angiography image are obtained; the centerline includes multiple vessel center points; the initial vessel contour includes an initial left contour and an initial right contour, the initial left contour includes multiple left contour points, and the initial right contour includes multiple right contour points; each vessel center point corresponds to one left contour point and one right contour point. Based on the vessel centerline and the initial vessel contour, the convex regions on the initial vessel contour are determined, including: For each vessel center point, the initial vessel contour length corresponding to the vessel center point is determined based on the vessel contour point corresponding to the vessel center point; wherein, the initial vessel contour length refers to the distance from the vessel center point to the corresponding vessel contour point. Based on the initial blood vessel contour length corresponding to each blood vessel center point, convolution calculation is used to determine the minimum blood vessel contour length corresponding to each blood vessel center point. Determine the length ratio between the minimum vessel contour length and the initial vessel contour length corresponding to each vessel center point; Based on the length ratio corresponding to each blood vessel center point and a preset ratio threshold, at least one protruding blood vessel center point is selected from the plurality of blood vessel center points. At least one protruding vessel center point is divided into at least one group of protruding vessel center points, and the protruding interval on the initial vessel contour is determined according to the vessel center point number of each group of protruding vessel center points; wherein, the vessel center point number of each group of protruding vessel center points is consecutive. For each convex region, based on the vessel contours before and after the convex region on the initial vessel contour, the convex vessel contour corresponding to the convex region is corrected to obtain the target vessel contour of the target vessel, including: Determine the first initial blood vessel contour length corresponding to the blood vessel center point of the previous sequence number in the protrusion interval with the smallest sequence number, and determine the second initial blood vessel contour length corresponding to the blood vessel center point of the next sequence number in the protrusion interval with the largest sequence number. Based on the first initial blood vessel contour length and the second initial blood vessel contour length, the initial blood vessel contour length corresponding to the center point of each protruding blood vessel in the protruding region is corrected using the following formula to obtain the target blood vessel contour length corresponding to the center point of each protruding blood vessel: In the formula, the convex region is represented as , and Indicates the sequence number of the center point of the blood vessel. This indicates the sequence number of the center point of the protruding blood vessel in the protruding area; Indicates the sequence number is The length of the target blood vessel contour corresponding to the center point of the protruding blood vessel; This represents the length of the first initial blood vessel contour; This represents the length of the second initial blood vessel contour; For each protruding blood vessel center point, the corrected coordinates of the protruding blood vessel contour point corresponding to the center point of the protruding blood vessel are determined based on the target blood vessel contour length corresponding to the center point of the protruding blood vessel. The target blood vessel contour is determined based on the corrected coordinates of the protruding blood vessel contour points and the coordinates of other blood vessel contour points. After acquiring the vessel centerline and initial vessel contour of the target vessel in the coronary angiography image, and before determining the convex region on the initial vessel contour based on the vessel centerline and the initial vessel contour, the generation method further includes: For each vessel center point, based on the coordinates of the left and right contour points of the vessel corresponding to that center point, determine whether the vessel center point has experienced contour disorder. For each disordered vessel center point with a disordered contour, the coordinates of the vessel contour point on the side with the disordered contour corresponding to the disordered vessel center point are corrected based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel.

2. The generation method according to claim 1, characterized in that, For each disordered vessel center point with a disordered contour, based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel, the coordinates of the vessel contour points on the side of the disordered vessel center point are corrected, including: For each disordered vessel center point with a distorted contour, based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel, the left contour direction, right contour direction, initial left contour length, and initial right contour length of the vessel center point are determined. Here, the contour direction on each side of the vessel refers to the direction from the vessel center point to the contour point on that side; the initial contour length on each side of the vessel refers to the distance from the vessel center point to the contour point on that side. Compare the initial left contour length and the initial right contour length of the disordered blood vessel center point to determine whether the side where the contour is disordered at the center point of the disordered blood vessel is the left or right contour. The direction of the blood vessel contour on the side where the center point of the disordered blood vessel is located is reversed, and the coordinates of the blood vessel contour point on the side where the contour is disordered are corrected based on the reversed contour direction.

3. The generation method according to claim 1, characterized in that, The step of determining the minimum blood vessel contour length corresponding to each blood vessel center point using convolution calculation based on the initial blood vessel contour length corresponding to each blood vessel center point includes: The initial blood vessel contour length sequence is composed of the initial blood vessel contour length corresponding to the center point of each blood vessel. On the initial blood vessel contour length sequence, convolution calculation is performed using a preset convolution kernel through sliding multiplication. The minimum blood vessel contour length corresponding to each blood vessel center point is determined by the following formula: In the formula, Indicates the first on the central line of the blood vessel The center point of each blood vessel Represents the first in the preset convolution kernel The number of kernels, M is the length of the preset convolutional kernel; Indicates the subscript in the initial blood vessel contour length sequence The corresponding initial blood vessel contour length; Indicates the subscript in the preset convolution kernel The corresponding numerical value.

4. A device for generating blood vessel contours, characterized in that, The generating apparatus includes: The acquisition module is used to acquire the vessel centerline and initial vessel contour of the target vessel in the coronary angiography image; the vessel centerline includes multiple vessel center points; the initial vessel contour includes an initial left vessel contour and an initial right vessel contour, the initial left vessel contour includes multiple left vessel contour points, and the initial right vessel contour includes multiple right vessel contour points; each vessel center point corresponds to one left vessel contour point and one right vessel contour point. The determination module is used to determine the convex regions on the initial blood vessel contour based on the blood vessel centerline and the initial blood vessel contour, including: For each vessel center point, the initial vessel contour length corresponding to the vessel center point is determined based on the vessel contour point corresponding to the vessel center point; wherein, the initial vessel contour length refers to the distance from the vessel center point to the corresponding vessel contour point. Based on the initial blood vessel contour length corresponding to each blood vessel center point, convolution calculation is used to determine the minimum blood vessel contour length corresponding to each blood vessel center point. Determine the length ratio between the minimum vessel contour length and the initial vessel contour length corresponding to each vessel center point; Based on the length ratio corresponding to each blood vessel center point and a preset ratio threshold, at least one protruding blood vessel center point is selected from the plurality of blood vessel center points. At least one protruding vessel center point is divided into at least one group of protruding vessel center points, and the protruding interval on the initial vessel contour is determined according to the vessel center point number of each group of protruding vessel center points; wherein, the vessel center point number of each group of protruding vessel center points is consecutive. The correction module is used to correct the protruding vessel contour corresponding to each protruding interval based on the vessel contours located before and after the protruding interval on the initial vessel contour, to obtain the target vessel contour of the target vessel, including: Determine the first initial blood vessel contour length corresponding to the blood vessel center point of the previous sequence number in the protrusion interval with the smallest sequence number, and determine the second initial blood vessel contour length corresponding to the blood vessel center point of the next sequence number in the protrusion interval with the largest sequence number. Based on the first initial blood vessel contour length and the second initial blood vessel contour length, the initial blood vessel contour length corresponding to the center point of each protruding blood vessel in the protruding region is corrected using the following formula to obtain the target blood vessel contour length corresponding to the center point of each protruding blood vessel: In the formula, the convex region is represented as , and Indicates the sequence number of the center point of the blood vessel. This indicates the sequence number of the center point of the protruding blood vessel in the protruding area; Indicates the sequence number is The length of the target blood vessel contour corresponding to the center point of the protruding blood vessel; This represents the length of the first initial blood vessel contour; This represents the length of the second initial blood vessel contour; For each protruding blood vessel center point, the corrected coordinates of the protruding blood vessel contour point corresponding to the center point of the protruding blood vessel are determined based on the target blood vessel contour length corresponding to the center point of the protruding blood vessel. The target blood vessel contour is determined based on the corrected coordinates of the protruding blood vessel contour points and the coordinates of other blood vessel contour points. After acquiring the vessel centerline and initial vessel contour of the target vessel in the coronary angiography image, and before determining the convex region on the initial vessel contour based on the vessel centerline and the initial vessel contour, the correction module is further configured to: For each vessel center point, based on the coordinates of the left and right contour points of the vessel corresponding to that center point, determine whether the vessel center point has experienced contour disorder. For each disordered vessel center point with a disordered contour, the coordinates of the vessel contour point on the side with the disordered contour corresponding to the disordered vessel center point are corrected based on the coordinates of the disordered vessel center point and the coordinates of the corresponding left and right contour points of the vessel.

5. An electronic device, characterized in that, include: The device includes a processor, a memory, and a bus. The memory stores machine-readable instructions executable by the processor. When the electronic device is running, the processor communicates with the memory via the bus. The machine-readable instructions are executed by the processor to perform the steps of a method for generating a blood vessel contour as described in any one of claims 1 to 3.

6. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, performs the steps of a method for generating a blood vessel contour as described in any one of claims 1 to 3.

Images