WO2005015495A2 - Online wavefront measurement and display - Google Patents

Online wavefront measurement and display Download PDF

Info

Publication number
WO2005015495A2
WO2005015495A2 PCT/EP2004/008205 EP2004008205W WO2005015495A2 WO 2005015495 A2 WO2005015495 A2 WO 2005015495A2 EP 2004008205 W EP2004008205 W EP 2004008205W WO 2005015495 A2 WO2005015495 A2 WO 2005015495A2
Authority
WO
WIPO (PCT)
Prior art keywords
centroid
image
algorithm
wavefront
centroids
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.)
Ceased
Application number
PCT/EP2004/008205
Other languages
English (en)
French (fr)
Other versions
WO2005015495A3 (en
Inventor
Hans-Joachim Polland
Gerhard Youssefi
Stefan Franzke
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Bausch and Lomb Inc
Original Assignee
Bausch and Lomb Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Bausch and Lomb Inc filed Critical Bausch and Lomb Inc
Priority to AU2004263957A priority Critical patent/AU2004263957B2/en
Priority to EP04763406.8A priority patent/EP1658587B1/en
Priority to ES04763406.8T priority patent/ES2689935T3/es
Priority to CA002539944A priority patent/CA2539944A1/en
Priority to US10/565,703 priority patent/US7708411B2/en
Priority to JP2006520794A priority patent/JP2006528499A/ja
Publication of WO2005015495A2 publication Critical patent/WO2005015495A2/en
Publication of WO2005015495A3 publication Critical patent/WO2005015495A3/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/0016Operational features thereof
    • A61B3/0025Operational features thereof characterised by electronic signal processing, e.g. eye models
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/66Analysis of geometric attributes of image moments or centre of gravity
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30041Eye; Retina; Ophthalmic

Definitions

  • online aberrometry information referred to herein as "online” aberrometry.
  • Hartman-Shack wavefront sensor typically includes a microlens array that
  • the image produced by the microlens array comprises an
  • centroid This aspect of the wavefront analysis process is known as centroid
  • Hartmann-Shack wavefront sensors and other well-known types such as
  • Tscherning typically measure single images of centroids or,
  • a single wavefront image requires about 400Kb of
  • aberration measurements e.g., sphere
  • Figure 1 is a flow chart schematic of a fast centroid detection
  • Figure 2 is a reproduced photographic of a Hartmann-Shack
  • Figure 3 is a drawing representation of a wavefront (centroid) image
  • Figure 4 is a flow chart schematic of a sorting algorithm according to
  • Figure 5 is another drawing representation of a wavefront (centroid)
  • Figure 6 is a schematic illustration representing another algorithm
  • Figure 7 is an illustration representing another algorithm process
  • Figure 8 is a reproduced photographic of part of a Hartmann-Shack
  • Figure 9 is a flow chart schematic directed to another embodiment of
  • Figure 10 is a block diagram of an apparatus embodiment of the invention.
  • the image may comprise, for illustration purposes only, centroids in a
  • centroids 32 in a CCD image can be located and sorted in
  • taken from a CCD camera usually consist of an array of pixels in which
  • every pixel is assigned a number proportional to the amount of charge
  • This number is referred to as the signal of the pixel.
  • the image is compressed from
  • next pixel is set to the average of the next (second) square, and so on,
  • the compressed image is then divided into square
  • one tile is a square of 64x64 pixels, but other
  • a signal-to-noise ratio of two was improved to a signal-to-noise ratio of 10 by
  • a maximum is defined as the highest signal in the
  • the maximum is determined by a scan through the
  • the first entry (highest signal value) is defined as the first rough structure
  • the pre-set is defined as rough structure points.
  • condition is that the position of the particular entry is farther away from all
  • the distance is 17 pixels.
  • step 108 can be repeated as shown at block 110, setting a new
  • the ultimate centroid positions are determined. Since, the
  • each square is smaller than 2x the
  • this value is larger than the centroid itself. In the exemplary embodiment this value is
  • step 112 can be repeated, for
  • structure point can also be assigned a quality factor depending on how much
  • this distance is five pixels.
  • an aspect 200 of the embodiment as illustrated in Figure 4 is directed to the
  • the desired sorting configuration is selected.
  • the configuration is a square grid based upon the
  • the centroid with the smallest nj value is assigned to
  • Row 1 and is stored in memory as the last centroid of Row 1.
  • the last centroids of the existing rows are stored in
  • a region 610 is defined that, in an exemplary embodiment, comprises an area to the right of the last centroid
  • search area 610 for the next centroid Any shape suited for detecting other
  • parameters include maximum angle 702, minimum distance 704, maximum
  • centroid of that row If no, than that centroid is assigned the last centroid of
  • Steps 214-216 are now repeated for all centroids. In this
  • This step facilitates marking of the top row as Row 1,
  • optional step 222 involves merging these rows.
  • the criteria for merging rows j and k is: If y j -y ⁇
  • f is a variable parameter in the range between about 0.1-0.7, that is set by
  • a is the mean distance between rows (see above);
  • P k ,&st is the x value of the first (left-most) centroid of the k row
  • P k jast is the x value of the last (right-most) centroid of the k row.
  • the rows are merged if they are much closer in the, y-position
  • row j is either completely to the
  • the process for sorting the columns begins at step 224 where the list
  • centroid with the smallest nj is assigned to Column 1 and is stored in memory as the last centroid of Column 1.
  • the last centroid of Column 1 is assigned to Column 1 and is stored in memory as the last centroid of Column 1.
  • centroids of an existing column are always stored in memory during step
  • a region is defined that, in the exemplary embodiment,
  • nj value is selected and that centroid is checked, with respect to all existing
  • centroid is assigned as the last centroid of that column. If no, than
  • This step facilitates marking ofthe left-most column as Column 1, the next
  • 236 involves merging these columns. This is accomplished by the following sub-steps: From the mean average x-position for each column
  • f is a variable parameter in the range between about 0.1-0,7, that is set by
  • a is the mean distance between columns
  • P k ,- Brst is the y value ofthe first (top-most) centroid ofthe k column
  • column j is either
  • Pupil coordinate location (optional): ⁇ 6-8ms;
  • centroid position data results in data storage requirements of only about 700Kb of memory for a 20 second measurement at 25Hz, which will yield
  • the method comprises at step 902 acquiring a plurality of wavefront
  • each ofthe images includes
  • the images are acquired at a rate
  • the images are
  • the display will be limited to second order aberrations (sphere
  • a pupil diameter value can be selected prior to measurement allowing
  • any Zernike order e.g., coma, spherical aberration, higher-orders
  • diameters is between about 2mm to 10mm.
  • the first 125 images might be obtained for pupil diameter, D
  • step 910 an average value of a selected Zernike order can be calculated
  • blinking periods can be determined which
  • step 914 a In a related aspect according to this embodiment, at step 914 a
  • sequence of pupil images corresponding to the wavefront images can also be
  • the apparatus includes an illumination component 1010 that delivers a
  • a detector 1040 adapted to acquire the
  • a display component 1060 operatively connected
  • pupil camera 1080 and pupilometer 1090 is also shown, where components 1015 are beam splitters or optical equivalents thereof.

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Public Health (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Biomedical Technology (AREA)
  • Ophthalmology & Optometry (AREA)
  • Geometry (AREA)
  • Signal Processing (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Image Analysis (AREA)
  • Eye Examination Apparatus (AREA)
  • Testing Of Optical Devices Or Fibers (AREA)
  • Image Processing (AREA)
  • Length Measuring Devices By Optical Means (AREA)
PCT/EP2004/008205 2003-07-24 2004-07-22 Online wavefront measurement and display Ceased WO2005015495A2 (en)

Priority Applications (6)

Application Number Priority Date Filing Date Title
AU2004263957A AU2004263957B2 (en) 2003-07-24 2004-07-22 Online wavefront measurement and display
EP04763406.8A EP1658587B1 (en) 2003-07-24 2004-07-22 Online wavefront measurement and display
ES04763406.8T ES2689935T3 (es) 2003-07-24 2004-07-22 Medición y visualización de frente de onda en línea
CA002539944A CA2539944A1 (en) 2003-07-24 2004-07-22 Online wavefront measurement and display
US10/565,703 US7708411B2 (en) 2003-07-24 2004-07-22 Online wavefront measurement and display
JP2006520794A JP2006528499A (ja) 2003-07-24 2004-07-22 オンライン波面測定および表示

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE10333813.6 2003-07-24
DE10333813A DE10333813A1 (de) 2003-07-24 2003-07-24 Online-Wellenfrontmessung und Anzeige

Publications (2)

Publication Number Publication Date
WO2005015495A2 true WO2005015495A2 (en) 2005-02-17
WO2005015495A3 WO2005015495A3 (en) 2005-12-22

Family

ID=34071891

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2004/008205 Ceased WO2005015495A2 (en) 2003-07-24 2004-07-22 Online wavefront measurement and display

Country Status (11)

Country Link
US (1) US7708411B2 (enExample)
EP (1) EP1658587B1 (enExample)
JP (1) JP2006528499A (enExample)
KR (1) KR20060080578A (enExample)
CN (1) CN1856803A (enExample)
AU (1) AU2004263957B2 (enExample)
CA (1) CA2539944A1 (enExample)
DE (1) DE10333813A1 (enExample)
ES (1) ES2689935T3 (enExample)
SG (1) SG145707A1 (enExample)
WO (1) WO2005015495A2 (enExample)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7771053B2 (en) 2003-07-24 2010-08-10 Bausch & Lomb Incorporated Method and apparatus for online contact lens evaluation
CN104103078A (zh) * 2014-07-31 2014-10-15 中国航天科工集团第三研究院第八三五七研究所 一种螺旋对称中心检测方法

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7789513B2 (en) 2007-04-20 2010-09-07 Alcon Refractivehorizons, Inc. Adaptive wavefront modulation system and method for refractive laser surgery
CN101055223B (zh) * 2007-04-26 2010-12-08 中国科学院光电技术研究所 哈特曼波前传感器质心测量精度优化方法
US8740381B2 (en) * 2007-06-27 2014-06-03 Bausch & Lomb Incorporated Method and apparatus for extrapolating diagnostic data
DE102008014294A1 (de) * 2008-03-14 2009-09-17 Bausch & Lomb Inc. Schneller Algorithmus für Wellenfrontdatenstrom
US8254724B2 (en) * 2008-11-06 2012-08-28 Bausch & Lomb Incorporated Method and apparatus for making and processing aberration measurements
US9504376B2 (en) * 2009-12-22 2016-11-29 Amo Wavefront Sciences, Llc Optical diagnosis using measurement sequence
CN104955148B (zh) * 2014-12-09 2019-03-19 文春明 一种利用电磁波对称传播特性的无线传感网络定位方法
US9782064B1 (en) 2016-04-08 2017-10-10 Clarity Medical Systems, Inc. Obtaining and displaying histogram and/or confidence of intra-operative refraction and/or IOL power recommendation
CN106530278B (zh) * 2016-10-14 2020-01-07 中国科学院光电技术研究所 用于点源哈特曼波前探测器的点光斑检测与背景噪声特征估计方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4222395A1 (de) 1992-07-08 1994-01-13 Amtech Ges Fuer Angewandte Mic Vorrichtung und Verfahren zur Messung der Augenrefraktion
US5777719A (en) 1996-12-23 1998-07-07 University Of Rochester Method and apparatus for improving vision and the resolution of retinal images
FR2788597B1 (fr) 1999-01-15 2001-02-23 Imagine Optic Sarl Procede et dispositif d'analyse de front d'onde a grande dynamique
US6234631B1 (en) 2000-03-09 2001-05-22 Lasersight Technologies, Inc. Combination advanced corneal topography/wave front aberration measurement
DE10154194A1 (de) 2001-11-07 2003-05-22 Asclepion Meditec Ag Verfahren und Vorrichtung zur Messung des Dynamischen Verhaltens eines optischen Systems

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
DROSTE D ET AL: "AN ASIC FOR HARTMANN-SHACK WAVEFRONT DETECTION" IEEE JOURNAL OF SOLID-STATE CIRCUITS, IEEE INC. NEW YORK, US, vol. 37, no. 2, February 2002 (2002-02), pages 173-182, XP001220695 ISSN: 0018-9200 *
HOFER H ET AL: "DYNAMICS OF THE EYE'S WAVE ABERRATION" JOURNAL OF THE OPTICAL SOCIETY OF AMERICA - A, OPTICAL SOCIETY OF AMERICA, WASHINGTON, US, vol. 18, no. 3, March 2001 (2001-03), pages 497-506, XP002338276 ISSN: 1084-7529 *
I. MIRO, N. LOPEZ-GIL, AND P. ARTAL: "Pupil meter and tracking system based in a fast image processing" OPHTHALMIC TECHNOLOGIES IX, 3591, 497-506, 1999, XP008050380 *
LIANG J ET AL: "OBJECTIVE MEASUREMENT OF WAVE ABERRATIONS OF THE HUMAN EYE WITH THE USE OF A HARTMANN-SHACK WAVE-FRONT SENSOR" JOURNAL OF THE OPTICAL SOCIETY OF AMERICA. A, OPTICS AND IMAGE SCIENCE, vol. 11, no. 7, July 1994 (1994-07), pages 1949-1957, XP001182893 ISSN: 0740-3232 *
NIRMAIER T ET AL: "Hartmann-Shack sensor ASIC's for real-time adaptive optics in biomedical physics" 6TH WORLD MULTICONFERENCE ON SYSTEMICS, CYBERNETICS AND INFORMATICS. PROCEEDINGS INT. INST. INF. & SYST ORLANDO, FL, USA, [Online] vol. 13, 2002, pages 280-284 vol.1, XP002338274 ISBN: 980-07-8150-1 Retrieved from the Internet: URL:http://wwwasic.kip.uni-heidelberg.de/wavefront/iii_final.pdf> [retrieved on 2005-07-26] *
PRIETO P M ET AL: "ANALYSIS OF THE PERFORMANCE OF THE HARTMANN-SHACK SENSOR IN THE HUMAN EYE" JOURNAL OF THE OPTICAL SOCIETY OF AMERICA - A, OPTICAL SOCIETY OF AMERICA, WASHINGTON, US, vol. 17, no. 8, August 2000 (2000-08), pages 1388-1398, XP001041251 ISSN: 1084-7529 *
SOPHIA I. PANAGOPOULOU; IOANNIS G. PALLIKARIS: "Wavefront Customized Ablations With the WASCA Asclepion Workstation" JOURNAL OF REFRACTIVE SURGERY, [Online] vol. 17, September 2001 (2001-09), XP002338275 Retrieved from the Internet: URL:http://www.carlomasci.it/biblio/secondo%20congresso/wasca.pdf> [retrieved on 2005-07-26] *
SUNG-HOON BAIK ET AL: "New centroid detection algorithm for the Shack-Hartmann wavefront sensor" PROCEEDINGS OF THE SPIE - THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING SPIE-INT. SOC. OPT. ENG USA, vol. 4926, 2002, pages 251-260, XP002348357 ISSN: 0277-786X *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7771053B2 (en) 2003-07-24 2010-08-10 Bausch & Lomb Incorporated Method and apparatus for online contact lens evaluation
CN104103078A (zh) * 2014-07-31 2014-10-15 中国航天科工集团第三研究院第八三五七研究所 一种螺旋对称中心检测方法

Also Published As

Publication number Publication date
JP2006528499A (ja) 2006-12-21
WO2005015495A3 (en) 2005-12-22
EP1658587A2 (en) 2006-05-24
AU2004263957B2 (en) 2010-05-13
US20070008491A1 (en) 2007-01-11
AU2004263957A1 (en) 2005-02-17
CN1856803A (zh) 2006-11-01
ES2689935T3 (es) 2018-11-16
KR20060080578A (ko) 2006-07-10
CA2539944A1 (en) 2005-02-17
DE10333813A1 (de) 2005-02-17
EP1658587B1 (en) 2018-09-19
SG145707A1 (en) 2008-09-29
US7708411B2 (en) 2010-05-04

Similar Documents

Publication Publication Date Title
CA2533068C (en) Method and apparatus for online contact lens evaluation
AU2004263957B2 (en) Online wavefront measurement and display
AU2007324081B2 (en) Focus assist system and method
US20160012575A1 (en) Image processing apparatus and image processing method
Nagin et al. The reproducibility of computerized boundary analysis for measuring optic disc pallor in the normal optic disc
CN114727755B (zh) 用于评估泪膜的稳定性的方法
CN113470015A (zh) 一种基于图像处理的水体晃动检测分析方法及系统
CN117617889A (zh) 一种基于遮盖测试的斜视度数检测系统
JP2002000567A (ja) 瞳孔中心位置計測方法及び視点位置検出方法
JP2006528499A5 (enExample)
JP2002102172A (ja) 視線検出方法及びシステム
CN118161124A (zh) 基于自适应非线性滤波器的双目瞳孔直径测量方法及系统
KR20200068062A (ko) 열화상이미지의 해상도 보정장치
CN114093018A (zh) 一种基于瞳孔定位的视力筛查设备及系统
CN111401181A (zh) 一种利于学生视力的防控筛查系统
WO2014179728A1 (en) Methods and apparatuses for determning contact lens intolerance in contact lens wearer patients based on dry eye tear film characteristic analysis and dry eye symptoms
JP2022548111A (ja) 涙液層破壊を検出するための装置および方法
JP7631780B2 (ja) 眼球運動計測装置、眼球運動計測方法及びプログラム
JP4712676B2 (ja) 画像解析装置、画像解析方法及びコンピュータプログラム
CN121033159A (zh) 一种铝电解槽测量的图像处理方法、装置及测量设备
CN117422758A (zh) 虹膜直径的获取方法、装置、电子设备及存储介质
JP2018504198A (ja) 角膜内皮の細胞の形態計測的解析のための方法および装置

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200480027469.9

Country of ref document: CN

AK Designated states

Kind code of ref document: A2

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A2

Designated state(s): BW GH GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

DPEN Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed from 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2004763406

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2006520794

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 2539944

Country of ref document: CA

WWE Wipo information: entry into national phase

Ref document number: 2004263957

Country of ref document: AU

ENP Entry into the national phase

Ref document number: 2004263957

Country of ref document: AU

Date of ref document: 20040722

Kind code of ref document: A

WWP Wipo information: published in national office

Ref document number: 2004263957

Country of ref document: AU

WWP Wipo information: published in national office

Ref document number: 2004763406

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2007008491

Country of ref document: US

Ref document number: 10565703

Country of ref document: US

WWP Wipo information: published in national office

Ref document number: 10565703

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 1020067001479

Country of ref document: KR