Можно ли сохранить векторное поле в центре тяжести тетраэдрической сетки, используя vtkStructuredGrid
? Я попробовал следующий код, но VTK (версия 8.1) жалуется на это, что, я думаю, связано с тем, что это векторное поле определено в центроиде клетки.
Warning: In c:\vtk\src\common\datamodel\vtkdataset.cxx, line 443
vtkUnstructuredGrid (0000020A4EC10D50): Point array with 3 components, has 137 tuples but there are only 64 points
Векторное поле определяется:
vtkSmartPointer<vtkUnstructuredGrid> uGrid = vtkSmartPointer<vtkUnstructuredGrid>::New();
// ... I already populated uGrid with points and tetrahedral information
// numberOfPoints = 64
// numberOfTetrahedra = 103
// Add a vector field at the centroid of each tetrahedral
vtkSmartPointer<vtkDoubleArray> vectors = vtkSmartPointer<vtkDoubleArray>::New();
vectors->SetNumberOfTuples(numberOfTetrahedra);
vectors->SetNumberOfComponents(3);
for (vtkIdType ielement = 0; ielement < numberOfTetrahedra; ielement++)
{
vectors->InsertNextValue(vec[3 * ielement]);
vectors->InsertNextValue(vec[3 * ielement + 1]);
vectors->InsertNextValue(vec[3 * ielement + 2]);
}
uGrid->GetPointData()->SetVectors(vectors);
// Write the data to file
vtkSmartPointer<vtkXMLUnstructuredGridWriter> writer = vtkSmartPointer<vtkXMLUnstructuredGridWriter>::New();
writer->SetFileName("vtk_test_write_unstructured_grid.vtu");
writer->SetInputData(uGrid);
writer->Write();
Буду признателен за любые советы / предложения для решения этой проблемы.
Для данных, связанных с клетками, вы должны использовать vtkCellData. Так что в вашем коде попробуйте заменить строку, где вы установите вектор
uGrid->GetCellData()->SetVectors(vectors);
Эти векторы связаны с клеткой. Вы можете интерпретировать их как связанные с клеточным центроидом, хотя в .vtu
В файле вы не можете найти никакой такой явной связи с центроидами клеток.
Чтобы отобразить ячейки с центрированием полей в точки сетки, рассмотрите возможность vtkCellDataToPointData учебный класс.
Других решений пока нет …