CN113034608A - Corneal surface morphology measuring device and method - Google Patents
Corneal surface morphology measuring device and method Download PDFInfo
- Publication number
- CN113034608A CN113034608A CN202110264197.XA CN202110264197A CN113034608A CN 113034608 A CN113034608 A CN 113034608A CN 202110264197 A CN202110264197 A CN 202110264197A CN 113034608 A CN113034608 A CN 113034608A
- Authority
- CN
- China
- Prior art keywords
- camera
- cornea
- imaging plane
- point
- coordinates
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 210000004087 cornea Anatomy 0.000 claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 39
- 238000003384 imaging method Methods 0.000 claims description 62
- 238000012545 processing Methods 0.000 claims description 26
- 238000012876 topography Methods 0.000 claims description 23
- 230000003287 optical effect Effects 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000005259 measurement Methods 0.000 claims description 11
- 238000000691 measurement method Methods 0.000 claims description 7
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 210000001508 eye Anatomy 0.000 description 14
- 210000005252 bulbus oculi Anatomy 0.000 description 5
- 238000011161 development Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 201000009310 astigmatism Diseases 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000001356 surgical procedure Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 210000000695 crystalline len Anatomy 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010252 digital analysis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 208000030533 eye disease Diseases 0.000 description 1
- 206010023365 keratopathy Diseases 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000002980 postoperative effect Effects 0.000 description 1
- 201000010041 presbyopia Diseases 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000002054 transplantation Methods 0.000 description 1
- 210000004127 vitreous body Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
- G06T7/85—Stereo camera calibration
-
- 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/107—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for determining the shape or measuring the curvature of the cornea
-
- 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
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/04—Texture mapping
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Surgery (AREA)
- Computer Graphics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ophthalmology & Optometry (AREA)
- Biophysics (AREA)
- Software Systems (AREA)
- Eye Examination Apparatus (AREA)
- Image Processing (AREA)
Abstract
The invention provides a corneal surface morphology measuring device and a method, which are characterized in that firstly, a binocular camera system is used for acquiring an image of a cornea containing dot matrix light spots, then, three-dimensional coordinates of each measuring point are calculated through coordinate conversion, all the three-dimensional coordinate points are used for obtaining an approximate spherical surface where the cornea is located by utilizing a curved surface fitting method, and finally, a Gaussian diopter at each measuring point is calculated, and a corneal topographic map is generated by combining operations such as a data field interpolation method, a color coding method and the like.
Description
Technical Field
The invention relates to the technical field of medical examination, in particular to a corneal surface morphology measuring device and method.
Background
In recent years, corneal topography measurement has become an important part of eye visual optics research and clinical application, and corneal topography is a systematic and comprehensive quantitative analysis means for studying corneal anterior surface morphology, wherein corneal topography is a digital analysis of corneal topography by a computer image processing system, and obtained information is represented by a color map with different characteristics, which is called corneal topography because of the apparent topographic surface height state in geography, and can accurately measure and analyze the curvature of any point on the total corneal anterior surface and detect corneal refractive power. The corneal topography is applied to diagnosis of corneal astigmatism and quantitative analysis of corneal properties, can also be used for preoperative examination and postoperative curative effect evaluation of corneal refractive surgery and corneal transplantation surgery, provides a large amount of corneal morphological information for researchers, and becomes an essential measuring tool for ophthalmology and visual optics.
The cornea topographer most commonly used at present mainly consists of a Placido plate projection system, a real-time image detection system and a computer image processing system, and is finally presented in the form of data or images. The cornea topographer has the advantages of large measuring range, high accuracy and the like. However, the measurement process is very complicated, and requires horizontal and vertical adjustment of the device, so that the first light-emitting point is aligned with the center of the mark and the center of the cornea, which requires a lot of time for adjustment, and the characteristic of which makes it difficult to perform measurement on a patient in a hand-held manner. Obtaining accurate corneal topography is difficult, especially for children who cannot maintain a fixed posture and coordinate measurements.
The newborn eyeball is small and is a physiological presbyopia, and in the growth and development process of the newborn eyeball, the axis of the eye is matched with refractive media such as cornea, crystalline lens, vitreous body and the like, and the mutual influence gradually reaches emmetropization. The cornea curvature variation degree of infants under three years old is large, and especially the cornea curvature degree and the astigmatism degree of infants less than three months are large. Therefore, the detection of the corneal topography is important for the development of eyes of infants, and if the eyes can be found and treated early, the ophthalmic diseases caused by the early finding can be prevented, and the situations of near, far, astigmatism, keratopathy and the like can be effectively reduced. Therefore, the cornea topographer is very important for detecting the condition of the eyeball of the child, and the cornea topographer is an indispensable instrument for detecting the condition of the eyeball.
However, since the conventional corneal topographer requires a complicated measurement process, requires a high level of eyeball condition for the subject, and requires a high level of fitting for the subject, it is difficult to obtain an accurate corneal topography for children who cannot maintain a fixed posture and perform fitting measurement.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a corneal surface morphology measuring device which is characterized by comprising a camera I, a camera II, a dot matrix light source, a calculation processing module and a display module, wherein the camera I, the camera II and the dot matrix light source are respectively and electrically connected with the calculation processing module, and the calculation processing module is electrically connected with the display module;
the lattice light source is used for emitting lattice light, and lattice light spots are formed on the surface of the cornea of the eye as target points;
the camera is used for shooting an eye cornea image with a dot matrix light spot and sending the image to the calculation processing module;
the calculation processing module is used for solving the three-dimensional coordinates of the target points, fitting the three-dimensional coordinates according to the target points to obtain the spherical surface where the cornea is located, then obtaining a corneal topographic map by solving the Gaussian diopter of each target point, and further controlling the on-off state of the dot matrix light source;
the display module is used for displaying the corneal topography.
A measurement method using the corneal topography measuring apparatus, comprising the steps of:
step 1: calibrating parameters and distortion coefficients of a binocular camera system consisting of a camera I and a camera II to obtain an internal parameter matrix K and distortion coefficients of the camera;
step 2: a dot matrix light source is used for emitting dot matrix light, and dot matrix light spots are formed on the surface of the cornea of the eye and serve as target points;
and step 3: respectively shooting the cornea of the eye by a camera I and a camera II which are arranged on the same horizontal plane and have a distance L;
and 4, step 4: performing distortion removal on images shot by the two cameras to obtain two-dimensional coordinates (m, n) of the target point on an imaging plane C1 and two-dimensional coordinates (m ', n') on an imaging plane C2, wherein the imaging plane C1 refers to an imaging plane of the camera I, and the optical center of the imaging plane C1 is recorded as C1The imaging plane C2 refers to the imaging plane of camera II, and the optical center of the imaging plane C2 is denoted as C2;
And 5: with c1As the origin of the three-dimensional coordinate system, with c1、c2Establishing X axis of coordinate system on the straight line, establishing Y axis in the direction perpendicular to the X axis in the imaging plane, establishing Z axis in the direction perpendicular to the normal vector of the imaging plane, and recording optical center c1Has coordinates of (0,0,0) and optical center c2The coordinates of (d) are (L,0, 0);
step 6: respectively calculating the three-dimensional coordinates of the imaging points of the target point on two imaging planes according to the internal parameter matrix K, wherein the imaging point g of the target point on the imaging plane C11The three-dimensional coordinates (X, Y,0) of (a) are calculated as follows:
imaging point g of target point on imaging plane C22The three-dimensional coordinates (X ', Y',0) of (a) are calculated as follows:
and 7: using optical centre c1Image point g1Establish the equation of a straight line L1 using the optical center c2Image point g2Establishing a linear equation L2, and solving an intersection G of the straight lines L1 and L20Then point of intersection G0Coordinate (x) oft,yt,zt) That is, the coordinates of the target point, where the calculated three-dimensional coordinates are respectively expressed as:
wherein F is the focal length of the camera;
and 8: fitting by using a least square method according to the obtained target point coordinates to obtain a spherical equation where the cornea is located;
and step 9: and solving the Gaussian diopter of each target point on the spherical surface to obtain a corneal topography.
The step 8 comprises the following steps:
step 8.1: according to the n target point coordinates { (x)1,y1,z1),(x2,y2,z2),…,(xt,yt,zt),…,(xn,yn,zn) The solving equation for constructing the spherical surface to be fitted is as follows:
in the formula (I), the compound is shown in the specification,a matrix of coefficients is represented by a matrix of coefficients,a representation of the observation vector is shown,representing parameters to be estimated, t is 1,2, …, n, r represents the radius of the sphere to be fitted, and a, b and c represent the parameters to be estimated of the sphere to be fitted;
step 8.2: converting the formula (4) into U ═ ATA)-1ATB, bringing intoAnd solving specific values of a, b, c and r by using the coordinates of the n target points by using a least square method, and determining a determined equation of the spherical surface to be fitted as a spherical surface equation of the cornea.
The step 9 comprises:
step 9.1: rotating a normal plane of the target point intersected with the spherical surface on the spherical surface by taking the normal as a symmetry axis, calculating curvature radiuses along all tangential directions, taking the maximum curvature and the minimum curvature, and converting the maximum curvature and the minimum curvature into diopter;
step 9.2: calculating the average value of the diopter of the maximum curvature and the diopter of the minimum curvature to obtain Gaussian diopter;
step 9.3: repeating the step 9.1 to the step 9.2, and calculating the Gauss diopter of all the target points;
step 9.4: and (3) carrying out grid division on a data field formed by all Gaussian refraction values, solving the data value of an unknown point by a data field interpolation method, and then carrying out color coding processing to finally form a corneal topography.
The invention has the beneficial effects that:
the invention provides a corneal surface morphology measuring device and a method, firstly, an eye cornea image containing lattice facula is obtained through a binocular camera system, then, the three-dimensional coordinates of each measuring point are calculated through coordinate conversion, all the three-dimensional coordinate points are used for obtaining an approximate spherical surface where the cornea is located by utilizing a curved surface fitting method, finally, the Gaussian diopter of each measuring point is calculated, and the corneal topographic map is generated by combining the data field interpolation method, the color coding and other operations. The device can be integrated into a handheld end, is suitable for carrying out the detection of cornea to the infant, conveniently screens early eye disease.
Drawings
FIG. 1 is a block diagram of a corneal topography measuring device according to the present invention;
FIG. 2 is a flow chart of a measurement method based on corneal topography measurement according to the present invention;
FIG. 3 is a schematic diagram of a measurement method according to the present invention;
in the figure, 1, an object point, 2, an optical center of a camera I, 3, an optical center of a camera II, 4, an imaging plane C1, 5 of the camera I, an imaging plane C2, 6 of the camera II, an imaging point of the object point on an imaging plane C1, 7, an imaging point of the object point on an imaging plane C2, 8, a center point of an imaging plane C1, 9, and a center point of an imaging plane C2.
Detailed Description
The invention is further described with reference to the accompanying drawings and the embodiments, which do not require corneal centring, and which are simple and quick to operate and have a lower fitting requirement for infants than the current methods.
As shown in fig. 1, a corneal topography measuring device comprises two cameras I, a camera II, a dot matrix light source, a calculation processing module and a display module, wherein the camera I, the camera II and the dot matrix light source are respectively electrically connected with the calculation processing module, and the calculation processing module is electrically connected with the display module;
the lattice light source is used for emitting lattice light, and lattice light spots are formed on the surface of the cornea of the eye as target points;
the camera is used for shooting an eye cornea image with a dot matrix light spot and sending the image to the calculation processing module;
the calculation processing module is used for solving the three-dimensional coordinates of the target points, fitting the three-dimensional coordinates according to the target points to obtain the spherical surface where the cornea is located, then obtaining a corneal topographic map by solving the Gaussian diopter of each target point, and further controlling the on-off state of the dot matrix light source;
the display module is used for displaying the corneal topography.
The model of the components used in this embodiment is as follows: the two cameras are DH-GPP1080P USB binocular camera modules; the lattice light source is a HW520AD35-16GD lattice laser emission module; the calculation processing module is a Raspberry Pi 3B + development board; the display module is a 4inch HDMI LCD display screen; the method comprises the steps that cornea images are collected from different directions through two cameras and then transmitted to a calculation processing module, the calculation processing module constructs an approximate spherical surface where a cornea is located through three-dimensional coordinates of an extraction target point, a Gaussian diopter without determining the position of the vertex of the cornea is adopted as a numerical form of a cornea topographic map, then color coding processing is carried out to form the cornea topographic map, finally the cornea topographic map is displayed through a display module, the cameras, a dot matrix light source, the calculation processing module and the display module can be integrated into a whole handheld instrument, and a user can complete measurement on a detected person only through the handheld instrument.
The specific measurement principle is as shown in fig. 3, a dot matrix light is emitted by a dot matrix light source, a dot matrix light spot is formed on the surface of a cornea, an image is shot by a camera I and a right camera II at the same time, a calculation processing module processes and calculates the image, the intersection point of the connecting line of a target point 1 and a left light center 2 and a left imaging plane 4 is an image 6 formed by the target point on the left imaging plane 4, the intersection point of the connecting line of the target point 1 and a right light center 3 and a right imaging plane 5 is an image 7 formed by the target point on the right imaging plane 5, and the three-dimensional coordinates of the target point are solved under the condition that the three-dimensional coordinates of the left imaging point 6, the right imaging point 7, the left camera light center 2 and the right camera light center 3 are known.
As shown in fig. 2, a measurement method using a corneal surface morphology measurement device is implemented by using a Raspberry Pi 3B + development board carrying an OpenCV vision development library and matlab software design, and includes:
step 1: calibrating parameters and distortion coefficients of a binocular camera system consisting of a camera I and a camera II to obtain an internal parameter matrix K and distortion coefficients of the camera; during calibration, checkerboards with a fixed format can be printed by matlab software, then the checkerboards are shot by a camera, and then the photos are guided into the matlab software to complete the calibration process through the calibration function.
Step 2: a dot matrix light source is used for emitting dot matrix light, and dot matrix light spots are formed on the surface of the cornea of the eye and serve as target points;
and step 3: respectively shooting the cornea of the eye by a camera I and a camera II which are arranged on the same horizontal plane and have a distance L;
because the detection object is the cornea surface of the eye, after the spot array is irradiated on the cornea surface, the obtained light spot is similar to a square spot array. In an image shot by a camera, firstly, finding out light spots closest to the pixels at the upper left and the lower right of the image according to the coordinates of the light spots in the image, marking the light spots as the angular points at the upper left and the lower right, calculating a regular lattice similar to a lattice light source by taking the two points as the upper left angular point and the lower right angular point of the lattice, traversing the point coordinates of the lattice in the image one by one, finding out the light spots closest to the currently traversed point coordinates, and recording the values of the light spots and the rows and columns in the regular lattice to be the same. The above steps are performed on the images of camera I and camera II, respectively. In the images of the camera I and the camera II, the light spots with the same row and column values are the results of the same target point respectively shot by the camera I and the camera II.
And 4, step 4: carrying out distortion removal processing on images shot by the two cameras by using an OpenCV (open CV) library distortion removal algorithm to obtain two-dimensional coordinates (m, n) of the target point on an imaging plane C1 and two-dimensional coordinates (m ', n') on an imaging plane C2, wherein the imaging plane C1 refers to an imaging plane of the camera I, and the optical center of the imaging plane C1 is recorded as C1The imaging plane C2 refers to the imaging plane of camera II, and the optical center of the imaging plane C2 is denoted as C2;
And 5: with c1As the origin of the three-dimensional coordinate system, with c1、c2Establishing X axis of coordinate system on the straight line, establishing Y axis in the direction perpendicular to the X axis in the imaging plane, establishing Z axis in the direction perpendicular to the normal vector of the imaging plane, and recording optical center c1Has coordinates of (0,0,0) and optical center c2The coordinates of (d) are (L,0, 0);
step 6: respectively calculating the three-dimensional coordinates of the imaging points of the target point on two imaging planes according to the internal parameter matrix K, wherein the imaging point g of the target point on the imaging plane C11The three-dimensional coordinates (X, Y,0) of (a) are calculated as follows:
imaging point g of target point on imaging plane C22The three-dimensional coordinates (X ', Y',0) of (a) are calculated as follows:
and 7: using optical centre c1Image point g1Establish the equation of a straight line L1 using the optical center c2Image point g2Establishing a linear equation L2, and solving an intersection G of the straight lines L1 and L20Then point of intersection G0Coordinate (x) oft,yt,zt) That is, the coordinates of the target point, where the calculated three-dimensional coordinates are respectively expressed as:
wherein, F is the focal length of the camera and is contained in an internal parameter matrix during calibration;
because the cornea surface is approximate to a spherical surface, the point cloud data of the cornea surface can be utilized to carry out spherical surface fitting, namely, the spherical center coordinates and the spherical radius are obtained according to a certain number of spherical points.
The known spherical equation is: r is2=(x-a)2+(y-b)2+(z-c)2
Wherein (a, b, c) are coordinates of the center of sphere, r is the radius of sphere, and (x, y, z) are coordinates of points on the sphere;
the spherical equation is developed as follows:
x2+y2+z2=2xa+2yb+2zc+r2-a2-b2-c2;
in the formula, a, b, c and r2-a2-b2-c2The parameters are regarded as the parameters to be estimated;
and 8: and obtaining a spherical equation of the cornea according to the obtained target point coordinates by using least square fitting, wherein the spherical equation comprises the following steps:
step 8.1: according to the n target point coordinates { (x)1,y1,z1),(x2,y2,z2),…,(xt,yt,zt),…,(xn,yn,zn) Construction of the object to be treatedThe solution equation for fitting the sphere is as follows:
in the formula (I), the compound is shown in the specification,a matrix of coefficients is represented by a matrix of coefficients,a representation of the observation vector is shown,the parameters to be estimated are represented, t is 1,2, …, n, r represents the radius of the spherical surface to be fitted and also the approximate curvature radius of the corneal surface, and a, b and c represent the parameters to be estimated of the spherical surface to be fitted;
step 8.2: converting the formula (4) into U ═ ATA)-1ATAnd B, substituting n target point coordinates, solving specific values of a, B, c and r by using a least square method, and determining a determined equation of the spherical surface to be fitted as a spherical surface equation of the cornea.
And step 9: solving the Gaussian diopter of each target point on the spherical surface to obtain a corneal topography, wherein the corneal topography comprises the following steps:
step 9.1: calculating a normal plane of the target point, which is intersected with the spherical surface on the spherical surface, rotating the normal plane by taking a normal as a symmetry axis, calculating curvature radiuses along all tangential directions, taking a maximum curvature and a minimum curvature, and converting the maximum curvature and the minimum curvature into diopter;
step 9.2: calculating the average value of the diopter of the maximum curvature and the diopter of the minimum curvature to obtain Gaussian diopter;
step 9.3: repeating the step 9.1 to the step 9.2, and calculating the Gauss diopter of all the target points;
step 9.4: and (3) carrying out grid division on a data field formed by all Gaussian refraction values, solving the data value of an unknown point by a data field interpolation method, carrying out color coding processing, identifying different heights of the topographic map, and finally forming the corneal topographic map.
Claims (4)
1. The cornea surface morphology measuring device is characterized by comprising a camera I, a camera II, a dot matrix light source, a calculation processing module and a display module, wherein the camera I, the camera II and the dot matrix light source are respectively and electrically connected with the calculation processing module;
the lattice light source is used for emitting lattice light, and lattice light spots are formed on the surface of the cornea of the eye as target points;
the camera is used for shooting an eye cornea image with a dot matrix light spot and sending the image to the calculation processing module;
the calculation processing module is used for solving the three-dimensional coordinates of the target points, fitting the three-dimensional coordinates according to the target points to obtain the spherical surface where the cornea is located, then obtaining a corneal topographic map by solving the Gaussian diopter of each target point, and further controlling the on-off state of the dot matrix light source;
the display module is used for displaying the corneal topography.
2. A measurement method using the corneal topography measurement apparatus according to claim 1, comprising:
step 1: calibrating parameters and distortion coefficients of a binocular camera system consisting of a camera I and a camera II to obtain an internal parameter matrix K and distortion coefficients of the camera;
step 2: a dot matrix light source is used for emitting dot matrix light, and dot matrix light spots are formed on the surface of the cornea of the eye and serve as target points;
and step 3: respectively shooting the cornea of the eye by a camera I and a camera II which are arranged on the same horizontal plane and have a distance L;
and 4, step 4: performing distortion removal on images shot by the two cameras to obtain two-dimensional coordinates (m, n) of the target point on an imaging plane C1 and two-dimensional coordinates (m ', n') on an imaging plane C2, wherein the imaging plane C1 refers to an imaging plane of the camera I, and the optical center of the imaging plane C1 is recorded as C1The imaging plane C2 refers to the imaging plane of camera II, the optical center of the imaging plane C2Is marked as c2;
And 5: with c1As the origin of the three-dimensional coordinate system, with c1、c2Establishing X axis of coordinate system on the straight line, establishing Y axis in the direction perpendicular to the X axis in the imaging plane, establishing Z axis in the direction perpendicular to the normal vector of the imaging plane, and recording optical center c1Has coordinates of (0,0,0) and optical center c2The coordinates of (d) are (L,0, 0);
step 6: respectively calculating the three-dimensional coordinates of the imaging points of the target point on two imaging planes according to the internal parameter matrix K, wherein the imaging point g of the target point on the imaging plane C11The three-dimensional coordinates (X, Y,0) of (a) are calculated as follows:
imaging point g of target point on imaging plane C22The three-dimensional coordinates (X ', Y',0) of (a) are calculated as follows:
and 7: using optical centre c1Image point g1Establish the equation of a straight line L1 using the optical center c2Image point g2Establishing a linear equation L2, and solving an intersection G of the straight lines L1 and L20Then point of intersection G0Coordinate (x) oft,yt,zt) That is, the coordinates of the target point, where the calculated three-dimensional coordinates are respectively expressed as:
wherein F is the focal length of the camera;
and 8: fitting by using a least square method according to the obtained target point coordinates to obtain a spherical equation where the cornea is located;
and step 9: and solving the Gaussian diopter of each target point on the spherical surface to obtain a corneal topography.
3. The measurement method according to claim 2, wherein the step 8 comprises:
step 8.1: according to the n target point coordinates { (x)1,y1,z1),(x2,y2,z2),…,(xt,yt,zt),…,(xn,yn,zn) The solving equation for constructing the spherical surface to be fitted is as follows:
in the formula (I), the compound is shown in the specification,a matrix of coefficients is represented by a matrix of coefficients,a representation of the observation vector is shown,representing parameters to be estimated, t is 1,2, …, n, r represents the radius of the sphere to be fitted, and a, b and c represent the parameters to be estimated of the sphere to be fitted;
step 8.2: converting the formula (4) into U ═ ATA)-1ATAnd B, substituting n target point coordinates, solving specific values of a, B, c and r by using a least square method, and determining a determined equation of the spherical surface to be fitted as a spherical surface equation of the cornea.
4. The measurement method according to claim 2, wherein the step 9 comprises:
step 9.1: rotating a normal plane of the target point intersected with the spherical surface on the spherical surface by taking the normal as a symmetry axis, calculating curvature radiuses along all tangential directions, taking the maximum curvature and the minimum curvature, and converting the maximum curvature and the minimum curvature into diopter;
step 9.2: calculating the average value of the diopter of the maximum curvature and the diopter of the minimum curvature to obtain Gaussian diopter;
step 9.3: repeating the step 9.1 to the step 9.2, and calculating the Gauss diopter of all the target points;
step 9.4: and (3) carrying out grid division on a data field formed by all Gaussian refraction values, solving the data value of an unknown point by a data field interpolation method, and then carrying out color coding processing to finally form a corneal topography.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110264197.XA CN113034608B (en) | 2021-03-11 | 2021-03-11 | Corneal surface morphology measuring device and method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110264197.XA CN113034608B (en) | 2021-03-11 | 2021-03-11 | Corneal surface morphology measuring device and method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113034608A true CN113034608A (en) | 2021-06-25 |
CN113034608B CN113034608B (en) | 2022-08-23 |
Family
ID=76470339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110264197.XA Active CN113034608B (en) | 2021-03-11 | 2021-03-11 | Corneal surface morphology measuring device and method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113034608B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116485918A (en) * | 2023-06-25 | 2023-07-25 | 天府兴隆湖实验室 | Calibration method, calibration system and computer readable storage medium |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102499633A (en) * | 2011-09-30 | 2012-06-20 | 中国科学院长春光学精密机械与物理研究所 | Large-visual-field liquid-crystal self-adapting optical eyeground imaging method |
CN102508944A (en) * | 2011-10-10 | 2012-06-20 | 沈阳理工大学 | Corneal surface shape fitting method |
CN104688178A (en) * | 2015-03-05 | 2015-06-10 | 温州医科大学 | Keratectasia measurement method based on optical CT |
WO2015093965A1 (en) * | 2013-12-20 | 2015-06-25 | I-Optics B.V. | Determining global shape parameters of a cornea using a two-fold fitting from measurement data to an aspheric toric surface model |
CN105608314A (en) * | 2015-12-17 | 2016-05-25 | 深圳市斯尔顿科技有限公司 | Corneal diopter calculation method and system |
CN108073861A (en) * | 2016-11-16 | 2018-05-25 | 东北大学秦皇岛分校 | A kind of new abnormal gait analysis method and system |
CN108634928A (en) * | 2018-04-23 | 2018-10-12 | 东北大学秦皇岛分校 | A kind of cornea measurement method and system |
CN108652583A (en) * | 2018-02-11 | 2018-10-16 | 温州眼视光发展有限公司 | Detect corneal thickness and the device and method of curvature |
CN109059761A (en) * | 2018-08-01 | 2018-12-21 | 西安交通大学 | A kind of hand-held target head calibration method based on EIV model |
CN109199689A (en) * | 2016-10-13 | 2019-01-15 | 清华大学深圳研究生院 | The cornea therapy device of absorption plant is instilled based on tunable optical component group and solution |
CN109285189A (en) * | 2018-07-12 | 2019-01-29 | 哈尔滨工程大学 | A kind of straight path quick calculation method synchronous without binocular |
CN109696954A (en) * | 2017-10-20 | 2019-04-30 | 中国科学院计算技术研究所 | Eye-controlling focus method, apparatus, equipment and storage medium |
WO2019193529A1 (en) * | 2018-04-06 | 2019-10-10 | Amo Wavefront Sciences, Llc | Methods and systems for corneal topography with in-focus scleral imaging |
CN110974150A (en) * | 2019-12-02 | 2020-04-10 | 四川大学 | Method for measuring human eye cornea topography |
CN111604598A (en) * | 2020-05-06 | 2020-09-01 | 大连理工大学 | Tool setting method of mechanical arm feeding type laser etching system |
CN112022408A (en) * | 2020-07-09 | 2020-12-04 | 常州飞诺医疗技术有限公司 | Method for obtaining human dental arch curve based on binocular vision |
CN112284290A (en) * | 2020-10-20 | 2021-01-29 | 湖南大学 | Autonomous measurement method and system for aero-engine blade robot |
-
2021
- 2021-03-11 CN CN202110264197.XA patent/CN113034608B/en active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102499633A (en) * | 2011-09-30 | 2012-06-20 | 中国科学院长春光学精密机械与物理研究所 | Large-visual-field liquid-crystal self-adapting optical eyeground imaging method |
CN102508944A (en) * | 2011-10-10 | 2012-06-20 | 沈阳理工大学 | Corneal surface shape fitting method |
WO2015093965A1 (en) * | 2013-12-20 | 2015-06-25 | I-Optics B.V. | Determining global shape parameters of a cornea using a two-fold fitting from measurement data to an aspheric toric surface model |
CN104688178A (en) * | 2015-03-05 | 2015-06-10 | 温州医科大学 | Keratectasia measurement method based on optical CT |
CN105608314A (en) * | 2015-12-17 | 2016-05-25 | 深圳市斯尔顿科技有限公司 | Corneal diopter calculation method and system |
CN109199689A (en) * | 2016-10-13 | 2019-01-15 | 清华大学深圳研究生院 | The cornea therapy device of absorption plant is instilled based on tunable optical component group and solution |
CN108073861A (en) * | 2016-11-16 | 2018-05-25 | 东北大学秦皇岛分校 | A kind of new abnormal gait analysis method and system |
CN109696954A (en) * | 2017-10-20 | 2019-04-30 | 中国科学院计算技术研究所 | Eye-controlling focus method, apparatus, equipment and storage medium |
CN108652583A (en) * | 2018-02-11 | 2018-10-16 | 温州眼视光发展有限公司 | Detect corneal thickness and the device and method of curvature |
WO2019193529A1 (en) * | 2018-04-06 | 2019-10-10 | Amo Wavefront Sciences, Llc | Methods and systems for corneal topography with in-focus scleral imaging |
CN108634928A (en) * | 2018-04-23 | 2018-10-12 | 东北大学秦皇岛分校 | A kind of cornea measurement method and system |
CN109285189A (en) * | 2018-07-12 | 2019-01-29 | 哈尔滨工程大学 | A kind of straight path quick calculation method synchronous without binocular |
CN109059761A (en) * | 2018-08-01 | 2018-12-21 | 西安交通大学 | A kind of hand-held target head calibration method based on EIV model |
CN110974150A (en) * | 2019-12-02 | 2020-04-10 | 四川大学 | Method for measuring human eye cornea topography |
CN111604598A (en) * | 2020-05-06 | 2020-09-01 | 大连理工大学 | Tool setting method of mechanical arm feeding type laser etching system |
CN112022408A (en) * | 2020-07-09 | 2020-12-04 | 常州飞诺医疗技术有限公司 | Method for obtaining human dental arch curve based on binocular vision |
CN112284290A (en) * | 2020-10-20 | 2021-01-29 | 湖南大学 | Autonomous measurement method and system for aero-engine blade robot |
Non-Patent Citations (4)
Title |
---|
MICHEL: "Effect of ageing on keratoconic corneas", 《JOURNAL OF OPTOMETRY》 * |
周红仙: "基于谱域光学相干层析的表面形貌成像", 《实验室研究与探索》 * |
应靖璐: "人眼角膜形态研究进展", 《中华眼视光学与视觉科学杂志》 * |
徐丹阳: "角膜地形图仪图像的边缘检测及形貌重构", 《激光与光电子学进展》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116485918A (en) * | 2023-06-25 | 2023-07-25 | 天府兴隆湖实验室 | Calibration method, calibration system and computer readable storage medium |
CN116485918B (en) * | 2023-06-25 | 2023-09-08 | 天府兴隆湖实验室 | Calibration method, calibration system and computer readable storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN113034608B (en) | 2022-08-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20220241056A1 (en) | Monitoring of dentition | |
RU2562690C2 (en) | Individual contact lenses with alignment marks | |
JP5158842B2 (en) | Eye movement measuring method and eye movement measuring apparatus | |
CN108846866B (en) | Optical imaging-based craniomaxillofacial soft tissue sagittal axial median plane determination method and system | |
CN107481228A (en) | Human body back scoliosis angle measurement method based on computer vision | |
JPH10507953A (en) | Corneal analyzer for compound camera | |
JPH01166740A (en) | Contour measuring method and apparatus | |
JP2013000223A (en) | Ophthalmologic apparatus, ophthalmologic system, controlling method for the ophthalmologic apparatus, and program for the controlling method | |
CN111105881B (en) | Database system for 3D measurement of human phenotype | |
CN110441901A (en) | It is a kind of can real-time tracing watch the optical microscope system and method for position attentively | |
CN105066901B (en) | A kind of in vitro survey calculation experimental animal eyeball shape and the method for axis oculi and curvature | |
CN113034608B (en) | Corneal surface morphology measuring device and method | |
US5838811A (en) | System for measuring curved surfaces | |
JP3711053B2 (en) | Line-of-sight measurement device and method, line-of-sight measurement program, and recording medium recording the program | |
CN106618580B (en) | Strabismus and nystagmus head position detection method, device and system | |
CN113012112B (en) | Evaluation system for thrombus detection | |
CN107014592B (en) | Safety goggles visual field detection system and detection method | |
CN105011899B (en) | A kind of laboratory animal eye position radius of curvature and the computational methods of pupil diameter | |
CN114189623B (en) | Light field-based refraction pattern generation method, device, equipment and storage medium | |
CN118078205A (en) | Image processing method, storage medium, and image processing apparatus | |
CN112656366B (en) | Method and system for measuring pupil size in non-contact manner | |
US6556691B1 (en) | System for measuring curved surfaces | |
JP2018166634A (en) | Ophthalmologic device | |
Coghill et al. | Stereo vision based optic nerve head 3D reconstruction using a slit lamp fitted with cameras: performance trial with an eye phantom | |
CN112587083B (en) | Visual processing method, device and computer storage medium |
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 |