メッシュ上のVTK/ITKダイス類似度係数

Question

私はVTKを初めて使用しており、2メッシュから始まるダイ類似度係数（DSC）を計算しようとしています。 DSCはVabメッシュAの間で重複ボリュームであり、メッシュ（3Dスライサー、https://www.slicer.orgを用い.vtk形式でエクスポートすなわち器官輪郭）を読み出すためにB.

メッシュIは、次のスニペットを使用2 Vab/(Va + Vb)、として計算することができます。 （私はVTKと3Dスライサーを用いて計算するものと計算されたものとの間にいくつかの違いを観察したが）

string inputFilename1 = "organ1.vtk"; 

// Get all data from the file 
vtkSmartPointer<vtkGenericDataObjectReader> reader1 = vtkSmartPointer<vtkGenericDataObjectReader>::New(); 
reader1->SetFileName(inputFilename1.c_str()); 
reader1->Update(); 

vtkSmartPointer<vtkPolyData> struct1 = reader1->GetPolyDataOutput(); 

IはvtkMassPropertiesを使用して、2つのメッシュの体積を計算することができます。

次に、2メッシュを交差させるために、vtkIntersectionPolyDataFilterを使用しようとしています。ただし、このフィルタの出力は、入力されたvtkPolyDataオブジェクトの交差点をマークし、閉じたサーフェスではありません。したがって、私はこれらの行から何らかの形でメッシュを生成し、そのボリュームを計算する必要があります。

このようなメッシュを生成し、それを行う方法はどれが良い、正確な方法でしょうか？

---

あるいは、私もITKを使用しようとしました。この問題（2010年12月http://www.insight-journal.org/browse/publication/762）を処理するパッケージが見つかりましたが、ITKの最新バージョンに対してコンパイルできません。 ITKは（現在廃止されている）ITK_USE_REVIEWフラグONでコンパイルする必要があります。言うまでもなく、私はModule_ITKReviewをONに設定してコンパイルしていましたが、下位互換性はありましたが運がなかったのです。

---

最後に、この問題を解決するための他の代替（スクリプト可能）ソフトウェア/ライブラリがある場合は、教えてください。私はこれらの計算を自動的に実行する必要があります。

Answer 1

あなたはあなたのデータは滑らかで、よく振る舞っている場合、このフィルタはかなり良い作品vtkBooleanOperationPolyDataFilter

http://www.vtk.org/doc/nightly/html/classvtkBooleanOperationPolyDataFilter.html

filter->SetOperationToIntersection(); 

を試みることができます。しかし、鋭い構造、例えば。バイナリイメージマーチングキューブアルゴリズムに由来するものは、それを問題にする可能性があります。つまり、vtkPolyDataToImageStencilは必ずしもこの点に関してはそれほど優れているとは限りません。

ポリゴンのブール演算が100kポリゴン以上のサイズの「オルガン」には本当に理想的ではないという印象を一度も持ちました。依存します。

Answer 2

ダイ類似度係数を計算する場合は、まずvtkPolyDataToImageStencilを使用してメッシュからボリューム（ラスタライズ）を生成することをお勧めします。 DSCを計算するのは簡単です。

幸運:)