分割三维膝盖骨
为了量化和测量体积中组分的属性,分割是必要的第一步. 为在 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= | |
![]() | |









