CN113855387B - Indirect compensation method for full-higher-order aberration in cornea refraction correction - Google Patents
Indirect compensation method for full-higher-order aberration in cornea refraction correction Download PDFInfo
- Publication number
- CN113855387B CN113855387B CN202111245102.6A CN202111245102A CN113855387B CN 113855387 B CN113855387 B CN 113855387B CN 202111245102 A CN202111245102 A CN 202111245102A CN 113855387 B CN113855387 B CN 113855387B
- Authority
- CN
- China
- Prior art keywords
- astigmatism
- cornea
- lens
- corn
- axial position
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 210000004087 cornea Anatomy 0.000 title claims abstract description 103
- 230000004075 alteration Effects 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title claims abstract description 43
- 201000009310 astigmatism Diseases 0.000 claims abstract description 141
- 210000000695 crystalline len Anatomy 0.000 claims abstract description 20
- 239000013078 crystal Substances 0.000 claims abstract description 13
- 201000009308 regular astigmatism Diseases 0.000 claims abstract description 11
- 201000000766 irregular astigmatism Diseases 0.000 claims abstract description 7
- 238000012876 topography Methods 0.000 claims description 13
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 abstract description 3
- 230000004438 eyesight Effects 0.000 description 9
- 238000007689 inspection Methods 0.000 description 3
- 230000019771 cognition Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 210000001525 retina Anatomy 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 206010010071 Coma Diseases 0.000 description 1
- 206010020675 Hypermetropia Diseases 0.000 description 1
- 238000002679 ablation Methods 0.000 description 1
- 230000004308 accommodation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013479 data entry Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 201000006318 hyperopia Diseases 0.000 description 1
- 230000004305 hyperopia Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000004377 improving vision Effects 0.000 description 1
- 230000000860 keratorefractive effect Effects 0.000 description 1
- 208000001491 myopia Diseases 0.000 description 1
- 230000004379 myopia Effects 0.000 description 1
- 230000002980 postoperative effect Effects 0.000 description 1
- 208000014733 refractive error Diseases 0.000 description 1
Classifications
-
- 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
- A61F9/013—Instruments for compensation of ocular refraction ; Instruments for use in cornea removal, for reshaping or performing incisions in 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/103—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for determining refraction, e.g. refractometers, skiascopes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/103—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for determining refraction, e.g. refractometers, skiascopes
- A61B3/1035—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for determining refraction, e.g. refractometers, skiascopes for measuring astigmatism
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/107—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for determining the shape or measuring the curvature of the cornea
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- 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
- A61F9/008—Methods or devices for eye surgery using laser
- A61F2009/00861—Methods or devices for eye surgery using laser adapted for treatment at a particular location
- A61F2009/00872—Cornea
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Ophthalmology & Optometry (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Vascular Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Eye Examination Apparatus (AREA)
Abstract
The application provides an indirect compensation method of full-higher-order aberration in cornea refraction correction, and belongs to the field of cornea correction. The indirect compensation method comprises the steps of firstly obtaining a subjective optometry sphere, a cylindrical lens and an axial position, and then obtaining the axial position of the flat-axis curvature, the steep-axis curvature and the flat-axis curvature of the cornea to be corrected; then, respectively obtaining the astigmatism and the axial position of the back surface of the cornea to be corrected and the regular astigmatism and the axial position of the front surface of the cornea; the cornea front surface truly diffuses light and axial position, and the residual astigmatism proportion of the natural crystal is selected; then calculating total astigmatism and axial position of cornea, astigmatism and axial position of crystalline lens, calculating all astigmatism and axial position of cornea except irregular astigmatism, then calculating target diopter, and finally calculating target astigmatism, axial position and spherical lens for dioptric correction. The application realizes the unified compensation of the full higher order aberration influence quantity, reduces input items, simplifies calculation, optimizes operation flow and reduces the introduction of human errors; meanwhile, the compensation precision is improved, and better visual quality is obtained.
Description
Technical Field
The application belongs to the field of vision imaging and cornea correction, and particularly relates to an indirect compensation (ZZ All-auto-compensation Refraction, ZZAR) method for full-higher-order aberration in cornea refraction correction.
Background
When the human eye is imaged, aberrations are formed because the human eye is imperfect. Aberrations formed by the human eye include lower-order aberrations and higher-order aberrations, wherein the lower-order aberrations generally refer to aberrations caused by refractive errors of the eye such as myopia, hyperopia, regular astigmatism and the like, and the higher-order aberrations generally refer to aberrations caused by defects of the eye itself, and include three-order coma, clover-like astigmatism, four-order spherical aberration, other higher-order aberrations and the like. To obtain a better visual effect, the cornea can be corrected.
Because of the different cornea conditions of each individual, a personalized guiding mode is generally adopted for cornea refraction correction, such as wavefront aberration guiding, cornea topography guiding, light ray tracing and the like. Currently, the higher order aberrations of the vision system are generally corrected and/or the increase in post-operative higher order aberrations is reduced by an optimized ablation pattern to increase the resolution and contrast sensitivity of the patient's retina, thereby improving vision quality. However, when the higher order aberration is corrected for the cornea refraction in the personalized guiding mode, the cornea morphology is changed, the lower order aberration of the whole eye is also caused to be changed, and the influence of the lower order aberration on the resolution and the contrast sensitivity of the retina is far greater than that of the higher order aberration. Therefore, when the cornea refraction guided by the individualization corrects the full higher order aberration, if the lower order aberration cannot be well compensated, the correction effect is affected.
In the prior art, when the cornea refraction of a personalized guiding mode corrects higher order aberration, a Zernike polynomial expression is generally used, and fig. 1 is a 3D schematic diagram of a 6 th order 27 th order Zernike polynomial expression of higher order aberration in the prior art. The existing high-order aberration compensation methods include the Phorcides method developed in the United states and the ZZ Vector-compensation Refraction (ZZ VR) method developed in China. However, the above-described compensation scheme for the higher order aberration has the following problems:
firstly, in clinical results, although the correction accuracy of the cylindrical lens of the ZZ VR method is obviously better than that of the Phorcides method, in theory, the ZZ VR method only compensates 4 items with highest weight in the influence of low-order aberration, and the Phorcides method only compensates for the area with high weight according to graphic analysis software, namely, the method of directly calculating the compensation is difficult to exhaust all influence factors and add one compensation;
secondly, the qualitative relation between different Zernike polynomials and low-order aberration is less internationally recognized, and the method can be applied to clinical quantitative relation compensation more individually;
third, the clinical operation of the existing refraction correction full-higher order aberration compensation method needs to be optimized, compared with the Phorcides method, the ZZ VR method has the advantages that the calculation efficiency is remarkably improved, the target result can be obtained more quickly, and the data entry item still needs to be further optimized.
Disclosure of Invention
In order to solve the problem of compensation of all higher-order aberration in cornea refractive correction in the prior art, the embodiment of the application provides an indirect compensation method ZZ AR method of all higher-order aberration in cornea refractive correction, which compensates all higher-order aberration affecting vision and vision quality, improves vision quality after cornea refractive correction, and optimizes operation flow.
In order to achieve the above purpose, the embodiment of the present application adopts the following technical scheme:
the embodiment of the application provides an indirect compensation method of full-higher-order aberration in cornea refraction correction, which comprises the following steps:
s0, obtaining a subjective refraction sphere Mani S, a cylindrical mirror Mani C and an axial position Mani Ax;
step S1, obtaining the flat axis curvature Kf, the steep axis curvature Ks and the flat axis curvature axis KfAx of the cornea to be corrected;
s2, obtaining back surface astigmatism Corn-P C and back surface astigmatism axial position Corn-P Ax of the cornea to be corrected;
step S3, recording cornea front surface regular astigmatism TMRC and regular astigmatism axis TMRAx in cornea topography
S4, selecting a residual astigmatism ratio Res-Lens C of the natural crystal within a range of 0-100%;
s5, calculating the real astigmatism Corn-A C and the real astigmatism axial position Corn-Aax of the front surface of the cornea through Kf, ks, kfAx;
s6, calculating total cornea astigmatism Corn C and total astigmatism axial position Corn Ax according to the real astigmatism and axial position of the front surface of the cornea and the astigmatism and axial position of the rear surface of the cornea;
s7, calculating Lens astigmatism Lens C and axial position LensAx according to the subjective light-detecting cylindrical Lens Mani C, the axial position Mani Ax and total cornea astigmatism Corn C;
s8, calculating all astigmatism C1 and axial position Ax1 except irregular astigmatism of the cornea according to the regular astigmatism TMR C and axial position TMR Ax of the front surface of the cornea, the astigmatism Corn-P C and axial position Corn-P Ax of the rear surface of the cornea, the astigmatism Lens C and axial position Lens Ax;
s9, calculating target diopter according to astigmatism and axial position of the crystalline Lens and the residual astigmatism proportion Res-Lens C of the selected natural crystal;
step S10, calculating refraction correction target astigmatism AR C, axis ARAx and sphere AR S according to all astigmatism C1 except irregular astigmatism and axis Ax1, target diopter and subjective light inspection sphere Mani S.
As a preferred embodiment of the present application, the flat-axis curvature Kf, steep-axis curvature Ks, and flat-axis curvature axis KfAx in the step S1 are obtained from data imported into an excimer apparatus according to a corneal topography.
As a preferred embodiment of the present application, the posterior surface astigmatism Corn-P C and the posterior surface astigmatism axial position Corn-P Ax of the cornea to be corrected are obtained from corneal topography data.
As a preferred embodiment of the present application, the residual astigmatism ratio Res-Lens c=30 to 50% of the natural crystal is selected.
As a preferred embodiment of the application, the back surface astigmatism Corn-P C and the back surface astigmatism axial position Corn-P Ax of the cornea to be corrected are obtained through fixed proportion estimation of the front surface astigmatism of the cornea and the front and back astigmatism of the cornea of a conventional population; when the estimation is performed, corn-P Ax is calculated by the formula (1) according to Kf, ks, kf Ax:
Corn-P Ax＝Kf Ax (1)；
Corn-P C is calculated by equation (2):
as a preferred embodiment of the present application, the front surface true astigmatism is calculated by the formula (3):
Corn-A Ax＝Kf Ax(3)；
the true astigmatism axis is calculated by equation (4):
the total corneal astigmatism axis Corn Ax is calculated by equation (5):
after solving the Corn Ax by the formula (5), converting the LensAx value beyond the range to 0-180 by the Corn Ax + -180 to obtain the final Corn Ax;
and calculating total corneal astigmatism Corn C by formula (6):
the Lens astigmatism axis Lens Ax is calculated by equation (7):
after solving the LensAx by the formula (7), converting the LensAx value beyond the range to 0-180 by the LensAx + -180 to obtain the final LensAx;
and calculate Lens astigmatism Lens C by equation (8):
all astigmatism axis positions Ax1 are calculated by formula (9):
after solving Ax1 through a formula (9), converting an Ax1 value which exceeds the range to 0-180 through Ax 1+/-180 to obtain the final Ax1;
and all astigmatism C1 are calculated by formula (10):
as a preferred embodiment of the present application, the target diopter includes a percentage of residual crystalline astigmatism and an axis.
As a preferred embodiment of the present application, the residual crystalline astigmatism axis is calculated by equation (11):
Target Ax＝Lens Ax (11)：
and the percent residual crystalline astigmatism is calculated by equation (12):
as a preferred embodiment of the present application, the refractive correction target astigmatism axis ARAx is calculated by formula (13):
solving the Lens Ax by the formula (13), and converting the Lens Ax value beyond the range to 0-180 by the Lens Ax + -180 to obtain a final Lens Ax value;
and calculating the refraction correction target astigmatism by the formula (14), and calculating the refraction correction target sphere AR S by the formula (15):
AR S＝Mani S-(AR C-Mani C)/2 (15)。
the ZZ AR method provided by the embodiment of the application has the following beneficial effects:
(1) the theoretical full-higher order aberration term compensation is realized;
(2) the surgical accuracy is obviously improved, and the method has statistical significance;
(3) the operation design is simple, and the operation design can be rapidly completed in a short time through the built-in calculator of the webpage, so that the operation efficiency is greatly improved;
(4) the method is characterized in that the individuation quantitative reservation is made for the unknown variation after the crystal astigmatism morphology operation for the first time.
Of course, it is not necessary for any one product or method of practicing the application to achieve all of the advantages set forth above at the same time.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present application, and other drawings may be obtained according to these drawings without inventive effort for a person of ordinary skill in the art.
FIG. 1 is a 3D schematic of a 6 th order 27 th order Zernike polynomial representing higher order aberrations of the prior art;
fig. 2 is a flowchart of an indirect compensation method for the full higher order aberration in the corneal refractive correction according to an embodiment of the present application.
Detailed Description
The technical problems, technical solutions and advantages of the present application will be elucidated by referring to exemplary embodiments. However, the present application is not limited to the exemplary embodiments disclosed below; this may be implemented in different forms. The essence of the description is merely to aid one skilled in the relevant art in comprehensively understanding the specific details of the application.
Aiming at the compensation problem of the full-higher-order aberration in the refraction correction in the prior art, the embodiment of the application provides an indirect compensation method of the full-higher-order aberration in the refraction correction according to the cognition of Zernike polynomials and ZZ VR and the rule cognition after the cornea topography is converted into a digital model and combining with physical optics and higher mathematical discipline theory, and aims to indirectly and automatically compensate the full-higher-order aberration when the refraction correction is carried out so as to further improve the naked eye vision and vision quality after the refraction correction, optimize the operation flow, and facilitate the simpler improvement of the cornea refraction correction effect according to the relevant standardized operation flow in the refraction correction process.
As shown in fig. 2, the indirect compensation method ZZ AR method for the full higher order aberration in the corneal refractive correction according to the embodiment of the application includes the following steps:
step S0, obtaining a subjective refraction sphere Mani S, a cylindrical mirror Mani C and an axial position Mani Ax.
In this step, the ball mirror Mani S, the cylindrical mirror Mani C and the axial Mani Ax data are obtained by subjective refraction of vision quality required to be subjected to cornea refraction correction. For example, in one specific embodiment, taking the cornea X to be corrected as an example, the result of the refraction is: subjective optometry sphere man s= -4.35, cylinder man c= -1.50, axis man ax=65.
And S1, obtaining a corneal topography of the cornea to be corrected, and obtaining corresponding flat-axis curvature Kf, steep-axis curvature Ks and flat-axis curvature axial position KfAx.
In this step, the flat-axis curvature Kf, steep-axis curvature Ks, and flat-axis curvature axial position kfax are obtained from data introduced into an excimer apparatus according to a corneal topography. For example, taking the cornea X to be corrected as an example, a corneal topography of the cornea to be corrected is introduced into an excimer apparatus, and the following results are obtained: kf=42.88, ks=43.60, kfax=40.
And S2, obtaining the back surface astigmatism Corn-P C and the back surface astigmatism axial position Corn-P Ax of the cornea to be corrected.
In the step, the back surface astigmatism Corn-P C and the back surface astigmatism axial position Corn-P Ax of the cornea to be corrected can be obtained through cornea topographic map data, and can also be obtained through calculation of the real astigmatism of the front surface of the cornea and the fixed proportion of the real astigmatism of the front cornea and the rear cornea of a conventional crowd.
When obtained by estimation, in one calculation mode, the calculation may be performed by Kf, ks, kfAx. Specifically, corn-Pax is calculated by the formula (1):
Com-P Ax＝KfAx (1)：
Corn-P C is calculated by equation (2):
taking cornea X to be corrected as an example, the following results are obtained through corneal topography data: cornPC=0.16, cornPax=24.
Step S3, recording the cornea front surface regular astigmatism tmrc and the regular astigmatism axis TMRAx in the cornea topography.
In this step, taking the cornea X to be corrected as an example, the following results are read from the corneal topography of the cornea X to be corrected: tmrc= -1.16, tmax=18.
And S4, selecting the residual astigmatism ratio Res-Lens C (%) of the natural crystal within the range of 0-100%.
Preferably, the residual astigmatism ratio Res-Lens c=30 to 50% of the natural crystal selected in this step. The selection of this astigmatism ratio is chosen according to the astigmatism parameters of the different cornea.
Taking the cornea X to be corrected as an example, res-Lens c=20% is selected.
Step S5, calculating the real astigmatism Corn-AC and the real astigmatism axial position Corn-Ax of the front surface of the cornea through Kf, ks, kfAx.
In this step, the front surface true astigmatism is calculated by formula (3):
Corn-A Ax＝Kf Ax (3)：
the true astigmatism axis is calculated by equation (4):
taking cornea X to be corrected as an example, the calculation result is as follows:
Com-A Ax＝Kf Ax＝40
and S6, calculating total cornea astigmatism Corn C and total astigmatism axial position Corn Ax according to the real astigmatism and axial position of the front surface of the cornea and the astigmatism and axial position of the rear surface of the cornea.
Specifically, this step calculates the total corneal astigmatism axis position Corn Ax by the formula (5):
after solving the Corn Ax by the formula (5), converting the LensAx value out of the range to 0-180 by the Corn Ax + -180, and obtaining the final Corn Ax.
And calculating total corneal astigmatism Corn C by formula (6):
taking cornea X to be corrected as an example, the calculation result is as follows:
Corn Ax＝43.59；
step S7, calculating Lens astigmatism Lens C and axis position Lens Ax according to the subjective light inspection cylindrical Lens Mani C, the axis position Mani Ax and the total cornea astigmatism Corn C.
Specifically, this step calculates the Lens astigmatism axis Lens Ax by equation (7):
after solving the Lens Ax by the formula (7), converting the Lens Ax value beyond the range to 0-180 by the Lens Ax + -180, and obtaining the final Lens Ax.
And calculate Lens astigmatism Lens C by equation (8):
taking cornea X to be corrected as an example, the calculation result is as follows:
Lens Ax＝-12.64
after the Lens Ax is obtained, converting the value out of the range to 0-180 through the Lens Ax + -180, wherein the Lens ax= 167.36;
step S8, calculating all astigmatism C1 and axis Ax1 except irregular astigmatism of the cornea according to the regular astigmatism TMR C and axis TMR Ax of the front surface of the cornea, the astigmatism Corn-P C and axis Corn-P Ax of the rear surface of the cornea, the astigmatism Lens C and axis Lens Ax.
Specifically, this step calculates all astigmatism axis positions Ax1 by the formula (9):
after solving Ax1 by the formula (9), the value of Ax1 out of the range is converted to 0-180 by Ax 1+/-180, and the final Ax1 is obtained.
And all astigmatism C1 are calculated by formula (10):
taking cornea X to be corrected as an example, the calculation result is as follows:
after Ax1 is obtained, the value of Ax1 out of range is converted to 0 to 180 by ax1±180, ax1= 139.57;
step S9, calculating the target diopter according to the astigmatism and the axial position of the crystalline Lens and the residual astigmatism proportion Res-Lens C of the selected natural crystal.
The target diopter includes a percentage of residual crystalline astigmatism and an axis. After refractive correction, the lens has a degree of accommodation relaxation due to the improved vision quality. In this step, this relaxation factor of the lens is also compensated by calculation of the percent of residual astigmatism of the lens, resulting in a better correction.
Specifically, this step calculates the residual crystal astigmatism axis by formula (11):
Target Ax＝Lens Ax (11)；
and the percent residual crystalline astigmatism is calculated by equation (12):
taking cornea X to be corrected as an example, the calculation result is as follows:
targetax=lenaxrotary= 167.36;
step S10, calculating refraction correction target astigmatism AR C, axis AR Ax and sphere AR S according to all astigmatism C1 and axis Ax1 except irregular astigmatism, target diopter and subjective light inspection sphere Mani S.
Specifically, this step calculates the refractive correction target astigmatism axis AR Ax by the formula (13):
after solving the Lens Ax by the formula (13), converting the Lens Ax value beyond the range to 0-180 by the Lens Ax + -180, and obtaining the final Lens Ax value.
And calculating the refraction correction target astigmatism by the formula (14), and calculating the refraction correction target sphere AR S by the formula (15):
AR S＝Mani S-(AR C-Mani C)/2 (15)。
taking cornea X to be corrected as an example, the calculation result is as follows:
ARAx＝43.99；
AR S＝Maini S-(AR C-Manini C)/2＝-4.63。
according to the indirect compensation method of the full-higher-order aberration in the keratorefractive correction of the embodiment of the application, the actual residual diopter is 0 after the cornea X to be corrected is corrected by the finally calculated AR S= -4.63, AR C= -0.94 and AR ax=43.99.
For the same cornea X to be corrected, the FDA method is adopted, the FDA S=Mani S= -4.36, the FDA C= -1.50 and the FDA ax=65 are calculated, and after the cornea X to be corrected is corrected according to the obtained parameters, the actual residual diopter is as follows: FDAS error = 0.79, fda C error = -1.02, fda Ax error = 174.
For the same cornea X to be corrected, a TMR S= -4.53, TMR C= -1.16 and TMR ax=18 are calculated by adopting a TMR method of the prior art, and after the cornea X to be corrected is corrected according to the obtained parameters, the actual residual diopter is as follows: TMR S error= -0.19, TMR C error = 0.94, TMR Ax error = 172.
For the same cornea X to be corrected, a prior art method Phorcides method is adopted, so that Phorcides S= -4.66 and Phorcides C= -1.18,Phorcides Ax =17 are calculated, and after the cornea X to be corrected is corrected according to the obtained parameters, the actual residual diopter is as follows: phorcides S error= -0.07, phorcides C error = 0.99,Phorcides Ax error = 172.
For the same cornea X to be corrected, a ZZ VR method is adopted, ZZ VR S= -4.66, ZZ VR C= -0.96 and ZZ VR ax=52 are calculated, and after the cornea X to be corrected is corrected according to the obtained parameters, the actual residual diopter is as follows: ZZVR S error = 0.44, zz VR C error = -0.27, zz VR Ax error = 1.
According to the technical scheme, the indirect compensation method ZZ AR method of the full-order aberration in the cornea refraction correction bypasses a direct method of calculating the influence quantity of a single term of a Zernike polynomial and compensating the influence quantity one by one, so that the calculation is simplified, and the unified compensation of the influence quantity of all the high-order aberration of the Zernike is realized; because of simplifying calculation, the operation flow is obviously optimized, input items are fewer, and the introduction of human errors is reduced; meanwhile, higher-order aberration of the human eye cornea refraction correction is better and fully compensated, so that compensation precision is improved, and better visual quality is obtained.
The above description is only of the preferred embodiments of the present application and the description of the technical principles applied is not intended to limit the scope of the application as claimed, but merely represents the preferred embodiments of the present application. It will be appreciated by persons skilled in the art that the scope of the application referred to in the present application is not limited to the specific combinations of the technical features described above, but also covers other technical features formed by any combination of the technical features described above or their equivalents without departing from the inventive concept. All other embodiments, which can be made by a person skilled in the art without making any inventive effort, are intended to be within the scope of the present application.
Claims (4)
1. An indirect compensation method of full higher order aberration in cornea refraction correction comprises the following steps:
s0, obtaining a subjective refraction sphere Mani S, a cylindrical mirror Mani C and an axial position Mani Ax;
step S1, obtaining the flat axis curvature Kf, the steep axis curvature Ks and the flat axis curvature axis KfAx of the cornea to be corrected;
s2, obtaining back surface astigmatism Corn-P C and back surface astigmatism axial position Corn-P Ax of the cornea to be corrected;
s3, recording a cornea front surface regular astigmatism TMR C and a regular astigmatism axis TMR Ax in a cornea topographic map;
s4, selecting a residual astigmatism ratio Res-Lens C of the natural crystal within a range of 0-100%;
s5, calculating real astigmatism Corn-A C and real astigmatism axial position Corn-AAx of the front surface of the cornea through Kf, ks, kfAx;
s6, calculating total cornea astigmatism Corn C and total astigmatism axial position Corn Ax according to the real astigmatism and axial position of the front surface of the cornea and the astigmatism and axial position of the rear surface of the cornea;
s7, calculating Lens astigmatism Lens C and axial Lens Ax according to the subjective light-detecting cylindrical Lens Mani C, the axial Mani Ax and total cornea astigmatism Corn C;
s8, calculating all astigmatism C1 and axial position Ax1 of the cornea except irregular astigmatism according to the regular astigmatism TMR C and axial position TMR Ax of the front surface of the cornea, the astigmatism Corn-PC and axial position Corn-P Ax of the rear surface of the cornea, the astigmatism Lens C and axial position Lens Ax;
s9, calculating target diopter according to astigmatism and axial position of the crystalline Lens and the residual astigmatism proportion Res-Lens C of the selected natural crystal;
step S10, calculating refraction correction target astigmatism AR C, axis ARAx and sphere AR S according to all astigmatism C1 except irregular astigmatism and axis Ax1, target diopter and subjective optometry sphere Mani S;
the method is characterized in that the back surface astigmatism Corn-P C and the back surface astigmatism axial position Corn-P Ax of the cornea to be corrected are obtained through fixed proportion calculation of front surface astigmatism of the cornea and front and back astigmatism of the cornea of the general population; when the estimation is performed, corn-P Ax is calculated by the formula (1) according to Kf, ks, kf Ax:
Com-P Ax＝Kf Ax (1)；
Corn-P C is calculated by equation (2):
the front surface true astigmatism is calculated by formula (3):
Corn-A Ax＝Kf Ax (3)；
the true astigmatism axis is calculated by equation (4):
the total corneal astigmatism axis Corn Ax is calculated by equation (5):
after solving the Corn Ax by the formula (5), converting the Lens Ax value beyond the range to 0-180 by the Corn Ax + -180 to obtain the final Corn Ax;
and calculating total corneal astigmatism Corn C by formula (6):
the Lens astigmatism axis Lens Ax is calculated by equation (7):
solving the Lens Ax by the formula (7), and converting the Lens Ax value beyond the range to 0-180 by the Lens Ax + -180 to obtain the final Lens Ax;
and calculate Lens astigmatism Lens C by equation (8):
all astigmatism axis positions Ax1 are calculated by formula (9):
after solving Ax1 through a formula (9), converting an Ax1 value which exceeds the range to 0-180 through Ax 1+/-180 to obtain the final Ax1;
and all astigmatism C1 are calculated by formula (10):
the target diopter comprises the percentage and the axial position of residual crystal astigmatism;
the residual crystallographic astigmatism axis is calculated by formula (11):
Target Ax＝Lens Ax (11)；
and the percent residual crystalline astigmatism is calculated by equation (12):
the refractive correction target astigmatism axis AR Ax is calculated by equation (13):
solving the Lens Ax by the formula (13), and converting the Lens Ax value beyond the range to 0-180 by the Lens Ax + -180 to obtain a final Lens Ax value;
and calculating the refraction correction target astigmatism by the formula (14), and calculating the refraction correction target sphere AR S by the formula (15):
AR S＝Mani S-(AR C-Mani C)/2 (15)。
2. the method according to claim 1, wherein the data of the flat curvature Kf, steep curvature Ks and flat curvature axis KfAx introduced into the excimer device according to the corneal topography is obtained in step S1.
3. The method of indirect compensation of all higher order aberrations in refractive corneal correction according to claim 1, wherein the posterior surface astigmatism Corn-P C and the posterior surface astigmatism axial position Corn-P Ax of the cornea to be corrected are obtained from corneal topography data.
4. The method of indirect compensation for all higher order aberrations in the refractive correction of the cornea according to claim 1, characterized in that the ratio Res-Lens c=30-50% of the residual astigmatism of the natural crystal selected.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111245102.6A CN113855387B (en) | 2021-10-25 | 2021-10-25 | Indirect compensation method for full-higher-order aberration in cornea refraction correction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111245102.6A CN113855387B (en) | 2021-10-25 | 2021-10-25 | Indirect compensation method for full-higher-order aberration in cornea refraction correction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113855387A CN113855387A (en) | 2021-12-31 |
| CN113855387B true CN113855387B (en) | 2023-09-22 |
Family
ID=78997519
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111245102.6A Active CN113855387B (en) | 2021-10-25 | 2021-10-25 | Indirect compensation method for full-higher-order aberration in cornea refraction correction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113855387B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004069044A1 (en) * | 2003-02-03 | 2004-08-19 | Kabushiki Kaisha Topcon | Ophthalmic data measurement device, ophthalmic data measurement program, and eye characteristic measurement device |
| CN102566085A (en) * | 2012-03-20 | 2012-07-11 | 天津宇光光学有限公司 | Wave-front technology-based method for designing aspheric surface eyeglasses |
| CN109491083A (en) * | 2019-01-03 | 2019-03-19 | 杭州明视康眼科医院有限公司 | A kind of higher order aberratons compensation method for cornea refractive correction |
| CN113412080A (en) * | 2018-11-20 | 2021-09-17 | 马克·洛巴诺夫 | Intelligent corneal topography surgical advisor |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SE0004829D0 (en) * | 2000-12-22 | 2000-12-22 | Pharmacia Groningen Bv | Methods of obtaining ophthalmic lenses providing the eye with reduced aberrations |
| CA2604776A1 (en) * | 2005-04-14 | 2006-10-19 | University Of Rochester | System and method for treating vision refractive errors |
-
2021
- 2021-10-25 CN CN202111245102.6A patent/CN113855387B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004069044A1 (en) * | 2003-02-03 | 2004-08-19 | Kabushiki Kaisha Topcon | Ophthalmic data measurement device, ophthalmic data measurement program, and eye characteristic measurement device |
| CN102566085A (en) * | 2012-03-20 | 2012-07-11 | 天津宇光光学有限公司 | Wave-front technology-based method for designing aspheric surface eyeglasses |
| CN113412080A (en) * | 2018-11-20 | 2021-09-17 | 马克·洛巴诺夫 | Intelligent corneal topography surgical advisor |
| CN109491083A (en) * | 2019-01-03 | 2019-03-19 | 杭州明视康眼科医院有限公司 | A kind of higher order aberratons compensation method for cornea refractive correction |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113855387A (en) | 2021-12-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6817714B2 (en) | Method and apparatus relating to the optical zone of an optical element | |
| EP1597623B1 (en) | Methods for designing custom lenses for improved vision and corresponding lenses | |
| Philip et al. | Total ocular, anterior corneal and lenticular higher order aberrations in hyperopic, myopic and emmetropic eyes | |
| JP5450406B2 (en) | Method, computer readable recording medium and apparatus for checking and / or determining user data | |
| US20060279699A1 (en) | Wavefront fusion algorithms for refractive vision correction and vision diagnosis | |
| EP1943984A1 (en) | Ophthalmic lenses providing the eye with reduced aberrations | |
| US11561413B2 (en) | Population of an eye model using measurement data in order to optimize spectacle lenses | |
| DE102012000390A1 (en) | Spectacle lens optimization with individual eye model | |
| Fernández-Sánchez et al. | Effect of 3rd-order aberrations on human vision | |
| US10613347B2 (en) | Population of an eye model for optimizing spectacle lenses with measurement data | |
| Langenbucher et al. | Theoretical image performance with customized aspheric and spherical IOLs–when do we get a benefit from customized aspheric design? | |
| JP7278302B2 (en) | Adjustment of subjective and objective refraction | |
| CN113855387B (en) | Indirect compensation method for full-higher-order aberration in cornea refraction correction | |
| Charman et al. | Variations in the local refractive correction of the eye across its entrance pupil | |
| CN105339833B (en) | For determining the method and system of eyeglass prescription | |
| WO2022127030A1 (en) | Method for designing contact lens | |
| Mierdel et al. | Automated ocular wavefront analyzer for clinical use | |
| CN110267629B (en) | Calculation of actual astigmatism correction and nomogram for corneal laser treatment | |
| CN113197543B (en) | Method and system for evaluating vision quality after refractive surgery based on vector aberration theory | |
| GB2453993A (en) | Soft contact lens for correcting abnormal corneal topography | |
| Li et al. | A method to design aspheric spectacles for correction of high-order aberrations of human eye | |
| DE102005022683A1 (en) | Lens such as contact lens for compensating for aberration, has freeform rear surface compensated by freeform front surface taking account of refractive indices of media | |
| Charman | Visual Optics | |
| CN108066006A (en) | Vision correction effect simulation system and method | |
| CN118844926A (en) | A diopter compensation method for light-tracking guided corneal refractive correction |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |