<div dir="ltr">Are your normals all in the right directions?<br><div><br></div><div>You could try the vtkPolyDataNormals or vktReverseSense filters to try to make them consistent.<br></div><div style><br></div><div style>Seth</div>
<div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Jun 6, 2013 at 9:00 AM, <a href="mailto:vtkusers-request@vtk.org">vtkusers-request@vtk.org</a> <span dir="ltr"><<a href="mailto:vtkusers-request@vtk.org" target="_blank">vtkusers-request@vtk.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Send vtkusers mailing list submissions to<br>
<a href="mailto:vtkusers@vtk.org">vtkusers@vtk.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:vtkusers-request@vtk.org">vtkusers-request@vtk.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:vtkusers-owner@vtk.org">vtkusers-owner@vtk.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of vtkusers digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. computing the volume of a simple multi triangles polydata<br>
(Michele)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 6 Jun 2013 07:33:22 -0700 (PDT)<br>
From: Michele <<a href="mailto:michele.conconi@unibo.it">michele.conconi@unibo.it</a>><br>
Subject: [vtkusers] computing the volume of a simple multi triangles<br>
polydata<br>
To: <a href="mailto:vtkusers@vtk.org">vtkusers@vtk.org</a><br>
Message-ID: <<a href="mailto:1370529201962-5721243.post@n5.nabble.com">1370529201962-5721243.post@n5.nabble.com</a>><br>
Content-Type: text/plain; charset=us-ascii<br>
<br>
Hi everybody,<br>
<br>
I am tring to build a polydata starting from an ordered set of points in<br>
space in order to evaluate its volume.<br>
<br>
I built an exapmple (see below) which generate a correct visual<br>
representation of the polydata, but the computed volume is wrong (6.8 rather<br>
then 48 mm^3).<br>
<br>
I tried aldo the vtkDelaunay2D, with no better results.<br>
<br>
What am I doing wrong?<br>
Thanks in advance for the help<br>
<br>
Michele<br>
<br>
------------------------------<br>
<br>
<br>
#include <vtkVersion.h><br>
#include <vtkSmartPointer.h><br>
#include <vtkCellArray.h><br>
#include <vtkCellData.h><br>
#include <vtkDoubleArray.h><br>
#include <vtkPoints.h><br>
#include <vtkPolyLine.h><br>
#include <vtkPolyData.h><br>
#include <vtkPolyDataMapper.h><br>
#include <vtkActor.h><br>
#include <vtkRenderWindow.h><br>
#include <vtkRenderer.h><br>
#include <vtkRenderWindowInteractor.h><br>
#include <vtkMassProperties.h><br>
#include <vtkTriangle.h><br>
<br>
<br>
<br>
int main( int argc, char *argv[] )<br>
{<br>
// build the points set to be triangulated<br>
int Np = 5;<br>
double **P1, **P2, **P3;<br>
<br>
P1 = new double*[Np];<br>
P2 = new double*[Np];<br>
P3 = new double*[Np];<br>
<br>
for(int i = 0; i < Np; i++)<br>
{<br>
P1[i] = new double[4];<br>
P2[i] = new double[4];<br>
P3[i] = new double[4];<br>
}<br>
<br>
for(int i=0;i<Np;i++)<br>
{<br>
P1[i][0] = 0+i;<br>
P1[i][1] = 1;<br>
P1[i][2] = 0;<br>
<br>
P2[i][0] = 0+i;<br>
P2[i][1] = 3;<br>
P2[i][2] = 2;<br>
<br>
P3[i][0] = 0+i;<br>
P3[i][1] = 5;<br>
P3[i][2] = 4;<br>
}<br>
<br>
<br>
vtkSmartPointer&lt;vtkPoints> points = vtkSmartPointer<vtkPoints>::New();<br>
<br>
for(int i=0;i<Np;i++)<br>
points->InsertNextPoint(P1[i]);<br>
<br>
for(int i=0;i<Np;i++)<br>
points->InsertNextPoint(P2[i]);<br>
<br>
<br>
for(int i=0;i<Np;i++)<br>
points->InsertNextPoint(P3[i]);<br>
<br>
for(int i=0;i<Np;i++)<br>
points->InsertNextPoint(P1[i][0],0,P1[i][2]);<br>
<br>
<br>
for(int i=0;i<Np;i++)<br>
points->InsertNextPoint(P2[i][0],0,P2[i][2]);<br>
<br>
for(int i=0;i<Np;i++)<br>
points->InsertNextPoint(P3[i][0],0,P3[i][2]);<br>
<br>
<br>
//build manually the triangle<br>
<br>
vtkSmartPointer<vtkCellArray> triangles =<br>
vtkSmartPointer<vtkCellArray>::New();<br>
vtkSmartPointer<vtkTriangle> triangle =<br>
vtkSmartPointer<vtkTriangle>::New();<br>
<br>
// ------------------------------------- base<br>
for(int j=3;j<5;j++)<br>
for(int i = (j*Np)+1; i < (j+1)*Np; i++)<br>
{<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,i-1);<br>
triangle->GetPointIds()->SetId(1,i);<br>
triangle->GetPointIds()->SetId(2,Np+i-1);<br>
triangles ->InsertNextCell(triangle);<br>
}<br>
<br>
for(int j=4;j<6;j++)<br>
for(int i = j*Np; i < (j+1)*Np-1; i++)<br>
{<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,i);<br>
triangle->GetPointIds()->SetId(1,i+1);<br>
triangle->GetPointIds()->SetId(2,i-Np+1);<br>
triangles ->InsertNextCell(triangle);<br>
}<br>
<br>
<br>
// ---------------------------------------- side 1<br>
<br>
for(int i = 1; i < Np; i++)<br>
{<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,i-1);<br>
triangle->GetPointIds()->SetId(1,i);<br>
triangle->GetPointIds()->SetId(2,3*Np+i-1);<br>
triangles ->InsertNextCell(triangle);<br>
}<br>
<br>
for(int i = 3*Np; i < 4*Np-1; i++)<br>
{<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,i);<br>
triangle->GetPointIds()->SetId(1,i+1);<br>
triangle->GetPointIds()->SetId(2,i-3*Np+1);<br>
triangles ->InsertNextCell(triangle);<br>
}<br>
<br>
// ---------------------------------------- side 2<br>
for(int i = 2*Np+1; i < 3*Np; i++)<br>
{<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,i-1);<br>
triangle->GetPointIds()->SetId(1,i);<br>
triangle->GetPointIds()->SetId(2,3*Np+i-1);<br>
triangles ->InsertNextCell(triangle);<br>
}<br>
<br>
for(int i = 5*Np; i < 6*Np-1; i++)<br>
{<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,i);<br>
triangle->GetPointIds()->SetId(1,i+1);<br>
triangle->GetPointIds()->SetId(2,i-3*Np+1);<br>
triangles ->InsertNextCell(triangle);<br>
}<br>
<br>
// ---------------------------------------- side 3<br>
<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,0);<br>
triangle->GetPointIds()->SetId(1,3*Np);<br>
triangle->GetPointIds()->SetId(2,4*Np);<br>
triangles ->InsertNextCell(triangle);<br>
<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,Np);<br>
triangle->GetPointIds()->SetId(1,4*Np);<br>
triangle->GetPointIds()->SetId(2,5*Np);<br>
triangles ->InsertNextCell(triangle);<br>
<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,4*Np);<br>
triangle->GetPointIds()->SetId(1,Np);<br>
triangle->GetPointIds()->SetId(2,0);<br>
triangles ->InsertNextCell(triangle);<br>
<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,5*Np);<br>
triangle->GetPointIds()->SetId(1,2*Np);<br>
triangle->GetPointIds()->SetId(2,Np);<br>
triangles ->InsertNextCell(triangle);<br>
<br>
// ---------------------------------------- side 4<br>
<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,Np-1);<br>
triangle->GetPointIds()->SetId(1,4*Np-1);<br>
triangle->GetPointIds()->SetId(2,5*Np-1);<br>
triangles ->InsertNextCell(triangle);<br>
<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,2*Np-1);<br>
triangle->GetPointIds()->SetId(1,5*Np-1);<br>
triangle->GetPointIds()->SetId(2,6*Np-1);<br>
triangles ->InsertNextCell(triangle);<br>
<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,5*Np-1);<br>
triangle->GetPointIds()->SetId(1,2*Np-1);<br>
triangle->GetPointIds()->SetId(2,Np-1);<br>
triangles ->InsertNextCell(triangle);<br>
<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,6*Np-1);<br>
triangle->GetPointIds()->SetId(1,3*Np-1);<br>
triangle->GetPointIds()->SetId(2,2*Np-1);<br>
triangles ->InsertNextCell(triangle);<br>
<br>
<br>
// -------------------- top<br>
<br>
for(int j=0;j<2;j++)<br>
for(int i = (j*Np)+1; i < (j+1)*Np; i++)<br>
{<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,i-1);<br>
triangle->GetPointIds()->SetId(1,i);<br>
triangle->GetPointIds()->SetId(2,Np+i-1);<br>
triangles ->InsertNextCell(triangle);<br>
}<br>
<br>
for(int j=1;j<3;j++)<br>
for(int i = j*Np; i < (j+1)*Np-1; i++)<br>
{<br>
triangle->Initialize();<br>
triangle->GetPointIds()->SetId(0,i);<br>
triangle->GetPointIds()->SetId(1,i+1);<br>
triangle->GetPointIds()->SetId(2,i-Np+1);<br>
triangles ->InsertNextCell(triangle);<br>
}<br>
<br>
<br>
// build the polydata<br>
<br>
vtkPolyData *polyData = vtkPolyData::New();<br>
<br>
polyData->SetPoints(points);<br>
polyData->SetPolys(triangles );<br>
<br>
<br>
//evaluate the volume : expected outcome for the current points set 48<br>
<br>
vtkMassProperties *mass = vtkMassProperties::New();<br>
<br>
mass->SetInput(polyData);<br>
cout<<"volume "<<mass->GetVolume()<<endl;<br>
<br>
// Setup actor and mapper<br>
vtkSmartPointer&lt;vtkPolyDataMapper> mapper =<br>
vtkSmartPointer<vtkPolyDataMapper>::New();<br>
mapper->SetInput(polyData);<br>
<br>
<br>
vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();<br>
actor->SetMapper(mapper);<br>
<br>
// Setup render window, renderer, and interactor<br>
vtkSmartPointer<vtkRenderer> renderer =<br>
vtkSmartPointer<vtkRenderer>::New();<br>
vtkSmartPointer<vtkRenderWindow> renderWindow =<br>
vtkSmartPointer<vtkRenderWindow>::New();<br>
renderWindow->AddRenderer(renderer);<br>
vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =<br>
vtkSmartPointer<vtkRenderWindowInteractor>::New();<br>
renderWindowInteractor->SetRenderWindow(renderWindow);<br>
renderer->AddActor(actor);<br>
<br>
renderWindow->Render();<br>
renderWindowInteractor->Start();<br>
<br>
return 0;<br>
}<br>
<br>
<br>
<br>
--<br>
View this message in context: <a href="http://vtk.1045678.n5.nabble.com/computing-the-volume-of-a-simple-multi-triangles-polydata-tp5721243.html" target="_blank">http://vtk.1045678.n5.nabble.com/computing-the-volume-of-a-simple-multi-triangles-polydata-tp5721243.html</a><br>
Sent from the VTK - Users mailing list archive at Nabble.com.<br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
vtkusers mailing list<br>
<a href="mailto:vtkusers@vtk.org">vtkusers@vtk.org</a><br>
<a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
<br>
<br>
End of vtkusers Digest, Vol 110, Issue 9<br>
****************************************<br>
</blockquote></div><br></div></div></div>