Radiologia Brasileira - Publicação Científica Oficial do Colégio Brasileiro de Radiologia

 

 

INTRODUÇÃO

O ser humano é uma espécie biológica dinâmica que sob ação de processos metabólicos sucessivos evolui, segundo os fatores individuais e ambientais, devido aos mecanismos de crescimento e desenvolvimento. O que se observa é que o corpo humano tem em sua formação e crescimento sequências e regras próprias, as quais vêm sendo observadas e estudadas sob vários aspectos⁽¹⁾.

Uma das formas de se estudar o crescimento é por meio da análise de radiografias de determinadas partes do corpo, e entre as mais utilizadas está a mão, por ser um local de fácil acesso e por não conter órgãos vitais, podendo ser submetida a doses de radiação com mais frequência⁽¹⁾.

Os métodos mais difundidos e utilizados são: Greulich e Pyle (GP), que faz uma avaliação inspecional dos ossos da mão através de um Atlas^(1,2); o método de Tanner e Whitehouse (TW), que analisa 20 ossos da mão e punho atribuindo um escore específico para cada osso e sexo, e pela soma desses escores obtém-se a idade óssea^(3-5) ; e o método de Eklof e Ringertz (ER), o qual se baseia em medidas de dimensão de dez centros de ossificação, formados por ossos do carpo, largura da epífise distal do rádio, comprimento e largura do capitato e do hamato, comprimento dos metacarpais II, III e IV, e comprimento das falanges proximais II e III^(2,4). Este método apresenta duas tabelas, uma para cada sexo, nas quais cada centro de ossificação possui um padrão mínimo e um máximo de medida. A idade óssea é calculada através de um valor médio das idades obtidas para cada centro de ossificação. Descrição mais detalhada desses métodos e dos processos de estimação pode ser obtida no portal de informações de idade óssea (www.carpal.eesc.usp.br).

Encontram-se alguns trabalhos nesta linha de pesquisa, dentre eles o de Taffarel et al.⁽⁵⁾, que desenvolveram um processo de estimação da idade óssea através da automatização do método de TW. O trabalho de Haiter et al.⁽⁶⁾ relata um estudo com o objetivo de verificar se os métodos de GP, TW e ER poderiam ser aplicados à população brasileira e qual destes métodos seria o mais confiável, quando comparados à idade cronológica do indivíduo.

De posse das características de cada um dos métodos de estimação citados, o de ER foi o explorado e implementado neste trabalho, por se tratar de um método baseado em medidas de dimensão, podendo ser implementado através de processamentos computacionais. Apresentam-se, neste trabalho, processos simplificados de estimação da maturidade óssea derivados de ER - denominados ER5 e ER3 -, que operam de forma automatizada.

 

MATERIAIS E MÉTODOS

Conjunto de imagens do estudo

Foram selecionadas duas bases de imagens radiográficas contituídas por pacientes na faixa etária de 6 a 16 anos. A primeira é formada por 685 pacientes, sendo 357 do sexo masculino e 328 do sexo feminino (Figura 1). As radiografias foram obtidas junto ao Departamento de Odontologia da Faculdade de Odontologia de Piracicaba - Universidade Estadual de Campinas, sendo classificadas pelos métodos de GP e TW, sob supervisão do Prof. Francisco Haiter Neto, e classificadas por ER com auxílio do software Anacarp⁽⁷⁾. A segunda base de imagens radiográficas é formada por 234 pacientes, sendo 130 do sexo masculino e 104 do sexo feminino. As imagens são oriundas da região de Bauru, SP, e foram cedidas pelo Dr. Orivaldo Tavano⁽⁴⁾. Esta segunda base de imagens foi incluída no estudo para que fosse possível a investigação de resultados considerados discrepantes em relação à média dos laudos médicos (ER, GP e TW) obtidos pelos métodos quando aplicados na primeira base de imagens.

 

 

Métodos

Para a implementação do método foi utilizada a plataforma Borland-Builder C++ e as análises estatísticas foram realizadas por meio do software BioEstat 3.0⁽⁸⁾.

1 -Propostas de simplificações para o método de ER

Uma das dificuldades encontradas no método ER está relacionada com a sobreposição dos ossos do carpo e punho, causando dificuldades no processo de obtenção de suas medidas, o que prejudica a obtenção precisa das medidas dos centros de interesse e pode, por consequência, produzir resultados pouco confiáveis. Para sanar este problema, buscou-se excluir estes ossos do processo de estimação, propondo-se, dessa forma, uma simplificação do método ER, denominada ER5.

a) A simplificação ER5

Nesta simplificação, a estimação da idade óssea é realizada com base na análise de apenas cinco centros de ossificação (ER5), compostos por duas falanges proximais e três ossos metacarpais, descartando os ossos do carpo e punho (Figura 2).

 

 

Posteriormente, realizou-se estudo estatístico objetivando uma nova simplificação para o método ER, sendo assim elaborado o método ER3.

b) A simplificação ER3

Buscou-se verificar o limite mínimo do número de centros ósseos que apresentassem resultados estatisticamente significantes, ou seja, próximos ao valor da média dos laudos médicos. Foram feitas todas as combinações possíveis entre os cinco ossos (tomados 2 a 2, 3 a 3 e média dos cinco ossos) utilizados no ER5, comparando-as com a média dos laudos médicos. A denominação para os ossos utilizada nas combinações foi: P1, P2 e P3 - comprimento dos metacarpais II, III e IV; P4 e P5 - comprimento das falanges proximais II e III. Em seguida, foram gerados gráficos para cada imagem do conjunto (Figura 3), analisando quais combinações mais se aproximaram do valor do laudo médico.

 

 

A linha tracejada (com quadrados) no gráfico da Figura 3 refere-se ao valor médio dos laudos médicos (ER, GP e TW). Os pontos selecionados no gráfico foram as combinações P1 P5 (círculo com rótulo 1) e P1 P3 P5 (círculo com rótulo 2). A média da combinação P1 P3 P5, formada pelos ossos em destaque na Figura 4, apresentou excelente aproximação em relação ao valor médio dos laudos médicos, indicando uma simplificação para o método ER.

 

 

2 -Método de estimação da idade óssea

Na Figura 5 são apresentados os passos para estimação da idade óssea.

 

 

1. Pré-processamento - No sentido de eliminar as variações da não uniformidade encontrada no fundo da imagem, inerente ao processo de obtenção das radiografias, foi aplicado um algoritmo desenvolvido por Nascimento et al.⁽⁹⁾ e exposto no trabalho de Olivete et al.⁽¹⁰⁾ para correção do efeito heel e um filtro passa-baixa para suavizar os ruídos⁽¹¹⁾.

2. Segmentação dos ossos - Uma forma simples de se isolar os ossos da mão pode ser conseguida com a utilização de um método eficiente de thresholding⁽¹¹⁾. Os métodos implementados foram: Otsu, Pun, Niblack e Rosenfeld^(12-14) .

3. Isolamento dos centros de ossificação - Para selecionar os dedos da mão, construiu-se um procedimento que realiza varreduras horizontais na imagem, guiando-se pelo pixel com intensidade 255 (cor do fundo). Posteriormente, inseriram-se marcadores para isolar os ossos metacarpais II, III e IV e as falanges proximais II e III.

4. Medidas dos centros de ossificação - Após a inserção dos marcadores, é calculado o comprimento, em pixels, de cada osso, usando a distância euclidiana. Para ser utilizado pelo método ER, esse valor é convertido para milímetros.

5. Estimação da idade óssea - Para cada medida encontrada, consulta-se a tabela de medidas para identificação do valor referente à idade. A idade final é obtida pela média da idade encontrada para todos os centros de ossificação.

3 -Processo de análise dos dados

As idades ósseas foram estimadas utilizando a ferramenta desenvolvida Anacarp⁽⁷⁾; posteriormente, foram analisadas utilizando o software BioEstat 3.0. Calcularam-se os valores dos desvios-padrão, média, coeficiente de correlação e teste t de Student para os valores obtidos pela média dos laudos (TW, GP e ER) e métodos simplificados (ER5 e ER3). Com o objetivo de verificar a correlação entre a idade estimada, utilizando ER5 e ER3, e a média dos laudos médicos, utilizou-se a correlação de Pearson. Por esta análise obtém-se um modelo de regressão linear do tipo Y = aX + b, onde X representa a média dos laudos, Y representa a idade óssea estimada, e a e b são coeficientes desta regressão. Para medir a taxa de concordância com a média dos laudos médicos, o coeficiente a deve ser igual a 1 e o coeficiente b deve ser igual a zero; assim, Y = X. No sentido de verificar se esta relação é estatisticamente significativa, aplicou-se o teste t de Student pareado para o coeficiente de correlação, com nível de significância de 5%.

 

RESULTADOS

Para melhorar a apresentação e análise dos resultados, as imagens foram classificadas de acordo com a média dos laudos médicos considerando grupos de faixa etária de 6 a 16 anos, em que cada faixa etária abrange uma quantidade de indivíduos; por exemplo, 6 anos compreende indivíduos de 6 anos até 6 anos e 11 meses. Na Tabela 1 são apresentados os resultados para os indivíduos do sexo masculino.

Analisando a Tabela 1, nota-se que a estimação realizada pelo método ER5, sexo masculino, obteve correlação forte e positiva com as médias dos laudos médicos, variando entre 0,58 e 0,96. Observa-se que somente os resultados obtidos para os indivíduos da faixa etária 11 anos não foram estatisticamente significantes, com valor p (teste t de Student) igual a 0,4512.

A fim de se identificar se o método simplificado ER5 não poderia ser utilizado para estimação da maturidade óssea para a faixa etária de 11 anos (sexo masculino), selecionaram-se 32 indivíduos de outro conjunto de imagens, oriundas de outra região brasileira, e aplicou-se o método. Os resultados obtidos foram estatisticamente significantes (valor de p [t de Student pareado] = 0,0076, com correlação de 72%), o que comprovou que a simplificação pode ser aplicada à faixa etária de 11 anos. Para compreender a divergência encontrada no primeiro banco, foram analisadas, individualmente, as imagens nas quais os resultados não foram significantes, e das 44 imagens, se for considerado um desvio-padrão de 6 meses - um valor aceitável entre os padrões apresentados pelos métodos clássicos, que variam entre 6 meses e 1 ano, dependendo da faixa etária -, apenas três indivíduos (7%) permaneceram com valores estatisticamente não significantes. Com isso, pode-se inferir que na produção dos laudos para a faixa de 11 anos do primeiro conjunto, com dois dos métodos clássicos (GP e TW) que são realizados por procedimentos de comparação, pode ter havido imprecisão suficiente para a geração de valores distintos daqueles esperados para os laudos.

No conjunto de imagens do sexo masculino, com o método ER5, os resultados mostraram fortes correlações, com valores chegando a 97% e dentro da faixa de significância (p < 0,05). Os resultados não foram significantes apenas para os grupos com faixas etárias de 9 e 14 anos. Avaliando os resultados do grupo de indivíduos com 9 anos (total de 28 imagens), duas imagens (7,14%) tiveram idades estimadas que diferiram 7 e 8 meses, respectivamente, do valor esperado (média dos laudos). No grupo de 14 anos, cerca de 10% dos laudos gerados ficaram fora do intervalo de desvio-padrão de 6 meses em relação à média dos laudos. Aqui também cabe o mesmo comentário feito para a faixa etária de 11 anos, ou seja, há grande possibilidade de erro na emissão dos laudos com dois dos métodos clássicos (GP e TW), que por serem realizados por procedimentos de comparação estão sujeitos às subjetividades de interpretação daqueles que emitiram os laudos.

Na Figura 6 é apresentado um gráfico ilustrando as correlações entre as médias dos laudos médicos (TW, GP e ER) e as simplificações (ER5 e ER3), considerando as faixas etárias e o sexo masculino. É possível notar que as simplificações obtiveram forte correlação (valores acima de 0,5) quase na totalidade das faixas etárias.

 

 

Uma análise similar foi realizada para os indivíduos do sexo feminino e os resultados estão apresentados na Tabela 2.

Analisando a Tabela 2, verifica-se que os resultados mostram correlações positivas e elevadas com a média dos laudos médicos, com valores chegando a aproximadamente 91%. O pior resultado foi o do grupo de 14 anos, com uma correlação de 0,5317 e valor de p = 0,2830 (acima do nível de significância de 5%). Neste grupo foram avaliadas 31 imagens e em 5 casos os laudos estimados (variando entre 7 meses e 1 ano e 2 meses) diferiram do valor do laudo médio, considerando um desvio-padrão de 6 meses.

Na Tabela 2 também se observam resultados não significantes para as faixas etárias de 9 e 10 anos (p > 0,05). Analisando individualmente os laudos para estes grupos etários, constatou-se que para a faixa de 9 anos, 7,5% dos laudos estimados diferiram do valor do laudo médio. Já para a faixa etária de 10 anos, esse erro foi de 3,7%. Para verificar essas diferenças estatísticas dentro dessas faixas etárias, repetiu-se a estimação para os indivíduos pertencentes ao conjunto de imagens 2 (Database II) - mesmas faixas etárias -, e os resultados avaliados a partir da correlação de Pearson ficaram acima de 79% e valores de p (t de Student) dentro da faixa de significância (< 0,003).

 

DISCUSSÃO

Os métodos que tradicionalmente têm sido utilizados como referências confiáveis no processo de estimação da idade óssea são baseados na análise da radiografia carpal: GP, TW e ER. Os dois primeiros necessitam de uma análise inspecional comparativa, inserindo significante dose de subjetividade no processo. O método de ER é um dos mais utilizados em clínicas especializadas que utilizam ferramentas computadorizadas para auxílio ao diagnóstico. Este método apresenta algumas dificuldades de utilização, dentre elas o processo manual de obtenção das medidas dos centros de ossificação, em que é necessária a inserção de 20 marcadores para a identificação de todos os ossos utilizados no processo^(7,15).

Optou-se por explorar o método de ER por ser um método padrão bastante usado na avaliação da idade óssea no âmbito odontológico⁽¹⁶⁾. Utilizando-se parâmetros originais do método, desenvolveu-se uma ferramenta - Anacarp⁽⁷⁾ - que realiza o processo de estimação de forma automática e simplificada, diferentemente dos softwares encontrados na literatura, que operam de forma manual^(6,16). Estes necessitam de um elevado tempo do operador, pois são necessários 20 marcadores para identificar todos os centros de ossificação. Na ferramenta Anacarp este processo de identificação é feito de forma automática e simplificada, oferecendo resultados de elevada concordância com o valor médio dos laudos médicos (obtidos por GP, TW e ER).

Os resultados foram analisados pela correlação de Pearson e com o teste t de Student pareado. Resultados iniciais não foram estatisticamente significantes (valores de p acima de 0,05) quando utilizado o método ER5 para indivíduos da faixa etária de 11 anos (sexo masculino) e 14 anos (sexo feminino), o que motivou a utilização de um segundo banco de imagens, o que permitiu inferir que os laudos emitidos com GP e TW para o primeiro banco estão com valores distintos daqueles esperados para os laudos. No método ER3, a mesma observação pode ser feita para os indivíduos incluídos nas faixas etárias de 9 e 14 anos (sexo masculino) e 9 e 10 anos (sexo feminino). Repetiram-se as análises para estas faixas etárias utilizando indivíduos de outra base de imagens e os laudos gerados por ambos os métodos e sexos foram coerentes com a média dos laudos médicos, considerando a faixa de significância de 5%. Foram feitas análises individuais nos grupos de imagens em que os laudos estimados ficaram fora da faixa de significância, tomando-se como referência um desvio-padrão de 6 meses em relação à média dos laudos, e o erro obtido em ambos os sexos variou entre 3,7% e 10%.

 

CONCLUSÕES

Com base nos resultados obtidos nos processos de estimação usando as simplificações (ER5 e ER3), em ambos os sexos e para todas as faixas etárias, observou-se correlações positivas e elevadas com os valores das médias dos laudos médicos, podendo ser observadas no item "Correlação" das Tabelas 1 e 2, com valores ficando acima de 0,5000 (50% de concordância). Observando os valores do teste t de Student, nota-se que os resultados foram estatisticamente significantes, com valores de p ficando abaixo de 0,05 (ver item t de Student nas Tabelas 1 e 2).

Com isso, conclui-se que as simplificações do método ER podem ser utilizadas com segurança no processo de estimação, trazendo resultados significantes quando comparados com a média dos laudos médicos, e com a vantagem de operar de forma automatizada. Isso contribui para a redução da subjetividade presente nos métodos clássicos, oferecendo também uma solução apropriada para estimação da idade óssea em grandes bases de dados.

 

REFERÊNCIAS

1. Haiter Neto F, Almeida SM, Leite CC. Estudo comparativo dos métodos de estimativa da idade óssea de Greulich & Pyle e Tanner & Whitehouse. Pesqui Odontol Bras. 2000;14:378-84.         [  ]

2. Oliveira HW, Veeck EB, Souza PHC, et al. Avaliação radiográfica da idade óssea em crianças infectadas pelo HIV por via vertical. Radiol Bras. 2006;39:27-31.         [  ]

3. Tanner JM, Whitehouse RH, Healy MJR. A new system for estimating skeletal maturity from the hand and wrist, with standards derived from a study of 2,600 healthy British children. Part II: The Scoring System. Paris: International Children's Centre; 1962.         [  ]

4. Tavano O. Estudo das principais tabelas de avaliação da idade biológica, através do desenvolvimento ósseo, visando sua aplicação em brasileiros leucodermas da região de Bauru [tese de doutorado]. Bauru: Faculdade de Odontologia de Bauru - Universidade de São Paulo; 1976.         [  ]

5. Taffarel B, Wangenhein AV, Silva PCA. Processamento de imagens para análise e estimativa de problemas de crescimento ósseo em crianças. III Workshop de Informática Aplicada à Saúde - CBComp 2003; Florianópolis, SC, Brasil; 2003.         [  ]

6. Haiter-Neto F, Kurita LM, Menezes AV, et al. Skeletal age assessment: a comparison of 3 methods. Am J Orthod Dentofacial Orthop. 2006;130:435.e15-20.         [  ]

7. Olivete CJ. Estimativa da idade óssea através da análise carpal baseada na simplificação do método de Eklof & Ringertz [dissertação de mestrado]. São Carlos: Departamento de Engenharia Elétrica - Universidade de São Paulo; 2005.         [  ]

8. Ayres M. BioEstat. Versão 3.0. Belém: Sociedade Civil Mamirauá, MCT - CNPq; 2005. [acessado em 10 de fevereiro de 2009]. Disponível em: www.mamiraua.org.br        [  ]

9. Nascimento MZ, Frere AF, Marques MA, et al. Subtraction of images for automatic determination of center of the radiation field. Engineering in Medicine and Biology Society. Proceedings of the 25th Annual International Conference of the IEEE. 2003;1:937-40.         [  ]

10. Olivete CJ, Rodrigues ELL, Nascimento MZ. O efeito da correção do "efeito heel" em imagens radiográficas da mão. Rev Bras Fís Méd. 2005; 1:38-51.         [  ]

11. Gonzalez RC, Woods RE. Digital image processing techniques. New York: Addison Wesley; 1993.         [  ]

12. Otsu N. A threshold selection method from gray-level histograms. IEEE Trans Systems Man Cybernet. 1997;9:62-6.         [  ]

13. Abutaleb AS. Automatic thresholding of gray-level pictures using two dimensional entropy. Comput Vision Graph Image Process. 1989;47:22-32.         [  ]

14. Kapur JN, Sahoo PK, Wong AKC. A new method for gray-level picture thresholding using the entropy of the histogram. Comput Vision Graph Image Process. 1985;29:273-85.         [  ]

15. Olivete CJ, Rodrigues ELL. Software automático para determinação da idade óssea baseado no método de Eklof & Ringertz. Rev Bras Fís Méd. 2009;2:15-9. [acessado em 25 de março de 2009]. Disponível em: http://www.abfm.org.br/rbfm/publicado/au_tra_72_p15-19.pdf        [  ]

16. Duarte HEM, Viek R, Siqueira DF, et al. Avaliação das idades dentária e óssea em indivíduos situados antes do surto de crescimento puberal. Rev Ortodontia SPO. 2008;41:95-100.         [  ]

 

 

Endereço para correspondência:
Dr. Celso Olivete Júnior
Departamento de Engenharia Elétrica - EESC-USP
Avenida Trabalhador São-carlense, 400, Centro
São Carlos, SP, Brasil, 13566-590
E-mail: olivete@usp.br

Recebido para publicação em 8/4/2009.
Aceito, após revisão, em 3/11/2009.

 

 

* Trabalho realizado na Escola de Engenharia de São Carlos - Universidade de São Paulo (EESC-USP), São Carlos, SP, Brasil.