工具之: 赤平投影 软件计算

背景:

在岩土工程等行业中,我们时常需要通过赤平投影来判断岩体的走向,已经为层移等因素做进一步参考工作;

尽管市面上有各种程序已经实现了这方面的功能.但是基于一些技术迭代等原因, 一方面程序失去维护,无法运行. 一方面,旧的程序已经满足不了现在日新月异的需求.所以,就打算写个小软件,来应对更复杂的场景;

要求的小软件:

1.可以导入各种模型;

2.导入激光点云数据;

3.拟合求解最佳平面;

4,交互式对赤平投影的各种操作,

5.实时修改倾向倾角,实时计算赤平投影; 

6.可以依据真实微缩地球,进行查验;

初步的j结果:

 

 过程:

        参考网上一些资料以及ppt,弄明白了岩石基坑的一些概念,并结合着实际工程的理解.

《构造地质学》第12章极射赤平投影的原理和应用.ppt-原创力文档 (book118.com)

 部分关键代码:

		static vtkPolyData* CreateCirclePointsFromRadian(double center[3], double radius, int resolution, double r = 1.0, double g = 1.0, double b = 1.0, double a = 1.0)
		{
			double pi = vtkMath::Pi();
			double L = 2 * pi * radius; //周长.
			if (resolution < 0) resolution = 1;
			double radianStep = L / resolution; //每次前进的弧度长.
			double angleStep = 360.0 / resolution;//每次前进的角度.
			auto Angle2Degree = [&](double angle) {

				return  pi / 180 * angle;
			};

			//Z角度都一样.
			vtkNew<vtkPoints> mPoints;
			vtkNew<vtkCellArray> pointIndexs;
			for (size_t i = 0; i < resolution; i++)
			{
				double value = Angle2Degree(i * angleStep);
				double x = center[0] + radius * std::cos(value);
				double y = center[1] + radius * std::sin(value);
				double z = center[2];
				vtkIdType id[1]{};
				id[0] = mPoints->InsertNextPoint(x, y, z);
				pointIndexs->InsertNextCell(1, id);
			}

			vtkSmartPointer<vtkPolyData> mPolydata = vtkPolyData::New();
			mPolydata->SetPoints(mPoints);
			mPolydata->SetVerts(pointIndexs);
			return mPolydata;
		}