已知空间两点组成的直线求线上某点的Z值,为什么会有这种看起来比较奇怪的求值需求呢?因为真正三维空间的几何计算是比较麻烦的,很多时候需要投影到二维,再反推到三维空间上去。
复习下空间直线方程:已知空间上一点(M0(x0,y0,z0))和方向向量(S(m,n,p)),则直线方程的点向式为:
[frac{X-x0}{m}=frac{Y-y0}{n}=frac{Z-z0}{p}
]
根据该公式可以解决该计算几何问题,具体实现代码如下:
#include
using namespace std;
//三维double矢量
struct Vec3d
{
double x, y, z;
Vec3d()
{
x = 0.0;
y = 0.0;
z = 0.0;
}
Vec3d(double dx, double dy, double dz)
{
x = dx;
y = dy;
z = dz;
}
void Set(double dx, double dy, double dz)
{
x = dx;
y = dy;
z = dz;
}
};
bool CalLinePointZ(const Vec3d & v1, const Vec3d & v2, Vec3d & vp)
{