Mapeo cuantitativo de la tensión de un aneurisma de la aorta abdominal

Biomedical Engineering

Your institution must subscribe to JoVE's Engineering collection to access this content.

Fill out the form below to receive a free trial or learn more about access:

 

Overview

Fuente: Hannah L. Cebull1, Arvin H. Soepriatna1, John J. Boyle2 y Craig J. Goergen1

1 Weldon School of Biomedical Engineering, Universidad Purdue, West Lafayette, Indiana

2 Ingeniería Mecánica y Ciencia de los Materiales, Universidad de Washington en St. Louis, St Louis, Missouri

El comportamiento mecánico de los tejidos blandos, como los vasos sanguíneos, la piel, los tendones y otros órganos, está fuertemente influenciado por su composición de elastina y colágeno, que proporcionan elasticidad y fuerza. La orientación de la fibra de estas proteínas depende del tipo de tejido blando y puede variar desde una sola dirección preferida hasta intrincadas redes malladas, que pueden alterarse en el tejido enfermo. Por lo tanto, los tejidos blandos a menudo se comportan anisotrópicamente a nivel celular y de órganos, creando una necesidad de caracterización tridimensional. Desarrollar un método para estimar de forma fiable los campos de deformación unitaria dentro de tejidos o estructuras biológicas complejas es importante para caracterizar y comprender mecánicamente la enfermedad. La tensión representa cómo el tejido blando se deforma relativamente con el tiempo, y se puede describir matemáticamente a través de varias estimaciones.

La adquisición de datos de imagen a lo largo del tiempo permite estimar la deformación y la tensión. Sin embargo, todas las modalidades de imágenes médicas contienen cierta cantidad de ruido, lo que aumenta la dificultad de estimar con precisión la tensión in vivo. La técnica descrita aquí supera con éxito estos problemas mediante el uso de un método de estimación de deformación directa (DDE) para calcular campos de deformación unitaria 3D que varían espacialmente a partir de datos de imagen volumétrica.

Los métodos actuales de estimación de deformación unitaria incluyen correlación de imágenes digitales (DIC) y correlación de volumen digital. Desafortunadamente, DIC sólo puede estimar con precisión la tensión de un plano 2D, limitando severamente la aplicación de este método. Si bien son útiles, los métodos 2D como DIC tienen dificultadpara cuantificar la tensión en regiones que se someten a deformación 3D. Esto se debe a que el movimiento fuera del plano crea errores de deformación. La correlación de volumen digital es un método más aplicable que divide los datos de volumen inicialen regiones y encuentra la región más similar del volumen deformado, lo que reduce el error fuera del plano. Sin embargo, este método demuestra ser sensible al ruido y requiere suposiciones sobre las propiedades mecánicas del material.

La técnica que se muestra aquí elimina estos problemas mediante el uso de un método DDE, por lo que es muy útil en el análisis de datos de imágenes médicas. Además, es robusta a cepa alta o localizada. Aquí describimos la adquisición de datos de ultrasonido 4D cerrados y volumétricos, su conversión en un formato analizable y el uso de un código Matlab personalizado para estimar la deformación 3D y las cepas de Green-Lagrange correspondientes, un parámetro que describe mejor las grandes deformaciones. El tensor de cepa Green-Lagrange se implementa en muchos métodos de estimación de deformación unitaria 3D porque permite calcular F a partir de un ajuste de mínimos cuadrados (LSF) de los desplazamientos. La ecuación siguiente representa el tensor de cepa Green-Lagrange, E, donde F y I representan el gradiente de deformación y el tensor de identidad de segundo orden, respectivamente.

Equation 1(1)

Cite this Video

JoVE Science Education Database. Ingeniería Biomédica. Mapeo cuantitativo de la tensión de un aneurisma de la aorta abdominal. JoVE, Cambridge, MA, (2020).

Principles

El ultrasonido 4D es un volumen dinámico que se adquiere utilizando un motor de traducción lineal unido a un transductor de ultrasonido, lo que permite la adquisición de bucles de vídeo secuenciales cerrados cardíacos y respiratorios a través de una región de interés. Este método es útil para visualizar estructuras complejas como el corazón, donde la hipertrofia o infarto causa geometrías únicas, o aneurismas aórticos, donde la expansión y disección de los vasos asimétricos a menudo ocurren en vasos tortuosos. Además, los datos 4D pueden proporcionar información espacial y temporal de alta resolución, que también es importante para las imágenes cardiovasculares.

El método DDE aplicado a los datos de ultrasonido 4D es superior a otros métodos porque utiliza el registro de imágenes no rígidas. Los tensores de gradiente de deformación se estiman tradicionalmente a partir de campos de desplazamiento después de la correlación de volumen digital. Por el contrario, el método DDE estima intrínsecamente los tensores de gradiente de deformación durante el registro del volumen optimizando una función de deformación que se elige cuidadosamente para ser directamente análoga al tensor de deformación. La función de deformación depende tanto de la posición espacial como del parámetro de deformación (p):

Equation 2(2)

Los tres primeros elementos de esta función representan el tensor de gradiente de deformación, F, lo que permite que el cálculo de la deformación se incorpore directamente a la función de deformación. Este método de deformación se ha demostrado para aumentar la precisión y precisión de la estimación de la tensión en comparación con técnicas anteriores similares, ya que permite deformaciones grandes o localizadas que se encuentran comúnmente en los tejidos blandos.

Procedure

1. Configuración por ultrasonido 4D

  1. Cuando utilice el software de imágenes, utilice un portátil capaz de ejecutar software informático matemático para automatizar el proceso de adquisición 4D. Conecte el portátil con este código personalizado al sistema de ultrasonido a través del puerto USB. Tenga en cuenta que el software de imágenes tiene una función de ultrasonido 4D integrada en el software.
  2. Después de encender el sistema de ultrasonido, configure la unidad de monitoreo fisiológico mientras se asegura de que los botones de frecuencia cardíaca y temperatura estén encendidos. Inicializar la etapa del motor 3D conectada al soporte del transductor.
  3. Utilice el estadio y el transductor de ultrasonido adecuados para la adquisición de imágenes. Asegúrese de que se realizan todas las conexiones adecuadas.
  4. Proceda con la anestesia y preparación del animal para la toma de imágenes. Añadir pomada oftálmica a los ojos para evitar la desecación corneal, asegurar las patas a los electrodos del escenario, e insertar una sonda de temperatura rectal lubricada. Retire el pelaje en el área de interés usando una crema depilatoria.
  5. Asegúrese de que la crema depilatoria se haya completado. A continuación, aplicar una cantidad generosa de gel de transducción ultrasónica calentado al animal. Esto es especialmente importante para crear una buena conexión en toda la región de interés para las imágenes 4D.

2. Adquisición de ultrasonido sin ultrasonido 4D

  1. Comience un nuevo estudio sobre el sistema de ultrasonido y abra la ventana de imágenes en modo B (modo de brillo). Baje el transductor sobre el animal y localice la región de interés utilizando las perillas de los ejes X e Y en el escenario, asegurándose de que la frecuencia respiratoria no disminuya sustancialmente. Monitoree esto en la parte inferior de la pantalla.
  2. Coloque el transductor dentro del centro de la región de interés. A partir de ahí, se aproxima la distancia necesaria para que el transductor se mueva hacia arriba y hacia abajo de tal forma que se incluya toda la región de interés.
  3. Introduzca las dimensiones aproximadas en el código de software informático, incluido un tamaño de paso que suele ser de 0,08 mm para las imágenes de aneurisma de la aorta abdominal. Comience a ejecutar el código después de asegurarse de que el corazón del animal y las tasas respiratorias estén estables. Esto es importante para reducir los errores al reconstruir imágenes.
  4. Después de completar la adquisición de imágenes, exporte los datos como archivos XML sin procesar.

3. Conversión de datos de ultrasonido 4D

  1. Introduzca los archivos XML sin procesar en software que pueda convertir los datos en el formato adecuado para el análisis de manchas 3D. Aquí usamos Matlab para convertir archivos XML en archivos MAT. El script completo de Matlab está disponible aquí.
  2. Para una conversión adecuada, también será necesario introducir el número de fotogramas, el tamaño del paso y la resolución de salida deseada.
  3. Después de volver a muestrear la matriz a través del plano, importe el nuevo archivo MAT en el código de análisis de deformación unitaria 3D.

4. Análisis de código de tensión 3D

  1. Comience el análisis ajustando correctamente el archivo MAT importado. Por ejemplo, es posible que sea necesario cambiar la escala del volumen de la imagen para reducir el tiempo de cálculo.
  2. Introduzca la región que se va a analizar y determine la plantilla de malla adecuada para segmentar los datos de imagen como cuadros simples o polígonos elegidos manualmente. Es posible que sea necesario modificar el tamaño del cuadro de las regiones y el espaciado entre los puntos centrales para cada conjunto de datos. Los números óptimos elegidos para el tamaño del cuadro estarán alrededor del orden de píxeles de la entidad que se está rastreando, que se puede aproximar mirando el número de píxeles en dos dimensiones en un sector. El espaciado de las cajas determinará la resolución de los campos de deformación unitaria. Más cajas aumentarán la resolución, pero también pueden aumentar sustancialmente el tiempo de cálculo.
  3. Comience a calcular los jacobinos y gradientes iterativamente dentro de cada una de estas regiones. Una vez completado el precálculo, aplique la función de deformación.
  4. Calcule el tensor de gradiente de deformación. Primero calcule la deformación unitaria y, a continuación, calcule los valores propios y los autovectores utilizando el método de estimación de deformación directa.
  5. Trazar estos resultados en los planos deseados utilizando una técnica como la asignación de color de un plano de corte para representar el campo de deformación unitaria sobre la región de interés.

La toma de imágenes tridimensionales por tensión se utiliza para estimar la deformación de los tejidos blandos a lo largo del tiempo y comprender la enfermedad. El comportamiento mecánico de los tejidos blandos, como la piel, los vasos sanguíneos, los tendones y otros órganos, está fuertemente influenciado por su composición extracelular, que puede verse alterada por el envejecimiento y la enfermedad. Dentro de los tejidos biológicos complejos, es importante caracterizar estos cambios, que pueden afectar significativamente las propiedades mecánicas y funcionales de un órgano.

El mapeo cuantitativo de deformación unitaria utiliza datos de imagen volumétrica y un método de estimación de deformación directa para calcular los campos de deformación unitaria que varían espacialmente. Este video ilustrará los principios del mapeo de deformación unitaria, demostrará cómo se utiliza el mapeo cuantitativo de cepas para estimar los campos de tensión dentro de tejidos biológicos complejos y discutirá otras aplicaciones.

Los tejidos biológicos están fuertemente influenciados por la composición y orientación de la elastina y el colágeno. La elastina proteica es un componente altamente elástico de los tejidos que se estiran y contraen continuamente, como los vasos sanguíneos y los pulmones. El colágeno es la proteína más abundante en el cuerpo, y se ensambla a partir de polímeros triple-helicales individuales que se agrupan en fibras más grandes que proporcionan integridad estructural a los tejidos que van desde la piel hasta los huesos.

La orientación de estas proteínas abarca desde fibras alineadas hasta redes de malla fibrosa, que afecta a las propiedades mecánicas del tejido. La tensión es una medida de la deformación relativa de los tejidos blandos a lo largo del tiempo, y se puede utilizar para visualizar lesiones y enfermedades. Se describe y mapea mediante estimaciones matemáticas.

Para mapear la tensión en órganos complejos, como el corazón, se pueden utilizar datos de ultrasonido en cuatro dimensiones, que proporcionan información de alta resolución, espacial y temporal. A continuación, se aplica el método de estimación de deformación directa, o DDE, a los datos. Un código se utiliza para estimar la deformación 3D y las cepas de Green-Lagrange correspondientes utilizando la siguiente ecuación.

El tensor de cepa Green-Lagrange depende del tensor gradiente de deformación y del tensor de identidad del segundo orden. Los tensores de gradiente de deformación se estiman tradicionalmente a partir de campos de desplazamiento. En el método DDE, una función de deformación está optimizada para ser directamente análoga al tensor de deformación. La función de deformación depende tanto de la posición espacial como del parámetro de deformación. El cálculo de la deformación se incorpora directamente a la función de deformación. Los primeros nueve elementos representan el tensor de gradiente de deformación.

Este método se utiliza para estimar deformaciones grandes y localizadas en tejidos blandos. Ahora que entendemos los principios de la asignación de cepas, veamos ahora cómo se realiza el mapeo de deformación unitaria para detectar aneurismas aórticos en ratones.

Para comenzar la configuración, abra el software Vivo 2100 y conecte el portátil al sistema de ultrasonido. Asegúrese de que la unidad de monitoreo fisiológico esté encendida para medir la frecuencia cardíaca y la temperatura. A continuación, inicialice la etapa del motor 3D.

Instale el transductor de ultrasonido y asegúrese de que se realizan todas las conexiones adecuadas. A continuación, anestesia al animal que será imageado usando 3% de isoflurano en una cámara derribada. Una vez anestesiado el ratón, muévalo a la etapa calentada y asegure un cono nasal para entregar 1-2% de isoflurano. Aplique pomada oftálmica a los ojos y fije las patas a los electrodos del escenario para monitorear la respiración y la frecuencia cardíaca del animal. A continuación, inserte una sonda de temperatura rectal. Aplicar crema depilatoria para eliminar el cabello de la zona de interés, y luego aplicar una generosa cantidad de gel de ultrasonido caliente en el área depilada.

Para iniciar la adquisición de imágenes, primero, abra la ventana de imágenes y seleccione el modo B. A continuación, baje el transductor sobre el animal y utilice las perillas del eje X e Y en el escenario para localizar el área de interés. Controle la frecuencia respiratoria para asegurarse de que no disminuya sustancialmente. Coloque el transductor en el centro de la región de interés. A continuación, aproxime la distancia necesaria para cubrir toda la región de interés.

Introduzca estas dimensiones en el código MATLAB y elija un tamaño de paso de 0,08 milímetros. Asegúrese de que el corazón del animal y las tasas respiratorias estén estables y, a continuación, ejecute el código MATLAB.

Después de la adquisición de imágenes, exporte los datos como archivos XML sin procesar y conviértalos en archivos MAT. Asegúrese de introducir el número de fotogramas, el tamaño del paso y la resolución de salida. A continuación, vuelva a muestrear la matriz en el plano pasante.

Importe el nuevo archivo MAT al código de análisis de deformación unitaria 3D. Puede ser necesario cambiar la escala del archivo para reducir el tiempo de cálculo. A continuación, introduzca la región que desea analizar. Proxime el número de píxeles en un sector bidimensional de la entidad rastreada y seleccione la plantilla de malla como un cuadro simple o polígonos elegidos manualmente. Elija el número de píxel óptimo para el tamaño de malla. Calcular los jacobianos y los gradientes. Repita el proceso para cada región. A continuación, aplique la función de deformación.

A continuación, utilizando deformaciones cartesianas calculadas a partir de DDE, determine los valores propios y los autovectores de la deformación. A continuación, seleccione los sectores para los que desea trazar los valores de deformación unitaria desplazándose por las vistas de eje largo, eje de ordenación y eje coronal.

Pulse Seleccionar colector para análisis. A continuación, utilice el cursor para colocar marcadores a lo largo de la pared aórtica, incluyendo el trombo, el aneurisma y partes sanas de la aorta. Repita el proceso para todas las vistas. Por último, utilice la asignación de color para trazar los resultados del campo de deformación unitaria sobre la región de interés.

Echemos un vistazo de cerca al ejemplo de un aneurisma suprarrenal inducido por angiotensina II adquirido de un ratón. En primer lugar, se obtienen varios bucles de visualización de kilohercios de kilohercio saque-zqui-eje corto con un tamaño de paso determinado a lo largo de la aorta y se combinan para crear datos 4D.

Después de realizar el cálculo de deformación unitaria 3D utilizando una función de deformación optimizada, se obtiene la gráfica de visualización de sectores 3D de la aorta infrarrenal. El mapa de color de la tensión verde principal se superpone para resaltar regiones de tensión de pared aórtica heterogénea. Además, las vistas de eje largo y eje corto revelan variaciones espaciales heterogéneas en la deformación unitaria, especialmente cuando hay un trombo presente.

Las gráficas de deformación unitaria correspondientes muestran valores de deformación unitaria más altos en regiones sanas de la aorta en el eje largo, mientras que la región aneurisma muestra una disminución de la tensión en el eje corto.

La visualización cuantitativa precisa de la tensión mediante estimación de deformación directa es una herramienta útil utilizada en diversas aplicaciones biomédicas.

Por ejemplo, se puede cuantificar la tensión cardíaca. Durante el ciclo cardíaco, el miocardio sufre deformación 3D. La cuantificación de la tensión en tres dimensiones es fundamental para caracterizar de forma fiable la dinámica de este tejido a lo largo del tiempo. Esto es útil para rastrear la progresión de la enfermedad en modelos animales.

Otra aplicación está en la caracterización del tejido intestinal. La imagen in vivo de los intestinos es difícil debido a los efectos de las estructuras circundantes. Sin embargo, calcular la tensión a partir de imágenes de fibrosis intestinal podría ser particularmente útil para proporcionar detección temprana de áreas problemáticas que requieren intervención quirúrgica.

A una escala mucho menor, este método DDE también se aplica al nivel celular mediante el uso de técnicas de imagen de mayor resolución, como la microscopía confocal. Sirve, por ejemplo, en la caracterización de la matriz extracelular para entender cómo las células se comunican bajo cambios mecánicos.

Acaba de ver la introducción de JoVE a la visualización cuantitativa de cepas. Ahora debe entender cómo medir la tensión tridimensional en los tejidos biológicos y cómo se utiliza en la detección temprana de enfermedades. ¡Gracias por mirar!

Results

Mediante el procedimiento descrito anteriormente, se adquirió una ecografía 4D de un aneurisma suprarrenal (AAA) de una angiotensina II. Se adquirieron varios bucles de vídeo EKV de eje corto a lo largo de la aorta y se combinaron para crear datos 4D, como se muestra en la Figura 1. Estos datos se convirtieron en un archivo MAT mediante un código personalizado, que luego se analizó en un código de cálculo de deformación unitaria 3D mediante una función de deformación. Después de optimizar los parámetros del código para un conjunto de datos específico, se produjo una vista representativa de eje largo con los valores de deformación unitaria correspondientes, así como un trazado de visualización de sectores 3D con un mapa de color de deformación unitaria superpuesta (Figura 2). Esta técnica DDE y los datos de deformación unitaria resaltan las variaciones espaciales heterogéneas en la tensión, especialmente cuando hay un trombo presente. Estos resultados se pueden correlacionar con la estructura del recipiente para determinar la relación entre la deformación in vivo y la composición del aneurisma.

Equation 3
Figura 1: Los bucles de visualización de kilohercios cerrados por ECG (EKV) de la aorta se adquieren de ubicaciones iniciales y finales introducidas manualmente, siguiendo un tamaño de paso de 0,2 mm.

Figure 1
Figura 2: Datos de ultrasonido de alta frecuencia 4D de un aneurisma de la aorta abdominal disuado murino representado en la sístole (A) con los principales campos de tensión estimados y superpuestos (B) (Scalebar a 5 mm). Vistas de eje largo y corto que representan las regiones aneurisma y sana que corresponden a la tensión principal en un ciclo cardíaco (sistolia: t-0,4) (C, D). Estos datos muestran niveles de tensión relativamente altos en regiones sanas y valores de tensión reducidos dentro del aneurisma de disección.

Applications and Summary

La caracterización mecánica localizada in vivo es una parte importante de la comprensión del crecimiento y la remodelación de los tejidos biológicos. En comparación con los enfoques existentes, el procedimiento de cuantificación de deformación unitaria descrito aquí utiliza un método mejorado para calcular con precisión la tensión 3D mediante la deformación óptima de la imagen no deformada antes de la correlación cruzada. Este método no utiliza ninguna suposición material para determinar las deformaciones unitarias dentro de los volúmenes de tejido. Desafortunadamente, la estimación de la tensión es confiable sólo hasta un tamaño de núcleo de 15x15x15 vóxeles cuando se utilizan datos de ultrasonido, lo que sugiere que este enfoque DDE puede no detectar características sutiles dentro de un campo de tensión. A pesar de esta limitación, sigue siendo una herramienta importante para investigar las respuestas mecánicas, diagnosticar la patología y mejorar los modelos de enfermedades.

Muchas áreas de investigación más allá de los aneurismas aórticos pueden beneficiarse de esta herramienta de medición de tensión. La tensión cardíaca también se puede cuantificar fácilmente con este método. Debido a que el miocardio sufre deformación 3D durante el ciclo cardíaco, cuantificar la tensión en tres dimensiones es integral para caracterizar de forma fiable la dinámica de este tejido. Los datos fiables de las cepas son especialmente importantes cuando se realiza un seguimiento de la progresión de la enfermedad en modelos animales.

El análisis de deformación unitaria 3D también se puede aplicar a las imágenes por ultrasonido intestinal. La caracterización mecánica del tejido intestinal se lleva a cabo más comúnmente in vitro. Sin embargo, esto no siempre es una verdadera representación del comportamiento real de los intestinos in vivo debido a los efectos de las estructuras circundantes. Como ejemplo de traducción clínica de este enfoque, el cálculo de la tensión a partir de imágenes de fibrosis intestinal debido a una presión luminal anormal podría proporcionar detección temprana de áreas problemáticas que requieren intervención quirúrgica.

Más allá de las aplicaciones a mayor escala, este método también se puede aplicar al nivel celular mediante el uso de técnicas de imagen de mayor resolución, como la microscopía confocal. Caracterizar la matriz extracelular es importante para entender cómo se comunican las células. Se ha llevado a cabo mucha investigación sobre la caracterización bioquímica, pero entender cómo se puede llevar a cabo la comunicación a través de respuestas mecánicas requiere una comprensión de la deformación y la tensión. La tensión a granel no es beneficiosa porque no hay manera de determinar el origen del cambio de deformación. La aplicación de un enfoque DDE de alta resolución podría revelar directamente cómo la matriz extracelular responde a los cambios mecánicos.

Agradecimientos

Nos gustaría reconocer a John Boyle, Guy Genin y Stavros Thomopoulos por la contribución del código Matlab personalizado DDE capaz de estimar directamente la cepa Lagrange-Green.

1. Configuración por ultrasonido 4D

  1. Cuando utilice el software de imágenes, utilice un portátil capaz de ejecutar software informático matemático para automatizar el proceso de adquisición 4D. Conecte el portátil con este código personalizado al sistema de ultrasonido a través del puerto USB. Tenga en cuenta que el software de imágenes tiene una función de ultrasonido 4D integrada en el software.
  2. Después de encender el sistema de ultrasonido, configure la unidad de monitoreo fisiológico mientras se asegura de que los botones de frecuencia cardíaca y temperatura estén encendidos. Inicializar la etapa del motor 3D conectada al soporte del transductor.
  3. Utilice el estadio y el transductor de ultrasonido adecuados para la adquisición de imágenes. Asegúrese de que se realizan todas las conexiones adecuadas.
  4. Proceda con la anestesia y preparación del animal para la toma de imágenes. Añadir pomada oftálmica a los ojos para evitar la desecación corneal, asegurar las patas a los electrodos del escenario, e insertar una sonda de temperatura rectal lubricada. Retire el pelaje en el área de interés usando una crema depilatoria.
  5. Asegúrese de que la crema depilatoria se haya completado. A continuación, aplicar una cantidad generosa de gel de transducción ultrasónica calentado al animal. Esto es especialmente importante para crear una buena conexión en toda la región de interés para las imágenes 4D.

2. Adquisición de ultrasonido sin ultrasonido 4D

  1. Comience un nuevo estudio sobre el sistema de ultrasonido y abra la ventana de imágenes en modo B (modo de brillo). Baje el transductor sobre el animal y localice la región de interés utilizando las perillas de los ejes X e Y en el escenario, asegurándose de que la frecuencia respiratoria no disminuya sustancialmente. Monitoree esto en la parte inferior de la pantalla.
  2. Coloque el transductor dentro del centro de la región de interés. A partir de ahí, se aproxima la distancia necesaria para que el transductor se mueva hacia arriba y hacia abajo de tal forma que se incluya toda la región de interés.
  3. Introduzca las dimensiones aproximadas en el código de software informático, incluido un tamaño de paso que suele ser de 0,08 mm para las imágenes de aneurisma de la aorta abdominal. Comience a ejecutar el código después de asegurarse de que el corazón del animal y las tasas respiratorias estén estables. Esto es importante para reducir los errores al reconstruir imágenes.
  4. Después de completar la adquisición de imágenes, exporte los datos como archivos XML sin procesar.

3. Conversión de datos de ultrasonido 4D

  1. Introduzca los archivos XML sin procesar en software que pueda convertir los datos en el formato adecuado para el análisis de manchas 3D. Aquí usamos Matlab para convertir archivos XML en archivos MAT. El script completo de Matlab está disponible aquí.
  2. Para una conversión adecuada, también será necesario introducir el número de fotogramas, el tamaño del paso y la resolución de salida deseada.
  3. Después de volver a muestrear la matriz a través del plano, importe el nuevo archivo MAT en el código de análisis de deformación unitaria 3D.

4. Análisis de código de tensión 3D

  1. Comience el análisis ajustando correctamente el archivo MAT importado. Por ejemplo, es posible que sea necesario cambiar la escala del volumen de la imagen para reducir el tiempo de cálculo.
  2. Introduzca la región que se va a analizar y determine la plantilla de malla adecuada para segmentar los datos de imagen como cuadros simples o polígonos elegidos manualmente. Es posible que sea necesario modificar el tamaño del cuadro de las regiones y el espaciado entre los puntos centrales para cada conjunto de datos. Los números óptimos elegidos para el tamaño del cuadro estarán alrededor del orden de píxeles de la entidad que se está rastreando, que se puede aproximar mirando el número de píxeles en dos dimensiones en un sector. El espaciado de las cajas determinará la resolución de los campos de deformación unitaria. Más cajas aumentarán la resolución, pero también pueden aumentar sustancialmente el tiempo de cálculo.
  3. Comience a calcular los jacobinos y gradientes iterativamente dentro de cada una de estas regiones. Una vez completado el precálculo, aplique la función de deformación.
  4. Calcule el tensor de gradiente de deformación. Primero calcule la deformación unitaria y, a continuación, calcule los valores propios y los autovectores utilizando el método de estimación de deformación directa.
  5. Trazar estos resultados en los planos deseados utilizando una técnica como la asignación de color de un plano de corte para representar el campo de deformación unitaria sobre la región de interés.

La toma de imágenes tridimensionales por tensión se utiliza para estimar la deformación de los tejidos blandos a lo largo del tiempo y comprender la enfermedad. El comportamiento mecánico de los tejidos blandos, como la piel, los vasos sanguíneos, los tendones y otros órganos, está fuertemente influenciado por su composición extracelular, que puede verse alterada por el envejecimiento y la enfermedad. Dentro de los tejidos biológicos complejos, es importante caracterizar estos cambios, que pueden afectar significativamente las propiedades mecánicas y funcionales de un órgano.

El mapeo cuantitativo de deformación unitaria utiliza datos de imagen volumétrica y un método de estimación de deformación directa para calcular los campos de deformación unitaria que varían espacialmente. Este video ilustrará los principios del mapeo de deformación unitaria, demostrará cómo se utiliza el mapeo cuantitativo de cepas para estimar los campos de tensión dentro de tejidos biológicos complejos y discutirá otras aplicaciones.

Los tejidos biológicos están fuertemente influenciados por la composición y orientación de la elastina y el colágeno. La elastina proteica es un componente altamente elástico de los tejidos que se estiran y contraen continuamente, como los vasos sanguíneos y los pulmones. El colágeno es la proteína más abundante en el cuerpo, y se ensambla a partir de polímeros triple-helicales individuales que se agrupan en fibras más grandes que proporcionan integridad estructural a los tejidos que van desde la piel hasta los huesos.

La orientación de estas proteínas abarca desde fibras alineadas hasta redes de malla fibrosa, que afecta a las propiedades mecánicas del tejido. La tensión es una medida de la deformación relativa de los tejidos blandos a lo largo del tiempo, y se puede utilizar para visualizar lesiones y enfermedades. Se describe y mapea mediante estimaciones matemáticas.

Para mapear la tensión en órganos complejos, como el corazón, se pueden utilizar datos de ultrasonido en cuatro dimensiones, que proporcionan información de alta resolución, espacial y temporal. A continuación, se aplica el método de estimación de deformación directa, o DDE, a los datos. Un código se utiliza para estimar la deformación 3D y las cepas de Green-Lagrange correspondientes utilizando la siguiente ecuación.

El tensor de cepa Green-Lagrange depende del tensor gradiente de deformación y del tensor de identidad del segundo orden. Los tensores de gradiente de deformación se estiman tradicionalmente a partir de campos de desplazamiento. En el método DDE, una función de deformación está optimizada para ser directamente análoga al tensor de deformación. La función de deformación depende tanto de la posición espacial como del parámetro de deformación. El cálculo de la deformación se incorpora directamente a la función de deformación. Los primeros nueve elementos representan el tensor de gradiente de deformación.

Este método se utiliza para estimar deformaciones grandes y localizadas en tejidos blandos. Ahora que entendemos los principios de la asignación de cepas, veamos ahora cómo se realiza el mapeo de deformación unitaria para detectar aneurismas aórticos en ratones.

Para comenzar la configuración, abra el software Vivo 2100 y conecte el portátil al sistema de ultrasonido. Asegúrese de que la unidad de monitoreo fisiológico esté encendida para medir la frecuencia cardíaca y la temperatura. A continuación, inicialice la etapa del motor 3D.

Instale el transductor de ultrasonido y asegúrese de que se realizan todas las conexiones adecuadas. A continuación, anestesia al animal que será imageado usando 3% de isoflurano en una cámara derribada. Una vez anestesiado el ratón, muévalo a la etapa calentada y asegure un cono nasal para entregar 1-2% de isoflurano. Aplique pomada oftálmica a los ojos y fije las patas a los electrodos del escenario para monitorear la respiración y la frecuencia cardíaca del animal. A continuación, inserte una sonda de temperatura rectal. Aplicar crema depilatoria para eliminar el cabello de la zona de interés, y luego aplicar una generosa cantidad de gel de ultrasonido caliente en el área depilada.

Para iniciar la adquisición de imágenes, primero, abra la ventana de imágenes y seleccione el modo B. A continuación, baje el transductor sobre el animal y utilice las perillas del eje X e Y en el escenario para localizar el área de interés. Controle la frecuencia respiratoria para asegurarse de que no disminuya sustancialmente. Coloque el transductor en el centro de la región de interés. A continuación, aproxime la distancia necesaria para cubrir toda la región de interés.

Introduzca estas dimensiones en el código MATLAB y elija un tamaño de paso de 0,08 milímetros. Asegúrese de que el corazón del animal y las tasas respiratorias estén estables y, a continuación, ejecute el código MATLAB.

Después de la adquisición de imágenes, exporte los datos como archivos XML sin procesar y conviértalos en archivos MAT. Asegúrese de introducir el número de fotogramas, el tamaño del paso y la resolución de salida. A continuación, vuelva a muestrear la matriz en el plano pasante.

Importe el nuevo archivo MAT al código de análisis de deformación unitaria 3D. Puede ser necesario cambiar la escala del archivo para reducir el tiempo de cálculo. A continuación, introduzca la región que desea analizar. Proxime el número de píxeles en un sector bidimensional de la entidad rastreada y seleccione la plantilla de malla como un cuadro simple o polígonos elegidos manualmente. Elija el número de píxel óptimo para el tamaño de malla. Calcular los jacobianos y los gradientes. Repita el proceso para cada región. A continuación, aplique la función de deformación.

A continuación, utilizando deformaciones cartesianas calculadas a partir de DDE, determine los valores propios y los autovectores de la deformación. A continuación, seleccione los sectores para los que desea trazar los valores de deformación unitaria desplazándose por las vistas de eje largo, eje de ordenación y eje coronal.

Pulse Seleccionar colector para análisis. A continuación, utilice el cursor para colocar marcadores a lo largo de la pared aórtica, incluyendo el trombo, el aneurisma y partes sanas de la aorta. Repita el proceso para todas las vistas. Por último, utilice la asignación de color para trazar los resultados del campo de deformación unitaria sobre la región de interés.

Echemos un vistazo de cerca al ejemplo de un aneurisma suprarrenal inducido por angiotensina II adquirido de un ratón. En primer lugar, se obtienen varios bucles de visualización de kilohercios de kilohercio saque-zqui-eje corto con un tamaño de paso determinado a lo largo de la aorta y se combinan para crear datos 4D.

Después de realizar el cálculo de deformación unitaria 3D utilizando una función de deformación optimizada, se obtiene la gráfica de visualización de sectores 3D de la aorta infrarrenal. El mapa de color de la tensión verde principal se superpone para resaltar regiones de tensión de pared aórtica heterogénea. Además, las vistas de eje largo y eje corto revelan variaciones espaciales heterogéneas en la deformación unitaria, especialmente cuando hay un trombo presente.

Las gráficas de deformación unitaria correspondientes muestran valores de deformación unitaria más altos en regiones sanas de la aorta en el eje largo, mientras que la región aneurisma muestra una disminución de la tensión en el eje corto.

La visualización cuantitativa precisa de la tensión mediante estimación de deformación directa es una herramienta útil utilizada en diversas aplicaciones biomédicas.

Por ejemplo, se puede cuantificar la tensión cardíaca. Durante el ciclo cardíaco, el miocardio sufre deformación 3D. La cuantificación de la tensión en tres dimensiones es fundamental para caracterizar de forma fiable la dinámica de este tejido a lo largo del tiempo. Esto es útil para rastrear la progresión de la enfermedad en modelos animales.

Otra aplicación está en la caracterización del tejido intestinal. La imagen in vivo de los intestinos es difícil debido a los efectos de las estructuras circundantes. Sin embargo, calcular la tensión a partir de imágenes de fibrosis intestinal podría ser particularmente útil para proporcionar detección temprana de áreas problemáticas que requieren intervención quirúrgica.

A una escala mucho menor, este método DDE también se aplica al nivel celular mediante el uso de técnicas de imagen de mayor resolución, como la microscopía confocal. Sirve, por ejemplo, en la caracterización de la matriz extracelular para entender cómo las células se comunican bajo cambios mecánicos.

Acaba de ver la introducción de JoVE a la visualización cuantitativa de cepas. Ahora debe entender cómo medir la tensión tridimensional en los tejidos biológicos y cómo se utiliza en la detección temprana de enfermedades. ¡Gracias por mirar!

JoVE Science Education is free through June 15th 2020.

RECOMMEND JoVE