ParaView热流图：如何计算单元格Cell的体积-卡核

1.Python Calculator过滤器

表达式(expression)：

volume(inputs[0])

Array Association:

Cell Data

在这里插入图片描述

2.vtkCellSizeFilter

@brief计算单元格大小。
计算VTK中所有类型单元格的单元格大小。对于三角形、四边形、四角形和六角形，使用vtkMeshQuality中的静态方法。
为了更高的准确性，这是通过Verdict来实现的。
其他类型的Cell在可能的情况下单独进行分析，在不可能的情况下分解成三角形或四角形。当细胞分裂成三角形或四角形时，精确度可能会降低。默认情况下，会计算所有大小（但可以选择忽略）：顶点计数、长度、面积和体积单元。每个请求的大小都使用单独的数组（例如，如果请求长度和体积，则将从该筛选器输出两个数组）。这4个数组可以单独命名，默认值为VertexCount、Length、Area和Volme。对于未计算大小的单元格尺寸，将给出0的值。
对于应该计算大小但无法计算的单元格，过滤器将返回-1。
ComputeSum选项将对单元格大小（不包括重影单元格）求和，并将值放入以相应单元格数据数组名称命名的vtkFieldData数组中。对于复合数据集，所有数据块的总和也将添加到顶层数据块的字段数据中进行总和。

相关过滤器：vtkPCellSizeFilter

std::string inputname = "test.vtu";
	vtkSmartPointer<vtkXMLUnstructuredGridReader> inputReader = vtkSmartPointer<vtkXMLUnstructuredGridReader>::New();
	inputReader->SetFileName(inputname.c_str());
	inputReader->Update();
	
	vtkNew<vtkCellSizeFilter> filter;	
	filter->SetInputConnection(inputReader->GetOutputPort());
	//filter->ComputeSumOn();
	//filter->ComputeVolumeOn();
	filter->Update();
	vtkDoubleArray* vertexCount = vtkDoubleArray::SafeDownCast(
		vtkUnstructuredGrid::SafeDownCast(filter->GetOutput())->GetCellData()->GetArray("VertexCount"));
	vtkDoubleArray* length = vtkDoubleArray::SafeDownCast(
		vtkUnstructuredGrid::SafeDownCast(filter->GetOutput())->GetCellData()->GetArray("Length"));
	vtkDoubleArray* area = vtkDoubleArray::SafeDownCast(
		vtkUnstructuredGrid::SafeDownCast(filter->GetOutput())->GetCellData()->GetArray("Area"));
	vtkDoubleArray* volume = vtkDoubleArray::SafeDownCast(
		vtkUnstructuredGrid::SafeDownCast(filter->GetOutput())->GetCellData()->GetArray("Volume"));

	vtkSmartPointer<vtkUnstructuredGrid> unstructuredGrid= vtkUnstructuredGrid::SafeDownCast(filter->GetOutput());
	vtkCellData* cell = unstructuredGrid->GetCellData();
	vtkDoubleArray* volume = vtkDoubleArray::SafeDownCast(cell->GetArray("Volume"));

	vtkSmartPointer <vtkXMLUnstructuredGridWriter> polyWriter = vtkSmartPointer <vtkXMLUnstructuredGridWriter>::New();
	polyWriter->SetFileName("mycellsize.vtu");
	polyWriter->SetInputData(unstructuredGrid);
	polyWriter->SetDataModeToAscii();
	polyWriter->Write();