背景:
在岩土工程等行业中,我们时常需要通过赤平投影来判断岩体的走向,已经为层移等因素做进一步参考工作;
尽管市面上有各种程序已经实现了这方面的功能.但是基于一些技术迭代等原因, 一方面程序失去维护,无法运行. 一方面,旧的程序已经满足不了现在日新月异的需求.所以,就打算写个小软件,来应对更复杂的场景;
要求的小软件:
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;
}