[vtkusers] How to visualize polydata after KmeansClustering?

Thu Aug 30 16:16:35 EDT 2018

Send me the working codr and I'll ad it as a VTKExample. You'll get credit
of course.

On Thu, Aug 30, 2018, 9:06 AM Polly Pui <polly_sukting at hotmail.com> wrote:

>
> Hi,
> The code works well now. I could visualize it.
> I made a mistake by changing into a different line.
> Thank you so much.
>
> Polly
> ------------------------------
> *From:* vtkusers <vtkusers-bounces at public.kitware.com> on behalf of Polly
> Pui <polly_sukting at hotmail.com>
> *Sent:* Thursday, August 30, 2018 12:22 PM
> *To:* kenichiro yoshimi
> *Cc:* vtkusers
> *Subject:* Re: [vtkusers] How to visualize polydata after
> KmeansClustering?
>
>
> Hi,
>
> It doesn't work.
>
> I changed the lines, but the identifier "point" is undefined.
> ------------------------------------------------
> And I have another question here,
> These codes could read polydata vtk, but failed to visualize.
> I have another dataset (approximate 2000 set) that I saved by using
> Paraview into .vtk (the data is saved into unstructured grid data instead
> of polydata)
> Then, I tried to open using the following codes but it failed to load for
> unstructured grid data.
> How should i modify the lines to read the unstructured grid dataset?
>
> Thank you so much.
>
> Polly
> ------------------------------
> *From:* kenichiro yoshimi <rccm.kyoshimi at gmail.com>
> *Sent:* Thursday, August 30, 2018 9:16 AM
> *To:* polly_sukting at hotmail.com
> *Cc:* vtkusers
> *Subject:* Re: [vtkusers] How to visualize polydata after
> KmeansClustering?
>
> Hi,
>
> Could you try to modify the lines of your code as below?
>
>   vtkSmartPointer<vtkPolyData> polydata =
>     vtkSmartPointer<vtkPolyData>::New();
> ==>
>   vtkSmartPointer<vtkPolyData> polydata = point;
>
> Regards
> 2018年8月26日(日) 1:59 Polly Pui <polly_sukting at hotmail.com>:
> >
> > Hi,
> >
> > I would like to show my polydata with points and mark with colours after
> the KmeansClustering process.
> >
> > However, I could only get the values of the points without the
> visualization of the data.
> >
> > Can anyone please advice?
> >
> > Thank you so much.
> >
> >
> > Polly
> >
> >
> --------------------------------------------------------------------------------------
> >
> > My code is as below:
> >
> > int main(int argc, char* argv[])
> > {
> > // Get the points into the format needed for KMeans
> > vtkSmartPointer<vtkPolyData> point =
> > vtkSmartPointer<vtkPolyData>::New();
> >
> > vtkSmartPointer<vtkPolyDataReader> reader =
> > vtkSmartPointer<vtkPolyDataReader>::New();
> > reader->SetFileName(argv[1]);
> > reader->Update();
> >
> > point = reader->GetOutput();
> >
> > vtkSmartPointer<vtkTable> inputData =
> > vtkSmartPointer<vtkTable>::New();
> >
> > for (int c = 0; c < 3; ++c)
> > {
> > std::stringstream colName;
> > colName << "coord " << c;
> > vtkSmartPointer<vtkDoubleArray> doubleArray =
> > vtkSmartPointer<vtkDoubleArray>::New();
> > doubleArray->SetNumberOfComponents(1);
> > doubleArray->SetName(colName.str().c_str());
> > doubleArray->SetNumberOfTuples(point->GetNumberOfPoints());
> >
> > for (int r = 0; r < point->GetNumberOfPoints(); ++r)
> > {
> > double p[3];
> > point->GetPoint(r, p);
> > doubleArray->SetValue(r, p[c]);
> > }
> > inputData->AddColumn(doubleArray);
> > }
> >
> > vtkSmartPointer<vtkKMeansStatistics> kMeansStatistics =
> > vtkSmartPointer<vtkKMeansStatistics>::New();
> >
> > #if VTK_MAJOR_VERSION <= 5
> > kMeansStatistics->SetInput(vtkStatisticsAlgorithm::INPUT_DATA,
> inputData);
> > #else
> > kMeansStatistics->SetInputData(vtkStatisticsAlgorithm::INPUT_DATA,
> inputData);
> > #endif
> > kMeansStatistics->SetColumnStatus(inputData->GetColumnName(0), 1);
> > kMeansStatistics->SetColumnStatus(inputData->GetColumnName(1), 1);
> > kMeansStatistics->SetColumnStatus(inputData->GetColumnName(2), 1);
> > kMeansStatistics->RequestSelectedColumns();
> > kMeansStatistics->SetAssessOption(true);
> > kMeansStatistics->SetDefaultNumberOfClusters(5);
> > //kMeansStatistics->SetMaxNumIterations(1);
> > kMeansStatistics->Update();
> > //double mean0[3];
> > //kMeansStatistics->GetMean(0, mean0);
> > // Display the results
> > kMeansStatistics->GetOutput()->Dump();
> > //Group the points according to ID number
> > vtkSmartPointer<vtkIntArray> clusterArray =
> > vtkSmartPointer<vtkIntArray>::New();
> > clusterArray->SetNumberOfComponents(1);
> > clusterArray->SetName("ClusterId");
> >
> > for (int r = 0; r < kMeansStatistics->GetOutput()->GetNumberOfRows();
> r++)
> > {
> > vtkVariant v = kMeansStatistics->GetOutput()->GetValue(r,
> kMeansStatistics->GetOutput()->GetNumberOfColumns() - 1);
> > //std::cout << "Point " << r << " is in cluster " << v.ToInt() <<
> std::endl;
> > clusterArray->InsertNextValue(v.ToInt());
> > }
> >
> > // Create a lookup table to map cell data to colors
> > vtkSmartPointer<vtkLookupTable> lut =
> > vtkSmartPointer<vtkLookupTable>::New();
> > int tableSize = (kMeansStatistics + 1, 10);
> > lut->SetNumberOfTableValues(tableSize);
> > lut->Build();
> >
> > // Fill in a few known colors, the rest will be generated if needed
> > //ColorCells tutorial
> https://lorensen.github.io/VTKExamples/site/Cxx/PolyData/ColorCells/
> ColorCells - GitHub Pages
> <https://lorensen.github.io/VTKExamples/site/Cxx/PolyData/ColorCells/>
> lorensen.github.io
> Download and Build ColorCells¶. Click here to download ColorCells and its
> CMakeLists.txt file. Once the tarball ColorCells.tar has been downloaded
> and extracted,
>
>
> > vtkSmartPointer<vtkNamedColors> colors =
> > vtkSmartPointer<vtkNamedColors>::New();
> > lut->SetTableValue(0, colors->GetColor4d("Black").GetData());
> > lut->SetTableValue(1, colors->GetColor4d("Banana").GetData());
> > lut->SetTableValue(2, colors->GetColor4d("Tomato").GetData());
> > lut->SetTableValue(3, colors->GetColor4d("Peacock").GetData());
> > lut->SetTableValue(4, colors->GetColor4d("Lavender").GetData());
> >
> > // Output the cluster centers
> >
> > vtkMultiBlockDataSet* outputMetaDS =
> vtkMultiBlockDataSet::SafeDownCast(kMeansStatistics->GetOutputDataObject(vtkStatisticsAlgorithm::OUTPUT_MODEL));
> > vtkSmartPointer<vtkTable> outputMeta =
> vtkTable::SafeDownCast(outputMetaDS->GetBlock(0));
> > //vtkSmartPointer<vtkTable> outputMeta = vtkTable::SafeDownCast(
> outputMetaDS->GetBlock( 1 ) );
> > vtkDoubleArray* coord0 =
> vtkDoubleArray::SafeDownCast(outputMeta->GetColumnByName("coord 0"));
> > vtkDoubleArray* coord1 =
> vtkDoubleArray::SafeDownCast(outputMeta->GetColumnByName("coord 1"));
> > vtkDoubleArray* coord2 =
> vtkDoubleArray::SafeDownCast(outputMeta->GetColumnByName("coord 2"));
> >
> > for (unsigned int i = 0; i < coord0->GetNumberOfTuples(); ++i)
> > {
> > std::cout << coord0->GetValue(i) << " " << coord1->GetValue(i) << " " <<
> coord2->GetValue(i) << std::endl;
> > }
> >
> > vtkSmartPointer<vtkPolyData> polydata =
> > vtkSmartPointer<vtkPolyData>::New();
> > //polydata->SetPoints(point);
> > polydata->GetPointData()->AddArray(clusterArray);
> > polydata->GetPointData()->SetScalars(clusterArray);
> >
> > // Display
> > vtkSmartPointer<vtkVertexGlyphFilter> glyphFilter =
> > vtkSmartPointer<vtkVertexGlyphFilter>::New();
> > #if VTK_MAJOR_VERSION <= 5
> > glyphFilter->SetInputConnection(polydata->GetProducerPort());
> > #else
> > glyphFilter->SetInputData(polydata);
> > #endif
> > glyphFilter->Update();
> >
> > // Create a mapper and actor
> > vtkSmartPointer<vtkPolyDataMapper> mapper =
> > vtkSmartPointer<vtkPolyDataMapper>::New();
> > mapper->SetInputConnection(glyphFilter->GetOutputPort());
> > mapper->SetScalarRange(0, tableSize - 1);
> > mapper->SetLookupTable(lut);
> >
> > vtkSmartPointer<vtkActor> actor =
> > vtkSmartPointer<vtkActor>::New();
> > actor->SetMapper(mapper);
> > actor->GetProperty()->SetPointSize(5);
> >
> > // Create a renderer, render window, and interactor
> > vtkSmartPointer<vtkRenderer> renderer =
> > vtkSmartPointer<vtkRenderer>::New();
> > vtkSmartPointer<vtkRenderWindow> renderWindow =
> > vtkSmartPointer<vtkRenderWindow>::New();
> > renderWindow->AddRenderer(renderer);
> > renderWindow->SetSize(800, 800);
> >
> > vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
> > vtkSmartPointer<vtkRenderWindowInteractor>::New();
> > renderWindowInteractor->SetRenderWindow(renderWindow);
> >
> > // Add the actor to the scene
> > renderer->AddActor(actor);
> > renderer->SetBackground(colors->GetColor3d("SlateGray").GetData());
> >
> > vtkSmartPointer<vtkInteractorStyleTrackballCamera> style =
> > vtkSmartPointer<vtkInteractorStyleTrackballCamera>::New();
> > renderWindowInteractor->SetInteractorStyle(style);
> >
> > // Render and interact
> > renderWindow->Render();
> >
> > vtkWindowToImageFilter *windowToImageFilter =
> vtkWindowToImageFilter::New();
> > windowToImageFilter->SetInput(renderWindow);
> > windowToImageFilter->Update();
> >
> > vtkJPEGWriter *writer = vtkJPEGWriter::New();
> > writer->SetInputConnection(windowToImageFilter->GetOutputPort());
> > writer->SetFileName("10nosepoints.jpg");
> > writer->Write();
> >
> > renderWindowInteractor->Start();
> >
> > return EXIT_SUCCESS;
> > }
> >
> > _______________________________________________
> > 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:
> > https://public.kitware.com/mailman/listinfo/vtkusers
> _______________________________________________
> 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:
> https://public.kitware.com/mailman/listinfo/vtkusers
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://public.kitware.com/pipermail/vtkusers/attachments/20180830/8e272796/attachment.html>