Summary
Presentamos un flujo de trabajo para segmentar y cuantificar los huesos trabeculares para imágenes 2D y 3D basados en el límite exterior del hueso utilizando un plugin de ImageJ. Este enfoque es más eficiente y preciso que el actual enfoque de contorno de mano manual y proporciona cuantificaciones capa por capa, que no están disponibles en software comercial actual.
Abstract
Tomografía computada de micro (micro-CT) se utiliza rutinariamente para evaluar la cantidad de hueso y trabeculares propiedades microestructurales en pequeños animales en condiciones de pérdida de hueso diferentes. Sin embargo, el enfoque estándar para análisis trabecular de imágenes micro-CT es rebanada por rebanada semiautomática mano contorno, uso intensivo de mano de obra y propenso a errores. Se describe aquí es un método eficaz para la segmentación automática de huesos trabeculares según límites exterior del hueso, donde huesos trabeculares pueden ser identificados y segmentados automáticamente con precisión con menos sesgo del operador cuando sea apropiado se establecen parámetros de segmentación. Para Perfil de parámetros de segmentación satisfactoria, se muestra una pila de imágenes de los resultados de la segmentación, donde todas las combinaciones posibles de los parámetros de segmentación son cambiados uno por uno en secuencia, y resultados de segmentación con parámetros asociados pueden fácilmente se comprueba visualmente. Como medida de control de calidad del plugin, se cuantifican objetos estándar simulados donde pueden comparar las cantidades medidas con valores teóricos. Capa por capa la cuantificación de propiedades trabeculares y el grosor trabecular están reportados por tal un plugin, y la distribución de tales propiedades dentro de las regiones seleccionadas puede ser perfilada fácilmente. Aunque la cuantificación capa por capa retiene más información sobre huesos trabeculares y facilita aún más el análisis estadístico de los cambios estructurales, tales medidas son de carácter de la salida de software comercial actual, donde solamente una sola valor cuantificado para cada parámetro se divulga para cada muestra. Por lo tanto, los flujos de trabajo descritas son mejores enfoques para el análisis de los huesos trabeculares con precisión y eficiencia.
Introduction
Micro-CT en el análisis de huesos trabeculares es el acercamiento estándar para el seguimiento de los cambios morfológicos de los huesos en animales pequeños bajo diferentes hueso pérdida condiciones1,2,3, donde diversas variables relacionadas con la estructuras de los huesos son reportados4. Sin embargo, estos parámetros no se distribuyen uniformemente en la metafísis de huesos largos5, y sólo un valor resumido o promedio es para cada variable estructural de cada muestra por comercial máquinas de micro-CT actual6,7 , aunque un solo valor no puede representar completamente las características del parámetro medido en la región de análisis. Cuantificación de capa por capa de huesos trabeculares no sólo retiene más información para cada variable, sino que también permite el perfilado de las distribuciones de variables en la región de análisis, facilitar el posterior análisis estadístico de estructural cambios bajo diferentes condiciones de5. Por lo tanto, el objetivo de este método es cuantificar los huesos trabeculares de micro-CT explora en cada nivel de sector, que no está disponible en cualquier paquete de análisis micro-CT comercialmente disponibles actualmente.
Eficientemente segmento huesos trabecular rebanada-de-rebanar, métodos de segmentación automática son deseables. Sin embargo, la técnica estándar actual para el análisis micro-CT se basa en contorno interactivo manual seguido por interpolación semiautomático para separar los huesos trabeculares de los compartimientos corticales, que es mano de obra, errores, y asociada a importante operador diagonal8,9,10. Se informaron los métodos de segmentación automática11,12 , pero estos métodos sólo son óptimos en regiones con buena separación entre huesos trabeculares y corticales huesos, pero no en las regiones sin separaciones claras. Por otra parte, se requieren para diferentes muestras12parámetros de segmentación diferente, y es tedioso seleccionar manualmente los parámetros de segmentación satisfactorio aplicables a grupos de muestras de huesos por probar distintos de combinaciones de parámetros12, a pesar de que el proceso de segmentación es automático cuando se establecen todos los parámetros relacionados. Como el límite exterior del hueso tiene el mayor contraste con el fondo de exploración y las cáscaras cortical metafisaria de los huesos largos muestran pocos cambios en el elegido analizar la región métodos de segmentación según el contorno del límite exterior de los huesos largos pueden forma fiable y precisa separar huesos trabeculares conchas corticales. La ventaja de este método de segmentación es que la segmentación se basa en la diferencia entre el fondo y el límite exterior del hueso, pero no en las diferencias entre los huesos trabeculares y corticales6,12, 13, por lo tanto es generalmente fácil de encontrar una combinación de parámetros de segmentación que es satisfactoria para un grupo de muestras de huesos, facilitando el análisis más fiable de trabeculares cambios entre diferentes grupos.
En cada rebanada de nivel, área, perímetro y espesor (2D) bidimensional se divulgan para análisis 2D, mientras que el volumen, superficie y grosor de (3D) tridimensional se divulgan en cuantificaciones 3D. Dicha información no se divulga generalmente por las actuales herramientas de análisis de imagen, lo que indica que los procedimientos reportados pueden aplicarse a imágenes general donde se desea dicha información.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
Procedimientos con sujetos animales fueron conducidos de acuerdo con la guía para el cuidado y uso de animales de laboratorio (publicación de NIH, 8ª edición, 2011) y han sido revisados y aprobada por el cuidado de Animal institucional y Comité uso de Wuhan Universidad.
1. instalación del software
- Instalar el software ImageJ. Descargar la versión para Windows del software ImageJ (versión p 1,51) liado con Java de 64 bits de https://imagej.nih.gov/ij/. Extraiga el software descargado en una carpeta, que se hará referencia posteriormente a como el "directorio de ImageJ".
Nota: Plugins Trabecular análisis requieren java 64 bits ejecutar tiempo (versión 1.8) y un sistema operativo de 64 bits de windows, preferiblemente un 64-bit sistema operativo Windows 7. - Instalar plugins de análisis trabecular. Solicitar análisis trabecular plugins de http://www.bomomics.com y descomprimir el contenido en el directorio de plugins de ImageJ, que es "ImageJ directorio/plugins".
Nota: El plugin puede obtenerse una versión libre, donde las mediciones de 5 especifican sectores adyacentes se divulgan, o como una versión comercial, donde se puede especifica y mide una gama de sectores.
2. prepare el conjunto de datos 3D para análisis Trabecular
- Exploración fémur de rata con un siguiente máquina micro-CT en el estándar de análisis de protocolo5 y luego guardar los datos en un formato que puede importarse en ImageJ, e.g., formato tiff. Si múltiples muestras de hueso fueron analizadas simultáneamente en un tubo de muestra, primero importar los datos en ImageJ y luego separar cada hueso individual por cultivo de otras muestras usando herramientas de procesamiento de imágenes de ImageJ. Posteriormente, guardar la imagen resultante en un formato que puede importarse en ImageJ más adelante.
Nota: Un archivo de imagen de muestra representativa en el análisis se incluye (1 archivo suplementario). -
Simulación de objetos 2D y 3D estándar.
- Abra el software ImageJ. Bajo los Plugins de | BoMomics | Simular objetos menú, haga clic en el botón círculo , escriba 200 el diámetro en la ventana emergente y haga clic en OK para generar el círculo simulado con un diámetro de 200 píxeles (figura 3B). Guardar el círculo generado en formato tiff.
- Abra el software ImageJ. Bajo los Plugins de | BoMomics | Simular objetos menú, haga clic en el botón cuadrado , introduzca 200 como la longitud del lado de la ventana emergente y haga clic en OK para generar el cuadrado simulado con una longitud lateral de 200 píxeles (figura 3B). Guardar la plaza generada en formato tiff.
- Abra el software ImageJ. Bajo los Plugins de | BoMomics | Simular objetos menú, haga clic en el botón de rectángulo , introduzca 200 como la anchura y 100 como la altura de la ventana emergente y haga clic en OK para generar el rectángulo simulado con un ancho de 200 píxeles y una altura de 100 píxeles (figura 3 B). Guardar el rectángulo generado en formato tiff.
Nota: Un círculo (diámetro: 200 píxeles), un cuadrado (longitud del lado: 200 píxeles) y un rectángulo (ancho: 200 píxeles, altura: 100 píxeles) se guardan para los análisis posteriores. - Abra el software ImageJ. Bajo los Plugins de | BoMomics | Simular objetos menú, haga clic en el botón del cubo , introducir 30 como la longitud del lado de la ventana emergente y haga clic en OK para generar el cubo simulado. Por último, haga clic en el Plugins | 3D | Visor de volumen a ver el cubo generado y guardar en formato tiff (figura 3C).
- Abra el software ImageJ. Bajo los Plugins de | BoMomics | Simular objetos menú, haga clic en el botón de rectángulo , introduzca 80 longitud 40 como el ancho y 30 como la altura de la ventana emergente, y haga clic en OK para generar el cubo simulado. Por último, haga clic en el Plugins | 3D | Visor de volumen a ver el cubo generado y guardar en formato tiff (figura 3C).
- Abra el software ImageJ. Bajo los Plugins de | BoMomics | Simular objetos menú, haga clic en el botón esfera , introducir 30 como el diámetro en la ventana emergente y haga clic en OK para generar el ámbito simulado. Por último, haga clic en el Plugins | 3D | Visor de volumen a ver el ámbito generado y guardar en formato tiff (figura 3C).
- Abra el software ImageJ. Bajo los Plugins de | BoMomics | Simular objetos menú, haga clic en el botón cilindro , introducir 30 como el diámetro y 100 como la altura de la ventana emergente y haga clic en OK para generar el cilindro simulado. Por último, haga clic en el Plugins | 3D | Visor de volumen a ver el cilindro generado y guardar en formato tiff (figura 3C).
Nota: Una esfera (diámetro: 30 píxeles), un cubo (longitud lateral: 30 píxeles), un ortoedro (longitud: 80 píxeles, ancho: 40 píxeles; Altura: 30 píxeles) y un cilindro (diámetro: 30 pixeles, altura: 100 píxeles) se guardan para los análisis posteriores.
3. parámetros de análisis de perfiles
- Abra el software ImageJ y abrir o importar una imagen escaneada.
- Deslice la barra de desplazamiento inferior para elegir una rebanada, luego haga clic en imagen | Ajustar | Umbral botón. En la ventana de umbral, ajustar los valores de umbral mínimo y máximo por inspección manual para asegurarse de que los huesos están bien separados del fondo y registre el valor de umbral mínimo como el valor de umbral de hueso cortical.
Nota: Como fabricantes de micro-CT en diferentes utilizan diferentes factores de escala para almacenar los coeficientes de atenuación de rayos x en las imágenes, el valor real debe determinarse empíricamente o por inspección manual o seguir el fabricante específico recomendaciones. En nuestra práctica, ajuste el umbral mínimo de ISQ archivos producidos por las máquinas de micro-CT a 6.000-7.000 puede confiablemente separar huesos desde el fondo de la exploración. - Haga clic en los Plugins de | BoMomics | Perfiles de trab Param botón. En la ventana emergente, establecer Índice de Slice a la posición del segmento representativo y la cortical ósea ("hueso de la Cort") ("gama") y valores ("paso") para el cálculo de un conjunto de umbrales corticales para Perfil de parámetros de segmentación, donde el valor de umbral de hueso cortical se adquiere al paso 3.2 (figura 1). Los valores por defecto gama y paso de 2.000 y 400 para la mayoría de las imágenes micro-CT. Establecer el Índice de corte a 5 para el conjunto de datos de ejemplo y mantener los valores por defecto para otros ajustes.
Nota: "Cort Bone" es el umbral para los huesos corticales y el umbral utilizado para cambios perfilado parámetro desde el valor umbral más bajo para el valor de umbral más alto con el incremento del valor del paso, donde el valor de umbral más bajo es el valor de Cort Hueso - gama y el mayor valor de umbral es Cort hueso + gama. Si los valores de set para hueso de Cort, gama y paso son 6.000, 1.000 y 500, respectivamente, entonces el valor de umbral más bajo para los huesos corticales es 6.000-1, 000 = 5, 000, el valor de umbral más alto es 6, 000 + 1, 000 = 7, 000 y los umbrales utilizan en los análisis son 5.000, 5.500 6.000, 6.500 y 7.000. - Ajustar diámetro de ruido ("diámetro del ruido"), paso ("paso") y rango de valores ("gama") para especificar un conjunto de valores de ruido en los análisis y el diámetro del agujero "(diámetro del agujero del) (" paso ") y valores (" gama ") para el cálculo de un conjunto de valores de agujero. En general, la configuración por defecto funciona para la mayoría de las muestras de hueso de micro-CT, donde, paso, el ruido y gama son 5, 5, 2, respectivamente, y gama, de paso y diámetro de agujero son 15, 5, 2, respectivamente. Mantenga la configuración predeterminada para el conjunto de datos de ejemplo.
Nota: "Diámetro de ruido" es el diámetro del filtro de supresión de ruido, y el "diámetro de agujero" es el diámetro de agujeros en la corteza de huesos. Los valores de ruido y agujero utilizados para perfiles de parámetro pueden ser calculados semejantemente como valores umbral cortical descritos anteriormente, utilizando el orificio de ruido especificado, gama y los valores de paso. Debido a las variaciones en el análisis parámetros al adquirir las imágenes, parámetros de ruido y el agujero deben ser determinados empíricamente según la calidad de imagen, como no general rangos de parámetros de ruido y el agujero están buenos para todos los escaneos. El valor más bajo umbral, ruido o agujeros debe ser mayor o igual a 0, y si el valor mínimo calculado es negativo mediante los parámetros proporcionados, el menor valor específico se establece en 0. - Haga clic en OK para realizar perfiles de parámetro. Visualmente Compruebe los resultados de la segmentación en la ventana de Resultados de perfiles de parámetro y seleccione una capa de la rebanada para que el límite exterior del hueso se describe con bastante precisión (figura 1B). Posteriormente, recuperar los parámetros de generación de perfiles a partir de la entrada en la tabla de Resultados de perfiles del parámetro correspondiente a la capa del segmento elegido (tabla 1).
4. trabecular análisis
-
Segmentación de los huesos trabeculares
- Abra el software ImageJ, luego abrir o importar una imagen escaneada.
- Haga clic en los Plugins de | BoMomics | Segmentación de trab botón y rellene los parámetros de análisis apropiado. Sistema de "Inicio", "Límite de contorno", "Trab. Huesos","Diámetro de reducción de ruido","Agujero de llenado de diámetro"y"Diámetro de espesor Cortical"5, 7.200, 7.000, 6, 12 y 25, respectivamente.
Nota: "Inicio" y "Fin" especifica el rango de división seleccionada para la segmentación de los huesos trabeculares, "Esquema límite" es correspondiente al perfilada parámetro "Cort Bone", "Ruido reducción diámetro" a "Diámetro ruido." Parámetro y "Agujero de llenado diámetro" parámetro "Diámetro de agujero". "Diámetro de grosor cortical" es el espesor especificado para excluir los huesos corticales externos. «Trab. Los huesos "es el umbral para la extracción de los huesos trabeculares (figura 2), donde los valores se determinan usando perfiles comando parámetros de análisis como se describe en el paso 3.5. - Haga clic en Aceptar para realizar la segmentación trabecular. Revise los resultados de la segmentación en la ventana de Resultados de segmentación de Trab (figura 2B). Guardar los huesos trabeculares extraídos que se muestra en la ventana de Huesos Trabecular segmentado en formato tiff (figura 2B), que puede ser más analizada por otro software.
-
Análisis de los huesos trabeculares.
- Abra el software ImageJ, luego abrir o importar una imagen escaneada.
- Haga clic en los Plugins de | BoMomics | Análisis de trab botón y rellene parámetros de análisis apropiado, como "Inicio", "Límite de contorno", "Trab. Los huesos","Diámetro de reducción de ruido","Agujero de llenado de diámetro"y"Diámetro de espesor Cortical"como describen anteriormente (figura 3A), donde se determinan los parámetros el comando perfiles análisis parámetros como se describe en paso 3.5. Sistema de "Inicio", "Límite de contorno", "Trab. Huesos","Diámetro de reducción de ruido","Agujero de llenado de diámetro"y"Diámetro de espesor Cortical"5, 7.200, 7.000, 6, 12 y 25, respectivamente.
Nota: En la versión gratuita del plugin, se seleccionan cinco rebanadas adyacentes desde el índice especificado del sector de "Inicio" para la medición, mientras que en la versión comercial, un número arbitrario de segmentos puede ser especificado por el usuario. - Seleccione una o más opciones en la sección de Informes de resultados para los parámetros a medir, donde el volumen de hueso trabecular (BV), volumen total (TV) de la región seleccionada y espesores medidos ya sea representan (2D) o tridimensional (3D) están disponibles para la selección a través de tres casillas de verificación, es decir, "BV sólo TV", "2D" y "3D". Seleccione las casillas de "2D" y "3D", haga clic en "Aceptar" para realizar análisis trabecular (figura 3, tabla 2).
Nota: Cuando "BV TV sólo", sin importar el estado de selección de "2D" y "3D", se reportan medidas crudas de BV, TV y la intensidad, y los huesos trabeculares segmentados son extraídos y aparecerá en una nueva ventana, que puede guardar y más analizada por otras software. Cuando está seleccionada "2D", se divulgan medidas crudas de BV, TV, intensidad y espesores medidos representan en cada nivel de segmento utilizando el modelo de placa. Si está marcado "3D", espesor tridimensional para cada voxel se calcula directamente sin cualquier asunción del modelo, luego crudas medidas de BV, TV, intensidad y espesor tridimensional muestreados en cada nivel de la rodaja se divulga. Si ninguna opción es seleccionada, huesos trabeculares están segmentados mediante el anterior conjunto de parámetros, pero no se informó de ninguna medición.
5. cuantificar objetos simulados
- Software ImageJ, y abra una imagen simulada. Aquí, abrir la esfera simulada con un diámetro de 30 píxeles por ejemplo.
- Seleccionar los Plugins de | BoMomics | Análisis de trab botón y rellene los parámetros adecuados de análisis como se describe anteriormente. Mantenga los valores predeterminados de "Inicio", "Fin", "Límite de contorno", "Trab. Huesos"y"Diámetro de reducción de ruido","Diámetro de reducción de ruido","Agujero de llenado de diámetro"y"Diámetro de espesor Cortical"a 0 (figura 3).
Nota: Para objetos simulados, no ruido de señal no deseados, y no pueden existir cáscaras corticales correspondientes. Por lo tanto, se deben establecer parámetros para tales valores apropiadamente (por defecto está a cero). Para cuantificaciones 3D, 30 rebanadas antes y después el segmento especificado "Start" se procesan utilizando la versión gratuita del plugin, mientras que sólo rebanadas en el rango especificado por rebanadas de "Inicio" y "Fin" son analizadas por la versión comercial del plugin. - En la sección de Informes de resultados , seleccione "2D" y "3D" para los parámetros a medir y haga clic en OK para realizar análisis de trabecular para el objeto simulado (tabla 3).
6. calibración de medidas Trabecular y presentación de datos: Perfil de la distribución de medidas Trabecular en la región analizando las
- Obtener la información de calibración del conjunto de datos escaneado, según las instrucciones del proveedor de micro-CT.
Nota: Solamente crudas medidas de BV, TV, intensidad y grosores son registradas por el plugin. Para obtener medidas comparables a los informes generados por otros programas, calibración del resultado es el deseado. - Abra Microsoft Excel y abra la tabla de resultados informados. Calcular calibrado hueso volumen (BV), total volumen (TV), contenido mineral óseo (BMC), fracción de volumen de hueso (BV/TV) y densidad mineral ósea (DMO) en nuevo excel columnas utilizando los valores reportados de BV, TV y la intensidad crudos según las siguientes expresiones.
- Nota: Escaneo resolución (resolución, μm), escala de grises (escala), unidades de densidad (mg HA/cm3), pendiente (slope) de la densidad y intercept (interceptar) de la densidad pueden ser extraídos de las micro-CT en imagen o meta archivos escaneados, como el archivo Scanco ISQ. Por lo tanto, las medidas calibradas se calculan como sigue:
Suponga BV, la TV y la intensidad son medidas crudas, BVc y TVc son valores calibrados, resolución, escala, pendiente, intercepto son de archivos de calibración de micro - CT imágenes.
BVc = resolución de × BV3 [μm3]
TVc = resolución de TV ×3 [μm3]
BV / TV = BV ÷ TV
BMC = (intensidad ÷ escala × pendiente - intercepción de × BV) resolución ×3 × 10-12 [mg HA]
DMO - BMC ÷ TVc × 1012 [mg HA / cm3] - Crear un diagrama XY (dispersión) medidas calibrado (Y) contra capas de corte (X) utilizando el software de Microsoft Excel (figura 4).
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
El plugin de análisis trabecular está diseñado para segmentar y cuantificar huesos trabeculares con exactitud automáticamente. Inicialmente, ósea límite externo es detectado y delineado seguida de una operación de relleno de agujero donde se llenan los orificios dentro de cáscaras cortical externa del hueso. Luego se realiza una operación de erosión para excluir los huesos corticales externos y los huesos trabeculares segmentados. Finalmente, se cuantifican las medidas de huesos trabeculares en la región segmentada.
Como micro-CT en imágenes es inherentemente ruidosos, segmentaciones con parámetros predefinidos de arbitrarios a menudo fallan identificar con precisión los límites exteriores del hueso. Es tedioso y mano de obra intensiva para probar un montón de combinaciones de parámetros para la selección de parámetros de segmentación satisfactoria. Por lo tanto, se proporciona un parámetro perfiles plugin para cambiar los parámetros uno a uno en la gama ajuste automáticamente para ayudar en la selección de combinaciones de parámetros satisfactorios, que facilita la selección de un conjunto común de parámetros de un grupo de hueso muestras. Figura 1 A muestra la configuración utilizada para perfilar los parámetros buena segmentación. Cuando los parámetros de umbral de hueso cortical (ósea Cort), rango de umbral a ser perfilados (gama), y se especifican el monto del incremento de umbral (paso) en cada paso, se generan una serie de umbrales a ser perfilados. Posteriormente, una serie de valores de ruido y el agujero se generan de manera similar mediante el establecimiento de los parámetros correspondientes. Finalmente, los huesos trabeculares están segmentados cambiando los parámetros de uno en uno para todas las combinaciones posibles del parámetro. Figura 1 B muestra los resultados de segmentación representativas para las combinaciones de diferentes parámetros. Obviamente, algunas combinaciones de parámetro son mejores que otros en delinear límites externa de un hueso, y más de una combinación de parámetros muestra resultados satisfactorios de la segmentación. Después de comprobar visualmente los resultados de segmentación, valores de los parámetros de segmentación satisfactoria se pueden recuperar de la tabla de resultados de generación de perfiles (tabla 1).
Para cuantificar las medidas de los huesos trabeculares, se realizan análisis y segmentación trabecular. Figura 2 A muestra el cuadro de diálogo Configuración de segmentación trabecular, donde huesos trabeculares en la región seleccionada son segmentados y extraídos (figura 2B), y pueden comprobarse resultados segmentación mediante los parámetros proporcionados rebanada por rebanada visualmente. Posteriormente, se analizan los huesos trabeculares con parámetros satisfactorios (figura 3A). Dependiendo del estado de selección de informes opciones raws cuantificaciones de volumen de hueso (BV), volumen (TV), suma de valores de gris (intensidad), total y espesores medición o dos dimensiones o tres dimensiones (tabla 2) se divulgan. Por último, información de calibración se extrae del conjunto de datos de micro-CT escaneado y calibrados medidas de BV, TV, BMC, BV/TV y DMO se calculan seguido perfilando sus distribuciones en el análisis región capa-por-capa seleccionada contra la capa de posiciones (figura 4).
Como medida de control de calidad del plugin, se admite cuantificación de objetos simulados. Simulación estándar objetos con dimensiones conocidas se cuantifican por el plugin para la comparación con los valores teóricos o las medidas de otro software, como software comercial suministrado con las máquinas de micro-CT o BoneJ14, un plugin gratuito de código abierto para el análisis de la imagen del hueso. Para analizar objetos simulados, el ruido para ajustar a cero, como imágenes simuladas se consideran imágenes de alta calidad sin ninguÌ n ruido. Se simularon objetos con diferentes grosores y se muestran los resultados (figura 3, tabla 3). Para simulados objetos 2D estándar, tales como círculos, cuadrados y rectángulos, exacto valores de área (TV o VB) y grueso son registrados (tabla 3). Para los objetos 3D, medidas de espesor exacto de cubos, esferas y cuboides se reportan, sin embargo, los espesores de cilindros no son precisos para voxels cerca de ambos extremos del cilindro, mientras que los espesores de voxel en rodajas de mediados de los cilindros son exactamente como predijo. Esta es una característica del algoritmo de medición de espesor subyacente, en espesor de voxel para cada objeto está determinado por el diámetro de la esfera más grande, o la longitud del lado del cubo más grande, que contiene este voxel y está totalmente dentro de la objeto. Por lo tanto, espesores de los objetos de diferentes esferas y cubos pueden medirse con precisión, mientras que los cilindros sólo se pueden medir exactamente en radio-distancia medio de la división de ambos extremos.
Figura 1 : Resultados representativos del parámetro perfiles análisis. (A) los parámetros de configuración de página. (B) resultados representativos del parámetro perfiles de análisis. Algunas combinaciones de parámetro son mejores que otros para la detección de las fronteras exteriores de los huesos. Haga clic aquí para ver una versión más grande de esta figura.
Figura 2 : Resultados representativos del análisis de segmentación trabecular. (A) los parámetros de configuración de página. (B) resultados de segmentación representante de huesos trabeculares en diferentes capas. Haga clic aquí para ver una versión más grande de esta figura.
Figura 3 : Análisis trabecular. (A) los parámetros de configuración de página. (B) resultados representativos de objetos 2D simulados. (C) resultados representativos de objetos 3D simulados visualizados por visor de volumen 3D. Haga clic aquí para ver una versión más grande de esta figura.
Figura 4 : Distribuciones de medidas de hueso trabecular en la región de análisis seleccionada. El eje horizontal representa la distancia relativa a la capa de corte inicial en la región de análisis. Los valores en el eje y son calibradas medidas trabeculares en la región analiza. Haga clic aquí para ver una versión más grande de esta figura.
| Slice | Límite de contorno | Diámetro de reducción de ruido | Diámetro de agujero relleno. |
| 4 | 5200 | 0 | 10 |
| 147 | 5200 | 2 | 20 |
| 361 | 7200 | 6 | 16 |
| 539 | 8000 | 10 | 20 |
Tabla 1: Resultados representativos del parámetro perfiles de análisis.
Tabla 2: Resultados representativos del análisis trabecular.
| Objeto | Dimensión | Volumenb | Superficiec | Espesor |
| Plaza | 200 X 200 | 40000 | 796 | 200 |
| Rectángulo de | 200 X 100 | 20000 | 596 | 100 |
| Círculo | Diámetro: 200 | 31428 | 796 | 200 |
| Cubo | 30 X 30 X 30 | 27000 | 5048 | 30 |
| Cuboid | 80 X 40 X 30 | 96000 | 13008 | 30 |
| Esfera | DIA: 30 | 14328 | 3944 | 30 |
| Cilindrod | DIA:30; H: 100 | 71600 | 12800 | 27.84 |
| Cilindroe | DIA:30; H: 100 | 51552 | 9552 | 30 |
| r: los resultados son en vóxeles crudo. Diámetro: diámetro; H: altura. | ||||
| b: volumen (3D) o área (2D). | ||||
| c: superficie (3D) o perímetro (2D). | ||||
| d: Measurs de segmento 1 segmento 100 se utilizan para el análisis. | ||||
| e: Measurs de laminar 15 cortar 85 se utilizan para el análisis. | ||||
Tabla 3: Resultados de la cuantificación de objetos simulados.
Archivo complementario 1. Muestra ósea. Haga clic aquí para descargar este archivo.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
Este estudio describe un plugin de ImageJ para el análisis de los huesos trabeculares, que es automático, eficiente, y fácil de usar. El plugin también puede utilizarse para cuantificar objetos 2D o 3D, capa por capa medidas de áreas, volúmenes y espesores. Actualmente, sólo un único valor medido para cada parámetro trabecular se divulga para cada muestra por análisis estándar de micro-CT, que no pueden representar completamente las características de la entidad medida en la región de análisis seleccionada. El plugin se describe informes de cantidades de capa por capa de cada parámetro para cada muestra, que conserva totalmente la información de la distribución del parámetro medido en la región de análisis seleccionada y por lo tanto más avanzada y sensible a la estadística planteamientos son aplicables para el análisis de dichos datos.
Con el semi automática, estándar mano-segmentación método, observaron una variación sustancial entre diferentes sectores que contornea la mano con operador diagonal8,9,10, garantizando una más uniforme y método automático para segmentación trabecular. Como el contraste es máximo entre el fondo de la exploración y el límite exterior del hueso, límite exterior asistida trabecular segmentación se realiza por el plugin, donde el límite exterior del hueso se puede detectar automáticamente con precisión cuando corresponda umbral, los ruidos y ajustes de agujero se proporcionan. Para facilitar la determinación de parámetros de segmentación, una gama de combinaciones de parámetros perfila uno por uno con una porción representativa de la región de análisis, y resultados satisfactorios de la segmentación se comprueban visualmente. Posteriormente, se analizan los huesos trabeculares con los parámetros de perfilado. Cuando se programan los parámetros adecuados de segmentación, y el plugin se ha aplicado con éxito para analizar huesos trabeculares de fémur distal de rata tratadas con hueso diferentes agentes anabolizantes5, se reportan resultados reproducibles para la misma imagen sin sesgo del operador.
El plugin en la actualidad sólo puede procesar una muestra en un momento. Si las muestras múltiples están presentes en la imagen, todas las muestras se cuantifican como un solo objeto sin discriminación. Por lo tanto, preprocesamiento es necesaria si las muestras múltiples se analizan simultáneamente en un tubo de muestra. Las salidas de la cuantificación son medidas crudas de pixel cuentas o valores de gris, que deben ser calibrados manualmente utilizando información de calibración apropiado.
El paso más crítico en este protocolo es la selección de parámetros adecuados para el análisis de trabecular. En general, grupos de muestras de huesos se analizan utilizando parámetros de segmentación idénticos para hacer los resultados comparables. En tal caso, debe tenerse cuidado para asegurar que resultados de segmentación para todas las muestras se examinan visualmente para posibles errores antes de que otros análisis se llevan a cabo.
Una limitación importante de esta técnica es que se aplican los ajustes de reducción de ruido a huesos corticales sólo para contornear el límite exterior de un hueso, pero no en huesos trabeculares extraídos. Como diversas estrategias de filtrado han sido reportados15 como óptimo para análisis trabecular bajo diferentes condiciones, no solo filtro es buena para todas las muestras de la imagen. Por lo tanto, no hay ningún paso de filtración incorporado antes de la cuantificación de los huesos trabeculares. Para que sea compatible con la salida por software con filtraciones de imagen incorporado, extraídos huesos trabeculares pueden ser procesado de manera similar primer uso varias herramientas de imagen, importen en ImageJ y cuantificaron por el plugin posteriormente. Además, formatos gráficos son generados por diferentes proveedores de micro-CT y la mayoría de ellos no puede ser importada o abierta de ImageJ, directamente, por lo tanto se desean que los módulos para importar diversos formatos de imagen en ImageJ para su posterior análisis. Otra limitación es que los espesores 3D en sectores cerca de los extremos de la gama de selección no son exactos si parte de los objetos está fuera del rango de selección, que también es válido para cualquier paquete de cuantificación 3D disponibles, por lo tanto, una gama más grande de rebanadas debe cuantificarse, y debe utilizarse sólo la parte media de las láminas para su posterior análisis. Este enfoque sólo puede estar disponible cuando cantidades medidas son divulgado capa por capa como se muestra en este plugin, pero no para las mediciones que se divulga solamente un único valor para cada parámetro cuantificado en el rango de selección.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
Se ha presentado una patente para el algoritmo subyacente de cuantificar grosor del objeto. Los autores colaboraron con otros registrados el bomomics.com sitio web hosting los plugins de la versión gratuita, donde consultoría y cuantificación de servicios de análisis de imágenes se proporcionan bajo petición.
Acknowledgments
Este trabajo fue parcialmente financiado por grant 81170806 de NFSC. Los autores desean dar las gracias a la instalación de micro-CT en base de la escuela de Estomatología, Universidad de Wuhan para ayudar a explorar y analizar el fémur de la rata.
Materials
| Name | Company | Catalog Number | Comments |
| ImageJ | NIH | imagej | Any version with a java 1.8 run time |
| trabecular analysis plugin | Bomomics | bomomics | free or commercial version |
| Micro CT scanner | Scanco | μ-50 | micro CT from any vendor |
| Computer System | Lenovo | any brand | |
| Windows Operating System | Microsoft | Windows 7 x64 | any 64-bit Windows operating system |
| Office Software | Microsoft | Office 2010 | any speadsheet software that has xy chart function |
References
- Ruegsegger, P., Koller, B., Muller, R. A microtomographic system for the nondestructive evaluation of bone architecture. Calcif Tissue Int. 58 (1), 24-29 (1996).
- Muller, R., Ruegsegger, P. Micro-tomographic imaging for the nondestructive evaluation of trabecular bone architecture. Stud Health Technol Inform. 40, 61-79 (1997).
- Clark, D. P., Badea, C. T. Micro-CT of rodents: state-of-the-art and future perspectives. Phys Med. 30 (6), 619-634 (2014).
- Bart, Z., Wallace, J. Microcomputed Tomography Applications in Bone and Mineral Research. Advances in Computed Tomography. 2, 121-127 (2013).
- Ji, Y., Ke, Y., Gao, S. Intermittent activation of notch signaling promotes bone formation. Am J Transl Res. 9 (6), 2933-2944 (2017).
- Jiang, Y., Zhao, J., White, D. L., Genant, H. K. Micro CT and Micro MR imaging of 3D architecture of animal skeleton. J Musculoskelet Neuronal Interact. 1 (1), 45-51 (2000).
- Laib, A., et al. 3D micro-computed tomography of trabecular and cortical bone architecture with application to a rat model of immobilisation osteoporosis. Med Biol Eng Comput. 38 (3), 326-332 (2000).
- Cole, H. A., Ichikawa, J., Colvin, D. C., O'Rear, L., Schoenecker, J. G. Quantifying intra-osseous growth of osteosarcoma in a murine model with radiographic analysis. J Orthop Res. 29 (12), 1957-1962 (2011).
- Jensen, M. M., Jorgensen, J. T., Binderup, T., Kjaer, A. Tumor volume in subcutaneous mouse xenografts measured by microCT is more accurate and reproducible than determined by 18F-FDG-microPET or external caliper. BMC Med Imaging. 8, 16 (2008).
- Soviero, V. M., Leal, S. C., Silva, R. C., Azevedo, R. B. Validity of MicroCT for in vitro detection of proximal carious lesions in primary molars. J Dent. 40 (1), 35-40 (2012).
- Kohler, T., Stauber, M., Donahue, L. R., Muller, R. Automated compartmental analysis for high-throughput skeletal phenotyping in femora of genetic mouse models. Bone. 41 (4), 659-667 (2007).
- Buie, H. R., Campbell, G. M., Klinck, R. J., MacNeil, J. A., Boyd, S. K. Automatic segmentation of cortical and trabecular compartments based on a dual threshold technique for in vivo micro-CT bone analysis. Bone. 41 (4), 505-515 (2007).
- Dougherty, G. Quantitative CT in the measurement of bone quantity and bone quality for assessing osteoporosis. Med Eng Phys. 18 (7), 557-568 (1996).
- Doube, M., et al. BoneJ: Free and extensible bone image analysis in ImageJ. Bone. 47 (6), 1076-1079 (2010).
- Bouxsein, M. L., et al. Guidelines for assessment of bone microstructure in rodents using micro-computed tomography. J Bone Miner Res. 25 (7), 1468-1486 (2010).
