Adding support for downrange-crossrange-above coordinate conversions with dca.py

Author: ryanpavlick

.github/contributors.md

@@ -2,7 +2,7 @@ Thanks to those who contributed code and ideas, including:
 
 ```
 @aldebaran1 (robustness)
-@rpavlick (multiple features + functions)
+@ryanpavlick (.rsphere, .rcurve, .latitude, .lox, .los, .dca, and more)
 @cchuravy (Ellipsoid parameters)
 @jprMesh  (more conversion functions)
 @Fil (docs)

src/pymap3d/__init__.py

@@ -32,6 +32,18 @@
 __version__ = "3.1.1"
 
 from .aer import aer2ecef, aer2geodetic, ecef2aer, geodetic2aer
+from .dca import (
+    enu2dca,
+    dca2enu,
+    dca2ned,
+    ned2dca,
+    ecef2dca,
+    dca2ecef,
+    geodetic2dca,
+    dca2geodetic,
+    aer2dca,
+    dca2aer,
+)
 from .ecef import (
     ecef2enu,
     ecef2enuv,

src/pymap3d/dca.py

@@ -0,0 +1,130 @@
+"""
+Transforms involving DCA (Downrange, Crossrange, Above)
+
+This module provides functions to convert coordinates between DCA and 
+other coordinate systems such as ENU (East, North, Up), NED (North, East, Down),
+ECEF (Earth-Centered, Earth-Fixed), and geodetic coordinates (latitude, longitude, altitude). 
+
+It also includes transformations to/from AER (Azimuth, Elevation, Range).
+
+Functions:
+----------
+- enu2dca: Convert ENU to DCA coordinates.
+- dca2enu: Convert DCA to ENU coordinates.
+- dca2ned: Convert DCA to NED coordinates.
+- ned2dca: Convert NED to DCA coordinates.
+- ecef2dca: Convert ECEF to DCA coordinates.
+- dca2ecef: Convert DCA to ECEF coordinates.
+- geodetic2dca: Convert geodetic to DCA coordinates.
+- dca2geodetic: Convert DCA to geodetic coordinates.
+- aer2dca: Convert AER to DCA coordinates.
+- dca2aer: Convert DCA to AER coordinates.
+"""
+
+from .mathfun import sin, cos, radians
+from .ecef import ecef2enu, enu2ecef
+from .enu import geodetic2enu, enu2geodetic, aer2enu, enu2aer
+from .ellipsoid import Ellipsoid
+
+__all__ = [
+    "enu2dca",
+    "dca2enu",
+    "dca2ned",
+    "ned2dca",
+    "ecef2dca",
+    "dca2ecef",
+    "geodetic2dca",
+    "dca2geodetic",
+    "aer2dca",
+    "dca2aer",
+]
+
+ELL = Ellipsoid.from_name("wgs84")  # Default reference ellipsoid (WGS84)
+
+
+def enu2dca(e, n, u, heading, deg=True):
+    """
+    Converts ENU (East, North, Up) coordinates to DCA (Downrange, Crossrange, Above).
+    """
+    if deg:
+        heading = radians(heading)
+    dr = e * sin(heading) + n * cos(heading)
+    cr = -e * cos(heading) + n * sin(heading)
+    a = u
+    return dr, cr, a
+
+
+def dca2enu(dr, cr, a, heading, deg=True):
+    """
+    Converts DCA (Downrange, Crossrange, Above) coordinates to ENU (East, North, Up).
+    """
+    if deg:
+        heading = radians(heading)
+    e = dr * sin(heading) - cr * cos(heading)
+    n = dr * cos(heading) + cr * sin(heading)
+    u = a
+    return e, n, u
+
+
+def dca2ned(dr, cr, a, heading, deg=True):
+    """
+    Converts DCA (Downrange, Crossrange, Above) coordinates to NED (North, East, Down).
+    """
+    e, n, u = dca2enu(dr, cr, a, heading, deg=deg)
+    return n, e, -u
+
+
+def ned2dca(n, e, d, heading, deg=True):
+    """
+    Converts NED (North, East, Down) coordinates to DCA (Downrange, Crossrange, Above).
+    """
+    dr, cr, a = enu2dca(e, n, -d, heading, deg=deg)
+    return dr, cr, a
+
+
+def ecef2dca(x, y, z, lat0, lon0, h0, heading, ell=ELL, deg=True):
+    """
+    Converts ECEF (Earth-Centered, Earth-Fixed) coordinates to DCA (Downrange, Crossrange, Above).
+    """
+    e, n, u = ecef2enu(x, y, z, lat0, lon0, h0, ell, deg=deg)
+    return enu2dca(e, n, u, heading, deg=deg)
+
+
+def dca2ecef(dr, cr, a, lat0, lon0, h0, heading, ell=ELL, deg=True):
+    """
+    Converts DCA (Downrange, Crossrange, Above) coordinates to ECEF (Earth-Centered, Earth-Fixed) coordinates.
+    """
+    e, n, u = dca2enu(dr, cr, a, heading, deg=deg)
+    return enu2ecef(e, n, u, lat0, lon0, h0, ell, deg=deg)
+
+
+def geodetic2dca(lat, lon, h, lat0, lon0, h0, heading, ell=ELL, deg=True):
+    """
+    Converts geodetic coordinates (latitude, longitude, altitude) to DCA (Downrange, Crossrange, Above) coordinates.
+    """
+    e, n, u = geodetic2enu(lat, lon, h, lat0, lon0, h0, ell, deg=deg)
+    return enu2dca(e, n, u, heading, deg=deg)
+
+
+def dca2geodetic(dr, cr, a, lat0, lon0, h0, heading, ell=ELL, deg=True):
+    """
+    Converts DCA (Downrange, Crossrange, Above) coordinates to geodetic coordinates (latitude, longitude, altitude).
+    """
+    e, n, u = dca2enu(dr, cr, a, heading, deg=deg)
+    return enu2geodetic(e, n, u, lat0, lon0, h0, ell, deg=deg)
+
+
+def aer2dca(az, el, srange, heading, deg=True):
+    """
+    Converts AER (Azimuth, Elevation, Range) coordinates to DCA (Downrange, Crossrange, Above).
+    """
+    e, n, u = aer2enu(az, el, srange, deg=deg)
+    return enu2dca(e, n, u, heading, deg=deg)
+
+
+def dca2aer(dr, cr, a, heading, deg=True):
+    """
+    Converts DCA (Downrange, Crossrange, Above) coordinates to AER (Azimuth, Elevation, Range).
+    """
+    e, n, u = dca2enu(dr, cr, a, heading, deg=deg)
+    return enu2aer(e, n, u, deg=deg)