VTK实现曲线等分踩坑总结

VTK实现曲线等分踩坑总结

设计VTK算法时需要用到等距的顶点，尝试了两种能够实现的等分的VTK函数，都发现各自的坑，或者是我以为结果会是等分，实际上并非等分。

---

踩坑一：利用vtkLineSource本身提供的函数做等分操作

vtkLineSource本身会提供一个类似等分操作的函数，SetResolution函数能够给用户去设置将线条分为多少段，预想结果是多个线段组成的线条，重新划分出来的线段步长是一样的。

	vtkSmartPointer<vtkLineSource> line = vtkSmartPointer<vtkLineSource>::New();
	//line->SetPoint1(p);
	//line->SetPoint2(q);
	line->SetPoints(ps);
	line->SetUseRegularRefinement(true);	//先开启RegularRefinement
	line->SetResolution(50);				//再设置分辨率，分辨率也理解为将线条划分为多少短
	line->Update();

但是实际结果：
（1）vtkLineSource只设置Point1和Point2：表示此时line是只有两个端点连接而成的一条直线段，此时划分结果是正确等分的，每个子线段是等距的。

（2）vtkLineSource设置多个顶点组成的Points：表示line由许多个顶点连接而成，宏观看是一条连续的曲线，实际上由多个子线段组成的。此时划分结果是不等分的，如下图。下图由四个端点连接而成的line，所以有三条子线段，绿色小球是划分后的顶点作为球心，顶点步长作为球体的直径，假如是等分的，那么小球是同样大小而且是紧挨的。可以看到三个原始子线段中，线段之间的小球大小是有差异的，线段内的小球是等分的，所以整体并非等分的。

---

踩坑二：利用vtkSplineFilter提供的函数做等分操作

VTK实现对曲线进行顶点等分处理

这是我之前写的等分处理方法，经过跟上图同样的验证后，发现结果也是不等分的。

vtkSplineFilter提供SetNumberOfSubdivisions函数能够对曲线细分，得到新顶点连接而成的同一条曲线，当初目测以为是等分的，就拿来用了，后面我的算法总是得不到正确的数据，我于是对vtkSplineFilter的细分操作做了检查，才发现实际上细分后的线段并非等距的，而更像是正态分布那样，两头的细分线段相对更窄，中间部分的细分线段更宽。

不过vtkSplineFilter经过细分后能够带来一个很大的好处是，访问顶点索引时是顶点有序的，不需要自己重新排序。

---

如何真正实现对曲线的等分呢？

答案是自己实现，不难，可以用最直观的方法去实现（以下是伪代码）：

input: 有序的曲线顶点集points（要求我按顶点id顺序访问，可以实现对三维顶点拓扑有序访问。可以用vtkSplineFilter对vtkLineSource细分处理一下就能实现排序）
input: 目标等分数量num
output: 新的等分顶点集ps

step1: //计算曲线本身的总长度
		for(int i = 1; i < points->GetNumberOfPoints(); i++)
		{
			//取出i和i-1的顶点计算距离
			p <- points[i - 1]
			q <- points[i]
			dist = distance(p, q)
			length += dist	//length用于计算曲线的总长度 
		}
step2: //执行等分
		value = 0	//当前行进的数值
		cur = 0		//遍历到当前路段的总长度
		index = 0	//当前顶点索引
		segment = 0	//当前所处线段的长度
		for(int i = 1; i <= num; i++)	//注意使用小于等于
		{
			value = length * i / numOfSub
			while(value >= cur)
			{
				index++;
				a1 <- points[index - 1]
				a2 <- points[index]
				segment = distance(a1, a2)
				cur = cur + segment
			}
			//先把p点和q点的三维坐标取出来，用于后续插值出等分点
			p <- points[index - 1]
			q <- points[index]
			pq = q - p
			alphe = (value - (cur - segment)) / segment
			t <- p + pq * alpha
			ps->InsertNextPoint(t)
		}
step3:	//返回等分结果
		return ps