diff --git a/src/Mod/Draft/CMakeLists.txt b/src/Mod/Draft/CMakeLists.txt index 69c2d5263d..df1df9f181 100644 --- a/src/Mod/Draft/CMakeLists.txt +++ b/src/Mod/Draft/CMakeLists.txt @@ -44,6 +44,7 @@ SET (Draft_geoutils draftgeoutils/cuboids.py draftgeoutils/circles.py draftgeoutils/circles_apollonius.py + draftgeoutils/circle_inversion.py ) SET(Draft_tests diff --git a/src/Mod/Draft/DraftGeomUtils.py b/src/Mod/Draft/DraftGeomUtils.py index 9e986d3459..dbe0b908fe 100644 --- a/src/Mod/Draft/DraftGeomUtils.py +++ b/src/Mod/Draft/DraftGeomUtils.py @@ -409,85 +409,12 @@ from draftgeoutils.circles import findRadicalAxis from draftgeoutils.circles import findRadicalCenter -def pointInversion(circle, point): - """Circle inversion of a point. - - pointInversion(Circle, Vector) - - Will calculate the inversed point an return it. - If the given point is equal to the center of the circle "None" will be returned. - - See also: - http://en.wikipedia.org/wiki/Inversive_geometry - """ - if (geomType(circle) == "Circle") and isinstance(point, FreeCAD.Vector): - cen = circle.Curve.Center - rad = circle.Curve.Radius - - if DraftVecUtils.equals(cen, point): - return None - - # Inverse the distance of the point - # dist(cen -> P) = r^2 / dist(cen -> invP) - - dist = DraftVecUtils.dist(point, cen) - invDist = rad**2 / d - - invPoint = Vector(0, 0, point.z) - invPoint.x = cen.x + (point.x - cen.x) * invDist / dist; - invPoint.y = cen.y + (point.y - cen.y) * invDist / dist; - - return invPoint - - else: - FreeCAD.Console.PrintMessage("debug: pointInversion bad parameters!\n") - return None +from draftgeoutils.circle_inversion import pointInversion -def polarInversion(circle, edge): - """Return the inversion pole of a line. - - polarInversion(circle, edge): - - edge ... The polar. - i.e. The nearest point on the line is inversed. - - http://mathworld.wolfram.com/InversionPole.html - """ - - if (geomType(circle) == "Circle") and (geomType(edge) == "Line"): - nearest = circle.Curve.Center.add(findDistance(circle.Curve.Center, edge, False)) - if nearest: - inversionPole = pointInversion(circle, nearest) - if inversionPole: - return inversionPole - else: - FreeCAD.Console.PrintMessage("debug: circleInversionPole bad parameters!\n") - return None +from draftgeoutils.circle_inversion import polarInversion -def circleInversion(circle, circle2): - """ - pointInversion(Circle, Circle) - - Circle inversion of a circle. - """ - if (geomType(circle) == "Circle") and (geomType(circle2) == "Circle"): - cen1 = circle.Curve.Center - rad1 = circle.Curve.Radius - - if DraftVecUtils.equals(cen1, point): - return None - - invCen2 = Inversion(circle, circle2.Curve.Center) - - pointOnCircle2 = Vector.add(circle2.Curve.Center, Vector(circle2.Curve.Radius, 0, 0)) - invPointOnCircle2 = Inversion(circle, pointOnCircle2) - - return Part.Circle(invCen2, norm, DraftVecUtils.dist(invCen2, invPointOnCircle2)) - - else: - FreeCAD.Console.PrintMessage("debug: circleInversion bad parameters!\n") - return None +from draftgeoutils.circle_inversion import circleInversion ## @} diff --git a/src/Mod/Draft/draftgeoutils/circle_inversion.py b/src/Mod/Draft/draftgeoutils/circle_inversion.py new file mode 100644 index 0000000000..c46dfdaf2b --- /dev/null +++ b/src/Mod/Draft/draftgeoutils/circle_inversion.py @@ -0,0 +1,119 @@ +# *************************************************************************** +# * Copyright (c) 2009, 2010 Yorik van Havre * +# * Copyright (c) 2009, 2010 Ken Cline * +# * * +# * This file is part of the FreeCAD CAx development system. * +# * * +# * This program is free software; you can redistribute it and/or modify * +# * it under the terms of the GNU Lesser General Public License (LGPL) * +# * as published by the Free Software Foundation; either version 2 of * +# * the License, or (at your option) any later version. * +# * for detail see the LICENCE text file. * +# * * +# * FreeCAD is distributed in the hope that it will be useful, * +# * but WITHOUT ANY WARRANTY; without even the implied warranty of * +# * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * +# * GNU Library General Public License for more details. * +# * * +# * You should have received a copy of the GNU Library General Public * +# * License along with FreeCAD; if not, write to the Free Software * +# * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 * +# * USA * +# * * +# *************************************************************************** +"""Provides various functions for inversive operations using a circle. + +http://en.wikipedia.org/wiki/Inversive_geometry +""" +## @package circle_inversion +# \ingroup DRAFTGEOUTILS +# \brief Provides various functions for inversive geometry operations. + +import lazy_loader.lazy_loader as lz + +import FreeCAD as App +import DraftVecUtils + +from draftgeoutils.general import NORM, geomType +from draftgeoutils.geometry import findDistance + +# Delay import of module until first use because it is heavy +Part = lz.LazyLoader("Part", globals(), "Part") + + +def pointInversion(circle, point): + """Return the circle inversion of a point. + + It will return `None` if the given point is equal to the center + of the circle. + """ + if geomType(circle) != "Circle" or isinstance(point, App.Vector): + print("debug: pointInversion bad parameters!") + return None + + cen = circle.Curve.Center + rad = circle.Curve.Radius + + if DraftVecUtils.equals(cen, point): + return None + + # Inverse the distance of the point + # dist(cen -> P) = r^2 / dist(cen -> invP) + + dist = DraftVecUtils.dist(point, cen) + invDist = rad**2 / dist + + invPoint = App.Vector(0, 0, point.z) + invPoint.x = cen.x + (point.x - cen.x) * invDist / dist + invPoint.y = cen.y + (point.y - cen.y) * invDist / dist + + return invPoint + + +def polarInversion(circle, edge): + """Return the inversion pole of a line. The edge is the polar. + + The nearest point on the line is inversed. + + http://mathworld.wolfram.com/InversionPole.html + """ + if geomType(circle) != "Circle" or geomType(edge) != "Line": + print("debug: circleInversionPole bad parameters! Must be a circle.") + return None + + nearest = circle.Curve.Center.add(findDistance(circle.Curve.Center, + edge, + False)) + if nearest: + inversionPole = pointInversion(circle, nearest) + if inversionPole: + return inversionPole + + return None + + +def circleInversion(circle, circle2): + """Circle inversion of a circle, inverting the center point. + + Returns the new circle created from the inverted center of circle2. + """ + if geomType(circle) != "Circle" or geomType(circle2) != "Circle": + print("debug: circleInversion bad parameters! Must be circles.") + return None + + cen1 = circle.Curve.Center + + cen2 = circle2.Curve.Center + rad2 = circle2.Curve.Radius + + if DraftVecUtils.equals(cen1, cen2): + return None + + invCen2 = pointInversion(circle, cen2) + + pointOnCircle2 = App.Vector.add(cen2, App.Vector(rad2, 0, 0)) + invPointOnCircle2 = pointInversion(circle, pointOnCircle2) + + return Part.Circle(invCen2, + NORM, + DraftVecUtils.dist(invCen2, invPointOnCircle2))