Sumario
NOTA TCNICA
Revista de Fsica Mdica 2001; 2(2): 90-97
Codificacin por sensibilidad SENSE1, Velocidad sin lmites en I.R.M.
H. Rademaker*, C. Santa Marta**
*Philips Ibrica S.A. Divisin Sistemas Mdicos **Profesora Asociada, UNED Departamento de Fsica Matemtica y Fluidos
Introduccin La imagen por resonancia magntica (IRM) combina tres ventajas fundamentales, la ausencia de radiacin ionizante, la gran variedad de contrastes entre los tejidos, y la posibilidad de obtener cortes tomogrficos en cualquier direccin del espacio. Por todo ello la IRM se ha convertido en muchos casos durante los ltimos aos, en la modalidad de diagnstico por excelencia. Los estudios por IRM ofrecen una gran variedad de informacin diagnstica valiosa; no obstante algunas aplicaciones presentes y muchas potenciales presentan dificultades por los largos tiempos de adquisicin, inherentes a la tcnica de RM. Por ejemplo en estudios cardacos o vasculares es necesario mantener apneas de unos 20 segundos, lo que resulta complicado para algunos pacientes con patologas cardiovasculares. La razn por la que los tiempos de adquisicin son largos, es que la obtencin de datos se realiza de manera secuencial, es decir, para formar la imagen es necesario adquirir varios ecos (generalmente 256) consecutivamente. El tiempo em-
Fig.1. Equipo moderno de RM Philips INTERA.
pleado en recolectar cada eco vara aproximadamente entre los 10 y los 3.000 ms, dependiendo de la secuencia. La resolucin espacial de la imagen de RM est limitada por la intensidad de los gradientes magnticos empleados descrito por Paul C. Lauterbur in 19732. La adquisicin de secuencias complejas requiere la aplicacin de un gran nmero de gradientes con diferentes intensidades o codificaciones de fase para poder reconstruir imgenes de buena calidad. La rapidez de adquisicin de la de imagen depende en gran manera de la velocidad a la que el hardware es capaz tanto de alcanzar la intensidad de trabajo del gradiente (tiempo de subida) como de conmutarlo. Los gradientes se crean mediante unas bobinas colocadas dentro del imn que permiten manipular el campo magntico estable. La intensidad de gradiente determina caractersticas geomtricas de la imagen como el grosor de corte, el tamao de la imagen o campo de visin y la resolucin espacial o nmero de pxeles y tambin limita el mnimo tiempo de eco (TE) que se puede utilizar en una determinada secuencia. Los ltimos aos se han convertido en una verdadera carrera entre los distintos fabricantes para desarrollar sistemas capaces de generar gradientes de gran intensidad y muy rpidos, cuyo tiempo de subida mnimo en este momento se han alcanzado intensidades de hasta 60 mT/m y tiempos de subida de 400 mT/m/ms SR3. Este continuo desarrollo se encuentra cada vez con ms problemas de implementacin. Por un lado, por los altos costos de dichos sistemas de gradientes debido a las complejas soluciones tecnolgicas necesarias para mejorar las especificaciones. Y por otro lado los gradientes comercialmente disponibles son ya tan potentes que pueden llegar a afectar a la electrofisiologa del enfermo debido a los rpidos tiempos de subida, que aumentan el riesgo de estimulacin nerviosa perifrica4. Por esta razn, la Organizacin de la Administracin Estadounidense de Alimentacin y Frmacos (USFDA, United States
Codificacin por sensibilidad SENSE, Velocidad sin lmites en I.R.M.
91
Food and Drug Administration5) impone un lmite en la tasa de los campos magnticos variables (dB/dt) inducidos en el pacientes. Otro lmite impuesto por la USFDA est relacionado con el ruido acstico que producen las bobinas de gradiente6.
Codificacin en sensibilidad (SENSE) Todos estos problemas ha llevado a Philips Sistemas Mdicos junto con sus socios tecnolgicos a desarrollar alternativas prcticas que permitan mejorar los tiempos de adquisicin sin comprometer ni el contraste ni la resolucin espacial de las imgenes y sin tener que emplear valores de gradiente cercanos a los lmites de seguridad recomendados por la USFDA. Fruto de esta investigacin es una nueva tecnologa revolucionaria llamado SENSE ("SENSitivity Encoding") [Pruessmann, 1999 #9]. El principio de SENSE consiste en utilizar las diferencias en sensibilidad de los distintos elementos (Fig. 2) que conforman una antena de recepcin multi elemento como mtodo de codificacin espacial, en lugar de la codificacin por los gradientes. La ventaja clave de este mtodo, sobre la codificacin por gradientes, es que el mecanismo de sensibilidad permite la codificacin simultnea de todo el FOV por los distintos segmentos que forman las bobinas multielemento. De este modo se ahorra tiempo de adquisicin por reemplazar codificaciones secuenciales de gradientes o de Fourier por codificaciones en sensibilidad. Para generar imgenes "convencionales" de RM se emplean dos tipos de codificaciones (ambas por gradien-
tes): la codificacin en fase y la codificacin en frecuencia. La codificacin en frecuencia se aplica al mismo tiempo que se lee el eco, y no supone un aumento del tiempo de adquisicin. La codificacin en fase se aplica antes de la lectura del eco, cada eco tiene un valor distinto del gradiente de codificacin de fase (Fig. 3a). El tiempo de adquisicin de la imagen es directamente proporcional al nmero de codificaciones de fase que se empleen. SENSE propone sustituir un cierto nmero de codificaciones de fase por codificaciones en sensibilidad, lo que supone aadir un tercer tipo de codificacin a los datos crudos. Para aplicar la codificacin en sensibilidad es necesario utilizar una antena multielemento o "phase array". La gran ventaja de este mtodo es que todos los elementos de la antena estn adquiriendo todo el FOV al mismo tiempo y que la seal que emite un punto de la muestra es recogido por los diferentes elementos de la antena con diferente intensidad (Fig. 3b), dependiendo del mapa de sensibilidad de cada elemento. De este modo se dice que se est codificando en sensibilidad. Aparece as un factor de aceleracin en la adquisicin de datos que es proporcional al nmero de elementos de la antena que estn trabajando. Para obtener una imagen a partir de datos crudos "convencionales" hay que deshacer las codificaciones en fase y en frecuencia mediante una transformada de Fourier. Si se utiliza una tercera codificacin, tambin es necesario deshacerla por los mtodos que se explicarn posteriormente. Debido a la independencia fsica del efecto de codificacin en sensibilidad de los mecanismos de contraste y obtencin de imagen, la tcnica SENSE no afecta al con-
Fig 2. Sensibilidad de las bobinas de superficie. (Se pierde seal alejndose de la superficie).
92
H. Rademaker, et al
Fig. 3a. En I.R.M. se muestrea el espacio k con una densidad homognea e igualmente espaciada. Cada lnea corresponde a una codificacin de fase.
traste y puede ser utilizado con todas las tcnicas de adquisicin de imagen actualmente en uso. La gran mayora de las tcnicas de adquisicin clnicamente relevantes utiliza un esquema de codificaciones de fase producidos por cambios en los gradientes (Fig. 3a). La codificacin en sensibilidad reduce la densidad y en consecuencia, el nmero de estas codificaciones. En el espacio k reduciendo la codificacin de fase de esta manera significa que se muestrea el mismo rea del espacio k con menos lneas y de manera ms espaciada (Fig. 3b). El factor de reduccin de las lneas de lectura se denomina factor R. En la reconstruccin convencional, este tipo de codificacin de fase resulta en una reduccin del campo de visin (FOV) en el mismo factor. Esto resulta
en una imagen con aliasing o doblamiento, donde partes de la imagen se superponen. La cantidad de veces en que se dobla la imagen depende del factor R y est limitado por el nmero de elementos de la bobina. En la adquisicin convencional con bobinas de un solo elemento, este efecto es irreversible, as que este tipo de muestreo no se puede usar. No obstante con SENSE, con bobinas multielemento, se permite reconstruir el campo de visin completo sin el artefacto presente. Con la reconstruccin SENSE, se reconstruye primero las imgenes de cada elemento de la bobina por separado usando la transformada de Fourier DFT tradicional. Cada imagen presenta los artefactos de doblamiento, pero cada elemento ver el artefacto de modo diferente, dependien-
Fig. 3b. En adquisicin SENSE en nmero de codificaciones de fase se reduce por un factor R (en este ejemplo R=2). Se incrementa el espacio entre lneas que resulta en una reduccin del FOV pero con artefacto del doblamiento.
Codificacin por sensibilidad SENSE, Velocidad sin lmites en I.R.M.
93
Fig. 4. Reversin del artefacto de doblamiento.
do de la situacin de los elementos sobre la zona a estudiar. El segundo paso es reconstruir el campo de visin completo usando la combinacin de imgenes, su diferente sensibilidad de recepcin y lgebra (Fig. 4).
un punto p2 en (x2,y2) que est fuera del FOV seleccionado pero que aparece en la imagen debido al aliasing. Se tiene por tanto en la imagen A una intensidad SA(x1, y1) = d1A S(p1) + d2A S(p2)
Adquisicin SENSE Para empezar se adquiere la imagen con un FOV reducido, menor que el tamao de la muestra, esto implica utilizar menos codificaciones de fase y, por lo tanto, reducir el tiempo de adquisicin. En condiciones normales (uso de una sola antena) esto da lugar a una imagen con aliasing, es decir, las zonas de la muestra que no han sido incluidas en el FOV aparecen en la imagen en una localizacin que no las corresponde. En estas condiciones, este artefacto de aliasing es irreversible. Si la imagen se adquiere con dos antenas la situacin es diferente: la imagen con aliasing captada por cada una de las antenas presenta diferentes intensidades en cada uno de los pxeles. Es decir, en la imagen con aliasing, un pixel situado en la posicin (x1,y1) tiene diferente intensidad en la imagen adquirida por la antena A (imagen A) que en la imagen adquirida por la antena B (imagen B). Se sabe que la intensidad del pxel (x1,y1) es la suma de las contribuciones del punto de la muestra p1 situado realmente en la posicin (x1,y1) ms la contribucin de
y en la imagen B una intensidad SB(x1, y1) = d1B S(p1) + d2B S(p2) Las intensidades SA y SB son conocidas, se miden sobre las imgenes A y B, y d1A, d2A, d1B y d2B dependen de la distancia de los puntos p1 y p2 a cada una de las antenas, vienen determinados por el mapa de sensibilidad de cada antena (tambin conocidos). S(p1) y S(p2) son las intensidades que tendran los puntos p1 y p2 si la imagen se hubiera obtenido slo con codificaciones de Fourier, sus valores se obtienen fcilmente resolviendo el sistema de dos ecuaciones con dos incgnitas planteado unas lneas ms arriba. Esta es la manera de deshacer la codificacin en sensibilidad.
Obtencin de mapas de sensibilidad Para la reconstruccin SENSE de imgenes sin artefactos residuales de doblamiento, es preciso conocer la
94
H. Rademaker, et al
Fig. 5. Determinacin de mapas de sensibilidad. Divisin de la imagen del primer elemento (a) por la imagen obtenido con la bobina de cuerpo (b) se muestra en (c). Se identifican las zonas que slo contienen ruido mediante un umbral (d), el cual es filtrado de la imagen (e). Se extrapola la zona mediante un crecimiento de regiones (f) que es filtrado para dar un mapa de sensibilidad refinada.
sensibilidad de los diferentes elementos con exactitud. Una dificultad aadida es que la sensibilidad, en cierta medida depende tambin de la secuencia de adquisicin a emplear, incluyendo el tamao y forma del paciente. Por eso, el conocimiento a priori de las distintas sensibilidades no suele ser lo suficientemente preciso. Para cada elemento de la bobina se obtiene un mapa de sensibilidad. Para eso se divide una imagen de referencia de baja resolucin con campo de visin completo por una hecha con la bobina de cuerpo integrada. Cualquier error en estos mapas se propagar hasta la imagen final utilizando reconstruccin SENSE. Por ello los mapas de sensibilidad se ajustan para eliminar ruido, incluso se realizar una interpolacin para asegurar una cobertura completa incluso con tejidos en movimiento (Fig. 5). El tiempo necesario para la obtencin de los mapas de sensibilidad es escaso por la baja resolucin de las referencias y por que los mapas sirven para todas las reconstrucciones SENSE aplicadas al mismo volumen de inters.
del tiempo de adquisicin. Lo mismo ocurre al reducir los tiempos de adquisicin con gradientes y ancho bandas ms grandes. Aunque recientes estudios demuestran que con la RSR el ltimo es significantemente mayor, sin mencionar el aumento del ruido acstico y costo que conlleva. No obstante en SENSE, existe el riesgo de resaltar el ruido en el caso que las relaciones geomtricas entre los distintos mapas de sensibilidad, los distintos elementos no son ideales. Este efecto especfico se describe de manera cuantitativo por el factor de geometra local g. La RSR de una imagen SENSE se define entonces como:
Relacin seal ruido (RSR) Para codificacin por sensibilidad, igual que con todas las secuencias rpidas, la relacin seal ruido es una preocupacin clave. El anlisis terico muestra que la RSR en imgenes SENSE es proporcional a la raz cuadrada
donde RSRconv es la obtenida con la misma bobina pero sin utilizar SENSE. El factor geomtrico g depende de la sensibilidad de los elementos, del factor R y vara dentro del campo de visin y por lo menos es uno, y tiende a aumentar con el factor R. Tericamente, la tcnica SENSE funciona hasta factores igual al nmero de elementos de la bobina. No obstante, el realce del ruido relacionado con el factor geomtrico g, limita el rango de los factores de reduccin hasta un mximo de tres o cuatro con bobinas con 6 elementos.
Codificacin por sensibilidad SENSE, Velocidad sin lmites en I.R.M.
95
Aplicaciones La Codificacion por sensibilidad permite mejorar en varias maneras el protocolo de imagen de un estudio de I.R.M. Adems de reducir el tiempo de adquisicin, es posible obtener ms resolucin espacial junto con una mejor cobertura anatmica. En algunos casos incluso es posible mejorar la RSR. La necesidad de aplicar SENSE en ciertas secuencias viene dada cuando el sistema haya agotado todas las posibilidades y trucos de la codificacin por gradientes, y cuando aun se quiere ir ms deprisa y algo se permitira sacrificar la RSR. Esto es especialmente cierto en las secuencias 3D. La RSR en secuencias 3D es en general alta, y en el caso de estudios vasculares con contraste aun ms, por la contribucin de la seal del gadolinio. Por otro lado, el tiempo total de las secuencias 3D suele ser largo, convirtindolas en candidatos idneos para la tecnologa SENSE.
Reduccin del tiempo de adquisicin Existen mltiples aplicaciones donde una reduccin del tiempo de adquisicin puede ser muy beneficiosa. En especial aquellas secuencias donde se requiere la colaboracin del enfermo para contener la respiracin durante la adquisicin. En la angiografa con contraste, el tiempo est tambin limitado por el proceso a estudiar. En general es preciso obtener la informacin en la fase arterial, que tiene una duracin limitada y es slo visible en el primer pase del contraste. Para evitar la contaminacin venosa la adquisicin se debe limitar a tan slo 15 - 20 segundos. SENSE permite estudiar las distintas fases sin la necesidad de calcular el tiempo exacto para cada una de ellas. Simplemente se dispara la secuencia dinmica con una resolucin temporal de 3 - 4 s. entre una y otra. A posteriori slo hay que identificar las distintas fases (Fig. 7). Otra aplicacin donde SENSE incrementa la resolucin temporal es en las secuencias de perfusin cerebral e imagen funcional, donde se podr obtener, en combina-
Fig. 6. Los experimentos con maniques muestran imgenes SENSE adquiridas con una bobina multielemento de 5 elementos. Dos elementos en la parte superior, y tres en la parte inferior. Se ha utilizado diferentes factores R y direcciones de fase. Imgenes SENSE del maniqu utilizando diferentes factores R. Izq: suma de las imgenes convencionales. Centro: Reconstruccin SENSE. Der: Mapas del ruido relativo segn la teora SENSE.
Est claro que el nmero y las especificaciones de los canales individuales junto con el diseo de las bobinas juegan un papel importante en SENSE.
Fig. 7. 4D Angiografa con contraste, 4 s. Por adquisicin 3D.
96
H. Rademaker, et al
adems se admite una ligera reduccin en la RSR, se convierte en tan slo 13 ms/imagen (Fig. 8). Despus de tantas secuencias innovadoras, casi se olvidan las secuencias que ms se aplican todos los das. stas tambin pueden beneficiarse de SENSE. Un buen ejemplo son las secuencias Espin-Eco o Vasculares cerebrales que de repente duran en lugar de 8 minutos, la mitad (Fig. 9).
Fig. 8. B-FFE con sense.
Aumentar la resolucin espacial No siempre es necesario aumentar la velocidad de adquisicin. En estos casos se puede aplicar SENSE para duplicar la resolucin espacial de la secuencia (Fig. 10). Otra forma de mejorar la resolucin espacial, especialmente en secuencias de difusin10 cerebral, es bajar la sensibilidad magntica en las secuencias ecoplaner, reduciendo el factor EPI (ETL) a la mitad. Este efecto reduce visiblemente los artefactos de inhomogeneidad en las zonas de los senos cabernosos y odos, mejorando as la sensibilidad a las pequeas lesiones en dichas zonas (Fig. 11).
Fig. 9. T.O.F. 3D 512, 100 cortes, 3.22 min.
cin con tcnicas 3D como PRESTO 7 ("Principle of Echo Shifting through a train of observation") hasta 32 imgenes o cortes en tan slo un segundo, mejorando as considerablemente la sensibilidad del proceso. Sin lugar a duda, la cardio-RM es la aplicacin con ms potencial para SENSE8 debido al carcter dinmico de casi todos los estudios de corazn. Actualmente apneas de unos 20 s. es prctica habitual en las distintas secuencias. Por un lado, SENSE permitir bajar este tiempo de apnea y por otro lado, aumentar el nmero de imgenes (fases del corazn) por latido visualizando con mucha ms precisin el proceso de contraccin del msculo cardaco. Incluso las secuencias en tiempo real9 son ya una realidad. La prioridad de dichas secuencias es que no haya movimiento apreciable entre una y otra imagen. Este concepto tiene mucho inters, porque permite la visualizacin del movimiento cardaco sin la colaboracin del enfermo ni ningn tipo de sincronismo, como respiratorio o cardaco. Y ser especialmente til para monitorizar el enfermo durante una prueba de esfuerzo, donde para lograr una calidad suficiente no se debe superar los 30 ms por imagen. Con tcnicas convencionales este tiempo est en torno a los 70 ms/imagen. Con la tecnologa SENSE, y con una reduccin de un factor 3, la adquisicin se reduce a unos 27 ms, y si
Fig. 10a. 1 x 2 x 4 mm.
Fig. 10b. 1 x 2 x 2 mm.
Codificacin por sensibilidad SENSE, Velocidad sin lmites en I.R.M.
97
Referencias
1. Philips Medica mundi, Vol 44/2 Nov 2000, pag.10-16 K.P. Pruessmann, M.Weiger and P.Boesiger. 2. Lanterbur PC, Image Formation by Induced Local Interaction: Examples Employing Nuclear Magnetic Resonance, Nature (London) 1973: 242: 190 - 191. 3. Especificaciones PHILIPS GYROSCAN INTERA 1.5T, gradientes explorer con SYNTRASCAN (SENSE) 1999-2000, Rel. 7.x. 4. Bourland JD, Nyenhuis JA, Mouchawar GA, et al. Physiologic indicators of high MRI gradient-induced fields. In, Book of Abstracts Society of Magnetic Resonance in Medicine 1990; 1276. 5. Guidance update for dB/dt for the "Guidance for content and review of a magnetic resonance diagnostic device 510(k) application" document. Office of Device Evaluation/Center for Devices and Radiological Health. Robert A. Phillips, Computed Imaging Devices Branch. Rockville, MD, Oct. 11, 1995. 6. McJury M, Blug A, Joerger C, Condon B, Wyper D. Acoustic noise levels during magnetic resonance imaging scanning at 1.5 T. Br J Radiol 1994;67:413-415. 7. Lui G, Sobering G, Duyn J, Moonen CTW. A Functional MRI Technique Combining Principles of Echo-Shifting with a Train of Observations (PRESTO). Magn Reson Med 1993;30: 764-768. 8. Pruessmann KP, Weiger M, Boesiger P. Sensitivity Encoded Cardiac MRI. J Cardiovasc Magn Reson 2001;3:1-9. 9. Pruessmann KP, Weiger M, Boesiger P. Cardiac Real Time imaging using SENSE Magn. Resonance Med 2000;43: 177 - 184. 10. Pruessmann KP, Weiger M, Muiswinkel AMC, Boesiger P, Sensitivity Encoding for Single Shot Diffusion Imaging. In: Proc ISMRM, 7th Annual Meeting, Philadelphia, 1999: 1815.
Fig. 11. Reduccin de susceptibilidad en secuencias SSh Difusin.