3Dでの膝の骨の分割
立体のコンポーネントの特性を数量化して測定するためには,分割は欠かせない第一ステップである.MRT立体で骨組織を分割するためには,クラスタリングアルゴリズムを使っておおまかな分割を実行して,grow-cutアルゴリズムを使って最終的な結果を得る.
Out[1]= |  |
ノイズを正規化するために,立体はMedianFilterを使って前処理される.最初の分割はClusteringComponents使ってボクセル値をクラスタリングすることで実行される.これによりデータはなにもない空間,筋肉,その他(骨,脂肪,皮膚等)の3つの領域に分けられる.
Out[2]= |  |
骨,脂肪,皮膚の組織を描写する最大の平均強度を持つ部分を抽出する.
Out[3]= |  |
Out[4]= |  |
GrowCutComponentsはきれいな最終的な分割を提供する.骨と骨ではない部分に2つのマーカーを生成するためには,モルフォロジ操作を行う.
骨のマーカーは半径4を使った部分を除去することにより計算することができる.これで薄い皮膚と脂肪層すべてがなくなり,骨の一部も除去される.
Out[5]= |  |
骨ではない部分のマーカーは,除去されたもの以上に骨の中心を拡張し,拡張された立体の周囲の面を抽出する.
Out[6]= |  |
分割を調整するために3Dのgrow-cutアルゴリズムを実行する.
Out[7]= |  |
分割を可視化する.
Out[8]= |  |
各部分について大きさと密度が計算できる.
Out[9]= |  |
Out[37]//TraditionalForm= |
|  |