VTK 面绘制三维数据源

最新推荐文章于 2024-10-06 21:43:01 发布

青青小禾

最新推荐文章于 2024-10-06 21:43:01 发布

阅读量666

点赞数

文章标签： VTK

本文通过使用VTK库和Marching Cubes算法实现了三维图像数据集的等值面提取及可视化。该程序创建了一个包含特定分布像素值的图像数据集，并通过Marching Cubes算法提取了这些像素值的等值面，最后利用VTK进行三维渲染。

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

// RawResolution.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"

#include<iostream>
using namespace std;

#include "vtkUnsignedCharArray.h"
#include "vtkPiecewiseFunction.h"
#include "vtkColorTransferFunction.h"
#include "vtkVolumeRayCastCompositeFunction.h"
#include "vtkVolumeRayCastMapper.h"
#include "vtkImageData.h"
#include "vtkVolumeProperty.h"
#include "vtkVolume.h"
#include "vtkRenderWindow.h"
#include "vtkRenderer.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkUnsignedShortArray.h"
#include "vtkPointData.h"
#include "vtkVolume16Reader.h"
#include "vtkMarchingCubes.h"

#include "vtkStripper.h"
#include "vtkPolyDataMapper.h"
#include "vtkProperty.h"
#include "vtkCamera.h"
#include "vtkBoxWidget.h"

#include "vtkImageReader.h"
#include "vtkVectorNorm.h"
#include "vtkDataSetMapper.h"

int _tmain(int argc, _TCHAR* argv[])
{

vtkImageData *image = vtkImageData::New();
image->SetDimensions(30, 30, 30);
image->SetSpacing(5, 5, 5);
image->SetOrigin(0, 0, 0);
image->SetScalarTypeToInt();
image->SetNumberOfScalarComponents(1);
image->AllocateScalars();
int *ptr = (int *)image->GetScalarPointer();

for(int i=0; i<30*30*30; i++)
{
if(i%9 == 0)
*ptr++ = 0;
else
*ptr++ = 2;
}

vtkMarchingCubes *iso = vtkMarchingCubes::New();
iso->SetInput(image);
iso->SetNumberOfContours(1);
iso->SetValue(0,1);
iso->ComputeGradientsOn();
iso->ComputeNormalsOff();
iso->ComputeScalarsOff();

vtkVectorNorm *gradiant = vtkVectorNorm::New();
gradiant->SetInputConnection(iso->GetOutputPort());

vtkDataSetMapper *cubeMapper = vtkDataSetMapper::New();
cubeMapper->SetInputConnection(gradiant->GetOutputPort());
vtkActor *cubeActor = vtkActor::New();
cubeActor->SetMapper(cubeMapper);

vtkCamera *camera = vtkCamera::New();
camera->SetPosition(1,1,1);
camera->SetFocalPoint(0,0,0);

vtkRenderer *renderer = vtkRenderer::New();
vtkRenderWindow *renWin = vtkRenderWindow::New();
renWin->AddRenderer(renderer);

vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
iren->SetRenderWindow(renWin);

renderer->AddActor(cubeActor);
renderer->SetActiveCamera(camera);
renderer->ResetCamera();
renderer->SetBackground(1,1,1);

renWin->SetSize(300,300);
renWin->Render();
iren->Start();

image->Delete();
cubeMapper->Delete();
cubeActor->Delete();
camera->Delete();
renderer->Delete();
renWin->Delete();
iren->Delete();

system("pause");
return 0;
}