+
+std::shared_ptr<GeomAPI_Sphere> GeomAPI_Face::getSphere() const
+{
+ GeomSpherePtr aSphere;
+
+ const TopoDS_Face& aFace = TopoDS::Face(impl<TopoDS_Shape>());
+ Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
+ if (aSurf->IsKind(STANDARD_TYPE(Geom_SphericalSurface))) {
+ gp_Sphere aSph = Handle(Geom_SphericalSurface)::DownCast(aSurf)->Sphere();
+ const gp_Pnt& aCenter = aSph.Location();
+ double aRadius = aSph.Radius();
+
+ GeomPointPtr aCenterPnt(new GeomAPI_Pnt(aCenter.X(), aCenter.Y(), aCenter.Z()));
+ aSphere = GeomSpherePtr(new GeomAPI_Sphere(aCenterPnt, aRadius));
+ }
+ return aSphere;
+}
+
+std::shared_ptr<GeomAPI_Cylinder> GeomAPI_Face::getCylinder() const
+{
+ GeomCylinderPtr aCylinder;
+
+ const TopoDS_Face& aFace = TopoDS::Face(impl<TopoDS_Shape>());
+ Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
+ if (aSurf->IsKind(STANDARD_TYPE(Geom_CylindricalSurface))) {
+ gp_Cylinder aCyl = Handle(Geom_CylindricalSurface)::DownCast(aSurf)->Cylinder();
+ gp_Pnt aLoc = aCyl.Location();
+ const gp_Dir& aDir = aCyl.Position().Direction();
+ double aRadius = aCyl.Radius();
+
+ double aUMin, aUMax, aVMin, aVMax;
+ BRepTools::UVBounds(aFace, aUMin, aUMax, aVMin, aVMax);
+ double aHeight = aVMax - aVMin;
+
+ aLoc.ChangeCoord() += aDir.XYZ() * aVMin;
+ GeomPointPtr aLocation(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z()));
+ GeomDirPtr aDirection(new GeomAPI_Dir(aDir.X(), aDir.Y(), aDir.Z()));
+ aCylinder = GeomCylinderPtr(new GeomAPI_Cylinder(aLocation, aDirection, aRadius, aHeight));
+ }
+ return aCylinder;
+}
+
+std::shared_ptr<GeomAPI_Cone> GeomAPI_Face::getCone() const
+{
+ GeomConePtr aCone;
+
+ const TopoDS_Face& aFace = TopoDS::Face(impl<TopoDS_Shape>());
+ Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
+ if (aSurf->IsKind(STANDARD_TYPE(Geom_ConicalSurface))) {
+ gp_Cone aCon = Handle(Geom_ConicalSurface)::DownCast(aSurf)->Cone();
+ gp_Pnt aLoc = aCon.Location();
+ gp_Dir aDir = aCon.Position().Direction();
+
+ double aUMin, aUMax, aVMin, aVMax;
+ BRepTools::UVBounds(aFace, aUMin, aUMax, aVMin, aVMax);
+
+ double aSemiAngle = Abs(aCon.SemiAngle());
+ double aRadius1 = Abs(aCon.RefRadius() + aVMin * Sin(aCon.SemiAngle()));
+ double aRadius2 = Abs(aCon.RefRadius() + aVMax * Sin(aCon.SemiAngle()));
+
+ aLoc.ChangeCoord() += aDir.XYZ() * aVMin * Cos(aCon.SemiAngle());
+
+ GeomPointPtr aLocation(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z()));
+ GeomDirPtr aDirection(new GeomAPI_Dir(aDir.X(), aDir.Y(), aDir.Z()));
+ aCone = GeomConePtr(new GeomAPI_Cone(aLocation, aDirection, aSemiAngle, aRadius1, aRadius2));
+ }
+ return aCone;
+}
+
+std::shared_ptr<GeomAPI_Torus> GeomAPI_Face::getTorus() const
+{
+ GeomTorusPtr aTorus;
+
+ const TopoDS_Face& aFace = TopoDS::Face(impl<TopoDS_Shape>());
+ Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
+ if (aSurf->IsKind(STANDARD_TYPE(Geom_ToroidalSurface))) {
+ gp_Torus aTor = Handle(Geom_ToroidalSurface)::DownCast(aSurf)->Torus();
+ const gp_Pnt& aLoc = aTor.Location();
+ const gp_Dir& aDir = aTor.Position().Direction();
+ double aMajorRadius = aTor.MajorRadius();
+ double aMinorRadius = aTor.MinorRadius();
+
+ GeomPointPtr aCenter(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z()));
+ GeomDirPtr aDirection(new GeomAPI_Dir(aDir.X(), aDir.Y(), aDir.Z()));
+ aTorus = GeomTorusPtr(new GeomAPI_Torus(aCenter, aDirection, aMajorRadius, aMinorRadius));
+ }
+ return aTorus;
+}
+
+GeomPointPtr GeomAPI_Face::middlePoint() const
+{
+ GeomPointPtr anInnerPoint;
+
+ const TopoDS_Face& aFace = impl<TopoDS_Face>();
+ if (aFace.IsNull())
+ return anInnerPoint;
+
+ BRepGProp_Face aProp(aFace);
+ double aUMin, aUMax, aVMin, aVMax;
+ aProp.Bounds(aUMin, aUMax, aVMin, aVMax);
+
+ Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
+ if (aSurf.IsNull())
+ return anInnerPoint;
+
+ gp_Pnt aPnt = aSurf->Value((aUMin + aUMax) * 0.5, (aVMin + aVMax) * 0.5);
+ anInnerPoint = GeomPointPtr(new GeomAPI_Pnt(aPnt.X(), aPnt.Y(), aPnt.Z()));
+ return anInnerPoint;
+}