FIX: When the user split, select the largest, or various surface, the vtkMassPro…

…perties because the resulting was not a homogenic triangle surface, so it is necessary to apply the vtkTriangleFilter

FIX: When the user split, select the largest, or various surface, the vtkMassPro…
…perties because the resulting was not a homogenic triangle surface, so it is necessary to apply the vtkTriangleFilter
tfmoraes
1 parent 8cbb2a51
Showing 1 changed file with 6 additions and 1 deletions Show diff stats
invesalius/data/surface.py
@@ -271,10 +271,15 @@ class SurfaceManager():
  
         # The following lines have to be here, otherwise all volumes disappear
         if not volume:
+            triangle_filter = vtk.vtkTriangleFilter()
+            triangle_filter.SetInput(polydata)
+            triangle_filter.Update()
+
             measured_polydata = vtk.vtkMassProperties()
-            measured_polydata.SetInput(polydata)
+            measured_polydata.SetInput(triangle_filter.GetOutput())
             volume =  measured_polydata.GetVolume()
             surface.volume = volume
+            print ">>>>", surface.volume
         else:
             surface.volume = volume
         self.last_surface_index = surface.index
...	...	@@ -271,10 +271,15 @@ class SurfaceManager():
271	271
272	272	# The following lines have to be here, otherwise all volumes disappear
273	273	if not volume:
	274	+ triangle_filter = vtk.vtkTriangleFilter()
	275	+ triangle_filter.SetInput(polydata)
	276	+ triangle_filter.Update()
	277	+
274	278	measured_polydata = vtk.vtkMassProperties()
275		- measured_polydata.SetInput(polydata)
	279	+ measured_polydata.SetInput(triangle_filter.GetOutput())
276	280	volume = measured_polydata.GetVolume()
277	281	surface.volume = volume
	282	+ print ">>>>", surface.volume
278	283	else:
279	284	surface.volume = volume
280	285	self.last_surface_index = surface.index
...	...