CN112353494A - Ascending aorta, aortic arch and three major branches parametric modeling method thereof - Google Patents
Ascending aorta, aortic arch and three major branches parametric modeling method thereof Download PDFInfo
- Publication number
- CN112353494A CN112353494A CN202011273546.6A CN202011273546A CN112353494A CN 112353494 A CN112353494 A CN 112353494A CN 202011273546 A CN202011273546 A CN 202011273546A CN 112353494 A CN112353494 A CN 112353494A
- Authority
- CN
- China
- Prior art keywords
- section
- point
- blood vessel
- target blood
- aorta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/101—Computer-aided simulation of surgical operations
- A61B2034/105—Modelling of the patient, e.g. for ligaments or bones
Abstract
The invention relates to a parametric modeling method for an ascending aorta, an aortic arch and three major branches thereof, belonging to the field of biological 3D printing modeling and relating to the parametric modeling method for the ascending aorta, the aortic arch and the three major branches thereof based on anatomical features. The method comprises the steps of firstly obtaining CT medical images of a patient, importing data into image analysis software, and adjusting various display parameters. A thresholding function and a region growing function of software are applied, and a target blood vessel section and a peripheral constraint organ tissue mask of a patient are drawn by combining manual correction; the peripheral organs and tissues of the artery of the patient including the sternum, the thoracic vertebra, the superior vena cava, the trachea and the pulmonary artery trunk are segmented, and a three-dimensional model of the target blood vessel segment and the peripheral constrained organ tissues is generated. The method improves the applicability and the service life of the implanted aorta, and can adjust the geometric model of the target blood vessel section by adjusting multiple groups of specific parameters, so that the blood flow design is personalized, the degree of freedom is higher, and the degree of recovery is high.
Description
Technical Field
The invention belongs to the field of biological 3D printing modeling, and relates to a parametric modeling method for ascending aorta, aortic arch and three major branches thereof based on anatomical features.
Background
The aorta is the main trunk of the systemic circulation artery, and serves as the largest blood vessel in the human body, and functions to transport blood to various parts throughout the body. Aortic disorders are largely divided into three types: congenital malformations, acquired physical injury and exacerbation, and inflammation of the aorta leading to injury. In the aspect of treatment, the medicine is often intervened, the stent is used for assisting treatment, and in the serious case, the replacement of the aorta partial section is needed. The available autologous source blood vessels are limited, and the artificial blood vessels on the market can ensure biocompatibility to a certain extent, but still have the defects of poor morphological goodness of fit and difficulty in realizing the co-growth with patients. The application of biological 3D printing technology in this aspect has great advantages, and firstly, the cells are from the patient, so the rejection reaction is avoided, and the problem that the cells cannot be well fused after the operation and cannot grow together with the patient is solved. Secondly, the applicability of the implanted aorta can be improved and the life of the implanted aorta can be prolonged by designing an artery model and applying haemodynamic simulation and optimization. Ascending aorta, aortic arch and its three major branches: the brachiocephalic trunk, left common cervical and left subclavian arteries are part of the aorta and are also the more frequent and replacement components.
There are currently two general categories of modeling for the aorta: an Aorta model with high reduction degree is constructed through CT images, for example, the model is constructed in the research of Aorta zero-stress state modeling with T-line characterization [ J ]. Computational Mechanics, 2019, of Sasaki T, Takizawa K and the like, and simulation analysis is carried out. The model has the advantages that the simulation degree is high, the blood flow condition in the body of a patient can be well restored through the reconstructed model, and doctors are helped to judge the focus position more accurately. But the defect is also obvious, because only the aorta of the patient body is reduced, the sizes and the relative position relation of all parts of the aorta are not adjusted, only the reduction is carried out, but no optimization design is carried out.
Another class of models the Aorta only roughly parametrically based on some basic data, such as M.Cilla and M.Casales, study of A Parametric Model for Studying the Aorta by means of the computing fluidic Dynamics [ J ]. Journal of biomemechanics, 2020:109691, since it is the Fluid analysis that is directed to a particular parameter, only three parameters of aortic inlet diameter, aortic arch width and angle between aortic arch and descending Aorta are involved in the modeling, and the complexity of the aortic arch itself and the constraints of the periphery are not taken into account.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, and provides a parametric modeling method of ascending aorta, aortic arch and three major branches thereof based on anatomical features. The method comprises the steps of firstly, acquiring a CT medical image of a patient, importing data into image analysis software Mimics, and adjusting various display parameters; the method comprises the steps of drawing a target blood vessel section and a peripheral constraint organ tissue mask of a patient by applying a thresholding function and a region growing function of a Mimics and combining manual correction, segmenting peripheral organs and tissues of an artery of the patient, including a sternum, a thoracic vertebra, a superior vena cava, a trachea and a pulmonary artery trunk, generating a three-dimensional model of the target blood vessel section and the peripheral constraint organ tissues, performing fairing treatment on the model by applying 3-matic, guiding the model back to the Mimics, fitting the shapes of a healthy ascending aorta and an aortic arch better and having higher reduction degree.
The technical scheme adopted by the invention is a parametric modeling method for ascending aorta, aortic arch and three major branches thereof, the method comprises the steps of firstly obtaining CT medical images of a patient, importing the data into image analysis software Mimics, and adjusting various display parameters; drawing a target blood vessel section and a peripheral constraint organ tissue mask of a patient by applying a thresholding function and a region growing function of the Mimics and combining manual correction, segmenting peripheral organs and tissues of arteries of the patient, including a sternum, a thoracic vertebra, a superior vena cava, a trachea and a pulmonary artery trunk, generating a three-dimensional model of the target blood vessel section and the peripheral constraint organ tissues, and performing photoaging treatment on the model by applying 3-matic and guiding the model back to the Mimics;
the method comprises the following specific steps:
step 1, importing the acquired CT image data of the patient into Mimics17.0, setting a display mode and an image reference direction, adjusting a display gray threshold range, and performing data processing;
a) analyzing the coronal plane in the image, finding the lowest point of the human scapula, connecting the coronal plane with the line, and analyzing the line in the sagittal plane to correspond to the spinous process part of the seventh thoracic vertebra; a fourth thoracic vertebra is positioned by a seventh thoracic vertebra, a connecting line of two transverse processes of the fourth thoracic vertebra is found, a line segment is taken as a diameter to make a circle to find two transverse process middle points, the middle point is taken as a positioning point to observe a sagittal plane, a bottom point of the fourth thoracic vertebra is found, and the bottom point is taken as the positioning point to position a horizontal plane, wherein the horizontal plane is a horizontal plane where an ascending aorta and an aortic arch junction B00 are located, and the horizontal plane where the aortic arch and a descending aorta junction D00 is located;
b) determining the upper limit and the lower limit of the target blood vessel section by a method of translating the horizontal plane image up and down; according to each positioning horizontal plane, a target blood vessel segment mask is generated by utilizing a thresholding function and a region growing function in combination with manual adjustment, and peripheral organs and tissues of the artery of the patient including a sternum, a thoracic vertebra, a superior vena cava, a trachea and a pulmonary artery trunk are segmented, so that a three-dimensional model of the target blood vessel segment and peripheral constrained organ tissues is generated; using 3-matic to carry out smooth treatment on the model and leading the model back to the Mimics;
step 2, drawing mark points on the target blood vessel section with the divided lesion to generate a central line, and performing cutting treatment;
respectively obtaining the cross-sectional areas of corresponding positions of an ascending aorta starting point mark point A00, an ascending aorta and aortic arch junction mark point B00, three segmented branch end point mark points E00, F00 and G00, three large branch starting mark points E01, F01 and G01, an aortic arch and descending aorta junction mark point C00, an aortic arch vertex mark point D00 and mark points A00, B00, C00, D00, E00, F00 and G00 according to the positioning horizontal plane obtained in the step 1;
the initial central line graph of the target blood vessel section comprises the central line of the target blood vessel section of the patient and the extracted mark points, the central line of the lesion blood vessel is converted into a point set, and the mark points and the central line are derived in an igs format;
generating an aortic arch center line by applying a curve-art spline in UG10.0 and a circular arc command; respectively determining corresponding model mark points A0, B0, C0, D0, E, F, G, E0, F0 and G0 from the initial centerline mark points and the initial centerline point set of the target blood vessel segment; drawing tangential direction line segments A0a2, C1C0, D1D0, E1E0, F1F0, G1G0 tangent to the respective endpoints; fitting by taking the central line of the initial target blood vessel section as a target by adjusting the distance between the control point and the end point; generating ascending aorta according to the tangent line segment A1A0 and the points A0 and B0
According to the ascending aorta
Generating a tangential direction line segment B1B 0; designating a junction BO of an ascending aorta and an aortic arch and a vertex CO of the aortic arch as spline endpoints, and finally obtaining a Bezier curve BOCO of an ascending section of the aortic arch; similarly, obtaining a Bezier curve CODO of the aortic arch descending segment;
step 4, constructing the central lines of three branch head-arm trunk branch sections EE0, a left neck total branch section FF0 and a left subclavian artery branch section GG 0;
taking the starting point E and the end point EO of the head-arm stem as end points, and taking a line segment E1E0 in the direction of the end point EO of the head-arm stem as a tangent to make a Bessel curve EE0 which is the central line of the head-arm stem segment; taking the left neck total starting point F and the end point FO as end points, and taking a line segment F1F0 in the direction of the left neck total end point as a tangent to make a Bezier curve FF0 which is a central line of the left neck total segment; taking a starting point G and a terminal point GO of the left subclavian artery as end points, and taking a line segment G1G0 in the terminal direction of the left subclavian artery as a tangent to make a Bessel curve GG0 which is the central line of the left subclavian artery segment;
respectively making reference planes which pass through the starting point of the ascending aorta, the junction point of the ascending aorta and the aortic arch, the vertex of the aortic arch, the junction point of the aortic arch and the descending aorta and the end points of the three branches and are vertical to the tangent direction line segments of the corresponding points; taking the ascending aorta starting point as an example, drawing a sketch on a reference plane at the ascending aorta starting point A0, taking the ascending aorta starting point as the center of a circle, and using the section information obtained in the step 2 as a circle to finish the sketch, so that the sketch is taken as an ascending aorta starting point A0 section circle; in the same way, the section circles of other endpoints can be obtained;
step 6, constructing a fluid domain model of the main section of the target blood vessel section;
extracting the central line of the newly constructed main section of the target blood vessel section by adopting an insertion-correlation replication-extraction geometric characteristic command: is composed of ascending aorta central line, aortic arch ascending section and descending section central line; constructing a main section fluid domain model of the target blood vessel section by adopting a curved surface-sweep command; sequentially selecting an ascending aorta starting point cross-sectional circle, an ascending aorta and aorta arch intersection point cross-sectional circle, an aorta arch vertex cross-sectional circle and an aorta arch and descending aorta intersection point cross-sectional circle as cross-sectional curves, and paying attention to selecting and adding a new set at a middle point; selecting the central line of the main section of the target blood vessel section as a guide line, setting the body type as an entity, clicking to determine to sweep, and obtaining a three-dimensional model of the fluid domain of the main section of the target blood vessel section;
step 7, constructing a fluid domain model of the branch section of the target blood vessel section, and finally obtaining a solid domain three-dimensional model of the target blood vessel section;
constructing a fluid domain model of a branch section of a target blood vessel section by adopting a curved surface-sweep command, selecting a central line of a head-arm trunk as a guide line, selecting a terminal point section circle of the head-arm trunk as a section curve, and setting the body type as an entity to sweep to obtain a three-dimensional fluid domain model of the head-arm trunk; obtaining the fluid domain three-dimensional models of the other two branches in the same way; combining the three-dimensional model of the fluid domain of the three branches and the three-dimensional model of the fluid domain of the main section of the target blood vessel section into a whole by applying a summation command;
and finally, thickening the non-end surface of the fluid domain to obtain a solid domain three-dimensional model of the target blood vessel section.
The invention has the beneficial effects that the modeling method is based on anatomical characteristics, and adopts hemodynamics simulation and optimization, thereby improving the applicability of the implanted aorta and prolonging the service life of the implanted object. The geometric model of the target blood vessel section can be adjusted by adjusting multiple groups of specific parameters, so that the blood flow design is personalized and the degree of freedom is high. The key parameters of the geometric model of the target blood vessel section are determined by introducing the constrained organ tissues and the peripheral blood vessel data, so that the model is more suitable for the form of the healthy ascending aorta on the premise of fully considering the individual specificity, and the reduction degree is higher.
Drawings
FIG. 1 is a flow chart of the parametric modeling method of the ascending aorta, the aortic arch and three major branches thereof.
FIG. 2 is a key cross-section applied during segmentation of a target vessel segment; in FIG. 2A, the center of the cross line is the sagittal plane intercept point corresponding to the line connecting the inferior angles of the two scapulae, from which the seventh thoracic vertebra can be determined; FIG. 2B is the fourth thoracic midsagittal plane defined by the transverse processes of the vertebral bodies, with the transverse line representing the fourth thoracic vertebral body nadir.
FIG. 3 is a diagram of a constrained tissue organ model; wherein, 1-sternum, 2-superior vena cava, 3-pulmonary trunk, 4-trachea, 5-thoracic vertebra.
FIG. 4 is an initial centerline map and marker points of a target vessel segment; wherein, A00-ascending aorta starting point mark point, B00-ascending aorta and aortic arch intersection point mark point, C00-aortic arch and descending aorta intersection point mark point, D00-apex mark point of aortic arch, E01-head arm stem starting point mark point, E00-head arm stem end point mark point, F01-left neck total starting point mark point, F00-left neck total end point mark point, G01-left subclavian artery starting point mark point, G00-left subclavian artery end point mark point.
Fig. 5 is a cross-sectional view of the centerline, tangent line segment and segment circle of the model target vessel segment, wherein,
-literAortic centerline, B0C0 bezier curve-aortic arch ascending segment centerline, C0D0 bezier curve-aortic arch descending segment centerline, bezier curve EE 0-brachiocephalic trunk centerline, bezier curve FF 0-left neck total centerline, bezier curve GG 0-left subclavian artery centerline; A0A 2-ascending aorta starting point tangent line segment, BOB 1-aortic arch starting point tangent line segment, EOE 1-brachiocephalic trunk terminal point tangent line segment, FOF 1-left neck total terminal point tangent line segment, GOG 1-left subclavian artery terminal point tangent line segment, COC 1-aortic arch descending segment starting point tangent line segment, and DOD 1-aortic arch descending segment terminal point direction segment tangent line segment.
FIG. 6 is a three-dimensional model of the ascending aorta, aortic arch and the internal fluid domains of its three major branches.
Fig. 7 is a three-dimensional model diagram of the ascending aorta, the aortic arch and the three major branches finally established. Wherein 7A is a front view of the three-dimensional model, 7B is a right view of the three-dimensional model, 7C is a top view of the three-dimensional model, and 7D is an isometric view of the three-dimensional model.
Detailed Description
The invention is further explained in detail with reference to the drawings and technical solutions.
The method is based on the following information of human anatomy:
A. the aorta rises from the left ventricle to the front, the upper part and the right part, turns to the left, the rear and the upper part, and descends along the front of the spine after turning to the left side of the lower edge of the 4 th thoracic vertebral body.
B. The peripheral constraints of the ascending aorta, the aortic arch and its three major branches are: sternum, thoracic vertebrae, superior vena cava, trachea, and pulmonary trunk.
C. The complexity of the aortic arch and the surrounding constraint conditions are fully considered, an ascending aorta, the aortic arch and three-branch parametric models thereof are constructed, a Bessel curve with high degree of freedom is adopted, the control degree of the central line of the aorta can be ensured, and the defects of the prior art can be effectively overcome.
In the examples, image analysis software micics 17.0 and UG10.0 were used.
Firstly, importing the obtained medical image data of the patient into Mimics17.0, setting a display mode, setting an image reference direction and adjusting a display gray threshold range, wherein the actual gray threshold range in the case is 1-65535;
then, drawing a target blood vessel section and a peripheral constraint organ tissue mask of a patient by applying a thresholding function and a region growing function of the Mimics and combining manual correction, generating a three-dimensional model of the target blood vessel section and the peripheral constraint organ tissue, and performing smooth treatment on the model by applying 3-matic and guiding the model back to the Mimics; a three-dimensional model of peripherally constrained organ tissue is shown in fig. 3;
the target blood vessel segmentation step is as follows:
step 1, analyzing a coronal plane in the image, finding the lowest point of the human scapula, and connecting lines. This connection line is used for analysis in the sagittal plane, corresponding to the spinous process portion of the seventh thoracic vertebra, as shown in fig. 2A; a fourth thoracic vertebra is positioned by a seventh thoracic vertebra, a connecting line of two transverse processes of the fourth thoracic vertebra is found, a circle is made by taking the line segment as the diameter to find two transverse process midpoints, the middle point is taken as a positioning point to observe a sagittal plane, and a bottom point of the fourth thoracic vertebra is found; the horizontal plane is located by using the point as a locating point, which is the horizontal plane where the junction of the ascending aorta and the aortic arch is located, as shown by the white line in fig. 2B, and which is also the horizontal plane where the junction of the aortic arch and the descending aorta is located.
And taking the horizontal plane corresponding to the bottom point of the fourth thoracic vertebra cone as a reference plane, continuously lowering the reference plane and observing the horizontal plane image of part of the ascending aorta until the part of the ascending aorta loses the arc smoothness of the part of the ascending aorta, and taking the part of the ascending aorta as the horizontal plane where the starting point of the ascending aorta is located. And taking the bottom point of the fourth thoracic vertebra cone corresponding to the horizontal plane as a reference plane, continuously lifting the reference plane and observing the horizontal plane image of the branch part of the head and arm trunk of the aorta until the horizontal plane image of the head and arm trunk begins to be divided into two parts, and taking the two parts as the horizontal planes of the end points of the three major branches.
And generating a target blood vessel segment mask by utilizing a thresholding function and a region growing function according to each positioning horizontal plane and combining manual adjustment and correction, and segmenting peripheral organs and tissues of the artery of the patient, including a sternum, a thoracic vertebra, a superior vena cava, a trachea and a pulmonary artery trunk, so as to generate a three-dimensional model of the target blood vessel segment and peripheral constrained organ tissues.
And 2, generating a central line according to the three-dimensional model of the target blood vessel section, and performing cutting treatment, as shown in fig. 4.
According to the positioning horizontal plane obtained in the step 1, starting point mark points A00 of the ascending aorta, section areas 1078.44 square millimeters, connecting point mark points B00 of the ascending aorta and the aortic arch and section areas 963.05 square millimeters are respectively obtained, head-arm trunk, left neck trunk and left subclavian artery three-large branch starting point mark points E01, F01 and G01, end point mark points E00, F00 and G00 are obtained through segmentation, the terminal point section areas of the head-arm trunk are 154.05 square millimeters, the total terminal point section areas of the left neck are 56.17 square millimeters, the section areas of the left subclavian artery are 95.45 square millimeters, the section areas of the aortic arch and descending aorta connecting point mark points C00 are 510.34 square millimeters, and the section areas of the top point mark points D00 and D00 of the aortic arch are 442.06 square millimeters.
Deriving the generated initial central line of the target blood vessel section and each mark point in an igs format; and (5) deriving the constrained tissue and organ three-dimensional models obtained by segmentation in an stl format. And converting the search and replacement functions of the exported. stl files and. igs files application notepad into a. dat point set format which can be accepted by UG, and importing UG 10.0.
Step 4 uses the arc command in UG10.0 with ends A0 and B0 and tangent line A0A2
According to
A tangential direction line segment B1B0 is generated. Selecting type according to the endpoint by using a curve-art spline command, wherein the selection times are 3 times; appointing the junction B0 of ascending aorta and aortic archThe vertex C0 is used as a spline endpoint, a continuous type is tangent G1 in a constraint and is appointed to be tangent to a respective direction line segment; at the moment, the control points of the spline curve are limited on the straight line corresponding to the corresponding direction line segment; and fitting with the central line of the initial target blood vessel section by adjusting the distance between the control point and the end point to construct an initial model. And then further adjustment can be carried out, and the finally obtained curve is a Bezier curve B0C0 of the ascending segment of the aortic arch. Similarly, a bezier curve C0D0 of the aortic arch descending segment can be obtained.
taking the starting point E and the end point EO of the head-arm stem as end points, and taking a line segment E1E0 in the direction of the end point EO of the head-arm stem as a tangent to make a Bessel curve EE0 which is the central line of the head-arm stem segment; taking the left neck total starting point F and the end point FO as end points, and taking a line segment F1F0 in the direction of the left neck total end point as a tangent to make a Bezier curve FF0 which is a central line of the left neck total segment; taking a starting point G and a terminal point GO of the left subclavian artery as end points, and taking a line segment G1G0 in the terminal direction of the left subclavian artery as a tangent to make a Bessel curve GG0 which is the central line of the left subclavian artery segment;
reference planes perpendicular to direction line segments of A0, B0, C0, D0, E0 and F0 and G0 and corresponding points are respectively drawn; taking the point a0 as an example, drawing a sketch on the reference plane at the point a0, taking the starting point of the ascending aorta on the line segment in the corresponding direction as the center of the circle, and using the cross-sectional information obtained in the step 2 to draw a circle, thereby completing the sketch, which is the cross-sectional circle at the point a 0. Similarly, the remaining cross-sectional circles are made, and the obtained set of the center line and the cross-sectional circle is shown in fig. 5;
step 6, constructing a three-dimensional model of the fluid domain of the main section of the target blood vessel section;
extracting newly constructed central lines A0B0, B0C0 and C0D0 of the main section of the target blood vessel section by adopting an insertion-correlation replication-extraction geometric feature command; constructing a main section fluid domain model of the target blood vessel section by adopting a curved surface-sweep command: sequentially selecting an AO point cross-sectional circle, a BO point cross-sectional circle, a CO point cross-sectional circle and a DO point cross-sectional circle as cross-sectional curves, and selecting a curve section AODO as a guide line; setting the type of the body as an entity, clicking to determine to sweep, and obtaining a three-dimensional model of a main section fluid domain of the target blood vessel section;
and 7, constructing a target blood vessel section branch section fluid domain model and a target blood vessel section blood three-dimensional model.
Constructing a fluid domain model of a branch section of a target blood vessel section by adopting a curved surface-sweep command; selecting an E0 point cross-section circle as a cross-section curve, selecting B0E0 as a guide line, setting the type of the body as an entity to sweep, and obtaining a three-dimensional model of a fluid domain of the head-arm stem; obtaining the fluid domain three-dimensional models of the other two branches in the same way; combining the three-dimensional model of the fluid domain of the three branches with the three-dimensional model of the blood of the main section of the target blood vessel section into a whole by applying a summation command; the resulting three-dimensional model of blood is shown in FIG. 6;
the obtained surface chain is thickened by applying a thickening command, and fig. 7 is a three-dimensional model diagram of the ascending aorta, the aortic arch and three branches thereof which are finally established. Wherein 7A is a front view of the three-dimensional model, 7B is a right view of the three-dimensional model, 7C is a top view of the three-dimensional model, and 7D is an isometric view of the three-dimensional model.
Claims (1)
1. A parametric modeling method for ascending aorta, aortic arch and three major branches thereof is characterized in that the method comprises the steps of firstly obtaining CT medical images of a patient, importing data into image analysis software Mimics, and adjusting various display parameters; drawing a target blood vessel section and a peripheral constraint organ tissue mask of a patient by applying a thresholding function and a region growing function of the Mimics and combining manual correction, segmenting peripheral organs and tissues of arteries of the patient, including a sternum, a thoracic vertebra, a superior vena cava, a trachea and a pulmonary artery trunk, generating a three-dimensional model of the target blood vessel section and the peripheral constraint organ tissues, and performing photoaging treatment on the model by applying 3-matic and guiding the model back to the Mimics;
the method comprises the following specific steps:
step 1, importing acquired CT image data of a patient into a Mimics17.0, setting a display mode and an image reference direction, adjusting a display gray threshold range, and performing data processing;
a) analyzing the coronal plane in the image, finding the lowest point of the human scapula, connecting the coronal plane with the line, and analyzing the line in the sagittal plane to correspond to the spinous process part of the seventh thoracic vertebra; a fourth thoracic vertebra is positioned by a seventh thoracic vertebra, a connecting line of two transverse processes of the fourth thoracic vertebra is found, a line segment is taken as a diameter to make a circle to find two transverse process middle points, the middle point is taken as a positioning point to observe a sagittal plane, a bottom point of the fourth thoracic vertebra is found, and the bottom point is taken as the positioning point to position a horizontal plane, wherein the horizontal plane is a horizontal plane where an ascending aorta and an aortic arch junction B00 are located, and the horizontal plane where the aortic arch and a descending aorta junction D00 is located;
b) determining the upper limit and the lower limit of the target blood vessel section by a method of translating the horizontal plane image up and down; according to each positioning horizontal plane, a target blood vessel segment mask is generated by utilizing a thresholding function and a region growing function in combination with manual adjustment, and peripheral organs and tissues of the artery of the patient including a sternum, a thoracic vertebra, a superior vena cava, a trachea and a pulmonary artery trunk are segmented, so that a three-dimensional model of the target blood vessel segment and peripheral constrained organ tissues is generated; using 3-matic to carry out smooth treatment on the model and leading the model back to the Mimics;
step 2, marking points are drawn on the target blood vessel section with the segmented lesion, an initial central line is generated, and cutting processing is carried out;
respectively obtaining the cross-sectional areas of corresponding positions of an ascending aorta starting point mark point A00, an ascending aorta and aortic arch junction mark point B00, three segmented branch end point mark points E00, F00 and G00, three large branch starting mark points E01, F01 and G01, an aortic arch and descending aorta junction mark point C00, an aortic arch vertex mark point D00 and mark points A00, B00, C00, D00, E00, F00 and G00 according to the positioning horizontal plane obtained in the step 1;
the initial central line graph of the target blood vessel section comprises the central line of the target blood vessel section of the patient and the extracted mark points, the central line of the lesion blood vessel is converted into a point set, and the mark points and the central line are derived in an igs format;
step 3, deriving the generated initial central line of the target blood vessel section and the constrained tissue and organ three-dimensional models obtained by dividing each mark point, and introducing the derived constrained tissue and organ three-dimensional models into modeling software UG in a point set form;
generating ascending aorta, aortic arch and three branch center lines by applying curve-art splines and arc commands in UG 10.0; respectively determining marker points A0, B0, C0, D0, E, F, G, E0, F0 and G0 corresponding to the model according to the initial centerline marker points and the initial centerline point set of the target blood vessel segment; drawing tangential direction line segments A0a2, C1C0, D1D0, E1E0, F1F0, G1G0 tangent to the respective endpoints; fitting by taking the central line of the initial target blood vessel section as a target by adjusting the distance between the control point and the end point; generating ascending aorta according to the tangent line segment A1A0 and the points A0 and B0
According to the ascending aorta
Generating a tangential direction line segment B1B 0; designating a junction BO of an ascending aorta and an aortic arch and a vertex CO of the aortic arch as spline endpoints, and finally obtaining a Bezier curve BOCO of an ascending section of the aortic arch; similarly, obtaining a Bezier curve CODO of the aortic arch descending segment;
step 4, constructing the center lines of three branch head-arm trunk branch sections B0E0, a left neck total branch section F0C0 and a left subclavian artery branch section G0C0 of the model;
taking the starting point E and the end point EO of the head-arm stem as end points, and taking a line segment E1E0 in the direction of the end point EO of the head-arm stem as a tangent to make a Bessel curve EE0 which is the central line of the head-arm stem segment; taking the total starting point F and the end point FO of the left neck as end points, and taking a line segment F1F0 in the direction of the total end point of the left neck as a tangent to make a Bezier curve FFO which is a central line of the total section of the left neck; taking a starting point G and a terminal point GO of the left subclavian artery as end points, and taking a line segment G1G0 in the terminal direction of the left subclavian artery as a tangent to make a Bessel curve GG0 which is the central line of the left subclavian artery segment;
step 5, constructing each cross-section circle;
respectively making reference planes which pass through the ascending aorta starting point, the ascending aorta and aorta arch junction, the aortic arch vertex, the aorta arch and descending aorta junction and the terminal points of the three branches and are vertical to the direction line segments of the corresponding points; taking the ascending aorta starting point as an example, drawing a sketch on a reference plane at the ascending aorta starting point A0, taking the ascending aorta starting point as the center of a circle, and using the section information obtained in the step 2 as a circle to finish the sketch, so that the sketch is taken as an ascending aorta starting point A0 section circle; in the same way, the section circles of other endpoints can be obtained;
step 6, constructing a fluid domain model of the main section of the target blood vessel section;
extracting the central line of the newly constructed main section of the target blood vessel section by adopting an insertion-correlation replication-extraction geometric characteristic command: is composed of ascending aorta central line, aortic arch ascending section and descending section central line; constructing a main section fluid domain model of the target blood vessel section by adopting a curved surface-sweep command; sequentially selecting an ascending aorta starting point cross-sectional circle, an ascending aorta and aorta arch intersection point cross-sectional circle, an aorta arch vertex cross-sectional circle and an aorta arch and descending aorta intersection point cross-sectional circle as cross-sectional curves, and paying attention to selecting and adding a new set at a middle point; selecting the central line of the main section of the target blood vessel section as a guide line, setting the body type as an entity, clicking to determine to sweep, and obtaining a three-dimensional model of the fluid domain of the main section of the target blood vessel section;
step 7, constructing a fluid domain model of the branch section of the target blood vessel section, and finally obtaining a solid domain three-dimensional model of the target blood vessel section;
constructing a fluid domain model of a branch section of a target blood vessel section by adopting a curved surface-sweep command, selecting a central line of a head-arm trunk as a guide line, selecting a terminal point section circle of the head-arm trunk as a section curve, and setting the body type as an entity to sweep to obtain a three-dimensional fluid domain model of the head-arm trunk; obtaining the fluid domain three-dimensional models of the other two branches in the same way; combining the three-dimensional model of the fluid domain of the three branches and the three-dimensional model of the fluid domain of the main section of the target blood vessel section into a whole by applying a summation command;
and finally, thickening the non-end surface of the fluid domain to obtain a solid domain three-dimensional model of the target blood vessel section.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011273546.6A CN112353494B (en) | 2020-11-14 | 2020-11-14 | Ascending aorta, aortic arch and three major branches parametric modeling method thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011273546.6A CN112353494B (en) | 2020-11-14 | 2020-11-14 | Ascending aorta, aortic arch and three major branches parametric modeling method thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112353494A true CN112353494A (en) | 2021-02-12 |
| CN112353494B CN112353494B (en) | 2021-11-19 |
Family
ID=74514829
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011273546.6A Active CN112353494B (en) | 2020-11-14 | 2020-11-14 | Ascending aorta, aortic arch and three major branches parametric modeling method thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112353494B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120082363A1 (en) * | 2010-09-30 | 2012-04-05 | Marcus Pfister | Method for Displaying a Vessel of a Particular Biological Subject |
| US20170027648A1 (en) * | 2012-05-14 | 2017-02-02 | Heartflow, Inc. | Method and system for providing information from a patient-specific model of blood flow |
| CN107292928A (en) * | 2017-06-16 | 2017-10-24 | 沈阳东软医疗系统有限公司 | A kind of method and device of blood vessel positioning |
| CN109700527A (en) * | 2019-01-03 | 2019-05-03 | 北京理工大学 | Method for establishing model, model and the simulation operation detection method of dissection of aorta |
-
2020
- 2020-11-14 CN CN202011273546.6A patent/CN112353494B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120082363A1 (en) * | 2010-09-30 | 2012-04-05 | Marcus Pfister | Method for Displaying a Vessel of a Particular Biological Subject |
| US20170027648A1 (en) * | 2012-05-14 | 2017-02-02 | Heartflow, Inc. | Method and system for providing information from a patient-specific model of blood flow |
| CN107292928A (en) * | 2017-06-16 | 2017-10-24 | 沈阳东软医疗系统有限公司 | A kind of method and device of blood vessel positioning |
| CN109700527A (en) * | 2019-01-03 | 2019-05-03 | 北京理工大学 | Method for establishing model, model and the simulation operation detection method of dissection of aorta |
Non-Patent Citations (5)
| Title |
|---|
| TAISIYA SIGAEVA: "Anisotropic residual stresses in arteries", 《ROYAL SOCIETY PUBLISHING》 * |
| TAKAFUMI SASAKI: "Aorta zero-stress state modeling with T-spline discretization", 《COMPUTATIONAL MECHANICS》 * |
| TAKAFUMI SASAKI: "Medical-image-based aorta modeling with zero-stress-state estimation", 《COMPUTATIONAL MECHANICS》 * |
| 付文宇: "基于CT图像的胸主动脉瘤模型数值模拟", 《科学导报》 * |
| 郭子义: "动态三维血管建模法分析Stanf_省略_型胸主动脉夹层四维相位对比MRI", 《第二军医大学学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112353494B (en) | 2021-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9672302B2 (en) | Producing a three-dimensional model of an implant | |
| CN105748177B (en) | A kind of personalized implantable spinal prosthesis and its manufacture method with bionic micropore | |
| CN110738681A (en) | automatic pedicle screw operation path planning method based on deep learning network | |
| US20180039754A1 (en) | Methods and systems for producing an implant | |
| CN103690278A (en) | 3D (three-dimensional) printing technology-based preparation method for personalized bionic vertebral column system | |
| EP2583255B1 (en) | Method for determining bone resection on a deformed bone surface from few parameters | |
| CN107274389A (en) | Femur and Acetabular dissection parameter acquiring method based on CT three-dimensional series images | |
| CN112102291A (en) | Method for obtaining reference data of middle-of-surface defect target by anatomical feature point matching | |
| CN109961436A (en) | A kind of median plane construction method based on artificial nerve network model | |
| CN109846580A (en) | Personalized range bone prosthese with biological function and preparation method thereof | |
| CN116982077A (en) | System and method for automatically segmenting patient-specific anatomy for pathology-specific measurements | |
| Tel et al. | Computer-guided in-house cranioplasty: establishing a novel standard for cranial reconstruction and proposal of an updated protocol | |
| CN112353494B (en) | Ascending aorta, aortic arch and three major branches parametric modeling method thereof | |
| CN206761800U (en) | A kind of human skull's defect repair hone lamella | |
| Mandolini et al. | A procedure for designing custom-made implants for forehead augmentation in people suffering from Apert syndrome | |
| CN110111354A (en) | Liver medical image arteriosomes dividing method and system | |
| CN114191075A (en) | Rapid construction method and system of personalized knee joint prosthesis model | |
| Ye et al. | Automatic design and fabrication of a custom ocular prosthesis using 3D volume difference reconstruction (VDR) | |
| CN108670362B (en) | Cutting guide plate prepared by digital space reconstruction and 3D printing technology | |
| CN109830157A (en) | A kind of production method of 5 sacrum of waist, 1 anterior stabilization model | |
| CN112017275B (en) | Auxiliary positioning method for face midsagittal reference plane | |
| CN108354651B (en) | Cutting guide for sizing a replacement flap for a soft tissue surface | |
| CN110390148A (en) | A kind of incidence soft tissue defects reconstruction operations template design method and surgical template based on curved surface features classification processing and expansion | |
| Medvedev | Construction of complex three-dimensional structures of the aorta of a particular patient using finite analytical formulas | |
| Ghosh et al. | Cranioplast fabrication in a comatose patient: A clinical report |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |