VTK—vtkQuadraticHexahedron细分

目录

1.提取单元格点和索引

2.用vtkUnstructuredGrid表达

3.显示效果,此时直边,细分无效

4.每个点位置添加球(字形)

5.折断直边,准备细分

6.细分过滤器vtkTessellatorFilter

7.使用vtkSliderWidget动态细分

8.完整代码


vtkQuadraticHexahedron是非线性的二次六面体单元格,就是普通的六面体每个边中间多了一个控制点。可以把点提取出来用vtkUnstructuredGrid显示。

1.提取单元格点和索引

cpp 复制代码
Create(vtkQuadraticHexahedron, hexadron);
    double *coord = hexadron->GetParametricCoords();
    int num = hexadron->GetNumberOfPoints();
    /*
     * 这两行无法获取数据
    hexadron->GetPoints();
    hexadron->GetPointIds();
    */

2.用vtkUnstructuredGrid表达

cpp 复制代码
for(int i = 0; i < num; i++)
    {
        double point[3] = {coord[i * 3 + 0], coord[i * 3 + 1], coord[i * 3 + 2]};
        points->InsertNextPoint(point);
        ids->InsertNextId(i);
    }
    Create(vtkUnstructuredGrid, grid);
    grid->SetPoints(points);
    grid->InsertNextCell(hexadron->GetCellType(), ids);

3.显示效果,此时直边,细分无效

4.每个点位置添加球(字形)

cpp 复制代码
Create(vtkPolyData, polydata);
    polydata->SetPoints(points);
    Create(vtkSphereSource, sphere);
    sphere->SetRadius(0.02);

    Create(vtkGlyph3D, glyph);
    glyph->SetInputData(polydata);
    glyph->SetSourceConnection(sphere->GetOutputPort());

    Create(vtkPolyDataMapper, gMapper);
    gMapper->SetInputConnection(glyph->GetOutputPort());

    Create(vtkActor, gActor);
    gActor->SetMapper(gMapper);

vtkGlyph3D有两个输入,setInputData()设置的其实点位置,setSourceConnection()设置的是要放在点位置的字形,此处是球。

5.折断直边,准备细分

cpp 复制代码
Create(vtkMinimalStandardRandomSequence, sequence);
    sequence->SetSeed(5);
    for(int i = 0; i < num; i++)
    {
        double point[3] = {coord[i * 3 + 0], coord[i * 3 + 1], coord[i * 3 + 2]};
        for(int j = 0; j < 3; j++)
        {
            point[j] = point[j] + sequence->GetRangeValue(-0.1, 0.1);
            sequence->Next();
        }
        points->InsertNextPoint(point);
        ids->InsertNextId(i);
    }

6.细分过滤器vtkTessellatorFilter

cpp 复制代码
Create(vtkTessellatorFilter, tessellator);
    tessellator->SetInputData(grid);
    tessellator->SetChordError(0.03);
    tessellator->Update();

    Create(vtkDataSetMapper, mapper);
    mapper->SetInputData(tessellator->GetOutput());

setChordError()能控制细分精细程度,目前还没理解参数含义(猜测是单元对角线的比例值)。

7.使用vtkSliderWidget动态细分

cpp 复制代码
Create(vtkSliderWidget, widget);
    widget->SetInteractor(inter);
    vtkSliderRepresentation3D *rep = (vtkSliderRepresentation3D *)widget->GetRepresentation();
    rep->GetPoint1Coordinate()->SetValue(-1, -1, 0);
    rep->GetPoint2Coordinate()->SetValue(2, -1, 0);
    rep->SetMaximumValue(0.1);
    rep->SetMinimumValue(0.001);
    rep->SetValue(tessellator->GetChordError());
    widget->SetEnabled(true);


    Create(myCallback, callback);
    callback->m_tessellator = tessellator;
    callback->m_mapper = mapper;
    callback->m_actor = actor;
    widget->AddObserver(vtkCommand::InteractionEvent, callback);

8.完整代码

cpp 复制代码
#include<vtkIdList.h>
#include<vtkPoints.h>
#include<vtkGlyph3D.h>
#include<vtkProperty.h>
#include<vtkSliderWidget.h>
#include<vtkSphereSource.h>
#include<vtkDataSetMapper.h>
#include<vtkUnstructuredGrid.h>
#include<vtkTessellatorFilter.h>
#include<vtkQuadraticHexahedron.h>
#include<vtkSliderRepresentation3D.h>
#include<vtkMinimalStandardRandomSequence.h>
#include<vtkCommand.h>


#include<vtkNew.h>
#include<vtkRenderer.h>
#include<vtkRenderWindow.h>
#include<vtkPolyDataMapper.h>
#include<vtkRenderWindowInteractor.h>

#define Create(type,name)    vtkNew<type> name


#include<vtkAutoInit.h>

VTK_MODULE_INIT(vtkRenderingOpenGL2)
VTK_MODULE_INIT(vtkInteractionStyle)
VTK_MODULE_INIT(vtkRenderingFreeType)


//用于交互回调的类
class myCallback: public vtkCommand
{
public:
    static myCallback *New() {
        return new myCallback;
    }

    vtkTypeMacro(myCallback, vtkCommand);

    void Execute(vtkObject *caller, unsigned long eventId, void *callData)override {
        vtkSliderWidget *slider = (vtkSliderWidget *)caller;
        m_tessellator->SetChordError(slider->GetSliderRepresentation()->GetValue());
        m_tessellator->Update();

    }

    vtkTessellatorFilter *m_tessellator;


private:
    myCallback() {
        m_tessellator = nullptr;

    }
    ~myCallback() {

    }
};

int main()
{

    Create(vtkQuadraticHexahedron, hexadron);
    double *coord = hexadron->GetParametricCoords();
    int num = hexadron->GetNumberOfPoints();
    /*
     * 这两行无法获取数据
    hexadron->GetPoints();
    hexadron->GetPointIds();
    */

    Create(vtkPoints, points);
    Create(vtkIdList, ids);
    
    //最常用的随机类
    Create(vtkMinimalStandardRandomSequence, sequence);
    sequence->SetSeed(5);
    for(int i = 0; i < num; i++)
    {
        double point[3] = {coord[i * 3 + 0], coord[i * 3 + 1], coord[i * 3 + 2]};
        for(int j = 0; j < 3; j++)
        {
            point[j] = point[j] + sequence->GetRangeValue(-0.1, 0.1);
            sequence->Next();
        }
        points->InsertNextPoint(point);
        ids->InsertNextId(i);
    }
    
    
    Create(vtkUnstructuredGrid, grid);
    grid->SetPoints(points);
    grid->InsertNextCell(hexadron->GetCellType(), ids);
    
    //用简式逼近非线性有限元单元
    Create(vtkTessellatorFilter, tessellator);
    tessellator->SetInputData(grid);
    tessellator->SetChordError(0.03);
    tessellator->Update();

    Create(vtkDataSetMapper, mapper);
    mapper->SetInputData(tessellator->GetOutput());


    Create(vtkActor, actor);
    actor->SetMapper(mapper);
    actor->GetProperty()->SetEdgeVisibility(true);
    actor->GetProperty()->SetRepresentationToSurface();

    Create(vtkPolyData, polydata);
    polydata->SetPoints(points);
    
    //放在每个点位置的字形(球)
    Create(vtkSphereSource, sphere);
    sphere->SetRadius(0.02);
    
    //将字形(球)放在每个点位置
    Create(vtkGlyph3D, glyph);
    glyph->SetInputData(polydata);
    glyph->SetSourceConnection(sphere->GetOutputPort());

    Create(vtkPolyDataMapper, gMapper);
    gMapper->SetInputConnection(glyph->GetOutputPort());

    Create(vtkActor, gActor);
    gActor->SetMapper(gMapper);

    Create(vtkRenderer, render);
    render->AddActor(actor);
    render->AddActor(gActor);

    Create(vtkRenderWindow, window);
    window->AddRenderer(render);

    Create(vtkRenderWindowInteractor, inter);
    inter->SetRenderWindow(window);
    
    
    //滑动条
    Create(vtkSliderWidget, widget);
    widget->SetInteractor(inter);
    vtkSliderRepresentation3D *rep = (vtkSliderRepresentation3D *)widget->GetRepresentation();
    rep->GetPoint1Coordinate()->SetValue(-1, -1, 0);
    rep->GetPoint2Coordinate()->SetValue(2, -1, 0);
    rep->SetMaximumValue(0.1);
    rep->SetMinimumValue(0.001);
    rep->SetValue(tessellator->GetChordError());
    widget->SetEnabled(true);
    
    
    //交互回调
    Create(myCallback, callback);
    callback->m_tessellator = tessellator;

    widget->AddObserver(vtkCommand::InteractionEvent, callback);

    inter->Start();


    return 0;
}
相关推荐
捕鲸叉1 小时前
创建线程时传递参数给线程
开发语言·c++·算法
A charmer1 小时前
【C++】vector 类深度解析:探索动态数组的奥秘
开发语言·c++·算法
Peter_chq1 小时前
【操作系统】基于环形队列的生产消费模型
linux·c语言·开发语言·c++·后端
青花瓷2 小时前
C++__XCode工程中Debug版本库向Release版本库的切换
c++·xcode
幺零九零零4 小时前
【C++】socket套接字编程
linux·服务器·网络·c++
捕鲸叉4 小时前
MVC(Model-View-Controller)模式概述
开发语言·c++·设计模式
Dola_Pan5 小时前
C++算法和竞赛:哈希算法、动态规划DP算法、贪心算法、博弈算法
c++·算法·哈希算法
yanlou2335 小时前
KMP算法,next数组详解(c++)
开发语言·c++·kmp算法
小林熬夜学编程5 小时前
【Linux系统编程】第四十一弹---线程深度解析:从地址空间到多线程实践
linux·c语言·开发语言·c++·算法
阿洵Rain6 小时前
【C++】哈希
数据结构·c++·算法·list·哈希算法