VTK对点进行采样

在VTK中对三维点进行采样，并保存采样后的结果，在同一个窗口显示采样前与采样后的点，采样前的点以红色显示，采样后的点以绿色显示

如果是二维点，则读取二维点坐标就是了

 

#include <iostream>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkPoints.h>
#include <vtkPolyData.h>
#include <vtkSmartPointer.h>
#include <vtkPolyDataMapper.h>
#include <vtkCleanPolyData.h>
#include <vtkPolyDataWriter.h>
#include <vtkActor.h>
#include <vtkProperty.h>
#include <vtkCellArray.h>
#include <vtkInteractorStyleTrackball.h>
void main()
{
	vtkSmartPointer<vtkPoints>m_points=vtkSmartPointer<vtkPoints>::New();
	vtkSmartPointer<vtkCellArray>m_vertices=vtkSmartPointer<vtkCellArray>::New();
	vtkSmartPointer<vtkRenderer>ren1=vtkSmartPointer<vtkRenderer>::New();
	vtkSmartPointer<vtkRenderWindow>renWin=vtkSmartPointer<vtkRenderWindow>::New();
	vtkSmartPointer<vtkRenderWindowInteractor>iren=vtkSmartPointer<vtkRenderWindowInteractor>::New();
	FILE *fp=NULL;
	fp=fopen("右心室浮点型.txt","r");
	if (!fp)
	{
		cout<<"打开文件失败"<<endl;
	}
	double x=0,y=0,z=0;
	int i=0;
	while(!feof(fp))
	{
		fscanf(fp,"%lf	%lf	%lf",&x,&y,&z);
		m_points->InsertPoint(i,x,y,z);
		m_vertices->InsertNextCell(1);
		m_vertices->InsertCellPoint(i);
		i++;
	}
	fclose(fp);
	vtkSmartPointer<vtkPolyData>m_polydata=vtkSmartPointer<vtkPolyData>::New();
	m_polydata->SetPoints(m_points);
	m_polydata->SetVerts(m_vertices);

	vtkSmartPointer<vtkCleanPolyData>m_cleanpoint=vtkSmartPointer<vtkCleanPolyData>::New();
	m_cleanpoint->SetInput(m_polydata);
	m_cleanpoint->SetTolerance(0.05);        //设置采样度
	m_cleanpoint->Update();

	vtkSmartPointer<vtkPolyDataWriter>pointcaiyang=vtkSmartPointer<vtkPolyDataWriter>::New();
	pointcaiyang->SetFileName("采样后的点.vtp");    //保存采样后的结果
	pointcaiyang->SetInput(m_cleanpoint->GetOutput());
	pointcaiyang->Write();

	vtkSmartPointer<vtkPolyDataMapper>caiyangqian=vtkSmartPointer<vtkPolyDataMapper>::New();
	caiyangqian->SetInputConnection(m_polydata->GetProducerPort());

	vtkSmartPointer<vtkPolyDataMapper>caiyanghou=vtkSmartPointer<vtkPolyDataMapper>::New();
	caiyanghou->SetInputConnection(m_cleanpoint->GetOutputPort());

	vtkSmartPointer<vtkActor>actor_qian=vtkSmartPointer<vtkActor>::New();
	actor_qian->SetMapper(caiyangqian);
	actor_qian->GetProperty()->SetColor(1,0,0);
//	actor_qian->GetProperty()->SetPointSize(2);

	vtkSmartPointer<vtkActor>actor_hou=vtkSmartPointer<vtkActor>::New();
	actor_hou->SetMapper(caiyanghou);
	actor_hou->GetProperty()->SetColor(0,1,0);
	actor_hou->GetProperty()->SetPointSize(4);

	ren1->AddActor(actor_qian);
	ren1->AddActor(actor_hou);

	ren1->SetBackground(1,1,1);
	renWin->AddRenderer(ren1);
	renWin->SetSize(600,600);
	iren->SetRenderWindow(renWin);
	renWin->Render();
	renWin->Start();
	getchar();
}

 

如果您觉得这篇博文有用，请访问我的个人站：http://www.stubbornhuang.com，更多博文干货等着您。