CN114224484B - Surface area localization method, device, computer equipment and storage medium for intracranial aneurysm - Google Patents
Surface area localization method, device, computer equipment and storage medium for intracranial aneurysm Download PDFInfo
- Publication number
- CN114224484B CN114224484B CN202111521343.9A CN202111521343A CN114224484B CN 114224484 B CN114224484 B CN 114224484B CN 202111521343 A CN202111521343 A CN 202111521343A CN 114224484 B CN114224484 B CN 114224484B
- 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.)
- Active
Links
- 201000008450 Intracranial aneurysm Diseases 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 48
- 230000004807 localization Effects 0.000 title claims description 5
- 206010002329 Aneurysm Diseases 0.000 claims abstract description 169
- 206010028980 Neoplasm Diseases 0.000 claims abstract description 37
- 238000007917 intracranial administration Methods 0.000 claims abstract description 35
- 238000002583 angiography Methods 0.000 claims abstract description 30
- 230000002792 vascular Effects 0.000 claims abstract description 28
- 230000017531 blood circulation Effects 0.000 claims abstract description 25
- 238000013507 mapping Methods 0.000 claims abstract description 22
- 210000004204 blood vessel Anatomy 0.000 claims description 36
- 238000004590 computer program Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 8
- 210000001367 artery Anatomy 0.000 claims description 5
- 230000008321 arterial blood flow Effects 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
- 238000011160 research Methods 0.000 description 3
- 230000000004 hemodynamic effect Effects 0.000 description 2
- 206010010356 Congenital anomaly Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 238000011144 upstream manufacturing Methods 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)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Image Analysis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention discloses a method, a device, a computer device and a storage medium for locating the surface area of an intracranial aneurysm, which comprise the following steps: acquiring an intracranial angiography image, and constructing a three-dimensional aneurysm vessel model based on the intracranial angiography image; dividing a three-dimensional aneurysm vascular model along a aneurysm neck to obtain a three-dimensional aneurysm model, and performing UV (ultraviolet) unfolding on the three-dimensional aneurysm model to obtain a UV plane model; determining the furthest point along the blood flow direction on a tumor neck boundary space curve and the highest point of the aneurysm on a three-dimensional aneurysm vascular model, taking the mapping point of the highest point on a UV plane model as a pole, taking a ray formed by the pole and the furthest point on the mapping point on the UV plane model as a polar axis, and establishing a polar coordinate system on the UV plane model; and positioning and dividing the surface area of the intracranial aneurysm according to the polar coordinates of the area division points of the three-dimensional aneurysm model at the mapping points of 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 surface area positioning method and device for intracranial aneurysms, computer equipment and a storage medium.
Background
Intracranial aneurysms are abnormal distensions of the intracranial arterial vessel due to congenital anomalies or acquired injuries and other factors. Aneurysms usually cause less symptoms when not ruptured, but after rupture, they bleed or bleed, are more fatal and often develop urgently. Therefore, intracranial medical imaging has great significance in three-dimensional reconstruction of aneurysms and nearby arterial blood vessels and studying surface area information of arterial aneurysms based on morphology and hemodynamics.
In the prior art, a three-dimensional aneurysm vascular model is reconstructed based on intracranial angiography images, the aneurysm model is segmented by using a tumor neck boundary space curve, the equal volume deformation is kept, the aneurysm is constructed into a sphere which is truncated by a tumor neck surface, finally, a plane is fitted by the tumor neck surface and is an XY plane, the direction perpendicular to the XY plane is a Z axis, and a three-dimensional Cartesian coordinate system is established to realize the surface positioning of the aneurysm. The method has certain disadvantages and specifically comprises the following steps: in the process of forming the aneurysm into a sphere in an equal volume, the surface characteristics of part of the aneurysm are lost, so that accurate positioning cannot be realized, and the method is not suitable for the aneurysm with complex morphology.
Patent document CN109447967a discloses a segmentation method and system of intracranial aneurysm image, comprising: intercepting a local three-dimensional image from an intracranial tumor-bearing blood vessel image 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; based on the central line and the radius of the intracranial aneurysm-carrying blood vessel image, the intracranial aneurysm image is segmented. The method can realize automatic segmentation of the intracranial aneurysm image, but can not realize positioning of intracranial aneurysm surface area information.
Patent document CN109584261a discloses a segmentation method and system of intracranial aneurysm image, comprising: intercepting a local three-dimensional image from an intracranial tumor-bearing blood vessel image 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-bearing blood vessel image; based on the central line and the radius of the intracranial aneurysm-carrying blood vessel image, the intracranial aneurysm image is segmented. The method can realize automatic segmentation of the intracranial aneurysm image, but can not realize positioning of intracranial aneurysm surface area information.
Disclosure of Invention
In view of the above technical problems, an object of the present invention is to provide a method, an apparatus, a computer device and a storage medium for locating a surface area of an intracranial aneurysm, so as to achieve accurate location of the surface area of the intracranial aneurysm.
To achieve the above object, according to a first aspect of the present invention, there is provided a method for locating a surface area of an intracranial aneurysm, comprising the steps of:
acquiring an intracranial angiography image, and constructing a three-dimensional aneurysm vessel model based on the intracranial angiography image;
dividing a three-dimensional aneurysm vascular model along a aneurysm neck to obtain a three-dimensional aneurysm model, and performing UV (ultraviolet) unfolding on the three-dimensional aneurysm model to obtain a UV plane model;
determining the furthest point along the blood flow direction on a tumor neck boundary space curve and the highest point of the aneurysm on a three-dimensional aneurysm vascular model, taking the mapping point of the highest point on a UV plane model as a pole, taking a ray formed by the pole and the furthest point on the mapping point on the UV plane model as a polar axis, and establishing a polar coordinate system on the UV plane model;
and carrying out positioning division on the surface area of the intracranial aneurysm according to the UV plane model by a polar coordinate system.
In one embodiment, the UV-unfolding of the three-dimensional aneurysm model comprises:
and taking the tumor neck boundary space curve of the three-dimensional aneurysm model as a free boundary, and performing UV expansion on the grid of the three-dimensional aneurysm model to obtain a UV plane model.
In one embodiment, determining the furthest point along the direction of blood flow on the tumor neck boundary spatial curve on the three-dimensional aneurysm vessel model comprises:
on the three-dimensional aneurysm vessel model, calculating a first straight line closest to the central line of the aneurysm-carrying artery at the proximal end of the aneurysm, wherein the direction of the first straight line represents the flowing direction of arterial blood flow flowing through the aneurysm; calculating the centroid of a tumor neck boundary space curve formed by dividing the three-dimensional aneurysm model and the nearest plane away from the tumor neck boundary space curve; and translating the first straight line 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 furthest point.
In one embodiment, a least squares method is used to calculate a first line that is closest to the centerline of the parent artery proximal to the aneurysm.
In one embodiment, the neck boundary space curve and the nearest plane to the 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 angiography images comprises:
extracting a blood vessel region image of interest from the intracranial angiography image, extracting data points from the blood vessel region image according to a preset threshold value, constructing a rough three-dimensional aneurysm blood vessel model, and 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 an embodiment of the present invention provides a surface area positioning device for intracranial aneurysms, including:
the acquisition module is used for acquiring intracranial angiography images;
the model construction module is used for constructing a three-dimensional aneurysm blood vessel model based on intracranial angiography images; the method is also used for dividing the three-dimensional aneurysm vascular model along the neck of the aneurysm to obtain a three-dimensional aneurysm model;
the UV unfolding module is used for carrying out UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model;
the system comprises a polar coordinate system construction module, a three-dimensional aneurysm vascular model, a model analysis module and a model analysis module, wherein the polar coordinate system construction module is used for determining the furthest point along the blood flow direction and the highest point of the aneurysm on a tumor neck boundary space curve on the three-dimensional aneurysm vascular model, taking the mapping point of the highest point on the UV plane model as a pole, taking a ray formed by the pole and the furthest point on the mapping point of the UV plane model as a polar axis, and constructing a polar coordinate system on the UV plane model;
and the positioning dividing module is used for carrying out positioning division on the surface area of the intracranial aneurysm on the UV plane model according to the polar coordinate system.
To achieve the above object, according to a third aspect of the embodiments of the present invention, there is provided a computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor executing the computer program to perform the steps of the method for positioning a surface area of an intracranial aneurysm according to the first aspect.
To achieve the above object, according to a fourth aspect of the embodiments of the present invention, there is provided a computer-readable storage medium having stored thereon a computer program which, when processed and executed, implements the steps of the method for positioning a surface area 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 positioning the surface area of the intracranial aneurysm have the beneficial effects that at least:
the three-dimensional aneurysm model is subjected to UV expansion to obtain a UV plane model, the surface characteristics of the aneurysm can be accurately reserved by the UV plane model, so that the surface area of the aneurysm can be accurately positioned, meanwhile, the UV expansion operation of the three-dimensional model is simple, and the method is suitable for aneurysms with complex forms; and then, determining the furthest point along the blood flow direction and the highest point of the aneurysm on a three-dimensional aneurysm vessel model, constructing a polar coordinate system by using the furthest point and the highest point at mapping points of a UV plane model, and carrying out positioning division of the surface area of the intracranial aneurysm on the UV plane model by using the polar coordinate system, wherein the furthest point and the highest point respectively consider the maximum stamping of the blood flow 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 that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for locating a surface area of an intracranial aneurysm, according to one embodiment;
FIG. 2 is a three-dimensional model of an aneurysm provided by an embodiment;
FIG. 3 is a UV planar model provided by an embodiment;
FIG. 4 is a schematic diagram of an embodiment of determining the furthest point along the direction of blood flow and the highest point of an aneurysm on a spatial curve of the neck boundary of a tumor;
FIG. 5 is a schematic diagram of constructing a polar coordinate system on a UV planar model according to one embodiment;
fig. 6 is a schematic structural view of a surface area positioning device for intracranial aneurysms, in accordance with an embodiment.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the detailed description is presented by way of example only and is not intended to limit the scope of the invention.
FIG. 1 is a flow chart of a method for locating a surface area of an intracranial aneurysm, according to one embodiment. As shown in fig. 1, the embodiment provides a method for locating a surface area of an intracranial aneurysm, which includes the following steps:
step 1, acquiring intracranial angiography images.
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 for various blood vessel researches.
And 2, dividing the three-dimensional aneurysm vascular model along the neck of the aneurysm based on the three-dimensional aneurysm vascular model constructed by the intracranial angiography image to obtain the three-dimensional aneurysm model.
In an embodiment, the process of constructing the three-dimensional aneurysm vessel model according to the intracranial angiography image is as follows:
extracting a blood vessel region image of interest from the intracranial angiography image, extracting data points from the blood vessel region image according to a preset threshold value, constructing a rough three-dimensional aneurysm blood vessel model, and 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 the roughness 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 angiography image is not limited to the above method.
The constructed three-dimensional aneurysm vascular model not only comprises an arterial blood vessel but also comprises an aneurysm carried on the arterial blood vessel, and the aneurysm is the research object, so that the aneurysm is separated from the three-dimensional aneurysm vascular model along the neck of the aneurysm (the joint of the aneurysm and the arterial blood vessel) to form an independent three-dimensional aneurysm model on the basis of obtaining the three-dimensional aneurysm vascular model. Upon segmentation, a tumor neck boundary space curve is formed at the three-dimensional aneurysm vessel model and the cut of the three-dimensional aneurysm model, which can be used as a basis for UV expansion and for determining the furthest point along the blood flow direction on the tumor neck boundary space curve.
And 3, performing UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model.
In the embodiment, on the basis of obtaining a three-dimensional aneurysm model, performing UV expansion on the three-dimensional aneurysm model, specifically, performing UV expansion on a grid of the three-dimensional aneurysm model by taking a tumor neck boundary space curve of the three-dimensional aneurysm model as a free boundary, so as to obtain a UV plane model.
For UV expansion, the parametric domain variable 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 parametric surface F (U, V), and UV coordinates are similar to the projection of a three-dimensional space, with two points to map points in a three-dimensional universe. The three-dimensional parametric surface is essentially the result of a two-dimensional parametric plane embedded in a three-dimensional space. For the parameter curved surface of the triangular mesh, each triangular vertex is built to be mapped with a parameter plane one by one, and the mapping is UV expansion.
Fig. 2 is a schematic diagram of an exemplary three-dimensional aneurysm model according to an embodiment, when the three-dimensional aneurysm model is UV-unfolded, a UV plane model as shown in fig. 3 is obtained, and when the three-dimensional aneurysm model is unfolded into a two-dimensional UV plane model, triangular meshes forming the surface of the three-dimensional aneurysm model are deformed, wherein the size of triangular meshes in a middle area is relatively reduced, the size of triangular meshes in an edge area is relatively increased, but the change of the size of triangular meshes does not affect the positioning of the surface area of the aneurysm.
And 4, determining the furthest point along the blood flow direction and the highest point of the aneurysm on the tumor neck boundary space curve on the three-dimensional aneurysm vascular model, and constructing a polar coordinate system according to the mapping points of the furthest point and the highest point on the UV plane model.
In the embodiment, in order to have universality on a polar coordinate system constructed on a UV planar model, the method is suitable for positioning the surface area of the aneurysm in any shape, and also in order to improve the convenience and accuracy of positioning the surface area of the aneurysm by using the polar coordinate system. Two most suitable points need to be found to determine the pole and polar axis to form a polar coordinate system. Through a great deal of researches, the highest point corresponding to the maximum height of the aneurysm is a point with physiological significance in the morphological index of the aneurysm, the mapping point corresponding to the highest point in the UV plane model is taken as a pole, and the constructed polar coordinate system is very convenient and accurate to position various surface areas of the intracranial aneurysm, wherein the surface areas comprise a tumor top area, a back flow area, a lateral area and an upstream area. When the characteristics of the aneurysm are studied, the position most susceptible to impact along the blood flow direction, namely the furthest point along the blood flow direction on the tumor neck boundary space curve, and the furthest point has important physiological significance due to the maximum blood flow stamping, so that the furthest point is selected to correspond to the mapping point and the pole of the UV plane model to determine the polar axis.
In an embodiment, determining the furthest point along the direction of blood flow on the tumor neck boundary spatial curve on the three-dimensional aneurysm vessel model comprises:
as shown in fig. 4, on the three-dimensional aneurysm vessel model, a first straight line L1 closest to the central line of the aneurysm-carrying artery at the proximal end of the aneurysm is calculated, and the direction of the first straight line L1 represents the flowing direction of arterial blood flow through the aneurysm; calculating a centroid P1 of a tumor neck boundary space curve C1 formed by dividing the three-dimensional aneurysm model and a nearest plane S1 away from the tumor neck boundary space curve C1; and translating the first straight line L1 through the centroid P1 to obtain a second straight line L2, projecting the second straight line L2 and the tumor neck boundary space curve C1 to the nearest plane S1 to obtain a corresponding projection straight line L3 and a projection curve C2, and taking the intersection point of the projection straight line L3 and the projection curve C2 at the far end in the blood flow direction as the farthest point P2.
After determining the furthest point P2, a corresponding point X' of the furthest point P2 on the three-dimensional model of the aneurysm, i.e. a point on the neck of the aneurysm corresponding to the direction of blood flow, is found according to the position of the furthest point P2 on the curve projection C2. While the point O' furthest from the nearest plane S1, i.e. the highest point representing the vertical height of the aneurysm, is determined on the three-dimensional model of the aneurysm. The points X 'and O' are also marked in the three-dimensional model of the aneurysm vessel.
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, where the mapping point corresponding to the point O 'on the UV plane model is the point O, and the mapping point corresponding to the point X' on the UV plane model is the point X. The point O is taken as a pole, the ray OX is taken as a polar axis, and a polar coordinate system is established, so that the points on the surface of the three-dimensional aneurysm model can be expressed by (rho, theta) in the polar coordinate system, rho represents the relative distance between the points O and the pole, and theta represents the angle between any point of OX and the pole.
And 5, carrying out positioning division on the surface area of the intracranial aneurysm according to the UV plane model by using a polar coordinate system.
In the embodiment, on the basis of obtaining a polar coordinate system, the surface of the aneurysm is partitioned according to the polar coordinates in the polar coordinate system, so that the surface of the aneurysm with different shapes can be partitioned into different areas according to a uniform standard, and meanwhile, the hemodynamic parameters of the surface of the aneurysm can be combined to analyze and evaluate each partition in a targeted manner.
According to the intracranial aneurysm surface area positioning method provided by the embodiment, the three-dimensional aneurysm model is subjected to UV expansion to obtain the UV planar model, the UV planar model can accurately retain the surface characteristics of the aneurysm so as to facilitate the accurate positioning of the aneurysm surface area, and meanwhile, the UV expansion operation of the three-dimensional model is simple and is suitable for aneurysms with complex forms; and then, determining the furthest point along the blood flow direction and the highest point of the aneurysm on a three-dimensional aneurysm vessel model, constructing a polar coordinate system by using the furthest point and the highest point at mapping points of a UV plane model, and carrying out positioning division of the surface area of the intracranial aneurysm on the UV plane model by using the polar coordinate system, wherein the furthest point and the highest point respectively consider the maximum stamping of the blood flow 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 view of a surface area positioning device for intracranial aneurysms, in accordance with an embodiment. As shown in fig. 6, the surface area positioning device provided in the embodiment includes:
the acquisition module is used for acquiring intracranial angiography images;
the model construction module is used for constructing a three-dimensional aneurysm blood vessel model based on intracranial angiography images; the method is also used for dividing the three-dimensional aneurysm vascular model along the neck of the aneurysm to obtain a three-dimensional aneurysm model;
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 construction module is used for determining the furthest point along the blood flow direction and the highest point of the aneurysm on the tumor neck boundary space curve on the three-dimensional aneurysm blood vessel model, and constructing a polar coordinate system according to the mapping points of the furthest point and the highest point on the UV plane model.
And the positioning dividing module is used for carrying out positioning division on the surface area of the intracranial aneurysm on the UV plane model according to the polar coordinate system.
It should be noted that, when the surface area positioning device for an intracranial aneurysm provided in the foregoing embodiment performs surface area positioning of the intracranial aneurysm, the division of the foregoing functional modules should be exemplified, and the foregoing functional allocation may be performed by different functional modules according to needs, that is, the internal structure of the terminal or the server is divided into different functional modules, so as to perform all or part of the functions described above. In addition, the surface area positioning device for intracranial aneurysm provided in the above embodiment and the surface area positioning method for intracranial aneurysm provided in the above embodiment belong to the same concept, and the specific implementation process of the surface area positioning device for intracranial aneurysm is detailed in the surface area positioning method for intracranial aneurysm, which is not described herein.
The embodiment also provides a computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the method for locating the surface area of the intracranial aneurysm when executing the computer program, namely the method comprises the following steps:
step 1, acquiring intracranial angiography images.
And 2, after constructing a three-dimensional aneurysm vascular model based on the intracranial angiography image, dividing the three-dimensional aneurysm vascular model along the neck of the aneurysm to obtain the three-dimensional aneurysm model.
And 3, performing UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model.
And 4, determining the furthest point along the blood flow direction and the highest point of the aneurysm on the tumor neck boundary space curve on the three-dimensional aneurysm vascular model, and constructing a polar coordinate system according to the mapping points of the furthest point and the highest point on the UV plane model.
And 5, carrying out positioning division on the surface area of the intracranial aneurysm according to the UV plane model by using a polar coordinate system.
In practical applications, the memory may be a volatile memory at the near end, such as a RAM, or a nonvolatile memory, such as a ROM, a FLASH, a floppy disk, a mechanical hard disk, or a remote storage cloud. The processor may be a Central Processing Unit (CPU), a Microprocessor (MPU), a Digital Signal Processor (DSP), or a Field Programmable Gate Array (FPGA), i.e. the surface area localization step of the intracranial aneurysm may be implemented by these processors.
The embodiment also provides a computer readable storage medium, on which a computer program is stored, the computer program being processed and executed to implement the above-mentioned method for locating the surface area of an intracranial aneurysm, namely comprising the following steps:
step 1, acquiring intracranial angiography images.
And 2, after constructing a three-dimensional aneurysm vascular model based on the intracranial angiography image, dividing the three-dimensional aneurysm vascular model along the neck of the aneurysm to obtain the three-dimensional aneurysm model.
And 3, performing UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model.
And 4, determining the furthest point along the blood flow direction and the highest point of the aneurysm on the tumor neck boundary space curve on the three-dimensional aneurysm vascular model, and constructing a polar coordinate system according to the mapping points of the furthest point and the highest point on the UV plane model.
And 5, carrying out positioning division on the surface area of the intracranial aneurysm according to the UV plane model by using a polar coordinate system.
The computer readable storage medium may be, among other things, ROM, random Access Memory (RAM), CD-ROM, magnetic tape, floppy disk, optical data storage device, etc.
The foregoing detailed description of the preferred embodiments and advantages of the invention will be appreciated that the foregoing description is merely illustrative of the presently preferred embodiments of the invention, and that no changes, additions, substitutions and equivalents of those embodiments are intended to be included within the scope of the invention.
Claims (9)
1. A method of locating a surface area of an intracranial aneurysm, comprising the steps of:
acquiring an intracranial angiography image, and constructing a three-dimensional aneurysm vessel model based on the intracranial angiography image;
dividing a three-dimensional aneurysm vascular model along a aneurysm neck to obtain a three-dimensional aneurysm model, and performing UV (ultraviolet) unfolding on the three-dimensional aneurysm model to obtain a UV plane model;
determining the furthest point along the blood flow direction and the highest point of the aneurysm on a three-dimensional aneurysm vascular model on a tumor neck boundary space curve, taking a mapping point of the highest point on a UV plane model as a pole, taking a ray formed by the pole and the furthest point on the mapping point of the UV plane model as a polar axis, and establishing a polar coordinate system on the UV plane model;
and carrying out positioning division on the surface area of the intracranial aneurysm according to the UV plane model by a polar coordinate system.
2. The method of claim 1, wherein UV-expanding the three-dimensional aneurysm model comprises:
and taking the tumor neck boundary space curve of the three-dimensional aneurysm model as a free boundary, and performing UV expansion on the grid of the three-dimensional aneurysm model to obtain a UV plane model.
3. The method of claim 1, wherein determining the furthest point along the direction of blood flow on the tumor neck boundary spatial curve on the three-dimensional aneurysm vessel model comprises:
on the three-dimensional aneurysm vessel model, calculating a first straight line closest to the central line of the aneurysm-carrying artery at the proximal end of the aneurysm, wherein the direction of the first straight line represents the flowing direction of arterial blood flow flowing through the aneurysm; calculating the centroid of a tumor neck boundary space curve formed by dividing the three-dimensional aneurysm model and the nearest plane away from the tumor neck boundary space curve; and translating the first straight line 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 furthest point.
4. A method of locating a surface area of an intracranial aneurysm as claimed in claim 3, wherein the first line is calculated by least squares to be closest to the centre line of the parent artery at the proximal end of the aneurysm.
5. A method of locating a surface area of an intracranial aneurysm as claimed in claim 3, wherein the neck boundary space curve and the nearest plane to the neck boundary space curve are calculated using a least squares method.
6. The method of claim 1, wherein constructing a three-dimensional aneurysm vessel model based on intracranial angiography images comprises:
extracting a blood vessel region image of interest from the intracranial angiography image, extracting data points from the blood vessel region image according to a preset threshold value, constructing a rough three-dimensional aneurysm blood vessel model, and 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 locating device for an intracranial aneurysm, comprising:
the acquisition module is used for acquiring intracranial angiography images;
the model construction module is used for constructing a three-dimensional aneurysm blood vessel model based on intracranial angiography images; the method is also used for dividing the three-dimensional aneurysm vascular model along the neck of the aneurysm to obtain a three-dimensional aneurysm model;
the UV unfolding module is used for carrying out UV unfolding on the three-dimensional aneurysm model to obtain a UV plane model;
the system comprises a polar coordinate system construction module, a three-dimensional aneurysm vascular model, a model analysis module and a model analysis module, wherein the polar coordinate system construction module is used for determining the furthest point along the blood flow direction and the highest point of the aneurysm on a tumor neck boundary space curve on the three-dimensional aneurysm vascular model, taking the mapping point of the highest point on the UV plane model as a pole, taking a ray formed by the pole and the furthest point on the mapping point of the UV plane model as a polar axis, and constructing a polar coordinate system on the UV plane model;
and the positioning dividing module is used for carrying out positioning division on the surface area of the intracranial aneurysm on 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, carries out the steps of the method for surface area localization of intracranial aneurysms as claimed in any of claims 1-6.
9. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being processed and executed, implements the steps of the method for localization of surface areas of intracranial aneurysms as claimed in any of claims 1-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 CN114224484A (en) | 2022-03-25 |
| CN114224484B true 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) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116172697B (en) * | 2023-05-04 | 2023-06-27 | 杭州脉流科技有限公司 | Method, device and equipment for estimating length of stent implanted in blood vessel after dragging |
| CN116503436B (en) * | 2023-06-26 | 2023-09-12 | 杭州脉流科技有限公司 | Method and device for automatically dividing aneurysm neck based on control points |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6905468B2 (en) * | 2002-09-18 | 2005-06-14 | Diagnostic Ultrasound Corporation | Three-dimensional system for abdominal aortic aneurysm evaluation |
-
2021
- 2021-12-13 CN CN202111521343.9A patent/CN114224484B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
|---|
| 浅谈数控铣加工工件外形侧面接刀痕的处理方法;向建平;;现代职业教育(第20期);98-99 * |
| 磁共振脑表面成像在颅内肿瘤手术定位中的应用;郑刚;中国优秀博硕士学位论文全文数据库 (硕士) 医药卫生科技辑(第3期);25-35 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114224484A (en) | 2022-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114224484B (en) | Surface area localization method, device, computer equipment and storage medium for intracranial aneurysm | |
| JP7112993B2 (en) | Laser Radar Internal Parameter Accuracy Verification Method and Its Apparatus, Equipment and Medium | |
| KR20190100005A (en) | Method and apparatus for generating blood vessel model | |
| US20160314587A1 (en) | Processing apparatus, processing method, and non-transitory computer-readable storage medium | |
| CN111340756B (en) | Medical image lesion detection merging method, system, terminal and storage medium | |
| US8218844B2 (en) | Method for generating analytical tool for analyzing abnormalities in human organs | |
| US20080317308A1 (en) | System and methods for image segmentation in N-dimensional space | |
| US20120207366A1 (en) | Method and system for segmenting a liver object in an image | |
| US20140044333A1 (en) | System and method for providing registration between breast shapes before and during surgery | |
| AU2020322893B2 (en) | Data processing method, apparatus, device, and storage medium | |
| CN107507212B (en) | Digital brain visualization method and device, computing equipment and storage medium | |
| CN110021053A (en) | A kind of image position method, device, storage medium and equipment based on coordinate conversion | |
| CN113706473A (en) | Method for determining long and short axes of lesion region in ultrasonic image and ultrasonic equipment | |
| US11481867B2 (en) | Device and method for registering three-dimensional data | |
| JPH09198490A (en) | Three-dimensional discrete data projector | |
| CN113506365A (en) | Image display method and related device, electronic equipment and storage medium | |
| CN113516758A (en) | Image display method and related device, electronic equipment and storage medium | |
| CN115063395A (en) | Ultrasonic image processing method, device, equipment and medium | |
| US20140032180A1 (en) | Method and apparatus for computing deformation of an object | |
| CN107945157B (en) | Method/system for measuring end plate on sacrum, computer readable storage medium and equipment | |
| CN114299096A (en) | Outline delineation method, device, equipment and storage medium | |
| CN113593036A (en) | Imaging system, method, device and program with reference to grid parameterized model | |
| CN112734830A (en) | Point cloud data-based pose control method, device, equipment and storage medium | |
| CN116681705B (en) | Surface morphology measurement method and processing equipment based on longitudinal structure of human brain hippocampus | |
| CN113408595B (en) | Pathological image processing method and device, electronic equipment and readable storage medium |
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 |