[vtkusers] Export mesh of cutting plane

Wed Oct 18 10:00:41 EDT 2017

Hello, sorry for coming back to this topic, but I have found some cases
where this doesn't work.

I attached a screenshot.
On top is the 3d mesh cut by a plane parallel to the screen. We can see
that are some elements (hexahedrons) that got cut and they were split in
sets of tetrahedrons.
At the bottom, we see the result of calling vtkCutter ->
vtkCleanUnstructuredGrid -> vtkXMLUnstructuredGridWriter. The result are
all quads and all the generated nodes by the cut are gone, and we see no
triangles.
I would expect to have the same result that we see on the top, but with 2d
elements instead of 3d.

Here's my code:

void cutModel(vtkSmartPointer<vtkUnstructuredGrid> model, double normal[3])
{
vtkSmartPointer<vtkPlane> plane = vtkSmartPointer<vtkPlane>::New();
plane->SetOrigin(0.0, 0.0, 0.0);
plane->SetNormal(normal);

vtkSmartPointer<vtkCutter> cutter = vtkSmartPointer<vtkCutter>::New();
cutter->SetCutFunction(plane);
cutter->SetInputData(model);
cutter->GenerateTrianglesOff();
cutter->Update();

vtkSmartPointer<vtkCleanUnstructuredGrid> cleaner =
vtkSmartPointer<vtkCleanUnstructuredGrid>::New();
cleaner->SetInputConnection(cutter->GetOutputPort());
cleaner->Update();
vtkSmartPointer<vtkXMLUnstructuredGridWriter> writer =
vtkSmartPointer<vtkXMLUnstructuredGridWriter>::New();
writer->SetFileName("planecut.vtu");
writer->SetInputConnection(cleaner->GetOutputPort());
writer->SetDataModeToAscii();
writer->Write();
}

Thank you.

On 29 September 2017 at 13:23, Fernando Nellmeldin <
f.nellmeldin at open-engineering.com> wrote:

> Hello, that worked as a charm, thank you very much for your time!!
>
> On 29 September 2017 at 11:31, kenichiro yoshimi <rccm.kyoshimi at gmail.com>
> wrote:
>
>> Hi Fernando,
>>
>> This will be accomplished by using vtkCleanUnstructuredGrid class
>> found in ParaView without triangulation. I have attached a simple
>> code.
>>
>> Thanks
>>
>> 2017-09-29 16:37 GMT+09:00 Fernando Nellmeldin
>> <f.nellmeldin at open-engineering.com>:
>> > Hello. That did the job! (although I get a triangle mesh while the
>> initial
>> > surface was full of sliced hexahedron (quads).
>> >
>> > Thank you!
>> >
>> > On 29 September 2017 at 04:05, kenichiro yoshimi <
>> rccm.kyoshimi at gmail.com>
>> > wrote:
>> >>
>> >> Hello,
>> >>
>> >> Using vtkDataSetTriangleFilter is a simple way which allows you to
>> >> convert polyData to unstructuredGrid. Something like this (VTK-8.0.0):
>> >> ---
>> >>   // Create cutter
>> >>   vtkSmartPointer<vtkCutter> cutter =
>> >>     vtkSmartPointer<vtkCutter>::New();
>> >>   cutter->SetCutFunction(plane);
>> >>   cutter->SetInputConnection(model->GetOutputPort());
>> >>
>> >>   // Convert polyData to unstructuredGrid
>> >>   vtkSmartPointer<vtkDataSetTriangleFilter> triFilter =
>> >>     vtkSmartPointer<vtkDataSetTriangleFilter>::New();
>> >>   triFilter->SetInputConnection(cutter->GetOutputPort());
>> >>   triFilter->Update();
>> >>
>> >>   vtkSmartPointer<vtkXMLUnstructuredGridWriter> writer =
>> >>     vtkSmartPointer<vtkXMLUnstructuredGridWriter>::New();
>> >>   writer->SetFileName("unstructuredGrid.vtu");
>> >>   writer->SetInputConnection(triFilter->GetOutputPort());
>> >>   writer->Write();
>> >> ---
>> >>
>> >> Regards
>> >>
>> >> 2017-09-28 23:43 GMT+09:00 Fernando Nellmeldin
>> >> <f.nellmeldin at open-engineering.com>:
>> >> > Hello.
>> >> > I have a method to cut a volume in a vtkUnstructuredGrid given a
>> >> > vtkPlane.
>> >> > I'm using vtkTableBasedClipDataSet to extract half of the model,
>> also as
>> >> > vtkUnstructuredGrid. This is working.
>> >> >
>> >> > Now, I would like to obtain ONLY the surface of the cut, than means:
>> I
>> >> > want
>> >> > a vtkUnstructuredGrid with the nodes and faces being on the surface
>> of
>> >> > the
>> >> > cut.
>> >> > How can this be achieved?
>> >> >
>> >> > I tried using vtkCutter with vtkStripper but the mesh I get doesn't
>> have
>> >> > any
>> >> > cells, only the points and their associated data.
>> >> >
>> >> > Here's how my code looks like.
>> >> >
>> >> > Thank You!
>> >> >
>> >> > vtkSmartPointer<vtkUnstructuredGrid> model =
>> >> > vtkSmartPointer<vtkUnstructuredGrid>::New();
>> >> > // fill the ugrid ...
>> >> >
>> >> > // Create The Plane to cut
>> >> > vtkSmartPointer<vtkPlane> plane = vtkSmartPointer<vtkPlane>::New();
>> >> > plane->SetOrigin(0.0, 0.0, 0.0);
>> >> > plane->SetNormal(1.0, 0, 0);
>> >> >
>> >> > // Code to cut the model in half
>> >> > vtkSmartPointer<vtkTableBasedClipDataSet> clipDataSet =
>> >> > vtkSmartPointer<vtkTableBasedClipDataSet>::New();
>> >> > clipDataSet->SetClipFunction(plane);
>> >> > clipDataSet->InsideOutOn();
>> >> > clipDataSet->GenerateClippedOutputOn();
>> >> > clipDataSet->SetInputConnection(model->GetProducerPort());
>> >> > clipDataSet->Update();
>> >> >
>> >> > // This has half the model, this is Ok but not what I want
>> >> > vtkSmartPointer<vtkUnstructuredGrid> halfModel =
>> >> > clipDataSet->GetOutput();
>> >> >
>> >> > ////****////
>> >> > // BEGINS code not working to extract the surface of the cut
>> >> > vtkSmartPointer<vtkCutter> cutter = vtkSmartPointer<vtkCutter>::Ne
>> w();
>> >> > cutter->SetCutFunction(plane);
>> >> > cutter->SetInputConnection(model->GetProducerPort());
>> >> > cutter->Update();
>> >> >
>> >> > vtkSmartPointer<vtkStripper> stripper =
>> >> > vtkSmartPointer<vtkStripper>::New();
>> >> > stripper->SetInputConnection(cutter->GetOutputPort());
>> >> > stripper->Update();
>> >> >
>> >> > vtkSmartPointer<vtkPolyData> pd = vtkSmartPointer<vtkPolyData>::
>> New();
>> >> > pd->SetPoints(stripper->GetOutput()->GetPoints());
>> >> > pd->SetPolys(stripper->GetOutput()->GetLines());
>> >> >
>> >> > vtkSmartPointer<vtkUnstructuredGrid> clipped =
>> >> > vtkSmartPointer<vtkUnstructuredGrid>::New();
>> >> > clipped->ShallowCopy(pd);
>> >> >
>> >> > // Here I see that there is no Cells
>> >> > clipped->Print(std::cout);
>> >> >
>> >> > // ENDS code to extract the surface of the cut
>> >> > ////****////
>> >> >
>> >> > // Write the unstructured grid
>> >> > vtkSmartPointer<vtkXMLUnstructuredGridWriter> writer =
>> >> > vtkSmartPointer<vtkXMLUnstructuredGridWriter>::New();
>> >> > writer->SetFileName("d:/unstructuredGrid.vtu");
>> >> > writer->SetInput(clipped);
>> >> > writer->SetDataModeToAscii();
>> >> > writer->Write();
>> >> > // file doesnt have any cell, only points and pointdata
>> >> >
>> >> > _______________________________________________
>> >> > Powered by www.kitware.com
>> >> >
>> >> > Visit other Kitware open-source projects at
>> >> > http://www.kitware.com/opensource/opensource.html
>> >> >
>> >> > Please keep messages on-topic and check the VTK FAQ at:
>> >> > http://www.vtk.org/Wiki/VTK_FAQ
>> >> >
>> >> > Search the list archives at: http://markmail.org/search/?q=vtkusers
>> >> >
>> >> > Follow this link to subscribe/unsubscribe:
>> >> > http://public.kitware.com/mailman/listinfo/vtkusers
>> >> >
>> >
>> >
>> >
>> >
>> > --
>> > Fernando NELLMELDIN
>> > Software Engineer
>> > _______________________________________________________________
>> >
>> > Open Engineering s.a.
>> >
>> > Rue Bois Saint-Jean 15/1
>> > B-4102 Seraing (Belgium)
>> > Tel: +32.4.353.30.34
>> >
>> > http://www.open-engineering.com
>> > https://www.linkedin.com/company/open-engineering?trk=biz-companies-cym
>> > ____________________________________________________________
>> _____________
>>
>
>
>
> --
> *Fernando NELLMELDIN*
> Software Engineer
> *_______________________________________________________________*
>
> *Open Engineering s.a.*
>
> Rue Bois Saint-Jean 15/1
> B-4102 Seraing (Belgium)
> Tel: +32.4.353.30.34
>
> http://www.open-engineering.com
> https://www.linkedin.com/company/open-engineering?trk=biz-companies-cym
>
> *_________________________________________________________________________*
>

-- 
*Fernando NELLMELDIN*
Software Engineer
*_______________________________________________________________*

*Open Engineering s.a.*

Rue Bois Saint-Jean 15/1
B-4102 Seraing (Belgium)
Tel: +32.4.353.30.34

http://www.open-engineering.com
https://www.linkedin.com/company/open-engineering?trk=biz-companies-cym

*_________________________________________________________________________*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/vtkusers/attachments/20171018/c9c42847/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: planecutcomparison.png
Type: image/png
Size: 74265 bytes
Desc: not available
URL: <http://public.kitware.com/pipermail/vtkusers/attachments/20171018/c9c42847/attachment-0001.png>