opencascade计算点到特征线的距离应该注意的问题

经常会遇到要求点到曲线的距离的问题，在运用opencascade计算点到曲线和曲面的距离常用的两个类是：

Extrema_ExtPC， Extrema_ExtPS

一般做法是先构造曲线特征：（构造了Wire拓扑边）

TopoDS_Wire w = BRepBuilderAPI_MakePolygon(
gp_Pnt(0,0,0),
gp_Pnt(10,0,0),
gp_Pnt(10,10,0),
gp_Pnt(0,10,0),Standard_True);


然后将其转化成曲线：
//create an adaptor curve
BRepAdaptor_CompCurve pseudoCurve(w);
//get the parameters
Standard_Real f = pseudoCurve.FirstParameter();
Standard_Real l = pseudoCurve.LastParameter();
//f = 0.0 ; l = 4.0

创建点，检验计算是否正确：
//create 2 points to project on the wire
gp_Pnt P1(5,0,0);
gp_Pnt P2(0,5,0);


gp_Pnt newP1,newP2;
Standard_Real newVal1,newVal2;
Extrema_ExtPC ext1(P1,pseudoCurve,f,l);
if ( ext1.IsDone() && ext1.NbExt() > 0 )
{
for ( Standard_Integer i = 1 ; i < ext1.NbExt() ; i++ )
{
if ( ext1.IsMin(i) )
{
newP1 = ext1.Point(i).Value();
cout << "points: " << newP1.X() << " " << newP1.Y() << " " << newP1.Z() << endl;
cout << "distance: " << sqrt((P1.X() - newP1.X() ) * (P1.X() - newP1.X() ) + (P1.Y() - newP1.Y() ) * (P1.Y() - newP1.Y() ) + (P1.Z() - newP1.Z() ) * (P1.Z() - newP1.Z() ) ) << endl;
//break;
}
}
}

输出结果：

最后有个非常关键的问题：

这里计算出了一系列的最短距离，开始看到一个函数

IsMin就让容易让人想到判断是不是就是最短距离的，

其实不是，这个类估计是分段计算最短距离的，

结果保存的是每一段中得到的最短距离，因此，每一个都是最短距离

非常值得注意的是：

1.计算到底哪个是最短距离 需要后续处理，这个非常关键，容易搞错。

2.给定任意的一个点，经常会出现，找不到最短距离的点，因此这种情况需要特殊处理。