I. INTRODUCCIÓN
Dada la importancia en el diagnóstico y tratamiento de diferentes dolencias, las técnicas de formación de imágenes médicas se encuentran en constante mejora y evolución, el espectro electromagnético brinda una serie de herramientas ópticas que pueden ser utilizadas a tal fin. Desde hace dos décadas, se ha desarrollo el estudio de una técnica alternativa y complementaria que emplea luz en el infrarrojo cercano (NIR), que ha captado especial interés por su característica de inocuidad y no invasividad.1-5
Esta luz NIR está comprendida entre los 600 y 900 nm de longitud de onda, en la llamada “ventana óptica o terapéutica”,6 y es poco absorbida por los componentes más significativos del tejido blando pudiendo así penetrarlos varios centímetros y explorarlos. La forma en que se produce esta propagación así como la cantidad de luz que es absorbida en cada longitud de onda depende de las características fisiológicas y metabólicas del tejido recorrido.7 El estudio de la propagación de luz de longitud de onda NIR en medios difusivos, es decir, con alto scattering, es de gran valor debido a sus potenciales aplicaciones en imágenes de tejidos tales como mamas, cerebro y articulaciones. Por lo tanto, resulta ser una importante herramienta para el diagnóstico de enfermedades así como, en algunos casos, para el tratamiento de las mismas. En este tipo de análisis, una de las técnicas que se utilizan es la tomografía óptica difusa (DOT) que permite la reconstrucción de la distribución espacial de la propiedades ópticas a partir de las medidas con NIR.2,4,8,9
Entre todas las aplicaciones de esta técnica en el campo de la biomedicina, una de las potenciales aplicaciones es el diagnóstico temprano y seguimiento de tratamiento del cáncer de mama. Esta enfermedad causa en Argentina 5400 muertes por año, y es el tipo de cáncer que mayor número de personas afecta.10
Si bien existen técnicas estandarizadas empleadas en el algoritmo bien establecido para el diagnóstico de cáncer de mama, presentan algunas desventajas: la mamografía tradicional de rayos x emplea radiación ionizante, limitando su uso en el tiempo, y si bien es una técnica de alta resolución, carece de especificidad. Por otro lado, la mamografía por ultrasonido, si bien no tiene ninguna contraindicación, es de muy baja resolución y técnico dependiente. Por este motivo, resulta de gran interés desarrollar una técnica diagnóstica complementaria, como lo es la mamografía óptica, de la que se puede obtener información metabólica, permitiendo caracterizar los tipos de lesiones presentes. Este aporte significativo en el diagnóstico, significaría la detección temprana, lo que aumenta la sobrevida de las mujeres con cáncer de mama, y permitiría contar además con una herramienta para evaluar la respuesta a un tratamiento de esta enfermedad.
Se presenta aquí un estudio de la plausibilidad del uso de la técnica de óptica difusa para sistemas difusivos que contienen regiones localizadas donde la aproximación difusiva no es válida, como lo son los quistes en mamas. La importancia de este de estudio no radica sólo en la detección de este tipo de alteraciones, sino además en el ajuste y calibración de esta técnica de diagnóstico por imágenes, pretendiendo incluir todos los posibles casos de estudio reales, donde pueden presentarse alteraciones de diferentes características ópticas; introduciéndose en un tipo de estudio poco tratado.11-13
A continuación se presenta un resumen de las consideraciones teóricas, la metodología experimental para la validación de los modelos teóricos empleados, los resultados y las principales conclusiones de este análisis.
II. CONSIDERACIONES TEÓRICAS
Existen diversas técnicas para estudiar medios turbios, y obtener información de importancia diagnósticas: pueden utilizarse diferentes tipos de fuentes y de detectores, acorde a la fuente elegida. Según los diferentes esquemas de medición, los sistemas pueden ser resueltos en el tiempo (TR), de onda continua (CW) o resueltos en el dominio de la frecuencias (FD). El primero requiere una fuente de onda continua y un detector tipo cámara CCD o medidor de potencia. Este esquema tiene la desventaja de que la información que brinda acerca de las propiedades ópticas es acoplada, es decir, no se pueden obtener los valores de los propiedades ópticas que describen al sistema, coeficiente de absorción y coeficiente de scattering reducido de manera independiente. Sin embargo, sirve para estudiar variaciones relativas de los mismos, determinando la presencia de inhomogeneidades localizadas en los medios de interés, además de ser una técnica de bajo costo y de muy rápida implementación, demostrando que puede ser utilizada para una primera localización.5
Para el segundo caso, se usan fuentes pulsadas, que constan de un pulso corto, típicamente del orden de los picosegundos. Al ingresar al medio, este pulso sufre un ensanchamiento, del que puede extraerse información acerca de las propiedades ópticas de dicho medio en estudio. A diferencia de las fuentes continuas, tienen la ventaja de que éstas pueden obtenerse de manera independiente a partir del perfil temporal de la distribuciones de tiempo de vuelo de los fotones difundidos.
En el caso trabajar en el dominio de frecuencias, el pulso de entrada se encuentra modulado en amplitud y fase, éstas sufren un cambio y un corrimiento respectivamente, luego de atravesar el medio turbio, de los que pueden establecerse de manera independiente los coeficientes de absorción y de scattering del medio de interés.
Pueden vincularse los tres esquemas, en particular es posible comparar el resultado de integrar sobre el tiempo los pulsos de TR con los valores de intensidad de CW.5,12,14 También, pueden conectarse los datos en el régimen resuelto en el tiempo con FD mediante una transformada de Fourier para pasar al dominio de frecuencias (es usual seleccionar amplitud y fase correspondientes a una frecuencia de 100 MHz8).
Modelos de transporte de la luz
Si bien la descripción más completa para la propagación de la luz en medios biológicos es a través de la ecuaciones de Maxwell, la complejidad matemática es muy elevada y, por lo tanto, su aplicabilidad práctica es limitada. El segundo enfoque posible es a través de la ecuación de transferencia radiativa (RTE), que es un acercamiento heurístico, que no comienza con la ecuaciones de onda, sino que se desprende de un análisis de balance energético. Si bien carece de la rigurosidad matemática de la teoría analítica, tiene la ventaja de ser aplicable a una gran variedad de problemas prácticos, dando resultados extremadamente precisos, convirtiéndose en una descripción muy aceptada y la más utilizada en problemas de transporte de luz en medios biológicos.15,16
En el caso más general dependiente del tiempo, la RTE puede escribirse como
donde es a intensidad, es la velocidad de la luz en el medio, es el término fuente, que es la potencia emitida en el tiempo por unidad de volumen, en la posición y de ángulo sólido a lo largo de la dirección , y es el elemento de ángulo sólido en la dirección y la función de scattering. Considerando un el balance de energía en un volumen , el primer término representa el cambio temporal total de energía, que se propaga a lo largo de en el elemento de volumen , en el ángulo sólido en . Es proporcional a la variación temporal total del número de fotones moviéndose en , y ; el segundo término representa el flujo de energía neto que se propaga a lo largo de a través del volumen en el ángulo en . El término es usualmente utilizado como ; el tercer término representa la fracción de energía propagándose en la dirección dentro de , y , debido a los fenómenos de scattering y absorción. Del segundo miembro de la ecuación para completar el balance de energía, el primer término representa la energía que viene de la dirección que dentro de , y es dispersada en la dirección y el último representa la energía generada a lo largo de en y por fuentes internas en .
Si bien, la RTE determina correctamente la propagación de la luz en los tejidos, teniendo en cuenta la absorción y los múltiples eventos de scattering debido a las heterogeneidades del medio en estudio, no posee una solución analítica para geometrías arbitrarias, y los métodos numéricos son muy costosos computacionalmente.
Para considerar soluciones, se emplean aproximaciones: el modelo de aproximación más usado es la aproximación difusiva (DA), Éste supone que la propagación de luz se vuelve isotrópica debido al alto scattering, cuyo valor es aproximadamente dos órdenes de magnitud mayor que la absorción. En el caso más general, se utilizan dos suposiciones simplificadoras. La primera es que la radiancia dentro del medio es casi isotrópica. La intensidad difusa se aproxima por los primeros dos términos de la expansión en serie en esféricos armónicos, es decir, los términos isotrópico y linealmente anisotrópico:
Esta expansión en armónicos esféricos truncada en el segundo término se llama aproximación .15 La Ec. (2) es una buena aproximación de la radiancia cuando los órdenes superiores de la expansión en esféricos armónicos son despreciables, lo que es usualmente cierto cuando el segundo término de la expansión es pequeño comparado con el primero, es decir, .
La segunda suposición que se hace es que la variación temporal del vector flujo en el tiempo es despreciable respecto al vector en sí mismo, es decir,
Dada esta limitación, la DA falla en describir correctamente la propagación de la luz en regiones donde el scattering es bajo o nulo. Ejemplos de regiones de bajo scattering en DOT son el fluido cerebro-espinal, fluido sinovial en articulaciones y quistes en mamas.
Sin embargo, es posible analizar estas situaciones utilizando la aproximación difusiva para modelar el transporte de la luz en todo el volumen del fantoma, pero compensando los errores de utilizar dicha aproximación en lugar de modelarlo con la RTE, particularmente para las regiones de bajo o nulo scattering.17 El Error de modelo corrige errores en la utilización de modelos físicos inadecuados, y puede emplearse para corregir errores de discretización de contornos físicos o de geometría.8,18,19 Este método ha demostrado ser útil para su aplicación en tomografía óptica.8,11,19
Para estudiar este método, se aplica la aproximación Bayesiana11 para compensar las discrepancias de la aproximación difusiva en regiones del volumen de estudio cercanas a fuentes o en casos de subdominios de bajo o nulo scattering. Se utiliza la RTE como el modelo exacto de transporte de luz, y las muestras de error son construidos en base a las predicciones de modelos difusivos (DA) y de transferencia radiativa. Puede hacerse una breve descripción como sigue: sean los parámetros ópticos e las medidas, vectores finitos.
En DOT, el modelo de ruido aditivo en el modelo de observación tiene la forma:
donde, son la medidas, los parámetros ópticos, es el modelo preciso que mapea la propiedades ópticas y el ruido. En general, el parámetro es aproximado a una base finita y la ecuación continua (4) es aproximada a una discreta de la forma:
donde es el parámetro discretizado, el parámetro de malla que controla el nivel de discretización, que se realiza para resolver el problema inverso. Así, si el modelo continuo puede discretizarse en un modelo de dimensión finita, el modelo discretizado es exacto dentro de la precisión de la medida, e , con pequeño, el modelo de observación será:
donde es el modelo reducido, en este caso la DA, es el modelo exacto de la RTE y el error de modelado. Bajo este esquema de aproximación de error, las reconstrucciones de baja dimensión DA son comparables a las de RTE. Por lo tanto es posible compensar el error de utilizar la DA en un esquema donde la RTE debe ser utilizada por las características ópticas de regiones presentes en el volumen de estudio.
III. METODOLOGÍA EXPERIMENTAL
Se estudiaron fantomas difusivos, fabricados ad hoc con características ópticas coincidentes con las del tejido mamamario,7 con inclusiones localizadas, de diferente niveles de turbidez. A continuación se describen las dos metodologías utilizadas, validas en fantomas con inclusiones difusivas,1,2,4,5 que fueron testeadas para este estudio con heterogeneidades localizadas traslúcidas.
Las propiedades ópticas del huésped y de las inclusiones fueron determinadas previamente, a través de un ajuste por cuadrados mínimos.20
Medidas topográficas
Este tipo de esquema se propone para el estudio en el plano 2d. En este trabajo se estudiaron dos técnicas: transmitancia de campo completo y transiluminancia.
Cada fantoma estudiado en el esquema de transmitancia es iluminado en una de las caras por una fuente láser y mediante la ayuda de una cámara CCD se toma la imagen de la otra cara. La lente de la cámara forma una imagen de campo completo en la cara de salida del fantoma sobre el arreglo de CCD. El láser es colocado en el eje óptico.
Como fuente de excitación se utiliza un láser de longitud de onda de 785 nm de Becker & Hickl GmbH, modelo BHLP-700, de potencia promedio 5 mW, que mediante una fibra óptica de 600 (m de diámetro y ilumina al medio de interés, de la cara opuesta a la cámara, sobre la linea del eje óptico. Como sistema de detección se emplea una cámara EMCCD (“Electron Multiplying Charge Coupled Device”), Andor Ixon Ultra 987, que toma imágenes con tiempos de adquisición de , tomada con una resolución de pixels y tamaño determinado por la geometría de estudio, es decir, la extensión del fantoma y del campo visual de la lente, formando una imagen de cm.
Las inhomogeneidades tienen una contribución muy sutil en la imagen, que no es posible de detectar en un simple análisis de la misma. La influencia de inclusiones de diversas características ópticas respecto del fondo o huésped puede hacerse evidente si se aplica un proceso de normalización adecuado, dividiendo a la imagen de un medio con inclusión por la correspondiente de un medio homogéneo (sin inclusiones) de las mismas características ópticas. Sin embargo, en el caso práctico real no es posible tomar la imagen del medio sin ellas, por lo que se toman múltiples exposiciones para obtener un fondo sintetizado del promedio de las mismas.1
Estas medidas son directamente contrastables con un estudio de transiluminancia. Para este se emplea un láser pulsado, en este caso el mismo láser de 785 nm, de Becker & Hickl GmbH, con duración de pulso de 120 ps, y colecta mediante una fibra óptica colineal a la incidencia de la otra cara, la intensidad temporal trasmitida, empleando como detector un fotomultiplicador, Hamamatsu PCM-100. Las medidas de intensidades se colectan en varias posiciones a lo largo de una línea o en una grilla en el plano , manteniendo colineales las fibras de fuente y de colección, obteniendo los valores de intensidades integradas punto a punto, resultado en los perfiles de transiluminancia en una y dos dimensiones (se emplea el mismo esquema experimental de la Fig. 1).
Medidas tomográficas
Para los experimentos tipo tomográficos, las fibras de colección y detección no se mueven colineales, sino que cada una barre un área: para la fuente se utilizan posiciones, distanciadas 10 mm, barriendo un área total de 900 mm2 de la cara entrante, y para el detector se emplean 81 posiciones, distanciadas 5 mm barriendo un área total de 1600 (40 x 40) mm2 de la cara saliente. Todas las posiciones se combinan, tomando un total de pares fuente-detector.
Los datos, que son adquiridos en el régimen resuelto en el tiempo, se transforman al dominio de las frecuencias mediante la transformada de Fourier, seleccionando la amplitud y fase correspondientes a MHz. Este valor de frecuencia es típicamente utilizado en este clase de procedimientos.8,17
La reconstrucción tomográfica de la que se obtiene una caracterización óptica de todo el volumen de estudio, aplicando la corrección del modelo para la recuperación de los parámetros ópticos. El proceso de reconstrucción tomográfica, que se encuentra desarrollado en Baez 2017,8 son reconstrucciones obtenidas a partir de la aplicación del Filtro de Kalman Extendido. Consiste en un problema de estimación de parámetro dada una medida compleja (se trabaja en el dominio de la frecuencias) , con el número de fuentes y el número de detectores.
IV. RESULTADOS
Se presentan algunos resultados representativos, de un fantoma con un huésped difusivo e inclusiones de dos tipo. En la Tabla 1, se indican las propiedades ópticas y posiciones de las inclusiones, de 1 cm de dimensión, están referidas al centro de de las imágenes y el respecto de la cara externa del fantoma. La inclusión 1 imita una lesión tipo tumor, de alto scattering del orden del huésped, y de mayor absorción, como sucede en el tejido mamario, mientras que la Inclusión 2 refiere a una lesión tipo quiste, de nulo scattering y absorción baja, como la del agua. Los experimentos realizados permiten dilucidar la validez de los esquemas experimentales empleados en los tipos de fantomas seleccionados: sistemas difusivos con heterogeneidades localizadas difusivas y no difusivas.
Medidas topográficas
En la Fig. 2 se observan dos imágenes: la obtenida para trasmitancia y la obtenida para transiluminancia de un barrido en 2d del mismo fantoma. Ambas imágenes, tienen el mismo tamaño.
De las figuras de contraste, puede asegurarse la presencia de inclusiones, debido a la diferencia relativa de absorción o scattering de las mismas respecto del huésped. Utilizando este análisis de normalización de imágenes no es posible determinar el valor de las propiedades ópticas de las inclusiones, pero sí es posible detectar variaciones de la atenuación efectiva del medio estudio, es decir, huésped e inhomogeneidades.
Si bien los valores de contraste son comparables, las posiciones de los valores extremos resultantes en los perfiles correspondientes pueden considerarse coincidentes, dentro de un margen, con las posiciones reales de las inclusiones. Los corrimientos de las manchas que representan las inclusiones en los respectivos fantomas heterogéneos, pueden atribuirse a que las inclusiones no se encuentran localizadas en el eje óptico del sistema fuente-cámara en todos los casos, como ya se indicara en los resultados de transmitancia de fantomas líquidos, debido a la distribución de trayectorias más probables de los fotones, corrimientos esquematizado en la Fig. 3.
Medidas tomográficas
En la Fig. 4 se presentan dos resultados de representativo de cortes en z del fantoma testigo, cuyos valores de propiedades ópticas indicadas en la escala de la figura, fuera reconstruidas con el modelo de DA, y con la corrección de modelo, EM.
De la reconstrucción basada en DA, el coeficiente de absorción, , muestra un valor aproximado del huésped de 0.004 mm-1, dicho valor es tres veces menor al valor determinado por el método de propiedades ópticas (Tabla 1), y sólo es posible distinguir algunas zonas de menor absorción localizadas, pero en las capas cercanas a las caras de entrada y salida, y cuyas posición no coincide con la posición nominal de ninguna de las inclusiones. Estas zonas pueden interpretarse como artefactos, que no pueden vincularse a ninguna inhomogeneidad. Por su parte, el coeficiente de scattering reducido muestra una sobreestimación de dicho coeficiente para todas las capas, y cuyo valor ronda en . En las capas centrales aparecen zonas de valores diferenciados en el centro, pero que no llegan a tener una forma distinguida (circular, vinculado a las inclusiones presentes).
Es evidente que en esta situación una reconstrucción basada en la DA requiere de algún tipo de consideración, es decir, debe emplearse la metodología de error de modelo. En la columna derecha, donde se considera esta reconstrucción, que presenta un mejoría, evidenciándose la presencia de ambas inclusiones.
Más específicamente, se obtiene del mapa de propiedades ópticas, se obtiene un huésped de mm-1, y dos zonas localizadas de con mayores valores y otra con menores valores, en coincidencia con la presencia de una inclusión tipo tumor y otra tipo quiste, respectivamente. La primera alcanza valores de hasta 0.005 mm-1, que es un poco mayor al coeficiente de absorción del agua; y la segunda valores de 0.022 mm-1, un poco por debajo del valor establecido mediante la técnica ajuste de DTOF por transmitancia resuelta en el tiempo, sin embargo, manteniendo la relación de veces entre huésped e inclusión, como se indica en la Tabla 1. Si bien logra reconstruir correctamente los valores, no localiza exactamente las inclusiones en la dirección , no obstante sí logra hacerlo correctamente en el plano , donde las distribuciones coinciden con las posiciones nominales de las heterogeneidades.
En el caso del coeficiente de scattering reducido recuperado, se aprecia una vez más que el valor del fondo, es decir, del huésped logra reconstruirse correctamente, presentando una fluctuación entre los valores 0.8 mm-1 y 1.2 mm-1. De particular interés, se observa una zona donde el es considerablemente menor, alcanzando valores de 0.4 mm-1 y hasta 0.3 mm-1. Esta región es donde se encuentra localizada la inclusión tipo quiste. Es interesante también destacar que el resultado de la reconstrucción de reducido refleja que la inclusión tipo tumor no se distingue, lo que es de esperar, ya que su valor es el mismo que el del huésped, como se indica en la Tabla 1.
Estos resultados muestran la aplicabilidad del error de modelo, demostrando que es posible utilizar una reconstrucción tomográfica basada en la aproximación difusiva en medios con heterogeneidades localizadas inmersas que no cumplen con las condiciones de dicha aproximación.
V. CONCLUSIONES
Ambos esquemas de medidas, permitieron validar la técnica de óptica difusa para sistemas difusivos, con diferente tipo de heterogeneidades, abarcando diferentes casos clínicos de lesiones mamarias recurrentes. Se llevó adelante la evaluación de experimentos de óptica biomédica, empleando técnicas de onda continua y resueltas en el tiempo. Las medidas experimentales se realizaron bajo dos esquemas: topográfico y tomográfico. El primero, brindó una localización en el plano de las heterogeneidades inmersas y un análisis cualitativo de las mismas, permitiendo probar arreglos experimentales ya establecidos. El segundo, proveyó una caracterización específica de los medios en estudio a través de una reconstrucción en tres dimensiones de las propiedades ópticas de un fantoma huésped difusivo con inclusiones de dos tipos: difusiva y no difusiva, lo que representa además una evaluación experimental del método de error de modelo.