billy_gisboy

VTK（Activiz）场景应用：将离散点构建三角网、等值线、颜色渲染、等值线值标注

1 概括

VTK将离散点构建三角网并进行颜色渲染，生成等值线，对等值线的值进行标注

现阶段了解下，对已有离散点进行三角剖分构建三角网，并进行具有实际含义（降雨量、含水量）的属性颜色渲染、等值线构建、等值线值的标注有如下思路。特别感谢前辈“硫酸亚铜”的经验分享，前辈的博客基于VTK的C++版本进行，个人开发采用的是VTK的C#版本（Activiz），现将实现总结如下。

前辈“硫酸亚铜”的Demo生成效果：

**
我的Demo实现效果C#版本（Activiz）：解释均在代码注释中
（1）以结点为单位构建标量思路

（2）以Cell为单位构建标量思路

（两个思路构建的颜色过渡情况不一致是因为标量设置的值不一样，前提数据都是一样的，可见下图的网格显示）

2 数据前提

**必备数据：**离散点数据（X1，Y1，Z1）、（X2，Y2，Z2）、（X3，Y3，Z3）…
**中间数据：**离散点三角剖分拓扑关系（即哪几个点构成了哪一个三角形），常用数据组成为（三角形序号，点1序号，点2序号，点3序号）

离散点数据就是我们最最基础的源数据了，这个是不可或缺的，也就是作为已知数据。对于中间数据（三角剖分拓扑关系数据）如果不是已知数据，我们有两种方案进行生成：
（1）基于Triangle工具进行三角剖分生成（我使用的Triangle工具）
（2）基于VTK自身的Delaunay多边形过滤器（vtkDelaunay2D/vtkDelaunay3D）对点集（VtkPoints）进行过滤生成
现阶段中间数据作为已知数据进行处理。

3 构建思路

针对于我们有离散点和离散点三角剖分拓扑关系的数据前提去构建等值线和标注，有两种构建模式：
（1）构建点集——构建Cell单元集（拓扑关系）——构建结点标量集——构建模型
（2）构建点集——构建Cell单元集（拓扑关系）——构建Cell标量集——构建模型
第一种方案是将每一个结点看作一个标量单位，适用于对点元模型的渲染，即一个区域由一堆点元控制，每个点元（结点）有一个数据属性，例如一口井的水位值、一座烟囱的排放量…
第二种则将每一个Cell单元看作一个标量单位，适用于面元模型的渲染，即一片区域被划分为了一个个的单元且每一个单元有一个数据属性，例如每个单元的降雨量、事故发生率…

4 构建过程

4.1 两种思路的公用数据构建部分

第一步：将离散点构建成点集 vtkPoints

			/*构建网格结点*/
            vtkPoints points = new vtkPoints();
            points.InsertPoint(1, 0, 0, 0);
            points.InsertPoint(2, 1, 0, 0);
            points.InsertPoint(3, 2, 0, 0);
            points.InsertPoint(4, 0, 1, 0);
            points.InsertPoint(5, 1, 1, 0);
            points.InsertPoint(6, 2, 1, 0);
            points.InsertPoint(7, 0, 2, 0);
            points.InsertPoint(8, 1, 2, 0);
            points.InsertPoint(9, 2, 2, 0);
            points.InsertPoint(10, 3, 2, 0);
            points.InsertPoint(11, 3, 1, 0);

**第二步：构建网格单元（即拓扑关系） **vtkCellArray

            /*构建网格单元*/
            vtkIdList idList1 = new vtkIdList();
            idList1.InsertNextId(1);
            idList1.InsertNextId(2);
            idList1.InsertNextId(4);
            vtkIdList idList2 = new vtkIdList();
            idList2.InsertNextId(4);
            idList2.InsertNextId(5);
            idList2.InsertNextId(7);
            vtkIdList idList3 = new vtkIdList();
            idList3.InsertNextId(2);
            idList3.InsertNextId(5);
            idList3.InsertNextId(4);
            vtkIdList idList4 = new vtkIdList();
            idList4.InsertNextId(3);
            idList4.InsertNextId(6);
            idList4.InsertNextId(5);
            vtkIdList idList5 = new vtkIdList();
            idList5.InsertNextId(2);
            idList5.InsertNextId(3);
            idList5.InsertNextId(5);
            vtkIdList idList6 = new vtkIdList();
            idList6.InsertNextId(6);
            idList6.InsertNextId(9);
            idList6.InsertNextId(8);
            vtkIdList idList7 = new vtkIdList();
            idList7.InsertNextId(5);
            idList7.InsertNextId(6);
            idList7.InsertNextId(8);
            vtkIdList idList8 = new vtkIdList();
            idList8.InsertNextId(5);
            idList8.InsertNextId(8);
            idList8.InsertNextId(7);
            vtkIdList idList9 = new vtkIdList();
            idList9.InsertNextId(9);
            idList9.InsertNextId(6);
            idList9.InsertNextId(11);
            vtkIdList idList10 = new vtkIdList();
            idList10.InsertNextId(9);
            idList10.InsertNextId(11);
            idList10.InsertNextId(10);

            vtkCellArray cellArray = new vtkCellArray();
            cellArray.InsertNextCell(idList1);
            cellArray.InsertNextCell(idList2);
            cellArray.InsertNextCell(idList3);
            cellArray.InsertNextCell(idList4);
            cellArray.InsertNextCell(idList5);
            cellArray.InsertNextCell(idList6);
            cellArray.InsertNextCell(idList7);
            cellArray.InsertNextCell(idList8);
            cellArray.InsertNextCell(idList9);
            cellArray.InsertNextCell(idList10);

4.2 两种思路的不同构建部分

4.2.1 结点标量构建思路

**第三步：设置结点标量 **vtkDataArray（根据标量数据类型选择其子类，此处选的Float）

            /*设置结点标量*/
            vtkFloatArray pointScalar = new vtkFloatArray();
            //这个插入方式必须和单元集（vtkCellArray）的插入方式一致，单元集插入Cell单元采用InsertNextCell自动递增插入
            pointScalar.InsertNextValue(1); //1
            pointScalar.InsertNextValue(2); //2
            pointScalar.InsertNextValue(3); //3
            pointScalar.InsertNextValue(2); //4
            pointScalar.InsertNextValue(3); //5
            pointScalar.InsertNextValue(4); //6
            pointScalar.InsertNextValue(3); //7
            pointScalar.InsertNextValue(4); //8
            pointScalar.InsertNextValue(5); //9
            pointScalar.InsertNextValue(4); //10
            pointScalar.InsertNextValue(5); //11

            double rangeMin = pointScalar.GetRange()[0];
            double rangeMax = pointScalar.GetRange()[1];

第四步：构建表面几何模型

            /*构建表面几何模型*/
            //PolyData
            vtkPolyData surfacePolydata = vtkPolyData.New();
            surfacePolydata.SetPoints(points);
            surfacePolydata.SetPolys(cellArray); //设置数据类型为多边形，并设置多边形单元数据
            surfacePolydata.GetPointData().SetScalars(pointScalar); //设置点集标量
            //Actor：表面1 构建无规则网格单元的表面模型，效果是网格单元中的一个个面片，有颜色过渡
            vtkDataSetMapper surfaceMapper = new vtkDataSetMapper();
            surfaceMapper.SetInput(surfacePolydata);
            surfaceMapper.ScalarVisibilityOn(); //开启标量颜色渲染
            surfaceMapper.SetScalarRange(rangeMin, rangeMax);
            //surfaceMapper.SetLookupTable(vtkLookupTable lookupTable); //此处可以根据标量范围设置自己需要的颜色映射表，其它的制图器Mapper也一样
            vtkActor surfaceActor = new vtkActor();
            surfaceActor.SetMapper(surfaceMapper);

第五步：构建等值线

            /*构建等值线*/
            //Fileter：等值线/面过滤器
            vtkContourFilter contourFilter = new vtkContourFilter();
            contourFilter.SetInput(surfacePolydata);
            contourFilter.GenerateValues(8, rangeMin, rangeMax);
            contourFilter.Update();
            //数据处理类：用于将离散的三角面片拼接为连续的等值面
            vtkStripper stripper = new vtkStripper();
            stripper.SetInput(contourFilter.GetOutput());
            stripper.Update();
            //Mapper
            vtkDataSetMapper contourMapper = new vtkDataSetMapper();
            contourMapper.SetInput(stripper.GetOutput());
            //Actor
            vtkActor isolineActor = new vtkActor();
            isolineActor.SetMapper(contourMapper);

第六步：构建等值线标注

            /*等值线值标注*/
            //构建等值线标注几何数据
            //long linesCount = stripper.GetOutput().GetNumberOfLines(); //等值线的数目
            vtkPoints isolinePoints = stripper.GetOutput().GetPoints(); //等值线的点集
            vtkCellArray isolineCells = stripper.GetOutput().GetLines(); //等值线的线单元数组
            vtkDataArray isolineScalars = stripper.GetOutput().GetPointData().GetScalars(); //等值线的点标量数组
            //等值线标注Label数据
            vtkPolyData labelPolyData = new vtkPolyData(); //标注点几何数据
            vtkPoints labelPoints = new vtkPoints(); //标注点集
            vtkDoubleArray labelScalars = new vtkDoubleArray(); //标注点标量数组
            labelScalars.SetNumberOfComponents(1);
            //labelScalars.SetName("Isovalues");
            //构建等值线标注Label数据，此处将每一个等值线进行遍历，获取线上的一点作为等值线的标注点
            vtkIdList isoCell = new vtkIdList();
            while (isolineCells.GetNextCell(isoCell) > 0)
            {
                long idCount = isoCell.GetNumberOfIds(); //获取这条等值线单元的点个数
                //long samplePointIndex = (long)vtkMath.Random(0, idCount); //随机取该等值线上一个点
                long labelPointId = isoCell.GetId(idCount / 2); //获取该点的id作为标记点。此方法是根据点数量顺序排序序号Index获取，获取到这个cell中的指定序号Index的点的vtkId（个人认为）

                double[] midPoint = isolinePoints.GetPoint(labelPointId); //获取该点的坐标

                labelPoints.InsertNextPoint(midPoint[0], midPoint[1], midPoint[2]); //插入到标注点集
                labelScalars.InsertNextTuple1(isolineScalars.GetTuple1(labelPointId)); //取该点的标量值插入到标注点标量数组
            }
            labelPolyData.SetPoints(labelPoints); //设置标注点集
            labelPolyData.GetPointData().SetScalars(labelScalars); //设置点数据的标量数组
            labelPolyData.Update();
            //Mapper
            vtkLabeledDataMapper labelMapper = new vtkLabeledDataMapper();
            //labelMapper.SetFieldDataName("Isovalues");
            labelMapper.SetInput(labelPolyData);
            labelMapper.SetLabelModeToLabelScalars();
            labelMapper.SetLabelFormat("%6.2f");
            labelMapper.GetLabelTextProperty().SetColor(0, 0, 1);
            //Actor
            vtkActor2D isoLabelActor = new vtkActor2D();
            isoLabelActor.SetMapper(labelMapper);

第七步：添加角色到场景

			SceneControl.AddActor(surfaceActor);
            SceneControl.AddActor(isolineActor);
            SceneControl.AddActor2D(isoLabelActor);
            SceneControl.RefreshScene(true);

这个SceneControl就是vtkSceneControl，我是基于C#版本的VTK进行开发的，为了便于开发工作，将其控件进行了封装。

4.2.2 单元标量构建思路

**第三步：设置单元标量 **vtkDataArray（根据标量数据类型选择其子类，此处选的Float）

            /*设置单元标量*/
            vtkFloatArray cellScalar = new vtkFloatArray();
            //这个插入方式必须和单元集（vtkCellArray）的插入方式一致，单元集插入cell单元采用InsertNextCell自动递增插入
            cellScalar.InsertNextValue(1); //1
            cellScalar.InsertNextValue(2); //2
            cellScalar.InsertNextValue(2); //3
            cellScalar.InsertNextValue(3); //4
            cellScalar.InsertNextValue(2); //5
            cellScalar.InsertNextValue(4); //6
            cellScalar.InsertNextValue(3); //7
            cellScalar.InsertNextValue(3); //8
            cellScalar.InsertNextValue(3); //9
            cellScalar.InsertNextValue(4); //10

            double rangeMin = cellScalar.GetRange()[0];
            double rangeMax = cellScalar.GetRange()[1];