CA2204692C - Neural network analysis for multifocal contact lens design - Google Patents

Neural network analysis for multifocal contact lens design Download PDF

Info

Publication number
CA2204692C
CA2204692C CA002204692A CA2204692A CA2204692C CA 2204692 C CA2204692 C CA 2204692C CA 002204692 A CA002204692 A CA 002204692A CA 2204692 A CA2204692 A CA 2204692A CA 2204692 C CA2204692 C CA 2204692C
Authority
CA
Canada
Prior art keywords
neural network
predicted
response
visual acuity
subjective
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.)
Expired - Fee Related
Application number
CA002204692A
Other languages
French (fr)
Other versions
CA2204692A1 (en
Inventor
Jeffrey H. Roffman
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Johnson and Johnson Vision Care Inc
Original Assignee
Johnson and Johnson Vision Care Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Johnson and Johnson Vision Care Inc filed Critical Johnson and Johnson Vision Care Inc
Publication of CA2204692A1 publication Critical patent/CA2204692A1/en
Application granted granted Critical
Publication of CA2204692C publication Critical patent/CA2204692C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02CSPECTACLES; SUNGLASSES OR GOGGLES INSOFAR AS THEY HAVE THE SAME FEATURES AS SPECTACLES; CONTACT LENSES
    • G02C7/00Optical parts
    • G02C7/02Lenses; Lens systems ; Methods of designing lenses
    • G02C7/04Contact lenses for the eyes
    • G02C7/041Contact lenses for the eyes bifocal; multifocal
    • G02C7/044Annular configuration, e.g. pupil tuned
    • GPHYSICS
    • G02OPTICS
    • G02CSPECTACLES; SUNGLASSES OR GOGGLES INSOFAR AS THEY HAVE THE SAME FEATURES AS SPECTACLES; CONTACT LENSES
    • G02C7/00Optical parts
    • G02C7/02Lenses; Lens systems ; Methods of designing lenses
    • G02C7/024Methods of designing ophthalmic lenses
    • G02C7/028Special mathematical design techniques
    • GPHYSICS
    • G02OPTICS
    • G02CSPECTACLES; SUNGLASSES OR GOGGLES INSOFAR AS THEY HAVE THE SAME FEATURES AS SPECTACLES; CONTACT LENSES
    • G02C7/00Optical parts
    • G02C7/02Lenses; Lens systems ; Methods of designing lenses
    • G02C7/04Contact lenses for the eyes
    • G02C7/041Contact lenses for the eyes bifocal; multifocal
    • G02C7/042Simultaneous type

Abstract

The present invention discloses a method for optimizing multifocal lens designs using neural network analysis. More specifically, a neural network is trained using data collected in clinical evaluations of various multifocal lens designs. The trained neural network is then used to predict optimal lens designs for large populations of patients.

Description

NEURAL NETWORK ANALYSIS FOR MULTIFOCAL
CONTACT LENS DESIGN
BACKGROUND OF THE INVENTION
This invention generally relates to multifocal opthamalic lens designs, and more particularly to a method for determining optimal designs for multifocal opthamalic lenses using neural network analysis.
The ability of an individual's eyes to focus on objects may be affected by various optical refractive conditions including presbyopia and cataracts. To elaborate, each eye contains a natural lens that is used to focus images on the retina of the eye. In a person with normal eyesight, the lens of the eye is naturally shaped to focus images of distant objects on the retina, and the eye lens is bent in order to focus images of near objects on the retina.
This adjustment of the eye lens to focus images of objects on the retina is referred to as accommodation.
With many people the accommodation of the eyes is inadequate. For example, as a person ages the natural crystalline lenses of the eyes tend to harden.
The hardening makes it difficult for the eye lens to change its surface curvature (i.e. bending). This inability of the eye lens to bend is a condition referred to as presbyopia. Presbyopia can be corrected with prescribed optical lenses.
Inadequate accommodation may also be caused by cataracts. Many people with cataracts have their natural eye lenses removed and replaced with artificial intraocular lenses which have no ability to 1 change or adjust inside the eye. People with this condition require corrective lenses to compensate for their optical condition.
Inadequate accommodation may be corrected by spectacles or other lenses having a number of different regions with different optical powers. A
person wearing spectacles must shift his or her line of vision so that the object being viewed is observed through a portion of the spectacle having the appropriate optical power. An alternative to spectacles is opthamalic contact lenses with multifocal capabilities. Multifocal contact lenses provide the wearer with the ability to focus on objects at both near and far distances without changing their line of vision.
Specific multifocal lens designs have a circular inner region surrounded by a plurality of concentric annular rings having alternating optical power. The lens is divided into two or more regions, including the circular inner region, having various optical power distributions. Each region, other than the circular inner region, may contain one or more annular rings.
A preferred multifocal lens design for the correction of presbyopia separates the lens into three regions: a circular inner region; an annular intermediate region consisting of multiple annular rings having alternating optical powers; and an annular outer region. This preferred design is described in detail in U.S. Patent No. 5,929,969.
In the preferred design, the inner region has distance optical corrective power. The inner region is most important in conditions of high luminance when the pupil of the eye is constricted to limit the amount of light that enters the eye. When the pupil is constricted the effective corrective coverage area of the multifocal lens is the inner region. High luminance conditions, when the pupil is constricted, usually occur in outdoor daylight activities which require distance vision correction.
Thus, it is advantageous to design multifocal lenses with distance optical power correction in the inner region.
In the preferred design, the intermediate region has both near and distance optical corrective powers. The intermediate region is most important in conditions of intermediate luminance when the pupil is in the middle of its dilation range. When the pupil is in the middle of its dilation range the boundary of the effective corrective coverage area of the multifocal lens is the intermediate region.
Intermediate luminance conditions, when the pupil is in the middle of its dilation range, can occur indoors under artificial light, or outdoors on cloudy days or at dusk. Activities under such conditions can range from reading a book, which requires near vision correction, to driving a vehicle, which requires distance vision correction. Thus, it is advantageous to design multifocal lenses with an intermediate region having both near and distance optical corrective power.
In the preferred design, the outer region is designed to have distance optical corrective power.
The outer region is most important in conditions of low luminance when the pupil of the eye is fully dilated. When the pupil is fully dilated the boundary of the effective corrective coverage area of the multifocal lens is the outer region. Low luminance conditions, when the pupil is fully dilated, can occur outdoors at night. Nighttime activities usually include driving a vehicle or walking down a street which require distance vision correction. Thus, it is advantageous to design multifocal lenses with an outer region having distance optical corrective power.
The multifocal lens design set forth above is a preferred design for correction of presbyopia, however, other designs are possible which may perform as well or better. Optimal designs vary from patient to patient and are dependent upon individual patient parameters.
A recent development in fitting contact lenses is the molding of extended wear and daily wear lenses that are discarded when removed from the eye.
This development requires a high volume automated facility to produce these disposable lenses at an acceptable cost/wear ratio. In this environment, it is impractical to keep small inventories of a large variety of multifocal lens designs that will accommodate small groups of patients. It is more desirable to keep larger inventories of a relatively 1 smaller number of designs that will accommodate large groups of patients. Thus, it is necessary to find optimal designs that will provide proper vision correction for large populations of patients. A
5 complicating factor in the design of lenses for large groups is that subjective performance evaluation does not always allow optimum optical correction. Finding optimal designs for large populations of patients i-s a primary objective of the present invention.
Designing optimal multifocal contact lenses for large groups of patients entails selecting optical design parameters based upon the common optical needs of the patients in the group. Relevant optical design parameters for fitting contact lenses to patients with presbyopia include, but are not limited to, the following: number of annular rings, spacing of annular rings, lens add power, monocular/binocular function, and pupil function. Each of the above-listed optical design parameters is described in the paragraphs that 2 0 follow.
As previously detailed, multifocal contact lenses for treatment of presbyopia generally are comprised of a circular inner region surrounded by one or more annular regions having various optical power 25 distributions. The optical power distribution of each region is defined by a parameter called pupil function. A near or distance optical power distribution indicates that a certain region is predominantly near or distance optical power, 30 respectfully. An equal optical power distribution indicates that a region is comprised of an equal amount of near and distance optical power (i.e.. the annular rings within a region are designed so that half of the total surface area of the region is near optical power and the other half is distance optical power ) .
Although multifocal lenses can be designed with a plurality of regions, practical designs usually have two or three regions with optical power distributions defined by the following pupil function types: distance/equal (d/q), wherein an inner region has distance optical power, and an outer region has an equal amount of near and distance optical power;
distance/equal/distance (d/q/d), wherein an inner region has distance optical power, an intermediate region has an equal amount of near and distance optical power, and an outer region has distance optical power; or near/equal/distance (n/q/d), wherein an inner region has near optical power, an intermediate region has an equal amount of near and distance optical power, and an outer region has distance optical power.
Another optical design parameter to be considered when designing multifocal contact lenses is monocular/binocular function, which is relevant to multifocal lenses as a pair. Monocular/binocular function indicates whether each lens in a pair is .identi.cal or slightly different. Monocular designs are those in which both eyes are fitted with identical lens designs. Binocular designs are those in which each eye is fitted with a slightly different lens designs. Binocular designs usually fit the distance dominant eye with a lens which has a higher proportion of distance optical power than does the lens for the non-dominant eye.
Two very important optical design parameters of a multifocal lens are the number of and spacing between the annular rings that surround the circular inner region. These two parameters are a function of patient pupil size under different illumination conditions. They also are affected by the optical power distribution requirements of each region.
Another important optical design parameter is lens add power which sets the overall add power correction of the multifocal lens. Lens add power is a function of a patients optical corrective needs.
In addition to the optical design parameters discussed above, multifocal lens designs are dependent upon several patient parameters. Relevant patient parameters include, but are not limited to, the following: patient age, patient refractive add and patient Hloss.
As discussed above, the accommodation of a person's eyes diminishes with that person's age. In fact, with advanced age, a person may lose that ability altogether. For this reason alone it is important to consider patient age when designing multifocal contact lenses.
Patient age has additional relevance to multifocal lens designs. The manner in which the size of a person's pupil varies is predictable, principally depending on the illumination level and the age of the person. For people of the same age, the size of their pupils at maximum and minimum dilation change, as a function of illumination level, in the same or substantially the same way. Thus, the size of a person's pupils at minimum and maximum dilation can be estimated based upon the age of that person.
The size of a patient's pupil at different illumination levels affects the effective corrective power ratio of the lens. For example, if a patient is performing a distance vision task under high illumination, the patient's pupil is constricted and the effective corrective area of the multifocal lens is the inner region, which must provide distance vision correction. Thus, since the effective corrective power ratio of a multifocal lens is a function of patient pupil size, and patient pupil size, as a function of illumination, changes with patient age, it is important to consider patient age when designing multifocal lenses.
The inability to focus images that are relatively near (approximately 18 inches) is measured by the amount of positive optical power that must be added to an individual's base distance correction, if any, in order to enable the individual to focus the image properly. The positive optical power that is provided for this reason is typically referred to as refractive add. Multifocal lens designs that will accommodate large groups of patients require that all patients in the particular group have approximately the same refractive add requirement.
A patient's Hloss is a measurement of how many lines of visual acuity a subject loses when 1 trying to perform near vision tasks with their distance Rx compared to their near Rx addition.
Actual clinical performance of specific lens designs are based upon clinical visual acuity measurements of patients while wearing the lenses.
Visual acuity (VA) can be measured in units defined by -10'LogMAR, where MAR is minimum angle of resolution.
In this system 20/20 corresponds to 0.00, values better than 20/20 are positive values, and values poorer than 20/20 are negative values. VA also can be measured in lines lost from a best spectacle correction. Such measurements are based upon number of letters on a line, assuming an eye chart with eight letters on a line.
VA measurements can be taken for near and distance optical performance. VA can also be measured at high, low or intermediate levels of luminance and contrast.
VA performance of specific lens designs can be determined subjectively. Subjective determinations are based upon subjective ratings by individual patients after being fitted with different lens designs. A subjective rating may be based upon a patient's optical performance rating of a lens as defined by a number selected from a predetermined range.
In clinical evaluations, it was found that the optical performance of specific lens designs based upon actual VA measurements did not always coincide with subjective determinations. Thus, it was determined that evaluations of specific lens designs would be more accurate if based upon both actual and subjective VA determinations.
It is not cost effective to inventory large varieties of daily wear or extended wear disposable 5 multifocal lenses having different designs to accommodate individual patients or small groups of patients. If the design is useful only to a small group of patients, then the costs expended to inventory the many different combinations of designs, 10 powers and adds will not be easily recouped. Thus, it is advantageous to design multifocal lenses that will accommodate large populations of patients having similar optical needs.
Finding optimal lens designs that will accommodate large populations of patients is extremely difficult because there are a multitude of design variables to consider. These variables include the optical design parameters and patient parameters discussed above. A complicating factor is that the optical performance of a particular lens design when worn by a particular patient is not linearly related to the optical performance of the lens when worn by other patients. Thus, the relationship between lens design and optical performance is not linear and can not be modeled by a linear equation which could be used to determine optimal designs.
There is a need to formulate a method which can accurately model the relationship between various lens designs and the optical performances of those lenses. Such a model could then be used to determine optimal multifocal lens designs for large populations of patients.

SUMMARY OF THE INVENTION
An object of the present invention is to optimize multifocal lens designs for large populations of patients.
Another object of the present invention is to utilize neural network analysis to determine optimal multifocal lens designs for large populations of patients.
A further object of the present invention is to analyze initial fit Visual Acuity (VA) data collected from clinical studies of various multifocal lens designs and to optimize mutlifocal lens designs based upon those studies.
These and other objectives are achieved by utilizing a neural network to learn true relationships between various multifocal lens designs and their associated visual acuity performance.
According to the above needs and objects of the present invention, there is provided a method of optimizing optical designs involving a plurality of design variables.
The method comprises: (a) identifying predetermined optical design parameters relevant to a predetermined optical refractive condition; (b) forming optical lenses utilizing one or more of said optical design parameters for use in clinical evaluations each such evaluation providing visual acuity data and subjective response ratings for a defined number of patients having said optical refractive condition;
(c) inputting said optical design parameters and related patient parameters as input components, and said visual acuity data and said subjective response ratings as output components into a neural network; (d) training said neural network to model significant relationships between said input components and said output components to thereby produce a trained neural network; (e) isolating and - lla -inputting one or more specific design parameters for evaluation by said trained neural network to predict visual acuity and subjective response as a function of said specific design parameter; and (f) integrating one or more of said predictions to determine the optimal optical design to correct said optical refractive condition.
Further benefits and advantages of the invention will become apparent from a consideration of the following detailed description, which specifies a preferred embodiment of the invention.

BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 shows the structure of a neural network.
Figure 2 shows a schematic representation of a processing element used in a neural network.

1 Figure 3 shows a graph of predicted actual distance VA at high luminance high contrast plotted against number of rings for five different lens designs.
Figure 4 shows a graph of predicted subjective distance response rating plotted against number of rings for five different lens designs.
Figure 5 shows a graph of predicted actual near VA at high luminance high contrast plotted against number of rings for five different lens designs.
Figure 6 shows a graph of predicted subjective near response rating plotted against number of rings for five different lens designs.
Figure 7 shows a graph of a normalized composite of predicted distance VA and predicted subjective distance response plotted against number of rings for five different lens designs.
Figure 8 shows a graph of a normalized composite of predicted near VA and subjective near response plotted against number of rings for five different lens designs.
Figure 9 shows a graph of a normalized composite of predicted near VA/subjective and distance VA/subjective plotted against number of rings for five different lens designs.
Figure 10 shows a graph of a normalized composite of five lens designs relative to predicted VA/subjective response plotted against number of rings.

Neural network analysis is a method of modeling non-linear relationships between independent and dependent variables. The analysis is performed by creating a neural network which accurately models the relationship between the independent and dependent variables. Once a valid neural network is created it can be used to predict values of unknown, dependent variables on the basis of known, independent variables. By convention, in neural network analysis, independent variables are called inputs and dependent variables are called outputs.
The power of a neural network is found in the non-linear equation that it uses to model the relationship between the inputs and the outputs. The equation is a complex function that is defined by a set of variables called connection weights. Connection weights are determined by a training algorithm which examines training data that is input into the neural network. The training data is a set of inputs and associated outputs that are representative of the non-linear relationship being modeled. The training algorithm processes the training data and finds a set of weights that minimize the error between the predicted output of the neural network and the training data output.
A neural network is structurally comprised of an input layer, one or more hidden layers, and an output layer. The output and hidden layers are comprised of interconnected processing elements which are the main building blocks of a neural network.

1 The main function of the input layer is to route input values to processing elements of the first hidden layer. Each processing element multiplies each input by a different weight value and adds the individual products. The results are passed through a non-linear transfer function to produce a processing element output. All processing element outputs of one layer are routed to processing element inputs of the next layer where the same processing is repeated.
The last layer is the output layer which contains either linear or non-linear processing elements. Non-linear processing elements process inputs in the same manner described above. Linear processing elements simply pass the input of the processing element to the output of the processing element. The outputs of the processing elements in the output layer produce the final output of the neural network.
Other neural network design considerations include whether the neural network is a fully connected and/or a feedforward design. A neural network is fully connected if all outputs from one layer are used as inputs to the next layer. A neural network is feedforward if there are no internal feedback loops (i.e. no outputs from one layer are used as inputs to a previous layer).
The first step in creating a neural network is to define what is to be predicted. These predictions will be the outputs of the neural network.
The next step is to identify all variables that could possibly influence the value of the predictions. These 1 variables will be the inputs to the neural network.
Once the inputs and outputs have been identified the remaining structure of the neural network, including the number of layers and the number of processing 5 elements in each layer, must be determined.
Once the structure of the neural network is determined, the neural network can be created. After creation, the neural network is trained using training data. Training data is a set of data, including input 10 variables and associated output variables, which represent the statistical relationship to be modeled by the neural network. The more training data collected and used the better, particularly if the relationship being modeled is statistical.
15 Training is accomplished by a training algorithm that is implemented by the neural network.
The training algorithm processes the training data and selects appropriate connection weights which will most closely model the relationship between the training data inputs and the training data outputs.
After training, the performance of the neural network can be evaluated using test data. Test data is gathered in the same manner as training data, in fact, for superior test results, test data should consist of a random ten percent of the training data.
Testing a neural network is accomplished as follows. Test data inputs are individually input into the neural network. The neural network is run and predicted outputs are generated for each test input.
The predicted outputs are compared to actual test data outputs to determine if the neural network is 1 performing properly. A neural network that performs poorly on test data should not be used.
After a neural network is trained it can be used to predict outputs based on various inputs. A set of variables can be used as inputs and the neural network can be run to predict outputs based upon those inputs. The resulting predictions then can be used for the purpose for which the neural network was designed.
The present invention comprises a method which utilizes a neural network to predict optimal multifocal lens designs for large populations of patients. The neural network is used to model the non-linear relationship between various optical design parameters and the associated visual acuity performance of lenses formed using those design parameters. A detailed description of the inventive method is set forth below.
The first step is to identify optical design parameters that are relevant to a specific optical refractive condition. The identified optical design parameters then are used to form optical lenses that are used in clinical evaluations to derive visual acuity performance data for a defined number of patients having the specified optical refractive condition.
Next, a neural network is created using the identified optical design parameters and various patients parameters as inputs, and the visual acuity data as outputs. After creation, the neural network is trained using the data collected in the clinical evaluations.

1 After training, specific optical design parameters are isolated and evaluated by the trained neural network to predict visual acuity as a function of the specified parameters. The resulting visual acuity predictions are analyzed with reference to the specified design parameters and optimal lens designs are determined for large populations of patients having the optical refractive condition.
In an embodiment of the present inventive method a neural network is created using @BRAIN Neural Network Development System (NNDS) from Talon Development Corporation. The neural network is used to model a non-linear relationship between various optical design parameters and associated V),/subjective performance of lenses formed using those design parameters. The neural network comprises several layers including an input layer, an output layer and one or more hidden layers having interconnected processing elements that are used by a training algorithm to determine a set of connection weights that accurately model the relationship between the optical design parameters and the VA/subjective performance of the associated lenses. The @BRAIN NNDS
utilizes Talon's training algorithm which defines each processing element as a sigmoid non-linear transfer function represented by the following equation:
out = 1.0/ (l+e-'II) - 0.5.
Where out is the output of the processing element, in is the collective sum of each processing element input multiplied by its associated weight, and e- is 1 x 10-iII.

1 In the present embodiment, a neural network was used to determine the optimal number of concentric annular rings for multifocal contact lenses designed for a large population of patients having presbyopia and requiring medium refractive add correction. Each step of the present embodiment of the disclosed inventive method is set forth below.
A first step was to identify optical design parameters for patients having presbyopia and requiring medium add correction. The optical design parameters included number of concentric annular rings, lens add power, pupil function, and monocular or binocular pair. Number of rings had a range from 1 to 12, lens add power had a value of either 1.5 or 2.0, pupil function was either distance/equal (d/q);
distance/equal/distance (d/q/d); or near/equal/distance (n/q/d), and monocular or binocular pair was either monocular pair or binocular pair.
Various combinations of design parameters were selected and labeled with various design codes Each design code indicated a pair of multifocal contact lenses having the same optical design parameters. The various design codes and the associated design parameters were as follows: design code 103 was a monocular design having 12 rings, a lens add power of 1.5 and a pupil function of d/q;
design code 104 was a monocular design having 7 rings, a lens add power of 1.5 and a pupil function of d/q/d;
design code 105/6 was a binocular design having 6.5 rings, a lens add power of 1.5 and a pupil function of d/q; design code 107/8 was a binocular design having 6 rings, a lens add power of 1.5 and a pupil function of d/q/d; design code 109/10 was a binocular design having 3 rings, a lens add power of 1.5 and a pupil function of d/q/d; design code 111/12 was a binocular design having 4.5 rings, a lens add power of 1.5 and a pupil function of d/q; design code 113 was a monocular design having 5 rings, a lens add power of 1.5 and a pupil function of d/q/d; design code 114 was a monocular design having 6 rings, a lens add power of 1.5 and a pupil function of n/q/d; design code 203 was a monocular design having 12 rings, a lens add power of 2 and a pupil function of d/q; and design code 204 was a monocular design having 7 rings, a lens add power of 2 and a pupil function of d/q/d.
Next, multifocal contact lenses were formed using the specific design parameters in each design group. Patients were fitted with the lenses so formed, and clinical evaluations were performed to collect VA
and subjective performance data for each lens design.
For each patient wearing a specified design, actual VA
for near and distance at high luminance and high contrast were measured in visual acuity units (v.a.u.) based upon a -10 Log MAR system. In addition, subjective response ratings for near, distance and overall lens performance were recorded for each patient based upon their satisfaction with the lens in question. Subjective response ratings were selected from a range between 0 and 50, where 50 was best spectacle correction and zero was no correction.

1 The optical design parameters and associated patient parameters were recorded in a spreadsheet along with the corresponding actual VA and subjective response data for each patient wearing each lens 5 design. Patient parameters included patient age, patient refractive add and patient Hloss. The spreadsheet consisted of the data listed in Table I.

v ~
EQ ~4 c v Q > W N O O CO N O tfl r-I O L!1 [ll l0 O O N O
Cl~ O M rn V~ V~ d~ d~ f='1 d~ d~ lll f='1 r1 a' M L(1 lzr r-I
N
N
Q ~4 ~ N c0 O ~ O 01 O Lf1 H V O Ol h l0 O O
-~ ~ f"1 d' M Lfl N L(1 .--I rn V' fn lf1 d' 61 J-1 U]

N C
=n=r-I 1-t UI

:5 m =rl m r+1 1D o lo lo o m N o o t~ O u1 C!] R: 'LS rn N d' d' d' M Ll1 d' m Ll1 m m O 111 S-I r-1 U N C N N c0 ao N CO N
f~ = N = = O = = = ~ = N = O O = O =
N-r-I -r-I O = O r-I = O r-I N = e--I = '-i = = O = r-I
E
~4 ~ UI J0 rl U O O O O O O O O N O O M M O O Ln O
1J fo W N V' O O O O W O = 1D lf1 r-I ~O O Ill L, t11 . . . . . . . . . . . . .
U OJ =.i =~I =rl rl ~i', N= r-I rl r-I O. O O O N e-I r-I ri O.
W
W V~ O O lD N ~ O OD l0 ~ O p l0 N O N ~O
~ . . . . . . . . . . . . . . .
O. N N N V M M l0 O O H
M N

H
S-I O O O L[1 Ll1 O O O lf1 lfl o O o lf1 O O (D
f6 'O O Lf1 tf1 [- [- O O O N N O O O N ul L[1 Lfl a (0 ~I ~-I e-i rl r-I N N N N N
J
C'.
(U
J-I 4) f6 O V~ N l0 Y~1 c0 f~1 (=1 O rl N ~D 111 61 r`'1 O O
2 0 a f6 tr [M d~ d~ lfl d~ Ln Ln Ln Ln w Ln Ln lo N M
II II
r--I '-1 'L3 -rl \
arrrrrr ~\\\
N a,'O '0 G ,~ ~I ~I ,-~ rl r=+ ~I rl ,~ ,~ ~I rl rl r+ ~ ~I ,~
>ti u) Q) U O II II
E
r.r.ar.
O=~I O=~I
E A E,A ~-I rl ri ,-i ,-I rl rl rl H ~i H rl ,-I rl r-A r-I '-+
q o rn ~
m N w v z 3 Ln ui Ln in Ln in In Ln In u) in Ln Ln In Ln Ln in Q) Ri ~
U O m =.~
GL O=.~ N N N N N N N N N N N N N N N N N
O '=ti' 34 rl ,-I f-1 rl rl rl ,-I ,-I ,-i f=I rl ,~ ,-I ,-I 1-1 1-1 r-I
= C q m w Ln ri co ,-I rn H co m rl rn ,-I m ~o m 'O =rl S=I O W H 11 r OD l0 [- l0 ri H N lf1 W H l0 lf1 G' i0 (O [- W L- Lfl f`1 O '-I H r-I [- r-i N m 00 d' f`'1 h c6 ~1 N . . . . . . . . . . .
N a~ ~ o 0 o O o 0 0 o p o O o O o O o 0 ~
=~u) 3~ tl] 'LS M M M (+1 <+1 r+1 M M f+1 (+1 M f`=7 Y~1 P-1 M r+1 4) Q O O O O O O O O O O O O O O O O O

1 r., v o in Ln Ln ,,,O Ln ,o Ln o In rl ~n w In M o o u, GL 0 M N V' M r-I M N M 01 N Q' V' d' <!' L!1 Lfl M

~, co , Ln O lD r \o O 10 Ln Ln o o o o O Ln O
ri r-I r-I m M H (N r-I H O lf1 lf1 m Lf1 lfl M
J W
N C
=n-r~ 1~
Q J-1 Ul ~(Q -rl O t!1 O l- tf1 tf1 d1 H O T O lf1 m lf1 O O O
cJ] a '[J cr M lfl M a' M I:r N M m h M V, V' C M Lfl Lf1 M
S-I ri U 0 O O 47 l0 O 1D N O CO l0 O
L (~ = = = = lp = = = = = = = O O O O O = O
N=~-1 =ri O r-I M O = O ri O N M O H
1 o= O= o= O= O= N =
r-I Sa U1 J-I r--I U O O m W L11 r M4 O L(7 O M W Lf1 LI1 O O O O O O
(~ U1 O O r-I W L~ = O N Lh l0 CD N N O O O O O O
U U) ri rl .,..I . . . . p . . . . . . . . . . . .
O O rl ~-I O O ri N N N N N (N N
N
(q d~ CO ~M O V~ O O l0 ~ CO N W T O
N O . . . . . . . . . . . . O = O O O = O
~4 a) i f ~ N D l0 M V i M L lD Li M d O N
i i i i O O O O i O
J.~

O O L(1 lf1 tf1 lf1 l71 O Lf1 lf1 u-f O O O Lf1 O O t.f) lf1 Lfl Lf1 t- l- [- L- [- O N N N O tf1 Lh L~ O O N N
a 2s . . .
N N N N M O O O r-I ===i r-I r=I
J
F.
-~
J-1 v (d tr1 O N M cM N O O 01 01 01 l0 0~ 00 O l, CD N ~0 lfl ~ L, L, U, L, L, Ln L, lzr ~D lzv M It m M IZV ~t :v a N M
II II
=.i II \ \
~\\\
GLO 'O .(; rl H rl H rl rl H rl H rl rl ri r+ H rl H

N
41 \ U.--IN
U O II II
O =~I O -~I
E 4 E 4 H H e-I H rl rl ,-I rl H H H 1-1 r-I ri ,-I rl 1-1 a rn -~ ~4 m m 0 q g in ui Ln ui Ln Ln Ln Ln Ln in Un Ln Ln Ln Ln Ln Ln Ln Ln Q rwl 0, ~
U O N
tr) 41 . C

O .f, }-I rl ri r-I rl ri r-I rl rl rl rl '-I rl rl rl ri ~-I ri r-1 rl .(.=' C,' M O CO '-I l0 l0 lD lD OJ O ~0 c+1 O 'IV L- r-1 lfl -IV r-I
=ri 7-1 N I'D m [- N m co I- M N OD H Lf1 lD O N lI M ri fo fo Ol h 01 L(7 N O 01 00 M V~ P=1 00 e-i p L(1 r-I Ol M M
~ ~I v Fd 1 . O. O . O O O O O . . . . O. O . O . O. O. O. O O . . O O . .
J r-I O O
L~.
-~ 6) (JJ O M M M M M M M M M M M M M M M M M M M
QJ O p p p p O O O o O o p O p O p o O p O
'O U rl i-i e-i rl ri ri rl rl rl rl rl r-I r-I rl ri ri rl ri rl ~
~
a) O > O lf1 O) O O O L(1 Lfl L(1 Lfl L!1 O O l0 rl T N N N
(,1+ O C C' a' V' a' LL1 N v -tv m m t!1 " m N d~
N
N
~, (Z
> N rn l11 r In t!1 O Ln o Ln Ln Ln O Ln Ln a~ rn m v~ lo =~ ('., M d' V' C' N ~ r-I M IZV N M ~ M N M N [i' y ~
N r.
-n='i JJ
A 4-3 ~l O O L- lf1 O O M GO tf1 O l~ CO 01 OD N O CM ~--I O
U) a ro d' M Lf1 C' C' Lf1 d' tl' d' V' m V' 11 m L(1 N S=1 r-I U O l0 O l0 O O l0 lD N N O
41 fp = O O O= = O = = 0 = = = 0 = = N = = d' ~ N -,y -,y r-I = = '-1 = e-I N = ~-i r-I O = ~--I O = O rl =
O O O i O i ~ O I I i O I ~ O I I O

Q) Sa tz ::5 ~ W J-1 rl U O O O O O O O O O O lf1 O (+1 (+1 O O O O O
~ !!1 O O O O O O O O O O N O H e-1 00 O 10 O O
U 4) -.i ==.i =~I
r-I r-i R(,' x N r-I N r-I r-I O M N r-I r-I O O r-I ~--I O r=i H

z (A O O O O O O V- O N C' O O -IV O O l0 N

p= = N p . N 10 M. m. lf1. '-I r-I . . . . .
y,~ O N M N N
O O O
~
E
~
}-I lfl C. O O O Lf1 Lf1 Lfl Lf) Lf1 Lf1 O L(1 O O O O 111 L(1 F (y~ ro N U1 u1 Lf1 lf1 [, I- h [, h l- O N L-I O lf1 L(7 [- L~
N r-I ~--I r-I ri rl '-I ~--I r-I r.f r-I r-I N N N rl f=f ~--I rl r-I
J
C.' N

16 Z71 00 L- r1 d' c0 tf1 OJ [, H r-I h Q1 L11 O O d' N l0 M
a ro c ~o ~r w w n m w n n n ~o n n a 11 d+ v n N m ,1 r O ZS
-rl II \ \
\\\
(1 `ZS ro G' r-I r--I H .--1 r-I H rl rl H 1-1 r=i H H H N N N N N
N
>y N
L \ U ~ N
U U O II II
E ~ G C L.
~ O = O H
~4 E A E H H H rl H ri H ~-.I rl '-I .-1 rl rl r-I r-I H r-I r-I r-I
ro ~ ro I~
ma) ~ q 3 n Ln in n Ln Ln Ln n In ~n n In n u n In in Ln in r-, ~ a 0 ~

U 0 m -H b) JJ = C
C~a O=rq N N N N N N N N N N N N N N
O q~4 r-I r-i rl rl .-I rl ,-1 H ,-I H H rl r-I ,-I [- r [- L- L, = C L. rl " C 00 Q1 L(1 Ll1 61 LIl y " N 07 l!7 m V, rl 61 l0 l0 -11 l4 01 -IV O CD O 61 L!1 O N N aD m O LI1 O O l, N l0 G," f6 (6 r L- O 01 m N 1~4 lf1 oJ l0 N lD l, m m O
rt ~4 a) 4 11 r-I O O O O O O O O O O O O O O O O O O O
r~
=~ v 0] "CS m M M M M m m M M M M M m m W V' V' U 0 O O O O O O O O O O O O 0 O O C. O O O
'O U rl 1-1 H ri H ri H Ii rl rl H ri lI rl rl H rl ri H

~
~
$4 ' ro 0 > 0 r1 if1 rn O w r - 01 J 0 0 0 OD tD r=1 01 ll1 l(1 CL 0 l!1 C -cr H ll1 'V' fn H C V N m H C N H V fn fA
U) ~, w ro O ~D r-I o O l0 r O LLl Ul 0 m 1D H U1 m Ol ll1 r1 a rl ul sr rl r ,-I C a ~-1 w .=~ C r~ .~ N rl .W N
U Ln UJ G =
n .I a.l 7f6 =.1 0 0 01 l11 O N r 0 O 0 0 if1 r L'1 '+1 r Lf) (p CJ1 a 'd l!1 lIl V' N ll1 m N H rl C' rl f+1 1-1 m p~ 1-1 H U (D O N N m iD N N V' a0 O ~D C
~ rt1 0 0 a~ o 0 RC N ~=~ =~ = N r=i ri r-1 . .-1 _ O = = N ri O = O N 0 O
C. .C .~ O i i i i e-I i i O r=1 i I i ei i i i i > E
4) ri Ya ro ~ M
:3 Cq 1=I r-I U 0 0 0 0 0 0 H O 0 _^. O O m Ln Lr) O 11 r, r-iJ 16 V] 0 O lp l0 0 tfl ~ 0 ~ O O CO r r l0 10 tD lO
u v . =~=~ o=
RC 2: O ~ ~ ~ o ~ o o . . ,-=I -, o N I ,-I 0 0 0 0 ,-I
m 0 ~ o m ~ ~ o 0 0 ~ a ao a o a 0 0 ~c w r=1 f'1 C m ('1 l0 0 = - 11 N
^'1 w 10 m C ri M lf1 l0 a) x I I ~ 0 N
E
rC
"-~ }I 0 0 0 t!1 l77 0 0 - O O O O l!1 6n LI1 LCt U1 O L!1 a 0 O O N N 0 0 _ U1 :7 l!1 lf1 r r r r r O N
f6 N N N N N H H N N

N
L ~
f0 Ol 00 f+1 m O r-1 m ~D = f~1 O O N M C' N O O 01 Q) (a w Ln Ln Ln U) a a ~ in ~ Ln in ui in :n Ln in v c~
N ~ 1 II II
~+~Ivro =~ n~~
rn CL'O 'O C. N N N N (N N N '. N :V N N N N N N N N N
~4 a) 41 U r=I N
u 0 u u E G C C C
ro ~4 0H 0H
m E A E A
rl '-1 .-I .=I rl rl e-I - ,-I -, ri ~ ri '-I r=1 r-I ~-i ri 'i a ,d c 'D
rn ~
.,, >-+
~ rn v r. a ~, ~ ~, ~, ~, ~, ~, _ L, ~ L, ~, L, ~, ~, L, ~, ~, ~, A ri GL .~ =-. ri .~ r-I .-+ ,~ - r-1 -i .-1 ri r1 r1 ri ri rl e-1 r-I
.~
ro w u O m =~
+.1 a 0 =r+
O ~ r r r r r r r r r r r r r r r r r r = G+ t".. N m 1-1 m l0 ll1 r'1 = Ol 01 (D V~ O N QN rl 01 ~-1 m 3 5 17 =.i 1 M V' 0) l0 N r Q~ C 0 L1 r=1 l0 N l0 r-I 01 O lD N
C f6 ro Qh d= P=1 w N M r f'1 :V O) Ltl r c+1 . 01 01 st' U1 V' ro v >4 ,J -4 o 0 0 0 0 0 0 _ o 0 0 0 0 0 0 0 0 0 0 C
~aJ
=o a IV c w a w a c a w er ~r er a c w w w a wo 0 0 0 0 0 0 o c 0 0 0 0 0 0 0 0 0 0 0 'O U ,-1 ri rl H rl H H - '-I =-I rl .=1 ri r-i rl rl ri ri 11 ~ ~4 C v o > o v o in Ln ui in Ln o 0 o r o 0 0 Ln r Ln ¾~ O ~I r M d' M It V~ V~ cM M d, ul lfl N
N
v ~4 ro > N r tf1 Lf1 O O O lf1 O O CO r o o O O O Lfl V' M H d' d' M Ll1 L[1 M d~ N L!1 N L(1 ~

U b1 U L,' =
-n =.I 1~
+1 N
=1 O (~ =r-I O N lfl M tf1 O N lf1 V~ lf7 N r O O O7 O 111 O
1 U) R: 'i~ (N r d' M d' C' C' d= M M dl d' Ln V~ d' C d' Lf) N
Sd r-I U M l0 l0 O O O O lD O
y f~ = = = O O O O O O O = = O O = O = = =
~ Q) =rl -.i O O O = = = = = = = r=I rl = = r-I = r-I N N
O O p o 0 0 o p I o > ~
r1 H
C ~
~ U2 17 ri U OD W Co o p o N o o in o o p p p o p o 0 l.l RS Ul a0 CO CO O O O M O O N O O O O O O O O O
-.I
U (1) .i -1 ~ O O O N N N '-I N N N N ~ N O rl r-I M r-I H
(A
Z OI W 00 CO C O o O O O O O ~
N 0 = = = O = O O O = O = = O O = =
.... S-i rl tf/ M d' = p = = = N = ri ~--I = = N N lD l0 M
O O O ~
~
ro Ln In p p o In 0 p n in n p o p Ln n In Ln F a ~ N N 0 ul ul r 0 O (N N N lfl ul Lf1 ll1 r r r r f6 N N M O O O r-I ~--I rl rl '-i H r-I r-I e--I rl r-I rl r-I
LJ
C
N
iJ N
b=1 61 l0 01 CO O r Oo N 0 l!1 W r r-I V' CO L!1 N r H
a RS d~ l0 d' M C M M 'cM d~ V' d' l0 d' [M [M L(1 M V' L(1 N M
~~vro =.i II
acrtr~r N iln 23 'LS G," N N N N N N N N N N N N N N N N N N N
~4 a) Ai \Ur-IN
a) U O II II
E C C C G -O -H 0 =~
E A E 11 rl r I rl ri rl H e=-I H H H rl H H H ri r-I rl H H
ro a ~ b rn rd -H
m v v q g In Ln Ln Ln in U) In In U) Ln In In In Ln Ln in in Ln Ln a) a ~
U 0 [n 4.) = G
a o =~
O G N r r r r r r r r r r r r r r r r r r r =[," [.; oJ O O M r M M r t0 O o l0 0 l0 l0 r 01 N H
'O =r1 FI Lf1 1=1 l0 1f1 0 01 r r (N w d' r I;v O M d' O l0 O
C, RS r6 V ri c0 OJ N CO N d' r-I r-I N [A l0 l0 ri O co M Sa N
}d 1J ~--I O O O O O O O O O O O O O O O O O O O
G
=rl N
mro ~ v w ~r ~ ~ a v d d v c ~ ~ ~ v ~ v c N 0 0 0 0 0 0 0 0 0 0 0 o p o 0 o p o 0 0 'O U ri .-i '=i ~-I .=i r1 f-I e-i ~-i '-I ~-i ri rl rl r~ ri r-I rl ,~

~
~
v ~
N SI
C O
0 > Ln r, Ln O [- m O m in Ln O ro ~v Ln m I, Ln m in a' 0 V' d' M M r-I C' d' N rl ri M ri N H H M
N
N
N ~
v O N O O~ Lf1 O N O N Oo Lf1 M h lD O O L(1 l0 >
-1 C L(1 d' d' tY' -1 d' M N H lf1 M H N 11 ri rl U 0) N ~
=n=i t~
A ~ N
~3 fo -.i O Lf1 ao o] 61 tf1 ,D t- M O O u'1 O O 1D O O N [, C!) (Y 'O lfl -;I, dl d' C' V' N r-I -zv V~ M N LI1 M rl M N N M
N S=1 rl U C' O O O 07 ~ O O Oo O d+ l0 N
L m Z N-rl =.i N ri .-i rl = ri = ri = = ri r-I O ~-I = O rl = M
FC N o ~ o i H O i i o O
> E
Q) 1-4 ~4 :~ O M
~ N 41 ri U O O O o O O lf1 O lf1 O O lf1 f+1 ri 1.f1 oJ O L!1 M
O o o O O o [- r o lll N r-I = [-OJ 11 h lD
U N -.i O
O O O N N O =--1 N r-I O r-I i H O N N O
G
F N
z N O N d' O O O O ~O N OJ l0 O V~ O O ~ d~ CO
O N p = o M lfl ~--I l0 r-I d~ O = O M l0 N M NI M ri lD ~0 Ul '~$' I i I I I I I O I I I I I I I I I I
a) r W
w ro CL' ~=I lfl lll O O lf1 O O 0 O O O O O O Ll1 Lf1 L() lf1 L!1 2U G[ d L5 L- !- O O N fl O O O LN lf1 f1 Lf1 Lf1 I- t~ t~ N N
a b fp ri r-I N N N
J- N
N
-~
J-~ N
(0 ZJ1 r-I L~ O~ V~ U1 O m lD Ilt O N O M O ~ O M 01 01 a f~ Ln Ln lo U) I.n lr, w d~ -zv m Ill ~O Ln Ln Ln Ln Ill c zv N M
II
-r II \\
a N a'O N N N N N N 11 r-I ri e-I rl r-I H r-I r-I rl rl .-I r-I
~4 N
U rl N
U O II II
E G] C C c 0 =~ O -~
E A E.A H H N N N N N N N N N N N N N
rt a ,d q 'i7 rn ~
-~, ~A
N N
~ C 3 Ln Ln Ill Ln Ln Ln Ln Ln Il n Ln n lfl Ln lf1 lf1 ln Ln Ln a ~
U O N
J~ = C Ln Un Un Ln un Ln un Ln In Ln Ln Ln Ir) a o.~
O C~+ r ~ ~ r r r ~o io ~o ~o ~o ~o o ~o ~o io ~o 0 0 = G' C d~ O ~--I d' d' d~ Ul N e-1 .-1 M O L~ N tfl Ol M r-I r-I
~ N O N U O H M L~ 0 O M 0 l0 CO 0 r CO r ,t'.. (t fo l0 O lfl N l0 CO d' H N N N H H Ol N lfl Ol d' O
~ N
R$ ~1 v }-I 1J ri O O O O O O O O O O O O O O O O O O O
~1 l0 1D l0 l0 l0 lD ~0 l0 l0 ~0 lD l0 l0 -ri N \ \ \ \ \ \ \ \ \ \ \ \ \
ro ' a' sM lf1 LI1 ll1 tf1 I.ft tf1 Ln l!1 lf1 t11 i!1 lf1 tf1 N a Q) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ,~

~
N ~
C v o > 1-4 a ,o Ln ~n o r-W m r Ln m Ln rn io W o Ln [L 0 N C r'1 c'1 N lfl m r'7 r'1 r+1 N l- e-1 m ~--I N r 1 N r1 ~
N
~-, v ft > N r'1 N h I~ O o rl fn O N O Ol Lf1 Ln ln tfl O
r-1 d' r+1 tf1 V' t+1 V' M N L~ .~ ri N N V' ri d' J-1 U]
U bl N ~ =
=n-.i t~
O W T
f0 -.i r+1 V' Lf1 l11 O O lf1 r+1 I~ [ U1 ~ 0, O 01 W
A~ N
U] 10 N a' M N N Lfl m f=1 l+1 M N (~ r-I IZV .-i N t+1 r-I N
71 r-I U ~O O l0 -IV O l0 W lO 'O
y~ ftl = O W N = = = N d+ O = = d~ = = O d+ = O
~. O .rl -rl ~y = = = ri r-I O = = = N O = rl N = = r~ =
N ~=~ ~ ~ o ri e-i >
~4 L1 r-I U m Lfl O Lf1 O M o l71 If1 f~ t.f) f+1 O Lf1 W c+1 tf1 O W
1J (~ U1 l0 N o h o ID O N N r-1 h r-I lfl h W M h O W
U N =rl =rl =~I
O ri r-I r-I r-I O N r-I r-I O r-I O rl ~
F q~
UJ W O ~O N W O d+ O O W W W W O N ~0 d' U N O
y,4 O O f+l lD M N ('1 cM r-I l0 M l0 Lf1 d~ O M N N
~ ],=' 1 I I I I I I I I I I I I I I I 1 I I

v ~
os Sd O O O O O O O U1 Lf1 lll Lf1 L(1 Lfl L(1 O O O O o E a b o 0 o 111 lll Ln In I, L, l- r N N N O o lf1 lfl L(7 r6 Yn r-i .--i ri r-I rl rl '--I rl ri rl N N N M rl ri ri rl G
N
-~
L N
f~ tT1 0~ W ll1 O N fn O 14, N O M l0 01 61 61 W o O O
a ro d~ d~ v lf, U, lf, L(, Ln Ln u, lf1 l0 dl ltv Ln l0 l!1 N ~'1 II II
r-I r-I 'd '~j =r II \\
~~~~
04 'LS 'O C.' '-I N N N N N N N N N N N N N N N N N N
?a N
\ U r-I N
U O II II
o-H o-H
E 4 E ,Q N N N N N N N N N N N N N N N N N N N
a ,0 rn ro ma~
N [~. 3 U) ln vl In Ln Ln Ln In Ln Ln Ln In Ln ln l!l Ln lfl Ln Ln N
Q) 0 ~
U O m =~ ~
y C Ln a o =~ =
O ~~+ n r h =['. .('. lf1 Lf1 N 61 U1 W M 01 W N Ol Ol d' I l0 (") H ~4 N rn M O lO M L, 110 M rl M H N O l, Ln O v [-~ f6 fC W v 01 O ID -0 (+1 V~ Ln N W M 01 W N N lf1 . . . . . . . . . . .
r6 Si N
Sd 1J ri O O O O O O O O O O O O O O O O O O O
C, O O O O
C$'1 lD W W W W W W W W W W W W W W H r-I rl ~-1 =ri v \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
m'LS In c~ L~ I~ h [~ !~ l~ L~ [~ l~ [~ h r L~ rn m rn m ~ o 0 o O o 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 'O U r', ~ ~ ~ ~ rl rl rl .-i rl rl rl rl rl rl ,-i rl rl r~

~
~
~' ro 0 m m o m O rl Ln O O H rl ri C ln L- I m rn Q.i 0 N M M ri rl N r-I V V N N m H M '-I ri r-I O rl m N
~4 N ro N m lf1 O m N O 01 Lf1 O 61 01 O O m C N M N H M N r-I M V~ Q, N N H d' r-I I~ H N N
11 Ul N C =
-n=.1 y A,1 IQ
::s ro -.4 rn m Ln r o H Ln o In o o M o m o rn rn N r-I
Cn fx '0 N rn N M ri N M N N M N M M H H

~4 ri U ~ O O V~ O V~ O er O l0 m = = O = =
O O = = = l0 O N = O C
C U.iy =ri N N H O= = = ri . = O=
r-I ri = rl N
4 N C.. f'. H O I H H i O= H O i > E
r-I la ro a C m 41 ri U m m m M lfl m M m m M m O m m O m O O m LJ f~ !A m m m H [- m l0 m m 1D m O m m O m O O m U N --i -.i =ri r-I r-I O e=i e-i ri .-i O rl N ri r-I O O N

W O O V~ O m V~ m O N lD V~ O C lzv O m ;Zp m m N 0 . . . . . . . . . . . . . . . . . .
~ ri M ri lD M M l0 LI1 O M N N M rl ~D M M l0 L() d~
r t~
N
E
ro C7 S-I tf1 U1 tf1 lf1 Lf1 Lfl L(1 O O O 0 tfl Lfl L(1 Ln Lfl lfl lf1 O
F a [- [- N N N O lf1 L(1 L~ L~
f6 ~I '--I '--I e--I N N N r-I rl '-I ~-i ri ~-i rl ri N N N M

C
N

fo O rn O l0 Ol 01 m O O O lzr O M O l0 01 dl Ol a fro lf, lr, lr, Lr, ~o C C lzv Ln \o lr, I11 lr, Ln Ln 0 v N M
II II
=ri II \ \
~5 a~~~
a'LS 'LS C N N N N N (N N r-I H H H H H H H H H H H
4) 41 U .1 N
w U 0 II II

ro O -HO-~I
E A E A N N N N N N N N N N N (N N N N N N N N
ro rn ro -~, ~, m a) ~ z 3 LIl ln L(1 l11 Ln Ln Ln In ln ln in In l11 In Ln ln Ln Ln Ln N
~ r-I a ~i rl rl ri ei rl rl ri ~i rl ~i ~i r1 e-i .1 .1 ri rl rl U 0 ri) =rl tr LJ = =~'. Ln Ln l() lfl tf) L71 lfl lf1 lfl ll1 Ln 11n O (:~ S-I M M M M M M M d~ V~ cM W q4 = C C V~ O L~ 01 M v Ol Ol [- l0 M m TV 0 N Lf1 Ol m N
'0 ri S=I V' l0 O tf7 Ol [- l0 [- O H H [, lf1 N m 61 N rl ~N
C ro ro O N N lIl m -V Lf7 " m I, cr d' O 01 m N m L~ O
ro Fi N . . . .
4 1J ri O O O O O O O O O O O O O O O O O O O
C O O O O O O O N N N N N N N N N N N N
b=1 r-I ri ri r=I r-I r-I ri ri ~--I r-I r-I ri ri r-I r-I ri r-I f-I f-I
=rl N \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
fA '1'~ Q~ 01 Q1 01 6~ 61 61 r-I r-i r-i e-i ri ri r-I r-I r-I ri r-I ri U 0 O 0 O O O O O r-I rl H '~ U ri 'i rl ri rl ri r=I ri e-i ri ~-1 ri rl rl rl ri ri r1 r1 ~
~ ~4 G O
0 > V' M l0 l0 O Oo 10 lD O Lfl m L(1 Ol O ~-i lD d' N
0 V' M C a' N M M M N M H N rl r-I N M M C, C, N
N
~A
N ro > 4) d' O l0 lf1 L- 00 lD Ill N ~0 O lfl O lfl lfl tf1 L(1 N L!1 M lzv ~T H M M M N N N N ~ rl r-I ri M "I' iJ N
N =
-ri=ri 1~
A ~ N
f6 -r-I d~ h r [- r-I CO N d' O 111 O r L(1 M I~ N I- h p C/) p.', M N M M M M M N H H N M H M a' S4 r-I U d' O OD n7 cr O ~ O N ~D
y~ f6 O ID 0 = O I'D [O = O = O l0 Q) .,.{ .,.{ . ~ . . ~ p . . . O O = rl O N rl rl =
~C U1 ~=~ =O rl rl rl i i ~-I O O i i o i i i i i rl O
>
~4 ro Ln Ln :3 Ul Jw ri U M N O n7 N O O OJ O lf1 o M c0 O CO O Lf1 c0 M
11 (Q U] l0 = N M = lfl O CO lf1 L- O l0 m O CO o M H
U N -~ -r-I -'-I O O O
1-1 N N i e-I O rl r-I e-i H O O N H r-I i N N
C
F N
2 W O 10 N l0 O 1~ O 10 fA ~ N OD lD N
p O

yõI r-I O O M N M N M d~ rl l0 Lf1 f`=1 ~O Lf1 V~ = O f~1 O x I o= I i , i r y N
ro ~ ~d 0 O 0 0 O O 0 Lh lll Lf1 lf1 O Lf1 lf1 t.f) o o O o H a ~ O O O tfl l11 l.f) lfl I~ [, l~ L- O N N N o O o L(1 N ~I rl ~i r-I r-I ~--I ri r-I rl r-I rl N N N N f`1 ri ri ri JJ
C
U) .~
~ 4) f0 a1 00 l0 Lll p o M o ~ N O M O\ l0 01 Ol 0) ~O Ln p a ro V' d~ Ln 1D L(1 Lfl L7, in Ln Ln lo d' -V v I:p ~ ln N M
II II
=~ u~~
au, rrtr z ~_ -- --GL'O Ra G, N N N N N N N N N N N N N N N N M M M
i-1 U O II II
E G G G C -ro 0 -H O -H
E A E A ~ r+ .1 r-I ,~ ~ ~=-i ~ ,~ ~ r-I rl rl ~ ~ ~ ~ ~ ,-i rt c b o~ ro -r ~
m v ~ G 3 n n in n n n n n in n n Ln n n in In n n In 0 p . . . . . . . . . . . . . . . . . . .
rl a .=~ rl r-I ~-I ,1 i-I rl ~-I ri rl rl rl rl ri rl rl rl ~i ~-I
r-i U O w -H CJl 41 = r~
a 0 -H
O G s4 in n n In n n n Ln n n n n in in ~n in ~o ~o ~o = f' f' l0 rl lD C - N e-i M H ~0 l0 Ol Ill M 01 d' l71 m -H Y-I O N ~O O H Lf1 IV m l0 lD lf1 Ol O N H l- r-I l0 r-I
z ro ro 111 lw a) O O h=7 l0 N [~ H lfl L(1 [M o l0 l0 CO IV O
fZ ~-1 v Sd iJ r-i O O O O O O O O O O O O O o O O O O O
G'.
-~ N
Q1 "d M M M M M t`=1 M M M M M Cl M M M M d~ d' cM
4) 0 'O U rl 'i ri rl ri rl rl e-I ~-i rl rl rl rl rl rl rl rl ri rl ~
~ ~
~ v 0 > 01 fD L(1 01 O O O co l!1 M lo ln ln O o rl O OJ
(L 0 r-f (l) f`-1 N m (N rH r-i N H M rl N V~ H N M ri M
@
N
~4 ~ ro > N r ~o m r rn a~ ~o r in Ln Ln o ~o Ln o co C. N M M N H r-I N r-I .--I [~ M r-I M ~ M e-i M =--I M
L VI
U bl N ~
r~- H i m N ~ c0 l0 O O Ol Il1 O) L71 lD O O M I11 01 A~ N
Cl~ fx 'O N M M N M C r-I e-I M N M 11 N M L(1 N N 01 M
Uj }-I ri U O N O O O CO N ~0 O l0 l0 y~ f{$ = 10 10 = = O O O l0 O = O
~ N -rl -rl N f-I = = ~--I N N O = .-I = rl = = = ri O = O
O O I O r-I O O O
> E
a) r-I ~4 ro 7 O
:~ QI 41 r-I U M lIl lfl N O CA O Ln 0 M LIl O O M O M M OJ O
1J f6 W ~ (N L- N O f0 O h O lD N C, l0 O ~o r-I OJ O
U N =.i =ri =rl rl H H H N O r-I rl '-I ri O O O O O rl O ri fA
z 01 ~D O O N cM CO Oo O N N OJ ~O d~ O O
0 p L' ~4 .--I N M N M l0 M M l0 lf1 ~ O = rl M 1D N N M r-I
U) 4: E
ro }.I O O 0 U1 lf1 U1 L(1 L(1 Ln 0 O O U1 O O O O ll1 t!1 ~ a ~ 11 f1 L(1 I~ h l~ N N N O O O N l(1 lfl L(1 lf1 [-H H r-I ri e-1 r-I (V N N M r-I r-I '-I ri e-i rl rl rl r-I
G
N
-~

f6 01 O M O c~ M O l0 01 61 01 co l0 0) O N O O -ZV O
a lD lr, lr, Ln lfl Ln lo -T ~v w Ln In \o lll lr, I11 N M
II II
r-i r-1 '17 '1~
=.I II \\
N M M M M M M M M M M H H H r-I H r-I r-I r-I ri >N
N
4.) U r-I N

ro 0 =.q O =~
E4 8 A r+ rl ,-I ,~ ri ~-1 H ,~ rl ri .-i r=I rl ,-I rl rl H rl H
ro a ,ti rn ro m a) ~
N G 3 n n Ln n in n n In u n ~ a) 0 Q r-I Qa ~-I ri ri rl ~--I r-I rl ri r=i ri N N N N N N N N N
RS W
U O m =~I b1 a~ = q Oa 0=~I N N N N N N (N N N
O r ~4 ~O ~c lo lo lo lo l0 lo lo lo e-i rl rl rl ri rl ri rl rl =.('. G 01 H tIl l0 OD rl M l0 !- H l0 V~ N M H O L- N N
~=ri 4 ==-1 [- (D 1-4 O l0 0 d' O] I- Ol ri m m N r-I t0 N p C 16 N ~ m M N Ill 61 r-I M V~ H L(1 O H lo 00 -W L- O L(1 . . . . . . . . . . . . . .
ro ~4 Q) N a~ - o O O p o 0 0 0 0 0 o p p O O O O O O
q -~ N
U2 d' d' d' di ~ lzr 11 lzv M M M M M M M M M
N 0 rH rl rH rl H rl 11 1+ H ri O p O O p O O o O
'O U H H H H H H H H e--I '-i N N N N N N (N N N
-ji I I I I I I I I I I I I I I I I I _j ~
~
G N
0 ,> r-I 61 lf1 o Lf1 L(1 M 111 01 0 0 0 0 61 a 0 I- 1-i H M V~ N N M t-i '-i ~ M N rl rl W
N
> v lf1 L(1 N lll 0 m lf1 o C) Lf1 o T tf1 H V' QD
H ~ C r-I N m r-I N Q1 rn -w N N M N H M N N ~0 H

U CSl N C
A +~ N
ru =H ~n lD Orn t~ lo In o M o -w M o t~ ~ o Q1 V] a 'O M M m O7 r-I h r-1 M M ~w N M H N M M ri rl r-I
~ LI ri U O N N o o l0 O l0 ~ O O
U ro o . . . m = o o co . . . v~ . . . o . o 0=r1 =ri = e-i O r-1 = rl = = = N O r=I = O '-i N = ri ~~ o o O o i > 5 r Sd ~5 Ln VI J-I rl U Ill o M M O r= 00 O M o O O ll1 O M o o 47 M
LJ I6 11 [- O ~D l0 tf1 OJ O l0 O O O N u1 lo o tf1 N H
U N =rl =~I =.i o ~C E b C C o 0 0 H H I o ~ .-{ o o N o 0 0 0 o rl 1 m z (!I cM O l0 n7 ~ W c~ o N c0 1D V~ o o ~ O 10 no V~
C 0 . . . . . . . . . . . . . . . . .
ri l0 M lfl M kO lfl V~ O M l0 N N M ~-1 lD M lf1 M l0 v G E
C
u1 t11 O lf1 u1 u1 O O O O o O t11 tf1 C1 u1 O Ln tI1 6 23 L~ L- O N N N o O u'1 lf1 Lf7 t71 [- I~ [- [- O N N
r a ~s ftf ~--I .--I N N N N M r-I e-I ~-i rl rl rl r-I r-I e-i N N N
G
N
J-I N
f6 01 M O 0) l0 Ol 01 01 [O o N o O ~ O M O 01 l0 01 a ro Ln lr, w ~O ~v ~T Ln lr, lo Ln Ln Ln Ln Lr, ~ l0 ~
N M
II II
r-I r-I 'O 'LS
=.i II \ \
arrrrrr ~ \\\
~ a'O 'O C --I r-I .-=I r-I rl H r-I N N N N N N N N N N N N

N
4j \ U .1 N

0 -~=i 0 -~I
EQ E.0 ~i ri ~I ~I ~-i rl rl ~-I ~-I ~i ri ri ,1 rl rl rl rl r-1 ri a ti rn ~
~4 v ~a r-I a N N N N N N N N N N N N N N N N N N N
U 0 m -~I tr ij = G
a 0=r-I N N N N N N N
O ('. }-I ,-I ri r-I .-1 ~-I rl rl L~ t~ t~ t~ t~ t~ t~ t~ t~ t~ t~ L~
=['. ~ c7 co [- N o lf1 LI1 Oo Ln M N H r-i O e-1 [~ l~ O rl =ri ~-I CO M lfl M d' M l0 V, h CO lzv l0 M H lD Lfl M o OJ
C: m r6 Oo L- lD d' m 61 01 h e--I M ul l0 t"1 [, O N ~0 r-I O
tt ~i U
S-1 J.J r/ O o O O O O O O O O O O O O O O O O O
~

W'j~ M M M M M M M
U 0 0 0 0 o o o o 0 o o o o 0 0 0 o 0 o o 'O U N N N N N N N N N N N N N N N N N N N

rn ~4 0 > Ln a a 0 m av ~, > ~ o -,f a a 4J rn U 0) N ~ =
-r~=rl 17 -H `"
cn x ~S ~r rn =~ rl U O
i ~
o=
~. Q) =rl -rl 9 ~
>N
ra a ~ W y rl U t+i U N -.i -.=~ -=a 'D s. ~
W
Ul 00 ~4 v x ' N
E
~
Yi Ln a ~ c'q r6 N

G
N
-~
J-I N
N
a ra v N Cn II II
r-I r-I '>~ 'i^f -.-I II \ \
a~rrrrr :3\\\
N aOro C N
>~
a \ U r-I N

E C, C. C C, -0 =H 0 ===1 EA EA .-I
ro i ro ~ ro rn co -,, >~

IH a N
~
U O m . 0) t~ = G
a o =~
o cs4 cc ~I
b11 ~4 ~
CZ ~4 C) o ~
02b ~

O U N

It is apparent from Table I that in some instances the subjective response ratings contradicted 1 the actual VA measurements. This indicates that patients believe that VA measurements do not always provide the optimal lens design for that particular patient. Thus, it is important for the neural network to be trained with both actual VA and subjective response data so that the resulting trained neural network correctly models the nonlinear relationship between the various lens designs and both actual VA
and subjective response ratings.
The next step was to utilize @BRAIN NNDS to create a fully connected feedforward neural network, as shown in Figure 1, using the identified optical design and patient parameters as inputs, and actual VA
and subjective response as outputs. Accordingly, the input layer has seven inputs (I1-I7) and the output layer has five processing elements (PE01-PE05) including five outputs (01-05). The remaining structure of the neural network consists of a first hidden layer having eleven processing elements (PEAl-PEA11) and a second hidden layer having eight processing elements (PEB1-PEB8).
The neural network is interconnected as shown in Figure 1. All seven inputs of the input layer are routed to the inputs of all eleven processing elements of the first hidden layer. The outputs of all eleven processing elements of the first hidden layer are input to all eight processing elements of the second hidden layer. The outputs of all eight processing elements of the second hidden layer are input to all five processing elements of the output layer. The outputs of the five processing elements of the output layer are the five outputs of the neural network which correspond to actual VA and subjective response.
After the neural network was created, the neural network was trained using the data collected during the clinical evaluations. Training was accomplished by inputting the spreadsheet containing the clinically collected data into the neural network 1 and allowing the training algorithm to learn the nonlinear relationship between the input design data and the output VA/subjective data.
The training process is accomplished as follows. The seven optical design and patient parameters are routed from the input layer to the inputs of all eleven processing elements of the first hidden layer. Each processing element, as shown in Figure 2, multiplies each input by a different connection weight and adds the individual products.
The results are passed through the sigmoid non-linear transfer function, as defined above, to produce eleven processing element outputs. All eleven processing element outputs of the first hidden layer are routed to the inputs of all eight processing elements of the second hidden layer. Again, each processing element multiplies each input by a different connection weight, adds the individual products and passes the results through the sigmoid non-linear transfer function to produce eight processing element outputs.
All eight processing element outputs of the second hidden layer are routed to the inputs of all five processing elements of the output layer. And again, each processing element multiplies each input by a different connection weight, adds the individual products and passes the results through the sigmoid non-linear transfer function to produce five processing element outputs. These five processing element outputs are the five predicted VA/subjective outputs of the neural network.
During the training process, the training algorithm dynamically compares the five predicted VA/subjective outputs with the associated training data VA/subjective outputs for each set of training data inputs and dynamically changes the connection weights to find a set which minimizes the error between the predicted outputs and the training data outputs. Thus, after all the training data has been processed by the neural network, the resulting set of connection weights should accurately model the relationship between all the training data inputs and the associated training data outputs. The resulting neural network is called a trained neural network.
5 After training, the trained neural network was tested using a random ten percent of the clinically collected data. This data is defined as test data. Test data should not consist of data that is used to train the neural network. This insures, 10 with successful test results, that the neural network has learned, not memorized, the relationship between the input training data and the output training data.
In the present embodiment, test data was distinguished from training data by utilizing a three 15 digit random number generator to assign a three digit number, from 0.000 to 0.999, to each design case. The three digit numbers were generated and integrated into the spreadsheet as random training/learning numbers (see Table I). The neural network was programed to use 20 any data associated with a number less than or equal to 0.100 as test data and any data associated with a number greater than 0.100 as training data.
Testing was accomplished by inputting test inputs, consisting of seven optical design and patient 25 parameters, into the input layer of the trained neural network. In the input layer, the seven test inputs are .routed from the input layer to the inputs of all eleven processing elements of the first hidden layer.
Each processing element multiplies each input by the 30 different connection weights that were determined during training and adds the individual products. The results are passed through the sigmoid non-linear transfer function to produce eleven processing element outputs. All eleven processing element outputs of the 35 first hidden layer are routed to the inputs of all eight processing elements of the second hidden layer.
Again, each processing element multiplies each input by the predetermined connection weights, adds the individual products and passes the results through the 1 sigmoid non-linear transfer function to produce eight processing element outputs. All eight processing element outputs of the second hidden layer are routed to the inputs of all five processing elements of the output layer. And again, each processing element multiplies each input by the predetermined connection weights, adds the individual products and passes the results through the sigmoid non-linear transfer function to produce five processing element outputs.
These five processing element outputs produce five predicted test outputs consisting of predicted VA and subjective test outputs.
The predicted test outputs were compared to the clinically collected VA/subjective outputs to determine if they were substantially the same. The predicted outputs were substantially the same for each patient in the test group, thus, the trained neural network was deemed valid.
A further sanity test was performed as follows. The clinically collected data was grouped by design code and averaged, thereby obtaining average values for all inputs and actual patient performance outputs of each design group (see Table II).

.., U
N v N N
U N=
47 N a C a~ rl Z U] O
rl i 1(1 i O C r+1 i cM
N
W
Q >~
a ~n r ~ o r ~n w aJ in o o ~n ~o w m m o m ro a~ m m ~o r a~ v .~ ~ c~ w H ,~
Z a ~ o o v W ~
L
a U

0 77 a f'1 \0 m Ot V m 01 00 H
,'7. C!] nl f'1 N m N N m N N (+1 ~
v U ~
-.+ C in w 0) r w H w w rn ~
'O SJ C 01 tD n1 V~ r=1 N lf1 47 f"1 t+1 v ro s~ y Q o 0 0 0 0 0 0 0 0 0 a z>

v iJ
a ~
y Q U
w o ., ro a y~ y Gl a A1 r H r 01 V7 0) l0 [O r ro U Cm m r1 N t+1 N N N N N N
a~
v ~ ^
m ~
~ o w r Ln w r m o .i , H N
O ~O N r ro ~ r r l0 ~f O
w ro ~ ; C~ o 0 0 0 0 0 0 0 0 0 ~ R m y 0 W $4 G) m N N N N N N N N N N ~
n7 ~ -.i W 1!1 U) In If1 l!1 1(1 I71 If1 Ifl I11 ~ aJ 0 rt r1 1'1 t~ r~ t~1 rn M ~n ~n m ~n at i i ~ G
y N
Q a N
x U
~ g -.0+ r r r r r r r r r r r ca s z oam w .-~ -u a q ro a~ w La Ze +i N N N N N N N N N N
L v \
LI N M1i QI ri ri r-i e-1 r-1 rl rl rl rl '-I L
a L a ro In l!1 Ifl Itl Ifl V1 tfl tfl t!1 lf1 m I

O u N'n 0 U
ro a =ri O! II \\ ~
cr a a~
a \\\ -~
~ H a ~ b ro a ~i N .-I N N H N m H C C) H tJl ro ro \ v 14 Ol Gl OO uN C.' 11 ~ E.q E A N ~ N N N N H H H H 0 M
A1C 3 in in in in in in in vi 'i ro .% m a N N a ~
N m ~ 0) Ln Ln ~
a W -, N N .C
q 0 ~+ ~ r ~o ~o r+i c ~c ~c ~ r q N

=a' a~ r~ w ~n r 0) ,-i rn c ~n ~ z m ro o 0 0 0 0 H ,1 0 0 41 0 ri rl H H H rl rl 11 N N
~O U > > > > > > > > > >

The average inputs for each individual design code were used as inputs to the tr-ained neural network 1 to obtain predicted near VA and surjective response outputs. These predicted outputs are listed in TABLE
II.
The predicted outputs were compared to the actual average outputs and the dif=erence between the VA and subjective outputs was calc~lated and recorded in TABLE II. The VA difference waS calculated in Letters gained or lost, which is calculated by multiplying the VA difference in V= units by eight, which is the number of Letters on a line. The subjective difference was based on a 0 to 50 unitless range.
As is shown in Table II, zhe VA difference ranged from -4.80 to 2.79 Letters and the subjective difference ranged from -4 to 5, bo-h of which are negligible. Thus, the sanity test further proved that the trained neural network was valid.
After training and testing, the trained neural network was used to predict distance and near VA and subjective response for several lens design as a function of the number of rings cf each lens. To elaborate, average input parameters for five lens designs were input into the trained neural network as constants, except for number of ri~gs which was varied from 1 to 12. Thus, a set of pred_cted VA and subjective outputs were obtained t:at were a function of number of rings.
The predicted VA and subjective output data was then plotted graphically for t:e following five lens designs: monocular d/q/d, monccular d/q, monocular n/q/d, binocular d/q, and binocular d/q/d.
Four graphs were created ;see Figures 3-6) for each of the following outputs:
(a) predicted actual distance VA at high luminance high contrast plotted against number of rings for the five lens designs defined above (Figure 3) ;

(b) predicted subjective distance response plotted against number of rings for the five lens designs defined above (Figure 4);
(c) predicted actual near VA at high luminance high contrast plotted against number of rings for the five lens designs defined above (Figure 5); and (d) predicted subjective near response plotted against number of rings for the five lens designs defined above (Figure 6).
It was then noted that in each case the subjective acuity results differed from the actual acuity results, for both distance vision and near vision. In each case, the actual acuity results indicated that a small number of rings would provide the best results, depicted as a negative slope, while the subjective results indicated a clear preference a larger number of rings, which was depicted as a positive slope. To determine the optimal number of rings, or peak value that would harmonize the conflicting results, each of the two sets of data was normalized and integrated as a composite graph, and these two composites integrations were integrated into a single composite which depicted a single optimal value for all of the integrated data sets.
This composite integration may be expressed in the following formula:

Qual(n) = Sub(n) x Obj(n) where Qual is quality merit function or composite integration;
Subj is subjective response;
Obj is objectively measured response (VA);
and n is the number of rings, or the clinical variation under evaluation.

To create the first composite integration, the predicted actual distance VA and subjective distance response were normalized to a value of 1Ø
The normalized distance VA values that corresponded to the same number of rings were multiplied by the subjective distance values and plotted against number 5 of rings for the five lens designs defined above. The resulting graph was a distance normalized composite as shown in Figure 7.
The same integration was then preformed on the near vision data, with the predicted actual near 10 VA and subjective near response normalized to a value of 1Ø The normalized near VA values that corresponded to the same number of rings were multiplied by the subjective near values and plotted against number of rings for the five lens designs 15 defined above. The resulting graph was a near normalized composite as shown in Figure 8.
The distance normalized composite was then multiplied by the near normalized composite and the result was plotted against number of rings for the 20 five lens designs defined above. The resulting graph was a distance/near normalized composite as shown in Figure 9.
The final integration was prepared multiplying each of the resulting distance/near 25 normalized values for each of the five lens designs and plotting them against the number of rings. The resulting graph was a distance/near normalized composite for all five lens types combined, as shown in Figure 10. This graph represents a single function 30 which defines the combination of all the predicted outputs for all five lens designs. This graph peaks between 5 and 7 rings. Thus, the optimal number of rings for a large population of patients having presbyopia and requiring medium refractive add 35 correction is between 5 and 7.
This composite integration technique could be used for harmonizing or peaking disparate data sets to select any single design element that is separately tested. It could also be used if zhe subjective and 1 actual results for a particular diagnostic remedy were in general harmony, as was the case with the integration of the results of Figure 9 to the single optimzation composite illustrated In Figure 10.
The preferred embodiment of the present inventive method described herein can be modified to determine the optimal number of rings for patients having other optical refractive cc-ditions including, but not limited to, astigmatism ar-d cataracts.
The preferred embodiment of the present inventive method described herein also can be modified to determine optimal values for other design parameters including, but not limi--ed to, monocular or binocular pair and pupil function.
While the present inventive method has been described herein with respect to a preferred embodiment, it will be understood by those skilled in the art that the foregoing and other changes in form and details may be made therein without departing from the spirit and scope of the invention.

Claims (40)

1. A method of optimizing optical designs involving a plurality of design variables, said method comprising:
(a) identifying predetermined optical design parameters relevant to a predetermined optical refractive condition;
(b) forming optical lenses utilizing one or more of said optical design parameters for use in clinical evaluations each such evaluation providing visual acuity data and subjective response ratings for a defined number of patients having said optical refractive condition;
(c) inputting said optical design parameters and related patient parameters as input components, and said visual acuity data and said subjective response ratings as output components into a neural network;
(d) training said neural network to model significant relationships between said input components and said output components to thereby produce a trained neural network;
(e) isolating and inputting one or more specific design parameters for evaluation by said trained neural network to predict visual acuity and subjective response as a function of said specific design parameter; and (f) integrating one or more of said predictions to determine the optimal optical design to correct said optical refractive condition.
2. The method according to claim 1, wherein said optical refractive condition is selected from a group consisting of presbyopia and astigmatism.
3. The method according to claim 2, wherein said optical refractive condition is presbyopia.
4. The method according to claim 3, wherein said optical lenses are formed as multifocal contact lenses having a plurality of concentric annular rings with alternating distance and near add power rings.
5. The method according to claim 4, wherein said multifocal contact lenses are formed as a monocular pair or a binocular pair wherein a monocular pair comprising identical lenses, and a binocular pair comprising slightly different lenses.
6. The method according to claim 5, wherein said optical design parameters are selected from a group consisting of:
(a) the number of concentric annular rings;
(b) a pupil function wherein said alternating rings are grouped into at least two regions having add power distributions defined by a pupil function type;
(c) a lens add power; and (d) a monocular or binocular pair.
7. The method according to claim 6, wherein said optical design parameters include said number of concentric annular rings and said pupil functions.
8. The method according to claim 6, wherein said optical design parameters include said number of concentric annular rings, said pupil functions, said lens add power and either said monocular or said binocular pair.
9. The method according to claim 5, wherein said optical refractive condition is presbyopia, and said predicted visual acuity and said predicted subjective response is selected from a group consisting of:

(a) predicted distance visual acuity;
(b) predicted near visual acuity;
(c) predicted subjective near response;
(d) predicted subjective distance response; and (e) predicted subjective overall response.
10. The method according to claim 9, wherein said predicted visual acuity include said predicted distance visual acuity and said predicted near visual acuity.
11. The method according to claim 9, wherein said predicted subjective response include said predicted subjective near response and said predicted subjective distance response.
12. The method according to claim 9, wherein said predicted visual acuity and said predicted subjective response include said predicted distance and near visual acuity and said predicted subjective distance and near response.
13. The method according to claim 12, wherein said predicted visual acuity and said predicted subjective response are normalized and combined to determine an optimum optical design for a large population of patients having said optical refractive condition.
14. The method according to claim 5, further comprising the steps of:

grouping said input components according to said optical design parameters thereby producing a plurality of design groups; and averaging said input components for each of said design groups thereby obtaining average input component values.
15. The method according to claim 14, wherein said specific design parameter is the number of concentric annular rings; said number of rings being varied and input, with said average input component values, into said trained neural network for evaluation by said trained neural network; said trained neural network outputting said predicted visual acuity and said predicted subjective response as a function of said number of concentric annular rings.
16. The method according to claim 15, wherein said predicted visual acuity and said predicted subjective response is normalized and combined to determine an optimum optical design for a large population of patients having said optical refractive condition.
17. The method according to claim 5, wherein the neural network used in said inputting, training and isolating steps has an input layer having at least one input for receiving said input components and at least one hidden layer having at least one first processing element for modeling complex functions.
18. The method according to claim 17, wherein the neural network further includes an output layer having at least one second processing element, said output layer having one output for receiving said output components and for outputting said predicted visual acuity and said predicted subjective response.
19. The method according to claim 18, wherein said second processing element may be either said first processing element or a linear processing element.
20. The method according to claim 19, wherein said neural network includes said input layer, said at least one hidden layer including a first hidden layer and a second hidden layer, and said output layer.
21. The method according to claim 20, wherein said input layer includes seven inputs, said first hidden layer includes eleven processing elements, said second hidden layer includes eight processing elements, and said output layer includes five processing elements and five outputs.
22. The method according to claim 21, wherein said optical refractive condition is presbyopia and said seven inputs are selected from a group consisting of:
(a) the number of concentric annular rings;
(b) a pupil function;
(c) a lens add power;
(d) monocular or binocular pair;
(e) patient age;
(f) patient refractive add; and (g) Hloss
23. The method according to claim 21, wherein said optical refractive condition is presbyopia, and said five outputs are selected from a group consisting of:
(a) actual distance visual acuity;

(b) actual near visual acuity;
(c) subjective near response;
(d) subjective distance response; and (e) subjective overall response.
24. The method according to claim 17, wherein said neural network utilizes Talon's training algorithm which defines each of said first processing elements as a sigmoid non-linear transfer function defined by an equation as follows:

out = 1 . 0/(1+e in) -0.5 where out is a processing element output, in is a processing element input and e-in is 1 × 10-in.
25. The method according to claim 17, wherein said neural network is a fully connected feedforward design.
26. The method according to claim 5, wherein said optical refractive condition is presbyopia and said defined number of patients are selected from a group consisting of:
(a) low refractive add patients;
(b) medium refractive add patients; and (c) high refractive add patients.
27. The method according to claim 26, wherein said defined number of patients are said medium refractive add patients.
28. The method according to claim 5, wherein said optical refractive condition is presbyopia, and said one or more specific design parameter is selected from a group consisting of:
(a) the number of concentric rings;
(b) a pupil function;
(c) a lens add power; and (d) a monocular or binocular pair.
29. The method according to claim 1, wherein said patient parameters are selected from a group consisting of:
(a) patient age;
(b) patient refractive add; and (c) Hloss.
30. The method according to claim 29, wherein said patient parameters include said patient age, said patient refractive add, and said Hloss.
31. The method according to claim 1, wherein said optical refractive condition is presbyopia.
32. The method according to claim 31, wherein said visual acuity data is obtained during a clinical evaluation which includes distance and near visual acuity at high luminance and high contrast measured in units of -10 LOG MAR, where MAR is minimum angle of resolution.
33. The method according to claim 32, wherein said visual acuity data is obtained from a clinical evaluation which includes distance and near visual acuity at high luminance and high contrast.
34. The method according to claims 32 or 33, wherein said visual acuity data is measured in lines lost from a patient's best spectacle correction.
35. The method according to claim 31, wherein said subjective response rating are obtained from a clinical evaluation which includes distance response, near response and overall response rated in unitless values from a predetermined range.
36. The method according to claim 35, wherein said subjective response ratings includes distance, near and overall subjective response ratings.
37. The method according to claim 1 which further comprises the step of testing said neural network with test data to verify that said trained neural network is valid.
38. The method according to claim 37, wherein said test data is a random 10% of said input data.
39. The method according to claim 38, wherein said specific design parameter is the number of concentric rings.
40. The method according to claim 39, wherein the number of concentric rings is varied and input into said trained neural network for evaluation, to thereby predict said visual acuity and said subject response as a function of the number of concentric rings.
CA002204692A 1996-05-09 1997-05-07 Neural network analysis for multifocal contact lens design Expired - Fee Related CA2204692C (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/647,005 US5724258A (en) 1996-05-09 1996-05-09 Neural network analysis for multifocal contact lens design
US08/647,005 1996-05-09

Publications (2)

Publication Number Publication Date
CA2204692A1 CA2204692A1 (en) 1997-11-09
CA2204692C true CA2204692C (en) 2009-08-04

Family

ID=24595329

Family Applications (1)

Application Number Title Priority Date Filing Date
CA002204692A Expired - Fee Related CA2204692C (en) 1996-05-09 1997-05-07 Neural network analysis for multifocal contact lens design

Country Status (8)

Country Link
US (1) US5724258A (en)
EP (1) EP0806694B1 (en)
JP (1) JP4018193B2 (en)
AT (1) ATE227854T1 (en)
AU (1) AU712104B2 (en)
CA (1) CA2204692C (en)
DE (1) DE69716994T2 (en)
SG (1) SG80571A1 (en)

Families Citing this family (55)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ATE336196T1 (en) * 1998-03-04 2006-09-15 Visx Inc SYSTEM FOR LASER TREATMENT OF PRESALES VISION
US6176580B1 (en) * 1999-04-02 2001-01-23 Johnson & Johnson Vision Care, Inc. Method of designing and fitting contact lenses taking into account material properties of the lenses
CA2798508A1 (en) * 2000-10-18 2002-04-25 Johnson & Johnson Consumer Companies, Inc. Intelligent performance-based product recommendation system
EP1205863A1 (en) * 2000-11-14 2002-05-15 Honda R&D Europe (Deutschland) GmbH Multi-objective optimization
AU2002252769A1 (en) 2000-11-30 2002-06-11 Sarver And Associates Advanced vision intervention algorithm
EP1235180A1 (en) 2001-02-26 2002-08-28 Honda R&D Europe (Deutschland) GmbH Parameter adaptation in evolution strategies
US7605649B2 (en) * 2001-03-13 2009-10-20 Marvell World Trade Ltd. Nested transimpedance amplifier
JP2004534964A (en) 2001-04-27 2004-11-18 ノバルティス アクチエンゲゼルシャフト Automatic lens design and manufacturing system
US7123985B2 (en) * 2002-04-12 2006-10-17 Johnson & Johnson Vision Care, Inc. Design build test cycle reduction
US6932808B2 (en) * 2002-11-19 2005-08-23 Visx, Incorporated Ablation shape for the correction of presbyopia
DE10255189B4 (en) * 2002-11-27 2005-12-22 Photeon Technologies Gmbh Method of designing waveguide geometries in integrated optical devices
US7896916B2 (en) 2002-11-29 2011-03-01 Amo Groningen B.V. Multifocal ophthalmic lens
SE0203564D0 (en) 2002-11-29 2002-11-29 Pharmacia Groningen Bv Multifocal opthalmic lens
US8342686B2 (en) 2002-12-06 2013-01-01 Amo Manufacturing Usa, Llc. Compound modulation transfer function for laser surgery and other optical applications
US8911086B2 (en) 2002-12-06 2014-12-16 Amo Manufacturing Usa, Llc Compound modulation transfer function for laser surgery and other optical applications
US7460288B2 (en) * 2002-12-06 2008-12-02 Amo Manufacturing Usa, Llc Methods for determining refractive corrections from wavefront measurements
EP1567907A4 (en) * 2002-12-06 2009-09-02 Amo Mfg Usa Llc Presbyopia correction using patient data
US7434936B2 (en) * 2002-12-06 2008-10-14 Amo Manufacturing Usa, Llc Residual accommodation threshold for correction of presbyopia and other presbyopia correction using patient data
US7320517B2 (en) 2002-12-06 2008-01-22 Visx, Incorporated Compound modulation transfer function for laser surgery and other optical applications
US6951391B2 (en) * 2003-06-16 2005-10-04 Apollo Optical Systems Llc Bifocal multiorder diffractive lenses for vision correction
AU2004292165B2 (en) 2003-11-14 2010-09-09 Essilor International (Compagnie Generale D'optique) System for manufacturing an optical lens
EP1598751B1 (en) * 2004-01-12 2014-06-25 Honda Research Institute Europe GmbH Estimation of distribution algorithm (EDA)
EP1557788B1 (en) * 2004-01-26 2008-04-16 Honda Research Institute Europe GmbH Reduction of fitness evaluations using clustering technique and neural network ensembles
JP4897497B2 (en) * 2004-02-20 2012-03-14 ヴィズイクス・インコーポレーテッド Volume point spread function for eye diagnosis and treatment
US20050200809A1 (en) * 2004-02-20 2005-09-15 Dreher Andreas W. System and method for analyzing wavefront aberrations
US20050261752A1 (en) * 2004-05-18 2005-11-24 Visx, Incorporated Binocular optical treatment for presbyopia
US7387387B2 (en) * 2004-06-17 2008-06-17 Amo Manufacturing Usa, Llc Correction of presbyopia using adaptive optics and associated methods
US7156516B2 (en) * 2004-08-20 2007-01-02 Apollo Optical Systems Llc Diffractive lenses for vision correction
US7025456B2 (en) * 2004-08-20 2006-04-11 Apollo Optical Systems, Llc Diffractive lenses for vision correction
US7922326B2 (en) 2005-10-25 2011-04-12 Abbott Medical Optics Inc. Ophthalmic lens with multiple phase plates
EP3480650A1 (en) * 2004-10-25 2019-05-08 Johnson & Johnson Surgical Vision, Inc. Ophthalmic lens with multiple phase plates
US20060229932A1 (en) * 2005-04-06 2006-10-12 Johnson & Johnson Services, Inc. Intelligent sales and marketing recommendation system
US7413566B2 (en) * 2005-05-19 2008-08-19 Amo Manufacturing Usa, Llc Training enhanced pseudo accommodation methods, systems and devices for mitigation of presbyopia
EP1768053A1 (en) * 2005-09-12 2007-03-28 Honda Research Institute Europe GmbH Evolutionary search for robust solutions
US8147063B2 (en) * 2006-04-21 2012-04-03 Otc Optics Llc Method for minimizing prism in over-the-counter eyeglasses and optical devices
SG148891A1 (en) * 2007-06-21 2009-01-29 Novartis Ag Engineering expert system
US7625086B2 (en) * 2007-08-28 2009-12-01 Johnson & Johnson Vision Care, Inc. Method of designing multifocal contact lenses
CA3123266A1 (en) 2012-08-31 2014-03-06 Amo Groningen B.V. Multi-ring lens, systems and methods for extended depth of focus
US9304332B2 (en) * 2013-08-22 2016-04-05 Bespoke, Inc. Method and system to create custom, user-specific eyewear
AU2017218680B2 (en) 2016-02-09 2021-09-23 Amo Groningen B.V. Progressive power intraocular lens, and methods of use and manufacture
US10018854B2 (en) 2016-06-22 2018-07-10 Indizen Optical Technologies of America, LLC Custom ophthalmic lens design derived from multiple data sources
EP3321831B1 (en) 2016-11-14 2019-06-26 Carl Zeiss Vision International GmbH Device for determining predicted subjective refraction data or predicted subjective correction data and computer program
AU2018235011A1 (en) 2017-03-17 2019-10-24 Amo Groningen B.V. Diffractive intraocular lenses for extended range of vision
US11523897B2 (en) 2017-06-23 2022-12-13 Amo Groningen B.V. Intraocular lenses for presbyopia treatment
CA3067116A1 (en) 2017-06-28 2019-01-03 Amo Groningen B.V. Diffractive lenses and related intraocular lenses for presbyopia treatment
CA3068351A1 (en) 2017-06-28 2019-01-03 Amo Groningen B.V. Extended range and related intraocular lenses for presbyopia treatment
US11327210B2 (en) 2017-06-30 2022-05-10 Amo Groningen B.V. Non-repeating echelettes and related intraocular lenses for presbyopia treatment
CN107976804B (en) * 2018-01-24 2020-11-24 苏州浪潮智能科技有限公司 Design method, device, equipment and storage medium of lens optical system
CN113056699B (en) * 2018-11-15 2023-07-21 依视路国际公司 Method and system for determining parameters for manufacturing an optical article and corresponding optical article
WO2020102762A1 (en) * 2018-11-16 2020-05-22 Arizona Board Of Regents On Behalf Of The University Of Arizona Optical system design
WO2021040078A1 (en) * 2019-08-27 2021-03-04 (주)비쥬웍스 Lens determination method and device using same
EP3798944A1 (en) * 2019-09-30 2021-03-31 Hoya Lens Thailand Ltd. Learning model generation method, computer program, eyeglass lens selection support method, and eyeglass lens selection support system
EP4085292A1 (en) 2019-12-30 2022-11-09 AMO Groningen B.V. Lenses having diffractive profiles with irregular width for vision treatment
DE102020004840A1 (en) * 2020-08-07 2022-02-10 Rodenstock Gmbh Improved calculation of ophthalmic lenses
WO2023135745A1 (en) * 2022-01-14 2023-07-20 オリンパス株式会社 Optical system design system, optical system design method, trained model, program, and information recording medium

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2573876A1 (en) * 1984-11-26 1986-05-30 Vinzia Francis MULTIFOCAL LENS, PROCESS FOR PREPARING SAME, AND USE AS A CONTACT LENS OR AN INTRAOCULAR IMPLANT TO REPLACE THE CRYSTALLINE
FI79619C (en) * 1984-12-31 1990-01-10 Antti Vannas Intraocular lens
EP0201231A3 (en) * 1985-05-03 1989-07-12 THE COOPER COMPANIES, INC. (formerly called CooperVision, Inc.) Method of treating presbyopia with concentric bifocal contact lenses
EP0227653A1 (en) * 1985-06-24 1987-07-08 BRONSTEIN, Leonard Contact lens
WO1987007496A1 (en) * 1986-06-05 1987-12-17 Precision-Cosmet Co., Inc. One-piece bifocal intraocular lens construction
US5225858A (en) * 1987-06-01 1993-07-06 Valdemar Portney Multifocal ophthalmic lens
US4838675A (en) * 1987-06-19 1989-06-13 Sola International Holdings, Ltd. Method for improving progressive lens designs and resulting article
US4869587A (en) * 1987-12-16 1989-09-26 Breger Joseph L Presbyopic contact lens
ES2081811T3 (en) * 1988-07-20 1996-03-16 Allen L Dr Cohen MULTIFOCAL DIFRACTIVE OPTICAL ELEMENT.
GB9008582D0 (en) * 1990-04-17 1990-06-13 Pilkington Diffractive Lenses Method and contact lenses for treating presbyobia
US5112351A (en) * 1990-10-12 1992-05-12 Ioptex Research Inc. Multifocal intraocular lenses
US5198844A (en) * 1991-07-10 1993-03-30 Johnson & Johnson Vision Products, Inc. Segmented multifocal contact lens
FR2688898A1 (en) * 1992-03-23 1993-09-24 Chemoul Alain Contact lens and method of producing it
US5430506A (en) * 1992-11-06 1995-07-04 Volk; Donald A. Indirect ophthalmoscopy lens for use with slit lamp biomicroscope
US5448312A (en) * 1992-12-09 1995-09-05 Johnson & Johnson Vision Products, Inc. Pupil-tuned multifocal ophthalmic lens
US5619289A (en) * 1993-03-31 1997-04-08 Permeable Technologies, Inc. Multifocal contact lens
US5526071A (en) * 1993-03-31 1996-06-11 Permeable Technologies, Inc. Multifocal contact lens and method for preparing
US5493350A (en) * 1993-03-31 1996-02-20 Seidner; Leonard Multipocal contact lens and method for preparing
US5404183A (en) * 1993-03-31 1995-04-04 Seidner; Leonard Multifocal contact lens and method for preparing
US5517260A (en) * 1994-03-28 1996-05-14 Vari-Site, Inc. Ophthalmic lens having a progressive multifocal zone and method of manufacturing same
AU7955794A (en) * 1994-09-16 1996-03-29 Permeable Technologies Inc. Multifocal contact lens and method for preparing
IL117935A0 (en) * 1995-05-04 1996-08-04 Johnson & Johnson Vision Prod Multifocal ophthalmic lens

Also Published As

Publication number Publication date
AU2006297A (en) 1997-11-27
DE69716994T2 (en) 2003-07-03
EP0806694A3 (en) 1999-02-24
SG80571A1 (en) 2001-05-22
US5724258A (en) 1998-03-03
EP0806694B1 (en) 2002-11-13
JPH1068913A (en) 1998-03-10
EP0806694A2 (en) 1997-11-12
MX9703455A (en) 1998-06-28
CA2204692A1 (en) 1997-11-09
JP4018193B2 (en) 2007-12-05
AU712104B2 (en) 1999-10-28
ATE227854T1 (en) 2002-11-15
DE69716994D1 (en) 2002-12-19

Similar Documents

Publication Publication Date Title
CA2204692C (en) Neural network analysis for multifocal contact lens design
RU2464604C2 (en) Method of making multifocal contact lenses
CN101686802B (en) Apparatus and method for determining necessary correction of defective vision of eye
US7954950B2 (en) System and method for analyzing wavefront aberrations
CN107358036A (en) A kind of child myopia Risk Forecast Method, apparatus and system
US20110270596A1 (en) Apparatus, System and Method for Predictive Modeling to Design, Evaluate and Optimize Ophthalmic Lenses
CN110235051B (en) Contact lens
MX2007005849A (en) Correction of higher order aberrations in intraocular lenses.
CN1849091A (en) Ocular aberration correction taking into account fluctuations dueto biophysical rhythms
US20110157547A1 (en) Method of designing progressive addition lenses
CN102460273A (en) Making of progressive spectacle lens customized based on blur perception
CN112313564A (en) Method for determining an ophthalmic device and associated system
KR20110042209A (en) Fitting method for multifocal lenses
CN113039479A (en) Method and apparatus for evaluating the efficacy of an ophthalmic lens in controlling vision disorders
KR101645953B1 (en) Lens design simplification process
MXPA97003455A (en) Neural network analysis for the design of multifour contact lenses
CN105319736A (en) Patient interactive fit tool and methodology for contact lens fitting
WO2023187089A1 (en) Device for determining a shift in a refraction value of an eye
KR20230002799A (en) Multi-lens system for presbyopia

Legal Events

Date Code Title Description
EEER Examination request
MKLA Lapsed

Effective date: 20170510