VTK Learning Twenty - RectilinearGridToTetrahedra

VTK Learning Twenty - RectilinearGridToTetrahedra

Description

規則不等間距格網體轉換成四面體不規則網格體

Code

#include <vtkSmartPointer.h>
#include <vtkRectilinearGrid.h>
#include <vtkRectilinearGridToTetrahedra.h>
#include <vtkUnstructuredGrid.h>
#include <vtkMath.h>
#include <vtkDoubleArray.h>
#include <vtkFloatArray.h>
#include <vtkPoints.h>
#include <vtkPointData.h>
#include <vtkDataSetMapper.h>
#include <vtkActor.h>
#include <vtkRenderWindow.h>
#include <vtkRenderer.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkProperty.h>

int main(int, char *[])
{
	// Create a grid
	vtkSmartPointer<vtkRectilinearGrid> grid =
		vtkSmartPointer<vtkRectilinearGrid>::New();

	grid->SetDimensions(5, 4, 3);

	vtkSmartPointer<vtkDoubleArray> xArray =
		vtkSmartPointer<vtkDoubleArray>::New();
	xArray->InsertNextValue(0.0);
	xArray->InsertNextValue(1.0);
	xArray->InsertNextValue(2.0);
	xArray->InsertNextValue(3.0);
	xArray->InsertNextValue(4.0);

	vtkSmartPointer<vtkDoubleArray> yArray =
		vtkSmartPointer<vtkDoubleArray>::New();
	yArray->InsertNextValue(0.0);
	yArray->InsertNextValue(1.0);
	yArray->InsertNextValue(2.0);
	yArray->InsertNextValue(3.0);

	vtkSmartPointer<vtkDoubleArray> zArray =
		vtkSmartPointer<vtkDoubleArray>::New();
	zArray->InsertNextValue(0.0);
	zArray->InsertNextValue(1.0);
	zArray->InsertNextValue(2.0);

	grid->SetXCoordinates(xArray);
	grid->SetYCoordinates(yArray);
	grid->SetZCoordinates(zArray);
	

	vtkSmartPointer<vtkRectilinearGridToTetrahedra> rectilinearGridToTetrahedra =
		vtkSmartPointer<vtkRectilinearGridToTetrahedra>::New();
	rectilinearGridToTetrahedra->SetInputData(grid);
	rectilinearGridToTetrahedra->Update();
	//設置屬性
	vtkSmartPointer<vtkPoints>points;
	points=rectilinearGridToTetrahedra->GetOutput()->GetPoints();
	vtkSmartPointer<vtkFloatArray>values = vtkSmartPointer<vtkFloatArray>::New();
	for (vtkIdType id = 0; id < points->GetNumberOfPoints(); id++) {
		values->InsertNextValue(vtkMath::Random(0, 2));
	}
	rectilinearGridToTetrahedra->GetOutput()->GetPointData()->SetScalars(values);
  // Create a mapper and actor
  vtkSmartPointer<vtkDataSetMapper> mapper = 
    vtkSmartPointer<vtkDataSetMapper>::New();
  mapper->SetInputConnection(rectilinearGridToTetrahedra->GetOutputPort());

  vtkSmartPointer<vtkActor> actor = 
    vtkSmartPointer<vtkActor>::New();
  actor->SetMapper(mapper);
  actor->GetProperty()->SetEdgeVisibility(true);
  // Visualize
  vtkSmartPointer<vtkRenderer> renderer = 
    vtkSmartPointer<vtkRenderer>::New();
  vtkSmartPointer<vtkRenderWindow> renderWindow = 
    vtkSmartPointer<vtkRenderWindow>::New();
  renderWindow->AddRenderer(renderer);
  vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor = 
    vtkSmartPointer<vtkRenderWindowInteractor>::New();
  renderWindowInteractor->SetRenderWindow(renderWindow);

  renderer->AddActor(actor);

  renderWindow->Render();
  renderWindowInteractor->Start();

  return EXIT_SUCCESS;
}

Result