Blood vessel contour splicing method and device, electronic equipment and storage medium

Abstract

The invention provides a blood vessel contour splicing method and device, electronic equipment and a storage medium, and the method comprises the steps: determining a relative space relation between a to-be-spliced branch blood vessel and a main and branch blood vessel according to the branch contour information of the to-be-spliced branch blood vessel and the main and branch contour information of the main and branch blood vessel; determining a target contour line in the main and branch blood vessels according to the relative spatial relationship; determining a target correction strategy of the to-be-spliced branch blood vessel according to the minimum distance between the to-be-spliced branch blood vessel and the target contour line; the to-be-spliced branch blood vessels are processed according to the target correction strategy, and an initial blood vessel contour map of the processed to-be-spliced branch blood vessels coincided with the target contour line is obtained; and identifying spliced cross points in the initial blood vessel contour map, and cutting off redundant contour lines exceeding the spliced cross points to obtain a target blood vessel contour map. In this way, the branch blood vessel and the main branch blood vessel are accurately spliced, and therefore the accurate blood vessel contour can be obtained.

Description

Technical Field

[0001] The present application relates to the technical field of medical image processing, and in particular, to a method, device, electronic device, and storage medium for splicing vascular contours. Background Art

[0002] The incidence rate of coronary heart disease in China has been remaining high, and coronary angiography interventional therapy has become one of the standard treatment methods for this disease. Among them, QFR (Quantitative Flow Ratio) is an advanced technology based on coronary angiography images, which can quickly calculate FFR (Fractional Flow Reserve) to assist PCI (Percutaneous Coronary Intervention) surgery. This technology can quickly calculate functional parameters such as pressure drop and FFR value at the lesion site, as well as structural indexes such as lesion length, proximal and distal diameters of blood vessels, and stenosis rate through three-dimensional reconstruction of blood vessels and hydrodynamic analysis. Therefore, it can well assist doctors to more quickly and accurately identify clinically significant vascular stenosis, so as to formulate a more precise interventional treatment plan.

[0003] However, the vascular contours determined by existing contour segmentation methods usually have the situation that the branches and the main branch do not completely coincide. And the accurate vascular contour is the key to calculating the above important parameters in the above technology. If the branches and the main branch do not completely coincide, it will affect the shunt of the main branch, resulting in inaccurate simulation results. Therefore, it is very important for the branch contour to be accurately generated and spliced onto the main branch, both in terms of visualization and simulation algorithm calculation. Summary of the Invention

[0004] In view of this, the purpose of the present application is to provide a method, device, electronic device, and storage medium for splicing vascular contours, which can obtain an accurate vascular contour by accurately splicing the branch blood vessel and the main blood vessel.

[0005] The embodiment of the present application provides a method for splicing vascular contours, and the splicing method includes:

[0006] Determine the branch contour information of the branch blood vessel to be spliced and the main branch contour information of the main blood vessel;

[0007] According to the branch contour information and the main branch contour information, determine the relative spatial relationship between the branch blood vessel to be spliced and the main blood vessel;

[0008] According to the relative spatial relationship, determine the target contour line in the main blood vessel;

[0009] According to the minimum distance between the branch blood vessel to be spliced and the target contour line, determine the target correction strategy of the branch blood vessel to be spliced;

[0010] Process the to-be-spliced branch vessel according to the target correction strategy to obtain an initial vessel contour map in which the processed to-be-spliced branch vessel coincides with the target contour line;

[0011] Identify the splicing intersection points in the initial vessel contour map, and truncate the redundant contour lines beyond the splicing intersection points to obtain the target vessel contour map.

[0012] Optionally, determine the target correction strategy for the to-be-spliced branch vessel according to the minimum distance between the to-be-spliced branch vessel and the target contour line, including:

[0013] Determine whether the minimum distance between the to-be-spliced branch vessel and the target contour line is less than a preset distance;

[0014] If it is less, the target correction strategy is to perform an outer extension process on the to-be-spliced branch vessel according to the minimum distance;

[0015] If it is not less, the target correction strategy is to update the to-be-spliced branch vessel according to the original medical image and perform an outer extension process on the updated to-be-spliced branch vessel.

[0016] Optionally, when the minimum distance between the to-be-spliced branch vessel and the target contour line is less than the preset distance, the process of processing the to-be-spliced branch vessel according to the target correction strategy to obtain an initial vessel contour map in which the processed to-be-spliced branch vessel coincides with the target contour line includes:

[0017] For each branch contour line of the to-be-spliced branch vessel, determine the splicing starting point of this branch contour line;

[0018] Determine the extension distance according to the shortest distance between the splicing starting point of this branch contour line and the target contour line;

[0019] According to the branch contour information of this branch contour line, starting from the splicing starting point of this branch contour line, use the interpolation method to extend this branch contour line by the extension distance;

[0020] When all branch contour lines are fully extended, obtain the initial vessel contour map.

[0021] Optionally, when the minimum distance between the to-be-spliced branch vessel and the target contour line is not less than the preset distance, the process of processing the to-be-spliced branch vessel according to the target correction strategy to obtain an initial vessel contour map in which the processed to-be-spliced branch vessel coincides with the target contour line includes:

[0022] Obtain the original medical image including the to-be-spliced branch vessel and the main branch vessel;

[0023] Perform skeletonization on the original medical image to determine the centerline of the entire blood vessel;

[0024] Extract the centerline corresponding to the branch blood vessel to be spliced from the centerline of the entire blood vessel to obtain the centerline of the branch to be processed;

[0025] According to the centerline of the entire blood vessel, use the depth - first traversal algorithm to complete the centerline of the branch to be processed, and determine the complete branch centerline corresponding to the branch blood vessel to be spliced that intersects the centerline of the main branch blood vessel;

[0026] According to the complete branch centerline and the original medical image, update the branch contour information of the branch blood vessel to be spliced to obtain the updated branch blood vessel to be spliced;

[0027] According to the shortest distance between the updated branch blood vessel to be spliced and the target contour line, use the interpolation method to perform an outer extension process on the updated branch blood vessel to be spliced to obtain the initial blood vessel contour map.

[0028] Optionally, the identifying the splicing intersection points in the initial blood vessel contour map and truncating the redundant contour lines beyond the splicing intersection points to obtain the target blood vessel contour map includes:

[0029] For each branch contour line in the processed branch blood vessel to be spliced in the initial blood vessel contour map, determine the extended contour point set of this branch contour line; wherein, the extended end point in the extended contour point set is the other end point of this branch contour line;

[0030] According to the extended contour point set of this branch contour line and the main support point set of the target contour line, determine the shortest distance between this branch contour line and the target contour line;

[0031] Determine the extended contour point in the shortest distance as the splicing intersection point between this branch contour line and the target contour line;

[0032] Truncate the contour line between the splicing intersection point and the extended end point on this branch contour line to obtain the target blood vessel contour map.

[0033] Optionally, the determining the relative spatial relationship between the branch blood vessel to be spliced and the main branch blood vessel according to the branch contour information and the main branch contour information includes:

[0034] According to the branch contour information and the main branch contour information, determine the interval distance between each main branch contour line in the main branch blood vessel and the branch blood vessel to be spliced;

[0035] According to the two interval distances, determine the relative spatial relationship between the branch blood vessel to be spliced and the main branch blood vessel.

[0036] Optionally, determining the target contour line in the main branch vessel according to the relative spatial relationship includes:

[0037] Determining the main branch contour line on the same side as the branch vessel to be spliced as the target contour line.

[0038] The embodiment of the present application also provides a splicing device for blood vessel contours, and the splicing device includes:

[0039] A first determination module, configured to determine the branch contour information of the branch vessel to be spliced and the main branch contour information of the main branch vessel;

[0040] A second determination module, configured to determine the relative spatial relationship between the branch vessel to be spliced and the main branch vessel according to the branch contour information and the main branch contour information;

[0041] A third determination module, configured to determine the target contour line in the main branch vessel according to the relative spatial relationship;

[0042] A fourth determination module, configured to determine the target correction strategy of the branch vessel to be spliced according to the minimum distance between the branch vessel to be spliced and the target contour line;

[0043] A processing module, configured to process the branch vessel to be spliced according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch vessel to be spliced coincides with the target contour line;

[0044] A truncation module, configured to identify the splicing intersection point in the initial blood vessel contour map and truncate the redundant contour line beyond the splicing intersection point to obtain a target blood vessel contour map.

[0045] Optionally, when the fourth determination module is used to determine the target correction strategy of the branch vessel to be spliced according to the minimum distance between the branch vessel to be spliced and the target contour line, the fourth determination module is used for:

[0046] Determining whether the minimum distance between the branch vessel to be spliced and the target contour line is less than a preset distance;

[0047] If it is less, the target correction strategy is to perform an outer extension process on the branch vessel to be spliced according to the minimum distance;

[0048] If it is not less, the target correction strategy is to update the branch vessel to be spliced according to the original medical image and perform an outer extension process on the updated branch vessel to be spliced.

[0049] Optionally, when the minimum distance between the branch blood vessel to be spliced and the target contour line is less than a preset distance, when the processing module is used to process the branch blood vessel to be spliced according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessel to be spliced coincides with the target contour line, the processing module is used to:

[0050] For each branch contour line of the branch blood vessel to be spliced, determine the splicing starting point of this branch contour line;

[0051] Determine the extension distance according to the shortest distance between the splicing starting point of this branch contour line and the target contour line;

[0052] According to the branch contour information of this branch contour line, starting from the splicing starting point of this branch contour line, use the interpolation method to extend this branch contour line by the extension distance;

[0053] When all branch contour lines are fully extended, an initial blood vessel contour map is obtained.

[0054] Optionally, when the minimum distance between the branch blood vessel to be spliced and the target contour line is not less than the preset distance, when the processing module is used to process the branch blood vessel to be spliced according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessel to be spliced coincides with the target contour line, the processing module is used to:

[0055] Obtain the original medical image including the branch blood vessel to be spliced and the main blood vessel;

[0056] Perform skeletonization processing on the original medical image to determine the full blood vessel centerline;

[0057] Extract the centerline corresponding to the branch blood vessel to be spliced from the full blood vessel centerline to obtain the branch centerline to be processed;

[0058] According to the full blood vessel centerline, use the depth-first traversal algorithm to perform complement processing on the branch centerline to be processed, and determine the complete branch centerline corresponding to the branch blood vessel to be spliced that intersects the centerline of the main blood vessel;

[0059] According to the complete branch centerline and the original medical image, update the branch contour information of the branch blood vessel to be spliced to obtain the updated branch blood vessel to be spliced;

[0060] According to the shortest distance between the updated branch blood vessel to be spliced and the target contour line, use the interpolation method to perform outer extension processing on the updated branch blood vessel to be spliced to obtain an initial blood vessel contour map.

[0061] Optionally, when the truncation module is used to identify the splicing intersection points in the initial blood vessel contour diagram and truncate the redundant contour lines beyond the splicing intersection points to obtain the target blood vessel contour diagram, the truncation module is used for:

[0062] For each branch contour line in the processed branch blood vessel to be spliced in the initial blood vessel contour diagram, determine the extended contour point set of this branch contour line; wherein, the extended end point in the extended contour point set is the other end point of this branch contour line;

[0063] According to the extended contour point set of this branch contour line and the main fulcrum set of the target contour line, determine the shortest distance between this branch contour line and the target contour line;

[0064] Determine the extended contour point in the shortest distance as the splicing intersection point between this branch contour line and the target contour line;

[0065] Truncate the contour line between the splicing intersection point and the extended end point on this branch contour line to obtain the target blood vessel contour diagram.

[0066] Optionally, when the second determination module is used to determine the relative spatial relationship between the branch blood vessel to be spliced and the main blood vessel according to the branch contour information and the main branch contour information, the second determination module is used for:

[0067] According to the branch contour information and the main branch contour information, determine the interval distance between each main branch contour line in the main blood vessel and the branch blood vessel to be spliced;

[0068] According to the two interval distances, determine the relative spatial relationship between the branch blood vessel to be spliced and the main blood vessel.

[0069] Optionally, when the third determination module is used to determine the target contour line in the main blood vessel according to the relative spatial relationship, the third determination module is used for:

[0070] Determine the main branch contour line on the same side as the branch blood vessel to be spliced as the target contour line.

[0071] An embodiment of the present application further 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 runs, the processor communicates with the memory through the bus. When the machine-readable instructions are executed by the processor, the steps of the splicing method as described above are executed.

[0072] An embodiment of the present application further provides a computer-readable storage medium. A computer program is stored on the computer-readable storage medium. When the computer program is run by a processor, the steps of the splicing method as described above are executed.

[0073] A method, device, electronic device, and storage medium for splicing blood vessel contours provided by an embodiment of the present application. The splicing method includes: determining branch contour information of a branch blood vessel to be spliced and main branch contour information of a main branch blood vessel; determining a relative spatial relationship between the branch blood vessel to be spliced and the main branch blood vessel according to the branch contour information and the main branch contour information; determining a target contour line in the main branch blood vessel according to the relative spatial relationship; determining a target correction strategy for the branch blood vessel to be spliced according to a minimum distance between the branch blood vessel to be spliced and the target contour line; processing the branch blood vessel to be spliced according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessel to be spliced coincides with the target contour line; identifying a splicing intersection point in the initial blood vessel contour map and truncating redundant contour lines beyond the splicing intersection point to obtain a target blood vessel contour map.

[0074] In this way, when generating blood vessel contours, the present application takes into account the connectivity between branch blood vessels and main branch blood vessels, and splices unconnected branch blood vessels and main branch blood vessels by combining segmentation methods. At the same time, during the splicing process, different splicing methods are used for branch blood vessels with different distances from the main branch blood vessel, so as to take into account the splicing accuracy of branch blood vessels far from the main branch, and further ensure the accuracy of the splicing result, the accuracy of the simulation result, and make the present solution have a wider applicability.

[0075] To make the above objects, features, and advantages of the present application more obvious and understandable, the following specifically gives preferred embodiments and, in conjunction with the accompanying drawings, makes a detailed description as follows. BRIEF DESCRIPTION OF THE DRAWINGS

[0076] To more clearly illustrate the technical solutions of the embodiments of the present application, the following briefly introduces the drawings required for the embodiments. It should be understood that the following drawings only show some embodiments of the present application, and therefore should not be regarded as limiting the scope. For those of ordinary skill in the art, without creative efforts, other related drawings can also be obtained based on these drawings.

[0077] Figure 1 A schematic diagram of a blood vessel contour generation result provided by the present application;

[0078] Figure 2 A flowchart of a method for splicing blood vessel contours provided by an embodiment of the present application;

[0079] Figure 3 A schematic diagram of a blood vessel centerline provided by the present application;

[0080] Figure 4 An example diagram of the determination process of a target blood vessel contour map provided by the present application;

[0081] Figure 5 This is a schematic structural diagram of a splicing device for a blood vessel contour provided by an embodiment of the present application;

[0082] Figure 6 This is a schematic structural diagram of an electronic device provided by an embodiment of the present application. Detailed implementation manners

[0083] To make the objectives, technical solutions, and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application. Apparently, the described embodiments are only some of the embodiments of the present application, rather than all the embodiments. Usually, the components of the embodiments of the present application described and illustrated in the accompanying drawings here can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present application provided in the accompanying drawings is not intended to limit the scope of the present application to be protected, but only represents the selected embodiments of the present application. Based on the embodiments of the present application, every other embodiment obtained by those skilled in the art without creative efforts belongs to the scope of protection of the present application.

[0084] The incidence rate of coronary heart disease in our country has been remaining high, and coronary angiography interventional therapy has become one of the standard treatment methods for this disease. Among them, QFR (Quantitative Flow Ratio) is an advanced technology based on coronary angiography images, which can quickly calculate FFR (Fractional Flow Reserve) to assist PCI (Percutaneous Coronary Intervention) surgery. This technology can quickly calculate functional parameters such as pressure drop and FFR value at the lesion site, as well as structural indexes such as lesion length, proximal and distal vessel diameters, and stenosis rate through three-dimensional reconstruction of blood vessels and hydrodynamic analysis. Therefore, it can well assist doctors to more quickly and accurately identify clinically significant vascular stenosis, so as to formulate a more precise interventional treatment plan.

[0085] However, the blood vessel contours determined by the existing contour segmentation methods usually have the situation that the branches and the main branch do not completely coincide. For example, please refer to Figure 1 , Figure 1 This is a schematic diagram of a blood vessel contour generation result provided by the present application. As Figure 1As shown in the left figure, when the doctor wants to analyze the branch indicated by the red arrow but the software does not provide this branch, this branch is newly added. However, as seen in the right figure, the branch added through the prior art has a certain distance from the main branch contour and does not join the main branch (in addition, there may also be a situation where the branch passes through the main branch contour. In short, it cannot be accurately joined to the main branch). In the QFR technology, an accurate blood vessel contour is the key to calculating important parameters such as FFR. If the branch does not completely coincide with the main branch, it will affect the shunt of the main branch, resulting in inaccurate simulation results. Therefore, it is very important that the branch contour can be accurately generated and joined to the main branch, both in terms of visualization and simulation algorithm calculation.

[0086] Based on this, the embodiments of the present application provide a method, device, electronic device, and storage medium for splicing blood vessel contours. By accurately splicing the branch blood vessel and the main blood vessel, an accurate blood vessel contour can be obtained.

[0087] Please refer to Figure 2 , Figure 2 which is a flowchart of a method for splicing blood vessel contours provided by the embodiments of the present application. As shown in Figure 2 , the splicing method provided by the embodiments of the present application includes:

[0088] S201. Determine the branch contour information of the branch blood vessel to be spliced and the main contour information of the main blood vessel.

[0089] S202. Determine the relative spatial relationship between the branch blood vessel to be spliced and the main blood vessel according to the branch contour information and the main contour information.

[0090] S203. Determine the target contour line in the main blood vessel according to the relative spatial relationship.

[0091] S204. Determine the target correction strategy of the branch blood vessel to be spliced according to the minimum distance between the branch blood vessel to be spliced and the target contour line.

[0092] S205. Process the branch blood vessel to be spliced according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessel to be spliced coincides with the target contour line.

[0093] S206. Identify the splicing intersection points in the initial blood vessel contour map and truncate the redundant contour lines beyond the splicing intersection points to obtain a target blood vessel contour map.

[0094] The following explains each step of the embodiments of the present application by way of example:

[0095] Regarding step S201, in this step, the contour information of the main blood vessel (specifically, the contour of the main blood vessel can be known) and the branch blood vessel to be spliced (specifically, the contour of the branch blood vessel to be spliced can be referred to) can be extracted from the image data.

[0096] Here, the contour information can usually be obtained through image processing techniques and may include geometric features such as the shape, boundary, curvature, and image position of the blood vessel.

[0097] Regarding step S202, in an implementation manner provided by the present application, determining the relative spatial relationship between the branch blood vessel to be spliced and the main blood vessel according to the branch contour information and the main branch contour information includes:

[0098] S2021. According to the branch contour information and the main branch contour information, determine the interval distance between each main branch contour line in the main blood vessel and the branch blood vessel to be spliced.

[0099] S2022. According to the two interval distances, determine the relative spatial relationship between the branch blood vessel to be spliced and the main blood vessel.

[0100] Regarding step S2021, in this step, for example, the interval distance can be determined through the following steps: for each main branch contour line in the main blood vessel, according to the main branch contour information of this main branch contour line and the branch contour information of at least one branch contour line in the branch blood vessel to be spliced, determine the minimum distance between this main branch contour line and at least one branch contour line; and according to at least one minimum distance, determine the interval distance between this main branch contour line and the branch blood vessel to be spliced.

[0101] Here, when determining the interval distance between this main branch contour line and the branch blood vessel to be spliced according to at least one minimum distance, it includes: determining any one of the minimum distances as the interval distance between this main branch contour line and the branch blood vessel to be spliced; or, determining the minimum distance between this main branch contour line and a pre-specified branch contour line as the interval distance between this main branch contour line and the branch blood vessel to be spliced; or, determining the sum value or average value of the two minimum distances as the interval distance between this main branch contour line and the branch blood vessel to be spliced.

[0102] Among them, when determining the minimum distance between this main branch contour line and a pre-specified branch contour line as the interval distance between this main branch contour line and the branch blood vessel to be spliced, only the distance between this main branch contour line and the pre-specified branch contour line needs to be calculated.

[0103] Regarding step S2022, the two interval distances refer to the interval distances between two main branch contour lines and the branch blood vessel to be spliced respectively.

[0104] When determining the relative spatial relationship between the branch blood vessel to be spliced and the main blood vessel according to two spacing distances, the relative position of the two blood vessels in space can be inferred by comparing the distances between them.

[0105] For step S203, in an embodiment provided by the present application, determining the target contour line in the main blood vessel according to the relative spatial relationship includes: determining the main branch contour line on the same side as the branch blood vessel to be spliced as the target contour line.

[0106] Here, when determining the main branch contour line on the same side as the branch blood vessel to be spliced, it can be determined according to the spacing distance. For example, the main branch contour line on the side with the minimum spacing distance is determined as the main branch contour line on the same side as the branch blood vessel to be spliced.

[0107] In addition, the following content is used to illustrate steps S202 - S203 again by way of example:

[0108] First, calculate the minimum distances between one side (s l , s r ) of the branch blood vessel to be spliced S = [s l , s r refers to the coordinate points of the left and right contours. For example, take the left s l for calculation) and the two side contours of the main blood vessel M = [m l , m r , and take the side with the smaller distance as the main branch contour (target contour line) of the branch blood vessel to be spliced, denoted as M t . The mathematical formula is as follows:

[0109]

[0110] In this way, the minimum distance between each side of the branch contour to be spliced and the target contour line can be determined, and the formula is as follows:

[0111] d = min(s l - M t )

[0112] For step S204, here, the target correction strategies corresponding to different minimum distances may also be different.

[0113] In an embodiment provided by the present application, determining the target correction strategy of the branch blood vessel to be spliced according to the minimum distance between the branch blood vessel to be spliced and the target contour line includes:

[0114] S2041. Determine whether the minimum distance between the branch blood vessel to be spliced and the target contour line is less than a preset distance.

[0115] S2042. The target correction strategy is to perform an outward extension process on the branch blood vessel to be spliced according to the minimum distance.

[0116] S2043. The target correction strategy is to update the branch blood vessel to be spliced according to the original medical image, and perform an outward extension process on the updated branch blood vessel to be spliced.

[0117] For step S2041, compare the minimum distance between the branch blood vessel to be spliced and the target contour line with a preset distance. If it is less than the preset distance, execute step S2042; if it is not less than the preset distance, execute step S2043.

[0118] Here, the preset distance can be set according to experience. For example, set 20 pixel distances as the preset distance. If it is not less than the preset distance, it means that the distance between the branch blood vessel to be spliced and the target contour line is relatively far, and it is not suitable for direct splicing.

[0119] For step S2043, when the minimum distance between the branch blood vessel to be spliced and the target contour line is not less than the preset distance, the target correction strategy can also be to update the branch blood vessel to be spliced according to the original medical image, and perform an outward extension process on the branch blood vessel to be spliced that has not been successfully spliced after the update.

[0120] Here, if the splicing is successful, an initial blood vessel contour map can be generated, and then step S206 can be executed.

[0121] For step S205, in this step, perform an extension or update process on the branch blood vessel to be spliced according to the target correction strategy, so that the processed branch blood vessel to be spliced is successfully spliced with the target contour line or spliced beyond it, and an initial blood vessel contour map is obtained.

[0122] In an embodiment provided by the present application, when the minimum distance between the branch blood vessel to be spliced and the target contour line is less than the preset distance, performing processing on the branch blood vessel to be spliced according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessel to be spliced coincides with the target contour line includes:

[0123] S20511. For each branch contour line of the branch blood vessel to be spliced, determine the splicing starting point of this branch contour line.

[0124] S20512. Determine the extension distance according to the shortest distance between the splicing starting point of this branch contour line and the target contour line.

[0125] S20513. According to the branch contour information of this branch contour line, starting from the splicing starting point of this branch contour line, use the interpolation method to extend this branch contour line by the extension distance.

[0126] S20514. When all the branch contour lines are completely extended, an initial blood vessel contour map is obtained.

[0127] Regarding step S20511, this step may include: for each branch contour line in the branch blood vessel to be spliced, determining the two end points of this branch contour line; and determining the end point closest to the target contour line as the splicing starting point of this branch contour line.

[0128] Regarding step S20512, according to the splicing starting point of this branch contour line and the main branch contour information of the target contour line, determining the shortest distance between the splicing starting point of this branch contour line and the target contour line.

[0129] Before determining the shortest distance between the splicing starting point of this branch contour line and the target contour line, the method further includes: performing interpolation processing on the main branch contour information of the target contour line. Interpolation is to estimate values between known data points.

[0130] The extension distance is a preset multiple of the shortest distance.

[0131] Exemplarily, the extension distance is 2 times the shortest distance.

[0132] Regarding step S20513, the interpolation method may adopt cubic spline interpolation. The interpolation method performed in this step is an extrapolation processing method. Extrapolation is to estimate values outside the range of known data points.

[0133] It should be noted that cubic spline is a piecewise-defined polynomial interpolation method. It fits a cubic polynomial between each adjacent pair of points and ensures that these polynomials have continuous first and second derivatives at the nodes, thereby ensuring the smoothness of the blood vessel.

[0134] Regarding step S20514, on the basis of the original blood vessel contour map (generally, there is a disconnection between branches and the main branch), after performing outward extension processing on each branch contour line in the branch blood vessel to be spliced, an initial blood vessel contour map is obtained.

[0135] In another implementation manner provided by the present application, when the minimum distance between the branch blood vessel to be spliced and the target contour line is not less than a preset distance, processing the branch blood vessel to be spliced according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessel to be spliced coincides with the target contour line includes:

[0136] S20521. Obtain the original medical image including the branch blood vessel to be spliced and the main branch blood vessel.

[0137] S20522. Perform skeletonization processing on the original medical image to determine the full blood vessel center line.

[0138] S20523. Extract the centerline corresponding to the branch blood vessel to be spliced from the full blood vessel centerline to obtain the branch centerline to be processed.

[0139] S20524. According to the full blood vessel centerline, use the depth-first traversal algorithm to complete the branch centerline to be processed, and determine the complete branch centerline corresponding to the branch blood vessel to be spliced and intersecting with the centerline of the main branch blood vessel.

[0140] S20525. According to the complete branch centerline and the original medical image, update the branch contour information of the branch blood vessel to be spliced to obtain the updated branch blood vessel to be spliced.

[0141] S20526. According to the shortest distance between the updated branch blood vessel to be spliced and the target contour line, use the interpolation method to perform an outer extension process on the updated branch blood vessel to be spliced to obtain an initial blood vessel contour map.

[0142] For step S20521, the original medical image can be a CT image, an angiogram image, etc.

[0143] For step S20522, performing a skeletonization process on the original medical image is to extract the blood vessel centerline in the original medical image. By way of example, please refer to Figure 3 , Figure 3 which is a schematic diagram of the blood vessel centerline provided by this application.

[0144] For step S20523, in this step, the centerline of the branch blood vessel to be spliced (branch centerline to be processed) is separated from the full blood vessel centerline.

[0145] Among them, the extracted branch centerline may be broken or incomplete, so subsequent completion is required.

[0146] By way of example, please continue to refer to Figure 3 , Figure 3 where the thick red line is the branch centerline to be processed.

[0147] For step S20524, by way of example, this step may include: selecting the point with the shortest blood vessel distance from the centerline coordinates to the main branch centerline as the target coordinate point; then, based on the full blood vessel centerline, using the depth-first traversal algorithm to search for the point coordinates in the 8-neighborhood. If the pixel value of the coordinate is 1, it is added to the path list until the traversal reaches the bifurcation point of the main branch (such as Figure 3 the diamond position); merging the path list and the original branch centerline list into a new centerline list, that is, determining the complete branch centerline.

[0148] For step S20525, based on the complete branch centerline and the original medical image, contour segmentation processing is performed again to obtain the newly generated branch blood vessels to be spliced; the originally to-be-spliced branch blood vessels are replaced with the newly generated branch blood vessels to be spliced to complete the update process.

[0149] For step S20526, the description of this step can refer to the description of S20511 to S20514 and can achieve the same technical effect, so it will not be elaborated here.

[0150] For step S206, in an implementation manner provided by the present application, identifying the splicing intersection points in the initial blood vessel contour diagram and truncating the redundant contour lines beyond the splicing intersection points to obtain the target blood vessel contour diagram includes:

[0151] S2061. For each branch contour line in the to-be-spliced branch blood vessels after processing in the initial blood vessel contour diagram, determine the extended contour point set of this branch contour line.

[0152] S2062. According to the extended contour point set of this branch contour line and the main fulcrum set of the target contour line, determine the shortest distance between this branch contour line and the target contour line.

[0153] S2063. Determine the splicing intersection point between this branch contour line and the target contour line as the extended contour point in the shortest distance.

[0154] S2064. Truncate the contour line between the splicing intersection point and the extended end point on this branch contour line to obtain the target blood vessel contour diagram.

[0155] It should be noted that whether it is a branch contour line, a main branch contour line, or even a blood vessel centerline, the corresponding data form is generally a discrete point set.

[0156] For step S2061, the extended contour point set is the point set determined after processing this branch contour line according to step S205. The extended end point in the extended contour point set is the other end point of this branch contour line after processing.

[0157] For step S2062, this step includes: randomly selecting a point from the extended contour point set and the main fulcrum set respectively, calculating the interval distance between this branch contour line and the target contour line; and determining the shortest distance between this branch contour line and the target contour line according to multiple interval distances.

[0158] Among them, when the number of contour points of the target contour line is small, the main fulcrum set can be a dense point set generated after interpolation processing.

[0159] For step S2063, this step includes: for each branch contour line, determining two points corresponding to the shortest distance between this branch contour line and the target contour line, and determining the point (expanded contour point) belonging to the expanded contour point set among them as the splicing intersection point.

[0160] For step S2064, this step includes: for each branch contour line, truncating the contour line from the splicing intersection point to the expanded end point on this branch contour line to obtain the target blood vessel contour map.

[0161] In addition, in another implementation manner provided by the present application, the redundant contour lines in the initial blood vessel contour map can also be truncated and deleted in the following manner to obtain the target blood vessel contour map. The specific implementation manner is as follows: Please refer to Figure 4 , Figure 4 which is an example diagram of the determination process of a target blood vessel contour map provided by the present application. As Figure 4 shown, with the starting point of the splicing point of each extrapolated branch as the center, search within a rectangular range of a set size (such as the orange frame in Figure 4 the extrapolated diagram (initial blood vessel contour map)), search for the point closest to the main branch dense point among the extrapolated contour points of the branch, and truncate at this point. Perform the same processing for each contour, and then the Figure 4 splicing diagram (target blood vessel contour map) can be obtained.

[0162] In this way, when generating the blood vessel contour in the present application, the connectivity between the branch blood vessels and the main branch blood vessels is considered, and the unconnected branch blood vessels and main branch blood vessels are spliced by combining the segmentation method. At the same time, during the splicing process, different splicing methods are adopted for the branch blood vessels with different distances from the main branch blood vessels, so as to take into account the splicing accuracy of the branch blood vessels far from the main branch, and further ensure the accuracy of the splicing result, the accuracy of the simulation result, and make the present solution have a wider applicability.

[0163] Based on the same inventive concept, a splicing device corresponding to the splicing method is also provided in the embodiments of the present application. Since the principle of solving problems by the device in the embodiments of the present application is similar to the above splicing method in the embodiments of the present application, the implementation of the device can refer to the implementation of the method, and the repeated parts will not be described again.

[0164] Please refer to Figure 5 , Figure 5 which is a schematic structural diagram of a blood vessel contour splicing device provided by the embodiments of the present application. As Figure 5 shown in, the splicing device 500 includes:

[0165] A first determination module 510, configured to determine the branch contour information of the branch blood vessel to be spliced and the main branch contour information of the main branch blood vessel;

[0166] The second determination module 520 is configured to determine the relative spatial relationship between the branch blood vessel to be spliced and the main blood vessel according to the branch contour information and the main branch contour information;

[0167] The third determination module 530 is configured to determine a target contour line in the main blood vessel according to the relative spatial relationship;

[0168] The fourth determination module 540 is configured to determine a target correction strategy for the branch blood vessel to be spliced according to the minimum distance between the branch blood vessel to be spliced and the target contour line;

[0169] The processing module 550 is configured to process the branch blood vessel to be spliced according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessel to be spliced coincides with the target contour line;

[0170] The truncation module 560 is configured to identify the splicing intersection point in the initial blood vessel contour map and truncate the redundant contour lines beyond the splicing intersection point to obtain a target blood vessel contour map.

[0171] Optionally, when the fourth determination module 540 is configured to determine a target correction strategy for the branch blood vessel to be spliced according to the minimum distance between the branch blood vessel to be spliced and the target contour line, the fourth determination module 540 is configured to:

[0172] Determine whether the minimum distance between the branch blood vessel to be spliced and the target contour line is less than a preset distance;

[0173] If it is less, the target correction strategy is to perform an outer extension process on the branch blood vessel to be spliced according to the minimum distance;

[0174] If it is not less, the target correction strategy is to update the branch blood vessel to be spliced according to the original medical image and perform an outer extension process on the updated branch blood vessel to be spliced.

[0175] Optionally, when the minimum distance between the branch blood vessel to be spliced and the target contour line is less than the preset distance, when the processing module 550 is configured to process the branch blood vessel to be spliced according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessel to be spliced coincides with the target contour line, the processing module 550 is configured to:

[0176] For each branch contour line of the branch blood vessel to be spliced, determine the splicing starting point of the branch contour line;

[0177] Determine the extension distance according to the shortest distance between the splicing starting point of the branch contour line and the target contour line;

[0178] According to the branch contour information of the branch contour line, starting from the splicing starting point of the branch contour line, the branch contour line is extended by the extension distance using an interpolation method;

[0179] When all branch contour lines are fully extended, an initial blood vessel contour map is obtained.

[0180] Optionally, when the minimum distance between the branch blood vessel to be spliced and the target contour line is not less than a preset distance, when the processing module 550 is used to process the branch blood vessel to be spliced according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessel to be spliced coincides with the target contour line, the processing module 550 is used to:

[0181] Obtain the original medical image including the branch blood vessel to be spliced and the main blood vessel;

[0182] Perform skeletonization processing on the original medical image to determine the full blood vessel centerline;

[0183] Extract the centerline corresponding to the branch blood vessel to be spliced from the full blood vessel centerline to obtain the branch centerline to be processed;

[0184] According to the full blood vessel centerline, use the depth-first traversal algorithm to perform complement processing on the branch centerline to be processed to determine the complete branch centerline corresponding to the branch blood vessel to be spliced that intersects the centerline of the main blood vessel;

[0185] According to the complete branch centerline and the original medical image, update the branch contour information of the branch blood vessel to be spliced to obtain the updated branch blood vessel to be spliced;

[0186] According to the shortest distance between the updated branch blood vessel to be spliced and the target contour line, perform an outer extension process on the updated branch blood vessel to be spliced using an interpolation method to obtain an initial blood vessel contour map.

[0187] Optionally, when the truncation module 560 is used to identify the splicing intersection points in the initial blood vessel contour map and truncate the redundant contour lines beyond the splicing intersection points to obtain the target blood vessel contour map, the truncation module 560 is used to:

[0188] For each branch contour line in the processed branch blood vessel to be spliced in the initial blood vessel contour map, determine the extended contour point set of the branch contour line; wherein, the extended end point in the extended contour point set is the other end point of the branch contour line;

[0189] According to the extended contour point set of the branch contour line and the main support point set of the target contour line, determine the shortest distance between the branch contour line and the target contour line;

[0190] Determine the extended contour points in the shortest distance as the splicing intersection points between this branch contour line and the target contour line;

[0191] Truncate the contour line between the splicing intersection point and the extended end point on this branch contour line to obtain the target blood vessel contour map.

[0192] Optionally, when the second determination module 520 is used to determine the relative spatial relationship between the branch blood vessel to be spliced and the main blood vessel according to the branch contour information and the main branch contour information, the second determination module 520 is used for:

[0193] Determine the interval distance between each main branch contour line in the main blood vessel and the branch blood vessel to be spliced according to the branch contour information and the main branch contour information;

[0194] Determine the relative spatial relationship between the branch blood vessel to be spliced and the main blood vessel according to the two interval distances.

[0195] Optionally, when the third determination module 530 is used to determine the target contour line in the main blood vessel according to the relative spatial relationship, the third determination module 530 is used for:

[0196] Determine the main branch contour line that belongs to the same side as the branch blood vessel to be spliced as the target contour line.

[0197] Please refer to Figure 6 , Figure 6 which is a schematic structural diagram of an electronic device provided by an embodiment of the present application. As Figure 6 shown in, the electronic device 600 includes a processor 610, a memory 620, and a bus 630.

[0198] The memory 620 stores machine-readable instructions executable by the processor 610. When the electronic device 600 runs, the processor 610 communicates with the memory 620 through the bus 630. When the machine-readable instructions are executed by the processor 610, the steps in the method embodiment as described above can be executed. The specific implementation manner can refer to the method embodiment and will not be elaborated here. Figures 1 to 4 shown, and the specific implementation manner can refer to the method embodiment and will not be elaborated here.

[0199] An embodiment of the present application further provides a computer-readable storage medium. A computer program is stored on the computer-readable storage medium. When the computer program is run by a processor, the steps in the method embodiment as described above can be executed. The specific implementation manner can refer to the method embodiment and will not be elaborated here. Figures 1 to 4 shown, and the specific implementation manner can refer to the method embodiment and will not be elaborated here.

[0200] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the systems, devices, and units described above can refer to the corresponding processes in the foregoing method embodiments and will not be elaborated herein.

[0201] In several embodiments provided in the present application, it should be understood that the disclosed systems, devices, and methods can be implemented in other ways. The device embodiments described above are merely illustrative. For example, the division of the units is only a logical function division, and there may be other division methods in actual implementation. For another example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed couplings or direct couplings or communication connections to each other can be through some communication interfaces. The indirect couplings or communication connections of the devices or units can be in electrical, mechanical, or other forms.

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

[0203] In addition, in each embodiment of the present application, the functional units can be integrated into one processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit.

[0204] If the function is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a non-volatile computer-readable storage medium executable by a processor. Based on this understanding, the technical solution of the present application, in essence, or the part that contributes to the prior art, or a part of this 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 for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in each embodiment of the present application. The foregoing storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical discs that can store program codes.

[0205] Finally, it should be noted that the above-described embodiments are only specific implementation manners of the present application, used to illustrate the technical solutions of the present application, rather than limiting it. The protection scope of the present application is not limited thereto. Although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that any person skilled in the art within the technical scope disclosed by the present application can still modify the technical solutions recorded in the foregoing embodiments, or can easily think of changes, or perform equivalent replacements on some of the technical features; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application, and should all be covered within the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims

1. A method for splicing blood vessel contours, characterized in that: The splicing method comprises: Determine the branch contour information of the branch blood vessels to be spliced ​​and the main branch contour information of the main blood vessels; Determining the relative spatial relationship between the branch blood vessel to be spliced ​​and the main blood vessel according to the branch outline information and the main blood vessel outline information; determining a target contour line in the main branch vessel according to the relative spatial relationship; Determining a target correction strategy for the branch blood vessels to be spliced ​​according to a minimum distance between the branch blood vessels to be spliced ​​and the target contour line; Processing the branch blood vessels to be spliced ​​according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessels to be spliced ​​coincide with the target contour line; The splicing intersection points in the initial blood vessel contour map are identified, and redundant contour lines beyond the splicing intersection points are cut off to obtain a target blood vessel contour map.

2. The splicing method according to claim 1, characterized in that: Determining a target correction strategy for the branch blood vessels to be spliced ​​according to the minimum distance between the branch blood vessels to be spliced ​​and the target contour line includes: Determining whether the minimum distance between the branch blood vessel to be spliced ​​and the target contour line is less than a preset distance; If it is less than, the target correction strategy is to perform an external extension process on the branch blood vessel to be spliced ​​according to the minimum distance; If it is not less than, the target correction strategy is to update the branch blood vessels to be spliced ​​according to the original medical image, and perform external extension processing on the updated branch blood vessels to be spliced.

3. The splicing method according to claim 2, characterized in that: When the minimum distance between the branch blood vessel to be spliced ​​and the target contour line is less than a preset distance, the branch blood vessel to be spliced ​​is processed according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessel to be spliced ​​coincides with the target contour line, including: For each branch contour line of the branch blood vessel to be spliced, determining a splicing starting point of the branch contour line; Determine the extension distance according to the shortest distance between the splicing starting point of the branch contour line and the target contour line; According to the branch contour information of the branch contour line, starting from the splicing starting point of the branch contour line, an interpolation method is used to extend the branch contour line by the extension distance; When all branch contour lines are extended, the initial blood vessel contour map is obtained.

4. The splicing method according to claim 2, characterized in that: When the minimum distance between the branch blood vessels to be spliced ​​and the target contour line is not less than a preset distance, the branch blood vessels to be spliced ​​are processed according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessels to be spliced ​​coincide with the target contour line, including: Acquiring an original medical image including the branch blood vessel to be spliced ​​and the main blood vessel; Performing skeleton processing on the original medical image to determine the centerline of the entire blood vessel; Extracting the centerline corresponding to the branch blood vessel to be spliced ​​from the centerline of the whole blood vessel to obtain the centerline of the branch blood vessel to be processed; According to the centerline of the whole blood vessel, a depth-first traversal algorithm is used to complete the centerline of the branch to be processed, and a complete branch centerline corresponding to the branch blood vessel to be spliced ​​and intersecting with the centerline of the main branch blood vessel is determined; updating the branch contour information of the branch blood vessel to be spliced ​​according to the complete branch centerline and the original medical image to obtain an updated branch blood vessel to be spliced; According to the shortest distance between the updated branch blood vessels to be spliced ​​and the target contour line, the updated branch blood vessels to be spliced ​​are extended outward by using an interpolation method to obtain an initial blood vessel contour map.

5. The splicing method according to claim 1, characterized in that: The step of identifying the splicing intersection points in the initial blood vessel contour map and truncating the redundant contour lines beyond the splicing intersection points to obtain the target blood vessel contour map includes: For each branch contour line in the branch blood vessel to be spliced ​​after processing in the initial blood vessel contour map, determine an extended contour point set of the branch contour line; wherein the extended end point in the extended contour point set is the other end point of the branch contour line; Determine the shortest distance between the branch contour line and the target contour line according to the extended contour point set of the branch contour line and the main support point set of the target contour line; Determine the extended contour point in the shortest distance as the splicing intersection point between the branch contour line and the target contour line; The contour line between the splicing intersection point and the extension end point on the branch contour line is truncated to obtain the target blood vessel contour map.

6. The splicing method according to claim 1, characterized in that: The determining, based on the branch outline information and the main branch outline information, the relative spatial relationship between the branch blood vessel to be spliced ​​and the main branch blood vessel comprises: Determine the interval distance between each main branch contour line in the main branch blood vessel and the branch blood vessel to be spliced ​​according to the branch contour information and the main branch contour information; The relative spatial relationship between the branch blood vessel to be spliced ​​and the main blood vessel is determined according to the two interval distances.

7. The splicing method according to claim 6, characterized in that: Determining the target contour line in the main branch vessel according to the relative spatial relationship includes: The main branch contour line that belongs to the same side as the branch blood vessel to be spliced ​​is determined as the target contour line.

8. A device for splicing blood vessel contours, characterized in that: The splicing device comprises: A first determination module is used to determine the branch contour information of the branch blood vessels to be spliced ​​and the main branch contour information of the main blood vessels; A second determination module is used to determine the relative spatial relationship between the branch blood vessel to be spliced ​​and the main branch blood vessel according to the branch outline information and the main branch outline information; A third determination module is used to determine a target contour line in the main branch vessel according to the relative spatial relationship; A fourth determination module, configured to determine a target correction strategy for the branch blood vessels to be spliced ​​according to a minimum distance between the branch blood vessels to be spliced ​​and the target contour line; A processing module, used for processing the branch blood vessels to be spliced ​​according to the target correction strategy to obtain an initial blood vessel contour map in which the processed branch blood vessels to be spliced ​​coincide with the target contour line; The truncation module is used to identify the splicing intersection points in the initial blood vessel contour map, and to truncate the redundant contour lines beyond the splicing intersection points to obtain the target blood vessel contour map.

9. An electronic device, characterized in that: include: A processor, a memory and a bus, wherein the memory stores machine-readable instructions executable by the processor, and when the electronic device is running, the processor and the memory communicate through the bus, and the machine-readable instructions are executed by the processor to execute the steps of the splicing method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the splicing method according to any one of claims 1 to 7 are executed.

Images