VTK读取外部数据并显示

最新推荐文章于 2024-08-15 09:49:55 发布

原创最新推荐文章于 2024-08-15 09:49:55 发布 · 2.9k 阅读

CC 4.0 BY-SA版权

文章标签：

vtk中用vtkTecplotReader类读取外部文件数据（该文件数据可以用Tecplot打开并显示）。
主要流程是：
（1）vtkTecplotReader读取外部文件，
（2）vtkAppendFilter数据过滤，
（3）vtkUnstructuredGrid过滤后数据转换非结构网格，
（4）读取存储在cell单元中的结果数据：
vtkSmartPointer ds = reader->GetOutput();
vtkDataSetAttributes* cellData = ds->GetBlock(0)->GetAttributes(vtkDataSet::CELL);
通过 int n = cellData->GetNumberOfArrays();可以获取cell单元中存储的标量结果数量，知道有多少个标量数据。
（5）将数据存入polyData
vtkSmartPointer polyData=vtkSmartPointer::New();
polyData->SetPoints(meshPoints);
polyData->SetPolys(meshCells);
（6）将标量数据设置到polyData：
vtkDataArray* scalar = cellData->GetScalars(arrayName);
polyData->GetCellData()->SetScalars(scalar);
（7）将cell单元标量数据转换为point点标量数据：
vtkSmartPointer cellDataToPointData = vtkSmartPointer::New();
cellDataToPointData->SetInputData(polyData);
cellDataToPointData->Update();
vtkDataArray* detail = cellDataToPointData->GetOutput()->GetPointData()->GetArray(0);
polyData->GetPointData()->SetScalars(detail);
（8）按照顺序定义和设置Mapper，Actor，Renderer，RenderWindow，RenderWindowInteractor进行调用vtk原生窗口进行显示。

#include <vtkAutoInit.h>
VTK_MODULE_INIT(vtkRenderingOpenGL2);// VTK was built with vtkRenderingOpenGL2
VTK_MODULE_INIT(vtkInteractionStyle);
VTK_MODULE_INIT(vtkRenderingFreeType)

#include "vtkSmartPointer.h"
#include "vtkMultiBlockDataSet.h"
#include "vtkTecplotReader.h"
#include "vtkPointData.h"
#include "vtkCellData.h"
#include "vtkDataSetMapper.h"
#include "vtkPolyDataMapper.h"
#include "vtkActor.h"
#include "vtkRenderWindow.h"
#include "vtkRenderer.h"
#include "vtkCellDataToPointData.h"
#include "vtkAppendFilter.h"

bool ::showContourData(QString& fileName, QString& selectAngle)
{
  QString filePath = areaFilePath + "/" + fileName;
  //支持中文路径
  QTextCodec::setCodecForLocale(QTextCodec::codecForName("GBK"));
  QByteArray ba = filePath.toLocal8Bit();
  const char* fname = ba.data();

  vtkSmartPointer<vtkPolyDataMapper> map1 = vtkSmartPointer<vtkPolyDataMapper>::New();
  vtkSmartPointer<vtkActor> actor1 = vtkSmartPointer<vtkActor>::New();
  vtkSmartPointer<vtkRenderer> ren = vtkSmartPointer<vtkRenderer>::New();
  vtkSmartPointer<vtkRenderWindow> win = vtkSmartPointer<vtkRenderWindow>::New();
  vtkSmartPointer<vtkTecplotReader> reader = vtkSmartPointer<vtkTecplotReader>::New();
  
  reader->SetFileName(fname);
  reader->Update();

  vtkSmartPointer<vtkMultiBlockDataSet> ds = reader->GetOutput();
  
  vtkSmartPointer<vtkAppendFilter> filter = vtkSmartPointer<vtkAppendFilter>::New();
  filter->AddInputData(reader->GetOutput()->GetBlock(0));
  filter->Update();
  
  vtkSmartPointer<vtkUnstructuredGrid> grid = vtkSmartPointer<vtkUnstructuredGrid>::New();
  grid=filter->GetOutput();
  
  int test1 = grid->GetNumberOfCells();
  int test2 = grid->GetNumberOfPoints();
  vtkPoints* meshPoints = grid->GetPoints();
  vtkCellArray* meshCells = grid->GetCells();

  vtkSmartPointer<vtkPolyData> polyData=vtkSmartPointer<vtkPolyData>::New();
  polyData->SetPoints(meshPoints);
  polyData->SetPolys(meshCells);

  vtkDataSetAttributes* cellData = ds->GetBlock(0)->GetAttributes(vtkDataSet::CELL);
  int n = cellData->GetNumberOfArrays();

  //获取cell单元中的标量数据
  const char* arrayName;
  QByteArray btArray = selectAngle.toUtf8();
  arrayName = btArray.data();
  vtkDataArray* scalar = cellData->GetScalars(arrayName);

  vtkIdType number = scalar->GetNumberOfValues();
  std::cout << " " << std::endl;
  polyData->GetCellData()->SetScalars(scalar);

  //将cell中的数据转换到point上
  vtkSmartPointer<vtkCellDataToPointData> cellDataToPointData = vtkSmartPointer<vtkCellDataToPointData>::New();
  cellDataToPointData->SetInputData(polyData);
  cellDataToPointData->Update();
  vtkDataArray* detail1 = cellDataToPointData->GetOutput()->GetPointData()->GetArray(0);
  polyData->GetPointData()->SetScalars(detail1);

  map1->SetInputData(polyData);
  map1->SetScalarModeToUsePointData();
  map1->SetScalarRange(scalar->GetRange());
  map1->Update();

  vtkNew<vtkLookupTable> lookupTable;
  lookupTable->SetTableRange(scalar->GetRange());
  lookupTable->Build();
  map1->SetLookupTable(lookupTable);

  vtkNew<vtkScalarBarActor> scalarBarActor;
  scalarBarActor->SetLookupTable(lookupTable);

  actor1->SetMapper(map1);
  actor1->GetProperty()->SetColor(255, 0, 0);

  win->AddRenderer(ren);

  vtkSmartPointer<vtkRenderWindowInteractor> iren = vtkSmartPointer<vtkRenderWindowInteractor>::New();
  iren->SetRenderWindow(win);
  ren->AddActor(actor1);
  ren->AddActor(scalarBarActor);
  win->SetSize(600, 600);
  win->Render();
  iren->Start();

  return true;
}

新的改变

我们对Markdown编辑器进行了一些功能拓展与语法支持，除了标准的Markdown编辑器功能，我们增加了如下几点新功能，帮助你用它写博客：

全新的界面设计 ，将会带来全新的写作体验；
在创作中心设置你喜爱的代码高亮样式，Markdown 将代码片显示选择的高亮样式 进行展示；
增加了 图片拖拽 功能，你可以将本地的图片直接拖拽到编辑区域直接展示；
全新的 KaTeX数学公式 语法；
增加了支持甘特图的mermaid语法¹ 功能；
增加了 多屏幕编辑 Markdown文章功能；
增加了 焦点写作模式、预览模式、简洁写作模式、左右区域同步滚轮设置 等功能，功能按钮位于编辑区域与预览区域中间；
增加了 检查列表 功能。

功能快捷键

撤销：Ctrl/Command + Z
重做：Ctrl/Command + Y
加粗：Ctrl/Command + B
斜体：Ctrl/Command + I
标题：Ctrl/Command + Shift + H
无序列表：Ctrl/Command + Shift + U
有序列表：Ctrl/Command + Shift + O
检查列表：Ctrl/Command + Shift + C
插入代码：Ctrl/Command + Shift + K
插入链接：Ctrl/Command + Shift + L
插入图片：Ctrl/Command + Shift + G
查找：Ctrl/Command + F
替换：Ctrl/Command + G

合理的创建标题，有助于目录的生成

直接输入1次#，并按下space后，将生成1级标题。
输入2次#，并按下space后，将生成2级标题。
以此类推，我们支持6级标题。有助于使用TOC语法后生成一个完美的目录。

如何改变文本的样式

强调文本 强调文本

加粗文本 加粗文本

标记文本

~~删除文本~~

引用文本

H₂O is是液体。

2¹⁰ 运算结果是 1024.

插入链接与图片

链接: link.

图片: Alt

带尺寸的图片:

居中的图片: Alt

居中并且带尺寸的图片:

当然，我们为了让用户更加便捷，我们增加了图片拖拽功能。

如何插入一段漂亮的代码片

去博客设置页面，选择一款你喜欢的代码片高亮样式，下面展示同样高亮的 代码片.

// An highlighted block
var foo = 'bar';

生成一个适合你的列表

项目
- 项目
  - 项目

项目1
项目2
项目3

计划任务
完成任务

创建一个表格

一个简单的表格是这么创建的：

项目	Value
电脑	$1600
手机	$12
导管	$1

设定内容居中、居左、居右

使用:---------:居中
使用:----------居左
使用----------:居右

第一列	第二列	第三列
第一列文本居中	第二列文本居右	第三列文本居左

SmartyPants

SmartyPants将ASCII标点字符转换为“智能”印刷标点HTML实体。例如：

TYPE	ASCII	HTML
Single backticks	`'Isn't this fun?'`	‘Isn’t this fun?’
Quotes	`"Isn't this fun?"`	“Isn’t this fun?”
Dashes	`-- is en-dash, --- is em-dash`	– is en-dash, — is em-dash