计算 volume of 3D .vtk points (triangle mesh, closed polydata)

问题描述：有一个.vtk点集，带triangle mesh，这些点构成一个闭合的三维物体表面。想计算这些点围成的这个closed object的内部体积。

解决思路：用vtkMassProperties就像这里对vtkMassProperties的提示: Currently only triangles are processed. Use vtkTriangleFilter to convert any strips or polygons to triangles. 也就是说，目前该类只能用于triangle mesh connectivity closed surface.

代码，polyDataVolume.cpp文件：

#include <vtkPolyDataReader.h>
#include "vtkPolyData.h"
#include "vtkSmartPointer.h"
#include <vtkMassProperties.h>


int main( int argc, char* argv[] ){

    // read in the .vtk file
    vtkSmartPointer<vtkPolyDataReader> reader = vtkSmartPointer<vtkPolyDataReader>::New();
    reader->SetFileName(argv[1]);
    reader->Update();

    vtkSmartPointer<vtkPolyData> polyData = vtkSmartPointer<vtkPolyData>::New();
    polyData = reader->GetOutput();

    // evaluate the volume
    vtkMassProperties *mass = vtkMassProperties::New();

    mass->SetInput(polyData);

    cout << "---- the volume of this vtk polydata is: " << mass->GetVolume() << endl;

    return 0;
}

CMakeLists.txt文件：

cmake_minimum_required(VERSION 2.8)
 
PROJECT(polyDataVolume)
 
find_package(VTK REQUIRED)
include(${VTK_USE_FILE})
 
add_executable(polyDataVolume MACOSX_BUNDLE polyDataVolume)
 
if(VTK_LIBRARIES)
  target_link_libraries(polyDataVolume ${VTK_LIBRARIES})
else()
  target_link_libraries(polyDataVolume vtkHybrid vtkWidgets)
endif()

编译：

ccmake ../

make

运行：

./polyDataVolume ../test.vtk

运行结果：

---- the volume of this vtk polydatais: 0.00339752

这个小实验的代码已打包，可点这里下载。

参考：

http://www.vtk.org/pipermail/vtkusers/2013-June/080004.html

http://www.vtk.org/Wiki/VTK/Examples/Cxx/PolyData/EmbedPointsIntoVolume

http://markmail.org/message/6u6psqha5ckz2ec3

http://www.paraview.org/Wiki/VTK/Examples/Cxx/PolyData/PolyDataToImageData

http://permalink.gmane.org/gmane.comp.lib.vtk.user/67470

http://public.kitware.com/pipermail/vtkusers/2011-January/065437.html

http://vtk.1045678.n5.nabble.com/Polydata-volume-calculation-problem-td3347493.html

http://public.kitware.com/pipermail/vtkusers/2014-July/084469.html

http://www.paraview.org/Wiki/VTK/Examples/Cxx/Modelling/Delaunay3D