WO2021098768A1 - Method and system for assessing aneurysm rupture risk - Google Patents

Method and system for assessing aneurysm rupture risk Download PDF

Info

Publication number
WO2021098768A1
WO2021098768A1 PCT/CN2020/130058 CN2020130058W WO2021098768A1 WO 2021098768 A1 WO2021098768 A1 WO 2021098768A1 CN 2020130058 W CN2020130058 W CN 2020130058W WO 2021098768 A1 WO2021098768 A1 WO 2021098768A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
aneurysm
labeled
triangle mesh
triangle
Prior art date
Application number
PCT/CN2020/130058
Other languages
French (fr)
Chinese (zh)
Inventor
马泽
宋凌
印胤
杨光明
秦岚
Original Assignee
强联智创(北京)科技有限公司
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 强联智创(北京)科技有限公司 filed Critical 强联智创(北京)科技有限公司
Publication of WO2021098768A1 publication Critical patent/WO2021098768A1/en

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Definitions

  • This specification relates to the field of medical imaging and computer technology, and in particular to a method and system for assessing the risk of aneurysm rupture.
  • Intracranial aneurysms are mostly abnormal bulges that occur on the walls of intracranial arteries. According to reports, the prevalence of unruptured intracranial aneurysms in Chinese adults can be as high as 7%. Subarachnoid hemorrhage can occur after rupture. Cause severe disability or death.
  • the treatment of intracranial aneurysm after rupture is mainly to remove intracranial hematoma and prevent blood from continuing to flow into the skull; intracranial unruptured aneurysms need to develop a personalized management plan based on the results of the aneurysm rupture risk assessment and the condition of the observed person, and be conservative Observation or surgical intervention. Therefore, the risk assessment of aneurysm rupture is of great significance.
  • the observer observes the aneurysm morphology based on the three-dimensional DSA (Digital subtraction angiography) image, MRA (Magnetic Resonance Angiography) or CTA (CT angiography, CT angiography) image, and Combining the condition of the observer to assess the risk of aneurysm rupture, this method often relies on the experience of the observer, is affected by subjective judgment, lacks objective support, and takes a long time.
  • DSA Digital subtraction angiography
  • MRA Magnetic Resonance Angiography
  • CTA CT angiography, CT angiography
  • a new method is needed that can eliminate or reduce the participation of human factors, shorten the time-consuming, realize simple and quick assessment of the risk of aneurysm rupture, and provide objective support for the assessment of the risk of aneurysm rupture.
  • the embodiments of this specification provide an aneurysm rupture risk assessment method and system, which are used to solve the following technical problems: in the prior art, the observer uses three-dimensional DSA (Digital subtraction angiography, digital subtraction angiography) angiography images or MRA (Magnetic Subtraction Angiography) images or MRA (Magnetic Subtraction Angiography) images according to the prior art.
  • DSA Digital subtraction angiography, digital subtraction angiography
  • MRA Magnetic Subtraction Angiography
  • MRA Magnetic Subtraction Angiography
  • the embodiment of this specification provides a method for assessing the risk of aneurysm rupture, which includes the following steps:
  • Acquiring data to be processed where the data to be processed includes first data and/or second data and/or third data;
  • the aneurysm risk assessment result includes the aneurysm rupture risk probability of the to-be-processed data and/or the vertex classification of the 3D triangular mesh data corresponding to the to-be-processed data, and the value of the 3D triangular mesh data corresponding to the to-be-processed data
  • the vertex is classified as whether the vertex of the 3D triangle mesh data corresponding to the to-be-processed data belongs to a blood vessel or an aneurysm;
  • the labeled 3D triangle mesh data is based on the known 3D triangle mesh data, the vertices belonging to the blood vessel in the known 3D triangle mesh data are marked as 0, and the known 3D triangle mesh data is marked as 0. The vertex belonging to the aneurysm in the 3D triangle mesh data is marked as 1.
  • labeled 3D triangle mesh data further includes:
  • the vertices belonging to the blood vessel and/or the vertices belonging to the aneurysm are randomly discarded for the labeled 3D triangle mesh data according to the preset ratio to obtain the coordinates of the vertices belonging to the aneurysm point and the points belonging to the aneurysm. Edge information connected between vertices.
  • the training of the aneurysm risk assessment model includes:
  • a one-dimensional vector of the labeled 3D triangle mesh data is obtained through an arithmetic block operation, where the one-dimensional vector represents the labeled 3D triangle mesh
  • the global feature vector of the data
  • the aneurysm risk assessment of the labeled 3D triangle grid data is obtained through a sigmoid function operation As a result, the aneurysm risk assessment model is obtained.
  • the obtaining a one-dimensional vector of the marked 3D triangular mesh data through an arithmetic block operation based on the marked 3D triangular mesh data specifically includes:
  • the labeled 3D triangle mesh data is subjected to the operations of the first operation block of (64, 64) dimensions and the second operation block of (64, 128, 1024) dimensions to obtain a matrix of (N, 1024), where N is The number of vertices in the labeled 3D triangle mesh data;
  • the second data and/or third data corresponding to the one-dimensional vector and the labeled 3D triangle mesh data are operated by a sigmoid function to obtain the labeled 3D triangle mesh data
  • the results of the aneurysm risk assessment include:
  • the matrix of (N, 2) is subjected to a sigmoid function operation to obtain the vertex classification of the labeled 3D triangle mesh data.
  • An input module to obtain data to be processed, where the data to be processed includes first data and/or second data and/or third data;
  • the evaluation module converts the first data into 3D triangle mesh data, and inputs the aneurysm risk assessment model based on the 3D triangle mesh data and/or the second data and/or the third data to obtain all The aneurysm risk assessment result of the to-be-processed data, wherein the aneurysm risk assessment model is based on the marked 3D triangle grid data and/or the second data and/or the number obtained by training on the third data Model, the aneurysm risk assessment result includes the aneurysm rupture risk probability of the data to be processed and/or the vertex classification of the 3D triangle mesh data corresponding to the data to be processed, and the 3D triangle mesh corresponding to the data to be processed The vertex of the grid data is classified as whether the vertex of the 3D triangle mesh data corresponding to the data to be processed belongs to a blood vessel or an aneurysm;
  • the output module outputs the aneurysm risk assessment result of the to-be-processed data.
  • the labeled 3D triangle mesh data is based on the known 3D triangle mesh data, the vertices belonging to the blood vessel in the known 3D triangle mesh data are marked as 0, and the known 3D triangle mesh data is marked as 0. The vertex belonging to the aneurysm in the 3D triangle mesh data is marked as 1.
  • labeled 3D triangle mesh data further includes:
  • the vertices belonging to the blood vessel and/or the vertices belonging to the aneurysm are randomly discarded for the labeled 3D triangle mesh data according to the preset ratio to obtain the coordinates of the vertices belonging to the aneurysm point and the points belonging to the aneurysm. Edge information connected between vertices.
  • the training of the aneurysm risk assessment model includes:
  • a one-dimensional vector of the labeled 3D triangle mesh data is obtained through an arithmetic block operation, where the one-dimensional vector represents the labeled 3D triangle mesh
  • the global feature vector of the data
  • the aneurysm risk assessment of the labeled 3D triangle grid data is obtained through a sigmoid function operation As a result, the aneurysm risk assessment model is obtained.
  • the obtaining a one-dimensional vector of the marked 3D triangular mesh data through an arithmetic block operation based on the marked 3D triangular mesh data specifically includes:
  • the labeled 3D triangle mesh data is subjected to the operations of the first operation block of (64, 64) dimensions and the second operation block of (64, 128, 1024) dimensions to obtain a matrix of (N, 1024), where N is The number of vertices in the labeled 3D triangle mesh data;
  • the second data and/or third data corresponding to the one-dimensional vector and the labeled 3D triangle mesh data are operated by a sigmoid function to obtain the labeled 3D triangle mesh data
  • the results of the aneurysm risk assessment include:
  • the matrix of (N, 2) is subjected to a sigmoid function operation to obtain the vertex classification of the labeled 3D triangle mesh data.
  • the embodiment of this specification acquires data to be processed, wherein the data to be processed includes first data and/or second data and/or third data; the first data is converted into 3D triangle mesh data based on the The 3D triangle grid data and/or the second data and/or the third data are input into an aneurysm risk assessment model to obtain an aneurysm risk assessment result of the to-be-processed data, wherein the aneurysm risk assessment model It is a digital model obtained by training based on the labeled 3D triangle grid data and/or the second data and/or the third data, and the aneurysm risk assessment result includes the aneurysm rupture risk of the data to be processed Probability and/or the vertex classification of the 3D triangle mesh data corresponding to the data to be processed, the vertex of the 3D triangle mesh data corresponding to the data to be processed is classified as the vertex of the 3D triangle mesh data corresponding to the data to be processed It belongs to a blood vessel or an aneurysm; outputting the aneur
  • Fig. 1 is a schematic flow chart of a method for assessing the risk of aneurysm rupture according to an embodiment of this specification
  • Figure 2 is a schematic diagram of the construction process of the aneurysm risk assessment model provided by the embodiments of this specification;
  • Fig. 3 is a schematic diagram of an aneurysm rupture risk assessment system provided by an embodiment of this specification.
  • the morphological parameters of intracranial aneurysms are of great significance for the diagnosis of intracranial aneurysms.
  • the morphological parameters of intracranial aneurysms can be obtained through medical imaging, including but not limited to DSA, MRA, and CTA.
  • DSA digitally input two frames of X-ray images taken before and after the injection of the contrast agent into the image computer. Through the process of subtraction, enhancement and re-imaging, a clear pure blood vessel image is obtained, and the vascular shadow is displayed in real time.
  • DSA has the advantages of high contrast resolution, short examination time, low contrast agent consumption, low concentration, significantly reduced patient X-ray absorption, and film saving. It is of great significance in the clinical diagnosis of vascular diseases. DSA has become the gold standard for the diagnosis of intracranial arterial malformations and aneurysms due to its imaging characteristics.
  • MRA The basic principle of MRA is based on saturation effect, inflow enhancement effect, and flow dephasing effect.
  • the pre-saturation zone is placed at the head end of the 3D slab to saturate the venous blood flow.
  • the arterial blood flowing in the reverse direction enters the 3D slab and generates MR signals because it is not saturated.
  • a thicker volume is divided into multiple thin layers for excitation, reducing the thickness of the excitation volume to reduce the inflow saturation effect, and can ensure the scanning volume range, obtain several layers of adjacent layers of thin-layer images, make the image clear, and the blood vessels are subtle
  • the structure is displayed well and the spatial resolution is improved. Due to its high-quality imaging characteristics, MRA is gradually used in the diagnosis of intracranial aneurysms.
  • the basic principle of CTA is that after intravenous injection of iodine-containing contrast agent, spiral CT or electron beam CT is used to continuously perform thin-slice scanning at the peak of the blood vessel filled with the contrast agent to quickly obtain a large number of thin-layer superimposed sections, which are reconstructed after computer image processing Stereoscopic images of blood vessels, clearly showing the anterior cerebral artery, middle artery, posterior artery and its main branches, Wi11is arterial ring, etc.
  • spiral CT or electron beam CT is used to continuously perform thin-layer scanning at the peak of the blood vessel filled with the contrast agent to quickly obtain a large number of thin-layer superimposed sections, and reconstruct the three-dimensional image of the blood vessel after computer image processing.
  • CTA has the characteristics of non-invasive, fast, simple operation and low price, it can mostly replace DSA in the diagnosis of clinical intracranial aneurysms.
  • the original image based on DSA and/or CTA and/or MRA can be digitally processed into 3D triangular grid data.
  • the estimated value can be realized Assessment of the risk of aneurysm rupture.
  • Fig. 1 is a schematic flow chart of a method for assessing the risk of aneurysm rupture according to an embodiment of this specification. The method specifically includes the following steps:
  • Step S101 Obtain data to be processed, where the data to be processed includes first data and/or second data and/or third data.
  • the first data and/or the second data and/or the third data are the data corresponding to each aneurysm sample
  • the first data is 3D image data, including but not limited to DSA image data, CTA image data, and MRA image data.
  • the second data includes but is not limited to: the age and/or gender of the observed person and/or whether it is a patient with multiple aneurysms and/or history of drinking and/or smoking and/or family history and/or arteries The location and/or symptoms of the tumor in the vascular segment.
  • the third data is the morphological parameters of the aneurysm of the observed person, including but not limited to: aneurysm volume, average diameter of tumor-bearing blood vessels, ratio of tumor body length to tumor neck diameter (SR), tumor body length and Ratio of aneurysm neck width (AR), aneurysm length, aneurysm height, aneurysm width, aneurysm neck width, inflow angle, aneurysm neck area, aneurysm neck diameter, tumor-carrying vessel length, fluctuation index (UI), non- Spherical Index (NSI).
  • aneurysm volume average diameter of tumor-bearing blood vessels
  • SR tumor body length
  • AR Ratio of aneurysm neck width
  • aneurysm length aneurysm height
  • aneurysm width aneurysm neck width
  • inflow angle aneurysm neck area
  • aneurysm neck diameter tumor-carrying vessel length
  • UI fluctuation index
  • the second data and the third data are stored in the form of a table for subsequent use.
  • the morphological parameters of the aneurysm of the observed person are obtained.
  • the method of threshold segmentation and region growth is adopted to convert the first data into a binary image;
  • the Marching Cube algorithm is used to construct 3D triangular mesh data;
  • the observed aneurysm is obtained Morphological parameters.
  • the 3D triangle mesh data is a topological structure composed of multiple triangles in a Cartesian coordinate system in space, where each triangle contains vertices and edges.
  • the 3D triangular mesh data adopts STL format data, specifically, it may be a binary text or a text file.
  • the data to be processed includes data corresponding to one or more aneurysms.
  • the observed person is the subject of the risk assessment of aneurysm rupture in this application.
  • Step S103 Convert the first data into 3D triangular mesh data, and input the aneurysm risk assessment model based on the 3D triangular mesh data and/or the second data and/or the third data to obtain all The aneurysm risk assessment result of the to-be-processed data, wherein the aneurysm risk assessment model is based on the marked 3D triangle grid data and/or the second data and/or the number obtained by training on the third data Model, the aneurysm risk assessment result includes the aneurysm rupture risk probability of the data to be processed and/or the vertex classification of the 3D triangle mesh data corresponding to the data to be processed, and the 3D triangle mesh corresponding to the data to be processed The vertex of the grid data is classified as whether the vertex of the 3D triangle mesh data corresponding to the data to be processed belongs to a blood vessel or an aneurysm.
  • the 3D triangle mesh data corresponding to the data to be processed is the method of threshold segmentation and region growth for the data to be processed, and the first data is converted into a binary image; the 3D triangle mesh data obtained by the Marching Cube algorithm .
  • FIG. 2 is a schematic diagram of the acquisition process of an aneurysm risk assessment model provided in an embodiment of this specification, which specifically includes:
  • Step S201 Annotate the 3D triangle mesh data to obtain the annotated 3D triangle mesh data.
  • the known 3D triangle mesh data Based on the known 3D triangle mesh data, the vertices belonging to the blood vessel in the known 3D triangle mesh data are marked as 0, and the vertices belonging to the aneurysm in the known 3D triangular mesh data are marked as 1 .
  • the known 3D triangle mesh data is based on the corresponding 3D triangle mesh data obtained as the first data of the training set
  • the labeled 3D triangle mesh data is further: randomly discarded the vertices belonging to the blood vessel and/or the aneurysm according to the preset ratio of the labeled 3D triangle mesh data, so as to obtain the vertices belonging to the aneurysm.
  • the preset ratio is comprehensively determined according to the amount of 3D triangle mesh data that has been marked and/or the number of vertex coordinates.
  • the labeled 3D triangle mesh data is further: randomly discarded vertices belonging to blood vessels and/or vertices belonging to aneurysms for the labeled 3D triangle mesh data according to a preset ratio.
  • the specific number of vertices 1024 or 2048 is only an exemplary description of this application, and does not constitute a limitation of this application.
  • the known 3D triangle mesh data can be input into the corresponding labeling software, so that the mesh data can be visualized, so that manual labeling can be performed point by point.
  • multi-person labeling can be used, and the number of labeling personnel shall be at least 3 people. The person making the labeling should have professional experience to ensure the accuracy of the labeling.
  • Step S203 Based on the labeled 3D triangular mesh data, a one-dimensional vector of the labeled 3D triangular mesh data is obtained through an arithmetic block operation, where the one-dimensional vector represents the labeled 3D The global feature vector of the triangle mesh data.
  • the labeled 3D triangle mesh data is subjected to the operations of the first operation block of (64, 64) dimensions and the second operation block of (64, 128, 1024) dimensions to obtain (N, 1024) ), where N is the number of vertices in the labeled 3D triangle mesh data;
  • Step S205 Based on the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle mesh data, perform a sigmoid function operation to obtain the artery of the labeled 3D triangle mesh data An aneurysm risk assessment result to obtain the aneurysm risk assessment model.
  • the second data and/or third data corresponding to the one-dimensional vector and the labeled 3D triangular mesh data are input into a multi-layer perceptron with a dimension of (512, 256, 2) to obtain Fourth data
  • the matrix of (N, 2) is subjected to a sigmoid function operation to obtain the vertex classification of the labeled 3D triangle mesh data.
  • the sigmoid function is also called the Logistic function.
  • the sigmoid classification function after the sigmoid classification function is processed, it can be ensured that the output data is within (0,1), and the risk of aneurysm rupture can be assessed.
  • the sigmoid function can also realize data classification to realize the vertex classification of the labeled 3D triangle mesh data.
  • the data used for the establishment of the aneurysm risk assessment model includes: training set and test set.
  • the training set is used to train the aneurysm risk assessment model
  • the test set is used to test the effect of the obtained aneurysm risk assessment model.
  • the aneurysm risk assessment model provided by the embodiments of this specification can automatically learn the local features of each vertex and surrounding vertices in the labeled 3D triangle mesh data, and further integrate the second data and/or third data of the observed person , So as to achieve a comprehensive and efficient prediction of aneurysm rupture risk and vertex classification.
  • Step S105 Output the aneurysm risk assessment result of the to-be-processed data.
  • the output of the aneurysm risk assessment result includes: the aneurysm rupture risk probability of the to-be-processed data and/or the vertex classification of the 3D triangle mesh corresponding to the to-be-processed data.
  • the probability of obtaining the aneurysm rupture risk of the observed person is 70%, which can be used as an auxiliary reference for subsequent diagnosis and treatment to determine Whether to perform surgical treatment.
  • a set of aneurysm data is data corresponding to an aneurysm, including the first data and/or the second data and/or the third data.
  • the method provided in the embodiments of this specification can eliminate or reduce the participation of human factors, shorten time-consuming, realize simple and quick assessment of the risk of aneurysm rupture, and provide objective support for the assessment of the risk of aneurysm rupture.
  • the evaluation method provided in the embodiments of this specification can be packaged as software in actual application to assist observers in making decisions about treatment of unruptured aneurysms to quickly obtain reasonable and reliable prediction results and/or vertex classification.
  • FIG. 3 is a schematic diagram of an aneurysm rupture risk assessment system provided by the embodiment of this specification, and the system includes:
  • the input module 301 obtains data to be processed, where the data to be processed includes first data and/or second data and/or third data;
  • the evaluation module 303 converts the first data into 3D triangle mesh data, and inputs the aneurysm risk assessment model based on the 3D triangle mesh data and/or the second data and/or the third data to obtain The aneurysm risk assessment result of the to-be-processed data, wherein the aneurysm risk assessment model is obtained by training based on the marked 3D triangle grid data and/or the second data and/or the third data A digital model, where the aneurysm risk assessment result includes the aneurysm rupture risk probability of the to-be-processed data and/or the vertex classification of the 3D triangle mesh data corresponding to the to-be-processed data, and the 3D triangle corresponding to the to-be-processed data The vertex of the mesh data is classified as whether the vertex of the 3D triangle mesh data corresponding to the data to be processed belongs to a blood vessel or an aneurysm;
  • the output module 305 outputs the aneurysm risk assessment result of the to-be-processed data.
  • the labeled 3D triangle mesh data is based on the known 3D triangle mesh data, the vertices belonging to the blood vessel in the known 3D triangle mesh data are marked as 0, and the known 3D triangle mesh data is marked as 0. The vertex belonging to the aneurysm in the 3D triangle mesh data is marked as 1.
  • labeled 3D triangle mesh data further includes:
  • the vertices belonging to the blood vessel and/or the vertices belonging to the aneurysm are randomly discarded for the labeled 3D triangle mesh data according to the preset ratio to obtain the coordinates of the vertices belonging to the aneurysm point and the points belonging to the aneurysm. Edge information connected between vertices.
  • the training of the aneurysm risk assessment model includes:
  • a one-dimensional vector of the labeled 3D triangle mesh data is obtained through an arithmetic block operation, where the one-dimensional vector represents the labeled 3D triangle mesh
  • the global feature vector of the data
  • the aneurysm risk assessment of the labeled 3D triangle grid data is obtained through a sigmoid function operation As a result, the aneurysm risk assessment model is obtained.
  • the obtaining a one-dimensional vector of the marked 3D triangular mesh data through an arithmetic block operation based on the marked 3D triangular mesh data specifically includes:
  • the labeled 3D triangle mesh data is subjected to the operations of the first operation block of (64, 64) dimensions and the second operation block of (64, 128, 1024) dimensions to obtain a matrix of (N, 1024), where N is The number of vertices in the labeled 3D triangle mesh data;
  • the second data and/or third data corresponding to the one-dimensional vector and the labeled 3D triangle mesh data are operated by a sigmoid function to obtain the labeled 3D triangle mesh data
  • the results of the aneurysm risk assessment include:
  • the matrix of (N, 2) is subjected to a sigmoid function operation to obtain the vertex classification of the labeled 3D triangular mesh data.
  • the device, electronic device, non-volatile computer storage medium and method provided in the embodiments of this specification correspond to each other. Therefore, the device, electronic device, and non-volatile computer storage medium also have beneficial technical effects similar to the corresponding method.
  • the beneficial technical effects of the method have been described in detail above, therefore, the beneficial technical effects of the corresponding device, electronic equipment, and non-volatile computer storage medium will not be repeated here.
  • the improvement of a technology can be clearly distinguished between hardware improvements (for example, improvements in circuit structures such as diodes, transistors, switches, etc.) or software improvements (improvements in method flow).
  • hardware improvements for example, improvements in circuit structures such as diodes, transistors, switches, etc.
  • software improvements improvements in method flow.
  • the improvement of many methods and processes of today can be regarded as a direct improvement of the hardware circuit structure.
  • Designers almost always get the corresponding hardware circuit structure by programming the improved method flow into the hardware circuit. Therefore, it cannot be said that the improvement of a method flow cannot be realized by the hardware entity module.
  • a programmable logic device for example, a Field Programmable Gate Array (Field Programmable Gate Array, FPGA)
  • PLD Programmable Logic Device
  • FPGA Field Programmable Gate Array
  • HDL Hardware Description Language
  • ABEL Advanced Boolean Expression Language
  • AHDL Altera Hardware Description Language
  • HDCal JHDL
  • Lava Lava
  • Lola MyHDL
  • PALASM RHDL
  • VHDL Very-High-Speed Integrated Circuit Hardware Description Language
  • Verilog Verilog
  • the controller can be implemented in any suitable manner.
  • the controller can take the form of, for example, a microprocessor or a processor and a computer-readable medium storing computer-readable program codes (such as software or firmware) executable by the (micro)processor. , Logic gates, switches, application specific integrated circuits (ASICs), programmable logic controllers and embedded microcontrollers. Examples of controllers include but are not limited to the following microcontrollers: ARC625D, Atmel AT91SAM, Microchip PIC18F26K20 and Silicon Labs C8051F320, the memory controller can also be implemented as part of the memory control logic.
  • controllers in addition to implementing the controller in a purely computer-readable program code manner, it is entirely possible to program the method steps to make the controller use logic gates, switches, application-specific integrated circuits, programmable logic controllers, and embedded logic.
  • the same function can be realized in the form of a microcontroller or the like. Therefore, such a controller can be regarded as a hardware component, and the devices included in it for realizing various functions can also be regarded as a structure within the hardware component. Or even, the device for realizing various functions can be regarded as both a software module for realizing the method and a structure within a hardware component.
  • a typical implementation device is a computer.
  • the computer may be, for example, a personal computer, a laptop computer, a cell phone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or Any combination of these devices.
  • the embodiments of this specification can be provided as a method, a system, or a computer program product. Therefore, the embodiments of this specification may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the embodiments of this specification may adopt the form of computer program products implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.
  • computer-usable storage media including but not limited to disk storage, CD-ROM, optical storage, etc.
  • These computer program instructions can also be stored in a computer-readable memory that can guide a computer or other programmable data processing equipment to work in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture including the instruction device.
  • the device implements the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
  • These computer program instructions can also be loaded on a computer or other programmable data processing equipment, so that a series of operation steps are executed on the computer or other programmable equipment to produce computer-implemented processing, so as to execute on the computer or other programmable equipment.
  • the instructions provide steps for implementing the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
  • the computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
  • processors CPUs
  • input/output interfaces network interfaces
  • memory volatile and non-volatile memory
  • the memory may include non-permanent memory in a computer readable medium, random access memory (RAM) and/or non-volatile memory, such as read-only memory (ROM) or flash memory (flash RAM). Memory is an example of computer readable media.
  • RAM random access memory
  • ROM read-only memory
  • flash RAM flash memory
  • Computer-readable media include permanent and non-permanent, removable and non-removable media, and information storage can be realized by any method or technology.
  • the information can be computer-readable instructions, data structures, program modules, or other data.
  • Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, CD-ROM, digital versatile disc (DVD) or other optical storage, Magnetic cassettes, magnetic tape storage or other magnetic storage devices or any other non-transmission media can be used to store information that can be accessed by computing devices. According to the definition in this article, computer-readable media does not include transitory media, such as modulated data signals and carrier waves.
  • program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types.
  • the instructions can also be practiced in distributed computing environments where tasks are performed by remote processing devices connected through a communication network.
  • program modules can be located in local and remote computer storage media including storage devices.

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Pathology (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

A method and a system for assessing an aneurysm rupture risk comprising: obtaining data to be processed; inputting 3D triangle grid data and/or second data and/or third data into an aneurysm risk assessment model, to obtain an aneurysm risk assessment result of the data to be processed, wherein the aneurysm risk assessment model is a digital model obtained by means of training on the basis of labeled 3D triangle grid data and/or the second data and/or the third data; and outputting the aneurysm risk assessment result of the data to be processed. The method can eliminate or reduce participation of human factors, reduce time-consumption, achieve simple and quick assessment of an aneurysm rupture risk, and provide objective support for assessment of the aneurysm rupture risk.

Description

一种动脉瘤破裂风险评估方法及系统A method and system for evaluating aneurysm rupture risk 技术领域Technical field
本说明书涉及医学影像和计算机技术领域,尤其涉及一种动脉瘤破裂风险评估方法及系统。This specification relates to the field of medical imaging and computer technology, and in particular to a method and system for assessing the risk of aneurysm rupture.
背景技术Background technique
颅内动脉瘤多为发生在颅内动脉管壁上的异常膨出,据报道,颅内未破裂动脉瘤在我国成人中患病率可高达7%,破裂后造成蛛网膜下腔出血,可导致严重残疾或死亡。颅内动脉瘤破裂后治疗主要是清除颅内血肿并阻止血液继续流向颅内;颅内未破裂动脉瘤需要根据动脉瘤破裂风险评估结果及被观察者的病情,制定个性化管理方案,进行保守观察或手术干预。因此,动脉瘤破裂风险评估具有重要的意义。Intracranial aneurysms are mostly abnormal bulges that occur on the walls of intracranial arteries. According to reports, the prevalence of unruptured intracranial aneurysms in Chinese adults can be as high as 7%. Subarachnoid hemorrhage can occur after rupture. Cause severe disability or death. The treatment of intracranial aneurysm after rupture is mainly to remove intracranial hematoma and prevent blood from continuing to flow into the skull; intracranial unruptured aneurysms need to develop a personalized management plan based on the results of the aneurysm rupture risk assessment and the condition of the observed person, and be conservative Observation or surgical intervention. Therefore, the risk assessment of aneurysm rupture is of great significance.
目前,由观察者根据三维DSA(Digital subtraction angiography,数字减影血管造影)造影图像或MRA(Magnetic Resonance Angiography,磁共振血管成像)或CTA(CT angiography,CT血管造影)图像观察动脉瘤形态,并结合被观察者的病情,评估动脉瘤破裂风险,这种方法往往依赖观察者的经验,受主观判断影响,缺少客观支持,耗时久。At present, the observer observes the aneurysm morphology based on the three-dimensional DSA (Digital subtraction angiography) image, MRA (Magnetic Resonance Angiography) or CTA (CT angiography, CT angiography) image, and Combining the condition of the observer to assess the risk of aneurysm rupture, this method often relies on the experience of the observer, is affected by subjective judgment, lacks objective support, and takes a long time.
因此,需要一种新的方法,能够排除或减少人为因素的参与,缩短耗时,实现简单快捷的进行动脉瘤破裂风险评估,为动脉瘤破裂风险评估提供客观支持。Therefore, a new method is needed that can eliminate or reduce the participation of human factors, shorten the time-consuming, realize simple and quick assessment of the risk of aneurysm rupture, and provide objective support for the assessment of the risk of aneurysm rupture.
发明内容Summary of the invention
本说明书实施例提供一种动脉瘤破裂风险评估方法及系统,用于解决以下技术问题:现有技术,由观察者根据三维DSA(Digital subtraction angiography,数字减影血管造影)造影图像或MRA(Magnetic Resonance Angiography,磁 共振血管成像)或CTA(CT angiography,CT血管造影)图像观察动脉瘤形态,并结合被观察者的病情,评估动脉瘤破裂风险,这种方法往往依赖观察者的经验,受主观判断影响,缺少客观支持,耗时久。The embodiments of this specification provide an aneurysm rupture risk assessment method and system, which are used to solve the following technical problems: in the prior art, the observer uses three-dimensional DSA (Digital subtraction angiography, digital subtraction angiography) angiography images or MRA (Magnetic Subtraction Angiography) images or MRA (Magnetic Subtraction Angiography) images according to the prior art. Resonance Angiography, magnetic resonance angiography) or CTA (CT angiography, CT angiography) image to observe the aneurysm morphology, combined with the condition of the subject to assess the risk of aneurysm rupture, this method often depends on the experience of the observer, subject to subjectivity It takes a long time to judge the impact without objective support.
本说明书实施例提供一种动脉瘤破裂风险评估方法,包括以下步骤:The embodiment of this specification provides a method for assessing the risk of aneurysm rupture, which includes the following steps:
获取待处理数据,其中,所述待处理数据包括第一数据和/或第二数据和/或第三数据;Acquiring data to be processed, where the data to be processed includes first data and/or second data and/or third data;
将所述第一数据转换为3D三角形网格数据,基于所述3D三角形网格数据和/或所述第二数据和/或所述第三数据输入动脉瘤风险评估模型,获得所述待处理数据的动脉瘤风险评估结果,其中,所述动脉瘤风险评估模型是基于已标注的3D三角形网格数据和/或所述第二数据和/或所述第三数据训练获得的数字模型,所述动脉瘤风险评估结果包括所述待处理数据的动脉瘤破裂风险概率和/或所述待处理数据对应的3D三角形网格数据的顶点分类,所述待处理数据对应的3D三角形网格数据的顶点分类为所述待处理数据对应的3D三角形网格数据的顶点是属于血管或者属于动脉瘤;Convert the first data into 3D triangle mesh data, and input the aneurysm risk assessment model based on the 3D triangle mesh data and/or the second data and/or the third data to obtain the to-be-processed Data based on the aneurysm risk assessment result, wherein the aneurysm risk assessment model is based on the labeled 3D triangle grid data and/or the second data and/or the digital model obtained by training on the third data, so The aneurysm risk assessment result includes the aneurysm rupture risk probability of the to-be-processed data and/or the vertex classification of the 3D triangular mesh data corresponding to the to-be-processed data, and the value of the 3D triangular mesh data corresponding to the to-be-processed data The vertex is classified as whether the vertex of the 3D triangle mesh data corresponding to the to-be-processed data belongs to a blood vessel or an aneurysm;
输出所述待处理数据的动脉瘤风险评估结果。Output the aneurysm risk assessment result of the to-be-processed data.
进一步地,所述已标注的3D三角形网格数据是基于已知的3D三角形网格数据,将所述已知的3D三角形网格数据中属于血管的顶点标注为0,将所述已知的3D三角形网格数据中属于动脉瘤的顶点标注为1。Further, the labeled 3D triangle mesh data is based on the known 3D triangle mesh data, the vertices belonging to the blood vessel in the known 3D triangle mesh data are marked as 0, and the known 3D triangle mesh data is marked as 0. The vertex belonging to the aneurysm in the 3D triangle mesh data is marked as 1.
进一步地,所述已标注的3D三角形网格数据进一步包括:Further, the labeled 3D triangle mesh data further includes:
对所述已标注的3D三角形网格数据按照预设比率随机舍弃属于血管的顶点和/或属于动脉瘤的顶点,以获得属于动脉瘤点周围的顶点坐标、以及所述属于动脉瘤点周围的顶点之间相连的边信息。The vertices belonging to the blood vessel and/or the vertices belonging to the aneurysm are randomly discarded for the labeled 3D triangle mesh data according to the preset ratio to obtain the coordinates of the vertices belonging to the aneurysm point and the points belonging to the aneurysm. Edge information connected between vertices.
进一步地,所述动脉瘤风险评估模型的训练包括:Further, the training of the aneurysm risk assessment model includes:
基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,其中,所述一维向量表示所述已标注的3D三角形网格数据的全局特征向量;Based on the labeled 3D triangle mesh data, a one-dimensional vector of the labeled 3D triangle mesh data is obtained through an arithmetic block operation, where the one-dimensional vector represents the labeled 3D triangle mesh The global feature vector of the data;
基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,以获得所述动脉瘤风险评估模型。Based on the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle grid data, the aneurysm risk assessment of the labeled 3D triangle grid data is obtained through a sigmoid function operation As a result, the aneurysm risk assessment model is obtained.
进一步地,所述基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,具体包括:Further, the obtaining a one-dimensional vector of the marked 3D triangular mesh data through an arithmetic block operation based on the marked 3D triangular mesh data specifically includes:
将所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块以及(64,128,1024)维度的第二运算块的运算,获得(N,1024)的矩阵,其中N为所述已标注的3D三角形网格数据中属于顶点的个数;The labeled 3D triangle mesh data is subjected to the operations of the first operation block of (64, 64) dimensions and the second operation block of (64, 128, 1024) dimensions to obtain a matrix of (N, 1024), where N is The number of vertices in the labeled 3D triangle mesh data;
将所述(N,1024)的矩阵经过最大值池化层降维,获得所述已标注的3D三角形网格数据的一维向量。The dimension reduction of the matrix of (N, 1024) through the maximum pooling layer, to obtain the one-dimensional vector of the labeled 3D triangular mesh data.
进一步地,所述基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,具体包括:Further, the second data and/or third data corresponding to the one-dimensional vector and the labeled 3D triangle mesh data are operated by a sigmoid function to obtain the labeled 3D triangle mesh data The results of the aneurysm risk assessment include:
将所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据输入维度为(512,256,2)的多层感知机,获得第四数据;Input the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle grid data into a multi-layer perceptron with a dimension of (512, 256, 2) to obtain fourth data;
将所述第四数据经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤破裂风险概率;Subjecting the fourth data to a sigmoid function operation to obtain the aneurysm rupture risk probability of the marked 3D triangle mesh data;
和/或and / or
将所述一维向量及所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块获得的新的矩阵,经过维度为(512,256,128)的第三运算块和维度为(128,2)的第四运算块,获得(N,2)的矩阵;A new matrix obtained by passing the one-dimensional vector and the labeled 3D triangle mesh data through the first operation block of (64, 64) dimensions, passing through the third operation block of dimensions (512, 256, 128), and the dimension is ( 128,2) the fourth operation block to obtain the matrix of (N,2);
将所述(N,2)的矩阵经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的顶点分类。The matrix of (N, 2) is subjected to a sigmoid function operation to obtain the vertex classification of the labeled 3D triangle mesh data.
本说明书实施例提供的一种动脉瘤破裂风险评估系统,包括:An aneurysm rupture risk assessment system provided by an embodiment of this specification includes:
输入模块,获取待处理数据,其中,所述待处理数据包括第一数据和/或第二数据和/或第三数据;An input module to obtain data to be processed, where the data to be processed includes first data and/or second data and/or third data;
评估模块,将所述第一数据转换为3D三角形网格数据,基于所述3D三角形网格数据和/或所述第二数据和/或所述第三数据输入动脉瘤风险评估模型,获得所述待处理数据的动脉瘤风险评估结果,其中,所述动脉瘤风险评估模型是基于已标注的3D三角形网格数据和/或所述第二数据和/或所述第三数据训练获得的数字模型,所述动脉瘤风险评估结果包括所述待处理数据的动脉瘤破裂风险概率和/或所述待处理数据对应的3D三角形网格数据的顶点分类,所述待处理数据对应的3D三角形网格数据的顶点分类为所述待处理数据对应的3D三角形网格数据的顶点是属于血管或者属于动脉瘤;The evaluation module converts the first data into 3D triangle mesh data, and inputs the aneurysm risk assessment model based on the 3D triangle mesh data and/or the second data and/or the third data to obtain all The aneurysm risk assessment result of the to-be-processed data, wherein the aneurysm risk assessment model is based on the marked 3D triangle grid data and/or the second data and/or the number obtained by training on the third data Model, the aneurysm risk assessment result includes the aneurysm rupture risk probability of the data to be processed and/or the vertex classification of the 3D triangle mesh data corresponding to the data to be processed, and the 3D triangle mesh corresponding to the data to be processed The vertex of the grid data is classified as whether the vertex of the 3D triangle mesh data corresponding to the data to be processed belongs to a blood vessel or an aneurysm;
输出模块,输出所述待处理数据的动脉瘤风险评估结果。The output module outputs the aneurysm risk assessment result of the to-be-processed data.
进一步地,所述已标注的3D三角形网格数据是基于已知的3D三角形网格数据,将所述已知的3D三角形网格数据中属于血管的顶点标注为0,将所述已知的3D三角形网格数据中属于动脉瘤的顶点标注为1。Further, the labeled 3D triangle mesh data is based on the known 3D triangle mesh data, the vertices belonging to the blood vessel in the known 3D triangle mesh data are marked as 0, and the known 3D triangle mesh data is marked as 0. The vertex belonging to the aneurysm in the 3D triangle mesh data is marked as 1.
进一步地,所述已标注的3D三角形网格数据进一步包括:Further, the labeled 3D triangle mesh data further includes:
对所述已标注的3D三角形网格数据按照预设比率随机舍弃属于血管的顶点和/或属于动脉瘤的顶点,以获得属于动脉瘤点周围的顶点坐标、以及所述属于动脉瘤点周围的顶点之间相连的边信息。The vertices belonging to the blood vessel and/or the vertices belonging to the aneurysm are randomly discarded for the labeled 3D triangle mesh data according to the preset ratio to obtain the coordinates of the vertices belonging to the aneurysm point and the points belonging to the aneurysm. Edge information connected between vertices.
进一步地,所述动脉瘤风险评估模型的训练包括:Further, the training of the aneurysm risk assessment model includes:
基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,其中,所述一维向量表示所述已标注的3D三角形网格数据的全局特征向量;Based on the labeled 3D triangle mesh data, a one-dimensional vector of the labeled 3D triangle mesh data is obtained through an arithmetic block operation, where the one-dimensional vector represents the labeled 3D triangle mesh The global feature vector of the data;
基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,以获得所述动脉瘤风险评估模型。Based on the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle grid data, the aneurysm risk assessment of the labeled 3D triangle grid data is obtained through a sigmoid function operation As a result, the aneurysm risk assessment model is obtained.
进一步地,所述基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,具体包括:Further, the obtaining a one-dimensional vector of the marked 3D triangular mesh data through an arithmetic block operation based on the marked 3D triangular mesh data specifically includes:
将所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块以及 (64,128,1024)维度的第二运算块的运算,获得(N,1024)的矩阵,其中N为所述已标注的3D三角形网格数据中属于顶点的个数;The labeled 3D triangle mesh data is subjected to the operations of the first operation block of (64, 64) dimensions and the second operation block of (64, 128, 1024) dimensions to obtain a matrix of (N, 1024), where N is The number of vertices in the labeled 3D triangle mesh data;
将所述(N,1024)的矩阵经过最大值池化层降维,获得所述已标注的3D三角形网格数据的一维向量。The dimension reduction of the matrix of (N, 1024) through the maximum pooling layer, to obtain the one-dimensional vector of the labeled 3D triangular mesh data.
进一步地,所述基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,具体包括:Further, the second data and/or third data corresponding to the one-dimensional vector and the labeled 3D triangle mesh data are operated by a sigmoid function to obtain the labeled 3D triangle mesh data The results of the aneurysm risk assessment include:
将所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据输入维度为(512,256,2)的多层感知机,获得第四数据;Input the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle grid data into a multi-layer perceptron with a dimension of (512, 256, 2) to obtain fourth data;
将所述第四数据经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤破裂风险概率;Subjecting the fourth data to a sigmoid function operation to obtain the aneurysm rupture risk probability of the marked 3D triangle mesh data;
和/或and / or
将所述一维向量及所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块获得的新的矩阵,经过维度为(512,256,128)的第三运算块和维度为(128,2)的第四运算块,获得(N,2)的矩阵;A new matrix obtained by passing the one-dimensional vector and the labeled 3D triangle mesh data through the first operation block of (64, 64) dimensions, passing through the third operation block of dimensions (512, 256, 128), and the dimension is ( 128,2) the fourth operation block to obtain the matrix of (N,2);
将所述(N,2)的矩阵经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的顶点分类。The matrix of (N, 2) is subjected to a sigmoid function operation to obtain the vertex classification of the labeled 3D triangle mesh data.
本说明书实施例采用的上述至少一个技术方案能够达到以下有益效果:The above at least one technical solution adopted in the embodiment of this specification can achieve the following beneficial effects:
本说明书实施例通过获取待处理数据,其中,所述待处理数据包括第一数据和/或第二数据和/或第三数据;所述第一数据转换为3D三角形网格数据,基于所述3D三角形网格数据和/或所述第二数据和/或所述第三数据输入动脉瘤风险评估模型,获得所述待处理数据的动脉瘤风险评估结果,其中,所述动脉瘤风险评估模型是基于已标注的3D三角形网格数据和/或所述第二数据和/或所述第三数据训练获得的数字模型,所述动脉瘤风险评估结果包括所述待处理数据的动脉瘤破裂风险概率和/或所述待处理数据对应的3D三角形网格数据的顶点分类,所述待处理数据对应的3D三角形网格数据的顶点分类为所述待处 理数据对应的3D三角形网格数据的顶点是属于血管或者属于动脉瘤;输出所述待处理数据的动脉瘤风险评估结果,能够排除或减少人为因素的参与,缩短耗时,实现简单快捷的进行动脉瘤破裂风险评估,为动脉瘤破裂风险评估提供客观支持。The embodiment of this specification acquires data to be processed, wherein the data to be processed includes first data and/or second data and/or third data; the first data is converted into 3D triangle mesh data based on the The 3D triangle grid data and/or the second data and/or the third data are input into an aneurysm risk assessment model to obtain an aneurysm risk assessment result of the to-be-processed data, wherein the aneurysm risk assessment model It is a digital model obtained by training based on the labeled 3D triangle grid data and/or the second data and/or the third data, and the aneurysm risk assessment result includes the aneurysm rupture risk of the data to be processed Probability and/or the vertex classification of the 3D triangle mesh data corresponding to the data to be processed, the vertex of the 3D triangle mesh data corresponding to the data to be processed is classified as the vertex of the 3D triangle mesh data corresponding to the data to be processed It belongs to a blood vessel or an aneurysm; outputting the aneurysm risk assessment result of the data to be processed can eliminate or reduce the participation of human factors, shorten time-consuming, and realize simple and quick assessment of aneurysm rupture risk, which is the risk of aneurysm rupture The evaluation provides objective support.
附图说明Description of the drawings
为了更清楚地说明本说明书实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本说明书中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly describe the technical solutions in the embodiments of this specification or the prior art, the following will briefly introduce the drawings that need to be used in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only These are some embodiments described in this specification. For those of ordinary skill in the art, other drawings can be obtained from these drawings without creative labor.
图1为本说明书实施例提供的一种动脉瘤破裂风险评估方法的流程示意图;Fig. 1 is a schematic flow chart of a method for assessing the risk of aneurysm rupture according to an embodiment of this specification;
图2为本说明书实施例提供的动脉瘤风险评估模型的构建流程示意图;Figure 2 is a schematic diagram of the construction process of the aneurysm risk assessment model provided by the embodiments of this specification;
图3为本说明书实施例提供的一种动脉瘤破裂风险评估系统的示意图。Fig. 3 is a schematic diagram of an aneurysm rupture risk assessment system provided by an embodiment of this specification.
具体实施方式Detailed ways
颅内动脉瘤形态学参数对于颅内动脉瘤的诊断具有重要意义,颅内动脉瘤形态学参数可以通过医学影像获得,包括但不限于DSA、MRA、CTA。The morphological parameters of intracranial aneurysms are of great significance for the diagnosis of intracranial aneurysms. The morphological parameters of intracranial aneurysms can be obtained through medical imaging, including but not limited to DSA, MRA, and CTA.
DSA的基本原理是将注入造影剂前后拍摄的两帧X线图像经数字化输入图像计算机,通过减影、增强和再成像过程来获得清晰的纯血管影像,同时实时地显现血管影。DSA具有对比度分辨率高、检查时间短、造影剂用量少,浓度低、患者X线吸收量明显降低以及节省胶片等优点,在血管疾患的临床诊断中,具有十分重要的意义。DSA因其成像特点成为颅内动脉血管畸形和动脉瘤诊断的金标准。The basic principle of DSA is to digitally input two frames of X-ray images taken before and after the injection of the contrast agent into the image computer. Through the process of subtraction, enhancement and re-imaging, a clear pure blood vessel image is obtained, and the vascular shadow is displayed in real time. DSA has the advantages of high contrast resolution, short examination time, low contrast agent consumption, low concentration, significantly reduced patient X-ray absorption, and film saving. It is of great significance in the clinical diagnosis of vascular diseases. DSA has become the gold standard for the diagnosis of intracranial arterial malformations and aneurysms due to its imaging characteristics.
MRA基本原理是基于饱和效应、流入增强效应、流动去相位效应。MRA是将预饱和带置于3D层块的头端以饱和静脉血流,反向流动的动脉血液进入3D层块,因未被饱和从而产生MR信号。扫描时将一个较厚容积分割成多个薄 层激发,减少激发容积厚度以减少流入饱和效应,且能保证扫描容积范围,获得数层相邻层面的薄层图像,使图像清晰,血管的细微结构显示好,空间分辨力提高。MRA因其高质量的成像特点,也逐步用于颅内动脉瘤的诊断。The basic principle of MRA is based on saturation effect, inflow enhancement effect, and flow dephasing effect. In MRA, the pre-saturation zone is placed at the head end of the 3D slab to saturate the venous blood flow. The arterial blood flowing in the reverse direction enters the 3D slab and generates MR signals because it is not saturated. During scanning, a thicker volume is divided into multiple thin layers for excitation, reducing the thickness of the excitation volume to reduce the inflow saturation effect, and can ensure the scanning volume range, obtain several layers of adjacent layers of thin-layer images, make the image clear, and the blood vessels are subtle The structure is displayed well and the spatial resolution is improved. Due to its high-quality imaging characteristics, MRA is gradually used in the diagnosis of intracranial aneurysms.
CTA的基本原理是静脉注射含碘造影剂后,利用螺旋CT或电子束CT在造影剂充盈的受检血管高峰期连续进行薄层扫描,快速获取大量薄层叠加断面,经计算机图像处理后重建血管立体影像,清晰显示大脑前动脉、中动脉、后动脉及其主要分支和Wi11is动脉环等。静脉注射含碘造影剂后,利用螺旋CT或电子束CT在造影剂充盈的受检血管高峰期连续进行薄层扫描,快速获取大量薄层叠加断面,经计算机图像处理后重建血管立体影像,清晰显示大脑前动脉、中动脉、后动脉及其主要分支和Wi11is动脉环等。由于CTA具有无创、快捷、操作简单、价格低廉等特点,在临床颅内动脉瘤的诊断中已可以大部分取代DSA。The basic principle of CTA is that after intravenous injection of iodine-containing contrast agent, spiral CT or electron beam CT is used to continuously perform thin-slice scanning at the peak of the blood vessel filled with the contrast agent to quickly obtain a large number of thin-layer superimposed sections, which are reconstructed after computer image processing Stereoscopic images of blood vessels, clearly showing the anterior cerebral artery, middle artery, posterior artery and its main branches, Wi11is arterial ring, etc. After intravenous injection of iodine-containing contrast agent, spiral CT or electron beam CT is used to continuously perform thin-layer scanning at the peak of the blood vessel filled with the contrast agent to quickly obtain a large number of thin-layer superimposed sections, and reconstruct the three-dimensional image of the blood vessel after computer image processing. Show the anterior cerebral artery, middle artery, posterior artery and its main branches, Wi11is arterial ring, etc. Because CTA has the characteristics of non-invasive, fast, simple operation and low price, it can mostly replace DSA in the diagnosis of clinical intracranial aneurysms.
本申请,基于DSA和/或CTA和/或MRA原始图像可以经数字化处理为3D三角形网格数据,基于被观察者的动脉瘤形态学参数及被观察者的信息参数,可以实现对被评估值的动脉瘤破裂风险的评估。In this application, the original image based on DSA and/or CTA and/or MRA can be digitally processed into 3D triangular grid data. Based on the morphological parameters of the aneurysm of the observed person and the information parameters of the observed person, the estimated value can be realized Assessment of the risk of aneurysm rupture.
为了使本技术领域的人员更好地理解本说明书中的技术方案,下面将结合本说明书实施例中的附图,对本说明书实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本说明书实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。In order to enable those skilled in the art to better understand the technical solutions in this specification, the following will clearly and completely describe the technical solutions in the embodiments of this specification in conjunction with the drawings in the embodiments of this specification. Obviously, the described The embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments of this specification, all other embodiments obtained by those of ordinary skill in the art without creative work shall fall within the protection scope of this application.
以下结合附图,详细说明本申请各实施例提供的技术方案。The technical solutions provided by the embodiments of the present application will be described in detail below with reference to the accompanying drawings.
图1为本说明书实施例提供的一种动脉瘤破裂风险评估方法的流程示意图。该方法具体包括以下步骤:Fig. 1 is a schematic flow chart of a method for assessing the risk of aneurysm rupture according to an embodiment of this specification. The method specifically includes the following steps:
步骤S101:获取待处理数据,其中,所述待处理数据包括第一数据和/或第二数据和/或第三数据。Step S101: Obtain data to be processed, where the data to be processed includes first data and/or second data and/or third data.
在本申请中,第一数据和/或第二数据和/或第三数据是每个动脉瘤样本对 应的数据In this application, the first data and/or the second data and/or the third data are the data corresponding to each aneurysm sample
在本申请中,第一数据为三维影像数据,包括但不限于DSA影像数据、CTA影像数据、MRA影像数据。In this application, the first data is 3D image data, including but not limited to DSA image data, CTA image data, and MRA image data.
在本申请中,第二数据包括但不限于:被观察者的年龄和/或性别和/或是否为多发动脉瘤患者和/或饮酒史和/或吸烟史和/或家族史和/或动脉瘤在血管分段中的位置和/或症状。In this application, the second data includes but is not limited to: the age and/or gender of the observed person and/or whether it is a patient with multiple aneurysms and/or history of drinking and/or smoking and/or family history and/or arteries The location and/or symptoms of the tumor in the vascular segment.
在本申请中,第三数据为被观察者的动脉瘤形态学参数,包括但不限于:动脉瘤体积、载瘤血管平均直径、瘤体长度与瘤颈直径比值(SR)、瘤体长度与瘤颈宽度的比值(AR)、动脉瘤长径、动脉瘤高度、动脉瘤宽度、动脉瘤颈宽度、流入角度、瘤颈面积、瘤颈直径、载瘤血管长度、波动指数(UI)、非球形指数(NSI)。In this application, the third data is the morphological parameters of the aneurysm of the observed person, including but not limited to: aneurysm volume, average diameter of tumor-bearing blood vessels, ratio of tumor body length to tumor neck diameter (SR), tumor body length and Ratio of aneurysm neck width (AR), aneurysm length, aneurysm height, aneurysm width, aneurysm neck width, inflow angle, aneurysm neck area, aneurysm neck diameter, tumor-carrying vessel length, fluctuation index (UI), non- Spherical Index (NSI).
在本申请实施例中,将第二数据和第三数据以表格的形式存储,以便于后续使用。In the embodiment of the present application, the second data and the third data are stored in the form of a table for subsequent use.
在本申请的一个实施例中,基于第一数据,获得被观察者的动脉瘤形态学参数。具体地,采用阈值分割和区域增长的方法,将第一数据转换成二值图;使用Marching Cube算法构建3D三角形网格数据;基于已经构建的3D三角形网格数据,获得被观察者的动脉瘤形态学参数。3D三角形网格数据是空间中笛卡尔坐标系下多个三角形组成的拓扑结构,其中每个三角形包含顶点和边。In an embodiment of the present application, based on the first data, the morphological parameters of the aneurysm of the observed person are obtained. Specifically, the method of threshold segmentation and region growth is adopted to convert the first data into a binary image; the Marching Cube algorithm is used to construct 3D triangular mesh data; based on the constructed 3D triangular mesh data, the observed aneurysm is obtained Morphological parameters. The 3D triangle mesh data is a topological structure composed of multiple triangles in a Cartesian coordinate system in space, where each triangle contains vertices and edges.
在本申请的一个实施例中,3D三角形网格数据采用STL格式数据,具体地,可以为二进制文本或者文本文件。In an embodiment of the present application, the 3D triangular mesh data adopts STL format data, specifically, it may be a binary text or a text file.
由于被观察者可能会存在一个或多个动脉瘤,因此,待处理数据包括一个或多个动脉瘤对应的数据。Since the observed person may have one or more aneurysms, the data to be processed includes data corresponding to one or more aneurysms.
在本申请中,被观察者是本申请中动脉瘤破裂风险评估的对象。In this application, the observed person is the subject of the risk assessment of aneurysm rupture in this application.
步骤S103:将所述第一数据转换为3D三角形网格数据,基于所述3D三角形网格数据和/或所述第二数据和/或所述第三数据输入动脉瘤风险评估模型,获得所述待处理数据的动脉瘤风险评估结果,其中,所述动脉瘤风险评估模型 是基于已标注的3D三角形网格数据和/或所述第二数据和/或所述第三数据训练获得的数字模型,所述动脉瘤风险评估结果包括所述待处理数据的动脉瘤破裂风险概率和/或所述待处理数据对应的3D三角形网格数据的顶点分类,所述待处理数据对应的3D三角形网格数据的顶点分类为所述待处理数据对应的3D三角形网格数据的顶点是属于血管或者属于动脉瘤。Step S103: Convert the first data into 3D triangular mesh data, and input the aneurysm risk assessment model based on the 3D triangular mesh data and/or the second data and/or the third data to obtain all The aneurysm risk assessment result of the to-be-processed data, wherein the aneurysm risk assessment model is based on the marked 3D triangle grid data and/or the second data and/or the number obtained by training on the third data Model, the aneurysm risk assessment result includes the aneurysm rupture risk probability of the data to be processed and/or the vertex classification of the 3D triangle mesh data corresponding to the data to be processed, and the 3D triangle mesh corresponding to the data to be processed The vertex of the grid data is classified as whether the vertex of the 3D triangle mesh data corresponding to the data to be processed belongs to a blood vessel or an aneurysm.
在本申请中,待处理数据对应的3D三角形网格数据是将待处理数据采用阈值分割和区域增长的方法,将第一数据转换成二值图;使用Marching Cube算法获得的3D三角形网格数据。In this application, the 3D triangle mesh data corresponding to the data to be processed is the method of threshold segmentation and region growth for the data to be processed, and the first data is converted into a binary image; the 3D triangle mesh data obtained by the Marching Cube algorithm .
在本申请中,动脉瘤风险评估模型是基于已标注的3D三角形网格数据训练获得的数字模型。为理解动脉瘤风险评估模型的获取,图2为本说明书实施例提供的一种动脉瘤风险评估模型的获取流程示意图,具体包括:In this application, the aneurysm risk assessment model is a digital model obtained by training based on labeled 3D triangle mesh data. In order to understand the acquisition of an aneurysm risk assessment model, FIG. 2 is a schematic diagram of the acquisition process of an aneurysm risk assessment model provided in an embodiment of this specification, which specifically includes:
步骤S201:对3D三角形网格数据进行标注,获得已标注的3D三角形网格数据。Step S201: Annotate the 3D triangle mesh data to obtain the annotated 3D triangle mesh data.
基于已知的3D三角形网格数据,将所述已知的3D三角形网格数据中属于血管的顶点标注为0,将所述已知的3D三角形网格数据中属于动脉瘤的顶点标注为1。在本申请中,已知的3D三角形网格数据是基于作为训练集的第一数据,获得的对应的3D三角形网格数据Based on the known 3D triangle mesh data, the vertices belonging to the blood vessel in the known 3D triangle mesh data are marked as 0, and the vertices belonging to the aneurysm in the known 3D triangular mesh data are marked as 1 . In this application, the known 3D triangle mesh data is based on the corresponding 3D triangle mesh data obtained as the first data of the training set
在具体实施过程中,对已标注的3D三角形网格数据进一步:对所述已标注的3D三角形网格数据按照预设比率随机舍弃属于血管的顶点和/或属于动脉瘤的顶点,以获得属于动脉瘤点周围的顶点坐标、以及所述属于动脉瘤点周围的顶点之间相连的边信息。需要特别说明的是,预设比率是根据已标注的3D三角形网格数据量和/或顶点坐标的数量综合确定的。In the specific implementation process, the labeled 3D triangle mesh data is further: randomly discarded the vertices belonging to the blood vessel and/or the aneurysm according to the preset ratio of the labeled 3D triangle mesh data, so as to obtain the vertices belonging to the aneurysm. The coordinates of the vertices around the aneurysm point and the edge information connected between the vertices belonging to the aneurysm point. It should be particularly noted that the preset ratio is comprehensively determined according to the amount of 3D triangle mesh data that has been marked and/or the number of vertex coordinates.
在本申请一个实施例中,对已标注的3D三角形网格数据进一步:对所述已标注的3D三角形网格数据按照预设比率随机舍弃属于血管的顶点和/或属于动脉瘤的顶点,以获得属于动脉瘤点周围的1024个或者2048个顶点坐标、以及所述属于动脉瘤点周围的1024个或者2048个顶点之间相连的边信息。顶点 的具体个数1024个或者2048个仅为本申请的一个示例性说明,并不构成对本申请的限定。In an embodiment of the present application, the labeled 3D triangle mesh data is further: randomly discarded vertices belonging to blood vessels and/or vertices belonging to aneurysms for the labeled 3D triangle mesh data according to a preset ratio. Obtain the coordinates of 1024 or 2048 vertices around the aneurysm point and the edge information connected between the 1024 or 2048 vertices around the aneurysm point. The specific number of vertices 1024 or 2048 is only an exemplary description of this application, and does not constitute a limitation of this application.
在本申请实施例中,可以将已知的3D三角形网格数据输入相应的标注软件中,使网格数据能够实现可视化,从而逐点进行人工标注。为了保证标注结果的准确性,可以采用多人标注,标注人员的数量至少3人。进行标注的人员应该具有专业经验,以保证标注的准确性。In the embodiment of the present application, the known 3D triangle mesh data can be input into the corresponding labeling software, so that the mesh data can be visualized, so that manual labeling can be performed point by point. In order to ensure the accuracy of the labeling results, multi-person labeling can be used, and the number of labeling personnel shall be at least 3 people. The person making the labeling should have professional experience to ensure the accuracy of the labeling.
在本申请实施例中,3D三角形网格数据是以N*F的矩阵形式存在的,其中N为顶点的数量(例如N=1024表示有1024个顶点),F为特征数量(例如F=3表示有x、y、z三个方向上的特征)。In the embodiment of this application, the 3D triangle mesh data is in the form of a N*F matrix, where N is the number of vertices (for example, N=1024 means there are 1024 vertices), and F is the number of features (for example, F=3 Indicates that there are features in the three directions of x, y, and z).
步骤S203:基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,其中,所述一维向量表示所述已标注的3D三角形网格数据的全局特征向量。Step S203: Based on the labeled 3D triangular mesh data, a one-dimensional vector of the labeled 3D triangular mesh data is obtained through an arithmetic block operation, where the one-dimensional vector represents the labeled 3D The global feature vector of the triangle mesh data.
在本申请实施例中,将所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块以及(64,128,1024)维度的第二运算块的运算,获得(N,1024)的矩阵,其中N为所述已标注的3D三角形网格数据中属于顶点的个数;In the embodiment of the present application, the labeled 3D triangle mesh data is subjected to the operations of the first operation block of (64, 64) dimensions and the second operation block of (64, 128, 1024) dimensions to obtain (N, 1024) ), where N is the number of vertices in the labeled 3D triangle mesh data;
将所述(N,1024)的矩阵经过最大值池化层降维,获得所述已标注的3D三角形网格数据的一维向量。The dimension reduction of the matrix of (N, 1024) through the maximum pooling layer, to obtain the one-dimensional vector of the labeled 3D triangular mesh data.
步骤S205:基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,以获得所述动脉瘤风险评估模型。Step S205: Based on the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle mesh data, perform a sigmoid function operation to obtain the artery of the labeled 3D triangle mesh data An aneurysm risk assessment result to obtain the aneurysm risk assessment model.
在本申请实施例中,将所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据输入维度为(512,256,2)的多层感知机,获得第四数据;In the embodiment of the present application, the second data and/or third data corresponding to the one-dimensional vector and the labeled 3D triangular mesh data are input into a multi-layer perceptron with a dimension of (512, 256, 2) to obtain Fourth data
将所述第四数据经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤破裂风险概率;Subjecting the fourth data to a sigmoid function operation to obtain the aneurysm rupture risk probability of the marked 3D triangle mesh data;
和/或and / or
将所述一维向量及所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块获得的新的矩阵,经过维度为(512,256,128)的第三运算块和维度为(128,2)的第四运算块,获得(N,2)的矩阵;A new matrix obtained by passing the one-dimensional vector and the labeled 3D triangle mesh data through the first operation block of (64, 64) dimensions, passing through the third operation block of dimensions (512, 256, 128), and the dimension is ( 128,2) the fourth operation block to obtain the matrix of (N,2);
将所述(N,2)的矩阵经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的顶点分类。The matrix of (N, 2) is subjected to a sigmoid function operation to obtain the vertex classification of the labeled 3D triangle mesh data.
在本申请实施例中,sigmoid函数也叫Logistic函数,在本申请中,经过sigmoid分类函数的处理,可以保证输出的数据位于(0,1)内,进行动脉瘤破裂风险的评估。sigmoid函数亦可实现数据分类,以实现已标注的3D三角形网格数据的顶点分类。In the embodiments of this application, the sigmoid function is also called the Logistic function. In this application, after the sigmoid classification function is processed, it can be ensured that the output data is within (0,1), and the risk of aneurysm rupture can be assessed. The sigmoid function can also realize data classification to realize the vertex classification of the labeled 3D triangle mesh data.
在动脉瘤风险评估模型的建立过程中,用于动脉瘤风险评估模型建立的数据包括:训练集和测试集。训练集用于训练动脉瘤风险评估模型,测试集用于测试得到的动脉瘤风险评估模型的效果。During the establishment of the aneurysm risk assessment model, the data used for the establishment of the aneurysm risk assessment model includes: training set and test set. The training set is used to train the aneurysm risk assessment model, and the test set is used to test the effect of the obtained aneurysm risk assessment model.
本说明书实施例提供的动脉瘤风险评估模型能够自动学习已标注的3D三角形网格数据中每个顶点及其周围顶点的局部特征,进一步通过整合被观察者的第二数据和/或第三数据,从而实现全面、高效的预测动脉瘤破裂风险和顶点分类。The aneurysm risk assessment model provided by the embodiments of this specification can automatically learn the local features of each vertex and surrounding vertices in the labeled 3D triangle mesh data, and further integrate the second data and/or third data of the observed person , So as to achieve a comprehensive and efficient prediction of aneurysm rupture risk and vertex classification.
步骤S105:输出所述待处理数据的动脉瘤风险评估结果。Step S105: Output the aneurysm risk assessment result of the to-be-processed data.
在本申请中,动脉瘤风险评估结果的输出包括:所述待处理数据的动脉瘤破裂风险概率和/或所述待处理数据对应的3D三角形网格的顶点分类。In this application, the output of the aneurysm risk assessment result includes: the aneurysm rupture risk probability of the to-be-processed data and/or the vertex classification of the 3D triangle mesh corresponding to the to-be-processed data.
在本申请的一个实施例中,利用本申请提供的动脉瘤风险评估方法,输入待处理数据,得到被观察者的动脉瘤破裂风险的概率为70%,可以作为后续诊断治疗的辅助参考,决定是否实施手术治疗。In one embodiment of the present application, using the aneurysm risk assessment method provided in the present application, inputting the data to be processed, the probability of obtaining the aneurysm rupture risk of the observed person is 70%, which can be used as an auxiliary reference for subsequent diagnosis and treatment to determine Whether to perform surgical treatment.
本申请实施例提供的方法,可以实现输入一组或多组动脉瘤数据,得到动脉瘤破裂风险的评估结果。在本申请中,一组动脉瘤数据为一个动脉瘤所对应的数据,包括第一数据和/或第二数据和/或第三数据。The method provided in the embodiments of the present application can realize the input of one or more sets of aneurysm data to obtain an assessment result of the risk of aneurysm rupture. In this application, a set of aneurysm data is data corresponding to an aneurysm, including the first data and/or the second data and/or the third data.
本说明书实施例提供的方法,能够排除或减少人为因素的参与,缩短耗时, 实现简单快捷的进行动脉瘤破裂风险评估,为动脉瘤破裂风险评估提供客观支持。The method provided in the embodiments of this specification can eliminate or reduce the participation of human factors, shorten time-consuming, realize simple and quick assessment of the risk of aneurysm rupture, and provide objective support for the assessment of the risk of aneurysm rupture.
本说明书实施例提供的评估方法,在实际应用时,可以封装为软件,用于辅助观察者在进行未破裂动脉瘤治疗决策时,快速得到合理可靠的预测结果和/或顶点分类。The evaluation method provided in the embodiments of this specification can be packaged as software in actual application to assist observers in making decisions about treatment of unruptured aneurysms to quickly obtain reasonable and reliable prediction results and/or vertex classification.
基于同样的思路,本说明书实施例还提供了一种动脉瘤破裂风险评估系统,图3为本说明书实施例提供的一种动脉瘤破裂风险评估系统的示意图,该系统包括:Based on the same idea, the embodiment of this specification also provides an aneurysm rupture risk assessment system. FIG. 3 is a schematic diagram of an aneurysm rupture risk assessment system provided by the embodiment of this specification, and the system includes:
输入模块301,获取待处理数据,其中,所述待处理数据包括第一数据和/或第二数据和/或第三数据;The input module 301 obtains data to be processed, where the data to be processed includes first data and/or second data and/or third data;
评估模块303,将所述第一数据转换为3D三角形网格数据,基于所述3D三角形网格数据和/或所述第二数据和/或所述第三数据输入动脉瘤风险评估模型,获得所述待处理数据的动脉瘤风险评估结果,其中,所述动脉瘤风险评估模型是基于已标注的3D三角形网格数据和/或所述第二数据和/或所述第三数据训练获得的数字模型,所述动脉瘤风险评估结果包括所述待处理数据的动脉瘤破裂风险概率和/或所述待处理数据对应的3D三角形网格数据的顶点分类,所述待处理数据对应的3D三角形网格数据的顶点分类为所述待处理数据对应的3D三角形网格数据的顶点是属于血管或者属于动脉瘤;The evaluation module 303 converts the first data into 3D triangle mesh data, and inputs the aneurysm risk assessment model based on the 3D triangle mesh data and/or the second data and/or the third data to obtain The aneurysm risk assessment result of the to-be-processed data, wherein the aneurysm risk assessment model is obtained by training based on the marked 3D triangle grid data and/or the second data and/or the third data A digital model, where the aneurysm risk assessment result includes the aneurysm rupture risk probability of the to-be-processed data and/or the vertex classification of the 3D triangle mesh data corresponding to the to-be-processed data, and the 3D triangle corresponding to the to-be-processed data The vertex of the mesh data is classified as whether the vertex of the 3D triangle mesh data corresponding to the data to be processed belongs to a blood vessel or an aneurysm;
输出模块305,输出所述待处理数据的动脉瘤风险评估结果。The output module 305 outputs the aneurysm risk assessment result of the to-be-processed data.
进一步地,所述已标注的3D三角形网格数据是基于已知的3D三角形网格数据,将所述已知的3D三角形网格数据中属于血管的顶点标注为0,将所述已知的3D三角形网格数据中属于动脉瘤的顶点标注为1。Further, the labeled 3D triangle mesh data is based on the known 3D triangle mesh data, the vertices belonging to the blood vessel in the known 3D triangle mesh data are marked as 0, and the known 3D triangle mesh data is marked as 0. The vertex belonging to the aneurysm in the 3D triangle mesh data is marked as 1.
进一步地,所述已标注的3D三角形网格数据进一步包括:Further, the labeled 3D triangle mesh data further includes:
对所述已标注的3D三角形网格数据按照预设比率随机舍弃属于血管的顶点和/或属于动脉瘤的顶点,以获得属于动脉瘤点周围的顶点坐标、以及所述属于动脉瘤点周围的顶点之间相连的边信息。The vertices belonging to the blood vessel and/or the vertices belonging to the aneurysm are randomly discarded for the labeled 3D triangle mesh data according to the preset ratio to obtain the coordinates of the vertices belonging to the aneurysm point and the points belonging to the aneurysm. Edge information connected between vertices.
进一步地,所述动脉瘤风险评估模型的训练包括:Further, the training of the aneurysm risk assessment model includes:
基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,其中,所述一维向量表示所述已标注的3D三角形网格数据的全局特征向量;Based on the labeled 3D triangle mesh data, a one-dimensional vector of the labeled 3D triangle mesh data is obtained through an arithmetic block operation, where the one-dimensional vector represents the labeled 3D triangle mesh The global feature vector of the data;
基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,以获得所述动脉瘤风险评估模型。Based on the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle grid data, the aneurysm risk assessment of the labeled 3D triangle grid data is obtained through a sigmoid function operation As a result, the aneurysm risk assessment model is obtained.
进一步地,所述基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,具体包括:Further, the obtaining a one-dimensional vector of the marked 3D triangular mesh data through an arithmetic block operation based on the marked 3D triangular mesh data specifically includes:
将所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块以及(64,128,1024)维度的第二运算块的运算,获得(N,1024)的矩阵,其中N为所述已标注的3D三角形网格数据中属于顶点的个数;The labeled 3D triangle mesh data is subjected to the operations of the first operation block of (64, 64) dimensions and the second operation block of (64, 128, 1024) dimensions to obtain a matrix of (N, 1024), where N is The number of vertices in the labeled 3D triangle mesh data;
将所述(N,1024)的矩阵经过最大值池化层降维,获得所述已标注的3D三角形网格数据的一维向量。The dimension reduction of the matrix of (N, 1024) through the maximum pooling layer, to obtain the one-dimensional vector of the labeled 3D triangular mesh data.
进一步地,所述基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,具体包括:Further, the second data and/or third data corresponding to the one-dimensional vector and the labeled 3D triangle mesh data are operated by a sigmoid function to obtain the labeled 3D triangle mesh data The results of the aneurysm risk assessment include:
将所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据输入维度为(512,256,2)的多层感知机,获得第四数据;Input the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle grid data into a multi-layer perceptron with a dimension of (512, 256, 2) to obtain fourth data;
将所述第四数据经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤破裂风险概率;Subjecting the fourth data to a sigmoid function operation to obtain the aneurysm rupture risk probability of the marked 3D triangle mesh data;
和/或and / or
将所述一维向量及所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块获得的新的矩阵,经过维度为(512,256,128)的第三运算块和维度为(128,2)的第四运算块,获得(N,2)的矩阵;A new matrix obtained by passing the one-dimensional vector and the labeled 3D triangle mesh data through the first operation block of (64, 64) dimensions, passing through the third operation block of dimensions (512, 256, 128), and the dimension is ( 128,2) the fourth operation block to obtain the matrix of (N,2);
将所述(N,2)的矩阵经过sigmoid函数运算,获得所述已标注的3D三角 形网格数据的顶点分类。The matrix of (N, 2) is subjected to a sigmoid function operation to obtain the vertex classification of the labeled 3D triangular mesh data.
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。The foregoing describes specific embodiments of this specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps described in the claims can be performed in a different order than in the embodiments and still achieve desired results. In addition, the processes depicted in the drawings do not necessarily require the specific order or sequential order shown in order to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置、电子设备、非易失性计算机存储介质实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。The various embodiments in this specification are described in a progressive manner, and the same or similar parts between the various embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the device, electronic equipment, and non-volatile computer storage medium embodiments, since they are basically similar to the method embodiments, the description is relatively simple, and for related parts, please refer to the part of the description of the method embodiments.
本说明书实施例提供的装置、电子设备、非易失性计算机存储介质与方法是对应的,因此,装置、电子设备、非易失性计算机存储介质也具有与对应方法类似的有益技术效果,由于上面已经对方法的有益技术效果进行了详细说明,因此,这里不再赘述对应装置、电子设备、非易失性计算机存储介质的有益技术效果。The device, electronic device, non-volatile computer storage medium and method provided in the embodiments of this specification correspond to each other. Therefore, the device, electronic device, and non-volatile computer storage medium also have beneficial technical effects similar to the corresponding method. The beneficial technical effects of the method have been described in detail above, therefore, the beneficial technical effects of the corresponding device, electronic equipment, and non-volatile computer storage medium will not be repeated here.
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable Gate Array,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且,如今,取代手工地制作集成 电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language)等,目前最普遍使用的是VHDL(Very-High-Speed Integrated Circuit Hardware Description Language)与Verilog。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。In the 1990s, the improvement of a technology can be clearly distinguished between hardware improvements (for example, improvements in circuit structures such as diodes, transistors, switches, etc.) or software improvements (improvements in method flow). However, with the development of technology, the improvement of many methods and processes of today can be regarded as a direct improvement of the hardware circuit structure. Designers almost always get the corresponding hardware circuit structure by programming the improved method flow into the hardware circuit. Therefore, it cannot be said that the improvement of a method flow cannot be realized by the hardware entity module. For example, a programmable logic device (Programmable Logic Device, PLD) (for example, a Field Programmable Gate Array (Field Programmable Gate Array, FPGA)) is such an integrated circuit whose logic function is determined by the user's programming of the device. It is programmed by the designer to "integrate" a digital system on a PLD, without requiring the chip manufacturer to design and manufacture a dedicated integrated circuit chip. Moreover, nowadays, instead of manually making integrated circuit chips, this kind of programming is mostly realized with "logic compiler" software, which is similar to the software compiler used in program development and writing, but before compilation The original code must also be written in a specific programming language, which is called Hardware Description Language (HDL), and there is not only one type of HDL, but many types, such as ABEL (Advanced Boolean Expression Language) , AHDL (Altera Hardware Description Language), Confluence, CUPL (Cornell University Programming Language), HDCal, JHDL (Java Hardware Description Language), Lava, Lola, MyHDL, PALASM, RHDL (Ruby Hardware Description), etc., currently most commonly used It is VHDL (Very-High-Speed Integrated Circuit Hardware Description Language) and Verilog. It should also be clear to those skilled in the art that only a little logic programming of the method flow in the above-mentioned hardware description languages and programming into an integrated circuit can easily obtain the hardware circuit that implements the logic method flow.
控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。The controller can be implemented in any suitable manner. For example, the controller can take the form of, for example, a microprocessor or a processor and a computer-readable medium storing computer-readable program codes (such as software or firmware) executable by the (micro)processor. , Logic gates, switches, application specific integrated circuits (ASICs), programmable logic controllers and embedded microcontrollers. Examples of controllers include but are not limited to the following microcontrollers: ARC625D, Atmel AT91SAM, Microchip PIC18F26K20 and Silicon Labs C8051F320, the memory controller can also be implemented as part of the memory control logic. Those skilled in the art also know that, in addition to implementing the controller in a purely computer-readable program code manner, it is entirely possible to program the method steps to make the controller use logic gates, switches, application-specific integrated circuits, programmable logic controllers, and embedded logic. The same function can be realized in the form of a microcontroller or the like. Therefore, such a controller can be regarded as a hardware component, and the devices included in it for realizing various functions can also be regarded as a structure within the hardware component. Or even, the device for realizing various functions can be regarded as both a software module for realizing the method and a structure within a hardware component.
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实 体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的,计算机例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。The systems, devices, modules, or units explained in the above embodiments may be implemented by computer chips or entities, or implemented by products with certain functions. A typical implementation device is a computer. Specifically, the computer may be, for example, a personal computer, a laptop computer, a cell phone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or Any combination of these devices.
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本说明书一个或多个实施例时可以把各单元的功能在同一个或多个软件和/或硬件中实现。For the convenience of description, when describing the above device, the functions are divided into various units and described separately. Of course, when implementing one or more embodiments of this specification, the functions of each unit may be implemented in the same one or more software and/or hardware.
本领域内的技术人员应明白,本说明书实施例可提供为方法、系统、或计算机程序产品。因此,本说明书实施例可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art should understand that the embodiments of this specification can be provided as a method, a system, or a computer program product. Therefore, the embodiments of this specification may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the embodiments of this specification may adopt the form of computer program products implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.
本说明书是参照根据本说明书实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。This specification is described with reference to flowcharts and/or block diagrams of methods, devices (systems), and computer program products according to the embodiments of this specification. It should be understood that each process and/or block in the flowchart and/or block diagram, and the combination of processes and/or blocks in the flowchart and/or block diagram can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing equipment to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing equipment are generated It is a device that realizes the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions can also be stored in a computer-readable memory that can guide a computer or other programmable data processing equipment to work in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture including the instruction device. The device implements the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处 理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded on a computer or other programmable data processing equipment, so that a series of operation steps are executed on the computer or other programmable equipment to produce computer-implemented processing, so as to execute on the computer or other programmable equipment. The instructions provide steps for implementing the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。In a typical configuration, the computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。The memory may include non-permanent memory in a computer readable medium, random access memory (RAM) and/or non-volatile memory, such as read-only memory (ROM) or flash memory (flash RAM). Memory is an example of computer readable media.
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带式磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。Computer-readable media include permanent and non-permanent, removable and non-removable media, and information storage can be realized by any method or technology. The information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, CD-ROM, digital versatile disc (DVD) or other optical storage, Magnetic cassettes, magnetic tape storage or other magnetic storage devices or any other non-transmission media can be used to store information that can be accessed by computing devices. According to the definition in this article, computer-readable media does not include transitory media, such as modulated data signals and carrier waves.
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。It should also be noted that the terms "include", "include" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or equipment including a series of elements not only includes those elements, but also includes Other elements that are not explicitly listed, or they also include elements inherent to such processes, methods, commodities, or equipment. If there are no more restrictions, the element defined by the sentence "including a..." does not exclude the existence of other identical elements in the process, method, commodity, or equipment that includes the element.
本说明书可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践 说明书,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。This specification may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types. The instructions can also be practiced in distributed computing environments where tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules can be located in local and remote computer storage media including storage devices.
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。The various embodiments in this specification are described in a progressive manner, and the same or similar parts between the various embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, as for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and for related parts, please refer to the part of the description of the method embodiment.
以上所述仅为本说明书实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。The above descriptions are only examples of this specification, and are not intended to limit this application. For those skilled in the art, this application can have various modifications and changes. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of this application shall be included in the scope of the claims of this application.

Claims (12)

  1. 一种动脉瘤破裂风险评估方法,其特征在于,所述评估方法包括:An aneurysm rupture risk assessment method, characterized in that the assessment method includes:
    获取待处理数据,其中,所述待处理数据包括第一数据和/或第二数据和/或第三数据;Acquiring data to be processed, where the data to be processed includes first data and/or second data and/or third data;
    将所述第一数据转换为3D三角形网格数据,基于所述3D三角形网格数据和/或所述第二数据和/或所述第三数据输入动脉瘤风险评估模型,获得所述待处理数据的动脉瘤风险评估结果,其中,所述动脉瘤风险评估模型是基于已标注的3D三角形网格数据和/或所述第二数据和/或所述第三数据训练获得的数字模型,所述动脉瘤风险评估结果包括所述待处理数据的动脉瘤破裂风险概率和/或所述待处理数据对应的3D三角形网格数据的顶点分类,所述待处理数据对应的3D三角形网格数据的顶点分类为所述待处理数据对应的3D三角形网格数据的顶点是属于血管或者属于动脉瘤;Convert the first data into 3D triangle mesh data, and input the aneurysm risk assessment model based on the 3D triangle mesh data and/or the second data and/or the third data to obtain the to-be-processed Data based on the aneurysm risk assessment result, wherein the aneurysm risk assessment model is based on the labeled 3D triangle grid data and/or the second data and/or the digital model obtained by training on the third data, so The aneurysm risk assessment result includes the aneurysm rupture risk probability of the to-be-processed data and/or the vertex classification of the 3D triangular mesh data corresponding to the to-be-processed data, and the value of the 3D triangular mesh data corresponding to the to-be-processed data The vertex is classified as whether the vertex of the 3D triangle mesh data corresponding to the to-be-processed data belongs to a blood vessel or an aneurysm;
    输出所述待处理数据的动脉瘤风险评估结果。Output the aneurysm risk assessment result of the to-be-processed data.
  2. 如权利要求1所述的方法,其特征在于,所述已标注的3D三角形网格数据是基于已知的3D三角形网格数据,将所述已知的3D三角形网格数据中属于血管的顶点标注为0,将所述已知的3D三角形网格数据中属于动脉瘤的顶点标注为1。The method of claim 1, wherein the labeled 3D triangle mesh data is based on known 3D triangle mesh data, and the vertices belonging to the blood vessel in the known 3D triangle mesh data It is marked as 0, and the vertex belonging to the aneurysm in the known 3D triangle mesh data is marked as 1.
  3. 如权利要求2所述的方法,其特征在于,所述已标注的3D三角形网格数据进一步包括:3. The method of claim 2, wherein the labeled 3D triangle mesh data further comprises:
    对所述已标注的3D三角形网格数据按照预设比率随机舍弃属于血管的顶点和/或属于动脉瘤的顶点,以获得属于动脉瘤点周围的顶点坐标、以及所述属于动脉瘤点周围的顶点之间相连的边信息。The vertices belonging to the blood vessel and/or the vertices belonging to the aneurysm are randomly discarded for the labeled 3D triangle mesh data according to the preset ratio to obtain the coordinates of the vertices belonging to the aneurysm point and the points belonging to the aneurysm. Edge information connected between vertices.
  4. 如权利要求1所述的方法,其特征在于,所述动脉瘤风险评估模型的训练包括:The method of claim 1, wherein the training of the aneurysm risk assessment model comprises:
    基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,其中,所述一维向量表示所述已标注的 3D三角形网格数据的全局特征向量;Based on the labeled 3D triangle mesh data, a one-dimensional vector of the labeled 3D triangle mesh data is obtained through an arithmetic block operation, where the one-dimensional vector represents the labeled 3D triangle mesh The global feature vector of the data;
    基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,以获得所述动脉瘤风险评估模型。Based on the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle grid data, the aneurysm risk assessment of the labeled 3D triangle grid data is obtained through a sigmoid function operation As a result, the aneurysm risk assessment model is obtained.
  5. 如权利要求4所述的方法,其特征在于,所述基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,具体包括:The method according to claim 4, wherein the obtaining a one-dimensional vector of the marked 3D triangular mesh data through an arithmetic block operation based on the marked 3D triangular mesh data specifically includes :
    将所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块以及(64,128,1024)维度的第二运算块的运算,获得(N,1024)的矩阵,其中N为所述已标注的3D三角形网格数据中属于顶点的个数;The labeled 3D triangle mesh data is subjected to the operations of the first operation block of (64, 64) dimensions and the second operation block of (64, 128, 1024) dimensions to obtain a matrix of (N, 1024), where N is The number of vertices in the labeled 3D triangle mesh data;
    将所述(N,1024)的矩阵经过最大值池化层降维,获得所述已标注的3D三角形网格数据的一维向量。The dimension reduction of the matrix of (N, 1024) through the maximum pooling layer, to obtain the one-dimensional vector of the labeled 3D triangular mesh data.
  6. 如权利要求4所述的方法,其特征在于,所述基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,具体包括:The method according to claim 4, wherein the second data and/or third data corresponding to the one-dimensional vector and the labeled 3D triangle mesh data are obtained through a sigmoid function operation. The aneurysm risk assessment result of the marked 3D triangle grid data specifically includes:
    将所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据输入维度为(512,256,2)的多层感知机,获得第四数据;Input the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle grid data into a multi-layer perceptron with a dimension of (512, 256, 2) to obtain fourth data;
    将所述第四数据经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤破裂风险概率;Subjecting the fourth data to a sigmoid function operation to obtain the aneurysm rupture risk probability of the marked 3D triangle mesh data;
    和/或and / or
    将所述一维向量及所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块获得的新的矩阵,经过维度为(512,256,128)的第三运算块和维度为(128,2)的第四运算块,获得(N,2)的矩阵;A new matrix obtained by passing the one-dimensional vector and the labeled 3D triangle mesh data through the first operation block of (64, 64) dimensions, passing through the third operation block of dimensions (512, 256, 128), and the dimension is ( 128,2) the fourth operation block to obtain the matrix of (N,2);
    将所述(N,2)的矩阵经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的顶点分类。The matrix of (N, 2) is subjected to a sigmoid function operation to obtain the vertex classification of the labeled 3D triangle mesh data.
  7. 一种动脉瘤破裂风险评估系统,其特征在于,所述评估系统包括:An aneurysm rupture risk assessment system, characterized in that the assessment system includes:
    输入模块,获取待处理数据,其中,所述待处理数据包括第一数据和/或第二数据和/或第三数据;An input module to obtain data to be processed, where the data to be processed includes first data and/or second data and/or third data;
    评估模块,将所述第一数据转换为3D三角形网格数据,基于所述3D三角形网格数据和/或所述第二数据和/或所述第三数据输入动脉瘤风险评估模型,获得所述待处理数据的动脉瘤风险评估结果,其中,所述动脉瘤风险评估模型是基于已标注的3D三角形网格数据和/或所述第二数据和/或所述第三数据训练获得的数字模型,所述动脉瘤风险评估结果包括所述待处理数据的动脉瘤破裂风险概率和/或所述待处理数据对应的3D三角形网格数据的顶点分类,所述待处理数据对应的3D三角形网格数据的顶点分类为所述待处理数据对应的3D三角形网格数据的顶点是属于血管或者属于动脉瘤;The evaluation module converts the first data into 3D triangle mesh data, and inputs the aneurysm risk assessment model based on the 3D triangle mesh data and/or the second data and/or the third data to obtain all The aneurysm risk assessment result of the to-be-processed data, wherein the aneurysm risk assessment model is based on the marked 3D triangle grid data and/or the second data and/or the number obtained by training on the third data Model, the aneurysm risk assessment result includes the aneurysm rupture risk probability of the data to be processed and/or the vertex classification of the 3D triangle mesh data corresponding to the data to be processed, and the 3D triangle mesh corresponding to the data to be processed The vertex of the grid data is classified as whether the vertex of the 3D triangle mesh data corresponding to the data to be processed belongs to a blood vessel or an aneurysm;
    输出模块,输出所述待处理数据的动脉瘤风险评估结果。The output module outputs the aneurysm risk assessment result of the to-be-processed data.
  8. 如权利要求7所述的系统,其特征在于,所述已标注的3D三角形网格数据是基于已知的3D三角形网格数据,将所述已知的3D三角形网格数据中属于血管的顶点标注为0,将所述已知的3D三角形网格数据中属于动脉瘤的顶点标注为1。The system of claim 7, wherein the labeled 3D triangle mesh data is based on known 3D triangle mesh data, and the vertices belonging to the blood vessel in the known 3D triangle mesh data It is marked as 0, and the vertex belonging to the aneurysm in the known 3D triangle mesh data is marked as 1.
  9. 如权利要求8所述的系统,其特征在于,所述已标注的3D三角形网格数据进一步包括:The system of claim 8, wherein the labeled 3D triangle mesh data further comprises:
    对所述已标注的3D三角形网格数据按照预设比率随机舍弃属于血管的顶点和/或属于动脉瘤的顶点,以获得属于动脉瘤点周围的顶点坐标、以及所述属于动脉瘤点周围的顶点之间相连的边信息。The vertices belonging to the blood vessel and/or the vertices belonging to the aneurysm are randomly discarded for the labeled 3D triangle mesh data according to the preset ratio to obtain the coordinates of the vertices belonging to the aneurysm point and the points belonging to the aneurysm. Edge information connected between vertices.
  10. 如权利要求7所述的系统,其特征在于,所述动脉瘤风险评估模型的训练包括:The system of claim 7, wherein the training of the aneurysm risk assessment model comprises:
    基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,其中,所述一维向量表示所述已标注的3D三角形网格数据的全局特征向量;Based on the labeled 3D triangle mesh data, a one-dimensional vector of the labeled 3D triangle mesh data is obtained through an arithmetic block operation, where the one-dimensional vector represents the labeled 3D triangle mesh The global feature vector of the data;
    基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,以获得所述动脉瘤风险评估模型。Based on the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle grid data, the aneurysm risk assessment of the labeled 3D triangle grid data is obtained through a sigmoid function operation As a result, the aneurysm risk assessment model is obtained.
  11. 如权利要求10所述的系统,其特征在于,所述基于所述已标注的3D三角形网格数据,经过运算块运算,获得所述已标注的3D三角形网格数据的一维向量,具体包括:The system according to claim 10, wherein the one-dimensional vector of the marked 3D triangular mesh data is obtained through an arithmetic block operation based on the marked 3D triangular mesh data, which specifically includes :
    将所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块以及(64,128,1024)维度的第二运算块的运算,获得(N,1024)的矩阵,其中N为所述已标注的3D三角形网格数据中属于顶点的个数;The labeled 3D triangle mesh data is subjected to the operations of the first operation block of (64, 64) dimensions and the second operation block of (64, 128, 1024) dimensions to obtain a matrix of (N, 1024), where N is The number of vertices in the labeled 3D triangle mesh data;
    将所述(N,1024)的矩阵经过最大值池化层降维,获得所述已标注的3D三角形网格数据的一维向量。The dimension reduction of the matrix of (N, 1024) through the maximum pooling layer, to obtain the one-dimensional vector of the labeled 3D triangular mesh data.
  12. 如权利要求10所述的系统,其特征在于,所述基于所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据,经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤风险评估结果,具体包括:The system according to claim 10, wherein the second data and/or third data corresponding to the one-dimensional vector and the labeled 3D triangle mesh data are obtained through a sigmoid function operation. The aneurysm risk assessment result of the marked 3D triangle grid data specifically includes:
    将所述一维向量及所述已标注的3D三角形网格数据对应的第二数据和/或第三数据输入维度为(512,256,2)的多层感知机,获得第四数据;Input the one-dimensional vector and the second data and/or third data corresponding to the labeled 3D triangle grid data into a multi-layer perceptron with a dimension of (512, 256, 2) to obtain fourth data;
    将所述第四数据经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的动脉瘤破裂风险概率;Subjecting the fourth data to a sigmoid function operation to obtain the aneurysm rupture risk probability of the marked 3D triangle mesh data;
    和/或and / or
    将所述一维向量及所述已标注的3D三角形网格数据经过(64,64)维度的第一运算块获得的新的矩阵,经过维度为(512,256,128)的第三运算块和维度为(128,2)的第四运算块,获得(N,2)的矩阵;A new matrix obtained by passing the one-dimensional vector and the labeled 3D triangle mesh data through the first operation block of (64, 64) dimensions, passing through the third operation block of dimensions (512, 256, 128), and the dimension is ( 128,2) the fourth operation block to obtain the matrix of (N,2);
    将所述(N,2)的矩阵经过sigmoid函数运算,获得所述已标注的3D三角形网格数据的顶点分类。The matrix of (N, 2) is subjected to a sigmoid function operation to obtain the vertex classification of the labeled 3D triangle mesh data.
PCT/CN2020/130058 2019-11-22 2020-11-19 Method and system for assessing aneurysm rupture risk WO2021098768A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201911154109.XA CN111081378B (en) 2019-11-22 2019-11-22 Aneurysm rupture risk assessment method and system
CN201911154109.X 2019-11-22

Publications (1)

Publication Number Publication Date
WO2021098768A1 true WO2021098768A1 (en) 2021-05-27

Family

ID=70311265

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/130058 WO2021098768A1 (en) 2019-11-22 2020-11-19 Method and system for assessing aneurysm rupture risk

Country Status (2)

Country Link
CN (1) CN111081378B (en)
WO (1) WO2021098768A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113744883A (en) * 2021-09-22 2021-12-03 皖南医学院第一附属医院(皖南医学院弋矶山医院) Construction method and device for predicting intracranial aneurysm rupture model

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111081378B (en) * 2019-11-22 2022-05-20 强联智创(北京)科技有限公司 Aneurysm rupture risk assessment method and system
CN113130078B (en) * 2021-05-11 2022-09-23 首都医科大学附属北京天坛医院 Method, device and equipment for predicting intracranial aneurysm occlusion
CN113130030B (en) * 2021-05-11 2022-09-23 首都医科大学附属北京天坛医院 Method, device and equipment for evaluating stability of intracranial aneurysm

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107273658A (en) * 2017-05-16 2017-10-20 哈尔滨医科大学 Rupture of intracranial aneurysm risk is estimated and its device that image is classified
CN109907732A (en) * 2019-04-09 2019-06-21 广州新脉科技有限公司 A kind of appraisal procedure and system of rupture of intracranial aneurysm risk
CN109961850A (en) * 2019-03-19 2019-07-02 肖仁德 A kind of method, apparatus, computer equipment for assessing rupture of intracranial aneurysm risk
CN110517780A (en) * 2019-09-02 2019-11-29 强联智创(北京)科技有限公司 A kind of aneurysm rupture methods of risk assessment and system
CN110534193A (en) * 2019-09-02 2019-12-03 强联智创(北京)科技有限公司 A kind of aneurysm rupture methods of risk assessment and system
CN111081378A (en) * 2019-11-22 2020-04-28 强联智创(北京)科技有限公司 Aneurysm rupture risk assessment method and system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102007026519A1 (en) * 2007-06-08 2008-12-18 Siemens Ag Method for determination of rupture risk of aneurysm of patient, involves determining rupture risk at side of calculation device depending on specific person related factor for patients, and anatomy related factor is also related

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107273658A (en) * 2017-05-16 2017-10-20 哈尔滨医科大学 Rupture of intracranial aneurysm risk is estimated and its device that image is classified
CN109961850A (en) * 2019-03-19 2019-07-02 肖仁德 A kind of method, apparatus, computer equipment for assessing rupture of intracranial aneurysm risk
CN109907732A (en) * 2019-04-09 2019-06-21 广州新脉科技有限公司 A kind of appraisal procedure and system of rupture of intracranial aneurysm risk
CN110517780A (en) * 2019-09-02 2019-11-29 强联智创(北京)科技有限公司 A kind of aneurysm rupture methods of risk assessment and system
CN110534193A (en) * 2019-09-02 2019-12-03 强联智创(北京)科技有限公司 A kind of aneurysm rupture methods of risk assessment and system
CN111081378A (en) * 2019-11-22 2020-04-28 强联智创(北京)科技有限公司 Aneurysm rupture risk assessment method and system

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113744883A (en) * 2021-09-22 2021-12-03 皖南医学院第一附属医院(皖南医学院弋矶山医院) Construction method and device for predicting intracranial aneurysm rupture model

Also Published As

Publication number Publication date
CN111081378A (en) 2020-04-28
CN111081378B (en) 2022-05-20

Similar Documents

Publication Publication Date Title
WO2021098768A1 (en) Method and system for assessing aneurysm rupture risk
US11875503B2 (en) Method and system for measuring morphological parameters of an intracranial aneurysm image
WO2020083374A1 (en) Method and system for measuring morphological parameters of an intracranial aneurysm image
CN109448004B (en) Centerline-based intracranial blood vessel image interception method and system
CN110517780A (en) A kind of aneurysm rupture methods of risk assessment and system
CN109584997B (en) Method and system for measuring morphological parameters of intracranial aneurysm image
Zeng et al. Automatic diagnosis based on spatial information fusion feature for intracranial aneurysm
CN109389637B (en) Method and system for measuring morphological parameters of intracranial aneurysm image
US20230139405A1 (en) Stenosis assessment method and device based on intracranial DSA imaging
CN110534193A (en) A kind of aneurysm rupture methods of risk assessment and system
CN109472823A (en) A kind of measurement method and system of the Morphologic Parameters of intracranial aneurysm image
CN109472780A (en) A kind of measurement method and system of the Morphologic Parameters of intracranial aneurysm image
CN110517242A (en) A kind of aneurysmal analysis method and device
CN109671066A (en) A kind of method and system of the cerebral infarction judgement based on head CT images
JP2017189308A (en) Medical image processing apparatus, medical image diagnostic apparatus and program
Eulzer et al. Vessel Maps: A Survey of Map‐Like Visualizations of the Cardiovascular System
CN109447967B (en) Method and system for segmenting intracranial aneurysm image
CN111223089B (en) Aneurysm detection method and device and computer readable storage medium
Spiegel et al. A 2D driven 3D vessel segmentation algorithm for 3D digital subtraction angiography data
Karmonik et al. Magnetic resonance imaging as a tool to assess reliability in simulating hemodynamics in cerebral aneurysms with a dedicated computational fluid dynamics prototype: preliminary results
CN109741339A (en) A kind of partition method and system
CN109377504A (en) A kind of entocranial artery blood-vessel image dividing method and system
CN109584261B (en) Method and system for segmenting intracranial aneurysm image
CN110739078B (en) Aneurysm rupture risk assessment method and system
CN112734726A (en) Typing method, device and equipment for angiography

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20890375

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20890375

Country of ref document: EP

Kind code of ref document: A1