ES2546076T3 - Método de calibración y corrección de la distorsión de barrido de un sistema de tomografía de coherencia óptica - Google Patents

Abstract

Procedimiento de calibración y corrección de la distorsión de barrido de cualquier sistema de tomografía de coherencia óptica, mediante el uso de patrones de referencia y la obtención de relaciones matemáticas entre las posiciones de los puntos de referencia en un patrón de referencia y las coordenadas locales de dichos puntos de referencia obtenidas mediante el citado sistema de tomografía de coherencia óptica.

Description

E12777359

25-08-2015

DESCRIPCIÓN

Método de calibración y corrección de la distorsión de barrido de un sistema de tomografía de coherencia óptica

5 Objeto de la invención

La presente invención, tal y como se expresa en el enunciado de esta memoria descriptiva se refiere a un método de calibración y corrección de la distorsión de barrido de sistemas de tomografía de coherencia óptica.

La presente invención supone una mejora de la técnica de tomografía de coherencia óptica, en cuanto a que la aplicación del método descrito en la invención permite la obtención de una topografía cuantitativa bidimensional y tridimensional de superficies a partir de imágenes de tomografía de coherencia óptica. El método es aplicable a cualquier sistema de tomografía de coherencia óptica cuyo sistema de barrido esté basado en un sistema de escáner de dos ejes.

15 La presente invención supone una mejora del estado de la técnica, al permitir la cuantificación de imágenes obtenidas mediante técnicas de tomografía de coherencia óptica, y la obtención de mapas topográficos de superficies en general, y de superficies oculares en particular, con ventajas sobre técnicas alternativas existentes.

En consecuencia, el objetivo de la presente invención es proporcionar un método para la corrección de la distorsión de cualquier sistema de tomografía de coherencia óptica, mediante el uso de patrones de referencia y la obtención de relaciones matemáticas entre las posiciones de los puntos de referencia en un patrón de referencia y las coordenadas locales de dichos puntos de referencia obtenidas mediante el citado sistema de tomografía de coherencia óptica.

25 La presente invención se refiere, en general, al campo de los sistemas de captura de imágenes, y en particular, a los sistemas de tomografía de coherencia óptica.

Antecedentes de la invención

La tomografía de coherencia óptica, OCT, (Huang, D. et al., 1991. Optical coherence tomography. Science 254:1178-1181) es una técnica interferométrica que permite obtener las diferencias de camino óptico entre superficies. Un escáner de barrido lateral de una muestra permite obtener una colección de interferogramas (Ascans) que forman una imagen de la sección transversal de la muestra (B-Scan). Un barrido en dos direcciones (x e

35 y) permite obtener una colección de B-scans y construir una imagen tridimensional de la muestra. La resolución axial de la técnica se encuentra en el orden de las micras (Povazay, B. et al. 2002. Submicrometer axial resolution optical coherence tomography. Opt. Lett. 27:1800) y viene determinada por el ancho de banda espectral de la fuente (típicamente se utilizan diodos superluminiscentes, aunque en también se utilizan láseres de femtosegundos o fuentes de barrido). Los interferogramas se pueden obtener en el dominio temporal, cambiando físicamente la longitud del brazo de referencia, o en el dominio de frecuencias: espacialmente (dominio frecuencial espacial (Fercher A. F. et al., “Measurement of Intraocular Distances by Backscattering Spectral Interferometry.” Optics Communications 1995,117:43-48)) o temporalmente (dominio frecuencial temporal (Chinn, S.R. et al. (1997). Optical coherence tomography using a frequency-tunable optical source. Opt. Lett., 22, 340-342)) codificadas mediante un espectrómetro o barriendo la frecuencia de la fuente.

45 El aumento de la velocidad de adquisición de datos en los sistemas OCT (de hasta 150.000 A-scans/s) ha permitido la captura de imágenes tridimensionales en tiempos inferiores a 1 segundo. La alta resolución axial (2-20 µm) y lateral (del orden de 100 µm) confiere a la tomografía de coherencia óptica un alto potencial para la caracterización topográfica y profilométrica de superficies y para la medida, in vivo, de la topografía corneal, entre otros.

En el estado actual de la técnica (relativo a las técnicas de tomografía de coherencia óptica, a los sistemas de topografía de superficies oculares, y a las técnicas de profilometría de superficies en general basadas en otros métodos) existe la necesidad de cuantificar los sistemas de topografía de coherencia óptica para mejorar la biometría ocular que se obtiene mediante estos sistemas, y conseguir así un nuevo y ventajoso método asociado al

55 uso de profilometría basada en tomografía de coherencia óptica. Además, existe la necesidad de un método general de calibración de sistemas de tomografía de coherencia óptica, con el fin de mejorar la información cuantitativa que se obtiene de estos sistemas. En general, el uso de la técnica de OCT como técnica topográfica viene limitada por la presencia de la distorsión de barrido asociada a la arquitectura del sistema de barrido (generalmente formada por un escáner de espejos bi-axial), y que produce distorsión de campo y astigmatismo en las imágenes. La principal contribución a dicha distorsión viene dada por la separación de los espejos en el escáner y la focal de la lente que colima el haz sobre la muestra, y en menor medida, por la planitud de los espejos y desalineamiento del haz de rotación de los espejos.

Hasta el momento, no existe ningún método general para la calibración y corrección de la distorsión de barrido

65 aplicable a cualquier sistema de tomografía de coherencia óptica, sin conocimiento previo de la configuración óptica y mecánica del sistema. La ausencia de calibración y corrección de la distorsión de barrido ha impedido la

E12777359

25-08-2015

generalización del uso cuantitativo de los sistemas de tomografía de coherencia óptica, o la correcta interpretación de datos topográficos. Uno de los objetivos principales de la presente invención es proporcionar un método para la calibración y corrección de la distorsión de barrido para la cuantificación de los datos topográficos obtenidos con cualquier sistema de Tomografía de Coherencia Óptica. La corrección de la distorsión óptica es relativamente

5 sencilla en sistemas de tomografía de coherencia óptica basados en un solo escáner, y con adquisición únicamente bidimensional de los datos. Sin embargo, en sistemas de dos escáneres, con adquisición tridimensional de datos, la distorsión es compleja, al ser no lineal, con dependencias entre las posiciones laterales y axiales, y dependiente de la configuración óptica y geométrica de cada equipo. Esta complejidad ha impedido la obtención general de datos topográficos tridimensionales cuantitativos.

Existen en el mercado varios sistemas de tomografía de coherencia óptica del segmento anterior del ojo. Estos sistemas proporcionan datos biométricos cuantitativos, generalmente en la dirección axial. Sin embargo, no se tiene constancia de la corrección de la distorsión de barrido en estos sistemas comerciales. De hecho, uno de los sistemas comerciales más extendidos (Visante, Zeiss) ha incorporado un sistema de topografía corneal de cara

15 anterior basado en anillos de Placido, a pesar de disponer de datos tridimensionales de elevación corneal. Algunos autores proporcionan configuraciones de barrido alternativas que minimizan la distorsión de barrido dependiendo de la configuración de espejos del mismo (Chinn, et al (1997). Optical coherence tomography using a frequency-tunable optical source. Opt. Lett., 22, 340-342), o en sistemas de escáneres orientados a máquinas de corte (Ireneusz Grulkowski et al “Anterior segment imaging with Spectral OCT system using a high-speed CMOS camera", Opt. Express 17, 4842-4858 (2009)). Sin embargo, estos sistemas siempre dejan distorsiones residuales que deben de ser corregidas para poder obtener las coordenadas tridimensionales de cada punto de una superficie.

Westphal et al. (Correction of geometric and refractive image distortions in optical coherence tomography applying Fermat’s principle, Opt. Express 10, 397–404 (2002)) proporciona una solución a la distorsión de barrido en sistemas

25 de OCT corneal en que el sistema de barrido es un sistema de escáner no lineal (espejos resonantes con barrido no telecéntrico), mediante la toma de imágenes axialmente alrededor de la posición axial, y aplicado únicamente a secciones bi-dimensionales de la muestra, y no imágenes tridimensionales.

Kim et al. (Automated analysis of OCT images of the crystalline lens, Proc. SPIE 7163, 716313 (2009)) utilizan un sistema telecéntrico para la adquisición de imágenes transversales (bi-dimensionales) libres de distorsión óptica. O’hara y Meyer (US7878651) proponen el uso de haces normales a la córnea para la obtención de la refracción de la córnea, sin embargo esto no produciría la pretendida corrección de la distorsión, si no más bien al contrario, dado que los haces tienen que recorrer caminos muy diferentes.

35 Ortiz et al. (Optical coherence tomography for quantitative surface topography, Appl. Opt. 48, 6708-6715 (2009)) propusieron un método para la optimización de la distorsión de barrido en un sistema propio de OCT de dominio temporal, y la corrección tri-dimensional de la distorsión de barrido residual, basado en la adquisición de imágenes axiales mediante un canal de imagen lateral confocal incorporado en el sistema de OCT. Sin embargo, este método requiere el uso de un canal confocal para la obtención de la distorsión de barrido. Dicho canal confocal no está disponible en general en instrumentos de Tomografía de Coherencia Óptica, por lo cual el método no es aplicable de manera general. Estos autores proporcionan además estimaciones teóricas de la distorsión de barrido que permiten predecir la distorsión de barrido medida experimentalmente, pero éstas exigen conocimiento preciso de la configuración óptica y geométrica del instrumento. Las estimaciones teóricas permiten obtener una configuración optimizada que permite minimizar estas distorsiones pero no eliminarlas, siendo necesario realizar el método

45 propuesto para las distorsiones residuales que permanecen en el sistema óptico de iluminación y recolección de luz.

El método de corrección de distorsión de barrido de la presente invención se puede aplicar a la obtención de profilometría de superficies en general, o a la topografía corneal en particular, mediante el uso de sistemas de tomografía de coherencia óptica.

La patente US7416300 ha descrito el uso de tomografía de coherencia óptica para la metrología de lentes y superficies, pero no alude a la corrección de la distorsión de barrido. Varias patentes US7878651, US716313 y US5491524 describen sistemas de mapeo topográfico corneal mediante tomografía de coherencia óptica, pero no aluden a la corrección de la distorsión de barrido. Generalmente, en estos estudios los mapas se obtienen a partir de

55 un conjunto de secciones transversales adquiridas para un conjunto de meridianos alrededor de un eje de rotación, centrado en el ápex corneal (de manera similar a los sistemas de Scheimpflug, o de barrido por rotación de rendija), limitando la resolución lateral en la dimensión radial.

La Solicitud de Patente US 2009/261240 divulga un modelo de prueba ("plantilla de calibración") que tiene características ópticas determinadas que se usa para generar una tabla de conversión para corregir las distorsiones que surjan, en OCT en el dominio de Fourier, a partir de las no linealidades en la relación entre el número de onda de la radiación y el desplazamiento espacial o el intervalo de tiempo de la señal de irradiación.

La Solicitud de Patente US 2010/256504 divulga un modelo de prueba, que por ejemplo comprende un número de

65 hilos finos de cruce de diámetro conocido, para calibrar una sonda de OCT. La calibración consiste en localizar cada pixel en la imagen OCT en el marco de referencia mundial de un dispositivo de medición de la posición que a su vez

E12777359

25-08-2015

forma imágenes del dispositivo de OCT así como de otros dispositivos de formación de imágenes (por ejemplo, dispositivo de formación de imágenes por ultrasonidos).

La técnica de tomografía de coherencia óptica, una vez corregida la distorsión de barrido, presenta ventajas frente a

5 la profilometría de contacto de superficies (por ejemplo Talysurf), incluyendo una mayor rapidez de adquisición de datos y la ausencia de contacto con la muestra. También presenta ventajas frente a la profilometría óptica basada en microscopía, incluyendo una mayor distancia de trabajo, mucha mayor rapidez de adquisición de datos en áreas más amplias y mayor independencia de las propiedades de reflexión especular de la muestra. La técnica de tomografía de coherencia óptica, una vez corregida la distorsión de barrido, presenta ventajas a la medida de la topografía corneal en pacientes frente a la videoqueratoscopía corneal basada en anillos de Placido, habitualmente utilizada en clínica, incluyendo una mayor resolución axial y lateral, en la dimensión radial, y la adquisición directa de datos de elevación, sin suposiciones derivadas de la presencia del “skew ray”. También presenta ventajas frente a la topografía corneal basada en Scheimpflug, incluyendo una mayor rapidez de adquisición, mayor resolución axial y lateral.

15

Descripción de la invención

Para lograr los objetivos y evitar los inconvenientes indicados anteriormente, la invención consiste en un método de calibración de cualquier tipo de sistema de tomografía de coherencia óptica.

El método de la presente invención se plantea como un protocolo de calibración de cualquier sistema de tomografía de coherencia óptica para la obtención de mapas topográficos cuantitativos a partir de las imágenes tridimensionales de tomografía de coherencia óptica de la muestra. El método permite corregir de distorsión de cualquier sistema de tomografía de coherencia óptica, independientemente de su configuración óptica y geométrica específica, a

25 diferencia de otros métodos descritos que se limitan a una configuración particular del sistema, a imágenes bidimensionales, o asumen hipótesis que exigen conocimiento previo del sistema, incluyendo la disposición de sus elementos, el alineamiento y respuesta del sistema. La presente invención aborda la medida empírica de la distorsión de barrido del sistema, determinando, mediante el uso de un patrón de referencia, la distorsión de las coordenadas espaciales del volumen de interés.

La presente invención hace uso de un patrón de referencia, que puede consistir en, aunque no está limitado a, una rejilla calibrada montada sobre un desplazador lineal axial; un cubo transparente con una rejilla tridimensional tallada, o un patrón en escalera.

35 La presente invención hace uso de imágenes de tomografía de coherencia óptica de dicho patrón de referencia, de donde se extraen las posiciones en la imagen de los puntos de referencia en el patrón, permitiendo establecer una correspondencia entre las coordenadas conocidas del patrón de referencia y las medidas resultantes del instrumento de tomografía de coherencia óptica.

En general, la correspondencia se establecerá entre puntos discretos del patrón, por ejemplo nodos en una rejilla para cada posición axial o nodos en una rejilla tridimensional y las posiciones de los puntos del patrón referencia en los datos obtenidos por el sistema de tomografía de coherencia óptica, expresadas en coordenadas locales del sistema. La correspondencia se generaliza a cualquier posición del volumen en la región calibrada mediante la interpolación –en general no lineal-de los datos entre los puntos muestreados, mediante funciones analíticas o 45 numéricas con un muestreo de puntos suficientemente denso. Por ejemplo, las líneas horizontales y verticales distorsionadas de una rejilla de calibración constituida por nodos equiespaciados unidos por líneas rectas se pueden ajustar a funciones cuádricas (coordenadas laterales) y las posiciones axiales de los nodos mediante una regresión lineal (cosenos directores). Dichas funciones analíticas permiten representar la posición de cualquier punto del volumen distorsionado en coordenadas locales. La presente invención establece una transformación matemática de coordenadas entre el volumen distorsionado, en las coordenadas locales del sistema y el volumen real del patrón de referencia, o equivalentemente, de la magnitud de la distorsión de barrido para cada punto, en general, en coordenadas angulares. La presente invención establece una transformación matemática de coordenadas entre el volumen distorsionado, en las coordenadas locales del sistema y el volumen real del patrón de referencia. La corrección de la distorsión de barrido se realiza mediante la aplicación de la transformación a todos los puntos de

55 una imagen cualquiera adquirida por el sistema de Tomografía de Coherencia Óptica.

En una realización de la invención el método se aplica directamente a imágenes de tomografía de coherencia óptica en niveles de gris. Estas imágenes se procesan mediante algoritmos digitales para la eliminación del ruido de la imagen y segmentación de superficies.

En otra realización de la invención el método se aplica a puntos de la superficie previamente segmentados.

En otra realización de la invención el método se aplica a funciones analíticas ajustadas a los bordes o superficies de la imagen.

65 En una realización de la invención el método se aplica a la imagen de cualquier superficie reflectante, obtenida

E12777359

25-08-2015

mediante cualquier sistema tomografía de coherencia óptica con escáner de barrido biaxial.

En otra realización de la invención el método se aplica a la superficie anterior corneal, obtenida mediante cualquier sistema tomografía de coherencia óptica con escáner de barrido biaxial de segmento anterior.

5 En otra realización de la invención el método se aplica a toda la imagen de segmento anterior del ojo, obtenida mediante cualquier sistema tomografía de coherencia óptica con escáner de barrido biaxial de segmento anterior.

En otra realización de la invención el método se aplica a toda la imagen de segmento anterior, obtenida mediante cualquier sistema tomografía de coherencia óptica con escáner de barrido biaxial de retina.

Así pues, la presente invención reivindica un método de calibración y corrección de la distorsión de barrido de un sistema de tomografía de coherencia óptica, que comprende las siguientes fases:

15 i) seleccionar un patrón de referencia que comprende unos puntos de referencia conocidos y que están descritos en unas coordenadas reales; ii) adquirir imágenes del patrón de referencia seleccionado en la fase i), estando situado el patrón de referencia en un espacio objeto del sistema de tomografía de coherencia óptica, mediante el sistema de tomografía de coherencia óptica; iii) identificar, en las imágenes adquiridas del patrón de referencia, unas posiciones de unos puntos de referencia descritas en unas coordenadas locales proporcionadas por el sistema de tomografía de coherencia óptica y que se corresponden a los puntos de referencia conocidos de la fase i); iv) obtener una relación matemática de distorsión que define una transformación entre las coordenadas locales proporcionadas por el sistema de tomografía de coherencia óptica y unas coordenadas reales, estando basada

25 la relación matemática en la comparación de las posiciones de los puntos de referencia conocidos en las coordenadas locales de la fase iii) y en las coordenadas reales de la fase i); y, v) corregir la distorsión mediante la aplicación de la relación matemática de distorsión obtenida en la fase iv) a los datos obtenidos por el sistema de tomografía de coherencia óptica.

En una realización preferente de la invención la relación matemática de distorsión descrita en la fase iv) comprende ser una interpolación de las posiciones de los puntos de referencia, mediante funciones seleccionadas entre funciones analíticas, funciones numéricas y una combinación de ambas.

En otra realización de la invención, el patrón de referencia de la fase i) está seleccionado entre cualquier estructura 35 bidimensional móvil con marcas localizadas en posiciones conocidas y cualquier estructura tridimensional con marcas localizadas en posiciones conocidas.

En otra realización de la invención el patrón de referencia de la fase i) está seleccionado entre:

 una rejilla de calibración tridimensional, y los puntos de referencia conocidos son los nodos de la rejilla de calibración tridimensional;  una rejilla de calibración bidimensional montada sobre una etapa lineal de desplazamiento calibrada, y los puntos de referencia conocidos son los nodos de la rejilla de calibración en distintas posiciones axiales;  un cubo con una rejilla de calibración tridimensional tallada, y los puntos de referencia conocidos son los nodos 45 de la rejilla de calibración tridimensional; y,  un patrón en escalera, y los puntos de referencia conocidos son las transiciones bruscas en profundidad entre escalones.

En otra realización de la invención, los puntos de referencia descritos mediante las coordenadas locales de la fase iii) se unen mediante líneas, estando la líneas definidas por funciones analíticas. Para llevar a cabo la calibración del sistema, los puntos identificados en la fase iii) del método se unen con líneas. Como debido a la distorsión de barrido las líneas que los unen, en cualquiera de los tres ejes del espacio, no son rectas sino curvas, para definir dichas rectas, se emplean funciones analíticas convencionales. Por ejemplo en una realización particular, las funciones analíticas corresponderían con funciones parabólicas ya que dichas parábolas se asemejarían a las curvas descritas

55 por las líneas que unirían los puntos identificados en la fase iii) del método.

Además, se ha previsto el uso del método descrito en la presente memoria para la obtención de datos calibrados y corregidos de distorsión de barrido con el sistema de tomografía de coherencia óptica seleccionados entre:

-datos de secciones bidimensionales; -datos de volúmenes tridimensionales; -datos de topografías corneales; -datos de topografías de la retina; -datos de superficies internas de un ojo en combinación con una compensación de una distorsión óptica y de

65 refracción; -datos de imagen de segmento anterior del ojo;

E12777359

25-08-2015

-datos de imagen de las capas de la retina; -datos de señales obtenidas en un fotodetector del sistema de tomografía de coherencia óptica; -datos de imágenes y volúmenes en intensidad; -datos de mapas de puntos correspondientes a unos bordes previamente extraídos de una imagen de

5 tomografía de coherencia óptica; y, -datos de superficies ajustadas a unos bordes previamente extraídos de una imagen de tomografía de coherencia óptica.

No obstante existen otros muchos usos posibles de la presente invención en diversos campos de la ciencia no especificados en la presente memoria descriptiva.

Obsérvese que la distorsión de barrido es un tipo de distorsión geométrica. Además la presente invención permite el que una vez corregida dicha distorsión de barrido mediante el proceso descrito, se puedan abordar otro tipo de correcciones como por ejemplo la corrección de la distorsión óptica, siendo dicha distorsión óptica la producida al ver

15 una superficie a través de otra superficie.

Breve descripción de los dibujos

Figura 1a.-Muestra una imagen integrada de una rejilla como patrón de referencia. Figura 1b.-Muestra la rejilla de la figura con los nodos marcados. Figura 1c.-Muestra la rejilla de la figura 1 con los nodos y la estimación de las coordenadas laterales reales marcada. Figura 2a.-Muestra la diferencia entre la topografía real (medida mediante profilometría) de una superficie esférica y la topografía basada en OCT sin corrección de la distorsión óptica.

25 Figura 2b.-Muestra la diferencia entre la topografía real (medida mediante profilometría) de una superficie esférica y la topografía basada en OCT tras la corrección de la distorsión geométrica. Figura 3a.-Muestra un mapa topográfico de la cara anterior de la córnea de un paciente obtenido mediante OCT antes de aplicar la corrección de distorsión de barrido. Figura 3b.-Muestra un mapa topográfico de la cara anterior de la córnea de un paciente obtenido mediante OCT después de aplicar la corrección de distorsión de barrido.

Descripción de varios ejemplos de realización de la invención

Seguidamente se realiza, con carácter ilustrativo y no limitativo, una descripción de varios ejemplos de realización de 35 la invención, haciendo referencia a la numeración adoptada en las figuras presentadas anteriormente.

Así pues a modo de ejemplo de realización, se describe a continuación un método preferido de calibración de un sistema de tomografía de coherencia óptica y su aplicación a la medida de una superficie de plástico y una córnea de un paciente.

El método se ha aplicado a las imágenes de OCT obtenidas mediante un Tomógrafo de Coherencia Óptica espectral de cámara anterior no comercial en un laboratorio de investigación. El dispositivo empleado incluye un sistema de escáneres x-y (8 mm de apertura y separación entre centros de los espejos de 13.6 mm) y una lente colimadorafocalizadora de 75 mm. La fuente de iluminación es un diodo superluminiscente (840 nm, 50 nm ancho espectral).

45 La línea de retardo está constituida por un espectrómetro con una red de difracción y una cámara CMOS. La velocidad de adquisición es de 25.000 A-scans (interferogramas por segundo).

El proceso descrito en la invención se aplicó siguiendo las siguientes etapas:

(1)

Fabricación de un patrón de referencia, consistente en una superficie opaca plana reflectante con interlineado conocido impreso o grabado, y en particular en una malla o rejilla de calibración milimetrada impresa en papel blanco con tinta negra;

(2)

Desplazamiento manual o motorizado con graduación, sobre la que se coloca dicha rejilla de calibración;

55 (3) Adquisición de volúmenes en 3D de la rejilla de calibración colocada en la posición de la muestra (espacio objeto), y en distintas posiciones axiales alrededor de dicha posición, en un rango de 7 mm y pasos de 0.5 mm. Las Figura 1a, 1b y 1c muestran las imágenes de la rejilla utilizada para la calibración de la distorsión de barrido en los presentes ejemplos de realización. La figura 1a muestra la rejilla de calibración empleada, la figura 1b muestra esa misma rejilla con los nodos conocidos identificados por puntos blancos. La figura 1c muestra la rejilla con los nodos identificados y el desplazamiento de dichos nodos debido a la distorsión de barrido. Como se ve de estas figuras la distorsión de barrido ha producido un desplazamiento de los nodos, estando dicho desplazamiento representado por las líneas que parten de las posiciones originales de los nodos hasta la posición en la que quedan tras la toma de imágenes en el espacio objeto del sistema de tomografía de coherencia óptica;

65 (4) Obtención de imágenes bidimensionales integradas de la rejilla de calibración, para cada posición axial, a partir de la suma de la señal de cada A-scan;

E12777359

25-08-2015

(5)

Eliminación del ruido de barrido de la imagen mediante filtros digitales;

(6)

Extracción de los bordes de las líneas de la rejilla mediante métodos de análisis de imagen, en particular una transformada de Hough;

(7)

Ajuste a funciones cuádricas de las líneas horizontales y verticales extraídas;

5 (8) Obtención de los puntos de intersección de las curvas horizontales y verticales y etiquetado de dichos nodos, para cada imagen de rejilla adquirida en una posición axial diferente. Por lo tanto, para cada nodo se obtiene un conjunto de puntos tridimensionales expresados en las coordenadas locales del dispositivo, como por ejemplo pixeles en las imágenes tridimensionales;

(9) Obtención de un factor de calibración entre coordenadas locales del dispositivo (calculadas según el paso 8)

10 y coordenadas euclídeas. La coordinada axial euclídea se obtiene a partir de las posiciones axiales de los desplazamientos lineales;

(10) Obtención de funciones analíticas que representan las posiciones de cualquier punto de la imagen de la rejilla de calibración mediante interpolación bicúbica entre las posiciones estimadas de los nodos (coordenadas laterales) y mediante una regresión lineal de las posiciones axiales de los nodos (cosenos directores);

15 (11) Obtención de las funciones de transformación entre las coordenadas reales y las coordenadas de la imagen,

o equivalentemente, de la magnitud de distorsión de barrido obtenida para cada punto, en coordenadas angulares; y,

(12)

Corrección de la distorsión de un punto genérico que se lleva a cabo mediante la aplicación de la función de transformación de coordenadas a cada punto de la imagen de una superficie obtenida mediante el sistema de

20 tomografía de coherencia óptica. La magnitud de distorsión de barrido para cada punto se sustrae (en coordenadas angulares) de la diferencia de camino óptico obtenida como señal de OCT para cada punto de la superficie detectada para obtener la posición actual de la superficie sin distorsión.

Para el sistema de OCT usado en el ejemplo, la distorsión de barrido medida es de 24 píxeles (86 m) en la

25 dirección horizontal y 7 píxeles (24 m) en la dirección vertical para un rango angular del escáner entre -7,5 y 7,5 grados.

En la aplicación del ejemplo descrito, el método se aplicó sobre datos tridimensionales de una superficie esférica de PMMA (Polimetilmetacrilato) y de una cornea de un paciente. Se utilizaron rutinas de procesado de imagen para la

30 eliminación de ruido y segmentación de las superficies. El método se utiliza para cada punto de la superficie detectado en la imagen tridimensional. Las imágenes se adquirieron sobre una zona de 10 x 10 mm, con una densidad de 200 x 200 A-scans para el caso de la superficie esférica d PMMA, y de 10 x 12 mm, con una densidad de 120x50 A-scans para el caso de la córnea del paciente.

35 La superficie de PMMA consistía en una superficie asférica tratada con un láser de cirugía refractiva (patrón de ablación miópico) lo que alteró la forma de la superficie incrementando su asfericidad corneal. Como referencia, la topografía de la superficie se evaluó con un profilómetro de no contacto basado en microscopía confocal (PLµ, Sensofar).

40 El método descrito en la invención, con los parámetros detallados en el ejemplo de aplicación, se aplicó sobre la imagen tridimensional de OCT adquirida sobre esta superficie. Los datos topográficos (del profilómetro y de OCT, datos brutos, y tras calibración/corrección de la distorsión de barrido) se ajustaron a funciones bicónicas (caracterizadas por el radio de curvatura y la asfericidad) y a polinomios de Zernike de orden 8. La discrepancia en el radio de curvatura ajustado a la superficie de OCT sin calibrar ni corregir con respecto al ajuste del perfil

45 profilométrico de no contacto fue de 4,6%, mientras que la discrepancia tras la calibración fue de 1,6%. La discrepancia en la asfericidad disminuyó de 130% al 5%. La Figura 2a muestra el mapa de diferencia entre la superficie y el ajuste del perfil profilométrico antes de la calibración y la figura 2b muestra el mapa de diferencia entre la superficie y el ajuste del perfil profilométrico después de la aplicación del método de calibración, objeto de la presente invención.

50 Para el otro ejemplo de realización en el que se aplicaba el método objeto de la invención sobre la córnea de un ser humano, el paciente medido era un sujeto joven normal. El sujeto fija la vista en un estímulo que permite alinear el eje pupilar keratométrico con el eje óptico del instrumento. El tiempo de adquisición de la imagen de OCT es de 0,72 segundos. El método descrito en la invención, con los parámetros detallados en el ejemplo de aplicación, se

55 aplicaron sobre la imagen tridimensional de OCT adquirida de la cara anterior corneal. Los datos topográficos se ajustaron a superficies bicónicas y a polinomios de Zernike de orden 8. La corrección de la distorsión de barrido disminuyó la asimetría de la superficie en 5,7 % (radio de curvatura) y 9,5 % (asfericidad). El radio de curvatura apical de la córnea antes y después de la corrección fue 7,38 y 7,59 mm respectivamente, y la asfericidad -0,38 y 0,42, respectivamente. La Figura 3a muestra el mapa topográfico de la cara anterior de la córnea (ajuste a

60 polinomios de Zernike) de un paciente obtenido mediante OCT antes de aplicar la corrección de la distorsión de barrido y la figura 3b muestra dicho mapa topográfico de la cara anterior de la córnea tras aplicar la corrección de la distorsión de barrido, en ambos casos, previa sustracción de la mejor esfera.

Claims (6)

1. REIVINDICACIONES

   1. Método de calibración y corrección de la distorsión de barrido de un sistema de tomografía de coherencia óptica, que comprende las siguientes fases:

   5 i) seleccionar un patrón de referencia que comprende unos puntos de referencia conocidos y que están descritos en unas coordenadas reales; ii) adquirir imágenes del patrón de referencia seleccionado en la fase i), estando situado el patrón de referencia en un espacio objeto del sistema de tomografía de coherencia óptica, mediante el sistema de tomografía de coherencia óptica; iii) identificar, en las imágenes adquiridas del patrón de referencia, unas posiciones de unos puntos de referencia descritas en unas coordenadas locales proporcionadas por el sistema de tomografía de coherencia óptica y que se corresponden a los puntos de referencia conocidos de la fase i); iv) obtener una relación matemática de distorsión que define una transformación entre las coordenadas locales

   15 proporcionadas por el sistema de tomografía de coherencia óptica y unas coordenadas reales, estando basada la relación matemática en la comparación de las posiciones de los puntos de referencia conocidos en las coordenadas locales de la fase iii) y en las coordenadas reales de la fase i); y, v) corregir la distorsión mediante la aplicación de la relación matemática de distorsión obtenida en la fase iv) a los datos obtenidos por el sistema de tomografía de coherencia óptica.
2. 2. Método de calibración y corrección de la distorsión de barrido de un sistema de tomografía de coherencia óptica, según la reivindicación 1, caracterizado por que la relación matemática de distorsión de la fase iv) es una interpolación de las posiciones de los puntos de referencia, mediante funciones seleccionadas entre funciones analíticas, funciones numéricas y una combinación de ambas.

   25

3. 3.

   Método de calibración y corrección de la distorsión de barrido de un sistema de tomografía de coherencia óptica, según la reivindicación 1, caracterizado por que el patrón de referencia de la fase i) está seleccionado entre cualquier estructura bidimensional móvil con marcas localizadas en posiciones conocidas y cualquier estructura tridimensional con marcas localizadas en posiciones conocidas.

4. 4.

   Método de calibración y corrección de la distorsión de barrido de un sistema de tomografía de coherencia óptica, según la reivindicación 1, caracterizado por que el patrón de referencia de la fase i) está seleccionado entre:

    una rejilla de calibración tridimensional, y los puntos de referencia conocidos son los nodos de la rejilla de

   35 calibración tridimensional;  una rejilla de calibración bidimensional montada sobre un elemento lineal de desplazamiento calibrado, y los puntos de referencia conocidos son los nodos de la rejilla de calibración en distintas posiciones axiales;  un cubo con una rejilla de calibración tridimensional tallada, y los puntos de referencia conocidos son los nodos de la rejilla de calibración tridimensional; y,  un patrón en escalera y los puntos de referencia conocidos son las transiciones bruscas en profundidad entre escalones.
5. 5. Método de calibración y corrección de la distorsión de barrido de un sistema de tomografía de coherencia óptica,

   según la reivindicación 1, caracterizado por que los puntos de referencia descritos mediante las coordenadas 45 locales de la fase iii) se unen mediante líneas, estando la líneas definidas por funciones analíticas.
6. 6. Uso del método descrito en las reivindicaciones 1 a 5 para la obtención de datos calibrados y corregidos de distorsión de barrido con el sistema de tomografía de coherencia óptica seleccionados entre:

   -datos de secciones bidimensionales; -datos de volúmenes tridimensionales; -datos de topografías corneales; -datos de topografías de la retina; -datos de superficies internas de un ojo en combinación con una compensación de una distorsión óptica y de

   55 refracción; -datos de imagen de segmento anterior del ojo; -datos de imagen de las capas de la retina; -datos de señales obtenidas en un fotodetector del sistema de tomografía de coherencia óptica; -datos de imágenes y volúmenes en intensidad; -datos de mapas de puntos correspondientes a unos bordes previamente extraídos de una imagen de tomografía de coherencia óptica; y, -datos de superficies ajustadas a unos bordes previamente extraídos de una imagen de tomografía de coherencia óptica.

   65

   8