CN113034608B - Corneal surface morphology measuring device and method - Google Patents
Corneal surface morphology measuring device and method Download PDFInfo
- Publication number
- CN113034608B CN113034608B CN202110264197.XA CN202110264197A CN113034608B CN 113034608 B CN113034608 B CN 113034608B CN 202110264197 A CN202110264197 A CN 202110264197A CN 113034608 B CN113034608 B CN 113034608B
- Authority
- CN
- China
- Prior art keywords
- camera
- imaging plane
- cornea
- point
- dot matrix
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 210000004087 cornea Anatomy 0.000 claims abstract description 45
- 239000011159 matrix material Substances 0.000 claims abstract description 45
- 238000003384 imaging method Methods 0.000 claims description 61
- 238000012545 processing Methods 0.000 claims description 27
- 238000012876 topography Methods 0.000 claims description 24
- 230000003287 optical effect Effects 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 20
- 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
- 238000005259 measurement Methods 0.000 description 11
- 210000005252 bulbus oculi Anatomy 0.000 description 5
- 238000011161 development Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 238000000691 measurement 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
- 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
-
- G06T5/80—
-
- 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
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 measurement 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 fit 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 also used for 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;
and 2, step: 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 3, 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 C 1 The imaging plane C2 refers to the imaging plane of camera II, and the optical center of the imaging plane C2 is denoted as C 2 ;
And 5: with c 1 As the origin of the three-dimensional coordinate system, with c 1 、c 2 The straight line establishes the X axis of the coordinate system for imagingThe direction perpendicular to the X axis in the plane establishes the Y axis, the direction perpendicular to the normal vector of the imaging plane establishes the Z axis, and the optical center c is recorded 1 Has coordinates of (0,0,0) and optical center c 2 The coordinates of (d) are (L,0, 0);
and 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 C1 1 The three-dimensional coordinates (X, Y,0) of (a) are calculated as follows:
imaging point g of target point on imaging plane C2 2 Is calculated as follows:
and 7: using optical centre c 1 Image point g 1 Establish the equation of a straight line L1 using the optical center c 2 Image point g 2 Establishing a linear equation L2, and solving an intersection G of the straight lines L1 and L2 0 Then point of intersection G 0 Coordinate (x) of (2) t ,y t ,z t ) 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 ,y 1 ,z 1 ),(x 2 ,y 2 ,z 2 ),…,(x t ,y t ,z t ),…,(x n ,y n ,z n ) The solution equation for constructing the sphere 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 ═ A T A) -1 A T And 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.
The step 9 comprises:
step 9.1: rotating a normal plane of the target point, which is intersected with the spherical surface on the spherical surface, by taking a normal line as a symmetry axis, calculating curvature radii 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 faculae 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, the method does not require corneal centring, and compared with the current methods, the method is simple and quick to operate and requires less cooperation with infants.
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 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 to serve 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 adopted 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, dot matrix light is emitted by a dot matrix light source, dot matrix light spots are formed on the surface of the cornea, images are simultaneously shot by a camera I and a camera II, the images are processed and calculated by a calculation processing module, the intersection point of the connecting line of a target point 1 and a left optical 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 optical center 3 and a right imaging plane 5 is an image 7 formed by the target point on a 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 optical center 2 and the right camera optical 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 pictures are led 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 captured 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 asc 1 The imaging plane C2 refers to the imaging plane of camera II, and the optical center of the imaging plane C2 is denoted as C 2 ;
And 5: with c 1 As the origin of the three-dimensional coordinate system, with c 1 、c 2 The X axis of the coordinate system is established on the straight line, the Y axis is established in the direction vertical to the X axis in the imaging plane, the Z axis is established in the direction vertical to the normal vector of the imaging plane, and the light recording center c 1 Has coordinates of (0,0,0) and optical center c 2 The coordinates of (d) are (L,0, 0);
and 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 C1 1 The three-dimensional coordinates (X, Y,0) of (a) are calculated as follows:
imaging point g of target point on imaging plane C2 2 The three-dimensional coordinates (X ', Y',0) of (a) are calculated as follows:
and 7: using optical centre c 1 Image point g 1 Establish the equation of a straight line L1 using the optical center c 2 An image point g 2 Establishing a linear equation L2, and solving an intersection G of the straight lines L1 and L2 0 Then point of intersection G 0 Coordinate (x) of t ,y t ,z t ) 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 is 2 =(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:
x 2 +y 2 +z 2 =2xa+2yb+2zc+r 2 -a 2 -b 2 -c 2 ;
in the formula, a, b, c and r 2 -a 2 -b 2 -c 2 Regarded as the parameter 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 ,y 1 ,z 1 ),(x 2 ,y 2 ,z 2 ),…,(x t ,y t ,z t ),…,(x n ,y n ,z n ) 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,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 ═ A T A) -1 A T And 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 (3)
1. A corneal surface morphology measuring method is realized by adopting a corneal surface morphology measuring device, wherein the measuring device comprises a camera I, a camera II, a dot matrix light source, a calculation processing module and a display module, 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 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 to serve 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; characterized in that the method comprises:
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 by a camera I and a camera II which are arranged on the same horizontal plane and have a distance L between each other;
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 C 1 The imaging plane C2 refers to the imaging plane of camera II, and the optical center of the imaging plane C2 is denoted as C 2 ;
And 5: with c 1 As the origin of the three-dimensional coordinate system, with c 1 、c 2 The X axis of the coordinate system is established on the straight line, the Y axis is established in the direction vertical to the X axis in the imaging plane, the Z axis is established in the direction vertical to the normal vector of the imaging plane, and the light recording center c 1 Has coordinates of (0,0,0) and optical center c 2 The 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 C1 1 The three-dimensional coordinates (X, Y,0) of (a) are calculated as follows:
imaging point g of target point on imaging plane C2 2 The three-dimensional coordinates (X ', Y',0) of (a) are calculated as follows:
and 7: using optical centre c 1 Image point g 1 Establish the equation of a straight line L1 using the optical center c 2 An image point g 2 Establishing a linear equation L2, and solving an intersection G of the straight lines L1 and L2 0 Then point of intersection G 0 Coordinate (x) of (2) t ,y t ,z t ) 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.
2. A corneal topography measuring method according to claim 1, wherein said step 8 comprises:
step 8.1: according to the n target point coordinates { (x) 1 ,y 1 ,z 1 ),(x 2 ,y 2 ,z 2 ),…,(x t ,y t ,z t ),…,(x n ,y n ,z n ) The solution equation for constructing the sphere 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 ═ A T A) -1 A T And 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.
3. A corneal topography measuring method according to claim 2, wherein said 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 CN113034608A (en) | 2021-06-25 |
CN113034608B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116485918B (en) * | 2023-06-25 | 2023-09-08 | 天府兴隆湖实验室 | Calibration method, calibration system and computer readable storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
CN110974150A (en) * | 2019-12-02 | 2020-04-10 | 四川大学 | Method for measuring human eye cornea topography |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102499633B (en) * | 2011-09-30 | 2013-09-25 | 中国科学院长春光学精密机械与物理研究所 | Large-visual-field liquid-crystal self-adapting optical eyeground imaging method |
NL2012004C2 (en) * | 2013-12-20 | 2015-06-26 | Optics B V I | Determining global shape parameters of a cornea using an aspheric toric surface model. |
CN105608314B (en) * | 2015-12-17 | 2018-03-20 | 深圳市斯尔顿科技有限公司 | Determine the method and system of corneal diopter |
CN106264855B (en) * | 2016-10-13 | 2018-11-06 | 清华大学深圳研究生院 | A kind of device for treating cornea tissue disease |
CN108073861B (en) * | 2016-11-16 | 2020-05-15 | 东北大学秦皇岛分校 | Novel abnormal gait analysis method and system |
CN109696954B (en) * | 2017-10-20 | 2021-05-07 | 中国科学院计算技术研究所 | Sight tracking method, device, equipment and storage medium |
CN108652583B (en) * | 2018-02-11 | 2020-12-04 | 温州眼视光发展有限公司 | Device and method for detecting thickness and curvature of cornea |
WO2019193529A1 (en) * | 2018-04-06 | 2019-10-10 | Amo Wavefront Sciences, Llc | Methods and systems for corneal topography with in-focus scleral imaging |
CN108634928B (en) * | 2018-04-23 | 2020-08-11 | 东北大学秦皇岛分校 | Corneal measurement method and system |
CN109285189B (en) * | 2018-07-12 | 2021-10-01 | 哈尔滨工程大学 | Method for quickly calculating straight-line track without binocular synchronization |
CN109059761B (en) * | 2018-08-01 | 2020-06-19 | 西安交通大学 | EIV model-based handheld target measuring head calibration method |
CN111604598B (en) * | 2020-05-06 | 2021-04-20 | 大连理工大学 | 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 |
CN112284290B (en) * | 2020-10-20 | 2021-09-28 | 湖南大学 | Autonomous measurement method and system for aero-engine blade robot |
-
2021
- 2021-03-11 CN CN202110264197.XA patent/CN113034608B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
CN110974150A (en) * | 2019-12-02 | 2020-04-10 | 四川大学 | Method for measuring human eye cornea topography |
Also Published As
Publication number | Publication date |
---|---|
CN113034608A (en) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107481228B (en) | Human back scoliosis angle measuring method based on computer vision | |
EP0843529B1 (en) | Multi-camera corneal analysis system | |
JP5728302B2 (en) | Ophthalmic apparatus, ophthalmic system, ophthalmic apparatus control method, and program for the control method | |
JP2010259605A (en) | Visual line measuring device and visual line measuring program | |
CN108846866B (en) | Optical imaging-based craniomaxillofacial soft tissue sagittal axial median plane determination method and system | |
JPH01166740A (en) | Contour measuring method and apparatus | |
KR20120084656A (en) | Optical coherence tomographic imaging method and optical coherence tomographic imaging apparatus | |
US20170100033A1 (en) | Ophthalmologic apparatus | |
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 | |
JPH08238222A (en) | Apparatus for recognizing line of gaze | |
CN113034608B (en) | Corneal surface morphology measuring device and method | |
CN105066901B (en) | A kind of in vitro survey calculation experimental animal eyeball shape and the method for axis oculi and curvature | |
US5838811A (en) | System for measuring curved surfaces | |
JP4944714B2 (en) | Anterior segment observation device | |
CN111989030B (en) | Image processing method, program, and image processing apparatus | |
JP5893248B2 (en) | Optical tomographic imaging method and optical tomographic imaging apparatus | |
CN114189623B (en) | Light field-based refraction pattern generation method, device, equipment and storage medium | |
CN105011899B (en) | A kind of laboratory animal eye position radius of curvature and the computational methods of pupil diameter | |
CN112656366B (en) | Method and system for measuring pupil size in non-contact manner | |
CN113012112A (en) | Evaluation method and system for thrombus detection | |
US6556691B1 (en) | System for measuring curved surfaces | |
JP2018166634A (en) | Ophthalmologic device | |
CN117598659A (en) | Digital speckle-based cornea measurement method, system, storage medium and processor | |
CN117333506B (en) | Three-dimensional eye movement tracking method based on TOF camera |
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 |