CN112669257A - Intracranial vascular lesion identification method based on transfer learning - Google Patents
Intracranial vascular lesion identification method based on transfer learning Download PDFInfo
- Publication number
- CN112669257A CN112669257A CN202011324149.7A CN202011324149A CN112669257A CN 112669257 A CN112669257 A CN 112669257A CN 202011324149 A CN202011324149 A CN 202011324149A CN 112669257 A CN112669257 A CN 112669257A
- Authority
- CN
- China
- Prior art keywords
- image
- blood
- bright
- black
- enhanced
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 174
- 238000007917 intracranial administration Methods 0.000 title claims abstract description 130
- 231100000216 vascular lesion Toxicity 0.000 title claims abstract description 24
- 238000013526 transfer learning Methods 0.000 title claims abstract description 12
- 210000004369 blood Anatomy 0.000 claims abstract description 529
- 239000008280 blood Substances 0.000 claims abstract description 529
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 187
- 238000005070 sampling Methods 0.000 claims abstract description 22
- 238000013508 migration Methods 0.000 claims abstract description 12
- 230000005012 migration Effects 0.000 claims abstract description 12
- 238000007781 pre-processing Methods 0.000 claims abstract description 10
- 230000002792 vascular Effects 0.000 claims description 49
- 238000002583 angiography Methods 0.000 claims description 44
- 238000010586 diagram Methods 0.000 claims description 33
- 230000011218 segmentation Effects 0.000 claims description 29
- 230000009466 transformation Effects 0.000 claims description 27
- 239000011159 matrix material Substances 0.000 claims description 25
- 238000012545 processing Methods 0.000 claims description 23
- 208000031481 Pathologic Constriction Diseases 0.000 claims description 21
- 230000036262 stenosis Effects 0.000 claims description 21
- 208000037804 stenosis Diseases 0.000 claims description 21
- 238000012360 testing method Methods 0.000 claims description 20
- 238000013528 artificial neural network Methods 0.000 claims description 16
- 230000007797 corrosion Effects 0.000 claims description 15
- 238000005260 corrosion Methods 0.000 claims description 15
- 238000005259 measurement Methods 0.000 claims description 13
- 238000007667 floating Methods 0.000 claims description 10
- 230000004927 fusion Effects 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 7
- 238000005282 brightening Methods 0.000 claims description 6
- 238000011524 similarity measure Methods 0.000 claims description 4
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 230000002708 enhancing effect Effects 0.000 claims description 2
- 230000008030 elimination Effects 0.000 abstract description 9
- 238000003379 elimination reaction Methods 0.000 abstract description 9
- 230000005291 magnetic effect Effects 0.000 description 40
- 230000000694 effects Effects 0.000 description 26
- 230000008569 process Effects 0.000 description 21
- 238000003384 imaging method Methods 0.000 description 20
- 238000004458 analytical method Methods 0.000 description 16
- 238000004422 calculation algorithm Methods 0.000 description 16
- 238000005516 engineering process Methods 0.000 description 12
- 238000011156 evaluation Methods 0.000 description 8
- 206010057469 Vascular stenosis Diseases 0.000 description 7
- 239000003086 colorant Substances 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 238000000605 extraction Methods 0.000 description 6
- 238000004445 quantitative analysis Methods 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 238000004451 qualitative analysis Methods 0.000 description 5
- 210000001519 tissue Anatomy 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 230000003628 erosive effect Effects 0.000 description 3
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 3
- 238000012805 post-processing Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 208000006011 Stroke Diseases 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010968 computed tomography angiography Methods 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 230000034994 death Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 238000003709 image segmentation Methods 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 201000001320 Atherosclerosis Diseases 0.000 description 1
- 206010008111 Cerebral haemorrhage Diseases 0.000 description 1
- 206010008190 Cerebrovascular accident Diseases 0.000 description 1
- 206010069729 Collateral circulation Diseases 0.000 description 1
- 206010028851 Necrosis Diseases 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 210000005013 brain tissue Anatomy 0.000 description 1
- 230000002308 calcification Effects 0.000 description 1
- 210000001715 carotid artery Anatomy 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 206010008118 cerebral infarction Diseases 0.000 description 1
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 230000004087 circulation Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000005802 health problem Effects 0.000 description 1
- 238000003711 image thresholding Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 201000005851 intracranial arteriosclerosis Diseases 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000017074 necrotic cell death Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000005298 paramagnetic effect Effects 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 210000005259 peripheral blood Anatomy 0.000 description 1
- 239000011886 peripheral blood Substances 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 210000002385 vertebral artery Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
- G06T5/30—Erosion or dilatation, e.g. thinning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
- G06T7/344—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Image Processing (AREA)
Abstract
The invention discloses an intracranial vascular lesion identification method based on transfer learning, which comprises the following steps: acquiring a bright blood image group, a black blood image group and an enhanced black blood image group; preprocessing to obtain a first bright blood image and a first black blood image; taking the corresponding first black blood image as a reference for each first bright blood image, and obtaining a bright blood image group after registration by using a mutual information based on Gaussian distribution sampling and an image pyramid registration method; eliminating flow space artifacts to obtain an artifact elimination enhanced black blood image group; subtracting corresponding images in the artifact removal enhanced black blood image group and the black blood image group to obtain K contrast enhanced images; and establishing a blood three-dimensional model by using the registered bright blood image group and a migration learning method, and establishing a blood boundary expanded blood vessel three-dimensional model by using the registered bright blood image group. The method of the invention can simply, conveniently, quickly and visually identify the focus area of the intracranial blood vessel in clinic.
Description
Technical Field
The invention belongs to the field of image processing, and particularly relates to an intracranial vascular lesion identification method based on transfer learning.
Background
With the rapid development of national economy, people pay more and more attention to health problems. A paper published by Lancet in 2019 in 6 months analyzes the death reasons of residents in 34 provinces (including Hongkong and Australia) in China from 1990 to 2017, and the first death reason of high-centrality Chinese people is found to be stroke. Cerebral apoplexy is a series of symptoms caused by brain tissue necrosis caused by rupture, stenosis or blockage of intracranial blood vessels, including cerebral hemorrhage, cerebral infarction and the like, and if treatment is not timely, a patient can die; even if the treatment is timely, the disability of the patient can be caused.
Currently, for clinically assessing the degree of intracranial vascular lesion and vascular stenosis, a lumen-based imaging method, such as Digital Subtraction Angiography (DSA), CT Angiography (CTA), Magnetic Resonance Angiography (MRA), and High-Resolution Magnetic Resonance Angiography (HRMRA), is commonly used. The intracranial artery blood vessel is connected with the carotid artery and the vertebral artery, and forms a ring structure at the bottom of the brain, and has special structural form, zigzag and extremely thin tube wall thickness. By means of the magnetic resonance blood vessel imaging technology, the path of the intracranial artery blood vessel can be clearly described.
The magnetic resonance angiography (MRA or HRMRA) is used as a non-invasive imaging method for a patient, the vascular wall structure of intracranial vessels can be clearly detected and analyzed, the magnetic resonance image obtained by scanning has high resolution ratio on soft tissues, no bone artifacts and good image quality, and the tissue structures with different imaging characteristics can be obtained by using multiple sequence scanning, so that the magnetic resonance angiography (MRA or HRMRA) has obvious superiority in displaying the intracranial vessels.
Because the images corresponding to the bright blood sequence and the black blood sequence obtained by the magnetic resonance blood vessel imaging technology are two-dimensional images, in clinic, doctors need to obtain the comprehensive condition of the intracranial blood vessels by combining the information of the two images through experience so as to analyze the intracranial vascular lesions. However, the two-dimensional image has limitations, which is not beneficial to simply and rapidly obtaining the real information of the intracranial blood vessel. Moreover, analysis data about the intracranial vascular stenosis degree cannot be intuitively and quickly obtained from the image obtained by the imaging method, and the positioning analysis of the intracranial vascular lesion area is not facilitated clinically.
Disclosure of Invention
In order to be used in clinical application, the real information of the intracranial blood vessel is simply, conveniently and quickly obtained for intracranial vascular lesion analysis. The embodiment of the invention provides an intracranial vascular lesion identification method based on transfer learning. The method comprises the following steps:
acquiring a bright blood image group, a black blood image group and an enhanced black blood image group of an intracranial vascular site; the bright blood image group, the black blood image group and the enhanced black blood image group respectively comprise K bright blood images, black blood images and enhanced black blood images; the images in the bright blood image group, the black blood image group and the enhanced black blood image group are in one-to-one correspondence; k is a natural number greater than 2;
taking each bright blood image and the corresponding enhanced black blood image as an image pair, and preprocessing each image pair to obtain a first bright blood image and a first black blood image of the image pair;
aiming at each first bright blood image, taking the corresponding first black blood image as a reference, and performing image registration by using a mutual information based on Gaussian distribution sampling and image pyramid registration method to obtain a registered bright blood image group comprising K registered bright blood images;
performing flow-space artifact removing operation on the enhanced black blood images in the enhanced black blood image group by using the registered bright blood image group to obtain an artifact-removed enhanced black blood image group comprising K target enhanced black blood images;
subtracting the corresponding image in the artifact removal enhanced black blood image group from the corresponding image in the black blood image group to obtain K contrast enhanced images;
establishing a blood three-dimensional model by using the registered bright blood image group and adopting a transfer learning method;
establishing a blood vessel three-dimensional model of blood boundary expansion by using the registered bright blood image group;
establishing a contrast enhanced three-dimensional model by using the K contrast enhanced images;
obtaining an intracranial angiography enhancement three-dimensional model based on the blood three-dimensional model, the blood vessel three-dimensional model and the angiography enhancement three-dimensional model;
and obtaining the numerical value of a target parameter representing the angiostenosis degree of each segment of blood vessel in the intracranial angiography enhanced three-dimensional model, and marking the intracranial angiography enhanced three-dimensional model by using the numerical value of the target parameter of each segment of blood vessel to obtain an intracranial vascular focus identification model.
In the scheme provided by the embodiment of the invention, firstly, the bright blood image and the enhanced black blood image which are obtained by scanning the magnetic resonance blood vessel imaging technology are subjected to image registration by adopting a mutual information and image pyramid registration method based on Gaussian distribution sampling, so that the registration efficiency can be improved, and the registration accuracy of the images is improved layer by layer from low resolution to high resolution. The bright blood image and the enhanced black blood image can be unified under the same coordinate system through the image registration. And secondly, the registered bright blood image is used for carrying out flow-space artifact elimination operation on the enhanced black blood image, so that more accurate and comprehensive blood vessel information can be displayed. And then, the flow-space artifact is eliminated from the angle of image post-processing without using a new imaging technology, an imaging mode or a pulse sequence, so that the flow-space artifact can be simply, accurately and quickly eliminated, and the method can be better popularized in clinical application. Thirdly, establishing a blood three-dimensional model by using the registered bright blood image and a migration learning method, establishing a blood three-dimensional model of the blood boundary expansion of the blood three-dimensional model of the blood boundary expansion by using the registered bright blood image, and subtracting the artifact removal enhanced black blood image and the black blood image to obtain a contrast enhanced three-dimensional model with a contrast enhancement effect; and finally, obtaining an intracranial angiography enhancement three-dimensional model corresponding to the vascular wall with an angiography enhancement effect based on the blood three-dimensional model, the vascular three-dimensional model and the angiography enhancement three-dimensional model. And finally, marking by using the numerical value of the target parameter for representing the angiostenosis degree in the intracranial angiography enhanced three-dimensional model to obtain the intracranial vascular lesion identification model. The intracranial blood vessel focus identification model realizes three-dimensional visualization of intracranial blood vessels, does not need a doctor to restore an intracranial blood vessel tissue structure, disease characteristics and the like through imagination, can provide vivid intracranial blood vessel three-dimensional space information, is convenient for visual observation, and is convenient for positioning and displaying a narrow focus area. The method can simply, conveniently, quickly and intuitively obtain the real information of the intracranial blood vessel and the analysis data about the intracranial blood vessel stenosis degree in clinical application.
Of course, not all of the advantages described above need to be achieved at the same time in the practice of any one product or method of the invention.
The present invention will be described in further detail with reference to the accompanying drawings and examples.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic flow chart of an intracranial vascular lesion identification method based on transfer learning according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a region to be registered of an intracranial vascular magnetic resonance image in accordance with an embodiment of the invention;
fig. 3(a) is a bright blood gaussian pyramid and a black blood gaussian pyramid of an intracranial vascular magnetic resonance image according to an embodiment of the invention; fig. 3(b) is a bright blood laplacian pyramid and a black blood laplacian pyramid of an intracranial vascular magnetic resonance image according to an embodiment of the present invention;
FIG. 4 is a result of registration of Laplacian pyramid images of intracranial vascular magnetic resonance images according to an embodiment of the invention;
FIG. 5 is a schematic diagram of a Gaussian pyramid image registration step based on mutual information for an intracranial vascular magnetic resonance image according to an embodiment of the invention;
FIG. 6 is a normalized mutual information for different iterations according to an embodiment of the present invention;
FIG. 7 is a registration result of an intracranial vascular magnetic resonance image including a plurality of registration methods of a mutual information pyramid method;
fig. 8 is a result of the registration of the mutual information and image pyramid based on gaussian distribution sampling and the intracranial vascular magnetic resonance image of the mutual information pyramid method in the embodiment of the present invention;
FIG. 9 is a schematic diagram of gray scale linear transformation and parameter setting according to an embodiment of the present invention;
FIG. 10 is a diagram of an image binarization result according to an embodiment of the present invention;
FIG. 11 is a flow-empty artifact removal result obtained by different methods for an intracranial blood vessel according to an embodiment of the present invention;
fig. 12 is an exemplary MIP diagram of an embodiment of the present invention;
FIG. 13 is an inverse graph, a characteristic MIP graph corresponding to the MIP graph of the present invention;
FIG. 14 is an effect diagram of an enhanced three-dimensional model of intracranial angiography of an embodiment of the present invention;
FIG. 15 is a graph showing the effect of the intracranial vascular lesion identification model according to the embodiment of the present invention;
fig. 16 is a diagram showing the effect of an intracranial vascular lesion recognition model and a sectional view according to an embodiment of the present invention.
Detailed Description
In order to be used in clinical application, the real information of the intracranial blood vessel is simply, conveniently and quickly obtained for intracranial vascular lesion analysis. The embodiment of the invention provides an intracranial vascular lesion identification method based on transfer learning.
As shown in fig. 1, fig. 1 is a schematic flow chart of an intracranial vascular lesion identification method based on migration learning according to an embodiment of the present invention, which may include the following steps:
s1, acquiring a bright blood image group, a black blood image group and an enhanced black blood image group of the intracranial vascular site;
the bright blood image group, the black blood image group and the enhanced black blood image group respectively comprise K bright blood images, black blood images and enhanced black blood images; the images in the bright blood image group, the black blood image group and the enhanced black blood image group are in one-to-one correspondence; k is a natural number greater than 2;
the bright blood image group is an image group obtained by performing bright blood sequence scanning on an intracranial vascular site by using a magnetic resonance vascular imaging technology. The black blood image group is an image group obtained by performing black blood sequence scanning on an intracranial vascular site by using a magnetic resonance vascular imaging technology. The enhanced black blood image group is an image group obtained by injecting paramagnetic contrast agent into a patient and then scanning a black blood sequence on an intracranial vascular part by using a magnetic resonance vascular imaging technology. In an embodiment of the invention, the magnetic resonance angiography technique is preferably HRMRA.
The K images in the group of bright blood images, the group of black blood images and the group of enhanced black blood images are in one-to-one correspondence in such a way that the images formed according to the scanning time are in the same order.
S2, taking each bright blood image and the corresponding enhanced black blood image as an image pair, and preprocessing each image pair to obtain a first bright blood image and a first black blood image of the image pair;
this step may be understood as a preprocessing process of the images, and in an alternative embodiment, each image pair is preprocessed to obtain a first bright blood image and a first black blood image of the image pair, which may include S21 and S22:
s21, aiming at each image pair, taking the enhanced black blood image as a reference, carrying out coordinate transformation and image interpolation on the bright blood image, using the similarity measurement based on mutual information, and adopting a preset search strategy to obtain a first bright blood image after pre-registration;
specifically, the enhanced black blood image and the bright blood image are to-be-registered images, and the enhanced black blood image is used as a reference image, the bright blood image is used as a floating image, and the bright blood image is subjected to coordinate transformation according to the orientation tag information in the DICOM file of the bright blood image, so that the purpose of rotating the bright blood image to the same coordinate system as the enhanced black blood image is achieved, and the scanning direction of the rotated bright blood image is also changed into a coronal plane.
And the two search strategies of a gradient descent optimizer and (1+1) -ES are used for respectively registering 160 bright blood images and 160 enhanced black blood images of corresponding scanning layers, wherein the enhanced black blood images are reference images, the bright blood images are floating images,
through the pre-registration of the step, the magnetic resonance images of the same scanning layer can be compared under the same coordinate system preliminarily, but because the scanning time of the bright blood sequence and the scanning time of the black blood sequence are different, and the patient possibly moves slightly before and after the scanning, the operation is only a rough coordinate transformation, the complete registration of the multi-mode magnetic resonance images can not be realized only through the pre-registration, but the step can omit unnecessary processing procedures for the subsequent accurate registration link, and the processing speed is improved.
S22, the same area content as the scanning range of the first bright blood image is extracted from the enhanced black blood image, and a first black blood image is formed.
Optionally, S22 may include the following steps:
1. obtaining edge contour information of a blood vessel in the first bright blood image;
specifically, the edge contour information may be obtained by using a Sobel edge detection method or the like. The edge profile information contains coordinate values of the respective edge points.
2. Extracting the minimum value and the maximum value of the abscissa and the ordinate from the edge profile information, and determining an initial extraction frame based on the obtained four coordinate values;
in other words, in the edge profile information, extracting a minimum abscissa value, a maximum abscissa value, a minimum ordinate value and a maximum ordinate value, and determining four vertexes of the square frame by using the four coordinate values, thereby obtaining an initial extracted frame;
3. in the size range of the first bright blood image, the size of the initial extraction frame is respectively enlarged by a preset number of pixels along four directions to obtain a final extraction frame;
wherein, the four directions are respectively the positive and negative directions of the horizontal and vertical coordinates; the preset number is reasonably selected according to the type of the blood vessel image, so as to ensure that the expanded final extraction frame does not exceed the size range of the first bright blood image, for example, the preset number may be 20.
4. And extracting the corresponding area content in the final extracted frame from the enhanced black blood image to form a first black blood image.
And extracting the content of the corresponding area in the enhanced black blood image according to the coordinate range defined by the final extraction frame, and forming the extracted content into a first black blood image. The step obtains the common scanning range of the magnetic resonance images under the two modes by extracting the region to be registered, thereby being beneficial to subsequent rapid registration.
Referring to fig. 2, fig. 2 is a schematic diagram of a region to be registered of an intracranial vascular magnetic resonance image according to an embodiment of the present invention, where the left image is a first bright blood image after pre-registration, the right image is an enhanced black blood image, and the box is a region to be extracted in the enhanced black blood image. The region contains the common scanning range of a bright blood sequence and a black blood sequence in an intracranial vascular magnetic resonance image, and useful information can be focused more quickly by determining the region to be extracted.
The two-step preprocessing process of the embodiment of the invention plays a very important role, the preprocessed image can pay more attention to useful information and exclude irrelevant information, and in actual use, the image preprocessing can be used for improving the reliability of intracranial blood vessel image registration and identification.
S3, aiming at each first bright blood image, taking the corresponding first black blood image as a reference, and carrying out image registration by using a mutual information based on Gaussian distribution sampling and an image pyramid registration method to obtain a registered bright blood image group comprising K registered bright blood images;
in an alternative embodiment, S3 may include S31-S34:
s31, sampling by adopting Gaussian distribution to select a part of preprocessed image pair as a test image pair;
the test image pair of the embodiment of the invention is the image pair to be registered, and the random selection of the image pair to be registered of the embodiment of the invention adopts Gaussian distribution sampling. The scanning directions of the bright blood image and the enhanced black blood image are different, and in the image preprocessing process, in order to observe the bright blood image and the enhanced black blood image on the same imaging layer, the bright blood image is subjected to coordinate transformation and interpolation, so that each bright blood image corresponds to the enhanced black blood image on the current layer; meanwhile, due to the fact that the scanning ranges of the bright blood image and the enhanced black blood image are different, data of the edge layer of the bright blood image may not be complete. In summary, the bright blood data and the enhanced black blood data of the scanned middle layer are most abundant, so that the gaussian mean μ is selected to be half of the total number of the images to be registered, and the probability of registration of the images in the middle layer selected by gaussian random is the largest.
S32, performing image registration on the first bright blood image and the first black blood image in each test image pair by adopting a registration method based on mutual information and an image pyramid to obtain a rotation matrix corresponding to the first bright blood image in the test image pair after registration; in an optional embodiment, S32 may specifically include steps S321 to S324:
s321, aiming at each test image pair, obtaining a bright blood Gaussian pyramid from the first bright blood image based on downsampling processing, and obtaining a black blood Gaussian pyramid from the first black blood image; the bright blood Gaussian pyramid and the black blood Gaussian pyramid comprise m images with resolution ratios which are sequentially reduced from bottom to top; m is a natural number greater than 3;
in order to improve the accuracy of image registration and avoid the convergence of an image to a local maximum value in the registration process, a multi-resolution strategy can be used for solving the problem of a local extreme value, and meanwhile, the multi-resolution strategy can improve the execution speed of the algorithm and increase the robustness under the condition of meeting the image registration precision. The method is an effective way to improve the registration accuracy and speed by increasing the complexity of the model, namely, in the registration process, the registration is performed in the order from coarse registration to fine registration, firstly, the registration is performed on the low-resolution image, and then, on the basis of the completion of the registration of the low-resolution image, the registration is performed on the high-resolution image.
In an alternative embodiment, the step S321 may include:
obtaining an input image of an ith layer, filtering the input image of the ith layer by using a Gaussian kernel, and deleting even rows and even columns of the filtered image to obtain an image G of the ith layer of the Gaussian pyramidiAnd the ith layer image GiObtaining an i +1 layer image G of a Gaussian pyramid as an i +1 layer input imagei+1;
Wherein i is 1, 2, …, m-1; when the gaussian pyramid is a bright blood gaussian pyramid, the input image of the 1 st layer is a first bright blood image, and when the gaussian pyramid is a black blood gaussian pyramid, the input image of the 1 st layer is a first black blood image.
Specifically, the multiple images in the gaussian pyramid are corresponding to the same original image with different resolutions. The Gaussian pyramid acquires an image through Gaussian filtering and downsampling, and each layer of construction steps can be divided into two steps: firstly, smoothing filtering is carried out on an image by using Gaussian filtering, namely filtering is carried out by using a Gaussian kernel; and then deleting even rows and even columns of the filtered image, namely reducing the width and height of the lower layer image by half to obtain the current layer image, so that the current layer image is one fourth of the size of the lower layer image, and finally obtaining the Gaussian pyramid by continuously iterating the steps.
In this step, the first bright blood image and the first black blood image after the preprocessing are subjected to the processing, so that a bright blood gaussian pyramid and a black blood gaussian pyramid can be obtained. As shown in fig. 3(a), fig. 3(a) is a bright blood gaussian pyramid and a black blood gaussian pyramid of an intracranial vascular magnetic resonance image according to an embodiment of the present invention.
These resolutions are gradually reduced, and the images from the same image combined at different resolutions are arranged to resemble a pyramid, and are therefore referred to as an image pyramid, where the highest resolution image is located at the bottom of the pyramid and the lowest resolution image is located at the top of the pyramid. In the aspect of image information processing, the multi-resolution images can more easily acquire the essential characteristics of the images compared with the traditional single-resolution images.
S322, based on the upsampling processing, obtaining a bright blood Laplacian pyramid by using a bright blood Gaussian pyramid, and obtaining a black blood Laplacian pyramid by using a black blood Gaussian pyramid; wherein the bright blood Laplacian pyramid and the black blood Laplacian pyramid comprise m-1 images with resolution which is sequentially reduced from bottom to top;
since the gaussian pyramid is downsampled, i.e., the image is reduced, a portion of the data of the image is lost. Therefore, in order to avoid data loss of the image in the zooming process and recover detailed data, the embodiment of the invention uses the Laplacian pyramid to realize image reconstruction together with the Gaussian pyramid, and details are highlighted on the basis of the Gaussian pyramid image.
In an alternative embodiment, the step S322 may include:
for the i +1 layer image G of the Gaussian pyramidi+1Performing upsampling, and filling the newly added rows and columns with data 0 to obtain a filled image;
performing convolution on the filling image by utilizing a Gaussian kernel to obtain an approximate value of the filling pixel to obtain an amplified image;
the ith layer image G of the Gaussian pyramidiSubtracting the amplified image to obtain the ith layer image L of the Laplacian pyramidi;
When the gaussian pyramid is the bright blood gaussian pyramid, the laplacian pyramid is the bright blood laplacian pyramid, and when the gaussian pyramid is the black blood laplacian pyramid, the laplacian pyramid is the black blood laplacian pyramid.
Since the laplacian pyramid is a residual between the image and the original image after downsampling, the laplacian pyramid is compared from bottom to top, and the laplacian pyramid has one layer of higher-level image less than the laplacian pyramid structure.
Specifically, the mathematical formula for generating the Laplacian pyramid structure is shown as (1), wherein LiIndicating the Laplacian pyramid (bright blood Laplacian pyramid or black blood Laplacian pyramid) of the i-th layer GiRepresenting the i-th layer of Gaussian pyramid (bright blood Gaussian pyramid or black blood height)Pyramid) and UP operates to UP-sample the enlarged image, symbolIs a sign of the convolution of the symbols,is the gaussian kernel used in constructing the gaussian pyramid.
Corresponding to the gaussian pyramid with 4 layers, the step can obtain a bright blood laplacian pyramid and a black blood laplacian pyramid with 3 image layers. As shown in fig. 3(b), fig. 3(b) is a bright blood laplacian pyramid and a black blood laplacian pyramid of an intracranial vascular magnetic resonance image according to an embodiment of the present invention. The image display uses gamma correction to achieve a clearer effect, and the gamma value is 0.5.
S323, registering images of corresponding layers in the bright blood Laplacian pyramid and the black blood Laplacian pyramid to obtain a registered bright blood Laplacian pyramid;
in an alternative embodiment, the step S323 may include:
aiming at each layer of the bright blood Laplacian pyramid and the black blood Laplacian pyramid, taking the corresponding black blood Laplacian image of the layer as a reference image, taking the corresponding bright blood Laplacian image of the layer as a floating image, using a similarity measure based on mutual information, and adopting a preset search strategy to realize image registration to obtain the registered bright blood Laplacian image of the layer;
forming a registered Laplacian pyramid of the bright blood from bottom to top according to the sequence of the sequential reduction of the resolution by the registered multilayer Laplacian images of the bright blood;
the black blood laplacian image is an image in the black blood laplacian pyramid, and the bright blood laplacian image is an image in the bright blood laplacian pyramid.
The registration process in this step is similar to the pre-registration process, and the image registration is realized by performing coordinate transformation and image interpolation on the bright blood laplacian image and using the similarity measure and the predetermined search strategy, so that the registered bright blood laplacian image can be obtained. The coordinate transformation, the image interpolation, the similarity measurement and the predetermined search strategy are not repeated.
As shown in fig. 4, fig. 4 is a registration result of laplacian pyramid images of an intracranial vascular magnetic resonance image according to an embodiment of the present invention, where the left image is a reference image in a black blood laplacian pyramid, the middle image is a registered image in a bright blood laplacian pyramid, the right image is an effect image obtained by directly superimposing the left and middle images, and the superimposed image displays a montage effect, and the black blood image and the bright blood image are enhanced by using pseudo-color transparency processing, where purple is the enhanced black blood laplacian pyramid image, and green is the bright blood laplacian pyramid image (the image is an image of an original image subjected to gray processing, and the color is not shown).
And S324, registering the images of all layers in the bright blood Gaussian pyramid and the black blood Gaussian pyramid from top to bottom by using the registered bright blood Laplacian pyramid as superposition information to obtain a registered bright blood Gaussian pyramid, and obtaining a rotation matrix corresponding to the first bright blood image in the test image pair after registration.
In the registration of the step, images with different resolutions in the Gaussian pyramid need to be registered, and the registration of the low-resolution images can easily hold essential characteristics of the images, so that the embodiment of the invention registers the high-resolution images on the basis of the registration of the low-resolution images, namely, the images of the Gaussian pyramid are registered from top to bottom, and the registration result of the image of the upper layer is used as the registration input of the image of the lower layer.
In an optional embodiment, in S324, the registering of the blood-brightening gaussian pyramid and the black blood gaussian pyramid from top to bottom is performed by using the registered blood-brightening laplacian pyramid as the overlay information, so as to obtain the registered blood-brightening gaussian pyramid, which may include:
for the j-th layer from top to bottom in the bright blood Gaussian pyramid and the black blood Gaussian pyramid, taking the black blood Gaussian image corresponding to the layer as a reference image, taking the bright blood Gaussian image corresponding to the layer as a floating image, using similarity measurement based on mutual information, and adopting a preset search strategy to realize image registration to obtain a registered j-th layer bright blood Gaussian image;
performing up-sampling operation on the registered jth layer of bright blood Gaussian image, adding the up-sampling operation to the registered corresponding layer of bright blood Laplacian image, and replacing the jth +1 layer of bright blood Gaussian image in the bright blood Gaussian pyramid by using the added image;
taking the black blood Gaussian image of the j +1 th layer as a reference image, taking the replaced bright blood Gaussian image of the j +1 th layer as a floating image, and using a preset similarity measure and a preset search strategy to realize image registration to obtain a registered bright blood Gaussian image of the j +1 th layer;
where j is 1, 2, …, m-1, the black blood gaussian image is an image in the black blood gaussian pyramid, and the bright blood gaussian image is an image in the bright blood gaussian pyramid.
And repeating the operations until the high-resolution registration of the bottom layer Gaussian pyramid image is completed to obtain the registered bright blood Gaussian pyramid. In the registration process, the pre-registration process is similar to the one described above. The coordinate transformation, the image interpolation, the similarity measurement and the predetermined search strategy are not repeated.
The specific steps of mutual information-based gaussian pyramid image registration are shown in fig. 5, and fig. 5 is a schematic diagram of mutual information-based gaussian pyramid image registration steps of an intracranial vascular magnetic resonance image in an embodiment of the present invention. Firstly, registering the low-resolution black blood Gaussian image of the top layer and the low-resolution bright blood Gaussian image of the top layer based on mutual information; then, performing up-sampling operation on the registered bright blood Gaussian image, and adding the up-sampled bright blood Gaussian image and the bright blood Laplacian image of the corresponding layer which retains high-frequency information and is registered according to the operation to be used as a next layer of bright blood Gaussian image; and then, taking the bright blood Gaussian image obtained by the operation as an input image, registering the input image with the black blood Gaussian image of the corresponding layer, and repeating the operation until the high-resolution registration of the bottom layer Gaussian pyramid image is completed.
In the registration of Gaussian pyramid images based on mutual information, the registration of each layer of bright blood Gaussian image and black blood Gaussian image is carried out by taking normalized mutual information as similarity measurement, and the NMI of the two images is calculated through loop iteration until the NMI reaches the maximum. Fig. 6 is normalized mutual information under different iteration times of the embodiment of the present invention, and when the registration of the first-layer image, that is, the bottom-layer image with the highest resolution in the gaussian pyramid reaches the maximum NMI value and the data is stable, the iteration is stopped.
Thus, a registered bright blood Gaussian pyramid is obtained, the coordinate system of the bright blood image is consistent with the coordinate system of the enhanced black blood image, and the images have high similarity, so that the blood vessel image registration process of the embodiment of the invention can be completed. And after registration, a rotation matrix corresponding to the first bright blood image in the test image pair after registration can be obtained.
In order to verify the effectiveness and the practicability of the mutual information and image pyramid-based registration method (referred to as the mutual information pyramid method for short) in the embodiment of the invention, a comparison experiment is also performed, and intracranial vascular magnetic resonance images of five patients are used together, wherein the enhanced black blood image and the bright blood image of the patient A, B, C, D are 160 respectively, and the enhanced black blood image and the bright blood image of the patient E are 150 respectively; meanwhile, an algorithm which only uses DICOM image orientation label information for registration and a registration algorithm based on mutual information measurement are selected and compared with the mutual information pyramid method in the embodiment of the invention, wherein the algorithm based on mutual information measurement is to search the optimal transformation between a reference image and a floating image by a multi-parameter optimization method, so that the mutual information value of the two images is maximum, and the image pyramid algorithm is not used.
The experimental platform was Matlab R2016 b. And combining qualitative analysis and quantitative analysis according to the image registration result of the experiment. In the aspect of qualitative analysis, because large gray scale difference exists between the multi-modal medical images, a difference image obtained by subtracting the registration image from the reference image cannot effectively reflect the registration result of the multi-modal medical images, the embodiment of the invention obtains a color overlapping image capable of reflecting the alignment degree of the registration image and the reference image by overlapping the registration image with the reference image, qualitatively analyzes the registration effect of the multi-modal registration algorithm by the color overlapping image, shows the registration result of the multi-modal intracranial vascular magnetic resonance image in fig. 7, and shows the registration result of the intracranial vascular magnetic resonance image of a plurality of registration methods including a mutual information pyramid method in fig. 7. Wherein, (a) is a reference image; (b) is a floating image; (c) is an overlay image based on image orientation label information; (d) is an overlay image based on a mutual information metric; (e) the invention discloses a superposed image of a mutual information pyramid method. The figures are gray scale images of the original image, not shown in color. In the aspect of quantitative analysis, since the root mean square error RMSE and the peak signal-to-noise ratio PSNR of the evaluation indexes are not suitable for evaluating images with large gray changes, in order to better evaluate the registration result of the multi-modal medical image, the normalized cross-correlation coefficient NCC is adopted, the normalized mutual information NMI is used as the evaluation index, when the values of the normalized cross-correlation coefficient NCC and the normalized mutual information NMI are larger, the higher the image registration accuracy is, and table 1 shows the evaluation index result analysis of different registration algorithms.
Table 1 analysis of the results of different registration methods
aThe value in (1) is the mean value of the evaluation index +/-mean square error based on the registration of a plurality of images of a patient
And (3) qualitative analysis: as is apparent from the overlaid images of fig. 7, the method based on mutual information metric has a large registration shift, and the analysis reason may be that it is easy to fall into a local optimum value rather than a global optimum value only using the method based on mutual information metric; the registration effect based on the image orientation label information is not good enough, and the images are partially not overlapped; the mutual information pyramid method has good image registration effect, the images are displayed more clearly, and the images are almost completely overlapped.
Quantitative analysis: as can be seen from table 1, from the two evaluation indexes NCC and NMI, compared with the registration algorithm using only the orientation tag information of the DICOM image and the registration algorithm based on the mutual information measurement, the mutual information pyramid method according to the embodiment of the present invention is improved in registration accuracy, which shows that the registration method based on the mutual information and the image pyramid according to the embodiment of the present invention can well process the registration of the multi-modal intracranial vascular magnetic resonance image.
S33, obtaining the mean value of the rotation matrix of all the test image pairs;
in the last step, the rotation matrix corresponding to the first bright blood image after registration in each test image pair can be obtained, and then the mean value of the rotation matrices of all the test image pairs can be calculated and obtained.
And S34, performing coordinate transformation on the first bright blood image in the other preprocessed image pairs except the test image pair by using the mean value of the rotation matrix, completing image registration, and obtaining a plurality of registered image pairs.
Considering that when a patient uses a magnetic resonance bright blood sequence scan, if slight movement occurs, the coordinate position of an intracranial blood vessel image obtained by the bright blood sequence scan slightly changes, and then each bright blood image needs to be subjected to a spatial coordinate transformation operation so as to keep the same coordinate position as that of an enhanced black blood image. Mutual information of every two images to be registered needs to be continuously and iteratively calculated in the registration process of the registration method based on the mutual information and the image pyramid, when the size and the number of the images are large, time consumption is overlarge, and if the registration method based on the mutual information and the image pyramid is used for all the images to be registered, the calculation speed is slow. The inventor considers that the intracranial blood vessel can be regarded as a rigid body, which is different from organs such as heart or lung and the like which change along with the movement of human breath and the like, so the space coordinate transformation operation of each bright blood image is almost consistent, namely almost the same rotation matrix is used. Referring to table 2, table 2 shows the mean and mean square error of the rotation matrix calculated by using the mutual information pyramid method, the mean square error of the rotation matrix is very small as can be seen from data, and the rotation matrix is obtained by registering 160 400 × 400 enhanced black blood images and 160 400 × 400 bright blood images derived from patient a by using the mutual information pyramid method. Then, the calculated rotation matrix can be registered by a few bright blood images of the patient to perform the same spatial coordinate transformation on all the bright blood images without the need of solving the rotation matrix for each bright blood image, thereby accelerating the image registration process.
TABLE 2 mean and mean square error of rotation matrix
In the step, the mean value of the rotation matrix, namely a new rotation matrix, is used for carrying out coordinate transformation on the first bright blood image in the rest of the preprocessed image pairs, so that the registration of all images can be completed quickly, and the registration speed is improved greatly. For the processes of implementing coordinate transformation of the bright blood image by using the rotation matrix, image interpolation, using the similarity measurement based on mutual information, and implementing image registration by using a predetermined search strategy, reference may be made to the foregoing, and it may also be understood by combining with the related art, and details are not described herein again.
After the step, a plurality of registered bright blood images can be obtained, and the bright blood images and the corresponding enhanced black blood images are in the same coordinate system and have higher similarity.
In order to verify the feasibility of the registration method based on the mutual information of gaussian distribution sampling and the image pyramid in the embodiment of the present invention, intracranial vascular magnetic resonance images of five patients were used together for registration, wherein the number of the enhanced black blood images and the number of the bright blood images of patient A, B, C, D were 160, and the number of the enhanced black blood images and the number of the bright blood images of patient E were 150, respectively. For the image registration result of an experiment, because the multi-mode magnetic resonance images are different in acquisition principle and different in presented information, a unified gold standard is not available at present to evaluate which registration algorithm is the best, and the quality of the registration result is evaluated according to a specific registration object and an application purpose, so that qualitative analysis and quantitative analysis are combined. In terms of qualitative analysis, the registration algorithm results are qualitatively analyzed by color overlay images that reflect the degree of alignment of the registered image and the reference image, and fig. 8 shows a comparison of the registration results of the multi-modality intracranial vascular magnetic resonance images. Fig. 8 is a result of the registration of the intracranial vascular magnetic resonance image based on the mutual information of gaussian distribution sampling and the image pyramid and the mutual information pyramid method in the embodiment of the present invention. Wherein (a) is a reference image; (b) is a floating image; (c) an overlay image for mutual information pyramid method; (d) for the overlapped images of the method of the invention, the standard deviation sigma is taken as 1; (e) for the overlapped images of the method of the invention, the standard deviation sigma is taken as 2; (f) for the overlapped images of the method of the invention, the standard deviation sigma is taken as 3; (g) for the overlapped images of the method of the invention, the standard deviation sigma is taken as 4; (h) for the overlapped images of the method of the invention, the standard deviation sigma is taken as 5; (i) for the overlaid images of the method of the invention, the standard deviation σ is taken to be 6. Each image is a processed gray scale image, not shown in color. In the aspect of quantitative analysis, a normalized cross-correlation coefficient NCC is adopted, normalized mutual information NMI is used as an evaluation index, and the higher the value of NCC and NMI is, the higher the image registration precision is. Table 3 shows the registration results of all image pairs of one patient using the mutual information and image pyramid-based registration method (abbreviated as mutual information pyramid method), and the mutual information and image pyramid-based registration method based on gaussian distribution sampling (abbreviated as the present invention method) proposed in the present embodiment (the rest of the patient data is not shown by space limitation), and the experimental platform is Matlab R2016 b. As the method of the invention does not need to register all images, only needs to randomly select a few images for registration, the experiment sets the Gaussian distribution mean value mu to be half of the total number of the images to be registered, and randomly extracts 20 enhanced black blood images to register with the corresponding bright blood images when the standard deviation sigma is 1, 2, 3, 4, 5 and 6 respectively.
TABLE 3 analysis of results of different registration algorithms
aThe value in (1) is based on the mean value of the evaluation indexes of the registration of 160 bright blood images and 160 enhanced black blood images +/-mean square error
As is apparent from the overlapped images of fig. 8, the registered images of the mutual information pyramid method and the registered images of the method of the present invention have good performance, and the images are almost completely overlapped together; as can be seen from table 3, from the two evaluation indexes of NCC and NMI, although the accuracy of the method of the present invention is lower than that of the mutual information pyramid method, the difference between them is not too large under different gaussian distribution function settings. When the size of the images participating in registration is large and the number of the images is large, the mutual information pyramid method has a large calculation amount. The method improves the algorithm aiming at the rigid space transformation characteristics of the intracranial blood vessels, accelerates the image registration by utilizing the similarity of the transformation matrix, and the experimental result proves that the time consumed by the improved algorithm registration is only one fifth of the registration time of the mutual information pyramid method, so that the registration speed can be greatly improved, and the registration of the multi-mode intracranial blood vessel magnetic resonance image can be well realized.
In the registration scheme provided by the embodiment of the invention, the intracranial blood vessel can be regarded as a rigid body, and the space coordinate transformation operation of each bright blood image is almost consistent, so that the same rotation matrix can be used. Therefore, on one hand, a small number of image pairs are selected for carrying out image registration of the bright blood images, the average value of the rotation matrix of the registered small number of bright blood images is utilized, the same space coordinate transformation is carried out on the bright blood images of the rest image pairs, and the rotation matrix of each of the rest bright blood images is not required to be solved, so that the image registration process can be accelerated. On the other hand, in the image registration process, mutual information is used as similarity measurement, and an image pyramid algorithm is adopted to increase the complexity of the model, so that the registration accuracy and speed can be improved. Compared with the prior art, the method needs doctors to observe the bright blood image and the black blood image of the intracranial blood vessel according to spatial imagination and subjective experience. The embodiment of the invention adopts the image registration method to unify the bright blood image and the enhanced black blood image in the same coordinate system, so that doctors can conveniently understand the intracranial blood vessel images corresponding to the black blood sequence and the bright blood sequence, the comprehensive information required by diagnosis can be simply, conveniently and quickly obtained, and accurate and reliable reference information can be provided for subsequent medical diagnosis, operation plan formulation, radiotherapy plan and the like. Meanwhile, the image registration is an important step of subsequent flow null artifact elimination. The registration scheme provided by the embodiment of the invention can provide a better reference mode for the registration of other medical images, and has great clinical application value. Meanwhile, the image registration process of the embodiment of the invention is an important basis for eliminating the flow-space artifact subsequently.
After image registration, flow and empty artifacts in the black blood image enhanced after registration can be eliminated, wherein the flow and empty artifacts occur because blood vessels are too small, the blood flow velocity at the tortuous part is slow, and peripheral blood and tissue fluid may have signal pollution and other problems during imaging of blood vessel walls, so that in the image obtained by scanning the black blood sequence, blood information which should be black is instead bright, thereby simulating wall thickening or plaque appearance of normal individuals and exaggerating the degree of blood vessel stenosis. The embodiment of the invention considers that the blood information in the bright blood image after registration is utilized to correct the blood information with incorrect signal display in the enhanced black blood image after registration, and the blood information in the bright blood image after registration is embedded into the enhanced black blood image after registration so as to achieve the effect of image fusion. The method can be realized by the following steps:
s4, carrying out flow-space artifact removing operation on the enhanced black blood image in the enhanced black blood image group by using the registered bright blood image group to obtain an artifact-removed enhanced black blood image group comprising K target enhanced black blood images;
in an alternative embodiment, S4 may include S41-S44:
s41, aiming at each post-registration bright blood image, improving the contrast of the post-registration bright blood image to obtain a contrast enhanced bright blood image;
for example, a specific gray scale linear transformation and parameter setting used for the registered bright blood image is shown in fig. 9, and fig. 9 is a schematic diagram of the gray scale linear transformation and parameter setting provided by the embodiment of the present invention. By using the gray scale linear transformation shown in fig. 9, the smaller gray scale value change interval in the original post-registration bright blood image f can be expanded to the larger gray scale value change interval in the new post-registration bright blood image f1 (contrast enhanced bright blood image), and the image gray scale range can be adjusted to achieve the purpose of improving the contrast of the post-registration bright blood image. By this step, a contrast-enhanced bright blood image can be obtained.
S42, extracting blood information from the contrast enhanced bright blood image to obtain a bright blood characteristic diagram;
in an alternative embodiment, S42 may include the following steps:
s421, determining a first threshold value by using a preset image binarization method;
s422, extracting blood information from the contrast enhanced bright blood image by using a first threshold value;
the method used in this step is called threshold segmentation.
S423, a bright blood feature map is obtained from the extracted blood information.
The preset image binarization method in the embodiment of the invention can comprise a maximum inter-class variance method OTSU, kittle and the like.
The formula for extracting blood information is shown in (2), where T (x, y) is the gray-level value of the contrast-enhanced bright blood image, F (x, y) is the gray-level value of the bright blood feature map, and T is the first threshold.
In an optional implementation manner, the maximum inter-class variance method OTSU is adopted, and the result is shown in fig. 10, where fig. 10 is an image binarization result diagram of the embodiment of the present invention, a left diagram is a contrast enhanced bright blood image, and a right diagram is blood information obtained after the contrast enhanced bright blood image is subjected to threshold extraction. It can be seen that the portion of the right image that appears bright is only blood related information.
S43, carrying out image fusion on the bright blood characteristic image and the enhanced black blood image corresponding to the bright blood image after registration according to a preset fusion formula to obtain a target enhanced black blood image with the flow space artifact eliminated corresponding to the enhanced black blood image;
in the step, firstly, a spatial mapping relation between the bright blood characteristic diagram and the corresponding enhanced black blood image is established, the bright blood characteristic diagram is mapped into the corresponding enhanced black blood image, and image fusion is performed according to a preset fusion formula, wherein the preset fusion formula is as follows:
wherein, F (x, y) is the gray value of the bright blood feature map, R (x, y) is the gray value of the corresponding enhanced black blood image, and g (x, y) is the gray value of the fused target enhanced black blood image.
Through the above operations, the gray value of the flow-space artifact which is supposed to be black but appears as bright color in the corresponding enhanced black blood image can be changed into black, so that the purpose of eliminating the flow-space artifact is achieved.
Referring to fig. 11, fig. 11 is a flow-empty artifact removal result obtained by different methods for intracranial vessels according to an embodiment of the present invention, wherein, the arrow shows the flowing and empty artifact, the left image is the original image of the encephalic blood vessel enhanced black blood image with the flowing and empty artifact, the left two images are the results obtained by the flowing and empty artifact eliminating method based on the mutual information and image pyramid registering method, the left three images are the results obtained by the flowing and empty artifact eliminating method based on the mutual information of Gaussian distribution sampling and the image pyramid registering method, the standard deviation sigma of the gaussian distribution is 3, and it can be seen that the elimination effect of the flow-space artifact elimination method based on the mutual information of the gaussian distribution sampling and the registration method of the image pyramid adopted in the embodiment of the present invention is better than the elimination effect of the flow-space artifact elimination method based on the mutual information and the registration method of the image pyramid. However, it can be understood that the flow-space artifact removing method based on the mutual information of gaussian distribution sampling and the registration method of the image pyramid adopted in the embodiment of the present invention has a faster processing speed than the flow-space artifact removing method based on the mutual information and the pyramid, and experiments prove that about 80% of time can be saved.
And S44, obtaining an artifact-eliminated enhanced black blood image group according to the target enhanced black blood images corresponding to the K enhanced black blood images.
After all the enhanced black blood images are subjected to the flow-space artifact elimination, an artifact eliminated enhanced black blood image group can be obtained.
In the scheme provided by the embodiment of the invention, the bright blood image and the enhanced black blood image obtained by scanning the magnetic resonance blood vessel imaging technology are subjected to image registration based on the mutual information of Gaussian distribution sampling and the registration method of the image pyramid, so that the registration precision and the registration speed can be improved. Blood information is extracted from the registered bright blood image through threshold segmentation and is fused into the registered enhanced black blood image, so that the blood information with incorrect signal display in the registered enhanced black blood image is corrected, and the gray value of the flow-space artifact expressed as bright color is changed into black, so that the aim of eliminating the flow-space artifact is fulfilled, and the intracranial blood vessel image which is displayed more accurately and comprehensively is obtained. The scheme provided by the embodiment of the invention is to eliminate the flow-space artifact from the angle of image post-processing without using a new imaging technology, an imaging mode or a pulse sequence, so that the flow-space artifact can be simply, accurately and quickly eliminated, and the better popularization can be realized in clinical application.
S5, subtracting the corresponding image in the artifact removal enhanced black blood image group and the black blood image group to obtain K contrast enhanced images;
subtracting the corresponding black blood image from each target enhanced black blood image to obtain a contrast enhanced image with a contrast enhanced effect, and subtracting the corresponding black blood image from all the target enhanced black blood images to obtain K contrast enhanced images.
S6, establishing a blood three-dimensional model by using the registered bright blood image group and adopting a transfer learning method;
in an alternative embodiment, S6 may include the following steps:
s61, projecting the registered bright blood image group in three preset directions by using a maximum intensity projection method to obtain MIP (maximum intensity projection) images in all directions;
the Maximum Intensity Projection (MIP) is one of the CT three-dimensional image reconstruction techniques, and is referred to as MIP. The MIP can reflect the X-ray attenuation value of the corresponding pixel, small density change can be displayed on the MIP image, and stenosis, expansion and filling defects of the blood vessel can be well displayed, and calcification on the blood vessel wall and contrast agents in the blood vessel cavity can be well distinguished.
For the MIP method, please refer to the related description of the prior art, which is not repeated herein, refer to fig. 12, and fig. 12 is an exemplary MIP diagram according to an embodiment of the present invention.
And S62, taking the MIP maps in all directions as target domains and the fundus blood vessel map as a source domain, and obtaining two-dimensional blood vessel segmentation maps corresponding to the MIP maps in all directions by using a migration learning method.
The inventor finds out through research that the MIP map of the intracranial vascular bright blood sequence has a distribution of a vascular tree similar to that of the fundus blood vessel. Therefore, the inventor considers that a model pre-trained by a fundus blood vessel (source domain) segmentation task is migrated into an intracranial blood vessel segmentation task by means of a migration learning method, particularly by means of feature migration. Feature based TL (Feature based TL) is to transform the features of the source domain and the target domain into the same space by Feature transformation, assuming that the source domain and the target domain contain some common cross features, so that the source domain data and the target domain data in the space have the same distributed data distribution, and then perform conventional machine learning.
For S62, an alternative embodiment may include S621 to S623:
s621, obtaining a target neural network pre-trained for the eye fundus blood vessel map segmentation task;
the target neural network is obtained by pre-training according to the fundus blood vessel map data set and the improved U-net network model.
As described above, the embodiment of the present invention intends to migrate a pre-trained model of a fundus blood vessel (source domain) segmentation task into an intracranial blood vessel segmentation task by means of a feature migration learning manner. Therefore, it is necessary to obtain a mature network model for the vessel segmentation of the fundus blood vessel map. Specifically, obtaining the target neural network may be performed by the following steps:
in the embodiment of the invention, the structure of the existing U-net network model can be improved, and each sub-module of the U-net network model is respectively replaced by a residual module with a residual connection form, so that the improved U-net network model is obtained. According to the embodiment of the invention, the residual error module is introduced into the U-net network model, so that the problem that the training error does not decrease or inversely increase due to the disappearance of the gradient caused by the deepening of the layer number of the neural network can be effectively solved.
Step 2, obtaining sample data of the fundus blood vessel map;
embodiments of the present invention acquire a fundus angiogram dataset, the DRIVE dataset, which is a dataset that has been labeled.
And 3, training the original network model by using the sample data of the fundus blood vessel map to obtain the trained target neural network.
The following summary describes some parameter characteristics of the target neural network of embodiments of the present invention:
the improved U-net network model in the embodiment of the invention has 5 levels, and a 2.5M parameter ladder network is formed. Each residual module uses 0.25 droout rate (droout means that the neural network unit is temporarily discarded from the network according to a certain probability in the training process of the deep learning network, generally, the droout rate can be set to be 0.3-0.5); and Batch Normalization (BN) is used, the variance size and the mean position are changed by optimization, so that the new distribution is more suitable for the real distribution of data, and the nonlinear expression capability of the model is ensured. The activating function adopts LeakyRelu; the last layer of the network model is activated using Softmax.
Moreover, because of the problem of uneven foreground and background distribution of the medical image sample, the loss function uses a common Dice coefficient (Dice coefficient) loss function for medical image segmentation, and specifically uses an improved Dice loss function, so as to solve the unstable condition of Dice loss function training.
In the aspect of neural network optimization, an Adam optimization algorithm and default parameters are adopted, and the batch size is 256. 250 epochs are trained using the "reduced learning rate" strategy, setting the learning rates at epochs 0, 20, and 150 to 0.01, 0.001, and 0.0001, respectively, and the total learning rate to 250. And the data enhancement is carried out by using a random clipping mode, and the training sample in the DRIVE data set is enlarged by 20000 times.
The process of obtaining the target neural network is briefly introduced, and the trained target neural network can realize the blood vessel segmentation of the fundus blood vessel map to obtain a corresponding two-dimensional blood vessel segmentation map.
S622, carrying out gray inversion processing and contrast enhancement processing on the MIP images in all directions respectively to obtain corresponding characteristic MIP images;
the realization of the feature transfer learning requires that a source domain (fundus blood vessel image) and a target domain (intracranial blood vessel bright blood sequence MIP image) have high similarity and realize the same data distribution.
Therefore, in step S622, the MIP map is subjected to the gradation inversion processing and the contrast enhancement processing, resulting in a characteristic MIP map so that the characteristic MIP map is closer to the fundus blood vessel image.
In an alternative embodiment, S622 may include S6221 and S6222:
s6221, carrying out pixel transformation on the MIP graph by utilizing a gray inversion formula to obtain an inversion graph; wherein, the grayscale inversion formula is t (x) 255-x, x is the pixel value in the MIP map, and t (x) is the pixel value in the inversion map;
the step can be understood in a popular way as grayscale inversion processing, since the pixel range of the MIP map is between 0 and 255, the original brighter region can be darkened and the original darker region can be lightened through the step, specifically, the pixel transformation can be performed through the grayscale inversion formula, the obtained inversion map please refer to the left map in fig. 13, and the left map in fig. 13 is the inversion map corresponding to the MIP map in the embodiment of the present invention.
S6222, the inverse map is used to enhance the contrast by using a contrast-limited adaptive histogram equalization method, so as to obtain a characteristic MIP map.
The main purpose of this step is to enhance the contrast of the inversion map to show a clearer vascularity. As for the way of enhancing the Contrast, any one of the prior arts can be used, and in an alternative embodiment, this step may employ a Contrast Limited Adaptive Histogram Equalization (CLAHE) to enhance the Contrast. For the CLAHE method, reference is made to the prior art for understanding, and no further description is given here. The obtained characteristic MIP map is shown in the right diagram of fig. 13, and the right diagram of fig. 13 is the characteristic MIP map corresponding to the MIP map of the embodiment of the present invention. It can be seen that the contrast of the characteristic MIP map is significantly enhanced and the blood vessels are clearer compared to the inversion map.
After S6222, corresponding characteristic MIP maps can be obtained for the MIP maps in each direction.
In the embodiment of the invention, the cross characteristics of the intracranial blood vessel bright blood sequence MIP and the fundus blood vessel image are considered, so that the MIP image characteristics are mapped to the fundus blood vessel image by adopting a characteristic migration learning method, and the intracranial blood vessel input sample and the fundus blood vessel input sample corresponding to the target neural network have the same sample distribution. Wherein, S621 and S622 may not be in sequence.
S623, respectively inputting the characteristic MIP images in all directions into a target neural network to obtain corresponding two-dimensional blood vessel segmentation images;
and respectively inputting the characteristic MIP images of all directions into a target neural network to obtain a two-dimensional blood vessel segmentation image corresponding to each direction, wherein the obtained two-dimensional blood vessel segmentation image is a binary image, namely pixels are only 0 and 255, white represents a blood vessel, and black represents a background.
S63, synthesizing the two-dimensional vessel segmentation maps in the three directions by using a back projection method to obtain first three-dimensional vessel volume data;
the principle of the back projection method is to evenly distribute measured projection values to each passing point according to the original projection path, back-project the projection values in all directions, and accumulate the back-projected images at all angles to estimate the original image. By synthesizing the two-dimensional vessel segmentation maps in the three directions by using a back projection method, three-dimensional volume data can be obtained, which is referred to as first three-dimensional vessel volume data in the embodiment of the invention. The back projection method in the embodiment of the present invention may be a direct back projection method, a filtered back projection method, a convolution back projection method, and the like, which is not limited herein.
In the embodiment of the present invention, the voxel value of the blood vessel portion in the obtained first three-dimensional blood vessel volume data is 0, and the voxel value of the non-blood vessel portion is minus infinity through the pixel control of the back projection method.
And S64, obtaining an intracranial blood vessel simulation three-dimensional model based on the first three-dimensional blood vessel volume data and the second three-dimensional blood vessel volume data corresponding to the registered bright blood image group.
In an alternative embodiment, S64 may include S641 and S642:
s641, adding the first three-dimensional blood vessel volume data and the second three-dimensional blood vessel volume data to obtain third three-dimensional blood vessel volume data;
the method can be used for directly correspondingly adding each voxel value in the first three-dimensional blood vessel volume data and the second three-dimensional blood vessel volume data to obtain third three-dimensional blood vessel volume data, and cerebrospinal fluid and fat signals with the same intracranial and blood vessel signal intensity can be eliminated through the step.
And S642, processing the third three-dimensional blood vessel volume data by using a threshold segmentation method to obtain an intracranial blood vessel simulation three-dimensional model.
The threshold segmentation method is an image segmentation technology based on regions, and the principle is to divide image pixel points into a plurality of classes. The purpose of image thresholding is to divide the set of pixels by gray level, each resulting subset forming a region corresponding to the real scene, each region having consistent properties within it, while adjacent regions do not have such consistent properties. Threshold segmentation is a method for processing an image into a high-contrast, easily recognizable image with a proper pixel value as a boundary.
The threshold segmentation method adopted by the embodiment of the invention comprises a maximum inter-class variance method, a maximum entropy, an iteration method, a self-adaptive threshold, a manual method, an iteration method, a basic global threshold method and the like. In an alternative implementation manner, the embodiment of the present invention may adopt a maximum inter-class variance method.
The maximum inter-class variance method (or referred to as "Otsu" for short) is a method for automatically calculating a threshold value suitable for a bimodal situation, and performing S642 by using OTSU may include the following steps:
firstly, calculating a first threshold corresponding to centered fourth three-dimensional blood vessel volume data in third three-dimensional blood vessel volume data by using the OTSU;
in this step, one threshold corresponding to a plurality of images in one small cube (referred to as fourth three-dimensional blood vessel volume data) located near the middle of the large three-dimensional cube of the third three-dimensional blood vessel volume data is determined as a first threshold by using the OTSU method. Because the blood information is substantially concentrated in the middle of the image in the third three-dimensional blood vessel volume data, the small cube data (fourth three-dimensional blood vessel volume data) in the middle is selected to determine the first threshold value in the third three-dimensional blood vessel volume data, so that the calculation amount of the threshold value can be reduced, the calculation speed can be improved, and the first threshold value can be accurately applied to all the blood information in the third three-dimensional blood vessel volume data.
For the size of the fourth three-dimensional blood vessel volume data, the central point of the third three-dimensional blood vessel volume data can be determined firstly, and then the preset side length extends in six directions corresponding to the cube, so that the size of the fourth three-dimensional blood vessel volume data is determined; the preset side length may be determined according to an empirical value including a Willis ring, such as 1/4 that is the side length of the cube of the third three-dimensional blood vessel volume data. The Willis loop is the most important collateral circulation pathway in the cranium, linking the bilateral hemisphere with the anterior and posterior circulation.
And then, threshold segmentation of the third three-dimensional blood vessel volume data is realized by utilizing the first threshold, and an intracranial blood vessel simulation three-dimensional model is obtained.
It can be understood by those skilled in the art that, by threshold segmentation, the gray-scale value of the point on the image corresponding to the third three-dimensional blood vessel volume data can be set to 0 or 255, that is, the whole image exhibits a distinct black-and-white effect, the blood information is highlighted as white, and the irrelevant information is displayed as black. For the processing procedure of threshold segmentation, please refer to the prior art, and will not be described herein. And finally, obtaining the intracranial blood vessel simulation three-dimensional model.
S7, establishing a blood vessel three-dimensional model of blood boundary expansion by using the registered bright blood image group;
the step S6 is to obtain a three-dimensional model of blood, which represents the flow direction and regional distribution of intracranial blood, and because there is a blood vessel wall at the periphery of blood in practice, the three-dimensional model of blood cannot represent the actual intracranial blood vessel condition completely.
Therefore, in step S7, the blood boundary in the registered bright blood image may be expanded so that the expanded blood boundary can cover the range of the intracranial blood vessel wall to form the effect of a hollow tube, and then the three-dimensional model is generated by the three-dimensional reconstruction method for the two-dimensional image after the blood boundary is expanded, so as to obtain the three-dimensional model of the blood vessel closer to the real intracranial blood vessel condition than the three-dimensional model of the blood in step S6.
In an alternative embodiment, S7 may include S71-S75:
s71, obtaining K bright blood characteristic graphs;
namely, the K bright blood feature maps obtained in step S42 are obtained.
S72, expanding the boundary of the blood in each bright blood characteristic map by utilizing an expansion operation to obtain an expanded bright blood characteristic map corresponding to the bright blood characteristic map;
the expansion operation is one of morphological operations, and the expansion operation can fill the holes in the image and expand the protruding points of the object at the edge outwards, so that the final expanded object has a larger area than the original one. The dilation operation may be denoted as a ≧ B, defined as a ≦ B ≦ x ≦ B (x) λ a ≠ Φ, where B is a structural element and a is the original. The original image a here is a bright blood feature map, in which only two pixel values, 0 and 255, are present, 0 corresponding to black and 255 corresponding to white.
S73, obtaining a difference characteristic diagram corresponding to the bright blood characteristic diagram by subtracting the expanded bright blood characteristic diagram corresponding to the bright blood characteristic diagram from the bright blood characteristic diagram;
the difference feature map obtained by this step for each bright blood feature map is a two-dimensional plan similar to a hollow blood vessel. Similarly, the pixel values of the difference feature map are only 0 and 255.
S74, determining a third threshold;
this step may select a pixel value as the third threshold value for all difference feature maps according to empirical values, for example, any value between 100 and 200, for example, 128, may be selected as the third threshold value.
And S75, taking the third threshold as an input threshold of the moving cube method, and processing the K difference feature maps by using the moving cube method to obtain the blood vessel three-dimensional model with the blood boundary expanded.
The moving cube method uses the third threshold as an input threshold, and a blood vessel three-dimensional model of blood boundary expansion can be obtained from the K difference feature maps. The specific implementation process of the method for moving cubes is not described herein.
S8, establishing a contrast enhanced three-dimensional model by using the K contrast enhanced images;
this step can be implemented by using a moving cube method, see S6 and S7, which are not described herein.
And S9, obtaining an intracranial angiography enhancement three-dimensional model based on the blood three-dimensional model, the blood vessel three-dimensional model and the angiography enhancement three-dimensional model.
In an alternative embodiment, S9 may include the following steps:
s91, reserving the overlapped part of the contrast enhanced three-dimensional model and the blood vessel three-dimensional model to obtain a reserved contrast enhanced three-dimensional model;
since the contrast-enhanced three-dimensional model obtained in S8 does not only include contrast enhancement of blood vessels, but also needs to exclude enhancement characteristics of unrelated tissues, the search range of the vascular wall contrast enhancement characteristics in the vascular three-dimensional model obtained in S7 is used to determine whether the contrast-enhanced three-dimensional model obtained in S8 is located in a vascular wall region near blood, that is, whether there is an overlapping portion with the vascular three-dimensional model in the contrast-enhanced three-dimensional model, and if so, it indicates that the overlapping portion is located within the search range, and the overlapping portion needs to be retained, so that the retained contrast-enhanced three-dimensional model is obtained.
And S92, fusing the reserved contrast enhanced three-dimensional model with the blood three-dimensional model to obtain the intracranial angiography enhanced three-dimensional model.
The reserved contrast enhanced three-dimensional model representing angiography enhancement is fused with the blood three-dimensional model representing blood information, so that the blood vessel wall with obvious contrast enhancement can be visually displayed, the contrast enhancement effect in which part range of the intracranial blood vessel is most obvious can be clearly seen, and then atherosclerosis or vulnerable plaques possibly appear in the region.
In an optional embodiment, a contrast enhancement quantitative analysis may be obtained in the intracranial angiography enhancement three-dimensional model, and specifically, a plaque enhancement index CE may be obtained for any one point on a blood vessel wall in the intracranial angiography enhancement three-dimensional model, where CE is defined as:
wherein S ispreBBMRAnd SpostBBMRSignal intensity in the black blood image and the contrast enhanced black blood image, respectively.
As will be understood by those skilled in the art, SpreBBMRAnd SpostBBMRThe information carried in the images after the black blood image and the contrast enhanced black blood image are taken, respectively. The embodiment of the invention obtains the plaque enhancement index CE of each point of the edge of the intracranial vascular wall by using the information, and the plaque enhancement index CE is embodied in the intracranial angiography enhanced three-dimensional model, so that a doctor can conveniently obtain more detailed vascular information, and particularly, when the CE is more than a plaque threshold value, such as 0.5, the plaque enhancement index CE indicates that plaque appears on the vascular wall, so that the measurement of the plaque enhancement index CE of the vascular wall area is helpful for identifying responsible intracranialArterial plaque, etc., may provide valuable diagnostic aid information. The fusion technique of the two three-dimensional models can be implemented by using the prior art, and is not described herein. Specific results referring to fig. 14, fig. 14 is a diagram illustrating the effect of the intracranial angiography enhanced three-dimensional model according to the embodiment of the present invention. Wherein the image is grey scale processed. In practice, different colors can be used for the distinction in fig. 14, for example, blue is the blood vessel part where no contrast enhancement appears, and red is the blood vessel part where contrast enhancement appears. In the attached drawings of the specification, the bright-colored part in the white coil is an intracranial vascular part with contrast enhancement, namely, intracranial atherosclerosis diseases or vulnerable plaques possibly appear at the part, the other part is a vascular part without contrast enhancement, and the angiography enhancement three-dimensional model can realize basic functions of rotation, amplification and reduction and the like, so that a doctor is assisted to position a focus area and make a more accurate judgment.
S10, obtaining the value of the target parameter representing the angiostenosis degree of each segment of blood vessel in the angiography enhanced three-dimensional model, and marking the angiography enhanced three-dimensional model by using the value of the target parameter of each segment of blood vessel to obtain the intracranial vascular lesion recognition model.
In an alternative embodiment, S10 may include S101 to S104:
s101, segmenting each section of blood vessel in the angiography enhanced three-dimensional model from three preset directions to obtain two-dimensional sectional views of each direction;
in this step, the blood vessels in the angiography enhanced three-dimensional model may be divided, and for each segment of blood vessel, the segmentation may be performed from three preset orientations to obtain two-dimensional sectional views in each orientation.
Wherein, three preset positions include: axial, coronal, and sagittal.
S102, carrying out corrosion operation on the blood vessel in the two-dimensional sectional diagram of each direction, and recording the target corrosion times when the blood vessel is corroded to a single pixel;
the corrosion operation is one of morphological operations, the corrosion operation can eliminate edge data of an object, the corroded object has a smaller area than the original area and even completely disappears, and corrosion can break some small and long communication areas. When the blood vessel is thick, a plurality of corrosion operations can be carried out, and when the blood vessel is thin, only a few corrosion operations can be carried out.
In step S102, performing erosion operation on the blood vessel in the axial two-dimensional sectional view, and recording the target erosion times n corresponding to the erosion of the blood vessel in the axial two-dimensional sectional view to a single pixel1(ii) a Carrying out corrosion operation on the blood vessel in the two-dimensional sectional diagram of the coronal position, and recording the corresponding target corrosion times n when the blood vessel in the two-dimensional sectional diagram of the azimuth corrodes to a single pixel2(ii) a Carrying out corrosion operation on the blood vessel in the two-dimensional sectional diagram of the sagittal position, and recording the corresponding target corrosion times n when the blood vessel in the two-dimensional sectional diagram of the azimuth corrodes to a single pixel3。
S103, obtaining a numerical value of a target parameter representing the stenosis degree of the section of the blood vessel according to the target corrosion times of the section of the blood vessel in the three directions respectively;
in an alternative embodiment, the target parameter includes stenosis rate and/or flatness; those skilled in the art will appreciate that both of these parameters may be indicative of the degree of vascular stenosis.
When the target parameter includes a stenosis rate, S103 may include:
according to n1、n2、n3Obtaining the value of the stenosis rate of the section of blood vessel by using a stenosis rate formula of the blood vessel; wherein, the stenosis rate formula is:
wherein, the resolution is the resolution of each azimuth two-dimensional sectional image (the resolution of the three azimuth two-dimensional sectional images is the same), and the smaller the numerical value of the stenosis rate is, the narrower the blood vessel is.
When the target parameter includes flatness, S103 may include:
according to n1、n2、n3Obtaining the value of the flatness of the section of the blood vessel by using a blood vessel flatness formula; wherein, the flatness formula is as follows:
a larger value of the degree of flattening indicates a narrower vessel.
And S104, marking the angiography enhanced three-dimensional model by using the numerical value of the target parameter of each section of blood vessel to obtain an intracranial vascular lesion identification model.
Through the steps, the numerical value of the target parameter of each segment of blood vessel can be obtained, and then the numerical values of each segment of blood vessel can be marked on the angiography enhanced three-dimensional model to obtain the intracranial blood vessel focus identification model. The numerical value of the target parameter of each point is embedded into the intracranial vascular lesion identification model, so that the numerical value of the target parameter of each point can be extracted and displayed when needed, a doctor can conveniently obtain the data of the vascular stenosis degree of each position in time when observing the overall three-dimensional vascular state, for example, when the intracranial vascular lesion identification model is displayed on a display screen of a computer, the numerical value of the stenosis rate and/or the flatness of the mouse position point can be displayed in a blank area of the model.
For convenience of visual display, different numerical values can be marked on the angiography enhanced three-dimensional model by different colors to obtain the angiography enhanced three-dimensional stenosis analysis model, for example, multiple colors from light to dark can be correspondingly marked for stenosis rate numerical values from small to large, and for flatness numerical values, because the numerical values are fewer and only 2 numerical values are possible, two color corresponding marks distinguished from the stenosis rate can be adopted. The narrowing degree of the blood vessel can be more intuitively shown by adopting the color display of different tones, so that the attention of a doctor can be attracted.
FIG. 15 is a graph showing the effect of an angiographic enhanced three-dimensional stenosis analysis model of intracranial vessels, in accordance with an embodiment of the present invention. Wherein the left graph is the stenosis rate marking effect and the right graph is the flatness marking effect. In practice, different colors are displayed on the model, so that the degree of narrowing can be distinguished, for example, a thinner part of a blood vessel is warm, the narrowest part is red, a thicker part of the blood vessel is cool, the thickest part is green, and the like, a white arrow indicates abrupt narrowing of the intracranial blood vessel, and color display with different colors can more intuitively show the narrowing of the blood vessel. In the figure are the effects of the grey scale processing, the colours not being shown.
Furthermore, since doctors are used to observe the two-dimensional medical images of the tangent plane, the embodiment of the invention can provide the simulated three-dimensional vascular stenosis analysis model and simultaneously provide the two-dimensional tangent plane images of three directions, i.e. the images of the coronal plane, the sagittal plane and the axial plane where the current point corresponding to each point in the simulated three-dimensional vascular stenosis analysis model is located can be displayed. Referring to fig. 16, fig. 16 is a diagram showing the effect of an angiographic enhanced three-dimensional stenosis analysis model and a sectional view of an intracranial vessel according to an embodiment of the present invention. In fig. 16, there may be a blood vessel narrowing at the warm tone of the blood vessel, there is no obvious blood vessel narrowing at the cold tone, and the three two-dimensional images on the right side of the image are respectively imaged from top to bottom on the axial plane, the sagittal plane and the coronal plane where the current point is located; when the simulated three-dimensional vascular stenosis analysis model is displayed, the functions of measuring the distance by two points and measuring the angle by three points can be realized by using the points with three colors such as red, green and blue, the three points are displayed on the left lower side of the display screen, and the volume size of the currently selected model is displayed on the right lower side of the display screen. So that the doctor can obtain more detailed data of the intracranial blood vessel.
In the scheme provided by the embodiment of the invention, firstly, the bright blood image and the enhanced black blood image which are obtained by scanning the magnetic resonance blood vessel imaging technology are subjected to image registration by adopting a mutual information and image pyramid registration method based on Gaussian distribution sampling, so that the registration efficiency can be improved, and the registration accuracy of the images is improved layer by layer from low resolution to high resolution. The bright blood image and the enhanced black blood image can be unified under the same coordinate system through the image registration. And secondly, the registered bright blood image is used for carrying out flow-space artifact elimination operation on the enhanced black blood image, so that more accurate and comprehensive blood vessel information can be displayed. And then, the flow-space artifact is eliminated from the angle of image post-processing without using a new imaging technology, an imaging mode or a pulse sequence, so that the flow-space artifact can be simply, accurately and quickly eliminated, and the method can be better popularized in clinical application. Thirdly, establishing a blood three-dimensional model by using the registered bright blood image and a migration learning method, establishing a blood vessel three-dimensional model with blood boundary expansion by using the registered bright blood image, and subtracting the artifact-removed enhanced black blood image and the black blood image to obtain a contrast enhanced three-dimensional model with a contrast enhancement effect; and finally, obtaining an intracranial angiography enhancement three-dimensional model corresponding to the vascular wall with an angiography enhancement effect based on the blood three-dimensional model, the vascular three-dimensional model and the angiography enhancement three-dimensional model. And finally, marking by using the numerical value of the target parameter for representing the angiostenosis degree in the intracranial angiography enhanced three-dimensional model to obtain an intracranial vascular lesion identification model. The intracranial blood vessel focus identification model realizes three-dimensional visualization of intracranial blood vessels, does not need a doctor to restore an intracranial blood vessel tissue structure, disease characteristics and the like through imagination, can provide vivid intracranial blood vessel three-dimensional space information, is convenient for visual observation, and is convenient for positioning and displaying a narrow focus area. The method can simply, conveniently, quickly and intuitively obtain the real information of the intracranial blood vessel and the analysis data about the intracranial blood vessel stenosis degree in clinical application.
Claims (10)
1. An intracranial vascular lesion identification method based on transfer learning is characterized by comprising the following steps:
acquiring a bright blood image group, a black blood image group and an enhanced black blood image group of an intracranial vascular site; the bright blood image group, the black blood image group and the enhanced black blood image group respectively comprise K bright blood images, black blood images and enhanced black blood images; the images in the bright blood image group, the black blood image group and the enhanced black blood image group are in one-to-one correspondence; k is a natural number greater than 2;
taking each bright blood image and the corresponding enhanced black blood image as an image pair, and preprocessing each image pair to obtain a first bright blood image and a first black blood image of the image pair;
aiming at each first bright blood image, taking the corresponding first black blood image as a reference, and performing image registration by using a mutual information based on Gaussian distribution sampling and image pyramid registration method to obtain a registered bright blood image group comprising K registered bright blood images;
performing flow-space artifact removing operation on the enhanced black blood images in the enhanced black blood image group by using the registered bright blood image group to obtain an artifact-removed enhanced black blood image group comprising K target enhanced black blood images;
subtracting the corresponding image in the artifact removal enhanced black blood image group from the corresponding image in the black blood image group to obtain K contrast enhanced images;
establishing a blood three-dimensional model by using the registered bright blood image group and adopting a transfer learning method;
establishing a blood vessel three-dimensional model of blood boundary expansion by using the registered bright blood image group;
establishing a contrast enhanced three-dimensional model by using the K contrast enhanced images;
obtaining an intracranial angiography enhancement three-dimensional model based on the blood three-dimensional model, the blood vessel three-dimensional model and the angiography enhancement three-dimensional model;
and obtaining the numerical value of a target parameter representing the angiostenosis degree of each segment of blood vessel in the intracranial angiography enhanced three-dimensional model, and marking the intracranial angiography enhanced three-dimensional model by using the numerical value of the target parameter of each segment of blood vessel to obtain an intracranial vascular focus identification model.
2. The method of claim 1, wherein preprocessing each image pair to obtain a first bright blood image and a first black blood image of the image pair comprises:
for each image pair, taking the enhanced black blood image as a reference, performing coordinate transformation and image interpolation on the bright blood image, using similarity measurement based on mutual information, and adopting a preset search strategy to obtain a pre-registered first bright blood image;
and extracting the same area content as the scanning range of the first bright blood image from the enhanced black blood image to form a first black blood image.
3. The method according to claim 1 or 2, wherein the performing, for each first bright blood image, image registration by using a mutual information based on gaussian distribution sampling and image pyramid registration method with reference to the corresponding first black blood image to obtain a registered bright blood image group including K registered bright blood images comprises:
adopting Gaussian distribution sampling to select a part of preprocessed image pair as a test image pair;
performing image registration on the first bright blood image and the first black blood image in each test image pair by adopting a registration method based on mutual information and an image pyramid to obtain a rotation matrix corresponding to the first bright blood image in the test image pair after registration;
obtaining the mean value of the rotation matrix of all the test image pairs;
and performing coordinate transformation on the first bright blood image in the other preprocessed image pairs except the test image pair by using the mean value of the rotation matrix to complete image registration to obtain a registered bright blood image group comprising K registered bright blood images.
4. The method according to claim 3, wherein the performing image registration on the first bright blood image and the first black blood image in each test image pair by using a registration method based on mutual information and an image pyramid to obtain a rotation matrix corresponding to the first bright blood image in the test image pair after registration comprises:
aiming at each test image pair, based on down-sampling processing, obtaining a bright blood Gaussian pyramid from the first bright blood image, and obtaining a black blood Gaussian pyramid from the first black blood image; the bright blood Gaussian pyramid and the black blood Gaussian pyramid comprise m images with resolution becoming smaller in sequence from bottom to top; m is a natural number greater than 3;
based on the upsampling processing, obtaining a bright blood Laplacian pyramid by using the bright blood Gaussian pyramid, and obtaining a black blood Laplacian pyramid by using the black blood Gaussian pyramid; the bright blood Laplacian pyramid and the black blood Laplacian pyramid comprise m-1 images with resolution which is sequentially reduced from bottom to top;
registering images of corresponding layers in the bright blood Laplacian pyramid and the black blood Laplacian pyramid to obtain a registered bright blood Laplacian pyramid;
and registering the images of all layers in the bright blood Gaussian pyramid and the black blood Gaussian pyramid from top to bottom by using the registered bright blood Laplacian pyramid as superposition information to obtain a registered bright blood Gaussian pyramid, and obtaining a rotation matrix corresponding to the first bright blood image in the registered test image pair.
5. The method according to claim 4, wherein the registering the images of the respective layers in the blood-brightening Gaussian pyramid and the black blood Gaussian pyramid from top to bottom by using the registered blood-brightening Gaussian pyramid as overlay information to obtain the registered blood-brightening Gaussian pyramid comprises:
for the j-th layer from top to bottom in the bright blood Gaussian pyramid and the black blood Gaussian pyramid, taking the black blood Gaussian image corresponding to the layer as a reference image, taking the bright blood Gaussian image corresponding to the layer as a floating image, using similarity measurement based on mutual information, and adopting a preset search strategy to realize image registration to obtain a registered j-th layer bright blood Gaussian image;
performing upsampling operation on the registered jth layer of bright blood Gaussian image, adding the upsampled operation to the registered corresponding layer of bright blood Laplacian image, and replacing the jth +1 layer of bright blood Gaussian image in the bright blood Gaussian pyramid by using the added image;
taking the black blood Gaussian image of the j +1 th layer as a reference image, taking the replaced bright blood Gaussian image of the j +1 th layer as a floating image, and using a preset similarity measure and a preset search strategy to realize image registration to obtain a registered bright blood Gaussian image of the j +1 th layer;
wherein j is 1, 2, …, m-1, the black blood gaussian image is an image in the black blood gaussian pyramid, and the bright blood gaussian image is an image in the bright blood gaussian pyramid.
6. The method according to claim 1 or 5, wherein said performing an empty artifact removing operation on the enhanced black blood image in the enhanced black blood image group by using the post-registration bright blood image group to obtain an artifact removed enhanced black blood image group including K target enhanced black blood images comprises:
for each post-registration bright blood image, improving the contrast of the post-registration bright blood image to obtain a contrast enhanced bright blood image;
extracting blood information from the contrast enhanced bright blood image to obtain a bright blood characteristic diagram;
carrying out image fusion on the bright blood characteristic graph and the enhanced black blood image corresponding to the registered bright blood image according to a preset fusion formula to obtain a target enhanced black blood image with the air artifact removed corresponding to the enhanced black blood image;
and enhancing the black blood image by using the targets corresponding to the K enhanced black blood images to obtain an artifact-eliminated enhanced black blood image group.
7. The method of claim 6, wherein extracting blood information from the contrast enhanced bright blood image to obtain a bright blood feature map comprises:
determining a first threshold value by using a preset image binarization method;
extracting blood information from the contrast-enhanced bright blood image using the first threshold;
and obtaining a bright blood characteristic map from the extracted blood information.
8. The method according to claim 1 or 7, wherein the establishing of the blood three-dimensional model by using the registered bright blood image group and a migration learning method comprises:
projecting the registered bright blood image group in three preset directions by using a maximum intensity projection method to obtain MIP (maximum intensity projection) images in all directions;
taking the MIP images in all directions as target domains and the fundus blood vessel images as source domains, and obtaining two-dimensional blood vessel segmentation images corresponding to the MIP images in all directions by using a migration learning method;
synthesizing the two-dimensional vessel segmentation maps in the three directions by using a back projection method to obtain first three-dimensional vessel volume data; wherein the voxel value of the blood vessel part in the first three-dimensional blood vessel volume data is 0, and the voxel value of the non-blood vessel part is minus infinity;
and obtaining a blood three-dimensional model based on the first three-dimensional blood vessel volume data and the second three-dimensional blood vessel volume data corresponding to the registered bright blood image group.
9. The method according to claim 8, wherein the MIP maps of all directions are used as target domains, the fundus blood vessel map is used as a source domain, and a migration learning method is used for obtaining two-dimensional blood vessel segmentation maps corresponding to the MIP maps of all directions; the method comprises the following steps:
obtaining a pre-trained target neural network aiming at the eye fundus blood vessel map segmentation task; the target neural network is obtained by pre-training according to the fundus blood vessel map data set and the improved U-net network model;
respectively carrying out gray level inversion processing and contrast enhancement processing on the MIP images in all directions to obtain corresponding characteristic MIP images; wherein the characteristic MIP map has the same sample distribution as the fundus blood vessel map;
and respectively inputting the characteristic MIP maps of all directions into the target neural network to obtain corresponding two-dimensional vessel segmentation maps.
10. The method according to claim 1 or 9, wherein obtaining values of target parameters representing the degree of stenosis of blood vessels of each segment of blood vessel in the intracranial angiography-enhanced three-dimensional model, and marking the intracranial angiography-enhanced three-dimensional model by using the values of the target parameters of each segment of blood vessel to obtain an intracranial vascular lesion recognition model comprises:
cutting each segment of blood vessel in the intracranial angiography enhanced three-dimensional model from three preset directions to obtain a two-dimensional sectional view of each direction;
carrying out corrosion operation on the blood vessel in the two-dimensional sectional diagram of each direction, and recording the target corrosion times when the blood vessel is corroded to a single pixel;
obtaining a numerical value of a target parameter representing the stenosis degree of the section of the blood vessel according to the target corrosion times of the section of the blood vessel in the three directions respectively;
and marking the intracranial angiography enhanced three-dimensional model by using the numerical value of the target parameter of each section of blood vessel to obtain an intracranial vascular lesion identification model.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011324149.7A CN112669257A (en) | 2020-11-23 | 2020-11-23 | Intracranial vascular lesion identification method based on transfer learning |
CN202111381503.4A CN114170151A (en) | 2020-11-23 | 2021-11-22 | Intracranial vascular lesion identification method based on transfer learning |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011324149.7A CN112669257A (en) | 2020-11-23 | 2020-11-23 | Intracranial vascular lesion identification method based on transfer learning |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112669257A true CN112669257A (en) | 2021-04-16 |
Family
ID=75403511
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011324149.7A Withdrawn CN112669257A (en) | 2020-11-23 | 2020-11-23 | Intracranial vascular lesion identification method based on transfer learning |
CN202111381503.4A Pending CN114170151A (en) | 2020-11-23 | 2021-11-22 | Intracranial vascular lesion identification method based on transfer learning |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111381503.4A Pending CN114170151A (en) | 2020-11-23 | 2021-11-22 | Intracranial vascular lesion identification method based on transfer learning |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN112669257A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022105623A1 (en) * | 2020-11-23 | 2022-05-27 | 西安科锐盛创新科技有限公司 | Intracranial vascular focus recognition method based on transfer learning |
CN115307582A (en) * | 2022-10-11 | 2022-11-08 | 南通电博士自动化设备有限公司 | Method for measuring sectional area of wire and cable |
-
2020
- 2020-11-23 CN CN202011324149.7A patent/CN112669257A/en not_active Withdrawn
-
2021
- 2021-11-22 CN CN202111381503.4A patent/CN114170151A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022105623A1 (en) * | 2020-11-23 | 2022-05-27 | 西安科锐盛创新科技有限公司 | Intracranial vascular focus recognition method based on transfer learning |
CN115307582A (en) * | 2022-10-11 | 2022-11-08 | 南通电博士自动化设备有限公司 | Method for measuring sectional area of wire and cable |
Also Published As
Publication number | Publication date |
---|---|
CN114170151A (en) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107909585B (en) | Intravascular intima segmentation method of intravascular ultrasonic image | |
WO2022105623A1 (en) | Intracranial vascular focus recognition method based on transfer learning | |
WO2022105647A1 (en) | Method for establishing enhanced three-dimensional model of intracranial angiography | |
CN112598619A (en) | Method for establishing intracranial vascular simulation three-dimensional narrowing model based on transfer learning | |
US11830193B2 (en) | Recognition method of intracranial vascular lesions based on transfer learning | |
CN112562058B (en) | Method for quickly establishing intracranial vascular simulation three-dimensional model based on transfer learning | |
CN112509075A (en) | Intracranial vascular lesion marking and three-dimensional display method based on intelligent medical treatment | |
CN112508873A (en) | Method for establishing intracranial vascular simulation three-dimensional narrowing model based on transfer learning | |
CN110648338B (en) | Image segmentation method, readable storage medium, and image processing apparatus | |
CN112634196A (en) | Medical image segmentation and display method based on intelligent medical treatment | |
CN114170152A (en) | Method for establishing simulated three-dimensional intracranial vascular stenosis analysis model | |
CN114170151A (en) | Intracranial vascular lesion identification method based on transfer learning | |
CN112509079A (en) | Method for establishing intracranial angiography enhanced three-dimensional narrowing analysis model | |
CN112509076A (en) | Intracranial vascular lesion marking and three-dimensional display system based on intelligent medical treatment | |
CN114170337A (en) | Method for establishing intracranial vascular enhancement three-dimensional model based on transfer learning | |
CN112669256B (en) | Medical image segmentation and display method based on transfer learning | |
CN112509077A (en) | Intracranial blood vessel image segmentation and display method based on intelligent medical treatment | |
CN112508874A (en) | Cerebrovascular lesion marking and three-dimensional display system based on intelligent medical treatment | |
CN112509080A (en) | Method for establishing intracranial vascular simulation three-dimensional model based on transfer learning | |
US20220164967A1 (en) | Method of establishing an enhanced three-dimensional model of intracranial angiography | |
CN114240841A (en) | Establishment method of simulated three-dimensional vascular stenosis analysis model | |
CN112634386A (en) | Method for establishing angiography enhanced three-dimensional narrowing analysis model | |
CN112509081A (en) | Method for establishing intracranial angiography enhanced three-dimensional narrowing analysis model | |
CN112669439B (en) | Method for establishing intracranial angiography enhanced three-dimensional model based on transfer learning | |
CN112508881A (en) | Intracranial blood vessel image registration method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20210416 |
|
WW01 | Invention patent application withdrawn after publication |