US20160346120A9 - Method and Apparatus for Calculating a Laser Shot File for Use in an Excimer Laser - Google Patents
Method and Apparatus for Calculating a Laser Shot File for Use in an Excimer Laser Download PDFInfo
- Publication number
- US20160346120A9 US20160346120A9 US14/134,541 US201314134541A US2016346120A9 US 20160346120 A9 US20160346120 A9 US 20160346120A9 US 201314134541 A US201314134541 A US 201314134541A US 2016346120 A9 US2016346120 A9 US 2016346120A9
- Authority
- US
- United States
- Prior art keywords
- shot
- laser
- ablation profile
- grid
- desired ablation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- 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
- 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
- A61F9/00825—Methods or devices for eye surgery using laser for photodisruption
- A61F9/0084—Laser features or special beam parameters therefor
-
- 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
-
- 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
- A61F9/00802—Methods or devices for eye surgery using laser for photoablation
- A61F9/00804—Refractive treatments
-
- 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
- A61F9/00802—Methods or devices for eye surgery using laser for photoablation
- A61F9/00812—Inlays; Onlays; Intraocular lenses [IOL]
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23K—SOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
- B23K26/00—Working by laser beam, e.g. welding, cutting or boring
- B23K26/36—Removing material
- B23K26/38—Removing material by boring or cutting
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/14—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range characterised by the material used as the active medium
- H01S3/22—Gases
- H01S3/223—Gases the active gas being polyatomic, i.e. containing two or more atoms
- H01S3/225—Gases the active gas being polyatomic, i.e. containing two or more atoms comprising an excimer or exciplex
-
- 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
-
- 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/00878—Planning
-
- 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/00897—Scanning mechanisms or algorithms
-
- 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
- A61F9/00802—Methods or devices for eye surgery using laser for photoablation
- A61F9/00814—Laser features or special beam parameters therefor
Definitions
- the present invention relates to a method and apparatus for calculating a laser shot file for use in an excimer laser in particular using a dithering algorithm.
- the invention is specifically suitable for applying the laser shot file when performing a laser treatment of an eye or when producing a customized contact lens or an intraocular lens (IOL) by laser ablation.
- IOL intraocular lens
- U.S. Pat. No. 6,090,100 relates to an excimer laser system for correction of vision with reduced thermal effects. It specifically relates to an apparatus and method for controlling the excimer laser system for removing tissue from the eye to perform various types of corrections, such as myopia, hyperopia, and astigmatism correction.
- the excimer laser system provides a relatively large spot size which provides a relatively large coverage of treatment area per shot. While using such large spot sizes, the shots are generally not “adjacent” to each other but instead overlap to generate the desired degree of ablation at a particular point. For calculating the result of the overlapping shots, an algorithm is used.
- a dithering algorithm is used. Specific reference is made to a rectangular dithering, circular dithering and a line-by-line oriented dithering.
- shot dithering methods an array of shots is created for a fixed spot size spread over a treatment area to correct to the desired degree of ablation.
- a grid is used with a constant grid width between individual grid positions.
- the shape of the desired ablation profile which usually is a continuous profile, has to be transferred into a whole-numbered discrete density distribution.
- the continuous profile represents a planned ablation and the whole-numbered discrete density distribution represents a series of ablating flying spot laser pulses.
- the residual structure i.e., the difference between the planned and the achieved profile, has to be minimised. Exact solutions can principally be found numerically but not in a reasonable time. Therefore, for this purpose, dither algorithms are used.
- the profile is discretised on a given grid. Using a cost function or merit function the algorithm decides for each position of the grid whether to place a shot or not. For this decision, usually only a few neighboring positions of the grid are taken into account. This dither algorithm saves calculation time without the need that the real size of the spot is taken into account. It is sufficient to know the shot volume which is ablated with one laser shot.
- the known dither algorithms produce artefacts in parts of the profile, e.g., in low-density regions where the next neighboring shot is too far away. Artefacts may also be produced in high-density regions where at nearly every position, a shot is placed. The positions with no shot also have too large a distance for the assumption that only a few neighbor positions are necessary.
- known laser printer systems are using a respective fixed resolution given as a number of dots per inch, i.e., a higher number of dots per inch results in a better resolution.
- a known laser printer has no problem with overlapping and touching dots because this does not result in an additional blackening when hitting a point twice or more often. Rather, to produce an image, a certain local area of the image having a certain grey level can be created by applying a corresponding number of dots in this local area.
- the object underlying the present invention is to provide a method and apparatus for calculating a laser shot file for use in a refractive excimer laser, wherein the difference between the planned and the achieved profile is minimised.
- a desired ablation profile for correcting for example myopia has a maximum shot density in the central part of the treatment zone whereas a minimum shot density is present along the circumferential border of the treatment zone.
- the number of laser shots to be applied to the central part of the treatment zone is higher than in other sub-areas in particular along the border of the treatment zone.
- the minimum shot density is present in the central part of the treatment zone.
- the ablation profile requires a higher number of laser shots along a circumferential border of the treatment zone.
- the invention is generally applicable for any ablation profile, wherein sub-areas having different shot densities are investigated in order to determine any sub-area having a maximum shot density and/or any sub-area having a minimum shot density.
- the general concept of the present invention is based on the idea to adapt the dither algorithm which is used for placing laser shots of the excimer laser when discretizing a given ablation profile on a given grid.
- the dither algorithm decides for each position of a grid whether to place a shot. More specifically, the shot density for obtaining a predetermined desired ablation profile is calculated first. Depending on the calculated shot density of the desired ablation profile, the dither algorithm is adapted by using a dynamic threshold value being used in a cost function for the shot calculation.
- the threshold value is selected from at least two different threshold values depending on the minimum shot density and/or maximum shot density of the desired ablation profile. Generally, for a desired ablation profile having low shot densities, a lower threshold value is used. For a desired ablation profile having high shot densities, a higher threshold value is used.
- a first threshold value is a value within the range of 0% to 20% of the maximum shot density of the ablation profile.
- a second threshold value is a value within the range of 20% to 80% of the maximum shot density.
- a third threshold value is a value within the range of 80% to 100% of the maximum shot density.
- threshold value “TV(x, y)” is related to the shot density “D(x, y)” according to the following equation (1):
- the threshold value is set for each grid position in correspondence to the density function. More preferably, the threshold value is set to a value equal to or nearby the value of the density function at a respective grid position.
- the threshold value is preferably at least a value within the range of 80% to 110%, more preferably of 90% to 100% of the value of the density function at a respective grid position.
- a local shot density D(x, y) within a sub-area around a grid position P(x, y) is calculated from an ablation profile z(x, y) within the respective sub-area using the ablation volume of a single laser shot V shot and a given width G using following equation:
- a dither algorithm is used for calculating the placement of the laser shots of the excimer laser on grid positions.
- the dither algorithm is adapted to the desired ablation profile by determining the optimised grid width for the grid to be used for the dither algorithm.
- optimised grid width for the grid to be used for the dither algorithm.
- the grid width is found for a maximum value of the Profile z max (x, y) and for a desired maximum density D max (x, y):
- the local shot density around the minimum of the desired profile is calculated with a given grid width.
- the grid width is calculated with equation 4.
- the influence of the dynamic threshold is explained using two examples.
- a treatment using a treatment zone of about 5.5 mm for a desired correction of +4 dpt is selected.
- This hyperopia correction has the maximum of the ablation along an annular portion circumfering the centre.
- the desired depth is approximately 26 ⁇ m.
- About 445 laser shots are necessary to reach a result with a typical excimer treatment laser.
- a grid width of 98 ⁇ m is chosen.
- an ablation is calculated using a constant threshold.
- the treatment zone is again 5.5 mm and the correction is +4 dpt.
- the desired maximum depth is also about 26 ⁇ m and about 445 laser shots are needed.
- a dynamic threshold is used.
- the second example shows the advantage to use a dynamic threshold when calculating the ablation.
- a desired ablation profile is divided into at least two ablation sub-profiles. Then for each ablation sub-profile, the respective shot density is calculated and a respective grid width based on the respective calculated density of the ablation sub-profile is determined. Each sub-profile is calculated using the dynamic threshold.
- the calculation of the laser shot file is made in two or more runs preferably using different grid constants or grid widths for each respective ablation sub-profile resulting in a corresponding laser shot file. Thereafter, the two or more laser shot files can be combined in one single laser shot file.
- the calculated, placed laser shots are processed in a further step of sorting to obtain a laser shot sequence.
- the sorting is performed taking into consideration that any thermal effects should be avoided, i.e., two consecutive laser shots are preferably placed on different grid positions in the treatment zone which are at a distance from each other.
- FIG. 1A is a diagram showing the location of laser spots for a first test using a constant threshold value
- FIG. 1B is a diagram showing the planned and the achieved profile as a cross-section along the horizontal axis of FIG. 1A ,
- FIG. 1C is a diagram showing the planned and the achieved profile as a cross-section along the vertical axis of FIG. 1A ,
- FIG. 2A is a diagram showing the location of laser spots for a second test using a dynamic threshold value according to a preferred embodiment of the present invention
- FIG. 2B is a diagram showing the planned and the achieved profile as a cross-section along the horizontal axis of FIG. 2A ,
- FIG. 2C is a diagram showing the planned and the achieved profile as a cross-section along the vertical axis of FIG. 2A ,
- FIG. 3 shows a flow diagram with a calculation of laser pulse patterns with a dither algorithm
- FIG. 4 shows an example of a sub-grid with weighting factors usable for weighting neighboring error values.
- FIGS. 1A, 1B and 1C show the simulated calculation of a laser shot file for use in an excimer laser for the correction of hyperopia with a value of about +4 dioptres, using a typical excimer laser for refractive treatments, within a treatment zone having a diameter of 5.5 mm and using a laser spot having a diameter of 1 mm.
- the grid width is 98 ⁇ m.
- the distance between two neighbor grid points is 98 ⁇ m.
- the grid points are arranged in rows and columns.
- 445 laser shots are used for achieving the ablation.
- the resulting treatment is expected to have a refraction of said about +4 dioptres.
- the diagram of FIG. 1A shows the respective centre position of each of the 445 laser shots which is related to one of the grid positions each marked with a “+”-sign.
- the grid is schematically shown having a grid width of 98 ⁇ m.
- Each of the shown laser shot centre positions are arranged on a grid point of this grid.
- the diagram of FIG. 1B shows, as a dashed line, the desired ablation profile, i.e., the ablation depth in ⁇ m with respect to a respective X-position.
- the ablation depth is approximately 26 ⁇ m in an annular portion of the treatment zone at about the x-positions ⁇ 2 and +2 and is smaller in the centre portion and at both sides.
- the ablation depth is almost zero in the centre portion. It further shows the simulated resulting ablation profile as a continuous line as a cross-section taken along the horizontal axis through the point 0-0 in FIG. 1A .
- FIG. 1C shows the desired ablation profile as a dashed line which is taken as a cross-section along the vertical axis through the point 0-0 in FIG. 1A .
- FIG. 1C further shows the resulting ablation profile as a continuous line taken as a cross-section along the vertical axis through point 0-0 of FIG. 1 A.
- the average shot density inside the treatment zone having in this example a diameter of 5.5 mm, is about 18% ( FIG. 1A ).
- the respective centre positions of the laser shots are placed within a range of ⁇ 2.7 mm in the X-direction and ⁇ 2.7 mm in the Y-direction.
- the use of a constant threshold value for the first test causes artefacts like the linewise arrangement of laser shot positions in the lower part of the ablation ( FIG. 1A ). For example, as shown, several laser shots are provided at grid positions which are arranged along a horizontal bottom line at a closer distance. Further laser shots are provided at grid positions which are arranged at a larger distance from this horizontal bottom line. Thus, the laser shots are not provided in an equal manner resulting in a deviation from the desired ablation profile (see FIG. 1C ).
- FIG. 1C shows that the resulting ablation profile deviates from the desired ablation profile, i.e., there is a shift with respect to the right part of the desired ablation profile.
- the dither algorithm produces artefacts in parts of the resulting ablation profile which may depend on the order of the calculation of the laser shots for respective grid positions. In regions with gradients in shot density, the shots are shifted. The shift depends on the depth of the desired ablation. Additionally, artefacts called worms may be introduced.
- the input parameters are the shot volume of a laser shot and the desired ablation profile. There is no need to take the beam diameter into account as the dither algorithm works independently therefrom.
- the dither algorithm provides a laser shot file as an output. More specifically, the dither algorithm is used for the placement of laser shots of the excimer laser on grid positions. Preferably, a cost function is used for deciding for each grid position whether a laser shot is placed or not. Herein, preferably the decision is made with regard to whether one or more laser shot(s) is (are) placed at a grid position(s) within the neighborhood of the given grid position.
- a dither algorithm is used as disclosed in U.S. Pat. No. 6,090,100.
- FIG. 3 shows a flow chart representing an example for the error diffusion.
- This dither algorithm is based on the concept of error diffusion.
- the desired ablation profile is calculated based, e.g., on the desired correction of a patient's eye or the modification of contact lenses or of IOLs.
- This profile is stored within a grid having a specific grid width. For example, such a grid has 256 ⁇ 256 values which covers an area of 15 2 mm 2 .
- the error diffusion may be started in one edge within that grid and follows it line by line.
- a first step S 1 the ablation profile and the dynamic threshold is determined using equation 1 and the active dither position is set to a point in one of the edges of the grid.
- a desired grid width is calculated.
- Said active dither position represents the actual position within the grid being processed.
- a desired ablation value for the active dither position is obtained.
- this desired ablation value is multiplied with a scaling factor f.
- the scaling factor f takes into account the different size of laser pulse and the positioning step, i.e., the grid width. More specifically, the scaling factor is calculated as follows to get the desired shot density at this position (see equation 3):
- the area of the smallest square the laser beam can be sent around is (58 ⁇ m) 2 .
- the number of calculated pulses are reduced in order to take into account for the overlapping of laser pulses.
- weighted neighboring errors are added to the scaled desired ablation value for the active dither position.
- These weighted neighboring errors are preferably the weighted sum of errors of adjacent grid points that have already been processed. An example will be described later.
- a predetermined threshold the sum of the value for the respective grid point and the weighted errors of adjacent grid points will be compared to this threshold value. If the value is not larger than the dynamic threshold T (x, y) step S 9 follows. If the value is larger than the threshold, a laser pulse is set for this grid position in step S 6 . One laser pulse is subtracted from said density value. Then in step S 7 it is determined whether a new value is still larger than the threshold. In case the new value is larger than the dynamic threshold in step S 8 , it is determined that an error of shot overflow occurred.
- this new value is stored as an error for this particular grid position. It will be used when processing neighboring positions for the calculation with respect to further dither positions.
- next step S 10 it is decided whether the line is complete; if not, in step S 11 a next point in the same line is selected as an active position and the before-mentioned processing is repeated.
- step S 12 a decision has to be made whether there is a new line; if yes, then in step S 13 a first point in the new line is selected as active position and the processing is repeated. Otherwise, if there is no new line, the processing ends with step S 14 .
- the before-mentioned grid point error represents the ablation error done at a particular grid point. For each grid point processed, this error is the sum of desired ablation value plus the weighted neighboring errors minus the laser pulse ablation depth (if a laser pulse has been set for that position).
- FIG. 4 shows an example for weighting of errors of neighboring grid points. More specifically, FIG. 4 shows a sub-grid of 7 ⁇ 7 grid points, wherein the active dither position is shown in the middle. In this case, the weighting function is determined as 8/distance with a distance measured in units of grid points. The sum of the errors will then be normalised by a division with 70.736 which is the sum of all weighting factors used. As apparent from FIG. 4 , the white positions indicate grid position not yet processed. Thus, before deciding whether a laser pulse has to be set at a given grid position, the error induced while processing adjacent grid points has to be added to the theoretical ablation value for that grid point.
- the errors of the neighboring grid points are not simply added but weighted due to their distance to the active grid point.
- the respective weighting factors are shown in FIG. 4 . It shall be noted that this is just one possible method for summing up the surrounding errors, which is working fine.
- a laser shot sequence may be determined thereafter by using a separate sorting algorithm.
- a sorting may be performed in order to avoid thermal effects.
- any two following laser shots should preferably be placed at two grid positions at a distance from each other.
- every four shots a laser shot is placed in the same region as the first shot.
Abstract
Description
- This is a continuation of U.S. patent application Ser. No. 12/375,428, filed on Apr. 21, 2009, and claims priority to PCT/EP07/57782, filed on Feb. 7, 2008.
- The present invention relates to a method and apparatus for calculating a laser shot file for use in an excimer laser in particular using a dithering algorithm. The invention is specifically suitable for applying the laser shot file when performing a laser treatment of an eye or when producing a customized contact lens or an intraocular lens (IOL) by laser ablation.
- U.S. Pat. No. 6,090,100 relates to an excimer laser system for correction of vision with reduced thermal effects. It specifically relates to an apparatus and method for controlling the excimer laser system for removing tissue from the eye to perform various types of corrections, such as myopia, hyperopia, and astigmatism correction. In one disclosed embodiment, the excimer laser system provides a relatively large spot size which provides a relatively large coverage of treatment area per shot. While using such large spot sizes, the shots are generally not “adjacent” to each other but instead overlap to generate the desired degree of ablation at a particular point. For calculating the result of the overlapping shots, an algorithm is used. In one method of calculating treatment patterns using large, fixed spot sizes distributed throughout the treatment area, a dithering algorithm is used. Specific reference is made to a rectangular dithering, circular dithering and a line-by-line oriented dithering. Using any variety of shot dithering methods, an array of shots is created for a fixed spot size spread over a treatment area to correct to the desired degree of ablation. For the respective array, a grid is used with a constant grid width between individual grid positions. With the known dither methods, the shape of the desired ablation profile, which usually is a continuous profile, has to be transferred into a whole-numbered discrete density distribution. Here, the continuous profile represents a planned ablation and the whole-numbered discrete density distribution represents a series of ablating flying spot laser pulses. The residual structure, i.e., the difference between the planned and the achieved profile, has to be minimised. Exact solutions can principally be found numerically but not in a reasonable time. Therefore, for this purpose, dither algorithms are used. The profile is discretised on a given grid. Using a cost function or merit function the algorithm decides for each position of the grid whether to place a shot or not. For this decision, usually only a few neighboring positions of the grid are taken into account. This dither algorithm saves calculation time without the need that the real size of the spot is taken into account. It is sufficient to know the shot volume which is ablated with one laser shot. However, under certain conditions, the known dither algorithms produce artefacts in parts of the profile, e.g., in low-density regions where the next neighboring shot is too far away. Artefacts may also be produced in high-density regions where at nearly every position, a shot is placed. The positions with no shot also have too large a distance for the assumption that only a few neighbor positions are necessary.
- As regards the general background of dithering algorithms, reference is made to U.S. Pat. No. 6,271,936 B1, which relates to the field of digital image processing. It particularly relates to a method for digitally multitoning a continuous-tone image using error diffusion, dithering and over-modulation methods. Reference is made to the problem that an artefact may occur like worms which are formed when the black or white output pixels appear to string together in an area that should be otherwise uniform. Wherein this U.S. patent gives a detailed description of these known methods, it is related to a completely different technical field. Among other differences, known laser printer systems are using a respective fixed resolution given as a number of dots per inch, i.e., a higher number of dots per inch results in a better resolution. Moreover, a known laser printer has no problem with overlapping and touching dots because this does not result in an additional blackening when hitting a point twice or more often. Rather, to produce an image, a certain local area of the image having a certain grey level can be created by applying a corresponding number of dots in this local area.
- The object underlying the present invention is to provide a method and apparatus for calculating a laser shot file for use in a refractive excimer laser, wherein the difference between the planned and the achieved profile is minimised. This object is solved with the features of the claims.
- A desired ablation profile for correcting for example myopia has a maximum shot density in the central part of the treatment zone whereas a minimum shot density is present along the circumferential border of the treatment zone. Thus, the number of laser shots to be applied to the central part of the treatment zone is higher than in other sub-areas in particular along the border of the treatment zone.
- For the correction of, for example, hyperopia the minimum shot density is present in the central part of the treatment zone. On the other hand, the ablation profile requires a higher number of laser shots along a circumferential border of the treatment zone.
- The invention is generally applicable for any ablation profile, wherein sub-areas having different shot densities are investigated in order to determine any sub-area having a maximum shot density and/or any sub-area having a minimum shot density.
- The general concept of the present invention is based on the idea to adapt the dither algorithm which is used for placing laser shots of the excimer laser when discretizing a given ablation profile on a given grid. Using a cost function the dither algorithm decides for each position of a grid whether to place a shot. More specifically, the shot density for obtaining a predetermined desired ablation profile is calculated first. Depending on the calculated shot density of the desired ablation profile, the dither algorithm is adapted by using a dynamic threshold value being used in a cost function for the shot calculation.
- According to a preferred embodiment of the present invention, the threshold value is selected from at least two different threshold values depending on the minimum shot density and/or maximum shot density of the desired ablation profile. Generally, for a desired ablation profile having low shot densities, a lower threshold value is used. For a desired ablation profile having high shot densities, a higher threshold value is used.
- According to a preferred embodiment of the present invention, a first threshold value is a value within the range of 0% to 20% of the maximum shot density of the ablation profile. Alternatively or in addition a second threshold value is a value within the range of 20% to 80% of the maximum shot density. Alternatively or in addition a third threshold value is a value within the range of 80% to 100% of the maximum shot density.
- According to a further preferred embodiment of the present invention, more than three different threshold values are used and more preferably the threshold value “TV(x, y)” is related to the shot density “D(x, y)” according to the following equation (1):
-
TV(x,y)=f(D(x,y)) (1) - More preferably, there is a linear relationship between the threshold value “TV(x, y)” and “D(x, y)” according to the following equation (2):
-
TV(x,y)=a·D(x,y) (2), - wherein “a” is a positive factor within the range of 0<a≦1.5 and wherein “x” and “y” are the coordinates of the grid position for which the calculation is made.
- Preferably, the threshold value is set for each grid position in correspondence to the density function. More preferably, the threshold value is set to a value equal to or nearby the value of the density function at a respective grid position.
- The threshold value is preferably at least a value within the range of 80% to 110%, more preferably of 90% to 100% of the value of the density function at a respective grid position. Thus, the factor “a” of equation (1) is a value preferably within the range of 0.8 to 1.1, more preferably of 0.9 to 1.0. Best results can be achieved with a=1.
- A local shot density D(x, y) within a sub-area around a grid position P(x, y) is calculated from an ablation profile z(x, y) within the respective sub-area using the ablation volume of a single laser shot Vshot and a given width G using following equation:
-
D(x,y)=z(xy)*G 2 /V Shot (3) - According to a preferred embodiment, a dither algorithm is used for calculating the placement of the laser shots of the excimer laser on grid positions. The dither algorithm is adapted to the desired ablation profile by determining the optimised grid width for the grid to be used for the dither algorithm. For a more detailed description of this aspect to optimize the grid width, reference is made to the co-pending patent application of the present applicant with the title “Method and apparatus for calculating a laser shot file for use in a refractive excimer laser”.
- According to a preferred embodiment, with following equation the grid width is found for a maximum value of the Profile zmax(x, y) and for a desired maximum density Dmax (x, y):
-
G=√{square root over (V shot *D max(x,y)/z max(x,y))} (4) - With equation 3 the local shot density around the minimum of the desired profile is calculated with a given grid width. Preferably the grid width is calculated with
equation 4. The influence of the dynamic threshold is explained using two examples. As a first example a treatment using a treatment zone of about 5.5 mm for a desired correction of +4 dpt is selected. This hyperopia correction has the maximum of the ablation along an annular portion circumfering the centre. The desired depth is approximately 26 μm. About 445 laser shots are necessary to reach a result with a typical excimer treatment laser. To get shot densities along the annular portion of about 18% a grid width of 98 μm is chosen. In this example an ablation is calculated using a constant threshold. In a second example of an ablation the treatment zone is again 5.5 mm and the correction is +4 dpt. The desired maximum depth is also about 26 μm and about 445 laser shots are needed. For the second example, a dynamic threshold is used. The second example shows the advantage to use a dynamic threshold when calculating the ablation. - According to a further preferred embodiment, a desired ablation profile is divided into at least two ablation sub-profiles. Then for each ablation sub-profile, the respective shot density is calculated and a respective grid width based on the respective calculated density of the ablation sub-profile is determined. Each sub-profile is calculated using the dynamic threshold. Thus, for a desired ablation profile where the contrast is too high, i.e., a difference between the maximum shot density and the minimum shot density is too high, the calculation of the laser shot file is made in two or more runs preferably using different grid constants or grid widths for each respective ablation sub-profile resulting in a corresponding laser shot file. Thereafter, the two or more laser shot files can be combined in one single laser shot file.
- According to the present invention, the calculated, placed laser shots are processed in a further step of sorting to obtain a laser shot sequence. The sorting is performed taking into consideration that any thermal effects should be avoided, i.e., two consecutive laser shots are preferably placed on different grid positions in the treatment zone which are at a distance from each other.
-
FIG. 1A is a diagram showing the location of laser spots for a first test using a constant threshold value, -
FIG. 1B is a diagram showing the planned and the achieved profile as a cross-section along the horizontal axis ofFIG. 1A , -
FIG. 1C is a diagram showing the planned and the achieved profile as a cross-section along the vertical axis ofFIG. 1A , -
FIG. 2A is a diagram showing the location of laser spots for a second test using a dynamic threshold value according to a preferred embodiment of the present invention, -
FIG. 2B is a diagram showing the planned and the achieved profile as a cross-section along the horizontal axis ofFIG. 2A , -
FIG. 2C is a diagram showing the planned and the achieved profile as a cross-section along the vertical axis ofFIG. 2A , -
FIG. 3 shows a flow diagram with a calculation of laser pulse patterns with a dither algorithm, and -
FIG. 4 shows an example of a sub-grid with weighting factors usable for weighting neighboring error values. -
FIGS. 1A, 1B and 1C show the simulated calculation of a laser shot file for use in an excimer laser for the correction of hyperopia with a value of about +4 dioptres, using a typical excimer laser for refractive treatments, within a treatment zone having a diameter of 5.5 mm and using a laser spot having a diameter of 1 mm. In this simulated first test, the grid width is 98 μm. Thus, the distance between two neighbor grid points is 98 μm. In this example, the grid points are arranged in rows and columns. In total, 445 laser shots are used for achieving the ablation. Depending on the ablated volume of a single shot the resulting treatment is expected to have a refraction of said about +4 dioptres. The diagram ofFIG. 1A shows the respective centre position of each of the 445 laser shots which is related to one of the grid positions each marked with a “+”-sign. In the upper right comer ofFIG. 1A , the grid is schematically shown having a grid width of 98 μm. Each of the shown laser shot centre positions are arranged on a grid point of this grid. The diagram ofFIG. 1B shows, as a dashed line, the desired ablation profile, i.e., the ablation depth in μm with respect to a respective X-position. The ablation depth is approximately 26 μm in an annular portion of the treatment zone at about the x-positions −2 and +2 and is smaller in the centre portion and at both sides. The ablation depth is almost zero in the centre portion. It further shows the simulated resulting ablation profile as a continuous line as a cross-section taken along the horizontal axis through the point 0-0 inFIG. 1A . Similarly,FIG. 1C shows the desired ablation profile as a dashed line which is taken as a cross-section along the vertical axis through the point 0-0 inFIG. 1A .FIG. 1C further shows the resulting ablation profile as a continuous line taken as a cross-section along the vertical axis through point 0-0 ofFIG. 1 A. InFIG. 1 the average shot density inside the treatment zone, having in this example a diameter of 5.5 mm, is about 18% (FIG. 1A ). The respective centre positions of the laser shots are placed within a range of ±2.7 mm in the X-direction and ±2.7 mm in the Y-direction. -
FIGS. 2A, 2B and 2C show the results of a similar second test as inFIGS. 1A, 1B and 1C except for using a dynamic threshold. More specifically in this test the shot density D (x, y) has been used as the threshold value TV (x, y). Thus, the factor “a” in the above equation (2) is selected as a=1. - The use of a constant threshold value for the first test causes artefacts like the linewise arrangement of laser shot positions in the lower part of the ablation (
FIG. 1A ). For example, as shown, several laser shots are provided at grid positions which are arranged along a horizontal bottom line at a closer distance. Further laser shots are provided at grid positions which are arranged at a larger distance from this horizontal bottom line. Thus, the laser shots are not provided in an equal manner resulting in a deviation from the desired ablation profile (seeFIG. 1C ). - A comparison of the Figures for the first test and the second test shows that the resulting ablation profile in the second test is better, i.e., the curve of the resulting ablation profile better follows to the curve of the desired ablation profile (see
FIGS. 2B and 2C ). In particular,FIG. 1C shows that the resulting ablation profile deviates from the desired ablation profile, i.e., there is a shift with respect to the right part of the desired ablation profile. The dither algorithm produces artefacts in parts of the resulting ablation profile which may depend on the order of the calculation of the laser shots for respective grid positions. In regions with gradients in shot density, the shots are shifted. The shift depends on the depth of the desired ablation. Additionally, artefacts called worms may be introduced. - By using a dither algorithm, the input parameters are the shot volume of a laser shot and the desired ablation profile. There is no need to take the beam diameter into account as the dither algorithm works independently therefrom. The dither algorithm provides a laser shot file as an output. More specifically, the dither algorithm is used for the placement of laser shots of the excimer laser on grid positions. Preferably, a cost function is used for deciding for each grid position whether a laser shot is placed or not. Herein, preferably the decision is made with regard to whether one or more laser shot(s) is (are) placed at a grid position(s) within the neighborhood of the given grid position. Preferably, a dither algorithm is used as disclosed in U.S. Pat. No. 6,090,100.
- In the following a preferred dithering algorithm will be described with reference to
FIG. 3 which shows a flow chart representing an example for the error diffusion. This dither algorithm is based on the concept of error diffusion. Prior to the step of error diffusion, the desired ablation profile is calculated based, e.g., on the desired correction of a patient's eye or the modification of contact lenses or of IOLs. This profile is stored within a grid having a specific grid width. For example, such a grid has 256×256 values which covers an area of 152 mm2. The error diffusion may be started in one edge within that grid and follows it line by line. - In a first step S1, the ablation profile and the dynamic threshold is determined using
equation 1 and the active dither position is set to a point in one of the edges of the grid. Optionally, a desired grid width is calculated. Said active dither position represents the actual position within the grid being processed. - In a next step S2, a desired ablation value for the active dither position is obtained. In step S3 this desired ablation value is multiplied with a scaling factor f. The scaling factor f takes into account the different size of laser pulse and the positioning step, i.e., the grid width. More specifically, the scaling factor is calculated as follows to get the desired shot density at this position (see equation 3):
-
- For the above-mentioned grid having 256×256 values covering an area of 152 mm2, the grid width is 15 mm/256=58 μm. Thus, the area of the smallest square the laser beam can be sent around is (58 μm)2. Thus, the number of calculated pulses are reduced in order to take into account for the overlapping of laser pulses.
- In a next step S4, weighted neighboring errors are added to the scaled desired ablation value for the active dither position. These weighted neighboring errors are preferably the weighted sum of errors of adjacent grid points that have already been processed. An example will be described later.
- In a further step S5, a decision is made whether the obtained value is larger than a predetermined threshold. Thus, the sum of the value for the respective grid point and the weighted errors of adjacent grid points will be compared to this threshold value. If the value is not larger than the dynamic threshold T (x, y) step S9 follows. If the value is larger than the threshold, a laser pulse is set for this grid position in step S6. One laser pulse is subtracted from said density value. Then in step S7 it is determined whether a new value is still larger than the threshold. In case the new value is larger than the dynamic threshold in step S8, it is determined that an error of shot overflow occurred. In other words, if at a grid position it would be necessary to set more than laser pulse, the algorithm has to stop with an error. By the use of grid width calculated with
equation 4 this error can be avoided. In this exemplary implementation of the error diffusion, a maximum of one laser pulse for each grid position is allowed. - On the other hand, if the new value is not larger than a dynamic threshold in step S9, this new value is stored as an error for this particular grid position. It will be used when processing neighboring positions for the calculation with respect to further dither positions.
- In the next step S10, it is decided whether the line is complete; if not, in step S11 a next point in the same line is selected as an active position and the before-mentioned processing is repeated. In case the line is complete, then in step S12 a decision has to be made whether there is a new line; if yes, then in step S13 a first point in the new line is selected as active position and the processing is repeated. Otherwise, if there is no new line, the processing ends with step S14. The before-mentioned grid point error represents the ablation error done at a particular grid point. For each grid point processed, this error is the sum of desired ablation value plus the weighted neighboring errors minus the laser pulse ablation depth (if a laser pulse has been set for that position).
-
FIG. 4 shows an example for weighting of errors of neighboring grid points. More specifically,FIG. 4 shows a sub-grid of 7×7 grid points, wherein the active dither position is shown in the middle. In this case, the weighting function is determined as 8/distance with a distance measured in units of grid points. The sum of the errors will then be normalised by a division with 70.736 which is the sum of all weighting factors used. As apparent fromFIG. 4 , the white positions indicate grid position not yet processed. Thus, before deciding whether a laser pulse has to be set at a given grid position, the error induced while processing adjacent grid points has to be added to the theoretical ablation value for that grid point. The errors of the neighboring grid points are not simply added but weighted due to their distance to the active grid point. The respective weighting factors are shown inFIG. 4 . It shall be noted that this is just one possible method for summing up the surrounding errors, which is working fine. - It shall be noted that the above described dither algorithm is only one example for using the present invention.
- A laser shot sequence may be determined thereafter by using a separate sorting algorithm. A sorting may be performed in order to avoid thermal effects. Thus, any two following laser shots should preferably be placed at two grid positions at a distance from each other. Preferably, every four shots a laser shot is placed in the same region as the first shot.
- The foregoing disclosure and description of the invention are illustrative and explanatory thereof and changes in the construction and method of operation may be made without departing from the scope of the invention.
Claims (15)
TV(x,y)=f(D(x,y)).
TV(x,y)=a·D(x,y),
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/134,541 US9636256B2 (en) | 2006-08-02 | 2013-12-19 | Method and apparatus for calculating a laser shot file for use in an excimer laser |
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102006036085.0 | 2006-08-02 | ||
DE102006036085A DE102006036085A1 (en) | 2006-08-02 | 2006-08-02 | Method and apparatus for calculating a laser shot file for use in an excimer laser |
PCT/EP2007/057782 WO2008015175A1 (en) | 2006-08-02 | 2007-07-27 | Method and apparatus for calculating a laser shot file for use in an excimer laser |
US37542809A | 2009-01-28 | 2009-01-28 | |
US14/134,541 US9636256B2 (en) | 2006-08-02 | 2013-12-19 | Method and apparatus for calculating a laser shot file for use in an excimer laser |
Related Parent Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12375428 Continuation | 2007-02-27 | ||
PCT/EP2007/057782 Continuation WO2008015175A1 (en) | 2006-08-02 | 2007-07-27 | Method and apparatus for calculating a laser shot file for use in an excimer laser |
US12/375,428 Continuation US20090306635A1 (en) | 2006-08-02 | 2008-02-07 | Method and apparatus for calculating a laser shot file for use in an excimer laser |
Publications (3)
Publication Number | Publication Date |
---|---|
US20150173946A1 US20150173946A1 (en) | 2015-06-25 |
US20160346120A9 true US20160346120A9 (en) | 2016-12-01 |
US9636256B2 US9636256B2 (en) | 2017-05-02 |
Family
ID=38486132
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/375,428 Abandoned US20090306635A1 (en) | 2006-08-02 | 2008-02-07 | Method and apparatus for calculating a laser shot file for use in an excimer laser |
US14/134,541 Active 2027-06-02 US9636256B2 (en) | 2006-08-02 | 2013-12-19 | Method and apparatus for calculating a laser shot file for use in an excimer laser |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/375,428 Abandoned US20090306635A1 (en) | 2006-08-02 | 2008-02-07 | Method and apparatus for calculating a laser shot file for use in an excimer laser |
Country Status (10)
Country | Link |
---|---|
US (2) | US20090306635A1 (en) |
EP (1) | EP2049053B1 (en) |
JP (1) | JP2009545349A (en) |
KR (1) | KR101401823B1 (en) |
CN (1) | CN101677873B (en) |
AU (1) | AU2007280515B2 (en) |
CA (1) | CA2660468C (en) |
DE (2) | DE102006036085A1 (en) |
ES (1) | ES2337523T3 (en) |
WO (1) | WO2008015175A1 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102006036085A1 (en) * | 2006-08-02 | 2008-02-07 | Bausch & Lomb Incorporated | Method and apparatus for calculating a laser shot file for use in an excimer laser |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6090100A (en) * | 1992-10-01 | 2000-07-18 | Chiron Technolas Gmbh Ophthalmologische Systeme | Excimer laser system for correction of vision with reduced thermal effects |
US20030023233A1 (en) * | 2001-07-30 | 2003-01-30 | Smith Michael J. | Technique for removal of material to achieve a desired shape with a laser |
Family Cites Families (61)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5425727A (en) * | 1988-04-01 | 1995-06-20 | Koziol; Jeffrey E. | Beam delivery system and method for corneal surgery |
US5984916A (en) * | 1993-04-20 | 1999-11-16 | Lai; Shui T. | Ophthalmic surgical laser and method |
US6325792B1 (en) | 1991-11-06 | 2001-12-04 | Casimir A. Swinger | Ophthalmic surgical laser and method |
DE4232915A1 (en) * | 1992-10-01 | 1994-04-07 | Hohla Kristian | Device for shaping the cornea by removing tissue |
WO1995027535A1 (en) | 1994-04-08 | 1995-10-19 | Summit Technology, Inc. | Profiling the intensity distribution of optical beams |
US6302877B1 (en) * | 1994-06-29 | 2001-10-16 | Luis Antonio Ruiz | Apparatus and method for performing presbyopia corrective surgery |
EP0697611B9 (en) | 1994-08-18 | 2003-01-22 | Carl Zeiss | Optical coherence tomography assisted surgical apparatus |
US6454761B1 (en) * | 1995-01-30 | 2002-09-24 | Philip D. Freedman | Laser surgery device and method |
US5891132A (en) * | 1996-05-30 | 1999-04-06 | Chiron Technolas Gmbh Opthalmologische Systeme | Distributed excimer laser surgery system |
DE19727573C1 (en) | 1996-10-26 | 1998-05-20 | Aesculap Meditec Gmbh | Device and method for shaping surfaces, in particular lenses |
US5777719A (en) * | 1996-12-23 | 1998-07-07 | University Of Rochester | Method and apparatus for improving vision and the resolution of retinal images |
WO1998048746A1 (en) | 1997-04-25 | 1998-11-05 | Technolas Gmbh Ophthalmologische Systeme | Dual mode ophthalmic laser ablation |
US5928221A (en) * | 1997-11-17 | 1999-07-27 | Coherent, Inc. | Fluence monitoring method for laser treatment of biological tissue |
US6132424A (en) * | 1998-03-13 | 2000-10-17 | Lasersight Technologies Inc. | Smooth and uniform laser ablation apparatus and method |
EP0947158B1 (en) * | 1998-03-31 | 2006-05-24 | Nidek Co., Ltd. | Ophthalmic apparatus |
US6203407B1 (en) * | 1998-09-03 | 2001-03-20 | Micron Technology, Inc. | Method and apparatus for increasing-chemical-polishing selectivity |
US6149643A (en) * | 1998-09-04 | 2000-11-21 | Sunrise Technologies International, Inc. | Method and apparatus for exposing a human eye to a controlled pattern of radiation |
AU5907299A (en) | 1998-09-04 | 2000-03-27 | Sunrise Technologies International, Inc. | Radiation treatment system and methods for correcting vision |
UA67870C2 (en) * | 2002-10-04 | 2004-07-15 | Сергій Васильович Молебний | Method for measuring wave aberrations of eyes |
US6271936B1 (en) * | 1998-12-11 | 2001-08-07 | Eastman Kodak Company | Combining error diffusion, dithering and over-modulation for smooth multilevel printing |
US6332216B1 (en) * | 1999-03-09 | 2001-12-18 | Hewlett-Packard Company | Hybrid just-in-time compiler that consumes minimal resource |
US6129722A (en) * | 1999-03-10 | 2000-10-10 | Ruiz; Luis Antonio | Interactive corrective eye surgery system with topography and laser system interface |
DE19938203A1 (en) * | 1999-08-11 | 2001-02-15 | Aesculap Meditec Gmbh | Method and device for correcting visual defects in the human eye |
ES2306667T3 (en) * | 1999-09-10 | 2008-11-16 | Haag-Streit Ag | DEVICE FOR THE PHOTO-ABLATION OF THE CORNEA THROUGH A LASER RAY. |
US6086204A (en) * | 1999-09-20 | 2000-07-11 | Magnante; Peter C. | Methods and devices to design and fabricate surfaces on contact lenses and on corneal tissue that correct the eye's optical aberrations |
US6322216B1 (en) * | 1999-10-07 | 2001-11-27 | Visx, Inc | Two camera off-axis eye tracker for laser eye surgery |
US7237898B1 (en) | 1999-10-21 | 2007-07-03 | Bausch & Lomb Incorporated | Customized corneal profiling |
CA2388014C (en) | 1999-10-21 | 2013-04-16 | Technolas Gmbh Ophthalmologische Systeme | Multi-step laser correction of ophthalmic refractive errors |
DE19958436B4 (en) * | 1999-12-03 | 2014-07-17 | Carl Zeiss Meditec Ag | Apparatus and method for active, physiologically evaluated, comprehensive correction of the aberrations of the human eye |
US6315771B1 (en) * | 1999-12-09 | 2001-11-13 | Nidek Co., Ltd. | Apparatus for corneal surgery |
US6550917B1 (en) * | 2000-02-11 | 2003-04-22 | Wavefront Sciences, Inc. | Dynamic range extension techniques for a wavefront sensor including use in ophthalmic measurement |
US6394999B1 (en) * | 2000-03-13 | 2002-05-28 | Memphis Eye & Cataract Associates Ambulatory Surgery Center | Laser eye surgery system using wavefront sensor analysis to control digital micromirror device (DMD) mirror patterns |
US6923802B2 (en) * | 2000-03-13 | 2005-08-02 | Memphis Eye & Cataract Assoc. | System for generating ablation profiles for laser refractive eye surgery |
JP4021136B2 (en) * | 2000-08-31 | 2007-12-12 | 株式会社ニデック | Cornea surgery device |
DE60121976D1 (en) * | 2000-10-10 | 2006-09-14 | Univ Rochester | DETERMINATION OF THE OKULAR REFECTION BY WAVE FRONT BRAKING DATA |
WO2002034178A1 (en) * | 2000-10-20 | 2002-05-02 | Bausch & Lomb Incorporated | Method and system for improving vision |
UA59488C2 (en) * | 2001-10-03 | 2003-09-15 | Василь Васильович Молебний | Method for measuring wave aberrations of eye and device for its realization (variants) |
SG154324A1 (en) * | 2001-04-18 | 2009-08-28 | Bausch & Lomb | Objective manifest refraction |
US6827442B2 (en) * | 2001-09-12 | 2004-12-07 | Denwood F. Ross | Ophthalmic wavefront measuring devices |
US6554429B1 (en) * | 2001-10-15 | 2003-04-29 | Alcon, Inc. | Method for determining accommodation |
JP4343699B2 (en) | 2002-02-11 | 2009-10-14 | ヴィズイクス・インコーポレーテッド | Closed loop system and method for cutting a lens with aberration |
US20050107775A1 (en) | 2002-03-04 | 2005-05-19 | The Cleveland Clinic Foundation | Method and apparatus for controlling ablation in refractive surgery |
US20040002697A1 (en) * | 2002-06-27 | 2004-01-01 | Gerhard Youssefi | Biconic ablation with controlled spherical aberration |
AU2003249884B2 (en) * | 2002-06-27 | 2008-08-14 | Technolas Gmbh Ophthalmologishe Systeme | Myopia correction enhancing biodynamic ablation |
US7133137B2 (en) * | 2002-06-27 | 2006-11-07 | Visx, Incorporated | Integrated scanning and ocular tomography system and method |
JP2004148074A (en) | 2002-09-06 | 2004-05-27 | Nidek Co Ltd | Cornea surgery apparatus |
ITMI20022332A1 (en) * | 2002-11-04 | 2004-05-05 | Univ Padova | LASER EQUIPMENT FOR THE TREATMENT OF HARD AND RELATED TISSUES |
ITTO20021007A1 (en) | 2002-11-19 | 2004-05-20 | Franco Bartoli | EXCIMER LASER EQUIPMENT AND DRIVING METHOD |
US7293873B2 (en) | 2002-12-06 | 2007-11-13 | Visx, Incorporated | Presbyopia correction using patient data |
CA2522787C (en) * | 2003-04-18 | 2012-09-18 | Visx, Incorporated | Systems and methods for correcting high order aberrations in laser refractive surgery |
DE10331792A1 (en) * | 2003-07-11 | 2005-02-17 | Medizinisches Laserzentrum Lübeck GmbH | Laser with dose rate control |
DE102005006897A1 (en) | 2005-02-15 | 2006-08-24 | Carl Zeiss Meditec Ag | Ablation program establishing method for correcting ametropia of human eye, involves establishing ablation program based on water concentration of surface to be ablated and pre-compensated set-point ablation profile |
ES2371407T3 (en) | 2005-05-02 | 2012-01-02 | Schwind Eye-Tech-Solutions Gmbh & Co. Kg | PROCEDURE FOR THE CONTROL OF A LASER FOR THE ABLATION OF A CORNEA LAYER. |
DE202005018911U1 (en) | 2005-06-02 | 2006-03-09 | Wavelight Laser Technologie Ag | Refractive ophthalmology surgery device for correcting defective vision, has two treatment stations in which patient is subjected to different treatments, and couch that is moved from one treatment station to other station |
US20070027438A1 (en) | 2005-07-26 | 2007-02-01 | Frieder Loesel | System and method for compensating a corneal dissection |
US20070282313A1 (en) | 2006-06-01 | 2007-12-06 | University Of Southern California | Method and apparatus to guide laser corneal surgery with optical measurement |
US10729585B2 (en) * | 2006-07-19 | 2020-08-04 | Alcon Inc. | Computer program for ophthalmological surgery |
DE102006036085A1 (en) * | 2006-08-02 | 2008-02-07 | Bausch & Lomb Incorporated | Method and apparatus for calculating a laser shot file for use in an excimer laser |
DE102006036086A1 (en) * | 2006-08-02 | 2008-02-07 | Bausch & Lomb Incorporated | Method and apparatus for calculating a laser shot file for use in a refractive excimer laser |
EP1886758B1 (en) * | 2006-08-07 | 2009-03-18 | WaveLight AG | Laser system for refractive surgery |
DE102008053827A1 (en) * | 2008-10-30 | 2010-05-12 | Technolas Perfect Vision Gmbh | Apparatus and method for providing a laser shot file |
-
2006
- 2006-08-02 DE DE102006036085A patent/DE102006036085A1/en not_active Ceased
-
2007
- 2007-07-27 CN CN2007800355308A patent/CN101677873B/en not_active Expired - Fee Related
- 2007-07-27 EP EP07787996A patent/EP2049053B1/en active Active
- 2007-07-27 ES ES07787996T patent/ES2337523T3/en active Active
- 2007-07-27 JP JP2009522237A patent/JP2009545349A/en not_active Withdrawn
- 2007-07-27 WO PCT/EP2007/057782 patent/WO2008015175A1/en active Application Filing
- 2007-07-27 AU AU2007280515A patent/AU2007280515B2/en not_active Ceased
- 2007-07-27 KR KR1020097004207A patent/KR101401823B1/en not_active IP Right Cessation
- 2007-07-27 DE DE602007003381T patent/DE602007003381D1/en active Active
- 2007-07-27 CA CA2660468A patent/CA2660468C/en active Active
-
2008
- 2008-02-07 US US12/375,428 patent/US20090306635A1/en not_active Abandoned
-
2013
- 2013-12-19 US US14/134,541 patent/US9636256B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6090100A (en) * | 1992-10-01 | 2000-07-18 | Chiron Technolas Gmbh Ophthalmologische Systeme | Excimer laser system for correction of vision with reduced thermal effects |
US20030023233A1 (en) * | 2001-07-30 | 2003-01-30 | Smith Michael J. | Technique for removal of material to achieve a desired shape with a laser |
Also Published As
Publication number | Publication date |
---|---|
US20150173946A1 (en) | 2015-06-25 |
ES2337523T3 (en) | 2010-04-26 |
US9636256B2 (en) | 2017-05-02 |
US20090306635A1 (en) | 2009-12-10 |
KR20090054433A (en) | 2009-05-29 |
CN101677873A (en) | 2010-03-24 |
AU2007280515A1 (en) | 2008-02-07 |
AU2007280515B2 (en) | 2013-07-04 |
WO2008015175A1 (en) | 2008-02-07 |
CA2660468A1 (en) | 2008-02-07 |
DE102006036085A1 (en) | 2008-02-07 |
EP2049053B1 (en) | 2009-11-18 |
CA2660468C (en) | 2013-03-19 |
DE602007003381D1 (en) | 2009-12-31 |
CN101677873B (en) | 2012-02-22 |
KR101401823B1 (en) | 2014-05-29 |
EP2049053A1 (en) | 2009-04-22 |
JP2009545349A (en) | 2009-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4396896B2 (en) | Threshold matrix generation method, threshold matrix generation apparatus, and recording medium | |
US9345620B2 (en) | Apparatus and method for providing a laser shot file | |
US6673062B2 (en) | Generating scanning spot locations for laser eye surgery | |
US9636256B2 (en) | Method and apparatus for calculating a laser shot file for use in an excimer laser | |
EP2046251B1 (en) | Method and apparatus for calculating a laser shot file for use in a refractive excimer laser | |
US6241356B1 (en) | Method for determining data for treating a surface | |
JP2012070269A (en) | Image processing apparatus and image processing method | |
JP4124581B2 (en) | Image pixel arrangement correction method | |
US20230263664A1 (en) | Treatment apparatus for an eye treatment, method, computer program and computer-readable medium | |
CN112863641A (en) | Radiation therapy system and offset determination method and device of radiation source thereof | |
EP3593858A1 (en) | Ripple filter unit for use in radiotherapy treatment, methods for radiotherapy treatment planning and delivery and computer program products |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: BARCLAYS BANK PLC, AS COLLATERAL AGENT, NEW YORK Free format text: SECURITY AGREEMENT;ASSIGNORS:TECHNOLAS PERFECT VISION GMBH;DR. GERHARD MANN CHEM-PHARM. FABRIK GMBH;REEL/FRAME:036400/0711 Effective date: 20150819 |
|
AS | Assignment |
Owner name: TECHNOLAS PERFECT VISION GMBH, GERMANY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:HEGELS, ERNST;REEL/FRAME:038257/0117 Effective date: 20160403 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
AS | Assignment |
Owner name: THE BANK OF NEW YORK MELLON, NEW YORK Free format text: SECURITY INTEREST;ASSIGNOR:TECHNOLAS PERFECT VISION GMBH;REEL/FRAME:043251/0910 Effective date: 20170717 |
|
AS | Assignment |
Owner name: THE BANK OF NEW YORK MELLON, AS COLLATERAL AGENT, NEW YORK Free format text: SECURITY INTEREST;ASSIGNORS:ATON PHARMA, INC.;BAUSCH & LOMB INCORPORATED;BAUSCH & LOMB PHARMA HOLDINGS CORP.;AND OTHERS;REEL/FRAME:045444/0634 Effective date: 20180213 Owner name: BARCLAYS BANK PLC, AS COLLATERAL AGENT, NEW YORK Free format text: SECURITY INTEREST;ASSIGNORS:ATON PHARMA, INC.;BAUSCH & LOMB INCORPORATED;BAUSCH & LOMB PHARMA HOLDINGS CORP.;AND OTHERS;REEL/FRAME:045444/0299 Effective date: 20180213 Owner name: THE BANK OF NEW YORK MELLON, AS COLLATERAL AGENT, Free format text: SECURITY INTEREST;ASSIGNORS:ATON PHARMA, INC.;BAUSCH & LOMB INCORPORATED;BAUSCH & LOMB PHARMA HOLDINGS CORP.;AND OTHERS;REEL/FRAME:045444/0634 Effective date: 20180213 |
|
AS | Assignment |
Owner name: THE BANK OF NEW YORK MELLON, AS COLLATERAL AGENT, Free format text: INTELLECTUAL PROPERTY SECURITY AGREEMENT;ASSIGNORS:BAUSCH HEALTH IRELAND LIMITED;BAUSCH HEALTH COMPANIES INC.;BAUSCH HEALTH, CANADA INC.;AND OTHERS;REEL/FRAME:049672/0652 Effective date: 20190701 Owner name: THE BANK OF NEW YORK MELLON, AS COLLATERAL AGENT, NEW YORK Free format text: INTELLECTUAL PROPERTY SECURITY AGREEMENT;ASSIGNORS:BAUSCH HEALTH IRELAND LIMITED;BAUSCH HEALTH COMPANIES INC.;BAUSCH HEALTH, CANADA INC.;AND OTHERS;REEL/FRAME:049672/0652 Effective date: 20190701 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 4 |
|
AS | Assignment |
Owner name: THE BANK OF NEW YORK MELLON, AS NOTES COLLATERAL AGENT, NEW YORK Free format text: SECURITY INTEREST;ASSIGNORS:BAUSCH & LOMB IRELAND LIMITED;BAUSCH HEALTH COMPANIES INC.;DR. GERHARD MANN CHEM.-PHARM. FABRIK GMBH;AND OTHERS;REEL/FRAME:057821/0800 Effective date: 20211004 |
|
AS | Assignment |
Owner name: TECHNOLAS PERFECT VISION GMBH, GERMANY Free format text: RELEASE OF SECURITY INTEREST IN SPECIFIED PATENTS (REEL/FRAME 036400/0711);ASSIGNOR:BARCLAYS BANK PLC;REEL/FRAME:061775/0826 Effective date: 20221019 Owner name: BAUSCH & LOMB INCORPORATED, NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN SPECIFIED PATENTS (REEL/FRAME 036400/0711);ASSIGNOR:BARCLAYS BANK PLC;REEL/FRAME:061775/0826 Effective date: 20221019 Owner name: LABORATOIRE CHAUVIN S.A.S., FRANCE Free format text: RELEASE OF SECURITY INTEREST IN SPECIFIED PATENTS (REEL/FRAME 045444/0299);ASSIGNOR:BARCLAYS BANK PLC;REEL/FRAME:061779/0001 Effective date: 20221019 Owner name: PF CONSUMER HEALTHCARE 1 LLC, DELAWARE Free format text: RELEASE OF SECURITY INTEREST IN SPECIFIED PATENTS (REEL/FRAME 045444/0299);ASSIGNOR:BARCLAYS BANK PLC;REEL/FRAME:061779/0001 Effective date: 20221019 Owner name: THE UNITED STATES OF AMERICA, AS REPRESENTED BY THE SECRETARY, DEPARTMENT OF HEALTH AND HUMAN SERVICES, MARYLAND Free format text: RELEASE OF SECURITY INTEREST IN SPECIFIED PATENTS (REEL/FRAME 045444/0299);ASSIGNOR:BARCLAYS BANK PLC;REEL/FRAME:061779/0001 Effective date: 20221019 Owner name: TECHNOLAS PERFECT VISION GMBH, GERMANY Free format text: RELEASE OF SECURITY INTEREST IN SPECIFIED PATENTS (REEL/FRAME 045444/0299);ASSIGNOR:BARCLAYS BANK PLC;REEL/FRAME:061779/0001 Effective date: 20221019 Owner name: BAUSCH & LOMB INCORPORATED, NEW YORK Free format text: RELEASE OF SECURITY INTEREST IN SPECIFIED PATENTS (REEL/FRAME 045444/0299);ASSIGNOR:BARCLAYS BANK PLC;REEL/FRAME:061779/0001 Effective date: 20221019 |
|
AS | Assignment |
Owner name: THE UNITED STATES OF AMERICA, AS REPRESENTED BY THE SECRETARY, DEPARTMENT OF HEALTH AND HUMAN SERVICES, MARYLAND Free format text: OMNIBUS PATENT SECURITY RELEASE AGREEMENT (REEL/FRAME 045444/0634);ASSIGNOR:THE BANK OF NEW YORK MELLON;REEL/FRAME:061872/0295 Effective date: 20221018 Owner name: TECHNOLAS PERFECT VISION GMBH, GERMANY Free format text: OMNIBUS PATENT SECURITY RELEASE AGREEMENT (REEL/FRAME 045444/0634);ASSIGNOR:THE BANK OF NEW YORK MELLON;REEL/FRAME:061872/0295 Effective date: 20221018 Owner name: LABORATOIRE CHAUVIN S.A.S., FRANCE Free format text: OMNIBUS PATENT SECURITY RELEASE AGREEMENT (REEL/FRAME 045444/0634);ASSIGNOR:THE BANK OF NEW YORK MELLON;REEL/FRAME:061872/0295 Effective date: 20221018 Owner name: BAUSCH & LOMB INCORPORATED, NEW YORK Free format text: OMNIBUS PATENT SECURITY RELEASE AGREEMENT (REEL/FRAME 045444/0634);ASSIGNOR:THE BANK OF NEW YORK MELLON;REEL/FRAME:061872/0295 Effective date: 20221018 Owner name: BAUSCH + LOMB IRELAND LIMITED, IRELAND Free format text: OMNIBUS PATENT SECURITY RELEASE AGREEMENT (REEL/FRAME 057821/0800);ASSIGNOR:THE BANK OF NEW YORK MELLON;REEL/FRAME:061884/0514 Effective date: 20221018 Owner name: TECHNOLAS PERFECT VISION GMBH, GERMANY Free format text: OMNIBUS PATENT SECURITY RELEASE AGREEMENT (REEL/FRAME 057821/0800);ASSIGNOR:THE BANK OF NEW YORK MELLON;REEL/FRAME:061884/0514 Effective date: 20221018 |