$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Se recolectaron imágenes longitudinales de μCT para este estudio en las semanas 0 y 6 de defectos radiales críticos de 3 mm en ratas hembras adultas Charles River SASCO-SD tratadas con un andamio a base de policaprolactona (PCL). Todo el uso de animales se realizó de acuerdo con los protocolos aprobados por el Comité de Recursos Animales (UCAR) de la Universidad de Rochester. La recolección de imágenes μCT se realizó con Scanco Medical VivaCT 40.
NOTA: Los pasos principales de este protocolo se dividen en segmentación de imágenes μCT, alineación, selección y recorte de ROI, y análisis y visualización (Figura 1). El protocolo para la segmentación de imágenes μCT es una adaptación de Kenney et al. (2022)11.

Figura 1: Diagrama de flujo de trabajo resumido. Los pasos del protocolo se dividen principalmente en segmentación de imágenes, alineación de modelos, selección de ROI y análisis de volumen. Haga clic aquí para ver una versión más grande de esta figura.
1. Segmentación de imágenes μCT
- Apertura de imágenes dentro del software
- Inicie el software Amira (en lo sucesivo, software de visualización) y seleccione Open Data. Dentro de la ventana emergente, navegue hasta el directorio de carpetas deseado para seleccionar el punto de tiempo inicial .dcm archivos. Presione CTRL+A para seleccionar todas las imágenes y haga clic en Abrir imágenes.
- Una ventana emergente mostrará información sobre el conjunto de datos de imagen abierto; haga clic en Aceptar para continuar.
- Aparecerá una advertencia sobre el escalado. Haga clic en Convertir en flotante (la opción predeterminada) para continuar.
- El icono del conjunto de datos aparecerá en la sección Vista del proyecto . Para cambiar el nombre del conjunto de datos, haga clic en el icono del conjunto de datos y presione F2; una vez que se ingresa el nombre deseado, haga clic en Aceptar. Establezca un patrón de nomenclatura para identificar el escaneo (por ejemplo, 'Material de andamio óseo' - 'WKX').
- Filtrado y umbral de la imagen μCT
- Con las imágenes cargadas en la sección Vista del proyecto , haga clic en el visor 2D predeterminado de Ortho Slice para ver las orientaciones y cortes 2D. Aparecerá una visualización 2D en la ventana de visualización de la derecha.
- Para visualizar las imágenes en 3D, haga clic con el botón derecho en el conjunto de datos, busque Representación de volumen y seleccione. Dentro de la sección Propiedades de la Representación de volúmenes, ajuste el umbral inferior del mapa de colores a 2500 y pulse Intro.
NOTA: 2D y 3D se pueden activar y desactivar haciendo clic en el cuadrado azul junto a Ortho Slice o Volume Rendering, respectivamente. Para la ventana de visualización de la derecha, hay varias opciones de visualización disponibles en la parte superior de la barra. Estos incluyen diferentes herramientas para mover, rotar o acercar/alejar la imagen; herramientas para cambiar la orientación de la imagen a las vistas XY, XZ e YZ, y herramientas para mediciones sencillas.
- Según sea necesario, algunos conjuntos de datos pueden requerir que se eliminen artefactos o datos fuera de la región de interés mediante una edición de volumen. Para agregar una edición de volumen, haga clic con el botón derecho en el conjunto de datos, busque Edición de volumen y seleccione.
- Haga clic en Editar volumen en la sección Vista del proyecto . Dentro de la sección Propiedades , seleccione TabBox en la primera ventana desplegable a la derecha de Herramienta.
- Al seleccionar las regiones que desea conservar, ajuste el cuadro de pestañas haciendo clic y arrastrando las esquinas verdes con el cursor de Interact (el icono del cursor apuntado a lo largo de la barra superior de la ventana de visualización). El ajuste de TabBox puede requerir el uso de vistas alternativas (XY, XZ o YZ) para identificar las regiones que se van a eliminar.
- Una vez que el TabBox rodea la región que se va a conservar, seleccione Cortar fuera para quitar el artefacto o los datos fuera del TabBox. Se generará un nuevo conjunto de datos modificado en la vista del proyecto.
- Haga clic y arrastre la cadena de representación de volumen del conjunto de datos original al conjunto de datos modificado. Esto cambiará la imagen que se muestra en la ventana de visualización al conjunto de datos modificado. Dentro de la sección Propiedades de la Representación de volúmenes, ajuste el umbral inferior del mapa de colores a 2500 y pulse Intro.
- Haga clic con el botón derecho en los datos modificados, busque Filtro de mediana y seleccione. En la sección Propiedades del filtro Mediana, seleccione 3D para Interpretación y haga clic en Aplicar. Esto creará un nuevo conjunto de datos filtrados (.filtered) en la sección Vista del proyecto .
- Haga clic con el botón derecho en los datos filtrados, busque Umbral interactivo y seleccione. En la sección Propiedades de Umbral interactivo, seleccione 3D para Tipo de vista previa, ajuste el valor inferior del Rango de intensidad a 2500 y haga clic en Aplicar. Esto creará un nuevo conjunto de datos con umbral (.thresholded) en la sección Vista del proyecto .
- Haga clic y arrastre la cadena de representación de volumen del conjunto de datos modificado al conjunto de datos con umbral. Desactive el umbral interactivo y active el aumento de volumen (si está desactivado).
- Segmentación de la imagen μCT
- Haga clic en el conjunto de datos filtrados y luego cambie a la pestaña Segmentación haciendo clic en Segmentación debajo de la barra de menú. Dentro de la pestaña Segmentación , los controles son similares a los de la pestaña Proyecto .
NOTA: Hay algunas herramientas adicionales relevantes para la ventana de visualización dentro de la pestaña Segmentación . Hay un botón de zoom +/- en la esquina superior de la ventana de visualización que se puede usar para que la imagen llene la ventana de visualización. La imagen se puede centrar en la ventana de visualización utilizando la barra de desplazamiento izquierda/derecha a lo largo de la parte inferior de la ventana de visualización. Encima de esta barra de desplazamiento, hay una segunda barra de desplazamiento para moverse a través de los segmentos de datos μCT.
- Para asegurarse de que la pestaña Segmentación esté configurada correctamente, en la sección Editor de segmentación , la imagen se establece en el conjunto de datos filtrados y el campo Etiqueta se establece en el conjunto de datos con umbral. En la sección Materiales , haga doble clic en el Material llamado Material y cámbiele el nombre a Hueso.
- Haga clic en Nuevo en el Editor de segmentación para crear un nuevo campo de etiqueta. En la ventana emergente, cambie el nombre utilizando una convención de nomenclatura similar (por ejemplo, 'Material de andamio óseo' - 'WKX' - Marcadores) como antes y como Marcadores y haga clic en Aceptar.
- Ajuste la ventana de visualización 2D utilizando vistas alternativas (XY, XZ o YZ) para mostrar una vista de plano sagital o coronal de la imagen.
- Utilice el pincel, el lazo, la varita y el umbral dentro de la sección Selección para identificar las áreas de interés. Para este protocolo, seleccione la herramienta pincel . Ajuste el tamaño de la herramienta de pincel ; Por lo general, un puntero de tamaño mediano funciona bien.
- Ajuste la pila de imágenes al primer corte en el que el hueso del radio comienza a ser visible. Esto puede requerir moverse a través de las rebanadas varias veces para determinar qué volúmenes pertenecen al radio y al cúbito y cuándo comienzan.
- Una vez identificado el corte, use el cursor del pincel para dibujar haciendo clic y arrastrando el cursor alrededor de este segmento de hueso. Aparecerá una línea roja donde se ha dibujado el cursor. Dibuje dentro de la sección del hueso, no fuera del hueso; Esto ayudará al software de visualización a diferenciar entre los dos huesos.
- Continúa moviéndote a través de las divisiones y dibujando alrededor del radio. Generalmente, a través de la sección media de la pila de imágenes, se pueden hacer dibujos cada 20-30 cortes de imagen. A lo largo del extremo de la pila de imágenes o alrededor del defecto de tamaño crítico, haga dibujos cada 5-10 cortes de imagen para ayudar al software de visualización a medida que el hueso se desvanece o se desvanece.
- Una vez que se haya completado la longitud completa del radio, haga clic en Selección, Relleno y luego en Todos los sectores en el menú desplegable de la barra de menú. Esto llenará todos los dibujos hechos con la herramienta pincel. Ahora aparecerán sombreados con rojo en lugar de contornos.
- Haga clic en Seleccionar e interpolar en el menú desplegable de la barra de menús. Esto rellenará las regiones sombreadas del hueso del radio en todas las rebanadas en función de los dibujos realizados.
NOTA: De vez en cuando se le escapan algunas secciones del hueso. Antes de continuar, desplácese por todas las rodajas para asegurarse de que generalmente haya identificado el hueso. Si se ha omitido una sección de hueso, use la herramienta de pincel para rodear esta región en algunas rodajas alrededor de la sección faltante y repita los dos pasos anteriores (1.3.9 y 1.3.10).
- Una vez que se haya identificado completamente el radio, en la sección Materiales , haga doble clic en el material llamado Interior y cámbiele el nombre a Radio.
- Haga clic en el símbolo Agregar (+) ubicado debajo de la sección Selección para agregar esto como un material. El radio ahora se delinea en la ventana de visualización. Haga clic en el icono de bloqueo para el material Radio en la sección Materiales antes de continuar; Esto asegurará que no se realicen cambios en este material.
- Completa el mismo proceso para el hueso cúbito. En la sección Materiales , haga clic en Agregar para crear un nuevo material. Haga doble clic en el nuevo material y cámbiele el nombre a Ulna.
- Repita los pasos 1.3.6 a 1.3.12 para el hueso cúbito.
- Una vez completados los pasos para los huesos Radio y Cúbito, desbloquee el Radio haciendo clic en el icono de Bloqueo para el material de Radio en la sección Materiales . Cambie a la pestaña Proyecto haciendo clic en Proyecto debajo de la barra de menú. Ahora hay un conjunto de datos denominado Marcadores aplicados al conjunto de datos filtrados en la sección Vista del proyecto.
- Haga clic con el botón derecho en el conjunto de datos filtrados y busque Cuenca hidrográfica basada en marcadores dentro de la marca (Segmentación de imagen) y seleccione . En la sección Propiedades , haciendo clic en los menús desplegables, establezca los Datos como el conjunto de datos filtrados, los Marcadores como conjuntos de datos de marcadores y la máscara binaria como el conjunto de datos con umbral y haga clic en aplicar.
- Se generará un archivo .grown dentro de la sección Vista del proyecto . Haga clic y arrastre la cadena de representación de volumen desde el conjunto de datos con umbrales hasta el conjunto de datos .grown para visualizar la segmentación del radio y el cúbito dentro de la ventana de visualización.
- Para convertir el archivo .grown a 8 bits, haga clic en el conjunto de datos .grown y busque Convertir tipo de imagen y seleccione. El archivo .grown se convertirá de 16 bits a 8 bits.
2. Alineación
- Alineación de ejes y extracción de radio y cúbito
NOTA: Si el modelo actual es el punto de tiempo inicial, se deben completar todos los pasos. De lo contrario, comience esta sección en el paso 2.1.4. Para permitir un corte transversal del radio, el modelo de punto de tiempo inicial debe alinearse perpendicularmente al plano de Ortho Slice; Este modelo servirá como línea de base para la alineación de las semanas siguientes.
- Dentro de la sección Vista de proyecto , active el Ortho Slice (si está desactivado) y, a continuación, haga clic en el conjunto de datos .grown. En la sección Propiedades , haga clic en el icono Editor de transformación . Con el cursor Interact, ajuste el ángulo del radio y el cúbito haciendo clic en los puntos de alineación del eje verde para que el corte ortogonal cree un corte transversal a través del radio.
- Una vez que se ha alineado, se debe crear un nuevo conjunto de datos para guardar esta transformación. Haga clic con el botón derecho en el conjunto de datos .grown, busque Resample Transformed Image y seleccione.
- Haga clic en Volver a muestrear imagen transformada y, dentro de la sección Propiedades , establezca los datos en el conjunto de datos .grown, la interpolación en el vecino más cercano, el modo en extendido, la opción conservar en el tamaño del vóxel y el valor de relleno en 0. Haga clic en Aplicar y se creará un nuevo conjunto de datos .transformado.
- Para extraer el radio y el cúbito del modelo de segmentación combinado, haga clic con el botón derecho en el conjunto de datos .transformed, busque Extraer etiqueta y seleccione.
- Haga clic en Extraer etiqueta, dentro de la sección Propiedades , establezca las etiquetas en el conjunto de datos .transformado , el ID de etiqueta en 1 y la opción Exportar a binario marcada. Haga clic en Aplicar y se creará un conjunto de datos de resultados.
- Haga clic y arrastre la cadena de representación de volumen desde el conjunto de datos .grown hasta el conjunto de datos Resultado . Haga clic en el conjunto de datos de resultados y presione F2 para cambiar el nombre del archivo (por ejemplo, 'Material de andamio óseo' - 'WKX' - Radio).
- Repita los pasos 2.1.4 a 2.1.6 para el cúbito estableciendo el ID de etiqueta en 2 y cambiando el nombre del archivo de resultados a Ulna.
- Guardado de archivos y apertura de imágenes longitudinales
NOTA: Si el modelo actual es el punto de tiempo inicial, se deben completar todos los pasos. De lo contrario, comience esta sección en el paso 2.2.2.
- Para guardar el archivo de punto de tiempo inicial, haga clic en Archivo, Guardar proyecto como y Establecer ubicación del directorio de archivos en el menú desplegable de la barra de menú. Establezca el tipo de guardado como un archivo de proyecto y archivos de datos (pack >) (*.hx) y asigne al archivo el nombre como la versión maestra (por ejemplo, 'Bone Scaffold Material'_ 'WKX' MASTER).
- Abra el archivo maestro de punto de tiempo inicial (si aún no está abierto) y guarde una copia haciendo clic en Archivo, Guardar proyecto como y Establecer ubicación del directorio de archivos en el menú desplegable de la barra de menú. Establezca el tipo de guardado como un archivo de proyecto y archivos de datos Amira (pack >) (*.hx) y asigne un nombre al archivo (por ejemplo, 'Material de andamio óseo'_ 'WK0' y 'WKX'). Este archivo se utilizará para comparar el punto de tiempo inicial con las semanas posteriores sin sobrescribir el archivo maestro.
- Dentro del nuevo archivo de comparación, abra las imágenes para el punto de tiempo de comparación haciendo clic en Abrir datos en la sección Vista del proyecto . Dentro de la ventana emergente, navegue hasta el directorio de carpetas deseado para seleccionar el punto de tiempo inicial .dcm archivos. Presione CTRL+A para seleccionar todas las imágenes y haga clic en Abrir imágenes.
- Una ventana emergente proporcionará información sobre el conjunto de datos de imagen abierto; haga clic en Aceptar para continuar.
- Aparecerá una advertencia sobre el escalado. Haga clic en Convertir en flotante (la opción predeterminada) para continuar.
- El icono del conjunto de datos aparecerá en la sección Vista del proyecto . Para cambiar el nombre del conjunto de datos, haga clic en el icono del conjunto de datos y presione F2; una vez que se ingresa el nombre deseado, haga clic en Aceptar. Se puede establecer un patrón de nomenclatura para identificar el escaneo (por ejemplo, 'Material de andamio óseo' - 'WKX').
- Alineación de modelos
NOTA: Este proceso discutirá la alineación de las segmentaciones óseas del radio. El mismo proceso se puede aplicar al cúbito, según sea necesario. No alinee ambos huesos al mismo tiempo; Esto causará problemas con la alineación adecuada.
- Desactive cualquiera de los puntos de tiempo iniciales Renderizado de volumen. No se muestra nada dentro de la ventana de visualización.
- Haga clic con el botón derecho en el conjunto de datos de comparación (semana siguiente), busque Ortho Slice y seleccione.
- Repita la sección 1 (comenzando en el paso 1.2), la sección 2 (paso 2.1) después de las notas adicionales para los números de paso proporcionados. Una vez finalizadas estas secciones, reanude en el paso 2.3.4. La finalización de estas secciones producirá la segmentación del conjunto de datos de la semana siguiente y extraerá los huesos del radio y el cúbito.
- Haga clic con el botón derecho en el conjunto de datos de comparación (semana siguiente) para el hueso de radio extraído, busque Asistente de registro de imágenes y seleccione. En la sección Propiedades , establezca Datos como el conjunto de datos de la semana de comparación para el hueso de radio extraído, Referencia como el conjunto de datos de punto de tiempo inicial para el hueso de radio extraído.
- En la sección Acciones del Asistente para registro de imágenes , haga clic en Omitir para el paso 1 de 4. Para los pasos 2 y 3 de 4, con el cursor de Interact, ajuste TabBox a la región común entre el punto de tiempo inicial y los conjuntos de datos de la semana de comparación haciendo clic en Aplicar en Acción después de cada paso. Para el paso 4 de 4, establezca Métrica como Correlación, Transformación como Rígida y Prealineación como Alinear ejes principales y haga clic en Aplicar en Acción.
- Una vez alineados los conjuntos de datos, haga clic con el botón derecho en el conjunto de datos de comparación (semana posterior) para el hueso de radio extraído y busque Remuestrear imagen transformada y seleccione. En la sección Propiedades , establezca los Datos como el conjunto de datos de la semana de comparación para el hueso de radio extraído, Interpolación al vecino más cercano, Modo en Extendido, Conservar en Tamaño de vóxel y Valor de relleno en 0 y haga clic en Aplicar. Se creará un nuevo conjunto de datos .transformed para el conjunto de datos de comparación.
3. Selección del ROI y cultivo
NOTA: Complete primero el recorte de ROI para determinar los números de corte alrededor de la fractura de tamaño crítico. Una vez que se determinan estos números de corte, estos se pueden usar en el mismo punto de tiempo para el cúbito, según sea necesario. Para este proceso, el recorte no se puede revertir una vez aplicado al modelo.
- Región de interés cultivo
- Haga clic en y active la división de ortofotos para el punto de tiempo inicial y establezca los datos en el conjunto de datos de punto de tiempo inicial para el radio extraído. Establezca la Orientación de modo que el plano produzca un corte transversal a través del hueso del radio.
- Con la barra deslizante Número de corte de la sección Propiedades , cambie el número de corte para determinar las ubicaciones de corte proximal y distal que rodean el defecto de tamaño crítico. Para ambos cortes, determine el corte más distal y proximal donde la fractura se encuentra con la diáfisis del hueso del radio. Documente el número de rebanada; Este número variará para cada modelo.
- Haga clic en y active el Ortho Slice para la semana de comparación y establezca los Datos en el conjunto de datos de punto de tiempo inicial para el radio extraído. Establezca la Orientación para que el plano produzca un corte transversal a través del hueso del radio.
- Con la barra deslizante Número de corte de la sección Propiedades y con el conjunto de datos de punto de tiempo inicial que muestra el corte ortodoxo distal, cambie el número de corte de la semana de comparación para que se alinee con el segmento distal del punto de tiempo inicial. El software de visualización permite una comprobación visual de la alineación; Cuando los sectores se superponen, aparecerán como un solo sector en la ventana de visualización.
- Anote el número de corte para el corte distal del conjunto de datos de la semana de comparación. Repita el paso 3.1.4 para el corte proximal.
NOTA: Como comprobación adicional, determine la diferencia entre el corte proximal y distal para los conjuntos de datos de la semana inicial y de comparación. Estos valores deben ser equivalentes, lo que ilustra un ROI de tamaño constante.
- Haga clic en el punto de tiempo inicial para el radio extraído. En la sección Propiedades , haga clic en la herramienta Editor de recortes .
- Dentro de la ventana emergente Editor de recortes , ingrese los valores mínimo y máximo en el campo X, Y o Z apropiado (la ventana de visualización alterará el ROI a medida que se ingresen). Una vez introducidos los valores, haga clic en Aceptar. El conjunto de datos se recortará al ROI de entrada.
- Repita el paso 3.1.7 para el conjunto de datos de la semana de comparación.
4. Análisis y visualización
- Análisis de volumen
- Para determinar el volumen del conjunto de datos de punto de tiempo inicial, haga clic con el botón derecho en el conjunto de datos de punto de tiempo inicial transformado para el radio extraído, busque Estadísticas de materiales y seleccione. Dentro de la sección Propiedades , establezca Datos como el punto de tiempo inicial transformado para el radio extraído y, defina Seleccionar como Materiales y haga clic en Aplicar. Un nuevo archivo . Se creará un conjunto de datos MaterialStatistics para el conjunto de datos de punto de tiempo inicial.
- Haga clic en el archivo . MaterialStatistics para el punto de tiempo inicial para el radio extraído y dentro de la ventana Propiedades , haga clic en Mostrar hoja de cálculo. Haga clic en la pestaña Tablas encima de la ventana Propiedades . Dentro de esta pestaña, se muestra el volumen del conjunto de datos de punto de tiempo inicial recortado para el radio extraído.
- Repita los pasos 4.1.1 y 4.1.2 para el conjunto de datos de la semana de comparación para el radio extraído. Dentro de la pestaña Tabla , el punto de tiempo inicial y la semana de comparación tendrán pestañas separadas en las que se puede hacer clic.
- Visualización de conjuntos de datos
- Para visualizar el cambio en el volumen óseo, haga clic con el botón derecho en el conjunto de datos transformados de la semana de comparación para el radio extraído, busque Aritmética y seleccione. Haga clic en Aritmética y dentro de la ventana Propiedades , establezca la Entrada A como el conjunto de datos transformados de la semana de comparación para el radio extraído, establezca la Entrada B como el conjunto de datos de punto de tiempo inicial para el radio extraído, establezca la Entrada C como SIN FUENTE, establezca el Tipo de resultado como Entrada A, deje la opción sin marcar, establezca los canales de resultado como como la entrada A, y establezca Expresión como A-B.
- Haga clic en el conjunto de datos de resultados y presione F2 para cambiar el nombre del archivo (por ejemplo, 'Material de andamio óseo' - 'WKX' - Cambio de hueso). Haga clic con el botón derecho en este conjunto de datos de resultados, busque Generar superficie y seleccione. Haga clic en Generar superficie y, dentro de la ventana Propiedades , haga clic en Aplicar y, en la ventana emergente, haga clic en Continuar. Se creará un nuevo conjunto de datos .surf.
- Haga clic con el botón derecho en el conjunto de datos .surf , busque Vista de superficie y seleccione. Se mostrará una vista de superficie del resultado aritmético en la ventana de visualización.
- Para cambiar el color de la vista de superficie, haga clic en la vista de superficie en la ventana Vista de proyecto y, a continuación, en la ventana de Propiedades , haga clic en el menú desplegable de Colores y seleccione Constante. Haga clic en el mapa de colores y asigne el color preferido.
- Para ver el cambio de volumen óseo en el conjunto de datos de la semana inicial, haga clic con el botón derecho en el conjunto de datos .transformado, busque Extraer etiqueta y seleccione. Haga clic en Extraer etiqueta y, dentro de la sección Propiedades , establezca las etiquetas en el conjunto de datos .transformado , el ID de etiqueta en 1 y la opción Exportar a binario activada. Haga clic en Aplicar y se creará un conjunto de datos de resultados . Haga clic en el conjunto de datos de resultados y presione F2 para cambiar el nombre del archivo (por ejemplo, 'Material de andamio óseo' - 'WKX' - Completo - Radio).
- Haga clic con el botón derecho en este conjunto de datos de resultados, busque Generar superficie y seleccione. Haga clic en Generar superficie y, dentro de la ventana Propiedades , haga clic en Aplicar y en la ventana emergente, haga clic en Continuar. Se creará un nuevo conjunto de datos .surf.
- Haga clic con el botón derecho en el conjunto de datos .surf, busque Vista de superficie y seleccione. Se mostrará una vista de superficie del resultado aritmético en la ventana de visualización.
- Para cambiar el color de la vista de superficie, haga clic en la vista de superficie en la ventana Vista de proyecto y, a continuación, en la ventana de Propiedades, haga clic en el menú desplegable de Colores y seleccione Constante. Haga clic en el mapa de colores y asigne el color preferido.
- Para continuar con semanas adicionales, vuelva al paso 2.2.2.