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.

In[1]:=
Click for copyable input
X
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.

In[2]:=
Click for copyable input
X
Out[2]=

Extraia o segmento com maior intensidade média, que retrata ossos, gordura e pele.

In[3]:=
Click for copyable input
X
Out[3]=
In[4]:=
Click for copyable input
X
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.

In[5]:=
Click for copyable input
X
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.

In[6]:=
Click for copyable input
X
Out[6]=

Execute o algoritmo grow-cut 3D para refinar o segmento.

In[7]:=
Click for copyable input
X
Out[7]=

Visualize a segmentação.

In[8]:=
Click for copyable input
X
Out[8]=

Medidas de tamanho e densidade podem ser calculadas para cada segmento.

In[9]:=
Click for copyable input
X
Out[9]=
Out[37]//TraditionalForm=
de en es ja zh