A coronary artery blood supply evaluation system based on parabolic three-dimensional coordinate system

Abstract

The present invention proposes a coronary artery blood supply evaluation system based on a parabola three-dimensional coordinate system, which includes locating the intersection point A of the parabola corresponding to the blood supply of the left anterior descending coronary artery determined by the coronary angiography results and the parabola corresponding to the blood supply of the right coronary artery in the X-axis plane in the parabola three-dimensional coordinate system. 1， The vertical axis marks the blood flow of the left coronary artery circumflex branch B, and the area S enclosed by the parabola F1(x) passing through the polar coordinates A1 and B is drawn RAC That is, the "coronary artery abundance area". Area S RAC The slope of the tangent line of polar coordinate A1 can intuitively reflect the coronary blood flow characteristics of the sample heart. By simulating different coronary revascularization schemes before implementing revascularization, the effects of different revascularization schemes can be intuitively demonstrated, allowing the optimal revascularization strategy to be selected. Furthermore, this parabolic three-dimensional coordinate system can simultaneously represent multiple sample data, making it suitable for studying coronary blood flow characteristics in large sample sizes.

Description

Technical Field

[0001] The present invention belongs to the field of statistics, and in particular relates to a coronary artery blood supply evaluation system based on a parabola three-dimensional coordinate system. Background Art

[0002] Since F. Mason Sones performed the first coronary angiography in 1960, the study of coronary artery blood flow has been limited to two-dimensional images such as tables and charts. The combination of intravascular ultrasound (IVUS) and optical coherence tomography (OCT) has significantly improved the quantification and characterization of coronary artery stenosis. However, due to the inherent technical limitations of IVUS, only one coronary artery can be observed at a time. Furthermore, the high cost of the equipment limits its clinical application.

[0003] The quantitative flow fraction (QFR) successfully developed in 2018 achieved the leap of digital angiography (DSA) from grayscale (Greyscale) to color (RGB color system), but its role in the global evaluation of the heart's coronary blood supply is limited. The analysis weights of each segment of the coronary artery cannot be adjusted, which limits the analysis of a certain segment of the coronary artery. Summary of the Invention

[0004] The purpose of the present invention is to provide intuitive visualization of coronary artery blood supply conditions, and at the same time, to provide a coronary artery blood supply evaluation system that is convenient for evaluating changes in blood supply to a certain segment or certain segments.

[0005] To achieve the above-mentioned object, the present invention adopts the following technical solution: a coronary artery blood supply evaluation system based on a parabolic three-dimensional coordinate system, comprising the parabolic three-dimensional coordinate system, and also comprising, in the parabolic three-dimensional coordinate system, an intersection point A1 of a parabola corresponding to the blood supply status of the left anterior descending coronary artery and a parabola corresponding to the blood supply status of the right coronary artery on the X-axis plane, and a vertical coordinate marking the blood supply status of the left circumflex coronary artery, B;

[0006] Also includes drawing the area S enclosed by the parabola F1(x) passing through the polar coordinates A1 and B RAC For the coronary artery abundance area, a polar coordinate A with A1 and A2 coordinate systems is established. The tangent slope of the polar coordinate A can intuitively reflect the coronary blood supply characteristics of the heart.

[0007] The intersection point A of the parabola corresponding to the left anterior descending coronary artery and the right coronary artery 1， and a virtual symmetrical position A2 which is symmetrical to the origin of A1;

[0008] The characteristics of coronary blood supply of the heart are expressed by the polar coordinates of the A1 and A2 coordinate systems and the slope of the tangent line between A1 and A2;

[0009] With the blood supply of the left coronary artery circumflex branch marked on the vertical axis B, draw the parabola F1 (x) passing through the polar coordinates A1 and B. RAC , S RAC Extracting the characteristics of the coronary blood supply of the heart; the slope of the parabola F1 (x) pointing from the polar coordinate A1 on the X-axis plane to the ordinate can extract the characteristics of the coronary blood supply of the heart;

[0010] By simulating the changes in the stenosis of each coronary artery, we can obtain the "coronary artery abundance area" corresponding to different F1 (x), that is, S RAC , thereby estimating the reconstruction results and prognosis before implementing revascularization, selecting the optimal revascularization strategy, and providing clinical revascularization decisions.

[0011] The system includes the following steps:

[0012] S1, express the coronary angiography results using the common origin equation:

[0013] The equation for the stenosis rate of the left anterior descending coronary artery is: X=-(N+1)Y 2 +(1-0.1N);

[0014] The equation for the right coronary artery stenosis rate is: X=(N+1)Y 2 -(1-0.1N);

[0015] N is determined based on the absolute value of the filling defect at two phases with an angle of 90 degrees on DSA angiography of the same heart's coronary arteries;

[0016] S2, the origin O is set as the Z axis, which is used to represent the stenosis rate B of the left circumflex coronary artery;

[0017] S3, complete the establishment of the parabolic three-dimensional coordinate system through S1 and S2;

[0018] S4, determine the intersection A1 of the left anterior descending coronary artery and the right coronary artery in the X-axis plane;

[0019] S5, virtual symmetrical position A2 of virtual A1 relative to the origin, A2 is used to assist in completing the establishment of F1 (X);

[0020] S6, set A1, A2 and B as the polar coordinates of the parabola F1(X), and B as the vertex coordinate of F1(X);

[0021] S7, the area enclosed by A1 and B on the parabola F1 (X) is defined as the coronary artery abundance area, that is, S RAC ;

[0022] S8, combined with the parabolic three-dimensional coordinate system to simulate the changes in different coronary artery segments.

[0023] In a preferred embodiment of the present invention, the coordinate system is established in step 1 as follows: X=(N+1)Y 2 -(1-0.1N); where N is a natural number greater than 0;

[0024] Then y can be expressed as ;

[0025] Draw the parabola plane, and then implement the coordinate Z perpendicular to the parabola plane at the origin of the parabola plane to complete the setting of the parabola three-dimensional coordinate system.

[0026] In a preferred embodiment of the present invention, the coordinate definition is also included:

[0027] Positive x-axis displacement (x) represents stenosis of the left anterior descending coronary artery;

[0028] Negative x-axis displacement (-x) represents stenosis of the right coronary artery;

[0029] Positive z-axis displacement (z) represents stenosis of the left circumflex coronary artery.

[0030] In a preferred embodiment of the present invention, in a parabolic coordinate system, any X and any -X have two intersection points located symmetrically on both sides of the X axis, wherein the left coronary dominant intersection point is within the first and fourth quadrants; and the right coronary dominant intersection point is within the second and third quadrants.

[0031] In another preferred embodiment of the present invention, in step S7, the area S of the image domain enclosed by A1 and B is RAC Calculation methods include ;

[0032] Where x1 represents the X coordinate of the intersection point A1 of the left anterior descending coronary artery and the right coronary artery in the X-axis plane; x2 represents the X coordinate of the virtual point A2 symmetrical with respect to the origin of A1, satisfying x2=-x1.

[0033] The above scheme achieved the following beneficial effects: 1. It provided an intuitive three-dimensional visualization of the coronary blood supply of the heart, showing the coronary blood supply of the research sample in an intuitive way, no longer simply in the form of percentages or rates;

[0034] 2. Through the polar coordinates and area S of the sample in the coordinate system RAC , polar coordinate slope, A1 in S RAC The tangent slope can be used to extract the coronary blood supply characteristics of a large sample of heart, and then deduce the clinical significance behind it, providing a new research method for the study of coronary blood supply of a large sample of heart.

[0035] 3. The parabolic coordinate system platform provides a global perspective for evaluating the coronary blood supply of the heart. Rather than studying one or two coronary arteries individually, all coronary branches are studied as a whole. Different weights are assigned to each segment based on the research purpose, fully exposing positive results and minimizing interference.

[0036] 4. Before implementing revascularization, by simulating the changes in different coronary artery segments, the results of different revascularization strategies can be displayed. Combined with previous clinical treatment effects, the most optimized revascularization strategy can be selected to provide clinical revascularization decision-making.

[0037] Additional aspects and advantages of the present invention will be set forth in part in the description which follows and, in part, will be obvious from the description which follows, or may be learned by practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] The above and / or additional aspects and advantages of the present invention will become apparent and readily understood from the following description of the embodiments with reference to the accompanying drawings, in which:

[0039] Figure 1 This is a diagram for establishing a parabolic three-dimensional coordinate system according to an embodiment of the present application.

[0040] Figure 2 yes Figure 1 Detailed description of the embodiment

[0041] The following describes embodiments of the present invention in detail. Examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended only to explain the present invention and are not to be construed as limiting the present invention.

[0042] In the description of the present invention, it should be understood that the terms "longitudinal", "transverse", "vertical", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside", etc., indicating the orientation or position relationship, are based on the orientation or position relationship shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention.

[0043] In the description of the present invention, unless otherwise specified and limited, it should be noted that the terms "installed", "connected" and "connected" should be understood in a broad sense. For example, it can be a mechanical connection or an electrical connection, or it can be the internal communication between two components. It can be a direct connection or an indirect connection through an intermediate medium. For ordinary technicians in this field, the specific meanings of the above terms can be understood according to the specific circumstances.

[0044] The present invention provides a coronary artery blood supply evaluation system based on a parabolic three-dimensional coordinate system, such as Figure 1 and Figure 2 As shown,

[0045] The key scientific issues to be addressed in this technical solution are:

[0046] 1. It provides an intuitive three-dimensional visualization of the coronary artery blood supply, offering a new research method for large-sample coronary artery blood supply studies;

[0047] 2. Express the coronary blood supply characteristics of the research sample through the coordinate system area and tangent slope;

[0048] 3. Provide research methods for coronary angiography results of specific segments;

[0049] 4. Quantify and verify the potential benefits of the upcoming revascularization strategy to assist in selecting the optimal revascularization strategy.

[0050] In view of the above problems that need to be solved, the parabolic three-dimensional coordinate system is first established.

[0051] The formula for drawing a parabola solid coordinate system is:

[0052] X=(N+1)Y 2 -(1-0.1N) (1)

[0053] Then y can be expressed as (2)

[0054] X=-(N+1)Y 2 +(1-0.1N) (3)

[0055] y can be expressed as (4)

[0056] Where N is a natural number greater than 0

[0057] The determination of N is based on the absolute value of the filling defects at two phases with an angle of 90 degrees in DSA angiography of the same heart's coronary arteries. Therefore, the semi-diameter of the parabola parallel to the Y axis corresponding to each DSA angiography result is different.

[0058] Draw the parabola plane, then implement the coordinate Z at the parabola origin perpendicular to the parabola plane to complete the setting of the parabola coordinate system. The position of a point in the parabola coordinate system is expressed as (x, y, z).

[0059] *The Y-axis is introduced here only to facilitate the calculation of the formula and does not participate in the final result expression

[0060] Coordinate definition:

[0061] Positive x-axis displacement (x) represents the stenosis of the left anterior descending coronary artery

[0062] Negative x-axis displacement (-x) represents stenosis of the right coronary artery

[0063] Positive z-axis displacement (z) represents the stenosis of the left circumflex coronary artery

[0064] In this coordinate system, any X and any -X have two intersection points located symmetrically on both sides of the X axis. The intersection points of the left-coronary dominant type are in the first and fourth quadrants, while the intersection points of the right-coronary dominant type are in the second and third quadrants.

[0065] Definition of Quadrant

[0066] The area enclosed by the positive x-axis and the positive y-axis is the first quadrant

[0067] The area enclosed by the negative x-axis and the positive y-axis is the second quadrant

[0068] The area enclosed by the negative directions of the x and y axes is the third quadrant

[0069] The area enclosed by the negative y-axis and the positive x-axis is the fourth quadrant

[0070] Definition of the image domain

[0071] The space enclosed by the positive x, y, and z axes is the first image domain

[0072] The space enclosed by the positive y and z axes and the negative x axis is the second image domain.

[0073] The space enclosed by the negative x, y axes and the positive z axis is the third image domain.

[0074] The space enclosed by the positive x, z and negative y directions is the fourth image domain.

[0075] Use of coordinate system:

[0076] According to the DSA results, the parabola of the left anterior descending coronary artery and the right coronary artery is determined to be the intersection point A1 (X1, 0) of the X-axis plane, and its virtual symmetrical position relative to the origin is A2 (X2, 0), and the parabola formula Y = aX 2 +n (5)

[0077] Then the distance from the coaxial origin O to the intersection point is recorded as OX1.

[0078] but (6)

[0079] The Laplace operator for a parabolic solid coordinate system is:

[0080] (7)

[0081] Then the slope calculation formula of (X1, 0) is:

[0082] Y ‘ (X1) = 2aX1 (7)

[0083] Then the area S of the image domain in this coordinate system is RAC The formula is

[0084] (8)

[0085] By simulating the changes in different coronary artery segments before implementing revascularization, the results of different revascularization strategies can be displayed. Combined with previous clinical treatment effects, the most optimized revascularization strategy can be selected to provide clinical revascularization decisions.

[0086] The three-dimensional visualization of coronary blood supply characteristics can provide an intuitive understanding of the coronary blood supply situation and provide a new method for the study of clinical revascularization strategies and coronary angiography results.

[0087] 1. It provides an intuitive three-dimensional visualization of the coronary blood supply of the heart, showing the coronary blood supply of the research sample in an intuitive way, no longer simply in the form of percentage or rate;

[0088] 2. The polar coordinates, area, polar coordinate slope, and tangent slope of the sample in the coordinate system can be used to extract the coronary blood supply characteristics of a large sample, and then deduce the meaning behind them, providing a new research method for large-sample coronary blood supply research;

[0089] 3. The parabolic coordinate system platform provides a global perspective for evaluating the coronary blood supply of the heart. Rather than studying one or two coronary arteries individually, all coronary branches are studied as a whole. Different weights are assigned to each segment based on the research purpose, fully exposing positive results and minimizing interference.

[0090] 4. Before implementing revascularization, by simulating the changes in different coronary artery segments, the results of different revascularization strategies can be displayed. Combined with previous clinical treatment effects, the most optimized revascularization strategy can be selected to provide clinical revascularization decision-making.

[0091] Application Example 1:

[0092] The coronary DSA results show that the left anterior descending coronary artery is stenotic by 70%, the right coronary artery is stenotic by 30%, and the left circumflex coronary artery is stenotic by 60%. Substituting these into formulas (1) and (3), we obtain:

[0093] X L =-(3+1)Y 2 +(1-0.1x3)

[0094] X R =(7+1)Y 2 -(1-0.1x7)

[0095] That is:

[0096] -8Y 2 +0.3=4Y 2 -0.7

[0097] y 0.2886

[0098] It can be seen that X 0.37;

[0099] If it is a left coronary dominant type, the intersection is located in the first quadrant and the first quadrant.

[0100] The distance from the coaxial origin O to the intersection point is recorded as OX, and then substituted into formula (6) to obtain:

[0101] 0.46;

[0102] LCx stenosis 60%, i.e. z = 0.6

[0103] Then, according to the parabola formula y=ax 2 +n, and substituting the coordinates of the two vertices (0, 0.6) and (0.46, 0) into the parabola formula, we get:

[0104] y=-2.8x 2 +0.6

[0105] Then (0.46, 0) is at y=-2.8x 2 The slope of +0.6 (7) is -2.58

[0106] The Ratio of Abundant Coronary flow (RAC) is OX and z is the area of ​​the fourth image domain in this coordinate system ( Figure 1 (medium blue transparent area), and substituting it into formula (8) yields:

[0107] S RAC=0.235

[0108] Then the area of ​​the right triangle with two bases 'S is:

[0109] 'S RAC = 0.138

[0110] Application Example 2:

[0111] The coronary DSA results show that the left anterior descending coronary artery is 40% narrowed, the right coronary artery is 90% narrowed, and the left circumflex coronary artery is 70% narrowed. Substituting into formulas (1) and (3) we get:

[0112] X L =-(6+1)Y 2 +(1-0.1x6)

[0113] X R =(1+1)Y 2 -(1-0.1x1)

[0114] That is:

[0115] -7Y 2 +0.4=2Y 2 -0.9

[0116] y 0.38

[0117] It can be seen that X -0.51;

[0118] If it is right coronary dominant type, the intersection is located in the second quadrant, the second quadrant

[0119] The distance from the coaxial origin O to the intersection point is recorded as OX, and then substituted into formula (6) to obtain:

[0120] 0.71;

[0121] LCx stenosis 70%, i.e. z = 0.7

[0122] Then, according to the parabola formula z=ax 2 +n, and substituting the coordinates of the two vertices (0, 0.7) and (0.71, 0) into the parabola formula, we get:

[0123] z=-1.38x 2 +0.7

[0124] Then (0.71, 0) at z=-1.38x 2 The slope of +0.7 (7) is 1.95

[0125] The Ratio of Abundant Coronary flow (RAC) is OX and z is the area of ​​the fourth image domain in this coordinate system ( Figure 2 Orange transparent area in the middle), substituting into formula (8) we get:

[0126] S RAC =0.663

[0127] Then the area of ​​the right triangle with two bases 'S is:

[0128] 'S RAC = 0.249

[0129] Throughout this specification, reference to terms such as "preferred embodiment," "one embodiment," "some embodiments," "example," "specific example," or "some examples" means that the specific features, structures, materials, or characteristics described in conjunction with that embodiment or example are included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

[0130] While embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions, and variations may be made to the embodiments without departing from the principles and spirit of the invention, and that the scope of the invention is defined by the claims and their equivalents.

Claims

1. A coronary artery blood supply evaluation system based on a parabolic three-dimensional coordinate system, characterized in that: It includes a parabola three-dimensional coordinate system, and also includes the intersection point A1 of the parabola corresponding to the blood supply of the left anterior descending coronary artery determined according to the coronary angiography results and the parabola corresponding to the blood supply of the right coronary artery in the X-axis plane, and the vertical coordinate marking the blood supply of the left circumflex coronary artery B; It also involves drawing the area S enclosed by the parabola F1(x) passing through the polar coordinates A1 and B. RAC For the coronary artery abundance area, a polar coordinate A with A1 and A2 coordinate systems is established. The tangent slope of the polar coordinate A can intuitively reflect the coronary blood supply characteristics of the heart. The intersection point A of the parabola corresponding to the left anterior descending coronary artery and the right coronary artery 1， and a virtual symmetrical position A2 of A1 symmetrical with respect to the origin; The characteristics of coronary blood supply of the heart are expressed by the polar coordinates of the A1 and A2 coordinate systems and the slope of the tangent line between A1 and A2; With the blood supply of the left coronary artery circumflex branch marked on the vertical axis B, draw the parabola F1 (x) passing through the polar coordinates A1 and B. RAC , S RAC Extracting the characteristics of the coronary blood supply of the heart; the slope of the parabola F1 (x) pointing from the polar coordinate A1 on the X-axis plane to the ordinate can extract the characteristics of the coronary blood supply of the heart; By simulating the changes in the stenosis of each coronary artery, we can obtain the "coronary artery abundance area" corresponding to different F1 (x), that is, S RAC , so as to estimate the reconstruction results and prognosis before implementing revascularization, select the most optimized revascularization strategy, and thus provide clinical revascularization decision-making; The system includes the following steps: S1, express the coronary angiography results using the common origin equation: The equation for the stenosis rate of the left anterior descending coronary artery is: X=-(N+1)Y 2 +(1-0.1N); The equation for the right coronary artery stenosis rate is: X=(N+1)Y 2 -(1-0.1N); N is determined based on the absolute value of the filling defect at two phases with an angle of 90 degrees on DSA angiography of the same heart's coronary arteries; S2, the origin O is set as the Z axis, which is used to represent the stenosis rate B of the left circumflex coronary artery; S3, complete the establishment of the parabolic three-dimensional coordinate system through S1 and S2; S4, determine the intersection A1 of the left anterior descending coronary artery and the right coronary artery in the X-axis plane; S5, virtual symmetric position A2 of virtual A1 relative to the origin, A2 assists in establishing F1 (X); S6, set A1, A2 and B as the polar coordinates of the parabola F1(X), and B as the vertex coordinate of F1(X); S7, the area enclosed by A1 and B on the parabola F1 (X) is defined as the coronary artery abundance area, which is S RAC ; S8, combined with the parabolic three-dimensional coordinate system to simulate the changes in different coronary artery segments.

2. The coronary artery blood supply evaluation system based on a parabolic three-dimensional coordinate system according to claim 1, characterized in that: The method for establishing the coordinate system in step S1 is as follows: X=(N+1)Y 2 -(1-0.1N); where N is a natural number greater than 0; Then y can be expressed as ; Draw the parabola plane, and then implement the coordinate Z perpendicular to the parabola plane at the origin of the parabola plane to complete the setting of the parabola three-dimensional coordinate system.

3. The coronary artery blood supply evaluation system based on a parabolic three-dimensional coordinate system according to claim 2, characterized in that: Also includes coordinate definitions: Positive x-axis displacement (x) represents stenosis of the left anterior descending coronary artery; Negative x-axis displacement (-x) represents stenosis of the right coronary artery; Positive z-axis displacement (z) represents stenosis of the left circumflex coronary artery.

4. The coronary artery blood supply evaluation system based on a parabolic three-dimensional coordinate system according to claim 3, characterized in that: In the parabolic coordinate system, any X and any -X have two intersection points located symmetrically on both sides of the X axis. The intersection points of the left coronal dominant type are in the first and fourth quadrants, while the intersection points of the right coronal dominant type are in the second and third quadrants.

5. The coronary artery blood supply evaluation system based on a parabolic three-dimensional coordinate system according to claim 4, characterized in that: In step S7, the area S of the image domain enclosed by A1 and B is RAC Calculation methods include ; Where x1 represents the X coordinate of the intersection point A1 of the left anterior descending coronary artery and the right coronary artery in the X-axis plane; x2 represents the X coordinate of the virtual point A2 symmetrical with respect to the origin of A1, satisfying x2=-x1.

Images