给VTU文件添加网格度量属性后在ParaView中无法显示的问题排查

我尝试给VTU文件添加一个网格度量属性，该属性对应网格中每个三角形的单个数值。当前VTU文件结构如下：

<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
    <UnstructuredGrid>
        <Piece NumberOfPoints="468" NumberOfCells="862">
            <Points>
                <DataArray NumberOfComponents="3" type="Float64" Format="binary">
                    *insert nodes coordinates*
                </DataArray>
            </Points>
            <Cells>
                <DataArray type="Int32" Name="connectivity" Format="binary">
                    *insert connectivity info*
                </DataArray>
                <DataArray type="Int32" Name="offsets" Format="binary">
                    * insert offset data*
                </DataArray>
                <DataArray type="Int32" Name="types" Format="binary">
                </DataArray>
                ****This is where I inserted my array
                <DataArray type="Float64" Name="meshmetric" Format="binary">
                    *insert mesh metric values*
                </DataArray>
            </Cells>
        </Piece>
    </UnstructuredGrid>
</VTKFile>

我将网格度量数组插入到<types>数组之后，但加载到ParaView中时该属性并未显示。怀疑可能的原因：

• 数组放置位置不正确
• 不清楚ParaView中查看该属性的方法
• 需要修改<offsets>或<types>数组的相关信息

注：未使用Python

问题解决方法 #

1. 修正数组放置位置 #

VTU格式中，<Cells>标签仅允许包含connectivity、offsets、types三个固定数组，单元格属性必须放在<CellData>标签内，该标签需置于<Cells>之后、</Piece>之前。

修改后的正确结构示例：

<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
    <UnstructuredGrid>
        <Piece NumberOfPoints="468" NumberOfCells="862">
            <Points>
                <DataArray NumberOfComponents="3" type="Float64" Format="binary">
                    *insert nodes coordinates*
                </DataArray>
            </Points>
            <Cells>
                <DataArray type="Int32" Name="connectivity" Format="binary">
                    *insert connectivity info*
                </DataArray>
                <DataArray type="Int32" Name="offsets" Format="binary">
                    * insert offset data*
                </DataArray>
                <DataArray type="Int32" Name="types" Format="binary">
                </DataArray>
            </Cells>
            <!-- 这里放置单元格属性 -->
            <CellData>
                <DataArray type="Float64" Name="meshmetric" Format="binary">
                    *insert mesh metric values*
                </DataArray>
            </CellData>
        </Piece>
    </UnstructuredGrid>
</VTKFile>

2. ParaView中查看属性的步骤 #

• 加载修改后的VTU文件，在左侧Pipeline Browser选中数据集
• 切换到右侧Properties面板，找到并展开Cell Data选项，即可看到meshmetric属性
• 右键点击该属性，选择Color by即可可视化网格度量值

3. 数据一致性检查 #

• 确保meshmetric数组的元素数量与<Piece>标签中的NumberOfCells值（此处为862）完全一致，每个三角形对应一个数值
• <offsets>和<types>数组无需修改，只需保证它们原本正确（types数组中每个元素应为5，对应VTK_TRIANGLE类型）

内容的提问来源于stack exchange，提问作者Alaa M. Ali