Hi,<br><br>I&#39;m trying to generate a mesh for export and analysis via finite element. <br><br>But the original mesh are very triangle. <br><br><div style="text-align: left;" id="result_box" dir="ltr">If you apply the filter mesh decimação is a little deformed.<br>
<br><div style="text-align: left;" id="result_box" dir="ltr">The Imagedata is generated from a CT 512x512 with 200 slices.</div><br>
<div style="text-align: left;" id="result_box" dir="ltr">What can I doing wrong ?</div>
<br><br>The Code:<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mcubes = vtk.vtkMarchingCubes()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mcubes.GetOutput().ReleaseDataFlagOff()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mcubes.SetInput(_imagedata)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mcubes.SetValue(0, 4)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mcubes.ComputeGradientsOn()<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mcubes.Update()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smoother = vtk.vtkSmoothPolyDataFilter()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smoother.GetOutput().ReleaseDataFlagOff()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smoother.SetInput(mcubes.GetOutput())<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smoother.SetNumberOfIterations(5)<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smoother.SetRelaxationFactor(0.70)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smoother.SetFeatureAngle(70)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smoother.FeatureEdgeSmoothingOn()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smoother.BoundarySmoothingOn()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smoother.Update()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; normals = vtk.vtkPolyDataNormals()<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; normals.GetOutput().ReleaseDataFlagOff()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; normals.SetInput(smoother.GetOutput())<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; normals.SetFeatureAngle(80)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; normals.Update()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; stripper = vtk.vtkStripper()<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; stripper.SetInput(normals.GetOutput())<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; stripper.GetOutput().ReleaseDataFlagOff()<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; stripper.Update()<br>&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mesh =&nbsp; normals.GetOutput()<br><br><br><br>Regards<br>Paul Jhin<br></div>