High speed topography measurement of semi-diffuse objects

• US 6,024,449 A
• Filed: 07/13/1998
• Issued: 02/15/2000
• Est. Priority Date: 07/13/1998
• Status: Expired due to Fees

• Current Assignee: Robert F Smith
• Original Assignee: Robert F Smith
• Inventors: Smith, Robert F.
• Primary Examiner(s): MANUEL, GEORGE C
• Application Number:
  US09/113,028
• Time in Patent Office: 582 Days
• Field of Search: 351/205, 351/211, 351/212, 606/12, 606/5, 606/108, 356/376
• US Class Current: 351/212
• CPC Class Codes: A61B 3/107 for determining the shape or measuring the curvature of the cornea

Abstract

Forming an image on a CCD of a grid pattern projected on a target surface and detecting the positional coordinates of the image by correlating voltage peaks in the CCD output voltage during pixel readout synchronized with the count of pulses from a master clock. Topography is calculated by knowing the equations of the rays comprising the projection pattern and the rays from the determined CCD positional coordinates and solving for the 3D coordinates formed by the intersections of the two sets of rays. Means are provided for: blocking the specular component of reflected light from the target surface while passing the diffuse component; forming the rays of the projected grid pattern so that fluctuations in the target surface position do not cause measurement error; adequate focussing of the full grid pattern on the CCD; differentiating the CCD output voltage and detecting a zero crossing for precise timing of count capture at the instant of a peak. In a preferred embodiment, the target surface is a cornea undergoing PRK and a laser projection source pulsed at a high enough rate to provide real time differential topographical display and feedback to the photoablation means. Alternate embodiments include: using the uv photoablating laser as also the projection source in conjunction with a uv sensitive CCD; an adjunct CCD type sensor in the 8 micron region for monitoring corneal sector temperatures during PRK concurrent with the topography measurement.

Claims

1. A system for rapidly measuring the topography of a semidiffuse target surface comprising:

a projection source means comprising a laser producing a pulsed plane polarized light beam having a substantially uniform intensity cross section, said light beam passing through a mask, said mask having apertures to form a grid pattern, said light beam through said grid pattern projecting sharply defined bright and dark areas on the semi-diffuse target surface;

a video camera means comprising a polarizing filter, image forming optics and a CCD (charge coupled device), said camera forming on the CCD an image of the grid pattern projected on said target surface, the polarizing filter blocking specular reflected light and passing diffuse reflected light from said image on the target surface;

a focussing means whereby said projection source and the video camera are mounted on a common platform, said platform being movable as to maintain focus and positioning of the grid pattern on the CCD;

a master clock producing low jitter constant frequency pulses;

a frequency divider receiving pulses from said master clock and delivering pulses at a reduced repetition rate to the CCD so that a single pixel readout occurs for each of the reduced repetition rate pulses, said readout resulting in an output voltage from the CCD proportional to pixel light intensity;

a master clock counter giving a count of all pulses from said master clock over a time period corresponding to a complete readout of the pixels of the CCD, said count precisely identifying the pixel location on the CCD;

an analog circuit assembly consisting of a differentiator, zero crossing detector, and gate circuit, said assembly enabling the capture of the master clock count at an instant when a voltage peak corresponding to a light intensity peak is outputted by the CCD and when said voltage peak is above a threshold value;

a digital register, into said register is stored the master clock count at the instant an enabling pulse is issued by said gate circuit;

an analog to digital converter producing a digital equivalent of the CCD output voltage, said digital equivalent being stored in said digital register simultaneous with storage of said master clock count, said converter being used to provide a feedback signal to automatically control movement of said platform;

a digital computing means, either in an application specific integrated circuit (ASIC) or computer based, performing in real time the functions of temporary storage, averaging, interpolation, error correction, topography calculation and displays.

2. The system of claim 1 where the CCD is a high efficiency back illuminated type with electronic shuttering on the order of microseconds.

3. The system of claim 1 where said laser produces pulse widths on the order of nanoseconds, and where pulsing is synchronized with the electronic shutter on the CCD.

4. The system of claim 1 where said source projection means comprises:

a concave lens to expand the beam of the laser producing a diverging beam;

an aspherical front surface convex lens with a spherical rear surface, said lens having a diameter greater than the diameter of the desired projected grid pattern, said convex lens capturing said diverging beam when said diverging beam has expanded to the diameter of the convex lens, said convex lens rapidly converging said diverging to form a rapidly converging beam;

a mask placed in close proximity to the rear surface of said convex lens, said mask apertured to form a pattern of vertical slits said pattern of vertical slits having a circular boundary producing a cross sectionally modulated converging beam;

a second concave lens placed at a distance from said mask capturing all of said modulated converging beam at a desired diameter, said second concave lens producing a slightly diverging beam, the cross section of said slightly diverging beam preserving the pattern of vertical slits sufficiently free from diffraction blurring at a maximized distance from the second concave lens.

5. The system of claim 1 where the source projection means comprises:

a concave lens to expand the beam of the laser producing a diverging beam;

a diffuser receiving said diverging beam, said diffuser producing disperse light beam;

a convex condensing lens capturing said disperse light beam, said condensing lens directing the disperse light beam through said mask;

an aspherical projection lens imaging the disperse light through said mask on said CCD.

6. The system according to claim 5 where the concave lens, the condensing lens, and the diffuser are replaced by a bundle of optical fibers, said fibers being bundled closely together at one end so that the normally narrow beam of the laser spans the bundle of fibers, each of the opposite ends of said fibers being separated and inserted into an alignment mask, said mask constraining said fibers to form a grid pattern directing each fiber towards said projection lens.

7. The system according to claim 1 where said semi-diffuse target surface is a deepithelized human cornea, said cornea reflecting diffuse light from said projected grid pattern, said cornea undergoing photoablation.

8. The system according to claim 1 where said image forming optics comprise an objective lens nominally placed 400 mm from said target surface, said objective lens being separated from the CCD by nominally 400 mm, said objective lens being biconvex with equal paracentral radii of 200 mm with a spherical rear surface and an aspherical front surface, the front surface being hyperbolic with an eccentricity of 2.12, said objective lens being 2 mm thick, 20 mm in diameter and having an index of refraction of 1.5;

a field flattening lens having a concave front surface with a paraxial radius of curvature of 3 mm and an eccentricity of 1.56, said field flattening lens having a central thickness of 1.8 mm and a flat rear surface, said field flattening lens having an index of refraction of 1.91.

9. The system of claim 1 where said focussing means comprises a digital to analog converter (DAC), said DAC receiving from said digital computing means values of errors from nominal positioning of said target surface along a set of three dimensional axes, said DAC outputting to a multiplexer a sampled set of three voltages proportional to said errors from nominal portioning, said multiplexer routing said sampled voltages to three motors, said motors driving said platform along said axes in directions to minimize said errors.

10. The system of claim 1 where said laser light source produces pulsed ultraviolet radiation of 150 to 200 nanometer wavelength, said laser being an excimer laser, said excimer laser being used to perform photoablation on a cornea, said excimer laser producing a beam of crosssectional diameter about equal to a desired projection diameter along an axis onto said cornea, said beam being diverted by a partially reflecting mirror to a second mirror to form an intensity attenuated beam, said second mirror displaced from said axis and reflecting said attenuated beam through said mask forming the grid pattern on the cornea, said laser beam being blocked by a modulator on alternate laser pulses, during said alternate pulses the grid pattern on the cornea is imaged by said video camera onto a CCD sensitive to said ultraviolet radiation wavelength.

11. In claim 10 the system where said video camera images a fluoresced projected grid pattern, said fluoresced pattern being caused by fluorescent material contained in or applied to said cornea.

12. The system of claim 1 where said video camera comprises:

an objective lens transparent to both laser and infrared radiation, said infrared radiation being of wavelength in the 8 nanometer range;

a mirror reflecting said infrared radiation but transmitting said laser beam;

a quantum well infrared photodetector (QWIP) receiving said infrared radiation, pixel readout of said QWIP being synchronized with the CCD readout, said infrared radiation emanating from the surface of a cornea undergoing photoablation, said QWIP being sensitive to a range of infrared radiation variations resulting from temperature increases induced in the cornea by said photoablation.

13. A method for performing high speed topography measurements of a semi diffuse target surface comprising the steps of;

aligning a light source to project through an apertured mask a cross sectionally modulated grid pattern;

projecting said grid pattern on said target surface;

imaging, using an objective lens, said grid pattern on a CCD with sufficient resolution to achieve a desired topography accuracy;

setting a pulse repetition rate of a master clock at a rate higher than the rate at which light intensity at successive pixels is measured;

reading out a voltage from the CCD proportional to light intensity of each pixel of the CCD, the step of pixel readout being performed at a submultiple of the pulse rate of the master clock, the step of pixel readout producing an output voltage as a function of time;

counting the pulses from the master clock to produce a count during the readout of the pixels of the CCD, resetting said count to zero upon completion of the readout of all of the pixels of the CCD;

capturing said counts precisely when peaks are detected in said output voltage, said peaks corresponding to peaks of light intensity;

storing said captured counts in memory, correlating the counts with two dimensional coordinates on the CCD corresponding to said detected peaks;

calculating three dimensional coordinate points corresponding to the topography of the target surface from said two dimensional coordinates.

14. The method of claim 13 where the step of aligning a light source also includes the steps of:

defining a projection axis from said light source through said target surface, defining an optical axis between said surface and said imaging means, defining a horizontal plane;

aligning the projection axis and the optical axis to both lie in said horizontal plane, aligning the projection axis and the optical axis to intersect on the target surface at a non-zero angle on said horizontal plane;

locating a focal point along said projection axis at a point corresponding to the center of said objective lens as measured relative to said optical axis, diverging light rays from said focal point through said mask resulting in an expanding grid pattern, projecting this pattern onto the target surface.

15. The method of claim 13 where the step of imaging includes the steps of:

converting to digital form the successive pixel intensity values readout from the CCD;

storing N most recent values of said successive pixel intensity values and after detecting a peak intensity value, storing N/2 more successive pixel values before;

forming a ratio of said peak intensity value to an average of remaining N--1 intensity values;

iterating the step of forming a ratio to produce a sequence of ratios;

forming an average of said ratios across a full row of pixels to produce a present row average;

checking said present row average with a previous row average, checking direction of motion of a focussing means;

modifying said direction of motion to increase said present row average depending on whether said previous row average is greater or less than the present row average.

16. The method of claim 13 where the steps of reading out pixel intensity, counting master clock pulses, and capturing the master clock count also involve the step of reading out all pixel rows of the CCD in the presence of a vertical slot grid pattern image with a circular boundary followed by the steps of:

selecting those rows of the CCD to produce the same approximate spacing as that of the vertical line image spacing resulting in selected rows;

adjusting the length of said vertical slots so that the number of light intensity peaks occurring across each selected row is a predefined number;

temporarily storing the count captured at intensity peaks over several consecutive rows of pixels;

checking the number of light intensity peaks actually detected in a given selected row and if fewer than said predefined number, back checking the positions in said given selected row to determine approximate positions of missing peaks, then

comparing peak positions in several previous adjacent rows with said approximate positions and if sufficiently close, extrapolating these previous positions to produce extrapolated peak position values, otherwise discarding all position values in this row containing missing peaks;

averaging the counts for each vertical line one or more rows on either side of each said selected row.

17. The method of claim 13 where the step of reading out a voltage is followed by the steps of:

differentiating with respect to time the voltage out of the CCD by means of a differentiator;

detecting a transition from a positive to a negative voltage out of the differentiator and triggering a gate circuit to produce a sharp voltage pulse if the voltage out of the CCD is above a threshold value, if so sending said pulse to cause capture of the master clock count value.

18. The method of claim 13 where the step of capturing said count includes the steps of:

dividing master clock pulse frequency by a factor producing a pulse interval, said interval corresponding to the interval between successive pixel readouts of said CCD;

setting said factor so that the number of master clock pulses between successive pixel readouts corresponds to the desired positional resolution of light intensity peaks on the CCD;

adjusting the master clock pulse frequency to select a desired CCD frame readout frequency;

compensating for a time delay between a pulse from the master clock corresponding to a fractional pixel position and the corresponding detected peak in the CCD output voltage.

19. The method of claim 13 where the steps of projecting and imaging a grid pattern also comprise the steps of:

pulsing the light source at a rate equal to a desired CCD frame rate;

selecting a short enough pulse width for said light source so that saccadic eye movement of an in-vivo semi-diffuse corneal target surface does not degrade measurement accuracy;

synchronizing CCD shutter action with the light source pulsing so as to capture all of the pulsing light while minimizing ambient light;

minimizing diffraction blurring caused by the projection system and mask.

20. A method for topography determination of a semi-diffuse surface comprising the steps of:

directing a plane polarized monochromatic pulsed through a mask to project a grid pattern on said surface;

placing a plane polarizing filter rotated at right angles to the plane of polarization of the monochromatic light source thereby blocking specular light while passing diffuse light emanating from said projected grid pattern;

imaging said pattern on a CCD of a video camera by means of an objective lens;

initiating a readout of pixels comprising the image on the CCD, said initiating also being the point in time of initiating a count of pulses from a master clock;

storing said count when a peak voltage out of the CCD corresponding to the readout of a light intensity peak occurring on a pixel or between two pixels is detected;

correlating said count with a two dimensional coordinate point on said CCD;

calculating the three dimensional topography of said surface by means of an algorithm that uses the fore knowledge of the point of origin and direction cosines of a first set of rays comprising the projected grid pattern, said algorithm using said two dimensional CCD coordinates to form a second set of rays, this second set of rays passing through a second point, said second point being the center of said objective lens, said second point and said two dimensional CCD coordinates forming a second set of direction cosines;

solving for the intersection points in three dimensional space, said intersection points being the desired topography.