分割三维膝盖骨
为了量化和测量体积中组分的属性,分割是必要的第一步. 为在 MRT 体积中分割骨组织,需要使用聚类算法来实现粗糙的分割,并通过运用生长切割(grow-cut)算法获取最终结果.
In[1]:= | ![]() X |
Out[1]= | ![]() |
用 MedianFilter 对体积进行预加工来调整噪声. 通过使用 ClusteringComponents 进行聚类立体像素强度来完成初始分割. 这将数据划分为三个区域:空隙、肌肉组织和骨头、脂肪和皮肤等.
In[2]:= | ![]() X |
Out[2]= | ![]() |
用最高平均强度提取分割,它描绘了骨骼、脂肪和皮肤组织.
In[3]:= | ![]() X |
Out[3]= | ![]() |
In[4]:= | ![]() X |
Out[4]= | ![]() |
GrowCutComponents 可提供精密的最终分割. 为了创建骨骼和非骨骼部分的标记,可使用形态学操作.
使用半径4腐蚀体节可计算骨骼标记, 它删除所有薄的皮肤和脂肪层并腐蚀部分骨骼.
In[5]:= | ![]() X |
Out[5]= | ![]() |
对于非骨骼标记,膨胀骨骼核心,大于被腐蚀部分,并提取被膨胀体积周围的表面.
In[6]:= | ![]() X |
Out[6]= | ![]() |
运行三维生长切割(grow-cut)算法细化分割.
In[7]:= | ![]() X |
Out[7]= | ![]() |
可视化分割部分.
In[8]:= | ![]() X |
Out[8]= | ![]() |
可对每个分割进行尺寸和密度测量计算.
In[9]:= | ![]() X |
Out[9]= | ![]() |
Out[37]//TraditionalForm= | |
![]() |