Summary
Apresentamos um fluxo de trabalho para segmentar e quantificação dos ossos trabecular para imagens 2D e 3D, com base no limite externo do osso, usando um plugin do ImageJ. Esta abordagem é mais eficiente e precisa do que a atual abordagem manual mão-contorno e fornece quantificações de camada por camada, que não estão disponíveis no software comercial atual.
Abstract
Microtomografia computadorizada (micro-CT) é usada rotineiramente para avaliar a quantidade óssea e trabecular Propriedades microestrutural em pequenos animais sob condições de perda óssea diferentes. No entanto, a abordagem padrão trabecular análise de imagens de micro-CT pedaço por pedaço semi-automática mão-contorno, que é mão de obra intensiva e sujeito a erros. Descrito aqui é um método eficiente para segmentação automática de ossos trabecular de acordo com os limites exteriores do osso, onde os ossos trabecular podem ser identificados e segmentados, automaticamente, com precisão, com menos preconceito do operador quando apropriado parâmetros de segmentação são definidos. Para o perfil de parâmetros de segmentação satisfatória, é exibida uma pilha de imagem dos resultados de segmentação, onde todas as combinações possíveis dos parâmetros de segmentação são alterados um por um, em sequência, e resultados de segmentação com parâmetros associados pode facilmente ser visualmente verificada. Como um recurso de controle de qualidade do plugin, objetos padrão simulados são quantificados onde as quantidades medidas podem ser comparadas com valores teóricos. Quantificação de camada por camada de propriedades trabecular e espessuras trabecular são relatados por este plugin, e as distribuições de tais propriedades dentro as regiões selecionadas podem ser perfiladas facilmente. Embora a quantificação de camada por camada retém mais informações sobre ossos trabecular e facilita ainda mais a análise estatística das mudanças estruturais, essas medidas são indisponíveis desde a saída do atual software comercial, onde apenas um único quantificado o valor para cada parâmetro é relatado para cada amostra. Portanto, os fluxos de trabalho descritos são melhores abordagens para a análise de ossos trabecular com precisão e eficiência.
Introduction
Análise de micro-CT dos ossos trabecular é a abordagem padrão para acompanhamento de alterações morfológicas dos ossos em pequenos animais sob diferentes osso perda condições1,2,3, onde diversas variáveis relacionadas com o estruturas dos ossos são relatados4. No entanto, esses parâmetros não estão distribuídos uniformemente na metáfise de ossos longos5, e somente um valor resumido ou média é relatado para cada variável estrutural de cada amostra pelo atual comercial máquinas de micro-CT6,7 , embora um único valor totalmente não pode representar as características do parâmetro medido na região analisando. Quantificação de camada por camada de ossos trabecular não só retém mais informações para cada variável, mas também permite a criação de perfil das distribuições das variáveis na região analisa, facilitando a posterior análise estatística da estrutural alterações sob diferentes condições de5. Portanto, o objetivo desse método é quantificar ossos trabecular de microtomografia computadorizada em cada nível de fatia, o que não está disponível em qualquer pacote de análise de micro-CT comercialmente disponíveis atualmente.
De forma eficiente segmento ossos trabecular fatia-por-fatia, métodos de segmentação automática são desejáveis. No entanto, a técnica padrão atual para a análise de micro-CT baseia-se no contorno interativo manual seguido por interpolação semi-automática para separar os ossos trabecular dos compartimentos corticais, que é trabalhoso, sujeito a erros, e associado com operador substancial viés8,9,10. Métodos de segmentação automática11,12 foram relatados, mas tais métodos são apenas ideais em regiões com boa separação entre os ossos trabecular e corticais ossos, mas não em regiões sem separações de cores claras. Além disso, são necessários para diferentes amostras de12parâmetros de segmentação diferente, e é tedioso selecionar manualmente os parâmetros de segmentação satisfatória aplicáveis aos grupos de amostras de osso tentando várias combinações parâmetro12, mesmo que o processo de segmentação é automático quando situam-se todos os parâmetros relacionados. Como o limite externo do osso tem o maior contraste com o fundo de digitalização e as conchas cortical metafisária de ossos longos mostram poucas mudanças no escolhido analisando a região, métodos de segmentação de acordo com o contorno limite externo de ossos longos pode confiável e com precisão separe os ossos trabecular conchas corticais. A vantagem de um tal método de segmentação é que a segmentação é baseada na diferença entre o fundo e o limite externo do osso, mas não sobre as diferenças entre os ossos trabecular e cortical6,12, 13, portanto, é geralmente fácil de encontrar uma combinação de parâmetros de segmentação que é satisfatória para um grupo de amostras de osso, facilitando a análise mais confiável de trabecular mudanças entre os diferentes grupos.
Em cada fatia de nível, área, perímetro e bidimensional (2D) de espessura são relatados para análise 2D, enquanto o volume, superfície e espessura de tridimensional (3D) são relatados em 3D quantificações. Tais informações geralmente não são relatadas por ferramentas de análise de imagem atual, indicando que os procedimentos relatados podem ser aplicados a imagens gerais onde tal informação for desejada.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
Procedimentos envolvendo assuntos animais foram conduzidos de acordo com o guia para o cuidado e o uso de animais de laboratório (publicação de NIH, 8ª edição, 2011) e já foi revisados e aprovados pelo cuidado institucional do Animal e uso Comitê de Wuhan Universidade.
1. software instalação
- Instale o software ImageJ. Baixe a versão Windows do software ImageJ (versão 1,51 p) junto com o Java 64-bit de https://imagej.nih.gov/ij/. Extraia o software baixado em uma pasta, que será posteriormente referida como o "diretório ImageJ".
Nota: Plugins Trabecular análise requerem tempo de 64-bit java executar (versão 1.8) e um sistema operacional de 64 bits do windows, de preferência um 64-bit sistema operacional Windows 7. - Instale plugins análise trabecular. Solicitar análise trabecular plugins do http://www.bomomics.com e descompacte o conteúdo para o diretório de plugins do ImageJ, que é "ImageJ diretório/plugins".
Nota: O plugin pode ser obtido como uma versão gratuita, onde as medições para 5 especificado fatias adjacentes são relatadas, ou como uma versão comercial, onde uma série de fatias pode ser especificada e medida.
2. preparar o conjunto de dados 3D para análise Trabecular
- Varredura de fêmures de ratos com um seguinte máquina micro-CT o padrão de varredura protocolo5 e então salvar os dados em um formato que pode ser importado no ImageJ, por exemplo, o formato tiff. Se várias amostras de osso foram escaneadas simultaneamente em um tubo de amostra única, primeiro importar os dados para o ImageJ e em seguida, separe cada osso individual por cortando fora outras amostras usando ferramentas de processamento de imagem do ImageJ. Posteriormente, salve a imagem resultante em um formato que pode ser importado para ImageJ mais tarde.
Nota: Um arquivo de imagem de amostra representativa usado na análise está incluído (arquivo complementar 1). -
Simule o padrão objetos 2D e 3D.
- Abra o software ImageJ. Sob os Plugins | BoMomics | Simular objetos menu, clique com o botão círculo e digite 200 como o diâmetro na janela pop-up, clique Okey para gerar simulado círculo com um diâmetro de 200 pixels (Figura 3B). Salve o círculo gerado no formato tiff.
- Abra o software ImageJ. Sob os Plugins | BoMomics | Simular objetos menu, clique no botão quadrado e digite 200 como o comprimento do lado na janela pop-up, clique Okey para gerar a Praça simulada com um comprimento lateral de 200 pixels (Figura 3B). Salve a Praça gerada no formato tiff.
- Abra o software ImageJ. Sob os Plugins | BoMomics | Simular objetos menu, clique no botão retângulo e digite 200 como a largura e 100 como a altura em janela pop-up, clique Okey para gerar o retângulo simulado com uma largura de 200 pixels e uma altura de 100 pixels (Figura 3 B). Salve o retângulo gerado no formato tiff.
Nota: Um círculo (diâmetro: 200 pixels), um quadrado (comprimento lateral: 200 pixels de) e um retângulo (largura: 200 pixels; Altura: 100 pixels) são salvos para as análises subsequentes. - Abra o software ImageJ. Sob os Plugins | BoMomics | Simular objetos menu, clique no botão do cubo e digite 30 como o comprimento do lado na janela pop-up, clique Okey para gerar o cubo simulado. Finalmente, clique no Plugins | 3D | Visualizador de volume para exibir o cubo gerado e salve-o no formato tiff (Figura 3C).
- Abra o software ImageJ. Sob os Plugins | BoMomics | Simular objetos menu, clique no botão cuboides e digite 80 comprimento 40 como a largura e 30, como a altura em janela pop-up, clique Okey para gerar o cuboide simulado. Finalmente, clique no Plugins | 3D | Visualizador de volume para exibir o cuboide gerado e salve-o no formato tiff (Figura 3C).
- Abra o software ImageJ. Sob os Plugins | BoMomics | Simular objetos menu, clique no botão de esfera e digite 30 como o diâmetro na janela pop-up, clique Okey para gerar a esfera simulada. Finalmente, clique no Plugins | 3D | Visualizador de volume ver a esfera gerada e salve-o no formato tiff (Figura 3C).
- Abra o software ImageJ. Sob os Plugins | BoMomics | Simular objetos menu, clique no botão do cilindro e digite 30 como o diâmetro e 100 como a altura em janela pop-up, clique Okey para gerar o cilindro simulado. Finalmente, clique no Plugins | 3D | Visualizador de volume para ver o cilindro gerado e salve-o no formato tiff (Figura 3C).
Nota: Uma esfera (diâmetro: 30 pixels), um cubo (comprimento lateral: 30 pixels), um cuboide (comprimento: 80 pixels; largura: 40 pixels; Altura: 30 pixels) e um cilindro (diâmetro: 30 pixels; Altura: 100 pixels) são salvos para as análises subsequentes.
3. parâmetros de análise de perfil
- Abra o software ImageJ e abrir ou importar uma imagem digitalizada.
- Deslize a barra de rolagem inferior para escolher uma fatia e, em seguida, clique o imagem | Ajustar | Limite botão. Na janela pop-up de limiar, ajustar os valores de limiar mínimo e máximo por inspeção manual para garantir que os ossos são bem separados do fundo e registre o valor de limiar mínimo como o valor de limite do osso cortical.
Nota: Como fabricantes de micro-CT diferentes usam diferentes fatores de escala para armazenar as imagens os coeficientes de atenuação de raios-x, o valor real limite deve ser determinado empiricamente por inspeção manual ou do fabricante específico a seguir recomendações. Em nossa prática, definir o limite mínimo para arquivos do CEI produzida por máquinas de micro-CT para 6.000-7.000 pode confiantemente separar ossos do fundo varredura. - Clique Plugins | BoMomics | Trab Param de perfil botão. Na janela pop-up, definir Índice de fatia para a posição da fatia representativa e definir o osso cortical ("osso de Cort"), gama ("gama") e valores ("passo") para o cálculo de um conjunto de limites corticais para perfil de parâmetros de segmentação, onde o valor de limiar de osso cortical é adquirido pelo passo 3.2 (Figura 1). Gama e passo valores padrão de trabalho 2.000 e 400 para a maioria das imagens de micro-CT. Definir o Índice de fatia de 5 para o conjunto de dados de exemplo e mantenha os valores padrão para outras configurações.
Nota: "Osso de Cort" é o limite para ossos corticais e o valor do limite utilizado para alterações de perfil parâmetro do valor limiar mais baixo para o valor de limiar mais alto com o incremento do valor do passo, onde o valor de limiar mais baixo é o valor de Cort Osso - gama e o maior valor de limiar é osso Cort + intervalo. Se os valores definidos para osso Cort, gama e passo são 6.000, 1.000 e 500, respectivamente, em seguida, o menor valor de limiar para ossos corticais é 6.000-1, 000 = 5, 000, valor limite máximo é de 6, 000 + 1, 000 = 7, 000 e os limiares utilizados em análises são 5.000, 5.500, 6.000, 6.500 e 7.000. - Defina o diâmetro de ruído ("diâmetro de ruído"), passo ("") e faixa de valores ("Range") para especificar um conjunto de valores de ruído nas análises e o diâmetro do furo ("diâmetro de furo") ("passo") e valores de gama ("gama") para o cálculo de um conjunto de valores de buraco. Em geral, a configuração padrão funciona para a maioria das amostras de osso de micro-CT, onde ruído, passo e gama são 5, 5, 2, respectivamente, e gama, passo e diâmetro do furo são 15, 5, 2, respectivamente. Mantenha as configurações padrão para o conjunto de dados de exemplo.
Nota: "Diâmetro de ruído" é o diâmetro do filtro de supressão de ruído, e o diâmetro de"buraco" é o diâmetro para buracos dentro dos ossos corticais. Os valores de ruído e o buraco usados para perfilação de parâmetro podem ser calculados da mesma forma como valores de limiar cortical descritos acima, usando o especificado ruído/buraco, gama e passo valores. Devido as variações nos parâmetros de digitalização aquando da aquisição de imagens, parâmetros de ruído e o buraco devem ser determinados empiricamente de acordo com a qualidade de imagem, como sem intervalos gerais dos parâmetros de ruído e buraco são bons para todos os exames. O menor valor de limiar, ruído ou furos deve ser maior ou igual a 0, e se o menor valor calculado é negativo usando os parâmetros fornecidos, o menor valor específico é definido como 0. - Clique em Okey para executar o parâmetro de criação de perfil. Visualmente a verificar os resultados da segmentação na janela de Resultados de criação de perfil de parâmetro e selecione uma camada de fatia, para que o limite externo do osso é delineado com bastante precisão (Figura 1B). Posteriormente, recupere os parâmetros de perfil a partir da entrada na tabela de Resultados de criação de perfil de parâmetro correspondente à camada de fatia escolhido (tabela 1).
4. trabecular análise
-
Segmentação dos ossos trabecular
- Abra o software ImageJ, em seguida, abrir ou importar uma imagem digitalizada.
- Clique Plugins | BoMomics | Segmentação de trab botão e preencher os parâmetros de análise adequada. Definir o "Iniciar", "Limite de contorno", "Trab. Ossos","Diâmetro de redução de ruído","Buraco de enchimento diâmetro"e"Diâmetro de espessura Cortical"para 5, 7.200, 7.000, 6, 12 e 25, respectivamente.
Nota: "Start" e "End" especifica o intervalo de fatia selecionada para segmentação dos ossos trabecular, "Limite de estrutura de tópicos" é correspondente ao parâmetro perfilado de "Cort-osso", "Diâmetro de redução de ruído" para "Diâmetro de ruído" Parâmetro e "Buraco preenchendo dia." para o parâmetro "Diâmetro de furo". "Espessura cortical diâmetro" é a espessura especificada para excluindo os ossos corticais externas. "Trab. Ossos "são o limite para a extração de ossos trabecular (Figura 2), onde as configurações são determinadas usando o comando de parâmetros de análise de perfil, conforme descrito na etapa 3.5. - Clique em Okey para realizar a segmentação trabecular. Verificar visualmente os resultados de segmentação na janela de Resultados de segmentação de Trab (Figura 2B). Salve os ossos trabecular extraídos mostrados na janela de Ossos Trabecular segmentado em formato tiff (Figura 2B), que pode ser mais analisada por outro software.
-
Análise dos ossos trabecular.
- Abra o software ImageJ, em seguida, abrir ou importar uma imagem digitalizada.
- Clique Plugins | BoMomics | Análise de trab botão e preencher os parâmetros de análise adequadas, tais como "Start", "Limite de contorno", "Trab. Ossos","Diâmetro de redução de ruído","Diâmetro de furo de enchimento"e"Diâmetro de espessura Cortical"como descreveram acima (Figura 3A), onde as configurações são determinadas usando o comando de parâmetros de análise de perfil, como descrito no passo 3.5. Definir o "Iniciar", "Limite de contorno", "Trab. Ossos","Diâmetro de redução de ruído","Buraco de enchimento diâmetro"e"Diâmetro de espessura Cortical"para 5, 7.200, 7.000, 6, 12 e 25, respectivamente.
Nota: Na versão gratuita do plugin, cinco fatias adjacentes, começando o índice especificado do fatia "Start" são selecionadas para medição, enquanto a versão comercial, um número arbitrário de fatias pode ser especificado pelo usuário. - Selecione uma ou mais opções na seção de Relatórios de resultados para os parâmetros a medir, onde o volume ósseo trabecular (BV), o volume total (TV) da região selecionada e espessuras medido qualquer forma bidimensional (2D) ou tridimensional (3D) estão disponíveis para seleção através de três caixas de seleção, ou seja "BV TV só", "2D" e "3D". Selecione as caixas de seleção "2D" e "3D", clique em "Okey" para executar a análise trabecular (Figura 3, tabela 2).
Nota: Quando "TV BV só" estiver marcada, independentemente do status de seleção de "2D" e "3D", crus medidas de BV, TV e intensidade são relatadas, e os ossos trabecular segmentados são extraídos e exibidos em uma nova janela, que pode ser salvos e mais analisada por outro software. Quando "2D" estiver marcada, crus medidas de BV, TV, intensidade e espessuras medidas de forma bidimensional em cada fatia de nível usando o modelo de placa são relatadas. Se "3D" estiver marcada, espessura tridimensional para cada voxel é calculada diretamente, sem qualquer hipótese de modelo, então crus medidas de BV, TV, intensidade e espessuras tridimensionais amostradas em cada nível da fatia são relatadas. No entanto, caso não esteja selecionada, ossos trabecular são segmentados usando o conjunto acima dos parâmetros, mas nenhuma medição é relatada.
5. quantificação dos objetos simulados
- Abra o software ImageJ e, em seguida, abrir uma imagem simulada. Aqui, abra a simulado esfera com um diâmetro de 30 pixels como exemplo.
- Selecione os Plugins | BoMomics | Análise de trab botão e preencher os parâmetros de análise apropriado conforme descrito anteriormente. Mantenha os valores padrão para "Iniciar", "Fim", "Limite de contorno", "Trab. Ossos"e definir"Diâmetro de redução de ruído","Diâmetro de redução de ruído","Buraco de enchimento diâmetro"e"Diâmetro de espessura Cortical"para 0 (Figura 3).
Nota: Para objetos simulados, não há nenhum ruído de sinal não intencionais, e pode haver sem conchas corticais correspondentes. Portanto, parâmetros para tais valores devem ser definidos apropriadamente (configurações padrão são a zero). Para 3D quantificações, 30 fatias de antes e depois a fatia de "Start" especificada são processadas usando a versão gratuita do plugin, enquanto apenas fatias no intervalo especificado pelo "Start" e "Fim" fatias são analisadas pela versão comercial do plugin. - Na seção de Relatórios de resultados , selecione "2D" e "3D" para os parâmetros a medir e clicar Okey para executar análise trabecular para o objeto simulado (tabela 3).
6. calibração de medidas Trabecular e apresentação de dados: perfil as distribuições das medidas Trabecular na região selecionada analisando
- Obter as informações de calibração do conjunto de dados digitalizado, de acordo com instruções do fornecedor micro-CT.
Nota: Apenas crus medidas de BV, TV, intensidade e espessuras são relatadas pelo plugin. Para obter medidas comparáveis aos relatórios gerados por outros softwares, calibração de resultado é desejada. - Abra o Microsoft Excel e abra a tabela de resultados relatados. Calcular calibrado volume ósseo (BV), total volume (TV), conteúdo mineral ósseo (BMC), fração de volume de osso (BV/TV) e densidade mineral óssea (DMO) no novo excel colunas usando os valores brutos relatados BV, TV e intensidade de acordo com as seguintes expressões.
- Nota: Digitalização resolução (resolução, µm), escala de cinza (escala), unidades de densidade (mg HA/cm3), inclinação de densidade (inclinação) e interceptação de densidade (intercepção) podem ser extraídos com os micro-CT imagem ou meta arquivos digitalizados, tais como o arquivo Scanco ISQ. Portanto, as medidas calibradas são calculadas como segue:
Suponha que BV, TV e intensidade são medidas crus, BVc e TVc são valores calibrados, resolução, escala, inclinação, intercepção são arquivos de calibração de micro - CT imagens.
BVc = BV × resolução3 [μm3]
TVc = resolução de TV ×3 [μm3]
BV / TV = BV ÷ TV
BMC = (intensidade ÷ dimensionamento × inclinação - interceptação de × BV) resolução ×3 × 10-12 [mg HA]
BMD - BMC ÷ TVc × 1012 [mg HA / cm3] - Crie uma trama de dispersão (XY) para medidas calibradas (Y) contra camadas de fatia (X) utilizando o software Microsoft Excel (Figura 4).
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
O plugin de análise trabecular é projetado para automaticamente segmentar e ossos trabecular de quantificar com precisão. Inicialmente, osso limite externo é detectado e delineado seguido por uma operação de preenchimento de buracos onde qualquer buracos dentro de conchas cortical externa do osso são preenchidos. Então, uma operação de erosão é realizada para excluir os ossos corticais externas e pegar os ossos trabecular segmentados. Finalmente, as medidas dos ossos trabecular na região segmentado são quantificadas.
Como micro-CT imagens são inerentemente barulhentas, segmentações usando parâmetros predefinidos arbitrários, muitas vezes não conseguem identificar com precisão os limites exteriores do osso. É tedioso e do trabalho intensivo para tentar um monte de combinações de parâmetros para a seleção de parâmetros de segmentação satisfatória. Portanto, um parâmetro de perfil plugin é fornecido para alterar os parâmetros de um por um no intervalo definido automaticamente para ajudar na seleção de combinações de parâmetros satisfatórios, que também facilita a seleção de um conjunto comum de parâmetros para um grupo de osso amostras. Figura 1 A mostra as configurações usadas para criação de perfil parâmetros de segmentação bom. Quando os parâmetros para o limite de osso cortical (osso Cort), intervalo do limiar para ser perfilado (gama), e o valor de incremento para o ponto inicial (etapa) em cada etapa são especificados, uma série de limites a ser perfilado são gerados. Posteriormente, uma série de valores de ruído e buraco são gerados do mesmo modo, definindo os parâmetros correspondentes. Finalmente, ossos trabecular são segmentados, alterando os parâmetros de um de cada vez para todas as combinações possíveis de parâmetro. Figura 1 B mostra os resultados representativos de segmentação para combinações de parâmetros diferentes. Obviamente, algumas combinações de parâmetro são melhores que outros em delinear limites exteriores do osso, e mais do que uma combinação de parâmetro mostra resultados satisfatórios de segmentação. Depois de verificar visualmente os resultados da segmentação, os valores de parâmetro para segmentações satisfatórias podem ser recuperados da tabela de resultados de perfil (tabela 1).
Para quantificar as medidas dos ossos trabecular, análise e segmentação trabecular são executadas. Figura 2 A mostra a caixa de diálogo configuração de segmentação trabecular, onde ossos trabecular na região selecionada são segmentados e extraídos (Figura 2B), e os resultados de segmentação usando os parâmetros fornecidos podem ser verificados, pedaço por pedaço visualmente. Posteriormente, os ossos trabecular são analisados com parâmetros satisfatórios(Figura 3). Dependendo do status de seleção de opções, crus quantificações de volume ósseo (BV), emissão de relatórios total volume (TV), a soma dos valores de cinza (intensidade), e espessuras medido duas dimensões ou três dimensionalmente (tabela 2) são relatados. Finalmente, informações de calibração são extraídas do conjunto de dados digitalizado micro-CT e calibrado medidas de BV, BV/TV TV, BMC e BMD são calculados seguido por suas distribuições na selecionado analisando região camada por camada contra a camada de criação de perfil posições (Figura 4).
Como um recurso de controle de qualidade do plugin, quantificação de objetos simulados é suportada. Simulado padrão objetos com dimensões conhecidas são quantificados pelo plugin para comparação com os valores teóricos ou medidas de outros softwares, tais como software comercial fornecido com máquinas de micro-CT ou BoneJ14, um fonte aberta livre plugin para a análise de imagem do osso. Para análise de objetos simulados, as configurações de ruído são definidas como zero, como imagens simuladas são consideradas imagens de alta qualidade, sem qualquer ruído. Objetos com várias espessuras foram simulados e os resultados são apresentados (Figura 3, tabela 3). Para objetos 2D padrão simulados, tais como círculos, quadrados e retângulos, exatos valores para área (TV ou BV) e espessura são relatados (tabela 3). Para objetos 3D, medidas de espessura exata para cubos, esferas e cuboides são relatadas, no entanto, as espessuras para cilindros não são precisas para voxels perto de ambas as extremidades do cilindro, enquanto as espessuras de voxel nas fatias médios dos cilindros são exatamente como previu. Esta é uma característica do algoritmo subjacente de medição espessura, onde a espessura do voxel para cada objeto é determinada pelo diâmetro da esfera maior, ou o comprimento do lado do cubo maior, que contém este voxel e é completamente dentro do objeto. Portanto, espessuras de objetos feitos de diferentes esferas e cubos podem ser medidos com precisão, enquanto cilindros só podem ser medidos com precisão, a fatia meio raio-distância longe de ambas as extremidades.
Figura 1 : Resultados representativos de parâmetro de análise de perfil. (A) os parâmetros de configuração de página. (B) resultados representativos de parâmetro de análise de perfil. Algumas combinações de parâmetro são melhores do que outros para detectar os limites exteriores dos ossos. Clique aqui para ver uma versão maior desta figura.
Figura 2 : Resultados representativos da análise de segmentação trabecular. (A) os parâmetros de configuração de página. (B) resultados de segmentação representativa dos ossos trabecular em diferentes camadas. Clique aqui para ver uma versão maior desta figura.
Figura 3 : Análise trabecular. (A) os parâmetros de configuração de página. (B) resultados representativos de objetos 2D simulados. (C) resultados representativos de objetos 3D simulados visualizados pelo Visualizador de volume 3D. Clique aqui para ver uma versão maior desta figura.
Figura 4 : Distribuições das medidas de osso trabecular na região analisa selecionada. O eixo horizontal representa a distância relativa à camada de fatia de partida na região analisa. Valores do eixo y são calibradas trabecular medidas na região analisa. Clique aqui para ver uma versão maior desta figura.
| Fatia | Limite de estrutura de tópicos | Diâmetro de redução de ruído | Diâmetro de furo enchimento |
| 4 | 5200 | 0 | 10 |
| 147 | 5200 | 2 | 20 |
| 361 | 7200 | 6 | 16 |
| 539 | 8000 | 10 | 20 |
Tabela 1: Resultados representativos de parâmetro de análise de perfil.
Tabela 2: Resultados representativos de análise trabecular.
| Objeto | Dimensão | Volumeb | Superfíciec | Espessura |
| Praça | 200 X 200 | 40000 | 796 | 200 |
| Retângulo | 200 X 100 | 20000 | 596 | 100 |
| Círculo | Diâmetro: 200 | 31428 | 796 | 200 |
| Cubo | 30 X 30 X 30 | 27000 | 5048 | 30 |
| Cuboide | 80 X 40 X 30 | 96000 | 13008 | 30 |
| Esfera de | Diâmetro: 30 | 14328 | 3944 | 30 |
| Cilindrod | Dia:30; H: 100 | 71600 | 12800 | 27,84 |
| Cilindroe | Dia:30; H: 100 | 51552 | 9552 | 30 |
| r: os resultados são em voxels cru. Diâmetro: diâmetro; H: altura. | ||||
| b: Volume (3D) ou área (2D). | ||||
| c: superfície (3D) ou perímetro (2D). | ||||
| d: Measurs de fatia 1 fatia 100 são utilizados para análise. | ||||
| e: Measurs de fatia 15 a fatia 85 são utilizados para análise. | ||||
Tabela 3: Resultados de quantificação de objetos simulados.
Complementar arquivo 1. Amostra de osso. Clique aqui para baixar este arquivo.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
Este estudo descreve um ImageJ plugin para análise de ossos trabecular, que é automático, eficiente, e amigável. O plugin também pode ser usado para quantificar qualquer objecto 2D ou 3D para medidas de camada por camada de áreas, volumes e espessuras. Atualmente, apenas um único valor medido para cada parâmetro trabecular é relatado para cada amostra por análise de micro-CT padrão, que não representa totalmente as características da entidade medida na região selecionada analisa. O plugin descrito relata quantidades de camada por camada de cada parâmetro para cada amostra, que totalmente retém as informações de distribuição do parâmetro medido na região analisa selecionada e, portanto, mais avançadas e sensíveis à estatística abordagens são aplicáveis para análise de tais dados.
Com o padrão, semiautomática mão-segmentação método, observaram-se uma variação substancial entre diferentes fatias mão-contorno com operador viés8,9,10, justificando um uniforme mais e o método automático de segmentação trabecular. Como o contraste é máximo entre o digitalização do fundo e o limite externo do osso, limite externo assistida segmentação trabecular é executada pelo plugin, onde limite externo do osso pode ser detectado automaticamente com precisão quando apropriado limiar, o ruído e as configurações do furo são fornecidas. Para facilitar a determinação de parâmetros de segmentação, uma gama de parâmetro são combinações perfilado usando uma fatia representativa da região analisa um por um, e resultados satisfatórios de segmentação são verificados visualmente. Posteriormente, os ossos trabecular são analisados com os parâmetros perfilados. Quando são definidos parâmetros de segmentação apropriada, e o plugin tem sido aplicado com sucesso para analisar ossos trabecular de fêmures distal de ratos tratados com osso diferentes agentes anabólicos5, resultados reprodutíveis são relatados para a mesma imagem sem viés de operador.
O plugin atualmente só pode processar uma amostra de cada vez. Se várias amostras estão presentes na imagem, todas as amostras são quantificadas como um único objeto sem discriminação. Portanto, pré-processamento é necessária se múltiplas amostras são verificadas simultaneamente em um tubo de amostra única. As saídas de quantificação são crus medidas de contagem de pixels ou valores de cinza, que devem ser calibrados manualmente usando informações de calibração adequados.
O passo mais crítico neste protocolo está selecionando os parâmetros adequados para análise trabecular. Geralmente, grupos de amostras de osso são analisados usando parâmetros idênticos segmentação para que os resultados sejam comparáveis. Nesse caso, deve ter cuidado para garantir que resultados de segmentação para todas as amostras são examinados visualmente por possíveis erros antes de mais análises são conduzidas.
Uma grande limitação dessa técnica é que as configurações de redução de ruído são aplicadas para ossos corticais, apenas para o limite externo de um osso de contorno, mas não em extraídos ossos trabecular. Como várias estratégias de filtragem foram relatados15 como ideal para análise trabecular sob diferentes condições, sem filtro único é bom para todas as amostras de imagem. Portanto, não há nenhuma etapa de filtragem interna antes de quantificação dos ossos trabecular. Para ser compatível com a saída de software com filtrações de imagem interna, extraído de ossos trabecular podem ser processados da mesma forma primeiro usando várias ferramentas de imagem e então importados ImageJ e quantificada pelo plugin posteriormente. Além disso, diferentes formatos de imagem são gerados por diferentes fornecedores de micro-CT e a maioria deles não pode ser importada ou aberta pelo ImageJ diretamente, assim, módulos para importar vários formatos de imagem em ImageJ para posterior análise são desejados. Outra limitação é que o 3D espessuras em fatias perto as extremidades da faixa de seleção não são precisas se parte dos objetos que estão fora do intervalo de seleção, que também é verdade para qualquer pacote de quantificação 3D disponível, portanto, uma maior gama de fatias devem ser quantificadas, e somente a parte média de fatias selecionadas deve ser usada para posterior análise. Tal abordagem só pode ser disponível quando quantidades medidas são relatada camada por camada, como mostrado neste plugin, mas não para as medições que apenas um único valor é relatado para cada parâmetro quantificado na faixa de seleção.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
Uma patente para o algoritmo subjacente de quantificar a espessura do objeto. Os autores colaboraram com outros registrado o bomomics.com site de hospedagem os plugins da versão gratuita, onde consultoria e quantificação dos serviços de análise de imagem são fornecidos mediante solicitação.
Acknowledgments
Este trabalho foi parcialmente financiado por grant 81170806 NFSC. Os autores gostaria de agradecer a facilidade do núcleo de micro-CT da escola de Estomatologia, Universidade de Wuhan para ajudar a digitalizar e analisar os fêmures de ratos.
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).
