Biofisica Imagenes
Biofisica Imagenes
Biofisica Imagenes
Biofísica
3.a edición
Biofísica
3.a edición
André Aurengo
Catedrático de Biofísica
Facultad Pierre et Marie Curie, París
Thierry Petitclerc
Catedrático de Biofísica
Facultad Pierre et Marie Curie, París
Equipo de traducción:
Prof. Carmela Calés
Profesora titular, Dpto. de Bioquímica y Biofísica
Facultad de Medicina, Universidad Autónoma de Madrid
Prof. Emilio Marco
Profesor titular, Dpto. de Fisiología
Facultad de Medicina, Universidad Autónoma de Madrid
Prof. Roberto Marco
Catedrático, Departamento de Bioquímica y Biofísica
Facultad de Medicina, Universidad Autónoma de Madrid
BIOFÍSICA
No está permitida la reproducción total o parcial de este libro, su tratamiento informático, ni la transmisión
de cualquier otra forma o por cualquier otro medio electrónico, mecánico, por fotocopia, por registro u otros
métodos, sin el permiso previo y por escrito de los titulares del Copyright.
Edificio Valrealty
c/Basauri, 17, 1.a planta
28023 Aravaca (Madrid)
ISBN: 978-84-486-0855-2
MEDIO INTERNO
Capítulo 1. Compartimientos líquidos del organismo ........................................................................................................... 3
Líquidos y soluciones ................................................................................................................................................................. 3
El agua y los solutos en el organismo ........................................................................................................................................ 9
Medida del contenido en agua y de los solutos ........................................................................................................................ 12
Preguntas de opción múltiple (POM) ........................................................................................................................................ 14
Ejercicios ...................................................................................................................................................................................... 15
Capítulo 5. Difusión y convección simultáneas del solvente a través de una membrana ............................................... 65
Presión osmótica ......................................................................................................................................................................... 65
Ultrafiltración .............................................................................................................................................................................. 73
Preguntas de opción múltiple (POM) ........................................................................................................................................ 79
Ejercicios ...................................................................................................................................................................................... 80
Capítulo 6. Difusión y migración eléctrica simultánea de los iones a través de una membrana ................................... 83
Generalidades .............................................................................................................................................................................. 84
Efecto Donnan ............................................................................................................................................................................ 85
Potencial de electrodo ................................................................................................................................................................ 90
Electrodos selectivos ................................................................................................................................................................... 94
VI CONTENIDO
ELECTROFISIOLOGÍA
Capítulo 7. Electrofisiología celular ......................................................................................................................................... 101
Potencial de reposo celular ........................................................................................................................................................ 101
Potencial de acción ..................................................................................................................................................................... 106
Preguntas de opción múltiple (POM) ........................................................................................................................................ 121
Ejercicios ...................................................................................................................................................................................... 121
BIOFÍSICA SENSORIAL
Capítulo 10. Biofísica de las funciones sensoriales ............................................................................................................... 151
Las diversas funciones sensoriales ............................................................................................................................................ 151
Cadena de transmisión sensorial ............................................................................................................................................... 151
Características de los receptores biológicos ............................................................................................................................. 152
Características del mensaje sensorial ....................................................................................................................................... 153
Características de la biofísica sensorial ..................................................................................................................................... 154
RADIACIONES
Capítulo 13. Radiaciones electromagnéticas ......................................................................................................................... 205
Ondas electromagnéticas ........................................................................................................................................................... 205
Espectro de la radiación electromagnética ............................................................................................................................... 206
Fotones ......................................................................................................................................................................................... 207
Clasificación de las radiaciones electromagnéticas ................................................................................................................. 208
Dualidad onda-corpúsculo ........................................................................................................................................................ 209
Ejercicios ...................................................................................................................................................................................... 210
Capítulo 15. Interacciones entre las radiaciones ionizantes y la materia .......................................................................... 219
Interacciones de las partículas cargadas con la materia ......................................................................................................... 219
Interacciones de los neutrones con la materia ......................................................................................................................... 222
Interacciones de los fotones con la materia ............................................................................................................................. 222
Ejercicios ...................................................................................................................................................................................... 226
Capítulo 19. Higiene y protección en el empleo de las radiaciones ionizantes ................................................................. 263
Irradiación del público (personas no expuestas por su profesión) ........................................................................................ 263
Irradiación profesional ............................................................................................................................................................... 265
Principios de protección radiológica ........................................................................................................................................ 265
Protección radiológica de los trabajadores ............................................................................................................................... 267
Protección radiológica del público ............................................................................................................................................ 267
Protección radiológica y medicina ............................................................................................................................................ 268
Conducta a seguir en el caso de una contaminación accidental ............................................................................................ 269
Ejercicios ...................................................................................................................................................................................... 269
IMÁGENES
Capítulo 20. Imágenes analógicas y digitales ......................................................................................................................... 273
Características de las imágenes analógicas .............................................................................................................................. 273
Características de las imágenes digitales .................................................................................................................................. 277
Obtención de imágenes digitales ............................................................................................................................................... 278
Visualización de una imagen digital .......................................................................................................................................... 279
Ventajas de la imagen digital ...................................................................................................................................................... 281
Inconvenientes de la imagen digital .......................................................................................................................................... 281
Ejercicios ...................................................................................................................................................................................... 282
Anexo 10. Teorema del perfil central de la transformada de Fourier ...................................................................................... 396
Anexo 11. Splines ......................................................................................................................................................................... 397
Anexo 12. Identificación de los parámetros de un modelo ..................................................................................................... 397
Ésta es la tercera vez que escribo un prólogo para un manual de biofísica. La primera fue en 1977. Lo firmaba junto a Jean Pe-
rrin, que nos ha dejado recientemente. Ambos éramos autores del tomo II de los Elementos de biofísica. Estos, en dos grandes
tomos rojos, trataban de dar una imagen casi exhaustiva de la biofísica. La segunda vez, en 1982, la obra constaba de un único
volumen, esta vez, blanco. Se titulaba Biofísica, sin la connotación restrictiva del término «elementos». En esta obra, yo era a la
vez autor y editor, en el sentido americano del término. Éramos seis autores, todos pertenecientes al departamento de biofísica
y biomatemática de la Facultad Pitié-Salpêtrière.
Hoy, la obra que presento es esencialmente la de dos de aquellos seis, y que son, para mí, dos alumnos muy queridos, si
me puedo atrever a usar este término, que implica ¡que les he enseñado realmente algo!
André Aurengo y Thierry Petitclerc me han pedido unir mi nombre a los suyos y aparecer de nuevo como coautor. La razón
que invocan, y de la que les cedo la responsabilidad, es que varios capítulos están directamente inspirados en los que había
escrito yo en 1982: cuestión de filiación, dicen. ¿Ha sido debilidad, vanidad o incluso deseo pueril de supervivencia lo que me
ha empujado a aceptar?
«Morir, de acuerdo, hay que hacerlo. Pero terminar de vivir, ¡no!», dijo, en alguna parte, Claude Roy. Mucho menos numerosos
que antaño son los que esperan, en la actualidad, la eternidad reservada a los bienaventurados. Todavía más escasos son los
que quieren sobrevivir en los «mañana que cantan», que a su vez habrán contribuido a iniciar. En el momento actual, parece
que la única esperanza legítima es dejar tras de sí una obra, un patrimonio, un nombre: algo parecido a una supervivencia
simbólica... y provisional. Y es desde esta perspectiva más modesta de una obra por hacer desde la que acepto hacer aparecer
mi nombre al comienzo de este libro.
¿Obra finalizada? No, más bien obra siempre por continuar. ¿Obra de uno solo? Seguramente no, más bien obra colectiva,
de un grupo que transmite y desarrolla una herencia. Escribiendo el tomo I de los Elementos, tomaba el relevo de A. Strohl,
A. Dognon y L. Gougerot. Hoy mismo, es posible que se encuentren, alrededor de A. Aurengo y T. Petitclerc, algunos jóvenes
asistentes, alentados por ellos, preparándose para tomar a su vez el testigo. De esta manera, la obra persistirá en la fidelidad y
la innovación, mientras los más ancianos se pierden en la «calle de las bodegas oscuras».
Con respecto a la obra de 1982, que a su vez era de extensión reducida en comparación con los Elementos de los años
setenta, el espectro de los dominios biofísicos que se tratan en este libro se ha reducido aún más. La biofísica de la circulación
y la de la respiración han desaparecido. Lo propio ha ocurrido con los fenómenos de superficies, de rayos ultravioletas, de
electrofisiología cerebral. En contraposición, el estudio del medio interno, la biofísica de las radiaciones y, en especial, de la
imagen, han tomado un espacio considerable y constituyen la esencia de esta obra.
¿Por qué estas variaciones? En primer lugar, porque André Aurengo y Thierry Petitclerc han querido llegar hoy a lo que
será esencial para la medicina del mañana. Restringiendo el campo cubierto, lo han hecho más profundo, exigiéndose a
ellos mismos y a sus lectores una precisión de conceptos sin contemplaciones y una articulación rigurosa del razonamiento.
Los amantes de la «canción gris en la que lo impreciso se une a lo indeciso» tienen el riesgo de penar. El trabajo que se nos
ofrece es espléndido, por su pureza y su exigencia. No es posible hacer «zapping», actividad, por otra parte, familiar a nuestros
estudiantes. Por ello, este libro es un gran educador del espíritu.
Pero hay tal vez otra razón, a saber, el estatus epistemológico de esta disciplina que es la biofísica. No tiene la misma
coherencia interna que la bioquímica o la inmunología, que representan campos del saber bien definidos. La biofísica no
tiene realmente un dominio propio. Entonces, ¿podría tratarse únicamente de una sub-sección del Consejo Nacional de las
Universidades, cuyo único fin fuera la supervivencia? No. Lo que hace la unidad de la biofísica, y explica al tiempo su diversidad,
es la luz, la visión del mundo que hace llegar a los diversos dominios de las ciencias de la vida. La física está considerada,
justamente, como el patrón de referencia de las ciencias experimentales, por las posibilidades inigualables de formalización
matemática, por la existencia de robustas teorías, y por las capacidades extraordinarias de precisión como herramienta teórica.
Cada vez que un sector de las ciencias de la vida puede acceder a una potente modelización que sobrepasa el estadio de la
modelización puramente lógica, la intervención de la (bio)física aporta a dichas ciencias una riqueza y un rigor excepcionales
para alcanzar el formalismo matemático. Por ejemplo, en electrofisiología, la electrocardiografía está en el dominio de la
biofísica, mientras que la electroencefalografía no lo está en realidad, pues apenas sobrepasa el estadio descriptivo. Hemos
deducido, a través de esta constatación, la elección que se imponía para esta obra.
XII PRÓLOGO A LA SEGUNDA EDICIÓN
Así, la biofísica aparece –con la teoría matemática de las probabilidades− como una de las canteras de la medicina
cuantitativa, y constituye para el facultativo un progreso considerable en la dirección de la «medicina basada en la evidencia».
La capacidad de justificar toda decisión sobre argumentos sólidos será, en adelante, el objetivo de toda educación médica de
calidad.
Puede ser oportuno precisar a los futuros facultativos que si bien la medicina cuantitativa es, desde este momento,
necesaria, no es, en manera alguna, suficiente para una buena práctica médica. Tratar a un cuerpo enfermo y, con más razón, a
un órgano enfermo, no debe hacer olvidar a la persona dueña de ese cuerpo, o mejor dicho, que es ese cuerpo. Uno no puede,
por tanto, darse por satisfecho con la mirada fría que lanza el ingeniero sobre el motor que está construyendo o reparando. Este
ingeniero no habla con su motor. O también, como recordaba R. Debré en su célebre discurso de 1973 (con el que ya apelaba a
una nueva revolución de la formación médica), la palabra y la presencia son tan importantes como el fármaco. Toda medicina
de calidad es, antes que nada, una medicina de la persona para la persona. Por supuesto, la segunda debe a la primera su
competencia técnica y científica. Es a esta parte esencial, pero no única, de la actividad médica, a la que contribuye este libro
debido, en lo esencial, a André Aurengo y Thierry Petitclerc. Son, el uno y el otro, ingenieros de formación, egresados de la más
ilustre de nuestras escuelas. Pero han sabido resistir a la tentación mortífera de la mirada helada, y tanto el uno como el otro
son también, y ante todo, médicos en el sentido más fuerte y más noble del término, pues no olvidan jamás que tratan a seres
humanos, sus semejantes, sus hermanos.
François GRÉMY
Prólogo a la primera edición
No le es fácil a uno renegar de sus orígenes… A pesar de haber abandonado desde hace ya varios años la enseñanza de la
biofísica, para dedicarme a la de las ciencias de la información y a la de la epidemiología, me veo de nuevo arrastrado a lanzar
un nuevo libro de Biofísica.
Para llegar a este punto, he debido, en primer lugar, ceder a las peticiones de la casa Flammarion, que me ha rogado dar
una continuidad o más bien un complemento a los «elementos de Biofísica». Más fuertes y más convincentes han resultado
las presiones de mis amigos A. Aurengo, B. Auvert, F. Leterrier y T. Petitclerc. Su talento, su entusiasmo por la enseñanza, su
gusto por el trabajo bien hecho han vencido mis dudas. Demasiado modestos, pensaban que necesitaban de un «antiguo» para
apoyarles en su empresa. Pidiéndome dicho apoyo, han halagado mi vanidad. En definitiva, cedí... Y no lo he lamentado por lo
agradable que ha sido trabajar con ellos.
Hemos estado de acuerdo en hacer sufrir a la enseñanza de nuestra disciplina una evolución significativa. En primer
lugar, este manual de biofísica quiere ser manifiestamente «bio». La física o la físico-química tradicionales han sido reducidas
a un mínimo y frecuentemente se suponen ya conocidas. Además, entre todas las aplicaciones biológicas, hemos querido
privilegiar aquellas más útiles a la fisiología, a la clínica y a las exploraciones funcionales o anatómicas.
No pensamos que esta determinación de encaminarse a las aplicaciones conlleve una devaluación de la enseñanza. La
lectura de algunos capítulos convencerá fácilmente al escéptico de que las ciencias aplicadas reclaman tanta inteligencia
como las fundamentales.
François GRÉMY
Imágenes analógicas
y digitales 20
La imagen médica nació en 1895 a raíz de los descubri- Señal física, imagen analógica e imagen visual
mientos de Roentgen, quien realizó la primera radiografía
utilizando la mano de su esposa. El progreso conseguido En el caso que nos ocupa, una imagen analógica es la
en un siglo descansa no sólo en avances tecnológicos, con traducción de las variaciones espaciales de una señal física
la generación de la gammacámara, la ecografía, el escáner, bajo una forma accesible a la observación visual. Para una
la resonancia magnética (IRM) y la tomografía por emisión imagen de dimensiones horizontal h y vertical v, respec-
de positrones (PET), sino también en el uso de la informá- tivamente, la señal física es una función de dos variables
tica. x ∈ [0,h], y ∈ [0,ν] ⎯→
S
⎯ S(x,y). Por ejemplo, en radiología clá-
Durante mucho tiempo, las imágenes se han obtenido sica, S(x,y) es el producto, llamado exposición, de la intensi-
de forma exclusivamente «analógica», en la pantalla fluo- dad del haz de rayos X por el tiempo de «posado», después de
rescente de los aparatos de radioscopia o sobre las películas atravesar los tejidos a radiografiar.
radiológicas. La informatización de la imagen ha producido La función S(x,y) puede ser caracterizada por sus valo-
imágenes «digitales» constituidas por tablas de cifras a las res mínimo Smín y máximo Smáx entre los cuales S(x,y) puede
que se debe asignar, a priori, un significado visual indepen- tomar un número infinito de valores. Se define el contraste
diente, del modo en el que la imagen numérica ha sido ob- de señal entre dos puntos P1 y P2, en los cuales la señal tiene
tenida. valores S1 y S2, como la razón:
Presentaremos las características principales de las imá- S1 − S2
genes analógicas, y a continuación las de las imágenes digi- CS = (20-1)
tales, explicando el modo de obtención, así como las venta- S1 + S2
jas e inconvenientes con respecto a las imágenes analógicas. La imagen analógica, si se supone de idénticas dimensio-
El tratamiento de las imágenes digitales es el objeto del ca- nes, se define por una función x ∈ [0,h], y ∈ [0,ν] ⎯→ A
⎯ A(x,y).
pítulo 26. En el caso de una película radiológica, A(x,y) representa la
transparencia de la película que debe ser situada sobre una
superficie luminosa uniforme. En el caso de una fotografía o
Características de las imágenes de papel térmico, A(x,y) representa el coeficiente de reflexión
analógicas de la luz que ilumina la imagen. Los valores máximo y míni-
mo de A(x,y) se representan Amín y Amáx y el contraste se define
Las imágenes analógicas de utilidad en Medicina son bien como para la función S(x,y) por:
imágenes reales, obtenidas mediante observación directa
A1 − A 2
(p. ej., en dermatología) o mediante técnicas endoscópicas, CA = (20-2)
o bien imágenes que traducen la intensidad de una señal fí- A1 + A 2
sica (p. ej., de un haz de rayos X o ultrasonidos). Limitaremos Finalmente, la imagen visual percibida por nuestro ojo
nuestro estudio al segundo caso, ya que el primero recoge está caracterizada por su luminancia en cada punto L(x,y). En
técnicas fotográficas en color, no específicas de la imagen radiología, escáner e IRM, se utilizan casi exclusivamente imá-
médica. genes en los niveles de gris y L(x,y) corresponde a la luminancia
273
de la película visualizada en un negatoscopio. A lo largo de este La función de dos variables F(u,v) se llama función de
capítulo, supondremos que se obtiene siempre de esta manera. dispersión puntual, o PSF (del inglés point spread function).
Así pues, en todo punto, la luminancia L(x,y) es proporcional a Por definición, el producto de convolución es la integral si-
la brillancia energética local B(x,y) (véase capítulo 12): guiente:
L(x,y) = k × V × B(x,y) (20-3) S(x,y) F = ∫∫DS(x – u,y – v).F(u,v).du.dv (20-7)
siendo V el coeficiente medio de eficiencia luminosa de la luz Esta formulación traduce el hecho de que el valor de
blanca utilizada. Por tanto, para la imagen visual, el contras- A en el punto (x,y) es una suma ponderada, calculada a par-
te puede ser definido utilizando indistintamente la razón de tir del valor de S en diferentes puntos, siendo el factor de
luminancias o la de brillancias energéticas: ponderación F(u,v). Por ejemplo, S(x,y) interviene con
el «peso» S(x,y) � F(0,0). Del mismo modo, un punto si-
L1 − L 2 B1 − B2
CL = = (20-4) tuado en la unidad de longitud a la izquierda de (x,y), por
L1 + L 2 B1 + B2 tanto de coordenadas (x – 1,y), interviene con el peso
El número de valores diferentes que puede tomar L(x,y) S(x – 1,y) � F(1,0)
es infinito, pero el número de valores percibidos como dife- Teóricamente, la integral doble (20-7) debe ser extendida
rentes por un observador no lo es, puesto que sólo distingui- a todo el espacio. En realidad, sólo los puntos más próximos
mos entre 200 a 250 tonalidades de grises. de (x,y) influyen sobre el valor A(x,y). Cuando |u| o |v| sobre-
pasan un determinado umbral, el punto (x – u,y – v) está
alejado del (x,y) y su influencia sobre A(x,y) es despreciable,
Correspondencia entre señal física por tanto F(u,v) = 0.
e imagen analógica La función de dispersión puntual F(u,v) presenta fre-
En una primera aproximación, a cada valor de S le co- cuentemente una simetría circular; sólo depende de la dis-
rresponde un valor de A, según la siguiente función f: tancia (u,v) en el origen ( )
u2 + v 2 .
f
S(x,y) ⎯→
⎯ A(x,y) = f(S(x,y)) (20-5) Por consiguiente, se puede escribir F(u, v) de la siguiente
En general, la función f no es lineal (no es del tipo manera:
f(x) = ax + b) y presenta fenómenos de umbral y de satura- F(u,v) = G ( )
u 2 + v 2 = G(ρ) , siendo ρ = u2 + v 2 (20-8)
ción. Cuando la señal es demasiado débil, inferior a un um-
bral Sinf, no puede ser detectada y, más allá de un valor límite El esparcimiento o dispersión de F se caracteriza enton-
Ssup, se produce un fenómeno de saturación (figura 20-1). ces por la «anchura a media altura» de la función G(ρ), que
Esto puede formularse: aquí denominaremos FWHM (del inglés full width at half
maximum) (figura 20-2).
S Sinf; f(S) Amín y S Ssup; f(S) Amáx
La función F produce un «desenfoque» de la función A,
Para dos puntos en los que se cumplan las condiciones tanto más acusado cuanto más esparcidos sean los valores
S Sinf o S Ssup, el contraste de la señal se conserva, pero no de F (se tienen en cuenta los valores más alejados) y mayores
el contraste de la imagen analógica. sean los valores de la FWHM de G(ρ). Como muestra la figu-
La correspondencia punto por punto, simbolizada por ra 20-3, este desenfoque se traduce en una mala restitución
la función f, sólo es aproximada. En realidad, debido a las de detalles y contornos de S sobre la imagen A.
imperfecciones de los sistemas físicos, A(x,y) no depende
únicamente de S(x,y), sino también de los valores que toma
S para los puntos cercanos de (x,y). Esta influencia de los
Función de transferencia de modulación
puntos vecinos se representa en forma de un producto de La degradación en el detalle de las informaciones conte-
convolución . nidas en S en el momento del paso a la imagen analógica A se
A(x,y) = f(S(x,y) F) (20-6) caracteriza enteramente por las dos funciones f (correspon-
A
Amáx
Amín
Sinf Ssup S
Figura 20-1. Correspondencia no lineal f, en la que se puede apreciar un fenómeno de umbral y de saturación.
274 IMÁGENES
G(ρ) 1 FTM
0.8
0.6
0.4 A
B C
FWHM 0.2
0 ϕ
0 50 100 150 200
Figura 20-2. Anchura a media altura (FWHM) de una función Para valores suficientemente pequeños como para que
de dispersión puntual y simetría circular. no se produzca saturación, la relación entre las amplitudes
de oscilación Aamp y Samp no depende más que de . Esta rela-
ción se denomina función de transferencia de modulación
� = (MTF, del inglés modulation transfer function)
A amp
MTF(ϕ ) = (20-9)
Samp
La MTF traduce la fidelidad con la que se conservan so-
� =
bre la imagen A las fluctuaciones, y por tanto los contrastes,
de la señal S. Si la MTF tiene un valor elevado, las fluctuacio-
nes correspondientes de S son claramente perceptibles en
A. Por el contrario, un valor bajo de MTF corresponde a una
� =
atenuación de las fluctuaciones de S que consecuentemente
no son visibles en A.
S G S�G Un ejemplo de MTF se representa en la figura 20-5. La
Figura 20-3. Convolución de una misma función S(x,y) por funciones
MTF disminuye con , lo que se traduce en la caída del con-
traste de la imagen cuanto mayor es la resolución de las os-
de simetría circular, F(u,v) = G ( )
u 2 + v 2 , de anchura a media altura
cilaciones (figura 20-6). Sobre la imagen, esta degradación
(FWHM) creciente.
conlleva un rendimiento deficiente, e incluso la desapari-
dencia de niveles) y F (degradación espacial). Se puede igual- ción de los detalles de pequeño tamaño y la atenuación de la
mente representar esta degradación examinando el modo en nitidez de los contornos.
el que un detalle de S se traduce en la imagen A.
Sea una señal S(x,y) formada por bandas verticales uni- Películas radiológicas analógicas
formes, con una variación sinusoidal de media Smed y ampli-
Actualmente, las películas radiológicas son el único tipo de
tud Samp en el sentido horizontal, o sea S(x,y) = Smed + Samp.
imagen analógica que traduce una señal física utilizada en ima-
sen(2πx) (figura 20-4). El parámetro (expresado en «osci-
gen médica; las otras imágenes, reproducidas sobre película o
laciones/cm») caracteriza la frecuencia espacial de imágenes
sobre papel de impresión apropiado, son obtenidas directamen-
cuyas oscilaciones son tanto más próximas cuanto mayor es .
te de modo digital o por digitalización de una señal de vídeo.
El período espacial de S es λ = −1 (cm).
El principio de la detección de las radiaciones por una
Se puede mostrar que la imagen A formada a partir de S
emulsión fotográfica y por una pantalla de cristal fluorescen-
tiene la siguiente formulación:
te se ha estudiado en el capítulo 16 (véase Emulsiones foto-
A(x,y) = Amed + Aamp sen (2πϕx) gráficas y Pantallas de cristal fluorescente).
Una película radiológica está compuesta de 4 ó 5 capas
principales (figura 20-7). Las películas «monocapa» se compo-
nen de una única capa de emulsión; se utilizan con una panta-
lla intensificadora que absorbe los rayos X y emite fotones de
fluorescencia que son los responsables de la exposición de la
λ emulsión de la película. Las películas «bicapa» se componen
de dos capas de emulsión, lo cual duplica su sensibilidad, y se
colocan entre dos pantallas intensificadoras fluorescentes.
Para retomar las anotaciones de los párrafos anteriores:
– la señal S(x,y) se denomina exposición. Es la energía
Figura 20-4. Señal de variaciones sinusoidales (λ = ϕ–1). depositada por el haz de rayos X tras atravesar la región ra-
pantalla intensificadora
Figura 20-7. Película radiográfica monocapa (A) y bicapa (B). La capa antiestrés antiestático protege la emulsión y permite separar
fácilmente la película de la pantalla sobre el que está puesta. La capa anti-halo evita el reflejo de la luz que produciría una imagen
con halo en la emulsión.
276 IMÁGENES
3 c
DO
saturación 1 2 3 C
2.5
1
2 pendiente γ 2
1.5
3
1 umbral
0.5
0
log(E)
Memoria ocupada por una imagen digital tiene como dimensiones N líneas y N columnas, cada
píxel representa una región cuadrada del espacio de lado
Si R es la resolución de la intensidad expresada en octetos ∆ = z/N, y el valor del píxel I(l,c) debe corresponder a los valo-
por píxel, el volumen de memoria ocupado por una imagen res de S en la región delimitada por los intervalos [∆ (l – 1),
N N es de N2 R octetos. Una exploración está frecuente- ∆ 1] de altura y [∆ (c – 1), ∆ c] de ancho. Se muestra que
mente constituida por una serie de P imágenes (p. ej., de 10 el valor I(l,c) se obtiene por convolución de los valores S(x,y) a
a 20 cortes anatómicos en IRM o una secuencia temporal de través de una función denominada de muestreo M(u,v):
120 imágenes en centellografía), que ocupan M = N2 R P
I(l,c) = ∫∫ DS(∆ (c 1/2) u,∆
octetos en la memoria.
(l + 1/2) v).M(u,v).du.dv
Para un muestreo perfecto, M(u,v) debería ser uniforme
Obtención de imágenes digitales en un cuadrado centrado en el origen y de lado ∆, y nulo fue-
ra de dicho cuadrado. En realidad M(u,v) no es uniforme en
En algunas modalidades, la imagen se obtiene por un el cuadrado asociado a un píxel y puede extenderse a las zo-
cálculo realizado sobre datos obtenidos por digitalización nas correspondientes a píxeles vecinos.
de una señal. Es el caso de las imágenes de reconstrucción En el caso de la imagen centellográfica (véase capítu-
tomográfica (escáner, IRM, tomo-centellografía) y las imá- lo 25), la gammacámara aporta la abscisa X y la ordenada Y
genes de sustracción en angiografía digital. Estas imágenes de cada fotón γ detectado. La imagen digital se obtiene
denominadas «intrínsecamente digitales» no derivan direc- adscribiendo dicho fotón al píxel (l,c) de tal manera que
tamente de una señal S(x,y). ∆ (c – 1) < X ∆ c y ∆ (l – 1) < Y ∆ l (el valor de I(l,c)
Por el contrario, en el caso de la radiología digitalizada se incrementa en una unidad).
(véase capítulo 22), la imagen digital se obtiene a partir de Una imagen digital puede también obtenerse por digita-
una señal S(x,y). La correspondencia entre los valores S(x,y) lización de una imagen analógica, por medio de un escáner
de la señal y los valores I(l,c) de la imagen digital, se denomi- de documentos. En este caso, los valores de la imagen digi-
na muestreo de S(x,y). Si las variaciones de la señal se hacen tal vienen dados por una ecuación de convolución similar
en una región cuadrada del espacio de lado z, y si la imagen a (20-8).
278 IMÁGENES
Visualización de una imagen digital El aspecto de la imagen obtenida depende en gran me-
dida de la función elegida y el contraste visual de la imagen
Imágenes aisladas será máximo para la gama de valores v para los cuales la pen-
diente de la función G es más elevada:
En el caso de una imagen de escáner, de IRM o de cen- – la función más sencilla es lineal, asocia el color negro
tellografía, el color de los píxeles de una imagen digitales es (G = 0) al valor I(l,c) = 0 y el blanco (tomaremos G = 100) al
una convención arbitraria, cuya meta es poner en evidencia, valor de I(l,c) = 255 en una correspondencia lineal entre es-
visualmente, las características pertinentes de la imagen. Se tos dos extremos: G(ν) = ν (100/255) (figura 20-12);
dice que la imagen obtenida está en pseudocolor. La prime- – una función de segmentación consiste en fijar un um-
ra etapa de la visualización de una imagen digital de este tipo bral Smín por encima del cual todos los píxeles serán repre-
consiste en la elección del color que haremos corresponder sentados en negro y un umbral Smáx por debajo del cual lo se-
con cada valor I(l,c) de la imagen. Nos limitaremos a la si- rán en blanco, con una correspondencia lineal en el intervalo
tuación en la que todos los colores elegidos se correspondan (figuras 20-13 y 20-14):
a diferentes niveles de gris; la generalización a otros colores ⎧⎪ν < Smín ⇒ G(ν) = 0
cualesquiera será inmediata. ⎪⎪
⎪⎪ ν − Smín
Para fijar los conceptos, supongamos que la imagen di- ⎨ ν ∈ [Smín , Smáx → G(ν) = ×100
gital está almacenada en 1 octeto por píxel y que toma va- ⎪⎪ Smáx − Smín
⎪⎪
lores comprendidos entre 0 y 255. Un nivel de gris se define ⎪⎩ν > Smáx ⇒ G(ν) = 100
por el valor idéntico de sus tres componentes rojo, verde y
La segmentación permite poner en evidencia partes es-
azul (RGB, del inglés red, green, blue) que denominaremos G.
pecíficas de la imagen y favorecer, con un contraste amplifi-
Supondremos que los valores posibles de G están compren-
cado, la visualización de píxeles cuyo valor está comprendi-
didos entre 0 (negro) y 100 (blanco).
do entre [Smín–Smáx];
Debemos, por tanto, elegir una función que proporcio-
ne un nivel de gris para cada valor de la imagen: ν = I(l,c)
[0.255] → G(ν).
nivel de gris
nivel de gris 100
100
0 0
0 255 0 Smáx Smín 255
valor numérico de la imagen valor numérico de la imagen
80
γ = 0.4
60
γ = 1.0
40
γ = 1.5
20
V
0
0 50 100 150 200 250
280 IMÁGENES
La representación visual puede consistir simplemente en las películas radiológicas que permiten conservar dichas
representar las imágenes unas al lado de las otras. En el caso imágenes.
de series espaciales, la imagen digital permite una represen- El archivo de las imágenes digitales es conceptualmente
tación global, tridimensional, de la información contenida sencillo y perfectamente adaptado a una gestión informati-
en varios cortes anatómicos, bajo la forma de imágenes 3D. zada en la cual la imagen no es más que un elemento más
Las técnicas utilizadas para realizar imágenes 3D serán ex- de la historia clínica. Las imágenes pueden ser indexadas en
puestas en el capítulo 26. En el caso de series dinámicas, se una base de datos y son fácilmente localizables según una
puede colocar en la pantalla de un ordenador, sucesivamen- gran diversidad de criterios, como por ejemplo, encontrar
te, con una cadencia rápida, las imágenes de la serie, produ- entre 12 000 pruebas todas las angiografías digitales realiza-
ciendo la impresión de movimiento. Otras técnicas de cálcu- das entre 1991 y 1993 en pacientes de menos de 40 años y
lo y visualización de la información contenida en una serie que muestren una malformación de la arteria silviana.
de imágenes serán estudiadas en el capítulo 26. La imagen multimodal y la transferencia de imágenes
plantean el problema de la multiplicidad de formatos de los
archivos informáticos que contienen imágenes. Debido a la
Ventajas de la imagen digital presión ejercida por los usuarios, el formato DICOM se ha
convertido en el estándar universalmente aceptado en ima-
Con respecto a las imágenes analógicas, las imágenes di- gen médica, lo que permite transmitir las imágenes entre to-
gitales tienen, en primer lugar, la ventaja de la universalidad: das las instalaciones de imagen. Este formato permite codifi-
que se trate de una fotografía personal, de una IRM o de una car no sólo la imagen, sino también la información asociada
prueba de centelleo, la imagen puede ser conservada en el (identificación del paciente, fecha, hora, lugar y modalidad
mismo tipo de soporte (disco magnético u óptico), visuali- de la prueba, etc.).
zado en el mismo tipo de pantalla de ordenador, reproduci- El transporte de la imágenes digitales es, teóricamente,
do por la misma impresora. Esta característica permite una sencillo: la red de Internet, por ejemplo, permite transferir
imagen multimodal, superponiendo sobre la misma imagen imágenes digitales a cualquier lugar del globo, si remiten-
informaciones obtenidas con técnicas diferentes (p. ej., escá- te y destinatario disponen de dirección electrónica y de un
ner + centelleo). formato de archivo común. Es posible consultar en Internet
El dinamismo de las imágenes digitales es intrínseca- varias y voluminosas bases de imágenes médicas digitales
mente más amplio que el de las imágenes analógicas. Debi- (p. ej., http://www.med.harvard.edu/AANLIB/home.html).
do a las limitaciones visuales, en una imagen analógica sólo
se pueden representar de 30 a 60 niveles de grises distintos.
Sin embargo, una imagen digital codificada sobre dos octe- Inconvenientes de la imagen digital
tos permite representar 65 536 niveles diferentes. No pueden
ser visualizados simultáneamente, pero sí a intervalos, gra- Con respecto a una imagen analógica, una imagen digi-
cias al método de segmentación. tal posee, a priori, una resolución más grosera, a menos que
La flexibilidad de visualización es una ventaja conside- se utilicen imágenes de muy elevada resolución espacial. Por
rable de las imágenes digitales, para las que los valores de ejemplo, una imagen digital de radiografía torácica de una
la señal física y la visualización propiamente dicha son dos calidad comparable a una película debería tener una resolu-
conceptos independientes. Una misma información puede ción 2048 2048, lo que requiere pantallas de visualización
ser presentada de diversas maneras según el objetivo de la especiales. La evolución rápida del material informático per-
prueba. La segmentación, por ejemplo, permite examinar mite progresivamente resolver problemas planteados por la
una misma prueba de escáner mostrando los detalles de las manipulación de gran número de imágenes de gran tamaño.
estructuras óseas, o bien los de los tejidos blandos. Esta flexi- La transferencia de imágenes a través de la red telefónica
bilidad se extiende a representaciones más complejas, como estándar es bastante rápida con las técnicas ADSL o el cable,
pueden ser las tridimensionales. pero las transferencias masivas en el interior de un hospital
La imagen digital permite una aproximación cuantitativa plantea problemas técnicos debidos a la velocidad limitada
de la imagen determinando, en ciertas áreas de la imagen, pa- de las redes, lo que hace necesario contar con instalaciones
rámetros pertinentes (p. ej., las dimensiones, la superficie y el específicas y costosas.
volumen de un órgano, la densidad media de un tejido). Estas Las imágenes pueden ser comprimidas antes de su envío
informaciones permiten una interpretación más objetiva de y descomprimidas a su llegada, pero por razones de respon-
las imágenes y son prácticamente imposibles de obtener con sabilidad médica, solo se pueden utilizar algoritmos de com-
imágenes analógicas. De un modo más general, las imágenes presión «conservadores», que den lugar a una restitución ad
digitales permiten numerosos tratamientos informáticos. integrum de la imagen, una vez descomprimida. Dichos al-
La conservación de las imágenes digitales es muy fiable, goritmos conservadores son mucho menos eficaces (tasas de
duradera, y poco costosa. Por ejemplo, un disco óptico digi- compresión de 2 a 3) que los algoritmos que conllevan una
tal de 600 Mb permite conservar 4500 imágenes 256 256 cierta degradación de las imágenes (tasas de compresión de
codificadas sobre dos octetos y cuesta 100 veces menos que hasta 30).
Ejercicio 20-1. Entre dos puntos de una imagen física en la Ejercicio 20-5. ¿Cómo se debe elegir el valor de γ en una fun-
cual la señal tiene valores S1 y S2, existe un contraste de señal ción de correspondencia, para acentuar el contraste en las
Cs. ¿Como variará Cs si las señales se multiplican por un mis- zonas de señal débil? ¿Y en las de señal fuerte?
282 IMÁGENES
Tomografía
computarizada 21
283
Y Retroproyección filtrada continua
La retroproyección de una proyección concreta θ(ρ),
asociada a una recta Dθ, consiste en «esparcir» sobre una
Dθ imagen inicialmente virgen, un valor igual a θ(ρ) para todos
los puntos que se proyectan sobre Dθ en el punto R definido
R
O ρ θ por OR = ρ (figura 21-4). Si se esparce θ , para todos los
X valores de θ comprendidos entre 0 y 2π, y con cada retropro-
yección que interviene con un peso infinitesimal dθ, se ob-
tiene una imagen B(x,y) llamada imagen de retroproyec-
W ción, que reproduce los principales aspectos de la imagen
inicial A(x,y) modificada por un desenfoque importante, que
borra los contornos (figura 21-5).
�θ(ρ) ρ
Se muestra que B(x,y) es la transformada de la imagen
inicial A(x,y) por convolución con una función de dispersión
puntual en l/r, donde r indica la distancia al origen referen-
cial (figura 21-6). Con las anotaciones explicadas en la pági-
na 274, y con r = u2 + v 2 :
Figura 21-2. Transformada de Radon de una imagen analógica. B(x,y) = A(x,y) (l/r) =
1
= ∫∫ A(x − u,y − v). .du.dv (21-2)
D
u2 + v 2
imagen analógica
Y La retroproyección permite, por tanto, reconstruir la
imagen B(x,y), próxima de A(x,y), a partir de proyecciones.
Para que la restitución sea perfecta, es necesario corregir el
efecto de la función de dispersión en l/r y por tanto invertir
X proyección el efecto de la convolución (l/r).
� � 90°
Esta corrección se obtiene a través de la operación inver-
sa de la convolución, llamada deconvolución. El principio se
proyección basa en utilizar la transformada de Fourier de las imágenes
� � 0° A, B y de la función l/r. La transformada de Fourier de una
imagen, presentada en el anexo 9, sobrepasa el marco de este
libro. Digamos, simplemente, que esta operación, invertible,
sinograma transforma una imagen A(x,y) en una función de dos varia-
�
bles a y b denominada F(A). Los valores de F(A) en dos puntos
próximos al origen traducen variaciones suaves (frecuencias
bajas) de A(x,y), y corresponde a zonas relativamente unifor-
mes de la imagen. Inversamente, los valores en puntos aleja-
� dos del origen traducen variaciones bruscas (frecuencias al-
tas), y corresponden a los detalles y los bordes de los objetos.
Figura 21-3. Sinograma de una imagen analógica
La transformada de Fourier tiene la propiedad de transfor-
mar la convolución en una sencilla multiplicación; por tanto,
tomando la transformada de Fourier de (21-2), se obtiene:
F(B(x,y)) = F(A(x,y) (l/r)) = F(A(x,y)) � F(l/r)
llamada sinograma, con ρ en abscisas y θ en ordenadas
(figura 21-3). Pero la transformada de Fourier de l/r es F(1/r) = 1/ a 2 + b2 =
F(1/r) = 1/ a 2 + b2 = 1/d (en la cual d indica la distancia al origen
Se demuestra que el conjunto de las transformadas de Ra-
don de una imagen A(x,y), es decir, su sinograma, contiene en el plano a, b), por lo que:
todas las informaciones necesarias para reconstruir A(x,y). F(B(x,y)) = F(A(x,y))/d
Dicho de otro modo, se puede reconstruir exactamente la es decir,
imagen inicial, si se conocen todas sus proyecciones θ(ρ), F(A(x,y)) = F(B(x,y)) � d
para cualquier valor de θ y de ρ. Esta demostración descan-
y, aplicando la transformación de Fourier inversa F–1:
sa sobre el teorema del perfil central de la transformada de
Fourier, cuyo principio se explica en el anexo 10. La deter- A(x,y) = F–1(F(B(x,y)) � d)
minación de A(x,y) a partir de θ(ρ), denominado problema La transformada de Fourier proporciona, de esta mane-
inverso, utiliza el método de retroproyección filtrada. ra, una herramienta teórica para reconstruir perfectamente
284 IMÁGENES
Y Y
Y
O F
X X X
Y
Retroproyección Retroproyección
en una dirección en tres direcciones
Y
Imagen inicial
286 IMÁGENES
1
0.8 1
0.6
0.4
2
0.2
3
4
0
0 20 40 60 80 100
Figura 21-11. Los filtros de reconstrucción y sus resultados. En abscisas: distancia d al centro en porcentaje de la distancia de corte.
En ordenadas: (d).
Calidad de la imagen reconstruida te elevado para utilizar los métodos clásicos de resolución de
ecuaciones lineales (65 536 para una imagen 256 � 256), y
La calidad de la imagen reconstruida depende, eviden-
que tiene que recurrir a métodos iterativos.
temente, del detector que proporciona las proyecciones dis-
El principio de los métodos iterativos se basa en partir
cretas. Depende igualmente del número de proyecciones gra-
de una imagen digital a priori (p. ej., uniforme) y calcular sus
badas (para reconstruir de manera satisfactoria una imagen
proyecciones en las mismas direcciones que la adquisición
N � N, teóricamente hay que obtener al menos πN/2 proyec-
real. Se comparan a continuación estas proyecciones calcu-
ciones). La degradación de la calidad de la reconstrucción,
ladas (simuladas) con las proyecciones obtenidas realmente
cuando el número de proyecciones disminuye, se ilustra en
y se modifica la imagen digital, de manera que se disminuya
la figura 21-12. La precisión de los cálculos y la naturaleza
la distancia entre la simulación y la realidad. Este procedi-
del filtro utilizado son, igualmente, elementos importantes.
miento se repite (de ahí su denominación «método iterati-
Se han propuesto numerosas mejoras de los algoritmos que
vo») hasta que la diferencia entre las proyecciones simuladas
hemos presentado.
y las reales sea suficientemente pequeña.
Los métodos iterativos requieren potencias de cálcu-
lo importantes, fácilmente disponibles en los ordenadores
Aproximación algebraica actuales. Tienen la gran ventaja de poder considerar de ma-
nera sencilla, en la simulación, todas las particularidades del
Se considera, en este segundo método, que la ecuación modo de adquisición (atenuación, correcciones geométri-
(21-3) es una relación lineal entre las N2 desconocidas I(l,c). cas, características del detector).
Se pueden calcular teóricamente todas estas incógnitas y
por consiguiente determinar exactamente la imagen digital
I(l,c), a condición de que se disponga de, como mínimo, un Reconstrucción en otras direcciones
mismo número de relaciones independientes que de incóg- en el espacio
nitas. Ahora bien, conocemos T � N valores de Fθ(i), por lo
que el problema tiene teóricamente solución, si T � N N2, Hemos descrito el principio de la reconstrucción digital
esto es, si T N. Esta aproximación se enfrenta con el hecho de cortes tomográficos transversales. A partir de una serie de
de que el número de incógnitas del sistema es excesivamen- cortes transversales ensamblados, reconstrucciones infor-
Imágenes filtradas
Figura 21-12. Calidad de la reconstrucción según el número de proyecciones. La imagen inicial (véase figura 21-9) ha sido proyectada
y posteriormente reconstruida por retroproyección y posterior filtrado por filtro de Hamming, para 4, 8, 16 y 32 proyecciones. La visualización
utiliza una correspondencia γ = 0.7 para poner mejor de manifiesto los errores de reconstrucción.
máticas permiten representar un corte en cualquier direc- La imagen de escáner, estudiada en el capítulo 22, per-
ción del espacio. Para un corte sagital, por ejemplo, se coloca mite establecer una cartografía tomográfica de los coeficien-
en una columna todas las imágenes transversales, una al lado tes de atenuación lineal de los tejidos para los rayos X. El
de la otra, obteniendo una nueva imagen (figura 21-13). principio del escáner de primera generación está más próxi-
mo a los conceptos generales que hemos estudiado en los
apartados anteriores. Los «escanógrafos» actuales utilizan
Aplicaciones a la imagen médica una batería de detectores y la «proyección» no es ortogonal,
sino en abanico. El algoritmo de reconstrucción es más com-
Las técnicas de reconstrucción tomográfica constituyen la
plicado que con la metodología de retroproyección filtrada,
base de tres modalidades fundamentales de la imagen médica.
aunque se fundamenta sobre los mismos conceptos: abarca
una etapa complementaria de formateado de los datos y uti-
liza filtros adaptados y una proyección en abanico.
La tomografía de centelleo, que se estudia en el capítu-
lo 25, permite construir una imagen tomográfica de la distri-
bución de fuentes radiactivas en el interior del organismo. La
formulación del problema inverso (construir la imagen co-
nociendo sus «proyecciones») es más compleja que la trans-
formada discreta de Radon, en razón de la disminución del
rendimiento de contaje de fuentes radiactivas cuando nos
alejamos del detector y los fenómenos de atenuación de la
emisión γ. Estas dificultades se resuelven más fácilmente por
técnicas algebraicas que por la aproximación de retropro-
yección filtrada. Ciertos colimadores llevan a cabo proyeccio-
nes cónicas o en abanico. En este caso, nos encontramos con
el mismo problema de reconstrucción que con los escanó-
grafos actuales y con el mismo tipo de método de recons-
Figura 21-13. Reconstrucción en otras direcciones en el espacio.
Las líneas homólogas de las secciones tomográficas transversales
trucción.
de una bifurcación vascular se reagrupan para formar una imagen En el caso de la IRM, la formulación del problema de
tomográfica frontal. reconstrucción de la imagen puede abordarse en el marco
288 IMÁGENES
general de la transformada de Radon, pero lo más frecuente Ejercicios
es que se utilice la codificación múltiple (véase capítulo 23)
que no contempla la reconstrucción tomográfica estudiada Ejercicio 21-1. Calcúlese la transformada de Radon θ(ρ)
en este capítulo. de la imagen analógica A(x,y) definida por: A(x,y) = 1 si x2 +
+ y2 1; A(x,y) = 1 si x2 + y2 1. Se demostrará que θ(ρ) no
depende de θ y que puede ser calculada para θ = 0°.
Límites teóricos y prácticos Ejercicio 21-2. Una imagen analógica es simétrica con res-
pecto al eje OY. ¿Tendrá su sinograma necesariamente un eje
de la imagen tomográfica de simetría?, y si es así, ¿cuál será dicho eje?
computarizada
Ejercicio 21-3. Demuéstrese que si una imagen tiene una
simetría de revolución (es decir, si es invariante por rotación
A pesar de sus éxitos espectaculares, la imagen tomo-
de un ángulo cualquiera alrededor del centro), su transfor-
gráfica computarizada tiene ciertos límites teóricos y prác-
mada de Radon no depende de θ. ¿Qué se puede deducir en
ticos:
cuanto a su sinograma? (Utilizando la intuición, este ejerci-
– el número de proyecciones y de puntos de adqui-
cio es fácil, pero su demostración rigurosa resulta difícil)
sición es limitado, lo que implica un límite teórico de re-
solución espacial de las imágenes que se pueden obtener. Ejercicio 21-4. Se proporciona la siguiente imagen digital
Este límite se va haciendo más pequeño con el progreso 4 � 4:
de detectores y de la informática, alcanzando la resolu-
2 5 2 1
ción espacial de los actuales escanógrafos alrededor de
0.3 mm; 5 30 5 1
– los cálculos son necesariamente aproximados, inclu- 2 5 2 1
so hechos con una precisión importante (errores de redon- 1 1 1 0
deo);
– el espectro de Fourier es limitado, con errores numé- Calcúlese su transformada de Radon en las direcciones
ricos de truncado (que llevan a oscilaciones y fenómenos de θ = 0°, θ = 90° y θ = 180°. Determinar los valores de la figura
«tornasoleado» de la imagen reconstruida); de retroproyección de las transformadas θ = 0° y 90°.
– las fuertes interrupciones de la señal (p. ej., en el escá- Ejercicio 21-5. Se quiere determinar la transformada de Ra-
ner, en caso de prótesis metálica) no son tenidas en cuenta don de la imagen del ejercicio 21-4 en la dirección θ = 45°.
de manera satisfactoria por las técnicas de filtrado que tra- Demostrar que el píxel de valor 30 (l = 2, c = 2) afecta los
bajan sobre un espectro limitado de frecuencias. Producen valores de F45(2) y F45 (3). ¿Con qué coeficientes se afectará
artefactos que degradan considerablemente la imagen re- el valor de dicho píxel en el cálculo F45 (2) y F45 (3)? (Estos
construida; coeficientes son H245(2,2) y H345(2,2); véase pág. 286).
– el muestreo no es uniforme debido a los fenómenos
de difusión y de atenuación de los rayos X o γ ; Ejercicio 21-6. De los diferentes filtros rampa, de Hamming,
– finalmente, los movimientos circulatorios y respirato- de Hanning o de Parzen: ¿cuál le parece más adecuado para
rios del paciente no son despreciables y son, de hecho, fuen- atenuar el ruido presente en una imagen? ¿Cuál es el más
te de artefactos. apropiado para obtener los contornos más nítidos?
La imagen radiológica se basa en el hecho de que los di- – interacciones con los núcleos de la diana que produ-
ferentes tejidos del organismo atenúan los rayos X de forma cen la radiación de frenado, cuyo espectro energético es
distinta según su composición. Cuando un haz homogéneo continuo;
de rayos X atraviesa el cuerpo humano, el haz emergente – interacciones con los electrones de la diana que llevan
exhibe cambios de homogeneidad que reflejan la forma y a la radiación característica cuyo espectro discontinuo, de
la naturaleza de los tejidos atravesados. La información re- rayas, es característico de la diana.
sultante, que constituye la denominada imagen radiante, En las fuentes utilizadas en radiodiagnóstico, los rayos X
puede ser utilizada para formar una imagen con un detector son emitidos en un dispositivo llamado tubo de rayos X,
adecuado. Se trata pues de una imagen obtenida por atenua- cuyo objetivo es obtener un haz de electrones homogéneo,
ción, y que está esencialmente representada por la radiogra- de energía ajustable, de modo que la dispersión de energía
fía estándar, el escáner y la radiología digital por sustracción. asociada a los fotones emitidos sea la menor posible, alrede-
En el momento actual, las principales evoluciones se orien- dor de un valor medio igualmente ajustable.
tan hacia una imagen enteramente digital y el desarrollo de El tubo de rayos X está constituido por un emisor de
las técnicas intervencionistas. electrones, un acelerador de electrones y una diana metáli-
Las imágenes son obtenidas pagando el precio de ca. Estos componentes se insertan dentro de una ampolla de
una irradiación del paciente y del personal que realiza el exa- vidrio transparente a los rayos X, y en la cual se produce el
men. Según los principios generales de Radio-protección, vacío (figura 22-1).
esta irradiación debe estar siempre médicamente justifi- El emisor de electrones o cátodo, es un filamento
cada, debe respetar las contraindicaciones, y debe ir acom- fino de tungsteno, que es calentado por una corriente que
pañada de todo el esfuerzo técnico razonable para redu- lo lleva a la incandescencia. Este calentamiento arrastra,
cirla. por efecto termo-iónico, la emisión de un flujo de electrones
que aumenta con la intensidad de la corriente de calenta-
miento.
La diana metálica o ánodo es una aleación de tungsteno
Producción de los rayos X y renio. La zona emisora de la diana se llama foco. Un foco
en radiodiagnóstico de área pequeña permite obtener imágenes muy precisas,
pero limita la intensidad del haz y viceversa. Se utiliza una
Los rayos X (RX) son fotones de gran energía en teoría placa o ánodo fijo para la instalaciones de potencia débil, y
superior a una decena, y en la práctica a un millar de elec- una pieza que gira o ánodo rotatorio para obtener potencias
tronvoltios. Los RX se generan cuando un haz de electrones elevadas con un foco de pequeñas dimensiones, evitando así
acelerados en el vacío percuten en una diana material. Se un calentamiento excesivamente importante de la diana (fi-
producen entonces dos tipos de interacciones entre los elec- gura 22-2).
trones acelerados y los átomos que componen la diana (véa- La aceleración de los electrones se obtiene por una alta
se capítulo 15, Interacción entre las radiaciones ionizantes y tensión U de varias decenas de kV, ajustable, aplicada entre
la materia): el cátodo y el ánodo. El flujo de electrones representa una
291
Generador
C A
Flujo
de electrones
Diafragma
RX
Electrones
Rotación rápida
Foco óptico
Rayos X
292 IMÁGENES
dΦ energía está comprendida entre 25 y 130 keV. La energía de
dE los fotones es regulable en función del tipo de examen ra-
KiZU diológico, ajustándose la tensión del tubo. Con una tensión
comprendida entre 25 y 30 kV, se llama sistema de baja ten-
sión. Entre 110 y 130 kV, se le llama de alta tensión.
En un examen radiológico, para una tensión determina-
da del tubo, la energía total emitida por el tubo, así como la
recibida por el paciente, es proporcional al número de elec-
trones que bombardean el ánodo. Este número es proporcio-
E nal al producto iT (mA.s) de la intensidad i (mA) del haz que
U bombardea el ánodo, por el tiempo T, denominado tiempo
de «posado», tiempo durante el cual el sujeto es sometido a
Figura 22-3. Espectro teórico de un tubo de rayos X.
la acción de los rayos X.
La fuente puede estar situada relativamente lejos del su-
jeto puesto que el aire, como cualquier otro gas, no atenúa
En realidad, los fotones de baja energía son absorbidos prácticamente los rayos X. El ánodo emite rayos X en todas
por la diana y el haz es filtrado, para posteriormente pasar las direcciones. El haz es por tanto divergente. Como el foco
por un diafragma para que sólo se irradie la zona del pa- emisor es por lo general muy pequeño con respecto a la su-
ciente que se desea estudiar. El filtrado sirve para eliminar perficie del receptor, la intensidad del haz a una distancia D
los fotones X de muy baja energía que serían absorbidos por de la fuente varía proporcionalmente a 1/D2; doblar la dis-
los tejidos superficiales y aumentarían la irradiación del pa- tancia entre la fuente y el paciente divide por cuatro la inten-
ciente sin valor añadido para la formación de la imagen ra- sidad del haz recibido por éste.
diológica. El rendimiento global real de los tubos de rayos X Intensidad del haz de rayos X y exposición. Se llama in-
es muy bajo. Alrededor del 1 por 1000 de la potencia utilizada tensidad de un haz de rayos X a la energía que atraviesa, por
se emite en forma de rayos X: el resto se disipa en forma de unidad de tiempo, una unidad de superficie perpendicular a
calor. la dirección del haz. Es pues el aclaramiento energético, tal y
La radiación característica del tungsteno y del renio se como lo hemos definido en el capítulo 17 (véase Parámetros
superponen en forma de rayas espectrales al espectro con- energéticos), expresándose en W.m–2.
tinuo de la radiación de frenado, pero su valor es muy pe- Se llama exposición (de la película radiológica por ejem-
queño con respecto a éste. Se hace despreciable cuando los plo), a la cantidad total de energía depositada por el haz
electrones incidentes son muy energéticos puesto que los durante toda la duración del examen. La exposición E es el
fotones de fluorescencia, que son precisamente los que la producto de la intensidad I por el tiempo de duración de la
constituyen, tienen energías bajas (del orden de 70KeV des- prueba T, es decir E = I T; se expresa en J.m–2.
pués de la expulsión de un electrón de la capa K del tungs-
teno). El espectro real de la emisión de un tubo de rayos X
queda esquematizado en la figura 22-4.
Tras la filtración, el haz real puede ser asimilado a un Imagen radiante
flujo de fotones X monoenergéticos cuya energía en keV es
igual a aproximadamente un tercio de la tensión del tubo ex- Formación de la imagen radiante
presada en kV. En general, se utilizan haces de fotones cuya
Cuando un haz monoenergético de rayos X paralelos
atraviesa un medio homogéneo, la intensidad I de los foto-
dΦ nes que han atravesado un espesor x sin sufrir interacción
dE decrece en función de x, siguiendo una ley exponencial:
I = I0 e–µx (22-3)
Serie L
en la cual Io es la intensidad del haz incidente en el punto de
entrada en el medio y µ (en m–1) el coeficiente lineal global de
atenuación (por difusión y por absorción). La diferencia de
Serie K
intensidad Io – I perdida por el haz incidente corresponde a
los fotones incidentes que han interaccionado con el medio
atravesado (figura 22-5).
E Cuando un haz de rayos X de intensidad incidente uni-
forme atraviesa un objeto no homogéneo desde el punto
U de vista de la atenuación de los rayos X, el haz transmitido
Figura 22-4. Espectro real de un tubo de rayos X. ya no es uniforme: su intensidad varía de un punto a otro
I1 − I 2
CR = (22-4)
I0 e–µx I1 + I 2
Esta relación está directamente ligada a las posibilidades
x
Distancia de discriminar visualmente estas dos zonas después del reve-
lado de la imagen radiante. Examinemos la manera de esti-
t mar el contraste de una estructura de espesor uniforme x, y de
d coeficiente de atenuación lineal µ, situada en un medio de es-
t pesor y con coeficiente de atenuación lineal ν (figura 22-7):
a
I1 = I0e–νy; I2 = I0e–ν(y – x)e–µx = I0e–ν(y – x) – µx
Si los exponentes son pequeños con respecto a 1, se pue-
Figura 22-5. Atenuación de un haz de rayos X. t : rayo X transmitido; de utilizar la aproximación eε 1 + ε, y por tanto:
d : rayo X difundido; a : rayo X absorbido. I1 I0 (1 – νy); I2 I0 (1 – ν(y – x) – µx)
y el contraste valdrá:
I1 − I 2 µx − νx 1
CR = µ − ν x (22-5)
de la sección perpendicular a la trayectoria del haz. Se de- I1 + I 2 2 − 2νy + x(ν − µ ) 2
nomina imagen radiante la distribución de la intensidad Para que una estructura produzca un contraste impor-
en la sección recta del haz transmitido (figura 22-6). Esta tante, es necesario pues que el valor absoluto de la diferencia
imagen radiante invisible a nuestro ojo es recogida por el |µ – ν| entre su coeficiente de atenuación y el del medio que
detector, cuyo papel es de transformarla en una «imagen la rodea sea grande, y que su espesor x en el eje del haz sea
visual» visible para el observador (véase capítulo 20, Señal importante.
física, imagen analógica e imagen visual). La presencia de El contraste parece, en este caso particular, indepen-
materia poco radiopaca (aire, por ejemplo) en un medio diente del espesor y del medio que la rodea. En realidad, la
más absorbente se traduce en la imagen radiante por una radiación difundida, que no se tiene en cuenta en el cálculo
zona cuya intensidad es relativamente mayor. En efecto, los anterior, y cuya importancia crece con el espesor atravesado,
rayos X que atraviesan esta materia poco radiopaca están produce una disminución del contraste. En efecto, si a I1 e I2
menos atenuados que los otros. Por el contrario, la presen- se superpone una intensidad difundida D, estas cantidades
cia de un cuerpo muy radiopaco (un hueso o una prótesis serán I1 + D e I2 + D, y el contraste vendrá dado por :
metálica) en un medio menos absorbente (agua o músculo)
corresponde a una zona de la imagen de relativamente me- (I1 + D) − (I 2 + D) I1 − I 2 I1 − I 2
C DR = = = C R (22-6)
nor intensidad. (I1 + D) + (I 2 + D) I1 + I 2 + 2D
D I1 + I 2
haz de haz de
RX incidentes RX transmitidos
1
2 detector
(pantalla luminiscente)
3
IO sujeto
imagen
radiante
1 aire
2 músculo imagen
3 hueso luminosa
294 IMÁGENES
En el efecto fotoeléctrico, prácticamente toda la energía
del fotón es absorbida por el medio. En el efecto Compton,
ν sólo una parte de esta energía es absorbida con emisión de
un fotón X difundido en una dirección al azar. Estos fotones
I0 I1
difundidos pueden atravesar el medio sin sufrir interacción,
o ser ellos mismos difundidos o absorbidos.
La probabilidad del efecto fotoeléctrico es mayor cuanto
µ
mayor es el número atómico Z de los átomos que componen
los tejidos. Así, los huesos, ricos en 40 Ca, tienen un número
I0 I2 20
Z medio mayor que el de los músculos donde predominan
1
1
H, 126C y 32
16
O, y son más absorbentes que estos últimos. Por
el contrario, la probabilidad del efecto Compton varía poco
x con la naturaleza de los tejidos.
La figura 22-8 y la tabla 22-I muestran las variaciones de
y µ en función de la energía de los fotones X del haz incidente y
de la naturaleza del medio atravesado. Por orden decreciente
Figura 22-7. Contraste de la imagen radiante. de atenuación encontramos los huesos, después el agua, los
músculos y las vísceras que tienen µ muy próximos, la grasa y
los pulmones (poca atenuación pues contienen mucho aire).
Por otra parte, cuando la energía de los fotones aumenta,
decrece y los valores de µ para los diferentes tejidos se acer-
Por tanto, el contraste disminuye si aumenta el espesor can unos a otros, lo que lleva a un relativo «aplastamiento»
de los tejidos atravesados por los rayos X por delante y por del contraste.
detrás de la estructura estudiada. Los rayos X obtenidos con tensiones elevadas (100-
130 kV ) son pues muy penetrantes: se denominan «duros»
y son en particular utilizados para las radiografías óseas.
Origen del contraste en la imagen radiante Por oposición a éstos, los rayos «blandos» (25-50kV ) per-
miten una mejor discriminación de los tejidos blandos
Teniendo en cuenta la energía de los fotones X utilizada con un mejor contraste de la imagen radiante. En este
(entre 25 y 130 keV), las interacciones de los rayos X con los caso, la intensidad recibida por el detector es menor (lo
tejidos son principalmente de tipo fotoeléctrico y Compton. que requiere aumentar el duración de la prueba), y la do-
Las interacciones de tipo fotoeléctrico predominan para los sis absorbida por el paciente es más elevada, dado que µ
rayos X de baja energía (25-50kV); las de tipo Compton lo ha- aumenta cuando la tensión del tubo, y por tanto la energía
cen, para los rayos X de alta energía (110-130 kV). de los fotones disminuye.
Así, se puede definir tres gamas de energía de los haces
de rayos X, según sus características y aplicaciones:
– rayos X blandos entre 25 y 50 kV. Son poco penetran-
tes y sus interacciones tienen lugar principalmente por efec-
cm–1
5
µ
2
1
Tabla 22-I Variaciones del coeficiente de atenuación (µ)
0.5
en función de la energía de los fotones X del haz incidente
0.3 hueso
y de la naturaleza del medio atravesado
0.2 agua
0.1 grasa Hueso Músculos Grasa Pulmones
0.05 pulmón ρ (g.cm–3) 1.8 1 0.9 0.3
0.03
E media de los fotones Z medio 20 7 6.5 7
µ (cm–1) a 50 keV 0.9 0.3 0.25 0.07
50 100 250 keV
tensión del tubo µ (cm–1) a 70 keV 0.5 0.2 0.15 0.06
Figura 22-8. Variaciones de en función de la energía de los fotones µ (cm ) a 120 keV
–1
0.35 0.18 0.13 0.05
en diferentes tejidos biológicos.
296 IMÁGENES
rayos X
objeto
Radiografía estándar transparencia del medio atravesado por los rayos X. Estas
zonas son por lo tanto oscuras (casi negras) en la película;
Principio es, por ejemplo, el caso de los pulmones. Por el contrario, las
«opacidades» hacen referencia a la radiopacidad y aparecen
La radiografía estándar utiliza una fuente de rayos X del
claras, es decir en blanco en la película; es por ejemplo el
tipo estudiado más arriba: un detector constituido por una
caso de las calcificaciones.
película fotográfica acoplada a una pantalla intensificadora
La imagen obtenida es una proyección cónica cuyo pun-
y a dispositivos de eliminación de la radiación difundida. El
to más elevado es el foco del generador de rayos X sobre la
paciente está situado entre la fuente de rayos X y el detector, lo
placa fotográfica. Los elementos anatómicos se superponen,
más cerca posible de este último (figura 22-10). La fuente emi-
ligeramente agrandados y deformados. El tamaño en la ima-
te durante un corto tiempo de exposición un haz uniforme,
gen es más grande que en la realidad. Las estructuras estarán
prácticamente monoenergético, dirigido hacia el paciente.
menos agrandadas cuanto más próximas se encuentren de
Este haz se atenúa de manera diferente según la composición
la película en el momento de la toma del negativo, y cuanto
y el espesor de los diferentes tejidos atravesados. Las partes
mayor sea la distancia fuente-paciente. En efecto, si se deno-
del haz más atenuadas dan lugar a un menor ennegrecimien-
minan F y S respectivamente la distancia del foco y de la es-
to de la película. Después del revelado de la película, estas re-
tructura anatómica al detector (figura 22-11) esta estructura
giones corresponderán a las regiones más transparentes del
se agrandará con un factor de homotecia (o amplificación
negativo. Así, los huesos, que son muy radiopacos, aparecerán
vectorial):
en blanco en el negativo observado en el negatoscopio.
¡Ojo! Cuando en una radiografía se habla de «claridad» F
r= (22-9)
o «hiperclaridad», se hace implícitamente referencia a una F−S
rayos X
fuente
radiografía
negatoscopio
placa
fotográfica
2 pendiente γ
1.5
F 1 umbral
0.5
0 log(E)
objeto
S
imagen
película
A B C
Figura 22-11. Ampliación homotética de los objetos. Figura 22-12. Ejemplo de curva de sensibilidad del par
película-pantalla.
⎛B ⎞⎟ ⎛ E + ∆E ⎞⎟ ⎛ (I + ∆I) + T ⎞⎟ ⎛ I + ∆I ⎞⎟
DO = log10 ⎜⎜ 0 ⎟⎟ = log⎜⎜ ⎟ = log⎜⎜ ⎟⎟ = log⎜⎜ ⎟
⎜⎝ B ⎟⎠ ⎜⎝ E ⎟⎟⎠ ⎝⎜ I×T ⎟
⎠ ⎝⎜ I ⎟⎟⎠
por lo tanto: ⎛ I + ∆I ⎞⎟ 1 ⎛ I + ∆I ⎞⎟
∆x = log⎜⎜ ⎟= ln ⎜ ⎟=
⎜⎝ I ⎟⎟⎠ M ⎜⎜⎝ I ⎟⎟⎠
B = B0.10–DO (22-10)
1 ⎛ ∆I ⎞⎟ 1 ∆I 2
DO depende de la exposición E = I T, producto de la = ln ⎜⎜1 + ⎟⎟ CR
M ⎝ ⎜ ⎟
I ⎠ M I M
intensidad I por la duración de la prueba T, y de la sensibili-
dad del par película-pantalla, caracterizado a su vez por la (Recordemos que la constante M = 2.3 permite pasar de lo-
curva de sensibilidad, de la cual la figura 22-12 muestra un garitmos neperianos a logaritmos decimales).
298 IMÁGENES
∆x corresponde por tanto a un contraste dado de la ima- (hiperclaridad), sobreexpuesta, con un muy deficiente con-
M traste de la imagen luminosa.
gen radiante C R ∆x .
2 Para que el contraste de una determinada región de la
Del mismo modo, un intervalo dado ∆DO de ordenadas imagen radiante se pueda detectar, una vez revelada la pe-
corresponde a densidades ópticas DO y DO + ∆DO, por tanto lícula, es preciso que la energía total correspondiente a esta
a un contraste de la imagen luminosa: parte, y recibida por la película durante la duración de la
prueba, esté en la parte rectilínea, de pendiente γ, de la curva
(B + ∆B) − B 10−DO − 10−DO−∆DO 1 − 10−∆DO M
CL = = ∆D de sensibilidad.
(B + ∆B) − B 10−DO− ∆DO
+ 10−DO
2 2
La curva de sensibilidad en coordenadas semilogarítmicas
establece, por lo tanto, una correspondencia inmediata entre el Influencia de la tensión, de la intensidad
contraste de la imagen radiante en abscisas y el de la imagen lu- y del tiempo de la prueba
minosa en ordenadas (véase figura 22-12). Se tiene, en efecto:
Cuando la región anatómica a examinar es muy hetero-
CL Μ.∆DO/2 ∆DO
= = (22-9) génea, el importante contraste presente de por sí en la imagen
CR M.∆x/2 ∆x radiante entre las zonas radiopacas y radiotransparentes se
La relación de los contrastes de la imagen luminosa y de refuerza incluso más en la utilización de rayos X blandos. En
la radiante es por tanto igual a la pendiente media de la curva este caso, no es posible, por lo general, explorar correctamen-
de sensibilidad (en coordenadas semilogarítmicas). La figura te y simultáneamente las partes radiopacas y las radiotranspa-
22-13 muestra para un mismo contraste de 10% de la imagen rentes. Por ejemplo, en el caso de una radiografía de tórax, si se
radiante, el contraste obtenido en la imagen radiológica para utilizan tensiones bajas, y se regula la intensidad y la duración
tres valores de E en particular. de la prueba para tener una buena visualización del parénqui-
La curva de sensibilidad presenta pues tres zonas dife- ma pulmonar, los huesos aparecen uniformemente blancos
rentes: y enmascaran los tejidos subyacentes. De forma general, se
– la primera parte, de débil pendiente, corresponde a in- puede decir que la ganancia de contraste conduce a una re-
tensidades demasiado débiles de la imagen radiante y a una ducción de la gama de atenuaciones que se pueden estudiar.
película infraexpuesta. La imagen luminosa es casi unifor- En la práctica, se utiliza el bajo voltaje para estudiar
memente blanca (opacidad). Para valores muy pequeños de estructuras poco radiopacas y relativamente homogéneas
E, el ennegrecimiento de la película no es nulo y se percibe (p. ej., la mama). Para las estructuras radiopacas (radiografía
un «velo de fondo» debido al revelado de granos de plata en ósea, o prueba con agente de contraste positivo), la utilización
el momento de la fabricación de la película y de la opacidad de bajo voltaje conduciría a una irradiación excesiva para el
del soporte de poliéster. El velo de fondo tiene una DO de paciente, para conseguir exponer correctamente las partes
alrededor 0.12 (atenuación 24 pag. 100); de la película correspondiente a las estructuras radiopacas y
– la parte media de la curva es rectilínea. Su pendiente radiotransparentes. Permiten estudiar simultáneamente, por
se denomina como la γ de la película. La prolongación de esta ejemplo en el caso de la radiografía del tórax, algunos ele-
parte rectilínea corta el eje de las abscisas en un punto que co- mentos óseos, el parénquima pulmonar y el mediastino.
rresponde a un valor de exposición U denominado «umbral». El aumento de la intensidad y del haz de electrones o la
En esta parte rectilínea, se tiene, según (22-15), y en este caso: duración de la prueba T, con igual tensión de tubo U, provoca
un aumento proporcional de la intensidad I0 del haz inciden-
CL γ.CR (22-16)
te y, proporcionalmente, de los haces transmitidos que cons-
El contraste de la imagen luminosa es por tanto igual al tituyen la imagen radiante. Si, por ejemplo, el producto iT se
de la imagen radiante, multiplicado por la γ de la película. multiplica por k, dos zonas de la imagen radiante I e I + ∆I
Esta zona corresponde a una exposición óptima de la pelícu- tienen nuevas intensidades k.I y k.(I + ∆I). Esto no cambia el
la y a una amplificación del contraste de la imagen radiante contraste de la imagen radiante que se convierte en:
por parte del par película-pantalla;
– finalmente, la última parte de la curva, de nuevo con k(I + ∆I) − kI (I + ∆I) − I ∆I
R =
C *k = = CR (22-17)
pendiente más suave, corresponde a una película casi negra k(I + ∆I) + kI (I + ∆I) + I 2I
300 IMÁGENES
foco foco
puntual no puntual
objeto objeto
película película
imagen imagen
zonas de penumbra
Figura 22-15. Desenfoque geométrico ligado a las dimensiones del soporte del tubo de rayos X.
cámara de televisión
fuente de
rayos X
haz de
objeto fotones de
a examinar fluorescencia
serie
de detectores pantalla
haz
de electrones de vídeo
cibida por el paciente con respecto a un examen radioscópico digital, que hace un muestreo del valor analógico con una fre-
clásico de idéntica duración. El registro continuo de las imá- cuencia determinada para cada línea. El valor se transmite a
genes permite estudios cinéticos en diferido y ralentizados. El un ordenador, que compone y almacena la imagen digital.
estudio de la circulación coronaria inyectando, por medio de Este proceso es muy rápido y permite adquisiciones di-
un catéter, un agente de contraste positivo en cada coronaria, námicas en tiempo real. Las principales ventajas de esta di-
es una de las aplicaciones más extendidas de esta técnica. gitalización son las posibilidades de visualización y de trata-
miento ulterior de las imágenes, y en particular por técnicas
de sustracción.
La sustracción permite visualizar la vascularización de
Radiografía por sustracción digital una región anatómica eliminando las estructuras parásitas,
Las técnicas de la imagen digital que hemos descrito en en particular las óseas. Se realiza una primera adquisición,
el capítulo 20 son cada vez más utilizadas en radiología. Ac- denominada máscara, que se conserva como imagen digital
tualmente, se usan principalmente en la tomodensitometría, M(l,c). En todo punto, M(l,c) es proporcional a la intensidad
pero están ganando terreno en radiología clásica, particular- I0 del haz incidente, multiplicada por el factor de atenuación
mente en el área de la imagen vascular y torácica. de los rayos X denominado A(l,c) en dicho punto de la ima-
La digitalización en un escáner de la película radioló- gen (ligado a los huesos, a los tejidos blandos, etc.):
gica es técnicamente posible, pero suministra imágenes de M(l, c) = k I0 A(l,c) (22-18)
calidad mediocre. En la práctica, las radiografías digitales se
A continuación, después de la inyección de un agente de
obtienen directamente, sin utilizar una película como sopor-
contraste, se realiza una segunda adquisición, estrictamente
te, por digitalización de la señal de vídeo obtenido con una
superponible. La imagen correspondiente N(l,c) traduce la
cámara. Dicha cámara puede estar acoplada, bien a la pan-
atenuación en cada punto (l,c) del haz: a) por el factor A(l,c)
talla secundaria de un amplificador de luminancia (técnica
que traduce la presencia de los huesos y de los tejidos blan-
utilizada en angiografía digital), bien a una placa fosfolumi-
dos; b) por un segundo factor V(l,c) que traduce la presencia
niscente, que sirve de soporte intermediario entre la imagen
eventual del agente de contraste:
radiante y la imagen captada por la cámara. Otra técnica,
muy prometedora, utiliza la cámara de hilos de Charpak. N(l,c) = k I0 A(l,c) V(l,c) (22-19)
Dividiendo para cada píxel (l,c) la imagen N(l,c) por el
Angiografía digital valor M(l,c), se obtiene una nueva imagen, enteramente arti-
ficial, cuyo valor en cada punto es:
Se utiliza un amplificador de luminancia cuya pantalla N(l,c) K × I0 × A(l,c) × V(l,c)
secundaria fluorescente de fósforo es captado por una cáma- I(l,c) = = = V(l,c)
M(l,c) k × I0 × A(l,c)
ra de vídeo. El receptor de la cámara está constituida por una
sustancia (trisulfuro de antimonio u óxido de plomo) que se Esta imagen digital obtenida por división representa,
vuelve conductora cuando le llegan los fotones luminosos que por tanto, la imagen exclusiva del agente de contraste en los
emergen del amplificador de brillancia. El captador es barri- vasos sanguíneos. En realidad, se utiliza una representación
do, línea por línea, por un haz de electrones que genera una logarítmica de las imágenes, que transforma la división en
corriente más intensa cuanto más conductor se vuelve. Se ob- sustracción, de ahí el nombre de la técnica. La figura 22-17
tiene de esta forma, línea por línea, una señal analógica que ilustra este método en el caso de una angiografía cerebral en
será función creciente de la luminancia de la pantalla. Esta la cual las estructuras óseas dificultan considerablemente la
señal es digitalizada por medio de un convertidor analógico- visualización de los vasos en una radiografía estándar.
302 IMÁGENES
A B C
Figura 22-17. Angiografía cerebral con sustracción. A) Máscara. B) Angiografía sin procesar. C) Imagen tras sustracción de la máscara. (negativos
amablemente cedidos por el Doctor Zouaoui, servicio de Neurorradiología, hospital de La Salpêtrière, Paris).
Radiografía digital por sistema de placas con menos pasos intermedios que en las técnicas preceden-
tes y con una resolución espacial de 0.6 0.6 mm. La mayor
Esta técnica utiliza la digitalización de la señal de una eficacia del elemento detector y la disminución de la radia-
cámara de vídeo acoplada a una placa fosfoluminiscente de ción difundida, obtenida por el sistema colimador utilizado,
memoria. Esta placa se expone a la imagen radiante, en el permite a la vez reducir considerablemente la irradiación de
lugar de la pareja película-pantalla utilizada en la técnica es- los pacientes (un factor de 8 a 10 para la radiología bidimen-
tándar. La imagen radiante convierte en luminiscentes a los sional y 800 a 1000 en la tomodensitometría tridimensional)
cristales de la placa, con una intensidad que aumenta con la y mejorar el contraste. Un sistema de tratamiento informáti-
del haz de rayos X, durante un tiempo suficiente para anali- co de adquisiciones de frente y de perfil permite reconstruir
zar la placa punto por punto en un momento posterior. Este imágenes osteoarticulares de calidad excelente.
análisis se realiza con una cámara de vídeo o con un disposi-
tivo de barrido óptico por un haz láser. En este último caso, la
señal grabada para cada píxel de la imagen digital se codifica Contraste de la imagen luminosa
en 10 bits (puede por tanto tomar valores comprendidos en- Una de las ventajas decisivas de la radiología digital es la
tre 0 y 1023), lo que representa unos valores dinámicos 100 flexibilidad con la que se pueden obtener toda una gama de
veces superior a una película radiográfica clásica, con una imágenes visuales, a partir de la adquisición inicial en la que
resolución espacial de 10 píxeles/mm. Se obtiene, de esta se convierte la simple digitalización de la imagen radiante.
manera, una imagen digital de alta calidad, que, además, En este caso, no se es tributario, como en el de la radiografía
puede ser tratada posteriormente, por ejemplo para mejorar clásica, de la curva de sensibilidad del par película-pantalla.
su contraste o compensar digitalmente la existencia de infra- Los métodos de segmentación o de realce de contraste (véase
o sobreexposición. capítulo 26) permiten descargar la imagen en una pantalla
Cada vez más utilizados para las imágenes torácicas, es- de ordenador, modificarla en tiempo real, y luego reproducir
tos sistemas tienen, no obstante, una resolución insuficiente en película o en papel, una parte arbitraria de los valores de
para sustituir todas las indicaciones de las radiografías es- la imagen digital con un contraste máximo. En la práctica, ya
tándar, pero su calidad está mejorando rápidamente. no existe la «radiografía fallida» que obliga a rehacer la prue-
ba. En definitiva, el contraste de la imagen luminosa no está
limitado por los fenómenos de ruido y por la resolución de
Cámara de hilos de Charpak
intensidad del sistema de detección.
Este sistema, llamado «EOS» está poco difundido por el
momento, pero sus propiedades le auguran un gran futuro.
La radiografía se realiza línea por línea, el emisor móvil de Tomodensitometría
los rayos X pasando delante de una cámara proporcional
de múltiples hilos de Charpak que comporta una batería de La tomodensitometría (también llamada TDM, escáner
640 elementos de detección (véase capítulo 16, Cámara de de rayos X o «escanografía») ha revolucionado la imagen mé-
hilos). La señal se obtiene directamente de manera digital, dica y supuso la concesión del premio Nóbel de Medicina en
304 IMÁGENES
Translación (ρ) Translación (ρ)
E E
RX
RX
Translación
Translación
E : emisor de rayos X
: detector
E
E
Rotación
Rotación continua
RX
Con los escanógrafos actuales, la proyección según Ra- partir del coeficiente de atenuación lineal del agua, µagua. Un
don no es ortogonal sino en abanico, sobre unos 40º. El algo- tejido de coeficiente de atenuación lineal µ tendrá una opa-
ritmo de reconstrucción es más complicado que el método cidad UH:
de retroproyección filtrado digital que hemos descrito, pero
se funda en los mismos principios. Comprende una etapa µ − µ agua
C UH = 1000 × (22-20)
adicional de puesta a punto de los datos; utiliza una repre- µ agua
sentación en abanico y filtros adaptados a esta geometría.
Los filtros de reconstrucción se eligen en función de los Esta escala se mueve entre un valor 0 para el agua, –1000
objetivos del examen: para el aire (puesto que µaire 0) y un valor cercano a +1000
– los filtros de alta frecuencia de corte (sharp) tienen para el hueso (en realidad un poco más para los huesos muy
una menor resolución en densidad, una mayor sensibilidad densos). Los UH de algunos tejidos se muestran en la tabla
al ruido, pero una excelente resolución espacial (0.5 mm). Se 22-II.
utilizan para la exploración de los huesos, pulmones, etc.; Las imágenes se descargan en la pantalla de un orde-
– los filtros de frecuencia baja de corte (smooth) tienen nador, con numerosas posibilidades de tratamiento y de
un efecto de alisamiento espacial y una resolución espacial visualización (en particular segmentación, zoom en una
mediocre (2 mm), pero una muy alta resolución en intensi- parte de la imagen con interpolación, reconstrucción 3D,
dad. Se utilizan para la exploración del hígado, de los discos medida directa de las longitudes, de densidades, de áreas
intervertebrales, etc. y ángulos) y pueden ser impresas en una película fotográfi-
Unidades Hounsfield. Las unidades de atenuación lineales ca o transmitidas a otro sistema informático bajo formato
se expresan en unidades Hounsfield (UH), escala definida a digital.
A B
Figura 22-20. Ejemplo de segmentación. A) Ventana de discriminación entre tejidos blandos (0-50UH). B) Ventana ancha (0-1000 UH).
Negativos amablemente cedidos por el Doctor Zouaoui, servicio de Neurorradiología, hospital de La Salpêtrière, Paris).
306 IMÁGENES
Riesgo de los exámenes radiológicos Tabla 22-III Coeficiente EDAP para algunas áreas
anatómicas
La irradiación del paciente y del personal y las complica-
Región explorada EDAP
ciones de las sustancias de contraste constituyen los princi-
pales riesgos de los exámenes radiológicos. Tórax de frente (130 kV) 0.33
Tórax de perfil (130 kV) 0.15
308 IMÁGENES
Imagen por resonancia
magnética 23
La imagen por resonancia magnética (IRM) se basa en to tiene una orientación y un módulo que expresa la «fuer-
el registro de las velocidades de retorno al equilibrio de pro- za» del imán. Un imán crea, en el espacio circundante, un
tones expuestos a un campo magnético, tras ser sometidos campo magnético, denominado B, magnitud vectorial cuya
a un estímulo «desestabilizador». En veinte años, los pro- intensidad, denominada inducción magnética, se expresa
gresos alcanzados por esta técnica no ionizante han sido en teslas (T). El gauss, la unidad utilizada antiguamente, es
espectaculares: la mejora considerable de la resolución es- igual a 10–4 T.
pacial y de la intensidad, la rápida secuencia de adquisición, Del mismo modo que un dipolo situado en un campo
la introducción de productos de contraste, la posibilidad de eléctrico está sometido a un par y se orienta en el sentido
obtención de imagen funcional y el desarrollo de protocolos del campo, un momento magnético M situado en el campo
intervencionistas, han hecho de la IRM un método costoso magnético B creado por otro imán está sometido a un par
pero rico en posibilidades reales y potenciales, ventajas no =
Γ M
B y se orienta en el mismo sentido del campo. Un
disponibles en ninguna de las modalidades concurrentes de ejemplo conocido desde hace siglos es la brújula, imán que
la imagen médica. se orienta en el campo magnético terrestre (cuya componen-
te horizontal está dirigida de Sur a Norte y su intensidad me-
dia es de 0.5 10–4 T).
Fenómeno de resonancia magnética Las propiedades magnéticas de la materia están liga-
das a la presencia de cargas eléctricas en movimiento. De
El fenómeno de resonancia magnética (RM) descansa
un modo general, toda carga en movimiento crea una in-
sobre la existencia del «espín nuclear», que confiere a algunos
ducción magnética y un hilo eléctrico recorrido por una co-
núcleos un momento magnético responsable de su interac-
rriente desvía la brújula. Una bobina recorrida por una co-
ción con un campo magnético externo. Descubierto en 1946
rriente crea un campo magnético orientado según el eje de la
de manera independiente por Félix Bloch y Edward Purcell,
bobina y cuya intensidad es proporcional a la de la corriente.
este fenómeno dio lugar, en primer lugar, a la espectroscopía
La IRM utiliza inducciones magnéticas muy elevadas
por RM, método de análisis químico potente y no destructi-
(0.5 a 4 T) creadas por grandes bobinas resistivas o super-
vo, y en 1972 a las primeras imágenes 2D obtenidas por Paul
conductoras.
Lauterbur a partir de una señal RM. La denominación de RM
proviene del hecho de que el núcleo de determinados áto-
mos posee las características de un imán, capaz de oscilar
entre dos estados posibles cuando es situado en el seno de Propiedades magnéticas de núcleo
un campo magnético y de resonar cuando está sometido a
El núcleo de los átomos está animado por un movimien-
una excitación de frecuencia adecuada.
to de rotación sobre sí mismo, llamado movimiento de es-
pín. Este movimiento de una carga eléctrica positiva confiere
Campo y momento magnéticos al núcleo, bajo determinadas condiciones de asimetría, un
Un imán es un sistema material dotado de un momento momento magnético llamado momento magnético nuclear
magnético, denominado
M, magnitud vectorial que por tan- y denominado µ.
309
Entre los núcleos cuantitativamente importantes en el B0
organismo, tan sólo el 1H posee un momento magnético nu-
momento
clear. Otros isótopos, los que poseen número impar de pro- paralelo
tones o de neutrones (13C, 19F, 23Na, 31P) tienen igualmente
un momento magnético, no así el 12C o el 16O. Por lo tanto,
la imagen RM actual utiliza esencialmente las propiedades �
del núcleo de hidrógeno. Describiremos el fenómeno de re-
sonancia magnética para este núcleo, siendo inmediata su
generalización para el resto.
Los momentos magnéticos nucleares de los núcleos
de hidrógeno contenidos en un volumen V de materia se �
suman vectorialmente para dar un momento magnético
macroscópico M
= ∑ µ.
Los momentos magnéticos ma-
V momento
croscópicos que describiremos tienen un módulo propor- antiparalelo
cional al número ρ de núcleos de hidrógeno por unidad de
volumen (expresado en átomos-gramos por m3).
Figura 23-2. Orientación de los momentos magnéticos elementales
En ausencia de un campo magnético externo, los mo-
en presencia de un campo magnético externo
B0 .
mentos elementales µ tienen orientaciones aleatorias que se
anulan en promedio, por lo que M =
0(figura 23-1). Z
Por el contrario, en presencia de un campo magnético
externo B0, se produce un fenómeno doble:
– los momentos magnéticos nucleares individuales µ se
orientan y forman un ángulo constante θ con B0 (figura 23-2).
Algunos momentos, mostrados como ⇑, se orientan en para-
lelo con respecto a
B0, esto es, de forma que (µ. B0 0), y otros,
mostrados como ⇓, se orientan en antiparalelo (µ.
B0 0);
– al igual que una peonza girando sobre sí misma, alre-
dedor de un eje vertical, los momentos magnéticos nuclea- B0
momento
res µ están animados por un movimiento de rotación, dicho paralelo
de precesión, alrededor del eje definido por B0 (figura 23-3).
Figura 23-3. Movimiento de precesión de los momentos magnéticos
Dicho movimiento se produce a la misma frecuencia y en el
elementales .
mismo sentido (positivo para la orientación definida por B0)
para los momentos paralelos o los antiparalelos. De ello re-
sulta que, asimismo, el momento resultante M gira alrededor La constante γ, característica del núcleo, se denomina
de
B0, a la misma frecuencia. En ambas orientaciones (⇑ y ⇓), relación o también constante giromagnética. La frecuencia
los momentos individuales no están en fase y se disponen de rotación de precesión, denominada frecuencia de Lar-
aleatoriamente alrededor de B0 (figura 23-4). mor, viene dada por:
La velocidad angular ω 0 del movimiento de rotación de
ω0 γ
precesión viene dada por la relación de Larmor: ω 0 = γ
B0. ν0 = = B0 (23-1)
2π 2π
ω0 γ
Los νvalores
0 = =
de ,Bexpresados
0 en revoluciones por se-
2π 2π
gundo y Tesla (o en hercios por Tesla) se indican en la tabla
23-I para algunos núcleos importantes.
Las posiciones paralelas y antiparalelas corresponden a
niveles de energía diferente. Un momento magnético parale-
γ hB0
lo tiene una energía individual negativa e ⇑ = − (sien-
4π
do h la constante de Planck = 6.62 10 J.s) y un mo-
–34
310 IMÁGENES
que denominaremos
individuales µ, M0 , no es nula, sino que
ML no nula, paralela a
tiene una componente longitudinal B0 y
con el mismo sentido. El módulo ML es igual a la diferencia en-
tre los módulos de los momentos magnéticos elementales ⇑ y ⇓;
ML
por consiguiente, después de proyectar en la dirección de
B0:
⎛ N( ) ⎞
ML = (N(⇑) − N(⇓))µ.cos θ = N(⇓)⎜⎜ ⇑ − 1⎟⎟⎟µ.cos θ =
⎜⎝ N(⇓) ⎟⎠
∆e
90°
MT Bobina
de detección 270°
MT
B0 t
Señal sinusoidal
180°
Fenómenos de relajación
Figura 23-6. Recorrido en el espacio de
ML en el curso del fenómeno
Relajación transversal
de resonancia magnética. Supongamos un pulso 90° interrumpido en el momen-
to que M se encuentra en el plano perpendicular a B . El
0
En el espacio, el extremo del vector
M describe una curva
momento magnético global M, que ha sido apartado de su
compleja en la superficie de una esfera y en forma de hélice, posición de equilibrio, va a retornar espontáneamente a su
debido a la rotación de la frecuencia de Larmor (figura 23-6). posición inicial: a este retorno se le llama relajación. En el
Si la estimulación «resonante» por el flujo de fotones de espacio, el extremo de M describe un movimiento comple-
frecuencia ν0 se mantiene, la inversión de los efectivos de las jo que lo acerca progresivamente a la dirección B0, girando
poblaciones ⇑ y ⇓ continúa hasta la anulación de la compo- a la frecuencia ν0 (figura 23-8). El módulo de la proyección
nente longitudinal
ML , siendo máxima la componente trans- transversal M (que gira a la frecuencia ν ) disminuye progre-
T 0
MT en esas condiciones. Más adelante,
versal ML se vuelve ne- sivamente de su valor máximo M0 hasta 0 (figura 23-9) según
MT disminuye hasta que
gativo y M sea el opuesto a su valor una ley exponencial: MT(t) = M0 .e–t/T2. Al parámetro T2 se le
inicial. Se observa a continuación un fenómeno simétrico, denomina tiempo de relajación transversal o espín-espín.
con incremento de M que pasa por valor nulo, para luego re-
T
Es el tiempo al cabo del cual la amplitud de MT se reduce por
tornar progresivamente al valor inicial positivo
M0 , mientras un factor 1/e = 0.37.
que
MT (en giro) pasa por un máximo, y después se anula. En En una bobina de eje perpendicular a B0 se genera
este momento, un ciclo idéntico puede empezar de nuevo. una corriente inducida de frecuencia ν0 y de intensidad I(t)
El ciclo descrito por
M en un sistema de referencia que proporcional a MT(t), por tanto progresivamente decre-
girara alrededor de B con la frecuencia de Larmor, puesto ciente (figura 23-10). Esta intensidad I(t) es el producto de
0
que en ese sistema, la rotación alrededor de B0 deja de ser cos(2πν0t) que traduce la rotación de MT a frecuencia ν0, por
312 IMÁGENES
��
M
��
M T
I(t)
I0
Imáx(t)
0
t
Imín(t)
Imáx(t)
I0/e
0 T2
I0/e t
Imín(t)
314 IMÁGENES
I(t) Imáx(t) ⇒ T2
I0
I*máx(t) ⇒ T2
Corriente inducida
Pulsos
90° t
Figura 23-12. Secuencia de saturación-recuperación.
I(t)
Imáx(t) ⇒T2
I0
I*máx(t) ⇒ T2*
Corriente inducida
TE/2 TE
14454244543
eco de espín
Pulsos
90° 180° t
Figura 23-13. Secuencia de eco de espín.
pero no de las diferencias de sus velocidades. Se obtiene una trada en el momento de los ecos sucesivos depende enton-
segunda acumulación de los corredores en la línea de llegada ces estrictamente de T1 (ponderación de T1).
haciéndoles que den media vuelta 1 hora y 30 minutos des- Estas propiedades se resumen en la tabla 23-II.
pués del inicio de la carrera (se supone que la velocidad de
cada corredor es constante). Secuencia de inversión-recuperación
En la práctica, la medida se efectúa en el momento del (inversion-recovery, IR)
eco (al tiempo TE) y las secuencias de eco de espín se repiten Esta secuencia tiene como objetivo obtener una señal
después de un tiempo TR (tiempo de repetición) y se suman. ligada a T1. Utiliza un primer pulso 180°, seguido de un pe-
Pueden considerarse dos casos: ríodo durante el cual la componente ML aumenta según la
– si el tiempo TR es largo frente a los tiempos T1 y T2, se función:
trata simplemente de la repetición del mismo experimento,
con el resultado explicado anteriormente. En ese caso, si TE ML(t) = M0(1 – 2e–t/T1) (23-7)
es corto frente a T2, la señal medida al cabo del tiempo TE no
disminuye excesivamente y refleja principalmente los valo- Tabla 23-II Secuencia de eco de espín
res de ρ (ponderación de ρ). Si por el contrario TE es largo, la
TE corto ( 30 ms) TE largo ( 60 ms)
señal refleja los valores de ρ, pero sobre todo los de T2 (pon-
deración en T2); TR corto ( 500ms) T1 T1, T2
– si el tiempo TR, así como el tiempo T2, son breves, se
TR largo ( 2 000ms) Densidad
repite la secuencia, mientras el sistema no haya regresado a
de protones T2
su equilibrio inicial. Se puede demostrar que la señal regis-
ML(t)
τ
t
I(t)
Corriente inducida
τ t
Pulsos
180° 90° t
Figura 23-14. Secuencia de inversión-recuperación.
316 IMÁGENES
Sustancia Sustancia
blanca gris
ML(t)
A B C t
Figura 23-15. Contraste en IRM. (La sustancia gris cerebral se muestra como una corona periférica rodeando la materia blanca central.)
318 IMÁGENES
B0
ν2
ν1
z
z1 z2
corte
seleccionado
.
Figura 23-17. Selección del corte estudiado por el gradiente de campo G z
se persiste durante la fase de recepción, lo que significa que rentes gradientes (selección del corte por
Gz ; codificación de
protones con idéntica abscisa (por tanto sometidos a idénti- la abscisa por Gx ; codificación de la ordenada por una serie
co campo Gx ) van a quedar desfasados con un mismo ángulo de
Gy ) y la evolución de la corriente inducida, durante una
que sólo depende de sus ordenadas respectivas. Las señales adquisición en «eco de espín» clásico. La escala del tiempo
emitidas por dos elementos de misma abscisa y ordenadas de las oscilaciones se ha ampliado para hacerlas visibles.
respectivas y1 e y2 aparecerán desfasadas con ángulos res-
pectivos θ1 = Θ(y1) y θ2 = Θ(y2) (se puede demostrar que la Imágenes IRM no transversales
función Θ es del tipo Θ(y) = γ B + (jy) ). Las corrientes in-
2
0
2
El principio de obtención de las imágenes IRM permite,
ducidas correspondientes I1(t) e I2(t) presentan desfases q1 y utilizando gradientes de campo magnético en una dirección
q2 que «codifican» por lo tanto el valor de la ordenada en el definida, obtener directamente una imagen tomográfica en
origen de la señal (figura 23-19). cualquier orientación (transversal, sagital, longitudinal u
La extracción de información de codificación de y, con- oblicua). En efecto, las bobinas de gradientes de campo
Gz ,
tenida en la señal en forma de desfases, necesita repetir la Gy y
Gx permiten obtener, a través de una suma vectorial,
misma secuencia de adquisición para valores crecientes del gradientes en cualquier dirección. Se puede pues generar un
módulo de Gy . Para reconstituir una imagen que contiene L primer gradiente de selección de corte, cualquiera que sea la
líneas (en la dirección y), se debe proceder a L adquisiciones dirección elegida, y a continuación un segundo gradiente de
para valores crecientes de j: j1, j2 … jL. Se obtiene de esta mane- codificación en fase, perpendicular al primero, y finalmente
ra una serie de L señales. La transformada de Fourier de cada un tercer gradiente de codificación por la frecuencia perpen-
señal lleva a separar las diferentes abscisas. Para cada absci- dicular a los dos primeros. Las secuencias de adquisición y
sa, se efectúa a continuación una transformada de Fourier las técnicas digitales de reconstrucción de imágenes no se
de los L valores obtenidos para las diferentes señales, lo que modifican. La figura 23-21 muestra un ejemplo de imagen
permite obtener los valores de las diferentes ordenadas. sagital de una cabeza obtenida en IRM.
Fourier (I(t))
Pico ligado a x1
Pico ligado a x2
0
B p.x 1 p.x 2
Figura 23-18. Codificación de la abscisa x ligada a la frecuencia de la señal: señal recogida. A) Señal compleja I(t) = I1(t) + I2(t).
B) Transformada de Fourier de la señal I(t).
I1(t) I2(t)
Desfase θ1 Desfase θ2
t t
Figura 23-19. Codificación de la ordenada ligada a la fase de la señal. (La figura ha sido realizada con θ1 π/4 y θ2 3π/4).
cias idénticas promediadas para aumentar la relación señal/ neadas, conducen a tiempos prolongados de examen, con
ruido o de las L secuencias para valores diferentes de
Gy . Si TR el doble inconveniente de reducir el número de exploracio-
es el tiempo de repetición (entre dos ecos o secuencias), la du- nes que se pueden practicar en la instalación y aumentar los
ración de la adquisición para cada corte es D = TR N L. efectos que el movimiento involuntario del paciente puedan
Para obtener una imagen 256 256, L = 256. Para las se- tener sobre la calidad de la imagen.
cuencias clásicas, TR está comprendido entre 300 y 2500 ms, Numerosas técnicas tienen por objetivo recortar el nú-
siendo N un número entre 1 y 4. El tiempo de adquisición mero de adquisiciones, a costa, frecuentemente, de un ligero
de cada corte está por tanto comprendido entre aproxima- deterioro de la calidad de imágenes. Se basan en lo siguiente:
damente 80 y 2500 s. Un examen que comporte un número – limitación del número de secuencias idénticas, lo que
generalmente importante de cortes (15 a 20), incluso si las conduce a reducir la relación señal/ruido o el número de lí-
adquisiciones de varios cortes pueden ser en parte simulta- neas (por tanto, la resolución espacial);
320 IMÁGENES
I(t)
Medida Medida
Corriente
t inducida
TR TR
TE/2 TE/2
Pulsos
90° 180° 90° 180° 90° t
��
Gradiente G z
t
��
Gradiente G x
t
��
Gradiente G y
t
El campo magnético principal B0 se obtiene con varios so- Resolución espacial y de intensidad
lenoides coaxiales, dispuestas según el método de Helmholtz
(figura 23-22). Estos solenoides son atravesados por una co- Los equipos de IRM permiten obtener rutinariamente
rriente de alta intensidad. Para limitar el consumo de electri- imágenes 256 256, 512 512 ó 1024 1024. Los cortes
cidad, se utilizan solenoides superconductores, mantenidos tienen unos milímetros de espesor. El campo explorado es
a muy baja temperatura ( 4 °K) en helio y nitrógeno líqui- variable, en función de la antena utilizada. Puede tratarse de
dos, que deben ser regularmente repuestos, debido a las pér- una sección transversal completa del cuerpo o de un cam-
didas por evaporación. Se utilizan también, aunque menos po reducido que explore con una gran resolución espacial
frecuentemente, imanes de resistencia (es decir, no super- una pequeña zona del organismo. Por ejemplo, una imagen
conductores) o imanes permanentes. 256 256 obtenida en un campo de 8 cm 8 cm tiene una
Los campos obtenidos tienen, según los equipos, inten- resolución espacial de 0.3 mm/píxel. Antenas especiales per-
sidades de 0.5 a 4 T, siendo los valores más corrientes 1 y 1.5 T. miten una imagen transcutánea de muy alta resolución.
En el centro del campo útil, la homogeneidad del campo al- La resolución en intensidad aumenta con la intensidad
canza ± 0.05 ppm (partes por millón) en un volumen esférico del campo magnético estático B0 y, para una misma secuen-
de 10 cm de diámetro y ± 5 ppm en uno de 50 cm de diáme- cia de adquisición, con el número de secuencias sumadas
tro. (en la práctica, de 1 a 4). La posibilidad de discriminar dos
tejidos descansa principalmente sobre las diferencias de los
valores T1 y T2, sobre la secuencia utilizada y sobre el instan-
Generación de los gradientes te en el que se ha realizado la medida:
del campo magnético – para obtener una imagen ponderada en T1, se utiliza
por ejemplo, bien una secuencia de inversión-recuperación,
Los gradientes de campo se obtienen con tres tipos de
aunque el tiempo de adquisición sea muy largo, bien una se-
dispositivos accesorios, cuya configuración geométrica per-
cuencia de eco de espín, con valores cortos de TE ( 30 ms)
Gz,
mite generar Gy y
Gx. Se trata de solenoides clásicos, cuya
y de TR ( 500 ms);
alimentación eléctrica permite un encendido muy rápido
– para obtener una imagen ponderada en T2, se utiliza
(estabilización en menos de un milisegundo). Los gradientes
por ejemplo una secuencia en eco de espín con un TE largo
obtenidos son del orden de 10 a 25 mT/m.
( 60 ms) y un TR muy largo ( 2000 ms).
��
B0
322 IMÁGENES
vascular y a la cuantificación del flujo sanguíneo. Se han menta localmente la intensidad de B0), mientras que la oxi-
utilizado varias técnicas, todas ellas fundadas sobre el des- hemoglobina (HbO2) no tiene este efecto. Por comparación
plazamiento de los protones del fluido durante la secuencia de las imágenes obtenidas antes y después de la activación,
excitación-adquisición. es posible identificar áreas corticales activadas. Reservada
Los métodos por «tiempo de vuelo» (o en inglés time- inicialmente para máquinas que disponían de un campo
of-flight) utilizan, por ejemplo, el hecho de que si entre la magnético estático muy intenso (4 T), en el momento actual
fase de excitación y la de adquisición los protones excitados únicamente se puede realizar con equipos 1.5 T. Un ejemplo
del fluido se desplazan y son reemplazados por protones no de IRM funcional se presenta en la figura 23-23.
excitados, la señal será más o menos intensa (según la se- La IRM funcional abre posibilidades considerables en el
cuencia utilizada) para el fluido en movimiento que para los campo de las ciencias cognitivas, de la neurofisiología y de la
tejidos circundantes inmóviles. exploración funcional del cerebro.
Los métodos por contraste de fase utilizan el desfase de
los espines producido por gradientes accesorios. Este des-
fase varía según la velocidad de los protones, pudiéndose
identificar la velocidad de desplazamiento de los elementos Espectrometría RM in vivo
del fluido. Se puede de esta manera cartografiar las velocida-
des de desplazamiento de un fluido en un vaso sanguíneo, Muy utilizada in vitro, la espectrometría RM se basa en
estudiar el régimen de flujo (laminar o turbulento) y cuanti- el hecho de que, en una secuencia de RM, la frecuencia de
ficar el caudal de los flujos laminares. Larmor de los diferentes protones depende de su entorno
Se pueden obtener, igualmente, imágenes de la ventila- químico pues éste modula localmente el valor del campo
ción pulmonar en tiempo real, gracias a una mezcla gaseo- magnético estático B0. Una muestra que contiene un lípido y
sa que contiene Helio-3, previamente hiperpolarizado por agua, por ejemplo, da lugar a una señal compleja resultante
bombeo óptico, y que da lugar a una señal IRM observable. de la superposición de la frecuencia de los H del agua y la de
Se trata todavía de técnicas en fase de investigación. los del lípido. Dichas frecuencias pueden ser separadas por
una transformación de Fourier.
El objetivo de la espectrometría RM es obtener informa-
ciones análogas in vivo. Dado que los desajustes de frecuen-
IRM funcional cia de Larmor según el entorno químico son muy pequeños
(unas cuantas partes por millón o ppm), es necesario obte-
La IRM funcional permite visualizar la activación de ner un campo B0 de elevada homogeneidad y una relación
áreas cerebrales, con una precisión topográfica del orden del señal/ruido excelente (por tanto un campo estático B0 inten-
milímetro, en el curso de tareas físicas (p. ej., mover los de- so, con un mínimo de 1.5 T) para poder evidenciarlos. Las
dos), o intelectuales (p. ej., contar) repetidas en un protocolo principales aplicaciones se refieren a la espectroscopía RM
de activación denominado por los psicólogos «paradigma». del fósforo, que permite el análisis del metabolismo energé-
Su principio descansa sobre la modificación del campo mag- tico del músculo. En efecto, es posible diferenciar y cuantifi-
nético local de los tejidos según la concentración local de car las formas químicas en las que se encuentran los núcleos
oxígeno, que aumenta con el riego sanguíneo de las zonas de fósforo en cada vóxel muscular: fosfato inorgánico, fos-
activadas. Esta modificación, llamada BOLD (blood oxygen focreatina, ATP. Un esfuerzo muscular se acompaña de una
level dependence) se origina por el hecho de que la hemo- bajada de los picos correspondientes a la fosfocreatina y al
globina no oxigenada (Hb) es paramagnética (por tanto au- ATP, con aumento de los picos del fosfato inorgánico.
D I D I
Figura 23-23. IRM funcional. Activación cerebral en el momento de la realización de movimientos voluntarios de los dedos de la mano izquierda
(parte derecha de la figura) o de la mano derecha (parte izquierda) para un diestro. Las áreas activadas se representan en blanco.
D: derecha; I: izquierda.
324 IMÁGENES
Imagen por ultrasonidos
24
Las ondas sonoras de alta frecuencia constituyen la base describen en un apartado del capítulo 11 (véase Señal física
de una de las modalidades de imagen médica más utilizada de la audición); recordemos algunas de las relaciones más
debido a su excelente poder de discriminación entre dos te- comunes entre estas magnitudes:
jidos diferentes y a su prácticamente total inocuidad. Estas – frecuencia del sonido (Hz): ν;
características han hecho de este tipo de prueba un instru- – velocidad del sonido (m.s–1): C;
mento irremplazable en obstetricia. Por su parte, el efecto – densidad del medio (Kg.m–3): ρ;
Doppler permite cuantificar la velocidad de desplazamiento – longitud de onda (m): λ = C/ν
de las células sanguíneas, por lo que ofrece un medio de in- ⎛ ⎛ x ⎞⎞
– elongación (m): u(x,t) = a cos⎜⎜2πν⎜⎜t − ⎟⎟⎟⎟⎟⎟ (24-1)
vestigar la patología cardiovascular de manera extremada- ⎜⎝ ⎜⎝ C ⎟⎠⎟⎠
mente eficiente. – presión acústica (Pa):
⎛ ⎛ x ⎞ π⎞
P(x,t) = p cos ⎜⎜2πν⎜⎜t − ⎟⎟⎟ + ⎟⎟⎟ (24-2)
⎜⎝ ⎜⎝ C ⎟⎠ 2 ⎟⎠
Propiedades físicas de los
ultrasonidos – impedancia acústica (kg.m–2.s–1 o rayls):
p
Los ultrasonidos (US) son ondas acústicas de frecuencia Z= = ρC (24-3)
2πνa
comprendida entre 20 kHz y 200 MHz. Las frecuencias inferio-
p2
res (de hasta 16 Hz) corresponden a los sonidos, audibles por – intensidad acústica (W.m–2): I = (24-4)
las personas jóvenes, y las más elevadas a los hipersonidos. 2Z
El interés de los ultrasonidos en imagen médica reside
en su poder de propagación en los tejidos y de reflexión en Reflexión y refracción de los ultrasonidos
las superficies de contacto entre dos medios de característi-
cas diferentes en cuanto a la propagación de los mismos. Cuando una onda US atraviesa la superficie de contac-
to (o «interfase») que separa dos medios materiales M1 y M2,
que transmiten el sonido a velocidades C1 y C2, e impedancias
Nociones fundamentales y nomenclatura acústicas Z1 y Z2, respectivamente (figura 24-1), se observa
que una parte de la energía incidente se refleja por M1 y una
En el marco de aplicación de los US en imagen, ninguna
parte se transmite en el medio M2. Si denominamos i, r y t a
de sus propiedades físicas los distingue de las ondas sono-
los ángulos que forman con la perpendicular a la interfase las
ras. Se trata pues de ondas elásticas longitudinales que sólo
direcciones de propagación de las ondas incidente, reflejada y
pueden propagarse en un medio elástico, sólido, líquido o
transmitida, respectivamente, se pueden demostrar las rela-
gaseoso. El estudio de las principales características de las
ciones de Snell, análogas a las leyes de Descartes en óptica:
ondas sonoras y la definición de los términos fundamentales
(onda acústica, onda de presión acústica, velocidad, longi- sen i sen t
i=r y = (24-5)
tud de onda, impedancia acústica, intensidad acústica) se C1 C2
Ii Ii
Ii
Ir It It
326 IMÁGENES
Tabla 24-II Coeficiente según la frecuencia de los ultra- + – + –
sonidos
+ + –
Interfase Frecuencia US (MHz) (dB.cm–1)
Tejidos blandos 1 1 – –
+ –
Tejidos blandos 3.5 3.5 +
+ + –
–
Pulmones 1 40
Pulmones 2 80 V=0
– +
Huesos 2,5 35
Figura 24-3. Efecto piezoeléctrico.
la longitud de onda. Este fenómeno, análogo a la dispersión efecto piezoeléctrico está ligado a la reorientación de las es-
(scattering) Rayleigh de las ondas electromagnéticas, se lla- tructuras dipolo bajo el efecto de un campo eléctrico creado
ma dispersión de la onda US. Este fenómeno atenúa la in- por la ddp aplicada (figura 24-3). Si se aplica una ddp varia-
tensidad de la onda US, pero permite obtener una señal de ble a lo largo del tiempo, las contracciones y dilataciones del
reflexión sobre las estructuras de tamaño muy pequeño. disco siguen las variaciones de la ddp.
Globalmente, la atenuación en función de la profundi- Si la ddp aplicada es sinusoidal, de amplitud fija V0 y de
dad de los tejidos atravesados sigue una ley exponencial que frecuencia variable ν, los desplazamientos de la cara anterior
se expresa en dB (decibelios) de atenuación. Si el haz de US del disco se producirán con la misma frecuencia ν, con una
tiene una intensidad acústica I0 a la entrada del material e Ix a amplitud a = V0.A(ν), proporcional a V0 (en tanto no sea de-
la profundidad x, su atenuación será, por definición: masiado elevada) y a una función de la frecuencia A(ν), de-
pendiente no sólo de la naturaleza y de la geometría del dis-
⎛I ⎞⎟
A x (db) = 10 log⎜⎜ 0 ⎟⎟ (24-10) co, sino también de la resonancia entre el disco piezoeléctri-
⎜⎝ I ⎟⎠ co y la frecuencia impuesta por la ddp. Los desplazamientos
x
Filamento
Fils z
de
d'excitation ρ
excitación
et recueil
y recogida Emission
Emisión US π
Lente acústica
Placa adaptadora
Bloc postérieur de impedancia Figura 24-7. Enfoque del haz emitido por una sonda
Bloque posterior rL
absorbant Elemento monotransductora. La distancia focal viene dada por f
absorbente piezoeléctrico 1 − C/C L
en la que rL es el radio de curvatura de la cara anterior de la lentilla, CL y
Figura 24-5. Sonda monotransductora. C las velocidades de los US en la lente acústica y en el medio estudiado.
328 IMÁGENES
Excitación
Transductor
aislado
Haz fuertemente
divergente
onda esférica
Sector
++ explorado
Haz débilmente Transductores
divergente agrupados
onda plana
eje de rotación
Sondas multitransductoras
Estas sondas, actualmente las más extendidas, están
constituidas por una batería de emisores, excitables de ma- Figura 24-10. Enfoque electrónico de una sonda multitransductora
lineal.
nera independiente. Los transductores se disponen en forma
de anillo o en línea recta (figura 24-9). El número de trans-
ductores es variable (influyendo, evidentemente, en el precio Características de la detección
de la sonda): por ejemplo 128 transductores en disposición
lineal. ecográfica
Cada transductor produce un haz rápidamente di-
vergente, pero, a pesar de que el ángulo de divergencia La ecografía se basa en la emisión breve de una onda US
θ = Arc sen (1.2λ0/D) es grande para un emisor pequeño, el y en el análisis de las ondas reflejadas por un medio material.
acoplamiento de varios transductores permite obtener un El transductor se utiliza para la emisión incidente y la recep-
haz de menor divergencia. Además, retrasando sucesiva- ción de la señal reflejada. Los US reflejados por una superfi-
mente las emisiones de los transductores acoplados se pue- cie de contacto retornan al transductor con un retraso que
den crear interferencias que lleven a cabo un enfoque del será mayor cuanto mayor sea la profundidad a la que se en-
haz a profundidades y direcciones variables (figura 24-10). cuentra la interfase. El análisis temporal de la señal reflejada
Como cada elemento es de pequeño tamaño, se trabaja en la permite por tanto conocer la localización de la interfase en
zona distal de Fraunhofer, obteniéndose así una gran homo- el eje del haz. La calidad de esta detección se cuantifica por
geneidad y un enfoque muy fino del haz de US. su resolución axial (en profundidad) y su resolución lateral,
Del mismo modo, se aplica idéntico retraso a los ecos re- su profundidad de penetración y su profundidad de campo,
cogidos por los diferentes transductores, antes de realizar el y finalmente, por su dinámica. Las otras características im-
sumatorio, lo que permite un enfoque de la recepción. Este portantes son la frecuencia de recurrencia y la ganancia de
enfoque en la recepción se puede modificar durante toda la profundidad.
fase de recepción para que, en el instante ∆t después del pul-
so, el receptor sea enfocado a una distancia C∆t/2, preparado Resolución
para recibir los ecos provenientes de esta zona. Se obtiene de
esta manera un enfoque dinámico de recepción. Resolución axial (en profundidad)
La resolución axial, o en profundidad, es el parámetro
que caracteriza la posibilidad de discriminar señales ecográ-
transductor ficas provenientes de estructuras de profundidades próxi-
mas z y z + ∆z. Si la velocidad del medio es C, el segundo
eco alcanzará el transductor con un retraso correspondiente
aislante
a un recorrido de ida y vuelta en la distancia ∆x con respecto
al primero: ∆t = 2∆x/C. Los dos ecos sólo se podrán distin-
guir si la duración del pulso incidente es inferior a ∆t (figura
24-11). La duración de la excitación eléctrica del transductor
es despreciable, pero éste continúa vibrando y la duración
Disposición anular Disposición lineal
real de la emisión de US está ligada a la rapidez con la que
Figura 24-9. Sondas multitransductoras. las vibraciones del transductor se amortiguan. Hemos visto
f 2λ 0 f 2C
π3 =3 2 (24-14)
ida Dz D 2
D ν0
gel vuelta
Debido a que el examen ecográfico se lleva a cabo so-
�t
bre estructuras de profundidad variable, no se puede elegir
Q elevado especialmente una profundidad dada, y la profundidad de
campo debe tener un valor importante. El enfoque electróni-
t co permite realizar sucesivamente varias exploraciones en la
misma dirección, con distancias focales distintas, y superpo-
Emisión Recepción
ner en la misma imagen los resultados obtenidos.
�t Profundidad de penetración
Q pequeño
La penetración de los US es el parámetro que caracteri-
t
za empíricamente el espesor para el que la intensidad acús-
Emisión Recepción tica es suficiente para obtener una imagen ecográfica. Direc-
tamente ligada a la atenuación, la penetración depende del
Figura 24-11. Resolución axial. La señal obtenida con un transductor
tejido explorado y decrece con la frecuencia. Para los tejidos
de factor de calidad Q elevado no permite separar los dos ecos
reflejados, contrariamente a un transductor de Q pequeño.
blandos, se obtiene la siguiente expresión aproximada:
40
que la intensidad emitida se reduce en un 85% al cabo de un penetración (cm) (24-15)
frecuencia (MHz)
tiempo τ. Si se desprecia la intensidad emitida tras ese retra-
so, la resolución axial viene dada por ∆t = τ, por lo que:
QC Compromiso resolución-profundidad
∆x = = Qλ 0 (24-12)
ν0 La elección de la frecuencia utilizada en ecografía resulta
La resolución axial es tanto más fina cuanto mayor es de un compromiso entre las resoluciones en profundidad y la-
la frecuencia de los US y el factor de calidad. Los ecógrafos teral, que mejoran cuando la frecuencia aumenta, y las profun-
actuales tienen una resolución axial de alrededor de 1 mm didades de penetración y de campo, que entonces disminuyen.
a 5 MHz. El compromiso más efectivo depende de los órganos examina-
dos. Se utiliza, por ejemplo, la gama de frecuencias siguientes:
Resolución lateral – 3.5 a 5 MHz: examen abdominal general, como cora-
zón, hígado, útero;
Se define como la distancia mínima que debe separar
– 5 a 10 MHz: tiroide, carótidas, mama, testículos, eco-
dos puntos situados a idéntica profundidad para que den
grafía pediátrica;
dos ecos separados (o para que uno dé lugar a un eco y el
– 10 a 15 MHz: ojo (de pequeño tamaño y bajo contraste
otro no, en el caso de un único transductor). Para un trans-
acústico).
ductor único no enfocado, la resolución lateral es máxima en
Las ecógrafos suelen estar equipados con más de una
la zona de Fresnel y aproximadamente igual al diámetro D
sonda, pudiendo trabajar cada una de ellas a una frecuencia
del transductor.
dada, adaptada a un tipo de examen.
Para una sonda de enfoque (de transductor único o de
enfoque electrónico), la resolución lateral es máxima a nivel
del foco, a distancia f y aproximadamente igual al diámetro Dinámica
del haz que viene dado por:
La dinámica es la característica que indica la posibilidad
fλ fC de detectar los ecos de muy baja amplitud y que correspon-
δ 2.44 0 = 2.44 (24-13)
D Dν0 den, por ejemplo, a las superficies de contacto, interfases,
entre dos medios de impedancias acústicas muy semejan-
tes. Se cuantifica en dB utilizando la relación de las inten-
Profundidad de campo
sidades del eco detectable más intenso Imáx y el del más bajo
Según (24-13), la obtención de una mejor resolución I
Imín: Dyn = 10.log máx . Para el transductor, esta relación es
lateral obliga a aumentar el diámetro D del transductor. I mín
Desafortunadamente, esta mejora se hace en detrimento de de alrededor de 70 a 80 dB (por tanto, una relación de in-
la profundidad de campo, que es la longitud π de la zona tensidades de 107 a 108). Después de diversos tratamientos
330 IMÁGENES
de las señales, la dinámica es de alrededor 40 a 50 dB. Los
sistemas de visualización tienen una dinámica de cerca de 20
a 30 dB, pero se benefician de posibilidades de segmentación
que permiten reproducir la gama de intensidades de la señal
que se desea estudiar.
Ganancia en profundidad t
nódulo
tiroides
vasos
sanguíneos
anchura tiempo
US
Figura 24-13. Diferentes modos de examen ecográfico. Se ha tomado como ejemplo la exploración de la región del tiroides (vista transversal, con
superposición de un nódulo y de la carótida cuyo modo TM muestra los latidos).
332 IMÁGENES
mero de líneas exploradas para una imagen de calidad es de Se obtienen haciendo girar una sonda multitransductora,
L = 100-200 líneas/imagen. bien sobre su propio eje, generando un barrido «circular»,
El factor limitante es la frecuencia de recurrencia: para bien alrededor del eje perpendicular al eje de la sonda (como
explorar una profundidad Z de un tejido que se mueve a una un abanico), generando un barrido «angular». En ambos ca-
velocidad C, el intervalo entre dos pulsos debe ser, como mí- sos, en cada momento de la rotación, se graba la imagen eco-
nimo, t = 2Z/C. Por tanto, se debe cumplir que 1/NL 2Z/C, gráfica de un plano. Estas diferentes imágenes se reagrupan
es decir, 2NLZ C. a continuación en formato de tabla numérica tridimensional
Por ejemplo, una sonda para una ecografía de tiroides que representará los valores del eco recogido en cada uno
(7.5 MHz) explora, a lo largo de una longitud de imagen de de los elementos del volumen explorado. Estos elementos,
8 cm barrida por L = 128 líneas, una profundidad de tejido de denominados vóxeles, se corresponden a los píxeles de una
Z = 6 cm. Como la velocidad es C = 1540 m.s–1, se deduce que: imagen digital bidimensional.
– la longitud real de cada línea explorada es l = 8 cm/ La utilización de los datos del volumen se puede llevar a
128 = 0.6 mm, orden de magnitud de la resolución lateral in- cabo de varias maneras:
trínseca del detector; – representación de uno o más planos de corte con
– la cadencia máxima de las imágenes es N = C/2LZ = orientación al azar (se puede, por ejemplo, explorar un plano
= 1540/ (2 128 0.06) 100 imágenes/segundo, lo que de corte perpendicular al eje principal del haz US, que es im-
resulta suficiente para esta aplicación. posible de conseguir en ecografía 2D clásica);
En una imagen ecográfica, las diferentes estructuras ex- – representación de una interfase en particular, lo que
ploradas se distinguen de varias maneras: resulta en una imagen que da la impresión de relieve (véase
– las superficies de separación entre dos zonas con dis- capítulo 26, Síntesis y resultado de imágenes tridimensiona-
tinta impedancia acústica generan un eco intenso y son per- les) o de las zonas cuyo eco se encuentra en un intervalo de
cibidas como contornos brillantes. El valor más pequeño del amplitudes determinado. Se puede de esta manera visuali-
coeficiente de reflexión que permite distinguir una interfase zar de manera muy realista la cara o las extremidades de un
es aproximadamente 0.01%; feto, o incluso las estructuras óseas, u órganos anecogénicos
– en función de su estructura celular, las estructuras (llenos de líquido). Por su lado, la utilización de la técnica
atravesadas generan una señal US retrodifundida (es decir, Doppler permite visualizar los vasos sanguíneos.
difundida hacia el transductor) más o menos intensa. Por El tiempo de adquisición de una imagen es de 3 a 10 se-
ejemplo, un quiste de tiroides lleno de líquido y sin células gundos aproximadamente, siendo el tiempo de su recons-
en su interior no genera ningún eco y aparece negro; se dice trucción de 10 a 30 segundos, dependiendo de la velocidad
que es anecogénico. Por el contrario, un adenoma sólido (tu- del ordenador utilizado.
mor benigno) genera una señal difusa, que puede ser más o La ecografía 3D permite, en algunos casos, un balance
menos intensa que el resto del tejido de la glándula (adeno- más preciso de malformaciones fetales. Por el momento,
ma hiperecogénico o hipoecogénico, respectivamente); sólo se aplica en obstetricia.
– la organización estructural a mayor escala no se tra-
duce directamente en la imagen pero aparece a menudo con
textura de mayor o menor definición; Utilización del efecto Doppler
– algunas estructuras son reconocidas gracias a su mo-
Las exploraciones ecográficas basadas en el efecto
vimiento (p. ej., pulso de las arterias) o por aparecer en varios
Doppler permiten cuantificar la velocidad de las células san-
cortes sucesivos.
guíneas en los vasos, utilizándose frecuentemente en la ex-
Con frecuencia aparecen artefactos inherentes a la técni-
ploración cardiovascular.
ca. Por ejemplo, una estructura de impedancia muy diferente
del resto de tejidos (p. ej., la tráquea) produce una reflexión
prácticamente total de la onda US y la desaparición de las Principio del efecto Doppler
estructuras situadas en la «zona de penumbra» por detrás de
la misma. Por el contrario, una estructura anecogénica que Cualquiera que se haya cruzado por la calle con una am-
no atenúe apenas los US (p. ej., un quiste) se traduce en un bulancia, ha experimentado el efecto Doppler al escuchar la
refuerzo de las estructuras subyacentes, cuyos ecos aparecen sirena haciéndose claramente más grave a medida que se ale-
amplificados por la ganancia de profundidad, puesto que los ja de nosotros. Supongamos un sonido de frecuencia ν0 refle-
haces incidente y reflejado sufren una menor atenuación.
u(x,t)
V
O x
Ecografía en 3D V
w(x,t)
La ecografía tridimensional utiliza sondas especiales
x X
que adquieren la señal ecográfica, no en dos dimensiones
(como las sondas multitransductoras clásicas), sino en tres. Figura 24-15. Principio del efecto Doppler.
334 IMÁGENES
modo B a tiempo real. Se puede, en la pantalla donde se vi- salvo en el caso de un examen prolongado de una misma
sualiza esta imagen, elegir una dirección y una distancia de zona. En el momento del pulso incidente US, la intensidad
examen Doppler, es decir, seleccionar finamente el origen de media llega a los 100 mW.mm–2 en la zona de enfoque, pero
la señal Doppler analizada. Este resultado se obtiene anali- solamente durante un tiempo muy breve. Se pueden produ-
zando únicamente la señal reflejada durante un intervalo de cir efectos de cavitación (que consiste en la creación de bur-
tiempo determinado, intervalo que corresponde a la dura- bujas de aire potencialmente muy peligrosas) producidos por
ción de ida y vuelta del pulso Doppler hasta la zona elegida. los US, pero para intensidades muy superiores a las utiliza-
La imagen ecográfica simultánea permite medir el ángu- das en imagen. Finalmente, no se ha detectado ningún efecto
lo θ, lo que brinda una medida absoluta de la velocidad cir- cancerígeno ni teratogénico. Las técnicas ecográficas pueden
culatoria, así como del diámetro del vaso sanguíneo, lo que ser consideradas, por tanto, como totalmente inocuas y por
permite estimar su caudal. ello son ampliamente utilizadas en Obstetricia. Se trata, sin
Es también posible superponer la imagen ecográfica embargo, de un análisis médico que debe estar debidamente
en gris a una «imagen funcional Doppler» que codifica para justificado, lo que no incluye algunas ecografías cuyo único
cada píxel la velocidad de un movimiento eventual detecta- objetivo es obtener para los padres una fotografía del feto.
do por el análisis Doppler. Los movimientos hacia la sonda se
representan en rojo, tanto más intenso cuanto más rápidos
son dichos movimientos, mientras que en azul se represen- Conclusión
tan los movimientos que se alejan del detector y en verde los
La ecografía ha alcanzado un lugar principal en nume-
inestables que corresponden a turbulencias. Los tejidos in-
rosos campos de la medicina, particularmente en obstetri-
móviles no se ven modificados y aparecen en gris. Se obtiene
cia, cardiología, gastroenterología y endocrinología. Su ino-
así, en tiempo real, una representación de los vasos de una
cuidad, su excelente resolución espacial, su poder de discri-
región que contendrá información morfológica y funcional,
minación de diferentes tejidos blandos, las posibilidades de
por ejemplo, la región cervical de la carótida o de la vascula-
análisis Doppler, se suman a su facilidad de puesta en mar-
rización de un nódulo tiroideo.
cha y al precio razonable de los aparatos (en comparación
Técnica muy espectacular y útil, el Doppler pulsado
con otras modalidades de imagen) y por todo ello su impor-
tiene sus límites. No permite, por ejemplo, analizar los mo-
tancia es evidente.
vimientos demasiado rápidos de los flujos (apareciendo, en
Hay que señalar, no obstante, que se trata de la más sub-
este caso, un fenómeno denominado aliasing, de naturaleza
jetiva de las modalidades de imagen médica. Una ecografía
estroboscópica, y que conduce a una estimación errónea de
sólo tiene el valor que ofrezca el ecografista que la hace (la
las velocidades). Sigue siendo modesta la calidad de la ima-
escala de competencia tiene muchos niveles) y, además, las
gen ecográfica asociada, en la que el número de barridos está
imágenes que acompañan a los informes no son más que un
necesariamente limitado para dejar tiempo para el análisis
resumen sucinto de lo que ha visto el facultativo, resumen
Doppler. No permite la exploración de los vasos sanguíneos
por lo general insuficiente para ser interpretado por otro
profundos, ya que aumentar la profundidad conduce a espa-
ecografista.
ciar los pulsos para dejar a los US el tiempo de ida y vuelta,
La tecnología de la ecografía está progresando rápida-
con aparición de artefactos de aliasing.
mente, con mejoras constantes de la resolución espacial, de
sistemas informáticos más eficientes y finalmente, de la eco-
Doppler en modo de energía grafía 3D, cuyo interés está por demostrar.
Después de eliminar los ecos provenientes de las estruc-
turas fijas, la función de autocorrelación de la señal Doppler
permite calcular su energía, independientemente de la ve-
Ejercicios
locidad de las partículas que la han originado, por tanto del Ejercicio 24-1. Se interpone una lámina de material de impe-
ángulo que forma la dirección del flujo con la del haz US in- dancia Z entre dos medios materiales M1 y M2 de impedan-
cidente. La representación en la pantalla de esta medida de cias acústicas Z1 y Z2. ¿Cuál tiene que ser el valor de Z para
la energía en cada punto codificada en pseudocolor muestra que el paso de US de M1 a M2 se produzca con intensidad
un conjunto de flujos circulantes con una elevada sensibili- máxima? Se tendrán en cuenta las reflexiones de los US, pero
dad (incluso se pueden ver algunos flujos capilares). consideraremos despreciable su absorción.
336 IMÁGENES
Imagen de centelleo
25
El principio de la imagen de centelleo se basa en la vi- tructuras anatómicas o funcionales en las que dicho reparto
sualización de una estructura anatómica o funcional del or- sea diferente.
ganismo por medio de la localización de un isótopo radiac- Las cámaras utilizadas en centellografía (cámara gamma,
tivo susceptible de una detección externa (figura 25-1). Es tomógrafo de emisión de positrones o PET) se encuentran,
por tanto una técnica de imagen por emisión. Los detectores cada vez con mayor frecuencia, acopladas a un escáner, lo que
utilizan el principio de centelleo sólido, de aquí el nombre permite disponer en la misma prueba de información funcio-
de «centellografía» dado a este tipo de imagen(1). nal junto con información morfológica de buena calidad.
Las imágenes de centelleo son de calidad morfológica
mediocre, comparadas a las que aportan la tomografía clá-
sica, la tomodensitometría o la IRM. Por su parte, la imagen
Trazadores y marcadores
de centelleo conlleva un poder diagnóstico considerable de- La imagen de centelleo utiliza radiofármacos consti-
bido a su carácter dinámico y funcional: permite hacer un tuidos por un trazador marcado con un isótopo radiactivo.
seguimiento de la evolución del reparto en el organismo de Muy frecuentemente, estas sustancias se introducen en el
una determinada sustancia, y por tanto diferenciar las es- organismo por vía venosa, aunque a veces se utiliza la vía
oral o la respiratoria.
órgano
a visualizar Trazadores
Se llama trazador a una sustancia que puede localizarse
de manera selectiva en una estructura del organismo en par-
ticular (p. ej., un órgano, un compartimiento fluido o una le-
sión). El trazador puede ser un elemento simple (yodo, xenón),
sujeto detector dispositivo
de formación
una molécula (albúmina, pirofosfato) o tener una estructura
de la imagen más compleja (células, macroagregados de albúmina).
Figura 25-1. Principio de la imagen gammagráfica. La presencia del trazador en el seno de la estructura que
se desea visualizar resulta de alguno de los siguientes fenó-
(1) El uso del término centellografía es limitado en la práctica mé-
menos:
dica. Asimismo, el anglicismo «escintigrafía» se utiliza esporádicamente – fenómenos pasivos puramente físicos, como bloqueo
para referirse a este tipo de imagen médica. Lo más habitual es hablar de macroagregados en los capilares pulmonares, o también
de gammagrafía, atendiendo al origen de la emisión y fundamento de visualización de las cavidades cardíacas por medio de un
la técnica. Sin embargo también en la práctica se distingue entre la trazador vascular;
gammagrafía clásica y la tomocentellografía de fotón único (SPECT) o la
– fenómenos metabólicos activos, como captación de
tomografía de emisión de positrones (PET). Aquí se usará centellografía
para las dos últimas, así como para la descripción de las características yodo por parte del tiroides, del pirofosfato por los osteoblas-
comunes a todas ellas, reservando el término de gammagrafía para las tos en las zonas de reconstrucción ósea o captación del talio,
pruebas clásicas. (N. del T.). análoga a la del potasio, por las células del miocardio;
337
– fenómenos excretores. Por ejemplo, numerosos tra-
Tabla 25-I Principales marcadores usados
zadores son filtrados por el riñón y concentrados en la ori-
en centellografía
na, mientras que otros son secretados en el túbulo renal. Se
puede de esta manera visualizar el parénquima renal y las Elemento Período Tipo de emisión y energía (keV)
vías urinarias; C
11
20 min β+ ⇒ 2ϕ – 511 KeV
– reacción antígeno-anticuerpo específica. Se utiliza
para este fin anticuerpos monoclonales marcados, dirigidos N
13
10 min β+ ⇒ 2ϕ – 511 KeV
contra antígenos expresados en la superficie de las células O
15
2 min β+ ⇒ 2ϕ – 511 KeV
tumorales como, por ejemplo, de cáncer de ovario, de apa-
F
18
112 min β+ ⇒ 2ϕ – 511 KeV
rato digestivo, de médula ósea o del tiroides.
Algunos trazadores se fijan selectivamente en las zonas Cr
51
27.8 d γ – 320 KeV (10%)
sanas y dan imágenes «negativas» de las zonas anormales Ga
67
78 h γ – 92 KeV (39%)
(hipofuncionales, mal vascularizadas o cancerosas, depen- γ – 184 KeV (21%)
diendo del trazador). Inversamente, algunos trazadores se fi- γ – 300 KeV (17%)
jan de manera más o menos específica a las zonas anormales e– – 84 KeV (29%)
y ofrecen imágenes «positivas», mucho menos contrastadas. 99m
Tc 6h γ – 140 KeV
La fijación del yodo es singular, ya que los niveles basa-
e– – 119 KeV (9%)
les en la glándula normal son de por sí elevados, disminu-
yendo, pero sin anularse, en el cáncer de tiroides y sus me-
111
In 2.8 d γ – 245 KeV (94%)
tástasis. Un nódulo canceroso de tiroides aparece por tanto γ – 171 KeV (90%)
como una zona de hipofijación (nódulo «frío»). Después de 113m
In 1.7 h γ – 390 KeV (64%)
una extracción quirúrgica total del tiroides, las posibles me- e– – 364 KeV (29%)
tástasis fijan el yodo débilmente. Aparecen por tanto con un 123
I 13 h γ – 160 KeV (84%)
contraste positivo en las zonas que no lo fijan normalmente,
e– – 127 KeV (14%)
como puede ser el caso de los pulmones.
El contraste entre las zonas de estudio puede resultar 131
I 8d γ – 364 KeV (82%)
igualmente de una distinta cinética de fijación entre las zo- β – 610 KeV (87%)
nas sanas y las patológicas. Estas interrupciones dinámi- 133
Xe 5.3 d γ – 81 KeV (38%)
cas de la homogeneidad se evidencian tras realizar análisis β – 340 KeV (100%)
a tiempos determinados, en los que la discriminación es 201
Tl 3d γ – 167 KeV (10%)
óptima, o bien realizando estudios dinámicos fundados en
X – 71 KeV (47%)
el análisis de series temporales de imágenes.
e – 84 KeV (16%)
Marcadores que se desee visualizar, lo que puede durar unas cuantas ho-
ras. Finalmente, el marcador no debe presentar toxicidad y
Al trazador se le acopla un marcador, un radionucleido, ser fácilmente disponible.
que permite hacer el seguimiento del destino del trazador en La tabla 25-I muestra las características físicas de los
el organismo así como de cuantificar su concentración local marcadores más frecuentemente utilizados en gammagrafía
para obtener índices funcionales. El marcador puede ser y otros tipos de centellografías. Sólo se indican las emisiones
el propio trazador (p. ej., 123I utilizado en las gammagrafías significativas. Se pueden agrupar en tres familias:
de tiroides) o encontrarse fijado al trazador por sustitución – los emisores (β–, γ) para los cuales la emisión β–, que
(p. ej., se marcan las proteínas con yodo radiactivo, sustitu- no puede ser detectada por un dispositivo externo, no tiene
yendo un H por un I en las cadenas laterales de tirosina), el más mínimo interés, pero aumenta inútilmente la dosis
o por quelación (asociación de un ligando con un radionu- de irradiación del paciente (131I, 133Xe);
cleido metálico como el 99mTc). – los emisores γ puros, por captura electrónica (67Ga,
Para someterse a una detección externa, el marcador 123
I), o metaestables (99mTc, 113mIn). Estas emisiones, conoci-
debe emitir fotones γ cuya energía esté dentro de los ran- das como «γ puras» se acompañan en realidad de una emi-
gos de los detectores habitualmente utilizados (entre 80 y sión no despreciable de electrones e–, que no provienen del
300 keV para gammacámaras, con un óptimo alrededor de núcleo, sino de la corteza electrónica por el fenómeno de
150 keV). El período de semidesintegración físico del mar- conversión interna (véase capítulo 14, Emisión β– y (β–, γ)).
cador debe ser lo suficientemente corto para que la activi- El 99mTc es el isótopo más utilizado en gammagrafía. Presen-
dad necesaria para obtener una imagen de buena calidad no ta, en efecto, unas propiedades físicas óptimas (período de
conlleve una dosis de irradiación excesiva (véase capítulo semidesintegración de 6 horas, energía igual a 140 keV) y
17, Dosimetría in vivo); debe ser lo suficientemente largo sus propiedades químicas permiten combinarlo con un nú-
para permitir la fijación del trazador en el órgano o la lesión mero importante de trazadores de interés biológico. A pesar
338 IMÁGENES
de su corto período de semidesintegración, se puede utili- Por ejemplo:
zar fácilmente en cualquier centro: en efecto, se prepara a 31 γ)
P ⎯(n,⎯→ 32 P
medida que se necesita, a partir del núcleo «padre», el 99Mo,
γ)
cuyo período de semidesintegración, de 67 horas, permite
130
Te ⎯(n,⎯→ 131 Te
tener una provisión suficiente para varios días (véase más El interés de la segunda reacción proviene de que el
adelante, Generadores de radionucleidos); 131
Te, emisor β–, se transforma en 131I, radionucleido muy uti-
– los emisores β+ (11C, 13N, 15O, 18F). Estos marcadores se lizado.
detectan indirectamente mediante los dos fotones de aniqui-
lación de 511 keV, emitidos en direcciones opuestas (véase Radionucleidos producidos en un acelerador de partículas
capítulo 14, Emision β+ y (β+, γ)). La utilización de los tres pri-
meros se reserva a centros altamente especializados, debido Los aceleradores de partículas de tipo ciclotrón permi-
a su muy breve período de semidesintegración que requiere ten obtener partículas pesadas cargadas (protones, deutero-
su producción en un ciclotrón en el mismo lugar de su uti- nes, etc.) con alta energía cinética asociada y susceptibles de
lización. Estos isótopos de los tres elementos mayoritarios producir reacciones nucleares cuando son dirigidas hacia
en los organismos vivos, permiten marcar un gran número una diana formada por un elemento puro. Las reacciones
de moléculas de interés biológico, como glucosa, aminoáci- nucleares que se producen entre las partículas incidentes y
dos o catecolaminas, así como diferentes medicamentos de los núcleos de los átomos de la diana dan lugar a nucleidos
cuyo metabolismo y distribución puede hacerse un segui- que presentan un exceso de cargas positivas y que son muy
miento in vivo. El período de semidesintegración del 18F (112 frecuentemente emisores β+, o emisores γ tras captura elec-
min) permite utilizarlo a cierta distancia del lugar de su pro- trónica.
ducción, como marcador de la desoxiglucosa. Los foto- Por ejemplo:
nes emitidos, muy energéticos (511 keV), son detectados 11
B ⎯(p,⎯
n)
→ 11 C (emisor β+)
con cámaras especiales (PET), que se estudiarán más ade- 122
Te ⎯(d,⎯
n)
→ 123 I (captura electrónica, emisor γ)
lante (véase Utilización de emisores de positrones: PET y
PET-CT).
Generadores de radionucleidos
Producción de radionucleidos utilizados Estos generadores utilizan el principio de las series ra-
en centellografía diactivas (véase capítulo 14, Cinética de las series radiacti-
vas). La mezcla radionucleido padre-radionucleido hijo está
Los radionucleidos utilizados en centellografía se pro- contenida en un dispositivo que permite la extracción re-
ducen artificialmente. Se originan bien en reactores nuclea- petida del radionucleido hijo, producido permanentemente
res o en aceleradores de partículas, bien como producto de por la transformación radiactiva del padre (figuras 25-2 y
una transformación radiactiva. 25-3). Los generadores permiten sacar un beneficio del pe-
ríodo de semidesintegración del nucleido padre, suficiente-
Radionucleidos producidos en reactores nucleares mente largo como para evitar la necesidad de provisión ex-
En un reactor nuclear, se producen reacciones en cesivamente frecuente, y del nucleido hijo, suficientemente
cadena en las que se generan neutrones, lo que permite breve como para limitar la irradiación del paciente.
obtener flujos de neutrones del orden de 1013 a 1014 partí-
culas. s–1.cm–2. Si se introduce en dicho flujo de neutrones
Frasco para la recogida
una diana constituida por un elemento puro AZ X, se pro-
de la solución de 99mTcO4Na
ducen reacciones nucleares entre el neutrón incidente y
el núcleo de uno de los átomos de la diana, reacción del
siguiente tipo:
A
X + 01 n → A +1
X+γ reacción
Z Z
columna de aluminio
denominada (n,γ)
o también:
A
Z X + 01 n → A
X + 11 p reacción
Z −1
denominada (n,p)
A
Estas reacciones dan lugar a isótopos nZ X→+A01+Zn
Z X + 0hijos
1 A 1
X→
+ó Zγ−A1 X ,+ 11 p
que presentan un exceso de neutrones y son, por lo general,
emisores β–. El rendimiento global de la reacción es muy pe- reservorio con la
queño; además, en el caso de la reacción (n,γ), el producto solución de elución
blindaje de plomo
final es químicamente idéntico al producto inicial, por lo
que no se pueden separar fácilmente. Figura 25-2. Esquema de un generador de 99mTc.
(λ1N1)0
actividad de 99Mo
actividad
de
99m
Tc ordeño
ordeño
ordeño
ordeño
12 24 36 48 60 72 84 96 108 120 132 144 156
Tiempo (horas)
Figura 25-3. Principio del ordeño de la «vaca» (obtención periódica) de 99mTc.
340 IMÁGENES
blindaje exterior
selector
de amplitud x
z e
cálculo x, y, z x
y
y
+1
fotomultiplicadores
guía de la luz imagen digital
cristal de centelleo
septum canal
colimador paralelo
emisión
Cabezal 1 C
ab
1
ez
al
al
ez
2
ab
C
1
Paciente
Camilla de examen
Soporte
Camilla Cabezal 2
de examen
Soporte
Figura 25-5. Disposición de las cámaras multidetectoras. (1) Barrido; (2) rotación para tomografía.
mum) de la función de dispersión obtenida para una fuente les y la distancia fuente-colimador (cuando se aleja la fuente
radiactiva puntual (véase capítulo 20, Correspondencia en- en el eje del canal, los fotones que emite pueden pasar por
tre la señal física y la imagen analógica, y figura 20-2). Para los canales vecinos, con una degradación de la resolución
un colimador paralelo (figura 25-6), se muestra que: espacial).
La sensibilidad SC o «eficacia geométrica» se cuantifica
d(l + b)
RC = FWHM (25-1) por la relación del flujo de fotones transmitidos por el coli-
l − 2µ −1
mador con respecto al flujo de los fotones emitidos. Para un
en la que d es el diámetro de los canales, l su longitud (gro- colimador paralelo, la sensibilidad viene dada por la ecua-
sor del colimador), b es la distancia fuente-colimador y µ el ción:
coeficiente lineal de atenuación del plomo (el término –2µ–1 2
⎛ Kd 2 ⎞⎟
toma en cuenta la penetración en los septos, su valor es de SC ⎜⎜ ⎟ (25-2)
–0.1 cm a 140 keV). RC aumenta con el diámetro de los cana- ⎜⎝ (1 − 2−1 )(d + t) ⎟⎟⎠
septum
canal F1 F2
l A B
F1 F2
cristal foco
d
b
cristal
F1 F2
interacción
Compton
foco
C D F1 F2
emisión
en la que K es un coeficiente dependiente de la forma de los nes ortogonales superpuestas, tanto más borrosas cuanto
orificios (orificios redondos: K = 0.24; orificios hexagonales: más alejada del colimador esté la fuente, siendo necesario
K = 0.26). SC aumenta con el diámetro de los canales y con el situar la zona estudiada lo más cerca posible del colimador.
área relativa que ocupan en el colimador. La sensibilidad no Como la sensibilidad no depende de la distancia fuente-co-
depende de la distancia de la fuente al colimador: si se aleja limador, se adaptan convenientemente a las pruebas cuan-
una fuente puntual en el eje de un canal, el ángulo sólido titativas (medida de la actividad contenida en una zona de
bajo el cual se ve el orificio de salida del canal disminuye, la imagen). La tabla 25-II muestra las características de al-
pero los fotones emitidos por la fuente pueden pasar por gunos colimadores de orificios paralelos.
los canales vecinos y esos dos fenómenos se compensan. Los colimadores convergentes (figura 25-7C) y, menos
La función de dispersión espacial de una fuente puntual se frecuentemente, los divergentes (figura 25-7D) permiten
amplía por tanto cuando la distancia al colimador aumenta, adaptar el tamaño de la zona de detección al del área explo-
pero el área bajo la curva permanece constante. rada. Los colimadores en abanico de enfoque lineal o fan
La elección de un colimador depende del radionucleido beam, sólo convergen en una dirección. Al utilizar óptima-
utilizado y de la información que se espera de la prueba. Di-
cha elección está basada en los siguientes parámetros: Tabla 25-II Colimadores de orificios paralelos
– la energía de los fotones a detectar, clasificados como
Tipo Energía Resolución Sensibilidad
de baja energía (BE ≤ 140 keV); energía media (140 < ME ≤
máxima a 10 cm
400 keV) y de alta energía (HE > 400 keV). El grosor de los de los fotones (mm)
septos, previsto para detener el 95% de los fotones parásitos, (keV)
va desde algunas décimas de milímetro de plomo para los
colimadores adaptados a las bajas energías, a varios milíme- BE-AR 150 7.4 1.84 10–4
tros para las energías más altas; Baja energía
Alta resolución
– un compromiso entre la sensibilidad y la resolución
espacial. Los colimadores de canales gruesos dejan pasar BE-MU 150 9.1 2.68 10–4
muchos fotones pero no permiten una localización muy Baja energía
precisa de su origen. Tienen por lo tanto una buena sensibi- Multiuso
lidad pero una resolución mediocre. Los colimadores de ca- BE-AS 150 13.2 5.74 10–4
nales finos, tienen, contrariamente, baja sensibilidad pero Baja energía
buena resolución. Alta sensibilidad
Los colimadores de canales paralelos (figura 25-7A)
ME-AS 400 13.4 1.72 10–4
son los más comúnmente utilizados. Los orificios son circu-
Media energía
lares o hexagonales. Los colimadores de alta calidad se pro-
Alta sensibilidad
ducen mediante fundición en moldes. Producen proyeccio-
342 IMÁGENES
mente el área de los detectores de gran tamaño, mejoran la
calidad de las gammagrafías cerebrales o cardíacas. Brindan
imágenes deformadas (en proyección «astigmática») pero
sólo son utilizados para llevar a cabo tomografías, por lo que
este problema no tiene excesiva importancia.
Los colimadores estenopeicos (pin-hole) (figura 25-7B)
están constituidos por un cono de plomo perforado con una
pupila de entrada en su vértice. Los colimadores pin-hole
producen imágenes de proyección cónica (producen erro-
res de paralaje), con efecto de agrandamiento y excelente re-
solución espacial. Su sensibilidad es proporcional al área del
orificio, más bien pequeña, disminuyendo muy rápidamen-
te a medida que aumenta la distancia fuente-colimador. Es-
tos colimadores se utilizan para la gammagrafía de órganos
superficiales de pequeño tamaño como el tiroides.
Figura 25-8. Disposición de los fotomultiplicadores (PM) de una
Cualquiera que sea el tipo de colimador utilizado, me- gammacámara circular. Los PM, cuya cara de entrada es hexagonal,
nos de un fotón de cada 1000 emitidos contribuye a la for- están dispuestos de manera ordenada según una red que desborda
mación de la imagen. Este desperdicio inevitable explica la en algunos centímetros la proyección del cristal de NaI (representado
baja resolución espacial y el pobre valor morfológico de las en gris).
imágenes centellográficas, cuyo interés es fundamental-
mente funcional.
– determinar la energía total transferida al cristal por
el fotón γ incidente, sumándose los resultados parciales
Detector de centelleo de cada PM. Un selector de amplitud (véase figura 16-12)
permite considerar únicamente los fotones cuya energía
El principio del detector de centelleo se ha estudiado transferida se sitúa en una ventana que rodea el pico de
en el capítulo 16 (véase Detectores de centelleo sólido). El absorción total, lo que permite rechazar la mayoría de los
de la gammacámara está constituido por un cristal (yodu- fotones difundidos. Con frecuencia, la anchura de la venta-
ro de sodio «dopado» con talio) de gran tamaño, circular o na sobrepasa el 10% de la energía del fotón incidente (p. ej.,
rectangular, pudiendo alcanzar 60 50 cm2 y 3/8 de pulgada 140 ± 14 keV para el 99mTc). Para algunos isótopos (67Ga, 111In),
(9.5 mm) de grosor. Este grosor permite la absorción comple- se pueden utilizar simultáneamente dos o tres ventanas de
ta y por tanto la detección de cerca del 90% de los fotones recuento para detectar los diferentes tipos de fotones γ emi-
140 keV emitidos por el 99mTc y de alrededor del 30% de los tidos. Se puede igualmente utilizar de manera simultánea
fotones de 364 keV emitidos por el 131I. Un cristal más grueso varias ventanas de recuento para pruebas efectuadas con
permite un mayor rendimiento de detección de los fotones dos marcadores diferentes (estudios con varios isótopos);
muy energéticos, pero a costa de una degradación de la re- – calcular la abscisa y la ordenada del punto de impacto
solución espacial. Los fotones de centelleo tienen una longi- del fotón γ incidente (figura 25-9). El método de determina-
tud de onda λ = 410 nm y una energía de alrededor 3 eV; son ción, enteramente digitalizado, se basa en la comparación
fotones luminosos, de color azul. El centelleo se atenúa con de las amplitudes de los pulsos generados en cada PM;
una constante de tiempo muy breve (250 ns), que permite
a priori la detección individual de fotones que llegan a ritmo
rápido. Un fotón de 140 keV crea alrededor de 4500 fotones
1 500 200 200
de centelleo (es decir, alrededor de 30 eV por cada fotón lu-
minoso creado).
1 500 1 500 200 200
Alrededor del 80% de los fotones de centelleo creados
por la interacción entre el fotón incidente y el cristal son re-
200 200
cogidos por intermedio de una «guía de luz» sobre una bate- 3300
ría de fotomultiplicadores o PM (hasta 90, dependiendo de 750 1 500
los diferentes modelos de cámara), dispuestos en una red
hexagonal (figura 25-8). Cada fotón de centelleo recibido en 1 500 750
un PM da lugar a un electrón que se «multiplica» por una
serie de dínodos, dando lugar a su vez a una avalancha de
Figura 25-9. Localización del punto de impacto detectado. La figura
electrones sobre el ánodo del PM (véase capítulo 16, Detec-
muestra tres ejemplos del número de fotones detectados por cada
tores). Esta cantidad de electricidad, que puede llegar a 1010 PM en función de los puntos de impacto de tres fotones incidentes
electrones por fotón de centelleo, es recogida y amplificada sucesivos (puntos negros). El análisis de las cantidades de electricidad
por un mecanismo electrónico que permite: resultantes permite localizar el punto de impacto.
e = 2.5 cm
20
e = 1.25 cm
10
5 e = 0.6 cm
foton
incidente 0 100 200 300 400 500 600 700
energía de los fotones (keV)
Figura 25-10. Error en la determinación del punto de impacto de un
fotón debido a una interacción Compton en el cristal de centelleo. Figura 25-11. Sensibilidad intrínseca de un detector de una
La suma de las energías transferidas en la interacción Compton (1) gammacámara en función del grosor del cristal y de la energía
y de la absorción total del fotón Compton de retroceso (2) es igual a de los fotones.
la energía del fotón incidente. El análisis espectrométrico acepta por
tanto dicha detección cuyo punto de impacto se localiza de manera cristal, de la guía de luz, del buen funcionamiento y equi-
errónea entre (1) y (2).
librio de los PM. Esta cualidad es esencial para la calidad
– efectuar una calibración automática e individual de de la prueba: se podría dar un diagnóstico de fijación baja
los PM. del trazador, cuando lo que existe es una avería del PM, por
El conjunto del detector de centelleo se caracteriza por su ejemplo. Controlada diariamente con una fuente ampliada,
resolución espacial RI denominada intrínseca, su sensibilidad la uniformidad se estima a partir de los valores mínimo y
intrínseca SI, su linearidad espacial, su uniformidad intrínse- máximo obtenidos, que deberían ser idénticos, si la unifor-
ca, su resolución de energía y su rendimiento en las tasas de midad fuera perfecta. Se utiliza el índice de uniformidad:
recuento.
1 máx − mín
La resolución intrínseca RI, cuantificada por la FWHM U= (25-3)
2 máx + mín
(«anchura a media altura») de la función de dispersión de
una fuente puntual, está limitada por dos fenómenos: U se encuentra normalmente entre el 2 al 4%.
– en el seno del cristal, un fotón puede ser desviado por La resolución de energía es la capacidad de cuantificar
efecto Compton antes de su absorción total por efecto foto- correctamente el valor de la energía total transferida al cris-
eléctrico (figura 25-10). Se localiza entonces de manera erró- tal, punto importante para diferenciar los fotones directos
nea por la cámara. El uso de cristales delgados limita este fenó- de los que han sufrido una interacción fuera del detector. Se
meno, pero entonces disminuye la sensibilidad del detector; evalúa a través del valor relativo de la menor diferencia de
– el reparto de fotones de centelleo entre los diferentes energía ∆E detectable, o por la anchura del pico fotoeléctrico
PM sufre fluctuaciones estadísticas tanto más importantes a media altura. La relación ∆E/E es de alrededor de un 10%.
(relativamente), cuanto menos numerosos sean los fotones. El rendimiento de la tasa de recuento es una caracte-
Los de energía superior a 100 keV permiten obtener una re- rística esencial para algunas aplicaciones que hacen uso de
solución intrínseca de 3 a 4 mm. La resolución espacial es actividades elevadas, en particular en cardiología.
peor si se utilizan fotones menos energéticos, que producen En los sistemas más comunes, cada detector «trata» cada
a su vez pocos fotones de centelleo. vez únicamente un evento, es decir, un fotón γ. Teniendo en
La sensibilidad intrínseca (o eficacia intrínseca) se cuenta el tiempo de extinción del centelleo (constante de
cuantifica por la relación SI del número de fotones detecta- tiempo 250 ns) y el tiempo de respuesta de los PM, el límite
dos con respecto al número que alcanza el cristal. Depende teórico es de alrededor 106 eventos por segundo (s–1), En rea-
de la malla del detector, pero sobre todo del grosor del cris- lidad, el número de eventos contados para una fuente de ac-
tal (figura 25-11) y de la anchura de la ventana de recuento. tividad creciente pasa por un máximo, y después disminuye.
La linearidad espacial es la capacidad de generar una Esta disminución proviene de la interacción simultánea de
imagen lineal a partir de fuentes lineales. Las gammacáma- varios fotones incidentes con el cristal: las energías de los di-
ras pueden presentar distorsiones geométricas cóncavas o ferentes fotones se suman en el sistema espectrométrico, so-
convexas, que dañan la uniformidad intrínseca y producen brepasando el límite superior de la ventana de recuento, por
errores de cuantificación. lo que ninguno de ellos es contabilizado. El rendimiento del
La uniformidad intrínseca es la capacidad de gene- detector se cuantifica por el número de eventos contados con
rar una imagen uniforme a partir de una fuente extendida una pérdida del 20% (alrededor de 100 000.s–1) y por el núme-
uniforme. Depende de la linearidad espacial, del estado del ro máximo de eventos contados (de 100 000 a 400 000.s–1).
344 IMÁGENES
Algunos sistemas permiten tomar en cuenta, en parale-
lo, varios impactos de fotones γ que se producen en puntos
suficientemente distantes del cristal y afectan a PM distin-
tos. A costa de una mayor complejidad del sistema electró-
nico de detección, se obtiene de esta manera una mejora
muy significativa de la tasa de recuento.
A B
La resolución espacial y la sensibilidad varían en senti-
do inverso cuando se modifica el grosor del cristal. Con 3/8 Figura 25-12. Ruido de origen estadístico. En las imágenes (A)
y (B), las elipses ofrecen el mismo contraste con el fondo, pero el
de pulgada (9.5 mm), obtenemos una situación de compro-
nivel de detección global es 8 veces superior para (A) que para (B).
miso satisfactorio entre un cristal demasiado fino (baja sen- La distribución de Poisson convierte a las imágenes en «ruidosas»,
sibilidad de detección) y un cristal demasiado grueso (me- circunstancia mucho más aparente en (B) cuyos contornos se
diocre resolución espacial). Los rendimientos intrínsecos visualizan pobremente.
de las cámaras de centelleo son óptimas para fotones cuya
energía esté comprendida entre 100 y 200 keV. Por debajo, la ble aleatoria que sigue una distribución de Poisson (véase
resolución es mala, y por encima, la sensibilidad es baja. anexo 7) de media M(l,c) = A τ H(l,c), en la que A es
la actividad inyectada, τ el tiempo de adquisición y H(l,c) la
Eficacia global del detector, calidad imagen «ideal» que se obtendría en ausencia de toda fluc-
de las imágenes tuación de la emisión radiactiva o de cada una de los eslabo-
Los rendimientos globales de una gammacámara se es- nes de la cadena de detección, esto es, si todo fuera perfecto.
tablecen en conjunto para el par colimador-detector especí- La varianza es igual a la media M(l,c) y la variación tipo a
ficamente utilizado para una imagen de centelleo. M(1,c) . Se cuantifica la importancia de las fluctuaciones y
La resolución espacial global RG, definida como la del desenfoque estadístico por la relación:
FWHM obtenida con una fuente puntual, está ligada a las variación tipo M(l,c) 1
resoluciones del colimador RC y la del detector RI por: σ(1,c) = = = (25-6)
media M(l,c) A × τ ×H(l,c)
R G R C2 + R I2 (25-4)
Se puede disminuir el desenfoque estadístico alargan-
Depende sobre todo de los rendimientos del colimador do la duración τ de la prueba siempre que su dinámica lo
y se degrada con la distancia a la fuente. A una distancia de permita, y administrando la actividad A máxima compatible
5 a 10 cm del colimador, que corresponde a las gammagra- con las limitaciones dosimétricas, de ahí el interés de utili-
fías más corrientes, la resolución global es de alrededor de zación de un marcador de vida corta.
7 a 12 mm con un colimador paralelo y de 5 a 8 mm con El desenfoque de origen geométrico proviene del he-
un colimador pin-hole. Estas cifras constituyen el límite cho de que la cámara no puede distinguir los fotones pro-
inferior del tamaño de las lesiones «frías» detectables, que venientes de dos fuentes puntuales muy próximas. La reso-
fijan menos trazador marcado que los tejidos circundantes. lución espacial es la distancia mínima que debe separar dos
Podemos visualizar las lesiones «calientes» más pequeñas, fuentes puntuales para que puedan ser diferenciadas. Las
pero su tamaño real es siempre sobrestimado. modificaciones susceptibles de mejorar la resolución espa-
La sensibilidad global SG es la probabilidad de que un cial, esencialmente la disminución del diámetro de los orifi-
fotón emitido alcance el detector y sea detectado efectiva- cios del colimador, se traducen en general por una disminu-
mente. Se muestra que: ción de la sensibilidad y por tanto por una agravamiento del
desenfoque de origen estadístico.
SG = SC SI (25-5)
La difusión de los fotones por efecto Compton en el pa-
Esta sensibilidad global depende principalmente del ciente arrastra una degradación de la información morfoló-
colimador. Puede expresarse en número de impactos detec- gica y cuantitativa contenida en la imagen centellográfica.
tados para una actividad dada medida durante un tiempo Se han propuesto técnicas eficaces para atenuar los efectos
dado. Es del orden de un fotón por cada 2000 emitidos. de la radiación dispersada. Se fundan por lo general en el
La calidad de las imágenes obtenidas resulta alterada estudio del espectro energético de los fotones detectados.
por el desenfoque de origen estadístico, el de origen geomé- La atenuación de los fotones emitidos por una zona
trico y la radiación difundida. Su valor cuantitativo se ve asi- del organismo cuando atraviesan el cuerpo es un obstácu-
mismo disminuido por los fenómenos de atenuación. lo para una cuantificación precisa, tanto más importante
El desenfoque de origen estático (figura 25-12) provie- cuanto más profunda es la emisión y menor energía tienen
ne del carácter aleatorio de la actividad detectada en cada los fotones γ, y por consiguiente son absorbidos en mayor
punto de una imagen. Un área de actividad homogénea número. Los métodos de corrección de la atenuación se
puede, de esta manera, dar lugar a una imagen falsamente basan en grabar dos imágenes mediante dos detectores
heterogénea. El número de «impactos por píxel», es decir, opuestos (véase ejercicio 25-7) o en el análisis del espectro
de fotones detectados en un píxel dado (l,c), es una varia- de emisión. En la tomogammagrafía de fotón único SPECT
346 IMÁGENES
una duración predefinida. Estos registros permiten hacer un
seguimiento de la dinámica de un trazador marcado en una
región dada, denominada región de interés (RDI, frecuente-
mente delimitada manualmente sobre las imágenes) y sa-
car información funcional. En la figura 25-16 se presenta
un ejemplo sacado de una gammagrafía dinámica de riñón
trasplantado. Este tipo de registro es el utilizado para llevar
a cabo las adquisiciones dinámicas sincronizadas con el
ECG.
Figura 25-15. Gammagrafía del esqueleto obtenida por barrido: cara anterior de frente. La gammagrafía real se ha reducido 4 veces.
Figura 25-16. Imágenes sacadas de una gammagrafía dinámica de riñón trasplantado, constituida por 60 imágenes obtenidas durante 20
segundos. La adquisición comienza desde el momento de la inyección del trazador, en este caso un diurético marcado con 99mTc y dura 20 minutos.
Se observa la actividad creciente y a continuación decreciente del riñón (captación y posterior excreción del trazador) y la actividad creciente
de la vejiga (abajo y a la izquierda de la imagen) en la que el trazador se acumula.
La determinación de la cartografía de la actividad por Figura 25-18. Corrección de la atenuación por una tomografía de
unidad de volumen V(x,y), conociendo los valores Fθ(ρ) en transmisión simultánea. Los cabezales E1 y E2 detectan clásicamente
un número finito T de direcciones θ, por un número finito R las emisiones de fotones γ que provienen del paciente (en el centro).
de valores de ρ, utiliza técnicas semejantes a las que hemos El tercer cabezal T está regulado para detectar únicamente los fotones
emitidos por la fuente radiactiva externa y que después de atravesar
presentado en el capítulo 20, pero es necesario introdu-
al paciente generan una imagen por transmisión que permite estimar
cir modificaciones para tener en cuenta las variaciones de y corregir la atenuación. El conjunto de tres cabezales y de la fuente
Ω (x,y,,) y de la forma de la zona Zθ(ρ). externa gira alrededor del paciente.
348 IMÁGENES
son importantes: los cabezales equipados de su colimador
Tabla 25-III Características de los principales emisores
son muy pesados y es imprescindible que el eje de rotación
de positrones
esté totalmente estabilizado. La resolución espacial de las
imágenes alcanza cerca de 9 mm. Isótopo Período de EC (keV) RLM (mm)
semidesintegración
Desarrollos futuros de la gammacámaras (min)
C
11
21 960 1.1
Destacan dos vías de mejora de las gammacámaras:
– el acoplamiento a un escáner de rayos X, para efectuar N
13
10 1200 1.4
en la misma prueba una imagen bimodal: el escáner aporta O
15
2 1730 1.5
una mejora considerable de la localización anatómica de las
anomalías gammagráficas y permite una buena corrección F
18
110 630 1,0
de la atenuación. Se comercializan modelos operacionales
de escáneres-gammacámaras. Es probable que esta mejora lulas donde, después de ser fosforilada, se acumula en forma
esencial llegue a generalizarse. de 2D-6P, puesto que no es sustrato de los enzimas de las
– la puesta a punto de nuevos detectores, más eficientes etapas siguientes de la glicólisis. La 2DG está marcada por la
que el par cristal de centelleo-batería de PM que alcanza sus sustitución de un grupo hidroxilo OH con 18F, generando de
límites en términos de tasa de recuento, resolución espacial esta manera el radiofármaco 18F-2-desoxiglucosa o 18FDG.
y energía. Las gammacámaras que utilizan semiconducto- La 18FDG es un trazador de las células con metabolismo
res como detectores presentan una vía muy prometedora: activado, lo que ocurre en los músculos durante el ejercicio, en
alta sensibilidad (30 000 cargas creadas por fotón γ frente a el corazón, en el cerebro, y, hecho importante, en las células
1500 con el par cristal-PM), muy buena resolución espacial, tumorales malignas, que se podrán visualizar con la 18FDG.
tasa de recuento elevada, posibilidad de utilizar el mismo La centellografía con 18FDG muestra una excelente
receptor para la imagen por emisión (centellografía) y por sensibilidad y una elevada especificidad para la detección
transmisión (escáner de rayos X). de tejido tumoral y sus indicaciones son muy numerosas
Una vía adicional de investigación muy activa en el mo- en oncología para diferenciar los tumores malignos de los
mento actual es la puesta a punto de nuevos radiofármacos benignos, hacer balances preterapeúticos, seguimientos de
(receptores de membrana, metabolitos específicos, etc.). un tratamiento, diferenciar cicatrices de tumores residuales,
buscar recidivas e incluso descubrir el lugar primitivo, hasta
ese momento desconocido, de un cáncer puesto en eviden-
Utilización de los emisores cia por sus metástasis.
El campo de aplicación es muy extenso: cáncer pulmo-
de positrones: PET y PET-CT nar, mamario, de tiroides, linfomas, melanomas, cáncer di-
gestivo, de cabeza y cuello (ORL), o bien cáncer de origen
Emisores de positrones y FDG
desconocido. Las otras aplicaciones, en cardiología o en
Las propiedades de los emisores β+ se han estudiado en neurología, están mucho menos desarrolladas.
el capítulo 14 (véase Emisión β+ y (β+, γ)). Recordemos que
su singularidad se debe al destino del positrón que, tras ha- Detector de positrones especializados (PET)
ber agotado su energía cinética EC en colisiones múltiples
en menos de 10–9 s y a distancia media del lugar de su emi- La localización espacial del 18FDG introducido en el or-
sión (denominada recorrido libre medio), se combina con ganismo (por vía intravenosa y en reposo) es difícil, debido
un electrón en una reacción de aniquilación. Esta reacción a la elevada energía de los fotones de aniquilación que hay
da lugar a dos fotones de 511 keV emitidos en direcciones que detectar (511 keV). Han llegado a utilizarse las gamma-
opuestas. Estos dos fotones de aniquilación pueden ser fá- cámaras clásicas, provistas de cristales de gran tamaño, pero
cilmente detectados. la centellografía con 18FDG utiliza actualmente, de forma
Las características de los principales emisores de posi- casi exclusiva, cámaras especiales denominadas cámaras de
trones (11C, 13N, 15O,18F) quedan recogidas en la tabla 25-III. tomografía de emisión de positrones o PET (del inglés posi-
El recorrido libre medio de cada uno de ellos se refiere a un tron emission tomography).
recorrido en el agua. Permiten marcar un gran número de La cámara PET está constituida por una batería de de-
moléculas de interés biológico o de medicamentos, pero tectores dispuestos en anillos de gran tamaño, por los que
sólo el 18F tiene un período de semidesintegración compati- avanza el paciente en una cama móvil, durante la duración
ble con su uso a distancia del lugar de producción (e incluso de la prueba. Cada anillo lleva alrededor de 500 detecto-
en este caso hay que asegurar dos o tres entregas diarias), y res independientes. La detección utiliza la propiedad de la
es el único marcador de este tipo utilizado rutinariamente. coincidencia: sólo se utiliza la información proveniente de
El trazador es casi siempre la 2-desoxiglucosa (2DG), los pares de fotones detectados simultáneamente en dos
compuesto análogo de la glucosa. La 2DG penetra en las cé- detectores y que pueden ser interpretados como fotones de
detector
de coincidencias
∆t < 12 ns
fotón
de aniquilación 2
anillo
de detectores
Figura 25-19. Principio de la cámara de positrones. Se representan los dos fotones de aniquilación que provienen de la misma fuente. El detector de
coincidencia interpreta como proveniente de la desintegración del 18FDG la detección sobre dos de los detectores del anillo de dos fotones de 511 keV
en un intervalo de tiempo de menos de 12 ns.
aniquilación del mismo positrón (figura 25-19). Se deduce Los cristales utilizados no son de yoduro de sodio como en
que la fuente está por tanto situada en la línea que une los las gammacámaras, sino de sales que difieren por su densi-
detectores. A partir de un gran número de coincidencias, se dad y su eficiencia luminosa (que condiciona la sensibilidad)
puede reconstruir una imagen por técnicas clásicas de retro- y el tiempo de disminución del centelleo (que debe ser lo más
proyección o por cálculo iterativo. El sinograma, representa- corto posible para asegurar una tasa de recuento elevada). Las
ción tridimensional de la señal, se obtiene reagrupando las características de los cristales BGO, GSO y LSO utilizados en
coincidencias registradas en direcciones paralelas (es decir, PET quedan recogidas en la tabla 25-IV donde se comparan
correspondientes a direcciones paralelas del trayecto de los con las del NaI utilizado por las gammacámaras clásicas.
fotones de aniquilación, vertical por ejemplo); la suma de
todos ellas conduce a una proyección tipo de Radón.
Tabla 25-IV Características de los cristales BGO, GSO
Este modo de detección no requiere colimador, lo que
y LSO utilizados en PET
confiere a la cámara PET una sensibilidad aproximadamen-
te 10 000 veces mayor que la gammacámara clásica. Cuando Nombre Composición (g.cm–3) Tiempo Rendi-
la cámara PET está equipada con varios anillos, es posible de rela- miento(1) (%)
utilizar un colimador constituido por septos que separan jación (ns)
dichos anillos e imposibilitan la detección de coincidencias BGO Germanato 7.1 300 22
entre anillos, por tanto realizando así una tomografía 2D por de bismuto
anillo. Sin colimador, se realiza una tomografía 3D (figura
GSO Oxi-orto-silicato 6.7 60 20
25-20). El PET 3D es más sensible que el 2D. Sus algoritmos
de gadolinio
de reconstrucción de imágenes son más complejos. Sobre dopado con cerio
todo, está alterado por un aumento de la proporción de ra-
diación dispersa, que pasa del 10-20% en el PET 2D a 40-60% LSO Oxi-orto-silicato 7.4 40 75
en el 3D. Este fenómeno aumenta la proporción de coinci- de lutecio dopado
dencias fortuitas (detección simultánea de dos fotones de con cerio
511 keV pero que no provienen del mismo positrón) y au- NaI Yoduro de sodio 3.7 230 100
menta el tiempo muerto de la cámara PET. dopado con talio
Los detectores dispuestos en anillo son del mismo tipo (1) El rendimiento luminoso relativo (número de fotones de centelleo
que los cristales de centelleo-PM. Los cristales son gruesos, emitidos por fotón de aniquilación detectado) se expresa en
bien adaptados a la detección de fotones muy energéticos. porcentajes con respecto a un cristal de NaI de dimensiones idénticas.
350 IMÁGENES
PET 2D PET 3D
septum
Figura 25-20. PET 2D y 3D. La figura está representada en un plano de corte que pasa por el eje de los anillos de detección. Se ha supuesto que el
detector estaba formado por tres anillos. Las flechas de dos puntas corresponden a los trayectos de los fotones de aniquilación. Los trayectos en
línea discontinua de la PET 2D corresponden a eventos no detectados, habiendo sido detenidos por los septos, o uno, o los dos fotones.
La resolución espacial de un PET usado con 18F es del El escáner se utiliza con una pequeña intensidad del
orden de 4 mm, correspondiente a la FWHM de la «mancha» haz de rayos X, con el fin de evitar una radiación excesiva
obtenida en el caso de la imagen de una fuente puntual. para el paciente, teniendo en cuenta que el escáner es justa-
mente el origen de la mayor parte de la dosis total recibida
en la prueba. Las imágenes resultantes son de una calidad
Acoplamiento PET-CT claramente más baja que las obtenidas con el escáner clá-
sico, pero su calidad es suficiente para la localización de las
Estas cámaras están cada vez más acopladas a un escá- anomalías gammagráficas y para corregir la atenuación.
ner helicoidal que comporta de 2 a 16 detectores, formando
un único aparato denominado PET-CT (figura 25-21), que
genera imágenes bimodales que agrupan informaciones Evolución de los PET-CT
funcionales aportadas por la PET e informaciones morfo-
lógicas aportadas por el escáner. El escáner permite igual- La tomografía de emisión de positrones tiene un impac-
mente la corrección de fenómenos de atenuación. to considerable en la admisión de los pacientes aquejados
de cáncer. Se encuentra en rápida evolución en cuatro di-
recciones de investigación:
– mejora del detector: utilización de cristales más efi-
cientes, detectores semiconductores, uso del time-of-flight,
esto es, de la diferencia en el tiempo de impacto de los dos
fotones de aniquilación en sus respectivos detectores, lo que
ofrecería indicaciones sobre la posición exacta del origen de
estos fotones, mejora de los algoritmos de reconstrucción y
del resultado de la fusión de las imágenes PET y CT;
– puesta a punto de nuevos trazadores o de fármacos
marcados cuyo PET permita estudiar el metabolismo. Cite-
mos, como ejemplo, los radiofármacos dopaminérgicos o
serotoninérgicos para el estudio de la neurotransmisión, la
18
F-fluordopamina para la búsqueda de un feocromocitoma,
el 11C-hidroxitriptófano, para la de los carcinoides o la 11C-
eritromicina para el estudio de este antibiótico;
– uso de trazadores emisores de positrones distintos
del 18F, lo que requiere la proximidad inmediata tanto del
Figura 25-21. La primera carcasa contiene el escáner y la segunda el ciclotrón que produce dichos trazadores con períodos de
tomógrafo de emisión de positrones. semidesintegración muy cortos, como de instalaciones de
352 IMÁGENES
Tabla 25-V Irradiación según el tipo de prueba centellográfica realizada
Órgano Trazador Actividad Cuerpo entero Médula ósea Ovario Testículo Otros
explorado marcado inyectada (MBq) (mSv) (mGy) (mGy) (mGy)
Tiroides 123
I 4 0.03 0.03 0.03 0.03 Tiroides
14 mGy
Miocardio 201
Tl 100 6 6 16 16 Tiroides
25 mGy
Esqueleto MDP 500 1 4 2 1.5 Vejiga
99m
Tc 30 mGy
Riñones DMSA 100 0.4 1 0.5 0.3 Riñones
99m
Tc 20 mGy
PET-CT 18
FDG 300 7 3 3 3 Cerebro
11 mGy
Escáner 9(1) 9 9 9 9 mGy
Total 16 12 12 12 320 mGy
(1) La dosis generada por el escáner depende de la calidad de las imágenes que se quiera obtener. Varía alrededor de 9mSv, para imágenes de baja
calidad, a 19mSv para imágenes de «alta calidad».
tumor, búsqueda sistemática de metástasis óseas), reuma- ficios de 0.3 cm de diámetro y 2 cm de altura, para fuentes
tología (diagnóstico de una fractura no detectable radioló- situadas a 1 cm y a 10 cm de distancia respectivamente. Se
gicamente) o neurología (diagnóstico de patologías vascu- recuerda que ρPb = 11.3 g.cm–3.
lares, demencias, epilepsias). La puesta a punto constante
de nuevos trazadores marcados mejora la especificidad y el Ejercicio 25-4. Calcúlese la sensibilidad de un colimador
valor diagnóstico de estas pruebas. de plomo con tres orificios redondos de diámetro 0.3 cm y
La difusión, lenta pero evidente, de las cámaras de emi- de altura 2 cm, a 140 keV (µ/ρPb = 1.9 cm2.g–1) y a 360 keV
sión de positrones, frecuentemente acopladas a un escáner, (µ/ρPb = 0.25 cm2.g–1); ρPb = 11.3 g.cm–3. Se utilizarán los re-
debería permitir responder a la demanda creciente de este sultados del ejercicio 25-1.
tipo de examen. Ejercicio 25-5. Calcúlese el ángulo máximo de cambio de di-
rección por efecto Compton de los fotones emitidos por una
fuente 99mTc y dispersados dentro de un paciente que no son
Ejercicios detenidos por el colimador, ni por el análisis espectrométri-
Ejercicio 25-1. Determínese para un colimador paralelo de co con una ventana + 10%.
plomo, el grosor t de los septos que serán necesarios para Ejercicio 25-6. Una fuente de 5 µCi situada a 20 cm de un
detener el 95% de los fotones que alcanza el detector según detector con su colimador genera en 10 minutos una ima-
el recorrido de la figura mostrada a continuación, con altura gen que corresponde a 15 000 fotones detectados. Calcúlese
l = 2 cm, diámetro de los orificios d = 0.3 cm. Se hará el cálcu- la sensibilidad global del detector.
lo para los fotones de 140 keV (µ/ρPb = 1.9 cm2.g–1) y 360 keV
(µ/ρPb = 0.25 cm2.g–1). Se recuerda que ρPb = 11.3 g.cm–3. Ejercicio 25-7. Tres fuentes de 99mTc, con idéntica actividad
A y que suponemos puntuales, se sitúan a diferentes pro-
Detector
fundidades p1, p2 y p3, cuyos valores se desconocen, en un
recipiente lleno de agua de profundidad total conocida h. Se
genera una imagen de estas tres fuentes con una gammacá-
mara equipada con dos cabezales situados respectivamente
encima y debajo del recipiente. Cada fuente i (i = 1, 2, 3)
l d
da lugar, en un tiempo T, en la imagen superior (o inferior,
t
fotón respectivamente), una imagen de forma circular más o me-
nos extendida, y cuya suma de píxeles vale Csup(i) (respecti-
Ejercicio 25-2. ¿Cuál es el inconveniente de hacer una vamente Cinf(i). Demostrar que los valores de Csup(i) o Cinf(i)
gammagrafía con 99mTc y con un colimador de energía me-
dependen de i, pero que el valor de Csup(i) C inf (i) es inde-
dia? ¿Y con 131I y un colimador de baja energía?
pendiente de la fuente considerada. ¿Qué se puede deducir
Ejercicio 25-3. Calcúlese la resolución espacial a 140keV de en cuanto a las posibilidades de corrección de la atenuación
un colimador de plomo (µ/ρPb = 1.9 cm2.g–1), que tenga ori- a partir de dos registros opuestos?
Dejando aparte la obtención de imágenes «intrínse- estas técnicas han sido el objeto de un número considerable
camente» digitales como las conseguidas en el escáner, la de publicaciones, ninguna que nosotros conozcamos, aparte
tomografía, la IRM o la imagen por sustracción, el procesa- de los métodos de segmentación, ha demostrado su utilidad
miento (también denominado tratamiento) de las imágenes en la rutina clínica y no serán presentadas en el marco res-
médicas digitales puede resumirse en los seis objetivos si- tringido de este libro.
guientes: Utilizaremos la nomenclatura introducida en el capítulo
– la visualización de las imágenes en pseudocolor con 20 (véase Características de las imágenes digitales). Recorde-
efectos de zoom y de interpolación adaptados a un propósito mos que una imagen digital es una tabla de números con L
diagnóstico; líneas y C columnas cuyo «elemento» (o píxel) situado en la
– la atenuación del ruido y la eliminación de artefactos línea l y la columna c se designa (l,c), asignándole su valor
por «filtrado»; I(l,c). El punto de partida, el origen, de la imagen es la esqui-
– la representación de una serie espacial de imágenes na superior izquierda, las líneas se numeran de arriba a aba-
por la síntesis y reproducción de imágenes tridimensionales jo, y las columnas de izquierda a derecha (véase figura 20-9).
o multimodales;
– la síntesis de la información de una serie temporal
de imágenes por métodos como la imagen paramétrica o el Visualización de las imágenes
análisis de Fourier temporal;
– la estimación cuantitativa de los parámetros prove- La representación de las imágenes digitales asignando,
nientes de imágenes individuales (morfometría, análisis de por convención, un color a cada valor numérico de la imagen
textura) o de secuencias de imágenes (regiones de interés, se ha presentado en el capítulo 20 (véase Visualización de
identificación de modelos); una imagen digital). Hemos utilizado una función de corres-
– la compresión de los datos para su transmisión y ar- pondencia entre el valor de un píxel y su color de visualiza-
chivo. ción, codificada en tres componentes rojo, verde y azul:
Es imposible ser exhaustivo en el marco de este libro y
presentaremos solamente las técnicas principales que per- I(l,c) = ν → (R(ν), V(ν), A(ν)) (26-1)
miten alcanzar esos objetivos. En el caso de la una imagen visualizada en niveles de gris
Por lo demás, se consagran numerosos trabajos a la pues- (G), esta relación se escribe más sencillamente, utilizando el
ta a punto de sistemas informáticos que permiten análisis de nivel de luminancia G(ν) = R(ν) = V(ν) = A(ν):
alto nivel, con el objetivo de aportar una ayuda directa a la
I(l,c) = ν → G(ν) (26-2)
interpretación de las imágenes. Dichos sistemas recurren a
técnicas de segmentación de las imágenes (recorte en zonas Para simplificar, y limitándonos a este último caso, he-
«similares» según el criterio elegido), de reconocimiento de mos demostrado que la función ν → G(ν) podría tomar for-
las formas, de representación de los conocimientos médicos mas muy diferentes en función del resultado deseado: fun-
(p. ej., bajo la forma de sistemas expertos o de redes neurona- ción lineal, segmentación, función gamma. Estos métodos
les) y de puesta en marcha de dichos conocimientos. Si todas permiten visualizar preferentemente una u otra gama de
355
A B C
Figura 26-3. Ampliación de una parte de una imagen
con o sin interpolación (C y B respectivamente).
Umbral 0-255 Umbral 0-140
356 IMÁGENES
dirección (véase figura 21-13). Si los cortes de la serie inicial Métodos lineales y no lineales
están demasiado espaciados, es necesario intercalar líneas
suplementarias en la imagen reconstruida. Se obtienen los Los métodos de corrección del ruido de las imágenes di-
valores correspondientes por interpolación a partir de los gitales tienen como propósito atenuar las fluctuaciones esta-
valores conocidos. dísticas (gaussianas o poissonianas) de las señales grabadas.
Este tipo de transformación se llama alisado (smoothing) o
filtrado. Es importante resaltar que el smoothing tiene como
Atenuación del ruido contrapartida la mayor o menor atenuación de la nitidez de
los contornos y de los detalles de la imagen inicial.
Naturaleza del ruido Un método se dice que es lineal si la suma y el producto
se conservan en forma de un escalar (es decir, es equivalente
El valor de un píxel en una imagen digital es el resulta-
tratar dos imágenes por separado y posteriormente sumar-
do de una medida. Teniendo en cuenta la naturaleza de los
las, o hacer la suma de las dos imágenes iniciales y aplicar el
fenómenos físicos puestos en juego y de las imperfecciones
tratamiento posteriormente). Si las imágenes se denominan
de los sistemas de detección, el resultado no es perfecto, sino
A y B, y el método M(), la linealidad se expresa como:
que corresponde a una estimación del «valor verdadero» que
no podemos conocer. Se designa como ruido la separación M(A + kB) = M(A) + k.M(B) (26-5)
entre el valor estimado y el valor real. Dependiendo de las
Un método se dice no lineal si esta relación no es válida
modalidades, la palabra ruido refleja diferentes fenómenos.
siempre.
En radiología digital, el ruido abarca las fluctuaciones
del emisor de rayos X, las de la atenuación por los tejidos
blandos, las del detector de rayos X y las del conjunto de la
cadena de medida de la señal provista por dicho detector
Convolución
(amplificación, conversión analógico-digital). El valor de la
señal es una variable aleatoria de tipo gaussiana; por tanto, La convolución es una técnica de filtrado lineal muy
si el valor real es ν, la probabilidad de que la medida resulte utilizada. Hace intervenir un núcleo de convolución, o más-
en una media de valor m es: cara de convolución, que no es sino una tabla de coeficien-
tes reales de dimensiones impares: (2K +1) (2K + 1). Este
1 ⎛ (m − ν)2 ⎞⎟ núcleo, de orden K, puede interpretarse como una imagen
Pν (m) = exp⎜⎜⎜− ⎟⎟ (26-3)
σ 2π ⎝ 2σ 2 ⎟⎠ G(l,c) cuyos índices de línea l y de columna c varían entre –K
y +K. Habitualmente se utilizan núcleos 3 3 (K = 1); 5 5
La desviación tipo σ es proporcional a ν y por consi- (K = 2) ó 7 7 (K = 3). La tabla 26-I muestra un ejemplo de
guiente la relación señal/ruido varía como 1/ ν . núcleo de convolución 3 3.
En gammagrafía, el ruido está esencialmente ligado a La convolución de una imagen I, aplicándole el núcleo
las fluctuaciones poissonianas de la emisión radiactiva cu- de convolución G (figura 26-5), representado como , con-
yos fotones γ hay que detectar. Si ν R es el valor real del siste en sustituir el valor de cada píxel I (l,c) por un valor pon-
número de fotones que se deberían detectar durante el tiem- derado entre los valores de los píxeles vecinos y él mismo,
po de adquisición, la probabilidad de detectar realmente un ponderación cuyos coeficientes son los valores del núcleo de
número m N fotones es (véase anexo 7): convolución. Si la suma s de los coeficientes no es nula, el re-
sultado se divide por s para conservar la suma de los valores
νm −ν de la imagen inicial. El nuevo valor I G(l,c) del píxel (l,c)
Pν (m) = e (26-4)
m! será entonces:
La desviación estándar σ es igual a ν y la relación se- 1 K K
Image
Imagen V = ra+sb+tc
ra + sb + tc Image
Imagen
inicial
initiale ++ud+ve+wf
ud + ve + wf de convolución
convoluée
++xg+yh+zi
xg + yh + zi
Es evidente que la convolución es una operación lineal. ponderado por los píxeles vecinos (dándosele más o menos
Si algunos coeficientes son negativos, ya que el resultado de peso al píxel central inicial, en función de que dicho píxel
la operación puede ser un valor negativo para algunos píxe- (0,0) tenga un valor preponderante o no). Estos filtros, lla-
les, se puede entonces asignarles un valor nulo o tomar su mados de paso de banda bajo debido a que respetan las fre-
valor absoluto. cuencias bajas de la imagen, atenúan las fluctuaciones con
Esta decisión plantea un problema para las K líneas y co- un desenfoque más o menos importante de los contornos y
lumnas que forman el cuadro exterior de la imagen I, para las de los detalles (figuras 26-6B y 26-6C). Un caso particular es
cuales I(l + x,c + y) no está siempre definido (la máscara se el filtro gaussiano cuyos coeficientes se basan en los valores
sale de los límites de la imagen). Este problema puede resol- de una distribución gaussiana de dimensión 2:
verse de distintas maneras, la más frecuente, extendiendo la ⎛ x2 + y2 ⎞⎟
1 ⎜⎜−
imagen inicial por simetría alrededor de sus bordes. G(x,y) = exp ⎟⎟ (26-7)
2πσ 2 ⎜⎝ 2σ 2 ⎟⎠
Según los valores de K y de los coeficientes G(x,y) del nú-
cleo de convolución, el resultado obtenido es muy diferente Según los valores de σ, la convolución por un núcleo
y la atenuación de las fluctuaciones de los píxeles de la ima- gaussiano de orden K produce un filtrado paso de banda
gen no es más que una de las posibilidades de aplicación de bajo tanto más importante cuanto mayor es el valor de σ (ex-
la convolución. presado en píxeles) y de K. Estos núcleos paso de banda bajo
Los núcleos de convolución cuyos coeficientes son to- son separables, lo que permite un cálculo rápido de la convo-
dos positivos sustituyen el valor del píxel central por un valor lución (véase ejercicio 25-5).
A B C D
1 1 1 1 1
1 1 1 1 1
1 1 1 1
1 1 1 1 1 1 –1 –1 –1
1 1 1 25
1 1 1 1 1 –1 9 –1
9
1 1 1 –1 –1 –1
1 1 1 1 1
Figura 26-6. Convolución con diferentes núcleos. Los núcleos aplicados se indican debajo de cada imagen. A) Imagen inicial. B) Filtro de paso
de banda bajo. C) Filtro de paso de banda bajo. D) Filtro de paso de banda alto.
358 IMÁGENES
x x
Imagen inicial UH UV 2
VH + VV
2
–1 0 1 –1 –2 –1
–2 0 2 0 0 0
–1 0 1 1 2 1 Conexión 4 Conexión 8
Imagen inicial Imagen filtrada Se utilizan dos técnicas principales para generar imáge-
nes 3D.
Figura 26-11. Filtro morfológico (DEEEDD). D: dilatación; E: erosión.
La técnica del nivel de profundidad consiste en generar
una imagen 3D en la que cada línea corresponde a cada una
de las imágenes de la serie. Se fija la dirección de observación
y un umbral U. Para cada línea de la imagen 3D, el valor de
cada píxel se calcula a partir de la profundidad p en la que se
encuentra en la imagen 2D correspondiente a dicha línea, un
valor que tiene que ser como mínimo igual a U (figura 26-13).
Se utiliza una relación exponencial de tal manera que los
píxeles más profundos se corresponden con las luminancias
más bajas. Un resultado obtenido por esta técnica se mues-
Imagen I Erosionada E(I) I – E(I) tra en la figura 26-14.
Figura 26-12. Determinación de los contornos por un método
La modelización de los contornos es más compleja pero
morfológico. La imagen erosionada, cuyos objetos se han ofrece numerosas posibilidades, por ejemplo, la representa-
empequeñecido, se sustrae de la imagen original. ción simultánea de varios órganos con efecto de transparen-
cia. Esta representación comporta tres etapas:
– la segmentación de las imágenes tiene como objeti-
vo localizar zonas de las imágenes que corresponden a las
al número total de dilataciones. Dos ejemplos en particu-
lar son el opening (apertura, erosión seguida de dilatación)
que elimina los «pelillos», y el closing (cierre, dilatación se- Corte n
guida de una erosión), que cierra los huecos (figura 26-10). S
Dirección de
La apertura y el cierre no pueden superponerse a sí mismos: la observación
por ejemplo, después de una primera apertura, una segunda
no tiene efecto. Un ejemplo de filtro morfológico más com-
plejo (1 dilatación, 3 erosiones, 2 dilataciones) denominada
P
(DEEEDD) se muestra en la figura 26-11. Un filtro de estas
3D
características permite eliminar completamente las barbas
imagen
y los huecos. n de la
línea
Las transformaciones morfológicas tienen múltiples
aplicaciones en el tratamiento de imágenes. Citemos, por Figura 26-13. Técnica del nivel de profundidad.
360 IMÁGENES
contorno spline
corte n
Síntesis de la información
de una serie temporal de imágenes
Las series temporales de imágenes, por ejemplo gam-
magrafías de un mismo órgano obtenidas a diferentes tiem-
pos, pueden visualizarse por simple yuxtaposición de las
mismas. Este modo de representación es generalmente di-
Figura 26-17. Imagen ósea tridimensional obtenida por la técnica fícil de interpretar; sin embargo existen numerosas técnicas
de los marching cubes. (Negativo cedido generosamente por el Doctor que ponen en evidencia detalles poco (o nada) aparentes en
Zouaoui, servicio de Neurorradiología, Hospital de La Salpêtrière). la serie primitiva. Daremos ejemplos de análisis de una serie
gammagráfica cardíaca (figura 26-18) obtenida por adqui-
valor diagnóstico está por demostrar. No permite ver, evi- sición sincronizada con un ECG (véase capítulo 25, Cons-
dentemente, los colores de las paredes exploradas ni hacer titución de las imágenes).
biopsias de las zonas sospechosas, y la irradiación producida
por el escáner no es despreciable.
Otras aplicaciones del procesamiento de imagen tridi- El tratamiento paramétrico de imágenes
mensional son, por ejemplo, la simulación de intervención
quirúrgica, la de los movimientos óseos (cadera) o la visuali- El tratamiento paramétrico de imágenes consiste en
zación de superficies «isodosis» en radioterapia (véase capí- construir una o más imágenes artificiales, en las que cada
tulo 28, Protocolos de tratamiento en radioterapia). píxel representa un parámetro calculado a partir de valores
de píxeles homólogos de la serie inicial. Consideremos por
ejemplo, una serie P de imágenes: I1, I2....IP. Un mismo píxel
Procesamiento multimodal de imágenes (l,c) tiene para esas P imágenes los valores I1(l,c), I2(l,c)....
IP(l,c). Una imagen paramétrica M puede, por ejemplo, cons-
El establecer una correspondencia entre imágenes saca- truirse tomando por valor en (l,c) el máximo de los valores
das de modalidades diferentes tiene como objetivo el reunir del píxel (l,c) de la serie inicial. La imagen M (figura 26-19:
en una misma imagen 2D o 3D, datos morfológicos de alta M) viene definida por:
definición (IRM o escáner) y datos funcionales (PET o tomo– M(l,c) = Supt=1…P(It(l,c))
gammagrafía). El tratamiento multimodal de imágenes hace
uso de referencias anatómicas o artificiales presentes en los Otro ejemplo es la imagen paramétrica R cuyo valor del
dos tipos de imágenes que se quiere hacer corresponder. La píxel (l,c) es el rango t de la imagen para el que el valor It(l,c)
imagen de la modalidad menos resolutiva desde el punto de es máximo (figura 26-19: R):
vista morfológico se somete a transformaciones geométricas R(l,c) = t
(translación, rotación, deformación regular) que permita su-
tal que:
perponer los puntos de referencia homólogos. Las imágenes
se representan entonces simultáneamente utilizando paletas (∀k; It(l,c) Ik(l,c)).
Figura 26-18. Gammagrafía cardíaca obtenida por adquisición sincronizada con el ECG. La serie completa está compuesta de 16 imágenes.
362 IMÁGENES
Con (l,c), se lleva a cabo una imagen de fase, que represen-
ta el retardo temporal de la cinética de los diferentes píxeles
(figura 26-20).
Por tanto, este método identifica para cada píxel el mo-
delo a priori (ecuación 26–9) y luego visualiza los valores ob-
tenidos para cada parámetro del modelo en forma de ima-
gen. Esta aproximación puede utilizarse para modelos muy
diversos, por ejemplo una disminución exponencial: It(l,c) =
= a(l,c) exp(–λ(l,c).t). En este último caso, se obtienen dos
M: máximo/píxel R: tiempo del máximo/píxel imágenes paramétricas: una con a(l,c) que representa la tasa
de fijación inicial del trazador y una con λ(l,c) que representa
Figura 26-19. Imágenes paramétricas construidas a partir de la la velocidad a la que se elimina.
serie 26-18 (véase el texto). Para la imagen R, se han representado
únicamente los píxeles cuya suma sobrepasa un determinado umbral.
Estimación cuantitativa
Análisis temporal de Fourier de parámetros
El análisis temporal de Fourier es un caso particular de La cuantificación de la información recogida en una
imagen paramétrica muy utilizado. Consiste en suponer que imagen o en una serie de imágenes se ha convertido en un
para cada píxel (l,c), los valores I1(l,c), I2(l,c)....IP(l,c) siguen objetivo esencial en el análisis de las imágenes médicas. Se
una evolución sinusoidal y que se pueden por tanto encon- trata de obtener índices cuantitativos objetivos que permitan
trar tres constantes m(l,c), a(l,c) y (l,c) que cumplan una clasificación nosológica más precisa y reproducible que
⎛ 2πt ⎞ la obtenida con una interpretación visual subjetiva, y así po-
It (l,c) m(l,c) + a(l,c) cos⎜⎜ + ϕ(l,c)⎟⎟⎟ (26-9) der comparar estos índices en los diferentes estadios de una
⎜⎝ P ⎠⎟
enfermedad para estimar su evolución o la eficacia de una
En esta ecuación, a(l,c) es denominada amplitud y terapia. Es imposible presentar las numerosísimas técnicas
(l,c), fase. Se trata pues de una modelización (a nivel de utilizadas (y por tanto, propuestas) de las que sólo daremos
cada píxel), de la evolución de la actividad registrada. El cál- dos ejemplos.
culo de los coeficientes m(l,c), a(l,c) y (l,c) para cada píxel
se produce determinando los coeficientes de Fourier (véase
anexo 8), cuyo cálculo es muy rápido. Con los valores a(l,c) se Morfometría
genera una imagen de amplitud, que representa para cada
La morfometría consiste en efectuar medidas sobre las
píxel la amplitud de las variaciones de la actividad registrada.
imágenes. Pueden ser medidas de longitud: habitualmente
manuales. Se realizan gracias a «utilidades» o «rutinas» de
los sistemas de procesamiento de imágenes. Este tipo de
AI medida es ampliamente utilizado en ecografía obstétrica,
AD por ejemplo. De la misma manera también se pueden ob-
tener los ángulos. Son muy utilizadas en ortopedia sobre
radiografías óseas simples, así como de áreas, perímetros o
volúmenes.
VI En estos últimos tres casos, es necesario trazar en la ima-
VD
gen el contorno de la región que se quiere medir, por ejem-
plo, el área; esta determinación se debe hacer necesariamen-
te en todas las imágenes de una serie para la determinación
Imagen de amplitud Imagen de fase de un volumen. El trazado del contorno se hace con la ayuda
Figura 26-20. Imágenes de Fourier obtenidas con la serie de la
de líneas continuas: polígonos o splines (véase anexo 11). El
figura 26-18. La imagen de amplitud es tanto más intensa cuanto área se obtiene contando los píxeles que se encuentran en el
más variable es la actividad de los píxeles, por tanto en los bordes interior del contorno o por cálculo directo del área delimita-
de los ventrículos (VD y VI) y de las aurículas (AD y AI). La imagen da por el spline. El perímetro de una región se estima por el
de fase tiene el mismo nivel de gris para todos los puntos cuyo ciclo cálculo del de los splines (figura 26-21).
está sincronizado. Muestra un retraso temporal en la actividad de
Las operaciones de morfometría sobre volumen requie-
una zona del ventrículo izquierdo (señalado por una flecha). Esta
zona corresponde a un bloqueo de la rama izquierda, totalmente
ren una intervención manual laboriosa, más o menos sub-
indetectable en la serie inicial. Normalmente, las imágenes de fase jetiva para trazar los contornos de las zonas examinadas en
utilizan una paleta de colores que permite una interpretación más fácil. una serie de cortes. El volumen se calcula a partir de las áreas
Figura 26-22. Gammagrafía dinámica de riñón trasplantado y dos RDI definidas en forma de spline y trazadas en una de las imágenes ampliada:
renal (R) y vesical (V). La serie completa se compone de 30 imágenes.
364 IMÁGENES
contenida en el ventrículo: hay sangre igualmente marcada 25:1), pero no permiten una reconstitución perfecta de la
en el miocardio, el tejido pulmonar, etc. Por tanto, los valores imagen inicial a partir de la imagen comprimida. Se produce
VTD y VTS sobrestiman el volumen de sangre intraventricular. una pérdida de información, en general imperceptible vi-
Se corrige esta sobrestimación utilizando una RDI acce- sualmente, pero que en el caso de documentación médica
soria, próxima al ventrículo izquierdo, denominada Z, y cuya puede plantear problemas diagnósticos y médico-legales.
actividad se supone que representa la actividad parásita que Esta pérdida de información tiene, además, efectos no con-
se sobrepone a la del ventrículo izquierdo. Se efectúa la co- trolables sobre el resultado de los procesamientos que pro-
rrección sustrayendo a VTD el valor ZTD que se obtiene en la ducen índices cuantitativos;
imagen al final de la diástole, teniendo en cuenta las áreas –los algoritmos no destructivos (p. ej., la codificación
respectivas AVI y AZ de las dos RDI utilizadas. Se obtiene, así, de Huffman o el algoritmo LZW) tienen una eficacia mucho
un valor corregido: menor (del orden de 2 a 3:1), pero realizan una restitución ad
integrum de la imagen inicial en el momento de la descom-
A VI presión. Además, estos algoritmos permiten la transmisión
WTD = VTD − Z TD
AZ de los valores numéricos de la imagen, independientemente
El mismo método permite corregir el valor telesistólico de la manera elegida para visualizarla en su destino, mien-
VTS en: tras que los algoritmos destructivos transmiten una imagen
cuyas características de visualización (paleta de colores, seg-
A VI mentación, función de correspondencia valor-nivel de gris,
WTS = VTS − Z TS
AZ zoom) vienen ya fijadas.
FE corregida =
WTD − WTS
(26-11)
Ejercicios
VI
WTD
Ejercicio 26-1. Las transformaciones siguientes, que hacen
Sencillas en su puesta a punto, al tiempo que robustas, corresponder la imagen J a la imagen I ¿son lineales? No se
las técnicas de RDI son ampliamente utilizadas, pero no hay tendrán en cuenta los problemas que puedan sobrevenir en
que subestimar sus inconvenientes. Su trazado es subjetivo los bordes de la imagen. Se explicará en cada caso lo que re-
y los fenómenos de superposición plantean al operador una presenta la transformación.
difícil elección: o trazar RDI grandes con el riesgo de mezclar
píxeles provenientes de estructuras diferentes, o trazar RDI a) J(l,c) = supi =−1,1(I(l + i,c + j)).
j=−1 ,1
pequeñas «puras», descartando una parte importante de la
l l
información disponible. b) J(l,c) = ∑ ∑ I(l + i,c + j).
i =−1 j=−1
L C
La transmisión de las imágenes médicas dentro de un d) J(l,c) = aI(l,c) + b; a, b R (discutir en función del
servicio de imagen (p. ej., entre el lugar de producción y el de valor de b).
interpretación), entre dos servicios de un mismo hospital (el
servicio que pide la prueba y el que la realiza) o a distancia Ejercicio 26-2. Sin considerar el problema de los bordes, cal-
(de uno a otro hospital o a la consulta de un especialista) no cular cuantas operaciones (sumas y multiplicaciones) son
plantea en todos los casos el mismo problema de entrega y necesarias para efectuar la convolución de una imagen N*N
confidencialidad. Si las dificultades en el primer caso son fá- por un núcleo (2K +1)*(2K +1) en el caso general. Aplicación
cilmente subsanables, con la instalación y uso de una línea numérica: N = 1024; K = 3.
especializada, en los otros dos casos el volumen de datos se
convierte en un factor limitante (ciertos exámenes represen- Ejercicio 26-3. ¿Cuál es el efecto de los núcleos de convolu-
tan varios centenares de megabytes) y la confidencialidad de ción siguientes?
la transmisión es esencial. El volumen considerable de los ⎡ 0 0 0 ⎤ ⎡1 0 0 ⎤ ⎡ 2 1 0 ⎤
datos plantea igualmente problemas difíciles (y costosos) ⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎢0 1 0⎥ , ⎢0 0 0⎥ , ⎢1 0 −1⎥ ,
para el archivo de las imágenes digitales. ⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎢0 0 0⎥ ⎢0 0 0⎥ ⎢0 −1 −2⎥
La compresión de las imágenes es, por tanto, una nece- ⎢⎣ ⎥⎦ ⎢⎣ ⎥⎦ ⎢⎣ ⎥⎦
sidad en la actualidad. Se pueden utilizar dos tipos de algo- ⎡1 1 1⎤ ⎡−1 −1 −1⎤
ritmos de compresión: 1 ⎢⎢ ⎥ 1⎢ ⎥
⎢ 1 5 1⎥⎥ , ⎢⎢−1 12 −1⎥⎥
–los algoritmos destructivos (p. ej., JPEG, muy utilizado 13 ⎢ ⎥ 4 ⎢−1 −1 −1⎥
en internet) tienen tasas de compresión elevadas (de hasta ⎢⎣1 1 1⎥⎦ ⎢⎣ ⎥⎦
366 IMÁGENES