VTK Learning Twenty - RectilinearGridToTetrahedra

最新推荐文章于 2024-08-23 21:47:54 发布

原创最新推荐文章于 2024-08-23 21:47:54 发布 · 263 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#RectilinearGrid #Tetrahedra

VTK 专栏收录该内容

72 篇文章

订阅专栏

本文介绍如何使用VTK将规则的不等间距三维格网转换为四面体不规则网格，并通过代码示例展示了从创建格网到可视化整个过程。关键步骤包括设置坐标、转换网格类型、添加随机属性以及渲染输出。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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

在这里插入图片描述