vtkQuadric创建二次曲面

在本实例中，我们将用到vtkQuadric、vtkSampleFunction、vtkContourFilter三个类，分别是二次曲面函数、函数曲面抽样和等高滤波。

vtkQuadric负责二次曲面基本参数的设置，vtkSampleFunction则是对二次曲面进行等间隔逐点采样，在本例中采样点数为50*50*50。vtkContourFilter负责将采集到空间点，转成vtkPolyData型对象。

#ifndef INITIAL_OPENGL

#define INITIAL_OPENGL

#include <vtkAutoInit.h>

VTK_MODULE_INIT(vtkRenderingOpenGL)

VTK_MODULE_INIT(vtkInteractionStyle)

VTK_MODULE_INIT(vtkRenderingFreeType)

#endif

#include <iostream>

using namespace std;

#include <vtkSmartPointer.h>

#include <vtkQuadric.h>

#include <vtkSampleFunction.h>

#include <vtkContourFilter.h>

#include <vtkActor.h>

#include <vtkProperty.h>

#include <vtkRenderWindow.h>

#include <vtkRenderer.h>

#include <vtkRenderWindowInteractor.h>

#include <vtkPolyDataMapper.h>

#include <vtkCellArray.h>

#include <vtkImageDataGeometryFilter.h>

#include <vtkImageData.h>

#include <vtkImageActor.h>

int main()

{

    //创建二次函数

    vtkSmartPointer<vtkQuadric>quadric=vtkSmartPointer<vtkQuadric>::New();

    quadric->SetCoefficients(0.5,,0.2,,0.1,,,0.2,,);

    //二次函数采样分辨率

    vtkSmartPointer<vtkSampleFunction>sample=vtkSmartPointer<vtkSampleFunction>::New();

    sample->SetSampleDimensions(,,);

    sample->SetImplicitFunction(quadric);

    double xmin=,xmax=,ymin=,ymax=,zmin=,zmax=;

    sample->SetModelBounds(xmin,xmax,ymin,ymax,zmin,zmax);

    vtkSmartPointer<vtkContourFilter> contourFilter=vtkSmartPointer<vtkContourFilter>::New();

    contourFilter->SetInputConnection(sample->GetOutputPort());

    contourFilter->GenerateValues(,,);

    contourFilter->Update();
/**在此处可以有两种方法，使得contourFilter得到的vtkPolyData型对象与mapper对接
1.用contourFilter->GetOutput();得到PolyData，然后输入mapper；
2.再有就是，直接调用contourFilter->GetOutputPort()输入到 mapper->SetInputConnection（）
******************************************************************************************/

vtkSmartPointer<vtkPolyData> polyData= contourFilter->GetOutput();

//可视化,设置mapper、actor

    vtkSmartPointer<vtkPolyDataMapper>mapper=vtkSmartPointer<vtkPolyDataMapper>::New();

//    mapper->SetInputConnection(contourFilter->GetOutputPort());

    mapper->SetInputData(polyData);

    vtkSmartPointer<vtkActor>actor=vtkSmartPointer<vtkActor>::New();

    actor->SetMapper(mapper);

    actor->SetOrigin(,,);

    //设置renderer、window和interactor

    vtkSmartPointer<vtkRenderer> renderer=vtkSmartPointer<vtkRenderer>::New();

    vtkSmartPointer<vtkRenderWindow> renWin=vtkSmartPointer<vtkRenderWindow>::New();

    vtkSmartPointer<vtkRenderWindowInteractor> iren=vtkSmartPointer<vtkRenderWindowInteractor>::New();

    renWin->AddRenderer(renderer);

    iren->SetRenderWindow(renWin);

    renderer->AddActor(actor);

//    renderer->AddActor(imageActor);

    renderer->SetBackground(0.3,0.6,0.3);

    renWin->Render();

    iren->Start();

    return ;

}

#ifndefINITIAL_OPENGL

#defineINITIAL_OPENGL

#include<vtkAutoInit.h>

VTK_MODULE_INIT(vtkRenderingOpenGL)

VTK_MODULE_INIT(vtkInteractionStyle)

VTK_MODULE_INIT(vtkRenderingFreeType)

#endif

#include<iostream>

usingnamespacestd;

#include<vtkSmartPointer.h>

#include<vtkQuadric.h>

#include<vtkSampleFunction.h>

#include<vtkContourFilter.h>

#include<vtkActor.h>

#include<vtkProperty.h>

#include<vtkRenderWindow.h>

#include<vtkRenderer.h>

#include<vtkRenderWindowInteractor.h>

#include<vtkPolyDataMapper.h>

#include<vtkCellArray.h>

#include<vtkImageDataGeometryFilter.h>

#include<vtkImageData.h>

#include<vtkImageActor.h>

intmain()

{

//创建二次函数

vtkSmartPointer<vtkQuadric>quadric=vtkSmartPointer<vtkQuadric>::New();

quadric->SetCoefficients(0.5,,0.2,,0.1,,,0.2,,);

//二次函数采样分辨率

vtkSmartPointer<vtkSampleFunction>sample=vtkSmartPointer<vtkSampleFunction>::New();

sample->SetSampleDimensions(,,);

sample->SetImplicitFunction(quadric);

doublexmin=,xmax=,ymin=,ymax=,zmin=,zmax=;

sample->SetModelBounds(xmin,xmax,ymin,ymax,zmin,zmax);

vtkSmartPointer<vtkContourFilter>contourFilter=vtkSmartPointer<vtkContourFilter>::New();

contourFilter->SetInputConnection(sample->GetOutputPort());

contourFilter->GenerateValues(,,);

contourFilter->Update();

vtkSmartPointer<vtkPolyData>polyData=contourFilter->GetOutput();

//可视化,设置mapper、actor

vtkSmartPointer<vtkPolyDataMapper>mapper=vtkSmartPointer<vtkPolyDataMapper>::New();

//mapper->SetInputConnection(contourFilter->GetOutputPort());

mapper->SetInputData(polyData);

vtkSmartPointer<vtkActor>actor=vtkSmartPointer<vtkActor>::New();

actor->SetMapper(mapper);

actor->SetOrigin(,,);

//设置renderer、window和interactor

vtkSmartPointer<vtkRenderer>renderer=vtkSmartPointer<vtkRenderer>::New();

vtkSmartPointer<vtkRenderWindow>renWin=vtkSmartPointer<vtkRenderWindow>::New();

vtkSmartPointer<vtkRenderWindowInteractor>iren=vtkSmartPointer<vtkRenderWindowInteractor>::New();

renWin->AddRenderer(renderer);

iren->SetRenderWindow(renWin);

renderer->AddActor(actor);

//renderer->AddActor(imageActor);

renderer->SetBackground(0.3,0.6,0.3);

renWin->Render();

iren->Start();

return;

}