Calibración y verificación de sistemas de imagen para guiado
Un sistema de guiado intervencionista no opera sobre imagen pura. Opera sobre correspondencias: entre lo que muestra la imagen y lo que ocurre en el espacio físico del procedimiento.
Esa correspondencia no está garantizada. Se establece, se verifica y se mantiene. La calibración no es un paso inicial de configuración que se hace una vez. Es una función continua de la que depende la exactitud del guiado.
1. El problema de la correspondencia geométrica
En diagnóstico, la imagen informa. En intervención, la imagen localiza. La diferencia operativa es que localizar exige una correspondencia cuantitativa entre la representación y el espacio físico.
Un sistema de imagen en intervención tiene que responder a preguntas del tipo: ¿dónde está exactamente este punto visible en imagen, en el espacio real del paciente? ¿A qué coordenadas del espacio quirúrgico corresponde este punto del detector?
Esas preguntas tienen respuesta solo si el sistema tiene un modelo geométrico válido y calibrado.
El problema surge porque los sistemas de imagen en intervención no tienen geometría fija:
- Un C-arm rota y su geometría cambia con cada posición del arco.
- Una cámara gamma portátil puede moverse libremente en el campo quirúrgico.
- Un sistema óptico acoplado a imagen de rayos X tiene una relación espacial relativa que debe ser caracterizada.
Sin un modelo geométrico calibrado que relacione la imagen con el espacio físico, la posición de un instrumento visible en pantalla es información cualitativa, no una localización cuantitativa.
2. Modelos geométricos de proyección
El modelo estándar para imagen proyectiva (rayos X, fluoroscopia, C-arm) es la proyección perspectiva central, conocido también como modelo pin-hole o modelo de cámara central.
Este modelo describe la relación entre un punto del espacio 3D y su proyección en el plano del detector mediante un conjunto de parámetros:
Parámetros intrínsecos: describen la geometría interna del sistema de adquisición.
- Distancia focal efectiva (en el modelo de rayos X: distancia foco-detector).
- Posición del punto principal (proyección del foco sobre el detector).
- Factor de escala entre coordenadas del detector (píxeles) y unidades físicas.
- Parámetros de distorsión (del detector, de la óptica en cámaras convencionales).
Parámetros extrínsecos: describen la posición y orientación del sistema respecto a un sistema de referencia externo.
- Matriz de rotación (3×3) que describe la orientación.
- Vector de traslación que describe la posición del foco en el sistema de referencia.
La calibración geométrica consiste en determinar estos parámetros con suficiente precisión para la aplicación concreta.
Ecuación de cámara y magnificación geométrica
La forma compacta de escribir el problema es la ecuación de cámara:
x ~ K [R | t] X
donde X es un punto del espacio 3D en coordenadas homogéneas, x es su proyección en el detector, K contiene los parámetros intrínsecos, y [R | t] describe la posición y orientación del sistema respecto al marco de referencia elegido.
En imagen óptica esta ecuación se interpreta como la proyección desde un centro de cámara sobre un sensor. En rayos X, la analogía física es directa: el foco actúa como centro de proyección y el detector como plano imagen. La diferencia es que aquí la magnificación tiene una lectura geométrica inmediata en el espacio clínico.
Si p_d es el tamaño de píxel en el detector, SID la distancia fuente-detector y SOD la distancia fuente-objeto, la magnificación geométrica es:
M = SID / SOD
y el tamaño de píxel efectivo en el plano del objeto es aproximadamente:
p_obj = p_d / M
Esto importa al informar resolución o error de localización. En un detector con p_d = 0,194 mm, un SID = 1200 mm y un objeto situado cerca del isocentro con SOD = 750 mm, la magnificación es M = 1,6 y el tamaño de píxel efectivo en ese plano es de unos 0,12 mm. Es decir: el tamaño de píxel expresado en el plano del isocentro no coincide con el tamaño físico del píxel en el detector.
El isocentro y su estabilidad
En sistemas de C-arm, el isocentro es el punto del espacio al que convergen todas las proyecciones durante una rotación. En un C-arm isocéntrico ideal, ese punto permanece fijo independientemente del ángulo del arco, lo que permite combinar proyecciones desde distintos ángulos para reconstruir volúmenes (tomografía cone-beam con C-arm).
En la práctica, la mecánica real introduce desviaciones respecto a la geometría nominal. El isocentro real puede no coincidir con el nominal, y puede desplazarse a lo largo de la rotación. Caracterizar esas desviaciones — la variación del isocentro en función del ángulo — forma parte de la calibración geométrica del sistema.
Para aplicaciones de reconstrucción tomográfica con C-arm, la exactitud de la geometría tiene un impacto directo en la calidad de la reconstrucción: errores geométricos no compensados generan artefactos y pérdida de resolución. En guiado 2D, el efecto es distinto pero igualmente crítico: la correspondencia entre imagen y espacio real depende de que los parámetros extrínsecos del sistema en cada posición del arco sean conocidos.
3. Fantomas de calibración
La calibración geométrica requiere objetos de referencia con posiciones conocidas en el espacio y detectables en imagen. Esos objetos se denominan fantomas de calibración.
Tipos de fantomas
Fantomas geométricos: contienen marcadores de posición conocida (esferas metálicas, hilos de tungsteno, rejillas) que son visibles en imagen de rayos X. La relación entre las posiciones de los marcadores en el espacio y sus proyecciones en el detector permite estimar los parámetros del modelo geométrico.
El diseño del fantoma condiciona la calidad de la calibración: la distribución y número de marcadores, su tamaño relativo al campo de visión, y su disposición tridimensional determinan la condicionabilidad del sistema de ecuaciones que se resuelve para obtener los parámetros.
Fantomas de tejido equivalente: no sirven para calibración geométrica directamente, pero permiten evaluar el rendimiento del sistema en condiciones próximas a las clínicas. Útiles para verificación dosimétrica y para evaluar contraste y resolución en condiciones de atenuación realista.
Fantomas multimodales: permiten verificar la correspondencia entre dos modalidades de imagen (por ejemplo, TAC y fluoroscopia, o imagen óptica e imagen de rayos X). Deben ser visibles en ambos sistemas con suficiente contraste y con marcadores localizables en cada uno.
Resolución y trazabilidad
El error de localización de los marcadores en imagen limita la exactitud de la calibración. Un marcador de 1 mm de diámetro localizado con una incertidumbre de 0,5 mm en el detector introduce un error de calibración que se propaga al modelo geométrico.
La trazabilidad metrológica de la posición de los marcadores en el fantoma (medida con incertidumbre conocida, referenciada a patrones físicos) es la base sobre la que se sostiene el proceso de calibración. Sin trazabilidad, la calibración es relativa pero no absoluta.
4. Calibración conjunta de sistemas RX y ópticos
Una línea de trabajo especialmente relevante en guiado intervencionista es la combinación de imagen de rayos X con imagen óptica visible (cámaras RGB o RGB-D).
La motivación es directa: los rayos X proporcionan visibilidad de estructuras internas, pero no aportan información superficial o de la escena quirúrgica. Una cámara óptica captura la geometría superficial y la posición de instrumentos, pero no ve a través del tejido. Combinar ambos canales amplía la información disponible y permite relacionar el espacio visual del cirujano con la imagen proyectiva interna.
La calibración conjunta de un sistema de rayos X y un sistema óptico requiere determinar la transformación rígida entre ambos sistemas de coordenadas: la rotación y traslación que relaciona las coordenadas del espacio óptico con las coordenadas del espacio físico del sistema de rayos X.
Este problema es análogo al de calibración extrínsecos entre dos cámaras ópticas (calibración estéreo), pero presenta complicaciones adicionales:
- Los dos sistemas operan en longitudes de onda muy distintas, por lo que los mismos objetos físicos pueden no ser visibles en ambas modalidades o no tener el mismo aspecto geométrico.
- La resolución espacial y las distorsiones son diferentes en cada sistema.
- El detector de rayos X puede no estar fijo en posición conocida respecto a la cámara óptica (como ocurre en un C-arm móvil).
La solución habitual es diseñar fantomas multimodales con marcadores específicamente diseñados para ser detectables en ambas modalidades, y resolver conjuntamente los parámetros geométricos de ambos sistemas usando las posiciones correspondientes de esos marcadores.
Restricción epipolar y triangulación
Cuando dos sistemas proyectivos observan el mismo punto 3D, las correspondencias entre ambas imágenes no son arbitrarias. Quedan restringidas por la geometría epipolar: un punto identificado en una imagen debe proyectarse en la otra sobre su recta epipolar, no en cualquier posición del plano.
En formulación matricial, si x y x' son dos observaciones correspondientes del mismo punto, la restricción ideal se expresa como:
x'^T F x = 0
donde F es la matriz fundamental. Si además los parámetros intrínsecos de ambos sistemas son conocidos, puede escribirse la relación equivalente en coordenadas normalizadas mediante la matriz esencial E.
La utilidad práctica es doble. Primero, reduce el problema de búsqueda de correspondencias: no hay que buscar en toda la imagen, sino sobre una recta. Segundo, permite triangular la posición 3D del punto como intersección, en el sentido geométrico o de mínimos cuadrados, de dos rayos de proyección calibrados. En sistemas RX-ópticos y en navegación radioscópica, esta restricción es una pieza central para pasar de dos imágenes 2D a una localización espacial cuantitativa.
Reconstrucción 3D a partir de proyecciones RX con referencia óptica
Una aplicación directa de la calibración conjunta es la posibilidad de recuperar información tridimensional a partir de proyecciones de rayos X usando la geometría óptica como referencia.
Dos proyecciones desde ángulos conocidos permiten, por triangulación, recuperar la posición 3D de un punto visible en ambas. Si la geometría del sistema (parámetros de cada proyección) está calibrada con suficiente exactitud, la posición 3D resultante puede relacionarse con el espacio óptico y, por tanto, con la escena quirúrgica.
Este principio está en la base de los sistemas de navegación radioscópica: usar múltiples proyecciones de rayos X con geometría calibrada para localizar instrumentos en 3D sin necesidad de TAC intraoperatorio o tomografía cone-beam.
5. Calibración, registro y localización: tres problemas distintos
En guiado intervencionista conviene separar tres problemas que a menudo aparecen encadenados pero no son equivalentes.
Calibración: estimar los parámetros del sistema de imagen o de la relación fija entre dos subsistemas. Responde a preguntas como: ¿cuál es la geometría real de este C-arm? ¿Cuál es la transformación rígida entre esta cámara RGB y este sistema de rayos X?
Registro: alinear imágenes, volúmenes o modalidades distintas en un mismo marco espacial. Responde a preguntas como: ¿cómo hago coincidir un TAC preoperatorio con una resonancia magnética? ¿Cómo relaciono una imagen gamma con una imagen anatómica obtenida por otra modalidad?
Localización: situar al paciente, al instrumento o al dispositivo en el espacio físico del procedimiento. Responde a preguntas como: ¿dónde está ahora mismo la aguja de biopsia? ¿Dónde está el extremo de un instrumento robotizado? ¿Cómo queda referido un aplicador de radioterapia intraoperatoria al campo real de intervención?
No son sinónimos. La calibración caracteriza el sistema. El registro relaciona representaciones. La localización conecta el procedimiento real con el sistema de coordenadas operativo.
Esta distinción aparece también en otros contextos del bloque, por ejemplo en radioterapia, donde el registro alinea datasets y la localización coloca al paciente real dentro del sistema geométrico del equipo. Véase Radioterapia.
Ejemplos típicos:
- TAC-RM: problema de registro entre datasets.
- Fluoroscopia + cámara RGB/RGB-D: problema de calibración conjunta entre sistemas, con posible registro posterior a otras imágenes.
- Aguja visible en fluoroscopia durante una biopsia: problema de localización en tiempo de intervención.
- Sistema operatorio o plataforma tipo da Vinci: problema de localización y seguimiento de instrumentos, normalmente apoyado sobre calibraciones previas y, en algunos flujos, sobre registro con imagen preoperatoria.
6. Registro y referenciación espacial de cámaras gamma intervencionistas
Las cámaras gamma portátiles e intervencionistas presentan un problema geométrico y operativo distinto al de los sistemas de rayos X.
Un sistema de rayos X forma imagen por proyección: cada punto del detector corresponde a una trayectoria espacial definida por la geometría de proyección. La calibración establece esa relación.
Una cámara gamma en configuración intervencionista no se organiza alrededor de un isocentro como un C-arm. Su problema operativo suele ser otro: cómo referir la información funcional al paciente, al campo quirúrgico o a otra modalidad anatómica. Eso puede resolverse mediante registro imagen-imagen, registro imagen-espacio o tracking externo del detector.
Además, no todas las cámaras gamma forman imagen del mismo modo:
- En colimación paralela, cada punto del detector corresponde a una región del espacio en la dirección de los canales del colimador. La relación es aproximadamente directa, pero la eficiencia y la resolución varían con la distancia.
- En colimación estenopeica (pinhole), la geometría es más próxima a una proyección perspectiva, pero con campo de visión muy limitado y fuerte dependencia de la distancia. Es un caso proyectivo particular, no el equivalente operativo de un C-arm isocéntrico.
- En cámaras Compton, no hay colimador físico: cada evento detectado define un cono de probabilidad en el espacio, y la imagen se reconstruye por acumulación estadística de muchos eventos.
Por eso, en gamma intervencionista la pregunta principal no suele ser “¿cuál es el isocentro del sistema?”, sino más bien: ¿cómo conecto la señal funcional detectada con la anatomía del paciente o con el espacio físico del procedimiento?
Para responderla hacen falta una o varias de estas operaciones:
- referenciación espacial del detector: conocer la posición y orientación de la cámara gamma en el campo quirúrgico mediante navegación óptica, marcadores electromagnéticos o rastreo mecánico;
- registro con otra modalidad: alinear la imagen gamma con una imagen anatómica previa o intraoperatoria, como TAC, resonancia magnética o fluoroscopia;
- localización clínica: usar esa cadena de correspondencias para decidir dónde está una fuente, una lesión marcada o un blanco funcional durante el procedimiento.
La relación fija entre la cámara gamma y su sistema de rastreo sí puede requerir una calibración extrínseca. Pero la tarea clínica relevante suele ser de registro y referenciación espacial, no de geometría isocéntrica.
En una gammacámara colimada portátil, el QA relevante no se limita a la geometría: incluye también estabilidad de sensibilidad, resolución espacial dependiente de la distancia, reproducibilidad del registro con el sistema anatómico y exactitud de la referenciación en el espacio de trabajo.
En una cámara Compton, además, hay una capa adicional de reconstrucción:
- Calibración energética de cada detector individual.
- Caracterización de la resolución espacial y energética de cada canal.
- Verificación de la estabilidad de la respuesta con temperatura y con tiempo de uso.
- Validación de la reconstrucción estadística y de la referenciación espacial mediante fuentes de posición conocida.
7. Dispositivos de uso dual: medida espacial e imagen de radiación
Una clase de dispositivos de especial interés en investigación y desarrollo es la que combina capacidad de imagen de radiación con medida espacial en un único sistema.
La motivación es que en intervención guiada, la imagen funcional (gamma, fluorescencia) y la imagen estructural (RX, óptica) aportan información complementaria. Un dispositivo que puede adquirir ambas — o que puede relacionar ambas mediante calibración — reduce la complejidad del sistema global y el número de transformaciones geométricas que hay que encadenar.
Ejemplos de esta filosofía de diseño:
- Sistemas que combinan una sonda gamma con un puntero óptico localizable por cámara externa: la señal gamma se atribuye a la posición 3D del puntero en el espacio quirúrgico.
- Sistemas que combinan una cámara gamma compacta con marcadores ópticos integrados: la cámara óptica externa localiza la gamma-cámara y permite referir su imagen al espacio quirúrgico.
- Sistemas de rayos X con cámara RGB-D integrada: la profundidad del objeto visible en óptica se usa para refinar la reconstrucción 3D de la proyección de rayos X.
En todos estos casos aparecen, encadenados, los tres niveles anteriores:
- calibración entre sensores o subsistemas;
- registro entre imágenes o espacios de referencia;
- localización del instrumento o del blanco durante el procedimiento.
La verificación de la exactitud del sistema integrado — frente a la exactitud de cada subsistema por separado — es un requisito de validación que va más allá de la caracterización individual de los componentes.
8. Controles numéricos y verificación periódica
Ni la calibración, ni el registro, ni la localización operativa son estados permanentes. Toda la cadena puede degradarse por:
- Cambios mecánicos en el sistema (colimador, detector, brazo portante).
- Temperatura (afecta a detectores, electrónica de adquisición, geometría mecánica).
- Uso acumulado (desgaste de colimadores, degradación de detectores).
- Accidentes mecánicos (golpes, transporte).
- Actualizaciones de software de reconstrucción o adquisición.
- Cambios en el algoritmo o en el flujo de registro entre modalidades.
- Desplazamientos del paciente, del instrumental o del sistema de rastreo durante el procedimiento.
El protocolo de controles numéricos periódicos tiene que verificar que cada capa sigue dentro de los límites de exactitud especificados para su uso clínico.
Indicadores de control
Los controles periódicos deben incluir al menos:
Error de calibración geométrica: discrepancia entre la proyección esperada por el modelo y la posición realmente observada de marcadores de referencia. Debe medirse con fantomas adecuados y en posiciones representativas del campo de trabajo.
Error de registro: discrepancia residual al alinear dos imágenes, volúmenes o modalidades. Debe evaluarse sobre estructuras o marcadores conocidos, y no solo sobre el valor que optimiza internamente el algoritmo.
Error de localización: diferencia entre la posición medida por el sistema y la posición real conocida de un instrumento, marcador o blanco durante una tarea representativa. Expresado en mm, es el indicador más directamente conectado con el uso clínico.
Reproducibilidad: variación del error de localización entre sesiones o entre operadores, bajo condiciones equivalentes. Expresa la precisión del sistema en condiciones reales de uso.
Resolución espacial: capacidad del sistema para distinguir estructuras próximas. Para sistemas de imagen proyectiva, se mide con fantomas de barra o de cuña. Para sistemas de reconstrucción, se mide mediante la función de transferencia de modulación (MTF) o el perfil de respuesta de un punto (PSF).
En sistemas proyectivos conviene indicar siempre en qué plano se expresa esa resolución. El valor en el detector puede darse en mm/píxel o en pares de líneas por milímetro del propio detector, pero el valor relevante para guiado suele ser el del plano de interés clínico, a menudo próximo al isocentro. Por magnificación geométrica, una misma matriz de detector puede corresponder a un tamaño de píxel efectivo menor en el isocentro que en el plano detector, y esa conversión debe explicitarse.
Como orden de magnitud, un panel plano con píxeles de 0,15-0,20 mm puede corresponder a tamaños efectivos del orden de 0,10-0,13 mm cerca del isocentro para geometrías de trabajo habituales. Eso no significa que la exactitud global del guiado sea de 0,10 mm: la exactitud final incorpora además localización subpíxel, distorsión residual, incertidumbre del fantoma, estabilidad mecánica, registro entre modalidades y errores del algoritmo de reconstrucción o navegación.
Estabilidad de la geometría bajo carga y posición: verificar que la relación geométrica entre componentes del sistema no varía de forma significativa con la posición del arco o del portante durante el uso.
Robustez del registro y del tracking: verificar que la alineación entre modalidades y la referenciación espacial externa se mantienen dentro de tolerancia cuando cambian operador, iluminación, posición del detector o condiciones del entorno.
Frecuencia de controles
La frecuencia depende de:
- El tipo de procedimiento (guiado de alta precisión exige controles más frecuentes).
- La intensidad de uso del sistema.
- Las especificaciones del fabricante.
- Los resultados históricos de los propios controles.
Los controles no son mera documentación. Su función es detectar desviaciones antes de que afecten al procedimiento. Un sistema que supera el límite de exactitud especificado debe ser retirado de uso clínico hasta que se identifique la causa y se corrija.
9. Evaluación de sistemas de guiado: exactitud, precisión y validación
La evaluación de un sistema de guiado intervencionista se realiza en capas progresivas:
Evaluación en banco
Verificación del modelo geométrico, de la calibración, del registro y de la exactitud de localización en condiciones controladas de laboratorio, usando fantomas de posición conocida y bajo condiciones estables de temperatura e iluminación.
Esta evaluación establece el límite inferior del error del sistema: el mejor rendimiento alcanzable en condiciones ideales.
Evaluación en condiciones de uso realistas
Verificación del rendimiento en condiciones próximas al uso clínico: temperatura de quirófano, iluminación variable, presencia de otros equipos, movimiento del operador, variación de posición del sistema y cambios en el flujo real de registro o rastreo.
La degradación del rendimiento respecto a la evaluación en banco cuantifica el impacto de las condiciones operativas reales.
Evaluación en fantoma clínico
Verificación del rendimiento con fantomas de tejido equivalente que simulan la tarea clínica concreta: biopsia de nódulo, localización de ganglio, guiado de catéter. La evaluación mide el error de ejecución de la tarea completa, no solo el error geométrico interno del sistema.
Validación clínica
Demostración de mejora en indicadores clínicamente relevantes (exactitud de localización en pacientes reales, tasa de complicaciones, tiempo de procedimiento, número de exposiciones adicionales necesarias) en comparación con el estándar de referencia.
La validación clínica solo tiene sentido cuando las capas anteriores han establecido que el sistema funciona como se espera técnicamente. Invertir el orden — validar clínicamente antes de verificar técnicamente — mezcla las fuentes de error y dificulta la interpretación de los resultados.
10. Resumen
La calibración de un sistema de imagen para guiado intervencionista es una parte del problema general de correspondencia entre la representación en imagen y el espacio físico del procedimiento.
En la práctica conviene separar tres operaciones:
- Calibración: estimar los parámetros del sistema o de la relación fija entre subsistemas.
- Registro: alinear imágenes, modalidades o espacios de representación cuando hace falta combinar información.
- Localización: situar al paciente, el instrumento o el dispositivo en el espacio real donde se actúa.
En sistemas proyectivos como el C-arm, la calibración geométrica y la estabilidad del isocentro son centrales. En sistemas gamma portátiles o multimodales, el problema suele desplazarse hacia el registro y la referenciación espacial. En todos los casos, la cadena clínica completa exige no solo calibrar, sino también registrar cuando corresponda y localizar con exactitud suficiente para actuar.
Para ampliar
Una parte de la bibliografía más directamente relacionada con este tema procede de trabajos sobre calibración RX-óptica, reconstrucción geométrica y medida 3D desarrollados en este mismo entorno de investigación.
- Richard Hartley y Andrew Zisserman - Multiple View Geometry in Computer Vision, 2nd ed.
- Zhengyou Zhang - A Flexible New Technique for Camera Calibration
- Evaluation of modern camera calibration techniques for conventional diagnostic X-ray imaging settings
- Geometrical Calibration of X-Ray Imaging With RGB Cameras for 3D Reconstruction
- Joint calibration of RGB and X-ray cameras
- 3D measurements from X-ray images and dense surface mappings
- 3D measurements in conventional X-ray imaging with RGB-D sensors
- A geometric calibration method for cone beam CT systems
- Geometric calibration of a mobile C-arm for intraoperative cone-beam CT
- Self-Calibration of Cone-Beam CT Geometry Using 3D-2D Image Registration: Development and Application to Task-Based Imaging with a Robotic C-Arm
- A versatile calibration procedure for portable coded aperture gamma cameras and RGB-D sensors
- MedPhysics.i-do.science - Publications
