CN113951813A - Retinal blood vessel branch angle calculation method and device and electronic equipment - Google Patents

Retinal blood vessel branch angle calculation method and device and electronic equipment Download PDF

Info

Publication number
CN113951813A
CN113951813A CN202111320203.5A CN202111320203A CN113951813A CN 113951813 A CN113951813 A CN 113951813A CN 202111320203 A CN202111320203 A CN 202111320203A CN 113951813 A CN113951813 A CN 113951813A
Authority
CN
China
Prior art keywords
black
white image
optic disc
retinal
determining
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.)
Pending
Application number
CN202111320203.5A
Other languages
Chinese (zh)
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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202111320203.5A priority Critical patent/CN113951813A/en
Publication of CN113951813A publication Critical patent/CN113951813A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/103—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for determining refraction, e.g. refractometers, skiascopes
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/14—Arrangements specially adapted for eye photography
    • 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
    • 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/30041—Eye; Retina; Ophthalmic
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Ophthalmology & Optometry (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Eye Examination Apparatus (AREA)

Abstract

The invention provides a method and a device for calculating a branch angle of a retinal blood vessel and electronic equipment, wherein the method for calculating the branch angle of the retinal blood vessel comprises the following steps: obtaining a retina fundus image, and performing binarization processing on the retina fundus image to obtain a black-and-white image; determining a central point of an optic disc based on the black-and-white image, and determining two target branch vessels of the retina based on the black-and-white image; and obtaining the branch angle of the retinal vessel based on the central point of the optic disc and the two target branch vessels. The method, the device and the electronic equipment for calculating the branch angle of the retinal vessel can accurately calculate the branch angle of the retinal vessel of the human eye.

Description

Retinal blood vessel branch angle calculation method and device and electronic equipment
Technical Field
The invention relates to the technical field of image processing, in particular to a method and a device for calculating a branch angle of a retinal blood vessel and electronic equipment.
Background
An eyeground camera belongs to the field of medical imaging and is used for acquiring retinal eyeground images so as to process and observe retinas of human eyes in the following process. Since the blood vessels of the fundus are the only blood vessels that can be directly observed by the human body through the body surface, the state of the optic nerve, retina, choroid and dioptric medium of the fundus can be checked, i.e. the state of the human eye can be determined, by the retinal fundus image acquired by the fundus camera.
When the state of human eyes is determined, the included angle between two main branches in the middle of a retinal blood vessel, namely the branch angle of the retinal blood vessel, generally needs to be calculated, but a scheme capable of accurately calculating the branch angle of the retinal blood vessel of the human eyes is still lacked at present.
Disclosure of Invention
The invention provides a method and a device for calculating a branch angle of a retinal blood vessel and electronic equipment, which can accurately calculate the branch angle of the retinal blood vessel of a human eye.
The invention provides a method for calculating a branch angle of a retinal blood vessel, which comprises the following steps:
obtaining a retina fundus image, and performing binarization processing on the retina fundus image to obtain a black-and-white image;
determining a central point of an optic disc based on the black-and-white image, and determining two target branch vessels of the retina based on the black-and-white image;
and obtaining the branch angle of the retinal vessel based on the central point of the optic disc and the two target branch vessels.
According to the method for calculating the branch angle of the retinal vessel provided by the invention, the determination of the central point of the optic disc based on the black-and-white image comprises the following steps:
determining a video disc based on the black and white image;
and determining the position of a target pixel point in the optic disc, and determining the central point of the optic disc based on the position of the target pixel point in the optic disc.
According to the method for calculating the branch angle of the retinal vessel provided by the invention, the determination of the optic disk based on the black-and-white image comprises the following steps:
traversing the black-and-white image based on a preset first traversal window; wherein the first traversal window is a circular window;
and taking the area with the maximum concentration of the target pixel points in the first traversal window as the video disc.
According to the method for calculating the retinal blood vessel branch angle provided by the invention, the determination of the central point of the optic disc based on the position of the target pixel point in the optic disc comprises the following steps:
creating a two-dimensional coordinate system on the black-and-white image;
under the two-dimensional coordinate system, the position of the central point of the optic disc is calculated and obtained based on the following formula:
Figure BDA0003345323850000021
wherein x isp、ypRespectively an abscissa and an ordinate of the central point of the optic disc,
Figure BDA0003345323850000022
respectively the mean value of the abscissa and the mean value of the ordinate, x, of all target pixel points in the black-and-white imagei、yiRespectively is the abscissa and the ordinate of the ith target pixel point in the black-and-white image.
According to the method for calculating the branch angle of the retinal vessel provided by the invention, two target branch vessels of the retina are determined based on the black-and-white image:
based on the central point of the optic disc, a polar coordinate system is established;
under the polar coordinate system, a preset second traversal window is rotated to traverse the black-and-white image, and two target branch vessels are determined;
the second traversal window is a trilateral, and the vertex of the trilateral is overlapped with the central point of the optic disc.
According to the method for calculating the branch angle of the retinal vessel provided by the invention, in the polar coordinate system, the black-and-white image is traversed by rotating a preset second traversal window to determine two target branch vessels, which comprises the following steps:
and under the polar coordinate system, traversing the black-and-white image by rotating a preset second traversal window, and taking two areas with the largest target pixel points in the second traversal window as two target branch vessels respectively.
The present invention also provides a retinal blood vessel branch angle calculation apparatus, including:
the image processing module is used for acquiring a retina fundus image and carrying out binarization processing on the retina fundus image to obtain a black-and-white image;
the characteristic determining module is used for determining the central point of the optic disc based on the black-and-white image and determining two target branch blood vessels of the retina based on the black-and-white image;
and the angle calculation module is used for obtaining the retinal blood vessel branch angle based on the central point of the optic disc and the two target branch blood vessels.
The invention also provides an electronic device, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the program to realize the steps of the retinal blood vessel branch angle calculation method.
The present invention also provides a non-transitory computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the retinal vessel branch angle calculation method as described in any one of the above.
The present invention also provides a computer program product comprising a computer program which, when executed by a processor, implements the steps of the method for calculating a retinal vessel branch angle as described in any one of the above.
According to the method, the device and the electronic equipment for calculating the retinal blood vessel branch angle, the corresponding black and white image is obtained by performing binarization processing on the retinal fundus image; on the basis of the black-white image, the target branch blood vessel corresponding to the central point of the optic disc and the retinal blood vessel can be determined, and then the branch angle of the retinal blood vessel can be accurately calculated.
Drawings
In order to more clearly illustrate the technical solutions of the present invention or the prior art, the drawings needed for the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
FIG. 1 is a schematic flow chart of a retinal vessel branch angle calculation method provided by the present invention;
FIG. 2 is a black-and-white image corresponding to the retinal vessel branch angle calculation method provided by the present invention;
FIG. 3 is a schematic view of a video disc in a black and white image provided by the present invention;
FIG. 4 is a schematic diagram of a black-and-white image divided into four regions according to the present invention;
FIG. 5 is a schematic structural diagram of a retinal vessel branch angle calculation apparatus provided by the present invention;
fig. 6 is a schematic structural diagram of an electronic device provided in the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The following describes a retinal vessel branch angle calculation method, device and electronic device according to the present invention with reference to fig. 1 to 6.
As shown in fig. 1, the present invention provides a retinal blood vessel branch angle calculation method, including:
and 110, acquiring a retina fundus image, and performing binarization processing on the retina fundus image to obtain a black-and-white image.
It can be understood that the process of performing binarization processing on the retinal fundus image, that is, setting the gray value of the pixel point on the retinal fundus image to be 0 or 255, and rendering the whole retinal fundus image to be obviously black or white, is performed.
And step 120, determining the central point of the optic disc based on the black-and-white image, and determining two target branch vessels of the retina based on the black-and-white image.
It can be understood that the statistical analysis can be performed on the pixel points in the black-and-white image, so as to determine the corresponding central point of the optic disc and two target branch vessels of the retina.
And step 130, obtaining a retinal vessel branch angle based on the central point of the optic disc and the two target branch vessels.
It is understood that the retinal vessel branch angle is the angle between two target branch vessels corresponding to the retinal vessel from the central point of the optic disc.
In some embodiments, said determining, based on said black and white image, a disk center point comprises:
determining a video disc based on the black and white image;
and determining the position of a target pixel point in the optic disc, and determining the central point of the optic disc based on the position of the target pixel point in the optic disc.
It is understood that since the eyeball is a nearly spherical body, the fundus camera is not vertically photographing either, and therefore the disk center position is not simply equivalent to the disk geometric center position. Therefore, it is necessary to read all target pixel positions in the video disc, perform statistical analysis on all target pixel positions, and determine the central point of the video disc.
It should be noted that the target pixel is also a pixel corresponding to the blood vessel. Further, in the black and white image, the target pixel may be a white pixel.
A black-and-white image corresponding to the retinal fundus image is shown in fig. 2, a optic disc in the black-and-white image is shown in fig. 3, a black circle region shown in fig. 3 is the optic disc, and a central point of the black circle region is a central point of the optic disc.
In some embodiments, said determining a disc based on said black and white image comprises:
traversing the black-and-white image based on a preset first traversal window; wherein the first traversal window is a circular window;
and taking the area with the maximum concentration of the target pixel points in the first traversal window as the video disc.
It is understood that the first traversal window may be a circular window with a preset size, and the relative center position of the optic disc, that is, the center position of the optic disc, may be determined by using a clustering method by reading all target pixel point positions in the first traversal window.
And taking the area with the maximum concentration of the target pixel points in the first traversal window as the video disc. Specifically, in the process of traversing the black-and-white image by using the first traversal window, the concentrations of the target pixels in the first traversal window are different, and when the target pixels in the first traversal window reach a certain moment, the concentration of the target pixels in the first traversal window is the moment with the maximum concentration of the target pixels in the traversal process, and then the area where the first traversal window is located at the moment is determined as the optic disc.
As shown in fig. 3, traversing the black-and-white image from the upper left corner to the lower right corner of the black-and-white image, the region of the black-and-white image having the maximum white pixel density in the black background within a specific window is determined as the root of the vessel tree or optic disc.
The location of the optic disc is determined using the attributes of the circle to color the vessel tree or the root of the optic disc in the black and white image. Setting the center of a circle as (x)1，y1) The radius r, stains the root of the vessel tree or optic disc located inside the circle, and the place where the stained points are most dense can be determined as the optic disc area.
In some embodiments, said determining a disc center point based on a target pixel point location in said disc comprises:
creating a two-dimensional coordinate system on the black-and-white image;
under the two-dimensional coordinate system, the position of the central point of the optic disc is calculated and obtained based on the following formula:
Figure BDA0003345323850000071
wherein x isp、ypRespectively an abscissa and an ordinate of the central point of the optic disc under a two-dimensional coordinate system,
Figure BDA0003345323850000072
respectively is the average value of the abscissa and the average value of the ordinate, x, of all target pixel points in the black-and-white image under a two-dimensional coordinate systemi、yiRespectively is the abscissa and the ordinate of the ith target pixel point in the black-and-white image under a two-dimensional coordinate system.
It can be understood that the above calculation formula adopts a clustering manner to obtain the position of the central point of the optic disc. When the position of the central point of the optic disc is calculated based on the formula, a two-dimensional coordinate system needs to be constructed first, the origin of coordinates of the two-dimensional coordinate system can be the central point of the optic disc, and the coordinate system can be established by taking the horizontal direction of the black-and-white image as the X axis and the vertical direction of the black-and-white image as the Y axis.
In some embodiments, said determining two target branch vessels of the retina based on said black and white image comprises:
based on the central point of the optic disc, a polar coordinate system is established;
and traversing the black-and-white image by rotating a preset second traversal window under the polar coordinate system to determine two target branch vessels.
The second traversal window is a trilateral, and the vertex of the trilateral is overlapped with the central point of the optic disc.
It will be appreciated that the optic disc center point may be the pole of the polar coordinate system. The second traversal window is a triangle, i.e. a shape composed of three sides, and the second traversal window may be an isosceles triangle window or a sector window.
Wherein, the angle between the first side and the second side of the second traversal window is 2-5 degrees, for example, it may be 3 degrees, and the intersection point of the first side and the second side coincides with the pole.
The length of the second traversal window is one third to two thirds of the width of the grayscale picture, for example, the length of the second traversal window can be one half, so as to avoid selecting a blood vessel branch which is too far away.
In some embodiments, the traversing the black-and-white image by rotating a preset second traversal window in the polar coordinate system to determine two target branch vessels includes:
and under the polar coordinate system, traversing the black-and-white image by rotating a preset second traversal window, and taking two areas with the largest target pixel points in the second traversal window as two target branch vessels respectively.
It will be appreciated that we link the relative center of the optic disc (i.e., the disc center) to all off-disc white pixels (i.e., the blood vessels) and calculate the slope of the line between the two points (i.e., the vector direction) at the same time.
In the polar coordinate system, traversal is performed from 0 degree based on the second traversal window, and in the polar coordinate system, traversal is performed counterclockwise with the vertical direction of the black-and-white image as the 0 degree direction. During traversal, the position with the most white pixels contained in the second traversal window is used as the position of the main retinal blood vessel, namely the position of the target branch blood vessel of the retina.
The angle of the straight line connecting the center of the target branch blood vessel and the central point of the optic disc is taken as the slope of the blood vessel. The calculation formula of the slope of the straight line is as follows:
Figure BDA0003345323850000081
wherein x isp、ypRespectively, the abscissa and the ordinate of the central point position of the visual disc under the two-dimensional coordinate system, and x and y are the abscissa and the ordinate of the corresponding center of the target branch blood vessel under the two-dimensional coordinate system.
Further, as shown in fig. 4, the black-and-white image may be divided into four regions along the polar coordinates, and the position of the main retinal blood vessel and the slope thereof are found in region 2 and region 4, respectively, for example, in region 2, the scanning is started from 0 degree and the slope is calculated from 1.5 degrees. In each area of the black-and-white image, the vertical direction is the 0 reading direction.
The angle tan θ between the two target branch vessels of the retina corresponds to the following calculation formula:
Figure BDA0003345323850000082
wherein if L1-L2|≥180°Error messages may be returned. L is1The slope L of a straight line connecting the center of the first target branch blood vessel and the center of the optic disk in a two-dimensional coordinate system2The slope of a straight line connecting the center of the second target branch blood vessel and the central point of the optic disc in a two-dimensional coordinate system, ". is the point multiplication operation.
In summary, the method for calculating the branch angle of the retinal blood vessel provided by the present invention includes: obtaining a retina fundus image, and performing binarization processing on the retina fundus image to obtain a black-and-white image; determining a central point of an optic disc based on the black-and-white image, and determining two target branch vessels of the retina based on the black-and-white image; and obtaining the branch angle of the retinal vessel based on the central point of the optic disc and the two target branch vessels.
Based on the retinal blood vessel branch angle calculation method provided by the invention, the retinal fundus image is subjected to binarization processing to obtain a corresponding black-and-white image; on the basis of the black-white image, a target branch blood vessel corresponding to the central point of the optic disc and the retinal blood vessel can be determined, and then the branch angle of the retinal blood vessel can be calculated. The state of the human eye can be determined by the branch angle of the retinal blood vessel.
For example, retinal blood vessels of a patient with severe myopia are generally straight, and the included angle between the middle two main branch blood vessels is smaller, because the axial length of the eye of a patient with myopia is generally higher, the axial length of the eye is longer, the blood vessels are straightened, and finally the included angle between the middle two main branch blood vessels is smaller.
The following describes the retinal blood vessel branch angle calculation device provided by the present invention, and the retinal blood vessel branch angle calculation device described below and the retinal blood vessel branch angle calculation method described above may be referred to in correspondence with each other.
As shown in fig. 5, the retinal blood vessel branch angle calculation apparatus 500 provided by the present invention includes: an image processing module 510, a feature determination module 520, and an angle calculation module 530.
The image processing module 510 is configured to obtain a retinal fundus image, and perform binarization processing on the retinal fundus image to obtain a black-and-white image.
The feature determination module 520 is configured to determine a central point of the optic disc based on the black-and-white image, and determine two target branch vessels of the retina based on the black-and-white image.
The angle calculation module 530 is configured to obtain a retinal vessel branch angle based on the central point of the optic disc and the two target branch vessels.
In some embodiments, the feature determination module 520 includes: a first determination unit and a second determination unit.
The first determining unit is used for determining the optic disc based on the black-and-white image.
The second determining unit is used for determining the position of a target pixel point in the video disc and determining the central point of the video disc based on the position of the target pixel point in the video disc.
In some embodiments, the first determination unit comprises: a first traversal unit and a disc determination unit.
The first traversal unit is used for traversing the black-and-white image based on a preset first traversal window; wherein the first traversal window is a circular window.
And the video disc determining unit is used for taking the area with the maximum concentration of the target pixel points in the first traversal window as the video disc.
In some embodiments, the second determination unit comprises: a first coordinate system creation unit and a disc center calculation unit.
The first coordinate system creating unit is used for creating a two-dimensional coordinate system on the black and white image;
the optic disc center calculating unit is used for calculating the position of the optic disc center point based on the following formula under the two-dimensional coordinate system:
Figure BDA0003345323850000101
wherein x isp、ypRespectively an abscissa and an ordinate of the central point of the optic disc,
Figure BDA0003345323850000102
respectively the mean value of the abscissa and the mean value of the ordinate, x, of all target pixel points in the black-and-white imagei、yiRespectively is the abscissa and the ordinate of the ith target pixel point in the black-and-white image.
In some embodiments, the feature determination module 520 further comprises: a second coordinate creating unit and a traversal calculating unit.
The second coordinate creating unit is used for creating a polar coordinate system based on the central point of the optic disc.
And the traversal calculation unit is used for traversing the black-and-white image by rotating a preset second traversal window under the polar coordinate system to determine two target branch vessels.
The second traversal window is a trilateral, and the vertex of the trilateral is overlapped with the central point of the optic disc.
In some embodiments, traversing the computing unit comprises: a second traversal unit and an angle calculation unit.
And the second traversal unit is used for traversing the black-and-white image by rotating a preset second traversal window under the polar coordinate system, and respectively taking two areas with the largest target pixel points in the second traversal window as two target branch vessels.
The angle calculation unit is used for determining the included angles between the two target branch vessels and the central point of the optic disc respectively based on the positions of the retinal vessels and the central point of the optic disc.
The electronic device, the computer program product, and the storage medium provided by the present invention are described below, and the electronic device, the computer program product, and the storage medium described below may be referred to in correspondence with the above-described retinal blood vessel branch angle calculation method.
Fig. 6 illustrates a physical structure diagram of an electronic device, which may include, as shown in fig. 6: a processor (processor)610, a communication Interface (Communications Interface)620, a memory (memory)630 and a communication bus 640, wherein the processor 610, the communication Interface 620 and the memory 630 communicate with each other via the communication bus 640. Processor 610 may invoke logic instructions in memory 630 to perform a retinal vessel branch angle calculation method comprising:
110, obtaining a retina fundus image, and performing binarization processing on the retina fundus image to obtain a black-and-white image;
step 120, determining a central point of the optic disc based on the black-and-white image, and determining two target branch vessels of the retina based on the black-and-white image;
and step 130, obtaining a retinal vessel branch angle based on the central point of the optic disc and the two target branch vessels.
In addition, the logic instructions in the memory 630 may be implemented in software functional units and stored in a computer readable storage medium when the logic instructions are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, the present invention also provides a computer program product, the computer program product including a computer program, the computer program being stored on a non-transitory computer-readable storage medium, wherein when the computer program is executed by a processor, a computer is capable of executing the method for calculating a retinal vessel branch angle provided by the above methods, the method including:
110, obtaining a retina fundus image, and performing binarization processing on the retina fundus image to obtain a black-and-white image;
step 120, determining a central point of the optic disc based on the black-and-white image, and determining two target branch vessels of the retina based on the black-and-white image;
and step 130, obtaining a retinal vessel branch angle based on the central point of the optic disc and the two target branch vessels.
In yet another aspect, the present invention also provides a non-transitory computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method for calculating a retinal vessel branch angle provided by the above methods, the method including:
110, obtaining a retina fundus image, and performing binarization processing on the retina fundus image to obtain a black-and-white image;
step 120, determining a central point of the optic disc based on the black-and-white image, and determining two target branch vessels of the retina based on the black-and-white image;
and step 130, obtaining a retinal vessel branch angle based on the central point of the optic disc and the two target branch vessels.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (10)

1. A retinal blood vessel branch angle calculation method is characterized by comprising the following steps:
obtaining a retina fundus image, and performing binarization processing on the retina fundus image to obtain a black-and-white image;
determining a central point of an optic disc based on the black-and-white image, and determining two target branch vessels of the retina based on the black-and-white image;
and obtaining the branch angle of the retinal vessel based on the central point of the optic disc and the two target branch vessels.
2. The method of claim 1, wherein the determining a disk center point based on the black-and-white image comprises:
determining a video disc based on the black and white image;
and determining the position of a target pixel point in the optic disc, and determining the central point of the optic disc based on the position of the target pixel point in the optic disc.
3. The retinal blood vessel branch angle calculation method according to claim 2, wherein the determining an optic disc based on the black-and-white image includes:
traversing the black-and-white image based on a preset first traversal window; wherein the first traversal window is a circular window;
and taking the area with the maximum concentration of the target pixel points in the first traversal window as the video disc.
4. The retinal vessel branch angle calculation method according to claim 2, wherein the determining a disc center point based on a target pixel point position in the disc comprises:
creating a two-dimensional coordinate system on the black-and-white image;
under the two-dimensional coordinate system, the position of the central point of the optic disc is calculated and obtained based on the following formula:
Figure FDA0003345323840000011
wherein x isp、ypRespectively an abscissa and an ordinate of the central point of the optic disc,
Figure FDA0003345323840000021
respectively the mean value of the abscissa and the mean value of the ordinate, x, of all target pixel points in the black-and-white imagei、yiRespectively is the abscissa and the ordinate of the ith target pixel point in the black-and-white image.
5. The retinal blood vessel branch angle calculation method according to any one of claims 1 to 4, wherein the two target branch blood vessels of the retina are determined based on the black-and-white image:
based on the central point of the optic disc, a polar coordinate system is established;
under the polar coordinate system, a preset second traversal window is rotated to traverse the black-and-white image, and two target branch vessels are determined;
the second traversal window is a trilateral, and the vertex of the trilateral is overlapped with the central point of the optic disc.
6. The method for calculating the branch angle of retinal blood vessel according to claim 5, wherein the traversing the black-and-white image by rotating a preset second traversal window in the polar coordinate system to determine two target branch blood vessels comprises:
and under the polar coordinate system, traversing the black-and-white image by rotating a preset second traversal window, and taking two areas with the largest target pixel points in the second traversal window as two target branch vessels respectively.
7. A retinal blood vessel branch angle calculation apparatus, comprising:
the image processing module is used for acquiring a retina fundus image and carrying out binarization processing on the retina fundus image to obtain a black-and-white image;
the characteristic determining module is used for determining the central point of the optic disc based on the black-and-white image and determining two target branch blood vessels of the retina based on the black-and-white image;
and the angle calculation module is used for obtaining the retinal blood vessel branch angle based on the central point of the optic disc and the two target branch blood vessels.
8. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of the retinal vessel branch angle calculation method according to any one of claims 1 to 6 when executing the program.
9. A non-transitory computer-readable storage medium having stored thereon a computer program, wherein the computer program, when executed by a processor, implements the steps of the retinal vessel branch angle calculation method according to any one of claims 1 to 6.
10. A computer program product comprising a computer program, wherein the computer program when executed by a processor implements the steps of the retinal vessel branch angle calculation method according to any one of claims 1 to 6.
CN202111320203.5A 2021-11-09 2021-11-09 Retinal blood vessel branch angle calculation method and device and electronic equipment Pending CN113951813A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111320203.5A CN113951813A (en) 2021-11-09 2021-11-09 Retinal blood vessel branch angle calculation method and device and electronic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111320203.5A CN113951813A (en) 2021-11-09 2021-11-09 Retinal blood vessel branch angle calculation method and device and electronic equipment

Publications (1)

Publication Number Publication Date
CN113951813A true CN113951813A (en) 2022-01-21

Family

ID=79469862

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111320203.5A Pending CN113951813A (en) 2021-11-09 2021-11-09 Retinal blood vessel branch angle calculation method and device and electronic equipment

Country Status (1)

Country Link
CN (1) CN113951813A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115100222A (en) * 2022-08-24 2022-09-23 首都医科大学附属北京朝阳医院 Image processing method and device for separating artery and vein blood vessels, storage medium and terminal
CN116309391A (en) * 2023-02-20 2023-06-23 依未科技(北京)有限公司 Image processing method and device, electronic equipment and storage medium

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115100222A (en) * 2022-08-24 2022-09-23 首都医科大学附属北京朝阳医院 Image processing method and device for separating artery and vein blood vessels, storage medium and terminal
CN116309391A (en) * 2023-02-20 2023-06-23 依未科技(北京)有限公司 Image processing method and device, electronic equipment and storage medium
CN116309391B (en) * 2023-02-20 2023-09-05 依未科技(北京)有限公司 Image processing method and device, electronic equipment and storage medium

Similar Documents

Publication Publication Date Title
CN113951813A (en) Retinal blood vessel branch angle calculation method and device and electronic equipment
CN109656373B (en) Fixation point positioning method and positioning device, display equipment and storage medium
US20230120295A1 (en) Methods and systems for ocular imaging, diagnosis and prognosis
CN112017187B (en) Method and device for locating center of macula lutea of fundus image, server and storage medium
WO2020019286A1 (en) Blepharoptosis detection method and system
Lee et al. 3-D segmentation of retinal blood vessels in spectral-domain OCT volumes of the optic nerve head
CN114937024A (en) Image evaluation method and device and computer equipment
CN114708280A (en) Cerebrovascular segmentation algorithm applicable to multiple modes
JP5935146B2 (en) Ophthalmic disease image analysis apparatus, ophthalmic image analysis method, and ophthalmic image analysis program
Setiawan et al. Robust pupil localization algorithm based on circular Hough transform for extreme pupil occlusion
EP3861524B1 (en) Method for automatic shape quantification of an optic nerve head
CN114980800A (en) Dioptric pattern generation method and device and computer-readable storage medium
CN111612713B (en) Method for removing occlusion of three-dimensional ultrasonic image
CN108231165A (en) A kind of auxiliary system of eye disease diagnosis
CN114373216A (en) Eye movement tracking method, device, equipment and storage medium for anterior segment OCTA
TWI673034B (en) Methods and system for detecting blepharoptosis
Janpongsri et al. Pseudo‐real‐time retinal layer segmentation for high‐resolution adaptive optics optical coherence tomography
Syga et al. A fully automated 3D in-vivo delineation and shape parameterization of the human lamina cribrosa in optical coherence tomography
Xu et al. Retinal vessel segmentation on SLO image
CN108269258B (en) Method and system for segmenting corneal structures from OCT corneal images
CN110610147A (en) Blood vessel image extraction method, related device and storage equipment
WO2023138273A1 (en) Image enhancement method and system
CN116309594B (en) OCT image processing method for anterior ocular segment
CN116503405B (en) Myopia fundus change visualization method and device, storage medium and electronic equipment
CN116309661B (en) Method for extracting OCT (optical coherence tomography) image contour of anterior segment of eye

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