WO2023280148A1 - Blood vessel segmentation method and apparatus, and electronic device and readable medium - Google Patents

Blood vessel segmentation method and apparatus, and electronic device and readable medium Download PDF

Info

Publication number
WO2023280148A1
WO2023280148A1 PCT/CN2022/103838 CN2022103838W WO2023280148A1 WO 2023280148 A1 WO2023280148 A1 WO 2023280148A1 CN 2022103838 W CN2022103838 W CN 2022103838W WO 2023280148 A1 WO2023280148 A1 WO 2023280148A1
Authority
WO
WIPO (PCT)
Prior art keywords
blood vessel
convolutional network
features
segmentation
scale
Prior art date
Application number
PCT/CN2022/103838
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 WO2023280148A1 publication Critical patent/WO2023280148A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20076Probabilistic image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Definitions

  • the present invention claims the priority of the application proposed by the applicant, the application date is July 7, 2021, the application number is CN2021107686178, and the application name is "a blood vessel segmentation method and device".
  • the entire content of the above application is hereby incorporated by reference in its entirety.
  • the invention belongs to the technical field of medical image blood vessel segmentation, and in particular relates to a blood vessel segmentation method, device, electronic equipment and readable medium.
  • the segmentation algorithm can realize automatic reconstruction of blood vessels (such as head and neck vessels, coronary arteries, etc.), which greatly improves the operating efficiency of the hospital while reducing the work pressure of technicians.
  • blood vessels such as head and neck vessels, coronary arteries, etc.
  • some external factors such as artifacts, noise, shooting technology, etc. will affect the quality of blood vessel imaging, so that the results of the segmentation algorithm may be missed or broken.
  • the existing blood vessel segmentation method mainly predicts the blood vessel area through the extraction of local features, and cannot model the overall structure of the blood vessel; the existing method mainly uses multi-scale local features for blood vessel prediction, and cannot predict the global structure of multiple scales. Features are learned and combined.
  • the present invention provides a blood vessel segmentation method, device, electronic equipment and readable medium.
  • the present invention adopts the following technical solutions.
  • the present invention provides a blood vessel segmentation method, including:
  • the pre-segmentation results and segmentation labels are used as input to construct a blood vessel distribution map.
  • the nodes of the graph are pixel areas with a high probability of blood vessels and consistent gray levels.
  • the shape of the nodes is consistent with the direction of the blood vessels;
  • the edges of the graph represent the correlation of connected nodes.
  • the length of the edge is less than the set threshold;
  • the CTA image is input to the second convolutional network with the same network structure and weight parameters as the first convolutional network, and the output of the second convolutional network is multiplied by the input CTA image and then input to the third convolutional network;
  • the third convolutional network and the graph convolutional network perform feature interaction based on bidirectional mapping, and perform global feature modeling based on the vascular distribution map to capture multi-scale local and global features, and realize vascular region prediction;
  • the method also includes the step of normalizing the input CTA image, normalizing the grayscale of the pixel to [0,255], and the calculation formula is as follows:
  • x min and x max are the minimum and maximum values of the gray value of the pixel before normalization, respectively.
  • the blood vessel pre-segmentation method specifically includes:
  • the multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size through upsampling or deconvolution, and combines them with the encoding features at the corresponding scale to obtain the segmentation result.
  • the method for predicting the vascular region specifically includes:
  • the third convolutional network projects multiple pixels of encoded features into a single node in the graph convolutional network through forward mapping
  • the graph convolutional network propagates the features on the blood vessel distribution map for global feature modeling
  • the graph convolutional network maps the node features back to the feature space of the third convolutional network through reverse mapping, and propagates the global features to each part for feature enhancement;
  • the method uses the Dice coefficient to measure the similarity between the blood vessel segmentation result and the gold standard marked by the doctor, and the formula is as follows:
  • a and B are the blood vessel pixel point set of the gold standard and the pixel point set of the blood vessel segmentation result respectively, and
  • are A, B and the intersection of A and B respectively The number of pixels.
  • the present invention provides a blood vessel segmentation device, comprising:
  • the pre-segmentation module is used to input the CTA image to the first convolutional network for pre-segmentation of blood vessels based on multi-scale feature extraction;
  • the distribution map building module is used to construct a blood vessel distribution map with pre-segmentation results and segmentation labels as input.
  • the nodes of the map are pixel areas with a high probability of blood vessels and consistent gray levels.
  • the shape of the nodes is consistent with the direction of the blood vessels;
  • the edges of the graph represent The relevance of the connected nodes, the length of the edge is less than the set threshold;
  • the first-stage segmentation module is used to input the CTA image to the second convolutional network with the same network structure and weight parameters as the first convolutional network, and to multiply the output of the second convolutional network with the input CTA image and then input it to The third convolutional network;
  • the second stage segmentation module is used for the third convolutional network and the graph convolutional network to perform feature interaction based on two-way mapping, and to perform global feature modeling based on the blood vessel distribution map to capture multi-scale local and global features, and realize blood vessel region prediction;
  • the feature fusion module is used to predict the blood vessel segmentation result by fusing multi-scale features.
  • the device also includes a normalization module for preprocessing the input CTA image, which is used to normalize the grayscale of the pixels to [0, 255], and the calculation formula is as follows:
  • x min and x max are the minimum and maximum values of the gray value of the pixel before normalization, respectively.
  • the blood vessel pre-segmentation method specifically includes:
  • the multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size through upsampling or deconvolution, and combines them with the encoding features at the corresponding scale to obtain the segmentation result.
  • the method for predicting the vascular region specifically includes:
  • the third convolutional network projects multiple pixels of encoded features into a single node in the graph convolutional network through forward mapping
  • the graph convolutional network propagates the features on the blood vessel distribution map for global feature modeling
  • the graph convolutional network maps the node features back to the feature space of the third convolutional network through reverse mapping, and propagates the global features to each part for feature enhancement;
  • the method uses the Dice coefficient to measure the similarity between the blood vessel segmentation result and the gold standard marked by the doctor, and the formula is as follows:
  • a and B are the blood vessel pixel point set of the gold standard and the pixel point set of the blood vessel segmentation result respectively, and
  • are A, B and the intersection of A and B respectively The number of pixels.
  • the present application provides a readable medium, including execution instructions, and when a processor of an electronic device executes the execution instructions, the electronic device executes the method described in any one of the first aspects.
  • the present application provides an electronic device, including a processor and a memory storing execution instructions, and when the processor executes the execution instructions stored in the memory, the processor executes the method described in the first aspect any of the methods described.
  • the present invention has the following beneficial effects:
  • the present invention models the global characteristics of blood vessels based on the construction of a blood vessel distribution network, uses the overall structure of blood vessels to enhance the features of segmented regions, and taps the potential direction of blood vessels.
  • the present invention performs vessel segmentation based on cross-network multi-scale feature fusion and fusion of local features and global features.
  • feature representation with stronger information can be obtained , to effectively prevent the model from predicting segmentation regions that do not conform to the vascular structure, thereby improving the accuracy of vascular segmentation.
  • FIG. 1 is a flowchart of a blood vessel segmentation method according to an embodiment of the present invention.
  • FIG. 2 is a schematic diagram of a network structure according to an embodiment of the present invention.
  • Fig. 3 is a block diagram of a blood vessel segmentation device according to an embodiment of the present invention.
  • FIG. 4 is a schematic structural diagram of an electronic device provided by an embodiment of the present invention.
  • Fig. 1 is a flowchart of a blood vessel segmentation method according to an embodiment of the present invention, including the following steps:
  • Step 101 input the CTA image to the first convolutional network, and perform blood vessel pre-segmentation based on multi-scale feature extraction;
  • Step 102 using the pre-segmentation result and the segmentation label as input to construct a blood vessel distribution map, the nodes of the map are pixel point areas with a high probability of blood vessels and consistent gray levels, and the shape of the nodes is consistent with the direction of the blood vessels; the edges of the graph represent the connected nodes. Relevance, the length of the edge is less than the set threshold;
  • Step 103 the CTA image is input to the second convolutional network with the same network structure and weight parameters as the first convolutional network, and the output of the second convolutional network is multiplied by the input CTA image and then input to the third convolutional network ;
  • Step 104 the third convolutional network and the graph convolutional network perform feature interaction based on bidirectional mapping, and perform global feature modeling based on the blood vessel distribution map to capture multi-scale local and global features, and realize blood vessel region prediction;
  • step 105 the vessel segmentation result is predicted by fusing the multi-scale features.
  • step 101 is mainly used for blood vessel pre-segmentation.
  • the blood vessel segmentation algorithm in this embodiment is mainly implemented by a convolutional neural network (CNN), including three convolutional networks (first to third convolutional networks) and a graph convolutional network, as shown in FIG. 2 .
  • CNN is a feedforward neural network, but unlike the general fully connected feedforward neural network, its convolutional layer has the characteristics of local connection and weight sharing, so it can greatly reduce the number of weight parameters, thereby reducing Model complexity and increased running speed.
  • a typical CNN is composed of convolutional layers, pooling layers (or pooling layers, downsampling layers), and fully connected layers cross-stacked.
  • the function of the convolution layer is to extract the features of a local area through the convolution operation of the convolution kernel and the input image. Different convolution kernels are equivalent to different feature extractors.
  • the role of the pooling layer is to perform feature selection and reduce the number of features, thereby further reducing the number of parameters. Generally, the maximum pooling method and the average pooling method are used.
  • the fully connected layer is used to fuse the different features obtained.
  • the CTA image is input to the first convolution network, and multiple convolution kernels of different sizes are used for multi-scale feature extraction, so as to pre-segment blood vessels.
  • step 102 is mainly used to construct a blood vessel distribution map.
  • the blood vessel pre-segmentation result obtained in the previous step (the probability that a pixel belongs to a blood vessel) and the segmentation label (the label marked on the image by a senior doctor is the gold standard) are used as input to construct a blood vessel distribution map.
  • the distribution graph consists of nodes and edges, as shown in Figure 2. The nodes in the figure have the following characteristics: first, the probability of blood vessels in the pixels in the nodes is high; second, the gray levels of the pixels in the nodes are consistent; third, the shape of the nodes conforms to the direction of blood vessels.
  • An edge in the graph is connected between two nodes, indicating the relevance of the connected nodes, and the length of the edge is less than the set threshold.
  • a pre-segmentation model is used to construct a blood vessel distribution map, which can model or fit the real structure of blood vessels, and is beneficial to improve the accuracy of blood vessel segmentation.
  • step 103 is mainly used for the first stage of blood vessel segmentation.
  • blood vessel segmentation is performed based on cross-network multi-scale feature fusion.
  • blood vessel segmentation consists of two stages: the first stage uses the second convolutional network (UNET-2); the second stage uses the third convolutional network (UNET-3) and graph convolutional network (UNET- G).
  • the network structure of the second convolutional network is the same as that of the pre-segmented first convolutional network, and is initialized with the weights of the pre-segmented model of the first convolutional network. Input the CTA image to the second convolutional network to obtain the prediction result of the first stage.
  • the prediction result of the first stage is multiplied by the original image, it is input to the third convolutional network of the second stage.
  • the prediction result of the first stage is not directly used as the input of the third convolutional network, but multiplied by the original image and then used as the input of the third convolutional network, in order to prevent loss of information of the original image.
  • step 104 is mainly used for the second-stage blood vessel segmentation of the cross-network model.
  • the second stage includes a third convolutional network and a graph convolutional network.
  • the graph convolutional network takes the blood vessel distribution map constructed in step 102 as input to perform global feature modeling, and performs two-way propagation and fusion of local features and global features through bidirectional (forward and reverse) mapping with the third convolutional network, In this way, the prediction of blood vessel area is realized.
  • step 105 is mainly used to predict the blood vessel segmentation result.
  • the network structure and weight parameters of the second convolutional network in the first stage are the same as those of the pre-segmented first convolutional network, so the second convolutional network also performs multi-scale feature extraction. Therefore, after the fusion of local features and global features, multi-scale features are fused to obtain the final vessel segmentation.
  • the method also includes the step of normalizing the input CTA image, normalizing the grayscale of the pixel to [0,255], the calculation formula is as follows:
  • x min and x max are the minimum and maximum values of the gray value of the pixel before normalization, respectively.
  • This embodiment provides an image preprocessing method.
  • normalization processing is performed on the CTA image before it is input, and the grayscale of the pixel is normalized to [0, 255], and the normalization formula is as the above formula.
  • the blood vessel pre-segmentation method specifically includes:
  • the multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size through upsampling or deconvolution, and combines them with the encoding features at the corresponding scale to obtain the segmentation result.
  • the first convolutional network (the same for the other two convolutional networks and the graph convolutional network) includes an encoding layer and a decoding layer.
  • the CTA sequence is first input to the encoding layer, and the encoding layer uses different convolution kernels to convolve the image.
  • Multi-scale feature extraction is realized after operation and pooling operation.
  • the size of extracted feature spaces of different scales constitutes a geometric series with a common ratio of 1/2, the space size of the i-th scale is 1/2 i of the original space size, and the number N of different scales is generally 4-5.
  • the extracted features contain the shape and spatial location information of vessel segmentation.
  • the multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size, and combines them with the encoding features at the corresponding scales (feature cascading or splicing) to obtain segmentation results.
  • the decoding layer maps the global features back to the original feature size through upsampling (interpolation) or deconvolution.
  • the method for predicting a blood vessel region specifically includes:
  • the third convolutional network projects multiple pixels of encoded features into a single node in the graph convolutional network through forward mapping
  • the graph convolutional network propagates the features on the blood vessel distribution map for global feature modeling
  • the graph convolutional network maps the node features back to the feature space of the third convolutional network through reverse mapping, and propagates the global features to each part for feature enhancement;
  • Cross-network mainly refers to feature interaction between the third convolutional network and the graph convolutional network through bidirectional mapping.
  • the third convolutional network performs forward mapping, and projects multiple pixels of encoded features into a single node in the graph convolutional network;
  • the graph convolutional network performs global feature modeling based on the blood vessel distribution map, and performs inverse To map the node features back to the feature space of the third convolutional network, propagate the global features to each local feature for feature enhancement; finally, capture the local features of the image through the third convolutional network, and use the graph convolution
  • the network captures the global features of blood vessels, and uses feature interaction to combine local features and global features (feature cascade) to realize the prediction of the area where blood vessels are located.
  • the method uses the Dice coefficient to measure the similarity between the blood vessel segmentation result and the gold standard marked by the doctor, and the formula is as follows:
  • a and B are the blood vessel pixel point set of the gold standard and the pixel point set of the blood vessel segmentation result respectively, and
  • are A, B and the intersection of A and B respectively The number of pixels.
  • the Dice coefficient is used to represent the similarity between the blood vessel segmentation result and the gold standard marked by doctors, and the formula is as above.
  • the numerator in the formula is twice the number of pixels where the segmentation result coincides with the gold standard, and the denominator is the sum of the pixel points of the segmentation result and the gold standard. According to the formula, the more pixels that overlap the segmentation result and the gold standard, the greater the similarity; when the two completely overlap, the maximum similarity is 1.
  • Fig. 3 is a schematic diagram of the composition of a blood vessel segmentation device according to an embodiment of the present invention, and the device includes:
  • the pre-segmentation module 11 is used to input the CTA image to the first convolutional network for pre-segmentation of blood vessels based on multi-scale feature extraction;
  • the distribution map construction module 12 is used to construct a blood vessel distribution map with pre-segmentation results and segmentation labels as input.
  • the nodes of the map are pixel areas with high blood vessel probability and consistent gray scale, and the shape of the nodes is consistent with the direction of the blood vessels; Indicates the relevance of the connected nodes, and the length of the edge is less than the set threshold;
  • the first stage segmentation module 13 is used to input the CTA image to the second convolutional network with the same network structure and weight parameters as the first convolutional network, and input the output after multiplying the output of the second convolutional network with the input CTA image to the third convolutional network;
  • the second-stage segmentation module 14 is used for feature interaction between the third convolutional network and the graph convolutional network based on bidirectional mapping, and global feature modeling based on the blood vessel distribution map to capture multi-scale local and global features and realize blood vessel region prediction;
  • the feature fusion module 15 is configured to predict the blood vessel segmentation result by fusing multi-scale features.
  • the device of this embodiment can be used to implement the technical solution of the method embodiment shown in FIG. 1 , and its implementation principle and technical effect are similar, and will not be repeated here. The same is true for the following embodiments, which will not be further described.
  • the device also includes a normalization module for preprocessing the input CTA image, which is used to normalize the grayscale of the pixels to [0,255], and the calculation formula is as follows:
  • x min and x max are the minimum and maximum values of the gray value of the pixel before normalization, respectively.
  • the blood vessel pre-segmentation method specifically includes:
  • the multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size through upsampling or deconvolution, and combines them with the encoding features at the corresponding scale to obtain the segmentation result.
  • the method for predicting a blood vessel region specifically includes:
  • the third convolutional network projects multiple pixels of encoded features into a single node in the graph convolutional network through forward mapping
  • the graph convolutional network propagates the features on the blood vessel distribution map for global feature modeling
  • the graph convolutional network maps the node features back to the feature space of the third convolutional network through reverse mapping, and propagates the global features to each part for feature enhancement;
  • the method uses the Dice coefficient to measure the similarity between the blood vessel segmentation result and the gold standard marked by the doctor, and the formula is as follows:
  • a and B are the blood vessel pixel point set of the gold standard and the pixel point set of the blood vessel segmentation result respectively, and
  • are A, B and the intersection of A and B respectively The number of pixels.
  • FIG. 4 is a schematic structural diagram of an electronic device provided by an embodiment of the present application.
  • the electronic device includes a processor, and optionally also includes an internal bus, a network interface, and a memory.
  • the memory may include a memory, such as a high-speed random-access memory (Random-Access Memory, RAM), and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.
  • RAM random-Access Memory
  • non-volatile memory such as at least one disk memory.
  • the electronic device may also include hardware required by other services.
  • the processor, the network interface and the memory can be connected to each other through an internal bus, which can be an ISA (Industry Standard Architecture, industry standard architecture) bus, a PCI (Peripheral Component Interconnect, peripheral component interconnection standard) bus or an EISA (Extended Industry Standard Architecture, extended industry standard architecture) bus, etc.
  • the bus can be divided into address bus, data bus, control bus and so on. For ease of representation, only one double-headed arrow is used in FIG. 4 , but it does not mean that there is only one bus or one type of bus.
  • Memory for storing execution instructions. Specifically, a computer program that can be executed by executing instructions.
  • the memory which can include internal memory and non-volatile memory, provides instructions and data to the processor for execution.
  • the processor reads the corresponding execution instructions from the non-volatile memory into the memory and then runs them. It can also obtain the corresponding execution instructions from other devices to form blood vessel segmentation at the logical level. device.
  • the processor executes the execution instructions stored in the memory, so as to implement the blood vessel segmentation method provided in any embodiment of the present application through the executed execution instructions.
  • the above-mentioned method performed by the blood vessel segmentation apparatus provided in the embodiment shown in FIG. 1 of the present application may be applied to a processor or implemented by the processor.
  • a processor may be an integrated circuit chip with signal processing capabilities.
  • each step of the above method can be completed by an integrated logic circuit of hardware in a processor or an instruction in the form of software.
  • the above-mentioned processor can be a general-purpose processor, including a central processing unit (Central Processing Unit, CPU), a network processor (Network Processor, NP), etc.; it can also be a digital signal processor (Digital Signal Processor, DSP), a dedicated integrated Circuit (Application Specific Integrated Circuit, ASIC), Field-Programmable Gate Array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.
  • DSP Digital Signal Processor
  • ASIC Application Specific Integrated Circuit
  • FPGA Field-Programmable Gate Array
  • a general-purpose processor may be a microprocessor, or the processor may be any conventional processor, or the like.
  • the steps of the method disclosed in connection with the embodiments of the present application may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor.
  • the software module can be located in a mature storage medium in the field such as random access memory, flash memory, read-only memory, programmable read-only memory or electrically erasable programmable memory, register.
  • the storage medium is located in the memory, and the processor reads the information in the memory, and completes the steps of the above method in combination with its hardware.
  • the embodiment of the present application also proposes a readable medium, the readable storage medium stores execution instructions, and when the stored execution instructions are executed by the processor of the electronic device, the electronic device can execute the electronic device provided in any embodiment of the present application.
  • the blood vessel segmentation method and is specifically used to implement the above blood vessel segmentation method.
  • the electronic equipment described in the foregoing embodiments may be a computer.
  • each embodiment in the present application is described in a progressive manner, the same and similar parts of each embodiment can be referred to each other, and each embodiment focuses on the differences from other embodiments.
  • the description is relatively simple, and for relevant parts, please refer to part of the description of the method embodiment.
  • a blood vessel segmentation method, device, electronic equipment, and readable medium provided by the present invention use image recognition technology to build a blood vessel distribution network to model the global characteristics of blood vessels, use the overall structure of blood vessels to enhance the features of the segmented area, and mine blood vessels
  • the potential direction of the structure makes full use of the processing basis of program automation in computer technology, which greatly improves the accuracy of the prediction and segmentation of vessel segmentation regions under the complex imaging quality environment of angiographic images.
  • the formed products can be mass-produced and quickly applied to systems or scenarios with high demand for vascular imaging.

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Computational Linguistics (AREA)
  • Software Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Medical Informatics (AREA)
  • Image Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

The present invention provides a blood vessel segmentation method and apparatus, and an electronic device and a readable medium. The method comprises: inputting a CTA image into a first convolutional network, and performing blood vessel pre-segmentation based on multi-scale feature extraction; constructing a blood vessel distribution graph; inputting the CTA image into a second convolutional network, multiplying an output of the second convolutional network by the input CTA image, and then inputting a multiplied result into a third convolutional network; the third convolutional network and a graph convolutional network performing feature interaction on the basis of bidirectional mapping, performing global feature modeling on the basis of the blood vessel distribution graph to capture multi-scale local and global features, and implementing blood vessel region prediction; and fusing the multi-scale features to predict a blood vessel segmentation result. According to the present invention, blood vessel segmentation is performed on the basis of cross-network multi-scale feature fusion and fusion of local features and global features, and compared with blood vessel segmentation based on local features in the prior art, the present invention can effectively prevent a model from predicting a segmentation region which does not conform to a blood vessel structure, thereby improving the accuracy of blood vessel segmentation.

Description

一种血管分割方法、装置、电子设备和可读介质A blood vessel segmentation method, device, electronic device and readable medium
本发明要求由申请人提出的,申请日为2021年7月7日,申请号为CN2021107686178,名称为“一种血管分割方法及装置”的申请的优先权。上述申请的全部内容通过整体引用结合于此。The present invention claims the priority of the application proposed by the applicant, the application date is July 7, 2021, the application number is CN2021107686178, and the application name is "a blood vessel segmentation method and device". The entire content of the above application is hereby incorporated by reference in its entirety.
技术领域technical field
本发明属于医学影像血管分割技术领域,具体涉及一种血管分割方法、装置、电子设备和可读介质。The invention belongs to the technical field of medical image blood vessel segmentation, and in particular relates to a blood vessel segmentation method, device, electronic equipment and readable medium.
发明背景Background of the invention
常见的血管造影技术已经被广泛应用于临床诊断和治疗中。分割算法可实现自动化的血管(如头颈部血管、冠脉等)重建,在减轻技师工作压力的同时,大幅提高了医院的运行效率。在实际场景中,一些外部因素(如伪影、噪声、拍摄技术等)会影响血管成像的质量,从而使得分割算法的结果存在遗漏或者断裂的情况。Common angiographic techniques have been widely used in clinical diagnosis and treatment. The segmentation algorithm can realize automatic reconstruction of blood vessels (such as head and neck vessels, coronary arteries, etc.), which greatly improves the operating efficiency of the hospital while reducing the work pressure of technicians. In the actual scene, some external factors (such as artifacts, noise, shooting technology, etc.) will affect the quality of blood vessel imaging, so that the results of the segmentation algorithm may be missed or broken.
现有血管分割方法,主要通过对局部特征的提取,来预测血管区域,无法对血管的整体结构进行建模;现有方法主要通过多尺度的局部特征进行血管预测,无法对多个尺度的全局特征进行学习和结合。The existing blood vessel segmentation method mainly predicts the blood vessel area through the extraction of local features, and cannot model the overall structure of the blood vessel; the existing method mainly uses multi-scale local features for blood vessel prediction, and cannot predict the global structure of multiple scales. Features are learned and combined.
发明内容Contents of the invention
为了解决现有技术中存在的上述问题,本发明提供一种血管分割方法、装置、电子设备和可读介质。In order to solve the above problems in the prior art, the present invention provides a blood vessel segmentation method, device, electronic equipment and readable medium.
为了实现上述目的,本发明采用以下技术方案。In order to achieve the above object, the present invention adopts the following technical solutions.
第一方面,本发明提供一种血管分割方法,包括:In a first aspect, the present invention provides a blood vessel segmentation method, including:
将CTA图像输入到第一卷积网络,进行基于多尺度特征提取的血管预分割;Input the CTA image to the first convolutional network for pre-segmentation of blood vessels based on multi-scale feature extraction;
以预分割结果和分割标签为输入构建血管分布图,图的节点为血管概率较高、 灰度一致的像素点区域,节点的形状与血管走向一致;图的边表示所连节点的相关性,边的长度小于设定的阈值;The pre-segmentation results and segmentation labels are used as input to construct a blood vessel distribution map. The nodes of the graph are pixel areas with a high probability of blood vessels and consistent gray levels. The shape of the nodes is consistent with the direction of the blood vessels; the edges of the graph represent the correlation of connected nodes. The length of the edge is less than the set threshold;
将CTA图像输入到网络结构和权重参数与第一卷积网络相同的第二卷积网络,并将第二卷积网络的输出与输入CTA图像相乘后输入到第三卷积网络;The CTA image is input to the second convolutional network with the same network structure and weight parameters as the first convolutional network, and the output of the second convolutional network is multiplied by the input CTA image and then input to the third convolutional network;
第三卷积网络与图卷积网络通过基于双向映射进行特征交互,基于血管分布图进行全局特征建模捕捉多尺度的局部和全局特征,实现血管区域预测;The third convolutional network and the graph convolutional network perform feature interaction based on bidirectional mapping, and perform global feature modeling based on the vascular distribution map to capture multi-scale local and global features, and realize vascular region prediction;
通过对多尺度特征进行融合,预测出血管分割结果。Through the fusion of multi-scale features, the result of blood vessel segmentation is predicted.
进一步地,所述方法还包括对输入CTA图像进行归一化的步骤,将像素点的灰度归一化为[0,255],计算公式如下:Further, the method also includes the step of normalizing the input CTA image, normalizing the grayscale of the pixel to [0,255], and the calculation formula is as follows:
Figure PCTCN2022103838-appb-000001
Figure PCTCN2022103838-appb-000001
式中,
Figure PCTCN2022103838-appb-000002
为任一像素点的灰度值x归一化后的灰度值,x min、x max分别为归一化前像素点灰度值的最小值和最大值。
In the formula,
Figure PCTCN2022103838-appb-000002
is the gray value x of any pixel after normalization, and x min and x max are the minimum and maximum values of the gray value of the pixel before normalization, respectively.
进一步地,血管预分割方法具体包括:Further, the blood vessel pre-segmentation method specifically includes:
将CTA序列输入到第一卷积网络的编码层;Input the CTA sequence into the encoding layer of the first convolutional network;
编码层通过卷积和池化操作进行多尺度特征提取,第i个尺度的空间大小为原始空间大小的1/2 i,i=1,2,…,N,N为尺度个数;所述特征包含血管分割的形态和空间位置信息; The encoding layer performs multi-scale feature extraction through convolution and pooling operations, the space size of the i-th scale is 1/2 i of the original space size, i=1,2,...,N, N is the number of scales; the Features include shape and spatial location information of vessel segmentation;
将编码层提取的多尺度特征输入到解码层,解码层将每个尺度的全局特征通过上采样或反卷积映射回原始特征大小,并与对应尺度下的编码特征进行结合得到分割结果。The multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size through upsampling or deconvolution, and combines them with the encoding features at the corresponding scale to obtain the segmentation result.
进一步地,所述血管区域预测的方法具体包括:Further, the method for predicting the vascular region specifically includes:
第三卷积网络通过前向映射,将编码特征的多个像素点投影为图卷积网络中的单一节点;The third convolutional network projects multiple pixels of encoded features into a single node in the graph convolutional network through forward mapping;
图卷积网络将特征在血管分布图上传播,进行全局特征建模;The graph convolutional network propagates the features on the blood vessel distribution map for global feature modeling;
图卷积网络通过反向映射,将节点特征映射回第三卷积网络的特征空间上,将全局特征传播到每个局部,进行特征增强;The graph convolutional network maps the node features back to the feature space of the third convolutional network through reverse mapping, and propagates the global features to each part for feature enhancement;
将第三卷积网络和图卷积网络的特征进行结合,通过第三卷积网络捕捉影像的局部特征,通过图卷积网络捕捉血管的全局特征,利用特征交互将局部特征和全局特征进行结合,完成血管所在区域的预测。Combine the features of the third convolutional network and the graph convolutional network, capture the local features of the image through the third convolutional network, capture the global features of blood vessels through the graph convolutional network, and use feature interaction to combine local features and global features , to complete the prediction of the area where the blood vessel is located.
进一步地,所述方法采用Dice系数度量血管分割结果与医生标注的金标准的相似度,公式如下:Further, the method uses the Dice coefficient to measure the similarity between the blood vessel segmentation result and the gold standard marked by the doctor, and the formula is as follows:
Figure PCTCN2022103838-appb-000003
Figure PCTCN2022103838-appb-000003
式中,A、B分别为金标准的血管像素点集合和血管分割结果的像素点集合,|A|、|B|和|A∩B|分别为A、B和A与B的交集中的像素点的数量。In the formula, A and B are the blood vessel pixel point set of the gold standard and the pixel point set of the blood vessel segmentation result respectively, and |A|, |B| and |A∩B| are A, B and the intersection of A and B respectively The number of pixels.
第二方面,本发明提供一种血管分割装置,包括:In a second aspect, the present invention provides a blood vessel segmentation device, comprising:
预分割模块,用于将CTA图像输入到第一卷积网络,进行基于多尺度特征提取的血管预分割;The pre-segmentation module is used to input the CTA image to the first convolutional network for pre-segmentation of blood vessels based on multi-scale feature extraction;
分布图构建模块,用于以预分割结果和分割标签为输入构建血管分布图,图的节点为血管概率较高、灰度一致的像素点区域,节点的形状与血管走向一致;图的边表示所连节点的相关性,边的长度小于设定的阈值;The distribution map building module is used to construct a blood vessel distribution map with pre-segmentation results and segmentation labels as input. The nodes of the map are pixel areas with a high probability of blood vessels and consistent gray levels. The shape of the nodes is consistent with the direction of the blood vessels; the edges of the graph represent The relevance of the connected nodes, the length of the edge is less than the set threshold;
第一阶段分割模块,用于将CTA图像输入到网络结构和权重参数与第一卷积网络相同的第二卷积网络,并将第二卷积网络的输出与输入CTA图像相乘后输入到第三卷积网络;The first-stage segmentation module is used to input the CTA image to the second convolutional network with the same network structure and weight parameters as the first convolutional network, and to multiply the output of the second convolutional network with the input CTA image and then input it to The third convolutional network;
第二阶段分割模块,用于第三卷积网络与图卷积网络通过基于双向映射进行特征交互,基于血管分布图进行全局特征建模捕捉多尺度的局部和全局特征,实现血管区域预测;The second stage segmentation module is used for the third convolutional network and the graph convolutional network to perform feature interaction based on two-way mapping, and to perform global feature modeling based on the blood vessel distribution map to capture multi-scale local and global features, and realize blood vessel region prediction;
特征融合模块,用于通过对多尺度特征进行融合,预测出血管分割结果。The feature fusion module is used to predict the blood vessel segmentation result by fusing multi-scale features.
进一步地,所述装置还包括对输入CTA图像进行预处理的归一化模块,用于将像素点的灰度归一化为[0,255],计算公式如下:Further, the device also includes a normalization module for preprocessing the input CTA image, which is used to normalize the grayscale of the pixels to [0, 255], and the calculation formula is as follows:
Figure PCTCN2022103838-appb-000004
Figure PCTCN2022103838-appb-000004
式中,
Figure PCTCN2022103838-appb-000005
为任一像素点的灰度值x归一化后的灰度值,x min、x max分别为归一化 前像素点灰度值的最小值和最大值。
In the formula,
Figure PCTCN2022103838-appb-000005
is the gray value x of any pixel after normalization, and x min and x max are the minimum and maximum values of the gray value of the pixel before normalization, respectively.
进一步地,血管预分割方法具体包括:Further, the blood vessel pre-segmentation method specifically includes:
将CTA序列输入到第一卷积网络的编码层;Input the CTA sequence into the encoding layer of the first convolutional network;
编码层通过卷积和池化操作进行多尺度特征提取,第i个尺度的空间大小为原始空间大小的1/2 i,i=1,2,…,N,N为尺度个数;所述特征包含血管分割的形态和空间位置信息; The encoding layer performs multi-scale feature extraction through convolution and pooling operations, the space size of the i-th scale is 1/2 i of the original space size, i=1,2,...,N, N is the number of scales; the Features include shape and spatial location information of vessel segmentation;
将编码层提取的多尺度特征输入到解码层,解码层将每个尺度的全局特征通过上采样或反卷积映射回原始特征大小,并与对应尺度下的编码特征进行结合得到分割结果。The multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size through upsampling or deconvolution, and combines them with the encoding features at the corresponding scale to obtain the segmentation result.
进一步地,所述血管区域预测的方法具体包括:Further, the method for predicting the vascular region specifically includes:
第三卷积网络通过前向映射,将编码特征的多个像素点投影为图卷积网络中的单一节点;The third convolutional network projects multiple pixels of encoded features into a single node in the graph convolutional network through forward mapping;
图卷积网络将特征在血管分布图上传播,进行全局特征建模;The graph convolutional network propagates the features on the blood vessel distribution map for global feature modeling;
图卷积网络通过反向映射,将节点特征映射回第三卷积网络的特征空间上,将全局特征传播到每个局部,进行特征增强;The graph convolutional network maps the node features back to the feature space of the third convolutional network through reverse mapping, and propagates the global features to each part for feature enhancement;
将第三卷积网络和图卷积网络的特征进行结合,通过第三卷积网络捕捉影像的局部特征,通过图卷积网络捕捉血管的全局特征,利用特征交互将局部特征和全局特征进行结合,完成血管所在区域的预测。Combine the features of the third convolutional network and the graph convolutional network, capture the local features of the image through the third convolutional network, capture the global features of blood vessels through the graph convolutional network, and use feature interaction to combine local features and global features , to complete the prediction of the area where the blood vessel is located.
进一步地,所述方法采用Dice系数度量血管分割结果与医生标注的金标准的相似度,公式如下:Further, the method uses the Dice coefficient to measure the similarity between the blood vessel segmentation result and the gold standard marked by the doctor, and the formula is as follows:
Figure PCTCN2022103838-appb-000006
Figure PCTCN2022103838-appb-000006
式中,A、B分别为金标准的血管像素点集合和血管分割结果的像素点集合,|A|、|B|和|A∩B|分别为A、B和A与B的交集中的像素点的数量。In the formula, A and B are the blood vessel pixel point set of the gold standard and the pixel point set of the blood vessel segmentation result respectively, and |A|, |B| and |A∩B| are A, B and the intersection of A and B respectively The number of pixels.
第三方面,本申请提供了一种可读介质,包括执行指令,当电子设备的处理器执行所述执行指令时,所述电子设备执行如第一方面中任一所述的方法。In a third aspect, the present application provides a readable medium, including execution instructions, and when a processor of an electronic device executes the execution instructions, the electronic device executes the method described in any one of the first aspects.
第四方面,本申请提供了一种电子设备,包括处理器以及存储有执行指令的 存储器,当所述处理器执行所述存储器存储的所述执行指令时,所述处理器执行如第一方面中任一所述的方法。In a fourth aspect, the present application provides an electronic device, including a processor and a memory storing execution instructions, and when the processor executes the execution instructions stored in the memory, the processor executes the method described in the first aspect any of the methods described.
与现有技术相比,本发明具有以下有益效果:Compared with the prior art, the present invention has the following beneficial effects:
本发明基于构建血管分布网络对血管全局特征进行建模,利用血管的整体结构对分割区域的特征进行增强,挖掘血管的潜在走向,相比于现有技术人为设定规则对断裂区域进行连接,更能适用实际中的复杂场景;本发明基于跨网络多尺度特征融合及局部特征和全局特征融合进行血管分割,相比于现有技术基于局部特征进行血管分割,可以获得信息更强的特征表示,有效避免模型预测出不符合血管结构的分割区域,从而提高血管分割的准确度。The present invention models the global characteristics of blood vessels based on the construction of a blood vessel distribution network, uses the overall structure of blood vessels to enhance the features of segmented regions, and taps the potential direction of blood vessels. Compared with the prior art, artificially setting rules to connect fractured regions, It is more applicable to complex scenes in practice; the present invention performs vessel segmentation based on cross-network multi-scale feature fusion and fusion of local features and global features. Compared with the prior art for vessel segmentation based on local features, feature representation with stronger information can be obtained , to effectively prevent the model from predicting segmentation regions that do not conform to the vascular structure, thereby improving the accuracy of vascular segmentation.
附图简要说明Brief description of the drawings
为了更清楚地说明本申请实施例或现有的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present application or existing technical solutions, the following will briefly introduce the accompanying drawings that need to be used in the description of the embodiments or prior art. Obviously, the accompanying drawings in the following description are only For some embodiments described in the application, those of ordinary skill in the art can also obtain other drawings based on these drawings without any creative effort.
图1为本发明实施例一种血管分割方法的流程图。FIG. 1 is a flowchart of a blood vessel segmentation method according to an embodiment of the present invention.
图2为本发明实施例的网络结构示意图。FIG. 2 is a schematic diagram of a network structure according to an embodiment of the present invention.
图3为本发明实施例一种血管分割装置的方框图。Fig. 3 is a block diagram of a blood vessel segmentation device according to an embodiment of the present invention.
图4为本发明实施例提供的一种电子设备的结构示意图。FIG. 4 is a schematic structural diagram of an electronic device provided by an embodiment of the present invention.
实施本发明的方式Modes of Carrying Out the Invention
为使本发明的技术目的、技术方案和有益效果更加清楚,下面结合附图对本发明的具体实施方式进行清楚、完整地描述,所描述的具体实施方式只是本发明的一部分实施例,而不是全部的实施例,基于本发明的具体实施方式,本领域技术人员在没有做出创造性劳动的前提下所获得的其他实施例,都属于本发明的保护范围。In order to make the technical objectives, technical solutions and beneficial effects of the present invention clearer, the specific embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings. The described specific embodiments are only part of the embodiments of the present invention, not all of them. Based on the specific implementation manners of the present invention, other embodiments obtained by those skilled in the art without creative efforts all belong to the protection scope of the present invention.
图1为本发明实施例一种血管分割方法的流程图,包括以下步骤:Fig. 1 is a flowchart of a blood vessel segmentation method according to an embodiment of the present invention, including the following steps:
步骤101,将CTA图像输入到第一卷积网络,进行基于多尺度特征提取的血管预分割; Step 101, input the CTA image to the first convolutional network, and perform blood vessel pre-segmentation based on multi-scale feature extraction;
步骤102,以预分割结果和分割标签为输入构建血管分布图,图的节点为血管概率较高、灰度一致的像素点区域,节点的形状与血管走向一致;图的边表示所连节点的相关性,边的长度小于设定的阈值; Step 102, using the pre-segmentation result and the segmentation label as input to construct a blood vessel distribution map, the nodes of the map are pixel point areas with a high probability of blood vessels and consistent gray levels, and the shape of the nodes is consistent with the direction of the blood vessels; the edges of the graph represent the connected nodes. Relevance, the length of the edge is less than the set threshold;
步骤103,将CTA图像输入到网络结构和权重参数与第一卷积网络相同的第二卷积网络,并将第二卷积网络的输出与输入CTA图像相乘后输入到第三卷积网络; Step 103, the CTA image is input to the second convolutional network with the same network structure and weight parameters as the first convolutional network, and the output of the second convolutional network is multiplied by the input CTA image and then input to the third convolutional network ;
步骤104,第三卷积网络与图卷积网络通过基于双向映射进行特征交互,基于血管分布图进行全局特征建模捕捉多尺度的局部和全局特征,实现血管区域预测; Step 104, the third convolutional network and the graph convolutional network perform feature interaction based on bidirectional mapping, and perform global feature modeling based on the blood vessel distribution map to capture multi-scale local and global features, and realize blood vessel region prediction;
步骤105,通过对多尺度特征进行融合,预测出血管分割结果。In step 105, the vessel segmentation result is predicted by fusing the multi-scale features.
本实施例中,步骤101主要用于进行血管预分割。本实施例的血管分割算法主要由卷积神经网络CNN实现,包括三个卷积网络(第一~第三卷积网络)和一个图卷积网络,如图2所示。CNN是一种前馈神经网络,但与一般的全连接前馈神经网络不同的是,它的卷积层具有局部连接和权重共享的特性,因此能够大大减小权重参数的数量,从而减小模型的复杂程度和提高运行速度。一个典型的CNN是由卷积层、汇聚层(或池化层、下采样层)、全连接层交叉堆叠而成的。卷积层的作用是通过卷积核与输入图像的卷积运算提取一个局部区域的特征,不同的卷积核相当于不同的特征提取器。汇聚层的作用是进行特征选择,降低特征数量,从而进一步减少参数数量。一般采用最大汇聚法和平均汇聚法。全连接层用于对得到的不同特征进行融合。本实施例是将CTA图像输入到第一卷积网络,采用多个不同尺寸的卷积核进行多尺度特征提取,从而血管的预分割。In this embodiment, step 101 is mainly used for blood vessel pre-segmentation. The blood vessel segmentation algorithm in this embodiment is mainly implemented by a convolutional neural network (CNN), including three convolutional networks (first to third convolutional networks) and a graph convolutional network, as shown in FIG. 2 . CNN is a feedforward neural network, but unlike the general fully connected feedforward neural network, its convolutional layer has the characteristics of local connection and weight sharing, so it can greatly reduce the number of weight parameters, thereby reducing Model complexity and increased running speed. A typical CNN is composed of convolutional layers, pooling layers (or pooling layers, downsampling layers), and fully connected layers cross-stacked. The function of the convolution layer is to extract the features of a local area through the convolution operation of the convolution kernel and the input image. Different convolution kernels are equivalent to different feature extractors. The role of the pooling layer is to perform feature selection and reduce the number of features, thereby further reducing the number of parameters. Generally, the maximum pooling method and the average pooling method are used. The fully connected layer is used to fuse the different features obtained. In this embodiment, the CTA image is input to the first convolution network, and multiple convolution kernels of different sizes are used for multi-scale feature extraction, so as to pre-segment blood vessels.
本实施例中,步骤102主要用于构建血管分布图。本实施例以上一步得到的血管预分割结果(像素点属于血管的概率)和分割标签(资深医生在影像上标注的标签即金标准)为输入构建血管分布图。分布图由节点和边组成,如图2所示。图中的节点具有以下特点:一是节点内像素点的血管概率较高;二是节点内像素 点的灰度一致;三是节点的形状符合血管走向。图中的边连接在两个节点之间,表示所连节点的相关性,边的长度小于设定的阈值。本实施例利用预分割模型构建血管分布图,能够建模或拟合血管的真实结构,有利于提高血管分割精度。In this embodiment, step 102 is mainly used to construct a blood vessel distribution map. In this embodiment, the blood vessel pre-segmentation result obtained in the previous step (the probability that a pixel belongs to a blood vessel) and the segmentation label (the label marked on the image by a senior doctor is the gold standard) are used as input to construct a blood vessel distribution map. The distribution graph consists of nodes and edges, as shown in Figure 2. The nodes in the figure have the following characteristics: first, the probability of blood vessels in the pixels in the nodes is high; second, the gray levels of the pixels in the nodes are consistent; third, the shape of the nodes conforms to the direction of blood vessels. An edge in the graph is connected between two nodes, indicating the relevance of the connected nodes, and the length of the edge is less than the set threshold. In this embodiment, a pre-segmentation model is used to construct a blood vessel distribution map, which can model or fit the real structure of blood vessels, and is beneficial to improve the accuracy of blood vessel segmentation.
本实施例中,步骤103主要用于第一阶段血管分割。本实施例基于跨网络多尺度特征融合进行血管分割。如图2所示,血管分割包括两个阶段:第一阶段采用第二卷积网络(UNET-2);第二阶段采用第三卷积网络(UNET-3)和图卷积网络(UNET-G)。第二卷积网络的网络结构与预分割的第一卷积网络相同,并使用第一卷积网络预分割模型的权重进行初始化。将CTA图像输入到第二卷积网络,得到一阶段的预测结果。一阶段的预测结果与原图进行相乘后,输入到二阶段的第三卷积网络。本实施例不将一阶段的预测结果直接作为第三卷积网络的输入,而是将其与原图相乘后作为第三卷积网络的输入,是为了防止原图信息丢失。In this embodiment, step 103 is mainly used for the first stage of blood vessel segmentation. In this embodiment, blood vessel segmentation is performed based on cross-network multi-scale feature fusion. As shown in Figure 2, blood vessel segmentation consists of two stages: the first stage uses the second convolutional network (UNET-2); the second stage uses the third convolutional network (UNET-3) and graph convolutional network (UNET- G). The network structure of the second convolutional network is the same as that of the pre-segmented first convolutional network, and is initialized with the weights of the pre-segmented model of the first convolutional network. Input the CTA image to the second convolutional network to obtain the prediction result of the first stage. After the prediction result of the first stage is multiplied by the original image, it is input to the third convolutional network of the second stage. In this embodiment, the prediction result of the first stage is not directly used as the input of the third convolutional network, but multiplied by the original image and then used as the input of the third convolutional network, in order to prevent loss of information of the original image.
本实施例中,步骤104主要用于跨网络模型的第二阶段血管分割。如图2所示,第二阶段包括第三卷积网络和图卷积网络。图卷积网络以步骤102构建的血管分布图为输入进行全局特征建模,并通过与第三卷积网络的双向(前向和反向)映射进行局部特征和全局特征的双向传播和融合,从而实现血管区域预测。In this embodiment, step 104 is mainly used for the second-stage blood vessel segmentation of the cross-network model. As shown in Figure 2, the second stage includes a third convolutional network and a graph convolutional network. The graph convolutional network takes the blood vessel distribution map constructed in step 102 as input to perform global feature modeling, and performs two-way propagation and fusion of local features and global features through bidirectional (forward and reverse) mapping with the third convolutional network, In this way, the prediction of blood vessel area is realized.
本实施例中,步骤105主要用于预测血管分割结果。第一阶段的第二卷积网络的网络结构和权重参数均与预分割的第一卷积网络相同,因此第二卷积网络同样是进行多尺度特征提取。所以,在进行局部特征和全局特征的融合后,再对多尺度特征进行融合,得到最终的血管分割。In this embodiment, step 105 is mainly used to predict the blood vessel segmentation result. The network structure and weight parameters of the second convolutional network in the first stage are the same as those of the pre-segmented first convolutional network, so the second convolutional network also performs multi-scale feature extraction. Therefore, after the fusion of local features and global features, multi-scale features are fused to obtain the final vessel segmentation.
作为一可选实施例,所述方法还包括对输入CTA图像进行归一化的步骤,将像素点的灰度归一化为[0,255],计算公式如下:As an optional embodiment, the method also includes the step of normalizing the input CTA image, normalizing the grayscale of the pixel to [0,255], the calculation formula is as follows:
Figure PCTCN2022103838-appb-000007
Figure PCTCN2022103838-appb-000007
式中,
Figure PCTCN2022103838-appb-000008
为任一像素点的灰度值x归一化后的灰度值,x min、x max分别为归一化前像素点灰度值的最小值和最大值。
In the formula,
Figure PCTCN2022103838-appb-000008
is the gray value x of any pixel after normalization, and x min and x max are the minimum and maximum values of the gray value of the pixel before normalization, respectively.
本实施例给出了一种图像预处理方法。本实施例在输入CTA图像前对其进行归一化处理,将像素点的灰度归一化为[0,255],归一化公式如上式。很显然,当 x=x min时,
Figure PCTCN2022103838-appb-000009
当x=x max时,
Figure PCTCN2022103838-appb-000010
因此,上式可以将原图像中的所有像素点的灰度值归一化到[0,255]。
This embodiment provides an image preprocessing method. In this embodiment, normalization processing is performed on the CTA image before it is input, and the grayscale of the pixel is normalized to [0, 255], and the normalization formula is as the above formula. Obviously, when x=x min ,
Figure PCTCN2022103838-appb-000009
When x=x max ,
Figure PCTCN2022103838-appb-000010
Therefore, the above formula can normalize the gray value of all pixels in the original image to [0,255].
作为一可选实施例,血管预分割方法具体包括:As an optional embodiment, the blood vessel pre-segmentation method specifically includes:
将CTA序列输入到第一卷积网络的编码层;Input the CTA sequence into the encoding layer of the first convolutional network;
编码层通过卷积和池化操作进行多尺度特征提取,第i个尺度的空间大小为原始空间大小的1/2 i,i=1,2,…,N,N为尺度个数;所述特征包含血管分割的形态和空间位置信息; The encoding layer performs multi-scale feature extraction through convolution and pooling operations, the space size of the i-th scale is 1/2 i of the original space size, i=1,2,...,N, N is the number of scales; the Features include shape and spatial location information of vessel segmentation;
将编码层提取的多尺度特征输入到解码层,解码层将每个尺度的全局特征通过上采样或反卷积映射回原始特征大小,并与对应尺度下的编码特征进行结合得到分割结果。The multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size through upsampling or deconvolution, and combines them with the encoding features at the corresponding scale to obtain the segmentation result.
本实施例给出了血管预分割的一种具体技术方案。如前述,血管预分割由第一卷积网络实现。具体地,第一卷积网络(其它两个卷积网络及图卷积网络也一样)包括编码层和解码层,CTA序列先输入到编码层,编码层采用不同的卷积核对图像进行卷积运算并经池化操作后实现多尺度特征提取。提取的不同尺度特征空间的大小构成公比为1/2的等比级数,第i个尺度的空间大小为原始空间大小的1/2 i,不同尺度的数量N一般取4~5。提取的特征包含了血管分割的形态和空间位置信息。然后将编码层提取的多尺度特征输入到解码层,解码层再将每个尺度的全局特征映射回原始特征大小,与对应尺度下的编码特征进行结合(特征级联即拼接)得到分割结果。解码层通过上采样(插值)或反卷积将全局特征映射回原始特征大小。 This embodiment provides a specific technical solution for blood vessel pre-segmentation. As mentioned above, blood vessel pre-segmentation is achieved by the first convolutional network. Specifically, the first convolutional network (the same for the other two convolutional networks and the graph convolutional network) includes an encoding layer and a decoding layer. The CTA sequence is first input to the encoding layer, and the encoding layer uses different convolution kernels to convolve the image. Multi-scale feature extraction is realized after operation and pooling operation. The size of extracted feature spaces of different scales constitutes a geometric series with a common ratio of 1/2, the space size of the i-th scale is 1/2 i of the original space size, and the number N of different scales is generally 4-5. The extracted features contain the shape and spatial location information of vessel segmentation. Then the multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size, and combines them with the encoding features at the corresponding scales (feature cascading or splicing) to obtain segmentation results. The decoding layer maps the global features back to the original feature size through upsampling (interpolation) or deconvolution.
作为一可选实施例,所述血管区域预测的方法具体包括:As an optional embodiment, the method for predicting a blood vessel region specifically includes:
第三卷积网络通过前向映射,将编码特征的多个像素点投影为图卷积网络中的单一节点;The third convolutional network projects multiple pixels of encoded features into a single node in the graph convolutional network through forward mapping;
图卷积网络将特征在血管分布图上传播,进行全局特征建模;The graph convolutional network propagates the features on the blood vessel distribution map for global feature modeling;
图卷积网络通过反向映射,将节点特征映射回第三卷积网络的特征空间上,将全局特征传播到每个局部,进行特征增强;The graph convolutional network maps the node features back to the feature space of the third convolutional network through reverse mapping, and propagates the global features to each part for feature enhancement;
将第三卷积网络和图卷积网络的特征进行结合,通过第三卷积网络捕捉影像 的局部特征,通过图卷积网络捕捉血管的全局特征,利用特征交互将局部特征和全局特征进行结合,完成血管所在区域的预测。Combine the features of the third convolutional network and the graph convolutional network, capture the local features of the image through the third convolutional network, capture the global features of blood vessels through the graph convolutional network, and use feature interaction to combine local features and global features , to complete the prediction of the area where the blood vessel is located.
本实施例给出了跨网络血管区域预测的一种技术方案。跨网络主要是指第三卷积网络和图卷积网络之间通过双向映射进行特征交互。首先是第三卷积网络进行前向映射,将编码特征的多个像素点投影为图卷积网络中的单一节点;然后,图卷积网络基于血管分布图进行全局特征建模,并进行反向映射,将节点特征映射回第三卷积网络的特征空间上,将全局特征传播到每个局部特征,进行特征增强;最后,通过第三卷积网络捕捉影像的局部特征,通过图卷积网络捕捉血管的全局特征,利用特征交互将局部特征和全局特征进行结合(特征级联),实现血管所在区域的预测。This embodiment provides a technical solution for prediction of cross-network blood vessel regions. Cross-network mainly refers to feature interaction between the third convolutional network and the graph convolutional network through bidirectional mapping. First, the third convolutional network performs forward mapping, and projects multiple pixels of encoded features into a single node in the graph convolutional network; then, the graph convolutional network performs global feature modeling based on the blood vessel distribution map, and performs inverse To map the node features back to the feature space of the third convolutional network, propagate the global features to each local feature for feature enhancement; finally, capture the local features of the image through the third convolutional network, and use the graph convolution The network captures the global features of blood vessels, and uses feature interaction to combine local features and global features (feature cascade) to realize the prediction of the area where blood vessels are located.
作为一可选实施例,所述方法采用Dice系数度量血管分割结果与医生标注的金标准的相似度,公式如下:As an optional embodiment, the method uses the Dice coefficient to measure the similarity between the blood vessel segmentation result and the gold standard marked by the doctor, and the formula is as follows:
Figure PCTCN2022103838-appb-000011
Figure PCTCN2022103838-appb-000011
式中,A、B分别为金标准的血管像素点集合和血管分割结果的像素点集合,|A|、|B|和|A∩B|分别为A、B和A与B的交集中的像素点的数量。In the formula, A and B are the blood vessel pixel point set of the gold standard and the pixel point set of the blood vessel segmentation result respectively, and |A|, |B| and |A∩B| are A, B and the intersection of A and B respectively The number of pixels.
本实施例给出了定量评价血管分割结果的一种技术方案。本实施例采用Dice系数表示血管分割结果与医生标注的金标准的相似度,公式如上式。公式中的分子是分割结果和金标准重合的像素点的数量的2倍,分母是分割结果和金标准的像素点的和。根据公式,分割结果和金标准重合的像素点越多,相似度越大;当二者完全重合时,相似度最大为1。This embodiment provides a technical solution for quantitatively evaluating blood vessel segmentation results. In this embodiment, the Dice coefficient is used to represent the similarity between the blood vessel segmentation result and the gold standard marked by doctors, and the formula is as above. The numerator in the formula is twice the number of pixels where the segmentation result coincides with the gold standard, and the denominator is the sum of the pixel points of the segmentation result and the gold standard. According to the formula, the more pixels that overlap the segmentation result and the gold standard, the greater the similarity; when the two completely overlap, the maximum similarity is 1.
图3为本发明实施例一种血管分割装置的组成示意图,所述装置包括:Fig. 3 is a schematic diagram of the composition of a blood vessel segmentation device according to an embodiment of the present invention, and the device includes:
预分割模块11,用于将CTA图像输入到第一卷积网络,进行基于多尺度特征提取的血管预分割;The pre-segmentation module 11 is used to input the CTA image to the first convolutional network for pre-segmentation of blood vessels based on multi-scale feature extraction;
分布图构建模块12,用于以预分割结果和分割标签为输入构建血管分布图,图的节点为血管概率较高、灰度一致的像素点区域,节点的形状与血管走向一致; 图的边表示所连节点的相关性,边的长度小于设定的阈值;The distribution map construction module 12 is used to construct a blood vessel distribution map with pre-segmentation results and segmentation labels as input. The nodes of the map are pixel areas with high blood vessel probability and consistent gray scale, and the shape of the nodes is consistent with the direction of the blood vessels; Indicates the relevance of the connected nodes, and the length of the edge is less than the set threshold;
第一阶段分割模块13,用于将CTA图像输入到网络结构和权重参数与第一卷积网络相同的第二卷积网络,并将第二卷积网络的输出与输入CTA图像相乘后输入到第三卷积网络;The first stage segmentation module 13 is used to input the CTA image to the second convolutional network with the same network structure and weight parameters as the first convolutional network, and input the output after multiplying the output of the second convolutional network with the input CTA image to the third convolutional network;
第二阶段分割模块14,用于第三卷积网络与图卷积网络通过基于双向映射进行特征交互,基于血管分布图进行全局特征建模捕捉多尺度的局部和全局特征,实现血管区域预测;The second-stage segmentation module 14 is used for feature interaction between the third convolutional network and the graph convolutional network based on bidirectional mapping, and global feature modeling based on the blood vessel distribution map to capture multi-scale local and global features and realize blood vessel region prediction;
特征融合模块15,用于通过对多尺度特征进行融合,预测出血管分割结果。The feature fusion module 15 is configured to predict the blood vessel segmentation result by fusing multi-scale features.
本实施例的装置,可以用于执行图1所示方法实施例的技术方案,其实现原理和技术效果类似,此处不再赘述。后面的实施例也是如此,均不再展开说明。The device of this embodiment can be used to implement the technical solution of the method embodiment shown in FIG. 1 , and its implementation principle and technical effect are similar, and will not be repeated here. The same is true for the following embodiments, which will not be further described.
作为一可选实施例,所述装置还包括对输入CTA图像进行预处理的归一化模块,用于将像素点的灰度归一化为[0,255],计算公式如下:As an optional embodiment, the device also includes a normalization module for preprocessing the input CTA image, which is used to normalize the grayscale of the pixels to [0,255], and the calculation formula is as follows:
Figure PCTCN2022103838-appb-000012
Figure PCTCN2022103838-appb-000012
式中,
Figure PCTCN2022103838-appb-000013
为任一像素点的灰度值x归一化后的灰度值,x min、x max分别为归一化前像素点灰度值的最小值和最大值。
In the formula,
Figure PCTCN2022103838-appb-000013
is the gray value x of any pixel after normalization, and x min and x max are the minimum and maximum values of the gray value of the pixel before normalization, respectively.
作为一可选实施例,血管预分割方法具体包括:As an optional embodiment, the blood vessel pre-segmentation method specifically includes:
将CTA序列输入到第一卷积网络的编码层;Input the CTA sequence into the encoding layer of the first convolutional network;
编码层通过卷积和池化操作进行多尺度特征提取,第i个尺度的空间大小为原始空间大小的1/2 i,i=1,2,…,N,N为尺度个数;所述特征包含血管分割的形态和空间位置信息; The encoding layer performs multi-scale feature extraction through convolution and pooling operations, the space size of the i-th scale is 1/2 i of the original space size, i=1,2,...,N, N is the number of scales; the Features include shape and spatial location information of vessel segmentation;
将编码层提取的多尺度特征输入到解码层,解码层将每个尺度的全局特征通过上采样或反卷积映射回原始特征大小,并与对应尺度下的编码特征进行结合得到分割结果。The multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size through upsampling or deconvolution, and combines them with the encoding features at the corresponding scale to obtain the segmentation result.
作为一可选实施例,所述血管区域预测的方法具体包括:As an optional embodiment, the method for predicting a blood vessel region specifically includes:
第三卷积网络通过前向映射,将编码特征的多个像素点投影为图卷积网络中的单一节点;The third convolutional network projects multiple pixels of encoded features into a single node in the graph convolutional network through forward mapping;
图卷积网络将特征在血管分布图上传播,进行全局特征建模;The graph convolutional network propagates the features on the blood vessel distribution map for global feature modeling;
图卷积网络通过反向映射,将节点特征映射回第三卷积网络的特征空间上,将全局特征传播到每个局部,进行特征增强;The graph convolutional network maps the node features back to the feature space of the third convolutional network through reverse mapping, and propagates the global features to each part for feature enhancement;
将第三卷积网络和图卷积网络的特征进行结合,通过第三卷积网络捕捉影像的局部特征,通过图卷积网络捕捉血管的全局特征,利用特征交互将局部特征和全局特征进行结合,完成血管所在区域的预测。Combine the features of the third convolutional network and the graph convolutional network, capture the local features of the image through the third convolutional network, capture the global features of blood vessels through the graph convolutional network, and use feature interaction to combine local features and global features , to complete the prediction of the area where the blood vessel is located.
作为一可选实施例,所述方法采用Dice系数度量血管分割结果与医生标注的金标准的相似度,公式如下:As an optional embodiment, the method uses the Dice coefficient to measure the similarity between the blood vessel segmentation result and the gold standard marked by the doctor, and the formula is as follows:
Figure PCTCN2022103838-appb-000014
Figure PCTCN2022103838-appb-000014
式中,A、B分别为金标准的血管像素点集合和血管分割结果的像素点集合,|A|、|B|和|A∩B|分别为A、B和A与B的交集中的像素点的数量。In the formula, A and B are the blood vessel pixel point set of the gold standard and the pixel point set of the blood vessel segmentation result respectively, and |A|, |B| and |A∩B| are A, B and the intersection of A and B respectively The number of pixels.
图4是本申请实施例提供的一种电子设备的结构示意图。在硬件层面,该电子设备包括处理器,可选地还包括内部总线、网络接口、存储器。其中,存储器可能包含内存,例如高速随机存取存储器(Random-Access Memory,RAM),也可能还包括非易失性存储器(non-volatile memory),例如至少1个磁盘存储器等。当然,该电子设备还可能包括其他业务所需要的硬件。FIG. 4 is a schematic structural diagram of an electronic device provided by an embodiment of the present application. At the hardware level, the electronic device includes a processor, and optionally also includes an internal bus, a network interface, and a memory. Wherein, the memory may include a memory, such as a high-speed random-access memory (Random-Access Memory, RAM), and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory. Of course, the electronic device may also include hardware required by other services.
处理器、网络接口和存储器可以通过内部总线相互连接,该内部总线可以是ISA(Industry Standard Architecture,工业标准体系结构)总线、PCI(Peripheral Component Interconnect,外设部件互连标准)总线或EISA(Extended Industry Standard Architecture,扩展工业标准结构)总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示,图4中仅用一个双向箭头表示,但并不表示仅有一根总线或一种类型的总线。The processor, the network interface and the memory can be connected to each other through an internal bus, which can be an ISA (Industry Standard Architecture, industry standard architecture) bus, a PCI (Peripheral Component Interconnect, peripheral component interconnection standard) bus or an EISA (Extended Industry Standard Architecture, extended industry standard architecture) bus, etc. The bus can be divided into address bus, data bus, control bus and so on. For ease of representation, only one double-headed arrow is used in FIG. 4 , but it does not mean that there is only one bus or one type of bus.
存储器,用于存放执行指令。具体地,执行指令即可被执行的计算机程序。存储器可以包括内存和非易失性存储器,并向处理器提供执行指令和数据。Memory for storing execution instructions. Specifically, a computer program that can be executed by executing instructions. The memory, which can include internal memory and non-volatile memory, provides instructions and data to the processor for execution.
在一种可能实现的方式中,处理器从非易失性存储器中读取对应的执行指令 到内存中然后运行,也可从其它设备上获取相应的执行指令,以在逻辑层面上形成血管分割装置。处理器执行存储器所存放的执行指令,以通过执行的执行指令实现本申请任一实施例中提供的血管分割方法。In a possible implementation, the processor reads the corresponding execution instructions from the non-volatile memory into the memory and then runs them. It can also obtain the corresponding execution instructions from other devices to form blood vessel segmentation at the logical level. device. The processor executes the execution instructions stored in the memory, so as to implement the blood vessel segmentation method provided in any embodiment of the present application through the executed execution instructions.
上述如本申请图1所示实施例提供的血管分割装置执行的方法可以应用于处理器中,或者由处理器实现。处理器可能是一种集成电路芯片,具有信号的处理能力。在实现过程中,上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器可以是通用处理器,包括中央处理器(Central Processing Unit,CPU)、网络处理器(Network Processor,NP)等;还可以是数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。The above-mentioned method performed by the blood vessel segmentation apparatus provided in the embodiment shown in FIG. 1 of the present application may be applied to a processor or implemented by the processor. A processor may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by an integrated logic circuit of hardware in a processor or an instruction in the form of software. The above-mentioned processor can be a general-purpose processor, including a central processing unit (Central Processing Unit, CPU), a network processor (Network Processor, NP), etc.; it can also be a digital signal processor (Digital Signal Processor, DSP), a dedicated integrated Circuit (Application Specific Integrated Circuit, ASIC), Field-Programmable Gate Array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components. Various methods, steps, and logic block diagrams disclosed in the embodiments of the present application may be implemented or executed. A general-purpose processor may be a microprocessor, or the processor may be any conventional processor, or the like.
结合本申请实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成,或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器,处理器读取存储器中的信息,结合其硬件完成上述方法的步骤。The steps of the method disclosed in connection with the embodiments of the present application may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module can be located in a mature storage medium in the field such as random access memory, flash memory, read-only memory, programmable read-only memory or electrically erasable programmable memory, register. The storage medium is located in the memory, and the processor reads the information in the memory, and completes the steps of the above method in combination with its hardware.
本申请实施例还提出了一种可读介质,该可读存储介质存储有执行指令,存储的执行指令被电子设备的处理器执行时,能够使该电子设备执行本申请任一实施例中提供的血管分割方法,并具体用于执行上述血管分割方法。The embodiment of the present application also proposes a readable medium, the readable storage medium stores execution instructions, and when the stored execution instructions are executed by the processor of the electronic device, the electronic device can execute the electronic device provided in any embodiment of the present application. The blood vessel segmentation method, and is specifically used to implement the above blood vessel segmentation method.
前述各个实施例中所述的电子设备可以为计算机。The electronic equipment described in the foregoing embodiments may be a computer.
本领域内的技术人员应明白,本申请的实施例可提供为方法或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例,或软件和硬件相结合的形式。Those skilled in the art should understand that the embodiments of the present application may be provided as methods or computer program products. Therefore, the present application may adopt an entirely hardware embodiment, an entirely software embodiment, or a combination of software and hardware.
本申请中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。Each embodiment in the present application is described in a progressive manner, the same and similar parts of each embodiment can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, as for the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and for relevant parts, please refer to part of the description of the method embodiment.
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。It should also be noted that the term "comprises", "comprises" or any other variation thereof is intended to cover a non-exclusive inclusion such that a process, method, article, or apparatus comprising a set of elements includes not only those elements, but also includes Other elements not expressly listed, or elements inherent in the process, method, commodity, or apparatus are also included. Without further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article or apparatus comprising said element.
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。The above descriptions are only examples of the present application, and are not intended to limit the present application. For those skilled in the art, various modifications and changes may occur in this application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included within the scope of the claims of the present application.
工业实用性Industrial Applicability
本发明提供的一种血管分割方法、装置、电子设备和可读介质利用图像识别技术结合构建血管分布网络对血管全局特征进行建模,利用血管的整体结构对分割区域的特征进行增强,挖掘血管结构的潜在走向,充分利用了计算机技术中程序自动化的处理基础,使得在血管造影影像的复杂成像质量环境下对血管分割区域预测和分割的准确性获得极大地提高。形成的产品可以批量生产,快速应用于对血管成像具有高需求的系统或场景。A blood vessel segmentation method, device, electronic equipment, and readable medium provided by the present invention use image recognition technology to build a blood vessel distribution network to model the global characteristics of blood vessels, use the overall structure of blood vessels to enhance the features of the segmented area, and mine blood vessels The potential direction of the structure makes full use of the processing basis of program automation in computer technology, which greatly improves the accuracy of the prediction and segmentation of vessel segmentation regions under the complex imaging quality environment of angiographic images. The formed products can be mass-produced and quickly applied to systems or scenarios with high demand for vascular imaging.

Claims (12)

  1. 一种血管分割方法,其特征在于,包括以下步骤:A blood vessel segmentation method, characterized in that, comprising the following steps:
    将CTA图像输入到第一卷积网络,进行基于多尺度特征提取的血管预分割;Input the CTA image to the first convolutional network for pre-segmentation of blood vessels based on multi-scale feature extraction;
    以预分割结果和分割标签为输入构建血管分布图,图的节点为血管概率较高、灰度一致的像素点区域,节点的形状与血管走向一致;图的边表示所连节点的相关性,边的长度小于设定的阈值;The pre-segmentation results and segmentation labels are used as input to construct a blood vessel distribution map. The nodes of the graph are pixel areas with high probability of blood vessels and consistent gray levels. The shape of the nodes is consistent with the direction of the blood vessels; the edges of the graph represent the correlation of connected nodes. The length of the edge is less than the set threshold;
    将CTA图像输入到网络结构和权重参数与第一卷积网络相同的第二卷积网络,并将第二卷积网络的输出与输入CTA图像相乘后输入到第三卷积网络;The CTA image is input to the second convolutional network with the same network structure and weight parameters as the first convolutional network, and the output of the second convolutional network is multiplied by the input CTA image and then input to the third convolutional network;
    第三卷积网络与图卷积网络通过基于双向映射进行特征交互,基于血管分布图进行全局特征建模捕捉多尺度的局部和全局特征,实现血管区域预测;The third convolutional network and the graph convolutional network perform feature interaction based on bidirectional mapping, and perform global feature modeling based on the vascular distribution map to capture multi-scale local and global features, and realize vascular region prediction;
    通过对多尺度特征进行融合,预测出血管分割结果。Through the fusion of multi-scale features, the result of blood vessel segmentation is predicted.
  2. 根据权利要求1所述的血管分割方法,其特征在于,所述方法还包括对输入CTA图像进行归一化的步骤,将像素点的灰度归一化为[0,255],计算公式如下:The blood vessel segmentation method according to claim 1, characterized in that the method also includes the step of normalizing the input CTA image, normalizing the grayscale of the pixels to [0,255], and the calculation formula is as follows:
    Figure PCTCN2022103838-appb-100001
    Figure PCTCN2022103838-appb-100001
    式中,
    Figure PCTCN2022103838-appb-100002
    为任一像素点的灰度值x归一化后的灰度值,x min、x max分别为归一化前像素点灰度值的最小值和最大值。
    In the formula,
    Figure PCTCN2022103838-appb-100002
    is the gray value x of any pixel after normalization, and x min and x max are the minimum and maximum values of the gray value of the pixel before normalization, respectively.
  3. 根据权利要求1所述的血管分割方法,其特征在于,血管预分割方法具体包括:The blood vessel segmentation method according to claim 1, wherein the blood vessel pre-segmentation method specifically comprises:
    将CTA序列输入到第一卷积网络的编码层;Input the CTA sequence into the encoding layer of the first convolutional network;
    编码层通过卷积和池化操作进行多尺度特征提取,第i个尺度的空间大小为原始空间大小的1/2 i,i=1,2,…,N,N为尺度个数;所述特征包含血管分割的形态和空间位置信息; The encoding layer performs multi-scale feature extraction through convolution and pooling operations, the space size of the i-th scale is 1/2 i of the original space size, i=1,2,...,N, N is the number of scales; the Features include shape and spatial location information of vessel segmentation;
    将编码层提取的多尺度特征输入到解码层,解码层将每个尺度的全局特征通过上采样或反卷积映射回原始特征大小,并与对应尺度下的编码特征进行结合得到分割结果。The multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size through upsampling or deconvolution, and combines them with the encoding features at the corresponding scale to obtain the segmentation result.
  4. 根据权利要求1所述的血管分割方法,其特征在于,所述血管区域预测的 方法具体包括:The blood vessel segmentation method according to claim 1, wherein the method for predicting the blood vessel region specifically comprises:
    第三卷积网络通过前向映射,将编码特征的多个像素点投影为图卷积网络中的单一节点;The third convolutional network projects multiple pixels of encoded features into a single node in the graph convolutional network through forward mapping;
    图卷积网络将特征在血管分布图上传播,进行全局特征建模;The graph convolutional network propagates the features on the blood vessel distribution map for global feature modeling;
    图卷积网络通过反向映射,将节点特征映射回第三卷积网络的特征空间上,将全局特征传播到每个局部,进行特征增强;The graph convolutional network maps the node features back to the feature space of the third convolutional network through reverse mapping, and propagates the global features to each part for feature enhancement;
    将第三卷积网络和图卷积网络的特征进行结合,通过第三卷积网络捕捉影像的局部特征,通过图卷积网络捕捉血管的全局特征,利用特征交互将局部特征和全局特征进行结合,完成血管所在区域的预测。Combine the features of the third convolutional network and the graph convolutional network, capture the local features of the image through the third convolutional network, capture the global features of blood vessels through the graph convolutional network, and use feature interaction to combine local features and global features , to complete the prediction of the area where the blood vessel is located.
  5. 根据权利要求1所述的血管分割方法,其特征在于,所述方法采用Dice系数度量血管分割结果与医生标注的金标准的相似度,公式如下:The blood vessel segmentation method according to claim 1, wherein the method uses the Dice coefficient to measure the similarity between the blood vessel segmentation result and the gold standard marked by a doctor, and the formula is as follows:
    Figure PCTCN2022103838-appb-100003
    Figure PCTCN2022103838-appb-100003
    式中,A、B分别为金标准的血管像素点集合和血管分割结果的像素点集合,|A|、|B|和|A∩B|分别为A、B和A与B的交集中的像素点的数量。In the formula, A and B are the blood vessel pixel point set of the gold standard and the pixel point set of the blood vessel segmentation result respectively, and |A|, |B| and |A∩B| are A, B and the intersection of A and B respectively The number of pixels.
  6. 一种血管分割装置,其特征在于,包括:A blood vessel segmentation device, characterized in that it comprises:
    预分割模块,用于将CTA图像输入到第一卷积网络,进行基于多尺度特征提取的血管预分割;The pre-segmentation module is used to input the CTA image to the first convolutional network for pre-segmentation of blood vessels based on multi-scale feature extraction;
    分布图构建模块,用于以预分割结果和分割标签为输入构建血管分布图,图的节点为血管概率较高、灰度一致的像素点区域,节点的形状与血管走向一致;图的边表示所连节点的相关性,边的长度小于设定的阈值;The distribution map building module is used to construct a blood vessel distribution map with pre-segmentation results and segmentation labels as input. The nodes of the map are pixel areas with a high probability of blood vessels and consistent gray levels. The shape of the nodes is consistent with the direction of the blood vessels; the edges of the graph represent The relevance of the connected nodes, the length of the edge is less than the set threshold;
    第一阶段分割模块,用于将CTA图像输入到网络结构和权重参数与第一卷积网络相同的第二卷积网络,并将第二卷积网络的输出与输入CTA图像相乘后输入到第三卷积网络;The first-stage segmentation module is used to input the CTA image to the second convolutional network with the same network structure and weight parameters as the first convolutional network, and to multiply the output of the second convolutional network with the input CTA image and then input it to The third convolutional network;
    第二阶段分割模块,用于第三卷积网络与图卷积网络通过基于双向映射进行特征交互,基于血管分布图进行全局特征建模捕捉多尺度的局部和全局特征,实现血管区域预测;The second stage segmentation module is used for the third convolutional network and the graph convolutional network to perform feature interaction based on two-way mapping, and to perform global feature modeling based on the blood vessel distribution map to capture multi-scale local and global features, and realize blood vessel region prediction;
    特征融合模块,用于通过对多尺度特征进行融合,预测出血管分割结果。The feature fusion module is used to predict the blood vessel segmentation result by fusing multi-scale features.
  7. 根据权利要求6所述的血管分割装置,其特征在于,所述装置还包括对输入CTA图像进行预处理的归一化模块,用于将像素点的灰度归一化为[0,255],计算公式如下:The blood vessel segmentation device according to claim 6, wherein the device further comprises a normalization module for preprocessing the input CTA image, which is used to normalize the grayscale of the pixel to [0, 255], and calculate The formula is as follows:
    Figure PCTCN2022103838-appb-100004
    Figure PCTCN2022103838-appb-100004
    式中,
    Figure PCTCN2022103838-appb-100005
    为任一像素点的灰度值x归一化后的灰度值,x min、x max分别为归一化前像素点灰度值的最小值和最大值。
    In the formula,
    Figure PCTCN2022103838-appb-100005
    is the gray value x of any pixel after normalization, and x min and x max are the minimum and maximum values of the gray value of the pixel before normalization, respectively.
  8. 根据权利要求6所述的血管分割装置,其特征在于,血管预分割方法具体包括:The blood vessel segmentation device according to claim 6, wherein the blood vessel pre-segmentation method specifically comprises:
    将CTA序列输入到第一卷积网络的编码层;Input the CTA sequence into the encoding layer of the first convolutional network;
    编码层通过卷积和池化操作进行多尺度特征提取,第i个尺度的空间大小为原始空间大小的1/2 i,i=1,2,…,N,N为尺度个数;所述特征包含血管分割的形态和空间位置信息; The encoding layer performs multi-scale feature extraction through convolution and pooling operations, the space size of the i-th scale is 1/2 i of the original space size, i=1,2,...,N, N is the number of scales; the Features include shape and spatial location information of vessel segmentation;
    将编码层提取的多尺度特征输入到解码层,解码层将每个尺度的全局特征通过上采样或反卷积映射回原始特征大小,并与对应尺度下的编码特征进行结合得到分割结果。The multi-scale features extracted by the encoding layer are input to the decoding layer, and the decoding layer maps the global features of each scale back to the original feature size through upsampling or deconvolution, and combines them with the encoding features at the corresponding scale to obtain the segmentation result.
  9. 根据权利要求6所述的血管分割装置,其特征在于,所述血管区域预测的方法具体包括:The blood vessel segmentation device according to claim 6, wherein the method for predicting the blood vessel region specifically comprises:
    第三卷积网络通过前向映射,将编码特征的多个像素点投影为图卷积网络中的单一节点;The third convolutional network projects multiple pixels of encoded features into a single node in the graph convolutional network through forward mapping;
    图卷积网络将特征在血管分布图上传播,进行全局特征建模;The graph convolutional network propagates the features on the blood vessel distribution map for global feature modeling;
    图卷积网络通过反向映射,将节点特征映射回第三卷积网络的特征空间上,将全局特征传播到每个局部,进行特征增强;The graph convolutional network maps the node features back to the feature space of the third convolutional network through reverse mapping, and propagates the global features to each part for feature enhancement;
    将第三卷积网络和图卷积网络的特征进行结合,通过第三卷积网络捕捉影像的局部特征,通过图卷积网络捕捉血管的全局特征,利用特征交互将局部特征和全局特征进行结合,完成血管所在区域的预测。Combine the features of the third convolutional network and the graph convolutional network, capture the local features of the image through the third convolutional network, capture the global features of blood vessels through the graph convolutional network, and use feature interaction to combine local features and global features , to complete the prediction of the area where the blood vessel is located.
  10. 根据权利要求6所述的血管分割装置,其特征在于,所述方法采用Dice系数度量血管分割结果与医生标注的金标准的相似度,公式如下:The blood vessel segmentation device according to claim 6, wherein the method uses the Dice coefficient to measure the similarity between the blood vessel segmentation result and the gold standard marked by the doctor, and the formula is as follows:
    Figure PCTCN2022103838-appb-100006
    Figure PCTCN2022103838-appb-100006
    式中,A、B分别为金标准的血管像素点集合和血管分割结果的像素点集合,|A|、|B|和|A∩B|分别为A、B和A与B的交集中的像素点的数量。In the formula, A and B are the blood vessel pixel point set of the gold standard and the pixel point set of the blood vessel segmentation result respectively, and |A|, |B| and |A∩B| are A, B and the intersection of A and B respectively The number of pixels.
  11. 一种可读介质,其特征在于,所述可读介质包括执行指令,当电子设备的处理器执行所述执行指令时,所述电子设备执行如权利要求1-5中任一所述的血管分割方法。A readable medium, characterized in that the readable medium includes execution instructions, and when the processor of the electronic device executes the execution instructions, the electronic device executes the blood vessel according to any one of claims 1-5. split method.
  12. 一种电子设备,其特征在于,所述电子设备包括处理器以及存储有执行指令的存储器,当所述处理器执行所述存储器存储的所述执行指令时,所述处理器执行如权利要求1-5中任一所述的血管分割方法。An electronic device, characterized in that the electronic device includes a processor and a memory storing execution instructions, and when the processor executes the execution instructions stored in the memory, the processor executes the method according to claim 1. - the blood vessel segmentation method described in any one of 5.
PCT/CN2022/103838 2021-07-07 2022-07-05 Blood vessel segmentation method and apparatus, and electronic device and readable medium WO2023280148A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202110768617.8 2021-07-07
CN202110768617.8A CN113554665A (en) 2021-07-07 2021-07-07 Blood vessel segmentation method and device

Publications (1)

Publication Number Publication Date
WO2023280148A1 true WO2023280148A1 (en) 2023-01-12

Family

ID=78131406

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/103838 WO2023280148A1 (en) 2021-07-07 2022-07-05 Blood vessel segmentation method and apparatus, and electronic device and readable medium

Country Status (2)

Country Link
CN (1) CN113554665A (en)
WO (1) WO2023280148A1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116523888A (en) * 2023-05-08 2023-08-01 北京天鼎殊同科技有限公司 Pavement crack detection method, device, equipment and medium
CN116342588B (en) * 2023-05-22 2023-08-11 徕兄健康科技(威海)有限责任公司 Cerebrovascular image enhancement method
CN116630386A (en) * 2023-06-12 2023-08-22 新疆生产建设兵团医院 CTA scanning image processing method and system thereof
CN117726633A (en) * 2024-02-07 2024-03-19 安徽大学 Segmentation method and system of double-branch coronary artery image based on feature fusion
CN118015017A (en) * 2024-02-06 2024-05-10 中国科学院宁波材料技术与工程研究所 Training method and device for segmentation model, electronic equipment and storage medium

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113554665A (en) * 2021-07-07 2021-10-26 杭州深睿博联科技有限公司 Blood vessel segmentation method and device
CN114723683B (en) * 2022-03-22 2023-02-17 推想医疗科技股份有限公司 Head and neck artery blood vessel segmentation method and device, electronic device and storage medium
CN115239744A (en) * 2022-07-25 2022-10-25 推想医疗科技股份有限公司 Blood vessel segmentation method, device, equipment and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112017190A (en) * 2020-08-06 2020-12-01 杭州深睿博联科技有限公司 Global network construction and training method and device for vessel segmentation completion
CN112329801A (en) * 2020-12-03 2021-02-05 中国石油大学(华东) Convolutional neural network non-local information construction method
US20210064959A1 (en) * 2019-08-27 2021-03-04 Nec Laboratories America, Inc. Flexible edge-empowered graph convolutional networks with node-edge enhancement
CN112766280A (en) * 2021-01-16 2021-05-07 北京工业大学 Remote sensing image road extraction method based on graph convolution
CN113554665A (en) * 2021-07-07 2021-10-26 杭州深睿博联科技有限公司 Blood vessel segmentation method and device

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111815598B (en) * 2020-06-30 2024-04-26 上海联影医疗科技股份有限公司 Calculation method, device and equipment of vascular parameters and storage medium
CN112053363B (en) * 2020-08-19 2023-12-15 苏州超云生命智能产业研究院有限公司 Retina blood vessel segmentation method, retina blood vessel segmentation device and model construction method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20210064959A1 (en) * 2019-08-27 2021-03-04 Nec Laboratories America, Inc. Flexible edge-empowered graph convolutional networks with node-edge enhancement
CN112017190A (en) * 2020-08-06 2020-12-01 杭州深睿博联科技有限公司 Global network construction and training method and device for vessel segmentation completion
CN112329801A (en) * 2020-12-03 2021-02-05 中国石油大学(华东) Convolutional neural network non-local information construction method
CN112766280A (en) * 2021-01-16 2021-05-07 北京工业大学 Remote sensing image road extraction method based on graph convolution
CN113554665A (en) * 2021-07-07 2021-10-26 杭州深睿博联科技有限公司 Blood vessel segmentation method and device

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116523888A (en) * 2023-05-08 2023-08-01 北京天鼎殊同科技有限公司 Pavement crack detection method, device, equipment and medium
CN116523888B (en) * 2023-05-08 2023-11-03 北京天鼎殊同科技有限公司 Pavement crack detection method, device, equipment and medium
CN116342588B (en) * 2023-05-22 2023-08-11 徕兄健康科技(威海)有限责任公司 Cerebrovascular image enhancement method
CN116630386A (en) * 2023-06-12 2023-08-22 新疆生产建设兵团医院 CTA scanning image processing method and system thereof
CN116630386B (en) * 2023-06-12 2024-02-20 新疆生产建设兵团医院 CTA scanning image processing method and system thereof
CN118015017A (en) * 2024-02-06 2024-05-10 中国科学院宁波材料技术与工程研究所 Training method and device for segmentation model, electronic equipment and storage medium
CN117726633A (en) * 2024-02-07 2024-03-19 安徽大学 Segmentation method and system of double-branch coronary artery image based on feature fusion
CN117726633B (en) * 2024-02-07 2024-04-19 安徽大学 Segmentation method and system of double-branch coronary artery image based on feature fusion

Also Published As

Publication number Publication date
CN113554665A (en) 2021-10-26

Similar Documents

Publication Publication Date Title
WO2023280148A1 (en) Blood vessel segmentation method and apparatus, and electronic device and readable medium
WO2019120110A1 (en) Image reconstruction method and device
EP4099220A1 (en) Processing apparatus, method and storage medium
CN114897779A (en) Cervical cytology image abnormal area positioning method and device based on fusion attention
WO2020168820A1 (en) Yolo convolutional neural network-based cholelithiasis ct medical image data enhancement method
Huang et al. Medical image segmentation using deep learning with feature enhancement
CN113378933A (en) Thyroid ultrasound image classification and segmentation network, training method, device and medium
CN113592060A (en) Neural network optimization method and device
CN115631112B (en) Building contour correction method and device based on deep learning
WO2023207778A1 (en) Data recovery method and device, computer, and storage medium
WO2021227058A1 (en) Text processing method and apparatus, and electronic device and storage medium
CN114612681A (en) GCN-based multi-label image classification method, model construction method and device
CN108932715A (en) A kind of optimization method of the coronarogram segmentation based on deep learning
WO2024179485A1 (en) Image processing method and related device thereof
Wu et al. Pneumonia detection based on RSNA dataset and anchor-free deep learning detector
WO2024011859A1 (en) Neural network-based face detection method and device
CN113011415A (en) Improved target detection method and system based on Grid R-CNN model
CN115496744B (en) Lung cancer image segmentation method, device, terminal and medium based on mixed attention
WO2023197857A1 (en) Model partitioning method and related device thereof
CN112884702A (en) Polyp identification system and method based on endoscope image
CN116703945A (en) Tumor segmentation method, system, medium and terminal for breast ultrasound image
CN115690143B (en) Image segmentation method, device, electronic equipment and storage medium
CN113610856B (en) Method and device for training image segmentation model and image segmentation
GB2611601A (en) Image processing method and device
CN115079818A (en) Hand capturing method and system

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: 22836899

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: 22836899

Country of ref document: EP

Kind code of ref document: A1