CN110956107A - Three-dimensional blood vessel type distinguishing method based on OCT (optical coherence tomography) scanning system - Google Patents

Three-dimensional blood vessel type distinguishing method based on OCT (optical coherence tomography) scanning system Download PDF

Info

Publication number
CN110956107A
CN110956107A CN201911145891.9A CN201911145891A CN110956107A CN 110956107 A CN110956107 A CN 110956107A CN 201911145891 A CN201911145891 A CN 201911145891A CN 110956107 A CN110956107 A CN 110956107A
Authority
CN
China
Prior art keywords
matrix
blood vessel
dimensional
diameter
scanning system
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911145891.9A
Other languages
Chinese (zh)
Other versions
CN110956107B (en
Inventor
韩定安
刘碧旺
廖锤
曾锟
张章
曾亚光
王雪花
王茗祎
熊红莲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Foshan University
Original Assignee
Foshan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Foshan University filed Critical Foshan University
Priority to CN201911145891.9A priority Critical patent/CN110956107B/en
Publication of CN110956107A publication Critical patent/CN110956107A/en
Application granted granted Critical
Publication of CN110956107B publication Critical patent/CN110956107B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/26Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
    • G06V10/267Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/30Noise filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03Recognition of patterns in medical or anatomical images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/14Vascular patterns
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses a three-dimensional blood vessel type distinguishing method based on an OCT (optical coherence tomography) scanning system, which comprises the steps of obtaining a three-dimensional stereo image group of a blood vessel in the time direction through the OCT scanning system; extracting a two-dimensional image corresponding to unit depth from the three-dimensional stereo image group, converting the two-dimensional image into a two-dimensional matrix, and forming the two-dimensional matrix and time corresponding to the two-dimensional matrix into a three-dimensional matrix T; carrying out noise reduction processing on the three-dimensional matrix T to obtain a matrix G; performing Fourier transform on the matrix G to convert the matrix G from a time domain space to a frequency domain space to obtain a matrix K; carrying out quantization processing on the matrix K through a parameter IR to obtain a matrix A, wherein an image represented by the matrix A is an artery distribution map; subtracting the matrix A from the matrix T to obtain a matrix B; performing convolution on the matrix B by taking a Gaussian function as a convolution kernel to obtain a matrix C; determining a vessel radial search direction; and obtaining the diameter of the blood vessel and obtaining a vein distribution diagram and a capillary distribution diagram by searching the radial search direction of the blood vessel. The method is mainly used for biological imaging.

Description

Three-dimensional blood vessel type distinguishing method based on OCT (optical coherence tomography) scanning system
Technical Field
The invention relates to the technical field of medical detection, in particular to a three-dimensional blood vessel type distinguishing method based on an OCT (optical coherence tomography) scanning system.
Background
Before methods for automatically distinguishing blood vessel types are proposed, the blood vessel types are mainly judged by doctors or biological laboratory technicians according to self experiences. With the development of optical imaging technology, some methods for automatically distinguishing blood vessel types based on optical imaging are also continuously proposed. In 1982 and 1998, Akita et al and Jusephsky-Haw Yu et al set threshold values to distinguish the type of retinal vessels in the fundus according to the difference of the reflection light intensity of arteriovenous vessels. However, in some cases of single-wavelength illumination, the light intensity of more artery and vein regions is similar, which is not enough to distinguish the blood vessel types according to the light intensity difference. In 2007, zhanhoboshi et al, university of texas a & M, used a photoacoustic microscope to observe changes in hemoglobin oxygen saturation (SO2), and discriminated arteriovenous types from measured blood oxygen saturation. In 2007, Indian scientist Narasimhaiyer et al identified blood vessel types by structural and functional features of fundus retinal images at 570nm and 600nm wavelengths. Furthermore, single wavelength based imaging methods are also used to distinguish vessel types. In 2010 and 2011, the brain cortex blood vessel type recognition is realized by the blood vessel anatomical features and the light intensity distribution in single-wavelength LSCI (SW-LSCI) images of the Shanghai traffic university nerve engineering laboratory, Roento doctor and the like and the Huazhong science and technology university Von Neng cloud doctor in sequence. With the application and development of MRA and CT techniques, they are also used for three-dimensional vessel classification to distinguish head, lung and heart.
At present, two-dimensional blood vessel classification technologies mainly comprise a multi-wavelength artery and vein classification method and a single-wavelength blood vessel classification method, but the two methods are complex and are not beneficial to rapidly distinguishing the types of blood vessels.
Disclosure of Invention
The present invention is directed to a method for distinguishing three-dimensional blood vessel types based on an OCT scanning system, which solves one or more of the problems of the prior art, and provides at least one of the advantages.
The solution of the invention for solving the technical problem is as follows: a three-dimensional blood vessel type distinguishing method based on an OCT scanning system comprises the following steps:
step 1, performing continuous C scanning on a region where a blood vessel is located through an OCT (optical coherence tomography) scanning system to obtain a three-dimensional map group of the blood vessel in the time direction;
step 2, extracting a two-dimensional image corresponding to unit depth from the three-dimensional stereo image group, converting the two-dimensional image into a two-dimensional matrix, and forming the two-dimensional matrix and time corresponding to the two-dimensional matrix into a three-dimensional matrix T;
step 3, carrying out noise reduction processing on the three-dimensional matrix T to obtain a matrix G;
step 4, carrying out Fourier transform on the matrix G to convert the matrix G from a time domain space to a frequency domain space to obtain a matrix K;
step 5, carrying out quantization processing on the matrix K through a parameter IR to obtain a matrix A, wherein an image represented by the matrix A is an artery distribution map;
6, subtracting the matrix A from the matrix T to obtain a matrix B;
step 7, performing convolution on the matrix B by taking a Gaussian function as a convolution kernel to obtain a matrix C;
step 8, calculating a Hessian matrix H (x, y) of the matrix C;
step 9, calculating eigenvalues and corresponding eigenvectors of a Hessian matrix H (x, y), taking the included angle direction of the two eigenvectors as the direction of the blood vessel, and rotating by 90 degrees to obtain the radial search direction of the blood vessel;
step 10, obtaining the diameter of the blood vessel by searching the radial search direction of the blood vessel;
step 11, obtaining a vein distribution map and a capillary distribution map through the diameter of a blood vessel;
wherein the expression of the parameter IR is:
Figure BDA0002282184730000031
Figure BDA0002282184730000032
is given as I (upsilon) at a specific frequency f0Intensity value of (v) of0) Is the intensity value of I (upsilon) at zero frequency, and is expressed as the intensity value of each point of the matrix K, and the specific frequency f0Is preset.
Further, in step 3, the noise reduction processing is eight-domain noise reduction processing.
Further, the eight-domain denoising specifically comprises: and (3) taking 3-by-3 windows with all values of 1 to convolute the matrix T in turn, and smoothing data by using Savitzky-Golay filtering operators for the eight neighborhoods in a single column in sequence along the time direction.
Further, in step 10, the method for determining the diameter of the blood vessel includes: performing blood vessel segmentation, then performing binarization, respectively starting to search blood vessel boundaries along the radial search direction and the opposite direction of the blood vessel, wherein the distance between adjacent blood vessel boundaries is the diameter of the blood vessel, and the diameter of the blood vessel is recorded as a first diameter; calculating gradient value change conditions along the radial search direction of the blood vessel and the opposite direction of the radial search direction of the blood vessel by adopting a gradient extreme value method, determining the edge of the blood vessel, wherein the gradient value of the edge of the blood vessel is increased from zero to gradient, the initial end and the termination end of the mutation are taken as two edge points of the blood vessel, the Euclidean distance between the two edge points is taken as the diameter of the blood vessel, and the diameter of the blood vessel is taken as a second diameter; the first diameter and the second diameter are averaged to obtain the vessel diameter.
Further, in step 11, the acquisition of the vein profile and the capillary profile includes: marking the blood vessels with the diameter of more than 15 mu m to obtain a vein distribution map, and marking the blood vessels with the diameter of less than or equal to 15 mu m to obtain a capillary distribution map.
The invention has the beneficial effects that: the three-dimensional map group is obtained by using an OCT scanning system, then the artery is separated out by using the difference of vibration signals, and then the vein and the capillary vessel are separated according to the difference between the diameters of the vein and the capillary vessel, so that the three-dimensional blood vessel types are distinguished with high precision.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings used in the description of the embodiments will be briefly described below. It is obvious that the described drawings are only a part of the embodiments of the invention, not all embodiments, and that a person skilled in the art will be able to derive other designs and drawings from these drawings without the exercise of inventive effort.
Fig. 1 is a flow chart of steps of a three-dimensional blood vessel type distinguishing method based on an OCT scanning system.
Detailed Description
Reference will now be made in detail to the present preferred embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to like elements throughout.
In the description of the present invention, it should be understood that the orientation or positional relationship referred to in the description of the orientation, such as up, down, front, rear, left, right, etc., is the orientation or positional relationship shown in the drawings, and is only for convenience of description and simplification of the description of the present invention, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus, should not be construed as limiting the present invention.
In the description of the invention, if words such as "a number" or the like are used, the meaning is one or more, the meaning of a plurality is two or more, more than, less than, more than, etc. are understood as not including the number, and more than, less than, more than, etc. are understood as including the number.
In the description of the present invention, unless otherwise explicitly defined, terms such as setup, installation, connection, and the like should be understood in a broad sense, and those skilled in the art can reasonably determine the specific meanings of the terms in the present invention in combination with the detailed contents of the technical solutions.
Embodiment 1, referring to fig. 1, a three-dimensional blood vessel type distinguishing method based on an OCT scanning system includes:
step 1, performing continuous C scanning on a region where a blood vessel is located through an OCT (optical coherence tomography) scanning system to obtain a three-dimensional map group of the blood vessel in the time direction;
step 2, extracting a two-dimensional image corresponding to unit depth from the three-dimensional stereo image group, converting the two-dimensional image into a two-dimensional matrix, and forming the two-dimensional matrix and time corresponding to the two-dimensional matrix into a three-dimensional matrix T;
step 3, carrying out noise reduction processing on the three-dimensional matrix T to obtain a matrix G;
step 4, carrying out Fourier transform on the matrix G to convert the matrix G from a time domain space to a frequency domain space to obtain a matrix K;
step 5, carrying out quantization processing on the matrix K through a parameter IR to obtain a matrix A, wherein an image represented by the matrix A is an artery distribution map;
6, subtracting the matrix A from the matrix T to obtain a matrix B;
step 7, performing convolution on the matrix B by taking a Gaussian function as a convolution kernel to obtain a matrix C;
step 8, calculating a Hessian matrix H (x, y) of the matrix C;
step 9, calculating eigenvalues and corresponding eigenvectors of a Hessian matrix H (x, y), taking the included angle direction of the two eigenvectors as the direction of the blood vessel, and rotating by 90 degrees to obtain the radial search direction of the blood vessel;
step 10, obtaining the diameter of the blood vessel by searching the radial search direction of the blood vessel;
step 11, obtaining a vein distribution map and a capillary distribution map through the diameter of a blood vessel;
wherein the expression of the parameter IR is:
Figure BDA0002282184730000061
Figure BDA0002282184730000062
is given as I (upsilon) at a specific frequency f0Intensity value of (v) of0) Is the intensity value of I (upsilon) at zero frequency, and is expressed as the intensity value of each point of the matrix K, and the specific frequency f0Is preset.
Wherein, in some optimized embodiments, in step 3, the denoising process is an eight-domain denoising process. The eight-field denoising treatment specifically comprises the following steps: and (3) taking 3-by-3 windows with all values of 1 to convolute the matrix T in turn, and smoothing data by using Savitzky-Golay filtering operators for the eight neighborhoods in a single column in sequence along the time direction. The Savitzky-Golay filter is a special low-pass filter, can be used for directly carrying out smooth denoising on data in a time domain, and has the characteristic of high fidelity of abnormal signals. The idea of this filter is to slide filter the data with a window width of 2m +1, and fit an nth order polynomial to the data within the window:
Figure BDA0002282184730000071
calculating expression of fitting polynomial P (x) based on least square principle, and substituting the horizontal coordinate of the central point in the window into P (x) to obtain function value P (x)0) The Savitzky-Golay filtering value at the point is sequentially and continuously slid from head to tail to complete the filtering of all data.
Since the vibration of the artery originates from the beating of the animal's heart, the I (upsilon) of the artery location is mainly distributed at a specific frequency f0(ii) a The veins and capillaries have no obvious and stable vibration signals, so that the I (upsilon) of the veins and capillaries are randomly distributed in the whole frequency domain. Based on the characteristics, the vibration noise at non-artery positions can be weakened through the parameter IR, so thatThe artery position is highlighted. Thereby obtaining an arterial map.
After obtaining the artery distribution map, it is necessary to distinguish veins and capillaries, and a gaussian function is used as a convolution kernel to reduce the interference of random noise, so as to obtain a matrix G, which is represented by T (x, y; s), that is:
T(x,y;s)=T(x,y)*G(x,y;s)
wherein, is convolution symbol, variance is s2Is defined as:
Figure BDA0002282184730000072
in order to accurately calculate the blood vessel diameter, it is necessary to calculate the direction of the blood vessel first, and to use the vertical direction of the blood vessel direction as the search direction in the calculation of the blood vessel diameter. The invention adopts eigenvalue analysis based on Hessian matrix to calculate the direction of the blood vessel. The Hessian matrix of a two-dimensional image is a symmetric square matrix composed of second-order partial derivatives of a bivariate real-valued function, and is defined as:
Figure BDA0002282184730000081
respectively convolving the original image T (X, Y) by calculating first and second derivatives of the Gaussian kernel function in the X direction and the Y direction, and finally obtaining a Hessian matrix of the matrix C at the point (X, Y):
Figure BDA0002282184730000082
and calculating the eigenvalue of the Hessian matrix H (x, y) and the corresponding eigenvector, taking the included angle direction of the two eigenvectors as the direction of the blood vessel, and rotating by 90 degrees to obtain the radial search direction of the blood vessel. And obtaining the diameter of the blood vessel by searching the radial search direction of the blood vessel. The specific method comprises the following steps: performing blood vessel segmentation, then performing binarization, respectively starting to search blood vessel boundaries along the radial search direction and the opposite direction of the blood vessel, wherein the distance between adjacent blood vessel boundaries is the diameter of the blood vessel, and the diameter of the blood vessel is recorded as a first diameter; calculating gradient value change conditions along the radial search direction of the blood vessel and the opposite direction of the radial search direction of the blood vessel by adopting a gradient extreme value method, determining the edge of the blood vessel, wherein the gradient value of the edge of the blood vessel is increased from zero to gradient, the initial end and the termination end of the mutation are taken as two edge points of the blood vessel, the Euclidean distance between the two edge points is taken as the diameter of the blood vessel, and the diameter of the blood vessel is taken as a second diameter; the first diameter and the second diameter are averaged to obtain the vessel diameter. The method for segmenting the blood vessel includes but is not limited to: frangi algorithm based on Hessian matrix, algorithm based on PCA, matched filtering algorithm and adaptive contrast enhancement algorithm. By averaging the first diameter and the second diameter, the diameter of the blood vessel at a certain position can be determined more accurately, so that veins and capillaries can be better distinguished. Wherein the method of distinguishing between veins and capillaries is by comparing the diameters of veins and capillaries. Specifically, blood vessels with the diameter of more than 15 mu m are marked to obtain a vein distribution map, and blood vessels with the diameter of less than or equal to 15 mu m are marked to obtain a capillary distribution map.
The invention utilizes an OCT scanning system to obtain a three-dimensional map group, then separates out the artery by utilizing the difference of vibration signals, and separates the vein and the capillary vessel by the difference of the diameters of the vein and the capillary vessel, thereby realizing the high-precision distinguishing of the three-dimensional blood vessel types.
While the preferred embodiments of the present invention have been described in detail, it will be understood by those skilled in the art that the invention is not limited to the details of the embodiments shown, but is capable of various modifications and substitutions without departing from the spirit of the invention.

Claims (5)

1. A three-dimensional blood vessel type distinguishing method based on an OCT scanning system is characterized by comprising the following steps:
step 1, performing continuous C scanning on a region where a blood vessel is located through an OCT (optical coherence tomography) scanning system to obtain a three-dimensional map group of the blood vessel in the time direction;
step 2, extracting a two-dimensional image corresponding to unit depth from the three-dimensional stereo image group, converting the two-dimensional image into a two-dimensional matrix, and forming the two-dimensional matrix and time corresponding to the two-dimensional matrix into a three-dimensional matrix T;
step 3, carrying out noise reduction processing on the three-dimensional matrix T to obtain a matrix G;
step 4, carrying out Fourier transform on the matrix G to convert the matrix G from a time domain space to a frequency domain space to obtain a matrix K;
step 5, carrying out quantization processing on the matrix K through a parameter IR to obtain a matrix A, wherein an image represented by the matrix A is an artery distribution map;
6, subtracting the matrix A from the matrix T to obtain a matrix B;
step 7, performing convolution on the matrix B by taking a Gaussian function as a convolution kernel to obtain a matrix C;
step 8, calculating a Hessian matrix H (x, y) of the matrix C;
step 9, calculating eigenvalues and corresponding eigenvectors of a Hessian matrix H (x, y), taking the included angle direction of the two eigenvectors as the direction of the blood vessel, and rotating by 90 degrees to obtain the radial search direction of the blood vessel;
step 10, obtaining the diameter of the blood vessel by searching the radial search direction of the blood vessel;
step 11, obtaining a vein distribution map and a capillary distribution map through the diameter of a blood vessel;
wherein the expression of the parameter IR is:
Figure FDA0002282184720000021
Figure FDA0002282184720000022
is given as I (upsilon) at a specific frequency f0Intensity value of (v) of0) Is the intensity value of I (upsilon) at zero frequency, and is expressed as the intensity value of each point of the matrix K, and the specific frequency f0Is preset.
2. The three-dimensional blood vessel type distinguishing method based on the OCT scanning system of claim 1, wherein in step 3, the noise reduction processing is eight-domain noise reduction processing.
3. The three-dimensional blood vessel type distinguishing method based on the OCT scanning system of claim 2, wherein the eight-domain noise reduction processing specifically comprises: and (3) taking 3-by-3 windows with all values of 1 to convolute the matrix T in turn, and smoothing data by using Savitzky-Golay filtering operators for the eight neighborhoods in a single column in sequence along the time direction.
4. The three-dimensional blood vessel type distinguishing method based on the OCT scanning system of claim 1, wherein in step 10, the method for determining the diameter of the blood vessel comprises: performing blood vessel segmentation, then performing binarization, respectively starting to search blood vessel boundaries along the radial search direction and the opposite direction of the blood vessel, wherein the distance between adjacent blood vessel boundaries is the diameter of the blood vessel, and the diameter of the blood vessel is recorded as a first diameter; calculating gradient value change conditions along the radial search direction of the blood vessel and the opposite direction of the radial search direction of the blood vessel by adopting a gradient extreme value method, determining the edge of the blood vessel, wherein the gradient value of the edge of the blood vessel is increased from zero to gradient, the initial end and the termination end of the mutation are taken as two edge points of the blood vessel, the Euclidean distance between the two edge points is taken as the diameter of the blood vessel, and the diameter of the blood vessel is taken as a second diameter; the first diameter and the second diameter are averaged to obtain the vessel diameter.
5. The method for distinguishing three-dimensional blood vessel types based on the OCT scanning system of claim 1, wherein in step 11, the obtaining manner of the vein profile and the capillary vessel profile comprises: marking the blood vessels with the diameter of more than 15 mu m to obtain a vein distribution map, and marking the blood vessels with the diameter of less than or equal to 15 mu m to obtain a capillary distribution map.
CN201911145891.9A 2019-11-21 2019-11-21 Three-dimensional blood vessel type distinguishing method based on OCT scanning system Active CN110956107B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911145891.9A CN110956107B (en) 2019-11-21 2019-11-21 Three-dimensional blood vessel type distinguishing method based on OCT scanning system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911145891.9A CN110956107B (en) 2019-11-21 2019-11-21 Three-dimensional blood vessel type distinguishing method based on OCT scanning system

Publications (2)

Publication Number Publication Date
CN110956107A true CN110956107A (en) 2020-04-03
CN110956107B CN110956107B (en) 2023-05-30

Family

ID=69977864

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911145891.9A Active CN110956107B (en) 2019-11-21 2019-11-21 Three-dimensional blood vessel type distinguishing method based on OCT scanning system

Country Status (1)

Country Link
CN (1) CN110956107B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111681242A (en) * 2020-08-14 2020-09-18 北京至真互联网技术有限公司 Retinal vessel arteriovenous distinguishing method, device and equipment

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105188550A (en) * 2013-03-12 2015-12-23 光学实验室成像公司 Vascular data processing and image registration systems, methods, and apparatuses
CN107468234A (en) * 2017-08-28 2017-12-15 佛山科学技术学院 A kind of vascular flow projection chromatography 3-D imaging system and imaging method
CN107485383A (en) * 2017-09-29 2017-12-19 佛山科学技术学院 A kind of speckle blood flow imaging method and apparatus based on constituent analysis
CN107595250A (en) * 2017-09-30 2018-01-19 浙江大学 The blood flow imaging method and system of contrast is mixed with figure based on motion
US20180020909A1 (en) * 2016-07-20 2018-01-25 Oregon Health & Science University Automated quantification of nonperfusion in the retina using optical coherence tomography angiography
CN108670239A (en) * 2018-05-22 2018-10-19 浙江大学 A kind of the three-dimensional flow imaging method and system in feature based space
CN109907731A (en) * 2019-01-31 2019-06-21 浙江大学 The three-dimensional flow angiographic method and system of optical coherence tomography based on feature space
CN110160440A (en) * 2019-03-07 2019-08-23 佛山科学技术学院 A kind of three-dimensional colour dynamic imaging device and method based on frequency domain OCT technology

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105188550A (en) * 2013-03-12 2015-12-23 光学实验室成像公司 Vascular data processing and image registration systems, methods, and apparatuses
US20180020909A1 (en) * 2016-07-20 2018-01-25 Oregon Health & Science University Automated quantification of nonperfusion in the retina using optical coherence tomography angiography
CN107468234A (en) * 2017-08-28 2017-12-15 佛山科学技术学院 A kind of vascular flow projection chromatography 3-D imaging system and imaging method
CN107485383A (en) * 2017-09-29 2017-12-19 佛山科学技术学院 A kind of speckle blood flow imaging method and apparatus based on constituent analysis
CN107595250A (en) * 2017-09-30 2018-01-19 浙江大学 The blood flow imaging method and system of contrast is mixed with figure based on motion
CN108670239A (en) * 2018-05-22 2018-10-19 浙江大学 A kind of the three-dimensional flow imaging method and system in feature based space
CN109907731A (en) * 2019-01-31 2019-06-21 浙江大学 The three-dimensional flow angiographic method and system of optical coherence tomography based on feature space
CN110160440A (en) * 2019-03-07 2019-08-23 佛山科学技术学院 A kind of three-dimensional colour dynamic imaging device and method based on frequency domain OCT technology

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
FULI ZHANG 等: "In vivo full-field functional optical hemocytometer", 《J. BIOPHOTONICS.》 *
RIWEI LIAO 等: "High-temporal-resolution, full-field optical angiography based on short time modulation depth for vascular occlusion tests", 《JOURNAL OF BIOMEDICAL OPTICS》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111681242A (en) * 2020-08-14 2020-09-18 北京至真互联网技术有限公司 Retinal vessel arteriovenous distinguishing method, device and equipment

Also Published As

Publication number Publication date
CN110956107B (en) 2023-05-30

Similar Documents

Publication Publication Date Title
Akbar et al. Arteriovenous ratio and papilledema based hybrid decision support system for detection and grading of hypertensive retinopathy
Wu et al. On the adaptive detection of blood vessels in retinal images
Zhao et al. Automated vessel segmentation using infinite perimeter active contour model with hybrid region information with application to retinal images
CN111340789A (en) Method, device, equipment and storage medium for identifying and quantifying eye fundus retinal blood vessels
US9940722B2 (en) Segmentation and identification of closed-contour features in images using graph theory and quasi-polar transform
Rossant et al. A morphological approach for vessel segmentation in eye fundus images, with quantitative evaluation
Moghimirad et al. Retinal vessel segmentation using a multi-scale medialness function
ES2921974T3 (en) Image processing method for the detection of glaucoma and associated software products
KR102095723B1 (en) Improvements in or relating to image processing
Khan et al. A generalized multi-scale line-detection method to boost retinal vessel segmentation sensitivity
Kafieh et al. An accurate multimodal 3-D vessel segmentation method based on brightness variations on OCT layers and curvelet domain fundus image analysis
Rehman et al. Microscopic retinal blood vessels detection and segmentation using support vector machine and K‐nearest neighbors
Amini et al. Classification of medical image modeling methods: A review
Liu et al. A multiscale ray-shooting model for termination detection of tree-like structures in biomedical images
Mukherjee et al. Tree2tree2: Neuron tracing in 3d
CN110956107A (en) Three-dimensional blood vessel type distinguishing method based on OCT (optical coherence tomography) scanning system
Morales et al. Segmentation and analysis of retinal vascular tree from fundus images processing
KR101601041B1 (en) Pattern Classification Apparatus and Method for fMRI
Bhattacharya et al. A new approach to automated retinal vessel segmentation using multiscale analysis
Li et al. Automated retinal vessel segmentation using multiscale analysis and adaptive thresholding
KR102223263B1 (en) Method and Apparatus for Interpretation of Magnifying Endoscopy Images of Esophageal Cancer using Computer Aided Diagnostic System
Singh et al. Automatic detection of hard exudates shadow region within retinal layers of OCT images
Hatanaka et al. Artery and vein diameter ratio measurement based on improvement of arteries and veins segmentation on retinal images
DelPozo‐Baños et al. Localized component filtering for electroencephalogram artifact rejection
CN112017132A (en) Vein image enhancement method based on maximum curvature method and multi-scale Hessian matrix

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant