基于vtk非结构网格的eclipse网格模型显示

基于vtk中的 vtkHexahedron网格，绘制空间上8个顶点组成的网格的三维模型，最整个模型需要绘制每个vtkHexhedron网格，组成整体模型，是vtkUnstructuredGrid的显示方法。这里只放了规则模型的框架图和颜色，可以用于显示eclipse grid格式的油藏地质三维模型。

public static void DrawCoordXYZGrid(Point3D[,] coordXYZGrid, ref RenderWindowControl myRenderWindowControl)

    {

        int node_ni = coordXYZGrid.GetLength(0);

        int node_nj = coordXYZGrid.GetLength(1);

        int node_nk = coordXYZGrid.GetLength(2);



        vtkPoints points = vtkPoints.New();

        vtkFloatArray scalars = vtkFloatArray.New();         



        points.SetNumberOfPoints(node_ni * node_nj * node_nk);

        int seq = 0;

        for (int k = 0; k < node_nk; k++)

            for (int j = 0; j < node_nj; j++)

                for (int i = 0; i < node_ni; i++)

                {

                    points.InsertPoint(seq, coordXYZGrid[i, j, k].X, coordXYZGrid[i, j, k].Y, coordXYZGrid[i, j, k].Z);

                    scalars.InsertTuple1(seq, k);

                    //Console.WriteLine(“{0} {1} {2}”, coordXYZGrid[i, j, k].x, coordXYZGrid[i, j, k].y, coordXYZGrid[i, j, k].z);

                    seq += 1;

                }



        vtkHexahedron hexahedron;

        vtkUnstructuredGrid coordXYZUG = vtkUnstructuredGrid.New();

        coordXYZUG.Allocate(1, 1);

        coordXYZUG.SetPoints(points);

        coordXYZUG.GetPointData().SetScalars(scalars);

        seq = 0;

        for (int k = 0; k < node_nk -1; k++)

            for (int j = 0; j < node_nj - 1; j++)

                for (int i = 0; i < node_ni - 1; i++)

                {

                    seq = k * node_nj * node_ni + j * node_ni + i;

                    hexahedron = vtkHexahedron.New();

                    hexahedron.GetPointIds().SetId(0, seq);

                    hexahedron.GetPointIds().SetId(1, seq + 1);

                    hexahedron.GetPointIds().SetId(2, seq + node_ni + 1);

                    hexahedron.GetPointIds().SetId(3, seq + node_ni);

                    hexahedron.GetPointIds().SetId(4, seq + node_nj * node_ni);

                    hexahedron.GetPointIds().SetId(5, seq + node_nj * node_ni + 1);

                    hexahedron.GetPointIds().SetId(6, seq + node_nj * node_ni + node_ni + 1);

                    hexahedron.GetPointIds().SetId(7, seq + node_nj * node_ni + node_ni);

                    coordXYZUG.InsertNextCell(hexahedron.GetCellType(), hexahedron.GetPointIds());

                }



        vtkDataSetMapper aHexahedronMapper = vtkDataSetMapper.New();

        aHexahedronMapper.SetInputData(coordXYZUG);

        aHexahedronMapper.SetScalarRange(0, node_nk);



        vtkActor aHexahedronActor = vtkActor.New();

        aHexahedronActor.SetMapper(aHexahedronMapper);

        aHexahedronActor.GetProperty().SetDiffuseColor(1, 1, 0);



        vtkRenderer render = vtkRenderer.New();

        vtkRenderWindow renWin = myRenderWindowControl.RenderWindow;



        renWin.AddRenderer(render);

        render.SetBackground(0, 0, 1);

        render.AddActor(aHexahedronActor);

    }

https://afengwind.files.wordpress.com/2019/02/eclipsegrid.jpg?w=584

https://afengwind.files.wordpress.com/2019/02/eclipsegrid2.jpg?w=584

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yanfeng1022/article/details/89553881