Draft: move functions to draftgeoutils.circle_inversion
This commit is contained in:
vocx-fc
@@ -44,6 +44,7 @@ SET (Draft_geoutils
|
||||
draftgeoutils/cuboids.py
|
||||
draftgeoutils/circles.py
|
||||
draftgeoutils/circles_apollonius.py
|
||||
draftgeoutils/circle_inversion.py
|
||||
)
|
||||
|
||||
SET(Draft_tests
|
||||
|
||||
@@ -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
|
||||
|
||||
## @}
|
||||
|
||||
119
src/Mod/Draft/draftgeoutils/circle_inversion.py
Normal file
119
src/Mod/Draft/draftgeoutils/circle_inversion.py
Normal file
@@ -0,0 +1,119 @@
|
||||
# ***************************************************************************
|
||||
# * Copyright (c) 2009, 2010 Yorik van Havre <yorik@uncreated.net> *
|
||||
# * Copyright (c) 2009, 2010 Ken Cline <cline@frii.com> *
|
||||
# * *
|
||||
# * 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))
|
||||
Reference in New Issue
Block a user