Как написать скалярное поле с VTK в C ++? - PullRequest
Новая
       12

Как написать скалярное поле с VTK в C ++?

1 голос
/ 23 марта 2019

Я начинаю узнавать о ВТК. Я хотел бы написать скалярное поле в файле .vts. Структурированная сетка создается с помощью кода, приведенного ниже.

Я пытался найти информацию в руководстве пользователя и руководствах, но мне не удалось это сделать. 

// Create a grid
vtkSmartPointer<vtkStructuredGrid> structuredGrid = vtkSmartPointer<vtkStructuredGrid>::New();

vtkSmartPointer<vtkPoints> points = vtkSmartPointer<vtkPoints>::New();
unsigned int numi = 20;
unsigned int numj = 30;
unsigned int numk = 1;

for (unsigned int k = 0; k < numk; k++) {
    for (unsigned int j = 0; j < numj; j++) {
        for (unsigned int i = 0; i < numi; i++) {
            points->InsertNextPoint(i, j, k);
        }
    }
}

// Specify the dimensions of the grid
structuredGrid->SetDimensions(numi, numj, numk);
structuredGrid->SetPoints(points);

// Write file
vtkSmartPointer<vtkXMLStructuredGridWriter> writer = vtkSmartPointer<vtkXMLStructuredGridWriter>::New();
writer->SetFileName("output.vts");
writer->SetInputData(structuredGrid);
writer->Write();

1 Ответ

0 голосов
/ 24 марта 2019

Скаляр может быть установлен для атрибута PointData объекта vtkStructuredGrid.
Затем объект vtkXMLStructuredGridWriter записывает structdGrid в файл vts следующим образом: 

    vtkFloatArray *scalars = vtkFloatArray::New();
    for (int i=0; i< numi*numj*numk; i++) scalars->InsertTuple1(i, 1.0 * i /(numi*numj*numk));
    structuredGrid->GetPointData()->SetScalars( scalars );

    // Write file
    vtkSmartPointer<vtkXMLStructuredGridWriter> writer = vtkSmartPointer<vtkXMLStructuredGridWriter>::New();
    writer->SetFileName("output.vts");
    writer->SetInputData(structuredGrid);
    writer->Write();
...