CN114224484A - Method, device, computer equipment and storage medium for locating surface area of intracranial aneurysm - Google Patents
Method, device, computer equipment and storage medium for locating surface area of intracranial aneurysm Download PDFInfo
- Publication number
- CN114224484A CN114224484A CN202111521343.9A CN202111521343A CN114224484A CN 114224484 A CN114224484 A CN 114224484A CN 202111521343 A CN202111521343 A CN 202111521343A CN 114224484 A CN114224484 A CN 114224484A
- Authority
- CN
- China
- Prior art keywords
- aneurysm
- model
- dimensional
- intracranial
- plane
- 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
- 201000008450 Intracranial aneurysm Diseases 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 43
- 206010002329 Aneurysm Diseases 0.000 claims abstract description 164
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 57
- 206010028980 Neoplasm Diseases 0.000 claims abstract description 36
- 238000007917 intracranial administration Methods 0.000 claims abstract description 35
- 230000017531 blood circulation Effects 0.000 claims abstract description 25
- 238000002583 angiography Methods 0.000 claims abstract description 24
- 238000013507 mapping Methods 0.000 claims abstract description 21
- 230000002792 vascular Effects 0.000 claims abstract description 12
- 238000004590 computer program Methods 0.000 claims description 10
- 210000001367 artery Anatomy 0.000 claims description 5
- 230000004807 localization Effects 0.000 claims description 4
- 239000008280 blood Substances 0.000 claims description 3
- 210000004369 blood Anatomy 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000000877 morphologic effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 230000000004 hemodynamic effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 208000032170 Congenital Abnormalities Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/501—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the head, e.g. neuroimaging or craniography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- 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
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/404—Angiography
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Physiology (AREA)
- Neurology (AREA)
- Neurosurgery (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Vascular Medicine (AREA)
- Robotics (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention discloses a method, a device, computer equipment and a storage medium for positioning the surface area of intracranial aneurysm, comprising the following steps: acquiring an intracranial angiography image, and constructing a three-dimensional aneurysm vascular model based on the intracranial angiography image; segmenting the three-dimensional aneurysm blood vessel model along the neck of the aneurysm to obtain a three-dimensional aneurysm model, and carrying out UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model; determining a farthest point and a highest point of the aneurysm along the blood flow direction on a curve of a tumor neck boundary space on a three-dimensional aneurysm blood vessel model, taking a mapping point of the highest point on a UV plane model as a pole, and establishing a polar coordinate system on the UV plane model by taking a ray formed by the pole and the farthest point at the mapping point of the UV plane model as a polar axis; and realizing the positioning division of the surface area of the intracranial aneurysm according to the polar coordinates of the mapping points of the area division points of the three-dimensional aneurysm model in the UV plane model. The method can realize accurate positioning of the surface area of the intracranial aneurysm.
Description
Technical Field
The invention belongs to the field of medical image data processing, and particularly relates to a method and a device for positioning a surface region of an intracranial aneurysm, computer equipment and a storage medium.
Background
Intracranial aneurysm is the abnormal expansion of intracranial arterial blood vessels due to congenital abnormality or acquired injury. Aneurysms, when unbroken, often cause less symptoms, and after rupture, the aneurysm bleeds or bleeds, with higher mortality and often acute morbidity. Therefore, the intracranial medical image has important significance for three-dimensional reconstruction of the aneurysm and the nearby arterial blood vessel and research of the surface area information of the aneurysm according to morphology and hemodynamics.
In the prior art, a three-dimensional aneurysm vascular model is reconstructed based on an intracranial angiography image, the aneurysm model is divided by a aneurysm neck boundary space curve, isometric deformation is kept, the aneurysm is constructed into a sphere which is cut off by a aneurysm neck surface, finally, the aneurysm neck surface is fitted to form a plane which is an XY plane, and the direction perpendicular to the XY plane is a Z axis, and a three-dimensional Cartesian coordinate system is established to realize surface positioning of the aneurysm. The method has certain disadvantages, and specifically comprises the following steps: in the process of changing the volume of the aneurysm into a sphere, part of the surface characteristics of the aneurysm is lost, so that accurate positioning cannot be realized, and the method is not suitable for aneurysms with complex shapes.
Patent document CN109447967A discloses a method and system for segmenting an intracranial aneurysm image, including: intercepting a local three-dimensional image from an image of an intracranial tumor-bearing blood vessel to be segmented; acquiring the maximum inscribed circle of the tumor-bearing blood vessel image in the local three-dimensional image, and calculating the central line and the radius of the intracranial tumor-bearing blood vessel image; and performing the segmentation of the intracranial aneurysm image based on the central line and the radius of the intracranial aneurysm vessel image. The method can realize automatic segmentation of images of the intracranial aneurysm, but cannot realize positioning of information of surface areas of the intracranial aneurysm.
Patent document CN109584261A discloses a method and system for segmenting an intracranial aneurysm image, including: intercepting a local three-dimensional image from an image of an intracranial tumor-bearing blood vessel to be segmented; acquiring a tree-shaped central line of the local three-dimensional image, and calculating the central line and the radius of the intracranial tumor-carrying blood vessel image; and performing the segmentation of the intracranial aneurysm image based on the central line and the radius of the intracranial aneurysm vessel image. The method can realize automatic segmentation of images of the intracranial aneurysm, but cannot realize positioning of information of surface areas of the intracranial aneurysm.
Disclosure of Invention
In view of the above technical problems, it is an object of the present invention to provide a method, an apparatus, a computer device and a storage medium for locating a surface region of an intracranial aneurysm, so as to achieve accurate locating of the surface region of the intracranial aneurysm.
To achieve the above object, a first aspect of the embodiments of the present invention provides a method for locating a surface region of an intracranial aneurysm, including the steps of:
acquiring an intracranial angiography image, and constructing a three-dimensional aneurysm vascular model based on the intracranial angiography image;
segmenting the three-dimensional aneurysm blood vessel model along the neck of the aneurysm to obtain a three-dimensional aneurysm model, and carrying out UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model;
determining a farthest point and a highest point of the aneurysm along the blood flow direction on a curve of a tumor neck boundary space on a three-dimensional aneurysm blood vessel model, taking a mapping point of the highest point on a UV plane model as a pole, and establishing a polar coordinate system on the UV plane model by taking a ray formed by the pole and the farthest point at the mapping point of the UV plane model as a polar axis;
and positioning and dividing the surface area of the intracranial aneurysm of the UV plane model according to the polar coordinate system.
In one embodiment, the UV unfolding the three-dimensional aneurysm model comprises:
and taking a tumor neck boundary space curve of the three-dimensional aneurysm model as a free boundary, and carrying out UV expansion on the grid of the three-dimensional aneurysm model to obtain a UV plane model.
In one embodiment, determining the farthest point in the direction of blood flow on the tumor neck boundary space curve on a three-dimensional aneurysm vessel model comprises:
calculating a first straight line closest to a central line of a parent artery of the aneurysm on the three-dimensional aneurysm vascular model, wherein the direction of the first straight line represents the flowing direction of arterial blood flowing through the aneurysm; calculating the centroid of a tumor neck boundary space curve formed by segmenting the three-dimensional aneurysm model and the nearest plane from the tumor neck boundary space curve; and translating the first straight line to pass through the centroid to obtain a second straight line, projecting the second straight line and the tumor neck boundary space curve to the nearest plane to obtain a corresponding projection straight line and a projection curve, and taking the intersection point of the projection straight line and the projection curve at the far end in the blood flow direction as the farthest point.
In one embodiment, a first straight line closest to the central line of the parent artery of the aneurysm is calculated by a least square method.
In one embodiment, the tumor neck boundary space curve and the nearest plane to the tumor neck boundary space curve are calculated using a least squares method.
In one embodiment, the constructing a three-dimensional aneurysm vessel model based on intracranial angiographic images comprises:
extracting an interested blood vessel region image from an intracranial angiography image, extracting data points from the blood vessel region image according to a preset threshold value to construct a rough three-dimensional aneurysm blood vessel model, and then obtaining an accurate three-dimensional aneurysm blood vessel model by using a level set method according to the rough three-dimensional aneurysm blood vessel model and the blood vessel region image.
To achieve the above object, a second aspect of embodiments of the present invention provides a surface area localization apparatus for an intracranial aneurysm, comprising:
an acquisition module for acquiring an intracranial angiographic image;
the model construction module is used for constructing a three-dimensional aneurysm vascular model based on the intracranial angiography image; the three-dimensional aneurysm model is obtained by segmenting the three-dimensional aneurysm blood vessel model along the neck of the aneurysm;
the UV unfolding module is used for carrying out UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model;
the polar coordinate system building module is used for determining a farthest point and a highest point of the aneurysm along the blood flow direction on a tumor neck boundary space curve on the three-dimensional aneurysm blood vessel model, taking a mapping point of the highest point on the UV plane model as a pole, and building a polar coordinate system on the UV plane model by taking a ray formed by the pole and the farthest point on the mapping point of the UV plane model as a polar axis;
and the positioning division module is used for positioning and dividing the surface area of the intracranial aneurysm of the UV plane model according to the polar coordinate system.
To achieve the above object, a third aspect of the embodiments of the present invention provides a computer device, including a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the method for locating a surface region of an intracranial aneurysm according to the first aspect when executing the computer program.
To achieve the above object, a fourth aspect of the embodiments of the present invention provides a computer-readable storage medium on which a computer program is stored, which, when being processed and executed, realizes the steps of the method for locating a surface region of an intracranial aneurysm according to the first aspect.
Compared with the prior art, the method, the device, the computer equipment and the storage medium for locating the surface area of the intracranial aneurysm provided by the embodiment of the invention have the beneficial effects that at least:
the UV unfolding operation of the three-dimensional model is simple, and the UV unfolding operation is suitable for aneurysms with complex forms; then, determining a farthest point and a highest point of the aneurysm along the blood flow direction on a curve of a tumor neck boundary space on a three-dimensional aneurysm vessel model, constructing a polar coordinate system by using the farthest point and the highest point as mapping points of a UV plane model, and performing positioning division on the surface area of the intracranial aneurysm on the UV plane model by using the polar coordinate system, wherein the farthest point and the highest point of the positioning respectively consider the maximum blood flow pressure and the morphological maximum height of the aneurysm, so that the constructed polar coordinate system has wider applicability, and the surface area of the aneurysm divided by using the polar coordinate system is more accurate.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a flow chart of a method for locating a surface region of an intracranial aneurysm according to one embodiment;
FIG. 2 is a three-dimensional aneurysm model provided by an embodiment;
FIG. 3 is a UV planar model provided by an embodiment;
FIG. 4 is a schematic illustration of the determination of the farthest point in the direction of blood flow on the curve of the boundary space of the neck of a tumor and the highest point of an aneurysm according to one embodiment;
FIG. 5 is a schematic diagram of constructing a polar coordinate system on a UV planar model according to an embodiment;
fig. 6 is a schematic structural diagram of a surface region locating device for an intracranial aneurysm according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be further described in detail with reference to the accompanying drawings and examples. It should be understood that the detailed description and specific examples, while indicating the scope of the invention, are intended for purposes of illustration only and are not intended to limit the scope of the invention.
Fig. 1 is a flow chart of a method for locating a surface region of an intracranial aneurysm according to one embodiment. As shown in fig. 1, an embodiment provides a method for locating a surface region of an intracranial aneurysm, comprising the steps of:
step 1, acquiring an intracranial angiography image.
An intracranial angiography image is obtained through an angiography technology, the intracranial angiography image can reflect the actual condition of blood vessels, and a three-dimensional blood vessel model can be obtained through three-dimensional reconstruction of the intracranial angiography image so as to be used for various blood vessel researches.
And 2, segmenting the three-dimensional aneurysm blood vessel model along the neck of the aneurysm to obtain the three-dimensional aneurysm model based on the three-dimensional aneurysm blood vessel model constructed by the intracranial angiography image.
In an embodiment, the process of constructing the three-dimensional aneurysm vascular model according to the intracranial angiography image comprises the following steps:
extracting an interested blood vessel region image from an intracranial angiography image, extracting data points from the blood vessel region image according to a preset threshold value to construct a rough three-dimensional aneurysm blood vessel model, and then obtaining an accurate three-dimensional aneurysm blood vessel model by using a level set method according to the rough three-dimensional aneurysm blood vessel model and the blood vessel region image. The preset threshold value determines the roughness of the three-dimensional aneurysm vascular model, and is determined according to the situation when the three-dimensional aneurysm vascular model is applied. Of course, the manner of constructing the three-dimensional aneurysm vessel model from the intracranial angiographic image is not limited to the above method.
Since the constructed three-dimensional aneurysm blood vessel model includes not only arterial blood vessels but also aneurysms carried on arterial blood vessels, which are the objects of study, the three-dimensional aneurysm blood vessel model is obtained and then the aneurysm is divided from the three-dimensional aneurysm blood vessel model along the neck of the aneurysm (the junction between the aneurysm and the arterial blood vessel) to become an independent three-dimensional aneurysm model. During segmentation, a tumor neck boundary space curve is formed at the cut of the three-dimensional aneurysm blood vessel model and the three-dimensional aneurysm model, and the tumor neck boundary space curve can be used as a basis for UV expansion and is used for determining the farthest point on the tumor neck boundary space curve along the blood flow direction.
And 3, carrying out UV expansion on the three-dimensional aneurysm model to obtain a UV plane model.
In the embodiment, on the basis of obtaining the three-dimensional aneurysm model, the three-dimensional aneurysm model is subjected to UV expansion, specifically, a tumor neck boundary space curve of the three-dimensional aneurysm model is taken as a free boundary, and a grid of the three-dimensional aneurysm model is subjected to UV expansion to obtain a UV plane model.
For the UV unfolding, the parameter domain variables of a parametric surface can be generally represented by the coordinates UV, where U represents the horizontal direction and V represents the vertical direction, such as the parametric surface F (U, V), with the UV coordinates being similar to the projection of a three-dimensional space with two points to map points in the three-dimensional full surface. So essentially a three-dimensional parametric surface is the result of a two-dimensional parametric plane embedded in a three-dimensional space. For the parametric surface of the triangular mesh, a one-to-one mapping is established between each triangular vertex and the parametric surface, and the mapping is UV expansion.
Fig. 2 is an exemplary three-dimensional aneurysm model provided by an embodiment, when the three-dimensional aneurysm model is UV-expanded, a UV planar model as shown in fig. 3 is obtained, and when the three-dimensional aneurysm model is expanded into a two-dimensional UV planar model, a triangle mesh forming the surface of the three-dimensional aneurysm model is deformed, wherein the size of the triangle mesh in the middle area is relatively reduced, and the size of the triangle mesh in the edge area is relatively increased, but the change of the size of the triangle mesh does not affect the positioning of the surface area of the aneurysm.
And 4, determining a farthest point and a highest point of the aneurysm along the blood flow direction on a curve of a tumor neck boundary space on the three-dimensional aneurysm blood vessel model, and constructing a polar coordinate system according to mapping points of the farthest point and the highest point on the UV plane model.
In the embodiment, the polar coordinate system constructed on the UV plane model has universality, is suitable for positioning the surface region of an aneurysm in any shape, and also improves the convenience and accuracy of positioning the surface region of the aneurysm by using the polar coordinate system. Two most suitable points need to be found to determine the poles and polar axes to form a polar coordinate system. Through a large number of researches, in the morphological indexes of the aneurysm, the highest point corresponding to the maximum height of the aneurysm is a point with a relatively physiological significance, and the mapping point of the highest point corresponding to the UV plane model is used as a pole, so that the constructed polar coordinate system can very conveniently and accurately position various surface areas of the intracranial aneurysm, wherein the surface areas comprise a top area, a back flow area, a lateral area and a head flow area. When studying the characteristics of the aneurysm, the position most susceptible to impact in the blood flow direction, that is, the farthest point in the blood flow direction on the curve of the tumor neck boundary space, which has a physiologically important significance because the blood flow has the largest impact, is selected to correspond to the mapping point of the UV plane model and the pole to determine the polar axis.
In an embodiment, determining the farthest point in the direction of blood flow on the tumor neck boundary space curve on a three-dimensional aneurysm vessel model comprises:
as shown in fig. 4, on the three-dimensional aneurysm blood vessel model, a first straight line L1 closest to the central line of a parent artery of the aneurysm is calculated, and the direction of the first straight line L1 represents the flowing direction of the arterial blood flowing through the aneurysm; calculating a centroid P1 of a neck boundary space curve C1 formed by segmenting the three-dimensional aneurysm model and a nearest plane S1 to the neck boundary space curve C1; the first straight line L1 is translated to pass through the centroid P1 to obtain a second straight line L2, the second straight line L2 and the tumor neck boundary space curve C1 are projected to the nearest plane S1 to obtain a corresponding projection straight line L3 and a corresponding projection curve C2, and the intersection point of the projection straight line L3 and the distal end of the projection curve C2 in the blood flow direction is taken as a farthest point P2.
After determining the farthest point P2, the corresponding point X' of the farthest point P2 on the three-dimensional model of the aneurysm, i.e. the point on the neck of the aneurysm corresponding to the direction of blood flow, is found from the position of the farthest point P2 on the curvilinear projection C2. And simultaneously determining the point O' farthest from the nearest plane S1 on the three-dimensional model of the aneurysm, namely the highest point representing the vertical height of the aneurysm. The points X 'and O' are also labeled in the three-dimensional aneurysm vessel model.
After determining the point X 'and the point O', a polar coordinate system may be established on the UV plane model, as shown in fig. 5, the mapped point corresponding to the point O 'on the UV plane model is the point O, and the mapped point corresponding to the point X' on the UV plane model is the point X. And (3) establishing a polar coordinate system by taking the point O as a pole and the ray OX as a polar axis, so that the points on the surface of the three-dimensional aneurysm model can be represented by (rho, theta) in the polar coordinate system, rho represents the relative distance from the pole O, and theta represents the angle of the line connecting any one point with the pole.
And 5, positioning and dividing the surface area of the intracranial aneurysm of the UV plane model according to the polar coordinate system.
In the embodiment, on the basis of obtaining the polar coordinate system, the aneurysm surface is partitioned according to the polar coordinates in the polar coordinate system, so that the aneurysm surfaces with different shapes can be partitioned into different areas according to a unified standard, and meanwhile, the targeted analysis and evaluation can be performed on each partition by combining the hemodynamic parameters of the aneurysm surface.
According to the intracranial aneurysm surface region positioning method provided by the embodiment, the three-dimensional aneurysm model is subjected to UV expansion to obtain the UV plane model, the UV plane model can accurately keep the surface characteristics of the aneurysm so as to facilitate accurate positioning of the aneurysm surface region, and meanwhile, the UV expansion operation of the three-dimensional model is simple, so that the method is suitable for aneurysms with complex shapes; then, determining a farthest point and a highest point of the aneurysm along the blood flow direction on a curve of a tumor neck boundary space on a three-dimensional aneurysm vessel model, constructing a polar coordinate system by using the farthest point and the highest point as mapping points of a UV plane model, and performing positioning division on the surface area of the intracranial aneurysm on the UV plane model by using the polar coordinate system, wherein the farthest point and the highest point of the positioning respectively consider the maximum blood flow pressure and the morphological maximum height of the aneurysm, so that the constructed polar coordinate system has wider applicability, and the surface area of the aneurysm divided by using the polar coordinate system is more accurate.
Fig. 6 is a schematic structural diagram of a surface region locating device for an intracranial aneurysm according to an embodiment. As shown in fig. 6, an embodiment provides a surface area positioning device, including:
an acquisition module for acquiring an intracranial angiographic image;
the model construction module is used for constructing a three-dimensional aneurysm vascular model based on the intracranial angiography image; the three-dimensional aneurysm model is obtained by segmenting the three-dimensional aneurysm blood vessel model along the neck of the aneurysm;
the UV unfolding module is used for carrying out UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model;
and the polar coordinate system building module is used for determining a farthest point and a highest point of the aneurysm along the blood flow direction on a curve of a tumor neck boundary space on the three-dimensional aneurysm blood vessel model, and building a polar coordinate system according to the mapping points of the farthest point and the highest point on the UV plane model.
And the positioning division module is used for positioning and dividing the surface area of the intracranial aneurysm of the UV plane model according to the polar coordinate system.
It should be noted that, when the surface region positioning device for an intracranial aneurysm provided in the above embodiment is used to position the surface region of the intracranial aneurysm, the division of the functional modules is taken as an example, and the function distribution may be completed by different functional modules according to needs, that is, the internal structure of the terminal or the server is divided into different functional modules to complete all or part of the functions described above. In addition, the surface area positioning device for the intracranial aneurysm provided by the above embodiment and the surface area positioning method for the intracranial aneurysm belong to the same concept, and the specific implementation process is described in detail in the embodiment of the surface area positioning method for the intracranial aneurysm, and is not described again here.
Embodiments also provide a computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor when executing the computer program implementing the above method for locating a surface region of an intracranial aneurysm, comprising the steps of:
step 1, acquiring an intracranial angiography image.
And 2, constructing a three-dimensional aneurysm blood vessel model based on the intracranial angiography image, and segmenting the three-dimensional aneurysm blood vessel model along the neck of the aneurysm to obtain the three-dimensional aneurysm model.
And 3, carrying out UV expansion on the three-dimensional aneurysm model to obtain a UV plane model.
And 4, determining a farthest point and a highest point of the aneurysm along the blood flow direction on a curve of a tumor neck boundary space on the three-dimensional aneurysm blood vessel model, and constructing a polar coordinate system according to mapping points of the farthest point and the highest point on the UV plane model.
And 5, positioning and dividing the surface area of the intracranial aneurysm of the UV plane model according to the polar coordinate system.
In practical applications, the memory may be a volatile memory at the near end, such as RAM, a non-volatile memory, such as ROM, FLASH, a floppy disk, a mechanical hard disk, etc., or a remote storage cloud. The processor may be a Central Processing Unit (CPU), a microprocessor unit (MPU), a Digital Signal Processor (DSP), or a Field Programmable Gate Array (FPGA), i.e. the surface region localization step of the intracranial aneurysm may be achieved by these processors.
Embodiments also provide a computer-readable storage medium having stored thereon a computer program which, when being processed and executed, realizes the above method for locating a surface region of an intracranial aneurysm, comprising the steps of:
step 1, acquiring an intracranial angiography image.
And 2, constructing a three-dimensional aneurysm blood vessel model based on the intracranial angiography image, and segmenting the three-dimensional aneurysm blood vessel model along the neck of the aneurysm to obtain the three-dimensional aneurysm model.
And 3, carrying out UV expansion on the three-dimensional aneurysm model to obtain a UV plane model.
And 4, determining a farthest point and a highest point of the aneurysm along the blood flow direction on a curve of a tumor neck boundary space on the three-dimensional aneurysm blood vessel model, and constructing a polar coordinate system according to mapping points of the farthest point and the highest point on the UV plane model.
And 5, positioning and dividing the surface area of the intracranial aneurysm of the UV plane model according to the polar coordinate system.
The computer readable storage medium may be, among others, ROM, Random Access Memory (RAM), CD-ROM, magnetic tape, floppy disk, optical data storage device, and the like.
The above-mentioned embodiments are intended to illustrate the technical solutions and advantages of the present invention, and it should be understood that the above-mentioned embodiments are only the most preferred embodiments of the present invention, and are not intended to limit the present invention, and any modifications, additions, equivalents, etc. made within the scope of the principles of the present invention should be included in the scope of the present invention.
Claims (9)
1. A method of locating a surface region of an intracranial aneurysm, comprising the steps of:
acquiring an intracranial angiography image, and constructing a three-dimensional aneurysm vascular model based on the intracranial angiography image;
segmenting the three-dimensional aneurysm blood vessel model along the neck of the aneurysm to obtain a three-dimensional aneurysm model, and carrying out UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model;
determining a farthest point and a highest point of the aneurysm along the blood flow direction on a curve of a tumor neck boundary space on a three-dimensional aneurysm blood vessel model, taking a mapping point of the highest point on a UV plane model as a pole, and establishing a polar coordinate system on the UV plane model by taking a ray formed by the pole and the farthest point at the mapping point of the UV plane model as a polar axis;
and positioning and dividing the surface area of the intracranial aneurysm of the UV plane model according to the polar coordinate system.
2. The method for locating the surface area of an intracranial aneurysm as in claim 1, wherein the UV-deploying the three-dimensional aneurysm model comprises:
and taking a tumor neck boundary space curve of the three-dimensional aneurysm model as a free boundary, and carrying out UV expansion on the grid of the three-dimensional aneurysm model to obtain a UV plane model.
3. The method for locating the surface area of an intracranial aneurysm as in claim 1, wherein determining the farthest point in the direction of blood flow on the contour of the tumor neck boundary space on the three-dimensional aneurysm vessel model comprises:
calculating a first straight line closest to a central line of a parent artery of the aneurysm on the three-dimensional aneurysm vascular model, wherein the direction of the first straight line represents the flowing direction of arterial blood flowing through the aneurysm; calculating the centroid of a tumor neck boundary space curve formed by segmenting the three-dimensional aneurysm model and the nearest plane from the tumor neck boundary space curve; and translating the first straight line to pass through the centroid to obtain a second straight line, projecting the second straight line and the tumor neck boundary space curve to the nearest plane to obtain a corresponding projection straight line and a projection curve, and taking the intersection point of the projection straight line and the projection curve at the far end in the blood flow direction as the farthest point.
4. The method of claim 3, wherein the first line closest to the centerline of the parent artery of the aneurysm is calculated using least squares.
5. The method for locating the surface area of an intracranial aneurysm as in claim 3, wherein the neck boundary space curve and the closest plane to the neck boundary space curve are calculated using a least squares method.
6. The method for locating the surface region of an intracranial aneurysm as in claim 1, wherein the constructing a three-dimensional aneurysm vessel model based on intracranial angiographic images comprises:
extracting an interested blood vessel region image from an intracranial angiography image, extracting data points from the blood vessel region image according to a preset threshold value to construct a rough three-dimensional aneurysm blood vessel model, and then obtaining an accurate three-dimensional aneurysm blood vessel model by using a level set method according to the rough three-dimensional aneurysm blood vessel model and the blood vessel region image.
7. A surface area localization device for an intracranial aneurysm, comprising:
an acquisition module for acquiring an intracranial angiographic image;
the model construction module is used for constructing a three-dimensional aneurysm vascular model based on the intracranial angiography image; the three-dimensional aneurysm model is obtained by segmenting the three-dimensional aneurysm blood vessel model along the neck of the aneurysm;
the UV unfolding module is used for carrying out UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model;
the polar coordinate system building module is used for determining a farthest point and a highest point of the aneurysm along the blood flow direction on a tumor neck boundary space curve on the three-dimensional aneurysm blood vessel model, taking a mapping point of the highest point on the UV plane model as a pole, and building a polar coordinate system on the UV plane model by taking a ray formed by the pole and the farthest point on the mapping point of the UV plane model as a polar axis;
and the positioning division module is used for positioning and dividing the surface area of the intracranial aneurysm of the UV plane model according to the polar coordinate system.
8. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor when executing the computer program realizes the steps of the method for surface area localization of intracranial aneurysms according to any one of claims 1-6.
9. A computer-readable storage medium, having stored thereon a computer program, characterized in that the computer program, when being processed and executed, is adapted to carry out the steps of the method for locating the surface area of an intracranial aneurysm as defined in any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111521343.9A CN114224484B (en) | 2021-12-13 | 2021-12-13 | Surface area localization method, device, computer equipment and storage medium for intracranial aneurysm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111521343.9A CN114224484B (en) | 2021-12-13 | 2021-12-13 | Surface area localization method, device, computer equipment and storage medium for intracranial aneurysm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114224484A true CN114224484A (en) | 2022-03-25 |
| CN114224484B CN114224484B (en) | 2023-05-30 |
Family
ID=80755353
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111521343.9A Active CN114224484B (en) | 2021-12-13 | 2021-12-13 | Surface area localization method, device, computer equipment and storage medium for intracranial aneurysm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114224484B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116172697A (en) * | 2023-05-04 | 2023-05-30 | 杭州脉流科技有限公司 | Method, device and equipment for estimating length of stent implanted in blood vessel after dragging |
| CN116503436A (en) * | 2023-06-26 | 2023-07-28 | 杭州脉流科技有限公司 | Method and device for automatically dividing aneurysm neck based on control points |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040054280A1 (en) * | 2002-09-18 | 2004-03-18 | Mcmorrow Gerald J. | Three-dimensional system for abdominal aortic aneurysm evaluation |
| JP2004201873A (en) * | 2002-12-25 | 2004-07-22 | Ge Medical Systems Global Technology Co Llc | Image processing method, image processing apparatus, and program |
| JP2005345373A (en) * | 2004-06-04 | 2005-12-15 | Aloka Co Ltd | Interface detecting method, and interface detector |
| CN101283929A (en) * | 2008-06-05 | 2008-10-15 | 华北电力大学 | Rebuilding method of blood vessel three-dimensional model |
| JP2009240543A (en) * | 2008-03-31 | 2009-10-22 | Kgt Inc | Aneurysm measuring method, its apparatus and computer program |
| CN107049484A (en) * | 2017-01-21 | 2017-08-18 | 绍兴市人民医院 | Intracranial aneurysm 3D model Method of printings based on DSA data |
| CN108261217A (en) * | 2018-03-27 | 2018-07-10 | 向建平 | A kind of blood flow guider of no cripetura |
| CN109260596A (en) * | 2018-05-29 | 2019-01-25 | 杨富云 | A kind of device of micro-wave therapeutic tumour |
| CN109924956A (en) * | 2019-04-19 | 2019-06-25 | 广州新脉科技有限公司 | A kind of measurement method and device of the Morphologic Parameters of intracranial aneurysm image |
| WO2020083375A1 (en) * | 2018-10-26 | 2020-04-30 | Union Strong (Beijing) Technology Co. Ltd. | Method and System for Measuring Morphological Parameters of An Intracranial Aneurysm Image |
| CN111652881A (en) * | 2020-07-01 | 2020-09-11 | 杭州脉流科技有限公司 | Coronary artery reconstruction and fractional flow reserve calculation method, device and equipment based on deep learning and readable storage medium |
| CN112807083A (en) * | 2020-06-01 | 2021-05-18 | 上海库欣医疗科技有限公司 | Craniocerebral puncture path establishing method and system for neurosurgical navigation |
| CN113066583A (en) * | 2021-03-30 | 2021-07-02 | 昆明同心医联科技有限公司 | Aneurysm rupture risk prediction method, aneurysm rupture risk prediction device and storage medium |
| WO2021153721A1 (en) * | 2020-01-30 | 2021-08-05 | 国立大学法人京都大学 | Flow pattern word representation device, word representation method, and program |
-
2021
- 2021-12-13 CN CN202111521343.9A patent/CN114224484B/en active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040054280A1 (en) * | 2002-09-18 | 2004-03-18 | Mcmorrow Gerald J. | Three-dimensional system for abdominal aortic aneurysm evaluation |
| JP2004201873A (en) * | 2002-12-25 | 2004-07-22 | Ge Medical Systems Global Technology Co Llc | Image processing method, image processing apparatus, and program |
| JP2005345373A (en) * | 2004-06-04 | 2005-12-15 | Aloka Co Ltd | Interface detecting method, and interface detector |
| JP2009240543A (en) * | 2008-03-31 | 2009-10-22 | Kgt Inc | Aneurysm measuring method, its apparatus and computer program |
| CN101283929A (en) * | 2008-06-05 | 2008-10-15 | 华北电力大学 | Rebuilding method of blood vessel three-dimensional model |
| CN107049484A (en) * | 2017-01-21 | 2017-08-18 | 绍兴市人民医院 | Intracranial aneurysm 3D model Method of printings based on DSA data |
| CN108261217A (en) * | 2018-03-27 | 2018-07-10 | 向建平 | A kind of blood flow guider of no cripetura |
| CN109260596A (en) * | 2018-05-29 | 2019-01-25 | 杨富云 | A kind of device of micro-wave therapeutic tumour |
| WO2020083375A1 (en) * | 2018-10-26 | 2020-04-30 | Union Strong (Beijing) Technology Co. Ltd. | Method and System for Measuring Morphological Parameters of An Intracranial Aneurysm Image |
| CN109924956A (en) * | 2019-04-19 | 2019-06-25 | 广州新脉科技有限公司 | A kind of measurement method and device of the Morphologic Parameters of intracranial aneurysm image |
| WO2021153721A1 (en) * | 2020-01-30 | 2021-08-05 | 国立大学法人京都大学 | Flow pattern word representation device, word representation method, and program |
| CN112807083A (en) * | 2020-06-01 | 2021-05-18 | 上海库欣医疗科技有限公司 | Craniocerebral puncture path establishing method and system for neurosurgical navigation |
| CN111652881A (en) * | 2020-07-01 | 2020-09-11 | 杭州脉流科技有限公司 | Coronary artery reconstruction and fractional flow reserve calculation method, device and equipment based on deep learning and readable storage medium |
| CN113066583A (en) * | 2021-03-30 | 2021-07-02 | 昆明同心医联科技有限公司 | Aneurysm rupture risk prediction method, aneurysm rupture risk prediction device and storage medium |
Non-Patent Citations (2)
| Title |
|---|
| 向建平;: "浅谈数控铣加工工件外形侧面接刀痕的处理方法", 现代职业教育 * |
| 郑刚: "磁共振脑表面成像在颅内肿瘤手术定位中的应用", 中国优秀博硕士学位论文全文数据库 (硕士) 医药卫生科技辑 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116172697A (en) * | 2023-05-04 | 2023-05-30 | 杭州脉流科技有限公司 | Method, device and equipment for estimating length of stent implanted in blood vessel after dragging |
| CN116172697B (en) * | 2023-05-04 | 2023-06-27 | 杭州脉流科技有限公司 | Method, device and equipment for estimating length of stent implanted in blood vessel after dragging |
| CN116503436A (en) * | 2023-06-26 | 2023-07-28 | 杭州脉流科技有限公司 | Method and device for automatically dividing aneurysm neck based on control points |
| CN116503436B (en) * | 2023-06-26 | 2023-09-12 | 杭州脉流科技有限公司 | Method and device for automatically dividing aneurysm neck based on control points |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114224484B (en) | 2023-05-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114224484B (en) | Surface area localization method, device, computer equipment and storage medium for intracranial aneurysm | |
| CN106504331B (en) | Tooth modeling method based on three-dimensional model retrieval | |
| CN108732582B (en) | Vehicle positioning method and device | |
| US8358819B2 (en) | System and methods for image segmentation in N-dimensional space | |
| CN103247073B (en) | Three-dimensional brain blood vessel model construction method based on tree structure | |
| US20090322749A1 (en) | Systems and methods for extracting a curve-skeleton from a volumetric image of a vessel | |
| US20080033302A1 (en) | System and method for semi-automatic aortic aneurysm analysis | |
| US20090278846A1 (en) | System and method for geometric modeling of tubular structures | |
| CN111340756B (en) | Medical image lesion detection merging method, system, terminal and storage medium | |
| WO2009022283A1 (en) | Imaging method for sampling a cross-section plane in a three-dimensional (3d) image data volume | |
| CN114419047B (en) | Method, device, equipment and storage medium for determining blood vessel morphological characteristics | |
| US20120265494A1 (en) | Method of Online Building-Model Reconstruction Using Photogrammetric Mapping System | |
| CN112161622B (en) | Robot footprint planning method and device, readable storage medium and robot | |
| CN113706473A (en) | Method for determining long and short axes of lesion region in ultrasonic image and ultrasonic equipment | |
| JPH09198490A (en) | Three-dimensional discrete data projector | |
| JP2005508680A (en) | Branch selection method for probe alignment | |
| CN112712590B (en) | Method and system for generating animal point cloud | |
| CN113506365A (en) | Image display method and related device, electronic equipment and storage medium | |
| JP2003502723A (en) | 3D shape measurement using statistical curvature analysis | |
| CN116228718A (en) | Detection method, detection device, electronic equipment and storage medium | |
| CN115018825B (en) | Coronary artery dominant type classification method, classification device and storage medium | |
| CN115063395A (en) | Ultrasonic image processing method, device, equipment and medium | |
| CN116188734A (en) | AR-based house graphic mapping method | |
| CN113506313A (en) | Image processing method and related device, electronic equipment and storage medium | |
| JP2003256871A (en) | Method and device for generating shape model, program for executing the same method and recording medium having its program recorded thereon |
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 |