itk中的花式数据切割（三）_itk extractfilter 斜切平面-优快云博客

本文链接：https://blog.youkuaiyun.com/rabbitbride/article/details/72626566

本文介绍了一种基于已知点和方向向量的三维图像斜切方法，并详细解释了该方法的工作原理及其优缺点。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

前两篇基本上都是中规中矩的沿着XYZ切割，本文呢，我们开始斜着切。

斜着切需要有两个输入：已知点，方向向量（垂直于切面）。

void complateplane(double* planeOrigin,double* planeNormal)
{
	std::auto_ptr<itkPlane> plane ( new itkPlane);
	plane->SetOrigin(planeOrigin);
	plane->SetNormal(planeNormal);

	ImageType::RegionType inputRegion = this_data->GetLargestPossibleRegion();
	ImageType::SizeType size = inputRegion.GetSize();

	for (int x=0;x<size[0];x++)
	{
		for (int y=0;y<size[1];y++)
		{
			for(int z=0;z<size[2];z++)
			{
				ImageType::IndexType  point_temp;
				point_temp[0] = x;
				point_temp[1] = y;
				point_temp[2] = z;

				double pworld[3];
				pworld[0]=point_temp[0];
				pworld[1]=point_temp[1];
				pworld[2]=point_temp[2];

				double val=plane->EvaluateFunction(pworld);

				if (val < 1 && val > -1)
				{
					this_data->SetPixel(point_temp,1);
				}
			}
		}
	}
}

斜着切的思路完全依赖一个数学公式：

double something = normal[0]*(seed[0]-origin[0]) + normal[1]*(seed[1]-origin[1]) + normal[2]*(seed[2]-origin[2])

其中：

normal：法向量

seed：已知点

origin：图像原点

通过判断something的值来选择对应的数据。

1.something = 0

此时得到的结果是一个面，即通过seed点并且垂直于normal方向向量的空间平面

2.something > 0

1中得到的空间平面会将图像数据分成两部分，依赖右手定则，左侧的部分就是大于0

3.something < 0

相对于2中描述的空间，另外一边的空间小于0

（汗，好拗口，凑合看，简直不会说中文）

优点：这种切割方法不只可以用来提取某一部分，而且可以在原始数据中划分区域。

缺点：虽然比起前两篇的切割方法要灵活，但记得一定要把法向量算对。

参考文献：http://www.vtk.org/doc/nightly/html/classvtkPlane.html