src/GEOMUtils/GEOMUtils.hxx

   1 // Copyright (C) 2007-2015  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   6 // This library is free software; you can redistribute it and/or
   7 // modify it under the terms of the GNU Lesser General Public
   8 // License as published by the Free Software Foundation; either
   9 // version 2.1 of the License, or (at your option) any later version.
  10 //
  11 // This library is distributed in the hope that it will be useful,
  12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 // Lesser General Public License for more details.
  15 //
  16 // You should have received a copy of the GNU Lesser General Public
  17 // License along with this library; if not, write to the Free Software
  18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22
  23 #ifndef _GEOMUtils_HXX_
  24 #define _GEOMUtils_HXX_
  25
  26 #include <Standard_Macro.hxx>
  27 #include <TopoDS_Shape.hxx>
  28 #include <TopoDS_Vertex.hxx>
  29
  30 #include <TopTools_ListOfShape.hxx>
  31
  32 #include <TopAbs.hxx>
  33
  34 #include <gp_Ax3.hxx>
  35 #include <gp_Vec.hxx>
  36
  37 #include <V3d_View.hxx>
  38
  39 #include <NCollection_DataMap.hxx>
  40
  41 #include <functional>
  42
  43 #include <map>
  44 #include <vector>
  45 #include <string>
  46 #include <utility>
  47
  48 class Bnd_Box;
  49
  50 inline Standard_Boolean IsEqual (const TopoDS_Shape& S1, const TopoDS_Shape& S2)
  51 {
  52   return S1.IsSame(S2);
  53 }
  54
  55 namespace GEOMUtils
  56 {
  57
  58   /**
  59    * This enumeration represents comparison conditions.
  60    */
  61   enum ComparisonCondition {
  62     CC_GT, ///< Greater then
  63     CC_GE, ///< Greater then or equal to
  64     CC_LT, ///< Less then
  65     CC_LE  ///< Less then or equal to
  66   };
  67
  68   typedef std::vector<std::string> NodeLinks;
  69   typedef std::map<std::string, NodeLinks> LevelInfo;
  70   typedef std::vector<LevelInfo> LevelsList;
  71   typedef std::map<std::string,std::pair<LevelsList,LevelsList> > TreeModel;
  72
  73   /*!
  74    * \brief Compute numerical functor for the shape.
  75    *
  76    * Resulting value can be used to sort out shapes according to some parameter.
  77    *
  78    * Returns a pair of two values (dist, functor) where
  79    * - \a dist is a some value that is computed according to the center of mass of given shape;
  80    * - \a functor is a numerical functor value
  81    *
  82    * The numerical functor is computed according to the shape's topological properties as follows:
  83    * - orientation for vertices
  84    * - length for edges and wires
  85    * - area for faces and shells
  86    * - volume for solids, compounds, compsolids
  87    *
  88    * If \a isOldSorting parameter is set to \c true, for all cases linear properties of the shape
  89    * are used (to support backward compatibility in some methods). By default, this parameter is
  90    * set to \c false.
  91    */
  92   Standard_EXPORT std::pair<double, double> ShapeToDouble (const TopoDS_Shape& theShape,
  93                                                            bool isOldSorting = false);
  94
  95   /*!
  96    * \brief Get Local Coordinate System, corresponding to the given shape.
  97    *
  98    * Origin of the LCS is situated at the shape's center of mass.
  99    * Axes of the LCS are obtained from shape's location or,
 100    * if the shape is a planar face, from position of its plane.
 101    */
 102   Standard_EXPORT gp_Ax3 GetPosition (const TopoDS_Shape& theShape);
 103
 104   /*!
 105    * \brief Get vector, defined by the given edge.
 106    * \param theShape The edge.
 107    * \param doConsiderOrientation If True, take into account the edge orientation.
 108    * \note It is recommended to use doConsiderOrientation=Standard_False, because
 109    *       the same edge can have different orientation depending on the way it was
 110    *       extracted from a shape.
 111    */
 112   Standard_EXPORT gp_Vec GetVector (const TopoDS_Shape& theShape,
 113                                     Standard_Boolean doConsiderOrientation);
 114
 115   /*!
 116    * \brief Sort shapes in the list by their coordinates.
 117    * \param SL The list of shapes to sort.
 118    */
 119   struct CompareShapes : public std::binary_function<TopoDS_Shape, TopoDS_Shape, bool>
 120   {
 121     CompareShapes (bool isOldSorting)
 122       : myIsOldSorting(isOldSorting) {}
 123
 124     bool operator() (const TopoDS_Shape& lhs, const TopoDS_Shape& rhs);
 125
 126     typedef NCollection_DataMap<TopoDS_Shape, std::pair<double, double> > GEOMUtils_DataMapOfShapeDouble;
 127     GEOMUtils_DataMapOfShapeDouble myMap;
 128     bool myIsOldSorting;
 129   };
 130
 131   /*!
 132    * \brief Sort shapes by their centers of mass, using formula X*999 + Y*99 + Z*0.9
 133    */
 134   Standard_EXPORT void SortShapes (TopTools_ListOfShape& SL,
 135                                    const Standard_Boolean isOldSorting = Standard_True);
 136
 137   /*!
 138    * \brief Convert TopoDS_COMPSOLID to TopoDS_COMPOUND.
 139    *
 140    * If the argument shape is not of type TopoDS_COMPSOLID, this method returns it as is.
 141    *
 142    * \param theCompsolid The compsolid to be converted.
 143    * \retval TopoDS_Shape Returns the resulting compound.
 144    */
 145   Standard_EXPORT TopoDS_Shape CompsolidToCompound (const TopoDS_Shape& theCompsolid);
 146
 147   /*!
 148    * \brief Recursively extract all shapes from compounds and compsolids of the given shape into theList.
 149    *
 150    * If theShape is not compound or compsolid, theList will contain only theShape itself.
 151    *
 152    * \param theShape The shape to be exploded.
 153    * \param theList Output parameter.
 154    */
 155   Standard_EXPORT void AddSimpleShapes (const TopoDS_Shape& theShape,
 156                                         TopTools_ListOfShape& theList);
 157
 158   /*!
 159    * \brief Build a triangulation on \a theShape if it is absent.
 160    * \param theShape The shape to check/build triangulation on.
 161    * \retval bool Returns false if the shape has no faces, i.e. impossible to build triangulation.
 162    */
 163   Standard_EXPORT bool CheckTriangulation (const TopoDS_Shape& theShape);
 164
 165   /*!
 166    * \brief Return type of shape for explode. In case of compound it will be a type of its first sub shape.
 167    * \param theShape The shape to get type of.
 168    * \retval TopAbs_ShapeEnum Return type of shape for explode.
 169    */
 170   Standard_EXPORT TopAbs_ShapeEnum GetTypeOfSimplePart (const TopoDS_Shape& theShape);
 171
 172   /*!
 173    * \brief Find an edge of theShape, closest to thePoint.
 174    *
 175    * \param theShape The shape to explore.
 176    * \param thePoint The point near the required edge.
 177    * \retval TopoDS_Shape Returns the found edge or an empty shape if multiple edges found.
 178    */
 179   Standard_EXPORT TopoDS_Shape GetEdgeNearPoint (const TopoDS_Shape&  theShape,
 180                                                  const TopoDS_Vertex& thePoint);
 181
 182   /*!
 183    * \brief Compute precise bounding box of the shape based on the rough bounding box.
 184    *
 185    * \param theShape the shape.
 186    * \param theBox rough bounding box on input; precise bounding box on output.
 187    * \retval Standard_True in case of success; Standard_False otherwise.
 188    */
 189   Standard_EXPORT Standard_Boolean PreciseBoundingBox(const TopoDS_Shape &theShape, Bnd_Box &theBox);
 190
 191   /*!
 192    * \brief Computes minumal distance between two shapes for singular cases
 193    *        (workaround for bugs 19899, 19908 and 19910 from Mantis).
 194    *
 195    * \param aSh1 the first shape
 196    * \param aSh2 the second shape
 197    * \param Ptmp1 the output result point on the first shape
 198    * \param Ptmp2 the output result point on the second shape
 199    * \retval negative value if it is not a singular case; actual distance for singular case.
 200    */
 201   Standard_EXPORT Standard_Real GetMinDistanceSingular(const TopoDS_Shape& aSh1,
 202                                                        const TopoDS_Shape& aSh2,
 203                                                        gp_Pnt& Ptmp1, gp_Pnt& Ptmp2);
 204
 205   /*!
 206    * \brief Computes minumal distance between two shapes.
 207    *
 208    * \param theShape1 the first shape
 209    * \param theShape2 the second shape
 210    * \param thePnt1 the output result point on the first shape
 211    * \param thePnt2 the output result point on the second shape
 212    * \retval negative value in case of failure; otherwise the real distance.
 213    */
 214   Standard_EXPORT Standard_Real GetMinDistance(const TopoDS_Shape& theShape1,
 215                                                const TopoDS_Shape& theShape2,
 216                                                gp_Pnt& thePnt1, gp_Pnt& thePnt2);
 217
 218   /*!
 219    * \brief Returns the point clicked in 3D view.
 220    *
 221    * \param x The X coordinate in the view.
 222    * \param y The Y coordinate in the view.
 223    * \param theView View where the given point takes place.
 224    * \retval gp_Pnt Returns the point clicked in 3D view
 225    */
 226   Standard_EXPORT gp_Pnt ConvertClickToPoint( int x, int y, Handle(V3d_View) theView );
 227
 228   /*!
 229    * \brief Convert dependency tree data to the string representation
 230    *
 231    * \param tree dependency tree data
 232    * \param dependencyStr output string
 233    */
 234   Standard_EXPORT void ConvertTreeToString( const TreeModel& tree,
 235                                             std::string& dependencyStr );
 236
 237   /*!
 238    * \brief Restore dependency tree data from the string representation
 239    *
 240    * \param dependencyStr string representation of tree data
 241    * \param tree output dependency tree data
 242    */
 243   Standard_EXPORT void ConvertStringToTree( const std::string& dependencyStr,
 244                                             TreeModel& tree );
 245
 246   /*!
 247    * \brief Check shape
 248    *
 249    * \param shape input shape object
 250    * \param checkGeometry when set to \c true, causes check of underlying geometry
 251    *        in addition to the topology
 252    * \return \c true if shape is valid or \c false otherwise
 253    */
 254   Standard_EXPORT bool CheckShape( TopoDS_Shape& shape, bool checkGeometry = false );
 255
 256   /*!
 257    * \brief Limit shape tolerance to the given value
 258    *
 259    * \param shape shape being fixed
 260    * \param type topology type which tolerance is to be limited; TopAbs_SHAPE means
 261    *             all types of topology
 262    * \param tolerance expected tolerance value (1e-7 by default)
 263    * \param checkGeometry check geometry validity of result
 264    * \return \c true if resulting shape is valid
 265    *
 266    * \note Resulting tolerance of the shape is not mandatory equal to requested value
 267    *       as it might be changed by fixshape operation in order to get valid shape where possible
 268    * \note By default, result only checked for topology validity; check of geometry can be done by
 269    *       passing \c true to \a checkGeometry parameter
 270    */
 271   Standard_EXPORT bool FixShapeTolerance( TopoDS_Shape& shape,
 272                                           TopAbs_ShapeEnum type,
 273                                           Standard_Real tolerance = Precision::Confusion(),
 274                                           bool checkGeometry = false );
 275
 276   /*!
 277    * \brief Limit shape tolerance to the given value
 278    * This is overloaded function, it behaves exactly as previous one
 279    */
 280   Standard_EXPORT bool FixShapeTolerance( TopoDS_Shape& shape,
 281                                           Standard_Real tolerance = Precision::Confusion(),
 282                                           bool checkGeometry = false );
 283
 284   /*!
 285    * \brief Limit shape tolerance to the given value
 286    * This is overloaded function, it behaves exactly as previous one
 287    */
 288   Standard_EXPORT bool FixShapeTolerance( TopoDS_Shape& shape,
 289                                           bool checkGeometry );
 290
 291   /*!
 292    * \brief Fix curves of the given shape
 293    *
 294    * The function checks each curve of the input shape in the following way:
 295    * - compute deviation of the curve from the underlying surface in a set of points
 296    *   computed with the certain discretization step value
 297    * - find maximum tolerance between computed deviation values
 298    * - limit tolerance of the curve with the computed maximum value
 299    *
 300    * \param shape shape being fixed
 301    * \return \c true if resulting shape is valid
 302    */
 303   Standard_EXPORT bool FixShapeCurves( TopoDS_Shape& shape );
 304
 305   /*!
 306    * \brief Write shape to the BREP file
 307    *
 308    * \param source shape
 309    * \return \c true if file was written or \c false otherwise
 310    */
 311   Standard_EXPORT bool Write( const TopoDS_Shape& shape,
 312                               const char* fileName );
 313
 314   /*!
 315    * \brief Extract single SOLID from COMPSOLID or COMPOUND.
 316    *
 317    * If the argument shape is a COMPUND or COMPSOLID and there's
 318    * only single simple-shape type inside, this sub-shape is returned as a result;
 319    * otherwise, the shape is not changed.
 320    *
 321    * \param shape compound or compsolid being processed.
 322    * \retval TopoDS_Shape resulting shape
 323    */
 324   Standard_EXPORT TopoDS_Shape ReduceCompound( const TopoDS_Shape& shape );
 325
 326   /*!
 327    * \brief Generate triangulation for the shape.
 328    *
 329    * \param shape shape being meshed
 330    * \param deflection deflection coefficient to be used
 331    * \param forced if \c true, causes generation of mesh regardless it is already present in the shape
 332    */
 333   Standard_EXPORT void MeshShape( const TopoDS_Shape shape,
 334                                   double deflection, bool forced = true );
 335
 336   /*!
 337    * \brief Get default deflection coefficient used for triangulation
 338    * \return default deflection value
 339    */
 340   Standard_EXPORT double DefaultDeflection();
 341
 342   /**
 343    * \brief Check if the shape is not a closed wire or edge.
 344    *
 345    * This function is used for pipe creation algorithm to test if
 346    * the pipe path is not closed. It returns false if theShape is a wire or
 347    * an edge with coincident end vertices. Otherwise it returns true.
 348    *
 349    * \param theShape the shape to be tested.
 350    * \return true if theShape is not a closed wire or edge.
 351    */
 352   Standard_EXPORT bool IsOpenPath(const TopoDS_Shape &theShape);
 353
 354   /**
 355    * This function compares two tolerances. The shape tolerance (the first
 356    * argument) is considered less than the reference tolerance (the second
 357    * argument) if theTolShape < theTolRef - Tolerance(theTolRef). theTolShape is
 358    * considered greater than theTolRef if theTolShape > theTolRef +
 359    * Tolerance(theTolRef). Otherwise these tolerances are equal.
 360    * Tolerance(theTolRef) = theTolRef*DEFAULT_TOLERANCE_TOLERANCE. But this value
 361    * should not be greated than DEFAULT_MAX_TOLERANCE_TOLERANCE.
 362    *
 363    * \param theTolShape the shape tolerance
 364    * \param theTolRef the reference tolerance
 365    * \return -1 if theTolShape is less than theTolRef; 1 if theTolShape is greater
 366    * than theTolRef; 0 if they are equal
 367    */
 368   Standard_EXPORT int CompareToleranceValues(const double theTolShape,
 369                                              const double theTolRef);
 370
 371   /**
 372    * Check if the comarison of tolerances fit the condition. The comparison of
 373    * tolerances is performed using the function CompareToleranceValues.
 374    *
 375    * \param theCondition the condition
 376    * \param theTolShape the shape tolerance
 377    * \param theTolRef the reference tolerance
 378    * \return true if the shape tolerance fits the condition; false otherwise.
 379    */
 380   Standard_EXPORT bool IsFitCondition(const ComparisonCondition theCondition,
 381                                       const double              theTolShape,
 382                                       const double              theTolRef);
 383
 384 };
 385
 386 #endif