pygeodesy 24.11.11__py2.py3-none-any.whl → 24.12.12__py2.py3-none-any.whl

pygeodesy/cartesianBase.py

@@ -24,32 +24,32 @@ from pygeodesy.fsums import fsumf_,  Fmt
 from pygeodesy.interns import _COMMASPACE_, _datum_, _no_, _phi_
 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.named import _name2__, _NamedLocal, _Pass
 from pygeodesy.namedTuples import LatLon4Tuple, _NamedTupleTo , Vector3Tuple, \
                                   Vector4Tuple,  Bearing2Tuple  # PYCHOK .sphericalBase
 # from pygeodesy.nvectorBase import _N_vector  # _MODS
 from pygeodesy.props import deprecated_method, Property, Property_RO, property_doc_, \
-                            property_RO, property_ROnce, _update_all
+                            property_RO, _update_all
 # from pygeodesy,resections import cassini, collins5, pierlot, pierlotx, \
 #                                  tienstra7  # _MODS
 # from pygeodesy.streprs import Fmt  # from .fsums
 # from pygeodesy.triaxials import Triaxial_  # _MODS
 from pygeodesy.units import Degrees, Height, _heigHt, _isMeter, Meter, Radians
-from pygeodesy.utily import acos1, sincos2d, sincos2_,  atan2, degrees, radians
+from pygeodesy.utily import acos1, atan2, sincos2d, sincos2_,  degrees, radians
 from pygeodesy.vector3d import Vector3d, _xyzhdlln4
 # from pygeodesy.vector3dBase import _xyz3  # _MODS
 # from pygeodesy import ltp  # _MODS
 
-# from math import atan2, degrees, fabs, radians, sqrt  # from .fmath, .utily
+# from math import degrees, fabs, radians, sqrt  # from .fmath, .utily
 
 __all__ = _ALL_LAZY.cartesianBase
-__version__ = '24.11.06'
+__version__ = '24.12.04'
 
 _r_     = 'r'
 _theta_ = 'theta'
 
 
-class CartesianBase(Vector3d):
+class CartesianBase(Vector3d, _NamedLocal):
     '''(INTERNAL) Base class for ellipsoidal and spherical C{Cartesian}.
     '''
     _datum  = None  # L{Datum}, to be overriden
@@ -211,12 +211,6 @@ class CartesianBase(Vector3d):
             r = Cartesian(*c, **kwds)
         return r.renamed(n) if n else r
 
-    @property_ROnce
-    def Ecef(self):
-        '''Get the ECEF I{class} (L{EcefKarney}), I{once}.
-        '''
-        return _MODS.ecef.EcefKarney
-
     @Property_RO
     def _ecef9(self):
         '''(INTERNAL) Helper for L{toEcef}, L{toLocal} and L{toLtp} (L{Ecef9Tuple}).
@@ -231,14 +225,14 @@ class CartesianBase(Vector3d):
 
     def hartzell(self, los=False, earth=None):
         '''Compute the intersection of a Line-Of-Sight from this cartesian Point-Of-View
-           (pov) and this cartesian's ellipsoid surface.
+           (pov) and this cartesian's C{datum} ellipsoid surface.
 
            @kwarg los: Line-Of-Sight, I{direction} to the ellipsoid (L{Los}, L{Vector3d}),
                        C{True} for the I{normal, plumb} onto the surface or I{False} or
                        C{None} to point to the center of the ellipsoid.
            @kwarg earth: The earth model (L{Datum}, L{Ellipsoid}, L{Ellipsoid2}, L{a_f2Tuple}
                          or C{scalar} radius in C{meter}), overriding this cartesian's
-                         C{datum} ellipsoid.
+                         datum.
 
            @return: The intersection (C{Cartesian}) with C{.height} set to the distance to
                     this C{pov}.
@@ -284,7 +278,8 @@ class CartesianBase(Vector3d):
         return r
 
     def height3(self, earth=None, height=None, **Cartesian_and_kwds):
-        '''Compute the cartesian at a height above or below this certesian's ellipsoid.
+        '''Compute the cartesian at a height above or below this certesian's
+           C{datum} ellipsoid surface.
 
            @kwarg earth: A datum, ellipsoid, triaxial ellipsoid or earth radius,
                          I{overriding} this cartesian's datum (L{Datum}, L{Ellipsoid},
@@ -411,12 +406,6 @@ class CartesianBase(Vector3d):
         '''
         return self.toEcef().latlonheightdatum
 
-    @Property_RO
-    def _Ltp(self):
-        '''(INTERNAL) Cache for L{toLtp}.
-        '''
-        return _MODS.ltp.Ltp(self._ecef9, ecef=self.Ecef(self.datum), name=self.name)
-
     @Property_RO
     def _N_vector(self):
         '''(INTERNAL) Get the (C{nvectorBase._N_vector_}).
@@ -696,8 +685,8 @@ class CartesianBase(Vector3d):
     def toEcef(self):
         '''Convert this cartesian to I{geodetic} (lat-/longitude) coordinates.
 
-           @return: An L{Ecef9Tuple}C{(x, y, z, lat, lon, height,
-                    C, M, datum)} with C{C} and C{M} if available.
+           @return: An L{Ecef9Tuple}C{(x, y, z, lat, lon, height, C, M, datum)}
+                    with C{C} and C{M} if available.
 
            @raise EcefError: A C{.datum} or an ECEF issue.
         '''
@@ -706,18 +695,16 @@ class CartesianBase(Vector3d):
     def toLatLon(self, datum=None, height=None, LatLon=None, **LatLon_kwds):  # see .ecef.Ecef9Tuple.toDatum
         '''Convert this cartesian to a I{geodetic} (lat-/longitude) point.
 
-           @kwarg datum: Optional datum (L{Datum}, L{Ellipsoid}, L{Ellipsoid2}
-                         or L{a_f2Tuple}).
-           @kwarg height: Optional height, overriding the converted height
-                          (C{meter}), only if C{B{LatLon} is not None}.
-           @kwarg LatLon: Optional class to return the geodetic point
-                          (C{LatLon}) or C{None}.
-           @kwarg LatLon_kwds: Optional, additional B{C{LatLon}} keyword
-                               arguments, ignored if C{B{LatLon} is None}.
+           @kwarg datum: Optional datum (L{Datum}, L{Ellipsoid}, L{Ellipsoid2} or L{a_f2Tuple}).
+           @kwarg height: Optional height, overriding the converted height (C{meter}), only if
+                          C{B{LatLon} is not None}.
+           @kwarg LatLon: Optional class to return the geodetic point (C{LatLon}) or C{None}.
+           @kwarg LatLon_kwds: Optional, additional B{C{LatLon}} keyword arguments, ignored if
+                               C{B{LatLon} is None}.
 
-           @return: The geodetic point (B{C{LatLon}}) or if C{B{LatLon}
-                    is None}, an L{Ecef9Tuple}C{(x, y, z, lat, lon, height,
-                    C, M, datum)} with C{C} and C{M} if available.
+           @return: The geodetic point (B{C{LatLon}}) or if C{B{LatLon}is None}, an
+                    L{Ecef9Tuple}C{(x, y, z, lat, lon, height, C, M, datum)} with C{C}
+                    and C{M} if available.
 
            @raise TypeError: Invalid B{C{datum}} or B{C{LatLon_kwds}}.
         '''
@@ -735,32 +722,6 @@ class CartesianBase(Vector3d):
         _xdatum(r.datum, d)
         return r
 
-    def toLocal(self, Xyz=None, ltp=None, **Xyz_kwds):
-        '''Convert this I{geocentric} cartesian to I{local} C{X}, C{Y} and C{Z}.
-
-           @kwarg Xyz: Optional class to return C{X}, C{Y} and C{Z} (L{XyzLocal},
-                       L{Enu}, L{Ned}) or C{None}.
-           @kwarg ltp: The I{local tangent plane} (LTP) to use, overriding this
-                       cartesian's LTP (L{Ltp}).
-           @kwarg Xyz_kwds: Optional, additional B{C{Xyz}} keyword arguments,
-                            ignored if C{B{Xyz} is None}.
-
-           @return: An B{C{Xyz}} instance or a L{Local9Tuple}C{(x, y, z, lat, lon,
-                    height, ltp, ecef, M)} if C{B{Xyz} is None} (with C{M=None}).
-
-           @raise TypeError: Invalid B{C{ltp}}.
-        '''
-        return _MODS.ltp._toLocal(self, ltp, Xyz, Xyz_kwds)  # self._ecef9
-
-    def toLtp(self, Ecef=None, **name):
-        '''Return the I{local tangent plane} (LTP) for this cartesian.
-
-           @kwarg Ecef: Optional ECEF I{class} (L{EcefKarney}, ...
-                        L{EcefYou}), overriding this cartesian's C{Ecef}.
-           @kwarg name: Optional C{B{name}=NN} (C{str}).
-        '''
-        return _MODS.ltp._toLtp(self, Ecef, self._ecef9, name)  # self._Ltp
-
     def toNvector(self, Nvector=None, datum=None, **name_Nvector_kwds):
         '''Convert this cartesian to C{n-vector} components, I{including height}.
 
@@ -1030,7 +991,7 @@ __all__ += _ALL_DOCS(CartesianBase)
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/clipy.py

@@ -683,7 +683,7 @@ def clipSH3(points, corners, closed=False, inull=False):
 
 # **) MIT License
 #
-# Copyright (C) 2018-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2018-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/constants.py

@@ -348,8 +348,8 @@ def isfinite(obj):
        @arg obj: Value (C{scalar}, C{complex}, an L{Fsum} or
                         L{Fsum2Tuple}).
 
-       @return: C{False} if B{C{obj}} is C{INF}, C{NINF}
-                or C{NAN}, C{True} otherwise.
+       @return: C{False} if B{C{obj}} is C{INF}, C{NINF} or
+                C{NAN}, C{True} otherwise.
 
        @raise TypeError: Non-scalar and non-complex B{C{obj}}.
     '''
@@ -515,7 +515,7 @@ if __name__ == '__main__':
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/css.py

@@ -642,7 +642,7 @@ def toCss(latlon, cs0=None, height=None, Css=Css, **name):
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/datums.py

@@ -733,7 +733,7 @@ if __name__ == '__main__':
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/__init__.py

@@ -27,7 +27,7 @@ __all__ = (_ALL_DEPRECATED.deprecated_bases +
            _ALL_DEPRECATED.deprecated_classes +
            _ALL_DEPRECATED.deprecated_consterns +
            _ALL_DEPRECATED.deprecated_functions)
-__version__ = '24.11.07'
+__version__ = '24.11.28'
 
 if _unLazy0:
     from pygeodesy.deprecated import bases, datum, nvector, rhumbBase, \
@@ -40,7 +40,7 @@ else:  # lazily import modules and exported attrs
 
 # **) MIT License
 #
-# Copyright (C) 2018-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2018-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/bases.py

@@ -19,7 +19,7 @@ class LatLonHeightBase(_LatLonBase):  # PYCHOK no cover
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/classes.py

@@ -275,7 +275,7 @@ __all__ += _ALL_DOCS(_Deprecated_NamedTuple)
 
 # **) MIT License
 #
-# Copyright (C) 2018-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2018-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/consterns.py

@@ -33,7 +33,7 @@ OK     = _Deprecated_Str(OK='OK')
 
 # **) MIT License
 #
-# Copyright (C) 2018-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2018-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/datum.py

@@ -19,7 +19,7 @@ assert _ALL_OTHER(Curvature2Tuple, Datum,  Ellipsoid,  Transform) + \
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/functions.py

@@ -11,10 +11,10 @@ from pygeodesy.interns import NN, _area_, _COMMASPACE_, _negative_, \
 from pygeodesy.lazily import _ALL_DEPRECATED, _ALL_MODS as _MODS
 from pygeodesy.props import deprecated_function
 # from pygeodesy.resections import TriAngle5Tuple as _TriAngle5Tuple  # from .classes
-from pygeodesy.units import Number_, Scalar_
+from pygeodesy.units import _1mm, Number_, Scalar_
 
 __all__ = _ALL_DEPRECATED.deprecated_functions
-__version__ = '24.11.07'
+__version__ = '24.11.28'
 
 _WGS84 = _UTM = object()
 
@@ -295,11 +295,21 @@ def scalar(value, low=EPS, high=1.0, name=_scalar_, Error=ValueError):  # PYCHOK
 
 
 @deprecated_function
-def simplify2(points, pipe, radius=R_M, shortest=False, indices=False, **options):  # PYCHOK no cover
-    '''DEPRECATED, use function L{pygeodesy.simplifyRW}.
-    '''
-    return _MODS.simplify.simplifyRW(points, pipe, radius=radius, shortest=shortest,
-                                                   indices=indices, **options)
+def simplify2(points, pipe=_1mm, **radius_shortest_indices_options):  # PYCHOK no cover
+    '''DEPRECATED, use function L{pygeodesy.simplifyRW}.'''
+    return _MODS.simplify.simplifyRW(points, pipe=pipe, **radius_shortest_indices_options)
+
+
+@deprecated_function
+def simplifyRDPm(points, distance=_1mm, **radius_shortest_indices_options):
+    '''DEPRECATED on 2024.11.28, use function L{pygeodesy.simplifyRDP}C{(points, modified=True, ...)}.'''
+    return _MODS.simplify.simplifyRDP(points, modified=True, distance=distance, **radius_shortest_indices_options)
+
+
+@deprecated_function
+def simplifyVWm(points, area=_1mm, **radius_attr_indices_options):
+    '''DEPRECATED on 2024.11.28, use function L{pygeodesy.simplifyVW}C{(points, modified=True, ...)}.'''
+    return _MODS.simplify.simplifyVW(points, modified=True, area=area, **radius_attr_indices_options)
 
 
 @deprecated_function
@@ -366,7 +376,7 @@ def utmZoneBand2(lat, lon):  # PYCHOK no cover
 
 # **) MIT License
 #
-# Copyright (C) 2018-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2018-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/nvector.py

@@ -27,7 +27,7 @@ assert (_ALL_OTHER(LatLonNvectorBase, Nvector, sumOf) +
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/rhumbBase.py

@@ -11,7 +11,7 @@ __version__ = '23.11.26'
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/rhumbaux.py

@@ -12,7 +12,7 @@ __version__ = '23.12.02'
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/rhumbsolve.py

@@ -12,7 +12,7 @@ __version__ = '23.11.26'
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/deprecated/rhumbx.py

@@ -12,7 +12,7 @@ __version__ = '23.11.30'
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/dms.py

@@ -941,7 +941,7 @@ def toDMS(deg, form=_F_DMS, prec=2, sep=S_SEP, ddd=2, neg=_MINUS_, pos=_PLUS_, *
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/ecef.py

@@ -73,20 +73,23 @@ from pygeodesy.interns import NN, _a_, _C_, _datum_, _ellipsoid_, _f_, _height_,
                              _lat_, _lon_, _M_, _name_, _singular_, _SPACE_, \
                              _x_, _xyz_, _y_, _z_
 from pygeodesy.lazily import _ALL_DOCS, _ALL_LAZY, _ALL_MODS as _MODS
-from pygeodesy.named import _name__, _name1__, _NamedBase, _NamedTuple, _Pass, _xnamed
+from pygeodesy.named import _name__, _name1__, _NamedBase, _NamedLocal, \
+                            _NamedTuple, _Pass, _xnamed
 from pygeodesy.namedTuples import LatLon2Tuple, LatLon3Tuple, \
                                   PhiLam2Tuple, Vector3Tuple, Vector4Tuple
-from pygeodesy.props import deprecated_method, Property_RO, property_ROver, property_doc_
+from pygeodesy.props import deprecated_method, Property_RO, property_RO, \
+                            property_ROver, property_doc_
 # from pygeodesy.streprs import Fmt, unstr  # from .fsums
 from pygeodesy.units import _isRadius, Degrees, Height, Int, Lam, Lat, Lon, Meter, \
                              Phi, Scalar, Scalar_
-from pygeodesy.utily import atan1, atan1d, atan2d, degrees90, degrees180, \
+from pygeodesy.utily import atan1, atan1d, atan2, atan2d, degrees90, degrees180, \
                             sincos2, sincos2_, sincos2d, sincos2d_
+# from pygeodesy.vector3d import Vector3d  # _MODS
 
-from math import atan2, cos, degrees, fabs, radians, sqrt
+from math import cos, degrees, fabs, radians, sqrt
 
 __all__ = _ALL_LAZY.ecef
-__version__ = '24.10.28'
+__version__ = '24.12.06'
 
 _Ecef_    = 'Ecef'
 _prolate_ = 'prolate'
@@ -266,7 +269,7 @@ class _EcefBase(_NamedBase):
         '''(INTERNAL) Get the valid geocentric classes. I{once}.
         '''
         return (Ecef9Tuple,  # overwrite property_ROver
-               _MODS.cartesianBase.CartesianBase)
+               _MODS.vector3d.Vector3d)  # _MODS.cartesianBase.CartesianBase
 
     @Property_RO
     def _isYou(self):
@@ -326,7 +329,8 @@ class EcefFarrell21(_EcefBase):
     def reverse(self, xyz, y=None, z=None, M=None, **lon00_name):  # PYCHOK unused M
         '''Convert from geocentric C{(x, y, z)} to geodetic C{(lat, lon, height)} using
            I{Farrell}'s U{Table 2.1<https://Books.Google.com/books?id=fW4foWASY6wC>},
-           page 29.
+           page 29, aka the I{Heikkinen application} of U{Ferrari's solution
+           <https://WikiPedia.org/wiki/Geographic_coordinate_conversion>}.
 
            @arg xyz: A geocentric (C{Cartesian}, L{Ecef9Tuple}) or C{scalar} ECEF C{x}
                      coordinate (C{meter}).
@@ -366,9 +370,9 @@ class EcefFarrell21(_EcefBase):
             p2 = p**2
             G  = p2 + ez - e2 * (a2 - b2)  # p2 + ez - e4 * a2
             F  = b2 * z2 * 54
-            t  = e4 * p2 * F / G**3
-            t  = cbrt(sqrt(t * (t + _2_0)) + t + _1_0)
-            G *= fsumf_(t , _1_0, _1_0 / t)
+            c  = e4 * p2 * F / G**3
+            s  = cbrt(sqrt(c * (c + _2_0)) + c + _1_0)
+            G *= fsumf_(s , _1_0, _1_0 / s)  # k
             P  = F / (G**2 * _3_0)
             Q  = sqrt(_2_0 * e4 * P + _1_0)
             Q1 = Q +  _1_0
@@ -376,7 +380,7 @@ class EcefFarrell21(_EcefBase):
                                                -P * ez / (Q * Q1),
                                                -P * p2 * _0_5))
             r = p + e2 * r0
-            v = b2 / (sqrt(r**2 + ez) * a)
+            v = b2 / (sqrt(r**2 + ez) * a)  # z0 / z
 
             h   = hypot(r, z) * (_1_0 - v)
             lat = atan1d((e2_ * v + _1_0) * z, p)
@@ -1014,12 +1018,12 @@ class EcefMatrix(_NamedTuple):
         return xyz_
 
 
-class Ecef9Tuple(_NamedTuple):
-    '''9-Tuple C{(x, y, z, lat, lon, height, C, M, datum)} with I{geocentric}
-       C{x}, C{y} and C{z} plus I{geodetic} C{lat}, C{lon} and C{height}, case
-       C{C} (see the C{Ecef*.reverse} methods) and optionally, the rotation
-       matrix C{M} (L{EcefMatrix}) and C{datum}, with C{lat} and C{lon} in
-       C{degrees} and C{x}, C{y}, C{z} and C{height} in C{meter}, conventionally.
+class Ecef9Tuple(_NamedTuple, _NamedLocal):
+    '''9-Tuple C{(x, y, z, lat, lon, height, C, M, datum)} with I{geocentric} C{x},
+       C{y} and C{z} plus I{geodetic} C{lat}, C{lon} and C{height}, case C{C} (see
+       the C{Ecef*.reverse} methods) and optionally, rotation matrix C{M} (L{EcefMatrix})
+       and C{datum}, with C{lat} and C{lon} in C{degrees} and C{x}, C{y}, C{z} and
+       C{height} in C{meter}, conventionally.
     '''
     _Names_ = (_x_,   _y_,   _z_,   _lat_, _lon_, _height_, _C_,  _M_,   _datum_)
     _Units_ = ( Meter, Meter, Meter, Lat,   Lon,   Height,   Int, _Pass, _Pass)
@@ -1035,6 +1039,10 @@ class Ecef9Tuple(_NamedTuple):
         '''DEPRECATED, use method L{toDatum}.'''
         return self.toDatum(datum2)
 
+    @property_RO
+    def _ecef9(self):
+        return self
+
     @Property_RO
     def lam(self):
         '''Get the longitude in C{radians} (C{float}).
@@ -1088,12 +1096,19 @@ class Ecef9Tuple(_NamedTuple):
     @Property_RO
     def lonVermeille(self):
         '''Get the longitude in C{degrees [-225..+225]} after U{Vermeille
-           <https://Search.ProQuest.com/docview/639493848>} (2004), p 95.
+           <https://Search.ProQuest.com/docview/639493848>} 2004, p 95.
 
            @see: Property C{lamVermeille}.
         '''
         return Lon(Vermeille=degrees(self.lamVermeille))
 
+    def _ltp_toLocal(self, ltp, Cls_and_kwds, **nam_Cls):  # for C{_NamedLocal}
+        '''(INTERNAL) Invoke C{ltp._xLtp(ltp)._ecef2local}.
+        '''
+        Cls, kwds = self._ltp_toLocal2(Cls_and_kwds, **nam_Cls)  # in ._NamedLocal
+        ltp = self._ltp._xLtp(ltp, self._Ltp)  # both in ._NamedLocal
+        return ltp._ecef2local(self, Cls, kwds)
+
     @Property_RO
     def phi(self):
         '''Get the latitude in C{radians} (C{float}).
@@ -1134,13 +1149,13 @@ class Ecef9Tuple(_NamedTuple):
                              L{ellipsoidalNvector.Cartesian}, L{ellipsoidalVincenty.Cartesian},
                              L{sphericalNvector.Cartesian} or L{sphericalTrigonometry.Cartesian})
                              or C{None}.
-           @kwarg Cartesian_kwds: Optional, additional B{C{Cartesian}} keyword arguments, ignored
+           @kwarg Cartesian_kwds: Optionally, additional B{C{Cartesian}} keyword arguments, ignored
                                   if C{B{Cartesian} is None}.
 
-           @return: A C{B{Cartesian}(x, y, z, **B{Cartesian_kwds})} instance or
-                    a L{Vector4Tuple}C{(x, y, z, h)} if C{B{Cartesian} is None}.
+           @return: A B{C{Cartesian}} instance or a L{Vector4Tuple}C{(x, y, z, h)} if C{B{Cartesian}
+                    is None}.
 
-           @raise TypeError: Invalid B{C{Cartesian}} or B{C{Cartesian_kwds}}.
+           @raise TypeError: Invalid B{C{Cartesian}} or B{C{Cartesian_kwds}} item.
         '''
         if Cartesian in (None, Vector4Tuple):
             r = self.xyzh
@@ -1174,17 +1189,15 @@ class Ecef9Tuple(_NamedTuple):
     def toLatLon(self, LatLon=None, **LatLon_kwds):
         '''Return the geodetic C{(lat, lon, height[, datum])} coordinates.
 
-           @kwarg LatLon: Optional class to return C{(lat, lon, height[, datum])}
-                          or C{None}.
-           @kwarg LatLon_kwds: Optional B{C{height}}, B{C{datum}} and other
-                               B{C{LatLon}} keyword arguments.
+           @kwarg LatLon: Optional class to return C{(lat, lon, height[, datum])} or C{None}.
+           @kwarg LatLon_kwds: Optional B{C{height}}, B{C{datum}} and other B{C{LatLon}}
+                               keyword arguments.
 
-           @return: An instance of C{B{LatLon}(lat, lon, **B{LatLon_kwds})}
-                    or if C{B{LatLon} is None}, a L{LatLon3Tuple}C{(lat, lon,
-                    height)} respectively L{LatLon4Tuple}C{(lat, lon, height,
-                    datum)} depending on whether C{datum} is un-/specified.
+           @return: A B{C{LatLon}} instance or if C{B{LatLon} is None}, a L{LatLon4Tuple}C{(lat,
+                    lon, height, datum)} or L{LatLon3Tuple}C{(lat, lon, height)} if C{datum} is
+                    specified or not.
 
-           @raise TypeError: Invalid B{C{LatLon}} or B{C{LatLon_kwds}}.
+           @raise TypeError: Invalid B{C{LatLon}} or B{C{LatLon_kwds}} item.
         '''
         lat, lon, D = self.lat, self.lon, self.datum  # PYCHOK Ecef9Tuple
         kwds = _name1__(LatLon_kwds, _or_nameof=self)
@@ -1202,33 +1215,17 @@ class Ecef9Tuple(_NamedTuple):
         _xdatum(_xattr(r, datum=D), D)
         return r
 
-    def toLocal(self, ltp, Xyz=None, **Xyz_kwds):
-        '''Convert this geocentric to I{local} C{x}, C{y} and C{z}.
-
-           @kwarg ltp: The I{local tangent plane} (LTP) to use (L{Ltp}).
-           @kwarg Xyz: Optional class to return C{x}, C{y} and C{z}
-                       (L{XyzLocal}, L{Enu}, L{Ned}) or C{None}.
-           @kwarg Xyz_kwds: Optional, additional B{C{Xyz}} keyword
-                            arguments, ignored if C{B{Xyz} is None}.
-
-           @return: An B{C{Xyz}} instance or if C{B{Xyz} is None},
-                    a L{Local9Tuple}C{(x, y, z, lat, lon, height,
-                    ltp, ecef, M)} with C{M=None}, always.
-
-           @raise TypeError: Invalid B{C{ltp}}.
-        '''
-        return _MODS.ltp._xLtp(ltp)._ecef2local(self, Xyz, Xyz_kwds)
-
     def toVector(self, Vector=None, **Vector_kwds):
-        '''Return the geocentric C{(x, y, z)} coordinates as vector.
+        '''Return these geocentric C{(x, y, z)} coordinates as vector.
+
+           @kwarg Vector: Optional vector class to return C{(x, y, z)} or C{None}.
+           @kwarg Vector_kwds: Optional, additional B{C{Vector}} keyword arguments,
+                               ignored if C{B{Vector} is None}.
 
-           @kwarg Vector: Optional vector class to return C{(x, y, z)} or
-                          C{None}.
-           @kwarg Vector_kwds: Optional, additional B{C{Vector}} keyword
-                               arguments, ignored if C{B{Vector} is None}.
+           @return: A B{C{Vector}} instance or a L{Vector3Tuple}C{(x, y, z)} if
+                    C{B{Vector} is None}.
 
-           @return: A C{Vector}C{(x, y, z, **Vector_kwds)} instance or a
-                    L{Vector3Tuple}C{(x, y, z)} if C{B{Vector} is None}.
+           @raise TypeError: Invalid B{C{Vector}} or B{C{Vector_kwds}} item.
 
            @see: Propertes C{xyz} and C{xyzh}
         '''
@@ -1316,7 +1313,7 @@ __all__ += _ALL_DOCS(_EcefBase)
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/elevations.py

@@ -254,7 +254,7 @@ if __name__ == '__main__':
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),

pygeodesy/ellipsoidalBase.py

@@ -35,7 +35,7 @@ from pygeodesy.units import Epoch, _isDegrees, Radius_, _1mm as _TOL_M
 # from math import fabs  # from .latlonBase
 
 __all__ = _ALL_LAZY.ellipsoidalBase
-__version__ = '24.11.07'
+__version__ = '24.12.04'
 
 
 class CartesianEllipsoidalBase(CartesianBase):
@@ -666,7 +666,7 @@ class LatLonEllipsoidalBase(LatLonBase):
                         B{C{point2}} (C{bool}).
            @kwarg equidistant: An azimuthal equidistant projection (I{class} or function
                                L{pygeodesy.equidistant}) or C{None} for this point's preferred
-                               C{Equidistant}, like L{Equidistant<LatLonEllipsoidalBase.Equistant>}.
+                               C{Equidistant}, like L{Equidistant<LatLonEllipsoidalBase.Equidistant>}.
            @kwarg tol: Convergence tolerance (C{meter}, conventionally).
 
            @return: Closest point (C{LatLon}).
@@ -1148,7 +1148,7 @@ __all__ += _ALL_DOCS(CartesianEllipsoidalBase, LatLonEllipsoidalBase)
 
 # **) MIT License
 #
-# Copyright (C) 2016-2024 -- mrJean1 at Gmail -- All Rights Reserved.
+# Copyright (C) 2016-2025 -- mrJean1 at Gmail -- All Rights Reserved.
 #
 # Permission is hereby granted, free of charge, to any person obtaining a
 # copy of this software and associated documentation files (the "Software"),