EP4078518A1 - System und verfahren zur erstellung des profils einer augenlinsenkapsel - Google Patents
System und verfahren zur erstellung des profils einer augenlinsenkapselInfo
- Publication number
- EP4078518A1 EP4078518A1 EP20743306.1A EP20743306A EP4078518A1 EP 4078518 A1 EP4078518 A1 EP 4078518A1 EP 20743306 A EP20743306 A EP 20743306A EP 4078518 A1 EP4078518 A1 EP 4078518A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- central
- equatorial
- posterior
- anterior
- controller
- 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
Links
- 239000002775 capsule Substances 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000007704 transition Effects 0.000 claims abstract description 92
- 238000003384 imaging method Methods 0.000 claims abstract description 32
- 210000001747 pupil Anatomy 0.000 claims abstract description 9
- 230000006870 function Effects 0.000 claims description 7
- 230000001131 transforming effect Effects 0.000 claims description 7
- 238000004891 communication Methods 0.000 description 5
- 230000000875 corresponding effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000002350 accommodative effect Effects 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 238000012014 optical coherence tomography Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 208000002177 Cataract Diseases 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001886 ciliary effect Effects 0.000 description 1
- 210000004087 cornea Anatomy 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013523 data management Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000001404 mediated effect Effects 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003319 supportive effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/50—Depth or shape recovery
-
- 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/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- 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/0016—Operational features thereof
- A61B3/0025—Operational features thereof characterised by electronic signal processing, e.g. eye models
-
- 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/1005—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for measuring distances inside the eye, e.g. thickness 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/117—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for examining the anterior chamber or the anterior chamber angle, e.g. gonioscopes
- A61B3/1173—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for examining the anterior chamber or the anterior chamber angle, e.g. gonioscopes for examining the eye lens
-
- 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
- A61B3/145—Arrangements specially adapted for eye photography by video means
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
- A61F2/00—Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
- A61F2/02—Prostheses implantable into the body
- A61F2/14—Eye parts, e.g. lenses, corneal implants; Implanting instruments specially adapted therefor; Artificial eyes
- A61F2/16—Intraocular lenses
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
- A61F9/00—Methods or devices for treatment of the eyes; Devices for putting-in contact lenses; Devices to correct squinting; Apparatus to guide the blind; Protective devices for the eyes, carried on the body or in the hand
- A61F9/007—Methods or devices for eye surgery
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
-
- 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/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- 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
Definitions
- the disclosure relates to a system and method for obtaining a profile of a lens capsule of an eye.
- intraocular lens may be implanted into the eye in a cataract procedure to replace a human lens that has become cloudy. Having a profile of the lens capsule of the eye, prior to the procedure, assists in the selection of the intraocular lens.
- a system having a controller with at least one processor and at least one non-transitory, tangible memory on which instructions are recorded for executing a method for obtaining a profile of a lens capsule of an eye.
- Execution of the instructions by the processor causes the controller to obtain imaging data for a portion of the lens capsule visible through a pupil of the eye.
- the imaging data includes posterior datapoints and anterior datapoints and is transformed to an adjusted frame of reference having a first axis (X) and a second axis (Y).
- X first axis
- Y second axis
- Also disclosed is a corresponding method for obtaining a profile of a lens capsule of an eye
- the profile is represented by respective central surfaces and respective equatorial surfaces separated by respective transition points.
- the controller is configured to fit the imaging data in the adjusted frame of reference to the respective central surfaces in a predefined central region of the lens capsule.
- the method includes obtaining a transition coordinate as a coordinate value of the respective transition points in a positive X-domain.
- the controller is configured to determine a set of fitting parameters for the respective central surfaces and the respective equatorial surfaces based on the transition coordinate and a plurality of constraints.
- the profile is obtained based on the set of fitting parameters for the respective central surfaces and the respective equatorial surfaces.
- the controller may be configured to select an intraocular lens based at least partially on the profile of the lens capsule.
- Transforming the imaging data to the adjusted frame of reference includes fitting the posterior datapoints and the anterior datapoints to a first circle and a second circle, respectively, and determining intersection points of the first circle and the second circle.
- Transforming the imaging data to the adjusted frame of reference includes transforming the posterior datapoints and the anterior datapoints such that they are centered about a respective center of the intersection points and rotated such that a tilt angle of rotation is zero.
- the controller is configured to fit the respective central surfaces in the adjusted frame of reference to respective conic equations for a conic surface.
- the controller is configured to determine a conic x-intercept as a coordinate on the first axis (X) where the respective conic equations intersect in a positive-x domain, the transition coordinate being a product of the conic x-intercept and a predefined constant, the predefined constant being less than 1.
- the respective central surfaces include a central anterior surface.
- the controller may be configured to represent the central anterior surface as an elliptical cone characterized by a first plurality of variables (Ka, Qa, Ra), the first plurality of variables (Ka, Qa, Ra) being obtained by fitting the imaging data to the central anterior surface in the predefined central region in the adjusted frame of reference.
- the respective central surfaces include a central posterior surface.
- the controller may be configured to represent the central posterior surface as an elliptical cone characterized by a second plurality of variables (Kp, Qp, Rp), the second plurality of parameters (Kp, Qp, Rp) being obtained by fitting the imaging data to the central posterior surface in the predefined central region in the adjusted frame of reference.
- the set of fitting parameters include a first anterior parameter (Ga), a second anterior parameter (Pa), a first posterior parameter (Gp), a second posterior parameter (Pp) and respective coordinates (Xe, Ye) of a vertex in the adjusted frame of reference.
- the respective equatorial surfaces include an equatorial anterior surface and an equatorial posterior surface meeting at the vertex.
- the equatorial anterior surface and the equatorial posterior surface are represented by respective skewed parabola functions.
- the equatorial anterior surface is based in part on the first anterior parameter (Ga), the second anterior parameter (Pa) and the respective coordinates (Xe, Ye) of the vertex.
- the equatorial posterior surface is based in part on the first posterior parameter ( Gp ), the second posterior parameter (Pp) and the respective coordinates ( Xe , Ye) of the vertex.
- the equatorial posterior surface Ep(x) may be defined as:
- the respective central surfaces include a central posterior surface and a central anterior surface.
- the respective equatorial surfaces include an equatorial anterior surface and an equatorial posterior surface.
- the plurality of constraints includes a first equation matching respective values of the central anterior surface and the equatorial anterior surface at the transition coordinate, and a second equation matching the respective values of the central posterior surface and the equatorial posterior surface at the transition coordinate.
- the plurality of constraints includes a third equation matching respective first derivatives of the central anterior surface and the equatorial anterior surface at the transition coordinate, and a fourth equation matching the respective first derivatives of the central posterior surface and the equatorial posterior surface at the transition coordinate.
- the plurality of constraints includes a fifth equation matching respective second derivatives of the central anterior surface and the equatorial anterior surface at the transition coordinate, and a sixth equation matching the respective second derivatives of the central posterior surface and the equatorial posterior surface at the transition coordinate.
- the profile is represented by respective central surfaces and respective equatorial surfaces separated by respective transition points.
- Execution of the instructions by the processor causes the controller to obtain a lens diameter and at least two variables from a set of variables, the set of variables including a lens thickness, a central anterior apex and a central posterior apex.
- a transition coordinate is set as a product of the lens diameter and a predefined constant, the predefined constant being less than 0.5.
- the controller is configured to obtain a first plurality of variables (Ka, Qa,
- the controller is configured to obtain a second plurality of variables (Kp, Qp, Rp) and a pair of posterior parameters (Gp, Pp) by simultaneously solving a second group of constraints based in part on the transition coordinate.
- the controller is configured to obtain the profile based on the first plurality of variables (Ka, Qa, Ra), the pair of anterior parameters (Ga, Pa), the second plurality of variables (Kp, Qp, Rp) and the pair of posterior parameter (Gp, Pp).
- An updated value of the lens diameter is obtained based on the profile, and an updated value of the transition coordinate is obtained based on the updated value of the lens diameter.
- the controller is configured to update the first plurality of variables (Ka, Qa, Ra) and the pair of anterior parameters (Ga, Pa) by simultaneously solving the first group of constraints based in part on the updated value of the transition coordinate.
- the controller is configured to update the second plurality of variables (Kp, Qp, Rp) and the pair of posterior parameter ( Gp, Pp) by simultaneously solving the second group of constraints based in part on the updated value of the transition coordinate.
- the respective central surfaces include a central anterior surface and a central posterior surface, the respective equatorial surfaces including an equatorial anterior surface and an equatorial posterior surface.
- the controller is configured to represent the central anterior surface and the central posterior surface as respective elliptical cones characterized by the first plurality of variables (Ka, Qa, Ra) and the second plurality of variables (Kp, Qp, Rp), respectively.
- the controller is configured to represent the equatorial anterior surface and the equatorial posterior surface as respective skewed parabolas characterized by the pair of anterior parameters (Ga, Pa) and the pair of posterior parameters (Gp, Pp), respectively.
- the first group of constraints includes: a first equation matching respective values of the central anterior surface and the equatorial anterior surface at the transition coordinate; a second equation matching respective first derivatives of the central anterior surface and the equatorial anterior surface at the transition coordinate; a third equation matching respective second derivatives of the central anterior surface and the equatorial anterior surface at the transition coordinate; and a fourth equation matching a respective coordinate of the central anterior surface to the central anterior apex.
- the second group of constraints includes: a fifth equation matching the respective values of the central posterior surface and the equatorial posterior surface at the transition coordinate; a sixth equation matching respective first derivatives of the central posterior surface and the equatorial posterior surface at the transition coordinate; a seventh equation matching respective second derivatives of the central posterior surface and the equatorial posterior surface at the transition coordinate; and an eighth equation matching a respective coordinate of the central posterior surface to the central posterior apex.
- FIG. 1 is a schematic illustration of a system for obtaining a profile of a lens capsule of an eye, the system having a controller;
- FIG. 2 is a schematic flowchart for a method executable by the controller of FIG. 1, in accordance with a first embodiment
- FIG. 3 is a schematic fragmentary example of a sectional image of an eye, the image having posterior datapoints and anterior datapoints of a lens capsule;
- FIG. 4 is a schematic illustration of the posterior datapoints and anterior datapoints of FIG. 3, in an original frame of reference;
- FIG. 5. is a schematic illustration of the posterior datapoints and anterior datapoints of FIG. 4, after transformation to an adjusted frame of reference;
- FIG. 6 is a schematic diagram of a profile of the lens capsule obtained by the system of FIG. 1, in the adjusted frame of reference;
- FIG. 7 is a schematic diagram of the profile of the lens capsule shown in
- FIG. 8 is a schematic flowchart for a method executable by the controller of FIG. 1, in accordance with a second embodiment.
- FIG. 1 schematically illustrates system 10 for obtaining a profile of a lens capsule of an eye.
- the system 10 includes a controller C having at least one processor 12 and at least one memory 14 (or non-transitory, tangible computer readable storage medium) on which are recorded instructions for executing one or more methods.
- Method 100 and Method 300 are shown in and described below with reference to FIGS. 2 and 8, respectively.
- the system 10 may include a user interface 16 for collecting user data from one or more clinical facilities or electronic medical record units.
- the system 10 may include a data management unit 18 for storing and/or facilitating transfer of the user data and other functions.
- the various components of the system 10 may be configured to communicate via a short-range network 20 and/or a long-range network 22.
- the controller C may be in communication with a remote server 24 and/or a cloud unit 26, which may include one or more servers hosted on the Internet to store, manage, and process data.
- the cloud unit 26 may be a private or public source of information maintained by an organization, such as for example, a research institute, a company, a university and/or a hospital.
- the short-range network 20 may be a bus implemented in various ways, such as for example, a serial communication bus in the form of a local area network.
- the local area network may include, but is not limited to, a Controller Area Network (CAN), a Controller Area Network with Flexible Data Rate (CAN-FD),
- CAN Controller Area Network
- CAN-FD Controller Area Network with Flexible Data Rate
- the long-range network 22 may be a Wireless Local Area Network (LAN) which links multiple devices using a wireless distribution method, a Wireless Metropolitan Area Networks (MAN) which connects several wireless LANs or a Wireless Wide Area Network (WAN) which covers large areas such as neighboring towns and cities. Other types of connections may be employed.
- LAN Wireless Local Area Network
- MAN Wireless Metropolitan Area Networks
- WAN Wireless Wide Area Network
- the controller C may be configured to receive and transmit wireless communication to the remote server 24 through a mobile application 28, shown in FIG. 1.
- the mobile application 28 may in communication with the controller C via the short- range network 20 such that it has access to the data in the controller C.
- the mobile application 28 is physically connected (e.g. wired) to the controller C.
- the mobile application 28 is embedded in the controller C.
- the circuitry and components of a remote server 24 and mobile application 28 (“apps”) available to those skilled in the art may be employed.
- the user interface 16 and/or the controller C may be configured to communicate with an imaging device 30, which may be an optical coherence tomography machine.
- the imaging device 30 may be an ultrasound machine, a magnetic resonance imaging machine or other imaging device available to those skilled in the art.
- the user interface 16 and/or the controller C may be in communication with a profile output module 32 and a lens selection module 34 for selecting an intraocular lens 36, as described below.
- a profile output module 32 and a lens selection module 34 for selecting an intraocular lens 36, as described below.
- FIG. 3 an example image of an eye E is shown.
- the controller C is configured to obtain imaging data for a portion 140 of a lens capsule 142 visible through the pupil 144 of the eye E.
- the lens capsule 142 has a lens thickness 146 and a lens diameter 148.
- the cornea 152 and iris 154 are shown in FIG. 3 . It is understood that FIG. 3 is not to scale.
- the method 100 includes obtaining imaging data for a portion 140 of the lens capsule 142 visible through the pupil 144 of the eye E.
- Method 100 enables a prediction of the entire shape of the lens capsule 142, including the portions of the shape that are obscured by the iris 154, using the portions 140 of the surface which is visible through the open pupil 144.
- the imaging data may be obtained via ultrasound bio microscopy, optical coherence tomography, magnetic resonance imaging or any other imaging modality available to those skilled in the art.
- the imaging data may be derived from a single image or from multiple images.
- the imaging data may be obtained from imaging device 30.
- the imaging data includes posterior datapoints 160 (at a posterior side P) and anterior datapoints 170 (at an anterior side A).
- the posterior datapoints 160 and anterior datapoints 170 are shown plotted in FIG. 4 in pixel units (referred to herein as the original frame of reference 150), with a horizontal axis Q and a vertical axis R.
- the image shown in FIG. 4 (with posterior side P being above the anterior side A) is flipped vertically compared to FIG. 3 (with anterior side A being above the posterior side P) since the coordinate value in the vertical axis R advances going downwards in FIG. 3.
- the method 100 includes transforming the posterior datapoints 160 and anterior datapoints 170 from the original frame of reference 150 (shown in FIG. 4) to an adjusted frame of reference 200 (shown in FIGS. 5 and 6).
- the posterior datapoints 160 are fitted to the equation for a first circle 162 at the posterior side P.
- the anterior datapoints 170 are fitted to the equation for a second circle 172 at the anterior side A.
- the first circle 162 and the second circle 172 intersect at intersection points 180, 182, shown in FIG. 4.
- the controller C is configured to determine a center 184 between the intersection points 180, 182.
- the coordinate-transformed points (in the adjusted frame of reference 200 of FIG. 5) are obtained by centering the posterior datapoints 160 and the anterior datapoints 170 about the center 184 and rotating the centered coordinates by an angle exactly opposite to a tilt angle 188, such that the resulting tilt angle 188 becomes zero.
- the tilt angle 188 may be obtained as an arctangent of a ratio of the increase in the respective vertical coordinate value divided by an increase in the respective horizontal coordinate value between the intersection points 180, 182.
- the adjusted frame of reference 200 is shown in FIGS. 5-6 and has an X axis and a Y axis.
- the controller C is configured to fit the posterior datapoints 160 and the anterior datapoints 170 in the adjusted frame of reference 200 to respective central surfaces 210 (see FIG. 5) in a predefined central region 205 of the lens capsule 142.
- the predefined central region 205 is in a range 3-7 mm.
- the predefined central region 205 is about 5 mm, corresponding to an average-sized open pupil.
- the predefined central region 205 may be correlated to the opening diameter of the pupil 144 of a specific patient.
- the respective central surfaces 210 are elliptical cones.
- the respective central surfaces 210 may be other types of conic surfaces. It is understood that the form of the respective central surfaces 210 may be varied.
- trace 202 and trace 212 show extrapolation of the respective central surfaces 210 at the posterior side P and the anterior side A, respectively.
- the respective central surfaces 210 are extrapolated to determine a positive conic intersection point 220 and a negative conic intersection point 222.
- the controller is configured to determine a conic x-intercept 224 as a coordinate on the X axis where the respective conic equations intersect in a positive-x domain.
- the conic x-intercept 224 is the X- coordinate of the positive conic intersection point 220.
- the profile L is represented by respective central surfaces 210 and respective equatorial surfaces 240.
- the respective equatorial surfaces 240 include equatorial posterior surfaces 252 and equatorial anterior surfaces 262.
- the central anterior surface 260 is flanked by the equatorial anterior surface 262 on either side.
- the central posterior surface 250 is flanked by the equatorial posterior surface 252 on either side.
- the controller C is configured to obtain a transition coordinate Xt (see FIG. 6) which is the coordinate value of the respective transition points 274 in a positive X-domain.
- the transition coordinate Xt is the positive coordinate on the X axis for the respective transition points 274.
- the predefined constant F is less than 1.
- the predefined constant F is within a range of 0.5 to 0.9.
- the predefined constant F is set to 0.7.
- the central posterior surface 250 and the central anterior surface 260 may be generated from expressions for the conic equations over the range -Xt ⁇ x ⁇ Xt, where Xt is the transition coordinate.
- the central anterior surface 260 is represented as an elliptical cone characterized by a first plurality of variables (Ka, Qa, Ra).
- the first plurality of variables (Ka, Qa, Ra) are obtained by fitting the anterior datapoints 170 (in the adjusted frame of reference 200) in the predefined central region 205.
- the central posterior surface 250 is represented as an elliptical cone characterized by a second plurality of variables (Kp, Qp, Rp).
- the equatorial posterior surfaces 252 and the equatorial anterior surfaces 262 may be represented by respective skewed parabola functions.
- the equatorial anterior surfaces 262 and the equatorial posterior surfaces 252 meet at a vertex 270 in the positive X-domain and another vertex 272 in the negative X- domain.
- the method 100 includes determining a set of fitting parameters for the respective central surfaces 210 and the respective equatorial surfaces 240.
- the set of fitting parameters include a first anterior parameter (Ga), a second anterior parameter (Pa), a first posterior parameter (Gp), a second posterior parameter (Pp) and respective coordinates ( Xe , Ye) of the vertex 270 in the positive X- domain in the adjusted frame of reference 200.
- the equatorial posterior surface 252 (see FIGS. 6 and 7) is based in part on the first posterior parameter (Gp), the second posterior parameter (Pp) and the respective coordinates (Xe, Ye) of the vertex 270 in the positive X-domain.
- the equatorial posterior surface 252 or Ep(x ) is defined as:
- the equatorial anterior surface 262 (see FIGS. 6 and 7) based in part on the first anterior parameter (Ga), the second anterior parameter (Pa) and respective coordinates (Xe, Ye) of the vertex 270.
- the equatorial anterior surface 262 or Ea(x ) is defined as:
- the set of fitting parameters are based on the transition coordinate Xt and a plurality of constraints.
- the plurality of constraints includes first, second, third, fourth, fifth and sixth equations.
- the six constraint equations may be solved numerically, for example, using the MATLAB function vo/ve, employing the trust-region algorithm.
- the first equation matches respective values of the central anterior surface
- the third equation matches respective first derivatives of the central anterior surface 260 and the equatorial anterior surface 262 at the transition coordinate Xt, dCa(x ) dEa(x ) as follows: dx d .
- controller C is configured to obtain the profile
- FIG. 6 shows the profile L after the set of fitting parameters have been obtained. Additionally, per block 112, the profile L may be re-adjusted to the original frame of reference 150 shown in FIG. 4.
- FIG. 7 is a schematic diagram of the profile L after the opposite transformation from that of block 104, including de-centering and re-tilting. In order to provide visualization for clinicians, the profile L in the original frame of reference 150 may be superimposed onto the image of the eye E in FIG. 3.
- the controller C may be configured to select an intraocular lens 36 based at least partially on the profile L of the lens capsule 142. Obtaining an accurate shape of the lens capsule L optimizes selection of the power of the intraocular lens 36. This effect is heightened where the intraocular lens 36 is an accommodative lens which may change its shape in response to external forces. In other words, the intraocular lens 36 may react differently to the same accommodative changes mediated by the ciliary muscles, depending on the geometric dimension and shape of the lens capsule 142.
- the controller C is configured to obtain a lens diameter 148 (see FIG. 3) and at least two variables from a set of variables.
- the set of variables includes a lens thickness 146 (see FIG. 3), respective coordinates on the Y-axis of the central anterior apex 264 (see FIG. 6) and the central posterior apex 254 (see FIG. 6).
- the set of variables may be obtained from the imaging device 30 of FIG. 1 or any other source.
- Method 300 enables a prediction of the entire shape of the lens capsule 142 using a handful of parameters.
- the method 300 includes representing the profile L with respective central surfaces 210 and respective equatorial surfaces 240 in the adjusted frame of reference 200, as shown in FIG. 6.
- the respective central surfaces 210 include a central anterior surface 260 and a central posterior surface 250, represented as elliptical cones characterized by a first plurality of variables (Ka, Qa, Ra) and a second plurality of variables (Kp, Qp, Rp), respectively.
- the respective equatorial surfaces 240 include an equatorial anterior surface 262 and an equatorial posterior surface 252, represented as skewed parabolas characterized by a pair of anterior parameters (Ga, Pa) and a pair of posterior parameters (Gp, Pp), respectively.
- the predefined constant J is less than 0.5. In one example, the predefined constant J is 0.35.
- the controller C is configured to obtain the first plurality of variables (Ka, Qa, Ra) and the pair of anterior parameters (Ga, Pa) by simultaneously solving a first group of constraints based in part on the transition coordinate Xt.
- the first group of constraints includes four equations which may be solved numerically, for example, using the MATLAB function fsolve. Since there are five unknowns and four equations, the Levenberg-Marquardt method may be used. Other numerical algorithms available to those skilled in the art may be employed.
- the first equation matches respective values of the central anterior surface
- the second equation matches respective first derivatives of the central anterior surface 260 and the equatorial anterior surface 262 at the transition dCa(x) dEa(x ) coordinate Xt, as follows: .
- the controller C is configured to obtain the second plurality of variables (Kp, Qp, Rp) and the pair of posterior parameter (Gp, Pp) by simultaneously solving a second group of constraints based in part on the transition coordinate Xt.
- the second group of constraints includes four equations (fifth through eighth equations) which may be solved numerically, for example, using the MATLAB function fsolve. Since there are five unknowns and four equations, the Levenberg- Marquardt method may be used. Other numerical algorithms available to those skilled in the art may be employed.
- the sixth equation matches the respective first derivatives of the central posterior surface 250 and the equatorial posterior surface 252 at the transition dCp(x ) dEp(x) coordinate Xt, as follows: .
- the output of block 306 may be used to obtain the central anterior surface
- the equatorial anterior surface 262 or Ea(x ) is defined as:
- the output of block 308 may be used to obtain the central posterior surface 250 and the equatorial posterior surface 252.
- controller C is configured to obtain the profile
- the controller C is configured to obtain an updated value of the lens diameter 148 based on the profile L, and an updated value of the transition coordinate Xt based on the updated value of the lens diameter 148.
- the controller C is configured to determine if a difference between the updated value of the transition coordinate and the transition coordinate is less than a predefined threshold, in other words, if the updated value of the transition coordinate and the transition coordinate converge to within the predefined threshold. If so, the method 300 is ended.
- the controller C may be configured to select an intraocular lens 36 based at least partially on the profile L of the lens capsule 142.
- the method 300 loops back to block 306 and the controller C is configured to update the first plurality of variables (Ka, Qa, Ra) and the pair of anterior parameters (Ga, Pa) by simultaneously solving the first group of constraints based in part on the updated value of the transition coordinate. Additionally, the controller C is configured to update the second plurality of variables (Kp, Qp, Rp) and the pair of posterior parameter (Gp, Pp) by simultaneously solving the second group of constraints based in part on the updated value of the transition coordinate. [0064] In summary, the system 10 (via execution of the method 100 and/or method 300) enables the prediction of a profile L of the lens capsule 142 with relatively high accuracy while requiring a relatively small number of parameters. The system 10 uses separate parameter values for the anterior side A and the posterior side P of the lens capsule 142, thus capturing the physiologic asymmetric nature of the shape of the lens capsule 142, which can be flatter on one side compared to the other.
- the controller C of FIG. 1 includes a computer-readable medium (also referred to as a processor-readable medium), including a non- transitory (e.g., tangible) medium that participates in providing data (e.g., instructions) that may be read by a computer (e.g., by a processor of a computer).
- a medium may take many forms, including, but not limited to, non-volatile media and volatile media.
- Non-volatile media may include, for example, optical or magnetic disks and other persistent memory.
- Volatile media may include, for example, dynamic random-access memory (DRAM), which may constitute a main memory.
- DRAM dynamic random-access memory
- Such instructions may be transmitted by one or more transmission media, including coaxial cables, copper wire and fiber optics, including the wires that comprise a system bus coupled to a processor of a computer.
- Some forms of computer-readable media include, for example, a floppy disk, a flexible disk, hard disk, magnetic tape, other magnetic medium, a CD-ROM, DVD, other optical medium, punch cards, paper tape, other physical medium with patterns of holes, a RAM, a PROM, an EPROM, a FLASH-EEPROM, other memory chip or cartridge, or other medium from which a computer can read.
- Look-up tables, databases, data repositories or other data stores described herein may include various kinds of mechanisms for storing, accessing, and retrieving various kinds of data, including a hierarchical database, a set of files in a file system, an application database in a proprietary format, a relational database management system (RDBMS), etc.
- Each such data store may be included within a computing device employing a computer operating system such as one of those mentioned above and may be accessed via a network in one or more of a variety of manners.
- a file system may be accessible from a computer operating system and may include files stored in various formats.
- An RDBMS may employ the Structured Query Language (SQL) in addition to a language for creating, storing, editing, and executing stored procedures, such as the PL/SQL language mentioned above.
- SQL Structured Query Language
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Ophthalmology & Optometry (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Vascular Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Geometry (AREA)
- Signal Processing (AREA)
- Multimedia (AREA)
- Cardiology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Transplantation (AREA)
- Prostheses (AREA)
- Eye Examination Apparatus (AREA)
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201962950412P | 2019-12-19 | 2019-12-19 | |
PCT/IB2020/055757 WO2021123926A1 (en) | 2019-12-19 | 2020-06-18 | System and method for obtaining profile of eye lens capsule |
Publications (1)
Publication Number | Publication Date |
---|---|
EP4078518A1 true EP4078518A1 (de) | 2022-10-26 |
Family
ID=71728811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP20743306.1A Pending EP4078518A1 (de) | 2019-12-19 | 2020-06-18 | System und verfahren zur erstellung des profils einer augenlinsenkapsel |
Country Status (7)
Country | Link |
---|---|
US (2) | US11688093B2 (de) |
EP (1) | EP4078518A1 (de) |
JP (1) | JP2023507148A (de) |
CN (1) | CN114846516A (de) |
AU (1) | AU2020405407A1 (de) |
CA (1) | CA3158550A1 (de) |
WO (1) | WO2021123926A1 (de) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4078518A1 (de) * | 2019-12-19 | 2022-10-26 | Alcon Inc. | System und verfahren zur erstellung des profils einer augenlinsenkapsel |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8048155B2 (en) * | 2002-02-02 | 2011-11-01 | Powervision, Inc. | Intraocular implant devices |
US10709610B2 (en) * | 2006-01-20 | 2020-07-14 | Lensar, Inc. | Laser methods and systems for addressing conditions of the lens |
US11185226B2 (en) * | 2008-07-25 | 2021-11-30 | Lensar, Inc. | System and method for measuring tilt in the crystalline lens for laser phaco fragmentation |
BRPI1016052A2 (pt) * | 2009-03-26 | 2020-07-28 | National Digital Research Centre Limited | métodos para modelar uma lente de um olho, e para determinar uma posição ótima para uma lente intraocular, aparelhos para modelar uma lente de um olho, e para determinar uma posição ótima para uma lente intraocular, e, meio de armazenamento legível por computador |
JP5894076B2 (ja) * | 2009-08-31 | 2016-03-23 | パワーヴィジョン・インコーポレーテッド | 水晶体嚢サイズ推定方法 |
JP6502322B2 (ja) * | 2013-04-18 | 2019-04-17 | オプティメディカ コーポレイション | 角膜外科的処置の角膜トポグラフィー測定及びアラインメント |
JP6436888B2 (ja) * | 2015-10-19 | 2018-12-12 | 株式会社トーメーコーポレーション | 眼内レンズ度数決定装置 |
WO2017187267A1 (en) * | 2016-04-29 | 2017-11-02 | Consejo Superior De Investigaciones Científicas (Csic) | Method of estimating a full shape of the crystalline lens from measurements taken by optic imaging techniques and method of estimating an intraocular lens position in a cataract surgery |
EP4078518A1 (de) * | 2019-12-19 | 2022-10-26 | Alcon Inc. | System und verfahren zur erstellung des profils einer augenlinsenkapsel |
-
2020
- 2020-06-18 EP EP20743306.1A patent/EP4078518A1/de active Pending
- 2020-06-18 CN CN202080088745.1A patent/CN114846516A/zh active Pending
- 2020-06-18 US US16/905,693 patent/US11688093B2/en active Active
- 2020-06-18 CA CA3158550A patent/CA3158550A1/en active Pending
- 2020-06-18 WO PCT/IB2020/055757 patent/WO2021123926A1/en unknown
- 2020-06-18 JP JP2022537059A patent/JP2023507148A/ja active Pending
- 2020-06-18 AU AU2020405407A patent/AU2020405407A1/en active Pending
-
2023
- 2023-05-11 US US18/315,591 patent/US20230281856A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
US20210192766A1 (en) | 2021-06-24 |
US11688093B2 (en) | 2023-06-27 |
CA3158550A1 (en) | 2021-06-24 |
AU2020405407A1 (en) | 2022-06-09 |
US20230281856A1 (en) | 2023-09-07 |
CN114846516A (zh) | 2022-08-02 |
JP2023507148A (ja) | 2023-02-21 |
WO2021123926A1 (en) | 2021-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11547484B2 (en) | Selection of intraocular lens based on a plurality of machine learning models | |
US20230281856A1 (en) | System and method for obtaining profile of eye lens capsule | |
WO2022166399A1 (zh) | 基于双模态深度学习的眼底疾病辅助诊断方法和装置 | |
CN108618749B (zh) | 基于便携式数字化眼底照相机的视网膜血管三维重建方法 | |
Zhang et al. | A novel deep learning method for nuclear cataract classification based on anterior segment optical coherence tomography images | |
Lockwood et al. | Lamina cribrosa microarchitecture in normal monkey eyes part 1: methods and initial results | |
Siddiqui et al. | Artificial intelligence in cornea, refractive, and cataract surgery | |
US20220189608A1 (en) | Selection of a preferred intraocular lens based on ray tracing | |
CA3180001A1 (en) | Selection of intraocular lens based on a predicted subjective outcome score | |
Masood et al. | Automatic choroid layer segmentation using normalized graph cut | |
Marsolo et al. | Spatial modeling and classification of corneal shape | |
US11931104B2 (en) | Vision quality assessment based on machine learning model and wavefront analysis | |
JP2023545091A (ja) | 位相シフト構造を有する眼科用レンズとその方法 | |
US12033321B2 (en) | Augmentation of optical coherence tomography image of eye based on learning modules | |
CN111899879A (zh) | 一种眼表疾病自动筛查方法与系统、区块链 | |
CN116823828B (zh) | 黄斑变性程度参数确定方法、装置、设备及存储介质 | |
JP2024151312A (ja) | 顕微鏡撮像モード移行方法、システムおよび記憶媒体 | |
CN116528743A (zh) | 假晶状体眼睛中的术后渐晕的预测 | |
Baroudy et al. | Computer-aided decision system for refractive surgeries with excimer laser |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: UNKNOWN |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE |
|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
|
17P | Request for examination filed |
Effective date: 20220530 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
DAV | Request for validation of the european patent (deleted) | ||
DAX | Request for extension of the european patent (deleted) | ||
P01 | Opt-out of the competence of the unified patent court (upc) registered |
Effective date: 20230508 |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: GRANT OF PATENT IS INTENDED |