Segmentação de um osso do joelho em 3D
Para quantificar e medir as propriedades de um componente em um volume, a segmentação é um primeiro passo necessário. Para segmentar o tecido ósseo em um volume MRT, um algoritmo de clusterização é usado para obter uma segmentação grosseira e aplique um algoritmo grow-cut para obter o resultado final.
Out[1]= | |
O volume é pré-processado usando MedianFilter para regularizar o ruído. Uma segmentação inicial é realizada através da clusterização das intensidades dos voxels usando ClusteringComponents. Os dados são assim divididos em três regiões: espaço vazio, tecido muscular e o resto, que inclui ossos, gordura e pele.
Out[2]= | |
Extraia o segmento com maior intensidade média, que retrata ossos, gordura e pele.
Out[3]= | |
Out[4]= | |
GrowCutComponents pode fornecer uma boa segmentação final. A fim de criar dois marcadores para as áreas ósseas e não ósseas, pode-se usar operações morfológicas.
O marcador de osso pode ser computado pela erosão do segmento usando um raio de tamanho 4, que exclui todas as camadas de pele fina e gordura, e também parte do osso.
Out[5]= | |
Para o marcador de áreas não ósseas, dilate o núcleo do osso mais do que foi erodido e extraia a superfície ao redor da área dilatada.
Out[6]= | |
Execute o algoritmo grow-cut 3D para refinar o segmento.
Out[7]= | |
Visualize a segmentação.
Out[8]= | |
Medidas de tamanho e densidade podem ser calculadas para cada segmento.
Out[9]= | |
Out[37]//TraditionalForm= |
| |