Disclosure of Invention
The present application mainly aims to provide a method and an apparatus for naming a coronary artery, a computer storage medium, and a processor, so as to solve the technical problem in the prior art that the accuracy of a scheme for naming a coronary artery is low.
In order to achieve the above object, according to one aspect of the present application, there is provided a method for naming a coronary artery of a heart, the method including: acquiring a 3D image of a heart coronary vessel, and extracting a central line tree of the heart coronary vessel from the 3D image; obtaining a heart atrium and ventricle segmentation model; determining the anatomical structure of the atrioventricular heart according to the heart atrioventricular segmentation model, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane; and naming the coronary artery blood vessel of the heart according to the relative position relation of the central line tree and the anatomical structure of the atrium and the ventricle.
Further, the heart, atrium and ventricle segmentation model comprises a right atrium and right ventricle model and a left atrium and left ventricle model, and the anatomical structure of the atrium and ventricle is determined according to the heart, atrium and ventricle segmentation model, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane, and the method comprises the following steps: expanding the right atrium and right ventricle model by a preset multiple, so that a first intersection is generated between the right atrium and the right ventricle; determining a right atrioventricular groove intersection line according to the first intersection; determining a right atrioventricular groove plane according to the right atrioventricular groove intersection line; expanding the left atrial-left ventricular model by the predetermined factor such that a second intersection is created between the left atrium and the left ventricle; determining a left atrioventricular groove intersection line according to the second intersection; and determining a left atrioventricular groove plane according to the left atrioventricular groove intersection line.
Further, the anatomical structure of the atrioventricular chamber is determined based on said cardiac atrioventricular segmentation model, the method further comprising: acquiring a first central point of the right atrioventricular groove intersecting line and a second central point of the left atrioventricular groove intersecting line; acquiring a first normal vector of a first connecting line of the first central point and the second central point; acquiring a third central point of the right atrioventricular groove plane and a fourth central point of the left atrioventricular groove plane; acquiring a midpoint of a second connecting line of the third central point and the fourth central point; and constructing a cross surface according to the first normal vector and the middle point.
Further, a first center point of the right atrioventricular groove intersection is obtained, the method comprising: acquiring a first average value of distances from any one first target point on the right atrioventricular groove intersecting line to all points except the first target point; and determining the first target point corresponding to the minimum first average value as the first central point.
Further, acquiring a third center point of the right atrioventricular groove plane, the method comprising: acquiring a second average value of distances from any one second target point on the right atrioventricular groove plane to all points except the second target point; and determining the second target point corresponding to the minimum second average value as the third central point.
Further, the anatomical structure of the atrioventricular chamber is determined based on said cardiac atrioventricular segmentation model, the method further comprising: and acquiring the interventricular sulcus intersection line of the left ventricle and the right ventricle.
Further, the centerline tree includes a left crown and a right crown, and the cardiac coronary vessels are named according to the relative position relationship between the centerline tree and the anatomical structure of the atrioventricular chamber, and the method includes: naming a centerline in the right crown that passes through the cruciform intersection and extends furthest rearward as primary PL; acquiring a bifurcation point of a centerline of the centerline tree, except the main PL, and the main PL; acquiring the average value of the minimum distances from all points on a line segment from a bifurcation point to an end point of a preset central line to the intersection line of the left ventricular chamber and the right ventricular chamber; acquiring the minimum value of the distance from all points on a line segment from a bifurcation point to an end point of the preset central line to the intersection line of the left ventricular chamber and the right ventricular chamber; the product of the average value and the minimum value is minimum, and the corresponding predetermined central line is PD.
Further, the coronary artery blood vessel of the heart is named according to the relative position relationship between the central line tree and the anatomical structure of the atrium and the ventricle, and the method further comprises the following steps: the part from the origin of the right corona to the bifurcation of the main PL and the PD is named RCA.
Further, the coronary artery blood vessel of the heart is named according to the relative position relationship between the central line tree and the anatomical structure of the atrium and the ventricle, and the method further comprises the following steps: dividing the RCA into a first front section, a middle section and a first rear section; a portion which is emitted from the first anterior segment and in which the dot product of a second normal vector of a vector composed of a bifurcation point and an end point of the RCA and a third normal vector of the right atrioventricular groove plane is smaller than a predetermined value is named CA; a part which is emitted from the middle section and the dot product of a fourth normal vector of a vector formed by a bifurcation point and an end point of the RCA and a fifth normal vector of the right atrioventricular groove plane is smaller than the preset value is named as Rv; and a portion which is emitted from the first posterior segment and in which the dot product of the sixth normal vector of the vector constituted by the bifurcation point and the end point of the RCA and the seventh normal vector of the right atrioventricular groove plane is smaller than the predetermined value is named AM.
Further, the coronary artery blood vessel of the heart is named according to the relative position relationship between the central line tree and the anatomical structure of the atrium and the ventricle, and the method further comprises the following steps: dividing the RCA into a second front section and a second rear section; a portion which is emitted from the second front segment and in which the dot product of the eighth normal vector of the vector constituted by the bifurcation point and the end point of the RCA and the ninth normal vector of the right atrioventricular groove plane is smaller than a predetermined value is named SAN; and designating a portion which is emitted from the second posterior segment and in which a dot product of a tenth normal vector of a vector composed of a bifurcation point and an end point of the RCA and an eleventh normal vector of the right atrioventricular groove plane is smaller than a predetermined value, as RA Branch.
Further, the centerline tree includes a left crown and a right crown, and the cardiac coronary vessels are named according to the relative position relationship between the centerline tree and the anatomical structure of the atrioventricular chamber, and the method includes: the center line with the minimum score with the left ventricular interventricular sulcus intersection line in the left crown is named as LAD, wherein the score is obtained by dividing the sum of the average value of the minimum distances from all points on the center line to the left ventricular interventricular sulcus intersection line and the minimum value of the distances from all points on the center line to the left ventricular interventricular sulcus intersection line by the length of the center line; obtaining a bifurcation point between a centerline in the left crown, other than the LAD, and the LAD; determining the point with the minimum bifurcation point order as an LM bifurcation point, wherein the bifurcation point order means that all centerlines except the LAD and bifurcation points of the centerline of the LAD are sequenced to obtain a bifurcation point order; and dividing the left crown into an LAD coronary artery tree and an LCX coronary artery tree according to the LM bifurcating point.
Further, the coronary artery blood vessel of the heart is named according to the relative position relationship between the central line tree and the anatomical structure of the atrium and the ventricle, and the method further comprises the following steps: fitting all points on the LAD from the LM bifurcation point to the LAD end point to an LAD plane; calculating the dot product of the vector formed by the central line in the LAD coronary artery tree and the bifurcation point and the tail point of the central line of the LAD and the normal vector of the LAD plane; the central line of which the dot product is smaller than a first preset value is named as D, and the central line of which the dot product is larger than or equal to the first preset value is named as S; and if the central line is D and the distance between the central line and the branch point of the LAD and the LM branch point is less than the preset distance, naming the central line as RM.
Further, the coronary artery blood vessel of the heart is named according to the relative position relationship between the central line tree and the anatomical structure of the atrium and the ventricle, and the method further comprises the following steps: and (3) naming the central line with the minimum score with the left atrioventricular groove in the LCX coronary tree as a main LCX, wherein the scoring rule is as follows: the sum of the average value of the minimum distances from all points on the central line behind the LM bifurcation point to the left atrioventricular groove intersection line and the minimum value of the distances from all points on the central line behind the LM bifurcation point to the left atrioventricular groove intersection line is divided by the length of the central line to obtain the score; calculating the dot product of a vector formed by a central line of the LCX coronary artery tree except the main LCX and a point from a bifurcation point of the main LCX to the tail point of the central line and a normal vector of the left atrioventricular groove plane; naming a centerline with the dot product smaller than a second predetermined value as OM, naming a centerline with the dot product larger than a third predetermined value larger than the second predetermined value as LA Branch, otherwise, naming LCX; in a case where the centerline is OM, and a distance between the centerline and a bifurcation point of the main LCX and the LM bifurcation point is less than a predetermined distance, the centerline is named RM.
According to another aspect of the present application, there is provided a naming apparatus for coronary vessels of the heart, the apparatus comprising: the device comprises a first acquisition unit, a second acquisition unit and a third acquisition unit, wherein the first acquisition unit is used for acquiring a 3D image of a cardiac coronary vessel and extracting a central line tree of the cardiac coronary vessel from the 3D image; the second acquisition unit is used for acquiring a heart atrium and ventricle segmentation model; a determining unit, configured to determine an anatomical structure of an atrioventricular ventricle according to the cardiac atrioventricular segmentation model, where the anatomical structure at least includes an atrioventricular groove intersection line and an atrioventricular groove plane; and the naming unit is used for naming the coronary artery blood vessel of the heart according to the relative position relation of the central line tree and the anatomical structure of the atrium and the ventricle.
According to another aspect of the present application, there is also provided a computer-readable storage medium, which includes a stored program, wherein the program, when executed, controls a device on which the computer-readable storage medium is located to perform the above-mentioned method for naming a coronary artery of a heart.
According to another aspect of the present application, there is also provided a processor for executing a program, wherein the program is executed to perform the above-mentioned naming method for coronary vessels of the heart.
By applying the technical scheme of the application, firstly, a 3D image of the coronary artery blood vessel of the heart is obtained, and a central line tree of the coronary artery blood vessel of the heart is extracted from the 3D image; and obtaining a cardiac atrial ventricular segmentation model; then according to the heart atrium and ventricle segmentation model, determining the anatomical structure of the atrium and ventricle, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane; finally, the coronary vessels of the heart are named according to the relative position relationship between the central line tree and the anatomical structure of the atrioventricular heart. According to the method, the heart coronary artery blood vessel is named according to the relative position relation of the central line tree and the anatomical structure of the atrium ventricle, so that the technical effect that the blood vessel information can be continuously and accurately extracted and named in the complex background of the 3D image of the heart coronary artery blood vessel is achieved, and the technical problem that the accuracy of the scheme for naming the heart coronary artery blood vessel in the prior art is low is solved.
Detailed Description
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present application will be described in detail below with reference to the embodiments with reference to the attached drawings.
In order to make the technical solutions better understood by those skilled in the art, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only partial embodiments of the present application, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
It should be noted that the terms "first," "second," and the like in the description and claims of this application and in the drawings described above are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It should be understood that the data so used may be interchanged under appropriate circumstances such that embodiments of the application described herein may be used. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
It will be understood that when an element such as a layer, film, region, or substrate is referred to as being "on" another element, it can be directly on the other element or intervening elements may also be present. Also, in the specification and claims, when an element is described as being "connected" to another element, the element may be "directly connected" to the other element or "connected" to the other element through a third element.
Example 1
According to an embodiment of the present application, a method for naming a coronary vessel of a heart is provided. It should be noted that the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer-executable instructions and that, although a logical order is illustrated in the flowcharts, in some cases, the steps illustrated or described may be performed in an order different than presented herein.
Fig. 1 is a flow chart of a naming method of coronary vessels of the heart according to an embodiment of the present application. As shown in fig. 1, the method comprises the steps of:
step S101, acquiring a 3D image of the coronary artery blood vessel of the heart, and extracting a center line tree of the coronary artery blood vessel of the heart from the 3D image.
The 3D image of the coronary vessels of the heart in the above step can be obtained by CTA (Computed tomography Angiography), which is a non-invasive examination method for evaluating the vascular system; the centerline tree of the coronary artery blood vessel of the heart may be extracted by using a threshold segmentation method, and further, after the centerline tree is obtained, the centerline in the centerline tree may be extracted by using a skeleton method, so as to obtain centerline information of the coronary artery blood vessel of the heart.
Step S102, a heart atrium and ventricle segmentation model is obtained.
The cardiac atrial ventricular segmentation model in the above steps may be obtained by performing cardiac atrial ventricular extraction and building a cardiac atrial ventricular segmentation model by using an AI segmentation method, and as an example, fig. 3 is a schematic diagram of an alternative cardiac atrial ventricular segmentation model.
Step S103, determining the anatomical structure of the atria and ventricles according to the heart atria and ventricular segmentation model, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane.
The heart, atrium and ventricle segmentation model in the above steps includes a right atrium and right ventricle model and a left atrium and left ventricle model, and the anatomical structure of the atrium and ventricle can be determined by determining a left atrioventricular groove intersection line, a left atrioventricular groove plane, a right atrioventricular groove intersection line, and a right atrioventricular groove plane, for example, as shown in fig. 4, a schematic diagram of an atrioventricular groove intersection line is shown, and as shown in fig. 5, a schematic diagram of an atrioventricular groove plane is shown.
And step S104, naming the coronary artery blood vessels of the heart according to the relative position relation of the central line tree and the anatomical structures of the atria and the ventricles.
The relative positional relationship in the above steps may include an extension line, a bifurcation point, a maximum value or a minimum value of the distance, and the like, but it should be noted that any positional relationship is explained based on the centerline tree and the anatomical structure of the atrioventricular chamber.
By applying the technical scheme of the application, firstly, a 3D image of the coronary artery blood vessel of the heart is obtained, and a central line tree of the coronary artery blood vessel of the heart is extracted from the 3D image; and obtaining a cardiac atrial ventricular segmentation model; then according to the heart atrium and ventricle segmentation model, determining the anatomical structure of the atrium and ventricle, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane; finally, the coronary vessels of the heart are named according to the relative position relationship between the central line tree and the anatomical structure of the atrioventricular heart. According to the method, the heart coronary artery blood vessel is named according to the relative position relation of the central line tree and the anatomical structure of the atrium ventricle, so that the technical effect that the blood vessel information can be continuously and accurately extracted and named in the complex background of the 3D image of the heart coronary artery blood vessel is achieved, and the technical problem that the accuracy of the scheme for naming the heart coronary artery blood vessel in the prior art is low is solved.
Further, the heart, atrium and ventricle segmentation model comprises a right atrium and right ventricle model and a left atrium and left ventricle model, and the anatomical structure of the atrium and ventricle is determined according to the heart, atrium and ventricle segmentation model, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane, and the method comprises the following steps: expanding the right atrium and right ventricle model by a preset multiple, so that a first intersection is generated between the right atrium and the right ventricle; determining a right atrioventricular groove intersection line according to the first intersection; determining a right atrioventricular groove plane according to a right atrioventricular groove intersection line; expanding the left atrium-left ventricle model by a predetermined factor such that a second intersection is created between the left atrium and the left ventricle; determining the left atrioventricular groove intersection line according to the second intersection; and determining the left atrioventricular groove plane according to the left atrioventricular groove intersection line.
Specifically, the right and left atrioventricular groove planes are oriented normal to the atrium.
In an optional embodiment, the right atrium and right ventricle model is first expanded by 1.1 times, so that the right atrium and right ventricle generate an intersection, that is, the first intersection, of the right atrium and right ventricle, the right atrioventricular groove intersection line is obtained, based on all points on the right atrioventricular groove intersection line, a correlation algorithm is used for fitting to determine a right atrioventricular groove plane, wherein the correlation algorithm can use a least square method, the same method is applied to the left atrium and left ventricle, the left atrium and right ventricle model is expanded by 1.1 times, the left atrioventricular groove intersection line is determined according to the obtained second intersection, and the left atrioventricular groove plane is determined by fitting.
Further, the anatomical structure of the atrioventricular chamber is determined based on a cardiac atrioventricular segmentation model, the method further comprising: acquiring a first central point of a right atrioventricular groove intersecting line and a second central point of a left atrioventricular groove intersecting line; acquiring a first normal vector of a first connecting line of the first central point and the second central point; acquiring a third central point of a right atrioventricular groove plane and a fourth central point of a left atrioventricular groove plane; acquiring a midpoint of a second connecting line of the third central point and the fourth central point; and constructing a cross surface according to the first normal vector and the midpoint.
In an alternative embodiment, a unique cross-plane may be determined based on the normal vector of the line connecting the midpoints of the left and right atrioventricular grooves, i.e., the first normal vector, and the midpoint of the line connecting the centers of the left and right atrioventricular grooves, i.e., the midpoint of the second line.
Further, a first center point of a right atrioventricular sulcus intersection is obtained, the method comprising: acquiring a first average value of distances from any first target point on the right atrioventricular groove intersecting line to all points except the first target point; and determining the first target point corresponding to the minimum first average value as a first central point.
In an alternative embodiment, the first central point in the above step, as shown in fig. 5, may be a right atrioventricular groove intersecting line which is an irregular closed curve, for further explanation, as shown in fig. 6, it is assumed that the trajectory shown in the figure is a certain atrioventricular groove intersecting line, theoretically, the atrioventricular groove intersecting line is composed of numerous points, that is, there are numerous first target points, for convenience of explanation, A, B, C points are taken to describe, that is, A, B, C points are the first target points, there are many distances to points other than the first target point at A, B, C points, three distances are respectively selected at A, B, C point in fig. 6 for explanation, the first average value determined by a point in the figure is obtained by calculating the lengths of three solid lines, and similarly, the first average values at B, C point are respectively calculated by three dotted lines starting from B, C, after the first average values of A, B, C are obtained, the average values of the three points are compared, and the smallest average value corresponds to the first central point of the right atrioventricular groove intersection line.
Further, acquiring a third center point of the right atrioventricular groove plane, the method comprising: acquiring a second average value of the distances from any one second target point on the right atrioventricular groove plane to all points except the second target point; and determining a second target point corresponding to the minimum second average value as a third central point.
The determination of the third central point in the above steps is similar to the determination of the first central point and the second central point, except that the atrioventricular groove plane is a plane, that is, when the third central point is determined, if the atrioventricular groove plane is a regular circle, the third central point may be located at the center of the plane, and similarly, the determination methods of the fourth central point and the third central point are the same.
Further, the anatomical structure of the atrioventricular chamber is determined based on a cardiac atrioventricular segmentation model, the method further comprising: and acquiring the interventricular sulcus intersection line of the left ventricle and the right ventricle.
In an alternative embodiment, a schematic view of the interventricular sulcus of the left and right ventricles is shown in FIG. 7.
Further, the centerline tree includes a left crown and a right crown, and the cardiac coronary vessels are named according to the relative position relationship between the centerline tree and the anatomical structures of the atrioventricular chambers, and the method includes: the center line which passes through the cross-shaped cross surface in the right crown and extends furthest backwards is named as a main PL; acquiring a bifurcation point of a central line except the main PL and the main PL in the central line tree; acquiring the average value of the minimum distances from all points on a line segment from a bifurcation point to an end point of a preset central line to a groove intersection line between a left ventricle chamber and a right ventricle chamber; acquiring the minimum value of the distance from all points on a line segment from a bifurcation point to a tail point of a preset central line to a groove intersection line between a left ventricle chamber and a right ventricle chamber; the product of the average and minimum values is the smallest and the corresponding predetermined centerline is PD.
In an alternative embodiment, as shown in fig. 8, fig. 8 is a diagram of the result of naming the coronary vessels in the heart, the cross plane is perpendicular to fig. 8, because of the mirror image relationship of the diagram, the right crown is the left half of the diagram, the main PL can be judged according to the cross plane, the central line of the right crown which passes through the cross plane and extends furthest backwards is the main PL, the central line except the main PL and the bifurcation point of the main PL are calculated, the part from the bifurcation point of the processed central line to the end point is calculated, and the intersection line with the interventricular sulcus scores are the smallest PD, and two bifurcation points, namely PD1 and PD2, can be seen in fig. 8, wherein the scoring rules are as follows: the product of the average value of the minimum distances from each point on the central line behind the bifurcation point to the interchamber groove intersecting line and the minimum value of the distances from all points on the central line to the interchamber groove intersecting line, namely the product of the average value and the minimum value is minimum.
In another alternative embodiment, the dot product of the vector from the centerline bifurcation point to the end point and the right atrioventricular groove plane vector is AVN if it is greater than 0.3, and PL if it is not.
Further, the method includes naming the coronary vessels of the heart according to the relative position relationship between the centerline tree and the anatomical structures of the atria and ventricles, and the method further includes: the part from the origin of the right corona to the bifurcation of the main PL and PD is named RCA.
In an alternative embodiment, based on fig. 8, the part from the origin of the right coronary to the bifurcation of the main PL and PD is named RCA, as shown in fig. 9, and fig. 9 is a schematic diagram of RCA naming the coronary vessels of the heart.
Further, the method includes naming the coronary vessels of the heart according to the relative position relationship between the centerline tree and the anatomical structures of the atria and ventricles, and the method further includes: dividing the RCA into a first front section, a middle section and a first rear section; a part which is emitted from the first front section and the dot product of a second normal vector of a vector formed by a bifurcation point and an end point of the RCA and a third normal vector of a right atrioventricular groove plane is smaller than a preset value is named as CA; a part which is emitted from the middle section and the dot product of the fourth normal vector of the vector formed by the bifurcation point and the end point of the RCA and the fifth normal vector of the right atrioventricular groove plane is smaller than a preset value is named as Rv; and a portion which is emitted from the first posterior segment and in which the dot product of the sixth normal vector of the vector constituted by the bifurcation point and the end point of the RCA and the seventh normal vector of the right atrioventricular groove plane is smaller than a predetermined value is named AM.
In an alternative embodiment, still taking fig. 9 as an example for description, the RCAs are divided into front, middle and rear sections, that is, the first front section, the middle section and the first rear section; the dot product of the normal direction of the vector from the first front section to the tail point of the RCA bifurcation and the normal direction of the right atrioventricular groove plane is less than or equal to 0.5 and is named as CA; the dot product of the normal direction of the vector from the bifurcation point of the middle section to the tail point of the RCA and the normal direction of the right atrioventricular groove plane is less than or equal to 0.5, and is named as Rv, and the diagram comprises Rv1 andRv 2; the dot product of the normal direction of the vector from the bifurcation point to the tail point of the RCA and the normal direction of the right atrioventricular groove plane, which is emitted from the first posterior segment, is less than or equal to 0.5, and the name is AM.
Further, the method includes naming the coronary vessels of the heart according to the relative position relationship between the centerline tree and the anatomical structures of the atria and ventricles, and the method further includes: dividing the RCA into a second front section and a second rear section; a portion which is emitted from the second front segment and the dot product of the eighth normal vector of the vector formed by the bifurcation point and the end point of the RCA and the ninth normal vector of the right atrioventricular groove plane is smaller than a preset value is named as SAN; a portion which is emitted from the second posterior segment and in which the dot product of the tenth normal vector of the vector constituted by the bifurcation point and the end point of the RCA and the eleventh normal vector of the right atrioventricular groove plane is smaller than a predetermined value is named RA Branch.
In an alternative embodiment, the RCA is divided into two segments, namely, the second front segment and the second back segment, and the dot product of the normal direction of the vector from the second front segment to the end point of the RCA bifurcation and the normal direction of the right atrioventricular groove plane is greater than 0.5, which is named as SAN; the dot product of the normal direction of the vector from the bifurcation point to the end point of the RCA and the normal direction of the right atrioventricular groove plane, which is emitted from the second posterior segment, is more than 0.5, and is named as RA Branch.
Further, the centerline tree includes a left crown and a right crown, and the cardiac coronary vessels are named according to the relative position relationship between the centerline tree and the anatomical structures of the atrioventricular chambers, and the method includes: the center line with the minimum score with the left ventricular interventricular sulcus intersection line in the left crown is named as LAD, wherein the score is obtained by dividing the sum of the average value of the minimum distances from all points on the center line to the left ventricular interventricular sulcus intersection line and the minimum value of the distances from all points on the center line to the left ventricular interventricular sulcus intersection line by the length of the center line; acquiring a center line of the center line in the left crown except the LAD and a bifurcation point of the LAD; determining the point with the minimum bifurcation point order as an LM bifurcation point, wherein the bifurcation point order means that all centerlines except the LAD and bifurcation points of the centerline of the LAD are sequenced to obtain a bifurcation point order; the left coronary is divided into LAD coronary tree and LCX coronary tree according to the LM bifurcate point.
In an alternative embodiment, the scoring rule is as follows according to the smallest LAD score between the left coronal centerline and the interventricular sulcus intersection: dividing the sum of the average value of the minimum distances from each point on the central line to the interchamber gully intersection line behind the bifurcation point and the minimum value of the distances from all points on the central line to the interchamber gully intersection line by the length of the central line to obtain a score; and then, obtaining bifurcation points according to the LAD and all the left crown center lines except the LAD, and obtaining the point with the minimum bifurcation point order as an LM bifurcation point, thereby obtaining an LAD coronary artery tree and an LCX coronary artery tree.
Further, the method includes naming the coronary vessels of the heart according to the relative position relationship between the centerline tree and the anatomical structures of the atria and ventricles, and the method further includes: fitting all points on the LAD from the LM bifurcation point to the LAD end point to form an LAD plane; calculating the dot product of the vector formed by the central line in the LAD coronary artery tree and the bifurcation point of the LAD and the tail point of the central line and the normal vector of the LAD plane; a central line with the dot product smaller than a first preset value is named as D, and a central line with the dot product larger than or equal to the first preset value is named as S; if the central line is D and the distance between the central line and the branch point of the LAD and the branch point of the LM is less than the preset distance, the central line is named as RM.
In an alternative embodiment, as shown in fig. 10, fig. 10 is a diagram illustrating the naming result of the coronary artery tree centerline LAD, and a least square method is used to fit a plane, called LAD plane for short, from the center line point of the part starting from the LM bifurcation to the center line point of the last point of the LAD on the LAD. Calculating the point product of the vector formed by the central line of the LAD coronary artery tree and the bifurcation point of the LAD and the tail point of the central line and the normal direction of the LAD plane, wherein D is calculated if the point product is less than 0.1, and D is calculated as D shown in figure 10, D1, D2, D3 and D4, otherwise S is calculated as S shown in figure 10, and S is calculated as S1 and S2, wherein D is calculated and the point product is named as RM if the point product is D and the distance from the bifurcation point of the LAD to the LM is less than 0.5 mm.
Further, the method includes naming the coronary vessels of the heart according to the relative position relationship between the centerline tree and the anatomical structures of the atria and ventricles, and the method further includes: the center line with the minimum score with the intersection line of the left atrioventricular groove in the LCX coronary artery tree is named as a main LCX, and the scoring rule is as follows: the sum of the average value of the minimum distances from all points on the central line behind the LM bifurcation point to the left atrioventricular groove intersection line and the minimum value of the distances from all points on the central line behind the LM bifurcation point to the left atrioventricular groove intersection line is divided by the length of the central line to obtain a score; calculating the dot product of a vector formed by a central line except the main LCX in the LCX coronary artery tree and a bifurcation point of the main LCX to the tail point of the central line and a normal vector of the left atrioventricular groove plane; the center line with the dot product smaller than the second preset value is named as OM, the center line with the dot product larger than a third preset value is named as LA Branch, otherwise, the center line is named as LCX, wherein the third preset value is larger than the second preset value; in the case where the center line is OM, and the distance between the center line and the branch point of the main LCX and the LM branch point is less than the predetermined distance, the center line is named RM.
In an alternative embodiment, as shown in fig. 11, fig. 11 is a diagram illustrating the naming result of the central line LCX of the coronary artery tree, and the scoring rule is as follows according to the main LCX with the smallest score of the intersection line of the central line and the left atrioventricular groove: and dividing the sum of the average value of the minimum distances from each point on the central line to the left atrioventricular groove intersection line behind the LM bifurcation point and the minimum value of the distances from all points on the central line to the left atrioventricular groove intersection line by the length of the central line to obtain a score, judging the bifurcation points of the rest central lines and the main LCX, and obtaining a vector from the bifurcation point to the tail point. The vector is calculated as the dot product of the vector normal to the left atrioventricular groove plane. If the dot product is less than-0.2, it is OM, such as OM1, OM2, OM3, OM4, and OM5 shown in FIG. 11; if the dot product is larger than 0.1, the LA Branch is obtained, otherwise, the LCX is obtained. Wherein, if the name is OM and the distance between the bifurcation point and the LM bifurcation point is less than 0.5mm, the name is RM.
Example 2
The embodiment of the present application further provides a device for naming a coronary artery blood vessel, and it should be noted that the device for naming a coronary artery blood vessel of the present application may be used to execute the method for naming a coronary artery blood vessel of the heart provided by the embodiment of the present application. The following describes a device for naming a coronary artery blood vessel of the heart provided by the embodiment of the present application.
Fig. 12 is a schematic diagram of a naming device for coronary vessels of the heart according to an embodiment of the present application. As shown in fig. 12, the apparatus includes:
a first obtaining unit 1201, configured to obtain a 3D image of a cardiac coronary artery blood vessel, and extract a centerline tree of the cardiac coronary artery blood vessel from the 3D image;
a second obtaining unit 1202 for obtaining a cardiac atrioventricular segmentation model;
a determining unit 1203, configured to determine an anatomical structure of the atrioventricular chamber according to the cardiac atrioventricular segmentation model, where the anatomical structure at least includes an atrioventricular groove intersection line and an atrioventricular groove plane;
a naming unit 1204 for naming the coronary vessels of the heart according to the relative position relationship of the centerline tree and the anatomical structure of the atrioventricular chambers.
By applying the technical scheme of the application, the 3D image of the coronary artery blood vessel of the heart is obtained through the first obtaining unit, and the central line tree of the coronary artery blood vessel of the heart is extracted from the 3D image; acquiring a heart atrium and ventricle segmentation model through a second acquisition unit; determining the anatomical structure of the atrioventricular chamber by a determining unit according to the heart atrioventricular segmentation model, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane; the coronary vessels of the heart are named by the naming unit according to the relative position relationship of the centerline tree and the anatomical structure of the atria and ventricles. According to the method, the heart coronary artery blood vessel is named according to the relative position relation of the central line tree and the anatomical structure of the atrium ventricle, so that the technical effect that the blood vessel information can be continuously and accurately extracted and named in the complex background of the 3D image of the heart coronary artery blood vessel is achieved, and the technical problem that the accuracy of the scheme for naming the heart coronary artery blood vessel in the prior art is low is solved.
Further, the determination unit includes: the first expansion module is used for expanding the right atrium and right ventricle model by a preset multiple, so that a first intersection is generated between the right atrium and the right ventricle; the first determining module is used for determining a right atrioventricular groove intersection line according to the first intersection; the second determining module is used for determining a right atrioventricular groove plane according to a right atrioventricular groove intersecting line; a second expansion module for expanding the left atrium-left ventricle model by a predetermined multiple such that a second intersection is generated between the left atrium and the left ventricle; a third determining module, configured to determine a left atrioventricular groove intersection line according to the second intersection; and the fourth determining module is used for determining the left atrioventricular groove plane according to the left atrioventricular groove intersecting line.
Further, the determination unit further includes: the first acquisition module is used for acquiring a first central point of a right atrioventricular groove intersecting line and a second central point of a left atrioventricular groove intersecting line; the second acquisition module is used for acquiring a first normal vector of a first connecting line of the first central point and the second central point; the third acquisition module is used for acquiring a third central point of the right atrioventricular groove plane and a fourth central point of the left atrioventricular groove plane; the fourth acquisition module is used for acquiring the midpoint of a second connecting line of the third central point and the fourth central point; and the construction module constructs a cross plane according to the first normal vector and the midpoint.
Further, the first obtaining module comprises: the first acquisition submodule is used for acquiring a first average value of distances from any one first target point on the right atrioventricular groove intersecting line to all points except the first target point; and the first determining submodule is used for determining the first target point corresponding to the minimum first average value as the first central point.
Further, the third obtaining module includes: the second acquisition submodule is used for acquiring a second average value of the distances from any one second target point on the right atrioventricular groove plane to all points except the second target point; and the second determining submodule is used for determining a second target point corresponding to the minimum second average value as a third central point.
Further, the determination unit further includes: and the fifth acquisition module is used for acquiring the interventricular sulcus intersection line of the left ventricle and the right ventricle.
Further, the naming unit includes: a first naming module for naming a centerline in the right crown passing through the cruciform baffle and extending furthest backwards as a primary PL; a sixth obtaining module, configured to obtain a centerline of the centerline tree, except for the primary PL, and a bifurcation point with the primary PL; a seventh obtaining module, configured to obtain an average value of minimum distances from all points on a line segment from a bifurcation point to a tail point of a predetermined centerline to a groove intersection line between left and right ventricles; the eighth acquisition module is used for acquiring the minimum value of the distance from all points on a line segment from a bifurcation point to a tail point of a preset central line to a groove intersection line between the left ventricle chamber and the right ventricle chamber; and the second naming module is used for minimizing the product of the average value and the minimum value, and the corresponding preset central line is PD.
Further, the naming unit further includes: and the third naming module is used for naming the part from the starting point of the right crown to the bifurcation point of the main PL and the PD as RCA.
Further, the naming unit further includes: the first segmentation module is used for dividing the RCA into a first front segment, a middle segment and a first rear segment; a fourth naming module for naming a portion, which is emitted from the first front segment and in which a dot product of a second normal vector of a vector formed by a bifurcation point and an end point of the RCA and a third normal vector of a right atrioventricular groove plane is smaller than a predetermined value, as CA; a fifth naming module for naming a portion, which is emitted from the middle section and in which a dot product of a fourth normal vector of a vector formed by a bifurcation point and an end point of the RCA and a fifth normal vector of the right atrioventricular groove plane is smaller than a predetermined value, as Rv; and a sixth naming module for naming a portion, which is emitted from the first posterior segment and in which a dot product of a sixth normal vector of a vector formed by a bifurcation point and an end point of the RCA and a seventh normal vector of the right atrioventricular groove plane is smaller than a predetermined value, as AM.
Further, the naming unit further includes: the second segmentation module is used for dividing the RCA into a second front segment and a second rear segment; a seventh naming module configured to name a portion, which is emitted from the second front segment and in which a dot product of an eighth normal vector of a vector formed by a bifurcation point and an end point of the RCA and a ninth normal vector of the right atrioventricular groove plane is smaller than a predetermined value; and an eighth naming module for naming a portion, which is emitted from the second posterior segment and in which a dot product of a tenth normal vector of a vector formed by a bifurcation point and an end point of the RCA and an eleventh normal vector of the right atrioventricular groove plane is smaller than a predetermined value, as RA Branch.
Further, the naming unit further includes: a ninth naming module, configured to name a center line in the left crown with the smallest score with respect to the left and right ventricular interventricular sulci lines as LAD, where the score is obtained by dividing the sum of the average value of the minimum distances from all points on the center line to the left and right ventricular interventricular sulci lines and the minimum value of the distances from all points on the center line to the left and right ventricular interventricular sulci lines by the length of the center line; the ninth acquisition module is used for acquiring a center line of the center line in the left crown except the LAD and a bifurcation point of the LAD; a fifth determining module, configured to determine a point with a minimum bifurcation point order as an LM bifurcation point, where the bifurcation point order refers to a bifurcation point order obtained by ordering all centerlines other than the LAD and the bifurcation point of the centerline of the LAD; and the classification module is used for dividing the left crown into the LAD coronary artery tree and the LCX coronary artery tree according to the LM branching point.
Further, the naming unit further includes: the fitting module is used for fitting all points from the LM bifurcation point to the LAD end point on the LAD into an LAD plane; the first calculation module is used for calculating a vector formed by a central line in the LAD coronary artery tree and a bifurcation point of the LAD and a tail point of the central line, and a dot product of the vector and a normal vector of an LAD plane; a ninth naming module, configured to name the centerline with the dot product smaller than the first predetermined value as D, and the centerline with the dot product greater than or equal to the first predetermined value as S; and the tenth naming module is used for naming the center line as RM under the condition that the center line is D and the distance between the center line and the branch point of the LAD and the branch point of the LM is less than the preset distance.
Further, the naming unit further includes: the eleventh naming module is used for naming the central line with the minimum score with the left atrioventricular groove in the LCX coronary tree as a main LCX, and the scoring rule is as follows: the sum of the average value of the minimum distances from all points on the central line behind the LM bifurcation point to the left atrioventricular groove intersection line and the minimum value of the distances from all points on the central line behind the LM bifurcation point to the left atrioventricular groove intersection line is divided by the length of the central line to obtain a score; the second calculation module is used for calculating the dot product of a vector formed by a central line of the LCX coronary artery tree except the main LCX and a bifurcation point of the main LCX to the tail point of the central line and a normal vector of the left atrioventricular groove plane; a twelfth naming module for naming a centerline with a dot product smaller than the second predetermined value as OM, and a centerline with a dot product larger than a third predetermined value as LA Branch, otherwise, as LCX, wherein the third predetermined value is larger than the second predetermined value; and the thirteenth naming module is used for naming the center line as RM under the condition that the center line is OM and the distance between the center line and the branching point of the main LCX and the LM branching point is less than the preset distance.
The device for naming the coronary artery blood vessel comprises a processor and a memory, wherein the first acquiring unit, the second acquiring unit, the determining unit, the naming unit and the like are stored in the memory as program units, and the processor executes the program units stored in the memory to realize corresponding functions.
The processor comprises a kernel, and the kernel calls the corresponding program unit from the memory. The kernel can be set to be one or more than one, and the naming of the coronary artery blood vessel of the heart is realized according to the relative position relation of the central line tree and the anatomical structure of the atrium and the ventricle by adjusting the kernel parameters, so that the technical effect that the blood vessel information can be continuously and accurately extracted and named in the complex background of the 3D image of the coronary artery blood vessel of the heart is achieved, and the technical problem that the naming scheme of the coronary artery blood vessel of the heart in the prior art is low in accuracy is solved.
The memory may include volatile memory in a computer readable medium, Random Access Memory (RAM) and/or nonvolatile memory such as Read Only Memory (ROM) or flash memory (flash RAM), and the memory includes at least one memory chip.
The embodiment of the invention provides a computer readable storage medium, which comprises a stored program, wherein when the program runs, a device where the computer readable storage medium is located is controlled to execute the naming method of the coronary artery blood vessel of the heart.
The embodiment of the invention provides a processor, which is used for running a program, wherein the program runs to execute the naming method of the coronary artery blood vessel of the heart.
The embodiment of the invention provides equipment, which comprises a processor, a memory and a program which is stored on the memory and can run on the processor, wherein when the processor executes the program, at least the following steps are realized: step S101, acquiring a 3D image of a cardiac coronary vessel, and extracting a central line tree of the cardiac coronary vessel from the 3D image; step S102, obtaining a heart atrium and ventricle segmentation model; step S103, determining an anatomical structure of an atrium and a ventricle according to the heart atrium and ventricle segmentation model, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane; and step S104, naming the coronary artery blood vessels of the heart according to the relative position relation of the central line tree and the anatomical structures of the atria and the ventricles. The device herein may be a server, a PC, a PAD, a mobile phone, etc.
The present application further provides a computer program product adapted to perform a program of initializing at least the following method steps when executed on a data processing device: step S101, acquiring a 3D image of a cardiac coronary vessel, and extracting a central line tree of the cardiac coronary vessel from the 3D image; step S102, obtaining a heart atrium and ventricle segmentation model; step S103, determining an anatomical structure of an atrium and a ventricle according to the heart atrium and ventricle segmentation model, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane; and step S104, naming the coronary artery blood vessels of the heart according to the relative position relation of the central line tree and the anatomical structures of the atria and the ventricles.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). The memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in the process, method, article, or apparatus that comprises the element.
From the above description, it can be seen that the above-described embodiments of the present application achieve the following technical effects:
1) by applying the technical scheme of the application, firstly, a 3D image of the coronary artery blood vessel of the heart is obtained, and a central line tree of the coronary artery blood vessel of the heart is extracted from the 3D image; and obtaining a cardiac atrial ventricular segmentation model; then according to the heart atrium and ventricle segmentation model, determining the anatomical structure of the atrium and ventricle, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane; finally, the coronary vessels of the heart are named according to the relative position relationship between the central line tree and the anatomical structure of the atrioventricular heart. According to the method, the heart coronary artery blood vessel is named according to the relative position relation of the central line tree and the anatomical structure of the atrium ventricle, so that the technical effect that the blood vessel information can be continuously and accurately extracted and named in the complex background of the 3D image of the heart coronary artery blood vessel is achieved, and the technical problem that the accuracy of the scheme for naming the heart coronary artery blood vessel in the prior art is low is solved.
2) By applying the technical scheme of the application, the 3D image of the cardiac coronary artery blood vessel is obtained through the first obtaining unit, and the central line tree of the cardiac coronary artery blood vessel is extracted from the 3D image; acquiring a heart atrium and ventricle segmentation model through a second acquisition unit; determining the anatomical structure of the atrioventricular chamber by a determining unit according to the heart atrioventricular segmentation model, wherein the anatomical structure at least comprises an atrioventricular groove intersection line and an atrioventricular groove plane; the coronary vessels of the heart are named by the naming unit according to the relative position relationship of the centerline tree and the anatomical structure of the atria and ventricles. According to the method, the heart coronary artery blood vessel is named according to the relative position relation of the central line tree and the anatomical structure of the atrium ventricle, so that the technical effect that the blood vessel information can be continuously and accurately extracted and named in the complex background of the 3D image of the heart coronary artery blood vessel is achieved, and the technical problem that the accuracy of the scheme for naming the heart coronary artery blood vessel in the prior art is low is solved.
The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.