Wie man das Volumen einer vtk-Datei berechnet

Question

Ich bin nur neu bei Python und ich kann nicht scheinen, eine Lösung für mein Problem zu finden, da es scheint, ziemlich einfach zu sein. Ich habe eine Geometrie auf Paraview, ich speichere es als vtk-Datei und ich versuche Python zu verwenden, um sein Volumen zu berechnen.

Dies ist der Code Ich verwende:

import vtk 
reader = vtk.vtkPolyDataReader() 
reader.SetFileName("C:\Users\Pauuu\Google Drive\2016-01\SURF\Sim Vascular\Modelos\apoE183 Day 14 3D\AAA.vtk") 
reader.Update() 
polydata = reader.GetOutput() 

Mass = vtk.vtkMassProperties() 
Mass.SetInputConnection(polydata.GetOutput()) 
Mass.Update() 

print "Volume = ", Mass.GetVolume() 
print "Surface = ", Mass.GetSurfaceArea() 

Ich denke, es könnte die Daten ein Problem mit der Art und Weise im Loding sein, und ich bekomme die AttributeError: GetOutput.

Wissen Sie, was passiert oder was ich falsch mache? Vielen Dank im Voraus.

Answer 1

Sie haben reader.GetOutput() in polydata zugewiesen. Von polydata, ich glaube, Sie tun müssen, polydata.GetOutputPort()

Answer 2

Je nach Ihrer Version von vtk Paket erhalten Sie die folgende Syntax, wenn Ihre Version < = 5 testen möchten:

Mass.SetInput(polydata.GetOutput()); 

Andernfalls wird die tatsächliche Syntax :

Mass.SetInputData(polydata.GetOutputPort()); 

PS: Sie können die python-wrapped überprüfen vtk Version durch Ausführen von:

import vtk 
print vtk.vtkVersion.GetVTKSourceVersion() 

Answer 3

Ich denke, Sie haben VTK 6, Sie als Eingabe für einen Filter entweder der Ausgangstor eines Filters oder ein vtkDataObject zur Verfügung stellen können:

Mass.SetInputConnection(reader.GetOutputPort()) 
Mass.SetInputData(polydata) #that is Mass.SetInputData(reader.GetOutput()) 

Für das Verständnis, warum diese Methode nicht gleichwertig ist, wenn eine Aktualisierung Pipeline und zum Vergleich mit der vorherigen Version, siehe http://www.vtk.org/Wiki/VTK/VTK_6_Migration/Removal_of_GetProducerPorthttp://www.vtk.org/Wiki/VTK/VTK_6_Migration/Replacement_of_SetInput