pygeodesy 25.5.5__py2.py3-none-any.whl → 25.5.28__py2.py3-none-any.whl

pygeodesy/auxilats/auxily.py

@@ -12,7 +12,7 @@ under the MIT/X11 License.  For more information, see the U{GeographicLib
 <https://GeographicLib.SourceForge.io>} documentation.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 # from pygeodesy import auxilats  # _MODS
 from pygeodesy.auxilats._CX_Rs import _Rkey
@@ -29,7 +29,7 @@ from pygeodesy.utily import atan1
 from math import asinh, copysign
 
 __all__ = ()
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 
 class Aux(object):

pygeodesy/azimuthal.py

@@ -40,7 +40,7 @@ altitude in Earth radii<https://WikiPedia.org/wiki/Azimuthal_equidistant_project
 #/media/File:Comparison_azimuthal_projections.svg>}.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 # from pygeodesy.basics import _isin, _xinstanceof  # from .ellipsoidalBase
 from pygeodesy.constants import EPS, EPS0, EPS1, NAN, isnon0, _umod_360, \
@@ -71,7 +71,7 @@ from pygeodesy.utily import asin1, atan1, atan2, atan2b, atan2d, \
 from math import acos, degrees, fabs, sin, sqrt
 
 __all__ = _ALL_LAZY.azimuthal
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _EPS_K         = _EPStol * _0_1  # Karney's eps_ or _EPSmin * _0_1?
 _over_horizon_ = 'over horizon'

pygeodesy/booleans.py

@@ -15,7 +15,7 @@ C{reverse-difference}, C{sum} and C{union}.
       <http://www.Inf.USI.CH/hormann/papers/Greiner.1998.ECO.pdf>}.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _isin, isodd, issubclassof, map2, \
                              _xscalar,  typename
@@ -45,7 +45,7 @@ from pygeodesy.utily import fabs, _unrollon, _Wrap
 # from math import fabs  # from .utily
 
 __all__ = _ALL_LAZY.booleans
-__version__ = '25.04.30'
+__version__ = '25.05.26'
 
 _0EPS  =  EPS  # near-zero, positive
 _EPS0  = -EPS  # near-zero, negative
@@ -199,8 +199,8 @@ class _LatLonBool(_Named):
 
     def __sub__(self, other):
         _other(self, other)
-        return self.__class__(self.y - other.y,  # classof
-                              self.x - other.x)
+        return type(self)(self.y - other.y,  # classof
+                          self.x - other.x)
 
     def _2A(self, p2, p3):
         # I{Signed} area of a triangle, I{doubled}.
@@ -314,7 +314,7 @@ class LatLonFHP(_LatLonBool):
 
     def __add__(self, other):
         _other(self, other)
-        return self.__class__(self.y + other.y, self.x + other.x)
+        return type(self)(self.y + other.y, self.x + other.x)
 
     def __mod__(self, other):  # cross product
         _other(self, other)
@@ -326,7 +326,7 @@ class LatLonFHP(_LatLonBool):
 
     def __rmul__(self, other):  # scalar product
         _xscalar(other=other)
-        return self.__class__(self.y * other, self.x * other)
+        return type(self)(self.y * other, self.x * other)
 
 #   def _edge2(self):
 #       # Return the start and end point of the
@@ -959,7 +959,7 @@ class _CompositeBase(_Named):
         # Return a new instance
         _g = kwds.get
         kwds = dict((n, _g(n, v)) for n, v in dflts.items())
-        return self.__class__(corners or (), **kwds)
+        return type(self)(corners or (), **kwds)
 
     @property_RO
     def _clipids(self):  # PYCHOK no cover
@@ -1549,7 +1549,7 @@ class _EdgeFHP(object):
             dq  = q2 - q1
             dq2 = dq * dq  # dot product, hypot2
             if dq2 > EPS2:  # like ._clip
-                T, _E  = None, _EdgeFHP  # self.__class__
+                T, _E  = None, _EdgeFHP  # type(self)
                 p1, p2 = self._p1_p2
                 ap1   = p1._2A(q1, q2)
                 ap2_1 = p2._2A(q1, q2) - ap1
@@ -1750,8 +1750,8 @@ class _BooleanBase(object):
 
     def _boolean4(self, other, op):
         # Set up a new C{Boolean[FHP|GH]}.
-        C = self.__class__
-        kwds = C._kwds(self, op)
+        C = type(self)
+        kwds = C._kwds(self, op)  # PYCHOK in _Composite[FHP|GH]
         a =  C(self, **kwds)
         b = _other(self, other)
         return a, b, C, kwds
@@ -1964,7 +1964,7 @@ def _min_max_eps2(*xs):
 def _other(this, other):
     '''(INTERNAL) Check for compatible C{type}s.
     '''
-    C = this.__class__
+    C = type(this)
     if isinstance(other, C):
         return other
     raise _IsnotError(C, other=other)

pygeodesy/cartesianBase.py

@@ -28,7 +28,7 @@ from pygeodesy.interns import _ellipsoidal_, _spherical_  # PYCHOK used!
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
 from pygeodesy.named import _name2__, _Pass
 from pygeodesy.namedTuples import LatLon4Tuple, _NamedTupleTo , Vector3Tuple, \
-                                  Vector4Tuple,  Bearing2Tuple  # PYCHOK .sphericalBase
+                                  Vector4Tuple
 # from pygeodesy.nvectorBase import _N_vector  # _MODS
 from pygeodesy.props import deprecated_method, Property, Property_RO, property_doc_, \
                             property_RO, _update_all
@@ -44,7 +44,7 @@ from pygeodesy.vector3d import Vector3d, _xyzhdlln4
 # from math import degrees, fabs, radians, sqrt  # from .fmath, .utily
 
 __all__ = _ALL_LAZY.cartesianBase
-__version__ = '25.04.28'
+__version__ = '25.05.07'
 
 _r_         = 'r'
 _resections = _MODS.into(resections=__name__)

pygeodesy/clipy.py

@@ -10,7 +10,7 @@ I{Greiner-Hormann} and L{clipSH} and L{clipSH3} I{Sutherland-Hodgeman}.
 .
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import len2,  typename
 from pygeodesy.constants import EPS, _0_0, _1_0
@@ -31,7 +31,7 @@ from pygeodesy.units import Bool, FIx, HeightX, Lat, Lon, Number_
 # from math import fabs  # from .fmath
 
 __all__ = _ALL_LAZY.clipy
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _fj_       = 'fj'
 _original_ = 'original'

pygeodesy/constants.py

@@ -1,14 +1,14 @@
 
 # -*- coding: utf-8 -*-
 
-u'''Single-instance C{float} and C{int} constants across C{pygeodesy} modules and
-related functions L{pygeodesy.float_}, L{pygeodesy.isclose}, L{pygeodesy.isfinite},
+u'''Single-instance C{float} and C{int} constants across C{pygeodesy} modules and related
+functions L{pygeodesy.float_}, L{pygeodesy.float0_}, L{pygeodesy.isclose}, L{pygeodesy.isfinite},
 L{pygeodesy.isinf}, L{pygeodesy.isint0}, L{pygeodesy.isnan}, L{pygeodesy.isnear0},
 L{pygeodesy.isnear1}, L{pygeodesy.isnear90}, L{pygeodesy.isneg0}, L{pygeodesy.isninf},
 L{pygeodesy.isnon0} and L{pygeodesy.remainder}.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _copysign, isbool, iscomplex, isint
 from pygeodesy.errors import _xError, _xError2, _xkwds_get1, _xkwds_item2
@@ -26,7 +26,7 @@ except ImportError:  # Python 2-
     _inf, _nan = float(_INF_), float(_NAN_)
 
 __all__ = _ALL_LAZY.constants
-__version__ = '25.04.14'
+__version__ = '25.05.26'
 
 
 def _copysign_0_0(y):
@@ -495,12 +495,12 @@ def _umod_PI2(rad):
 
 if __name__ == _DMAIN_:
 
-    def _main(locals):
+    def _main(globalocals):
         from pygeodesy import itemsorted, printf
         from pygeodesy.interns import _DALL_, _UNDER_
 
         t = n = v = []
-        for n, v in itemsorted(locals):
+        for n, v in itemsorted(globalocals):
             if isinstance(v, (Float, Int, Radius)):
                 printf('%9s: %r', n, v.toRepr(std=False))
                 if v.name != n:
@@ -512,7 +512,7 @@ if __name__ == _DMAIN_:
         t.append(typename(float_))
         printf('%s = %r', _DALL_, tuple(t))
 
-    _main(locals())
+    _main(globals())  # or locals()
 
 # **) MIT License
 #

pygeodesy/css.py

@@ -470,7 +470,7 @@ class Css(_NamedBase):
 #                                h=self.height, cs0=self.cs0,
 #                                name=_name__(name, _or_nameof(self)))
 #       args, kwds = _args_kwds(**kwds)
-#       return self.__class__(*args, **kwds)  # .classof
+#       return type(self)(*args, **kwds)  # .classof
 
     @Property_RO
     def easting(self):

pygeodesy/datums.py

@@ -65,7 +65,7 @@ datum, q.v. U{"A Guide to Coordinate Systems in Great Britain", Section 6
 @var Transforms.WGS84: Transform(name='WGS84', tx=0.0, ty=0.0, tz=0.0, s1=1.0, rx=0.0, ry=0.0, rz=0.0, s=0.0, sx=0.0, sy=0.0, sz=0.0)
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 from pygeodesy.basics import _isin, islistuple, map2, neg, _xinstanceof, _zip
 from pygeodesy.constants import R_M, _float as _F, _0_0, _1_0, _2_0, _8_0, _3600_0
@@ -94,7 +94,7 @@ from pygeodesy.units import _isRadius, Radius_,  radians
 # import operator as _operator  # from .fmath
 
 __all__ = _ALL_LAZY.datums
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 _a_ellipsoid_ = _UNDER_(_a_, _ellipsoid_)
 _BD72_        = 'BD72'

pygeodesy/deprecated/__init__.py

@@ -10,13 +10,13 @@ bases} or C{from pygeodesy.deprecated import bases}.  Likewise for C{datum} and
 C{nvector}.
 '''
 
-from pygeodesy.deprecated.bases import *  # PYCHOK not pygeodesy.__init__
-from pygeodesy.deprecated.datum import *  # PYCHOK not pygeodesy.__init__
-from pygeodesy.deprecated.nvector import *  # PYCHOK not pygeodesy.__init__
+from pygeodesy.deprecated.bases import *  # noqa: F403 not pygeodesy.__init__
+from pygeodesy.deprecated.datum import *  # noqa: F403 not pygeodesy.__init__
+from pygeodesy.deprecated.nvector import *  # noqa: F403 not pygeodesy.__init__
 
-from pygeodesy.deprecated.classes import *  # PYCHOK expected
-from pygeodesy.deprecated.consterns import *  # PYCHOK expected
-from pygeodesy.deprecated.functions import *  # PYCHOK expected
+from pygeodesy.deprecated.classes import *  # noqa: F403
+from pygeodesy.deprecated.consterns import *  # noqa: F403
+from pygeodesy.deprecated.functions import *  # noqa: F403
 
 from pygeodesy.lazily import _ALL_ATTRS, _ALL_DEPRECATED, _lazy_import_as, _unLazy0  # _lazy_import_star
 
@@ -27,11 +27,11 @@ __all__ = (_ALL_DEPRECATED.deprecated_bases +
            _ALL_DEPRECATED.deprecated_classes +
            _ALL_DEPRECATED.deprecated_consterns +
            _ALL_DEPRECATED.deprecated_functions)
-__version__ = '24.12.31'
+__version__ = '25.05.12'
 
 if _unLazy0:
     from pygeodesy.deprecated import bases, datum, nvector, rhumbBase, \
-                                     rhumbaux, rhumbsolve, rhumbx  # PYCHOK expected
+                                     rhumbaux, rhumbsolve, rhumbx  # noqa: F401 PYCHOK expected
     __all__ += _ALL_ATTRS(_ALL_DEPRECATED.deprecated)  # DEPRECATED modules
 
 else:  # lazily import modules and exported attrs

pygeodesy/deprecated/bases.py

@@ -4,12 +4,12 @@
 u'''DEPRECATED on 2021.02.10, use (INTERNAL) module L{pygeodesy.latlonBase} instead.
 '''
 
-from pygeodesy.iters import points2  # PYCHOK exported
+from pygeodesy.iters import points2  # noqa: F401
 from pygeodesy.latlonBase import LatLonBase as _LatLonBase
 from pygeodesy.lazily import _ALL_DEPRECATED
 
 __all__ = _ALL_DEPRECATED.deprecated_bases
-__version__ = '25.04.09'
+__version__ = '25.05.12'
 
 
 class LatLonHeightBase(_LatLonBase):  # PYCHOK no cover

pygeodesy/deprecated/rhumbBase.py

@@ -4,10 +4,10 @@
 u'''DEPRECATED on 2023.11.26, use module L{pygeodesy.rhumb.bases} instead.
 '''
 
-from pygeodesy.rhumb.bases import *  # PYCHOK expected
+from pygeodesy.rhumb.bases import *  # noqa: F403
 
 __all__ = ()  # _ALL_LAZY.rhumb._bases
-__version__ = '23.11.26'
+__version__ = '25.05.12'
 
 # **) MIT License
 #

pygeodesy/deprecated/rhumbaux.py

@@ -4,11 +4,11 @@
 u'''DEPRECATED on 2023.11.26, use module L{pygeodesy.rhumb.aux_} instead.
 '''
 
-from pygeodesy.rhumb.aux_ import *  # PYCHOK expected
+from pygeodesy.rhumb.aux_ import *  # noqa: F403
 from pygeodesy.lazily import _ALL_LAZY
 
 __all__ = _ALL_LAZY.rhumb_aux_  # attrs NOT DEPRECATED
-__version__ = '23.12.02'
+__version__ = '25.05.12'
 
 # **) MIT License
 #

pygeodesy/deprecated/rhumbsolve.py

@@ -4,11 +4,11 @@
 u'''DEPRECATED on 2023.11.26, use module L{pygeodesy.rhumb.solve} instead.
 '''
 
-from pygeodesy.rhumb.solve import *  # PYCHOK expected
+from pygeodesy.rhumb.solve import *  # noqa: F403
 from pygeodesy.lazily import _ALL_LAZY
 
 __all__ = _ALL_LAZY.rhumb_solve  # attrs NOT DEPRECATED
-__version__ = '23.11.26'
+__version__ = '25.05.12'
 
 # **) MIT License
 #

pygeodesy/deprecated/rhumbx.py

@@ -4,11 +4,11 @@
 u'''DEPRECATED on 2023.11.26, use module L{pygeodesy.rhumb.ekx} instead.
 '''
 
-from pygeodesy.rhumb.ekx import *  # PYCHOK expected
+from pygeodesy.rhumb.ekx import *  # noqa: F403
 from pygeodesy.lazily import _ALL_LAZY
 
 __all__ = _ALL_LAZY.rhumb_ekx  # attrs NOT DEPRECATED
-__version__ = '23.11.30'
+__version__ = '25.05.12'
 
 # **) MIT License
 #

pygeodesy/ecef.py

@@ -60,9 +60,8 @@ plane} as opposed to I{geocentric} (ECEF) ones.
 from pygeodesy.basics import copysign0, _isin, isscalar, issubclassof, neg, map1, \
                             _xinstanceof, _xsubclassof,  typename  # _args_kwds_names
 from pygeodesy.constants import EPS, EPS0, EPS02, EPS1, EPS2, EPS_2, INT0, PI, PI_2, \
-                               _0_0, _0_0001, _0_01, _0_5, _1_0, _1_0_1T, _N_1_0, \
-                               _2_0, _N_2_0, _3_0, _4_0, _6_0, _60_0, _90_0, _N_90_0, \
-                               _100_0, _copysign_1_0,  isnon0  # PYCHOK used!
+                               _0_0, _0_5, _1_0, _1_0_1T, _2_0, _N_2_0, _3_0, _4_0, \
+                               _6_0, _90_0, _N_90_0, _copysign_1_0,  isnon0  # PYCHOK used!
 from pygeodesy.datums import _ellipsoidal_datum, _WGS84,  a_f2Tuple, _EWGS84
 from pygeodesy.ecefLocals import _EcefLocal
 # from pygeodesy.ellipsoids import a_f2Tuple, _EWGS84  # from .datums
@@ -90,7 +89,7 @@ from pygeodesy.utily import atan1, atan1d, atan2, atan2d, degrees90, degrees180,
 from math import cos, degrees, fabs, radians, sqrt
 
 __all__ = _ALL_LAZY.ecef
-__version__ = '25.04.28'
+__version__ = '25.05.07'
 
 _Ecef_    = 'Ecef'
 _prolate_ = 'prolate'

pygeodesy/ellipsoidalBase.py

@@ -10,7 +10,7 @@ and published under the same MIT Licence**, see for example U{latlon-ellipsoidal
 <https://www.Movable-Type.co.UK/scripts/geodesy/docs/latlon-ellipsoidal.js.html>}.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 # from pygeodesy.azimuthal import EquidistantExact, EquidistantKarney  # _MODS
 from pygeodesy.basics import _isin, _xinstanceof
@@ -46,7 +46,7 @@ from pygeodesy.units import Epoch, _isDegrees, Radius_, _1mm as _TOL_M
 # from math import fabs  # from .latlonBase
 
 __all__ = _ALL_LAZY.ellipsoidalBase
-__version__ = '25.04.14'
+__version__ = '25.05.12'
 
 
 class CartesianEllipsoidalBase(CartesianBase):

pygeodesy/ellipsoidalBaseDI.py

@@ -5,7 +5,7 @@ u'''(INTERNAL) Private, ellipsoidal Direct/Inverse geodesy base
 class C{LatLonEllipsoidalBaseDI} and functions.
 '''
 # make sure int/int division yields float quotient, see .basics
-from __future__ import division as _; del _  # PYCHOK semicolon
+from __future__ import division as _; del _  # noqa: E702 ;
 
 # from pygeodesy.azimuthal import _Equidistants  # _MODS
 from pygeodesy.basics import isLatLon, _xsubclassof,  typename
@@ -36,7 +36,7 @@ from pygeodesy.utily import m2km, unroll180, _unrollon, _unrollon3, _Wrap, wrap3
 from math import degrees, radians
 
 __all__ = _ALL_LAZY.ellipsoidalBaseDI
-__version__ = '25.04.21'
+__version__ = '25.05.23'
 
 _polar__ = 'polar?'
 _TRIPS   =  33   # _intersect3, _intersects2, _nearestOn interations, 6..9 sufficient?
@@ -115,14 +115,15 @@ class LatLonEllipsoidalBaseDI(LatLonEllipsoidalBase):
             r = self._Direct2Tuple(LL, height, r)
         return r
 
-    def _Direct2Tuple(self, LL, height, r):
+    def _Direct2Tuple(self, LL, height, r):  # in .ellipsoidalVincenty.py
         '''(INTERNAL) Helper for C{._Direct} result L{Destination2Tuple}.
         '''
         h =  self._heigHt(height)
-        d = _xkwds_not(None, datum=self.datum, name=self.name,
+        n =  self.name
+        d = _xkwds_not(None, datum=self.datum, height=h, name=n,
                              epoch=self.epoch, reframe=self.reframe)
-        d = LL(*_Wrap.latlon(r.lat, r.lon), height=h, **d)
-        return Destination2Tuple(d, wrap360(r.final), name=self.name)
+        d = LL(*_Wrap.latlon(r.lat, r.lon), **d)
+        return Destination2Tuple(d, wrap360(r.final), name=n)
 
     def distanceTo(self, other, wrap=False, **unused):  # radius=R_M
         '''Compute the distance between this and an other point along

pygeodesy/ellipsoidalExact.py

@@ -11,12 +11,12 @@ L{GeodesicExact}, L{GeodesicAreaExact} and L{GeodesicLineExact}.
 
 # from pygeodesy.datums import _WGS84  # from .ellipsoidalBase
 from pygeodesy.ellipsoidalBase import CartesianEllipsoidalBase, \
-                                     _nearestOn, _WGS84
+                                     _nearestOn,  _WGS84
 from pygeodesy.ellipsoidalBaseDI import LatLonEllipsoidalBaseDI, \
                                        _intersection3, _intersections2, \
                                        _TOL_M, intersecant2
 # from pygeodesy.errors import _xkwds  # from .karney
-from pygeodesy.karney import fabs, _polygon,  Property_RO, _xkwds
+from pygeodesy.karney import _polygon,  fabs, Property_RO, _xkwds
 from pygeodesy.lazily import _ALL_LAZY, _ALL_MODS as _MODS, _ALL_OTHER
 from pygeodesy.points import _areaError, ispolar  # PYCHOK exported
 # from pygeodesy.props import Property_RO  # from .karney
@@ -24,7 +24,7 @@ from pygeodesy.points import _areaError, ispolar  # PYCHOK exported
 # from math import fabs  # from .karney
 
 __all__ = _ALL_LAZY.ellipsoidalExact
-__version__ = '24.11.06'
+__version__ = '25.05.28'
 
 
 class Cartesian(CartesianEllipsoidalBase):
@@ -88,17 +88,19 @@ class LatLon(LatLonEllipsoidalBaseDI):
         return LatLonEllipsoidalBaseDI.toCartesian(self, **kwds)
 
 
-def areaOf(points, datum=_WGS84, wrap=True):
+def areaOf(points, datum=_WGS84, wrap=True, polar=False):
     '''Compute the area of an (ellipsoidal) polygon or composite.
 
-       @arg points: The polygon points (L{LatLon}[], L{BooleanFHP} or
-                    L{BooleanGH}).
+       @arg points: The polygon points (L{LatLon}[], L{BooleanFHP} or L{BooleanGH}).
        @kwarg datum: Optional datum (L{Datum}).
-       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the
-                    B{C{points}} (C{bool}).
+       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the B{C{points}}
+                    (C{bool}).
+       @kwarg polar: Use C{B{polar}=True} if the polygon encloses a pole (C{bool}), see
+                     function L{ispolar<pygeodesy.points.ispolar>} and U{area of a polygon
+                     enclosing a pole<https://GeographicLib.SourceForge.io/C++/doc/
+                     classGeographicLib_1_1GeodesicExact.html#a3d7a9155e838a09a48dc14d0c3fac525>}.
 
-       @return: Area (C{meter} I{squared}, same units as the B{C{datum}}'s
-                ellipsoid axes).
+       @return: Area (C{meter} I{squared}, same units as the B{C{datum}}'s ellipsoid axes).
 
        @raise PointsError: Insufficient number of B{C{points}}.
 
@@ -110,12 +112,8 @@ def areaOf(points, datum=_WGS84, wrap=True):
        @see: Functions L{pygeodesy.areaOf}, L{ellipsoidalGeodSolve.areaOf},
              L{ellipsoidalKarney.areaOf}, L{sphericalNvector.areaOf} and
              L{sphericalTrigonometry.areaOf}.
-
-       @note: The U{area of a polygon enclosing a pole<https://GeographicLib.SourceForge.io/
-              C++/doc/classGeographicLib_1_1GeodesicExact.html#a3d7a9155e838a09a48dc14d0c3fac525>}
-              can be found by adding half the datum's ellipsoid surface area to the polygon's area.
     '''
-    return fabs(_polygon(datum.ellipsoid.geodesicx, points, True, False, wrap))
+    return fabs(_polygon(datum.ellipsoid.geodesicx, points, True, False, wrap, polar))
 
 
 def intersection3(start1, end1, start2, end2, height=None, wrap=False,  # was=True
@@ -216,13 +214,14 @@ def intersections2(center1, radius1, center2, radius2, height=None, wrap=False,
                                              equidistant=equidistant, tol=tol, **kwds)
 
 
-def isclockwise(points, datum=_WGS84, wrap=True):
+def isclockwise(points, datum=_WGS84, wrap=True, polar=False):
     '''Determine the direction of a path or polygon.
 
        @arg points: The path or polygon points (C{LatLon}[]).
        @kwarg datum: Optional datum (L{Datum}).
-       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the
-                    B{C{points}} (C{bool}).
+       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the B{C{points}} (C{bool}).
+       @kwarg polar: Use C{B{polar}=True} if the C{B{points}} enclose a pole (C{bool}),
+                     see function U{ispolar<pygeodeys.points.ispolar>}.
 
        @return: C{True} if B{C{points}} are clockwise, C{False} otherwise.
 
@@ -234,7 +233,7 @@ def isclockwise(points, datum=_WGS84, wrap=True):
 
        @see: L{pygeodesy.isclockwise}.
     '''
-    a = _polygon(datum.ellipsoid.geodesicx, points, True, False, wrap)
+    a = _polygon(datum.ellipsoid.geodesicx, points, True, False, wrap, polar)
     if a < 0:
         return True
     elif a > 0:
@@ -288,29 +287,25 @@ def nearestOn(point, point1, point2, within=True, height=None, wrap=False,
 def perimeterOf(points, closed=False, datum=_WGS84, wrap=True):
     '''Compute the perimeter of an (ellipsoidal) polygon or composite.
 
-       @arg points: The polygon points (L{LatLon}[], L{BooleanFHP} or
-                    L{BooleanGH}).
+       @arg points: The polygon points (L{LatLon}[], L{BooleanFHP} or L{BooleanGH}).
        @kwarg closed: Optionally, close the polygon (C{bool}).
        @kwarg datum: Optional datum (L{Datum}).
-       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the
-                    B{C{points}} (C{bool}).
+       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the B{C{points}} (C{bool}).
 
-       @return: Perimeter (C{meter}, same units as the B{C{datum}}'s
-                ellipsoid axes).
+       @return: Perimeter (C{meter}, same units as the B{C{datum}}'s ellipsoid axes).
 
        @raise PointsError: Insufficient number of B{C{points}}.
 
        @raise TypeError: Some B{C{points}} are not L{LatLon}.
 
-       @raise ValueError: Invalid C{B{wrap}=False}, unwrapped, unrolled
-                          longitudes not supported or C{B{closed}=False}
-                          with C{B{points}} a composite.
+       @raise ValueError: Invalid C{B{wrap}=False}, unwrapped, unrolled longitudes not
+                          supported or C{B{closed}=False} with C{B{points}} a composite.
 
        @see: Functions L{pygeodesy.perimeterOf}, L{ellipsoidalGeodSolve.perimeterOf},
              L{ellipsoidalKarney.perimeterOf}, L{sphericalNvector.perimeterOf} and
              L{sphericalTrigonometry.perimeterOf}.
     '''
-    return _polygon(datum.ellipsoid.geodesicx, points, closed, True, wrap)
+    return _polygon(datum.ellipsoid.geodesicx, points, closed, True, wrap, False)
 
 
 __all__ += _ALL_OTHER(Cartesian, LatLon,  # classes

pygeodesy/ellipsoidalGeodSolve.py

@@ -24,7 +24,7 @@ from pygeodesy.points import _areaError, ispolar  # PYCHOK exported
 # from math import fabs  # from .karney
 
 __all__ = _ALL_LAZY.ellipsoidalGeodSolve
-__version__ = '24.08.13'
+__version__ = '25.05.28'
 
 
 class Cartesian(CartesianEllipsoidalBase):
@@ -87,12 +87,16 @@ class LatLon(LatLonEllipsoidalBaseDI):
         return LatLonEllipsoidalBaseDI.toCartesian(self, **kwds)
 
 
-def areaOf(points, datum=_WGS84, wrap=True):
+def areaOf(points, datum=_WGS84, wrap=True, polar=False):
     '''Compute the area of an (ellipsoidal) polygon or composite.
 
        @arg points: The polygon points (L{LatLon}[], L{BooleanFHP} or L{BooleanGH}).
        @kwarg datum: Optional datum (L{Datum}).
        @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the B{C{points}} (C{bool}).
+       @kwarg polar: Use C{B{polar}=True} if the polygon encloses a pole (C{bool}), see
+                     function L{ispolar<pygeodesy.points.ispolar>} and U{area of a polygon
+                     enclosing a pole<https://GeographicLib.SourceForge.io/C++/doc/
+                     classGeographicLib_1_1GeodesicExact.html#a3d7a9155e838a09a48dc14d0c3fac525>}.
 
        @return: Area (C{meter}, same as units of the B{C{datum}}'s ellipsoid axes, I{squared}).
 
@@ -105,7 +109,7 @@ def areaOf(points, datum=_WGS84, wrap=True):
        @see: Functions L{pygeodesy.areaOf}, L{ellipsoidalExact.areaOf}, L{ellipsoidalKarney.areaOf},
              L{sphericalNvector.areaOf} and L{sphericalTrigonometry.areaOf}.
     '''
-    return fabs(_polygon(datum.ellipsoid.geodsolve, points, True, False, wrap))
+    return fabs(_polygon(datum.ellipsoid.geodsolve, points, True, False, wrap, polar))
 
 
 def intersection3(start1, end1, start2, end2, height=None, wrap=False,  # was=True
@@ -209,13 +213,14 @@ def intersections2(center1, radius1, center2, radius2, height=None, wrap=False,
                            equidistant=equidistant, tol=tol, LatLon=LatLon, **LatLon_kwds)
 
 
-def isclockwise(points, datum=_WGS84, wrap=True):
+def isclockwise(points, datum=_WGS84, wrap=True, polar=False):
     '''Determine the direction of a path or polygon.
 
        @arg points: The path or polygon points (C{LatLon}[]).
        @kwarg datum: Optional datum (L{Datum}).
-       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the
-                    B{C{points}} (C{bool}).
+       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the B{C{points}} (C{bool}).
+       @kwarg polar: Use C{B{polar}=True} if the C{B{points}} enclose a pole (C{bool}),
+                     see function U{ispolar<pygeodeys.points.ispolar>}.
 
        @return: C{True} if B{C{points}} are clockwise, C{False} otherwise.
 
@@ -223,12 +228,11 @@ def isclockwise(points, datum=_WGS84, wrap=True):
 
        @raise TypeError: Some B{C{points}} are not C{LatLon}.
 
-       @raise ValueError: The B{C{points}} enclose a pole or zero
-                          area.
+       @raise ValueError: The B{C{points}} enclose a pole or zero area.
 
        @see: L{pygeodesy.isclockwise}.
     '''
-    a = _polygon(datum.ellipsoid.geodsolve, points, True, False, wrap)
+    a = _polygon(datum.ellipsoid.geodsolve, points, True, False, wrap, polar)
     if a < 0:
         return True
     elif a > 0:
@@ -282,29 +286,25 @@ def nearestOn(point, point1, point2, within=True, height=None, wrap=False,
 def perimeterOf(points, closed=False, datum=_WGS84, wrap=True):
     '''Compute the perimeter of an (ellipsoidal) polygon or composite.
 
-       @arg points: The polygon points (L{LatLon}[], L{BooleanFHP} or
-                    L{BooleanGH}).
+       @arg points: The polygon points (L{LatLon}[], L{BooleanFHP} or L{BooleanGH}).
        @kwarg closed: Optionally, close the polygon (C{bool}).
        @kwarg datum: Optional datum (L{Datum}).
-       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the
-                    B{C{points}} (C{bool}).
+       @kwarg wrap: If C{True}, wrap or I{normalize} and unroll the B{C{points}} (C{bool}).
 
-       @return: Perimeter (C{meter}, same as units of the B{C{datum}}'s
-                ellipsoid axes).
+       @return: Perimeter (C{meter}, same as units of the B{C{datum}}'s ellipsoid axes).
 
        @raise PointsError: Insufficient number of B{C{points}}.
 
        @raise TypeError: Some B{C{points}} are not L{LatLon}.
 
-       @raise ValueError: Invalid C{B{wrap}=False}, unwrapped, unrolled
-                          longitudes not supported or C{B{closed}=False}
-                          with C{B{points}} a composite.
+       @raise ValueError: Invalid C{B{wrap}=False}, unwrapped, unrolled longitudes not
+                          supported or C{B{closed}=False} with C{B{points}} a composite.
 
        @see: Functions L{pygeodesy.perimeterOf}, L{ellipsoidalExact.perimeterOf},
              L{ellipsoidalKarney.perimeterOf}, L{sphericalNvector.perimeterOf}
              and L{sphericalTrigonometry.perimeterOf}.
     '''
-    return _polygon(datum.ellipsoid.geodsolve, points, closed, True, wrap)
+    return _polygon(datum.ellipsoid.geodsolve, points, closed, True, wrap, False)
 
 
 __all__ += _ALL_OTHER(Cartesian, LatLon,  # classes