CN112164046B - Automatic analysis method for ocular microvascular hemodynamic parameters - Google Patents

Automatic analysis method for ocular microvascular hemodynamic parameters Download PDF

Info

Publication number
CN112164046B
CN112164046B CN202011017868.4A CN202011017868A CN112164046B CN 112164046 B CN112164046 B CN 112164046B CN 202011017868 A CN202011017868 A CN 202011017868A CN 112164046 B CN112164046 B CN 112164046B
Authority
CN
China
Prior art keywords
image
ocular
microvasculature
diameter
central line
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011017868.4A
Other languages
Chinese (zh)
Other versions
CN112164046A (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.)
Meishi Optical Technology Guangdong Co ltd
Original Assignee
Sun Yat Sen University
Zhongshan Ophthalmic Center
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 Sun Yat Sen University, Zhongshan Ophthalmic Center filed Critical Sun Yat Sen University
Priority to CN202011017868.4A priority Critical patent/CN112164046B/en
Publication of CN112164046A publication Critical patent/CN112164046A/en
Application granted granted Critical
Publication of CN112164046B publication Critical patent/CN112164046B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration by the use of histogram techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30041Eye; Retina; Ophthalmic
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Abstract

The invention provides an automatic analysis method for ocular microvascular hemodynamics, which comprises the following steps: step 1: inputting an image set I which needs to carry out analysis on ocular microvascular hemodynamic parameters in a video file containing ocular microvascular information o (ii) a Step 2: for image set I o Carrying out image registration on the continuous images; and step 3: carrying out image segmentation on the registered image to obtain a segmentation mask image set S of the ocular microvasculature r (ii) a And 4, step 4: obtaining a segmented mask image set S r The central line of the middle eye microvasculature; and 5: fitting the central line, extracting the edge of the ocular microvasculature in the segmentation mask image, and calculating the average diameter, the inner diameter and the outer diameter of the ocular microvasculature; and 6: forming a spatiotemporal image; and 7: carrying out histogram equalization processing on the obtained space-time image; and step 8: carrying out linear detection processing on the spatio-temporal image, carrying out statistics on the slope of the obtained straight line, and obtaining the slope k with the most dense distribution; and step 9: meterCalculating blood flow velocity V of ocular microvasculature a (ii) a Step 10: and calculating blood flow and tube wall shear rate.

Description

Automatic analysis method for ocular microvascular hemodynamic parameters
Technical Field
The invention relates to the field of parameter analysis of ocular microvasculature, in particular to an automatic analysis method for ocular microvasculature hemodynamic parameters.
Background
The eyeball and the brain share the blood supply of the internal carotid artery, and the function of the ocular microvasculature is a potential index for reflecting ocular surface inflammation and systemic microvascular pathological diseases, so that the method can provide important basis for researching the pathogenesis of cerebrovascular, cardiovascular and diabetes and other diseases, disease progress, judging the treatment effect and the like. With the rapid development of ophthalmic imaging technology in recent years, the emergence and development of new imaging technologies such as functional imaging OCT and ocular microvascular functional analysis provide important information for diagnosis and treatment of ophthalmic diseases, and the morphological measurement function of the imaging technology is applied to the precise analysis of eyeball tissue functions, which is called a milestone result of ophthalmic diagnostic technology. The existing eye imaging research technology has the defects that the matching analysis software is incomplete, and the processing and intelligent analysis of batch images can not meet the scientific research and clinical requirements and the like.
The existing system for collecting the ocular microvasculature is realized by integrating a high-speed camera, a video capture card and other digital image collecting devices into a standard slit lamp platform, and the slight shake of a patient can generate larger offset in an image in the collecting process, so that the image which can be used for blood flow analysis in the whole video image only accounts for less than 20% of the whole video image, the existing system adopts manual labeling of doctors, and the workload is very huge; in addition, blood vessels displayed in the eye microvascular image are discontinuous, and the situation that the blood vessels are discontinuous exists in a single image, so that the blood vessels acquired in the single image by adopting the traditional method for analysis are very limited, and the usability of the blood vessels is greatly reduced; the subsequent analysis of the vascular mechanics in the prior art is completed by a manual method, which causes different errors due to different individuals or fatigue degrees and consumes time.
Disclosure of Invention
The invention mainly aims to provide an automatic analysis method for ocular microvascular hemodynamic parameters, which realizes full-automatic hemodynamic analysis and greatly improves the speed of the hemodynamic analysis and the accuracy of core parameter calculation.
In order to achieve the above purposes, the technical scheme adopted by the invention is as follows: an analysis method of blood flow velocity of ocular microvasculature, comprising the steps of:
step 1: inputting an image set I which needs to carry out analysis on ocular microvascular hemodynamic parameters in a video file containing ocular microvascular information o
Step 2: for image set I o The image registration is carried out on the continuous images, the first image is taken as a fixed image, the rest images are aligned with the first image, and a registered image set I is obtained r
And 3, step 3: carrying out image segmentation on the registered image to obtain a segmentation mask image set S of the ocular microvasculature r
And 4, step 4: obtaining a segmented mask image set S r The central line of the middle eye microvasculature;
and 5: fitting the central line, extracting the edge of the eye microvessels in the segmentation mask image, and calculating the average observation diameter, the inner diameter and the outer diameter of the eye microvessels according to the fitted central line and the edge of the eye microvessels in the segmentation mask image;
step 6: extracting the first fraction in step 5Cutting each fitted central line in the mask image, taking the coordinates of the central line as the column of the space-time diagram, and taking the image set I r The number of the images in (1) is line, and a space-time image is formed;
and 7: carrying out histogram equalization processing on the obtained spatio-temporal image;
and 8: carrying out linear detection processing on the spatio-temporal image, carrying out statistics on the slope of the obtained straight line, and obtaining the slope k with the most dense distribution;
and step 9: calculating blood flow velocity of ocular microvasculature
Figure DEST_PATH_IMAGE001
The following formula is adopted:
Figure 932614DEST_PATH_IMAGE002
preferably, in step 3, the contour in the segmentation result needs to be calculated, and the region with small contour is removed.
Preferably, in step 5, the method specifically includes the following steps:
step 51: sampling at equal intervals from the starting point of the central line of each blood vessel;
step 52: sequentially calculating the coordinate Pi and the normal vector of each sampling point on the blood vessel center line, and obtaining the observation diameter of the blood vessel at each sampling point according to the edge of the eye microvessel of the segmentation mask image along the normal vector direction;
step 53: counting the diameters of all the sampling points, and obtaining a normal distribution parameter by adopting a normal distribution approximation mode, wherein the highest probability is the inner diameter B of the ocular capillary i The diameter corresponding to the confidence coefficient of 0.94 is the outer diameter B of the ocular microvasculature o The mean of the observed diameters is the mean observed diameter B m
The invention also provides an analysis method of blood flow volume of ocular microvasculature, which adopts the following formula:
Figure 69328DEST_PATH_IMAGE003
q is blood flow, B m Is the average observed diameter obtained in step 53 above.
The invention also provides an analysis method of the tube wall shear rate of the eyeball blood vessel, which comprises the following steps: the following formula is adopted:
Figure 133886DEST_PATH_IMAGE004
and WSR is the pipe wall shear rate.
The invention also provides an analysis method of the diameter of the ocular microvasculature, which specifically comprises the following steps:
step 1: inputting an image set I which needs to carry out analysis on ocular microvascular hemodynamic parameters in a video file containing ocular microvascular information o
Step 2: for image set I o Performing image registration on the continuous images, taking the first image as a fixed image, aligning the residual images with the first image to obtain a registered image set I r
And step 3: carrying out image segmentation on the registered image to obtain a segmentation mask image set S of the ocular microvasculature r
And 4, step 4: obtaining a segmented mask image set S r The centerline of the middle vessel;
and 5: fitting the central line, extracting the edge of the ocular microvasculature in the segmentation mask image, and calculating the average diameter, the inner diameter and the outer diameter of the ocular microvasculature according to the fitted central line and the edge of the ocular microvasculature in the segmentation mask image.
Preferably, in step 3, the contour in the segmentation result needs to be calculated, and the region with small contour is removed.
Preferably, in step 5, the method specifically includes the following steps:
step 51: sampling at equal intervals from the starting point of the central line of each blood vessel;
step 52: sequentially calculating the coordinates Pi and normal vector of each sampling point on the blood vessel central line, and obtaining the observation diameter of the blood vessel at each sampling point according to the edge of the ocular microvasculature of the segmentation mask image along the normal vector direction;
step 53: counting the diameters of all the sampling points, and obtaining a normal distribution parameter by adopting a normal distribution approximation mode, wherein the highest probability is the inner diameter B of the ocular capillary i The diameter corresponding to the confidence coefficient of 0.94 is the outer diameter B of the ocular microvasculature o The mean of the observed diameters is the mean observed diameter B m
Compared with the prior art, the invention has the following beneficial effects:
1) The processing speed is high, unnecessary human-computer interaction is reduced, and automatic analysis of the eye microvascular hemodynamics can be completed quickly; 2) The calculated ocular microvascular kinetic parameters have high precision, curve fitting is carried out on the center line of the blood vessel, and the image of blood vessel wall red blood cells is reduced by adopting a Gaussian distribution method; 3) A statistical method is introduced into the space-time image, so that the accuracy of blood flow parameters is improved.
Drawings
FIG. 1 is a flow chart in accordance with a preferred embodiment of the present invention;
FIG. 2 (a) is an original spatiotemporal image;
FIG. 2 (b) is a spatiotemporal image that has been histogram equalized;
fig. 2 (c) is a spatiotemporal image subjected to the line detection process.
Detailed Description
The following description is provided to disclose the invention so as to enable any person skilled in the art to practice the invention. The preferred embodiments in the following description are given by way of example only, and other obvious variations will occur to those skilled in the art.
An automatic analysis method for hemodynamics parameters of ocular microvasculature, wherein the dynamics parameters mainly comprise blood flow velocity, blood flow volume and tube wall shear rate, and the method specifically comprises the following steps:
step 1: image set I which is input into video file and needs to carry out analysis on ocular microvascular hemodynamic parameters o And input the relevant parameters.
The video file is obtained by shooting the eye microvasculature through a camera and the likeAnd (4) obtaining. For the image set I o The video file is composed of a plurality of successive pictures, a picture set I o Is the collection of these images.
Step 2: for image set I o Performing image registration on the continuous images, taking the first image as a fixed image, aligning the residual images with the first image to obtain a registered image set I r . The registration algorithm may employ a mutual information based rigid registration or other configuration methods in the prior art.
And step 3: carrying out image segmentation on the registered image by adopting a U-Net network, calculating the outline in the segmentation result, removing the area with small outline, and obtaining a segmentation mask image set S of the blood vessel r
And 4, step 4: obtaining a segmentation image set S by adopting a skeleton extraction method of an image r The center line of the middle vessel.
And 5: and (3) spline fitting is carried out on the central line, the edge of the blood vessel in the segmentation image is extracted, the diameter of the blood vessel is calculated, and the average observed diameter, the inner diameter and the outer diameter of the blood vessel are obtained by adopting a Gaussian fitting method.
Specifically, the central line is processed by using a fixed mask, the starting point and the end point of each blood vessel central line are detected and recorded, equal-interval sampling is performed from the starting point of each blood vessel central line, the interval size is the number of pixels, the number of pixels between two adjacent sampling points is the same, the specific number of pixels can be adjusted according to the resolution of the image, and the interval can be larger as the resolution is higher until the end point is reached. B-Spline is adopted to carry out blood vessel center line fitting, after fitting, coordinates Pi and normal vector of each sampling point on the blood vessel center line are calculated in sequence, the diameter of the blood vessel at the sampling point is obtained according to the boundary of a blood vessel mask image along the direction of the normal vector, the observation diameters of all the sampling points in the blood vessel are obtained in the same mode, the observation diameters refer to two times of the value from the sampling point to the boundary of the blood vessel along the normal vector, but when the boundary of the blood vessel image is actually obtained, the diameters of the blood vessel corresponding to the images at different positions are different due to the problem of the images, so that the blood vessel sampling points are subjected to blood vessel samplingAll diameters are subjected to statistical analysis, a normal distribution parameter is obtained by adopting a normal distribution approximation mode, and the highest probability is the inner diameter B of the blood vessel i The position with the confidence coefficient of 0.94 is the outer diameter B of the blood vessel o The mean of all observed diameters is the mean observed diameter B m Wall thickness of tube = B o -B m
Step 6: extracting each fitted central line in the first blood vessel mask image in the step 5, taking the coordinates of the central line as the column of the space-time image, and collecting the images I r Wherein all the image gray values on the central line are extracted to form a group, and an image set I r The number of images in (a) is a row value, that is, a spatiotemporal image is obtained according to the blood flow situation at the central line at different times as shown in fig. 2 (a). In other words, the video is composed of a set of images I r In the image composition of (1), the number of images in the videos with different time lengths is different, the time of the videos corresponding to the number of different images is different from 0 for the same video, and if the number of images in the video of 1 minute is 200, the number of images in the video of two minutes is 400, conversely, 200 images correspond to the video within 0-1 minute, and 400 images correspond to the video within 0-2 minutes, so that the number of images corresponds to the time in the video, and the number of images is taken as a line, which is actually the time in the video.
And 7: the obtained spatio-temporal image is subjected to histogram equalization processing, as shown in fig. 2 (b). Of course, other methods known in the art may be used for equalization.
And 8: and (3) performing linear detection processing on the spatio-temporal image processed in the step (7), counting the obtained linear, counting the distribution of the slope, and taking the slope k with the most dense distribution as the basis of blood flow calculation, as shown in fig. 2 (c). When the straight line is counted, it is necessary to remove the case where the slope of the straight line is 0 and the vertical abscissa.
And step 9: calculating the blood flow velocity of the blood vessel by adopting the following formula:
Figure 217510DEST_PATH_IMAGE005
,
wherein S represents the distance between image pixels, in um, and is input after the resolution is determined by the camera,
Figure 890937DEST_PATH_IMAGE006
is the cross-sectional flow velocity of the blood vessel;
Figure 336569DEST_PATH_IMAGE007
wherein the content of the first and second substances,
Figure 524973DEST_PATH_IMAGE008
the axial flow velocity of the blood vessel is the blood flow velocity of the blood vessel.
Step 10: calculating the flow of the blood vessel by adopting the following formula:
flow rate (Q) (pl/s):
Figure 520611DEST_PATH_IMAGE009
,B m is the average observed diameter obtained in step 5.
Step 11, calculating the vessel wall shear rate of the blood vessel by adopting the following formula:
wall Shear Rate (WSR) (s-1):
Figure 753272DEST_PATH_IMAGE010
the foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are merely illustrative of the principles of the invention, but that various changes and modifications may be made without departing from the spirit and scope of the invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.

Claims (5)

1. An analysis method of blood flow velocity of ocular microvasculature, comprising the steps of:
step 1: inputting an image set which needs to analyze the blood flow dynamic parameters of the ocular microvasculature in a video file containing the information of the ocular microvasculature
Figure 651662DEST_PATH_IMAGE001
Step 2: for image set
Figure 785840DEST_PATH_IMAGE002
Performing image registration on the continuous images, taking the first image as a fixed image, aligning the residual images with the first image to obtain a registered image set
Figure 85103DEST_PATH_IMAGE003
And step 3: carrying out image segmentation on the registered image to obtain a segmentation mask image set of the ocular microvasculature
Figure 797101DEST_PATH_IMAGE004
And 4, step 4: obtaining a segmented mask image set
Figure 853918DEST_PATH_IMAGE004
The central line of the middle-eye microvasculature;
and 5: fitting the central line, extracting the edge of the eye microvessels in the segmentation mask image, and calculating the average observation diameter, the inner diameter and the outer diameter of the eye microvessels according to the fitted central line and the edge of the eye microvessels in the segmentation mask image;
step 6: extracting each fitted central line in the first segmentation mask image in the step 5, taking the coordinates of the central lines as the columns of the space-time diagram, and taking the image set
Figure 276941DEST_PATH_IMAGE003
The number of the images in (1) is line, and a space-time image is formed;
and 7: carrying out histogram equalization processing on the obtained space-time image;
and 8: carrying out linear detection processing on the spatio-temporal image, carrying out statistics on the slope of the obtained straight line, and obtaining the slope k with the most dense distribution;
and step 9: calculating blood flow velocity of ocular microvasculature
Figure 617661DEST_PATH_IMAGE005
The following formula is adopted:
Figure 733516DEST_PATH_IMAGE006
wherein S represents the distance between image pixels, the unit is um, and the resolution is determined by the camera and then input.
2. The method of analyzing an ocular microvascular blood flow velocity according to claim 1, wherein in step 3, it is further necessary to calculate a contour in the segmentation result and remove a region with a small contour.
3. The method for analyzing an ocular microvascular blood flow velocity according to claim 1, further comprising, in step 5, the steps of:
step 51: sampling at equal intervals from the starting point of the central line of each blood vessel;
step 52: sequentially calculating the coordinate Pi and the normal vector of each sampling point on the blood vessel center line, and obtaining the observation diameter of the blood vessel at each sampling point according to the edge of the eye microvessel of the segmentation mask image along the normal vector direction;
step 53: counting the diameters of all the sampling points, and obtaining a normal distribution parameter by adopting a normal distribution approximation mode, wherein the highest probability is the inner diameter of the ocular capillary
Figure 338283DEST_PATH_IMAGE007
The diameter corresponding to the confidence coefficient of 0.94 is the ocular microvasculatureOuter diameter of
Figure 627313DEST_PATH_IMAGE008
The mean of the observed diameters is the average observed diameter
Figure 619278DEST_PATH_IMAGE009
4. An analysis method of blood flow volume of ocular microvasculature, characterized by comprising the following steps:
step 1: inputting image set required to analyze blood flow dynamic parameters of ocular microvasculature in video file containing ocular microvasculature information
Figure 437192DEST_PATH_IMAGE001
Step 2: for image set
Figure 986378DEST_PATH_IMAGE002
Performing image registration on the continuous images, taking the first image as a fixed image, aligning the rest images with the first image to obtain a registered image set
Figure 485624DEST_PATH_IMAGE003
And step 3: carrying out image segmentation on the registered image to obtain a segmentation mask image set of the ocular microvasculature
Figure 286090DEST_PATH_IMAGE004
And 4, step 4: obtaining a segmented mask image set
Figure 711124DEST_PATH_IMAGE004
The central line of the middle eye microvasculature;
and 5: fitting the central line, extracting the edge of the eye microvessels in the segmentation mask image, and calculating the average observation diameter, the inner diameter and the outer diameter of the eye microvessels according to the fitted central line and the edge of the eye microvessels in the segmentation mask image;
step 6: extracting each fitted central line in the first segmentation mask image in the step 5, taking the coordinates of the central lines as columns of a space-time diagram, and taking the image set
Figure 980562DEST_PATH_IMAGE003
The number of the images in (1) is line, and a space-time image is formed;
and 7: carrying out histogram equalization processing on the obtained space-time image;
and 8: carrying out linear detection processing on the spatio-temporal image, carrying out statistics on the slope of the obtained straight line, and obtaining the slope k with the most dense distribution;
and step 9: calculating the blood flow volume of ocular microvasculature by using the following formula:
Figure 798346DEST_PATH_IMAGE010
q is the blood flow, S represents the distance between image pixels, in um, input after the resolution is determined by the camera,
Figure 236674DEST_PATH_IMAGE009
the average observation diameter is obtained by adopting the following steps:
step 51: sampling at equal intervals from the starting point of the central line of each blood vessel;
step 52: sequentially calculating the coordinates Pi and normal vector of each sampling point on the blood vessel central line, and obtaining the observation diameter of the blood vessel at each sampling point according to the edge of the ocular microvasculature of the segmentation mask image along the normal vector direction;
step 53: counting the diameters of all the sampling points, and obtaining a normal distribution parameter by adopting a normal distribution approximation mode, wherein the highest probability is the inner diameter of the ocular capillary
Figure 130811DEST_PATH_IMAGE007
The diameter corresponding to the confidence coefficient of 0.94 is the outer diameter of the ocular microvasculature
Figure 120502DEST_PATH_IMAGE008
The mean of the observed diameters is the average observed diameter
Figure 210818DEST_PATH_IMAGE009
5. The method for analyzing the tube wall shear rate of the eyeball blood vessel is characterized by comprising the following steps:
step 1: inputting image set required to analyze blood flow dynamic parameters of ocular microvasculature in video file containing ocular microvasculature information
Figure 736608DEST_PATH_IMAGE001
Step 2: for image set
Figure 771953DEST_PATH_IMAGE002
Performing image registration on the continuous images, taking the first image as a fixed image, aligning the residual images with the first image to obtain a registered image set
Figure 265252DEST_PATH_IMAGE003
And step 3: carrying out image segmentation on the registered image to obtain a segmentation mask image set of ocular microvasculature
Figure 441149DEST_PATH_IMAGE004
And 4, step 4: obtaining a segmented mask image set
Figure 851140DEST_PATH_IMAGE004
The central line of the middle-eye microvasculature;
and 5: fitting the central line, extracting the edge of the eye microvessels in the segmentation mask image, and calculating the average observation diameter, the inner diameter and the outer diameter of the eye microvessels according to the fitted central line and the edge of the eye microvessels in the segmentation mask image;
and 6: extracting each fitted central line in the first segmentation mask image in the step 5, taking the coordinates of the central lines as the columns of the space-time diagram, and taking the image set
Figure 389876DEST_PATH_IMAGE003
The number of the images in (1) is line, and a space-time image is formed;
and 7: carrying out histogram equalization processing on the obtained space-time image;
and 8: carrying out linear detection processing on the spatio-temporal image, carrying out statistics on the slope of the obtained straight line, and obtaining the slope k with the most dense distribution;
and step 9: calculating the tube wall shear rate of the eyeball blood vessel by adopting the following formula:
Figure 652361DEST_PATH_IMAGE011
WSR is the pipe wall shear rate, S represents the distance between image pixels, the unit is um, the resolution is determined by a camera and then input,
Figure 130484DEST_PATH_IMAGE009
the average observation diameter is obtained by adopting the following steps:
step 51: sampling at equal intervals from the starting point of the central line of each blood vessel;
step 52: sequentially calculating the coordinate Pi and the normal vector of each sampling point on the blood vessel center line, and obtaining the observation diameter of the blood vessel at each sampling point according to the edge of the eye microvessel of the segmentation mask image along the normal vector direction;
step 53: counting the diameters of all the sampling points, obtaining a normal distribution parameter by adopting a normal distribution approximation mode, wherein the position with the highest probability is the inner diameter of the ocular microvasculature
Figure 162025DEST_PATH_IMAGE007
The diameter corresponding to the confidence coefficient of 0.94 is the outer diameter of the ocular microvasculature
Figure 335911DEST_PATH_IMAGE008
The mean of the observed diameters is the mean observed diameter
Figure 272643DEST_PATH_IMAGE009
CN202011017868.4A 2020-09-24 2020-09-24 Automatic analysis method for ocular microvascular hemodynamic parameters Active CN112164046B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011017868.4A CN112164046B (en) 2020-09-24 2020-09-24 Automatic analysis method for ocular microvascular hemodynamic parameters

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011017868.4A CN112164046B (en) 2020-09-24 2020-09-24 Automatic analysis method for ocular microvascular hemodynamic parameters

Publications (2)

Publication Number Publication Date
CN112164046A CN112164046A (en) 2021-01-01
CN112164046B true CN112164046B (en) 2022-12-16

Family

ID=73863722

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011017868.4A Active CN112164046B (en) 2020-09-24 2020-09-24 Automatic analysis method for ocular microvascular hemodynamic parameters

Country Status (1)

Country Link
CN (1) CN112164046B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113012198B (en) * 2021-03-22 2022-04-01 复旦大学 Blood flow automatic quantitative analysis method in intracranial and extracranial vascular anastomosis

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110060275A (en) * 2019-04-11 2019-07-26 霍尔果斯奇妙软件科技有限公司 A kind of method and system of microcirculation in human body blood flow velocity detection

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101474083A (en) * 2009-01-15 2009-07-08 西安交通大学 System and method for super-resolution imaging and multi-parameter detection of vascular mechanical characteristic
EP3179902B1 (en) * 2014-08-11 2020-10-14 The Board of Trustees of the University of Illionis Epidermal device for analysis of temperature and thermal transport characteristics
US10176408B2 (en) * 2015-08-14 2019-01-08 Elucid Bioimaging Inc. Systems and methods for analyzing pathologies utilizing quantitative imaging
US11094058B2 (en) * 2015-08-14 2021-08-17 Elucid Bioimaging Inc. Systems and method for computer-aided phenotyping (CAP) using radiologic images
AU2016379175A1 (en) * 2015-12-21 2018-07-05 The Regents Of The University Of California Perfusion digital subtraction angiography
CN109069114B (en) * 2017-02-16 2021-10-22 深圳迈瑞生物医疗电子股份有限公司 Ultrasonic medical detection equipment, imaging method, imaging system and display terminal
CN110313943A (en) * 2018-03-30 2019-10-11 佳能医疗系统株式会社 Medical diagnostic apparatus, medical image-processing apparatus and image processing method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110060275A (en) * 2019-04-11 2019-07-26 霍尔果斯奇妙软件科技有限公司 A kind of method and system of microcirculation in human body blood flow velocity detection

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Co-Occurrence and Cyclical Growth Law Analysis of User Innovation Knowledge Map Based on Temporal-Weighted Network;Qiaojiayu Wang et al.;《IEEE Access ( Volume: 7)》;20190501;第60026 - 60041页 *
吕菲 ; 赵兴群.一种基于ST图轨迹方向的微循环流速测量算法.《生物医学工程学杂志》.2015, *
基于图像分析的人体微循环参数测量及识别研究;陈远;《中国博士学位论文全文数据库 信息科技辑》;20140615;第I138-22页 *
基于视频图像的微循环参数测量算法研究;王钰;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180415;第I138-2794页 *
基于视频图像的微循环参数自动化测量研究;王钰 等;《全国第17次光纤通信暨第18届集成光学学术会议——可实施光元部件、光纤光缆、子系统、系统、网络与应用最新进展专辑》;20151218;第1-5页 *
基于视频图像的微循环血流分析系统的研制;吕菲;《中国优秀硕士学位论文全文数据库 信息科技辑》;20160815;第I138-894页 *

Also Published As

Publication number Publication date
CN112164046A (en) 2021-01-01

Similar Documents

Publication Publication Date Title
WO2021169128A1 (en) Method and apparatus for recognizing and quantifying fundus retina vessel, and device and storage medium
CN109544512B (en) Multi-mode-based embryo pregnancy result prediction device
CN108765392B (en) Digestive tract endoscope lesion detection and identification method based on sliding window
CN111798425B (en) Intelligent detection method for mitotic image in gastrointestinal stromal tumor based on deep learning
WO2020019757A1 (en) Intelligent intracranial plaque and carotid plaque segmentation method for magnetic resonance image
CN107292835B (en) Method and device for automatically vectorizing retinal blood vessels of fundus image
CN111667456A (en) Method and device for detecting vascular stenosis in coronary artery X-ray sequence radiography
CN112164046B (en) Automatic analysis method for ocular microvascular hemodynamic parameters
CN109754388B (en) Carotid artery stenosis degree calculation method and device and storage medium
CN108670297A (en) Parkinson's disease assistant diagnosis system and method based on multi-modal TCD,transcranial Doppler
CN110610480B (en) MCASPP neural network eyeground image optic cup optic disc segmentation model based on Attention mechanism
Jian et al. Dual-Branch-UNet: A Dual-Branch Convolutional Neural Network for Medical Image Segmentation.
Caliskan et al. Automated bolus detection in videofluoroscopic images of swallowing using mask-RCNN
Patwari et al. Classification and calculation of retinal blood vessels parameters
CN112700409A (en) Automatic retinal microaneurysm detection method and imaging method
CN112419246B (en) Depth detection network for quantifying esophageal mucosa IPCLs blood vessel morphological distribution
CN110033861B (en) Blood vessel and macular avascular region quantitative analysis method and system suitable for OCTA image
CN115797329A (en) Fundus image optic disk optic cup segmentation method based on improved Unet network
CN111292285A (en) Automatic screening method for diabetes mellitus based on naive Bayes and support vector machine
CN111798426B (en) Deep learning and detecting system for mitotic image in gastrointestinal stromal tumor of moving end
Huang et al. Automatic quantitative analysis of bioresorbable vascular scaffold struts in optical coherence tomography images using region growing
CN114266742A (en) Method for detecting narrow region of cerebral blood vessel CTA image
CN114494177A (en) IVOCT (in-vivo visual optical coherence tomography) branch blood vessel identification method by utilizing longitudinal section and withdrawal property
Sathishkumar Segmentation of optic nerve head for glaucoma detection using fundus images
TW201309264A (en) Automatic nailfold microcirculation analysis system

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
TA01 Transfer of patent application right

Effective date of registration: 20210819

Address after: 510000 No. 7, Jinsui Road, Tianhe District, Guangzhou City, Guangdong Province

Applicant after: ZHONGSHAN OPHTHALMIC CENTER, SUN YAT-SEN University

Applicant after: SUN YAT-SEN University

Address before: Zhongshan Ophthalmic Center, Sun Yat sen University, No.7, Jinsui Road, Tianhe District, Guangzhou, Guangdong 510000

Applicant before: ZHONGSHAN OPHTHALMIC CENTER, SUN YAT-SEN University

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20230426

Address after: Room 406, No. 1, Yichuang Street, Huangpu District (Zhongxin Guangzhou Knowledge City), Guangzhou City, Guangdong Province, 510000

Patentee after: Meishi Optical Technology (Guangdong) Co.,Ltd.

Address before: 510000 No. 7, Jinsui Road, Tianhe District, Guangzhou City, Guangdong Province

Patentee before: ZHONGSHAN OPHTHALMIC CENTER, SUN YAT-SEN University

Patentee before: SUN YAT-SEN University

TR01 Transfer of patent right