src/ProjLib/ProjLib_ProjectOnSurface.cxx

   1 // Created on: 1994-09-15
   2 // Created by: Bruno DUMORTIER
   3 // Copyright (c) 1994-1999 Matra Datavision
   4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   5 //
   6 // This file is part of Open CASCADE Technology software library.
   7 //
   8 // This library is free software; you can redistribute it and/or modify it under
   9 // the terms of the GNU Lesser General Public License version 2.1 as published
  10 // by the Free Software Foundation, with special exception defined in the file
  11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  12 // distribution for complete text of the license and disclaimer of any warranty.
  13 //
  14 // Alternatively, this file may be used under the terms of Open CASCADE
  15 // commercial license or contractual agreement.
  16
  17
  18 #include <Adaptor3d_HCurve.hxx>
  19 #include <Adaptor3d_HSurface.hxx>
  20 #include <AppCont_Function.hxx>
  21 #include <AppParCurves_MultiCurve.hxx>
  22 #include <Approx_FitAndDivide.hxx>
  23 #include <BSplCLib.hxx>
  24 #include <Extrema_ExtPS.hxx>
  25 #include <Extrema_POnSurf.hxx>
  26 #include <Geom_BSplineCurve.hxx>
  27 #include <PLib.hxx>
  28 #include <Precision.hxx>
  29 #include <ProjLib_ProjectOnSurface.hxx>
  30 #include <Standard_NoSuchObject.hxx>
  31 #include <TColgp_Array1OfPnt.hxx>
  32 #include <TColStd_Array1OfInteger.hxx>
  33 #include <TColStd_Array1OfReal.hxx>
  34
  35 //=======================================================================
  36 //function : OnSurface_Value
  37 //purpose  : Evaluate current point of the projected curve
  38 //=======================================================================
  39 static gp_Pnt OnSurface_Value(const Standard_Real U,
  40                               const Handle(Adaptor3d_HCurve)& myCurve,
  41                               Extrema_ExtPS * myExtPS)
  42 {
  43   // on essaie de rendre le point solution le plus proche.
  44   myExtPS->Perform(myCurve->Value(U));
  45
  46   Standard_Real Dist2Min = RealLast();
  47   Standard_Integer Index = 0;
  48
  49   for ( Standard_Integer i = 1; i <= myExtPS->NbExt(); i++) {
  50     if ( myExtPS->SquareDistance(i) < Dist2Min) {
  51       Index = i;
  52       Dist2Min = myExtPS->SquareDistance(Index);
  53     }
  54   }
  55   if ( Index == 0 ) {
  56     cout << " Extrema non trouve pour U = " << U << endl;
  57     return gp_Pnt(0.,0.,0.);
  58   }
  59   else {
  60     return (myExtPS->Point(Index)).Value();
  61   }
  62 }
  63
  64 //=======================================================================
  65 //function : OnSurface_D1
  66 //purpose  :
  67 //=======================================================================
  68
  69 static Standard_Boolean OnSurface_D1(const Standard_Real , // U,
  70                                      gp_Pnt& ,       // P,
  71                                      gp_Vec&       , // V,
  72                                      const  Handle(Adaptor3d_HCurve)& , //  myCurve,
  73                                      Extrema_ExtPS *) // myExtPS)
  74 {
  75   return Standard_False;
  76 }
  77
  78
  79 //=======================================================================
  80 //  class  : ProjLib_OnSurface
  81 //purpose  : Use to approximate the projection on a plane
  82 //=======================================================================
  83
  84 class ProjLib_OnSurface : public AppCont_Function
  85
  86 {
  87 public:
  88
  89   ProjLib_OnSurface(const Handle(Adaptor3d_HCurve)   & C,
  90                     const Handle(Adaptor3d_HSurface) & S)
  91  : myCurve(C)
  92   {
  93     myNbPnt = 1;
  94     myNbPnt2d = 0;
  95     Standard_Real U = myCurve->FirstParameter();
  96     gp_Pnt P = myCurve->Value(U);
  97     Standard_Real Tol = Precision::PConfusion();
  98     myExtPS = new Extrema_ExtPS(P,S->Surface(),Tol,Tol);
  99   }
 100
 101   ~ProjLib_OnSurface() { delete myExtPS; }
 102
 103   Standard_Real FirstParameter() const
 104     {return myCurve->FirstParameter();}
 105
 106   Standard_Real LastParameter() const
 107     {return myCurve->LastParameter();}
 108
 109   Standard_Boolean Value(const Standard_Real   theT,
 110                          NCollection_Array1<gp_Pnt2d>& /*thePnt2d*/,
 111                          NCollection_Array1<gp_Pnt>&   thePnt) const
 112   {
 113       thePnt(1) = OnSurface_Value(theT, myCurve, myExtPS);
 114       return Standard_True;
 115   }
 116
 117   Standard_Boolean D1(const Standard_Real   theT,
 118                       NCollection_Array1<gp_Vec2d>& /*theVec2d*/,
 119                       NCollection_Array1<gp_Vec>&   theVec) const
 120   {
 121     gp_Pnt aPnt;
 122     return OnSurface_D1(theT, aPnt, theVec(1), myCurve, myExtPS);
 123   }
 124
 125 private:
 126   ProjLib_OnSurface (const ProjLib_OnSurface&);
 127   ProjLib_OnSurface& operator= (const ProjLib_OnSurface&);
 128
 129 private:
 130   Handle(Adaptor3d_HCurve)       myCurve;
 131   Extrema_ExtPS*                 myExtPS;
 132 };
 133
 134
 135 //=====================================================================//
 136 //                                                                     //
 137 //  D E S C R I P T I O N   O F   T H E   C L A S S  :                 //
 138 //                                                                     //
 139 //         P r o j L i b _ A p p r o x P r o j e c t O n P l a n e     //
 140 //                                                                     //
 141 //=====================================================================//
 142
 143
 144 //=======================================================================
 145 //function : ProjLib_ProjectOnSurface
 146 //purpose  :
 147 //=======================================================================
 148
 149 ProjLib_ProjectOnSurface::ProjLib_ProjectOnSurface() :
 150 myIsDone(Standard_False)
 151 {
 152 }
 153
 154 //=======================================================================
 155 //function : ProjLib_ProjectOnSurface
 156 //purpose  :
 157 //=======================================================================
 158
 159 ProjLib_ProjectOnSurface::ProjLib_ProjectOnSurface
 160 (const Handle(Adaptor3d_HSurface)& S ) :
 161 myIsDone(Standard_False)
 162 {
 163   mySurface = S;
 164 }
 165
 166 void ProjLib_ProjectOnSurface::Load(const Handle(Adaptor3d_HCurve)& C,
 167                                     const Standard_Real  Tolerance)
 168 {
 169   myTolerance = Tolerance ;
 170   myCurve = C;
 171   myIsDone = Standard_False ;
 172   if (!mySurface.IsNull()) {
 173
 174     ProjLib_OnSurface F(myCurve, mySurface);
 175
 176     Standard_Integer Deg1, Deg2;
 177     Deg1 = 8; Deg2 = 8;
 178
 179     Approx_FitAndDivide Fit(F,Deg1,Deg2,Precision::Approximation(),
 180                             Precision::PApproximation(),Standard_True);
 181     Standard_Integer i;
 182     Standard_Integer NbCurves = Fit.NbMultiCurves();
 183     Standard_Integer MaxDeg = 0;
 184
 185     // Pour transformer la MultiCurve en BSpline, il faut que toutes
 186     // les Bezier la constituant aient le meme degre -> Calcul de MaxDeg
 187     Standard_Integer NbPoles  = 1;
 188     for (i = 1; i <= NbCurves; i++) {
 189       Standard_Integer Deg = Fit.Value(i).Degree();
 190       MaxDeg = Max ( MaxDeg, Deg);
 191     }
 192     NbPoles = MaxDeg * NbCurves + 1;               //Poles sur la BSpline
 193     TColgp_Array1OfPnt  Poles( 1, NbPoles);
 194
 195     TColgp_Array1OfPnt TempPoles( 1, MaxDeg + 1);  //pour augmentation du degre
 196
 197     TColStd_Array1OfReal Knots( 1, NbCurves + 1);  //Noeuds de la BSpline
 198
 199     Standard_Integer Compt = 1;
 200     for (i = 1; i <= Fit.NbMultiCurves(); i++) {
 201       Fit.Parameters(i, Knots(i), Knots(i+1));
 202
 203       AppParCurves_MultiCurve MC = Fit.Value( i);   //Charge la Ieme Curve
 204       TColgp_Array1OfPnt LocalPoles( 1, MC.Degree() + 1);//Recupere les poles
 205       MC.Curve(1, Poles);
 206
 207       //Augmentation eventuelle du degre
 208       Standard_Integer Inc = MaxDeg - MC.Degree();
 209       if ( Inc > 0) {
 210         BSplCLib::IncreaseDegree( Inc, LocalPoles, BSplCLib::NoWeights(),
 211                                  TempPoles, BSplCLib::NoWeights());
 212           //mise a jour des poles de la PCurve
 213           for (Standard_Integer j = 1 ; j <= MaxDeg + 1; j++) {
 214             Poles.SetValue( Compt, TempPoles( j));
 215             Compt++;
 216           }
 217         }
 218       else {
 219         //mise a jour des poles de la PCurve
 220         for (Standard_Integer j = 1 ; j <= MaxDeg + 1; j++) {
 221           Poles.SetValue( Compt, LocalPoles( j));
 222           Compt++;
 223         }
 224       }
 225
 226       Compt--;
 227     }
 228
 229     //mise a jour des fields de ProjLib_Approx
 230
 231     Standard_Integer NbKnots = NbCurves + 1;
 232
 233     TColStd_Array1OfInteger  Mults( 1, NbKnots);
 234     Mults.SetValue( 1, MaxDeg + 1);
 235     for ( i = 2; i <= NbCurves; i++) {
 236       Mults.SetValue( i, MaxDeg);
 237     }
 238     Mults.SetValue(NbKnots, MaxDeg + 1);
 239     myResult =
 240       new Geom_BSplineCurve(Poles,
 241                             Knots,
 242                             Mults,
 243                             MaxDeg,
 244                             Standard_False) ;
 245     myIsDone = Standard_True ;
 246   }
 247 }
 248
 249 //=======================================================================
 250 //function : ~ProjLib_ProjectOnSurface
 251 //purpose  :
 252 //=======================================================================
 253
 254 ProjLib_ProjectOnSurface::~ProjLib_ProjectOnSurface()
 255 {}
 256
 257 //=======================================================================
 258 //function : BSpline
 259 //purpose  :
 260 //=======================================================================
 261
 262 Handle(Geom_BSplineCurve) ProjLib_ProjectOnSurface::BSpline() const
 263 {
 264   Standard_NoSuchObject_Raise_if
 265     (!myIsDone,
 266      "ProjLib_ProjectOnSurface:BSpline");
 267   return myResult ;
 268 }
 269
 270 //=======================================================================
 271 //function : IsDone
 272 //purpose  :
 273 //=======================================================================
 274
 275 Standard_Boolean ProjLib_ProjectOnSurface::IsDone() const
 276 {
 277   return myIsDone;
 278 }